Natural selection at major histocompatibility complex loci of vertebrates

The dynamics of reciprocal selective sweeps of host resistance and a parasite counter-adaptation in Drosophila

Host-parasite coevolution can result in consecutive selective sweeps of host resistance alleles and parasite counter-adaptations. To illustrate the dynamics of this important but little studied form of coevolution, we have modeled an ongoing arms race between Drosophila melanogaster and the vertically transmitted sigma virus, using parameters we estimated in the field. We integrate these results with previous work showing that the spread of a resistance allele of the ref(2)P gene in the host was followed by the spread of a virus genotype, which overcomes this resistance. In line with these observations, our model predicts that there can be rapid selective sweeps in both the host and parasite, which can drive large changes in the prevalence of infection. The virus will tend to be ahead in the arms race, as incomplete dominance slows down host adaptation and selection for host resistance is weaker than selection for parasites to overcome resistance--the "life-dinner" principle. This asymmetry in the adaptation rates results in a partial sweep of the host resistance allele, as it loses its advantage part way through the selective sweep. This well-understood natural system illustrates how the outcome of host-parasite coevolution is determined by different population genetic parameters in the field.

Cardiac structural and sarcomere genes associated with cardiomyopathy exhibit marked intolerance of genetic variation

BACKGROUND

The clinical significance of variants in genes associated with inherited cardiomyopathies can be difficult to determine because of uncertainty regarding population genetic variation and a surprising amount of tolerance of the genome even to loss-of-function variants. We hypothesized that genes associated with cardiomyopathy might be particularly resistant to the accumulation of genetic variation.

METHODS AND RESULTS

We analyzed the rates of single nucleotide genetic variation in all known genes from the exomes of >5000 individuals from the National Heart, Lung, and Blood Institute's Exome Sequencing Project, as well as the rates of structural variation from the Database of Genomic Variants. Most variants were rare, with over half unique to 1 individual. Cardiomyopathy-associated genes exhibited a rate of nonsense variants, about 96.1% lower than other Mendelian disease genes. We tested the ability of in silico algorithms to distinguish between a set of variants in MYBPC3, MYH7, and TNNT2 with strong evidence for pathogenicity and variants from the Exome Sequencing Project data. Algorithms based on conservation at the nucleotide level (genomic evolutionary rate profiling, PhastCons) did not perform as well as amino acid-level prediction algorithms (Polyphen-2, SIFT). Variants with strong evidence for disease causality were found in the Exome Sequencing Project data at prevalence higher than expected.

CONCLUSIONS

Genes associated with cardiomyopathy carry very low rates of population variation. The existence in population data of variants with strong evidence for pathogenicity suggests that even for Mendelian disease genetics, a probabilistic weighting of multiple variants may be preferred over the single gene causality model.

Sequencing of the core MHC region of black grouse (Tetrao tetrix) and comparative genomics of the galliform MHC

BACKGROUND

The MHC, which is regarded as the most polymorphic region in the genomes of jawed vertebrates, plays a central role in the immune system by encoding various proteins involved in the immune response. The chicken MHC-B genomic region has a highly streamlined gene content compared to mammalian MHCs. Its core region includes genes encoding Class I and Class IIB molecules but is only ~92Kb in length. Sequences of other galliform MHCs show varying degrees of similarity as that of chicken. The black grouse (Tetrao tetrix) is a wild galliform bird species which is an important model in conservation genetics and ecology. We sequenced the black grouse core MHC-B region and combined this with available data from related species (chicken, turkey, gold pheasant and quail) to perform a comparative genomics study of the galliform MHC. This kind of analysis has previously been severely hampered by the lack of genomic information on avian MHC regions, and the galliformes is still the only bird lineage where such a comparison is possible.

RESULTS

In this study, we present the complete genomic sequence of the MHC-B locus of black grouse, which is 88,390 bp long and contains 19 genes. It shows the same simplicity as, and almost perfect synteny with, the corresponding genomic region of chicken. We also use 454-transcriptome sequencing to verify expression in 17 of the black grouse MHC-B genes. Multiple sequence inversions of the TAPBP gene and TAP1-TAP2 gene block identify the recombination breakpoints near the BF and BLB genes. Some of the genes in the galliform MHC-B region also seem to have been affected by selective forces, as inferred from deviating phylogenetic signals and elevated rates of non-synonymous nucleotide substitutions.

CONCLUSIONS

We conclude that there is large synteny between the MHC-B region of the black grouse and that of other galliform birds, but that some duplications and rearrangements have occurred within this lineage. The MHC-B sequence reported here will provide a valuable resource for future studies on the evolution of the avian MHC genes and on links between immunogenetics and ecology of black grouse.

Association between the polymorphism of CfPGRP-S1 gene and disease susceptibility/resistance of zhikong scallop (Chlamys farreri) to Listonella anguillarum challenge

