Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequencing data annotation”

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 181 records · Page 10

Galba: genome annotation with miniprot and AUGUSTUS

The Earth Biogenome Project has rapidly increased the number of available eukaryotic genomes, but most released genomes continue to lack annotation of protein-coding genes. In addition, no transcriptome data is available for some genomes. Various gene annotation tools have been developed but each has its limitations. Here, we introduce GALBA, a fully automated pipeline that utilizes miniprot, a rapid protein-to-genome aligner, in combination with AUGUSTUS to predict genes with high accuracy. Accuracy results indicate that GALBA is particularly strong in the annotation of large vertebrate genomes. We also present use cases in insects, vertebrates, and a land plant. GALBA is fully open source and available as a docker image for easy execution with Singularity in high-performance computing environments. Our pipeline addresses the critical need for accurate gene annotation in newly sequenced genomes, and we believe that GALBA will greatly facilitate genome annotation for diverse organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Single-molecule Real-time (SMRT) Isoform Sequencing (Iso-Seq) in Plants: The Status of the Bioinformatics Tools to Unravel the Transcriptome Complexity

Background: The advent of the Single-Molecule Real-time (SMRT) Isoform Sequencing(Iso-Seq) has paved the way to obtain longer full-length transcripts. This method was found tobe much superior in identifying full-length splice variants and other post-transcriptional events ascompared to the Next Generation Sequencing (NGS)-based short read sequencing (RNA-Seq).Several different bioinformatics tools to analyze the Iso-Seq data have been developed and someof them are still being refined to address different aspects of transcriptome complexity. However, acomprehensive summary of the available tools and their utility is still lacking. Objective: Here, we summarized the existing Iso-Seq analysis tools and presented an integratedbioinformatics pipeline for Iso-Seq analysis, which overcomes the limitations of NGS and generateslong contiguous Full-Length Non-Chimeric (FLNC) reads for the analysis of posttranscriptionalevents. Results: In this review, we summarized recent applications of Iso-Seq in plants, which include improvedgenome annotations, identification of novel genes and lncRNAs, identification of fulllengthsplice isoforms, detection of novel Alternative Splicing (AS) and Alternative Polyadenylation(APA) events. In addition, we also discussed the bioinformatics pipeline for comprehensiveIso-Seq data analysis, including how to reduce the error rate in the reads and how to identify andquantify post-transcriptional events. Furthermore, the visualization approach of Iso-Seq was discussedas well. Finally, we discussed methods to combine Iso-Seq data with RNA-Seq for transcriptomequantification. Conclusion: Overall, this review demonstrates that the Iso-Seq is pivotal for analyzing transcriptomecomplexity and this new method offers unprecedented opportunities to comprehensively understandtranscripts diversity.

Biochemistry & Molecular Biology↗

