Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis

The genetic association of the major histocompatibility complex (MHC) to rheumatoid arthritis risk has commonly been attributed to alleles in HLA-DRB1. However, debate persists about the identity of the causal variants in HLA-DRB1 and the presence of independent effects elsewhere in the MHC. Using existing genome-wide SNP data in 5,018 individuals with seropositive rheumatoid arthritis (cases) and 14,974 unaffected controls, we imputed and tested classical alleles and amino acid polymorphisms in HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1 and HLA-DRB1, as well as 3,117 SNPs across the MHC. Conditional and haplotype analyses identified that three amino acid positions (11, 71 and 74) in HLA-DRβ1 and single-amino-acid polymorphisms in HLA-B (at position 9) and HLA-DPβ1 (at position 9), which are all located in peptide-binding grooves, almost completely explain the MHC association to rheumatoid arthritis risk. This study shows how imputation of functional variation from large reference panels can help fine map association signals in the MHC.

High HLA-DP Expression and Graft-versus-Host Disease

BACKGROUND

Transplantation of hematopoietic cells from unrelated donors can cure blood disorders but carries a significant risk of acute graft-versus-host disease (GVHD). The risk is higher when the recipient and donor are HLA-DPB1-mismatched, but the mechanisms leading to GVHD are unknown. The HLA-DPB1 regulatory region variant rs9277534 is associated with HLA-DPB1 expression. We tested the hypothesis that the GVHD risk correlates with the rs9277534 allele linked to the mismatched HLA-DPB1 in the recipient.

METHODS

We genotyped rs9277534 in 3505 persons to define rs9277534-DPB1 haplotypes. Among 1441 recipients of transplants from HLA-A,B,C,DRB1,DQB1-matched unrelated donors with only one HLA-DPB1 mismatch, linkage of the rs9277534 A and G alleles to the mismatched HLA-DPB1 was determined. HLA-DPB1 expression was assessed by means of a quantitative polymerase-chain-reaction assay. The risk of acute GVHD among recipients whose mismatched HLA-DPB1 allele was linked to rs9277534G (high expression) was compared with the risk among recipients whose mismatched HLA-DPB1 allele was linked to rs9277534A (low expression).

RESULTS

The mean HLA-DPB1 expression was lower with rs9277534A than with rs9277534G. Among recipients of transplants from donors with rs9277534A-linked HLA-DPB1, the risk of acute GVHD was higher for recipients with rs9277534G-linked HLA-DPB1 mismatches than for recipients with rs9277534A-linked HLA-DPB1 mismatches (hazard ratio, 1.54; 95% confidence interval [CI], 1.25 to 1.89; P<0.001), as was the risk of death due to causes other than disease recurrence (hazard ratio, 1.25; 95% CI, 1.00 to 1.57; P=0.05).

CONCLUSIONS

The risk of GVHD associated with HLA-DPB1 mismatching was influenced by the HLA-DPB1 rs9277534 expression marker. Among recipients of HLA-DPB1-mismatched transplants from donors with the low-expression allele, recipients with the high-expression allele had a high risk of GVHD. (Funded by the National Institutes of Health and others.).

Fine-mapping the genetic association of the major histocompatibility complex in multiple sclerosis: HLA and non-HLA effects

The major histocompatibility complex (MHC) region is strongly associated with multiple sclerosis (MS) susceptibility. HLA-DRB1*15:01 has the strongest effect, and several other alleles have been reported at different levels of validation. Using SNP data from genome-wide studies, we imputed and tested classical alleles and amino acid polymorphisms in 8 classical human leukocyte antigen (HLA) genes in 5,091 cases and 9,595 controls. We identified 11 statistically independent effects overall: 6 HLA-DRB1 and one DPB1 alleles in class II, one HLA-A and two B alleles in class I, and one signal in a region spanning from MICB to LST1. This genomic segment does not contain any HLA class I or II genes and provides robust evidence for the involvement of a non-HLA risk allele within the MHC. Interestingly, this region contains the TNF gene, the cognate ligand of the well-validated TNFRSF1A MS susceptibility gene. The classical HLA effects can be explained to some extent by polymorphic amino acid positions in the peptide-binding grooves. This study dissects the independent effects in the MHC, a critical region for MS susceptibility that harbors multiple risk alleles.

