Next Article in Journal
Taxonomic, Functional, and Phylogenetic Diversity of Bats in Urban and Suburban Environments in Southern México
Next Article in Special Issue
Coahuilasaurus lipani, a New Kritosaurin Hadrosaurid from the Upper Campanian Cerro Del Pueblo Formation, Northern Mexico
Previous Article in Journal
Environmental DNA as Early Warning for Alien Species in Mediterranean Coastal Lagoons: Implications for Conservation and Management
Previous Article in Special Issue
Molecular Phylogeny and Historical Biogeography of Byrsonima (Malpighiaceae) Corroborates the Mid-Miocene Origins of Neotropical Savannas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flourishing in Darkness: Protist Communities of Water Sites in Shulgan-Tash Cave (Southern Urals, Russia)

by
Natalia E. Gogoleva
1,2,3,*,
Marina A. Nasyrova
2,
Alexander S. Balkin
2,3,
Olga Ya. Chervyatsova
4,
Lyudmila Yu. Kuzmina
5,
Elena I. Shagimardanova
3,6,7,
Yuri V. Gogolev
3,8 and
Andrey O. Plotnikov
9,10,*
1
Research Department for Limnology, Mondsee, Universität Innsbruck, 5310 Mondsee, Austria
2
Laboratory of Biomedical Technologies, Institute for Cellular and Intracellular Symbiosis, Ural Branch of Russian Academy of Sciences, 460000 Orenburg, Russia
3
Institute of Fundamental Medicine and Biology, Kazan Federal University, 420111 Kazan, Russia
4
State Nature Reserve “Shulgan-Tash”, 453585 Irgyzly, Russia
5
Ufa Institute of Biology, Ufa Federal Research Center, Russian Academy of Sciences, 450054 Ufa, Russia
6
SBHI Moscow Clinical Scientific Center Named after Loginov MHD, 111123 Moscow, Russia
7
Life Improvement by Future Technologies Center, Skolkovo, 143025 Moscow, Russia
8
Kazan Institute of Biochemistry and Biophysics, Federal Research Center “Kazan Scientific Center of the Russian Academy of Sciences”, 420111 Kazan, Russia
9
Faculty of Biology, Shenzhen MSU-BIT University, Shenzhen 518172, China
10
“Persistence of Microorganisms” Science Resource Center, Institute for Cellular and Intracellular Symbiosis, Ural Branch of Russian Academy of Sciences, 460000 Orenburg, Russia
*
Authors to whom correspondence should be addressed.
Diversity 2024, 16(9), 526; https://doi.org/10.3390/d16090526
Submission received: 7 July 2024 / Revised: 22 August 2024 / Accepted: 26 August 2024 / Published: 1 September 2024
(This article belongs to the Special Issue Diversity in 2024)

Abstract

:
Karst caves, formed by the erosion of soluble carbonate rocks, provide unique ecosystems characterized by stable temperatures and high humidity. These conditions support diverse microbial communities, including wall microbial fouling, aquatic biofilms, and planktonic communities. This study discloses the taxonomic diversity of protists in aquatic biotopes of Shulgan-Tash Cave, a culturally significant site and popular tourist destination, by 18S rRNA gene metabarcoding. Our findings reveal the rich protist communities in the cave’s aquatic biotopes, with the highest diversity observed in Blue Lake at the cave entrance. In contrast, Distant Lake in the depth of the cave was inhabited by specific communities of plankton, mats, and pool fingers, which exhibited lower richness and evenness, and were adapted to extreme conditions (cold, darkness, and limited nutrients). High-rank taxa including Opisthokonta, Stramenopiles, and Rhizaria dominated all biotopes, aligning with observations from other subterranean environments. Specific communities of biotopes inside the cave featured distinct dominant taxa: amoeboid stramenopile (Synchromophyceae) and flagellates (Choanoflagellatea and Sandona) in mats; flagellates (Choanoflagellatea, Bicoecaceae, Ancyromonadida) and amoeboid protists (Filasterea) in pool fingers; flagellates (Ochromonadales, Glissomonadida, Synchromophyceae), fungi-like protists (Peronosporomycetes), and fungi (Ustilaginomycotina) in plankton. The specificity of the communities was supported by LEfSe analysis, which revealed enriched or differentially abundant protist taxa in each type of biotope. The predominance of Choanoflagellatea in the communities of cave mats and pool fingers, as well as the predominance of Synchromophyceae in the cave mats, appears to be a unique feature of Shulgan-Tash Cave. The cold-tolerant yeast Malassezia recorded in other caves was present in both plankton and biofilm communities, suggesting its resilience to low temperatures. However, no potentially harmful fungi were detected, positioning this research as a baseline for future monitoring. Our results emphasize the need for ongoing surveillance and conservation efforts to protect the fragile ecosystems of Shulgan-Tash Cave from human-induced disturbances and microbial invasions.

1. Introduction

Karst caves are often formed by the erosion of soluble carbonate rocks. On the one hand, the stable temperature and high humidity make karst caves attractive to various organisms. On the other hand, the cold environment and lack of nutrients present challenges that require special adaptations. Wall microbial fouling, aquatic biofilms, and planktonic communities are commonly found in karst caves [1,2,3,4]. The invasion of underground cavities can occur by microorganisms living on the surface and in the epikarst zone [5]. However, most cave-dwelling species are not found on the surface, suggesting that underground ecosystems form cave-specific communities [6]. To date, several prokaryotic cave communities that support subterranean ecosystems and participate in the primary production of organic matter have been described in detail [4,7,8,9]. As for protists, which are a crucial part of the food chains [10,11], data on their diversity in karst caves are very limited. The first descriptions of protist diversity in cave communities have been presented by Gittleson and Hoover [12,13] and Hill et al. [14], using culture methods followed by microscopy. In recent years, in addition to classical methods [15,16,17], the use of DNA metabarcoding and metagenomics has become increasingly popular [18,19,20] for studying environmental communities. Although the taxonomic resolution of protists remains challenging in metabarcoding, the V4 region of the 18S rRNA gene is considered a rather reliable taxonomic marker for the discrimination of many protist taxa, especially high-rank ones [21,22].
Historical artifacts and the exhilarating atmosphere of the caves make them the subject of tourist attention. However, tourism-related anthropization leads to changes in microclimate and intense growth of allochthonous microorganisms forming spots and biofilms, e.g., in Petralona Cave (Chalkidiki, Greece) and Lascaux Cave (France) [20,23]. Moreover, the equipment of tourist routes with artificial lighting encourages the development of lampenflora, which is uncharacteristic for caves [20]. The term “lampenflora” (also known as lamp-flora or lamp flora) was introduced by a botanist Klaus Dobat in the 1960s [24] and refers to all phototrophic organisms (cyanobacteria, algae, plants) growing in permanently illuminated areas. Frequent visits to caves cause an increase in air temperature and the introduction and development of extraneous microorganisms, including potential pathogens. The development of bacteria and fungi in Paleolithic caves can result in stained spots formation on the cave walls and damage Paleolithic art due to microbial alterations [23,25]. The conservation of these fragile, unique cave ecosystems requires their close study and monitoring.
Shulgan-Tash Cave, located in the Southern Urals (the Republic of Bashkortostan, Russia), is a well-known cultural heritage site and a popular tourist destination [26,27]. Studying and monitoring the microbial communities of Shulgan-Tash Cave is an important task not only to characterize the underground ecosystems but also to assess human-induced threats to them. We have previously described the biodiversity of prokaryotic communities of wall fouling in Shulgan-Tash Cave [4,28,29]. In this work, we applied the method of metabarcoding based on the 18S rRNA gene to describe the general taxonomic diversity of protists in aquatic biotopes of the cave. We also performed a comparative analysis of the described subterranean communities with communities of Blue Lake, a similar water site on the surface, in the immediate vicinity of the cave.

2. Materials and Methods

2.1. Cave Site Description and Sampling

