Macro and Microelements Drive Diversity and Composition of Prokaryotic and Fungal Communities in Hypersaline Sediments and Saline-Alkaline Soils.

Understanding the effects of environmental factors on microbial communities is critical for microbial ecology, but it remains challenging. In this study, we examined the diversity (alpha diversity) and community compositions (beta diversity) of prokaryotes and fungi in hypersaline sediments and salinized soils from northern China. Environmental variables were highly correlated, but they differed significantly between the sediments and saline soils. The compositions of prokaryotic and fungal communities in the hypersaline sediments were different from those in adjacent saline-alkaline soils, indicating a habitat-specific microbial distribution pattern. The macroelements (S, P, K, Mg, and Fe) and Ca were, respectively, correlated closely with the alpha diversity of prokaryotes and fungi, while the macronutrients (e.g., Na, S, P, and Ca) were correlated with the prokaryotic and fungal beta-diversity (P ≤ 0.05). And, the nine microelements (e.g., Al, Ba, Co, Hg, and Mn) and micronutrients (Ba, Cd, and Sr) individually shaped the alpha diversity of prokaryotes and fungi, while the six microelements (e.g., As, Ba, Cr, and Ge) and only the trace elements (Cr and Cu), respectively, influenced the beta diversity of prokaryotes and fungi (P < 0.05). Variation-partitioning analysis (VPA) showed that environmental variables jointly explained 55.49% and 32.27% of the total variation for the prokaryotic and fungal communities, respectively. Together, our findings demonstrate that the diversity and community composition of the prokaryotes and fungi were driven by different macro and microelements in saline habitats, and that geochemical elements could more widely regulate the diversity and community composition of prokaryotes than these of fungi.

Introduction

Microbial communities are ubiquitous and can even thrive in extreme environments such as hypersaline lakes, marine habitats and saline–alkaline soils (Jiang et al., 2007; Andrei et al., 2015; Yakimov et al., 2015; Xie et al., 2017). They are key components of these extreme ecosystems, and thus play central roles in their geochemical cycles and ecological stability (Nunoura et al., 2013; Liang et al., 2014; Sorokin et al., 2015). Microbial communities in extreme environments have formed unique community structures (Oren, 2013). Moreover, the species diversity and distribution patterns of microbial assemblages vary, together with dynamic ecological variables such as pH, geochemical elements, moisture and terrestrial locations (Xiong et al., 2012; Liu K. et al., 2014; Bryanskaya et al., 2016; Fernandez et al., 2016; Zhong et al., 2016). Revealing the linkage between the diversity and structures of microbial communities and environmental driving forces will provide significant insight into the mechanisms that microbial communities use to adapt to extreme habitats, community succession and ecological functioning.

Accumulating evidence suggests that geochemical elements could be major factors influencing microbial assemblages (Oren, 2013). For example, macroelements such as K+, Ca2+, and Mg2+significantly contributed to the species richness and compositions of prokaryotic communities in soils and hypersaline sediments (Podell et al., 2014; Bryanskaya et al., 2016; Xia et al., 2016; Zhong et al., 2016) and Ca and P significantly affect fungal community compositions in forest soils (Sun et al., 2016). Although some microelements are regarded as toxic to most forms of life, they regulated microbial communities in extreme conditions. For instance, they are significantly correlated with bacterial community composition in hot springs (Jiang et al., 2016) and certain elements (Co, Ni, and Mn) explained variations of prokaryotic and fungal assemblages in polymetallic mining areas (Reith et al., 2015). These previous studies have enriched our knowledge of microbial biogeography, but the respective reports are mainly focused on various macroelements. In fact, microbial communities are jointly driven by various factors. What’s more, these two types of elements drive different physiological processes of microbial life (Stolz et al., 2006; Wintsche et al., 2016). Therefore, it is important to investigate the effects of these multiple variables on the diversity and composition of microbial communities. However, the relevant knowledge is still unavailable especially as regards saline habitats.

Salt lakes and saline–alkaline lands are widely distributed across the earth, where they are characterized by high osmotic pressure and high ionic concentrations (Oren, 2002; Hallsworth et al., 2007). Microbial communities in saline habitats are adapted to high ionic conditions, and salt ions play an important role in driving the diversity and composition of microbial communities. In this study, we collected 36 samples (including hypersaline sediments and saline soils) from salt lakes and saline–alkaline regions of northwest China. The diversity and composition of prokaryotic and fungal communities in these saline samples were determined based on the sequencing data of 16S rRNA genes and fungal ITS regions on an Illumina HiSeq. Specifically, this study has two major objectives: (1) evaluate the influence of macro- and microelements on the diversity of prokaryotic and fungal communities; (2) assess the relative importance of macro- and microelements in driving the distribution of prokaryotic and fungal communities.

