Julie A. Lane

HLA Class I and Genetic Susceptibility to Type 1 Diabetes Results From the Type 1 Diabetes Genetics Consortium

OBJECTIVE We report here genotyping data and type 1 diabetes association analyses for HLA class I loci (A, B, and C) on 1,753 multiplex pedigrees from the Type 1 Diabetes Genetics Consortium (T1DGC), a large international collaborative study. RESEARCH DESIGN AND METHODS Complete eight-locus HLA genotyping data were generated. Expected patient class I (HLA-A, -B, and -C) allele frequencies were calculated, based on linkage disequilibrium (LD) patterns with observed HLA class II DRB1-DQA1-DQB1 haplotype frequencies. Expected frequencies were compared to observed allele frequencies in patients. RESULTS Significant type 1 diabetes associations were observed at all class I HLA loci. After accounting for LD with HLA class II, the most significantly type 1 diabetes–associated alleles were B*5701 (odds ratio 0.19; P = 4 × 10−11) and B*3906 (10.31; P = 4 × 10−10). Other significantly type 1 diabetes–associated alleles included A*2402, A*0201, B*1801, and C*0501 (predisposing) and A*1101, A*3201, A*6601, B*0702, B*4403, B*3502, C*1601, and C*0401 (protective). Some alleles, notably B*3906, appear to modulate the risk of all DRB1-DQA1-DQB1 haplotypes on which they reside, suggesting a class I effect that is independent of class II. Other class I type 1 diabetes associations appear to be specific to individual class II haplotypes. Some apparent associations (e.g., C*1601) could be attributed to strong LD to another class I susceptibility locus (B*4403). CONCLUSIONS These data indicate that HLA class I alleles, in addition to and independently from HLA class II alleles, are associated with type 1 diabetes.

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−16). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10−12), *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−13. 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.

HLA genotyping in the international Type 1 Diabetes Genetics Consortium

Background Although human leukocyte antigen (HLA) DQ and DR loci appear to confer the strongest genetic risk for type 1 diabetes, more detailed information is required for other loci within the HLA region to understand causality and stratify additional risk factors. The Type 1 Diabetes Genetics Consortium (T1DGC) study design included high-resolution genotyping of HLA-A, B, C, DRB1, DQ, and DP loci in all affected sibling pair and trio families, and cases and controls, recruited from four networks worldwide, for analysis with clinical phenotypes and immunological markers. Purpose In this article, we present the operational strategy of training, classification, reporting, and quality control of HLA genotyping in four laboratories on three continents over nearly 5 years. Methods Methods to standardize HLA genotyping at eight loci included: central training and initial certification testing; the use of uniform reagents, protocols, instrumentation, and software versions; an automated data transfer; and the use of standardized nomenclature and allele databases. We implemented a rigorous and consistent quality control process, reinforced by repeated workshops, yearly meetings, and telephone conferences. Results A total of 15,246 samples have been HLA genotyped at eight loci to four-digit resolution; an additional 6797 samples have been HLA genotyped at two loci. The genotyping repeat rate decreased significantly over time, with an estimated unresolved Mendelian inconsistency rate of 0.21%. Annual quality control exercises tested 2192 genotypes (4384 alleles) and achieved 99.82% intra-laboratory and 99.68% inter-laboratory concordances. Limitations The chosen genotyping platform was unable to distinguish many allele combinations, which would require further multiple stepwise testing to resolve. For these combinations, a standard allele assignment was agreed upon, allowing further analysis if required. Conclusions High-resolution HLA genotyping can be performed in multiple laboratories using standard equipment, reagents, protocols, software, and communication to produce consistent and reproducible data with minimal systematic error. Many of the strategies used in this study are generally applicable to other large multi-center studies. Clinical Trials 2010; 7: S75—S87. http:// ctj.sagepub.com

HLA Class II Genotyping of African American Type 1 Diabetic Patients Reveals Associations Unique to African Haplotypes

HLA genotyping was performed in African American type 1 diabetic patients (n = 772) and controls (n = 1,641) in the largest study of African Americans and type 1 diabetes reported to date. Cases were from Children’s Hospital and Research Center Oakland and from existing collections (Type 1 Diabetes Genetics Consortium [T1DGC], Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications [DCCT/EDIC], and Genetics of Kidneys in Diabetes [GoKinD]). Controls were from the T1DGC and from newborn bloodspot cards. The diversity of HLA DRB1-DQA1-DQB1 haplotypes and genotypes is far greater than that found in Europeans and European Americans. Association analyses replicated many type 1 diabetes risk effects of European-derived haplotypes but also revealed novel effects for African-derived haplotypes. Notably, the African-specific “DR3” haplotype DRB1*03:02-DQA1*04:01-DQB1*04:02 is protective for type 1 diabetes, in contrast to the common and highly-susceptible DR3 DRB1*03:01-DQA1*05:01-DQB1*02:01. Both DRB1*07:01 and DRB1*13:03 haplotypes are predisposing when they include DQA1*03:01-DQB1*02:01g but are protective with DQA1*02:01-DQB1*02:01g. The heterozygous DR4/DR9 genotype, containing the African-derived “DR9” haplotype DRB1*09:01-DQA1*03:01-DQB1*02:01g, exhibits extremely high risk (odds ratio = 30.88), approaching that for DR3/DR4 in European populations. Disease risk assessment for African Americans differs greatly from risk assessment in European populations. This has profound implications on risk screening programs and underscores the need for high-resolution genotyping of multiple populations for the rational design of screening programs with tests that will fairly represent the population being screened.

