Deciphering microbial communities of three Savoyard raw milk cheeses along ripening and regarding the cheese process.

Different Savoyard cheeses are granted with PDO (Protected Designation or Origin) and PGI (Protected Geographical Indication) which guarantees consumers compliance with strict specifications. The use of raw milk is known to be crucial for specific flavor development. To unravel the factors influencing microbial ecosystems across cheese making steps, according to the seasonality (winter and summer) and the mode of production (farmhouse and dairy factory ones), gene targeting on bacteria and fungus was used to have a full picture of 3 cheese making technologies, from the raw milk to the end of the ripening. Our results revealed that Savoyard raw milks are a plenteous source of biodiversity together with the brines used during the process, that may support the development of specific features for each cheese. It was shown that rinds and curds have very contrasted ecosystem diversity, composition, and evolution. Ripening stage was selective for some bacterial species, whereas fungus were mainly ubiquitous in dairy samples. All ripening stages are impacted by the type of cheese technologies, with a higher impact on bacterial communities, except for fungal rind communities, for which the technology is the more discriminant. The specific microorganism's abundance for each technology allow to see a real bar-code, with more or less differences regarding bacterial or fungal communities. Bacterial structuration is shaped mainly by matrices, differently regarding technologies while the influence of technology is higher for fungi. Production types showed 10 differential bacterial species, farmhouses showed more ripening taxa, while dairy factory products showing more lactic acid bacteria. Meanwhile, seasonality looks to be a minor element for the comprehension of both microbial ecosystems, but the uniqueness of each dairy plant is a key explicative feature, more for bacteria than for fungus communities.

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.

Differential effects of tree species identity on rhizospheric bacterial and fungal community richness and composition across multiple trace element-contaminated sites.

Soil microbial communities play a key role in plant nutrition and stress tolerance. This is particularly true in sites contaminated by trace metals, which often have low fertility and stressful conditions for woody plants in particular. However, we have limited knowledge of the abiotic and biotic factors affecting the richness and composition of microbial communities inhabiting the rhizosphere of plants in contaminated sites. Using high-throughput amplicon sequencing, we studied the rhizospheric bacterial and fungal community structures of 14 woody plant families planted in three contrasting sites contaminated by metals (Pb, Cd, Zn, Mn, Fe, S). The rhizospheric bacterial communities in the given sites showed no significant difference between the various woody species but did differ significantly between sites. The Proteobacteria phylum was dominant, accounting for over 25 % of the overall relative abundance, followed by Actinobacteria, Bacteroidetes and Gemmatimonadetes. Site was also the main driver of fungal community composition, yet unlike bacteria, tree species identity significantly affected fungal communities. The Betulaceae, Salicaceae and Fagaceae families had a high proportion of Basidiomycota, particularly ectomycorrhizal fungi, and the lowest diversity and richness. The other tree families and the unplanted soil harboured a greater abundance of Ascomycota and Mucoromycota. Consequently, for both bacteria and fungi, the site effect significantly impacted their community richness and composition, while the influence of plants on the richness and composition of rhizospheric microbial communities stayed consistent across sites and was dependent on the microbial kingdom. Finally, we highlighted the importance of considering this contrasting response of plant rhizospheric microbial communities in relation to their host identity, particularly to improve assisted revegetation efforts at contaminated sites.

Individual genotypes from environmental DNA: Fingerprinting snow tracks of three large carnivore species.

Continued advancements in environmental DNA (eDNA) research have made it possible to access intraspecific variation from eDNA samples, opening new opportunities to expand non-invasive genetic studies of wildlife populations. However, the use of eDNA samples for individual genotyping, as typically performed in non-invasive genetics, still remains elusive. We present successful individual genotyping of eDNA obtained from snow tracks of three large carnivores: brown bear (Ursus arctos), European lynx (Lynx lynx) and wolf (Canis lupus). DNA was extracted using a protocol for isolating water eDNA and genotyped using amplicon sequencing of short tandem repeats (STR), and for brown bear a sex marker, on a high-throughput sequencing platform. Individual genotypes were obtained for all species, but genotyping performance differed among samples and species. The proportion of samples genotyped to individuals was higher for brown bear (5/7) and wolf (7/10) than for lynx (4/9), and locus genotyping success was greater for brown bear (0.88). The sex marker was typed in six out of seven brown bear samples. Results for three species show that reliable individual genotyping, including sex identification, is now possible from eDNA in snow tracks, underlining its vast potential to complement the non-invasive genetic methods used for wildlife. To fully leverage the application of snow track eDNA, improved understanding of the ideal species- and site-specific sampling conditions, as well as laboratory methods promoting genotyping success, is needed. This will also inform efforts to retrieve and type nuclear DNA from other eDNA samples, thereby advancing eDNA-based individual and population-level studies.

Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints.

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

