The Haplotype Structure of the Human Major Histocompatibility Complex

Exploring the HLA complex in autoimmunity: From the risk haplotypes to the modulation of expression

The genes mapping at the HLA region show high density, strong linkage disequilibrium and high polymorphism, which affect the association of HLA class I and class II genes with autoimmunity. We focused on the HLA haplotypes, genomic structures consisting of an array of specific alleles showing some degrees of genetic association with different autoimmune disorders. GWASs in many pathologies have identified variants in either the coding loci or the flanking regulatory regions, both in linkage disequilibrium in haplotypes, that are frequently associated with increased risk and may influence gene expression. We discuss the relevance of the HLA gene expression because the level of surface heterodimers determines the number of complexes presenting self-antigen and, thus, the strength of pathogenic autoreactive T cells immune response.

Integrative HLA typing of tumor and adjacent normal tissue can reveal insights into the tumor immune response

BACKGROUND

The HLA complex is the most polymorphic region of the human genome, and its improved characterization can help us understand the genetics of human disease as well as the interplay between cancer and the immune system. The main function of HLA genes is to recognize "non-self" antigens and to present them on the cell surface to T cells, which instigate an immune response toward infected or transformed cells. While sequence variation in the antigen-binding groove of HLA may modulate the repertoire of immunogenic antigens presented to T cells, alterations in HLA expression can significantly influence the immune response to pathogens and cancer.

METHODS

RNA sequencing was used here to accurately genotype the HLA region and quantify and compare the level of allele-specific HLA expression in tumors and patient-matched adjacent normal tissue. The computational approach utilized in the study types classical and non-classical Class I and Class II HLA alleles from RNA-seq while simultaneously quantifying allele-specific or personalized HLA expression. The strategy also uses RNA-seq data to infer immune cell infiltration into tumors and the corresponding immune cell composition of matched normal tissue, to reveal potential insights related to T cell and NK cell interactions with tumor HLA alleles.

RESULTS

The genotyping method outperforms existing RNA-seq-based HLA typing tools for Class II HLA genotyping. Further, we demonstrate its potential for studying tumor-immune interactions by applying the method to tumor samples from two different subtypes of breast cancer and their matched normal breast tissue controls.

CONCLUSIONS

The integrative RNA-seq-based HLA typing approach described in the study, coupled with HLA expression analysis, neoantigen prediction and immune cell infiltration, may help increase our understanding of the interplay between a patient's tumor and immune system; and provide further insights into the immune mechanisms that determine a positive or negative outcome following treatment with immunotherapy such as checkpoint blockade.

HLA Genetics for the Human Diseases

Highly polymorphic human leukocyte antigen (HLA) molecules (alleles) expressed by different classical HLA class I and class II genes have crucial roles in the regulation of innate and adaptive immune responses, transplant rejection and in the pathogenesis of numerous infectious and autoimmune diseases. To date, over 35,000 HLA alleles have been published from the IPD-IMGT/HLA database, and specific HLA alleles and HLA haplotypes have been reported to be associated with more than 100 different diseases and phenotypes. Next generation sequencing (NGS) technology developed in recent years has provided breakthroughs in various HLA genomic/gene studies and transplant medicine. In this chapter, we review the current information on the HLA genomic structure and polymorphisms, as well as the genetic context in which numerous disease associations have been identified in this region.

The Path to Conserved Extended Haplotypes: Megabase-Length Haplotypes at High Population Frequency

This minireview describes the history of the conceptual development of conserved extended haplotypes (CEHs): megabase-length haplotypes that exist at high (≥0.5%) population frequency. My career began in internal medicine, shifted to pediatrics, and clinical practice changed to research. My research interest was initially in hematology: on plasma proteins, their metabolism, synthesis, and function. This narrowed to a focus on proteins of the human complement system, their role in immunity and their genetics, beginning with polymorphism and deficiency of C3. My group identified genetic polymorphisms and/or inherited deficiencies of C2, C4, C6, and C8. After defining glycine-rich beta glycoprotein as factor B (Bf) in the properdin system, we found that the genes for Bf (CFB), C2, C4A, and C4B were inherited as a single haplotypic unit which we named the "complotype." Complotypes are located within the major histocompatibility complex (MHC) between HLA-B and HLA-DRB1 and are designated (in arbitrary order) by their CFB, C2, C4A, and C4B types. Pedigree analysis revealed long stretches (several megabases) of apparently fixed DNA within the MHC that we referred to as "extended haplotypes" (later as "CEHs"). About 10 to 12 common CEHs constitute at least 25 - 30% of MHC haplotypes among European Caucasian populations. These CEHs contain virtually all the most common markers of MHC-associated diseases. In the case of type 1 diabetes, we have proposed a purely genetic and epigenetic model (with a small number of Mendelian recessive disease genes) that explains all the puzzling features of the disease, including its rising incidence.

Haplotypic Associations and Differentiation of MHC Class II Polymorphic Alu Insertions at Five Loci With HLA-DRB1 Alleles in 12 Minority Ethnic Populations in China

The analysis of polymorphic variations in the human major histocompatibility complex (MHC) class II genomic region on the short-arm of chromosome 6 is a scientific enquiry to better understand the diversity in population structure and the effects of evolutionary processes such as recombination, mutation, genetic drift, demographic history, and natural selection. In order to investigate associations between the polymorphisms of HLA-DRB1 gene and recent Alu insertions (POALINs) in the HLA class II region, we genotyped HLA-DRB1 and five Alu loci (AluDPB2, AluDQA2, AluDQA1, AluDRB1, AluORF10), and determined their allele frequencies and haplotypic associations in 12 minority ethnic populations in China. There were 42 different HLA-DRB1 alleles for ethnic Chinese ranging from 12 alleles in the Jinuo to 28 in the Yugur with only DRB1^∗08:03, DRB1^∗09:01, DRB1^∗12:02, DRB1^∗14:01, DRB1^∗15:01, and DRB1^∗15:02 present in all ethnic groups. The POALINs varied in frequency between 0.279 and 0.514 for AluDPB2, 0 and 0.127 for AluDQA2, 0.777 and 0.995 for AluDQA1, 0.1 and 0.455 for AluDRB1 and 0.084 and 0.368 for AluORF10. By comparing the data of the five-loci POALIN in 13 Chinese ethnic populations (including Han-Yunnan published data) against Japanese and Caucasian published data, marked differences were observed between the populations at the allelic or haplotypic levels. Five POALIN loci were in significant linkage disequilibrium with HLA-DRB1 in different populations and AluDQA1 had the highest percentage association with most of the HLA-DRB1 alleles, whereas the nearby AluDRB1 indel was strongly haplotypic for only DRB1^∗01, DRB1^∗10, DRB1^∗15 and DRB1^∗16. There were 30 five-locus POALIN haplotypes inferred in all populations with H5 (no Alu insertions except for AluDQA1) and H21 (only AluDPB2 and AluDQA1 insertions) as the two predominant haplotypes. Neighbor joining trees and principal component analyses of the Alu and HLA-DRB1 polymorphisms showed that genetic diversity of these genomic markers is associated strongly with the population characteristics of language family, migration and sociality. This comparative study of HLA-DRB1 alleles and multilocus, lineage POALIN frequencies of Chinese ethnic populations confirmed that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic and evolutionary markers for the identification of allele and haplotype lineages and genetic variations within the same and/or different populations.

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.

