Linkage disequilibrium and haplotype blocks in the MHC vary in an HLA haplotype specific manner assessed mainly by DRB1*03 and DRB1*04 haplotypes

Predicting the molecular mechanism-driven progression of breast cancer through comprehensive network pharmacology and molecular docking approach

Identification of key regulators is a critical step toward discovering biomarker that participate in BC. A gene expression dataset of breast cancer patients was used to construct a network identifying key regulators in breast cancer. Overexpressed genes were identified with BioXpress, and then curated genes were used to construct the BC interactome network. As a result of selecting the genes with the highest degree from the BC network and tracing them, three of them were identified as novel key regulators, since they were involved at all network levels, thus serving as the backbone. There is some evidence in the literature that these genes are associated with BC. In order to treat BC, drugs that can simultaneously interact with multiple targets are promising. When compared with single-target drugs, multi-target drugs have higher efficacy, improved safety profile, and are easier to administer. The haplotype and LD studies of the FN1 gene revealed that the identified variations rs6707530 and rs1250248 may both cause TB, and endometriosis respectively. Interethnic differences in SNP and haplotype frequencies might explain the unpredictability in association studies and may contribute to predicting the pharmacokinetics and pharmacodynamics of drugs using FN1.

MICA Polymorphism and Genetic Predisposition to T1D in Jordanian Patients: A Case-Control Study

Type 1 diabetes (T1D) is an autoimmune disorder whose etiology includes genetic and environmental factors. The non-classical Major Histocompatibility Complex (MHC) class I chain-related gene A (MICA) gene has been associated with increased susceptibility to T1D as the interaction of MICA to the Natural Killer Group 2D (NK2GD) receptors found on the cell surface of natural killer (NK) cells and T cells is responsible for inducing immune responses. MICA polymorphisms were reported in association with T1D among different ethnic groups. However, data from different populations revealed conflicting results, so the association of MICA polymorphisms with predisposition to T1D remains uncertain. The aim of this sequencing-based study was to identify, for the first time, the possible MICA alleles and/or genotypes that could be associated with T1D susceptibility in the Jordanian population. Polymorphisms in exons 2–4 and the short tandem repeats (STR) in exon 5 of the highly polymorphic MICA gene were analyzed. No evidence for association between T1D and MICA alleles/genotypes was found in this study, except for the MICA*011 allele which was found to be negatively associated with T1D (p = 0.023, OR = 0.125). In conclusion, MICA polymorphisms seem not to be associated with increasing T1D susceptibility in Jordanian patients.

Identification, analysis of deleterious SNPs of the human GSR gene and their effects on the structure and functions of associated proteins and other diseases

Hereditary glutathione reductase deficiency, caused by mutations of the GSR gene, is an autosomal recessive disorder characterized by decreased glutathione disulfide (GSSG) reduction activity and increased thermal instability. This study implemented computational analysis to screen the most likely mutation that might be associated with hereditary glutathione reductase deficiency and other diseases. Using ten online computational tools, the study revealed four nsSNPs among the 17 nsSNPs identified as most deleterious and disease associated. Structural analyses and evolutionary confirmation study of native and mutant GSR proteins using the HOPE project and ConSruf. HOPE revealed more flexibility in the native GSR structure than in the mutant structure. The mutation in GSR might be responsible for changes in the structural conformation and function of the GSR protein and might also play a significant role in inducing hereditary glutathione reductase deficiency. LD and haplotype studies of the gene revealed that the identified variations rs2978663 and rs8190955 may be responsible for obstructive heart defects (OHDs) and hereditary anemia, respectively. These interethnic differences in the frequencies of SNPs and haplotypes might help explain the unpredictability that has been reported in association studies and can contribute to predicting the pharmacokinetics and pharmacodynamics of drugs that make use of GSR.

MS-Based HLA-II Peptidomics Combined With Multiomics Will Aid the Development of Future Immunotherapies

