Pathogen burden, co-infection and major histocompatibility complex variability in the European badger (Meles meles)

Changes in capture availability due to infection can lead to detectable biases in population-level infectious disease parameters

Correctly identifying the strength of selection that parasites impose on hosts is key to predicting epidemiological and evolutionary outcomes of host-parasite interactions. However, behavioral changes due to infection can alter the capture probability of infected hosts and thereby make selection difficult to estimate by standard sampling techniques. Mark-recapture approaches, which allow researchers to determine if some groups in a population are less likely to be captured than others, can be used to identify infection-driven capture biases. If a metric of interest directly compares infected and uninfected populations, calculated detection probabilities for both groups may be useful in identifying bias. Here, we use an individual-based simulation to test whether changes in capture rate due to infection can alter estimates of three key metrics: 1) reduction in the reproductive success of infected parents relative to uninfected parents, 2) the relative risk of infection for susceptible genotypes compared to resistant genotypes, and 3) changes in allele frequencies between generations. We explore the direction and underlying causes of the biases that emerge from these simulations. Finally, we argue that short series of mark-recapture sampling bouts, potentially implemented in under a week, can yield key data on detection bias due to infection while not adding a significantly higher burden to disease ecology studies.

Immunogenetics, sylvatic plague and its vectors: insights from the pathogen reservoir Mastomys natalensis in Tanzania

Yersinia pestis is a historically important vector-borne pathogen causing plague in humans and other mammals. Contemporary zoonotic infections with Y. pestis still occur in sub-Saharan Africa, including Tanzania and Madagascar, but receive relatively little attention. Thus, the role of wildlife reservoirs in maintaining sylvatic plague and spillover risks to humans is largely unknown. The multimammate rodent Mastomys natalensis is the most abundant and widespread rodent in peri-domestic areas in Tanzania, where it plays a major role as a Y. pestis reservoir in endemic foci. Yet, how M. natalensis' immunogenetics contributes to the maintenance of plague has not been investigated to date. Here, we surveyed wild M. natalensis for Y. pestis vectors, i.e., fleas, and tested for the presence of antibodies against Y. pestis using enzyme-linked immunosorbent assays (ELISA) in areas known to be endemic or without previous records of Y. pestis in Tanzania. We characterized the allelic and functional (i.e., supertype) diversity of the major histocompatibility complex (MHC class II) of M. natalensis and investigated links to Y. pestis vectors and infections. We detected antibodies against Y. pestis in rodents inhabiting both endemic areas and areas considered non-endemic. Of the 111 nucleotide MHC alleles, only DRB*016 was associated with an increased infestation with the flea Xenopsylla. Surprisingly, we found no link between MHC alleles or supertypes and antibodies of Y. pestis. Our findings hint, however, at local adaptations towards Y. pestis vectors, an observation that more exhaustive sampling could unwind in the future.

MHC-B Diversity and Signs of Respiratory Illness in Wild, East African Chimpanzees ( Pan troglodytes schweinfurthii )

Many traits, intrinsic and extrinsic to an organism, contribute to interindividual variation in immunity in wild habitats. The vertebrate Major Histocompatibility Complex (MHC) includes genes encoding antigen-presenting molecules that are highly variable, and that variation often predicts susceptibility/resistance to and recovery from pathogen infection. I compare MHC-B variation at two long-term chimpanzee research sites, Kibale National Park in Uganda and Gombe National Park in Tanzania. Using decades of respiratory health data available for these chimpanzees, I test hypotheses associated with maintenance of diversity at MHC loci, including heterozygote, divergent allele, and rare allele advantage hypotheses, and predictions for unique function of MHC-B in great apes. I found, despite confirmation of recent shared ancestry between Kibale and Gombe chimpanzees, including an overlapping MHC-B allele repertoire and similar MHC-B phenotype compositions, chimpanzees from the two research sites experienced differences in the occurrence of respiratory signs and had different associations of MHC-B diversity with signs of respiratory illness. Kibale chimpanzees with heterozygous genotypes and different peptide-binding supertypes were observed less often with respiratory signs than those homozygous or possessing the same supertypes, but this same association was not observed among Gombe chimpanzees. Gombe chimpanzees with specific MHC-B phenotypes that enable engagement of Natural Killer (NK) cells were observed more often with respiratory signs than chimpanzees with other phenotypes, but this was not observed at Kanyawara. This study emphasizes local adaptation in shaping genetic and phenotypic traits in different infectious disease contexts, even among close genetic relatives of the same subspecies, and highlights utility for continued and simultaneous tracking of host immune genes and specific pathogens in wild species.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Assessing chronic stress in wild mammals using claw-derived cortisol: a validation using European badgers (Meles meles)