High-Resolution HLA Typing of HLA-A, -B, -C, -DRB1, and -DQB1 in Kinh Vietnamese by Using Next-Generation Sequencing

Human leukocyte antigen (HLA) genotyping displays the particular characteristics of HLA alleles and haplotype frequencies in each population. Although it is considered the current gold standard for HLA typing, high-resolution sequence-based HLA typing is currently unavailable in Kinh Vietnamese populations. In this study, high-resolution sequence-based HLA typing (3-field) was performed using an amplicon-based next-generation sequencing platform to identify the HLA-A, -B, -C, -DRB1, and -DQB1 alleles of 101 unrelated healthy Kinh Vietnamese individuals from southern Vietnam. A total of 28 HLA-A, 41 HLA-B, 21 HLA-C, 26 HLA-DRB1, and 25 HLA-DQB1 alleles were identified. The most frequently occurring HLA alleles were A^∗11:01:01, B^∗15:02:01, C^∗07:02:01, DRB1^∗12:02:01, and DQB1^∗03:01:01. Haplotype calculation showed that A^∗29:01:01∼B^∗07:05:01, DRB1^∗12:02:01∼DQB1^∗3:01:01, A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01, A^∗33:03:01∼B^∗58:01:01∼DRB1^∗03:01:01, and A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01∼DRB1^∗10:01:01∼DQB1^∗05:01:01 were the most common haplotypes in the southern Kinh Vietnamese population. Allele distribution and haplotype analyses demonstrated that the Vietnamese population shares HLA features with South-East Asians but retains unique characteristics. Data from this study will be potentially applicable in medicine and anthropology.

MHC Genomics and Disease: Looking Back to Go Forward

Ancestral haplotypes are conserved but extremely polymorphic kilobase sequences, which have been faithfully inherited over at least hundreds of generations in spite of migration and admixture. They carry susceptibility and resistance to diverse diseases, including deficiencies of CYP21 hydroxylase (47.1) and complement components (18.1), as well as numerous autoimmune diseases (8.1). The haplotypes are detected by segregation within ethnic groups rather than by SNPs and GWAS. Susceptibility to some other diseases is carried by specific alleles shared by multiple ancestral haplotypes, e.g., ankylosing spondylitis and narcolepsy. The difference between these two types of association may explain the disappointment with many GWAS. Here we propose a pathway for combining the two different approaches. SNP typing is most useful after the conserved ancestral haplotypes have been defined by other methods.

A New Pedigree-Based SNP Haplotype Method for Genomic Polymorphism and Genetic Studies

Single nucleotide polymorphisms (SNPs) are usually the most frequent genomic variants. Directly pedigree-phased multi-SNP haplotypes provide a more accurate view of polymorphic population genomic structure than individual SNPs. The former are, therefore, more useful in genetic correlation with subject phenotype. We describe a new pedigree-based methodology for generating non-ambiguous SNP haplotypes for genetic study. SNP data for haplotype analysis were extracted from a larger Type 1 Diabetes Genetics Consortium SNP dataset based on minor allele frequency variation and redundancy, coverage rate (the frequency of phased haplotypes in which each SNP is defined) and genomic location. Redundant SNPs were eliminated, overall haplotype polymorphism was optimized and the number of undefined haplotypes was minimized. These edited SNP haplotypes from a region containing HLA-DRB1 (DR) and HLA-DQB1 (DQ) both correlated well with HLA-typed DR,DQ haplotypes and differentiated HLA-DR,DQ fragments shared by three pairs of previously identified megabase-length conserved extended haplotypes. In a pedigree-based genetic association assay for type 1 diabetes, edited SNP haplotypes and HLA-typed HLA-DR,DQ haplotypes from the same families generated essentially identical qualitative and quantitative results. Therefore, this edited SNP haplotype method is useful for both genomic polymorphic architecture and genetic association evaluation using SNP markers with diverse minor allele frequencies.

Common and well-documented HLA alleles of German stem cell donors by haplotype frequency estimation

We present a catalog of common and well-documented (CWD) alleles of the German population for the six HLA loci A, B, C, DRB1, DQB1, and DPB1. This study is based on a sample of over 5 million volunteer adult hematopoietic stem cell donors from the 26 German donor centers. To establish the catalog, allele and haplotype frequencies were estimated with a validated implementation of the expectation-maximization algorithm. CWD criteria similar to existing CWD catalogs were applied in order to be able to put our findings into the context of relevant existing references. Overall, 2155 HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1 alleles were identified as CWD in the German donor population representing about 20% of the HLA alleles at two-field resolution in the IPD-IMGT/HLA Database release v3.25.0 from July 2016 for these six loci. We found a substantial concordance of CWD alleles between the three catalogs and showed the contribution of the German donor population to the CWD alleles domain. In conclusion, the definition of CWD criteria that allow interoperability, scalability, and flexibility will be crucial for the development of a worldwide CWD catalog.

The MHC in the era of next-generation sequencing: Implications for bridging structure with function

The MHC continues to have the most disease-associations compared to other regions of the human genome, even in the genome-wide association study (GWAS) and single nucleotide polymorphism (SNP) era. Analysis of non-coding variation and their impact on the level of expression of HLA allotypes has shed new light on the potential mechanisms underlying HLA disease associations and alloreactivity in transplantation. Next-generation sequencing (NGS) technology has the capability of delineating the phase of variants in the HLA antigen-recognition site (ARS) with non-coding regulatory polymorphisms. These relationships are critical for understanding the qualitative and quantitative implications of HLA gene diversity. This article summarizes current understanding of non-coding region variation of HLA loci, the consequences of regulatory variation on HLA expression, the role for evolution in shaping lineage-specific expression, and the impact of HLA expression on disease susceptibility and transplantation outcomes. A role for phased sequencing methods for the MHC, and perspectives for future directions in basic and applied immunogenetic studies of the MHC are presented.

A stochastic epigenetic Mendelian oligogenic disease model for type 1 diabetes

