Identification of photosynthetic plankton communities using sedimentary ancient DNA and their response to late-Holocene climate change on the Tibetan Plateau

Ancient DNA reveals potentially toxic cyanobacteria increasing with climate change

Cyanobacterial blooms in freshwater systems are a global threat to human and aquatic ecosystem health, exhibiting particularly harmful effects when toxin-producing taxa are present. While climatic change and nutrient over-enrichment control the global expansion of total cyanobacterial blooms, it remains unknown to what extent this expansion reflected cyanobacterial assemblage due to the scarcity of long-term monitoring data. Here we use high-throughput sequencing of sedimentary DNA to track ∼100 years of changes in cyanobacterial community in hyper-eutrophic Lake Taihu, China's third largest freshwater lake and the key water source for ∼30 million people. A steady increase in the abundance of Microcystis (as potential toxin producers) during the past thirty years was correlated with increasing temperatures and declining wind speeds, but not with temporal trends in lakewater nutrient concentrations, highlighting recent climate effects on potentially increasing toxin-producing taxa. The socio-environmental repercussions of these findings are worrisome as continued anthropogenic climate change may counteract nutrient amelioration efforts in this critical freshwater resource.

Distinct strategies of the habitat generalists and specialists in sediment of Tibetan lakes

Microbial metacommunities normally comprise generalists and specialists. Uncovering the mechanisms underlying the diversity patterns of these two sub-communities is crucial for aquatic biodiversity maintenance. However, little is known about the ecological assembly processes and co-occurrence patterns of the habitat generalists and specialists across large spatial scales in plateau lake sediments, particularly regarding their environmental adaptations. Here, we investigated assembly processes of the habitat generalists and specialists in sediment of Tibetan lakes and their role in the stability of metacommunity co-occurrence network. Our results showed that the habitat generalists exhibited broader environmental thresholds and closer phylogenetic clustering than specialist counterparts. In contrast, the specialists exhibited stronger phylogenetic signals of ecological preferences compared with the habitat generalists. Stochastic processes dominated the habitat generalist (63.2%) and specialist (81.3%) community assembly. Sediment pH was the major factor mediating the balance between stochastic and deterministic processes in the habitat generalists and specialists. In addition, revealed by network analysis, the habitat specialists played a greater role in maintaining the stability of metacommunity co-occurrence network. The insights gained from this study can be helpful to understand the mechanisms underlying maintenance of sediment microbial diversity in plateau lakes. This article is protected by copyright. All rights reserved.

Sediment archives reveal irreversible shifts in plankton communities after World War II and agricultural pollution

To evaluate the stability and resilience¹ of coastal ecosystem communities to perturbations that occurred during the Anthropocene,² pre-industrial biodiversity baselines inferred from paleoarchives are needed.^3,4 The study of ancient DNA (aDNA) from sediments (sedaDNA)⁵ has provided valuable information about past dynamics of microbial species^6-8 and communities^9-18 in relation to ecosystem variations. Shifts in planktonic protist communities might significantly affect marine ecosystems through cascading effects,^19-21 and therefore the analysis of this compartment is essential for the assessment of ecosystem variations. Here, sediment cores collected from different sites of the Bay of Brest (northeast Atlantic, France) allowed ca. 1,400 years of retrospective analyses of the effects of human pollution on marine protists. Comparison of sedaDNA extractions and metabarcoding analyses with different barcode regions (V4 and V7 18S rDNA) revealed that protist assemblages in ancient sediments are mainly composed of species known to produce resting stages. Heavy-metal pollution traces in sediments were ascribed to the World War II period and coincided with community shifts within dinoflagellates and stramenopiles. After the war and especially from the 1980s to 1990s, protist genera shifts followed chronic contaminations of agricultural origin. Community composition reconstruction over time showed that there was no recovery to a Middle Ages baseline composition. This demonstrates the irreversibility of the observed shifts after the cumulative effect of war and agricultural pollutions. Developing a paleoecological approach, this study highlights how human contaminations irreversibly affect marine microbial compartments, which contributes to the debate on coastal ecosystem preservation and restoration.

