Localization of new genes and markers to the distal part of the human major histocompatibility complex (MHC) region and comparison with the mouse: new insights into the evolution of mammalian genomes

RNA-sequencing study of peripheral blood monocytes in chronic periodontitis

BACKGROUND

Monocytes are an important cell type in chronic periodontitis (CP) by interacting with oral bacteria and mediating host immune response. The aim of this study was to reveal new functional genes and pathways for CP at monocyte transcriptomic level.

METHODS

We performed an RNA-sequencing (RNA-seq) study of peripheral blood monocytes (PBMs) in 5 non-smoking moderate to severe CP (case) individuals vs. 5 controls. We took advantage of a microarray study of periodontitis to support our findings. We also performed pathway-based analysis on the identified differentially expressed (DEx) transcripts/isoforms using DAVID (Database for Annotation, Visualization and Integrated Discovery).

RESULTS

Through differential expression analyses at both whole gene (or whole non-coding RNA) and isoform levels, we identified 380 DEx transcripts and 5955 DEx isoforms with a PPEE (posterior probability of equal expression) of <0.05. Pervasive up-regulation of transcripts at isoform level in CP vs. control individuals was observed, suggesting a more functionally active monocyte transcriptome for CP. By comparing with the microarray dataset, we identified several CP-associated novel genes (e.g., FACR and CUX1) that have functions to interact with invading microorganisms or enhance TNF production on lipopolysaccharide stimulation. DAVID analysis of both the RNA-seq and the microarray datasets leads to converging evidence supporting "endocytosis", "cytokine production" and "apoptosis" as significant biological processes in CP.

CONCLUSIONS

As the first RNA-seq study of PBMs for CP, this study provided novel findings at both gene (e.g., FCAR and CUX1) and biological process level. The findings will contribute to better understanding of CP disease mechanisms.

Genetic study of the myelin oligodendrocyte glycoprotein (MOG) gene in schizophrenia

Schizophrenia (SCZ) is a neuropsychiatric disorder that affects approximately 1% of the general population. The human leukocyte antigen (HLA) system has been implicated in several genetic studies of SCZ. The myelin oligodendrocyte glycoprotein (MOG) gene, which is located close to the HLA region, is considered a candidate for SCZ due to its association with white matter abnormalities and its importance in mediating the complement cascade. Four polymorphisms in the MOG gene (CA)n (TAAA)n, and two intronic polymorphisms, C1334T and C10991T, were investigated for the possibility of association with SCZ using 111 SCZ proband and their families. We examined the transmission of the alleles of each of these polymorphisms with the transmission disequilibrium test. We did not observe significant evidence for biased transmission of alleles at the (CA)n (chi2=2.430, 6 df, P=0.876) (TAAA)n (chi2=3.550, 5 df, P=0.616), C1334T (chi2=0.040, 1 df, P=0.841) and C10991T (chi2=0.154, 1 df, P=0.695) polymorphisms. Overall haplotype analysis using the TRANSMIT program was also not significant (chi2=7.954, 9 df, P=0.539). Furthermore, our results comparing mean age at onset in the genotype groups using the Kruskal-Wallis Test were not significant. Our case-control analyses (182 cases age-, sex- and ethnicity-matched with healthy controls) and combined z-score [(CA)n: z-score=-1.126, P=0.130; (TAAA)n: z-score=-0.233, P=0.408; C1334T: z-score=0.703, P=0.241; C10991T: z-score=0.551, P=0.291] were also not significant. Although our data are negative, the intriguing hypothesis for MOG in SCZ may warrant further investigation of this gene.

Co-duplication of olfactory receptor and MHC class I genes in the mouse major histocompatibility complex

We report the 897 kb sequence of a cluster of olfactory receptor (OR) genes located at the distal end of the major histocompatibility complex (MHC) class I region on mouse chromosome 17 of strain 129/SvJ (H2bc). With additional information from the mouse genome draft sequence, we identified 59 OR loci (approximately 20% pseudogenes) in contrast to only 25 OR loci (approximately 50% pseudogenes) in the corresponding centromeric OR cluster that is part of the 'extended MHC class I region' on human chromosome 6. Comparative analysis leads to three major observations: (i) most of the OR subfamilies have evolved independently in the two species, expanding more in the mouse, and resulting in co-orthologs--subfamilies of highly similar paralogs that keep orthologous relationships with their human counterparts; (ii) three of the mouse OR subfamilies have no orthologs in humans; and (iii) MHC class I loci are interspersed in the OR cluster in mouse but not in human, and were subjected to co-duplication with OR genes. Screening of our sequence against the available sequences of other strains/haplotypes revealed that most of the OR loci are polymorphic and that the number of OR loci may vary among strains/haplotypes. Our findings that MHC-linked OR loci share duplication with MHC class I loci, have duplicated extensively and are polymorphic revives questions about potential reciprocal influences acting on the dynamics and evolution of the H2 region and the H2-linked OR loci.

