Like Wings of a Bird: Functional Divergence and Complementarity between HLA-A and HLA-B Molecules

Human leukocyte antigen (HLA) genes are among the most polymorphic of our genome, as a likely consequence of balancing selection related to their central role in adaptive immunity. HLA-A and HLA-B genes were recently suggested to evolve through a model of joint divergent asymmetric selection conferring all human populations, including those with severe loss of diversity, an equivalent immune potential. However, the mechanisms by which these two genes might undergo joint evolution while displaying very distinct allelic profiles in populations are still unknown. To address this issue, we carried out extensive data analyses (among which factorial correspondence analysis and linear modeling) on 2,909 common and rare HLA-A, HLA-B, and HLA-C alleles and 200,000 simulated pathogenic peptides by taking into account sequence variation, predicted peptide-binding affinity and HLA allele frequencies in 123 populations worldwide. Our results show that HLA-A and HLA-B (but not HLA-C) molecules maintain considerable functional divergence in almost all populations, which likely plays an instrumental role in their immune defense. We also provide robust evidence of functional complementarity between HLA-A and HLA-B molecules, which display asymmetric relationships in terms of amino acid diversity at both inter- and intraprotein levels and in terms of promiscuous or fastidious peptide-binding specificities. Like two wings of a flying bird, the functional complementarity of HLA-A and HLA-B is a perfect example, in our genome, of duplicated genes sharing their capacity of assuming common vital functions while being submitted to complex and sometimes distinct environmental pressures.     

Common chimpanzees have greater diversity than humans at two of the three highly polymorphic MHC class I genes

 MHC class I polymorphism improves the defense of vertebrate species against viruses and other intracellular pathogens. To see how polymorphism at the same class I genes can evolve in different species we compared the MHC-A, MHC-B, and MHC-C loci of common chimpanzees and humans. Diversity in 23 Patr-A, 32 Patr-B, and 18 Patr-C alleles obtained from study of 48 chimpanzees was compared to diversity in 66 HLA-A, 149 HLA-B, and 41 HLA-C alleles obtained from a study of over 1 million humans. At each locus, alleles group hierarchically into families and then lineages. No alleles or families are shared by the two species, commonality being seen only at the lineage level. The overall nucleotide sequence diversity of MHC class I is estimated to be greater for modern chimpanzees than humans. Considering the numbers of lineages, families, and alleles, Patr-B and Patr-C have greater diversity than the HLA-B and HLA-C, respectively. In contrast, Patr-A has less polymorphism than HLA-A, due to the absence of A2 lineage alleles. The results are consistent with ancestral humans having passed through a narrower population bottleneck than chimpanzees, and with pathogen-mediated selection having favored either preservation of A2 lineage alleles on the human line and/or their extinction on the chimpanzee line. Received: 8 December 1999 / Accepted: 30 December 1999     

Evidence for an HLA-C-like locus in the orangutan Pongo pygmaeus

 HLA-B and C are related class I genes which are believed to have arisen by duplication of a common ancestor. Previous study showed the presence of orthologues for both HLA-B and C in African apes but only for HLA-B in Asian apes. These observations suggested that the primate C locus evolved subsequent to the divergence of the Pongidae and Hominidae. From an analysis of orangutan Tengku two HLA-C-like alleles (Popy C*0101 and Popy C*0201) were defined as well as three HLA-B-like (Popy-B) alleles. By contrast, no Popy-C alleles were obtained from orangutan Hati, although three Popy-B alleles were defined. Thus an HLA-C-like locus exists in the orangutan (as well as a duplicated B locus), implying that the primate C locus evolved prior to the divergence of the Pongidae and Hominidae and is at least 12–13 million years old. Uncertain is whether all orangutan MHC haplotypes contain a C locus, as the failure to find C alleles in some individuals could be due to a mispairing of HLA-C-specific primers with certain Popy-C alleles. These results raise the possibilities that other primate species have a C locus and that the regulation of natural killer cells by C allotypes evolved earlier in primate evolution than has been thought. Received: 18 January 1999 / Revised: 23 March 1999     