Immunotherapies have emerged to treat diseases by selectively modulating a patient's immune response. Although the roles of T and B cells in adaptive immunity have been well studied, it remains difficult to select targets for immunotherapeutic strategies. Because human leukocyte antigen class II (HLA-II) peptides activate CD4+ T cells and regulate B cell activation, proliferation, and differentiation, these peptide antigens represent a class of potential immunotherapy targets and biomarkers. To better understand the molecular basis of how HLA-II antigen presentation is involved in disease progression and treatment, systematic HLA-II peptidomics combined with multiomic analyses of diverse cell types in healthy and diseased states is required. For this reason, MS-based innovations that facilitate investigations into the interplay between disease pathologies and the presentation of HLA-II peptides to CD4+ T cells will aid in the development of patient-focused immunotherapies.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

HLA Heterozygote Advantage against HIV-1 Is Driven by Quantitative and Qualitative Differences in HLA Allele-Specific Peptide Presentation

Pathogen-mediated balancing selection is regarded as a key driver of host immunogenetic diversity. A hallmark for balancing selection in humans is the heterozygote advantage at genes of the human leukocyte antigen (HLA), resulting in improved HIV-1 control. However, the actual mechanism of the observed heterozygote advantage is still elusive. HLA heterozygotes may present a broader array of antigenic viral peptides to immune cells, possibly resulting in a more efficient cytotoxic T-cell response. Alternatively, heterozygosity may simply increase the chance to carry the most protective HLA alleles, as individual HLA alleles are known to differ substantially in their association with HIV-1 control. Here, we used data from 6,311 HIV-1-infected individuals to explore the relative contribution of quantitative and qualitative aspects of peptide presentation in HLA heterozygote advantage against HIV. Screening the entire HIV-1 proteome, we observed that heterozygous individuals exhibited a broader array of HIV-1 peptides presented by their HLA class I alleles. In addition, viral load was negatively correlated with the breadth of the HIV-1 peptide repertoire bound by an individual's HLA variants, particularly at HLA-B. This suggests that heterozygote advantage at HLA-B is at least in part mediated by quantitative peptide presentation. We also observed higher HIV-1 sequence diversity among HLA-B heterozygous individuals, suggesting stronger evolutionary pressure from HLA heterozygosity. However, HLA heterozygotes were also more likely to carry certain HLA alleles, including the highly protective HLA-B*57:01 variant, indicating that HLA heterozygote advantage ultimately results from a combination of quantitative and qualitative effects in antigen presentation.

What has GWAS done for HLA and disease associations?

The major histocompatibility complex (MHC) is located in chromosome 6p21 and contains crucial regulators of immune response, including human leucocyte antigen (HLA) genes, alongside other genes with nonimmunological roles. More recently, a repertoire of noncoding RNA genes, including expressed pseudogenes, has also been identified. The MHC is the most gene dense and most polymorphic part of the human genome. The region exhibits haplotype-specific linkage disequilibrium patterns, contains the strongest cis- and trans-eQTLs/meQTLs in the genome and is known as a hot spot for disease associations. Another layer of complexity is provided to the region by the extreme structural variation and copy number variations. While the HLA-B gene has the highest number of alleles, the HLA-DR/DQ subregion is structurally most variable and shows the highest number of disease associations. Reliance on a single reference sequence has complicated the design, execution and analysis of GWAS for the MHC region and not infrequently, the MHC region has even been excluded from the analysis of GWAS data. Here, we contrast features of the MHC region with the rest of the genome and highlight its complexities, including its functional polymorphisms beyond those determined by single nucleotide polymorphisms or single amino acid residues. One of the several issues with customary GWAS analysis is that it does not address this additional layer of polymorphisms unique to the MHC region. We highlight alternative approaches that may assist with the analysis of GWAS data from the MHC region and unravel associations with all functional polymorphisms beyond single SNPs. We suggest that despite already showing the highest number of disease associations, the true extent of the involvement of the MHC region in disease genetics may not have been uncovered.

