Domestication does not narrow MHC diversity in Sus scrofa

The Major Histocompatibility Complex (MHC) is a multigene family of outstanding polymorphism. MHC molecules bind antigenic peptides in the peptide-binding region (PBR) that consists of five binding pockets (P). In this study, we compared the genetic diversity of domestic pigs to that of the modern representatives of their wild ancestors, the wild boar, in two MHC loci, the oligomorphic DQA and the polymorphic DRB1. MHC nucleotide polymorphism was compared with the actual functional polymorphism in the PBR and the binding pockets P1, P4, P6, P7, and P9. The analysis of approximately 200 wild boars collected throughout Europe and 120 domestic pigs from four breeds (three pureblood, Pietrain, Leicoma, and Landrace, and one mixed Danbred) revealed that wild boars and domestic pigs share the same levels of nucleotide and amino acid polymorphism, allelic richness, and heterozygosity. Domestication did not appear to act as a bottleneck that would narrow MHC diversity. Although the pattern of polymorphism was uniform between the two loci, the magnitude of polymorphism was different. For both loci, most of the polymorphism was located in the PBR region and the presence of positive selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous substitutions in the PBR. P4 and P6 were the most polymorphic binding pockets. Functional polymorphism, i.e., the number and the distribution of pocket variants within and among populations, was significantly narrower than genetic polymorphism, indicative of a hierarchical action of selection pressures on MHC loci.     

Major histocompatibility complex variation at class II DQA locus in the brown hare (Lepus europaeus)

The major histocompatability complex (MHC) is a multigene family of receptors that bind and present antigenic peptides to T-cells. Genes of the MHC are characterized by an outstanding genetic polymorphism, which is considered to be maintained by positive selection. Sites involved in peptide binding form binding pockets (P) that are collectively termed the peptide-binding region (PBR). In this study, we examined the level of MHC genetic diversity within and among natural populations of brown hare ( Lepus europaeus ) from Europe and Anatolia choosing for analysis of the second exon of the DQA locus, one of the most polymorphic class II loci. We aimed at an integrated population genetic analysis of L. europeaus by (i) correlating MHC polymorphism to genetic variability and phylogenetic status estimated previously from maternally (mtDNA) and biparentally (allozymes, microsatellites) inherited loci; and (ii) comparing full-length exon amino acid polymorphism with functional polymorphism in the PBR and the binding pockets P1, P6 and P9. A substantial level of DQA exon 2 polymorphism was detected with two completely different set of alleles between the Anatolian and European populations. However, the phylogeny of full-length exon 2 Leeu-DQA alleles did not show a strong phylogeographic signal. The presence of balancing selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous in the PBR and a trans-species pattern of evolution detected after phylogenetic reconstruction. The differentiating patterns detected between genetic and functional polymorphism, i.e. the number and the distribution of pocket variants within and among populations, indicated a hierarchical action of selection pressures.     

Major histocompatibility complex class I loci from the gray whale (Eschrichtius robustus)

Sequences from exons encoding the peptide binding region of MHC class I (MHC-I) molecules were isolated from California gray whale (Eschrichtius robustus) genomic DNA to initiate an investigation of variation in these genes in a cetacean. These represent the first mysticete MHC-I sequences to be reported. The analysis of gray whale MHC-I sequences suggests the presence of at least three loci, which share greatest similarity to MHC-I in the ungulates, consistent with current views on cetacean phylogenetics. The peptide binding region of MHC is the most polymorphic part of the molecule and analysis of the variation and synonymous to nonsynonymous substitution ratios in gray whale sequences found these genes to display polymorphism characteristics similar to that attributed to selection in other species.     

Circumsporozoite Protein Genes of Malaria Parasites (Plasmodium Spp.): Evidence for Positive Selection on Immunogenic Regions

