Validation of picogram- and femtogram-input DNA libraries for microscale metagenomics

Metagenomic insights into taxonomic and functional patterns in shallow coastal and deep subseafloor sediments in the Western Pacific

Marine sediments are vast, underexplored habitats and represent one of the largest carbon deposits on our planet. Microbial communities drive nutrient cycling in these sediments, but the full extent of their taxonomic and metabolic diversity remains to be explored. Here, we analysed shallow coastal and deep subseafloor sediment cores from 0.01 to nearly 600 metres below the seafloor, in the Western Pacific Region. Applying metagenomics, we identified several taxonomic clusters across all samples, which mainly aligned with depth and sediment type. Inferring functional patterns provided insights into possible ecological roles of the main microbial taxa. These included Chloroflexota, the most abundant phylum across all samples, whereby the classes Dehalococcoida and Anaerolineae dominated deep-subsurface and most shallow coastal sediments, respectively. Thermoproteota and Asgardarchaeota were the most abundant phyla among Archaea, contributing to high relative abundances of Archaea reaching over 50% in some samples. We recovered high-quality metagenome-assembled genomes for all main prokaryotic lineages and proposed names for three phyla, i.e. Tangaroaeota phyl. nov. (former RBG-13-66-14), Ryujiniota phyl. nov. (former UBA6262) and Spongiamicota phyl. nov. (former UBA8248). Metabolic capabilities across all samples ranged from aerobic respiration and photosynthesis in the shallowest sediment layers to heterotrophic carbon utilization, sulphate reduction and methanogenesis in deeper anoxic sediments. We also identified taxa with the potential to be involved in nitrogen and sulphur cycling and heterotrophic carbon utilization. In summary, this study contributes to our understanding of the taxonomic and functional diversity in benthic prokaryotic communities across marine sediments in the Western Pacific Region.

Metagenomic study of the microbiome and key geochemical potentials associated with architectural heritage sites: a case study of the Song Dynasty city wall in Shou County, China

Historical cultural heritage sites are valuable for all of mankind, as they reflect the material and spiritual wealth of by nations, countries, or specific groups during the development of human civilization. The types and functions of microorganisms that form biofilms on the surfaces of architectural heritage sites influence measures to preserve and protect these sites. These microorganisms contribute to the biocorrosion of architectural heritage structures through the cycling of chemical elements. The ancient city wall of Shou County is a famous architectural and cultural heritage site from China's Song Dynasty, and the protection and study of this site have substantial historical and cultural significance. In this study, we used metagenomics to study the microbial diversity and taxonomic composition of the Song Dynasty city wall in Shou County, a tangible example of Chinese cultural heritage. The study covered three main topics: (1) examining the distribution of bacteria in the biofilm on the surfaces of the Song Dynasty city wall in Shou County; (2) predicting the influence of bacteria involved in the C, N, and S cycles on the corrosion of the city wall via functional gene analysis; and (3) discussing cultural heritage site protection measures for biocorrosion-related bacteria to investigate the impact of biocorrosion on the Song Dynasty city wall in Shou County, a tangible example of Chinese cultural heritage. The study revealed that (1) the biofilm bacteria mainly belonged to Proteobacteria, Actinobacteria, Cyanobacteria, Bacteroidetes, and Firmicutes, which accounted for more than 70% of the total bacteria in the biofilms. The proportion of fungi in the microbial community of the well-preserved city wall was greater than that in the damaged city wall. The proportion of archaea was low-less than 1%. (2) According to the Shannon diversity index, the well-preserved portion of the ancient city wall had the highest diversity of bacteria, fungi, and archaea, and bacterial diversity on the good city wall was greater than that on the corroded city wall. (3) Bray-Curtis distances revealed that the genomes of the two good city walls were similar and that the genomes of the corroded city wall portions were similar. Researchers also detected human intestine-related bacteria in four locations on the city walls, with the proportion of these bacteria in the microbial community being greater on good city walls than on bad city walls. (4) KEGG functional analysis revealed that the energy metabolism and inorganic ion transport activities of the bacterial community on the corroded city wall were greater than those of the good city wall. (5) In the carbon cycle, the absence of active glycolysis, the ED pathway, and the TCA cycle played significant roles in the collapse of the east city wall. (6) The nitrogen cycling processes involved ammonia oxidation and nitrite reduction to nitrate. (7) In the sulfur cycle, researchers discovered a crucial differential functional gene, SoxY, which facilitates the conversion of thiosulfate to sulfate. This study suggests that, in the future, biological approaches can be used to help cultural heritage site protectors achieve targeted and precise protection of cultural relics through the use of microbial growth inhibition technology. The results of this study serve as a guide for the protection of cultural heritage sites in other parts of China and provide a useful supplement to research on the protection of world cultural heritage or architectural heritage sites.

Metagenomic airborne resistome from urban hot spots through the One Health lens

Human activities are a significant contributor to the spread of antibiotic resistance genes (ARGs), which pose a serious threat to human health. These ARGs can be transmitted through various pathways, including air, within the context of One Health. This study used metagenomics to monitor the resistomes in urban air from two critical locations: a wastewater treatment plant and a hospital, both indoor and outdoor. The presence of cell-like structures was confirmed through fluorescence microscopy. The metagenomic analysis revealed a wide variety of ARGs and a high diversity of antibiotic-resistant bacteria in the airborne particles collected. The wastewater treatment plant showed higher relative abundances with 32 ARG hits per Gb and m3, followed by the main entrance of the hospital (indoor) with ≈5 ARG hits per Gb and m3. The hospital entrance exhibited the highest ARG richness, with a total of 152 different ARGs classified into nine categories of antibiotic resistance. Common commensal and pathogenic bacteria carrying ARGs, such as Moraxella, Staphylococcus and Micrococcus, were detected in the indoor airborne particles of the hospital. Interestingly, no ARGs were shared among all the samples analysed, indicating a highly variable dynamic of airborne resistomes. Furthermore, the study found no ARGs in the airborne viral fractions analysed, suggesting that airborne viruses play a negligible role in the dissemination of ARGs.

A genome-centric view of the role of the Acropora kenti microbiome in coral health and resilience

Microbial diversity has been extensively explored in reef-building corals. However, the functional roles of coral-associated microorganisms remain poorly elucidated. Here, we recover 191 bacterial and 10 archaeal metagenome-assembled genomes (MAGs) from the coral Acropora kenti (formerly A. tenuis) and adjacent seawater, to identify microbial functions and metabolic interactions within the holobiont. We show that 82 MAGs were specific to the A. kenti holobiont, including members of the Pseudomonadota, Bacteroidota, and Desulfobacterota. A. kenti-specific MAGs displayed significant differences in their genomic features and functional potential relative to seawater-specific MAGs, with a higher prevalence of genes involved in host immune system evasion, nitrogen and carbon fixation, and synthesis of five essential B-vitamins. We find a diversity of A. kenti-specific MAGs encode the biosynthesis of essential amino acids, such as tryptophan, histidine, and lysine, which cannot be de novo synthesised by the host or Symbiodiniaceae. Across a water quality gradient spanning 2° of latitude, A. kenti microbial community composition is correlated to increased temperature and dissolved inorganic nitrogen, with corresponding enrichment in molecular chaperones, nitrate reductases, and a heat-shock protein. We reveal mechanisms of A. kenti-microbiome-symbiosis on the Great Barrier Reef, highlighting the interactions underpinning the health of this keystone holobiont.

Strong chemotaxis by marine bacteria towards polysaccharides is enhanced by the abundant organosulfur compound DMSP

