Proxy comparison in ancient peat sediments: pollen, macrofossil and plant DNA

Using ancient sedimentary DNA to forecast ecosystem trajectories under climate change

Ecosystem response to climate change is complex. In order to forecast ecosystem dynamics, we need high-quality data on changes in past species abundance that can inform process-based models. Sedimentary ancient DNA (sedaDNA) has revolutionised our ability to document past ecosystems' dynamics. It provides time series of increased taxonomic resolution compared to microfossils (pollen, spores), and can often give species-level information, especially for past vascular plant and mammal abundances. Time series are much richer in information than contemporary spatial distribution information, which have been traditionally used to train models for predicting biodiversity and ecosystem responses to climate change. Here, we outline the potential contribution of sedaDNA to forecast ecosystem changes. We showcase how species-level time series may allow quantification of the effect of biotic interactions in ecosystem dynamics, and be used to estimate dispersal rates when a dense network of sites is available. By combining palaeo-time series, process-based models, and inverse modelling, we can recover the biotic and abiotic processes underlying ecosystem dynamics, which are traditionally very challenging to characterise. Dynamic models informed by sedaDNA can further be used to extrapolate beyond current dynamics and provide robust forecasts of ecosystem responses to future climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

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.

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.

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.

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.

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.

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.

Integrative analysis of DNA, macroscopic remains and stable isotopes of dog coprolites to reconstruct community diet

Paleofeces or coprolites are often used to reconstruct diet at archaeological sites, usually using macroscopic analyses or targeted DNA amplification and sequencing. Here we present an integrative analysis of dog coprolites, combining macroscopic analyses, stable isotope measurements, and DNA shotgun sequencing to examine diet and health status. Dog coprolites used in this study were recovered from the Janey B. Goode and East Saint Louis archaeological sites, both of which are located in the American Bottom, an extensive Mississippi River floodplain in Southwestern Illinois. Based on the context of recovery, coprolites are assigned to the Late Woodland and Terminal Late Woodland periods (ca. 600-1050 AD). Given the scarcity of human remains from this time period, these dog coprolites can be useful as a proxy for understanding human diet during the Late Woodland period. We find that the Late Woodland dogs consumed a variety of fish as well as bird and plant taxa, possibly including maize, and also harbored intestinal parasites and pathogenic bacteria. By sequencing the fecal microbiome of the coprolites, we find some similarities to modern dog microbiomes, as well as specific taxa that can be used to discriminate between modern and ancient microbiomes, excluding soil contaminants. As dogs are often used as a surrogate to assess human diet, humans living with these dogs likely had a similar diet and were affected by similar parasites. These analyses, when integrated, show a more comprehensive view of ancient dog and human diet and health in the region during the initial expansion of maize agriculture than any individual method could alone.

Multiscale climate change impacts on plant diversity in the Atacama Desert

Comprehending ecological dynamics requires not only knowledge of modern communities but also detailed reconstructions of ecosystem history. Ancient DNA (aDNA) metabarcoding allows biodiversity responses to major climatic change to be explored at different spatial and temporal scales. We extracted aDNA preserved in fossil rodent middens to reconstruct late Quaternary vegetation dynamics in the hyperarid Atacama Desert. By comparing our paleo-informed millennial record with contemporary observations of interannual variations in diversity, we show local plant communities behave differentially at different timescales. In the interannual (years to decades) time frame, only annual herbaceous expand and contract their distributional ranges (emerging from persistent seed banks) in response to precipitation, whereas perennials distribution appears to be extraordinarily resilient. In contrast, at longer timescales (thousands of years) many perennial species were displaced up to 1,000 m downslope during pluvial events. Given ongoing and future natural and anthropogenically induced climate change, our results not only provide baselines for vegetation in the Atacama Desert, but also help to inform how these and other high mountain plant communities may respond to fluctuations of climate in the future.

From microbial eukaryotes to metazoan vertebrates: Wide spectrum paleo-diversity in sedimentary ancient DNA over the last ~14,500 years

Most studies that utilize ancient DNA have focused on specific groups of organisms or even single species. Instead, the whole biodiversity of eukaryotes can be described using universal phylogenetic marker genes found within well-preserved sediment cores that cover the post-glacial period. Sedimentary ancient DNA samples from Lake Lielais Svētiņu, eastern Latvia, at a core depth of 1,050 cm in ~150 year intervals were used to determine phylotaxonomy in domain Eukaryota. Phylotaxonomic affiliation of >1,200 eukaryotic phylotypes revealed high richness in all major eukaryotic groups-Alveolata, Stramenopiles, Cercozoa, Chlorophyta, Charophyta, Nucletmycea, and Holozoa. The share of organisms that originate from terrestrial ecosystems was about one third, of which the most abundant molecular operational taxonomic units were Fungi and tracheal/vascular plants, which demonstrates the usefulness of using lake sediments to reconstruct the terrestrial paleoecosystems that surround them. Phylotypes that originate from the lake ecosystem belonged to various planktonic organisms; phyto-, proto,- and macrozooplankton, and vascular aquatic plants. We observed greater richness of several planktonic organisms that can be associated with higher trophic status during the warm climate period between 4,000 and 8,000 years ago and an increase in eukaryotic richness possibly associated with moderate human impact over the last 2,000 years.

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.

