Standardized and accessible multi-omics bioinformatics workflows through the NMDC EDGE resource
Not Available
59 BASIC BIOLOGICAL SCIENCES↗
Engineering topics
Shakya, Migun
Publications and source records attributed to Shakya, Migun.
Exabiome: Advancing Microbial Science through Exascale Computing
The Exabiome project seeks to improve the understanding of microbiomes through the development of methods for accelerating metagenomic science using exascale computing. This article gives an overview of scientific impact of the three components of the project: metagenome assembly, protein family detection, and comparative analysis of metagenomes. Exabiome developed MetaHipMer, the only metagenome assembler capable of scaling to full exascale systems. MetaHipMer has enabled ground-breaking assemblies on the Frontier supercomputer, with many scientific benefits, such as the discovery of rare species and viral genomes. To investigate protein families, Exabiome developed two exascale tools, PASTIS and HipMCL. Together, these can utilize exascale resources to understand the functional diversity of billions of dark matter proteins and novel protein families. For comparative analysis, Exabiome developed kmerprof, a tool that can be used to compare huge metagenomes for many different scientific purposes, for example, grouping human microbiomes according to body location.
59 BASIC BIOLOGICAL SCIENCES↗
Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows
Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.
60 APPLIED LIFE SCIENCES↗
Migun Shakya, Bio for UC-Davis
Explore the source record for details and available documents.
Biological Science↗
Recombinant And Mix-Infection Finder for SARS-CoV-2 sample
The scientific and public health communities responded to the COVID-19 pandemic with sample acquisition and genome sequencing on a scale that eclipsed all prior sequencing efforts. While this can only be characterized as a resounding success story that has cemented the use of genomics for epidemiological investigations for any future infectious disease outbreak, several retrospective studies are cataloging an array of lessons learned and issues that have yet to be addressed in order to realize the full potential of genomics as a routine biosurveillance tool. We have been both developing methods to accurately assess SARS-CoV-2 genomes from complex samples, and analyzing the large volumes of international data, both at the consensus level and the raw sequencing data. During the course of our investigations and similar to other groups, we have examined COVID-19 samples with signatures from multiple lineages of SARS-CoV-2 and will describe some of our findings during the development of a novel workflow that incorporates detection and reporting of potential co-infection within samples and also highlights any evidence of within-host recombination.
Lo, Chien-Chi↗
Intrahost SARS-CoV-2 k-mer Identification Method (iSKIM) for Rapid Detection of Mutations of Concern Reveals Emergence of Global Mutation Patterns
Despite unprecedented global sequencing and surveillance of SARS-CoV-2, timely identification of the emergence and spread of novel variants of concern (VoCs) remains a challenge. Several million raw genome sequencing runs are now publicly available. We sought to survey these datasets for intrahost variation to study emerging mutations of concern. We developed iSKIM (“intrahost SARS-CoV-2 k-mer identification method”) to relatively quickly and efficiently screen the many SARS-CoV-2 datasets to identify intrahost mutations belonging to lineages of concern. Certain mutations surged in frequency as intrahost minor variants just prior to, or while lineages of concern arose. The Spike N501Y change common to several VoCs was found as a minor variant in 834 samples as early as October 2020. This coincides with the timing of the first detected samples with this mutation in the Alpha/B.1.1.7 and Beta/B.1.351 lineages. Using iSKIM, we also found that Spike L452R was detected as an intrahost minor variant as early as September 2020, prior to the observed rise of the Epsilon/B.1.429/B.1.427 lineages in late 2020. iSKIM rapidly screens for mutations of interest in raw data, prior to genome assembly, and can be used to detect increases in intrahost variants, potentially providing an early indication of novel variant spread.
60 APPLIED LIFE SCIENCES↗
Benchmark datasets for SARS-CoV-2 surveillance bioinformatics
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of coronavirus disease 2019 (COVID-19), has spread globally and is being surveilled with an international genome sequencing effort. Surveillance consists of sample acquisition, library preparation, and whole genome sequencing. This has necessitated a classification scheme detailing Variants of Concern (VOC) and Variants of Interest (VOI), and the rapid expansion of bioinformatics tools for sequence analysis. These bioinformatic tools are means for major actionable results: maintaining quality assurance and checks, defining population structure, performing genomic epidemiology, and inferring lineage to allow reliable and actionable identification and classification. Additionally, the pandemic has required public health laboratories to reach high throughput proficiency in sequencing library preparation and downstream data analysis rapidly. However, both processes can be limited by a lack of a standardized sequence dataset. We identified six SARS-CoV-2 sequence datasets from recent publications, public databases and internal resources. In addition, we created a method to mine public databases to identify representative genomes for these datasets. Using this novel method, we identified several genomes as either VOI/VOC representatives or non-VOI/VOC representatives. To describe each dataset, we utilized a previously published datasets format, which describes accession information and whole dataset information. Additionally, a script from the same publication has been enhanced to download and verify all data from this study.
60 APPLIED LIFE SCIENCES↗
Experimental evidence for the impact of soil viruses on carbon cycling during surface plant litter decomposition
To date, the potential impact of viral communities on biogeochemical cycles in soil has largely been inferred from correlational evidence, such as virus-driven changes in microbial abundances, viral auxiliary metabolic genes, and links with soil physiochemical properties. To more directly test the impact of soil viruses on carbon cycling during plant litter decomposition, we added concentrated viral community suspensions to complex litter decomposer communities in 40-day microcosm experiments. Microbial communities from two New Mexico alpine soils, Pajarito (PJ) and Santa Fe (SF), were inoculated onto grass litter on sand, and three treatments were applied in triplicate to each set of microcosms: addition of buffer (no added virus), live virus (+virus), or killed-virus (+killed-virus) fractions extracted from the same soil. Significant differences in respiration were observed between the +virus and +killed-virus treatments in the PJ, but not the SF microcosms. Bacterial and fungal community composition differed significantly by treatment in both PJ and SF microcosms. Combining data across both soils, viral addition altered links between bacterial and fungal diversity, dissolved organic carbon and total nitrogen. Overall, we demonstrate that increasing viral pressure in complex microbial communities can impact terrestrial biogeochemical cycling but is context-dependent.
54 ENVIRONMENTAL SCIENCES↗
The National Microbiome Data Collaborative Data Portal: an integrated multi-omics microbiome data resource
Explore the source record for details and available documents.
Eloe-Fadrosh, Emiley A↗
NCBI’s Virus Discovery Codeathon: Building “FIVE” —The Federated Index of Viral Experiments API Index
Viruses represent important test cases for data federation due to their genome size and the rapid increase in sequence data in publicly available databases. However, some consequences of previously decentralized (unfederated) data are lack of consensus or comparisons between feature annotations. Unifying or displaying alternative annotations should be a priority both for communities with robust entry representation and for nascent communities with burgeoning data sources. To this end, during this three-day continuation of the Virus Hunting Toolkit codeathon series (VHT-2), a new integrated and federated viral index was elaborated. This Federated Index of Viral Experiments (FIVE) integrates pre-existing and novel functional and taxonomy annotations and virus–host pairings. Variability in the context of viral genomic diversity is often overlooked in virus databases. As a proof-of-concept, FIVE was the first attempt to include viral genome variation for HIV, the most well-studied human pathogen, through viral genome diversity graphs. As per the publication of this manuscript, FIVE is the first implementation of a virus-specific federated index of such scope. FIVE is coded in BigQuery for optimal access of large quantities of data and is publicly accessible. Many projects of database or index federation fail to provide easier alternatives to access or query information. To this end, a Python API query system was developed to enhance the accessibility of FIVE.
59 BASIC BIOLOGICAL SCIENCES↗
A public website for the automated assessment and validation of SARS-CoV-2 diagnostic PCR assays
Abstract Summary Polymerase chain reaction-based assays are the current gold standard for detecting and diagnosing SARS-CoV-2. However, as SARS-CoV-2 mutates, we need to constantly assess whether existing PCR-based assays will continue to detect all known viral strains. To enable the continuous monitoring of SARS-CoV-2 assays, we have developed a web-based assay validation algorithm that checks existing PCR-based assays against the ever-expanding genome databases for SARS-CoV-2 using both thermodynamic and edit-distance metrics. The assay-screening results are displayed as a heatmap, showing the number of mismatches between each detection and each SARS-CoV-2 genome sequence. Using a mismatch threshold to define detection failure, assay performance is summarized with the true-positive rate (recall) to simplify assay comparisons. Availability and implementation The assay evaluation website and supporting software are Open Source and freely available at https://covid19.edgebioinformatics.org/#/assayValidation, https://github.com/jgans/thermonucleotide BLAST and https://github.com/LANL-Bioinformatics/assay_validation. Supplementary information Supplementary data are available at Bioinformatics online.
59 BASIC BIOLOGICAL SCIENCES↗
The National Microbiome Data Collaborative: Enabling Microbiome Science
To harness the potential of microbiome science across the broad range of relevant disciplines, new approaches to data infrastructure and transdisciplinary collaboration are necessary. The National Microbiome Data Collaborative (NMDC) is a new initiative to support microbiome data exploration and discovery through a collaborative integrative data science ecosystem.
Wood-Charlson, Elisha M.↗
A Gene Cluster That Encodes Histone Deacetylase Inhibitors Contributes to Bacterial Persistence and Antibiotic Tolerance in Burkholderia thailandensis
Persister cells are genetically identical variants in a bacterial population that have phenotypically modified their physiology to survive environmental stress. In bacterial pathogens, persisters are able to survive antibiotic treatment and reinfect patients in a frustrating cycle of chronic infection. To better define core persistence mechanisms for therapeutics development, we performed transcriptomics analyses of Burkholderia thailandensis populations enriched for persisters via three methods: flow sorting for low proton motive force, meropenem treatment, and culture aging. Although the three persister-enriched populations generally displayed divergent gene expression profiles that reflect the multimechanistic nature of stress adaptations, there were several common gene pathways activated in two or all three populations. These include polyketide and nonribosomal peptide synthesis, Clp proteases, mobile elements, enzymes involved in lipid metabolism, and ATP binding cassette (ABC) transporter systems. In particular, identification of genes that encode polyketide synthases (PKSs) and fatty acid catabolism factors indicates that generation of secondary metabolites, natural products, and complex lipids could be part of the metabolic program that governs the persistence state. We also found that loss-of-function mutations in the PKS-encoding gene locus BTH_I2366, which plays a role in biosynthesis of histone deacetylase (HDAC) inhibitors, resulted in increased sensitivity to antibiotics targeting DNA replication. Furthermore, treatment of multiple bacterial pathogens with a fatty acid synthesis inhibitor, CP-640186, potentiated the efficacy of meropenem against the persister populations. Altogether, our results suggest that bacterial persisters may exhibit an outwardly dormant physiology but maintain active metabolic processes that are required to maintain persistence.
59 BASIC BIOLOGICAL SCIENCES↗
Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life
There is growing interest in reconstructing phylogenies from the copious amounts of genome sequencing projects that target related viral, bacterial or eukaryotic organisms. To facilitate the construction of standardized and robust phylogenies for disparate types of projects, we have developed a complete bioinformatic workflow, with a web-based component to perform phylogenetic and molecular evolutionary (PhaME) analysis from sequencing reads, draft assemblies or completed genomes of closely related organisms. Furthermore, the ability to incorporate raw data, including some metagenomic samples containing a target organism (e.g. from clinical samples with suspected infectious agents), shows promise for the rapid phylogenetic characterization of organisms within complex samples without the need for prior assembly.
59 BASIC BIOLOGICAL SCIENCES↗
Machine-Learning Classification Suggests That Many Alphaproteobacterial Prophages May Instead Be Gene Transfer Agents
Abstract Many of the sequenced bacterial and archaeal genomes encode regions of viral provenance. Yet, not all of these regions encode bona fide viruses. Gene transfer agents (GTAs) are thought to be former viruses that are now maintained in genomes of some bacteria and archaea and are hypothesized to enable exchange of DNA within bacterial populations. In Alphaproteobacteria, genes homologous to the “head–tail” gene cluster that encodes structural components of the Rhodobacter capsulatus GTA (RcGTA) are found in many taxa, even if they are only distantly related to Rhodobacter capsulatus. Yet, in most genomes available in GenBank RcGTA-like genes have annotations of typical viral proteins, and therefore are not easily distinguished from their viral homologs without additional analyses. Here, we report a “support vector machine” classifier that quickly and accurately distinguishes RcGTA-like genes from their viral homologs by capturing the differences in the amino acid composition of the encoded proteins. Our open-source classifier is implemented in Python and can be used to scan homologs of the RcGTA genes in newly sequenced genomes. The classifier can also be trained to identify other types of GTAs, or even to detect other elements of viral ancestry. Using the classifier trained on a manually curated set of homologous viruses and GTAs, we detected RcGTA-like “head–tail” gene clusters in 57.5% of the 1,423 examined alphaproteobacterial genomes. We also demonstrated that more than half of the in silico prophage predictions are instead likely to be GTAs, suggesting that in many alphaproteobacterial genomes the RcGTA-like elements remain unrecognized.
Evolutionary Biology↗
Novel Insights Into the Spread of Enteric Pathogens Using Genomics
Diarrheal diseases continue to be a leading cause of mortality and morbidity across the globe. The use of next-generation sequencing data has provided major insights into understanding epidemiological spread, antimicrobial resistance, vaccine effects, and basic mechanisms of bacterial and viral evolution. Furthermore, we discuss how genomic data have informed investigations across enteric pathogens, such as Shigella, Vibrio cholerae, Salmonella Typhi, rotavirus, and norovirus.
60 APPLIED LIFE SCIENCES↗
EDGE Bioinformatics
This report gives an overview of Empowering the Development of Genomics Expertise (EDGE) Bioinformatics software.
59 BASIC BIOLOGICAL SCIENCES↗