Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 649 records · Page 36

Methanococcus jannaschii genome: revisited

Analysis of genomic sequences is necessarily an ongoing process. Initial gene assignments tend (wisely) to be on the conservative side (Venter, 1996). The analysis of the genome then grows in an iterative fashion as additional data and more sophisticated algorithms are brought to bear on the data. The present report is an emendation of the original gene list of Methanococcus jannaschii (Bult et al., 1996). By using a somewhat more updated database and more relaxed (and operator-intensive) pattern matching methods, we were able to add significantly to, and in a few cases amend, the gene identification table originally published by Bult et al. (1996).

Non-NASA Center↗

Basin-scale study of CO 2 storage in stacked sequence of geological formations

Commercial scale decarbonization through carbon capture and storage may likely involve many CO 2 storage projects located in close proximity. The close proximity could raise concerns over caprock integrity associated with reservoir pressure buildup and interference among adjacent projects. Commercial-scale injection will also require large prospective CO 2 storage resource and high injectivity in the targeted storage formations. To accommodate the need for both large resource and high injectivity, project operators could consider injecting CO 2 into a stacked sequence of formations. This analysis investigates the benefits of injecting CO 2 into a vertically stacked sequence of saline formations, over injecting the same amount of CO 2 into a single saline formation, in addressing these challenges. Our analysis shows that injecting into the stacked sequence mitigates the extent of pressure buildup among the stacked formations, while still achieving the same or greater target CO 2 storage volumes. Among cases modeled, the resulting pressure buildup front is most reduced when each storage site distributes injection volumes over several wells, each of which injects a portion of the total CO 2 mass across the stacked sequence. This favorable case not only results in the smallest CO 2 aerial footprint, but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 46% compared to the single-formation storage), the result of which is crucial to maintain caprock integrity. This analysis provides insights into required decision-making when considering multi-project deployment in a shared basin.

42 ENGINEERING↗

CRITICA: coding region identification tool invoking comparative analysis

Gene recognition is essential to understanding existing and future DNA sequence data. CRITICA (Coding Region Identification Tool Invoking Comparative Analysis) is a suite of programs for identifying likely protein-coding sequences in DNA by combining comparative analysis of DNA sequences with more common noncomparative methods. In the comparative component of the analysis, regions of DNA are aligned with related sequences from the DNA databases; if the translation of the aligned sequences has greater amino acid identity than expected for the observed percentage nucleotide identity, this is interpreted as evidence for coding. CRITICA also incorporates noncomparative information derived from the relative frequencies of hexanucleotides in coding frames versus other contexts (i.e., dicodon bias). The dicodon usage information is derived by iterative analysis of the data, such that CRITICA is not dependent on the existence or accuracy of coding sequence annotations in the databases. This independence makes the method particularly well suited for the analysis of novel genomes. CRITICA was tested by analyzing the available Salmonella typhimurium DNA sequences. Its predictions were compared with the DNA sequence annotations and with the predictions of GenMark. CRITICA proved to be more accurate than GenMark, and moreover, many of its predictions that would seem to be errors instead reflect problems in the sequence databases. The source code of CRITICA is freely available by anonymous FTP (rdp.life.uiuc.edu in/pub/critica) and on the World Wide Web (http:/(/)rdpwww.life.uiuc.edu).

Non-NASA Center↗

Short-Read Assembly of Full-Length 16S Amplicons Reveals Bacterial Diversity in Subsurface Sediments