Measuring stress experienced by wild mammals is increasingly important in the context of human-induced rapid environmental change and initiatives to mitigate human-wildlife conflicts. Glucocorticoids (GC), such as cortisol, mediate responses by promoting physiological adjustments during environmental perturbations. Measuring cortisol is a popular technique; however, this often reveals only recent short-term stress such as that incurred by restraining the animal to sample blood, corrupting the veracity of this approach. Here we present a protocol using claw cortisol, compared with hair cortisol, as a long-term stress bio-indicator, which circumvents this constraint, where claw tissue archives the individual's GC concentration over preceding weeks. We then correlate our findings against detailed knowledge of European badger life history stressors. Based on a solid-phase extraction method, we assessed how claw cortisol concentrations related to season and badger sex, age and body-condition using a combination of generalized linear mixed models (GLMM) (n = 668 samples from 273 unique individuals) followed by finer scale mixed models for repeated measures (MMRM) (n = 152 re-captured individuals). Claw and hair cortisol assays achieved high accuracy, precision and repeatability, with similar sensitivity. The top GLMM model for claw cortisol included age, sex, season and the sex*season interaction. Overall, claw cortisol levels were significantly higher among males than females, but strongly influenced by season, where females had higher levels than males in autumn. The top fine scale MMRM model included sex, age and body condition, with claw cortisol significantly higher in males, older and thinner individuals. Hair cortisol was more variable than claw; nevertheless, there was a positive correlation after removing 34 outliers. We discuss strong support for these stress-related claw cortisol patterns from previous studies of badger biology. Given the potential of this technique, we conclude that it has broad application in conservation biology.

Infection with a Recently Discovered Gammaherpesvirus Variant in European Badgers, Meles meles, is Associated with Higher Relative Viral Loads in Blood

Herpesviruses are ubiquitous pathogens infecting most animals. Although host immunity continually coevolves to combat virulence, viral variants with enhanced transmissibility or virulence occasionally emerge, resulting in disease burdens in host populations. Mustelid gammaherpesvirus 1 (MusGHV-1) is the only herpesvirus species identified thus far in European badgers, Meles meles. No MusGHV-1 associated pathomorbidity has been reported, but reactivation of MusGHV-1 in genital tracts is linked to impaired female reproductive success. An analysis of a short sequence from the highly conserved DNA polymerase (DNApol) gene previously identified two variants in a single host population. Here we compared genetic variance in blood samples from 66 known individuals of this same free-ranging badger population using a partial sequence comprising 2874 nucleotides of the DNApol gene, among which we identified 15 nucleotide differences resulting in 5 amino acid differences. Prevalence was 86% (59/66) for the common and 17% (11/66) for the novel variant, with 6% (4/66) of badgers presenting with coinfection. MusGHV-1 variants were distributed unevenly across the population, with individuals infected with the novel genotype clustered in 3 of 25 contiguous social groups. Individuals infected with the novel variant had significantly higher MusGHV-1 viral loads in their blood (p = 0.002) after adjusting for age (juveniles > adults, p < 0.001) and season (summer > spring and autumn, p = 0.005; mixed-effect linear regression), likely indicating higher virulence of the novel variant. Further genome-wide analyses of MusGHV-1 host resistance genes and host phenotypic variations are required to clarify the drivers and sequelae of this new MusGHV-1 variant.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Adverse weather during in utero development is linked to higher rates of later-life herpesvirus reactivation in adult European badgers, Meles meles

Maternal immune and/or metabolic conditions relating to stress or nutritional status can affect in utero development among offspring with subsequent implications for later-life responses to infections. We used free-ranging European badgers as a host-pathogen model to investigate how prenatal weather conditions affect later-life herpesvirus genital tract reactivation. We applied a sliding window analysis of weather conditions to 164 samples collected in 2018 from 95 individuals born between 2005-2016. We test if the monthly mean and variation in rainfall and temperature experienced by their mother during the 12 months of delayed implantation and gestation prior to parturition subsequently affected individual herpes reactivation rates among these offspring. We identified four influential prenatal seasonal weather windows that corresponded with previously identified critical climatic conditions affecting badger survival, fecundity and body condition. These all occurred during the pre-implantation rather than the post-implantation period. We conclude that environmental cues during the in utero period of delayed implantation may result in changes that affect an individual's developmental programming against infection or viral reactivation later in life. This illustrates how prenatal adversity caused by environmental factors, such as climate change, can impact wildlife health and population dynamics-an interaction largely overlooked in wildlife management and conservation programmes.

Stress-Related Herpesvirus Reactivation in Badgers Can Result in Clostridium Proliferation

Clostridium perfringens is an important food-borne zoonotic pathogen and a member of the commensal gut microbiome of many mammals. Predisposing factors such as coinfection with other pathogens or diet change can, however, cause overgrowth and subsequent disease development. Here we investigated the occurrence of C. perfringens in a free-ranging badger population with up to 100% prevalence of herpesvirus infection. Herpesvirus reactivation is known to be associated with increased susceptibility bacterial infections. PCR screening of rectal swabs from 69 free-ranging badgers revealed 15.9% (11/69, 95% CI = 9.1-26.3%) prevalence of detectable C. perfringens (Type A) DNA in the digestive tracts of assymptomatic animals. The results of Fisher's exact test revealed C. perfringens detection was not biased by age, sex and seasons. However, badgers with genital tract gammaherpesvirus (MusGHV-1) reactivation (p = 0.007) and infection with a specific MusGHV-1 genotype (p = 0.019) were more prone to of C. perfringens proliferation, indicating coinfection biased dynamics of intestinal C. perfringens. An inclusion pattern analysis further indicated that, causally, MusGHV-1 reactivation potentiated C. perfringens detection. Whether or not specific MusGHV-1 genotype infection or reactivation plays a role in C. perfringens overgrowth or disease development in badgers will require further investigation. Nevertheless, a postmortem examination of a single badger that died of fatal disease, likely associated with C. perfringens, revealed MusGHV-1 detection in the small intestine.

