Anthropogenic impact on the historical phytoplankton community of Lake Constance reconstructed by multimarker analysis of sediment-core environmental DNA

Sediment DNA Records the Critical Transition of Bacterial Communities in the Arid Lake

It is necessary to predict the critical transition of lake ecosystems due to their abrupt, non-linear effects on social-economic systems. Given the promising application of paleolimnological archives to tracking the historical changes of lake ecosystems, it is speculated that they can also record the lake's critical transition. We studied Lake Dali-Nor in the arid region of Inner Mongolia because of the profound shrinking the lake experienced between the 1300 s and the 1600 s. We reconstructed the succession of bacterial communities from a 140-cm-long sediment core at 4-cm intervals and detected the critical transition. Our results showed that the historical trajectory of bacterial communities from the 1200 s to the 2010s was divided into two alternative states: state1 from 1200 to 1300 s and state2 from 1400 to 2010s. Furthermore, in the late 1300 s, the appearance of a tipping point and critical slowing down implied the existence of a critical transition. By using a multi-decadal time series from the sedimentary core, with general Lotka-Volterra model simulations, local stability analysis found that bacterial communities were the most unstable as they approached the critical transition, suggesting that the collapse of stability triggers the community shift from an equilibrium state to another state. Furthermore, the most unstable community harbored the strongest antagonistic and mutualistic interactions, which may imply the detrimental role of interaction strength on community stability. Collectively, our study showed that sediment DNA can be used to detect the critical transition of lake ecosystems.

Relation between the relative abundance and collapse of Aphanizomenon flos-aquae and microbial antagonism in Upper Klamath Lake, Oregon

Aphanizomenon flos-aquae (AFA) is the dominant filamentous cyanobacterium that develops into blooms in Upper Klamath Lake, Oregon, each year. During AFA bloom and collapse, ecosystem conditions for endangered Lost River and shortnose suckers deteriorate, thus motivating the need to identify processes that limit AFA abundance and decline. Here, we investigate the relations between AFA and other members of the microbial community (photosynthetic and nonphotosynthetic bacteria and archaea), how those relations impact abundance and collapse of AFA, and the types of microbial conditions that suppress AFA. We found significant spatial variation in AFA relative abundance during the 2016 bloom period using 16S rRNA sequencing. The Pelican Marina site had the lowest AFA relative abundance, and this was coincident with increased relative abundance of Candidatus Sericytochromatia, Flavobacterium, and Rheinheimera, some of which are known AFA antagonists. The AFA collapse coincided with phosphorus limitation relative to nitrogen and the increased relative abundance of Cyanobium and Candidatus Sericytochromatia, which outcompete AFA when dissolved inorganic nitrogen is available. The data collected in this study indicate the importance of dissolved inorganic nitrogen combined with microbial community structure in suppressing AFA abundance.

Characterization of lacustrine harmful algal blooms using multiple biomarkers: Historical processes, driving synergy, and ecological shifts

Harmful algal blooms (HABs) producing toxic metabolites are increasingly threatening environmental and human health worldwide. Unfortunately, long-term process and mechanism triggering HABs remain largely unclear due to the scarcity of temporal monitoring. Retrospective analysis of sedimentary biomarkers using up-to-date chromatography and mass spectrometry techniques provide a potential means to reconstruct the past occurrence of HABs. By combining aliphatic hydrocarbons, photosynthetic pigments, and cyanotoxins, we quantified herein century-long changes in abundance, composition, and variability of phototrophs, particularly toxigenic algal blooms, in China's third largest freshwater Lake Taihu. Our multi-proxy limnological reconstruction revealed an abrupt ecological shift in the 1980s characterized by elevated primary production, Microcystis-dominated cyanobacterial blooms, and exponential microcystin production, in response to nutrient enrichment, climate change, and trophic cascades. The empirical results from ordination analysis and generalized additive models support climate warming and eutrophication synergy through nutrient recycling and their feedback through buoyant cyanobacterial proliferation, which sustain bloom-forming potential and further promote the occurrence of increasingly-toxic cyanotoxins (e.g., microcystin-LR) in Lake Taihu. Moreover, temporal variability of the lake ecosystem quantified using variance and rate of change metrics rose continuously after state change, indicating increased ecological vulnerability and declined resilience following blooms and warming. With the persistent legacy effects of lake eutrophication, nutrient reduction efforts mitigating toxic HABs probably be overwhelmed by climate change effects, emphasizing the need for more aggressive and integrated environmental strategies.

Environmental DNA metabarcoding revealed the impacts of anthropogenic activities on phytoplankton diversity in Dianchi Lake and its three inflow rivers

