Natural reservoirs of zoonotic pathogens generally seem to be capable of tolerating infections. Tolerance and its underlying mechanisms remain difficult to assess using experiments or wildlife surveys. High-throughput sequencing technologies give the opportunity to investigate the genetic bases of tolerance, and the variability of its mechanisms in natural populations. In particular, population genomics may provide preliminary insights into the genes shaping tolerance and potentially influencing epidemiological dynamics. Here, we addressed these questions in the bank vole Myodes glareolus, the specific asymptomatic reservoir host of Puumala hantavirus (PUUV), which causes nephropathia epidemica (NE) in humans. Despite the continuous spatial distribution of M. glareolus in Sweden, NE is endemic to the northern part of the country. Northern bank vole populations in Sweden might exhibit tolerance strategies as a result of coadaptation with PUUV. This may favor the circulation and maintenance of PUUV and lead to high spatial risk of NE in northern Sweden. We performed a genome-scan study to detect signatures of selection potentially correlated with spatial variations in tolerance to PUUV. We analyzed six bank vole populations from Sweden, sampled from northern NE-endemic to southern NE-free areas. We combined candidate gene analyses (Tlr4, Tlr7, and Mx2 genes) and high-throughput sequencing of restriction site-associated DNA (RAD) markers. Outlier loci showed high levels of genetic differentiation and significant associations with environmental data including variations in the regional number of NE human cases. Among the 108 outliers that matched to mouse protein-coding genes, 14 corresponded to immune-related genes. The main biological pathways found to be significantly enriched corresponded to immune processes and responses to hantavirus, including the regulation of cytokine productions, TLR cascades, and IL-7, VEGF, and JAK-STAT signaling. In the future, genome-scan replicates and functional experimentations should enable to assess the role of these biological pathways in M. glareolus tolerance to PUUV.

Infectious pathogens are major selective forces acting on individuals. The recent advent of high-throughput sequencing technologies now enables to investigate the genetic bases of resistance/susceptibility to infections in non-model organisms. From an evolutionary perspective, the analysis of the genetic diversity observed at these genes in natural populations provides insight into the mechanisms maintaining polymorphism and their epidemiological consequences. We explored these questions in the context of the interactions between Puumala hantavirus (PUUV) and its reservoir host, the bank vole Myodes glareolus . Despite the continuous spatial distribution of M. glareolus in Europe, PUUV distribution is strongly heterogeneous. Different defence strategies might have evolved in bank voles as a result of co-adaptation with PUUV, which may in turn reinforce spatial heterogeneity in PUUV distribution. We performed a genome scan study of six bank vole populations sampled along a North/South transect in Sweden, including PUUV endemic and non-endemic areas. We combined candidate gene analyses ( Tlr4 , Tlr7 , Mx2 genes) and high throughput sequencing of RAD (Restriction-site Associated DNA) markers. We found evidence for outlier loci showing high levels of genetic differentiation. Ten outliers among the 52 that matched to mouse protein-coding genes corresponded to immune related genes and were detected using ecological associations with variations in PUUV prevalence. One third of the enriched pathways concerned immune processes, including platelet activation and TLR pathway. In the future, functional experimentations should enable to confirm the role of these these immune related genes with regard to the interactions between M. glareolus and PUUV.

Crassostrea gigas originated from the Pacific coast of Asia, but was introduced into several European countries in the early 1970s. Natural populations have now spread across the length of the western seaboard of Europe. To elucidate the demographic and selective processes at play during this rapid expansion, genome-scan analysis was performed on different populations. High diversities and low differentiation were observed overall, but significant genetic differentiation was found among newly established populations and between the newly established northern group and a nearly panmictic group composed of southern European populations and a population from Japan. Loss of genetic diversity was also seen in the north, likely caused by founder events during colonization. The few strongly supported outlier loci revealed a genetic structure uncorrelated with the north/south differentiation, but grouping two samples from the Danish fjords (northern group) and one from the Dutch Scheldt estuary (southern group) with the one from Japan. These findings might reflect the following: (i) parallel adaptation to similar environmental pressures (fjord-like environment) within each of the two groups or (ii) a footprint of a secondary introduction of an alternative genomic background maintained by multifarious isolation factors. Our results call for a closer examination of adaptive genetic structure in the area of origin.

Le campagnol roussâtre Myodes glareolus est le réservoir de l’hantavirus Puumala (PUUV), responsable chez l’Homme d’une forme atténuée de Fièvre Hémorragique à Syndrome Rénal (FHSR), la Néphropathie Épidémique (NE). En Europe, l’incidence de la NE présente, malgré la distribution continue du réservoir, une forte variabilité géographique dont les causes ne sont à ce jour pas identifiées. Aux hypothèses climatiques et paysagères, nous proposons que des facteurs intrinsèques aux campagnols puissent également être impliqués. Une plus forte tolérance à l’infection par le virus PUUV, chez certains campagnols roussâtres, favoriserait la persistance et la transmission de ce virus, ce qui devrait accroître le risque de NE chez l’Homme. Nous avons testé cette hypothèse en étudiant les polymorphismes et/ou les niveaux d''expression de six gènes candidats impliqués dans la réponse immunitaire à PUUV (DRB-CMH, promoteur du TNF-, TLR4, TLR7, Mx2, intégrine 3) chez dix populations de campagnols échantillonnées le long d’un axe nord/ sud dans les Ardennes françaises, couvrant des zones endémiques et non endémiques à PUUV. Des signatures de sélection ont été détectées pour TNF-et Mx2 grâce à des approches de génétique des populations (scan FST) et d’associations génotypes/phénotypes. Ces gènes codent des protéines dont les propriétés antivirales sont connues, mais qui induisent des coûts immunopathologiques importants. Ils pourraient donc jouer un rôle central dans une balance de tolérance/résistance à PUUV. De plus, dans les zones d’endémie, les gènes TNF-et Mx2 sont sous-exprimés, ce qui suggère l’évolution d’une plus forte tolérance à PUUV, potentiellement au bénéfice d’un moindre coût immunopathologique. Certains de ces résultats ont été confirmés à l’échelle européenne.

Hybrid zones are noteworthy systems for the study of environmental adaptation to fast-changing environments, as they constitute reservoirs of polymorphism and are key to the maintenance of biodiversity. They can move in relation to climate fluctuations, as temperature can affect both selection and migration, or remain trapped by environmental and physical barriers. There is therefore a very strong incentive to study the dynamics of hybrid zones subjected to climate variations. The infaunal bivalve Macoma balthica emerges as a noteworthy model species, as divergent lineages hybridize, and its native NE Atlantic range is currently contracting to the North. To investigate the dynamics and functioning of hybrid zones in M. balthica, we developed new molecular markers by sequencing the collective transcriptome of 30 individuals. Ten individuals were pooled for each of the three populations sampled at the margins of two hybrid zones. A single 454 run generated 277 Mb from which 17K SNPs were detected. SNP density averaged 1 polymorphic site every 14 to 19 bases, for mitochondrial and nuclear loci, respectively. An scan detected high genetic divergence among several hundred SNPs, some of them involved in energetic metabolism, cellular respiration and physiological stress. The high population differentiation, recorded for nuclear-encoded ATP synthase and NADH dehydrogenase as well as most mitochondrial loci, suggests cytonuclear genetic incompatibilities. Results from this study will help pave the way to a high-resolution study of hybrid zone dynamics in M. balthica, and the relative importance of endogenous and exogenous barriers to gene flow in this system.