Shared evolutionary origin of major histocompatibility complex polymorphism in sympatric lemurs

Genes of the major histocompatibility complex (MHC) play a central role in adaptive immune responses of vertebrates. They exhibit remarkable polymorphism, often crossing species boundaries with similar alleles or allelic motifs shared across species. This pattern may reflect parallel parasite‐mediated selective pressures, either favouring the long maintenance of ancestral MHC allelic lineages across successive speciation events by balancing selection (“trans‐species polymorphism”), or alternatively favouring the independent emergence of functionally similar alleles post‐speciation via convergent evolution. Here, we investigate the origins of MHC similarity across several species of dwarf and mouse lemurs (Cheirogaleidae). We examined MHC class II variation in two highly polymorphic loci (DRB, DQB) and evaluated the overlap of gut–parasite communities in four sympatric lemurs. We tested for parasite‐MHC associations across species to determine whether similar parasite pressures may select for similar MHC alleles in different species. Next, we integrated our MHC data with those previously obtained from other Cheirogaleidae to investigate the relative contribution of convergent evolution and co‐ancestry to shared MHC polymorphism by contrasting patterns of codon usage at functional vs. neutral sites. Our results indicate that parasites shared across species may select for functionally similar MHC alleles, implying that the dynamics of MHC‐parasite co‐evolution should be envisaged at the community level. We further show that balancing selection maintaining trans‐species polymorphism, rather than convergent evolution, is the primary mechanism explaining shared MHC sequence motifs between species that diverged up to 30 million years ago.     

Evidence of gene orthology and trans‐species polymorphism,but not of parallel evolution,despite high levels of concerted evolution in the major histocompatibility complex of flamingo species

The major histocompatibility complex (MHC) is a cornerstone in the study of adaptive genetic diversity. Intriguingly, highly polymorphic MHC sequences are often not more similar within species than between closely related species. Divergent selection of gene duplicates, balancing selection maintaining trans‐species polymorphism (TSP) that predate speciation and parallel evolution of species sharing similar selection pressures can all lead to higher sequence similarity between species. In contrast, high rates of concerted evolution increase sequence similarity of duplicated loci within species. Assessing these evolutionary models remains difficult as relatedness and ecological similarities are often confounded. As sympatric species of flamingos are more distantly related than allopatric species, flamingos represent an ideal model to disentangle these evolutionary models. We characterized MHC Class I exon 3, Class IIB exon 2 and exon 3 of the six extant flamingo species. We found up to six MHC Class I loci and two MHC Class IIB loci. As all six species shared the same number of MHC Class IIB loci, duplication appears to predate flamingo speciation. However, the high rate of concerted evolution has prevented the divergence of duplicated loci. We found high sequence similarity between all species regardless of codon position. The latter is consistent with balancing selection maintaining TSP, as under this mechanism amino acid sites under pathogen‐mediated selection should be characterized by fewer synonymous codons (due to their common ancestry) than under parallel evolution. Overall, balancing selection maintaining TSP appears to result in high MHC similarity between species regardless of species relatedness and geographical distribution.     

Multiple parasites are driving major histocompatibility complex polymorphism in the wild

Parasite mediated selection may result in arms races between host defence and parasite virulence. In particular, simultaneous infections from multiple parasite species should cause diversification (i.e. balancing selection) in resistance genes both at the population and the individual level. Here, we tested these ideas in highly polymorphic major histocompatibility complex (MHC) genes from three-spined sticklebacks (Gasterosteus aculeatus L.). In eight natural populations, parasite diversity (15 different species), and MHC class IIB diversity varied strongly between habitat types (lakes vs. rivers vs. estuaries) with lowest values in rivers. Partial correlation analysis revealed an influence of parasite diversity on MHC class IIB variation whereas general genetic diversity assessed at seven microsatellite loci was not significantly correlated with parasite diversity. Within individual fish, intermediate, rather than maximal allele numbers were associated with minimal parasite load, supporting theoretical models of self-reactive T-cell elimination. The optimal individual diversity matched those values female fish try to achieve in their offspring by mate choice. We thus present correlative evidence supporting the 'allele counting' strategy for optimizing the immunocompetence in stickleback offspring.     

