Allelic genealogy under overdominant and frequency-dependent selection and polymorphism of major histocompatibility complex loci

Maintenance of multiallelic polymorphism at the MHC region

Models that purport to explain the maintenance of MHC polymorphism must be able to explain a variety of phenomena. (1) The range of MHC allele frequencies at some of the loci is very large, with some alleles quite common and many others rare, while at others the range of allele frequencies is far narrower. (2) MHC alleles and their frequencies often have long persistence times, in some cases tens of millions of years. (3) Random-mating populations appear to be in Hardy-Weinberg equilibrium for MHC. (4) There is no obvious, strong and consistent selection pressure yet detected that acts differentially on different MHC genotypes. (5) Because the allelic composition of the MHC polymorphism does change over evolutionary time, the MHC system must be capable of accommodating new alleles with similar properties without destruction of the equilibria that permit the maintenance of the older alleles. In this review I examined the degree to which a large number of models that have been proposed fit these criteria. These include heterosis, marginal overdominance, conditional heterosis, assortative mating, maternal-fetal incompatibility, molecular mimicry, minority advantage, pathogen adaptation, and optimum allele frequency models. Most of the models do poorly at accounting for a number of the above phenomena. The last class, optimum allele frequency models, have the most satisfactory set of properties. However, optimum allele frequency models require mechanisms that somehow "feed back" from the frequency of an allele in the population to the fitness of an organism carrying that allele. Thus, these models require that MHC polymorphisms be maintained by some type of group selection. Evidence for an against optimum allele frequency selection, and ways in which this type of selection might be detected experimentally, are presented.

Excess nonsynonymous substitution of shared polymorphic sites among self-incompatibility alleles of Solanaceae

The function of the self-incompatibility locus (S locus) of many plant species dictates that natural selection will favor high levels of protein diversity. Pairwise sequence comparisons between S alleles from four species of Solanaceae reveal remarkably high sequence diversity and evidence for shared polymorphism. The level of amino acid constraint was found to be significantly heterogeneous among different regions of the gene, with some regions being highly constrained and others appearing to be virtually unconstrained. In some regions of the protein, there was an excess of nonsynonymous over synonymous substitution, consistent with the strong diversifying selection that must operate on this locus. These hypervariable regions are candidates for the sites that determine functional allelic identity. Simple contingency table tests show that sites that have polymorphism shared between species have more nonsynonymous substitution than polymorphic sites that do not exhibit shared polymorphism. This is consistent with the idea that adaptive evolution favoring amino acid replacement is occurring at sites with shared polymorphism. Tests of clustered polymorphism reveal that an unusually low rate of recombination must be occurring in this locus, allowing very ancient alleles to preserve their identity.

Origins of H-2 polymorphism in the house mouse. II. Characterization of a model population and evidence for heterozygous advantage

Comparison of the rate of synonymous and nonsynonymous nucleotide substitutions suggests that certain regions of the functional H-2 genes, which are part of the mouse major histocompatibility complex (Mhc), are under strong positive selection pressure. Thus far, however, little evidence has been provided for the existence of such pressure in natural mouse populations. We have, therefore, initiated experiments designed to test the hypothesis of positive selection acting on H-2 loci. The experiments are being carried out on two natural mouse populations in Jerusalem, Israel. One population occupies a space of about 100 m2 in a chicken coop, the other lives in a nearby field in which "mouse stations" providing food and shelter have been set up. Extensive typing of these two populations revealed the presence of only four H-2 haplotypes. Mice in the two populations breed continually all year around, yet population size varies seasonally, with population maxima in winter and minima in summer. The population in the chicken coop contains a relatively stable nucleus which may be organized in demes with an excess of females over males and limited territorial mobility. The rest of the mice stay in the population for a short time only and then either die or emigrate. The field population is smaller and more loosely organized than the chicken-coop population, with demes probably forming only during population maxima. For the rest of the time breeding in this population is probably panmictic. At a population minimum in the summer of 1984, H-2 homozygotes happened to predominate over heterozygotes. This situation, however, lasted for a short time only and thereafter there was a continuous, statistically highly significant increase in the proportion of H-2 heterozygotes of one or two types. The increase occurred in both populations but was more apparent in the chicken-coop population. This observation provides the first experimental evidence that heterozygous advantage might be one of the mechanisms maintaining high H-2 polymorphism in natural populations of the house mouse.

Common west African HLA antigens are associated with protection from severe malaria

A large case-control study of malaria in West African children shows that a human leucocyte class I antigen (HLA-Bw53) and an HLA class II haplotype (DRB1*1302-DQB1*0501), common in West Africans but rare in other racial groups, are independently associated with protection from severe malaria. In this population they account for as great a reduction in disease incidence as the sickle-cell haemoglobin variant. These data support the hypothesis that the extraordinary polymorphism of major histocompatibility complex genes has evolved primarily through natural selection by infectious pathogens.