Expression, genomic structure and mapping of the thymus specific protease prss16: a candidate gene for insulin dependent diabetes mellitus susceptibility

PRSS16 is a serine protease specifically expressed by epithelial cells in the thymic cortex. The human gene is encoded on 6p21.3-p22 where recent linkage analysis has identified an association with insulin dependent diabetes mellitus (IDDM) susceptibility independent of HLA-DR3. To further investigate its potential role in autoimmunity, we characterized the mouse orthologue, Prss16. The genomic structure of Prss16 shows conservation with the human gene in size, number of exons and chromosomal location. Mapping of Prss16 places it on mouse chromosome 13 centromeric of thesatin locus. This region is comparable to the PRSS16 region on human chromosome 6 and has also been linked to quantitative trait locus for IDDM in the nonobese diabetic mouse. Similar to the human gene, Prss16 expression is highly specific in the mouse with expression limited to the cortical thymic epithelium. Notably, embryonic expression coincides with population of the thymic anlage with T-cell precursors and initiation of T-cell development. We also show that NOD and New Zealand Black mice, which have a disrupted thymic architecture and autoimmune phenotype, have lower levels of Prss16 expression compared to C57BL/6 mice. These findings support the role of Prss16 in T-cell development and susceptibility to autoimmunity in the mouse.

Hemochromatosis gene in leukemia and lymphoma

The gene causing hereditary hemochromatosis (HH), HFE is an HLA class I-like gene with no known immunological function but indirectly related to the immune functions because of its role in iron transport. It is located 6.5 Mb telomeric to HLA-A. The most common mutation of HFE, C282Y, has a Celtic origin and most patients with HH are homozygous for it in Northern European populations. While there is an enormously increased risk for hepatocellular cancer in hemochromatosis that is attributed to the toxic effects of iron, the risk for extra-hepatic cancers is also increased slightly. Recent studies have found genetic associations between several cancers and C282Y but only in the presence of a particular allele of the transferrin receptor gene. This suggests that the increased cancer risk is more likely due to the effects of iron. In childhood acute lymphoblastic leukemia (ALL), however, there is a strong association of C282Y with a gender effect in two different Celtic populations. This association does not require homozygosity for C282Y or an interaction with the transferrin receptor gene, and is male-specific. The other HFE mutation H63D does not confer increased risk to childhood ALL. Acute myeloblastic leukemia and Hodgkin's disease in adults do not have an association with HFE. Its male-specificity, occurrence in childhood and the lack of a gene-dosage effect suggest that the C282Y association in childhood ALL may reflect the involvement of another HLA-linked gene in leukemia susceptibility.

Immunopathogenic and clinical relevance of antibodies against myelin oligodendrocyte glycoprotein (MOG) in Multiple Sclerosis

Recent neuropathological findings identified four distinct immunopathogenic pathways of demyelination and tissue destruction in the most common inflammatory demyelinating central nervous system disorder, Multiple Sclerosis. One of this neuropathological subtypes is characterised by features of antibody-mediated demyelination. A role of anti-myelin antibodies in the disease evolution of multiple sclerosis has been suggested already for a long time, however, their pathogenetic and clinical relevance is not understood yet. This present article will discuss recently published and some preliminary data on the immunopathogenic role of antibodies against myelin oligodendrocyte glycoprotein (MOG) and other myelin/nonmyelin targets in multiple sclerosis, as well as possible clinical implications for prognosis and therapy in the future.

Comparative genomic aspects of rat, mouse and human MHC class I gene regions

In this review a particular aspect of the genomic structure of the major histocompatibility complex (MHC), the organization of MHC class I regions, will be discussed for the rat in comparison to mouse and human.

Characterization of clustered MHC-linked olfactory receptor genes in human and mouse