Shulgan-Tash Cave (also known as Kapova Cave) is located within the territory of Shulgan-Tash Nature Reserve (53°02′ N, 57°03′ E) in the Tirmentau massif (the height of the massif is 420 m.a.s.l., and the entrance to the cave is 280 m.a.s.l.) (Figure 1).
Distant Lake (Russian name “Bepxнee Дaльнee”, spelled “Verhnee Dalnee”) is a low-flow groundwater body located 700 m from the cave entrance, at an elevation of 305 m.a.s.l. and lying 90 m below the earth’s surface (Figure 1A). The lake extends 29 m from south to north, is 11 m wide, and has a maximum depth of 1.65 m [30]. The lake is recharged by karst water from a vadose zone. A permanent inflow comes from a dead-end fissure passage located 10 m from the northern shore (Figure 1B). The magnitude of the inflow ranges from 50 to 200 mL/sec during dry periods and from 1 to 2 L/sec during the snowmelt season. A small creek flows out of the lake at the southern shore and is absorbed by loose sediments after 20 m (Figure 1B). The vadose water of Distant Lake, when infiltrated by the outflowing stream, recharges the phreatic water of the Blue Lake (Russian name “Гoлyбoe”, spelled “Goluboe”) hydrosystem located at the entrance of the cave [31]. Vaucluse-type karst spring Blue Lake is characterized by a flow rate ranging from 230–240 L/sec during low-water periods to more than 3000 L/sec during rain floods and snowmelt. The dimensions of this lake-spring are 37 by 10 m, and the maximum depth of the phreatic zone channels reached by divers exceeds 80 m.
Figure 1. Map of the research area and Shulgan-Tash Cave sampling sites: (A) Overview map showing the global position of the research area. (B) The position of the cave in relation to the Southern Urals. (C) Digital model of the cave combined with the relief of the Tirmentau massif. Visualization was created using Cloud Compare 2020 software with laser scanning point clouds [32]. The cavities of the phreatic zone are shown schematically (by Snetkov E., unpublished data). (D) Sampling scheme of Distant Lake and its watercourses. PF—pool finger biofilm; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Figure 1. Map of the research area and Shulgan-Tash Cave sampling sites: (A) Overview map showing the global position of the research area. (B) The position of the cave in relation to the Southern Urals. (C) Digital model of the cave combined with the relief of the Tirmentau massif. Visualization was created using Cloud Compare 2020 software with laser scanning point clouds [32]. The cavities of the phreatic zone are shown schematically (by Snetkov E., unpublished data). (D) Sampling scheme of Distant Lake and its watercourses. PF—pool finger biofilm; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Diversity 16 00526 g001
The waters of Distant Lake and its watercourses have HCO3-Ca major ions, a stable TDS (Total dissolved solids) of 4000–4450 ppm, a pH close to neutral (6.5–7.5), and an almost constant temperature of 6.2–6.4 °C. A notable feature of the atmosphere over the lake is the elevated CO2 concentration, ranging from 0.8–1.2% near the northern shore to 2.2% near the southern shore throughout the year. The waters of the Blue Lake karst spring are characterized by HCO3-Ca composition, TDS 0.26–3.16 g/L, and a slightly alkaline reaction (pH 6.5–7.9). The water temperature varies from 2 °C to 9 °C (unpublished data).
Active biogenic-mineral subaquatic stalactoids, known in the literature as pool fingers [33], are found 5–10 cm below the water surface along the west and east shores. Additionally, cascades of calcite gyre dams are formed along the outflowing creek, within which filamentous littoral microbial mats develop (Figure 2).
The samples for metabarcoding were collected in July 2018. Sampling sites are marked in Figure 1 and included the following locations: PF—biofilm on the terminal part of pool fingers (stalactoids immersed in water); MDL—littoral microbial mats of Distant Lake; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; BL—water from Blue Lake (Figure 2). MDL biofilm samples were collected using a mechanical single-channel pipette (100–1000 µL) with sterile filtered pipet tips and placed in sterile tubes (1.5 mL) with 1 mL RNAlater solution (Invitrogen, Carlsbad, CA, USA) (Figure 2b). PF samples were collected by plucking the tip of subaquatic stalactoids. Biofilm samples were collected in quadruplicate. All water samples were collected from a depth of 10 cm. Water samples of 1.5–3 L were collected in sterile bottles, brought to the surface, and filtered through a 0.22 µm membrane immediately after collection. The membranes were cut into four pieces equal in size, which were then placed in RNAlater solution. Two pieces of each filter were used for downstream analysis. All 16 samples (4 PF, 4 MDL, 2 CEW, 2 WPF, 2 WDL, 2 BL) were transported to the laboratory and stored at −24 °C until DNA extraction.

2.2. DNA Extraction, PCR Amplification, and Sequencing

DNA was extracted using the DNeasy PowerBiofilm Kit (Qiagen, Germany) with a FastPrep-24 homogenizer (MP Biomedicals, Solon, OH, USA), according to the protocol provided by the manufacturer. The concentration of DNA was assessed using Qubit (Invitrogen, Carlsbad, CA, USA) and the Qubit dsDNA High Sensitivity Kit (Invitrogen, Carlsbad, CA, USA). The V4 region of the 18S rRNA gene was amplified with the universal eukaryotic UNonMet primers EK-565F (5′-GCAGTTAAAAAGCTCGTAGT-3′) [34] and EUK-1134R-UNonMet (5′-TTTAAGTTTCAGCCTTGCG-3′), which are biased against Metazoa [35] with Illumina i5 and i7 adapters (https://support-docs.illumina.com/SHARE/AdapterSequences/Content/SHARE/AdapterSeq/Nextera/SequencesNextera_Illumina.htm (accessed on 15 may 2024). A negative sterile water control was always performed to confirm the absence of DNA contamination. The libraries were sequenced using a MiSeq v3 Reagent Kit 600 cycles (Illumina, San Diego, CA, USA) on the MiSeq platform (Illumina, San Diego, CA, USA). All DNA manipulations and sequencing were performed at the Joint KFU-Riken Laboratory, Kazan Federal University (Kazan, Russia).

2.3. Bioinformatic Treatment and Data Analysis

Bioinformatic treatment and data analysis were performed at the “Persistence of Microorganisms” Science Resource Center of the Institute for Cellular and Intracellular Symbiosis UrB RAS (Orenburg, Russia).
The quality of sequenced reads was checked using FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/, accessed on 1 December 2022). Primer sequences were removed using Cutadapt v. 4.1 [36]. DADA2 v. 1.28 [37] was used for further sequence analysis. Reads were sequentially filtered by quality with a maximum expected error threshold of 4.0 (maxEE = 4) and a minimum length of 250 bp, then merged with a minimum length of 12 bp overlap. Amplicon sequence variants (ASVs) were generated with the pseudo-pooling option, and chimera sequences were removed. ASVs were additionally filtered for a minimum length of 400 bp and a maximum length of 520 bp. The obtained ASVs were taxonomically classified using a pre-trained naive Bayes classifier, which was trained on the PR2 database v. 5.0.0 [38]. Taxa with low support were additionally manually checked against the NCBI nucleotide database using the BLAST algorithm. Non-eukaryotic ASVs and those belonging to Embryophyta, chloroplasts, mitochondria, and Metazoa were removed before further analysis.
All statistical analyses and visualizations were performed within the R environment v. 4.4.0 using phyloseq v. 1.48.0 [39], ggplot2 v. 3.5.1 [40], MicrobiotaProcess v. 1.16.0 [41], and microbiome v.1.26.0 [42] packages. Indices of alpha diversity were calculated using the alpha function from the microbiome R package [42]. The Bray–Curtis dissimilarities and weighted UniFrac distances [43] were used for β-diversity comparison of protist communities based on the ASV distribution and phylogenetic differences, respectively. Permutational multivariate analyses of variance (PERMANOVA) of Bray–Curtis dissimilarities and weighted UniFrac distances were conducted using the “mp_adonis” function in the MicrobiotaProcess [41] package with 999 permutations. Principal coordinate analysis (PCoA) plots of Bray–Curtis dissimilarity or weighted UniFrac distances (based on Hellinger-transformed ASV relative abundances) were constructed using the “mp_plot_ord” function in the MicrobiotaProcess [41] package.
Marker taxa were identified using LEfSe (Linear discriminant analysis Effect Size) analysis [44] on CPM (counts per million)-normalized ASV abundance, with multiple testing correction options (“multigrp_strat”), as implemented in the microbiomeMarker v. 1.6.0 [45] package. This process was set with the Kruskal (p = 0.05) test based on linear discriminant analysis (LDA) effect size (LEfSe) and Wilcox (p = 0.05), corrected with False Discovery Rate (FDR).

2.4. Functional Annotation of ASVs

In this study, we used previously published works [46,47] to assign ASVs into three major functional groups: consumers (Ancyromonadida, Ciliophora, Rhizaria, Obazoa non-Ichthyosporea, CRuMs (Collodictyonida, Rigifilida, and Mantamonadida), Amoebozoa, non-Xanthophyceae and non-Peronosporomycetes Stramenopiles, Centroplasthelida, Telonemia, Malawimonadida, and Diplonemidae), phototrophs (Archaeplastida, Cryptophyceae, and Xanthophyceae), and parasites (Apicomplexa, Peronosporomycetes, Phytomyxea, Perkinsida, and Rozellomycota). Besides, we considered as parasites exclusively predatory protists such as Orciraptor, Viridiraptor, Allantion, and Colpodellida. All other groups of protists, for which the function could not be clearly defined, were classified as “unknown”.

3. Results

In this study, we conducted a comparative evaluation of the diversity, taxonomic composition, and functional composition of protist communities in Distant Lake inside Shulgan-Tash Cave in contrast to Blue Lake located outside the cave. Two types of biotopes, plankton and biofilm (mats and pool fingers), were investigated in Distant Lake. The rarefaction curves formed based on the Chao1 index calculated with ASVs reached a plateau except for BL communities, indicating that the sequencing depth of each sample inside the cave was sufficient (Figure S1).