The incidence of type 1 diabetes (T1D) and some other complex diseases is increasing. The cause has been attributed to an undefined changing environment. We examine the role of the environment (or any changing non-genetic mechanism) in causing the rising incidence, and find much evidence against it: 1) Dizygotic twin T1D concordance is the same as siblings of patients in general; 2) If the environment is responsible for both the discordance among identical twins of patients with T1D and its rising incidence, the twin concordance rate should be rising, but it is not; 3) Migrants from high-to low-incidence countries continue to have high-incidence children; 4) TID incidence among the offspring of two T1D parents is identical to the monozygotic twin rate. On the other hand, genetic association studies of T1D have revealed strong susceptibility in the major histocompatibility complex and many optional additive genes of small effect throughout the human genome increasing T1D risk. We have, from an analysis of previously published family studies, developed a stochastic epigenetic Mendelian oligogenic (SEMO) model consistent with published observations. The model posits a few required recessive causal genes with incomplete penetrance explaining virtually all of the puzzling features of T1D, including its rising incidence and the specific low T1D incidence rates among first-degree relatives of patients. Since historic selection against any causal gene could prevent T1D, we postulate that the rising incidence is because of increasing population mixing of parents from some previously isolated populations that had selected against different causal genes.

Haplotypes for Type, Degree, and Rate of Marbling in Cattle Are Syntenic with Human Muscular Dystrophy

Traditional analyses of a QTL on Bota 19 implicate a surfeit of candidates, but each is of marginal significance in explaining the deposition of healthy, low melting temperature fat within marbled muscle of Wagyu cattle. As an alternative approach, we have used genomic, multigenerational segregation to identify 14 conserved, ancestral 20 Mb haplotypes. These determine the degree and rate of marbling in Wagyu and other breeds of cattle. The melting temperature of intramuscular fat is highly heritable and traceable by haplotyping. Fortunately, for the production of healthy beef, some of these haplotypes are sufficiently penetrant to be expressed in heterozygous crossbreds, thereby allowing selection of sires which will improve the healthiness of beef produced under even harsh climatic conditions. The region of Bota 19 is syntenic to a region of Hosa 17 known to be important in muscle metabolism and in determining susceptibility to a form of human muscular dystrophy.

Pedigree-Defined Haplotypes and Their Applications to Genetic Studies

A haplotype is a string of nucleotides or alleles at nearby loci on one chromosome, usually inherited as a unit. Within the major histocompatibility complex (MHC) region on human chromosome 6p, independent population studies of multiple families have identified conserved extended haplotypes (CEHs) that segregate as long stretches (≥1 megabase) of essentially identical DNA sequence at relatively high (≥0.5 %) population frequency ("genetic fixity"). CEHs were first identified through segregation analysis in the early 1980s. In European Caucasian populations, the most frequent 30 CEHs account for at least one-third of all MHC haplotypes. These CEHs provide all of the known individual MHC susceptibility and protective genetic markers within those populations for several complex genetic diseases. Haplotypes are rigorously determined directly by sequencing single chromosomes or by Mendelian segregation analysis using families with informative genotypes. Four parental haplotypes are assigned unambiguously using genotypes from the two parents and from two of their haploidentical (to each other) children. However, the most common current technique to phase haplotypes is probabilistic statistical imputation, using unrelated subjects. Such probabilistic techniques have failed to detect CEHs and are thus of questionable value in identifying long-range haplotype structure and, consequently, genetic structure-function relationships. Finally, with haplotypes rigorously defined, association studies can determine frequencies of alleles among unrelated patient haplotypes vs. those among only unaffected family members (i.e., control alleles/haplotypes). Such studies reduce, as much as possible, the confounding effects of population stratification common to all genetic studies.

Ensembl core software resources: storage and programmatic access for DNA sequence and genome annotation

The Ensembl software resources are a stable infrastructure to store, access and manipulate genome assemblies and their functional annotations. The Ensembl 'Core' database and Application Programming Interface (API) was our first major piece of software infrastructure and remains at the centre of all of our genome resources. Since its initial design more than fifteen years ago, the number of publicly available genomic, transcriptomic and proteomic datasets has grown enormously, accelerated by continuous advances in DNA-sequencing technology. Initially intended to provide annotation for the reference human genome, we have extended our framework to support the genomes of all species as well as richer assembly models. Cross-referenced links to other informatics resources facilitate searching our database with a variety of popular identifiers such as UniProt and RefSeq. Our comprehensive and robust framework storing a large diversity of genome annotations in one location serves as a platform for other groups to generate and maintain their own tailored annotation. We welcome reuse and contributions: our databases and APIs are publicly available, all of our source code is released with a permissive Apache v2.0 licence at http://github.com/Ensembl and we have an active developer mailing list ( http://www.ensembl.org/info/about/contact/index.html ).

Database URL

http://www.ensembl.org.

Genome-wide autozygosity is associated with lower general cognitive ability

Inbreeding depression refers to lower fitness among offspring of genetic relatives. This reduced fitness is caused by the inheritance of two identical chromosomal segments (autozygosity) across the genome, which may expose the effects of (partially) recessive deleterious mutations. Even among outbred populations, autozygosity can occur to varying degrees due to cryptic relatedness between parents. Using dense genome-wide single-nucleotide polymorphism (SNP) data, we examined the degree to which autozygosity associated with measured cognitive ability in an unselected sample of 4854 participants of European ancestry. We used runs of homozygosity-multiple homozygous SNPs in a row-to estimate autozygous tracts across the genome. We found that increased levels of autozygosity predicted lower general cognitive ability, and estimate a drop of 0.6 s.d. among the offspring of first cousins (P=0.003-0.02 depending on the model). This effect came predominantly from long and rare autozygous tracts, which theory predicts as more likely to be deleterious than short and common tracts. Association mapping of autozygous tracts did not reveal any specific regions that were predictive beyond chance after correcting for multiple testing genome wide. The observed effect size is consistent with studies of cognitive decline among offspring of known consanguineous relationships. These findings suggest a role for multiple recessive or partially recessive alleles in general cognitive ability, and that alleles decreasing general cognitive ability have been selected against over evolutionary time.

HLA-A*33-DR3 and A*33-DR9 haplotypes enhance the risk of type 1 diabetes in Han Chinese

Aims/Introduction

To investigate the typing for human leukocyte antigen (HLA) class I in Chinese patients with type 1 diabetes as a complement screening for HLA class II.

Materials and Methods

A total of 212 type 1 diabetic patients and 200 healthy controls were enrolled. The genetic polymorphisms of HLA class I and II were examined with a high‐resolution polymerase chain reaction sequence‐based typing method.