Hybridization with mountain hares increases the functional allelic repertoire in brown hares

Brown hares (Lepus europaeus Pallas) are able to hybridize with mountain hares (L. timidus Linnaeus) and produce fertile offspring, which results in cross-species gene flow. However, not much is known about the functional significance of this genetic introgression. Using targeted sequencing of candidate loci combined with mtDNA genotyping, we found the ancestral genetic diversity in the Finnish brown hare to be small, likely due to founder effect and range expansion, while gene flow from mountain hares constitutes an important source of functional genetic variability. Some of this variability, such as the alleles of the mountain hare thermogenin (uncoupling protein 1, UCP1), might have adaptive advantage for brown hares, whereas immunity-related MHC alleles are reciprocally exchanged and maintained via balancing selection. Our study offers a rare example where an expanding species can increase its allelic variability through hybridization with a congeneric native species, offering a route to shortcut evolutionary adaptation to the local environmental conditions.

Borrelia Infection in Bank Voles Myodes glareolus Is Associated With Specific DQB Haplotypes Which Affect Allelic Divergence Within Individuals

The high polymorphism of Major Histocompatibility Complex (MHC) genes is generally considered to be a result of pathogen-mediated balancing selection. Such selection may operate in the form of heterozygote advantage, and/or through specific MHC allele–pathogen interactions. Specific MHC allele–pathogen interactions may promote polymorphism via negative frequency-dependent selection (NFDS), or selection that varies in time and/or space because of variability in the composition of the pathogen community (fluctuating selection; FS). In addition, divergent allele advantage (DAA) may act on top of these forms of balancing selection, explaining the high sequence divergence between MHC alleles. DAA has primarily been thought of as an extension of heterozygote advantage. However, DAA could also work in concert with NFDS though this is yet to be tested explicitly. To evaluate the importance of DAA in pathogen-mediated balancing selection, we surveyed allelic polymorphism of MHC class II DQB genes in wild bank voles (Myodes glareolus) and tested for associations between DQB haplotypes and infection by Borrelia afzelii, a tick-transmitted bacterium causing Lyme disease in humans. We found two significant associations between DQB haplotypes and infection status: one haplotype was associated with lower risk of infection (resistance), while another was associated with higher risk of infection (susceptibility). Interestingly, allelic divergence within individuals was higher for voles with the resistance haplotype compared to other voles. In contrast, allelic divergence was lower for voles with the susceptibility haplotype than other voles. The pattern of higher allelic divergence in individuals with the resistance haplotype is consistent with NFDS favouring divergent alleles in a natural population, hence selection where DAA works in concert with NFDS.

Towards a more healthy conservation paradigm: integrating disease and molecular ecology to aid biological conservation†

Parasites, and the diseases they cause, are important from an ecological and evolutionary perspective because they can negatively affect host fitness and can regulate host populations. Consequently, conservation biology has long recognized the vital role that parasites can play in the process of species endangerment and recovery. However, we are only beginning to understand how deeply parasites are embedded in ecological systems, and there is a growing recognition of the important ways in which parasites affect ecosystem structure and function. Thus, there is an urgent need to revisit how parasites are viewed from a conservation perspective and broaden the role that disease ecology plays in conservation-related research and outcomes. This review broadly focusses on the role that disease ecology can play in biological conservation. Our review specifically emphasizes on how the integration of tools and analytical approaches associated with both disease and molecular ecology can be leveraged to aid conservation biology. Our review first concentrates on disease-mediated extinctions and wildlife epidemics. We then focus on elucidating how host–parasite interactions has improved our understanding of the eco-evolutionary dynamics affecting hosts at the individual, population, community and ecosystem scales. We believe that the role of parasites as drivers and indicators of ecosystem health is especially an exciting area of research that has the potential to fundamentally alter our view of parasites and their role in biological conservation. The review concludes with a broad overview of the current and potential applications of modern genomic tools in disease ecology to aid biological conservation.

Pathogen-mediated selection favours the maintenance of innate immunity gene polymorphism in a widespread wild ungulate