Materials and Methods

Site Description and Sample Collection

A total of 27 hypersaline sediment samples were collected from three salt lakes (Gouchi, Yuncheng, and Luyang lakes), and nine salinized soil samples were collected from areas adjacent to these lakes (Figure ). Gouchi and Luyang lakes represent typical chloride-type lakes and are located in Dingbian County and Weinan City of Shaanxi Province in Northern China. Lake Yuncheng, known as one of the world’s three largest sodium-sulfate-type inland saline lakes, is located in Yuncheng City, in Shaanxi Province. Sediment samples from Gouchi, Yuncheng and Luyang lakes were numbered with the prefixes G, Y and W, respectively, and saline soil samples surrounding these three salt lakes were assigned the prefixes GT, YT and WT, respectively. Each sample was obtained by mixing five 10-g subsamples from 1 m × 1 m quadrants, which were taken from a depth of approximately 10–15 cm from a sediment or soil surface. Each set of five subsamples was combined as one sample, giving three biological replicates per plot. All samples were collected in sterile 50-ml plastic screw-top tubes, transported to the laboratory on ice, and stored at -80°C until the extraction of DNA.

Geographic location of the Gouchi salt lake, Yuncheng salt lake, and Luyang salt lake. Location map of the sampling sites was created using ArcGIS software by Esri.

Measurements of Environmental Factors

Macroelements and trace elements in the soil and sediment samples were measured using an inductively coupled plasma mass spectrometer (ICP-MS; Perkin-Elmer Sciex Elan+ 5000). Briefly, sediment and soil samples were dried in a hot-air oven at 40°C until a stable weight was achieved, and sieved through a 0.25-μm filter. Each test sample (200 mg) was fully digested in a microwave digestion system with the addition of a mixture of HNO3 and HF (Arslan and Tyson, 2008), and the resulting solutions were transferred to 25-ml volumetric flasks and diluted to the fixed volume with 3% HNO3. The digestion procedure was performed in triplicate for each test sample. The calibration curve was prepared using a working standard solution with concentrations ranging from 0.01 to 10 μg ml-1 for all macroelements (Na, Mg, K, Ca, S, P, and Fe) and microelements (Al, As, B, Ba, Cd, Co, Cr, Cu, Ga, Ge, Hg, Li, Mn, Ni, Pt, Sr, Ti, Zn, and Rb). The concentrations of macroelements and trace elements in the soil and saline sediments were determined by ICP-MS (Sandroni and Smith, 2002). In addition, moisture was reported based on the water content of the samples. Geographic locations of each site were recorded in situ, and pH values of the samples were measured using pH strips with 1:5 (wt/vol) sediment to water.

DNA Extraction, PCR Amplification, and Illumina-Based Sequencing

DNA was extracted from 0.5-g samples of soil or saline sediments using a FastDNATM SPIN Kit for soil (MP Biomedicals, LLC) according to the manufacturer’s instructions. Triplicate DNA extracts from the same sample were combined and analyzed using a NanoDrop® TM ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE, United States). The 16S rRNA genes and the ITS regions of the extracted DNA were amplified using general bacterial primers 515F/806R (Caporaso et al., 2011) and fungal-specific primers ITS1F (Gardes and Bruns, 1993)/ITS4 (White et al., 1990). PCR was performed in triplicate for each sample with the following thermal cycling: an initial denaturation at 94°C for 5 min, followed by 25 cycles of 30 s at 94°C, 30 s at 50°C (ITS) or 55°C (16S rRNA) and 30 s at 72°C. Replicates were further pooled, and the resulting amplicons from each sample were sequenced on the HiSeq2500 system (Illumina Inc., San Diego, CA, United States).

Sequence Processing

