Caught in the middle with multiple displacement amplification: the myth of pooling for avoiding multiple displacement amplification bias in a metagenome

Overcoming donor variability and risks associated with fecal microbiota transplants through bacteriophage-mediated treatments

BACKGROUND

Fecal microbiota transplantation (FMT) and fecal virome transplantation (FVT, sterile filtrated donor feces) have been effective in treating recurrent Clostridioides difficile infections, possibly through bacteriophage-mediated modulation of the gut microbiome. However, challenges like donor variability, costly screening, coupled with concerns over pathogen transfer (incl. eukaryotic viruses) with FMT or FVT hinder their wider clinical application in treating less acute diseases.

METHODS

To overcome these challenges, we developed methods to broaden FVT's clinical application while maintaining efficacy and increasing safety. Specifically, we employed the following approaches: (1) chemostat-fermentation to reproduce the bacteriophage FVT donor component and remove eukaryotic viruses (FVT-ChP), (2) solvent-detergent treatment to inactivate enveloped viruses (FVT-SDT), and (3) pyronin-Y treatment to inhibit RNA virus replication (FVT-PyT). We assessed the efficacy of these processed FVTs in a C. difficile infection mouse model and compared them with untreated FVT (FVT-UnT), FMT, and saline.

RESULTS

FVT-SDT, FVT-UnT, and FVT-ChP reduced the incidence of mice reaching the humane endpoint (0/8, 2/7, and 3/8, respectively) compared to FMT, FVT-PyT, and saline (5/8, 7/8, and 5/7, respectively) and significantly reduced the load of colonizing C. difficile cells and associated toxin A/B levels. There was a potential elimination of C. difficile colonization, with seven out of eight mice treated with FVT-SDT testing negative with qPCR. In contrast, all other treatments exhibited the continued presence of C. difficile. Moreover, the results were supported by changes in the gut microbiome profiles, cecal cytokine levels, and histopathological findings. Assessment of viral engraftment following FMT/FVT treatment and host-phage correlations analysis suggested that transfer of phages likely were an important contributing factor associated with treatment efficacy.

CONCLUSIONS

This proof-of-concept study shows that specific modifications of FVT hold promise in addressing challenges related to donor variability and infection risks. Two strategies lead to treatments significantly limiting C. difficile colonization in mice, with solvent/detergent treatment and chemostat propagation of donor phages emerging as promising approaches. Video Abstract.

Nanopore and Illumina sequencing reveal different viral populations from human gut samples

The advent of viral metagenomics, or viromics, has improved our knowledge and understanding of global viral diversity. High-throughput sequencing technologies enable explorations of the ecological roles, contributions to host metabolism, and the influence of viruses in various environments, including the human intestinal microbiome. However, bacterial metagenomic studies frequently have the advantage. The adoption of advanced technologies like long-read sequencing has the potential to be transformative in refining viromics and metagenomics. Here, we examined the effectiveness of long-read and hybrid sequencing by comparing Illumina short-read and Oxford Nanopore Technology (ONT) long-read sequencing technologies and different assembly strategies on recovering viral genomes from human faecal samples. Our findings showed that if a single sequencing technology is to be chosen for virome analysis, Illumina is preferable due to its superior ability to recover fully resolved viral genomes and minimise erroneous genomes. While ONT assemblies were effective in recovering viral diversity, the challenges related to input requirements and the necessity for amplification made it less ideal as a standalone solution. However, using a combined, hybrid approach enabled a more authentic representation of viral diversity to be obtained within samples.

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Host DNA depletion methods and genome-centric metagenomics of bovine hindmilk microbiome

Bovine mastitis is a multi-etiological and complex disease, resulting in serious economic consequences for dairy farmers and industry. In recent years, the microbiological evaluation of raw milk has been investigated in-depth using next-generation sequencing approaches such as metataxonomic analysis. Despite this, host DNA is a major concern in the shotgun metagenomic sequencing of microbial communities in milk samples, and it represents a big challenge. In this study, we aimed to evaluate different methods for host DNA depletion and/or microbial DNA enrichment and assess the use of PCR-based whole genome amplification in milk samples with high somatic cell count (SCC) by using short- and long-read sequencing technologies. Our results evidenced that DNA extraction performed differently in terms of host DNA removal, impacting metagenome composition and functional profiles.. Moreover, the ratio of SCC/bacteria ultimately impacts microbial DNA yield, and samples with low SCC (SCC below 100,000 cells/mL) are the most problematic. When milk samples with high SCC (SCC above 200,000 cells/mL) underwent multiple-displacement amplification (MDA), we successfully recovered high-quality metagenome-assembled genomes (MAGs), and long-read sequencing was feasible even for samples with low DNA concentration. By associating MDA and short-read sequencing, we recovered two times more MAGs than in untreated samples, and an ongoing co-infection not reported by traditional methods was detected for mastitis pathogen. Overall, this new approach will improve the detection of mastitis-associated microorganisms and make it possible to examine host-microbiome interactions in bovine mastitis.IMPORTANCENext-generation sequencing technologies have been widely used to gain new insights into the diversity of the microbial community of milk samples and dairy products for different purposes such as microbial safety, profiling of starter cultures, and host-microbiome interactions. Milk is a complex food matrix, and additionally, the presence of host nucleic acid sequences is considered a contaminant in untargeted high-throughput sequencing studies. Therefore, genomic-centric metagenomic studies of milk samples focusing on the health-disease status in dairy cattle are still scarce, which makes it difficult to evaluate the microbial ecophysiology of bovine hindmilk. This study provides an alternative method for genome-centric metagenome studies applied to hindmilk samples with high somatic cell content, which is indispensable to examining host-microbiome interactions in bovine mastitis.

The Diversity of Bacteriophages in the Human Gut

Bacteriophages, commonly referred to as phages, are viruses that infect bacteria and are among the most numerous microorganisms on the planet. They occur throughout nature occupying every habitat where their bacterial hosts can be found. Within these communities, phages are responsible for shaping the bacterial community structure and function through their interactions. Phages shape the community structure and function within the human gut but are also able to influence the human host. As such, there is increased interest in understanding the composition and activity of the gastrointestinal phages, although these studies have been hindered by the difficulties accompanying the study of the human gut. Here, we summarize the methods and findings pertaining to the diversity of the human gastrointestinal phages.

Single-cell transcriptomics and data analyses for prokaryotes-Past, present and future concepts

Transcriptomics, or more specifically mRNA sequencing, is a powerful tool to study gene expression at the single-cell level (scRNA-seq) which enables new insights into a plethora of biological processes. While methods for single-cell RNA-seq in eukaryotes are well established, application to prokaryotes is still challenging. Reasons for that are rigid and diverse cell wall structures hampering lysis, the lack of polyadenylated transcripts impeding mRNA enrichment, and minute amounts of RNA requiring amplification steps before sequencing. Despite those obstacles, several promising scRNA-seq approaches for bacteria have been published recently, albeit difficulties in the experimental workflow and data processing and analysis remain. In particular, bias is often introduced by amplification which makes it difficult to distinguish between technical noise and biological variation. Future optimization of experimental procedures and data analysis algorithms are needed for the improvement of scRNA-seq but also to aid in the emergence of prokaryotic single-cell multi-omics. to help address 21st century challenges in the biotechnology and health sector.

Soil viral diversity, ecology and climate change

Soil viruses are highly abundant and have important roles in the regulation of host dynamics and soil ecology. Climate change is resulting in unprecedented changes to soil ecosystems and the life forms that reside there, including viruses. In this Review, we explore our current understanding of soil viral diversity and ecology, and we discuss how climate change (such as extended and extreme drought events or more flooding and altered precipitation patterns) is influencing soil viruses. Finally, we provide our perspective on future research needs to better understand how climate change will impact soil viral ecology.

Conservation of Genomic Information in Multiple Displacement Amplified Low-Quantity Metagenomic Material from Marine Invertebrates

