Search NASA⌕ Search

SEARCH · Search NASA

Results for “genomic prediction”

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 73 records · Page 4

Towards a robust out-of-the-box neural network model for genomic data

The accurate prediction of biological features from genomic data is paramount for precision medicine and sustainable agriculture. For decades, neural network models have been widely popular in fields like computer vision, astrophysics and targeted marketing given their prediction accuracy and their robust performance under big data settings. Yet neural network models have not made a successful transition into the medical and biological world due to the ubiquitous characteristics of biological data such as modest sample sizes, sparsity, and extreme heterogeneity. Here, we investigate the robustness, generalization potential and prediction accuracy of widely used convolutional neural network and natural language processing models with a variety of heterogeneous genomic datasets. Mainly, recurrent neural network models outperform convolutional neural network models in terms of prediction accuracy, overfitting and transferability across the datasets under study. While the perspective of a robust out-of-the-box neural network model is out of reach, we identify certain model characteristics that translate well across datasets and could serve as a baseline model for translational researchers.

59 BASIC BIOLOGICAL SCIENCES↗

Genome-wide functional screens enable the prediction of high activity CRISPR-Cas9 and -Cas12a guides in Yarrowia lipolytica

Abstract Genome-wide functional genetic screens have been successful in discovering genotype-phenotype relationships and in engineering new phenotypes. While broadly applied in mammalian cell lines and in E. coli , use in non-conventional microorganisms has been limited, in part, due to the inability to accurately design high activity CRISPR guides in such species. Here, we develop an experimental-computational approach to sgRNA design that is specific to an organism of choice, in this case the oleaginous yeast Yarrowia lipolytica . A negative selection screen in the absence of non-homologous end-joining, the dominant DNA repair mechanism, was used to generate single guide RNA (sgRNA) activity profiles for both SpCas9 and LbCas12a. This genome-wide data served as input to a deep learning algorithm, DeepGuide, that is able to accurately predict guide activity. DeepGuide uses unsupervised learning to obtain a compressed representation of the genome, followed by supervised learning to map sgRNA sequence, genomic context, and epigenetic features with guide activity. Experimental validation, both genome-wide and with a subset of selected genes, confirms DeepGuide’s ability to accurately predict high activity sgRNAs. DeepGuide provides an organism specific predictor of CRISPR guide activity that with retraining could be applied to other fungal species, prokaryotes, and other non-conventional organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Novosphingobium aromaticivorans LigR coordinates transcription of genes involved in metabolism of multiple types of aromatics

Aromatic compounds are a ubiquitous and diverse family of chemicals with functions as biomolecules, natural products, industrial chemicals, and pollutants. Novosphingobium aromaticivorans DSM 12444 uses multiple inducible pathways to catabolize H-, G-, and S-type aromatics that contain zero, one, or two methoxy groups, respectively. Here, we obtain a systems-level view of the transcriptional control of its aromatic metabolic pathways. Several in vitro analyses found that a N. aromaticivorans homolog of the Sphingobium lignivorans SYK-6 transcription factor LigR bound genomic DNA upstream of genes involved in metabolism of multiple aromatic types. We found that a ΔLigR mutant had growth defects on all three types of aromatics as sole carbon sources. Transcriptomic analysis revealed that LigR was required to increase expression of gene products that function in metabolism of all three aromatic types. We also found that, in media containing both glucose and an aromatic carbon source, the ΔLigR mutant directed intermediates through alternative aromatic metabolic pathways. Protein-DNA binding assays showed that N. aromaticivorans LigR binds immediately upstream of promoters of genes involved in aromatic metabolism. We found that N. aromaticivorans LigR coordinates the expression of enzymes that function in the catabolism of H-, G-, and S-type aromatics, and that there are differences in the role of LigR in N. aromaticivorans and S. lignivorans. A comparative genomic analysis predicted that LigR homologs and the aromatic-metabolizing genes that it directly regulates are often co-localized in the genomes of Sphingomonadales, but often not found in this arrangement in many other known aromatic metabolizing bacteria.

Aromatic Compound Degradation↗

Predicting metabolic modules in incomplete bacterial genomes with MetaPathPredict

The reconstruction of complete microbial metabolic pathways using ‘omics data from environmental samples remains challenging. Computational pipelines for pathway reconstruction that utilize machine learning methods to predict the presence or absence of KEGG modules in incomplete genomes are lacking. Here, we present MetaPathPredict, a software tool that incorporates machine learning models to predict the presence of complete KEGG modules within bacterial genomic datasets. Using gene annotation data and information from the KEGG module database, MetaPathPredict employs deep learning models to predict the presence of KEGG modules in a genome. MetaPathPredict can be used as a command line tool or as a Python module, and both options are designed to be run locally or on a compute cluster. Benchmarks show that MetaPathPredict makes robust predictions of KEGG module presence within highly incomplete genomes.