Toll-like receptors (TLR) play a central role in recognition and host frontline defence against a wide range of pathogens. A number of recent studies have shown that TLR genes (Tlrs) often exhibit large polymorphism in natural populations. Yet, there is little knowledge on how this polymorphism is maintained and how it influences disease susceptibility in the wild. In previous work, we showed that some Tlrs exhibit similarly high levels of genetic diversity as genes of the Major Histocompatibility Complex (MHC), and signatures of contemporary balancing selection in roe deer (Capreolus capreolus), the most abundant cervid species in Europe. Here, we investigated the evolutionary mechanisms by which pathogen-mediated selection could shape this innate immunity genetic diversity by examining the relationships between Tlr (Tlr2, Tlr4 and Tlr5) genotypes (heterozygosity status and presence of specific alleles) and infections with Toxoplasma and Chlamydia, two widespread intracellular pathogens known to cause reproductive failure in ungulates. We showed that Toxoplasma and Chlamydia exposures vary significantly across years and landscape features with few co-infection events detected and that the two pathogens exert antagonistic selection on Tlr2 polymorphism. By contrast, we found limited support for Tlr heterozygote advantage. Our study confirmed the importance of looking beyond Mhc genes in wildlife immunogenetic studies. It also emphasized the necessity to consider multiple pathogen challenges and their spatiotemporal variation to improve our understanding of vertebrate defence evolution against pathogens.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.

Patterns of Genital Tract Mustelid Gammaherpesvirus 1 (Musghv-1) Reactivation Are Linked to Stressors in European Badgers (Meles Meles)

Gammaherpesvirus reactivation can promote diseases or impair reproduction. Understanding reactivation patterns and associated risks of different stressors is therefore important. Nevertheless, outside the laboratory or captive environment, studies on the effects of stress on gammaherpesvirus reactivation in wild mammals are lacking. Here we used Mustelid gammaherpesvirus 1 (MusGHV-1) infection in European badgers (Meles meles) as a host-pathogen wildlife model to study the effects of a variety of demographic, physiological and environmental stressors on virus shedding in the genital tract. We collected 251 genital swabs from 150 free-ranging individuals across three seasons and screened them for the presence of MusGHV-1 DNA using PCR targeting the DNA polymerase gene. We explored possible links between MusGHV-1 DNA presence and seven variables reflecting stressors, using logistic regression analysis. The results reveal different sets of risk factors between juveniles and adults, likely reflecting primary infection and reactivation. In adults, virus shedding was more likely in badgers in poorer body condition and younger than 5 years or older than 7; while in juveniles, virus shedding is more likely in females and individuals in better body condition. However, living in social groups with more cubs was a risk factor for all badgers. We discuss possible explanations for these risk factors and their links to stress in badgers.

Patterns of MHC-dependent sexual selection in a free-living population of sheep

The MHC is one of the most polymorphic gene clusters in vertebrates and play an essential role in adaptive immunity. Apart from pathogen-mediated selection, sexual selection can also contribute to the maintenance of MHC diversity. MHC-dependent sexual selection could occur via several mechanisms but at present there is no consensus as to which of these mechanisms are involved and their importance. Previous studies have often suffered from limited genetic and behavioural data and small sample size, and were rarely able to examine all the mechanisms together, determine whether signatures of MHC-based non-random mating are independent of genomic effects or differentiate whether MHC-dependent sexual selection takes place at the pre- or post-copulatory stage. In this study, we use Monte Carlo simulation to investigate evidence for non-random MHC-dependent mating patterns by all three mechanisms in a free-living population of Soay sheep. Using 1710 sheep diplotyped at the MHC class IIa region and genome-wide SNPs, together with field observations of consorts, we found sexual selection against a particular haplotype in males at the pre-copulatory stage and sexual selection against female MHC heterozygosity during the rut. We also found MHC-dependent disassortative mating at the post-copulatory stage, along with strong evidence of inbreeding avoidance at both stages. However, results from generalized linear mixed models suggest that the pattern of MHC-dependent disassortative mating could be a by-product of inbreeding avoidance. Our results therefore suggest that while multiple apparent mechanisms of non-random mating with respect to the MHC may occur, some of them have alternative explanations.

Autochthonous Trypanosoma spp. in European Mammals: A Brief Journey amongst the Neglected Trypanosomes

The genus Trypanosoma includes flagellated protozoa belonging to the family Trypanosomatidae (Euglenozoa, Kinetoplastida) that can infect humans and several animal species. The most studied species are those causing severe human pathology, such as Chagas disease in South and Central America, and the human African trypanosomiasis (HAT), or infections highly affecting animal health, such as nagana in Africa and surra with a wider geographical distribution. The presence of these Trypanosoma species in Europe has been thus far linked only to travel/immigration history of the human patients or introduction of infected animals. On the contrary, little is known about the epidemiological status of trypanosomes endemically infecting mammals in Europe, such as Trypanosomatheileri in ruminants and Trypanosomalewisi in rodents and other sporadically reported species. This brief review provides an updated collection of scientific data on the presence of autochthonous Trypanosoma spp. in mammals on the European territory, in order to support epidemiological and diagnostic studies on Trypanosomatid parasites.

Spillover events of rabbit haemorrhagic disease virus 2 (recombinant GI.4P-GI.2) from Lagomorpha to Eurasian badger

