Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genomic methods”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

Long-reads↗
DOE OSTI2021

Novel Gloeobacterales spp. from Diverse Environments across the Globe

Photosynthetic Cyanobacteria and their descendants are the only known organisms capable of oxygenic photosynthesis. Their metabolism permanently changed the Earth’s surface and the evolutionary trajectory of life, but little is known about their evolutionary history. Genomes of the Gloeobacterales, an order of deeply divergent photosynthetic Cyanobacteria, may hold clues about the evolutionary process. However, there are only three published genomes within this order, and it is difficult to make broad inferences based on such little data. Here, I describe five species within the Gloeobacterales retrieved from publicly available databases and examine their photosynthetic gene content and the environments in which Gloeobacterales genomes and 16S rRNA gene sequences are found. The Gloeobacterales contain reduced photosystems and inhabit cold, wet-rock, and low-light environments. They are likely present in low abundances due to their low growth rate. Future searches for Gloeobacterales should target these environments, and samples should be deeply sequenced to capture the low-abundance taxa. Publicly available databases contain undescribed taxa within the Gloeobacterales. However, searching through all available data with current methods is computationally expensive. Therefore, new methods must be developed to search for these and other evolutionarily important taxa. Once identified, these novel photosynthetic Cyanobacteria will help illuminate the origin and evolution of oxygenic photosynthesis. Early branching photosynthetic Cyanobacteria such as the Gloeobacterales may provide clues into the evolutionary history of oxygenic photosynthesis, but there are few genomes or cultured taxa from this order. Five new metagenome-assembled genomes suggest that members of the Gloeobacterales all contain reduced photosystems and lack genes associated with thylakoids and circadian rhythms. Their distribution suggests that they may thrive in environments that are marginal for other species, including wet-rock and cold environments. These traits may aid in the discovery and cultivation of novel species in this clade.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Natural transformation as a tool in Acinetobacter baylyi : Evolution by amplification of gene copy number

For many years, the natural competency of Acinetobacter baylyi ADP1 facilitated studies of bacterial metabolism, biochemistry, and physiology. With the advent of synthetic biology, new opportunities arise to exploit the remarkable transformability and chromosomal plasticity of this model organism. In this chapter, we describe a recently developed method, “Evolution by Amplification and Synthetic Biology” (EASy). EASy allows the targeted amplification of chromosomal segments that give rise to new phenotypes. Increased gene dosage regulates protein expression in a rudimentary fashion by establishing a chromosomal array in which copy number adjusts via recombination between repeated DNA sequences. Selective conditions enrich for cells within the population that confer a growth advantage. Under continuous selective pressure, beneficial mutations may accumulate in any genomic region. Such mutations favor decreases in the average copy number of the target region. Thus, the genetic flexibility afforded by transient copy number variation helps accelerate the selection of engineered strains with desired traits during laboratory evolution. As a result, thanks to the extremely simple genetic manipulation of Acinetobacter baylyi ADP1, the EASy method can be readily implemented by researchers without the need for advanced instrumentation or complex cloning techniques.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Hypergraph Models of Biological Networks to Identify Genes Critical to Pathogenic Viral Response

Motivation: Representing biological networks as graphs is a powerful approach to reveal underlying patterns, signatures, and critical components from high-throughput biomolecular data. However, graphs do not natively capture the multi-way relationships present among genes and proteins in biological systems such as protein complexes, metabolic reactions, and signal transduction pathways. Hypergraphs are generalizations of graphs that naturally model multi-way interactions in data, and we therefore seek to understand how they can more faithfully identify, and potentially predict, complex relationships in genomic expression data sets. Results: We compiled a novel data set of transcriptional host response to pathogenic viral infections and formulated relationships between genes as a hypergraph where hyperedges are differentially expressed genes and vertices represent conditions. We find that hypergraph betweenness centrality is a superior method for identification of genes important to viral response when compared with graph centrality. Our results demonstrate the utility of using hypergraphs to represent complex biological systems, and highlight potentially interesting biological results about host response to highly pathogenic viruses.

