High rates of gene duplication and the highest levels of functional allelic diversity in vertebrate genomes are the main hallmarks of the major histocompatibility complex (MHC), a multigene family with a primordial role in pathogen recognition. The usual tight linkage among MHC gene duplicates may provide an opportunity for the evolution of haplotypes that associate functionally divergent alleles and thus grant the transmission of optimal levels of diversity to coming generations. Even though such associations may be a crucial component of disease resistance, this hypothesis has been given little attention in wild populations. Here, we leveraged pedigree data from a barn owl (Tyto alba) population to characterize MHC haplotype structure across two MHC class I (MHC-I) and two MHC class IIB (MHC-IIB) duplicates, in order to test the hypothesis that haplotypes’ genetic diversity is higher than expected from randomly associated alleles. After showing that MHC loci are tightly linked within classes, we found limited evidence for shifts towards MHC haplotypes combining high diversity. Neither amino acid nor functional within-haplotype diversity were significantly higher than in random sets of haplotypes, regardless of MHC class. Our results therefore provide no evidence for selection towards high-diversity MHC haplotypes in barn owls. Rather, high rates of convergent evolution may constrain the evolution of high-diversity haplotypes at MHC-I, while, in contrast, for MHC-IIB, fixed differences among loci may provide barn owls with already optimized functional diversity. This suggests that at the MHC-I and MHC-IIB, respectively, different evolutionary dynamics may govern the evolution of within-haplotype diversity.

In the last few years, the study of environmental DNA (eDNA) has drawn attention for many reasons, including its advantages for monitoring and conservation purposes. So far, in aquatic environments, most of eDNA research has focused on the detection of single species using species-specific markers. Recently, species inventories based on the analysis of a single generalist marker targeting a larger taxonomic group (eDNA metabarcoding) have proven useful for bony fish and amphibian biodiversity surveys. This approach involves in situ filtering of large volumes of water followed by amplification and sequencing of a short discriminative fragment from the 12S rDNA mitochondrial gene. In this study, we went one step further by investigating the spatial representativeness (i.e. ecological reliability and signal variability in space) of eDNA metabarcoding for large-scale fish biodiversity assessment in a freshwater system including lentic and lotic environments. We tested the ability of this approach to characterize large-scale organization of fish communities along a longitudinal gradient, from a lake to the outflowing river. First, our results confirm that eDNA metabarcoding is more efficient than a single traditional sampling campaign to detect species presence, especially in rivers. Second, the species list obtained using this approach is comparable to the one obtained when cumulating all traditional sampling sessions since 1995 and 1988 for the lake and the river, respectively. In conclusion, eDNA metabarcoding gives a faithful description of local fish biodiversity in the study system, more specifically within a range of a few kilometers along the river in our study conditions, i.e. longer than a traditional fish sampling site.

With their direct link to individual fitness, genes of the major histocompatibility complex (MHC) are a popular system to study the evolution of adaptive genetic diversity. However, owing to the highly dynamic evolution of the MHC region, the isolation, characterization and genotyping of MHC genes remain a major challenge. While high-throughput sequencing technologies now provide unprecedented resolution of the high allelic diversity observed at the MHC, in many species, it remains unclear (i) how alleles are distributed among MHC loci, (ii) whether MHC loci are linked or segregate independently and (iii) how much copy number variation (CNV) can be observed for MHC genes in natural populations. Here, we show that the study of allele segregation patterns within families can provide significant insights in this context. We sequenced two MHC class I (MHC-I) loci in 1267 European barn owls (Tyto alba), including 590 offspring from 130 families using Illumina MiSeq technology. Coupled with a high per-individual sequencing coverage (~3000×), the study of allele segregation patterns within families provided information on three aspects of the architecture of MHC-I variation in barn owls: (i) extensive sharing of alleles among loci, (ii) strong linkage of MHC-I loci indicating tandem architecture and (iii) the presence of CNV in the barn owl MHC-I. We conclude that the additional information that can be gained from high-coverage amplicon sequencing by investigating allele segregation patterns in families not only helps improving the accuracy of MHC genotyping, but also contributes towards enhanced analyses in the context of MHC evolutionary ecology.