3.1. Diversity and Taxonomic Composition of Protist Communities

3.1.1. Alpha Diversity

Diversity estimation based on sequencing the V4 region of the 18S rRNA gene revealed fairly rich and diverse communities of protists in water sites of Shulgan-Tash Cave: 1621 ASVs in total, with 450–474 ASVs in the richest community of Blue Lake and 41–60 ASVs in the least rich community of pool fingers (Table S1). The taxonomic richness of communities, assessed using the Chao1 index, in Blue Lake, located in the illuminated site of the cave entrance, was about twice as high as in the inner cave sites (Wilcoxon rank-sum test, p ≤ 0.05; Figure 3). Similar ratios were observed using the Shannon index (diversity assessed with a focus on rare taxa, which depends more on richness than evenness). The Gini–Simpson index (diversity assessed with a focus on common taxa, which depends more on evenness than richness) and Pielou’s evenness (evenness of relative abundances for all ASVs) demonstrated similar patterns of alpha diversity in the protist communities, which were substantially lower even in relatively richer biofilm communities (MDL and PF of Distant Lake) than in more even water communities (Figure 3). A pairwise comparison of the alpha diversity indices for protist communities of the cave biotopes did not reveal significant differences, except for the Chao1 richness (Wilcoxon rank-sum test, p ≤ 0.05) between the communities of MDL and PF sites. Nevertheless, the Chao1 richness values of the communities in CEW, MDL, and WPF sites were higher than those in PF and WDL sites.

3.1.2. Beta Diversity

Principal coordinate analysis based on the Bray–Curtis and weighted UniFrac distances grouped the cave protist communities in three clusters according to the main cave biotopes: inner and outer water sites, cave aquatic biofilms, and pool fingers (Figure 4). Protist communities were confirmed to be associated with specific biotopes by both Bray–Curtis dissimilarity (PERMANOVA; F = 5.03, R2 = 0.71, p ≤ 0.001) and the weighted Unifrac algorithm (PERMANOVA; F = 9.46, R2 = 0.82, p ≤ 0.001). In general, all protist communities inside Shulgan-Tash Cave differed significantly from the communities of Blue Lake (weighted UniFrac; F = 4.93, R2 = 0.26, p ≤ 0.01). Communities of different biotopes in Distant Lake such as plankton (CEW, WDL, WPF), mats (MDL), and pool fingers (PF) also differed significantly from each other (weighted UniFrac; F = 8.97, R2 = 0.62, p ≤ 0.001).
Blue Lake protist community featured ASV11 and ASV22, which belonged to an unclassified phototrophic chrysophyte alga Ochromonadales. The cave mat community was characterized by specific ASV2 (unclassified amoeboid photosynthetic or heterotrophic stramenopile of the class Synchromophyceae), ASV7, and ASV18 (both belonged to a heterotrophic cercomonad flagellate Sandona). The latter two ASVs have homologous sequences in GenBank assigned to Sandona hexamutans HQ918174.1 [48] and Sandona pentamutans EU709147.1 [49] with 99% and 100% identity, respectively. ASV3 was a feature of the pool finger protist community and was identified as an unclassified heterotrophic flagellate from the family Bicoecaceae, whereas ASV4 characterized the communities of both pool fingers and mats and belonged to an unclassified heterotrophic flagellate Choanoflagellatea.
Although the protist communities of the water samples from Distant Lake and Blue Lake formed isolated clusters based on the Bray–Curtis dissimilarity, PCoA ordination using phylogenetically weighted UniFrac distances placed all water samples in a common cluster (Figure 4), which confirmed, at the supraspecies level, the close phylogenetic relatedness of the protists living in the water inside and outside the cave.

3.1.3. Taxonomic Composition

A total of 9 supergroups of eukaryotes, 16 divisions, 23 subdivisions, 61 classes, 133 orders, 186 families, and 272 genera were identified, except for unidentified taxa (Figures S2–S8). However, the following discussion focuses only on protist taxa (ASVs) represented by more than 1% of the reads, namely 5 supergroups, 8 divisions, 12 subdivisions, 24 classes, 31 orders, 33 families, and 32 genera (Table S2). The taxonomic structure of the protist communities at the genus level and at higher taxonomic levels demonstrated similar patterns. In all biotopes, the most abundant supergroups were Obazoa, represented primarily by Opisthokonta, and TSAR, mainly represented by Stramenopiles, Alveolata, and Rhizaria (Figures S2 and S3; Table S2). Blue Lake (BL) had a significantly higher proportion of Archaeplastida compared to the other sites (PF, MDL, CEW, WPF, WDL) (Figure S2; Table S2).
Supergroups Haptista, Cryptista, and CRuMs were present in small numbers, each representing less than 1% of the total abundance in all communities (Figure S2; Table S2).
Stramenopiles dominated in all biotopes studied, occupying 26–49% of the total abundance in the cave lake communities and up to 78–81% in Blue Lake (Figure S3; Table S2). Blue Lake was dominated by the mixotrophic flagellates unclassified Ochromonadales (28–32%) and heterotrophic flagellates Spumella (22–23%), whereas in pool fingers, unclassified Bicoecales were the most abundant (24–37%) (Figure 5 and Figures S6 and S8; Table S2). In MDL, Stramenopiles were represented by unclassified Bicoecales, Synchromophyceae, and unclassified Gyrista (Figure 5 and Figures S5 and S6). Diatoms, mainly Melosirales, were found in different samples (Figure 5 and Figure S6; Table S2).
Although Opisthokonta were found in all samples, their proportion was the largest in the communities of the oligotrophic Distant Lake, ranging from 10 to 39% (Figure S3). Moreover, the biofilm communities of mats and pool fingers were characterized by the predominance of sessile heterotrophic flagellates, which belonged to unclassified Choanoflagellatea (7–18%) (Figure 5 and Figures S5 and S6; Table S2).
Fungi were represented by Basidiomycota (Ustilaginomycotina, Pucciniomycotina, and Agaricomycotina), Chytridiomycota, and Ascomycota (Saccharomycotina), occupying a significant proportion in all aquatic communities (11–15% in WDL; 15–20% in CEW; 32–35% in WPF). Surprisingly, fungi-like saprotrophic and parasitic protists Peronosporomycetes (4–14%) and fungi Ustilaginomycotina (6–23% in CEW and WPF only) were the abundant orders in aquatic communities (Figures S5 and S6; Table S2).
The pool fingers community was the only one in which representatives of Filasterea (heterotrophic amoeboid protists) were present as the dominant taxon (Figure S4). Interestingly, Ancyromonadida also constituted a large proportion (4–13%) in the pool finger community (Figure S3; Table S2).
Rhizaria, another representative of TSAR, was predominant in the inner cave communities and differed from the community of Blue Lake in both diversity and genera representation. The presence of Sandona (up to 45%) was a feature of the mat community. Rhizaria were also represented by the genera Neoheteromita, Paracercomonas, and Teretomonas (Figure 5 and Figure S3).

3.2. Protist Taxa with Differential Abundance between Main Groups of Samples

LEfSe analysis was used to identify biomarkers (enriched or differentially abundant taxa) associated with the cave biotopes and an outside Blue Lake (Table S3). In Blue Lake communities, 71 genera had significantly larger proportions in contrast to the cave sites (p < 0.01, FDR, LDA > 2). Obviously, the most represented genera belonged to the photosynthetic members of classes Chrysophyceae (Hydrurus, heterotrophic flagellates Paraphysomonas), Xanthophyceae (Vaucheria), Chlorophyceae (Planophila), and Bacillariophyceae (Planothidium, Cocconeis, Diatoma) (Table S3).
Among cave communities, as expected, planktonic and biofilm (mat and pool finger) communities differed in numerous markers. In the aquatic samples CEW, WDL, and WPF, 70 differentially abundant taxa were identified, among which Alveolata, Cryptophyta, and Discoba were significantly enriched (p < 0.05, FDR, LDA > 4) (Table S3). Ancyromonadida and Evosia were enriched in the mat and pool finger communities. The enriched genera in the cave planktonic communities included Spumella, Segregatospumella, Oikomonas, an unidentified genus of Saprolegniales, Melosira, Wickerhamomyces, Rhodotorula, Colpidium, Malassezia (Figure S9; Table S3). At the same time, taxa enriched in the MDL and PF together included mainly taxa of suprageneric levels, namely divisions Ancyromonadida and Stramenopiles, subdivision Filasterea, classes Choanoflagellatea and Variosea, and family Bicoecaceae (Figure S9; Table S3).
In the communities of the inner cave biotopes (CEW, WDL, WPF, PF, and MDL), 88 taxa with significant differential abundance (p < 0.05, FDR, LDA > 2) were found (Figure S10; Table S3). Generally, in every cave biotope, the specifically enriched division was distinguished: Rhizaria (LDA = 5.1) in MDL, Ancyromonadida (LDA = 4.9) in PF, and Alveolata (LDA = 5.1) in plankton, respectively. At the genus level, the communities of mats were enriched with genera Sandona, Flectomonas, Cercomonas (belonging to the Rhizaria clade), Bracteacoccus, and Archaeorhizomyces. Pool finger communities were enriched with only one taxon of the genus level, Schizoplasmodium; 21 enriched taxa of higher taxonomic levels were found. Segregatospumella, Spumella, Paracercomonas, Cryptococcus, and Pseudocyrtolophosis genera were enriched in the communities of the cave plankton.

