Nonpareil: a redundancy-based approach to assess the level of coverage in metagenomic datasets

Recovering complete and draft population genomes from metagenome datasets

Assembly of metagenomic sequence data into microbial genomes is of fundamental value to improving our understanding of microbial ecology and metabolism by elucidating the functional potential of hard-to-culture microorganisms. Here, we provide a synthesis of available methods to bin metagenomic contigs into species-level groups and highlight how genetic diversity, sequencing depth, and coverage influence binning success. Despite the computational cost on application to deeply sequenced complex metagenomes (e.g., soil), covarying patterns of contig coverage across multiple datasets significantly improves the binning process. We also discuss and compare current genome validation methods and reveal how these methods tackle the problem of chimeric genome bins i.e., sequences from multiple species. Finally, we explore how population genome assembly can be used to uncover biogeographic trends and to characterize the effect of in situ functional constraints on the genome-wide evolution.

Metagenomic Comparison of Antibiotic Resistance Genes Associated with Liquid and Dewatered Biosolids

Municipal biosolids (MBs) that are land-applied in North America are known to possess an active microbial population that can include human pathogens. Activated sludge is a hotspot for the accumulation of antibiotics and has been shown to be a selective environment for microorganisms that contain antibiotic resistance genes (ARGs); however, the prevalence of ARGs in MBs is not well characterized. In this study, we enriched the plasmid metagenome from raw sewage sludge and two CP2 MBs, a mesophilic anaerobic digestate and a dewatered digestate, to evaluate the presence of ARGs in mobile genetic elements. The CP2-class biosolids are similar to Class B biosolids in the United States. The CP2 biosolids must meet a microbiological cut off of 2 × 10 colony-forming units (CFU) per dry gram or 100 mL of biosolids. The enriched plasmid DNA was sequenced (Illumina MiSeq). Sequence matching against databases, including the Comprehensive Antibiotic Resistance Database (CARD), MG-RAST, and INTEGRALL, identified potential genes of interest related to ARGs and their ability to transfer. The presence and abundance of different ARGs varied between treatments with heterogeneity observed among the same sample types. The MBs plasmid-enriched metagenomes contained ARGs associated with resistance to a variety of antibiotics, including β-lactams, rifampicin, quinolone, and tetracycline as well as the detection of extended spectrum β-lactamase genes. Cultured bacteria from CP2 MBs possessed antibiotic resistances consistent with the MBs metagenome data including multiantibiotic-resistant isolates. The results from this study provide a better understanding of the ARG and MGE profile of the plasmid-enriched metagenome of CP2 MBs.

Draft genome sequence of Halomonas meridiana R1t3 isolated from the surface microbiota of the Caribbean Elkhorn coral Acropora palmata

Members of the gammaproteobacterial genus Halomonas are common in marine environments. Halomonas and other members of the Oceanospirillales have recently been identified as prominent members of the surface microbiota of reef-building corals. Halomonas meridiana strain R1t3 was isolated from the surface mucus layer of the scleractinian coral Acropora palmata in 2005 from the Florida Keys. This strain was chosen for genome sequencing to provide insight into the role of commensal heterotrophic bacteria in the coral holobiont. The draft genome consists of 290 scaffolds, totaling 3.5 Mbp in length and contains 3397 protein-coding genes.

Strong spurious transcription likely contributes to DNA insert bias in typical metagenomic clone libraries

BACKGROUND

Clone libraries provide researchers with a powerful resource to study nucleic acid from diverse sources. Metagenomic clone libraries in particular have aided in studies of microbial biodiversity and function, and allowed the mining of novel enzymes. Libraries are often constructed by cloning large inserts into cosmid or fosmid vectors. Recently, there have been reports of GC bias in fosmid metagenomic libraries, and it was speculated to be a result of fragmentation and loss of AT-rich sequences during cloning. However, evidence in the literature suggests that transcriptional activity or gene product toxicity may play a role.

RESULTS