Mating patterns in seminatural populations of mice influenced by MHC genotype

Because of the central role of major histocompatibility complex (MHC) genes in immune recognition, it is often assumed that parasite-driven selection maintains the unprecendented genetic diversity of these genes. But associations between MHC genotype and specific infectious diseases have been difficult to identify with a few exceptions such as Marek's disease and malaria. Alternatively, MHC-related reproductive mechanisms such as selective abortion and mating preferences could be responsible for the diversity. To determine both the nature and strength of selection operating on MHC genes by we have studied components of selection in seminatural populations of mice (Mus musculus domesticus). Here we assess MHC-related patterns of reproduction and early (preweaning) mortality by analysing 1,139 progeny born in nine populations, and 662 progeny from laboratory matings. Reproductive mechanisms, primarily mating preferences, result in 27% fewer MHC-homozygous offspring than expected from random mating. MHC genotype had no detectable influence on neonatal (preweaning) mortality. These mating preferences are strong enough to account for most of the MHC genetic diversity found in natural populations of Mus.

Role of diversifying selection and gene conversion in evolution of major histocompatibility complex loci

Genes at the major histocompatibility complex (MHC) in mammals are known to have exceptionally high polymorphism and linkage disequilibrium. In addition, these genes form highly complicated gene families that have evolved through gene conversion and unequal crossing-over. It has been shown recently that amino acid substitution at the antigen recognition site (ARS) is more rapid than synonymous substitution, suggesting some kind of positive natural selection working at the ARS. It is highly desirable to know the interactive effect of gene conversion and natural selection on the evolution and variation of MHC gene families. A population genetic model is constructed that incorporates both selection and gene conversion. Diversifying selection is assumed in which sequence diversity is enhanced not only between alleles at the same locus but also between duplicated genes. Expressed and nonexpressed loci are assumed as in the class I gene family of MHC, with gene conversion occurring among all loci. Extensive simulation studies reveal that very weak selection at individual amino acid sites in combination with gene conversion can explain the unusual pattern of evolution and polymorphisms. Here both gene conversion and natural selection contribute to enhancing polymorphism.

HLA population genetics

The HLA system has been extensively studied from an evolutionary perspective. The region contains a number of closely linked genes whose products control a variety of functions concerned with the regulation of immune responses. In addition, the genetic predisposition to over 40 diseases maps to this region. A number of observations indicate that strong selection is acting on the HLA region, including its extensive polymorphism with very even allele frequencies, the preferential occurrence of high levels of variability at positions critical to antigen recognition, the great age of alleles and the patterns of linkage disequilibrium among loci. The form of the selection is unknown. Although balancing selection is a strong candidate, it seems unlikely that only one selective mechanism is operating in this complex multigene family region. Mutation, recombination and gene conversion all contribute to the generation of HLA variability. The apparent great age of many HLA alleles revealed by phylogenetic analysis suggests that the absolute rate of production of new variants is not high. Detailed studies of population and evolutionary features of the HLA region are necessary for an informed discussion of the evolution of disease predisposing genes and epitopes, and of complex multigene families.

Statistical models of the overdispersed molecular clock

The most commonly used statistical model to describe the rate constancy of molecular evolution (molecular clock) is a simple Poisson process in which the variance of the number of amino acid or nucleotide substitutions in a particular gene should be equal to the mean and henceforth the dispersion index, the ratio of the variance to the mean, should be equal to one. Recent sequence data, however, have shown that the substitutional process in molecular evolution is often considerably overdispersed and have called into question the generality of using a simple Poisson process. Several efforts have been made to develop more realistic models of molecular evolution. In this paper, I will show that the spatial (site-specific) variation in the rate of molecular evolution is an improbable cause of the overdispersion and then review various statistical models which take the temporal variation into account. Although these models do not immediately specify what the mechanisms of molecular evolution might be, they do make qualitatively different predictions and give some insight into their inference. One way to distinguish them is suggested. In addition, effects of selected substitutions that presumably occur after a major change in a molecule are quasi-quantitatively examined. It is most likely that the overdispersion of molecular clock is due either to a major molecular reconfiguration (fluctuating neutral space) led by a series of subliminal neutral changes or to selected substitutions fine-tuning a molecule after a major molecular change. Although the latter possibility, of course, violates the simplest neutrality assumption, it would not impair the neutral theory as a whole.

Evolution of MHC polymorphism: extensive sharing of polymorphic sequence motifs between human and bovine DRB alleles

