1
|
Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis. NATURE PLANTS 2024; 10:673-688. [PMID: 38589485 PMCID: PMC11035146 DOI: 10.1038/s41477-024-01666-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Accepted: 03/06/2024] [Indexed: 04/10/2024]
Abstract
The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.
Collapse
|
2
|
Quality MAGnified. Nat Rev Microbiol 2023; 21:771. [PMID: 37789075 DOI: 10.1038/s41579-023-00981-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
|
3
|
Multicellular magnetotactic bacterial consortia are metabolically differentiated and not clonal. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.27.568837. [PMID: 38076927 PMCID: PMC10705294 DOI: 10.1101/2023.11.27.568837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/24/2023]
Abstract
Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing eight new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nano-scale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal non-canonical amino acid tagging (BONCAT) we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.
Collapse
|
4
|
Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea. PLoS Biol 2023; 21:e3002292. [PMID: 37747940 PMCID: PMC10553843 DOI: 10.1371/journal.pbio.3002292] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Revised: 10/05/2023] [Accepted: 08/08/2023] [Indexed: 09/27/2023] Open
Abstract
Sulfate-coupled anaerobic oxidation of methane (AOM) is performed by multicellular consortia of anaerobic methanotrophic archaea (ANME) in obligate syntrophic partnership with sulfate-reducing bacteria (SRB). Diverse ANME and SRB clades co-associate but the physiological basis for their adaptation and diversification is not well understood. In this work, we used comparative metagenomics and phylogenetics to investigate the metabolic adaptation among the 4 main syntrophic SRB clades (HotSeep-1, Seep-SRB2, Seep-SRB1a, and Seep-SRB1g) and identified features associated with their syntrophic lifestyle that distinguish them from their non-syntrophic evolutionary neighbors in the phylum Desulfobacterota. We show that the protein complexes involved in direct interspecies electron transfer (DIET) from ANME to the SRB outer membrane are conserved between the syntrophic lineages. In contrast, the proteins involved in electron transfer within the SRB inner membrane differ between clades, indicative of convergent evolution in the adaptation to a syntrophic lifestyle. Our analysis suggests that in most cases, this adaptation likely occurred after the acquisition of the DIET complexes in an ancestral clade and involve horizontal gene transfers within pathways for electron transfer (CbcBA) and biofilm formation (Pel). We also provide evidence for unique adaptations within syntrophic SRB clades, which vary depending on the archaeal partner. Among the most widespread syntrophic SRB, Seep-SRB1a, subclades that specifically partner ANME-2a are missing the cobalamin synthesis pathway, suggestive of nutritional dependency on its partner, while closely related Seep-SRB1a partners of ANME-2c lack nutritional auxotrophies. Our work provides insight into the features associated with DIET-based syntrophy and the adaptation of SRB towards it.
Collapse
|
5
|
A standardized quantitative analysis strategy for stable isotope probing metagenomics. mSystems 2023; 8:e0128022. [PMID: 37377419 PMCID: PMC10469821 DOI: 10.1128/msystems.01280-22] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Accepted: 04/19/2023] [Indexed: 06/29/2023] Open
Abstract
Stable isotope probing (SIP) facilitates culture-independent identification of active microbial populations within complex ecosystems through isotopic enrichment of nucleic acids. Many DNA-SIP studies rely on 16S rRNA gene sequences to identify active taxa, but connecting these sequences to specific bacterial genomes is often challenging. Here, we describe a standardized laboratory and analysis framework to quantify isotopic enrichment on a per-genome basis using shotgun metagenomics instead of 16S rRNA gene sequencing. To develop this framework, we explored various sample processing and analysis approaches using a designed microbiome where the identity of labeled genomes and their level of isotopic enrichment were experimentally controlled. With this ground truth dataset, we empirically assessed the accuracy of different analytical models for identifying active taxa and examined how sequencing depth impacts the detection of isotopically labeled genomes. We also demonstrate that using synthetic DNA internal standards to measure absolute genome abundances in SIP density fractions improves estimates of isotopic enrichment. In addition, our study illustrates the utility of internal standards to reveal anomalies in sample handling that could negatively impact SIP metagenomic analyses if left undetected. Finally, we present SIPmg, an R package to facilitate the estimation of absolute abundances and perform statistical analyses for identifying labeled genomes within SIP metagenomic data. This experimentally validated analysis framework strengthens the foundation of DNA-SIP metagenomics as a tool for accurately measuring the in situ activity of environmental microbial populations and assessing their genomic potential. IMPORTANCE Answering the questions, "who is eating what?" and "who is active?" within complex microbial communities is paramount for our ability to model, predict, and modulate microbiomes for improved human and planetary health. These questions can be pursued using stable isotope probing to track the incorporation of labeled compounds into cellular DNA during microbial growth. However, with traditional stable isotope methods, it is challenging to establish links between an active microorganism's taxonomic identity and genome composition while providing quantitative estimates of the microorganism's isotope incorporation rate. Here, we report an experimental and analytical workflow that lays the foundation for improved detection of metabolically active microorganisms and better quantitative estimates of genome-resolved isotope incorporation, which can be used to further refine ecosystem-scale models for carbon and nutrient fluxes within microbiomes.
Collapse
|
6
|
Genome-wide fitness profiling reveals molecular mechanisms that bacteria use to interact with Trichoderma atroviride exometabolites. PLoS Genet 2023; 19:e1010909. [PMID: 37651474 PMCID: PMC10516422 DOI: 10.1371/journal.pgen.1010909] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2023] [Revised: 09/22/2023] [Accepted: 08/07/2023] [Indexed: 09/02/2023] Open
Abstract
Trichoderma spp. are ubiquitous rhizosphere fungi capable of producing several classes of secondary metabolites that can modify the dynamics of the plant-associated microbiome. However, the bacterial-fungal mechanisms that mediate these interactions have not been fully characterized. Here, a random barcode transposon-site sequencing (RB-TnSeq) approach was employed to identify bacterial genes important for fitness in the presence of Trichoderma atroviride exudates. We selected three rhizosphere bacteria with RB-TnSeq mutant libraries that can promote plant growth: the nitrogen fixers Klebsiella michiganensis M5aI and Herbaspirillum seropedicae SmR1, and Pseudomonas simiae WCS417. As a non-rhizosphere species, Pseudomonas putida KT2440 was also included. From the RB-TnSeq data, nitrogen-fixing bacteria competed mainly for iron and required the siderophore transport system TonB/ExbB for optimal fitness in the presence of T. atroviride exudates. In contrast, P. simiae and P. putida were highly dependent on mechanisms associated with membrane lipid modification that are required for resistance to cationic antimicrobial peptides (CAMPs). A mutant in the Hog1-MAP kinase (Δtmk3) gene of T. atroviride showed altered expression patterns of many nonribosomal peptide synthetase (NRPS) biosynthetic gene clusters with potential antibiotic activity. In contrast to exudates from wild-type T. atroviride, bacterial mutants containing lesions in genes associated with resistance to antibiotics did not show fitness defects when RB-TnSeq libraries were exposed to exudates from the Δtmk3 mutant. Unexpectedly, exudates from wild-type T. atroviride and the Δtmk3 mutant rescued purine auxotrophic mutants of H. seropedicae, K. michiganensis and P. simiae. Metabolomic analysis on exudates from wild-type T. atroviride and the Δtmk3 mutant showed that both strains excrete purines and complex metabolites; functional Tmk3 is required to produce some of these metabolites. This study highlights the complex interplay between Trichoderma-metabolites and soil bacteria, revealing both beneficial and antagonistic effects, and underscoring the intricate and multifaceted nature of this relationship.
Collapse
|
7
|
BONCAT-FACS-Seq reveals the active fraction of a biocrust community undergoing a wet-up event. Front Microbiol 2023; 14:1176751. [PMID: 37434715 PMCID: PMC10330726 DOI: 10.3389/fmicb.2023.1176751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 06/05/2023] [Indexed: 07/13/2023] Open
Abstract
Determining which microorganisms are active within soil communities remains a major technical endeavor in microbial ecology research. One promising method to accomplish this is coupling bioorthogonal non-canonical amino acid tagging (BONCAT) with fluorescence activated cell sorting (FACS) which sorts cells based on whether or not they are producing new proteins. Combined with shotgun metagenomic sequencing (Seq), we apply this method to profile the diversity and potential functional capabilities of both active and inactive microorganisms in a biocrust community after being resuscitated by a simulated rain event. We find that BONCAT-FACS-Seq is capable of discerning the pools of active and inactive microorganisms, especially within hours of applying the BONCAT probe. The active and inactive components of the biocrust community differed in species richness and composition at both 4 and 21 h after the wetting event. The active fraction of the biocrust community is marked by taxa commonly observed in other biocrust communities, many of which play important roles in species interactions and nutrient transformations. Among these, 11 families within the Firmicutes are enriched in the active fraction, supporting previous reports indicating that the Firmicutes are key early responders to biocrust wetting. We highlight the apparent inactivity of many Actinobacteria and Proteobacteria through 21 h after wetting, and note that members of the Chitinophagaceae, enriched in the active fraction, may play important ecological roles following wetting. Based on the enrichment of COGs in the active fraction, predation by phage and other bacterial members, as well as scavenging and recycling of labile nutrients, appear to be important ecological processes soon after wetting. To our knowledge, this is the first time BONCAT-FACS-Seq has been applied to biocrust samples, and therefore we discuss the potential advantages and shortcomings of coupling metagenomics to BONCAT to intact soil communities such as biocrust. In all, by pairing BONCAT-FACS and metagenomics, we are capable of highlighting the taxa and potential functions that typifies the microbes actively responding to a rain event.
Collapse
|
8
|
Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.06.543983. [PMID: 37333348 PMCID: PMC10274767 DOI: 10.1101/2023.06.06.543983] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/20/2023]
Abstract
Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.
Collapse
|
9
|
HT-SIP: a semi-automated stable isotope probing pipeline identifies cross-kingdom interactions in the hyphosphere of arbuscular mycorrhizal fungi. MICROBIOME 2022; 10:199. [PMID: 36434737 PMCID: PMC9700909 DOI: 10.1186/s40168-022-01391-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Accepted: 10/04/2022] [Indexed: 06/16/2023]
Abstract
BACKGROUND Linking the identity of wild microbes with their ecophysiological traits and environmental functions is a key ambition for microbial ecologists. Of many techniques that strive for this goal, Stable-isotope probing-SIP-remains among the most comprehensive for studying whole microbial communities in situ. In DNA-SIP, actively growing microorganisms that take up an isotopically heavy substrate build heavier DNA, which can be partitioned by density into multiple fractions and sequenced. However, SIP is relatively low throughput and requires significant hands-on labor. We designed and tested a semi-automated, high-throughput SIP (HT-SIP) pipeline to support well-replicated, temporally resolved amplicon and metagenomics experiments. We applied this pipeline to a soil microhabitat with significant ecological importance-the hyphosphere zone surrounding arbuscular mycorrhizal fungal (AMF) hyphae. AMF form symbiotic relationships with most plant species and play key roles in terrestrial nutrient and carbon cycling. RESULTS Our HT-SIP pipeline for fractionation, cleanup, and nucleic acid quantification of density gradients requires one-sixth of the hands-on labor compared to manual SIP and allows 16 samples to be processed simultaneously. Automated density fractionation increased the reproducibility of SIP gradients compared to manual fractionation, and we show adding a non-ionic detergent to the gradient buffer improved SIP DNA recovery. We applied HT-SIP to 13C-AMF hyphosphere DNA from a 13CO2 plant labeling study and created metagenome-assembled genomes (MAGs) using high-resolution SIP metagenomics (14 metagenomes per gradient). SIP confirmed the AMF Rhizophagus intraradices and associated MAGs were highly enriched (10-33 atom% 13C), even though the soils' overall enrichment was low (1.8 atom% 13C). We assembled 212 13C-hyphosphere MAGs; the hyphosphere taxa that assimilated the most AMF-derived 13C were from the phyla Myxococcota, Fibrobacterota, Verrucomicrobiota, and the ammonia-oxidizing archaeon genus Nitrososphaera. CONCLUSIONS Our semi-automated HT-SIP approach decreases operator time and improves reproducibility by targeting the most labor-intensive steps of SIP-fraction collection and cleanup. We illustrate this approach in a unique and understudied soil microhabitat-generating MAGs of actively growing microbes living in the AMF hyphosphere (without plant roots). The MAGs' phylogenetic composition and gene content suggest predation, decomposition, and ammonia oxidation may be key processes in hyphosphere nutrient cycling. Video Abstract.
Collapse
|
10
|
A centimeter-long bacterium with DNA contained in metabolically active, membrane-bound organelles. Science 2022; 376:1453-1458. [PMID: 35737788 DOI: 10.1126/science.abb3634] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Cells of most bacterial species are around 2 micrometers in length, with some of the largest specimens reaching 750 micrometers. Using fluorescence, x-ray, and electron microscopy in conjunction with genome sequencing, we characterized Candidatus (Ca.) Thiomargarita magnifica, a bacterium that has an average cell length greater than 9000 micrometers and is visible to the naked eye. These cells grow orders of magnitude over theoretical limits for bacterial cell size, display unprecedented polyploidy of more than half a million copies of a very large genome, and undergo a dimorphic life cycle with asymmetric segregation of chromosomes into daughter cells. These features, along with compartmentalization of genomic material and ribosomes in translationally active organelles bound by bioenergetic membranes, indicate gain of complexity in the Thiomargarita lineage and challenge traditional concepts of bacterial cells.
Collapse
|
11
|
Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution. THE ISME JOURNAL 2022; 16:1337-1347. [PMID: 34969995 PMCID: PMC9039060 DOI: 10.1038/s41396-021-01178-4] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/29/2021] [Revised: 11/29/2021] [Accepted: 12/10/2021] [Indexed: 02/07/2023]
Abstract
With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.
Collapse
|
12
|
Subsurface hydrocarbon degradation strategies in low- and high-sulfate coal seam communities identified with activity-based metagenomics. NPJ Biofilms Microbiomes 2022; 8:7. [PMID: 35177633 PMCID: PMC8854433 DOI: 10.1038/s41522-022-00267-2] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Accepted: 01/05/2022] [Indexed: 11/18/2022] Open
Abstract
Environmentally relevant metagenomes and BONCAT-FACS derived translationally active metagenomes from Powder River Basin coal seams were investigated to elucidate potential genes and functional groups involved in hydrocarbon degradation to methane in coal seams with high- and low-sulfate levels. An advanced subsurface environmental sampler allowed the establishment of coal-associated microbial communities under in situ conditions for metagenomic analyses from environmental and translationally active populations. Metagenomic sequencing demonstrated that biosurfactants, aerobic dioxygenases, and anaerobic phenol degradation pathways were present in active populations across the sampled coal seams. In particular, results suggested the importance of anaerobic degradation pathways under high-sulfate conditions with an emphasis on fumarate addition. Under low-sulfate conditions, a mixture of both aerobic and anaerobic pathways was observed but with a predominance of aerobic dioxygenases. The putative low-molecular-weight biosurfactant, lichysein, appeared to play a more important role compared to rhamnolipids. The methods used in this study—subsurface environmental samplers in combination with metagenomic sequencing of both total and translationally active metagenomes—offer a deeper and environmentally relevant perspective on community genetic potential from coal seams poised at different redox conditions broadening the understanding of degradation strategies for subsurface carbon.
Collapse
|
13
|
Metapangenomics reveals depth-dependent shifts in metabolic potential for the ubiquitous marine bacterial SAR324 lineage. MICROBIOME 2021; 9:172. [PMID: 34389059 PMCID: PMC8364033 DOI: 10.1186/s40168-021-01119-5] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Accepted: 06/22/2021] [Indexed: 06/01/2023]
Abstract
BACKGROUND Oceanic microbiomes play a pivotal role in the global carbon cycle and are central to the transformation and recycling of carbon and energy in the ocean's interior. SAR324 is a ubiquitous but poorly understood uncultivated clade of Deltaproteobacteria that inhabits the entire water column, from ocean surface waters to its deep interior. Although some progress has been made in elucidating potential metabolic traits of SAR324 in the dark ocean, very little is known about the ecology and the metabolic capabilities of this group in the euphotic and twilight zones. To investigate the comparative genomics, ecology, and physiological potential of the SAR324 clade, we examined the distribution and variability of key genomic features and metabolic pathways in this group from surface waters to the abyss in the North Pacific Subtropical Gyre, one of the largest biomes on Earth. RESULTS We leveraged a pangenomic ecological approach, combining spatio-temporally resolved single-amplified genome, metagenomic, and metatranscriptomic datasets. The data revealed substantial genomic diversity throughout the SAR324 clade, with distinct depth and temporal distributions that clearly differentiated ecotypes. Phylogenomic subclade delineation, environmental distributions, genomic feature similarities, and metabolic capacities revealed strong congruence. The four SAR324 ecotypes delineated in this study revealed striking divergence from one another with respect to their habitat-specific metabolic potentials. The ecotypes living in the dark or twilight oceans shared genomic features and metabolic capabilities consistent with a sulfur-based chemolithoautotrophic lifestyle. In contrast, those inhabiting the sunlit ocean displayed higher plasticity energy-related metabolic pathways, supporting a presumptive photoheterotrophic lifestyle. In epipelagic SAR324 ecotypes, we observed the presence of two types of proton-pumping rhodopsins, as well as genomic, transcriptomic, and ecological evidence for active photoheterotrophy, based on xanthorhodopsin-like light-harvesting proteins. CONCLUSIONS Combining pangenomic and both metagenomic and metatranscriptomic profiling revealed a striking divergence in the vertical distribution, genomic composition, metabolic potential, and predicted lifestyle strategies of geographically co-located members of the SAR324 bacterial clade. The results highlight the utility of metapangenomic approaches employed across environmental gradients, to decipher the properties and variation in function and ecological traits of specific phylogenetic clades within complex microbiomes. Video abstract.
Collapse
|
14
|
Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design. mSystems 2020; 5:e00151-20. [PMID: 32694124 PMCID: PMC7566279 DOI: 10.1128/msystems.00151-20] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2020] [Accepted: 06/29/2020] [Indexed: 12/14/2022] Open
Abstract
Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets-from soil and seawater microbiomes-to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon's abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.IMPORTANCE One of the biggest challenges in microbial ecology is correlating the identity of microorganisms with the roles they fulfill in natural environmental systems. Studies of microbes in pure culture reveal much about their genomic content and potential functions but may not reflect an organism's activity within its natural community. Culture-independent studies supply a community-wide view of composition and function in the context of community interactions but often fail to link the two. Quantitative stable isotope probing (qSIP) is a method that can link the identity and functional activity of specific microbes within a naturally occurring community. Here, we explore how the resolution of density gradient fractionation affects the error and precision of qSIP results, how they may be improved via additional experimental replication, and discuss cost-benefit balanced scenarios for SIP experimental design.
Collapse
|
15
|
Ecological and genomic analyses of candidate phylum WPS-2 bacteria in an unvegetated soil. Environ Microbiol 2020; 22:3143-3157. [PMID: 32372527 DOI: 10.1111/1462-2920.15054] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2019] [Revised: 04/28/2020] [Accepted: 04/29/2020] [Indexed: 12/01/2022]
Abstract
Members of the bacterial candidate phylum WPS-2 (or Eremiobacterota) are abundant in several dry, bare soil environments. In a bare soil deposited by an extinct iron-sulfur spring, we found that WPS-2 comprised up to 24% of the bacterial community and up to 108 cells per g of soil based on 16S rRNA gene sequencing and quantification. A single genus-level cluster (Ca. Rubrimentiphilum) predominated in bare soils but was less abundant in adjacent forest. Nearly complete genomes of Ca. Rubrimentiphilum were recovered as single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs). Surprisingly, given the abundance of WPS-2 in bare soils, the genomes did not indicate any capacity for autotrophy, phototrophy, or trace gas metabolism. Instead, they suggest a predominantly aerobic organoheterotrophic lifestyle, perhaps based on scavenging amino acids, nucleotides, and complex oligopeptides, along with lithotrophic capacity on thiosulfate. Network analyses of the entire community showed that some species of Chloroflexi, Actinobacteria, and candidate phylum AD3 (or Dormibacterota) co-occurred with Ca. Rubrimentiphilum and may represent ecological or metabolic partners. We propose that Ca. Rubrimentiphilum act as efficient heterotrophic scavengers. Combined with previous studies, these data suggest that the phylum WPS-2 includes bacteria with diverse metabolic capabilities.
Collapse
|
16
|
Niche differentiation is spatially and temporally regulated in the rhizosphere. THE ISME JOURNAL 2020; 14:999-1014. [PMID: 31953507 PMCID: PMC7082339 DOI: 10.1038/s41396-019-0582-x] [Citation(s) in RCA: 69] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/28/2019] [Revised: 10/28/2019] [Accepted: 12/18/2019] [Indexed: 01/14/2023]
Abstract
The rhizosphere is a hotspot for microbial carbon transformations, and is the entry point for root polysaccharides and polymeric carbohydrates that are important precursors to soil organic matter (SOM). However, the ecological mechanisms that underpin rhizosphere carbohydrate depolymerization are poorly understood. Using Avena fatua, a common annual grass, we analyzed time-resolved metatranscriptomes to compare microbial functions in rhizosphere, detritusphere, and combined rhizosphere-detritusphere habitats. Transcripts were binned using a unique reference database generated from soil isolate genomes, single-cell amplified genomes, metagenomes, and stable isotope probing metagenomes. While soil habitat significantly affected both community composition and overall gene expression, the succession of microbial functions occurred at a faster time scale than compositional changes. Using hierarchical clustering of upregulated decomposition genes, we identified four distinct microbial guilds populated by taxa whose functional succession patterns suggest specialization for substrates provided by fresh growing roots, decaying root detritus, the combination of live and decaying root biomass, or aging root material. Carbohydrate depolymerization genes were consistently upregulated in the rhizosphere, and both taxonomic and functional diversity were highest in the combined rhizosphere-detritusphere, suggesting coexistence of rhizosphere guilds is facilitated by niche differentiation. Metatranscriptome-defined guilds provide a framework to model rhizosphere succession and its consequences for soil carbon cycling.
Collapse
|
17
|
A pipeline for targeted metagenomics of environmental bacteria. MICROBIOME 2020; 8:21. [PMID: 32061258 PMCID: PMC7024552 DOI: 10.1186/s40168-020-0790-7] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Accepted: 01/19/2020] [Indexed: 05/19/2023]
Abstract
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.
Collapse
|
18
|
Abstract
The ability to link soil microbial diversity to soil processes requires technologies that differentiate active microbes from extracellular DNA and dormant cells. Here, we use BONCAT (bioorthogonal non-canonical amino acid tagging) to measure translationally active cells in soils. We compare the active population of two soil depths from Oak Ridge (Tennessee, USA) and find that a maximum of 25-70% of the extractable cells are active. Analysis of 16S rRNA sequences from BONCAT-positive cells recovered by fluorescence-activated cell sorting (FACS) reveals that the phylogenetic composition of the active fraction is distinct from the total population of extractable cells. Some members of the community are found to be active at both depths independently of their abundance rank, suggesting that the incubation conditions favor the activity of similar organisms. We conclude that BONCAT-FACS is effective for interrogating the active fraction of soil microbiomes in situ and provides a new approach for uncovering the links between soil processes and specific microbial groups.
Collapse
|
19
|
Optimizing de novo genome assembly from PCR-amplified metagenomes. PeerJ 2019; 7:e6902. [PMID: 31119088 PMCID: PMC6511391 DOI: 10.7717/peerj.6902] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2018] [Accepted: 04/03/2019] [Indexed: 12/13/2022] Open
Abstract
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.
Collapse
|
20
|
Freshwater carbon and nutrient cycles revealed through reconstructed population genomes. PeerJ 2018; 6:e6075. [PMID: 30581671 PMCID: PMC6292386 DOI: 10.7717/peerj.6075] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Accepted: 11/05/2018] [Indexed: 02/01/2023] Open
Abstract
Although microbes mediate much of the biogeochemical cycling in freshwater, the categories of carbon and nutrients currently used in models of freshwater biogeochemical cycling are too broad to be relevant on a microbial scale. One way to improve these models is to incorporate microbial data. Here, we analyze both genes and genomes from three metagenomic time series and propose specific roles for microbial taxa in freshwater biogeochemical cycles. Our metagenomic time series span multiple years and originate from a eutrophic lake (Lake Mendota) and a humic lake (Trout Bog Lake) with contrasting water chemistry. Our analysis highlights the role of polyamines in the nitrogen cycle, the diversity of diazotrophs between lake types, the balance of assimilatory vs. dissimilatory sulfate reduction in freshwater, the various associations between types of phototrophy and carbon fixation, and the density and diversity of glycoside hydrolases in freshwater microbes. We also investigated aspects of central metabolism such as hydrogen metabolism, oxidative phosphorylation, methylotrophy, and sugar degradation. Finally, by analyzing the dynamics over time in nitrogen fixation genes and Cyanobacteria genomes, we show that the potential for nitrogen fixation is linked to specific populations in Lake Mendota. This work represents an important step towards incorporating microbial data into ecosystem models and provides a better understanding of how microbes may participate in freshwater biogeochemical cycling.
Collapse
|
21
|
Abstract
Known giant virus diversity is currently skewed towards viruses isolated from aquatic environments and cultivated in the laboratory. Here, we employ cultivation-independent metagenomics and mini-metagenomics on soils from the Harvard Forest, leading to the discovery of 16 novel giant viruses, chiefly recovered by mini-metagenomics. The candidate viruses greatly expand phylogenetic diversity of known giant viruses and either represented novel lineages or are affiliated with klosneuviruses, Cafeteria roenbergensis virus or tupanviruses. One assembled genome with a size of 2.4 Mb represents the largest currently known viral genome in the Mimiviridae, and others encode up to 80% orphan genes. In addition, we find more than 240 major capsid proteins encoded on unbinned metagenome fragments, further indicating that giant viruses are underexplored in soil ecosystems. The fact that most of these novel viruses evaded detection in bulk metagenomes suggests that mini-metagenomics could be a valuable approach to unearth viral giants.
Collapse
|
22
|
Single-cell genomics of co-sorted Nanoarchaeota suggests novel putative host associations and diversification of proteins involved in symbiosis. MICROBIOME 2018; 6:161. [PMID: 30223889 PMCID: PMC6142677 DOI: 10.1186/s40168-018-0539-8] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/19/2018] [Accepted: 08/26/2018] [Indexed: 05/20/2023]
Abstract
BACKGROUND Nanoarchaeota are obligate symbionts of other Archaea first discovered 16 years ago, yet little is known about this largely uncultivated taxon. While Nanoarchaeota diversity has been detected in a variety of habitats using 16S rRNA gene surveys, genome sequences have been available for only three Nanoarchaeota and their hosts. The host range and adaptation of Nanoarchaeota to a wide range of environmental conditions has thus largely remained elusive. Single-cell genomics is an ideal approach to address these questions as Nanoarchaeota can be isolated while still attached to putative hosts, enabling the exploration of cell-cell interactions and fine-scale genomic diversity. RESULTS From 22 single amplified genomes (SAGs) from three hot springs in Yellowstone National Park, we derived a genome-based phylogeny of the phylum Nanoarchaeota, linking it to global 16S rRNA gene diversity. By exploiting sequencing of co-sorted tightly attached cells, we associated Nanoarchaeota with 6 novel putative hosts, 2 of which were found in multiple SAGs, and showed that the same host species may associate with multiple species of Nanoarchaeota. Comparison of single nucleotide polymorphisms (SNPs) within a population of Nanoarchaeota SAGs indicated that Nanoarchaeota attached to a single host cell in situ are likely clonal. In addition to an overall pattern of purifying selection, we found significantly higher densities of non-synonymous SNPs in hypothetical cell surface proteins, as compared to other functional categories. Genes implicated in interactions in other obligate microbe-microbe symbioses, including those encoding a cytochrome bd-I ubiquinol oxidase and a FlaJ/TadC homologue possibly involved in type IV pili production, also had relatively high densities of non-synonymous SNPs. CONCLUSIONS This population genetics study of Nanoarchaeota greatly expands the known potential host range of the phylum and hints at what genes may be involved in adaptation to diverse environments or different hosts. We provide the first evidence that Nanoarchaeota cells attached to the same host cell are clonal and propose a hypothesis for how clonality may occur despite diverse symbiont populations.
Collapse
|
23
|
Single cell genomic and transcriptomic evidence for the use of alternative nitrogen substrates by anammox bacteria. ISME JOURNAL 2018; 12:2706-2722. [PMID: 29991764 DOI: 10.1038/s41396-018-0223-9] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2017] [Revised: 03/20/2018] [Accepted: 03/29/2018] [Indexed: 11/09/2022]
Abstract
Anaerobic ammonium oxidation (anammox) contributes substantially to ocean nitrogen loss, particularly in anoxic marine zones (AMZs). Ammonium is scarce in AMZs, raising the hypothesis that organic nitrogen compounds may be ammonium sources for anammox. Biochemical measurements suggest that the organic compounds urea and cyanate can support anammox in AMZs. However, it is unclear if anammox bacteria degrade these compounds to ammonium themselves, or rely on other organisms for this process. Genes for urea degradation have not been found in anammox bacteria, and genomic evidence for cyanate use for anammox is limited to a cyanase gene recovered from the sediment bacterium Candidatus Scalindua profunda. Here, analysis of Ca. Scalindua single amplified genomes from the Eastern Tropical North Pacific AMZ revealed genes for urea degradation and transport, as well as for cyanate degradation. Urease and cyanase genes were transcribed, along with anammox genes, in the AMZ core where anammox rates peaked. Homologs of these genes were also detected in meta-omic datasets from major AMZs in the Eastern Tropical South Pacific and Arabian Sea. These results suggest that anammox bacteria from different ocean regions can directly access organic nitrogen substrates. Future studies should assess if and under what environmental conditions these substrates contribute to the ammonium budget for anammox.
Collapse
|
24
|
Corrigendum: Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol 2018; 36:660. [PMID: 29979671 PMCID: PMC7608355 DOI: 10.1038/nbt0718-660a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
25
|
Corrigendum: Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol 2018; 36:196. [PMID: 29406516 PMCID: PMC7609277 DOI: 10.1038/nbt0218-196a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|
26
|
Genome-wide identification of bacterial plant colonization genes. PLoS Biol 2017; 15:e2002860. [PMID: 28938018 PMCID: PMC5627942 DOI: 10.1371/journal.pbio.2002860] [Citation(s) in RCA: 115] [Impact Index Per Article: 16.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2017] [Revised: 10/04/2017] [Accepted: 09/05/2017] [Indexed: 01/09/2023] Open
Abstract
Diverse soil-resident bacteria can contribute to plant growth and health, but the molecular mechanisms enabling them to effectively colonize their plant hosts remain poorly understood. We used randomly barcoded transposon mutagenesis sequencing (RB-TnSeq) in Pseudomonas simiae, a model root-colonizing bacterium, to establish a genome-wide map of bacterial genes required for colonization of the Arabidopsis thaliana root system. We identified 115 genes (2% of all P. simiae genes) with functions that are required for maximal competitive colonization of the root system. Among the genes we identified were some with obvious colonization-related roles in motility and carbon metabolism, as well as 44 other genes that had no or vague functional predictions. Independent validation assays of individual genes confirmed colonization functions for 20 of 22 (91%) cases tested. To further characterize genes identified by our screen, we compared the functional contributions of P. simiae genes to growth in 90 distinct in vitro conditions by RB-TnSeq, highlighting specific metabolic functions associated with root colonization genes. Our analysis of bacterial genes by sequence-driven saturation mutagenesis revealed a genome-wide map of the genetic determinants of plant root colonization and offers a starting point for targeted improvement of the colonization capabilities of plant-beneficial microbes. Plants fix carbon to create an abundance of sugars and amino acids, thus providing an enticing environment for microorganisms that reside in soil. Once these microorganisms have colonized the root environment, they can dramatically influence plant growth and development. We set out to identify a comprehensive set of microbial genes that control or influence root colonization, using a genome-wide transposon mutagenesis approach (randomly barcoded transposon sequencing [RB-TnSeq]). By using this method, we identified several hundred genes that, when mutated, affect the ability of the bacterium P. simiae to competitively colonize the root system of the model plant A. thaliana. These included many genes purported to be involved in carbohydrate metabolism, cell wall biosynthesis, and motility, underscoring the notion that sugar metabolism, defense, and motility are all key features of a root-colonizing microbe. We also identified several amino acid transport and metabolism genes with mutations that confer a fitness advantage in root colonization. Lastly, we identified several genes with no known function that significantly alter root colonization ability when mutated. These findings suggest novel engineering strategies to improve biological product development, and will facilitate the mechanistic exploration of the root colonization process.
Collapse
|
27
|
Ecophysiology of Freshwater Verrucomicrobia Inferred from Metagenome-Assembled Genomes. mSphere 2017; 2:e00277-17. [PMID: 28959738 PMCID: PMC5615132 DOI: 10.1128/msphere.00277-17] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Accepted: 09/05/2017] [Indexed: 11/20/2022] Open
Abstract
Microbes are critical in carbon and nutrient cycling in freshwater ecosystems. Members of the Verrucomicrobia are ubiquitous in such systems, and yet their roles and ecophysiology are not well understood. In this study, we recovered 19 Verrucomicrobia draft genomes by sequencing 184 time-series metagenomes from a eutrophic lake and a humic bog that differ in carbon source and nutrient availabilities. These genomes span four of the seven previously defined Verrucomicrobia subdivisions and greatly expand knowledge of the genomic diversity of freshwater Verrucomicrobia. Genome analysis revealed their potential role as (poly)saccharide degraders in freshwater, uncovered interesting genomic features for this lifestyle, and suggested their adaptation to nutrient availabilities in their environments. Verrucomicrobia populations differ significantly between the two lakes in glycoside hydrolase gene abundance and functional profiles, reflecting the autochthonous and terrestrially derived allochthonous carbon sources of the two ecosystems, respectively. Interestingly, a number of genomes recovered from the bog contained gene clusters that potentially encode a novel porin-multiheme cytochrome c complex and might be involved in extracellular electron transfer in the anoxic humus-rich environment. Notably, most epilimnion genomes have large numbers of so-called "Planctomycete-specific" cytochrome c-encoding genes, which exhibited distribution patterns nearly opposite to those seen with glycoside hydrolase genes, probably associated with the different levels of environmental oxygen availability and carbohydrate complexity between lakes/layers. Overall, the recovered genomes represent a major step toward understanding the role, ecophysiology, and distribution of Verrucomicrobia in freshwater. IMPORTANCE Freshwater Verrucomicrobia spp. are cosmopolitan in lakes and rivers, and yet their roles and ecophysiology are not well understood, as cultured freshwater Verrucomicrobia spp. are restricted to one subdivision of this phylum. Here, we greatly expanded the known genomic diversity of this freshwater lineage by recovering 19 Verrucomicrobia draft genomes from 184 metagenomes collected from a eutrophic lake and a humic bog across multiple years. Most of these genomes represent the first freshwater representatives of several Verrucomicrobia subdivisions. Genomic analysis revealed Verrucomicrobia to be potential (poly)saccharide degraders and suggested their adaptation to carbon sources of different origins in the two contrasting ecosystems. We identified putative extracellular electron transfer genes and so-called "Planctomycete-specific" cytochrome c-encoding genes and identified their distinct distribution patterns between the lakes/layers. Overall, our analysis greatly advances the understanding of the function, ecophysiology, and distribution of freshwater Verrucomicrobia, while highlighting their potential role in freshwater carbon cycling.
Collapse
|
28
|
Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol 2017; 35:725-731. [PMID: 28787424 PMCID: PMC6436528 DOI: 10.1038/nbt.3893] [Citation(s) in RCA: 991] [Impact Index Per Article: 141.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2016] [Accepted: 04/27/2017] [Indexed: 12/20/2022]
Abstract
Standards for sequencing the microbial 'uncultivated majority', namely bacterial and archaeal single-cell genome sequences, and genome sequences from metagenomic datasets, are proposed. We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.
Collapse
|
29
|
Abstract
Bacteria of the SAR11 clade constitute up to one half of all microbial cells in the oxygen-rich surface ocean. SAR11 bacteria are also abundant in oxygen minimum zones (OMZs), where oxygen falls below detection and anaerobic microbes have vital roles in converting bioavailable nitrogen to N2 gas. Anaerobic metabolism has not yet been observed in SAR11, and it remains unknown how these bacteria contribute to OMZ biogeochemical cycling. Here, genomic analysis of single cells from the world's largest OMZ revealed previously uncharacterized SAR11 lineages with adaptations for life without oxygen, including genes for respiratory nitrate reductases (Nar). SAR11 nar genes were experimentally verified to encode proteins catalysing the nitrite-producing first step of denitrification and constituted ~40% of OMZ nar transcripts, with transcription peaking in the anoxic zone of maximum nitrate reduction activity. These results link SAR11 to pathways of ocean nitrogen loss, redefining the ecological niche of Earth's most abundant organismal group.
Collapse
|
30
|
Genome-wide selective sweeps and gene-specific sweeps in natural bacterial populations. THE ISME JOURNAL 2016; 10:1589-601. [PMID: 26744812 PMCID: PMC4918448 DOI: 10.1038/ismej.2015.241] [Citation(s) in RCA: 142] [Impact Index Per Article: 17.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2015] [Revised: 10/22/2015] [Accepted: 11/18/2015] [Indexed: 12/30/2022]
Abstract
Multiple models describe the formation and evolution of distinct microbial phylogenetic groups. These evolutionary models make different predictions regarding how adaptive alleles spread through populations and how genetic diversity is maintained. Processes predicted by competing evolutionary models, for example, genome-wide selective sweeps vs gene-specific sweeps, could be captured in natural populations using time-series metagenomics if the approach were applied over a sufficiently long time frame. Direct observations of either process would help resolve how distinct microbial groups evolve. Here, from a 9-year metagenomic study of a freshwater lake (2005-2013), we explore changes in single-nucleotide polymorphism (SNP) frequencies and patterns of gene gain and loss in 30 bacterial populations. SNP analyses revealed substantial genetic heterogeneity within these populations, although the degree of heterogeneity varied by >1000-fold among populations. SNP allele frequencies also changed dramatically over time within some populations. Interestingly, nearly all SNP variants were slowly purged over several years from one population of green sulfur bacteria, while at the same time multiple genes either swept through or were lost from this population. These patterns were consistent with a genome-wide selective sweep in progress, a process predicted by the 'ecotype model' of speciation but not previously observed in nature. In contrast, other populations contained large, SNP-free genomic regions that appear to have swept independently through the populations prior to the study without purging diversity elsewhere in the genome. Evidence for both genome-wide and gene-specific sweeps suggests that different models of bacterial speciation may apply to different populations coexisting in the same environment.
Collapse
|
31
|
Large-scale maps of variable infection efficiencies in aquatic Bacteroidetes phage-host model systems. Environ Microbiol 2016; 18:3949-3961. [PMID: 27235779 DOI: 10.1111/1462-2920.13392] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2015] [Accepted: 05/21/2015] [Indexed: 01/21/2023]
Abstract
Microbes drive ecosystem functioning and their viruses modulate these impacts through mortality, gene transfer and metabolic reprogramming. Despite the importance of virus-host interactions and likely variable infection efficiencies of individual phages across hosts, such variability is seldom quantified. Here, we quantify infection efficiencies of 38 phages against 19 host strains in aquatic Cellulophaga (Bacteroidetes) phage-host model systems. Binary data revealed that some phages infected only one strain while others infected 17, whereas quantitative data revealed that efficiency of infection could vary 10 orders of magnitude, even among phages within one population. This provides a baseline for understanding and modeling intrapopulation host range variation. Genera specific host ranges were also informative. For example, the Cellulophaga Microviridae, showed a markedly broader intra-species host range than previously observed in Escherichia coli systems. Further, one phage genus, Cba41, was examined to investigate nonheritable changes in plating efficiency and burst size that depended on which host strain it most recently infected. While consistent with host modification of phage DNA, no differences in nucleotide sequence or DNA modifications were detected, leaving the observation repeatable, but the mechanism unresolved. Overall, this study highlights the importance of quantitatively considering replication variations in studies of phage-host interactions.
Collapse
|
32
|
Abstract
DNA methylation acts in concert with restriction enzymes to protect the integrity of prokaryotic genomes. Studies in a limited number of organisms suggest that methylation also contributes to prokaryotic genome regulation, but the prevalence and properties of such non-restriction-associated methylation systems remain poorly understood. Here, we used single molecule, real-time sequencing to map DNA modifications including m6A, m4C, and m5C across the genomes of 230 diverse bacterial and archaeal species. We observed DNA methylation in nearly all (93%) organisms examined, and identified a total of 834 distinct reproducibly methylated motifs. This data enabled annotation of the DNA binding specificities of 620 DNA Methyltransferases (MTases), doubling known specificities for previously hard to study Type I, IIG and III MTases, and revealing their extraordinary diversity. Strikingly, 48% of organisms harbor active Type II MTases with no apparent cognate restriction enzyme. These active ‘orphan’ MTases are present in diverse bacterial and archaeal phyla and show motif specificities and methylation patterns consistent with functions in gene regulation and DNA replication. Our results reveal the pervasive presence of DNA methylation throughout the prokaryotic kingdoms, as well as the diversity of sequence specificities and potential functions of DNA methylation systems. DNA methylation is a chemical modification of DNA present in many prokaryotic genomes. The best-known role of DNA methylation is as a component of restriction-modification systems. In these systems, restriction enzymes target foreign DNA for cleavage, while DNA methylation protects the host genome from destruction. Studies in a handful of organisms show that DNA methylation may also act independently of restriction systems and function in genome regulation. However, a lack of technologies has limited the study of DNA methylation to a small number of organisms, and the broader patterns and functions of DNA methylation remain unknown. Here we use SMRT-sequencing to determine the genome wide DNA methylation patterns of more than 200 diverse bacteria and archaea. We show that DNA methylation is pervasive and present in more than 90% of studied organisms. Analysis of this data enabled annotation of the specific DNA binding sites of more than 600 restriction systems, revealing their extraordinary diversity. Strikingly, we observed widespread DNA methylation in the absence of restriction systems. Analyses of these patterns reveal that they are conserved through evolution, and likely function in genome regulation. Thus DNA methylation may play a far wider function in prokaryotic genome biology than was previously supposed.
Collapse
|
33
|
Abstract
Extensive genomic diversity within coexisting members of a microbial species has been revealed through selected cultured isolates and metagenomic assemblies. Yet, the cell-by-cell genomic composition of wild uncultured populations of co-occurring cells is largely unknown. In this work, we applied large-scale single-cell genomics to study populations of the globally abundant marine cyanobacterium Prochlorococcus. We show that they are composed of hundreds of subpopulations with distinct "genomic backbones," each backbone consisting of a different set of core gene alleles linked to a small distinctive set of flexible genes. These subpopulations are estimated to have diverged at least a few million years ago, suggesting ancient, stable niche partitioning. Such a large set of coexisting subpopulations may be a general feature of free-living bacterial species with huge populations in highly mixed habitats.
Collapse
|
34
|
Ecology of uncultured Prochlorococcus clades revealed through single-cell genomics and biogeographic analysis. ISME JOURNAL 2012; 7:184-98. [PMID: 22895163 DOI: 10.1038/ismej.2012.89] [Citation(s) in RCA: 81] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Prochlorococcus is the numerically dominant photosynthetic organism throughout much of the world's oceans, yet little is known about the ecology and genetic diversity of populations inhabiting tropical waters. To help close this gap, we examined natural Prochlorococcus communities in the tropical Pacific Ocean using a single-cell whole-genome amplification and sequencing. Analysis of the gene content of just 10 single cells from these waters added 394 new genes to the Prochlorococcus pan-genome--that is, genes never before seen in a Prochlorococcus cell. Analysis of marker genes, including the ribosomal internal transcribed sequence, from dozens of individual cells revealed several representatives from two uncultivated clades of Prochlorococcus previously identified as HNLC1 and HNLC2. While the HNLC clades can dominate Prochlorococcus communities under certain conditions, their overall geographic distribution was highly restricted compared with other clades of Prochlorococcus. In the Atlantic and Pacific oceans, these clades were only found in warm waters with low Fe and high inorganic P levels. Genomic analysis suggests that at least one of these clades thrives in low Fe environments by scavenging organic-bound Fe, a process previously unknown in Prochlorococcus. Furthermore, the capacity to utilize organic-bound Fe appears to have been acquired horizontally and may be exchanged among other clades of Prochlorococcus. Finally, one of the single Prochlorococcus cells sequenced contained a partial genome of what appears to be a prophage integrated into the genome.
Collapse
|
35
|
Abstract
Background Different high-throughput nucleic acid sequencing platforms are currently available but a trade-off currently exists between the cost and number of reads that can be generated versus the read length that can be achieved. Methodology/Principal Findings We describe an experimental and computational pipeline yielding millions of reads that can exceed 200 bp with quality scores approaching that of traditional Sanger sequencing. The method combines an automatable gel-less library construction step with paired-end sequencing on a short-read instrument. With appropriately sized library inserts, mate-pair sequences can overlap, and we describe the SHERA software package that joins them to form a longer composite read. Conclusions/Significance This strategy is broadly applicable to sequencing applications that benefit from low-cost high-throughput sequencing, but require longer read lengths. We demonstrate that our approach enables metagenomic analyses using the Illumina Genome Analyzer, with low error rates, and at a fraction of the cost of pyrosequencing.
Collapse
|
36
|
Abstract
Background Single-cell genome sequencing has the potential to allow the in-depth exploration of the vast genetic diversity found in uncultured microbes. We used the marine cyanobacterium Prochlorococcus as a model system for addressing important challenges facing high-throughput whole genome amplification (WGA) and complete genome sequencing of individual cells. Methodology/Principal Findings We describe a pipeline that enables single-cell WGA on hundreds of cells at a time while virtually eliminating non-target DNA from the reactions. We further developed a post-amplification normalization procedure that mitigates extreme variations in sequencing coverage associated with multiple displacement amplification (MDA), and demonstrated that the procedure increased sequencing efficiency and facilitated genome assembly. We report genome recovery as high as 99.6% with reference-guided assembly, and 95% with de novo assembly starting from a single cell. We also analyzed the impact of chimera formation during MDA on de novo assembly, and discuss strategies to minimize the presence of incorrectly joined regions in contigs. Conclusions/Significance The methods describe in this paper will be useful for sequencing genomes of individual cells from a variety of samples.
Collapse
|
37
|
Assimilation of polysaccharides and glucose by major bacterial groups in the Delaware Estuary. Appl Environ Microbiol 2005; 71:7799-805. [PMID: 16332754 PMCID: PMC1317420 DOI: 10.1128/aem.71.12.7799-7805.2005] [Citation(s) in RCA: 99] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2005] [Accepted: 08/09/2005] [Indexed: 11/20/2022] Open
Abstract
The contribution of major bacterial groups to the assimilation of extracellular polymeric substances (EPS) and glucose in the Delaware Estuary was assessed using microautoradiography and fluorescence in situ hybridization. Bacterial groups contributed to EPS and glucose assimilation in part according to their distribution in the estuary. Abundance of the phylogenetic groups explained 35% and 55% of the variation in EPS and glucose assimilation, respectively. Actinobacteria contributed 70% to glucose assimilation in freshwater, while Alphaproteobacteria assimilated 60% of this compound in saline water. In contrast, various bacterial groups dominated the assimilation of EPS. Actinobacteria and Betaproteobacteria contributed the most in the freshwater section, whereas Cytophaga-like bacteria and Alpha- and Gammaproteobacteria participated in EPS assimilation in the lower part of the estuary. In addition, we examined the fraction of bacteria in each group that assimilated glucose or EPS. Overall, the fraction of bacteria in all groups that assimilated glucose was higher than the fraction that assimilated EPS (15 to 30% versus 5 to 20%, respectively). We found no correlation between the relative abundance of a group in the estuary and the fraction of bacteria actively assimilating glucose or EPS; the more active groups were often less abundant. Our results imply that the bacterial community in the Delaware Estuary is not controlled solely by "bottom-up" factors such as dissolved organic matter.
Collapse
|
38
|
Biomass production and assimilation of dissolved organic matter by SAR11 bacteria in the Northwest Atlantic Ocean. Appl Environ Microbiol 2005; 71:2979-86. [PMID: 15932993 PMCID: PMC1151852 DOI: 10.1128/aem.71.6.2979-2986.2005] [Citation(s) in RCA: 121] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2004] [Accepted: 12/21/2004] [Indexed: 11/20/2022] Open
Abstract
Members of the SAR11 clade often dominate the composition of marine microbial communities, yet their contribution to biomass production and the flux of dissolved organic matter (DOM) is unclear. In addition, little is known about the specific components of the DOM pool utilized by SAR11 bacteria. To better understand the role of SAR11 bacteria in the flux of DOM, we examined the assimilation of leucine (a measure of biomass production), as well as free amino acids, protein, and glucose, by SAR11 bacteria in the Northwest Atlantic Ocean. We found that when SAR11 bacteria were >25% of total prokaryotes, they accounted for about 30 to 50% of leucine incorporation, suggesting that SAR11 bacteria were major contributors to bacterial biomass production and the DOM flux. Specific growth rates of SAR11 bacteria either equaled or exceeded growth rates for the total prokaryotic community. In addition, SAR11 bacteria were typically responsible for a greater portion of amino acid assimilation (34 to 61%) and glucose assimilation (45 to 57%) than of protein assimilation (< or = 34%). These data suggest that SAR11 bacteria do not utilize various components of the DOM pool equally and may be more important to the flux of low-molecular-weight monomers than to that of high-molecular-weight polymers.
Collapse
|
39
|
Contribution of SAR11 bacteria to dissolved dimethylsulfoniopropionate and amino acid uptake in the North Atlantic ocean. Appl Environ Microbiol 2004; 70:4129-35. [PMID: 15240292 PMCID: PMC444831 DOI: 10.1128/aem.70.7.4129-4135.2004] [Citation(s) in RCA: 126] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
SAR11 bacteria are abundant in marine environments, often accounting for 35% of total prokaryotes in the surface ocean, but little is known about their involvement in marine biogeochemical cycles. Previous studies reported that SAR11 bacteria are very small and potentially have few ribosomes, indicating that SAR11 bacteria could have low metabolic activities and could play a smaller role in the flux of dissolved organic matter than suggested by their abundance. To determine the ecological activity of SAR11 bacteria, we used a combination of microautoradiography and fluorescence in situ hybridization (Micro-FISH) to measure assimilation of (3)H-amino acids and [(35)S]dimethylsulfoniopropionate (DMSP) by SAR11 bacteria in the coastal North Atlantic Ocean and the Sargasso Sea. We found that SAR11 bacteria were often abundant in surface waters, accounting for 25% of all prokaryotes on average. SAR11 bacteria were typically as large as, if not larger than, other prokaryotes. Additionally, more than half of SAR11 bacteria assimilated dissolved amino acids and DMSP, whereas about 40% of other prokaryotes assimilated these compounds. Due to their high abundance and activity, SAR11 bacteria were responsible for about 50% of amino acid assimilation and 30% of DMSP assimilation in surface waters. The contribution of SAR11 bacteria to amino acid assimilation was greater than would be expected based on their overall abundance, implying that SAR11 bacteria outcompete other prokaryotes for these labile compounds. These data suggest that SAR11 bacteria are highly active and play a significant role in C, N, and S cycling in the ocean.
Collapse
|
40
|
Contribution of SAR11 bacteria to dissolved dimethylsulfoniopropionate and amino acid uptake in the North Atlantic ocean. Appl Environ Microbiol 2004. [PMID: 15240292 DOI: 10.1128/aem.70.7.4129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023] Open
Abstract
SAR11 bacteria are abundant in marine environments, often accounting for 35% of total prokaryotes in the surface ocean, but little is known about their involvement in marine biogeochemical cycles. Previous studies reported that SAR11 bacteria are very small and potentially have few ribosomes, indicating that SAR11 bacteria could have low metabolic activities and could play a smaller role in the flux of dissolved organic matter than suggested by their abundance. To determine the ecological activity of SAR11 bacteria, we used a combination of microautoradiography and fluorescence in situ hybridization (Micro-FISH) to measure assimilation of (3)H-amino acids and [(35)S]dimethylsulfoniopropionate (DMSP) by SAR11 bacteria in the coastal North Atlantic Ocean and the Sargasso Sea. We found that SAR11 bacteria were often abundant in surface waters, accounting for 25% of all prokaryotes on average. SAR11 bacteria were typically as large as, if not larger than, other prokaryotes. Additionally, more than half of SAR11 bacteria assimilated dissolved amino acids and DMSP, whereas about 40% of other prokaryotes assimilated these compounds. Due to their high abundance and activity, SAR11 bacteria were responsible for about 50% of amino acid assimilation and 30% of DMSP assimilation in surface waters. The contribution of SAR11 bacteria to amino acid assimilation was greater than would be expected based on their overall abundance, implying that SAR11 bacteria outcompete other prokaryotes for these labile compounds. These data suggest that SAR11 bacteria are highly active and play a significant role in C, N, and S cycling in the ocean.
Collapse
|