Analysis of the gene-dense major histocompatibility complex class III region and its comparison to mouse

High-density SNP screening of the major histocompatibility complex in systemic lupus erythematosus demonstrates strong evidence for independent susceptibility regions

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10⁻¹⁶). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10⁻¹²), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10⁻¹³. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and involvement of multiple organ systems. Although the cause of SLE remains unknown, several lines of evidence underscore the importance of genetic factors. As is true for most autoimmune diseases, a substantial genetic contribution to disease risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6. This region of the genome contains a large number of genes that participate in the immune response. However, the full contribution of this genomic region to SLE risk has not yet been defined. In the current study we characterize a large number of SLE patients and family members for approximately 2,000 MHC region variants to identify the specific genes that influence disease risk. Our results, for the first time, implicate four different MHC regions in SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

An extensive analysis of the hereditary hemochromatosis gene HFE and neighboring histone genes: associations with childhood leukemia

The most common mutation of the HFE gene C282Y has shown a risk association with childhood acute lymphoblastic leukemia (ALL) in Welsh and Scottish case-control studies. This finding has not been replicated outside Britain. Here, we present a thorough analysis of the HFE gene in a panel of HLA homozygous reference cell lines and in the original population sample from South Wales (117 childhood ALL cases and 414 newborn controls). The 21 of 24 variants analyzed were from the HFE gene region extending 52 kb from the histone gene HIST1H1C to HIST1H1T. We identified the single-nucleotide polymorphism (SNP) rs807212 as a tagging SNP for the most common HFE region haplotype, which contains wild-type alleles of all HFE variants examined. This intergenic SNP rs807212 yielded a strong male-specific protective association (per allele OR = 0.38, 95% CI = 0.22-0.64, P (trend) = 0.0002; P = 0.48 in females), which accounted for the original C282Y risk association. In the HapMap project data, rs807212 was in strong linkage disequilibrium with 25 other SNPs spanning 151 kb around HFE. Minor alleles of these 26 SNPs characterized the most common haplotype for the HFE region, which lacked all disease-associated HFE variants. The HapMap data suggested positive selection in this region even in populations where the HFE C282Y mutation is absent. These results have implications for the sex-specific associations observed in this region and suggest the inclusion of rs807212 in future studies of the HFE gene and the extended HLA class I region.

The human Major Histocompatibility Complex as a paradigm in genomics research

Since its discovery more than 50 years ago, the human Major Histocompatibility Complex (MHC) on chromosome 6p21.3 has been at the forefront of human genetic research. Here, we review from a historical perspective the major advances in our understanding of the nature and consequences of genetic variation which have involved the MHC, as well as highlighting likely future directions. As a consequence of its particular genomic structure, its remarkable polymorphism and its early implication in numerous diseases, the MHC has been considered as a model region for genomics, being the first substantial region to be sequenced and establishing fundamental concepts of linkage disequilibrium, haplotypic structure and meiotic recombination. Recently, the MHC became the first genomic region to be entirely re-sequenced for common haplotypes, while studies mapping gene expression phenotypes across the genome have strongly implicated variation in the MHC. This review shows how the MHC continues to provide new insights and remains in the vanguard of contemporary research in human genomics.

Identification and characterization of human LYPD6, a new member of the Ly-6 superfamily

The Ly-6 protein superfamily is usually identified as a group of proteins with a LU protein domain. LU domain is about 80 amino acids long and characterized by a conserved pattern of 10 cysteine residues. Here we report the cloning and characterization of a novel human LU domain containing gene, LYPD6, isolated from human testis cDNA library, and mapped to 2q23.1-23.2 by searching the UCSC genomic database. The LYPD6 cDNA sequence of 3,501 base pairs contains an open reading frame encoding 171 amino acids. Subcellular localization of LYPD6 demonstrated that the protein was localized in the cytoplasm when overexpressed in COS-7 cells. RT-PCR analysis showed that LYPD6 was widely expressed in human tissues and the expression levels in brain and heart were relatively high. Furthermore, the subsequent analysis based on reporter gene assays suggested that overexpression of LYPD6 in HEK 293T cells was able to suppress the transcriptional activities of AP1.

Regulation of expression of two LY-6 family genes by intron retention and transcription induced chimerism

BACKGROUND

Regulation of the expression of particular genes can rely on mechanisms that are different from classical transcriptional and translational control. The LY6G5B and LY6G6D genes encode LY-6 domain proteins, whose expression seems to be regulated in an original fashion, consisting of an intron retention event which generates, through an early premature stop codon, a non-coding transcript, preventing expression in most cell lines and tissues.