The quality filtering, chimera removal and operational taxonomic unit (OTU)-based clustering of the sequences were performed using the Mothur software V 1.35.1 (Quast et al., 2013). Briefly, sequences with a barcode ambiguity at one end and more than one ambiguous nucleotide (N), and those shorter than 200 bp in length after removal of the barcode and primer were discarded. Chimeric sequences were removed using the UCHIME modules in Mothur. The quality-filtered sequences were truncated to a constant length and were then clustered into OTUs with a 3% dissimilarity cutoff using the UPARSE pipeline (Edgar, 2013). The OTUs that contained fewer than two reads were excluded from further analysis. A representative sequence from each OTU was assigned to bacteria, archaea, and fungi using the Ribosomal Database Project (RDP) classifier (Wang et al., 2007). At the phylum-taxonomy level, OTUs in each of the three samples from the same saline pools or soil plots were combined and assigned, for example, (G1, G2, and G3) to G_1, (G4, G5, and G6) to G_2, (G7, G8, and G9) to G_3, (Y11, Y12, and Y13) to Y_1, (Y21, Y22, and Y23) to Y_2, (Y31, Y32, and Y33) to Y_3, (GT1, GT2, and GT3) to GT, (YT1, YT2, and YT3) to YT and (WT1, WT2, and WT3) to WT.

Multivariate Statistical Analysis

All statistical analyses in this study were performed using the vegan package in R software, except where otherwise noted. Prior to data analysis, the number of sequences was normalized by randomly sub-sampling reads of prokaryotes (30,126 sequences) and fungi (30,198 sequences) for each sample to avoid a potential bias caused by different sequencing depths. The alpha and beta diversities were calculated using species-level OTUs in QIIME (Caporaso et al., 2010). Relationships of environmental variables and alpha-diversity indices (species richness and Shannon) with factors were tested based on Pearson’s correlation. Principal component analyses (PCoA), based on unweighted UniFrac metrics, were performed to reveal the spatial patterns of the microbial communities within the samples. Canonical correspondence analysis (CCA) was used to identify the significant factors (P ≤ 0.05) to the community composition of prokaryotes and fungi. The permutational multivariate analysis of variance (PerMANOVA) was performed using the distance matrices (R: adonis) in vegan with 999 permutations, and variation partitioning analyses (VPA) were performed by the vegan package in R 2.14.0 to evaluate the relative contribution of significant factors (P ≤ 0.05, variance inflation coefficient < 10) to the prokaryotic and fungal community structures. Due to the microbial community shaped by multiple factors, we also analyzed the joint effects of macroelements, microelements and other factors on the community composition.

Accession Numbers

The raw sequencing data generated in this study were deposited into the National Center for Biotechnology Information (NCBI) database under accession numbers SAMN07764381-SAMN07764416 (16S rRNA) and SAMN07764417-SAMN07764452 (ITS).

Results

Characteristics of Environmental Variables

The macroelemental contents differed significantly among the 36 samples of salt-lake sediments and saline–alkaline soils (Supplementary Table S1). The Na content in all samples was variable. The concentrations of P, Ca, K, and Fe were low in samples G1, G2, and G3 compared to the other 33 samples, as expected. Moreover, the average P concentration in the saline sediments from Gouchi Lake (G4 to G9) was higher than that in the adjacent soils, whereas the average P concentrations in the hypersaline sediments from Yuncheng Lake and Luyang Lake were lower than those in the adjacent soils. In contrast, the average concentrations of S and Mg in the sediments from the salt lakes were significantly higher than in the neighboring saline soils.

As shown in Supplementary Table S1, the microelemental contents (Al, Co, Cr, Ga, Hg, Mn, Ni, Sr, and Ti) of the 33 samples were higher than those of G1, G2, and G3. The concentrations of Al and Ba fluctuated greatly among all samples. Moreover, the average Al content was found to be greater in the saline soils surrounding Yuncheng and Luyang lakes than in the sediments from those salt lakes. The Sr content varied significantly in saline sediments from each lake but was relatively constant in the soil samples. Mn and Ti concentrations were relatively constant in the all samples, except for samples G1, G2, and G3, in which they were exceptionally low. The contents of Co, Cr, Ga, Hg, and Ni were relatively constant in all samples. In general, moisture content in saline sediments was greater than in their corresponding surrounding soils.

The Pearson’s correlations showed significant relationships between many environmental variables of the investigated sites (Figure and Supplementary Table S2). For example, the presence of S positively correlated with that of Cr (P = 0.005, r = 0.327) and Ga (P = 0.004, r = 0.310). The P content was significantly correlated with Na, K, Ca, Fe, Al, Co, Cr, Ga, Hg, Li, Mn, Ni, Rb, Ti, latitude and longitude, with P = 0.001 and r ranging from 0.213 to 0.740 for all correlations. The concentrations of Ca were correlated with Fe, Al, Co, Cr, Hg, Li, Mn, Rb, Sr, latitude and longitude (P = 0.001), with r values ranging from 0.420 to 0.541. The K content was correlated with Fe, Al, Co, Cr, Ga, Li, Mn, Ni, Rb and Ti, with P = 0.001 and r values of 0.491–0.809. The presence of Ba was positively related to Cd (P = 0.001 and r = 0.782). The Cr content had an even higher correlation to Na, Fe, Al, Co, Ga, Hg, Li, Mn, Ni, Rb, Ti, longitude and latitude (P = 0.001 and r values in the range of 0.258–0.959). The Ga content was correlated with Fe, Al, Co, Hg, Li, Mn, Ni, Rb and Ti at (P = 0.001 and r values of 0.616–0.959). The Ti content correlated with Na, Fe, Al, Co, Hg, Li, Mn, Ni, and Rb, with P = 0.001 and r values of 0.555 to 0.847. Moisture correlated with S (P = 0.015, r = 0.174) and Ba (P = 0.005, r = 0.256).