Marine invertebrate microbiomes have been a rich source of bioactive compounds and interesting genomic features. In cases where the achievable amounts of metagenomic DNA are too low for direct sequencing, multiple displacement amplification (MDA) can be used for whole genome amplification. However, MDA has known limitations which can affect the quality of the resulting genomes and metagenomes. In this study, we evaluated the conservation of biosynthetic gene clusters (BGCs) and enzymes in MDA products from low numbers of prokaryotic cells (estimated 2–850). Marine invertebrate microbiomes collected from Arctic and sub-Arctic areas served as source material. The cells were separated from the host tissue, lysed, and directly subjected to MDA. The MDA products were sequenced by Illumina sequencing. Corresponding numbers of bacteria from a set of three reference strains were treated the same way. The study demonstrated that useful information on taxonomic, BGC, and enzyme diversities was obtainable from such marginal quantities of metagenomic material. Although high levels of assembly fragmentation resulted in most BGCs being incomplete, we conclude that this genome mining approach has the potential to reveal interesting BGCs and genes from hard-to-reach biological sources.

Interpersonal variability of the human gut virome confounds disease signal detection in IBD

Viruses are increasingly recognised as important components of the human microbiome, fulfilling numerous ecological roles including bacterial predation, immune stimulation, genetic diversification, horizontal gene transfer, microbial interactions, and augmentation of metabolic functions. However, our current view of the human gut virome is tainted by previous sequencing requirements that necessitated the amplification of starting nucleic acids. In this study, we performed an original longitudinal analysis of 40 healthy control, 19 Crohn's disease, and 20 ulcerative colitis viromes over three time points without an amplification bias, which revealed and highlighted the interpersonal individuality of the human gut virome. In contrast to a 16 S rRNA gene analysis of matched samples, we show that α- and β-diversity metrics of unamplified viromes are not as efficient at discerning controls from patients with inflammatory bowel disease. Additionally, we explored the intrinsic properties of unamplified gut viromes and show there is considerable interpersonal variability in viral taxa, infrequent longitudinal persistence of intrapersonal viruses, and vast fluctuations in the abundance of temporal viruses. Together, these properties of unamplified faecal viromes confound the ability to discern disease associations but significantly advance toward an unbiased and accurate representation of the human gut virome.

Long-Read Metagenomics and CAZyme Discovery

Microorganisms play a primary role in regulating biogeochemical cycles and are a valuable source of enzymes that have biotechnological applications, such as carbohydrate-active enzymes (CAZymes). However, the inability to culture the majority of microorganisms that exist in natural ecosystems restricts access to potentially novel bacteria and beneficial CAZymes. While commonplace molecular-based culture-independent methods such as metagenomics enable researchers to study microbial communities directly from environmental samples, recent progress in long-read sequencing technologies are advancing the field. We outline key methodological stages that are required as well as describe specific protocols that are currently used for long-read metagenomic projects dedicated to CAZyme discovery.

Strain Tracking of 'Candidatus Liberibacter asiaticus', the Citrus Greening Pathogen, by High-Resolution Microbiome Analysis of Asian Citrus Psyllids

The Asian citrus psyllid, Diaphorina citri, is an invasive insect and a vector of 'Candidatus Liberibacter asiaticus' (CLas), a bacterium whose growth in Citrus species results in huanglongbing (HLB), also known as citrus greening disease. Methods to enrich and sequence CLas from D. citri often rely on biased genome amplification and nevertheless contain significant quantities of host DNA. To overcome these hurdles, we developed a simple pretreatment DNase and filtration (PDF) protocol to remove host DNA and directly sequence CLas and the complete, primarily uncultivable microbiome from D. citri adults. The PDF protocol yielded CLas abundances upward of 60% and facilitated direct measurement of CLas and endosymbiont replication rates in psyllids. The PDF protocol confirmed our lab strains derived from a progenitor Florida CLas strain and accumulated 156 genetic variants, underscoring the utility of this method for bacterial strain tracking. CLas genetic polymorphisms arising in lab-reared psyllid populations included prophage-encoding regions with key functions in CLas pathogenesis, putative antibiotic resistance loci, and a single secreted effector. These variants suggest that laboratory propagation of CLas could result in different phenotypic trajectories among laboratories and could confound CLas physiology or therapeutic design and evaluation if these differences remain undocumented. Finally, we obtained genetic signatures affiliated with Citrus nuclear and organellar genomes, entomopathogenic fungal mitochondria, and commensal bacteria from laboratory-reared and field-collected D. citri adults. Hence, the PDF protocol can directly inform agricultural management strategies related to bacterial strain tracking, insect microbiome surveillance, and antibiotic resistance screening.

Characterization of the human gut virome in metabolic and autoimmune diseases

The intestinal microbiome is dominated by bacteria and plays a pivotal role in the occurrence and development of disease, including several metabolic and autoimmune disorders. While intestinal viral communities, primarily made up of bacteriophages, are also thought to play a role in disease pathogenesis in the gastrointestinal tract, they have received much less attention than intestinal bacteria. Thus, there is limited information about the relationship between bacteriophages and disease. This review explores a potential role for the intestinal viral microbiome in various metabolic and autoimmune diseases.

How clear is our current view on microbial dark matter? (Re-)assessing public MAG & SAG datasets with MDMcleaner

As of today, the majority of environmental microorganisms remain uncultured and is therefore referred to as 'microbial dark matter' (MDM). Hence, genomic insights into these organisms are limited to cultivation-independent approaches such as single-cell- and metagenomics. However, without access to cultured representatives for verifying correct taxon-assignments, MDM genomes may cause potentially misleading conclusions based on misclassified or contaminant contigs, thereby obfuscating our view on the uncultured microbial majority. Moreover, gradual database contaminations by past genome submissions can cause error propagations which affect present as well as future comparative genome analyses. Consequently, strict contamination detection and filtering need to be applied, especially in the case of uncultured MDM genomes. Current genome reporting standards, however, emphasize completeness over purity and the de facto gold standard genome assessment tool, checkM, discriminates against uncultured taxa and fragmented genomes. To tackle these issues, we present a novel contig classification, screening, and filtering workflow and corresponding open-source python implementation called MDMcleaner, which was tested and compared to other tools on mock and real datasets. MDMcleaner revealed substantial contaminations overlooked by current screening approaches and sensitively detects misattributed contigs in both novel genomes and the underlying reference databases, thereby greatly improving our view on 'microbial dark matter'.

Global Composition of the Bacteriophage Community in Honey Bees

The microbial communities in animal digestive systems are critical for host development and health. They stimulate the immune system during development, synthesize important chemical compounds like hormones, aid in digestion, competitively exclude pathogens, etc. Compared to the bacterial and fungal components of the microbiome, we know little about the temporal and spatial dynamics of bacteriophage communities in animal digestive systems. Recently, the bacteriophages of the honey bee gut were characterized in two European bee populations. Most of the bacteriophages described in these two reports were novel, harbored many metabolic genes in their genomes, and had a community structure that suggests coevolution with their bacterial hosts. To describe the conservation of bacteriophages in bees and begin to understand their role in the bee microbiome, we sequenced the virome of Apis mellifera from Austin, TX, and compared bacteriophage compositions among three locations around the world. We found that most bacteriophages from Austin are novel, sharing no sequence similarity with anything in public repositories. However, many bacteriophages are shared among the three bee viromes, indicating specialization of bacteriophages in the bee gut. Our study, along with the two previous bee virome studies, shows that the bee gut bacteriophage community is simple compared to that of many animals, consisting of several hundred types of bacteriophages that primarily infect four of the dominant bacterial phylotypes in the bee gut. IMPORTANCE Viruses that infect bacteria (bacteriophages) are abundant in the microbial communities that live on and in plants and animals. However, our knowledge of the structure, dynamics, and function of these viral communities lags far behind our knowledge of their bacterial hosts. We sequenced the first bacteriophage community of honey bees from the United States and compared the U.S. honey bee bacteriophage community to those of samples from Europe. Our work is an important characterization of an economically critical insect species and shows how bacteriophage communities can contain highly conserved individuals and be highly variable in composition across a wide geographic range.

