Environmental drivers of plankton protist communities along latitudinal and vertical gradients in the oldest and deepest freshwater lake

Characterization of bacterial community dynamics dominated by salinity in lakes of the Inner Mongolian Plateau, China

Microorganisms in lakes are sensitive to salinity fluctuations. Despite extensive prior research on bacterial communities, our understanding of their characteristics and assembly mechanisms in lakes, especially in desert lakes with different salinities. To address this issue, we collected three samples from freshwater lakes, six from brackish lakes, and five from salt lakes in the Badanjilin Desert. The 16S rRNA gene sequencing was applied to investigate the bacterial interactions with rising salinity, community coexistence patterns, and assembly mechanisms. Our findings suggested that the increased lake salinity significantly reduces the bacterial community diversity and enhanced the community differentiation. Significant variations were observed in the contribution of biomarkers from Cyanobacteria, Chloroflexi, and Halobacterota to the composition of the lake bacterial communities. The bacterial communities in the salt lakes exhibited a higher susceptibility to salinity limitations than those in the freshwater and brackish lakes. In addition, the null modeling analyses confirmed the quantitative biases in the stochastic assembly processes of bacterial communities across freshwater, brackish, and saline lakes. With the increasing lake salinity, the significance of undominated and diffusion limitation decreased slightly, and the influence of homogenizing dispersal on community assembly increased. However, the stochasticity remained the dominant process across all lakes in the Badanjilin Desert. The analysis of co-occurring networks revealed that the rising salinity reduced the complexity of bacterial network structures and altered the interspecific interactions, resulting in the increased interspecies collaboration with increasing salinity levels. Under the influence of salinity stress, the key taxon Cyanobacteria in freshwater lakes (Schizothrix_LEGE_07164) was replaced by Proteobacteria (Thalassobaculum and Polycyclovorans) in brackish lakes, and Thermotogota (SC103) in salt lakes. The results indicated the symbiotic patterns of bacterial communities across varying salinity gradients in lakes and offer insights into potential mechanisms of community aggregation, thereby enhancing our understanding of bacterial distribution in response to salinity changes.

Metagenomic 18S rDNA reads revealed zonation of eukaryotic communities in the Yongle blue hole

The Yongle blue hole (YBH), situated in the South China Sea, represents a compelling subject of study in marine microbiology due to its unique redox-layered microbial ecosystems. However, the diversity and ecology of microbial eukaryotes within the YBH remains underexplored. This study endeavors to bridge this gap through the application of the in situ microbial filtration and fixation (ISMIFF) device to collect 0.22-30 μm microbial samples from 21 water layers of YBH. Subsequent extraction of 18S rRNA metagenomic reads of 21 metagenomes and 10 metatranscriptomes facilitated a comprehensive analysis of community structures. Findings revealed a pronounced superiority in the diversity and richness of eukaryotic microorganisms in the oxic zone compared to its suboxic and anoxic counterparts. Notably, Dinophyceae and Maxillopoda emerged as the predominant taxa based on the analysis of the 18S rRNA reads for the V4 and V9 regions, which showed stratification In their relative abundance and suggested their potential role in the thermo-halocline boundaries and oxic-anoxic interface. Specifically, In these eukaryotic microbial communities, Dinophyceae exhibited significant abundance at 20 m (20.01%) and 105 m (26.13%) water depths, while Maxillopoda was prevalent at 40 m (22.84%), 80 m (23.19%), and 100 m (15.42%) depths. A part of these organisms, identified as larvae and protists, were likely attracted by swarming chemosynthetic bacterial prey prevailing at the thermo-halocline boundaries and oxic-anoxic interface. Furthermore, the phylogenetic relationships of the major 18S operational taxonomic units (OTUs) showed a close adjacency to known species, except for three Dinophyceae OTUs. In conclusion, this study provides critical insights into the vertical distribution and transcriptional activity of <30-μm eukaryotic microbes, shedding light on the taxonomic novelty of eukaryotic microorganisms within the semi-enclosed blue holes.

