SeqScreen: accurate and sensitive functional screening of pathogenic sequences via ensemble learning.

The COVID-19 pandemic has emphasized the importance of accurate detection of known and emerging pathogens. However, robust characterization of pathogenic sequences remains an open challenge. To address this need we developed SeqScreen, which accurately characterizes short nucleotide sequences using taxonomic and functional labels and a customized set of curated Functions of Sequences of Concern (FunSoCs) specific to microbial pathogenesis. We show our ensemble machine learning model can label protein-coding sequences with FunSoCs with high recall and precision. SeqScreen is a step towards a novel paradigm of functionally informed synthetic DNA screening and pathogen characterization, available for download at www.gitlab.com/treangenlab/seqscreen .

Seasonal dynamics in taxonomy and function within bacterial and viral metagenomic assemblages recovered from a freshwater agricultural pond.

BACKGROUND: Ponds are important freshwater habitats that support both human and environmental activities. However, relative to their larger counterparts (e.g. rivers, lakes), ponds are understudied, especially with regard to their microbial communities. Our study aimed to fill this knowledge gap by using culture-independent, high-throughput sequencing to assess the dynamics, taxonomy, and functionality of bacterial and viral communities in a freshwater agricultural pond. RESULTS: Water samples (n = 14) were collected from a Mid-Atlantic agricultural pond between June 2017 and May 2018 and filtered sequentially through 1 and 0.2 µm filter membranes. Total DNA was then extracted from each filter, pooled, and subjected to 16S rRNA gene and shotgun sequencing on the Illumina HiSeq 2500 platform. Additionally, on eight occasions water filtrates were processed for viral metagenomes (viromes) using chemical concentration and then shotgun sequenced. A ubiquitous freshwater phylum, Proteobacteria was abundant at all sampling dates throughout the year. However, environmental characteristics appeared to drive the structure of the community. For instance, the abundance of Cyanobacteria (e.g. Nostoc) increased with rising water temperatures, while a storm event appeared to trigger an increase in overall bacterial diversity, as well as the relative abundance of Bacteroidetes. This event was also associated with an increase in the number of antibiotic resistance genes. The viral fractions were dominated by dsDNA of the order Caudovirales, namely Siphoviridae and Myovirdae. CONCLUSIONS: Overall, this study provides one of the largest datasets on pond water microbial ecology to date, revealing seasonal trends in the microbial taxonomic composition and functional potential.

Metagenomic analysis of bacterial and viral assemblages from a freshwater creek and irrigated field reveals temporal and spatial dynamics.

Lotic surface water sites (e.g. creeks) are important resources for localized agricultural irrigation. However, there is concern that microbial contaminants within untreated surface water may be transferred onto irrigated soil and crops. To evaluate this issue, water samples were collected between January 2017 and August 2018 from a freshwater creek used to irrigate kale and radish plants on a small farm in the Mid-Atlantic, United States. In addition, on one sampling date, a field survey was conducted in which additional water (creek source and point-of-use) and soil samples were collected to assess the viral and bacterial communities pre- and post- irrigation. All samples were processed for DNA extracts and shotgun sequenced on the Illumina HiSeq platform. The resulting metagenomic libraries were assembled de novo and taxonomic and functional features were assigned at the contig and peptide level. From these data, we observed that Betaproteobacteria (e.g. Variovorax) dominated the water, both at the source and point-of-use, and Alphaproteobacteria (e.g. Streptomyces) dominated both pre- and post-irrigated soil. Additionally, in the creek source water there were variations in the abundance of the dominant bacterial genera and functional annotations associated with seasonal characteristics (e.g. water temperature). Antibiotic resistance genes and virulence factors were also identified in the creek water and soil, with the majority specific to their respective habitat. Moreover, an analysis of clustered regularly interspaced short palindromic repeat (CRISPR) arrays showed the persistence of certain spacers through time in the creek water, as well as specific interactions between creek bacteriophages and their hosts. Overall, these findings provide a more holistic picture of bacterial and viral composition, dynamics, and interactions within a freshwater creek that can be utilized to further our knowledge on its suitability and safety for irrigation.

Comparative metagenomic analysis of microbial taxonomic and functional variations in untreated surface and reclaimed waters used in irrigation applications.