A novel variant marking HLA-DP expression levels predicts recovery from hepatitis B virus infection

Variants near the HLA-DP gene show the strongest genome-wide association with chronic hepatitis B virus (HBV) infection and HBV recovery/persistence in Asians. To test the effect of the HLA-DP region on outcomes to HBV infection, we sequenced the polymorphic HLA-DPB1 and DPA1 coding exons and the corresponding 3' untranslated regions (3'UTRs) in 662 individuals of European-American and African-American ancestry. The genome-wide association study (GWAS) variant (rs9277535; 550A/G) in the 3'UTR of the HLA-DPB1 gene that associated most significantly with chronic hepatitis B and outcomes to HBV infection in Asians had a marginal effect on HBV recovery in our European- and African-American samples (odds ratio [OR] = 0.39, P = 0.01, combined ethnic groups). However, we identified a novel variant in the HLA-DPB1 3'UTR region, 496A/G (rs9277534), which associated very significantly with HBV recovery in both European and African-American populations (OR = 0.37, P = 0.0001, combined ethnic groups). The 496A/G variant distinguishes the most protective HLA-DPB1 allele (DPB1*04:01) from the most susceptible (DPB1*01:01), whereas 550A/G does not. 496A/G has a stronger effect than any individual HLA-DPB1 or DPA1 allele and any other HLA alleles that showed an association with HBV recovery in our European-American cohort. The 496GG genotype, which confers recessive susceptibility to HBV persistence, also associates in a recessive manner with significantly higher levels of HLA-DP surface protein and transcript level expression in healthy donors, suggesting that differences in expression of HLA-DP may increase the risk of persistent HBV infection.

Genome-wide association study confirming association of HLA-DP with protection against chronic hepatitis B and viral clearance in Japanese and Korean

Hepatitis B virus (HBV) infection can lead to serious liver diseases, including liver cirrhosis (LC) and hepatocellular carcinoma (HCC); however, about 85-90% of infected individuals become inactive carriers with sustained biochemical remission and very low risk of LC or HCC. To identify host genetic factors contributing to HBV clearance, we conducted genome-wide association studies (GWAS) and replication analysis using samples from HBV carriers and spontaneously HBV-resolved Japanese and Korean individuals. Association analysis in the Japanese and Korean data identified the HLA-DPA1 and HLA-DPB1 genes with P(meta) = 1.89×10⁻¹² for rs3077 and P(meta) = 9.69×10⁻¹⁰ for rs9277542. We also found that the HLA-DPA1 and HLA-DPB1 genes were significantly associated with protective effects against chronic hepatitis B (CHB) in Japanese, Korean and other Asian populations, including Chinese and Thai individuals (P(meta) = 4.40×10⁻¹⁹ for rs3077 and P(meta) = 1.28×10⁻¹⁵ for rs9277542). These results suggest that the associations between the HLA-DP locus and the protective effects against persistent HBV infection and with clearance of HBV were replicated widely in East Asian populations; however, there are no reports of GWAS in Caucasian or African populations. Based on the GWAS in this study, there were no significant SNPs associated with HCC development. To clarify the pathogenesis of CHB and the mechanisms of HBV clearance, further studies are necessary, including functional analyses of the HLA-DP molecule.

