Dynamics and drivers of mycorrhizal fungi after glacier retreat

The development of terrestrial ecosystems depends greatly on plant mutualists such as mycorrhizal fungi. The global retreat of glaciers exposes nutrient-poor substrates in extreme environments and provides a unique opportunity to study early successions of mycorrhizal fungi by assessing their dynamics and drivers. We combined environmental DNA metabarcoding and measurements of local conditions to assess the succession of mycorrhizal communities during soil development in 46 glacier forelands around the globe, testing whether dynamics and drivers differ between mycorrhizal types. Mycorrhizal fungi colonized deglaciated areas very quickly (< 10 yr), with arbuscular mycorrhizal fungi tending to become more diverse through time compared to ectomycorrhizal fungi. Both alpha- and beta-diversity of arbuscular mycorrhizal fungi were significantly related to time since glacier retreat and plant communities, while microclimate and primary productivity were more important for ectomycorrhizal fungi. The richness and composition of mycorrhizal communities were also significantly explained by soil chemistry, highlighting the importance of microhabitat for community dynamics. The acceleration of ice melt and the modifications of microclimate forecasted by climate change scenarios are expected to impact the diversity of mycorrhizal partners. These changes could alter the interactions underlying biotic colonization and belowground-aboveground linkages, with multifaceted impacts on soil development and associated ecological processes.

Transient social-ecological dynamics reveal signals of decoupling in a highly disturbed Anthropocene landscape

Understanding the transient dynamics of interlinked social-ecological systems (SES) is imperative for assessing sustainability in the Anthropocene. However, how to identify critical transitions in real-world SES remains a formidable challenge. In this study, we present an evolutionary framework to characterize these dynamics over an extended historical timeline. Our approach leverages multidecadal rates of change in socioeconomic data, paleoenvironmental, and cutting-edge sedimentary ancient DNA records from China's Yangtze River Delta, one of the most densely populated and intensively modified landscapes on Earth. Our analysis reveals two significant social-ecological transitions characterized by contrasting interactions and feedback spanning several centuries. Initially, the regional SES exhibited a loosely connected and ecologically sustainable regime. Nevertheless, starting in the 1950s, an increasingly interconnected regime emerged, ultimately resulting in the crossing of tipping points and an unprecedented acceleration in soil erosion, water eutrophication, and ecosystem degradation. Remarkably, the second transition occurring around the 2000s, featured a notable decoupling of socioeconomic development from ecoenvironmental degradation. This decoupling phenomenon signifies a more desirable reconfiguration of the regional SES, furnishing essential insights not only for the Yangtze River Basin but also for regions worldwide grappling with similar sustainability challenges. Our extensive multidecadal empirical investigation underscores the value of coevolutionary approaches in understanding and addressing social-ecological system dynamics.

The importance of species addition 'versus' replacement varies over succession in plant communities after glacier retreat

The mechanisms underlying plant succession remain highly debated. Due to the local scope of most studies, we lack a global quantification of the relative importance of species addition 'versus' replacement. We assessed the role of these processes in the variation (β-diversity) of plant communities colonizing the forelands of 46 retreating glaciers worldwide, using both environmental DNA and traditional surveys. Our findings indicate that addition and replacement concur in determining community changes in deglaciated sites, but their relative importance varied over time. Taxa addition dominated immediately after glacier retreat, as expected in harsh environments, while replacement became more important for late-successional communities. These changes were aligned with total β-diversity changes, which were more pronounced between early-successional communities than between late-successional communities (>50 yr since glacier retreat). Despite the complexity of community assembly during plant succession, the observed global pattern suggests a generalized shift from the dominance of facilitation and/or stochastic processes in early-successional communities to a predominance of competition later on.

Local climate modulates the development of soil nematode communities after glacier retreat

The worldwide retreat of glaciers is causing a faster than ever increase in ice-free areas that are leading to the emergence of new ecosystems. Understanding the dynamics of these environments is critical to predicting the consequences of climate change on mountains and at high latitudes. Climatic differences between regions of the world could modulate the emergence of biodiversity and functionality after glacier retreat, yet global tests of this hypothesis are lacking. Nematodes are the most abundant soil animals, with keystone roles in ecosystem functioning, but the lack of global-scale studies limits our understanding of how the taxonomic and functional diversity of nematodes changes during the colonization of proglacial landscapes. We used environmental DNA metabarcoding to characterize nematode communities of 48 glacier forelands from five continents. We assessed how different facets of biodiversity change with the age of deglaciated terrains and tested the hypothesis that colonization patterns are different across forelands with different climatic conditions. Nematodes colonized ice-free areas almost immediately. Both taxonomic and functional richness quickly increased over time, but the increase in nematode diversity was modulated by climate, so that colonization started earlier in forelands with mild summer temperatures. Colder forelands initially hosted poor communities, but the colonization rate then accelerated, eventually leveling biodiversity differences between climatic regimes in the long term. Immediately after glacier retreat, communities were dominated by colonizer taxa with short generation time and r-ecological strategy but community composition shifted through time, with increased frequency of more persister taxa with K-ecological strategy. These changes mostly occurred through the addition of new traits instead of their replacement during succession. The effects of local climate on nematode colonization led to heterogeneous but predictable patterns around the world that likely affect soil communities and overall ecosystem development.

