Next Article in Journal
Reproductive Biology of the Golden Cuttlefish Sepia esculenta (Cephalopoda, Sepiida)
Next Article in Special Issue
Free-Living and Particle-Associated Microbial Communities of Lake Baikal Differ by Season and Nutrient Intake
Previous Article in Journal
Space Use and Movements of Southeastern Breeding Double-Crested Cormorants (Nannopterum auritum) in the United States
Previous Article in Special Issue
Under-Ice Development of Silica-Scaled Chrysophytes with Different Trophic Mode in Two Ultraoligotrophic Lakes of Yakutia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microeukaryotic Communities of the Long-Term Ice-Covered Freshwater Lakes in the Subarctic Region of Yakutia, Russia

1
Limnological Institute, Siberian Branch of the Russian Academy of Sciences, 3 Ulan-Batorskaya Street, Irkutsk 664033, Russia
2
Institute for Biological Problems of Cryolithozone, Siberian Branch of the Russian Academy of Sciences, 41 Lenin Ave, Yakutsk 677980, Russia
*
Author to whom correspondence should be addressed.
Diversity 2023, 15(3), 454; https://doi.org/10.3390/d15030454
Submission received: 29 January 2023 / Revised: 10 March 2023 / Accepted: 16 March 2023 / Published: 18 March 2023
(This article belongs to the Special Issue Diversity and Ecology of Freshwater Lake Microbial Communities)

Abstract

:
Currently, microeukaryotic communities of the freshwater arctic and subarctic ecosystems are poorly studied. Still, these are of considerable interest due to the species biogeography and autecology as well as global climate change. Here, we used high-throughput 18S rRNA amplicon sequencing to study the microeukaryotic communities of the large subarctic freshwater lakes Labynkyr and Vorota in Yakutia, Russia, during the end of the ice cover period, from April to June. By applying the statistical methods, we coupled the microeukaryotic community structure profiles with available discrete factor variables and hydrophysical, hydrochemical, and environmental parameters. The sub-ice layer and the water column communities were differentiated due to the temporal change in environmental conditions, particularly temperature regime and electric conductivity. Additionally, the community composition of unicellular eukaryotes in lakes Labynkyr and Vorota was changing due to seasonal environmental factors, with these alterations having similar patterns in both sites. We suggest the community developed in the sub-ice layer in April serves as a primer for summer freshwater microeukaryotes. Our results extend the current knowledge on the community composition and seasonal succession of unicellular eukaryotes within subarctic freshwater ecosystems.

1. Introduction

Arctic and subarctic freshwater ecosystems are apparently insufficiently studied with respect to microeukaryotic community composition, as well as its spatial variability and temporal dynamics [1]. However, this problem is not only of significant interest with regard to biogeography, autecology and paleohistory of individual species, but is also important for environmental monitoring, for the study of under-ice water ecosystems and in the context of global climate changes [2]. In this respect, two large freshwater subarctic lakes, Labynkyr (LL) and Vorota (LV), which are located in close proximity to the cold pole of Eurasia, are very valuable objects. Until recently, communities of unicellular eukaryotes in these lakes were not analyzed by metagenomics or high-throughput sequencing (HTS). A typical subarctic freshwater ecosystem is characterized by a low temperature, not exceeding 10 °C even in the warmest summer period, decreased concentrations of dissolved inorganic nutrients, low primary production, and oligotrophic state [3,4]. LL and LV are covered by ice for nearly two-thirds of the year; therefore, changes in the community composition between the ice cover and the open water periods are of particular interest. Recent studies of these lakes have shown high diversity and unique composition of microalgae by optical and electron microscopy [5,6,7,8,9] and bacteria by high-throughput sequencing [10].
Extensive knowledge of the under-ice environmental conditions in lake ecosystems has already been accumulated. This might help us to understand how a particular lake ecosystem will evolve over the next season. The under-ice environmental monitoring and evaluation of seasonal connectivity can improve our understanding of the environmental response to global climate change [11]. However, only a small number of such studies consider microeukaryotic community structure as a proxy for continuous environmental monitoring in freshwater lakes. For instance, the long-term study of phytoplankton in oligotrophic Lake Tovel, Brenta Dolomites, during the ice-covered and open water periods was recently reported [12]. However, investigation of the microeukaryotic communities in oligotrophic mountain, subarctic and arctic lakes is often limited by late spring and summer periods due to the difficulty of accessing the sampling sites in these regions.
Within the ice covered period, the environmental conditions are different between the sub-ice layer adjacent to the bottom ice surface and the upper layer of the water column. Further, conditions in both environments are changing during the ice cover period, which results in continuous differentiation of the under-ice communities of microeukaryotes and bacteria [12,13]. Over the last 10 years, Lake Baikal, the largest freshwater lake on Earth, located 10° to the south from LL and LV in latitudinal direction, has been actively studied with respect to microeukaryotic community features by applying a metagenomics approach and HTS techniques. Particularly, the under-ice phytoplankton communities of Lake Baikal were characterized using microscopy and HTS of SSU amplicons [14,15,16,17,18]. These results provide a valuable background that could be used to explore the microeukaryotic communities of LL/LV during the ice-covered and open water periods.
In this work, we compiled an initial inventory of unicellular eukaryotes in LL and LV during the ice cover period by using the metabarcoding of 18S rRNA and investigated the compositional differences between the sub-ice layer and the water column communities in these lakes. Furthermore, we evaluated the impact of seasonal changes on the community structure. Below, we discuss the consistency of the obtained results and hypothesize on the environmental factors influencing changes in community structure.

2. Materials and Methods

2.1. Sampling and Environmental Parameters

Detailed descriptions of the geographical position, climatic characteristics, physical and chemical environmental parameters of the Lake Vorota and Lake Labynkyr sampling sites were presented by Zakharova et al. [10]. Water samples were collected from the sub-ice layer of the bottom ice surface (SI) and from the water column (WC) during April and May 2016 and June 2017. For WC, two liters of water were sampled from 5 m depth using a Niskin bottle. The SI water samples were collected manually from the bottom ice surface (water layer thickness ~1–2 cm) by divers holding the syringe in the outstretched arm.

2.2. Amplicon Library Preparation, Sequencing, Raw Data Processing and Quality Control

