Next Article in Journal
SOX11, SOX10 and MITF Gene Interaction: A Possible Diagnostic Tool in Malignant Melanoma
Next Article in Special Issue
Stachybotrys musae sp. nov., S. microsporus, and Memnoniella levispora (Stachybotryaceae, Hypocreales) Found on Bananas in China and Thailand
Previous Article in Journal
Slr0320 Is Crucial for Optimal Function of Photosystem II during High Light Acclimation in Synechocystis sp. PCC 6803
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Soil Fungal and Bacterial Communities in Tianchi Volcano Crater, Northeast China

School of Pharmaceutical Science and Technology, Tianjin University, 92 Weijin Road, Nankai District, Tianjin 300072, China
*
Author to whom correspondence should be addressed.
Life 2021, 11(4), 280; https://doi.org/10.3390/life11040280
Submission received: 16 March 2021 / Revised: 23 March 2021 / Accepted: 24 March 2021 / Published: 26 March 2021

Abstract

:
High-altitude volcanoes, typical examples of extreme environments, are considered of particular interest in biology as a possible source of novel and exclusive microorganisms. We analyzed the crater soil microbial diversity of Tianchi Volcano, northeast China, by combining molecular and morphological analyses of culturable microbes, and metabarcoding based on Illumina sequencing, in order to increase our understanding of high-altitude volcanic microbial community structure. One-hundred and seventeen fungal strains belonging to 51 species and 31 genera of Ascomycota, Basidiomycota and Mucoromycota were isolated. Penicillium, Trichoderma, Cladosporium, Didymella, Alternaria and Fusarium dominated the culturable fungal community. A considerable number of isolated microbes, including filamentous fungi, such as Aureobasidium pullulans and Epicoccum nigrum, yeasts (Leucosporidium creatinivorum), and bacteria (Chryseobacterium lactis and Rhodococcus spp.), typical of high-altitude, cold, and geothermal extreme environments, provided new insights in the ecological characterization of the investigated environment, and may represent a precious source for the isolation of new bioactive compounds. A total of 1254 fungal and 2988 bacterial operational taxonomic units were generated from metabarcoding. Data analyses suggested that the fungal community could be more sensitive to environmental and geographical change compared to the bacterial community, whose network was characterized by more complicated and closer associations.

1. Introduction

Microorganisms are considered an essential component of natural environments [1]. They are ubiquitous in nature based on their characteristics, such as small size, flexible capability to exploit nutrients and adaptability to unfavorable and extreme environmental conditions [1]. Microbial communities in volcanic environments are of particular interest for research on the emergence and evolution of life due to the unique extreme conditions that characterize these natural habitats [2]. Extremophilic microorganisms inhabiting volcanic environments can be a source of important biotechnological products like antibiotics, bio-fertilizers or bio-control agents [3,4,5]. New extremophiles have been isolated from different, very peculiar environments, which are often considered “empty” in the biological sense [6]. Novel microbial species have been recently described from barren alpine environments [7,8] and high-altitude volcanic soils [4,5].
Both morphological culture-based and molecular approaches are generally used to characterize environmental microbial communities. Culturable microorganisms provide more detailed information for their identification at species level, and for understanding their role in environmental processes, but only less than 5% of the microorganisms on Earth are culturable in lab condition, due to the selectivity of media and culture conditions [9,10]. To overcome this limitation, diverse methods have been developed to obtain genetic information on the presence of microorganisms living in natural environments. These methods allow the culture-independent analysis of the total microbial genomes called “metagenome” in a particular environment [2,11,12,13]. The development of molecular biological techniques (particularly high-throughput sequencing) represents a powerful tool to access a much larger proportion of microbial communities, compare to traditional culture-based methods, by means of environmental sample total DNA extraction [14].
Changbai Mountain Nature Reserve is located in the northeast of China. It is considered one of the most well-protected and conserved natural ecosystems in China [15]. Many new microbial species have been described in this area, which is particularly predisposed to harbor uncommon microorganisms due to its distinctive landforms and complex climatic conditions [16]. For instance, a novel bacterial strain Bacillus methylotrophicus was isolated from the soil near the roots of Pinus koraiensis plants at an altitude of 2749 m a.s.l. [4], while a novel thermophilic anaerobic bacterium Fervidobacterium changbaicum was collected from the mixture of water and mud from a hot spring, in Changbai Mountain [17]. A number of novel macro-mycete taxa were also described from the highly diverse fungal community of this Nature Reserve [18,19]. Tianchi Volcano is located at the highest point of Changbai Mountain ranges [20]. It is an active volcano with the highest potential eruption risk in China [21,22]. Water accumulates at the top of the volcanic cone and forms Tianchi lake, the highest volcanic lake in China, at an altitude of 2189 m a.s.l. The lake is surrounded by 16 peaks, the highest peak reaching the altitude of 2749.5 m a.s.l. [23,24]. The crater of Tianchi Volcano is the accumulation area of volcanic lavas and debris, which are mainly composed by Ti, Fe, Mn, Si, Al, Ca, Na, K and Mg oxides, trace elements, including high concentration of rare-earth elements, and heavy metals, such as Zn, Pb and Cu [23,25].
To date, a few studies on microbial diversity have been performed in Changbai Mountain Nature Reserve, mainly focusing on the structure and function of the soil microbial communities in vertical vegetation zones, with little attention to microbial communities in the high crater area of Tianchi Volcano. For instance, the fungal and bacterial diversity and community composition along an elevation gradient (2000–2500 m a.s.l.) on the northern alpine tundra belt of Changbai Mountain were analyzed using Illumina sequencing [26,27]. Other studies focused on specific soil microbial diversity associated with vegetation zones [28,29,30,31]. Zhao et al. characterized microbes in rhizosphere soils of Cowskin Azalea (Rhododendron aureum) using culture-dependent methods [32]. The diversity of culturable forest micro-fungi in different vegetational belts from 700 to 2600 m a.s.l. was described by Yang et al. [16]. In this study, we focused on the mountaintop area of Tianchi Volcano, which was expected to harbor a peculiar microbial community adapted to the particularly extreme environmental conditions characterizing the highest elevation montane zones, including colder temperatures, and higher exposure to wind and solar radiation. We performed a comprehensive analysis of the soil fungal and bacterial diversity colonizing the crater margin of Tianchi Volcano, using a combination of molecular and morphological analyses of culturable microbes, and metabarcoding based on Illumina sequencing, in order to shed light on the structure of the microbial community inhabiting this unexplored extreme environment and increase our understanding of high-altitude volcanic microorganism network.

2. Materials and Methods

2.1. Study Area and Sampling

The study plots were located along the crater margin of Tianchi Volcano, which is part of the Changbai Mountain Nature Reserve (41°41′49′′–42°25′18′′ N, 127°42′55′′–128°16′48′′ E) in Jilin Province, northeast China [32]. This Nature Reserve, at an altitude ranging from 500 m to 2691 m a.s.l., is characterized by a typical mountain climate, with low temperature and heavy precipitation, being the mean annual temperature 4.9–7.3 °C, and the mean annual precipitation over 800 mm [33]. The Changbai Mountain high altitude zone environment is constituted by a tundra belt, which is distributed between 1950 m and 2650 m a.s.l., and includes volcanic, glacial and periglacial landforms. The specific mean annual temperature (−4.8 °C) and precipitation (1154 mm) of this area are typical of tundra-periglacial climate. The tundra belt is covered with snow for about 6 months every year, from mid-October to mid-May [26]. The Tianchi Volcano crater area, around Tianchi lake, belongs to the tundra belt.
During 26–27 October 2019, soil stone mixture samples from northern and western crater area of Tianchi Volcano were collected. Six plots of 20 m diam were randomly selected at the very top of the mountain, along the crater margin, on each of the two analyzed sides of the crater to collect a total of 12 soil stone mixture samples. All plots were completely devoid of vegetation. The number of plots was chosen based on the accessibility of the crater site and following the sampling restrictions in the protected area. Each sample was collected from the center of the plot at shallow depth (1–5 cm), after removing the surface layer (ca. 1 cm). Figure 1 shows the geographical distribution of sampling locations, while specific information of coordinates and altitudes are listed in Table 1. The samples were immediately stored on ice in insulated containers. After returning to the laboratory, each sample was divided into two subsamples, one stored at 4 °C for isolation of culturable microbes, the other one stored at −80 °C for metabarcoding analysis. The pH of soil samples was measured adding distilled water to ground material at a ratio of 1:2.5 (w/v).

2.2. Cultivable Microbe Isolation, Microscopy and Molecular Analyses