To explore possible mechanisms responsible for sequence bias in clone libraries, we constructed a cosmid library from a human microbiome sample and sequenced DNA from different steps during library construction: crude extract DNA, size-selected DNA, and cosmid library DNA. We confirmed a GC bias in the final cosmid library, and we provide evidence that the bias is not due to fragmentation and loss of AT-rich sequences but is likely occurring after DNA is introduced into Escherichia coli. To investigate the influence of strong constitutive transcription, we searched the sequence data for promoters and found that rpoD/σ(70) promoter sequences were underrepresented in the cosmid library. Furthermore, when we examined the genomes of taxa that were differentially abundant in the cosmid library relative to the original sample, we found the bias to be more correlated with the number of rpoD/σ(70) consensus sequences in the genome than with simple GC content.

CONCLUSIONS

The GC bias of metagenomic libraries does not appear to be due to DNA fragmentation. Rather, analysis of promoter sequences provides support for the hypothesis that strong constitutive transcription from sequences recognized as rpoD/σ(70) consensus-like in E. coli may lead to instability, causing loss of the plasmid or loss of the insert DNA that gives rise to the transcription. Despite widespread use of E. coli to propagate foreign DNA in metagenomic libraries, the effects of in vivo transcriptional activity on clone stability are not well understood. Further work is required to tease apart the effects of transcription from those of gene product toxicity.

Microbial community successional patterns in beach sands impacted by the Deepwater Horizon oil spill

Although petroleum hydrocarbons discharged from the Deepwater Horizon (DWH) blowout were shown to have a pronounced impact on indigenous microbial communities in the Gulf of Mexico, effects on nearshore or coastal ecosystems remain understudied. This study investigated the successional patterns of functional and taxonomic diversity for over 1 year after the DWH oil was deposited on Pensacola Beach sands (FL, USA), using metagenomic and 16S rRNA gene amplicon techniques. Gamma- and Alphaproteobacteria were enriched in oiled sediments, in corroboration of previous studies. In contrast to previous studies, we observed an increase in the functional diversity of the community in response to oil contamination and a functional transition from generalist populations within 4 months after oil came ashore to specialists a year later, when oil was undetectable. At the latter time point, a typical beach community had reestablished that showed little to no evidence of oil hydrocarbon degradation potential, was enriched in archaeal taxa known to be sensitive to xenobiotics, but differed significantly from the community before the oil spill. Further, a clear succession pattern was observed, where early responders to oil contamination, likely degrading aliphatic hydrocarbons, were replaced after 3 months by populations capable of aromatic hydrocarbon decomposition. Collectively, our results advance the understanding of how natural benthic microbial communities respond to crude oil perturbation, supporting the specialization-disturbance hypothesis; that is, the expectation that disturbance favors generalists, while providing (microbial) indicator species and genes for the chemical evolution of oil hydrocarbons during degradation and weathering.

Isolation and metagenomic characterization of bacteria associated with calcium carbonate and struvite precipitation in a pure moving bed biofilm reactor-membrane bioreactor

A bench-scale pure moving bed bioreactor-membrane bioreactor (MBBR-MBR) used for the treatment of urban wastewater was analyzed for the identification of bacterial strains with the potential capacity for calcium carbonate and struvite biomineral formation. Isolation of mineral-forming strains on calcium carbonate and struvite media revealed six major colonies with a carbonate or struvite precipitation capacity in the biofouling on the membrane surface and showed that heterotrophic bacteria with the ability to precipitate calcium carbonate and struvite constituted ~7.5% of the total platable bacteria. These belonged to the genera Lysinibacillus, Trichococcus, Comamomas and Bacillus. Pyrosequencing analysis of the microbial communities in the suspended cells and membrane biofouling showed a high degree of similarity in all the samples collected with respect to bacterial assemblage. The study of operational taxonomic units (OTUs) identified through pyrosequencing suggested that ~21% of the total bacterial community identified in the biofouling could potentially form calcium carbonate or struvite crystals in the pure MBBR-MBR system used for the treatment of urban wastewater.

Strengths and limitations of 16S rRNA gene amplicon sequencing in revealing temporal microbial community dynamics