Phytoplankton diversity is closely related to environmental variables and has been widely used in ecological health assessment of rivers and lakes. Combining advantages of DNA-based identification and high-throughput sequencing technology, environmental DNA (eDNA) metabarcoding permits a new measurement for biodiversity monitoring in aquatic ecosystems. However, it had rarely been used to explore the variability and similarity of phytoplankton diversity between lake and its inflow rivers and the effects of environmental variables on phytoplankton. This study utilized eDNA metabarcoding to investigate the spatial distribution of phytoplankton and the impacts of environmental variables on the phytoplankton diversity in Dianchi Lake (one of the most polluted urban lakes in China) and its main inflow rivers (Panlong River, Baoxiang River, and Chai River). A total of 243 distinct phytoplankton taxa were detected, covering 9 phyla, 30 classes, 84 orders, and 132 families, and the taxonomic richness of rivers was higher than that of Dianchi Lake. Distinct biodiversity patterns (e.g., community structure, dominant taxon, ɑ-diversity) were exhibited among Dianchi Lake and its three inflow rivers, but similar biodiversity patterns were also observed in Dianchi Lake and the estuarine sites. The patterns of phytoplankton diversity were closely related to environmental variables, which were associated with pollution sources from different anthropogenic activities (e.g., urbanization, water diversion, industrial and agricultural activities). The primary environmental variables correlated with phytoplankton diversity varied in different habitats. The total phosphorus (TP) and chemical oxygen demand (COD) positively correlated with the phytoplankton community structures in Dianchi Lake, whereas negatively correlated in Panlong River and Baoxiang River. The total nitrogen (TN) positively correlated with the phytoplankton community structures in Baoxiang River and Chai River but negatively correlated in Dianchi Lake. Overall, this study provides insights on the phytoplankton diversity monitoring and the conservation of its diversity and healthy management of Dianchi Lake.

Paleoreconstructions of ciliate communities reveal long-term ecological changes in temperate lakes

Ciliates are unicellular heterotrophic organisms that play a key role in aquatic planktonic and benthic food webs. Advances in sedimentary DNA (sed-DNA) analysis offer the possibility to integrate these bioindicators in paleoenvironmental reconstructions. In this study, we used the top–bottom paleolimnological approach and metabarcoding techniques applied to sed-DNA to compare the recent and past (i.e. prior to major anthropogenic impacts) ciliate communities of 48 lakes located along an elevation gradient. Our results show an overall decline in the β-diversity in recent time, especially in lowland lakes, which are more strongly exposed to local human pressures. Analyses of the functional groups indicate important restructuration of the food web, including the recent increase in mixotrophs. Moreover, changes in the benthic ciliates were consistent with the widespread increase in deep water anoxia. Our results provided evidence that sed-DNA can uncover information about past ciliate communities on a wide variety of lakes. Overall, our study demonstrates the potential of using ciliates as new paleoindicators, integrating information from the pelagic to the benthic zones, and providing valuable insights into ecosystem functioning through a trait-based functional community approach. As paleoindicator, they thus offer a more holistic view on the long-term changes of aquatic ecosystems.

Century-Long Homogenization of Algal Communities Is Accelerated by Nutrient Enrichment and Climate Warming in Lakes and Reservoirs of the North Temperate Zone

Anthropogenic pressures can threaten lake and reservoir ecosystems, leading to harmful algal blooms that have become globally widespread. However, patterns of phytoplankton diversity change and community assembly over long-term scales remain unknown. Here, we explore biodiversity patterns in eukaryotic algal (EA) and cyanobacterial (CYA) communities over a century by sequencing DNA preserved in the sediment cores of seven lakes and reservoirs in the North Temperate Zone. Comparisons within lakes revealed temporal algal community homogenization in mesotrophic lakes, eutrophic lakes, and reservoirs over the last century but no systematic losses of α-diversity. Temporal homogenization of EA and CYA communities continued into the modern day probably due to time-lags related to historical legacies, even if lakes go through a eutrophication phase followed by a reoligotrophication phase. Further, algal community assembly in lakes and reservoirs was mediated by both deterministic and stochastic processes, while homogeneous selection played a relatively important role in recent decades due to intensified anthropogenic activities and climate warming. Overall, these results expand our understanding of global change effects on algal community diversity and succession in lakes and reservoirs that exhibit different successional trajectories while also providing a baseline framework to assess their potential responses to future environmental change.

Shifts in DNA yield and biological community composition in stored sediment: implications for paleogenomic studies

