Next Article in Journal
Resistance to Boscalid, Fluopyram and Fluxapyroxad in Blumeriella jaapii from Michigan (U.S.A.): Molecular Characterization and Assessment of Practical Resistance in Commercial Cherry Orchards
Previous Article in Journal
Propionibacterium acnes-Derived Circulating Immune Complexes in Sarcoidosis Patients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 11
10.3390/microorganisms9112195
microorganisms-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

Microbial Diversity in Subarctic Biocrusts from West Iceland following an Elevation Gradient

by
Ekaterina Pushkareva
1,2,*,
Israel Barrantes
3,
Peter Leinweber
4 and
Ulf Karsten
1
1
Department of Applied Ecology and Phycology, Institute of Biological Sciences, University of Rostock, 18059 Rostock, Germany
2
Department of Biology, Botanical Institute, University of Cologne, 50674 Cologne, Germany
3
Research Group Translational Bioinformatics, Institute for Biostatistics and Informatics in Medicine and Ageing Research, Rostock University Medical Center, 18057 Rostock, Germany
4
Department Soil Sciences, Faculty of Agricultural and Environmental Sciences, University of Rostock, Justus-von-Liebig-Weg 6, 18051 Rostock, Germany
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(11), 2195; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9112195
Submission received: 30 August 2021 / Revised: 15 October 2021 / Accepted: 17 October 2021 / Published: 21 October 2021
(This article belongs to the Section Environmental Microbiology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Biological soil crusts (biocrusts) are essential communities of organisms in the Icelandic soil ecosystem, as they prevent erosion and cryoturbation and provide nutrients to vascular plants. However, biocrust microbial composition in Iceland remains understudied. To address this gap in knowledge, we applied high-throughput sequencing to study microbial community composition in biocrusts collected along an elevation gradient (11–157 m a.s.l.) stretching away perpendicular to the marine coast. Four groups of organisms were targeted: bacteria and cyanobacteria (16S rRNA gene), fungi (transcribed spacer region), and other eukaryotes (18S rRNA gene). The amplicon sequencing of the 16S rRNA gene revealed the dominance of Proteobacteria, Bacteroidetes, and Actinobacteria. Within the cyanobacteria, filamentous forms from the orders Synechococcales and Oscillatoriales prevailed. Furthermore, fungi in the biocrusts were dominated by Ascomycota, while the majority of reads obtained from sequencing of the 18S rRNA gene belonged to Archaeplastida. In addition, microbial photoautotrophs isolated from the biocrusts were assigned to the cyanobacterial genera Phormidesmis, Microcoleus, Wilmottia, and Oscillatoria and to two microalgal phyla Chlorophyta and Charophyta. In general, the taxonomic diversity of microorganisms in the biocrusts increased following the elevation gradient and community composition differed among the sites, suggesting that microclimatic and soil parameters might shape biocrust microbiota.

Graphical Abstract

1. Introduction

Biological soil crusts (biocrusts) are assemblages of organisms, such as bacteria, cyanobacteria, fungi, eukaryotic microalgae, lichens, and mosses, living on the soil surface. They can survive long periods of droughts and live under extreme conditions with strong temperature and light gradients. Biocrust organisms are poikilohydric and, thus, can be rapidly reactivated once water is available. Moreover, biocrusts play an important role in primary production, nutrient cycling, mineralization, weathering, and soil stabilization [1]. Microbial carbon dioxide (CO2) fixation is mostly provided by photoautotrophic microorganisms, which are prominent members of biocrusts. Cyanobacteria, for example, usually participate in biogeochemical cycles such as carbon (C) or nitrogen (N) fixation [2]. Furthermore, cyanobacteria promote biocrust development forming a cohesive layer with soil particles, which increases surface stability and inhibits erosion [3]. Non-photoautotrophic bacterial communities in biocrusts typically consist of a number of ubiquitous phyla including Actinobacteria, Bacteroidetes, Proteobacteria, and Verrucomicrobia [4,5]. Proteobacteria and Actinobacteria that typically exhibit wide metabolic ranges, competitive advantages, and multiple UV repair mechanisms are functionally important in nutrient-limited arid environments [6]. Eukaryotic microbial members of biocrust communities are usually represented by fungi and microalgae. Fungi, including Basidiomycota and Ascomycota, are ecologically important, as they can release bioavailable nutrients and form symbiotic associations with microbial photoautotrophs into lichens [7]. In general, fungi produce vegetative filaments (hyphae) that bind soil particles together and help to consolidate the surface of biocrust communities. Eukaryotic microalgae, similar to cyanobacteria, contribute to the C cycle through their photosynthesis. Due to the structural and functional differences in photosynthetic organelles, eukaryotic microalgae have higher rates of photosynthesis than cyanobacteria but a lower tolerance toward freeze–thaw cycles [8]. While most knowledge about biocrust community composition and ecological functions originates from drylands, much less is known from higher latitudes [9].
Iceland, located right below the Arctic Circle, is a geologically young island with frequent volcanic eruptions. Icelandic soils, known as Andosols, are derived from volcanic ash and, therefore, are very distinct from other soils around the world [10]. They are susceptible to weathering which, in turn, results in the release of aluminum (Al), iron (Fe), and silicon (Si) [11]. Moreover, thin Andosols together with a harsh climate slow down vascular plant colonization in Iceland and, hence, biocrusts occupy considerable areas located on free spaces between the higher plants [12].
Despite the importance of biocrusts in high latitude areas, particularly Iceland, there are a limited number of studies about their microbial community composition. Most of the available literature focuses on the microbial communities in the Icelandic soil environment in general. For example, bacteria were characterized in volcanic deposits in Northern Iceland using high-throughput sequencing (HTS) which revealed the dominance of Proteobacteria, Acidobacteria, Actinobacteria, and Firmicutes [13]. Likewise, HTS was applied to study bacterial diversity in the deep soil (below 145 m depth) of the volcanic island Surtsey, and the majority of reads belonged to Methanobacteriales and Archaeoglobales [14]. Furthermore, arbuscular mycorrhizal fungi [15,16] and other terrestrial fungi [17,18] were described in Icelandic soils using morphological methods. Similarly, the diversity of cyanobacteria and eukaryotic microalgae was only reported in papers from the last century [19,20]. However, our recent study [12] based on a morphological approach revealed unexpectedly rich communities of microbial photoautotrophs in Icelandic biocrusts. Eukaryotic microalgae, such as green algae (Chlorella, Coccomyxa, and Stichococcus) and diatoms (e.g., Pinnularia and Eunotia) dominated these biocrusts [12].
Given the lack of information on microbial community composition inhabiting Icelandic biocrusts, the objective of this study was to obtain an overview of bacteria, cyanobacteria, fungi, and other eukaryotes in the biocrusts of west Iceland using HTS. Moreover, we hypothesized that different elevations (from 11 to 156 m a.s.l.) and distance from the sea (from 0 to 80 km) as well as soil properties would affect microbiota in these biocrusts.

2. Materials and Methods

2.1. Site Description

Litla-Skard (64°43’36” N, 21°37’48” W) is a research station in the west of Iceland approximately 100 km farther north from Reykjavik and belongs to the National Forestry Service (Icelandic Forest Service). The area has a subarctic climate with an annual average temperature of 3.1 °C. The coldest and warmest months are February and July with mean temperatures of −1.8 and 10.8 °C, respectively. The average annual precipitation is 930 mm in the form of rain and 100 mm in the form of snow (meteorological data are provided from the research station Litla-Skard; no online data are available). The vegetation is represented by shrub birches, moss heaths, marsh grass, and grassland. The soils in this area are Andosols, which have coarse texture, high porosity, low cohesion, and the dispersal of many layers of tephra caused by frequent volcanic eruptions.
Five sampling sites were chosen for this study that represented a catena of approximately 80 km in distance from the sea into inland and were located at different elevations from 11 to 157 m a.s.l.: Krákunes (KRA; 11 m a.s.l., at the sea coast), Borgarfjarðarbraut (BOR; 25 m a.s.l., 60 km from the sea), Fiflholt (FIF; 46 m a.s.l., 11 km from the sea), Litla-Skard station (LSK; 106 m a.s.l., 50 km from the sea), and Giljar (GIL; 157 m a.s.l., 80 km from the sea).

2.2. Sampling and Soil Characterization

Sampling was conducted in July 2014. Biocrust samples (up to 1 cm in depth) were taken in 4 replicates from each site and kept moist and cool until they were delivered to the laboratory where they were stored in a freezer.
The description of the samples (soil chemistry, sequential phosphorus fractionation, microbial P, and potential activity of phosphatase as well as community composition of algae and fungi assessed by microscope) are reported in Pushkareva et al. [12]. In brief, the biocrusts were acidic (pH 5.5–5.8) with a high P content (982–1571 mg kg−1) and moderately labile P was dominant. Microphotoautotrophic communities were dominated by eukaryotic microalgae including diatoms. Green algae, such as Chlorella vulgaris, Chlamydomonas sp., and Coccomyxa sp., and diatoms, such as Pinnularia sp. and Eunotia sp., were present in all sites. The majority of cyanobacteria belonged to filamentous forms such as Microcoleus sp., Leptolyngbya sp., Phormidesmis sp., and Stenomitos sp. In addition, the majority of fungal strains isolated from the biocrusts were assigned to Penicillium sp. and Fusarium sp. [12]. Moreover, mosses and lichens were not observed in the samples.

2.3. Cultivation and Molecular Analyses of Microbial Photoautotrophs

Cyanobacteria and eukaryotic microalgae were isolated from the enrichment cultures used in our previous study [12]. Single colonies were transferred to the Petri dishes with Bold’s Basal Medium, and unialgal cultures were established. The cultures were then kept at 15 °C under 30 μmol photons m−2 s−1 (Osram Lumilux Cool White lamps L36W/840) with a light/dark regime of 16/8 h.
The DNA of cyanobacterial and microalgal strains was extracted using the NucleoSpin Plant II mini kit (Macherey Nagel, Düren, Germany) according to the manufacturer’s instructions. The 16S (for cyanobacteria) and 18S rRNA (for eukaryotic microalgae) genes were amplified using primers BS-1F/CPL-10R and EAF3/ITS055R, respectively (Supplementary Materials Table S1). PCRs for both sets of primers were performed as follows: an initial denaturation step at 95 °C for 3 min, followed by 35 cycles of DNA denaturation at 95 °C for 15 s, primer annealing at 55 °C for 15 s, strand extension at 72 °C for 1 min 30 s, and a final extension step at 72 °C for 2 min. Sanger sequencing was performed at GATC Sequencing Services (Eurofins Genomics Germany, Ebersberg, Germany) with the primers listed in the Supplementary Materials Table S1.

2.4. Environmental DNA Extraction and Sequencing

Total DNA was extracted from the biocrust samples using the PowerSoil DNA Isolation Kit (MOBIO, Carlsbad, CA, USA) according to the manufacturer’s instructions. DNA concentrations were quantified using a Qubit3.0 Fluorometer according to the manufacturer’s protocol. Extracted DNAs were sent to the Microsynth AG (Balgach, Switzerland), where PCR and sequencing using the Illumina MiSeq platform were performed.
Different sets of primers were used to target four groups of organisms: bacteria, cyanobacteria, fungi, and other eukaryotes (Supplementary Materials Table S1). Universal bacterial primers (341F and 802R) and cyanobacteria specific primers (CYA_F and CYA_R) were used to amplify the V3–V4 region of the 16S rRNA gene, while primers tarEuk_F and tarEuk_F were used for the 18S rRNA gene amplification. To assess fungal community composition, the internal transcribed spacer region (ITS) was amplified with primers Fungi_F and Fungi_R. The raw reads were submitted to the European Nucleotide Archive (ENA) under the project PRJEB45587.

2.5. Bioinformatic and Statistical Analyses

Sequence contigs of isolates were generated using the software Geneious 8.1.9 (Biomatters Ltd., Auckland, New Zealand) and then deposited in GenBank with accession numbers MZ020203–MZ020221 for eukaryotic microalgae and MZ020189–MZ020202 for cyanobacteria. The sequences of the most closely related isolates from GenBank were retrieved using BLAST.
The amplicons were reconstructed from the Illumina sequencing runs using the pandaseq program (version 2.11; [21]). Operational taxonomic units (OTUs) were then identified from these amplicons using USEARCH (version 6.1.544; [22]) with the pick_open_reference_otus.py script of QIIME 1.9.1 [23]. The following databases were used: 16S rRNA Greengenes (version 13.8; [24]) for the bacteria and cyanobacteria; UNITE (version 12_11; [25]) for the ITS amplicons; and SILVA (version 132; [26]) for the 18S rRNA data. In all cases, a cutoff of 97% identity was applied. Afterwards, low confidence OTUs were excluded via the remove_low_confidence_otus.py script [27]. OTUs classified as chloroplasts and mitochondria were removed from the bacterial data set, while OTUs assigned to bacteria were deleted from the 18S rRNA results. Likewise, non-cyanobacterial OTUs were removed from the data set obtained with cyanobacteria-specific primers with additional manual inspection.
One replicate of site KRA (KRA-2) revealed only 31 reads with cyanobacteria-specific primers and, thus, was removed from further analyses. Likewise, this sample did not produce any reads using primers for 18S rRNA amplification. Moreover, two replicates were not included to alpha and beta diversity analyses of eukaryotes (18S rRNA data set) due to the low number of reads (LSK-2 and GIL-1 with 255 and 814 total reads, respectively).
Alpha diversity indices (OTU richness and Shannon’s diversity) were calculated in R (version 1.3.1073) using the package vegan [28]. The differences in parameters among sampling sites were tested with one-way analysis of variance (ANOVA) and Tukey’s HSD post hoc test using the software JMP 13.0.1 (SAS Institute Inc., Cary, NC, USA). A p-value < 0.05 was considered significant in all statistical tests. Normality of variance was assessed using the Shapiro–Wilk’s test. If necessary, data were Log or Box–Cox transformed. The relationships between diversity indices and soil parameters were tested by Pearson’s correlation coefficient (PCC). Community dissimilarities among the sites based on number of reads obtained by amplicon sequencing were assessed by non-metrical multidimensional scaling (NMDS) using the package vegan in R. In addition, permutational multivariate analysis of variance (PERMANOVA) was performed, and homogeneity of variances was tested by the betadisper command prior this analysis. Various taxonomic levels were used for both multivariate analyses: bacterial phyla, cyanobacterial orders, fungal phyla, and eukaryotic clades. Furthermore, Venn diagrams for each group of organisms were created using an online platform (http://bioinformatics.psb.ugent.be/webtools/Venn, accessed on 15 November 2020).
In addition, using FUNGuild database, fungal OTUs were assigned to functional (trophic) groups according to their nutrient acquisition strategy of the fungus: pathotrophs, saprotrophs, symbiotrophs, or a combination of these [29]. Furthermore, a PICRUSt analysis was carried out to predict the relative abundance of functional genes and pathways in the communities (version 2.11; [30]) for bacteria and cyanobacteria. For this, the OTUs were re-assigned with a closed-reference OTU picking protocol (from QIIME 1.9.1) and the Greengenes database (version 13.5), preclustered at 97% identity. The obtained OTU tables were then normalized by a 16S rRNA copy number, and functional genes were predicted from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [31]. Finally, the results from the PICRUSt analysis and KEGG predictions were processed within STAMP (version 2.1.3; [32]).

3. Results

3.1. Overall Diversity of Microorganisms in the Icelandic Biocrusts

3.1.1. Bacteria

Next generation sequencing with general bacterial primers resulted in 369,016 bacterial reads for 20 samples, which were clustered into 3155 OTUs. Of these, the majority of reads and OTUs belonged to the orders Proteobacteria (38% of total reads, 999 OTUs).
Bacteroidetes (15% of total reads, 368 OTUs), Actinobacteria (12% of total reads, 395 OTUs), Planctomycetes (9% of total reads, 396 OTUs), Acidobacteria (7% of total reads, 299 OTUs), Verrucomicrobia (6% of total reads, 242 OTUs), Chloroflexi (4% of total reads, 150 OTUs), Gemmatimonadetes (3% of total reads, 75 OTUs), and Cyanobacteria (1% of total reads, 42 OTUs) (Figure 1a). Within Proteobacteria, Alphaproteobacteria was the dominant class and constituted 18% of total bacterial reads (475 OTUs).
In addition, community function analysis revealed only one pathway in the 16S data set which significantly (p < 0.05) differed among the studied sites (Supplementary Materials Figure S1). A total of 74 bacterial genera were related to the glutathione metabolism pathway (K00480), and the majority of sequences assigned to it were recorded at site BOR.

3.1.2. Cyanobacteria

Sequencing of the 16S rRNA gene with cyanobacteria-specific primers produced 835,441 reads (603 OTUs). After removal of non-cyanobacterial OTUs (other bacteria and chloroplasts), 244 OTUs (527,772 reads) assigned to cyanobacterial taxa remained. Six orders were distinguished according to the Greengenes database: Synechococcales (63% of total reads, 129 OTUs), Oscillatoriales (20% of total reads, 37 OTUs), Nostocales (13% of total reads, 53 OTUs), Chroococcales (2% of total reads, 12 OTUs), Gloeobacterales (<1% of total reads, 5 OTUs), and Stogonematales (<1% of total reads, 1 OTU) (Figure 1b). Seven OTUs (1% of total reads) could not be assigned to any order.
Ten strains of filamentous cyanobacteria isolated from the biocrusts were successfully sequenced (Table 1a). Sequences of six isolates were assigned to the order Oscillatoriales, particularly to the genera Microcoleus and Wilmottia (>94% identity). Further, sequences of four other isolates corresponded to Phormidesmis sp. from the order Synechococcales (>93% identity).

3.1.3. Fungi

A total of 364,233 reads and 624 OTUs assigned to Fungi were obtained by sequencing the ITS region. Five fungal phyla were identified: Ascomycota (37% of total reads, 294 OTUs), Basidiomycota (17% of total reads, 89 OTUs), Zygomycota (7% of total reads, 10 OTUs), Glomeromycota (<1% of total reads, 3 OTUs), and Chytridiomycota (<1% of total reads, 2 OTUs) (Figure 1c). However, a large number of OTUs (38% of total reads, 226 OTUs) could not be assigned to any phyla. The majority of identified Ascomycota belonged to the classes Leotiomycetes, Sordariomycetes, and Eurotiomycetes (8%, 6%, and 5% of total fungal reads, respectively). Sequences affiliated with the class Agaricomycetes comprised the majority of those from Basidiomycota with 11% of total ITS reads.
Furthermore, 67% of fungal OTUs were identified by the FUNGuild database, and only 22% of identified OTUs were assigned to different functional groups. The majority of fungal reads belonged to symbiotrophs and saprotrophs (Figure 2).
In addition, fungal OTUs (132 OTUs) obtained by amplicon sequencing of the 18S rRNA constituted 10% of total eukaryotic reads and 58% of Opisthokonta clade. Phyla Mucoromycota and Cryptomycota exhibited the majority of reads within fungi (38 and 30% of total fungal reads, respectively).

3.1.4. Eukaryotes

Sequencing of the 18S rRNA gene produced 850,183 reads and 1033 OTUs assigned to Eukaryota. The majority of reads affiliated with the clade Archaeplastida (57% of total reads, 325 OTUs) followed by clades SAR (23% of total reads, 425 OTUs), Opisthokonta (17% of total reads, 224 OTUs), Amoebozoa (3% of total reads, 56 OTUs), Centrohelida (<1% of total reads, 56 OTUs) and Excavata (<1% of total reads, 2 OTUs) (Figure 1d). Archaeplastida was represented by the clade Chloroplastida only which, in turn, included two algal phyla, Charophyta (80 OTUs) and Chlorophyta (245 OTUs). The SAR groups were almost evenly distributed within the clade: 30% (105 OTUs) affiliated with Stramenopiles, 34 % (109 OTUs) with Alveolata, and 36% (211 OTUs) with Rhizaria. The majority of reads within Stramenopiles were assigned to Bacillariophyceae (33% of total Stramenopiles reads; 31 OTUs).
In addition, 20 strains of eukaryotic microalgae were isolated and sequenced as previously described (Table 1b). Of these, the majority of the strains were from the order Chlorophyta, and only one strain isolated from site FIF affiliated with the order Charophyta (class Klebsormidiophyceae). Sites KRA and GIL had more distinct isolates in comparison to other sites and each other. For example, Gonium sp. was isolated only from site KRA (KRA-4a), while Stichococcus sp. (GIL-3a), Nannochloris sp. (GIL-3d), and Parietochloris sp. (GIL-4a) only from site GIL. The genera Xerochlorella and Chlorella were common for sites BOR, FIF, and LSK (Table 1b).

3.2. Alpha and Beta Diversities

In general, two diversity indices (OTU richness and Shannon index) showed a trend to increase following the elevation gradient from the sea, with some exceptions (Figure 3). A significant difference between the lowest and highest topographic elevations was observed in the diversity indices for bacteria, cyanobacteria, and fungi. However, a significantly higher OTU richness of eukaryotes and fungi was observed at site FIF. For the cyanobacteria, site BOR had significantly lower OTU richness and Shannon index than other sites. Further, community structure at the BOR site appeared distinct from those at the other sites (Figure 1). For example, significantly more reads from BOR affiliated with Proteobacteria, cyanobacterial order Synechococcales, and clade Archaeplastida than in biocrusts from the other sites. Similarly, most reads from the BOR site did not affiliate with a fungal phylum, and the highest relative abundance of pathotrophs was found here (Figure 2). In addition, Venn diagrams showed that more OTUs were shared among sites KRA, FIF, LSK, and GIL rather than with site BOR (Supplementary Materials Figure S2).
The NMDS plot revealed that microbial communities were distinct among the five sites and overall grouped by the site of origin, except site KRA, which exhibited a rather heterogeneous community composition (Figure 4). This was confirmed by PERMANOVA analysis (Table 2 and Table S2). Moreover, community dissimilarities for each microbial group were explained by different environmental parameters, measured in our previous publication [12]. For example, significant correlations (p < 0.05) were observed between manganese (Mn), Fe, and Al contents and fungal community composition, while dissimilarities in the composition of bacterial communities were explained by C, N, and S contents, C:N ratio, and the potential activity of acid phosphatase. Furthermore, differences in cyanobacterial and eukaryotic community composition were related to Na and microbial P (only for eukaryotes) contents. In addition, Pearson’s correlation test revealed no significant correlations between alpha diversity indices and soil chemistry parameters.

4. Discussion

4.1. Overall Diversity of Microbiota in Icelandic Biocrusts

Most available studies about microbial biodiversity in biocrusts of the Northern Hemisphere are from Arctic regions (e.g., [33,34]). These investigations often reveal a high diversity of microorganisms in polar biocrusts and indicate the striking ecological importance of their microbiota in such high latitude areas. The subarctic zone, located right below the Arctic Circle, has received far less attention. Overall, the bacterial community in the studied Icelandic biocrusts comprised largely Proteobacteria, Bacteroidetes, and Actinobacteria (Figure 1a), which is consistent with previous observations from Arctic and other subarctic environments [35]. These bacteria play an important role in biocrust formation and maintenance due to the fact of their production of exopolysaccharides (Bacteroidetes), pigments (Proteobacteria), or mycelia (Actinobacteria) [36]. They are often observed in permafrost habitats and might be well adapted to cold conditions [37,38]. Moreover, soils with high organic matter content promote growth of Proteobacteria [39]. Andosols, derived from volcanic parent material, accumulate more organic matter than other soil types [40]. Other important bacteria, Planctomycetes, observed in the Icelandic biocrusts and also previously in Arctic soils, strongly participate in the nitrogen cycle by performing an anaerobic ammonium oxidation [35] and, therefore, might play an important role in biocrusts particularly from cold environments.
The assessment of cyanobacterial community composition using amplicon sequencing showed that, in general, the distribution of cyanobacterial orders in Icelandic biocrusts was similar to those from Svalbard archipelago using the same set of primers (Figure 1b; [41,42]). The biocrusts were dominated by filamentous forms from the order Synechococcales, and that is also in agreement with microscopic observations [12] and culture isolation (Table 1a). Moreover, Leptolyngbya sp., which was the major representative of the Synechococcales, is common in biocrusts of cold environments [3]. Interestingly, the strains of filamentous cyanobacteria from this order isolated from the Icelandic biocrusts were assigned to Phormidesmis sp. (>97% identity) and not Leptolyngbya sp. (Table 1a). This might indicate that some sequences stored in the database used for the identification of the reads obtained by amplicon sequencing have incorrect taxonomic names on the genus level. In addition, most of these strains were closely related to Phormidesmis arctica reported previously from Central Svalbard [43] and, furthermore, corresponded to the higher elevations of the Icelandic and Svalbard sampling sites. For example, a strain isolated from the Icelandic site GIL (GIL-3b) was similar to Phormidesmis arctica recorded in Mumien Peak (Svalbard) at the altitude of 442 m a.s.l. and strains from Icelandic sites BOR (BOR-4c) and FIF (FIF-1a) corresponded to one strain isolated from a glacier foreland at the lower altitude of Svalbard [43]. Such data reflect a high degree of genetic relatedness and, hence, it is possible to speculate that, for example, migrating birds, such as Arctic terns (Sterna paradisea), could act as a vector for transportation of cyanobacteria between Iceland and Svalbard. A considerable number of reads in the biocrusts obtained with amplicon sequencing were assigned to another filamentous order of cyanobacteria (Oscillatoriales), and the majority of cultured isolates belonged indeed to the Oscillatoriales. Similar to Synechococcales, members of this order can survive extreme conditions due to the array of adaptive strategies such as dormancy, motility, formation of associations, and production of mucilage [44].
The dominance of the fungal phyla Ascomycota, Basidiomycota, and Zygomycota has been described in different Arctic soils [35], and it was also observed in this study (Figure 1c). Several members of these phyla have evolved ecophysiological traits that permit their survival in cold environments such as production of antifreeze proteins (AFPs) as well as protective sugars, polyols, lipids, and fatty acids [45]. In addition, the high number of OTUs that could not be assigned to any phyla might indicate yet undescribed fungal species in Icelandic biocrusts. Therefore, more studies focusing on fungal isolation and following morphological and molecular determinations are required to confirm this fact. The widely used fungal database UNITE often identifies fungal sequences to a kingdom level only [46], which might also lead to a high number of non-assigned OTUs. In addition, higher relative abundance of saprotrophs occurred in the sites located further from the sea (i.e., BOR, LSK and GIL), which could be related to the higher amount of litter in the inland. These fungi play an important role in decomposition and stabilization of soil organic matter [47].
Archaeplastida, represented here only by Chlorophyta and Charophyta, numerically dominated the eukaryotic communities in terms of the number of sequences in each biocrust (Figure 1d). In agreement with the amplicon sequencing, the strains isolated from the biocrust samples belonged to these two phyla (Table 1b), which are very common in northern latitudes [3]. Moreover, all sequences of cultured isolates were found in the amplicon sequencing data set but were underrepresented in number, pointing to a large part of the microalgae remaining uncultivable under the conditions we used in this study. Diversity and abundance of green algae in polar biocrusts usually exceed those of yellow-green algae [48], which seems to be similar in Icelandic biocrusts. In addition, only one isolate and three OTUs of Klebsormidium sp. were reported in the Icelandic biocrusts. Despite their ability to live in cold environments due to the increased desiccation tolerance [49], they were also detected in low abundance in Svalbard biocrusts [4]. Diatoms constituted only 2% of total eukaryotic reads, which did not correspond to the results of microscopic observations that showed their dominance in the biocrusts. This could be due to the low primer coverage efficiency towards different eukaryotic groups. To the best of our knowledge, there are no primers available with 100% coverage of all eukaryotic microorganisms. In addition, the phylum Cercozoa (Rhizaria), a typical soil protist group, was the dominant representative of the SAR clade in the biocrusts. These single-celled eukaryotes have already been described here, with members of the order Glissomonadida dominating [50].

4.2. Shift of Microbial Community Composition in Biocrust Depending on Environmental Parameters

The microbial diversity in Icelandic biocrusts increased from the lowest to the highest topographic elevations (Figure 3), and several previous studies reached similar observations [42,51,52]. Perhaps, there is more water available in the soil at higher altitudes, which would promote microbial growth [53]. Furthermore, the samples grouping according to the site of origin (except site KRA) showed overall community homogeneity at each elevation (Figure 4). Site KRA located at the coast, however, exhibited a more heterogeneous community composition in comparison with all other sites, and this is in agreement with microscopic observations of microbial photoautotrophs [12]. Moreover, there was a lower bacterial diversity than at the other sites, which could be explained by frequent disturbance induced by marine water or sea spray. This supports a previous finding that salt in the coastal area negatively impacted soil microbial community composition [52]. Moreover, an increase of Bacteroidetes at site KRA is consistent with previous studies showing a positive correlation between electrical conductivity and abundance of Bacteroidetes in soils [54]. Likewise, the increase in the fungi Ascomycota at site KRA is similar to a previous study from India indicating a higher number of this phylum in sand dunes than in other soil environments [55]. In addition, sequences of two microalgal strains isolated from the site KRA shared >99% sequence identity with Chlamydomonas hedleyi, a typical marine species, and Gonium pectorale, a freshwater species (Table 1a). The latter is usually found in nutrient-rich lakes or ponds [56] but to the best of our knowledge has never been described in biocrusts.
Interestingly, community composition at site BOR, located further from the sea but at lower elevation, was distinct from other sites (Figure 1). This site showed a higher abundance of typical biocrust representatives such as Proteobacteria and filamentous forms of Synechococcales [33,41,52] than the other sites. Furthermore, the prevalent clade Archaeplastida was dominated by algae from the phylum Charophyta at site BOR unlike the other sites, where members of the phylum Chlorophyta prevailed. Perhaps, distance from the sea and low altitude played a role in community shaping rather than the chemical composition, as there were no prominent statistical correlations between alpha diversity indices and soil chemistry data. Nevertheless, there were positive correlations between two prokaryotic groups and fungi.
In addition, a glutathione metabolism pathway was revealed in several bacterial genera, and most of them were recorded at site BOR (Supplementary Materials Figure S1). Glutathione is involved in bacterial redox-regulation and acclimation to various stressors [57], which could explain higher proportions of sequences responsible for glutathione metabolism at BOR, with the lowest carbon and nitrogen contents and C/N ratio [12].

5. Conclusions

In this study, we showed that the biocrusts from west Iceland harbor diverse communities of various prokaryotic and eukaryotic microorganisms that differ in their composition depending on their microenvironment. Furthermore, we observed an increase in microbial diversity with increasing elevation.
This work provides a biodiversity baseline for understanding the structure and functions of biocrusts in this unique subarctic environment.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms9112195/s1, Figure S1: Proportion of bacterial sequences potentially involved in glutathione metabolism estimated by matching of 16S rRNA gene amplicons with genomes from reference database (KEGG), Figure S2: Venn diagram showing the relationships between the studied sites based on shared OTUs of bacteria (a), cyanobacteria (b), fungi (c) and eukaryotes (d), Table S1: List of primers used for the (a) amplicon and (b) sanger sequencing, Table S2: Differences of microbial community composition between the five studied sites based on pairwise PERMANOVA (9999 permutations).

Author Contributions

Data curation, E.P. and I.B.; Funding acquisition, P.L. and U.K.; Investigation, E.P.; Resources, P.L.; Software, I.B.; Supervision, U.K.; Visualization, E.P.; Writing—original draft, E.P. All authors have read and agreed to the published version of the manuscript.

Funding

P.L. was funded by the European Union Seventh Framework Program (FP7/2007–2013) under grant agreement no. 262693 (INTERACT). E.P. and U.K. were supported by the BiodivERsA CLIMARCTIC project (KA899/33-1) funded by DFG (Germany).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Sequence contigs of isolates were deposited in GenBank with accession numbers MZ020203–MZ020221 for eukaryotic microalgae and MZ020189–MZ020202 for cyanobacteria. The raw reads from amplicon sequencing were submitted to the European Nucleotide Archive (ENA) under the project PRJEB45587.

Acknowledgments

We would like to acknowledge the staff of the research station Litla-Skard (Iceland) for their kind support during the sampling. We also thank Martin Albrecht for his assistance in the molecular work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Belnap, J.; Lange, O.L. Biological Soil Crusts: Structure, Function, and Management; Springer: Berlin, Germany, 2003; ISBN 9783540437574. [Google Scholar]
  2. Stewart, K.J.; Lamb, E.G.; Coxson, D.S.; Siciliano, S.D. Bryophyte-cyanobacterial associations as a key factor in N2-fixation across the Canadian Arctic. Plant Soil 2011, 344, 335–346. [Google Scholar] [CrossRef]
  3. Büdel, B.; Dulić, T.; Darienko, T.; Rybalka, N.; Friedl, T. Cyanobacteria and Algae of Biological Soil Crusts. In Biological Soil Crusts: An Organizing Principle in Drylands. Ecological Studies (Analysis and Synthesis); Weber, B., Büdel, B., Belnap, J., Eds.; Springer: Cham, Switzerland, 2016; Volume 226, pp. 55–80. [Google Scholar]
  4. Rippin, M.; Borchhardt, N.; Williams, L.; Colesie, C.; Jung, P.; Büdel, B.; Karsten, U.; Becker, B. Genus richness of microalgae and Cyanobacteria in biological soil crusts from Svalbard and Livingston Island: Morphological versus molecular approaches. Polar Biol. 2018, 41, 909–923. [Google Scholar] [CrossRef]
  5. Steven, B.; Lionard, M.; Kuske, C.R.; Vincent, W.F. High Bacterial Diversity of Biological Soil Crusts in Water Tracks over Permafrost in the High Arctic Polar Desert. PLoS ONE 2013, 8, e71489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rosenberg, E.; DeLong, E.F.; Lory, S.; Stackebrandt, E.; Thompson, F. (Eds.) The Prokaryotes: Prokaryotic Communities and Ecophysiology; Springer: Berlin/Heidelberg, Germany, 2012; pp. 1–528. [Google Scholar] [CrossRef]
  7. Makhalanyane, T.P.; Valverde, A.; Velázquez, D.; Gunnigle, E.; Van Goethem, M.W.; Quesada, A.; Cowan, D.A. Ecology and biogeochemistry of cyanobacteria in soils, permafrost, aquatic and cryptic polar habitats. Biodivers. Conserv. 2015, 24, 819–840. [Google Scholar] [CrossRef]
  8. Pichrtová, M.; Kulichová, J.; Holzinger, A. Nitrogen limitation and slow drying induce desiccation tolerance in conjugating green algae (Zygnematophyceae, Streptophyta) from polar habitats. PLoS ONE 2014, 9, e113137. [Google Scholar] [CrossRef]
  9. Weber, B.; Büdel, B.; Belnap, J. Biological Soil Crusts: An Organizing Principle in Drylands; Springer International Publishing: Cham, Switzerland, 2016. [Google Scholar]
  10. Arnalds, Ó. Icelandic soils. Dev. Quat. Sci. 2005, 5, 309–318. [Google Scholar] [CrossRef]
  11. Opfergelt, S.; Williams, H.M.; Cornelis, J.T.; Guicharnaud, R.A.; Georg, R.B.; Siebert, C.; Gislason, S.R.; Halliday, A.N.; Burton, K.W. Iron and silicon isotope behaviour accompanying weathering in Icelandic soils, and the implications for iron export from peatlands. Geochim. Cosmochim. Acta 2017, 217, 273–291. [Google Scholar] [CrossRef]
  12. Pushkareva, E.; Baumann, K.; Van, A.T.; Mikhailyuk, T.; Baum, C.; Hrynkiewicz, K.; Demchenko, E.; Thiem, D.; Köpcke, T.; Karsten, U.; et al. Diversity of microbial phototrophs and heterotrophs in Icelandic biocrusts and their role in phosphorus-rich Andosols. Geoderma 2021, 386. [Google Scholar] [CrossRef]
  13. Byloos, B.; Monsieurs, P.; Mysara, M.; Leys, N.; Boon, N.; Van Houdt, R. Characterization of the bacterial communities on recent Icelandic volcanic deposits of different ages 04 Earth Sciences 0403 Geology. BMC Microbiol. 2018, 18, 1–11. [Google Scholar] [CrossRef]
  14. Marteinsson, V.; Klonowski, A.; Reynisson, E.; Vannier, P.; Sigurdsson, B.D.; Ólafsson, M. Microbial colonization in diverse surface soil types in Surtsey and diversity analysis of its subsurface microbiota. Biogeosciences 2015, 12, 1191–1203. [Google Scholar] [CrossRef] [Green Version]
  15. Greipsson, S.; El-Mayas, H.; Vestberg, M.; Walker, C. Arbuscular mycorrhizal fungi in sandy soils in Iceland. Arct. Antarct. Alp. Res. 2002, 34, 419–427. [Google Scholar] [CrossRef]
  16. Greipsson, S.; El-Mayas, H. Arbuscular mycorrhizae of Leymus arenarius on coastal sands and reclamation sites in Iceland and response to inoculation. Restor. Ecol. 2000, 8, 144–150. [Google Scholar] [CrossRef]
  17. Eyjolfsdottir, G.G. Soil fungi isolated from Icelandic farmland. Acta Bot. Isl. 1995, 12, 53–62. [Google Scholar]
  18. Henriksson, L.E.; Henriksson, E. Occurrence of fungi on the volcanic island of Surtsey, Iceland. Acta Bot. Isl. 1974, 3, 82–88. [Google Scholar]
  19. Broady, P.A. The terrestrial algae of Glerardalur Akureyri, Iceland. Acta Bot. Islandica 1978, 5, 3–60. [Google Scholar]
  20. Broady, P. Green and yellow-green terrestrial algae from Surtsey (Iceland) in 1978. Surtsey Res. Prog. Rep. 1982, 9, 13–32. [Google Scholar]
  21. Masella, A.P.; Bartram, A.K.; Truszkowski, J.M.; Brown, D.G.; Neufeld, J.D. PANDAseq: Paired-end assembler for illumina sequences. BMC Bioinform. 2012, 13, 31. [Google Scholar] [CrossRef] [Green Version]
  22. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  23. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high- throughput community sequencing data. Nat. Publ. Gr. 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  24. McDonald, D.; Price, M.N.; Goodrich, J.; Nawrocki, E.P.; DeSantis, T.Z.; Probst, A.; Andersen, G.L.; Knight, R.; Hugenholtz, P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6, 610–618. [Google Scholar] [CrossRef]
  25. Abarenkov, K.; Nilsson, R.H.; Larsson, K.H.; Alexander, I.J.; Eberhardt, U.; Erland, S.; Høiland, K.; Kjøller, R.; Larsson, E.; Pennanen, T.; et al. The UNITE database for molecular identification of fungi-recent updates and future perspectives. New Phytol. 2010, 186, 281–285. [Google Scholar] [CrossRef]
  26. 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. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  27. Comeau, A.M.; Douglas, G.M.; Langille, M.G.I. Microbiome Helper: A Custom and Streamlined Workflow for Microbiome Research. mSystems 2017, 2, 1–11. [Google Scholar] [CrossRef] [Green Version]
  28. Oksanen, J. Multivariate analysis of ecological communities in R:vegan tutorial. R Doc. 2013, 1–43. [Google Scholar] [CrossRef]
  29. Nguyen, N.H.; Song, Z.; Bates, S.T.; Branco, S.; Tedersoo, L.; Menke, J.; Schilling, J.S.; Kennedy, P.G. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 2016, 20, 241–248. [Google Scholar] [CrossRef]
  30. Langille, M.G.I.; Zaneveld, J.; Caporaso, J.G.; McDonald, D.; Knights, D.; Reyes, J.A.; Clemente, J.C.; Burkepile, D.E.; Vega Thurber, R.L.; Knight, R.; et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 2013, 31, 814–821. [Google Scholar] [CrossRef] [PubMed]
  31. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36, D480–D484. [Google Scholar] [CrossRef]
  32. Parks, D.H.; Tyson, G.W.; Hugenholtz, P.; Beiko, R.G. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics 2014, 30, 3123–3124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Rippin, M.; Lange, S.; Sausen, N.; Becker, B. Biodiversity of biological soil crusts from the Polar Regions revealed by metabarcoding. FEMS Microbiol. Ecol. 2018, 94, 36. [Google Scholar] [CrossRef] [PubMed]
  34. McCann, C.M.; Wade, M.; Gray, N.D.; Roberts, J.A.; Hubert, C.R.J.; Graham, D.W. Microbial communities in a High Arctic polar desert landscape. Front. Microbiol. 2016, 7, Art 419. [Google Scholar] [CrossRef]
  35. Malard, L.A.; Pearce, D.A. Minireview Microbial diversity and biogeography in Arctic soils. Environ. Microbiol. Rep. 2018, 10, 611–625. [Google Scholar] [CrossRef]
  36. Gundlapally, S.R.; Garcia-Pichel, F. The community and phylogenetic diversity of biological soil crusts in the Colorado Plateau studied by molecular fingerprinting and intensive cultivation. Microb. Ecol. 2006, 52, 345–357. [Google Scholar] [CrossRef] [PubMed]
  37. Mackelprang, R.; Waldrop, M.P.; DeAngelis, K.M.; David, M.M.; Chavarria, K.L.; Blazewicz, S.J.; Rubin, E.M.; Jansson, J.K. Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw. Nature 2011, 480, 368–371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Frey, B.; Rime, T.; Phillips, M.; Stierli, B.; Hajdas, I.; Widmer, F.; Hartmann, M. Microbial diversity in European alpine permafrost and active layers. FEMS Microbiol. Ecol. 2016, 92, 1–17. [Google Scholar] [CrossRef] [Green Version]
  39. Bell, T.H.; Callender, K.L.; Whyte, L.G.; Greer, C.W. Microbial Competition in Polar Soils: A Review of an Understudied but Potentially Important Control on Productivity. Biology 2013, 2, 533–554. [Google Scholar] [CrossRef]
  40. Arnalds, O.; Hallmark, C.T.; Wilding, L.P. Andisols from four different regions of Iceland. Soil Sci. Soc. Am. J. 1995, 59, 161–169. [Google Scholar] [CrossRef]
  41. Pushkareva, E.; Pessi, I.S.; Wilmotte, A.; Elster, J. Cyanobacterial community composition in Arctic soil crusts at different stages of development. FEMS Microbiol. Ecol. 2015, 91, fiv143. [Google Scholar] [CrossRef] [PubMed]
  42. Pushkareva, E.; Wilmotte, A.; Láska, K.; Elster, J. Comparison of Microphototrophic Communities Living in Different Soil Environments in the High Arctic. Front. Ecol. Evol. 2019, 7, 393. [Google Scholar] [CrossRef] [Green Version]
  43. Raabová, L.; Kovacik, L.; Elster, J.; Strunecký, O. Review of the genus Phormidesmis (Cyanobacteria) based on environmental, morphological, and molecular data with description of a new genus Leptodesmis. Phytotaxa 2019, 395, 1–16. [Google Scholar] [CrossRef]
  44. Tashyreva, D.; Elster, J. Production of Dormant Stages and Stress Resistance of Polar Cyanobacteria. In Life on Earth and other Planetary Bodies, Cellular Origin, Life in Extreme Habitats and Astrobiology; Hanslmeier, A., Kempe, S., Seckbach, J., Eds.; Springer Science + Business Media: Dordrecht, The Netherlands, 2012; Volume 24, pp. 367–386. ISBN 978-94-007-4965-8. [Google Scholar]
  45. Robinson, C.H. Cold adaptation in Arctic and Antarctic fungi. New Phytol. 2001, 151, 341–353. [Google Scholar] [CrossRef]
  46. Nilsson, R.H.; Tedersoo, L.; Abarenkov, K.; Ryberg, M.; Kristiansson, E.; Hartmann, M.; Schoch, C.L.; Nylander, J.A.A.; Bergsten, J.; Porter, T.M.; et al. Five simple guidelines for establishing basic authenticity and reliability of newly generated fungal ITS sequences. MycoKeys 2012, 4, 37–63. [Google Scholar] [CrossRef] [Green Version]
  47. Marañón-Jiménez, S.; Radujković, D.; Verbruggen, E.; Grau, O.; Cuntz, M.; Peñuelas, J.; Richter, A.; Schrumpf, M.; Rebmann, C. Shifts in the Abundances of Saprotrophic and Ectomycorrhizal Fungi With Altered Leaf Litter Inputs. Front. Plant Sci. 2021, 12, 1–17. [Google Scholar] [CrossRef] [PubMed]
  48. Borchhardt, N.; Baum, C.; Mikhailyuk, T.; Karsten, U. Biological soil crusts of arctic svalbard-water availability as potential controlling factor for microalgal biodiversity. Front. Microbiol. 2017, 8, 1–12. [Google Scholar] [CrossRef]
  49. Rippin, M.; Borchhardt, N.; Karsten, U.; Becker, B. Cold acclimation improves the desiccation stress resilience of polar strains of Klebsormidium (Streptophyta). Front. Microbiol. 2019, 10, 1730. [Google Scholar] [CrossRef] [Green Version]
  50. Roshan, S.K.; Dumack, K.; Bonkowski, M.; Leinweber, P.; Karsten, U.; Glaser, K. Taxonomic and functional diversity of heterotrophic protists (Cercozoa and endomyxa) from biological soil crusts. Microorganisms 2021, 9, 205. [Google Scholar] [CrossRef]
  51. Castillo-Monroy, A.P.; Benítez, Á.; Reyes-Bueno, F.; Donoso, D.A.; Cueva, A. Biocrust structure responds to soil variables along a tropical scrubland elevation gradient. J. Arid Environ. 2016, 124, 31–38. [Google Scholar] [CrossRef]
  52. Abed, R.M.M.; Tamm, A.; Hassenrück, C.; Al-Rawahi, A.N.; Rodríguez-Caballero, E.; Fiedler, S.; Maier, S.; Weber, B. Habitat-dependent composition of bacterial and fungal communities in biological soil crusts from Oman. Sci. Rep. 2019, 9, 1–16. [Google Scholar] [CrossRef]
  53. Elster, J. Ecological Classification of Terrestrial Algal Communities of Polar Environment; Springer: Berlin, Germany, 2002. [Google Scholar]
  54. Canfora, L.; Bacci, G.; Pinzari, F.; Lo Papa, G.; Dazzi, C.; Benedetti, A. Salinity and bacterial diversity: To what extent does the concentration of salt affect the bacterial community in a saline soil? PLoS ONE 2014, 9, e106662. [Google Scholar] [CrossRef] [Green Version]
  55. Rao, S.; Chan, Y.; Bugler-Lacap, D.C.; Bhatnagar, A.; Bhatnagar, M.; Pointing, S.B. Microbial Diversity in Soil, Sand Dune and Rock Substrates of the Thar Monsoon Desert, India. Indian J. Microbiol. 2016, 56, 35–45. [Google Scholar] [CrossRef] [Green Version]
  56. Pentecost, A. Order Volvocales. In The Freshwater Algal Flora of the British Isles. An Identification Guide to Freshwater and Terrestrial Algae; John, D., Whitton, B., Brook, A., Eds.; Cambridge University Press: Cambridge, MA, USA, 2002; pp. 303–327. [Google Scholar]
  57. Smirnova, G.V.; Oktyabrsky, O.N. Glutathione in bacteria. Biokhimiya 2005, 70, 1459–1473. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 09 02195 g001 550
Figure 1. General microbial profile of Icelandic biocrusts. The figures show the relative abundances based on the read numbers obtained by amplicon sequencing: (a) bacterial phyla; (b) cyanobacterial orders, (c) fungal phyla, and (d) eukaryotic clades. KRA—Krákunes, BOR—Borgarfjarðarbraut, FIF—Fiflholt, LSK—Litla-Skard, GIL—Giljar.
Figure 1. General microbial profile of Icelandic biocrusts. The figures show the relative abundances based on the read numbers obtained by amplicon sequencing: (a) bacterial phyla; (b) cyanobacterial orders, (c) fungal phyla, and (d) eukaryotic clades. KRA—Krákunes, BOR—Borgarfjarðarbraut, FIF—Fiflholt, LSK—Litla-Skard, GIL—Giljar.
Microorganisms 09 02195 g001
Microorganisms 09 02195 g002 550
Figure 2. Proportion of different functional groups of biocrust fungi identified by the FUNGuild database.
Figure 2. Proportion of different functional groups of biocrust fungi identified by the FUNGuild database.
Microorganisms 09 02195 g002
Microorganisms 09 02195 g003 550
Figure 3. OTU richness and Shannon index in Icelandic biocrusts. The sites are listed in order of elevation gradient from 11 to 157 m a.s.l., respectively. KRA—Krákunes, BOR—Borgarfjarðarbraut, FIF—Fiflholt, LSK—Litla-Skard, and GIL—Giljar. Same letters indicate no statistical difference between the samples according to one-way ANOVA followed by Tukey HSD post hoc test.
Figure 3. OTU richness and Shannon index in Icelandic biocrusts. The sites are listed in order of elevation gradient from 11 to 157 m a.s.l., respectively. KRA—Krákunes, BOR—Borgarfjarðarbraut, FIF—Fiflholt, LSK—Litla-Skard, and GIL—Giljar. Same letters indicate no statistical difference between the samples according to one-way ANOVA followed by Tukey HSD post hoc test.
Microorganisms 09 02195 g003
Microorganisms 09 02195 g004 550
Figure 4. Non-metric multidimensional scaling (NMDS) plots based on relative abundance of the OTUs in the biocrusts. Arrows indicate significant correlations (p < 0.05) with environmental variables. Phosphatase (acid), Pmic, C, and N refer to the potential activity of acid phosphatase, microbial P, total organic carbon, and total organic nitrogen, respectively.
Figure 4. Non-metric multidimensional scaling (NMDS) plots based on relative abundance of the OTUs in the biocrusts. Arrows indicate significant correlations (p < 0.05) with environmental variables. Phosphatase (acid), Pmic, C, and N refer to the potential activity of acid phosphatase, microbial P, total organic carbon, and total organic nitrogen, respectively.
Microorganisms 09 02195 g004
Table
Table 1. Cyanobacteria (a) and eukaryotic microalgae (b) isolated from the localities and their most closely related isolates retrieved from the NCBI database.
Table 1. Cyanobacteria (a) and eukaryotic microalgae (b) isolated from the localities and their most closely related isolates retrieved from the NCBI database.
(a)
CodeSiteNearest Cultivated Neighbor from NCBI Database (ID%)
KRA-1aKRAMK211225.1 Microcoleus vaginatus KZ-2-2-5 (94.8%)
KRA-4bMK211225.1 Microcoleus vaginatus KZ-2-2-5 (95.1%)
BOR-4cBORKU219728.1 Phormidesmis arctica HOR 11-2 (99.9%)
FIF-1aFIFKU219728.1 Phormidesmis arctica HOR 11-2 (99.8%)
FIF-1dMK211225.1 Microcoleus vaginatus KZ-2-2-5 (95.3%)
LSK-4cLSKMK211225.1 Microcoleus vaginatus KZ-2-2-5 (97.9%)
GIL-2aGILKJ939033.1 Phormidesmis sp. WJT36-NPBG15 (97.9%)
GIL-2bMN158649.1 Wilmottia murrayi ACKU582 (98.3%)
GIL-2cMK211225.1 Microcoleus vaginatus KZ-2-2-5 (98.1%)
GIL-3bKU219738.1 Phormidesmis arctica MUM 11-8 (97.5%)
(b)
CodeSiteNearest Cultivated Neighbor from NCBI Database (ID%)
KRA-2bKRAJQ315503.1 Chlamydomonas hedleyi KMMCC 188 (99.9%)
KRA-4aMK541759.1 Gonium pectorale CCAP 32/23 (99.3%)
BOR-1cBORAB936289.1 Chlorococcum lobatum (99.0%)
BOR-3bAY846370.1 Ankistrodesmus fusiformis Itas 8/18 M-7w (99.0%)
BOR-3cMN267184.1 Xerochlorella olmae UTEX B 2993 (99.3%)
BOR-4aMN248530.1 Chlorella vulgaris CCAP 211/21B (100%)
FIF-1bFIFMN267184.1 Xerochlorella olmae UTEX B2993 (100%)
FIF-1cMK541792.1 Chlorella vulgaris CCAP 211/19 (100%)
FIF-2aMH703776.1 Chlamydomonas callunae SAG 68.81 (99.8%)
FIF-3aMK262904.1 Klebsormidium elegans ACSSI 187 (100%)
FIF-3bMN248530.1 Chlorella vulgaris CCAP 211/21B (99.9%)
FIF-4aLT560349.1 Elliptochloris subsphaerica SAG 2202 (100%)
LSK-1aLSKMN267184.1 Xerochlorella olmae UTEX B 2993 (100%)
LSK-1bKM020100.1 Chlorococcum citriforme SAG 62.80 (100%)
LSK-1cKF673372.1 Chlorella emersonii SAG 2337 (99.9%)
LSK-4aMN248530.1 Chlorella vulgaris CCAP 211/21B (100%)
LSK-4bMT425956.1 Heterochlamydomonas sp. ACSSI 328 (99.9%)
GIL-3aGILKJ756841.1 Stichococcus bacillaris CCAP 379/5 (99.0%)
GIL-3dMG696565.1 Nannochloris sp. ACSSI 144 (100%)
GIL-4aEU878373.1 Parietochloris alveolaris UTEX 836 (99.3%)
Table
Table 2. Differences in microbial community composition among five sites revealed by PERMANOVA (9999 permutations). Asterisks corresponds to the level of significance (*** p < 0.001). Pairwise comparisons are shown in Supplementary Materials Table S2.
Table 2. Differences in microbial community composition among five sites revealed by PERMANOVA (9999 permutations). Asterisks corresponds to the level of significance (*** p < 0.001). Pairwise comparisons are shown in Supplementary Materials Table S2.
r2F
Bacteria0.453.07 ***
Cyanobacteria0.584.75 ***
Fungi0.412.59 ***
Eukaryotes0.543.47 ***
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pushkareva, E.; Barrantes, I.; Leinweber, P.; Karsten, U. Microbial Diversity in Subarctic Biocrusts from West Iceland following an Elevation Gradient. Microorganisms 2021, 9, 2195. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9112195

AMA Style

Pushkareva E, Barrantes I, Leinweber P, Karsten U. Microbial Diversity in Subarctic Biocrusts from West Iceland following an Elevation Gradient. Microorganisms. 2021; 9(11):2195. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9112195

Chicago/Turabian Style

Pushkareva, Ekaterina, Israel Barrantes, Peter Leinweber, and Ulf Karsten. 2021. "Microbial Diversity in Subarctic Biocrusts from West Iceland following an Elevation Gradient" Microorganisms 9, no. 11: 2195. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9112195

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