The Pearson’s correlation of geochemical variables of the study sites at the Gouchi salt lake, Yuncheng salt lake, and Luyang salt lake. (A,B) Represent P-value and correlation coefficient, respectively.

Distribution Patterns of Microbial Communities

Illumina sequencing yielded 1,152,000 prokaryotic reads and 1,156,000 fungal sequences from the 36 samples, after removing low-quality tags and chimeras, resulting in a total of 43,579 OTUs and 3,182 OTUs, respectively (Supplementary Table S3). We found that the Chao1 index and observed OTUs for prokaryotes were significantly higher than these for fungi, and that the diversity of prokaryotes and fungi was higher in saline soils than in hypersaline samples (Supplementary Table S4). We normalized the sequence numbers to 30,124 and 30,198 reads per sample for prokaryotes and fungi, respectively, based on the fewest reads available among the 36 samples. These normalized reads were used in the subsequent analyses.

The phylum-level distribution patterns of prokaryotes differed significantly between the hypersaline sediments and saline–alkaline soils (Figure ). A total of 74 prokaryotic phyla, including 13 abundant groups, were found across all samples. Proteobacteria were quite abundant in hypersaline sediments, accounting for 57–95% of the total reads per sediment sample (other than G_1, which was dominated by Euryarchaeota), while they exhibited lower richness in the soils. By contrast, Actinobacteria were the most abundant in saline soil samples (14–39%) but represented a relatively small proportion of prokaryotic reads in saline samples (<4%). This phylum was not detected in G_1. Bacteroidetes were more abundant in saline sediments than in soil samples, with the exception of sample Y_2, and Firmicutes were only present in some saline samples. However, Chloroflexi and Acidobacteria were richer in the soils than in the several salt-lake sediment samples.

Relative abundance of the dominant phyla in the prokaryotic (A) and fungal (B) communities from the soil-saline sediment regions. Observed OTUs in three samples from the same saline pool (or soil plot) were combined and further assigned to G_1, G_2, G_3, Y_1, Y_2, Y_3, W_1, W_2, W_3, GT, YT and WT. Only the phyla with >1% mean abundance were shown.

Similarly, the community composition of fungi at the phylum level was distinct between the saline sediments and soils (Figure ). The relative abundances of Ascomycota varied greatly in the saline sediment samples (8–87%), while they were more consistent in the different soil samples (66–86%). Basidiomycota was more abundant in hypersaline sediments (9–60%) than in salinized soils (<3%), with the exception of sample W_1, which did not contain this group. Eukaryota inhabited all investigated sites but could not be assigned to any known group at the phylum level.

The PCoA analyses showed that prokaryotic (Figure ) and fungal (Figure ) communities in all samples were sorted by habitat, which was demonstrated by the community distributions at the phylum-level (Figures ). For example, the prokaryotic communities in hypersaline sediments were better separated along the PC1 axis, while these groups in saline soils were partly arrayed along the PC2 axis. The fungal communities in the sediments and saline soils exhibited the similar separation trend. Moreover, prokaryotic and fungal communities from the same saline pools or plots tended to be assembled together.

Principal coordinates analysis (PCoA) of the prokaryotic (A) and fungal (B) communities in all the samples along a salinity gradient. Microbial communities in three samples from the same saline pool (or soil plot) were represented by G_1, G_2, G_3, Y_1, Y_2, Y_3, W_1, W_2, W_3, GT, YT, and WT, respectively.

Effects of Environmental Factors on Microbial Alpha Diversity

The rarefaction curves for all samples approached the saturation plateau (Supplementary Figure S1), indicating that the sequencing depths for each sample were sufficient to cover the microbial diversity. The species diversity of the prokaryotic community was higher in the soil samples than in the corresponding hypersaline sediments; the fungal community showed similar trends, with the exception of samples G3, Y13, and Y31 (Supplementary Table S4). In addition, species diversity of prokaryotic and fungal microbes was highly variable in the hypersaline sediment samples, with OTUs ranging from 74 to 2,394 and 12 to 366, respectively.

