A comparison of sedimentary DNA and pollen from lake sediments in recording vegetation composition at the Siberian treeline

A 1 Ma sedimentary ancient DNA (sedaDNA) record of catchment vegetation changes and the developmental history of tropical Lake Towuti (Sulawesi, Indonesia)

Studying past ecosystems from ancient environmental DNA preserved in lake sediments (sedaDNA) is a rapidly expanding field. This research has mainly involved Holocene sediments from lakes in cool climates, with little known about the suitability of sedaDNA to reconstruct substantially older ecosystems in the warm tropics. Here, we report the successful recovery of chloroplast trnL (UAA) sequences (trnL-P6 loop) from the sedimentary record of Lake Towuti (Sulawesi, Indonesia) to elucidate changes in regional tropical vegetation assemblages during the lake's Late Quaternary paleodepositional history. After the stringent removal of contaminants and sequence artifacts, taxonomic assignment of the remaining genuine trnL-P6 reads showed that native nitrogen-fixing legumes, C3 grasses, and shallow wetland vegetation (Alocasia) were most strongly associated with >1-million-year-old (>1 Ma) peats and silts (114-98.8 m composite depth; mcd), which were deposited in a landscape of active river channels, shallow lakes, and peat-swamps. A statistically significant shift toward partly submerged shoreline vegetation that was likely rooted in anoxic muddy soils (i.e., peatland forest trees and wetland C3 grasses (Oryzaceae) and nutrient-demanding aquatic herbs (presumably Oenanthe javanica)) occurred at 76 mcd (~0.8 Ma), ~0.2 Ma after the transition into a permanent lake. This wetland vegetation was most strongly associated with diatom ooze (46-37 mcd), thought to be deposited during maximum nutrient availability and primary productivity. Herbs (Brassicaceae), trees/shrubs (Fabaceae and Theaceae), and C3 grasses correlated with inorganic parameters, indicating increased drainage of ultramafic sediments and laterite soils from the lakes' catchment, particularly at times of inferred drying. Downcore variability in trnL-P6 from tropical forest trees (Toona), shady ground cover herbs (Zingiberaceae), and tree orchids (Luisia) most strongly correlated with sediments of a predominantly felsic signature considered to be originating from the catchment of the Loeha River draining into Lake Towuti during wetter climate conditions. However, the co-correlation with dry climate-adapted trees (i.e., Castanopsis or Lithocarpus) plus C4 grasses suggests that increased precipitation seasonality also contributed to the increased drainage of felsic Loeha River sediments. This multiproxy approach shows that despite elevated in situ temperatures, tropical lake sediments potentially comprise long-term archives of ancient environmental DNA for reconstructing ecosystems, which warrants further exploration.

Earthquake, floods and changing land use history: A 200-year overview of environmental changes in Selenga River basin as indicated by n-alkanes and related proxies in sediments from shallow lakes

The Selenga River basin, located in southern Siberia, is an important component of the Lake Baikal ecosystem, and comprises approximately 80 % of the Baikal watershed. Within the Selenga River basin, two localized study regions were chosen. The first, the Selenga Delta, is one of the largest inland freshwater floodplains in the world and plays an important role in the ecosystem functioning of Baikal. It purifies the river waters before they enter the lake and acts as a refuge for many of Baikal's endemic species. The second location, the Gusinoozersk region, is southwest of Lake Baikal and the Selenga Delta, and was chosen as a more heavily industrialized region within the Selenga River basin. Anthropogenic activities, including industry, urban settlements, aquaculture and agriculture, have historically increased ecological damage within this area. We assessed possible drivers of changes in sedimentary organic matter (OM) composition within two shallow lakes (SLNG04 and Black Lake), located in the Selenga Delta and the Selenga watershed, respectively. We focused on individual n-alkanes, one of the most abundant and common lipids used to provide information on past vegetation and used multivariate statistics to disentangle changes in the sources of sedimentary OM over time. The depositional OM history of SLNG04B core can be divided in four zones: (i) major influence of non-emergent vascular plants, typically found in transitional environments (ca. 1835 to ca. 1875); (ii) increased influence of grasses/herbs (ca. 1880 to ca. 1910); (iii) transition from non-emergent vascular plants and grasses/herbs to submerged and floating macrophytes and phytoplankton (ca. 1915 to ca. 1945); (iv) maintenance of autochthonous OM from submerged and floating macrophytes and phytoplankton (ca. 1945 to ca. 2014). The depositional OM history of the Black Lake core can be divided in two main zones: (i) major influence of non-emergent vascular plants and submerged and floating macrophytes (ca. 1915 to ca. 1980); (ii) increased influence of grasses/herbs and phytoplankton (ca. 1980 to ca. 2010). Natural events (e.g., an earthquake in 1862 caused flooding and subsidence of much of the land surrounding SLNG04 lake and a further catastrophic flood event in 1897) and anthropogenic activities (e.g., nutrient pollution from expansion of agricultural and livestock population) changed the composition of sedimentary OM resulting in ecological shifts across trophic levels in the Selenga River basin.