Metabarcoding data reveal vertical multitaxa variation in topsoil communities during the colonization of deglaciated forelands

Ice-free areas are expanding worldwide due to dramatic glacier shrinkage and are undergoing rapid colonization by multiple lifeforms, thus representing key environments to study ecosystem development. It has been proposed that the colonization dynamics of deglaciated terrains is different between surface and deep soils but that the heterogeneity between communities inhabiting surface and deep soils decreases through time. Nevertheless, tests of this hypothesis remain scarce, and it is unclear whether patterns are consistent among different taxonomic groups. Here, we used environmental DNA metabarcoding to test whether community diversity and composition of six groups (Eukaryota, Bacteria, Mycota, Collembola, Insecta, and Oligochaeta) differ between the surface (0-5 cm) and deeper (7.5-20 cm) soil at different stages of development and across five Alpine glaciers. Taxonomic diversity increased with time since glacier retreat and with soil evolution. The pattern was consistent across groups and soil depths. For Eukaryota and Mycota, alpha-diversity was highest at the surface. Time since glacier retreat explained more variation of community composition than depth. Beta-diversity between surface and deep layers decreased with time since glacier retreat, supporting the hypothesis that the first 20 cm of soil tends to homogenize through time. Several molecular operational taxonomic units of bacteria and fungi were significant indicators of specific depths and/or soil development stages, confirming the strong functional variation of microbial communities through time and depth. The complexity of community patterns highlights the importance of integrating information from multiple taxonomic groups to unravel community variation in response to ongoing global changes.

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

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

Cascading effects of moth outbreaks on subarctic soil food webs

The increasing severity and frequency of natural disturbances requires a better understanding of their effects on all compartments of biodiversity. In Northern Fennoscandia, recent large-scale moth outbreaks have led to an abrupt change in plant communities from birch forests dominated by dwarf shrubs to grass-dominated systems. However, the indirect effects on the belowground compartment remained unclear. Here, we combined eDNA surveys of multiple trophic groups with network analyses to demonstrate that moth defoliation has far-reaching consequences on soil food webs. Following this disturbance, diversity and relative abundance of certain trophic groups declined (e.g., ectomycorrhizal fungi), while many others expanded (e.g., bacterivores and omnivores) making soil food webs more diverse and structurally different. Overall, the direct and indirect consequences of moth outbreaks increased belowground diversity at different trophic levels. Our results highlight that a holistic view of ecosystems improves our understanding of cascading effects of major disturbances on soil food webs.

Altitudinal Zonation of Green Algae Biodiversity in the French Alps

Mountain environments are marked by an altitudinal zonation of habitat types. They are home to a multitude of terrestrial green algae, who have to cope with abiotic conditions specific to high elevation, e.g., high UV irradiance, alternating desiccation, rain and snow precipitations, extreme diurnal variations in temperature and chronic scarceness of nutrients. Even though photosynthetic green algae are primary producers colonizing open areas and potential markers of climate change, their overall biodiversity in the Alps has been poorly studied so far, in particular in soil, where algae have been shown to be key components of microbial communities. Here, we investigated whether the spatial distribution of green algae followed the altitudinal zonation of the Alps, based on the assumption that algae settle in their preferred habitats under the pressure of parameters correlated with elevation. We did so by focusing on selected representative elevational gradients at distant locations in the French Alps, where soil samples were collected at different depths. Soil was considered as either a potential natural habitat or temporary reservoir of algae. We showed that algal DNA represented a relatively low proportion of the overall eukaryotic diversity as measured by a universal Eukaryote marker. We designed two novel green algae metabarcoding markers to amplify the Chlorophyta phylum and its Chlorophyceae class, respectively. Using our newly developed markers, we showed that elevation was a strong correlate of species and genus level distribution. Altitudinal zonation was thus determined for about fifty species, with proposed accessions in reference databases. In particular, Planophila laetevirens and Bracteococcus ruber related species as well as the snow alga Sanguina genus were only found in soil starting at 2,000 m above sea level. Analysis of environmental and bioclimatic factors highlighted the importance of pH and nitrogen/carbon ratios in the vertical distribution in soil. Capacity to grow heterotrophically may determine the Trebouxiophyceae over Chlorophyceae ratio. The intensity of freezing events (freezing degree days), proved also determinant in Chlorophyceae distribution. Guidelines are discussed for future, more robust and precise analyses of environmental algal DNA in mountain ecosystems and address green algae species distribution and dynamics in response to environmental changes.

Ancient DNA, lipid biomarkers and palaeoecological evidence reveals construction and life on early medieval lake settlements

Direct evidence of ancient human occupation is typically established through archaeological excavation. Excavations are costly and destructive, and practically impossible in some lake and wetland environments. We present here an alternative approach, providing direct evidence from lake sediments using DNA metabarcoding, steroid lipid biomarkers (bile acids) and from traditional environmental analyses. Applied to an early Medieval Celtic settlement in Ireland (a crannog) this approach provides a site chronology and direct evidence of human occupation, crops, animal farming and on-site slaughtering. This is the first independently-dated, continuous molecular archive of human activity from an archeological site, demonstrating a link between animal husbandry, food resources, island use. These sites are under threat but are impossible to preserve in-situ so this approach can be used, with or without excavation, to produce a robust and full site chronology and provide direct evidence of occupation, the use of plants and animals, and activities such as butchery.

