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

Charline Giguet-Covex, Gentile Fransesco Ficetola, Kevin Walsh, Jérôme Poulenard, Manon Bajard, Laurent Fouinat, Pierre Sabatier, Ludovic Gielly, Erwan Messager, Anne Lise Develle, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Charline Giguet-Covex, Gentile Fransesco Ficetola, Kevin Walsh, Jérôme Poulenard, Manon Bajard, et al. New insights on lake sediment DNA from the catchment: importance of taphonomic and analytical issues on the record quality. Scientific Reports, 2019, 9 (14676), pp.1-21. 10.1038/s41598-019-50339-1. hal-02320430

HAL Id: hal-02320430

https://amu.hal.science/hal-02320430

Submitted on 18 Oct 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
OPEN

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

Received: 5 March 2019

Accepted: 12 August 2019

Published online: 11 October 2019

C. Giguet-Covex<sup>1,2</sup>, G. F. Ficetola<sup>3,4</sup>, K. Walsh<sup>1</sup>, J. Poulenard<sup>2</sup>, M. Bajard<sup>2</sup>, L. Fouinat<sup>2</sup>, P. Sabatier<sup>2</sup>, L. Gielly<sup>3</sup>, E. Messager<sup>2</sup>, A. L. Develle<sup>2</sup>, F. David<sup>5</sup>, P. Taberlet<sup>3</sup>, E. Brisset<sup>6,7,8</sup>, F. Guiter<sup>6</sup>, R. Sinet<sup>6</sup> & F. Arnaud<sup>2</sup>

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").

**History and potential of the lake sediment DNA (sedDNA).** The earliest studies of ancient DNA (aDNA) from lake sediment archives (lake sedDNA) date to the late-1990s<fnref n="1" />. However, it is only during the last eight years there has been the extensive application of molecular biological techniques in lake sediment analyses, with a significant contribution over recent years²⁻⁷. Lake sediments accumulate through time, as both autochthonous (in-lake biological production and chemical precipitation) and allochthonous (particles brought from the catchment and beyond) materials, which can contain DNA, and are preserved within lacustrine environments. The study of such sediments using molecular biological techniques has enormous potential for the identification of organisms present within the "lake's sediment source area" (i.e., the lake itself, its catchment area as well as the atmosphere) and thus to trace changes of biodiversity over time, from the scale of the population to that of the ecosystem. This approach allows us to address a wide range of questions, especially in ecology and about the nature of human-environment interactions over time⁴,⁸,⁹. Before 2008, there had been relatively few studies of lake sedDNA from terrestrial organisms, with all of these studies focussing on pollen DNA¹⁰,¹¹. On the contrary, many

<footnote n="1">¹BioArch-Department of Archaeology, University of York, York, YO10 5DD, UK. ²EDYTEM, UMR 5204 CNRS, Univ. Savoie Mont Blanc, Pôle Montagne, 73376, Le Bourget du Lac, France. ³Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, 38000, Grenoble, France. ⁴Department of Environmental Science and Policy, Università degli Studi di Milano, Via Celoria 26, 20133, Milano, Italy. ⁵CEREGE, UMR CNRS 7330, IRD 161-Marseille Université, Technopôle de l'Arbois Méditerranée, BP 80, 13545, Aix en Provence cedex 4, France. ⁶Aix-Marseille Univ, Avignon Univ, CNRS, IRD, IMBE, Aix-en-Provence, France. ⁷Institut Català de Paleocologia Humana i Evolució Social (IPHES), Tarragona, Spain. ⁸Área de Prehistòria, Universitat Rovira i Virgili, Tarragona, Spain. Correspondence and requests for materials should be addressed to C.G.-C. (email: charline.giguet-covex@univ-smb.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Flow chart of taphonomic and analytical processes likely to affect reconstructions of the past, especially reconstructions of landscapes and agricultural activities.</figcaption></figure>

studies have focused on aquatic organisms<sup>1,12–18</sup>. This may be due to the perception that the DNA from organisms within the lake would be accumulated in higher quantities in the sediments compared to the DNA derived from the catchment area. However, since 2008, researchers have successfully tracked organisms derived from terrestrial environments, using bulk sediments and focusing on the analysis of extracellular or total DNA<sup>19</sup>. These studies on bulk sediments targeted plants<sup>20–29</sup>, mammals<sup>24,30</sup>, humans and/or animal specific faecal bacteria<sup>30–35</sup>, and more recently other eukaryotes, such as fungi and worms<sup>5,36</sup>. The aim of such research includes the reconstruction of past biodiversity, vegetation cover, landscape, and climate changes, the nature of agro-pastoral activities and the relationships between humans and landscapes. Such studies have demonstrated the great potential of this tool in providing new knowledge for palaeoecology and archaeology. However, as discussed below, there have been several studies that have interrogated the reliability and accuracy of lake sedDNA results. To our knowledge, the data and analyses presented here constitute the first study that explicitly engages with a range of pre and post-depositional taphonomic process from a range of very different lake catchments. Here, our aim is not to reconstruct a particular ecological dynamic or environmental history but instead offer a model or framework that others might employ and develop when they plan, execute and interpret lake sedDNA research and data. While this pilot study is only based on three lakes, it represents a first step and is a timely contribution for a community increasingly using lake sedDNA.

**Issues and limits: taphonomic and analytical considerations.** *Plant DNA records.* Along with research that highlights the reliability of lake sedDNA in tracking vegetation changes, many studies have questioned the interpretation of some data, raising concerns over analytical and/or taphonomic processes (modified from<sup>37</sup>; Fig. 1). Taphonomic processes refer to all processes that govern the production, transfer, incorporation, and preservation of the lake sedDNA. For instance, Pedersen *et al.*<sup>38</sup> did not detect a substantial proportion of DNA from the local flora, which was independently identified by macrofossils. They proposed multiple, non-exclusive explanations, such as the high abundance of some taxa that may overwhelm the rarest taxa. The taxonomic resolution and assignment rate could have been limited by the degradation of DNA sequences, the sequencing depth or the incompleteness of the reference database. Indeed, in a more recent study, also from the Arctic, the authors obtained superior taxonomic recovery between data from DNA and macrofossils analyses, probably due to the use of an almost complete reference library, as well as optimised extraction protocols (sediment quantity) and sequencing conditions<sup>20</sup>. Several studies also revealed discrepancies between records of plant DNA, pollen and macrofossils, which may reflect differences in the source (production, origin), modes of transfer, and preservation conditions for these different vegetation-cover proxies<sup>22,28,38</sup>. Whereas taphonomic processes are relatively well-understood for pollen and plant macro-remains, their understanding for lake sedDNA is still limited<sup>2,8,37</sup>, especially for extracellular DNA, which by definition excludes the DNA from pollen and plant macro-remains. However, a recent review<sup>2,19</sup>, along with two studies that compared modern vegetation with pollen and DNA analyses from surface sediments from a large set of lakes in different vegetation environments (tundra to forest-tundra environments<sup>39</sup>; boreal and alpine<sup>40</sup>) have suggested the following: (1) pollen does not significantly contribute to DNA records, (2) the DNA has a local origin and probably has a similar source as the macrofossils, (3) aquatic plants are well-represented, (4) taxa detection seems to depend on the distance to lakeshore, the relief and its abundance (biomass) in the vegetation, (5) different physical and chemical sediment characteristics might have an impact on the DNA preservation. These studies targeted both intra and extracellular DNA, of which the respective contribution to the sediments remains unclear<sup>2</sup> while the taphonomic processes affecting each of these DNA pools can be expected to differ.
**Mammal DNA records.** Some studies have also raised questions regarding the taphonomic processes that affect the mammal DNA record. For example, Giguet-Covex and co-authors²⁴ did not find sheep extracellular DNA in modern sediments from a small subalpine lake (Lake Anterne, 2063 m a.s.l, Northern French Alps), despite the fact that sheep flocks are present today in the catchment. Here, low stocking-rates and scattered distributions of domestic animals, triggering a low biomass and a "diffuse source" of DNA, have been proposed as an explanation for the non-detection of DNA. On the contrary, high stocking-rates and/or the existence of areas used for the herding or flocking of animals (e.g. enclosures or folds), which represent a "point source" of DNA because of the « concentration effect » of animals, might explain the enhanced supply of mammal DNA in the sediments during earlier periods²⁴. Moreover, urine and faeces, the two primary sources of animal DNA⁴¹,⁴² - are produced mainly during the night within the enclosures or folds⁴³. The presence of enclosures within a catchment is thus expected to significantly favour the detection of domestic animal DNA. Another study that aimed to identify the presence of humans in a catchment using human-specific bacteria DNA also proposed potential biases in the record due to taphonomic issues³¹. In fact, the absence of human-specific bacteria DNA from a core where the pollen data suggests the presence of humans might be explained by DNA concentrations that were below detection limits. For example, situations where human camps/villages are at some distance from the lake or the inlet are expected to limit the DNA transfer to the lake. Moreover, situations when the population density is low are expected to limit the DNA production and biomass. An alternative explanation might also be that pollen reflects human activity from the wider region.

**“Time shifts”**. Several studies have raised the issue of potential “time shifts” in lake sediment DNA records, due to DNA leaching through sediment layers²⁸ or DNA storage in soils and its delayed release into the environment several centuries after its production²⁴. The delayed release into the environment of molecules stored in soils for several decades has already been observed for pesticides, which are persistent molecules⁴⁴. In a study of alpine soils, it was demonstrated that DNA from crops cultivated over 50 years ago could be detected⁴⁵. This study also shows a significant correlation between the proportion of DNA in soils and the proportion of above ground biomass for different functional plant groups, suggesting that the DNA transported by soil erosion will primarily reflect the ecosystem established at the time of the erosion event, and the influence of older DNA on the overall DNA signal will be minimal. This is supported by recent studies in which DNA accurately recorded the timing of changes in a vegetation cover and mammal distribution, in accordance with detailed evidence from historical and other sedimentological sources²⁸,²⁹. This high level of concordance with an independent approach highlights not only the low levels of older DNA stored in soils but also suggests limited DNA leaching through the sediment layers.

**DNA degradation/preservation processes.** DNA degradation and preservation processes within the lake water-column and sediments also have to be considered. DNA preservation and degradation are the most extensively studied taphonomic processes, as these issues are important for several research communities. These communities are interested in nutrient cycles, gene transfer, palaeoenvironmental reconstructions, or genetic studies from archaeological remains such as bones. DNA degradation is triggered by both abiotic and biotic mechanisms.