The circumsporozoite (CS) protein is a cell surface protein of the sporozoite, the stage of the life cycle of malaria parasites (Plasmodium spp.) that infects the vertebrate host. Analysis of DNA sequences supports the hypothesis that in Plasmodium falciparum, positive Darwinian selection favors diversity in the T-cell epitopes (peptides presented to T cells by host MHC molecules) of the CS protein. In gene regions encoding T cell epitopes of P. falciparum, the rate of nonsynonymous nucleotide substitution is significantly higher than that of synonymous substitution, whereas this is not true of other gene regions. Furthermore nonsynonymous nucleotide substitutions in these regions cause a change of amino acid residue charge significantly more frequently than expected by chance. By contrast, in Plasmodium cynomolgi, the same regions show no evidence of positive selection, and residue charge is conserved. The CS protein has a central repeat region, which is the target of host antibodies. In P. falciparum, the amino acid sequence of the repeat region is conserved within and between alleles. In P. cynomolgi, on the other hand, there is evidence that positive selection has favored evolution of two different repeat types within a given allele.     

Polymorphism and signature of selection in the MHC class I genes of the three-spined stickleback Gasterosteus aculeatus

The role and intensity of positive selection maintaining the polymorphism of major histocompatibility complex (MHC) class I genes in the three-spined stickleback Gasterosteus aculeatus was investigated. The highly polymorphic set of MHC class I genes found was organized in a single linkage group. Between 5 and 14 sequence variants per individual were identified by single-stranded conformation polymorphism (SSCP) analysis. Segregation analysis studied in 10 three-spined stickleback families followed the expected pattern of Mendelian inheritance. The gamete fusion in three-spined stickleback thus seems to be random with respect to the MHC class I genes. The DNA sequence analyses showed that the expressed MHC class I loci are under strong selection pressure, possibly mediated by parasites. Codons that were revealed to be under positive selection are potentially important in antigen binding. MHC class I sequences did not form significant supported clusters within a phylogenetic tree. Analogous to MHC class II genes, it was not possible to assign the class I sequences to a specific locus, suggesting that the class I genes may have been generated by recent gene duplication.     

Characterization of MHC class I and II genes in a subantarctic seabird,the blue petrel, <Emphasis Type="Italic">Halobaena caerulea</Emphasis> (Procellariiformes)

The great polymorphism observed in the major histocompatibility complex (MHC) genes is thought to be maintained by pathogen-mediated selection possibly combined with MHC-disassortative mating, guided by MHC-determined olfactory cues. Here, we partly characterize the MHC class I and II B of the blue petrel, Halobaena caerulea (Procellariiformes), a bird with significant olfactory abilities that lives under presumably low pathogen burdens in Subantarctica. Blue petrels are long-lived, monogamous birds which suggest the necessity of an accurate mate choice process. The species is ancestral to songbirds (Passeriformes; many MHC loci), although not to gamefowls (Galliformes; few MHC loci). Considering the phylogenetic relationships and the low subantarctic pathogen burden, we expected few rather than many MHC loci in the blue petrel. However, when we analysed partial MHC class I and class II B cDNA and gDNA sequences we found evidence for as many as at least eight MHC class I loci and at least two class II B loci. These class I and II B sequences showed classical MHC characteristics, e.g. high nucleotide diversity, especially in putative peptide-binding regions where signatures of positive selection was detected. Trans-species polymorphism was found between MHC class II B sequences of the blue petrel and those of thin-billed prion, Pachyptila belcheri, two species that diverged ∼25 MYA. The observed MHC allele richness in the blue petrel may well serve as a basis for mate choice, especially since olfactory discrimination of MHC types may be possible in this species.     

Evolution of a new nonclassical MHC class I locus in two Old World primate species

 HLA-G is a nonclassical major histocompatibility complex (MHC) class I molecule that is expressed only in the human placenta, suggesting that it plays an important role at the fetal-maternal interface. In rhesus monkeys, which have similar placentation to humans, the HLA-G orthologue is a pseudogene. However, rhesus monkeys express a novel placental MHC class I molecule, Mamu-AG, which has HLA-G-like characteristics. Phylogenetic analysis of AG alleles in two Old World primate species, the baboon and the rhesus macaque, revealed limited diversity characteristic of a nonclassical MHC class I locus. Gene trees constructed using classical and nonclassical primate MHC class I alleles demonstrated that the AG locus was most closely related to the classical A locus. Interestingly, gene tree analyses suggested that the AG alleles were most closely related to a subset of A alleles which are the products of an ancestral interlocus recombination event between the A and B loci. Calculation of the rates of synonymous and nonsynonymous substitution at the AG locus revealed that positive selection was not acting on the codons encoding the peptide binding region. In exon 4, however, the rate of nonsynonymous substitution was significantly lower than the rate of synonymous substitution, suggesting that negative selection was acting on these codons. Received: 22 April 1998 / Revised: 15 July 1998     