Topsoil organic matter build-up in glacier forelands around the world

Since the last glacial maximum, soil formation related to ice-cover shrinkage has been one major sink of carbon accumulating as soil organic matter (SOM), a phenomenon accelerated by the ongoing global warming. In recently deglacierized forelands, processes of SOM accumulation, including those that control carbon and nitrogen sequestration rates and biogeochemical stability of newly sequestered carbon, remain poorly understood. Here, we investigate the build-up of SOM during the initial stages (up to 410 years) of topsoil development in 10 glacier forelands distributed on four continents. We test whether the net accumulation of SOM on glacier forelands (i) depends on the time since deglacierization and local climatic conditions (temperature and precipitation); (ii) is accompanied by a decrease in its stability and (iii) is mostly due to an increasing contribution of organic matter from plant origin. We measured total SOM concentration (carbon, nitrogen), its relative hydrogen/oxygen enrichment, stable isotopic (¹³ C, ¹⁵ N) and carbon functional groups (C-H, C=O, C=C) compositions, and its distribution in carbon pools of different thermal stability. We show that SOM content increases with time and is faster on forelands experiencing warmer climates. The build-up of SOM pools shows consistent trends across the studied soil chronosequences. During the first decades of soil development, the low amount of SOM is dominated by a thermally stable carbon pool with a small and highly thermolabile pool. The stability of SOM decreases with soil age at all sites, indicating that SOM storage is dominated by the accumulation of labile SOM during the first centuries of soil development, and suggesting plant carbon inputs to soil (SOM depleted in nitrogen, enriched in hydrogen and in aromatic carbon). Our findings highlight the potential vulnerability of SOM stocks from proglacial areas to decomposition and suggest that their durability largely depends on the relative contribution of carbon inputs from plants.

Effects of soil preservation for biodiversity monitoring using environmental DNA

Environmental DNA (eDNA) metabarcoding is becoming a key tool for biodiversity monitoring over large geographical or taxonomic scales and for elusive taxa such as soil organisms. Increasing sample sizes and interest in remote or extreme areas often require the preservation of soil samples and thus deviations from optimal standardized protocols. However, we still ignore the impact of different methods of soil sample preservation on the results of metabarcoding studies and there is no guideline for best practices so far. Here, we assessed the impact of four methods of soil sample preservation that can be conveniently used also in metabarcoding studies targeting remote or difficult to access areas. Tested methods include: preservation at room temperature for 6 hr, preservation at 4°C for 3 days, desiccation immediately after sampling and preservation for 21 days, and desiccation after 6 hr at room temperature and preservation for 21 days. For each preservation method, we benchmarked resulting estimates of taxon diversity and community composition of three different taxonomic groups (bacteria, fungi and eukaryotes) in three different habitats (forest, river bank and grassland) against results obtained under ideal conditions (i.e., extraction of eDNA immediately after sampling). Overall, the different preservation methods only marginally impaired results and only under certain conditions. When rare taxa were considered, we detected small but significant changes in molecular operational taxonomic units (MOTU) richness of bacteria, fungi and eukaryotes across treatments, but MOTU richness was similar across preservation methods if rare taxa were not considered. All the approaches were able to identify differences in community structure among habitats, and the communities retrieved using the different preservation conditions were extremely similar. We propose guidelines on the selection of the optimal soil sample preservation conditions for metabarcoding studies, depending on the practical constraints, costs and ultimate research goals.

Ecological specialization and niche overlap of subterranean rodents inferred from DNA metabarcoding diet analysis

Knowledge of how animal species use food resources available in the environment can increase our understanding of many ecological processes. However, obtaining this information using traditional methods is difficult for species feeding on a large variety of food items in highly diverse environments. We amplified the DNA of plants for 306 scat and 40 soil samples, and applied an environmental DNA metabarcoding approach to investigate food preferences, degree of diet specialization and diet overlap of seven herbivore rodent species of the genus Ctenomys distributed in southern and midwestern Brazil. The metabarcoding approach revealed that these species consume more than 60% of the plant families recovered in soil samples, indicating generalist feeding habits of ctenomyids. The family Poaceae was the most common food resource retrieved in scats of all species as well in soil samples. Niche overlap analysis indicated high overlap in the plant families and molecular operational taxonomic units consumed, mainly among the southern species. Interspecific differences in diet composition were influenced, among other factors, by the availability of resources in the environment. In addition, our results provide support for the hypothesis that the allopatric distributions of ctenomyids allow them to exploit the same range of resources when available, possibly because of the absence of interspecific competition.

Persistence of arctic-alpine flora during 24,000 years of environmental change in the Polar Urals

Plants adapted to extreme conditions can be at high risk from climate change; arctic-alpine plants, in particular, could "run out of space" as they are out-competed by expansion of woody vegetation. Mountain regions could potentially provide safe sites for arctic-alpine plants in a warmer climate, but empirical evidence is fragmentary. Here we present a 24,000-year record of species persistence based on sedimentary ancient DNA (sedaDNA) from Lake Bolshoye Shchuchye (Polar Urals). We provide robust evidence of long-term persistence of arctic-alpine plants through large-magnitude climate changes but document a decline in their diversity during a past expansion of woody vegetation. Nevertheless, most of the plants that were present during the last glacial interval, including all of the arctic-alpines, are still found in the region today. This underlines the conservation significance of mountain landscapes via their provision of a range of habitats that confer resilience to climate change, particularly for arctic-alpine taxa.

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

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