From the moment of cell death, DNA repair mechanisms cease, and DNA starts to degrade via several chemical reactions (oxidation, hydrolysis, alkylation and Maillard reaction). These processes act both inside and outside of the cells after their lysis, thus affecting both intracellular and extracellular DNA⁴⁶,⁴⁷. The rate of chemically-induced degradation is controlled by several environmental factors. Low temperature, high salt concentration (high ionic strength) and high pH, limit the hydrolysis and thus favour DNA preservation⁴⁸⁻⁵⁰. Environments protected from ultraviolet (UV) radiation also favour DNA preservation as this radiation causes DNA damage⁵⁰. The extracellular DNA is also affected by microbial activity. In fact, the degradation by DNases produced by bacteria is considered the primary mechanism responsible for extracellular DNA degradation⁵¹. However, DNA can be protected from this process when it is adsorbed onto charged surfaces (clays and humic substances), or absorbed into the crystal lattice of fine particles, amorphous crystals and particulate organic compounds⁵¹,⁵². This protection can also result from the inactivation of DNases via their binding onto particles⁵³. The binding of extracellular DNA onto particles, as well as the degree of protection, are complex processes as they are dependent on; the mineralogy of the sorbent, the presence of organic material, pH conditions, and the ionic strength and length of the DNA molecules⁵⁴,⁵⁵. In soils, nucleic acids released from cells bind quickly to particles⁵¹,⁵³,⁵⁵,⁵⁶, which delays DNA degradation and might explain the detection of a few sequences of crop DNA in the alpine soils 50 years after arable agriculture had come to an end⁴⁵. Inside the lake, bacterial activity, oxygenation, salt concentration, organic and mineral particles, UV penetration and pH conditions can vary through time and thus differentially affect the DNA preservation. When sediments are deposited in the lake bottom, they quickly become anoxic after burial, which limits microbial activity and thus favours long-term DNA preservation. However, the uppermost sediments often represent an active layer that can significantly modify the concentration and composition of DNA⁷. With burial, DNA is entirely protected from UV radiation. In marine sediments, it has also been shown that a high proportion of extracellular DNA is bound to minerals or humic substances⁵⁶⁻⁵⁸. Given the mechanisms of DNA protection afforded by binding, the absence of oxygen and UV radiation, aquatic sediments are, *a priori*, suitable environments for DNA preservation¹⁹. However, low bacterial activity, along with the binding of DNA on particles, does not prevent chemically induced DNA degradation, especially hydrolysis. DNA degradation should result in a decrease of the DNA pool with time, and a decrease in the size of DNA fragments still present within the sediment. A time-dependent DNA decrease was reported in a study of dinoflagellate DNA from fjord sediments in Antarctica⁵⁹, and several studies reported the loss of long fragments with age⁴,²⁷,⁶⁰. Ageing also triggers cytosine to thymine substitutions at the single-stranded ends of the DNA fragments, which can be used to discriminate between ancient DNA sequences and contaminations from modern DNA²⁷,⁶⁰,⁶¹. DNA
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Presentation of the study sites. (A) Location of sites (map background purchased, along with copyright, from Map Resources [Extended Use Electronic/Broadcast License]). (B) Presentation of the characteristics of each catchment-lake system (pastoral pressure, physical characteristics and plant cover).</figcaption></figure>

preservation can also vary among different groups of organisms as well as among different species of the same group⁵⁹,⁶²⁻⁶⁴.

**Challenges ahead and aim of the study.** In the light of the preceding discussion and the issues raised therein, it is apparent that there is a need to investigate the potential distortions within lake sediment DNA records due to the taphonomic processes (production, transfer, preservation of DNA) that affect DNA in sediments. Then, we are obliged to consider biases or issues with our analytical procedures (extraction/amplification/identification)⁹,⁸,³⁷ (Fig. 1). Without a good understanding of these potential biases, the potential of lake sediment DNA analyses cannot be fully realised. Of particular importance is the issue as to whether the DNA archived within sediment represents a reliable diachronic signal. More specifically, we have to ask if the following characteristics or processes are constant over time: 1) the source of DNA, 2) the processes and efficiency of DNA transfer, and 3) the preservation conditions of DNA?

Our review of the literature demonstrates that our understanding of DNA preservation processes is improving. However, few studies have focused on identifying terrestrial DNA sources and transfer processes from catchments to lakes. In order to engage with all of these processes and questions, we present the empirical analysis of temporal lake sedDNA datasets from three mountainous lake-catchment systems characterised by a range of erosion dynamics that are the product of different geological formations, topographical characteristics and vegetation and soil cover characteristics (Fig. 2A,B). Our review of these different contexts and concomitant processes allows us to elucidate complex, catchment specific, taphonomic (i.e. source and transfer) processes. Both plant and mammal extracellular DNA were investigated using the DNA *metabarcoding* approach, which is the amplification and sequencing of DNA molecules found in the environment using universal markers⁶⁵. This extracellular DNA may represent the main DNA pool in sediments⁵⁷,⁶⁶ and is of great interest as it may provide the most integrated view of aquatic, sedimentary and terrestrial biodiversity⁵⁸. Here, we only focused on this particular DNA pool to avoid the extraction of DNA from plant macro-remains, which might lead to an overrepresentation of these taxa and limit the detection of the other, rarer taxa. Sedimentological and geochemical data were also acquired, and provide information on the processes of sediment production, transfer, and deposit, as well as of lake water physicochemical conditions. Pollen or coprophilous fungi data were included in the study as complementary evidence of vegetation cover changes and domestic herd presence. All these data are essential for our understanding of the processes that drive the production and preservation of DNA records. In addition,
<table><caption>**Table 1.** Synthesis of plant DNA results for the three lakes. La Thuile and Muzelle were analysed in the same sequencing run.</caption><thead><tr><th>Lakes</th><th>La Thuile</th><th>Muzelle</th><th>Serre de l'Homme</th></tr></thead><tbody><tr><td>number of samples</td><td>50</td><td>30</td><td>41</td></tr><tr><td>replicate number performed</td><td>4</td><td>4</td><td>8</td></tr><tr><td>Illumina Hi-seq run numero</td><td>run 1</td><td>run 1</td><td>run 2</td></tr><tr><td rowspan="2">number of MOTU</td><td>Terrestrial</td><td>96</td><td>83</td><td>12</td></tr><tr><td>Aquatic</td><td>11</td><td>0</td><td>7</td></tr><tr><td rowspan="4">number of reads</td><td>Terrestrial</td><td>796266</td><td>1836110</td><td>1205395</td></tr><tr><td>Aquatic</td><td>326988</td><td>0</td><td>4517931</td></tr><tr><td>Terrestrial (%)</td><td>70,90</td><td>100</td><td>21,10</td></tr><tr><td>Aquatic (%)</td><td>29,10</td><td>0</td><td>78,90</td></tr><tr><td rowspan="10">% of samples with x positive replicates</td><td rowspan="5">Terrestrial</td><td>0</td><td>2</td><td>0</td><td>58,5</td></tr><tr><td>1</td><td>4</td><td>0,3</td><td>26,8</td></tr><tr><td>2</td><td>6</td><td>0</td><td>4,9</td></tr><tr><td>3</td><td>12</td><td>10</td><td>9,8</td></tr><tr><td>4</td><td>76</td><td>86,7</td><td>0</td></tr><tr><td>&gt;4</td><td>No analyses</td><td>No analyses</td><td>0</td></tr><tr><td rowspan="5">Aquatic</td><td>0</td><td>28</td><td>0</td><td>44</td></tr><tr><td>1</td><td>22</td><td>0</td><td>22</td></tr><tr><td>2</td><td>2</td><td>0</td><td>22</td></tr><tr><td>3</td><td>14</td><td>0</td><td>0</td></tr><tr><td>4</td><td>34</td><td>0</td><td>0</td></tr><tr><td></td><td>&gt;4</td><td>No analyses</td><td>No analyses</td><td>12</td></tr></tbody></table>

these data allow us to consider how changes in taphonomic conditions over time can affect the quality of the DNA record and thus facilitate the evaluation of our landscape and land-use reconstructions.

## Results and Interpretations

**Plant DNA detected in the three lakes.** The analysis of sequencing data was realised using appropriate bioinformatic tools (OBITOOLS software) and a filtering procedure based on the percentage of similarity with sequences in the reference database, the identification of potential contaminations and exotic taxa, and the consideration of stochastic detections. The detailed procedure is described in the “Materials and methods” section and is discussed in the Supplementary materials section 2.2. After this filtering procedure, 107 and 83 MOTUs (Molecular Operational Taxonomic Units) of plants were detected in lakes La Thuile and Muzelle, respectively, while only 19 MOTUs were found in Lake Serre de l'Homme. In Lake Muzelle, we exclusively detected DNA from terrestrial plants (100% of the reads; a read is a sequence of base pairs corresponding to a single DNA fragment and produced by a sequencer). Lake La Thuile presents a mixed recording, but most of the DNA reads are of terrestrial origin (71% of reads distributed in 96 MOTUs, Table 1). Conversely, most of the DNA reads detected in Lake Serre de l'Homme are aquatic in origin (79% of reads distributed in 7 MOTUs but probably only representing 3 different taxa, Table 1 and Supplementary Fig. 3).

Based on the comparison between the proportions of samples in which terrestrial plants are detected in 0, 1, 2, 3 or 4 replicates, it is clear that the low terrestrial plant richness detected in Lake Serre de l'Homme also corresponds to very low quantities of DNA extracted from the samples compared to the two other lakes. In fact, terrestrial plants were never detected in more than three replicates over eight. In addition, in 85% of the samples, either terrestrial plants were not detected, or they were only detected in one replicate over eight (Table 1). However, in 12% of the samples, aquatic plants were detected in more than four replicates (34% of the samples contain aquatic plant DNA in more than one replicate). Conversely to Serre de l'Homme, terrestrial plant DNA was detected in all four of the replicates in most of the samples from lakes Muzelle and La Thuile (87% and 76%, respectively) (Table 1). Some samples from Lake La Thuile also lacked plant DNA (Table 1, Fig. 3). The three lake-catchment systems are thus characterised by different plant DNA records in terms of quality. The absence of terrestrial plant DNA in most samples from Serre de l'Homme, and in some samples from La Thuile, clearly raises the question of taphonomic issues. In order to assess the potential existence of other taphonomic biases and understand their origin, it is necessary to study each record in detail, integrating other proxies for vegetation cover (pollen data) and data that characterises sedimentary dynamics (sedimentological and geochemical data).

**La Thuile: evidence for the roles of erosion processes and vegetation cover in the production of the DNA record.** The diachronic record of terrestrial plant DNA concentration (Fig. 3) can be divided into seven phases ((a) from 0 to 1000 cal. BP, (b) from 1000 to 1400 cal. BP, (c) from 1400 to 2500 cal. BP, (d) from 2500 to 3600 cal. BP, phase (e) 3600 to 4500 cal. BP, (f) from 4500 to 5200 cal. BP and (g) from 5200 to 6400 cal. BP). These phases correspond to changes in environmental conditions inferred from the sedimentological and geochemical proxies⁶⁷. In most
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Comparison between global terrestrial plant DNA and the sedimentological/geochemical properties of sediments in Lake La Thuile over the last 6500 years. To study the behaviour of land plant extracellular DNA we focused on the proxies of the richness (mean and standard deviations of the number of MOTU) and the DNA contents in the samples (mean and standard deviations of the log(number of DNA reads + 1)/ dry mass of sediment). These variables were compared to several selected sedimentological and geochemical data: the organic matter content (LOI<sub>550°C</sub>) and origin, the non-carbonate mineral matter (LOI residue), the carbonates (LOI<sub>950°C</sub>) and the total sediment flux (g/cm<sup>2</sup>/yr). The organic matter origin is determined from the combination of data from pyrolysis Rock Eval analyses (Hydrogen Index in mg HC/g TOC and Oxygen Index in mg O<sub>2</sub>/g TOC<sup>69</sup>), X-Ray fluorescence core scanner analyses (Si/Ti as a proxy of biogenic silica<sup>67</sup>), the lithological description and the aquatic plant DNA analyses (Supplementary materials, Figs 2 and 4). Seven specific phases of changes in DNA content were defined and discussed in the text (purple shaded areas **a–g**). They correspond to different sedimentological and geochemical characteristics, which inform hypotheses explaining the behaviour of the extracellular DNA from the catchment.</figcaption></figure>