In microbial ecology, a fundamental question relates to how community diversity and composition change in response to perturbation. Most studies have had limited ability to deeply sample community structure (e.g. Sanger-sequenced 16S rRNA libraries), or have had limited taxonomic resolution (e.g. studies based on 16S rRNA hypervariable region sequencing). Here, we combine the higher taxonomic resolution of near-full-length 16S rRNA gene amplicons with the economics and sensitivity of short-read sequencing to assay the abundance and identity of organisms that represent as little as 0.01% of sediment bacterial communities. We used a new version of EMIRGE optimized for large data size to reconstruct near-full-length 16S rRNA genes from amplicons sheared and sequenced with Illumina technology. The approach allowed us to differentiate the community composition among samples acquired before perturbation, after acetate amendment shifted the predominant metabolism to iron reduction, and once sulfate reduction began. Results were highly reproducible across technical replicates, and identified specific taxa that responded to the perturbation. All samples contain very high alpha diversity and abundant organisms from phyla without cultivated representatives. Surprisingly, at the time points measured, there was no strong loss of evenness, despite the selective pressure of acetate amendment and change in the terminal electron accepting process. However, community membership was altered significantly. The method allows for sensitive, accurate profiling of the "long tail" of low abundance organisms that exist in many microbial communities, and can resolve population dynamics in response to environmental change.

59 BASIC BIOLOGICAL SCIENCES↗

Numerical simulations of vertical growth of hydraulic fractures and brine migration in geological formations above the Marcellus shale

One of the critical environmental questions about hydraulic fracturing in shales is the potential for contamination of ground and surface water. There are two specific concerns arising from hydraulic treatments: 1) whether hydraulic fractures extend upward through overlying strata to reach overlying aquifers containing drinking water, and 2) whether injected fluids push native fluids upward into these overlying aquifers. In this work, the extent of likely fracture growth through overlying layers during hydraulic treatment of the Marcellus shale was estimated using a hydraulic fracture model. A wide range of material and fluid flow properties in a multi-layered geologic model was considered. The model was based on conditions and characteristics applicable to the Marcellus shale in that part of the Appalachian basin within southwestern Pennsylvania. Predictions of vertical termination frequencies for hydraulic fractures were used in a multi-layer model of the strata and natural fractures for studying brine migration through the natural and induced fracture network. NFFLOW, the software for explicitly modeling flow within networks of fractures, was utilized to compute transient flow rates according to the schedule of injected fluid during hydraulic fracturing. To aid our analysis, the modeled sequence of geologic strata was capped with a fictitious unfractured, but moderately-permeable layer, which serves as a monitoring zone. The analysis assumes one well lateral was placed in the middle of the Marcellus shale with hydraulic fractures penetrating layers in the model. The newly-developed geomechanical module within NFFLOW was used to represent stress-sensitivity of the fractures. This allows the opening and closing of fracture apertures with changes in fluid pressures within fracture segments. Pressure increases in the formations overlying the Tully limestone, indicating fluid flow, was observed due to the hydraulic stimulation; and the impact of these increased pressures on brine migration towards the surface was considered.

03 NATURAL GAS↗

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.

alternative polyadenylation↗

Technical Report on Modeling for Quasispecies Abundance Inference with Confidence Intervals from Metagenomic Sequence Data

The overall aim of this project is to develop a software package, called MetaQuant, that can determine the constituents of a complex microbial sample and estimate their relative abundances by analysis of metagenomic sequencing data. The goal for Task 1 is to create a generative model describing the stochastic process underlying the creation of sequence read pairs in the data set. The stages in this generative process include the selection of a source genome sequence for each read pair, with probability dependent on its abundance in the sample. The other stages describe the evolution of the source genome from its nearest common ancestor with a reference genome, breakage of the source DNA into short fragments, and the errors in sequencing the ends of the fragments to produce read pairs.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of the myosins encoded in the recently completed Arabidopsis thaliana genome sequence