RESULTS

The MHC LY-6 non-coding transcripts have shown to be stable and very abundant in the cell, and not subject to Nonsense Mediated Decay (NMD). This retention event appears not to be solely dependent on intron features, because in the case of LY6G5B, when the intron is inserted in the artificial context of a luciferase expression plasmid, it is fully spliced but strongly stabilises the resulting luciferase transcript. In addition, by quantitative PCR we found that the retained and spliced forms are differentially expressed in tissues indicating an active regulation of the non-coding transcript. EST database analysis revealed that these genes have an alternative expression pathway with the formation of Transcription Induced Chimeras (TIC). This data was confirmed by RT-PCR, revealing the presence of different transcripts that would encode the chimeric proteins CSNKbeta-LY6G5B and G6F-LY6G6D, in which the LY-6 domain would join to a kinase domain and an Ig-like domain, respectively.

CONCLUSION

In conclusion, the LY6G5B and LY6G6D intron-retained transcripts are not subjected to NMD and are more abundant than the properly spliced forms. In addition, these genes form chimeric transcripts with their neighbouring same orientation 5' genes. Of interest is the fact that the 5' genes (CSNKbeta or G6F) undergo differential splicing only in the context of the chimera (CSNKbeta-LY6G5B or G6F-LY6G6C) and not on their own.

Chromatin profiling across the human tumour necrosis factor gene locus reveals a complex, cell type-specific landscape with novel regulatory elements

The TNF locus on chromosome 6p21 encodes a family of proteins with key roles in the immune response whose dysregulation leads to severe disease. Transcriptional regulation is important, with cell type and stimulus-specific enhancer complexes involving the proximal TNF promoter. We show how quantitative chromatin profiling across a 34 kb region spanning the TNF locus has allowed us to identify a number of novel DNase hypersensitive sites and characterize more distant regulatory elements. We demonstrate DNase hypersensitive sites corresponding to the lymphotoxin alpha (LTA) and tumour necrosis factor (TNF) promoter regions, a CpG island in exon 4 of lymphotoxin beta (LTB), the 3′ end of nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor-like 1 (NFKBIL1) and 3.4 kb upstream of LTA. These sites co-localize to highly conserved DNA sequences and show evidence of cell type specificity when lymphoblastoid, Jurkat, U937, HeLa and HEK293T cell lines are analysed using Southern blotting. For Jurkat T cells, we define histone modifications across the locus. Peaks of acetylated histone H3 and H4, together with tri-methyl K4 of histone H3, correspond to hypersensitive sites, notably in exon 4 of LTB. We provide evidence of a functional role for an intergenic DNase I hypersensitive site distal to LTA in Jurkat cells based on reporter gene analysis, with evidence of recruitment of upstream stimulatory factors (USF) transcription factors.

Identification and hormonal regulation of a novel form of NKp30 in human endometrial epithelium

This study reports the discovery and steroid hormonal regulation of a novel glycosylated form of NKp30 in human endometrial epithelium. NKp30 is a member of the immunoglobulin superfamily and one of three existing natural cytotoxicity-triggering receptors. NKp30 is a glycosylated protein and is thought to be selectively expressed in resting and activated natural killer cells. The aims of the present study were to fully characterize NKp30 mRNA and protein in human endometrium during the menstrual cycle, and to investigate the hormonal regulation of NKp30. NKp30 mRNA was significantly up-regulated in fresh tissues during late secretory phase of the menstrual cycle. Interestingly, NKp30 mRNA was also present in clonally derived endometrial epithelial cells (CEE) in comparable amounts to fresh tissue. NKp30 protein was predominantly found in the endometrial glands and luminal epithelia of the secretory phase endometrium. Western blotting and de-glycosylation studies show that a novel glycosylated form of NKp30 is present in endometrial epithelium and that it can dimerize. Further phenotyping of CEE by flow cytometry revealed that they are CK8(+)CD49f(+)NKp30(+)CD45(-)CD56(-). The data also show that transcription and translation of the novel form of NKp30 can be induced by progesterone treatment after 48 h in endometrial explants in vitro. This is the first report to show the presence of both NKp30 mRNA and a novel glycosylated form of NKp30 protein in endometrial epithelial cells.

Activation-dependent intrachromosomal interactions formed by the TNF gene promoter and two distal enhancers