Genetics of primary sclerosing cholangitis and pathophysiological implications

Primary sclerosing cholangitis (PSC) is a chronic disease leading to fibrotic scarring of the intrahepatic and extrahepatic bile ducts, causing considerable morbidity and mortality via the development of cholestatic liver cirrhosis, concurrent IBD and a high risk of bile duct cancer. Expectations have been high that genetic studies would determine key factors in PSC pathogenesis to support the development of effective medical therapies. Through the application of genome-wide association studies, a large number of disease susceptibility genes have been identified. The overall genetic architecture of PSC shares features with both autoimmune diseases and IBD. Strong human leukocyte antigen gene associations, along with several susceptibility genes that are critically involved in T-cell function, support the involvement of adaptive immune responses in disease pathogenesis, and position PSC as an autoimmune disease. In this Review, we survey the developments that have led to these gene discoveries. We also elaborate relevant interpretations of individual gene findings in the context of established disease models in PSC, and propose relevant translational research efforts to pursue novel insights.

Genome-wide Association Study Identifies 27 Loci Influencing Concentrations of Circulating Cytokines and Growth Factors

Circulating cytokines and growth factors are regulators of inflammation and have been implicated in autoimmune and metabolic diseases. In this genome-wide association study (GWAS) of up to 8,293 Finns we identified 27 genome-widely significant loci (p < 1.2 × 10^-9) for one or more cytokines. Fifteen of the associated variants had expression quantitative trait loci in whole blood. We provide genetic instruments to clarify the causal roles of cytokine signaling and upstream inflammation in immune-related and other chronic diseases. We further link inflammatory markers with variants previously associated with autoimmune diseases such as Crohn disease, multiple sclerosis, and ulcerative colitis and hereby elucidate the molecular mechanisms underpinning these diseases and suggest potential drug targets.

Human transcriptome response to immunization with live-attenuated Venezuelan equine encephalitis virus vaccine (TC-83): Analysis of whole blood

Venezuelan equine encephalitis virus (VEEV) is an important human and animal alphavirus pathogen transmitted by mosquitoes. The virus is endemic in Central and South America, but has also caused equine outbreaks in southwestern areas of the United States. In an effort to better understand the molecular mechanisms of the development of immunity to this important pathogen, we performed transcriptional analysis from whole, unfractionated human blood of patients who had been immunized with the live-attenuated vaccine strain of VEEV, TC-83. We compared changes in the transcriptome between naïve individuals who were mock vaccinated with saline to responses of individuals who received TC-83. Significant transcriptional changes were noted at days 2, 7, and 14 following vaccination. The top canonical pathways revealed at early and intermediate time points (days 2 and 7) included the involvement of the classic interferon response, interferon-response factors, activation of pattern recognition receptors, and engagement of the inflammasome. By day 14, the top canonical pathways included oxidative phosphorylation, the protein ubiquitination pathway, natural killer cell signaling, and B-cell development. Biomarkers were identified that differentiate between vaccinees and control subjects, at early, intermediate, and late stages of the development of immunity as well as markers which were common to all 3 stages following vaccination but distinct from the sham-vaccinated control subjects. The study represents a novel examination of molecular processes that lead to the development of immunity against VEEV in humans and which may be of value as diagnostic targets, to enhance modern vaccine design, or molecular correlates of protection.

Mapping the genetic diversity of HLA haplotypes in the Japanese populations

