Next Article in Journal
Application of Freeze Concentration Technologies to Valorize Nutrient-Rich Effluents Generated from the Anaerobic Digestion of Agro-Industrial Wastes
Previous Article in Journal
Dynamics and Economics of Shallow Lakes: A Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity Variation of Silica-Scaled Chrysophytes Related to Differences in Physicochemical Variables in Estuaries of Rivers in an Arctic Watershed

1
Limnological Institute Siberian Branch of the Russian Academy of Sciences, 3 Ulan-Batorskaya, 664033 Irkutsk, Russia
2
Institute for Biological Problems of Cryolithozone, Siberian Branch of the Russian Academy of Science, 41 Lenina Avenue, 677980 Yakutsk, Russia
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(24), 13768; https://doi.org/10.3390/su132413768
Submission received: 10 November 2021 / Revised: 3 December 2021 / Accepted: 6 December 2021 / Published: 13 December 2021

Abstract

:
The present study examined the taxonomic diversity of silica-scaled chrysophytes in the estuaries of the Arctic watershed of Yakutia in the context of global climate change, as these aquatic organisms are highly sensitive to environmental changes. Previously, 41 species of silica-scaled chrysophytes were recorded in the waters of Yakutia. In the present study, we supplemented this list with 55 species. We observed a high species richness (82 taxa) of silica-scaled chrysophytes in the study area. Of these, eight species were recorded in the waters of Russia for the first time. At present, the study area represents the northernmost habitat for most species observed during the study. The diversity of silica-scaled chrysophytes in Arctic rivers is comparable to or even greater than the diversity of chrysophycean flora at the middle latitudes. Most of these chrysophytes are represented by polyzonal and ubiquitous species. During 2008–2010, we noted increased relative numbers of boreal species in northern waters compared with data obtained during the previous 30 years. Overall, the species richness of silica-scaled chrysophytes in the studied rivers increased with increase in water transparency, temperature, and magnesium ion concentration, but decreased with increase in surfactant concentration. The origin of samples from a certain river; distance of the sampling site from the estuary; water temperature, transparency, and colour; and concentration of carbon dioxide, magnesium ions, total iron, surfactants, and oil products affect the species composition of silica-scaled chrysophytes in the study area. Our findings are fundamental to evaluate the current status of Arctic aquatic microflora and its further monitoring in the context of anthropogenic and climatic impacts.

Graphical Abstract

1. Introduction

Representatives of the class Chrysophyceae have long been recognised as important indicators of environmental conditions [1,2,3]. They inhabit all climatic zones, including the tropics [4,5,6,7] and the Arctic [8,9,10]. However, their diversity and abundance are the highest at the middle latitudes [11,12,13], decreasing northward [14]. Few studies have explored the diversity of silica-scaled chrysophytes in northern polar waters (~66°33′ N). To date, 4 species in a small lake in Alaska [15], 17 species in the lakes of the Mackenzie Delta in Canada [16], 21 species in the Disko Island and Søndre Strømfjord region of West Greenland [17], 5 species in Central Norway [18], 32 species in Swedish Lapland [19,20], 6 species in Finnish Lapland [21], 75 species in the Bolshezemelskaya and Vorkuta tundra in Russia [22,23], 42 species in the Polar Urals of Russia [24,25], 44 species in the gulfs and bays of the Kara Sea in Russia [26,27], 25 species in the Khantay Reservoir in Russia [28], 26 species in the waters of the Taymyr Peninsula in Russia [8,9], and 8 species in the estuary of the Anabar River, Yakutia, Russia [10] have been recorded.
In the face of the current trends of global climate change [29,30] and rising atmospheric CO2 concentrations [31,32,33], silica-scaled chrysophytes at latitudes must be studied, considering that these organisms are highly sensitive to environmental changes. During the 20th century, the mean annual temperature of the Earth’s surface increased by 0.6 °C. In the Arctic, the temperature has risen to almost two times the earlier value and is expected to further increase by 4–7 °C in the present century [34]. The duration of open-water period, that is, the duration of the growing season [35], is an important factor that relies upon ambient temperature [35]. From 1980 to 2014, the duration of open water period in various water bodies in Asian Russia was extended by 4.63 to 11 days on average every decade [36], and this trend has been predicted to persist in the coming decades [34]. This prolongation of growing season favours the northward dispersal of southern species, leading to the expansion of their biogeographical ranges [37,38,39]. In a previous study [40], we analysed factors affecting the abundance and diversity of chrysophyceans in northern waters and hypothesised a scenario of their dispersal. Recent studies have frequently reported increased abundance of silica-scaled chrysophytes in boreal and Arctic lakes [29,30]. In their study in the Subarctic Region, Siver et al. [22] confirmed this trend for the first time based on the high diversity (75 species) of silica-scaled chrysophytes in lakes and ponds in the Bolshezemelskaya and Vorkuta tundra; the authors recorded some representatives of boreal flora, such as species of the genus Lepidochromonas (Lepidochromonas homolepis (Preisig & D.J. Hibberd) Kapustin & Guiry, Lepidochromonas quadrispina (Thomsen & Kristiansen) Kapustin & Guiry, Lepidochromonas subquadrangularis (Preisig & D.J. Hibberd) Kapustin & Guiry, and Lepidochromonas undulata (Preisig & D.J. Hibberd) Kapustin & Guiry) [22], at higher latitudes for the first time.
The water bodies we investigated in the present study are located north of the Arctic Circle (68°–72° parallel) in the Arctic zone. Samples collected from this region remain relatively understudied, particularly using electron microscopy. In previous studies using light microscopy, 32 species of silica-scaled chrysophytes were identified in the waters of Yakutia, including 21 species in the Lena River delta [41,42,43,44,45,46,47]. However, these records warrant a taxonomic revision given the limitations of light microscopy. Subsequent studies using electron microscopy revealed 41 species in the waters of lakes Labynkyr and Vorota [48], lake Toko and its surrounding area, and the estuary of Anabar [10].
Given water temperature is a factor determining the diversity of silica-scaled chrysophytes in northern reservoirs, we aimed to establish the biogeographic structures of communities of these organisms in the rivers of the Arctic region of Yakutia and compare these structures with those in other previously studied northern reservoirs. In this study, we tested the hypothesis of the introduction of species from temperate latitudes into Arctic reservoirs in the context of global climate change.

2. Site Description

2.1. Habitat Characteristics

The republic of Sakha (Yakutia) is the largest region in Russia (Figure 1), covering a total area of approximately 31,03,200 km2. Over 40% of Yakutia lies north of the Arctic Circle. The climate is extremely continental, with long winters and short summers [49].
The study sites are located north of the Polar Circle downstream the rivers, which flow in the Arctic geographic zone over three terrestrial ecosystems, including tundra, wooded tundra, and taiga. The southern boundary of the tundra in the region runs west-to-east along the right shores of the Bur River (left tributary of the Olenyok River), crossing the Lena River in the north of Kyusyur.
It runs parallel to the Indigirka River for up to 70 km and then turns southward toward the Kolyma River, crossing it near Nizhnekolymsk. The study area is situated in the permafrost zone.
Olenyok is one of the largest rivers, where the anthropogenic load is low. Amongst rivers in Asia, Olenyok is the fifteenth longest (2270 km) and the nineteenth largest by the area of watershed (219,000 km2) [50]. There are only two small settlements without industrial facilities, namely Olenyok and Taymylyr, on the banks of this large river. The river has never known navigation. Indigirka is a large river of the Arctic watershed in East Siberia; it is 1900 km long, and its drainage basin expands over 362,000 km2. There are almost no human activities along the river, and its ecosystems remain pristine. Yana is another large river of the Arctic watershed in East Siberia; it is 872 km long, and its drainage basin expands over 230,000 km2. The river begins at the confluence of the rivers Sartang and Dulgalakh. As the river flows into the Laptev Sea, it forms a huge river delta spanning 10,200 km2 [51]. Kolyma is a large river of the Arctic watershed in East Siberia; it is 2600 km long, and its drainage basin expands over 665,000 km2 [50]. Large fields of mineral resources are found in the drainage basins of these rivers; therefore, intense development of these areas is anticipated in near future. Mineral resources of Yana, Indigirka, and Kolyma are being developed since the second half of last century. The Kolyma Hydroelectric Station regulates the flow of the Kolyma River. The study sites are situated far from large scientific centres and remain hardly accessible. There is no developed route network in the study area, and the mode of transportation is mainly by air; hence, the field samples were collected by rafting.

2.2. Climate

The climate of the region is under a great influence of the Siberian High, which develops in winter at the centre of Asia; a large offset of this high covers the entirety of East Siberia. The climate is also strongly affected by frequent air intrusions from the Arctic Ocean, which are poor in aqueous vapour in summer. According to the data on WorldClim (www.worldclim.org), the mean annual ambient temperature at the sampling sites varies between −12.4 and −14.8 °C, with the maximum temperature recorded in summer (14.9–18.6 °C); mean annual precipitation is 160–256 mm. The mean duration of open water, which determines the length of the growing season, is 115, 118, 122, and 135 days in Olenyok, Yana, Indigirka, and Kolyma, respectively [52,53]. Based on our observations in summer when the water flow was low, the median water temperature of the rivers was 14.8 °C.

3. Materials and Methods

3.1. Field and Laboratory Methods

Samples were collected on 10–14 August 2008, 26 June–2 July 2009, 17–18 July 2009, and 5–7 August 2010, either from the littoral zone or from the midstream in the surface water layer (0–0.3 m) of Olenyok, Yana, Indigirka, and Kolyma, respectively, (see Figure 1).
Phytoplankton samples (n = 33) were collected using the Apstein net [54,55] with a SEFAR NITEX filter cloth (30 μm mesh). The collected material was fixed in 4% formaldehyde solution at the 1:10 ratio.
Transmission electron microscopy was used for the identification of chrysophytes based on the structure of their scales. Samples were transferred to 1.5 mL Eppendorf tubes and centrifuged using the MiniSpin centrifuge (Eppendorf, Germany) at 13,400 rpm for 10 min. The supernatant was removed using a pipette, and deionised water was added to the precipitate. The procedure was performed three times. Next, 30% H2O2 was added to the precipitate, and the sample was heated in a thermostat at 80 °C for 5 h. Thereafter, the precipitate was washed with deionised water and centrifuged. The procedure was performed five times. The washed samples were pipetted onto a grid (diameter, 3 mm) covered with Formvar support, dried at room temperature, and examined under a transmission electron microscope (LEO 906E; Carl Zeiss, Germany).
Chemical analyses of 25 water samples were performed following standard methods [56,57]. Gaseous components of water (O2 and CO2) and specific physical parameters were measured in situ, and the remaining chemical analyses were performed in vitro. A Secchi disc was used to measure water transparency. Water colour was determined using a photometric method with the SF-26 tool. pH was measured using a potentiometric method with the Multitest IPL-101 device. A gravimetric method was used to determine the suspended solid content [57]. Oxygen and carbon dioxide concentrations were measured using a titration method with iodimetric determination and phenolphthalein, respectively. Oxygen saturation was determined using a calculation method, and the hardness of water was measured using a titration method. Water salinity was calculated as the sum of ions using the following methods: turbidimetry for sulphate anions; flame spectrophotometry for potassium and sodium cations; mercurimetry for chloride ions; and titration for calcium, magnesium, and bicarbonate ions. A photometric method with the SF-26 tool was applied to determine nutrient concentrations. Nessler’s reagent, Griess reagent, salicylic acid, ammonium molybdate, and sulfosalicylic acid were used for the measurement of ammonium ion, nitrite ion, nitrate ion, phosphate and silicon ions, and total iron, respectively. A combined reagent composed of ammonium molybdate and ascorbic acid was used to determine total phosphorus content. Chemical oxygen demand (COD) was measured with the Fluorat-02 device using a luminescent method. A titration method with iodimetric determination was used to measure biological oxygen demand (BOD). A fluorimetric method with the Fluorat-02 device was used to measure the levels of oil products, phenols, and anionic surfactants. The maximum allowable concentrations (MACs) in the waters of commercial fishing [58] are used in this article.

3.2. Geographical Distribution