Environmental DNA provides information on sediment sources: A study in catchments affected by Fukushima radioactive fallout

An excessive supply of sediment is observed in numerous rivers across the world where it leads to deleterious impacts. Information on the sources delivering this material to waterbodies is required to design effective management measures, and sediment tracing or fingerprinting techniques are increasingly used to quantify the amount of sediment derived from different sources. However, the current methods used to identify the land use contributions to sediment have a limited discrimination power. Here, we investigated the potential of environmental DNA (eDNA) to provide more detailed information on the plant species found in sediment source areas as a next generation fingerprint. To this end, flood sediment deposits (n = 12) were collected in 2017 in two catchments impacted by the Fukushima radioactive fallout along differing river sections draining forests, cropland or a mix of both land uses. Conventional fingerprints (i.e. fallout radionuclides and organic matter properties) were also measured in these samples. The conventional fingerprint model results showed that most sediment samples contained a dominant proportion of subsoil material. Nevertheless, the eDNA information effectively discriminated the three above-mentioned groups of sediment, with the dominance of tree, shrub and fern species in sediment sampled in rivers draining forests versus a majority of grass, algae and cultivated plant species in sediment collected in rivers draining cropland. Based on these encouraging results, future research should examine the potential of eDNA in mixed land use catchments where the contribution of topsoil to sediment dominates and where the cultivation of land has not been abandoned in order to better characterize the memory effect of eDNA in soils and sediment.

DNA from lake sediments reveals long-term ecosystem changes after a biological invasion

What are the long-term consequences of invasive species? After invasion, how long do ecosystems require to reach a new equilibrium? Answering these questions requires long-term, high-resolution data that are vanishingly rare. We combined the analysis of environmental DNA extracted from a lake sediment core, coprophilous fungi, and sedimentological analyses to reconstruct 600 years of ecosystem dynamics on a sub-Antarctic island and to identify the impact of invasive rabbits. Plant communities remained stable from AD 1400 until the 1940s, when the DNA of invasive rabbits was detected in sediments. Rabbit detection corresponded to abrupt changes of plant communities, with a continuous decline of a dominant plant species. Furthermore, erosion rate abruptly increased with rabbit abundance. Rabbit impacts were very fast and were stronger than the effects of climate change during the 20th century. Lake sediments can allow an integrated temporal analysis of ecosystems, revealing the impact of invasive species over time and improving our understanding of underlying mechanisms.

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

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

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

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

Understanding the evolution of holoparasitic plants: the complete plastid genome of the holoparasite Cytinus hypocistis (Cytinaceae)

Background and Aims Plant plastid genomes are highly conserved in size, gene content and structure; however, parasitic plants are a noticeable exception to this evolutionary stability. Although the evolution of parasites could help to better understand plastome evolution in general, complete plastomes of parasites have been sequenced only for some lineages so far. Here we contribute to filling this gap by providing and analysing the complete plastome sequence of Cytinus hypocistis, the first parasite sequenced for Malvales and a species suspected to have an extremely small genome. Methods We sequenced and assembled de novo the plastid genome of Cytinus hypocistis using a shotgun approach on genomic DNA. Phylogenomic analyses based on coding regions were performed on Malvidae. For each coding region present in Cytinus, we tested for relaxation or intensification of selective pressures in the Cytinus lineage compared with autotrophic Malvales. Key Results Cytinus hypocistis has an extremely divergent genome that is among the smallest sequenced to date (19·4 kb), with only 23 genes and no inverted repeat regions. Phylogenomic analysis confirmed the position of Cytinus within Malvales. All coding regions of Cytinus plastome presented very high substitution rates compared with non-parasitic Malvales. Conclusions Some regions were inferred to be under relaxed negative selection in Cytinus, suggesting that further plastome reduction is occurring due to relaxed purifying selection associated with the loss of photosynthetic activity. On the other hand, increased selection intensity and strong positive selection were detected for rpl22 in the Cytinus lineage, which might indicate an evolutionary role in the host-parasite arms race, a point that needs further research.

Long-lasting modification of soil fungal diversity associated with the introduction of rabbits to a remote sub-Antarctic archipelago

During the late nineteenth century, Europeans introduced rabbits to many of the sub-Antarctic islands, environments that prior to this had been devoid of mammalian herbivores. The impacts of rabbits on indigenous ecosystems are well studied; notably, they cause dramatic changes in plant communities and promote soil erosion. However, the responses of fungal communities to such biotic disturbances remain unexplored. We used metabarcoding of soil extracellular DNA to assess the diversity of plant and fungal communities at sites on the sub-Antarctic Kerguelen Islands with contrasting histories of disturbance by rabbits. Our results suggest that on these islands, the simplification of plant communities and increased erosion resulting from the introduction of rabbits have driven compositional changes, including diversity reductions, in indigenous soil fungal communities. Moreover, there is no indication of recovery at sites from which rabbits were removed 20 years ago. These results imply that introduced herbivores have long-lasting and multifaceted effects on fungal biodiversity as well as highlight the low resiliency of sub-Antarctic ecosystems.