Rabbit haemorrhagic disease (RHD) is a major threat to domestic and wild European rabbits. Presently, in Europe, the disease is caused mainly by Rabbit haemorrhagic disease virus 2 (RHDV2/b or Lagovirus europaeus GI.2), the origin of which is still unclear, as no RHDV2 reservoir hosts were identified. After the RHDV2 emergence in 2010, viral RNA was detected in a few rodent species. Furthermore, RHDV2 was found to cause disease in some hare species resembling the disease in rabbits, evidencing the ability of the virus to cross the species barrier. In this study, through molecular, histopathologic, antigenic and morphological evidences, we demonstrate the presence and replication of RHDV2 in Eurasian badgers (Meles meles) found dead in the district of Santarém, Portugal, between March 2017 and January 2020. In these animals, we further classify the RHDV2 as a Lagovirus europaeus recombinant GI.4P-GI.2. Our results indicate that Meles meles is susceptible to RHDV2, developing systemic infection, and excreting the virus in the faeces. Given the high viral loads seen in several organs and matrices, we believe that transmission to the wild rabbit is likely. Furthermore, transmission electron microscopy data show the presence of calicivirus compatible virions in the nucleus of hepatocytes, which constitutes a paradigm shift for caliciviruses' replication cycle.

Negative density-dependent parasitism in a group-living carnivore

Animals living at high population densities commonly experience greater exposure to disease, leading to increased parasite burdens. However, social animals can benefit immunologically and hygienically from cooperation, and individuals may alter their socio-spatial behaviour in response to infection, both of which could counteract density-related increases in exposure. Consequently, the costs and benefits of sociality for disease are often uncertain. Here, we use a long-term study of a wild European badger population (Meles meles) to investigate how within-population variation in host density determines infection with multiple parasites. Four out of five parasite taxa exhibited consistent spatial hotspots of infection, which peaked among badgers living in areas of low local population density. Combined movement, survival, spatial and social network analyses revealed that parasite avoidance was the likely cause of this negative density dependence, with possible roles for localized mortality, encounter-dilution effects, and micronutrient-enhanced immunity. These findings demonstrate that animals can organize their societies in space to minimize parasite infection, with important implications for badger behavioural ecology and for the control of badger-associated diseases.

Estimation of environmental, genetic and parental age at conception effects on telomere length in a wild mammal

Understanding individual variation in fitness-related traits requires separating the environmental and genetic determinants. Telomeres are protective caps at the ends of chromosomes that are thought to be a biomarker of senescence as their length predicts mortality risk and reflect the physiological consequences of environmental conditions. The relative contribution of genetic and environmental factors to individual variation in telomere length is, however, unclear, yet important for understanding its evolutionary dynamics. In particular, the evidence for transgenerational effects, in terms of parental age at conception, on telomere length is mixed. Here, we investigate the heritability of telomere length, using the 'animal model', and parental age at conception effects on offspring telomere length in a wild population of European badgers (Meles meles). Although we found no heritability of telomere length and low evolvability (<0.001), our power to detect heritability was low and a repeatability of 2% across individual lifetimes provides a low upper limit to ordinary narrow-sense heritability. However, year (32%) and cohort (3%) explained greater proportions of the phenotypic variance in telomere length, excluding qPCR plate and row variances. There was no support for cross-sectional or within-individual parental age at conception effects on offspring telomere length. Our results indicate a lack of transgenerational effects through parental age at conception and a low potential for evolutionary change in telomere length in this population. Instead, we provide evidence that individual variation in telomere length is largely driven by environmental variation in this wild mammal.

Effects of Mustelid gammaherpesvirus 1 (MusGHV-1) Reactivation in European Badger (Meles meles) Genital Tracts on Reproductive Fitness

Reactivation of latent Gammaherpesvirus in the genital tract can lead to reproductive failure in domestic animals. Nevertheless, this pathophysiology has not received formal study in wild mammals. High prevalence of Mustelid gammaherpesvirus 1 (MusGHV-1) DNA detected in the genital tracts of European badgers (Meles meles) implies that this common pathogen may be a sexual transmitted infection. Here we used PCR to test MusGHV-1 DNA prevalence in genital swabs collected from 144 wild badgers in Ireland (71 males, 73 females) to investigate impacts on male fertility indicators (sperm abundance and testes weight) and female fecundity (current reproductive output). MusGHV-1 reactivation had a negative effect on female reproduction, but not on male fertility; however males had a higher risk of MusGHV-1 reactivation than females, especially during the late-winter mating season, and genital MusGHV-1 reactivation differed between age classes, where 3-5 year old adults had significantly lower reactivation rates than younger or older ones. Negative results in foetal tissues from MusGHV-1 positive mothers indicated that cross-placental transmission was unlikely. This study has broader implications for how wide-spread gammaherpesvirus infections could affect reproductive performance in wild Carnivora species.

Neoplasms and novel gammaherpesviruses in critically endangered captive European minks (Mustela lutreola)

