Next Article in Journal
The Amylolytic Regulator AmyR of Aspergillus niger Is Involved in Sucrose and Inulin Utilization in a Culture-Condition-Dependent Manner
Previous Article in Journal
Races CYR34 and Suwon11-1 of Puccinia striiformis f. sp. tritici Played an Important Role in Causing the Stripe Rust Epidemic in Winter Wheat in Yili, Xinjiang, China
Previous Article in Special Issue
Development and Determinants of Topsoil Bacterial and Fungal Communities of Afforestation by Aerial Sowing in Tengger Desert, China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 9
Issue 4
10.3390/jof9040437
jof-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrative Study of Mycobiome in Different Habitats from a High Arctic Region: Diversity, Distribution, and Functional Role

by
Xiufei Chen
1,
Dong Yan
2,
Liyan Yu
1,* and
Tao Zhang
1,*
1
China Pharmaceutical Culture Collection, Institute of Medicinal Biotechnology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100050, China
2
Xinxiang Key Laboratory of Pathogenic Biology, Department of Pathogenic Biology, School of Basic Medical Sciences, Xinxiang Medical University, Xinxiang 453003, China
*
Authors to whom correspondence should be addressed.
J. Fungi 2023, 9(4), 437; https://doi.org/10.3390/jof9040437
Submission received: 27 January 2023 / Revised: 24 March 2023 / Accepted: 28 March 2023 / Published: 3 April 2023
(This article belongs to the Special Issue Fungal Diversity in Various Environments)
Browse Figures
Versions Notes

Abstract

:
In the Arctic ecosystems, fungi are crucial for interactions between soil and plants, the cycling of nutrients, and the transport of carbon. To date, no studies have been conducted to thoroughly examine the mycobiome and its functional role in various habitats of the High Arctic region. The aim was to unravel the mycobiome in the nine habitats (i.e., soil, lichen, vascular plant, moss, freshwater, seawater, marine sediment, dung, and marine alga) in the Ny-Ålesund Region (Svalbard, High Arctic) using a high-throughput sequencing approach. A total of 10,419 ASVs were detected. Among them, 7535 ASVs were assigned to unidentified phyla, while the remaining 2884 ASVs were assigned to 11 phyla, 33 classes, 81 orders, 151 families, 278 genera, and 261 species that were known. The distribution of the mycobiome was driven by habitat specificity, indicating that habitat filtering is a crucial factor influencing the fungal assemblages at a local scale in this High Arctic region. Six growth forms and 19 fungal guilds were found. The ecological guild (e.g., lichenized, ectomycorrhizal) and growth form (e.g., yeast, thallus photosynthetic) varied significantly among various habitats. In addition, the occurrence of 31 fungal species that are considered to be potential pathogens was determined. These results will increase our understanding of fungal diversity and its functional significance in this distinctive High Arctic area and thereby establish the groundwork for prediction about how the mycobiome will alter in various environments as a result of anticipated climate change.

1. Introduction

Global climate change is causing the Arctic to warm at unprecedented rates, which has a significant impact on the Arctic ecosystems by altering weather patterns, lowering sea ice cover, melting massive areas of permafrost soil, and changes in Arctic vegetation [1,2,3]. The Arctic tundra, which makes up around 5% of the planet’s land area, is one of the harshest habitats on earth due to its low temperatures, repeated freeze-thaw cycles, wet-dry cycles, and lack of organic matter. Up to now, understanding the effects of warming on the Arctic ecosystem has been a difficult task. The Arctic has a variety of habitats that are home to numerous potentially viable, mainly uncharacterized microorganisms. Warming has been demonstrated to change the composition of the microbial communities in Arctic soil, including bacteria [4] and fungi [4,5,6]. Therefore, to predict how Arctic ecosystems will respond to rising temperatures, it is crucial to have a thorough understanding of the region’s microbiome [7].
The majority of earlier studies used conventional isolation techniques to investigate the diversity of cultivated fungi in Arctic habitats (e.g., soils [8,9], lichens [10], freshwater [11], and cryoconite holes [12]). Due to their selectivity, the cultivation-based approaches, however, probably do not accurately reveal the diversity of the mycobiome. Numerous studies have recently used targeted amplicon sequencing (i.e., Roche 454 pyrosequencing, Illumina sequencing) to characterize the diversity and composition of the mycobiome in the Arctic, including those found in the soils [13,14,15,16], marine sediments [17], lichens [18], mosses [16,19], vascular plants [16,20,21,22], seawater [23,24], and freshwater [16,25,26]. However, there are no comprehensive studies that unravel and compare the diversity and distribution of the mycobiome in various habitats from the High Arctic region. The fungal assemblages at the local scale of the High Arctic may be impacted by habitat overlap or specialization, but this is still unclear.
In the Arctic, shifts in the composition of fungal functional groups are considered to have an impact on permafrost carbon storage, ecosystem respiration, and nutrient cycling [27]. The ecological role of fungi in the Arctic ecosystem, however, is only partially understood (e.g., soil [13,27,28,29], vegetation [19], and seawater [30]). In addition, some environmental fungi may evolve pathogenic potential as a result of climate warming, and the presence of potential fungal pathogens in Arctic environments (such as snow, subglacial ice, glacial meltwater, and air) has been observed [31]. However, there is still a considerable information gap between the diversity of fungi and their ecological functions in the Arctic ecosystem.
To address this knowledge gap, we used high-throughput sequencing to analyze the mycobiome associated with the nine habitat types (i.e., soil, lichen, vascular plant, moss, freshwater, seawater, sediment, dung, and marine alga) in the Arctic region. The aims of the present study were to (1) comprehensively reveal the diversity and functional role of the mycobiome in various habitats from a High Arctic region; (2) assess how and to what extent the habitat types influence the mycobiome at the local scale. We hypothesized that habitat selection instead of habitat overlap would determine the mycobiome in different habitats of the High Arctic. Knowledge of the diversity and functional role of the mycobiome in the High Arctic will be enhanced by the findings of the present study.