Olfactory receptor (OR) loci frequently cluster and are present on most human chromosomes. They are members of the seven transmembrane receptor (7-TM) superfamily and, as such, are part of one of the largest mammalian multigene families, with an estimated copy number of up to 1000 ORs per haploid genome. As their name implies, ORs are known to be involved in the perception of odors and possibly also in other, nonolfaction-related, functions. Here, we report the characterization of ORs that are part of the MHC-linked OR clusters in human and mouse (partial sequence only). These clusters are of particular interest because of their possible involvement in olfaction-driven mate selection. In total, we describe 50 novel OR loci (36 human, 14 murine), making the human MHC-linked cluster the largest sequenced OR cluster in any organism so far. Comparative and phylogenetic analyses confirm the cluster to be MHC-linked but divergent in both species and allow the identification of at least one ortholog that will be useful for future regulatory and functional studies. Quantitative feature analysis shows clear evidence of duplications of blocks of OR genes and reveals the entire cluster to have a genomic environment that is very different from its neighboring regions. Based on in silico transcript analysis, we also present evidence of extensive long-distance splicing in the 5'-untranslated regions and, for the first time, of alternative splicing within the single coding exon of ORs. Taken together with our previous finding that ORs are also polymorphic, the presented data indicate that the expression, function, and evolution of these interesting genes might be more complex than previously thought.

Characterization of Clustered MHC-Linked Olfactory Receptor Genes in Human and Mouse

Olfactory receptor (OR) loci frequently cluster and are present on most human chromosomes. They are members of the seven transmembrane receptor (7-TM) superfamily and, as such, are part of one of the largest mammalian multigene families, with an estimated copy number of up to 1000 ORs per haploid genome. As their name implies, ORs are known to be involved in the perception of odors and possibly also in other, nonolfaction-related, functions. Here, we report the characterization of ORs that are part of the MHC-linked OR clusters in human and mouse (partial sequence only). These clusters are of particular interest because of their possible involvement in olfaction-driven mate selection. In total, we describe 50 novel OR loci (36 human, 14 murine), making the human MHC-linked cluster the largest sequenced OR cluster in any organism so far. Comparative and phylogenetic analyses confirm the cluster to be MHC-linked but divergent in both species and allow the identification of at least one ortholog that will be useful for future regulatory and functional studies. Quantitative feature analysis shows clear evidence of duplications of blocks of OR genes and reveals the entire cluster to have a genomic environment that is very different from its neighboring regions. Based on in silico transcript analysis, we also present evidence of extensive long-distance splicing in the 5′-untranslated regions and, for the first time, of alternative splicing within the single coding exon of ORs. Taken together with our previous finding that ORs are also polymorphic, the presented data indicate that the expression, function, and evolution of these interesting genes might be more complex than previously thought.

[The sequence data described in this paper have been submitted to the EMBL nucleotide data library under accession nos.Z84475, Z98744, Z98745, AL021807, AL021808, AL022723, AL022727,AL031893, AL035402, AL035542, AL050328, AL050339, AL078630, AL096770,AL121944, AL133160, and AL133267.]

Comparative mapping of two adjacent regions of MMU19 with their human counterpart on HSA11q13

High resolution physical maps of two adjacent regions of MMU19 were constructed in order to establish a comparative map between the pericentromeric region of MMU19 and its human counterpart on HSA11q13. These two physical maps span 2.5 and 0.5 megabases on MMU19. Long range restriction analysis and YAC contigs have been built, five genes were located on MMU19 and eight new STSs were generated. The 0.5-Mb map which has been positioned close to the centromere of MMU19, based on dual-color FISH experiments and genetic data, includes eight genes (Type I markers), three microsatellites (Type II markers) and five new STSs. The 2.5-Mb map is located more telomeric and contains seven genes, four microsatellites and four new STSs. Gene order and physical distances appear to be similar in human and in mouse in this 2.5-Mb region. Strikingly, the 0.5-Mb region has a similar size in human but gene order is shuffled. The overall comparative map shows that these two regions are inverted on MMU19 when compared with HSA11q13.

A 356-Kb sequence of the subtelomeric part of the MHC Class I region

The subtelomeric part of the MHC Class I region contains 11 of the 21 genes described on chromosome 6 at position 6p21.3. The general organization of those and other genes resident in the region was revealed by determining a 356,376 bp sequence. Potential exons for new genes were identified by computer analysis and a large number of ESTs were selected by testing the sequence by the BLAST algorithm against the GenBank nonredundant and EST databases. Most of the ESTs are clustered in two regions. In contrast, the whole HLA-gene region is crammed with LINE and SINE repeats, fragments of genes and microsatellites, which tends to hinder the identification of new genes.

Dissection of the 5.5 Mbp region directly telomeric of HLA-B including a long range restriction map, YAC and PAC contigs

