Nucleotide sequence analysis of the HLA class I region spanning the 237-kb segment around the HLA-B and -C genes

Fibrosis and Hepatocarcinogenesis: Role of Gene-Environment Interactions in Liver Disease Progression

The liver is a complex organ that performs vital functions in the body. Despite its extraordinary regenerative capacity compared to other organs, exposure to chemical, infectious, metabolic and immunologic insults and toxins renders the liver vulnerable to inflammation, degeneration and fibrosis. Abnormal wound healing response mediated by aberrant signaling pathways causes chronic activation of hepatic stellate cells (HSCs) and excessive accumulation of extracellular matrix (ECM), leading to hepatic fibrosis and cirrhosis. Fibrosis plays a key role in liver carcinogenesis. Once thought to be irreversible, recent clinical studies show that hepatic fibrosis can be reversed, even in the advanced stage. Experimental evidence shows that removal of the insult or injury can inactivate HSCs and reduce the inflammatory response, eventually leading to activation of fibrolysis and degradation of ECM. Thus, it is critical to understand the role of gene-environment interactions in the context of liver fibrosis progression and regression in order to identify specific therapeutic targets for optimized treatment to induce fibrosis regression, prevent HCC development and, ultimately, improve the clinical outcome.

First characterization of major histocompatibility complex class II DRB3 diversity in cattle breeds raised in Egypt

Genes encoding bovine leukocyte antigen (BoLA) enable the immune system to identify pathogens. Therefore, these genes have been used as genetic markers for infectious and autoimmune diseases as well as for immunological traits in cattle. Although BoLA polymorphisms have been reported in various cattle breeds worldwide, they have not been studied in cattle populations in Egypt. In this study, we characterized BoLA-DRB3 in two local Egyptian populations and one foreign population using polymerase chain reaction-sequence-based typing (PCR-SBT) method. Fifty-four previously reported BoLA-DRB3 alleles and eight new alleles (BoLA-DRB3*005:08, *015:07, *016:03, *017:04, *020:02:02, *021:03, *164:01, and *165:01) were identified. Alignment analysis of the eight new alleles revealed 90.7-98.9 %, and 83.1-97.8 % nucleotide and amino acid identities, respectively, with the BoLA-DRB3 cDNA clone NR-1. Interestingly, BoLA-DRB3 in Egyptian cattle showed a high degree of allelic diversity in native (na = 28, hE > 0.95), mixed (na = 61, hE > 0.96), and Holstein (na = 18, hE > 0.88) populations. BoLA-DRB3*002:01 (14.3 %), BoLA-DRB3*001:01 (8.5 %), and BoLA-DRB3*015:01 (20.2 %) were the most frequent alleles in native, mixed, and Holstein populations, respectively, indicating that the genetic profiles differed in each population. Based on the allele frequencies of BoLA-DRB3, genetic variation among Egyptian, Asian, African, and American breeds was examined using Nei's distances and principal component analysis. The results suggested that native and mixed cattle populations were most closely associated with African breeds in terms of their gene pool, whereas Holstein cattle were more distinct from the other breeds and were closely related to Holstein cattle populations from other countries.

Potential Molecular Markers Related to Lymph Node Metastasis and Stalk Resection Margins in Pedunculated T1 Colorectal Cancers Using Digital Spatial Profiling: A Pilot Study with a Small Case Series

There is a debate regarding the prediction of lymph node metastasis (LNM) in pedunculated T1 colorectal cancer (CRC). In this study with four cases of pedunculated T1 CRCs, we aimed to investigate gene expression variations based on the distance from the Haggitt line (HL) and identify potential molecular risk factors for LNM. By leveraging the Cancer Transcriptome Atlas and digital spatial profiling technology, we meticulously analyzed discrete regions, including the head, HL, proximal stalk region (300-1000 μm from HL), and distal stalk region (1500-2000 μm from HL) to identify spatially sequential molecular changes. Our findings showed significant overall gene expression variations among the head, proximal stalk, and distal stalk regions of pedunculated T1 CRCs compared to the control adenoma. Compared to LNM-negative T1 CRCs, LNM-positive T1 CRC showed that the expression of genes involved in immune-related pathways such as B2M, HLA-B, and HLA-E were significantly downregulated in the distal stalk region compared to the proximal stalk region. In summary, our results may tentatively suggest considering endoscopic resection of the stalk with a minimum 2000 μm margin from the HL, taking into account the gene expression alterations related to immune-related pathways. However, we acknowledge the limitations of this pilot study, notably the small case series, which may restrict the depth of interpretation. Further validation is imperative to substantiate these findings.

Multiomics Integration at Single-Cell Resolution Using Bayesian Networks: A Case Study in Hepatocellular Carcinoma

Multiomics data integration is one of the leading frontiers of complex disease research and integrative biology. The advances in single-cell sequencing technologies offer yet another crucial dimension in multiomics research. The single-cell studies enable the study and integration of multiomics data simultaneously in the same cell. We report in this study multiomics data integration in single-cell resolution using Bayesian networks (BNs) in a case study of hepatocellular carcinoma (HCC). A BN encodes the conditional dependencies/independencies of variables using a graphical model with an accompanying joint probability. RNA-seq and Reduced Representation Bisulfite Sequencing data were analyzed separately, and copy number variations were estimated by the hidden Markov model method. Several BN models were constructed to reveal omics' causal and associational relationships. These methods were subjected to a validation study using an independent data set. We show the heterogeneity of the multiple cellular layers of HCC at single-cell omics resolution by identifying best-fitted BN models of 295 genes. We also provide novel insights into the multiomics mechanistic relationships in the human lymphocyte antigen class I genes in HCC. To the best of our knowledge, this is the first study to focus on integrating omics data using a machine learning algorithm, BNs, at the single-cell resolution using a case study of HCC.

Establishment of HLA class I and MICA/B null HEK-293T panel expressing single MICA alleles to detect anti-MICA antibodies

Pre- and post-transplantation anti-MICA antibody detection development are associated with an increased rejection risk and low graft survival. We previously generated HLA class I null HEK-293T using CRISPR/Cas9, while MICA and MICB genes were removed in this study. A panel of 11 cell lines expressing single MICA alleles was established. Anti-MICA antibody in the sera of kidney transplant patients was determined using flow cytometric method (FCM) and the Luminex method. In the 44 positive sera, the maximum FCM value was 2879 MFI compared to 28,135 MFI of Luminex method. Eleven sera (25%) were determined as positive by FCM and 32 sera (72%) were positive by the Luminex method. The sum of total MICA antigens, MICA*002, *004, *009, *019, and *027 correlation showed a statistically significant between the two methods (P = 0.0412, P = 0.0476, P = 0.0019, P = 0.0098, P = 0.0467, and P = 0.0049). These results demonstrated that HEK-293T-based engineered cell lines expressing single MICA alleles were suitable for measuring specific antibodies against MICA antigens in the sera of transplant patients. Studies of antibodies to MICA antigens may help to understand responses in vivo and increase clinical relevance at the cellular level such as complement-dependent cytotoxicity.