59 BASIC BIOLOGICAL SCIENCES↗

The state of algal genome quality and diversity

The genomic era of biology has created unprecedented opportunity to study how life works, including understanding evolutionary principles, bioprospecting for novel antibiotics, and genetic manipulation of bioeconomically-relevant species. While genome sequencing and genomic analysis was previously restricted to large-scale projects and consortia, sequencing has become democratized. However, as genomics has become commonplace, the cataloging of sequenced organisms has become challenging, and standardized practices for sequencing, assembly, and annotation have not been adopted. This is equally true in research fields such as algal biology, despite the growing importance of algae in a bio-based economy. Here, in this study, we provide a comprehensive review of the state of eukaryotic algal genomics, explore the quality of algal genome assemblies, and identify the biases and gaps in the current species distribution to inform the development of future genome projects. Overall, we find a trend of declining quality of genomic resources, including a reduction in assembly quality, gene annotation quality, and genome completeness. Potential solutions to improve genome quality include widespread utilization of long read and scaffolding technologies, implementation of standards for assembly quality, evidence-based gene annotation and requisite publication, and support of continued development of gene prediction and genome assessment software.

59 BASIC BIOLOGICAL SCIENCES↗

Integrating chromatin conformation information in a self-supervised learning model improves metagenome binning

Metagenome binning is a key step, downstream of metagenome assembly, to group scaffolds by their genome of origin. Although accurate binning has been achieved on datasets containing multiple samples from the same community, the completeness of binning is often low in datasets with a small number of samples due to a lack of robust species co-abundance information. In this study, we exploited the chromatin conformation information obtained from Hi-C sequencing and developed a new reference-independent algorithm, Metagenome Binning with Abundance and Tetra-nucleotide frequencies—Long Range (metaBAT-LR), to improve the binning completeness of these datasets. This self-supervised algorithm builds a model from a set of high-quality genome bins to predict scaffold pairs that are likely to be derived from the same genome. Then, it applies these predictions to merge incomplete genome bins, as well as recruit unbinned scaffolds. We validated metaBAT-LR’s ability to bin-merge and recruit scaffolds on both synthetic and real-world metagenome datasets of varying complexity. Benchmarking against similar software tools suggests that metaBAT-LR uncovers unique bins that were missed by all other methods.

59 BASIC BIOLOGICAL SCIENCES↗

Application of the metabolic modeling pipeline in KBase to categorize reactions, predict essential genes, and predict pathways in an isolate genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

DOE knowledgebase↗

Application of the Metabolic Modeling Pipeline in KBase to Categorize Reactions, Predict Essential Genes, and Predict Pathways in an Isolate Genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

Allen, Benjamin↗

Genome-wide identification and functional prediction of silicon (Si) transporters in poplar (Populus trichocarpa)