Japan has often been viewed as an Asian country that possesses a genetically homogenous community. The basis for partitioning the country into prefectures has largely been geographical, although cultural and linguistic differences still exist between some of the districts/prefectures, especially between Okinawa and the mainland prefectures. The Major Histocompatibility Complex (MHC) region has consistently emerged as the most polymorphic region in the human genome, harbouring numerous biologically important variants; nevertheless the presence of population-specific long haplotypes hinders the imputation of SNPs and classical HLA alleles. Here, we examined the extent of genetic variation at the MHC between eight Japanese populations sampled from Okinawa, and six other prefectures located in or close to the mainland of Japan, specifically focusing at the haplotypes observed within each population, and what the impact of any variation has on imputation. Our results indicated that Okinawa was genetically farther to the mainland Japanese than were Gujarati Indians from Tamil Indians, while the mainland Japanese from six prefectures were more homogeneous than between northern and southern Han Chinese. The distribution of haplotypes across Japan was similar, although imputation was most accurate for Okinawa and several mainland prefectures when population-specific panels were used as reference.

Intrahaplotypic Variants Differentiate Complex Linkage Disequilibrium within Human MHC Haplotypes

Distinct regions of long-range genetic fixation in the human MHC region, known as conserved extended haplotypes (CEHs), possess unique genomic characteristics and are strongly associated with numerous diseases. While CEHs appear to be homogeneous by SNP analysis, the nature of fine variations within their genomic structure is unknown. Using multiple, MHC-homozygous cell lines, we demonstrate extensive sequence conservation in two common Asian MHC haplotypes: A33-B58-DR3 and A2-B46-DR9. However, characterization of phase-resolved MHC haplotypes revealed unique intra-CEH patterns of variation and uncovered 127 single nucleotide variants (SNVs) which are missing from public databases. We further show that the strong linkage disequilibrium structure within the human MHC that typically confounds precise identification of genetic features can be resolved using intra-CEH variants, as evidenced by rs3129063 and rs448489, which affect expression of ZFP57, a gene important in methylation and epigenetic regulation. This study demonstrates an improved strategy that can be used towards genetic dissection of diseases.

Genotyping of human leukocyte antigen (HLA) ancestral haplotypes as prognostic marker in cancer using PCR analysis

The major histocompatibility complex (MHC) comprises a set of genes that are essential to immunity and surveillance against neoplastic transformation. MHC antigens not only regulate antitumor immune responses in experimental animal models but also directly correlate with survival and prognosis of patients with various types of cancers. Effective recognition of tumor cells by effector T cells may be affected by the genotype and the extent of expression of human leukocyte antigen (HLA)-peptide complexes. Therefore, MHC antigens may serve as potential biomarkers for prognosis and allow selection of cancer patients for specific therapy. We describe PCR-based method to determine the HLA genotype in healthy individuals and patients using blood and tumor tissue as DNA source.

Population-specific recombination sites within the human MHC region

Genetic rearrangement by recombination is one of the major driving forces for genome evolution, and recombination is known to occur in non-random, discreet recombination sites within the genome. Mapping of recombination sites has proved to be difficult, particularly, in the human MHC region that is complicated by both population variation and highly polymorphic HLA genes. To overcome these problems, HLA-typed individuals from three representative populations: Asian, European and African were used to generate phased HLA haplotypes. Extended haplotype homozygosity (EHH) plots constructed from the phased haplotype data revealed discreet EHH drops corresponding to recombination events and these signatures were observed to be different for each population. Surprisingly, the majority of recombination sites detected are unique to each population, rather than being common. Unique recombination sites account for 56.8% (21/37 of total sites) in the Asian cohort, 50.0% (15/30 sites) in Europeans and 63.2% (24/38 sites) in Africans. Validation carried out at a known sperm typing recombination site of 45 kb (HLA-F-telomeric) showed that EHH was an efficient method to narrow the recombination region to 826 bp, and this was further refined to 660 bp by resequencing. This approach significantly enhanced mapping of the genomic architecture within the human MHC, and will be useful in studies to identify disease risk genes.

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.

Reproducible association with type 1 diabetes in the extended class I region of the major histocompatibility complex

