Nucleotide substitution at major histocompatibility complex class II loci: evidence for overdominant selection

SNP profile within the human major histocompatibility complex reveals an extreme and interrupted level of nucleotide diversity

The human major histocompatibility complex (MHC) is characterized by polymorphic multicopy gene families, such as HLA and MIC (PERB11); duplications; insertions and deletions (indels); and uneven rates of recombination. Polymorphisms at the antigen recognition sites of the HLA class I and II genes and at associated neutral sites have been attributed to balancing selection and a hitchhiking effect, respectively. We, and others, have previously shown that nucleotide diversity between MHC haplotypes at non-HLA sites is unusually high (>10%) and up to several times greater than elsewhere in the genome (0.08%-0.2%). We report here the most extensive analysis of nucleotide diversity within a continuous sequence in the genome. We constructed a single nucleotide polymorphism (SNP) profile that reveals a pattern of extreme but interrupted levels of nucleotide diversity by comparing a continuous sequence within haplotypes in three genomic subregions of the MHC. A comparison of several haplotypes within one of the genomic subregions containing the HLA-B and -C loci suggests that positive selection is operating over the whole subgenomic region, including HLA and non-HLA genes. [The sequence data for the multiple haplotype comparisons within the class I region have been submitted to DDBJ/EMBL/GenBank under accession nos. AF029061, AF029062, and AB031005-AB031010. Additional sequence data have been submitted to the DDBJ data library under accession nos. AB031005-AB03101 and AF029061-AF029062.]

Evolutionary diversification of multigene families: allelic selection of toxins in predatory cone snails

In order to investigate the evolution of conotoxin multigene families among two closely related vermivorous CONUS: species, we sequenced 104 four-loop conotoxin mRNAs from two individuals of CONUS: ebraeus and compared these with sequences already obtained from CONUS: abbreviatus. In contrast to the diversity of conotoxin sequences obtained from C. abbreviatus, only two common sequence variants were recovered from C. ebraeus. Segregation patterns of the variants in these two individuals and restriction digests of four-loop conotoxin amplification products from nine additional individuals suggest that the common variants are alleles from a single locus. These two putative alleles differ at nine positions that occur nonrandomly in the toxin-coding region of the sequences. Moreover, all substitutions are at nonsynonymous sites and are responsible for seven amino acid differences among the predicted amino acid sequences of the alleles. These results imply that conotoxin diversity is driven by strong diversifying selection and some form of frequency-dependent or overdominant selection at conotoxin loci, and they suggest that diverse conotoxin multigene families can originate from duplications at polymorphic loci. Furthermore, none of the sequences recovered from C. ebraeus appeared to be orthologs of loci from C. abbreviatus, and attempts to amplify orthologous sequences with locus-specific primers were unsuccessful among these species. These patterns suggest that venoms of closely related CONUS: species may differ due to the differential expression of conotoxin loci.

A 39-kb sequence around a blackbird Mhc class II gene: ghost of selection past and songbird genome architecture

To gain an understanding of the evolution and genomic context of avian major histocompatibility complex (Mhc) genes, we sequenced a 38.8-kb Mhc-bearing cosmid insert from a red-winged blackbird (Agelaius phoeniceus). The DNA sequence, the longest yet retrieved from a bird other than a chicken, provides a detailed view of the process of gene duplication, divergence, and degeneration ("birth and death") in the avian Mhc, as well as a glimpse into major noncoding features of a songbird genome. The peptide-binding region (PBR) of the single Mhc class II B gene in this region, Agph-DAB2, is almost devoid of polymorphism, and a still-segregating single-base-pair deletion and other features suggest that it is nonfunctional. Agph-DAB2 is estimated to have diverged about 40 MYA from a previously characterized and highly polymorphic blackbird Mhc gene, Aph-DAB1, and is therefore younger than most mammalian Mhc paralogs and arose relatively late in avian evolution. Despite its nonfunctionality, Agph-DAB2 shows very high levels of nonsynonymous divergence from Agph-DAB1 and from reconstructed ancestral sequences in antigen-binding PBR codons-a strong indication of a period of adaptive divergence preceding loss of function. We also found that the region sequenced contains very few other unambiguous genes, a partial Mhc- class II gene fragment, and a paucity of simple-sequence and other repeats. Thus, this sequence exhibits some of the genomic streamlining expected for avian as compared with mammalian genomes, but is not as densely packed with functional genes as is the chicken Mhc.

Do infectious diseases drive MHC diversity?

The primary function of the major histocompatibility complex (MHC) is to allow the immune system to identify infectious pathogens and eliminate them. Infectious diseases are now thought to be the main selection force that drives and maintains the extraordinary diversity of the MHC.

Sojourn times and substitution rate at overdominant and linked neutral loci

The sojourn times until fixation of an overdominant allele were investigated based on the diffusion equation. Furthermore, the rate of accumulation of mutations, or the substitution rate, was predicted from the mean extinction time of a common overdominant allele. The substitution rate calculated theoretically agreed well with that determined by computer simulation. Overdominant selection enhances the polymorphism at linked loci, while its effect on the sojourn times and the substitution rate at linked loci has not been studied yet. To solve these problems, a model that assumed two linked loci, each with infinite alleles, was examined by computer simulation. A decrease in the recombination rate between two loci markedly changed the distribution of sojourn times of a neutral allele. Although overdominant selection obviously increased the sojourn times and the polymorphism at a linked locus, the rate of nucleotide substitution at the neutral locus was not influenced significantly even if complete linkage was assumed. These results suggest that, in regions containing overdominant genes, linked neutral loci will exhibit elevated levels of polymorphism, but their rate of molecular evolution remains that predicted by neutral theory.

Evolutionary lineages of RT1.Ba in the Australian Rattus

In this study, the evolutionary history of the variable second exon of RT1.Ba and its adjoining intron b are compared across a number of species and subspecies of the Australian RATTUS: Three lineages are identified in the second intron across a range of Rattus species. Two of these lineages, separated by the insertion of a probable rodent short interspersed nucleotide element and by point mutations outside the indel region, are both found in each of the major clades of the endemic Australian RATTUS: This pattern of ancestral polymorphism is reflected in the adjoining exon 2 sequences, although phylogenetic constraints confirm that the clustering is not identical to that of the associated intron sequences. In addition, the coding sequences show evidence of the retention of ancestral polymorphism, with identical exon sequences found in two divergent species, and some indication of gene conversion detected for the exon sequences.

Alu elements support independent origin of prosimian, platyrrhine, and catarrhine Mhc-DRB genes