Patterns of reticulate evolution for the classical class I and II HLA loci

 Some alleles of the major histocompatibility complex (MHC) genes have a reticulate pattern of evolution, probably resulting from the exchange of segments by gene conversion or recombination. Here we compare the extent and patterns of reticulate evolution among the classical class I and class II loci of the human MHC using the recently developed compatibility and partition matrix methods. A complex pattern is revealed with substantial differences among loci in the extent and pattern of reticulation. Extremely high levels of reticulation are observed at HLA-B and HLA-DPB1, high levels at HLA-A and HLA-DRB1, moderate levels at HLA-C and HLA-DQB1, and low levels at HLA-DQA1. The reticulate events are concentrated in the exons encoding the highly variable, peptide-binding domains, suggesting that the sequence combinations produced by these events are maintained by natural selection. Received: 3 December 1997 / Revised: 30 March 1998     

Molecular mapping of a new public HLA class I epitope shared by all HLA-B and HLA-C antigens and defined by a monoclonal antibody

It has previously been shown that a mouse monoclonal antibody, designated 4E, reacts with an epitope common to all HLA-B and -C antigens and those of the HLA-Aw19 cross-reactive group, namely, HLA-A29,-A30, -A31, -A32, -Aw33, and -Aw74. In order to pinpoint the amino acid residues which comprise the public specificity recognized by 4E, an HLA-A29 cDNA clone was isolated and its predicted amino acid sequence compared with those of other clonedHLA class I genes. The isolated HLA-A29 cDNA corresponded to the rarer of the twoA29 variant alleles,A29.1. Two amino acid residues of HLA-A29.1, gln-144 and arg-151, were found in all 24HLA-B andHLA-C alleles examined but were present in only one of 15HLA-A alleles for which sequence data are available. Importantly, this exceptional allele wasHLA-A32, another member of the HLA-Aw19 cross-reactive group. Gln-144 and arg-151 should be capable of jointly contributing to the binding site for 4E, as they are situated in successive alpha-helical subregions and are predicted to be juxtaposed in the three-dimensional HLA molecule. Four other residues in the first or second external domains of HLA-A29.1 (thr-9, leu-62, gln-63, and his-102) were unique among theHLA-A alleles, but none of these was found in corresponding positions ofHLA-B or-C alleles and thus failed to correlate with presence or absence of the 4E determinant. These observations are consistent with the notion that gln-144 and arg-151 define a determinant common to HLA-B, HLA-C, and the HLA-Awl9 cross-reactive group and the binding site of the monoclonal antibody 4E.     

The molecular origin and consequences of escape from miRNA regulation by HLA-C alleles

Differential expression of human leukocyte antigen C (HLA-C) allotypes is mediated by the binding of a microRNA, miR-148a, to the 3′ untranslated region of some, but not all, HLA-C alleles. The binding results in lower levels of HLA-C expression, which is associated with higher levels of HIV-1 viral load among infected individuals. The alternative set of HLA-C alleles has several substitutions in the miR-148a binding site that prevent binding and HLA-C downregulation; these high-expression alleles associate with control of HIV-1 viral load. We show that the common ancestor of all extant HLA-C alleles was suppressed by miR-148a. Substitutions that prevent miR-148a binding arose by a sequence exchange event between an HLA-C allele and an HLA-B (MIM 142830) allele of a B^∗07-like lineage. The event occurred 3–5 million years ago, resulting in an HLA-C variant that escape from miR-148a downregulation. We present evidence suggesting that selection played a role in the successful spread of the HLA-C escape alleles, giving rise to 7 of the 14 extant HLA-C lineages. Notably, critical peptide and KIR binding residues of the escape variants have selectively converged to resemble the sequence of their inhibited counterparts, such that the inhibited and escape groupings differ primarily by their levels of expression.     

Restimulation in secondary MLC by a non-D-locus determinant within the MHC

By testing a family in which one offspring had inherited a chromosome including a recombination between theHLA-B andD loci, we sought to obtain some insight into whether determinants other than those atHLA-D are capable of restimulating in a secondary MLC. Results obtained in the primary MLC indicated that the recombinant child did not express products of theHLA-D region normally associated with this haplotype. When sensitization in a primary MLC was to the entire haplotype, everyone who had the sensitizing haplotype restimulated strongly and specifically. In addition, however, weak to moderate stimulation was obtained when the cells of the recombinant child were used to restimulate. Presumably these cells possessed determinants coded for within theHLA-A toB chromosomal segments of the sensitizing haplotype, but not those coded for by theHLA-D locus. Our results indicate that a simple structure atHLA-D is not the sole factor in the secondary MLC. Either the products of loci outside theHLA-D locus can also restimulate or the recombination occurred within theD locus, suggesting that theD region contains more than one locus coding for restimulating determinants.     

Parakeratosis Mibelli and skin carcinoma

