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.

SourceResearch connectionFrom yearThrough year

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

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Machine learning analysis of RB-TnSeq fitness data predicts functional gene modules in Pseudomonas putida KT2440

ABSTRACT There is growing interest in engineering Pseudomonas putida KT2440 as a microbial chassis for the conversion of renewable and waste-based feedstocks, and metabolic engineering of P. putida relies on the understanding of the functional relationships between genes. In this work, independent component analysis (ICA) was applied to a compendium of existing fitness data from randomly barcoded transposon insertion sequencing (RB-TnSeq) of P. putida KT2440 grown in 179 unique experimental conditions. ICA identified 84 independent groups of genes, which we call fModules (“functional modules”), where gene members displayed shared functional influence in a specific cellular process. This machine learning-based approach both successfully recapitulated previously characterized functional relationships and established hitherto unknown associations between genes. Selected gene members from fModules for hydroxycinnamate metabolism and stress resistance, acetyl coenzyme A assimilation, and nitrogen metabolism were validated with engineered mutants of P. putida . Additionally, functional gene clusters from ICA of RB-TnSeq data sets were compared with regulatory gene clusters from prior ICA of RNAseq data sets to draw connections between gene regulation and function. Because ICA profiles the functional role of several distinct gene networks simultaneously, it can reduce the time required to annotate gene function relative to manual curation of RB-TnSeq data sets. IMPORTANCE This study demonstrates a rapid, automated approach for elucidating functional modules within complex genetic networks. While Pseudomonas putida randomly barcoded transposon insertion sequencing data were used as a proof of concept, this approach is applicable to any organism with existing functional genomics data sets and may serve as a useful tool for many valuable applications, such as guiding metabolic engineering efforts in other microbes or understanding functional relationships between virulence-associated genes in pathogenic microbes. Furthermore, this work demonstrates that comparison of data obtained from independent component analysis of transcriptomics and gene fitness datasets can elucidate regulatory-functional relationships between genes, which may have utility in a variety of applications, such as metabolic modeling, strain engineering, or identification of antimicrobial drug targets.

09 BIOMASS FUELS↗

DOE OSTI2016

The role of structural parameters in DNA cyclization

Abstract Background The intrinsic bendability of DNA plays an important role with relevance for myriad of essential cellular mechanisms. The flexibility of a DNA fragment can be experimentally and computationally examined by its propensity for cyclization, quantified by the Jacobson-Stockmayer J factor. In this study, we use a well-established coarse-grained three-dimensional model of DNA and seven distinct sets of experimentally and computationally derived conformational parameters of the double helix to evaluate the role of structural parameters in calculating DNA cyclization. Results We calculate the cyclization rates of 86 DNA sequences with previously measured J factors and lengths between 57 and 325 bp as well as of 20,000 randomly generated DNA sequences with lengths between 350 and 4000 bp. Our comparison with experimental data is complemented with analysis of simulated data. Conclusions Our data demonstrate that all sets of parameters yield very similar results for longer DNA fragments, regardless of the nucleotide sequence, which are in agreement with experimental measurements. However, for DNA fragments shorter than 100 bp, all sets of parameters performed poorly yielding results with several orders of magnitude difference from the experimental measurements. Our data show that DNA cyclization rates calculated using conformational parameters based on nucleosome packaging data are most similar to the experimental measurements. Overall, our study provides a comprehensive large-scale assessment of the role of structural parameters in calculating DNA cyclization rates.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2012

Livermore Metagenomics Analysis Toolkit

LMAT is designed to take as input a collection of raw metagenomic sequencer reads, and search each read against a reference genome database and assign a taxonomic label and confidence value to each read and report a summary of the predicted taxonomic contents of the metagenomic sample.

Source record↗

DOE OSTI2016

A novel analysis method for paired-sample microbial ecology experiments

Many microbial ecology experiments use sequencing data to measure a community s response to an experimental treatment. In a common experimental design, two units, one control and one experimental, are sampled before and after the treatment is applied to the experimental unit. The four resulting samples contain information about the dynamics of organisms that respond to the treatment, but there are no analytical methods designed to extract exactly this type of information from this configuration of samples. Here we present an analytical method specifically designed to visualize and generate hypotheses about microbial community dynamics in experiments that have paired samples and few or no replicates. The method is based on the Poisson lognormal distribution, long studied in macroecology, which we found accurately models the abundance distribution of taxa counts from 16S rRNA surveys. To demonstrate the method s validity and potential, we analyzed an experiment that measured the effect of crude oil on ocean microbial communities in microcosm. Our method identified known oil degraders as well as two clades, Maricurvus and Rhodobacteraceae, that responded to amendment with oil but do not include known oil degraders. Furthermore, our approach is sensitive to organisms that increased in abundance only in the experimental unit but less sensitive to organisms that increased in both control and experimental units, thus mitigating the role of bottle effects .

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2012

