Genomic Organization Around the Centromeric End of the HLA Class I Region: Large-Scale Sequence Analysis

We previously sequenced two regions around the centromeric end of HLA class I and the boundary between class I and class III. In this paper we analyze the two regions of about 385 kb and confirm, giving a new line of evidence, that the following two pairs of the genomic segments were duplicated in evolution: (i) a 43-kb genomic segment including the HLA-B gene showing the highest polymorphism among the classical HLA class I loci (class Ia) and a 40-kb segment including the HLA-C locus showing the lowest polymorphism and (ii) a 52-kb segment including the MIC (MHC class I chain related gene) B and a 35-kb segment including MICA. We also found that repetitive elements such as SINEs, LINEs, and LTRs occupy as much as 47% of nucleotides in this 385-kb region. This unusually high content of repetitive elements indicates that repeat-mediated rearrangements have frequently occurred in the evolutionary history of the HLA class Ia region. Analysis of LINE compositions within the two pairs of duplicated segments revealed that (i) LINEs in these regions had been dispersed prior to both the duplication of the HLA-B and -C loci and the duplication of the MICB and MICA loci, and (ii) the divergence of the HLA-B and -C loci occurred prior to the duplication of the MICA and MICB loci. To find novel genes responsible for HLA class I-associated or other diseases, we performed computer analysis applying GenScan and GRAIL to GenBank's dbEST. As a result, at least five as yet uncharacterized genes were newly mapped on the HLA class I centromeric region studied. These novel genes should be analyzed further to determine their relationships to diseases associated with this region. Received: 16 June 1998 / Accepted: 18 August 1998     

Refined Linkage Map of Chromosome 5 in the Region of the Spinal Muscular Atrophy Gene

The genetic map in the region of human chromosome 5 that harbors the gene for autosomal recessive forms of spinal muscular atrophy (SMA) has been refined by a multilocus linkage study in 50 SMA-segregating families. Among six markers spanning 8 cM for combined sexes, four were shown to be tightly linked to the SMA locus. Multipoint linkage analysis was used to establish the best estimate of the SMA gene location. Our data suggest that the most likely location for the SMA locus is between blocks AFM114ye7 (D5S465)/EF5.15 (D5S125) and MAP-1B/JK53 (D5S112) at a sex-combined genetic distance of 2.4 and 1.7 cM, respectively. Thus the SMA gene lies in the 4-cM region between these two blocks. This information is of primary importance for designing strategies for isolating the SMA gene.     

Coevolution of PERB11 (MIC) and HLA Class I Genes with HERV-16 and Retroelements by Extended Genomic Duplication

The recent availability of genomic sequence information for the class I region of the MHC has provided an opportunity to examine the genomic organization of HLA class I (HLAcI) and PERB11/MIC genes with a view to explaining their evolution from the perspective of extended genomic duplications rather than by simple gene duplications and/or gene conversion events. Analysis of genomic sequence from two regions of the MHC (the alpha- and beta-blocks) revealed that at least 6 PERB11 and 14 HLAcI genes, pseudogenes, and gene fragments are contained within extended duplicated segments. Each segment was searched for the presence of shared (paralogous) retroelements by RepeatMasker in order to use them as markers of evolution, genetic rearrangements, and evidence of segmental duplications. Shared Alu elements and other retroelements allowed the duplicated segments to be classified into five distinct groups (A to E) that could be further distilled down to an ancient preduplication segment containing a HLA and PERB11 gene, an endogenous retrovirus (HERV-16), and distinctive retroelements. The breakpoints within and between the different HLAcI segments were found mainly within the PERB11 and HLA genes, HERV-16, and other retroelements, suggesting that the latter have played a major role in duplication and indel events leading to the present organization of PERB11 and HLAcI genes. On the basis of the features contained within the segments, a coevolutionary model premised on tandem duplication of single and multipartite genomic segments is proposed. The model is used to explain the origins and genomic organization of retroelements, HERV-16, DNA transposons, PERB11, and HLAcI genes as distinct segmental combinations within the alpha- and beta-blocks of the human MHC. Received: 5 December 1998 / Accepted: 27 January 1999     

Human trophoblast and the choriocarcinoma cell line BeWo express a truncated HLA Class I molecule

Extravillous trophoblast from normal human placenta has been shown previously to express an unusual form of HLA class I molecule, as does a choriocarcinoma cell line, BeWo. This molecule has a H chain of approximately 40 kDa and appears to be nonpolymorphic. We have isolated and sequenced a HLA class I cDNA clone, which probably corresponds to this molecule, from a library derived from BeWo. The nucleotide sequence shows a high degree of homology with the published sequence of a genomic clone, HLA 6.0, which is the product of a class I locus other than A, B, or C, (provisionally designated "HLA G"). The expressed product of this locus has not previously been localized. We have used the polymerase chain reaction to demonstrate the presence of similar HLA class I sequences in cDNA from normal extravillous trophoblast. Although there is some nucleotide sequence polymorphism the amino acid sequence of this molecule is conserved. It is therefore unlikely to provoke immune responses even though it is found at the fetal-maternal interface.     