The ability of marine bacteria to direct their movement in response to chemical gradients influences inter-species interactions, nutrient turnover, and ecosystem productivity. While many bacteria are chemotactic towards small metabolites, marine organic matter is predominantly composed of large molecules and polymers. Yet, the signalling role of these large molecules is largely unknown. Using in situ and laboratory-based chemotaxis assays, we show that marine bacteria are strongly attracted to the abundant algal polysaccharides laminarin and alginate. Unexpectedly, these polysaccharides elicited stronger chemoattraction than their oligo- and monosaccharide constituents. Furthermore, chemotaxis towards laminarin was strongly enhanced by dimethylsulfoniopropionate (DMSP), another ubiquitous algal-derived metabolite. Our results indicate that DMSP acts as a methyl donor for marine bacteria, increasing their gradient detection capacity and facilitating their access to polysaccharide patches. We demonstrate that marine bacteria are capable of strong chemotaxis towards large soluble polysaccharides and uncover a new ecological role for DMSP in enhancing this attraction. These navigation behaviours may contribute to the rapid turnover of polymers in the ocean, with important consequences for marine carbon cycling.

Sequencing-based analysis of the microbiomes of Spanish food processing facilities reveals environment-specific variation in the dominant taxa and antibiotic resistance genes

In the last years, advances in high throughput sequencing technologies have opened the possibility to broaden environmental monitoring activities in facilities processing food, offering expanded opportunities for characterizing in an untargeted manner the microbiome and resistome of foods and food processing environments (FPE) with huge potential benefits in food safety management systems. Here the microbiome and resistome of FPE from slaughterhouses (n = 3), dairy (n = 12) and meat (n = 10) processing plants were assessed through whole metagenome sequencing of 2 composite samples for each facility, comprising 10 FPE swabs taken from food contact surfaces and 10 FPE samples from non-food contact surfaces, respectively. FPE from slaughterhouses had more diverse microbiomes and resistomes, while FPE from dairy processing plants showed the highest β-dispersion, consistent with a more heterogeneous microbiome and resistome composition. The predominant bacterial genera depended on the industry type, with Pseudomonas and Psychrobacter being highly dominant in surfaces from slaughterhouses and meat industries, while different lactic acid bacteria predominated in dairy industries. The most abundant antimicrobial resistance genes (ARG) found were associated with resistance to aminoglycosides, tetracyclines and quaternary ammonium compounds (QAC). ARGs relating to resistance to aminoglycosides and tetracyclines were significantly more prevalent in slaughterhouses than in food processing plants, while QAC resistance genes were particularly abundant in some food contact surfaces from dairy and meat processing plants, suggesting that daily sanitation under suboptimal conditions may be selecting for persistent microbiota tolerant to these biocides in some facilities. The taxonomic mapping of ARG pointed to specific bacterial genera, such as Escherichia, Bacillus, or Staphylococcus, as carriers of the most relevant resistance determinants. About 63% of all ARG reads were assigned to contigs classified as plasmid-associated, indicating that the resistome of FPE may be strongly shaped through the spread of mobile genetic elements. Overall, the relevance of FPE as reservoirs of ARG was confirmed and it was demonstrated that next generation sequencing technologies allowing a deep characterisation of sources and routes of spread of microorganisms and antimicrobial resistance determinants in food industry settings hold promise to be integrated in monitoring and food safety management programmes.

Microbiota and Resistome Analysis of Colostrum and Milk from Dairy Cows Treated with and without Dry Cow Therapies

This study investigated the longitudinal impact of methods for the drying off of cows with and without dry cow therapy (DCT) on the microbiota and resistome profile in colostrum and milk samples from cows. Three groups of healthy dairy cows (n = 24) with different antibiotic treatments during DCT were studied. Colostrum and milk samples from Month 0 (M0), 2 (M2), 4 (M4) and 6 (M6) were analysed using whole-genome shotgun-sequencing. The microbial diversity from antibiotic-treated groups was different and higher than that of the non-antibiotic group. This difference was more evident in milk compared to colostrum, with increasing diversity seen only in antibiotic-treated groups. The microbiome of antibiotic-treated groups clustered separately from the non-antibiotic group at M2-, M4- and M6 milk samples, showing the effect of antibiotic treatment on between-group (beta) diversity. The non-antibiotic group did not show a high relative abundance of mastitis-causing pathogens during early lactation and was more associated with genera such as Psychrobacter, Serratia, Gordonibacter and Brevibacterium. A high relative abundance of antibiotic resistance genes (ARGs) was observed in the milk of antibiotic-treated groups with the Cephaguard group showing a significantly high abundance of genes conferring resistance to cephalosporin, aminoglycoside and penam classes. The data support the use of non-antibiotic alternatives for drying off in cows.

Metagenomic Methods for Addressing NASA's Planetary Protection Policy Requirements on Future Missions: A Workshop Report

Molecular biology methods and technologies have advanced substantially over the past decade. These new molecular methods should be incorporated among the standard tools of planetary protection (PP) and could be validated for incorporation by 2026. To address the feasibility of applying modern molecular techniques to such an application, NASA conducted a technology workshop with private industry partners, academics, and government agency stakeholders, along with NASA staff and contractors. The technical discussions and presentations of the Multi-Mission Metagenomics Technology Development Workshop focused on modernizing and supplementing the current PP assays. The goals of the workshop were to assess the state of metagenomics and other advanced molecular techniques in the context of providing a validated framework to supplement the bacterial endospore-based NASA Standard Assay and to identify knowledge and technology gaps. In particular, workshop participants were tasked with discussing metagenomics as a stand-alone technology to provide rapid and comprehensive analysis of total nucleic acids and viable microorganisms on spacecraft surfaces, thereby allowing for the development of tailored and cost-effective microbial reduction plans for each hardware item on a spacecraft. Workshop participants recommended metagenomics approaches as the only data source that can adequately feed into quantitative microbial risk assessment models for evaluating the risk of forward (exploring extraterrestrial planet) and back (Earth harmful biological) contamination. Participants were unanimous that a metagenomics workflow, in tandem with rapid targeted quantitative (digital) PCR, represents a revolutionary advance over existing methods for the assessment of microbial bioburden on spacecraft surfaces. The workshop highlighted low biomass sampling, reagent contamination, and inconsistent bioinformatics data analysis as key areas for technology development. Finally, it was concluded that implementing metagenomics as an additional workflow for addressing concerns of NASA's robotic mission will represent a dramatic improvement in technology advancement for PP and will benefit future missions where mission success is affected by backward and forward contamination.

Novel complete methanogenic pathways in longitudinal genomic study of monogastric age-associated archaea

BACKGROUND

Archaea perform critical roles in the microbiome system, including utilizing hydrogen to allow for enhanced microbiome member growth and influencing overall host health. With the majority of microbiome research focusing on bacteria, the functions of archaea are largely still under investigation. Understanding methanogenic functions during the host lifetime will add to the limited knowledge on archaeal influence on gut and host health. In our study, we determined lifelong archaea dynamics, including detection and methanogenic functions, while assessing global, temporal and host distribution of our novel archaeal metagenome-assembled genomes (MAGs). We followed 7 monogastric swine throughout their life, from birth to adult (1-156 days of age), and collected feces at 22 time points. The samples underwent gDNA extraction, Illumina sequencing, bioinformatic quality and assembly processes, MAG taxonomic assignment and functional annotation. MAGs were utilized in downstream phylogenetic analysis for global, temporal and host distribution in addition to methanogenic functional potential determination.

RESULTS