3.3. Functional Composition of Protist Communities

We investigated the functional diversity of protists in Blue Lake and inner cave sites, based on the relative abundance of consumers, phototrophic organisms, and parasites (Figure 6 and Figure S11). Out of all, 1118 (69%) ASVs were classified as consumers, 250 (15.5%) as phototrophs, 139 (8.5%) as parasites (including some predatory protists), and 114 ASVs (7%) were assigned to unknown categories due to their unclear functional roles. Consumers were prevalent (84.8%) in the Distant Lake communities, whereas phototrophic and parasitic groups had minor proportions (5.3% and 4.5%, respectively) (Figure 6). Phototrophic protists were prevalent in the plankton community of Blue Lake (55%), although consumers were also widely represented (37.2%). Parasites were equally abundant in communities inside the cave (4.4%) and in Blue Lake (4.5%) and were mostly represented by parasitic Peronosporomycetes (Gyrista) and Rozellomycota (Fungi). Consumers were represented mainly by the supergroup TSAR and its division Stramenopiles. Phototrophic protists were mostly from the division Stramenopiles (TSAR) and from the division Chlorophyta (supergroup Archaeplastida) with a smaller proportion.
The functional diversity of protists in Distant Lake was analyzed in three biotopes, including plankton (CEW, DLL, WPF), pool fingers (PF), and mats (MDL). Consumers were dominant in all biotopes, accounting for 94.5%, 89%, and 72% in PF, MDL, and cave plankton, respectively (Figure S11). Parasitic (11.9%) and phototrophic (11.6%) protists comprised a fairly significant part of cave plankton communities, while in PF and MDL, they constituted only small proportions (0.08% and 0.58% for parasites; 2.2% and 1.1% for phototrophs, respectively) (Figure S11).

4. Discussion

DNA metabarcoding has revealed quite rich protist communities in the water sites of Shulgan-Tash Cave. Among the studied cave biotopes, the richest community was found in Blue Lake located at the cave entrance, whereas the least even communities were in the cave mats (MDL) and pool fingers (PF) (Figure 3). Distant Lake is located in a remote dead-end gallery of the cave, where, apart from the absence of photosynthesis, there is almost no direct influx of air and water from the surface, except for scant infiltration through the thickness of the epikarst and condensate from the upper gallery. Therefore, the higher alpha diversity metrics (richness and Shannon index) for protist communities of Blue Lake compared to the cave sites can be easily explained by the intensive growth of photosynthetic protists (Figure 6), many of which are absent in Distant Lake, and the intake of organic matter and nutrients, as well as other protists, into Blue Lake from surrounding terrestrial habitats, as noted generally for subsurface aquatic habitats [50]. At the same time, water from Distant Lake most likely flows into the underground Shulgan River and then into Blue Lake (Figure 1). This connection suggests the ingress of microorganisms from Distant Lake into Blue Lake, but a comparison of taxonomic lists for both locations revealed only 51 common ASVs, which may be attributed to both transfer through water flow and a poor discriminant power of the V4 region of the 18S rRNA for some closely related taxa of protists (Figure S12). This finding clearly indicates that in the depth of the cave, the specific communities of Distant Lake adapted to extreme conditions (cold, darkness, and substantially limited nutrients) are formed, as noted for other caves [6,19,51]. According to our preliminary data (unpublished), a unique microbial community represented by mats has formed in Distant Lake (Figure 2), with some members possibly acting as primary chemosynthetic producers.
The protist communities in pool fingers (PF) and mats (MDL) were characterized by the least evenness, with a high proportion of a few dominant taxa. Dominant taxa in MDL communities included Sandona, unclassified Synchromophyceae, and Choanoflagellatea, whereas, in PF, they were represented mostly by unclassified Choanoflagellatea and Filasterea, families Bicoecaceae and Sandonidae. All plankton communities inside the cave and outside (Blue Lake) were more even than the biofilm communities because they were richer, and their dominant taxa occupied lower proportions in the communities. In all sites and biotopes, three high-rank taxa of protists, Stramenopiles, Opisthokonta, and Rhizaria, dominated (Figure 4 and Figure S3). Stramenopiles is a rather diverse taxon, including heterotrophic flagellates, chrysophyte and diatom algae, fungi-like protists such as oomycetes and labyrinthulomycetes, etc. Stramenopiles have been revealed to prosper in other caves and groundwaters [3,20,50,52]. Opisthokonta, another megataxon of eukaryotes including heterotrophic organisms such as choanoflagellates, fungi, and multicellular animals (Metazoa), is also known to be dominant in caves and groundwaters [20,23,50]. However, in all mentioned ecosystems, the dominant Opisthokonta were represented mainly by different fungi, such as Ascomycota, Basidiomycota, Mucoromycota, and Chytridiomycota, but not Choanoflagellatea [20,23,50]. The predominance of Choanoflagellatea in communities of the lake mats and pool fingers seems to be a unique feature of Shulgan-Tash Cave. Choanoflagellates have been mentioned as dominants in the cave biotopes only once, namely in the pristine cave Allas in France [23]. However, a single finding of choanoflagellate Salpingoeca amphoridium is known from a plankton sample of the Echo River, a subterranean stream in Mammoth Cave (Kentucky, USA) [53].
Dominant Rhizaria were represented by bacteriotrophic (Sandona, Flectomonas, etc.) and predatory (Allantion, Orciraptor) cercozoan flagellates, which agrees with other observations of Cercozoa predominance in caves and subsurface waters [3,20,23,50,52]. Representatives of the genus Sandona (subdivision Cercozoa) have previously been described as small, rounded, gliding cells that are found ubiquitously in soils [48,49], likely under conditions similar to those of biofilms attached to underground stream beds. Interestingly, Cercozoa dominated ice communities, reaching 72% of total sequences, which confirms their ability to grow under low temperatures in caves [54,55]. Cercozoa also occupied a significant proportion of sediments and biofilms in the hot springs [56,57].
Pool fingers in Shulgan-Tash Cave, unique stalactoid speleothems formed in the cave water, were inhabited by specific communities [30]. The biofilms on the terminal parts of pool fingers were dominated mainly by sessile and gliding heterotrophic flagellates, including Choanoflagellatea, Ancyromonadida, cercomonad Sandona, and heterotrophic amoeboid protists Filasterea. Currently, it is the first study of the protist community of pool fingers, and further investigations are necessary to fully characterize the unique and enigmatic communities of pool fingers and the role of protists in their formation and functioning. Another interesting finding, noted for the first time for caves, seems to be a specific feature of Shulgan-Tash Cave: the predominance of Synchromophyceae in the cave mats. Synchromophyceae, amoeboid stramenopiles, include only three genera inhabiting surface waters: heterotrophic Leukarachnion and algae Chlamydomyxa labyrinthuloides and Synchroma grande [58].
The presence of diatom alga Melosira and green alga Bracteacoccus is not a specific trait of Shulgan-Tash Cave. Previously, algae in the cave were recorded using microscopic methods [59]. Algoflora has been found in almost all available habitats, including caves [15,51]. A number of studies show the possibility of diatom mixotrophy. Particularly, Claus [60] and Abdullin [61] performed growth experiments with algae in a cave and found that some species are able to survive in mineral-rich environments, even in complete darkness. Alternatively, some algae may be found in the caves due to their anthropogenic invasion [23]. At the moment, we have limited suggestions about the algoflora in caves, but as information is accumulating and the most predominant taxa are being identified, it would be promising to search for algae, especially diatoms, in the underground world.
The parasite representation was quite similar to that of other freshwater bodies, where the percentage of parasites was approximately 5% [46,47]. Interestingly, the taxonomic composition of the parasites differed between the communities of Blue Lake and the communities inside the cave. Namely, members of the Rosellomycota class, which infect fungi and green algae [62], dominated among the parasites of Blue Lake. The parasites of Distant Lake communities were mainly presented by predatory members of Glissomonadida and fungi-like protists Peronosporomycetes. The high relative abundances of parasitic fungi Ustilaginomycotina and fungi-like protists Peronosporomycetes seem to be a specific trait of water communities inside Shulgan-Tash Cave, as previous observations have not found these protists in large proportions in other caves [3,20,23,63,64,65,66], although Peronosporomycetes were predominant in groundwaters [50].
Fungi Malassezia and Mortierella (Tables S1 and S2) found in Shulgan-Tash Cave have been recorded previously in other caves [64,65,66]. Cold-resistant Malassezia is associated with animals and is capable of growing at low temperatures [65,67]. It has been shown that fungal outbreaks in caves may be associated with the development of Mucor circinelloides and Cladosporium microsporum in Castañar Cave, Spain [66], and Fusarium solani in Lascaux Cave, France. [68]. These fungal outbreaks in caves clearly demonstrate the fragility of cave ecosystems. In this regard, caves as specific ecosystems require a more reverent attitude. No potentially dangerous fungi were found in Shulgan-Tash Cave in our study, indicating that it can serve as a reference point for future monitoring of the cave.
In conclusion, it is worth emphasizing that the specific composition of the protist communities is obviously determined by different environmental conditions in the water sites inside and outside Shulgan-Tash Cave. Unfortunately, environmental physicochemical parameters were not measured during sampling, but some features of the biotopes studied might be identified easily. The main difference of Distant Lake from Blue Lake is full darkness, constant cold, and substantial limitation of nutrients. However, a high concentration of CO2 may contribute to the growth of autotrophic prokaryotes capable of chemosynthesis, especially in mats (unpublished data). These prokaryotes, together with others brought by karst water from a vadose zone, are the first component of the simplest food chain in the cave community and serve as food for abundant bacterivorous protists, mainly flagellates and, more rarely, amoebae. Bacteriotrophic protists with different types of swimming behavior and movement were revealed in the cave protist communities, as shown in a microscopic study of protists in aquifers [69]. In mats and pool fingers, they were represented by numerous creeping or gliding flagellates (Sandona, Cercomonas, Flectomonas), sessile flagellates (Bicoecaceae, Choanoflagellatea), and amoeboid protists (Copromyxa, Filamoeba, Synchromophyceae). In the water communities, there was a large abundance of actively swimming flagellates (Spumella elongata) and ciliates (Colpoda, Halteria), as well as flagellates with ambiguous behavior (Segregatospumella, some Bicoecaceae). The bacterivorous protists are in turn consumed by parasitic (Apicomplexa, Peronosporomycetes, Phytomyxea, Perkinsida, Rozellomycota) and predatory (Orciraptor, Viridiraptor, Allantion, Colpodellida) protists. Finally, the remains of organic matter are consumed by fungi (Basidiomycota, Chytridiomycota, Ascomycota) and fungi-like protists (Labyrinthulomycetes). Similar food chains have been found in groundwater and aquifers [50,69]. The food chain of Blue Lake, located in the illuminated area, differed from the food chain of Distant Lake, with a large proportion of phototrophic protists, mainly Ochromonadales and green algae.
We believe that this study enriches our knowledge of protist communities in karst caves. However, it has some limitations that should be noted here. First of all, the method of DNA metabarcoding based on the high-throughput sequencing amplicons of SSU rRNA genes is the best choice for the evaluation of total richness, phylogenetic diversity, and community dynamics of protists [70]. However, this approach provides information on short target regions of the gene that have insufficient discrimination power for the identification of protists at the species and genus level [70]. The identification of protist species and genera using the method of alignment against gene databases and the search for the closest homolog is often complicated due to blurry limits between genetic intra- and interspecific variation for most taxa [71,72]. The deposition of sequences with unreliable and often even unpublished taxonomic descriptions significantly complicates the systematic classification of protist sequences. For instance, some genera with high diversity and complex taxonomic structure, such as Spumella, cannot be identified reliably using current databases due to numerous deposited organisms formerly identified as Spumella, but after a thorough examination, they were described as new genera, e.g., Segregatospumella, Poteriospumella, etc., and species [73,74]. Similar issues occur with certain genera, e.g., Oikomonas, which have never been characterized according to the current taxonomic standards in protistology.
Another limitation of DNA metabarcoding is insufficient compliance of its results with estimates of the number and biomass of living protists [70]. Metabarcoding provides data on the relative abundance of organisms only, whereas for ecological research, the number and biomass of protists are usually the most crucial. These parameters may currently be estimated with direct microscopic observation and size measurement. The inability to distinguish active vegetative cells from dormant cysts, which are especially widespread in soils and aquifers [69], is an issue of metabarcoding that cannot be resolved with genetic methods. However, the most obvious way to resolve the main limitations of DNA metabarcoding involves a combination of this method with RNA metabarcoding, which provides information on metabolically active, rather than dormant, parts of the community, as well as the use of classical microscopic examination and cultivation methods to reveal the number, biomass, and morphological characteristics of protist species [70].

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/d16090526/s1, Figure S1: Rarefaction curves of sequenced DNA libraries; Figure S2: Relative abundances of the top 20 protist supergroups in the communities of Shulgan-Tash Cave; Figure S3: Relative abundances of the top 20 protist divisions in the communities of Shulgan-Tash Cave; Figure S4. Relative abundances of the top 20 protist subdivisions in the communities of Shulgan-Tash Cave; Figure S5. Relative abundances of the top 20 protist classes in the communities of Shulgan-Tash Cave; Figure S6. Relative abundances of the top 20 protist orders in the communities of Shulgan-Tash Cave; Figure S7. Relative abundances of the top 20 protist families in the communities of Shulgan-Tash Cave; Figure S8. Relative abundances of the top 20 protist genera in the communities of Shulgan-Tash Cave; Figure S9. Top 40 of protist taxonomic groups with significant differential abundance; Figure S10. Results of linear discriminant analysis effect size (LEfSe) analysis; Figure S11. The relative abundance of the protist functional groups; Figure S12. Venn diagram of the total amount of ASVs in Distant Lake and Blue Lake and those shared by both sample groups; Table S1: ASV table with taxonomy for ASVs remaining after bioinformatic treatment of 18S rRNA amplicons obtained from samples of Shulgan-Tash Cave biotopes and sequenced on MiSeq (Illumina); Table S2: Dominant protist taxa (ASVs) represented by more than 1% of reads; Table S3: Results of linear discriminant analysis effect size (LEfSe) analysis.