of these phases (a, b, c, e and g), the terrestrial plant DNA concentration is positively correlated with the organic
matter content (Pearson correlation coefficient r = 0.82, p < 0.001 excluding phases d and f; Fig. 3). This relation-
ship probably reflects the significant role of biomass production that has been described in previous studies<sup>2,40</sup>.
However, phases (d) and (f) do not follow this pattern. They are, respectively, impoverished and enriched in DNA,
compared to the organic content. Phase (d) is also characterised by very low carbonate content (<4%) (Fig. 3),
which might indicate acidic conditions in the water column. Acidic conditions are not favourable for DNA pres-
ervation<sup>48,49,50</sup>. Moreover, our method of DNA extraction might not be efficient enough to unbind organically
(humic substances)-complexed DNA<sup>58</sup>, which might be an important pool of extracellular DNA in this part of
the sedimentary sequence mostly made of leaves and needles<sup>67</sup>. Humic substances are also known to inhibit the
PCR reaction<sup>68</sup>. The low-DNA content in phase (d) might thus be due to unfavourable preservation conditions
and/or analytical limits. Phase (f) contains as much organic matter as phase (g), but the DNA content is higher.
However, phase (f) contains much more organic matter of terrestrial origin (vs aquatic; cf Fig. 3), derived from
the erosion of forest litter and/or the direct input of the upper parts of plants within the lake<sup>67</sup>. Very high levels
of organic matter derived from forest litter are also recorded in phase (e), but the DNA content does not sig-
nificantly increase relative to phase (f). This result is probably due to the presence of humic substances and the
acidic conditions suggested by the low carbonate content, as seen in phase (d). Phase (b) has a slightly lower DNA
content than phase (a), while organic matter values are similar in these phases. Moreover, this phase presents a
very low number of taxa, especially compared to those detected by pollen analyses (Fig. 4). However, this phase is
dominated by inputs of material from deep soils, i.e. mineral soil horizons, while phase (a) is dominated by inputs
from the soil surface, i.e. organo-mineral soil horizons (Fig. 3<sup>69</sup>). Consequently, sediments are enriched in terres-
trial plant DNA when soil-surface horizons (litter and organo-mineral horizons) are transported into the lake.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Community composition of terrestrial plants provided by the DNA and pollen analyses from Lake La Thuile. The temporal evolution of the richness, the percentage of arboreal taxa, and some taxa are presented for both methods. The taxa presented here were selected in order to show examples of temporally consistent and inconsistent records, which can have different explanations. The other taxa are presented in Supplementary Fig. 8. For the richness and the percentage of arboreal taxa determined from the terrestrial plant DNA dataset, we present the mean values and standard deviations of the four replicates. The maximum richness, i.e. cumulating all the replicates is also presented as it provides a more pertinent absolute value to compare with that of pollen. For each taxon in DNA, the size of circles is proportional to the number of reads (see scale on the top of the figure). The purple shaded areas underline the periods (b,d,e), when no or very few DNA was detected.</figcaption></figure>

However, this enrichment is lower or non-existent if the lake water is acidic and/or contains humic substances; conditions that do not favour DNA preservation/recovery. On the contrary, the sediments are impoverished in terrestrial plant DNA when the erosion strongly affects the deep soil horizons. This pattern is also evidenced by a PCA analysis (see Supplementary Fig. 8). Consequently, the erosion processes (e.g. sheet erosion, gully erosion or bank undercutting), controlling the origin of the organic matter, are key processes driving terrestrial plant DNA concentration in sediments.

A comparison of the pollen and DNA records allows us to identify other potential taphonomic processes and thereby assess the ways in which these affect the vegetation cover reconstructions derived from the DNA analysis. We consider the floristic richness, the proportion of trees and shrubs, and the 34 taxa detected by both DNA and pollen analyses (sometimes associated at different taxonomic resolution; Fig. 4 and Supplementary Fig. 9). Any temporal differences between the proxies are considered a direct consequence of taphonomic biases.

Both pollen and DNA records show an increase in floristic richness from 2500 cal. BP, i.e. from phase (c) (Fig. 4). However, the results from the DNA analyses reveal a more significant increase (especially in phase (a)). In fact, before 2500 cal. BP, 31 and 11 taxa on average are detected by pollen and molecular analyses, respectively (without taking into account phases d and e of lower DNA detection). From 2500 cal. BP, the number of taxa detected with pollen increases to 34 on average for phase (c) and to 38 for phase (a), whereas with the molecular
analyses, the mean number of taxa in phases (c) and (a) reach 19 and 30 on average, respectively. Therefore, the effectiveness of plant-community detection via the DNA analyses may well be higher after 2500 cal. BP. Moreover, from 2500 to 1400 cal. BP (i.e. in phase (c)), an increase of the proportion of arboreal taxa is recorded by DNA whereas pollen data suggest deforestation. The significant increase of the erosion from 2500 cal. BP (see LOI resi- due and total flux of sediments in Fig. 3<sup>67,69</sup>), supports the argument that deforestation occurred during this phase; a possible consequence of soil instability due to human activity. Consequently, the higher detection of trees (for instance, *Quercus* sp., *Acer* sp., Betulaceae, Ulmaceae and to a lesser extent *Viburnum opulus* and *lantana*, Fig. 4 and Supplementary Fig. 9) and the increase in species richness revealed via the DNA dataset when compared with the pollen data might be explained by the higher erosion rate. In fact, erosion increases the degree of connectivity across the catchment area (i.e. creates new connexions between patches of the catchment and the hydrographic web, including the lake). On the contrary, before 2500 cal. BP, in the forested landscape, there is a probable bias towards the recording of plants that were growing on the lakeshore and the riverside, as suggested by the domi- nance riparian *Alnus* sp., (*Alnus glutinosa* and *incana*), and by the presence of *Frangula* sp. (Supplementary Fig. 9). This would be due to the transfer of proximal litter or tree leaves into the sedimentary system.

Temporal inconsistencies are recorded between *Cannabis sativa*, detected via DNA analyses, and *Cannabis sativa* or *Humulus lupulus* (from the Cannabaceae family), detected via pollen analyses (Fig. 4). The high values of these pollen, comprising 10–15% of total pollen, suggest that they originate from retting activity. In this case, both pollen and DNA are directly transferred to the lake. Consequently, high quantities of DNA from *Cannabis sativa* can be transferred to the sediments and this might explain the high detection levels during phase (b), i.e. when erosion affects the deep soil horizons and dilutes DNA inputs from other terrestrial plants (Figs 3 and 4). Conversely, in phases (a) and (c), when the erosion predominantly affects soil surface horizons, the DNA from *Cannabis sativa* may be diluted by the DNA from other plants in the catchment. As the DNA from this spe- cies becomes rarer, it competes with other more abundant DNA fragments and is therefore no longer amplified. Nevertheless, we can see that for many taxa, DNA and pollen signals are the same (excluding phases b and d). Such trends are particularly coherent for tree taxa such as *Taxus* sp., *Tilia* sp., *Abies* sp., *Alnus* sp., *Fagus* sp., *Cupressaceae* (*Juniperus* with pollen) and *Juglandaceae* (*Juglans* with pollen). Herbaceous plants, like *Rumex* sp., *Plantago* sp., *Mentha* sp./ *Mentheae*, *Helianthemum nummularium* (*Helianthemum* with pollen) and others also record the same history (Fig. 4 and Supplementary Fig. 9).

**Serre de l'Homme: evidence for the roles of topographic and physicochemical preservation conditions in DNA recording.** Very little terrestrial plant DNA (low DNA concentration and richness) is recorded in Lake Serre de l'Homme (Fig. 5). The sediments mostly comprised non-carbonate mineral matter (35.5–78%) of clastic and biogenic (diatoms) origin as well as organic matter (20.4–62%). The C:N atomic ratio fluctuates from 9.3 to 15.4, i.e. between a pure aquatic end-member and a mixed terrestrial/aquatic end-member<sup>70-74</sup> (Fig. 5). The sediments contain terrestrial plant macrofossils. The lake catchment is flat, and the "lake surface: catchment surface" ratio is high, which explains the low terrigenous inputs reflected by the low total flux of sediments (between 1 and 20 mg/cm<sup>2</sup>/yr). In these topographical conditions, only the most easily erodable materials are mobilised. These materials may be plant remains that have fallen onto the soils (constituting the source of terrestrial plant macrofossils) as well as the bare soils on sandstones (Fig. 2), which contribute to the non-carbonate mineral matter. These materi- als are not expected to bear extracellular DNA from plants, which probably contribute to poor detection levels for terrestrial plant DNA. Moreover, poor-DNA preservation conditions may be triggered by soil acidity: pH values of 4.3–5.3 have been recorded on soils developed on the same geological substratum close to the catchment. In addition, poor DNA preservation might also be explained by low water depth that engenders high temperatures and oxygenation at the lake bottom. Higher detection probability of taxa has been demonstrated in studies of deeper lakes in boreal to alpine environments in Northern Norway<sup>40</sup>. In Lake Serre de l'Homme, better in-lake preservation conditions are assumed to have existed from 300–100 cal. BP, due to the higher organic matter pro- duction that favoured the establishment of anoxic conditions, thus reducing bacterial activity. These good preser- vation conditions may contribute to the detection of high quantities of aquatic plant DNA, which is otherwise in agreement with the decrease of the C:N atomic ratio (Fig. 5).

The poor quality of the terrestrial flora reconstruction is characterised by the stochastic detection of only eight different taxa (Fig. 6). At least four of these plants live in wet environments (*Athyrium* sp., *Caltha* sp., *Saliceae* and *Filipendula ulmaria*). The proximity or good connection between these wet environments and the lake might have favoured DNA transfer of plants that grow in these environments, such as DNA from the aquatic plants<sup>28,40</sup> that appear throughout the Serre de l'Homme record (successions of *Myriophyllum* sp., *Sparganium* sp. and *Potamogeton* sp. as well as *Potamogetonaceae*, Fig. 6). On the contrary, the very poor spatial representativeness of catchment-scale flora at Serre de l'Homme probably reflects the low connectivity between the whole catchment and the lake. This is due to the absence of a well-developed hydrographic network and low erosion, both a consequence of flat topography. The role of catchment relief on catchment flora reconstructions has also been proposed in two recent studies, in Arctic and African environments<sup>40,75</sup>.

**Muzelle: evidence for the roles of glaciers and the hydrographic web in DNA recording.** The analyses of the sed- iments from Lake Muzelle reveal significant variations in terrestrial-plant DNA concentration (from 0.28 to 2.10, Fig. 7A). However, throughout the core, these sediments possess almost homogeneous concentrations of non-carbonate mineral matter (93.6% +/- 0.8) and total organic matter (4.2% +/- 0.6) and carbonates (2.2% +/- 0.4). The sedimentological dynamic of this lake is dominated by significant changes in grain size<sup>76</sup>. The concentration of terrestrial-plant DNA tends to decrease with the increase in clay content (Pearson correlation coefficient r = -0.72, p < 0.0001; Fig. 7B). These inputs of clays increase substantially during two phases: 750–625 and 310–50 cal. BP (Fig. 7A), i.e. during the Little Ice Age (LIA)<sup>77</sup>. In this context, and given the presence of a glacier in the catchment, clays are interpreted as representing a proxy for glacier sediment input (glacial flour) to
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Comparison between plant DNA (terrestrial and aquatic) and the sedimentological/geochemical properties of sediments from Lake Serre de L'Homme over the last 3800 years. To study the behaviour of plant extracellular DNA we focused on the proxies of the richness (mean and standard deviations of the number of MOTU) and the DNA content (mean and standard deviations of the log(number of DNA reads + 1)/dry mass of sediment). These variables were compared to the organic matter content (LOI$_{550°C}$) and origin (C/N atomic ratio), the content in non-carbonate mineral matter (LOI residue) and the total sediment flux (g/cm$^2$/yr). The ranges of C/N values of land plants (green shaded area), soils (brown shaded area) and algae and aquatic plants (blue shaded area) come from the literature$^{70-74}$. The main change in sediment composition is characterised by an increase in aquatic organic matter accumulation corresponding to an increase in aquatic plant DNA.</figcaption></figure>

the lake. In fact, glacier advances triggered by colder and/or wetter conditions produce more glacial flour, which increases the input of clays into the lake, especially during high precipitation events as shown by the increase of the flood frequency$^{76}$. Because these clays do not come from soils covered by plants, no extracellular DNA fragments from terrestrial plants are expected to be bound to these clays. Thus, the inputs of these DNA-free clays might dilute the DNA coming from vegetated-soil erosion and thereby explain the decreases in DNA content when clays increase (Fig. 7A).