Microbes were isolated following a modified dilution plate method [35]. Soil stone mixture suspensions were prepared as follows: samples were first evenly grounded in sterilized mortars, then 10 g of each sample were suspended in 40 mL sterile water to obtain the diluted suspension. 1 mL of uniform diluted solution was taken and further diluted in 9 mL sterile water. All sample diluted suspensions were finally shaken at 220 rpm for 5 min. 100 μL of suspension from each of the two dilutions were plated on Potato Dextrose Agar (PDA) amended with 100 mg/L penicillin and streptomycin to reduce the growth of bacteria to a minimum level and favor the isolation of fungi. Petri dishes were sealed and incubated at room temperature (25 °C), in the darkness. After microbial colonies appeared, the different morphotypes were accurately selected for strain isolation, based on characters such as texture and pigmentation. Single colonies were picked and transferred to new Petri dishes containing the same medium to obtain pure cultures.
Fungal morphological traits (such as branched septate hyphae, pseudohyphae, conidiophores, conidia, poroconidia, arthroconidia and sporangiospores) were examined under Nikon ECLIPSE Ci microscope for identification of isolates following the taxonomic keys for different taxa described by Gilman [36], Alexopoulos and Mims [37].
The modified CTAB (cetyl trimethyl ammonium bromide) method was used to extract fungal DNA from mycelia of pure cultures [38]. A total of 500 mg of fungal mycelia was added to a 2.0 mL microcentrifuge tube with 0.2 g zirconia beads, then 100 μL of sterile lysis buffer (1 M Tris-HCl (pH 8.0) 10 mL, 0.5 M EDTA (pH 8.0) 4 mL, 5 M NaCl 28 mL, CTAB 2 g, add ddH2O to 100 mL) was added. The mixture was homogenized for 30 s, 60 Hz, 2 cycles using a high-throughput tissue grinder (SCIENTZ-48, SCIENTZ, Ningbo, China). 400 μL lysis buffer was added after homogenization, then the mixture was incubated at 60 °C for 30–60 min in a water basin by inverting the tubes at intervals of 15 min. The tubes were added with 500 μL of SEVAG (chloroform:isoamyl alcohol, 24:1), mixed evenly, then centrifuged at 12,000 rpm for 15 min. The supernatant was transferred to a new clean 2.0 mL tube, and an equal volume of SEVAG was added. The resulting solution was mixed evenly and centrifuged at 12,000 rpm for 15 min. The last step was repeated twice, then the supernatant (225–300 μL) was taken and transferred to a new clean 1.5 mL microcentrifuge tube, 600 μL isopropyl alcohol were added and the tubes inverted gently several times, thus making DNA “ropes” precipitate visible. The tubes were put at −20 °C for 1 h and then centrifuged at 12,000 rpm for 5 min. Supernatant was removed carefully from each tube and the pellets were washed with 400 μL of 70% ethanol, then centrifuged at 12,000 rpm for 3 min. After repeating the last step twice, the aqueous phase was discarded and the tubes placed in super clean bench or fume cupboard to dry. The DNA was redissolved in 50 μL of sterile ddH2O for further study. For bacteria, DNA was obtained by boiling method [39]. Bacterial samples were put into 100 μL sterile ddH2O in 1.5 mL microcentrifuge tubes and heated at 100 °C boiling water for 15 min, DNA was dissolved in the liquid phase. This step was followed by centrifugation (14,000 rpm/1 min) to recover supernatant, then the tubes containing supernatant were stored at −20 °C.
The internal transcribed spacer (ITS) region of fungal rDNA was amplified using primers ITS1 and ITS4, whereas the V3 to V4 hypervariable region of bacterial 16S rRNA was amplified using the primers U341F and U806R (Table 2). PCR amplifications were performed in 25 mL reaction mixtures that were prepared with 1 μL of genomic DNA as template DNA, 12.5 mL of 2× Vazyme Rapid Taq Mater Mix, 1 μL of each primer (10 μM), and double-distilled H2O to a total volume of 25 μL, using the following reaction conditions: Temperature of 95 °C for 3 min, 30 cycles at 95 °C for 15 s, 55 °C for 15 s, 72 °C for 15 s, and a final extension at 72 °C for 5 min. PCR products were Sanger sequenced at Tianjin Tsingke Biological Technology Co., Ltd. (Tianjin, China). Sequences were assembled using BioEdit [40]. All obtained sequences were BLASTn searched in NCBI and assigned to potential genera and species. The nomenclature followed Index Fungorum (indexfungorum.org) [41]. Sequences were deposited in GenBank under the accession numbers MW582315–MW582428 (fungi) and MW577449–MW577451 (bacteria).
All fungal colonies isolated in this study were inoculated, both on slanted medium in glass tubes and in strain preservation tubes containing double-sterilized ultra-pure water, and stored at 4 °C for long term preservation. Bacterial strains were preserved in 30% glycerol at –80 °C. All strains were deposited in the LP Culture Collection (personal culture collection held in the laboratory of Prof. Lorenzo Pecoraro), at the School of Pharmaceutical Science and Technology, Tianjin University, Tianjin, China.

2.3. Assessment of Microbial Community Using Illumina Sequencing

The total genomic DNA was extracted from 0.5 g of soil stone mixture for each sample by using a FastDNA® Spin Kit for soil (MP Biomedicals, Solon, OH, USA) according to the manufacturer’s instructions. The products were checked in 1% agarose gel and quantified with NanoDrop 2000 UV-vis spectrophotometer (Thermo Scientific, Waltham, MA, USA). The fungal ITS2 rDNA region was amplified using the fungus-specific primer pair ITS3F (5′-GCATCGATGAAGAACGCAGC-3′)–ITS4R (5′-TCCTCCGCTTATTGATATGC-3′). The V3–V4 hypervariable region of the bacterial 16S rRNA gene was amplified using the primers 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTA AT-3′). PCR reactions were carried out in triplicate 20 μL mixture containing 4 μL of 5× TransStart FastPfu Buffer, 2 μL of dNTPs (2.5 mM), 0.8 μL of each primer (5 μM), 0.4 μL of TransStart FastPfu Polymerase, 0.2 μL of bovine serum albumin (Takara Bio Inc., Kusatsu, Shiga, Japan) and 10 ng of template DNA. The PCR reaction conditions were as follows: an initial denaturation of 95 °C for 3 min, followed by 35 cycles (for the amplification of fungal ITS region, for bacteria was 27 cycles) of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s and extension at 72 °C for 45 s, and a final extension at 72 °C for 10 min. PCR products were checked in 2% agarose gel and purified with the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using QuantiFluor-ST Fluorometer (Promega, Madison, WI, USA). Amplicon libraries were generated using the NEXTFLEX® Rapid DNA-Seq Kit (Bioo Scientific, Austin, TX, USA) following the manufacturer’s recommendations to add index codes. The library quality was assessed by using the QuantiFluorTM-ST Blue fluorescence quantitative system (Promega) according to the manual. Equimolar amounts of purified amplicons from different samples were pooled and paired-end sequenced (2 × 300 bases) on an Illumina MiSeq PE300 platform (Illumina Inc., San Diego, CA, USA) at Shanghai Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). The MiSeq raw sequencing fastq files were deposited in the Sequence Read Archive (SRA) at the National Center for Biotechnology Information (NCBI) as BioProject ID PRJNA700061.

2.4. Bioinformatics of Fungal and Bacterial Sequences

Raw sequences were quality-filtered by Fastp v0.19.6 [42] according to the following rules: (1) Filter the bases in the tail of the reads with quality score below 20; (2) set 50 bp sliding windows on the reads, cut off the back-end bases from the window if the average quality value in the window is lower than 20, and the reads below 50 bp after quality control were removed; (3) reads containing ambiguous nucleotide (N) were eliminated; (4) the barcode mismatches were not allowed and the maximum number of primer mismatches was 2. Paired-end reads were subsequently merged by FLASH v1.2.11 [43] with a minimum overlap of 10 bp and maximum allowable mismatch ratio 0.2. Quality-controlled sequences were clustered in operational taxonomic units (OTUs) at a 97% similarity threshold using UPARSE v7.1 [44]. Chimeric sequences were identified and removed using UCHIME [45]. The sequence with the highest abundance were selected as representative sequence for each OTU. Taxonomic assignments of the representative sequence of each OTU were based on RDP Classifier v2.2 [46] using the UNITE Database v8.0 (for fungal sequences) [47], and SILVA 16S rRNA Database (release 138, for bacterial sequences) [48,49]. Volcanic fungal and bacterial OTU table were rarefied to counts up to 43,002 and 36,117 reads per sample, respectively, (these were the lowest sequence depths obtained from all samples) and used for downstream analyses, including calculation of relative abundance of different taxonomic groups, α-diversity and β-diversity analyses, as well as network analyses.

2.5. Statistical Analysis

α-diversity indexes, including richness (number of OTUs), Shannon, Chao1, and Good’s Coverage were calculated in Mothur used to estimate the richness, diversity and coverage of microbial communities. Rarefaction curves [50] and Rank abundance curve [51] were drawn in R package for evaluation of sufficient sequence depth, as well as assessment of community richness and evenness. Between-groups Venn diagram was plotted using R to identify unique and common OTUs. Principal co-ordinate analysis (PCoA) and Analysis of Similarity (ANOSIM) with 999 permutations were performed for the assessment of dissimilarities between samples based on taxonomic Bray–Curtis [52], phylogenetic Weighted and Unweighted Unifrac distance matrices [53].

2.6. Network Analysis

The interactions between microbial taxa were analyzed through a network structure to assess the complexity of fungal and bacterial communities in Tianchi crater. Rarefied fungal and bacterial OTU tables were used for network analyses of crater soils. Only the OTUs in the top 100 abundance level were retained for the analysis to reduce the complexity of the data sets. Spearman’s rank correlation coefficients (r) between the OTUs with a magnitude ≥0.8 or ≤−0.8 and statistically significant (p < 0.05) were further included for network construction. Gephi (version 0.9.2) was used for the topological properties estimation, visualization and modular analysis of the network [54]. The nodes in the constructed network represent the OTUs indicating different taxa, whereas the edges correspond to significant positive or negative correlations between taxa.