Historical and contemporary selection of teleost MHC genes: did we leave the past behind?

The extreme polymorphism of antigen‐presenting genes of the major histocompatibility complex (MHC) has spurred intense research unparalleled for any other gene family. This applies also to teleosts where sequence information is available for 3559 MHC class I and class II allelic variants from 137 species. This review summarizes current knowledge on the origin and maintenance of diversity at classical MHC loci. Most studies identified positive selection (i.e. elevated rates of non‐synonymous over synonymous substitutions, dN/dS) as a sign of balancing selection. A meta‐analysis on nine species with sufficient numbers of class I and class II sequences revealed that recombination rate and intensity of positive selection were positively correlated, suggesting that recombination and gene conversion played a significant role in shaping the allelic repertoire. Processes that create diversity over long timescales need to be complemented by contemporary balancing selection, either through overdominance or frequency‐dependent selection, in order to explain the high allelic diversity observed today. While some evidence for overdominance exists for a few taxa (mainly salmonids) by correlating parasite infection data or survival to MHC genotypes, field or experimental data on negative frequency‐dependent selection are lacking altogether, even though some fish species are particularly suitable as model systems. Theoretical predictions suggest that negative frequency‐dependent selection is necessary to maintain the existing polymorphism. Hence, future empirical studies should focus on detecting signals that differentiate between mechanisms of contemporary selection rather than repeatedly showing historical selection events.     

The role of major histocompatibility complex diversity in vigour of fish males (Abramis brama L.) and parasite selection

The major histocompatibility complex (MHC) presents a group of genes with highly polymorphic loci involved in specific immune responses. The factors maintaining extensive MHC polymorphism have been questioned, considering three possible hypotheses of parasite‐mediated selection driving an extensive MHC diversity (i.e. heterozygote advantage, rare‐allele advantage, and favouring optimal MHC diversity). The patterns of MHC diversity of class IIB genes were investigated following two noncontradicting hypotheses, parasite‐driven selection and MHC‐based mating preferences, using males of common bream collected in the spawning period. Two allelic groups DAB1 and DAB3 were recognized from the phylogenetic analyses. Individuals expressed one or two alleles of the same or different allelic groups. Several individuals shared identical alleles; however, the presence of parasite species was not associated with the occurrence of a particular allele. The presence of different allelic groups (only DAB1, only DAB3, or both DAB1 and DAB3) in individuals was not associated with parasite presence or diversity. The expression of two DAB1 alleles was associated with higher endoparasite abundance. Moreover, nucleotide diversity in individuals expressing a single type of alleles (DAB1 or DAB3) increased with the abundance of ectoparasitic Dactylogyrus spp. (Monogenea) and Ergasilus sp. (Crustacea). This suggests that the expression of two alleles of a single allelic type is related to high metazoan parasite infection whereas no significant influence of parasitism on the combined allelic form (the presence of both DAB1 and DAB3 alleles) was found. Moreover, the expression of two alleles of a single allelic type was related to decreased immunocompetence measured by spleen size. The condition factor was higher in fish expressing the combined allelic type. Thus, the presence of alleles of different lineages in individuals appears to be advantageous for individual male fitness. The expression of a single allelic type was related to higher sexual ornamentation, which could support the role of MHC in the hypothesis of the sexual selection of ‘good genes’. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society, 2007, 90 , 525–538.     

Selection from parasites favours immunogenetic diversity but not divergence among locally adapted host populations