Results

The haplotype, A*33:03‐B*58:01‐C*03:02(A33), was associated with type 1 diabetes (P = 1.0 × 10⁻⁴, odds ratio 3.2 [1.738–5.843]). The A33‐DR3 and A33‐DR9 haplotypes significantly enhanced the risk of type 1 diabetes (A33‐DR3, odds ratio 5.1 [2.40–10.78], P = 4.0 × 10⁻⁶; A33‐DR9, odds ratio 13.0 [1.69–100.32], P = 0.004). In type 1 diabetic patients, compared with A33‐DR3‐negative carriers, A33‐DR3‐positive carriers had significantly lower percentages of CD3⁺CD4⁺ T cells (42.5 ± 7.72 vs 37.0 ± 8.35%, P = 0.023), higher percentages of CD3⁺CD8⁺ T cells (27.4 ± 7.09 vs 32.8 ± 5.98%, P = 0.005) and T‐cell receptor α/β T cells (70.0 ± 7.00 vs 73.6 ± 6.25%, P = 0.031), and lower CD4/CD8 ratios (1.71 ± 0.75 vs 1.16 ± 0.35, P = 0.003).

Conclusions

It is the first time that the haplotypes A33‐DR3 and A33‐DR9 were found with an enhanced predisposition to type 1 diabetes in Han Chinese. A33‐DR3 was associated with a reduction in the helper‐to‐cytotoxic cell ratio and preferential increase of T‐cell receptor α/β T cell. The typing for HLA class I and its immunogenetic effects are important for more accurate HLA class II haplotype risk prediction and etiology research in type 1 diabetic patients.

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.

Soma-to-germline feedback is implied by the extreme polymorphism at IGHV relative to MHC: The manifest polymorphism of the MHC appears greatly exceeded at Immunoglobulin loci, suggesting antigen-selected somatic V mutants penetrate Weismann's Barrier

Soma-to-germline feedback is forbidden under the neo-Darwinian paradigm. Nevertheless, there is a growing realization it occurs frequently in immunoglobulin (Ig) variable (V) region genes. This is a surprising development. It arises from a most unlikely source in light of the exposure of co-author EJS to the haplotype data of RL Dawkins and others on the polymorphism of the Major Histocompatibility Complex, which is generally assumed to be the most polymorphic region in the genome (spanning ∼4 Mb). The comparison between the magnitude of MHC polymorphism with estimates for the human heavy chain immunoglobulin V locus (spanning ∼1 Mb), suggests IGHV could be many orders of magnitude more polymorphic than the MHC. This conclusion needs airing in the literature as it implies generational churn and soma-to-germline gene feedback. Pedigree-based experimental strategies to resolve the IGHV issue are outlined.

No Evidence of Human Leukocyte Antigen Gene Association With Rheumatic Fever Among Children in Samoa

Human leukocyte antigens (HLAs) have been implicated in rheumatic fever pathogenesis. This pilot whole genome association study compares genotypes of Samoan children with rheumatic fever to unaffected siblings and unrelated healthy controls. No risk-related genotypes were associated with HLA genes. Thirteen Regions of Interest were identified as candidates for further study.

HLA-B allele and haplotype diversity among Thai patients identified by PCR-SSOP: evidence for high risk of drug-induced hypersensitivity

Background: There are 3 classes of HLA molecules; HLA class I, II and III, of which different classes have different functions. HLA-B gene which belongs to HLA class I play an important role predicting drug hypersensitivity.

Materials and Methods: Nine hundred and eighty-six Thai subjects who registered at a pharmacogenomics laboratory were determined for HLA-B genotype using a two-stage sequence-specific oligonucleotide probe system (PCR-SSOP).

Results: In this study, HLA-B alleles did not deviate from Hardy-Weinberg equilibrium (P > 0.05). The most common HLA-B alleles observed in this population were HLA-B^*46:01 (11.51%), HLA-B^*58:01 (8.62%), HLA-B^*40:01 (8.22%), HLA-B^*15:02 (8.16%) and HLA-B^*13:01 (6.95%). This finding revealed that HLA-B allele frequency in the Thai population was consistent with the Chinese population (p > 0.05), however, differed from the Malaysian population (p < 0.05). The top five HLA-B genotypes were HLA-B^*40:01/46:01 (2.13%), HLA-B^*46:01/46:01 (2.03%), HLA-B^*40:01/58:01 (2.03%), HLA-B^*46:01/58:01 (1.93%) and HLA-B^*15:02/46:01 (1.83%). This study found that 15.92% of Thai subjects carry HLA-B^*15:02, which has been associated with carbamazepine-induced severe cutaneous adverse drug reactions (SCARs). Moreover, 16.33% of Thai subjects carry the HLA-B^*58:01 allele, which has been associated with allopurinol-induced SCARs.

Conclusion: This study demonstrates a high diversity of HLA-B polymorphisms in this Thai population. The high frequency of HLA-B pharmacogenomic markers in the population emphasizes the importance of such screening to predict/avoid drug hypersensitivity.

Human herpesvirus 6B immediate-early I protein contains functional HLA-A*02, HLA-A*03, and HLA-B*07 class I restricted CD8(+) T-cell epitopes

Human herpesvirus 6B (HHV-6B) is a ubiquitous pathogen with frequent reactivation observed in immunocompromised patients such as BM transplant (BMT) recipients. Adoptive immunotherapy is a promising therapeutic avenue for the treatment of opportunistic infections, including herpesviruses. While T-cell immunotherapy can successfully control CMV and EBV reactivations in BMT recipients, such therapy is not available for HHV-6 infections, in part due to a lack of identified protective CD8(+) T-cell epitopes. Our goal was to identify CD8(+) T-cell viral epitopes derived from the HHV-6B immediate-early protein I and presented by common human leukocyte Ag (HLA) class I alleles including HLA-A*02, HLA-A*03, and HLA-B*07. These epitopes were functionally tested for their ability to induce CD8(+) T-cell expansion and kill HHV-6-infected autologous cells. Cross-reactivity of specific HHV-6B-expanded T cells against HHV-6A-infected cells was also confirmed for a conserved epitope presented by HLA-A*02 molecule. Our findings will help push forward the field of adoptive immunotherapy for the treatment and/or the prevention of HHV-6 reactivation in BMT recipients.

Dominant sequences of human major histocompatibility complex conserved extended haplotypes from HLA-DQA2 to DAXX