The evolution of MHC polymorphism has been studied by comparing the amino acid and nucleotide sequences of 14 bovine and 32 human DRB alleles. The comparison revealed an extensive sharing of polymorphic sequence motifs in the two species. Almost identical sets of residues were found at several highly polymorphic amino acid positions in the putative antigen recognition site. Consequently, certain bovine alleles were found to be more similar to certain human alleles than to other bovine alleles. In contrast, the frequencies of silent nucleotide substitutions were found to be much higher in comparisons between species than within species implying that none of the human or bovine DRB alleles originated before the divergence of these distantly related species. The results suggest that the observed similarity in DRB polymorphism is due to convergent evolution and possibly the sharing of short ancestral sequence motifs. However, the relative role of the latter mechanism is difficult to assess due to the biased base composition in the first domain exon of polymorphic class II beta genes. The frequency of silent substitutions between DRB alleles was markedly lower in cattle than in man suggesting that the DRB diversity has evolved more rapidly in the former species.

Estimation of gene diversity at the b locus of the constant region of the immunoglobulin light chain in natural populations of European rabbit (Oryctolagus cuniculus) in Portugal, Andalusia and on the Azorean Islands

The minimal gene diversity at a locus of the antibody constant region, as estimated in natural populations of rabbit, revealed levels of heterozygosity similar to those reported for the major histocompatibility complex in human and murine populations. Sera of 416 wild rabbits were collected on the Iberian peninsula and on three islands of the Azorean archipelago and analyzed for the occurrence of the serological markers of the b locus of the immunoglobulin light chain. All four serotypes present in domestic rabbits were found in Portugal. They represented less than 50% of the gene pool. In Andalusia this was less than 15% and on the Azorean islands less than 10%. The pronounced and systematic hierarchy in allele frequencies, previously found in populations from the more recent distribution area of the species, was not observed. On the peninsula, the frequencies of the "domestic" alleles were similar, averaging 10%. The Portuguese sample revealed a total heterozygosity of at least 87%. This high value was supported by at least 11 serologically different alleles, none of them occurring at frequencies above 20%. These data are in agreement with an Iberian origin of the European rabbit and strongly suggest the coalescence of b locus allelic lines drawn from Iberian and western populations. The role of balancing selection in the evolution of the b locus polymorphism was further emphasized.

Overdispersed molecular clock at the major histocompatibility complex loci

The extent of amino acid differences of major histocompatibility complex molecules within species is unusually high, consistent with the finding that some pairs of alleles have persisted for more than ten million years and the view that the polymorphism has been maintained by natural selection. The disparity between synonymous and non-synonymous substitutions in the antigen recognition site, however, suggests that some non-synonymous sites have undergone a number of substitutions whereas others have little or none. To describe statistically such an overdispersed underlying process, commonly used Poisson processes are inadequate. An alternative process leads to the surprising conclusion that each non-synonymous site has accumulated as many as 2.6 substitutions, on the average, in the two lineages leading to humans and mice. The standard deviation is also very large (6.6) and the dispersion index (the ratio of the variance to the mean) is at least 17. The substitution process thus inferred qualitatively agrees with the disposition (a boomerang pattern) of substitutions between HLA-A2 and Aw68 alleles, and quantitatively agrees well with that expected where the evolution of major histocompatibility complex molecules has long been driven mostly by balancing selection.

Polymorphism at the self-incompatibility locus in Solanaceae predates speciation

Sequences of 11 alleles of the gametophytic self-incompatibility locus (S locus) from three species of the Solanaceae family have recently been determined. Pairwise comparisons of these alleles reveal two unexpected observations: (i) amino acid sequence similarity can be as low as 40% within species and (ii) some interspecific similarities are higher than intraspecific similarities. The gene genealogy clearly illustrates this unusual pattern of relationships. The data suggest that some of the polymorphism at the S locus existed prior to the divergence of these species and has been maintained to the present. In support of this hypothesis, the number of shared polymorphic sites was found to exceed the number found in simulations with independent accumulation of mutations. Strictly neutral evolution is exceedingly unlikely to maintain the polymorphism for such a long time. The allele multiplicity and extreme age of the alleles is consistent with Wright's classic one-locus population genetic model of gametophytic self-incompatibility. Similarities between the plant S locus and the mammalian major histocompatibility complex are discussed.

Alternative fitness models with the same allele frequency dynamics

For any set of one- or two-locus genotypic fitnesses there are alternative sets, usually frequency-dependent and often with quite different biological meanings, that give rise to the same equations for change of allele or haplotype frequencies. Therefore, it is not possible to distinguish among alternative fitness models from allele or haplotype frequency trajectories or equilibrium distributions. For a single locus and for two loci when linkage equilibrium can be assumed, a simple procedure generates some of the alternative fitness sets.