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.

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.

Heterogeneous changes of soil microclimate in high mountains and glacier forelands.

Landscapes nearby glaciers are disproportionally affected by climate change, but we lack detailed information on microclimate variations that can modulate the impacts of global warming on proglacial ecosystems and their biodiversity. Here, we use near-subsurface soil temperatures in 175 stations from polar, equatorial and alpine glacier forelands to generate high-resolution temperature reconstructions, assess spatial variability in microclimate change from 2001 to 2020, and estimate whether microclimate heterogeneity might buffer the severity of warming trends. Temporal changes in microclimate are tightly linked to broad-scale conditions, but the rate of local warming shows great spatial heterogeneity, with faster warming nearby glaciers and during the warm season, and an extension of the snow-free season. Still, most of the fine-scale spatial variability of microclimate is one-to-ten times larger than the temporal change experienced during the past 20 years, indicating the potential for microclimate to buffer climate change, possibly allowing organisms to withstand, at least temporarily, the effects of warming.

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.

Nitrogen and carbon stable isotope analysis sheds light on trophic competition between two syntopic land iguana species from Galápagos.

Coexistence between closely related species can lead to intense competition for resources. Stable isotope analysis (SIA) is a reliable tool to estimate the extent of species competition. We employed SIA to evaluate niche partitioning among two syntopic species of Galápagos land iguanas: Conolophus subcristatus and C. marthae. Samples were collected on Wolf Volcano, Isabela Island, where C. marthae is endemic and syntopic with C. subcristatus. We determined δ13C and δ15N ratios and described the isotopic niche of each species using corrected standard ellipse area (SEAc). We tested for differentiation between the isotopic niches, while controlling for sex, body size, spatial location of samples and mean annual primary productivity at capture points, using bivariate linear models. Despite the extensive overlap of the isotopic niches, we found species and sex to be a significant, interacting predictor of a sample's location in the δ13C, δ15N space, indicating the existence of niche partitioning mechanisms acting between species and sexes. We also found that body size and productivity at the capture points, compounded with yet undetermined spatial effects, explain ca. 75% of the differences observed between species and sexes, providing evidence for differential microhabitat and food-items usage. Our study provides essential baselines for evaluating conservation actions for C. marthae, such as the potential translocation to a sanctuary area free of competition from C. subcristatus.

Condition- and context-dependent variation of sexual dimorphism across lizard populations at different spatial scales.

The evolution of sexual dimorphism (SD) is driven by intricate interplays between sexual and natural selection. When it comes to SD variation within populations, however, environmental factors play a major role. Sexually selected traits are expected to be strongly dependent on individual body condition, which is influenced by the local environment that individuals experience. As a consequence, the degree of SD may also depend on resource availability. Here, we investigated the potential drivers of SD expression at two sexually dimorphic morphometric traits, body size (snout vent length) and head shape (head geometric morphometrics), in the Italian wall lizard (Podarcis siculus). We assessed the existence of condition- and context-dependent SD across ten islands of the Aeolian archipelago (southern Italy), at within- and among-population scales. We observed strong geographical variation of SD among islands, and tested three potential SD predictors related to resource availability (individual body condition, ecosystem productivity, temperature). Body condition and ecosystem productivity were the main drivers of body size SD variation, and body condition was also the main driver for head shape SD. Our results highlight that the expression of SD in the Italian wall lizard is both condition- and context-dependent. These results are congruent at within- and among-populations scales highlighting that spatial multi-scale analysis represents a useful approach to understand patterns of SD expression.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout.

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen-mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC-DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human-mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

Climate and land-use changes drive biodiversity turnover in arthropod assemblages over 150 years.

Long-term studies are essential to understand the impacts of global changes on the multiple facets of biological diversity. Here, we use distribution data for over 600 species of arthropods collected over 150 years from locations across Italy and test how multiple environmental stressors (climate, land use and human population density) influenced assemblage composition and functionality. By carefully reconstructing the temporal changes in these stressors, we explicitly tested how environmental changes can determine the observed changes in taxonomic and functional diversity. We found that rapid changes in precipitation destabilize the assemblages and maximize colonization and extinction rates, especially when coupled with changes in human population density (for taxonomy) or temperature (for functionality). Higher microclimatic heterogeneity increases the stability of biodiversity by reducing taxonomic and functional loss. Finally, changes in natural habitats increased colonization, influencing taxonomic nestedness and functional replacement. The integration of long-term datasets combining distributions, climate and traits may deepen our understanding of the processes underlying biodiversity responses to global-scale drivers.

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.

ClimCKmap, a spatially, temporally and climatically explicit distribution database for the Italian fauna.

Understanding and counteracting biodiversity losses requires quantitative knowledge on species distribution and abundance across space and time, as well as integrated and interoperable information on climate conditions and climatic changes. In this paper we developed a new biodiversity-climate database for Italy, ClimCKmap, based on the critical analysis, quality estimation and subsequent integration of the CKmap database with several high-resolution climate datasets. The original database was quality-checked for errors in toponym, species name and dating; the retained records were georeferenced and their distribution polygonised via Voronoi tessellation. We then integrated the species distribution information with several high-resolution climatic datasets: average monthly minimum and maximum temperature and total monthly precipitation were reconstructed for each Voronoi cell and year. The resulting database contains 268,977 occurrence records from 8,445 binomials and 16,332 localities, dating between 1680 and 2006 CE. This dataset, fully available at https://doi.org/10.6084/m9.figshare.7906739.v4 and http://hdl.handle.net/21.11125/a91f85cb-befd-4e14-8e83-24f17c4a0491 , represents the largest, fully quality-checked, spatially, temporally and climatically explicit distribution database ever assembled for the Italian fauna, now ready for scientific exploitation.