3. Results

3.1. Diversity of Culturable Microbes

In total, 117 fungal strains were isolated from the 12 soil samples collected along the margin of Tianchi crater (Table S1). The isolated fungi belonged to 51 species in 31 genera, based on combined ITS sequencing and microscopic identification (Table 3). Among them, 92.16% (i.e., 47 species, 93 strains) belonged to 27 genera of Ascomycota, 5.88% (3 species, 6 strains) belonged to three genera of Basidiomycota, while Mucoromycota was represented by a single species, Mucor hiemalis, which yielded 18 strains (Figure S1, Table S2). The most common genera included Penicillium (9.80%), Trichoderma (9.80%), Cladosporium (7.84%), Didymella (7.84%), Alternaria (5.88%) and Fusarium (5.88%) (Figure S2, Table S2). The most common species was Mucor hiemalis (15.38%), followed by Fusarium tricinctum (9.40%) and Cladosporium cladosporioides (6.84%). Three resistant bacterial strains were also isolated in the media supplemented with antibiotics, and were identified as Rhodococcus degradans, R. qingshengii and Chryseobacterium lactis (Table S1).

3.2. Microbial Diversity by Illumina Sequencing

From metabarcoding analyses, a total of 4,371,622 raw reads with 772,186 effective sequences were obtained by fungal ITS2 gene sequencing, using the Illumina MiSeq PE300 platform. The number of sequences from the samples ranged from 43,002 to 70,583. After normalization of sequences by resampling 43,002 reads every sample, 1254 OTUs could be classified into the fungal kingdom (Table S3). From bacterial 16S V3–V4 region sequencing, 3,105,944 raw reads with 507,467 effective sequences were obtained. 2988 bacterial OTUs were clustered subsequently after the OTU table was rarefied to counts up to 36,117 reads per sample (the number of sequences ranged from 36,117 to 49,230) (Table S4).
OTU richness, Shannon and Chao1 indexes were calculated for both fungal and bacterial communities (Table 4). Tianchi crater margin soil contained more diverse and richer bacterial than fungal communities as shown by the significantly higher α-diversity indexes measured according to Mean ± SD values (Table 4). Good’s coverage for all samples was > 99.0% (Table 4), which indicated the small probability that sequences were not detected in each sample. This result was consistent with the rarefaction curves related to the OTUs that showed all the curves approaching saturation plateau (Figure 2). The rank abundance curve representing richness and evenness of fungal and bacterial communities showed that the highest richness and evenness was found from northern soil samples, for fungi, while from western samples, for bacteria (Figure 3).
From the rarefied fungal OTUs table, 1077 out of the 1254 retrieved OTUs were identified into eight phyla, 38 classes, 112 orders, 250 families and 455 genera (Table S3), whereas 165 OTUs remained unidentified and 12 Incertae sedis at phylum level (14.11%). Among the 1077 identified fungal OTUs, 804 were assigned to Ascomycota, 213 to Basidiomycota, 32 Chytridiomycota and 16 Mortierellomycota. Other phyla, including Monoblepharomycota, Glomeromycota, Rozellomycota and Mucoromycota were recovered in small proportions (12 OTUs) (Table S3). For the bacterial community, except 52 OTUs that remained unidentified at phylum level (1.74%), 2936 OTUs were assigned to 32 phyla, 91 classes, 217 orders, 353 families and 654 genera (Table S4). The phyla Proteobacteria (557OTUs), Chloroflexi (375 OTUs), Actinobacteriota (334 OTUs), Bacteroidota (314 OTUs), Acidobacteriota (217 OTUs), Myxococcota (184 OTUs), Patescibacteria (158 OTUs), Bdellovibrionota (151 OTUs), Gemmatimonadota (112 OTUs) and Planctomycetota (106 OTUs) were the most dominant in the crater soil (Table S4).
In terms of the relative abundance, the class Leotiomycetes (49.13%), followed by Dothieomycetes (12.90%), dominated the Tianchi crater soil fungal community (Figure 4A), while for the bacterial community, Gammaproteobacteria (26.30%), Actinobacteria (26.04%) and Alphaproteobacteria (17.29%) were dominant (Figure 4B). At genus level, the relative abundance of Leucosporidium was found to be the highest (2.30%) among Basidiomycota, while the most abundant ascomycetes genera remained not clearly discriminated as belonging to the class Leotiomycetes (Figure 5A). The genera from other phyla were with very low relative abundances (<1%) and merged as “others” in the bar plot (Figure 5A). Among bacteria, the genus Sphingomonas (5.50%) was the most abundant in the analyzed crater soil (Figure 5B).
From a comparison of OTUs diversity between northern and western crater soil samples, a total of 566 fungal and 1816 bacterial OTUs, accounting for 45.14% and 60.78% of the total OTUs, respectively, were found to be common to both sides of the crater (Figure 6). Some of the common OTUs exhibited a relatively high abundance among all OTUs. For instance, fungal OTU1055 (unclassified Leotiomycetes) and OTU1230 (unclassified Helotiales) were the most abundant OTUs detected, showing a relative abundance of 11.89% and 11.87%, respectively. Similarly, common bacterial OTU11 (unclassified Erwiniaceae), with relative abundance 5.07%, and OTU541 (unclassified Pseudarthrobacter, 4.62%) dominated the bacterial community.
The variations of fungal and bacterial OTUs were analyzed by principal co-ordinate analysis (PCoA). Fungal communities in northern and western crater soil were clearly distinguished based on Bray–Curtis distance (ANOSIM p = 0.011) (Figure 7A), weighted Unifrac distance (ANOSIM p = 0.035) and unweighted Unifrac distance (ANOSIM p = 0.003) (Figure 7B,C). The separation of communities was more significant using unweighted Unifrac distance, which is more sensitive to taxa presence or absence regardless of their abundances compared with Bray–Curtis distance and weighted Unifrac distance. For bacterial communities, on the contrary, no clear separation was found based on the different dissimilarity matrices (Figure 8).

3.3. Network Complexity for Bacterial and Fungal Communities

Separated co-occurrence network analyses were performed for both fungal and bacterial communities to explore and evaluate the complex interactions between the microbial taxa detected in the analyzed crater soils. All interactions comprised in the network were strongly correlated and statistically significant (r ≥ 0.8 or ≤ −0.8 and p < 0.05). In general, there were considerable differences in network topology and structure between fungal and bacterial communities (Figure 9, Table 5). The more compact and complicated bacterial network than the fungal network was reflected by the average path lengths (2.661 vs. 4.982) and average degrees (10.143 vs. 3.487) (Figure 9, Table 5). Contrasting correlations between OTUs within fungal and bacterial communities were observed. Seventy-seven nodes were included in the network of the fungal community, with 136 edges, 108 positive interactions and 28 negative interactions, while 84 nodes were included in the bacterial network. The network metrics of fungal community were significantly lower than those for bacterial community (426 edges, 302 positive and 124 negative interactions) (Figure 9, Table 5). Both networks were predominated by positive interactions (79.41% and 70.89%), indicating that the mutual effects dominate the microbial communities in the crater. For fungal community, Leucosporidium sp. (OTU348), Paraphoma fimeti (OTU398), Alatospora sp. (OTU952), unclassified Chaetothyriales (OTU955) and unclassified Eurotiomycetes (OTU987) were found to have the most abundant interactions with other nodes (degree: 9) (Figure 9A). Pedobacter sp. (OTU2266) was found to have the most abundant interactions in the network of bacterial community (degree = 26) (Figure 9B).

4. Discussion