The primate major histocompatibility complex (Mhc) genes fall into two classes and each of the classes into several families. Of the class II families, the DRB family has a long and complex evolutionary history marked by gene turnover, rearrangement, and molecular convergence. Because the history is not easily decipherable from sequences alone, Alu element insertions were used as cladistic markers to support the surmised phylogenetic relationships among the DRB genes. Intron 1 segments of 24 DRB genes from five platyrrhine species and five DRB genes from three prosimian species were amplified by PCR and cloned, and the amplification products were sequenced or PCR-typed for Alu repeats. Three Alu elements were identified in the platyrrhine and four in the prosimian DRB genes. One of the platyrrhine elements (Alu50J) is also found in the Catarrhini, whereas the other two (Alu62Sc, Alu63Sc) are restricted to the New World monkeys. Similarly, the four prosimian elements are found only in this taxon. This distribution of Alu elements is consistent with the phylogeny of the DRB genes as determined from their intron 1 sequences in an earlier and the present study. It contradicts the exon 2-based phylogeny and thus corroborates the conclusion that the evolution of DRB exon 2 sequences is, to some extent, shaped by molecular convergence. Taken together, the data indicate that each of the assemblages of DRB genes in prosimians, platyrrhines, and catarrhines is derived from a separate ancestral gene.

Primitive synteny of vertebrate major histocompatibility complex class I and class II genes

Major histocompatibility complex (MHC) class I and class II molecules bind to and display peptidic antigens acquired from pathogens that are recognized by lymphocytes coordinating and executing adaptive immune responses. The two classes of MHC proteins have nearly identical tertiary structures and were derived from a common ancestor that probably existed not long before the emergence of the cartilaginous fish. Class I and class II genes are genetically linked in tetrapods but are not syntenic in teleost fish, a phylogenetic taxon derived from the oldest vertebrate ancestor examined to date. Cartilaginous fish (sharks, skates, and rays) are in the oldest taxon of extant jawed vertebrates; we have carried out segregation analyses in two families of nurse sharks and one family of the banded houndshark that revealed a close linkage of class IIalpha and beta genes both with each other and with the classical class I (class Ia) gene. These results strongly suggest that the primordial duplication giving rise to classical class I and class II occurred in cis, and the close linkage between these two classes of genes has been maintained for at least 460 million years in representatives of most vertebrate taxa.

Gene structure of mouse BIT/SHPS-1

BIT/SHPS-1/SIRPalpha/P84 is a unique molecule with a high degree of homology with immune antigen recognition molecules (immunoglobulin, T-cell receptor and MHC), and is highly expressed in the brain. The extracellular region contains three immunoglobulin-like domains (V-type, C1-type and C1-type), and the intracellular region contains two signalling motifs that interact with SHP-2 protein tyrosine phosphatase. BIT-coated plates support cell-substrate adhesion and neurite extension of neurons, and BIT participates in neuronal signal transduction. Diversity of the V-type domain sequences of human BIT has been reported. In the present study we analysed the structure of the mouse BIT gene (Bit). The protein coding region consists of eight exons corresponding to a signal peptide, a V-type domain, a C1-type domain, a C1-type domain, a transmembrane region and three parts of one cytoplasmic region. The two signalling motifs are encoded in one exon. Four splicing forms of mouse BIT were revealed. We also found the sequence diversity in three mouse strains, namely BALB/c, 129/Sv and C57BL/6. The substitution patterns of amino acids and nucleotides indicate positive pressure to alter the amino acids in the V-type domain in evolution. Immunoblot analyses showed that mouse BIT and human BITalpha are predominantly expressed in the brain. On the bases of these findings we discuss the possibility that BIT contributes to the genetic individuality and diversity of the brain.

Trans-speciation maintenance in the MHC region of a polymorphism which includes a polymorphic dinucleotide locus, and the de novo arisal of a polymorphic tetranucleotide microsatellite

Alleles and the surrounding regions of DQCAR, a dinucleotide repeat tightly linked to HLA-DQB1, were sequenced in a range of primate species including man. Three polymorphic regions can usefully be defined in the description of these sequences: the dinucleotide GT repeat itself, the anonymous region 5' of this repeat, and a variable CTGT repeat in the 3' region. The 5' sequence displayed six alleles in the individuals studied. One of these alleles was invariably associated with substitutions in the GT repeat and absence of the CTGT repeat, the others with pure, polymorphic GT repeats and variation in the numbers of CTGT repeats. Haplotypes can be classified by the allele in the 5' region. Those carrying allele 1 were only found in man, those with allele 2 in man, chimpanzee and gorilla. The third haplotype (indicated by the presence of allele 3) was found in chimpanzee, gorilla and orang-utan, the fourth in chimpanzee and gibbon, the fifth in baboon, guenon and mangabey and the sixth in guenon and macaque. The alleles in the 5' region, but from different species, are thus often more similar than alleles from the same species, a phenomenon already shown for some HLA genes. This suggests that major histocompatibility sequences and surrounding sequences shared a correlated evolutionary history. The new polymorphic tetranucleotide microsatellite (CTGT, 3rd region) has possibly arisen de novo from the pre-existing dinucleotide GT. This study provides information not only on the molecular evolution of this particular microsatellite but also of the trans-speciation maintenance of polymorphism of its surrounding sequences.

Variation on islands: major histocompatibility complex (Mhc) polymorphism in populations of the Australian bush rat

Loss of genetic variation in small, isolated populations is commonly observed at neutral or nearly neutral loci. In this study, the loss of genetic variation was assessed in island populations for a locus of major histocompatibility complex (Mhc), a locus shown to be under the influence of balancing selection. A total of 36 alleles was found at the second exon of RT1.Ba in 14 island and two mainland populations of Rattus fuscipes greyii. Despite this high overall diversity, a substantial lack of variation was observed in the small island populations, with 13 islands supporting only one to two alleles. Two populations, Waldegrave and Williams Islands, showed moderately high levels of heterozygosity (52-56%) which were greater than expected under neutrality, suggesting the action of balancing selection. However, congruence between the level of variation at this Mhc locus and in previous allozyme electrophoresis and mitochondrial DNA studies highlights the dominant influence of genetic drift and population factors, such as bottlenecks and structuring in the founding population, in the loss of genetic variation in these small, isolated populations.

Evolution of HLA-class I compared to HLA-class II polymorphism in Terena, a South-American Indian tribe