Macrophyte coverage drives microbial community structure and interactions in a shallow sub-tropical lake

Shallow lakes are typically dominated by macrophytes, which have important functional roles regulating trophic conditions and creating biological habitat. Macrophytes have been shown to strongly influence water chemistry and shape microbial communities in shallow lakes. In Florida, many large, shallow lakes are dominated by alien invasive, submersed macrophytes, such as hydrilla (Hydrilla verticillata [L.F.] Royle) and are intensively managed to reduce infestations and contain the spread of these alien invasive macrophytes. In this study, we investigated the effects of large (40 ha) herbicidal and mechanical control treatments on a large lake located in Central Florida that resulted in the reduction of Hydrilla and concomitant changes in water chemistry and microbial communities (both bacteria and protists [microbial eukaryotes]). We observed a considerable decrease in macrophyte coverage associated with plant control treatments as well as a temporal change in macrophyte coverage in Lake Tohopekaliga. We found that changes in macrophyte coverage, regardless of treatment type, significantly affected the water chemistry of the lake, resulting in a sharp increase of chlorophyll a concentration as well as an increase in turbidity with the decrease of macrophyte coverage. Moreover, the decline in macrophytes led to decreases in microbial community diversity with over-representation of phototrophic functional groups. Specifically, we observed an increase in cyanobacteria with the decrease in macrophyte coverage. Our study highlights the advantages and disadvantages of macrophyte control. Although there was an initial decrease in macrophyte coverage associated with the chemical and mechanical control of aquatic plants, after a few months, we found a considerable increase in coverage. In addition, the increase of cyanobacterial relative abundance demonstrates the possible consequences of aquatic plant control such as cyanobacterial blooms if there is a continued decline of macrophytes.

Spatiotemporal assembly and functional composition of planktonic microeukaryotic communities along productivity gradients in a subtropical lake

Microeukaryotes play crucial roles in the microbial loop of freshwater ecosystems, functioning both as primary producers and bacterivorous consumers. However, understanding the assembly of microeukaryotic communities and their functional composition in freshwater lake ecosystems across diverse environmental gradients remains limited. Here, we utilized amplicon sequencing of 18S rRNA gene and multivariate statistical analyses to examine the spatiotemporal and biogeographical patterns of microeukaryotes in water columns (at depths of 0.5, 5, and 10 m) within a subtropical lake in eastern China, covering a 40 km distance during spring and autumn of 2022. Our results revealed that complex and diverse microeukaryotic communities were dominated by Chlorophyta (mainly Chlorophyceae), Fungi, Alveolata, Stramenopiles, and Cryptophyta lineages. Species richness was higher in autumn than in spring, forming significant hump-shaped relationships with chlorophyll a concentration (Chl-a, an indicator of phytoplankton biomass). Microeukaryotic communities exhibited significant seasonality and distance-decay patterns. By contrast, the effect of vertical depth was negligible. Stochastic processes mainly influenced the assembly of microeukaryotic communities, explaining 63, 67, and 55% of community variation for spring, autumn, and both seasons combined, respectively. Trait-based functional analysis revealed the prevalence of heterotrophic and phototrophic microeukaryotic plankton with a trade-off along N:P ratio, Chl-a, and dissolved oxygen (DO) gradients. Similarly, the mixotrophic proportions were significantly and positively correlated with Chl-a and DO concentrations. Overall, our findings may provide useful insights into the assembly patterns of microeukaryotes in lake ecosystem and how their functions respond to environmental changes.

Lakes-scale pattern of eukaryotic phytoplankton diversity and assembly process shaped by electrical conductivity in central Qinghai-Tibet Plateau