A large number of diseases are associated with the human major histocompatibility (HLA) complex located in 6p21.3. The underlying defect of most of these has not yet been determined even after detailed analysis of the HLA region. Due to the extended haplotypes found in this area, several of the HLA-linked disease genes may be located also telomeric of the class I region. In order to analyse the area covering the 4 megabases directly telomeric of HLA-F in close detail, we have generated 50 new markers. These and other markers have been used to establish a SalI restriction map from 46 YACs. A subset of 42 markers was applied to construct a genomic long range restriction map from an HLA-A2/B13 haplotype. Both maps have been compared revealing the presence of additional 150 kb in the HLA-A2 haplotype close to the RFP locus. Additionally, 47 PACs have been selected mapping to this region and grouped into 7 contigs. Sequencing of these PAC contigs has already been initiated.

Physical mapping of two histone gene clusters on human chromosome 6p22.1-22.2

Histones are basic proteins which are responsible for the assembly and maintenance of the nucleosomal structure within the chromosomal fiber in eukaryotes. Two clusters of these genes have previously been mapped to the region 6p21.1-p22.2. We describe here a radiation hybrid map, a long range restriction map and a YAC contig covering and linking these two clusters and giving the precise localisation with respect to the HLA complex. The large cluster contains five H1 histone genes in the 6p22.2 region, the smaller only one, H1F5 (H1.5), in 6p22.1. In both clusters, each H1 locus is accompanied by several core histone genes. The large cluster has additionally been covered by a sequence ready PAC contig and three probably unrelated genes (TRMI2, BTN and SSADH) have been accurately localized within the 6p22.2-p22.1 region.

Ermap, a gene coding for a novel erythroid specific adhesion/receptor membrane protein

Ermap (erythroid membrane-associated protein), a gene coding for a novel transmembrane protein produced exclusively in erythroid cells, is described. It is mapped to murine Chromosome 4, 57 cM distal to the centromere. The initial cDNA clone was isolated from a day 9 murine embryonic erythroid cell cDNA library. The predicted peptide sequence suggests that ERMAP is a transmembrane protein with two extracellular immunoglobulin folds, as well as a highly conserved B30.2 domain and several phosphorylation consensus sequences in the cytoplasmic region. ERMAP shares a high homology throughout the entire peptide with butyrophilin, a glycoprotein essential for milk lipid droplet formation and release. A GFP-ERMAP fusion protein was localized to the plasma membrane and cytoplasmic vesicles in transiently transfected 293T cells. Northern blot analysis and in-situ hybridization demonstrated that Ermap expression was restricted to fetal and adult erythroid tissues. ERMAP is likely a novel adhesion/receptor molecule specific for erythroid cells.

Structure and evolution of the extended B7 family

Here, Joëlle Henry and colleagues explore structural and evolutionary relationships between the B7 costimulator molecules and a growing number of molecules encoded within the major histocompatibility complex. They propose that B7 and MHC genes are derived from a common ancestor, with several members of this large gene family possibly having pivotal influences on T-cell activation.

The chromosomal duplication model of the major histocompatibility complex

The major histocompatibility complex (MHC) is a genetic region that has been extensively studied by immunologists, molecular biologists, and evolutionary biologists. Nevertheless, our knowledge of how the MHC acquired its present-day organization is quite limited. The recent discovery that the mammalian genome contains regions paralogous to the MHC has led us to the proposal that the MHC region of jawed vertebrates arose as a result of ancient chromosomal duplications. Here, I review the current status of this proposal.

Comparative organization and function of the major histocompatibility complex of domesticated cattle

This review focuses on recent advances in research on the bovine major histocompatibility complex (BoLA), with specific reference to the genetic organization, polymorphism and function of the class II genes. The BoLA region is unlike the MHC of humans and mice in that a large inversion has moved several class II genes, including the TAP/LMP cluster, close to the centromere of bovine chromosome 23. Therefore, close linkage of MHC genes and other genes associated with the MHC in humans and mice does not appear to be required for normal immunological function. In cattle, polymorphism in the class IIa genes influences both the magnitude and the epitope specificity of antigen-specific T-cell responses to foot-and-mouth disease virus peptides. Disease association studies have demonstrated that BoLA alleles affect the subclinical progression of bovine leukemia virus (BLV) infection. This association is strongly correlated with the presence of specific amino acid motifs within the DRB3 antigen-binding domain. In addition to the practical significance of these findings, the association between BoLA and BLV provides a unique model to study host resistance to retrovirus infection in a non-inbred species. These studies contribute to our understanding of the evolution of the MHC in mammals, to the development of broadly effective vaccines, and to breeding strategies aimed at improving resistance to infectious diseases.

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.

Comparative gene mapping: a fine-scale survey of chromosome rearrangements between ruminants and humans

