Rare Species Shift the Structure of Bacterial Communities Across Sphagnum Compartments in a Subalpine Peatland

Sphagnum-associated microbiomes are crucial to Sphagnum growth and peatland ecological functions. However, roles of rare species in bacterial communities across Sphagnum compartments are poorly understood. Here the structures of rare taxa (RT) and conditionally abundant and rare taxa (CART) from Sphagnum palustre peat (SP), S. palustre ectosphere (Ecto) and S. palustre endosphere (Endo) were investigated in the Dajiuhu Peatland, central China. Our results showed that plant compartment effects significantly altered the diversities and structures of bacterial communities. The Observed species and Simpson indices of RT and CART in alpha diversity significantly increased from Endo to SP, with those of Ecto in-between. The variations of community dissimilarities of RT and CART among compartments were consistent with those of whole bacterial communities (WBC). Network analysis indicated a non-random co-occurrence pattern of WBC and all keystone species are affiliated with RT and CART, indicating their important role in sustaining the WBC. Furthermore, the community structures of RT and CART in SP were significantly shaped by water table and total nitrogen content, which coincided with the correlations between WBC and environmental factors. Collectively, our results for the first time confirm the importance of rare species to bacterial communities through structural and predicted functional analyses, which expands our understanding of rare species in Sphagnum-associated microbial communities in subalpine peatlands.

Introduction

Peatlands are the largest terrestrial carbon pool and important nitrogen reservoir, covering an area of 4 × 10⁶ km² on the Earth (Limpens et al., 2008). Globally, northern peatlands store 473 ∼ 621 Gt C and 8 ∼ 15 Gt N, which approximately account for 33% of the soil organic C and 9 ∼ 16% of the terrestrial N (Limpens et al., 2006; Yu et al., 2010), respectively. These abundant organic matters mainly originate from Sphagnum mosses compiling about 300 species around the world (Daniels, 1990; Shaw et al., 2003). Sphagnum mosses are the prevalent vegetation in acidic peatlands and greatly contribute to peatland development. Sphagnum mosses host diverse bacterial communities either endosymbiotic or ectosymbiotic with the plant (Opelt et al., 2007a, b; Bragina et al., 2012, 2014; Tian et al., 2019). The bacterial communities can produce antifungal compounds to improve pathogen defense of hosts and provide nitrogen for Sphagnum mosses via nitrogen-fixing bacteria such as Burkholderia (Opelt et al., 2007c; Shcherbakov et al., 2013, 2015). Besides the above-mentioned functions, endosymbiotic methanotrophs are demonstrated to provide carbon dioxide for the photosynthesis of Sphagnum mosses via methane oxidation (Raghoebarsing et al., 2005; Larmola et al., 2010; Kip et al., 2011). Therefore these Sphagnum-associated microorganisms play important roles in Sphagnum growth and elemental cycling in peatland ecosystems (Bragina et al., 2014; Weston et al., 2015; Kostka et al., 2016).

Generally, microbial communities in peatlands are comprised of a few dominant taxa (relative abundance >1%) with a long tail of low abundance species (Gilbert and Mitchell, 2006; Andersen et al., 2013; Tian et al., 2019), as well as in other habitats (Delgado-Baquerizo et al., 2018). These abundant taxa largely contribute to carbon cycling and energy flow due to their broad niche, high competitive ability and fast growth rate (Campbell et al., 2011). However, little knowledge is known about the low abundance populations, designated as rare biosphere (Sogin et al., 2006), or opportunistic species (Jousset et al., 2017), with high diversity. They may be part of the microbial “seed bank” and contribute disproportionately to community diversities and variations (Pedrós-Alió, 2007; Shade et al., 2014). They may become dominant and conduct essential functions in nutrient cycling, which compensate for the function deficiency of abundant species (Jousset et al., 2017) under favorable environmental conditions (Pedrós-Alió, 2007; Aanderud et al., 2015). Additionally, the responses of abundant and rare communities to environmental factors are different (Chen et al., 2017; Ruiz-González et al., 2019) due to the differences in niche breadth, competitiveness of resources and growth rate (Wu et al., 2017; Xue et al., 2018). Previous studies about Sphagnum-associated microbial communities mainly concentrate on abundant taxa, whereas the distribution of rare species among different compartments of Sphagnum (surface peat, ectosphere and endosphere of Sphagnum) still remains largely unknown. Therefore, we hypothesize that (1) rare species are major contributors to bacterial community variations across different compartments of Sphagnum palustre; (2) rare species may play an important role in the whole co-occurrence network and potential functions of bacterial communities; (3) rare species and abundant taxa may response differently to environmental factors in surface peat.

To test our hypotheses, bacterial communities were investigated via 16S rRNA gene sequencing from different compartments of S. palustre (S. palustre peat, S. palustre ectosphere, and S. palustre endosphere) in the Dajiuhu Peatland, Central China. In present study, we aim to address the following questions: (1) compare the diversity and composition of abundant and rare species, and their relative contributions to the variation of whole bacterial communities; (2) reveal the co-occurrence patterns and potential functional variations of abundant and rare species in Sphagnum compartments; (3) explore the response in abundant and rare species to environmental factors in S. palustre peat. Altogether, this study will enhance our understanding of Sphagnum-associated bacterial communities in different compartments, with emphasis on rare species.

Materials and Methods

Depiction of Site and Sampling

Dajiuhu (31°24′ – 31°33′ N, 109°56′ – 110°11′ E) is a typical subalpine peatland located in Shennongjia, Hubei province, central China (Supplementary Figure S1A). Peat has developed in a closed intermontane basin (1,730 m, above sea level), with an area around 16 km² and a depth of 2–3 m (Huang et al., 2018). The mean annual temperature and precipitation are 7.2°C and 1,560 mm, respectively. The present vegetation types are mainly Carex spp., S. palustre, Sanguisorba officinalis, and Euphorbia esula accompanying shrubs. More information about the Dajiuhu Peatland was given by Huang et al. (2012) and Li et al. (2016).

In this study, S. palustre and underlying peat samples (0–5 cm) were collected with sterile plastic bags and centrifuge tubes, respectively. Samples were collected in triplicate, in July 2016. Sampling sites were located at the peatland core area are Erhaoba (EHB), Niangniangfen (NNF) and two sites at Yangluchang (YLC) where S. palustre was actively growing (Supplementary Figure S1A). Water table at the four sites was determined in situ at the time of sampling. Visible plant root, litter and debris were removed from S. palustre peat samples. All the samples were transported to the geomicrobiology laboratory at China University of Geosciences (Wuhan) with ice box within 12 h.

Sample Processing

