Major histocompatibility complex (MHC) heterozygote superiority to natural multi-parasite infections in the water vole (Arvicola terrestris)

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.

Gene duplication and divergence produce divergent MHC genotypes without disassortative mating

Genes of the major histocompatibility complex (MHC) exhibit heterozygote advantage in immune defence, which in turn can select for MHC-disassortative mate choice. However, many species lack this expected pattern of MHC-disassortative mating. A possible explanation lies in evolutionary processes following gene duplication: if two duplicated MHC genes become functionally diverged from each other, offspring will inherit diverse multilocus genotypes even under random mating. We used locus-specific primers for high-throughput sequencing of two expressed MHC Class II B genes in Leach's storm-petrels, Oceanodroma leucorhoa, and found that exon 2 alleles fall into two gene-specific monophyletic clades. We tested for disassortative vs. random mating at these two functionally diverged Class II B genes, using multiple metrics and different subsets of exon 2 sequence data. With good statistical power, we consistently found random assortment of mates at MHC. Despite random mating, birds had MHC genotypes with functionally diverged alleles, averaging 13 amino acid differences in pairwise comparisons of exon 2 alleles within individuals. To test whether this high MHC diversity in individuals is driven by evolutionary divergence of the two duplicated genes, we built a phylogenetic permutation model. The model showed that genotypic diversity was strongly impacted by sequence divergence between the most common allele of each gene, with a smaller additional impact of monophyly of the two genes. Divergence of allele sequences between genes may have reduced the benefits of actively seeking MHC-dissimilar mates, in which case the evolutionary history of duplicated genes is shaping the adaptive landscape of sexual selection.

Development of novel associations between MHC alleles and susceptibility to parasitic infections in an isolated population of an endangered mammal

The role of pathogens in dynamics of endangered species is not fully understood, and the effect of infection often interacts with other processes affecting those species, such as fragmentation and isolation or loss of genetic variation. Small, isolated populations are prone to losing functional alleles due to demographic processes and genetic drift, which may diminish their ability to resist infection if immune genes are affected. Demographic processes may also alter the selective pressure exerted by a parasite, as they influence the rate of parasite transmission between individuals. In the present paper we studied changes in parasite infection levels and genetic variability in an isolated population of spotted suslik (Spermophillus suslicus). Over a three-year period (approx. three generations), when the population size remained relatively stable, we observed a considerable increase in parasite prevalence and infection intensity, followed by the development of novel associations between MHC DRB alleles and parasite burden. Contrary to expectations, the change in MHC allele frequency over time was not consistent with the effect of the allele - for instance, Spsu-DRB*07, associated with higher intensity of infection with a nematode Capillaria sp., increased in frequency from 11.8 to 20.2%. Yet, we found no signatures of selection in the studied loci. Our results show that an isolated, stable population may experience a sudden increase in parasitic infections, resulting in a development of novel associations between MHC alleles and parasite susceptibility/resistance, even though no signatures of selection can be found.

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.

Evidence for selection maintaining MHC diversity in a rodent species despite strong density fluctuations

Strong spatiotemporal variation in population size often leads to reduced genetic diversity limiting the adaptive potential of individual populations. Key genes of adaptive variation are encoded by the immune genes of the major histocompatibility complex (MHC) playing an essential role in parasite resistance. How MHC variation persists in rodent populations that regularly experience population bottlenecks remains an important topic in evolutionary genetics. We analysed the consequences of strong population fluctuations on MHC class II DRB exon 2 diversity in two distant common vole (Microtus arvalis) populations in three consecutive years using a high-throughput sequencing approach. In 143 individuals, we detected 25 nucleotide alleles translating into 14 unique amino acid MHC alleles belonging to at least three loci. Thus, the overall allelic diversity and amino acid distance among the remaining MHC alleles, used as a surrogate for the range of pathogenic antigens that can be presented to T-cells, are still remarkably high. Both study populations did not show significant population differentiation between years, but significant differences were found between sites. We concluded that selection processes seem to be strong enough to maintain moderate levels of MHC diversity in our study populations outcompeting genetic drift, as the same MHC alleles were conserved between years. Differences in allele frequencies between populations might be the outcome of different local parasite pressures and/or genetic drift. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in cyclic populations.

MHC, parasites and antler development in red deer: no support for the Hamilton & Zuk hypothesis

The Hamilton-Zuk hypothesis proposes that the genetic benefits of preferences for elaborated secondary sexual traits have their origins in the arms race between hosts and parasites, which maintains genetic variance in parasite resistance. Infection, in turn, can be reflected in the expression of costly sexual ornaments. However, the link between immune genes, infection and the expression of secondary sexual traits has rarely been investigated. Here, we explored whether the presence and identity of functional variants (supertypes) of the highly polymorphic major histocompatibility complex (MHC), which is responsible for the recognition of parasites, predict the load of lung and gut parasites and antler development in the red deer (Cervus elaphus). While we found MHC supertypes to be associated with infection by a number of parasite species, including debilitating lung nematodes, we did not find support for the Hamilton-Zuk hypothesis. On the contrary, we found that lung nematode load was positively associated with antler development. We also found that the supertypes that were associated with resistance to certain parasites at the same time cause susceptibility to others. Such trade-offs may undermine the potential genetic benefits of mate choice for resistant partners.