2. Materials and Methods

2.1. Sampling Area

The sampling area was situated in the Ny-Ålesund Region (78°55′ N, 11°56′ E), which lies in the north-western region of Spitsbergen, Svalbard. Svalbard is an Arctic archipelago that lies completely within the High Arctic Zone and is geographically separated from mainland Eurasia. Polar cyclones and the circulation of the North Atlantic Ocean both have an impact on the Svalbard archipelago. According to Maturilli et al. [32], the mean temperatures of the Ny-Ålesund Region range from −17.0 to −3.8 °C in January to 4.6 to 6.9 °C in July. There is no sunset from 18 April to 24 August, and there is no sunrise from 25 October to 17 February (polar night). A total of 114 samples from nine habitats were collected during China’s Arctic expedition in July 2013 (Figure 1, Table S1), comprising 69 samples in the present study (including lichen, marine alga, seawater, marine sediment, dung, and moss) and 45 samples in a previous study (including soil, freshwater, moss, and vascular plant) [16]. The Svalbard Science Forum, the Chinese Arctic and Antarctic Administration (CAA), and the State Oceanic Administration (SOA) of China all provided their permission for the sample collection. Before being flown to our laboratory in China, the samples were frozen at –20 °C in the Yellow River Station (China) for approximately 20 days.

2.2. DNA Extraction, PCR Amplification, and Illumina Sequencing

(1) For the samples of marine sediment and dung, a PowerSoil DNA Isolation Kit (MO BIO Laboratories Inc., Carlsbad, CA, USA) was used to extract the DNA from 0.25 g of each sample according to the manufacturer’s instructions. (2) For the samples of freshwater and seawater (1000 mL of each sample), a PowerWater DNA Isolation Kit (MO BIO Laboratories Inc., Carlsbad, CA, USA) was used to extract the DNA from the membrane-bound biomass according to the manufacturer’s instructions. (3) For the samples of lichen, moss, and marine alga, before DNA extraction, the tissues were surface sterilized, as described by Zhang et al. [12]. A PowerSoil DNA Isolation Kit (MO BIO Laboratories Inc., Carlsbad, CA, USA) was then used to extract the DNA from the sterilized tissues, as described by Zhang et al. [18].
The nuclear ribosomal internal transcribed spacer 1 (ITS1) of the fungi was amplified from the acquired DNA using ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) [33] and ITS2 (5′-GCTGCGTTCTTCATCGATGC-3′) [34] with a 10-nucleotide barcode. PCR amplification was then performed, and the reaction mixture (50 μL) contained the template DNA (about 10 ng), 10× EasyTaq buffer (5 μL), dNTPs (2.5 mM, 2 μL), EasyTaq DNA polymerase (5 units, 1 μL), primers (10 μM, 1 μL), and nuclease-free water. The PCR program began with an initial denaturation step at 95 °C for 2 min, followed by 30 cycles of 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s, and a final extension step at 72 °C for 5 min [18]. Amplicons were extracted, purified, and paired-end sequenced (2 × 300 bp) on an Illumina MiSeq platform by Shanghai Majorbio Bio-pharm Technology Co., Ltd. (Shanghai, China). The negative control samples were subjected to both DNA extraction and PCR. Since no quantifiable DNA was discovered in these negative controls, further sequencing was not performed on them.

2.3. Sequencing Data Processing

FLASH software was used to combine the paired-end reads [35]. Subsequently, based on the unique barcodes of each sample, reads were assigned to each one. QIIME 2 version 2022.01 was used to process the unprocessed raw demultiplexed sequences [36]. The DADA2 plugin, a component of QIIME 2, was used to denoise, dereplicate, and filter (for chimeras) the paired-end reads [37]. The bases with quality scores below 35 were removed from the raw reads. The first 26 nucleotides of the 5′ end were removed from both the forward and reverse reads. The forward and reverse sequences both had their 3′ ends cut off at locations 240 and 230, respectively. The read number was set to 100,000 for training the error model. Denovo clustering (100% similarity) was conducted using VSEARCH [38], which is a part of QIIME 2. Amplicon sequence variations (ASVs) were assigned to specific taxa using the classify-sklearn classifier [39] against the UNITE Fungi version 8.3 reference database [40], which has been trained to recognize the fungal ITS1 gene. ASVs with read counts under five were eliminated. The data were then subsampled, with each sample containing 14,866 sequences. Using the FFT-NS-1 method, MAFFT version 7 was used to align with typical ASV sequences, and the Average linkage (UPGMA) method was utilized to build a phylogenetic tree [41].

2.4. Statistical Analyses

MicrobiomeAanalyst was used to analyze the data obtained [42]. Non-metric multidimensional scaling (NMDS) ordination showing the mycobiome composition (at the ASV level) in different samples was performed using unweighted UniFrac distance. An analysis of similarity (ANOSIM) was also used to reveal the dissimilarity of the mycobiome among nine habitat types at the ASV level (unweighted UniFrac distances). Dendrogram analysis of the mycobiome in the 114 samples from the nine habitats was performed (average clustering algorithm and unweighted UniFrac distance). The number of fungal ASVs in each of the nine habitats was represented by a Venn diagram. To investigate how different fungal genera interact with one another in different habitats, correlation network analysis was performed using Spearman’s rank correlation test (permutation: 100, p-value threshold: 0.01, correlation threshold: 0.70). The functional compositions of the mycobiome were identified using FungalTraits version 1.2 [43]. Using a factorial Kruskal-Wallis test, the distinct taxonomic units (phylum, class, order, family, and genus) and functional roles (guild and growth form) among the nine habitat types were identified using the linear discriminant analysis effect size (LEfSe) analysis. The Atlas of Clinical Fungi was used to determine potential human pathogens [44].