The Pearson’s analyses showed that the alpha-diversity of prokaryotic and fungal communities was influenced differently by macroelements (Table ). The macroelements S (P = 0.001, r = 0.562 and 0.610, respectively), Mg (P = 0.002, r = 0.156 and 0.350, respectively) and Fe (P ≤ 0.047, r = 0.238 and 0.136, respectively) were significantly correlated with the Chao1 (species richness) and Shannon (diversity and evenness) indices of the prokaryotic community. Meanwhile, P (P = 0.021, r = 0.100) and K (P = 0.007, r = 0.118) were positively correlated with the Chao1 index of these groups. In contrast, Ca (P = 0.004, r = 0.228) correlated with the Shannon index of the fungal community.

As shown in Table , the microelements had a greater effect on the species diversity of prokaryotic assemblages than on those of fungal communities. The microelements (Al, Ba, Co, Cr, Ga, Ni, and Ti) were significantly correlated with the Chao1 (P ≤ 0.004, r = 0.231 ± 0.064) and Shannon indices (P ≤ 0.025, r = 0.192 ± 0.033) of the prokaryotic community. Concentrations of Hg and Mn significantly influenced the Chao1 of prokaryotic taxa (P = 0.001, r = 0.167 and 0.218, respectively). However, for the fungal community, Ba and Cd were positively related to Chao1 (P ≤ 0.017, r = 0.467 and 0.362, respectively), and Sr was significantly correlated with Shannon index (P = 0.001, r = 0.433) as well.

Of other environmental variables, moisture was significantly correlated with the Chao1 index of the prokaryotic and fungal communities (P ≤ 0.005, r = 0.158 and 0.275, respectively). And, longitude and latitude were linked to the Chao1 (P ≤ 0.020, r = 0.129 and 0.160, respectively) and Shannon indices (P ≤ 0.047, r = 0.162 and 0.107, respectively) of the prokaryotic assemblages.

Effects of Environmental Factors on Microbial Beta Diversity

The CCA results showed significant correlations (P ≤ 0.05) between macro and microelements (e.g., Na, S, P, Ca, Ba, Cr, and Ge) and the prokaryotic microbial assemblages (Figure ), indicating that these geochemical parameters have strong impacts on the composition of the prokaryotic communities. In contrast, the CCA analysis also indicated that the trends of the factors, e.g., S, P, Ca, Cr, Cu, and moisture, correlated with the two axes and that these nutrient components significantly influenced the fungal community structures (P ≤ 0.05) (Figure ).

Canonical correspondence analysis (CCA) of the prokaryotic (A) and fungal (B) communities with the significant variables (P ≤ 0.05). G, Y, W, GT, YT, and WT represent microbial communities from Gouchi, Yuncheng, Luyang lakes, and the surrounding soils.

The PerMANOVA results indicated that the compositions of the prokaryotic communities in all samples were strongly driven by macroelements, including Na, S, P and Ca, with P-values ranging from 0.001 to 0.023 and r2-values ranging from 0.036 to 0.192. They were influenced by microelements As, Ba, Cr, Cu, Ge and Zn, with P-values varying from 0.001 to 0.048 and r2-values varying from 0.030 to 0.075. In addition, pH (P = 0.049, r2 = 0.031) correlated significantly with the compositions of the prokaryotic communities (Table ). However, the fungal community structures were significantly influenced by the macroelements, including Na (P = 0.007, r2 = 0.041), S (P = 0.001, r2 = 0.074), P (P = 0.004, r2 = 0.043), Ca (P = 0.016, r2 = 0.036), and K (P = 0.003, r2 = 0.043), and by the microelements Cr (P = 0.006, r2 = 0.039) and Cu (P = 0.008, r2 = 0.039). Latitude shaped the compositions of fungal communities (P = 0.038, r2 = 0.033) (Table ).

Determinants of the Microbial Community Composition

Variance partitioning analyses were used to further quantify the relative effects of environmental factors on prokaryotic and fungal community structures. Macroelements (Na, S, P, and Ca), microelements (As, Ba, Cr, Cu, Ge and Zn) and pH contributed 19.34, 21.28, and 3.23% to the prokaryotic community variation, respectively. The combination of the significant geochemical variables (Na, S, P, Ca, As, Ba, Cr, Cu, Ge, Zn, and pH) (P ≤ 0.05) explained 55.49% of the prokaryotic community variation, leaving 44.51% of the variation unexplained (Figure ). In contrast, macronutrients (Na, S, P, Ca, and K), trace elements (Cr and Cu), other variables (latitude and moisture) accounted for 17.13%, 6.93% and 6.79% of the prokaryotic community compositions, respectively. The combination of the significant environmental factors (Na, S, P, Ca, K, C, Cu, moisture and latitude) (P ≤ 0.05) accounted for 32.27% of the observed variation in the fungal community, leaving 67.73% of the variation unexplained (Figure ).