We have studied the HLA alleles of 60 unrelated healthy Terena and 10 Terena families. They are members of an isolated Brazilian tribe located in Mato Grosso do Sul (South Central Brazil). Six novel alleles were found in this population: HLA-A*0219 (gf = 0.02), A*0222 (gf = 0.15), HLA-B* 3520 (gf = 0.01), B*3521 (gf = 0.03), B*3912 (gf = 0.03) and B*4803 (gf = 0.16). Five of the six novel alleles differ from their putative progenitors by amino acid replacements in residues that contribute to the pockets of the peptide-binding site. Many of the variants defined by molecular methods were not identified correctly by serological typing. We calculated heterozygosity values (H) for HLA-A, -B, -C, DRB1, DQB1 and DPB . The highest values were observed at the HLA-B locus, followed by HLA-A, -DRB1 and DQB1. Residue positions 9, 24, 45, 62, 67, 95, 114, 116, 156, and 163 of HLA class I showed heterozygosity values greater than 0.50. Nine of them contribute to the peptide-binding specificity pockets and one to the T cell receptor binding site. If HLA antigens are useful for defense against pathogenic agents, heterozygosity would offer an advantage by allowing binding of a larger repertoire of peptides to the class I molecules. Individuals that are heterozygous at these positions would probably have a wider repertoire of peptide presentation to T cells. The observed results including the presence of novel alleles in the class I HLA loci suggest a functionally significant, more rapid evolution of class I compared to class II loci in this South American isolated population.

Fifty-million-year-old polymorphism at an immunoglobulin variable region gene locus in the rabbit evolutionary lineage

The Ig heavy chain variable region (V(H)) genes encode the antigen-binding regions of antibodies. The rabbit genome contains more than 100 V(H) genes, but only one (V(H)1) is preferentially used in the VDJ gene rearrangement. Three highly divergent alleles occur at this V(H)1 locus in most rabbit populations. These three V(H) alleles are also present in snowshoe hare populations, indicating that the polymorphism of the V(H)1 alleles is trans-specific. Here we report the results of a phylogenetic analysis of rabbit Ig germ-line V(H) genes (alleles) together with V(H) genes from humans and mice. We have found that all rabbit V(H) genes belong to one mammalian V(H) group (group C), which also includes various human and mouse V(H) genes. Using the rate of nucleotide substitution obtained from human and mouse V(H) sequences, we have estimated that the V(H)1 polymorphism in the rabbit lineage has been maintained for about 50 million years. This extremely long persistence of V(H)1 polymorphism is apparently caused by overdominant selection, though the real mechanism is unclear.

Major histocompatibility complex variation in the endangered Przewalski's horse

The major histocompatibility complex (MHC) is a fundamental part of the vertebrate immune system, and the high variability in many MHC genes is thought to play an essential role in recognition of parasites. The Przewalski's horse is extinct in the wild and all the living individuals descend from 13 founders, most of whom were captured around the turn of the century. One of the primary genetic concerns in endangered species is whether they have ample adaptive variation to respond to novel selective factors. In examining 14 Przewalski's horses that are broadly representative of the living animals, we found six different class II DRB major histocompatibility sequences. The sequences showed extensive nonsynonymous variation, concentrated in the putative antigen-binding sites, and little synonymous variation. Individuals had from two to four sequences as determined by single-stranded conformation polymorphism (SSCP) analysis. On the basis of the SSCP data, phylogenetic analysis of the nucleotide sequences, and segregation in a family group, we conclude that four of these sequences are from one gene (although one sequence codes for a nonfunctional allele because it contains a stop codon) and two other sequences are from another gene. The position of the stop codon is at the same amino-acid position as in a closely related sequence from the domestic horse. Because other organisms have extensive variation at homologous loci, the Przewalski's horse may have quite low variation in this important adaptive region.

Allelic and haplotype variation of major histocompatibility complex class II DRB1 and DQB loci in the St Lawrence beluga (Delphinapterus leucas)

In order to assess levels of major histocompatibility complex (Mhc) variation within the St Lawrence beluga (Delphinapterus leucas) the variation at the beluga Mhc DRB1 class II locus was assessed by single-strand conformation polymorphism (SSCP) analysis of the peptide-binding region for 313 whales collected from 13 sampling locations across North America. In addition, samples from west Greenland and the St Lawrence were also typed at the DQB locus, allowing comparison to a previous study and assessment of linkage disequilibrium of alleles at the two loci. Comparisons of DRB1 and DQB allele frequencies among all sampling locations indicated genetic structure (alpha < 0.005). Most of this structure resulted from differences between the different wintering groups. Significant genetic structure (alpha = 0.05) exists among each pair of the following groups at both the DRB1 and DQB loci; St Lawrence, Hudson Strait, Bering Sea, Cunningham Inlet, and Davis Strait (minus Cunningham Inlet), except the St Lawrence and Hudson Strait for the DQB locus. In the St Lawrence population, six of the eight DRB1 alleles are present representing all five known allelic lineages. Evidence of linkage disequilibrium between the DRB1 and DQB is present in two sampling locations, the St Lawrence and Nuussuaq (alpha = 0.05). Analysis of probable DRB1-DQB haplotypes among groups of beluga suggests a haplotype reduction in the St Lawrence.

Predictability of mutant sequences. Relationships between mutational mechanisms and mutant specificity

Spontaneous mutations are rare and are produced by multiple biochemical mechanisms. Nonetheless, studies of these mechanisms have revealed striking examples in which mutational specificity can be regularly related to a characteristic of the surrounding DNA sequence and/or the enzymes participating in mutagenesis. Thus, to an increasing degree the DNA sequences of mutants are "predictable." This report considers some examples of predictable sequence changes, evidence for their contribution to mutagenesis in populations, and how the predictability of mutant sequences may be useful to improve our interpretation of the molecular course of evolution from DNA sequence comparisons.

Evolution of HLA class II molecules: Allelic and amino acid site variability across populations

Analysis of the highly polymorphic beta1 domains of the HLA class II molecules encoded by the DRB1, DQB1, and DPB1 loci reveals contrasting levels of diversity at the allele and amino acid site levels. Statistics of allele frequency distributions, based on Watterson's homozygosity statistic F, reveal distinct evolutionary patterns for these loci in ethnically diverse samples (26 populations for DQB1 and DRB1 and 14 for DPB1). When examined over all populations, the DQB1 locus allelic variation exhibits striking balanced polymorphism (P < 10(-4)), DRB1 shows some evidence of balancing selection (P < 0.06), and while there is overall very little evidence for selection of DPB1 allele frequencies, there is a trend in the direction of balancing selection (P < 0.08). In contrast, at the amino acid level all three loci show strong evidence of balancing selection at some sites. Averaged over polymorphic amino acid sites, DQB1 and DPB1 show similar deviation from neutrality expectations, and both exhibit more balanced polymorphic amino acid sites than DRB1. Across ethnic groups, polymorphisms at many codons show evidence for balancing selection, yet data consistent with directional selection were observed at other codons. Both antigen-binding pocket- and non-pocket-forming amino acid sites show overall deviation from neutrality for all three loci. Only in the case of DRB1 was there a significant difference between pocket- and non-pocket-forming amino acid sites. Our findings indicate that balancing selection at the MHC occurs at the level of polymorphic amino acid residues, and that in many cases this selection is consistent across populations.