Record-Breaking Rain Event Altered Estuarine Viral Assemblages

Viruses are the dominant biological entity in the ocean, play a vital role in biogeochemical cycles, and provide their hosts with novel metabolic capabilities through auxiliary metabolic genes (AMGs). Hurricane Harvey was a category 4 hurricane that made landfall on the Texas coast in 2017 and lashed the Houston area with 1.4–1.7 × 1010 m3 of rainfall. In this paper, we aim to characterize how the changes in abiotic conditions brought by Hurricane Harvey altered the viral assemblages of Galveston Bay at the taxonomic level and determine how viral ecosystem functions were altered. Metagenomes of the viruses and their hosts were sequenced from a transect in Galveston Bay over the five weeks following the storm. Our results show that the viral assemblages of Galveston Bay dramatically changed following Hurricane Harvey’s landfall. Of the abiotic parameters measured, salinity had the strongest effect on shaping the viral assemblages. In the five weeks following Hurricane Harvey, there was a steady increase of metabolic genes and putative viral infections. Our study provides the first in-depth look at how marine viral assemblages respond and recover from extreme rainfall events, which models predict will become more frequent and intense with climate change.

RNA Viruses in Aquatic Ecosystems through the Lens of Ecological Genomics and Transcriptomics

Massive amounts of data from nucleic acid sequencing have changed our perspective about diversity and dynamics of marine viral communities. Here, we summarize recent metatranscriptomic and metaviromic studies targeting predominantly RNA viral communities. The analysis of RNA viromes reaffirms the abundance of lytic (+) ssRNA viruses of the order Picornavirales, but also reveals other (+) ssRNA viruses, including RNA bacteriophages, as important constituents of extracellular RNA viral communities. Sequencing of dsRNA suggests unknown diversity of dsRNA viruses. Environmental metatranscriptomes capture the dynamics of ssDNA, dsDNA, ssRNA, and dsRNA viruses simultaneously, unravelling the full complexity of viral dynamics in the marine environment. RNA viruses are prevalent in large size fractions of environmental metatranscriptomes, actively infect marine unicellular eukaryotes larger than 3 µm, and can outnumber bacteriophages during phytoplankton blooms. DNA and RNA viruses change abundance on hourly timescales, implying viral control on a daily temporal basis. Metatranscriptomes of cultured protists host a diverse community of ssRNA and dsRNA viruses, often with multipartite genomes and possibly persistent intracellular lifestyles. We posit that RNA viral communities might be more diverse and complex than formerly anticipated and that the influence they exert on community composition and global carbon flows in aquatic ecosystems may be underestimated.

Metagenomics Approaches for Improving Food Safety: A Review

ABSTRACT

Advancements in next-generation sequencing technology have dramatically reduced the cost and increased the ease of microbial whole genome sequencing. This approach is revolutionizing the identification and analysis of foodborne microbial pathogens, facilitating expedited detection and mitigation of foodborne outbreaks, improving public health outcomes, and limiting costly recalls. However, next-generation sequencing is still anchored in the traditional laboratory practice of the selection and culture of a single isolate. Metagenomic-based approaches, including metabarcoding and shotgun and long-read metagenomics, are part of the next disruptive revolution in food safety diagnostics and offer the potential to directly identify entire microbial communities in a single food, ingredient, or environmental sample. In this review, metagenomic-based approaches are introduced and placed within the context of conventional detection and diagnostic techniques, and essential considerations for undertaking metagenomic assays and data analysis are described. Recent applications of the use of metagenomics for food safety are discussed alongside current limitations and knowledge gaps and new opportunities arising from the use of this technology.

HIGHLIGHTS

Alterations in the Gut Microbiota and Their Metabolites in Colorectal Cancer: Recent Progress and Future Prospects

Colorectal cancer (CRC) is a leading cause of cancer morbidity and mortality worldwide. The etiology and pathogenesis of CRC remain unclear. A growing body of evidence suggests dysbiosis of gut bacteria can contribute to the occurrence and development of CRC by generating harmful metabolites and changing host physiological processes. Metabolomics, a systems biology method, will systematically study the changes in metabolites in the physiological processes of the body, eventually playing a significant role in the detection of metabolic biomarkers and improving disease diagnosis and treatment. Metabolomics, in particular, has been highly beneficial in tracking microbially derived metabolites, which has substantially advanced our comprehension of host-microbiota metabolic interactions in CRC. This paper has briefly compiled recent research progress of the alterations of intestinal flora and its metabolites associated with CRC and the application of association analysis of metabolomics and gut microbiome in the diagnosis, prevention, and treatment of CRC; furthermore, we discuss the prospects for the problems and development direction of this association analysis in the study of CRC. Gut microbiota and their metabolites influence the progression and causation of CRC, and the association analysis of metabolomics and gut microbiome will provide novel strategies for the prevention, diagnosis, and therapy of CRC.

Viral Characteristics of the Warm Atlantic and Cold Arctic Water Masses in the Nordic Seas

Nordic Seas are the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean with complex water masses, experiencing an abrupt climate change. Though knowledge of the marine virosphere has expanded rapidly, the diversity of viruses and their relationships with host cells and water masses in the Nordic Seas remain to be fully revealed. Here, we establish the Nordic Sea DNA virome (NSV) data set of 55,315 viral contigs including 1,478 unique viral populations from seven stations influenced by both the warm Atlantic and cold Arctic water masses. Caudovirales dominated in the seven NSVs, especially in the warm Atlantic waters. The major giant nucleocytoplasmic large DNA viruses (NCLDVs) contributed a significant proportion of the classified viral contigs in the NSVs (32.2%), especially in the cold Arctic waters (44.9%). The distribution patterns of Caudovirales and NCLDVs were a reflection of the community structure of their hosts in the corresponding water masses and currents. Latitude, pH, and flow speed were found to be key factors influencing the microbial communities and coinfluencing the variation of viral communities. Network analysis illustrated the tight coupling between the variation of viral communities and microbial communities in the Nordic Seas. This study suggests a probable linkage between viromes, host cells, and surface water masses from both the cool Arctic and warm Atlantic Oceans. IMPORTANCE This is a systematic study of Nordic Sea viromes using metagenomic analysis. The viral diversity, community structure, and their relationships with host cells and the complex water masses from both the cool Arctic and the warm Atlantic oceans were illustrated. The NCLDVs and Caudovirales are proposed as the viral characteristics of the cold Arctic and warm Atlantic waters, respectively. This study provides an important background for the viromes in the subarctic seas connecting the Arctic Ocean and North Atlantic Ocean and sheds light on their responses to abrupt climate change in the future.

The intestinal virome: lessons from animal models

Mucosal surfaces in contact with the environment host specific microbiota. The intestinal tract harbours the most abundant and diverse bacterial and viral populations interacting with each other as well as with the host. Viruses of the microbiota are important components of this ecosystem, as shown by viral alterations associated with various pathologies. However, practical and ethical constraints limit functional studies of the virome in humans, making animal models invaluable experimental tools to understand its impact on intestinal physiology. In this review, we present the recent advances in the study of virome in animal models. We focus on the strategies used to characterise viral changes in disease models and approaches to modulate the microbiota using viruses. In reviewing the interplay between viruses, bacteria, and the animal host, we highlight the potential and limitations of these models in elucidating the role of the virome in determining human health and disease.

Next-Generation Food Research: Use of Meta-Omic Approaches for Characterizing Microbial Communities Along the Food Chain

Microorganisms exist along the food chain and impact the quality and safety of foods in both positive and negative ways. Identifying and understanding the behavior of these microbial communities enable the implementation of preventative or corrective measures in public health and food industry settings. Current culture-dependent microbial analyses are time-consuming and target only specific subsets of microbes. However, the greater use of culture-independent meta-omic approaches has the potential to facilitate a thorough characterization of the microbial communities along the food chain. Indeed, these methods have shown potential in contributing to outbreak investigation, ensuring food authenticity, assessing the spread of antimicrobial resistance, tracking microbial dynamics during fermentation and processing, and uncovering the factors along the food chain that impact food quality and safety. This review examines the community-based approaches, and particularly the application of sequencing-based meta-omics strategies, for characterizing microbial communities along the food chain. Expected final online publication date for the Annual Review of Food Science and Technology, Volume 13 is March 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Comparison of PCR versus PCR-Free DNA Library Preparation for Characterising the Human Faecal Virome