Microbial survival strategies in ancient permafrost: insights from metagenomics

In permafrost (perennially frozen ground) microbes survive oligotrophic conditions, sub-zero temperatures, low water availability and high salinity over millennia. Viable life exists in permafrost tens of thousands of years old but we know little about the metabolic and physiological adaptations to the challenges presented by life in frozen ground over geologic time. In this study we asked whether increasing age and the associated stressors drive adaptive changes in community composition and function. We conducted deep metagenomic and 16 S rRNA gene sequencing across a Pleistocene permafrost chronosequence from 19 000 to 33 000 years before present (kyr). We found that age markedly affected community composition and reduced diversity. Reconstruction of paleovegetation from metagenomic sequence suggests vegetation differences in the paleo record are not responsible for shifts in community composition and function. Rather, we observed shifts consistent with long-term survival strategies in extreme cryogenic environments. These include increased reliance on scavenging detrital biomass, horizontal gene transfer, chemotaxis, dormancy, environmental sensing and stress response. Our results identify traits that may enable survival in ancient cryoenvironments with no influx of energy or new materials.

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.

Lake sedimentary DNA accurately records 20th Century introductions of exotic conifers in Scotland

Sedimentary DNA (sedDNA) has recently emerged as a new proxy for reconstructing past vegetation, but its taphonomy, source area and representation biases need better assessment. We investigated how sedDNA in recent sediments of two small Scottish lakes reflects a major vegetation change, using well-documented 20^th Century plantations of exotic conifers as an experimental system. We used next-generation sequencing to barcode sedDNA retrieved from subrecent lake sediments. For comparison, pollen was analysed from the same samples. The sedDNA record contains 73 taxa (mainly genus or species), all but one of which are present in the study area. Pollen and sedDNA shared 35% of taxa, which partly reflects a difference in source area. More aquatic taxa were recorded in sedDNA, whereas taxa assumed to be of regional rather than local origin were recorded only as pollen. The chronology of the sediments and planting records are well aligned, and sedDNA of exotic conifers appears in high quantities with the establishment of plantations around the lakes. SedDNA recorded other changes in local vegetation that accompanied afforestation. There were no signs of DNA leaching in the sediments or DNA originating from pollen.

First physical evidence for forested environment in the Arctic during MIS 3

Old sedimentological and geochronological records can be preserved underneath the central parts of the continental ice sheets under non-erosive, cold-based subglacial conditions. Organic deposits that predate the last deglaciation are of particular value for the information held on glacial-time climate and environmental conditions. In this study, we present multiproxy data derived from a well-preserved MIS 3 interstadial (55–25 ka ago) organic layer from inside the Arctic Circle in the Finnish Lapland. Biological proxy evidence, namely coming from aquatic plant species, indicates July temperatures as high as 14.4 °C, i.e. higher than those of today for the study site. Macrofossil evidence demonstrates for the first time the presence of pines accompanied by tree birch during the MIS 3 interstadial in northern Fennoscandia. These results concur with contemporary insolation model outcomes but contradict with the previous proxy-based view of open tundra conditions during the MIS 3. The data suggest that there are highly dynamic interstadial continental ice-sheet dynamics following changes in orbital forcing. Warm climate enabled the establishment of forests on exposed landscape. Moreover, we suggest that in the light of these new data, previous MIS 3 pollen data could be re-interpreted.

How have studies of ancient DNA from sediments contributed to the reconstruction of Quaternary floras?

499 I. 499 II. 500 III. 500 IV. 500 V. 500 VI. 501 VII. 502 VIII. 504 504 References 505 SUMMARY: Ancient DNA (aDNA) from lake sediments, peats, permafrost soils, preserved megafaunal gut contents and coprolites has been used to reconstruct late-Quaternary floras. aDNA is either used alone for floristic reconstruction or compared with pollen and/or macrofossil results. In comparative studies, aDNA may complement pollen and macrofossil analyses by increasing the number of taxa found. We discuss the relative contributions of each fossil group to taxon richness and the number of unique taxa found, and situations in which aDNA has refined pollen identifications. Pressing problems in aDNA studies are contamination and ignorance about taphonomy (transportation, incorporation, and preservation in sediments). Progress requires that these problems are reduced to allow aDNA to reach its full potential contribution to reconstructions of Quaternary floras.