We generated 1130 non-redundant MAGs, representing 588 unique taxa at the species level, with 8 classified as methanogenic archaea. The taxonomic classifications were as follows: orders Methanomassiliicoccales (5) and Methanobacteriales (3); genera UBA71 (3), Methanomethylophilus (1), MX-02 (1), and Methanobrevibacter (3). We recovered the first US swine Methanobrevibacter UBA71 sp006954425 and Methanobrevibacter gottschalkii MAGs. The Methanobacteriales MAGs were identified primarily during the young, preweaned host whereas Methanomassiliicoccales primarily in the adult host. Moreover, we identified our methanogens in metagenomic sequences from Chinese swine, US adult humans, Mexican adult humans, Swedish adult humans, and paleontological humans, indicating that methanogens span different hosts, geography and time. We determined complete metabolic pathways for all three methanogenic pathways: hydrogenotrophic, methylotrophic, and acetoclastic. This study provided the first evidence of acetoclastic methanogenesis in archaea of monogastric hosts which indicated a previously unknown capability for acetate utilization in methanogenesis for monogastric methanogens. Overall, we hypothesized that the age-associated detection patterns were due to differential substrate availability via the host diet and microbial metabolism, and that these methanogenic functions are likely crucial to methanogens across hosts. This study provided a comprehensive, genome-centric investigation of monogastric-associated methanogens which will further improve our understanding of microbiome development and functions.

Back to Basics: A Simplified Improvement to Multiple Displacement Amplification for Microbial Single-Cell Genomics

Microbial single-cell genomics (SCG) provides access to the genomes of rare and uncultured microorganisms and is a complementary method to metagenomics. Due to the femtogram-levels of DNA in a single microbial cell, sequencing the genome requires whole genome amplification (WGA) as a preliminary step. However, the most common WGA method, multiple displacement amplification (MDA), is known to be costly and biased against specific genomic regions, preventing high-throughput applications and resulting in uneven genome coverage. Thus, obtaining high-quality genomes from many taxa, especially minority members of microbial communities, becomes difficult. Here, we present a volume reduction approach that significantly reduces costs while improving genome coverage and uniformity of DNA amplification products in standard 384-well plates. Our results demonstrate that further volume reduction in specialized and complex setups (e.g., microfluidic chips) is likely unnecessary to obtain higher-quality microbial genomes. This volume reduction method makes SCG more feasible for future studies, thus helping to broaden our knowledge on the diversity and function of understudied and uncharacterized microorganisms in the environment.

Chemotaxis may assist marine heterotrophic bacterial diazotrophs to find microzones suitable for N2 fixation in the pelagic ocean

Heterotrophic bacterial diazotrophs (HBDs) are ubiquitous in the pelagic ocean, where they have been predicted to carry out the anaerobic process of nitrogen fixation within low-oxygen microenvironments associated with marine pelagic particles. However, the mechanisms enabling particle colonization by HBDs are unknown. We hypothesized that HBDs use chemotaxis to locate and colonize suitable microenvironments, and showed that a cultivated marine HBD is chemotactic toward amino acids and phytoplankton-derived DOM. Using an in situ chemotaxis assay, we also discovered that diverse HBDs at a coastal site are motile and chemotactic toward DOM from various phytoplankton taxa and, indeed, that the proportion of diazotrophs was up to seven times higher among the motile fraction of the bacterial community compared to the bulk seawater community. Finally, three of four HBD isolates and 16 of 17 HBD metagenome assembled genomes, recovered from major ocean basins and locations along the Australian coast, each encoded >85% of proteins affiliated with the bacterial chemotaxis pathway. These results document the widespread capacity for chemotaxis in diverse and globally relevant marine HBDs. We suggest that HBDs could use chemotaxis to seek out and colonize low-oxygen microenvironments suitable for nitrogen fixation, such as those formed on marine particles. Chemotaxis in HBDs could therefore affect marine nitrogen and carbon biogeochemistry by facilitating nitrogen fixation within otherwise oxic waters, while also altering particle degradation and the efficiency of the biological pump.

Toward efficient and high-fidelity metagenomic data from sub-nanogram DNA: evaluation of library preparation and decontamination methods

BACKGROUND

Shotgun metagenomic sequencing has greatly expanded the understanding of microbial communities in various biological niches. However, it is still challenging to efficiently convert sub-nanogram DNA to high-quality metagenomic libraries and obtain high-fidelity data, hindering the exploration of niches with low microbial biomass.

RESULTS

To cope with this challenge comprehensively, we evaluated the performance of various library preparation methods on 0.5 pg-5 ng synthetic microbial community DNA, characterized contaminants, and further applied different in silico decontamination methods. First, we discovered that whole genome amplification prior to library construction led to worse outcomes than preparing libraries directly. Among different non-WGA-based library preparation methods, we found the endonuclease-based method being generally good for different amounts of template and the tagmentation-based method showing specific advantages with 0.5 pg template, based on evaluation metrics including fidelity, proportion of designated reads, and reproducibility. The load of contaminating DNA introduced by library preparation varied from 0.01 to 15.59 pg for different kits and accounted for 0.05 to 45.97% of total reads. A considerable fraction of the contaminating reads were mapped to human commensal and pathogenic microbes, thus potentially leading to erroneous conclusions in human microbiome studies. Furthermore, the best performing in silico decontamination method in our evaluation, Decontam-either, was capable of recovering the real microbial community from libraries where contaminants accounted for less than 10% of total reads, but not from libraries with heavy and highly varied contaminants.

CONCLUSIONS

This study demonstrates that high-quality metagenomic data can be obtained from samples with sub-nanogram microbial DNA by combining appropriate library preparation and in silico decontamination methods and provides a general reference for method selection for samples with varying microbial biomass.

Historical contingencies and phage induction diversify bacterioplankton communities at the microscale

In many natural environments, microorganisms decompose microscale resource patches made of complex organic matter. The growth and collapse of populations on these resource patches unfold within spatial ranges of a few hundred micrometers or less, making such microscale ecosystems hotspots of heterotrophic metabolism. Despite the potential importance of patch-level dynamics for the large-scale functioning of heterotrophic microbial communities, we have not yet been able to delineate the ecological processes that control natural populations at the microscale. Here, we address this challenge by characterizing the natural marine communities that assembled on over 1,000 individual microscale particles of chitin, the most abundant marine polysaccharide. Using low-template shotgun metagenomics and imaging, we find significant variation in microscale community composition despite the similarity in initial species pools across replicates. Chitin-degrading taxa that were rare in seawater established large populations on a subset of particles, resulting in a wide range of predicted chitinolytic abilities and biomass at the level of individual particles. We show, through a mathematical model, that this variability can be attributed to stochastic colonization and historical contingencies affecting the tempo of growth on particles. We find evidence that one biological process leading to such noisy growth across particles is differential predation by temperate bacteriophages of chitin-degrading strains, the keystone members of the community. Thus, initial stochasticity in assembly states on individual particles, amplified through ecological interactions, may have significant consequences for the diversity and functionality of systems of microscale patches.

Turnover in Life-Strategies Recapitulates Marine Microbial Succession Colonizing Model Particles

Particulate organic matter (POM) in the ocean sustains diverse communities of bacteria that mediate the remineralization of organic complex matter. However, the variability of these particles and of the environmental conditions surrounding them present a challenge to the study of the ecological processes shaping particle-associated communities and their function. In this work, we utilize data from experiments in which coastal water communities are grown on synthetic particles to ask which are the most important ecological drivers of their assembly and associated traits. Combining 16S rRNA amplicon sequencing with shotgun metagenomics, together with an analysis of the full genomes of a subset of isolated strains, we were able to identify two-to-three distinct community classes, corresponding to early vs. late colonizers. We show that these classes are shaped by environmental selection (early colonizers) and facilitation (late colonizers) and find distinctive traits associated with each class. While early colonizers have a larger proportion of genes related to the uptake of nutrients, motility, and environmental sensing with few pathways enriched for metabolism, late colonizers devote a higher proportion of genes for metabolism, comprising a wide array of different pathways including the metabolism of carbohydrates, amino acids, and xenobiotics. Analysis of selected pathways suggests the existence of a trophic-chain topology connecting both classes for nitrogen metabolism, potential exchange of branched chain amino acids for late colonizers, and differences in bacterial doubling times throughout the succession. The interpretation of these traits suggests a distinction between early and late colonizers analogous to other classifications found in the literature, and we discuss connections with the classical distinction between r- and K-strategists.