Approximately 2 L out of each of 20 water samples (Table 1) was first prefiltered through a 27 μm mesh and then filtered through 0.2 μm polycarbonate membrane (Whatman Part of GE HealthCare, Chicago, Illinois, USA). Biomass collected on 0.2 μm membrane was washed out of the filters into sterile flasks filled with 5 mL of TE buffer (10 mM Tris-HCl, 1 mM EDTA; pH 7.5) and stored at −20 °C until samples were transported to the laboratory. Then, samples were stored at −80 °C until further analysis. DNA was extracted using lysozyme (1 mg/mL), proteinase K, 10% SDS, and phenol:chloroform:isoamyl alcohol mixture (25:24:1) according to the protocol based on Rusch et al. [19].
The purified genomic DNA was sent to HTS provider for amplicon library preparation and sequencing. Amplicon libraries were prepared with primers V8f (5′-ATAACAGGTCTGTGATGCCCT-3′) and 1510R (5′-CCTTCYGCAGGTTCACCTAC-3′) [20], targeting a broad set of eukaryotic taxa. Amplicon library preparation was performed as described in “Illumina 16S Metagenomic Library Preparation Guide #15044223 Rev. B”. Libraries were analyzed by the Core Centrum “Genomic Technologies, Proteomics and Cell Biology”, ARRIAM (Saint-Petersburg, Russia) using Illumina MiSeq Sequencer with Reagent Kit v3 to obtain 300 bp paired-end reads. DNA samples L4SI05 and L4WC06 were used to generate technical replicates L4SI05-02 and L4WC06-02 with independent amplicon library preparation and sequencing procedures (Table 2).
Analysis of sequencing data was performed in Usearch v.10 (https://www.drive5.com/usearch/ (accessed on 25 January 2023)) [21], Vsearch v.2.9.1 (https://github.com/torognes/vsearch (accessed on 25 January 2023)) [22], and mothur v.1.43.10 (https://mothur.org/ (accessed on 25 January 2023)) [23]. Raw reads were merged, primer sequences were truncated, and contigs were filtered by the expected error threshold 1.0 with usearch (option -fastq_maxee 1.0). Next, sequences were clustered with vsearch command “-cluster_size” at 0.97 identity threshold, and OTUs with less than two reads were discarded. OTU centroids were subjected to chimera filtering by UCHIME-denovo followed by UCHIME-reference of vsearch. Finally, the community composition matrix was generated by remapping of merged and quality-filtered reads to a chimera-free set of OTUs with identity threshold 0.97 (vsearch command “-usearch_global”).

2.3. Bioinformatics, Statistical Analyses, and Data Visualization

OTU sequences were taxonomically classified using Silva v.138 (Bremen, Germany) [24] in mothur with probability cutoff set to 80 and further filtered by taxonomy to drop chloroplast-specific, mitochondria-specific and multicellular eukaryotic OTUs as well as OTUs with Kingdom = “unknown”. Additionally, the taxonomy mapping was performed with the online ACT service https://www.arb-silva.de/aligner (25 February 2023). Results of taxonomic assignment and sequence similarity search are listed in Table S3. To generate the summary statistics on the abundance of phylotypes (Table S2), the raw community composition matrix was standardized by a median of sample sequencing depth, OTUs with total abundance less than 10 were filtered out as low-covered, phylotype counts were aggregated up to the “Class” taxonomic level, which was manually assigned using the results of mothur/ACT taxonomic pipelines (Table S4), and the top-13 abundant classes were selected to draw the histogram.
Statistical analyses were performed in R. The rarefaction curves, ACE index (non-parametric species richness estimators), and Shannon, Simpson, and inverse Simpson indices were calculated to measure the alpha diversity. OTUs with a total abundance below 10 were excluded from the community composition matrix prior to these calculations. Exploratory analyses of community composition were performed using vegan v.2.5-6 [25], phyloseq [26], and pvclust [27] (https://www.R-project.org/ (accessed on 25 January 2023)).
One-way ANOVA and Kruskal–Wallis statistical tests were used to examine the impact of independent factor variables on the discrete/continuous environmental variables (Figure S1) and alpha-diversity metrics (Figure S2) associated with the microeukaryotic community profiles. Two-way ANOVA test was used to estimate the main effect of factors “Month” and “Layer” and the effect of their interaction. We calculated two-way ANOVA statistics with R function “Anova” from package “car”, using the “type-III” test for experiments with unbalanced design. The p-values, which were obtained by testing the independent factor variable (or two of them) versus a set of dependent variables (14 environmental variables and five alpha-diversity metrics), were adjusted using the false discovery rate (FDR) procedure. Both raw and adjusted p-values were reported. May and June groups were merged into “MJ” category when the alpha-diversity metrics were tested by ANOVA and Kruskal–Wallis procedures for the following reasons: (i) the SI June group was presented by a single profile with values falling into the range of SI May profiles; (ii) WC May group was presented by two profiles with values within the range of WC June profiles (Figure S2(C1–D5)). R-package ggpubr v. 0.6.0 (https://cran.r-project.org/web/packages/ggpubr/index.html (accessed on 25 January 2023)) was used to visualize the group-wise distribution of dependent variables.
For exploratory analyses, OTUs with a relative abundance above 0.1% or standardized abundance above 50 were selected to produce a count matrix. OTU counts were transformed with log(x + 1) and subjected to the transformation-based principal coordinate analysis (tb-PCoA). Linear regression of explanatory variables was performed by the envfit function of package vegan, followed by adjustment of permutation-based regression p-values by the Holm procedure. Environmental factors having an adjusted p-value threshold below 0.05 were drawn on the ordination plane. The pairwise distance matrix computed with Bray–Curtis similarity index was used to cluster communities by UPGMA. A robustness of clustering dendrogram was accessed by bootstrap- and au-values in the pvclust package.
The transformation-based redundancy analysis (tb-RDA) was used to evaluate the impact of environmental factors on species composition. OTUs with relative abundance above 0.1% were selected to produce a count matrix. OTU counts were transformed with log(x + 1) and subjected to tb-RDA. Explanatory variables were chosen by a “forward selection” approach, followed by further filtering on p-value and the impact of the model ability to explain the total variance (Table S1).
DESeq2 package [28] was used to evaluate the significance of OTU abundance differences between the groups of samples. The raw community composition matrix of the non-transformed OTU counts, which was generated by the quality control/clustering pipeline described above (usearch-vsearch-taxonomic filtering), was used as DESeq2 input. The grouping of samples for differential abundance testing was established considering the results of exploratory analyses. It included five subsets of microeukaryotic community profiles: (i) LL SI April; (ii) LL SI May; (iii) LL WC April; (iv) LL WC May–June; (v) LV SI/WC June (all LV profiles except V1SI05). Ten contrasts corresponding to pairwise combinations of these five subsets were used for OTU differential abundance analysis. A phyloseq object was transformed to DESeq2 object, followed by estimation of size factors with geometric means of OTU counts, “local” dispersion estimate, and computation of Wald test statistics. The FDR-adjusted p-value and log-fold change (LFC) thresholds were set to 0.05 and log2(1.5), respectively. OTUs indicating significant difference between groups of samples are listed in Table S5.

3. Results

3.1. Environmental Parameters at Sampling Sites

The lake ice was covered by snow at the time of sampling, with its thickness decreasing from 30–35 cm in April to 0.5–5 cm in June. Samples from the LL ice bottom surface were not collected in June, as the process of active ice degradation started by this time, and vertical channels filled with water were formed in the ice. During the sampling period, the water temperature of the SI layer varied between 0.4 and 1.2 °C, while the WC temperature range was 1.2–5.6 °C. The pH value varied between 6.8 and 9.34, electric conductivity was between 8.5 and 55 μS/cm2, and mineral nitrogen concentration ranged between 0.02 and 0.13 mg/L (Table 1).

3.2. Quality Control of the Amplicon Clustering Results and Rarefaction Analysis

Sample V1WC05 was excluded from further analysis due to the lack of coverage for targeted microeukaryotic taxa: almost all OTUs of V1WC05 were classified as Arthropoda-specific sequences. For the remaining 20 samples, 22 amplicon data bunches were generated, including pairs of technical replicates for L4WC06 and L4SI05 (Table 2). After quality control steps, the coverage per sample varied between 7032 and 51,984 microeukaryote-specific reads. Processing of the amplicon sequencing data produced 459 distinct OTUs (3% dissimilarity threshold, minimal OTU coverage above 10 reads), ranging from 73 to 337 OTUs per sample with mean sequence length of 330 bp. Rarefaction curves of the amplicons reached saturation, indicating the diversity of the microeukaryotic V8-V9 SSU fragments was decently recovered, at least under current experimental settings (Figure 1). The exception was V1SI06, which had minimal sample coverage for microeukaryote-specific taxa because of contamination by Arthropoda-specific sequences.

3.3. Factors Influencing the Hydrophysical and, Hydrochemical Characteristics of Samples and the Alpha-Diversity Metrics of Microeukaryotic Communities

ANOVA revealed the phosphate and nitrate anion concentrations, and Nmin were significantly different between LL and LV samples (Table 3, Figure S1(A5,A8,A9)). The phosphate anion concentration was higher in LV, while the nitrate anion and Nmin concentrations were higher in LL. The concentrations of TOC and nitrite anion by LL/LV showed a barely significant response (Figure S1(A7,A10)). A number of hydrochemical, hydrophysical and environmental parameters, specifically pH, phosphate, ammonium, nitrite, Nmin, TOC, TMA and the snow thickness, varied from April to June (Table 3, Figure S1(B1–B14)). The phosphate and nitrite anion concentrations, as well as TMA, gradually increased from April to June (Figure S1(B5,B7,B11)). The ammonium cation, Nmin and TOC concentrations decreased from April to June (Figure S1(B6,B9,B10)). However, the only variable differing between layers was water temperature (Table 3, Figure S1(C1)), which obviously was lower in the SI layer than in the WC layer.
Importantly, the two-factor sorting of environmental variables by the sampling layer and month also produced a set of distinct spatiotemporal patterns. In the WC layer, the pH, concentration of ammonium cation, concentration of nitrite anion, and TMA differed significantly over time (Table 3, Figure S1(D2,D6,D7,D11)), reinforcing the simple temporal response (Figure S1(B2,B6,B7,B11). Similarly, in SI layer, the concentrations of nitrate anion and Nmin decreased over time (Figure S1(D8,D9)). The layer-wise temporal concentrations of TOC basically mirrored its general temporal dynamics (Table 3, Figure S1(B10,D10)), with some lag in the WC layer.
The alpha-diversity metrics were also detected to be different over the sampling layer and month factors: the Shannon, Simpson and inverse Simpson indices were higher in WC profiles (Table 3, Figure S2(C3,C5)), while the species richness and ACE were gradually increased from April to June (Table 3, Figure S2(B1,B2)). The two-factor sorting of alpha-diversity metrics by the sampling layer and month also produced discrete patterns: the species richness and ACE were significantly increased in May–June WC profiles, while the Shannon and inverse Simpson indices were higher in April SI profiles (Table 3, Figure S2(D1,D5)).
With respect to interaction of the sampling layer and month factor variables, two-way ANOVA produced barely significant results for two out of 14 environmental variables: concentration of TOC and TMA (Table 4). At the same time, Shannon and Simpson alpha-diversity metrics displayed good support for the two-factor interaction model, which included the sampling layer and month factors. Additionally, the species richness, ACE, and the inverse Simpson metrics revealed that the “Month:Layer” interaction term of the two-factor model was barely significant (Table 4).

3.4. Overall Community Composition and Comparison of Microeukaryotic Profiles

The generated OTUs were taxonomically assigned to 34 classes of microeukaryotes. This set represented a wide spectrum of unicellular protists. The most abundant were (in decreasing order) the classes Dinophyceae, Chrysophyceae, Ciliophora, Trebouxiophyceae (green algae), Bacillariophyceae, Basidio- and Ascomycota (Fungi), unclassified Chlorophyta, Cercozoa, Cryptophyceae, Prymnesiophyceae (Haptophyta), and Telonena (SAR group) (Figure 2). Ciliophora were abundant in April LL profiles, with decreasing counts in June. Diatoms were predominantly abundant in WC communities of May–June. A distinctive feature of LV profiles was the increased number of Trebouxiophyceae, which were also present in June LL WC, while their abundance was lower. Peaks of Basidiomycota and Cercozoa in LL occurred in April, with Basidiomycota predominantly abundant in the WC layer. On the other hand, Ascomycota predominantly proliferated in May communities.
Exploratory analyses revealed that the major factor variables influencing the profile clustering were the month of sampling (April/May–June) and the sampling layer (SI/WC) (Figure 3). In tb-PCoA, the LL SI April communities were similar to the LL WC April ones (Figure 3A,C, quadrant III). The LL SI May profiles were completely different from others, forming a distinct group (Figure 3A,C, quadrant II). The LL WC microeukaryotic communities of May–June clustered together (Figure 3A,C, quadrant IV).
Analysis of alpha diversity (Table 2 and Table 3) showed that the LL WC microeukaryotic communities of May–June had the highest richness. Noticeably, the diversity of the May–June WC microeukaryotic community profiles was also increased in comparison with the May–June SI profiles. Using the results of differential abundance analysis (Table S5) and OTU abundance heatmap (Figure 4), one may identify the key OTUs that are responsible for this effect. Particularly, OTU146 (family Suessiaceae), OTU153 (Chlorophyta), OTU305 (Chrysophyceae), OTU442 (order Perkinsidae), OTU115 (Chrysophyceae), and OTU576 (order Eustigmatales) were detected in the LL WC May–June profiles, while the rest of the communities lacked these OTUs. Typically, the LL WC May–June profiles were the most similar to those of LL SI April. The latter group was highly similar to LL WC April (Figure 3B,D). The common feature of these groups was the increased abundance of OTUs in the middle part of Figure 4 (top-down from OTU45 to OTU221). Species classified as Perkinsidae, Chrysophyceae, Choanoflagellata, Chlorophyta, Intramacronucleata (orders Spirotrichea and Conthreep), Dinophyceae and Cercozoa were presented in this subset. At the same time, several OTUs from this subset (top-down from OTU56 to OTU221), which were classified as Dinophyceae, Chlorophyta, Chlorodendrophyceae, Intramacronucleata, Prymnesiophyceae, MAST-2, and MAST-12, had increased abundance in the LV June profiles.
OTUs that were present in almost all profiles at top-abundant positions were taxonomically assigned to dinoflagellates (OTU3, OTU50, order Peridiniophyceae—OTU71), diatoms (OTU1), Cryptophyceae (OTU17), Ciliophora (order Spirotrichea—OTU42), Chrysophyceae (OTU49), Cercozoa (OTU15), and various fungi: Cystobasidiomycetes (OTU7), Leotiomycetes (OTU163), Dothideomycetes (OTU97, OTU117), and Tremellomycetes (OTU47).

3.5. Differentially Abundant OTUs of LL Communities in April

A comparatively small subset of differentially abundant OTUs with the relative abundance above 0.1% was found between the April SI and WC microeukaryotic community profiles (Figure 4). Eight OTUs were abundant in SI, while three were abundant in WC (Table S5).
The following OTUs were abundant in the LL SI profiles: OTU3, OTU65, OTU12, OTU158, OTU13, OTU27, OTU42, and OTU109. OTU3 was one of the most abundant sequences in the dataset. This dinoflagellate-specific OTU sequence (Table 5), which was classified as Peridiniphycidae, was abundant in the SI April and May profiles, as well as in the WC May–June communities. Dinoflagellate-specific OTU65, chrysophycean OTU12 and OTU158 had a high similarity with sequences found in freshwater ecosystems (Table 5).
The Ciliophora-specific OTUs from the SI April profiles were classified as free-living ciliates of classes Oligohymenophorea (OTU13), Oligotrichia (OTU27), uncultured ciliates (OTU42) and Haptotia (OTU109). These OTUs were similar to environmental sequences from freshwater ecosystems, including those from cold habitats (Table 5). The LL WC April profiles were abundant in OTU22, OTU48 and OTU97, possessing similarity with Ciliophora-, chrysophyte- and fungal-specific sequences (Table 5).

3.6. Differentially Abundant OTUs of LL Communities in May–June

The species diversity of LL SI communities was decreased in comparison with that of WC (Table 2, Figure 3A,B and Figure 4). Out of 52 differentially abundant OTUs (relative abundance above 0.1%), 45 had increased abundance in the WC profiles, while seven OTUs were more abundant in SI (Table S5). A majority of OTUs, which were differentially abundant in the WC May profiles, possessed small absolute counts. Thus, the WC May profiles were characterized by increased species richness and diversity at the cost of appearance of a number of minor species. The exceptions from this pattern were the diatom-specific OTU1 and the chrysophyte-specific OTU21 (Table 5). Diatoms started massive development in late May in LL WC, while they were also abundant in the LL SI April profiles. Chrysophytes had their peak abundance in WC during May–June in both LL and LV, while they also were present in the SI May profiles.
Among the OTUs abundant in SI during May–June, there were sequences that were taxonomically assigned to dinoflagellates, chrysophytes, Cercosoa, ciliates and Fungi. The most abundant of these were Dinophyta-specific (OTU3 and OTU50) and chrysophyte-specific OTU404. As already noted, the LL SI April profiles had considerable similarity with the LL WC May–June ones. For a major number of microeukaryotic communities, similarity between the LL SI April and LL WC May–June groups was even more prominent, than that between the LL SI April and LL WC May (Figure 3B and Figure 4) groups. One may argue for a scenario where the communities of microeukaryotes, which were formed in the SI layer in April, migrate into the water column once the ice starts to degrade. Thus, the SI April communities might represent the primer stage for development of the summer microeukaryotic communities in Labynkyr.

3.7. Microeukaryotic Communities of LV and Their Comparison with LL Communities

There were four microeukaryotic community profiles in our dataset sampled from LV (Table 1). V1SI05 was collected in May 2015, while the three remaining samples (V1SI06, V1WC06 and V2WC06) were collected in June 2017. Although the LV communities were insufficiently sampled for in-depth statistical analysis in terms of their spatial and temporal variability, it seemed appropriate to use these results to compare the LV profiles between each other and with those of LL. The most prominent feature of the LV community profiles was the increased abundance of several OTUs, which were absent or very minor in the LL profiles (Figure 4, the upper part of heatmap). OTU1131 (Trebouxiophyceae), OTU444 (Chrysophyceae, genus Paraphysomonas), and OTU926 (Chlorophyceae) of this group were differentially abundant, as revealed by DESeq2 results (Table S5). These OTUs had sequences similar to those of chrysophytes and green algae from cold freshwater environments (Table 5).

3.8. Environmental Factors as Quantitative Explanatory Variables of Microeukaryotic Community Structure

We used the following continuous environmental variables to build a model of factors influencing the community structure: thickness of snow cover at sampling site; water temperature, pH, and electric conductivity (EC); concentrations of dissolved oxygen, NH4+, PO43−, NO3, NO2 ions, total mineral nitrogen (Nmin), total organic carbon (TOC); total microeukaryote abundance (TMA), and biomass (TMB) (Table 1). In tb-PCA, variability explained by community composition was 67.5%. Searching for the most important environmental predictors with forward selection, the variability explained was 49.2%, and variables selected were the snow thickness, mineral nitrogen, electric conductivity, and water temperature (Table S1, Figure 5A). The ordination pattern obtained with the constrained approach was similar to that with the unconstrained technique (Figure 3A,C). The results of variation partitioning between the model variables showed partial overlap between the snow thickness, mineral nitrogen concentration and water temperature; these factors seemed to have a strong conditional effect on each other (Figure 5B). Nevertheless, every model variable had an individual effect on the explained variation in the range between 5 and 10%: snow thickness, mineral nitrogen concentration, electric conductivity, and water temperature variables yielded 9, 8, 10 and 5%, respectively. These individual effects summed up to 32% of the total 49.2% of variation explained by the model. Thus, the remaining 17% could be thought of as a shared fraction of variation, where the effect of these four explanatory variables could not be disentangled.

3.9. Consistency of Results from the Point of View of Technical Replicates

The paired technical replicates, L4SI05/L4SI05-2 and L4WC06/L4WC06-2, could be treated as internal control for consistency of DNA library preparation, amplicon sequencing and data processing framework. The L4WC06/L4WC06-2 and L4SI05/L4SI05-2 were positioned very close on the ordination planes (Figure 3A,C and Figure 5A), had the smallest Bray–Curtis distances (Figure 3B,D), and revealed very similar abundance profiles in the heatmap of major OTUs (Figure 4). Such a robustness was observed in spite of rarefaction curves of L4SI05/L4SI05-2 profiles that were substantially different (Figure 1): the L4SI05 curve was more similar to that of V1SI05 than L4SI05-2. Additionally, the L4SI05-2 curve assumed lower sequencing depth and alpha-diversity indices (Table 2) than that of L4SI05. On the other hand, the rarefaction curves of L4WC06/L4WC06-2 were mainly similar to each other, as also confirmed by comparable alpha-diversity values. Thus, we argue that the replicates are consistent and expect that our results have no issues from a technical standpoint.

4. Discussion

4.1. Environmental Factors Influencing the Spatial Differences and Temporal Changes in LL Microeukaryotic Community Composition

Exploratory analysis suggested that during the ice cover period, LL SI and WC profiles are different (Figure 3 and Figure 5). This difference is relatively small in April, while WC and SI profiles are much less similar to each other in May–June. Electric conductivity and temperature were found to be explanatory variables that separated the SI and WC profiles (Table 4, Figure 5). As reverse stratification was expected during the ice cover period [44,45], electric conductivity and, consequently, mineralization in the sub-ice water layer were greater than in the water column [13,46]. This fact can be clearly demonstrated by the difference in EC between LL April samples, when the ice cover is stable (Figure 5A, quadrant IV). In mid May, the ice cover is starting to degrade because of the increase in average air temperature, thinning of the snow cover (Table 2), and the increase in light flux incoming into the under-ice environment. Degradation of the bottom surface of the ice results in demineralization of the sub-ice water layer. This moment is characterized by a sharp drop in the electric conductivity in sub-ice layer, while that of the water column remains constant (Figure 5A, quadrants I and III, Table 2).
Interestingly, an important factor influencing the community structure and revealed by the constrained ordination was the concentration of mineral nitrogen (Figure 5 and Figure S1(D10), Table S1). It was high in the April samples and decreased in May–June. This fact can be explained by the overall increase in ecosystem productivity linked with the increment in light flux [13]. The growth in primary production is associated with the increase in TMA, which was observed in WC May samples (Figure S1(D11)) and resulted in a decrease in free inorganic nutrients within the ecosystem.

4.2. Seasonal Dynamics of Microeukaryotic Community Structure

Considering the LV June samples, it should be noted that V1SI06/V1WC06 were collected after the process of active ice degradation began. It is not surprising that all three June LV profiles, including V1SI06, were grouped together (Figure 3 and Figure 5). Importantly, the LV May SI profile (V1SI05) and the corresponding LL May SI ones were very similar. Additionally, although good support was found for LL and LV ecosystems to be different in a set of hydrochemical parameters (see results of ANOVA), no clear differentiation of LL/LV microeukaryotic community profiles was observed. Taking this fact into consideration, one may suggest the changes in the microeukaryotic community composition in LL and LV were driven by analogous forces and followed the general compositional trend in May–June.
Summarizing the findings of the differential abundance analysis of LL April microeukaryotic communities, during the ice cover period, these were targeted to degrade the organic matter synthesized in the sub-ice layer. The sub-ice layer community was mainly responsible for primary production mediated by diatoms and dinoflagellates. Bacterivorous flagellates, dinoflagellates, ciliates, and Cercosoa are facultative heterotrophs that are able to compensate for a lack of nutrients and energy by feeding on bacterioplankton [47,48]. The LL WC microeukaryotic community profiles of May–June revealed an increased species richness. On the other hand, the SI May microeukaryotic communities were less diverse and mainly dominated by dinoflagellates. The peak of development of phytoplankton species, such as diatoms and green algae, suggests an increment in primary production at the end of May. Additionally, the primary production was facilitated by photosynthesis of mixotrophic species such as Chrysophyceae and Dinophyceae [35,36]. Considering the taxonomic affiliation of OTUs in LL and LV microeukaryotic communities on a global scale, we found a number of good OTU sequence matches with 18S rRNA barcodes from a wide variety of natural freshwater/brackish, arctic/boreal/temperate environments in Europe, Asia, North America, and even Antarctica (Table 5). This fact suggests that lacustrine ecosystems around the Earth are not isolated, and the exchange of species between geographically distant sites took place recently or is under way right now [49].

4.3. Peculiarities of the Specific Microeukaryotic High-Rank Phylotypes in LL and LV

HTS revealed a number of OTUs assigned to taxonomic groups of microeukaryotes that were never observed in LL and LV using microscopy. This allowed us to extend the knowledge on the species richness and diversity of the insufficiently explored phylotypes. While dinoflagellates, Chrysophyceae, Bacillariophyceae, Chlorophyta, and Prymnesiophyceae were observed in LL/LV by previous studies [10], the Ciliophora-specific OTUs, fungal OTUs assigned to Ascomycota and Basidiomycota, as well as many other eukaryotic phylotypes, were not previously reported in these lakes.
Dinoflagellate phylotypes were constantly present in the under-ice microeukaryotic communities of LL/LV. Their relative abundance decreased in June profiles and was relatively more constant in those of LL than of LV. The closest relatives of the most abundant dinoflagellate-specific OTUs were Peridinium aciculiferum/Scrippsiella hangoei (OTU3, which was constantly present in all profiles) [29], and Woloszynskia pascheri (OTU50, which was more abundant in SI profiles) [30]. These species are the typical members of spring phytoplankton blooms in under-ice freshwater ecosystems [39,40,50]. The massive bloom of dinoflagellates was reported in both marine [51] and freshwater environments [14,15,52], and likely accounts for a considerable proportion of the primary production in LL/LV ecosystems under high snow-cover thickness conditions. For instance, the sub-ice communities of Franklin Bay (Beaufort Sea, Canada) were shown to be dominated by dinoflagellates if the snow mantle thickness was considerable [53]. Furthermore, many dinoflagellate species are adapted to switch to a mixotrophic diet in low light conditions by feeding on bacterioplankton [49,50,54].
The composition of the chrysophyte-specific phylotypes was also changed from April to June. In April, the abundant chrysophyte-specific OTUs were observed in the SI layer; they proliferated in both SI and WC layers in May, while being abundant in WC in June. Chrysophytes are nanoflagellates, which are characterized by a small cell size and motile flagella. Most of the chrysophyte-specific OTUs were revealed to be similar with different species from a variety of freshwater environments (Table S3). The species richness of chrysophytes in LL and LV was studied previously by microscopy. Nineteen species of chrysophytes were classified as Chrysosphaerella, Paraphysomonas, Lepidochromonas, Spiniferomonas, and Mallomonas genera [9]. Silica-scaled chrysophytes are believed to have an advantage in ice-covered environments by covering a wide spectrum of trophic modes. The major factor influencing the change in abundance of phototrophic and mixotrophic chrysophytes is hypothesized to be the light influx into the under-ice environment [9]. Obviously, the available light influx depends on the snow cover thickness and the ice properties [13]. Notably, we reported the presence of chrysophyte-specific phylotypes (OTU12, OTU444, Table S3), which have high similarity with Chrysolepidomonas dendrolepidota from Lake Superior (USA) [31], Spumella-like isolate from Lake Hallstatt (Austria) [32], and Paraphysomonas foraminifera from the lake Mondsee (Austria) [32] (Table 5). These phylotypes were not previously detected in LL and LV by microscopy methods.
Ciliophora-specific phylotypes were not previously reported in LL and LV by microscopy. However, a variety of ciliate species were reported to inhabit cold freshwater environments. We detected several Intramacronucleata-specific OTUs, which were classified as Spirotrichea and had the closest sequence hits with SSU of free-living Oligohymenophorea [34], Cyrtostrombidium longisomum [35], and Limnostrombidium viride SW2012122001 [36] (Table 5). Ciliophora species are known to eventually bloom during the ice-cover period in Lake Baikal, as revealed by microscopy [55] and HTS [56]. Most likely, the seasonal dynamics of ciliate-specific phylotypes result from the change in the nutrient supply [57]. Ciliates are able to effectively feed on phytoplankton during the early spring blooms, and their diversity is enhanced when the abundance of phytoplankton species changes due to succession [58].
Bacillariophyta were one of the most abundant microeukaryotic phylotypes, which were present in both the SI and WC layers. According to previous study by optical microscopy, diatom species were shown to be quite specific and diverse in LL and LV [6]. Interestingly, more than half of the described species were detected in the under-ice environment. The abundance of diatoms was increased during the ice-cover period [10], with more than 80% of cold-adapted diatom species, such as phytoplankton and periphyton, continuing to grow in the open water conditions in summer [6]. Such a scenario of under-ice proliferation of diatom species is in line with recent reports of studies on phytoplankton in the long-term ice-covered freshwater lakes [14,59,60,61].
HTS allowed us to detect cercozoan phylotypes in the LL/LV microeukaryotic community profiles. Cercozoa are a diverse group of protists that have no common morphological features required for effective detection by microscopy methods. The Thecofilosea phylotype was abundant in all April microeukaryotic profiles. This taxon of amoeboid protists is known to contain mainly freshwater algivores, which are able to feed on unicellular algae and fungal cells [62]. Thus, Thecofilosea should be thought of as a predator taxon with a wide spectrum of dietary habits.
Fungal-specific phylotypes were not previously reported in LL and LV. At the same time, members of Basidiomycota/Cystobasidiomycetes were detected during all three sampling months, with the peak of abundance in LL in April and LV in June. Cystobasidiomycetes were reported to grow in psychrophylic environments, as demonstrated by the strains isolated from East Ongul Island, East Antarctica [63]. Cystobasidium tubakii and Cystobasidium ongulense, which are highly similar to OTU7, were reported to grow in below-zero temperatures on a vitamin-free medium. Such features were apparently evolved in these species during their adaptation to the cold and nutrient-depleted environments of subarctic lakes near the Pole of Cold in the northern hemisphere. Ascomycota-specific phylotypes mainly belonged to Thelebolaceae. These phylotypes mainly proliferated in the SI layer, which is also in line with their known features. Several species in the Thelebolaceae family were reported to express antifreeze and ice-binding proteins [64].
Cryptophyte-specific phylotypes were evenly distributed across all of the groups of the microeukaryotic community profiles. This can be explained by the fact that Katablepharidophyta, detected as one of the major OTUs with relative abundance above 0.01%, are characterized as a bacterivore species [65], and this phylotype does not strongly depend on dietary habits. Some Cryptophyte species were reported to be psychrophylic and are sensitive to light conditions. Thus, Chyptophytes are expected to be abundant during the under-ice period with a thick snow mantle. Cryptophytes were reported to proliferate during the whole winter season in high-mountain and subarctic/arctic freshwater environments. Furthermore, Cryptophyte species dominate the phytoplankton in the ice-covered freshwater Antarctic reservoirs, where they represent up to 70% of phytoplankton biomass [66].

5. Conclusions

In this study, the first inventory of the microeukaryotic communities of two large subarctic freshwater lakes in Yakutia was compiled by using HTS of 18S rRNA amplicons. As revealed by HTS, dinoflagellates, chrysophytes, Bacillariophyta, green algae, ciliates and Cercosoa were the major classes of phytoplankton and unicellular heterotrophs in the microeukaryotic communities. During the end of the ice cover period, from April to June, the complex communities of microeukaryotes developed in these lakes under harsh environmental conditions characterized by low temperature and decreased inorganic nutrient concentrations. We confirmed that the overall pattern of seasonal changes was shared between the lakes and likely was similar across different years. Importantly, the stratification of the sub-ice layer and the water column communities was constantly observed during the ice cover period, emphasizing the substantial role of temperature gradient and solar irradiance in the microeukaryotic community differentiation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d15030454/s1, Table S1: Impact of exploratory variables evaluated by redundancy analysis and “forward selection” approach; Table S2: Standardized abundance of top-ten classes of LL and LV microeukaryotic communities; Table S3: Taxonomic assignment and the closest GenBank matches of OTU sequences; Table S4: Manually curated taxonomic assignments of OTUs with standardized abundance above 10; Table S5: Differential abundance analysis of LL and LV microeukaryotic communities; Figure S1: One-way ANOVA of environmental variables with categorical variables “Lake”, “Month” and “Layer”; Figure S2: One-way ANOVA of alpha-diversity metrics with categorical variables “Lake”, “Month” and “Layer”.

Author Contributions

Conceptualization, Y.L., Y.G. and Y.Z.; methodology, Y.G.; software, Y.G.; validation, Y.G., Y.Z. and D.P.; formal analysis, Y.G.; investigation, Y.G., Y.Z., M.B., L.K. and D.P.; writing—original draft preparation, Y.G.; writing—review and editing, Y.G., Y.Z., Y.L. and M.B.; visualization, Y.G.; supervision, Y.L.; project administration, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the State Assignments of the Institute for Biological Problems of Cryolithozone, (FWSR-2021-0023, 121012190038-0; expeditions) and of the Limnological Institute (FWSR-2021-0008, 121032300186-9; analysis) of the Siberian Branch of the Russian Academy of Sciences.

Data Availability Statement

The sequence data were submitted to the Sequence Read Archive database (https://www.ncbi.nlm.nih.gov/sra, 19 December 2022) of the National Center for Biotechnology Information under Bio Project accession number PRJNA913696.

Acknowledgments

The authors acknowledge A.A. Dolzhenkov and I.E. Zhullyarov for organizing the expedition and technical support, and the divers A.S. Gubin, M.V. Astakhov, S.V. Bulochkin, and V.I. Chernykh for their assistance during the field expeditions. We gratefully acknowledge Irkutsk Supercomputer Center of Siberian Branch of the Russian Academy of Sciences for providing the access to HPC-cluster “Akademik V.M. Matrosov” and Shared Equipment Center for Integrated information and computing network of Irkutsk Research and Educational Complex for the data storage infrastructure. We would like to thank Reviewers for taking the time and effort necessary to review the manuscript. We sincerely appreciate all comments and suggestions of Reviewers, especially anonymous Reviewer #1, which helped us to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no competing interests.

References

  1. Debroas, D.; Domaizon, I.; Humbert, J.-F.; Jardillier, L.; Lepère, C.; Oudart, A.; Taïb, N. Overview of freshwater microbial eukaryotes diversity: A first analysis of publicly available metabarcoding data. FEMS Microbiol. Ecol. 2017, 93, fix023. [Google Scholar] [CrossRef] [Green Version]
  2. Schartau, A.K.; Mariash, H.L.; Christoffersen, K.S.; Bogan, D.; Dubovskaya, O.P.; Fefilova, E.B.; Hayden, B.; Ingvason, H.R.; Ivanova, E.A.; Kononova, O.N.; et al. First circumpolar assessment of Arctic freshwater phytoplankton and zooplankton diversity: Spatial patterns and environmental factors. Freshw. Biol. 2022, 67, 141–158. [Google Scholar] [CrossRef]
  3. Larsen, A.S.; O’Donnell, J.A.; Schmidt, J.H.; Kristenson, H.J.; Swanson, D.K. Physical and chemical characteristics of lakes across heterogeneous landscapes in arctic and subarctic Alaska. J. Geophys. Res. Biogeosci. 2017, 122, 989–1008. [Google Scholar] [CrossRef]
  4. Arzhakova, S.K.; Zhirkov, I.I.; Kusatov, K.I.; Androsov, I.M. A Quick-Reference Book on Rivers and Lakes of Yakutia; Bichik: Yakutsk, Russia, 2007; pp. 1–136. (In Russian) [Google Scholar]
  5. Bessudova, A.U.; Tomberg, I.V.; Firsova, A.D.; Kopyrina, L.I.; Likhoshway, Y.V. Silica-scaled chrysophytes in lakes Labynkyr and Vorota, of the Sakha (Yakutia) Republic, Russia. Nova Hedwigua 2019, 148, 35–48. [Google Scholar] [CrossRef]
  6. Kopyrina, L.; Firsova, A.; Rodionova, E.; Zakharova, Y.; Bashenkhaeva, M.; Usoltseva, M.; Likhoshway, Y. The insight into diatom diversity, ecology, and biogeography of an extreme cold ultraoligotrophic Lake Labynkyr at the Pole of Cold in the Northern Hemisphere. Extremophiles 2020, 24, 603–623. [Google Scholar] [CrossRef]
  7. Usoltseva, M.; Kopyrina, L.; Titova, L.; Morozov, A.; Firsova, A.; Zakharova, Y.; Bashenkhaeva, M.; Maslennikova, M.; Likhoshway, Y. Finding of a putative Lake Baikal endemic, Lindavia minuta, in distant lakes near the Arctic pole in Yakutia (Russia). Diatom Res. 2020, 35, 141–153. [Google Scholar] [CrossRef]
  8. Zakharova, Y.R.; Bedoshvili, Y.D.; Petrova, D.P.; Marchenkov, A.M.; Volokitina, N.A.; Bashenkhaeva, M.V.; Kopyrina, L.I.; Grachev, M.A.; Likhoshway, Y.V. Morphological description and molecular phylogeny of two diatom clones from the genus Ulnaria (Kützing) Compère isolated from an ultraoligotrophic lake at the Pole of Cold in the Northern Hemisphere, Republic of Sakha (Yakutia), Russia. Cryptogam. Algol. 2020, 41, 37–45. [Google Scholar] [CrossRef]
  9. Bessudova, A.; Firsova, A.; Bukin, Y.; Kopyrina, L.; Zakharova, Y.; Likhoshway, Y. Under-ice development of scaled chrysophytes with different trophic mode in two ultraoligotrophic lakes of Yakutia. Diversity 2023, 15, 326. [Google Scholar] [CrossRef]
  10. Zakharova, Y.; Bashenkhaeva, M.; Galachyants, Y.; Petrova, D.; Tomberg, I.; Marchenkov, A.; Kopyrina, L.; Likhoshway, Y. Variability of microbial communities in two long-term ice-covered freshwater lakes in the subarctic region of Yakutia, Russia. Microb. Ecol. 2022, 84, 958–973. [Google Scholar] [CrossRef]
  11. Hampton, S.E.; Galloway, A.W.E.; Powers, S.M.; Ozersky, T.; Woo, K.H.; Batt, R.D.; Labou, S.G.; O’Reilly, C.M.; Sharma, S.; Lottig, N.R.; et al. Ecology under lake ice. Ecol. Lett. 2017, 20, 98–101. [Google Scholar] [CrossRef]
  12. Obertegger, U.; Flaim, G.; Corradini, S.; Cerasino, L.; Zohary, T. Multi-annual comparisons of summer and under-ice phytoplankton communities of a mountain lake. Hydrobiologia 2022, 849, 4613–4635. [Google Scholar] [CrossRef]
  13. Jansen, J.; MacIntyre, S.; Barrett, D.C.; Chin, Y.-P.; Cortés, A.; Forrest, A.L.; Hrycik, A.R.; Martin, R.; McMeans, B.C.; Rautio, M.; et al. Winter limnology: How do hydrodynamics and biogeochemistry shape ecosystems under ice? J. Geophys. Res. Biogeosci. 2021, 126, e2020JG006237. [Google Scholar] [CrossRef]
  14. Bashenkhaeva, M.V.; Zakharova, Y.R.; Petrova, D.P.; Khanaev, I.V.; Galachyants, Y.P.; Likhoshway, Y.V. Sub-ice microalgal and bacterial communities in freshwater Lake Baikal, Russia. Microb. Ecol. 2015, 70, 751–765. [Google Scholar] [CrossRef]
  15. Bashenkhaeva, M.V.; Zakharova, Y.R.; Galachyants, Y.P.; Khanaev, I.V.; Likhoshway, Y.V. Bacterial communities during the period of massive under-ice dinoflagellate development in Lake Baikal. Microbiology 2017, 86, 524–532. [Google Scholar] [CrossRef]
  16. Bondarenko, N.A.; Timoshkin, O.A.; Ropstorf, P.; Melnik, N.G. The under-ice and bottom periods in the life cycle of Aulacoseira baicalensis (K. Meyer) Simonsen, a principal Lake Baikal alga. Hydrobiologia 2006, 568, 107–109. [Google Scholar] [CrossRef]
  17. Bondarenko, N.A.; Belykh, O.I.; Golobokova, L.P.; Artemyeva, O.V.; Logacheva, N.F.; Tikhonova, I.V.; Lipko, I.A.; Kostornova, T.Y.; Parfenova, V.V.; Khodzher, T.V.; et al. Stratified distribution of nutrients and extremophile biota within freshwater ice covering the surface of Lake Baikal. J. Microbiol. 2012, 50, 8–16. [Google Scholar] [CrossRef]
  18. Bashenkhaeva, M.V.; Galachyants, Y.P.; Khanaev, I.V.; Sakirko, M.V.; Petrova, D.P.; Likhoshway, Y.V.; Zakharova, Y.R. Comparative analysis of free-living and particle-associated bacterial communities of Lake Baikal during the ice-covered period. J. Great Lakes Res. 2020, 46, 508–518. [Google Scholar] [CrossRef]
  19. Rusch, D.B.; Halpern, A.L.; Sutton, G.; Heidelberg, K.B.; Williamson, S.; Yooseph, S.; Wu, D.; Eisen, J.A.; Hoffman, J.M.; Remington, K.; et al. The Sorcerer II Global Ocean Sampling Expedition: Northwest Atlantic through Eastern Tropical Pacific. PLoS Biol. 2007, 5, e77. [Google Scholar] [CrossRef] [PubMed]
  20. Bradley, I.M.; Pinto, A.J.; Guest, J.S. Design and evaluation of Illumina MiSeq-compatible, 18S rRNA gene-specific primers for improved characterization of mixed phototrophic communities. Appl. Environ. Microbiol. 2016, 82, 5878–5891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahe, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ 2016, 4, e2584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Pruesse, E.; Peplies, J.; Glöckner, F.O. SINA: Accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 2012, 28, 1823–1829. [Google Scholar] [CrossRef] [Green Version]
  25. Oksanen, J.; Simpson, G.L.; Blanchet, G.F.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’Hara, R.B.; Solymos, P.; Stevens, M.H.H.; Szoecs, E.; et al. Vegan: Community Ecology Package. R Package Version 2.5-6. 2022. Available online: https://CRAN.R-project.org/package=vegan (accessed on 24 November 2022).
  26. 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] [Green Version]
  27. Suzuki, R.; Terada, Y.; Shimodaira, H. Pvclust: Hierarchical Clustering with p-Values via Multiscale Bootstrap Resampling. R Package Version 2.2-0. 2019. Available online: https://CRAN.R-project.org/package=pvclust (accessed on 25 January 2023).
  28. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Logares, R.; Daugbjerg, N.; Boltovskoy, A.; Kremp, A.; Laybourn-Parry, J.; Rengefors, K. Recent evolutionary diversification of a protist lineage. Environ. Microbiol. 2008, 10, 1231–1243. [Google Scholar] [CrossRef]
  30. Lefèvre, E.; Bardot, C.; Noël, C.; Carrias, J.-F.; Viscogliosi, E.; Amblard, C.; Sime-Ngando, T. Unveiling fungal zooflagellates as members of freshwater picoeukaryotes: Evidence from a molecular diversity study in a deep meromictic lake. Environ. Microbiol. 2007, 9, 61–71. [Google Scholar] [CrossRef]
  31. Andersen, R.A.; de Peer, Y.V.; Potter, D.; Sexton, J.P.; Kawachi, M.; LaJeunesse, T. Phylogenetic Analysis of the SSU rRNA from Members of the Chrysophyceae. Protist 1999, 150, 71–84. [Google Scholar] [CrossRef]
  32. Boenigk, J.; Pfandl, K.; Stadler, P.; Chatzinotas, A. High diversity of the “Spumella-like” flagellates: An investigation based on the SSU rRNA gene sequences of isolates from habitats located in six different geographic regions. Environ. Microbiol. 2005, 7, 685–697. [Google Scholar] [CrossRef]
  33. Nakamura, K.; Sudo, M. Analysis of bacterial predatory protozoa inhabiting natural environments. J. Jpn. Soc. Civ. Eng. 2012, 68, III_31–III_40. [Google Scholar] [CrossRef] [Green Version]
  34. Monchy, S.; Sanciu, G.; Jobard, M.; Rasconi, S.; Gerphagnon, M.; Chabé, M.; Cian, A.; Meloni, D.; Niquil, N.; Christaki, U.; et al. Exploring and quantifying fungal diversity in freshwater lake ecosystems using rDNA cloning/sequencing and SSU tag pyrosequencing. Environ. Microbiol. 2011, 13, 1433–1453. [Google Scholar] [CrossRef]
  35. Tsai, S.-F.; Chen, W.-T.; Chiang, K.-P. Phylogenetic position of the genus Cyrtostrombidium, with a description of Cyrtostrombidium paralongisomum nov. spec. and a redescription of Cyrtostrombidium longisomum Lynn & Gilron, 1993 (Protozoa, Ciliophora) based on live observation, protargol impregnation, and 18S r DNA sequences. J. Eukaryot. Microbiol. 2014, 62, 239–248. [Google Scholar] [CrossRef]
  36. Gao, F.; Li, J.; Song, W.; Xu, D.; Warren, A.; Yi, Z.; Gao, S. Multi-gene-based phylogenetic analysis of oligotrich ciliates with emphasis on two dominant groups: Cyrtostrombidiids and strombidiids (Protozoa, Ciliophora). Mol. Phylogenet. Evol. 2016, 105, 241–250. [Google Scholar] [CrossRef] [Green Version]
  37. Nakai, R.; Abe, T.; Baba, T.; Imura, S.; Kagoshima, H.; Kanda, H.; Kohara, Y.; Koi, A.; Niki, H.; Yanagihara, K.; et al. Eukaryotic phylotypes in aquatic moss pillars inhabiting a freshwater lake in East Antarctica, based on 18S rRNA gene analysis. Polar Biol. 2012, 35, 1495–1504. [Google Scholar] [CrossRef] [Green Version]
  38. Motohashi, K.; Inaba, S.; Anzai, K.; Takamatsu, S.; Nakashima, C. Phylogenetic analyses of japanese species of Phyllosticta sensu stricto. Mycoscience 2009, 50, 291–302. [Google Scholar] [CrossRef]
  39. Charvet, S.; Vincent, W.F.; Lovejoy, C. Effects of light and prey availability on Arctic freshwater protist communities examined by high-throughput DNA and RNA sequencing. FEMS Microbiol. Ecol. 2014, 88, 550–564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Nicholls, K.H. Introduction to the biology and ecology of the freshwater cryophilic dinoflagellate Woloszynskia pascheri causing red ice. Hydrobiologia 2017, 784, 305–319. [Google Scholar] [CrossRef]
  41. Wolf, M.; Buchheim, M.; Hegewald, E.; Krienitz, L.; Hepperle, D. Phylogenetic position of the Sphaeropleaceae (Chlorophyta). Plant Syst. Evol. 2002, 230, 161–171. [Google Scholar] [CrossRef]
  42. Pröschold, T.; Darienko, T. Choricystis and Lewiniosphaera gen. nov. (Trebouxiophyceae, Chlorophyta), two different green algal endosymbionts in freshwater sponges. Symbiosis 2020, 82, 175–188. [Google Scholar] [CrossRef] [PubMed]
  43. Kulakova, N.V.; Kaskin, S.A.; Bukin, Y.S. The genetic diversity and phylogeny of green microalgae in the genus Choricystis (Trebouxiophyceae, Chlorophyta) in Lake Baikal. Limnology 2020, 21, 15–24. [Google Scholar] [CrossRef]
  44. Votintsev, K.K. Hydrochemistry of Lake Baikal. Proc. Baikal Limnol. St. 1961, 20, 1–312. [Google Scholar]
  45. Jewson, D.H.; Granin, N.G.; Zhdanov, A.A.; Gnatovsky, R.Y. Effect of snow depth on under-ice irradiance and growth of Aulacoseira baicalensis in Lake Baikal. Aquat. Ecol. 2009, 43, 673–679. [Google Scholar] [CrossRef]
  46. Granin, N.G.; Jewson, D.H.; Gnatovsky, R.Y.; Levin, L.A.; Zhdanov, A.A.; Gorbunova, L.A.; Tsekhanovsky, V.V.; Doroschenko, L.M.; Mogilev, N.Y. Turbulent mixing under ice and the growth of diatoms in Lake Baikal. SIL Proc. 1922–2010 2000, 275, 2812–2814. [Google Scholar] [CrossRef]
  47. Dagenais-Bellefeuille, S.; Morse, D. Putting the N in dinoflagellates. Front. Microbiol. 2013, 4, 369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Terrado, R.; Pasulka, A.; Lie, A.Y.; Orphan, V.J.; Heidelberg, K.B.; Caron, D.A. Autotrophic and heterotrophic acquisition of carbon and nitrogen by a mixotrophic chrysophyte established through stable isotope analysis. ISME J. 2017, 11, 2022–2034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Bessudova, A.; Bukin, Y.; Likhoshway, Y. Dispersal of Silica-Scaled Chrysophytes in Northern Water Bodies. Diversity 2021, 13, 284. [Google Scholar] [CrossRef]
  50. Charvet, S.; Vincent, W.F.; Comeau, A.; Lovejoy, C. Pyrosequencing analysis of the protist communities in a High Arctic meromictic lake: DNA preservation and change. Front. Microbiol. 2012, 3, 422. [Google Scholar] [CrossRef] [Green Version]
  51. Spilling, K. Dense sub-ice bloom of dinoflagellates in the Baltic Sea, potentially limited by high pH. J. Plankton Res. 2007, 29, 895–901. [Google Scholar] [CrossRef] [Green Version]
  52. Phillips, K.A.; Fawley, M.W. Winter phytoplankton blooms under ice associated with elevated oxygen levels. J. Phycol. 2002, 38, 1068–1073. [Google Scholar] [CrossRef]
  53. Różańska, M.; Gosselin, M.; Poulin, M.; Wiktor, J.M.; Michel, C. Influence of environmental factors on the development of bottom ice protist communities during the winter–spring transition. Mar. Ecol. Prog. Ser. 2009, 386, 43–59. [Google Scholar] [CrossRef] [Green Version]
  54. Kirchner, M.; Sahling, G.; Uhlig, G.; Gunkel, W.; Klings, K.-W. Does the red tide-forming dinoflagellate Noctiluca scintillans feed on bacteria? Sarsia 1996, 81, 45–55. [Google Scholar] [CrossRef]
  55. Obolkina, L.A. Planktonic ciliates of Lake Baikal. Hydrobiologia 2006, 568, 193–199. [Google Scholar] [CrossRef]
  56. Mikhailov, I.S.; Galachyants, Y.P.; Bukin, Y.S.; Petrova, D.P.; Bashenkhaeva, M.V.; Sakirko, M.V.; Blinov, V.V.; Titova, L.A.; Zakharova, Y.R.; Likhoshway, Y.V. Seasonal succession and coherence among bacteria and microeukaryotes in Lake Baikal. Microb. Ecol. 2022, 84, 404–422. [Google Scholar] [CrossRef]
  57. Posch, T.; Eugster, B.; Pomati, F.; Pernthaler, J.; Pitsch, G.; Eckert, E.M. Network of interactions between ciliates and phytoplankton during spring. Front. Microbiol. 2015, 6, 1289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Tirok, K.; Gaedke, U. Regulation of planktonic ciliate dynamics and functional composition during spring in Lake Constance. Aquat. Microb. Ecol. 2007, 49, 87–100. [Google Scholar] [CrossRef] [Green Version]
  59. Twiss, M.R.; McKay, R.M.L.; Bourbonniere, R.A.; Bullerjahn, G.S.; Carrick, H.J.; Smith, R.E.H.; Winter, J.G.; D’souza, N.A.; Furey., P.C.; Lashaway, A.R.; et al. Diatoms abound in ice-covered Lake Erie: An investigation of offshore winter limnology in Lake Erie over the period 2007 to 2010. J. Great Lakes Res. 2012, 38, 18–30. [Google Scholar] [CrossRef]
  60. Bertilsson, S.; Burgin, A.; Carey, C.C.; Fey, S.B.; Grossart, H.P.; Grubisic, L.M.; Jones, I.D.; Kirillin, G.; Lennon, J.T.; Shade, A.; et al. The under-ice microbiome of seasonally frozen lakes. Limnol. Oceanogr. 2013, 58, 1998–2012. [Google Scholar] [CrossRef]
  61. D’souza, N.A.; Kawarasaki, Y.; Ganyz, J.D., Jr.; Lee, R.E.; Beall, B.F.N.; Shtarkman, Y.M.; Koçer, Z.A.; Rogers, S.O.; Wildschutte, H.; Bullerjahn, G.S.; et al. Diatom assemblages promote ice formation in large lakes. ISME J. 2013, 7, 1632–1640. [Google Scholar] [CrossRef] [Green Version]
  62. Dumack, K.; Pundt, J.; Bonkowski, M. Food choice experiments indicate selective fungivorous predation in fisculla terrestris (Thecofilosea, Cercozoa). J. Eukaryot. Microbiol. 2019, 66, 525–527. [Google Scholar] [CrossRef]
  63. Tsuji, M.; Tsujimoto, M.; Imura, S. Cystobasidium tubakii and Cystobasidium ongulense, new basidiomycetous yeast species isolated from East Ongul Island, East Antarctica. Mycoscience 2017, 58, 103–110. [Google Scholar] [CrossRef]
  64. Batista, T.M.; Hilário, H.O.; de Brito, G.A.M.; Moreira, R.G.; Furtado, C.; de Menezes, G.C.A.; Rosa, C.A.; Rosa, L.H.; Franco, G.R. Whole-genome sequencing of the endemic Antarctic fungus Antarctomyces pellizariae reveals an ice-binding protein, a scarce set of secondary metabolites gene clusters and provides insights on Thelebolales phylogeny. Genomics 2020, 112, 2915–2921. [Google Scholar] [CrossRef] [PubMed]
  65. Grujcic, V.; Nuy, J.K.; Salcher, M.M.; Shabarova, T.; Kasalicky, V.; Boenigk, J.; Jensen, M.; Simek, K. Cryptophyta as major bacterivores in freshwater summer plankton. ISME J. 2018, 12, 1668–1681. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Henshaw, T.; Laybourn-Parry, J. The annual patterns of photosynthesis in two large, freshwater, ultra-oligotrophic Antarctic lakes. Polar Biol. 2002, 25, 744–752. [Google Scholar] [CrossRef]
Figure 1. Rarefaction curves of microeukaryotic communities. Curves are grouped by the month of sampling.
Figure 1. Rarefaction curves of microeukaryotic communities. Curves are grouped by the month of sampling.
Diversity 15 00454 g001
Figure 2. Hellinger-transformed relative abundances of the top-13 microeukaryotic classes in LL and LV communities. The bar stacks are arranged into six facets by two factor variables: “Layer” and “Month”. The profiles are horizontally sorted by sampling site, month and sampling layer. Taxa presented in the color legend are sorted by total abundance in decreasing order, except the last category “Other”. The order of Classes in each bar is the same as in the color legend.
Figure 2. Hellinger-transformed relative abundances of the top-13 microeukaryotic classes in LL and LV communities. The bar stacks are arranged into six facets by two factor variables: “Layer” and “Month”. The profiles are horizontally sorted by sampling site, month and sampling layer. Taxa presented in the color legend are sorted by total abundance in decreasing order, except the last category “Other”. The order of Classes in each bar is the same as in the color legend.
Diversity 15 00454 g002
Figure 3. Exploratory analysis of LL and LV microeukaryotic community profiles. Upper and bottom rows of the image were generated using different minimum OTU abundance thresholds: OTU relative abundance ≥0.1% (“major OTUs”) (A,B) and OTU standardized total abundance ≥50 (C,D). The former OTU count matrix contains 95 top-abundant OTUs with minimum OTU standardized total abundance ≥361. The latter count matrix consists of 253 OTUs. (A,C) tb-PCoA ordination biplots. Point shape: SI—circle, WC—square. Point color: LL—grey, LV—yellow. Blue arrows—linear regression of explanatory variables, showing the direction and range of their impact. (B,D) UPGMA-assisted community profile clustering dendrograms based on Bray–Curtis pairwise distance matrices. Red and green numbers indicate the “approximately unbiased” p-values and “bootstrap probability” values, respectively. Color annotation above the sample IDs encodes the spatial and temporal sampling parameters.
Figure 3. Exploratory analysis of LL and LV microeukaryotic community profiles. Upper and bottom rows of the image were generated using different minimum OTU abundance thresholds: OTU relative abundance ≥0.1% (“major OTUs”) (A,B) and OTU standardized total abundance ≥50 (C,D). The former OTU count matrix contains 95 top-abundant OTUs with minimum OTU standardized total abundance ≥361. The latter count matrix consists of 253 OTUs. (A,C) tb-PCoA ordination biplots. Point shape: SI—circle, WC—square. Point color: LL—grey, LV—yellow. Blue arrows—linear regression of explanatory variables, showing the direction and range of their impact. (B,D) UPGMA-assisted community profile clustering dendrograms based on Bray–Curtis pairwise distance matrices. Red and green numbers indicate the “approximately unbiased” p-values and “bootstrap probability” values, respectively. Color annotation above the sample IDs encodes the spatial and temporal sampling parameters.
Diversity 15 00454 g003
Figure 4. Heatmap of LL and LV microeukaryotic community profiles generated with a set of 95 top-abundant OTUs (see Figure 3A,B). Order of communities (columns) as in Figure 3B. UPGMA-assisted clustering of the Bray–Curtis pairwise distance matrices was used to produce the sample-wise and OTU-wise clustering trees. Color annotations above the heatmap describe the spatial and temporal categories of communities. Color annotation on the left denotes the OTU taxonomic affiliation on the “Class” level and is the same as in Figure 2.
Figure 4. Heatmap of LL and LV microeukaryotic community profiles generated with a set of 95 top-abundant OTUs (see Figure 3A,B). Order of communities (columns) as in Figure 3B. UPGMA-assisted clustering of the Bray–Curtis pairwise distance matrices was used to produce the sample-wise and OTU-wise clustering trees. Color annotations above the heatmap describe the spatial and temporal categories of communities. Color annotation on the left denotes the OTU taxonomic affiliation on the “Class” level and is the same as in Figure 2.
Diversity 15 00454 g004
Figure 5. Constrained ordination of count matrix with 95 top-abundant OTUs (see Figure 3A,B). (A) Ordination plot computed with transformation-based RDA. Exploratory variables were chosen by “forward selection” approach and are indicated by the blue vectors. Point shape: SI—circle, WC—square. Point color: LL—grey, LV—yellow. Gradients of electric conductivity and temperature are drawn by orange and green isolines. Gradients of the first two model variables—Snow and Nmin—are not shown, as the are nearly linear. (B) Venn diagram describing the variation partitioning.
Figure 5. Constrained ordination of count matrix with 95 top-abundant OTUs (see Figure 3A,B). (A) Ordination plot computed with transformation-based RDA. Exploratory variables were chosen by “forward selection” approach and are indicated by the blue vectors. Point shape: SI—circle, WC—square. Point color: LL—grey, LV—yellow. Gradients of electric conductivity and temperature are drawn by orange and green isolines. Gradients of the first two model variables—Snow and Nmin—are not shown, as the are nearly linear. (B) Venn diagram describing the variation partitioning.
Diversity 15 00454 g005
Table 1. Hydrophysical and hydrochemical characteristics of samples.
Table 1. Hydrophysical and hydrochemical characteristics of samples.
MonthSampleTemppHECDOPO43−NH4+NO2NO3NminTDSTOCTMATMBSnowIce
AprilL1SI040.48.284910.60.0010.0210.0020.350.133.193.2619.570.073091
L1WC041.27.8443.2210.60.0030.0590.0020.30.1229.292.5813.70.033091
L2SI041.27.7145.987.50.0010.0620.0020.340.1331.822.8824.10.063286
L2WC042.67.9339.497.50.0020.0580.0010.290.1126.672.2135.950.013286
L3SI040.47.5544.928.90.0010.0140.0020.330.0830.723.4520.950.0230109
L3WC041.97.9339.498.90.0020.0580.0010.290.1126.672.2128.850.0330109
L4SI040.47.4850.18.20.0040.0190.0010.420.1133.963.0440.40.043586
L4WC043.77.739.858.20.0010.0690.0020.280.1226.543.0439.550.133586
MayL1SI050.58.7838.418.50.0090.0130.0020.330.0825.331.4515.750.015111
L1WC052.69.3445.728.50.0030.0120.0020.340.0930.32.14222.230.085111
L3SI051.27.8123.349.80.0060.0070.0030.190.0517.161.9934.40.030.5110
L3WC053.49.2640.769.80.0020.0160.0010.20.0927.712.66313.670.030.5110
L4SI051.36.978.5179.600.0090.0020.170.056.561.552.560.060.586
V1SI050.48.9441.488.50.0080.0010.0030.10.0231.741.24192.70.031100
V1WC053.18.8251.218.50.0120.00400.080.0242.470.8679.550.081100
JuneL1WC063.66.830.898.30.0030.0170.0030.330.09-1.05173.30.025110
L3WC062.56.9737.717.30.0030.0120.0030.380.09-0.86127.10.03180
L4WC065.66.9841.078.40.0160.0120.0050.450.11-0.86121.30.050.586
V1SI060.47.1549.556.70.0140.0150.0030.050.02-0.8355.670.025110
V1WC063.27.2154.756.70.0230.0160.0060.060.02-1.8495.130.045110
V2WC062.27.133.88.30.0140.0160.0030.040.03-0.7584.910.02583
Column legend: Month—month of sampling. Sample—Sample id: the sampling site is encoded by the first two symbols, with the first letter denoting Labynkyr or Vorota and the digit denoting the sampling point; the fourth and fifth letters denote the sampling layer: sub-ice (SI) or water column (WC); month is encoded by the last two digits. Temp—water temperature, °C; EC—water electric conductivity, μS/cm2; DO—dissolved oxygen, mg/L; PO43−—concentration of phosphate anion, mg/L; NH4+—concentration of ammonium cation, mg/L; NO2—concentration of nitrite anion, mg/L; NO3—concentration of nitrate anion, mg/L; Nmin—total concentration of mineral nitrogen from ammonium, nitrate and nitrite, mg/L; TDS—total dissolved solids, mg/L; TOC—total organic carbon, mg/L; TMA—total microalgae abundance, ×103 cells/L; TMB—total microalgae biomass, g/m3; Snow—snow cover thickness at the moment of sampling, cm; Ice—ice thickness at the moment of sampling, cm.
Table 2. Summary NGS-data statistics and alpha-diversity indices of LL and LV microeukaryotic communities.
Table 2. Summary NGS-data statistics and alpha-diversity indices of LL and LV microeukaryotic communities.
SampleReadsRichnessACEShannonSimpsonInverse Simpson
L1SI0414,781155168.373.180.9010.11
L1WC0451,9849898.002.540.805.05
L2SI0410,12399102.693.210.9211.95
L2WC0445,9147373.002.740.856.67
L3SI0412,274124133.402.690.867.34
L3WC0411,225113114.463.050.899.26
L4SI0410,601123124.473.100.909.53
L4WC0410,2589393.703.300.9415.55
L1SI0520,9797577.400.990.311.45
L1WC0526,057191191.913.640.9314.74
L3SI0511,460106114.761.620.542.19
L3WC0514,548206212.883.480.9111.04
L4SI05 *35,709161177.542.250.815.15
L4SI05-2 *21,7349698.732.230.825.41
V1SI0517,205244257.792.900.825.64
V1WC05531-----
L1WC0617,685213224.493.310.909.74
L3WC0627,683337340.354.160.9521.76
L4WC06 *22,401326331.594.050.9417.57
L4WC06-2 *24,580153171.892.420.815.35
V1SI067032112134.982.270.753.99
V1WC0614,788133134.422.910.825.58
V2WC0615,535131135.603.090.9211.78
Column legend: Sample—sample id; Reads—number of reads after quality control; Richness—observed number of OTUs0.03; ACE—Abundance-based coverage estimator index; Shannon—Shannon index; Simpson—Simpson index; Inverse Simpson—inverse Simpson index; *—L4SI05/L4SI05-2 and L4WC06/L4WC06-2 samples are technical replicates.
Table 3. One-way ANOVA computed for environmental variables and alpha-diversity metrics.
Table 3. One-way ANOVA computed for environmental variables and alpha-diversity metrics.
Indep * Var #1Indep Var #2Dep VarANOVAKruskal–Wallis
ppadjpadj Signppadjpadj Sign
Environmental variables
Lake-PO42.64 × 10−41.23 × 10−3**7.56 × 10−33.53 × 10−2*
Lake-NO21.81 × 10−26.33 × 10−2.1.48 × 10−25.19 × 10−2.
Lake-NO36.52 × 10−65.49 × 10−5***2.44 × 10−31.71 × 10−2*
Lake-totN7.84 × 10−65.49 × 10−5***2.25 × 10−31.71 × 10−2*
Lake-TOC3.39 × 10−29.49 × 10−2.2.32 × 10−26.49 × 10−2.
Month-pH7.36 × 10−42.49 × 10−3**5.73 × 10−31.61 × 10−2*
Month-DO3.76 × 10−25.84 × 10−2.2.71 × 10−24.21 × 10−2*
Month-PO42.92 × 10−36.81 × 10−3**1.22 × 10−22.13 × 10−2*
Month-NH46.06 × 10−42.49 × 10−3**1.86 × 10−37.58 × 10−3**
Month-NO28.88 × 10−42.49 × 10−3**2.16 × 10−37.58 × 10−3**
Month-Nmin6.54 × 10−31.31 × 10−2*7.66 × 10−31.79 × 10−2*
Month-TOC7.00 × 10−64.90 × 10−5***6.49 × 10−44.54 × 10−3**
Month-TMA2.20 × 10−23.86 × 10−2*1.21 × 10−22.13 × 10−2*
Month-Snow4.15 × 10−155.82 × 10−14***6.03 × 10−44.54 × 10−3**
Layer-Temp2.81 × 10−53.93 × 10−4***2.69 × 10−43.77 × 10−3**
WCMonthpH8.98 × 10−82.52 × 10−6***1.34 × 10−29.87 × 10−2.
SIMonthPO42.89 × 10−27.35 × 10−2.1.66 × 10−12.74 × 10−1.
WCMonthNH41.87 × 10−72.53 × 10−6***2.52 × 10−29.87 × 10−2.
WCMonthNO21.64 × 10−24.72 × 10−2*1.89 × 10−29.87 × 10−2.
SIMonthNO31.54 × 10−24.72 × 10−2*4.50 × 10−29.87 × 10−2.
SIMonthNmin1.69 × 10−24.72 × 10−2*5.00 × 10−21.00 × 10−1.
SIMonthTOC2.52 × 10−41.17 × 10−3**3.57 × 10−29.87 × 10−2.
WCMonthTOC1.76 × 10−37.06 × 10−3**2.20 × 10−29.87 × 10−2.
WCMonthTMA1.41 × 10−47.90 × 10−4***1.36 × 10−29.87 × 10−2.
SIMonthSnow3.99 × 10−62.80 × 10−5***3.80 × 10−29.87 × 10−2.
WCMonthSnow2.72 × 10−72.53 × 10−6***2.39 × 10−29.87 × 10−2.
Alpha-diversity metrics **
Month-Richness3.95 × 10−29.87 × 10−2.1.54 × 10−23.84 × 10−2*
Month-ACE2.98 × 10−29.87 × 10−2.1.16 × 10−23.84 × 10−2*
Layer-Shannon3.16 × 10−31.58 × 10−2*6.86 × 10−33.43 × 10−2*
Layer-Simpson3.74 × 10−26.24 × 10−2.2.50 × 10−24.16 × 10−2*
Layer-Inverse Simpson2.00 × 10−24.99 × 10−2*2.50 × 10−24.16 × 10−2*
WCMonthRichness2.71 × 10−26.78 × 10−2.6.58 × 10−32.10 × 10−2*
WCMonthACE2.30 × 10−26.78 × 10−2.6.58 × 10−32.10 × 10−2*
SIMonthShannon6.84 × 10−33.42 × 10−2*1.05 × 10−22.10 × 10−2*
SIMonthInverse Simpson1.03 × 10−31.03 × 10−2*1.05 × 10−22.10 × 10−2*
Column legend: Indep Var #1/2—independent variables used to divide the dataset into subgroups; Dep Var—tested dependent variable; ANOVA/Kruskal–Wallis p-values; pp-value; padj—FDR-adjusted p-value; padj sign—the adjusted p-value significance code: *** ≤ 0.001; 0.001 ≤ ** ≤ 0.01; 0.01 ≤ * ≤ 0.05; 0.05 ≤ . ≤ 0.1; *—only rows with ANOVA padj ≤ 0.1 are shown, i.e., cases where the factor-wise distributions of dependent variable differ significantly; **—May and June groups were merged when the alpha-diversity metrics were tested by ANOVA.
Table 4. Two-way ANOVA with interaction of factors computed for environmental variables and alpha-diversity metrics.
Table 4. Two-way ANOVA with interaction of factors computed for environmental variables and alpha-diversity metrics.
Factor *Dep Varppadjpadj Sign
Environmental variables
LayerTemp1.83 × 10−29.61 × 10−2.
MonthEC1.47 × 10−28.84 × 10−2.
LayerNH49.94 × 10−41.39 × 10−2*
MonthNmin1.01 × 10−27.05 × 10−2.
MonthTOC1.53 × 10−53.22 × 10−4***
Month:LayerTOC8.00 × 10−36.72 × 10−2.
Month:LayerTMA7.64 × 10−36.72 × 10−2.
MonthSnow1.58 × 10−106.63 × 10−9***
Alpha-diversity metrics **
Month:LayerRichness3.61 × 10−28.30 × 10−2.
Month:LayerACE3.87 × 10−28.30 × 10−2.
MonthShannon2.16 × 10−31.62 × 10−2*
Month:LayerShannon8.96 × 10−41.34 × 10−2*
MonthSimpson8.41 × 10−34.20 × 10−2*
Month:LayerSimpson2.57 × 10−28.19 × 10−2.
MonthInverse Simpson4.54 × 10−28.51 × 10−2.
Month:LayerInverse Simpson2.73 × 10−28.19 × 10−2.
Column legend: Factor (Month, Layer or Month:Layer)—independent variable or their interaction; Dep Var—tested dependent variable; pp-value; padj—FDR-adjusted p-value; padj sign—the adjusted p-value significance code: *** ≤0.001; 0.001 ≤ ** ≤ 0.01; 0.01 ≤ * ≤ 0.05; 0.05 ≤ . ≤ 0.1; *—only rows with padj ≤ 0.1 are shown, i.e., cases where the two-factor interaction model has the significant impact to the dependent variable; **—May and June groups were merged when the alpha-diversity metrics were tested by ANOVA.
Table 5. Environmental sources of selected SSU sequences having high-similarity matches with OTUs from LL and LV.
Table 5. Environmental sources of selected SSU sequences having high-similarity matches with OTUs from LL and LV.
OUT
Number
Nearest MatchGenBank AccessionSequence Identity, %Source
OTU3Scrippsiella hangoei strain SHTV6
Peridinium aciculiferum strain PAER-2
EF417316
EF417314
99The Gulf of Finland, Baltic Sea
Lake Erken, Sweden [29]
OTU65Uncultured alveolateDQ24401999Lake Pavin, France [30]
OTU12Chrysolepidomonas dendrolepidota CCMP293Spumella-like isolate JBM19AF123297
FR865768
98Lake Superior, Keeweenaw Country, MI, USA [31]
Lake Hallstatt, Austria [32]
OTU158Bacterivorous protozoaAB74914997Hirose River, Japan [33]
OTU13Free-living ciliate
of class Oligohymenophorea
LR025746
HQ219368
98Lake Zurich, Switzerland
Lake Aydat, France [34]
OTU27Cyrtostrombidium longisomumKJ53458297Coastal waters of northeastern Taiwan [35]
OTU42Uncultured ciliate
Uncultured ciliate
Limnostrombidium viride SW2012122001
GU067975 LC165025
KU525754
99Esch-sur-Sure, Luxembourg
Biwa, JapanZhanjiang, Guangdong province, China [36]
OTU109Uncultured eukaryoteAB69550599Freshwater lake, East Antarctica [37]
OTU1Stephanodiscus sp.AB43059499Japan and South Korea [38]
OTU21Chrysophycean isolateAY08297094High-latitude Arctic lakes,
Ellesmere Island [39]
OTU50Woloszynskia pascheriEF05825396Freshwater environment [40]
OTU444Paraphysomonas foraminifera strain TPC2AY651096 99Lake Mondsee, Austria [32]
OTU926Ankyra lanceolata strain Hg 1998-5AF30276998Channel of Danube, Hungary [41]
OTU1131Choricystis sp. AS-29,
green alga symbiont of sponges
AY19597299Freshwater, free-living [42]
Lake Baikal, sponge symbiont [43]
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

Galachyants, Y.; Zakharova, Y.; Bashenkhaeva, M.; Petrova, D.; Kopyrina, L.; Likhoshway, Y. Microeukaryotic Communities of the Long-Term Ice-Covered Freshwater Lakes in the Subarctic Region of Yakutia, Russia. Diversity 2023, 15, 454. https://doi.org/10.3390/d15030454

AMA Style

Galachyants Y, Zakharova Y, Bashenkhaeva M, Petrova D, Kopyrina L, Likhoshway Y. Microeukaryotic Communities of the Long-Term Ice-Covered Freshwater Lakes in the Subarctic Region of Yakutia, Russia. Diversity. 2023; 15(3):454. https://doi.org/10.3390/d15030454

Chicago/Turabian Style

Galachyants, Yuri, Yulia Zakharova, Maria Bashenkhaeva, Darya Petrova, Liubov Kopyrina, and Yelena Likhoshway. 2023. "Microeukaryotic Communities of the Long-Term Ice-Covered Freshwater Lakes in the Subarctic Region of Yakutia, Russia" Diversity 15, no. 3: 454. https://doi.org/10.3390/d15030454

APA Style

Galachyants, Y., Zakharova, Y., Bashenkhaeva, M., Petrova, D., Kopyrina, L., & Likhoshway, Y. (2023). Microeukaryotic Communities of the Long-Term Ice-Covered Freshwater Lakes in the Subarctic Region of Yakutia, Russia. Diversity, 15(3), 454. https://doi.org/10.3390/d15030454

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