Subseafloor Archaea reflect 139 kyrs of paleodepositional changes in the northern Red Sea

The vertical distribution of subseafloor archaeal communities is thought to be primarily controlled by in situ conditions in sediments such as the availability of electron acceptors and donors, although sharp community shifts have also been observed at lithological boundaries suggesting that at least a subset of vertically stratified Archaea form a long-term genetic record of coinciding environmental conditions that occurred at the time of sediment deposition. To substantiate this possibility, we performed a highly resolved 16S rRNA gene survey of vertically stratified archaeal communities paired with paleo-oceanographic proxies in a sedimentary record from the northern Red Sea spanning the last glacial-interglacial cycle (i.e., marine isotope stages 1-6; MIS1-6). Our results show a strong significant correlation between subseafloor archaeal communities and drastic paleodepositional changes associated with glacial low vs. interglacial high stands (ANOSIM; R = .73; p = .001) and only a moderately strong correlation with lithological changes. Bathyarchaeota, Lokiarchaeota, MBGA, and DHVEG-1 were the most abundant identified archaeal groups. Whether they represented ancient cell lines from the time of deposition or migrated to the specific sedimentary horizons after deposition remains speculative. However, we show that the majority of sedimentary archaeal tetraether membrane lipids were of allochthonous origin and not produced in situ. Slow post-burial growth under energy-limited conditions would explain why the downcore distribution of these dominant archaeal groups still indirectly reflect changes in the paleodepositional environment that prevailed during the analyzed marine isotope stages. In addition, archaea seeded from the overlying water column such as Thaumarchaeota and group II and III Euryarchaeota, which were likely not have been able to subsist after burial, were identified from a lower abundance of preserved sedimentary DNA signatures, and represented direct markers of paleoenvironmental changes in the Red Sea spanning the last six marine isotope stages.

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.

An optimized method for the extraction of ancient eukaryote DNA from marine sediments

Marine sedimentary ancient DNA (sedaDNA) provides a powerful means to reconstruct marine palaeo-communities across the food web. However, currently there are few optimized sedaDNA extraction protocols available to maximize the yield of small DNA fragments typical of ancient DNA (aDNA) across a broad diversity of eukaryotes. We compared seven combinations of sedaDNA extraction treatments and sequencing library preparations using marine sediments collected at a water depth of 104 m off Maria Island, Tasmania, in 2018. These seven methods contrasted frozen versus refrigerated sediment, bead-beating induced cell lysis versus ethylenediaminetetraacetic acid (EDTA) incubation, DNA binding in silica spin columns versus in silica-solution, diluted versus undiluted DNA in shotgun library preparations to test potential inhibition issues during amplification steps, and size-selection of low molecular-weight (LMW) DNA to increase the extraction efficiency of sedaDNA. Maximum efficiency was obtained from frozen sediments subjected to a combination of EDTA incubation and bead-beating, DNA binding in silica-solution, and undiluted DNA in shotgun libraries, across 45 marine eukaryotic taxa. We present an optimized extraction protocol integrating these steps, with an optional post-library LMW size-selection step to retain DNA fragments of ≤500 base pairs. We also describe a stringent bioinformatic filtering approach for metagenomic data and provide a comprehensive list of contaminants as a reference for future sedaDNA studies. The new extraction and data-processing protocol should improve quantitative paleo-monitoring of eukaryotes from marine sediments, as well as other studies relying on the detection of highly fragmented and degraded eukaryote DNA in sediments.

High Diversity of Myocyanophage in Various Aquatic Environments Revealed by High-Throughput Sequencing of Major Capsid Protein Gene With a New Set of Primers