The use of irrigation water sourced from reclamation facilities and untreated surface water bodies may be a practical solution to attenuate the burden on diminishing groundwater aquifers. However, comprehensive microbial characterizations of these water sources are generally lacking, especially with regard to variations through time and across multiple water types. To address this knowledge gap we used a shotgun metagenomic approach to characterize the taxonomic and functional variations of microbial communities within two agricultural ponds, two freshwater creeks, two brackish rivers, and three water reclamation facilities located in the Mid-Atlantic, United States. Water samples (n = 24) were collected from all sites between October and November 2016, and filtered onto 0.2 µm membrane filters. Filters were then subjected to total DNA extraction and shotgun sequencing on the Illumina HiSeq platform. From these data, we found that Betaproteobacteria dominated the majority of freshwater sites, while Alphaproteobacteria were abundant at times in the brackish waters. One of these brackish sites was also host to a greater abundance of the bacterial genera Gimesia and Microcystis. Furthermore, predicted microbial features (e.g. antibiotic resistance genes (ARGs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) arrays) varied based on specific site and sampling date. ARGs were found across samples, with the diversity and abundance highest in those from a reclamation facility and a wastewater-impacted freshwater creek. Additionally, we identified over 600 CRISPR arrays, containing ∼2600 unique spacers, suggestive of a diverse and often site-specific phage community. Overall, these results provide a better understanding of the complex microbial community in untreated surface and reclaimed waters, while highlighting possible environmental and human health impacts associated with their use in agriculture.

Zero-valent iron sand filtration reduces concentrations of virus-like particles and modifies virome community composition in reclaimed water used for agricultural irrigation.

OBJECTIVE: Zero-valent iron sand filtration can remove multiple contaminants, including some types of pathogenic bacteria, from contaminated water. However, its efficacy at removing complex viral populations, such as those found in reclaimed water used for agricultural irrigation, has not been fully evaluated. Therefore, this study utilized metagenomic sequencing and epifluorescent microscopy to enumerate and characterize viral populations found in reclaimed water and zero-valent iron-sand filtered reclaimed water sampled three times during a larger greenhouse study. RESULTS: Zero-valent iron-sand filtered reclaimed water samples had significantly less virus-like particles than reclaimed water samples at all collection dates, with the reclaimed water averaging between 108 and 109 and the zero-valent iron-sand filtered reclaimed water averaging between 106 and 107 virus-like particles per mL. In addition, for both sample types, viral metagenomes (viromes) were dominated by bacteriophages of the order Caudovirales, largely Siphoviridae, and genes related to DNA metabolism. However, the proportion of sequences homologous to bacteria, as well as the abundance of genes possibly originating from a bacterial host, was higher in the viromes of zero-valent iron-sand filtered reclaimed water samples. Overall, zero-valent iron-sand filtered reclaimed water had a lower total concentration of virus-like particles and a different virome community composition compared to unfiltered reclaimed water.

CRISPR Spacers Indicate Preferential Matching of Specific Virioplankton Genes.

Viral infection exerts selection pressure on marine microbes, as virus-induced cell lysis causes 20 to 50% of cell mortality, resulting in fluxes of biomass into oceanic dissolved organic matter. Archaeal and bacterial populations can defend against viral infection using the clustered regularly interspaced short palindromic repeat (CRISPR)-associated (Cas) system, which relies on specific matching between a spacer sequence and a viral gene. If a CRISPR spacer match to any gene within a viral genome is equally effective in preventing lysis, no viral genes should be preferentially matched by CRISPR spacers. However, if there are differences in effectiveness, certain viral genes may demonstrate a greater frequency of CRISPR spacer matches. Indeed, homology search analyses of bacterioplankton CRISPR spacer sequences against virioplankton sequences revealed preferential matching of replication proteins, nucleic acid binding proteins, and viral structural proteins. Positive selection pressure for effective viral defense is one parsimonious explanation for these observations. CRISPR spacers from virioplankton metagenomes preferentially matched methyltransferase and phage integrase genes within virioplankton sequences. These virioplankton CRISPR spacers may assist infected host cells in defending against competing phage. Analyses also revealed that half of the spacer-matched viral genes were unknown, some genes matched several spacers, and some spacers matched multiple genes, a many-to-many relationship. Thus, CRISPR spacer matching may be an evolutionary algorithm, agnostically identifying those genes under stringent selection pressure for sustaining viral infection and lysis. Investigating this subset of viral genes could reveal those genetic mechanisms essential to virus-host interactions and provide new technologies for optimizing CRISPR defense in beneficial microbes.IMPORTANCE The CRISPR-Cas system is one means by which bacterial and archaeal populations defend against viral infection which causes 20 to 50% of cell mortality in the ocean. We tested the hypothesis that certain viral genes are preferentially targeted for the initial attack of the CRISPR-Cas system on a viral genome. Using CASC, a pipeline for CRISPR spacer discovery, and metagenome data from oceanic microbes and viruses, we found a clear subset of viral genes with high match frequencies to CRISPR spacers. Moreover, we observed a many-to-many relationship of spacers and viral genes. These high-match viral genes were involved in nucleotide metabolism, DNA methylation, and viral structure. It is possible that CRISPR spacer matching is an evolutionary algorithm pointing to those viral genes most important to sustaining infection and lysis. Studying these genes may advance the understanding of virus-host interactions in nature and provide new technologies for leveraging CRISPR-Cas systems in beneficial microbes.