Towards exhaustive community ecology via DNA metabarcoding.

Exhaustive biodiversity data, covering all the taxa in an environment, would be fundamental to understand how global changes influence organisms living at different trophic levels, and to evaluate impacts on interspecific interactions. Molecular approaches such as DNA metabarcoding are boosting our ability to perform biodiversity inventories. Nevertheless, even though a few studies have recently attempted exhaustive reconstructions of communities, holistic assessments remain rare. The majority of metabarcoding studies published in the last years used just one or two markers and analysed a limited number of taxonomic groups. Here, we provide an overview of emerging approaches that can allow all-taxa biological inventories. Exhaustive biodiversity assessments can be attempted by combining a large number of specific primers, by exploiting the power of universal primers, or by combining specific and universal primers to obtain good information on key taxa while limiting the overlooked biodiversity. Multiplexes of primers, shotgun sequencing and capture enrichment may provide a better coverage of biodiversity compared to standard metabarcoding, but still require major methodological advances. Here, we identify the strengths and limitations of different approaches, and suggest new development lines that might improve broad scale biodiversity analyses in the near future. More holistic reconstructions of ecological communities can greatly increase the value of metabarcoding studies, improving understanding of the consequences of ongoing environmental changes on the multiple components of biodiversity.

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.

Focal vs. fecal: Seasonal variation in the diet of wild vervet monkeys from observational and DNA metabarcoding data.

Assessing the diet of wild animals reveals valuable information about their ecology and trophic relationships that may help elucidate dynamic interactions in ecosystems and forecast responses to environmental changes. Advances in molecular biology provide valuable research tools in this field. However, comparative empirical research is still required to highlight strengths and potential biases of different approaches. Therefore, this study compares environmental DNA and observational methods for the same study population and sampling duration. We employed DNA metabarcoding assays targeting plant and arthropod diet items in 823 fecal samples collected over 12 months in a wild population of an omnivorous primate, the vervet monkey (Chlorocebus pygerythrus). DNA metabarcoding data were subsequently compared to direct observations. We observed the same seasonal patterns of plant consumption with both methods; however, DNA metabarcoding showed considerably greater taxonomic coverage and resolution compared to observations, mostly due to the construction of a local plant DNA database. We found a strong effect of season on variation in plant consumption largely shaped by the dry and wet seasons. The seasonal effect on arthropod consumption was weaker, but feeding on arthropods was more frequent in spring and summer, showing overall that vervets adapt their diet according to available resources. The DNA metabarcoding assay outperformed also direct observations of arthropod consumption in both taxonomic coverage and resolution. Combining traditional techniques and DNA metabarcoding data can therefore not only provide enhanced assessments of complex diets and trophic interactions to the benefit of wildlife conservationists and managers but also opens new perspectives for behavioral ecologists studying whether diet variation in social species is induced by environmental differences or might reflect selective foraging behaviors.

The generality of cryptic dietary niche differences in diverse large-herbivore assemblages.

Ecological niche differences are necessary for stable species coexistence but are often difficult to discern. Models of dietary niche differentiation in large mammalian herbivores invoke the quality, quantity, and spatiotemporal distribution of plant tissues and growth forms but are agnostic toward food plant species identity. Empirical support for these models is variable, suggesting that additional mechanisms of resource partitioning may be important in sustaining large-herbivore diversity in African savannas. We used DNA metabarcoding to conduct a taxonomically explicit analysis of large-herbivore diets across southeastern Africa, analyzing ∼4,000 fecal samples of 30 species from 10 sites in seven countries over 6 y. We detected 893 food plant taxa from 124 families, but just two families-grasses and legumes-accounted for the majority of herbivore diets. Nonetheless, herbivore species almost invariably partitioned food plant taxa; diet composition differed significantly in 97% of pairwise comparisons between sympatric species, and dissimilarity was pronounced even between the strictest grazers (grass eaters), strictest browsers (nongrass eaters), and closest relatives at each site. Niche differentiation was weakest in an ecosystem recovering from catastrophic defaunation, indicating that food plant partitioning is driven by species interactions, and was stronger at low rainfall, as expected if interspecific competition is a predominant driver. Diets differed more between browsers than grazers, which predictably shaped community organization: Grazer-dominated trophic networks had higher nestedness and lower modularity. That dietary differentiation is structured along taxonomic lines complements prior work on how herbivores partition plant parts and patches and suggests that common mechanisms govern herbivore coexistence and community assembly in savannas.