Introduction to the principles and methods underlying the recovery of metagenome-assembled genomes from metagenomic data

The rise of metagenomics offers a leap forward for understanding the genetic diversity of microorganisms in many different complex environments by providing a platform that can identify potentially unlimited numbers of known and novel microorganisms. As such, it is impossible to imagine new major initiatives without metagenomics. Nevertheless, it represents a relatively new discipline with various levels of complexity and demands on bioinformatics. The underlying principles and methods used in metagenomics are often seen as common knowledge and often not detailed or fragmented. Therefore, we reviewed these to guide microbiologists in taking the first steps into metagenomics. We specifically focus on a workflow aimed at reconstructing individual genomes, that is, metagenome-assembled genomes, integrating DNA sequencing, assembly, binning, identification and annotation.

Comprehensive human amniotic fluid metagenomics supports the sterile womb hypothesis

As metagenomic approaches for detecting infectious agents have improved, each tissue that was once thought to be sterile has been found to harbor a variety of microorganisms. Controversy still exists over the status of amniotic fluid, which is part of an immunologically privileged zone that is required to prevent maternal immune system rejection of the fetus. Due to this privilege, the exclusion of microbes has been proposed to be mandatory, leading to the sterile womb hypothesis. Since nucleic acid yields from amniotic fluid are very low, contaminating nucleic acid found in water, reagents and the laboratory environment frequently confound attempts to address this hypothesis. Here we present metagenomic criteria for microorganism detection and a metagenomic method able to be performed with small volumes of starting material, while controlling for exogenous contamination, to circumvent these and other pitfalls. We use this method to show that human mid-gestational amniotic fluid has no detectable virome or microbiome, supporting the sterile womb hypothesis.

Chemotaxis shapes the microscale organization of the ocean's microbiome

The capacity of planktonic marine microorganisms to actively seek out and exploit microscale chemical hotspots has been widely theorized to affect ocean-basin scale biogeochemistry^1-3, but has never been examined comprehensively in situ among natural microbial communities. Here, using a field-based microfluidic platform to quantify the behavioural responses of marine bacteria and archaea, we observed significant levels of chemotaxis towards microscale hotspots of phytoplankton-derived dissolved organic matter (DOM) at a coastal field site across multiple deployments, spanning several months. Microscale metagenomics revealed that a wide diversity of marine prokaryotes, spanning 27 bacterial and 2 archaeal phyla, displayed chemotaxis towards microscale patches of DOM derived from ten globally distributed phytoplankton species. The distinct DOM composition of each phytoplankton species attracted phylogenetically and functionally discrete populations of bacteria and archaea, with 54% of chemotactic prokaryotes displaying highly specific responses to the DOM derived from only one or two phytoplankton species. Prokaryotes exhibiting chemotaxis towards phytoplankton-derived compounds were significantly enriched in the capacity to transport and metabolize specific phytoplankton-derived chemicals, and displayed enrichment in functions conducive to symbiotic relationships, including genes involved in the production of siderophores, B vitamins and growth-promoting hormones. Our findings demonstrate that the swimming behaviour of natural prokaryotic assemblages is governed by specific chemical cues, which dictate important biogeochemical transformation processes and the establishment of ecological interactions that structure the base of the marine food web.

The Chronic Wound Phageome: Phage Diversity and Associations with Wounds and Healing Outcomes

Two leading impediments to chronic wound healing are polymicrobial infection and biofilm formation. Recent studies have characterized the bacterial fraction of these microbiomes and have begun to elucidate compositional correlations to healing outcomes. However, the factors that drive compositional shifts are still being uncovered. The virome may play an important role in shaping bacterial community structure and function. Previous work on the skin virome determined that it was dominated by bacteriophages, viruses that infect bacteria. To characterize the virome, we enrolled 20 chronic wound patients presenting at an outpatient wound care clinic in a microbiome survey, collecting swab samples from healthy skin and chronic wounds (diabetic, venous, arterial, or pressure) before and after a single, sharp debridement procedure. We investigated the virome using a virus-like particle enrichment procedure, shotgun metagenomic sequencing, and a k-mer-based, reference-dependent taxonomic classification method. Taxonomic composition, diversity, and associations with covariates are presented. We find that the wound virome is highly diverse, with many phages targeting known pathogens, and may influence bacterial community composition and functionality in ways that impact healing outcomes.

IMPORTANCE Chronic wounds are an increasing medical burden. These wounds are known to be rich in microbial content, including both bacteria and bacterial viruses (phages). The viruses may play an important role in shaping bacterial community structure and function. We analyzed the virome and bacterial composition of 20 patients with chronic wounds. The viruses found in wounds are highly diverse compared to normal skin, unlike the bacterial composition, where diversity is decreased. These data represent an initial look at this relatively understudied component of the chronic wound microbiome and may help inform future phage-based interventions.

Accurate prediction of metagenome-assembled genome completeness by MAGISTA, a random forest model built on alignment-free intra-bin statistics

BACKGROUND

Although the total number of microbial taxa on Earth is under debate, it is clear that only a small fraction of these has been cultivated and validly named. Evidently, the inability to culture most bacteria outside of very specific conditions severely limits their characterization and further studies. In the last decade, a major part of the solution to this problem has been the use of metagenome sequencing, whereby the DNA of an entire microbial community is sequenced, followed by the in silico reconstruction of genomes of its novel component species. The large discrepancy between the number of sequenced type strain genomes (around 12,000) and total microbial diversity (10⁶-10¹² species) directs these efforts to de novo assembly and binning. Unfortunately, these steps are error-prone and as such, the results have to be intensely scrutinized to avoid publishing incomplete and low-quality genomes.

RESULTS

We developed MAGISTA (metagenome-assembled genome intra-bin statistics assessment), a novel approach to assess metagenome-assembled genome quality that tackles some of the often-neglected drawbacks of current reference gene-based methods. MAGISTA is based on alignment-free distance distributions between contig fragments within metagenomic bins, rather than a set of reference genes. For proper training, a highly complex genomic DNA mock community was needed and constructed by pooling genomic DNA of 227 bacterial strains, specifically selected to obtain a wide variety representing the major phylogenetic lineages of cultivable bacteria.

CONCLUSIONS

MAGISTA achieved a 20% reduction in root-mean-square error in comparison to the marker gene approach when tested on publicly available mock metagenomes. Furthermore, our highly complex genomic DNA mock community is a very valuable tool for benchmarking (new) metagenome analysis methods.

Dual-layered hydrogels allow complete genome recovery with nucleic acid cytometry

Targeting specific cells for sequencing is important for applications in cancer, microbiology, and infectious disease. Nucleic acid cytometry is a powerful approach for accomplishing this because it allows specific cells to be isolated based on sequence biomarkers that are otherwise impossible to detect. However, existing methods require specialized microfluidic devices, limiting adoption. Here, we describe a modified workflow that uses particle-templated emulsification and flow cytometry to conduct the essential steps of cell detection and sorting normally accomplished by microfluidics. Our microfluidic-free workflow allows facile isolation and sequencing of cells, viruses, and nucleic acids and thus provides a powerful enrichment approach for targeted sequencing applications. This article is protected by copyright. All rights reserved.

Microbial colonization and resistome dynamics in food processing environments of a newly opened pork cutting industry during 1.5 years of activity

BACKGROUND