Predicting variable gene content in Escherichia coli using conserved genes

Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%–90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943–0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including “hypothetical proteins” was accurately predicted (F1 = 0.902 [0.898–0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876–0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data.

59 BASIC BIOLOGICAL SCIENCES↗

Show me your secret(ed) weapons: a multifaceted approach reveals a wide arsenal of type III‐secreted effectors in the cucurbit pathogenic bacterium Acidovorax citrulli and novel effectors in the Acidovorax genus

Summary The cucurbit pathogenic bacterium Acidovorax citrulli requires a functional type III secretion system (T3SS) for pathogenicity. In this bacterium, as with Xanthomonas and Ralstonia spp . , an AraC‐type transcriptional regulator, HrpX, regulates expression of genes encoding T3SS components and type III‐secreted effectors (T3Es) . The annotation of a sequenced A. citrulli strain revealed 11 T3E genes. Assuming that this could be an underestimation, we aimed to uncover the T3E arsenal of the A. citrulli model strain, M6. Thorough sequence analysis revealed 51 M6 genes whose products are similar to known T3Es. Furthermore, we combined machine learning and transcriptomics to identify novel T3Es. The machine‐learning approach ranked all A. citrulli M6 genes according to their propensity to encode T3Es. RNA‐Seq revealed differential gene expression between wild‐type M6 and a mutant defective in HrpX: 159 and 28 genes showed significantly reduced and increased expression in the mutant relative to wild‐type M6, respectively. Data combined from these approaches led to the identification of seven novel T3E candidates that were further validated using a T3SS‐dependent translocation assay. These T3E genes encode hypothetical proteins that seem to be restricted to plant pathogenic Acidovorax species. Transient expression in Nicotiana benthamiana revealed that two of these T3Es localize to the cell nucleus and one interacts with the endoplasmic reticulum. This study places A. citrulli among the ‘richest’ bacterial pathogens in terms of T3E cargo. It also revealed novel T3Es that appear to be involved in the pathoadaptive evolution of plant pathogenic Acidovorax species.

Jiménez‐Guerrero, Irene↗

Whole-genome sequencing and single nucleotide polymorphisms in multidrug-resistant clinical isolates of Mycobacterium tuberculosis from the Philippines

Thousands of cases of multidrug-resistant tuberculosis (TB) have been observed in the Philippines, but studies on the Mycobacterium tuberculosis (MTB) genotypes that underlie the observed drug resistance profiles are lacking. This study aimed to analyse the whole genomes of clinical MTB isolates representing various resistance profiles to identify single nucleotide polymorphisms (SNPs) in resistance-associated genes. The genomes of ten MTB isolates cultured from banked sputum sources were sequenced. Bioinformatics analysis consisted of assembly, annotation and SNP identification in genes reported to be associated with resistance to isoniazid (INH), rifampicin (RIF), ethambutol (ETH), streptomycin, pyrazinamide (PZA) and fluoroquinolones (FQs). The draft assemblies covered an average of 97.08% of the expected genome size. Seven of the ten isolates belonged to the Indo-Oceanic lineage/EA12-Manila clade. Two isolates were classified into the Euro-American lineage, whilst the pre-XDR (pre-extensively drug-resistant) isolate was classified under the East Asian/Beijing clade. The SNPs katG Ser315Thr, rpoB Ser450Leu and embB Met306Val were found in INH- (4/7), RIF- (3/6) and ETH-resistant (2/6) isolates, respectively, but not in susceptible isolates. Mutations in the inhA promoter and in the pncA and gyrA genes known to be involved in resistance to INH, PZA and FQs, respectively, were also identified. This study represents the first effort to investigate the whole genomes of Philippine clinical strains of MTB exhibiting various multidrug resistance profiles. Whole-genome data can provide valuable insights to the mechanistic and epidemiological qualities of TB in a high-burden setting such as the Philippines.

60 APPLIED LIFE SCIENCES↗

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↗

CONSTAX2: improved taxonomic classification of environmental DNA markers

Abstract Summary CONSTAX—the CONSensus TAXonomy classifier—was developed for accurate and reproducible taxonomic annotation of fungal rDNA amplicon sequences and is based upon a consensus approach of RDP, SINTAX and UTAX algorithms. CONSTAX2 extends these features to classify prokaryotes as well as eukaryotes and incorporates BLAST-based classifiers to reduce classification errors. Additionally, CONSTAX2 implements a conda-installable command-line tool with improved classification metrics, faster training, multithreading support, capacity to incorporate external taxonomic databases and new isolate matching and high-level taxonomy tools, replete with documentation and example tutorials. Availability and implementation CONSTAX2 is available at https://github.com/liberjul/CONSTAXv2, and is packaged for Linux and MacOS from Bioconda with use under the MIT License. A tutorial and documentation are available at https://constax.readthedocs.io/en/latest/. Data and scripts associated with the manuscript are available at https://github.com/liberjul/CONSTAXv2_ms_code. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

A scaffolded and annotated reference genome of giant kelp (Macrocystis pyrifera)

Abstract Macrocystis pyrifera (giant kelp), is a brown macroalga of great ecological importance as a primary producer and structure-forming foundational species that provides habitat for hundreds of species. It has many commercial uses (e.g. source of alginate, fertilizer, cosmetics, feedstock). One of the limitations to exploiting giant kelp’s economic potential and assisting in giant kelp conservation efforts is a lack of genomic tools like a high quality, contiguous reference genome with accurate gene annotations. Reference genomes attempt to capture the complete genomic sequence of an individual or species, and importantly provide a universal structure for comparison across a multitude of genetic experiments, both within and between species. We assembled the giant kelp genome of a haploid female gametophyte de novo using PacBio reads, then ordered contigs into chromosome level scaffolds using Hi-C. We found the giant kelp genome to be 537 MB, with a total of 35 scaffolds and 188 contigs. The assembly N50 is 13,669,674 with GC content of 50.37%. We assessed the genome completeness using BUSCO, and found giant kelp contained 94% of the BUSCO genes from the stramenopile clade. Annotation of the giant kelp genome revealed 25,919 genes. Additionally, we present genetic variation data based on 48 diploid giant kelp sporophytes from three different Southern California populations that confirms the population structure found in other studies of these populations. This work resulted in a high-quality giant kelp genome that greatly increases the genetic knowledge of this ecologically and economically vital species.

60 APPLIED LIFE SCIENCES↗

Leveraging computational genomics to understand the molecular basis of metal homeostasis

Genome-based data is helping to reveal the diverse strategies plants and algae use to maintain metal homeostasis. In addition to acquisition, distribution and storage of metals, acclimating to feast or famine can involve a wealth of genes that we are just now starting to understand. The fast-paced acquisition of genome-based data, however, is far outpacing our ability to experimentally characterize protein function. Computational genomic approaches are needed to fill the gap between what is known and unknown. To avoid misconstruing bioinformatically derived data, which is the root cause of the inaccurate functional annotations that plague databases, functional inferences from diverse sources and contextualization of that evidence with a robust understanding of protein family evolution is needed. Phylogenomic- and comparative-genomic-based studies can aid in the interpretation of experimental data or provide a spark for the discovery of a new function. These analyses not only lead to novel insight into a target protein's function but can generate thought-provoking insights across protein families.

59 BASIC BIOLOGICAL SCIENCES↗

General Biology 2 Sugar Beet Lab - SU2025 - Final

This module uses the Department of Energy Systems Biology Knowledgebase (KBase) platform to explore topics such as genome assembly, metagenomics, and phylogenomics. Here students will use sequences from DNA that they collected to compare the metagenomes of microbial communities from the rhizosphere of plants grown in fertilized vs. unfertilized soils. Using these data students will evaluate the impact of fertilizer on these communities and how these microbial communities influence soil health and plant growth.

59 BASIC BIOLOGICAL SCIENCES↗

Poplar Interactome: Project Final Report

The feedstock plant Poplar has many advantages over traditional crop plants. Not only Poplar needs low energy input and off season storage as compared to feedstocks such as corn, in the winter season Poplar biomass is stored on the stem/trunk, and Poplar plantations serve as large carbon sink. A key constraint to the expansion of cellulosic bioenergy sources such as in Poplar however, is the negative consequence of converting land use from food crops to energy crops. Therefore in order for Poplar to become a viable energy crop it needs to be grown mostly on marginal land unsuitable agricultural crops. For this we need a better understanding of abiotic stress and adaptation response in poplar. In the process we expected to find new and existing poplar genes and their function that respond to sustain abiotic stress. We carried out an extensive gene expression study on the control untreated and stress (drought, salinity, cold and heat) treated poplar plants. The samples were collected from the stem, leaf and root tissues. The RNA of protein coding genes and regulatory smallRNA genes were sequenced generating more than a billion reads. This is the first such known study in Poplar plants. These were used for quantification and genomic analysis to identify stress responsive genes in poplar. Based on the quantification and genomic analysis, a select set of genes were studied for gene-gene interactions to find their association to stress response. The data was also used to find novel stress responsive genes in poplar that were previously not identified in the Poplar reference genome. The data is made available to the public through the national and international genomic data archives.

59 BASIC BIOLOGICAL SCIENCES↗

Data from: Candidate variants for additive and interactive effects on bioenergy traits in switchgrass (Panicum virgatum L.) identified by genome-wide association analyses

Switchgrass is a promising herbaceous energy crop, but further gains in biomass yield and quality must be achieved to enable a viable bioenergy industry. Developing DNA markers can contribute to such progress, but depiction of genetic bases should be reliable, involving not only simple additive marker effects but also interactions with genetic backgrounds, e.g., ecotypes, or synergies with other markers. We analyzed plant height, carbon content, nitrogen content, and mineral concentration in a diverse panel consisting of 512 genotypes of upland and lowland ecotype. We performed association analyses based on exome capture sequencing and tested 439,170 markers for marginal effects, but also 83,290 markers for marker-by-ecotype interactions and up to 311,445 marker pairs for pairwise interactions. Analyses of pairwise interactions focused on subsets of marker pairs preselected based on marginal marker effects, gene ontology annotation, and pairwise marker associations. Our tests identified 12 significant effects. Homology and gene expression information corroborated seven effects and indicated plausible causal pathways: flowering time and lignin synthesis for plant height; plant growth and senescence for carbon content and mineral concentration. Four pairwise interactions were detected, including three interactions preselected based on pairwise marker correlations. Furthermore, one marker-by-ecotype interaction and one pairwise interaction were confirmed in an independent switchgrass panel. Our analyses identified reliable candidate variants for important bioenergy traits in switchgrass. Moreover, they exemplified the importance of interactive effects for the depiction of genetic bases, and illustrated the usefulness of preselection of marker pairs for identifying pairwise marker interactions in association testing.

09 BIOMASS FUELS↗

Proteome-scale Deployment of Protein Structure Prediction Workflows on the Summit Supercomputer

Deep learning has contributed to major advances in the prediction of protein structure from sequence, a fundamental problem in structural bioinformatics. With predictions now approaching the accuracy of crystallographic experiments, and with accelerators like GPUs and TPUs making inference using large models rapid, genome-level structure prediction becomes an obvious aim. Leadership-class computing resources can be used to perform genome-scale protein structure prediction using state-of-the-art deep learning models, providing a wealth of new data for systems biology applications. Here we describe our efforts to efficiently deploy the AlphaFold v.2 program, for full-proteome structure prediction, at scale on the Oak Ridge Leadership Computing Facility's resources, including the Summit supercomputer. We performed inference to produce the predicted structures for 40,526 protein sequences, corresponding to four prokaryotic proteomes and one plant proteome, using under 4,400 total Summit node hours, equivalent to using the majority of the supercomputer for a little over one hour. We also designed an optimized structure refinement that reduced the time for the relaxation stage of the AlphaFold pipeline by over 10X for longer sequences. We demonstrate the types of analyses that can be performed on proteome-scale collections of sequences, including a search for novel quaternary structures and implications for functional annotation.

Gao, Mu↗

Compilation and utilization of a sorghum transcriptome compendium for gene regulatory network analysis and crop trait engineering

Sorghum bicolor (Sorghum) is a drought and heat tolerant C4 grass crop used to produce grain, forage, biofuels, and other bioproducts. Genetic improvement of sorghum hybrid crops is aided by a large and diverse germplasm, sorghum's diploid inbreeding genetics, and a relatively small genome that has facilitated genomic research. Over the past 20 years, the sorghum research community characterized the cytogenetic and recombinant landscapes of sorghum's 10 chromosomes, sequenced and annotated the sorghum genome, and used that information to identify genes/alleles that modulate flowering time, plant height, seed shattering, and other important traits. More recently, >1000 RNA-seq transcriptome profiles were collected from 15 sorghum genotypes to help understand the genetic basis of variation in growth and development of sorghum stems, tillers, roots, and leaves, and the regulation of biosynthetic pathways that produce epicuticular wax, dhurrin, and RFOs, compounds that contribute to sorghum's resilience. Transcriptome studies were designed to identify differentially expressed genes that are co-expressed during development or in response to a treatment to enable construction of gene regulatory networks. Co-expression and network analysis identified transcription factors and their cognate binding sites in target gene promoters and signaling pathways that modulate gene regulatory networks providing gene editing targets for further trait optimization. RNA-seq data from >20 experiments targeting sorghum organs, tissues, cell types, developmental stages, and responses to environmental conditions (i.e., diel, day-length, shading, water-deficit, temperature) has been compiled in a sorghum transcriptome compendium. The goal of this resource paper is to describe compendium content, accessibility, and a compendium data analysis pipeline and to illustrate the types of information that can be derived from the compendium with a focus on the elucidation of gene regulatory networks useful for guiding the improvement of sorghum traits through gene editing.

RNA-seq↗

Decoding the effects of synonymous variants

Synonymous single nucleotide variants (sSNVs) are common in the human genome but are often overlooked. However, sSNVs can have significant biological impact and may lead to disease. Existing computational methods for evaluating the effect of sSNVs suffer from the lack of gold-standard training/evaluation data and exhibit over-reliance on sequence conservation signals. We developed synVep (synonymous Variant effect predictor), a machine learning-based method that overcomes both of these limitations. Our training data was a combination of variants reported by gnomAD (observed) and those unreported, but possible in the human genome (generated). We used positive-unlabeled learning to purify the generated variant set of any likely unobservable variants. We then trained two sequential extreme gradient boosting models to identify subsets of the remaining variants putatively enriched and depleted in effect. Our method attained 90% precision/recall on a previously unseen set of variants. Furthermore, although synVep does not explicitly use conservation, its scores correlated with evolutionary distances between orthologs in cross-species variation analysis. synVep was also able to differentiate pathogenic vs. benign variants, as well as splice-site disrupting variants (SDV) vs. non-SDVs. Thus, synVep provides an important improvement in annotation of sSNVs, allowing users to focus on variants that most likely harbor effects.

Zishuo Zeng↗

Evolution and comparative genomics of the most common Trichoderma species

Background: The growing importance of the ubiquitous fungal genus Trichoderma (Hypocreales, Ascomycota) requires understanding of its biology and evolution. Many Trichoderma species are used as biofertilizers and biofungicides and T. reesei is the model organism for industrial production of cellulolytic enzymes. In addition, some highly opportunistic species devastate mushroom farms and can become pathogens of humans. A comparative analysis of the first three whole genomes revealed mycoparasitism as the innate feature of Trichoderma. However, the evolution of these traits is not yet understood. Results: We selected 12 most commonly occurring Trichoderma species and studied the evolution of their genome sequences. Trichoderma evolved in the time of the Cretaceous-Palaeogene extinction event 66 (±15) mya, but the formation of extant sections (Longibrachiatum, Trichoderma) or clades (Harzianum/Virens) happened in Oligocene. The evolution of the Harzianum clade and section Trichoderma was accompanied by significant gene gain, but the ancestor of section Longibrachiatum experienced rapid gene loss. The highest number of genes gained encoded ankyrins, HET domain proteins and transcription factors. We also identified the Trichoderma core genome, completely curated its annotation, investigated several gene families in detail and compared the results to those of other fungi. Eighty percent of those genes for which a function could be predicted were also found in other fungi, but only 67% of those without a predictable function. Conclusions: Our study presents a time scaled pattern of genome evolution in 12 Trichoderma species from three phylogenetically distant clades/sections and a comprehensive analysis of their genes. The data offer insights in the evolution of a mycoparasite towards a generalist.

59 BASIC BIOLOGICAL SCIENCES↗

IsoForma: An R Package for Quantifying and Visualizing Positional Isomers in Top-Down LC-MS/MS Data

Proteoforms, the different forms of a protein with sequence variations including post-translational modifications (PTMs), execute vital functions in biological systems such as cell signaling and epigenetic regulation. Precisely defining the stoichiometry of PTMs has been challenging because, in the widely used bottom-up proteomics methods, the detection occurs at the peptide level and thus the link between peptides and their specific modification site is lost, resulting in proteoform ambiguity. Advances in top-down mass spectrometry (MS) technology have permitted the direct characterization of intact proteoforms and their exact number of modification sites, allowing for the relative quantification of positional isomers (PI). Proteins with positional isomers refers to proteoforms with identical total mass and set of modifications but varying PTM site combinations. The relative abundance of PI can be estimated by matching proteoform-specific fragment ions to top-down tandem MS (MS2) data to localize and quantify modifications. However, current approaches heavily rely on manual annotation. Here, we present IsoForma, an open-source R package for relative quantification of PI within a single tool. We benchmarked IsoForma’s performance against two existing workflows and highlight the similarity of the results and improvements in speed. Overall, IsoForma provides a streamlined process, reduces the time of conducting isoform-based analyses, and offers an essential framework for developing customized proteoform analysis workflows. Finally, the software is open source and available at https://github.com/EMSL-Computing/isoforma-lib.

59 BASIC BIOLOGICAL SCIENCES↗

Development of a Knowledgebase (MetRxn) of Metabolites, Reactions and Atom Mappings to Accelerate Discovery and Redesign

With advances in DNA sequencing and genome annotation techniques, the breadth of metabolic knowledge across all kingdoms of life is increasing. The construction of genome-scale models (GSMs) facilitates this distillation of knowledge by systematically accounting for reaction stoichiometry and directionality, gene to protein to reaction relationships, reaction localization among cellular organelles, metabolite transport costs and routes, transcriptional regulation, and biomass composition. Genome-scale reconstructions available now span across all kingdoms of life, from microbes to whole-plant models, and have become indispensable for driving informed metabolic designs and interventions. A key barrier to the pace of this development is our inability to utilize metabolite/reaction information from databases such as BRENDA [1], KEGG [2], MetaCyc [3], etc. due to incompatibilities of representation, duplications, and errors. Duplicate entries constitute a major impediment, where the same metabolite is found with multiple names across databases and models, which significantly slows downs the collating of information from multiple data sources. This can also lead to serious modeling errors such as charge/mass imbalances [4,5] which can thwart model predictive abilities such as identifying synthetic lethal gene pairs and quantifying metabolic flows. Hence, we created the MetRxn database [6] that takes the next step in integrating data from multiple sources and formats to automatically create a standardized knowledgebase. We subsequently deployed this resource to bring about new paradigms in genome-scale metabolic model reconstruction, metabolic flux elucidation through MFA, modeling of microbial communities, and pathway prospecting. This research has enabled the PI’s group to continue building upon research milestones and reach new ones (see list of MetRxn-related publications below).

59 BASIC BIOLOGICAL SCIENCES↗