Evolution and biological characterization of human immunodeficiency virus type 1 subtype E gp120 V3 sequences following horizontal and vertical virus transmission in a single family

It has been suggested that immune-pressure-mediated positive selection operates to maintain the antigenic polymorphism on the third variable (V3) loop of the gp120 of human immunodeficiency virus type 1 (HIV-1). Here we present evidence, on the basis of sequencing 147 independently cloned env C2/V3 segments from a single family (father, mother, and their child), that the intensity of positive selection is related to the V3 lineage. Phylogenetic analysis and amino acid comparison of env C2/V3 and gag p17/24 regions indicated that a single HIV-1 subtype E source had infected the family. The analyses of unique env C2/V3 clones revealed that two V3 lineage groups had evolved in the parents. Group 1 was maintained with low variation in all three family members regardless of the clinical state or the length of infection, whereas group 2 was only present in symptomatic individuals and was more positively charged and diverse than group 1. Only virus isolates carrying the group 2 V3 sequences infected and induced syncytia in MT2 cells, a transformed CD4(+)-T-cell line. A statistically significant excess of nonsynonymous substitutions versus synonymous substitutions was demonstrated only for the group 2 V3 region. The data suggest that HIV-1 variants, possessing the more homogeneous group 1 V3 element and exhibiting the non-syncytium-inducing phenotype, persist in infected individuals independent of clinical status and appear to be more resistant to positive selection pressure.

Identification of antigenic escape variants in an immunodominant epitope of hepatitis C virus

Numerous investigators have postulated that one mechanism by which hepatitis C virus (HCV) may evade the immune system is through the formation of escape mutants. This hypothesis is based largely on the observed mutability of the viral genome resulting in evolution of diverse quasispecies over the course of infection. That such diversification is a product of viral RNA polymerase infidelity, immune-driven selection or a combination of the two processes has not been addressed. We have examined sequence variability in a specific segment of HCV RNA encoding a known immunodominant region of the viral helicase, amino acids 358-375 of the non-structural 3 protein. Using sequence-specific oligonucleotide probe hybridization and automated DNA sequencing, we report a high frequency of mutations, essentially all of which result in amino acid replacements. To assess the biological impact of such mutations, corresponding chemically synthesized peptides were compared to wild-type peptide in T cell proliferation assays. We observed that a sizeable fraction of such peptides stimulated attenuated or negligible levels of proliferation by peripheral T cells from a chronically infected patient. This observation is consistent with expectations for immune-mediated selection of escape variants at the epitope level. We postulate that such a mechanism may be important in the immunopathogenesis of HCV infections.

Mutations in Immunodominant T Cell Epitopes Derived from the Nonstructural 3 Protein of Hepatitis C Virus Have the Potential for Generating Escape Variants That May Have Important Consequences for T Cell Recognition

One of the most disturbing features of hepatitis C virus (HCV) is its long-term persistence in the host. One hypothesis to explain this phenomenon is that HCV escapes immune recognition through its intrinsic hypermutability. To determine whether immunodominant T cell epitopes derived from HCV nonstructural 3 (NS3) protein might be subject to sequence variations leading to escape mutants, we examined sequence variations of one IL-2-producing epitope, NS3358–375, and one IL-10-producing epitope, NS3505–521. By PCR amplification, cloning, and sequencing, we observed significant sequence variations in the two epitopes, although the selection intensity for each epitope was different. For NS3358–375, more variants were observed, and for NS3505–521, fewer mutations were observed. Moreover, functional studies revealed that three NS3358–375 and one NS3505–521 variants failed to stimulate T cell proliferation, and two other NS3358–375 and NS3505–521 variants weakly stimulated T cell responses. Our results are consistent with immune selection of viral variants at the epitope level, which may enable HCV to evade host defenses over time.

Natural selection at major histocompatibility complex loci of vertebrates

The loci of the vertebrate major histocompatibility complex encode cell-surface glycoproteins that present peptides to T cells. Certain of these loci are highly polymorphic, and the mechanisms responsible for this polymorphism have been intensely debated. Four independent lines of evidence support the hypothesis that MHC polymorphisms are selectively maintained: (a) The distribution of allelic frequencies does not fit the neutral expectation. (b) The rate of nonsynonymous nucleotide substitution significantly exceeds the rate of synonymous substitution in the codons encoding the peptide-binding region of the molecule. (c) Polymorphisms have been maintained for long periods of time ("trans-species polymorphism"). (d) Introns have been homogenized relative to exons over evolutionary time, suggesting that balancing selection acts to maintain diversity in the latter, in contrast to the former.

Effect of gene conversion on polymorphic patterns at major histocompatibility complex loci

The pattern of polymorphisms at MHC loci was studied by computer simulations and by DNA sequence analyses, with special reference to synonymous and non-synonymous divergence among alleles. A model of overdominance plus short-term selection was simulated for a gene family with allelic and non-allelic gene conversion. It was found that the ratio of non-synonymous-to-synonymous divergence among alleles decreases as the non-allelic conversion becomes more frequent. On the other hand, the ratio of synonymous divergence at the antigen recognition site to that at non-coding region increases by non-allelic conversion. By comparing such a pattern with the observed values of mammalian genes, it is suggested that non-allelic conversion is less frequent in human genes than in genes of other mammals. A quantitative analysis on the patchwork pattern was attempted by using identity excess among nucleotide sites of a gene. Comparison of simulation results with DNA data analysis shows that the rate of allelic conversion is high.

Evolution of the mammalian MHC: natural selection, recombination, and convergent evolution

The genes that encode molecules involved in antigen presentation within the class I and class II regions of the mammalian major histocompatibility complex (MHC) include several that are highly polymorphic. There is evidence that this polymorphism is maintained by positive selection, most likely overdominant selection, relating to their role in presenting foreign peptides to T cells. This selection can maintain allelic lineages for much longer periods of time than neutral polymorphisms are expected to last, but sharing of polymorphic amino acid motifs among species of different mammalian orders is due to independent (or convergent) evolution rather than common ancestry. It has been suggested that interallelic recombination (gene conversion) plays a role in enhancing polymorphism, but there is evidence of striking differences among loci with respect to the rate at which such recombination has contributed to current polymorphism. Recent attempts to interpret linkage relationships in the MHC region as evidence of ancient genomic duplications are not supported by phylogenetic analysis. Rather, natural selection may have played a role in the linkage of other genes to those of the MHC.

Monomorphism and polymorphism at Mhc DRB loci in domestic and wild ruminants