Class I HLA antigens in two long-separated populations: Melanesians and South Amerinds

Class I HLA antigens have been compared in 5,835 Melanesians of Papua New Guinea and 2,028 Amerindians of South America. The sample includes 50 PNGMel ethnolinguistic groups and 22 SAmlnd groups. Both carry 15 serologically defined antigens and an undefined C allele. Except for A2 in Papua New Guinea and Cwl in South America, these antigens are widely distributed in their respective populations. Nine (A2 and A24, B39, B60 and B62, and Cwl, Cw3, Cw4, and Cw7) are common to both. This commonality suggests that these two populations derive from an ancestral population with less polymorphism than modern East Asians. In both populations several theoretically possible haplotypes were absent, and other haplo-types were in positive disequilibrium in both. The parallels in disequilibria suggest that haplotypes are subject to selective forces acting on the level of allelic interaction. Based on three locus haplotype frequencies, the PNGMel groups form five clusters with internally typical linguistic and geographic characteristics and a miscellaneous category, but Samlnd groups show no cluster. © 1995 Wiley-Liss, Inc.     

HLA Class I and Class II Alleles and Haplotypes Confirm the Berber Origin of the Present Day Tunisian Population

In view of its distinct geographical location and relatively small area, Tunisia witnessed the presence of many civilizations and ethnic groups throughout history, thereby questioning the origin of present-day Tunisian population. We investigated HLA class I and class II gene profiles in Tunisians, and compared this profile with those of Mediterranean and Sub-Sahara African populations. A total of 376 unrelated Tunisian individuals of both genders were genotyped for HLA class I (A, B) and class II (DRB1, DQB1), using reverse dot-blot hybridization (PCR-SSO) method. Statistical analysis was performed using Arlequin software. Phylogenetic trees were constructed by DISPAN software, and correspondence analysis was carried out by VISTA software. One hundred fifty-three HLA alleles were identified in the studied sample, which comprised 41, 50, 40 and 22 alleles at HLA-A,-B,-DRB1 and -DQB1 loci, respectively. The most frequent alleles were HLA-A*02:01 (16.76%), HLA-B*44:02/03 (17.82%), HLA-DRB1*07:01 (19.02%), and HLA-DQB1*03:01 (17.95%). Four-locus haplotype analysis identified HLA-A*02:01-B*50:01-DRB1*07:01-DQB1*02:02 (2.2%) as the common haplotype in Tunisians. Compared to other nearby populations, Tunisians appear to be genetically related to Western Mediterranean population, in particular North Africans and Berbers. In conclusion, HLA genotype results indicate that Tunisians are related to present-day North Africans, Berbers and to Iberians, but not to Eastern Arabs (Palestinians, Jordanians and Lebanese). This suggests that the genetic contribution of Arab invasion of 7^th-11^th century A.D. had little impact of the North African gene pool.     

Immunogenicity of HLA Class I and II Double Restricted Influenza A-Derived Peptides

The aim of the present study was to identify influenza A-derived peptides which bind to both HLA class I and -II molecules and by immunization lead to both HLA class I and class II restricted immune responses. Eight influenza A-derived 9-11mer peptides with simultaneous binding to both HLA-A*02:01 and HLA-DRB1*01:01 molecules were identified by bioinformatics and biochemical technology. Immunization of transgenic HLA-A*02:01/HLA-DRB1*01:01 mice with four of these double binding peptides gave rise to both HLA class I and class II restricted responses by CD8 and CD4 T cells, respectively, whereas four of the double binding peptides did result in HLA-A*02:01 restricted responses only. According to their cytokine profile, the CD4 T cell responses were of the Th2 type. In influenza infected mice, we were unable to detect natural processing in vivo of the double restricted peptides and in line with this, peptide vaccination did not decrease virus titres in the lungs of intranasally influenza challenged mice. Our data show that HLA class I and class II double binding peptides can be identified by bioinformatics and biochemical technology. By immunization, double binding peptides can give rise to both HLA class I and class I restricted responses, a quality which might be of potential interest for peptide-based vaccine development.     

The Transcription Map of the 13q14 Region Frequently Deleted in B-cell Chronic Lymphocytic Leukemia