Phytoplankton are the main primary producers in aquatic ecosystems and play an important role in food web and geochemical cycles. Its diversity, community structure, and assembly process are influenced by several factors. Alpine lake ecosystems are relatively weak and extremely sensitive to global climate change. However, the impact of climate change on phytoplankton in Qinghai-Tibet Plateau lakes and their responses are still unclear. In this study, we analyzed the diversity, environmental drivers, and assembly process of phytoplankton community in the central QTP lakes. The phytoplankton of these lakes can be primarily distinguished into freshwater and brackish types, with significant differences in species diversity and community dissimilarity. Both shared nearly same key environmental factors that significantly affecting phytoplankton such as EC, and brackish lakes were also positively correlative with TN. Stochastic process was predominant in phytoplankton assembly. Additionally, freshwater and brackish lakes were dominated by dispersal limitation and heterogeneous selection respectively. Alpine lakes had significant EC thresholds, and their diversity and assembly processes changed significantly around the thresholds. The present findings have important implications for understanding and predicting the response of lake phytoplankton communities to climate change and for making decisions to protect the ecological resources of alpine lakes.

Genomic Tools in Biological Invasions: Current State and Future Frontiers

Human activities are accelerating rates of biological invasions and climate-driven range expansions globally, yet we understand little of how genomic processes facilitate the invasion process. Although most of the literature has focused on underlying phenotypic correlates of invasiveness, advances in genomic technologies are showing a strong link between genomic variation and invasion success. Here, we consider the ability of genomic tools and technologies to (i) inform mechanistic understanding of biological invasions and (ii) solve real-world issues in predicting and managing biological invasions. For both, we examine the current state of the field and discuss how genomics can be leveraged in the future. In addition, we make recommendations pertinent to broader research issues, such as data sovereignty, metadata standards, collaboration, and science communication best practices that will require concerted efforts from the global invasion genomics community.

Molecular and morphological characterization of four new ancyromonad genera and proposal for an updated taxonomy of the Ancyromonadida

Ancyromonads are small biflagellated protists with a bean-shaped morphology. They are cosmopolitan in marine, freshwater, and soil environments, where they attach to surfaces while feeding on bacteria. These poorly known grazers stand out by their uncertain phylogenetic position in the tree of eukaryotes, forming a deep-branching "orphan" lineage that is considered key to a better understanding of the early evolution of eukaryotes. Despite their ecological and evolutionary interest, only limited knowledge exists about their true diversity. Here, we aimed to characterize ancyromonads better by integrating environmental surveys with behavioral observation and description of cell morphology, for which sample isolation and culturing are indispensable. We studied 18 ancyromonad strains, including 14 new isolates and seven new species. We described three new and genetically divergent genera: Caraotamonas, Nyramonas, and Olneymonas, together encompassing four species. The remaining three new species belong to the already-known genera Fabomonas and Ancyromonas. We also raised Striomonas, formerly a subgenus of Nutomonas, to full genus status, on morphological and phylogenetic grounds. We studied the morphology of diverse ancyromonads under light and electron microscopy and carried out molecular phylogenetic analyses, also including 18S rRNA gene sequences from several environmental surveys. Based on these analyses, we have updated the taxonomy of Ancyromonadida.

The effect of metabarcoding 18S rRNA region choice on diversity of microeukaryotes including phytoplankton

Metabarcoding using high throughput sequencing of amplicons of the 18S rRNA gene is one of the widely used methods for assessing the diversity of microeukaryotes in various ecosystems. We investigated the effectiveness of the V4 and V8-V9 regions of the 18S rRNA gene by comparing the results of metabarcoding microeukaryotic communities using the DADA2 (ASV), USEARCH-UNOISE3 (ZOTU), and USEARCH-UPARSE (OTU with 97% similarity) algorithms. Both regions showed similar levels of genetic variability and taxa identification accuracy. Richness for DADA2 datasets of both regions was lower than for UNOISE3 and UPARSE datasets, which is due to more accurate error correction in amplicons. Microeukaryotic communities (autotrophs and heterotrophs) structure identified using both regions showed a significant relationship with phytoplankton (autotrophs) communities structure based on microscopy in a seasonal freshwater sample series. The strongest relationship was found between the phytoplankton species and V8-V9 ASVs produced by DADA2.

Mechanistic insights into dissolved organic matter-driven protistan and bacterial community dynamics influenced by vegetation restoration