The unprecedented polymorphism in the major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection from parasites. However, do parasites also drive divergence at MHC loci between host populations, or do the effects of balancing selection maintain similarities among populations? We examined MHC variation in populations of the livebearing fish Poecilia mexicana and characterized their parasite communities. Poecilia mexicana populations in the Cueva del Azufre system are locally adapted to darkness and the presence of toxic hydrogen sulphide, representing highly divergent ecotypes or incipient species. Parasite communities differed significantly across populations, and populations with higher parasite loads had higher levels of diversity at class II MHC genes. However, despite different parasite communities, marked divergence in adaptive traits and in neutral genetic markers, we found MHC alleles to be remarkably similar among host populations. Our findings indicate that balancing selection from parasites maintains immunogenetic diversity of hosts, but this process does not promote MHC divergence in this system. On the contrary, we suggest that balancing selection on immunogenetic loci may outweigh divergent selection causing divergence, thereby hindering host divergence and speciation. Our findings support the hypothesis that balancing selection maintains MHC similarities among lineages during and after speciation (trans‐species evolution).     

Genetic variation in MHC class II expression and interactions with MHC sequence polymorphism in three-spined sticklebacks

Genes of the major histocompatibility complex (MHC) have been studied for several decades because of their pronounced allelic polymorphism. Structural allelic polymorphism is, however, not the only source of variability subjected to natural selection. Genetic variation may also exist in gene expression patterns. Here, we show that in a natural population of three-spined sticklebacks (Gasterosteus aculeatus) the expression of MHC class IIB genes was positively correlated with parasite load, which indicates increased immune activation of the MHC when infections are frequent. To experimentally study MHC expression, we used laboratory-bred sticklebacks that were exposed to three naturally occurring species of parasite. We found strong differences in MHC class IIB expression patterns among fish families, which were consistent over two generations, thus demonstrating a genetic component. The average number of MHC class IIB sequence variants within families was negatively correlated to the MHC expression level suggesting compensatory up-regulation in fish with a low (i.e. suboptimal) MHC sequence variability. The observed differences among families and the negative correlation with individual sequence diversity imply that MHC expression is evolutionary relevant for the onset and control of the immune response in natural populations.     

Host association influences variation at salivary protein genes in the bat ectoparasite Cimex adjunctus

Parasite–host relationships create strong selection pressures that can lead to adaptation and increasing specialization of parasites to their hosts. Even in relatively loose host–parasite relationships, such as between generalist ectoparasites and their hosts, we may observe some degree of specialization of parasite populations to one of the multiple potential hosts. Salivary proteins are used by blood‐feeding ectoparasites to prevent hemostasis in the host and maximize energy intake. We investigated the influence of association with specific host species on allele frequencies of salivary protein genes in Cimex adjunctus, a generalist blood‐feeding ectoparasite of bats in North America. We analysed two salivary protein genes: an apyrase, which hydrolyses ATP at the feeding site and thus inhibits platelet aggregation, and a nitrophorin, which brings nitrous oxide to the feeding site, inhibiting platelet aggregation and vasoconstriction. We observed more variation at both salivary protein genes among parasite populations associated with different host species than among populations from different spatial locations associated with the same host species. The variation in salivary protein genes among populations on different host species was also greater than expected under a neutral scenario of genetic drift and gene flow. Finally, host species was an important predictor of allelic divergence in genotypes of individual C. adjunctus at both salivary protein genes. Our results suggest differing selection pressures on these two salivary protein genes in C. adjunctus depending on the host species.     

Coevolutionary relationship between helminth diversity and MHC class II polymorphism in rodents

Parasite-mediated selection on major histocompatibility complex (MHC) genes has mainly been explored at the intraspecific level, although many molecular studies have revealed trans-species polymorphism. Interspecific patterns of MHC diversity might reveal factors responsible for the long-term evolution of MHC polymorphism. We hypothesize that host taxa harbouring high parasite diversity should exhibit high levels of MHC genetic diversity. We test this assumption using data on rodent species and their helminth parasites compiled from the literature. Controlling for similarity due to common descent, we present evidence indicating that high helminth species richness in rodent species is associated with increased MHC class II polymorphism. Our results are consistent with the idea that parasites sharing a long-term coevolutionary history with their hosts are the agents of selection explaining MHC polymorphism.     