The high-risk human leukocyte antigen (HLA)-DRB1, DQA1 and DQB1 alleles cannot explain the entire type 1 diabetes (T1D) association observed within the extended major histocompatibility complex. We have earlier identified an association with D6S2223, located 2.3 Mb telomeric of HLA-A, on the DRB1(*)03-DQA1(*)0501-DQB1(*)0201 haplotype, and this study aimed to fine-map the associated region also on the DRB1(*)0401-DQA1(*)03-DQB1(*)0302 haplotype, characterized by less extensive linkage disequilibrium. To exclude associations secondary to DRB1-DQA1-DQB1 haplotypes, 205 families with at least one parent homozygous for these loci, were genotyped for 137 polymorphisms. We found novel associations on the DRB1(*)0401-DQA1(*)03-DQB1(*)0302 haplotypic background with eight single nucleotide polymorphisms (SNPs) located within or near the PRSS16 gene. In addition, association at the butyrophilin (BTN)-gene cluster, particularly the BTN3A2 gene, was observed by multilocus analyses. We replicated the associations with SNPs in the PRSS16 region and, albeit weaker, to the BTN3A2 region, in an independent material of 725 families obtained from the Type 1 Diabetes Genetics Consortium. It is important to note that these associations were independent of the HLA-DRB1-DQA1-DQB1 genes, as well as of associations observed at HLA-A, -B and -C. Taken together, our results identify PRSS16 and BTN3A2, two genes thought to play important roles in regulating the immune response, as potentially novel susceptibility genes for T1D.

Three microsatellites from the T1DGC MHC data set show highly significant association with type 1 diabetes, independent of the HLA-DRB1, -DQA1 and -DQB1 genes

AIM

The aim of this study was to test the microsatellites in the Type 1 Diabetes Genetics Consortium major histocompatibility complex (MHC) data set for association with type 1 diabetes (T1D) independent of the HLA-DRB1, -DQA1 and -DQB1 genes.

METHODS

The data set was edited to contain only one affected child per family, and broad ethnic subgroups were defined. Genotypes for HLA-DRB1, -DQA1 and -DQB1 were replaced by a haplotype code spanning all three loci, with phase inferred based on common haplotypes. The final data set contained 8190 samples in 2301 families, 59 microsatellites and the DRB1-DQA1-DQB1 haplotype code. Statistical analyses consisted of conditional logistic regression and haplotype estimations and linkage disequilibrium calculations.

RESULTS

The data set was screened using a main effects test approach adjusted for DRB1-DQA1-DQB1, and significant results tested for validity. After these procedures, four markers remained significant at the Bonferroni-corrected threshold: D6S2773 (p = 0.00014), DG6S185 (p = 0.00015), DG6S398 (p = 0.00043) and D6S2998 (p = 0.00015). These results were supported by allelic tests conditioned on DRB1-DQA1-DQB1 haplotypes, except for DG6S185, which may contain artefacts.

CONCLUSIONS

We have identified three microsatellites that mark additional risk factors for T1D at highly significant levels in the MHC. Further analyses are needed to establish the relationship with other possible genetic determinants in this region.

The HLA genomic loci map: expression, interaction, diversity and disease

The human leukocyte antigen (HLA) super-locus is a genomic region in the chromosomal position 6p21 that encodes the six classical transplantation HLA genes and at least 132 protein coding genes that have important roles in the regulation of the immune system as well as some other fundamental molecular and cellular processes. This small segment of the human genome has been associated with more than 100 different diseases, including common diseases, such as diabetes, rheumatoid arthritis, psoriasis, asthma and various other autoimmune disorders. The first complete and continuous HLA 3.6 Mb genomic sequence was reported in 1999 with the annotation of 224 gene loci, including coding and non-coding genes that were reviewed extensively in 2004. In this review, we present (1) an updated list of all the HLA gene symbols, gene names, expression status, Online Mendelian Inheritance in Man (OMIM) numbers, including new genes, and latest changes to gene names and symbols, (2) a regional analysis of the extended class I, class I, class III, class II and extended class II subregions, (3) a summary of the interspersed repeats (retrotransposons and transposons), (4) examples of the sequence diversity between different HLA haplotypes, (5) intra- and extra-HLA gene interactions and (6) some of the HLA gene expression profiles and HLA genes associated with autoimmune and infectious diseases. Overall, the degrees and types of HLA super-locus coordinated gene expression profiles and gene variations have yet to be fully elucidated, integrated and defined for the processes involved with normal cellular and tissue physiology, inflammatory and immune responses, and autoimmune and infectious diseases.