To determine eventual changes in the context of global climate change, we analysed the geographical distribution of silica-scaled chrysophytes in waters located north of the Arctic Circle over the long-term. Considering the scarcity of such investigations, with no studies conducted in some regions since 1981, we undertook preliminary revision of species composition to match the current taxonomy of silica-scaled chrysophytes (Table S1). We ascribed the taxa of silica-scaled chrysophytes to the formerly proposed latitudinal and longitudinal groups of geographical distribution [59,60] according to their preferred distribution in a given geographic zone and their geographic range over the continents. The following types of geographic distribution were included within the latitudinal group: P: polyzonal (species found in all climatic zones); A-Bor: Arctic-boreal (species found in the Northern Hemisphere in the temperate and/or Arctic zones); Bor: boreal (species only found in the boreal zone); and Arc: Arctic (species only found in the Arctic zone). The following types of geographic distribution were included within the longitudinal group: K: cosmopolitan (species found on all six continents); W: widespread (species absent from one or two continents); R: species with scattered distribution (species rarely found across different latitudes); BP: bipolar species (species occurring at the middle latitudes of the Northern and Southern Hemispheres); and End: endemic (species restricted to a certain geographical range in a continent). We used the tag ‘unknown’ i.e., with unknown geographical characteristic) for new species, bipolar species, or species identified only to the genus level. Since some species from the literature we analysed were no longer valid, species composition in Table S1 is presented following the taxonomic revision. Species complexes within the groups Paraphysomonas vestita (Stokes) De Saedeeler and Paraphysomonas imperforata Lucas as well as Synura petersenii f. kufferathii Petersen & Hansen, which could not be synonymised, are simply designated as species.

3.3. Statistical Analysis

Only those samples for which all physicochemical parameters were known, and which harboured at least one chrysophycean species, were included in the statistical analysis (n = 18).
PERMANOVA [61] in the R package ‘Vegan’ [62] was used for analysing the effects of abiotic physical and chemical parameters on the species composition of chrysophytes. Species composition data in samples were represented as a binary table, with rows and columns describing samples and species, respectively. The presence of a species in a sample was coded as 1 and its absence as 0. The Jaccard distance matrix was calculated for PERMANOVA based on the table of species composition. p values were calculated based on 10,000 permutations. The effect of a factor was assumed true when the p value was ≤0.05. Furthermore, variance coefficients (R2) were obtained using PERMANOVA, which represented the extent of variation in dependent variables (differences in species composition) driven by the independent variables.
Differences in species composition amongst samples based on the Jaccard distances were visualised using non-metric multidimensional scaling (NMDS). Gradient vectors of abiotic physical and chemical parameters that truly affected (according to PERMANOVA) the species composition of samples were plotted on an NMDS biplot. Additionally, a gradient vector of the species richness of samples was plotted on the NMDS biplot. All analyses were performed in R using the ‘Vegan’ [62] package following the tutorial [63].

3.4. Data Availability

All raw data, including the binary table of species per sample and physicochemical indicators of water per sample, are provided in Table S2.

4. Results

4.1. Diversity and Geographic Distribution of Chrysophytes

In total, 82 species belonging to six genera, namely Chrysosphaerella (n = 5), Paraphysomonas (n = 12), Lepidochromonas (n = 10), Spiniferomonas (n = 10), Mallomonas (n = 20), and Synura (n = 25) (Table 1, Figures S1–S7), were identified in the estuaries of Olenyok, Yana, Indigirka, and Kolyma. Of these, eight species, including Lepidochromonas eiffelii, Lepidochromonas elegantissima, Lepidochromonas subquadrangularis, Mallomonas lychenensis, Synura cf. americana, Synura cf. cornuta, Synura soroconopea, and Synura cf. vinlandica, were described in the waters of Russia for the first time. No silica-scaled chrysophytes were detected in samples collected from stations 10 (estuary of Yana), 13, 15–19, 22–25 (estuaries of Indigirka), 28, 29, and 33 (estuaries of Kolyma). Species composition differed amongst the four rivers.
The highest diversity of chrysophytes was recorded in Olenyok, with 75 species belonging to Chrysosphaerella (n = 5), Paraphysomonas (n = 12), Lepidochromonas (n = 10), Spiniferomonas (n = 10), Mallomonas (n = 16) and Synura (22) (Table 2). Of these, the cosmopolitan species Chrysosphaerella coronacircumspina, Paraphysomonas gladiata, Spiniferomonas bourrellyi, Spiniferomonas serrata, Spiniferomonas trioralis, Mallomonas akrokomos, Mallomonas crassisquama and Synura echinulate were the most abundant and most frequently detected at the study site. Station 1 differed from the other stations in the estuary in terms of hydrochemical parameters. At station 1, total ion content was 270.14 mg·L−1, and total phosphorus content was 186 μg·L−1 (Table 2). However, the species richness at station 1 (14 species) was lower than that at the other stations situated downstream. Further downstream in the north toward station 5, species composition became more diverse and species richness was the maximum (52 species). Amongst all stations in Olenyok, stations 4 and 5, located ~3 km apart, harboured the most diverse species composition. The observed differences in the species composition and hydrochemical parameters across stations of the same river can be explained by the effects of the Bur River, which flows close to station 4. At station 4, total ion content was 67.46 mg·L−1 and pH was 6.78; meanwhile, at station 5, these values were, respectively, 180.67 mg·L−1 and 7.52. Moreover, six species, namely Paraphysomonas acuminata, Spiniferomonas minuta, Mallomonas lychenensis, Mallomonas papillosa, Mallomonas cf. pumilio, and Synura sp. 3, were detected only at station 4. Meanwhile, 36 species were detected only at station 5. Furthermore, station 5 harboured high richness of species belonging to the genera Paraphysomonas and Lepidochromonas, which are atypical in the waters at high latitudes. Specifically, Paraphysomonas cf. punctata ssp. simplicior, Paraphysomonas vulgaris, Lepidochromonas eiffelii, and Lepidochromonas elegantissima, which are typical of the waters at middle latitudes, were identified in the Arctic Region for the first time. Water temperature at station 5 (16.8 °C) is favourable for the development of boreal species. Closer to the estuary, water temperature falls to 15 °C, drastically reducing species richness (21 species).
In the Yana River, 40 species of the genera Mallomonas (n = 11), Synura (n = 19), Chrysosphaerella (n = 3), Paraphysomonas (n = 3), and Spiniferomonas (n = 4) were identified (Table 1). Hydrochemical parameters were relatively stable across the sampling stations of this river. Only water temperature dropped from 14.6 °C at the southernmost station 8 to 14 °C at the most downstream sampling station 12. In contrast, the species diversity of silica-scaled chrysophytes increased from the south (6 species) to north (31 species). This trend may be attributed to the progression of vegetation peak, which typically occurs at high latitudes between mid- and late summer. Paraphysomonas cf. punctata ssp. simplicior, Spiniferomonas cornuta, Mallomonas actinoloma, Synura praefracta, Synura nygaardii, and Synura leptorhabda showed scattered distribution (including rare species). Moreover, Mallomonas parvula, Synura leptorhabda, and Synura sp. 1 were detected only in Yana.
The species composition in the Indigirka River was the least diverse, with only seven species of Mallomonas (2) and Synura (5). Interestingly, no species of the genus Paraphysomonas and Lepidochromonas spp. were found in the river (Table 1). The widest ranges of some hydrochemical parameters were recorded in this habitat (Table 2). In this river, water temperature and pH values are favourable for the growth of chrysophytes. However, low transparency and high concentrations of suspended solids limit the development of phytoplankton, including silica-scaled chrysophytes, at the majority of the stations. High turbidity might be a greater problem for non-photosynthetic species, such as Paraphysomonas spp. In the present study, silica-scaled chrysophytes were detected at only 3 of the 13 stations. All species detected in this river, except Synura praefracta, are cosmopolitan and ubiquitous.
In the Kolyma River, 31 species belonging to the genera Chrysosphaerella (n = 2), Paraphysomonas (n = 2), Spiniferomonas (n = 2), Mallomonas (n = 9), and Synura (n = 16) (Table 1) were detected. Of these, Mallomonas caudata, Synura spinosa, and Synura spinosa f. longispina were the most frequently detected species at the study sites. Similarly, to the trends in Indigirka, the same hydrochemical parameters showed wide ranges in Kolyma. For instance, water temperature varied between 12 °C and 16.7 °C, and no chrysophyceans were found in the regions with the lowest temperatures (12–14.4 °C). Similarly, total phosphorus content also showed a wide range (40 to 380 μg·L−1), and no chrysophyceans were found in regions with the lowest or highest concentration. The distribution of silica-scaled chrysophytes from the southernmost sampling site to the northernmost one was uneven. Specifically, 16 species were identified at the southernmost sampling site, which disappeared from stations 28 and 29 with a drop in temperature and a rise in total phosphorus content. Farther northward, where the hydrochemical parameters remained stable, the species richness of silica-scaled chrysophytes increased, reaching the maximum for this river (24 species at station 32). Only three species, namely Paraphysomonas vulgaris, Synura nygaardii, and Synura praefracta, showed scattered distribution, and the rest were cosmopolitan and ubiquitous.
Paraphysomonasuniformis ssp. hemiradia, Mallomonas acaroides, Mallomonas alpina, Mallomonas crassisquama, Mallomonas tonsurata, Synura borealis, Synura petersenii, Synura praefracta, Synura cf. americana, and Synura conopea were the most frequently detected species. All these species, except Synura praefracta, are cosmopolitan and ubiquitous. The prevalent taxa in estuaries can be divided into two groups based on their response to pH gradient [3]: pH-indifferent species, such as Mallomonas crassisquama and Synura petersenii, which occur in waters with a pH of ~7, and alkaliphilic species, such as Mallomonas acaroides, Mallomonas alpina, Mallomonas crassisquama, and Mallomonas tonsurata, which prefer waters with a high pH (>7).

4.2. Addition of Autecology of Rare Species

