Cloning the human major histocompatibility complex in YACs

Humanized transgenic mice expressing HLA DR4-DQ3 haplotype: reconstitution of phenotype and HLA-restricted T-cell responses

Many autoimmune conditions have close genetic linkages to particular human histocompatibility leukocyte antigen (HLA) class II genes. With the aim of establishing a murine model of autoimmune disease, we have generated an HLA DR4-DQ3 haplotype transgenic (Tg) mouse that expresses a 440-kb yeast artificial chromosome harbouring DRA, DRB1*040101, DRB4*010301, DQA1*030101, DQB1*0302 and all the internal regulatory segments. This Tg mouse line was crossed to human CD4 (hCD4) Tg mice and endogenous class II knockout mice (I-A(o/o) and I-E(o/o)) lines to generate a DR4-DQ3.hCD4.IAE(o/o) Tg line. The Tg DR and DQ molecules are expressed on the physiological cell types in these animals, i.e. on most B cells (>85%), dendritic cells (DCs) and macrophages but not on T cells, with levels of expression comparable with those of human B cells (where DR > DQ expression). The DR4/DQ3 transgenes fully reconstituted the CD4 T-cell compartment, in both the thymus and the periphery, and the analysis of the T-cell receptor repertoire in the Tg mice confirmed that these class II molecules were able to mediate thymic selection of a broad range of Vbeta families. HLA DR4- and DQ3-restricted T-cell responses were elicited following immunization with known T-cell determinants presented by these molecules. Furthermore, the DR4-DQ3-restricted CD4(+) T cells conferred protective antibody-mediated immunity against an otherwise lethal infection with Salmonella enterica var. typhimurium. These new DR4-DQ3 Tg mice should prove to be valuable tools for dissecting the importance of this class II haplotype in autoimmune disorders like rheumatoid arthritis.

Genome-wide transcriptome mapping analysis identifies organ-specific gene expression patterns along human chromosomes

The Human Genome Project has revealed that there about 32,000 protein-encoding genes, which are distributed throughout the genome. It is unclear, however, whether genes are distributed on the chromosomes according to patterns linked to organ specificity. To explore the relationship between genes actively transcribed in normal tissues and their chromosomal locations, we analyzed serial analysis of gene expression libraries of normal human liver, brain, breast, and colon tissues. Transcriptome mapping analysis revealed that transcriptional activity in each tissue varied according to the chromosomal domains, and a weak positive correlation was observed between transcription density and gene density. We identified six liver-related and five colon-related chromosomal domains highly transcribed in each tissue, whereas no brain-related or breast-related chromosomal domains were identified. Representative genes located on these chromosomal domains were associated with the function of each organ and were highly conserved in both mouse and rat genomes. These data revealed that the transcriptional activities of normal human tissues are well orchestrated at chromosomal levels, suggesting that highly expressed genes may share physical proximity.

Chromosomal distribution of the human cardiovascular transcriptome

On the basis of previous observations in chromosomes 21 and 22, we hypothesize that there is a tissue-specific organization of cardiovascular gene transcripts in the human genome. To examine the distribution of heart-derived transcripts, we assigned a nonredundant set of 4628 fetal and 3574 adult known and uncharacterized cardiovascular expressed-sequence tags (cvESTs) to 5-Mb chromosomal 'windows' on the basis of publicly available sequence mapping data. On a whole-genome level (36,617 genes), chromosome 17 (19.2% in fetal, 16.5% in adult) contained the highest proportion of cvESTs, whereas chromosome Y (2.0% in fetal and adult) contained the lowest. In total, 50 of the 639 windows contained a significantly higher proportion of cvESTs (P < 0.003) compared with the genome-wide cvEST gene density, particularly on gene-dense chromosomes (that is, 17, 19, 22) as opposed to gene-rich chromosomes (for example, 1, 2, 11). This report provides insight into a possible role for complex tissue-specific gene regulation in the human genome.

Molecular dynamics of MHC genesis unraveled by sequence analysis of the 1,796,938-bp HLA class I region