Conditional analyses on the T1DGC MHC dataset: novel associations with type 1 diabetes around HLA-G and confirmation of HLA-B

The major histocompatibility complex (MHC) is known to harbour genetic risk factors for type 1 diabetes (T1D) additional to the class II determinants HLA-DRB1, -DQA1 and -DQB1, but strong linkage disequilibrium (LD) has made efforts to establish their location difficult. This study utilizes a dataset generated by the T1D genetics consortium (T1DGC), with genotypes for 2965 markers across the MHC in 2321 T1D families of multiple (mostly Caucasian) ethnicities. Using a comprehensive approach consisting of complementary conditional methods and LD analyses, we identified three regions with T1D association, independent both of the known class II determinants and of each other. A subset of polymorphisms that could explain most of the association in each region included single nucleotide polymorphisms (SNPs) in the vicinity of HLA-G, particular HLA-B and HLA-DPB1 alleles, and SNPs close to the COL11A2 and RING1 genes. Apart from HLA-B and HLA-DPB1, all of these represent novel associations, and subpopulation analyses did not indicate large population-specific differences among Caucasians for our findings. On account of the unusual genetic complexity of the MHC, further fine mapping is demanded, with the possible exception of HLA-B. However, our results mean that these efforts can be focused on narrow, defined regions of the MHC.

Genomic imprinting in the development and evolution of psychotic spectrum conditions

I review and evaluate genetic and genomic evidence salient to the hypothesis that the development and evolution of psychotic spectrum conditions have been mediated in part by alterations of imprinted genes expressed in the brain. Evidence from the genetics and genomics of schizophrenia, bipolar disorder, major depression, Prader-Willi syndrome, Klinefelter syndrome, and other neurogenetic conditions support the hypothesis that the etiologies of psychotic spectrum conditions commonly involve genetic and epigenetic imbalances in the effects of imprinted genes, with a bias towards increased relative effects from imprinted genes with maternal expression or other genes favouring maternal interests. By contrast, autistic spectrum conditions, including Kanner autism, Asperger syndrome, Rett syndrome, Turner syndrome, Angelman syndrome, and Beckwith-Wiedemann syndrome, commonly engender increased relative effects from paternally expressed imprinted genes, or reduced effects from genes favouring maternal interests. Imprinted-gene effects on the etiologies of autistic and psychotic spectrum conditions parallel the diametric effects of imprinted genes in placental and foetal development, in that psychotic spectrum conditions tend to be associated with undergrowth and relatively-slow brain development, whereas some autistic spectrum conditions involve brain and body overgrowth, especially in foetal development and early childhood. An important role for imprinted genes in the etiologies of psychotic and autistic spectrum conditions is consistent with neurodevelopmental models of these disorders, and with predictions from the conflict theory of genomic imprinting.

Genetic epidemiology of primary sclerosing cholangitis

The aetiology of primary sclerosing cholangitis (PSC) is not known. A more than 80-fold increased risk of PSC among first-degree relatives emphasizes the importance of genetic factors. Genetic associations within the human leukocyte antigen (HLA) complex on chromosome 6p21 were detected in PSC 25 years ago. Subsequent studies have substantiated beyond doubt that one or more genetic variants located within this genetic region are important. The true identities of these variants, however, remain to be identified. Several candidate genes at other chromosomal loci have also been investigated. However, according to strict criteria for what may be denominated a susceptibility gene in complex diseases, no such gene exists for PSC today. This review summarises present knowledge on the genetic susceptibility to PSC, as well as genetic associations with disease progression and clinical subsets of particular interest (inflammatory bowel disease and cholangiocarcinoma).