Vegetation restoration projects can not only improve water quality by absorbing and transferring pollutants and nutrients from non-vegetation sources, but also protect biodiversity by providing habitat for biological growth. However, the mechanism of the protistan and bacterial assembly processes in the vegetation restoration project were rarely explored. To address this, based on 18 S rRNA and 16 S rRNA high-throughput sequencing, we investigated the mechanism of protistan and bacterial community assembly processes, environmental conditions, and microbial interactions in the rivers with (out) vegetation restoration. The results indicated that the deterministic process dominated the protistan and bacterial community assembly (94.29% and 92.38%), influenced by biotic and abiotic factors. For biotic factors, microbial network connectivity was higher in the vegetation zone (average degree = 20.34) than in the bare zone (average degree = 11.00). For abiotic factors, the concentration of dissolved organic carbon ([DOC]) was the most important environmental factor affecting the microbial community composition. [DOC] was lower significantly in vegetation zone (18.65 ± 6.34 mg/L) than in the bare zone (28.22 ± 4.82 mg/L). In overlying water, vegetation restoration upregulated the protein-like fluorescence components (C1 and C2) by 1.26 and 1.01-folds and downregulated the terrestrial humic-like fluorescence components (C3 and C4) by 0.54 and 0.55-folds, respectively. The different DOM components guided bacteria and protists to select different interactive relationships. The protein-like DOM components led to bacterial competition, whereas the humus-like DOM components resulted in protistan competition. Finally, the structural equation model was established to explain that DOM components can affect protistan and bacterial diversity by providing substrates, facilitating microbial interactions, and promoting nutrient input. In general, our study provides insights into the responses of vegetation restored ecosystems to the dynamics and interactives in the anthropogenically influenced river and evaluates the ecological restoration performance of vegetation restoration from a molecular biology perspective.

Local-Scale Damming Impact on the Planktonic Bacterial and Eukaryotic Assemblages in the upper Yangtze River

Dam construction and impoundment cause discontinuities in the natural biophysical gradients in rivers. These discontinuities may alter distinctive habitats and different microbial community assembly mechanisms upstream and downstream of dams, which reflect the potential impacts of damming on riverine aquatic ecosystems. In this study, we investigated the planktonic microbial assemblages of three large dams in the upper Yangtze River by using high-throughput sequencing. The results revealed that the alpha diversity indexes increased downstream of the dams. In addition, more eukaryotic ASVs solely occurred downstream of the dams, which indicated that a large proportion of eukaryotes appeared downstream of the dams. The nonmetric multidimensional scaling analysis indicated that there was no obvious geographic clustering of the planktonic microbial assemblages among the different locations or among the different dams. However, the dam barriers changed dam-related variables (maximum dam height and water level) and local environmental variables (water temperature, DOC, etc.) that could possibly affect the assembly of the planktonic microbial communities that are closest to the dams. A co-occurrence network analysis demonstrated that the keystone taxa of the planktonic bacteria and eukaryotes decreased downstream of the dams. In particular, the keystone taxa of the eukaryotes disappeared downstream of the dams. The robustness analysis indicated that the natural connectivity of the microbial networks decreased more rapidly upstream of the dams, and the downstream eukaryotic network was more stable. In conclusion, damming has a greater impact on planktonic eukaryotes than on bacteria in near-dam areas, and planktonic microbial assemblages were more susceptible to the environmental changes. Our study provides a better understanding of the ecological effects of river damming.

Revealing the Differences in Ulnaria acus and Fragilaria radians Distribution in Lake Baikal via Analysis of Existing Metabarcoding Data

Two diatom species, Ulnaria acus and Fragilaria radians, are morphologically very similar and often coexist, which makes it difficult to compare their abundances. However, they are easily separated by molecular data; thus, in this work, we attempted to estimate the differences in their spatial and temporal distribution from existing metabarcoding datasets. Reanalyzing published sequences with an ASV-based pipeline and ad hoc classification routine allowed us to estimate the relative abundances of the two species, increasing the precision compared to usual OTU-based analyses. Existing data permit qualitative comparisons between two species that cannot be differentiated by other methods, detecting the distinct seasonal peaks and spatial distributions of F. radians and U. acus.