The microorganisms that inhabit food processing environments (FPE) can strongly influence the associated food quality and safety. In particular, the possibility that FPE may act as a reservoir of antibiotic-resistant microorganisms, and a hotspot for the transmission of antibiotic resistance genes (ARGs) is a concern in meat processing plants. Here, we monitor microbial succession and resistome dynamics relating to FPE through a detailed analysis of a newly opened pork cutting plant over 1.5 years of activity.

RESULTS

We identified a relatively restricted principal microbiota dominated by Pseudomonas during the first 2 months, while a higher taxonomic diversity, an increased representation of other taxa (e.g., Acinetobacter, Psychrobacter), and a certain degree of microbiome specialization on different surfaces was recorded later on. An increase in total abundance, alpha diversity, and β-dispersion of ARGs, which were predominantly assigned to Acinetobacter and associated with resistance to certain antimicrobials frequently used on pig farms of the region, was detected over time. Moreover, a sharp increase in the occurrence of extended-spectrum β-lactamase-producing Enterobacteriaceae and vancomycin-resistant Enterococcaceae was observed when cutting activities started. ARGs associated with resistance to β-lactams, tetracyclines, aminoglycosides, and sulphonamides frequently co-occurred, and mobile genetic elements (i.e., plasmids, integrons) and lateral gene transfer events were mainly detected at the later sampling times in drains.

CONCLUSIONS

The observations made suggest that pig carcasses were a source of resistant bacteria that then colonized FPE and that drains, together with some food-contact surfaces, such as equipment and table surfaces, represented a reservoir for the spread of ARGs in the meat processing facility. Video Abstract.

Host-associated microbe PCR (hamPCR) enables convenient measurement of both microbial load and community composition

The ratio of microbial population size relative to the amount of host tissue, or 'microbial load', is a fundamental metric of colonization and infection, but it cannot be directly deduced from microbial amplicon data such as 16S rRNA gene counts. Because existing methods to determine load, such as serial dilution plating, quantitative PCR, and whole metagenome sequencing add substantial cost and/or experimental burden, they are only rarely paired with amplicon sequencing. We introduce host-associated microbe PCR (hamPCR), a robust strategy to both quantify microbial load and describe interkingdom microbial community composition in a single amplicon library. We demonstrate its accuracy across multiple study systems, including nematodes and major crops, and further present a cost-saving technique to reduce host overrepresentation in the library prior to sequencing. Because hamPCR provides an accessible experimental solution to the well-known limitations and statistical challenges of compositional data, it has far-reaching potential in culture-independent microbiology.

Glacier ice archives nearly 15,000-year-old microbes and phages

BACKGROUND

Glacier ice archives information, including microbiology, that helps reveal paleoclimate histories and predict future climate change. Though glacier-ice microbes are studied using culture or amplicon approaches, more challenging metagenomic approaches, which provide access to functional, genome-resolved information and viruses, are under-utilized, partly due to low biomass and potential contamination.

RESULTS

We expand existing clean sampling procedures using controlled artificial ice-core experiments and adapted previously established low-biomass metagenomic approaches to study glacier-ice viruses. Controlled sampling experiments drastically reduced mock contaminants including bacteria, viruses, and free DNA to background levels. Amplicon sequencing from eight depths of two Tibetan Plateau ice cores revealed common glacier-ice lineages including Janthinobacterium, Polaromonas, Herminiimonas, Flavobacterium, Sphingomonas, and Methylobacterium as the dominant genera, while microbial communities were significantly different between two ice cores, associating with different climate conditions during deposition. Separately, ~355- and ~14,400-year-old ice were subject to viral enrichment and low-input quantitative sequencing, yielding genomic sequences for 33 vOTUs. These were virtually all unique to this study, representing 28 novel genera and not a single species shared with 225 environmentally diverse viromes. Further, 42.4% of the vOTUs were identifiable temperate, which is significantly higher than that in gut, soil, and marine viromes, and indicates that temperate phages are possibly favored in glacier-ice environments before being frozen. In silico host predictions linked 18 vOTUs to co-occurring abundant bacteria (Methylobacterium, Sphingomonas, and Janthinobacterium), indicating that these phages infected ice-abundant bacterial groups before being archived. Functional genome annotation revealed four virus-encoded auxiliary metabolic genes, particularly two motility genes suggest viruses potentially facilitate nutrient acquisition for their hosts. Finally, given their possible importance to methane cycling in ice, we focused on Methylobacterium viruses by contextualizing our ice-observed viruses against 123 viromes and prophages extracted from 131 Methylobacterium genomes, revealing that the archived viruses might originate from soil or plants.

CONCLUSIONS

Together, these efforts further microbial and viral sampling procedures for glacier ice and provide a first window into viral communities and functions in ancient glacier environments. Such methods and datasets can potentially enable researchers to contextualize new discoveries and begin to incorporate glacier-ice microbes and their viruses relative to past and present climate change in geographically diverse regions globally. Video Abstract.

Mosaic Ends Tagmentation (METa) Assembly for Highly Efficient Construction of Functional Metagenomic Libraries

Functional metagenomic libraries, physical bacterial libraries which allow the high-throughput capture and expression of microbiome genes, have been instrumental in the sequence-naive and cultivation-independent exploration of metagenomes. However, preparation of these libraries is often limited by their high DNA input requirement and their low cloning efficiency. Here, we describe a new method, mosaic ends tagmentation (METa) assembly, for highly efficient functional metagenomic library preparation. We applied tagmentation to metagenomic DNA from soil and gut microbiomes to prepare DNA inserts for high-throughput cloning into functional metagenomic libraries. The presence of mosaic end sequences in the resulting DNA fragments synergized with homology-based assembly cloning to result in a 300-fold increase in cloning efficiency compared to traditional blunt-cloning-based protocols. We show that compared to published libraries prepared by state-of-the-art protocols, METa assembly is on average ca. 20- to 200-fold more efficient and can prepare gigabase-sized libraries with as little as 200 ng of input DNA. We show the usefulness of METa assembly first by using a normative 5-μg mass of soil metagenomic DNA to prepare a 700-Gb library that allowed us to discover novel nourseothricin resistance genes and a potentially new mode of resistance to this antibiotic and second by using only 300 ng of goose fecal metagenomic DNA to prepare a 27-Gb library that captured numerous tetracycline and colistin resistance genes. METa assembly provides a streamlined, flexible, and efficient method for preparing functional metagenomic libraries, enabling new avenues of genetic and biochemical research into low-biomass or scarce microbiomes. IMPORTANCE Medically and industrially important genes can be recovered from microbial communities by high-throughput sequencing, but precise annotation is often limited to characterized genes and their relatives. Cloning a metagenome en masse into an expression host to produce a functional metagenomic library, directly connecting genes to functions, is a sequence-naive and cultivation-independent method to discover novel genes. The process of preparing these libraries is DNA greedy and inefficient, however. Here, we describe a library preparation method that is an order of magnitude more efficient and less DNA greedy. This method is consistently efficient across libraries prepared from cultures, a soil microbiome, and a goose fecal microbiome and allowed us to discover new antibiotic resistance genes and mechanisms. This library preparation method will potentially allow the functional metagenomic exploration of microbiomes that were previously off limits due to their rarity or low microbial biomass, such as biomedical swabs or exotic samples.

Metaviromes Reveal the Dynamics of Pseudomonas Host-Specific Phages Cultured and Uncultured by Plaque Assay