systems biology, hypergraph, viral infection, biol↗
DOE OSTI2024

The Defects Genome of Janus Transition Metal Dichalcogenides

2D Janus Transition Metal Dichalcogenides (TMDs) have attracted much interest due to their exciting quantum properties arising from their unique two-faced structure, broken-mirror symmetry, and consequent colossal polarization field within the monolayer. While efforts are made to achieve high-quality Janus monolayers, the existing methods rely on highly energetic processes that introduce unwanted grain-boundary and point defects with still unexplored effects on the material's structural and excitonic properties Through high-resolution scanning transmission electron microscopy (HRSTEM), density functional theory (DFT), and optical spectroscopy measurements; this work introduces the most encountered and energetically stable point defects. It establishes their impact on the material's optical properties. HRSTEM studies show that the most energetically stable point defects are single (V S and V Se ) and double chalcogen vacancy (V S −V Se ), interstitial defects (Mi), and metal impurities (MW) and establish their structural characteristics. DFT further establishes their formation energies and related localized bands within the forbidden band. Cryogenic excitonic studies on h-BN-encapsulated Janus monolayers offer a clear correlation between these structural defects and observed emission features, which closely align with the results of the theory. Finally, the overall results introduce the defect genome of Janus TMDs as an essential guideline for assessing their structural quality and device properties.

2D materials↗
DOE OSTI2024

Machine learning and deep learning tools for the automated capture of cancer surveillance data

The National Cancer Institute and the Department of Energy strategic partnership applies advanced computing and predictive machine learning and deep learning models to automate the capture of information from unstructured clinical text for inclusion in cancer registries. Applications include extraction of key data elements from pathology reports, determination of whether a pathology or radiology report is related to cancer, extraction of relevant biomarker information, and identification of recurrence. With the growing complexity of cancer diagnosis and treatment, capturing essential information with purely manual methods is increasingly difficult. These new methods for applying advanced computational capabilities to automate data extraction represent an opportunity to close critical information gaps and create a nimble, flexible platform on which new information sources, such as genomics, can be added. This will ultimately provide a deeper understanding of the drivers of cancer and outcomes in the population and increase the timeliness of reporting. These advances will enable better understanding of how real-world patients are treated and the outcomes associated with those treatments in the context of our complex medical and social environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A rapid thioacidolysis method for biomass lignin composition and tricin analysis

Abstract Background Biomass composition varies from plant to plant and greatly affects biomass utilization. Lignin is a heterogeneous phenolic polymer derived mainly from p -coumaryl, coniferyl, and sinapyl alcohols and makes up to 10–25% of lignocellulosic biomass. Recently, tricin, an O -methylated flavone, was identified as a lignin monomer in many grass species. Tricin may function as a nucleation site for lignification and is advocated as a novel target for lignin engineering to reduce lignin content and improve biomass digestibility in grasses. Thioacidolysis is an analytical method that can be adapted to analyze both lignin monomeric composition and tricin content in the lignin polymer. However, the original thioacidolysis procedure is complex, laborious, and time consuming, making it difficult to be adopted for large-scale screening in biomass research. In this study, a modified, rapid higher throughput thioacidolysis method was developed. Results In combination with gas chromatography–mass spectrometry (GC–MS) and liquid chromatography–mass spectrometry (LC–MS), the modified thioacidolysis method can be used to simultaneously characterize the lignin composition and tricin content using 2–5 mg of dry samples. The modified method eliminates the solvent extraction and drastically improves the throughput; 80 samples can be processed in one day per person. Our results indicate that there is no significant difference in the determination of lignin S/G ratio and tricin content between the original and modified methods. Conclusions A modified thioacidolysis protocol was established. The results demonstrate that the modified method can be used for rapid, high-throughput, and reliable lignin composition and tricin content analyses for screening transgenic plants for cell wall modifications or in large-scale genome-wide association studies (GWAS).

09 BIOMASS FUELS↗
DOE OSTI2020

Data from A Genetic Toolbox for Metabolic Engineering of Issatchenkia orientalis