To get the microbes from the ectosphere, S. palustre was rinsed three times for 5 min each time with sterile deionized water to remove adherent debris. Then approx. 5 g of S. palustre samples were placed in 20 mL sterile phosphate buffered saline (PBS), and shaken and sonicated to separate the tightly adhering microbes (Bulgarelli et al., 2015; Edwards et al., 2015). The S. palustre plant material was transferred into another tube, and PBS solution was centrifuged at 10,000 g for 5 min and the supernatant was discarded (Supplementary Figure S1B). The remaining PBS fraction was stored at −80°C until DNA extraction.

The retrieved S. palustre designated for endosphere underwent further sonication and shake twice following the above-depicted approach. The S. palustre samples were then placed into 75% ethanol and 20% hydrogen peroxide for 1 min, respectively (Bulgari et al., 2009). The S. palustre samples were rinsed three times for 5 min each time with sterile deionized water and the final sterile water was inoculated on the R2A medium. The S. palustre plant materials were stored at −80°C until DNA extraction.

Additional S. palustre peat samples (approx. 20 g) were dried at 50°C for 24 h and subject to the organic matter and total nitrogen analyses.

DNA Extraction, Amplification, and Sequencing

Total DNA for each approximately 0.5 g freeze-dried S. palustre peat, ectosphere and endosphere fractions were extracted by PowerSoil DNA Isolation Kit (QIAGEN, Düsseldorf, Germany). The concentration and quality of DNA were checked with NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States) and 1% agarose gel stained with ethidium bromide. The V3-V4 regions of bacterial 16S rRNA gene were amplified with the primer set of 347F (5′-CCTACGGRRBGCASCAGKVRVGAAT-3′) and 802R (5′-GGACTACNVGGGTWTCTAATCC-3′) (Ren et al., 2019). PCR reactions contained 25 μL 2× Premix Taq (Takara Biotechnology, Dalian Co., Ltd., China), 1 μL each primer (10 mM) and 3 μL DNA (20 ng/μL) template in a volume of 50 μL. PCR products were then purified and quantified with Quit 2.0 fluorometer (Invitrogen, Carlsbad, CA, United States). DNA libraries were generated using NEBNext Ultra DNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, United States) following manufacturer’s protocols and the quality of libraries was validated by the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, United States) and the Qubit 2.0 Fluorometer. DNA libraries were sequenced on an Illumina MiSeq platform with PE300 paired-end reads at GENEWIZ, Inc. (Suzhou, China).

Sequence Analysis

Raw sequences were processed in Quantitative Insight Into Microbial Ecology (QIIME v. 1.9.1) platform (Caporaso et al., 2010). The sequences were assigned to each sample based on their unique barcode, followed by the removal of primers and barcodes. The forward and reverse sequences were paired-end assembly with overlaps more than 20 bp. The high quality reads were determined according to without any ambiguous base “N”, length >200 bp and base quality score >20. The clean sequences were clustered into operational taxonomic units (OTUs) by UCLUST (Edgar, 2010) based on 3% sequence dissimilarity and singletons (OTU with only one sequence in all samples) were eliminated. The mitochondria and chloroplast were removed and the chimera sequences were identified by UCHIME algorithm (Edgar et al., 2011). Represent sequences for each OTU were classified using Ribosomal Database Project classifier at a threshold of 0.8 (Cole et al., 2009) and annotated in the SILVA 119 database. All the samples were resampled to 29,577 sequences per sample by mothur (Schloss et al., 2009).

Raw 16S rRNA gene data was available at the NCBI Sequence Read Archive database¹ with the accession number of PRJNA512496.

Definition of OTU Classification

We set relative abundance thresholds as 0.01% for rare taxa (RT) and 1% for abundant taxa (AT) and classified all OTUs (archaea and unclassified OTUs excluded) into six subcategories (Dai et al., 2016; Xue et al., 2018) to compare the roles of taxa with different abundance in the whole bacterial communities (WBC): always abundant taxa (AAT) with a relative abundance >1% in all samples; conditionally abundant taxa (CAT) with a relative abundance ≥1% in some samples and always ≥0.01% in all samples; always rare taxa (ART) with a relative abundance <0.01% in all samples; conditionally rare taxa (CRT) with a relative abundance <0.01% in some samples but never abundant (≥1%) in any samples; moderate taxa (MT) with a relative abundance between 0.01% and 1% in all samples; conditionally abundant and rare taxa (CART) with a relative abundance varying from rare (<0.01%) to abundant (≥1%). The detailed depictions of OTU category were listed in Supplementary Table S1.

Date Processing and Analysis

The contribution of each subcategory to community variations was determined by similarity percentage (SIMPER) analysis in PAST 3.26². The alpha diversity indices (Observed species, Chao1, ACE, Simpson, Pielou, and Shannon), NMDS, PERMANOVA, and ANOSIM analyses were performed using “vegan” package in R 3.5.1³. Functional profiles were predicted using Tax4Fun package (Aßhauer et al., 2015) and principal component analysis (PCA) was employed to reveal compartment differences of potential functional communities based on Euclidean distance of the relative abundance of overall Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthologs (KOs) with “ggord” package in R 3.5.1. Differences in alpha diversity and potential function were tested with one-way analysis of variance (ANOVA) in SPSS 18⁴. Beta diversity analysis was performed based on Bray-Curtis distance of the relative abundance of OTUs in non-metric multidimensional scaling (NMDS). Analysis of similarity (ANOSIM) was used to test the differences among groups. Monte–Carlo tests of redundancy analysis determined the relationships between environmental factors and bacterial communities of S. palustre peat in Canoco 5⁵.

The OTUs with relative abundance above 0.01% were selected, and the correlation matrix between two OTUs was calculated in “psych” R package. Correlation coefficient (r > 0.8 or r < −0.8) and significance (FDR-corrected P < 0.01) were integrated into network analysis (Jiao et al., 2016). Topology attributes of the network (numbers of node and edge, average degree, network diameter, graph density, modularity, clustering coefficient and average of path length) were calculated using “igraph” package in R 3.5.1. Network visualization was conducted with Gephi 0.9.2⁶. Meanwhile, 1,000 Erdös-Réyni random networks were constructed to compare with the real network (Erdös and Réyni, 1960). Nodes with high degree (>100) and low betweenness centrality (<5,000) were considered as the keystone species in co-occurrence networks (Ma et al., 2016).

Results

Diversity of Bacterial Communities