MetaDamage tool: Examining post-mortem damage in sedaDNA on a metagenomic scale

The use of metagenomic datasets to support ancient sedimentary DNA (sedaDNA) for paleoecological reconstruction has been demonstrated to be a powerful tool to understand multi-organism responses to climatic shifts and events. Authentication remains integral to the ancient DNA discipline, and this extends to sedaDNA analysis. Furthermore, distinguishing authentic sedaDNA from contamination or modern material also allows for a better understanding of broader questions in sedaDNA research, such as formation processes, source and catchment, and post-depositional processes. Existing tools for the detection of damage signals are designed for single-taxon input, require a priori organism specification, and require a significant number of input sequences to establish a signal. It is therefore often difficult to identify an established cytosine deamination rate consistent with ancient DNA across a sediment sample. In this study, we present MetaDamage, a tool that examines cytosine deamination on a metagenomic (all organisms) scale for multiple previously undetermined taxa and can produce a damage profile based on a few hundred reads. We outline the development and testing of the MetaDamage tool using both authentic sedaDNA sequences and simulated data to demonstrate the resolution in which MetaDamage can identify deamination levels consistent with the presence of ancient DNA. The MetaDamage tool offers a method for the initial assessment of the presence of sedaDNA and a better understanding of key questions of preservation for paleoecological reconstruction.

High resolution ancient sedimentary DNA shows that alpine plant diversity is associated with human land use and climate change

The European Alps are highly rich in species, but their future may be threatened by ongoing changes in human land use and climate. Here, we reconstructed vegetation, temperature, human impact and livestock over the past ~12,000 years from Lake Sulsseewli, based on sedimentary ancient plant and mammal DNA, pollen, spores, chironomids, and microcharcoal. We assembled a highly-complete local DNA reference library (PhyloAlps, 3923 plant taxa), and used this to obtain an exceptionally rich sedaDNA record of 366 plant taxa. Vegetation mainly responded to climate during the early Holocene, while human activity had an additional influence on vegetation from 6 ka onwards. Land-use shifted from episodic grazing during the Neolithic and Bronze Age to agropastoralism in the Middle Ages. Associated human deforestation allowed the coexistence of plant species typically found at different elevational belts, leading to levels of plant richness that characterise the current high diversity of this region. Our findings indicate a positive association between low intensity agropastoral activities and precipitation with the maintenance of the unique subalpine and alpine plant diversity of the European Alps.

Assessing the lysis of diverse pollen from bulk environmental samples for DNA metabarcoding

Pollen is ubiquitous year-round in bulk environmental samples and can provide useful information on previous and current plant communities. Characterization of pollen has traditionally been completed based on morphology, requiring significant time and expertise. DNA metabarcoding is a promising approach for characterizing pollen from bulk environmental samples, but accuracy hinges on successful lysis of pollen grains to free template DNA. In this study, we assessed the lysis of morphologically and taxonomically diverse pollen from one of the most common bulk environmental sample types for DNA metabarcoding, surface soil. To achieve this, a four species artificial pollen mixture was spiked into surface soils collected from Colorado, North Carolina, and Pennsylvania, and subsequently subjected to DNA extraction using both the PowerSoil and PowerSoil Pro Kits (Qiagen) with a heated incubation (either 65 °C or 90 °C). Amplification and Illumina sequencing of the internal transcribed spacer subunit 2 (ITS2) was completed in duplicate for each sample (total n, 76), and the resulting sequencing reads taxonomically identified using GenBank. The PowerSoil Pro Kit statistically outperformed the PowerSoil Kit for total DNA yield. When using either kit, incubation temperature (65 °C or 90 °C) used had no impact on the recovery of DNA, plant amplicon sequence variants (ASVs), or total plant ITS2 reads. This study highlighted that lysis of pollen in bulk environmental samples is feasible using commercially available kits, and downstream DNA metabarcoding can be used to accurately characterize pollen DNA from such sample types.

Modern Pollen Assemblages From Lake Sediments and Soil in East Siberia and Relative Pollen Productivity Estimates for Major Taxa

