Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT

Phage family classification under Caudoviricetes: A review of current tools using the latest ICTV classification framework

Bacteriophages, which are viruses infecting bacteria, are the most ubiquitous and diverse entities in the biosphere. There is accumulating evidence revealing their important roles in shaping the structure of various microbiomes. Thanks to (viral) metagenomic sequencing, a large number of new bacteriophages have been discovered. However, lacking a standard and automatic virus classification pipeline, the taxonomic characterization of new viruses seriously lag behind the sequencing efforts. In particular, according to the latest version of ICTV, several large phage families in the previous classification system are removed. Therefore, a comprehensive review and comparison of taxonomic classification tools under the new standard are needed to establish the state-of-the-art. In this work, we retrained and tested four recently published tools on newly labeled databases. We demonstrated their utilities and tested them on multiple datasets, including the RefSeq, short contigs, simulated metagenomic datasets, and low-similarity datasets. This study provides a comprehensive review of phage family classification in different scenarios and a practical guidance for choosing appropriate taxonomic classification pipelines. To our best knowledge, this is the first review conducted under the new ICTV classification framework. The results show that the new family classification framework overall leads to better conserved groups and thus makes family-level classification more feasible.

Insights into the human gut virome by sampling a population from the Indian subcontinent

Gut virome plays an important role in human physiology but remains poorly understood. This study reports an investigation of the human gut DNA-virome of a previously unexplored ethnic population through metagenomics of faecal samples collected from individuals residing in Northern India. Analysis shows that, similar to the populations investigated earlier, majority of the identified virome belongs to bacteriophages and a smaller fraction (<20 %) consists of viruses that infect animals, archaea, protists, multiple domains or plants. However, crAss-like phages, in this population, are dominated by the genera VI, VII and VIII. Interestingly, it also reveals the presence of a virus family, Sphaerolipoviridae, which has not been detected in the human gut earlier. Viral families, Siphoviridae, Myoviridae, Podoviridae, Microviridae, Herelleviridae and Phycodnaviridae are detected in all of the analysed individuals, which supports the existence of a core virome. Lysogeny-associated genes were found in less than 10 % of the assembled genomes and a negative correlation was observed in the richness of bacterial and free-viral species, suggesting that the dominant lifestyle of gut phage is not lysogenic. This is in contrast to some of the earlier studies. Further, several hundred high-quality viral genomes were recovered. Detailed characterization of these genomes would be useful for understanding the biology of these viruses and their significance in human physiology.

Utilizing Big Data to Identify Tiny Toxic Components: Digitalis

The botanical genus Digitalis is equal parts colorful, toxic, and medicinal, and its bioactive compounds have a long history of therapeutic use. However, with an extremely narrow therapeutic range, even trace amounts of Digitalis can cause adverse effects. Using chemical methods, the United States Food and Drug Administration traced a 1997 case of Digitalis toxicity to a shipment of Plantago (a common ingredient in dietary supplements marketed to improve digestion) contaminated with Digitalis lanata. With increased accessibility to next generation sequencing technology, here we ask whether this case could have been cracked rapidly using shallow genome sequencing strategies (e.g., genome skims). Using a modified implementation of the Site Identification from Short Read Sequences (SISRS) bioinformatics pipeline with whole-genome sequence data, we generated over 2 M genus-level single nucleotide polymorphisms in addition to species-informative single nucleotide polymorphisms. We simulated dietary supplement contamination by spiking low quantities (0-10%) of Digitalis whole-genome sequence data into a background of commonly used ingredients in products marketed for "digestive cleansing" and reliably detected Digitalis at the genus level while also discriminating between Digitalis species. This work serves as a roadmap for the development of novel DNA-based assays to quickly and reliably detect the presence of toxic species such as Digitalis in food products or dietary supplements using genomic methods and highlights the power of harnessing the entire genome to identify botanical species.

Monitoring Microbial Populations and Antibiotic Resistance Gene Enrichment Associated with Arctic Waste Stabilization Ponds

Given that the microbial communities of Arctic waste stabilization ponds (WSPs) are poorly studied to date, our characterization of multiple WSP systems and time points provides important baseline data that will assist with ongoing monitoring of effluent impacts on downstream aquatic ecosystems in the Arctic. This research also identifies indicator amplicon sequence variants (ASVs) of WSPs that will be helpful for future monitoring for WSP effluent attenuation and demonstrates that WSP microbial communities are enriched in antibiotic resistance genes.