Haplotype-specific pattern of association of human major histocompatibility complex with non-Hodgkin's lymphoma outcome

In the previous studies, some human major histocompatibility complex (MHC) genes such as TNF, LTA and human leukocyte antigen (HLA)-DR2 genes and A1-B8-TNF(-308A) haplotype were implied in non-Hodgkin's lymphoma (NHL) outcome. In the current study, we have assigned most probable six-locus haplotypes determined by HLA-A, -Cw, -B and -DRB1 highly polymorphic genes and non-HLA LTA(+252) and TNF(-308) single nucleotide polymorphisms (SNPs) in 152 NHL Caucasian French patients. We have broadly mapped the MHC region by its component blocks and tagging alleles. Ten frequent (with haplotype frequency >1%) six-locus extended haplotypes (EHs) were revealed in NHL patients. The only two adjacent locus fragment of 8.1 EH associated with shortened freedom from progression (FFP) was B*08-LTA(+252G) (P= 0.0084, RR = 2.45). Interestingly, 305-kbp-long, four-locus fragment of 8.1 EH, Cw*07-B*08-LTA(+252G)-TNF(-308A) block was much strongly associated with shortened FFP (P= 0.00045, RR = 3.26). The analysis of further extended haploblocks comprising five or six loci showed weaker association with outcome measures, suggesting linkage disequilibrium to be the cause of DRB1*03 and A*01 allele associations. In contrast, all fragments of 7.1 EH influenced FFP favorably with top association of TNF(-308G) allele. In multivariate analysis, only Cw*07-B*08-LTA(+252G)-TNF(-308A) and TNF(-308G)-DRB1*01 haplotypes remained predictive for shortened FFP (P= 0.024 and 0.027, respectively) and independent of International Prognostic Index (P= 0.00044). This study reveals that the block composition of EHs may cause important functional differences for NHL outcomes. Further study will be required in NHL patients by fine mapping with dense microsatellite or SNP tags to define susceptibility genes in associating regions.

Association analysis in type 1 diabetes of the PRSS16 gene encoding a thymus-specific serine protease

We have previously mapped a separate type 1 diabetes (T1D) association in the extended MHC class I region, marked by D6S2223, on the DRB1*03-DQA1*0501-DQB1*0201 haplotype. The associated region encompasses a gene encoding a thymus-specific serine protease (PRSS16), presumably involved in positive selection of T cells or in T-cell regulation. Fourteen PRSS16 polymorphisms were genotyped in two steps using a total of six T1D family data sets, as well as case-control materials for both T1D and celiac disease (CD). An association with a 15 base-pair deletion in exon 12 of PRSS16 was found on the DRB1*03-DQA1*0501-DQB1*0201 haplotype for both T1D and CD, but it could not explain the more pronounced disease associations observed at marker D6S2223. We compared the performance of the 14 tested PRSS16 polymorphisms, selected after our previous comprehensive screen, against HapMap selected tag SNPs. Use of a HapMap based SNP selection strategy would result in loss of a large proportion of the genetic variation in PRSS16. Our data suggest that it is unlikely that polymorphisms within the PRSS16 gene are involved in the predisposition to T1D. However, we cannot rule out that regulatory polymorphisms located some distance away from the gene may be involved.

14th International HLA and Immunogenetics Workshop: report on the HLA component of type 1 diabetes