In total, 2,413,649 high quality sequences and 3,144 operational taxonomic units (OTUs) were obtained. After the removal of archaea and unclassified OTUs, 2,848 OTUs were subject to further analysis. Among these, 13 OTUs representing 9.01% of all sequences were abundant taxa (AT) and present in all samples, whereas 2,706 OTUs occupying 39.2% of all sequences belonged to rare taxa (RT). Conditionally abundant and rare taxa (CART) included 126 OTUs representing 51.13% of all sequences (Supplementary Table S1).

Significant differences were observed in alpha diversities of whole bacterial communities (WBC) and RT among compartments (Supplementary Table S2). Observed species, Simpson, Pielou, and Shannon indices of WBC and CART significantly decreased from Sphagnum palustre peat (SP) to S. palustre endosphere (Endo), with those of S. palustre ectosphere (Ecto) in-between (Figure 1). In contrast, there were no significant differences in alpha diversity across sampling sites (Supplementary Table S2). Venn diagrams indicated over 50% OTUs in WBC were shared in three compartments, and the distributions of shared and specific OTUs in RT were similar to those of WBC (Figure 2B).

FIGURE 1

Figure 1. Changes in alpha diversity (Observed species, Simpson, Pielou, and Shannon indices) across Sphagnum palustre compartments. One-way analysis of variance (ANOVA) is used to compare the differences among S. palustre compartments. Cross sign represents the mean value (n = 12). WBC, whole bacterial communities; AT, abundant taxa; RT, rare taxa; CART, conditionally abundant and rare taxa; SP, S. palustre peat; Ecto, S. palustre ectosphere; Endo, S. palustre endosphere. The F and P-values of Observed species cannot be calculated in AT. Colored dots represent SP, Ecto, and Endo samples, respectively. Lines are fitted curves.

FIGURE 2

Figure 2. Structural variations of bacterial communities among S. palustre compartments. (A) Structure of bacterial communities based on Bray-Curtis distance in NMDS. (B) The number of shared and unique OTUs across different compartments of S. palustre. (C) The community dissimilarities of bacterial communities between S. palustre compartments. Each box indicates the maximum, minimum, 75th and 25th quartiles, respectively, and line within each box represents the median (n = 144). Lowercase letters indicate significances at the 95% confidence interval according to ANOVA with Tukey’s HSD post hoc test. The abbreviations are the same as indicated in Figure 1.

The dissimilarities in bacterial communities across compartments were visualized by NMDS analysis based on Bray-Curtis matrix. The communities of the same compartment clustered together (Figure 2A). The pairwise comparison of between group distances (Figure 2C), and the results of ANOSIM (Table 1) indicated that bacterial communities were significantly different among compartments, rather than among sampling sites.

TABLE 1

Table 1. Dissimilarity test showing differences of bacterial communities across compartments and sampling sites.

Composition of Bacterial Communities

As for WBC, Acidobacteria and Alphaproteobacteria dominated SP, whereas the predominant phyla in Ecto were Cyanobacteria, Alphaproteobacteria, and Gammaproteobacteria, followed by Acidobacteria. In Endo, Cyanobacteria and Alphaproteobacteria dominated the communities, followed by Acidobacteria (Figure 3A). Alphaproteobacteria was commonly found across all compartments and Acidobacteria presented frequently in SP, whereas Cyanobacteria was abundant in Ecto and Endo.

FIGURE 3

Figure 3. Composition of bacterial communities in WBC (A), AT (B), RT (C), and CART (D). E1, the first site of Erhaoba; N1, the first site of Niangniangfen; Y2, the second site of Yangluchang; Y6, the sixth site of Yangluchang. Others abbreviations are described in Figure 1.

In total, 13 OTUs were classified into AT, which were affiliated to Actinobacteria, Alphaproteobacteria and Gammaproteobacteria phyla (Figure 3B), and accounted for 0.46% of the OTU richness (Supplementary Table S3). As many as 2,706 OTUs were identified as RT, which were affiliated with 30 phyla (Figure 3C) and represented 95.01% of the OTU richness (Supplementary Table S3). Furthermore, 126 OTUs were recognized as CART belonging to 13 phyla (Figure 3D), and represented 4.42% of the OTU richness (Supplementary Table S3).

Co-occurrence Network

A co-occurrence network was constructed based on correlation relationships among OTU relative abundances. The resulting network contained 674 nodes linked by 8,455 edges, where positive correlations were commonly observed compared to negative ones (Figure 4C). The scale-free degree distribution of the real network (power-law: R² = 0.945, P < 0.0001, Supplementary Figure S2) indicated a scale-free and non-random network structure. Meanwhile the network had small-world features and modular structure as indicated by higher values of the modularity index (0.762) and clustering coefficient (0.521) compared with those of identical size Erdös-Réyni random networks (Erdös and Réyni, 1960; Table 2).

FIGURE 4

Figure 4. Co-occurrence network of bacterial communities sorted in color by phylum (A), modularity (B) and OTU category (C) rank, respectively. A connection between two nodes stands for a significant correlation (r > 0.8 or r < –0.8, P < 0.01). The numbers inside parentheses indicate the number of nodes in (A,B,C). Each node represents a single OTU. Positive edge and negative edge numbers among taxa with various abundances are put outside and inside the parentheses in a small panel at the bottom right corner. The size of each node is proportional to its degree; the thickness of edge linking two nodes is proportional to the corresponding correlation coefficient. RT, rare taxa; CART, conditionally abundant and rare taxa; AT, abundant taxa; MT, moderate taxa.

TABLE 2

Table 2. Topological properties of co-occurrence networks of microbial communities in the Dajiuhu Peatland.

The top 10 phyla in the network accounted for 93% of all nodes (Figure 4A). The whole network was grouped into 4 major modules (Figure 4B). Module I included 36.35% of the nodes in the entire network, which mainly fell into Proteobacteria. Module II and module III were mainly comprised of Proteobacteria and Acidobacteria occupying 41.99% of the nodes in the whole network, whereas module IV mainly consisted of Proteobacteria and Bacteroidetes (Supplementary Figure S3A). Ternary plots showed the specificity of different modules to a particular compartment (Supplementary Figures S3B–E). For example, most OTUs in module I showed high abundances in Ecto and Endo, whereas the majority of OTUs in modules II, III, and IV showed high proportions in SP.

Node-level topological features of microbial groups with different abundances (RT, CART, AT, and MT) indicated that the degree, betweenness and eigenvector centrality values of CART were significantly higher than that of RT (Supplementary Figures S4A,B,D). The closeness centrality values of AT were the highest among 4 microbial groups (Supplementary Figure S4C). Based on high degree (>100) and low betweenness centrality (<5,000) in the co-occurrence networks, 5 OTUs were identified as keystone species including Acidobacteria (1 OTU), Betaproteobacteria (1 OTU), Alphaproteobacteria (1 OTU), Ignavibacteriae (1 OTU) and Spirochaetae (1 OTU), which belonged to RT and CART (Supplementary Table S4).

