The reconstruction of human-driven, Earth-shaping dynamics is important for understanding past human/environment interactions and for helping human societies that currently face global changes. However, it is often challenging to distinguish the effects of the climate from human activities on environmental changes. Here we evaluate an approach based on DNA metabarcoding used on lake sediments to provide the first high-resolution reconstruction of plant cover and livestock farming history since the Neolithic Period. By comparing these data with a previous reconstruction of erosive event frequency, we show that the most intense erosion period was caused by deforestation and overgrazing by sheep and cowherds during the Late Iron Age and Roman Period. Tracking plants and domestic mammals using lake sediment DNA (lake sedDNA) is a new, promising method for tracing past human practices, and it provides a new outlook of the effects of anthropogenic factors on landscape-scale changes.

Although it is generally agreed that the arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we additionally explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr BP (before present). For much of the period investigated, arctic vegetation consisted of dry steppe tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25–15 kyr BP), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr BP, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such our findings question the predominance of a late Quaternary graminoid-dominated arctic “mammoth steppe”.

Closely related sympatric species commonly develop different ecological strategies to avoid competition. Ctenomys minutus and C. flamarioni are subterranean rodents parapatrically distributed in the southern Brazilian coastal plain, showing a narrow sympatric zone. To gain understanding on food preferences and possible competition for food resources, we evaluated their diet composition performing DNA metabarcoding analyzes of 67 C. minutus and 100 C. flamarioni scat samples, collected along the species geographical ranges. Thirteen plant families, mainly represented by Poaceae, Araliaceae, Asteraceae and Fabaceae, were identified in the diet of C. minutus. For C. flamarioni, 10 families were recovered, with a predominance of Poaceae, Araliaceae and Asteraceae. A significant correlation between diet composition and geographical distance was detected in C. minutus, whereas the diet of C. flamarioni was quite homogeneous throughout its geographical distribution. No significant differences were observed between males and females of each species. However, differences in diet composition between species were evident according to multivariate analysis. Our results suggest some level of diet partitioning between C. flamarioni and C. minutus in the sympatric region. While the first species is more specialized on few plant items, the second showed a more varied and heterogeneous diet pattern among individuals. These differences might have been developed to avoid competition in the region of co-occurrence. Resource availability in the environment also seems to influence food choices. Our data indicate that C. minutus and C. flamarioni are generalist species, but that some preference for Poaceae, Asteraceae and Araliaceae families can be suggested for both rodents.

Mosquito control is often the main method used to reduce mosquito-transmitted diseases. In order to investigate the genetic basis of resistance to the bio-insecticide Bacillus thuringiensis subsp. israelensis (Bti), we used information on polymorphism obtained from cDNA tag sequences from pooled larvae of laboratory Bti-resistant and susceptible Aedes aegypti mosquito strains to identify and analyse 1520 single nucleotide polymorphisms (SNPs). Of the 372 SNPs tested, 99.2% were validated using DNA Illumina GoldenGate® array, with a strong correlation between the allelic frequencies inferred from the pooled and individual data (r = 0.85). A total of 11 genomic regions and five candidate genes were detected using a genome scan approach. One of these candidate genes showed significant departures from neutrality in the resistant strain at sequence level. Six natural populations from Martinique Island were sequenced for the 372 tested SNPs with a high transferability (87%), and association mapping analyses detected 14 loci associated with Bti resistance, including one located in a putative receptor for Cry11 toxins. Three of these loci were also significantly differentiated between the laboratory strains, suggesting that most of the genes associated with resistance might differ between the two environments. It also suggests that common selected regions might harbour key genes for Bti resistance.

Ecosystems across the globe are threatened by climate change and human activities. New rapid survey approaches for monitoring biodiversity would greatly advance assessment and understanding of these threats. Taking advantage of next-generation DNA sequencing, we tested an approach we call metabarcoding: high-throughput and simultaneous taxa identification based on a very short (usually less than 100 base pairs) but informative DNA fragment. Short DNA fragments allow the use of degraded DNA from environmental samples. All analyses included amplification using plant-specific versatile primers, sequencing and estimation of taxonomic diversity. We tested in three steps whether degraded DNA from dead material in soil has the potential of efficiently assessing biodiversity in different biomes. First, soil DNA from eight boreal plant communities located in two different vegetation types (meadow and heath) was amplified. Plant diversity detected from boreal soil was highly consistent with plant functional and structural diversity estimated from conventional above-ground surveys. Second, we assessed DNA persistence using samples from formerly cultivated soils in temperate environments. We found that number of crop DNA sequences retrieved strongly varied with years since last cultivation, and crop sequences were absent from nearby, uncultivated plots. Third, we assessed the universal applicability of DNA metabarcoding using soil samples from tropical environments: a large proportion of species and families from the study site was efficiently recovered. The results open unprecedented opportunities for large-scale DNA-based biodiversity studies across a range of taxonomic groups using standardized metabarcoding approaches.