Current progress and future opportunities in applications of bioinformatics for biodefense and pathogen detection: report from the Winter Mid-Atlantic Microbiome Meet-up, College Park, MD, January 10, 2018.

The Mid-Atlantic Microbiome Meet-up (M3) organization brings together academic, government, and industry groups to share ideas and develop best practices for microbiome research. In January of 2018, M3 held its fourth meeting, which focused on recent advances in biodefense, specifically those relating to infectious disease, and the use of metagenomic methods for pathogen detection. Presentations highlighted the utility of next-generation sequencing technologies for identifying and tracking microbial community members across space and time. However, they also stressed the current limitations of genomic approaches for biodefense, including insufficient sensitivity to detect low-abundance pathogens and the inability to quantify viable organisms. Participants discussed ways in which the community can improve software usability and shared new computational tools for metagenomic processing, assembly, annotation, and visualization. Looking to the future, they identified the need for better bioinformatics toolkits for longitudinal analyses, improved sample processing approaches for characterizing viruses and fungi, and more consistent maintenance of database resources. Finally, they addressed the necessity of improving data standards to incentivize data sharing. Here, we summarize the presentations and discussions from the meeting, identifying the areas where microbiome analyses have improved our ability to detect and manage biological threats and infectious disease, as well as gaps of knowledge in the field that require future funding and focus.

RefSeq database growth influences the accuracy of k-mer-based lowest common ancestor species identification.

In order to determine the role of the database in taxonomic sequence classification, we examine the influence of the database over time on k-mer-based lowest common ancestor taxonomic classification. We present three major findings: the number of new species added to the NCBI RefSeq database greatly outpaces the number of new genera; as a result, more reads are classified with newer database versions, but fewer are classified at the species level; and Bayesian-based re-estimation mitigates this effect but struggles with novel genomes. These results suggest a need for new classification approaches specially adapted for large databases.

Agricultural Freshwater Pond Supports Diverse and Dynamic Bacterial and Viral Populations.

Agricultural ponds have a great potential as a means of capture and storage of water for irrigation. However, pond topography (small size, shallow depth) leaves them susceptible to environmental, agricultural, and anthropogenic exposures that may influence microbial dynamics. Therefore, the aim of this project was to characterize the bacterial and viral communities of pond water in the Mid-Atlantic United States with a focus on the late season (October-December), where decreasing temperature and nutrient levels can affect the composition of microbial communities. Ten liters of freshwater from an agricultural pond were sampled monthly, and filtered sequentially through 1 and 0.2 µm filter membranes. Total DNA was then extracted from each filter, and the bacterial communities were characterized using 16S rRNA gene sequencing. The remaining filtrate was chemically concentrated for viruses, DNA-extracted, and shotgun sequenced. Bacterial community profiling showed significant fluctuations over the sampling period, corresponding to changes in the condition of the pond freshwater (e.g., pH, nutrient load). In addition, there were significant differences in the alpha-diversity and core bacterial operational taxonomic units (OTUs) between water fractions filtered through different pore sizes. The viral fraction was dominated by tailed bacteriophage of the order Caudovirales, largely those of the Siphoviridae family. Moreover, while present, genes involved in virulence/antimicrobial resistance were not enriched within the viral fraction during the study period. Instead, the viral functional profile was dominated by phage associated proteins, as well as those related to nucleotide production. Overall, these data suggest that agricultural pond water harbors a diverse core of bacterial and bacteriophage species whose abundance and composition are influenced by environmental variables characteristic of pond topology and the late season.

Family A DNA Polymerase Phylogeny Uncovers Diversity and Replication Gene Organization in the Virioplankton.