Environmental Factors Related to Bacterial Communities

Monte–Carlo tests indicated that three physicochemical parameters, i.e., water table, total nitrogen and organic matter contents were significantly correlated with the variations of WBC, RT, CART, and AT in SP samples (Table 3). The WBC was significantly affected by water table and total nitrogen content (P < 0.05), which were also identified as the significant determinants (P < 0.05) in RT and CART. In contrast, organic matter content rather than water table and total nitrogen significantly (P < 0.05) shaped the structure of AT.

TABLE 3

Table 3. Correlations between bacterial communities and environmental factors in S. palustre peat samples in the Dajiuhu Peatland.

Potential Functions of Bacterial Communities

Similar to the structural variations of communities, the potential functions of communities showed an obvious separation among WBC, RT, and CART (Supplementary Figures S5A,C,D) compared to that in AT across compartments (Supplementary Figure S5B) as indicated by principal component analysis (PCA) based on Euclidean distance. The dissimilarity comparisons in each two compartments and sampling sites were further quantified using PERMANOVA (Supplementary Table S5), which showed significant (α = 0.05) differences in functional potentials across compartments rather than those among sampling sites.

Among WBC, RT, and CART, the relative abundance of potential nitrogen-fixing functional groups was the highest in SP, followed by those in Endo, and the lowest was in Ecto (Figure 5A). The relative abundance of predicted methane-oxidizing functional groups significantly decreased from Endo to SP, with that of Ecto in-between (Figure 5B). Together, specific potential functional groups such as N₂-fixers and CH₄-oxidizers in RT and CART can reflect consistent variation patterns with those in WBC across compartments, whereas AT did not.

FIGURE 5

Figure 5. The relative abundances of predicted nitrogen-fixing (A) and methane-oxidizing groups (B) in WBC, AT, RT, and CART. Abbreviations are the same as those in Figure 1. Significant differences (ANOVA tests) are labeled as follows: ^∗∗P < 0.01; ^∗∗∗P < 0.001. Lowercase letters represent significant differences at 95% confidence interval.

Discussion

The main objective of this study was to explore the roles of taxa with various abundances on the shift of bacterial communities in Sphagnum compartments at a subalpine peatland. Our results revealed that rare species (referring to RT and CART and hereafter) were the major contributors to the shift of Sphagnum-associated bacterial communities in different compartments. Rare species performed a consistent variation pattern with WBC in structures, predicted functional profiles and response to environmental factors. Additionally, they exerted large effects on the WBC co-occurrence network. Thus, our study provides new insights about the roles of rare species on Sphagnum-associated microbial communities.

Changes in Bacterial Communities

Our results validated that compartment effects drove the shift of bacterial communities in S. palustre as indicated by alpha and beta diversities. Alpha diversity decreased from SP to Ecto and was the lowest in Endo (Figure 1). This result agrees with the view of plant-associated microbial colonization, which shows a reduction in the species diversity from rhizosphere to rhizoplane to endosphere (Edwards et al., 2015; Beckers et al., 2017). The changes of alpha diversity from SP to Ecto to Endo may closely relate to the specific tissues of Sphagnum comprised of chlorocytes and hyalocytes. Microbial colonization in Sphagnum is strictly selected by chlorocytes (Opelt et al., 2007b; Bragina et al., 2012), and reduced the bacterial diversity in Sphagnum compartments (Figure 1). The significant increase of community dissimilarities (beta diversity) between compartments from Ecto vs. Endo to SP vs. Ecto and further to SP vs. Endo (Figure 2C and Table 1) also confirmed plant compartment effects on bacterial communities. Each compartment is a unique and heterogeneous microenvironment (Coleman-Derr et al., 2016; Stacie et al., 2016), which selects their specific microbial communities (Beckers et al., 2017). The plant ectosphere is usually oligotrophic with frequent fluctuations in temperature, precipitation, moisture and radiation, which is challenging for microbial survival (Müller et al., 2016). Therefore microorganisms living on the surface of plants have to deal with hard conditions, which resulted in a dispersive distribution of bacterial communities in Ecto (Figure 2A). In contrast, plant cells can act as oases for microbial survival in peatlands. Usually, microorganisms with specific traits, e.g., chemotaxis, motility, and quorum sensing are recruited by the host-dependent selection or passive transport, and thus lead to the formation of a set of specific plant-associated microbial communities (Compant et al., 2010; Kardol and Wardle, 2010; Bulgarelli et al., 2012; Philippot et al., 2013; Müller et al., 2016). Unlike vascular plants, rootless Sphagnum mosses intercept and hold water and nutrients efficiently through non-photosynthetic cells and adjacent photosynthetic cells (Weston et al., 2015; Kostka et al., 2016), which can facilitate the capture of their endosymbiotic microorganisms.

Our results for the first time indicated that microbial taxa with various abundances accounted for different variations of microbial communities in a subalpine peatland. The richness of RT was more than 60 times higher than that of AT (Figure 1), suggesting their important contributions to microbial diversity (Lynch and Neufeld, 2015; Zhan and Macisaac, 2015). Variation patterns of microbial taxa with different abundances in alpha and beta diversities among S. palustre compartments (Figures 1, 2) supported the hypothesis 1 that rare species mainly contributed to bacterial community variations compared to abundant taxa.

Co-occurrence Pattern of Communities in Sphagnum Compartments

The results of co-occurrence network analysis demonstrated keystone species were crucial to the whole bacterial communities and might exert large effects on other community components. Co-occurrence network can provide profound and unique insights into microbial interactions and community assembly (Barberán et al., 2012). Our results indicated that more positive interactions were observed between RT and CART (Figure 4C) than negative ones, which suggested a possible cooperation among them (Ju et al., 2014). Keystone species are characterized by high degree and low betweenness centrality based on scale-free feature (Barabási, 2009; Berry and Widder, 2014), and recognized as the initial components and affect other community components in networks. In our study, all the keystone species identified were affiliated to RT and CART, which confirmed that rare species might play a critical role in community co-occurrence network. Among keystone species, Candidatus Nitrotoga are cold-adapted nitrite-oxidizing bacteria, which involved in nitrogen cycling (Alawi et al., 2007). Spirochaeta are anaerobic bacteria, and some members are capable of cellulose-derivative degradation (glucose and cellobiose) and methane production in peat soils (Schmidt et al., 2015; Juottonen et al., 2017). These keystone species may benefit Sphagnum development under extreme conditions and mediate the recycle of macroelements (C, N) (Bragina et al., 2014) and sustain ecosystem network (Layeghifard et al., 2017). Therefore, rare species may be important in supporting peatland ecosystem functioning.