Summary Bf allele frequencies in a material of 172 unrelated Norwegians are given. Bf/HLA linkage relations in 49 informative matings with 178 children, and Bf/HLA association data of a material of 212 Bf-HLA haplotypes are presented. Of 171 informative meioses, there were no Bf-HLA-B recombinations, while 3 out of 158 Bf-HLA-A informative meioses showed recombination. There is significant association between the Bf ^F and the HLA-BW 35 allele.It is concluded that the Bf locus is situated on the HLA-B side of HLA-A within the HLA region, in very close proximity to HLA-B.     

Frequent segmental sequence exchanges and rapid gene duplication characterize the<Emphasis Type="Italic"> MHC</Emphasis> class I genes in lemurs

Major histocompatibility complex (MHC) class I genes have complicated and profound evolutionary histories. To reconstruct and better understand their histories, partial class I genes (exon 2–intron 2–exon 3) were sequenced in a sampling of prosimians (Strepsirhini, Primates). In total, we detected 117 different sequences from 36 Malagasy prosimians (lemurs) and 1 non-Malagasy prosimian (galago) representing 4 families, 7 genera, and 13 species. Unlike the MHC class II genes (MHC-DRB), MHC class I genes show a generally genus-specific mode of evolution in lemurs. Additionally, no prosimian class I loci were found to be orthologous to HLA genes, even at highly conserved loci (such as HLA-E, HLA-F). Phylogenetic analysis indicates that nucleotide diversity among loci was very small and the persistence time of the polymorphisms was short, suggesting that the origin of the lemur MHC class I genes detected in this study was relatively recent. The evolutionary mode of these genes is similar to that of classic HLA genes, HLA-A, HLA-B, and HLA-C, in terms of their recent origin and rarity of pseudogenes, and differs from them with respect to the degree of gene duplications. From the viewpoint of MHC genes evolution, some interlocus sequence exchanges were apparently observed in the lemur lineage upon phylogenetic and amino acid motif analyses. This is also in contrast to the evolutionary mode of HLA genes, where intralocus exchanges have certainly occurred but few interlocus exchanges have taken place. Consequently, the gene conversion model for explaining the generation of the MHC diversity among different loci can be thought to play more important roles in the evolution of lemur MHC class I genes than in that of HLA genes.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Genetic Epidemiology of Glioblastoma Multiforme: Confirmatory and New Findings from Analyses of Human Leukocyte Antigen Alleles and Motifs

Background

Human leukocyte antigen (HLA) class I genes mediate cytotoxic T-lymphocyte responses and natural killer cell function. In a previous study, several HLA-B and HLA-C alleles and haplotypes were positively or negatively associated with the occurrence and prognosis of glioblastoma multiforme (GBM).

Methodology/Principal Findings

As an extension of the Upper Midwest Health Study, we have performed HLA genotyping for 149 GBM patients and 149 healthy control subjects from a non-metropolitan population consisting almost exclusively of European Americans. Conditional logistic regression models did not reproduce the association of HLA-B*07 or the B*07-Cw*07 haplotype with GBM. Nonetheless, HLA-A*32, which has previously been shown to predispose GBM patients to a favorable prognosis, was negatively associated with occurrence of GBM (odds ratio = 0.41, p = 0.04 by univariate analysis). Other alleles (A*29, A*30, A*31 and A*33) within the A19 serology group to which A*32 belongs showed inconsistent trends. Sequencing-based HLA-A genotyping established that A*3201 was the single A*32 allele underlying the observed association. Additional evaluation of HLA-A promoter and exon 1 sequences did not detect any unexpected single nucleotide polymorphisms that could suggest differential allelic expression. Further analyses restricted to female GBM cases and controls revealed a second association with a specific HLA-B sequence motif corresponding to Bw4-80Ile (odds ratio = 2.71, p = 0.02).

Conclusions/Significance

HLA-A allelic product encoded by A*3201 is likely to be functionally important to GBM. The novel, sex-specific association will require further confirmation in other representative study populations.     

Quantitative estimation of HLA-A and HLA-B antigens carrying the Bw4 supertypic specificity in human peripheral blood lymphocytes