Wastewater management in the Canadian Arctic is challenging due to climate extremes, small population sizes, and lack of conventional infrastructure for wastewater treatment. Although many northern communities use waste stabilization ponds (WSPs) as their primary form of wastewater treatment, few studies have explored WSP microbial communities and assessed effluent impacts on receiving waters from a microbiological perspective. Here, we used 16S rRNA gene and metagenome sequencing to characterize WSP and receiving water microbial communities for two time points bracketing the spring WSP thaw in Baker Lake (Nunavut) and compared these results to other Nunavut WSPs in Cambridge Bay and Kugluktuk. Most amplicon sequence variants (ASVs) recovered from these WSP samples belonged to the phylum Proteobacteria, with considerable variation between the three locations and only six ASVs shared among the WSPs at >0.2% relative abundance. Wastewater indicator ASVs for the Baker Lake WSP were identified, and few indicator ASVs were detected in samples originating from other upstream or downstream sites. The metagenomic data revealed a strong enrichment of antibiotic resistance genes for WSP samples relative to downstream and reference samples, especially for genes associated with macrolide resistance. Together, our results provide a baseline characterization for WSP microbial communities, demonstrate how indicator ASVs can be used to monitor attenuation and dilution of effluent microorganisms, and reveal that WSPs can serve as hot spots for antibiotic resistance genes.

IMPORTANCE Given that the microbial communities of Arctic waste stabilization ponds (WSPs) are poorly studied to date, our characterization of multiple WSP systems and time points provides important baseline data that will assist with ongoing monitoring of effluent impacts on downstream aquatic ecosystems in the Arctic. This research also identifies indicator amplicon sequence variants (ASVs) of WSPs that will be helpful for future monitoring for WSP effluent attenuation and demonstrates that WSP microbial communities are enriched in antibiotic resistance genes. Given operational and infrastructure changes anticipated for wastewater treatment systems in the Arctic, baseline data such as these are essential for further development of safe and effective wastewater treatment systems.

Metabolic Contributions of an Alphaproteobacterial Endosymbiont in the Apicomplexan Cardiosporidium cionae

Apicomplexa is a diverse protistan phylum composed almost exclusively of metazoan-infecting parasites, including the causative agents of malaria, cryptosporidiosis, and toxoplasmosis. A single apicomplexan genus, Nephromyces, was described in 2010 as a mutualist partner to its tunicate host. Here we present genomic and transcriptomic data from the parasitic sister species to this mutualist, Cardiosporidium cionae, and its associated bacterial endosymbiont. Cardiosporidium cionae and Nephromyces both infect tunicate hosts, localize to similar organs within these hosts, and maintain bacterial endosymbionts. Though many other protists are known to harbor bacterial endosymbionts, these associations are completely unknown in Apicomplexa outside of the Nephromycidae clade. Our data indicate that a vertically transmitted α-proteobacteria has been retained in each lineage since Nephromyces and Cardiosporidium diverged. This α-proteobacterial endosymbiont has highly reduced metabolic capabilities, but contributes the essential amino acid lysine, and essential cofactor lipoic acid to C. cionae. This partnership likely reduces resource competition with the tunicate host. However, our data indicate that the contribution of the single α-proteobacterial endosymbiont in C. cionae is minimal compared to the three taxa of endosymbionts present in the Nephromyces system, and is a potential explanation for the virulence disparity between these lineages.

Integrated microbiota and metabolite profiles link Crohn's disease to sulfur metabolism

Gut microbial and metabolite alterations have been linked to the pathogenesis of inflammatory bowel diseases. Here we perform a multi-omics microbiome and metabolite analysis of a longitudinal cohort of Crohn's disease patients undergoing autologous hematopoietic stem cell transplantation, and investigational therapy that induces drug free remission in a subset of patients. Via comparison of patients who responded and maintained remission, responded but experienced disease relapse and patients who did not respond to therapy, we identify shared functional signatures that correlate with disease activity despite the variability of gut microbiota profiles at taxonomic level. These signatures reflect the disease state when transferred to gnotobiotic mice. Taken together, the integration of microbiome and metabolite profiles from human cohort and mice improves the predictive modelling of disease outcome, and allows the identification of a network of bacteria-metabolite interactions involving sulfur metabolism as a key mechanism linked to disease activity in Crohn's disease.