In this study, we provided a comprehensive characterization of fungal and bacterial communities along the crater margin of Tianchi Volcano, in Changbai Mountain Nature Reserve, by combing molecular and morphological analyses of culturable microbes, and metabarcoding analyses based on Illumina sequencing. A few previous studies have been conducted on microbial diversity in relatively lower altitudes of Changbai Mountain northern slope, mostly using high-throughput sequencing methods alone [14,33,55]. To the best of our knowledge, this study represents the first microbial community analysis in different sides of Tianchi Volcano crater margin, using integrated culture-dependent and metabarcoding approaches.
Many microbial species isolated from this study have been frequently reported from high-altitude or cold environments, and have been extensively applied in biotechnological fields, showing particular research significance due to their peculiar ecology. For instance, Aureobasidium strains isolated from Tianchi northern crater sediments, highly matched with A. pullulans isolated from saxicolous lichens growing in the north side of Taibai Mountain in China, at the altitude of 2614 m a.s.l. [56], and A. pullulans found from the glacier surface snow of Tibetan Plateau in China (>99% similarity) by Shao [57]. Besides, A. pullulans have been previously isolated from high Italian glaciers [7], high Arctic glaciers [58] and Antarctica soils [59]. It has been well documented that A. pullulans is a psychrophilic fungus able to produce various useful bioproducts (enzymes, pullulan, single-cell proteins and siderophore) for waste treatments, chemical industry materials, food industries and biocontrol [59]. Our finding of Candida tropicalis in the western crater of Tianchi Volcano is in agreement with previous studies in which Candida is regarded as a common Antarctic yeast genus [8]. Candida tropicalis has been described capable of degrading high-concentration phenol [60], and tolerant to copper [61]. This species may play an important functional role in the studied crater area of Tianchi volcano characterized by high heavy-metal concentration. The ascomycete Epicoccum nigrum, isolated from both sides of the analyzed crater, has been described as a xerotolerant species in a study conducted in Annapurna Mountains of Nepal [62]. This fungal species has been found to produce antibacterial compounds [3] and described as a biocontrol agent [63]. Among the few bacterial isolates belonging to Chryseobacterium and Rhodococcus genera, Chryseobacterium strain showed high similarity with C. lactis (accession number = MT065804) isolated from Qinling high altitude area in China by Men X. in an unpublished study, while Rhodococcus strains highly matched with R. degradans (accession number: LR216744) isolated from sodic soil in Hungary by Krett G., and R. qingshengii (accession number: MT632489) isolated from oil well in Russia by Borzenkov et al., in unpublished studies. Species from Rhodococcus have been previously found to be psychrotolerant in a study conducted in Antarctica [64]. This considerable number of microbes isolated from Tianchi Volcano, which have previously shown preference for high-altitude, cold, and geothermal extreme environments, and have been described as psychrotolerant and metal-tolerant microorganisms, provide new insights in the characterization of the ecological features of the investigated environment, and may represent a precious source for the isolation of new bioactive compounds [65,66].
The use of molecular and high-sensitive metabarcoding methods allowed the detection and identification of both common and rare microbial taxa, thus providing detailed and accurate information on fungal and bacterial communities in the investigated high-altitude volcanic habitat. The high relative abundances of unclassified genera detected in this study may indicate the presence of a significant number of undiscovered and possibly endemic taxa in the analyzed communities. Recent studies have described the fungal diversity in apparently barren (plant-free) zones of high alpine environments, including the Himalayas (Nepal, 5146 and 5509 m a.s.l.) [67], Colorado (USA, 3660–3800 m a.s.l.) [68,69], Llullaillaco (6030 and 6330 m a.s.l.) [8] and Socompa (5235 m a.s.l.) volcanoes [70] in the Andes, along the Chilean-Argentinian border. A major group of basidiomycetes found from these barren high-altitude soils (mostly Colorado) was constituted by members of the Microbotryomycetes, including Leucosporidium antarcticum. Accordingly, our results revealed that Leucosporidium was the most abundant genus in Basidiomycota from high-throughput sequencing, and two sequences obtained from isolated strains showed 100% similarity with L. creatinivorum isolated from the soil of King George Island, in the sub-Antarctic region. Leucosporidium creatinivorum is known as a psychrotolerant yeast, which has been frequently isolated from cold extreme environments, such as glaciers in Argentina, Russia, Iceland and Italian Alps [71,72], and has been also reported from soil, marine sponges and lichens in Antarctica [73,74,75]. Members of Leucosporidium have been described as sources of cold-active enzymes [73,76,77] and anti-freeze proteins [78], and have been shown to possess phenol and phenol-related compounds degradation ability [79]. The high relative abundance of Leucosporidium in high-altitude volcanic soil stone mixtures of Tianchi volcano confirmed the peculiar preference of this genus for extreme cold environments. Up to now, to our knowledge, the existence of Leucosporidium in China was only detected by metagenomic analysis [80,81]. Our results represent the first isolation of L. creatinivorum in China. These isolated strains certainly deserve further attention as potential producers of metabolites for biotechnological application. Besides, in our study, Leucosporidium was also found to have the most complicated interactions in the fungal network, which may indicate a fundamental ecological role of this taxon in maintaining the stability of the fungal community in high elevation volcanic environments.
The bacterial community was dominated by Proteobacteria, result in agreement with previous studies where this taxon has been described as the most dominant soil bacterial phylotype in different continents across the globe [82]. Gammaproteobacteria was the most abundant class within Proteobacteria, followed by Alphaproteobacteria, which is consistent with previous findings from Deception Island Volcano in Antarctica [83]. Sphingomonas, the most abundant bacterial genus in Tianchi Volcano, has been previously reported from other high-altitude volcanic areas, including the active volcano Socompa in Argentina [84], and El Chichón volcano in Mexico [85]. Members of this genus play prominent functional roles in the remediation of contaminated environments. For instance, several Sphingomonas species have been extensively associated to the chelation of heavy metals, and degradation/bio-transformation of aromatic hydrocarbons [86]. Sphingomonas strains have been also found to produce highly beneficial phytohormones [86]. The dominant presence of Sphingomonas in high-altitude, metal-enriched, and carbon poor sediments of Tianchi crater, suggested an important contribution of this bacterial group to the functionality of the studied ecosystem, and confirmed the metal-tolerant ability of these microorganisms, which probably possess special nutritional strategies to adapt and survive in barren habitats, under extreme conditions. Further studies focusing on the metabolic activity and nutritional mode of Sphingomonas species may clarify the mechanisms allowing extremotolerant microbes to strive in environments which are unsuitable for the majority of life forms.
The significant clear variation of soil fungal communities in the two sides of the Tianchi crater suggested that geographical and environmental factors may have a stronger influence on the diversity of fungi than bacteria, which, instead, did not show any remarkable community difference between North and West sampling areas. Different community diversity patterns may be the result of mixed effects of multiple factors, such as solar radiation, temperature, and soil properties that may be more relevant to fungi than bacteria [87]. However, our analysis of samples pH showed highly homogeneous values, with average pH at 7.76, which could not provide useful information to explain the observed community diversity variation. Further studies are needed in order to analyze environmental factors affecting bacterial and fungal community diversity in high elevation volcanic extreme environments. The results of PCoA were consistent with co-occurrence network analyses. Compare to the fungal network, the bacterial network was composed by more complicated and closer associations, in terms of numbers of nodes, edges, average degree and average path length, thus indicating that the bacterial communities possess rapid responsiveness to perturbation [88], and much stronger resistance to environmental changes than fungal communities. However, these related environmental factors remain to be studied.
The detailed characterization of Tianchi Volcano microbial diversity provided in this study may constitute an important reference for future long-term monitoring, aiming at tracing the effects of global warming on this delicate environment. Indeed, previous studies have proved that increasing global average temperatures can drive successions of microbial communities, for instance by changing the quality and quantity of genes potentially available for horizontal gene transfer, and can lead to increasingly divergent succession, with possibly higher impact on fungi than bacteria [89,90]. Besides, the impact of climate changes is predicted to be very pronounced at high-elevations and in tundra ecosystems [91,92]. The area along the margin of the Tianchi crater represents a typical example of a high-altitude tundra ecosystem in China. It has been reported that the mean annual temperature in Changbai Mountain Nature Reserve increased at a rate of 0.392 °C/10 years from 1958 to 2015, which is significantly higher than the national warming rate (0.22 °C/10 years) [93]. Because global warming may particularly alter the biogeochemistry and ecology of cold soil ecosystems [94,95,96], the increasing warming registered in Changbai Mountain region deserves attention, and further studies are needed to verify the response to environmental changes of the particularly sensitive and fragile high-altitude tundra microbes unveiled from our study in Tianchi volcano.

5. Conclusions

Our results represent an unprecedent comprehensive microbial community analysis along the high-altitude crater margin of Tianchi Volcano, by combing culture-dependent and metabarcoding analyses. We observed that Tianchi Volcano hosted a combination of taxonomic groups characteristic of high-altitude, cold and geothermal environments, with a considerable number of isolated microbes being of particular research significance, due to their rarity and peculiar ecology. Our study suggested that the structure and diversity of fungal community was more sensitive to environmental and geographical changes compare to bacterial community in the analyzed area. Our findings may represent an important starting point for future studies to explore the precious metabolite resources of the isolated microbes, and elucidate the effect of different environmental factors on community structure and dynamics in high-altitude volcanic environments.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/life11040280/s1, Figure S1: Number of fungal genera, species and strains in different phyla obtained from Tianchi crater soil. Figure S2: Number of fungal species and strains from most abundant genera observed in Tianchi crater soil. Table S1: Fungal and bacterial diversity molecularly detected in Tianchi Volcano soil, from DNA extracted from isolated microbes. Table S2: An overview of fungal strains isolated from the crater soil samples. Table S3: Rarefied fungal OTUs table based on fungal ITS2 rDNA region Illumina sequencing. Table S4: Rarefied bacterial OTUs table based on V3–V4 hypervariable region of the bacterial 16S rRNA gene Illumina sequencing.

Author Contributions

X.W. and L.P. conceived the study; samples were collected by X.W. and L.P.; the experiments were designed and supervised by L.P.; laboratory experiments and analysis were performed by X.W.; results were analyzed by X.W. and L.P.; X.W. prepared the original draft while L.P. critically revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The Fungal and bacterial DNA sequences amplified during this study are available in GenBank under accessions MW582315–MW582428 (fungi) and MW577449–MW577451 (bacteria), and in the Sequences Read Archive of NCBI as BioProject ID PRJNA700061.

Acknowledgments