Peptidoglycan recognition protein (PGRP) is a pattern recognition receptor, playing important roles in the innate immune response against invasive pathogens. The single nucleotide polymorphism (SNP) loci in scallop PGRP gene (CfPGRP) were screened from Chlamys farreri to investigate their association with disease resistance of scallop against Listonella anguillarum. Thirteen SNP sites were identified in PGRP domain of CfPGRP, and two of them at positions 4407 and 4408 which are located in the same codon resulted in a nonsynonymous substitution. The genotype frequency of CG/CG in the resistant stock was significantly lower than that in susceptible stock (0% vs 32.4%), while that of CG/TA in the resistant stock was significantly higher than that in susceptible stock (P < 0.01). The pathogen-associated molecular patterns (PAMP) binding activity of two recombinant proteins, rCfPGRP-S1 (R) with CG variant in 4407-4408 site, rCfPGRP-S1 (Y) with TA variant in 4407-4408 site, were elucidated by examining their P/N value at 405 nm with ELISA assay. The in vitro binding activities of the two rCfPGRP-S1 variants to both lipopolysaccharide (LPS) and peptidoglycan (PGN) varied (P < 0.05) in a dose-dependent manner, and rCfPRPP-S1(Y) exhibited significantly higher affinity to PGN and LPS than that of rCfPGRP-S1(R) (P < 0.05). The growth inhibition assay was conducted to find the antibacterial activities of the two variants. Both rCfPGRP-S1(R) and rCfPGRP-S1 (Y) displayed obvious activity to suppress the growth of Escherichia coli, but there was no significant difference in suppression activity of two variants (P > 0.05). The results suggested that the polymorphism at locus 4407-4408 of CfPGRP-S1 considerably affected its PAMP binding activity, and the SNP locus 4407-4408 CG/TA was associated with disease resistance of scallop against L. anguillarum infection, which could be used as a candidate marker for future selection in zhikong scallop breeding program.

Long-term and short-term evolutionary impacts of transposable elements on Drosophila

Transposable elements (TEs) are considered to be genomic parasites and their interactions with their hosts have been likened to the coevolution between host and other nongenomic, horizontally transferred pathogens. TE families, however, are vertically inherited as integral segments of the nuclear genome. This transmission strategy has been suggested to weaken the selective benefits of host alleles repressing the transposition of specific TE variants. On the other hand, the elevated rates of TE transposition and high incidences of deleterious mutations observed during the rare cases of horizontal transfers of TE families between species could create at least a transient process analogous to the influence of horizontally transmitted pathogens. Here, we formally address this analogy, using empirical and theoretical analysis to specify the mechanism of how host-TE interactions may drive the evolution of host genes. We found that host TE-interacting genes actually have more pervasive evidence of adaptive evolution than immunity genes that interact with nongenomic pathogens in Drosophila. Yet, both our theoretical modeling and empirical observations comparing Drosophila melanogaster populations before and after the horizontal transfer of P elements, which invaded D. melanogaster early last century, demonstrated that horizontally transferred TEs have only a limited influence on host TE-interacting genes. We propose that the more prevalent and constant interaction with multiple vertically transmitted TE families may instead be the main force driving the fast evolution of TE-interacting genes, which is fundamentally different from the gene-for-gene interaction of host-pathogen coevolution.

Identity by descent between distant relatives: detection and applications

Short segments of identity by descent (IBD) between individuals with no known relationship can be detected using genome-wide single nucleotide polymorphism data and recently developed statistical methodology. Emerging applications for the detected IBD segments include IBD mapping, haplotype phase inference, genotype imputation, and inference of population structure. In this review, we explain the principles behind methods for IBD segment detection, describe recently developed methods, discuss approaches to comparing methods, and give an overview of applications.

Evolutionary conflicts between viruses and restriction factors shape immunity

Host restriction factors are potent, widely expressed intracellular blocks to viral replication that are an important component of the innate immune response to viral infection. However, viruses have evolved mechanisms that antagonize restriction factors. Through evolutionary pressure for both host survival and virus replication, an evolutionary 'arms race' has developed that drives continuous rounds of selection for beneficial mutations in the genes encoding restriction factors and their viral antagonists. Because viruses can evolve faster than their hosts, the innate immune system of modern-day vertebrates is for the most part optimized to defend against ancient viruses, rather than newer viral threats. Thus, the evolutionary history of restriction factors might, in part, explain why humans are susceptible or resistant to the viruses present in the modern world.

Variation at innate immunity Toll-like receptor genes in a bottlenecked population of a New Zealand robin

Toll-like receptors (TLRs) are an ancient family of genes encoding transmembrane proteins that bind pathogen-specific molecules and initiate both innate and adaptive aspects of the immune response. Our goal was to determine whether these genes show sufficient genetic diversity in a bottlenecked population to be a useful addition or alternative to the more commonly employed major histocompatibility complex (MHC) genotyping in a conservation genetics context. We amplified all known avian TLR genes in a severely bottlenecked population of New Zealand's Stewart Island robin (Petroica australis rakiura), for which reduced microsatellite diversity was previously observed. We genotyped 17-24 birds from a reintroduced island population (including the 12 founders) for nine genes, seven of which were polymorphic. We observed a total of 24 single-nucleotide polymorphisms overall, 15 of which were non-synonymous, representing up to five amino-acid variants at a locus. One locus (TLR1LB) showed evidence of past directional selection. Results also confirmed a passerine duplication of TLR7. The levels of TLR diversity that we observe are sufficient to justify their further use in addressing conservation genetic questions, even in bottlenecked populations.

Evolution of MHC class I genes in the European badger (Meles meles)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune system and provides a good model with which to understand the evolutionary processes underlying functional genes. Trans-species polymorphism and orthology are both commonly found in MHC genes; however, mammalian MHC class I genes tend to cluster by species. Concerted evolution has the potential to homogenize different loci, whereas birth-and-death evolution can lead to the loss of orthologs; both processes result in monophyletic groups within species. Studies investigating the evolution of MHC class I genes have been biased toward a few particular taxa and model species. We present the first study of MHC class I genes in a species from the superfamily Musteloidea. The European badger (Meles meles) exhibits moderate variation in MHC class I sequences when compared to other carnivores. We identified seven putatively functional sequences and nine pseudogenes from genomic (gDNA) and complementary (cDNA) DNA, signifying at least two functional class I loci. We found evidence for separate evolutionary histories of the α1 and α2/α3 domains. In the α1 domain, several sequences from different species were more closely related to each other than to sequences from the same species, resembling orthology or trans-species polymorphism. Balancing selection and probable recombination maintain genetic diversity in the α1 domain, evidenced by the detection of positive selection and a recombination event. By comparison, two recombination breakpoints indicate that the α2/α3 domains have most likely undergone concerted evolution, where recombination has homogenized the α2/α3 domains between genes, leading to species-specific clusters of sequences. Our findings highlight the importance of analyzing MHC domains separately.

Large-scale MHC class II genotyping of a wild lemur population by next generation sequencing

The critical role of major histocompatibility complex (MHC) genes in disease resistance, along with their putative function in sexual selection, reproduction and chemical ecology, make them an important genetic system in evolutionary ecology. Studying selective pressures acting on MHC genes in the wild nevertheless requires population-wide genotyping, which has long been challenging because of their extensive polymorphism. Here, we report on large-scale genotyping of the MHC class II loci of the grey mouse lemur (Microcebus murinus) from a wild population in western Madagascar. The second exons from MHC-DRB and -DQB of 772 and 672 individuals were sequenced, respectively, using a 454 sequencing platform, generating more than 800,000 reads. Sequence analysis, through a stepwise variant validation procedure, allowed reliable typing of more than 600 individuals. The quality of our genotyping was evaluated through three independent methods, namely genotyping the same individuals by both cloning and 454 sequencing, running duplicates, and comparing parent-offspring dyads; each displaying very high accuracy. A total of 61 (including 20 new) and 60 (including 53 new) alleles were detected at DRB and DQB genes, respectively. Both loci were non-duplicated, in tight linkage disequilibrium and in Hardy-Weinberg equilibrium, despite the fact that sequence analysis revealed clear evidence of historical selection. Our results highlight the potential of 454 sequencing technology in attempts to investigate patterns of selection shaping MHC variation in contemporary populations. The power of this approach will nevertheless be conditional upon strict quality control of the genotyping data.

MHC diversity and mate choice in the magellanic penguin, Spheniscus magellanicus

We estimated levels of diversity at the major histocompatibility complex (MHC) class II DRß1 gene in 50 breeding pairs of the Magellanic penguin and compared those to estimates from Humboldt and Galapagos penguins. We tested for positive selection and 2 conditions required for the evolution of MHC-based disassortative mating: 1) greater MHC diversity between breeding pairs compared to random mating, and 2) associations between MHC genotype and fitness. Cloning and sequencing of the DRß1 gene showed that Magellanic penguins had higher levels of genetic variation than Galapagos and Humboldt penguins. Sequence analysis revealed 45 alleles with 3.6% average proportion of nucleotide differences, nucleotide diversity of 0.030, and observed heterozygosity of 0.770. A gene phylogeny showed 9 allelic lineages with interspersed DRß1 sequences from Humboldt and Galapagos penguins, indicating ancestral polymorphisms. d (N)/d (S) ratios revealed evidence for positive selection. Analysis of breeding pairs showed no disassortative mating preferences. Significant MHC genotype/fitness associations in females suggest, however, that selection for pathogen resistance plays a more important role than mate choice in maintaining diversity at the MHC in the Magellanic penguin. The differential effect of MHC heterozygosity on fitness between the sexes is likely associated with the relative role of hatching and fledging rates as reliable indicators of overall fitness in males and females.