The nonconventional yeast Issatchenkia orientalis can grow under highly acidic conditions and has been explored for production of various organic acids. However, its broader application is hampered by the lack of efficient genetic tools to enable sophisticated metabolic manipulations. We recently constructed an episomal plasmid based on the autonomously replicating sequence (ARS) from Saccharomyces cerevisiae (ScARS) in I. orientalis and developed a CRISPR/Cas9 system for multiplexed gene deletions. Here we report three additional genetic tools including: (1) identification of a 0.8 kb centromere-like (CEN-L) sequence from the I. orientalis genome by using bioinformatics and functional screening; (2) discovery and characterization of a set of constitutive promoters and terminators under different culture conditions by using RNA-Seq analysis and a fluorescent reporter; and (3) development of a rapid and efficient in vivo DNA assembly method in I. orientalis , which exhibited ~100% fidelity when assembling a 7 kb-plasmid from seven DNA fragments ranging from 0.7 kb to 1.7 kb. As proof of concept, we used these genetic tools to rapidly construct a functional xylose utilization pathway in I. orientalis .

Conversion↗
DOE OSTI2023

Detecting macroevolutionary genotype–phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype–phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω C , which measures the error-corrected convergence rate of protein evolution. While ω C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype–phenotype associations, even in lineages that diverged for hundreds of millions of years.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Distant-acting Enhancers in Development and Disease

Multiple converging lines of evidence from experimental and human genetic studies indicate that both common and rare sequence variants involving enhancers play major roles in Mendelian and complex human disease phenotypes. We use a combination of sequence-based molecular approaches including chromatin mapping and synthesis-enabled modification of enhancer sequences, large-scale transgenic mouse studies (http://enhancer.lbl.gov), and CRISPR genome editing in the mouse model to study the in vivofunction of enhancers in developmental, evolutionary, and disease-related processes. Using examples from our ongoing work, I will illustrate how these methods provide insight into the function and evolution of distant-acting regulatory sequences and offer a starting point for understanding their role in human disease. I will discuss how improved in vivo mouse assays can be used to study candidate pathogenic human mutations in enhancers at scale, revealing the extent to which human genetic studies alone can be insufficient to identify truly causal mutations in enhancers. I will also show how such assays can be used to explore the sensitivity of enhancers that are extremely conserved ("ultraconserved") in evolution to mutations, revealing a remarkable functional robustness to sequence changes that contrasts with apparent evolutionary constraints.

Visel, A↗
DOE OSTI2025

Data imbalance in drug response prediction: multi-objective optimization approach in deep learning setting

Abstract Drug response prediction (DRP) methods tackle the complex task of associating the effectiveness of small molecules with the specific genetic makeup of the patient. Anti-cancer DRP is a particularly challenging task requiring costly experiments as underlying pathogenic mechanisms are broad and associated with multiple genomic pathways. The scientific community has exerted significant efforts to generate public drug screening datasets, giving a path to various machine learning models that attempt to reason over complex data space of small compounds and biological characteristics of tumors. However, the data depth is still lacking compared to application domains like computer vision or natural language processing domains, limiting current learning capabilities. To combat this issue and improves the generalizability of the DRP models, we are exploring strategies that explicitly address the imbalance in the DRP datasets. We reframe the problem as a multi-objective optimization across multiple drugs to maximize deep learning model performance. We implement this approach by constructing Multi-Objective Optimization Regularized by Loss Entropy loss function and plugging it into a Deep Learning model. We demonstrate the utility of proposed drug discovery methods and make suggestions for further potential application of the work to achieve desirable outcomes in the healthcare field.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Populus_trichocarpa_Breeding_Population_SNPs

These data are from the manuscript “Application of Genomic Prediction in a Populus trichocarpa Breeding Program”, by Brian J. Stanton, David Macaya-Sanz, Chanaka Roshan Abeyratne, David Kainer, Kathy Haiby, Austin Himes, Carlos Gantz, Gerald A. Tuskan, and Stephen P. DiFazio. The data are based on genome resequencing to approximately 10X depth on two collections of Populus trichocarpa trees from Oregon, Washington, California, and British Columbia. The first collection consists of 293 genets collected by Poplar Innovations LLC for a breeding program. The second collection consists of 961 trees collected for the purpose of genome-wide association studies. These genets were sequenced using short, paired-end Illumina sequence reads (Chhetri et al. 2019). Reads were aligned to the P. trichocarpa ′Stettler-14′ reference (Hofmeister et al. 2020), with minor modifications to correct mis-assemblies (Zhou et al. 2020), and variants were called as per methods described in (Abeyratne et al. 2023). Identified variants were filtered using GATK’s VariantFiltration tool (DePristo et al. 2011), with filter expression flag set to “AF < 0.01 || AF > 0.99 || QD < 10.0 || ExcessHet > 20.0 || FS > 10.0 || MQ < 58.0”. SNPs with severe departures from Hardy−Weinberg expectations (exact-test p< 0.01) were also removed using vcftools --hwe flag (Danecek et al. 2011), resulting in 15,627,211 bi-allelic SNPs. The data included here consist of 141,903 high quality bi-allelic genome-wide SNPs obtained by further filtering the original SNP dataset using vcftools with flags --maf 0.05, --max-maf 0.95, --max-missing 0.95, --min-meanDP 10.75, --max-meanDP 43.00, --thin 2000. Collectively, these filtering parameters removed SNPs with 1) a minor allele frequency ≤ 0.05; 2) proportion of missing data for individual loci exceeding 5%; 3) sequencing depth more than 2X mean-depth or less than 0.5X mean-depth; or 4) a distance of