A monoclonal antibody specific for HLA-Bw4 was employed for the quantitative estimation of class I antigens on human peripheral blood lymphocytes and on the human macrophage cell line U 937. The epitopes reactive with the HLA-Bw4-specific antibody, which are present on different antigens coded by the HLA-A or the HLA-B locus, were characterized in terms of equilibrium and kinetic binding parameters. The level of expression of class I antigens on human lymphocytes was found to be in direct proportion to the gene dose. Variations between donors of the same phenotype were of minor importance. Estimation of the association constant, association rate constant, and half-life of dissociation for the interaction of the antibody with lymphocytes heterozygous for either HLA-A24, HLA-A32, HLA-A9, or HLA-Bw4 strongly suggested that the public antigenic determinant shared by these antigens is identical.     

Revision of the bulbous Iridaceae of North America

The native bulbous Iridaceae of North America excluding Mexico comprise seven species in four genera. The type species ofAlophia Herb. (1836) is shown to belong toEustylis Engelm. & Gray (1836), hence the latter genus is reduced to synonymy. The well known genusAlophia (sensu auct. non Herb.) (Herbertia Sweet) now takes the nameTrifurcia. As a result of the new interpretation ofAlophia, new combinations are made for the North American and Latin American species ofTrifurcia andAlophia. New chromosome counts are given forNemastylis floridana (2n = 56) andTrifurcia lahue (2n = 56).     

Mapping SB in relation to HLA and GL01 using cells from first-cousin marriage offspring

A newly defined leukocyte antigen system, the secondary B-cell (SB) system, was shown to be linked to HLA. The SB marker has been investigated in lymphocyte donors presumed to be genetically homozygous for HLA-A through HLA-D/DR by virtue of descent from a first-cousin marriage and of phenotypic homozygosity for these HLA markers. Of 19 donors, 3 were found to be heterozygous for SB. Studies of the families of these three donors could not distinguish with certainty whether the heterozygosity resulted from SB/DR recombination or from “pseudohomozygosity” for HLA-A through -D/DR by inheritance of two genetically unrelated but similar haplotypes. However, our data favored the occurrence of SB/DR recombination with a meiotic distance perhaps as large as 3.3 cMorgan. Recombinations were identified which mapped SB between HLA-B and GL01. These studies demonstrate the usefulness of cells from first-cousin marriage offspring in mapping a polymorphic genetic system.     

FISH and GISH analysis of the genomic relationships amongPanax species

Ginseng is a well-known medicinal plant that has been used as an anti-aging agent for many years in East Asia. In the genusPanax, only three species,P. ginseng (Oriental ginseng),P. quinquefolius (American ginseng) andP. notoginseng (Chinese ginseng), are currently considered to be important medicinal herbs. Despite the increase in their breeding value, molecular cytogenetic information on the species is not available. To analyze the genomic relationships among thePanax species, FISH (fluorescencein situ hybridization) and GISH (genomicin situ hybridization) techniques were applied. FISH analysis revealed that the 45S and 5S rRNA genes ofP. notoginseng (2n=2x=24) andP. ginseng (2n=4x=48) cluster on a single locus on different chromosomes, whileP. quinquefolius (2n=4x=48),P. japonicus (2n=4x=48), and Korean wild ginseng (2n =4x= 48) had one locus of the 45S rRNA gene and two loci of the 5S rRNA gene, respectively. GISH analysis using genomic DNA as a probe detected strong cross-hybridization of genomes betweenP. ginseng andP. quinquefolius. GISH analysis of other species showed weak or no distinct signals on the chromosomes. Our data revealed thatP. ginseng andP. quinquefolius showed the highest degree of homology, indicating that these species diverged in most recent years.     

Chromosome numbers of Orthocarpus and related monotypic genera (Scrophulariaceae: Subtribe castillejinae)

A chromosome number ofn=12 is reported for the three monotypic genera of subtribe Castillejinae:Clevelandia beldingii, Gentrya racemosa, andOpicopephalus angustifolius. Chromosome numbers ofOrthocarpus correspond mostly with current infrageneric classification. SubgenusTriphysaria hasn=11.Orthocarpus sectionsCastillejoides andCordylanthoides, which are closely related toCastilleja (x=12) and the three monotypic genera above, haven=12 with aneuploid reductions ton=10 inO. linearilobus andn=11 inO. lacerus (a species also withn=12). Tetraploids are found in two species.O. brevistylus (n=24) andO. hispidus (n=12, 24). The polyploid.O. laciniatus (n=36, 48) of Peru is postulated to be of hybrid origin between a species ofCastilleja andOrthocarpus attenuatus. SubgenusOrthocarpus sectionOrthocarpus, which hasn=14 in all species except.O. bracteosus (n=15), stands apart both morphologically and in chromosome number from the remainder of the genus.     