Notably, some rare species were observed in the estuaries of the studied rivers, and majority of them were detected at station 5 in Olenyok. The eight species described for the first time in the waters of Russia were found only at this station. We classified the species as rare if they had been described a long time ago in few studies, or if they had been described only recently because of their poorly known autecology, as in the case of Synura cf. vinlandica, Synura cf. cornuta, and Synura cf. soroconopea.
Synura cf. vinlandica (Figure S6o–q) was recently described based on molecular and morphological studies of isolated colonies in samples collected from Newfoundland, Canada [64]. We classified our samples as Synura cf. vinlandica based on the morphological features of scales. In the present study, we detected this species at station 7 in Olenyok (pH = 7.80, T = 15.0 °C) (Table 2).
Synura cf. cornuta (Figure S6m,n) was recently described based on molecular and morphological studies of isolated colonies in samples collected from Newfoundland, Canada [64]. We classified our samples as Synura cf. cornuta based on the most distinct morphological feature of scales possessing a ‘narrow and prominent tip, resembling a horn’ [64]. In the present study, we detected this species at stations 1 and 3 in Olenyok River (pH = 7.89–7.93, T = 16–16.4 °C) (Table 2).
Synura cf. soroconopea (Figure S7a,b) Jo, Shin, Kim & Siver was described based on molecular and morphological studies of isolated colonies in samples collected from small ponds in South Korea [65]. In the present study, we frequently detected this species at different stations in Olenyok, Yana, and Kolyma (pH = 6.69–7.93, T = 14–17.4 °C) (Table 2).
Chrysosphaerella rotundata (Figure S1f) was described from a water body in Finland [66]. Before its description, chrysophytes with rounded scales similar to Chrysosphaerella rotundata were found in Germany (classified by the author as Chrysosphaerella brevispina in [67]: Figure 76). The species has also been recorded in the Selenga Delta of Russia [13]. It has been reported to inhabit waters with temperatures between 5 °C [66] and 19 °C [67], pH between 7 [66] and 7.85 [13], and conductivity between 40 [67] and 340 μS·cm−1 [67]. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
Paraphysomonas circumvallata (Figure S1q) was described from a water body in Danish [68]. It has also been recorded in Greece [69], North America [70,71], Germany [72], New Zealand [73], and the Pojo Bay of Finland [74]. It was recorded for the last time in the Gulf of Finland, Baltic Sea [75]. However, only one previous study has reported the ideal hydrochemical parameters for the growth of this species in a lake in Germany (pH = 7.8, conductivity 136 μS·cm−1) [72]. In the present study, this species was detected at stations 3, 4, and 5 in Olenyok (pH = 6.78–7.93, T = 16.4–17.4 °C) (Table 2).
Lepidochromonas eiffelii (Figure S2f,g) was described from a water body in Danish [68]. It has also been recorded in England [76], Chile [77], North America [78], Germany [72], Hungary [79], and Austria [80], China [81]. It has been reported to inhabit waters with temperatures between 0 °C [68] and 20 °C [81], pH between 6.2 [72] and 8.7 [81], and conductivity between 102 [80] and 490 μS·cm−1 [72]. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
Lepidochromonaselegantissima (Figure S2h) was described from a water body in North America [70]. It has also been recorded in North America [71], Denmark [82], and Baltic Sea ice stations in Finland [83]. It has been reported to inhabit waters with temperatures between 0 °C [83] and 3 °C [82], pH between 5.6 [82] and 7.7 [83], and conductivity of 70 μS·cm−1 [82]. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
Lepidochromonas quadrispina (Figure S2e) was described from a water body in Danish [68]. It has also been recorded in Greece [69], Finland [11,84,85], and the Vorkuta tundra of Russia [22]. It has been reported to inhabit waters with temperatures between 10.7 °C and 13 °C [22], pH between 6.7 and 6.8 [84], and conductivity between 23 [22] and 47 μS·cm−1 [84]. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
Lepidochromonasundulata (Figure S2k) was described from a water body in England [76]. It has also been recorded in Finland [85], Hungary [79], and the Vorkuta tundra of Russia [22]. It has been reported to inhabit waters with temperatures between 10.5 °C ([79]) and 17 °C [22], pH between 6.2 and 6.8 [79], and conductivity between 24 [22] and 240 μS·cm−1 [79]. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
Spiniferomonas conica (Figure S3j) was described from a water body in Japan [86]. It has also been recorded in North America [87], Greenland [88], and the Lake Labynkyr in Russia [48]. It has been reported to inhabit waters with temperatures between 10 °C [88] and 15.8 °C [48], pH between 5.9 [88] and 7.63 [48], and conductivity between 11 [88] and 38 μS·cm−1 [48]. In the present study, this species was detected at station 2 in Olenyok (T = 16.4 °C) (Table 2).
Spiniferomonas minuta (Figure S3k) was described from a water body in North America [89,90]. It has also been recorded in North America (as Spiniferomonas takahashii in Nicholls, [87]: Figure 38; [89]), Finland [84], the Gulf of Bothnia, Baltic Sea [91], and small lakes of the Lower Yenisei Basin in Russia [27]. It has been reported to inhabit waters with temperatures between 4.5 °C [27] and 13.1 °C [84], pH between 5.6 [91], and 8.4 [27], conductivity between 20 [90] and 150 μS·cm−1 [27]. In the present study, this species was detected at station 4 in Olenyok (pH = 6.78, T = 17.4 °C) (Table 2).
Spiniferomonas takahashii (Figure S3l) was described from a water body in North America [87]. It has also been recorded in North America [90,92], the Gulf of Bothnia, Baltic Sea [91], and the Lake Baikal in Russia [93]. It has been reported to inhabit waters with temperatures between 8 °C and 20 °C [90], pH between 5.8 [91] and 8 [93], and conductivity between 34 [90] and 110 μS·cm−1 [93]. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
Synurapetersenii f. taymyrensis (Figure S7f) was described from a water body in the Taymyr Peninsula of Russia [9]. In the present study, Synura petersenii f. taymyrensis was detected at station 4 in Olenyok (pH = 6.78, T = 17.4 °C) (Table 2), representing the first record of this species since its description.

4.3. Undetermined Species

Eleven species showing scale morphology different from all known species were found in areas with high diversity of silica-scaled chrysophytes. These are likely to be new species, and their scales and spines are described below.
We found Chrysosphaerella sp. (Figure S1c,d) bearing two types of scales, including round (2.3–2.9 μm diameter) or oval (2.4–2.6 μm long and 1.4–1.7 μm wide) ones. Both types of scales were ornamented with series of elongated, round ridges, forming a serrated pattern. The considerably large, round ridges differentiate these scales from structures in other known species. In some oval scales, a secondary pattern was observed under the ridges (Figure S1d). In the present study, this species was detected at stations 4 and 5 in Olenyok (pH = 6.78–7.52, T = 16.8–17.4 °C) (Table 2).
We found Paraphysomonas sp. 1 (Figure S1h) bearing a round to oval, rimless baseplate (diameter, 0.9–1.1 μm). The plate bore 2.4–2.8-μm-long spines, tapering in a pinched tip. The proximal part of the spine (0.5 μm from the baseplate) was darker in colour, suggesting that this part was denser than the rest of the spine. In the present study, this species was detected at stations 4, 5, and 6 in Olenyok (pH = 6.78–7.51, T = 16.8–17.4 °C) (Table 2).
We found Paraphysomonas sp. 2 (Figure S1k) bearing a round baseplate (diameter, 0.6 μm) with a highly incurved border (width 0.04μm). The margins of the baseplate possessed small perforations. A spine (length 0.9 μm) protruded from the centre of the round baseplate, and it was rounded at the tip. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
We found Paraphysomonas sp. 3 (Figure S2a–c) bearing elliptical, deeply undulate plate-scales (0.9–1.5 × 1.3–1.9 μm) with a deep crater-like depression (crater-scales). The rim was 0.2–0.5 μm thick, with a ring of 70–74 small pegs along its inner edge on the distal side. Five to seven large apertures (diameter 0.08–0.2 μm) were located in the crater-like depression. In the present study, this species was detected at stations 4 and 5 in Olenyok (pH = 6.78–7.52, T = 16.8–17.4 °C) (Table 2).
We found Paraphysomonas sp. 4 (Figure S2d) bearing an elliptical baseplate (1.4–1.5 × 1.6–1.8 μm). The inner edge of the baseplate was ornamented with five rows of 62–85 min marginal pegs on the distal side. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
We found Lepidochromonas sp. 1 (Figure S2l) bearing perforated, elliptical plate-scales (1 × 1.2 μm). The plate-scale perforations varied in size and were randomly distributed. In the central part of the scale, there were four large apertures (diameter, 0.2 μm). In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
We found Lepidochromonas sp. 2 (Figure S2n) bearing crown-like scales with an elliptical base, which was wider on one side of the scale (0.6 μm) but slightly narrower on the other (0.5 μm). The base of the crown-like scales possessed six long, perpendicular shafts, distally connected with the arches. The arch structures were perforated. In the present study, this species was detected at station 5 in Olenyok (pH = 7.52, T = 16.8 °C) (Table 2).
We found Synura sp. 1 (Figure S5o) bearing oval body scales (2.8 × 2.3 μm) with an obtuse (flat) tip. Spine was unknown. There was no upturned edge around the periphery of the scale. The basal plate showed regular pore ornamentation. The vermiform ornamentation area on the distal portion was absent. In the distal flange area, 15 ribs (0.27 μm long) were present, bearing large pores between them. In the present study, this species was detected at station 12 in Yana River (pH = 7.13, T = 14 °C) (Table 2).
We found Synura sp. 2 (Figure S7c,d) bearing body scales (3.8–4.3 μm long and 2.3–2.4 μm wide) possessing a basal plate with a centrally raised keel protruding into an acute tip. The keel was anteriorly widened and had ornamented pores (diameter, 80–95 nm). Numerous small pores ornamented the basal plate. Numerous struts (34–36) regularly extended from the keel toward the scale perimeter. Short struts were present between the primary ones. One to two short struts spurred off from one side of the keel, and a short strut was present on the other side at the scale perimeter. A secondary layer running along the scale perimeter was present between the struts in the upper part of the scale. There were randomly arranged pores (diameter, 80–95 nm) close to the secondary layer. In the present study, this species was detected at stations 4, 6, and 7 Olenyok, at station 9 in Yana, and at stations 27 and 32 in Kolyma (pH = 6.78–7.80, T = 14.5–17.4 °C) (Table 2).
We found Synura sp. 3 (Figure S6i) possessing body scales (4.3 μm long and 1.8 μm wide) bearing a basal plate with a centrally raised keel protruding into an acute tip. The keel, ornamented with small pores, widened anteriorly, then tapered and widened again. Transverse folds interconnected numerous struts (38). In the present study, this species was detected at station 4 in Olenyok (pH = 6.78, T = 17.4 °C) (Table 2).
We found Synura sp. 4 (Figure S7i–l) possessing elongated body scales (3.0–3.8 μm long and 1.1–1.5 μm wide). The keel of the body scales was much reduced and ornamented with small pores (diameter, 34–60 nm). Struts (26–38), interconnected with numerous transverse ribs, regularly extended outward from the keel toward the scale perimeter. Short struts were present between the primary ones. One to two short struts spurred off from one side of the keel, and a short strut was present on the other side at the scale perimeter. In the present study, this species was detected at stations 3, 4, 5, and 7 in Olenyok and at stations 8, 11, and 12 of Yana (pH = 6.78–7.93, T = 14.0–17.4 °C) (Table 2).

4.4. Physicochemical Parameters and Species Composition and Richness