Deletions in the region located between the STS markers D13S1168 and D13S25 on chromosome 13 are the most frequent genomic changes in patients with B-cell chronic lymphocytic leukemia (B-CLL). After sequencing of this region, two novel candidate genes were identified: C13orf1(chromosome 13 open reading frame 1) and PLCC (putative large CLL candidate). Analysis of the repeat distribution revealed two subregions differing in composition of repetitious DNA and gene organization. The interval D13S1168–D13S319 contains 131 Alu repeats accounting for 24.8% of its length, whereas the interval GCT16C05–D13S25, which is no more than 180 kb away from the former one is extremely poor in Alu repeats (4.1% of the total length). Both intervals contain almost the same amount of the LINE-type repeats L1 and L2 (20.3 and 21.24%, respectively). In the chromosomal region studied, 29 Alu repeats were found to belong to the evolutionary young subfamily Y, which is still capable of amplifying. A considerable proportion of repeats of this type with similar nucleotide sequences may contribute to the recombinational activity of the chromosomal region 13q14.3, which is responsible for its rearrangements in some tumors in humans.     

HLA-DPB1 and HLA Class I Confer Risk of and Protection from Narcolepsy

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1^∗01:02-DQB1^∗06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1^∗01:03-DPB1^∗04:02 (DP0402; odds ratio [OR] = 0.51 [0.38–0.67], p = 1.01 × 10⁻⁶) and HLA-DPA1^∗01:03-DPB1^∗04:01 (DP0401; OR = 0.61 [0.47–0.80], p = 2.07 × 10⁻⁴) and predisposing effects of HLA-DPB1^∗05:01 in Asians (OR = 1.76 [1.34–2.31], p = 4.71 × 10⁻⁰⁵). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38–0.55] p = 8.99 × 10⁻¹⁷) and DP0501 (OR = 1.38 [1.18–1.61], p = 7.11 × 10⁻⁵). HLA-class-II-independent associations with HLA-A^∗11:01 (OR = 1.32 [1.13–1.54], p = 4.92 × 10⁻⁴), HLA-B^∗35:03 (OR = 1.96 [1.41–2.70], p = 5.14 × 10⁻⁵), and HLA-B^∗51:01 (OR = 1.49 [1.25–1.78], p = 1.09 × 10⁻⁵) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.     

Genomic and Phylogenetic Analysis of the Human CD1 and HLA Class I Multicopy Genes

The human CD1 proteins belong to a lipid-glycolipid antigen-presenting gene family and are related in structure and function to the MHC class I molecules. Previous mapping and DNA hybridization studies have shown that five linked genes located within a cluster on human chromosome 1q22-23 encode the CD1 protein family. We have analyzed the complete genomic sequence of the human CD1 gene cluster and found that the five active genes are distributed over 175,600 nucleotides and separated by four expanded intervening genomic regions (IGRs) ranging in length between 20 and 68 kb. The IGRs are composed mostly of retroelements including five full-length L1 PA sequences and various pseudogenes. Some L1 sequences have acted as receptors for other subtypes or families of retroelements. Alu molecular clocks that have evolved during primate history are found distributed within the HLA class I duplicated segments (duplicons) but not within the duplicons of CD1. Phylogeny of the alpha3 domain of the class I-like superfamily of proteins shows that the CD1 cluster is well separated from HLA class I by a number of superfamily members including MIC (PERB11), HFE, Zn-alpha2-GP, FcRn, and MR1. Phylogenetically, the human CD1 sequences are interspersed by CD1 sequences from other mammalian species, whereas the human HLA class I sequences cluster together and are separated from the other mammalian sequences. Genomic and phylogenetic analyses support the view that the human CD1 gene copies were duplicated prior to the evolution of primates and the bulk of the HLA class I genes found in humans. In contrast to the HLA class I genomic structure, the human CD1 duplicons are smaller in size, they lack Alu clocks, and they are interrupted by IGRs at least 4 to 14 times longer than the CD1 genes themselves. The IGRs seem to have been created as "buffer zones" to protect the CD1 genes from disruption by transposable elements.     

Nucleotide Sequencing Analysis of the 146-Kilobase Segment around theIkBLandMICAGenes at the Centromeric End of the HLA Class I Region

To elucidate the complete gene structure and to identify new genes involved in the development ofHLAclass I antigen-associated diseases in the class I region of the human major histocompatibility complex on chromosome 6, a YAC clone (745D12) covering the 146-kb segment around theIkBLandMICAloci was isolated from a YAC library constructed from the B-cell line, BOLETH. A physical map of this region was constructed by isolation of overlapping cosmid clones derived from 745D12. Of these, five contiguous cosmids were chosen for DNA sequencing by the shotgun strategy to give a single contig of 146,601 bp from 2.8 kb telomeric of theIkBLgene to exon 6 ofMICA.This region was confirmed to contain five known genes,IkBL, BAT1, MICB, P5-1,andHLA-X(class I fragment), from centromere to telomere, and their exon–intron organizations were determined. The3.8-1homologue gene (3.8-1-hom) showing 99.7% identity with the3.8-1cDNA clone, which was originally isolated using the 3.8-kbEcoRI fragment between theHLA-54/Hand theHLA-Ggenes, was detected betweenMICAandMICBand was suggested to represent the cognate3.8-1genomic sequence from which the cDNA clone was derived. No evidence for the presence of expressed new genes could be obtained in this region by homology and EST searches or coding and exon prediction analyses. One TA microsatellite repeat spanning 2545 bases with as many as 913 repetitions was found on the centromeric side of theMICAgene and was indicated to be a potential hot spot for genetic recombination. The two segments of approximately 35 kb upstream of theMICAandMICBgenes showed high sequence homology (about 85%) to each other, suggesting that segmental genome duplication including theMICAandMICBgenes must have occurred during the evolution of the humanMHC.     