The type 1 diabetes (T1D) component of the 13th International Histocompatibility Workshop (IHW) obtained microsatellite (msat) and human leukocyte antigen (HLA)-DR/DQ data on case/control and family samples through an international collaboration. The aim was to detect the effects of susceptibility loci on the HLA complex independent of the primary determinants in the class II region (HLA-DR/DQ). As part of the activity of the 14th International HLA and Immunogenetics Workshop (14th IHIWS), a T1D workshop was held to present analyses of the 13th IHW data and to discuss the current status of knowledge about the genetics of T1D. These data are now available online through dbMHC, a web-based resource established by the National Center for Biotechnology. Continuing work since the 13th IHW has resulted in published work showing heterogeneity of DR3 haplotypes in data sets from the 13th IHW and Human Biological Data Interchange (HBDI). In addition, we identified markers that define DRB1*1501 DQB1*0602 haplotypes conferring reduced protection from diabetes in a Swedish 13th IHW data set. Further analyses of the 13th IHW data set not only showed some significant results but also demonstrated extensive heterogeneity reminiscent of non-HLA genes. The haplotype analysis in HBDI families identified two msats with significant effects on susceptibility and statistically significant age of onset effects at class III markers that are not because of linkage disequilibrium, with class I alleles known to affect age of onset. The above studies underscore the importance of refining our understanding of susceptibility associated with genes in the HLA complex.

Evaluation of association of SNPs in the TNF alpha gene region with schizophrenia

The association of the tumor necrosis factor alpha (TNFalpha) -G308A promoter polymorphism with schizophrenia has complemented clinical findings of increased levels of the TNFalpha cytokine in schizophrenic patients, with some support for a functional consequence of the variant. Our previous studies of genetic causes in schizophrenia supported findings of linkage to the major histocompatibility complex (MHC) region where the TNFalpha gene is located as well as association with the -G308A promoter polymorphism. While the common G-allele shows association in our sample, association with the A-allele has been reported by other groups. This suggests linkage disequilibrium (LD) rather than direct involvement in the disorder. In order to define LD of DNA variants with the disorder in this area, we analyzed 36 SNPs in a 165-kb region around this polymorphism. We detected nominally significant associations (P < 0.05) of three markers (including the -G308A promoter polymorphism) and multiple haplotypes with schizophrenia in our sample of 204 families (79 sib-pairs and 125 trios). The association is largely restricted to a 30 kb high LD region/block and should assist in the identification of a schizophrenia susceptibility gene within the block or elsewhere in the MHC.

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.

Conserved extended haplotypes discriminate HLA-DR3-homozygous Basque patients with type 1 diabetes mellitus and celiac disease

The major susceptibility locus for type 1 diabetes mellitus (T1D) maps to the human lymphocyte antigen (HLA) class II region in the major histocompatibility complex on chromosome 6p21. In southern European populations, like the Basques, the greatest risk to T1D is associated with DR3 homo- and heterozygosity and is comparable to that of DR3/DR4, the highest risk genotype in northern European populations. Celiac disease (CD) is another DR3-associated autoimmune disorder showing certain overlap with T1D that has been explained by the involvement of common genetic determinants, a situation more frequent in DR3-rich populations, like the Basques. As both T1D- and CD-associated HLA alleles are part of conserved extended haplotypes (CEH), we compared DR3-homozygous T1D and CD patients to determine whether CEHs were equally distributed between both disorders or there was a differential contribution of different haplotypes. We observed a very pronounced distribution bias (P<10(-5)) of the two major DR3 CEHs, with DR3-B18 predominating in T1D and DR3-B8 in CD. Additionally, high-density single nucleotide polymorphism (SNP) analysis of the complete CEH [A*30-B*18-MICA*4-F1C30-DRB1*0301-DQB1*0201-DPB1*0202] revealed extraordinary conservation throughout the 4.9 Mbp analyzed supporting the existence of additional diabetogenic variants (other than HLA-DRB1*0301-DQB1*0201), conserved within the DR3-B18 CEH (but not in other DR3 haplotypes) that could explain its enhanced diabetogenicity.