Author Contributions

Conceptualization, Y.V.G. and N.E.G.; methodology, N.E.G.; software, A.S.B.; formal analysis, N.E.G., M.A.N. and A.S.B.; investigation, N.E.G., M.A.N., O.Y.C., L.Y.K. and E.I.S.; resources, Y.V.G., E.I.S. and A.O.P.; data curation, A.S.B.; writing—original draft preparation, N.E.G., M.A.N., A.S.B., O.Y.C., Y.V.G. and A.O.P.; writing—review and editing, N.E.G., M.A.N., A.S.B., O.Y.C., L.Y.K., E.I.S., Y.V.G. and A.O.P.; visualization, A.S.B. and O.Y.C.; supervision, Y.V.G.; project administration, N.E.G. and A.O.P.; funding acquisition, E.I.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The 18S rRNA gene sequencing data have been submitted to NCBI under BioProject PRJNA1136168 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1136168, accessed on 1 July 2024), Sequence Read Archive (SRA) study SRP520052 (https://trace.ncbi.nlm.nih.gov/Traces?study=SRP520052, accessed on 1 July 2024), and BioSamples SAMN42504430-SAMN42504445 accessed on July 2024.

Acknowledgments

The authors appreciate the Strategic Academic Leadership Program (PRIORITY-2030) for the support, Kazan Federal University. The authors thank the State Nature Reserve “Shulgan-Tash”, represented by Director Michael Kosarev, and the Center for Protection and Use of Objects of Cultural Heritage, represented by Director Dannir Gaynullin, for their long-term support of research in Shulgan-Tash Cave. The authors also express their sincere gratitude to Denis Karimov for his technical assistance with sampling. Y.V.G. expresses gratitude for administrative support within the framework of the state assignment to the Federal Research Center “Kazan Scientific Center of the Russian Academy of Sciences”.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sallstedt, T.; Ivarsson, M.; Lundberg, J.; Sjöberg, R.; Romaní, J.R.V. Speleothem and Biofilm Formation in a Granite/dolerite Cave, Northern Sweden. Int. J. Speleol. 2014, 43, 7. [Google Scholar] [CrossRef]
  2. Brar, A.K.; Bergmann, D. Culture-Based Analysis of “Cave Silver” Biofilms on Rocks in the Former Homestake Mine in South Dakota, USA. Int. J. Speleol. 2019, 48, 3. [Google Scholar] [CrossRef]
  3. Reboul, G.; Moreira, D.; Bertolino, P.; Hillebrand-Voiculescu, A.M.; López-García, P. Microbial Eukaryotes in the Suboxic Chemosynthetic Ecosystem of Movile Cave, Romania. Environ. Microbiol. Rep. 2019, 11, 464–473. [Google Scholar] [CrossRef] [PubMed]
  4. Gogoleva, N.; Chervyatsova, O.; Balkin, A.; Kuzmina, L.; Shagimardanova, E.; Kiseleva, D.; Gogolev, Y. Microbial Tapestry of the Shulgan-Tash Cave (Southern Ural, Russia): Influences of Environmental Factors on the Taxonomic Composition of the Cave Biofilms. Environ. Microbiome 2023, 18, 82. [Google Scholar] [CrossRef]
  5. Morse, K.V.; Richardson, D.R.; Brown, T.L.; Vangundy, R.D.; Cahoon, A.B. Longitudinal Metabarcode Analysis of Karst Bacterioplankton Microbiomes Provide Evidence of Epikarst to Cave Transport and Community Succession. PeerJ 2021, 9, e10757. [Google Scholar] [CrossRef] [PubMed]
  6. Thompson, B.; Richardson, D.; Vangundy, R.; Cahoon, A.B. Metabarcoding Comparison of Prokaryotic Microbiomes from Appalachian Karst Caves to Surface Soils in Southwest Virginia, USA. J. Caves Karst Stud. 2019, 81, 244–253. [Google Scholar] [CrossRef]
  7. Stomeo, F.; Portillo, M.C.; Gonzalez, J.M.; Laiz, L.; Saiz-Jimenez, C. Pseudonocardia in White Colonizations in Two Caves with Paleolithic Paintings. Int. Biodeterior. Biodegrad. 2008, 62, 483–486. [Google Scholar] [CrossRef]
  8. Spilde, M.N.; Northup, D.E.; Caimi, N.A.; Boston, P.J.; Stone, F.D.; Smith, S. Microbial mat communities in Hawaiian lava caves. Geol. Soc. Am. Abstr. Programs 2016, 48, 283965. [Google Scholar] [CrossRef]
  9. Gonzalez-Pimentel, J.L.; Martin-Pozas, T.; Jurado, V.; Miller, A.Z.; Caldeira, A.T.; Fernandez-Lorenzo, O.; Sanchez-Moral, S.; Saiz-Jimenez, C. Prokaryotic Communities from a Lava Tube Cave in La Palma Island (Spain) Are Involved in the Biogeochemical Cycle of Major Elements. PeerJ 2021, 9, e11386. [Google Scholar] [CrossRef]
  10. Sherr, E.B.; Sherr, B.F. Phagotrophic Protists: Central Roles in Microbial Food Webs. In Aquatic Microbial Ecology and Biogeochemistry: A Dual Perspective; Glibert, P., Kana, T., Eds.; Springer: Cham, Switzerland, 2016. [Google Scholar] [CrossRef]
  11. Garner, R.E.; Kraemer, S.A.; Onana, V.E.; Huot, Y.; Gregory-Eaves, I.; Walsh, D.A. Protist Diversity and Metabolic Strategy in Freshwater Lakes Are Shaped by Trophic State and Watershed Land Use on a Continental Scale. mSystems 2022, 7, e0031622. [Google Scholar] [CrossRef]
  12. Gittleson, S.M.; Hoover, R.L. Cavernicolous Protozoa: Review of the Literature and New Studies in Mammoth Cave, Kentucky. Ann. Speleol. 1969, 24, 737–776. [Google Scholar]
  13. Gittleson, S.M.; Hoover, R.L. Protozoa of Underground Waters in Caves. Ann. Econ. Soc. Civilis. 1970, 25, 91–106. [Google Scholar]
  14. Hill, B.F.; Small, E.B.; Iliffe, T.M. Euplotes Iliffei n. Sp.: A New Species of Euplotes (Ciliophora, Hypotrichida) from the Marine Caves of Bermuda. J. Wash. Acad. Sci. 1986, 76, 244–249. [Google Scholar]
  15. Walochnik, J.; Mulec, J. Free-Living Amoebae in Carbonate Precipitating Microhabitats of Karst Caves and a New Vahlkampfiid Amoeba, Allovahlkampfia Spelaea Gen. Nov., Sp. Nov. Acta Protozool. 2009, 48, 25–33. [Google Scholar]
  16. Sigala-Regalado, I.; Mayen-Estrada, R.; Morales-Malacara, J. Spatial and Temporal Distribution of Protozoa at Cueva de Los Riscos, Querétaro, México. J. Caves. Karst Stud. 2011, 73, 55–62. [Google Scholar] [CrossRef]
  17. Bharti, D.; Kumar, S.; Buonanno, F.; Ortenzi, C.; Montanari, A.; Quintela-Alonso, P.; La Terza, A. Free Living Ciliated Protists from the Chemoautotrophic Cave Ecosystem of Frasassi (Italy). SBC 2022, 44, 167–198. [Google Scholar] [CrossRef]
  18. Zepeda Mendoza, M.L.; Lundberg, J.; Ivarsson, M.; Campos, P.; Nylander, J.A.A.; Sallstedt, T.; Dalen, L. Metagenomic Analysis from the Interior of a Speleothem in Tjuv-Ante’s Cave, Northern Sweden. PLoS ONE 2016, 11, e0151577. [Google Scholar] [CrossRef] [PubMed]
  19. Kajan, K.; Cukrov, N.; Cukrov, N.; Bishop-Pierce, R.; Orlić, S. Microeukaryotic and Prokaryotic Diversity of Anchialine Caves from Eastern Adriatic Sea Islands. Microb. Ecol. 2022, 83, 257–270. [Google Scholar] [CrossRef]
  20. Tsouggou, N.; Oikonomou, A.; Papadimitriou, K.; Skandamis, P.N. 16S and 18S rDNA Amplicon Sequencing Analysis of Aesthetically Problematic Microbial Mats on the Walls of the Petralona Cave: The Use of Essential Oils as a Cleaning Method. Microorganisms 2023, 11, 2681. [Google Scholar] [CrossRef]
  21. del Campo, J.; Massana, R. Emerging Diversity within Chrysophytes, Choanoflagellates and Bicosoecids Based on Molecular Surveys. Protist 2011, 162, 435–448. [Google Scholar] [CrossRef]
  22. Harder, C.B.; Rønn, R.; Brejnrod, A.; Bass, D.; Al-Soud, W.A.; Ekelund, F. Local Diversity of Heathland Cercozoa Explored by in-Depth Sequencing. ISME J. 2016, 10, 2488–2497. [Google Scholar] [CrossRef] [PubMed]
  23. Alonso, L.; Pommier, T.; Kaufmann, B.; Dubost, A.; Chapulliot, D.; Doré, J.; Douady, C.J.; Moënne-Loccoz, Y. Anthropization Level of Lascaux Cave Microbiome Shown by Regional-Scale Comparisons of Pristine and Anthropized Caves. Mol. Ecol. 2019, 28, 3383–3394. [Google Scholar] [CrossRef] [PubMed]
  24. Dobàt, K. “Höhlenalgen” Bedrohen Die Eiszeitmalereien von Lascaux. Die Höhle 1963, 14, 41–45. [Google Scholar]
  25. Bontemps, Z.; Alonso, L.; Pommier, T.; Hugoni, M.; Moënne-Loccoz, Y. Microbial Ecology of Tourist Paleolithic Caves. Sci. Total Environ. 2022, 816, 151492. [Google Scholar] [CrossRef] [PubMed]
  26. Lyakhnitsky, Y.; Minnikov, O.; Yushko, A. Drawings and Signs of Shulgan-Tash (Kapova) Cave, Catalogue of Images; Kitap: Ufa, Russia, 2013; ISBN 9785295058479. [Google Scholar]
  27. Dublyansky, Y.; Moseley, G.E.; Lyakhnitsky, Y.; Cheng, H.; Edwards, L.R.; Scholz, D.; Koltai, G.; Spötl, C. Late Palaeolithic Cave Art and Permafrost in the Southern Ural. Sci. Rep. 2018, 8, 12080. [Google Scholar] [CrossRef]
  28. Kuzmina, L.Y.; Galimzyanova, N.F.; Chervyatsova, O.Y.; Sailfullina, N.M.; Kapralov, S.A.; Ryabova, A.S. Biogenous Fouling in Shulgan-Tash Cave (kapova, Southern Urals) and Factors Influencing on Their Expansion. Èkobioteh 2019, 2, 128–142. [Google Scholar] [CrossRef]
  29. Galimzianova, N.F.; Gilvanova, E.A.; Ryabova, A.S.; Guvatova, Z.G.; Kudryavtseva, A.V.; Melentiev, A.I. Phylogenetic Diversity of Prokaryotes in Microbial Communities Inhabiting Rock Surfaces of Shulgan-Tash (kapova) Cave, Southern Urals. Èkobioteh 2020, 3, 298–304. [Google Scholar] [CrossRef]
  30. Chervyatsova, O.Y.; Potapov, S.S.; Kuz’mina, L.Y.; Leonova, L.V. Subaqueous stalactoids in the Dal’nee Verkhnee Lake of the Shilgan-Tash Cave (Southern Urals). News Ural State Min. Univ. 2018, 2, 20–25. [Google Scholar] [CrossRef]
  31. Chervyatsova, O.; Grigoryev, N.; Kuzmina, L.; Dublyansky, Y.; Akhmedyanov, R.; Potapov, S. Documentations of the Sudden Karst Lake Drainage Event in the Shulgan River Valley (Southern Ural, Russia). J. Cave Karst Stud. 2024; in press. [Google Scholar]
  32. Svoyskiy, Y.M.; Romanenko, E.V.; Grigoryev, N.N.; Levanova, E.S. Experience in Documenting Shul’gan-Tash (Kapova) Cave and the Surrounding Landscape with Modern Methods. Brief Commun. Inst. Archaeol. 2020, 261, 67–81. [Google Scholar] [CrossRef]
  33. Melim, A.; Shinglman, K.M.; Boston, P.J.; Northup, D.E.; Spilde, M.N.; Michael Queen, J. Evidence for Microbial Involvement in Pool Finger Precipitation, Hidden Cave, New Mexico. Geomicrobiol. J. 2001, 18, 311–329. [Google Scholar] [CrossRef]
  34. Simon, M.; Jardillier, L.; Deschamps, P.; Moreira, D.; Restoux, G.; Bertolino, P.; López-García, P. Complex Communities of Small Protists and Unexpected Occurrence of Typical Marine Lineages in Shallow Freshwater Systems. Environ. Microbiol. 2015, 17, 3610–3627. [Google Scholar] [CrossRef] [PubMed]
  35. Bower, S.M.; Carnegie, R.B.; Goh, B.; Jones, S.R.; Lowe, G.J.; Mak, M.W. Preferential PCR Amplification of Parasitic Protistan Small Subunit rDNA from Metazoan Tissues. J. Eukaryot. Microbiol. 2004, 51, 325–332. [Google Scholar] [CrossRef] [PubMed]
  36. Martin, M. Cutadapt Removes Adapter Sequences from High-Throughput Sequencing Reads. EMBnet J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  37. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-Resolution Sample Inference from Illumina Amplicon Data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef]
  38. Guillou, L.; Bachar, D.; Audic, S.; Bass, D.; Berney, C.; Bittner, L.; Boutte, C.; Burgaud, G.; de Vargas, C.; Decelle, J.; et al. The Protist Ribosomal Reference Database (PR2): A Catalog of Unicellular Eukaryote Small Sub-Unit rRNA Sequences with Curated Taxonomy. Nucleic Acids Res. 2013, 41, D597–D604. [Google Scholar] [CrossRef] [PubMed]
  39. McMurdie, P.J.; Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef]
  40. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; 2197-5736; Springer: Cham, Switzerland, 2016; ISBN 9783319242774. [Google Scholar]
  41. Xu, S.; Zhan, L.; Tang, W.; Wang, Q.; Dai, Z.; Zhou, L.; Feng, T.; Chen, M.; Wu, T.; Hu, E.; et al. MicrobiotaProcess: A Comprehensive R Package for Deep Mining Microbiome. Innovation 2023, 4, 100388. [Google Scholar] [CrossRef]
  42. Lahti, L.; Shetty, S. Tools for Microbiome Analysis in R. Available online: https://www.bioconductor.org/packages/devel/bioc/vignettes/microbiome/inst/doc/vignette.html (accessed on 19 July 2023).
  43. Lozupone, C.; Knight, R. UniFrac: A New Phylogenetic Method for Comparing Microbial Communities. Appl. Environ. Microbiol. 2005, 71, 8228–8235. [Google Scholar] [CrossRef]
  44. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic Biomarker Discovery and Explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef]
  45. Cao, Y.; Dong, Q.; Wang, D.; Zhang, P.; Liu, Y.; Niu, C. microbiomeMarker: An R/Bioconductor Package for Microbiome Marker Identification and Visualization. Bioinformatics 2022, 38, 4027–4029. [Google Scholar] [CrossRef]
  46. Singer, D.; Seppey, C.V.W.; Lentendu, G.; Dunthorn, M.; Bass, D.; Belbahri, L.; Blandenier, Q.; Debroas, D.; de Groot, G.A.; de Vargas, C.; et al. Protist Taxonomic and Functional Diversity in Soil, Freshwater and Marine Ecosystems. Environ. Int. 2021, 146, 106262. [Google Scholar] [CrossRef] [PubMed]
  47. Carvalho da Silva, V.; Fernandes, N. Protist Taxonomic and Functional Diversity in Aquatic Ecosystems of the Brazilian Atlantic Forest. PeerJ 2023, 11, e15762. [Google Scholar] [CrossRef] [PubMed]
  48. Howe, A.T.; Bass, D.; Chao, E.E.; Cavalier-Smith, T. New Genera, Species, and Improved Phylogeny of Glissomonadida (Cercozoa). Protist 2011, 162, 710–722. [Google Scholar] [CrossRef] [PubMed]
  49. Howe, A.T.; Bass, D.; Vickerman, K.; Chao, E.E.; Cavalier-Smith, T. Phylogeny, Taxonomy, and Astounding Genetic Diversity of Glissomonadida Ord. Nov., the Dominant Gliding Zooflagellates in Soil (Protozoa: Cercozoa). Protist 2009, 160, 159–189. [Google Scholar] [CrossRef] [PubMed]
  50. Groult, B.; St-Jean, V.; Lazar, C.S. Linking Groundwater to Surface Discharge Ecosystems: Archaeal, Bacterial, and Eukaryotic Community Diversity and Structure in Quebec (Canada). Microorganisms 2023, 11, 1674. [Google Scholar] [CrossRef]
  51. Baković, N.; Siemensma, F.; Puljas, S.; Baković, R.; Ozimec, R.; Ostojić, A.; Mesić, Z. First Data on Testate Amoebae Associated with the Endemic Cave Bivalve Congeria Jalzici Morton & Bilandžija, 2013 with a Description of Psammonobiotus Dinarica Sp. Nov. SBC 2023, 45, 53–74. [Google Scholar] [CrossRef]
  52. Cahoon, A.B.; VanGundy, R.D. Alveolates (dinoflagellates, Ciliates and Apicomplexans) and Rhizarians Are the Most Common Microbial Eukaryotes in Temperate Appalachian Karst Caves. Environ. Microbiol. Rep. 2022, 14, 538–548. [Google Scholar] [CrossRef]
  53. Kofoid, C.A. The Plankton of Echo River, Mammoth Cave. Trans. Am. Microsc. Soc. 1900, 21, 113–126. [Google Scholar] [CrossRef]
  54. Jungblut, A.D.; Vincent, W.F.; Lovejoy, C. Eukaryotes in Arctic and Antarctic Cyanobacterial Mats. FEMS Microbiol. Ecol. 2012, 82, 416–428. [Google Scholar] [CrossRef]
  55. Vimercati, L.; Darcy, J.L.; Schmidt, S.K. The Disappearing Periglacial Ecosystem atop Mt. Kilimanjaro Supports Both Cosmopolitan and Endemic Microbial Communities. Sci. Rep. 2019, 9, 10676. [Google Scholar] [CrossRef]
  56. Takishita, K.; Miyake, H.; Kawato, M.; Maruyama, T. Genetic Diversity of Microbial Eukaryotes in Anoxic Sediment around Fumaroles on a Submarine Caldera Floor Based on the Small-Subunit rDNA Phylogeny. Extremophiles 2005, 9, 185–196. [Google Scholar] [CrossRef]
  57. Malygina, A.; Balkin, A.; Polyakova, E.; Stefanov, S.; Potekhin, A.; Gogoleva, N. Taxonomic Diversity of the Microbial Biofilms Collected along the Thermal Streams on Kunashir Island. Ecologies 2023, 4, 106–123. [Google Scholar] [CrossRef]
  58. Grant, J.; Tekle, Y.I.; Anderson, O.R.; Patterson, D.J.; Katz, L.A. Multigene Evidence for the Placement of a Heterotrophic Amoeboid Lineage Leukarachnion Sp. among Photosynthetic Stramenopiles. Protist 2009, 160, 376–385. [Google Scholar] [CrossRef]
  59. Abdullin, S.R.; Sharipova, M.Y. Algoflora of the Shulgan-Tash Cave (Kapova) in Various Years. Bull. Bashkir Univ. 2005, 3, 49–50. [Google Scholar]
  60. Claus, G. Algae and Their Mode of Life in the Baradla Cave at Aggtelek II. Int. J. Speleol. 1964, 1, 2. [Google Scholar] [CrossRef]
  61. Abdullin, S.R.; Bagmet, V.B. Mixotrophy of Cyanobacteria and Algae in Cave Conditions. J. Gen. Biol. 2016, 77, 54–62. [Google Scholar]
  62. Letcher, P.M.; Powell, M.J. A Taxonomic Summary and Revision of Rozella (Cryptomycota). IMA Fungus 2018, 9, 383–399. [Google Scholar] [CrossRef]
  63. Jurado, V.; del Rosal, Y.; Gonzalez-Pimentel, J.L.; Hermosin, B.; Saiz-Jimenez, C. Biological Control of Phototrophic Biofilms in a Show Cave: The Case of Nerja Cave. NATO Adv. Sci. Inst. Ser. E Appl. Sci. 2020, 10, 3448. [Google Scholar] [CrossRef]
  64. Karunarathna, S.C.; Dong, Y.; Karasaki, S.; Tibpromma, S.; Hyde, K.D.; Lumyong, S.; Xu, J.; Sheng, J.; Mortimer, P.E. Discovery of Novel Fungal Species and Pathogens on Bat Carcasses in a Cave in Yunnan Province, China. Emerg. Microbes Infect. 2020, 9, 1554–1566. [Google Scholar] [CrossRef]
  65. Biagioli, F.; Coleine, C.; Piano, E.; Nicolosi, G.; Poli, A.; Prigione, V.; Zanellati, A.; Varese, C.; Isaia, M.; Selbmann, L. Microbial Diversity and Proxy Species for Human Impact in Italian Karst Caves. Sci. Rep. 2023, 13, 689. [Google Scholar] [CrossRef]
  66. Martin-Pozas, T.; Nováková, A.; Jurado, V.; Cuezva, S.; Fernandez-Cortes, A.; Saiz-Jimenez, C.; Sanchez-Moral, S. A Second Fungal Outbreak in Castañar Cave, Spain, Discloses the Fragility of Subsurface Ecosystems. Microb. Ecol. 2024, 87, 53. [Google Scholar] [CrossRef]
  67. Lorch, J.M.; Palmer, J.M.; Vanderwolf, K.J.; Schmidt, K.Z.; Verant, M.L.; Weller, T.J.; Blehert, D.S. Malassezia Vespertilionis Sp. Nov.: A New Cold-Tolerant Species of Yeast Isolated from Bats. Persoonia 2018, 41, 56–70. [Google Scholar] [CrossRef]
  68. Martin-Sanchez, P.M.; Nováková, A.; Bastian, F.; Alabouvette, C.; Saiz-Jimenez, C. Use of Biocides for the Control of Fungal Outbreaks in Subterranean Environments: The Case of the Lascaux Cave in France. Environ. Sci. Technol. 2012, 46, 3762–3770. [Google Scholar] [CrossRef] [PubMed]
  69. Novarino, G.; Warren, A.; Butler, H.; Lambourne, G.; Boxshall, A.; Bateman, J.; Kinner, N.E.; Harvey, R.W.; Mosse, R.A.; Teltsch, B. Protistan Communities in Aquifers: A Review. FEMS Microbiol. Rev. 1997, 20, 261–275. [Google Scholar] [CrossRef] [PubMed]
  70. Santoferrara, L.; Burki, F.; Filker, S.; Logares, R.; Dunthorn, M.; McManus, G.B. Perspectives from Ten Years of Protist Studies by High-Throughput Metabarcoding. J. Eukaryot. Microbiol. 2020, 67, 612–622. [Google Scholar] [CrossRef]
  71. Bachy, C.; Dolan, J.R.; López-García, P.; Deschamps, P.; Moreira, D. Accuracy of Protist Diversity Assessments: Morphology Compared with Cloning and Direct Pyrosequencing of 18S rRNA Genes and ITS Regions Using the Conspicuous Tintinnid Ciliates as a Case Study. ISME J. 2013, 7, 244–255. [Google Scholar] [CrossRef] [PubMed]
  72. Decelle, J.; Romac, S.; Sasaki, E.; Not, F.; Mahé, F. Intracellular Diversity of the V4 and V9 Regions of the 18S rRNA in Marine Protists (radiolarians) Assessed by High-Throughput Sequencing. PLoS ONE 2014, 9, e104297. [Google Scholar] [CrossRef] [PubMed]
  73. Grossmann, L.; Bock, C.; Schweikert, M.; Boenigk, J. Small but Manifold—Hidden Diversity in “Spumella-like Flagellates”. J. Eukaryot. Microbiol. 2016, 63, 419–439. [Google Scholar] [CrossRef]
  74. Jeong, M.; Kim, J.I.; Nam, S.W.; Shin, W. Molecular Phylogeny and Taxonomy of the Genus Spumella (Chrysophyceae) Based on Morphological and Molecular Evidence. Front. Plant Sci. 2021, 12, 758067. [Google Scholar] [CrossRef]