Association of granulomatosis with polyangiitis (Wegener's) with HLA-DPB1*04 and SEMA6A gene variants: evidence from genome-wide analysis

OBJECTIVE

To identify genetic determinants of granulomatosis with polyangiitis (Wegener's) (GPA).

METHODS

We carried out a genome-wide association study (GWAS) of 492 GPA cases and 1,506 healthy controls (white subjects of European descent), followed by replication analysis of the most strongly associated signals in an independent cohort of 528 GPA cases and 1,228 controls.

RESULTS

Genome-wide significant associations were identified in 32 single-nucleotide polymorphic (SNP) markers across the HLA region, the majority of which were located in the HLA-DPB1 and HLA-DPA1 genes encoding the class II major histocompatibility complex (MHC) DPβ chain 1 and DPα chain 1 proteins, respectively. Peak association signals in these 2 genes, emanating from SNPs rs9277554 (for DPβ chain 1) and rs9277341 (DPα chain 1) were strongly replicated in an independent cohort (in the combined analysis of the initial cohort and the replication cohort, P = 1.92 × 10(-50) and 2.18 × 10(-39) , respectively). Imputation of classic HLA alleles and conditional analyses revealed that the SNP association signal was fully accounted for by the classic HLA-DPB1*04 allele. An independent single SNP, rs26595, near SEMA6A (the gene for semaphorin 6A) on chromosome 5, was also associated with GPA, reaching genome-wide significance in a combined analysis of the GWAS and replication cohorts (P = 2.09 × 10(-8) ).

CONCLUSION

We identified the SEMA6A and HLA-DP loci as significant contributors to risk for GPA, with the HLA-DPB1*04 allele almost completely accounting for the MHC association. These two associations confirm the critical role of immunogenetic factors in the development of GPA.

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.

Risk alleles for chronic hepatitis B are associated with decreased mRNA expression of HLA-DPA1 and HLA-DPB1 in normal human liver

A genome-wide association study identified single nucleotide polymorphisms (SNPs) rs3077 and rs9277535 located in the 3' untranslated regions of human leukocyte antigen (HLA) class II genes HLA-DPA1 and HLA-DPB1, respectively, as the independent variants most strongly associated with chronic hepatitis B. We examined whether these SNPs are associated with mRNA expression of HLA-DPA1 and HLA-DPB1. We identified gene expression-associated SNPs (eSNPs) in normal liver samples obtained from 651 individuals of European ancestry by integrating genotype (~650 000 SNPs) and gene expression (>39 000 transcripts) data from each sample. We used the Kruskal-Wallis test to determine associations between gene expression and genotype. To confirm findings, we measured allelic expression imbalance (AEI) of complementary DNA compared with DNA in liver specimens from subjects who were heterozygous for rs3077 and rs9277535. On a genome-wide basis, rs3077 was the SNP most strongly associated with HLA-DPA1 expression (p=10(-48)), and rs9277535 was strongly associated with HLA-DPB1 expression (p=10(-15)). Consistent with these gene expression associations, we observed AEI for both rs3077 (p=3.0 × 10(-7); 17 samples) and rs9277535 (p=0.001; 17 samples). We conclude that the variants previously associated with chronic hepatitis B are also strongly associated with mRNA expression of HLA-DPA1 and HLA-DPB1, suggesting that expression of these genes is important in control of HBV.

Effect of host and viral factors on hepatitis B e antigen-positive chronic hepatitis B patients receiving pegylated interferon-α-2a therapy

BACKGROUND

Pegylated interferon (PEG-IFN)-α-2a improves the hepatitis B e antigen (HBeAg) seroconversion rate in HBeAg-positive chronic hepatitis B patients. However, baseline factors predicting favourable responses to PEG-IFN-α-2a remain largely unknown.

METHODS

A total of 115 HBeAg-positive chronic hepatitis B patients who had a pre-therapy serum alanine aminotransferase (ALT) level over two times the upper limit of normal and received PEG-IFN-α-2a for 6-12 months were consecutively enrolled according to the local reimbursed guidelines. HBeAg seroconversion and combined response defined as HBeAg seroconversion, HBV-DNA level <20,000 IU/ml as well as ALT normalization at 6 months off therapy were primary and secondary therapeutic end points, respectively. Baseline viral factors, including viral load, genotype and major sequences of precore stop codon/basal core promoter (BCP), and host factors, including three single nucleotide polymorphisms among the HLA-DPA1, HLA-DPB1 and IL28B regions, were determined to correlate with therapeutic end points.

RESULTS

HBeAg seroconversion and combined response rates were 26.1% and 18.3%, respectively. By multivariate analysis, BCP mutation (OR 8.04, 95% CI 2.00-32.28) and rs3077 G/G genotype (OR 3.49, 95% CI 1.12-10.84) were associated with a higher HBeAg seroconversion rate; BCP mutation (OR 9.28, 95% CI 1.92-44.99) and baseline viral load <2 × 10(6) IU/ml (OR 4.78, 95% CI 1.37-16.69) were associated with a higher combined response rate.

CONCLUSIONS

BCP mutation is associated with higher HBeAg seroconversion and combined response rates at 6 months off therapy in HBeAg-positive chronic hepatitis B patients treated with PEG-IFN-α-2a. Genetic variants in the HLA-DPA1 region may also affect treatment-induced HBeAg seroconversion.

Classical HLA-DRB1 and DPB1 alleles account for HLA associations with primary biliary cirrhosis

Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with human leukocyte antigen (HLA)-region polymorphisms. To determine if associations can be explained by classical HLA determinants, we studied Italian, 676 cases and 1440 controls, genotyped with dense single-nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (*08, *11 and *14) account for most of the DRB1 association signal. Consistent with previous studies, DRB1*08 (P=1.59 × 10(-11)) was the strongest predisposing allele, whereas DRB1*11 (P=1.42 × 10(-10)) was protective. Additionally, DRB1*14 and the DPB1 association (DPB1*03:01; P=9.18 × 10(-7)) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC.

Genome-wide association study of aspirin-exacerbated respiratory disease in a Korean population

Aspirin-exacerbated respiratory disease (AERD) is a nonallergic clinical syndrome characterized by a severe decline in forced expiratory volume in one second (FEV1) following the ingestion of non-steroidal anti-inflammatory drugs (NSAIDs) such as aspirin. The effects of genetic variants have not fully explained all of the observed individual differences to an aspirin challenge despite previous attempts to identify AERD-related genes. In the present study, we performed genome-wide association study (GWAS) and targeted association study in Korean asthmatics to identify new genetic factors associated with AERD. A total of 685 asthmatic patients without AERD and 117 subjects with AERD were used for the GWAS of the first stage, and 996 asthmatics without AERD and 142 subjects with AERD were used for a follow-up study. A total of 702 SNPs were genotyped using the GoldenGate assay with the VeraCode microbead. GWAS revealed the top-ranked variants in 3' regions of the HLA-DPB1 gene. To investigate the detailed genetic effects of an associated region with the risk of AERD, a follow-up targeted association study with the 702 single nucleotide polymorphisms (SNPs) of 14 genes was performed on 802 Korean subjects. In a case-control analysis, HLA-DPB1 rs1042151 (Met105Val) shows the most significant association with the susceptibility of AERD (p = 5.11 × 10(-7); OR = 2.40). Moreover, rs1042151 also shows a gene dose for the percent decline of FEV1 after an aspirin challenge (p = 2.82 × 10(-7)). Our findings show that the HLA-DPB1 gene polymorphism may be the most susceptible genetic factor for the risk of AERD in Korean asthmatics and confirm the importance of HLA-DPB1 in the genetic etiology of AERD.

Transancestral mapping of the MHC region in systemic lupus erythematosus identifies new independent and interacting loci at MSH5, HLA-DPB1 and HLA-G

OBJECTIVES

Systemic lupus erythematosus (SLE) is a chronic multisystem genetically complex autoimmune disease characterised by the production of autoantibodies to nuclear and cellular antigens, tissue inflammation and organ damage. Genome-wide association studies have shown that variants within the major histocompatibility complex (MHC) region on chromosome 6 confer the greatest genetic risk for SLE in European and Chinese populations. However, the causal variants remain elusive due to tight linkage disequilibrium across disease-associated MHC haplotypes, the highly polymorphic nature of many MHC genes and the heterogeneity of the SLE phenotype.

METHODS

A high-density case-control single nucleotide polymorphism (SNP) study of the MHC region was undertaken in SLE cohorts of Spanish and Filipino ancestry using a custom Illumina chip in order to fine-map association signals in these haplotypically diverse populations. In addition, comparative analyses were performed between these two datasets and a northern European UK SLE cohort. A total of 1433 cases and 1458 matched controls were examined.

RESULTS

Using this transancestral SNP mapping approach, novel independent loci were identified within the MHC region in UK, Spanish and Filipino patients with SLE with some evidence of interaction. These loci include HLA-DPB1, HLA-G and MSH5 which are independent of each other and HLA-DRB1 alleles. Furthermore, the established SLE-associated HLA-DRB1*15 signal was refined to an interval encompassing HLA-DRB1 and HLA-DQA1. Increased frequencies of MHC region risk alleles and haplotypes were found in the Filipino population compared with Europeans, suggesting that the greater disease burden in non-European SLE may be due in part to this phenomenon.

CONCLUSION

These data highlight the usefulness of mapping disease susceptibility loci using a transancestral approach, particularly in a region as complex as the MHC, and offer a springboard for further fine-mapping, resequencing and transcriptomic analysis.

Association of human leukocyte antigen with interstitial lung disease in rheumatoid arthritis: a protective role for shared epitope

Introduction

Interstitial Lung Disease (ILD) is frequently associated with Rheumatoid Arthritis (RA) as one of extra-articular manifestations. Many studies for Human Leukocyte Antigen (HLA) allelic association with RA have been reported, but few have been validated in an RA subpopulation with ILD. In this study, we investigated the association of HLA class II alleles with ILD in RA.

Methods

An association study was conducted on HLA-DRB1, DQB1, and DPB1 in 450 Japanese RA patients that were or were not diagnosed with ILD, based on the findings of computed tomography images of the chest.

Results

Unexpectedly, HLA-DRB1*04 (corrected P [Pc] = 0.0054, odds ratio [OR] 0.57), shared epitope (SE) (P = 0.0055, OR 0.66) and DQB1*04 (Pc = 0.0036, OR 0.57) were associated with significantly decreased risk of ILD. In contrast, DRB1*16 (Pc = 0.0372, OR 15.21), DR2 serological group (DRB1*15 and *16 alleles) (P = 0.0020, OR 1.75) and DQB1*06 (Pc = 0.0333, OR 1.57, respectively) were significantly associated with risk of ILD.

Conclusion

HLA-DRB1 SE was associated with reduced, while DR2 serological group (DRB1*15 and *16) with increased, risk for ILD in Japanese patients with RA.

Genetic and infectious profiles of Japanese multiple sclerosis patients

BACKGROUND

Nationwide surveys conducted in Japan over the past thirty years have revealed a four-fold increase in the estimated number of multiple sclerosis (MS) patients, a decrease in the age at onset, and successive increases in patients with conventional MS, which shows an involvement of multiple sites in the central nervous system, including the cerebrum and cerebellum. We aimed to clarify whether genetic and infectious backgrounds correlate to distinct disease phenotypes of MS in Japanese patients.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed HLA-DRB1 and -DPB1 alleles, and IgG antibodies specific for Helicobacter pylori, Chlamydia pneumoniae, varicella zoster virus, and Epstein-Barr virus nuclear antigen (EBNA) in 145 MS patients and 367 healthy controls (HCs). Frequencies of DRB1*0405 and DPB1*0301 were significantly higher, and DRB1*0901 and DPB1*0401 significantly lower, in MS patients as compared with HCs. MS patients with DRB1*0405 had a significantly earlier age of onset and lower Progression Index than patients without this allele. The proportion and absolute number of patients with DRB1*0405 successively increased with advancing year of birth. In MS patients without DRB1*0405, the frequency of the DRB1*1501 allele was significantly higher, while the DRB1*0901 allele was significantly lower, compared with HCs. Furthermore, DRB1*0405-negative MS patients were significantly more likely to be positive for EBNA antibodies compared with HCs.

CONCLUSIONS

Our study suggests that MS patients harboring DRB1*0405, a genetic risk factor for MS in the Japanese population, have a younger age at onset and a relatively benign disease course, while DRB1*0405-negative MS patients have features similar to Western-type MS in terms of association with Epstein-Barr virus infection and DRB1*1501. The recent increase of MS in young Japanese people may be caused, in part, by an increase in DRB1*0405-positive MS patients.

A weighted genetic risk score using all known susceptibility variants to estimate rheumatoid arthritis risk

BACKGROUND

There is currently great interest in the incorporation of genetic susceptibility loci into screening models to identify individuals at high risk of disease. Here, we present the first risk prediction model including all 46 known genetic loci associated with rheumatoid arthritis (RA).

METHODS

A weighted genetic risk score (wGRS) was created using 45 RA non-human leucocyte antigen (HLA) susceptibility loci, imputed amino acids at HLA-DRB1 (11, 71 and 74), HLA-DPB1 (position 9) HLA-B (position 9) and gender. The wGRS was tested in 11 366 RA cases and 15 489 healthy controls. The risk of developing RA was estimated using logistic regression by dividing the wGRS into quintiles. The ability of the wGRS to discriminate between cases and controls was assessed by receiver operator characteristic analysis and discrimination improvement tests.

RESULTS

Individuals in the highest risk group showed significantly increased odds of developing anti-cyclic citrullinated peptide-positive RA compared to the lowest risk group (OR 27.13, 95% CI 23.70 to 31.05). The wGRS was validated in an independent cohort that showed similar results (area under the curve 0.78, OR 18.00, 95% CI 13.67 to 23.71). Comparison of the full wGRS with a wGRS in which HLA amino acids were replaced by a HLA tag single-nucleotide polymorphism showed a significant loss of sensitivity and specificity.

CONCLUSIONS

Our study suggests that in RA, even when using all known genetic susceptibility variants, prediction performance remains modest; while this is insufficiently accurate for general population screening, it may prove of more use in targeted studies. Our study has also highlighted the importance of including HLA variation in risk prediction models.

Host genetic factors affecting spontaneous HBsAg seroclearance in chronic hepatitis B patients

Spontaneous clearance of hepatitis B surface antigen (HBsAg) in chronic hepatitis B (CHB) patients usually indicates a remission of hepatitis activity and a favorable outcome. Two single nucleotide polymorphisms (SNP), rs3077 near HLA-DPA1 region and rs9277535 near HLA-DPB1 region, have been shown to be associated with HBV persistence after acute HBV infection. However, little is known about the impact of these 2 SNPs on spontaneous HBsAg clearance in CHB patients. In this case-control study, a total of 100 male HBeAg-negative chronic HBV carriers who cleared HBsAg spontaneously (case group) and 100 age-matched HBeAg-negative male patients with persistent HBsAg positivity (control group) were enrolled. We investigated the relationship between these 2 SNPs and HBsAg clearance. There was a higher frequency of rs9277535 non-GG genotype in the case group (57% vs. 42%). Patients with rs9277535 non-GG genotype had a higher chance to clear HBsAg [Odds ratio (OR): 1.83, 95% confidence interval (CI): 1.04∼3.21, P = 0.034]. Compared to GG haplotype of rs3077 and rs9277535, GA haplotype had a higher chance of achieving spontaneous HBsAg loss (OR: 2.17, 95% CI: 1.14∼4.16, P = 0.030). In conclusion, rs9277535 non-GG genotype is associated with a higher likelihood of spontaneous HBsAg seroclearance in CHB patients.

Association between GWAS-identified lung adenocarcinoma susceptibility loci and EGFR mutations in never-smoking Asian women, and comparison with findings from Western populations

To evaluate associations by EGFR mutation status for lung adenocarcinoma risk among never-smoking Asian women, we conducted a meta-analysis of 11 loci previously identified in genome-wide association studies (GWAS). Genotyping in an additional 10,780 never-smoking cases and 10,938 never-smoking controls from Asia confirmed associations with eight known single nucleotide polymorphisms (SNPs). Two new signals were observed at genome-wide significance (P < 5 × 10-8), namely, rs7216064 (17q24.3, BPTF), for overall lung adenocarcinoma risk, and rs3817963 (6p21.3, BTNL2) which is specific to cases with EGFR mutations. In further sub-analyses by EGFR status, rs9387478 (ROS1/DCBLD1) and rs2179920 (HLA-DPB1) showed stronger estimated associations in EGFR-positive compared to EGFR-negative cases. Comparison of the overall associations with published results in Western populations revealed that the majority of these findings were distinct, underscoring the importance of distinct contributing factors for smoking and non-smoking lung cancer. Our results extend the catalogue of regions associated with lung adenocarcinoma in non-smoking Asian women and highlight the importance of how the germline could inform risk for specific tumour mutation patterns, which could have important translational implications.

Predicting an HLA-DPB1 expression marker based on standard DPB1 genotyping: Linkage analysis of over 32,000 samples

The risk of acute graft-versus-host disease (GvHD) after hematopoietic stem cell transplantation is increased with donor-recipient HLA-DPB1 allele mismatching. The single-nucleotide polymorphism (SNP) rs9277534 within the 3' untranslated region (UTR) correlates with HLA-DPB1 allotype expression and serves as a marker for permissive HLA-DPB1 mismatches. Since rs9277534 is not routinely typed, we analyzed 32,681 samples of mostly European ancestry to investigate if the rs9277534 allele can be reliably imputed from standard DPB1 genotyping. We confirmed the previously-defined linkages between rs9277534 and 18 DPB1 alleles and established additional linkages for 46 DPB1 alleles. Based on these linkages, the rs9277534 allele could be predicted for 99.6% of the samples based on DPB1 genotypes (99.99% concordance). We demonstrate that 100% prediction accuracy could be achieved if the prediction utilized exon 3 sequence information. DPB1 genotyping based on exon 2 data alone allows no unambiguous rs9277534 allele prediction but was estimated to maintain 99% accuracy for samples of European descent. We conclude that DPB1 genotyping is sufficient to infer the DPB1 expression marker rs9277534 with high accuracy. This information could be used to select donors with permissive HLA-DPB1 mismatches without directly screening for rs9277534.

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.

Genome-wide association study of classical Hodgkin lymphoma identifies key regulators of disease susceptibility

Several susceptibility loci for classical Hodgkin lymphoma have been reported. However, much of the heritable risk is unknown. Here, we perform a meta-analysis of two existing genome-wide association studies, a new genome-wide association study, and replication totalling 5,314 cases and 16,749 controls. We identify risk loci for all classical Hodgkin lymphoma at 6q22.33 (rs9482849, P = 1.52 × 10-8) and for nodular sclerosis Hodgkin lymphoma at 3q28 (rs4459895, P = 9.43 × 10-17), 6q23.3 (rs6928977, P = 4.62 × 10-11), 10p14 (rs3781093, P = 9.49 × 10-13), 13q34 (rs112998813, P = 4.58 × 10-8) and 16p13.13 (rs34972832, P = 2.12 × 10-8). Additionally, independent loci within the HLA region are observed for nodular sclerosis Hodgkin lymphoma (rs9269081, HLA-DPB1*03:01, Val86 in HLA-DRB1) and mixed cellularity Hodgkin lymphoma (rs1633096, rs13196329, Val86 in HLA-DRB1). The new and established risk loci localise to areas of active chromatin and show an over-representation of transcription factor binding for determinants of B-cell development and immune response.

Correlation between human leukocyte antigen class II alleles and HAI titers detected post-influenza vaccination

Influenza is a major cause of morbidity and mortality. Despite vaccination, many elderly recipients do not develop a protective antibody response. To determine whether Human Leukocyte Antigen (HLA) alleles modulate seroprotection to influenza, a cohort of HLA class II-typed high-risk vaccine recipients was investigated. Haemagglutinin inhibition (HAI) titres were measured 14–40 days post-subunit vaccination. Seroprotection was defined as HAI titres reaching 40 or greater for all three vaccine strains. HLA-DRB1*04∶01 and HLA-DPB1*04∶01 alleles were detected at higher frequencies in seroprotected compared with non-seroprotected individuals. Thus, the presence of certain HLA class II alleles may determine the magnitude of antibody responses to influenza vaccination.

Association of HLA-DPB1 with scleroderma and its clinical features in Chinese population

Human leukocyte antigen DPB1 was reported to contain singly nucleotide polymorphisms conferring the strongest susceptibility to systemic sclerosis in Korean population. However, associations of specific DPB1 alleles with SSc vary in different ethnic populations. The aim of this study was to profile DPB1 alleles in Chinese population and to identify specific DPB1 alleles in association with SSc and clinical and serological features of SSc in Han Chinese. A cohort containing 338 patients with SSc and 480 gender-matched and unrelated controls were examined in the study. The HLA-DPB1 genotyping was performed with sequence-based typing method. Exact p-values were obtained (Fisher's test) from 2×2 tables of allele counts or allele carriers and disease status. Thirty eight DPB1 alleles were found in the cohort. DPB1*05:01 was the most common allele in this cohort. DPB1*03:01 and *13:01 were significantly increased in SSc. DPB1*13:01 association had already been described in other ethnic populations, whereas DPB1*03:01 was specific to Han Chinese patients with SSc. In addition, comparisons between SSc subsets indicated that patients carrying DPB1*03:01 were more likely to develop pulmonary fibrosis, DPB1*04 carriers were increased in SSc patients with anti-centromere autoantibodies and in contrast, SSc patients with homozygous DPB1*05:01 showed an opposite association with marginal significance.

HLA-DPB1*04:01 Protects Genetically Susceptible Children from Celiac Disease Autoimmunity in the TEDDY Study

OBJECTIVES

Tissue transglutaminase autoantibodies (tTGAs) represent the first evidence of celiac disease (CD) development. Associations of HLA-DR3-DQA1*05:01-DQB1*02:01 (i.e., DR3-DQ2) and, to a lesser extent, DR4-DQA1*03:01-DQB1*03:02 (i.e., DR4-DQ8) with the risk of CD differ by country, consistent with additional genetic heterogeneity that further refines risk. Therefore, we examined human leukocyte antigen (HLA) factors other than DR3-DQ2 for their contribution to developing tTGAs.

METHODS

The Environmental Determinants of Diabetes in the Young (TEDDY) study enrolled 8,676 infants at an increased HLA-DR-DQ risk for type 1 diabetes and CD into a 15-year prospective surveillance follow-up. Of those followed up, 21% (n=1,813) carried DR3-DQ2/DR3-DQ2, 39% (n=3,359) carried DR3-DQ2/DR4-DQ8, 20% (n=1701) carried DR4-DQ8/DR4-DQ8, and 17% (n=1,493) carried DR4-DQ8/DQ4. Within TEDDY, a nested case-control design of 248 children with CD autoimmunity (CDA) and 248 matched control children were genotyped for HLA-B, -DRB3, -DRB4, -DPA1, and -DPB1 genes, and the entire cohort was genotyped for single-nucleotide polymorphisms (SNPs) using the Illumina ImmunoChip. CDA was defined as a positive tTGA test at two consecutive clinic visits, whereas matching in those with no evidence of tTGAs was based on the presence of HLA-DQ2, country, and sex.

RESULTS

After adjustment for DR3-DQ2 and restriction to allele frequency (AF) ≥5%, HLA-DPB1*04:01 was inversely associated with CDA by conditional logistic regression (AF=44%, odds ratio=0.71, 95% confidence interval (CI)=0.53-0.96, P=0.025). This association of time to CDA and HLA-DPB1*04:01 was replicated with statistical significance in the remainder of the cohort using imputation for specific HLA alleles based on SNP genotyping (hazard ratio=0.84, 95% CI=0.73-0.96, P=0.013).

CONCLUSIONS

HLA-DPB1*04:01 may reduce the risk of tTGAs, an early marker of CD, among DR3-DQ2 children, confirming that additional variants in the HLA region influence the risk for CDA.

Genetic association of human leukocyte antigens with chronicity or resolution of hepatitis B infection in thai population

BACKGROUND

Previous studies showed that single nucleotide polymorphisms (SNPs) in the HLA-DP, TCF19 and EHMT2 genes may affect the chronic hepatitis B (CHB). To predict the degree of risk for chronicity of HBV, this study determined associations with these SNPs.

METHODS

The participants for this study were defined into 4 groups; HCC (n = 230), CHB (n = 219), resolved HBV infection (n = 113) and HBV uninfected subjects (n = 123). The HLA-DP SNPs (rs3077, rs9277378 and rs3128917), TCF19 SNP (rs1419881) and EHMT2 SNP (rs652888) were genotyped.

RESULTS

Due to similar distribution of genotype frequencies in HCC and CHB, we combined these two groups (HBV carriers). The genotype distribution in HBV carriers relative to those who resolved HBV showed that rs3077 and rs9277378 were significantly associated with protective effects against CHB in minor dominant model (OR = 0.45, p<0.001 and OR = 0.47, p<0.001). The other SNPs rs3128917, rs1419881 and rs652888 were not associated with HBV carriers.

CONCLUSIONS

Genetic variations of rs3077 and rs9277378, but not rs3128917, rs1419881 and rs652888, were significantly associated with HBV carriers relative to resolved HBV in Thai population.