Myocyanophages, a group of viruses infecting cyanobacteria, are abundant and play important roles in elemental cycling. Here we investigated the particle-associated viral communities retained on 0.2 μm filters and in sediment samples (representing ancient cyanophage communities) from four ocean and three lake locations, using high-throughput sequencing and a newly designed primer pair targeting a gene fragment (∼145-bp in length) encoding the cyanophage gp23 major capsid protein (MCP). Diverse viral communities were detected in all samples. The fragments of 142-, 145-, and 148-bp in length were most abundant in the amplicons, and most sequences (>92%) belonged to cyanophages. Additionally, different sequencing depths resulted in different diversity estimates of the viral community. Operational taxonomic units obtained from deep sequencing of the MCP gene covered the majority of those obtained from shallow sequencing, suggesting that deep sequencing exhibited a more complete picture of cyanophage community than shallow sequencing. Our results also revealed a wide geographic distribution of marine myocyanophages, i.e., higher dissimilarities of the myocyanophage communities corresponded with the larger distances between the sampling sites. Collectively, this study suggests that the newly designed primer pair can be effectively used to study the community and diversity of myocyanophage from different environments, and the high-throughput sequencing represents a good method to understand viral diversity.

Specific Chemical and Genetic Markers Revealed a Thousands-Year Presence of Toxic Nodularia spumigena in the Baltic Sea

In the Baltic Sea, diazotrophic cyanobacteria have been present for thousands of years, over the whole brackish water phase of the ecosystem. However, our knowledge about the species composition of the cyanobacterial community is limited to the last several decades. In the current study, the presence of species-specific chemical and genetic markers in deep sediments were analyzed to increase the existing knowledge on the history of toxic Nodularia spumigena blooms in the Baltic Sea. As chemical markers, three cyclic nonribosomal peptides were applied: the hepatotoxic nodularin, which in the sea was detected solely in N. spumigena, and two anabaenopeptins (AP827 and AP883a) characteristic of two different chemotypes of this species. From the same sediment samples, DNA was isolated and the gene involved in biosynthesis of nodularin, as well as the phycocyanin intergenic spacer region (PC-IGS), were amplified. The results of chemical and genetic analyses proved for the first time the thousands-year presence of toxic N. spumigena in the Baltic Sea. They also indicated that through all this time, the same two sub-populations of the species co-existed.

Evolutionary Patterns and Processes: Lessons from Ancient DNA

Ever since its emergence in 1984, the field of ancient DNA has struggled to overcome the challenges related to the decay of DNA molecules in the fossil record. With the recent development of high-throughput DNA sequencing technologies and molecular techniques tailored to ultra-damaged templates, it has now come of age, merging together approaches in phylogenomics, population genomics, epigenomics, and metagenomics. Leveraging on complete temporal sample series, ancient DNA provides direct access to the most important dimension in evolution—time, allowing a wealth of fundamental evolutionary processes to be addressed at unprecedented resolution. This review taps into the most recent findings in ancient DNA research to present analyses of ancient genomic and metagenomic data.

Localized high abundance of Marine Group II archaea in the subtropical Pearl River Estuary: implications for their niche adaptation

Marine Group II archaea are widely distributed in global oceans and dominate the total archaeal community within the upper euphotic zone of temperate waters. However, factors controlling the distribution of MGII are poorly delineated and the physiology and ecological functions of these still-uncultured organisms remain elusive. In this study, we investigated the planktonic MGII associated with particles and in free-living forms in the Pearl River Estuary (PRE) over a 10-month period. We detected high abundance of particle-associated MGII in PRE (up to ∼10⁸ 16S rRNA gene copies/l), which was around 10-fold higher than the free-living MGII in the same region, and an order of magnitude higher than previously reported in other marine environments. 10‰ salinity appeared to be a threshold value for these MGII because MGII abundance decreased sharply below it. Above 10‰ salinity, the abundance of MGII on the particles was positively correlated with phototrophs and MGII in the surface water was negatively correlated with irradiance. However, the abundances of those free-living MGII showed positive correlations with salinity and temperature, suggesting the different physiological characteristics between particle-attached and free-living MGIIs. A nearly completely assembled metagenome, MGIIa_P, was recovered using metagenome binning methods. Compared with the other two MGII genomes from surface ocean, MGIIa_P contained higher proportions of glycoside hydrolases, indicating the ability of MGIIa_P to hydrolyse glycosidic bonds in complex sugars in PRE. MGIIa_P is the first assembled MGII metagenome containing a catalase gene, which might be involved in scavenging reactive oxygen species generated by the abundant phototrophs in the eutrophic PRE. Our study presented the widespread and high abundance of MGII in the water columns of PRE, and characterized the determinant abiotic factors affecting their distribution. Their association with heterotrophs, preference for particles and resourceful metabolic traits indicate MGII might play a significant role in metabolising organic matters in the PRE and other temperate estuarine systems.