Neutral and selective processes shape MHC gene diversity and expression in stocked brook charr populations (Salvelinus fontinalis)

The capacity of an individual to battle infection is an important fitness determinant in wild vertebrate populations. The major histocompatibility complex (MHC) genes are crucial for a host's adaptive immune system to detect pathogens. However, anthropogenic activities may disrupt natural cycles of co‐evolution between hosts and pathogens. In this study, we investigated the dynamic sequence and expression variation of host parasite interactions in brook charr (Salvelinus fontinalis) in a context of past human disturbance via population supplementation from domestic individuals. To do so, we developed a new method to examine selection shaping MHC diversity within and between populations and found a complex interplay between neutral and selective processes that varied between lakes that were investigated. We provided evidence for a lower introgression rate of domestic alleles and found that parasite infection increased with domestic genomic background of individuals. We also documented an association between individual MHC alleles and parasite taxa. Finally, longer cis‐regulatory minisatellites were positively correlated with MHC II down‐regulation and domestic admixture, suggesting that inadvertent selection during domestication resulted in a lower immune response capacity, through a trade‐off between growth and immunity, which explained the negative selection of domestic alleles at least under certain circumstances.     

Polymorphic microsatellite loci for the trematode Diplostomum pseudospathaceum

We present primers for five polymorphic microsatellite loci in the eye fluke Diplostomum pseudospathaceum (Trematoda), a widely distributed parasite with a complex life cycle used as a model for parasitology and fish immunology. The loci were identified using a GA/CT‐enriched genomic library by subtractive hybridization with magnetic particles. All five loci were highly polymorphic, displaying 17 to 61 alleles and heterozygosities ranging from 0.53 to 0.92. We isolated populations of parasites within the first (snail) and second (fish) intermediate host and found small but significant genetic differentiation (F_ST = 0.012) between the two life stages of the parasite.     

Evolution of MHC class II SLA‐DRB1 locus in the Croatian wild boar (Sus scrofa) implies duplication and weak signals of positive selection

The wild boar is an ancestor of the domestic pig and an important game species with the widest geographical range of all ungulates. Although a large amount of data are available on major histocompatibility complex (MHC) variability in domestic pigs, only a few studies have been performed on wild boars. Due to their crucial role in appropriate immune responses and extreme polymorphism, MHC genes represent some of the best candidates for studying the processes of adaptive evolution. Here, we present the results on the variability and evolution of the entire MHC class II SLA‐DRB1 locus exon 2 in 133 wild boars from Croatia. Using direct sequencing and cloning methods, we identified 20 SLA‐DRB1 alleles, including eight new variants, with notable divergence. In some individuals, we documented functional locus duplication, and SLA‐DRB1*04:10 was identified as the allele involved in the duplication. The expression of a duplicated locus was confirmed by cloning and sequencing cDNA‐derived amplicons. Based on individual genotypes, we were able to assume that alleles SLA‐DRB1*04:10 and SLA‐DRB1*06:07 are linked as an allelic combination that co‐evolves as a two‐locus haplotype. Our investigation of evolutionary processes at the SLA‐DRB1 locus confirmed the role of intralocus recombination in generating allelic variability, whereas tests of positive selection based on the dN/dS (non‐synonymous/synonymous substitution rate ratio) test revealed atypically weak and ambiguous signals.     

Footprints of ancient‐balanced polymorphisms in genetic variation data from closely related species