We resequenced and phased 27 kb of DNA within 580 kb of the MHC class II region in 158 population chromosomes, most of which were conserved extended haplotypes (CEHs) of European descent or contained their centromeric fragments. We determined the single nucleotide polymorphism and deletion-insertion polymorphism alleles of the dominant sequences from HLA-DQA2 to DAXX for these CEHs. Nine of 13 CEHs remained sufficiently intact to possess a dominant sequence extending at least to DAXX, 230 kb centromeric to HLA-DPB1. We identified the regions centromeric to HLA-DQB1 within which single instances of eight “common” European MHC haplotypes previously sequenced by the MHC Haplotype Project (MHP) were representative of those dominant CEH sequences. Only two MHP haplotypes had a dominant CEH sequence throughout the centromeric and extended class II region and one MHP haplotype did not represent a known European CEH anywhere in the region. We identified the centromeric recombination transition points of other MHP sequences from CEH representation to non-representation. Several CEH pairs or groups shared sequence identity in small blocks but had significantly different (although still conserved for each separate CEH) sequences in surrounding regions. These patterns partly explain strong calculated linkage disequilibrium over only short (tens to hundreds of kilobases) distances in the context of a finite number of observed megabase-length CEHs comprising half a population's haplotypes. Our results provide a clearer picture of European CEH class II allelic structure and population haplotype architecture, improved regional CEH markers, and raise questions concerning regional recombination hotspots.

The human major histocompatibility complex (MHC) is a gene-dense region highly enriched in immune response genes. MHC genetic variation is among the highest in the human genome and is associated with both tissue transplant compatibility and many genetic disorders. Long-range (1–3 Mb) MHC haplotypes of essentially identical DNA sequence at relatively high (≥0.5%) population frequency (“genetic fixity”), called conserved extended haplotypes (CEHs), comprise roughly half of all European population haplotypes. We sequenced an aggregate of 27 kb over 580 kb in the MHC class II region from HLA-DQA2 to DAXX in 158 European haplotypes to quantify the breakdown of this genetic fixity in the centromeric portion of the MHC and to determine the representative nature within that region of eight previously fully or nearly fully sequenced “common” European haplotypes. We identified the dominant sequences of 13 European CEHs and determined where the “common” sequences did (or did not) represent related CEHs. We found patterns of shared sequence identity among different CEHs surrounded by fixed (for each CEH) but differing sequence. Our direct observational results for population haplotypes explain the mutual occurrence of CEHs and short (5–200 kb) blocks of fixed sequence detected by the statistical measure of linkage disequilibrium.

JC polyomavirus infection is strongly controlled by human leucocyte antigen class II variants

JC polyomavirus (JCV) carriers with a compromised immune system, such as in HIV, or subjects on immune-modulating therapies, such as anti VLA-4 therapy may develop progressive multifocal leukoencephalopathy (PML) which is a lytic infection of oligodendrocytes in the brain. Serum antibodies to JCV mark infection occur only in 50–60% of infected individuals, and high JCV-antibody titers seem to increase the risk of developing PML. We here investigated the role of human leukocyte antigen (HLA), instrumental in immune defense in JCV antibody response. Anti-JCV antibody status, as a surrogate for JCV infection, were compared to HLA class I and II alleles in 1621 Scandinavian persons with MS and 1064 population-based Swedish controls and associations were replicated in 718 German persons with MS. HLA-alleles were determined by SNP imputation, sequence specific (SSP) kits and a reverse PCR sequence-specific oligonucleotide (PCR-SSO) method. An initial GWAS screen displayed a strong HLA class II region signal. The HLA-DRB1*15 haplotype was strongly negatively associated to JCV sero-status in Scandinavian MS cases (OR = 0.42, p = 7×10⁻¹⁵) and controls (OR = 0.53, p = 2×10⁻⁵). In contrast, the DQB1*06:03 haplotype was positively associated with JCV sero-status, in Scandinavian MS cases (OR = 1.63, p = 0.006), and controls (OR = 2.69, p = 1×10⁻⁵). The German dataset confirmed these findings (OR = 0.54, p = 1×10⁻⁴ and OR = 1.58, p = 0.03 respectively for these haplotypes). HLA class II restricted immune responses, and hence CD4+ T cell immunity is pivotal for JCV infection control. Alleles within the HLA-DR1*15 haplotype are associated with a protective effect on JCV infection. Alleles within the DQB1*06:03 haplotype show an opposite association. These associations between JC virus antibody response and human leucocyte antigens supports the notion that CD4+ T cells are crucial in the immune defence to JCV and lays the ground for risk stratification for PML and development of therapy and prevention.

JC virus infection can lead to progressive multifocal leukoencephalopathy in individuals with a compromised immune system, such as during HIV infections or when treated with immunosuppressive or immunomodulating therapies. Progressive multifocal leukoencephalopathy is a rare but potentially fatal disease characterized by progressive damage of the brain white matter at multiple locations. It is therefore of importance to understand the host genetic control of response to JC virus in order to identify patients that can be treated with immunomodulating therapies, common treatments for autoimmune diseases, without increased risk for progressive multifocal leukoencephalopathy. This may also lead to development of preventative or curative anti-JC virus therapies. We here identify genetic variants being associated with JC virus antibody development; a negative association with the human leucocyte antigen DRB1*15-DQA1*01:02-DQB1*06:02 haplotype and a positive association with the DRB1*13-DQA1*01:03-DQB1*06:03 haplotype among controls and patients with multiple sclerosis from Scandinavia. We confirmed the associations in patients with multiple sclerosis from Germany. These associations between JC virus antibody response and human leucocyte antigens imply that CD4+ T cells are crucial in the immune defence and lay the ground for development of therapy and prevention.

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.

Human leukocyte antigen alleles, genotypes and haplotypes frequencies in renal transplant donors and recipients from West Central India

BACKGROUND:

Human leukocyte antigen (HLA) is comprised of a highly polymorphic set of genes which determines the histocompatibility of organ transplantation. The present study was undertaken to identify HLA class I and class II allele, genotype and haplotype frequencies in renal transplant recipients and donors from West Central India.

MATERIALS AND METHODS:

HLA typing was carried out using Polymerase Chain Reaction-Sequence Specific Primer in 552 live related and unrelated renal transplant recipients and donors.

RESULTS:

The most frequent HLA class I and class II alleles and their frequencies in recipients were HLA-AFNx0101 (0.1685) and AFNx0102 (0.1649), HLA-BFNx0135 (0.1322), and HLA-DR beta 1 (DRB 1)FNx0115 (0.2192), whereas in donors, these were HLA-AFNx0102 (0.1848) and AFNx0101 (0.1667), HLA-BFNx0135 (0.1359), and HLA-DRB1FNx0115 (0.2409). The two-locus haplotype statistical analysis revealed HLA-AFNx0102-B61 as the most common haplotype with the frequency of 0.0487 and 0.0510 in recipients and donors, respectively. Further, among the three locus haplotypes HLA-AFNx0133-BFNx0144-DRB1FNx0107 and HLA-AFNx0102-BFNx0161-DRB1FNx0115 were the most common haplotypes with frequencies 0.0362 and 0.0326, respectively in recipients and 0.0236 and 0.0323, respectively in donors. Genotype frequency revealed a high prevalence of genotype HLA-AFNx0102/AFNx0124 in recipients (0.058) compared to donors (0.0109) whereas low prevalence of HLA-AFNx0101/AFNx0102 in recipients (0.0435) than in donors (0.0797). The phylogenetic and principal component analysis of HLA allele and haplotype frequency distribution revealed genetic similarities of various ethnic groups. Further, case control analysis provides preliminary evidence of association of HLA-A genotype (P < 0.05) with renal failure.

CONCLUSION:

This study will be helpful in suitable donor search besides providing valuable information for population genetics and HLA disease association analysis.

SNP association mapping across the extended major histocompatibility complex and risk of B-cell precursor acute lymphoblastic leukemia in children

The extended major histocompatibility complex (xMHC) is the most gene-dense region of the genome and harbors a disproportionately large number of genes involved in immune function. The postulated role of infection in the causation of childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) suggests that the xMHC may make an important contribution to the risk of this disease. We conducted association mapping across an approximately 4 megabase region of the xMHC using a validated panel of single nucleotide polymorphisms (SNPs) in childhood BCP-ALL cases (n=567) enrolled in the Northern California Childhood Leukemia Study (NCCLS) compared with population controls (n=892). Logistic regression analyses of 1,145 SNPs, adjusted for age, sex, and Hispanic ethnicity indicated potential associations between several SNPs and childhood BCP-ALL. After accounting for multiple comparisons, one of these included a statistically significant increased risk associated with rs9296068 (OR=1.40, 95% CI=1.19-1.66, corrected p=0.036), located in proximity to HLA-DOA. Sliding window haplotype analysis identified an additional locus located in the extended class I region in proximity to TRIM27 tagged by a haplotype comprising rs1237485, rs3118361, and rs2032502 (corrected global p=0.046). Our findings suggest that susceptibility to childhood BCP-ALL is influenced by genetic variation within the xMHC and indicate at least two important regions for future evaluation.

HLA polymorphisms influence the development of skin rash arising from treatment with EGF receptor inhibitors

AIM

Development of a skin rash under treatment with EGF receptor (EGFR) inhibitors (EGFRIs) has been linked to a favorable prognosis in some studies, suggesting a possible immunological role for EGFRIs in addition to direct antagonistic downstream effects. The present study aimed to investigate whether particular HLA polymorphisms found in cancer patients treated with EGFRIs are associated with the development of skin rash and overall survival rates.

PATIENTS & METHODS

HLA typing was performed on 105 cancer patients and the development of skin rash was rated during the first 4 weeks of therapy with EGFRIs.

RESULTS

A significantly lower incidence of skin rash was found in patients carrying the HLA-A*02:01 or HLA-A*03:01 alleles (hazard ratio: 0.277; 95% CI: 0.121-0.634; p = 0.002 and hazard ratio: 0.292; 95% CI: 0.113-0.752; p = 0.011, respectively); however, no association with worse survival was seen.

CONCLUSION

The chances of developing a skin rash in patients treated with EGFRIs may be lower in patients that carry the HLA-A*02:01 or HLA-A*03:01 alleles, while the antitumor efficacy of EGFRIs does not seem to be significantly impaired in these patients.

Genomic evaluation of HLA-DR3+ haplotypes associated with type 1 diabetes

We have defined three sets of HLA-DR3(+) haplotypes that provide maximum risk of type 1 disease development in Indians: (1) a diverse array of B8-DR3 haplotypes, (2) A33-B58-DR3 haplotype, and (3) A2-B50-DR3 occurring most predominantly in this population. Further analysis has revealed extensive diversity in B8-DR3 haplotypes, particularly at the HLA-A locus, in contrast to the single fixed HLA-A1-B8-DR3 haplotype (generally referred to as AH8.1) reported in Caucasians. However, the classical AH8.1 haplotype was rare and differed from the Caucasian counterpart at multiple loci. In our study, HLA-A26-B8-DR3 (AH8.2) was the most common B8-DR3 haplotype constituting >50% of the total B8-DR3 haplotypes. Further, A2-B8-DR3 contributed the maximum risk (RR = 48.7) of type 1 diabetes, followed by A2-B50-DR3 (RR = 9.4), A33-B58-DR3 (RR = 6.6), A24-B8-DR3 (RR = 4.5), and A26-B8-DR3 (RR = 4.2). Despite several differences, the disease-associated haplotypes in Indian and Caucasian populations share a frozen DR3-DQ2 block, suggesting a common ancestor from which multiple haplotypes evolved independently.

Evolutionary analysis of classical HLA class I and II genes suggests that recent positive selection acted on DPB1*04:01 in Japanese population

The human leukocyte antigen (HLA) genes exhibit the highest degree of polymorphism in the human genome. This high degree of variation at classical HLA class I and class II loci has been maintained by balancing selection for a long evolutionary time. However, little is known about recent positive selection acting on specific HLA alleles in a local population. To detect the signature of recent positive selection, we genotyped six HLA loci, HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQB1, and HLA-DPB1 in 418 Japanese subjects, and then assessed the haplotype homozygosity (HH) of each HLA allele. There were 120 HLA alleles across the six loci. Among the 80 HLA alleles with frequencies of more than 1%, DPB1*04∶01, which had a frequency of 6.1%, showed exceptionally high HH (0.53). This finding raises the possibility that recent positive selection has acted on DPB1*04∶01. The DPB1*04∶01 allele, which was present in the most common 6-locus HLA haplotype (4.4%), A*33∶03-C*14∶03-B*44∶03-DRB1*13∶02-DQB1*06∶04-DPB1*04∶01, seems to have flowed from the Korean peninsula to the Japanese archipelago in the Yayoi period. A stochastic simulation approach indicated that the strong linkage disequilibrium between DQB1*06∶04 and DPB1*04∶01 observed in Japanese cannot be explained without positive selection favoring DPB1*04∶01. The selection coefficient of DPB1*04∶01 was estimated as 0.041 (95% credible interval 0.021–0.077). Our results suggest that DPB1*04∶01 has recently undergone strong positive selection in Japanese population.

Predicting the onset of Addison's disease: ACTH, renin, cortisol and 21-hydroxylase autoantibodies

CONTEXT