Shotgun metagenomics, which allows for broad sampling of viral diversity, has uncovered genes that are widely distributed among virioplankton populations and show linkages to important biological features of unknown viruses. Over 25% of known dsDNA phage carry the DNA polymerase I (polA) gene, making it one of the most widely distributed phage genes. Because of its pivotal role in DNA replication, this enzyme is linked to phage lifecycle characteristics. Previous research has suggested that a single amino acid substitution might be predictive of viral lifestyle. In this study Chesapeake Bay virioplankton were sampled by shotgun metagenomic sequencing (using long and short read technologies). More polA sequences were predicted from this single viral metagenome (virome) than from 86 globally distributed virome libraries (ca. 2,100, and 1,200, respectively). The PolA peptides predicted from the Chesapeake Bay virome clustered with 69% of PolA peptides from global viromes; thus, remarkably the Chesapeake Bay virome captured the majority of known PolA peptide diversity in viruses. This deeply sequenced virome also expanded the diversity of PolA sequences, increasing the number of PolA clusters by 44%. Contigs containing polA sequences were also used to examine relationships between phylogenetic clades of PolA and other genes within unknown viral populations. Phylogenic analysis revealed five distinct groups of phages distinguished by the amino acids at their 762 (Escherichia coli IAI39 numbering) positions and replication genes. DNA polymerase I sequences from Tyr762 and Phe762 groups were most often neighbored by ring-shaped superfamily IV helicases and ribonucleotide reductases (RNRs). The Leu762 groups had non-ring shaped helicases from superfamily II and were further distinguished by an additional helicase gene from superfamily I and the lack of any identifiable RNR genes. Moreover, we found that the inclusion of ribonucleotide reductase associated with PolA helped to further differentiate phage diversity, chiefly within lytic podovirus populations. Altogether, these data show that DNA Polymerase I is a useful marker for observing the diversity and composition of the virioplankton and may be a driving factor in the divergence of phage replication components.

The In-Feed Antibiotic Carbadox Induces Phage Gene Transcription in the Swine Gut Microbiome.

Carbadox is a quinoxaline-di-N-oxide antibiotic fed to over 40% of young pigs in the United States that has been shown to induce phage DNA transduction in vitro; however, the effects of carbadox on swine microbiome functions are poorly understood. We investigated the in vivo longitudinal effects of carbadox on swine gut microbial gene expression (fecal metatranscriptome) and phage population dynamics (fecal dsDNA viromes). Microbial metagenome, transcriptome, and virome sequences were annotated for taxonomic inference and gene function by using FIGfam (isofunctional homolog sequences) and SEED subsystems databases. When the beta diversities of microbial FIGfam annotations were compared, the control and carbadox communities were distinct 2 days after carbadox introduction. This effect was driven by carbadox-associated lower expression of FIGfams (n = 66) related to microbial respiration, carbohydrate utilization, and RNA metabolism (q < 0.1), suggesting bacteriostatic or bactericidal effects within certain populations. Interestingly, carbadox treatment caused greater expression of FIGfams related to all stages of the phage lytic cycle 2 days following the introduction of carbadox (q ≤0.07), suggesting the carbadox-mediated induction of prophages and phage DNA recombination. These effects were diminished by 7 days of continuous carbadox in the feed, suggesting an acute impact. Additionally, the viromes included a few genes that encoded resistance to tetracycline, aminoglycoside, and beta-lactam antibiotics but these did not change in frequency over time or with treatment. The results show decreased bacterial growth and metabolism, prophage induction, and potential transduction of bacterial fitness genes in swine gut bacterial communities as a result of carbadox administration.IMPORTANCE FDA regulations on agricultural antibiotic use have focused on antibiotics that are important for human medicine. Carbadox is an antibiotic not used in humans but frequently used on U.S. pig farms. It is important to study possible side effects of carbadox use because it has been shown to promote bacterial evolution, which could indirectly impact antibiotic resistance in bacteria of clinical importance. Interestingly, the present study shows greater prophage gene expression in feces from carbadox-fed animals than in feces from nonmedicated animals 2 days after the initiation of in-feed carbadox treatment. Importantly, the phage genetic material isolated in this study contained genes that could provide resistance to antibiotics that are important in human medicine, indicating that human-relevant antibiotic resistance genes are mobile between bacteria via phages. This study highlights the collateral effects of antibiotics and demonstrates the need to consider diverse antibiotic effects whenever antibiotics are being used or new regulations are considered.

Novel chaperonins are prevalent in the virioplankton and demonstrate links to viral biology and ecology.

Chaperonins are protein-folding machinery found in all cellular life. Chaperonin genes have been documented within a few viruses, yet, surprisingly, analysis of metagenome sequence data indicated that chaperonin-carrying viruses are common and geographically widespread in marine ecosystems. Also unexpected was the discovery of viral chaperonin sequences related to thermosome proteins of archaea, indicating the presence of virioplankton populations infecting marine archaeal hosts. Virioplankton large subunit chaperonin sequences (GroELs) were divergent from bacterial sequences, indicating that viruses have carried this gene over long evolutionary time. Analysis of viral metagenome contigs indicated that: the order of large and small subunit genes was linked to the phylogeny of GroEL; both lytic and temperate phages may carry group I chaperonin genes; and viruses carrying a GroEL gene likely have large double-stranded DNA (dsDNA) genomes (>70 kb). Given these connections, it is likely that chaperonins are critical to the biology and ecology of virioplankton populations that carry these genes. Moreover, these discoveries raise the intriguing possibility that viral chaperonins may more broadly alter the structure and function of viral and cellular proteins in infected host cells.