The intensely studied MHC has become the paradigm for understanding the architectural evolution of vertebrate multigene families. The 4-Mb human MHC (also known as the HLA complex) encodes genes critically involved in the immune response, graft rejection, and disease susceptibility. Here we report the continuous 1,796,938-bp genomic sequence of the HLA class I region, linking genes between MICB and HLA-F. A total of 127 genes or potentially coding sequences were recognized within the analyzed sequence, establishing a high gene density of one per every 14.1 kb. The identification of 758 microsatellite provides tools for high-resolution mapping of HLA class I-associated disease genes. Most importantly, we establish that the repeated duplication and subsequent diversification of a minimal building block, MIC-HCGIX-3.8-1-P5-HCGIV-HLA class I-HCGII, engendered the present-day MHC. That the currently nonessential HLA-F and MICE genes have acted as progenitors to today's immune-competent HLA-ABC and MICA/B genes provides experimental evidence for evolution by "birth and death," which has general relevance to our understanding of the evolutionary forces driving vertebrate multigene families.

Construction of a swine BAC library: application to the characterization and mapping of porcine type C endoviral elements

A porcine bacterial artificial chromosome (BAC) library was constructed using the pBeloBAC11 vector. It comprised 107,520 clones with an average insert size of 135 kb, representing an almost fivefold coverage of the swine haploid genome. Screening of the library allowed recovery of one to eight clones for 142 unique markers located all over the genome, while it failed for only one marker. About 4% chimeric clones were found. The library was also screened for the protease gene of type C porcine endoviral sequences (PERVs), and 62 clones were recovered, all but two of which contained one protease gene. We found 20 protease sequences (PERV-1 to PERV-20) which, despite differing by point mutations, were all coding sequences. The most frequent sequence, PERV-2, was 100% similar to a protease sequence expressed in the porcine PK-15 cell line. Most of the clones harbored envelope genes. Thirty-three BAC clones were mapped by fluorescence in situ hybridization to 22 distinct locations on 14 chromosomes, including the X and Y chromosomes. These overall results indicate that there is generally one PERV copy per integration site. Although PERV sequences were not tandemly arranged, clusters of integration sites were observed at positions 3p1.5 and 7p1.1. Southern blot experiments revealed 20-30 PERV copies in the Large White pig genome studied here, and variations in PERV content among pigs of different breeds were observed. In conclusion, this BAC collection represents a significant contribution to the swine large genomic DNA cloned insert resources and provides the first detailed map of PERV sequences in the swine genome. This work is the first step toward identification of potential active sites of PERV elements.

Sequence organisation of the class II region of the human MHC

We present the genomic organisation of the extended class II region of the human MHC. This initial sequence, which is nearing completion, spans about 1.2 Mbp and is at present a composite of more than one haplotype. The sequencing of single haplotypes is planned for the future. The current sequence encompasses all of the known class II genes at the DP, DO, DM, DQ and DR loci as well as the transporter associated with antigen processing (TAP)/low molecular weight protein (LMP) antigen processing genes and the Tapasin locus, at the extended centromeric end.

Genome sequencing analysis of the 1.8 Mb entire human MHC class I region

The human MHC class I region spans 1.8 Mb from the MICB gene to the HLA-F gene at the telomeric end of the HLA region. There are fewer genes recognized in this region than in the class II or class III region, probably because this region remained uncharacterized for genomic organization. Based on the 1,796,938 bp genomic sequence of the entire class I region determined in our laboratory, the complete gene structure of this region has finally emerged. This region embraces as many as 118 genes (73 known and 45 new genes) with a gene density of one gene every 15.2 kb, which is comparable to that of the gene-rich class III region. The GC content is fairly uniform throughout the class I region, being 45.8% on average, which corresponds to the isochore H1. By investigation of genetic polymorphisms in 26 out of 758 microsatellite repeats identified in the class I region, we could reduce the critical region for Behçet's disease (associated with B51) and psoriasis vulgaris (associated with Cw6) to approximately 50 kb segments, between MICA and HLA-B and between TCF19 and S, respectively. Thus, systematic large-scale genomic sequencing provides an efficient way of identifying genes and of mapping disease-susceptible genes in the genome.

The human major histocompatibility complex: 42,221 bp of genomic sequence, high-density sequence-tagged site map, evolution, and polymorphism for HLA class I