Does the parasite-mediated selection drive the MHC class IIB diversity in wild populations of European chub (Squalius cephalus)?

The genes of major histocompatibility complex (MHC) provide an excellent opportunity to study host-parasite relationships because they are expected to evolve in response to parasites and variation in parasite communities. In this study, we investigated the potential role of parasite-mediated selection acting on MHC class IIB (DAB) genes in European chub (Squalius cephalus) natural populations. We found significant differences between populations in metazoan parasites, neutral and adaptive genetic diversities. The analyses based on pairwise data revealed that populations with dissimilar MHC allelic profiles were geographically distant populations with significantly different diversity in microsatellites and a dissimilar composition of parasite communities. The results from the generalized estimating equations method (GEE) on the level of individuals revealed that metazoan parasite load in European chub was influenced by the diversity of DAB alleles as well as by the diversity of neutral genetic markers and host traits reflecting condition and immunocompetence. The multivariate co-inertia analysis showed specific associations between DAB alleles and parasite species. DAB1-like alleles were more involved in associations with ectoparasites, while DAB3-like alleles were positively associated with endoparasites which could suggest potential differences between DAB genes caused by different selection pressure. Our study revealed that parasite-mediated selection is not the only variable affecting MHC diversity in European chub; however, we strongly support the role of neutral processes as the main driver of DAB diversity across populations. In addition, our study contributes to the understanding of the evolution of MHC genes in wild living fish.

Tolerance to noninherited maternal antigens, reproductive microchimerism and regulatory T cell memory: 60 years after 'Evidence for actively acquired tolerance to Rh antigens'

Compulsory exposure to genetically foreign maternal tissue imprints in offspring sustained tolerance to noninherited maternal antigens (NIMA). Immunological tolerance to NIMA was first described by Dr. Ray D. Owen for women genetically negative for erythrocyte rhesus (Rh) antigen with reduced sensitization from developmental Rh exposure by their mothers. Extending this analysis to HLA haplotypes has uncovered the exciting potential for therapeutically exploiting NIMA-specific tolerance naturally engrained in mammalian reproduction for improved clinical outcomes after allogeneic transplantation. Herein, we summarize emerging scientific concepts stemming from tolerance to NIMA that includes postnatal maintenance of microchimeric maternal origin cells in offspring, expanded accumulation of immune suppressive regulatory T cells with NIMA-specificity, along with teleological benefits and immunological consequences of NIMA-specific tolerance conserved across mammalian species.

Lack of Spatial Immunogenetic Structure among Wolverine (Gulo gulo) Populations Suggestive of Broad Scale Balancing Selection

Elucidating the adaptive genetic potential of wildlife populations to environmental selective pressures is fundamental for species conservation. Genes of the major histocompatibility complex (MHC) are highly polymorphic, and play a key role in the adaptive immune response against pathogens. MHC polymorphism has been linked to balancing selection or heterogeneous selection promoting local adaptation. However, spatial patterns of MHC polymorphism are also influenced by gene flow and drift. Wolverines are highly vagile, inhabiting varied ecoregions that include boreal forest, taiga, tundra, and high alpine ecosystems. Here, we investigated the immunogenetic variation of wolverines in Canada as a surrogate for identifying local adaptation by contrasting the genetic structure at MHC relative to the structure at 11 neutral microsatellites to account for gene flow and drift. Evidence of historical positive selection was detected at MHC using maximum likelihood codon-based methods. Bayesian and multivariate cluster analyses revealed weaker population genetic differentiation at MHC relative to the increasing microsatellite genetic structure towards the eastern wolverine distribution. Mantel correlations of MHC against geographical distances showed no pattern of isolation by distance (IBD: r = -0.03, p = 0.9), whereas for microsatellites we found a relatively strong and significant IBD (r = 0.54, p = 0.01). Moreover, we found a significant correlation between microsatellite allelic richness and the mean number of MHC alleles, but we did not observe low MHC diversity in small populations. Overall these results suggest that MHC polymorphism has been influenced primarily by balancing selection and to a lesser extent by neutral processes such as genetic drift, with no clear evidence for local adaptation. This study contributes to our understanding of how vulnerable populations of wolverines may respond to selective pressures across their range.

Genetic structure in insular and mainland populations of house sparrows (Passer domesticus) and their hemosporidian parasites

Small and isolated populations usually exhibit low levels of genetic variability, and thus, they are expected to have a lower capacity to adapt to changes in environmental conditions, such as exposure to pathogens and parasites. Comparing the genetic variability of selectively neutral versus functional loci allows one to assess the evolutionary history of populations and their future evolutionary potential. The genes of the major histocompatibility complex (MHC) control immune recognition of parasites, and their unusually high diversity is genes which is likely driven by parasite-mediated balancing selection. Here, we examined diversity and differentiation of neutral microsatellite loci and functional MHC class I genes in house sparrows (Passer domesticus), living in six insular and six mainland populations, and we aimed to determine whether their diversity or differentiation correlates with the diversity and the prevalence of infection of hemosporidian parasites. We found that island bird populations tended to have lower neutral genetic variability, whereas MHC variability gene was similar between island and mainland populations. Similarly, island populations tended to show greater genetic differentiation than mainland populations, especially at microsatellite markers. The maintenance of MHC genetic diversity and its less marked structure in the island populations could be attributed to balancing-selection. The greater MHC differentiation among populations was negatively correlated with similarity in blood parasites (prevalence and diversity of parasite strains) between populations. Even at low prevalence and small geographical scale, haemosporidian parasites might contribute to structure the variability of immune genes among populations of hosts.