Here we provide a mechanism for specific, efficient transcription of the TNF gene and, potentially, other genes residing within multigene loci. We identify and characterize highly conserved noncoding elements flanking the TNF gene, which undergo activation-dependent intrachromosomal interactions. These elements, hypersensitive site (HSS)-9 and HSS+3 (9 kb upstream and 3 kb downstream of the TNF gene, respectively), contain DNase I hypersensitive sites in naive, T helper 1, and T helper 2 primary T cells. Both HSS-9 and HSS+3 inducibly associate with acetylated histones, indicative of chromatin remodeling, bind the transcription factor nuclear factor of activated T cells (NFAT)p in vitro and in vivo, and function as enhancers of NFAT-dependent transactivation mediated by the TNF promoter. Using the chromosome conformation capture assay, we demonstrate that upon T cell activation intrachromosomal looping occurs in the TNF locus. HSS-9 and HSS+3 each associate with the TNF promoter and with each other, circularizing the TNF gene and bringing NFAT-containing nucleoprotein complexes into close proximity. TNF gene regulation thus reveals a mode of intrachromosomal interaction that combines a looped gene topology with interactions between enhancers and a gene promoter.

Evolution and comparative analysis of the MHC Class III inflammatory region

BACKGROUND

The Major Histocompatibility Complex (MHC) is essential for immune function. Historically, it has been subdivided into three regions (Class I, II, and III), but a cluster of functionally related genes within the Class III region has also been referred to as the Class IV region or "inflammatory region". This group of genes is involved in the inflammatory response, and includes members of the tumour necrosis family. Here we report the sequencing, annotation and comparative analysis of a tammar wallaby BAC containing the inflammatory region. We also discuss the extent of sequence conservation across the entire region and identify elements conserved in evolution.

RESULTS

Fourteen Class III genes from the tammar wallaby inflammatory region were characterised and compared to their orthologues in other vertebrates. The organisation and sequence of genes in the inflammatory region of both the wallaby and South American opossum are highly conserved compared to known genes from eutherian ("placental") mammals. Some minor differences separate the two marsupial species. Eight genes within the inflammatory region have remained tightly clustered for at least 360 million years, predating the divergence of the amphibian lineage. Analysis of sequence conservation identified 354 elements that are conserved. These range in size from 7 to 431 bases and cover 15.6% of the inflammatory region, representing approximately a 4-fold increase compared to the average for vertebrate genomes. About 5.5% of this conserved sequence is marsupial-specific, including three cases of marsupial-specific repeats. Highly Conserved Elements were also characterised.

CONCLUSION

Using comparative analysis, we show that a cluster of MHC genes involved in inflammation, including TNF, LTA (or its putative teleost homolog TNF-N), APOM, and BAT3 have remained together for over 450 million years, predating the divergence of mammals from fish. The observed enrichment in conserved sequences within the inflammatory region suggests conservation at the transcriptional regulatory level, in addition to the functional level.

Characterization of the five novel Ly-6 superfamily members encoded in the MHC, and detection of cells expressing their potential ligands

Lymphocyte Antigen 6 (Ly-6) superfamily members are cysteine-rich, generally GPI-anchored cell surface proteins, which have definite or putative immune related roles. There are 27 members of this family described so far in the human genome and 37 in the mouse. Five of them are clustered in the class III region of the human and mouse MHCs. Following computational analyses, we functionally characterized the encoded proteins by creating epitope-tagged fusion constructs to determine molecular weight, complex formation, subcellular localization, post-translational modifications and ligand binding. We found that all human and mouse proteins were glycosylated, and most could form part of larger complexes. Human and mouse Ly6G6c and Ly6G6d, and mouse Ly6g6e were found to be GPI-anchored cell surface proteins, highly expressed at the leading edges of cells, on filopodia, which are normally involved in cell adhesion and migration. However, analysis of Ly6G5c and Ly6G5b indicated that they are potentially secreted proteins. Our results indicate that there are two subclusters of related Ly-6 proteins in this region of the MHC, with Ly6G6c, Ly6G6d, and Ly6G6e forming one and Ly6G5c and Ly6G5b forming another. In addition, by FACS analysis we have found that the potential ligands for human LY6G6C, LY6G6D, and LY6G5C are expressed on K562 cells, an undifferentiated megakaryocyte cell line, indicating a potential role in hematopoietic cell differentiation. This characterization of the five MHC class III region Ly-6 family members is of great relevance, as they represent 18% of the human Ly-6 protein family and 50% of the secreted ones.

YAC/BAC contig spanning the MHC class III region of cattle