Variance partitioning analysis (VPA) demonstrating the effects of environmental factors on the variation of the prokaryotic (A) and fungal (B) community compositions.

Discussion

In this study, we found that the average diversity index (Shannon index) of prokaryotes and fungi in the hypersaline sediments of salt lakes was consistently lower than in the adjacent saline–alkaline soils, and examined a noticeable decreasing trend in average species richness (Chao1 index) of the microbial communities in the hypersaline sediments (Supplementary Table S4), suggesting that the extreme saline environments harbor the lowest alpha-diversity of prokaryotes and fungi. It was agreed well with the general principle for ecology that microbial diversity is low in extreme habitats (Benlloch et al., 2002; Pavloudi et al., 2017). Furthermore, we found that the species richness and diversity of prokaryotes and fungi in the sediments drastically changed but remained relatively constant in the soils. This fluctuation of alpha-diversity indices might be caused by geochemical variations in the hypersaline sediments (Supplementary Table S1).

As expected, prokaryotic and fungal communities in hypersaline sediments and saline soils exhibited a habitat-dependent distribution pattern. Proteobacteria dominated hypersaline sediments (except for G_1), whereas many classes, such as Proteobacteria, Actinobacteria, Chloroflexi and Acidobacteria, were relatively evenly distributed in the saline soils (Figure ). Some groups of Proteobacteria widely occupied hypersaline habitats (Hollister et al., 2010; Quaiser et al., 2011; Kambura et al., 2016). Similarly, Ascomycota dominated the fungal communities in the soil samples, while its abundance varied obviously in the sediments. Basidiomycota were more abundant in the various hypersaline sediments and less in the soils (Figure ). Notably, prokaryotic and fungal communities were similar to each other when their habitats (hypersaline sediments versus saline soils) were similar, even if geographic distance between sample sites was greater, suggesting that habitat types play an important role in driving the microbial beta diversity.

The findings in this study demonstrate that the macroelements S, P, K, Mg, and Fe were drivers of the alpha-diversity of the prokaryotic communities in the salt-lake sediments and saline soils (P ≤ 0.05) (Table ). Previous reports have mainly addressed that elements K+, S2-, and Mg2+ determine the species richness of archaeal communities in a hypersaline lake (Podell et al., 2014), and that K+, Mg2+, and total P affect the OTU richness of soil bacterial groups such as Chloroflexi, Actinobacteria, Proteobacteria, and Bacteroidetes (Kim et al., 2016; Xia et al., 2016). Aside from the macronutrients (S, P, K, and Mg) associated with the diversity, we also explored the importance of Fe driving the prokaryotic alpha diversity. Here, we only analyzed the correlation of Ca (P = 0.004) with the Shannon index of the fungi in saline habitats (Table ). Ca as well was found to regulate global soil fungal diversity (Tedersoo et al., 2014). It might be because Ca2+-mediated signals pathways are indispensable for eukaryotes’ life processes including salt-stress tolerance (Shi et al., 2011; Asano et al., 2012; An et al., 2014). In contrast, K+ is crucial to the survival of halotolerant and halophilic prokaryotes in saline environments because these microbes concentrate K+ inside cells to maintain their abundance (Edbeib et al., 2016). This could partially explain why the alpha-diversity of prokaryotes and fungi in saline habitats are driven by different macroelements.

Studies have showed that the geochemical conditions of habitats govern microbial community compositions (Wang et al., 2013, 2017; Hazard et al., 2014; Oloo et al., 2016). Here, we observed that macroelements (Na, S, P, and Ca) shaped prokaryotic community compositions across the studied sites (P ≤ 0.05) (Table ). These macroelements explained 19.34% of the total variation in prokaryotic community structures (Figure ). The significance of Na, S, and Ca concentrations for bacterial community structures has been suggested by research conducted on the mine tailings and salt lakes of the Tibetan Plateau (Valentín-Vargas et al., 2014; Zhong et al., 2016). We showed that Na, S, P, Ca, and K were significantly correlated with fungal community compositions (Table ) and that these macronutrients explained 17.13% of the total variation (Figure ). S, P, and Ca were found to regulate fungal community compositions in a variety of soil habitats (Jie et al., 2012; Tedersoo et al., 2014), however, K and Na as drivers of fungal beta diversity were rarely reported in previous studies. Macroelements are responsible for multiple microbial processes, such as the synthesis of biological macromolecules, signal transduction and osmotic balance (Dominguez, 2004; Jung and Bahn, 2009; Madigan et al., 2009; Edbeib et al., 2016); thus, they play a crucial role in driving the community composition of prokaryotes and fungi.