Isolating single phages using plaque assays is a laborious and time-consuming process. Whether single isolated phages are the most lyse-effective, the most abundant in viromes, or those with the highest ability to make plaques in solid media is not well known. With the increasing accessibility of high-throughput sequencing, metaviromics is often used to describe viruses in environmental samples. By extracting and sequencing metaviromes from organic waste with and without exposure to a host-of-interest, we show a host-related phage community's shift, as well as identify the most enriched phages. Moreover, we isolated plaque-forming single phages using the same virome-host matrix to observe how enrichments in liquid media correspond to the metaviromic data. In this study, we observed a significant shift (p = 0.015) of the 47 identified putative Pseudomonas phages with a minimum twofold change above zero in read abundance when adding a Pseudomonas syringae DC3000 host. Surprisingly, it appears that only two out of five plaque-forming phages from the same organic waste sample, targeting the Pseudomonas strain, were highly abundant in the metavirome, while the other three were almost absent despite host exposure. Lastly, our sequencing results highlight how long reads from Oxford Nanopore elevates the assembly quality of metaviromes, compared to short reads alone.

CheckV assesses the quality and completeness of metagenome-assembled viral genomes

Millions of new viral sequences have been identified from metagenomes, but the quality and completeness of these sequences vary considerably. Here we present CheckV, an automated pipeline for identifying closed viral genomes, estimating the completeness of genome fragments and removing flanking host regions from integrated proviruses. CheckV estimates completeness by comparing sequences with a large database of complete viral genomes, including 76,262 identified from a systematic search of publicly available metagenomes, metatranscriptomes and metaviromes. After validation on mock datasets and comparison to existing methods, we applied CheckV to large and diverse collections of metagenome-assembled viral sequences, including IMG/VR and the Global Ocean Virome. This revealed 44,652 high-quality viral genomes (that is, >90% complete), although the vast majority of sequences were small fragments, which highlights the challenge of assembling viral genomes from short-read metagenomes. Additionally, we found that removal of host contamination substantially improved the accurate identification of auxiliary metabolic genes and interpretation of viral-encoded functions.

STROBE-metagenomics: a STROBE extension statement to guide the reporting of metagenomics studies

The term metagenomics refers to the use of sequencing methods to simultaneously identify genomic material from all organisms present in a sample, with the advantage of greater taxonomic resolution than culture or other methods. Applications include pathogen detection and discovery, species characterisation, antimicrobial resistance detection, virulence profiling, and study of the microbiome and microecological factors affecting health. However, metagenomics involves complex and multistep processes and there are important technical and methodological challenges that require careful consideration to support valid inference. We co-ordinated a multidisciplinary, international expert group to establish reporting guidelines that address specimen processing, nucleic acid extraction, sequencing platforms, bioinformatics considerations, quality assurance, limits of detection, power and sample size, confirmatory testing, causality criteria, cost, and ethical issues. The guidance recognises that metagenomics research requires pragmatism and caution in interpretation, and that this field is rapidly evolving.

Detection of mixed-strain infections by FACS and ultra-low input genome sequencing

The epidemiological tracking of a bacterial outbreak may be jeopardized by the presence of multiple pathogenic strains in one patient. Nevertheless, this fact is not considered in most of the epidemiological studies and only one colony per patient is sequenced. On the other hand, the routine whole genome sequencing of many isolates from each patient would be costly and unnecessary, because the number of strains in a patient is never known a priori. In addition, the result would be biased by microbial culture conditions. Herein we propose an approach for detecting mixed-strain infection, providing C. difficile infection as an example. The cells of the target pathogenic species are collected from the bacterial suspension by the fluorescence activated cell sorting (FACS) and a shallow genome sequencing is performed. A modified sequencing library preparation protocol for low-input DNA samples can be used for low prevalence gut pathogens (< 0.1% of the total microbiome). This FACS-seq approach reduces diagnostics time (no culture is needed) and may promote discoveries of novel strains. Methodological details, possible issues and future directions for the sequencing of these natural pan-genomes are herein discussed.

Morphology of obligate ectosymbionts reveals Paralaxus gen. nov.: A new circumtropical genus of marine stilbonematine nematodes

Stilbonematinae are a subfamily of conspicuous marine nematodes, distinguished by a coat of sulphur-oxidizing bacterial ectosymbionts on their cuticle. As most nematodes, the worm hosts have a relatively simple anatomy and few taxonomically informative characters, and this has resulted in numerous taxonomic reassignments and synonymizations. Recent studies using a combination of morphological and molecular traits have helped to improve the taxonomy of Stilbonematinae but also raised questions on the validity of several genera. Here, we describe a new circumtropically distributed genus Paralaxus (Stilbonematinae) with three species: Paralaxus cocos sp. nov., P. bermudensis sp. nov. and P. columbae sp. nov. We used single worm metagenomes to generate host 18S rRNA and cytochrome c oxidase I (COI) as well as symbiont 16S rRNA gene sequences. Intriguingly, COI alignments and primer matching analyses suggest that the COI is not suitable for PCR-based barcoding approaches in Stilbonematinae as the genera have a highly diverse base composition and no conserved primer sites. The phylogenetic analyses of all three gene sets, however, confirm the morphological assignments and support the erection of the new genus Paralaxus as well as corroborate the status of the other stilbonematine genera. Paralaxus most closely resembles the stilbonematine genus Laxus in overlapping sets of diagnostic features but can be distinguished from Laxus by the morphology of the genus-specific symbiont coat. Our re-analyses of key parameters of the symbiont coat morphology as character for all Stilbonematinae genera show that with amended descriptions, including the coat, highly reliable genus assignments can be obtained.

A pipeline for targeted metagenomics of environmental bacteria

BACKGROUND

Metagenomics and single cell genomics provide a window into the genetic repertoire of yet uncultivated microorganisms, but both methods are usually taxonomically untargeted. The combination of fluorescence in situ hybridization (FISH) and fluorescence activated cell sorting (FACS) has the potential to enrich taxonomically well-defined clades for genomic analyses.

METHODS

Cells hybridized with a taxon-specific FISH probe are enriched based on their fluorescence signal via flow cytometric cell sorting. A recently developed FISH procedure, the hybridization chain reaction (HCR)-FISH, provides the high signal intensities required for flow cytometric sorting while maintaining the integrity of the cellular DNA for subsequent genome sequencing. Sorted cells are subjected to shotgun sequencing, resulting in targeted metagenomes of low diversity.

RESULTS

Pure cultures of different taxonomic groups were used to (1) adapt and optimize the HCR-FISH protocol and (2) assess the effects of various cell fixation methods on both the signal intensity for cell sorting and the quality of subsequent genome amplification and sequencing. Best results were obtained for ethanol-fixed cells in terms of both HCR-FISH signal intensity and genome assembly quality. Our newly developed pipeline was successfully applied to a marine plankton sample from the North Sea yielding good quality metagenome assembled genomes from a yet uncultivated flavobacterial clade.

CONCLUSIONS

With the developed pipeline, targeted metagenomes at various taxonomic levels can be efficiently retrieved from environmental samples. The resulting metagenome assembled genomes allow for the description of yet uncharacterized microbial clades. Video abstract.

Optimizing sequencing protocols for leaderboard metagenomics by combining long and short reads

As metagenomic studies move to increasing numbers of samples, communities like the human gut may benefit more from the assembly of abundant microbes in many samples, rather than the exhaustive assembly of fewer samples. We term this approach leaderboard metagenome sequencing. To explore protocol optimization for leaderboard metagenomics in real samples, we introduce a benchmark of library prep and sequencing using internal references generated by synthetic long-read technology, allowing us to evaluate high-throughput library preparation methods against gold-standard reference genomes derived from the samples themselves. We introduce a low-cost protocol for high-throughput library preparation and sequencing.

GraftM: a tool for scalable, phylogenetically informed classification of genes within metagenomes