We are grateful to the Staff of Changbai Mountain Nature Reserve for the assistance during sampling.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rawat, N.; Joshi, G.K. Bacterial community structure analysis of a hot spring soil by next generation sequencing of ribosomal RNA. Genomics 2019, 111, 1053–1058. [Google Scholar] [CrossRef]
  2. Herrera, A.; Cockell, C.S. Exploring microbial diversity in volcanic environments: A review of methods in DNA extraction. J. Microbiol. Methods 2007, 70, 1–12. [Google Scholar] [CrossRef]
  3. Marie-Antoinette, B.; Gerard, D.; Robert, B.; Arlette, N. New antibiotics from the fungus Epicoccum nigrum. I. Fermentation, iso-lation and antibacterial properties. J. Antibiot. 1978, 31, 1099–1101. [Google Scholar]
  4. Ge, B.; Liu, B.; Nwet, T.T.; Zhao, W.; Shi, L.; Zhang, K. Bacillus methylotrophicus Strain NKG-1, Isolated from Changbai Mountain, China, Has Potential Applications as a Biofertilizer or Biocontrol Agent. PLoS ONE 2016, 11, e0166079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Tapia-Vázquez, I.; Sánchez-Cruz, R.; Arroyo-Domínguez, M.; Lira-Ruan, V.; Sánchez-Reyes, A.; Sánchez-Carbente, M.D.R.; Padilla-Chacón, D.; Batista-García, R.A.; Folch-Mallol, J.L. Isolation and characterization of psychrophilic and psychrotolerant plant-growth promoting microorganisms from a high-altitude volcano crater in Mexico. Microbiol. Res. 2020, 232, 126394. [Google Scholar] [CrossRef] [PubMed]
  6. Prins, R. Adaptation of microorganisms to extreme environments. FEMS Microbiol. Lett. 1990, 75, 103–104. [Google Scholar] [CrossRef]
  7. Turchetti, B.; Buzzini, P.; Goretti, M.; Branda, E.; Diolaiuti, G.; D’Agata, C.; Smiraglia, C.; Vaughan-Martini, A. Psychrophilic yeasts in glacial environments of Alpine glaciers. FEMS Microbiol. Ecol. 2008, 63, 73–83. [Google Scholar] [CrossRef] [Green Version]
  8. Schmidt, S.; Naff, C.; Lynch, R. Fungal communities at the edge: Ecological lessons from high alpine fungi. Fungal Ecol. 2012, 5, 443–452. [Google Scholar] [CrossRef]
  9. Alexander, M. Introduction to Soil Microbiology, 2nd ed.; John Wiley and Sons: New York, NY, USA, 1977; pp. 1–467. ISBN 978-00-6054-646-5. [Google Scholar]
  10. Amann, R.I.; Ludwig, W.; Schleifer, K.H. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol. Rev. 1995, 59, 143–169. [Google Scholar] [CrossRef] [PubMed]
  11. Mirete, S.; Morgante, V.; González-Pastor, J.E. Functional metagenomics of extreme environments. Curr. Opin. Biotechnol. 2016, 38, 143–149. [Google Scholar] [CrossRef] [PubMed]
  12. Xu, W.; Gao, Y.-H.; Gong, L.-F.; Li, M.; Pang, K.-L.; Luo, Z.-H. Fungal diversity in the deep-sea hadal sediments of the Yap Trench by cultivation and high throughput sequencing methods based on ITS rRNA gene. Deep. Sea Res. Part I Oceanogr. Res. Pap. 2019, 145, 125–136. [Google Scholar] [CrossRef]
  13. Rao, M.P.N.; Dong, Z.-Y.; Luo, Z.-H.; Li, M.-M.; Liu, B.-B.; Guo, S.-X.; Hozzein, W.N.; Xiao, M.; Li, W.-J. Physicochemical and Microbial Diversity Analyses of Indian Hot Springs. Front. Microbiol. 2021, 12, 627200. [Google Scholar] [CrossRef]
  14. Han, D.; Wang, N.; Sun, X.; Hu, Y.; Feng, F. Biogeographical distribution of bacterial communities in Changbai Mountain, Northeast China. Microbiology 2018, 7, e00529. [Google Scholar] [CrossRef] [Green Version]
  15. Liu, Y.G.; Zhang, M.L.; Guan, D.X.; Song, X.F.; Yuan, F.H.; Wu, J.B.; Wang, A.Z. Natural science research in Changbai Mountain during 1956–2018: A review. J. Appl. Ecol. 2019, 30, 1783–1796. [Google Scholar]
  16. Yang, H.; Lü, G.; Jiang, H.; Shi, D.-N.; Liu, Z. Diversity and distribution of soil micro-fungi along an elevation gradient on the north slope of Changbai Mountain. J. For. Res. 2017, 28, 831–839. [Google Scholar] [CrossRef]
  17. Cai, J.; Wang, Y.; Liu, D.; Zeng, Y.; Xue, Y.; Ma, Y.; Feng, Y. Fervidobacterium changbaicum sp. nov., a novel thermophilic anaerobic bacterium isolated from a hot spring of the Changbai Mountains, China. Int. J. Syst. Evol. Microbiol. 2007, 57, 2333–2336. [Google Scholar] [CrossRef] [PubMed]
  18. Zhou, L.-W.; Qin, W.-M. Sistotrema subconfluens sp. nov. (Cantharellales, Basidiomycota) from Changbaishan Nature Reserve, northeastern China. Mycoscience 2013, 54, 178–182. [Google Scholar] [CrossRef]
  19. He, X.-L.; Wang, D.; Peng, W.-H.; Gan, B.-C. Two new Entoloma s.L. species with serrulatum-type lamellar edge from Changbai Mountains, Northeast China. Mycol. Prog. 2017, 16, 761–768. [Google Scholar] [CrossRef]
  20. Wei, H.; Wang, Y.; Jin, J.; Gao, L.; Yun, S.-H.; Jin, B. Timescale and evolution of the intracontinental Tianchi volcanic shield and ignimbrite-forming eruption, Changbaishan, Northeast China. Lithos 2007, 96, 315–324. [Google Scholar] [CrossRef]
  21. Liu, J.-Q.; Chen, S.-S.; Guo, Z.-F.; Guo, W.-F.; He, H.-Y.; You, H.-T.; Kim, H.-M.; Sung, G.-H.; Kim, H. Geological background and geodynamic mechanism of Mt. Changbai volcanoes on the China–Korea border. Lithos 2015, 236–237, 46–73. [Google Scholar] [CrossRef]
  22. Yan, B.; Qiu, S.; Liu, Z.; ChangLai, X. Characteristics of the geothermal water in Changbai Mountain volcanic region, northeast of China. Arab. J. Geosci. 2017, 10, 261. [Google Scholar] [CrossRef]
  23. Sun, S. Volcanic Geology and Secondary Volcanic Hazards in the Changbaishan-Tianchi Volcano and its Adjacent Areas. Ph.D. Thesis, Jilin University, Changchun, Jilin, China, 2018. [Google Scholar]
  24. Yang, Q.F.; Yuan, X.J.; Wu, C.Z.; Pan, X.D.; Zhang, Y. The multibeam sounding exploration of the Tianchi caldera lakebed to-pography at the China-North Korea border. Acta Petrol. Sin. 2018, 34, 185–193. [Google Scholar]
  25. Andreevaa, O.A.; Yarmolyuka, V.V.; Andreevaa, I.A.; Ji, J.Q.; Li, W.R. The Composition and Sources of Magmas of Chang-baishan Tianchi Volcano (China-North Korea). Dokl. Earth Sci. 2014, 456, 572–578. [Google Scholar] [CrossRef]
  26. Shen, C.; Ni, Y.; Liang, W.; Wang, J.; Chu, H. Distinct soil bacterial communities along a small-scale elevational gradient in alpine tundra. Front. Microbiol. 2015, 6, 582. [Google Scholar] [CrossRef] [Green Version]
  27. Ni, Y.; Yang, T.; Zhang, K.; Shen, C.; Chu, H. Fungal Communities Along a Small-Scale Elevational Gradient in an Alpine Tundra Are Determined by Soil Carbon Nitrogen Ratios. Front. Microbiol. 2018, 9, 1815. [Google Scholar] [CrossRef] [PubMed]
  28. Shen, C.; Xiong, J.; Zhang, H.; Feng, Y.; Lin, X.; Li, X.; Liang, W.; Chu, H. Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem. 2013, 57, 204–211. [Google Scholar] [CrossRef]
  29. Shen, C.C.; Shi, Y.; Ni, Y.Y.; Deng, Y.; Van Nostrand, J.D.; He, Z.; Zhou, J.; Chu, H. Dramatic Increases of Soil Microbial Func-tional Gene Diversity at the Treeline Ecotone of Changbai Mountain. Front. Microbiol. 2016, 7, 1184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Wang, X.; Li, L.; Zhao, W.; Zhao, J.; Chen, X. Rhododendron aureum Georgi formed a special soil microbial community and competed with above-ground plants on the tundra of the Changbai Mountain, China. Ecol. Evol. 2017, 7, 7503–7514. [Google Scholar] [CrossRef]
  31. Tang, M.; Li, L.; Wang, X.; You, J.; Li, J.; Chen, X. Elevational is the main factor controlling the soil microbial community structure in alpine tundra of the Changbai Mountain. Sci. Rep. 2020, 10, 1–15. [Google Scholar] [CrossRef]
  32. Zhao, W.; Qi, X.; Lyu, J.; Yu, Z.; Chen, X. Characterization of microbial community structure in rhizosphere soils of Cowskin Azalea (Rhododendron aureum Georgi) on northern slope of Changbai Mountains, China. Chin. Geogr. Sci. 2015, 26, 78–89. [Google Scholar] [CrossRef]
  33. Ping, Y.; Han, D.; Wang, N.; Hu, Y.; Mu, L.; Feng, F. Vertical zonation of soil fungal community structure in a Korean pine forest on Changbai Mountain, China. World J. Microbiol. Biotechnol. 2016, 33, 12. [Google Scholar] [CrossRef] [PubMed]
  34. ArcGIS Online. Available online: https://maps.arcgis.com/index.html (accessed on 1 January 2021).
  35. Zhang, Y.; Liu, F.; Wu, W.; Cai, L. A phylogenetic assessment and taxonomic revision of the thermotolerant hyphomycete genera Acrophialophora and Taifanglania. Mycologia 2015, 107, 768–779. [Google Scholar] [CrossRef] [PubMed]
  36. Gilman, J.C. A Manual of Soil Fungi; The Iowa State University Press: Ames, IA, USA, 1957; pp. 1–450. ISBN 978-70-0001-234-5. [Google Scholar]
  37. Alexopoulos, C.J.; Mims, C.W. Introductory Mycology, 3rd ed.; John Wiley & Sons: New York, NY, USA, 1979; pp. 1–632. ISBN 978-71-2456-557-2. [Google Scholar]
  38. Lee, S.B.; Milgroom, M.G.; Taylor, J.W. A rapid, high yield mini-prep method for isolation of total genomic DNA from fungi. Fungal Genet. Rep. 1988, 35, 23. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, R.-Q.; Li, G.-X.; Li, X.-N.; Shen, X.-X.; Gao, Y.; Wang, L.; Fan, T.; Duan, Q.-X.; Wang, Y.-K.; Wang, J.; et al. A rapid and sensitive recombinase aided amplification assay incorporating competitive internal control to detect Bordetella pertussis using the DNA obtained by boiling. Int. J. Infect. Dis. 2019, 86, 108–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Hall, T. Bioedit: A user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95e98. [Google Scholar]
  41. Index Fungorum. Available online: indexfungorum.org (accessed on 15 September 2020).
  42. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  43. Magoč, T.; Magoč, M.; Salzberg, S.L. FLASH: Fast Length Adjustment of Short Reads to Improve Genome Assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef] [PubMed]
  44. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef] [PubMed]
  45. Edgar, R.C.; Haas, B.J.; Clemente, J.C.; Quince, C.; Knight, R. UCHIME Improves Sensitivity and Speed of Chimera Detection. Bioinformatics 2011, 27, 2194–2200. [Google Scholar] [CrossRef] [Green Version]
  46. Cole, J.R.; Wang, Q.; Fish, J.A.; Chai, B.; McGarrell, D.M.; Sun, Y.; Brown, C.T.; Porras-Alfaro, A.; Kuske, C.R.; Tiedje, J.M. Ribosomal Database Project: Data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014, 42, D633–D642. [Google Scholar] [CrossRef] [Green Version]
  47. Nilsson, R.H.; Larsson, K.-H.; Taylor, A.F.S.; Bengtsson-Palme, J.; Jeppesen, T.S.; Schigel, D.; Kennedy, P.; Picard, K.; Glöckner, F.O.; Tedersoo, L.; et al. The UNITE database for molecular identification of fungi: Handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 2019, 47, D259–D264. [Google Scholar] [CrossRef]
  48. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Opens external link in new window. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  49. Glöckner, F.O.; Yilmaz, P.; Quast, C.; Gerken, J.; Beccati, A.; Ciuprina, A.; Bruns, G.; Yarza, P.; Peplies, J.; Westram, R.; et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J. Biotechnol. 2017, 261, 169–176. [Google Scholar] [CrossRef] [PubMed]
  50. Gotelli, N.J.; Colwell, R.K. Quantifying biodiversity: Procedures and pitfalls in the measurement and comparison of species richness. Ecol. Lett. 2001, 4, 379–391. [Google Scholar] [CrossRef] [Green Version]
  51. Izsák, J.; Pavoine, S. Links between the species abundance distribution and the shape of the corresponding rank abundance curve. Ecol. Indic. 2012, 14, 1–6. [Google Scholar] [CrossRef]
  52. Sørensen, T. A method of establishing groups of equal amplitude in plant sociology based on similarity of species content and its application to analyses of the vegetation on Danish commons. Biol. Skr. 1948, 5, 1–34. [Google Scholar]
  53. Lozupone, C.; Knight, R. UniFrac: A New Phylogenetic Method for Comparing Microbial Communities. Appl. Environ. Microbiol. 2005, 71, 8228–8235. [Google Scholar] [CrossRef] [Green Version]
  54. Bastian, M.; Heymann, S.; Jacomy, M. Gephi: An Open Source Software for Exploring and Manipulating Networks. In Proceedings of the Third International Conference on Weblogs and Social Media (ICWSM 2009), San Jose, CA, USA, 17–20 May 2009; pp. 361–362. [Google Scholar]
  55. Liu, M.; Sui, X.; Hu, Y.; Feng, F. Microbial community structure and the relationship with soil carbon and nitrogen in an original Korean pine forest of Changbai Mountain, China. BMC Microbiol. 2019, 19, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Wang, D.S. Study of Microbial Flora and Antagonistic Actinomycete Resources in Five Habitats on North Side of Taibai Mountain, the Main Peak of Qin Mountains. Ph.D. Thesis, Northwest A&F University, Yangling, Shanxi, China, 2014. [Google Scholar]
  57. Shao, S.Y. Diversity of Culturable Yeast in the Snowpack from Different Glaciers on Tibetan Plateau and its Adjacent Regions. Master’s Thesis, Lanzhou University, Lanzhou, Gansu, China, 2012. [Google Scholar]
  58. Gostinčar, C.; Uršič, V.; De Hoog, S.; Gunde-Cimerman, N. Local evolution of black yeast A. pullulans in subglacial Arctic ice. In Proceedings of the International Conference Alpine Polar Microbiology, Innsbruck, Austria, 27–31 March 2006; p. 19. [Google Scholar]
  59. Chi, Z.M.; Wang, F.; Chi, Z. Bioproducts from Aureobasidium pullulans, a biotechnologically important yeast. Appl. Microbiol. Biotechnol. 2009, 82, 793–804. [Google Scholar] [CrossRef] [PubMed]
  60. Wang, L.; Lin, Y.; Yang, L.; Yu, P.; Xie, Z.; Luo, Y. Candida tropicalis: Characterization of a strain capable of degrading high concentrations of phenol. Biotechnol. Lett. 2011, 33, 943–946. [Google Scholar] [CrossRef] [PubMed]
  61. Rehman, A.; Farooq, H.; Hasnain, S. Biosorption of copper by yeast, Loddermyces elongisporus, isolated from industrial effluents: Its potential use in wastewater treatment. J. Basic Microbiol. 2008, 48, 195–201. [Google Scholar] [CrossRef]
  62. Petrovič, U.; Gunde-Cimerman, N.; Zalar, P. Xerotolerant mycobiota from high altitude Anapurna soil, Nepal. FEMS Microbiol. Lett. 2000, 182, 339–342. [Google Scholar] [CrossRef] [PubMed]
  63. Ogórek, R.; Plaskowska, E. Epicoccum nigrum for biocontrol agents in vitro of plant fungal pathogens. Commun. Agric. Appl. Biol. Sci. 2011, 76, 691–697. [Google Scholar]
  64. Silva, L.J.; Souza, D.T.; Genuario, D.B.; Hoyos, H.A.V.; Santos, S.N.; Rosa, L.H.; Zucchi, T.D.; Melo, I.S. Rhodococcus psychrotolerans sp. nov., isolated from rhizosphere of Deschampsia antarctica. Antonie van Leeuwenhoek 2018, 111, 629–636. [Google Scholar] [CrossRef] [PubMed]
  65. Niu, S.; Liu, D.; Hu, X.; Proksch, P.; Shao, Z.; Lin, W. Spiromastixones A-O, antibacterial chlorodepsidones from a deep-sea-derived Spiromastix sp. fungus. J. Nat. Prod. 2014, 77, 1021–1030. [Google Scholar] [CrossRef] [PubMed]
  66. Figueroa, L.; Jiménez, C.; Rodríguez, J.; Areche, C.; Chávez, R.; Henríquez, M.; de la Cruz, M.; Díaz, C.; Segade, Y.; Vaca, I. 3-Nitroasterric acid derivatives from an Antarctic sponge-derived Pseudogymnoascus sp. fungus. J. Nat. Prod. 2015, 78, 919–923. [Google Scholar] [CrossRef]
  67. King, A.J.; Karki, D.; Nagy, L.; Racoviteanu, A.; Schmidt, S.K. Microbial biomass and activity in high elevation (>5100 meters) soils from the Annapurna and Sagarmatha regions of the Nepalese Himalayas. Himal. J. Sci. 2011, 6, 11–18. [Google Scholar] [CrossRef] [Green Version]
  68. King, A.; Meyer, A.; Schmidt, S.K. High levels of microbial biomass and activity in unvegetated tropical and temperate alpine soils. Soil Biol. Biochem. 2008, 40, 2605–2610. [Google Scholar] [CrossRef]
  69. King, A.J.; Freeman, K.R.; McCormick, K.F.; Lynch, R.C.; Lozupone, C.; Knight, R.; Schmidt, S.K. Biogeography and habitat modelling of high-alpine bacteria. Nat. Commun. 2010, 1, 53. [Google Scholar] [CrossRef] [Green Version]
  70. Costello, E.K.; Halloy, S.R.P.; Reed, S.C.; Sowell, P.; Schmidt, S.K. Fumarole-Supported Islands of Biodiversity within a Hyperarid, High-Elevation Landscape on Socompa Volcano, Puna de Atacama, Andes. Appl. Environ. Microbiol. 2008, 75, 735–747. [Google Scholar] [CrossRef] [Green Version]
  71. De Garcã­a, V.; Brizzio, S.; Libkind, D.; Buzzini, P.; Van Broock, M. Biodiversity of cold-adapted yeasts from glacial meltwater rivers in Patagonia, Argentina. FEMS Microbiol. Ecol. 2007, 59, 331–341. [Google Scholar] [CrossRef]
  72. De García, V.; Coelho, M.A.; Maia, T.M.; Rosa, L.H.; Vaz, A.M.; Rosa, C.A.; Sampaio, J.P.; Gonçalves, P.; Van Broock, M.; Libkind, D. Sex in the cold: Taxonomic reorganization of psychrotolerant yeasts in the order Leucosporidiales. FEMS Yeast Res. 2015, 15, fov019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Carrasco, M.; Rozas, J.M.; Barahona, S.; Alcaíno, J.; Cifuentes, V.; Baeza, M. Diversity and extracellular enzymatic activities of yeasts isolated from King George Island, the sub-Antarctic region. BMC Microbiol. 2012, 12, 251. [Google Scholar] [CrossRef] [Green Version]
  74. Vaca, I.; Faúndez, C.; Maza, F.; Paillavil, B.; Hernández, V.; Acosta, F.; Levicán, G.; Martínez, C.; Chávez, R. Cultivable psy-chrotolerant yeasts associated with Antarctic marine sponges. World J. Microbiol. Biotechnol. 2013, 29, 183–189. [Google Scholar] [CrossRef] [PubMed]
  75. Santiago, I.F.; Soares, M.A.; Rosa, C.A.; Rosa, L.H. Lichensphere: A protected natural microhabitat of the non-lichenised fungal communities living in extreme environments of Antarctica. Extremophiles 2015, 19, 1087–1097. [Google Scholar] [CrossRef]
  76. Pazgier, M.; Turkiewicz, M.; Kalinowska, H.; Bielecki, S. The unique cold-adapted extracellular subtilase from psychrophilic yeast Leucosporidium antarcticum. J. Mol. Catal. B Enzym. 2003, 21, 39–42. [Google Scholar] [CrossRef]
  77. Matsui, M.; Kawamata, A.; Kosugi, M.; Imura, S.; Kurosawa, N. Diversity of proteolytic microbes isolated from Antarctic freshwater lakes and characteristics of their cold-active proteases. Polar Sci. 2017, 13, 82–90. [Google Scholar] [CrossRef]
  78. Park, K.S.; Do, H.; Lee, J.H.; Park, S.I.; Kim, E.J.; Kim, S.-J.; Kang, S.-H.; Kim, H.J. Characterization of the ice-binding protein from Arctic yeast Leucosporidium sp. AY30. Cryobiology 2012, 64, 286–296. [Google Scholar] [CrossRef] [PubMed]
  79. Bergauer, P.; Fonteyne, P.-A.; Nolard, N.; Schinner, F.; Margesin, R. Biodegradation of phenol and phenol-related compounds by psychrophilic and cold-tolerant alpine yeasts. Chemosphere 2005, 59, 909–918. [Google Scholar] [CrossRef] [PubMed]
  80. Ma, C.; Hu, B.; Wei, M.-B.; Zhao, J.-H.; Zhang, H.-Z. Influence of matured compost inoculation on sewage sludge composting: Enzyme activity, bacterial and fungal community succession. Bioresour. Technol. 2019, 294, 122165. [Google Scholar] [CrossRef]
  81. Zhang, M.; O’Connor, P.J.; Zhang, J.; Ye, X. Linking soil nutrient cycling and microbial community with vegetation cover in riparian zone. Geoderma 2021, 384, 114801. [Google Scholar] [CrossRef]
  82. Delgado-Baquerizo, M.; Oliverio, A.M.; Brewer, T.E.; Benavent-González, A.; Eldridge, D.J.; Bardgett, R.D.; Maestre, F.T.; Singh, B.K.; Fierer, N. A global atlas of the dominant bacteria found in soil. Science 2018, 359, 320–325. [Google Scholar] [CrossRef] [Green Version]
  83. Bendia, A.G.; Signori, C.N.; Franco, D.C.; Duarte, R.T.D.; Bohannan, B.J.M.; Pellizari, V.H. A Mosaic of Geothermal and Marine Features Shapes Microbial Community Structure on Deception Island Volcano, Antarctica. Front. Microbiol. 2018, 9, 899. [Google Scholar] [CrossRef] [Green Version]
  84. Farias, M.E.; Revale, S.; Mancini, E.; Ordoñez, O.; Turjanski, A.; Cortez, N.; Vazquez, M.P. Genome Sequence of Sphingomonas sp. S17, Isolated from an Alkaline, Hyperarsenic, and Hypersaline Volcano-Associated Lake at High Altitude in the Argentinean Puna. J. Bacteriol. 2011, 193, 3686–3687. [Google Scholar] [CrossRef] [Green Version]
  85. Rincón-Molina, C.I.; Martínez-Romero, E.; Ruiz-Valdiviezo, V.M.; Velázquez, E.; Ruiz-Lau, N.; Rogel-Hernández, M.A.; Villalobos-Maldonado, J.J.; Rincón-Rosales, R. Plant growth-promoting potential of bacteria associated to pioneer plants from an active volcanic site of Chiapas (Mexico). Appl. Soil Ecol. 2020, 146, 103390. [Google Scholar] [CrossRef]
  86. Asaf, S.; Numan, M.; Khan, A.L.; Al-Harrasi, A. Sphingomonas: From diversity and genomics to functional role in environmental remediation and plant growth. Crit. Rev. Biotechnol. 2020, 40, 138–152. [Google Scholar] [CrossRef]
  87. Shen, C.; Gunina, A.; Luo, Y.; Wang, J.; He, J.; Kuzyakov, Y.; Hemp, A.; Classen, A.T.; Ge, Y. Contrasting patterns and drivers of soil bacterial and fungal diversity across a mountain gradient. Environ. Microbiol. 2020, 22, 3287–3301. [Google Scholar] [CrossRef] [PubMed]
  88. Watts, D.J.; Strogatz, S.H. Collective dynamics of ‘small-world’ networks. Nature 1998, 393, 440–442. [Google Scholar] [CrossRef]
  89. Martínez-Rosales, C.; Fullana, N.; Musto, H.; Castro-Sowinski, S. Antarctic DNA moving forward: Genomic plasticity and bio-technological potential. FEMS Microbiol. Lett. 2012, 331, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Guo, X.; Feng, J.; Shi, Z.; Zhou, X.; Yuan, M.; Tao, X.; Hale, L.; Yuan, T.; Wang, J.; Qin, Y.; et al. Climate warming leads to divergent succession of grassland microbial communities. Nat. Clim. Chang. 2018, 8, 813–818. [Google Scholar] [CrossRef] [Green Version]
  91. Mertens, S.; Nijs, I.; Heuer, M.; Kockelbergh, F.; Beyens, L.; Van Kerckvoorde, A.; Impens, I. Influence of High Temperature on End-of-Season Tundra CO2 Exchange. Ecosystems 2001, 4, 226–236. [Google Scholar] [CrossRef]
  92. Mountain Research Initiative EDW Working Group; Pepin, N.; Bradley, R.S.; Diaz, H.F.; Baraer, M.; Caceres, E.B.; Forsythe, N.; Fowler, H.J.; Greenwood, G.; Hashmi, M.Z.; et al. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 2015, 5, 424–430. [Google Scholar] [CrossRef] [Green Version]
  93. Jia, X.; Jin, H.; Zhao, Y.; Yin, H.; Liu, L.J.; Wang, C.; Cui, K.F.; Gu, X.P. Characteristics of Climate Change in Changbai Mountain from 1958 to 2015. J. Beihua Univ. 2017, 18, 727–731. [Google Scholar]
  94. Davidson, E.A.; Janssens, I.A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nat. Cell Biol. 2006, 440, 165–173. [Google Scholar] [CrossRef] [PubMed]
  95. Pritchard, S.G. Soil organisms and global climate change. Plant Pathol. 2011, 60, 82–99. [Google Scholar] [CrossRef]
  96. Salazar, A.; Rousk, K.; Jónsdóttir, I.S.; Bellenger, J.; Andrésson, Ó.S. Faster nitrogen cycling and more fungal and root biomass in cold ecosystems under experimental warming: A meta-analysis. Ecology 2020, 101, e02938. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of Tianchi Volcano in China (A) and sampling area in northern and western margin of Tianchi crater (B). The map was generated from ArcGIS Online (https://maps.arcgis.com/index.html, accessed on 1 January 2021) [34].
Figure 1. Location of Tianchi Volcano in China (A) and sampling area in northern and western margin of Tianchi crater (B). The map was generated from ArcGIS Online (https://maps.arcgis.com/index.html, accessed on 1 January 2021) [34].
Life 11 00280 g001
Figure 2. Rarefaction curve of fungal OTUs (A) and bacterial OTUs (B) for α-diversity analyses.
Figure 2. Rarefaction curve of fungal OTUs (A) and bacterial OTUs (B) for α-diversity analyses.
Life 11 00280 g002
Figure 3. Rank abundance curve of fungal OTUs (A) and bacterial OTUs (B) for α-diversity analyses.
Figure 3. Rank abundance curve of fungal OTUs (A) and bacterial OTUs (B) for α-diversity analyses.
Life 11 00280 g003
Figure 4. The relative abundance of fungi (A) and bacteria (B) at class level.
Figure 4. The relative abundance of fungi (A) and bacteria (B) at class level.
Life 11 00280 g004
Figure 5. The relative abundance of fungi (A) and bacteria (B) at genus level.
Figure 5. The relative abundance of fungi (A) and bacteria (B) at genus level.
Life 11 00280 g005
Figure 6. The Venn diagram displaying the distribution of shared and unique OTUs of the fungal community (A) and bacterial community (B) in two sides of the crater. The number of OTUs in each side is indicated in the respective circle.
Figure 6. The Venn diagram displaying the distribution of shared and unique OTUs of the fungal community (A) and bacterial community (B) in two sides of the crater. The number of OTUs in each side is indicated in the respective circle.
Life 11 00280 g006
Figure 7. Principal co-ordinate analysis (PCoA) plots of fungal communities in northern and western crater soil based on Bray–Curtis distance (A), weighted Unifrac distance (B) and unweighted Unifrac distance (C). The r and p-values of analysis of similarity (ANOSIM) were shown respectively in each figure (p < 0.05).
Figure 7. Principal co-ordinate analysis (PCoA) plots of fungal communities in northern and western crater soil based on Bray–Curtis distance (A), weighted Unifrac distance (B) and unweighted Unifrac distance (C). The r and p-values of analysis of similarity (ANOSIM) were shown respectively in each figure (p < 0.05).
Life 11 00280 g007
Figure 8. Principal co-ordinate analysis (PCoA) plots of bacterial communities in northern and western crater soil based on Bray–Curtis distance (A), weighted Unifrac distance (B) and unweighted Unifrac distance (C). The r and p-values of analysis of similarity (ANOSIM) were shown respectively in each figure (p > 0.05).
Figure 8. Principal co-ordinate analysis (PCoA) plots of bacterial communities in northern and western crater soil based on Bray–Curtis distance (A), weighted Unifrac distance (B) and unweighted Unifrac distance (C). The r and p-values of analysis of similarity (ANOSIM) were shown respectively in each figure (p > 0.05).
Life 11 00280 g008
Figure 9. Soil OTUs network analysis in the Tianchi crater soil (Fruchterman Reingold layout). (A) Network of fungal community; (B) network of bacterial community. Each node represents an OTU indicating an individual species. Color codes for nodes belonging to different dominant phyla. The node size is proportional to the degree (degree = number of direct correlations to a node). Positive interactions are displayed as red edges and negative interactions are displayed as blue edges.
Figure 9. Soil OTUs network analysis in the Tianchi crater soil (Fruchterman Reingold layout). (A) Network of fungal community; (B) network of bacterial community. Each node represents an OTU indicating an individual species. Color codes for nodes belonging to different dominant phyla. The node size is proportional to the degree (degree = number of direct correlations to a node). Positive interactions are displayed as red edges and negative interactions are displayed as blue edges.
Life 11 00280 g009
Table 1. Details of sampling sites.
Table 1. Details of sampling sites.
Sampling SitesLatitudeLongitudeAltitude (m a.s.l)Sample Type
Northern margin of Tianchi crater42°1′33″ N128°4′8″ E26206 Soil stone mixture samples
Western margin of Tianchi crater41°59′48″ N128°1′44″ E24606 Soil stone mixture samples
Table 2. Primers used for culturable microbe identification.
Table 2. Primers used for culturable microbe identification.
MicrobesLocusPrimer NameDirectionSequenceTarget Region
FungiInternal transcribed spacer (ITS)ITS1Forward5′-TCCGTAGGTGAACCTGCGG-3′ITS
ITS4Reverse5′-TCCTCCGCTTATTGATATGC-3′
BacteriaV3 to V4 hypervariable regionU341FForward5′-ACTCCTACGGGAGGCAGCAG-3′16S rRNA
U806RReverse5′-GGACTACHVGGGTWTCTAAT-3′
Table 3. Taxonomic groups of fungi isolated from Tianchi Crater of Changbai Mountain.
Table 3. Taxonomic groups of fungi isolated from Tianchi Crater of Changbai Mountain.
PhylumClassOrderGenus (Number of Species)
AscomycotaDothideomycetesCapnodialesCladosporium (4)
ConiosporialesConiosporium (1)
DothidealesAureobasidium (1)
PleosporalesAlternaria (3), Curvularia (1), Didymella (4), Epicoccum (1), Juxtiphoma (1), Leptosphaeria (1), Paraphaeosphaeria (1), Paraphoma (1), Phaeosphaeria (2), Phoma (2)
EurotiomycetesChaetothyrialesExophiala (1)
EurotialesByssochlamys (1), Penicillium (5)
Incertae sedisIncertae sedisTricellula (1)
LeotiomycetesHelotialesXenochalara (1)
HelotialesHelotiales sp. (1)
SaccharomycetesSaccharomycetalesCandida (1)
SordariomycetesAmphisphaerialesMicrodochium (1)
HypocrealesFusarium (3), Nectria (1), Tangerinosporium (1), Trichoderma (5)
MicroascalesCephalotrichum (1)
XylarialesXylaria (1)
BasidiomycotaMicrobotryomycetesLeucosporidialesLeucosporidium (1)
TremellomycetesCystofilobasidialesTausonia (1)
TremellalesNaganishia (1)
MucoromycotaMucoromycetesMucoralesMucor (1)
Table 4. α-diversity indexes (OTU richness, Shannon, Chao1 and Coverage) data of fungal and bacterial communities from crater soil samples from North (N) and West (W) slopes. OTU = Operational Taxonomic Unit.
Table 4. α-diversity indexes (OTU richness, Shannon, Chao1 and Coverage) data of fungal and bacterial communities from crater soil samples from North (N) and West (W) slopes. OTU = Operational Taxonomic Unit.
SampleFungal CommunityBacterial Community
OTU NumbersShannonChaoCoverage (%)OTU NumbersShannonChaoCoverage (%)
Mean ± SD368 ± 1082.95 ± 0.74446.14 ± 10899.8 ± 0.11009 ± 3654.72 ± 0.931240.98 ± 393.6099.3 ± 0.2
N13242.36395.6799.84622.59607.0199.6
N24803.45526.3699.85083.38707.9499.5
N33082.44426.3299.88254.781061.3599.4
N43553.35412.3799.87824.621049.1499.4
N53183.42397.4499.810964.841297.2799.3
N65244.04555.5299.913015.581481.0699.3
W14993.20588.2599.715585.751828.5699.1
W24043.01496.6499.710685.221326.7499.3
W34513.86554.2699.713065.441569.1899.2
W43582.34477.2899.715325.471840.9599.0
W51961.54255.1199.98904.491101.6899.3
W62022.39268.599.97804.521020.999.4
Table 5. Network topological properties of fungal and bacterial communities in the crater soil of the Tianchi Volcano.
Table 5. Network topological properties of fungal and bacterial communities in the crater soil of the Tianchi Volcano.
Network MetricsCommunity
FungiBacteria
Number of nodes7784
Total number of edges136426
Number and percentage of positive correlations108 (79.41%)302 (70.89%)
Number and percentage of negative correlations28 (20.59%)124 (29.11%)
Average degree3.48710.143
Network diameter127
Modularity0.7090.341
Number of communities1410
Average clustering coefficient0.5210.565
Average path length4.9822.661
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, X.; Pecoraro, L. Analysis of Soil Fungal and Bacterial Communities in Tianchi Volcano Crater, Northeast China. Life 2021, 11, 280. https://doi.org/10.3390/life11040280

AMA Style

Wang X, Pecoraro L. Analysis of Soil Fungal and Bacterial Communities in Tianchi Volcano Crater, Northeast China. Life. 2021; 11(4):280. https://doi.org/10.3390/life11040280

Chicago/Turabian Style

Wang, Xiao, and Lorenzo Pecoraro. 2021. "Analysis of Soil Fungal and Bacterial Communities in Tianchi Volcano Crater, Northeast China" Life 11, no. 4: 280. https://doi.org/10.3390/life11040280

APA Style

Wang, X., & Pecoraro, L. (2021). Analysis of Soil Fungal and Bacterial Communities in Tianchi Volcano Crater, Northeast China. Life, 11(4), 280. https://doi.org/10.3390/life11040280

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop