Hereditary haemochromatosis (HH), which affects some 1 in 400 and has an estimated carrier frequency of 1 in 10 individuals of Northern European descent, results in multi-organ dysfunction caused by increased iron deposition, and is treatable if detected early. Using linkage-disequilibrium and full haplotype analysis, we have identified a 250-kilobase region more than 3 megabases telomeric of the major histocompatibility complex (MHC) that is identical-by-descent in 85% of patient chromosomes. Within this region, we have identified a gene related to the MHC class I family, termed HLA-H, containing two missense alterations. One of these is predicted to inactivate this class of proteins and was found homozygous in 83% of 178 patients. A role of this gene in haemochromatosis is supported by the frequency and nature of the major mutation and prior studies implicating MHC class I-like proteins in iron metabolism.

In the process of positionally cloning a candidate gene responsible for hereditary hemochromatosis (HH), we constructed a 1.1-Mb transcript map of the region of human chromosome 6p that lies 4.5 Mb telomeric to HLA-A. A combination of three gene-finding techniques, direct cDNA selection, exon trapping, and sample sequencing, were used initially for a saturation screening of the 1.1-Mb region for expressed sequence fragments. As genetic analysis further narrowed the HH candidate locus, we sequenced completely 0.25 Mb of genomic DNA as a final measure to identify all genes. Besides the novel MHC class 1-like HH candidate gene HLA-H, we identified a family of five butyrophilin-related sequences, two genes with structural similarity to a type 1 sodium phosphate transporter, 12 novel histone genes, and a gene we named RoRet based on its strong similarity to the 52-kD Ro/SSA lupus and Sjogren’s syndrome auto-antigen and the RET finger protein. Several members of the butyrophilin family and the RoRet gene share an exon of common evolutionary origin called B30-2. The B30-2 exon was originally isolated from the HLA class 1 region, yet has apparently “shuffled” into several genes along the chromosome telomeric to the MHC. The conservation of the B30-2 exon in several novel genes and the previously described amino acid homology of HLA-H to MHC class 1 molecules provide further support that this gene-rich region of 6p21.3 is related to the MHC. Finally, we performed an analysis of the four approaches for gene finding and conclude that direct selection provides the most effective probes for cDNA screening, and that as much as 30% of ESTs in this 1.1-Mb region may be derived from noncoding genomic DNA.[The sequence data described in this paper have been submitted to GenBank under accession nos. U90543–U90548,U90550–U90552, and U91328.]

Hereditary hemochromatosis is a recessive disease of iron metabolism widely distributed among people of European descent. Most patients have inherited the causative mutation from a single ancestor. In the course of cloning the hemochromatosis gene, genotypes were generated for these samples at 43 microsatellite repeat markers that span the 6.5-Mb hemochromatosis gene region. The data used to reconstruct the ancestral haplotype across the hemochromatosis gene region are presented in this paper. Portions of the ancestral haplotype were present on 85% of patient chromosomes in this sample and ranged in size from approximately 500 kb to greater than 6.5 Mb. Only one marker, D6S2239, was identical by descent on all of the patient chromosomes containing the ancestral mutation. In contrast, only 3 of the 128 control chromosomes, or 2.3%, carried the ancestral mutation and the surrounding ancestral haplotype. To test new methods for gene finding using linkage disequilibrium we analyzed the genotypic data with a multilocus maximum likelihood method (DISMULT) and a single point method (DISLAMB), both written to analyze data generated from multi-allelic markers. The maximum value from DISLAMB analysis occurred at marker D6S2239, which is less than 20 kb from the hemochromatosis gene HFE, consistent with the haplotype analysis. The peak of the multi-point analysis was 700 kb from HFE, possibly due to the nonuniform recombination rates within this large region. The recombination rate appears to be lower than expected centromeric of the HFE gene.

No abstract