Autoantibodies to 21-hydroxylase (21OH-AA) precede onset of autoimmune Addison's disease (AD). Progression to AD can take months to years, and early detection of metabolic decompensation may prevent morbidity and mortality.

OBJECTIVE

To define optimal methods of predicting progression to overt AD (defined by subnormal peak cortisol response to Cosyntropin) in 21OH-AA+ individuals.

DESIGN, SETTING AND PARTICIPANTS

Individuals were screened for 21OH-AA at the Barbara Davis Center from 1993 to 2011. Subjects positive for 21OH-AA (n = 87) were tested, and the majority prospectively followed for the development of Addison's disease, including seven diagnosed with AD upon 21OH-AA discovery (discovered), seven who progressed to AD (progressors) and 73 nonprogressors.

MAIN OUTCOME MEASURED

Plasma renin activity (PRA), ACTH, baseline cortisol, peak cortisol and 21OH-AA were measured at various time points relative to diagnosis of AD or last AD-free follow-up.

RESULTS

Compared with nonprogressors, in the time period 2 months-2 years prior to the onset of AD, progressors were significantly more likely to have elevated ACTH (11-22 pM, P < 1E-4), with no significant differences in mean PRA (P = 0·07) or baseline cortisol (P = 0·08), and significant but less distinct differences seen with 21OH-AA levels (P < 1E-4) and poststimulation cortisol levels (P = 6E-3).

CONCLUSION

Moderately elevated ACTH is a more useful early indicator of impending AD than 21OH-AA, PRA or peak cortisol, in the 2 months-2 years preceding the onset of AD.

Runs of homozygosity implicate autozygosity as a schizophrenia risk factor

Autozygosity occurs when two chromosomal segments that are identical from a common ancestor are inherited from each parent. This occurs at high rates in the offspring of mates who are closely related (inbreeding), but also occurs at lower levels among the offspring of distantly related mates. Here, we use runs of homozygosity in genome-wide SNP data to estimate the proportion of the autosome that exists in autozygous tracts in 9,388 cases with schizophrenia and 12,456 controls. We estimate that the odds of schizophrenia increase by ~17% for every 1% increase in genome-wide autozygosity. This association is not due to one or a few regions, but results from many autozygous segments spread throughout the genome, and is consistent with a role for multiple recessive or partially recessive alleles in the etiology of schizophrenia. Such a bias towards recessivity suggests that alleles that increase the risk of schizophrenia have been selected against over evolutionary time.

Congruence as a measurement of extended haplotype structure across the genome

Background

Historically, extended haplotypes have been defined using only a few data points, such as alleles for several HLA genes in the MHC. High-density SNP data, and the increasing affordability of whole genome SNP typing, creates the opportunity to define higher resolution extended haplotypes. This drives the need for new tools that support quantification and visualization of extended haplotypes as defined by as many as 2000 SNPs. Confronted with high-density SNP data across the major histocompatibility complex (MHC) for 2,300 complete families, compiled by the Type 1 Diabetes Genetics Consortium (T1DGC), we developed software for studying extended haplotypes.

Methods

The software, called ExHap (Extended Haplotype), uses a similarity measurement we term congruence to identify and quantify long-range allele identity. Using ExHap, we analyzed congruence in both the T1DGC data and family-phased data from the International HapMap Project.

Results

Congruent chromosomes from the T1DGC data have between 96.5% and 99.9% allele identity over 1,818 SNPs spanning 2.64 megabases of the MHC (HLA-DRB1 to HLA-A). Thirty-three of 132 DQ-DR-B-A defined haplotype groups have > 50% congruent chromosomes in this region. For example, 92% of chromosomes within the DR3-B8-A1 haplotype are congruent from HLA-DRB1 to HLA-A (99.8% allele identity). We also applied ExHap to all 22 autosomes for both CEU and YRI cohorts from the International HapMap Project, identifying multiple candidate extended haplotypes.

Conclusions

Long-range congruence is not unique to the MHC region. Patterns of allele identity on phased chromosomes provide a simple, straightforward approach to visually and quantitatively inspect complex long-range structural patterns in the genome. Such patterns aid the biologist in appreciating genetic similarities and differences across cohorts, and can lead to hypothesis generation for subsequent studies.

The past, present, and future of genetic associations in type 1 diabetes

Type 1 diabetes mellitus (T1DM) is an autoimmune disease affecting approximately one in 300 individuals in the United States. The majority of genetic research to date has focused on the heritability that predisposes to islet autoimmunity and T1DM. The evidence so far points to T1DM being a polygenic, common, complex disease with major susceptibility lying in the major histocompatibility complex (MHC) on chromosome 6 with other smaller effects seen in loci outside of the MHC. With recent advances in technology, novel means of exploring the human genome have given way to new information in the development of T1DM. The newest technologies, namely high-throughput polymorphism typing and sequencing, have led to a paradigm shift in studying common diseases such as T1DM. In this review we highlight the advances in genetic associations in T1DM in the last several decades and how they have led to a better understanding of T1DM pathogenesis.

Almost total protection from age-related macular degeneration by haplotypes of the Regulators of Complement Activation

Age-related macular degeneration (AMD) is the leading cause of blindness in developed countries. It has been proposed that the polymorphism encoding Y402H (T1277C) in the complement factor H gene (CFH) is one of the main determinants of disease. We genotyped the polymorphism at a number of loci in the region encompassing the Regulators of Complement Activation (RCA) on chromosome 1, including T1277C SNP, in 187 patients and 146 controls. Haplotypes have been classified as protective (P) or susceptible (S) with respect to AMD. This included the identification of an S haplotype with a T at 1277. The results show that no single locus should be assumed to be directly responsible for AMD, but rather argue for the existence of RCA haplotypes, which can be assigned meaningful predictive values for AMD. We conclude that the critical sequences are within a region 450 kb centromeric to 128 kb telomeric of CFH.

The HLA-B 3906 allele imparts a high risk of diabetes only on specific HLA-DR/DQ haplotypes

AIMS/HYPOTHESIS

We investigated the risk associated with HLA-B*39 alleles in the context of specific HLA-DR/DQ haplotypes.

METHODS

We studied a readily available dataset from the Type 1 Diabetes Genetics Consortium that consists of 2,300 affected sibling pair families genotyped for both HLA alleles and 2,837 single nucleotide polymorphisms across the major histocompatibility complex region.

RESULTS

The B*3906 allele significantly enhanced the risk of type 1 diabetes when present on specific HLA-DR/DQ haplotypes (DRB1 0801-DQB1 0402: p = 1.6 × 10(-6), OR 25.4; DRB1 0101-DQB1 0501: p = 4.9 × 10(-5), OR 10.3) but did not enhance the risk of DRB1 0401-DQB1 0302 haplotypes. In addition, the B 3901 allele enhanced risk on the DRB1 1601-DQB1 0502 haplotype (p = 3.7 × 10(-3), OR 7.2).