The European mink (Mustela lutreola) is a riparian mustelid, considered one of the most endangered carnivores in the world. Alpha, beta and gammaherpesviruses described in mustelids have been occasionally associated with different pathological processes. However, there is no information about the herpesviruses species infecting European minks. In this study, 141 samples of swabs (oral, conjunctival, anal), faeces and tissues from 23 animals were analysed for herpesvirus (HV) using a pan-HV-PCR assay. Two different, potentially novel, gammaherpesvirus species were identified in 12 samples from four animals (17.3%), and tentatively named Mustelid gammaherpesvirus-2 (MUGHV-2) and MuGHV-3. Gross examination was performed on dead minks (n = 11), while histopathology was performed using available samples from HV-positive individuals (n = 2), identifying several neoplasms, including B-cell lymphoma (identified by immunohistochemistry) with intralesional syncytia and intranuclear inclusion bodies characteristic of HV (n = 1), pulmonary adenocarcinoma (n = 1), and biliary (n = 1) and preputial (n = 1) cystadenomas, as well as other lesions (e.g., axonal vacuolar degeneration [n = 2] and neuritis [n = 1]). Viral particles, consistent with HVs, were observed by electron microscopy in the mink with neural lymphoma and inclusion bodies. This is the first description of neoplasms and concurrent gammaherpesvirus infection in European minks. The pathological, ultrastructural and PCR findings (MuGHV-2) in the European mink with lymphoma strongly suggest a potential role for this novel gammaherpesvirus in its pathogenesis, as it has been reported in other HV-infected species with lymphoma. The occurrence of neural lymphoma with intralesional syncytia and herpesviral inclusions is, however, unique among mammals. Further research is warranted to elucidate the potential oncogenic properties of gammaherpesviruses in European mink and their epidemiology in the wild population.

Social effects on age-related and sex-specific immune cell profiles in a wild mammal

Evidence for age-related changes in innate and adaptive immune responses is increasing in wild populations. Such changes have been linked to fitness, and knowledge of the factors driving immune response variation is important for understanding the evolution of immunity. Age-related changes in immune profiles may be owing to factors such as immune system development, sex-specific behaviour and responses to environmental conditions. Social environments may also contribute to variation in immunological responses, for example, through transmission of pathogens and stress arising from resource and mate competition. Yet, the impact of the social environment on age-related changes in immune cell profiles is currently understudied in the wild. Here, we tested the relationship between leukocyte cell composition (proportion of neutrophils and lymphocytes [innate and adaptive immunity, respectively] that were lymphocytes) and age, sex and group size in a wild population of European badgers (Meles meles). We found that the proportion of lymphocytes in early life was greater in males in smaller groups compared to larger groups, but with a faster age-related decline in smaller groups. By contrast, the proportion of lymphocytes in females was not significantly related to age or group size. Our results provide evidence of sex-specific age-related changes in immune cell profiles in a wild mammal, which are influenced by the social environment.

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

Genetic Diversity and Differentiation at Structurally Varying MHC Haplotypes and Microsatellites in Bottlenecked Populations of Endangered Crested Ibis

Investigating adaptive potential and understanding the relative roles of selection and genetic drift in populations of endangered species are essential in conservation. Major histocompatibility complex (MHC) genes characterized by spectacular polymorphism and fitness association have become valuable adaptive markers. Herein we investigate the variation of all MHC class I and II genes across seven populations of an endangered bird, the crested ibis, of which all current individuals are offspring of only two pairs. We inferred seven multilocus haplotypes from linked alleles in the Core Region and revealed structural variation of the class II region that probably evolved through unequal crossing over. Based on the low polymorphism, structural variation, strong linkage, and extensive shared alleles, we applied the MHC haplotypes in population analysis. The genetic variation and population structure at MHC haplotypes are generally concordant with those expected from microsatellites, underlining the predominant role of genetic drift in shaping MHC variation in the bottlenecked populations. Nonetheless, some populations showed elevated differentiation at MHC, probably due to limited gene flow. The seven populations were significantly differentiated into three groups and some groups exhibited genetic monomorphism, which can be attributed to founder effects. We therefore propose various strategies for future conservation and management.

Both candidate gene and neutral genetic diversity correlate with parasite resistance in female Mediterranean mouflon

BACKGROUND

Parasite infections can have substantial impacts on population dynamics and are accordingly a key challenge for wild population management. Here we studied genetic mechanisms driving parasite resistance in a large herbivore through a comprehensive approach combining measurements of neutral (16 microsatellites) and adaptive (MHC DRB1 exon 2) genetic diversity and two types of gastrointestinal parasites (nematodes and coccidia).

RESULTS

While accounting for other extrinsic and intrinsic predictors known to impact parasite load, we show that both neutral genetic diversity and DRB1 are associated with resistance to gastrointestinal nematodes. Intermediate levels of multi-locus heterozygosity maximized nematodes resistance, suggesting that both in- and outbreeding depression might occur in the population. DRB1 heterozygosity and specific alleles effects were detected, suggesting the occurrence of heterozygote advantage, rare-allele effects and/or fluctuating selection. On the contrary, no association was detected between genetic diversity and resistance to coccidia, indicating that different parasite classes are impacted by different genetic drivers.