3. Results

3.1. Sequence Data

From 114 samples, a total of 4,265,338 raw reads were obtained. The raw reads per sample ranged from 19,836 to 105,159 (Table S1). The trimmed reads per sample ranged from 14,896 to 97,477 (Table S2). Following the removal of the uncommon ASVs (<five reads), the total number of reads for 10,453 ASVs in 114 samples was decreased to 3,491,225. After each sequence library was rarefied to 14,866 reads, 1,694,724 reads from 114 samples and 10,419 ASVs were preserved. The ASV number per sample ranged from 2 to 681. Overall, 7535 ASVs were assigned to unidentified phyla, while the remaining 2884 ASVs were assigned to 11 phyla, 33 classes, 81 orders, 151 families, 278 genera, and 261 species that were known (Table S3).

3.2. Diversity and Composition of the Mycobiome in the Nine Habitats

The nine habitats have different mycobiome compositions, as shown by the stacked bar plots (Figure 2). Reads assigned to unknown phyla predominated in most samples of the nine habitats, especially in seawater and marine sediment. The eleven known phyla included Ascomycota, Basidiomycota, Chytridiomycota, Rozellomycota, Mortierellomycota, Olpidiomycota, Monoblepharomycota, Zoopagomycota, Aphelidiomycota, Mucoromycota, and Kickxellomycota (Figure 2a). At the order level, Thelobolales was more abundant than other orders in the habitat of animal dung, whereas Hypocreales was dominant in the seawater habitat. In addition, the order Agaricales was abundant in the soil habitat and the order Helotiales was abundant in the habitats of moss and vascular plant (Figure 2b).

3.3. Dissimilarity of the Mycobiome among the Nine Habitats

The NMDS plot revealed that the mycobiome compositions at the ASV level were different among the habitat types (Figure 3a). The results of the ANOSIM test also revealed that there were significant differences in the mycobiome compositions between the nine distinct habitats (R = 0.79507, p < 0.001). The results from the pairwise ANOSIM comparisons revealed differences in the mycobiome compositions among various habitats (Table 1). For example, seawater and marine sediment harbored significantly different mycobiome compositions (R = 0.90845, p < 0.001), and marine sediment and vascular plant harbored significantly different mycobiome compositions (R = 1, p < 0.001). Only one ASV was shared by all nine habitats, according to the Venn diagram (Figure 3b). The linkages between the mycobiome and the nine habitat types were also revealed using Dendrogram analysis, showing that the mycobiome in the samples clustered by habitat type (Figure S1). For example, 24 seawater samples grouped and were isolated from those in the other eight habitats; 8 marine sediment samples clustered together and were also isolated from those in the other samples; and 12 dung samples grouped and were also isolated from those in the other habitats.
The co-occurrence networks demonstrated that there was a strong relationship between the fungal genera and habitat types (Figure 4). For instance, six fungal genera were found in the same module and were present in the sediment, including Pyxidiophora, Trichosporon, Libertasomyces, Infundichalara, Lasiodiplodia, and Spissiomyces (Figure 4a). In addition, nine fungal genera (i.e., Tausonia, Curvularia, Clarireedia, Stagonosporopsis, Beauveria, Candida, Nigrospora, Kurtzmaniella, and Paraboeremia) were in the same module and were harbored in the marine alga (Figure 4b).
The LEfSe analyses showed that some fungal phyla, classes, orders, families, and genera were significantly different among the different habitat types (Figure 5, Figures S2–S4). For example, at the phylum level, Ascomycota was significantly abundant in the habitat of vascular plant, Basidiomycota was significantly abundant in the soil habitat, and Chytridiomycota was significantly abundant in the habitats of marine alga and freshwater (Figure S2). At the order level, Theleobolales was significantly abundant in the dung habitat, and Helotiales was significantly abundant in the habitat of vascular plant. Moreover, Agaricales, Sebacinales, and Chaetothyriales were significantly abundant in the soil habitat, while Capnodiales and Dothideales were significantly abundant in the moss habitat (Figure 5a). At the genus level, Xenoacremonium and Fontanospora were significantly abundant in the habitat of vascular plant; Mrakia and Glaciozyma were significantly abundant in the habitat of marine alga; Fusarium was significantly abundant in the habitat of seawater; Cortinarius and Inocybe were significantly abundant in the habitat of soil (Figure 5b).

3.4. Ecological Role of the Mycobiome in the Nine Habitats

A total of 19 functional guilds and 6 growth forms were discovered in the 114 samples. Table S3 showed the functional assignments of the fungal ASVs identified at the general level. Saprophytic fungi (989 ASVs) dominated the mycobiome in these samples, followed by plant pathogens (185 ASVs), lichenized fungi (116 ASVs), endophytic fungi (45 ASVs), and ectomycorrhizal fungi (41 ASVs).
The results of the ANOSIM tests revealed that the fungal guilds (R = 0.50281, p < 0.001) and growth forms (R = 0.50285, p < 0.001) among the nine different habitats differed significantly. Many significant differences in the fungal functional role were also observed among different habitats using LEfSe analysis (Figure 6). In terms of growth forms, dimorphic yeast and yeast were significantly abundant in the marine alga habitat; thallus photosynthetic was significantly abundant in the lichen habitat; and filamentous mycelium was significantly abundant in the vascular plant habitat (Figure 6a). In terms of the ecological guild, lichenized fungi were significantly abundant in the lichen habitat; ectomycorrhizal fungi and soil saprotroph were significantly abundant in the soil habitat; and root endophyte was significantly abundant in the moss habitat (Figure 6b).
The occurrence of 31 fungal species that were considered to be potential human pathogens was determined. The common species included Alternaria tenuissima (41 samples in 7 habitats), Malassezia restricta (41 samples in 7 habitats), Aspergillus sydowii (17 samples in 8 habitats), Didymella glomerata (12 samples in 5 habitats), Filobasidium magnum (12 samples in 4 habitats), Cladosporium cladosporioides (11 samples in 6 habitats), and Naganishia albida (10 samples in 3 habitats) (Table S4).