Figure 2. Distant Lake sampling sites: (a) Distant Lake; (b) the outflowing creek with filamentous microbial mats; (c) microbial mats; (d) pool fingers.
Figure 2. Distant Lake sampling sites: (a) Distant Lake; (b) the outflowing creek with filamentous microbial mats; (c) microbial mats; (d) pool fingers.
Diversity 16 00526 g002
Figure 3. Alpha diversity indices of the protist communities in the sites of Shulgan-Tash Cave and the outer Blue Lake: Chao 1, Shannon, Gini-Simpson, and Pielou’s evenness. The points on the plot correspond to the individual samples and are colored according to the different sites: PF—pool finger biofilm; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Figure 3. Alpha diversity indices of the protist communities in the sites of Shulgan-Tash Cave and the outer Blue Lake: Chao 1, Shannon, Gini-Simpson, and Pielou’s evenness. The points on the plot correspond to the individual samples and are colored according to the different sites: PF—pool finger biofilm; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Diversity 16 00526 g003
Figure 4. 2D plots based on the results of principal coordinate analysis (PCoA) of the protist communities using Bray–Curtis (left) and the weighted UniFrac (right) distances. PCo1 (Axis 1) and PCo2 (Axis 2) explained 33.65% and 25.8% of the protist community variance at the ASV level, respectively. PF—pool finger biofilms; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Figure 4. 2D plots based on the results of principal coordinate analysis (PCoA) of the protist communities using Bray–Curtis (left) and the weighted UniFrac (right) distances. PCo1 (Axis 1) and PCo2 (Axis 2) explained 33.65% and 25.8% of the protist community variance at the ASV level, respectively. PF—pool finger biofilms; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Diversity 16 00526 g004
Figure 5. The most abundant protist taxa in the communities of water sites in Shulgan-Tash Cave. Circle size indicates the inferred relative abundance based on amplicon numbers (in %). PF—pool finger biofilm; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Figure 5. The most abundant protist taxa in the communities of water sites in Shulgan-Tash Cave. Circle size indicates the inferred relative abundance based on amplicon numbers (in %). PF—pool finger biofilm; CEW—lake water at the point of inflow into Distant Lake; WPF—water surrounding pool fingers; WDL—water in the middle of Distant Lake; MDL—mats of Distant Lake; BL—water from Blue Lake.
Diversity 16 00526 g005
Figure 6. The relative abundances of the protist functional groups in communities of Blue Lake and Distant Lake (Shulgan-Tash Cave). The colors correspond to the supergroup names.
Figure 6. The relative abundances of the protist functional groups in communities of Blue Lake and Distant Lake (Shulgan-Tash Cave). The colors correspond to the supergroup names.
Diversity 16 00526 g006
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