Eukaryotic diversity in late Pleistocene marine sediments around a shallow methane hydrate deposit in the Japan Sea

Marine sediments contain eukaryotic DNA deposited from overlying water columns. However, a large proportion of deposited eukaryotic DNA is aerobically biodegraded in shallow marine sediments. Cold seep sediments are often anaerobic near the sediment-water interface, so eukaryotic DNA in such sediments is expected to be preserved. We investigated deeply buried marine sediments in the Japan Sea, where a methane hydrate deposit is associated with cold seeps. Quantitative PCR analysis revealed the reproducible recovery of eukaryotic DNA in marine sediments at depths up to 31.0 m in the vicinity of the methane hydrate deposit. In contrast, the reproducible recovery of eukaryotic DNA was limited to a shallow depth (8.3 m) in marine sediments not adjacent to the methane hydrate deposit in the same area. Pyrosequencing of an 18S rRNA gene variable region generated 1,276-3,307 reads per sample, which was sufficient to cover the biodiversity based on rarefaction curves. Phylogenetic analysis revealed that most of the eukaryotic DNA originated from radiolarian genera of the class Chaunacanthida, which have SrSO4 skeletons, the sea grass genus Zostera, and the seaweed genus Sargassum. Eukaryotic DNA originating from other planktonic fauna and land plants was also detected. Diatom sequences closely related to Thalassiosira spp., indicative of cold climates, were obtained from sediments deposited during the last glacial period (MIS-2). Plant sequences of the genera Alnus, Micromonas, and Ulmus were found in sediments deposited during the warm interstadial period (MIS-3). These results suggest the long-term persistence of eukaryotic DNA from terrestrial and aquatic sources in marine sediments associated with cold seeps, and that the genetic information from eukaryotic DNA from deeply buried marine sediments associated with cold seeps can be used to reconstruct environments and ecosystems from the past.

Ancient plant DNA in lake sediments

Contents 924 I. 925 II. 925 III. 927 IV. 929 V. 930 VI. 930 VII. 931 VIII. 933 IX. 935 X. 936 XI. 938 938 References 938 SUMMARY: Recent advances in sequencing technologies now permit the analyses of plant DNA from fossil samples (ancient plant DNA, plant aDNA), and thus enable the molecular reconstruction of palaeofloras. Hitherto, ancient frozen soils have proved excellent in preserving DNA molecules, and have thus been the most commonly used source of plant aDNA. However, DNA from soil mainly represents taxa growing a few metres from the sampling point. Lakes have larger catchment areas and recent studies have suggested that plant aDNA from lake sediments is a more powerful tool for palaeofloristic reconstruction. Furthermore, lakes can be found globally in nearly all environments, and are therefore not limited to perennially frozen areas. Here, we review the latest approaches and methods for the study of plant aDNA from lake sediments and discuss the progress made up to the present. We argue that aDNA analyses add new and additional perspectives for the study of ancient plant populations and, in time, will provide higher taxonomic resolution and more precise estimation of abundance. Despite this, key questions and challenges remain for such plant aDNA studies. Finally, we provide guidelines on technical issues, including lake selection, and we suggest directions for future research on plant aDNA studies in lake sediments.

Sedimentary DNA Reveals Cyanobacterial Community Diversity over 200 Years in Two Perialpine Lakes