A total of 202 genes were cytogenetically mapped to goat chromosomes, multiplying by five the total number of regional gene localizations in domestic ruminants (255). This map encompasses 249 and 173 common anchor loci regularly spaced along human and murine chromosomes, respectively, which makes it possible to perform a genome-wide comparison between three mammalian orders. Twice as many rearrangements as revealed by ZOO-FISH were observed. The average size of conserved fragments could be estimated at 27 and 8 cM with humans and mice, respectively. The position of evolutionary breakpoints often correspond with human chromosome sites known to be vulnerable to rearrangement in neoplasia. Furthermore, 75 microsatellite markers, 30 of which were isolated from gene-containing bacterial artificial chromosomes (BACs), were added to the previous goat genetic map, achieving 88% genome coverage. Finally, 124 microsatellites were cytogenetically mapped, which made it possible to physically anchor and orient all the linkage groups. We believe that this comprehensive map will speed up positional cloning projects in domestic ruminants and clarify some aspects of mammalian chromosomal evolution.

Physical map of human 6p21.2-6p21.3: region flanking the centromeric end of the major histocompatibility complex

We have physically mapped and cloned a 2.5-Mb chromosomal segment flanking the centromeric end of the major histocompatibility complex (MHC). We characterized in detail 27 YACs, 144 cosmids, 51 PACs, and 5 BACs, which will facilitate the complete genomic sequencing of this region of chromosome 6. The contig contains the genes encoding CSBP, p21, HSU09564 serine kinase, ZNF76, TCP-11, RPS10, HMGI(Y), BAK, and the human homolog of Tctex-7 (HSET). The GLO1 gene was mapped further centromeric in the 6p21.2-6p21.1 region toward TCTE-1. The gene order of the GLO1-HMGI(Y) segment in respect to the centromere is similar to the gene order in the mouse t-chromosome distal inversion, indicating that there is conservation in gene content but not gene order between humans and mice in this region. The close linkage of the BAK and CSBP genes to the MHC is of interest because of their possible involvement in autoimmune disease.

The haemochromatosis candidate gene HFE (HLA-H) of man and mouse is located in syntenic regions within the histone gene cluster

The HFE (HLA-H) gene is a strong candidate gene for hereditary haemochromatosis and was localized on the short arm of chromosome 6 to 6p21.3-p22. In addition, the sequence of the homologous mouse and rat cDNA and a partial sequence from the mouse gene have been reported recently. In this report, we describe the location of the human and the mouse HFE (HLA-H) gene within the histone gene clusters on the human chromosome 6 and the mouse chromosome 13. Both the human and the murine gene were located on syntenic regions within the histone gene clusters in the vicinity of the histone H1t gene. The genomic sequence of the human HFE (HLA-H) gene and the 3' portion of the homologous mouse gene were determined. Comparison of the genomic sequences from man and mouse and the cDNA sequence from rat shows significant similarities, also beyond the transcribed region of the mouse gene.

Structure and chromosomal localization of mouse G protein subunit gamma 4 gene

The G protein gamma subunits are members of a multigene family and are implicated in determining the specificity of receptor-G protein interaction. The gene structures for many of the gamma subunits remain to be determined. Here, we report the gene structure for the brain-specific gamma 4 subunit and its map position on a mouse chromosome. The gene (Gng4) comprises at least three exons spanning over 20 kb. The 225-bp coding region, which spans two exons, is interrupted by a large 18.2-kb intron whose position is conserved in other gamma subunit genes. There is a putative additional intron in the 5' untranslated region just upstream of the translation initiation codon. Introns are present in most of the other gamma subunits at this position. The mouse Gng4 gene is mapped to chromosome 13.

BAC clones and STS markers near the distal breakpoint of the fourth t-inversion, In(17)4d, in the H2-M region on mouse chromosome 17

The H2-M region is the most distal part of the mouse major histocompatibility complex (Mhc) and is likely to include the distal breakpoint of the fourth t-inversion, In(17)4d. The conserved synteny breakpoint between mouse and human is located in the H2-M region between D17Leh89, a putative olfactory receptor gene, and Pgk2 (phosphoglycerate kinase 2). To analyze the H2-M region, we screened a mouse bacterial artificial chromosome (BAC) library, using the D17Mit64, D17Tu49, D17Leh89, D17Leh467, and Pgk2 markers. Thirty-eight BAC clones were obtained and mapped in five clusters, and 25 sequence-tagged site (STS) markers were newly developed. The regions surrounding D17Tu49 and D17Leh467 are abundant in L1 repeat sequences and may, therefore, be candidates for the breakpoints of conserved synteny and t-inversion. D17Leh89 was linked to D17Mit64 by two contiguous BAC clones. The Aeg1 (acidic epididymal glycoprotein 1) and Aeg2 genes were mapped close to Pgk2, on the same BAC clones. The genetic length between D17Leh89-D17Mit64 and Pgk2-Aeg can be estimated as 0.5-0.7 centiMorgan (cM), and the most distal class I gene, H2-M2, can be placed 0.3-1.0 cM proximal to the t-inversion breakpoint. A recombinational hotspot is suggested to be located between Aeg and Tpxl in an interspecific cross of (C57BL/6J x Mus spretus).