A contig of the class III region of the bovine major histocompatibility complex (MHC) was established from bacterial and yeast artificial chromosomes using PCR and BAC-end sequencing. The marker content of individual clones was determined by gene and BAC-end specific PCR, and the location of genes and BAC-ends was confirmed analyzing somatic hybrid cells. A comparative analysis indicated that the content and order of MHC class III genes is strongly conserved between cattle and other mammalian species. Fluorescence in situ hybridization localized the bovine class III region to BTA23q21-->q22. The results show that the collection of sequenced BAC-ends is a powerful resource for generating high-resolution comparative chromosome maps.

The genomic sequence and analysis of the swine major histocompatibility complex

We describe the generation and analysis of an integrated sequence map of a 2.4-Mb region of pig chromosome 7, comprising the classical class I region, the extended and classical class II regions, and the class III region of the major histocompatibility complex (MHC), also known as swine leukocyte antigen (SLA) complex. We have identified and manually annotated 151 loci, of which 121 are known genes (predicted to be functional), 18 are pseudogenes, 8 are novel CDS loci, 3 are novel transcripts, and 1 is a putative gene. Nearly all of these loci have homologues in other mammalian genomes but orthologues could be identified with confidence for only 123 genes. The 28 genes (including all the SLA class I genes) for which unambiguous orthology to genes within the human reference MHC could not be established are of particular interest with respect to porcine-specific MHC function and evolution. We have compared the porcine MHC to other mammalian MHC regions and identified the differences between them. In comparison to the human MHC, the main differences include the absence of HLA-A and other class I-like loci, the absence of HLA-DP-like loci, and the separation of the extended and classical class II regions from the rest of the MHC by insertion of the centromere. We show that the centromere insertion has occurred within a cluster of BTNL genes located at the boundary of the class II and III regions, which might have resulted in the loss of an orthologue to human C6orf10 from this region.

An investigation of transmission ratio distortion in the central region of the human MHC

Transmission ratio distortion (TRD) describes a significant departure from expected Mendelian inheritance ratios that is fundamental to both the biology of reproduction and statistical genetics. The relatively high fetal wastage in humans, with consequent selection of alleles in utero, makes it likely that TRD is prevalent in the human genome. The central region of the human major histocompatibility complex (MHC) is a strong TRD candidate, as it houses a number of immune and regulatory genes that may be important in pregnancy outcome. We used a nonhaplotype-based method to select 13 tagging SNPs from three central MHC candidate regions, and analysed their transmission in 380 newborns and their parents (1138 individuals). A TRD of 54:46 was noted in favour of the common allele of a promoter SNP in the CLIC1 gene (P = 0.025), with a similar distortion using haplotypes across the same gene region (P = 0.016). We also found evidence that markers in the CLIC1 gene region may have been subject to recent selection (P < 0.001). The study illustrates the potential benefits of screening for TRD and highlights the difficulties encountered therein.

Epistatic suppression of systemic lupus erythematosus: fine mapping of Sles1 to less than 1 mb

Sle is a susceptibility locus for systemic autoimmunity derived from the lupus-prone NZM2410 mouse. The New Zealand White-derived suppressive modifier Sles1 was identified as a specific modifier of Sle1 and prevents the development of IgG anti-chromatin autoantibodies mediated by Sle1 on the C57BL/6 (B6) background. Fine mapping of Sles1 with truncated congenic intervals localizes it to a approximately 956-kb segment of mouse chromosome 17. Sles1 completely abrogates the development of activated T and B cell populations in B6.Sle1. Despite this suppression of the Sle1-mediated cell surface activation phenotypes, B6.Sle1 Sles1 splenic B cells still exhibit intrinsic ERK phosphorylation. Classic genetic complementation tests using the nonautoimmmune 129/SvJ mouse suggests that this strain possesses a Sles1 allele complementary to that of New Zealand White, as evidenced by the lack of glomerulonephritis, splenomegaly, and antinuclear autoantibody production seen in (129 x B6.Sle1 Sles1)F(1)s. These findings localize and characterize the suppressive properties of Sles1 and implicate 129 as a useful strain for aiding in the identification of this elusive epistatic modifier gene.

Features of MHC and NK gene clusters

When comparing the immune genome to the genome in general, a higher prevalence for association with disease is the only genetic feature significant in immune genes as a group. However, some genetic features, such as marked levels of polymorphism and gene duplication, are present in subsets of immune genes, namely the Major Histocompatibility Complex (MHC) and Natural Killer (NK) cell receptor gene complexes. In this review, we discuss features of MHC and NK receptor gene clusters, their epistatic interactions, and the impact of both on association to disease.

Genetics of the central MHC