Hypomethylation within gene promoter regions and type 1 diabetes in discordant monozygotic twins

Genetic susceptibility to type 1 diabetes (T1D) is well supported by epidemiologic evidence; however, disease risk cannot be entirely explained by established genetic variants identified so far. This study addresses the question of whether epigenetic modification of the inherited DNA sequence may contribute to T1D susceptibility. Using the Infinium HumanMethylation450 BeadChip array (450k), a total of seven long-term disease-discordant monozygotic (MZ) twin pairs and five pairs of HLA-identical, disease-discordant non-twin siblings (NTS) were examined for associations between DNA methylation (DNAm) and T1D. Strong evidence for global hypomethylation of CpG sites within promoter regions in MZ twins with TID compared to twins without T1D was observed. DNA methylation data were then grouped into three categories of CpG sites for further analysis, including those within: 1) the major histocompatibility complex (MHC) region, 2) non-MHC genes with reported T1D association through genome wide association studies (GWAS), and 3) the epigenome, or remainder of sites that did not include MHC and T1D associated genes. Initial results showed modest methylation differences between discordant MZ twins for the MHC region and T1D-associated CpG sites, BACH2, INS-IGF2, and CLEC16A (DNAm difference range: 2.2%-5.0%). In the epigenome CpG set, the greatest methylation differences were observed in MAGI2, FANCC, and PCDHB16, (DNAm difference range: 6.9%-16.1%). These findings were not observed in the HLA-identical NTS pairs. Targeted pyrosequencing of five candidate CpG loci identified using the 450k array in the original discordant MZ twins produced similar results using control DNA samples, indicating strong agreement between the two DNA methylation profiling platforms. However, findings for the top five candidate CpG loci were not replicated in six additional T1D-discordant MZ twin pairs. Our results indicate global DNA hypomethylation within gene promoter regions may contribute to T1D; however, findings do not support the involvement of large DNAm differences at single CpG sites alone in T1D.

The rs4774 CIITA missense variant is associated with risk of systemic lupus erythematosus

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor required for human leukocyte antigens (HLA) class II MHC-restricted antigen presentation. MHC genes, including the HLA class II DRB1*03:01 allele, are strongly associated with systemic lupus erythematosus (SLE). Recently the rs4774 CIITA missense variant (+1632G/C) was reported to be associated with susceptibility to multiple sclerosis. In the current study, we investigated CIITA, DRB1*03:01 and risk of SLE using a multi-stage analysis. In stage 1, 9 CIITA variants were tested in 658 cases and 1363 controls (N=2021). In stage 2, rs4774 was tested in 684 cases and 2938 controls (N=3622). We also performed a meta-analysis of the pooled 1342 cases and 4301 controls (N=5643). In stage 1, rs4774*C was associated with SLE (odds ratio (OR)=1.24, 95% confidence interval (95% CI)=1.07–1.44, P=4.2 × 10−3). Similar results were observed in stage 2 (OR=1.16, 95% CI=1.02–1.33, P=8.5 × 10−3) and the meta-analysis of the combined data set (OR=1.20, 95% CI=1.09–1.33, Pmeta=2.5 × 10−4). In all three analyses, the strongest evidence for association between rs4774*C and SLE was present in individuals who carried at least one copy of DRB1*03:01 (Pmeta=1.9 × 10−3). Results support a role for CIITA in SLE, which appears to be stronger in the presence of DRB1*03:01.

Type 1 diabetes risk for human leukocyte antigen (HLA)-DR3 haplotypes depends on genotypic context: association of DPB1 and HLA class I loci among DR3- and DR4-matched Italian patients and controls.

Patients with high-risk human leukocyte antigen (HLA)-DR-DQ genotypes for type 1 diabetes (T1D) were compared with HLA-matched controls to evaluate T1D risk for other HLA loci, including HLA-A, -B, -Cw, and DPB1. Patients (n = 133) with high-risk genotypes (DR3/DR3, DR3/DR4, DR4/DR4) were selected from the Lazio (Rome) region of Italy. Screening of more than 9000 patients from the Lazio region and northern Italy yielded 162 controls with high-T1D-risk haplotypes. Although the overall distributions did not differ significantly, allele frequency differences were discovered between the controls from Lazio and controls from northern Italy for some alleles previously determined to affect T1D risk, such as A*3002, DPB1*0301, and DPB1*0402. Therefore, Lazio patient data were compared both with the Lazio subset of controls (n = 53) and with the entire group of controls for association analyses. Significant allele frequency differences between patients and DR-DQ-matched controls existed for specific alleles at all loci. Data for the DR3/DR3 subset of patients and controls demonstrated an increase of Cw*0702 in patients. Compared with controls, reduced patient frequencies were seen for several alleles, including A*0101, B*0801, and Cw*0701, all on the highly conserved, extended DR3 haplotype known as 8.1 in DR3/DR3, but not DR3/DR4, subgroup. DPB1*0101, often reported on 8.1 haplotypes, was also less frequent in DR3/DR3 patients than controls. Analysis of family-based data from the HBDI repository was consistent with the observed results from the Italian patients, indicating the presence of a T1D-protective locus at or near A*0101 and a second T1D-protective locus at or near DPB1*0101. These data indicate that T1D risk conferred by the 8.1 haplotype is genotype dependent.