Modern pollen–vegetation–climate relationships underpin palaeovegetation and palaeoclimate reconstructions from fossil pollen records. East Siberia is an ideal area for investigating the relationships between modern pollen assemblages and near natural vegetation under cold continental climate conditions. Reliable pollen-based quantitative vegetation and climate reconstructions are still scarce due to the limited number of modern pollen datasets. Furthermore, differences in pollen representation of samples from lake sediments and soils are not well understood. Here, we present a new pollen dataset of 48 moss/soil and 24 lake surface-sediment samples collected in Chukotka and central Yakutia in East Siberia. The pollen–vegetation–climate relationships were investigated by ordination analyses. Generally, tundra and taiga vegetation types can be well distinguished in the surface pollen assemblages. Moss/soil and lake samples contain generally similar pollen assemblages as revealed by a Procrustes comparison with some exceptions. Overall, modern pollen assemblages reflect the temperature and precipitation gradients in the study areas as revealed by constrained ordination analysis. We estimate the relative pollen productivity (RPP) of major taxa and the relevant source area of pollen (RSAP) for moss/soil samples from Chukotka and central Yakutia using Extended R-Value (ERV) analysis. The RSAP of the tundra-forest transition area in Chukotka and taiga area in central Yakutia are ca. 1300 and 360 m, respectively. For Chukotka, RPPs relative to both Poaceae and Ericaceae were estimated while RPPs for central Yakutia were relative only to Ericaceae. Relative to Ericaceae (reference taxon, RPP = 1), Larix, Betula, Picea, and Pinus are overrepresented while Alnus, Cyperaceae, Poaceae, and Salix are underrepresented in the pollen spectra. Our estimates are in general agreement with previously published values and provide the basis for reliable quantitative reconstructions of East Siberian vegetation.

Two Millennia of Complexity and Variability in a Perialpine Socioecological System (Savoie, France): The Contribution of Palynology and sedaDNA Analysis

Over the last two millennia, European Alpine ecosystems have experienced major changes in response to the important, yet fluctuating, impact of human activities. This study aims to reconstruct the environmental history of the last 1800 years on the western edge of the Alps by analyzing sediments from Lake Aiguebelette, a large lake located in the perialpine area. We have combined analyses of pollen and other palynomorphs, such as coprophilous fungal spores, together with sedimentary DNA (from plants and mammals) in order to reconstruct both vegetation and land-use histories. A sedimentological and geochemical analysis was also conducted in order to gain an understanding of changes in erosion dynamics in response to landscape modifications that were influenced by climate and human activities. This work highlights alternating phases of anthropization and agricultural abandonment allowing forest recovery. While pollen reflects the major phases of regional deforestation and afforestation related to the dynamic of farming activities, plant DNA provides precise information on the plants cultivated in fields, orchards and vegetable gardens over the past centuries. The combination of mammal DNA and coprophilous fungal spores completes this work by documenting the history of pastoral practices.

Collapse of the mammoth-steppe in central Yukon as revealed by ancient environmental DNA

The temporal and spatial coarseness of megafaunal fossil records complicates attempts to to disentangle the relative impacts of climate change, ecosystem restructuring, and human activities associated with the Late Quaternary extinctions. Advances in the extraction and identification of ancient DNA that was shed into the environment and preserved for millennia in sediment now provides a way to augment discontinuous palaeontological assemblages. Here, we present a 30,000-year sedimentary ancient DNA (sedaDNA) record derived from loessal permafrost silts in the Klondike region of Yukon, Canada. We observe a substantial turnover in ecosystem composition between 13,500 and 10,000 calendar years ago with the rise of woody shrubs and the disappearance of the mammoth-steppe (steppe-tundra) ecosystem. We also identify a lingering signal of Equus sp. (North American horse) and Mammuthus primigenius (woolly mammoth) at multiple sites persisting thousands of years after their supposed extinction from the fossil record.

Plant Sedimentary Ancient DNA From Far East Russia Covering the Last 28,000 Years Reveals Different Assembly Rules in Cold and Warm Climates