The MHC, primarily known for its antigen-presenting class I and II molecules, harbours, within a central segment of less than 1 Mb, a dense collection of genes involved in various biological functions. Although MHC I and MHC II are principal players of adaptive immunity, several loci within this central (still called class III) MHC region encode members of the innate immune system. These include the long known factors of the complement system--potentially inhibitory and triggering natural killer receptors as well as stress proteins. Whether this physical proximity is fortuitous or functionally advantageous is an important question for the future of MHC genetics.

Comparative genomics of major histocompatibility complexes

The major histocompatibility complex (MHC) is a gene dense region found in all jawed vertebrates examined to date. The MHC contains a high percentage of immune genes, in particular genes involved in antigen presentation, which are generally highly polymorphic. The region plays an important role in disease resistance. The clustering of MHC genes could be advantageous for co-evolution or regulation, and its study in many species is desirable. Even though some linkage of MHC genes is apparent in all gnathostomes, the genomic organization can differ greatly by species, suggesting rapid evolution of MHC genes after divergence from a common ancestor. Previous reviews of comparative MHC organization have been written when relatively fragmentary sequence and mapping data were available on many species. This review compares maps of MHC gene orders in commonly studied species, where extensive sequencing has been performed.

Gene map of the extended human MHC

The major histocompatibility complex (MHC) is the most important region in the vertebrate genome with respect to infection and autoimmunity, and is crucial in adaptive and innate immunity. Decades of biomedical research have revealed many MHC genes that are duplicated, polymorphic and associated with more diseases than any other region of the human genome. The recent completion of several large-scale studies offers the opportunity to assimilate the latest data into an integrated gene map of the extended human MHC. Here, we present this map and review its content in relation to paralogy, polymorphism, immune function and disease.

Complete complement components C4A and C4B deficiencies in human kidney diseases and systemic lupus erythematosus

Although a heterozygous deficiency of either complement component C4A or C4B is common, and each has a frequency of approximately 20% in a Caucasian population, complete deficiencies of both C4A and C4B proteins are extremely rare. In this paper the clinical courses for seven complete C4 deficiency patients are described in detail, and the molecular defects for complete C4 deficiencies are elucidated. Three patients with homozygous HLA A24 Cw7 B38 DR13 had systemic lupus erythematosus, mesangial glomerulonephritis, and severe skin lesions or membranous nephropathy. Immunofixation, genomic restriction fragment length polymorphisms, and pulsed field gel electrophoresis experiments revealed the presence of monomodular RP-C4-CYP21-TNX (RCCX) modules, each containing a solitary, long C4A mutant gene. Sequencing of the mutant C4A genes revealed a 2-bp, GT deletion in exon 13 that leads to protein truncation. The other four patients with homozygous HLA A30 B18 DR7 had SLE, severe kidney disorders including mesangial or membranoproliferative glomerulonephritis, and/or Henoch Schoenlein purpura. Molecular genetic analyses revealed an unusual RCCX structure with two short C4B mutant genes, each followed by an intact gene for steroid 21-hydroxylase. Nine identical, intronic mutations were found in each mutant C4B. In particular, the 8127 g-->a mutation present at the donor site of intron 28 may cause an RNA splice defect. Analyses of 12 complete C4 deficiency patients revealed two hot spots of deleterious mutations: one is located at exon 13, the others within a 2.6-kb genomic region spanning exons 20-29. Screening of these mutations may facilitate epidemiologic studies of C4 in infectious, autoimmune, and kidney diseases.

Complete MHC haplotype sequencing for common disease gene mapping

The future systematic mapping of variants that confer susceptibility to common diseases requires the construction of a fully informative polymorphism map. Ideally, every base pair of the genome would be sequenced in many individuals. Here, we report 4.75 Mb of contiguous sequence for each of two common haplotypes of the major histocompatibility complex (MHC), to which susceptibility to >100 diseases has been mapped. The autoimmune disease-associated-haplotypes HLA-A3-B7-Cw7-DR15 and HLA-A1-B8-Cw7-DR3 were sequenced in their entirety through a bacterial artificial chromosome (BAC) cloning strategy using the consanguineous cell lines PGF and COX, respectively. The two sequences were annotated to encompass all described splice variants of expressed genes. We defined the complete variation content of the two haplotypes, revealing >18,000 variations between them. Average SNP densities ranged from less than one SNP per kilobase to >60. Acquisition of complete and accurate sequence data over polymorphic regions such as the MHC from large-insert cloned DNA provides a definitive resource for the construction of informative genetic maps, and avoids the limitation of chromosome regions that are refractory to PCR amplification.