The human intestinal microbiota is abundant in viruses, comprising mainly bacteriophages, occasionally outnumbering bacteria 10:1 and is termed the virome. Due to their high genetic diversity and the lack of suitable tools and reference databases, the virome remains poorly characterised and is often referred to as "viral dark matter". However, the choice of sequencing platforms, read lengths and library preparation make study design challenging with respect to the virome. Here we have compared the use of PCR and PCR-free methods for sequence-library construction on the Illumina sequencing platform for characterising the human faecal virome. Viral DNA was extracted from faecal samples of three healthy donors and sequenced. Our analysis shows that most variation was reflecting the individually specific faecal virome. However, we observed differences between PCR and PCR-free library preparation that affected the recovery of low-abundance viral genomes. Using three faecal samples in this study, the PCR library preparation samples led to a loss of lower-abundance vOTUs evident in their PCR-free pairs (vOTUs 128, 6202 and 8364) and decreased the alpha-diversity indices (Chao1 p-value = 0.045 and Simpson p-value = 0.044). Thus, differences between PCR and PCR-free methods are important to consider when investigating "rare" members of the gut virome, with these biases likely negligible when investigating moderately and highly abundant viruses.

Viromes outperform total metagenomes in revealing the spatiotemporal patterns of agricultural soil viral communities

Viruses are abundant yet understudied members of soil environments that influence terrestrial biogeochemical cycles. Here, we characterized the dsDNA viral diversity in biochar-amended agricultural soils at the preplanting and harvesting stages of a tomato growing season via paired total metagenomes and viral size fraction metagenomes (viromes). Size fractionation prior to DNA extraction reduced sources of nonviral DNA in viromes, enabling the recovery of a vaster richness of viral populations (vOTUs), greater viral taxonomic diversity, broader range of predicted hosts, and better access to the rare virosphere, relative to total metagenomes, which tended to recover only the most persistent and abundant vOTUs. Of 2961 detected vOTUs, 2684 were recovered exclusively from viromes, while only three were recovered from total metagenomes alone. Both viral and microbial communities differed significantly over time, suggesting a coupled response to rhizosphere recruitment processes and/or nitrogen amendments. Viral communities alone were also structured along an 18 m spatial gradient. Overall, our results highlight the utility of soil viromics and reveal similarities between viral and microbial community dynamics throughout the tomato growing season yet suggest a partial decoupling of the processes driving their spatial distributions, potentially due to differences in dispersal, decay rates, and/or sensitivities to soil heterogeneity.

Quantifying genome DNA during whole-genome amplification via quantitative real-time multiple displacement amplification

DNA quantification is important in the research of life sciences. In an independent quantification process, the extracted part of a DNA sample is usually difficult to be recycled for further use while the widely used real-time PCR is used to count the copies with certain sequences. Based on the popular multiple displacement amplification (MDA), we proposed and performed quantitative real-time MDA to obtain the information of template amount based on fluorescence signals while amplifying whole-genome DNA. The detection limit of real-time MDA was as low as 0.5 pg μl-1 (5 pg DNA input), offering the whole-genome research a promising tool to quantify the entire DNA during amplification without sacrificing sample completeness or introducing redundant steps.

Faecal virome transplantation decreases symptoms of type 2 diabetes and obesity in a murine model

OBJECTIVE

Development of obesity and type 2 diabetes (T2D) are associated with gut microbiota (GM) changes. The gut viral community is predominated by bacteriophages (phages), which are viruses that attack bacteria in a host-specific manner. The antagonistic behaviour of phages has the potential to alter the GM. As a proof-of-concept, we demonstrate the efficacy of faecal virome transplantation (FVT) from lean donors for shifting the phenotype of obese mice into closer resemblance of lean mice.

DESIGN

The FVT consisted of viromes with distinct profiles extracted from the caecal content of mice from different vendors that were fed a low-fat (LF) diet for 14 weeks. Male C57BL/6NTac mice were divided into five groups: LF (as diet control), high-fat (HF) diet, HF+ampicillin (Amp), HF+Amp+FVT and HF+FVT. At weeks 6 and 7 of the study, the HF+FVT and HF+Amp+FVT mice were treated with FVT by oral gavage. The Amp groups were treated with Amp 24 hours prior to first FVT treatment.

RESULTS

Six weeks after first FVT, the HF+FVT mice showed a significant decrease in weight gain compared with the HF group. Further, glucose tolerance was comparable between the LF and HF+FVT mice, while the other HF groups all had impaired glucose tolerance. These observations were supported by significant shifts in GM composition, blood plasma metabolome and expression levels of genes associated with obesity and T2D development.

CONCLUSIONS

Transfer of caecal viral communities from mice with a lean phenotype into mice with an obese phenotype led to reduced weight gain and normalised blood glucose parameters relative to lean mice. We hypothesise that this effect is mediated via FVT-induced GM changes.

Characterization of the Vaginal DNA Virome in Health and Dysbiosis

Bacterial vaginosis (BV) is characterized by a reduction in Lactobacillus (L.) spp. abundance and increased abundance of facultative anaerobes, such as Gardnerella spp. BV aetiology is not fully understood; however, bacteriophages could play a pivotal role in the perturbation of the vaginal bacterial community. We investigated the vaginal viral community, including bacteriophages and the association to the bacterial community and BV-status. Vaginal samples from 48 patients undergoing IVF treatment for non-female factor infertility were subjected to metagenomic sequencing of purified virus-like particles. The vaginal viral community was characterized and correlated with the BV-status by Nugent score, bacterial community, structure, and the presence of key vaginal bacterial species. The majority of identified vaginal viruses belonged to the class of double-stranded DNA bacteriophages, with eukaryotic viruses constituting 4% of the total reads. Clear links between the viral community composition and BV (q = 0.006, R = 0.26) as well as the presence of L. crispatus (q = 0.001, R = 0.43), L. iners, Gardnerella spp., and Atopobium vaginae were found (q < 0.002, R > 0.15). The eukaryotic viral community also correlated with BV-status (q = 0.018, R = 0.20). In conclusion, the vaginal virome was clearly linked with bacterial community structure and BV-status.

Shining Light on Human Gut Bacteriophages

The human gut is a complex environment that contains a multitude of microorganisms that are collectively termed the microbiome. Multiple factors have a role to play in driving the composition of human gut bacterial communities either toward homeostasis or the instability that is associated with many disease states. One of the most important forces are likely to be bacteriophages, bacteria-infecting viruses that constitute by far the largest portion of the human gut virome. Despite this, bacteriophages (phages) are the one of the least studied residents of the gut. This is largely due to the challenges associated with studying these difficult to culture entities. Modern high throughput sequencing technologies have played an important role in improving our understanding of the human gut phageome but much of the generated sequencing data remains uncharacterised. Overcoming this requires database-independent bioinformatic pipelines and even those phages that are successfully characterized only provide limited insight into their associated biological properties, and thus most viral sequences have been characterized as “viral dark matter.” Fundamental to understanding the role of phages in shaping the human gut microbiome, and in turn perhaps influencing human health, is how they interact with their bacterial hosts. An essential aspect is the isolation of novel phage-bacteria host pairs by direct isolation through various screening methods, which can transform in silico phages into a biological reality. However, this is also beset with multiple challenges including culturing difficulties and the use of traditional methods, such as plaquing, which may bias which phage-host pairs that can be successfully isolated. Phage-bacteria interactions may be influenced by many aspects of complex human gut biology which can be difficult to reproduce under laboratory conditions. Here we discuss some of the main findings associated with the human gut phageome to date including composition, our understanding of phage-host interactions, particularly the observed persistence of virulent phages and their hosts, as well as factors that may influence these highly intricate relationships. We also discuss current methodologies and bottlenecks hindering progression in this field and identify potential steps that may be useful in overcoming these hurdles.