We report the isolation and characterization of newly identified yeast artificial chromosome (YAC) and bacterial artificial chromosome (BAC) clones spanning the HLA class I region between HLA-C and HLA-E and of YACs extending telomeric of HLA-F. When included with previously characterized HLA class I YACs, a contiguous stretch of over 2.4 Mb pairs including the entire class I region has been isolated as a series of overlapping YAC and BAC clones. Evidence that the cloned DNA faithfully represents the source genomic DNA was obtained by extensive characterization of the YACs and by independent isolation of two or more overlapping YACs or BACs spanning the entire region. As a result of this work, over 80 unique sequence probes were identified, the majority of which were sequenced to yield 42,221 bp of new major histocompatibility complex (MHC)-derived sequence. Some of these data were reduced to sequenced tagged site primer sets, facilitating the isolation of all or nearly all of HLA class I from a variety of genomic libraries. The sequence data were analyzed for protein coding capacity and homology to existing expressed tagged sites and tested for conservation of sequences in other mammalian genomes. These results indicated that large portions of the HLA class I region are conserved among mammals. Measurements of polymorphism within non-HLA class I loci generated additional data pointing toward information of potential relevance to MHC-associated diseases. The combined data and clones presented here set the stage for the determination of the complete nucleotide sequence of HLA class I.

Genomic analysis of the Tapasin gene, located close to the TAP loci in the MHC

The Tapasin molecule is a member of the immunoglobulin (Ig) superfamily required for the association of TAP transporters and MHC class I heterodimers in the endoplasmic reticulum. In this study, the Tapasin gene was precisely mapped in relation to the MHC. The gene was centromeric of the HLA-DP locus between the HSET and HKE1.5 genes and within 500 kbp of the TAP1 and TAP2 genes. A homologous mouse EST was mapped to a syntenic position on chromosome 17, centromeric of the H-2 K locus. Similarly, the rat Tapasin gene was shown to be in an equivalent location with respect to the RT1.A locus. The localization of Tapasin, TAP, LMP and class I genes within such a short distance of each other on the chromosome implies some regulatory or functional significance. We determined the Tapasin gene sequence for comparison of its structure to that of other Ig superfamily members, such as MHC class I genes. The IgC domain was encoded by a separate exon. However, the positions of the other introns were not characteristic of other Ig superfamily genes, indicating that Tapasin has a distinct phylogeny.

Physical mapping of the mouse genome

This review is intended to provide an overview of techniques and a source of reagents for physical mapping of the mouse genome. It focuses on those applications, methods, or resources unique to the mouse and on the generation of comparative physical maps. The reference list is not comprehensive; rather, recent reviews on each topic and selected representative examples are given.

Nucleotide sequencing analysis of the 146-kilobase segment around the IkBL and MICA genes 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 of HLA class 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 the IkBL and MICA loci 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 the IkBL gene to exon 6 of MICA. This region was confirmed to contain five known genes, IkBL, BAT1, MICB, P5-1, and HLA-X (class I fragment), from centromere to telomere, and their exon-intron organizations were determined. The 3.8-1 homologue gene (3.8-1-hom) showing 99.7% identity with the 3.8-1 cDNA clone, which was originally isolated using the 3.8-kb EcoRI fragment between the HLA-54/H and the HLA-G genes, was detected between MICA and MICB and was suggested to represent the cognate 3.8-1 genomic 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 the MICA gene and was indicated to be a potential hot spot for genetic recombination. The two segments of approximately 35 kb upstream of the MICA and MICB genes showed high sequence homology (about 85%) to each other, suggesting that segmental genome duplication including the MICA and MICB genes must have occurred during the evolution of the human MHC.

Haplotype analysis of hemochromatosis: evaluation of different linkage-disequilibrium approaches and evolution of disease chromosomes