Recurrent positive selection at bgcn, a key determinant of germ line differentiation, does not appear to be driven by simple coevolution with its partner protein bam

Surveys of nucleotide sequence polymorphism in Drosophila melanogaster and Drosophila simulans were performed at 2 interacting loci crucial for gametogenesis: bag-of-marbles (bam) and benign gonial cell neoplasm (bgcn). At the polymorphism level, both loci appear to be evolving under the expectations of the neutral theory. However, ratios of polymorphism and divergence for synonymous and nonsynonymous mutations depart significantly from neutral expectations for both loci consistent with a previous observation of positive selection at bam. The deviations suggest either an excess of synonymous polymorphisms or an excess of nonsynonymous fixations at both loci. Synonymous evolution appears to conform to neutrality at bam. At bgcn, there is evidence of positive selection affecting preferred synonymous mutations along the D. simulans lineage. However, there is also a significantly higher rate of nonsynonymous fixations at bgcn within D. simulans. Thus, the deviation from neutrality detected by the McDonald-Kreitman test at these 2 loci is likely due to the selective acceleration of nonsynonymous fixations. Differences in the pattern of amino acid fixations between these 2 interacting proteins suggest that the detected positive selection is not due to a simple model of coevolution.     

Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber)

Major histocompatibility complex (MHC) genes, coding molecules which play an important role in immune response, are the most polymorphic genes known in vertebrates. However, MHC polymorphism in some species is limited. MHC monomorphism at several MHC class I and II loci was previously reported for two neighbouring northern European populations of the Eurasian beaver (Castor fiber) and reduced selection for polymorphism has been hypothesized. Here, we analysed a partial sequence of the second exon of the MHC II DRB locus from seven relict European and Asian beaver populations. We detected 10 unique alleles among 76 beavers analysed. Only a western Siberian population was polymorphic, with four alleles detected in 10 individuals. Each of the remaining populations was fixed for a different allele. Sequences showed considerable divergence, suggesting the long persistence of allelic lineages. A significant excess of nonsynonymous substitutions was detected at the antigen binding sites, indicating that sequence evolution of beaver DRB was driven by positive selection. Current MHC monomorphism in the majority of populations may be the result of the superimposition of the recent bottleneck on pre-existing genetic structure resulting from population subdivision and differential pathogen pressure.     

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.     

Codon usage bias and base composition in MHC genes in humans and common chimpanzees

 Codon bias and base composition in major histocompatibility complex (MHC) sequences have been studied for both class I and II loci in Homo sapiens and Pan troglodytes. There is low to moderate codon bias for the MHC of humans and chimpanzees. In the class I loci, the same level of moderate codon bias is seen for HLA-B, HLA-C, Patr-A, Patr-B, and Patr-C, while at HLA-A the level of codon bias is lower. There is a correlation between codon usage bias and G+C content in the A and B loci in humans and chimps, but not at the C locus. To examine the effect of diversifying selection on codon bias, we subdivided class I alleles into antigen recognition site (ARS) and non-ARS codons. ARS codons had lower bias than non-ARS codons. This may indicate that the constraint of codon bias on nucleotide substitution may be selected against in ARS codons. At the class II loci, there are distinct differences between alpha and beta chain genes with respect to codon usage, with the beta chain genes being much more biased. Species-specific differences in base composition were seen in exon 2 at the DRB1 locus, with lower GC content in chimpanzees. Considering the complex evolutionary history of MHC genes, the study of codon usage patterns provides us with a better understanding of both the evolutionary history of these genes and the evolution of synonymous codon usage in genes under natural selection. Received: 2 April 1998 / Revised: 2 September 1998     