CONCLUSIONS

This study provides important insights for large herbivores and wild sheep pathogen management, and in particular suggests that factors likely to impact genetic diversity and allelic frequencies, including global changes, are also expected to impact parasite resistance.

Evolution of MHC class I genes in Eurasian badgers, genus Meles (Carnivora, Mustelidae)

Because of their role in immune defense against pathogens, major histocompatibility complex (MHC) genes are useful in evolutionary studies on how wild vertebrates adapt to their environments. We investigated the molecular evolution of MHC class I (MHCI) genes in four closely related species of Eurasian badgers, genus Meles. All four species of badgers showed similarly high variation in MHCI sequences compared to other Carnivora. We identified 7-21 putatively functional MHCI sequences in each of the badger species, and 2-7 sequences per individual, indicating the existence of 1-4 loci. MHCI exon 2 and 3 sequences encoding domains α1 and α2 exhibited different clade topologies in phylogenetic networks. Non-synonymous nucleotide substitutions at codons for antigen-binding sites exceeded synonymous substitutions for domain α1 but not for domain α2, suggesting that the domains α1 and α2 likely had different evolutionary histories in these species. Positive selection and recombination seem to have shaped the variation in domain α2, whereas positive selection was dominant in shaping the variation in domain α1. In the separate phylogenetic analyses for exon 2, exon 3, and intron 2, each showed three clades of Meles alleles, with rampant trans-species polymorphism, indicative of the long-term maintenance of ancestral MHCI polymorphism by balancing selection.

MHC-mediated sexual selection on birdsong: Generic polymorphism, particular alleles and acoustic signals

Several hypotheses predict that the major histocompatibility complex (MHC) drives mating preference in females. Olfactory, colour or morphological traits are often found as reliable signals of the MHC profile, but the role of avian song mediating MHC-based female choice remains largely unexplored. We investigated the relationship between several MHC and acoustic features in the collared flycatcher (Ficedula albicollis), a European passerine with complex songs. We screened a fragment of the class IIB second exon of the MHC molecule, of which individuals harbour 4-15 alleles, while considerable sequence diversity is maintained at the population level. To make statistical inferences from a large number of comparisons, we adopted both null-hypothesis testing and effect size framework in combination with randomization procedures. After controlling for potential confounding factors, neither MHC allelic diversity nor the presence of particular alleles was associated remarkably with the investigated qualitative and quantitative song traits. Furthermore, genetic similarity among males based on MHC sequences was not reflected by the similarity in their song based on syllable content. Overall, these results suggest that the relationship between features of song and the allelic composition and diversity of MHC is not strong in the studied species. However, a biologically motivated analysis revealed that individuals that harbour an MHC allele that impairs survival perform songs with broader frequency range. This finding suggests that certain aspects of the song may bear reliable information concerning the MHC profile of the individuals, which can be used by females to optimize mate choice.

Experimental manipulation of population-level MHC diversity controls pathogen virulence evolution in Mus musculus

The virulence levels attained by serial passage of pathogens through similar host genotypes are much higher than observed in natural systems; however, it is unknown what keeps natural virulence levels below these empirically demonstrated maximum levels. One hypothesis suggests that host diversity impedes pathogen virulence, because adaptation to one host genotype carries trade-offs in the ability to replicate and cause disease in other host genotypes. To test this hypothesis, with the simplest level of population diversity within the loci of the major histocompatibility complex (MHC), we serially passaged Friend virus complex (FVC) through two rounds, in hosts with either the same MHC genotypes (pure passage) or hosts with different MHC genotypes (alternated passage). Alternated passages showed a significant overall reduction in viral titre (31%) and virulence (54%) when compared to pure passages. Furthermore, a resistant host genotype initially dominated any effects due to MHC diversity; however, when FVC was allowed to adapt to the resistant host genotype, predicted MHC effects emerged; that is, alternated lines show reduced virulence. These data indicate serial exposure to diverse MHC genotypes is an impediment to pathogen adaptation, suggesting genetic variation at MHC loci is important for limiting virulence in a rapidly evolving pathogen and supports negative frequency-dependent selection as a force maintaining MHC diversity in host populations.

GENITAL TRACT SCREENING FINDS WIDESPREAD INFECTION WITH MUSTELID GAMMAHERPESVIRUS 1 IN THE EUROPEAN BADGER ( MELES MELES)

:  Sexually transmitted diseases (STDs) can be important drivers of population dynamics because of their negative effects on reproduction. However, screening for STDs, especially in wildlife populations, is widely neglected. Using the promiscuous, polygynandrous European badger ( Meles meles) as a model, we investigated the presence and prevalence of herpesviruses (HVs) in a wild, high-density population and assessed potential differences in somatic fitness and female reproductive condition between infected and uninfected individuals. We collected n=98 genital swabs from 71 females (51 adults and 20 cubs) and 27 males (26 adults and 1 cub) during spring and summer 2015. Using a PCR specific for a mustelid α-HV, all genital-swab samples tested negative. In a panherpes PCR, a γ-HV was found in 55% (54/98; 39 adults and 15 cubs), identified as mustelid gammaherpesvirus 1 (MusGHV-1) using DNA sequencing. This contrasts with the results of a previous study, which reported MusGHV-1 in 98% (354/361) of blood samples taken from 218 badgers in the same population using PCR. The detection of MusHV-1 in the female reproductive tract strongly indicates the potential for a horizontal and, likely also a vertical, route of transmission. Our results suggest a potential linkage of genital HVs and impaired future reproductive success in females, but because reproductive failure can have many reasons in badgers, the causative link of this negative relationship remains to be investigated.