We applied several types of linkage-disequilibrium calculations to analyze the hereditary hemochromatosis (hh) locus. Twenty-four polymorphic markers in the major histocompatibility complex (MHC) class I region were used to haplotype hh and normal chromosomes. A total of 169 hh and 161 normal chromosomes were analyzed. Disequilibrium values were found to be high over an unusually large region beginning 150 kb centromeric of HLA-A and extending nearly 5 Mb telomeric of it. Recombination in this region was approximately 28% of the expected value. This low level of recombination contributes to the unusually broad region of linkage disequilibrium found with hh. The strongest disequilibrium was found at locus HLA-H (delta = .84) and at locus D6S2239 (delta = .85), a marker approximately 10 kb telomeric to HLA-H. All disequilibrium methods employed in this study found peak disequilibrium at HLA-H or D6S2239. The cys282tyr mutation in HLA-H, a candidate gene for hh, was found in 85% of disease chromosomes. A haplotype phylogeny for hh chromosomes was constructed and suggests that the mutation associated with the most common haplotype occurred relatively recently. The age of the hh mutation was estimated to be approximately 60-70 generations. Disequilibrium was maintained over a greater distance for hh-carrying chromosomes, consistent with a recent mutation for hh. Our data provide a reasonable explanation for previous difficulties in localizing the hh locus and provide an evolutionary history for disease chromosomes.

Clone-contig and STS maps of the hereditary hemochromatosis region on human chromosome 6p21.3-p22

YAC-based and bacterial-clone based STS-content maps were constructed that served as the framework physical maps for the positional cloning of a candidate gene for hereditary hemochromatosis. The YAC-based map comprises 43 YACs and 86 STS and spans approximately 8 Mb of DNA between the class I region of the major histocompatibility complex on human chromosome 6p21.3 and D6S276 in 6p22. Comparison with published maps revealed a hole in the MIT/Whitehead and CEPH YAC maps that includes the immediate region around the hemochromatosis gene itself. Approximately 3 Mb of DNA was covered by a bacterial clone contig that consists of 38 BACs, 45 PACs, 26 PI clones and one lambda phage. The bacterial clone-based STS map comprises 153 STSs. A contiguous block of 8 STSs could be amplified from both human chromosome 6 and 5. Further characterization of selected STSs and bacterial clones by radiation hybrid mapping and fluorescence in situ hybridization, respectively, revealed the presence of a multicopy DNA segment, more than one bacterial clone length in size, which is duplicated near the chromosome-6 centromere and part of which is present in multiple copies on chromosome 5. Possible implications of the incomplete public YAC-contig map and of the multicopy segment for physical mapping and linkage disequilibrium studies of the hemochromatosis candidate region are discussed.

Construction of a swine YAC library allowing an efficient recovery of unique and centromeric repeated sequences

A swine DNA genomic library was constructed in yeast artificial chromosome (YAC) using the pYAC4 vector and the AB1380 strain. The DNA prepared from two Large White males was partially digested with EcoRI and size selected after both digestion and ligation. The YAC library contained 33792 arrayed clones with an average size of 280 kb as estimated by analysis of 2% of the clones, thus representing a threefold coverage of the swine haploid genome. The library was organized in pools to facilitate the PCR screening. The complexity of the library was tested both for unique and centromeric repeated sequences. In all, 20 out of 22 primer sets allowed the characterization of one to six clones containing specific unique sequences. These sequences are known to be on Chromosomes (Chrs) 1, 2, 5, 6, 7, 8, 13, 14, 15, 17, and X. Eight additional clones carrying centromeric repeat units were also isolated with a single primer set. The sequencing of 37 distinct repeat units of about 340 bp subcloned from these eight YACs revealed high sequence diversity indicating the existence of numerous centromeric repeat unit subfamilies in swine. Furthermore, the analysis of the restriction patterns with selected enzymes suggested a higher order organization of the repeat units. According to preliminary FISH experiments on a small number of randomly chosen YACs and YACs carrying specific sequences, the chimerism appeared to be low. In addition, primed in situ labeling experiments favored the idea that the YACs with centromeric repeat sequences were derived from a subset of metacentric and submetacentric chromosomes.

The gene G13 in the class III region of the human MHC encodes a potential DNA-binding protein