4. Discussion

Many studies have indicated that changes in environmental conditions have a significant impact on the fungal communities in the Arctic [4,5,6], and thereby it is urgent for an integrative study to unravel the diversity and functional role of the mycobiome in the Arctic that exists now. Studies on the mycobiome have, to date, largely focused on limited Arctic habitats, and there hasn’t been much research conducted to investigate different habitats at a local scale.
The most prevalent fungal ASVs in the present study were categorized into unidentified phyla, which suggests the dominance of undescribed fungal species. The common phyla identified were Ascomycota, Basidiomycota, Chytridiomycota, Mortierellomycota, Mucoromycota, and Rozellomycota, which were reported in previous 454 pyrosequencing studies [15,17,18,19,22,26]. ASVs belonging to the infrequently reported phyla Olpidiomycota, Monoblepharomycota, Zoopagomycota, Aphelidiomycota, and Kickxellomycota, were also discovered using Illumina sequencing (Figure 2a). Furthermore, a total of 81 orders, 278 genera, and 261 species were found in this area’s nine habitats. In earlier researches, 28 orders and 95 genera were detectable in Arctic soils using 454 pyrosequencing [15], while 24 orders, 90 genera, and 76 species were detectable in freshwater [26], 28 orders, 51 genera, and 31 species were detectable in vascular plants [22], 40 orders, 78 genera, and 55 species were detectable in lichens [18], and 25 orders, 30 genera, and 30 species were detectable in marine sediments [17]. The results provided a more integrative understanding of the fungal taxonomic diversity in a High Arctic region.
The abundance of some fungal phyla, classes, orders, families, and genera differed significantly among the nine habitats. For example, the order Agariales was abundant in the soil habitat, and the order Helotiales was abundant in the moss and vascular plant habitats (Figure 2b). The roots of the Arctic vascular plant (Ericaceae) host a diverse mycobiome dominated by the order Helotiales [45]. In previous 454 pyrosequencing analyses of fungal communities in Arctic mosses, the order Helotiales was the dominant component of the communities inhabiting Arctic mosses and vascular plants [19]. In our previous study, a total of 21 macrofungal species belonging to the order Agaricales were identified in the Ny-Ålesund Region [46].
In the present study, 19 functional guilds and 6 growth forms were identified, suggesting that the mycobiome in the High Arctic may serve a range of functional roles in various habitats. We found saprophytic fungi (989 ASVs) dominated the fungal community in this High Arctic region. Saprotrophic fungi increase the availability of soil nutrients by actively decomposing organic matter through the secretion of a range of hydrolases, such as proteinases, cellulases, and laccases [47]. Interestingly, many fungal functional groups (i.e., guilds and growth forms) were also significantly different among the nine habitats (Figure 6), which was inconsistent with the results of the taxonomic groups (i.e., order and genera) (Figure 5). For example, ectomycorrhizal fungi (order Agaricales, genera Cortinarius and Inocybe) were more abundant in the soil habitat, lichenized fungi (order Lecanorales, Umbiliariales, and Pertusariales, genera Flavocetraria, Cetrariella, Umbilicaria, and Ochrolechia) were more abundant in the lichen habitat, root endophytes (genera Serendipita) were more abundant in the moss habitat, and yeasts (order Cystofilobasidiales, genera Glaciozyma and Mrakia) were more abundant in the marine alga habitat. Ectomycorrhizal fungi, which provide plants with mineral nutrients and water from the soil, are necessary for Arctic vascular plants to survive (e.g., Betula, Dryas, and Salix) [48]. The root endophytic fungus Serendipita belongs to the order Sebacinales and can alter the root architecture and improve phosphorus (P) and zinc (Zn) uptake [49]. Therefore, habitat filtering may affect the colonization of fungi, with some fungal species being chosen based on their functional characteristics [13].
The habitat specificity of the mycobiome at the local scale may be attributed to variations in environmental factors (such as the physicochemical proprieties, nutritional levels, antagonistic factors, etc.). There are many unmeasured environmental factors in the nine habitats of the Ny-Ålesund Region, which may affect the diversity and composition of the mycobiome in this region. In Western Greenland (Low Arctic), varied vegetation types and environmental conditions, such as pH and water content, were shown to be associated with changes in the fungal community composition and functional role [28].
The composition of the mycobiome may be affected by climate change directly or indirectly [4]. As a result of climate warming, the melting of permafrost can alter aquatic fungal diversity and fungal community composition through interactions with carbon released from sediments [25]. Some fungi can be utilized as bioindicators of changes taking place in the Arctic environment, such as the genus Cladosporium [50]. In addition, permafrost may be a reservoir of potential fungal pathogens [51], which may be released as a result of climate warming. Previous research has shown the presence of potential fungal pathogens in Arctic environments, such as snow, ice, meltwater, and air. These pathogens include Aureobasidium melanogenum, Naganishia albida, and Rhodotorula mucilaginosa, which can grow at 37 °C, exhibit β-hemolytic activity, and are resistant to azoles and echinocandins [31]. Potential human fungal pathogens (31 species) were discovered in the present study, according to the Atlas of Clinical Fungi [44]. Among the potential fungal pathogens (31 species) observed in the present study, seven fungal species are listed as Biosafety Level 2 (BSL-2), including Aspergillus fumigatus, Fusarium solani, Malassezia globosa, Rhizopus arrhizus, Sarocladium kiliense, Trichosporon asahii, and Wickerhamomyces anomalus. To determine the fungal virulence detected in the Arctic habitats, which is essential for the pathogenicity of fungi, additional analyses of fungal isolates are required in further studies.