Does filter pore size introduce bias in DNA sequence-based plankton community studies?

The cell size of microbial eukaryotic plankton normally ranges from 0.2 to 200 μm. During the past decade, high-throughput sequencing of DNA has been revolutionizing their study on an unprecedented scale. Nonetheless, it is currently unclear whether we can accurately, effectively, and quantitatively depict the microbial eukaryotic plankton community using size-fractionated filtration combined with environmental DNA (eDNA) molecular methods. Here we assessed the microbial eukaryotic plankton communities with two filtering strategies from two subtropical reservoirs, that is one-step filtration (0.2–200 μm) and size-fractionated filtration (0.2–3 and 3–200 μm). The difference of 18S rRNA gene copy abundance between the two filtering treatments was less than 50% of the 0.2–200 μm microbial eukaryotic community for 95% of the total samples. Although the microbial eukaryotic plankton communities within the 0.2–200 μm and the 0.2–3 and 3–200 μm size fractions had approximately identical 18S rRNA gene copies, there were significant differences in their community composition. Furthermore, our results demonstrate that the systemic bias introduced by size-fractionation filtration has more influence on unique OTUs than shared OTUs, and the significant differences in abundance between the two eukaryotic plankton communities largely occurred in low-abundance OTUs in specific seasons. This work provides new insights into the use of size-fractionation in molecular studies of microbial eukaryotes populating the plankton.

Seasonal Succession and Coherence Among Bacteria and Microeukaryotes in Lake Baikal

Microorganisms exhibit seasonal succession governed by physicochemical factors and interspecies interactions, yet drivers of this process in different environments remain to be determined. We used high-throughput sequencing of 16S rRNA and 18S rRNA genes to study seasonal dynamics of bacterial and microeukaryotic communities at pelagic site of Lake Baikal from spring (under-ice, mixing) to autumn (direct stratification). The microbial community was subdivided into distinctive coherent clusters of operational taxonomic units (OTUs). Individual OTUs were consistently replaced during different seasonal events. The coherent clusters change their contribution to the microbial community depending on season. Changes of temperature, concentrations of silicon, and nitrates are the key factors affected the structure of microbial communities. Functional prediction revealed that some bacterial or eukaryotic taxa that switched with seasons had similar functional properties, which demonstrate their functional redundancy. We have also detected specific functional properties in different coherent clusters of bacteria or microeukaryotes, which can indicate their ability to adapt to seasonal changes of environment. Our results revealed a relationship between seasonal succession, coherency, and functional features of freshwater bacteria and microeukaryotes.

Diversity and ecology of protists revealed by metabarcoding

Protists are the dominant eukaryotes in the biosphere where they play key functional roles. While protists have been studied for over a century, it is the high-throughput sequencing of molecular markers from environmental samples - the approach of metabarcoding - that has revealed just how diverse, and abundant, these small organisms are. Metabarcoding is now routine to survey environmental diversity, so data have rapidly accumulated from a multitude of environments and at different sampling scales. This mass of data has provided unprecedented opportunities to study the taxonomic and functional diversity of protists, and how this diversity is organised in space and time. Here, we use metabarcoding as a common thread to discuss the state of knowledge in protist diversity research, from technical considerations of the approach to important insights gained on diversity patterns and the processes that might have structured this diversity. In addition to these insights, we conclude that metabarcoding is on the verge of an exciting added dimension thanks to the maturation of high-throughput long-read sequencing, so that a robust eco-evolutionary framework of protist diversity is within reach.

Microbial eukaryotes assemblages and potential novel diversity in four tropical East-African Great Lakes