G13 is a single-copy gene lying approx. 75 kb centromeric of the complement gene cluster in the class III region of the human MHC. The gene spans approx. 17 kb of DNA and has been shown to encode mRNA of approx. 2.7 kb that is present in cell lines representing lymphoid and non-lymphoid tissues, indicating that it is ubiquitously expressed. The complete nucleotide sequence of the 2.7 kb mRNA has been derived from cDNA and genomic clones. The longest open reading frame obtained for G13 codes for a 703 amino acid protein of approx. 77 kDa in molecular mass. Comparison of the putative G13 amino acid sequence with the protein databases revealed significant similarities with DNA-binding proteins of the leucine zipper class, including a human cAMP response element binding protein. G13 contains a bZIP motif, a region rich in basic amino acids adjacent to a coiled-coil leucine zipper domain, common to this class of proteins that is known to be involved in dimerization and DNA binding. Antibodies raised against a fragment encoding the C-terminal half of the putative G13 protein recognized a major polypeptide of approx. 86 kDa and a minor polypeptide of approx. 78 kDa on immunoblotting of U937 cell extracts; this has been confirmed by immunoprecipitation experiments. Even though it contained at least one potential bipartite nuclear localization signal, the G13 protein was present both in the cytoplasm and the nucleus of the fibroblast cells. Thus G13 might be a novel DNA-binding protein that is perhaps translocated to the nucleus in a regulated manner.

A new non-HLA multigene family associated with the PERB11 family within the MHC class I region

In an effort to initiate steps designed to characterize the idiopathic hemochromatosis disease gene, the HLA-A/HLA-F region where this gene is in disequilibrium linkage with some polymorphic markers has been overlapped by a yeast artificial chromosome (YAC) contig. In order to achieve the physical mapping of these YACs and of the corresponding genomic region, we subcloned one of the YACs involved. A computer-assisted analysis of the sequence of one subclone led to the isolation of a potential exon that proved to belong to a new expressed messenger named HCGIX. After Southern blot analysis, the corresponding cDNA clone was found to belong to a new multigene family whose members are dispersed throughout the HLA class I region and are closely associated with members of another recently described multigene family designated PERB11. The data reported here suggest that these two multigene families form a cluster that have been dispersed together throughout the telomeric part of the major histocompatibility complex and have been involved in the genesis of this human class I region.

Physical map of the HLA-A/HLA-F subregion and identification of two new coding sequences

As part of an effort to characterize the hemochromatosis gene, we selected three non-chimeric yeast artificial chromosomes (YACs) overlapping with the YAC B30 previously described and forming an 800 kilobase contig covering the HLA-A/HLA-F region. The precise physical map of these YACs and of the corresponding genomic region were established. Nine concentrated sites of CpG cutter elements, potentially HTF islands, were mapped. In addition, several probes have been generated as tools for mapping and examining transcripts produced in the region. This allowed for the characterization and localization of two new coding sequences, provisionally named HCG (for hemochromatosis candidate gene) and numbered VIII and IX.

Allelic association of microsatellites of 6p in Italian hemochromatosis patients

Hemochromatosis (HC) is an inherited disorder of iron metabolism and is frequently seen in Caucasians. The biochemical defect and the responsible gene are unknown, but the HC locus is closely linked to HLA-A on human chromosome 6 in the region 6p21.3. Although extensive studies have been performed in several populations, the precise location of the gene is still undefined. Linkage disequilibrium with HC has been detected for loci that are 3 cM apart: HLA class I and D6S105, which is located on the telomeric side of HLA-A. We have analyzed the inheritance of several multi-allele polymorphisms that map to 6p (D6S265, Y52, HLA-F, D6S306, D6S105, D6S464, D6S299) in 34 Italian HC families and in 17 unrelated patients. Significant association with HC was shown for alleles of multiple markers in the HLA-A region, for the distant marker D6S105, but not for the D6S299 marker at 4 cM from HLA-A on the telomeric side. HC status was unambiguously assigned to 70 affected and 63 unaffected chromosomes from family studies. Thirty five different haplotypes were found in 70 HC chromosomes when considering four markers most tightly associated with the disease. A predominant haplotype comprising alleles 1-3-1-8 (marker order D6S265, HLA-A, Y52, D6S105) accounted for 30% of the HC chromosomes and was absent in normals. A minority of other HC haplotypes could be related to the major haplotype by assuming single crossover events. Results of haplotype studies suggest a founder effect in the Italian population, as previously shown in Australian patients, and a possible common mutation shared with affected individuals of Celtic origin.

A long-range physical map of human chromosome 21q22.1 band from the YAC continuum