CONCLUSIONS/INTERPRETATION

These associations indicate that the B 39 alleles significantly increase risk when present on specific HLA-DR/DQ haplotypes, and HLA-B typing in concert with specific HLA-DR/DQ genotypes should facilitate genetic prediction of type 1 diabetes, particularly in a research setting.

Dominant suppression of Addison's disease associated with HLA-B15

CONTEXT

Autoimmune Addison's disease (AD) is the major cause of primary adrenal failure in developed nations. Autoantibodies to 21-hydroxylase (21OH-AA) are associated with increased risk of progression to AD. Highest genetic risk is associated with the Major Histocompatibility region (MHC), specifically human leukocyte antigen (HLA)-DR3 haplotypes (containing HLA-B8) and HLA-DR4.

OBJECTIVE

The objective of the study was the further characterization of AD risk associated with MHC alleles.

DESIGN, SETTING, AND PARTICIPANTS

MHC genotypes were determined for HLA-DRB1, DQA1, DQB1, MICA, HLA-B, and HLA-A in 168 total individuals with 21OH-AA (85 with AD at referral and 83 with positive 21OH-AA but without AD at referral).

MAIN OUTCOME MEASURE(S)

Genotype was evaluated in 21OH-AA-positive individuals. Outcomes were compared with general population controls and type 1 diabetes patients.

RESULTS

In HLA-DR4+ individuals, HLA-B15 was found in only one of 55 (2%) with AD vs. 24 of 63 (40%) 21OH-AA-positive nonprogressors (P = 2 × 10(-7)) and 518 of 1558 (33%) HLA-DR4 patients with type 1 diabetes (P = 1 × 10(-8)). On prospective follow-up, none of the HLA-B15-positive, 21-hydroxylase-positive individuals progressed to AD vs. 25% non-HLA-B15 autoantibody-positive individuals by life table analysis (P = 0.03). Single nucleotide polymorphism analysis revealed the HLA-DR/DQ region associated with risk and HLA-B15 were separated by multiple intervening single-nucleotide polymorphism haplotypes.

CONCLUSIONS

HLA-B15 is not associated with protection from 21OH-AA formation but is associated with protection from progression to AD in 21OH-AA-positive individuals. To our knowledge, this is one of the most dramatic examples of genetic disease suppression in individuals who already have developed autoantibodies and of novel dominant suppression of an autoimmune disease by a class I HLA allele.

Genomic evolution in domestic cattle: ancestral haplotypes and healthy beef

We have identified numerous Ancestral Haplotypes encoding a 14-Mb region of Bota C19. Three are frequent in Simmental, Angus and Wagyu and have been conserved since common progenitor populations. Others are more relevant to the differences between these 3 breeds including fat content and distribution in muscle. SREBF1 and Growth Hormone, which have been implicated in the production of healthy beef, are included within these haplotypes. However, we conclude that alleles at these 2 loci are less important than other sequences within the haplotypes. Identification of breeds and hybrids is improved by using haplotypes rather than individual alleles.

Conditional meta-analysis stratifying on detailed HLA genotypes identifies a novel type 1 diabetes locus around TCF19 in the MHC

The human leukocyte antigen (HLA) class II genes HLA-DRB1, -DQA1 and -DQB1 are the strongest genetic factors for type 1 diabetes (T1D). Additional loci in the major histocompatibility complex (MHC) are difficult to identify due to the region's high gene density and complex linkage disequilibrium (LD). To facilitate the association analysis, two novel algorithms were implemented in this study: one for phasing the multi-allelic HLA genotypes in trio families, and one for partitioning the HLA strata in conditional testing. Screening and replication were performed on two large and independent datasets: the Wellcome Trust Case-Control Consortium (WTCCC) dataset of 2,000 cases and 1,504 controls, and the T1D Genetics Consortium (T1DGC) dataset of 2,300 nuclear families. After imputation, the two datasets have 1,941 common SNPs in the MHC, of which 22 were successfully tested and replicated based on the statistical testing stratifying on the detailed DRB1 and DQB1 genotypes. Further conditional tests using the combined dataset confirmed eight novel SNP associations around 31.3 Mb on chromosome 6 (rs3094663, p = 1.66 × 10(-11) and rs2523619, p = 2.77 × 10(-10) conditional on the DR/DQ genotypes). A subsequent LD analysis established TCF19, POU5F1, CCHCR1 and PSORS1C1 as potential causal genes for the observed association.

Haplotype analysis discriminates genetic risk for DR3-associated endocrine autoimmunity and helps define extreme risk for Addison's disease

CONTEXT

Multiple autoimmune disorders (e.g. Addison's disease, type 1 diabetes, celiac disease) are associated with HLA-DR3, but it is likely that alleles of additional genes in linkage disequilibrium with HLA-DRB1 contribute to disease.

OBJECTIVE

The objective of the study was to characterize major histocompatability complex (MHC) haplotypes conferring extreme risk for autoimmune Addison's disease (AD).

DESIGN, SETTING, AND PARTICIPANTS

Eighty-six 21-hydroxylase autoantibody-positive, nonautoimmune polyendocrine syndrome type 1, Caucasian individuals collected from 1992 to 2009 with clinical AD from 68 families (12 multiplex and 56 simplex) were genotyped for HLA-DRB1, HLA-DQB1, MICA, HLA-B, and HLA-A as well as high density MHC single-nucleotide polymorphism (SNP) analysis for 34.

MAIN OUTCOME MEASURES

AD and genotype were measured.

RESULT

Ninety-seven percent of the multiplex individuals had both HLA-DR3 and HLA-B8 vs. 60% of simplex AD patients (P = 9.72 × 10(-4)) and 13% of general population controls (P = 3.00 × 10(-19)). The genotype DR3/DR4 with B8 was present in 85% of AD multiplex patients, 24% of simplex patients, and 1.5% of control individuals (P = 4.92 × 10(-191)). The DR3-B8 haplotype of AD patients had HLA-A1 less often (47%) than controls (81%, P = 7.00 × 10(-5)) and type 1 diabetes patients (73%, P = 1.93 × 10(-3)). Analysis of 1228 SNPs across the MHC for individuals with AD revealed a shorter conserved haplotype (3.8) with the loss of the extended conserved 3.8.1 haplotype approximately halfway between HLA-B and HLA-A.

CONCLUSION

Extreme risk for AD, especially in multiplex families, is associated with haplotypic DR3 variants, in particular a portion (3.8) but not all of the conserved 3.8.1 haplotype.