MHC structuring and divergent allele advantage in a urodele amphibian: a hierarchical multi-scale approach.

Proteins encoded by extraordinarily polymorphic major histocompatibility complex (MHC) genes are involved in the adaptive immune response. Balancing selection is believed to maintain MHC polymorphism in the long term, although neutral processes also play a role in shaping MHC diversity. However, the relative contribution of these processes is poorly understood. Here we characterized MHC class II variation of a low-dispersal, pond-breeding newt (Triturus carnifex) over a restricted, geographically structured area. We aimed to (1) evaluate the contribution of selection and neutral processes to shaping MHC diversity at two geographic scales, and (2) test for signatures of divergent allele advantage (DAA), which is a potentially important mechanism of balancing selection. The dominant role of selection in shaping MHC variation was suggested by the lack of correlation between MHC and neutral (microsatellite) variation. Although most variation occurred within populations for both types of markers, they differed in the extent of structuring at the two spatial scales. MHC structuring was more pronounced at local scales, suggesting the role of local selection, while structuring was not detectable at a larger scale, possibly due to the effect of balancing selection. Microsatellites showed the opposite pattern. As expected under DAA, the observed genotypes combined more sequence diversity than expected under a random association of alleles. Thus, DAA may contribute to maintaining MHC polymorphism, which is ancient, as supported by signatures of historical positive selection and trans-species polymorphism. Our results point to the importance of a multi-scale approach in studying MHC variation, especially in low-dispersal taxa, which are genetically structured at fine spatial scales.

Which Latitudinal Gradients for Genetic Diversity?

A recent global analysis of GenBank DNA sequences from amphibians and mammals indicated consistent poleward decrease of intraspecific genetic diversity in both classes. We highlight that this result was biased by not accounting for distance decay of similarity and reanalyse the datasets, revealing distinct latitudinal gradients in mammals and amphibians.

From pest data to abundance-based risk maps combining eco-physiological knowledge, weather, and habitat variability.

Noxious species, i.e., crop pest or invasive alien species, are major threats to both natural and managed ecosystems. Invasive pests are of special importance, and knowledge about their distribution and abundance is fundamental to minimize economic losses and prioritize management activities. Occurrence models are a common tool used to identify suitable zones and map priority areas (i.e., risk maps) for noxious species management, although they provide a simplified description of species dynamics (i.e., no indication on species density). An alternative is to use abundance models, but translating abundance data into risk maps is often challenging. Here, we describe a general framework for generating abundance-based risk maps using multi-year pest data. We used an extensive data set of 3968 records collected between 2003 and 2013 in Wisconsin during annual surveys of soybean aphid (SBA), an exotic invasive pest in this region. By using an integrative approach, we modelled SBA responses to weather, seasonal, and habitat variability using generalized additive models (GAMs). Our models showed good to excellent performance in predicting SBA occurrence and abundance (TSS = 0.70, AUC = 0.92; R2  = 0.63). We found that temperature, precipitation, and growing degree days were the main drivers of SBA trends. In addition, a significant positive relationship between SBA abundance and the availability of overwintering habitats was observed. Our models showed aphid populations were also sensitive to thresholds associated with high and low temperatures, likely related to physiological tolerances of the insects. Finally, the resulting aphid predictions were integrated using a spatial prioritization algorithm ("Zonation") to produce an abundance-based risk map for the state of Wisconsin that emphasized the spatiotemporal consistency and magnitude of past infestation patterns. This abundance-based risk map can provide information on potential foci of pest outbreaks where scouting efforts and prophylactic measures should be concentrated. The approach we took is general, relatively simple, and can be applied to other species, habitats and geographical areas for which species abundance data and biotic and abiotic data are available.

Testing Classical Species Properties with Contemporary Data: How "Bad Species" in the Brassy Ringlets (Erebia tyndarus complex, Lepidoptera) Turned Good.

All species concepts are rooted in reproductive, and ultimately genealogical, relations. Genetic data are thus the most important source of information for species delimitation. Current ease of access to genomic data and recent computational advances are blooming a plethora of coalescent-based species delimitation methods. Despite their utility as objective approaches to identify species boundaries, coalescent-based methods (1) rely on simplified demographic models that may fail to capture some attributes of biological species, (2) do not make explicit use of the geographic information contained in the data, and (3) are often computationally intensive. In this article, we present a case of species delimitation in the Erebia tyndarus species complex, a taxon regarded as a classic example of problematic taxonomic resolution. Our approach to species delimitation used genomic data to test predictions rooted in the biological species concept and in the criterion of coexistence in sympatry. We (1) obtained restriction-site associated DNA (RAD) sequencing data from a carefully designed sample, (2) applied two genotype clustering algorithms to identify genetic clusters, and (3) performed within-clusters and between-clusters analyses of isolation by distance as a test for intrinsic reproductive barriers. Comparison of our results with those from a Bayes factor delimitation coalescent-based analysis, showed that coalescent-based approaches may lead to overconfident splitting of allopatric populations, and indicated that incorrect species delimitation is likely to be inferred when an incomplete geographic sample is analyzed. While we acknowledge the theoretical justification and practical usefulness of coalescent-based species delimitation methods, our results stress that, even in the phylogenomic era, the toolkit for species delimitation should not dismiss more traditional, biologically grounded, approaches coupling genomic data with geographic information.