Livermore Metagenomics Analysis Toolkit

Source record↗

DOE OSTI2016

From DNA to FBA: How to Build Your Own Genome-Scale Metabolic Model

Microbiological studies are increasingly relying on in silico methods to perform exploration and rapid analysis of genomic data, and functional genomics studies are supplemented by the new perspectives that genome-scale metabolic models offer. A mathematical model consisting of a microbe’s entire metabolic map can be rapidly determined from whole-genome sequencing and annotating the genomic material encoded in its DNA. Flux-balance analysis (FBA), a linear programming technique that uses metabolic models to predict the phenotypic responses imposed by environmental elements and factors, is the leading method to simulate and manipulate cellular growth in silico. However, the process of creating an accurate model to use in FBA consists of a series of steps involving a multitude of connections between bioinformatics databases, enzyme resources, and metabolic pathways. We present the methodology and procedure to obtain a metabolic model using PyFBA, an extensible Pythonbased open-source software package aimed to provide a platform where functional annotations are used to build metabolic models (http://linsalrob.github.io/PyFBA). Backed by the Model SEED biochemistry database, PyFBA contains methods to reconstruct a microbe’s metabolic map, run FBA upon different media conditions, and gap-fill its metabolism. The extensibility of PyFBA facilitates novel techniques in creating accurate genome-scale metabolic models.

genome annotation↗

DOE OSTI2015

New signatures of the Milky Way formation in the local halo and inner-halo streamers in the era of Gaia

We explore the vicinity of the Milky Way through the use of spectrophotometric data from the Sloan Digital Sky Survey and high-quality proper motions derived from multi-epoch positions extracted from the Guide Star Catalog II database. In order to identify and characterize streams as relics of the Milky Way formation, we start with classifying, selecting, and studying 2417 subdwarfs with [Fe/H] < −1.5 up to 3 kpc away from the Sun as tracers of the local halo system. Then, through phase-space analysis, we find statistical evidence of five discrete kinematic overdensities among 67 of the fastest-moving stars and compare them to high-resolution N-body simulations of the interaction between a Milky Way–like galaxy and orbiting dwarf galaxies with four representative cases of merging histories. The observed overdensities can be interpreted as fossil substructures consisting of streamers torn from their progenitors; such progenitors appear to be satellites on prograde and retrograde orbits on different inclinations. In particular, of the five detected overdensities, two appear to be associated, yielding 21 additional main-sequence members, with the stream of Helmi et al. that our analysis confirms is on a high-inclination prograde orbit. The three newly identified kinematic groups could be associated with the retrograde streams detected by Dinescu and Kepley et al.; whatever their origin, the progenitor(s) would be on retrograde orbit(s) and inclination(s) within the range 10{sup ∘}÷60{sup ∘}. Finally, we use our simulations to investigate the impact of observational errors and compare the current picture to the promising prospect of highly improved data expected from the Gaia mission.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

HiFiAdapterFilt, a memory efficient read processing pipeline, prevents occurrence of adapter sequence in PacBio HiFi reads and their negative impacts on genome assembly

Abstract Background Pacific Biosciences HiFi read technology is currently the industry standard for high accuracy long-read sequencing that has been widely adopted by large sequencing and assembly initiatives for generation of de novo assemblies in non-model organisms. Though adapter contamination filtering is routine in traditional short-read analysis pipelines, it has not been widely adopted for HiFi workflows. Results Analysis of 55 publicly available HiFi datasets revealed that a read-sanitation step to remove sequence artifacts derived from PacBio library preparation from read pools is necessary as adapter sequences can be erroneously integrated into assemblies. Conclusions Here we describe the nature of adapter contaminated reads, their consequences in assembly, and present HiFiAdapterFilt, a simple and memory efficient solution for removing adapter contaminated reads prior to assembly.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Radioisotope Identification with List-Mode Gamma-Ray Data

This work explores the potential of utilizing temporal data from gamma-ray detectors, known as list-mode data, to enhance radioisotope identification. Traditional identification methods, which rely on full gamma-ray spectrum analysis, often require long dwell times and struggle with spectra containing similarly spaced spectral peaks. We hypothesize that by leveraging the probabilistic nature of nuclear decay and the time-encoded information from decay sequences and interactions with surrounding materials, we can improve classification accuracy over static spectral analysis. This research examines the temporal content of list-mode data through exploratory data analysis via correlation discovery and qualitative distribution analysis. Additionally, we propose a probabilistic classification model that can utilize spectral data, temporal data, or both to determine if the incorporation of temporal information improves radioisotope identification. Our findings suggest that the temporal information present in list-mode gamma-ray data has merit and should be further investigated to develop more robust and optimal methods for utilizing this temporal information in applications requiring radioisotope identification.

List-mode data↗

DOE OSTI2015

RADIAL VELOCITY OBSERVATIONS AND LIGHT CURVE NOISE MODELING CONFIRM THAT KEPLER-91b IS A GIANT PLANET ORBITING A GIANT STAR

Kepler-91b is a rare example of a transiting hot Jupiter around a red giant star, providing the possibility to study the formation and composition of hot Jupiters under different conditions compared to main-sequence stars. However, the planetary nature of Kepler-91b, which was confirmed using phase-curve variations by Lillo-Box et al., was recently called into question based on a re-analysis of Kepler data. We have obtained ground-based radial velocity observations from the Hobby-Eberly Telescope and unambiguously confirm the planetary nature of Kepler-91b by simultaneously modeling the Kepler and radial velocity data. The star exhibits temporally correlated noise due to stellar granulation which we model as a Gaussian Process. We hypothesize that it is this noise component that led previous studies to suspect Kepler-91b to be a false positive. Our work confirms the conclusions presented by Lillo-Box et al. that Kepler-91b is a 0.73 ± 0.13 M {sub Jup} planet orbiting a red giant star.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

IMG/M: integrated genome and metagenome comparative data analysis system

The Integrated Microbial Genomes with Microbiome Samples (IMG/M: https://img.jgi.doe.gov/m/) system contains annotated DNA and RNA sequence data of (i) archaeal, bacterial, eukaryotic and viral genomes from cultured organisms, (ii) single cell genomes (SCG) and genomes from metagenomes (GFM) from uncultured archaea, bacteria and viruses and (iii) metagenomes from environmental, host associated and engineered microbiome samples. Sequence data are generated by DOE's Joint Genome Institute (JGI), submitted by individual scientists, or collected from public sequence data archives. Structural and functional annotation is carried out by JGI's genome and metagenome annotation pipelines. A variety of analytical and visualization tools provide support for examining and comparing IMG/M's datasets. IMG/M allows open access interactive analysis of publicly available datasets, while manual curation, submission and access to private datasets and computationally intensive workspace-based analysis require login/password access to its expert review(ER) companion system (IMG/M ER: https://img.jgi.doe.gov/ mer/). Since the last report published in the 2014 NAR Database Issue, IMG/M's dataset content has tripled in terms of number of datasets and overall protein coding genes, while its analysis tools have been extended to cope with the rapid growth in the number and size of datasets handled by the system.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Ecogenomics of Groundwater Phages Suggests Niche Differentiation Linked to Specific Environmental Tolerance

To our knowledge, this is the first study to identify the bacteriophage distribution in a groundwater ecosystem shedding light on their prevalence and distribution across metal-contaminated and background sites. Our study is uniquely based on selective sequencing of solely the extrachromosomal elements of a microbiome followed by analysis for viral signatures, thus establishing a more focused approach for phage identifications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2013

From molecules to dynamic biological communities

Microbial ecology is flourishing, and in the process, is making contributions to how the ecology and biology of large organisms is understood. Ongoing advances in sequencing technology and computational methods have enabled the collection and analysis of vast amounts of molecular data from diverse biological communities. While early studies focused on cataloguing microbial biodiversity in environments ranging from simple marine ecosystems to complex soil ecologies, more recent research is concerned with community functions and their dynamics over time. Models and concepts from traditional ecology have been used to generate new insight into microbial communities, and novel system-level models developed to explain and predict microbial interactions. The process of moving from molecular inventories to functional understanding is complex and challenging, and never more so than when many thousands of dynamic interactions are the phenomena of interest. We outline the process of how epistemic transitions are made from producing catalogues of molecules to achieving functional and predictive insight, and show how those insights not only revolutionize what is known about biological systems but also about how to do biology itself. Examples will be drawn primarily from analyses of different human microbiota, which are the microbial consortia found in and on areas of the human body, and their associated microbiomes (the genes of those communities). Molecular knowledge of these microbiomes is transforming microbiological knowledge, as well as broader aspects of human biology, health and disease.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2016

27-Hydroxycholesterol increases Myc protein stability via suppressing PP2A, SCP1 and FBW7 transcription in MCF-7 breast cancer cells

27-hydroxycholesterol (27-HC), the most abundant metabolite of cholesterol, is a risk factor for breast cancer. It can increase the proliferation of breast cancer cells and promote the metastasis of breast tumours in mouse models. Myc is a critical oncoprotein overexpressed in breast cancer. However, whether 27-HC affects Myc expression has not been reported. In the current study, we aimed to investigate the effects of 27-HC on Myc and the underlying mechanisms in MCF-7 breast cancer cells. Our data demonstrated that 27-HC activated Myc via increasing its protein stability. Three key negative modulators of Myc protein stability, PP2A, SCP1 and FBW7, were suppressed by 27-HC at the transcriptional level. We performed a data-mining analysis of the chromatin immunoprecipitation with next-generation DNA sequencing (ChIP-Seq) data in the ChIPBase, and discovered that a number of putative transcription factors (TFs), including Myc itself, were involved in the transcriptional regulation of PP2A, SCP1 and FBW7. Our results provide a novel mechanistic insight into the activation of Myc by 27-HC via transcriptional repression of PP2A, SCP1 and FBW7 to increase Myc protein stability in breast cancer cells. - Highlights: • 27-Hydroxycholesterol (27-HC) activates Myc via increasing its protein stability. • 27-HC inhibits PP2A and SCP1 transcription to block pS62-Myc dephosphorylation. • 27-HC suppresses FBW7 transcription to prevent pT58-Myc degradation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

TOUGH-UDEC: A simulator for coupled multiphase fluid flows, heat transfers and discontinuous deformations in fractured porous media

A numerical simulator entitled TOUGH-UDEC is presented for the assessment of coupled thermal-hydraulic-mechanical processes in fractured porous media. Two existing well-established codes, TOUGH2 and UDEC, are coupled to model multiphase fluid flows, heat transfers, and discontinuous deformations in fractured porous media by means of discrete fracture representation. TOUGH2 is widely used for the modeling of heat transfers and multiphase multicomponent fluid flows, and UDEC is a well-known distinct element code for rock mechanics. The two codes are solved sequentially, with coupling parameters passed to each equation at specific intervals. After solving thermal-hydraulic equations within the TOUGH2 code, pressure and temperature information is imported into the UDEC model. Upon solving the mechanical equation within the UDEC code the calculated fracture apertures are converted to the equivalent permeability and porosity values for a TOUGH2 flow analysis. The solution is calculated by iteratively following an explicit sequence for numerical efficiency. Verifications are presented to demonstrate the capabilities of the coupled TOUGH-UDEC simulator. Three application examples of (1) shear dilation due to increased pore pressure, (2) thermal stress and (3) CO 2 injection, show that the new simulator can be an effective tool for geoengineering applications involving shear activation of fractures and faults.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Aftershocks of the 2012 Off-Coast of Sumatra Earthquake Sequence

Aftershocks of the 2012 Off-Coast of Sumatra Earthquake Sequence exhibit a complex and diffuse spatial distribution. The first-order complexity in aftershock distribution is clear and well beyond the influence of typical earthquake location uncertainty. The sequence included rupture of multiple fault segments, spatially separated. We use surface-wave based relative centroid locations to examine whether, at the small scale, the distribution of the aftershocks was influenced by location errors. Surface-wave based relative location has delineated precise oceanic transform fault earthquake locations in multiple regions. However, the relocated aftershocks off the coast of Sumatra seldom align along simple linear trends that are compatible with the corresponding fault strikes as estimated for the GCMT catalog. The relocation of roughly 60 moderate-earthquake epicentroids suggests that the faulting involved in the 2012 earthquake aftershock sequence included strain release along many short fault segments. Statistical analysis and temporal variations of aftershocks show a typical decay of the aftershocks but a relatively low number of aftershocks, as is common for intraplate oceanic earthquakes. Coulomb stress calculations indicate that most of the moderate-magnitude aftershocks are compatible with stress changes predicted by the large-event slip models. Furthermore the patterns in the aftershocks suggest that the formation of the boundary and eventual localization of deformation between the Indian and Australian plate is a complicated process.

2012 Sumatra Earthquake↗

DOE OSTI2023

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins—including gene expression, chromatin occupancy and three-dimensional conformation—we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

59 BASIC BIOLOGICAL SCIENCES↗