Woody plants are expanding into the Arctic in response to the warming climate. The impact on arctic plant communities is not well understood due to the limited knowledge about plant assembly rules. Records of past plant diversity over long time series are rare. Here, we applied sedimentary ancient DNA metabarcoding targeting the P6 loop of the chloroplast trnL gene to a sediment record from Lake Ilirney (central Chukotka, Far Eastern Russia) covering the last 28 thousand years. Our results show that forb-rich steppe-tundra and dwarf-shrub tundra dominated during the cold climate before 14 ka, while deciduous erect-shrub tundra was abundant during the warm period since 14 ka. Larix invasion during the late Holocene substantially lagged behind the likely warmest period between 10 and 6 ka, where the vegetation biomass could be highest. We reveal highest richness during 28–23 ka and a second richness peak during 13–9 ka, with both periods being accompanied by low relative abundance of shrubs. During the cold period before 14 ka, rich plant assemblages were phylogenetically clustered, suggesting low genetic divergence in the assemblages despite the great number of species. This probably originates from environmental filtering along with niche differentiation due to limited resources under harsh environmental conditions. In contrast, during the warmer period after 14 ka, rich plant assemblages were phylogenetically overdispersed. This results from a high number of species which were found to harbor high genetic divergence, likely originating from an erratic recruitment process in the course of warming. Some of our evidence may be of relevance for inferring future arctic plant assembly rules and diversity changes. By analogy to the past, we expect a lagged response of tree invasion. Plant richness might overshoot in the short term; in the long-term, however, the ongoing expansion of deciduous shrubs will eventually result in a phylogenetically more diverse community.

Harnessing the Power of Metabarcoding in the Ecological Interpretation of Plant-Pollinator DNA Data: Strategies and Consequences of Filtering Approaches

Although DNA metabarcoding of pollen mixtures has been increasingly used in the field of pollination biology, methodological and interpretation issues arise due to its high sensitivity. Filtering or maintaining false positives, contaminants, and rare taxa or molecular features could lead to different ecological results. Here, we reviewed how this choice has been addressed in 43 studies featuring pollen DNA metabarcoding, which highlighted a very high heterogeneity of filtering methods. We assessed how these strategies shaped pollen assemblage composition, species richness, and interaction networks. To do so, we compared four processing methods: unfiltering, filtering with a proportional 1% of sample reads, a fixed threshold of 100 reads, and the ROC approach (Receiver Operator Characteristic). The results indicated that filtering impacted species composition and reduced species richness, with ROC emerging as a conservative approach. Moreover, in contrast to unfiltered networks, filtering decreased network Connectance and Entropy, and it increased Modularity and Connectivity, indicating that using cut-off thresholds better describes interactions. Overall, unfiltering might compromise reliable ecological interpretations, unless a study targets rare species. We discuss the suitability of each filtering type, plead for justifying filtering strategies on biological or methodological bases and for developing shared approaches to make future studies more comparable.

Sedimentary ancient DNA reveals a threat of warming-induced alpine habitat loss to Tibetan Plateau plant diversity

Studies along elevational gradients worldwide usually find the highest plant taxa richness in mid-elevation forest belts. Hence, an increase in upper elevation diversity is expected in the course of warming-related treeline rise. Here, we use a time-series approach to infer past taxa richness from sedimentary ancient DNA from the south-eastern Tibetan Plateau over the last ~18,000 years. We find the highest total plant taxa richness during the cool phase after glacier retreat when the area contained extensive and diverse alpine habitats (14-10 ka); followed by a decline when forests expanded during the warm early- to mid-Holocene (10-3.6 ka). Livestock grazing since 3.6 ka promoted plant taxa richness only weakly. Based on these inferred dependencies, our simulation yields a substantive decrease in plant taxa richness in response to warming-related alpine habitat loss over the next centuries. Accordingly, efforts of Tibetan biodiversity conservation should include conclusions from palaeoecological evidence.

Vegetation Reconstruction From Siberia and the Tibetan Plateau Using Modern Analogue Technique–Comparing Sedimentary (Ancient) DNA and Pollen Data

To reconstruct past vegetation from pollen or, more recently, lake sedimentary DNA (sedDNA) data is a common goal in palaeoecology. To overcome the bias of a researcher’s subjective assessment and to assign past assemblages to modern vegetation types quantitatively, the modern analogue technique (MAT) is often used for vegetation reconstruction. However, a rigorous comparison of MAT-derived pollen-based and sedDNA-based vegetation reconstruction is lacking. Here, we assess the dissimilarity between modern taxa assemblages from lake surface-sediments and fossil taxa assemblages from four lake sediment cores from the south-eastern Tibetan Plateau and northern Siberia using receiver operating characteristic (ROC) curves, ordination methods, and Procrustes analyses. Modern sedDNA samples from 190 lakes and pollen samples from 136 lakes were collected from a variety of vegetation types. Our results show that more modern analogues are found with sedDNA than pollen when applying similarly derived thresholds. In particular, there are few modern pollen analogues for open vegetation such as alpine or arctic tundra, limiting the ability of treeline shifts to be clearly reconstructed. In contrast, the shifts in the main vegetation communities are well captured by sedimentary ancient DNA (sedaDNA). For example, pronounced shifts from late-glacial alpine meadow/steppe to early–mid-Holocene coniferous forests to late Holocene Tibetan shrubland vegetation types are reconstructed for Lake Naleng on the south-eastern Tibetan Plateau. Procrustes and PROTEST analyses reveal that intertaxa relationships inferred from modern sedaDNA datasets align with past relationships generally, while intertaxa relationships derived from modern pollen spectra are mostly significantly different from fossil pollen relationships. Overall, we conclude that a quantitative sedaDNA-based vegetation reconstruction using MAT is more reliable than a pollen-based reconstruction, probably because of the more straightforward taphonomy that can relate sedDNA assemblages to the vegetation surrounding the lake.