Genomic evolution of the distal Mhc class I region on mouse Chr 17

A 5-Mb YAC contig, partly supplemented with BAC contigs, was created from the distal Mhc class I region on mouse Chr 17. The gene order of Znf173-Tctex5-Mog-D17Tu42-D17Leh 89 is conserved between mouse and human but not the physical distance, supporting the independent expansion of Mhc class I genes in the so-called accordion model of Mhc evolution. The distal H2-M region includes the breakpoint of conserved synteny between mouse and human as well as the In(17)4 t-inversion. The H2-M region is rich in L1 repeats, implying that the insertion of L1 repeats may be associated with the evolutionary flexibility to break a chromosome.

Zooming in on the human-mouse comparative map: genome conservation re-examined on a high-resolution scale

Over the past decade, conservation of genetic linkage groups has been shown in mammals and used to great advantage, fueling significant exchanges of gene mapping and functional information especially between the genomes of humans and mice. As human physical maps increase in resolution from chromosome bands to nucleotide sequence, comparative alignments of mouse and human regions have revealed striking similarities and surprising differences between the genomes of these two best-mapped mammalian species. Whereas, at present, very few mouse and human regions have been compared on the physical level, existing studies provide intriguing insights to genome evolution, including the observation of recent duplications and deletions of genes that may play significant roles in defining some of the biological differences between the two species. Although high-resolution conserved marker-based maps are currently available only for human and mouse, a variety of new methods and resources are speeding the development of comparative maps of additional organisms. These advances mark the first step toward establishment of the human genome as a reference map for vertebrate species, providing evolutionary and functional annotation to human sequence and vast new resources for genetic analysis of a variety of commercially, medically, and ecologically important animal models.

Identification of a Val 145 Ile substitution in the human myelin oligodendrocyte glycoprotein: lack of association with multiple sclerosis. The Réseau de Recherche Clinique INSERM sur la Susceptibilité Génétique à la Sclérose en Plaques

Myelin/oligodendrocyte glycoprotein (MOG) is a major target antigen in experimental autoimmune encephalomyelitis and it has been suggested that it may as well play a key role in the demyelination process in multiple sclerosis (MS). As MOG variants could be pathogenic in autoimmune demyelinating diseases of the central nervous system, we analysed the coding sequence of MOG in MS patients and described a G-->A transition occurring in exon 3 of the human MOG gene. The mutation predicts that isoleucine substitutes for a valine at codon 145 (Val 145 Ile) in the transmembrane region of the protein. This is the first aminoacid substitution reported in human MOG. The polymorphism can be detected by restriction enzyme digestion of genomic DNA or reverse-transcribed PCR amplified products, making it a simple tool to detect a potential implication of MOG alleles in susceptibility to MS by association study. The analysis of 83 unrelated MS patients and 82 unrelated healthy controls showed that the polymorphism is found in similar proportions in MS patients (18%) and controls (14.6%). It is therefore unlikely that the MOG Val 145 Ile variant is responsible for genetic susceptibility to MS.

A 1.1-Mb transcript map of the hereditary hemochromatosis locus

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.

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.

Identification of a mouse homolog for the human hereditary haemochromatosis candidate gene

Recently, a novel human major histocompatibility complex (MHC) class I-like gene (HLA-H) was reported as a candidate gene for human hereditary haemochromatosis, a recessive disease of iron metabolism with a remarkably high incidence in northern Europeans. Independently we have isolated this gene in the course of a search for new human MHC class I-related genes and named it MR2. Here we report a mouse homolog of this human gene. The mouse MR2 gene is similar to the human counterpart with an overall predicted amino acid sequence similarity of approximately 66% and it is expressed in various tissues as in human. The extra eight amino acid residues between the alpha1 and the alpha2 domains in the mouse molecule compared to the human counterpart can be explained by the creation of the coding sequence from the intron. While the human gene is located at the site telomeric to the MHC region on human chromosome 6, our study indicated the translocation of the mouse homolog from the site telomeric to the MHC on mouse chromosome 17 to chromosome 13 along with other genes. This mouse gene should be important in clarifying a possible role in iron metabolism.