When long‐lasting, balancing selection can lead to “trans‐species” polymorphisms that are shared by two or more species identical by descent. In such cases, the gene genealogy at the selected site clusters by allele instead of by species, and nearby neutral sites also have unusual genealogies because of linkage. While this scenario is expected to leave discernible footprints in genetic variation data, the specific patterns remain poorly characterized. Motivated by recent findings in primates, we focus on the case of a biallelic polymorphism under ancient balancing selection and derive approximations for summaries of the polymorphism data from two species. Specifically, we characterize the length of the segment that carries most of the footprints, the expected number of shared neutral single nucleotide polymorphisms (SNPs), and the patterns of allelic associations among them. We confirm the accuracy of our approximations by coalescent simulations. We further show that for humans and chimpanzees—more generally, for pairs of species with low genetic diversity levels—these patterns are highly unlikely to be generated by neutral recurrent mutations. We discuss the implications for the design and interpretation of genome scans for ancient balanced polymorphisms in primates and other taxa.     

Blood parasites shape extreme major histocompatibility complex diversity in a migratory passerine

Pathogens are one of the main forces driving the evolution and maintenance of the highly polymorphic genes of the vertebrate major histocompatibility complex (MHC). Although MHC proteins are crucial in pathogen recognition, it is still poorly understood how pathogen‐mediated selection promotes and maintains MHC diversity, and especially so in host species with highly duplicated MHC genes. Sedge warblers (Acrocephalus schoenobaenus) have highly duplicated MHC genes, and using data from high‐throughput MHC genotyping, we were able to investigate to what extent avian malaria parasites explain temporal MHC class I supertype fluctuations in a long‐term study population. We investigated infection status and infection intensities of two different strains of Haemoproteus, that is avian malaria parasites that are known to have significant fitness consequences in sedge warblers. We found that prevalence of avian malaria in carriers of specific MHC class I supertypes was a significant predictor of their frequency changes between years. This finding suggests that avian malaria infections partly drive the temporal fluctuations of the MHC class I supertypes. Furthermore, we found that individuals with a large number of different supertypes had higher resistance to avian malaria, but there was no evidence for an optimal MHC class I diversity. Thus, the two studied malaria parasite strains appear to select for a high MHC class I supertype diversity. Such selection may explain the maintenance of the extremely high number of MHC class I gene copies in sedge warblers and possibly also in other passerines where avian malaria is a common disease.     

Adaptive evolution of lateral plates in three‐spined stickleback Gasterosteus aculeatus: a case study in functional analysis of natural variation

The three‐spined stickleback Gasterosteus aculeatus is a model species for studying questions in ecology and evolution. The rapid diversification of G. aculeatus in post‐glacial freshwater environments, combined with recently developed molecular tools, provides a unique opportunity to study the functional basis of fitness variation in natural populations. In derived freshwater populations, a number of morphological traits have diverged in parallel from the marine ancestral state, including the number of lateral armour plates. Evolution of reduced armour in freshwater populations is due to positive selection from both abiotic and biotic mechanisms. The major effect gene (ectodysplasin‐A or Eda), along with several minor effect genetic regions, has recently been shown to control lateral plate variation. Field experiments have further determined the fitness consequences of allelic variation at the major effect locus. This work helps elucidate the mechanisms connecting genetic variation with phenotypic variation and fitness in the wild, a synthesis that should be applicable to many other phenotypic traits and species of fishes.     

Weak link between dispersal and parasite community differentiation or immunogenetic divergence in two sympatric cichlid fishes

Geographical isolation, habitat variation and trophic specialization have contributed to a large extent to the astonishing diversity of cichlid fishes in the Great East African lakes. Because parasite communities often vary across space and environments, parasites can accompany and potentially enhance cichlid species diversification. However, host dispersal may reduce opportunities for parasite‐driven evolution by homogenizing parasite communities and allele frequencies of immunity genes. To test for the relationships between parasite community variation, host dispersal and parasite‐induced host evolution, we studied two sympatric cichlid species with contrasting dispersal capacities along the shores of southern Lake Tanganyika. Whereas the philopatric Tropheus moorii evolved into several genetically differentiated colour morphs, Simochromis diagramma is phenotypically rather uniform across its distribution range and shows only weak population structure. Populations of both species were infected with divergent parasite communities and harbour differentiated variant pools of an important set of immune genes, the major histocompatibility complex (MHC). The overall extent of geographical variation of parasites and MHC genes was similar between host species. This indicates that immunogenetic divergence among populations of Lake Tanganyika cichlids can occur even in species that are strongly dispersing. However, because this also includes species that are phenotypically uniform, parasite‐induced evolution may not represent a key factor underlying species diversification in this system.     