Targeted Cell Sorting Combined With Single Cell Genomics Captures Low Abundant Microbial Dark Matter With Higher Sensitivity Than Metagenomics

Rare members of environmental microbial communities are often overlooked and unexplored, primarily due to the lack of techniques capable of acquiring their genomes. Chloroflexi belong to one of the most understudied phyla, even though many of its members are ubiquitous in the environment and some play important roles in biochemical cycles or biotechnological applications. We here used a targeted cell-sorting approach, which enables the selection of specific taxa by fluorescent labeling and is compatible with subsequent single-cell genomics, to enrich for rare Chloroflexi species from a wastewater-treatment plant and obtain their genomes. The combined workflow was able to retrieve a substantially higher number of novel Chloroflexi draft genomes with much greater phylogenetical diversity when compared to a metagenomics approach from the same sample. The method offers an opportunity to access genetic information from rare biosphere members which would have otherwise stayed hidden as microbial dark matter and can therefore serve as an essential complement to cultivation-based, metagenomics, and microbial community-focused research approaches.

Comparison of Illumina MiSeq and the Ion Torrent PGM and S5 platforms for whole-genome sequencing of picornaviruses and caliciviruses

Next-generation sequencing is a powerful tool for virological surveillance. While Illumina® and Ion Torrent® sequencing platforms are used extensively for generating viral RNA genome sequences, there is limited data comparing different platforms. The Illumina MiSeq, Ion Torrent PGM and Ion Torrent S5 platforms were evaluated using a panel of sixteen specimens containing picornaviruses and human caliciviruses (noroviruses and sapoviruses). The specimens were processed, using combinations of three library preparation and five sequencing kits, to assess the quality and completeness of assembled viral genomes, and an estimation of cost per sample to generate the data was calculated. The choice of library preparation kit and sequencing platform was found to impact the breadth of genome coverage and accuracy of consensus viral genomes. The Ion Torrent S5 510 chip runs produced more reads at a lower cost per sample than the highest output Ion Torrent PGM 318 chip run, and generated the highest proportion of reads for enterovirus D68 samples. However, indels at homopolymer regions impacted the accuracy of consensus genome sequences. For lower throughput sequencing runs (i.e., Ion Torrent 510 and Illumina MiSeq Nano V2), the cost per sample was lower on the MiSeq platform, whereas with higher throughput runs (Ion Torrent 530 and Illumina MiSeq V2) there is less of a difference in the cost per sample between the two sequencing platforms ($5.47-$10.25 more per sample for an Ion Torrent 530 chip run when multiplexing 24 samples). These findings suggest that the Ion Torrent S5 and Illumina MiSeq platforms are both viable options for genomic sequencing of RNA viruses, each with specific advantages and tradeoffs.

Diversity and potential biogeochemical impacts of viruses in bulk and rhizosphere soils

Viruses can affect microbial dynamics, metabolism and biogeochemical cycles in aquatic ecosystems. However, viral diversity and functions in agricultural soils are poorly known, especially in the rhizosphere. We used virome analysis of eight rhizosphere and bulk soils to study viral diversity and potential biogeochemical impacts in an agro-ecosystem. The order Caudovirales was the predominant viral type in agricultural soils, with Siphoviridae being the most abundant family. Phylogenetic analysis of the terminase large subunit of Caudovirales identified high viral diversity and three novel groups. Viral community composition differed significantly between bulk and rhizosphere soils. Soil pH was the main environmental driver of the viral community structure. Remarkably, abundant auxiliary carbohydrate-active enzyme (CAZyme) genes were detected in viromes, including glycoside hydrolases, carbohydrate esterases and carbohydrate-binding modules. These results demonstrate that virus-encoded putative auxiliary metabolic genes or metabolic genes that may change bacterial metabolism and indirectly contribute to biogeochemical cycling, especially carbon cycling, in agricultural soil.

Metagenome Mining Reveals Hidden Genomic Diversity of Pelagimyophages in Aquatic Environments

SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 μm), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear.

IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

Metagenomic Analysis of the Diversity of DNA Viruses in the Surface and Deep Sea of the South China Sea

A metagenomic analysis of the viral community from five surface and five deep sea water (>2000 m below the surface, mbs) samples collected from the central basin of the South China Sea and adjacent Northwest Pacific Ocean during July-August 2017 was conducted herein. We builded up a South China Sea DNA virome (SCSV) dataset of 29,967 viral Operational Taxonomic Units (vOTUs), which is comparable to the viral populations from the original Tara Ocean and Malaspina expeditions. The most abundant and widespread viral populations were from the uncultivated viruses annotated from the viral metagenomics. Only 74 and 37 vOTUs have similarity with the reported genomes from the cultivated viruses and the single-virus genomics, respectively. The community structures of deep sea viromes in the SCSV were generally different from the surface viromes. The carbon flux and nutrients (PO₄ and NOx) were related to the surface and deep sea viromes in the SCSV, respectively. In the SCSV, the annotated vOTUs could be affiliated to the cultivated viruses mainly including Pelagibacter (SAR11) phage HTVC010P, Prochlorococcus phages (P-GSP1, P-SSM4, and P-TIM68), Cyanophages (MED4-184 and MED4-117) and Mycobacterium phages (Sparky and Squirty). It indicated that phage infection to the SAR11 cluster may occur ubiquitously and has significant impacts on bathypelagic SAR11 communities in the deep sea. Meanwhile, as Prochlorococcus is prominently distributed in the euphotic ocean, the existence of their potential phages in the deep sea suggested the sedimentation mechanism might contribute to the formation of the deep sea viromes. Intriguingly, the presence of Mycobacterium phages only in the deep sea viromes, suggests inhabitance of endemic viral populations in the deep sea viromes in the SCSV. This study provided an insight of the viral community in the South China Sea and for the first time uncovered the deep sea viral diversity in the central basin of the South China Sea.

A simple, reproducible and cost-effective procedure to analyse gut phageome: from phage isolation to bioinformatic approach

The microbiota of the human gut is a complex and rich community where bacteria and their viruses, the bacteriophages, are dominant. There are few studies on the phage community and no clear standard for isolating them, sequencing and analysing their genomes. Since this makes comparisons between studies difficult, we aimed at defining an easy, low-cost, and reproducible methodology. We analysed five different techniques to isolate phages from human adult faeces and developed an approach to analyse their genomes in order to quantify contamination and classify phage contigs in terms of taxonomy and lifestyle. We chose the polyethylene glycol concentration method to isolate phages because of its simplicity, low cost, reproducibility, and of the high number and diversity of phage sequences that we obtained. We also tested the reproducibility of this method with multiple displacement amplification (MDA) and showed that MDA severely decreases the phage genetic diversity of the samples and the reproducibility of the method. Lastly, we studied the influence of sequencing depth on the analysis of phage diversity and observed the beginning of a plateau for phage contigs at 20,000,000 reads. This work contributes to the development of methods for the isolation of phages in faeces and for their comparative analysis.

Towards optimized viral metagenomes for double-stranded and single-stranded DNA viruses from challenging soils