09 BIOMASS FUELS↗
DOE OSTI2023

Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids

This dataset contains all data used in the paper "Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids". The dataset includes genotypes and phenotypic data from two autotetraploid species Miscanthus sacchariflorus and Vaccinium corymbosum that was used used for genome wide association studies and genomic prediction and the scripts used in the analysis. In this V2, 2 files have the raw data are added: "Miscanthus_sacchariflorus_RADSeq.vcf" is the VCF file with the raw SNP calls of the Miscanthus sacchariflorus data used for genotype calling using the 6 genotype calling methods. "Blueberry_data_read_depths.RData" is the a RData file with the read depth data that was used for genotype calling in the Blueberry dataset.

allelic dosage↗
DOE OSTI2026

wavess 1.2: presenting an HLA-aware within-host virus sequence simulation framework

Motivation Understanding how virus sequences are shaped by selection can inform vaccine design and transmission inference. Modeling within-host evolution to interrogate these questions requires a detailed mechanistic framework that accurately captures sequence diversification. The CD8 + cytotoxic T-lymphocyte (CTL) response plays an important role in immune-mediated selection and can leave strong signatures in virus sequences; however, existing sequence-based within-host virus modeling frameworks do not explicitly include a human leukocyte antigen (HLA)-aware CTL response. Results We extended our previously published within-host sequence evolution simulator, wavess, to include an explicit CTL response, and share a method for identifying HLA-specific CTL epitopes given a founder virus sequence. We also updated the model to permit a variable recombination rate, which allows for modeling non-adjacent genes, segmented genomes, and recombination hotspots. These extensions to wavess allow for more accurate simulation of viruses and virus genes, particularly in regions of the genome where the immune response is dominated by CTLs (rather than antibodies). It also provides the foundation for investigations of how these newly-added biological mechanisms influence within-host evolution. Availability and implementation The core of wavess is written in Python 3, with helper functions written in R. It is available at https://github.com/MolEvolEpid/wavess.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Machine learning approaches for integrating multi-omics data to expand microbiome annotation