Computer simulation analysis suggests weak balancing selection operative at the MICA locus

A high degree of polymorphism has been reported at the major histocompatibility class I chain-related gene A (MICA) locus, which is located 46 kb away from HLA-Bin the human major histocompatibility complex (MHC) class I region. Although it is known that the polymorphisms at the conventional MHC class I loci have been maintained by balancing selection, it is unclear whether positive natural selection is also operative in maintaining the polymorphism at the MICA locus. In order to explain the degree of polymorphism at the MICA locus, a computer simulation study was carried out. The high degree of polymorphism at the MICA locus (heterozygosity and number of polymorphic residues) could not be explained solely by balancing selection at the HLA-B locus even if no recombination was assumed between MICA and HLA-B. Although there is no definite evidence indicating that balancing selection is operative at the MICA locus, our results suggest that the MICA gene is subject to weak balancing selection.     

MHC class I genes in a New World primate, the cotton-top tamarin (Saguinus oedipus), have evolved by an active process of loci turnover

 Lymphocytes of a New World primate, the cotton-top tamarin (Saguinus oedipus), express classical G–related major histocompatibility complex (MHC) class I molecules with unusually limited polymorphism and variability. Three G-related loci, an F locus, an E locus, and two pseudogenes (So-N1 and So-N3) have been identified by cDNA library screening and extensive PCR analysis of both cDNA and genomic DNA from the cotton-top tamarin. Furthermore, each genus of the subfamily Callitrichinae (tamarins and marmosets) appears to express its own unique set of MHC class I genes, likely due to a rapid turnover of loci. The rapid emergence of unique MHC class I genes in the Callitrichinae genera, resulting from an active process of duplication and inactivation of loci, may account for the limited diversity of the MHC class I genes in the cotton-top tamarin. To determine the nature of the entire complement of MHC class I genes in the cotton-top tamarin, we synthesized a genomic DNA library and screened it with MHC class I-specific probes. We isolated nine new MHC class I pseudogenes from this library. These newly isolated tamarin G–related MHC class I pseudogenes are not closely related to any of their functional counterparts in the tamarin, suggesting that they do not share a recent common ancestral gene with the tamarin's currently expressed MHC class I loci. In addition, these tamarin sequences display a high rate of nonsynonymous substitutions in their putative peptide binding region. This indicates that the genes from which they have derived were likely subject to positive selection and, therefore, were once functional. Our data support the notion that an extremely high rate of loci turnover is largely responsible for the limited diversity of the MHC class I genes in the cotton-top tamarin. Received: 15 September 1997 / Revised: 2 July 1998     

Trans-species polymorphism and evidence of selection on class II MHC loci in tuco-tucos (Rodentia: Ctenomyidae)

Balancing selection acting over the evolutionary history of a lineage can result in the retention of alleles among species for longer than expected under neutral evolution. The associated pattern of trans-species polymorphism, in which similar or even identical alleles are shared among species, is often used to infer that balancing selection has occurred. The genes of the major histocompatibility complex (MHC) are thought to be subject to balancing selection that maintains alleles associated with response to specific pathogens. To explore the role of balancing selection in shaping MHC diversity in ctenomyid rodents, we examined allelic variability at the class II DRB and DQA loci in 18 species in the genus Ctenomys. Previous studies of four of these species had revealed significant within-population evidence of positive selection on MHC loci. The current study expands upon these analyses to (1) evaluate among-species evidence of positive selection and (2) explore the potential for balancing selection on MHC genes. Interspecific nucleotide sequence variation revealed significant evidence of positive selection on the DRB and DQA loci. At the same time, comparisons of phylogenetic trees for these MHC loci with a putative species tree based on mitochondrial sequence data revealed multiple examples of trans-specific polymorphism, including sharing of identical DRB and DQA alleles among distantly related species of Ctenomys. These findings suggest that MHC genes in these animals have historically been subject to balancing selection and yield new insights into the complex suite of forces shaping MHC diversity in free-living vertebrates.     

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake,the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long‐term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.     