Counts and sequences, observations that continue to change our understanding of viruses in nature.

The discovery of abundant viruses in the oceans and on land has ushered in a quarter century of groundbreaking advancements in our understanding of viruses within ecosystems. Two types of observations from environmental samples--direct counts of viral particles and viral metagenomic sequences--have been critical to these discoveries. Accurate direct counts have established ecosystem-scale trends in the impacts of viral infection on microbial host populations and have shown that viral communities within aquatic and soil environments respond to both short term and seasonal environmental change. Direct counts have been critical for estimating viral production rate, a measurement essential to quantifying the implications of viral infection for the biogeochemical cycling of nutrients within ecosystems. While direct counts have defined the magnitude of viral processes; shotgun sequences of environmental viral DNA--virome sequences--have enabled researchers to estimate the diversity and composition of natural viral communities. Virome-enabled studies have found the virioplankton to contain thousands of viral genotypes in communities where the most dominant viral population accounts for a small fraction of total abundance followed by a long tail of diverse populations. Detailed examination of long virome sequences has led to new understanding of genotype-to-phenotype connections within marine viruses and revealed that viruses carry metabolic genes that are important to maintaining cellular energy during viral replication. Increased access to long virome sequences will undoubtedly reveal more genetic secrets of viruses and enable us to build a genomics rulebook for predicting key biological and ecological features of unknown viruses.

Ribonucleotide reductases reveal novel viral diversity and predict biological and ecological features of unknown marine viruses.

Virioplankton play a crucial role in aquatic ecosystems as top-down regulators of bacterial populations and agents of horizontal gene transfer and nutrient cycling. However, the biology and ecology of virioplankton populations in the environment remain poorly understood. Ribonucleotide reductases (RNRs) are ancient enzymes that reduce ribonucleotides to deoxyribonucleotides and thus prime DNA synthesis. Composed of three classes according to O2 reactivity, RNRs can be predictive of the physiological conditions surrounding DNA synthesis. RNRs are universal among cellular life, common within viral genomes and virioplankton shotgun metagenomes (viromes), and estimated to occur within >90% of the dsDNA virioplankton sampled in this study. RNRs occur across diverse viral groups, including all three morphological families of tailed phages, making these genes attractive for studies of viral diversity. Differing patterns in virioplankton diversity were clear from RNRs sampled across a broad oceanic transect. The most abundant RNRs belonged to novel lineages of podoviruses infecting α-proteobacteria, a bacterial class critical to oceanic carbon cycling. RNR class was predictive of phage morphology among cyanophages and RNR distribution frequencies among cyanophages were largely consistent with the predictions of the "kill the winner-cost of resistance" model. RNRs were also identified for the first time to our knowledge within ssDNA viromes. These data indicate that RNR polymorphism provides a means of connecting the biological and ecological features of virioplankton populations.

VIROME: a standard operating procedure for analysis of viral metagenome sequences.

One consistent finding among studies using shotgun metagenomics to analyze whole viral communities is that most viral sequences show no significant homology to known sequences. Thus, bioinformatic analyses based on sequence collections such as GenBank nr, which are largely comprised of sequences from known organisms, tend to ignore a majority of sequences within most shotgun viral metagenome libraries. Here we describe a bioinformatic pipeline, the Viral Informatics Resource for Metagenome Exploration (VIROME), that emphasizes the classification of viral metagenome sequences (predicted open-reading frames) based on homology search results against both known and environmental sequences. Functional and taxonomic information is derived from five annotated sequence databases which are linked to the UniRef 100 database. Environmental classifications are obtained from hits against a custom database, MetaGenomes On-Line, which contains 49 million predicted environmental peptides. Each predicted viral metagenomic ORF run through the VIROME pipeline is placed into one of seven ORF classes, thus, every sequence receives a meaningful annotation. Additionally, the pipeline includes quality control measures to remove contaminating and poor quality sequence and assesses the potential amount of cellular DNA contamination in a viral metagenome library by screening for rRNA genes. Access to the VIROME pipeline and analysis results are provided through a web-application interface that is dynamically linked to a relational back-end database. The VIROME web-application interface is designed to allow users flexibility in retrieving sequences (reads, ORFs, predicted peptides) and search results for focused secondary analyses.