Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases

Comparing newly obtained and previously known nucleotide and amino-acid sequences underpins modern biological research. BLAST is a well-established tool for such comparisons but is challenging to use on new data sets. We combined a user-centric design philosophy with sustainable software development approaches to create Sequenceserver, a tool for running BLAST and visually inspecting BLAST results for biological interpretation. Sequenceserver uses simple algorithms to prevent potential analysis errors and provides flexible text-based and visual outputs to support researcher productivity. Our software can be rapidly installed for use by individuals or on shared servers.

59 BASIC BIOLOGICAL SCIENCES↗

Identifying RNA contacts from SHAPE-MaP by partial correlation analysis

In a recent paper Siegfried et al. published a new sequence-based structural RNA assay that utilizes mutational profiling to detect base pairing (MaP). Output from MaP provides information about both pairing (via reactivities) and contact (via correlations). Reactivities can be coupled to partition function folding models for structural inference, while correlations can reveal pairs of sites that may be in structural proximity. The possibility for inference of 3D contacts via MaP suggests a novel approach to structural prediction for RNA analogous to covariance structural prediction for proteins. We explore this approach and show that partial correlation analysis outperforms na\"ive correlation analysis. Our results should be applicable to a wide range of high-throughput sequencing based RNA structural assays that are under development.

Tambe, Akshay↗

Simulation-Based Recovery Action Analysis Using the EMRALD Dynamic Risk Assessment Tool

A recovery action is defined as the action that prevents deviant conditions from producing unwanted effects. It generally indicates a kind of countermeasure performed in response to a failure of human action. The recovery actions especially play an important role in complex systems like nuclear power plants (NPPs), which consist of highly sophisticated controllers to ensure that desired performance and safety must be achieved and maintained. This is because a combination of human error and its recovery failure may be able to cause a catastrophic effect on a system. Analyzing recovery actions has been a critical part of HRA, which is a technique to evaluate human errors and provide human error probabilities (HEPs) for application in probabilistic safety assessment (PSA). If recovery actions are not adequately analyzed and applied to PSA models, the PSA results may be under-estimated or be not able to reasonably account for the failure of human actions in the context of PSA. For this reason, some regulatory documents such as ASME/ANS RA-Sb-2013 by the American Society for Mechanical Engineers and the American Nuclear Society and NUREG-1792 by U.S. Nuclear Regulatory Commission have emphasized the importance of recovery analysis within the HRA. A couple of existing HRA methods, such as the Technique for Human Error-Rate Prediction (THERP), the Cause-Based Decision Tree (CBDT), and the Korean Standard HRA (K-HRA), have respectively suggested their own approaches to the HRA recovery analysis. However, there are a couple of limitations to treating recovery actions using only the current HRA methods available. The biggest limitation is that the existing recovery analysis does not explicitly consider a variety of recovery action types and recovery sequences as they occur in actual NPPs. To handle the limitations of existing recovery analysis, this study proposes a simulation-based recovery analysis method using the Event Modeling Risk Assessment Using Linked Diagram (EMRALD) software. The EMRALD software is a dynamic simulation tool for PSA. It supports realistic and dynamic modeling of human actions as they would be performed at NPPs. It is also favorable to simultaneously model the specific moment at which an action is performed, the time it takes to perform the action, and the failure probability of that action. In this paper, a detailed methodology for modeling recovery actions in the simulation platform is proposed with a couple of examples. Then, outputs from the simulation are discussed as reviewing if this novel approach can complement the challenges of existing recovery analyses.

99 GENERAL AND MISCELLANEOUS↗

Biomass yields, reproductive fertility, compositional analysis, and genetic diversity of newly developed triploid giant miscanthus hybrids