Microbial community analysis was carried out on ruminal digesta obtained directly via rumen fistula and buccal fluid, regurgitated digesta (bolus) and faeces of dairy cattle to assess if non-invasive samples could be used as proxies for ruminal digesta. Samples were collected from five cows receiving grass silage based diets containing no additional lipid or four different lipid supplements in a 5 x 5 Latin square design. Extracted DNA was analysed by qPCR and by sequencing 16S and 18S rRNA genes or the fungal ITS1 amplicons. Faeces contained few protozoa, and bacterial, fungal and archaeal communities were substantially different to ruminal digesta. Buccal and bolus samples gave much more similar profiles to ruminal digesta, although fewer archaea were detected in buccal and bolus samples. Bolus samples overall were most similar to ruminal samples. The differences between both buccal and bolus samples and ruminal digesta were consistent across all treatments. It can be concluded that either proxy sample type could be used as a predictor of the rumen microbial community, thereby enabling more convenient large-scale animal sampling for phenotyping and possible use in future animal breeding programs aimed at selecting cattle with a lower environmental footprint.

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

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

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

In tropical regions, most primary ecosystems have been replaced by mosaic landscapes in which species must cope with a large shift in the distribution of their habitat and associated food resources. Primates are particularly vulnerable to habitat modifications. Most species persist in small fragments surrounded by complex human-mediated matrices whose structure and connectivity may strongly influence their dispersal and feeding behavior. Behavioral plasticity appears to be a crucial parameter governing the ability of organisms to exploit the resources offered by new matrix habitats and thus to persist in fragmented habitats. In this study, we were interested in the dietary plasticity of the golden-crowned sifaka (Propithecus tattersalli), an endangered species of lemur, found only in the Daraina region in north-eastern Madagascar. We used a DNA-based approach combining the barcoding concept and Illumina next-generation sequencing to (i) describe the species diet across its entire range and (ii) evaluate the influence of landscape heterogeneity on diet diversity and composition. Faeces from 96 individuals were sampled across the entire species range and their contents were analyzed using the trnL metabarcoding approach. In parallel, we built a large DNA reference database based on a checklist of the plant species of the Daraina region. Our results suggest that golden-crowned sifakas exhibit remarkable dietary diversity with at least 130 plant species belonging to 80 genera and 49 different families. We highlighted an influence of both habitat type and openness on diet composition suggesting a high flexibility of foraging strategies. Moreover, we observed the presence of numerous cultivated and naturalized plants in the faeces of groups living in forest edge areas. Overall, our findings support our initial expectation that P. tattersalli is able to cope with the current level of alteration of the landscape and confirm our previous results on the distribution and the dispersal ability of this species.

Species range shifts in response to climate and land use change are commonly forecasted with species distribution models based on species occurrence or abundance data. Although appealing, these models ignore the genetic structure of species, and the fact that different populations might respond in different ways due to adaptation to their environment. Here, we introduced ancestry distribution models, i.e., statistical models of the spatial distribution of ancestry proportions, for forecasting intra-specific changes based on genetic admixture instead of species occurrence data. Using multi-locus genotypes and extensive geographic coverage of distribution data across the European Alps, we applied this approach to 20 alpine plant species considering a global increase in temperature from 0.25°C to 4°C. We forecasted the magnitudes of displacement of contact zones between plant populations potentially adapted to warmer environments and other populations. While a global trend of movement in a northeast direction was predicted, the magnitude of displacement was species-specific. For a temperature increase of 2°C, contact zones were predicted to move by 92 km on average (minimum of 5 km, maximum of 212 km), and by 188 km for an increase of 4°C (minimum of 11 km, maximum of 393 km). Intra-specific turnover – measuring the extent of change in global population genetic structure – was generally found to be moderate for 2°C of temperature warming. For 4°C of warming, however, the models indicated substantial intra-specific turnover for ten species. These results illustrate that, in spite of unavoidable simplifications, ancestry distribution models open new perspectives to forecast population genetic changes within species, and complement more traditional distribution-based approaches.