DNA metabarcoding diet analysis for species with parapatric vs sympatric distribution: a case study on subterranean rodents

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.

Reconstructing long-term human impacts on plant communities: an ecological approach based on lake sediment DNA

Paleoenvironmental studies are essential to understand biodiversity changes over long timescales and to assess the relative importance of anthropogenic and environmental factors. Sedimentary ancient DNA (sedaDNA) is an emerging tool in the field of paleoecology and has proven to be a complementary approach to the use of pollen and macroremains for investigating past community changes. SedaDNA-based reconstructions of ancient environments often rely on indicator taxa or expert knowledge, but quantitative ecological analyses might provide more objective information. Here, we analysed sedaDNA to investigate plant community trajectories in the catchment of a high-elevation lake in the Alps over the last 6400 years. We combined data on past and present plant species assemblages along with sedimentological and geochemical records to assess the relative impact of human activities through pastoralism, and abiotic factors (temperature and soil evolution). Over the last 6400 years, we identified significant variation in plant communities, mostly related to soil evolution and pastoral activities. An abrupt vegetational change corresponding to the establishment of an agropastoral landscape was detected during the Late Holocene, approximately 4500 years ago, with the replacement of mountain forests and tall-herb communities by heathlands and grazed lands. Our results highlight the importance of anthropogenic activities in mountain areas for the long-term evolution of local plant assemblages. SedaDNA data, associated with other paleoenvironmental proxies and present plant assemblages, appear to be a relevant tool for reconstruction of plant cover history. Their integration, in conjunction with classical tools, offers interesting perspectives for a better understanding of long-term ecosystem dynamics under the influence of human-induced and environmental drivers.

Highly overlapping winter diet in two sympatric lemming species revealed by DNA metabarcoding

Sympatric species are expected to minimize competition by partitioning resources, especially when these are limited. Herbivores inhabiting the High Arctic in winter are a prime example of a situation where food availability is anticipated to be low, and thus reduced diet overlap is expected. We present here the first assessment of diet overlap of high arctic lemmings during winter based on DNA metabarcoding of feces. In contrast to previous analyses based on microhistology, we found that the diets of both collared (Dicrostonyx groenlandicus) and brown lemmings (Lemmus trimucronatus) on Bylot Island were dominated by Salix while mosses, which were significantly consumed only by the brown lemming, were a relatively minor food item. The most abundant plant taxon, Cassiope tetragona, which alone composes more than 50% of the available plant biomass, was not detected in feces and can thus be considered to be non-food. Most plant taxa that were identified as food items were consumed in proportion to their availability and none were clearly selected for. The resulting high diet overlap, together with a lack of habitat segregation, indicates a high potential for resource competition between the two lemming species. However, Salix is abundant in the winter habitats of lemmings on Bylot Island and the non-Salix portion of the diets differed between the two species. Also, lemming grazing impact on vegetation during winter in the study area is negligible. Hence, it seems likely that the high potential for resource competition predicted between these two species did not translate into actual competition. This illustrates that even in environments with low primary productivity food resources do not necessarily generate strong competition among herbivores.

Replication levels, false presences and the estimation of the presence/absence from eDNA metabarcoding data

Environmental DNA (eDNA) metabarcoding is increasingly used to study the present and past biodiversity. eDNA analyses often rely on amplification of very small quantities or degraded DNA. To avoid missing detection of taxa that are actually present (false negatives), multiple extractions and amplifications of the same samples are often performed. However, the level of replication needed for reliable estimates of the presence/absence patterns remains an unaddressed topic. Furthermore, degraded DNA and PCR/sequencing errors might produce false positives. We used simulations and empirical data to evaluate the level of replication required for accurate detection of targeted taxa in different contexts and to assess the performance of methods used to reduce the risk of false detections. Furthermore, we evaluated whether statistical approaches developed to estimate occupancy in the presence of observational errors can successfully estimate true prevalence, detection probability and false-positive rates. Replications reduced the rate of false negatives; the optimal level of replication was strongly dependent on the detection probability of taxa. Occupancy models successfully estimated true prevalence, detection probability and false-positive rates, but their performance increased with the number of replicates. At least eight PCR replicates should be performed if detection probability is not high, such as in ancient DNA studies. Multiple DNA extractions from the same sample yielded consistent results; in some cases, collecting multiple samples from the same locality allowed detecting more species. The optimal level of replication for accurate species detection strongly varies among studies and could be explicitly estimated to improve the reliability of results.

Sources of variation in small rodent trophic niche: new insights from DNA metabarcoding and stable isotope analysis

Intraspecific competition for food is expected to increase the trophic niche width of consumers, defined here as their diet diversity, but this process has been little studied in herbivores. Population densities of small rodents fluctuate greatly, providing a good study model to evaluate effects of competition on trophic niche. We studied resource use in five arctic small rodent populations of four species combining DNA metabarcoding of stomach contents and stable isotope analysis (SIA). Our results suggest that for small rodents, the most pronounced effect of competition on trophic niche is due to increased use of secondary habitats and to habitat-specific diets, rather than an expansion of trophic niche in primary habitat. DNA metabarcoding and SIA provided complementary information about the composition and temporal variation of herbivore diets. Combing these two approaches requires caution, as the underlying processes causing observed patterns may differ between methodologies due to different spatiotemporal scales.

