The development of terrestrial ecosystems emerging after glacier retreat.

The global retreat of glaciers is dramatically altering mountain and high-latitude landscapes, with new ecosystems developing from apparently barren substrates1-4. The study of these emerging ecosystems is critical to understanding how climate change interacts with microhabitat and biotic communities and determines the future of ice-free terrains1,5. Here, using a comprehensive characterization of ecosystems (soil properties, microclimate, productivity and biodiversity by environmental DNA metabarcoding6) across 46 proglacial landscapes worldwide, we found that all the environmental properties change with time since glaciers retreated, and that temperature modulates the accumulation of soil nutrients. The richness of bacteria, fungi, plants and animals increases with time since deglaciation, but their temporal patterns differ. Microorganisms colonized most rapidly in the first decades after glacier retreat, whereas most macroorganisms took longer. Increased habitat suitability, growing complexity of biotic interactions and temporal colonization all contribute to the increase in biodiversity over time. These processes also modify community composition for all the groups of organisms. Plant communities show positive links with all other biodiversity components and have a key role in ecosystem development. These unifying patterns provide new insights into the early dynamics of deglaciated terrains and highlight the need for integrated surveillance of their multiple environmental properties5.

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.

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.

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.

Vegetation structure and climate shape mountain arthropod distributions across trophic levels.

Arthropods play a vital role in ecosystems; yet, their distributions remain poorly understood, particularly in mountainous regions. This study delves into the modelling of the distribution of 31 foliar arthropod genera in the French Alps, using a comprehensive approach encompassing multi-trophic sampling, community DNA metabarcoding and random forest models. The results underscore the significant importance of vegetation structure, such as herbaceous vegetation density, and forest density and heterogeneity, along with climate, in shaping the distributions of most arthropods. These responses to environmental gradients are consistent across trophic groups, with the exception of nectarivores, whose distributions are more sensitive to landscape structure and water availability. By leveraging community DNA metabarcoding, this study sheds light on the understudied drivers of arthropod distributions, emphasizing the importance of modelling across diverse trophic groups to anticipate arthropod responses to global change.

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.

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.

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 (13 C, 15 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.

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.

Fifty thousand years of Arctic vegetation and megafaunal diet.

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 also 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.

Body size determines soil community assembly in a tropical forest.

Tropical forests shelter an unparalleled biological diversity. The relative influence of environmental selection (i.e., abiotic conditions, biotic interactions) and stochastic-distance-dependent neutral processes (i.e., demography, dispersal) in shaping communities has been extensively studied for various organisms, but has rarely been explored across a large range of body sizes, in particular in soil environments. We built a detailed census of the whole soil biota in a 12-ha tropical forest plot using soil DNA metabarcoding. We show that the distribution of 19 taxonomic groups (ranging from microbes to mesofauna) is primarily stochastic, suggesting that neutral processes are prominent drivers of the assembly of these communities at this scale. We also identify aluminium, topography and plant species identity as weak, yet significant drivers of soil richness and community composition of bacteria, protists and to a lesser extent fungi. Finally, we show that body size, which determines the scale at which an organism perceives its environment, predicted the community assembly across taxonomic groups, with soil mesofauna assemblages being more stochastic than microbial ones. These results suggest that the relative contribution of neutral processes and environmental selection to community assembly directly depends on body size. Body size is hence an important determinant of community assembly rules at the scale of the ecological community in tropical soils and should be accounted for in spatial models of tropical soil food webs.

Mapping the imprint of biotic interactions on ß-diversity.

Investigating how trophic interactions influence the ß-diversity of meta-communities is of paramount importance to understanding the processes shaping biodiversity distribution. Here, we apply a statistical method for inferring the strength of spatial dependencies between pairs of species groups. Using simulated community data generated from a multi-trophic model, we showed that this method can approximate biotic interactions in multi-trophic communities based on ß-diversity patterns across groups. When applied to soil multi-trophic communities along an elevational gradient in the French Alps, we found that fungi make a major contribution to the structuring of ß-diversity across trophic groups. We also demonstrated that there were strong spatial dependencies between groups known to interact specifically (e.g. plant-symbiotic fungi, bacteria-nematodes) and that the influence of environment was less important than previously reported in the literature. Our method paves the way for a better understanding and mapping of multi-trophic communities through space and time.

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 20th 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.

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.

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.

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.

Genetic diversity in widespread species is not congruent with species richness in alpine plant communities.

The Convention on Biological Diversity (CBD) aims at the conservation of all three levels of biodiversity, that is, ecosystems, species and genes. Genetic diversity represents evolutionary potential and is important for ecosystem functioning. Unfortunately, genetic diversity in natural populations is hardly considered in conservation strategies because it is difficult to measure and has been hypothesised to co-vary with species richness. This means that species richness is taken as a surrogate of genetic diversity in conservation planning, though their relationship has not been properly evaluated. We tested whether the genetic and species levels of biodiversity co-vary, using a large-scale and multi-species approach. We chose the high-mountain flora of the Alps and the Carpathians as study systems and demonstrate that species richness and genetic diversity are not correlated. Species richness thus cannot act as a surrogate for genetic diversity. Our results have important consequences for implementing the CBD when designing conservation strategies.

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.

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.