Heterozygote advantage at MHC DRB may influence response to infectious disease epizootics

The effect of MHC polymorphism on individual fitness variation in the wild remains equivocal; however, much evidence suggests that heterozygote advantage is a major determinant. To understand the contribution of MHC polymorphism to individual disease resistance or susceptibility in natural populations, we investigated two MHC class II B loci, DQB and DRB, in the New Zealand sea lion (NZSL, Phocarctos hookeri). The NZSL is a threatened species which is unusually susceptible to death by bacterial infection at an early age; it has suffered three bacterial induced epizootics resulting in high mortality levels of young pups since 1997. The MHC DQB and DRB haplotypes of dead NZSL pups with known cause of death (bacteria, enteritis or trauma) were sequenced and reconstructed, compared to pups that survived beyond 2 months of age, and distinct MHC DRB allele frequency and genotype differences were identified. Two findings were striking: (i) one DRB allele was present only in dead pups, and (ii) one heterozygous DRB genotype, common in live pups, was absent from dead pups. These results are consistent with some functional relationship with these variants and suggest heterozygote advantage is operating at DRB. We found no association between heterozygosity and fitness at 17 microsatellite loci, indicating that general heterozygosity is not responsible for the effect on fitness detected here. This result may be a consequence of recurrent selection by multiple pathogen assault over recent years and highlights the importance of heterozygote advantage at MHC as a potential mechanism for fitness differences in wild populations.

Isolation over 35 years in a heated biotest basin causes selection on MHC class IIß genes in the European perch (Perca fluviatilis L.)

Genes that play key roles in host immunity such as the major histocompatibility complex (MHC) in vertebrates are expected to be major targets of selection. It is well known that environmental conditions can have an effect on host–parasite interactions and may thus influence the selection on MHC. We analyzed MHC class IIß variability over 35 years in a population of perch (Perca fluviatilis) from the Baltic Sea that was split into two populations separated from each other. One population was subjected to heating from cooling water of a nuclear power plant and was isolated from the surrounding environment in an artificial lake, while the other population was not subjected to any change in water temperature (control). The isolated population experienced a change of the allelic composition and a decrease in allelic richness of MHC genes compared to the control population. The two most common MHC alleles showed cyclic patterns indicating ongoing parasite–host coevolution in both populations, but the alleles that showed a cyclic behavior differed between the two populations. No such patterns were observed at alleles from nine microsatellite loci, and no genetic differentiation was found between populations. We found no indications for a genetic bottleneck in the isolated population during the 35 years. Additionally, differences in parasitism of the current perch populations suggest that a change of the parasite communities has occurred over the isolation period, although the evidence in form of in-depth knowledge of the change of the parasite community over time is lacking. Our results are consistent with the hypothesis of a selective sweep imposed by a change in the parasite community.

Host-parasite network structure is associated with community-level immunogenetic diversity

Genes of the major histocompatibility complex (MHC) encode proteins that recognize foreign antigens and are thus crucial for immune response. In a population of a single host species, parasite-mediated selection drives MHC allelic diversity. However, in a community-wide context, species interactions may modulate selection regimes because the prevalence of a given parasite in a given host may depend on its prevalence in other hosts. By combining network analysis with immunogenetics, we show that host species infected by similar parasites harbour similar alleles with similar frequencies. We further show, using a Bayesian approach, that the probability of mutual occurrence of a functional allele and a parasite in a given host individual is nonrandom and depends on other host-parasite interactions, driving co-evolution within subgroups of parasite species and functional alleles. Therefore, indirect effects among hosts and parasites can shape host MHC diversity, scaling it from the population to the community level.

Major histocompatibility complex alleles associated with parasite susceptibility in wild giant pandas

Major histocompatibility complex (MHC) polymorphism is thought to be driven by antagonistic coevolution between pathogens and hosts, mediated through either overdominance or frequency-dependent selection. However, investigations under natural conditions are still rare for endangered mammals which often exhibit depleted variation, and the mechanism of selection underlying the maintenance of characteristics remains a considerable debate. In this study, 87 wild giant pandas were used to investigate MHC variation associated with parasite load. With the knowledge of the MHC profile provided by the genomic data of the giant panda, seven DRB1, seven DQA1 and eight DQA2 alleles were identified at each single locus. Positive selection evidenced by a significantly higher number of non-synonymous substitutions per non-synonymous codon site relative to synonymous substitutions per synonymous codon site could only be detected at the DRB1 locus, which leads to the speculation that DRB1 may have a more important role in dealing with parasite infection for pandas. Coprological analyses revealed that 55.17% of individuals exhibited infection with 1-2 helminthes and 95.3% of infected pandas carried Baylisascaris shroederi. Using a generalized linear model, we found that Aime-DRB1*10 was significantly associated with parasite infection, but no resistant alleles could be detected. MHC heterozygosity of the pandas was found to be uncorrelated with the infection status or the infection intensity. These results suggested that the possible selection mechanisms in extant wild pandas may be frequency dependent rather than being determined by overdominance selection. Our findings could guide the candidate selection for the ongoing reintroduction or translocation of pandas.