Ancient plant DNA reveals High Arctic greening during the Last Interglacial

Summer warming is driving a greening trend across the Arctic, with the potential for large-scale amplification of climate change due to vegetation-related feedbacks [Pearson et al., Nat. Clim. Chang. (3), 673-677 (2013)]. Because observational records are sparse and temporally limited, past episodes of Arctic warming can help elucidate the magnitude of vegetation response to temperature change. The Last Interglacial ([LIG], 129,000 to 116,000 y ago) was the most recent episode of Arctic warming on par with predicted 21st century temperature change [Otto-Bliesner et al., Philos. Trans. A Math. Phys. Eng. Sci. (371), 20130097 (2013) and Post et al., SciAdv (5), eaaw9883 (2019)]. However, high-latitude terrestrial records from this period are rare, so LIG vegetation distributions are incompletely known. Pollen-based vegetation reconstructions can be biased by long-distance pollen transport, further obscuring the paleoenvironmental record. Here, we present a LIG vegetation record based on ancient DNA in lake sediment and compare it with fossil pollen. Comprehensive plant community reconstructions through the last and current interglacial (the Holocene) on Baffin Island, Arctic Canada, reveal coherent climate-driven community shifts across both interglacials. Peak LIG warmth featured a ∼400-km northward range shift of dwarf birch, a key woody shrub that is again expanding northward. Greening of the High Arctic-documented here by multiple proxies-likely represented a strong positive feedback on high-latitude LIG warming. Authenticated ancient DNA from this lake sediment also extends the useful preservation window for the technique and highlights the utility of combining traditional and molecular approaches for gleaning paleoenvironmental insights to better anticipate a warmer future.

Holocene chloroplast genetic variation of shrubs (Alnus alnobetula, Betula nana, Salix sp.) at the siberian tundra-taiga ecotone inferred from modern chloroplast genome assembly and sedimentary ancient DNA analyses

Climate warming alters plant composition and population dynamics of arctic ecosystems. In particular, an increase in relative abundance and cover of deciduous shrub species (shrubification) has been recorded. We inferred genetic variation of common shrub species (Alnus alnobetula, Betula nana, Salix sp.) through time. Chloroplast genomes were assembled from modern plants (n = 15) from the Siberian forest-tundra ecotone. Sedimentary ancient DNA (sedaDNA; n = 4) was retrieved from a lake on the southern Taymyr Peninsula and analyzed by metagenomics shotgun sequencing and a hybridization capture approach. For A. alnobetula, analyses of modern DNA showed low intraspecies genetic variability and a clear geographical structure in haplotype distribution. In contrast, B. nana showed high intraspecies genetic diversity and weak geographical structure. Analyses of sedaDNA revealed a decreasing relative abundance of Alnus since 5,400 cal yr BP, whereas Betula and Salix increased. A comparison between genetic variations identified in modern DNA and sedaDNA showed that Alnus variants were maintained over the last 6,700 years in the Taymyr region. In accordance with modern individuals, the variants retrieved from Betula and Salix sedaDNA showed higher genetic diversity. The success of the hybridization capture in retrieving diverged sequences demonstrates the high potential for future studies of plant biodiversity as well as specific genetic variation on ancient DNA from lake sediments. Overall, our results suggest that shrubification has species-specific trajectories. The low genetic diversity in A. alnobetula suggests a local population recruitment and growth response of the already present communities, whereas the higher genetic variability and lack of geographical structure in B. nana may indicate a recruitment from different populations due to more efficient seed dispersal, increasing the genetic connectivity over long distances.

Hybridization capture of larch (Larix Mill.) chloroplast genomes from sedimentary ancient DNA reveals past changes of Siberian forest