Main insights of genome wide association studies into HCV-related HCC

Background

Hepatocellular carcinoma (HCC) is one of the most common causes of cancer-mortality globally. Hepatocarcinogenesis is a complex multifactorial process. Host genetic background appeared to play a crucial role in the progression of HCC among chronic hepatitis C patients, especially in the era of Genome Wide Association Studies (GWAS) which allowed us to study the association of millions of single nucleotide polymorphisms (SNPs) with different complex diseases. This article aimed to review the discovered SNPs associated with the risk of HCV-related HCC development which was reported in the published GWA studies and subsequent validation studies and also try to explain the possible functional pathways.

Main text

We reviewed the recent GWA studies which reported several new loci associated with the risk of HCV-related HCC, such as (SNPs) in MHC class I polypeptide-related sequence A (MICA), DEP domain-containing 5 (DEPDC5), Tolloid-like protein 1 (TLL1), and human leukocyte antigen (HLA) genes. We also explained the possible underlying biological mechanisms that affect the host immune response pathways. Additionally, we discussed the controversial results reported by the subsequent validation studies of different ethnicities.

Conclusions

Although GWA studies reported strong evidence of the association between the identified SNPs and the risk of HCV-related HCC development, more functional experiments are necessary to confirm the defined roles of these genetic mutations for the future clinical application in different populations.

The role of Human leukocyte antigen class I on patient survival in Gastrointestinal cancers: a systematic review and meta- analysis