Introgression from domestic goat generated variation at the major histocompatibility complex of Alpine ibex

The major histocompatibility complex (MHC) is a crucial component of the vertebrate immune system and shows extremely high levels of genetic polymorphism. The extraordinary genetic variation is thought to be ancient polymorphisms maintained by balancing selection. However, introgression from related species was recently proposed as an additional mechanism. Here we provide evidence for introgression at the MHC in Alpine ibex (Capra ibex ibex). At a usually very polymorphic MHC exon involved in pathogen recognition (DRB exon 2), Alpine ibex carried only two alleles. We found that one of these DRB alleles is identical to a DRB allele of domestic goats (Capra aegagrus hircus). We sequenced 2489 bp of the coding and non-coding regions of the DRB gene and found that Alpine ibex homozygous for the goat-type DRB exon 2 allele showed nearly identical sequences (99.8%) to a breed of domestic goats. Using Sanger and RAD sequencing, microsatellite and SNP chip data, we show that the chromosomal region containing the goat-type DRB allele has a signature of recent introgression in Alpine ibex. A region of approximately 750 kb including the DRB locus showed high rates of heterozygosity in individuals carrying one copy of the goat-type DRB allele. These individuals shared SNP alleles both with domestic goats and other Alpine ibex. In a survey of four Alpine ibex populations, we found that the region surrounding the DRB allele shows strong linkage disequilibria, strong sequence clustering and low diversity among haplotypes carrying the goat-type allele. Introgression at the MHC is likely adaptive and introgression critically increased MHC DRB diversity in the genetically impoverished Alpine ibex. Our finding contradicts the long-standing view that genetic variability at the MHC is solely a consequence of ancient trans-species polymorphism. Introgression is likely an underappreciated source of genetic diversity at the MHC and other loci under balancing selection.

The major histocompatibility complex (MHC), a crucial component of the defense against pathogens, contains the most polymorphic functional genes in vertebrate genomes. The extraordinary genetic variation is generally considered to be ancient. We investigated whether a previously neglected mechanism, introgression from related species, provides an additional source of MHC variation. We show that introgression from domestic goat dramatically increased genetic variation at the MHC of Alpine ibex, a species that had nearly gone extinct during the 18^th century, but has been restored to large numbers since. We show that Alpine ibex share one of only two alleles at a generally highly polymorphic MHC locus with domestic goats and that the chromosomal region containing the goat-type allele has a signature of recent introgression. Our finding contradicts the long-standing view that ancient trans-species polymorphism is the sole source of the extraordinary genetic variability at the MHC. Instead, we show that in Alpine ibex introgression generated genetic diversity at a MHC locus. Our study supports the view that loci favoring genetic polymorphism may be susceptible to adaptive introgression from related species and will encourage future research to identify unexpected signatures of introgression.

Nonsynonymous substitution rate heterogeneity in the peptide-binding region among different HLA-DRB1 lineages in humans

An extraordinary diversity of amino acid sequences in the peptide-binding region (PBR) of human leukocyte antigen [HLA; human major histocompatibility complex (MHC)] molecules has been maintained by balancing selection. The process of accumulation of amino acid diversity in the PBR for six HLA genes (HLA-A, B, C, DRB1, DQB1, and DPB1) shows that the number of amino acid substitutions in the PBR among alleles does not linearly correlate with the divergence time of alleles at the six HLA loci. At these loci, some pairs of alleles show significantly less nonsynonymous substitutions at the PBR than expected from the divergence time. The same phenomenon was observed not only in the HLA but also in the rat MHC. To identify the cause for this, DRB1 sequences, a representative case of a typical nonlinear pattern of substitutions, were examined. When the amino acid substitutions in the PBR were placed with maximum parsimony on a maximum likelihood tree based on the non-PBR substitutions, heterogeneous rates of nonsynonymous substitutions in the PBR were observed on several branches. A computer simulation supported the hypothesis that allelic pairs with low PBR substitution rates were responsible for the stagnation of accumulation of PBR nonsynonymous substitutions. From these observations, we conclude that the nonsynonymous substitution rate at the PBR sites is not constant among the allelic lineages. The deceleration of the rate may be caused by the coexistence of certain pathogens for a substantially long time during HLA evolution.

Parasite-mediated selection drives an immunogenetic trade-off in plains zebras (Equus quagga)

Pathogen evasion of the host immune system is a key force driving extreme polymorphism in genes of the major histocompatibility complex (MHC). Although this gene family is well characterized in structure and function, there is still much debate surrounding the mechanisms by which MHC diversity is selectively maintained. Many studies have investigated relationships between MHC variation and specific pathogens, and have found mixed support for and against the hypotheses of heterozygote advantage, frequency-dependent or fluctuating selection. Few, however, have focused on the selective effects of multiple parasite types on host immunogenetic patterns. Here, we examined relationships between variation in the equine MHC gene, ELA-DRA, and both gastrointestinal (GI) and ectoparasitism in plains zebras (Equus quagga). Specific alleles present at opposing population frequencies had antagonistic effects, with rare alleles associated with increased GI parasitism and common alleles with increased tick burdens. These results support a frequency-dependent mechanism, but are also consistent with fluctuating selection. Maladaptive GI parasite 'susceptibility alleles' were reduced in frequency, suggesting that these parasites may play a greater selective role at this locus. Heterozygote advantage, in terms of allele mutational divergence, also predicted decreased GI parasite burden in genotypes with a common allele. We conclude that an immunogenetic trade-off affects resistance/susceptibility to parasites in this system. Because GI and ectoparasites do not directly interact within hosts, our results uniquely show that antagonistic parasite interactions can be indirectly modulated through the host immune system. This study highlights the importance of investigating the role of multiple parasites in shaping patterns of host immunogenetic variation.