Metabarcoding approaches use total and typically degraded DNA from environmental samples to analyse biotic assemblages and can potentially be carried out for any kinds of organisms in an ecosystem. These analyses rely on specific markers, here called metabarcodes, which should be optimized for taxonomic resolution, minimal bias in amplification of the target organism group and short sequence length. Using bioinformatic tools, we developed metabarcodes for several groups of organisms: fungi, bryophytes, enchytraeids, beetles and birds. The ability of these metabarcodes to amplify the target groups was systematically evaluated by (1) in silico PCRs using all standard sequences in the EMBL public database as templates, (2) in vitro PCRs of DNA extracts from surface soil samples from a site in Varanger, northern Norway, and (3) in vitro PCRs of DNA extracts from permanently frozen sediment samples of late-Pleistocene age (~ 16 000–50 000 yr BP) from two Siberian sites, Duvanny Yar and Main River. Comparison of the results from the in silico PCR with those obtained in vitro showed that the in silico approach offered a reliable estimate of the suitability of a marker. All target groups were detected in the environmental DNA, but we found large variation in the level of detection among the groups and between modern and ancient samples. Success rates for the Pleistocene samples were highest for fungal DNA, whereas bryophyte, beetle and bird sequences could also be retrieved, but to a much lesser degree. The metabarcoding approach has considerable potential for biodiversity screening of modern samples and also as a paleoecological tool.

Diet analysis is a prerequisite to fully understand the biology of a species and the functioning of ecosystems. For carnivores, traditional diet analyses mostly rely upon the morphological identification of undigested remains in the feces. Here, we developed a methodology for carnivore diet analyses based on next generation sequencing. We applied this approach to the analysis of the vertebrate component of leopard cat diet in two ecologically distinct regions in northern Pakistan. Despite being a relatively common species with a wide distribution in Asia, little is known about this elusive predator. We analyzed a total of 38 leopard cat feces. After a classical DNA extraction, the DNA extracts were amplified using primers for vertebrates targeting about 100 bp of the mitochondrial 12S rRNA gene, with and without a blocking oligonucleotide specific to the predator sequence. The amplification products were then sequenced on a next generation sequencer. We identified a total of 18 prey taxa, including eight mammals, eight birds, one amphibian, and one fish. In general, our results confirmed that the leopard cat has a very eclectic diet, and feeds mainly on rodents, and particularly on the Muridae family. The DNA-based approach we propose here represents a valuable complement to current conventional methods. It can be applied to other carnivore species with only a slight adjustment relating to the design of the blocking oligonucleotide. It is robust, simple to implement, and allows the possibility of very large-scale analyses.

Earthworms are known for their important role within the functioning of an ecosystem, and their diversity can be used as an indicator of ecosystem health. To date, earthworm diversity has been investigated through conventional extraction methods such as handsorting, soil washing or the application of a mustard solution. Such techniques are time-consuming and often difficult to apply. We showed that combining DNA metabarcoding and next generation sequencing facilitates the identification of earthworm species from soil samples. The first step of our experiments was to create a reference database of mitochondrial DNA (mtDNA) 16S gene for 14 earthworm species found in the French Alps. Using this database, we designed two new primer pairs targeting very short and informative DNA sequences (about 30 bp and 70 bp) that allow unambiguous species identification. Finally, we analyzed extracellular DNA taken from soil samples in two localities (two plots per locality, eight samples per plot). The two short metabarcode regions led to the identification of a total of eight earthworm species. The earthworm communities identified by the DNA-based approach appeared to be well-differentiated between the two localities, and are consistent with results derived from inventories collected using the handsorting method. The possibility of assessing earthworm communities from hundreds or even thousands of localities through the use of extracellular soil DNA will undoubtedly stimulate further ecological research on these organisms. Using the same DNA extracts, our study also illustrates the potential of environmental DNA as a tool to assess the diversity of other soil-dwelling animal taxa.

DNA metabarcoding corresponds to the DNA-based identification of multiple species from a single complex and degraded environmental sample. We developed new sampling and extraction protocols suitable for DNA metabarcoding analyses, targeting soil extracellular DNA. The proposed sampling protocol has been designed to reduce as much as possible the influence of the local heterogeneity by processing large amount of soil, resulting from the mixing of many different cores. The DNA extraction is based on the use of saturated phosphate buffer. The sampling and extraction protocols were validated first by analyzing plant DNA from a set of 12 plots corresponding to four plant communities in alpine meadows, and second by conducting pilot experiments on fungi and earthworms. The results of the validation experiments clearly demonstrated that sound biological information can be retrieved when following these sampling and extraction procedures. Such a protocol can be implemented at any time of the year without any preliminary knowledge of specific types of organisms during the sampling. It offers the opportunity to analyze all groups of organisms using a single sampling/extraction procedure and opens the possibility to fully standardize biodiversity surveys.

Although ancient DNA from sediments (sedaDNA) has been used to investigate past ecosystems, the approach has never been directly compared to the traditional methods of pollen and macrofossil analysis. We conducted a comparative survey of 18 ancient permafrost samples spanning the Late Pleistocene (46–12.5 thousand years ago), from the Taymyr Peninsula in northern Siberia. The results show that pollen, macrofossils and sedaDNA are complementary rather than overlapping, and in combination reveal more detailed information on plant palaeocommunities than can be achieved by each individual approach. SedaDNA and macrofossils share greater overlap in plant identifications than with pollen, suggesting that sedaDNA is local in origin. These two proxies also permit identification to lower taxonomic levels than pollen, enabling investigation of temporal changes in species composition and the determination of indicator species to describe environmental changes. Combining data from all three proxies, reveals an area continually dominated by a mosaic vegetation of tundra-steppe, pioneer and wet-indicator plants. Such vegetational stability is unexpected, given the severe climate changes taking place in the northern hemisphere during this time, with changes in average annual temperatures of > 22ºC. This may explain the abundance of ice-age mammals such as horse and bison in Taymyr Peninsula during the Pleistocene, and why it acted as a refugium for the last mainland woolly mammoth. Our finding reveals the benefits of combining sedaDNA, pollen and macrofossil for palaeovegetational reconstruction and add to the increasing evidence suggesting large areas of the northern hemisphere remained ecologically stable during the Late Pleistocene.