5. Conclusions

This integrative study reveals the diversity and distribution of the mycobiome in the nine distinct habitats and highlights their functional role in the High Arctic region. We also found that, while the mycobiome in distinct habitats was related locally through dispersal, habitat specialization rather than overlap drove the mycobiome distribution in the High Arctic region. This suggests that habitat filtering is a significant determinant of fungal assemblages in this Arctic region. Our data may lay a foundation for speculation on the changes in the mycobiome in different habitats due to climate warming in the High Arctic regions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jof9040437/s1. Figure S1: Dendrogram showing fungal communities in the 114 samples collected from the nine habitats; Figure S2: LEfSe analysis showing the fungal phyla that are significantly different among the nine habitats of the Ny-Ålesund Region; Figure S3: LEfSe analysis showing the fungal classes that are significantly different among the nine habitats of the Ny-Ålesund Region; Figure S4: LEfSe analysis showing the fungal families that are significantly different among the nine habitats of the Ny-Ålesund Region; Table S1: Information on the 114 samples collected from the Ny-Ålesund Region; Table S2: Information on the 11,850 fungal ASVs in the 114 samples collected from the Ny-Ålesund Region; Table S3: Information on the 10,419 fungal ASVs in the 114 samples collected from the Ny-Ålesund Region; Table S4: An overview of potentially pathogenic fungi found in the nine habitats from the Ny-Ålesund Region.

Author Contributions