Siberian larch (Larix Mill.) forests dominate vast areas of northern Russia and contribute important ecosystem services to the world. It is important to understand the past dynamics of larches in order to predict their likely response to a changing climate in the future. Sedimentary ancient DNA extracted from lake sediment cores can serve as archives to study past vegetation. However, the traditional method of studying sedimentary ancient DNA-metabarcoding-focuses on small fragments, which cannot resolve Larix to species level nor allow a detailed study of population dynamics. Here, we use shotgun sequencing and hybridization capture with long-range PCR-generated baits covering the complete Larix chloroplast genome to study Larix populations from a sediment core reaching back to 6700 years from the Taymyr region in northern Siberia. In comparison with shotgun sequencing, hybridization capture results in an increase in taxonomically classified reads by several orders of magnitude and the recovery of complete chloroplast genomes of Larix. Variation in the chloroplast reads corroborates an invasion of Larix gmelinii into the range of Larix sibirica before 6700 years ago. Since then, both species have been present at the site, although larch populations have decreased with only a few trees remaining in what was once a forested area. This study demonstrates for the first time that hybridization capture applied directly to ancient DNA of plants extracted from lake sediments can provide genome-scale information and is a viable tool for studying past genomic changes in populations of single species, irrespective of a preservation as macrofossil.

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.

Plant diversity in sedimentary DNA obtained from high-latitude (Siberia) and high-elevation lakes (China)

Background

Plant diversity in the Arctic and at high altitudes strongly depends on and rebounds to climatic and environmental variability and is nowadays tremendously impacted by recent climate warming. Therefore, past changes in plant diversity in the high Arctic and high-altitude regions are used to infer climatic and environmental changes through time and allow future predictions. Sedimentary DNA (sedDNA) is an established proxy for the detection of local plant diversity in lake sediments, but still relationships between environmental conditions and preservation of the plant sedDNA proxy are far from being fully understood. Studying modern relationships between environmental conditions and plant sedDNA will improve our understanding under which conditions sedDNA is well-preserved helping to a.) evaluate suitable localities for sedDNA approaches, b.) provide analogues for preservation conditions and c.) conduct reconstruction of plant diversity and climate change. This study investigates modern plant diversity applying a plant-specific metabarcoding approach on sedimentary DNA of surface sediment samples from 262 lake localities covering a large geographical, climatic and ecological gradient. Latitude ranges between 25°N and 73°N and longitude between 81°E and 161°E, including lowland lakes and elevated lakes up to 5168 m a.s.l. Further, our sampling localities cover a climatic gradient ranging in mean annual temperature between -15°C and +18°C and in mean annual precipitation between 36­ and 935 mm. The localities in Siberia span over a large vegetational gradient including tundra, open woodland and boreal forest. Lake localities in China include alpine meadow, shrub, forest and steppe and also cultivated areas. The assessment of plant diversity in the underlying dataset was conducted by a specific plant metabarcoding approach.

New information

We provide a large dataset of genetic plant diversity retrieved from surface sedimentary DNA from lakes in Siberia and China spanning over a large environmental gradient. Our dataset encompasses sedDNA sequence data of 259 surface lake sediments and three soil samples originating from Siberian and Chinese lakes. We used the established chloroplastidal P6 loop trnL marker for plant diversity assessment. The merged, filtered and assigned dataset includes 15,692,944 read counts resulting in 623 unique plant DNA sequence types which have a 100% match to either the EMBL or to the specific Arctic plant reference database. The underlying dataset includes a taxonomic list of identified plants and results from PCR replicates, as well as extraction blanks (BLANKs) and PCR negative controls (NTCs), which were run along with the investigated lake samples. This collection of plant metabarcoding data from modern lake sediments is still ongoing and additional data will be released in the future.

Arctic shrub colonization lagged peak postglacial warmth: Molecular evidence in lake sediment from Arctic Canada

Arctic shrubification is an observable consequence of climate change, already resulting in ecological shifts and global-scale climate feedbacks including changes in land surface albedo and enhanced evapotranspiration. However, the rate at which shrubs can colonize previously glaciated terrain in a warming world is largely unknown. Reconstructions of past vegetation dynamics in conjunction with climate records can provide critical insights into shrubification rates and controls on plant migration, but paleoenvironmental reconstructions based on pollen may be biased by the influx of exotic pollen to tundra settings. Here, we reconstruct past plant communities using sedimentary ancient DNA (sedaDNA), which has a more local source area than pollen. We additionally reconstruct past temperature variability using bacterial cell membrane lipids (branched glycerol dialkyl glycerol tetraethers) and an aquatic productivity indicator (biogenic silica) to evaluate the relative timing of postglacial ecological and climate changes at a lake on southern Baffin Island, Arctic Canada. The sedaDNA record tightly constrains the colonization of dwarf birch (Betula, a thermophilous shrub) to 5.9 ± 0.1 ka, ~3 ka after local deglaciation as determined by cosmogenic ¹⁰ Be moraine dating and >2 ka later than Betula pollen is recorded in nearby lake sediment. We then assess the paleovegetation history within the context of summer temperature and find that paleotemperatures were highest prior to 6.3 ka, followed by cooling in the centuries preceding Betula establishment. Together, these molecular proxies reveal that Betula colonization lagged peak summer temperatures, suggesting that inefficient dispersal, rather than climate, may have limited Arctic shrub migration in this region. In addition, these data suggest that pollen-based climate reconstructions from high latitudes, which rely heavily on the presence and abundance of pollen from thermophilous taxa like Betula, can be compromised by both exotic pollen fluxes and vegetation migration lags.