Targeted assemblies of cas1 suggest CRISPR-Cas's response to soil warming

There is an increasing interest in the clustered regularly interspaced short palindromic repeats CRISPR-associated protein (CRISPR-Cas) system to reveal potential virus–host dynamics. The universal and most conserved Cas protein, cas1 is an ideal marker to elucidate CRISPR-Cas ecology. We constructed eight Hidden Markov Models (HMMs) and assembled cas1 directly from metagenomes by a targeted-gene assembler, Xander, to improve detection capacity and resolve the diverse CRISPR-Cas systems. The eight HMMs were first validated by recovering all 17 cas1 subtypes from the simulated metagenome generated from 91 prokaryotic genomes across 11 phyla. We challenged the targeted method with 48 metagenomes from a tallgrass prairie in Central Oklahoma recovering 3394 cas1. Among those, 88 were near full length, 5 times more than in de-novo assemblies from the Oklahoma metagenomes. To validate the host assignment by cas1, the targeted-assembled cas1 was mapped to the de-novo assembled contigs. All the phylum assignments of those mapped contigs were assigned independent of CRISPR-Cas genes on the same contigs and consistent with the host taxonomies predicted by the mapped cas1. We then investigated whether 8 years of soil warming altered cas1 prevalence within the communities. A shift in microbial abundances was observed during the year with the biggest temperature differential (mean 4.16 °C above ambient). cas1 prevalence increased and even in the phyla with decreased microbial abundances over the next 3 years, suggesting increasing virus–host interactions in response to soil warming. This targeted method provides an alternative means to effectively mine cas1 from metagenomes and uncover the host communities.

Integrating Computational Methods to Investigate the Macroecology of Microbiomes

Studies in microbiology have long been mostly restricted to small spatial scales. However, recent technological advances, such as new sequencing methodologies, have ushered an era of large-scale sequencing of environmental DNA data from multiple biomes worldwide. These global datasets can now be used to explore long standing questions of microbial ecology. New methodological approaches and concepts are being developed to study such large-scale patterns in microbial communities, resulting in new perspectives that represent a significant advances for both microbiology and macroecology. Here, we identify and review important conceptual, computational, and methodological challenges and opportunities in microbial macroecology. Specifically, we discuss the challenges of handling and analyzing large amounts of microbiome data to understand taxa distribution and co-occurrence patterns. We also discuss approaches for modeling microbial communities based on environmental data, including information on biological interactions to make full use of available Big Data. Finally, we summarize the methods presented in a general approach aimed to aid microbiologists in addressing fundamental questions in microbial macroecology, including classical propositions (such as “everything is everywhere, but the environment selects”) as well as applied ecological problems, such as those posed by human induced global environmental changes.

Isopod holobionts as promising models for lignocellulose degradation

BACKGROUND

Isopods have colonized all environments, partly thanks to their ability to decompose the organic matter. Their enzymatic repertoire, as well as the one of their associated microbiota, has contributed to their colonization success. Together, these holobionts have evolved several interesting life history traits to degrade the plant cell walls, mainly composed of lignocellulose. It has been shown that terrestrial isopods achieve lignocellulose degradation thanks to numerous and diverse CAZymes provided by both the host and its microbiota. Nevertheless, the strategies for lignocellulose degradation seem more diversified in isopods, in particular in aquatic species which are the least studied. Isopods could be an interesting source of valuable enzymes for biotechnological industries of biomass conversion.

RESULTS

To provide new features on the lignocellulose degradation in isopod holobionts, shotgun sequencing of 36 metagenomes of digestive and non-digestive tissues was performed from several populations of four aquatic and terrestrial isopod species. Combined to the 15 metagenomes of an additional species from our previous study, as well as the host transcriptomes, this large dataset allowed us to identify the CAZymes in both the host and the associated microbial communities. Analyses revealed the dominance of Bacteroidetes and Proteobacteria in the five species, covering 36% and 56% of the total bacterial community, respectively. The identification of CAZymes and new enzymatic systems for lignocellulose degradation, such as PULs, cellulosomes and LPMOs, highlights the richness of the strategies used by the isopods and their associated microbiota.

CONCLUSIONS

Altogether, our results show that the isopod holobionts are promising models to study lignocellulose degradation. These models can provide new enzymes and relevant lignocellulose-degrading bacteria strains for the biotechnological industries of biomass conversion.