Samples were collected during summer when water temperature of the northern rivers is relatively high (12.0–17.7 °C) (Table 2). Chrysophyceans were not found at the minimum recorded temperature (station 28: 12.0 °C). The maximum number of species was recorded in Kolyma at temperatures of 14.7 °C (station 32, n = 24), 15.9 °C (station 27, n = 16), and 16.7 °C (station 26, n = 16). Similarly, few species were found (n = 3 at station 20 and n = 1 at station 21 in Indigirka) at the maximum recorded temperature (17.7 °C). Amongst the studied rivers, the highest number of species (n = 52) was found in Olenyok at the temperature of 16.8 °C, followed by Yana (n = 31) at the temperature of 14.0 °C, Kolyma (n = 24) at the temperature of 14.7 °C, and Indigirka (n = 5) at the temperature of 14.9 °C. Therefore, the temperature regime of 14.0–16.8 °C may be considered optimal for silica-scaled chrysophytes, under otherwise equal parameters.
During the present study, we recorded neutral to slightly alkalescent pH values in all rivers except Yana. The highest species richness of chrysophyceans was recorded at closely situated stations 4 and 5 in Olenyok at both neutral and slightly acidic pH (6.78, the lowest value at which species were found) or slightly alkaline (7.52) values. Of note, the species composition of silica-scaled chrysophytes at these two stations was different (Table 2).
Low water transparency is typical of Indigirka and Yana, because easily erodible sandy clay and Upper Paleozoic and Mesozoic deposits [94] dominate the drainage basins at the sampling sites in both rivers. Permafrost outcrops and intense thawing of permafrost in summer enhance the processes of bank erosion downstream the rivers, thus increasing the concentration of suspended solids and decreasing water transparency. The highest concentration of suspended solids was recorded in Indigirka (130.40 mg·L−1), and no silica-scaled chrysophytes were found even at turbidity of 21.60 mg·L−1 (Table 2). Water transparency of Kolyma was higher, being comparable to that at some stations in Olenyok. The highest water transparency was recorded in the most species rich location, Olenyok River.
In general, the studied rivers have low salt content, with the exception of Olenyok, which shows moderate salinity. The highest species richness was recorded at salt concentrations ranging from 67.46 (station 4, n = 40) to 228.15 mg·L−1 (station 3, n = 26). This range of salt concentration does not limit the development of silica-scaled chrysophytes.
Silicon concentration showed a wide range (0.5 to 2.5 mg·L−1). Algal blooms are known to depend on silicon concentration in water. The lowest silicon concentration (1.34 mg·L−1) at which chrysophycean species were found was recorded in Kolyma (station 27), whilst the highest concentration (2.48 mg·L−1) was recorded in Olenyok (station 7). Almost all stations in Indigirka showed comparable silicon concentrations (1.9–2.5 mg·L−1), except station 19 (0.5 mg·L−1); however, very low numbers of silica-scaled chrysophytes were recorded in this river (stations 20, n = 3 and station 21, n = 1), and none station 19 with the minimum concentration of silicon (0.5 mg·L−1).
Some samples from Indigirka, Olenyok, and Kolyma showed remarkably high total phosphorus content (168, 186, and 380 μg·L−1, respectively). The total phosphorus content of the majority of the remaining samples was low and did not exceed 100 μg·L−1. Despite the higher phosphorus content at some sites, the biological availability of this element was apparently low. In the present study, chrysophyceans were found at total phosphorus concentrations of between 13 μg·L−1 in Indigirka (stations 20, n = 3 and station 21, n = 1) and 80 μg·L−1 in Kolyma (station 27, n = 15). The maximum number of species (23–52) was recorded at total phosphorus concentrations within 26–180 μg·L−1, which may be optimal for the development of silica-scaled chrysophytes in the study area. The highest total phosphorus concentration at which chrysophyceans were found was recorded in Kolyma (180 μg·L−1; station 32, n = 24) and Olenyok (186 μg·L−1; station 1, n = 14) (Table 2).
There were no considerable variations in the concentration of dissolved oxygen (9.19–10.63 mg·L−1) amongst the sampling stations, indicating a favourable oxygen regimen and good aeration of waters. The concentration of dissolved oxygen within the above limits did not affect the total number of species in the sample. Meanwhile, the concentration of carbon dioxide was low (2.64–7.48 mg·L−1) (Table 2).
Overall, the maximum number of species (n = 52) was recorded at station 5 in Olenyok, where the values of all physical and chemical parameters of the river environment were average (Table 2).
Furthermore, we statistically evaluated the effects of the studied physical and chemical factors of the river environment on species richness and composition of silica-scaled chrysophytes. PERMANOVA (Table 3) revealed that factors including the provenance of a sample from a certain river; distance of the sampling site from the estuary (species composition changed toward the sea); water temperature, transparency, and colour; and concentration of carbon dioxide (CO2), magnesium ions (Mg2+), total iron (Fetotal), surfactants, and oil products, truly affected the species composition of chrysophycean algae in all rivers. Regarding the coefficient of variation (R2), the most important factor determining the species composition in the samples was the provenance of a sample from a certain site of the river (R2 = 0.29). The importance of other factors, even though they reliably affected the species composition in samples, was low (R2 = 0.1).
In NMDS scatter plot (see biplot in Figure 2), the sampling sites in Olenyok formed a close cluster in the right part of the reference frame. This cluster also included three sampling sites in Yana. The species composition of chrysophyceans was similar at all sites in the cluster, but different from that at the remaining sites. The remaining sampling sites were scattered in the left and upper parts of the reference frame. These sites differed in terms of species composition, both from one another and from sampling sites in Olenyok and Yana that formed the above-mentioned cluster.
Three groups of environmental parameters reliably affecting the species composition in samples could be identified based on the direction of gradient vectors (see biplot in Figure 2). The first group included codirectional vectors of water temperature, water transparency, and magnesium ion (Mg2+) concentration, which were positively correlated with one another. The second group included surfactant concentration and distance of a sampling site from the estuary, which were positively correlated with each other. The third group included codirectional vectors of carbon dioxide (CO2), oil product, and total iron (Fetotal) concentration and water colour, which were positively correlated with one another. The vectors of the first and second group were contradirectional and constituted an assembly of negatively correlated parameters. Water temperature, water transparency, and magnesium ion (Mg2+) concentration increased but surfactant concentration decreased as the sampling sites approached the river estuary. Meanwhile, the environmental parameters in the first and second groups (water temperature and transparency, magnesium ion (Mg2+) and surfactant concentration, and distance of a sampling site from the estuary) constituted an assembly of factors that complementarily affected the species composition in samples. Parameters in the third group (carbon dioxide (CO2), oil product, and total iron (Fetotal) concentration and water colour) exerted another type of effect, which was not related to the mechanisms of effect of parameters in the first and second groups.
The species richness in samples increased with increasing water temperature, water transparency, and magnesium ion (Mg2+) concentration, as reflected by the direction of gradient vector on the NMDS scatter plot (Figure 2). Similarly, the species richness in samples increased as the sampling site approached the estuary (i.e., the sea) but decreased as surfactant concentration increased. The other physical and chemical parameters, within the limits of their variation range, did not significantly affect species composition in the samples.

4.5. Specific Morphology of Some Chrysophytes Found in the Study Area

The present investigation of the species richness of silica-scaled chrysophytes revealed some peculiarities in the structure of scales in Spiniferomonas bourrellyi, Synura mammillosa, and Synura echinulate.
Synura echinulate is a cosmopolitan species that mainly occurs in oligotrophic and slightly acidic waters. In Synura echinulate collected in the present study, the scales occasionally lacked the labyrinth pattern (Figure S5i, station 12) on the distal third and instead presented a dotted pattern, which was only slightly similar to the former pattern (Figure S5j, station 5). In some other samples, the scales completely lacked the papillar ornament (Figure S5k, station 32) and possessed a weakly pronounced, upturned rim surrounding one-half to two-thirds of the scale (Figure S5j,k). Given the sufficient concentration of silicon in the waters of these three stations (1.44–1.90 mg·L−1) (Table 1), such alterations may be explained by variations in the degree of silification of Synura echinulate scales. These alterations have also been reported to occur in vitro [95].
Synura mammillosa is a boreal species. In Synura mammillosa collected in the present study, some scales bore transverse ribs on the marginal part (Figure S5m) and some present a pattern of contradirectional ribs (Figure S5l). The presence of such ornamentation on the scales of Synura mammillosa has been previously described in Aquitaine, France ([96]; Figure 52).
Spiniferomonas bourrellyi is a cosmopolitan species. Interestingly, we noted differences in the constructions and length of spines on the scales of this species collected in the present study. First, some samples possessed a 14–16.5-μm-long spine with a conical base and tapering toward the tip (Figure S3c). This type of spine was described by Green [97] in a new species of Spiniferomonas (Spiniferomonas andersonii Green). However, Nicholls [87] synonymised Spiniferomonas andersonii as Spiniferomonas bourrellyi. Second, in some samples, the spine was only 2.3–3.0 μm long and possessed a saucer-shaped base (diameter, 2.2–2.5 μm); this spine slightly widened at the base and then tapered toward the tip (Figure S3e). This type of spine has been described in samples collected from lakes in Connecticut [90] and the Labynkyr and Vorota lakes in Yakutia, Russia [48]. However, in Spiniferomonas bourrellyi collected from Olenyok in the present study, the spine was considerably smaller than that described previously.

5. Discussion

The considerable diversity (82 taxa) of silica-scaled chrysophytes in the Arctic rivers we investigated is comparable to the diversity in another Arctic region–the Bolshezemelskaya tundra of Russia (75 taxa) [22]. The species richness of silica-scaled chrysophytes is often high in waters with slightly acidic pH and low dissolved salt content, and in oligotrophic to mesotrophic waters with moderate content of dissolved humic compounds [98,99,100,101]. The geochemical context of the development of waters of the Bolshezemelskaya tundra is, however, considerably different. The waterbodies in this tundra are small lakes and ponds, that is, closed natural systems, with different levels of trophicity (both oligotrophic lakes and lakes with anthropogenic load), different acidic to alkaline pH values (5.5–8), and wide ranges of conductivity (19–146 μS·cm−1) and temperature (8–19.7 °C) [22]. Waters of the rivers studied here showed comparable values of hydrochemical parameters (Table 2). The studied rivers of Yakutia, despite the increased phosphorus content, can be classified as oligotrophic, according to the classification proposed by Dodds et al. [102]. In oxygen/iron-rich waters, biologically available phosphorus accounts for only 8% of total phosphorus [102]. The pH and temperature at different stations range within 6.69–7.93 (neutral to slightly alkaline) and 12–17 °C, respectively. Likewise, a high species richness of silica-scaled chrysophytes (79 species) has been previously recorded in the Baikal Region with similar primary hydrochemical parameters [13,103]. pH is one of the major factors driving the development and distribution of silica-scaled chrysophytes [3]. In the present study, however, we obtained data that water pH at different stations, ranging from neutral to alkalescent or alkaline values, did not affect the species composition and richness of silica-scaled chrysophytes in the rivers of Yakutia. Summarising data obtained by the statistical analysis of 28 parameters of the habitat, the major factors governing species composition and the richness of silica-scaled chrysophytes were the provenance of the sample from a specific river; the distance of the sampling site from the estuary (i.e., the sea); water temperature, transparency, and colour; and concentration of carbon dioxide, magnesium ions, and total iron. The factor related to sampling from a specific river confirms the previously identified regional nature of the diversity of silica-scaled chrysophytes in northern reservoirs [40]. However, the mouths of the rivers we examined are located in the same region, and therefore these differences may have been driven by the diversity inherited by the rivers, which flow over a vast territory. Further, water pollution from surfactants and petroleum products have a considerable impact on species diversity. While settlements and anthropogenic factors other than poor navigation do not exist in the region, active development of the region and increased water pollution with surfactants and petroleum products can have a negative impact on the diversity of silica-scaled chrysophytes. The reaction of silica-scaled chrysophytes to changes in the concentration of magnesium and iron ions was intriguing, as this reaction was previously undetected. We have formerly discussed the importance of transparency and colour to the development of silica-scaled chrysophytes and phytoplankton in general [103,104,105]. Increases in water turbidity as a result of warming and the removal of terrigenous material by meltwater play significant roles in inhibiting the development of chrysophytes. Accordingly, we did not observe chrysophytes in rivers with high turbidity in this study.
The silica-scaled chrysophytes are known to produce different phenotypes of scales in response to different environmental conditions [95,96,106]. However, some species exhibit such phenotypic plasticity more often than others do. Previous geometric morphometric analyses have revealed the phenotypic plasticity of scales in Synura echinulate grown with different combinations of light, temperature, and pH [95]. Variations in the shapes of spine bases have been described in the Spiniferomonas bourrellii complex [87]. Despite the observed morphological differences in the structure of individual scales, taxonomically relevant traits of species remained intact. Therefore, such morphological alterations are a result of phenotypic plasticity, likely induced by one or many environmental factors.
Analysing the species composition of silica-scaled chrysophytes in other northern reservoirs located above the Arctic Circle on a time scale, we observed that cosmopolitans and widespread species predominate in the longitude group (Figure 3).
As such, the proportion of cosmopolitan species amongst all chrysophytes has changed from 68% in the Khantay Reservoir of Russia during 1981–1982 [28] to 30.5% in the Arctic rivers of Yakutia in the present study. Simultaneously, however, the proportion of ubiquitous species has increased from 20% to 30.5%. Species with scattered distribution, including rare species, became abundant during 1981–1982. The maximum number of bipolar species (14.3%) has been recorded in the Vorkuta and Bolshezemelskaya tundra of Russia [22]. Synura obesa Němcova, Novakova & Řezačova-Škaloudova remains restricted to Abisko in Swedish Lapland [19], which is the only species that has retained its endemicity amongst all other endemics. In the present study, for the first time since its description, we found Synura petersenii f. taymyrensis, which was suggested to be endemic to the Taymyr Peninsula [14]. Owing to a high number of potentially new species (17%), the proportion of species with unidentified geographic characteristics was the highest in the present study.
Polyzonal species were predominant in the latitudinal group. The proportion of these species amongst all chrysophytes varies between 84.6% in the Taymyr Peninsula of Russia in 1993 [9] and 51.2% in the Arctic rivers of Yakutia in the present study. The maximum number of Arctic-boreal species (32.7%) was recorded in the Bolshezemelskaya tundra of Russia in 2003 [22], and the minimum (8%) was recorded in the Khantay Reservoir of Russia [28]. We noted a decrease in the proportion of Arctic-boreal species and an increase in the proportion of boreal species of silica-scaled chrysophytes in the waters of the studied Arctic Region (Figure 3) compared with previously published data. Boreal species found in the studied rivers included Chrysosphaerella rotundata, Paraphysomonas cf. punctata ssp. simplicior, Paraphysomonas vulgaris Lepidochromonas eiffelii, Lepidochromonas elegantissima, Spiniferomonas takahashii, Synura curtispina f. reticulata, Synura glabra, and Synura mammillosa. In addition, two species, namely Synura petersenii f. taymyrensis and Synura obesa, can be considered Arctic. A high proportion of species with unidentified geographic characteristics (14.2%) was recorded in the Vorkuta tundra of Russia during 2002–2003 [22] due to the high number of bipolar species. In the present study, this proportion was 19.5% due to the high number of potentially new species.