The taxonomic richness increases significantly from 550 cal. BP, i.e. when tree-shrub cover % decreases (Fig. 8). From this period, plant communities with different ecological preferences are recorded. In fact, heathland plants, characteristic of well-developed acid soils (e.g. Vaccinium uliginosum), are detected together with plants of calcareous meadow (Myosotis alpestris), siliceous screes, snow beds or moraines (Oxyria digyna, Veronica alpina), siliceous rocks (Eritrichium sp.), calcareous rocks (Saxifraga paniculata), nutrient-rich soils (Rumex sp., most of Mentheae sp.) and wet environments (Bartsia alpina) (Fig. 8). This record of a mosaic landscape may have been favoured by the well-developed hydrographic network connecting different parts of the catchment to the lake (Fig. 2), and the important erosion dynamic as represented by the high total sediment flux (14–77 mg/cm$^2$/yr) and by the high contribution of non-carbonate mineral matter (Fig. 7). This mosaic landscape is probably the result of the landscape opening caused by the development of pastoral activities, as suggested by the presence of plants with a preference for nutrient-rich soils. Mammal DNA analyses, considered below, can help test this hypothesis.

**Mammal DNA detection and indirect evidence for pastoral activities.** Mammal DNA is only detected in the sediments from Lake La Thuile (Table 2), while herds/flocks of domestic animals currently graze
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Community composition of terrestrial and aquatic plants provided by the DNA analyses from Lake Serre de L'Homme. For each taxon, the size of circles is proportional to the number of reads (see scale on the top of the figure). Four over eight terrestrial taxa are specific of wet environments. The detection of terrestrial taxa is relatively stochastic and only three taxa are detected in more than one replicate but in one sample (*Filipendula ulmaria*, *Caltha* and *Apiaceae*). However, each aquatic taxon is more frequently detected and often in at least two replicates. Moreover, their detections are clustered in specific periods highlighted by the green areas: the periods 3800–2950 and 2250–700 cal. BP are mostly characterised by *Myriophyllum* sp., the period 700–10 cal. BP by *Sparganium* sp. and the period from 10 to –59 cal. BP the three taxa.</figcaption></figure>

on all study sites, with high pastoral pressure around lakes Serre de l'Homme and Muzelle (Fig. 2). In the first run of sequencing (four PCR replicates per sample), only cattle are detected in La Thuile (Table 2), and always in only one replicate. In the second run of sequencing (twelve PCR replicates), the number of positive replicates where mammals were detected increases to four, and we detected two additional taxa (*Ovis* sp. and *Canis* sp. in addition to *Bos* sp.) (Table 2). More mammal DNA is detected from the sediments covering the last thousand years, which correlates nicely with the detection of *Rumex* sp. (Fig. 8A; 9), a nitrophilous plant commonly associated with animal stalls. *Plantago* sp., generally associated with grazing activity because it is resistant to trampling and not eaten by animals (especially *P. alpina* and *P. Lanceolata*), is also detected in previous periods (DNA and pollen, Fig. 4), e.g. from the Late Iron Age to the Early Medieval Period. Its occurrence suggests that herds/flocks of domestic animals might have been present in the catchment during these periods, although the mammal DNA analyses do
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Comparison between terrestrial plants DNA archived in Lake Muzelle sediments and the sedimentological/geochemical properties of sediments. (A) Evolutions of the richness (mean values and standard deviations of the four replicates), the contents in DNA reads in the samples (mean number of DNA reads normalised by the dry mass of sediment and standard deviations of the four replicates), the organic matter content (LOI 550 °C), the clay content, the total sediment flux (g/cm2/yr) and the flood frequency over the last 1600 years. Blue areas highlight phases of high inputs of clays and high flood frequency, which corresponds to low DNA concentration in the sediment samples. (B) Relationship between the DNA content in the samples and the clay content.</figcaption></figure>

not reveal the presence of domesticates. This absence of animal DNA before 1000 cal. BP might be due to several factors: (1) a real low number of animals, (2) a dominance of sheep or goats relative to cattle, (3) an absence of areas of animal stalls, (4) an absence of connexion between the area of animal stalls and the hydrographic web, (5) scattered distributions of animals, (6) a predominance of deep-soil horizon erosion between 1400 and 1000 yr cal. BP (Fig. 3), and (7) a combination of these factors. The first two factors lead to low DNA quantity due to low biomass and refer to the notion of Livestock Unit (LU) in agronomy. The third, fourth and fifth factors refer to the notion of stock density (LU/surface unit/time unit) in agronomy. In another alpine lake (Anterne), sheep DNA was detected in only one of eight replicates during the Late Bronze Age, whereas *Plantago* sp. DNA was regularly recorded during this period²⁴,²⁵. In this case, the low DNA content may also be explained by a dilution triggered by the significant increase in deep soil-horizon erosion²⁴,²⁷. Furthermore, as observed at Lake La Thuile (Fig. 3), this period was also characterised by the detection of few terrestrial plant taxa²⁶.

The absence of mammal DNA in sediments from Lake Muzelle is quite unexpected. Indeed, in addition to high pastoral pressure around the lake today, DNA from *Rumex* sp. and spores from coprophilous fungi (*Sporomiella sp.*)⁷⁶ are found in the sediments dated to the last few centuries (Fig. 9B). These results strongly suggest the presence of domestic flocks/herds at least during the recent past. Coprophilous fungi spores, as well as extracellular DNA from both *Rumex* sp. and domestic animals, are supposed to share the same area of production. *Sporomiella* spores mainly come from the faeces of herbivores, mammal DNA is assumed to be largely derived from dung, and urine⁴² and DNA from *Rumex* comes from places of high nutrient accumulation, such as domestic animal stalls where faeces accumulate (hence the good correspondence with the mammal DNA observed for La Thuile). However, the production (and thus concentration) of each of these proxies, as well as their distribution through the soil profiles, can be different. Consequently, the non-detection of mammal DNA in the sediments from Lake Muzelle might be due to low production/concentration of mammal DNA compared to DNA from *Rumex* sp. as well as *Sporomiella sp. spores*. Another explanation might be the different detection limits between these distinct proxies. The difficulty of detecting mammal DNA is well illustrated by the repeated amplification of DNA from sediments from Lake La Thuile. In fact, improved detection (higher numbers of positive replicates and more taxa) of mammal DNA is attained when the number of DNA replicates is increased (Lake La Thuile Table 2 and Fig. 9A), as this increases the detection probability of “rare” taxa⁷⁹,⁸⁰. In particular, *Ovis* sp. is consistently detected in Lake La Thuile only after the performance of numerous PCR replicates (Table 2). Even if these taxa are not “rare” in the catchment, because of contaminations of samples by human DNA during the lab work (still high even with the use of blocking primers, see Supplementary Fig. 6), these taxa have to be considered as “rare” in our samples. In fact, the DNA sequences from the domestic animals are diluted by the high amount of sequences from the human contamination. Moreover, these human DNA sequences are not degraded and thus preferentially amplified during the PCR. Consequently, the low number of replicates analysed in Lake Muzelle (only four), could contribute to the non-detection of the domestic animals.
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Community composition of terrestrial plants provided by the DNA analyses from Lake Muzelle. The richness (mean and maximum), the percentage of arboreal taxa and several taxa (species and genus) of different ecological preferences (mentioned on the right side of the figure) were selected to document the landscape and environmental changes. *Alchemilla* sp. and *Potentilla* sp. can have different ecological preferences according to the species. However, these pollen types were frequently observed in overgrazed and trampling sites⁹². A study on lake sediments DNA also observed these taxa during phases when pastoral activities with sheep and/or cow were recognised²⁵. For each taxon, the size of circles is proportional to the number of reads (see scale on the top of the figure).</figcaption></figure>

<table><caption>**Table 2.** Synthesis of mammal DNA results from the three lake sediment cores.</caption><thead><tr><th rowspan="2">Primer</th><th rowspan="2">PCR replicate number</th><th rowspan="2">Illumina Hi-se run</th><th colspan="3">detected taxa</th></tr><tr><th>La Thuile</th><th>Muzelle</th><th>Serre de l'Homme</th></tr></thead><tbody><tr><td>Mam-P007</td><td>4</td><td>1</td><td>Bos sp.</td><td>No DNA</td><td>No analyses</td></tr><tr><td>Mam-P007</td><td>12</td><td>2</td><td>Bos sp., Ovis sp., (Canis sp.)</td><td>No analyses</td><td>No analyses</td></tr><tr><td>Mam-P007</td><td>8</td><td>3</td><td>No analyses</td><td>No analyses</td><td>No DNA</td></tr></tbody></table>

The absence of mammal DNA in the sediments from Lake Serre de l'Homme, where spores of *Sporomiella* sp. are also detected and 8 replicates were performed, is probably due to the low detrital input combined with poor-DNA preservation conditions as was hypothesised for terrestrial plants.
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Comparison of proxies of the presence of domestic animals in the aim of studying the taphonomic processes and analytical biases affecting mammal DNA. (A) Comparison for Lake La Thuile between the mammal DNA results obtained from the same primer “mam P007”, but not with the same replicate numbers (4 vs 12). The DNA from *Rumex sp.* is also presented as a proxy of high animal stocking rate or DNA. (B) Comparison on Lake Muzelle between the DNA from *Rumex sp.* and spores of coprophilous fungi (*Sporomiella sp.* in Fouinat *et al.*$^{76}$).</figcaption></figure>

Discussion

Based on the case studies and the literature review presented above, we propose a model summarising the archiv-
ing of the extracellular DNA from the catchment around a lake (Fig. 10). This model can be used to guide the
choice of lakes most suitable for the reconstruction of the catchment history (landscape changes, agropastoral
activities, biodiversity).

The model integrates three equations. The first is a mixing equation between the different materials affected by erosion in the catchment and transferred to the lake. This equation can be written as follows, for one taxon (Eq. 1) and several taxa (Eq. 2):

$$
[\text{DNA}_{\text{Taxa } j, \text{TERRinit}}] = \sum_{i=0}^{x_1} [\text{DNA}_{\text{Taxa } j, \text{ Source } i}] [\text{Source } i], \quad (1)
$$

where [DNA<sub>Taxa j, TERRinit</sub>] is the concentration of one taxon “j” in the different terrigenous materials affected by the erosion (log(N reads + 1)/g of terrigenous materials) and *Source i* represents the proportion of each of these terrigenous sources affected by the erosion (the number of sources varies between 0 and x1).

$$
[\text{DNA}_{\text{TERRinit}}] = \sum_{i=0}^{x_1} \left( \sum_{j=0}^{x_2} [\text{DNA}_{\text{Taxa } j, \text{ Source } i}] \right) [\text{Source } i], \quad (2)
$$