We reconstructed cyanobacterial community structure and phylogeny using DNA that was isolated from layers of stratified sediments spanning 200 years of lake history in the perialpine lakes Greifensee and Lake Zurich (Switzerland). Community analysis based on amplification and sequencing of a 400-nucleotide (nt)-long 16S rRNA fragment specific to Cyanobacteria revealed operational taxonomic units (OTUs) capturing the whole phylum, including representatives of a newly characterized clade termed Melainabacteria, which shares common ancestry with Cyanobacteria and has not been previously described in lakes. The reconstruction of cyanobacterial richness and phylogenetic structure was validated using a data set consisting of 40 years of pelagic microscopic counts from each lake. We identified the OTUs assigned to common taxa known to be present in Greifensee and Lake Zurich and found a strong and significant relationship (adjusted R² = 0.89; P < 0.001) between pelagic species richness in water and OTU richness in the sediments. The water-sediment richness relationship varied between cyanobacterial orders, indicating that the richness of Chroococcales and Synechococcales may be underestimated by microscopy. PCR detection of the microcystin synthetase gene mcyA confirmed the presence of potentially toxic cyanobacterial taxa over recent years in Greifensee and throughout the last century in Lake Zurich. The approach presented in this study demonstrates that it is possible to reconstruct past pelagic cyanobacterial communities in lakes where the integrity of the sedimentary archive is well preserved and to explore changes in phylogenetic and functional diversity over decade-to-century timescales.

IMPORTANCE

Cyanobacterial blooms can produce toxins that affect water quality, especially under eutrophic conditions, which are a consequence of human-induced climate warming and increased nutrient availability. Lakes worldwide have suffered from regular cyanobacterial blooms over the last century. The lack of long-term data limits our understanding of how these blooms form. We successfully reconstructed the past diversity of whole cyanobacterial communities over two hundred years by sequencing genes preserved in the sediments of two perialpine lakes in Switzerland. We identified changes in diversity over time and validated our results using existing data collected in the same two lakes over the past 40 years. This work shows the potential of our approach for addressing important ecological questions about the effects of a changing environment on lake ecology.

Temporal Succession of Ancient Phytoplankton Community in Qinghai Lake and Implication for Paleo-environmental Change

Tibetan lake sediments in NW China are sensitive recorders of climate change. However, many important plankton members do not leave any microscopic features in sedimentary records. Here we used ancient DNA preserved in Qinghai Lake sediments to reconstruct the temporal succession of plankton communities in the past 18,500 years. Our results showed that seven classes and sixteen genera of phytoplankton in the lake underwent major temporal changes, in correlation with known climatic events. Trebouxiophyceae and Eustigmatophyceae were predominant during the cold periods, whereas Chlorophyceae, Phaeophyceae, Xanthophyceae, Bacillariophyceae, and Cyanophyceae were abundant during the warm periods. The inferred changes in temperature, nutrients, precipitation, and salinity, as driven by the Westerlies and summer Monsoon strength, likely contributed to these observed temporal changes. Based on these correlations, we propose the phytoplankton index as a proxy to reconstruct the stadial versus interstadial climate change history in Qinghai Lake. This taxon-specific index is free of terrestrial contamination, sensitive to short-term climatic oscillations, and continuous in recording all climatic events in the lake. The validity of this index and its applicability to other lakes is demonstrated by its good correlations with multiple climate records of Qinghai Lake and another lake on the Tibetan Plateau, Kusai Lake.

Sedimentary archaeal amoA gene abundance reflects historic nutrient level and salinity fluctuations in Qinghai Lake, Tibetan Plateau

Integration of DNA derived from ancient phototrophs with their characteristic lipid biomarkers has been successfully employed to reconstruct paleoenvironmental conditions. However, it is poorly known that whether the DNA and lipids of microbial functional aerobes (such as ammonia-oxidizing archaea: AOA) can be used for reconstructing past environmental conditions. Here we identify and quantify the AOA amoA genes (encoding the alpha subunit of ammonia monooxygenases) preserved in a 5.8-m sediment core (spanning the last 18,500 years) from Qinghai Lake. Parallel analyses revealed that low amoA gene abundance corresponded to high total organic carbon (TOC) and salinity, while high amoA gene abundance corresponded to low TOC and salinity. In the Qinghai Lake region, TOC can serve as an indicator of paleo-productivity and paleo-precipitation, which is related to historic nutrient input and salinity. So our data suggest that temporal variation of AOA amoA gene abundance preserved in Qinghai Lake sediment may reflect the variations of nutrient level and salinity throughout the late Pleistocene and Holocene in the Qinghai Lake region.