Patterns of nucleotide substitutions inferred from the phylogenies of the class I major histocompatibility complex genes

Summary Patterns of nucleotide substitutions in human major histocompatibility complex (MHC) class I genes were estimated by using phylogenetic trees of DNA sequences. The pattern is defined as a set of 12 parameters, each of which represents the relative frequency of substitutions from a particular nucleotide to another. The pattern at the antigen recognition sites (ARS) in functional MHC genes was remarkably different from that at the remaining coding region (non-ARS). In particular, the proportion of transitions among all the nucleotide substitutions (P _s) was extremely low at the third codon positions of ARS. In the HLA-A genes, P _s at the third codon positions was only 6% in ARS, whereas it was 69% in non-ARS. In HLA-B, the corresponding values were 30% in ARS and 80% in non-ARS, respectively. On the other hand, P _s in a class I pseudogene (HLA-H) was 57%, which was in good agreement with P _s in other pseudogenes. Because pseudogenes are selectively neutral, the pattern in pseudogenes is regarded as the pattern of spontaneous substitution mutations. In general, the pattern in functional genes that are subject to selective forces deviates from the pattern in pseudogenes. At the third codon positions in coding regions, transitions scarcely cause amino acid replacements, whereas about half of transversions do cause replacements. Accordingly, P _s at the third codon positions decreases if amino acid replacements are accelerated by natural selection but increases if amino acids are conserved by functional constraint. Our observations imply that the ARS region is subject to natural selection favoring amino acid replacements, whereas the non-ARS region is subject to functional constraint. Offprint requests to: T. Gojobori     

Adaptive evolution in the rat olfactory receptor gene family

Summary Comparison of DNA sequences of the rat (Rattus norvegicus) olfactory receptor gene family revealed an unusual pattern of nucleotide substitution in the gene region encoding the second extracellular domain (E2) of the protein. In this domain, nonsynonymous nucleotide differences between members of this subfamily that caused a change in amino acid residue polarity were over four times more frequent than nonsynonymous differences that did not cause a polarity change. This nonrandom pattern of nucleotide substitution is evidence of past directional selection favoring diversification of the E2 domain among members of this subfamily. This in turn suggests that E2 may play some important role in the functions unique to each member of the olfactory receptor family, and that it may perhaps be an odorant binding domain.Offprint requests to: A.L. Hughes     

Maintenance of DQB1 polymorphisms in primates.

To understand the evolution of the class II major histocompatibility complex (MHC) DQB1 locus in primates, the second exons of seven DQB1 alleles from five non-human primate species were amplified by polymerase chain reaction. Comparisons of these and other primate sequences show that no between-species diversity is greater than within-species diversity, suggesting maintenance of DQB1 alleles through the history of Old-World primates. There is a preponderance of nonsynonymous nucleotide substitutions at antigen-binding-site codons; this pattern is in marked contrast to what is seen at the closely related, presumably nonfunctional DQB2 gene. The results support the hypothesis that DQB1 polymorphism is maintained by overdominant selection relating to antigen presentation.     

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.     

Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline Mhc-DRB peptide-binding region

The long-term evolution of major histocompatibility complex (MHC) involves the birth-and-death process and independent divergence of loci during episodes punctuated by natural selection. Here, we investigated the molecular signatures of natural selection at exon-2 of MHC class II DRB gene which includes a part of the peptide-binding region (PBR) in seven of eight putative extant Felidae lineages. The DRB alleles in felids can be mainly divided into five lineages. Signatures of trans-species polymorphism among major allelic lineages indicate that balancing selection has maintained the MHC polymorphism for a long evolutionary time. Analysis based on maximum likelihood models of codon substitution revealed overall purifying selection acting on the feline DRB. Sites that have undergone positive selection and those that are under divergent selective pressure among lineages were detected and found to fall within the putative PBR. This study increased our understanding of the nature of selective forces acting on DRB during feline radiation.