Insight into gene evolution within Cervidae and Bovidae through genetic variation in MHC-DQA in the black muntjac (Muntiacus crinifrons)

The critical role that the major histocompatibility complex plays in the immune recognition of parasites and pathogens makes its evolutionary dynamics exceptionally relevant to ecology, population biology, and conservation studies. The black muntjac is a rare deer endemic to a small mountainous region in eastern China. We found that this species has two DQA loci through cDNA expression and sequence variation analysis. The level of variation at both DQA loci was found to be extremely low (three alleles for DQA1 and four alleles for DQA2), possibly because of past bottlenecks and the species' relatively solitary behavior pattern. The ratio of d(N)/d(S) in the putative peptide binding region of the DQA2 locus (13.36, P = 0.012) was significantly larger than one but not that of DQA1 (0.94, P = 0.95), suggesting strong positive selection at the DQA2 but not at the DQA1 locus. This difference might reflect different sets of evolutionary selection pressures acting on the two loci. The phylogenetic tree showed that DQA1 alleles from two species of Cervidae and two of Bovidae grouped together, as did the DQA2 alleles. However, different genes from the four species were located in separate branches. These results lead us to suggest that these DQA alleles are derived from primordial DQA genes from a common ancestor and are maintained in Cervidae and Bovidae since their divergence around 25.5-27.8 million years ago.

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.

The Père David's deer MHC class I genes show unexpected diversity patterns, with monomorphic classical genes but polymorphic nonclassical genes and pseudogenes

Père David's deer (Elaphurus davidianus) is a highly inbred species that arose from 11 founders but now comprises a population of about 3,000 individuals, making it interesting to investigate the adaptive variation of this species from the major histocompatibility complex (MHC) perspective. In this study, we isolated Elda-MHC class I loci using magnetic bead-based cDNA hybridization, and examined the molecular variations of these loci using single-strand conformation polymorphism (SSCP) and sequence analysis. We obtained seven MHC class I genes, which we designated F1, F12, G2, I7, AF, I8, and C1. Our analyses of stop codons, phylogenetic trees, amino acid conservation, and G+C content revealed that F1, F12, G2, and I7 were classical genes, AF was a nonclassical gene, and I8 and C1 were pseudogenes. Our subsequent molecular examinations showed that the diversity pattern in the Père David's deer was unusual. Most mammals have more polymorphic classical class I loci vs. the nonclassical and neutral genes. In contrast, the Père David's deer was found to be monomorphic at classical genes F1, F12, G2, and I7, dimorphic at the nonclassical AF gene, dimorphic at pseudogene I8, and tetramorphic at pseudogene C1. The adverse polymorphism patterns of Elda-I genes might provide evidence for selection too faster deplete MHC variation than drift in the bottlenecked populations, while the postbottleneck tetramorphism of the C1 pseudogene appears to be evidence of strong historical balancing selection.

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.

Polymorphisms at MHC class II DRB1 exon 2 locus in Pyrenean chamois (Rupicapra pyrenaica pyrenaica)

Chamois (Rupicapra spp.) are mountain ungulates from Southern and Central Europe and the Near East. A newly reported border disease virus (BDV) has affected the easternmost populations of Pyrenean chamois, leading to a dramatic population decrease that may drive to genetic variability loss. The Major Histocompatibility Complex (MHC) is a sensitive marker for genetic variation of populations: polymorphism on the MHC genes is affected both by pathogens and population dynamics and it is ecologically relevant, as depending on host-pathogen relationships and life history features. In the present study MHC class II DRB1 exon 2 variation was investigated in 81 Pyrenean chamois (Rupicapra pyrenaica pyrenaica) belonging to four populations. Haplotype analysis, population genetics statistics and network analysis were carried out, in order to analyze variability, phylogeography and genealogy, and the effects of geography and demographic trend. Twenty-nine haplotypes were identified, 26 of them newly described, with high Gene diversity (Gd). The variability observed in the easternmost populations of Pyrenean chamois showed a higher genetic diversity than that previously reported for other populations of Pyrenean and Cantabrian chamois (Rupicapra pyrenaica parva). The most frequent allele was RupyDRB*15, previously undetected, which seems to play a significant role in genotyping the variability, suggesting a possible effect of positive selection.

Analytical methods for immunogenetic population data