The study was designed by T.Z. and L.Y. T.Z. collected the samples, conducted the lab work and data analysis, and wrote the manuscript; X.C. conducted parts of the lab work and data analysis; D.Y. conducted parts of the data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (NSFC) (No. 32270013); the projects of the Chinese Arctic and Antarctic Administration and State Oceanic Administration (Nos. 2013YR06006 and 2013YR05005); the CAMS Innovation Fund for Medical Sciences (Grant no. 2021-I2M-1-055); the National Microbial Resource Center (Grant no. NMRC-2022-3); and the Non-profit Central Research Institute Fund of CAMS (Grant no. 2021-PT350-001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the raw sequence files are available in the National Center for Biotechnology Information under the BioProject ID PRJNA746944.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Biskaborn, B.K.; Smith, S.L.; Noetzli, J.; Matthes, H.; Vieira, G.; Streletskiy, D.A.; Schoeneich, P.; Romanovsky, V.E.; Lewkowicz, A.G.; Abramov, A.; et al. Permafrost is warming at a global scale. Nat. Commun. 2019, 10, 264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hassett, B.T.; Gradinger, R. Chytrids dominate arctic marine fungal communities. Environ. Microbiol. 2016, 18, 2001–2009. [Google Scholar] [CrossRef] [PubMed]
  3. Tape, K.E.N.; Sturm, M.; Racine, C. The evidence for shrub expansion in Northern Alaska and the Pan-Arctic. Glob. Chang. Biol. 2006, 12, 686–702. [Google Scholar] [CrossRef]
  4. Deslippe, J.R.; Hartmann, M.; Simard, S.W.; Mohn, W.W. Long-term warming alters the composition of Arctic soil microbial communities. FEMS Microbiol. Ecol. 2012, 82, 303–315. [Google Scholar] [CrossRef] [Green Version]
  5. Morgado, L.N.; Semenova, T.A.; Welker, J.M.; Walker, M.D.; Smets, E.; Geml, J. Summer temperature increase has distinct effects on the ectomycorrhizal fungal communities of moist tussock and dry tundra in Arctic Alaska. Glob. Chang. Biol. 2015, 21, 959–972. [Google Scholar] [CrossRef] [Green Version]
  6. Timling, I.; Taylor, D.L. Peeking through a frosty window: Molecular insights into the ecology of Arctic soil fungi. Fungal Ecol. 2012, 5, 419–429. [Google Scholar] [CrossRef]
  7. Blaud, A.; Lerch, T.Z.; Phoenix, G.K.; Osborn, A.M. Arctic soil microbial diversity in a changing world. Res. Microbiol. 2015, 166, 796–813. [Google Scholar] [CrossRef]
  8. Hafizah, S.; Alias, S.; Siang, H.; Smykla, J.; Pang, K.-L.; Guo, S.-Y.; Convey, P. Studies on diversity of soil microfungi in the Hornsund area, Spitsbergen. Pol. Polar Res. 2013, 34, 39–54. [Google Scholar] [CrossRef]
  9. Singh, S.; Singh, S.; Yadav, L.S.; Singh, P.; Ravindra, R. Filamentous Soil Fungi from Ny-Ålesund, Spitsbergen, and Screening for Extracellular Enzymes. ARCTIC 2012, 65, 45–55. [Google Scholar] [CrossRef]
  10. Zhang, T.; Wei, X.L.; Wei, Y.Z.; Liu, H.Y.; Yu, L.Y. Diversity and distribution of cultured endolichenic fungi in the Ny-Alesund Region, Svalbard (High Arctic). Extremophiles 2016, 20, 461–470. [Google Scholar] [CrossRef]
  11. Zhang, T.; Zhao, L.; Yu, C.; Wei, T.; Yu, L. Diversity and bioactivity of cultured aquatic fungi from the High Arctic region. Adv. Polar Sci. 2017, 28, 29–42. [Google Scholar]
  12. Borzęcka, J.; Suchodolski, J.; Dudek, B.; Matyaszczyk, L.; Spychała, K.; Ogórek, R. The first comprehensive biodiversity study of culturable fungal communities inhabiting cryoconite holes in Werenskiold Glacier on Spitsbergen (Svalbard Archipelago, Arctic). Biology 2022, 11, 1224. [Google Scholar] [CrossRef]
  13. Grau, O.; Geml, J.; Pérez-Haase, A.; Ninot, J.M.; Semenova-Nelsen, T.A.; Peñuelas, J. Abrupt changes in the composition and function of fungal communities along an environmental gradient in the high Arctic. Mol. Ecol. 2017, 26, 4798–4810. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, T.; Wang, N.; Yu, L. Soil fungal community composition differs significantly among the Antarctic, Arctic, and Tibetan Plateau. Extremophiles 2020, 24, 821–829. [Google Scholar] [CrossRef]
  15. Zhang, T.; Wang, N.F.; Liu, H.Y.; Zhang, Y.Q.; Yu, L.Y. Soil pH is a Key Determinant of Soil Fungal Community Composition in the Ny-Alesund Region, Svalbard (High Arctic). Front. Microbiol. 2016, 7, 227. [Google Scholar] [CrossRef]
  16. Zhang, T.; Wang, N.F.; Yu, L.Y. Geographic Distance and Habitat Type Influence Fungal Communities in the Arctic and Antarctic Sites. Microb. Ecol. 2021, 82, 224–232. [Google Scholar] [CrossRef]
  17. Zhang, T.; Wang, N.F.; Zhang, Y.Q.; Liu, H.Y.; Yu, L.Y. Diversity and distribution of fungal communities in the marine sediments of Kongsfjorden, Svalbard (High Arctic). Sci. Rep. 2015, 5, 14524. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, T.; Wei, X.L.; Zhang, Y.Q.; Liu, H.Y.; Yu, L.Y. Diversity and distribution of lichen-associated fungi in the Ny-Alesund Region (Svalbard, High Arctic) as revealed by 454 pyrosequencing. Sci. Rep. 2015, 5, 14850. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, T.; Wang, N.; Yu, L. Host-specificity of moss-associated fungal communities in the Ny-Ålesund region (Svalbard, High Arctic) as revealed by amplicon pyrosequencing. Fungal Ecol. 2021, 53, 101092. [Google Scholar] [CrossRef]
  20. Blaalid, R.; Davey, M.L.; Kauserud, H.; Carlsen, T.; Halvorsen, R.; Hoiland, K.; Eidesen, P.B. Arctic root-associated fungal community composition reflects environmental filtering. Mol. Ecol. 2014, 23, 649–659. [Google Scholar] [CrossRef]
  21. Botnen, S.S.; Thoen, E.; Eidesen, P.B.; Krabberod, A.K.; Kauserud, H. Community composition of arctic root-associated fungi mirrors host plant phylogeny. FEMS Microbiol. Ecol. 2020, 96, fiaa185. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, T.; Yao, Y.F. Endophytic Fungal Communities Associated with Vascular Plants in the High Arctic Zone Are Highly Diverse and Host-Plant Specific. PLoS ONE 2015, 10, e0130051. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hassett, B.T.; Ducluzeau, A.L.; Collins, R.E.; Gradinger, R. Spatial distribution of aquatic marine fungi across the western Arctic and sub-arctic. Environ. Microbiol. 2017, 19, 475–484. [Google Scholar] [CrossRef]
  24. Rämä, T.; Hassett, B.T.; Bubnova, E. Arctic marine fungi: From filaments and flagella to operational taxonomic units and beyond. Bot. Mar. 2017, 60, 433–452. [Google Scholar] [CrossRef]
  25. Kluge, M.; Wauthy, M.; Clemmensen, K.E.; Wurzbacher, C.; Hawkes, J.A.; Einarsdottir, K.; Rautio, M.; Stenlid, J.; Peura, S. Declining fungal diversity in Arctic freshwaters along a permafrost thaw gradient. Glob. Chang. Biol. 2021, 27, 5889–5906. [Google Scholar] [CrossRef]
  26. Zhang, T.; Wang, N.F.; Zhang, Y.Q.; Liu, H.Y.; Yu, L.Y. Diversity and Distribution of Aquatic Fungal Communities in the Ny-Alesund Region, Svalbard (High Arctic): Aquatic Fungi in the Arctic. Microb. Ecol. 2016, 71, 543–554. [Google Scholar] [CrossRef]
  27. Mundra, S.; Halvorsen, R.; Kauserud, H.; Bahram, M.; Tedersoo, L.; Elberling, B.; Cooper, E.J.; Eidesen, P.B. Ectomycorrhizal and saprotrophic fungi respond differently to long-term experimentally increased snow depth in the High Arctic. Microbiologyopen 2016, 5, 856–869. [Google Scholar] [CrossRef]
  28. Canini, F.; Zucconi, L.; Pacelli, C.; Selbmann, L.; Onofri, S.; Geml, J. Vegetation, pH and Water Content as Main Factors for Shaping Fungal Richness, Community Composition and Functional Guilds Distribution in Soils of Western Greenland. Front. Microbiol. 2019, 10, 2348. [Google Scholar] [CrossRef] [Green Version]
  29. Geml, J.; Semenova, T.A.; Morgado, L.N.; Welker, J.M. Changes in composition and abundance of functional groups of arctic fungi in response to long-term summer warming. Biol. Lett. 2016, 12, 20160503. [Google Scholar] [CrossRef] [Green Version]
  30. Hassett, B.T.; Borrego, E.J.; Vonnahme, T.R.; Rama, T.; Kolomiets, M.V.; Gradinger, R. Arctic marine fungi: Biomass, functional genes, and putative ecological roles. ISME J. 2019, 13, 1484–1496. [Google Scholar] [CrossRef] [Green Version]
  31. Perini, L.; Mogrovejo, D.C.; Tomazin, R.; Gostincar, C.; Brill, F.H.H.; Gunde-Cimerman, N. Phenotypes Associated with Pathogenicity: Their Expression in Arctic Fungal Isolates. Microorganisms 2019, 7, 600. [Google Scholar] [CrossRef] [Green Version]
  32. Maturilli, M.; Herber, A.; König-Langlo, G. Climatology and time series of surface meteorology in Ny-Ålesund, Svalbard. Earth Syst. Sci. Data 2013, 5, 155–163. [Google Scholar] [CrossRef] [Green Version]
  33. Gardes, M.; Bruns, T.D. ITS Primers with enhanced specificity for Basidiomycetes—Application to the identification of mycorrhizae and rusts. Mol. Ecol. 1993, 2, 113–118. [Google Scholar] [CrossRef]
  34. White, T.J.; Bruns, T.D.; Lee, S.B.; Taylor, J.W. Amplification and Direct Sequencing of Fungal Ribosomal RNA Genes for Phy-logenetics. In PCR Protocols: A Guide to Methods and Applications; Innis, M.A., Gelfand, D.H., Sninsky, J.J., White, T.J., Eds.; Academic Press: New York, NY, USA, 1990; pp. 315–322. [Google Scholar]
  35. Magoc, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef] [Green Version]
  36. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  37. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  38. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahe, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ 2016, 4, e2584. [Google Scholar] [CrossRef] [Green Version]
  39. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2′s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef]
  40. Abarenkov, K.; Zirk, A.; Piirmann, T.; Pöhönen, R.; Ivanov, F.; Nilsson, R.H.; Kõljalg, U. UNITE QIIME release for Fungi 2. Version 10.05.2021. UNITE Community 2021, 7, 1264763. [Google Scholar] [CrossRef]
  41. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 2019, 20, 1160–1166. [Google Scholar] [CrossRef] [Green Version]
  42. Dhariwal, A.; Chong, J.; Habib, S.; King, I.L.; Agellon, L.B.; Xia, J. MicrobiomeAnalyst: A web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res. 2017, 45, W180–W188. [Google Scholar] [CrossRef] [PubMed]
  43. Põlme, S.; Abarenkov, K.; Nilsson, R.H.; Lindahl, B.; Clemmensen, K.; Kauserud, H.; Nguyen, N.; Kjøller, R.; Bates, S.; Baldrian, P.; et al. FungalTraits: A user-friendly traits database of fungi and fungus-like stramenopiles. Fungal Divers. 2020, 105, 1–16. [Google Scholar] [CrossRef]
  44. De Hoog, G.S.; Guarro, J.; Gené, J.; Ahmed, S.A.; Al-Hatmi, A.M.S.; Figueras, M.J.; Vitale, R.G. Atlas of Clinical Fungi, 4th ed.; CBS: Utrecht, The Netherlands, 2020. [Google Scholar]
  45. Walker, J.; Aldrich-Wolfe, L.; Riffel, A.; Barbare, H.; Simpson, N.; Trowbridge, J.; Jumpponen, A. Diverse Helotiales associated with the roots of three species of Arctic Ericaceae provide no evidence for host specificity. New Phytol. 2011, 191, 515–527. [Google Scholar] [CrossRef] [PubMed]
  46. He, C.Y.; Feng, J.J.; Zhang, T.; Yu, L.Y. Species diversity and floristic analysis of macrofungi in Ny-Ålesund Region, High Arctic. Chin. J. Polar Res. 2022, 34, 139–148. [Google Scholar] [CrossRef]
  47. Baldrian, P.; Valaskova, V. Degradation of cellulose by basidiomycetous fungi. FEMS Microbiol. Rev. 2008, 32, 501–521. [Google Scholar] [CrossRef] [Green Version]
  48. Geml, J.; Timling, I.; Robinson, C.; Lennon, N.; Nusbaum, C.; Brochmann, C.; Noordeloos, M.; Taylor, D.L. An arctic community of symbiotic fungi assembled by long-distance dispersers: Phylogenetic diversity of ectomycorrhizal basidiomycetes in Svalbard based on soil and sporocarp DNA. J. Biogeogr. 2011, 39, 74–88. [Google Scholar] [CrossRef]
  49. Prasad, R.; Chhabra, S.; Gill, S.; Kumar Singh, P.; Tuteja, N. The Microbial Symbionts: Potential for Crop Improvement in Changing Environments. In Advancement in Crop Improvement Techniques; Tuteja, N., Tuteja, R., Passricha, N., Saifi, S.K., Eds.; Woodhead Publishing: London, UK, 2020; pp. 233–240. [Google Scholar] [CrossRef]
  50. Pusz, W.; Urbaniak, J. Airborne fungi in Longyearbyen area (Svalbard, Norway)—Case study. Environ. Monit. Assess. 2021, 193, 290. [Google Scholar] [CrossRef]
  51. Wu, R.; Trubl, G.; Taş, N.; Jansson, J.K. Permafrost as a potential pathogen reservoir. One Earth 2022, 5, 351–360. [Google Scholar] [CrossRef]