Arctic Small Rodents Have Diverse Diets and Flexible Food Selection

The ecology of small rodent food selection is poorly understood, as mammalian herbivore food selection theory has mainly been developed by studying ungulates. Especially, the effect of food availability on food selection in natural habitats where a range of food items are available is unknown. We studied diets and selectivity of grey-sided voles (Myodes rufocanus) and tundra voles (Microtus oeconomus), key herbivores in European tundra ecosystems, using DNA metabarcoding, a novel method enabling taxonomically detailed diet studies. In order to cover the range of food availabilities present in the wild, we employed a large-scale study design for sampling data on food availability and vole diets. Both vole species had ingested a range of plant species and selected particularly forbs and grasses. Grey-sided voles also selected ericoid shrubs and tundra voles willows. Availability of a food item rarely affected its utilization directly, although seasonal changes of diets and selection suggest that these are positively correlated with availability. Moreover, diets and selectivity were affected by availability of alternative food items. These results show that the focal sub-arctic voles have diverse diets and flexible food preferences and rarely compensate low availability of a food item with increased searching effort. Diet diversity itself is likely to be an important trait and has previously been underrated owing to methodological constraints. We suggest that the roles of alternative food item availability and search time limitations for small rodent feeding ecology should be investigated.

Nomenclature

Annotated Checklist of the Panarctic Flora (PAF), Vascular plants. Available at: http://nhm2.uio.no/paf/, accessed 15.6.2012.

Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.

DNA from soil mirrors plant taxonomic and growth form diversity

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 <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 taxonomic and growth form diversity estimated from conventional above-ground surveys. Second, we assessed DNA persistence using samples from formerly cultivated soils in temperate environments. We found that the 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 were efficiently recovered. The results open unprecedented opportunities for large-scale DNA-based biodiversity studies across a range of taxonomic groups using standardized metabarcoding approaches.

AFLP markers reveal high clonal diversity and extreme longevity in four key arctic-alpine species

We investigated clonal diversity, genet size structure and genet longevity in populations of four arctic-alpine plants (Carex curvula, Dryas octopetala, Salix herbacea and Vaccinium uliginosum) to evaluate their persistence under past climatic oscillations and their potential resistance to future climate change. The size and number of genets were determined by an analysis of amplified fragment length polymorphisms and a standardized sampling design in several European arctic-alpine populations, where these species are dominant in the vegetation. Genet age was estimated by dividing the size by the annual horizontal size increment from in situ growth measurements. Clonal diversity was generally high but differed among species, and the frequency distribution of genet size was strongly left-skewed. The largest C. curvula genet had an estimated minimum age of c. 4100 years and a maximum age of c. 5000 years, although 84.8% of the genets in this species were <200 years old. The oldest genets of D. octopetala, S. herbacea and V. uliginosum were found to be at least 500, 450 and 1400 years old, respectively. These results indicate that individuals in the studied populations have survived pronounced climatic oscillations, including the Little Ice Age and the postindustrial warming. The presence of genets in all size classes and the dominance of presumably young individuals suggest repeated recruitment over time, a precondition for adaptation to changing environmental conditions. Together, persistence and continuous genet turnover may ensure maximum ecosystem resilience. Thus, our results indicate that long-lived clonal plants in arctic-alpine ecosystems can persist, despite considerable climatic change.

Tracking genes of ecological relevance using a genome scan in two independent regional population samples of Arabis alpina

Understanding the genetic basis of adaptation in response to environmental variation is fundamental as adaptation plays a key role in the extension of ecological niches to marginal habitats and in ecological speciation. Based on the assumption that some genomic markers are correlated to environmental variables, we aimed to detect loci of ecological relevance in the alpine plant Arabis alpina L. sampled in two regions, the French (99 locations) and the Swiss (109 locations) Alps. We used an unusually large genome scan [825 amplified fragment length polymorphism loci (AFLPs)] and four environmental variables related to temperature, precipitation and topography. We detected linkage disequilibrium among only 3.5% of the considered AFLP loci. A population structure analysis identified no admixture in the study regions, and the French and Swiss Alps were differentiated and therefore could be considered as two independent regions. We applied generalized estimating equations (GEE) to detect ecologically relevant loci separately in the French and Swiss Alps. We identified 78 loci of ecological relevance (9%), which were mainly related to mean annual minimum temperature. Only four of these loci were common across the French and Swiss Alps. Finally, we discuss that the genomic characterization of these ecologically relevant loci, as identified in this study, opens up new perspectives for studying functional ecology in A. alpina, its relatives and other alpine plant species.

Selection criteria for scoring amplified fragment length polymorphisms (AFLPs) positively affect the reliability of population genetic parameter estimates