Regulatory role of a monomorphic determinant of HLA Class I antigens in T cell proliferation

The role of HLA Class I antigens in T cell proliferation was investigated by using the anti-HLA Class I monoclonal antibodies (MoAb) CR10-215, CR10-325, and CR11-115. MoAb CR10-215 and CR11-115 recognize the same (or spatially close) monomorphic determinant, which is distinct and spatially distant from that reacting with MoAb CR10-325. Addition of MoAb CR10-215 and CR11-115 to cultures of peripheral blood mononuclear cells stimulated with MoAb OKT3, MoAb Pan T2, PHA, or PPD inhibited cell proliferation. The blocking is specific in that the anti-HLA Class I MoAb CR10-325 and the Pan T MoAb Pan T1 had no effect on the proliferation. The inhibitory activity of MoAb CR10-215 and CR11-115 does not reflect i) toxic effects, ii) induction of suppressor cells and factors, iii) blocking of the binding of mitogens to lymphocytes, iv) inhibition of the production of interleukin 1 (IL 1) and interleukin 2 (IL 2), or v) function of IL 2 receptor. Anti-HLA Class I MoAb were able to inhibit the proliferation of purified, Tac-, T cells. The inhibited cells did not express Tac antigen, as assayed by direct immunofluorescence, with MoAb anti-Tac, but released a normal amount of IL 2 in culture medium. These results indicate that monomorphic determinants of the HLA Class I complex are involved in the regulation of T cell proliferation. The effect appears to occur at the level of IL 2 receptor expression.     

A Large-Scale Genetic Analysis Reveals a Strong Contribution of the HLA Class II Region to Giant Cell Arteritis Susceptibility

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10⁻⁴⁰, OR = 1.73). A multivariate model including class II amino acids of HLA-DRβ1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1^∗04. An omnibus test on polymorphic amino acid positions highlighted DRβ1 13 (p = 4.08 × 10⁻⁴³) and HLA-DQα1 47 (p = 4.02 × 10⁻⁴⁶), 56, and 76 (both p = 1.84 × 10⁻⁴⁵) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10⁻⁶, OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10⁻⁶, OR = 1.20), and REL (rs115674477, p = 1.10 × 10⁻⁵, OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function.     

Identification of Disease-Promoting HLA Class I and Protective Class II Modifiers in Japanese Patients with Familial Mediterranean Fever

ObjectivesThe genotype-phenotype correlation of MEFV remains unclear for the familial Mediterranean fever (FMF) patients, especially without canonical MEFV mutations in exon 10. The risk of FMF appeared to be under the influence of other factors in this case. The contribution of HLA polymorphisms to the risk of FMF was examined as strong candidates of modifier genes.MethodsGenotypes of HLA-B and -DRB1 loci were determined for 258 mutually unrelated Japanese FMF patients, who satisfied modified Tel-Hashomer criteria, and 299 healthy controls. The effects of carrier status were evaluated for the risk of FMF by odds ratio (OR). The HLA effects were also assessed for clinical forms of FMF, subsets of FMF with certain MEFV genotypes and responsiveness to colchicine treatment.ResultsThe carriers of B*39:01 were increased in the patients (OR = 3.25, p = 0.0012), whereas those of DRB1*15:02 were decreased (OR = 0.45, p = 0.00050), satisfying Bonferroni’s correction for multiple statistical tests (n = 28, p<0.00179). The protective effect of DRB1*15:02 was completely disappeared in the co-existence of B*40:01. The HLA effects were generally augmented in the patients without a canonical MEFV variant allele M694I, in accordance with the notion that the lower penetrance of the mutations is owing to the larger contribution of modifier genes in the pathogenesis, with a few exceptions. Further, 42.9% of 14 colchicine-resistant patients and 13.5% of 156 colchicine-responders possessed B*35:01 allele, giving OR of 4.82 (p = 0.0041).ConclusionsThe differential effects of HLA class I and class II polymorphisms were identified for Japanese FMF even in those with high-penetrance MEFV mutations.