Soils impact global carbon cycling and their resident microbes are critical to their biogeochemical processing and ecosystem outputs. Based on studies in marine systems, viruses infecting soil microbes likely modulate host activities via mortality, horizontal gene transfer, and metabolic control. However, their roles remain largely unexplored due to technical challenges with separating, isolating, and extracting DNA from viruses in soils. Some of these challenges have been overcome by using whole genome amplification methods and while these have allowed insights into the identities of soil viruses and their genomes, their inherit biases have prevented meaningful ecological interpretations. Here we experimentally optimized steps for generating quantitatively-amplified viral metagenomes to better capture both ssDNA and dsDNA viruses across three distinct soil habitats along a permafrost thaw gradient. First, we assessed differing DNA extraction methods (PowerSoil, Wizard mini columns, and cetyl trimethylammonium bromide) for quantity and quality of viral DNA. This established PowerSoil as best for yield and quality of DNA from our samples, though ∼1/3 of the viral populations captured by each extraction kit were unique, suggesting appreciable differential biases among DNA extraction kits. Second, we evaluated the impact of purifying viral particles after resuspension (by cesium chloride gradients; CsCl) and of viral lysis method (heat vs bead-beating) on the resultant viromes. DNA yields after CsCl particle-purification were largely non-detectable, while unpurified samples yielded 1–2-fold more DNA after lysis by heat than by bead-beating. Virome quality was assessed by the number and size of metagenome-assembled viral contigs, which showed no increase after CsCl-purification, but did from heat lysis relative to bead-beating. We also evaluated sample preparation protocols for ssDNA virus recovery. In both CsCl-purified and non-purified samples, ssDNA viruses were successfully recovered by using the Accel-NGS 1S Plus Library Kit. While ssDNA viruses were identified in all three soil types, none were identified in the samples that used bead-beating, suggesting this lysis method may impact recovery. Further, 13 ssDNA vOTUs were identified compared to 582 dsDNA vOTUs, and the ssDNA vOTUs only accounted for ∼4% of the assembled reads, implying dsDNA viruses were dominant in these samples. This optimized approach was combined with the previously published viral resuspension protocol into a sample-to-virome protocol for soils now available at protocols.io, where community feedback creates ‘living’ protocols. This collective approach will be particularly valuable given the high physicochemical variability of soils, which will may require considerable soil type-specific optimization. This optimized protocol provides a starting place for developing quantitatively-amplified viromic datasets and will help enable viral ecogenomic studies on organic-rich soils.

Beyond cells - The virome in the human holobiont

Viromics, or viral metagenomics, is a relatively new and burgeoning field of research that studies the complete collection of viruses forming part of the microbiota in any given niche. It has strong foundations rooted in over a century of discoveries in the field of virology and recent advances in molecular biology and sequencing technologies. Historically, most studies have deconstructed the concept of viruses into a simplified perception of viral agents as mere pathogens, which demerits the scope of large-scale viromic analyses. Viruses are, in fact, much more than regular parasites. They are by far the most dynamic and abundant entity and the greatest killers on the planet, as well as the most effective geo-transforming genetic engineers and resource recyclers, acting on all life strata in any habitat. Yet, most of this uncanny viral world remains vastly unexplored to date, greatly hindered by the bewildering complexity inherent to such studies and the methodological and conceptual limitations. Viromic studies are just starting to address some of these issues but they still lag behind microbial metagenomics. In recent years, however, higher-throughput analysis and resequencing have rekindled interest in a field that is just starting to show its true potential. In this review, we take a look at the scientific and technological developments that led to the advent of viral and bacterial metagenomics with a particular, but not exclusive, focus on human viromics from an ecological perspective. We also address some of the most relevant challenges that current viral studies face and ponder on the future directions of the field.

Mouse Vendor Influence on the Bacterial and Viral Gut Composition Exceeds the Effect of Diet

Often physiological studies using mice from one vendor show different outcome when being reproduced using mice from another vendor. These divergent phenotypes between similar mouse strains from different vendors have been assigned to differences in the gut microbiome. During recent years, evidence has mounted that the gut viral community plays a key role in shaping the gut microbiome and may thus also influence mouse phenotype. However, to date inter-vendor variation in the murine gut virome has not been studied. Using a metavirome approach, combined with 16S rRNA gene sequencing, we here compare the composition of the viral and bacterial gut community of C57BL/6N mice from three different vendors exposed to either a chow-based low-fat diet or high-fat diet. Interestingly, both the bacterial and the viral component of the gut community differed significantly between vendors. The different diets also strongly influenced both the viral and bacterial gut community, but surprisingly the effect of vendor exceeded the effect of diet. In conclusion, the vendor effect is substantial not only on the gut bacterial community but also strongly influences viral community composition. Given the effect of GM on mice phenotype, this is essential to consider for increasing reproducibility of mouse studies.

Optimizing de novo genome assembly from PCR-amplified metagenomes

BACKGROUND

Metagenomics has transformed our understanding of microbial diversity across ecosystems, with recent advances enabling de novo assembly of genomes from metagenomes. These metagenome-assembled genomes are critical to provide ecological, evolutionary, and metabolic context for all the microbes and viruses yet to be cultivated. Metagenomes can now be generated from nanogram to subnanogram amounts of DNA. However, these libraries require several rounds of PCR amplification before sequencing, and recent data suggest these typically yield smaller and more fragmented assemblies than regular metagenomes.

METHODS

Here we evaluate de novo assembly methods of 169 PCR-amplified metagenomes, including 25 for which an unamplified counterpart is available, to optimize specific assembly approaches for PCR-amplified libraries. We first evaluated coverage bias by mapping reads from PCR-amplified metagenomes onto reference contigs obtained from unamplified metagenomes of the same samples. Then, we compared different assembly pipelines in terms of assembly size (number of bp in contigs ≥ 10 kb) and error rates to evaluate which are the best suited for PCR-amplified metagenomes.

RESULTS

Read mapping analyses revealed that the depth of coverage within individual genomes is significantly more uneven in PCR-amplified datasets versus unamplified metagenomes, with regions of high depth of coverage enriched in short inserts. This enrichment scales with the number of PCR cycles performed, and is presumably due to preferential amplification of short inserts. Standard assembly pipelines are confounded by this type of coverage unevenness, so we evaluated other assembly options to mitigate these issues. We found that a pipeline combining read deduplication and an assembly algorithm originally designed to recover genomes from libraries generated after whole genome amplification (single-cell SPAdes) frequently improved assembly of contigs ≥10 kb by 10 to 100-fold for low input metagenomes.

CONCLUSIONS

PCR-amplified metagenomes have enabled scientists to explore communities traditionally challenging to describe, including some with extremely low biomass or from which DNA is particularly difficult to extract. Here we show that a modified assembly pipeline can lead to an improved de novo genome assembly from PCR-amplified datasets, and enables a better genome recovery from low input metagenomes.

Archaeal Sources of Intact Membrane Lipid Biomarkers in the Oxygen Deficient Zone of the Eastern Tropical South Pacific

Archaea are ubiquitous in the modern ocean where they are involved in the carbon and nitrogen biogeochemical cycles. However, the majority of Archaea remain uncultured. Archaeal specific membrane intact polar lipids (IPLs) are biomarkers of the presence and abundance of living cells. They comprise archaeol and glycerol dibiphytanyl glycerol tetraethers (GDGTs) attached to various polar headgroups. However, little is known of the IPLs of uncultured marine Archaea, complicating their use as biomarkers. Here, we analyzed suspended particulate matter (SPM) obtained in high depth resolution from a coastal and open ocean site in the eastern tropical South Pacific (ETSP) oxygen deficient zone (ODZ) with the aim of determining possible biological sources of archaeal IPL by comparing their composition by Ultra High Pressure Liquid Chromatography coupled to high resolution mass spectrometry with the archaeal diversity by 16S rRNA gene amplicon sequencing and their abundance by quantitative PCR. Thaumarchaeotal Marine Group I (MGI) closely related to Ca. Nitrosopelagicus and Nitrosopumilus dominated the oxic surface and upper ODZ water together with Marine Euryarchaeota Group II (MGII). High relative abundance of hexose phosphohexose- (HPH) crenarchaeol, the specific biomarker for living Thaumarchaeota, and HPH-GDGT-0, dihexose- (DH) GDGT-3 and -4 were detected in these water masses. Within the ODZ, DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaea) of the Woesearchaeota DHVE-6 group and Marine Euryarchaeota Group III (MGIII) were present together with a higher proportion of archaeol-based IPLs, which were likely made by MGIII, since DPANN archaea are supposedly unable to synthesize their own IPLs and possibly have a symbiotic or parasitic partnership with MGIII. Finally, in deep suboxic/oxic waters a different MGI population occurred with HPH-GDGT-1, -2 and DH-GDGT-0 and -crenarchaeol, indicating that here MGI synthesize membranes with IPLs in a different relative abundance which could be attributed to the different detected population or to an environmental adaptation. Our study sheds light on the complex archaeal community of one of the most prominent ODZs and on the IPL biomarkers they potentially synthesize.