Jof 09 00437 g001 550
Figure 1. The location of the sampling site and views of various habitats (e.g., soil, lichen, vascular plant, moss, freshwater, seawater, dung, and marine alga) in the Ny-Ålesund Region (Svalbard, High Arctic).
Figure 1. The location of the sampling site and views of various habitats (e.g., soil, lichen, vascular plant, moss, freshwater, seawater, dung, and marine alga) in the Ny-Ålesund Region (Svalbard, High Arctic).
Jof 09 00437 g001
Jof 09 00437 g002 550
Figure 2. Bar plots showing the abundance of (a) fungal phyla and (b) fungal orders in the 114 samples collected from nine habitats of the Ny-Ålesund Region (Svalbard, High Arctic).
Figure 2. Bar plots showing the abundance of (a) fungal phyla and (b) fungal orders in the 114 samples collected from nine habitats of the Ny-Ålesund Region (Svalbard, High Arctic).
Jof 09 00437 g002
Jof 09 00437 g003 550
Figure 3. (a) NMDS ordination plot showing the spatial pattern of the mycobiome in the 114 samples from the nine habitats, (b) Venn diagram showing the number of ASVs in the nine habitats.
Figure 3. (a) NMDS ordination plot showing the spatial pattern of the mycobiome in the 114 samples from the nine habitats, (b) Venn diagram showing the number of ASVs in the nine habitats.
Jof 09 00437 g003
Jof 09 00437 g004 550
Figure 4. Correlation network analyses showing the correlations of the fungal genera in (a) four habitat types (i.e., freshwater, seawater, sediment, and soil) and (b) five habitat types (i.e., dung, lichen, marine alga, moss, and vascular plant). Each node represents taxa at the genus level and the edges represent correlations between taxa pairs. The nodes are colored based on habitat types.
Figure 4. Correlation network analyses showing the correlations of the fungal genera in (a) four habitat types (i.e., freshwater, seawater, sediment, and soil) and (b) five habitat types (i.e., dung, lichen, marine alga, moss, and vascular plant). Each node represents taxa at the genus level and the edges represent correlations between taxa pairs. The nodes are colored based on habitat types.
Jof 09 00437 g004
Jof 09 00437 g005 550
Figure 5. LEfSe analysis showing (a) the fungal orders and (b) the genera that differed significantly among the nine habitat types. The mini heatmap to the right of the plot indicates whether the taxa are higher (red) or lower (blue) in each group.
Figure 5. LEfSe analysis showing (a) the fungal orders and (b) the genera that differed significantly among the nine habitat types. The mini heatmap to the right of the plot indicates whether the taxa are higher (red) or lower (blue) in each group.
Jof 09 00437 g005
Jof 09 00437 g006 550
Figure 6. LEfSe analysis showing (a) the growth forms and (b) the functional guilds that differed significantly among the nine habitat types. The mini heatmap to the right of the plot indicates whether the taxa are higher (red) or lower (blue) in each group.
Figure 6. LEfSe analysis showing (a) the growth forms and (b) the functional guilds that differed significantly among the nine habitat types. The mini heatmap to the right of the plot indicates whether the taxa are higher (red) or lower (blue) in each group.
Jof 09 00437 g006
Table
Table 1. ANOSIM results for the pairwise comparisons of the similarity among the nine habitats.
Table 1. ANOSIM results for the pairwise comparisons of the similarity among the nine habitats.
Habitat TypeSoilFreshwaterSeawaterLichenSedimentDungMarine algaMossVascular Plant
Soil-0.67664 ***0.84645 ***0.48733 ***0.97115 ***0.97516 ***0.47797 **0.37873 ***0.58004 ***
Freshwater0.67664 ***-0.9078 ***0.77907 ***0.97385 ***1 ***0.76711 **0.65965 ***0.94234 ***
Seawater0.84645 ***0.9078 ***-0.84193 ***0.90845 ***1 ***0.94132 ***0.92049 ***0.99311 ***
Lichen0.48733 ***0.77907 ***0.84193 ***-0.82897 ***0.40703 ***0.43981 **0.50916 ***0.48796 ***
Sediment0.97115 ***0.97385 ***0.90845 ***0.82897 ***-1 ***0.94868 ***0.96906 ***1 ***
Dung0.97516 ***1 ***1 ***0.40703 ***1 ***-0.90614 **0.95968 ***0.92961 ***
Marine alga0.47797 **0.76711 **0.94132 ***0.43981 **0.94868 ***0.90614 **-0.61455 ***0.42719 **
Moss0.37873 ***0.65965 ***0.92049 ***0.50916 ***0.96906 ***0.95968 ***0.61455 ***-0.60594 ***
Vascular plant0.58004 ***0.94234 ***0.99311 ***0.48796 ***1 ***0.92961 ***0.42719 **0.60594 ***-
** 0.001 < p ≤ 0.01, *** p ≤ 0.001.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, X.; Yan, D.; Yu, L.; Zhang, T. An Integrative Study of Mycobiome in Different Habitats from a High Arctic Region: Diversity, Distribution, and Functional Role. J. Fungi 2023, 9, 437. https://doi.org/10.3390/jof9040437

AMA Style

Chen X, Yan D, Yu L, Zhang T. An Integrative Study of Mycobiome in Different Habitats from a High Arctic Region: Diversity, Distribution, and Functional Role. Journal of Fungi. 2023; 9(4):437. https://doi.org/10.3390/jof9040437

Chicago/Turabian Style

Chen, Xiufei, Dong Yan, Liyan Yu, and Tao Zhang. 2023. "An Integrative Study of Mycobiome in Different Habitats from a High Arctic Region: Diversity, Distribution, and Functional Role" Journal of Fungi 9, no. 4: 437. https://doi.org/10.3390/jof9040437

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