BACKGROUND: Three types of molecular motors play an important role in the organization, dynamics and transport processes associated with the cytoskeleton. The myosin family of molecular motors move cargo on actin filaments, whereas kinesin and dynein motors move cargo along microtubules. These motors have been highly characterized in non-plant systems and information is becoming available about plant motors. The actin cytoskeleton in plants has been shown to be involved in processes such as transportation, signaling, cell division, cytoplasmic streaming and morphogenesis. The role of myosin in these processes has been established in a few cases but many questions remain to be answered about the number, types and roles of myosins in plants. RESULTS: Using the motor domain of an Arabidopsis myosin we identified 17 myosin sequences in the Arabidopsis genome. Phylogenetic analysis of the Arabidopsis myosins with non-plant and plant myosins revealed that all the Arabidopsis myosins and other plant myosins fall into two groups - class VIII and class XI. These groups contain exclusively plant or algal myosins with no animal or fungal myosins. Exon/intron data suggest that the myosins are highly conserved and that some may be a result of gene duplication. CONCLUSIONS: Plant myosins are unlike myosins from any other organisms except algae. As a percentage of the total gene number, the number of myosins is small overall in Arabidopsis compared with the other sequenced eukaryotic genomes. There are, however, a large number of class XI myosins. The function of each myosin has yet to be determined.

NASA Discipline Plant Biology↗

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Sequencing at sea: challenges and experiences in Ion Torrent PGM sequencing during the 2013 Southern Line Islands Research Expedition

Genomics and metagenomics have revolutionized our understanding of marine microbial ecology and the importance of microbes in global geochemical cycles. However, the process of DNA sequencing has always been an abstract extension of the research expedition, completed once the samples were returned to the laboratory. During the 2013 Southern Line Islands Research Expedition, we started the first effort to bring next generation sequencing to some of the most remote locations on our planet. We successfully sequenced twenty six marine microbial genomes, and two marine microbial metagenomes using the Ion Torrent PGM platform on the Merchant Yacht Hanse Explorer. Onboard sequence assembly, annotation, and analysis enabled us to investigate the role of the microbes in the coral reef ecology of these islands and atolls. This analysis identified phosphonate as an important phosphorous source for microbes growing in the Line Islands and reinforced the importance of Lserine in marine microbial ecosystems. Sequencing in the field allowed us to propose hypotheses and conduct experiments and further sampling based on the sequences generated. By eliminating the delay between sampling and sequencing, we enhanced the productivity of the research expedition. By overcoming the hurdles associated with sequencing on a boat in the middle of the Pacific Ocean we proved the flexibility of the sequencing, annotation, and analysis pipelines.

59 BASIC BIOLOGICAL SCIENCES↗

Clinical, Molecular, and Computational Analysis in two cases with mitochondrial encephalomyopathy associated with SUCLG1 mutation in a consanguineous family

Highlights: • We reported 2 patients with mitochondrial encepahlomyopathy disorder. • The homoplasmic mutation c.41T > C (p.M > 14) in SUCLG1 gene was reported. • A decrease of mtDNA copy number in both patients was detected. • These are the first cases of mitochondrial encepahlomyopathy with succinyl coA deficiency reported in Africa. Deficiency of the mitochondrial enzyme succinyl COA ligase (SUCL) is associated with encephalomyopathic mtDNA depletion syndrome and methylmalonic aciduria. This disorder is caused by mutations in both SUCL subunits genes: SUCLG1 (α subnit) and SUCLA2 (β subnit). We report here, two Tunisian patients belonging to a consanguineous family with mitochondrial encephalomyopathy, hearing loss, lactic acidosis, hypotonia, psychomotor retardation and methylmalonic aciduria. Mutational analysis of SUCLG1 gene showed, for the first time, the presence of c.41T > C in the exon 1 at homozygous state. In-silico analysis revealed that this mutation substitutes a conserved methionine residue to a threonine at position 14 (p.M14T) located at the SUCLG1 protein mitochondrial targeting sequence. Moreover, these analysis predicted that this mutation alter stability structure and mitochondrial translocation of the protein. In Addition, a decrease in mtDNA copy number was revealed by real time PCR in the peripheral blood leukocytes in the two patients compared with controls.

60 APPLIED LIFE SCIENCES↗

Assessment of Critical Experiment Benchmark Applicability to a Large-Capacity HALEU Transportation Package Concept

This work presents an assessment of the applicability of existing benchmark critical experiments to the criticality safety code validation for a large-capacity high-assay low-enriched uranium (HALEU) transportation package concept. Numerous next-generation nuclear reactor designs require HALEU fuel, which is characterized by an enrichment between 5 and 20 wt% 235 U. The U.S. Department of Energy (DOE) has proposed to recover and downblend highly enriched uranium from DOE-owned used nuclear fuel to accelerate the demonstration of commercially viable microreactor technologies. One element of the infrastructure needed to demonstrate HALEU-fueled reactors is the ability to safely transport enriched product to be used for fuel fabrication. There is uncertainty as to whether existing critical benchmark experiment data are sufficient to support criticality safety code validation for HALEU transportation applications. The anticipated chemical form of the HALEU in the proposed transportation concept is UO 2 with 20 wt% 235 U/U. The concept uses a combination of an existing transportation packaging design and a novel basket design, including borated aluminum flux traps. The basket provides space for 18 reusable, stainless steel canisters that contain the HALEU. In 10 CFR 71, normal conditions of transport (NCTs) and hypothetical accident conditions (HACs) are defined for fissile material transportation packages. NCT and HAC KENO-VI models of the transportation package were developed using the Standardized Computer Analyses for Licensing Evaluation (SCALE) 6.2.3 computer code package, and optimum moderation conditions were determined using the SCALE SAMPLER sequence. The SCALE Tools for Sensitivity and Uncertainty Analysis Methodology Implementation (TSUNAMI) sequences were then used to compare the neutronic characteristics of 1584 International Criticality Safety Benchmark Evaluation Project benchmark critical experiments with the NCT and HAC HALEU transportation models. The TSUNAMI integral correlation coefficient c k was the criterion used to rank neutronic similarity. Thirty-four experiments were identified as similar (c k ≥ 0.9) to the NCT model, and 55 experiments were identified as similar to the HAC model. Hundreds of experiments were also identified as at least marginally similar (c k ≥ 0.8) to both models. The results indicate that additional critical experiments are unlikely to be needed to support HALEU transportation criticality safety analyses for package concepts similar to the concept package analyzed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Automated Identification of Nucleotide Sequences

STITCH is a computer program that processes raw nucleotide-sequence data to automatically remove unwanted vector information, perform reverse-complement comparison, stitch shorter sequences together to make longer ones to which the shorter ones presumably belong, and search against the user s choice of private and Internet-accessible public 16S rRNA databases. ["16S rRNA" denotes a ribosomal ribonucleic acid (rRNA) sequence that is common to all organisms.] In STITCH, a template 16S rRNA sequence is used to position forward and reverse reads. STITCH then automatically searches known 16S rRNA sequences in the user s chosen database(s) to find the sequence most similar to (the sequence that lies at the smallest edit distance from) each spliced sequence. The result of processing by STITCH is the identification of the most similar well-described bacterium. Whereas previously commercially available software for analyzing genetic sequences operates on one sequence at a time, STITCH can manipulate multiple sequences simultaneously to perform the aforementioned operations. A typical analysis of several dozen sequences (length of the order of 103 base pairs) by use of STITCH is completed in a few minutes, whereas such an analysis performed by use of prior software takes hours or days.

Osman, Shariff↗

Effects of different cultivation media on root bacterial community characteristics of greenhouse tomatoes