Abstract Miscanthus × giganteus (giant miscanthus), first found as a naturally occurring hybrid, has shown promise as a bioenergy/biomass crop throughout much of the temperate world. This allotriploid (2 n = 3 x = 57) hybrid resulted from a cross between tetraploid Miscanthus sacchariflorus (2 n = 4 x = 76) and diploid Miscanthus sinensis (2 n = 2 x = 38) and is particularly desirable due to its low fertility that minimizes reseeding and potential invasiveness. However, there is limited genetic diversity in commonly grown cultivars of triploid M. × giganteus and breeding and development efforts to improve and domesticate this crop have been minimal. Here, we report on newly developed M. × giganteus hybrids compared with the industry standard M. × giganteus '1993‐1780'. Dry biomass yields of new hybrids ranged from 19.5 to 32.4 Mg/ha/year for the fourth growing season, compared with 21.0 Mg/ha/year for M. × giganteus '1993‐1780'. Plant reproductive fertility remained low for all accessions with overall fertility [(seed set × seed germination)/100] ranging from 0.3% to 4.5% for new hybrids compared to 0.4% for M. × giganteus '1993‐1780'. Culm density and height varied among accessions and were positively correlated with increased biomass. Based on compositional analyses, theoretical ethanol yields ranged from 9, 740 to 16,278 L/ha/year for new hybrids compared to 10,406 L/ha/year for M. × giganteus '1993‐1780'. Relative feed value indices were low overall and ranged between 66.0 and 72.8 for new hybrids compared to M. × giganteus '1993‐1780' with 71.3. The genetic diversity of new hybrids, compared with existing cultivars, was characterized using whole genome sequences. Based on pair‐wise distances, cluster analysis clearly showed increased diversity of new hybrids compared with earlier selections. These results document new triploid hybrids of M. × giganteus with enhanced biomass and theoretical ethanol yields in combination with broader genetic diversity and lowreproductive fertility.

Touchell, Darren H.↗

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

59 BASIC BIOLOGICAL SCIENCES↗

KBase Narrative - Supplemental Data for "Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes"

Here we provide the complete set of circular elements produced by SCAPP from metagenomic data that may represent mobile genetic elements (MGEs). For the publication, these sequences were de-replicated within the high [U] and low [U] sample sets. Additionally, 98 of the circular elements were removed for downstream analysis due to high suspicion of being (1) chloroplast or mitochondrial sequences or (2) erroneously circularized long repeat regions. The complete list of circular elements used for the analysis along with metadata can be found in Tables S5 and S6.

Goff, Jennifer↗

Single-Cell Genomics of Novel Actinobacteria With the Wood–Ljungdahl Pathway Discovered in a Serpentinizing System

Serpentinite-hosted systems represent modern-day analogs of early Earth environments. In these systems, water-rock interactions generate highly alkaline and reducing fluids that can contain hydrogen, methane, and low-molecular-weight hydrocarbons-potent reductants capable of fueling microbial metabolism. In this study, we investigated the microbiota of Hakuba Happo hot springs (~50°C; pH~10.5–11), located in Nagano (Japan), which are impacted by the serpentinization process. Analysis of the 16S rRNA gene amplicon sequences revealed that the bacterial community comprises Nitrospirae (47%), “Parcubacteria” (19%), Deinococcus-Thermus (16%), and Actinobacteria (9%), among others. Notably, only 57 amplicon sequence variants (ASV) were detected, and fifteen of these accounted for 90% of the amplicons. Among the abundant ASVs, an early-branching, uncultivated actinobacterial clade identified as RBG-16-55-12 in the SILVA database was detected. Ten single-cell genomes (average pairwise nucleotide identity: 0.98–1.00; estimated completeness: 33–93%; estimated genome size: ~2.3 Mb) that affiliated with this clade were obtained. Taxonomic classification using single copy genes indicates that the genomes belong to the actinobacterial class-level clade UBA1414 in the Genome Taxonomy Database. Based on metabolic pathway predictions, these actinobacteria are anaerobes, capable of glycolysis, dissimilatory nitrate reduction and CO 2 fixation via the Wood–Ljungdahl (WL) pathway. Several other genomes within UBA1414 and two related class-level clades also encode the WL pathway, which has not yet been reported for the Actinobacteria phylum. For the Hakuba actinobacterium, the energy metabolism related to the WL pathway is likely supported by a combination of the Rnf complex, group 3b and 3d [NiFe]-hydrogenases, [FeFe]-hydrogenases, and V-type (H + /Na + pump) ATPase. The genomes also harbor a form IV ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) complex, also known as a RubisCO-like protein, and contain signatures of interactions with viruses, including clustered regularly interspaced short palindromic repeat (CRISPR) regions and several phage integrases. This is the first report and detailed genome analysis of a bacterium within the Actinobacteria phylum capable of utilizing the WL pathway. The Hakuba actinobacterium is a member of the clade UBA1414/RBG-16-55-12, formerly within the group “OPB41.” We propose to name this bacterium ‘Candidatus Hakubanella thermoalkaliphilus.'