The prognostic role of Human leukocyte antigen class I (HLA- I) in gastrointestinal cancers has been remained controversial. We performed a meta-analysis to determine the role of classical HLA-I in predicting survival of patients. In addition, the relationship between HLA- I and some clinicopathological factors was evaluated. Published studies investigated HLA-I expression effect on gastrointestinal cancers were evaluated to determine association between HLA- I and overall survival (OS) and recurrence-free survival (RFS) in patients. The used effect sizes were hazard ratio (HR) and Odds ratio (OR) with 95% confidence interval (CI). A total of ten studies included 1307 patients were analyzed. The pooled results revealed that HLA- I overexpression was positively related to OS (HR: 0.72; 95% CI: 0.53–0.96) and demonstrated little association for RFS (HR: 0.70; 95% CI: 0.46–1.08). HLA-I overexpression is negative associated with poorer differentiation of tumor (OR: 0.53; 95% CI (0.43–0.81) and also higher stages of cancer (OR: 0.29; 95% CI (0.13–0.64). HLA- I overexpression was related to a better prognosis on OS and probably had little impact on RFS.

Novel genetic locus at MHC region for esophageal squamous cell carcinoma in Chinese populations

Background

Our previous genome-wide association study (GWAS) identified three independent single nucleotide polymorphisms (SNPs) in human major histocompatibility complex (MHC) region showing association with esophageal squamous cell carcinoma (ESCC). In this study, we increased GWAS sample size on MHC region and performed validation in an independent ESCC cases and normal controls with aim to find additional loci at MHC region showing association with an increased risk to ESCC.

Methods

The 1,077 ESCC cases and 1,733 controls were genotyped using Illumina Human 610-Quad Bead Chip, and 451 cases and 374 controls were genotyped using Illumina Human 660W-Quad Bead Chip. After quality control, the selected SNPs were replicated by TaqMan genotyping assay on another 2,026 ESCC cases and 2,384 normal controls.

Results

By excluding low quality SNPs in primary GWAS screening, we selected 2,533 SNPs in MHC region for association analysis, and identified 5 SNPs with p <10⁻⁴. Further validation analysis in an independent case-control cohort confirmed one of the 5 SNPs (rs911178) that showed significant association with ESCC. rs911178 (P_GWAS = 6.125E-04, OR = 0.644 and P_replication = 1.406E-22, OR = 0.489) was located at upstream of SCAND3.

Conclusion

The rs911178 (SCAND3 gene) in MHC region is significantly associated with high risk of ESCC. This study not only reveal the potential role of MHC region for the pathogenesis of ESCC, but also provides important clues for the establishment of tools and methods for screening high risk population of ESCC.

Mutation at positively selected positions in the binding site for HLA-C shows that KIR2DL1 is a more refined but less adaptable NK cell receptor than KIR2DL3

Through recognition of HLA class I, killer cell Ig-like receptors (KIR) modulate NK cell functions in human immunity and reproduction. Although a minority of HLA-A and -B allotypes are KIR ligands, HLA-C allotypes dominate this regulation, because they all carry either the C1 epitope recognized by KIR2DL2/3 or the C2 epitope recognized by KIR2DL1. The C1 epitope and C1-specific KIR evolved first, followed several million years later by the C2 epitope and C2-specific KIR. Strong, varying selection pressure on NK cell functions drove the diversification and divergence of hominid KIR, with six positions in the HLA class I binding site of KIR being targets for positive diversifying selection. Introducing each naturally occurring residue at these positions into KIR2DL1 and KIR2DL3 produced 38 point mutants that were tested for binding to 95 HLA- A, -B, and -C allotypes. Modulating specificity for HLA-C is position 44, whereas positions 71 and 131 control cross-reactivity with HLA-A*11:02. Dominating avidity modulation is position 70, with lesser contributions from positions 68 and 182. KIR2DL3 has lower avidity and broader specificity than KIR2DL1. Mutation could increase the avidity and change the specificity of KIR2DL3, whereas KIR2DL1 specificity was resistant to mutation, and its avidity could only be lowered. The contrasting inflexibility of KIR2DL1 and adaptability of KIR2DL3 fit with C2-specific KIR having evolved from C1-specific KIR, and not vice versa. Substitutions restricted to activating KIR all reduced the avidity of KIR2DL1 and KIR2DL3, further evidence that activating KIR function often becomes subject to selective attenuation.

Comparative genome analysis of the major histocompatibility complex (MHC) class I B/C segments in primates elucidated by genomic sequencing in common marmoset (Callithrix jacchus)

Common marmoset monkeys (Callithrix jacchus) have emerged as important animal models for biomedical research, necessitating a more extensive characterization of their major histocompatibility complex (MHC) region. However, the genomic information of the marmoset MHC (Caja) is still lacking. The MHC-B/C segment represents the most diverse MHC region among primates. Therefore, in this paper, to elucidate the detailed gene organization and evolutionary processes of the Caja class I B (Caja-B) segment, we determined two parts of the Caja-B sequences with 1,079 kb in total, ranging from H6orf15 to BAT1 and compared the structure and phylogeny with that of other primates. This segment contains 54 genes in total, nine Caja-B genes (Caja-B1 to Caja-B9), two MIC genes (MIC1 and MIC2), eight non-MHC genes, two non-coding genes, and 33 non-MHC pseudogenes that have not been observed in other primate MHC-B/C segments. Caja-B3, Caja-B4, and Caja-B7 encode proper MHC class I proteins according to amino acid structural characteristics. Phylogenetic relationships based on 48 MHC-I nucleotide sequences in primates suggested (1) species-specific divergence for Caja, Mamu, and HLA/Patr/Gogo lineages, (2) independent generation of the "seven coding exon" type MHC-B genes in Mamu and HLA/Patr/Gogo lineages from an ancestral "eight coding exon" type MHC-I gene, (3) parallel correlation with the long and short segmental duplication unit length in Caja and Mamu lineages. These findings indicate that the MHC-B/C segment has been under permanent selective pressure in the evolution of primates.

Methylation of the promoter of human leukocyte antigen class I in human esophageal squamous cell carcinoma and its histopathological characteristics

OBJECTIVE

The downregulation of human leukocyte antigen class I (HLA-I) has been proposed to contribute to the immune evasion of cancer cells in some cancers. Meanwhile, transcriptional silencing by means of promoter methylation is now believed to be an important mechanism of carcinogenesis. The aim of this study was (1) to examine the expression of HLA-I antigen and the antigen-processing machinery components in patients with esophageal squamous cell carcinoma lesions and (2) to detect the methylation pattern of the HLA-I gene in patients with esophageal squamous cell carcinoma and assess its association with histopathological characteristics.

METHODS

A total of 87 formalin-fixed and paraffin-embedded esophageal squamous cell carcinoma lesions were collected. HLA-I and antigen-processing machinery component expression was investigated by means of immunohistochemistry with anti-HLA-I monoclonal antibody, and methylation changes in the promoter region of HLA-1 genes were determined by using methylation-specific polymerase chain reaction.

RESULTS

HLA-I, transporter associated with antigen processing 1, and low molecular weight protein were lost or downregulated in 67%, 29.8%, and 47.0% of the esophageal squamous cell carcinoma lesions, respectively. The positive rates of gene promoter hypermethylation of HLA-I was 70.1% (61/87) in tumor tissues compared with 3.6% in normal tissue (P < .01). Also, the higher methylation rates and the HLA-I expression were significantly associated with tumor grade, including lymph node metastasis (P < .05).

CONCLUSIONS

HLA-I promoter hypermethylation was associated with loss of HLA-I antigen, which frequently occurred in primary tumors, especially in metastatic lymph node lesions, and was associated with patients' prognoses.

Role of gene methylation in antitumor immune response: implication for tumor progression

Cancer immunosurveillance theory has emphasized the role of escape mechanisms in tumor growth. In this respect, a very important factor is the molecular characterization of the mechanisms by which tumor cells evade immune recognition and destruction. Among the many escape mechanisms identified, alterations in classical and non-classical HLA (Human Leucocyte Antigens) class I and class II expression by tumor cells are of particular interest. In addition to the importance of HLA molecules, tumor-associated antigens and accessory/co-stimulatory molecules are also involved in immune recognition. The loss of HLA class I antigen expression and of co-stimulatory molecules can occur at genetic, transcriptional and post-transcriptional levels. Epigenetic defects are involved in at least some mechanisms that preclude mounting a successful host-antitumor response involving the HLA system, tumor-associated antigens, and accessory/co-stimulatory molecules. This review summarizes our current understanding of the role of methylation in the regulation of molecules involved in the tumor immune response.

Extended gene map reveals tripartite motif, C-type lectin, and Ig superfamily type genes within a subregion of the chicken MHC-B affecting infectious disease

MHC haplotypes have a remarkable influence on whether tumors form following infection of chickens with oncogenic Marek's disease herpesvirus. Although resistance to tumor formation has been mapped to a subregion of the chicken MHC-B region, the gene or genes responsible have not been identified. A full gene map of the subregion has been lacking. We have expanded the MHC-B region gene map beyond the 92-kb core previously reported for another haplotype revealing the presence of 46 genes within 242 kb in the Red Jungle Fowl haplotype. Even though MHC-B is structured differently, many of the newly revealed genes are related to loci typical of the MHC in other species. Other MHC-B loci are homologs of genes found within MHC paralogous regions (regions thought to be derived from ancient duplications of a primordial immune defense complex where genes have undergone differential silencing over evolutionary time) on other chromosomes. Still others are similar to genes that define the NK complex in mammals. Many of the newly mapped genes display allelic variability and fall within the MHC-B subregion previously shown to affect the formation of Marek's disease tumors and hence are candidates for genes conferring resistance.

MICA marks additional risk factors for Type 1 diabetes on extended HLA haplotypes: an association and meta-analysis

The association of the HLA complex on chromosome 6 does not explain total linkage of the HLA region to Type 1 Diabetes (T1D), leading to the hypothesis that there may be additional causal genes in the HLA region for immune-related disorders. Reports on the MHC Class I chain-related A (MICA) gene as candidate for association with T1D are contradicting. We investigated whether variation in MICA is associated to T1D in a cohort of 350 unrelated individuals with juvenile-onset T1D and 540 control subjects, followed by a meta-analysis of 14 studies. We also investigated an HLA-independent association for MICA with T1D. In our case-control study, we found that the MICA*A5 variant was significantly associated with an increased risk for T1D, while MICA*A6 was significantly associated with a decreased risk that was confirmed by our meta-analysis. However, the meta-analysis did not show an association of MICA*A5 T1D. Analysis of MICA alleles conditional on T1D-associated high-risk MHC class II haplotypes revealed that MICA*A6 was associated with an increased risk for T1D when this marker co-occurred with HLA DQ2DR17 T1D-risk-haplotypes. In contrast, MICA*A6 reduced the risk from the HLA DQ8DR4 T1D-risk haplotype. Moreover, MICA*A9 showed a significant association to increased risk for T1D on DQ8DR4 haplotypes. Co-inheritance of MICA*A6 with the HLA DQ2DR17 haplotype in T1D indicates this haplotype may carry the additional genetic factors for T1D, but our study does not support an independent association between MICA variants and T1D.

Construction of a nurse shark (Ginglymostoma cirratum) bacterial artificial chromosome (BAC) library and a preliminary genome survey

BACKGROUND

Sharks are members of the taxonomic class Chondrichthyes, the oldest living jawed vertebrates. Genomic studies of this group, in comparison to representative species in other vertebrate taxa, will allow us to theorize about the fundamental genetic, developmental, and functional characteristics in the common ancestor of all jawed vertebrates.

AIMS

In order to obtain mapping and sequencing data for comparative genomics, we constructed a bacterial artificial chromosome (BAC) library for the nurse shark, Ginglymostoma cirratum.

RESULTS

The BAC library consists of 313,344 clones with an average insert size of 144 kb, covering ~4.5 x 1010 bp and thus providing an 11-fold coverage of the haploid genome. BAC end sequence analyses revealed, in addition to LINEs and SINEs commonly found in other animal and plant genomes, two new groups of nurse shark-specific repetitive elements, NSRE1 and NSRE2 that seem to be major components of the nurse shark genome. Screening the library with single-copy or multi-copy gene probes showed 6-28 primary positive clones per probe of which 50-90% were true positives, demonstrating that the BAC library is representative of the different regions of the nurse shark genome. Furthermore, some BAC clones contained multiple genes, making physical mapping feasible.

CONCLUSION

We have constructed a deep-coverage, high-quality, large insert, and publicly available BAC library for a cartilaginous fish. It will be very useful to the scientific community interested in shark genomic structure, comparative genomics, and functional studies. We found two new groups of repetitive elements specific to the nurse shark genome, which may contribute to the architecture and evolution of the nurse shark genome.

MICA is associated with type 1 diabetes in the Belgian population, independent of HLA-DQ

To ascertain association of MICA with type 1 diabetes (T1D) in the Belgian population, well-characterized antibody-positive patients were analyzed for MICA transmembrane gene polymorphism in both an association study and a nuclear family study. The frequency of MICA5 was significantly increased in the T1D patient group (18%) compared with the control population (12%, OR=1.6, pc<10(-3)), whereas MICA9 was decreased (11% versus 16%, OR=0.7, pc<0.01). A p value<10(-3) for the association of MICA conditional on HLA class II and p=0.01 for the conditional extended transmission disequilibrium test were obtained, indicating that MICA is associated with type 1 diabetes, independent of HLA-DQ. Analysis of estimated extended HLA-DQ-MICA haplotypes revealed individual effects of MICA alleles. The most significant effect was seen for MICA5 on the HLA-DQA1*03-DQB1*0302-MICA haplotype (OR=2.5, p<10(-3)). A significant protective effect was seen for the combination of DQA1*01-DQB1*0602/3 and MICA5.1 (OR=0.3, p<10(-3)). However, patients stratified according to the presence or absence of the different MICA alleles did not differ in terms of age at onset, sex, or other diabetes-related clinical and epidemiological data. In conclusion, MICA is associated with type 1 diabetes in the Belgian population and the observed association does not result from the HLA-DQ associated risk.

Genomic sequence analysis of the 238-kb swine segment with a cluster of TRIM and olfactory receptor genes located, but with no class I genes, at the distal end of the SLA class I region

Continuous genomic sequence has been previously determined for the swine leukocyte antigen (SLA) class I region from the TNF gene cluster at the border between the major histocompatibility complex (MHC) class III and class I regions to the UBD gene at the telomeric end of the classical class I gene cluster (SLA-1 to SLA-5, SLA-9, SLA-11). To complete the genomic sequence of the entire SLA class I genomic region, we have analyzed the genomic sequences of two BAC clones carrying a continuous 237,633-bp-long segment spanning from the TRIM15 gene to the UBD gene located on the telomeric side of the classical SLA class I gene cluster. Fifteen non-class I genes, including the zinc finger and the tripartite motif (TRIM) ring-finger-related family genes and olfactory receptor genes, were identified in the 238-kilobase (kb) segment, and their location in the segment was similar to their apparent human homologs. In contrast, a human segment (alpha block) spanning about 375 kb from the gene ETF1P1 and from the HLA-J to HLA-F genes was absent from the 238-kb swine segment. We conclude that the gene organization of the MHC non-class I genes located in the telomeric side of the classical SLA class I gene cluster is remarkably similar between the swine and the human segments, although the swine lacks a 375-kb segment corresponding to the human alpha block.

DNA sequence based typing of swine leukocyte antigens in Yucatan Miniature Pigs

BACKGROUND

We have established a breeding program to develop additional lines of swine leukocyte antigens (SLA) homozygous miniature pigs derived from the Yucatan Miniature Pig. Yucatan pigs have been used extensively in biomedical research since the 1970s and are known for their docile nature and small size. The breed has no consanguinity with the Hormel or Pittman-Moore breeds used to produce the NIH Miniature Pigs, the only other SLA homozygous miniature pigs in the USA.

METHODS

SLA typing was initially done by restriction fragment length polymorphism. Then a cDNA library was constructed using spleen cells from these pigs, from which we cloned and sequenced nearly all alleles for the SLA-1, 3, 2, 6, DRB1, DQA and DQB1 loci.

RESULTS

Four SLA homozygous lines were established with haplotypes that we have designated 'w, x, y and z'. An SLA class I/II crossover haplotype, designated 'q', with the SLA class I alleles of 'w' and the SLA class II alleles of 'z' was discovered and used to establish a fifth line. The cDNA sequences were used to develop locus specific primers for each locus and an reverse transcription-polymerase chain reaction sequence based typing (SBT) method. We have used this method to perform SBT on SLA homozygous Yucatan pigs with these haplotypes and on the NIH pig 'a, c and d' haplotypes.

CONCLUSIONS

These pig lines represent a new resource for transplantation research and the methods we describe can be used to SLA type any herd of pigs.

ERVK9, transposons and the evolution of MHC class I duplicons within the alpha-block of the human and chimpanzee

The genomic sequences within the alpha-block (approximately 288-310 kb) of the human and chimpanzee MHC class I region contains ten MHC class I genes and three MIC gene fragments grouped together within alternating duplicated genomic segments or duplicons. In this study, the chimpanzee and human genomic sequences were analyzed in order to determine whether the remnants of the ERVK9 and other retrotransposon sequences are useful genomic markers for reconstructing the evolutionary history of the duplicated MHC gene families within the alpha-block. A variety of genes, pseudogenes, autologous DNA transposons and retrotransposons such as Alu and ERVK9 were used to categorize the ten duplicons into four distinct structural groups. The phylogenetic relationship of the ten duplicons was examined by using the neighbour joining method to analyze transposon sequence topologies of selected Alu members, LTR16B and Charlie9. On the basis of these structural groups and the phylogeny of the duplicated transposon sequences, a duplication model was reconstructed involving four multipartite tandem duplication steps to explain the organization and evolution of the ten duplicons within the alpha-block of the chimpanzee and human. The phylogenetic analysis and inferred duplication history suggests that the Patr/HLA-F was the first MHC class I gene to have been fixed and not required as a precursor for further duplication within the alpha-block of the ancestral species.

Genomic evolution of MHC class I region in primates

To elucidate the origins of the MHC-B-MHC-C pair and the MHC class I chain-related molecule (MIC)A-MICB pair, we sequenced an MHC class I genomic region of humans, chimpanzees, and rhesus monkeys and analyzed the regions from an evolutionary stand-point, focusing first on LINE sequences that are paralogous within each of the first two species and orthologous between them. Because all the long interspersed nuclear element (LINE) sequences were fragmented and nonfunctional, they were suitable for conducting phylogenetic study and, in particular, for estimating evolutionary time. Our study has revealed that MHC-B and MHC-C duplicated 22.3 million years (Myr) ago, and the ape MICA and MICB duplicated 14.1 Myr ago. We then estimated the divergence time of the rhesus monkey by using other orthologous LINE sequences in the class I regions of the three primate species. The result indicates that rhesus monkeys, and possibly the Old World monkeys in general, diverged from humans 27-30 Myr ago. Interestingly, rhesus monkeys were found to have not the pair of MHC-B and MHC-C but many repeated genes similar to MHC-B. These results support our inference that MHC-B and MHC-C duplicated after the divergence between apes and Old World monkeys.

Nucleotide sequencing analysis of the swine 433-kb genomic segment located between the non-classical and classical SLA class I gene clusters

Genome analysis of the swine leukocyte antigen ( SLA) region is needed to obtain information on the MHC genomic sequence similarities and differences between the swine and human, given the possible use of swine organs for xenotransplantation. Here, the genomic sequences of a 433-kb segment located between the non-classical and classical SLA class I gene clusters were determined and analyzed for gene organization and contents of repetitive sequences. The genomic organization and diversity of this swine non-class I gene region was compared with the orthologous region of the human leukocyte antigen ( HLA) complex. The length of the fully sequenced SLA genomic segment was 433 kb compared with 595 kb in the corresponding HLA class I region. This 162-kb difference in size between the swine and human genomic segments can be explained by indel activity, and the greater variety and density of repetitive sequences within the human MHC. Twenty-one swine genes with strong sequence similarity to the corresponding human genes were identified, with the gene order from the centromere to telomere of HCR - SPR1 - SEEK1 - CDSN - STG - DPCR1 - KIAA1885 - TFIIH - DDR - IER3 - FLOT1 - TUBB - KIAA0170 - NRM - KIAA1949 - DDX16 - FLJ13158 - MRPS18B - FB19 - ABCFI - CAT56. The human SEEK1 and DPCR1 genes are pseudogenes in swine. We conclude that the swine non-class I gene region that we have sequenced is highly conserved and therefore homologous to the corresponding region located between the HLA-C and HLA-E genes in the human.

Genomic organization of the mammalian MHC

The Human Genome Project transformed the quest of more than 50 years to understand the major histocompatibility complex (Mhc). The sequence of the Mhc from human and mouse, together with a large amount of sequence and mapping information from several other species, allows us to draw general conclusions about the organization and origin of this crucial part of the immune system. The Mhc is a mosaic of stretches formed by conserved and nonconserved genes. Surprisingly, of the approximately 3.6-Mb Mhc, the stretches that encode the class I and class II genes, which epitomize the Mhc, are the least conserved part, whereas the approximately 1.7-Mb stretches that encode at least 115 other genes are highly conserved. We summarize the available data to answer the questions (a) What is the Mhc? and (b) How can we define it in a general, not species-specific, way? Knowing what is essential and what is incidental helps us understand the fundamentals of the Mhc, and defining the species differences makes the model organisms more useful.

Comparative sequencing of human and chimpanzee MHC class I regions unveils insertions/deletions as the major path to genomic divergence

Despite their high degree of genomic similarity, reminiscent of their relatively recent separation from each other ( approximately 6 million years ago), the molecular basis of traits unique to humans vs. their closest relative, the chimpanzee, is largely unknown. This report describes a large-scale single-contig comparison between human and chimpanzee genomes via the sequence analysis of almost one-half of the immunologically critical MHC. This 1,750,601-bp stretch of DNA, which encompasses the entire class I along with the telomeric part of the MHC class III regions, corresponds to an orthologous 1,870,955 bp of the human HLA region. Sequence analysis confirms the existence of a high degree of sequence similarity between the two species. However, and importantly, this 98.6% sequence identity drops to only 86.7% taking into account the multiple insertions/deletions (indels) dispersed throughout the region. This is functionally exemplified by a large deletion of 95 kb between the virtual locations of human MICA and MICB genes, which results in a single hybrid chimpanzee MIC gene, in a segment of the MHC genetically linked to species-specific handling of several viral infections (HIV/SIV, hepatitis B and C) as well as susceptibility to various autoimmune diseases. Finally, if generalized, these data suggest that evolution may have used the mechanistically more drastic indels instead of the more subtle single-nucleotide substitutions for shaping the recently emerged primate species.

Hyperkeratosis and leukocytosis in transgenic mice carrying MHC class I chain-related gene B (MICB)

Major histocompatibility complex (MHC) class I chain-related gene A and B (MICA and MICB) are located very close to HLA-B. MICA is reported to be strongly associated with Behçet's disease (BD), a multisysytemic inflammation disorder characterized by oral apthous ulcers, skin lesions and genital ulcers. These two molecules are highly conserved at the amino acid levels. To determine the function of MICB in vivo and the relationship between the expression of MICB and BD experimentally, we produced several transgenic mouse lines (termed CAG-MICB) expressing human MICB cDNA under a ubiquitous promoter. They exhibited a 50% increase in the number of white blood cells compared with their non-transgenic littermates, and also exhibited a 10-20% reduction in body weight compared with non-transgenic littermates. Exfoliation of the skin first appeared around 7 days after birth and disappeared after 2 weeks of age. This was repeatedly observed in the transgenic offspring of two independent CAG-MICB lines examined. Histopathological analysis of skin of young mice exhibiting skin abnormalities revealed hyperkeratosis of the epidermis and thickening of the granular layer with slight infiltration of inflammatory cells in the dermis without any vasculitis. Other remarkable abnormalities associated with BD have not been observed in the CAG-MICB lines. Furthermore, fluorescein angiography of eyes of the CAG-MICB lines was performed, but there were no marked changes of BD-related uveitis in the ocular fundus. These findings suggest that (i) MICB expression is related to temporary skin inflammation, and (ii) expression of MICB is not directly associated with BD.

Comparative genomic analysis of the MHC: the evolution of class I duplication blocks, diversity and complexity from shark to man

The major histocompatibility complex (MHC) genomic region is composed of a group of linked genes involved functionally with the adaptive and innate immune systems. The class I and class II genes are intrinsic features of the MHC and have been found in all the jawed vertebrates studied so far. The MHC genomic regions of the human and the chicken (B locus) have been fully sequenced and mapped, and the mouse MHC sequence is almost finished. Information on the MHC genomic structures (size, complexity, genic and intergenic composition and organization, gene order and number) of other vertebrates is largely limited or nonexistent. Therefore, we are mapping, sequencing and analyzing the MHC genomic regions of different human haplotypes and at least eight nonhuman species. Here, we review our progress with these sequences and compare the human MHC structure with that of the nonhuman primates (chimpanzee and rhesus macaque), other mammals (pigs, mice and rats) and nonmammalian vertebrates such as birds (chicken and quail), bony fish (medaka, pufferfish and zebrafish) and cartilaginous fish (nurse shark). This comparison reveals a complex MHC structure for mammals and a relatively simpler design for nonmammalian animals with a hypothetical prototypic structure for the shark. In the mammalian MHC, there are two to five different class I duplication blocks embedded within a framework of conserved nonclass I and/or nonclass II genes. With a few exceptions, the class I framework genes are absent from the MHC of birds, bony fish and sharks. Comparative genomics of the MHC reveal a highly plastic region with major structural differences between the mammalian and nonmammalian vertebrates. Additional genomic data are needed on animals of the reptilia, crocodilia and marsupial classes to find the origins of the class I framework genes and examples of structures that may be intermediate between the simple and complex MHC organizations of birds and mammals, respectively.

Analysis of microsatellite polymorphism around the HLA-B locus in Iranian patients with Behçet's disease

We have previously suggested that in a Japanese population the susceptible locus for Behçet's disease (BD) is HLA-B51 itself. To confirm this finding in another population, we performed HLA class I typing using the PCR-SSP method and analyzed eight polymorphic markers distributed within 1100 kb around the HLA-B gene using automated sequencer and subsequent automated fragment detection by fluorescent-based technology with the DNA samples of 84 Iranian patients with BD and 87 healthy ethnically matched controls. As a result, three microsatellite alleles (MICA-A6, MIB-348, C1-4-1-217) and HLA-B51 were found to be strongly associated with BD. Of these alleles HLA-B51 is the most strongly associated allele. There were no alleles that were increased in allele frequency at any microsatellite loci centromeric of MICA or telomeric of HLA-B51. Therefore, HLA-B51 was confirmed to be by far the most strongly associated gene with BD in an Iranian population.

The absence of disease-specific polymorphisms within the HLA-B51 gene that is the susceptible locus for Behçet's disease

Behçet's disease is known to be associated with HLA-B51 in many different populations. Genetic evidence supports that the susceptible gene for Behçet's disease is the HLA-B51 allele at the HLA-B locus. This study was aimed to determine the HLA-B51 nucleotide sequence variation in three Behçet's disease patients and three healthy controls in order to elucidate if any disease specific mutations or polymorphisms may exist in the HLA-B51 gene of patients. Long-range polymerase chain reaction (PCR) was first carried out to give a PCR-amplified product of 9.5 kb which was then used as a template for nested PCR to give a final amplified product of 4.2 kb. This final product containing the 1.3-kb promoter/enhancer region and the entire HLA-B gene except for a 363-bp 3' terminal end segment encoding the 3' untranslated region was subcloned by the BP cloning technique and sequenced. The sequencing results showed that all the patients possessed the HLA-B*51011 allele, and there were no differences in the exonic nucleotide sequences between the three Behçet's disease patients and the three healthy controls. The HLA-B*51011 intronic and promoter/enhancer nucleotide sequences from the three patients had 22 single nucleotide polymorphisms (SNPs), a single insertion of 6 bp and a single deletion of 2 bp. On the other hand, the three healthy controls had 24 SNPs in their intronic and promoter/enhancer regions. However, none of these polymorphisms in the patients were specific for the disease. Therefore, these results clearly demonstrate that the HLA-B exonic sequence that encodes the HLA-B51 allele is the real pathogenic factor in Behçet's disease.

MICA and MICB genes: can the enigma of their polymorphism be resolved?

The human MHC class I chain-related genes (MICA and MICB) are located within the HLA class I region of chromosome 6. Their organization, expression and products differ considerably from classical HLA class I genes. MIC proteins are considered to be markers of "stress" in the epithelia, and act as ligands for cells expressing a common activatory natural killer-cell receptor (NKG2D). Molecular models are now available for the MICA protein, both bound and complexed with NKG2D. MICA molecules appear to be highly flexible and polymorphic, although the functional relevance and implications of their polymorphism have yet to be fully discerned.

Identification of MICA alleles with a long Leu-repeat in the transmembrane region and no cytoplasmic tail due to a frameshift-deletion in exon 4

MHC class I chain-related gene A (MICA) is located close to HLA-B gene and expressed in epithelial cells. The MICA gene is reported to be highly polymorphic as are the classical class I genes. To further assess the polymorphism in the MICA gene, we analyzed a total of 60 HLA-homozygous cells for the sequences spanning exons 2-6. In the analysis, four new MICA alleles were identified and six variations were recognized in exon 6. MICA*017, which was identified in three HLA-B57 homozygous cells (DBB, DEM and WIN), differed from MICA*002 in exon 3 and had a guanine deletion at the 3' end of exon 4. MICA*015 identified in an HLA-B45 homozygous cell (OMW) also had the same deletion that causes a frameshift mutation resulting in complete change of the transmembrane region and premature termination in the cytoplasmic tail; these alleles have a long hydrophobic leucine-rich region instead of the alanine repeat in the transmembrane region and terminate at the second position in the cytoplasmic domain. The frameshift deletion was found only in HLA-B45- or -B57-positive panels tested, suggesting a strong linkage disequilibrium between the deletion and B45 or B57. MICA*048, which was different in exon 5 from MICA*008, was identified in an HLA-B61 homozygous cell (TA21), while MICA*00901 identified in HLA-B51 homozygous cells (LUY and KT2) was distinguished from MICA*009 by exon 6.

Genomic anatomy of a premier major histocompatibility complex paralogous region on chromosome 1q21-q22

Human chromosomes 1q21-q25, 6p21.3-22.2, 9q33-q34, and 19p13.1-p13.4 carry clusters of paralogous loci, to date best defined by the flagship 6p MHC region. They have presumably been created by two rounds of large-scale genomic duplications around the time of vertebrate emergence. Phylogenetically, the 1q21-25 region seems most closely related to the 6p21.3 MHC region, as it is only the MHC paralogous region that includes bona fide MHC class I genes, the CD1 and MR1 loci. Here, to clarify the genomic structure of this model MHC paralogous region as well as to gain insight into the evolutionary dynamics of the entire quadriplication process, a detailed analysis of a critical 1.7 megabase (Mb) region was performed. To this end, a composite, deep, YAC, BAC, and PAC contig encompassing all five CD1 genes and linking the centromeric +P5 locus to the telomeric KRTC7 locus was constructed. Within this contig a 1.1-Mb BAC and PAC core segment joining CD1D to FCER1A was fully sequenced and thoroughly analyzed. This led to the mapping of a total of 41 genes (12 expressed genes, 12 possibly expressed genes, and 17 pseudogenes), among which 31 were novel. The latter include 20 olfactory receptor (OR) genes, 9 of which are potentially expressed. Importantly, CD1, SPTA1, OR, and FCERIA belong to multigene families, which have paralogues in the other three regions. Furthermore, it is noteworthy that 12 of the 13 expressed genes in the 1q21-q22 region around the CD1 loci are immunologically relevant. In addition to CD1A-E, these include SPTA1, MNDA, IFI-16, AIM2, BL1A, FY and FCERIA. This functional convergence of structurally unrelated genes is reminiscent of the 6p MHC region, and perhaps represents the emergence of yet another antigen presentation gene cluster, in this case dedicated to lipid/glycolipid antigens rather than antigen-derived peptides.

A second susceptibility gene for developing rheumatoid arthritis in the human MHC is localized within a 70-kb interval telomeric of the TNF genes in the HLA class III region

Rheumatoid arthritis (RA) is a chronic inflammatory joint disease with a multifactorial genetic basis. However, pathogenic genes for RA other than the human leukocyte antigen (HLA)-DRB1 gene have yet to be identified. Here, we investigated whether there is a second susceptibility locus for RA within the human major histocompatibility complex using 18 microsatellite markers distributed from the centromeric (HSET) to the telomeric end (P5-15) of the 3.6-Mb HLA region. Statistical studies of associated alleles on each microsatellite locus showed that one pathogenic gene for RA in the HLA region is localized in the DRB1 gene, as expected. Further, a second susceptibility gene of RA was suggested to be present in the HLA class III region, narrowed to 70 kb, that is just telomeric of the TNF gene cluster (TNFA and LTA) and that is located between the microsatellites TNFa and C1-2-A. In this critical segment, four expressed genes have been thus far identified, NFKBIL1 (IkappaBL), ATP6G, BAT1, and MICB, all of which are candidate genes for determining susceptibility to RA. These results exclude the possibility of involvement of the TNFA genes (TNF-alpha) in the development of RA, which was suggested previously to be a strong candidate for RA in the class III region.

Microsatellite mapping of a susceptible locus within the HLA region for Behçet's disease using Jordanian patients

Behçet's disease (BD) has been established to be associated with HLA-B51. However, it has not been revealed whether the HLA-B51 gene itself or another gene located near the HLA-B gene is directly involved in the pathogenesis of BD. Previously, using Japanese BD patients, our group has narrowed down a BD-causative gene to 46 kb between the MICA and HLA-B genes by means of fine mapping analysis with eight microsatellite markers distributed within a 1100 kb segment around the HLA-B gene. To know whether this mapping result is generally observed in BD of another population we have investigated repeat polymorphisms of the same microsatellite markers in Jordanian BD patients. Furthermore, we have evaluated these data by Mantel-Haenzel stratified analysis to find out a primarily associated locus for BD. As a result, HLA-B51 was found to be the most strongly and primarily associated marker. This result suggests that the pathogenic gene of BD is HLA-B51 itself, but unlikely to be other genes located in the vicinity of HLA-B.

Wide distribution of the MICA-MICB null haplotype in East Asians

Stress-inducible MICA (MHC class I chain-related A) is known to bind to NKG2D, which is one of the natural killer (NK) cell receptors, and plays a role in immune surveillance. We have reported that a MICA-MICB null haplotype is in linkage disequilibrium with HLA-B*4801 in the Japanese population. In the haplotype, an approximately 100-kb deletion, including the entire MICA gene, was observed and MICB possessed a premature stop codon. In this study, a multiplex polymerase chain reaction (PCR) method was developed for detecting the MICA deletion. MICB alleles were typed by PCR-single-strand conformation polymorphism (SSCP) method and direct sequencing. The frequency of the MICA-MICB null haplotype was 3.7% on the average, and was strongly associated with HLA-B48 in seven East Asian populations. It was presumed that the stop codon of MICB gene generated after the large-scale deletion. The wide distribution of the null haplotype at polymorphic frequencies suggests that the haplotype has been conservatively maintained because of some selective advantage.

Genomic organization, chromosomal localization, and the complete 22 kb DNA sequence of the human GCMa/GCM1, a placenta-specific transcription factor gene

The genomic sequence of the human GCMa/GCM1 gene, a mammalian homologue of Drosophila melanogaster GCM, was determined. Drosophila GCM is a neural transcription factor that regulates glial cell fate. The mammalian homolog however, is a placenta-specific transcription factor that is necessary for placental development. The 22 kb DNA sequence spanning the GCMa gene contains six exons and five introns, encoding a 2.8 kb cDNA. Overall genomic organization is similar for the human and mouse. Several potential binding sites for transcription factors like GATA, Oct-1, and bHLH proteins were found in the 5'-flanking region of the human gene. A DNA motif for GCM protein binding exists in the 5'-flanking region that is highly homologous with that of the mouse gene. The location of this gene was mapped to chromosome 6 using fluorescence in situ hybridization.

The implications of intergenic polymorphism for major histocompatibility complex evolution

A systematic survey of six intergenic regions flanking the human HLA-B locus in eight haplotypes reveals the regions to be up to 20 times more polymorphic than the reported average degree of human neutral polymorphism. Furthermore, the extent of polymorphism is directly related to the proximity to the HLA-B locus. Apparently linkage to HLA-B locus alleles, which are under balancing selection, maintains the neutral polymorphism of adjacent regions. For these linked polymorphisms to persist, recombination in the 200-kb interval from HLA-B to TNF must occur at a low frequency. The high degree of polymorphism found distal to HLA-B suggests that recombination is uncommon on both sides of the HLA-B locus. The least-squares estimate is 0.15% per megabase with an estimated range from 0.02 to 0.54%. These findings place strong restrictions on possible recombinational mechanisms for the generation of diversity at the HLA-B.

SNP profile within the human major histocompatibility complex reveals an extreme and interrupted level of nucleotide diversity

The human major histocompatibility complex (MHC) is characterized by polymorphic multicopy gene families, such as HLA and MIC (PERB11); duplications; insertions and deletions (indels); and uneven rates of recombination. Polymorphisms at the antigen recognition sites of the HLA class I and II genes and at associated neutral sites have been attributed to balancing selection and a hitchhiking effect, respectively. We, and others, have previously shown that nucleotide diversity between MHC haplotypes at non-HLA sites is unusually high (>10%) and up to several times greater than elsewhere in the genome (0.08%-0.2%). We report here the most extensive analysis of nucleotide diversity within a continuous sequence in the genome. We constructed a single nucleotide polymorphism (SNP) profile that reveals a pattern of extreme but interrupted levels of nucleotide diversity by comparing a continuous sequence within haplotypes in three genomic subregions of the MHC. A comparison of several haplotypes within one of the genomic subregions containing the HLA-B and -C loci suggests that positive selection is operating over the whole subgenomic region, including HLA and non-HLA genes. [The sequence data for the multiple haplotype comparisons within the class I region have been submitted to DDBJ/EMBL/GenBank under accession nos. AF029061, AF029062, and AB031005-AB031010. Additional sequence data have been submitted to the DDBJ data library under accession nos. AB031005-AB03101 and AF029061-AF029062.]

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.

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.

Association analysis using refined microsatellite markers localizes a susceptibility locus for psoriasis vulgaris within a 111 kb segment telomeric to the HLA-C gene

The HLA-Cw6 antigen has been associated with psoriasis vulgaris despite racial and ethnic differences. However, it remains unclear whether it is the HLA-Cw6 antigen itself or a closely linked, hitherto unidentified, locus that predisposes to the disease. Here, in order to map the susceptibility locus for psoriasis vulgaris precisely within the HLA class I region, 11 polymorphic microsatellite markers distributed throughout a 1060 kb segment surrounding the HLA-C locus were subjected to association analysis in Japanese psoriasis vulgaris patients. Statistical analyses of the distribution and deviation from Hardy-Weinberg equilibrium of the allelic frequency at each micro-satellite locus revealed that the pathogenic gene for psoriasis vulgaris is located within a reduced interval of 111 kb spanning 89-200 kb telomeric of the HLA-C gene. In addition to three known genes, POU5F1, TCF19 and S, this 111 kb fragment contains four new, expressed genes identified in the course of our genomic sequencing of the entire HLA class I region. Therefore, these seven genes are the potential candidates for susceptibility to psoriasis vulgaris.

Extensive nucleotide variability within a 370 kb sequence from the central region of the major histocompatibility complex

The recent availability of the genomic sequence spanning the central and telomeric end of the major histocompatibility complex (MHC) has allowed a detailed study of its organisation, gene content and level of nucleotide variability. Previous analyses of nucleotide variability in the MHC have focused on the coding regions of the human leukocyte antigen (HLA) Class I and II genes. Non-coding nucleotide variability has been considered a by-product of exonic diversity. However, with the advent of genomic sequencing, the extent of non-coding nucleotide variability within the MHC has just begun to be appreciated. In this study, we compared different human haplotypes in 370 kb of sequence in the central region of the MHC to show the following: 1. unusually high levels of non-coding nucleotide variability, up to 80 times greater than elsewhere in the genome; 2. non-coding nucleotide variability greater than 1% at nucleotide sites distant to the Class I genes; 3. nucleotide variability greater than 1% maintained over regions containing highly linked loci; and 4. distinct troughs and peaks in the level of nucleotide variability. We will discuss these observations in relation to a possible role of nucleotide variability in the organisation of the MHC.

New polymorphic microsatellite markers in the human MHC class I region

The human major histocompatibility complex (MHC) class I region is believed to contain a large number of genes encoding susceptible factors for diseases such as Behcet's disease, Graves disease and psoriasis vulgaris. To identify the causative genes of those diseases, we have conducted large-scale genomic sequencing and determined the 1.8 Mb entire HLA class I region from the MICB gene to the HLA-F gene. During the course of genomic sequencing, a total of 731 microsatellite sequences with dinucleotide to pentanucleotide repeats were found in this region. Previously, we reported that 26 microsatellites between MICB and S on the most centromeric side of the class I region, and between HSR1 and HLA-92/L in the midst of the class I region were highly polymorphic, and served as excellent genetic markers. In this paper, in order to fill the gaps with no known polymorphic microsatellites available in the HLA class I region, 12 new polymorphic microsatellite markers were recruited from the 1.8 Mb region including the remaining class I segments, namely between S and HSR1, and between HLA-92/L and HLA-F The average number of alleles at these new microsatellite loci was 8.2 with a polymorphism content value (PIC) of 0.63. These 38 markers in total almost uniformly interspersed in the HLA class I region will enable us to search precisely for the location of disease susceptible loci within the HLA class I region by association and for linkage analyses.

Genetic mapping and physical cloning of UVB susceptibility region in mice

One of the most important goals of cancer research is to identify environmental and host factors that contribute to the malignant state. Human skin cancers are among the few tumor types for which the predominant environmental carcinogen is known. Ultraviolet light, a component of sunlight, is an important cause of skin cancer in humans. In humans and mice, ultraviolet B radiation induces systematic and local immunosuppression. A consequence of that is inappropriate immune surveillance of somatic tissues for evidence of malignantly transformed cells. The impairment of contact hypersensitivity, as it develops early and correlates well with tumor frequency in various mouse strains, has been used for over 15 y as a model of immunologic events occurring in photocarcinogenesis. In mice, as well as in humans, ultraviolet B radiation induced impairment of contact hypersensitivity is not uniform in all individuals; some individuals are susceptible to the deleterious effects of ultraviolet B, whereas others are resistant to ultraviolet B. We have defined the genetic locus responsible for ultraviolet B susceptibility and resistance in mice within the Bat5 and H-2D segment of the mouse chromosome 17.

Sequencing based typing for genetic polymorphisms in exons, 2, 3 and 4 of the MICA gene

We have established a sequencing based typing (SBT) method for detection of genetic polymorphism in the exon 2 to 4 domains of the major histocompatibility complex (MHC) class I chain-related gene A (MICA) and applied it to allele typing of 130 healthy Japanese individuals. A 2.2-kb segment including exons 2, 3 and 4 of the MICA gene was amplified by a pair of generic primers followed by cycle sequencing using exon-specific nested primers. In total, 8 alleles were observed in a Japanese population and the most frequent allele was MICA008 with the gene frequency of 30.8%. MICA009 was the second most frequent (16.5%), while the rarest one was MICA007 (1.2%). MICA alleles displayed strong linkage equilibria with HLA-B antigens (i.e. MICA008 with B7, B48, B60 and B61; MICA009 with B51 and B52; MICA002 with B35, B39, B58 and B67; MICA004 with B44, MICA007 with B13 and B27; MICA010 with B46, B62 and B48, MICA012 with B54, B55, B56 and B59; MICA019 and B70, B71 and B62). Recently, the B48 haplotype has been reported to lack the entire MICA gene by a large-scale deletion in a Japanese population. Among 8 serologically B48 homozygous individuals, 4 were found to represent this MICA null allele as assessed by no polymerase chain reaction (PCR) amplification using MICA-specific primers, while the remaining four possessed the intact MICA gene with MICA008 or MICA010.

The critical region for Behçet disease in the human major histocompatibility complex is reduced to a 46-kb segment centromeric of HLA-B, by association analysis using refined microsatellite mapping

The HLA-B51 allele is known to be associated with Behçet disease. Recently, we found a higher risk for Behçet disease in the MICA gene, 46 kb centromeric of HLA-B, by investigation of GCT repetitive polymorphism within exon 5 of MICA. The pathogenic gene causing Behçet disease, however, has remained uncertain. Here, eight polymorphic microsatellite markers, distributed over a 900-kb region surrounding the HLA-B locus, were subjected to association analysis for Behçet disease. Statistical studies of associated alleles detected on each microsatellite locus showed that the pathogenic gene for Behçet disease is most likely found within a 46-kb segment between the MICA and HLA-B genes. The results of this mapping study, and the results of an earlier study of ours, suggest that MICA is a strong candidate gene for the development of Behçet disease.

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.