Human cytomegalovirus genomics and transcriptomics through the lens of next-generation sequencing: revision and future challenges

The human cytomegalovirus (HCMV) genome was sequenced by hierarchical shotgun almost 30 years ago. Over these years, low and high passaged strains have been sequenced, improving, albeit still far from complete, the understanding of the coding potential, expression dynamics and diversity of wild-type HCMV strains. Next-generation sequencing (NGS) platforms have enabled a huge advancement, facilitating the comparison of differentially passaged strains, challenging diagnostics and research based on a single or reduced gene set genotyping. In addition, it allowed to link genetic features to different viral phenotypes as for example, correlating large genomic re-arrangements to viral attenuation or different mutations to antiviral resistance and cell tropism. NGS platforms provided the first high-resolution experiments to HCMV dynamics, allowing the study of intra-host viral population structures and the description of rare transcriptional events. Long-read sequencing has recently become available, helping to identify new genomic re-arrangements, partially accounting for the genetic variability displayed in clinical isolates, as well as, in changing the understanding of the HCMV transcriptome. Better knowledge of the transcriptome resulted in a vast number of new splicing events and alternative transcripts, although most of them still need additional validation. This review summarizes the sequencing efforts reached so far, discussing its approaches and providing a revision and new nuances on HCMV sequence variability in the sequencing field.

A combined lipidomic and 16S rRNA gene amplicon sequencing approach reveals archaeal sources of intact polar lipids in the stratified Black Sea water column

Archaea are important players in marine biogeochemical cycles, and their membrane lipids are useful biomarkers in environmental and geobiological studies. However, many archaeal groups remain uncultured and their lipid composition unknown. Here, we aim to expand the knowledge on archaeal lipid biomarkers and determine the potential sources of those lipids in the water column of the euxinic Black Sea. The archaeal community was evaluated by 16S rRNA gene amplicon sequencing and by quantitative PCR. The archaeal intact polar lipids (IPLs) were investigated by ultra-high-pressure liquid chromatography coupled to high-resolution mass spectrometry. Our study revealed both a complex archaeal community and large changes with water depth in the IPL assemblages. In the oxic/upper suboxic waters (<105 m), the archaeal community was dominated by marine group (MG) I Thaumarchaeota, coinciding with a higher relative abundance of hexose phosphohexose crenarchaeol, a known marker for Thaumarchaeota. In the suboxic waters (80-110 m), MGI Nitrosopumilus sp. dominated and produced predominantly monohexose glycerol dibiphytanyl glycerol tetraethers (GDGTs) and hydroxy-GDGTs. Two clades of MGII Euryarchaeota were present in the oxic and upper suboxic zones in much lower abundances, preventing the detection of their specific IPLs. In the deep sulfidic waters (>110 m), archaea belonging to the DPANN Woesearchaeota, Bathyarchaeota, and ANME-1b clades dominated. Correlation analyses suggest that the IPLs GDGT-0, GDGT-1, and GDGT-2 with two phosphatidylglycerol (PG) head groups and archaeol with a PG, phosphatidylethanolamine, and phosphatidylserine head groups were produced by ANME-1b archaea. Bathyarchaeota represented 55% of the archaea in the deeper part of the euxinic zone and likely produces archaeol with phospho-dihexose and hexose-glucuronic acid head groups.

Viral Diversity and Its Relationship With Environmental Factors at the Surface and Deep Sea of Prydz Bay, Antarctica

A viral metagenomic analysis of five surface and two bottom water (878 meters below surface, mbs, and 3,357 mbs) samples from Prydz Bay, was conducted during February-March 2015. The results demonstrated that most of the DNA viruses were dsDNA viruses (79.73-94.06%, except at PBI1, 37.51%). Of these, Caudovirales (Siphoviridae, Myoviridae, and Podoviridae) phages were most abundant in surface seawater (67.67-71.99%), while nucleocytoplasmic large DNA viruses (NCLDVs) (Phycodnaviridae, Mimiviridae, and Pandoraviridae accounted for >30% of dsDNA viruses) were most abundant in the bottom water (3,357 mbs). Of the ssDNA viruses, Microviridae was the dominant family in PBI2, PBI3, PBOs, and PBI4b (57.09-87.55%), while Inoviridae (58.16%) was the dominant family in PBI1. Cellulophaga phages (phi38:1 and phi10:1) and Flavobacterium phage 11b, infecting the possible host strains affiliated with the family Flavobacteriaceae of the phylum Bacteroidetes, were abundant in surface water dsDNA viromes. The long contig (PBI2_1_C) from the viral metagenomes were most similar to the genome architectures of Cellulophaga phage phi10:1 and Flavobacterium phage 11b from the Arctic Ocean. Comparative analysis showed that the surface viral community of Prydz Bay could be clearly separated from other marine and freshwater environments. The deep sea viral community was similar to the deep sea viral metagenome at A Long-term Oligotrophic Habitat Assessment Station (ALOHA, at 22°45'N, 158°00'W). The multivariable analysis indicated that nutrients probably played an important role in shaping the local viral community structure. This study revealed the preliminary characteristics of the viral community in Prydz Bay, from both the surface and the deep sea. It provided evidence of the relationships between the virome and the environment in Prydz Bay and provided the first data from the deep sea viral community of the Southern Ocean.

The smallest ssDNA phage infecting a marine bacterium

In the marine environment, only a few lytic single-stranded DNA (ssDNA) phages have been isolated and characterized, despite the fact that diverse ssDNA bacteriophages have been discovered via metagenomic studies. In this study, we isolated and characterized a new ssDNA phage, vB_RpoMi-Mini, which infects a marine bacterium Ruegeria pomeroyi DSS-3. With a genome size of 4248 bp and only four putative open reading frames (ORF), vB_RpoMi-Mini becomes the smallest ssDNA phage among the known ssDNA phage isolates and represents the DNA bacteriophage with the least number of ORFs. Genome-wide analysis reveals that bacteriophage Mini is distantly related to the known ssDNA phages and belongs to an unclassified ssDNA phage within the Microviridae family. The presence of peptidase in vB_RpoMi-Mini genome further implies that horizontal gene transfer could be an important driving force in the evolution of ssDNA phages. Bacteriophage Mini seems to have lost the spike protein commonly seen in ssDNA phages, suggesting that ssDNA phage can be more diverse than previously thought. Metagenomic analysis indicates that Mini-like phages are widely distributed in the environments. The discovery of vB_RpoMi-Mini expands our understanding of ssDNA phages in nature, and also indicates our dearth of knowledge regarding of ssDNA phages.

Soil Viruses Are Underexplored Players in Ecosystem Carbon Processing