A reliable data set is a fundamental prerequisite for consistent results and conclusions in population genetic studies. However, marker scoring of genetic fingerprints such as amplified fragment length polymorphisms (AFLPs) is a highly subjective procedure, inducing inconsistencies owing to personal or laboratory-specific criteria. We applied two alternative marker selection algorithms, the newly developed script scanAFLP and the recently published AFLPScore, to a large AFLP genome scan to test how population genetic parameters and error rates were affected. These results were confronted with replicated random selections of marker subsets. We show that the newly developed marker selection criteria reduced the mismatch error rate and had a notable influence on estimates of genetic diversity and differentiation. Both effects are likely to influence biological inference. For example, genetic diversity (HS) was 29% lower while genetic differentiation (FST) was 8% higher when applying scanAFLP compared with AFLPScore. Likewise, random selections of markers resulted in substantial deviations of population genetic parameters compared with the data sets including specific selection criteria. These randomly selected marker sets showed surprisingly low variance among replicates. We conclude that stringent marker selection and phenotype calling reduces noise in the data set while retaining patterns of population genetic structure.

Analysing diet of small herbivores: the efficiency of DNA barcoding coupled with high-throughput pyrosequencing for deciphering the composition of complex plant mixtures

BACKGROUND

In order to understand the role of herbivores in trophic webs, it is essential to know what they feed on. Diet analysis is, however, a challenge in many small herbivores with a secretive life style. In this paper, we compare novel (high-throughput pyrosequencing) DNA barcoding technology for plant mixture with traditional microhistological method. We analysed stomach contents of two ecologically important subarctic vole species, Microtus oeconomus and Myodes rufocanus, with the two methods. DNA barcoding was conducted using the P6-loop of the chloroplast trnL (UAA) intron.

RESULTS

Although the identified plant taxa in the diets matched relatively well between the two methods, DNA barcoding gave by far taxonomically more detailed results. Quantitative comparison of results was difficult, mainly due to low taxonomic resolution of the microhistological method, which also in part explained discrepancies between the methods. Other discrepancies were likely due to biases mostly in the microhistological analysis.

CONCLUSION

We conclude that DNA barcoding opens up for new possibilities in the study of plant-herbivore interactions, giving a detailed and relatively unbiased picture of food utilization of herbivores.

Exonuclease activity of proofreading DNA polymerases is at the origin of artifacts in molecular profiling studies

CE fingerprint methods are commonly used in microbial ecology. We have previously noticed that the position and number of peaks in CE-SSCP (single-strand conformation polymorphism) profiles depend on the DNA polymerase used in PCR [1]. Here, we studied the fragments produced by Taq polymerase as well as four commercially available proofreading polymerases, using the V3 region of the Escherichia coli rss gene as a marker. PCR products rendered multiple peaks in denaturing CE; Taq polymerase was observed to produce the longest fragments. Incubation of the fragments with T4 DNA polymerase indicated that the 3'-ends of the proofreading polymerase amplicons were recessed, while the Taq amplicon was partially +A tailed. Treatment of the PCR product with proofreading DNA polymerase rendered trimmed fragments. This was due to the 3'-5' exonuclease activity of these enzymes, which is essential for proofreading. The nuclease activity was reduced by increasing the concentration of dNTP. The Platinum Pfx DNA polymerase generated very few artifacts and could produce 85% of blunted PCR products. Nevertheless, despite the higher error rate, we recommend the use of Taq polymerase rather than proofreading in the framework for molecular fingerprint studies. They are more cost-effective and therefore ideally suited for high-throughput analysis; the +A tail artifact rate can be controlled by modifying the PCR primers and the reaction conditions.

CE-SSCP and CE-FLA, simple and high-throughput alternatives for fungal diversity studies

Fungal communities are key components of soil, but the study of their ecological significance is limited by a lack of appropriated methods. For instance, the assessment of fungi occurrence and spatio-temporal variation in soil requires the analysis of a large number of samples. The molecular signature methods provide a useful tool to monitor these microbial communities and can be easily adapted to capillary electrophoresis (CE) allowing high-throughput studies. Here we assess the suitability of CE-FLA (Fragment Length Polymorphism, denaturing conditions) and CE-SSCP (Single-Stranded Conformation Polymorphism, native conditions) applied to environmental studies since they require a short molecular marker and no post-PCR treatments. We amplified the ITS1 region from 22 fungal strains isolated from an alpine ecosystem and from total genomic DNA of alpine and infiltration basin soils. The CE-FLA and CE-SSCP separated 17 and 15 peaks respectively from a mixture of 19 strains. For the alpine soil-metagenomic DNA, the FLA displayed more peaks than the SSCP and the converse result was found for infiltration basin sediments. We concluded that CE-FLA and CE-SSCP of ITS1 region provided complementary information. In order to improve CE-SSCP sensitivity, we tested its resolution according to migration temperature and found 32 degrees C to be optimal. Because of their simplicity, quickness and reproducibility, we found that these two methods were promising for high-throughput studies of soil fungal communities.

Improvements of polymerase chain reaction and capillary electrophoresis single-strand conformation polymorphism methods in microbial ecology: toward a high-throughput method for microbial diversity studies in soil