Extraordinary MHC class II B diversity in a non-passerine, wild bird: the Eurasian Coot Fulica atra (Aves: Rallidae)

The major histocompatibility complex (MHC) hosts the most polymorphic genes ever described in vertebrates. The MHC triggers the adaptive branch of the immune response, and its extraordinary variability is considered an evolutionary consequence of pathogen pressure. The last few years have witnessed the characterization of the MHC multigene family in a large diversity of bird species, unraveling important differences in its polymorphism, complexity, and evolution. Here, we characterize the first MHC class II B sequences isolated from a Rallidae species, the Eurasian Coot Fulica atra. A next-generation sequencing approach revealed up to 265 alleles that translated into 251 different amino acid sequences (β chain, exon 2) in 902 individuals. Bayesian inference identified up to 19 codons within the presumptive peptide-binding region showing pervasive evidence of positive, diversifying selection. Our analyses also detected a significant excess of high-frequency segregating sites (average Tajima's D = 2.36, P < 0.05), indicative of balancing selection. We found one to six different alleles per individual, consistent with the occurrence of at least three MHC class II B gene duplicates. However, the genotypes comprised of three alleles were by far the most abundant in the population investigated (49.4%), followed by those with two (29.6%) and four (17.5%) alleles. We suggest that these proportions are in agreement with the segregation of MHC haplotypes differing in gene copy number. The most widespread segregating haplotypes, according to our findings, would contain one single gene or two genes. The MHC class II of the Eurasian Coot is a valuable system to investigate the evolutionary implications of gene copy variation and extensive variability, the greatest ever found, to the best of our knowledge, in a wild population of a non-passerine bird.

Growth and fluctuating asymmetry of human newborns: influence of inbreeding and parental education

Historically, medical concerns about the deleterious effects of closely inbred marriages have focused on the risk posed by recessive Mendelian disease, with much less attention to developmental instability. We studied the effects of inbreeding (first-cousin marriage) on growth and fluctuating asymmetry of 200 full-term infants (101 inbred and 99 outbred) whose parents were of similar socioeconomic status in Sivas Province, Turkey. In addition to differences in their mean inbreeding coefficients (f = 1/16 for first cousins and f < 1/1,024 for unrelated parents), the consanguineous parents were less well educated (3 years, on average for both husbands and wives). We measured weight, height, head circumference, and chest circumference of the newborns, as well as four bilateral traits (ear width, ear length, and second and fourth digit lengths). After taking education into account, none of the measures of size (weight, height, head circumference, and chest circumference) and fluctuating asymmetry differed between the inbred and outbred groups. Male children of well-educated parents, however, were larger and had less fluctuating asymmetry. Female children of well-educated parents weighed more than those of less well-educated parents, but were otherwise indistinguishable for height, head circumference, chest circumference, and fluctuating asymmetry. We conclude that inbreeding depression causes neither an increase in fluctuating asymmetry of full-term newborns, nor a decrease in body size. Unmeasured variables correlated with education appear to have an effect on fluctuating asymmetry and size of male children and only a weak effect on size (weight) of female children.

MHC-I affects infection intensity but not infection status with a frequent avian malaria parasite in blue tits

Host resistance against parasites depends on three aspects: the ability to prevent, control and clear infections. In vertebrates the immune system consists of innate and adaptive immunity. Innate immunity is particularly important for preventing infection and eradicating established infections at an early stage while adaptive immunity is slow, but powerful, and essential for controlling infection intensities and eventually clearing infections. Major Histocompatibility Complex (MHC) molecules are central in adaptive immunity, and studies on parasite resistance and MHC in wild animals have found effects on both infection intensity (parasite load) and infection status (infected or not). It seems MHC can affect both the ability to control infection intensities and the ability to clear infections. However, these two aspects have rarely been considered simultaneously, and their relative importance in natural populations is therefore unclear. Here we investigate if MHC class I genotype affects infection intensity and infection status with a frequent avian malaria infection Haemoproteus majoris in a natural population of blue tits Cyanistes caeruleus. We found a significant negative association between a single MHC allele and infection intensity but no association with infection status. Blue tits that carry a specific MHC allele seem able to suppress H. majoris infection intensity, while we have no evidence that this allele also has an effect on clearance of the H. majoris infection, a result that is in contrast with some previous studies of MHC and avian malaria. A likely explanation could be that the clearance rate of avian malaria parasites differs between avian malaria lineages and/or between avian hosts.

An offer you cannot refuse: down-regulation of immunity in response to a pathogen's retaliation threat

According to the Red Queen hypothesis, hosts and pathogens are engaged in an escalating coevolutionary arms race between resistance and virulence. However, the vast majority of symbionts colonize their hosts' mucosal compartments without triggering any immune response, resulting in durable commensal associations. Here, I propose a simple extension of previous mathematical models for antagonistic coevolution in which the host can mount a delayed immune response; in response, the symbiont can change its virulence following this activation. Even though the levels of virulence in both phases are assumed to be genetically determined, this simple form of plasticity can select for commensal associations. In particular, coevolution can result in hosts that do not activate their immune response, thus preventing phenotypically plastic pathogens from switching to a higher virulence level. I argue that, from the host's point of view, this state is analogous to the mafia behaviour previously described in avian brood parasites. More importantly, this study provides a new hypothesis for the maintenance of a commensal relationship through antagonistic coevolution.

Parasite-mediated selection of major histocompatibility complex variability in wild brandt's voles (Lasiopodomys brandtii) from Inner Mongolia, China

Background

Genes of the major histocompatibility complex (MHC) exhibit high levels of variability, which is believed to have arisen through pathogen-mediated selection. We investigated the relationship between parasite load and genetic diversity at selectively neutral, non-coding markers (microsatellites) and adaptive genetic variation at a functionally important part of the MHC in six independent natural populations of Brandt’s voles (Lasiopodomys brandtii) from two regions of the Xilingol Grassland area of Inner Mongolia.

Results

Two-hundred and fifty-two voles were screened for gastrointestinal parasites, and were assessed for genetic variation. Parasite screening was done through non-invasive fecal egg counts, while allelic diversity was determined via single-stranded conformation polymorphism and DNA sequencing. We detected eight distinct helminth egg morphotypes. A total of 10 microsatellite loci were genotyped and 19 unique MHC class II B alleles were isolated. The rate of nonsynonymous substitutions (d_N) exceeded the rate of synonymous substitutions (d_S) at putative antigen binding sites of DRB. Neutral and adaptive genetic diversity differed between the six vole populations. To test the main pathogen-driven selection hypotheses for the maintenance of host MHC diversity and parasite species-specific co-evolutionary effects, multivariate approaches (generalized linear mixed models) were used to test for associations between the MHC class II DRB genotype and infections with nematodes. We found no evidence for heterozygote advantage, and overall heterozygosity was lower than expected in the MHC alleles. We identified an association between the parasite load and specific MHC alleles in the voles, and this pattern varied between geographic regions.

Conclusions

The results suggest that MHC variability in Brandt’s voles is maintained by rare allele advantage and fluctuating selection, but the data failed to show any heterozygote advantage effect. Our results add to a growing body of evidence showing that the mode and relative strength of pathogen-driven selection acting on MHC diversity varies within specific wild populations. In addition, our study contributes to the understanding of what maintains MHC diversity, of host-pathogen coevolution and of how genetic diversity is maintained in voles.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

Identification and characterization of the major histocompatibility complex class II DQB (MhcMath-DQB1) alleles in Tibetan macaques (Macaca thibetana)

Tibetan macaque (Macaca thibetana), an endangered primate species endemic to China, have been used as experimental animal model for various human diseases. Major histocompatibility complex (MHC) genes play a crucial role in the susceptibility and/or resistance to many human diseases, but little is known about Tibetan macaques. To gain an insight into the MHC background and to facilitate the experimental use of Tibetan macaques, the second exon of Mhc-DQB1 gene was sequenced in a cohort of wild Tibetan macaques living in the Sichuan province of China. A total of 23 MhcMath-DQB1 alleles were identified for the first time, illustrating a marked allelic polymorphism at the DQB1 locus for these macaques. Most of the sequences (74%) observed in this study belong to DQB1*06 (9 alleles) and DQB1*18 (8 alleles) lineages, and the rest (26%) belong to DQB1*15 (3 alleles) and DQB1*17 (3 alleles) lineages. The most frequent alleles detected among these macaques were MhcMath-DQB1*15:02:02 (17.9%), followed by Math-DQB1*06:06, 17:03 and 18:01, which were detected in 9 (16.1%) of the monkeys, respectively. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the peptide-binding region, suggesting balancing selection for maintaining polymorphisms at the MHC class II DQB1 locus. Phylogenetic analyses confirms the trans-species model of evolution of the Mhc-DQB1 genes in non-human primates, and in particular, the extensive allele sharing is observed between Tibetan and other macaque species.

Mate choice for genetic compatibility in the house mouse

In house mice, genetic compatibility is influenced by the t haplotype, a driving selfish genetic element with a recessive lethal allele, imposing fundamental costs on mate choice decisions. Here, we evaluate the cost of genetic incompatibility and its implication for mate choice in a wild house mice population. In laboratory reared mice, we detected no fertility (number of embryos) or fecundity (ability to conceive) costs of the t, and yet we found a high cost of genetic incompatibility: heterozygote crosses produced 40% smaller birth litter sizes because of prenatal mortality. Surprisingly, transmission of t in crosses using +/t males was influenced by female genotype, consistent with postcopulatory female choice for + sperm in +/t females. Analysis of paternity patterns in a wild population of house mice showed that +/t females were more likely than +/+ females to have offspring sired by +/+ males, and unlike +/+ females, paternity of their offspring was not influenced by +/t male frequency, further supporting mate choice for genetic compatibility. As the major histocompatibility complex (MHC) is physically linked to the t, we investigated whether females could potentially use variation at the MHC to identify male genotype at the sperm or individual level. A unique MHC haplotype is linked to the t haplotype. This MHC haplotype could allow the recognition of t and enable pre- and postcopulatory mate choice for genetic compatibility. Alternatively, the MHC itself could be the target of mate choice for genetic compatibility. We predict that mate choice for genetic compatibility will be difficult to find in many systems, as only weak fertilization biases were found despite an exceptionally high cost of genetic incompatibility.

Evolution of MHC class I genes in the endangered loggerhead sea turtle (Caretta caretta) revealed by 454 amplicon sequencing

Background

In evolutionary and conservation biology, parasitism is often highlighted as a major selective pressure. To fight against parasites and pathogens, genetic diversity of the immune genes of the major histocompatibility complex (MHC) are particularly important. However, the extensive degree of polymorphism observed in these genes makes it difficult to conduct thorough population screenings.

Methods

We utilized a genotyping protocol that uses 454 amplicon sequencing to characterize the MHC class I in the endangered loggerhead sea turtle (Caretta caretta) and to investigate their evolution at multiple relevant levels of organization.

Results

MHC class I genes revealed signatures of trans-species polymorphism across several reptile species. In the studied loggerhead turtle individuals, it results in the maintenance of two ancient allelic lineages. We also found that individuals carrying an intermediate number of MHC class I alleles are larger than those with either a low or high number of alleles.

Conclusions

Multiple modes of evolution seem to maintain MHC diversity in the loggerhead turtles, with relatively high polymorphism for an endangered species.

Mhc supertypes confer both qualitative and quantitative resistance to avian malaria infections in a wild bird population

Major histocompatibility complex (Mhc) genes are believed to play a key role in the genetic basis of disease control. Although numerous studies have sought links between Mhc and disease prevalence, many have ignored the ecological and epidemiological aspects of the host-parasite interaction. Consequently, interpreting associations between prevalence and Mhc has been difficult, whereas discriminating alleles for qualitative resistance, quantitative resistance and susceptibility remains challenging. Moreover, most studies to date have quantified associations between genotypes and disease status, overlooking the complex relationship between genotype and the properties of the Mhc molecule that interacts with parasites. Here, we address these problems and demonstrate avian malaria (Plasmodium) parasite species-specific associations with functional properties of Mhc molecules (Mhc supertypes) in a wild great tit (Parus major) population. We further show that correctly interpreting these associations depends crucially on understanding the spatial variation in risk of infection and the fitness effects of infection. We report that a single Mhc supertype confers qualitative resistance to Plasmodium relictum, whereas a different Mhc supertype confers quantitative resistance to Plasmodium circumflexum infections. Furthermore, we demonstrate common functional properties of Plasmodium-resistance alleles in passerine birds, suggesting this is a model system for parasite-Mhc associations in the wild.

Tick infestation patterns in free ranging African buffalo (Syncercus caffer): Effects of host innate immunity and niche segregation among tick species

Ticks are of vast importance to livestock health, and contribute to conflicts between wildlife conservation and agricultural interests; but factors driving tick infestation patterns on wild hosts are not well understood. We studied tick infestation patterns on free-ranging African buffalo (Syncercus caffer), asking (i) is there evidence for niche segregation among tick species?; and (ii) how do host characteristics affect variation in tick abundance among hosts? We identified ticks and estimated tick burdens on 134 adult female buffalo from two herds at Kruger National Park, South Africa. To assess niche segregation, we evaluated attachment site preferences and tested for correlations between abundances of different tick species. To investigate which host factors may drive variability in tick abundance, we measured age, body condition, reproductive and immune status in all hosts, and examined their effects on tick burdens.

Two tick species were abundant on buffalo, Amblyomma hebraeum and Rhipicephalus evertsi evertsi. A. hebraeum were found primarily in the inguinal and axillary regions; R. e. evertsi attached exclusively in the perianal area. Abundances of A. hebraeum and R. e. evertsi on the host were unrelated. These results suggest spatial niche segregation between A. hebraeum and R. e. evertsi on the buffalo. Buffalo with stronger innate immunity, and younger buffalo, had fewer ticks. Buffalo with low body condition scores, and pregnant buffalo, had higher tick burdens, but these effects varied between the two herds we sampled.

This study is one of the first to link ectoparasite abundance patterns and immunity in a free-ranging mammalian host population. Based on independent abundances of A. hebraeum and R. e. evertsi on individual buffalo, we would expect no association between the diseases these ticks transmit. Longitudinal studies linking environmental variability with host immunity are needed to understand tick infestation patterns and the dynamics of tick-borne diseases in wildlife.

A sequence-based approach demonstrates that balancing selection in classical human leukocyte antigen (HLA) loci is asymmetric

Balancing selection has maintained human leukocyte antigen (HLA) allele diversity, but it is unclear whether this selection is symmetric (all heterozygotes are comparable and all homozygotes are comparable in terms of fitness) or asymmetric (distinct heterozygote genotypes display greater fitness than others). We tested the hypothesis that HLA is under asymmetric balancing selection in populations by estimating allelic branch lengths from genetic sequence data encoding peptide-binding domains. Significant deviations indicated changes in the ratio of terminal to internal branch lengths. Such deviations could arise even if no individual alleles present a strikingly altered branch length (e.g. if there is an overall distortion, with all or many terminal branches being longer than expected). DQ and DP loci were also analyzed as haplotypes. Using allele frequencies for 419 distinct populations in 10 geographical regions, we examined population differentiation in alleles within and between regions, and the relationship between allelic branch length and frequency. The strongest evidence for asymmetrical balancing selection was observed for HLA-DRB1, HLA-B and HLA-DPA1, with significant deviation (P ≤ 1.1 × 10(-4)) in about half of the populations. There were significant results at all loci except HLA-DQB1/DQA1. We observed moderate genetic variation within and between geographic regions, similar to the rest of the genome. Branch length was not correlated with allele frequency. In conclusion, sequence data suggest that balancing selection in HLA is asymmetric (some heterozygotes enjoy greater fitness than others). Because HLA polymorphism is crucial for pathogen resistance, this may manifest as a frequency-dependent selection with fluctuation in the fitness of specific heterozygotes over time.

Variation in MHC genotypes in two populations of house sparrow (Passer domesticus) with different population histories

Small populations are likely to have a low genetic ability for disease resistance due to loss of genetic variation through inbreeding and genetic drift. In vertebrates, the highest genetic diversity of the immune system is located at genes within the major histocompatibility complex (MHC). Interestingly, parasite-mediated selection is thought to potentially maintain variation at MHC loci even in populations that are monomorphic at other loci. Therefore, general loss of genetic variation in the genome may not necessarily be associated with low variation at MHC loci. We evaluated inter- and intrapopulation variation in MHC genotypes between an inbred (Aldra) and a relatively outbred population (Hestmannøy) of house sparrows (Passer domesticus) in a metapopulation at Helgeland, Norway. Genomic (gDNA) and transcribed (cDNA) alleles of functional MHC class I and IIB loci, along with neutral noncoding microsatellite markers, were analyzed to obtain relevant estimates of genetic variation. We found lower allelic richness in microsatellites in the inbred population, but high genetic variation in MHC class I and IIB loci in both populations. This suggests that also the inbred population could be under balancing selection to maintain genetic variation for pathogen resistance.

Ancestral polymorphism at the major histocompatibility complex (MHCIIß) in the Nesospiza bunting species complex and its sister species (Rowettia goughensis)

BACKGROUND

The major histocompatibility complex (MHC) is an important component of the vertebrate immune system and is frequently used to characterise adaptive variation in wild populations due to its co-evolution with pathogens. Passerine birds have an exceptionally diverse MHC with multiple gene copies and large numbers of alleles compared to other avian taxa. The Nesospiza bunting species complex (two species on Nightingale Island; one species with three sub-species on Inaccessible Island) represents a rapid adaptive radiation at a small, isolated archipelago, and is thus an excellent model for the study of adaptation and speciation. In this first study of MHC in Nesospiza buntings, we aim to characterize MHCIIß variation, determine the strength of selection acting at this gene region and assess the level of shared polymorphism between the Nesospiza species complex and its putative sister taxon, Rowettia goughensis, from Gough Island.

RESULTS

In total, 23 unique alleles were found in 14 Nesospiza and 2 R. goughensis individuals encoding at least four presumably functional loci and two pseudogenes. There was no evidence of ongoing selection on the peptide binding region (PBR). Of the 23 alleles, 15 were found on both the islands inhabited by Nesospiza species, and seven in both Nesospiza and Rowettia; indications of shared, ancestral polymorphism. A gene tree of Nesospiza MHCIIß alleles with several other passerine birds shows three highly supported Nesospiza-specific groups. All R. goughensis alleles were shared with Nesospiza, and these alleles were found in all three Nesospiza sequence groups in the gene tree, suggesting that most of the observed variation predates their phylogenetic split.

CONCLUSIONS

Lack of evidence of selection on the PBR, together with shared polymorphism across the gene tree, suggests that population variation of MHCIIß among Nesospiza and Rowettia is due to ancestral polymorphism rather than local selective forces. Weak or no selection pressure could be attributed to low parasite load at these isolated Atlantic islands. The deep divergence between the highly supported Nesospiza-specific sequence Groups 2 and 3, and the clustering of Group 3 close to the distantly related passerines, provide strong support for preserved ancestral polymorphism, and present evidence of one of the rare cases of extensive ancestral polymorphism in birds.

Evidence for the 'good genes' model: association of MHC class II DRB alleles with ectoparasitism and reproductive state in the neotropical lesser bulldog bat, Noctilio albiventris

The adaptive immune system has a major impact on parasite resistance and life history strategies. Immunological defence is costly both in terms of immediate activation and long-term maintenance. The 'good genes' model predicts that males with genotypes that promote a good disease resistance have the ability to allocate more resources to reproductive effort which favours the transmission of good alleles into future generations. Our study shows a correlation between immune gene constitution (Major Histocompatibility Complex, MHC class II DRB), ectoparasite loads (ticks and bat flies) and the reproductive state in a neotropical bat, Noctilio albiventris. Infestation rates with ectoparasites were linked to specific Noal-DRB alleles, differed among roosts, increased with body size and co-varied with reproductive state particularly in males. Non-reproductive adult males were more infested with ectoparasites than reproductively active males, and they had more often an allele (Noal-DRB*02) associated with a higher tick infestation than reproductively active males or subadults. We conclude that the individual immune gene constitution affects ectoparasite susceptibility, and contributes to fitness relevant trade-offs in male N. albiventris as suggested by the 'good genes' model.