H2-M3, a full-service class Ib histocompatibility antigen

H2-M3 is an MHC class Ib molecule of the mouse with a unique preference for N-formylated peptides, which may come from the N-termini of endogenous, mitochondrial proteins or foreign, bacterial proteins. The crystal structure of M3 revealed a hydrophobic peptide-binding groove with an occluded A pocket and the peptide shifted one residue relative to class Ia structures. The formyl group is held by a novel hydrogen bonding network, involving His9 on the bottom of the groove, and the side chain of the P1 methionine is lodged in the B pocket. M3 is a full-service histocompatibility (H) antigen, i.e. self-M3 can present endogenous peptides as minor H antigens and foreign, bacterial antigens in a defensive immune response to infection; and foreign M3 complexed with endogenous self-peptides.

Structural organization and mammary-specific expression of the butyrophilin gene

Butyrophilin is a glycoprotein of the immunoglobulin superfamily that is secreted in association with the milk-fat-globule membrane from mammary epithelial cells. As a first step towards determining the possible function(s) of this protein in lactation, the mouse butyrophilin gene (Btn) has been cloned from a 129-ES cell genomic library. Over 14 kb of DNA was sequenced, including the entire transcriptional unit of the gene, and 4.6 kb and 1.1 kb of the 5' and 3' flanking region, respectively. In addition, the overall structure of the bovine gene (BTN) was determined by amplification of genomic DNA by the polymerase chain reaction. Both Btn and BTN comprise seven exons and six introns. The signal sequence and two immunoglobulin-like folds of the exoplasmic domain and the membrane anchor are encoded by separate exons, and the cytoplasmic domain is encoded by two short exons and a large terminal exon that also includes 3' untranslated sequence. The butyrophilin gene appears to have evolved from a subset of genes in the immunoglobulin superfamily and genes encoding the B30.2 domain, which is conserved in a family of zinc-finger proteins. Murine butyrophilin mRNA was detected specifically in the mammary gland by RNase protection analysis. Expression increased during the last half of pregnancy and was maximal during lactation. The 5' flanking region of Btn was analyzed for putative regulatory elements and is different from the promoters of other mammary-specific genes. Btn should be useful for determining the mechanisms underlying mammary-specific gene expression and potentially for the production of heterologous proteins in the milk of transgenic animals.

Comparative maps: adding pieces to the mammalian jigsaw puzzle

Comparative maps display the chromosomal location of homologous genes in different species and highlight genetic segments that are conserved in evolution. These maps are used to study chromosomal changes that occurred during the divergence of mammalian lineages, to identify candidates for hereditary disease genes, and to facilitate mapping in other species. Recently, physical mapping in regions of known conserved linkage has revealed previously undetected chromosomal changes that may provide clues to understanding chromosomal structure and function and evolutionary processes. The availability of these data in electronically accessible formats is critical to the growth and analysis of comparative maps.

Cloning of a human homologue of the mouse Tctex-5 gene within the MHC class I region

Using a positional cloning strategy to identify the hemochromatosis gene (HFE), we isolated seven cDNAs by cDNA selection from a region of 400 kilobases (kb) located near the HLA-A and HLA-F loci. In this paper, we report the study of one of the corresponding genes, referred to as HCG V (hemochromatosis candidate gene), localized 150 kb centromeric to HLA-A. This gene was found to be expressed ubiquitously in the form of a 1.8 kb transcript, and to be apparently well conserved during evolution. The gene spanned 3.1 kb and is organized in three exons and two introns. The cDNA of 1620 base pairs (bp) showed an open reading frame of 378 bp, encoding for a 126 amino acid polypeptide which displayed a strong identity with the predicted product of a mouse Tctex-5 gene (t complex, testis expressed) localized in the t complex on chromosome 17. The HCG V gene was assessed as a potential candidate for hemochromatosis in regard to its localization in the linkage disequilibrium area between HFE and polymorphic markers. The study of deletions and point mutations in hemochromatosis patients revealed a single bp polymorphism within the coding region; however, no associated disease changes were found. Therefore we conclude that HCG V is unlikely to be involved in the pathogenesis of hemochromatosis.

Structure and content of the major histocompatibility complex (MHC) class I regions of the great anthropoid apes