Fine mapping of the psoriasis susceptibility locus PSORS1 supports HLA-C as the susceptibility gene in the Han Chinese population

PSORS1 (psoriasis susceptibility gene 1) is a major susceptibility locus for psoriasis. Several fine-mapping studies have highlighted a 300-kb candidate region of PSORS1 where multiple biologically plausible candidate genes were suggested. The most recent study has indicated HLA-Cw6 as the primary PSORS1 risk allele within the candidate region in a Caucasian population. In this study, a family-based association analysis of the PSORS1 locus was performed by analyzing 10 polymorphic microsatellite markers from the PSORS1 region as well as HLA-B, HLA-C and CDSN loci in 163 Chinese families of psoriasis. Five marker loci show strong evidence (P<10⁻³), and one marker locus shows weak evidence (P=0.04) for association. The haplotype cluster analysis showed that all the risk haplotypes are Cw6 positive and share a 369-kb region of homologous marker alleles which carries all the risk alleles, including HLA-Cw6 and CDSN*TTC, identified in this study. The recombinant haplotype analysis of the HLA-Cw6 and CDSN*TTC alleles in 228 Chinese families showed that the HLA-Cw6⁻/CDSN*TTC⁺ recombinant haplotype is clearly not associated with risk for psoriasis (TNT=29:57, p=0.0025) in a Chinese population, suggesting that the CDSN*TTC allele itself does not confer risk without the presence of the HLA-Cw6 allele. The further exclusion analysis of the non-risk HLA-Cw6⁻/CDSN*TTC⁺ recombinant haplotypes with common recombination breakpoints has allowed us to refine the location of PSORS1 to a small candidate region. Finally, we performed a conditional linkage analysis and showed that the HLA-Cw6 is a major risk allele but does not explain the full linkage evidence of the PSORS1 locus in a Chinese population. By performing a series of family-based association analyses of haplotypes as well as an exclusion analysis of recombinant haplotypes, we were able to refine the PSORS1 gene to a small critical region where HLA-C is a strong candidate to be the PSORS1 susceptibility gene.     

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.     

Cloning and analysis of HLA class I cDNA encoding a new HLA-C specificity Cx52

HLA-C loci frequently have an unclassifiable blank (CwBL) specificity. It is unclear whether HLA-C specificities associated with the haplotypes of A24 Bw52 CwBL DR2 DQw1 and Aw33 B44 CwBL DRw13 DQw1 in Japanese (tentatively named Cx52 and Cx44, respectively) really exist. Southern hybridization experiments revealed that restriction enzyme-cleaved genomic DNA from AKIBA, consanguineous HLA homozygote, two other homozygotes with the former haplotype, and three homozygous cells with the latter haplotype hybridized strongly with an HLA-C-specific probe. We have screened the cDNA library constructed from AKIBA to isolate cDNA clones encoding the putative Cx52 antigen, and picked up 103 cDNA clones with HLA-class I DNA probes as possible candidates. By restriction enzyme mapping and Southern hybridization of selected clones, we identified three isotypes of cDNA clones, pA01, pB55, and pC68, which appeared to encode A24, Bw52, and Cx52, respectively. The nucleotide sequence of pC68 showed higher homology with exons of the HLA-C gene than with those of the HLA-A and HLA-B genes, especially in exons 6–8 which include the HLA-C-specific region. Comparison of amino acid sequences showed more than 86% homology among Cw1, Cw2, Cw3, and new pC68-encoded Cx52 proteins. These results support the notion that the inability to define C antigens serologically in this Cx52 haplotype is not due to a HLA-C gene deletion or mutation, but to the absence of typing sera.     

HLA-D typing in 111 multiple sclerosis patients: Distribution of four HLA-D alleles

Multiple sclerosis patients-111 of them-were typed for theHLA-D allelesw1, w2, andw3 and a new determinant, EI. Statistically significant increase was obtained for thew2 andw3 frequencies as compared to healthy controls. The distribution of HLA-D phenotypes among MS patients revealed a good fit according to the Hardy-Weinberg law. By calculating linkage disequilibrium parameters, theHLA-B7,Dw2 allele combination was found to be more closely associated than in normal controls, whereas forHLA-Bw35, Dw1, no linkage disequilibrium could be detected in the patients' group. From these data we conclude that in multiple sclerosis, the disturbance affects the frequency ofDw3 and the gametic association between alleles of theHLA-B andD loci, as well as the already known increase ofHLA-B7 andDw2.