The human Chromosome (Chr) 21q22.1 region contains several genes for cytokines and neurotransmitters and the gene for superoxide dismutase (mutant forms of which can cause familial amyotrophic lateral sclerosis). A region of approximately 5.8 Mb encompassing D21S82 and the glycinamide ribonucleotide transformylase (GART) loci was covered by overlapping YAC clones, which were contiguously ordered by clone walking with sequence-tagged site (STSs). A total of 76 markers, including 29 YAC end-specific STSs, were unambiguously ordered in this 5.8-Mb region, and the average interval between markers was 76 kb. Restriction maps of the YAC clones with rare-cutting enzymes were simultaneously prepared, and the restriction sites were aligned to obtain a consensus restriction map of the proximal region of the 21q22.1 band. The restriction map made from 44 overlapping YACs contains 54 physically assigned STSs. By integrating the consensus map of the adjacent 1.8-Mb region, we obtained a fine physical map spanning 6.5 Mb of human Chr 21q22.1. This map contains 24 precisely positioned end-specific STSs and 12 NotI-linking markers. More than 39 potential CpG islands were identified in this region and were found to be unevenly distributed. This physical map and the YACs should be useful as a reference map and as a resource for further structural analysis of the Giemsa-negative band (R-band) of Chr 21q22.1.

A new highly polymorphic marker in the 5' untranslated region of HLA-F shows strong allelic association with haemochromatosis

The 5' untranslated region of HLA-F contains a polypurine tract comprising repeats of tri- and hexa-nucleotide motifs. We have recently demonstrated that this polypurine tract is highly polymorphic by using the polymerase chain reaction. Here, we demonstrate that some of the alleles can be explained by a deletion of approximately 100 bp DNA and show that alleles of this novel, highly polymorphic locus are as strongly associated with haemochromatosis as HLA-A3 or D6S105-8. The observed frequency of heterozygosity at HLA-RF is extremely high (95%) and this locus has been found to be informative in pedigrees that are non-informative at HLA-A and D6S105. We also show an example of replication slippage at HLA-F in one pedigree.

1.8-megabases fine physical map encompassing IFNAR and AML1 loci on human chromosome 21q22.1

A long-range restriction map of the 1.8-megabases (mb) region encompassing the area between the interferon-alpha receptor and the acute myelogenous leukemia loci on human chromosome 21q22.1 was constructed after analysis of both the contiguous yeast artificial chromosome (YAC) clones and genomic DNA. Analysis of pulsed-field gel electrophoresis of lymphoblastoid DNA digested with three rare-cutting enzymes, Not I, Mlu I, and Nru I, revealed the positions of 17 markers on each restriction map. The 1.8-mb YAC contig that covers this region was obtained through YAC walking mediated by sequence-tagged sites (STSs), with 29 STSs including 12 newly generated YAC end-specific STSs. The consensus restriction map from 15 overlapping YACs and the positioning of the STS markers on each clone allowed 24 markers including 4 Not I-linking STSs to be ordered and mapped physically. Comparison of the maps revealed that the proximal region contains more unmethylated CpG islands than the distal region, which suggests that many expressed genes are in the proximal region. This fine consensus physical map will be informative and useful for construction of contigs of cosmid, P1, or BAC clones for further large-scale sequencing in this gene-rich region.

Genomic organization of a mouse MHC class II region including the H2-M and Lmp2 loci

The region encompassing the Ma, Mb1, Mb2, and Lmp2 genes of the mouse class II major histocompatibility complex (MHC) was sequenced. Since this region contains clusters of genes required for efficient class I and class II antigen presentation, it was interesting to search for putative additional genes in the 21 kilobase gap between the Mb1 and Lmp2 genes. Computer predictions of coding regions and CpG islands, exon trapping experiments, and cross-species comparison with the corresponding human sequence indicate that no additional functional gene is present in that stretch. However, computer analysis revealed the possible existence of an alternative 3' exon for Mb1. Except for the fact that the mouse MHC contains two Mb genes, the genomic organization of the H2-M loci was found to be almost identical to the organization of the human HLA-DM genes. The promoter regions of the Ma and Mb genes also resemble classical class II promoters, containing typical S, X, and Y boxes. Like the human genes, the three H2-M genes displayed very limited polymorphism when we compared the cDNA sequences from six haplotypes. Finally, comparison of DMB with Mb1 and Mb2, both at the genomic level and in their coding regions, suggests that the Mb gene was recently duplicated, probably only in certain rodents.