Abstract Silicon (Si) enhances plant tolerance to various biotic and abiotic stressors such as salinity, drought, and heat. In addition, Si can be biomineralized within plants to form organic carbon-containing phytoliths that can have ecosystem-level consequences by contributing to long-term carbon sequestration. Si is taken up and transported in plants via different transporter proteins such as influx transporters (e.g., Lsi1, Lsi6) and efflux transporters (e.g., Lsi2). Additionally, the imported Si can be deposited in plant leaves via silicification process using the Siliplant 1 (e.g., Slp1) protein. Functional homologs of these proteins have been reported in different food crops. Here, we performed a genome-wide analysis to identify different Si transporters and Slp1 homologs in the bioenergy crop poplar ( Populus trichocarpa Torr. and A. Gray ex W. Hook). We identified one channel-type Si influx transporter (PtLsi1; Potri.017G083300), one Si efflux transporter (PtLsi2; Potri.012G144000) and two proteins like Slp1 (PtSlp1a; Potri.004G168600 and PtSlp1b; Potri.009G129900 ) in the P. trichocarpa genome. We found a unique sequence (KPKPPVFKPPPVPI) in PtSlp1a which is repeated six times. Repeated presence of this sequence in PtSlp1a indicates that this protein might be important for silicification processes in P. trichocarpa. The mutation profiles of different Si transporters in a P. trichocarpa genome-wide association study population identified significant and impactful mutations in Potri.004G168600 and Potri.009G129900 . Using a publically accessible database ( http://bar.utoronto.ca/eplant_poplar/ ), digital expression analysis of the putative Si transporters in P. trichocarpa found low to moderate expression in the anticipated tissues, such as roots and leaves. Subcellular localization analysis found that PtLsi1/PtLsi2 are localized in the plasma membrane, whereas PtSlp1a/PtSlp1b are found in the extracellular spaces. Protein–Protein interaction analysis of PtLsi1/PtLsi2 identified Delta-1-pyrroline-5-carboxylate synthase (P5CS) as one of the main interacting partners of PtLsi2, which plays a key role in proline biosynthesis. Proline is a well-known participant in biotic and abiotic stress tolerance in plants. These findings will reinforce future efforts to modify Si accumulation for enhancing plant stress tolerance and carbon sequestration in poplar.

59 BASIC BIOLOGICAL SCIENCES↗

Decoding the chemical language of Suillus fungi: genome mining and untargeted metabolomics uncover terpene chemical diversity

ABSTRACT Ectomycorrhizal fungi establish mutually beneficial relationships with trees, trading nutrients for carbon. Suillus are ectomycorrhizal fungi that are critical to the health of boreal and temperate forest ecosystems. Comparative genomics has identified a high number of non-ribosomal peptide synthetase and terpene biosynthetic gene clusters (BGC) potentially involved in fungal competition and communication. However, the functionality of these BGCs is not known. This study employed co-culture techniques to activate BGC expression and then used metabolomics to investigate the diversity of metabolic products produced by three Suillus species ( Suillus hirtellus EM16, Suillus decipiens EM49, and Suillus cothurnatus VC1858), core members of the pine microbiome. After 28 days of growth on solid media, liquid chromatography–tandem mass spectrometry identified a diverse range of extracellular metabolites (exometabolites) along the interaction zone between Suillus co-cultures. Prenol lipids were among the most abundant chemical classes. Out of the 62 unique terpene BGCs predicted by genome mining, 41 putative prenol lipids (includes 37 putative terpenes) were identified across the three Suillus species using metabolomics. Notably, some terpenes were significantly more abundant in co-culture conditions. For example, we identified a metabolite matching to isomers isopimaric acid, sandaracopimaric acid, and abietic acid, which can be found in pine resin and play important roles in host defense mechanisms and Suillus spore germination. This research highlights the importance of combining genomics and metabolomics to advance our understanding of the chemical diversity underpinning fungal signaling and communication. IMPORTANCE Using a combination of genomics and metabolomics, this study’s findings offer new insights into the chemical diversity of Suillus fungi, which serve a critical role in forest ecosystems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GAL08, an Uncultivated Group of Acidobacteria, Is a Dominant Bacterial Clade in a Neutral Hot Spring

GAL08 are bacteria belonging to an uncultivated phylogenetic cluster within the phylum Acidobacteria . We detected a natural population of the GAL08 clade in sediment from a pH-neutral hot spring located in British Columbia, Canada. To shed light on the abundance and genomic potential of this clade, we collected and analyzed hot spring sediment samples over a temperature range of 24.2–79.8°C. Illumina sequencing of 16S rRNA gene amplicons and qPCR using a primer set developed specifically to detect the GAL08 16S rRNA gene revealed that absolute and relative abundances of GAL08 peaked at 65°C along three temperature gradients. Analysis of sediment collected over multiple years and locations revealed that the GAL08 group was consistently a dominant clade, comprising up to 29.2% of the microbial community based on relative read abundance and up to 4.7 × 10 5 16S rRNA gene copy numbers per gram of sediment based on qPCR. Using a medium quality threshold, 25 single amplified genomes (SAGs) representing these bacteria were generated from samples taken at 65 and 77°C, and seven metagenome-assembled genomes (MAGs) were reconstructed from samples collected at 45–77°C. Based on average nucleotide identity (ANI), these SAGs and MAGs represented three separate species, with an estimated average genome size of 3.17 Mb and GC content of 62.8%. Phylogenetic trees constructed from 16S rRNA gene sequences and a set of 56 concatenated phylogenetic marker genes both placed the three GAL08 bacteria as a distinct subgroup of the phylum Acidobacteria , representing a candidate order ( Ca. Frugalibacteriales) within the class Blastocatellia. Metabolic reconstructions from genome data predicted a heterotrophic metabolism, with potential capability for aerobic respiration, as well as incomplete denitrification and fermentation. In laboratory cultivation efforts, GAL08 counts based on qPCR declined rapidly under atmospheric levels of oxygen but increased slightly at 1% (v/v) O 2 , suggesting a microaerophilic lifestyle.

59 BASIC BIOLOGICAL SCIENCES↗

Tandem repeats in giant archaeal Borg elements undergo rapid evolution and create new intrinsically disordered regions in proteins

Borgs are huge, linear extrachromosomal elements associated with anaerobic methane-oxidizing archaea. Striking features of Borg genomes are pervasive tandem direct repeat (TR) regions. Here, we present six new Borg genomes and investigate the characteristics of TRs in all ten complete Borg genomes. We find that TR regions are rapidly evolving, recently formed, arise independently, and are virtually absent in host Methanoperedens genomes. Flanking partial repeats and A-enriched character constrain the TR formation mechanism. TRs can be in intergenic regions, where they might serve as regulatory RNAs, or in open reading frames (ORFs). TRs in ORFs are under very strong selective pressure, leading to perfect amino acid TRs (aaTRs) that are commonly intrinsically disordered regions. Proteins with aaTRs are often extracellular or membrane proteins, and functionally similar or homologous proteins often have aaTRs composed of the same amino acids. We propose that Borg aaTR-proteins functionally diversify Methanoperedens and all TRs are crucial for specific Borg–host associations and possibly cospeciation.

59 BASIC BIOLOGICAL SCIENCES↗

Using online tools at the Bovine Genome Database to manually annotate genes in the new reference genome

With the availability of a new highly contiguous Bos taurus reference genome assembly (ARS-UCD1.2), it is the opportune time to upgrade the bovine gene set by seeking input from researchers. Furthermore, advances in graphical genome annotation tools now make it possible for researchers to leverage sequence data generated with the latest technologies to collaboratively curate genes. For many years the Bovine Genome Database (BGD) has provided tools such as the APOLLO genome annotation editor to support manual bovine gene curation. The goal of this paper is to explain the reasoning behind the decisions made in the manual gene curation process while providing examples using the existing BGD tools. We will describe the sources of gene annotation evidence provided at the BGD, including RNAseq and Iso-Seq data. We will also explain how to interpret various data visualizations when curating gene models, and will demonstrate the value of manual gene annotation. The process described here can be applied to manual gene curation for other species with similar tools. With a better understanding of manual gene annotation, researchers will be encouraged to edit gene models and contribute to the enhancement of livestock gene sets.

59 BASIC BIOLOGICAL SCIENCES↗

Catabolism of β-5 linked aromatics by Novosphingobium aromaticivorans

ABSTRACT Aromatic compounds are an important source of commodity chemicals traditionally produced from fossil fuels. Aromatics derived from plant lignin can potentially be converted into commodity chemicals through depolymerization followed by microbial funneling of monomers and low molecular weight oligomers. This study investigates the catabolism of the β-5 linked aromatic dimer dehydrodiconiferyl alcohol (DC-A) by the bacterium Novosphingobium aromaticivorans . We used genome-wide screens to identify candidate genes involved in DC-A catabolism. Subsequent in vivo and in vitro analyses of these candidate genes elucidated a catabolic pathway composed of four required gene products and several partially redundant dehydrogenases that convert DC-A to aromatic monomers that can be funneled into the central aromatic metabolic pathway of N. aromaticivorans . Specifically, a newly identified γ-formaldehyde lyase, PcfL, opens the phenylcoumaran ring to form a stilbene and formaldehyde. A lignostilbene dioxygenase, LsdD, then cleaves the stilbene to generate the aromatic monomers vanillin and 5-formylferulate (5-FF). We also showed that the aldehyde dehydrogenase FerD oxidizes 5-FF before it is decarboxylated by LigW, yielding ferulic acid. We found that some enzymes involved in the β-5 catabolism pathway can act on multiple substrates and that some steps in the pathway can be mediated by multiple enzymes, providing new insights into the robust flexibility of aromatic catabolism in N. aromaticivorans . A comparative genomic analysis predicted that the newly discovered β-5 aromatic catabolic pathway is common within the order Sphingomonadales. IMPORTANCE In the transition to a circular bioeconomy, the plant polymer lignin holds promise as a renewable source of industrially important aromatic chemicals. However, since lignin contains aromatic subunits joined by various chemical linkages, producing single chemical products from this polymer can be challenging. One strategy to overcome this challenge is using microbes to funnel a mixture of lignin-derived aromatics into target chemical products. This approach requires strategies to cleave the major inter-unit linkages of lignin to release monomers for funneling into valuable products. In this study, we report newly discovered aspects of a pathway by which the Novosphingobium aromaticivorans DSM12444 catabolizes aromatics joined by the second most common inter-unit linkage in lignin, the β-5 linkage. This work advances our knowledge of aromatic catabolic pathways, laying the groundwork for future metabolic engineering of this and other microbes for optimized conversion of lignin into products.

59 BASIC BIOLOGICAL SCIENCES↗

Apomixis in Farmers’ Fields: Overview, Case Studies from Forage Grasses and Considerations for Future Apomictic Crops

Apomixis occurs naturally in several commercially important species from diverse plant families. While in some of these species apomixis is yet to be exploited in breeding schemes aimed at fixing heterosis, genetic progress and cultivar development, in other species apomixis has been integrated at different stages of breeding. Some of the most relevant examples come from the subfamily Panicoideae, the second largest subfamily of the Poaceae, and are the main focus of this review. The subfamily encompasses many tropical and sub-tropical grasses and grains of worldwide economic importance. Apomictic tropical forages are prime examples of how apomixis can be used and exploited in the development of marketable cultivars, which are essential to the meat and milk production industries globally. The main commercial forages used as grass pastures covering millions of hectares in tropical and sub-tropical regions are polyploids exhibiting gametophytic apomixis that belong to the genus Urochloa spp. (brachiariagrasses) and to the species Megathyrsus maximus (guineagrass). Buffel grass (Cenchrus ciliaris) and Paspalum spp. are other important apomictic forages bred and used in these regions. Breeding involves large germplasm collections from the centers of origin of the species, and for most of them, sexually reproducing diploid plants have been found. Chromosomically duplicated plants that maintain sexual reproduction are used in crosses with apomictic genotypes for the development and selection of cultivars to be marketed or used as progenitors in subsequent breeding cycles. The peculiarities of each genus/species breeding programs, the cultivars obtained from these programs, and the impact of use of marker assisted selection in cultivar development are presented. In addition, the test or implementation of new technologies such as high throughput phenotyping, and the use of machine learning methods for trait prediction and genomic selection are positively impacting the selection and speed of development of new polyploid apomictic cultivars. Furthermore, genetic transformation techniques, including genome editing, provide an additional layer for design of tailor-made, customer-oriented cultivars.

Cenchrus↗

Phylogenetic distribution and experimental characterization of corrinoid production and dependence in soil bacterial isolates

Abstract Soil microbial communities impact carbon sequestration and release, biogeochemical cycling, and agricultural yields. These global effects rely on metabolic interactions that modulate community composition and function. However, the physicochemical and taxonomic complexity of soil and the scarcity of available isolates for phenotypic testing are significant barriers to studying soil microbial interactions. Corrinoids—the vitamin B12 family of cofactors—are critical for microbial metabolism, yet they are synthesized by only a subset of microbiome members. Here, we evaluated corrinoid production and dependence in soil bacteria as a model to investigate the ecological roles of microorganisms involved in metabolic interactions. We isolated and characterized a taxonomically diverse collection of 161 soil bacteria from a single study site. Most corrinoid-dependent bacteria in the collection prefer B12 over other corrinoids, while all tested producers synthesize B12, indicating metabolic compatibility between producers and dependents in the collection. Furthermore, a subset of producers release B12 at levels sufficient to support dependent isolates in laboratory culture at estimated ratios of up to 1000 dependents per producer. Within our isolate collection, we did not find strong phylogenetic patterns in corrinoid production or dependence. Upon investigating trends in the phylogenetic dispersion of corrinoid metabolism categories across sequenced bacteria from various environments, we found that these traits are conserved in 47 out of 85 genera. Together, these phenotypic and genomic results provide evidence for corrinoid-based metabolic interactions among bacteria and provide a framework for the study of nutrient-sharing ecological interactions in microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗

Atlas of Transcription Factor Binding Sites from ENCODE DNase Hypersensitivity Data across 27 Tissue Types

Characterizing the tissue-specific binding sites of transcription factors (TFs) is essential to reconstruct gene regulatory networks and predict functions for non-coding genetic variation. DNase-seq footprinting enables the prediction of genome-wide binding sites for hundreds of TFs simultaneously. Despite the public availability of high-quality DNase-seq data from hundreds of samples, a comprehensive, up-to-date resource for the locations of genomic footprints is lacking. Here, we develop a scalable footprinting workflow using two state-of-the-art algorithms: Wellington and HINT. We apply our workflow to detect footprints in 192 ENCODE DNase-seq experiments and predict the genomic occupancy of 1,515 human TFs in 27 human tissues. We validate that these footprints overlap true-positive TF binding sites from ChIP-seq. We demonstrate that the locations, depth, and tissue specificity of footprints predict effects of genetic variants on gene expression and capture a substantial proportion of genetic risk for complex traits.

59 BASIC BIOLOGICAL SCIENCES↗