Genetic polymorphism at Mhc class II DRB loci was investigated in samples of musk-ox from Canada and Greenland; moose from Sweden, Norway, Canada, and Alaska; roe deer from Norway and Sweden; reindeer from Svalbard and Norway; fallow deer from Norway and Sweden; and red deer from Norway. The results were compared with published data on cattle, bison, goat, sheep, and red deer. Cattle-specific primers amplified a single DRB locus in all species except fallow deer and red deer, in which two loci were found. Single strand conformation polymorphism analysis and DNA sequence analysis were employed to detect genetic polymorphism. Complete monomorphism was found in musk-ox and fallow deer. Limited polymorphism was found in the moose, roe deer, and reindeer from Svalbard, whereas intermediate to extensive DRB diversity was present in reindeer from Norway and in bison, sheep, goat, cattle, and red deer. The restricted Mhc diversity in moose, roe deer, and fallow deer is notable in relation to the dramatic population expansion of moose and roe deer in Sweden during this century and since fallow deer is used for meat and game production with good results and without any marked disease problems. The results question the view that species or populations with restricted Mhc diversity have poor resistance to infectious diseases. A phylogenetic tree analysis revealed a clustering of DRB sequences within species rather than within allelic lineages across species. The results suggest trans-species persistence of polymorphic sequence motifs rather than of allelic lineages.

Comparative genome organization of the major histocompatibility complex: lessons from the Felidae

The mammalian major histocompatibility complex (MHC) has taught both immunologists and evolutionary biologists a great deal about the patterns and processes that have led to immune defenses. Driven principally by human and mouse studies, comparative MHC projects among other mammalian species offer certain advantages in connecting MHC genome characters to natural situations. We have studied the MHC in the domestic cat and in several wild species of Felidae. Our observations affirm class I and class II homology with other mammalian orders, derivative gene duplications during the Felidae radiation, abundant persistent trans-species allele polymorphism, recombination-derived amino acid motifs, and inverted ratios of non-synonymous to silent substitutions in the MHC peptide-binding regions, consistent with overdominant selection in class I and II genes. MHC diversity as quantified in population studies is a powerful barometer of historic demographic reduction for several endangered species including cheetahs, Asiatic lions, Florida panthers and tigers. In two cases (Florida panther and cheetah), reduced MHC variation may be contributing to uniform population sensitivity to emerging infectious pathogens. The Felidae species, nearly all endangered and monitored for conservation concerns, have allowed a glimpse of species adaptation, mediated by MHC divergence, using comparative inferences drawn from human and mouse models.

The mouse major histocompatibility complex: some assembly required

We have assembled a contig of 81 yeast artificial chromosome clones that spans 8 Mb and contains the entire major histocompatibility complex (Mhc) from mouse strain C57BL/6 (H2b), and we are in the process of assembling an Mhc contig of bacterial artificial chromosome (BAC) clones from strain 129 (H2bc), which differs from C57BL/6 in the H2-Q and H2-T regions. The current BAC contig extends from Tapasin to D17Leh89 with gaps in the class II, H2-Q, and distal H2-M regions. Only four BAC clones were required to link the class I genes of the H2-Q and H2-T regions, and no new class I gene was found in the previous gap. The proximal 1 Mb of the H2-M region has been analyzed in detail and is ready for sequencing; it includes 21 class I genes or fragments, at least 14 olfactory receptor-like genes, and a number of non-class I genes that clearly establish a conserved synteny with the class I regions of the human and rat Mhc.

Molecular clock and recombination in primate Mhc genes

To set an accurate chronological framework to the evolution of primate class I and II genes in the major histocompatibility complex (Mhc), the rate of silent nucleotide substitutions in exons and introns is examined for various cDNA and genome sequences currently available. The rate is sensitive to the GC content and correlates negatively with increased GC biases at the third codon positions of Mhc genes. The intergenic recombination rate in the HLA region is estimated from the synonymous nucleotide differences at 37 linked loci. Any HLA subregion is recombined more or less at the ordinary rate of 1 cM per 1 Mb, although the rate may be reduced in some subregions. This information is used to discuss HLA haplotypes when they are applied to studies of human demography. The unusual polymorphism in the alpha-helix of HLA-DRB1 is also revisited in relation to intragenic recombination, but the molecular mechanism and the evolutionary cause both remain enigmatic.

Recombination within the human MHC

Recombination (crossing over) in the human MHC is thought to have played a role in generation of novel alleles at various HLA loci. It is also responsible for the diversity observed at the haplotype level, although the functional consequences of this activity are not clear. Historic and family studies of recombination have provided estimations of recombination fractions across the MHC and identified potential hotspots for recombination in the class II region. Other characteristics of recombination in the human MHC such as haplotype specificity in recombination frequency and localized sequence motifs involved in recombination have been considered, but have been difficult to address given the constraints of human population studies. Single-sperm typing holds promise in overcoming some of the limitations inherent in the study of recombination in human populations. Both family-based and sperm typing analyses of recombination, along with our knowledge of linkage disequilibrium patterns in the MHC, may provide novel information regarding the evolution of HLA haplotypes that will be difficult to obtain by other means.

Identification of regions in which positive selection may operate in S-RNase of Rosaceae: implication for S-allele-specific recognition sites in S-RNase

A stylar S-RNase is associated with gametophytic self-incompatibility in the Rosaceae, Solanaceae, and Scrophulariaceae. This S-RNase is responsible for S-allele-specific recognition in the self-incompatible reaction, but how it functions in specific discrimination is not clear. Window analysis of the numbers of synonymous (dS) and non-synonymous (dN) substitutions in rosaceous S-RNases detected four regions with an excess of dN over dS in which positive selection may operate (PS regions). The topology of the secondary structure of the S-RNases predicted by the PHD method is very similar to that of fungal RNase Rh whose tertiary structure is known. When the sequences of S-RNases are aligned with the sequence of RNase Rh based on the predicted secondary structures, the four PS regions correspond to two surface sites on the tertiary structure of RNase Rh. These findings suggest that in S-RNases the PS regions also form two sites and are candidates for the recognition sites for S-allele-specific discrimination.

What sharks can tell us about the evolution of MHC genes

Similarity in structural features would argue that sharks possess class I, class IIA and class IIB genes, coding for classical peptide-presenting molecules, as well as non-classical class I genes. Some aspects of shark major histocompatibility complex genes are similar to teleost genes and others are similar to tetrapod genes. Shark class I genes form a monophyletic group, as also seen for tetrapods, but the classical and nonclassical genes form two orthologous clades, as seen for teleosts. Teleost class I genes arose independently at least four different times with the nonclassical genes of ray-finned fishes and all the shark and lobe-finned fish class I genes forming 1 clade. The ray-finned fish classical class I genes arose separately. In phylogenetic trees of class II alpha 2 and beta 2 domains, the shark and tetrapod genes cluster more closely than the teleost genes and, unlike the teleost sequences, the class II alpha 1 domains of sharks and tetrapods lack cysteines. On the other hand, both shark and teleost genes display sequence motifs in the antigen-binding cleft that have persisted over very long time periods. The similarities may reflect common selective pressures on species in aqueous environments while differences may be due to different evolutionary rates.