In this chapter, we describe analyses commonly applied to immunogenetic population data, along with software tools that are currently available to perform those analyses. Where possible, we focus on tools that have been developed specifically for the analysis of highly polymorphic immunogenetic data. These analytical methods serve both as a means to examine the appropriateness of a dataset for testing a specific hypothesis, as well as a means of testing hypotheses. Rather than treat this chapter as a protocol for analyzing any population dataset, each researcher and analyst should first consider their data, the possible analyses, and any available tools in light of the hypothesis being tested. The extent to which the data and analyses are appropriate to each other should be determined before any analyses are performed.

Major histocompatibility complex class IIA and IIB genes of the spotted halibut Verasper variegatus: genomic structure, molecular polymorphism, and expression analysis

The major histocompatibility complex (MHC) is a large genomic region characterized by extremely high polymorphism and its association with resistance/susceptibility to disease in vertebrates. In this study, the full lengths of MHC IIA and IIB cDNA were obtained from spotted halibut (Verasper variegates) by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The genomic structure, molecular polymorphism, and expression patterns were examined to study MHC II gene functions in fish. As in other teleosts, the genomic structure of the spotted halibut MHC IIA contained 4 exons and 3 introns. The deduced amino acid sequence of the class IIA molecule shared 28-79% similarity with those of teleosts and mammals. Nine class IIA alleles were identified from five individuals. Three alleles originating from a single individual suggested the existence of at least two class IIA loci in the genome. Six exons and 5 introns were identified from spotted halibut MHC IIB, and the deduced amino acid sequence shared 33-79% similarity with those of teleosts and mammals. Twelve alleles were identified, among which five were observed in a single individual, which suggested at least three class IIB loci. Quantitative real-time PCR analysis revealed the presence of class IIA and IIB transcripts in nine normal tissues with high expression level in kidney and gill. Furthermore, MHC IIA and IIB are probably two candidates of immune molecules involved in the acute-phase response in spotted halibut, because their transcriptional levels were significantly up-regulated in blood and liver after bacterial challenge.

Candidate targets of balancing selection in the genome of Staphylococcus aureus

Signatures of balancing selection can highlight polymorphisms and functions that are important to the long-term fitness of a species. We performed a first genome-wide scan for balancing selection in a bacterial species, Staphylococcus aureus, which is a common cause of serious antimicrobial-resistant infections of humans. Using a sliding window approach, the genomes of 16 strains of S. aureus, including 5 new genome sequences presented here, and 1 outgroup strain of S. epidermidis were scanned for signatures of balancing selection. A total of 195 short windows were investigated based on their extreme values of both Tajima's D (>2.03) and π/K ratios (>0.12) relative to the rest of the genome. To test the unusualness of these windows, an Approximate Bayesian Computation framework was used to select a null demographic model that better accounted for the observed data than did the standard neutral model. A total of 186 windows were demonstrated to be unusual under the null model and, thus, represented candidate loci under balancing selection. These 186 candidate windows were located within 99 candidate genes that were spread across 62 different loci. Nearly all the signal (97.2%) was located within coding sequences; balancing selection on gene regulation apparently occurs through the targeting of global regulators such as agr and gra/aps. The agr locus had some of the strongest signatures of balancing selection, which provides new insight into the causes of diversity at this locus. The list of candidate genes included multiple virulence-associated genes and was significantly enriched for functions in amino acid and inorganic ion transport and metabolism and in defense mechanisms against innate immunity and antimicrobials, highlighting these particular functions as important to the fitness of this pathogen.

Heligmosomoides polygyrus infection is associated with lower MHC class II gene expression in Apodemus flavicollis: indication for immune suppression?

Due to their key role in recognizing foreign antigens and triggering the subsequent immune response the genes of the major histocompatibility complex (MHC) provide a potential target for parasites to attack in order to evade detection and expulsion from the host. A diminished MHC gene expression results in less activated T cells and might serve as a gateway for pathogens and parasites. Some parasites are suspected to be immune suppressors and promote co-infections of other parasites even in other parts of the body. In our study we found indications that the gut dwelling nematode Heligmosomoides polygyrus might exert a systemic immunosuppressive effect in yellow-necked mice (Apodemus flavicollis). The amount of hepatic MHC class II DRB gene RNA transcripts in infected mice was negatively associated with infection intensity with H. polygyrus. The hepatic expression of immunosuppressive cytokines, such as transforming growth factor β and interleukin 10 was not associated with H. polygyrus infection. We did not find direct positive associations of H. polygyrus with other helminth species. But the prevalence and infection intensity of the nematodes Syphacia stroma and Trichuris muris were higher in multiple infected individuals. Furthermore, our data indicated antagonistic effects in the helminth community of A. flavicollis as cestode infection correlated negatively with H. polygyrus and helminth species richness. Our study shows that expression analyses of immune relevant genes can also be performed in wildlife, opening new aspects and possibilities for future ecological and evolutionary research.

MHC genotypes associate with resistance to a frog-killing fungus

