- Research article
- Open Access
Co-occurrence of planktonic bacteria and archaea affects their biogeographic patterns in China’s coastal wetlands
Environmental Microbiome volume 16, Article number: 19 (2021)
Planktonic bacteria and archaea play a key role in maintaining ecological functions in aquatic ecosystems; however, their biogeographic patterns and underlying mechanisms have not been well known in coastal wetlands including multiple types and at a large space scale. Therefore, planktonic bacteria and archaea and related environmental factors were investigated in twenty-one wetlands along China’s coast to understand the above concerns. The results indicated that planktonic bacteria had different biogeographic pattern from planktonic archaea, and both patterns were not dependent on the wetland's types. Deterministic selection shapes the former’s community structure, whereas stochastic processes regulate the latter’s, being consistent with the fact that planktonic archaea have a larger niche breadth than planktonic bacteria. Planktonic bacteria and archaea co-occur, and their co-occurrence rather than salinity is more important in shaping their community structure although salinity is found to be a main environmental deterministic factor in the coastal wetland waters. This study highlights the role of planktonic bacteria-archaea co-occurrence on their biogeographic patterns, and thus provides a new insight into studying underlying mechanisms of microbial biogeography in coastal wetlands.
Understanding microbial biogeographic patterns is a fundamental task of aquatic microbial ecology [1,2,3,4,5,6]. The biogeographic studies on aquatic microorganisms have experienced from arguing whether there are biogeographic patterns to concerning community assembly mechanisms [7,8,9]. Synchronously, the highlight is from a distance-decay relationship to a relative influence of deterministic and stochastic processes on present-day distribution patterns [10,11,12]. Deterministic processes are based on traditional niche theory including biotic interactions and environmental filtering, while stochastic processes are based on neutral theory including drift and extinction [13, 14]. Although four basic conceptual processes (i.e., selection, drift, dispersal, and mutation) are well known to create and maintain microbial biogeographic patterns [15,16,17], it is still a challenge to clarify their actual control mechanisms, especially with respect to biological interplay.
Coastal wetlands, located in the transition zone between the land and the sea, include multiple types such as estuary, delta, saltmarsh, mangrove, and seagrass wetland . They provide many valuable ecosystem services (e.g., flood protection, erosion control, and wildlife habitat) and are thus very important for human well-being . Planktonic bacteria and archaea play an important role in driving nutrient biogeochemical cycles and maintaining ecological functions in aquatic ecosystems [20,21,22]. However, to the best of our knowledge, their biogeographic patterns are largely unknown in coastal wetlands including multiple types and at a large space scale. The existing studies are confined to a certain estuary , saltmarsh , or mangrove, and the subjects are mainly focused on soil microorganisms. Therefore, it is necessary to understand biogeographic patterns of planktonic bacteria and archaea and the processes shaping their community assembly in the coastal wetlands.
Twenty-one China’s coastal wetlands were investigated to understand the biogeographic patterns and control mechanisms of their planktonic bacteria and archaea. The wetlands stretch across nearly 20° of latitude from the northernmost of Liaoning Province to the southernmost of Hainan Province, including reservoir, river, estuary, delta, lagoon, bay, marsh, and mangrove wetland. The inlet and outlet, open water area, and vegetation area of wetland are considered to understand microbial regional spatial difference. Our hypothesis is that planktonic bacteria and archaea co-occurrence to shape their biogeographic patterns in coastal wetlands. Using 16S rRNA amplicon sequencing and quantitative PCR (qPCR), the community composition and abundance of planktonic bacteria and archaea were investigated. In addition, the physical and chemical parameters of wetland waters were measured, and null model and correlation-based network are carried out to understand the underlying processes and the bacteria-archaea co-occurrence relationship. The main aims of this study are: (1) to understand planktonic bacteria and archaea community composition and controlling factors; (2) to discern the relative importance of deterministic and stochastic processes on the biogeographic patterns, and (3) to discover planktonic bacteria-archaea co-occurrence patterns in various coastal wetlands. Therefore, the study helps to improve our understanding of aquatic microbial ecology in coastal wetlands.
Study area and sampling
China’s coastal wetlands have subtropical and temperate monsoon climates. The average annual temperature increases gradually from north to south, but in summer the temperature difference is small. Twenty-one coastal wetlands were selected for this study (Fig. 1; Additional file 1: Table S1). Their types are diverse and include four reservoir, one lagoon, four estuary, five river, two marsh, one bay, one delta, and three mangrove wetlands. Except the mangrove wetlands, the main vegetation of wetlands are Phragmites communis and Acorus calamus. Samples were collected from the inlet, outlet, and wetland waters. For the wetlands, surface and bottom waters were collected in open water area, while only surface water was collected in vegetated area. As such, a total of 101 water samples were collected during June and August 2019. Water temperature (WT), dissolved oxygen (DO), chlorophyll a (Chl a), salinity, and pH were measured directly using professional plus multi-parameter probe (YSI, EXO1, America) with precorrection in situ. For analyses of nutrients, the water samples were filtered through 0.45 μm filters (GF/F-Whatman, UK) within 12 h after sampling, and the filtered waters were stored in 100 ml polyethylene bottles at 4 °C. Particulate organic carbon (POC) and nitrogen (PON) were collected with pre-combusted (450 °C, 8 h) 0.70 μm filters (GF/F-Whatman, UK). Bacterial and archaeal assemblages were collected by 0.22 μm sterilized filter membrane (MF-Millipore, USA) with diaphragm vacuum pump (GM-1.0A, China), and then stored at − 20 °C until DNA extraction.
Total nitrogen (TN), ammonium nitrogen (NH4-N), nitrate nitrogen (NO3-N), nitrite nitrogen (NO2-N), total phosphorus (TP), phosphate phosphorus (PO4-P), and dissolved silicon (DSi) were quantified by Skalar (San++, Netherlands) with detection limit of 0.01 mg L−1 for TN and NO3-N, 0.005 mg L−1 for NO2-N, 0.02 mg L−1 for NH4-N, TP, PO4-P, and DSi, respectively. POC and PON were measured by Elemental analyzer (Vario EL III, Germany). Dissolved organic carbon (DOC) was determined by a TOC analyzer (Aurora 1030, USA).
DNA extraction and quantitative PCR
Water microbial DNA extraction was carried out using the E.Z.N.ATM Water DNA Kit (OMEGA, USA) according to manufacturer’s protocols. The primer sets 906F/1062R  and 1106F/1378R  were respectively used for bacterial and archaeal 16S rRNA gene amplification. 16S rRNA gene abundance of planktonic bacteria and archaea in each sample was determined by qPCR using Bio-rad Fax (CFX ConnectTM, USA). Triplicate amplifications were conducted in a 20 μl reaction system containing 12.5 μl of Ultra SYBR Mixture (2 ×), 1.0 μl of each primer (10 nM), 10 ng of DNA template, and 4.5 μl ddH2O. The thermal cycling was performed in two steps: initial denaturation was at 95 °C for 10 min, 40 cycles of 95 °C for 15 s, 60 °C for 1 min, and then melting curve program was 95 °C for 15 s, 60 °C for 1 min, 95 °C for 15 s, and 60 °C for 15 s. The standard curve was generated by tenfold dilution of plasmids (containing single-copy bacteria or archaea 16S rRNA gene, respectively) from 1010 to 104 copies μl−1. The quantitative standard curve was built as bacteria 16S rRNA gene copies = − 2 × 1010 Ct + 5 × 1011 (R2 = 0.99) and archaea 16S rRNA gene copies = − 9 × 108 Ct + 3 × 1010 (R2 = 0.82), respectively, where Ct was the cycle threshold. 16S rRNA gene abundance was calculated based on the Ct value together with plasmid standard curve.
16S rRNA sequencing denoising, community composition and α-diversity analysis
Sequencing was performed on Illumina NovaSeq 6000 platform at Meige Technology Co., Ltd., Guangzhou, China. Raw fastq data were quality-filtered using Trimmomatic  to remove contaminating adaptors, low-quality ends of reads, and short length reads (< 200 bp). The overlapping paired-end reads with sequence mismatching < 5 bp and alignment similarity > 90% were merged using FLASH . Operational taxonomic units (OTUs) were clustered in UPARSE software at 97% consistency level . Singleton and doubleton OTUs representing sequencing errors were removed, and the representative OTUs with the highest occurrence frequency were assigned using Silva, RDP and Greengene . To equalize sequencing depth, each sample was rarefied to the minimum sequencing depth, and sequence normalization was performed using MOTHUR v.1.33.3 . Bacterial and archaeal raw data obtained have been deposited in NCBI SRA database with the accession numbers of PRJNA681135 and PRJNA674461. Phylogenetic analysis of community composition was performed in the R package, α-diversity including shannon_e was conducted with the “vegan” package v 2.0-4 in R .
Quantification of ecological processes using null model
Null model can integrate phylogenetic and taxonomic β-diversity metrics to quantify the relative importance of stochastic and deterministic processes on microbial community assembly [33, 34]. Beta mean nearest taxon metric (βMNTD) was calculated using the picante R package iCAMP , and then null model was implemented to calculate the beta nearest taxon index (βNTI) [33, 36]. |βNTI| > 2 is interpreted as variable selection of deterministic processes, and |βNTI| < 2 indicates that community composition difference is the result of stochastic processes . Further, the Bray–Curtis dissimilarities based Raup-Cricck metric (RCbary) was calculated. When |βNTI| < 2, RCbary < − 0.95 and RCbary > 0.95 indicate that community assembly is dominated by homogenizing dispersal and dispersal limitation respectively, and |RCbray| < 0.95 indicate the influence of the “non-dominant” fraction . Niche breadth of bacteria and archaea calculated using the “spaa” libraries in R.
Correlation-based network analysis has been extensively used for microbial interactions [38, 39], and the co-occurrence patterns illustrated by a network can capture important information in microbial ecology . Co-occurrence networks were constructed using the “WGCNA” libraries in R 4.0.0, and were visualized with “Gephi 0.9.2” software. To reduce complexity, only OTUs with the total number above 62 across all samples and occurring in more than 30% of all samples were retained, meanwhile, removing the relationship by themselves. The pairwise Spearman’s correlations between OTUs were calculated, with a correlation coefficient > 0.7 and a P value < 0.05 (Benjamini and Hochberg adjusted) being considered as a valid relationship. The network-level topological features of a network was calculated.
Spearman’s correlation analyses between physicochemical factors were performed using the core function in the ‘corrplot’ package. Mantel test was used for correlational analysis between environmental factors and microbial community composition according to the r and significance level P-value of the two matrices. All statistical analyses were performed using R software (version 3.6.3). Variation partitioning analysis (VPA) was determined the proportion of specific environmental factors to explain the changes of community structure by using “vegan” libraries in R.
The WT ranged from 21.0 to 36.0 °C with an average of 29.0 °C and showed a weak increasing trend from the north to the south (Fig. 2). Salinity ranged from 0.00 to 36.00‰ with an average of 5.81‰ and was usually high in the north wetlands. Trophic level-related parameters ( including TN, TP, Chl a, POC, and PON) usually showed high values in the middle wetlands, whereas DSi seemed decreasing from the north to the south. Other physicochemical parameters had disorderly spatial variation. All these parameters did not show regular patterns based on the wetland types (Fig. 2).
Planktonic bacteria and archaea 16S rRNA gene abundance, α-diversity and community composition
Bacteria 16S rRNA gene abundance was from 0.67 to 31.95 × 1012 copies L−1 with an average of 5.4 × 1012 copies L−1, and that of archaea was from 0.71 to 12.92 × 1012 copies L−1 with an average of 1.8 × 1012 copies L−1. They were usually high in the middle wetlands (Fig. 2). Shannon_e index of bacteria was higher than that of archaea, and both of them showed high values in the south wetlands.
A total of 18,819,251 and 6,736,335 reads were obtained after quality control and then clustered into 53,503 and 21,350 OTUs in bacteria and archaea, respectively. All bacteria OTUs yielded by high-throughput sequencing were clustered into 63 phyla, whereas archaea OTUs only assembled 9 phyla. For bacteria, Proteobacteria and Bacteroidetes were top 2 phyla, and the next phyla were followed by Actinobacteria and Cyanobacteria with relative abundance more than 10% (Fig. 3a). For archaea, Nanoarchaeaota, Thaumarchaeaota, Crenarchaeota and Euryarchaeota were dominant phyla (Fig. 3b). From the north to the south, archaea community structure showed larger variation than that of bacteria (Fig. 3, Additional file 1: Table S1).
Planktonic bacteria and archaea co-occurrence networks
The co-occurrence networks consisted of 357 nodes and 1687 edges (Fig. 4). Most of them were positive correlation. Planktonic archaea were in a core position and had more contact with planktonic bacteria. Sub-networks of each water sample were generated, and the nodes of sub-networks of each wetland decreased with the increasing latitude (Fig. 2).
Bate mean nearest taxon metric
Only 19.4% of the βNTI scores for 16S rRNA OTUs derived from bacteria community were in the range of − 2 to + 2, indicating that deterministic processes dominated planktonic bacteria community assembly in the wetlands (Fig. 6a). However, planktonic archaea βNTI scores in the range of − 2 to + 2 accounted for 59.3% of the total, which was indicative of stochastic processes to be greater than deterministic processes in their community assembly in the wetlands (Fig. 6b). For stochastic processes, dispersal (homogenizing), drift (acting alone), and dispersal (limitation combined with drift) was divided based on the RCbary, and the dispersal (limitation combined with drift) was prominent in planktonic bacteria (19.0%) and archaea (53.4%) community assembly.
The βNTI were grouped according to the spatial distance, salinity, and nodes, respectively (Fig. 5). The planktonic bacteria βNTI showed a decrease with an increase in those parameters, indicating an enhancement of deterministic processes on their community assembly. However, this phenomenon was not found for planktonic archaea in the wetlands. For the niche breadth, planktonic archaea showed greater value than planktonic bacteria (Fig. 6a).
Relationship between community composition and environmental factors
Variation partitioning analysis indicated the unexplained of bacterial community composition and archaeal community composition was accounted for 67% and 56%, respectively (Fig. 7b, c). Salinity and node were found to be key influencing factors based on mantel test (Fig. 7a) and showed significant negative correlation with shannon_e index of bacteria and archaea, respectively (Fig. 8).
Mechanisms underlying geographical patterns of planktonic bacteria and archaea in the coastal wetlands
There are different contributions of stochastic and deterministic processes to planktonic bacteria and archaea community assembly in the coastal wetlands along geographical distance. The deterministic selection shapes planktonic bacteria community, while stochastic processes regulate planktonic archaea community assembly (Fig. 5a, b). For the latter, the dispersal combined with drift was suggested to be a dominance in the stochastic processes, and the relative importance of dispersal and ecological drift is stronger than deterministic selection (Fig. 6b). Few studies are about planktonic archaea community assembly across complex coastal wetlands, and a stochastic-dominated archaea assembly in coastal waters of Northern Zhejiang has been reported recently . The contrasting mechanisms for planktonic bacteria and archaea assembly may be attributed to their different niche breadth. The wider niche breadth is less affected by environmental selection , and planktonic archaea are found to be wider niche breadth than planktonic bacteria (Fig. 6a) and thus tend to be stochastic regulation. With an increase in geographical distance, deterministic-dominated processes underlying the planktonic bacteria and archaea assembly tend to strengthen (Fig. 5a, b), probably because an excessive increase in geographical distance would weaken the ability of microbial dispersal and drift , and local deterministic selection gradually dominates community assembly.
In addition, there were two clusters with significant habitat differences: the inlet and outlet waters and the wetland waters (Additional file 1: Fig. S1), suggesting that local ecological niche plays a crucial role in shaping planktonic microbial community composition over a short distance. The inlet and outlet waters are lotic, while the wetland waters have a higher hydraulic retention time. This can result in their different microbial communities because hydrologic conditions have been reported to control nutrient biogeochemical cycling  and shape microbial composition . In addition, wetland vegetation maintain ecosystem structure and function by recycling nutrients, attenuating flow velocities, releasing oxygen and organic carbon, and stabilizing the sediment [46, 47]. This helps to shape the unique environmental conditions of wetlands, where bacteria community structure in the surface waters is similar to that in the bottom waters due to the shallow vertical depth (< 2 m). Different from planktonic bacteria, the similarity of planktonic archaea community structure is generally low, but cluster analysis still show dynamic difference (Additional file 1: Fig. S1). All these differences are attributed to their different environmental factors (Additional file 1: Fig. S2).
Deterministic environmental factors for planktonic bacteria and archaea community in the coastal wetlands
In China’s coastal wetlands, salinity is a key factor in shaping planktonic bacteria community composition; whereas a key determinant for planktonic archaea community assembly is not found (Fig. 7a), which is consistent with its stochastic-dominanted regulation. Salinity has been considered as a stress factor of bacteria community and can affect bacteria respiratory activity, membrance polarization and integrity, and DNA and RNA contents [48, 49]; therefore, salinity has been reported to be a major determinant affecting the biogeographic pattern of bacterial communities across lake and marine ecosystems . For planktonic archaea, although Crenarchaeota, Euryarchaeota and Thaumarchaeota are well known as the dominant archaea in brackish and pond sediment [51, 52], Nanoarchaeaeota are the top one in this study (Fig. 3), indicating their pervasiveness in the coastal wetland waters. In general, deterministic environmental factors only explain less than 44% of planktonic bacteria and archaea community composition in the coastal wetlands (Fig. 7b, c), and other factors such as biological interaction might play an important role in their community assembly.
The abundance of planktonic bacteria and archaea is mainly determined by TP, Chl a, POC and PON in the coastal wetlands (Fig. 7a). It is reasonable because they are nutrient-related environmental factors, and planktonic bacteria and archaea need them to support their growths. These factors are usually positively correlated from each other (Fig. 7a), indicating the importance of phytoplankton as a primary producer to some degree.
Effect of planktonic bacteria-archaea co-occurrence on their biogeographic patterns in the coastal wetlands
To the best of our knowledge, this is the first study that used network analysis to explore planktonic bacteria-archaea co-occurrence in the coastal wetlands. The large number of correlations imply that bacteria and archaea coexist in the same habitat and are functional interdependencies (Fig. 4). Nanoarchaeaeota are obligate symbiont with chemolithotrophic Crenarchaeon Ignicoccus hospitalis that is a sulfur reducer . As such, Nanoarchaeaeota should be involved in sulfur cycling along with Proteobacteria. In addition, Thaumarchaeota and Crenarchaeota are ammonia-oxidizing archaea containing amoA enzyme and can oxidize ammonia to nitrite during nitrification [54, 55]; as such, they will play an important role in nitrogen cycling together with Proteobacteria and Bacteroidetes. These interactions among the specific species finnally result in the tight coorelations in the network analysis although variations in co-occurrence patterns that are derived from a topology-based system approach have been reported not to reflect true inter taxa correlations and only to reveal partial complex interaction within a microbial community [56,57,58]. In addition, the nodes of sub-networks rather than salinity showed greater significant relationship with the shannon_e index of planktonic bacteria and archaea (Fig. 8), and their community composition also exhibited significant correlation with their biological parameters (i.e., shannon_e index, nodes of sub-networks, and 16S rRNA gene abundance) (Fig. 7a). All these demonstrate that planktonic bacteria-archaea co-occurrence play an important role in shaping their biogeographic patterns in the coastal wetlands.
Coastal wetlands are special and complex ecosystems. Planktonic bacteria and archaea in China’s coastal wetlands have different biogeographic patterns and control mechanisms. Deterministic selection shapes the planktonic bacteria community structure, and salinity is a main controlling factor for their community assembly. Stochastic processes regulate the planktonic archaea community structure, being consistent with the fact that planktonic archaea have a larger niche breadth than planktonic bacteria. Planktonic bacteria and archaea co-occur, and their co-occurrence play an important role in their community assembly. This study helps to improve an understanding about biogeographic patterns of planktonic microbes, and thus provides new insight into studying underlying mechanisms of microbial biogeography in coastal wetlands.
Availability of data and materials
The datasets analysed during the current study are available in NCBI SRA database with the accession numbers of PRJNA681135 and PRJNA674461.
Nogales B, Aguilό-Ferretjans MM, Martín-Cardona C, Lalucat J, Bosch R. Bacterial diversity, composition and dynamics in and around recreational coastal areas. Environ Microbiol. 2007;9:19–29.
Fortunato CS, Herfort L, Zuber P, Baptista AM, Crump BC. Spatial variability overwhelms seasonal patterns in bacterioplankton communities across a river to ocean gradient. ISME J. 2012;6:554–63.
Zhang L, Shen T, Cheng Y, Zhao T, Li L, Qi P. Temporal and spatial variations in the bacterial community composition in Lake Bosten, a large, brackish lake in China. Sci Rep-UK. 2020;10:304.
Hanson CA, Fuhrman JA, Horner-Devine MC, Martiny JBH. Beyond biogeographic patterns: processes shaping the microbial landscape. Nat Rev Microbiol. 2012;10:497–506.
Ramette A, Tiedje JM. Biogeography: an emerging cornerstone for understanding prokaryotic diversity, ecology, and evolution. Microb Ecol. 2007;53:197–207.
Martiny JBH, Bohannan BJM, Brown JH, Colwell RK, Fuhrman JA, Green JL, et al. Microbial biogeography: putting microorganisms on the map. Nat Rev Microbiol. 2006;4:102–12.
Du Y, Yang W, Ding X, Zhang J, Zheng Z, Zhu J. High heterogeneity of bacterioplankton community shaped by spatially structured environmental factors in west lake, a typical urban lake in eastern China. Environ Sci Pollut R. 2020;27:42283–93.
Guo L, Wang G, Sheng Y, Sun X, Shi Z, Xu Q, et al. Temperature governs the distribution of hot spring microbial community in three hydrothermal fields, eastern tibetan plateau geothermal belt, western China. Sci Total Environ. 2020;720:137574.
Nio-García JP, Ruiz-González C, Giorgio P. Interactions between hydrology and water chemistry shape bacterioplankton biogeography across boreal freshwater networks. ISME J. 2016;10:1755–66.
Zhou J, Ning D. Stochastic community assembly: does it matter in microbial ecology? Microbiol Mol Biol R. 2017;81:e00002-17.
Ofiţeru ID, Lunn M, Curtis TP, Wells GF, Criddle CS, Francis CA, et al. Combined niche and neutral effects in a microbial wastewater treatment community. P Natl Acad Sci USA. 2010;107:15345–50.
Stegen JC, Lin X, Konopka AE, Fredrickson JK. Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 2012;6:1653–64.
Chave J. Neutral theory and community ecology. Ecol Lett. 2004;7:241–53.
Stegen JC, Lin X, Fredrickson JK, Chen X, Kennedy DW, Murray CJ, et al. Quantifying community assembly processes and identifying features that impose them. ISME J. 2013;7:2069–79.
Vellend M. Conceptual synthesis in community ecology. Q Rev Biol. 2010;85:183–206.
Hubbell SP, Borda-De-Agua L. The unified neutral theory of biodiversity and biogeography: reply. Ecol. 2004;85:3175–8.
Chase JM, Myers JA. Disentangling the importance of ecological niches from stochastic processes across scales. Phil Trans R Soc B. 2011;366:2351–63.
Liu L, Wang H, Yue Q. China’s coastal wetlands: ecological challenges, restoration, and management suggestions. Reg Stud Mar Sci. 2020;37:101337.
Xu W, Fan X, Ma J, Pimm SL, Kong L, Zeng Y, et al. Hidden loss of wetlands in China. Cur Biol. 2019;29:1–7.
Wang C, Tang S, He X, Ji G. The abundance and community structure of active ammonia-oxidizing archaea and ammonia-oxidizing bacteria shape their activities and contributions in coastal wetlands. Water Res. 2020;171:115464.
Bunse C, Pinhassi J. Marine bacterioplankton seasonal succession dynamics. Trends Microbiol. 2017;25:494–505.
Hutchins DA, Fu F. Microorganisms and ocean global change. Nat Microbiol. 2017;2:17058.
Campbell BJ, Kirchman DL. Bacterial diversity, community structure and potential growth rates along an estuarine salinity gradient. ISME J. 2013;7:210–20.
Barreto CR, Morrissey EM, Wykoff DD, Chapman SK. Co-occurring mangroves and salt marshes differ in microbial community composition. Wetlands. 2018;38:497–508.
Karamipour N, Fathipour Y, Mehrabadi M. Gammaproteobacteria as essential primary symbionts in the striped shield bug, Graphosoma Lineatum (Hemiptera: Pentatomidae). Sci Rep. 2016;6:33168.
Zu Q, Zhong L, Deng Y, Shi Y, Wang B, Jia Z, et al. Geographical distribution of methanogenic archaea in nine representative paddy soils in China. Front Microbiol. 2016;7:1447.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.
Lan Y, Wang Q, Cole JR, Rosen GL. Using the RDP classifier to predict taxonomic novelty and reduce the search space for finding novel organisms. PLoS ONE. 2012;7:e32491.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microb. 2009;75:7537–41.
Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14:927–30.
Jiao S, Yang Y, Xu Y, Zhang J, Lu Y. Balance between community assembly processes mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 2020;14:202–16.
Zhang W, Lei M, Li Y, Wang P, Wang C, Gao Y, et al. Determination of vertical and horizontal assemblage drivers of bacterial community in a heavily polluted urban river. Water Res. 2019;161:98–107.
Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, et al. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4.
Xun W, Li W, Xiong W, Ren Y, Liu Y, Miao Y, et al. Diversity-triggered deterministic bacterial assembly constrains community functions. Nat Commun. 2019;10:3833.
Dini-Andreote F, Stegen JC, van Elsas JD, Salles JF. Disentangling mechanisms that mediate the balance between stochastic and deterministic processes in microbial succession. P Natl Acad Sci USA. 2015;112:e1326-1332.
Wang K, Razzano M, Mou X. Cyanobacterial blooms alter the relative importance of neutral and selective processes in assembling freshwater bacterioplankton community. Sci Total Environ. 2020;706:135724.
Zhou L, Zhou Y, Hu Y, Cai J, Liu X, Bai C, et al. Microbial production and consumption of dissolved organic matter in glacial ecosystems on the Tibetan Plateau. Water Res. 2019;160:18–28.
Liu J, Zhu S, Liu X, Yao P, Ge T, Zhang X-H. Spatiotemporal dynamics of the archaeal community in coastal sediments: assembly process and co-occurrence relationship. ISME J. 2020;14:1463–78.
Wang K, Hu H, Yan H, Hou D, Wang Y, Dong P, et al. Archaeal biogeography and interactions with microbial community across complex subtropical coastal waters. Mol Ecol. 2019;28:3101–18.
Wu W, Lu H-P, Sastri A, Yeh Y-C, Gong G-C, Chou W-C, et al. Contrasting the relative importance of species sorting and dispersal limitation in shaping marine bacterial versus protist communities. ISME J. 2018;12:485–94.
Heino J, Melo AS, Siqueira T, Soininen J, Valanko S, Bini LM. Metacommunity organisation, spatial extent and dispersal in aquatic systems: patterns, processes and prospects. Freshwater Biol. 2015;60:845–69.
Liang X, Xing T, Li J, Wang B, Wang F, He C, et al. Control of the hydraulic load on nitrous oxide emissions from cascade reservoirs. Environ Sci Technol. 2019;53:11745–54.
Yang M, Shi J, Wang B, Xiao J, Li W, Liu C-Q. Control of hydraulic load on bacterioplankton diversity in cascade hydropower reservoirs, southwest China. Microb Ecol. 2020;80:537–45.
Moustafa M, Wang N. Assessment of wind and vegetation interactions in constructed wetlands. Water. 2020;12:1937.
Rejmánková E. The role of macrophytes in wetland ecosystems. J Ecol Biol. 2011;34:333–45.
del Giorgio PA, Bouvier TC. Linking the physiologic and phylogenetic successions in free-living bacterial communities along an estuarine salinity gradient. Limnol Oceanogr. 2002;47:471–86.
Eswaran R, Khandeparker L. Influence of salinity stress on bacterial community composition and beta-glucosidase activity in a tropical estuary: elucidation through microcosm experiments. Mar Environ Res. 2020;159:104997.
Logares R, Lindström ES, Langenheder S, Logue JB, Paterson H, Laybourn-Parry J, et al. Biogeography of bacterial communities exposed to progressive long-term environmental change. ISME J. 2013;7:937–48.
Behera P, Mohapatra M, Kimb JY, Rastogi G. Benthic archaeal community structure and carbon metabolic profiling of heterotrophic microbial communities in brackish sediments. Sci Total Environ. 2020;706:135709.
Silveira R, Silveira Sartori Silva MR, Bandeira de Mello TR, Cunha Carvalho Alvim EA, Santos-Marques NC, Kruger RH, et al. Bacteria and archaea communities in cerrado natural pond sediments. Microb Ecol. 2020;81:563–78.
Huber H, Hohn MJ, Rachel R, Fuchs T, Wimmer VC, Stetter KO. A new phylum of archaea represented by a nanosized hyperthermophilic symbiont. Nature. 2002;417:63–7.
Auguet J-C, Triadό-Margarit X, Nomokonova N, Camarero L, Casamayor EO. Vertical segregation and phylogenetic characterization of ammonia-oxidizing Archaea in a deep oligotrophic lake. ISME J. 2012;6:1786–97.
Yakimov MM, Cono VL, Smedile F, DeLuca TH, Juárez S, Ciordia S, et al. Contribution of crenarchaeal autotrophic ammonia oxidizers to the dark primary production in Tyrrhenian deep waters (Central Mediterranean Sea). ISME J. 2011;5:945–61.
Ma B, Wang H, Dsouza M, Lou J, He Y, Dai Z, et al. Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. ISME J. 2016;10:1891–901.
Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol. 2012;10:538–50.
Blanchet FG, Cazelles K, Gravel D. Co-occurrence is not evidence of ecological interactions. Ecol Lett. 2020;23:1050–63.
This work was financially supported by the National Natural Science Foundation of China (U1612441) and the National Key R & D Program of China (2016YFA0601001).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Fig. S1. Cluster analysis of planktonic bacteria and archaea based on community composition similarity. Fig. S2. Correlations between physicochemical factors and community characteristic parameter of planktonic bacteria and archaea. Table S1. Sampling information of twenty-one China's coastal wetlands.
About this article
Cite this article
Wang, B., Liu, N., Yang, M. et al. Co-occurrence of planktonic bacteria and archaea affects their biogeographic patterns in China’s coastal wetlands. Environmental Microbiome 16, 19 (2021). https://doi.org/10.1186/s40793-021-00388-9
- Planktonic bacteria and archaea
- Coastal wetlands