Major histocompatibility genes in cyprinid fishes: theory and practice

The first teleostean MHC sequences were described for carp. Subsequent studies in a number of cyprinid fishes showed that the class I sequences of these fishes are of particular interest. Two distinct lineages (Cyca-Z and Cyca-U) are found in the common and ginbuna crucian carp, but only the U lineage is present in zebrafish and other non-cyprinid species. The presence of the Z lineage is hypothesised to be the result of an allotetraploidisation event. Both phylogenetic analyses and amino acid sequence comparisons suggest that Cyca-Z sequences are non-classical class I sequences, probably similar to CD1. The comprehensive phylogenetic analyses of these sequences revealed different phylogenetic histories of the exons encoding the extracellular domains. The MHC genes were studied in laboratory and natural models. The natural model addressed the evolution of MHC genes in a Barbus species flock. Sequence analysis of class I and class II supported the species designation of the morphotypes present in the lake, and as a consequence the trans-species hypothesis of MHC polymorphism. The laboratory model involves the generation of gynogenetic clones, which can be divergently selected for traits such as high and low antibody response. The role of MHC molecules can be investigated further by producing a panel of isogenic lines.

Evidence for balancing selection operating at the het-c heterokaryon incompatibility locus in a group of filamentous fungi

In filamentous fungi, het loci (for heterokaryon incompatibility) are believed to regulate self/nonself-recognition during vegetative growth. As filamentous fungi grow, hyphal fusion occurs within an individual colony to form a network. Hyphal fusion can occur also between different individuals to form a heterokaryon, in which genetically distinct nuclei occupy a common cytoplasm. However, heterokaryotic cells are viable only if the individuals involved have identical alleles at all het loci. One het locus, het-c, has been characterized at the molecular level in Neurospora crassa and encodes a glycine-rich protein. In an effort to understand the role of this locus in filamentous fungi, we chose to study its evolution by analyzing het-c sequence variability in species within Neurospora and related genera. We determined that the het-c locus was polymorphic in a field population of N. crassa with close to equal frequency of each of the three allelic types. Different species and even genera within the Sordariaceae shared het-c polymorphisms, indicating that these polymorphisms originated in an ancestral species. Finally, an analysis of the het-c specificity region shows a high occurrence of nonsynonymous substitution. The persistence of allelic lineages, the nearly equal allelic distribution within populations, and the high frequency of nonsynonymous substitutions in the het-c specificity region suggest that balancing selection has operated to maintain allelic diversity at het-c. Het-c shares this particular evolutionary characteristic of departing from neutrality with other self/nonself-recognition systems such as major histocompatibility complex loci in mammals and the S (self-incompatibility) locus in angiosperms.

Point mutation of a rubella virus E1 protein T-cell epitope by substitution of single amino acid reversed the restrictive HLA-DR polymorphism: a possible mechanism maintaining HLA polymorphism

The influence of single amino acid substitutions within a rubella E1 protein T-cell epitope, E1(273-284) on T-cell recognition was studied. Substitutions of an uncharged amino acid A for an E or for a T and substitution of a T for S were found to not significantly reduce the T-cell responses. However, substitution of a charged residue such as E for hydrophobic residues (I, V, or W); D for Q; or a relatively larger size amino acid for polar residues completely abolished the cytotoxicities mediated by E1(273-284)-specific T-cell clone. A set of single amino acid-substituted peptide analogs of E1(273-284) not eliciting cytotoxicity of the T-cell clone was used to test the influence of point mutation of the epitope on HLA DR restrictions. A panel of B-cell lines with different DR4 subtypes was used as targets in cytotoxicity assays to determine the restrictive HLA molecules. Results showed that modification of the T-cell epitope by point mutation could reverse the HLA DR restriction from one allele to other alleles. A model based on these results has been proposed to explain the mechanism balancing major histocompatibility complex (MHC) polymorphism in outbred populations.

Characterization and distribution of Mhc-DPB1 alleles in chimpanzee and rhesus macaque populations

Allelic diversity at the nonhuman primate Mhc-DPB1 locus was studied by determining exon 2 nucleotide sequences. This resulted in the detection of 17 chimpanzee (Pan troglodytes), 2 orangutan (Pongo pygmaeus) and 16 rhesus macaque (Macaca mulatta) alleles. These were compiled with primate Mhc-DPB1 nucleotide sequences that were published previously. Based upon the results, a sequence specific oligotyping method was developed allowing us to investigate the distribution of Mhc-DPB1 alleles in distinct chimpanzee and rhesus macaque colonies. Like found in humans, chimpanzee and rhesus macaque populations originating from different geographic backgrounds appear to be characterized by the presence of a few dominant Mhc-DPB1 alleles.

A procedure for detecting selection in highly variable viral genomes: evidence of positive selection in antigenic regions of capsid protein VP1 of foot-and-mouth disease virus

A new procedure is described for the detection of positive selection among sequences of viral proteins from highly variable viruses. The approach is based on the estimation of the rates of nonsynonymous to synonymous (ns/s) mutations to the overall genetic distances amongst the sequences compared. Rates of ns/s substitutions were calculated, and the individual profiles were arranged as a function of the genetic distance observed between the complete sequences. The resulting surfaces allowed identification of protein regions whose rates of ns/s substitutions were consistent with the existence of positive selection. This procedure has been applied to the study of a highly variable antigenic protein, VP1, a protein present in foot-and-mouth disease virus (FMDV). The analysis of groups of VP1 sequences corresponding to FMDV serotypes A, O and C, resulted in the identification of two regions, which contribute to an important antigenic site, where positive selection appears to operate.

Differences in the evolutionary pattern of feline panleukopenia virus and canine parvovirus