The origins of the functional class I genes predated human speciation, a phenomenon known as trans-speciation. The retention of class Ia orthologues within the great apes, however, has not been paralleled by studies designed to examine the pseudogene content, organization, and structure of their class I regions. Therefore, we have begun the systematic characterization of the Old World primate MHCs. The numbers and sizes of fragments harboring class I sequences were similar among the chimpanzee, gorilla, and human genomes tested. Both of the gorillas included in our study possessed genomic fragments carrying orthologues of the recently evolved HLA-H pseudogene identical to those found in the human. The overall megabase restriction fragment patterns of humans and chimpanzees appeared slightly more similar to each other, although the HLA-A subregional megabase variants may have been generated following the emergence of Homo sapiens. Based on the results of this initial study, it is difficult to generate a firm species tree and to determine human's closest evolutionary neighbor. Nevertheless, an analysis of MHC subregional similarities and differences in the hominoid apes may ultimately aid in localizing and identifying MHC haplotype-associated disease genes such as idiopathic hemochromatosis.

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.

The bovine butyrophilin gene maps to chromosome 23

Chromosomal assignment of the bovine butyrophilin gene (BTN) was performed by analysis of DNA from somatic hybrid cell lines using the polymerase chain reaction. The gene was assigned to bovine chromosome 23 using two sets of primers specific for bovine BTN.

Restriction fragment length polymorphism in amplification products of the bovine butyrophilin gene: assignment of bovine butyrophilin to bovine chromosome 23

A polymorphism was identified in the bovine butyrophilin (BTN) gene by digesting polymerase chain reaction products with the restriction enzyme HaeIII. This polymorphism was segregating in a Holstein-Friesian sire selected as part of an ongoing study directed towards the identification of quantitative trait loci affecting milk composition. Screening of a half-sib family established for the heterozygous sire allowed the localization of BTN to bovine chromosome 23 (BTA23).

Identification of seven new human MHC class I region genes around the HLA-F locus

Using cDNA hybridization selection techniques, we identified seven new genes in a 280 kilobase YAC covering the HLA-F locus. The new genes were mapped back to the YAC by a combination of optical restriction mapping and pulse field gel electrophoresis. Northern analysis of individual clones demonstrated the presence of either different mRNA sizes or different expression patterns. Two of the cDNA clones were expressed only in lymphoid cell lines: one in Jurkat cells (T cell) and another in JY cells (B cell). All the genes lacked sequence similarity to any known classical and non-classical major histocompatibility complex (MHC) class I genes, indicating that the MHC class I region has more functions than anticipated. Of the seven new genes, one is highly similar (97%) to mouse 60S ribosomal protein, and another is homologous to diubiquitin proteins. Of the two G-coupled receptor-like cDNAs, one was fully sequenced and found to be an olfactory receptor-like gene. The study strengthens evidence that the MHC complex not only plays a key role in the immune system, but also contributes to non-immunological functions.

Polymorphic tri- and tetranucleotide repeats in exons 1 and 8 of the myelin oligodendrocyte glycoprotein (MOG) gene

Polymorphic (CTC)n and (TAAA)n sequences were identified in exons 1 and 8 of the myelin oligodendrocyte glycoprotein (MOG) gene. The different alleles were detected by a method combining fluorescence labeling of polymerase chain reaction (PCR) products and use of an automated DNA sequencer. Although only two alleles differing by the number of leucine residues encoded by the (CTC)n array were detected at the first locus, seven alleles were identified at the second. The high degree of polymorphism (75%) of the tetranucleotide repeat makes this marker informative for association or linkage studies with diseases such as hemochromatosis or multiple sclerosis.

Three highly polymorphic microsatellites at the human myelin oligodendrocyte glycoprotein locus, 100 kb telomeric to HLA-F. Characterization and relation to HLA haplotypes

The MOG locus, located on chromosomal bands 6p21.3-p22 and mapped about 100 kb telomeric to HLA-F, was isolated from cosmid ICRFc109A2434 and shown to contain three microsatellites. These CA-repeat polymorphic markers were characterized in a sample of 173 healthy unrelated individuals and 84 DNAs from the HLA Workshop reference panel, by a method combining fluorescence labeling of PCR products and use of an automated DNA sequencer. For the three markers, frequencies of heterozygotes are well predicted from allele frequencies by the Hardy-Weinberg rule, which suggests that problems of allele nonamplification are unlikely. Typing of cell lines homozygous in the HLA region allowed unambiguous definition of 81 HLA-MOG haplotypes and showed that several HLA ancestral haplotypes extended to the MOG region. The high degree of polymorphism (59%, 51%, and 81% at the three loci, respectively, and 87% at the haplotype level) makes these new markers informative for association or linkage studies with diseases such as hemochromatosis or multiple sclerosis, and for studies aimed at precisely delineating the site of crossover in chromosomes in which recombination occurred in the distal part of the HLA class I region.