New insights on lake sediment DNA from the catchment: importance of taphonomic and analytical issues on the record quality

Over the last decade, an increasing number of studies have used lake sediment DNA to trace past landscape changes, agricultural activities or human presence. However, the processes responsible for lake sediment formation and sediment properties might affect DNA records via taphonomic and analytical processes. It is crucial to understand these processes to ensure reliable interpretations for “palaeo” studies. Here, we combined plant and mammal DNA metabarcoding analyses with sedimentological and geochemical analyses from three lake-catchment systems that are characterised by different erosion dynamics. The new insights derived from this approach elucidate and assess issues relating to DNA sources and transfer processes. The sources of eroded materials strongly affect the “catchment-DNA” concentration in the sediments. For instance, erosion of upper organic and organo-mineral soil horizons provides a higher amount of plant DNA in lake sediments than deep horizons, bare soils or glacial flours. Moreover, high erosion rates, along with a well-developed hydrographic network, are proposed as factors positively affecting the representation of the catchment flora. The development of open and agricultural landscapes, which favour the erosion, could thus bias the reconstructed landscape trajectory but help the record of these human activities. Regarding domestic animals, pastoral practices and animal behaviour might affect their DNA record because they control the type of source of DNA (“point” vs. “diffuse”).

Chloroplast and mitochondrial genetic variation of larches at the Siberian tundra-taiga ecotone revealed by de novo assembly

Larix populations at the tundra-taiga ecotone in northern Siberia are highly under-represented in population genetic studies, possibly due to the remoteness of these regions that can only be accessed at extraordinary expense. The genetic signatures of populations in these boundary regions are therefore largely unknown. We aim to generate organelle reference genomes for the detection of single nucleotide polymorphisms (SNPs) that can be used for paleogenetic studies. We present 19 complete chloroplast genomes and mitochondrial genomic sequences of larches from the southern lowlands of the Taymyr Peninsula (northernmost range of Larix gmelinii (Rupr.) Kuzen.), the lower Omoloy River, and the lower Kolyma River (both in the range of Larix cajanderi Mayr). The genomic data reveal 84 chloroplast SNPs and 213 putatively mitochondrial SNPs. Parsimony-based chloroplast haplotype networks show no spatial structure of individuals from different geographic origins, while the mitochondrial haplotype network shows at least a slight spatial structure with haplotypes from the Omoloy and Kolyma populations being more closely related to each other than to most of the haplotypes from the Taymyr populations. Whole genome alignments with publicly available complete chloroplast genomes of different Larix species show that among official plant barcodes only the rcbL gene contains sufficient polymorphisms, but has to be sequenced completely to distinguish the different provenances. We provide 8 novel mitochondrial SNPs that are putatively diagnostic for the separation of L. gmelinii and L. cajanderi, while 4 chloroplast SNPs have the potential to distinguish the L. gmelinii/L. cajanderi group from other Larix species. Our organelle references can be used for a targeted primer and probe design allowing the generation of short amplicons. This is particularly important with regard to future investigations of, for example, the biogeographic history of Larix by screening ancient sedimentary DNA of Larix.

Temporal and spatial patterns of mitochondrial haplotype and species distributions in Siberian larches inferred from ancient environmental DNA and modeling

Changes in species' distributions are classically projected based on their climate envelopes. For Siberian forests, which have a tremendous significance for vegetation-climate feedbacks, this implies future shifts of each of the forest-forming larch (Larix) species to the north-east. However, in addition to abiotic factors, reliable projections must assess the role of historical biogeography and biotic interactions. Here, we use sedimentary ancient DNA and individual-based modelling to investigate the distribution of larch species and mitochondrial haplotypes through space and time across the treeline ecotone on the southern Taymyr peninsula, which at the same time presents a boundary area of two larch species. We find spatial and temporal patterns, which suggest that forest density is the most influential driver determining the precise distribution of species and mitochondrial haplotypes. This suggests a strong influence of competition on the species' range shifts. These findings imply possible climate change outcomes that are directly opposed to projections based purely on climate envelopes. Investigations of such fine-scale processes of biodiversity change through time are possible using paleoenvironmental DNA, which is available much more readily than visible fossils and can provide information at a level of resolution that is not reached in classical palaeoecology.