Allelic associations and homozygosity at loci from HLA-B to D6S299 in genetic haemochromatosis

Haemochromatosis (GH) is an autosomal recessive disorder in which increased iron absorption causes iron overload. The gene (HFE) is closely linked to HLA-A on chromosome 6 (6p21.3) but has not yet been identified. We have examined eight polymorphic loci, HLA-B (most centromeric), I82, D6S265, HLA-A, D6S128, HLA-F, D6S105, and D6S299 (most telomeric) in 37 unrelated patients and 60 control subjects. There are also significant positive associations between GH and alleles at all loci except D6S299. Analysis of 48 GH chromosomes in which haplotypes could be established showed that the most common haplotype was I82-2:D6S265-1:HLA-A3:D6S128-2:HLA-F1:D6S105-8. This was present in 28 of 48 chromosomes. In 14 the haplotype included HLA-B7 but only in seven did this extend beyond the telomere to D6S299-2 (the most common allele on GH chromosomes at this locus). In 36 out of 48 chromosomes the two locus haplotype, F1:D6S105-8 was present. Since haemochromatosis appears to originate from a founder mutation we have examined linkage disequilibrium between these various loci and GH using calculations of pexcess. The maximum value (0.72, 95% CI 0.55-0.85) is given by D6S105-8 but is not significantly different from values for HLA-A3 and HLA-F1 (0.50, 95% CI 0.34-0.61 and 0.49, 0.25-0.66 respectively). However, both HLA-A and D6S105 give a value for pexcess which is significantly higher than that for the most centromeric marker, HLA-B (0.17, 95% CI 0.02-0.30). We have counted the number of patients who are homozygous for the common allele at each locus. At D6S105, 22 patients are homozygous for allele 8, with 18 homozygous for HLA-F1 and 10 homozygous for A3. The pattern of cumulative homozygosity suggests a gene location closer to D6S105 than HLA-A. We have also analysed our data for divergence from the apparent founder haplotype (A3:F1:105-8) and have calculated the theoretical frequencies of crossovers between loci. These data suggest a location telomeric to D6S105. A more precise localisation of the gene may be possible with the identification of new markers around D6S105.

A new marker in the HLA class I region is associated with the age at onset of IDDM

The (MHC) class II association with insulin-dependent diabetes mellitus (IDDM) is well documented. However, it is likely that genes within the MHC class III and the class I region also play a role in determining susceptibility to IDDM. In this study we have used a novel molecular probe to investigate the class I P3A and P3B loci of 179 patients with IDDM and 142 normal control subjects. A highly significant increase in the frequency of the class I P3 4.0;1.5 kilobase (kb) and 4.0;1.8;1.5 kb genotypes was found in patients compared to the control subjects (chi 2 46.8, 6 df, p < 0.0001). The association with the P3B 1.5 kb allele was strongly associated with the age at onset of diabetes, being present in 96.2% of subjects who developed diabetes between the age of 10-20 years compared to 55.0 and 74.6% who developed diabetes before 10 years or after 20 years, respectively (chi 2 31.4, p < 0.0001). There was no evidence for linkage disequilibrium between the DQA1 and DQB1 loci and P3B suggesting that this is an independent association. In conclusion, these results suggest that genes in both the MHC class I and II regions confer susceptibility to IDDM and are related to the age at onset of the disease.

Structure and evolution of a member of a new subfamily of GTP-binding proteins mapping to the human MHC class I region

A gene coding for a putative GTP-binding protein, MMR1, has been localized on band C of the murine Chr 17 within or close to the MHC (Denizot et al. 1992). Its human homolog, HSR1, localized to the human MHC class I region, is described in this paper. Its sequence, compared with MMR1, shows that the conceptual proteins encoded by these genes are highly homologous and have thus been subjected to high constraints during evolution. Furthermore, a detailed databank search with HSR1 leads to the characterization of a new subfamily of GTP-binding proteins, of which HSR1 and MMR1 are the only eukaryotic members. The precise localization of HSR1 within the human MHC class I region is also presented.