Tomato, as a typical greenhouse crop, is commonly first planted as seedlings in a variety of substrates before being transplanted into soil. However, there is rare research on the characteristics of the bacterial community in tomato roots under this planting mode. In this study, tomatoes were planted in pots containing three different cultivation media, including soil and two types of substrates in a greenhouse, followed by a transplanting treatment. After collecting tomato root samples, high-throughput sequencing and bioinformatic analysis were used to compare the differences in bacterial diversity and functions between tomato roots before and after transplanting in different cultivation media. In total, 702776 sequences were obtained, and the OTUs were belonging to 109 genera, 58 families, 41 orders, 14 classes, and 12 phyla. Among the three cultivation media, the β -diversity was significant, and there was a slight difference in bacterial species diversity along with a large difference in their abundance at the genus level. Soil and both substrates had 79 bacterial genera in common, these genera accounted for 68.70%, 76.70%, and 71.17% of the total genera found in the soil, substrate 1, and substrate 2, respectively. After being transplanted from the two substrates to the soil, the bacterial community structure and abundance exhibited similarities with those found in the soil. Furthermore, based on microbial function prediction, the microbial communities in the two-substrate environment demonstrated a greater potential for promoting growth, while the microbial communities in the soil exhibited a greater tendency to exert their antibacterial potential. Our findings offer theoretical support for the creation of artificially reconstructed microbial communities in greenhouse cultivation.

Zhang, Xinjian↗

Prochlorococcus phage ferredoxin: structural characterization and electron transfer to cyanobacterial sulfite reductases

Marine cyanobacteria are infected by phages whose genomes encode ferredoxin (Fd) electron carriers. These Fds are thought to redirect the energy harvested from light to phage-encoded oxidoreductases that enhance viral fitness, but it is unclear how the biophysical properties and partner specificities of phage Fds relate to those of photosynthetic organisms. Here, results of a bioinformatics analysis using a sequence similarity network revealed that phage Fds are most closely related to cyanobacterial Fds that transfer electrons from photosystems to oxidoreductases involved in nutrient assimilation. Structural analysis of myovirus P-SSM2 Fd (pssm2-Fd), which infects the cyanobacterium Prochlorococcus marinus, revealed high levels of similarity to cyanobacterial Fds (root mean square deviations of ≤0.5 Å). Additionally, pssm2-Fd exhibited a low midpoint reduction potential (–336 mV versus a standard hydrogen electrode), similar to other photosynthetic Fds, although it had lower thermostability (Tm = 28 °C) than did many other Fds. When expressed in an Escherichia coli strain deficient in sulfite assimilation, pssm2-Fd complemented bacterial growth when coexpressed with a P. marinus sulfite reductase, revealing that pssm2-Fd can transfer electrons to a host protein involved in nutrient assimilation. The high levels of structural similarity with cyanobacterial Fds and reactivity with a host sulfite reductase suggest that phage Fds evolved to transfer electrons to cyanobacterially encoded oxidoreductases.

cyanobacteria↗

The promising role of proteomes and metabolomes in defining the single-cell landscapes of plants

The plant community has a strong track-record of RNA sequencing technology deployment, which combined with the recent advent of spatial platforms (e.g., 10x genomics), has resulted in an explosion of outstanding single cell and nuclei datasets that can be put in an in situ context within tissues (e.g., a cell atlas)1. In the genomics era, application of proteomics technologies in the plant sciences has always trailed behind that of RNA sequencing technologies, largely due to accessibility, ease-of-use and access to expertise along with depth of analysis benefits. On the other hand, the use of early analytical tools for characterizing small molecules (metabolites) from plant systems predates nucleic acid sequencing and proteomics analysis2, as the search for plant-based natural products has played a significant role in improving human health throughout history. However, the employment of proteomics and metabolomics assays for characterizing plant cell processes now remains significantly behind transcriptional approaches, even though both provide a direct functional readout of cell states and phenotypes.

Anderton, Christopher R. [BATTELLE (PACIFIC NW LAB↗

pnnl/metagenomic_analysis

Metagenomic analysis is a computationally intensive task, this problem is compounded when many hundreds of samples need to be analyzed. We needed a program that would process the metagenomic sequencing data, manage computing resources, and combine the results into a unified output. We developed a script that prepares raw metagenomic sequencing data for analysis, configures analysis parameters, launches analysis jobs, keeps track of analysis progress, and combines analysis output. The analysis itself itself is done by the functionality provided by mothur (https://www.mothur.org/).

Avila, Andrew↗