Lake sediments hold a wealth of information from past environments that is highly valuable for paleolimnological reconstructions. These studies increasingly apply modern molecular tools targeting sedimentary DNA (sedDNA). However, sediment core sampling can be logistically difficult, making immediate subsampling for sedDNA challenging. Sediment cores are often refrigerated (4 °C) for weeks or months before subsampling. We investigated the impact of storage time on changes in DNA (purified or as cell lysate) concentrations and shifts in biological communities following storage of lake surface sediment at 4 °C for up to 24 weeks. Sediment samples (~ 0.22 g, in triplicate per time point) were spiked with purified DNA (100 or 200 ng) or lysate from a brackish water cyanobacterium that produces the cyanotoxin nodularin or non-spiked. Samples were analysed every 1–4 weeks over a 24-week period. Droplet digital PCR showed no significant decrease in the target gene (nodularin synthetase – subunit F; ndaF) over the 24-week period for samples spiked with purified DNA, while copy number decreased by more than half in cell lysate-spiked samples. There was significant change over time in bacteria and eukaryotic community composition assessed using metabarcoding. Amongst bacteria, the cyanobacterial signal became negligible after 5 weeks while Proteobacteria increased. In the eukaryotic community, Cercozoa became dominant after 6 weeks. These data demonstrate that DNA yields and community composition data shift significantly when sediments are stored chilled for more than 5 weeks. This highlights the need for rapid subsampling and appropriate storage of sediment core samples for paleogenomic studies.

Investigating variability in microbial community composition in replicate environmental DNA samples down lake sediment cores

Lake sediments are natural archives that accumulate information on biological communities and their surrounding catchments. Paleolimnology has traditionally focussed on identifying fossilized organisms to reconstruct past environments. In the last decade, the application of molecular methodologies has increased in paleolimnological studies, but further research investigating factors such as sample heterogeneity and DNA degradation are required. In the present study we investigated bacterial community heterogeneity (16S rRNA metabarcoding) within depth slices (1-cm width). Sediment cores were collected from three lakes with differing sediment compositions. Samples were collected from a variety of depths which represent a period of time of approximately 1,200 years. Triplicate samples were collected from each depth slice and bacterial 16S rRNA metabarcoding was undertaken on each sample. Accumulation curves demonstrated that except for the deepest (oldest) slices, the combination of three replicate samples were insufficient to characterise the entire bacterial diversity. However, shared Amplicon Sequence Variants (ASVs) accounted for the majority of the reads in each depth slice (max. shared proportional read abundance 96%, 86%, 65% in the three lakes). Replicates within a depth slice generally clustered together in the Non-metric multidimensional scaling analysis. There was high community dissimilarity in older sediment in one of the cores, which was likely due to the laminae in the sediment core not being horizontal. Given that most paleolimnology studies explore broad scale shifts in community structure rather than seeking to identify rare species, this study demonstrates that a single sample is adequate to characterise shifts in dominant bacterial ASVs.

Lake Sedimentary DNA Research on Past Terrestrial and Aquatic Biodiversity: Overview and Recommendations

The use of lake sedimentary DNA to track the long-term changes in both terrestrial and aquatic biota is a rapidly advancing field in paleoecological research. Although largely applied nowadays, knowledge gaps remain in this field and there is therefore still research to be conducted to ensure the reliability of the sedimentary DNA signal. Building on the most recent literature and seven original case studies, we synthesize the state-of-the-art analytical procedures for effective sampling, extraction, amplification, quantification and/or generation of DNA inventories from sedimentary ancient DNA (sedaDNA) via high-throughput sequencing technologies. We provide recommendations based on current knowledge and best practises.

Landscape Setting Drives the Microbial Eukaryotic Community Structure in Four Swedish Mountain Lakes over the Holocene

On the annual and interannual scales, lake microbial communities are known to be heavily influenced by environmental conditions both in the lake and in its terrestrial surroundings. However, the influence of landscape setting and environmental change on shaping these communities over a longer (millennial) timescale is rarely studied. Here, we applied an 18S metabarcoding approach to DNA preserved in Holocene sediment records from two pairs of co-located Swedish mountain lakes. Our data revealed that the microbial eukaryotic communities were strongly influenced by catchment characteristics rather than location. More precisely, the microbial communities from the two bedrock lakes were largely dominated by unclassified Alveolata, while the peatland lakes showed a more diverse microbial community, with Ciliophora, Chlorophyta and Chytrids among the more predominant groups. Furthermore, for the two bedrock-dominated lakes-where the oldest DNA samples are dated to only a few hundred years after the lake formation-certain Alveolata, Chlorophytes, Stramenopiles and Rhizaria taxa were found prevalent throughout all the sediment profiles. Our work highlights the importance of species sorting due to landscape setting and the persistence of microbial eukaryotic diversity over millennial timescales in shaping modern lake microbial communities.