This study explored the short-term planktonic microbial community structure and resilience in Lake Lanier (GA, USA) while simultaneously evaluating the technical aspects of identifying taxa via 16S rRNA gene amplicon and metagenomic sequence data. 16S rRNA gene amplicons generated from four temporally discrete samples were sequenced with 454 GS-FLX-Ti yielding ∼40,000 rRNA gene sequences from each sample and representing ∼300 observed OTUs. Replicates obtained from the same biological sample clustered together but several biases were observed, linked to either the PCR or sequencing-preparation steps. In comparisons with companion whole-community shotgun metagenome datasets, the estimated number of OTUs at each timepoint was concordant, but 1.5 times and ∼10 times as many phyla and genera, respectively, were identified in the metagenomes. Our analyses showed that the 16S rRNA gene captures broad shifts in community diversity over time, but with limited resolution and lower sensitivity compared to metagenomic data. We also identified OTUs that showed marked shifts in abundance over four close timepoints separated by perturbations and tracked these taxa in the metagenome vs. 16S rRNA amplicon data. A strong summer storm had less of an effect on community composition than did seasonal mixing, which revealed a distinct succession of organisms. This study provides insights into freshwater microbial communities and advances the approaches for assessing community diversity and dynamics in situ.

Microbial metabolic potential for carbon degradation and nutrient (nitrogen and phosphorus) acquisition in an ombrotrophic peatland

This study integrated metagenomic and nuclear magnetic resonance (NMR) spectroscopic approaches to investigate microbial metabolic potential for organic matter decomposition and nitrogen (N) and phosphorus (P) acquisition in soils of an ombrotrophic peatland in the Marcell Experimental Forest (MEF), Minnesota, USA. This analysis revealed vertical stratification in key enzymatic pathways and taxa containing these pathways. Metagenomic analyses revealed that genes encoding laccases and dioxygenases, involved in aromatic compound degradation, declined in relative abundance with depth, while the relative abundance of genes encoding metabolism of amino sugars and all four saccharide groups increased with depth in parallel with a 50% reduction in carbohydrate content. Most Cu-oxidases were closely related to genes from Proteobacteria and Acidobacteria, and type 4 laccase-like Cu-oxidase genes were >8 times more abundant than type 3 genes, suggesting an important and overlooked role for type 4 Cu-oxidase in phenolic compound degradation. Genes associated with sulfate reduction and methanogenesis were the most abundant anaerobic respiration genes in these systems, with low levels of detection observed for genes of denitrification and Fe(III) reduction. Fermentation genes increased in relative abundance with depth and were largely affiliated with Syntrophobacter. Methylocystaceae-like small-subunit (SSU) rRNA genes, pmoA, and mmoX genes were more abundant among methanotrophs. Genes encoding N2 fixation, P uptake, and P regulons were significantly enriched in the surface peat and in comparison to other ecosystems, indicating N and P limitation. Persistence of inorganic orthophosphate throughout the peat profile in this P-limiting environment indicates that P may be bound to recalcitrant organic compounds, thus limiting P bioavailability in the subsurface. Comparative metagenomic analysis revealed a high metabolic potential for P transport and starvation, N2 fixation, and oligosaccharide degradation at MEF relative to other wetland and soil environments, consistent with the nutrient-poor and carbohydrate-rich conditions found in this Sphagnum-dominated boreal peatland.

Soil microbial community responses to a decade of warming as revealed by comparative metagenomics

Soil microbial communities are extremely complex, being composed of thousands of low-abundance species (<0.1% of total). How such complex communities respond to natural or human-induced fluctuations, including major perturbations such as global climate change, remains poorly understood, severely limiting our predictive ability for soil ecosystem functioning and resilience. In this study, we compared 12 whole-community shotgun metagenomic data sets from a grassland soil in the Midwestern United States, half representing soil that had undergone infrared warming by 2°C for 10 years, which simulated the effects of climate change, and the other half representing the adjacent soil that received no warming and thus, served as controls. Our analyses revealed that the heated communities showed significant shifts in composition and predicted metabolism, and these shifts were community wide as opposed to being attributable to a few taxa. Key metabolic pathways related to carbon turnover, such as cellulose degradation (∼13%) and CO2 production (∼10%), and to nitrogen cycling, including denitrification (∼12%), were enriched under warming, which was consistent with independent physicochemical measurements. These community shifts were interlinked, in part, with higher primary productivity of the aboveground plant communities stimulated by warming, revealing that most of the additional, plant-derived soil carbon was likely respired by microbial activity. Warming also enriched for a higher abundance of sporulation genes and genomes with higher G+C content. Collectively, our results indicate that microbial communities of temperate grassland soils play important roles in mediating feedback responses to climate change and advance the understanding of the molecular mechanisms of community adaptation to environmental perturbations.