Canine parvovirus (CPV) suddenly appeared in the late 1970s after which it showed continuous antigenic changes. Virological and molecular genetic analyses mainly focused on feline panleukopenia virus (FPLV) were conducted in this study because FPLV is the suspected ancestor of CPV; the way in which FPLV evolves may help to explain the emergence of CPV. Analysis of escape mutants against FPLV-specific monoclonal antibody showed that viruses possessing CPV-like properties were not easily detected in FPLV virus stocks. Phylogenetic analysis revealed that the nonstructural protein 1 (NS1) and capsid protein 2 (VP2) genes of FPLV changed with time. A similar tendency, however, was not observed in the FPLV VP2 proteins. In contrast, the topology of the phylogenetic tree of VP2 proteins of CPV basically concurred with that of the VP2 genes. Analysis of the ratio of nonsynonymous and synonymous substitutions revealed that synonymous substitutions exceeded nonsynonymous substitutions in both the NS1 and VP2 genes of FPLV, even when the analysis focused on specific regions in the VP2 gene that are known to be located on the capsid surface. Comparison of the CPV VP2 genes revealed that nonsynonymous substitution was found to dominate over synonymous substitution in one specific region in the VP2 gene. These results suggested that FPLV has changed mainly by random genetic drift. In contrast, after the appearance of CPV, changes in the CPV VP2 gene appear to be partly selected by certain positive selection forces. CPV and FPLV are known to be closely related viruses genetically and biologically, but the evolutionary mechanisms of the two viruses appeared to be different.

Large-scale sequence comparisons reveal unusually high levels of variation in the HLA-DQB1 locus in the class II region of the human MHC

Comparison of genomic sequences flanking the HLA-DQB1 locus in the human MHC class II region reveals local sequence variation of up to 10%, which is the highest level of sequence variation found in the human genome so far. The variation is haplotype-specific and extends far beyond the transcriptional unit of the DQB1 gene, suggesting hitch-hiking along with functionally selected alleles as the most likely mechanism. All major insertions/deletions (indels) were found to be of retroviral origin and in the immediate upstream region of DQB1. Possible cis-acting effects of these indels on the transcriptional regulation of DQB1 are discussed.

New perspectives on mate choice and the MHC

A long series of studies on mice have shown that mate choice decisions can be made on the basis of individual genotype at the major histocompatibility complex (MHC), which accords well with the importance of immunocompetence in some theories of sexual selection. Recent work on other vertebrate species, including humans, indicates that MHC-based mate choice is not restricted to the genus Mus. However, its importance may vary among species as a result of differences in social and mating system structure, and perhaps genome structure. There appears to be a general preference expressed for MHC-dissimilar mates, and such MHC-disassortative mating may be involved in maintaining MHC and/or genome-wide diversity in natural populations. The strength and direction of MHC-based mating preference can vary, and may be modulated by factors such as genetic background, sex, and early life experience.

Molecular analysis and polymorphism of the DLA-DQB genes

Partial-length cDNA clones and full-length genomic clones corresponding to a complete canine DQB class II gene were isolated. Southern analyses suggested the presence of two DQB genes--one of which appeared to be a pseudogene lacking exon 2 called DQB2. The other DQB gene, called DQB1, was isolated from a genomic phage clone and contained six exons. The DQB1 clone was restriction mapped, and exon 2 was sequenced from 70 dogs. Twenty alleles were found. Most of the amino acid substitutions occurred at putative positions in the peptide binding site. Inheritance of these sequences showed Mendelian segregation with one or two alleles per dog. Cluster analysis of the nucleotide and predicted amino acid sequences subdivided the canine DQB1 alleles into four major allelic groups. The number of nonsynonymous changes was higher than the number of synonymous changes in the putative antigen recognition sites suggestive of positive selection.

New perspectives on mate choice and the MHC

A long series of studies on mice has shown that mate choice decisions can be made on the basis of individual genotype at the major histocompatibility complex (MHC), which accords well with the importance of immunocompetence in some theories of sexual selection. Recent work on other vertebrate species, including humans, indicates that MHC-based mate choice is not restricted to the genus Mus. However, its importance may vary among species as a result of differences in social and mating system structure, and perhaps genome structure. There appears to be a general preference expressed for MHC-dissimilar mates, and such MHC-disassortative mating may be involved in maintaining MHC and/or genome-wide diversity in natural populations. The strength and direction of MHC-based mating preference can vary, and may be modulated by factors such as genetic background, sex, and early life experience.

Flowering plant self-incompatibility: the molecular population genetics of Brassica S-loci

Self-incompatibility systems in different angiosperm families are reviewed, and the evidence that incompatibility has arisen several times is outlined. New data on the sequence polymorphism of self-incompatibility loci from two different angiosperm families are compared with results from other highly polymorphic loci, particularly MHC loci. We discuss what molecular genetic analyses of these sequences can tell us about the nature and maintenance of the polymorphism at self-incompatibility loci. We suggest that there is evidence for recombination at the Brassica self-incompatibility loci, so that it may be possible to discern regions that are particularly functionally important in the recognition reaction, even though the long-term maintenance of polymorphisms in these amino acid residues has caused the evolution of many other sequence differences between alleles.

The nature of polymorphism of the HLA-DRB intron sequences is lineage specific

The sequence database of HLA-DRB genes is mainly derived from mRNA analysis or has focused exclusively on the polymorphism of the 2nd exon. Little is known about the non-coding sequences of the different DRB alleles which represent about 94% of the genes. In this study we have determined the sequence of the 3' 500 bp intron 1 fragment adjacent to exon 2 in all serologically defined HLA-DRB genes and their most frequent allelic subtypes. The intron sequences turned out to be highly polymorphic. Similar to the class I introns, this variability was not characterized by random point mutations but by a highly systematic diversity reflecting the lineage-specific relationship of the HLA-DR alleles. With a few exceptions in DRBI*15, 13 and 08 as well as DRB4 and 5, the variability mirrors the serological diversity. As well as delivering insight into the genetic relationship between the different DRB alleles, these sequences will provide an extremely valuable basis for developing advanced DRB sequencing strategies for clinical purposes.

Cloning and characterization of a new swine MHC (SLA) class II DQB allele

Major histocompatibility complex (MHC) of pigs is known as swine leukocyte antigen (SLA). The cDNA encoding a new allele of SLA class II DQ beta-chain was successfully isolated from a CSK miniature pig (derived from Göttingen strain) and characterized by sequence analyses. SLA-DQB cDNA fragment encoding beta 2-domain was amplified by reverse transcriptase-polymerase chain reaction using the sequences preserved in a various vertebrates as primers. Using non-radioisotope technique with the PCR product as a probe, cDNA clone G01 was isolated from a spleen cDNA library, and nucleotide sequence of this clone was determined. This clone encompassed a whole SLA-DQ beta-chain coding region, containing a total length of 1161 nucleotides with an open reading frame (ORF) of 786 nucleotides, 5' untranslated region of 15 nucleotides, and 3' untranslated region of 360 nucleotides ending with a canonical polyadenylation signal, followed by a poly A tail. Sequence comparisons of the ORF of this clone with those of known SLA-DQB genes confirmed that this clone is a new allele (SLA-DQB*G01). Phylogenetic analysis of the nucleotide sequences of swine, human, and murine MHC class II genes indicated that SLA-DQB was more similar to HLA-DQB1 than H-2A beta. Comparison of the nucleotide and deduced amino acid sequences among SLA-DQB alleles showed that the SLA-DQ beta-chain polymorphism was found almost in beta 1-domain which contains the antigenic peptide binding sites.