Interspecific coprophagia by wild red foxes: DNA metabarcoding reveals a potentially widespread form of commensalism among animals.

Vertebrate animals are known to consume other species' faeces, yet the role of such coprophagy in species dynamics remains unknown, not least due to the methodological challenges of documenting it. In a large-scale metabarcoding study of red fox and pine marten scats, we document a high occurrence of domestic dog DNA in red fox scats and investigate if it can be attributed to interspecific coprophagia. We tested whether experimental artifacts or other sources of DNA could account for dog DNA, regressed dog occurrence in the diet of fox against that of the fox's main prey, short-tailed field voles, and consider whether predation or scavenging could explain the presence of dog DNA. Additionally, we determined the calorific value of dog faeces through calorimetric explosion. The high occurrence of dog DNA in the diet of fox, the timing of its increase, and the negative relationship between dog and the fox's main prey, point to dog faeces as the source of DNA in fox scats. Dog faeces being highly calorific, we found that foxes, but not pine martens, regularly exploit them, seemingly as an alternative resource to fluctuating prey. Scattered accounts from the literature may suggest that interspecific coprophagia is a potentially frequent and widespread form of interaction among vertebrates. However, further work should address its prevalence in other systems and the implications for ecological communities. Tools such as metabarcoding offer a way forward.

Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros.

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Analysis of complex trophic networks reveals the signature of land-use intensification on soil communities in agroecosystems.

Increasing evidence suggests that agricultural intensification is a threat to many groups of soil biota, but how the impacts of land-use intensity on soil organisms translate into changes in comprehensive soil interaction networks remains unclear. Here for the first time, we use environmental DNA to examine total soil multi-trophic diversity and food web structure for temperate agroecosystems along a gradient of land-use intensity. We tested for response patterns in key properties of the soil food webs in sixteen fields ranging from arable crops to grazed permanent grasslands as part of a long-term management experiment. We found that agricultural intensification drives reductions in trophic group diversity, although taxa richness remained unchanged. Intensification generally reduced the complexity and connectance of soil interaction networks and induced consistent changes in energy pathways, but the magnitude of management-induced changes depended on the variable considered. Average path length (an indicator of food web redundancy and resilience) did not respond to our management intensity gradient. Moreover, turnover of network structure showed little response to increasing management intensity. Our data demonstrates the importance of considering different facets of trophic networks for a clearer understanding of agriculture-biodiversity relationships, with implications for nature-based solutions and sustainable agriculture.

Assessing changes in stream macroinvertebrate communities across ecological gradients using morphological versus DNA metabarcoding approaches.

Freshwater macroinvertebrates provide valuable indicators for biomonitoring ecosystem change in relation to natural and anthropogenic drivers. DNA metabarcoding is an efficient approach for estimating such indicators, but its results may differ from morphotaxonomic approaches traditionally used in biomonitoring. Here we test the hypothesis that despite differences in the number and identity of taxa recorded, both approaches may retrieve comparable patterns of community change, and detect similar ecological gradients influencing such changes. We compared results obtained with morphological identification at family level of macroinvertebrates collected at 80 streams under a Water Framework Directive biomonitoring program in Portugal, with results obtained with metabarcoding from the ethanol preserving the bulk samples, using either single (COI-M19BR2, 16S-Inse01, 18S-Euka02) or multiple markers. Metabarcoding recorded less families and different communities compared to morphotaxonomy, but community sensitivities to disturbance estimated with the IASPT index were more similar across approaches. Spatial variation in local community metrics and the factors influencing such variation were significantly correlated between morphotaxonomy and metabarcoding. After reducing random noise in the dissimilarity matrices, the spatial variation in community composition was also significantly correlated across methods. A dominant gradient of community change was consistently retrieved, and all methods identified a largely similar set of anthropogenic stressors strongly influencing such gradient. Overall, results confirm our initial hypothesis, suggesting that morphotaxonomy and metabarcoding can estimate consistent spatial patterns of community variation and their main drivers. These results are encouraging for macroinvertebrate biomonitoring using metabarcoding approaches, suggesting that they can be intercalibrated with morphotaxonomic approaches to recover equivalent spatial and temporal gradients of ecological change.