To be noted, rare species widely spread in all modules, which may affect the functioning of host-associated environments (Jousset et al., 2017). For example, nitrogen-fixing and methane-oxidizing microbes play essential roles in nitrogen supply for Sphagnum (Larmola et al., 2014; Vile et al., 2014) and decomposition of organic matter (Chen and Murrell, 2010). In the predicted functional communities, the variations in relative abundance of genes encoding nitrogenase and methane monooxygenase in RT and CART were consistent with those in WBC among compartments (Figure 5). Therefore, we surmised that rare species were important in the co-occurrence network and specific functions of bacterial communities, and the hypothesis 2 was supported. Although network analysis can give us some insight of microbial interactions, it only shows a statistical correlation among microbial groups and does not directly prove microbial interactions. Combination of microscopy-based or co-culture experiments with co-occurrence analysis will provide more direct evidences of microbial interactions.

Environmental Factors Shaping Bacterial Communities

Environmental conditions strongly impact the structure of microbial communities (Bahram et al., 2018). The effects of water table and total nitrogen content on bacterial communities in SP samples (Table 3) matches well with previous studies in different types of peatlands (Mishra et al., 2014; Urbanová and Bárta, 2016; Zhong et al., 2017). Interestingly we further found that RT and CART also responded to water table and total nitrogen content. In contrast, AT was sensitive to organic matter content. This indicated that microbial taxa with various abundances responded differently to environmental factors. It has been demonstrated that the content and form of organic matter affect the bioavailability and the quality of organic carbon, and further affect the priority of microbial metabolism (Stres et al., 2008). Our results further demonstrated that organic matter content mainly affect AT in WBC. Therefore, the hypothesis 3 that rare species and abundant taxa responding differently to environmental factors was supported.

Microbial groups with different abundances showed distinct environmental responses, which may result from the environmental filtering or dispersal limitation (Chen et al., 2017; Wu et al., 2017; Xue et al., 2018). Our results (Table 3) indicated RT and CART were more sensitive to environmental filtering than AT did. Usually RT and CART showed weak competition capability and low growth rate (Campbell et al., 2011; Logares et al., 2015), whereas AT can compete for more substrates and well adapt to an ecosystem via active growth and proliferation in a wide niche (Pedrós-Alió, 2006; Jousset et al., 2017; Xue et al., 2018). This suggests that the composition of microbial communities is strongly affected by environmentally induced species sorting (Roberto et al., 2018). Therefore, taxa with various abundances in bacterial communities from SP were impacted by unequal environmental filtering.

CONCLUSION

Our results demonstrated that the structures and predicted functions of rare species largely shifted across different compartments of S. palustre in accordance with those of whole bacterial communities. These results expanded our knowledge about the compartment effects on bacterial communities in the peatlands by emphasizing the importance of rare species. Co-occurrence network analysis revealed that rare species may play central roles in bacterial communities. Like whole bacterial communities, rare species were also significantly shaped by water table and total nitrogen content in SP, whereas abundant taxa correlated more closely with organic matter content. Collectively, our results emphasize the important impact of rare species on bacterial communities. Rare species mainly contribute to the variations both in structures and potential functions of Sphagnum-associated microbial communities among compartments in subalpine peatland ecosystems.

Data Availability Statement

The datasets generated for this study can be found in the NCBI Sequence Read Archive database (PRJNA512496).

Author Contributions

WT conducted the experiments, analyzed the data, and drafted the manuscript. XX helped the data analyses. LM, SE, and XQ gave suggestions about the manuscript. RW helped the sample collection. HW designed the research, provided the funding, and revised the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (41572325), the Fundamental Research Funds for the Central Universities, the China University of Geosciences (Wuhan) (CUGCJ1703, CUGQY1922), and the Natural Science Foundation of Hubei Province (2019CFB245).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Baiying Man and Dr. Yuan Yun for their help on experiment process. We also thank the editor HC-Q, and two reviewers for their insightful comments in improving the clarity of this article.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.03138/full#supplementary-material

Footnotes

1. ^ https://www.ncbi.nlm.nih.gov/sra
2. ^ http://folk.uio.no/ohammer/past/
3. ^ https://www.r-project.org/
4. ^ https://spss.en.softonic.com/
5. ^ http://canoco5.com/
6. ^ https://gephi.org/

References

Aanderud, Z. T., Jones, S. E., Fierer, N., and Lennon, J. T. (2015). Resuscitation of the rare biosphere contributes to pulses of ecosystem activity. Front. Microbiol. 6:24. doi: 10.3389/fmicb.2015.00024

PubMed Abstract | CrossRef Full Text | Google Scholar

Alawi, M., Lipski, A., Sanders, T., Pfeiffer, E. M., and Spieck, E. (2007). Cultivation of a novel cold-adapted nitrite oxidizing betaproteobacterium from the Siberian Arctic. ISME J. 1, 256–264. doi: 10.1038/ismej.2007.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersen, R., Chapman, S. J., and Artz, R. R. E. (2013). Microbial communities in natural and disturbed peatlands: a review. Soil Biol. Biochem. 57, 979–994. doi: 10.1016/j.soilbio.2012.10.003

CrossRef Full Text | Google Scholar

Aßhauer, K. P., Wemheuer, B., Daniel, R., and Meinicke, P. (2015). Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics 31, 2882–2884. doi: 10.1093/bioinformatics/btv287

PubMed Abstract | CrossRef Full Text | Google Scholar

Bahram, M., Hildebrand, F., Forslund, S. K., Anderson, J. L., Soudzilovskaia, N. A., Bodegom, P. M., et al. (2018). Structure and function of the global topsoil microbiome. Nature 560, 233–237. doi: 10.1038/s41586-018-0386-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Barabási, A.-L. (2009). Scale-free networks: a decade and beyond. Science 325, 412–413. doi: 10.1126/science.1173299

PubMed Abstract | CrossRef Full Text | Google Scholar

Barberán, A., Bates, S. T., Casamayor, E. O., and Fierer, N. (2012). Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 6, 343–351. doi: 10.1038/ismej.2013.236

PubMed Abstract | CrossRef Full Text | Google Scholar

Beckers, B., Beeck, M. O. D., Weyens, N., Boerjan, W., and Vangronsveld, J. (2017). Structural variability and niche differentiation in the rhizosphere and endosphere bacterial microbiome of field-grown poplar trees. Microbiome 5, 25–41. doi: 10.1186/s40168-017-0241-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Berry, D., and Widder, S. (2014). Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front. Microbiol. 5:219. doi: 10.3389/fmicb.2014.00219