The emerging amphibian disease chytridiomycosis is caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Amphibian populations and species differ in susceptibility to Bd, yet we know surprisingly little about the genetic basis of this natural variation. MHC loci encode peptides that initiate acquired immunity in vertebrates, making them likely candidates for determining disease susceptibility. However, MHC genes have never been characterized in the context of chytridiomycosis. Here, we performed experimental Bd infections in laboratory-reared frogs collected from five populations that show natural variation in Bd susceptibility. We found that alleles of an expressed MHC class IIB locus associate with survival following Bd infection. Across populations, MHC heterozygosity was a significant predictor of survival. Within populations, MHC heterozygotes and individuals bearing MHC allele Q had a significantly reduced risk of death, and we detected a significant signal of positive selection along the evolutionary lineage leading to allele Q. Our findings demonstrate that immunogenetic variation affects chytridiomycosis survival under controlled experimental conditions, confirming that host genetic polymorphisms contribute to chytridiomycosis resistance.

Characterization and evolution of MHC class II B genes in Ardeid birds

Major histocompatibility complex (MHC) is a multi-gene family that is very suitable to investigate a wide range of open questions in evolutionary ecology. In this study, we characterized two expressed MHC class II B genes (DAB1 and DAB2) in the Grey Heron (Aves: Ardea cinerea). We further developed the primer pairs to amplify and sequence two MHC class II B loci in ten ardeid birds. Phylogenetic analysis revealed that different parts of the genes showed different evolutionary patterns. The exon 2 sequences tended to cluster two gene-specific lineages. In each lineage, exon 2 sequences from several species showed closer relationships than sequences within species, and two shared identical alleles were found between species (Egretta sacra and Nycticorax nycticorax; Egretta garzetta and Bubulcus ibis), supporting the hypothesis of trans-species polymorphism. In contrast, the species-specific intron 2 plus partial exon 3 tree suggested that DAB1 and DAB2 were subject to concerted evolution. GENECONV analyses showed the gene exchange played an important role in the ardeid MHC evolution.

Natural genetic variation caused by small insertions and deletions in the human genome

Human genetic variation is expected to play a central role in personalized medicine. Yet only a fraction of the natural genetic variation that is harbored by humans has been discovered to date. Here we report almost 2 million small insertions and deletions (INDELs) that range from 1 bp to 10,000 bp in length in the genomes of 79 diverse humans. These variants include 819,363 small INDELs that map to human genes. Small INDELs frequently were found in the coding exons of these genes, and several lines of evidence indicate that such variation is a major determinant of human biological diversity. Microarray-based genotyping experiments revealed several interesting observations regarding the population genetics of small INDEL variation. For example, we found that many of our INDELs had high levels of linkage disequilibrium (LD) with both HapMap SNPs and with high-scoring SNPs from genome-wide association studies. Overall, our study indicates that small INDEL variation is likely to be a key factor underlying inherited traits and diseases in humans.

Contrasting responses to selection in class I and class IIα major histocompatibility-linked markers in salmon

Comparison of levels and patterns of genetic variation in natural populations either across loci or against neutral expectation can yield insight into locus-specific differences in the strength and direction of evolutionary forces. We used both approaches to test the hypotheses on patterns of selection on major histocompatibility (MH)-linked markers. We performed temporal analyses of class I and class IIα MH-linked markers and eight microsatellite loci in two Atlantic salmon populations in Ireland on two temporal scales: over six decades and 9 years in the rivers Burrishoole and Delphi, respectively. We also compared contemporary Burrishoole and Delphi samples with nearby populations for the same loci. On comparing patterns of temporal and spatial differentiation among classes of loci, the class IIα MH-linked marker was consistently identified as an outlier compared with patterns at the other microsatellite loci or neutral expectation. We found higher levels of temporal and spatial heterogeneity in heterozygosity (but not in allelic richness) for the class IIα MH-linked marker compared with microsatellites. Tests on both within- and among-population differentiation are consistent with directional selection acting on the class IIα-linked marker in both temporal and spatial comparisons, but only in temporal comparisons for the class I-linked marker. Our results indicate a complex pattern of selection on MH-linked markers in natural populations of Atlantic salmon. These findings highlight the importance of considering selection on MH-linked markers when using these markers for management and conservation purposes.

Inference on the strength of balancing selection for epistatically interacting loci