We revealed a series of microelements (Al, Ba, Co, Cr, Ga, Hg, Mn, Ni, and Ti) influencing the alpha-diversity of prokaryotes (Table ). As shown in previous studies, metal elements Al, Cr, Hg, Mn, and Ni shape the species richness of bacteria in soils or in marine sediments (Faoro et al., 2010; Liu Y.R. et al., 2014; Pereira et al., 2014; Quero et al., 2015; Zhang et al., 2015), but correlations between microchemicals (Ba, Co, Ga, and Ti) and prokaryotic alpha diversity had not been previously discerned. Although some microelements (e.g., Co, Cr, Hg, and Mn) are toxic to life, microchemicals with a suitable content serve as electron donors or acceptors and enzymatic activators of bacterial cells (Stolz et al., 2006; Huang et al., 2015; Wintsche et al., 2016). What’ more, it is known that diverse prokaryotes tolerate the effects of certain metals, and that microbial assemblages are significantly more resistant to heavy metals than pure cultures (Oregaard and Sørensen, 2007; Mejias Carpio et al., 2018). Therefore, microelements might be not toxic and instead be drivers of alpha-diversity of prokaryotic groups. Moreover, one study has shown that deep-sea hydrothermal sediments, rich in metal elements such as Ba and Sr, harbored the highest level of species diversity and phylogenetic uniqueness for eukaryotes (López-García et al., 2003), suggesting that some microelements could modulate fungal alpha-diversity. Here, our results showed that Ba correlated with the species richness of fungi (P = 0.001), while Sr correlated with the Shannon diversity, highlighting the importance of alkali-metal elements for the fungal alpha-diversity.

Similarly, microelements explained the variations in prokaryotic community composition (Table ). For example, the elements As, Ba, Cr, Cu, Ge, and Zn significantly contributed to the differences in prokaryotic community structures in both hypersaline sediments and saline soils. Our findings are supported by research showing that microelements, including As, Cr, Cu, and Zn, affected bacterial community structures in coastal sediments and metal-rich soils (Quero et al., 2015; Reith et al., 2015). The combined microelements (As, Ba, Cr, Cu, Ge, and Zn) and variables (Cr and Cu) individually explained 21.28% and 6.93% of the total variation in community structures of prokaryotes and fungi (Figure ), indicating the relative importance of micronutrients in modulating prokaryotic composition. Together, the alpha and beta-diversities of prokaryotes were more widely influenced by macro and microelements as compared to fungi, suggesting that geochemical element profiles could better predict the diversity and community compositions of prokaryotes.

Furthermore, we observed that variables including moisture, longitude and latitude had relationships with the species richness and community compositions of prokaryotes and fungi, supported by previously published accounts (Talley et al., 2002; Tedersoo et al., 2014; Ding et al., 2015; Shay et al., 2015). However, it should be noted that given terrestrial coordinates determine the specific geochemical characteristics, which could have direct effects on microbial diversity and community composition. Moreover, latitudes, on a broad scale, are related to temperature, which is a well-known controller of the primary productivities and nutrient uptake of cells (Fuhrman et al., 2008). The diversity and community compositions of prokaryotes and fungi were driven by different variables, suggesting that distinct mechanisms are involved in the maintenance and succession of these two communities.

Conclusion

We found that the prokaryotic and fungal community compositions in hypersaline sediments and saline–alkaline soils exhibited habitat-dependent patterns. We explored the concept that the diversity and composition of prokaryotic and fungal communities were driven by different macro and microelements. We also found that the diversity and composition of prokaryotes were more widely influenced by geochemical elements than those of fungal communities. Geochemistry selects unique microbial communities and, conversely, microbes drive geochemical cycles. Further studies could focus on the functional microbial communities involved in geochemical recycling in hypersaline environments.

Author Contributions

KL and XD collected the samples and designed the research. KL, XD, and XT produced the data. KL, XD, JW, WL, QY, and ZL performed the data analysis. KL and XD wrote the manuscript.

Conflict of Interest Statement

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.