East-African Great Lakes are old and unique natural resources heavily utilized by their bordering countries. In those lakes, ecosystem functioning is dominated by pelagic processes, where microorganisms are key components, however protistan diversity is barely known. We investigated the community composition of small eukaryotes (< 10 µm) in surface waters of four African Lakes (Kivu, Edward, Albert and Victoria) by sequencing the 18S rRNA gene. Moreover, in the meromictic Lake Kivu, two stations were vertically studied. We found high protistan diversity distributed in 779 operational taxonomic units (OTUs), spanning in eleven high-rank lineages, being Alveolata (31%), Opisthokonta (20%) and Stramenopiles (17%) the most represented supergroups. Surface protistan assemblage were associated to conductivity and productivity gradients; whereas depth, had a strong effect on protistan community in Kivu, with higher contribution of heterotrophic organisms. Approximately 40% of OTUs had low similarity (< 90%) with reported sequences in public databases, these were mostly coming from deep anoxic waters of Kivu, suggesting a high extent of novel diversity. We also detected several taxa so far considered exclusive of marine ecosystems. Our results unveiled a complex and largely undescribed protistan community, in which several lineages have adapted to different niches after crossing the salinity boundary.

Changes in protist communities in drainages across the Pearl River Delta under anthropogenic influence

Drainages in the Pearl River Delta urban agglomeration (PRDUA) host vital aquatic ecosystems and face enormous pressures from human activities in one of the largest urban agglomerations in the world. Despite being crucial components of aquatic ecosystems, the interactions and assembly processes of the protistan community are rarely explored in areas with serious anthropogenic disturbance. To elucidate the mechanisms of these processes, we used environmental DNA sequencing of 18S rDNA to investigate the influence of environmental factors and species interactions on the protistan community and its assembly in drainages of the PRDUA during summer. The protistan community showed a high level of diversity and a marked spatial pattern in this region. Community assembly was driven primarily by stochastic processes based on the Sloan neutral community model, explaining 74.28%, 75.82%, 73.67%, 74.40% and 51.24% of community variations in the BJ (Beijiang), XJ (Xijiang), PRD (Pearl River Delta), PRE (Pearl River Estuary) areas and in total, respectively. Meanwhile, environmental variables including temperature, pH, dissolved oxygen, transparency, nutrients and land use were strongly correlated with the composition and assembly of the protistan community, explaining 40.40% of variation in the protistan community. Furthermore, the bacterial community was simultaneously analysed by the 16S rDNA sequencing. Co-occurrence network analysis revealed that species interactions within bacteria (81.41% positive) or protists (82.80% positive), and those between bacteria and protists (50% positive and 50% negative) impacted the protistan community assembly. In summary, stochastic processes dominated, whereas species interactions and environmental factors also played important roles in shaping the protistan communities in drainages across the PRDUA. This study provides insights into the ecological patterns, assembly processes and species interactions underlying protistan dynamics in urban aquatic ecosystems experiencing serious anthropogenic disturbance.

Marine signature taxa and core microbial community stability along latitudinal and vertical gradients in sediments of the deepest freshwater lake

Lake Baikal is the deepest (~1.6 km) and most voluminous freshwater reservoir on Earth. Compared to plankton, its benthos remains poorly explored. Here, we ask whether latitude and/or depth determine benthic microbial community structure and how Baikal communities compare to those of other freshwater, brackish and marine sediments. To answer, we collected sediment upper layers (0-1 cm) across a ~600 km North-South transect covering the three basins of the lake and from littoral to bathybenthic depths (0.5-1450  m). Analysis of 16S and 18S rRNA gene amplicon sequences revealed communities with high richness and evenness where rare operational taxonomic units (OTUs) collectively dominated. Archaea represented up to 25% or prokaryotic sequences. Baikal sediments harbored typically marine eukaryotic and prokaryotic OTUs recently identified in some lakes (diplonemids, Bolidophyceae, Mamiellales, SAR202, marine-like Synechococcus, Pelagibacterales) but also SAR324, Syndiniales and Radiolaria. We hypothesize that, beyond the salinity barrier, adaptation to oligotrophy explains the presence of these otherwise typically marine lineages. Baikal core benthic communities were relatively stable across sites and seemed not determined by depth or latitude. Comparative analyses with other freshwater, brackish and marine prokaryotic sediment communities confirmed the distinctness of Baikal benthos, which include elements of similarity to marine and hydrothermally influenced systems.