Polymorphism of Mhc-DRB alleles in Cercopithecus aethiops (green monkey): generation and functionality

DRB genes have been studied for the first time in green monkeys (Cercopithecus aethiops). Eleven new DRB alleles (exon 2, exon 3) have been obtained and sequenced from cDNA. A limited number of lineages have been identified: DRB1*03 (4 alleles), DRB1*07 (3 alleles), DRB5 (1 allele), DRB*w6 (1 allele), and DRB*w7 (2 alleles). The existence of Ceae-DRB1 duplications is supported by the finding of 3 DRB1 alleles in 3 different individuals. Ceae-DRB1*0701 may be non-functional because it bears serine at position 82, which hinders molecule surface expression in mice; the allele is only found in Ceae-DRB duplicated haplotypes. Base changes in cDNA Ceae-DRB alleles are consistent with the generation of polymorphism by point mutations or short segment exchanges between alleles. The eleven green monkey DRB alleles meet the requirements for functionality as antigen-presenting molecules (perhaps, excluding DRB1*0701), since: 1) they have been isolated from cDNA and do not present deletions, insertions or stop codons: 2) structural motifs necessary for a correct folding of the molecule, for the formation of DR/DR dimers and for CD4 interactions are conserved, and 3) the number of non-synonymous substitutions is higher than the number of synonymous substitutions in the peptide binding region (PBR), while the contrary holds true for the non-PBR region.

Major histocompatibility complex variation associated with juvenile survival and parasite resistance in a large unmanaged ungulate population

Antagonistic coevolution between hosts and parasites has been proposed as a mechanism maintaining genetic diversity in both host and parasite populations. In particular, the high levels of genetic diversity widely observed at the major histocompatibility complex (MHC) of vertebrate hosts are consistent with the hypothesis of parasite-driven balancing selection acting to maintain MHC genetic diversity. To date, however, empirical evidence in support of this hypothesis, especially from natural populations, has been lacking. A large unmanaged population of Soay sheep (Ovis aries L.) is used to investigate associations between MHC variation, juvenile survival, and parasite resistance. We show in an unmanaged, nonhuman population that allelic variation within the MHC is significantly associated with differences in both juvenile survival and resistance to intestinal nematodes. Certain MHC alleles are associated with low survivorship probabilities and high levels of parasitism or vice versa. We conclude that parasites are likely to play a major role in the maintenance of MHC diversity in this population.

Hypervariability generated by natural selection in an extracellular complement-inhibiting protein of serotype M1 strains of group A Streptococcus

In many countries, M1 strains of the human pathogenic bacterium group A Streptococcus are the most common serotype recovered from patients with invasive disease episodes. Strains of this serotype express an extracellular protein that inhibits complement [streptococcal inhibitor of complement (Sic)] and is therefore believed to be a virulence factor. Comparative sequence analysis of the 915-bp sic gene in 165 M1 organisms recovered from diverse localities and infection types identified 62 alleles. Inasmuch as multilocus enzyme electrophoresis and pulsed-field gel electrophoresis previously showed that most M1 organisms represent a distinct streptococcal clone, the extent of sic gene polymorphism was unexpected. The level of polymorphism greatly exceeds that recorded for all other genes examined in serotype M1 strains. All insertions and deletions are in frame, and virtually all nucleotide substitutions alter the amino acid sequence of the Sic protein. These molecular features indicate that structural change in Sic is mediated by natural selection. Study of 70 strains recovered from two temporally distinct epidemics of streptococcal infections in the former East Germany found little sharing of Sic variants among strains recovered in the different time periods. Taken together, the data indicate that sic is a uniquely variable gene and provide insight into a potential molecular mechanism contributing to fluctuations in streptococcal disease frequency and severity.

Recent origin of HLA-DRB1 alleles and implications for human evolution

The HLA class I and class II loci are the most highly polymorphic coding regions in the human genome. Based on the similarity of the coding sequences of alleles between species, it has been claimed that the HLA polymorphism is ancient and predates the separation of human (Homo) and chimpanzee (Pan), 4-7.4 Myr ago. Analysis of intron sequences, however, provides support for a more recent origin and for rapid generation of alleles at the HLA class II DRB1 locus. The human DRB1 alleles can be divided into groups (allelic lineages); most of these lineages have diverged from each other before the separation of Homo and Pan. Alleles within such a lineage, however, appear to be, on average, 250,000 years old, implying that the vast majority (greater than 90%) of the more than 135 contemporary human DRB1 alleles have been generated after the separation of Homo and Pan. The coalescence time of alleles within allelic lineages indicates that the effective population size (Ne) for early hominids (over the last 1 Myr) was approximately 10(4) individuals, similar to estimates based on other nuclear loci and mitochondrial DNA. With a single exception, the genetic mechanisms (gene conversion and point mutation) that have diversified the exon-2 sequences do not appear to extend into the adjacent intron sequences. The part of exon 2 encoding the beta-sheet evolves in concert with the surrounding introns, while the alpha-helix appears to have been subjected to gene conversion-like events, suggesting that such exchange events are highly localised and occur over extremely short sequence tracts.

Major histocompatibility complex class II alleles in an Uygur population in the Silk Route of Northwest China

HLA class II (DRB1, DQA1, DQB1 and DPB1) genotyping was performed in 57 unrelated Uygur individuals inhabiting the northwestern China area by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Among 98 DRB1 alleles tested, 23 alleles were detected, and DRB1*0701 (16.7%) and DRB1*0301 (14.0%) were the most and the second most common alleles, respectively. In 8 DQA1 alleles detected, DQA1*05 (26.3%), DQA1*03 (21.9%) and DQA1*0201 (21.1%) were very frequent. Of 21 DQB1 alleles tested, 13 were observed. Among them, DQB1*02 was highly predominant with the gene frequency of 32.5%. Of 46 DPB1 alleles tested, 15 were detected, among which DPB1*0401 (31.6%) was the most frequent. Two haplotypes predominate clearly; DRB1*0701-DQA1*0201-DQB1*02 (15.5%) and DRB1*0301-DQA1*05-DQB1*02 (12.6%). The dendrogram constructed by the neighbour-joining (NJ) method based on the allele frequencies of the DRB1, DQA1, DQB1 and DPB1 genes of 13 representative populations all over the world suggested that Uygur belonged to the Asian group and lay at the closest genetic distance to a Kazak population inhabiting the same area.