PubMed Abstract | CrossRef Full Text | Google Scholar

Bragina, A., Berg, C., Cardinale, M., Shcherbakov, A., Chebotar, V., and Berg, G. (2012). Sphagnum mosses harbour highly specific bacterial diversity during their whole lifecycle. ISME J. 6, 802–813. doi: 10.1038/ismej.2011.151

PubMed Abstract | CrossRef Full Text | Google Scholar

Bragina, A., Oberaunerwappis, L., Zachow, C., Halwachs, B., Thallinger, G. G., Müller, H., et al. (2014). The Sphagnum microbiome supports bog ecosystem functioning under extreme conditions. Mol. Ecol. 23, 4498–4510. doi: 10.1111/mec.12885

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulgarelli, D., Garrido-Oter, R., Münch, P. C., Weiman, A., Dröge, J., Pan, Y., et al. (2015). Structure and function of the bacterial root microbiota in wild and domesticated barley. Cell Host Microbe 17, 392–403. doi: 10.1016/j.chom.2015.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulgarelli, D., Rott, M., Schlaeppi, K., Ver, L. V. T. E., Ahmadinejad, N., Assenza, F., et al. (2012). Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 488, 91–95. doi: 10.1038/nature11336

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulgari, D., Casati, P., Brusetti, L., Quaglino, F., Brasca, M., Daffonchio, D., et al. (2009). Endophytic bacterial diversity in grapevine (Vitis vinifera L.) leaves described by 16S rRNA gene sequence analysis and length heterogeneity-PCR. J. Microbiol. 47, 393–401. doi: 10.1007/s12275-009-0082-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, B. J., Liying, Y., Heidelberg, J. F., and Kirchman, D. L. (2011). Activity of abundant and rare bacteria in a coastal ocean. Proc. Natl. Acad. Sci. U.S.A. 108, 12776–12781. doi: 10.1073/pnas.1101405108

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., and Murrell, J. C. (2010). Geomicrobiology: methanotrophs in moss. Nat. Geosci. 3, 595–596. doi: 10.1038/ngeo952

CrossRef Full Text | Google Scholar

Chen, Y., Pan, Y., Yu, L., Yang, J., and Zhang, W. (2017). Patterns and processes in marine microeukaryotic community biogeography from Xiamen coastal waters and intertidal sediments, southeast China. Front. Microbiol. 8:1912. doi: 10.3389/fmicb.2017.01912

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., et al. (2009). The Ribosomal database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37, D141–D145. doi: 10.1093/nar/gkn879

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman-Derr, D., Desgarennes, D., Fonseca-Garcia, C., Gross, S., Clingenpeel, S., Woyke, T., et al. (2016). Plant compartment and biogeography affect microbiome composition in cultivated and native Agave species. New Phytol. 209, 798–811. doi: 10.1111/nph.13697

PubMed Abstract | CrossRef Full Text | Google Scholar

Compant, S., Clément, C., and Sessitsch, A. (2010). Plant growth-promoting bacteria in the rhizo- and endosphere of plants: their role, colonization, mechanisms involved and prospects for utilization. Soil Biol. Biochem. 42, 669–678. doi: 10.1016/j.soilbio.2009.11.024

CrossRef Full Text | Google Scholar

Dai, T., Zhang, Y., Tang, Y., Bai, Y., Tao, Y., Huang, B., et al. (2016). Identifying the key taxonomic categories that characterize microbial community diversity using full-scale classification: a case study of microbial community in the sediments of Hangzhou Bay. FEMS Microbiol. Ecol. 92:fiw150. doi: 10.1093/femsec/fiw150

PubMed Abstract | CrossRef Full Text | Google Scholar

Daniels, R. E. (1990). “Key to European species of Sphagnum,” in Handbook of European Sphagna, 2nd Edn, ed. A. Eddy, (London: HMSO), 32–42.

Google Scholar

Delgado-Baquerizo, M., Oliverio, A. M., Brewer, T. E., Benavent-González, A., Eldridge, D. J., Bardgett, R. D., et al. (2018). A global atlas of the dominant bacteria found in soil. Science 359, 320–325. doi: 10.1126/science.aap9516

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

PubMed Abstract | CrossRef Full Text | Google Scholar

Edwards, J., Johnson, C., Santos-Medellín, C., Lurie, E., Podishetty, N. K., Bhatnagar, S., et al. (2015). Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. U.S.A. 112, E911–E920. doi: 10.1073/pnas.1414592112

PubMed Abstract | CrossRef Full Text | Google Scholar

Erdös, P., and Réyni, A. (1960). On the evolution of random graphs. Am. Math. Soc. 286, 257–274.

Google Scholar

Gilbert, D., and Mitchell, E. A. D. (2006). “Microbial diversity in Sphagnum peatlands,” in Developments in Earth Surface Processes, eds I. P. Martini, A. Martínez Cortizas, and W. Chesworth, (Amsterdam: Elsevier), 287–318. doi: 10.1016/s0928-2025(06)09013-4

CrossRef Full Text | Google Scholar

Huang, X., Meyers, P., Yu, J., Wang, X., Huang, J., Jin, F., et al. (2012). Moisture conditions during the Younger Dryas and the early Holocene in the middle reaches of the Yangtze River, central China. Holocene 22, 1473–1479. doi: 10.1177/0959683612450202

CrossRef Full Text | Google Scholar

Huang, X., Pancost, R. D., Xue, J., Gu, Y., Evershed, R. P., and Xie, S. (2018). Response of carbon cycle to drier conditions in the mid-Holocene in central China. Nat. Commun. 9, 1369–1377. doi: 10.1038/s41467-018-03804-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiao, S., Liu, Z., Lin, Y., Yang, J., Chen, W., and Wei, G. (2016). Bacterial communities in oil contaminated soils: biogeography and co-occurrence patterns. Soil Biol. Biochem. 98, 64–73. doi: 10.1016/j.soilbio.2016.04.005

CrossRef Full Text | Google Scholar

Jousset, A., Bienhold, C., Chatzinotas, A., Gallien, L., Gobet, A., Kurm, V., et al. (2017). Where less may be more: how the rare biosphere pulls ecosystems strings. ISME J. 11, 853–862. doi: 10.1038/ismej.2016.174

PubMed Abstract | CrossRef Full Text | Google Scholar