No Compensatory Relationship between the Innate and Adaptive Immune System in Wild-Living European Badgers

The innate immune system provides the primary vertebrate defence system against pathogen invasion, but it is energetically costly and can have immune pathological effects. A previous study in sticklebacks found that intermediate major histocompatibility complex (MHC) diversity correlated with a lower leukocyte coping capacity (LCC), compared to individuals with fewer, or many, MHC alleles. The organization of the MHC genes in mammals, however, differs to the highly duplicated MHC genes in sticklebacks by having far fewer loci. Using European badgers (Meles meles), we therefore investigated whether innate immune activity, estimated functionally as the ability of an individual’s leukocytes to produce a respiratory burst, was influenced by MHC diversity. We also investigated whether LCC was influenced by factors such as age-class, sex, body condition, season, year, neutrophil and lymphocyte counts, and intensity of infection with five different pathogens. We found that LCC was not associated with specific MHC haplotypes, MHC alleles, or MHC diversity, indicating that the innate immune system did not compensate for the adaptive immune system even when there were susceptible MHC alleles/haplotypes, or when the MHC diversity was low. We also identified a seasonal and annual variation of LCC. This temporal variation of innate immunity was potentially due to physiological trade-offs or temporal variation in pathogen infections. The innate immunity, estimated as LCC, does not compensate for MHC diversity suggests that the immune system may function differently between vertebrates with different MHC organizations, with implications for the evolution of immune systems in different taxa.

Major histocompatibility complex class II DAB alleles associated with intestinal parasite load in the vulnerable Chinese egret (Egretta eulophotes)

The maintenance of major histocompatibility complex (MHC) polymorphism has been hypothesized to result from many mechanisms such as rare-allele advantage, heterozygote advantage, and allele counting. In the study reported herein, 224 vulnerable Chinese egrets (Egretta eulophotes) were used to examine these hypotheses as empirical results derived from bird studies are rare. Parasite survey showed that 147 (65.63%) individuals were infected with 1-3 helminths, and 82.31% of these infected individuals carried Ascaridia sp. Using asymmetric polymerase chain reaction technique, 10 DAB1, twelve DAB2, and three DAB3 exon 2 alleles were identified at each single locus. A significant association of the rare allele Egeu-DAB2*05 (allele frequency: 0.022) with helminth resistance was found for all helminths, as well as for the most abundant morphotype Ascaridia sp. in the separate analyses. Egeu-DAB2*05 occurred frequently in uninfected individuals, and individuals carrying Egeu-DAB2*05 had significantly lower helminth morphotypes per individual (HMI) (the number of HMI) and the fecal egg count values. Further, the parasite infection measurements were consistently lower in individuals with an intermediate number of different alleles in the duplicated DAB loci. Significantly, heterozygosity within each DAB locus was not correlated with any parasite infection measurements. These results indicate that the diversity in MHC Egeu-DAB gene is associated with intestinal parasite load and maintained by pathogen-driven selection that probably operate through both the rare-allele advantage and the allele counting strategy, and suggest that Egeu-DAB2*05 might be a valuable indicator of better resistance to helminth diseases in the vulnerable Chinese egret.

Mate choice for neutral and MHC genetic characteristics in Alpine marmots: different targets in different contexts?

Sexual selection through female mate choice for genetic characteristics has been suggested to be an important evolutionary force maintaining genetic variation in animal populations. However, the genetic targets of female mate choice are not clearly identified and whether female mate choice is based on neutral genetic characteristics or on particular functional loci remains an open question. Here, we investigated the genetic targets of female mate choice in Alpine marmots (Marmota marmota), a socially monogamous mammal where extra‐pair paternity (EPP) occurs. We used 16 microsatellites to describe neutral genetic characteristics and two MHC loci belonging to MHC class I and II as functional genetic characteristics. Our results reveal that (1) neutral and MHC genetic characteristics convey different information in this species, (2) social pairs show a higher MHC class II dissimilarity than expected under random mate choice, and (3) the occurrence of EPP increases when social pairs present a high neutral genetic similarity or dissimilarity but also when they present low MHC class II dissimilarity. Thus, female mate choice is based on both neutral and MHC genetic characteristics, and the genetic characteristics targeted seem to be context dependent (i.e., the genes involved in social mate choice and genetic mate choice differ). We emphasize the need for empirical studies of mate choice in the wild using both neutral and MHC genetic characteristics because whether neutral and functional genetic characteristics convey similar information is not universal.