Gogoleva, N.E.; Nasyrova, M.A.; Balkin, A.S.; Chervyatsova, O.Y.; Kuzmina, L.Y.; Shagimardanova, E.I.; Gogolev, Y.V.; Plotnikov, A.O. Flourishing in Darkness: Protist Communities of Water Sites in Shulgan-Tash Cave (Southern Urals, Russia). Diversity 2024, 16, 526. https://doi.org/10.3390/d16090526

AMA Style

Gogoleva NE, Nasyrova MA, Balkin AS, Chervyatsova OY, Kuzmina LY, Shagimardanova EI, Gogolev YV, Plotnikov AO. Flourishing in Darkness: Protist Communities of Water Sites in Shulgan-Tash Cave (Southern Urals, Russia). Diversity. 2024; 16(9):526. https://doi.org/10.3390/d16090526

Chicago/Turabian Style

Gogoleva, Natalia E., Marina A. Nasyrova, Alexander S. Balkin, Olga Ya. Chervyatsova, Lyudmila Yu. Kuzmina, Elena I. Shagimardanova, Yuri V. Gogolev, and Andrey O. Plotnikov. 2024. "Flourishing in Darkness: Protist Communities of Water Sites in Shulgan-Tash Cave (Southern Urals, Russia)" Diversity 16, no. 9: 526. https://doi.org/10.3390/d16090526

APA Style

Gogoleva, N. E., Nasyrova, M. A., Balkin, A. S., Chervyatsova, O. Y., Kuzmina, L. Y., Shagimardanova, E. I., Gogolev, Y. V., & Plotnikov, A. O. (2024). Flourishing in Darkness: Protist Communities of Water Sites in Shulgan-Tash Cave (Southern Urals, Russia). Diversity, 16(9), 526. https://doi.org/10.3390/d16090526

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