where [DNA<sub>TERRinit</sub>] corresponds to the sum of the concentrations of each taxa “j” (e.g. of plants) in the different terrigenous materials affected by the erosion (where the number of taxa “j” vary between 0 and x2). We hypothesise that the terrigenous materials contain different concentrations of DNA of each taxa j ([DNA<sub>Taxa j, Source i</sub>] due to variations in; (1) spatial distribution of the taxa in the catchment, (2) DNA distribution within soil profiles, (3) soil type, and (4) biomass produced by each taxon. For instance, according to our interpretations from Lake La Thuile, the soil litter is the richest source of plant extracellular DNA (humic substances-bound DNA; Fig. 10). However, we anticipate variations in DNA content in different types of litter (for instance, forest vs meadow), due, in particular, to different biomass production levels, litter turnover, and pH conditions. The role of these characteristics has been proposed in a study of boreal environments, although this research assessed total DNA<sup>45</sup>. Data
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Proposition of a model describing the processes driving the archiving of extracellular DNA from plants and mammals in the lake sediments. Taphonomic processes acting at the source and driving the transfer, deposit and preservation of the DNA in the lake sediments are summarised.</figcaption></figure>

from La Thuile also suggests that organo-mineral soil horizons contain less extracellular plant DNA (clay-bound DNA) than the litter, but more than mineral (deep) soil horizons. The distribution of extracellular plant DNA in soil profiles should thus have a decreasing trend from top to bottom (Fig. 10). A lower total extracellular DNA concentration has been observed in deeper horizons (B) compared to upper horizons (A) within Inceptisols (forest soils from Mediterranean regions)⁸¹. In instances of buried palaeosols⁸², higher DNA content might be expected in the “palaeo” soil surface horizon. Very acidic soils and bare soils would comprise little or no extracellular plant DNA. This may, in part, explain the poor DNA record from Lake Serre de l'Homme. Moreover, glacial flour is free of extracellular plant DNA, as exemplified by the data from Muzelle.

The content of animal extracellular DNA in soil profiles can easily differ to that derived from plants. Total DNA was shown to be strongly related to animal biomass (which is much lower than the biomass of plants) as well as to soil texture, with significant leaching in sandy soils and for larger animals⁴¹. For the livestock, this biomass in a given space and per unit of time depends on the stock density, which is driven by animal behaviour and pastoral practices (Fig. 10). These factors will also create spatial variations in mammal DNA distribution across the catchment. However, as with plants and microbes, the highest animal DNA quantities are found in top soils⁴¹.

The concentration of the different sources of terrigenous materials ([Source i]) will depend on their erodibility (capacity to be mobilised), the slope and the connections between the sources and the lake (direct or via runoff waters and tributaries). A well-developed hydrographic web should provide terrigenous inputs from the different parts of the catchment and thus afford a more reliable reconstruction of the floristic richness at the catchment scale, as exemplified by the records of a landscape mosaic in the sediments from Lake Muzelle as well as another mountain lake, Anterne²⁵. Moreover, open landscapes, with a higher erosion dynamic triggered by greater soil erodibility, should yield better spatial representativeness, such as data on the range of plants in the catchment. This process is well exemplified by the data from Lake La Thuile. However, erosion should preferentially affect the upper parts of the soils, as noted earlier. This also means that any significant developments in agricultural activities should be clearly recorded by the lake sedDNA. On the contrary, spatially extended practices, such as unmanaged grazing without stockading or animal enclosures, with less impact on the erosion dynamic, might be more difficult to detect.

Previous studies have proposed that biomass, distance from source to sink along with relief determine the terrestrial plant DNA record within sediments²,⁴⁰. Here, our model goes further, integrating, in a more explicit manner, the mechanisms behind the production and transfer of extracellular DNA in lake sediments. In fact,
our data demonstrate that the nature of erosion processes, such as sheet erosion, gully erosion, bank under-cutting or glacial erosion, should be taken into consideration as they control the sources and quantities of catchment-derived extracellular DNA inputs to the lake. Furthermore, the concept of “catchment connectivity”, combining the hydrographic web and the catchment erodibility, reflects the features and processes controlling the spatial representativeness of the DNA record, which is key for high quality, reliable reconstructions, especially when landscapes possess high habitat diversity (i.e. comprise a plant metacommunity).

The second equation of the model reflects the dilution by autochthonous production (lake production):

$$[DNA_{TERRSED}] = [DNA_{TERRinit}][TERR_{SED}] \text{ or } [DNA_{TERRinit}](1 - [AquaMat_{SED}]) \quad (3)$$

where [DNA$_{TERRSED}$] is the concentration of terrestrial DNA in the sediments (log(N reads + 1)/g dry sediments), [TERR$_{SED}$] is the concentration of terrigenous materials in the sediments (g of terrigenous materials/g of dry sediments) and [AquaMat$_{SED}$] represents the concentration of the aquatic production.

The aquatic end-member of the sediments can include organic matter from microalgae, and aquatic plants as well as mineral matter produced or induced by aquatic organisms or chemical reactions. The dilution effect by the aquatic end-member is illustrated by the records from phases (a), (c) and (g) at Lake La Thuile and probably contributes to the poor terrestrial DNA record in lake Serre de L'Homme. In the dilution equation, we did not consider atmospheric materials as these constitute very low quantities in comparison to aquatic and terrestrial materials.

Finally, the third equation integrates the DNA degradation process within the lake water column and the sediments into the model.

$$[DNA_{TERRSED}] = (1 - \alpha) \left[ \sum_{i=0}^{x1} \left( \sum_{j=0}^{x2} [DNA_{Taxa j Source i}] \right) [Source i] \right] [TERR_{SED}] \quad (4)$$

where α is a factor of degradation (if α = 1 all the DNA is degraded and if α = 0 all the DNA is preserved). Theoretically,

$$\alpha = f(pH, T^{\circ}, UV, O_2, microbial activity, salinity, sediment composition, time)$$

In the example of Lake La Thuile, we were able to recognise a probable negative impact of acidic conditions within the water column on DNA preservation (or on the capacity of our method to detect DNA due to the presence of humic substances). One hypothesis explaining DNA degradation in Lake Serre de l'Homme is the shallow water, favouring warm conditions, along with oxygenation. Interestingly, our data do not provide any clear evidence for a significant effect of DNA degradation over time. Indeed, the DNA concentration is not significantly higher towards the top of our cores, and all changes in DNA content occur abruptly and are always associated with sedimentological and/or geochemical changes.

Some of the factors influencing the quantity and the spatial representativeness of the DNA archived in the lake sediments are relatively constant over time (catchment slopes, lake surface to catchment surface ratio and the hydrographic web present during the Holocene). Therefore, these factors can be employed as an initial guide when choosing lakes that are suitable for the reconstruction of the catchment history (landscape and agropastoral activities). However, as the other factors could change over time (especially soil erodibility), a temporally continuous high-quality DNA record cannot be guaranteed. Therefore, the complex, temporally variable processes also require assessment. In fact, changes in the quality of the DNA record over time will limit reliable inter-period comparisons. This assessment is particularly important as the palaeosciences are primarily concerned with the identification and understanding of diachronic processes in socio-ecosystem trajectories, including tipping points and resilience. We demonstrate that the integration of data from sedimentary geology, geochemistry and soil science is a powerful approach for the assessment of potential taphonomic biases in DNA records. Similar approaches, integrating the context of sediment formation, should be more routinely adopted as interpretative tools.

The model that we propose is based on the study of only three lake-catchment systems. Therefore, a similar empirical field-study on modern sediments from a larger collection of lakes located across diverse geological and ecological environments, in order to avoid confounding variables, would be useful. Studies on a range of soils integrating the different soil horizons would also be informative and complementary. Moreover, there would be a need for experimental projects that recreate a series of different taphonomic scenarios. These projects will thus test and enhance the model proposed in the manuscript.

Lake sediment DNA is often considered as a biological/ecological proxy because it gives information about organisms. However, lake sediment DNA should also be considered as a bio-geological proxy for the following reasons; (1) the investigation and interpretation of the record must involve earth scientists who understand taphonomic approaches, and (2) it might be used to answer questions about the evolution of geological processes within a study area. Indeed, we feel that there is a potential to use the terrestrial DNA composition detected in lake sediments as a signature of the sources mobilised in a catchment in order to determine areas affected by erosion, today⁸³ as well as in the past.

## Material and Methods

**Regional setting and site presentation.** All three study sites are located in the French Alps, although in different ecological zones (Fig. 2A,B). The catchment of Lake La Thuile (874 m above sea level (asl)) is situated within the mountainous belt of a pre-alpine massif (the Bauges Massif, Northern French Alps). The catchment of lakes Muzelle (2105 m asl) and Serre de l'Homme (2235 m asl) are within the Ecrins massif (central French Alps),
i.e. in a more internal position relative to the alpine range. These sites are at a higher altitude than Lake La Thuile. Lake Muzelle's catchment area comprises several ecological zones/ecotones: the upper subalpine zone, the alpine zone, and the nival zone, with the presence of a relict glacier in the catchment (Fig. 2B). Serre de l'Homme is in the subalpine zone. The subalpine belt comprises the so-called "alpages" areas with high-altitude pastoral units used during the summer following the growth of grass. Given the range of altitudes covered by the sites, they cover zones that can support different types of agricultural activity. Until recently, the Lake La Thuile catchment hosted pastoral activities (including the presence of permanent farms), and multiple crops. The two other sites only support pastoral activity nowadays (Fig. 2B).

**Sites topography/geology.** Each of the catchment areas studied possesses different physical characteristics (Fig. 2B). The Lake Muzelle catchment area has the highest proportion of steep slopes of the three sites, a well-developed hydrographic network, highly erodible rocks, including schist, and partial meadow vegetation, with some bare soils exposed to erosion. The lake surface constitutes <2% of the catchment, which implies there is an important "concentration effect" of sediments derived from the catchment. Combined, these characteristics lead to significant terrigenous inputs to the lake. Furthermore, the catchment comprises a glacier. Thus, a part of these terrigenous inputs is derived from glacial erosion. This type of erosion provides glacial clayey materials ("glacial flour")⁷⁶. At Lake La Thuile, the lake surface to catchment surface ratio is 4.7%, i.e. 2.4 times higher than for Muzelle. This implies that in Lake La Thuile, the "concentration effect" is lower than in Lake Muzelle. The slopes are also less steep, the hydrographic network is poorly developed, and the vegetation cover greater (meadows, some agricultural and forested areas) than in the catchment of Lake Muzelle. However, the presence of agricultural activities triggers significant soil erosion and thus terrigenous inputs to the lake⁶⁷,⁶⁹. The physical characteristics of Serre de l'Homme's catchment are the opposite of those at Muzelle: high lake to catchment surface ratio (12.9%), gentle slopes, and no hydrographic network. These characteristics are not favourable for detrital supplies into the lake. However, rocks around the lake are easily erodible (sandstones), and there are some small barren/exposed areas (bare soils), which are susceptible to provide a few terrigenous (and more precisely clastic) inputs.

**Vegetation cover.** Around Lake La Muzelle, the vegetation cover is dominated by subalpine and alpine meadows with herbs such as grasses (Poaceae), wormwood (Artemisia), sedges (Cyperaceae) and creeping willows (Salix)⁸⁴. Lake Serre de l'Homme is surrounded by a eutrophic subalpine meadow with goosefoot (Chenopodium bonus henricus), yellow gentian (Gentiana lutea) and docks (Rumex sp.) (H. Cortot, Pers. Com.). Lake La Thuile (in a mountainous area) is surrounded by meadows and pastures. According to the exhaustive floristic survey undertaken around the lake (M. Pienne, T. Delahaye, S. Henriquet; Conservatoire Naturel de Savoie, 1999 and 2000), two types of meadows are present: a meadow with orchard grass (Dactylis glomerata) and heath false brome (Brachypodium pinnatum), which is sometimes grazed, and a mesophytic meadow dominated by grasses such as crested dogstail (Cynosurus cristatus), and ryegrass (Lolium perenne) used for grazing and mowing. Artificial grassland and kitchen garden are found in the northwest and southeast extremities of the lake. White willow (Salix alba), ashy willow (Salix cinerea), black poplar (Populus nigra), ash tree (Fraxinus excelsior) were also described at the edge of the lake. In the higher part of the catchment, there are coniferous forests comprised of spruce (Picea abies) on the north side, and of deciduous forest on the east side.

**Coring and dating.** All lake sediment cores were taken in the deepest part of the lakes, which are located approximately in the centre of the lakes (Fig. 2). For Lake La Thuile, cores were taken using a UWITEC platform and coring devices. The sediment sequence comprises two core sites. Sections from the second core are shifted by one meter in depth in order to have overlapping sections and create a continuous sequence (THU10, N45 31.813, E6 03.394. IGSN: IEFRA00BB - IGSN codes refer to an open international database. www.geosamples.org). Cores from Lake Muzelle (MUZ12, N44 57.037, E6 05.845, IGSN: IEFRA00A4) and two from Lake Serre de l'Homme (SDH-09-P1 and P2, N44 77.459, E6 23.772, IGSN: IEFRA00AW and IEFRA00AV, respectively) were taken using a UWITEC gravity corer. Core diameters are 90 mm for La Thuile and Serre de l'Homme and 63 and 90 mm for Muzelle. Another core on Lake Serre de l'Homme (SDH-1) was also taken with a Russian corer close to the shoreline. After coring, sediment cores were stored at 4°C.

The lake sediment cores used for DNA analyses as well as sedimentological/geochemical analyses measured 283.5 cm at Muzelle (core MUZ-12, 90 mm diameter from 0 to 130 cm depth and 63 mm from 130 to 183.5 cm depth), 549 cm at La Thuile (upper part of the core THU-10) and 81.5 cm (core SDH-09-P1) and 93 cm (core SDH-09-P2) at Serre de L'Homme. These cores cover different periods: 1700 years for Muzelle, 6450 years for La Thuile and 4000 years for Serre de L'Homme. Depending on the lakes, age-depth models are based on ¹⁴C dates, geomagnetic field secular variations, short-lived radionuclide measurements and known lead-pollution levels. All age-depth models were generated using the *R software* and the *R-code package 'Clam' version 2.2*⁸⁵. Details about sediment lithology and the age-depth models are provided in the Supplementary Materials (Section 1.). For Lake Serre de l'Homme, several cores were used. Thus, core correlations are also presented in detail in the "sediment lithology and dating" section of the Supplementary Materials. Age-depth models were used to estimate the sedimentation rate for each lake (cm/yr).

**Sedimentological, geochemical and microfossils analyses.** The cores were longitudinally cut, and a half-core was subsampled for DNA analyses (the heart of the slices, see section 4.6.1.) and for basic sedimentological analyses (edges of the slices). Samples reserved for DNA analyses were weighed wet. Edges of the sediment slices were weighed wet (Wet weight<sub>Edge</sub>; g) and dry (dried at 60 °C, Dry weight<sub>Edge</sub>; g) to determine the water
content (WC) and be able to calculate the total dry weight of the sediments (Dry weight<sub>Total</sub>; g) and finally the total flux of sediments (Flux<sub>Totsed</sub>; g/cm²/yr), as follow:

$$Flux_{Totsed} = (Dry weight_{Total} * Sedimentation rate)/(Half core surface * Sample thickness) \quad (5)$$

Where,

$$Dry weight_{Total} = Dryweight_{Edge} + Wet weight_{Heart} - (WC * Wetweight_{Heart});$$

and

$$WC = (Wet weight_{Edge} - Dryweight_{Edge})/Wetweight_{Edge}$$

The edge samples were then used for Loss on Ignition (LOI) analyses, except for Lake Serre de l'Homme for which the analyses were performed on another core (SDH-09-P2). Samples were firstly ground in an agate mortar, and then the standardised procedure proposed by⁸⁶ was applied. The LOI at 550 °C and then at 950 °C burns the organic matter and carbonate particles, respectively. The contributions (%) of these two components can thus be estimated. The residue of these two successive ignitions provides an estimation of the content in non-carbonate mineral matter (%) and corresponds to alumina and silica-rich particles, i.e. clastic particles and/or biogenic silica.

In Lake Muzelle, where the sediments are dominated by the mineral terrigenous fraction, grain size measurements were also undertaken at the same sampling resolution as that employed for DNA analyses (on the other half of the core). Particle size analyses were carried out on bulk sediments using a Malvern Mastersizer S, which operates on the laser diffraction principle. Only the proportion of clays (<2 µm), will be used in this study.

Complementary information about organic matter quality is used for lakes La Thuile and Serre de L'Homme (i.e. for which sediments are the richest in organic matter). In the case of Lake La Thuile, pyrolysis Rock Eval and XRF core scanner analyses from a previous study provide indices (Hydrogen Index, HI mgHC/gTOC, Oxygen Index, OI mgO₂/gTOC and Si/Ti as proxy of biogenic silica production) allowing us to distinguish the aquatic organic matter, the organic matter produced in the litter, the soil surface organo-mineral horizons, and the deep mineral soil horizons⁶⁷,⁶⁹. For Serre de l'Homme, the C/N atomic ratio was used as an indicator of aquatic organic matter and organic matter derived from soils and land plant macroremains⁷⁰,⁸⁷. The carbon (C) and nitrogen (N) contents were measured with an elemental analyser (CEREGE, Aix en Provence).

Pollen analyses from Lake La Thuile and spores of coprophilous fungi from Lake Muzelle were already published in⁶⁷ and⁷⁶, respectively. For Lake La Thuile, samples do not correspond to those used for the lake sediment DNA analyses. For Lake Muzelle, samples analysed for coprophilous fungi are the same as those for DNA.

**DNA metabarcoding approach.** Lake sediment core sub-sampling. To avoid contamination, the sampling of the three half-cores was performed in a room dedicated to sedimentological analyses at the EDYTEM laboratory (Université Savoie Mont Blanc, Le Bourget du Lac, France), where no DNA analyses were previously performed. Sediment core slices were taken using sterilised metal plates. The edges of slices were removed using sterile scalpel as the surface of the half-core was in contact with the air, and the concave edge was in contact with water that circulates along the coring tubes. For each lake, samples were cut in two parts to perform two extractions by sediment slices. Fifty, 30 and 41 samples were taken from the cores corresponding to lakes - La Thuile, Muzelle and Serre de l'Homme, respectively. The thicknesses of sediment slices are 1 cm for lakes Muzelle and Serre de l'Homme but 0.5 or 1 cm for Lake La Thuile due to substantial variations in the sedimentation rate (greater than 10-fold variations) and thus to avoid high differences in time covered by the different samples. Sample wet weights were between 2.22 and 13.04 g for Lake La Thuile, between 4.08 and 15.63 g for Lake La Muzelle and 10.49 and 23.92 g for Lake Serre de l'Homme. These significant differences are due to different water content values, particle densities (organic vs mineral) and, in cases of lakes La Thuile and Muzelle, also due to the changes in sample thickness and core diameters, respectively. In dry weights, these differences are higher because of the wide variability of the water content, especially between the top and bottom sediments (0.58 to 9.46 g for Lake La Thuile, 1.97 to 10.88 g for Lake La Muzelle and 0.76 to 14.3 g for Lake Serre de l'Homme.

**DNA extraction.** To limit artefacts and biases that can occur in metabarcoding studies, we followed strict laboratory conditions; we performed multiple controls at the different steps of laboratory work (extraction and PCR), we analysed samples in several replicates⁸⁸. DNA extractions were performed in the Laboratoire d'Ecologie Alpine (University Grenoble-Alpes, France), in a room dedicated to ancient DNA extraction. Eleven extraction controls were performed (3 for lakes Muzelle and La Thuile and 8 for Lake Serre de l'Homme).