59 BASIC BIOLOGICAL SCIENCES↗

Portable Rabies Virus Sequencing in Canine Rabies Endemic Countries Using the Oxford Nanopore MinION

As countries with endemic canine rabies progress towards elimination by 2030, it will become necessary to employ techniques to help plan, monitor, and confirm canine rabies elimination. Sequencing can provide critical information to inform control and vaccination strategies by identifying genetically distinct virus variants that may have different host reservoir species or geographic distributions. However, many rabies testing laboratories lack the resources or expertise for sequencing, especially in remote or rural areas where human rabies deaths are highest. We developed a low-cost, high throughput rabies virus sequencing method using the Oxford Nanopore MinION portable sequencer. A total of 259 sequences were generated from diverse rabies virus isolates in public health laboratories lacking rabies virus sequencing capacity in Guatemala, India, Kenya, and Vietnam. Phylogenetic analysis provided valuable insight into rabies virus diversity and distribution in these countries and identified a new rabies virus lineage in Kenya, the first published canine rabies virus sequence from Guatemala, evidence of rabies spread across an international border in Vietnam, and importation of a rabid dog into a state working to become rabies-free in India. Taken together, our evaluation highlights the MinION’s potential for low-cost, high volume sequencing of pathogens in locations with limited resources.

59 BASIC BIOLOGICAL SCIENCES↗

Kalman filter based range estimation for autonomous navigation using imaging sensors

Rotorcraft operating in high-threat environments fly close to the surface of the earth to utilize surrounding terrain, vegetation, or man-made objects to minimize the risk of being detected by the enemy. Two basic requirements for obstacle avoidance are detection and range estimation of the object from the current rotorcraft position. There are many approaches to the estimation of range using a sequence of images. The approach used in this analysis differes from previous methods in two significant ways: an attempt is not made to estimate the rotorcraft's motion from the images; and the interest lies in recursive algorithms. The rotorcraft parameters are assumed to be computed using an onboard inertial navigation system. Given a sequence of images, using image-object differential equations, a Kalman filter (Sridhar and Phatak, 1988) can be used to estimate both the relative coordinates and the earth coordinates of the objects on the ground. The Kalman filter can also be used in a predictive mode to track features in the images, leading to a significant reduction of search effort in the feature extraction step of the algorithm. The purpose is to summarize early results obtained in extending the Kalman filter for use with actual image sequences. The experience gained from the application of this algorithm to real images is very valuable and is a necessary step before proceeding to the estimation of range during low-altitude curvilinear flight. A simple recursive method is presented to estimate range to objects using a sequence of images. The method produces good range estimates using real images in a laboratory set up and needs to be evaluated further using several different image sequences to test its robustness. The feature generation part of the algorithm requires further refinement on the strategies to limit the number of features (Sridhar and Phatak, 1989). The extension of the work reported here to curvilinear flight may require the use of the extended Kalman filter.

Sridhar, Banavar↗

Unraveling the functional dark matter through global metagenomics

Metagenomes encode an enormous diversity of proteins, reflecting a multiplicity of functions and activities1,2. Exploration of this vast sequence space has been limited to a comparative analysis against reference microbial genomes and protein families derived from those genomes. Here, to examine the scale of yet untapped functional diversity beyond what is currently possible through the lens of reference genomes, we develop a computational approach to generate reference-free protein families from the sequence space in metagenomes. We analyse 26,931 metagenomes and identify 1.17 billion protein sequences longer than 35 amino acids with no similarity to any sequences from 102,491 reference genomes or the Pfam database3. Using massively parallel graph-based clustering, we group these proteins into 106,198 novel sequence clusters with more than 100 members, doubling the number of protein families obtained from the reference genomes clustered using the same approach. We annotate these families on the basis of their taxonomic, habitat, geographical and gene neighbourhood distributions and, where sufficient sequence diversity is available, predict protein three-dimensional models, revealing novel structures. Overall, our results uncover an enormously diverse functional space, highlighting the importance of further exploring the microbial functional dark matter.

54 ENVIRONMENTAL SCIENCES↗

germs-lab/LAMPS-miscanthus-microbiome