Minimizing polymerase biases in metabarcoding

DNA metabarcoding is an increasingly popular method to characterize and quantify biodiversity in environmental samples. Metabarcoding approaches simultaneously amplify a short, variable genomic region, or "barcode," from a broad taxonomic group via the polymerase chain reaction (PCR), using universal primers that anneal to flanking conserved regions. Results of these experiments are reported as occurrence data, which provide a list of taxa amplified from the sample, or relative abundance data, which measure the relative contribution of each taxon to the overall composition of amplified product. The accuracy of both occurrence and relative abundance estimates can be affected by a variety of biological and technical biases. For example, taxa with larger biomass may be better represented in environmental samples than those with smaller biomass. Here, we explore how polymerase choice, a potential source of technical bias, might influence results in metabarcoding experiments. We compared potential biases of six commercially available polymerases using a combination of mixtures of amplifiable synthetic sequences and real sedimentary DNA extracts. We find that polymerase choice can affect both occurrence and relative abundance estimates and that the main source of this bias appears to be polymerase preference for sequences with specific GC contents. We further recommend an experimental approach for metabarcoding based on results of our synthetic experiments.

Plant DNA metabarcoding of lake sediments: How does it represent the contemporary vegetation

Metabarcoding of lake sediments have been shown to reveal current and past biodiversity, but little is known about the degree to which taxa growing in the vegetation are represented in environmental DNA (eDNA) records. We analysed composition of lake and catchment vegetation and vascular plant eDNA at 11 lakes in northern Norway. Out of 489 records of taxa growing within 2 m from the lake shore, 17–49% (mean 31%) of the identifiable taxa recorded were detected with eDNA. Of the 217 eDNA records of 47 plant taxa in the 11 lakes, 73% and 12% matched taxa recorded in vegetation surveys within 2 m and up to about 50 m away from the lakeshore, respectively, whereas 16% were not recorded in the vegetation surveys of the same lake. The latter include taxa likely overlooked in the vegetation surveys or growing outside the survey area. The percentages detected were 61, 47, 25, and 15 for dominant, common, scattered, and rare taxa, respectively. Similar numbers for aquatic plants were 88, 88, 33 and 62%, respectively. Detection rate and taxonomic resolution varied among plant families and functional groups with good detection of e.g. Ericaceae, Roseaceae, deciduous trees, ferns, club mosses and aquatics. The representation of terrestrial taxa in eDNA depends on both their distance from the sampling site and their abundance and is sufficient for recording vegetation types. For aquatic vegetation, eDNA may be comparable with, or even superior to, in-lake vegetation surveys and may therefore be used as an tool for biomonitoring. For reconstruction of terrestrial vegetation, technical improvements and more intensive sampling is needed to detect a higher proportion of rare taxa although DNA of some taxa may never reach the lake sediments due to taphonomical constrains. Nevertheless, eDNA performs similar to conventional methods of pollen and macrofossil analyses and may therefore be an important tool for reconstruction of past vegetation.

The History of Tree and Shrub Taxa on Bol'shoy Lyakhovsky Island (New Siberian Archipelago) since the Last Interglacial Uncovered by Sedimentary Ancient DNA and Pollen Data

Ecosystem boundaries, such as the Arctic-Boreal treeline, are strongly coupled with climate and were spatially highly dynamic during past glacial-interglacial cycles. Only a few studies cover vegetation changes since the last interglacial, as most of the former landscapes are inundated and difficult to access. Using pollen analysis and sedimentary ancient DNA (sedaDNA) metabarcoding, we reveal vegetation changes on Bol'shoy Lyakhovsky Island since the last interglacial from permafrost sediments. Last interglacial samples depict high levels of floral diversity with the presence of trees (Larix, Picea, Populus) and shrubs (Alnus, Betula, Ribes, Cornus, Saliceae) on the currently treeless island. After the Last Glacial Maximum, Larix re-colonised the island but disappeared along with most shrub taxa. This was probably caused by Holocene sea-level rise, which led to increased oceanic conditions on the island. Additionally, we applied two newly developed larch-specific chloroplast markers to evaluate their potential for tracking past population dynamics from environmental samples. The novel markers were successfully re-sequenced and exhibited two variants of each marker in last interglacial samples. SedaDNA can track vegetation changes as well as genetic changes across geographic space through time and can improve our understanding of past processes that shape modern patterns.