The molecular signature of bacteria from soil ecosystems is an important tool for studying microbial ecology and biogeography. However, a high-throughput technology is needed for such studies. In this article, we tested the suitability of available methods ranging from soil DNA extraction to capillary electrophoresis single-strand conformation polymorphism (CE-SSCP) for high-throughput studies. Our results showed that the extraction method does not dramatically influence CE-SSCP profiles, and that DNA extraction of a 0.25 g soil sample is sufficient to observe overall bacterial diversity in soil matrices. The V3 region of the 16S rRNA gene was amplified by PCR, and the extension time was found to be critical. We have also found that proofreading DNA polymerases generate a better signal in CE-SSCP profiles. Experiments performed with different soil matrices revealed the repeatability, efficiency, and consistency of CE-SSCP. Studies on PCR and CE-SSCP using single-species genomic DNA as a matrix showed that several ribotypes may migrate at the same position, and also that single species can produce double peaks. Thus, the extrapolation between number of peaks and number of species remains difficult. Additionally, peak detection is limited by the analysis software. We conclude that the presented method, including CE-SSCP and the analyzing step, is a simple and effective technique to obtain the molecular signature of a given soil sample.

Power and limitations of the chloroplast trnL (UAA) intron for plant DNA barcoding

DNA barcoding should provide rapid, accurate and automatable species identifications by using a standardized DNA region as a tag. Based on sequences available in GenBank and sequences produced for this study, we evaluated the resolution power of the whole chloroplast trnL (UAA) intron (254–767 bp) and of a shorter fragment of this intron (the P6 loop, 10–143 bp) amplified with highly conserved primers. The main limitation of the whole trnL intron for DNA barcoding remains its relatively low resolution (67.3% of the species from GenBank unambiguously identified). The resolution of the P6 loop is lower (19.5% identified) but remains higher than those of existing alternative systems. The resolution is much higher in specific contexts such as species originating from a single ecosystem, or commonly eaten plants. Despite the relatively low resolution, the whole trnL intron and its P6 loop have many advantages: the primers are highly conserved, and the amplification system is very robust. The P6 loop can even be amplified when using highly degraded DNA from processed food or from permafrost samples, and has the potential to be extensively used in food industry, in forensic science, in diet analyses based on feces and in ancient DNA studies.

Authenticated DNA from ancient wood remains

BACKGROUND

The reconstruction of biological processes and human activities during the last glacial cycle relies mainly on data from biological remains. Highly abundant tissues, such as wood, are candidates for a genetic analysis of past populations. While well-authenticated DNA has now been recovered from various fossil remains, the final 'proof' is still missing for wood, despite some promising studies.

SCOPE

The goal of this study was to determine if ancient wood can be analysed routinely in studies of archaeology and palaeogenetics. An experiment was designed which included blind testing, independent replicates, extensive contamination controls and rigorous statistical tests. Ten samples of ancient wood from major European forest tree genera were analysed with plastid DNA markers.

CONCLUSIONS

Authentic DNA was retrieved from wood samples up to 1,000 years of age. A new tool for real-time vegetation history and archaeology is ready to use.

Impact of ice ages on circumpolar molecular diversity: insights from an ecological key species

We address the impact of the ice age cycles on intraspecific cpDNA diversity, for the first time on the full circumboreal-circumarctic scale. The bird-dispersed bog bilberry (or arctic blueberry, Vaccinium uliginosum) is a key component of northern ecosystems and is here used to assess diversity in previously glaciated vs. unglaciated areas and the importance of Beringia as a refugium and source for interglacial expansion. Eighteen chloroplast DNA haplotypes were observed in and among 122 populations, grouping into three main lineages which probably diverged before, and thus were affected more or less independently by, all major glaciations. The boreal 'Amphi-Atlantic lineage' included one haplotype occurring throughout northern Europe and one occurring in eastern North America, suggesting expansion from at least two bottlenecked, glacial refugium populations. The boreal 'Beringian lineage' included seven haplotypes restricted to Beringia and the Pacific coast of USA. The 'Arctic-Alpine lineage' included nine haplotypes, one of them fully circumpolar. This lineage was unexpectedly diverse, also in previously glaciated areas, suggesting that it thrived on the vast tundras during the ice ages and recolonized deglaciated terrain over long distances. Its largest area of persistence during glaciations was probably situated in the north, stretching from Beringia and far into Eurasia, and it probably also survived the last glaciation in southern mountain ranges. Although Beringia apparently was important for the initial divergence and expansion of V. uliginosum as well as for continuous survival of both the Beringian and Arctic-Alpine lineages during all ice ages, this region played a minor role as a source for later interglacial expansions.

Genetic structure of Rhododendron ferrugineum at a wide range of spatial scales

Rhododendron ferrugineum L. (Ericaceae) is a subalpine shrub found throughout the Pyrenees and Alps at elevations of 1600-2200 m. We examined relationships between genetic and geographic distance, using 115 dominant amplified fragment length polymorphism (AFLP) markers to assess genetic structure over a wide range of spatial scales. We sampled 17 sites with distances of 4 km to more than 1000 km between them. At these scales we detected no association between geographic distance and genetic distance between populations. This suggests that genetic drift and gene flow are not in equilibrium for these populations. This pattern could have resulted from recent and rapid postglacial colonization, from more recent human disturbance, or as a function of frequent and random "natural" long-distance colonization. At two of our sites we used transects (two horizontal and two vertical with respect to slope at each site) to sample at distances ranging from 10 m to more than 5000 m. At this scale we observed a positive relationship between genetic and spatial distance along two vertical transects, one at each site. We hypothesize that isolation-by-distance at this smaller scale is a function of restricted gene flow via seed dispersal.