DNA extraction was performed by mixing the sediment with 20 mL of saturated phosphate buffer (0.12 M Na₂HPO₄; pH ≈ 8) for 15 minutes. Then, the mixture was centrifuged (10 minutes at 10000 g) to recover 400 µL of the resulting supernatant. DNA was extracted from the supernatant using the NucleoSpin Soil commercial kit (Macherey-Nagel, Düren, Germany), following the manufacturer's instructions but omitting the lysis step. The DNA extract was eluted in 100 µL of SE buffer. This method of extraction allows the retrieval of the extracellular DNA pool that is dissolved in pore water and adsorbed onto mineral surfaces. It is unlikely that organically/inorganically complexed DNA is released by DNA-desorbing phosphate buffer⁵⁸.

**DNA amplification and high-throughput sequencing.** DNA amplification was realised in a second room of the ancient DNA laboratory using PCR. For the amplification of plants, we used the primers g-h, targeting the P6 loop region of the chloroplast trnL (UAA) intron⁸⁹. For the amplification of mammals, we used universal primer MamP007 amplifying 60–84 bp fragment of the mitochondrial 16S gene²⁴. To limit the
amplification of human DNA, we used a human-specific blocking oligonucleotide (MamP007_B_Hum1, 5'-GGAGCTTTAATTTATTAATGCAAACAGTACC-C3-3'). A unique combination of 8 bp long sequence of nucleotides (tag) was added at the 5'end of each primer, in order to recognise each sample after the parallel sequencing of multiple samples⁹⁰.

To improve the reliability of the detection/ non-detection pattern, we performed multiple PCR replicates on each DNA extract⁷⁹. For Lake Serre de l'Homme, we performed four PCR replicates on two DNA extraction replicates, yielding eight analyses replicates. For Muzelle and La Thuile samples, we performed four PCR replicates on one single extraction replicate using the g-h and Mam-P007 primers. For mammals in the La Thuile samples, we performed 12 additional PCR replicates per sample (33 over 50 selected samples) on a second extract obtained from the same samples (which were divided into two parts).

All DNA amplifications were carried out at a final volume of 30 µL containing 2.5 µL of DNA template. The amplification mixture contained 1 U of AmpliTaq Gold® DNA polymerase (Applied Biosystems), 15 mM Tris-HCl, 50 mM KCl, 2.5 mM MgCl₂, 0.2 mM of each dNTP, 0.1 µM of each primer and 4.8 µg of bovine serum albumin (Roche Diagnostic). We added 2 µM of the human-specific blocking oligonucleotide to the PCR mixture in mammal analyses. For all primer pairs, the PCR mixture was denatured at 95 °C for 10 minutes, followed by 45 cycles of 30 s at 95 °C also for the denaturation, 30 s at 50 °C for the hybridation and 1 min at 72 °C for the elongation. A final elongation step was applied for 7 min at 72 °C. The PCR products were then purified and mixed (equivolume mixes) before sequencing. Seventy-two PCR controls were included for each primer.

Sequencing was carried out using the Illumina Hi-seq technology (2*100 bp, paired-end reads), in three separate runs, one comprising four PCR replicates for plants and mammals from La Thuile and Muzelle samples; one for the additional 12 replicates of mammals in La Thuile samples and one for mammals and plants in Serre de l'Homme samples.