6. Conclusions

By establishing the biogeographic structure of silica-scaled chrysophyte communities in the mouths of rivers in the Arctic region of Yakutia, and comparing these with similar structures in other previously studied northern reservoirs, we detected in the temporal aspect a decrease in the proportion of cosmopolitans and an increase in the proportion of widespread species and those with scattered distribution. Specifically, we noted a decrease in the proportion of polyzonal species and an increase in the percentage of boreal and arcto-boreal species. This confirms our hypothesis that species have been introduced from temperate latitudes into Arctic reservoirs in the context of global climate change. The record level of species diversity of these microeukaryotes and the change in the structure of their communities reflect the beginning of manifestations of climate warming. The results of the study are important for the present and further monitoring of an untouched, environmentally sensitive area of the world, and are relevant for assessing the future consequences of local anthropogenic impacts and climate change.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/su132413768/s1, Table S1: Species composition of silica-scaled chrysophytes north of the Arctic Circle, Table S2: Raw data table for statistical analysis, Figure S1: Silica-scaled chrysophytes of the genera Chrysosphaerella, Clathromonas and Paraphysomonas identified in this study: a, b, Chrysosphaerella brevispina (site 4); c, d, Chrysosphaerella sp. (c site 5; d site 12); e, C. coronacircumspina (site 6); f, C. rotundata (site 5); g, C. longispina (site 4); h, Paraphysomonas sp. 1 (site 6); i, P. vulgaris (site 2); j, P. acuminata acuminata (site 4); k, Paraphysomonas sp. 2 (site 5); l, P. bandaiensis (site 5); m, P. gladiata (site 12); n, P. uniformis hemiradia (site 6); o, P. cf. punctata ssp. simplicior (site 12); p, P. punctata (site 5); q, L. circumvallata (site 5). Scale bars: j–l–500 nµ; a–e, h, i, m–q–1 µm; f, g–2 µm, Figure S2: Silica-scaled chrysophytes of the genera Paraphysomonas and Lepidochromonas identified in this study: a–c, Paraphysomonas sp. 3 (a, b convex side of scale proximal, c distal side of scale (a site 4; b, c site 5)); d, Paraphysomonas sp. 4 (site 5); e, Lepidochromonas quadrispina (site 5); f, g, L. eiffelii (site 5); h, L. elegantissima (site 5); i, L. butcheri (site 5); j, L. diadernifera (site 5); k, L. undulata (site 5); l, Lepidochromonas sp. 1 (site 5); m, L. subquadrangularis (site 5); n, Lepidochromonas sp. 2 (site 5); o, L. coronata (site 5). Scale bars: e, m–200 nµ; f, g, i–l, n–500 nµ; a–d, h, o–1 µm, Figure S3: Silica-scaled chrysophytes of the genus Spiniferomonas identified in this study: a, Spiniferomonas trioralis (site 7); b, S. trioralis f. cuspidata (site 4); c, e, S. bourrellyi (c site 5; e site 6); d, f, S. serrata (site 3); g, S. abei (site 5); h, S. cornuta (site 12); i, S. silverensis (site 2); j, S. conica (site 2); k, S. minuta (site 4); l, S. takahashii (site 5). Scale bars: h–500 nµ; l–1 µm; a–g, i–k–2 µm, Figure S4: Silica-scaled chrysophytes of the genus Mallomonas identified in this study: a, Mallomonas acaroides (site 3); b, M. akrokomos (site 12); c, M. actinoloma (site 12); d, M. alata f. hualvensis (site 5); e, M. alpina (site 3); f, g, M. kuzminii (f site 3; g site 11); h, M. tonsurata (site 32); i, M. elongata (site 12); j, M. areolata (site 2); k, M. crassisquama (site 26); l, M. cratis (site 3); m, n, M. caudata (site 26); o, M. sf. pumilio (site 4); p, M. eoa (site 5). Scale bars: b–1 µm; a, c–p–2 µm, Figure S5: Silica-scaled chrysophytes of the genera Mallomonas and Synura identified in this study: a, Mallomonas lychenensis (site 4); b, M. multiunca (site 2); c, d, M. punctifera (site 4); e, M. parvula (site 11); f, M. papillosa (site 4); g, M. striata (site 26); h, Synura uvella (site 26); i–k, S. echinulate (i site 12; j site 5; k site 32); l, m, S. mammillosa (site 11); n, S. leptorhabda (site 12); o, Synura sp. 1 (site 12); p, S. nygaardii (site 3). Scale bars: e, f–1 µm; a–d, g–p–2 µm, Figure S6: Silica-scaled chrysophytes of the genus Synura identified in this study: a, Synura sphagnicola (site 5); b, S. spinosa f. longispina (site 30); c, S. spinosa (site 7); d, S. curtispina f. reticulata (site 3); e, S. curtispina (site 5); f, S. petersenii (site 11); g, h, S. borealis (site 2); i, Synura sp. 3 (site 4); j–l, S. americana (j site 2; k, l site 7); m, n, S. sf. cornuta (m site 1; n site 3); o–q, S. sf. vinlandica (site 7); r, s, S. conopea (r site 5; s site 31). Scale bars 2 µm, Figure S7: Silica-scaled chrysophytes of the genus Synura identified in this study: a, b, Synura sf. soroconopea (site 3); c, d, Synura sp. 2 (c site 9; d site 7); e, S. macropora (site 11); f, S. petersenii f. taymyrensis (site 5); g, S. praefracta (site 5); h, S. glabra (site 4), i–l, Synura sp. 4 (i site 4; g site 5; k, l site 11). Scale bars 2 µm.

Author Contributions

Conceptualisation, A.B. and Y.L.; data curation, V.G., O.G., A.B., Y.B. and A.F. writing—original draft preparation, A.B. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

The work was conducted within the State Project 0297-2021-0023 ‘Vegetation of the permafrost taiga zone of Yakutia: biodiversity, environment-forming functions, protection and rational use’ (Institute for Biological Problems of Cryolithozone, organisation and execution of expeditions, water chemistry analyses), and the State Project 0279-2021-0008 (Limnological Institute, Siberian Branch of the Russian Academy of Science, study of chrysophytes and statistical analyses).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

Acknowledgments