Ju, F., Xia, Y., Guo, F., Wang, Z., and Zhang, T. (2014). Taxonomic relatedness shapes bacterial assembly in activated sludge of globally distributed wastewater treatment plants. Environ. Microbiol. 16, 2421–2432. doi: 10.1111/1462-2920.12355

PubMed Abstract | CrossRef Full Text | Google Scholar

Juottonen, H., Eiler, A., Biasi, C., Tuittila, E. S., Yrjälä, K., and Fritzec, H. (2017). Distinct anaerobic bacterial consumers of cellobiose-derived carbon in boreal fens with different CO₂ vs. CH₄ production ratio. Appl. Environ. Microb. 83, e02533–e02548. doi: 10.1128/AEM.02533-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Kardol, P., and Wardle, D. A. (2010). How understanding aboveground-belowground linkages can assist restoration ecology. Trends Ecol. Evol. 25, 670–679. doi: 10.1016/j.tree.2010.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Kip, N., Ouyang, W., van Winden, J., Raghoebarsing, A., van Niftrik, L., Pol, A., et al. (2011). Detection, isolation, and characterization of acidophilic methanotrophs from Sphagnum mosses. Appl. Environ. Microbiol. 77, 5643–5654. doi: 10.1128/AEM.05017-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Kostka, J. E., Weston, D. J., Glass, J. B., Lilleskov, E. A., Shaw, A. J., and Turetsky, M. R. (2016). The Sphagnum microbiome: new insights from an ancient plant lineage. New Phytol. 211, 57–64. doi: 10.1111/nph.13993

PubMed Abstract | CrossRef Full Text | Google Scholar

Larmola, T., Tuittila, E. S., Aarva, M., Fritze, H., and Tiirola, M. (2014). Methanotrophy induces nitrogen fixation during peatland development. Proc. Natl. Acad. Sci. U.S.A. 111, 734–739. doi: 10.1073/pnas.1314284111

PubMed Abstract | CrossRef Full Text | Google Scholar

Larmola, T., Tuittila, E.-S., Tiirola, M., Nykänen, H., Martikainen, P. J., Yrjälä, K., et al. (2010). The role of Sphagnum mosses in the methane cyclinag of a boreal mire. Ecology 91, 2356–2365. doi: 10.1890/09-1343.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Layeghifard, M., Hwang, D. M., and Guttman, D. S. (2017). Disentangling interactions in the microbiome: a network perspective. Trends Microbiol. 25, 217–228. doi: 10.1016/j.tim.2016.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Ma, C., Zhu, C., Huang, R., and Zheng, C. (2016). Historical anthropogenic contributions to mercury accumulation recorded by a peat core from Dajiuhu montane mire, central China. Environ. Pollut. 216, 332–339. doi: 10.1016/j.envpol.2016.05.083

PubMed Abstract | CrossRef Full Text | Google Scholar

Limpens, J., Berendse, F., Blodau, C., Canadel, J. G., Freeman, C., Holden, J., et al. (2008). Peatlands and the carbon cycle: from local processes to global implications - a synthesis. Biogeosciences 5, 1379–1419. doi: 10.5194/bg-5-1475-2008

CrossRef Full Text | Google Scholar

Limpens, J., Heijmans, M. M., and Berendse, F. (2006). “The nitrogen cycle in boreal peatlands,” in Boreal Peatland Ecosystems, eds R. K. Wieder, and D. H. Vitt, (Berlin: Springer), 195–230. doi: 10.1007/978-3-540-31913-9_10

CrossRef Full Text | Google Scholar

Logares, R., Mangot, J.-F. O., and Massana, R. (2015). Rarity in aquatic microbes: placing protists on the map. Res. Microbiol. 166, 831–841. doi: 10.1016/j.resmic.2015.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Lynch, M. D. J., and Neufeld, J. D. (2015). Ecology and exploration of the rare biosphere. Nat. Rev. Microbiol. 13, 217–229. doi: 10.1038/nrmicro3400

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, B., Wang, H., Dsouza, M., Lou, J., He, Y., Dai, Z., et al. (2016). Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. ISME J. 10, 1891–1901. doi: 10.1038/ismej.2015.261

PubMed Abstract | CrossRef Full Text | Google Scholar

Mishra, S., Lee, W. A., Hooijer, A., Reuben, S., Sudiana, I. M., Idris, A., et al. (2014). Microbial and metabolic profiling reveal strong influence of water table and land-use patterns on classification of degraded tropical peatlands. Biogeosciences 11, 14009–14042. doi: 10.5194/bg-11-1727-2014

CrossRef Full Text | Google Scholar

Müller, D. B., Vogel, C., Yang, B., and Vorholt, J. A. (2016). The plant microbiota: systems-level insights and perspectives. Annu. Rev. Genet. 50, 211–234. doi: 10.1146/annurev-genet-120215-034952

PubMed Abstract | CrossRef Full Text | Google Scholar

Opelt, K., Berg, C., and Berg, G. (2007a). The bryophyte genus Sphagnum is a reservoir for powerful and extraordinary antagonists and potentially facultative human pathogens. FEMS Microbiol. Ecol. 61, 38–53. doi: 10.1111/j.1574-6941.2007.00323.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Opelt, K., Berg, C., Schönmann, S., Eberl, L., and Berg, G. (2007b). High specificity but contrasting biodiversity of Sphagnum-associated bacterial and plant communities in bog ecosystems independent of the geographical region. ISME J. 1, 502–516. doi: 10.1038/ismej.2007.58

PubMed Abstract | CrossRef Full Text | Google Scholar

Opelt, K., Chobot, V., Hadacek, F., Schönmann, S., Eberl, L., and Berg, G. (2007c). Investigations of the structure and function of bacterial communities associated with Sphagnum mosses. Environ. Microbiol. 9:2795. doi: 10.1111/j.1462-2920.2007.01391.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedrós-Alió, C. (2006). Marine microbial diversity: can it be determined? Trends Microbiol. 14, 257–263. doi: 10.1016/j.tim.2006.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedrós-Alió, C. (2007). Ecology: dipping into the rare biosphere. Science 315, 192–193. doi: 10.1126/science.1135933

PubMed Abstract | CrossRef Full Text | Google Scholar

Philippot, L., Raaijmakers, J. M., Lemanceau, P., and Wh, V. D. P. (2013). Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799. doi: 10.1038/nrmicro3109

PubMed Abstract | CrossRef Full Text | Google Scholar