*Data treatment and representation.* The analysis of sequences and the taxonomic assignment were realised using the OBITOOLS software (http://www.grenoble.prabi.fr/trac/OBITools)⁹¹. The forward and reverse reads corresponding to the same DNA fragment were aligned and merged applying the *IlluminaPairEnd* function that takes into account the quality of merging. An "ngsfilter" file containing the list of samples and their associated combination of primer and tag was created and then used to assign each sequence to the relevant sample applying the *ngsfilter* function. Only sequences containing perfect tags and primers with a maximum of three errors were considered. The next step was to identify and merge the identical sequences for each sample using the *obiunig* function. Afterwards, the *obigrep* function allowed the filtering of sequences based on two parameters, (1) the sequence length and (2) the sequence occurrence in the entire dataset. For plants, sequences shorter than 10 bp and sequences detected less than 100 times were removed. The same filters were applied for mammals, but we only retained sequences longer than 60 bp. *Obiclean* was then used to determine the status of each sequence in each PCR product: "head", "internal" or "singleton"⁹¹. Only sequences that were more often "head" and "singleton" than "internal" in the global dataset were retained for the subsequent steps. Reference databases were built from the EMBL database with the *ecoPCR* program (gh-database-r113, mamP007-database-r113) and then used to assign a taxon to each unique sequence with the *ecoTag* function (the % of sequence similarity was calculated and specified in the final file).

For the subsequent analyses, only the sequences with a similarity >95% to taxa in the reference database were selected. We considered a sequence as present in a PCR replicate when at least five reads were counted²⁵. In each lake dataset, we did not consider taxa that were only detected in one sample, or stochastically in less than two replicates (i.e. taxa always detected in only one replicate but with detections in consecutive samples were kept). To remove contaminants, we excluded taxa frequently present in extraction and PCR negative controls (in more than 5 controls, where the total number of reads was greater than 10000), and taxa allochthonous in the Alps (like *Actinidia* sp.) (see Supplementary Section 2.1 as well as Supplementary Figs 3, 6 and Table 2 for more details on contamination and on the data filtering steps). Potential impacts of the filtering procedure on the main results of the study are also presented and discussed in the Supplementary Materials (Supplementary Section 2.2. and Supplementary Figs 4 and 5).

For each PCR replicate, we summed the total number of reads corresponding to terrestrial plants, aquatic plants and mammals separately. Then, we determined the mean and standard deviation of the log-transformed total number of reads across PCR replicates, as well as the number of replicates where more than 20 reads were detected. These two parameters are positively correlated (see Supplementary Section 3), which supports the assumption that the number of reads is correlated to the DNA quantity available for amplification as suggested by previous studies on soils and lake sediments²⁸,⁴⁵. We normalised the log-transformed number of reads by the dry weight of sediments used for the extractions in order to obtain a proxy of the DNA concentration that we can compare with the concentrations of the main sediment components. The log-transformation helps to correct the exponential DNA amplification during the PCR. We also determined a proxy of the richness (number of MOTUs: Molecular Operational Taxonomic Units) of mammals and plants, considering the presence of the taxa (more than 5 reads). As part of this process, for terrestrial plants, the mean value and standard deviation across replicates were calculated. We also determined a "maximum richness" from the sum of reads obtained in all the replicates for each detected taxon. For Lake La Thuile, we also calculated the pollen taxon richness to compare it with the proxy of the plant DNA richness, as that had already been carried out on another lake, but with plant macroremain data²⁰. For mammals, we only determined the maximum richness from the sum of reads obtained in all the replicates for each detected taxon.
Data Availability

Sequences for plant and mammal DNA (obitools output), post-filtering plant and mammal DNA datasets and sedimentological/geochemical data are deposited and available in the figshare database (10.6084/m9. figshare.9792644).

References

1. Coolen, M. J. L. & Overmann, J. R. Analysis of Subfossil Molecular Remains of Purple Sulfur Bacteria in a Lake Sediment. *APPL ENV. MICROBIOL* **64**, 9 (1998).

2. Parducci, L. *et al*. Ancient plant DNA in lake sediments. *New Phytol*. **214**, 924–942 (2017).

3. Ahmed, E. *et al*. Archaeal community changes in Lateglacial lake sediments: Evidence from ancient DNA. *Quat. Sci. Rev.* **181**, 19–29 (2018).

4. Domaizon, I., Winegardner, A., Capo, E., Gauthier, J. & Gregory-Eaves, I. DNA-based methods in paleolimnology: new opportunities for investigating long-term dynamics of lacustrine biodiversity. *J. Paleolimnol*. **58**, 1–21 (2017).

5. Kisand, V. *et al*. From microbial eukaryotes to metazoan vertebrates: Wide spectrum paleo-diversity in sedimentary ancient DNA over the last ~14,500 years. *Geobiology* **16**, 628–639 (2018).

6. Olajos, F. *et al*. Estimating species colonization dates using DNA in lake sediment. *Methods Ecol. Evol.* **9**, 535–543 (2018).

7. Vuillemin, A. *et al*. Preservation and Significance of Extracellular DNA in Ferruginous Sediments from Lake Towuti, Indonesia. *Front. Microbiol.* **8** (2017).

8. Bálint, M. *et al*. Environmental DNA Time Series in Ecology. *Trends Ecol. Evol.* **33**, 945–957 (2018).

9. Seddon, A. W. R. *et al*. Looking forward through the past: identification of 50 priority research questions in palaeoecology. *J. Ecol.* **102**, 256–267 (2014).

10. Bennett, K. D. & Parducci, L. DNA from pollen: principles and potential. *The Holocene* **16**, 1031–1034 (2006).

11. Parducci, L., Suyama, Y., Lascoux, M. & Bennett, K. D. Ancient DNA from pollen: a genetic record of population history in Scots pine. *Mol. Ecol.* **14**, 2873–2882 (2005).

12. Bissett, A., Gibson, J. A. E., Jarman, S. N., Swadling, K. M. & Cromer, L. Isolation, amplification, and identification of ancient copepod DNA from lake sediments: Ancient copepod DNA. *Limnol. Oceanogr. Methods* **3**, 533–542 (2005).

13. Coolen, M. Combined DNA and lipid analyses of sediments reveal changes in Holocene haptophyte and diatom populations in an Antarctic lake. *Earth Planet. Sci. Lett.* **223**, 225–239 (2004).

14. Coolen, M. J. L. *et al*. Evolution of the methane cycle in Ace Lake (Antarctica) during the Holocene: response of methanogens and methanotrophs to environmental change. *Org. Geochem.* **35**, 1151–1167 (2004).

15. D'Andrea, W. J. *et al*. Alkenone producers inferred from well-preserved 18S rDNA in Greenland lake sediments. *J. Geophys. Res.* **111** (2006).

16. Limburg, P. A. & Weider, L. J. 'Ancient' DNA in the resting egg bank of a microcrustacean can serve as a palaeolimnological database. *Proc. R. Soc. Lond. B Biol. Sci.* **269**, 281–287 (2002).

17. Marková, S., Černý, M., Rees, D. & Stuchlík, E. Are they still viable? Physical conditions and abundance of Daphnia pulicaria resting eggs in sediment cores from lakes in the Tatra Mountains. *Biologia (Bratisl.).* **61** (2006).

18. Reid, V. A., Carvalho, G. R. & George, D. G. Molecular genetic analysis of Daphnia in the English Lake District: species identity, hybridisation and resting egg banks. *Freshw. Biol.* **44**, 247–253 (2000).

19. Parducci, L., Nota, K. & Wood, J. Reconstructing Past Vegetation Communities Using Ancient DNA from Lake Sediments. In, https://doi.org/10.1007/13836_2018_38 (Springer International Publishing, 2018).

20. Alsos, I. G. *et al*. Sedimentary ancient DNA from Lake Skartjøra, Svalbard: Assessing the resilience of arctic flora to Holocene climate change. *The Holocene* **26**, 627–642 (2016).

21. Anderson-Carpenter, L. L. *et al*. Ancient DNA from lake sediments: Bridging the gap between paleoecology and genetics. *BMC Evol. Biol.* **11** (2011).

22. Boessenkool, S. *et al*. Use of Ancient Sedimentary DNA as a Novel Conservation Tool for High-Altitude Tropical Biodiversity: Reconstructing Afro-Alpine Biodiversity. *Conserv. Biol.* **28**, 446–455 (2014).

23. Epp, L. S. *et al*. Lake sediment multi-taxon DNA from North Greenland records early post-glacial appearance of vascular plants and accurately tracks environmental changes. *Quat. Sci. Rev.* **117**, 152–163 (2015).

24. Giguet-Covex, C. *et al*. Long livestock farming history and human landscape shaping revealed by lake sediment DNA. *Nat. Commun.* **5** (2014).

25. Pansu, J. *et al*. Reconstructing long-term human impacts on plant communities: an ecological approach based on lake sediment DNA. *Mol. Ecol.* **24**, 1485–1498 (2015).

26. Parducci, L. *et al*. Molecular- and pollen-based vegetation analysis in lake sediments from central Scandinavia. *Mol. Ecol.* **22**, 3511–3524 (2013).

27. Pedersen, M. W. *et al*. Postglacial viability and colonization in North America's ice-free corridor. *Nature* **537**, 45–49 (2016).

28. Sjögren, P. *et al*. Lake sedimentary DNA accurately records 20th Century introductions of exotic conifers in Scotland. *New Phytol.* **213**, 929–941 (2017).

29. Ficetola, G. F. *et al*. DNA from lake sediments reveals long-term ecosystem changes after a biological invasion. *Sci. Adv.* **4**, eaar4292 (2018).

30. Etienne, D. *et al*. Two thousand-year reconstruction of livestock production intensity in France using sediment-archived fecal Bacteroidales and source-specific mitochondrial markers. *The Holocene* **25**, 1384–1393 (2015).

31. Madeja, J. A new tool to trace past human presence from lake sediments: the human-specific molecular marker *Bacteroides* strain HF 183: Past human presence from lake sediments. *J. Quat. Sci.* **30**, 349–354 (2015).

32. Madeja, J., Wacnik, A., Wypasek, E., Chandran, A. & Stankiewicz, E. Integrated palynological and molecular analyses of late Holocene deposits from Lake Milkowskie (NE Poland): Verification of local human impact on environment. *Quat. Int.* **220**, 147–152 (2010).

33. Madeja, J. *et al*. Bacterial ancient DNA as an indicator of human presence in the past: its correlation with palynological and archaeological data. *J. Quat. Sci.* **24**, 317–321 (2009).

34. Mahaney, W. C. *et al*. Biostratigraphic Evidence Relating to the Age-Old Question of Hannibal's Invasion of Italy, II: Chemical Biomarkers and Microbial Signatures: *Hannibal's invasion of Italy, II. Archaeometry* **59**, 179–190 (2017).

35. Matisoo-Smith, E. *et al*. Recovery of DNA and pollen from New Zealand lake sediments. *Quat. Int.* **184**, 139–149 (2008).

36. Lammers, Y. *et al*. Clitellate worms (Annelida) in lateglacial and Holocene sedimentary DNA records from the Polar Urals and northern Norway. *Boreas* **48**, 317–329 (2018).

37. Birks, H. J. B. & Birks, H. H. How have studies of ancient DNA from sediments contributed to the reconstruction of Quaternary floras? *New Phytol.* **209**, 499–506 (2016).

38. Pedersen, M. W. *et al*. A comparative study of ancient environmental DNA to pollen and macrofossils from lake sediments reveals taxonomic overlap and additional plant taxa. *Quat. Sci. Rev.* **75**, 161–168 (2013).

39. Niemeyer, B., Epp, L. S., Stoof-Leichsenring, K. R., Pestrakova, L. A. & Herzschuh, U. A comparison of sedimentary DNA and pollen from lake sediments in recording vegetation composition at the Siberian treeline. *Mol. Ecol. Resour.* **17**, e46–e62 (2017).
40. Alsos, I. G. et al. Metabarcoding lake sediments: taphonomy and representation of contemporary vegetation in environmental DNA (eDNA) records. *bioRxiv*, https://doi.org/10.1101/264903 (2018).

41. Andersen, K. et al. Meta-barcoding of 'dirt' DNA from soil reflects vertebrate biodiversity: META-BARCODING OF 'DIRT' DNA FROM SOIL. *Mol. Ecol.* **21**, 1966-1979 (2012).

42. Lydolph, M. C. et al. Beringian Paleoecology Inferred from Permafrost-Preserved Fungal DNA. *Appl. Environ. Microbiol.* **71**, 1012–1017 (2005).

43. Dorioz, J.-M. Alpages, prairies et pâturages d'altitude: l'exemple du Beaufortain. *Courr. Environ. INRA Paris Inst. Natl. Rech. Agron. Délégation Perm. À Environ.* **35**, 33-42 (1998).

44. Sabatier, P. et al. Long-term relationships among pesticide applications, mobility, and soil erosion in a vineyard watershed. *Proc. Natl. Acad. Sci.* **111**, 15647–15652 (2014).

45. Yoccoz, N. G. et al. DNA from soil mirrors plant taxonomic and growth form diversity: DNA FROM SOIL MIRRORS PLANT DIVERSITY. *Mol. Ecol.* **21**, 3647–3655 (2012).

46. Hofreiter, M., Serre, D., Poinar, H. N., Kuch, M. & Pääbo, S. Ancient DNA. *Nat. Rev. Genet.* **2**, 353 (2001).

47. Willerslev, E. & Cooper, A. Review Paper. Ancient DNA. *Proc. R. Soc. B Biol. Sci.* **272**, 3–16 (2005).

48. Lindahl, T. & Nyberg, B. Heat-induced deamination of cytosine residues in deoxyribonucleic acid. *Biochemistry* **13**, 3405–3410 (1974).

49. Lindahl, T. & Nyberg, B. Rate of depurination of native deoxyribonucleic acid. *Biochemistry* **11**, 3610–3618 (1972).

50. Strickler, K. M., Fremier, A. K. & Goldberg, C. S. Quantifying effects of UV-B, temperature, and pH on eDNA degradation in aquatic microcosms. *Biol. Conserv.* **183**, 85–92 (2015).

51. Blum, S. A. E., Lorenz, M. G. & Wackernagel, W. Mechanism of Retarded DNA Degradation and Prokaryotic Origin of DNases in Nonsterile Soils. *Syst. Appl. Microbiol.* **20**, 513–521 (1997).

52. Levy-Booth, D. J. et al. Cycling of extracellular DNA in the soil environment. *Soil Biol. Biochem.* **39**, 2977–2991 (2007).

53. Pietramellara, G. et al. Extracellular DNA in soil and sediment: fate and ecological relevance. *Biol. Fertil. Soils* **45**, 219–235 (2009).

54. Ogram, A., Sayler, G. S., Gustin, D. & Lewis, R. J. DNA adsorption to soils and sediments. *Environ. Sci. Technol.* **22**, 982–984 (1988).

55. Crecchio, C. & Stotzky, G. Binding of DNA on humic acids: Effect on transformation of Bacillus subtilis and resistance to DNase. *Soil Biol. Biochem.* **30**, 1061–1067 (1998).

56. Lorenz, M. G. & Wackernagel, W. Bacterial Gene Transfer by Natural Genetic Transformation in the Environment. *MICROBIOL REV* **58**, 40 (1994).

57. Dell'Anno, A., Stefano, B. & Danovaro, R. Quantification, base composition, and fate of extracellular DNA in marine sediments. *Limnol. Oceanogr.* **47**, 899–905 (2002).

58. Torti, A., Lever, M. A. & Jørgensen, B. B. Origin, dynamics, and implications of extracellular DNA pools in marine sediments. *Mar. Genomics* **24**, 185–196 (2015).

59. Boere, A. C. et al. Source-specific variability in post-depositional DNA preservation with potential implications for DNA based paleoecological records. *Organic Geochemistry* **42**, 1216–1225 (2011).

60. Capo, E. et al. Long-term dynamics in microbial eukaryotes communities: a palaeolimnological view based on sedimentary DNA. *Mol. Ecol.* **25**, 5925–5943 (2016).

61. Slon, V. et al. Neandertal and Denisovan DNA from Pleistocene sediments. *Science* **356**, 605–608 (2017).

62. Capo, E. et al. Is Planktonic Diversity Well Recorded in Sedimentary DNA? Toward the Reconstruction of Past Protistan Diversity. *Microb. Ecol.* **70**, 865–875 (2015).

63. Epp, L. S. et al. New environmental metabarcodes for analysing soil DNA: potential for studying past and present ecosystems: METABARCODES TO ANALYSE SOIL DNA. *Mol. Ecol.* **21**, 1821–1833 (2012).

64. Pedersen, M. W. et al. Ancient and modern environmental DNA. *Philos. Trans. R. Soc. B Biol. Sci.* **370**, 20130383 (2015).

65. Taberlet, P., Coissac, E., Pompanon, F., Brochmann, C. & Willerslev, E. Towards next-generation biodiversity assessment using DNA metabarcoding: NEXT-GENERATION DNA METABARCODING. *Mol. Ecol.* **21**, 2045–2050 (2012).

66. Corinaldesi, C., Danovaro, R. & Dell'Anno, A. Simultaneous Recovery of Extracellular and Intracellular DNA Suitable for Molecular Studies from Marine Sediments. *Appl. Environ. Microbiol.* **71**, 46–50 (2005).

67. Bajard, M. et al. Erosion record in Lake La Thuile sediments (Prealps, France): Evidence of montane landscape dynamics throughout the Holocene. *The Holocene* **26**, 350–364 (2016).

68. Albers, C. N., Jensen, A., Bælum, J. & Jacobsen, C. S. Inhibition of DNA Polymerases Used in Q-PCR by Structurally Different Soil-Derived Humic Substances. *Geomicrobiol. J.* **30**, 675–681 (2013).

69. Bajard, M. et al. Progressive and regressive soil evolution phases in the Anthropocene. *CATENA* **150**, 39–52 (2017).

70. Bertrand, S. et al. Bulk organic geochemistry of sediments from Puyehue Lake and its watershed (Chile, 40°S): Implications for paleoenvironmental reconstructions. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **294**, 56–71 (2010).

71. Duarte, C. M. Nutrient concentration of aquatic plants: Patterns across species. *Limnol. Oceanogr.* **37**, 882–889 (1992).

72. Li, W. et al. Effects of water depth on carbon, nitrogen and phosphorus stoichiometry of five submersed macrophytes in an *in situ* experiment. *Ecol. Eng.* **61**, 358–365 (2013).

73. Meyers, P. A. Organic geochemical proxies of paleoceanographic, paleolimnologic, and paleoclimatic processes. *Org. Geochem.* **27**, 213–250 (1997).

74. Thevenon, F., Adatte, T., Spangenberg, J. E. & Anselmetti, F. S. Elemental (C/N ratios) and isotopic (δ15N, δ13Corg) compositions of sedimentary organic matter from a high-altitude mountain lake (Meidsee, 2661 m a.s.l., Switzerland): Implications for Lateglacial and Holocene Alpine landscape evolution. *The Holocene* **22**, 1135–1142 (2012).

75. Bremond, L. et al. Five thousand years of tropical lake sediment DNA records from Benin. *Quat. Sci. Rev.* **170**, 203–211 (2017).

76. Fouinat, L. et al. One thousand seven hundred years of interaction between glacial activity and flood frequency in proglacial Lake Muzelle (western French Alps). *Quat. Res.* **87**, 407–422 (2017).

77. Buntgen, U. et al. 2500 Years of European Climate Variability and Human Susceptibility. *Science* **331**, 578–582 (2011).

78. Giguet-Covex, C. et al. Changes in erosion patterns during the Holocene in a currently treeless subalpine catchment inferred from lake sediment geochemistry (Lake Anterne, 2063 m a.s.l., NW French Alps): The role of climate and human activities. *The Holocene* **21**, 651–665 (2011).

79. Ficetola, G. F. et al. Replication levels, false presences and the estimation of the presence/absence from eDNA metabarcoding data. *Mol. Ecol. Resour.* **15**, 543–556 (2015).

80. Furlan, E. M., Gleeson, D., Hardy, C. M. & Duncan, R. P. A framework for estimating the sensitivity of eDNA surveys. *Mol. Ecol. Resour.* **16**, 641–654 (2016).

81. Agnelli, A. et al. Purification and isotopic signatures (δ13C, δ15N, Δ14C) of soil extracellular DNA. *Biol. Fertil. Soils* **44**, 353–361 (2007).

82. Bajard, M. et al. Long-term changes in alpine pedogenetic processes: Effect of millennial agro-pastoralism activities (French-Italian Alps). *Geoderma* **306**, 217–236 (2017).

83. Evrard, O. et al. Environmental DNA provides information on sediment sources: A study in catchments affected by Fukushima radioactive fallout. *Sci. Total Environ.* **665**, 873–881 (2019).

84. Coûteaux, M. Fluctuations glaciaires de la fin du Würm dans les Alpes françaises, établies par des analyses polliniques. *Boreas* **12**, 35–56 (1983).

85. Blaauw, M. Methods and code for 'classical' age-modelling of radiocarbon sequences. *Quat. Geochronol.* **5**, 512–518 (2010).
86. Heiri, O., Lotter, A. F. & Lemcke, G. Loss on ignition as a method for estimating organic and carbonate content in sediments: reproducibility and comparability of results. *J. Paleolimnol.* **25**, 101–110 (2001).

87. Meyers, P. A. & Ishiwatar, R. Lacustrine organic geochemistryman overview of indicators of organic matter sources and diagenesis in lake sediments. *Org. Geochem.* **20**, 867–900 (1993).

88. Zinger, L. et al. DNA metabarcoding - need for robust experimental designs to draw sound ecological conclusions. *Mol. Ecol*. In press.

89. Taberlet, P. et al. Power and limitations of the chloroplast trnL (UAA) intron for plant DNA barcoding. *Nucleic Acids Res.* **35**, e14–e14 (2007).

90. Binladen, J. et al. The Use of Coded PCR Primers Enables High-Throughput Sequencing of Multiple Homolog Amplification Products by 454 Parallel Sequencing. *PLoS ONE* **2**, e197 (2007).

91. Boyer, F. et al. OBITOOLS:aUNIX-inspired software package for DNAmetabarcoding. *Mol. Ecol. Resour.* **16**, 176–182 (2016).

92. Court-Picon, M., Buttler, A. & de Beaulieu, J.-L. Modern pollen–vegetation relationships in the Champsaur valley (French Alps) and their potential in the interpretation of fossil pollen records of past cultural landscapes. *Rev. Paleobot. Palyno.* **135**, 13–39 (2005).

## Acknowledgements

This work was realised in the framework of several programs funded by the following organisms: DiPEE Chambéry-Grenoble (University Grenoble Alpes and University Savoie Mont Blanc), the Zone Atelier Alpes (ZAA), the PEPS de site Grenoble and the Conseil Général of Isère (AAP on biodiversity). These researches were also supported by the Parc National des Ecrins to obtain the funding as well as for the lakes coring. C. Giguet-Covex's salary was first supported by the University of Grenoble and the University Savoie Mont Blanc and then by the European Commission in the framework of a Marie Sklodowska Curie fellowship (PALEO-AGRI project, EC grant number 655331). We have also benefited of previous researches led by the pluridisciplinary team of the FRECHALP research program and conducted by F. Mocci. We would also like to thank Fransesco Carrer for the production of the maps of the lake catchments with their slope distributions.

## Author Contributions

C.G.-C., J.P. and F.A. and K.J.W., contributed to the concept and designed the study. C.G.-C., G.F.F., P.T. and L.G. performed the DNA analyses and bioinformatic treatments. M.B., L.F., A.-L.D., P.S., E.B., R.S., F.G. and F.D. created the sedimentological, geochemical and pollen datasets. C.G.-C. analysed the data with the help of G.F.F. and P.S. F.A. and J.P. contributed their expertise on the reconstruction of soil erosion and dynamics. E.M. provided expertise on the reconstructions of plant cover based on pollen analyses. K.J.W. provided its expertise on taphonomic processes in archaeological contexts and edited the text. C.G.-C. wrote the manuscript with the contributions of all co-authors. C.G.-C. and K.J.W. coordinated the Marie Sklodowska Curie Individual fellowship (PALEO-AGRI project, EC Grant Number 655331) and the program funded by the Parc national des Ecrins. C.G.-C. and P.T. coordinated the Amaryllis program funded by the Conseil Général of Isère.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-50339-1.

**Competing Interests:** L.G. and P.T. are co-inventors of patents related to the gh primers and the use of the P6 loop of the chloroplast trnL (UAA) intron for plant identification using degraded template DNA. These patents only restrict commercial applications and have no impact on the use of this locus by academic researchers. The other authors declare no competing interests.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_13.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019