The microscopy studies were performed in the Electron Microscopy Center of the Shared Research Facilities ‘Ultramicroanalysis’ of Limnological Institute (http://www.lin.irk.ru/copp/).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wujek, D.E. Freshwater silica-scaled heterotrophic protista: Heliozoa, thaumatomonad flagellates, amoebae, and bicosoecids, from the Lake Itasca Region, Minnesota. Minn. Acad. Sci. J. 2015, 79, 1–13. [Google Scholar]
  2. Siver, P.A. The Biology of Mallomonas: Morphology, Taxonomy and Ecology; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1991; p. 248. [Google Scholar]
  3. Siver, P.A. Synurophyte algae. In Freshwater Algae of North America: Ecology and Classification, 2nd ed.; Academic Press: Boston, MA, USA, 2015; pp. 607–651. [Google Scholar]
  4. Dürrschmidt, M.; Croome, R. Mallomonadaceae (Chrysophyceae) from Malaysia and Australia. Nord. J. Bot. 1985, 5, 285–298. [Google Scholar] [CrossRef]
  5. Neustupa, J.; Řezáčová, M. The genus Mallomonas (Mallomonadales, Synurophyceae) in several Southeast Asian urban water bodies—The biogeographical implications. Nova Hedwig. 2007, 84, 249–259. [Google Scholar] [CrossRef]
  6. Gusev, E.S.; Thant, N.T.H. Silica-scaled chrysophytes (Chrysophyceae and Synurophyceae) from Vietnam (Khanh Hoa and Quang Nam provinces). Nova Hedwig. 2011, 93, 191–199. [Google Scholar] [CrossRef]
  7. Gusev, E.S.; Kapustin, D.A.; Martynenko, N.A.; Guseva, E.E.; Kulikovskiy, M.S. Mallomonas gusakovii sp. nov. (Chrysophyceae, Synurales), a new species from Phu Quoc Island, Vietnam. Phytotaxa 2019, 406, 199–205. [Google Scholar] [CrossRef]
  8. Duff, K.E. Chrysophyte microfossils in arctic Siberian lakes. Chrysophytes: Progress and new horizons. Nova Hedwig. Beih. 1996, 114, 253–263. [Google Scholar]
  9. Kristiansen, J.; Düwel, L.; Wegeber, S. Silica-scaled chrysophytes from the Taymyr Peninsula, Northern Siberia. Nova Hedwig. 1997, 65, 337–351. [Google Scholar] [CrossRef]
  10. Gusev, E.S.; Guseva, E.E.; Gabyshev, V.A. Taxonomic composition of silica-scaled chrysophytes in rivers and lakes of Yakutia and Magadanskaya oblast (Russia). Nova Hedwig. Beih. 2018, 147, 105–117. [Google Scholar] [CrossRef]
  11. Hällfors, G.; Hällfors, S. Records of chrysophytes with siliceous scales (Mallomonadaceae and Paraphysomonadaceae) from Finnish inland waters. Flagellates in freshwater ecosystems. Hydrobiologia 1988, 161, 1–29. [Google Scholar] [CrossRef]
  12. Němcová, Y.; Kreidlová, J.; Kosová, A.; Neustupa, J. Lakes and pools of Aquitaine region (France)—A biodiversity hotspot of Synurales in Europe. Nova Hedwig. 2012, 95, 1–24. [Google Scholar] [CrossRef] [PubMed]
  13. Bessudova, A.Y.; Sorokovikova, L.M.; Tomberg, I.V.; Likhoshway, Y.V. Silica-scaled chrysophytes in large tributaries of Lake Baikal. Cryptogam. Algol. 2018, 39, 145–165. [Google Scholar] [CrossRef]
  14. Siver, P.A. The synurophyceae. In Freshwater Algae of North America; Wehr, J., Sheath, B., Eds.; Academic Press: Thornwood, NY, USA, 2003; pp. 523–558. [Google Scholar]
  15. Asmund, B.; Hilliard, D.K. Studies on Chrysophyceae from some ponds and lakes in Alaska. Hydrobiolology 1961, 7, 237–258. [Google Scholar] [CrossRef]
  16. McKenzie, C.; Kling, H. Scale-bearing Chrysophyceae (Mallomonadaceae and Paraphysomonadaceae) from Mackenzie Delta area lakes, Northwest Territories, Canada. Nord. J. Bot. 1989, 9, 103–112. [Google Scholar] [CrossRef]
  17. Kristiansen, J. Silica-scaled chrysophytes from west Greenland: Disko island and the Søndre Strømfjord region. Nord. J. Bot. 1992, 12, 525–536. [Google Scholar] [CrossRef]
  18. Forsström, L.; Sorvari, S.; Korhola, A.; Rautio, M. Seasonality of phytoplankton in subarctic Lake Saanajärvi in NW Finnish Lapland. Pol. Biol. 2005, 28, 846–861. [Google Scholar] [CrossRef]
  19. Němcová, Y.; Nováková, S.; Řezáčová-Škaloudová, M. Synura obesa sp. nov. (Synurophyceae) and other silica scaled chrysophytes from Abisko (Swedish Lapland). Nova Hedwig. 2008, 86, 243–254. [Google Scholar] [CrossRef] [Green Version]
  20. Pichrtová, M.; Janatková, K.; Němcová, Y. Silica-scaled chrysophytes from Abisko (Swedish Lapland). Nord. J. Bot. 2011, 29, 112–118. [Google Scholar] [CrossRef]
  21. Kristiansen, J. Flagellates from Finnish Lappland. Bot. Tidsskr. 1964, 59, 315–333. [Google Scholar]
  22. Siver, P.A.; Voloshko, L.N.; Gavrilova, O.V.; Getsen, M.V. The scaled chrysophyte flora of the Bolshezemelskaya tundra (Russia). Nova Hedwig. Beih. 2005, 128, 125–150. [Google Scholar]
  23. Voloshko, L.N. A new species of the genus Mallomonas (Chrysophyceae, Synurophyceae) from lakes of the Vorkuta tundra. Bot. Zhurn. 2012, 97, 1226–1234. (In Russian) [Google Scholar]
  24. Voloshko, L.N. New taxa of the genus Mallomonas (Chrysophyceae, Synurophyceae) from lakes of the Polar Urals. Bot. Zhurn. 2009, 94, 1068–1076. (In Russian) [Google Scholar]
  25. Voloshko, L.N. The chrysophycean algae from glacial lakes of Polar Ural (Russia). Nova Hedwig. Beih. 2010, 136, 191–211. [Google Scholar] [CrossRef]
  26. Firsova, A.D.; Bessudova, A.Y.; Sorokovikova, L.M.; Tomberg, I.V.; Likhoshway, Y.V. The diversity of chrysophycean algae in an arctic zone of river and sea water mixing, Russia. Am. J. Plant Sci. 2015, 6, 2439–2452. [Google Scholar] [CrossRef] [Green Version]
  27. Bessudova, A.Y.; Bukin, Y.S.; Sorokovikova, L.M.; Firsova, A.D.; Tomberg, I.V. Silica-scaled Chrysophytes in small lakes of the Lower Yenisei basin, the Arctic. Nova Hedwig. 2018, 107, 315–336. [Google Scholar] [CrossRef]
  28. Balonov, I.M.; Kuzmina, A.E. Golden algae. In Proceedings of Limnological Institute, RAS of the Academy of Sciences of USSRP; Votintsev, K.K., Ed.; Hydrochemical and Hydrobiological Studies of the Khantay Reservoir, Nauka, Siberian Branch: Novosibirsk, Russia, 1986; pp. 59–70. (In Russian) [Google Scholar]
  29. Wolfe, A.P.; Siver, P.A. A hypothesis linking chrysophyte microfossils to lake carbon dynamics on ecological and evolutionary time scales. Glob. Planetary Chang. 2013, 111, 189–198. [Google Scholar] [CrossRef]
  30. Mushet, G.R.; Laird, K.R.; Das, B.; Hesjedal, B.; Leavitt, P.R.; Scott, K.A.; Simpson, G.L.; Wissel, B.; Wolfe, J.D.; Cumming, B.F. Regional climate changes drive increased scaled-chrysophyte abundance in lakesdownwind of Athabasca Oil Sands nitrogen emissions. J. Paleolimnol. 2017, 58, 419–435. [Google Scholar] [CrossRef] [Green Version]
  31. Schindler, D.W. The cumulative effects of climate warming and other human stresses on Canadian freshwaters in the new millennium. Can. J. Fish. Aquat. Sci. 2001, 58, 18–29. [Google Scholar] [CrossRef]
  32. Paterson, A.M.; Winter, J.G.; Nicholls, K.H.; Clarks, B.J.; Ramcharan, C.W.; Yan, N.D.; Somers, K.M. Long-term changes in phytoplankton composition in seven Canadian Sheet lakes in response to multiple anthropogenic stressors. Can. J. Fish. Aquat. Sci. 2008, 65, 846–861. [Google Scholar] [CrossRef]
  33. Rühland, K.M.; Paterson, A.M.; Smol, J.P. Hemispheric-scale patterns of climate-related shifts in planktonic diatoms from North American and European lakes. Glob. Chang. Biol. 2008, 14, 2740–2754. [Google Scholar] [CrossRef]
  34. Field, C.B.; Barros, V.R.; Dokken, D.J.; Mach, K.J.; Mastrandrea, M.D.; Bilir, T.E.; Chatterjee, M.; Ebi, K.L.; Estrada, Y.O.; Genova, R.C.; et al. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. In Climate Change 2014: Impacts, Adaptation, and Vulnerability; Cambridge University Press: New York, NY, USA, 2014; pp. 1–32. [Google Scholar]
  35. Magnuson, J.J.; Robertson, D.M.; Benson, B.J.; Wynne, R.H.; Livingstone, D.M.; Arai, R.A.; Barry, R.G.; Card, V.; Kuusisto, E.; Granin, N.G.; et al. Historical trends in lake and river ice cover in the Northern Hemisphere. Science 2000, 289, 1743–1746. [Google Scholar] [CrossRef] [Green Version]
  36. Vuglinsky, V.; Valatin, D. Changes in ice cover duration and maximum ice thickness for rivers and lakes in the Asian part of Russia. Nat. Res. 2018, 9, 73–87. [Google Scholar] [CrossRef] [Green Version]
  37. Korneva, L.G. Recent invasion of planktonic diatom algae in the Volga River basin and their causes. Inland Wat. Boil. 2007, 1, 30–39. (In Russian) [Google Scholar]
  38. Korneva, L.G. Invasions of alien species of planktonic microalgae into the fresh waters of Holarctic (Review). Russ. J. Biol. Invasions 2014, 1, 9–37. [Google Scholar] [CrossRef]
  39. Korneva, L.G. Phytoplankton of Reservoirs of the Volga River Basin; Kostroma: Dom Pechati, Russia, 2015; p. 284. (In Russian) [Google Scholar]
  40. Bessudova, A.; Bukin, Y.; Likhoshway, Y. Dispersal of silica-scaled Chrysophytes in Northern Water bodies. Diversity 2021, 13, 284. [Google Scholar] [CrossRef]
  41. Kirillov, F.N.; Labutina, T.M. Biology of Vilyuy Reservoir; Nauka: Novosibirks, Russia, 1979; p. 272. (In Russian) [Google Scholar]
  42. Vasilyeva, I.I.; Remigaylo, P.A. Algae of Vilyuy Reservoir; Yakutia Scientific Center USSR Academy of Sciences Publishing House: Yakutsk, Russia, 1982; p. 115. (In Russian) [Google Scholar]
  43. Remigaylo, P.A. Phytoplankton of the Vilyuy River and Vilyuy Reservoir. Ph.D. Thesis, Novosibirsk University, Novosibirsk, Russia, 1995; p. 16. (In Russian). [Google Scholar]
  44. Vasilyeva-Kralina, I.I.; Remigaylo, P.A.; Gabushev, V.A.; Pshennikova, E.V.; Ivanova, A.P.; Kopyrina, L.I. Algae. In Diversity of Plants of Yakutia; Publishing House SB RAS: Novpsibirsk, Russia, 2005; pp. 150–272. (In Russian) [Google Scholar]
  45. Gabyshev, V.A.; Gabysheva, O.I. Phytoplankton of Big Rivers of Yakutia and Adjacent Territories of East Siberia; Korneva, L.G., Ed.; SibAK Press: Novosibirks, Russia, 2018; p. 414. (In Russian) [Google Scholar]
  46. Gabyshev, V.A.; Tsarenko, P.M.; Ivanova, A.P. Diversity and features of the spatial structure of algal communities of water bodies and watercourses in the Lena River estuary. Inland Water Biol. 2019, 12, 1–9. [Google Scholar] [CrossRef]
  47. Gabyshev, V.A.; Tsarenko, P.M.; Ivanova, A.P. Algae of mouth area of Lena River. In The Biological Resources of the Ust-Lensky Nature Reserve: Fungi, Algae, Vegetation, Fishes, Birds, Muskoxen; Science: Novosibirsk, Russia, 2019; pp. 14–35. [Google Scholar]
  48. Bessudova, A.Y.; 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 Hedwig. Beih. 2019, 148, 35–48. [Google Scholar] [CrossRef]
  49. Izyumenko, S.A. (Ed.) Climate of the Yakut ASSR. In Atlas; Gidrometeoizdat: Leningrad, Russia, 1968; p. 33. (In Russian) [Google Scholar]
  50. Chistyakov, G.E. River Water Resources of Yakutia; Nauka: Moskow, Russia, 1964; p. 255. (In Russian) [Google Scholar]
  51. Chistyakov, G.E. Hydropower Resources of Drainage Basins of the Yana River; Chistyakov, G.E., Nogovitsyn, D.D., Yakushev, M.V., Eds.; Nauka: Moskow, Russia, 1970; p. 214. (In Russian) [Google Scholar]
  52. Resources of Inland Waters of the USSR. Lena-Indigirka Region; Gidrometeoizdat: Leningrad, Russia, 1972; Volume 17, p. 651. (In Russian)
  53. Resources of Inland Waters of the USSR. North-East; Gidrometeoizdat: Leningrad, Russia, 1966; Volume 19, p. 602. (In Russian)
  54. Guidelines for Collection and Processing of Samples at Hydrobiological Studies of Fresh Waters; Nauka: Leningrad, Russia, 1981; p. 32. (In Russian)
  55. Wiebe, P.H.; Benfield, M.C. From the Hensen net toward four-dimensional biological oceanography. Prog. Oceanogr. 2003, 56, 7–136. [Google Scholar] [CrossRef]
  56. Alyokin, O.A.; Semyonov, A.D.; Skopintsev, B.A. Guidelines for Chemical Analysis of Inland Waters; Gidrometeoizdat: Leningrad, Russia, 1973; p. 269. (In Russian) [Google Scholar]
  57. Semyonov, A.D. Guidelines for Chemical Analysis of Inland Water; Gidrometeizdat: Leningrad, Russia, 1977; p. 540. (In Russian) [Google Scholar]
  58. Sokolova, S.A.; Anisova, S.N. (Eds.) List of Commercial Fishing Standards: Maximum Allowable Concentrations (MACs) and Safe Reference Levels of Contaminants in Commercial Fishing Waters; VNIRO: Moskow, Russia, 1999; p. 304. (In Russian) [Google Scholar]
  59. Kristiansen, J. Dispersal and biogeography of silica-scaled chrysophytes. Biodiv. Conserv. 2008, 17, 419–426. [Google Scholar] [CrossRef]
  60. Voloshko, L.N. Golden Algae of Waters of North Russia; Komarov Botanical Institute of the Russian Academy of Sciences: Sankt-Petersburg, Russia, 2017; p. 378. (In Russian) [Google Scholar]
  61. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Aust. Ecol. 2001, 26, 32–46. [Google Scholar] [CrossRef]
  62. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; Stevens, M.H.H.; Wagner, H. Vegan: Community Ecology Package, 2014. R Package Version 2.2-0. Available online: http://CRAN.Rproject.org/package=vegan (accessed on 13 August 2021).
  63. Oksanen, J. Vegan: An Introduction to Ordination. 2015. Available online: http://cran.r-project.org/web/packages/vegan/vignettes/introvegan.pdf (accessed on 13 August 2021).
  64. Škaloud, P.; Škaloudová, M.; Jadrná, I.; Bestová, H.; Pusztai, M.; Kapustin, D.A.; Siver, P.A. Comparing morphological and molecular estimates of species diversity in the freshwater genus Synura (stramenopiles): A model for understanding diversity of eukaryotic microorganisms. J. Phycol. 2020, 56, 574–591. [Google Scholar] [CrossRef] [PubMed]
  65. Jo, B.Y.; Kim, J.I.; Škaloud, P.; Siver, P.A.; Shin, W. Multigene phylogeny of Synura (Synurophyceae) and descriptions of four new species based on morphological and DNA evidence. Europ. J. Phycol. 2016, 51, 413–430. [Google Scholar] [CrossRef] [Green Version]
  66. Škaloudová, M.; Škaloud, P. A new species of Chrysosphaerella (Chrysophyceae: Chromulinales), Chrysosphaerella rotundata, sp. nov., from Finland. Phytotaxa 2013, 130, 34–42. [Google Scholar] [CrossRef]
  67. Dürrschmidt, M. Studies on scale-bearing Chrysophyceae from the Giessen area, Federal Republic of Germany. Nord. J. Bot. 1984, 4, 123–143. [Google Scholar] [CrossRef]
  68. Thomsen, H.A.; Zimmermann, B.; Moestrup, O.; Kristiansen, J. Some new freshwater species of Paraphysomonas (Chrysophyceae). Nord. J. Bot. 1981, 1, 559–581. [Google Scholar] [CrossRef]
  69. Kristiansen, J. On the species of Paraphysomonas (Chrysophyceae) in some Greek lakes. Nova Hedwig. 1983, 38, 65–72. [Google Scholar]
  70. Kling, H.J.; Kristiansen, J. Scale-bearing Chrysophyceae (Mallomonadaceae) from Central and Northern Canada. Nord. J. Bot. 1983, 3, 269–290. [Google Scholar] [CrossRef]
  71. Nicholls, K.H. Five Paraphysomonas species (Chrysophyceae) new to North America, with notes on three other rarely reported species. Can. J. Bot. 1985, 63, 1208–1212. [Google Scholar] [CrossRef]
  72. Hickel, B.; Maass, I. Scaled chrysophytes, including heterotrophic nanoflagellates from the lake district in Holstein, northern Germany. Nova Hedwig. Beih. 1989, 95, 233–257. [Google Scholar]
  73. Wujek, D.E.; O’Kelly, C.J. Silica-scaled Chrysophyceae (Mallomonadaceae and Paraphysomonadaceae) from New Zealand freshwaters. N. Z. J. Bot. 1992, 30, 405–414. [Google Scholar] [CrossRef] [Green Version]
  74. Ikävalko, J. Observations on silica-scaled flagellates (Chrysophyceae and Synurophyceae) in the brackish water of Pojo Bay, SW coast of Finland. Ann. Bot. Fen. 1994, 31, 1–27. [Google Scholar]
  75. Hällfors, G. Checklist of Baltic Sea Phytoplankton species (including some heterotrophic protistan groups). Balt. Sea Environ. Proc. 2004, 95, 1–208. [Google Scholar]
  76. Preisig, H.R.; Hibberd, D.J. Ultrastructure and taxonomy of Paraphysomonas (Chrysophyceae) and related genera 2. Nord. J. Bot. 1982, 2, 601–638. [Google Scholar] [CrossRef]
  77. Dürrschmidt, M. Studies on the Chrysophyceae from South Chilean inland waters by means of scanning and transmission electron microscopy, II. Archiv für Hydrobiologie 63. Algol. Stud. 1982, 63, 121–163. [Google Scholar]
  78. Nicholls, K.H. Eight chrysophyceae new to North America. Phycologia 1984, 23, 213–221. [Google Scholar] [CrossRef]
  79. Barreto, S. The silica-scaled chrysophyte flora of Hungary. Nova Hedwig. Beih. 2005, 128, 11–41. [Google Scholar]
  80. Pichrtová, M.; Němcová, Y.; Škaloud, P.; Rott, E. Silica-scaled chrysophytes from North Tyrol (Austria) including a description of Mallomonas tirolensis sp. nov. Nova Hedwig. Beih. 2013, 142, 69–85. [Google Scholar]
  81. Wei, Y.-X.; Yuan, X.-P. Studies on silica-scaled chrysophytes from the Daxinganling mountains and Wudalianchi Lake regions, China. Nova Hedwig. 2015, 101, 299–312. [Google Scholar] [CrossRef]
  82. Nielsen, Y. Spiniferomonas abrupta, sp. nov. and some rare species of Synurophyceae and Chrysophyceae, not formerly recorded from Denmark. Nord. J. Bot. 1994, 14, 473–480. [Google Scholar] [CrossRef]
  83. Ikävalko, J.; Thomsen, H.A. Scale-covered and loricate flagellates (Chrysophyceae and Synurophyceae) from Baltic Sea ice. Nova Hedwig. Beih. 1994, 114, 147–160. [Google Scholar]
  84. Ikävalko, J. Contribution to the flora of silica-scaled flagellates in Mikkeli, central Finland. Nova Hedwig. 1994, 58, 475–505. [Google Scholar]
  85. Finlay, B.J.; Clarke, K.J. Apparent global ubiquity of species in the protist genus Paraphysomonas. Protist 1999, 150, 419–430. [Google Scholar] [CrossRef]
  86. Takahashi, E. Studies on genera Mallomonas and Synura, and other plankton in fresh water with the electron microscope. VII. New genus Spiniferomonas of the Synuraceae (Chrysophyceae). Bot. Mag. Tok. 1973, 86, 75–88. [Google Scholar] [CrossRef]
  87. Nicholls, K.H. Spiniferomonas (Chrysophyceae) in Ontario lakes including a revision and descriptions of two new species. Can. J. Bot. 1981, 59, 107–117. [Google Scholar] [CrossRef]
  88. Jacobsen, V.A. Scale-bearing Chrysophyceae (Mallomonadaceae and Paraphysomonadaceae) from West Greenland. Nord. J. Bot. 1985, 5, 381–398. [Google Scholar] [CrossRef]
  89. Nicholls, K.H. Descriptions of Spiniferomonas silverensis sp. nov. and S. minuta sp. nov. and an assessment of form variation in their closest relative, S. trioralis (Chrysophyceae). Can. J. Bot. 1984, 62, 2329–2335. [Google Scholar] [CrossRef]
  90. Siver, P.A. The distribution and ecology of Spiniferomonas (Chrysophyceae) in Connecticut (U.S.A). Nord. J. Bot. 1988, 8, 205–212. [Google Scholar] [CrossRef]
  91. Němcová, Y.; Pusztai, M.; Škaloudová, M.; Neustupa, J. Silica-scaled chrysophytes (Stramenopiles, Ochrophyta) along a salinity gradient: A case study from the Gulf of Bothnia western shore (Northern Europe). Hydrobiologia 2016, 764, 187–197. [Google Scholar] [CrossRef]
  92. Siver, P.A.; Wujek, D.E. Scaled Chrysophyceae and Synurophyceae from Florida, U.S.A.: VI. Observations on the flora from waterbodies in the Ocala National Forest. Nova Hedwig. 1999, 68, 75–92. [Google Scholar] [CrossRef]
  93. Bessudova, A.Y.; Domysheva, V.M.; Firsova, A.D.; Likhoshway, Y.V. Silica-scaled chryso phytes of Lake Baikal. Acta Biol. Sib. 2017, 3, 47–56. [Google Scholar] [CrossRef] [Green Version]
  94. Babich, D.B.; Korotaev, V.N.; Magritsky, D.V.; Mikhailov, V.N. Lower Indigirka: Estuarine and Riverbed Processes; GEOS: Moskow, Russia, 2001; p. 202. [Google Scholar]
  95. Němcová, Y.; Neustupa, J.; Kvíderová, J.; Řezáčová-Škaloudová, M. Morphological plasticity of silica scales of Synura echinulate (Synurophyceae) in crossed gradients of light and temperature—A geometric morphometric approach. Nova Hedwig. Beih. 2010, 136, 21–32. [Google Scholar] [CrossRef]
  96. Němcová, Y.; Pichrtová, M. Shape dynamics of silica scales (Chrysophyceae, Stramenopiles) associated with pH. Fottea Olomouc 2012, 12, 281–291. [Google Scholar] [CrossRef] [Green Version]
  97. Green, R.B. A new species of Spiniferomonas (Chrysophyceae) from an Alberta Lake. Can. J. Bot. 1979, 57, 557–560. [Google Scholar] [CrossRef]
  98. Eloranta, P. Biogeography of chrysophytes in Finnish lakes. In Chrysophyte Algae: Ecology, Phylogeny and Development; Sandgren, C.D., Smol, J.P., Kristiansen, J., Eds.; Cambridge University Press: Cambridge, UK, 1995; pp. 214–231. [Google Scholar] [CrossRef]
  99. Siver, P.A. The distribution of chrysophytes along environmental gradients: Their use as biologicalindicators. In Chrysophyte Algae; Sandgren, C.D., Smol, J.P., Kristiansen, J., Eds.; Cambridge University Press: Cambridge, UK, 1995; pp. 232–268. [Google Scholar]
  100. Kristiansen, J. Golden Algae: A Biology of Chrysophytes; A.R.G. Gantner Verlag: Königstein, Germany, 2005; p. 167. [Google Scholar]
  101. Siver, P.A.; Lott, A.M. The scaled chrysophyte flora in freshwater ponds and lakes from Newfoundland, Canada, and their relationship to environmental variables. Cryptogam. Algol. 2017, 38, 325–347. [Google Scholar] [CrossRef]
  102. Dodds, W.K.; Jones, J.R.; Welch, E.B. Suggested classification for stream trophic state: Distributions of temperate stream types by chlorophyll, total nitrogen and Phosphorus. Water Res. 1998, 32, 1455–1462. [Google Scholar] [CrossRef]
  103. Bessudova, A.Y.; Sorokovikova, L.M.; Sinyukovich, V.N.; Firsova, A.D.; Tomberg, I.V.; Likhoshway, Y.V. Effects of water levels on species diversity of silica-scaled chrysophytes in large tributaries of Lake Baikal. Acta Biol. Sib. 2020, 95, 1–24. [Google Scholar] [CrossRef]
  104. Rojo, C.; Mesquita-Joanes, F.; Monrós, J.S.; Armengol, J.; Sasa, M.; Bonilla, F.; Rueda, R.; Benavent-Corai, J.; Piculo, R.; Seguraet, M.M. Hydrology affects environmental and spatial structuring of microalgal metacommunities in tropical Pacific coast wetlands. PLoS ONE 2016, 11, e0149505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Talbot, C.J.; Bennett, K.M.; Cassell, K.; Hanes, D.M.; Minor Hans, E.C.; Paerl, H.; Raymond, P.A.; Vargas, R.; Vidon, P.G.; Wollheim, W.; et al. The impact of flooding on aquatic ecosystem services. Biogeochemistry 2018, 141, 439–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Řezáčová-Škaloudová, M.; Neustupa, J.; Němcová, Y. Effect of temperature on the variability of silicate structures in Mallomonas kalinae and Synura curtispina (Synurophyceae). Nova Hedwig. Beih. 2010, 136, 55–69. [Google Scholar] [CrossRef]
  107. Kristiansen, J. Cosmopolitan chrysophytes. Syst. Geogr. Plants 2000, 70, 291–300. [Google Scholar] [CrossRef]
Figure 1. Schematic map of the study region and locations of sampling sites. Broken line marks the boundary of the Arctic Circle.
Figure 1. Schematic map of the study region and locations of sampling sites. Broken line marks the boundary of the Arctic Circle.
Sustainability 13 13768 g001
Figure 2. Scatter plot of sampling sites done on results of multidimensional scaling based on the Jaccard distances showing difference by species lists. Arrows designate gradients of quantitative factors reliably affecting the species composition in samples (Table 3). Broken lines show the species richness gradient vector in samples.
Figure 2. Scatter plot of sampling sites done on results of multidimensional scaling based on the Jaccard distances showing difference by species lists. Arrows designate gradients of quantitative factors reliably affecting the species composition in samples (Table 3). Broken lines show the species richness gradient vector in samples.
Sustainability 13 13768 g002
Figure 3. Comparison of species composition of silica-scaled chrysophytes by their grouping in the longitudinal (A) and latitudinal (B) groups (according to [59,60,107] for waters north of the Arctic Circle from 1981 to 2010.
Figure 3. Comparison of species composition of silica-scaled chrysophytes by their grouping in the longitudinal (A) and latitudinal (B) groups (according to [59,60,107] for waters north of the Arctic Circle from 1981 to 2010.
Sustainability 13 13768 g003
Table 1. List of chrysophytes and their distribution in the estuaries. Site numbers correspond to those in Figure 1 and Table 1. Species found in Russia for the first time are in bold.
Table 1. List of chrysophytes and their distribution in the estuaries. Site numbers correspond to those in Figure 1 and Table 1. Species found in Russia for the first time are in bold.
No Site12345678911121420212627303132
Number Species
1.Chrysosphaerella brevispina Korshikov +++ ++ + +
2.Ch. coronacircumspina Wujek & Kristiansen in Wujek, Gretz & Wujek+++++++ +
3.Ch. longispina Lauterborn ++
4.Ch. rotundata Škaloudová & Škaloud +
5.Chrysosphaerella sp. ++ + +
6.Paraphysomonas acuminata acuminata Scoble & Cavalier-Smith +
7.P. bandaiensis Takahashi +
8.P. circumvallata Preisig & Hibberd + ++
9.P. gladiata Preisig & Hibberd ++ ++ +
10.P. punctata Preisig & Hibberd +++
11.P. cf. punctata ssp. simplicior Preisig & Hibberd + +
12.P. vulgaris Scoble & Cavalier-Smith + + +
13.P. uniformis hemiradia Scoble & Cavalier-Smith + +++ + + ++
14.Paraphysomonas sp. 1 +++
15.Paraphysomonas sp. 2 +
16.Paraphysomonas sp. 3 ++
17.Paraphysomonas sp. 4 +
18.Lepidochromonas butcheri (Pennick & Clarke) Kapustin & Guiry +
19.L. coronata (Moestrup & Zimmerman) Kapustin & Guiry +
20.L. diadernifera (Takahashi) Kristiansen +
21.L. eiffelii (Thomsen) Kapustin & Guiry +
22.L. elegantissima (Kling & Kristiansen) Kapustin & Guiry +
23.L. quadrispina (Thomsen & Kristiansen) Kapustin & Guiry +
24.L. quadrispina (Preisig & Hibberd) Kapustin & Guiry +
25.L. undulata (Preisig & Hibberd) Kapustin & Guiry +
26.Lepidochromonas sp. 1 +
27.Lepidochromonas sp. 2 +
28.Spiniferomonas abei Takahashi +
29.S. bourrellyi Takahashi + +++ +
30.S. conica Takahashi +
31.S. cornuta Balonov + ++ +
32.S. minuta Nicholls +
33.S. serrata Nicholls ++++ + + +
34.S. silverensis Nicholls +
35.S. takahashii Nicholls +
36.S. trioralis Takahashi +++++ + +
37.S. trioralis f. cuspidata Balonov+ ++
38.Mallomonas acaroides Perty ++ + +++ ++ +
39.M. actinoloma Takahashi ++
40.M. akrokomos Ruttner ++++ + +
41.M. alpina Pascher & Ruttner+++++++ +++ ++ ++
42.M. alata f. hualvensis Asmund, Cronberg & Dürrschmidt ++
43.M. areolata Nygaard + +
44.M. caudata Iwanoff ++ ++ ++
45.M. cratis Harris & Bradley +
46.M. crassisquama (Asmund) Fott+++++ + ++++ ++ ++
47.M. elongata Reverdin + ++ +
48.M. eoa Takahashi in Asmund & Takahashi +
49.M. kuzminii Gusev & Kulikovskiy ++ + ++ +
50.M. lychenensis Conrad +
51.M. multiunca Asmund +
52.M. papillosa Harris & Bradley +
53.M. parvula Dürrschmidt +
54.M. cf. pumilio Harris & Bradley +
55.M. punctifera Korshikov +++ +++
56.M. striata Asmund +
57.M. tonsurata Teiling+++++++ ++++ ++ +
58.Synura cf. americana Kynclová & Škaloud+++++++++++++ +++++
59.S. borealis Škaloud & Škaloudová + +++++++++ ++++ +
60.S. cf. cornuta Škaloud, Škaloudová & Siver+ +
61.S. conopea Kynclová & Škaloud +++++++++++ + + +
62.S. curtispina (Petersen & Hansen) Asmund + +
63.S. curtispina f. reticulata Asmund +++ ++
64.S. echinulate Korshikov + ++ + ++ + +
65.S. glabra (Korshikov) Škaloud & Kynclová+ + ++ ++ +
66.S. leptorhabda Nicholls in Nicholls & Gerrath +
67.S. macropora Škaloud & Kynclová+ + ++ +
68.S. mammillosa Takahashi + +++ +
69.S. nygaardii (Petersen & Hansen) Kristiansen + + +
70.S. petersenii (Korshikov) Škaloud & Kynclová+++++++ + + + +
71.S. petersenii f. taymyrensis Kristiansen +
72.S. praefracta (Asmund) Škaloud & Škaloudová++++++++ +++ +++++
73.S. cf. soroconopea Jo, Shin, Kim & Siver+ +++ +++++ ++ +
74.S. sphagnicola Korshikov ++ + +
75.S. spinosa Korshikov + + + + +++ +
76.S. spinosa f. longispina Petersen & Hansen ++ + + + +++
77.S. uvella Ehrenberg + +
78.S. cf. vinlandica Škaloud, Škaloudová & Siver +
79.Synura sp. 1 +
80.Synura sp. 2 + ++ + + +
81.Synura sp. 3 +
82.Synura sp. 4 +++ ++ ++
Total14222640522121615233153116167924
Table 2. A listing of the sampling sites from the rivers of Yakutia area included in this study. Site names and physicochemical parameters of the river mouths (site numbers in Figure 1).
Table 2. A listing of the sampling sites from the rivers of Yakutia area included in this study. Site names and physicochemical parameters of the river mouths (site numbers in Figure 1).
Site NumberSampling DateCoordinatesRiverT, °CpH, EpiWater Transparency, mSuspended Matter, mg/LΣ Ions, mg/LSi, mg/LPtotal, μg/LO2, mg/LCO2, mg/LNumber of SpeciesTotal Number of Species in the River
110 August 200870°62′12′′ N,
122°55′654′′ E
Olenyok16.07.894.07.20270.142.44186.009.403.521475
212 August 200871°30′749′′ N,
122°59′67′′ E
Olenyok16.4 22
313 August 200871°78′533′′ N,
123°73′743′′ E
Olenyok16.47.934.011.20228.152.3026.009.193.9626
413 August 200871°78′935′′ N,
123°69′794′′ E
Olenyok17.46.784.06.4067.461.3433.009.812.6440
513 August 200871°81′098′′ N,
123°68′292′′ E
Olenyok16.87.524.07.20180.671.8038.009.264.4052
613 August 200872°13′966′′ N,
123°45′577′′ E
Olenyok17.1 21
714 August 200872°55′621′′ N,
122°04′922′′ E
Olenyok15.07.804.05.60205.522.4826.009.833.9621
817 July 200970°63′703′′ N,
135°17′208′′ E
Yana14.6 640
918 July 200970°75′302′′ N,
136°17′058′′ E
Yana14.57.130.814.80101.001.9040.0010.332.6415
1018 July 200970°74′704′′ N,
136°20′373′′ E
Yana15.16.690.94.8045.961.1020.0010.602.640
1118 July 200970°76′075′′ N,
136°23′337′′ E
Yana14.27.130.98.4092.061.6050.0010.282.6423
1218 July 200971°14′201′′ N,
136°12′924′′ E
Yana147.130.914.00104.261.9050.0010.302.6431
1326 June 200968°33′803′′ N,
146°03′186′′ E
Indigirka14.27.770.279.2097.652.4069.009.863.5207
1426 June 200968°34′601′′ N,
146°01′735′′ E
Indigirka14.97.240.18.0052.212.0033.009.455.725
1526 June 200968°35′042′′ N,
145°99′79′′ E
Indigirka14.37.290.121.6060.912.2033.009.474.400
1627 June 200968°56′854′′ N,
146°74′313′′ E
Indigirka14.2 0
1728 June 200968°85′087′′ N,
147°38′221′′ E
Indigirka14.77.430.110.0090.121.9033.009.613.520
1829 June 200969°19′623′′ N,
147°49′045′′ E
Indigirka15.1 0
1929 June 200969°55′775′′ N,
147°62′09′′ E
Indigirka15.57.550.144.8090.650.5033.009.004.400
2029 June 200969°57′149′′ N,
147°62′26′′ E
Indigirka17.77.090.16.4051.142.5013.009.607.483
2129 June 200969°58′76′′ N,
147°59′92′′ E
Indigirka17.76.990.17.6052.602.5013.009.366.161
2230 June 200970°52′957′′ N,
147°76′007′′ E
Indigirka15.5 0
232 July 200970°51′058′′ N,
147°68′349′′ E
Indigirka12.57.490.1130.4097.222.50100.009.592.640
242 July 200970°53′011′′ N,
147°70′653′′ E
Indigirka14.17.470.188.4080.662.40168.009.423.080
252 July 200970°52′957′′ N,
147°76′007′′ E
Indigirka14.17.440.184.8080.712.50168.009.363.080
265 August 201068°31′505′′ N,
157°74′07′′ E
Kolyma16.7 1631
276 August 201068°71′893′′ N,
158°66′36′′ E
Kolyma15.97.601.95.20126.061.3480.009.845.2816
286 August 201068°70′169′′ N,
158°70′129′′ E
Kolyma12.07.643.05.2068.351.6440.0010.523.960
296 August 201068°70′217′′ N,
158°72′184′′ E
Kolyma14.47.481.95.2092.002.04380.0010.633.520
306 August 201068°57′6098′′ N 159°60′2341′′ EKolyma 7
317 August 201068°51′189′′ N,
160° 88′912′′ E
Kolyma13.97.631.45.60124.011.34120.0010.305.289
327 August 201068°46′522′′ N,
160°80′042′′ E
Kolyma14.77.552.56.0094.491.44180.0010.603.9624
337 August 201068°50′719′′ N,
160°97′508′′ E
Kolyma14.47.261.45.6090.321.44180.009.893.960
Table 3. PERMANOVA analysis results on reliability of impact of abiotic environmental parameters on species richness in samples. The reliable impact of a factor is seen at p value ≤ 0.05. All reliably affecting factors are in bold.
Table 3. PERMANOVA analysis results on reliability of impact of abiotic environmental parameters on species richness in samples. The reliable impact of a factor is seen at p value ≤ 0.05. All reliably affecting factors are in bold.
Factor NamesR2 Covariance Coefficientp Value
Sampling river0.2900.001
Distance to the mouth, km 0.0980.049
T, °C0.0900.042
Water transparency, m0.0970.031
pH, Epi0.0650.275
Suspended matter, mg/L0.0620.373
Color, Pt-Co0.1050.019
O2, mg/L0.0680.228
Oxygen saturation, %0.0680.232
CO2, mg/L0.0970.048
Σions, mg/L0.0650.300
Water hardness, mg-equ/L0.0690.249
Ca2+, mg/L0.0560.505
Mg2+, mg/L0.0990.025
Na+, mg/L0.0520.564
K+, mg/L0.0600.415
HCO3, mg/L0.0630.342
Cl, mg/L0.0450.727
SO42−, mg/L0.0620.364
NH4, mg/L0.0750.204
NO2, mg/L0.0620.199
NO3, mg/L0.0770.138
PO4, mg/L0.0710.205
Ptotal, μg/L0.0550.515
Si, mg/L0.0590.419
COD, mg/L0.0740.177
BOD5, mg/L0.0630.328
Fetotal, mg/L0.1020.025
Phenols, mg/L0.0660.298
Surfactants, mg/L0.1330.007
Petrochemicals, mg/L0.1210.009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bessudova, A.; Gabyshev, V.; Firsova, A.; Gabysheva, O.; Bukin, Y.; Likhoshway, Y. Diversity Variation of Silica-Scaled Chrysophytes Related to Differences in Physicochemical Variables in Estuaries of Rivers in an Arctic Watershed. Sustainability 2021, 13, 13768. https://doi.org/10.3390/su132413768

AMA Style

Bessudova A, Gabyshev V, Firsova A, Gabysheva O, Bukin Y, Likhoshway Y. Diversity Variation of Silica-Scaled Chrysophytes Related to Differences in Physicochemical Variables in Estuaries of Rivers in an Arctic Watershed. Sustainability. 2021; 13(24):13768. https://doi.org/10.3390/su132413768

Chicago/Turabian Style

Bessudova, Anna, Viktor Gabyshev, Alena Firsova, Olga Gabysheva, Yurij Bukin, and Yelena Likhoshway. 2021. "Diversity Variation of Silica-Scaled Chrysophytes Related to Differences in Physicochemical Variables in Estuaries of Rivers in an Arctic Watershed" Sustainability 13, no. 24: 13768. https://doi.org/10.3390/su132413768

APA Style

Bessudova, A., Gabyshev, V., Firsova, A., Gabysheva, O., Bukin, Y., & Likhoshway, Y. (2021). Diversity Variation of Silica-Scaled Chrysophytes Related to Differences in Physicochemical Variables in Estuaries of Rivers in an Arctic Watershed. Sustainability, 13(24), 13768. https://doi.org/10.3390/su132413768

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