Variety matters: adaptive genetic diversity and parasite load in two mouse opossums from the Brazilian Atlantic forest

The adaptive potential of a species to a changing environment and in disease defence is primarily based on genetic variation. Immune genes, such as genes of the major histocompatibility complex (MHC), may thereby be of particular importance. In marsupials, however, there is very little knowledge about natural levels and functional importance of MHC polymorphism, despite their key role in the mammalian evolution. In a previous study, we discovered remarkable differences in the MHC class II diversity between two species of mouse opossums (Gracilinanus microtarsus, Marmosops incanus) from the Brazilian Atlantic forest, which is one of the most endangered hotspots for biodiversity conservation. Since the main forces in generating MHC diversity are assumed to be pathogens, we investigated in this study gastrointestinal parasite burden and functional associations between the individual MHC constitution and parasite load. We tested two contrasting scenarios, which might explain differences in MHC diversity between species. We predicted that a species with low MHC diversity would either be under relaxed selection pressure by low parasite diversity (‘Evolutionary equilibrium’ scenario), or there was a recent loss in MHC diversity leading to a lack of resistance alleles and increased parasite burden (‘Unbalanced situation’ scenario). In both species it became apparent that the MHC class II is functionally important in defence against gastrointestinal helminths, which was shown here for the first time in marsupials. On the population level, parasite diversity did not markedly differ between the two host species. However, we did observe considerable differences in the individual parasite load (parasite prevalence and infection intensity): while M. incanus revealed low MHC DAB diversity and high parasite load, G. microtarsus showed a tenfold higher population wide MHC DAB diversity and lower parasite burden. These results support the second scenario of an unbalanced situation.     

Spatio‐temporal variation in parasite communities maintains diversity at the major histocompatibility complex class IIβ in the endangered Rio Grande silvery minnow

Climate change will strongly impact aquatic ecosystems particularly in arid and semi‐arid regions. Fish–parasite interactions will also be affected by predicted altered flow and temperature regimes, and other environmental stressors. Hence, identifying environmental and genetic factors associated with maintaining diversity at immune genes is critical for understanding species’ adaptive capacity. Here, we combine genetic (MHC class IIβ and microsatellites), parasitological and ecological data to explore the relationship between these factors in the remnant wild Rio Grande silvery minnow (Hybognathus amarus) population, an endangered species found in the southwestern United States. Infections with multiple parasites on the gills were observed and there was spatio‐temporal variation in parasite communities and patterns of infection among individuals. Despite its highly endangered status and chronically low genetic effective size, Rio Grande silvery minnow had high allelic diversity at MHC class IIβ with more alleles recognized at the presumptive DAB1 locus compared to the DAB3 locus. We identified significant associations between specific parasites and MHC alleles against a backdrop of generalist parasite prevalence. We also found that individuals with higher individual neutral heterozygosity and higher amino acid divergence between MHC alleles had lower parasite abundance and diversity. Taken together, these results suggest a role for fluctuating selection imposed by spatio‐temporal variation in pathogen communities and divergent allele advantage in maintenance of high MHC polymorphism. Understanding the complex interaction of habitat, pathogens and immunity in protected species will require integrated experimental, genetic and field studies.     

Gene duplication,allelic diversity,selection processes and adaptive value of MHC class II <Emphasis Type="Italic">DRB</Emphasis> genes of the bank vole, <Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.