KIR haplotypes are associated with late-onset type 1 diabetes in European–American families

Classical human leukocyte antigens (HLA) genes confer the strongest, but not the only, genetic susceptibility to type 1 diabetes. Killer cell immunoglobulin-like receptors (KIR), on natural killer (NK) cells, bind ligands including class I HLA. We examined presence or absence, with copy number, of KIR loci in 1698 individuals, from 339 multiplex type 1 diabetes families, from the Human Biological Data Interchange, previously genotyped for HLA. Combining family data with KIR copy number information allowed assignment of haplotypes using identity by descent. This is the first disease study to use KIR copy number typing and unambiguously define haplotypes by gene transmission. KIR A1 haplotypes were positively associated with T1D in the subset of patients without the high T1D risk HLA genotype, DR3/DR4 (odds ratio=1.29, P=0.0096). The data point to a role for KIR in type 1 diabetes risk in late-onset patients. In the top quartile (age of onset>14), KIR A2 haplotype was overtransmitted (63.4%, odds ratio=1.73, P=0.024) and KIR B haplotypes were undertransmitted (41.1%, odds ratio=0.70, P=0.0052) to patients. The data suggest that inhibitory ‘A’ haplotypes are predisposing and stimulatory ‘B’ haplotypes confer protection in both DR3/DR4-negative and late-onset patient groups.

Maximizing deoxyribonucleic acid yield from dried blood spots.

Background: One source of deoxyribonucleic acid (DNA) for genetic studies is the utilization of dried blood spots stored on paper cards (Guthrie cards) collected shortly after birth. These cards represent an important source of material for epidemiologic and population-based genetic studies. Extraction of DNA from these cards can lead to variable amounts of recovered DNA. We report here results of our efforts to maximize yield from this valuable, but nonrenewable, resource. Method: Commercial methods of DNA extraction from blood cards were used, and protocol modifications were introduced that enhanced DNA yield. Results: Use of a commercial solvent prior to DNA extraction steps gave greater yields than extraction without the solvent. Modification of the elution step by use of prewarmed extraction buffer and a soaking step at an elevated temperature increased yield by 6- to 10-fold. Conclusions: The modified DNA extraction method yielded as much as 660 ng of DNA from a single 5-mm-diameter punch of a blood spot card. The DNA performed well in downstream, polymerase chain reaction-based applications.

Analysis of Maternal–Offspring HLA Compatibility, Parent-of-Origin Effects, and Noninherited Maternal Antigen Effects for HLA–DRB1 in Systemic Lupus Erythematosus

OBJECTIVE Genetic susceptibility to systemic lupus erythematosus (SLE) is well established, with the HLA class II DRB1 and DQB1 loci demonstrating the strongest association. However, HLA may also influence SLE through novel biologic mechanisms in addition to genetic transmission of risk alleles. Evidence for increased maternal-offspring HLA class II compatibility in SLE and differences in maternal versus paternal transmission rates (parent-of-origin effects) and nontransmission rates (noninherited maternal antigen [NIMA] effects) in other autoimmune diseases have been reported. Thus, we investigated maternal-offspring HLA compatibility, parent-of-origin effects, and NIMA effects at DRB1 in SLE. METHODS The cohort comprised 707 SLE families and 188 independent healthy maternal-offspring pairs (total of 2,497 individuals). Family-based association tests were conducted to compare transmitted versus nontransmitted alleles (transmission disequilibrium test) and both maternally versus paternally transmitted (parent-of-origin) and nontransmitted alleles (using the chi-square test of heterogeneity). Analyses were stratified according to the sex of the offspring. Maternally affected offspring DRB1 compatibility in SLE families was compared with paternally affected offspring compatibility and with independent control maternal-offspring pairs (using Fishers test) and was restricted to male and nulligravid female offspring with SLE. RESULTS As expected, DRB1 was associated with SLE (P < 1 x 10(-4)). However, mothers of children with SLE had similar transmission and nontransmission frequencies for DRB1 alleles when compared with fathers, including those for the known SLE risk alleles HLA-DRB1*0301, *1501, and *0801. No association between maternal-offspring compatibility and SLE was observed. CONCLUSION Maternal-offspring HLA compatibility, parent-of-origin effects, and NIMA effects at DRB1 are unlikely to play a role in SLE.