Existing inference methods for estimating the strength of balancing selection in multi-locus genotypes rely on the assumption that there are no epistatic interactions between loci. Complex systems in which balancing selection is prevalent, such as sets of human immune system genes, are known to contain components that interact epistatically. Therefore, current methods may not produce reliable inference on the strength of selection at these loci. In this paper, we address this problem by presenting statistical methods that can account for epistatic interactions in making inference about balancing selection. A theoretical result due to Fearnhead (2006) is used to build a multi-locus Wright-Fisher model of balancing selection, allowing for epistatic interactions among loci. Antagonistic and synergistic types of interactions are examined. The joint posterior distribution of the selection and mutation parameters is sampled by Markov chain Monte Carlo methods, and the plausibility of models is assessed via Bayes factors. As a component of the inference process, an algorithm to generate multi-locus allele frequencies under balancing selection models with epistasis is also presented. Recent evidence on interactions among a set of human immune system genes is introduced as a motivating biological system for the epistatic model, and data on these genes are used to demonstrate the methods.

Allelic polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of Cynoglossus semilaevis (Cynoglossidae)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and evolutionary mechanism of MHC genes have been investigated in many mammals, but much less is known in teleosts. We examined the polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of the half-smooth tongue sole (Cynoglossus semilaevis); 23 alleles were found in this species. Gene duplication manifested as three to six distinct sequences at each domain in the same individuals. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the PBR domain, suggesting balancing selection for maintaining polymorphisms at the MHC II DAB locus. Many positive selection sites were found to act very intensely on antigen-binding sites of MHC class II DAB gene.

Nonhuman primate models of human immunology

Nonhuman primates have been used for biomedical research for several decades. The high level of genetic homology to humans coupled with their outbred nature has made nonhuman primates invaluable preclinical models. In this review, we summarize recent advances in our understanding of the nonhuman primate immune system, with special emphasis on studies carried out in rhesus macaque (Macaca mulatta). We highlight the utility of nonhuman primates in the characterization of immune senescence and the evaluation of new interventions to slow down the aging of the immune system.

Putative causes and consequences of MHC variation within and between locally adapted stickleback demes

Genes of the major histocompatibility complex (MHC) have been a source of considerable research interest, owing in large part to the growing body of evidence that they may be subject to both natural and sexual selection. However, much remains to be learned about the dynamics of MHC genes in subdivided populations, particularly those characterized by divergent ecological pressures. In this study, we attempt to disentangle the relative roles of both parasite-mediated selection and MHC-mediated mate choice in an open estuarine system inhabited by two parapatric, adaptively divergent threespine stickleback (Gasterosteus aculeatus) demes. We sequenced the putative peptide-binding region (PBR) of an estimated four Class IIβ loci from 127 individuals, identifying 329 sequence variants (276 translated amino acid sequences). Demes differed significantly both in the frequency of MHC alleles and in the communities of helminth parasites infecting resident sticklebacks. Strong signatures of natural selection were inferred from analyses of codon substitutions, particularly in the derived (freshwater) rather than the ancestral (marine) deme. Relationships between parasite load and MHC diversity were indicative of balancing selection, but only within the freshwater deme. Signals of MHC-mediated mate choice were weak and differed significantly between demes. Moreover, MHC-mediated mate choice was significantly influenced by environmental salinity and appeared of secondary importance to tendencies towards assortative mating. We discuss the implications of these findings in respect to ecological adaptation and the potential demographic consequences of possible outcomes of MHC-mediated mate choice.

Parasite-mediated evolution of the functional part of the MHC in primates

The major histocompatibility complex (MHC) is a key model of genetic polymorphism, but the mechanisms underlying its extreme variability are debated. Most hypotheses for MHC diversity focus on pathogen-driven selection and predict that MHC polymorphism evolves under the pressure of a diverse parasite fauna. Several studies reported that certain alleles offer protection against certain parasites, yet it remains unclear whether variation in parasite pressure more generally covaries with allelic diversity and rates of molecular evolution of MHC across species. We tested this prediction in a comparative study of 41 primate species. We characterized polymorphism of the exon 2 of DRB region of the MHC class II. Our phylogenetic analyses controlled for the potential effects of neutral mutation rate, population size, geographic origin and body mass and revealed that nematode species richness associates positively with nonsynonymous nucleotide substitution rate at the functional part of the molecule. We failed to find evidence for allelic diversity being strongly related to parasite species richness. Continental distribution was a strong predictor of both allelic diversity and substitution rate, with higher values in Malagasy and Neotropical primates. These results indicate that parasite pressure can influence the different estimates of MHC polymorphism, whereas geography plays an independent role in the natural history of MHC.

Balancing selection maintains a form of ERAP2 that undergoes nonsense-mediated decay and affects antigen presentation