Raghoebarsing, A. A., Smolders, A. J., Schmid, M. C., Rijpstra, W. I., Woltersarts, M., Derksen, J., et al. (2005). Methanotrophic symbionts provide carbon for photosynthesis in peat bogs. Nature 436, 1153–1156. doi: 10.1038/nature03802

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, Z., Qu, X., Zhang, M., Yu, Y., and Peng, W. (2019). Distinct bacterial communities in wet and dry seasons during a seasonal water level fluctuation in the largest freshwater lake (Poyang Lake) in China. Front. Microbiol. 10:1167. doi: 10.3389/fmicb.2019.01167

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberto, A. A., Van Gray, J. B., and Leff, L. G. (2018). Sediment bacteria in an urban stream: spatiotemporal patterns in community composition. Water Res. 134, 353–369. doi: 10.1016/j.watres.2018.01.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruiz-González, C., Logares, R., Sebastián, M., Mestre, M., Rodríguez-Martínez, R., Galí Tàpias, M., et al. (2019). Higher contribution of globally rare bacterial taxa reflects environmental transitions across the surface ocean. Mol. Ecol. 28, 1930–1945. doi: 10.1111/mec.15026

PubMed Abstract | CrossRef Full Text | Google Scholar

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microb. 75, 7537–7541. doi: 10.1128/AEM.01541-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, O., Horn, M. A., Kolb, S., and Drake, H. L. (2015). Temperature impacts differentially on the methanogenic food web of cellulose-supplemented peatland soil. Environ. Microbiol. 17, 720–734. doi: 10.1111/1462-2920.12507

PubMed Abstract | CrossRef Full Text | Google Scholar

Shade, A., Jones, S. E., Caporaso, J. G., Handelsman, J., Knight, R., Fierer, N., et al. (2014). Conditionally rare taxa disproportionately contribute to temporal changes in microbial diversity. mBio 5, e01371–e01379. doi: 10.1128/mBio.01371-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaw, A. J., Cox, C. J., and Boles, S. B. (2003). Global patterns in peatmoss biodiversity. Mol. Ecol. 12, 2553–2570. doi: 10.1046/j.1365-294X.2003.01929.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shcherbakov, A. V., Bragina, A. V., Kuzmina, E. Y., Berg, C., Muntyan, A. N., Makarova, N. M., et al. (2013). Endophytic bacteria of Sphagnum mosses as promising objects of agricultural microbiology. Microbiology 82, 306–315. doi: 10.1134/S0026261713030107

CrossRef Full Text | Google Scholar

Shcherbakov, A. V., Kuzmina, E. Y., Lapshina, E. D., Shcherbakova, E. N., Gonchar, L. N., and Chebotar, V. K. (2015). Taxonomic diversity of bacterial populations inhabiting gametophytes of Sphagnum mosses from different geographic regions of Russia. Agron. Res. 13, 192–201.

Google Scholar

Sogin, M. L., Morrison, H. G., Huber, J. A., Welch, D. M., Huse, S. M., Neal, P. R., et al. (2006). Microbial diversity in the deep sea and the underexplored “rare biosphere”. Proc. Natl. Acad. Sci. U.S.A. 103, 12115–12120. doi: 10.1073/pnas.0605127103

PubMed Abstract | CrossRef Full Text | Google Scholar

Stacie, T., Yergeau, É, Tremblay, J., Legendre, P., Whyte, L. G., and Greer, C. W. (2016). The willow Microbiome is influenced by soil petroleum-hydrocarbon concentration with plant compartment-specific effects. Front. Microbiol. 7:1363. doi: 10.3389/fmicb.2016.01363

PubMed Abstract | CrossRef Full Text | Google Scholar

Stres, B., Danevcic, T., Pal, L., Fuka, M. M., Resman, L., Leskovec, S., et al. (2008). Influence of temperature and soil water content on bacterial, archaeal and denitrifying microbial communities in drained fen grassland soil microcosms. FEMS Microbiol. Ecol. 66, 110–122. doi: 10.1111/j.1574-6941.2008.00555.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, W., Wang, H., Xiang, X., Wang, R., and Xu, Y. (2019). Structural variations of bacterial community driven by Sphagnum microhabitat differentiation in a subalpine peatland. Front. Microbiol. 10:1661. doi: 10.3389/fmicb.2019.01661

PubMed Abstract | CrossRef Full Text | Google Scholar

Urbanová, Z., and Bárta, J. (2016). Effects of long-term drainage on microbial community composition vary between peatland types. Soil Biol. Biochem. 92, 16–26. doi: 10.1016/j.soilbio.2015.09.017

CrossRef Full Text | Google Scholar

Vile, M. A., Wieder, R. K., Živković, T., Scott, K. D., Vitt, D. H., Hartsock, J. A., et al. (2014). N₂-fixation by methanotrophs sustains carbon and nitrogen accumulation in pristine peatlands. Biogeochemistry 121, 317–328. doi: 10.1007/s10533-014-0019-6

CrossRef Full Text | Google Scholar

Weston, D. J., Timm, C. M., Walker, A. P., Gu, L., Muchero, W., Schmutz, J., et al. (2015). Sphagnum physiology in the context of changing climate: emergent influences of genomics, modelling and host-microbiome interactions on understanding ecosystem function. Plant Cell Environ. 38, 1737–1751. doi: 10.1111/pce.12458

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, W., Logares, R., Huang, B., and Hsieh, C. H. (2017). Abundant and rare picoeukaryotic sub-communities present contrasting patterns in the epipelagic waters of marginal seas in the northwestern Pacific Ocean. Environ. Microbiol. 19, 287–300. doi: 10.1111/1462-2920.13606

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, Y., Chen, H., Yang, J. R., Min, L., Huang, B., and Yang, J. (2018). Distinct patterns and processes of abundant and rare eukaryotic plankton communities following a reservoir cyanobacterial bloom. ISME J. 12, 2263–2277. doi: 10.1038/s41396-018-0159-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Z. C., Loisel, J., Brosseau, D. P., Beilman, D. W., and Hunt, S. J. (2010). Global peatland dynamics since the last glacial maximum. Geophys. Res. Lett. 37, 69–73. doi: 10.1029/2010GL043584

CrossRef Full Text | Google Scholar

Zhan, A., and Macisaac, H. J. (2015). Rare biosphere exploration using high-throughput sequencing: research progress and perspectives. Conserv. Genet. 16, 513–522. doi: 10.1007/s10592-014-0678-9

CrossRef Full Text | Google Scholar

Zhong, Q., Chen, H., Liu, L., He, Y., Zhu, D., Jiang, L., et al. (2017). Water table drawdown shapes the depth-dependent variations in prokaryotic diversity and structure in Zoige peatlands. FEMS Microbiol. Ecol. 93:fix049. doi: 10.1093/femsec/fix049

PubMed Abstract | CrossRef Full Text | Google Scholar