Table 1

The Pearson correlations of alpha-diversity indices of prokaryotic and fungal communities with environmental factors.

Factor	Prokaryotic community (P/R)			Fungal community (P/R)
	Chao1 (Richness)	Shannon		Chao1 (Richness)	Shannon
Na	0.489/-0.006	0.592/-0.033		0.280/0.042	0.067/0.156
S	0.001/0.526	0.001/0.610		0.320/0.012	0.319/0.029
P	0.021/0.100	0.214/0.049		0.067/0.181	0.188/0.080
Ca	0.223/0.026	0.085/0.088		0.690/-0.056	0.004/0.228
K	0.007/0.118	0.150/0.080		0.340/0.053	0.431/0.006
Mg	0.002/0.156	0.002/0.350		0.646/-0.075	0.158/0.103
Fe	0.001/0.238	0.047/0.136		0.494/-0.027	0.404/0.007
Al	0.001/0.225	0.020/0.165		0.347/0.029	0.356/0.023
As	0.878/-0.049	0.627/-0.042		0.143/0.146	0.776/-0.078
B	0.549/-0.012	0.269/0.033		0.299/0.031	0.784/-0.079
Ba	0.001/0.287	0.012/0.196		0.001/0.467	0.261/0.044
Cd	0.461/0.001	0.553/-0.030		0.017/0.362	0.514/-0.019
Co	0.001/0.245	0.012/0.187		0.414/0.013	0.491/-0.011
Cr	0.001/0.295	0.011/0.201		0.421/-0.012	0.455/-0.003
Cu	0.317/0.017	0.627/-0.037		0.591/-0.071	0.971/-0.156
Ga	0.001/0.270	0.009/0.206		0.342/0.008	0.215/0.064
Ge	0.200/0.033	0.072/0.144		0.322/-0.008	0.969/-0.125
Hg	0.001/0.167	0.159/0.072		0.417/-0.009	0.326/0.024
Li	0.259/0.022	0.475/-0.002		0.332/0.035	0.148/0.101
Mn	0.001/0.218	0.062/0.121		0.337/0.026	0.266/0.047
Ni	0.004/0.167	0.025/0.158		0.462/-0.021	0.216/0.072
Pt	0.918/-0.056	0.454/-0.025		0.524/-0.044	0.574/-0.047
Sr	0.417/0.004	0.177/0.064		0.712/-0.064	0.001/0.433
Ti	0.001/0.282	0.009/0.225		0.225/0.079	0.195/0.062
Zn	0.813/-0.041	0.931/-0.092		0.979/-0.117	0.837/-0.110
Rb	0.136/0.051	0.398/0.01		0.094/0.178	0.139/0.097
TCE	0.985/-0.076	0.899/-0.084		0.177/0.099	0.091/0.122
pH	0.779/-0.035	0.849/-0.07		0.396/0.015	0.390/0.016
Moisture	0.004/0.158	0.057/0.101		0.005/0.275	0.281/0.043
Longitude	0.020/0.129	0.006/0.162		0.930/-0.075	0.287/0.021
Latitude	0.007/0.160	0.047/0.107		0.836/-0.065	0.205/0.046

Table 2

The PerMANOVA of the community compositions of prokaryotes and fungi with geochemical factors.

Factor	R2	P
Prokaryotic community
Na	0.092	0.001
S	0.191	0.001
P	0.047	0.006
Ca	0.036	0.023
Mg	0.007	0.846
As	0.040	0.008
B	0.019	0.251
Ba	0.031	0.048
Cr	0.036	0.029
Cu	0.030	0.048
Ge	0.034	0.014
Ni	0.017	0.277
Zn	0.075	0.001
pH	0.031	0.049
Moisture	0.031	0.066
Fungal community
Na	0.041	0.007
S	0.074	0.001
P	0.043	0.004
Ca	0.036	0.016
K	0.043	0.003
Mg	0.027	0.117
Fe	0.029	0.088
Al	0.021	0.533
As	0.025	0.267
B	0.020	0.592
Ba	0.030	0.081
Cr	0.039	0.006
Cu	0.039	0.008
Ga	0.022	0.401
Ge	0.025	0.251
Hg	0.024	0.306
Li	0.021	0.534
Mn	0.027	0.171
Ni	0.030	0.073
Sr	0.026	0.194
Ti	0.026	0.169
Zn	0.021	0.531
Rb	0.026	0.170
TCE	0.020	0.593
pH	0.022	0.430
Moisture	0.032	0.051
Longitude	0.030	0.061
Latitude	0.033	0.038