Rapidly thawing permafrost harbors ∼30 to 50% of global soil carbon, and the fate of this carbon remains unknown. Microorganisms will play a central role in its fate, and their viruses could modulate that impact via induced mortality and metabolic controls. Because of the challenges of recovering viruses from soils, little is known about soil viruses or their role(s) in microbial biogeochemical cycling. Here, we describe 53 viral populations (viral operational taxonomic units [vOTUs]) recovered from seven quantitatively derived (i.e., not multiple-displacement-amplified) viral-particle metagenomes (viromes) along a permafrost thaw gradient at the Stordalen Mire field site in northern Sweden. Only 15% of these vOTUs had genetic similarity to publicly available viruses in the RefSeq database, and ∼30% of the genes could be annotated, supporting the concept of soils as reservoirs of substantial undescribed viral genetic diversity. The vOTUs exhibited distinct ecology, with different distributions along the thaw gradient habitats, and a shift from soil-virus-like assemblages in the dry palsas to aquatic-virus-like assemblages in the inundated fen. Seventeen vOTUs were linked to microbial hosts (in silico), implicating viruses in infecting abundant microbial lineages from Acidobacteria, Verrucomicrobia, and Deltaproteobacteria, including those encoding key biogeochemical functions such as organic matter degradation. Thirty auxiliary metabolic genes (AMGs) were identified and suggested virus-mediated modulation of central carbon metabolism, soil organic matter degradation, polysaccharide binding, and regulation of sporulation. Together, these findings suggest that these soil viruses have distinct ecology, impact host-mediated biogeochemistry, and likely impact ecosystem function in the rapidly changing Arctic. IMPORTANCE This work is part of a 10-year project to examine thawing permafrost peatlands and is the first virome-particle-based approach to characterize viruses in these systems. This method yielded >2-fold-more viral populations (vOTUs) per gigabase of metagenome than vOTUs derived from bulk-soil metagenomes from the same site (J. B. Emerson, S. Roux, J. R. Brum, B. Bolduc, et al., Nat Microbiol 3:870-880, 2018, https://doi.org/10.1038/s41564-018-0190-y). We compared the ecology of the recovered vOTUs along a permafrost thaw gradient and found (i) habitat specificity, (ii) a shift in viral community identity from soil-like to aquatic-like viruses, (iii) infection of dominant microbial hosts, and (iv) carriage of host metabolic genes. These vOTUs can impact ecosystem carbon processing via top-down (inferred from lysing dominant microbial hosts) and bottom-up (inferred from carriage of auxiliary metabolic genes) controls. This work serves as a foundation which future studies can build upon to increase our understanding of the soil virosphere and how viruses affect soil ecosystem services.

Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses

BACKGROUND

Viruses are key players regulating microbial ecosystems. Exploration of viral assemblages is now possible thanks to the development of metagenomics, the most powerful tool available for studying viral ecology and discovering new viruses. Unfortunately, several sources of bias lead to the misrepresentation of certain viruses within metagenomics workflows, hindering the shift from merely descriptive studies towards quantitative comparisons of communities. Therefore, benchmark studies on virus enrichment and random amplification protocols are required to better understand the sources of bias.

RESULTS

We assessed the bias introduced by viral enrichment on mock assemblages composed of seven DNA viruses, and the bias from random amplification methods on human saliva DNA viromes, using qPCR and deep sequencing, respectively. While iodixanol cushions and 0.45 μm filtration preserved the original composition of nuclease-protected viral genomes, low-force centrifugation and 0.22 μm filtration removed large viruses. Comparison of unamplified and randomly amplified saliva viromes revealed that multiple displacement amplification (MDA) induced stochastic bias from picograms of DNA template. However, the type of bias shifted to systematic using 1 ng, with only a marginal influence by amplification time. Systematic bias consisted of over-amplification of small circular genomes, and under-amplification of those with extreme GC content, a negative bias that was shared with the PCR-based sequence-independent, single-primer amplification (SISPA) method. MDA based on random priming provided by a DNA primase activity slightly outperformed those based on random hexamers and SISPA, which may reflect differences in ability to handle sequences with extreme GC content. SISPA viromes showed uneven coverage profiles, with high coverage peaks in regions with low linguistic sequence complexity. Despite misrepresentation of certain viruses after random amplification, ordination plots based on dissimilarities among contig profiles showed perfect overlapping of related amplified and unamplified saliva viromes and strong separation from unrelated saliva viromes. This result suggests that random amplification bias has a minor impact on beta diversity studies.

CONCLUSIONS

Benchmark analyses of mock and natural communities of viruses improve understanding and mitigate bias in metagenomics surveys. Bias induced by random amplification methods has only a minor impact on beta diversity studies of human saliva viromes.

Bioreactor virome metagenomics sequencing using DNA spike-ins

With the emergence of Next Generation Sequencing, major advances were made with regard to identifying viruses in natural environments. However, bioinformatical research on viruses is still limited because of the low amounts of viral DNA that can be obtained for analysis. To overcome this limitation, DNA is often amplified with multiple displacement amplification (MDA), which may cause an unavoidable bias. Here, we describe a case study in which the virome of a bioreactor is sequenced using Ion Torrent technology. DNA-spiking of samples is compared with MDA-amplified samples. DNA for spiking was obtained by amplifying a bacterial 16S rRNA gene. After sequencing, the 16S rRNA gene reads were removed by mapping to the Silva database. Three samples were tested, a whole genome from Enterobacteria P1 Phage and two viral metagenomes from an infected bioreactor. For one sample, the new DNA-spiking protocol was compared with the MDA technique. When MDA was applied, the overall GC content of the reads showed a bias towards lower GC%, indicating a change in composition of the DNA sample. Assemblies using all available reads from both MDA and the DNA-spiked samples resulted in six viral genomes. All six genomes could be almost completely retrieved (97.9%–100%) when mapping the reads from the DNA-spiked sample to those six genomes. In contrast, 6.3%–77.7% of three viral genomes was covered by reads obtained using the MDA amplification method and only three were nearly fully covered (97.4%–100%). This case study shows that DNA-spiking could be a simple and inexpensive alternative with very low bias for sequencing of metagenomes for which low amounts of DNA are available.

Insights into the metabolic functioning of a multipartner ciliate symbiosis from oxygen-depleted sediments

Symbioses between anaerobic or microaerophilic protists and prokaryotes are common in anoxic and oxygen-depleted habitats ranging from marine sediments to gastrointestinal tracts. Nevertheless, little is known about the mechanisms of metabolic interaction between partners. In these putatively syntrophic associations, consumption of fermentative end products (e.g., hydrogen) by the prokaryotic symbionts is thought to facilitate protistan anaerobic metabolism. Here, we employed metagenomic and metatranscriptomic sequencing of a microaerophilic or anaerobic karyorelictid ciliate and its prokaryotic symbionts from oxygen-depleted Santa Barbara Basin (CA, USA) sediments to assess metabolic coupling within this consortium. This sequencing confirmed the predominance of deltaproteobacterial symbionts from the Families Desulfobacteraceae and Desulfobulbaceae and suggested active symbiont reduction of host-provided sulphate, transfer of small organic molecules from host to symbionts and hydrogen cycling among the symbionts. In addition, patterns of gene expression indicated active cell division by the symbionts, their growth via autotrophic processes and nitrogen exchange with the ciliate host. Altogether, this research underscores the importance of symbiont metabolism to host fermentative metabolism and, thus, likely its success in anoxic and low-oxygen habitats, but also suggests ciliate-associated prokaryotes play a role in important biogeochemical processes.

Sequencing of Genomes from Environmental Single Cells

Sequencing of single bacterial and archaeal cells is an important methodology that provides access to the genetic makeup of uncultivated microorganisms. We here describe the high-throughput fluorescence-activated cell sorting-based isolation of single cells from the environment, their lysis and strand displacement-mediated whole genome amplification. We further outline 16S rRNA gene sequence-based screening of single-cell amplification products, their preparation for Illumina sequencing libraries, and finally propose computational methods for read and contig level quality control of the resulting sequence data.

Metagenomic assessment of the interplay between the environment and the genetic diversification of Acinetobacter

Most bacteria have poorly characterized environmental reservoirs and unknown closely related species. This hampers the study of bacterial evolutionary ecology because both the environment and the genetic background of ancestral lineages are unknown. We combined metagenomics, comparative genomics and phylogenomics to overcome this limitation, to identify novel taxa and to propose environments where they can be isolated. We applied this method to characterize the ecological distribution of known and novel lineages of Acinetobacter spp. We observed two major environmental transitions at deep phylogenetic levels, splitting the genus into three ecologically differentiated clades. One of these has rapidly shifted towards host‐association by acquiring genes involved in bacteria‐eukaryote interactions. We show that environmental perturbations affect species distribution in predictable ways: bovines have very diverse communities of Acinetobacter, unless they were administered antibiotics, in which case they show highly uniform communities of Acinetobacter spp. that resemble those of humans. Our results uncover the diversity of bacterial lineages, overpassing the limitations of classical cultivation methods and highlight the role of the environment in shaping their evolution.