A remarkable characteristic of the human major histocompatibility complex (MHC) is its extreme genetic diversity, which is maintained by balancing selection. In fact, the MHC complex remains one of the best-known examples of natural selection in humans, with well-established genetic signatures and biological mechanisms for the action of selection. Here, we present genetic and functional evidence that another gene with a fundamental role in MHC class I presentation, endoplasmic reticulum aminopeptidase 2 (ERAP2), has also evolved under balancing selection and contains a variant that affects antigen presentation. Specifically, genetic analyses of six human populations revealed strong and consistent signatures of balancing selection affecting ERAP2. This selection maintains two highly differentiated haplotypes (Haplotype A and Haplotype B), with frequencies 0.44 and 0.56, respectively. We found that ERAP2 expressed from Haplotype B undergoes differential splicing and encodes a truncated protein, leading to nonsense-mediated decay of the mRNA. To investigate the consequences of ERAP2 deficiency on MHC presentation, we correlated surface MHC class I expression with ERAP2 genotypes in primary lymphocytes. Haplotype B homozygotes had lower levels of MHC class I expressed on the surface of B cells, suggesting that naturally occurring ERAP2 deficiency affects MHC presentation and immune response. Interestingly, an ERAP2 paralog, endoplasmic reticulum aminopeptidase 1 (ERAP1), also shows genetic signatures of balancing selection. Together, our findings link the genetic signatures of selection with an effect on splicing and a cellular phenotype. Although the precise selective pressure that maintains polymorphism is unknown, the demonstrated differences between the ERAP2 splice forms provide important insights into the potential mechanism for the action of selection.

It has long been known that the extremely high levels of genetic diversity present in the major histocompatibility locus (MHC) are due to balancing selection, a type of natural selection that maintains advantageous genetic diversity in populations. The MHC encodes for molecules required for a type of antigen presentation that mediates detection of infected and cancerous cells by the immune system; the genetic diversity of the MHC thus ensures an adequate response to the wide variety of pathogens that humans encounter. Here, we show that other genes involved in the same antigen-presentation pathway are also subject to balancing selection in humans. Specifically, we show that balancing selection acts to maintain two forms of the endoplasmic reticulum aminopeptidase 2 gene (ERAP2), which encodes a protein also involved in antigen presentation. Although the two ERAP2 forms are present in a similar frequency (close to 0.5), they are associated with differences with respect to the levels of MHC molecules on the cell surface of immune cells. In summary, our findings show that natural selection maintains variants of ERAP2 that affect immune surveillance; they also establish ERAP2 as one of the few examples of balancing selection in humans where the selected variant, its functional consequences, and its influence in interpersonal diversity are known.

Characterisation and comparative analysis of MHC-DPA1 exon 2 in the owl monkey (Aotus nancymaae)

The Aotus nancymaae (owl monkey) is an important animal model in biomedical research, particularly for the preclinical evaluation of vaccine candidates against Plasmodium falciparum and Plasmodium vivax, which require a precisely typed major histocompatibility complex. The exon 2 from A. nancymaae MHC-DPA1 gene was characterised in order to infer its allelic diversity and evolutionary history. Aona-DPA1 shows no polymorphism and is related to other primate DPA alleles (including Catarrhini and Platyrrhini), constituting an ancient trans-specific and strongly supported lineage with different variability and selective patterns when compared to other primate-MHC-DPA1 lineages. A. nancymaae monkeys have thus a smaller MHC-DP polymorphism than MHC-DQ or MHC-DR.

Genetic diversity at endoplasmic reticulum aminopeptidases is maintained by balancing selection and is associated with natural resistance to HIV-1 infection

Human ERAP1 and ERAP2 encode two endoplasmic reticulum aminopeptidases. These enzymes trim peptides to optimal size for loading onto major histocompatibility complex class I molecules and shape the antigenic repertoire presented to CD8(+) T cells. Therefore, ERAP1 and ERAP2 may be considered potential selection targets and modulators of infection susceptibility. We resequenced two genic regions in ERAP1 and ERAP2 in three HapMap populations. In both cases, we observed high levels of nucleotide variation, an excess of intermediate-frequency alleles, and reduced population genetic differentiation. The genealogy of ERAP1 and ERAP2 haplotypes was split into two major branches with deep coalescence times. These features suggest that long-standing balancing selection has acted on these genes. Analysis of the Lys528Arg (rs30187 in ERAP1) and Asn392Lys (rs2549782 in ERAP2) variants in an Italian population of HIV-1-exposed seronegative (ESN) individuals and a larger number of Italian controls indicated that rs2549782 significantly deviates from Hardy-Weinberg equilibrium (HWE) in ESN but not in controls. Technical errors were excluded and a goodness-of-fit test indicated that a recessive model with only genetic effects adequately explains HWE deviation. The genotype distribution of rs2549782 is significantly different in the two cohorts (P = 0.004), mainly as the result of an over-representation of Lys/Lys genotypes in the ESN sample (P-value for a recessive model: 0.00097). Our data suggest that genetic diversity in ERAP1 and ERAP2 has been maintained by balancing selection and that variants in ERAP2 confer resistance to HIV-1 infection possibly via the presentation of a distinctive peptide repertoire to CD8(+) T cells.

Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.