Preliminary: This final report corresponds to a grant (DE-SC0021216) that was awarded to the University of Montana. Mid-way through the grant period, I relocated from the University of Montana to the University of Arizona. The grant was ended at University of Montana in late 2022, with all efforts concluding on 08/26/22; the remaining funds supporting the project were relinquished by University of Montana, and were later awarded to University of Arizona under a new grant, with start date 04/01/23. This report focuses on results of research efforts at UMontana through 08/26/22. Results: We made progress in each of the three aims of the proposal. We released software that identifies and fills gaps in the annotation of metabolic proteins within bacterial genomes. We made substantial progress in developing software for alignment-based annotation of protein coding DNA, allowing for coding frameshifts caused by sequencing error. Finally, we made notable progress in developing AI methods (specifically: a neural embedding model) for identifying similarities between protein sequences based on amino-wise latent vectors. These efforts were supplemented by development of methods for protein modeling in support of predicting protein-drug binding activity, and by my leadership of a team in the NIH/DOE 2021 Petabyte-Scale Sequence Search hack-a-thon.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Development of male-sterile lines of Setaria viridis to accelerate C 4 model plant genetics

Setaria viridis is a diploid C 4 grass in the Poaceae family, notable for its rapid life cycle of 6–8 weeks from sowing to seed—much shorter than the 4–5 months required by crops such as Zea mays and Sorghum bicolor . This fast growth makes S. viridis a valuable model for C 4 crop research. Genetic crosses are essential for studying gene function, but manual crossing is labor-intensive and time-consuming. Here, to address this, we developed a male-sterile line by targeting the S. viridis ortholog of Setaria italica NO POLLEN 1 ( SiNP1 ), which encodes a glucose–methanol–choline oxidoreductase required for pollen exine formation. Using Cas9 and TREX2 -mediated genome editing, we generated SiNP1 knockouts in both the S. viridis ME034V and A10.1 backgrounds that were fully male-sterile. Backcrossing T 0 male-sterile plants to ME034V wild-type followed by selfing yielded a stable BC 1 F 2 line homozygous for a 59 bp deletion in the S. viridis NO POLLEN 1 gene, easily genotyped by PCR and maintained by heterozygous siblings. Using this line, we developed a simple and efficient crossing protocol that eliminates the need for emasculation. This method enables a single person to perform up to 100 crosses per day—compared to 15 using traditional methods—and yields 20–32 F 1 hybrid seeds per panicle with 100% genetic purity. We also quantified pollen flow and outcrossing frequencies under greenhouse conditions to develop optimal bagging strategies and prevent unintended pollination. This resource accelerates genetic research in S. viridis , enhancing its utility as a premier C 4 model for mapping and functional genomics.

C4 research↗
DOE OSTI2023

Topological structure of complex predictions

Abstract Current complex prediction models are the result of fitting deep neural networks, graph convolutional networks or transducers to a set of training data. A key challenge with these models is that they are highly parameterized, which makes describing and interpreting the prediction strategies difficult. We use topological data analysis to transform these complex prediction models into a simplified topological view of the prediction landscape. The result is a map of the predictions that enables inspection of the model results with more specificity than dimensionality-reduction methods such as tSNE and UMAP. The methods scale up to large datasets across different domains. We present a case study of a transformer-based model previously designed to predict expression levels of a piece of DNA in thousands of genomic tracks. When the model is used to study mutations in the BRCA1 gene, our topological analysis shows that it is sensitive to the location of a mutation and the exon structure of BRCA1 in ways that cannot be found with tools based on dimensionality reduction. Moreover, the topological framework offers multiple ways to inspect results, including an error estimate that is more accurate than model uncertainty. Further studies show how these ideas produce useful results in graph-based learning and image classification.

Computer Science↗
DOE OSTI2024

Enhancers in Plant Development, Adaptation and Evolution

Understanding plant responses to developmental and environmental cues is crucial for studying morphological divergence and local adaptation. Gene expression changes, governed by cis-regulatory modules (CRMs) including enhancers, are a major source of plant phenotypic variation. However, while genome-wide approaches have revealed thousands of putative enhancers in mammals, far fewer have been identified and functionally characterized in plants. This review provides an overview of how enhancers function to control gene regulation, methods to predict DNA sequences that may have enhancer activity, methods utilized to functionally validate enhancers and the current knowledge of enhancers in plants, including how they impact plant development, response to environment and evolutionary adaptation.

59 BASIC BIOLOGICAL SCIENCES↗