Large-scale metagenomic datasets enable the recovery of hundreds of population genomes from environmental samples. However, these genomes do not typically represent the full diversity of complex microbial communities. Gene-centric approaches can be used to gain a comprehensive view of diversity by examining each read independently, but traditional pairwise comparison approaches typically over-classify taxonomy and scale poorly with increasing metagenome and database sizes. Here we introduce GraftM, a tool that uses gene specific packages to rapidly identify gene families in metagenomic data using hidden Markov models (HMMs) or DIAMOND databases, and classifies these sequences using placement into pre-constructed gene trees. The speed and accuracy of GraftM was benchmarked with in silico and in vitro mock communities using taxonomic markers, and was found to have higher accuracy at the family level with a processing time 2.0–3.7× faster than currently available software. Exploration of a wetland metagenome using 16S rRNA- and methyl-coenzyme M reductase (McrA)-specific gpkgs revealed taxonomic and functional shifts across a depth gradient. Analysis of the NCBI nr database using the McrA gpkg allowed the detection of novel sequences belonging to phylum-level lineages. A growing collection of gpkgs is available online (https://github.com/geronimp/graftM_gpkgs), where curated packages can be uploaded and exchanged.

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.

Single-Cell Genomics of Microbial Dark Matter

Single-cell genomics allows bypassing the culturing step and to directly access environmental microbes one cell at a time. The method has been successfully applied to explore archaeal and bacterial candidate phyla, referred to as microbial dark matter. Here I summarize the single-cell genomics workflow, including sample preparation and preservation, high-throughput fluorescence-activated cell sorting, cell lysis and amplification of environmental samples. Furthermore I describe phylogenetic screening based on 16S rRNA genes and suggest a suitable library preparation and sequencing approach.

Evaluation of Sequencing Library Preparation Protocols for Viral Metagenomic Analysis from Pristine Aquifer Groundwaters

Viral ecology of terrestrial habitats is yet-to be extensively explored, in particular the terrestrial subsurface. One problem in obtaining viral sequences from groundwater aquifer samples is the relatively low amount of virus particles. As a result, the amount of extracted DNA may not be sufficient for direct sequencing of such samples. Here we compared three DNA amplification methods to enrich viral DNA from three pristine limestone aquifer assemblages of the Hainich Critical Zone Exploratory to evaluate potential bias created by the different amplification methods as determined by viral metagenomics. Linker amplification shotgun libraries resulted in lowest redundancy among the sequencing reads and showed the highest diversity, while multiple displacement amplification produced the highest number of contigs with the longest average contig size, suggesting a combination of these two methods is suitable for the successful enrichment of viral DNA from pristine groundwater samples. In total, we identified 27,173, 5,886 and 32,613 viral contigs from the three samples from which 11.92 to 18.65% could be assigned to taxonomy using blast. Among these, members of the Caudovirales order were the most abundant group (52.20 to 69.12%) dominated by Myoviridae and Siphoviridae. Those, and the high number of unknown viral sequences, substantially expand the known virosphere.

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.

High-Throughput Miniaturized 16S rRNA Amplicon Library Preparation Reduces Costs while Preserving Microbiome Integrity

Next-generation sequencing technologies have enabled many advances across biology, with microbial ecology benefiting primarily through expanded sample sizes. Although the cost of running sequencing instruments has decreased substantially over time, the price of library preparation methods has largely remained unchanged. In this study, we developed a low-cost miniaturized (5-µl volume) high-throughput (384-sample) amplicon library preparation method with the Echo 550 acoustic liquid handler. Our method reduces costs of library preparation to $1.42 per sample, a 58% reduction compared to existing automated methods and a 21-fold reduction from commercial kits, without compromising sequencing success or distorting the microbial community composition analysis. We further validated the optimized method by sampling five body sites from 46 Pacific chub mackerel fish caught across 16 sampling events over seven months from the Scripps Institution of Oceanography pier in La Jolla, CA. Fish microbiome samples were processed with the miniaturized 5-µl reaction volume with 0.2 µl of genomic DNA (gDNA) and the standard 25-µl reaction volume with 1 µl of gDNA. Between the two methods, alpha diversity was highly correlated (R 2 > 0.95), while distances of technical replicates were much lower than within-body-site variation (P < 0.0001), further validating the method. The cost savings of implementing the miniaturized library preparation (going from triplicate 25-µl reactions to triplicate 5-µl reactions) are large enough to cover a MiSeq sequencing run for 768 samples while preserving accurate microbiome measurements. IMPORTANCE Reduced costs of sequencing have tremendously impacted the field of microbial ecology, allowing scientists to design more studies with larger sample sizes that often exceed 10,000 samples. Library preparation costs have not kept pace with sequencing prices, although automated liquid handling robots provide a unique opportunity to bridge this gap while also decreasing human error. Here, we take advantage of an acoustic liquid handling robot to develop a high-throughput miniaturized library preparation method of a highly cited and broadly used 16S rRNA gene amplicon reaction. We evaluate the potential negative effects of reducing the PCR volume along with varying the amount of gDNA going into the reaction. Our optimized method reduces sample-processing costs while continuing to generate a high-quality microbiome readout that is indistinguishable from the original method.

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.

Synthetic microbe communities provide internal reference standards for metagenome sequencing and analysis

The complexity of microbial communities, combined with technical biases in next-generation sequencing, pose a challenge to metagenomic analysis. Here, we develop a set of internal DNA standards, termed “sequins” (sequencing spike-ins), that together constitute a synthetic community of artificial microbial genomes. Sequins are added to environmental DNA samples prior to library preparation, and undergo concurrent sequencing with the accompanying sample. We validate the performance of sequins by comparison to mock microbial communities, and demonstrate their use in the analysis of real metagenome samples. We show how sequins can be used to measure fold change differences in the size and structure of accompanying microbial communities, and perform quantitative normalization between samples. We further illustrate how sequins can be used to benchmark and optimize new methods, including nanopore long-read sequencing technology. We provide metagenome sequins, along with associated data sets, protocols, and an accompanying software toolkit, as reference standards to aid in metagenomic studies.

Complex microbial communities pose a challenge to metagenomic analysis. Here the authors develop ‘sequins’, internal DNA standards that represent a synthetic community of artificial genomes.

Enrichment of Clinically Relevant Organisms in Spontaneous Preterm-Delivered Placentas and Reagent Contamination across All Clinical Groups in a Large Pregnancy Cohort in the United Kingdom

Preterm birth is associated with both psychological and physical disabilities and is the leading cause of infant morbidity and mortality worldwide. Infection is known to be an important cause of spontaneous preterm birth, and recent research has implicated variation in the “placental microbiome” in the risk of preterm birth. Consistent with data from previous studies, the abundances of certain clinically relevant species differed between spontaneous preterm- and nonspontaneous preterm- or term-delivered placentas. These results support the view that a proportion of spontaneous preterm births have an intrauterine-infection component. However, an additional observation from this study was that a substantial proportion of sequenced reads were contaminating reads rather than DNA from endogenous, clinically relevant species. This observation warrants caution in the interpretation of sequencing outputs from low-biomass samples such as the placenta.

In this study, differences in the placental microbiota from term and preterm deliveries in a large pregnancy cohort in the United Kingdom were studied by using 16S-targeted amplicon sequencing. The impacts of contamination from DNA extraction, PCR reagents, and the delivery itself were also examined. A total of 400 placental samples from 256 singleton pregnancies were analyzed, and differences between spontaneous preterm-, nonspontaneous preterm-, and term-delivered placentas were investigated. DNA from recently delivered placentas was extracted, and screening for bacterial DNA was carried out by using targeted sequencing of the 16S rRNA gene on the Illumina MiSeq platform. Sequenced reads were analyzed for the presence of contaminating operational taxonomic units (OTUs) identified via sequencing of negative extraction and PCR-blank samples. Differential abundances and between-sample (beta) diversity metrics were then compared. A large proportion of the reads sequenced from the extracted placental samples mapped to OTUs that were also found for negative extractions. Striking differences in the compositions of samples were also observed, according to whether the placenta was delivered abdominally or vaginally, providing strong circumstantial evidence for delivery contamination as an important contributor to observed microbial profiles. When OTU- and genus-level abundances were compared between the groups of interest, a number of organisms were enriched in the spontaneous preterm-delivery cohort, including organisms that have been associated previously with adverse pregnancy outcomes, specifically Mycoplasma spp. and Ureaplasma spp. However, analyses of the overall community structure did not reveal convincing evidence for the existence of a reproducible “preterm placental microbiome.”

IMPORTANCE Preterm birth is associated with both psychological and physical disabilities and is the leading cause of infant morbidity and mortality worldwide. Infection is known to be an important cause of spontaneous preterm birth, and recent research has implicated variation in the “placental microbiome” in the risk of preterm birth. Consistent with data from previous studies, the abundances of certain clinically relevant species differed between spontaneous preterm- and nonspontaneous preterm- or term-delivered placentas. These results support the view that a proportion of spontaneous preterm births have an intrauterine-infection component. However, an additional observation from this study was that a substantial proportion of sequenced reads were contaminating reads rather than DNA from endogenous, clinically relevant species. This observation warrants caution in the interpretation of sequencing outputs from low-biomass samples such as the placenta.

The Life Aquatic at the Microscale

There are more than one million microbial cells in every drop of seawater, and their collective metabolisms not only recycle nutrients that can then be used by larger organisms but also catalyze key chemical transformations that maintain Earth’s habitability. Understanding how these microbes interact with each other and with multicellular hosts is critical to reliably quantify any functional aspect of their metabolisms and to predicting their outcomes on larger scales.

There are more than one million microbial cells in every drop of seawater, and their collective metabolisms not only recycle nutrients that can then be used by larger organisms but also catalyze key chemical transformations that maintain Earth’s habitability. Understanding how these microbes interact with each other and with multicellular hosts is critical to reliably quantify any functional aspect of their metabolisms and to predict their outcomes on larger scales. Following a large body of literature pioneered by Farooq Azam and colleagues more than 30 years ago, I emphasize the importance of studying microbial interactions at the appropriate scale if we want to fully decipher the roles that they play in oceanic ecosystems.

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.

Library Construction from Subnanogram DNA for Pelagic Sea Water and Deep-Sea Sediments

Shotgun metagenomics is a low biased technology for assessing environmental microbial diversity and function. However, the requirement for a sufficient amount of DNA and the contamination of inhibitors in environmental DNA leads to difficulties in constructing a shotgun metagenomic library. We herein examined metagenomic library construction from subnanogram amounts of input environmental DNA from subarctic surface water and deep-sea sediments using two library construction kits: the KAPA Hyper Prep Kit and Nextera XT DNA Library Preparation Kit, with several modifications. The influence of chemical contaminants associated with these environmental DNA samples on library construction was also investigated. Overall, shotgun metagenomic libraries were constructed from 1 pg to 1 ng of input DNA using both kits without harsh library microbial contamination. However, the libraries constructed from 1 pg of input DNA exhibited larger biases in GC contents, k-mers, or small subunit (SSU) rRNA gene compositions than those constructed from 10 pg to 1 ng DNA. The lower limit of input DNA for low biased library construction in this study was 10 pg. Moreover, we revealed that technology-dependent biases (physical fragmentation and linker ligation vs. tagmentation) were larger than those due to the amount of input DNA.

Microbial ecology of mountain glacier ecosystems: biodiversity, ecological connections and implications of a warming climate

Glacier ecosystems are teeming with life on, beneath, and to a lesser degree, within their icy masses. This conclusion largely stems from polar research, with less attention paid to mountain glaciers that overlap environmentally and ecologically with their polar counterparts in some ways, but diverge in others. One difference lies in the susceptibility of mountain glaciers to the near-term threat of climate change, as they tend to be much smaller in both area and volume. Moreover, mountain glaciers are typically steeper, more dependent upon basal sliding for movement, and experience higher seasonal precipitation. Here, we provide a modern synthesis of the microbial ecology of mountain glacier ecosystems, and particularly those at low- to mid-latitudes. We focus on five ecological zones: the supraglacial surface, englacial interior, subglacial bedrock-ice interface, proglacial streams and glacier forefields. For each, we discuss the role of microbiota in biogeochemical cycling and outline ecological and hydrological connections among zones, underscoring the interconnected nature of these ecosystems. Collectively, we highlight the need to: better document the biodiversity and functional roles of mountain glacier microbiota; describe the ecological implications of rapid glacial retreat under climate change and resolve the relative contributions of ecological zones to broader ecosystem function.

Zooming in on the phycosphere: the ecological interface for phytoplankton-bacteria relationships

By controlling nutrient cycling and biomass production at the base of the food web, interactions between phytoplankton and bacteria represent a fundamental ecological relationship in aquatic environments. Although typically studied over large spatiotemporal scales, emerging evidence indicates that this relationship is often governed by microscale interactions played out within the region immediately surrounding individual phytoplankton cells. This microenvironment, known as the phycosphere, is the planktonic analogue of the rhizosphere in plants. The exchange of metabolites and infochemicals at this interface governs phytoplankton-bacteria relationships, which span mutualism, commensalism, antagonism, parasitism and competition. The importance of the phycosphere has been postulated for four decades, yet only recently have new technological and conceptual frameworks made it possible to start teasing apart the complex nature of this unique microbial habitat. It has subsequently become apparent that the chemical exchanges and ecological interactions between phytoplankton and bacteria are far more sophisticated than previously thought and often require close proximity of the two partners, which is facilitated by bacterial colonization of the phycosphere. It is also becoming increasingly clear that while interactions taking place within the phycosphere occur at the scale of individual microorganisms, they exert an ecosystem-scale influence on fundamental processes including nutrient provision and regeneration, primary production, toxin biosynthesis and biogeochemical cycling. Here we review the fundamental physical, chemical and ecological features of the phycosphere, with the goal of delivering a fresh perspective on the nature and importance of phytoplankton-bacteria interactions in aquatic ecosystems.

MetaSort untangles metagenome assembly by reducing microbial community complexity

Most current approaches to analyse metagenomic data rely on reference genomes. Novel microbial communities extend far beyond the coverage of reference databases and de novo metagenome assembly from complex microbial communities remains a great challenge. Here we present a novel experimental and bioinformatic framework, metaSort, for effective construction of bacterial genomes from metagenomic samples. MetaSort provides a sorted mini-metagenome approach based on flow cytometry and single-cell sequencing methodologies, and employs new computational algorithms to efficiently recover high-quality genomes from the sorted mini-metagenome by the complementary of the original metagenome. Through extensive evaluations, we demonstrated that metaSort has an excellent and unbiased performance on genome recovery and assembly. Furthermore, we applied metaSort to an unexplored microflora colonized on the surface of marine kelp and successfully recovered 75 high-quality genomes at one time. This approach will greatly improve access to microbial genomes from complex or novel communities.

Impact of sample processing on human airways microbial metagenomes

Whole metagenome shotgun sequencing provides information about the gene content and the composition of microbial communities provided that the processing of the samples does not introduce a methodology-driven bias. We tested the impact of DNA isolation and storage period on the metagenome profile. Deep throat swabs were collected from healthy adults and an infected infant. DNA was isolated by sonification or enzymatic lysis either immediately or after 24h storage in agar gel Amies transport medium at room temperature. Disruption of cells and subsequent fragmentation of DNA by sonification was as suitable as the common enzymatic lysis to generate high-quality metagenomes particularly for low total DNA input of less than ten nanograms. Conversely, storage of samples for 24h produced severely distorted metagenomes. The majority of species became less abundant or even extinct, whereas a few Streptococcus, Neisseria and Haemophilus spp. proliferated so that the total number of bacterial reads increased at the expense of human reads. We recommend that samples for metagenome analysis should be immediately processed or frozen at -80°C.