Data and analysis to accompany our submitted paper characterizing the LAMPS Miscanthus Microbiome comparing 16S rRNA gene DNA sequencing of Miscanthus in a staggered-start experiment. Authors for this analysis include Fernando Igne Rocha, Lanying Ma, and Adina Howe.

Igne Rocha, Fernando↗

Search for evidence of a clock related to the solar 154 day complex of periodicities

Evidence that has recently been compiled (Bai and Sturrock 1991) indicates that the enigmatic 154-day periodicity in solar activity may be viewed as part of a complex of periodicities that are approximate multiples of 25.8 days, suggesting that the Sun contains a 'clock' with frequency in the range 440 to 463 nano Hz. The clock may comprise either an oscillator or a rotator, each of which may be either real or virtual. We have reconsidered a previous spectrum analysis of the Zurich sunspot-number sequence by Knight, Schatten, and Sturrock (1979) which revealed a sharp, persistent and significant periodicity with a period of 12.072 days, corresponding to a frequency of about 958.8 nano Hz. This periodicity may be regarded as the (second) upper sideband of the second harmonic (2nu(sub R) + 2nu(sub E)) of a fundamental frequency of 447.7 nano Hz that is clearly within the search band. In this expression, nu(sub R) is the sidereal frequency of the hypothetical rotator and nu(sub E) is the frequency (31.69 nano Hz) of the Earth in its orbital motion around the Sun. In analyzing sunspot area data derived from the Greenwich data set, and on noting that any frequency is defined only to within the Nyquist frequency, we find clear evidence not only for the upper sideband of the second harmonic, but also for the second harmonic (2nu(sub R)) and the lower sideband of the second harmonic (2nu(sub R) - 2nu(sub E)). There is no strong peak at the fundamental frequency in the Greenwich data, but there is in the Zurich sunspot data. The effect of a linear oscillator is, to the lowest order in the amplitude, the same as the combined effect of two rotators of opposite polarities. A rotator that has arbitrary orientation with respect to the ecliptic may influence the outer layers of the Sun and thereby modulate the occurrence of solar activity such as sunspots. By analyzing a simple model, we find that such a rotator would influence surface activity in such a way that the spectrum of a 'signal' (such as the record of sunspots), as seen from the Earth, would contain components with frequencies that are certain integral combinations of nu(sub R) and nu(sub E). The amplitudes of the various components depend sensitively on theta, the angle between the axis of the rotator and the axis of the Earth's orbital motion. This simple model therefore offers a kinematical (but not dynamical) interpretation of the sunspot spectrum. The present analysis, while offering support of our conjectures that the Sun contains a clock that regulates the 154-day complex of periodicities, cannot distinguish between an osillator or a rotator (that might be a traveling wave), nor between a real rotator or a virtual rotator (that might be an apparent traveling wave due to the aliasing effect of an oscillator in a rotating system). Further analysis of sunspot and other data sets will be required to confirm the existence of such clock and (if it is real) to determine its physical nature.

Sturrock, P. A.↗

Sequence data - Magnitude and implications of some ambiguities.

A stochastic model is applied to the divergence of the horse-pig lineage from a common ansestor in terms of the alpha and beta chains of hemoglobin and fibrinopeptides. The results are compared with those based on the minimum mutation distance model of Fitch (1972). Buckwheat and cauliflower cytochrome c sequences are analyzed to demonstrate their ambiguities. A comparative analysis of evolutionary rates for various proteins of horses and pigs shows that errors of considerable magnitude are introduced by Glx and Asx ambiguities into evolutionary conclusions drawn from sequences of incompletely analyzed proteins.

Holmquist, R.↗

Reclassification of Methanogenium tationis and Methanogenium liminatans as Methanofollis tationis gen. nov., comb. nov. and Methanofollis liminatans comb. nov. and description of a new strain of Methanofollis liminatans

Sequencing of 16S rRNA genes and phylogenetic analysis of Methanogenium tationis DSM 2702T (OCM 43T) (T = type strain) and Methanogenium liminatans GKZPZT (= DSM 4140T) as well as other members of the family Methanomicrobiaceae revealed that both species belong to a separate line of descent within this family. In addition, a new strain of Methanogenium liminatans, strain BM1 (= DSM 10196), was isolated from a butyrate-degrading, fluidized bed reactor and characterized. Cells of both species are mesophilic, highly irregular cocci that use H2/CO2 and formate for growth and methanogenesis. In addition, Methanogenium liminatans strains GKZPZT and BM1 used 2-propanol/CO2, 2-butanol/CO2 and cyclopentanol/CO2. Both species contained diether and tetraether lipids. The polar lipids comprised amino-phosphopentanetetrol derivatives, which appear to be characteristic lipids within the family Methanomicrobiaceae. The pattern of glycolipids, phosphoglycolipids and amino-phosphoglycolipids was consistent with the assignment of these two species to a taxon within the family Methanomicrobiaceae, but also permitted them to be distinguished from other higher taxa within this family. The G+C contents of the DNA of Methanogenium tationis and Methanogenium liminatans were 54 and 60 mol% (Tm and HPLC), respectively. On the basis of the data presented, the transfer of Methanogenium tationis and Methanogenium liminatans to the genus Methanofollis gen. nov. as Methanofollis tationis comb. nov. and Methanofollis liminatans comb. nov., respectively, is proposed, with Methanofollis tationis as the type species.

NASA Discipline Exobiology↗

RNAseq-based transcriptome assembly of Clostridium acetobutylicum for functional genome annotation and discovery

Accurate genome annotations are essential in modern biology and biotechnology, yet they are still largely based on genome sequencing and comparative analyses. We show that the Clostridium acetobutylicum genome annotation can be markedly improved by integrating bioinformatic predictions with RNA sequencing (RNAseq) data. Samples were acquired under butanol, butyrate, and unstressed treatments across various growth conditions. Analysis of an initial assembly revealed errors due to background signals and limitations of assembly algorithms. Hurdles for RNAseq transcriptome mapping include optimizing library complexity and sequencing depth, yet most studies report low sequencing depth and ignore the effect of ribosomal RNA abundance. An integrative analysis was developed to combine motif predictions, single-nucleotide resolution sequencing depth, and library complexity to resolve difficulties in assembly curation. This minimized false positive error and determined gene boundaries, in some cases, to the exact base-pair of prior studies. This will be the first strand-specific transcriptome assembly in a Clostridium organism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Isolation and characterization of a novel mycovirus from Penicillium digitatum

A novel double-stranded RNA virus designated Penicillium digitatum virus 1 (PdV1) was isolated from the citrus fruit rot pathogen P. digitatum (HS-RH1). The full-length cDNA sequence of the dsRNA/PdV1 (5211 bp) possesses two partially overlapping open reading frames, which encode a coat protein (CP) and a putative RNA-dependent RNA polymerase (RdRp), respectively. Phylogenetic analysis based on multiple alignments of the amino acid sequences of the RdRp and CP indicated that PdV1 tentatively belongs to the genus Victorivirus in the Totiviridae family. Electron micrographs of negatively stained viral particles purified from the peak fraction of sucrose density gradient centrifugation showed spherical particles ~35 nm in diameter. Transfection experiments with purified virions indicated that PdV1 could reduce the vegetative growth and virulence of P. digitatum strain HS-F6. In summary, we report the first isolation and characterization of a mycovirus from P. digitatum that contributes to the hypovirulence phenotypes of the host strain. - Highlights: • A novel victorivirus designated Penicillium digitatum virus 1 (PdV1) was isolated fromP. digitatum. • The dsRNA genome of PdV1 are 5211 bp long, two ORFs encoding CP and RdRp, and are encased in virions of ~35 nm in diameter. • PdV1 infection led to hypovirulent effect on P. digitatum. PdV1 may potentially be used for citrus green mold biocontrol. • Our study provides a research basis for establishing a model system for the study of P. digitatum–mycovirus interactions.

60 APPLIED LIFE SCIENCES↗

Riemann solutions for counterflow combustion in light porous foam

The paper is motivated by a model for the injection of air into a porous medium that contains a solid fuel. In previous work, a system of three evolutionary partial differential equations that models combustion of light porous foam under air injection was considered. The existence and uniqueness of traveling waves were studied and the wave sequences appearing in Riemann solutions were identified. This analysis was done under assumption that the combustion wave velocity is positive. In the present work, this hypothesis was neglected and the results generalized including the case of negative combustion wave speed. The existence of such waves was proved and the uniqueness investigated for some particular cases using Melnikov’s integral. The wave sequences appearing in Riemann solutions were identified and the numerical examples using finite difference scheme were presented.

97 MATHEMATICS AND COMPUTING↗