Abstract:Psoriasis is a complex, multigenic, immune-mediated skin disease characterized by heterogeneity across diverse ancestries; its prevalence varies from 0.09% to 11.4%, depending on the population of origin (Gibbs et al., 1996Gibbs S. Skin disease and socioeconomic conditions in rural Africa: Tanzania.Int J Dermatol. 1996; 35: 633-639Crossref PubMed Scopus (137) Google Scholar, Danielsen et al., 2013Danielsen K. Olsen A.O. Wilsgaard T. Furberg A.S. Is the prevalence of psoriasis increasing? A 30-year follow-up of a population-based cohort.Br J Dermatol. 2013; 168: 1303-1310Crossref PubMed Scopus (159) Google Scholar). In China, it is estimated that 0.47% of the population suffering from psoriasis (Ding et al., 2012Ding X. Wang T. Shen Y. Wang X. Zhou C. Tian S. et al.Prevalence of psoriasis in China: a population-based study in six cities.Eur J Dermatol. 2012; 22: 663-667Crossref PubMed Scopus (154) Google Scholar), and the prevalence of psoriasis in northern China is higher than in southern China. Recent genome-wide association studies have identified more than 80 regions associated with psoriasis, and a large transethnic genome-wide association study meta-analysis showed different genetic architecture of psoriasis susceptibility among populations (Tsoi et al., 2017Tsoi L.C. Stuart P.E. Tian C. Gudjonsson J.E. Das S. Zawistowski M. et al.Large scale meta-analysis characterizes genetic architecture for common psoriasis associated variants.Nat Commun. 2017; 8: 15382Crossref PubMed Scopus (162) Google Scholar, Yin et al., 2015Yin X. Low H.Q. Wang L. Li Y. Ellinghaus E. Han J. et al.Genome-wide meta-analysis identifies multiple novel associations and ethnic heterogeneity of psoriasis susceptibility.Nat Commun. 2015; 6: 6916Crossref PubMed Scopus (127) Google Scholar). Studies on the genetic structure of population showed one-dimensional north-south stratification in Han Chinese and showed large differences in allele frequency between north and south Han Chinese populations (Chen et al., 2009Chen J. Zheng H. Bei J.-X. Sun L. Jia W-h Li T. et al.Genetic structure of the Han Chinese population revealed by genome-wide SNP variation.Am J Hum Genet. 2009; 85: 775-785Abstract Full Text Full Text PDF PubMed Scopus (248) Google Scholar, Xu et al., 2009Xu S. Yin X. Li S. Jin W. Lou H. Yang L. et al.Genomic dissection of population substructure of Han Chinese and its implication in association studies.Am J Hum Genet. 2009; 85: 762-774Abstract Full Text Full Text PDF PubMed Scopus (273) Google Scholar, Zhou et al., 2016Zhou F. Cao H. Zuo X. Zhang T. Zhang X. Liu X. et al.Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease.Nat Genet. 2016; 48: 740-746Crossref PubMed Scopus (67) Google Scholar). Genetic structural difference across two subpopulations may be the result of evolutionary adaptations, as a consequence of environmental influences, including the effects of seasonality and climate, agricultural distribution, or varying prevalence of infectious diseases (Suo et al., 2012Suo C. Xu H. Khor C.C. Ong R.T. Sim X. Chen J. et al.Natural positive selection and north-south genetic diversity in East Asia.Eur J Hum Genet. 2012; 20: 102-110Crossref PubMed Scopus (37) Google Scholar). To further uncover the distinct genetic architecture of psoriasis across two populations in China, we included a total of 8,150 psoriasis patients and 9,906 matched control individuals and classified them into north Han Chinese ancestry (3,657 cases vs. 5,257 controls) and south Han Chinese ancestry (4,493 cases vs. 4,649 controls). Major histocompatibility complex (MHC) targeted sequencing and stringent quality control have been described elsewhere (Zhou et al., 2016Zhou F. Cao H. Zuo X. Zhang T. Zhang X. Liu X. et al.Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease.Nat Genet. 2016; 48: 740-746Crossref PubMed Scopus (67) Google Scholar). Then, we calculated the heritabilities of psoriasis in northern China (PsN) and psoriasis in southern China (PsS) using the genome-wide complex trait analysis method (Yang et al., 2011Yang J. Lee S.H. Goddard M.E. Visscher P.M. GCTA: a tool for genome-wide complex trait analysis..Am J Hum Genet. 2011; 88: 76-82Abstract Full Text Full Text PDF PubMed Scopus (0) Google Scholar) and performed a fine-mapping study to identify region-specific HLA associations in psoriasis. We observed that HLA alleles collectively explained a higher heritability of PsN (20.2%) than PsS (17.1%), which implies that there is different genetic architecture between PsN and PsS susceptibility. To further explain this result, we performed a stratification analysis to search for psoriatic region-specific risk alleles using stepwise analysis. We identified HLA-C*06:02 (P = 1.00 × 10–604,), HLA-C*07:04 (P = 3.12 × 10–25), B-Y67C (P = 1.21 × 10–23), B-Y116S (P = 7.71 × 10–18), and HLA-DPB1*05:01 (P = 1.21 × 10–10) were independently associated with PsN (see Supplementary Figure S1 and Supplementary Table S1 online). HLA-C*06:02 (P = 1.00 × 10–604), HLA-C*07:04 (P = 9.55 × 10–19), HLA-C*01:02 (P = 2.45 × 10–7), B-Y67C (P = 3.00 × 10–39), B-Q32L (P = 1.62 × 10–12), HLA-A*02:07 (corresponding to A-Y99C, P = 2.49 × 10–11), and DPB1-F35L (P = 2.30 × 10–7) were independently associated with PsS (see Supplementary Figure S2 and Supplementary Table S2 online). Most of these loci were consistent with those of our previous association study in the overall Han Chinese population (Zhou et al., 2016Zhou F. Cao H. Zuo X. Zhang T. Zhang X. Liu X. et al.Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease.Nat Genet. 2016; 48: 740-746Crossref PubMed Scopus (67) Google Scholar). HLA-C*06:02, HLA-C*07:04, and B-Y67C are shared, independent of psoriasis risk variants, across the two subpopulations. Moreover, HLA-C*06:02 and HLA-C*07:04 showed similarities in peptide binding preference (Pavlos et al., 2017Pavlos R. McKinnon E.J. Ostrov D.A. Peters B. Buus S. Koelle D. et al.Shared peptide binding of HLA class I and II alleles associate with cutaneous nevirapine hypersensitivity and identify novel risk alleles.Sci Rep. 2017; 7: 8653Google Scholar), and B-Y67C was located at functional pockets of the HLA-B molecule (Zhou et al., 2016Zhou F. Cao H. Zuo X. Zhang T. Zhang X. Liu X. et al.Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease.Nat Genet. 2016; 48: 740-746Crossref PubMed Scopus (67) Google Scholar), which suggest that these markers confer susceptibility to psoriasis in the Han Chinese population, because they are related to the disease itself and not to geographic location. However, other independent HLA variants (HLA-C*01:02, HLA-A*02:07, B-Y116S, and B-Q32L) showed region-specific associations with psoriasis. To further confirm the region-specific susceptibility loci, we directly comparing the PsN group with the PsS group and validated two independent region-specific risk loci for HLA-C*01:02 (Figure 1, Table 1) (P = 3.24 × 10–19) and HLA-A*02:07 (P = 4.11 × 10–8). HLA-C*01:02 and HLA-A*02:07 were independently associated with PsS (P = 2.45 × 10–7 and P = 2.49 × 10–11, respectively) and not with PsN (P = 0.003 and P = 0.024, respectively). We also identified that A-L180X and A-F33T were significantly associated with geographic locations, but they were not associated with psoriasis. DPB1-F35, a synonymous variant with high linkage to DPB1-F35L (r2 = 0.83) and DBP1*05:01 (r2 = 0.63), was found to meet the significance threshold in PsN-PsS analysis. The frequencies among these three variations were significantly higher in southern than in northern Han Chinese (see Supplementary Table S3 online). Therefore, we speculate that a signal marker within HLA-DPB1 is independently associated with psoriasis and contributes a higher risk to the PsS than the PsN population.Table 1Association of HLA variants with geographic location in Chinese psoriasis populationHLA VariantFrequencyPsN vs PsSFreq_PsNFreq_PsSOR (95% CI)PHLA-C AllelesHLA-C*01:020.0820.1270.61 (0.54–0.68)3.24 × 10–19HLA-A AllelesA-L156X0.2230.2740.77 (0.71–0.83)1.06 × 10–11A-F9T0.0740.0880.71 (0.63–0.81)6.74 × 10–8A-Y99C0.0530.0830.6546 (0.56–0.76)4.11 × 10–8HLA-DPB1 AllelesDPB1-F350.4260.4840.84 (0.79–0.90)6.28 × 10–7Abbreviations: CI, confidence interval; Freq_PsN, allele frequency in psoriatic patients from northern China; Freq_PsS: Allele frequency in psoriatic patients from southern China; OR, odds ratio. Open table in a new tab Abbreviations: CI, confidence interval; Freq_PsN, allele frequency in psoriatic patients from northern China; Freq_PsS: Allele frequency in psoriatic patients from southern China; OR, odds ratio. It makes sense that C*01:02 and A*02:07 showed PsS-specific risk loci, because their frequencies are much higher in southern Han Chinese than those in northern Han Chinese (see Supplementary Table S3). In addition, HLA-A*02 and HLA-C*01 have very similar peptide binding motifs (Andersen et al., 1999Andersen M.H. Sondergaard I. Zeuthen J. Elliott T. Haurum J.S. An assay for peptide binding to HLA-Cw∗0102.Tissue Antigens. 1999; 54: 185-190Crossref PubMed Scopus (20) Google Scholar), suggesting that specific autoantigen peptides presented by HLA-A*02:07 and HLA-C*01:02 could share a common structure and be involved in the pathogenesis of PsS. HLA-A*02:07 was reported to be associated with psoriasis in Japanese, and A-Y99C was located at functional pockets of the HLA-A molecule (Hirata et al., 2018Hirata J. Hirota T. Ozeki T. Kanai M. Sudo T. Tanaka T. et al.Variants at HLA-A, HLA-C, and HLA-DQB1 confer risk of psoriasis vulgaris in Japanese.J Invest Dermatol. 2018; 138: 542-548Abstract Full Text Full Text PDF PubMed Scopus (31) Google Scholar). Our previous study, in which samples were enrolled mainly from individuals in southern China, also showed that HLA-A*02:07 was associated with psoriasis in Han Chinese (Yin et al., 2015Yin X. Low H.Q. Wang L. Li Y. Ellinghaus E. Han J. et al.Genome-wide meta-analysis identifies multiple novel associations and ethnic heterogeneity of psoriasis susceptibility.Nat Commun. 2015; 6: 6916Crossref PubMed Scopus (127) Google Scholar). Moreover, we noted that the A*02:07-C*01:02-B*46:01 haplotype with a frequency of greater than 5% was reported in Han Chinese from the Zhejiang and Yunnan provinces, located in southern China (Chen et al., 2018Chen N. Wang W. Wang F. Dong L. Zhao S. Zhang W. et al.The distributions of HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1 allele and haplotype at high-resolution level in Zhejiang Han population of China.Int J Immunogenet. 2018; 45: 201-209Google Scholar, Shen et al., 2014Shen Y. Cao D. Li Y. Kulski J.K. Shi L. Jiang H. et al.Distribution of HLA-A, -B, and -C alleles and HLA/KIR combinations in Han population in China.J Immunol Res. 2014; 2014565296Crossref PubMed Scopus (14) Google Scholar). In our data set, the frequency of A*02:07-C*01:02-B*46:01 in healthy southern Han Chinese was significantly higher than in healthy northern Han Chinese (5.2% and 3.1%, respectively; P = 3.85 × 10–7). We did not identify any independent PsN-specific risk variants in HLA, despite a higher heritability. Possible reasons for this might be as follows. First, considering the statistical power, we removed variants with a minor allele frequency (MAF) of less than 1% for association analysis, although we used all sequencing data to calculate the heritability of psoriasis. There is a possibility that very-low-frequency variants (MAF < 1%) that confer risk for PsN were missed. Second, we did not have sufficient power (<80%) to detect low-frequency variants (MAF < 5%) with odds ratios of less than 1.6. Therefore, variants with MAF between 1% and 5% have a small effect and were also lost. In addition, statistical power of association analysis in PsN was slightly lower than PsS because of relatively small sample size and unbalanced distribution for northern Han Chinese. Large studies with sequencing data would clarify this in the future. All participating individuals provided written informed consent and were recruited according to the protocols approved by the institutional ethics committee of the First Affiliated Hospital of Anhui Medical University. The raw sequencing data from samples evaluated in the Han-MHC project and Psoriasis Association Analysis have been deposited in the Sequence Read Archive with the accession number SRA205317. The constructed Han-MHC database is available through the GigaScience database at https://doi.org/10.5524/100156. Minglong Cai: http://orcid.org/0000-0002-5038-3703 He Huang: http://orcid.org/0000-0003-2195-3209 Delin Ran: http://orcid.org/0000-0002-0153-0930 Xiaodong Zheng: http://orcid.org/0000-0002-6658-6492 Leilei Wen: http://orcid.org/0000-0002-8943-1549 Zhengwei Zhu: http://orcid.org/0000-0001-9343-8287 Lu Liu: http://orcid.org/0000-0002-4997-3415 Chuanliang Zhang: http://orcid.org/0000-0002-9246-7462 Xiaojie Hong: http://orcid.org/0000-0002-1746-0500 Jiaqi Hong: http://orcid.org/0000-0002-6877-2838 Wenjuan Wu: http://orcid.org/0000-0002-9436-474X Jie Ma: http://orcid.org/0000-0002-9412-0438 Mingshun Wu: http://orcid.org/0000-0002-9370-2890 Danfeng Qian: http://orcid.org/0000-0001-6252-4406 Yujun Sheng: http://orcid.org/0000-0002-6022-3359 The authors state no conflict of interest. This study was supported financially by the National Natural Science Foundation of China (81872527, 81130031, 81803117) and the Natural Science Foundation of Anhui Province (1808085QH284). We are grateful to all participants, their families, and healthy control individuals who donated blood samples in this study. MC: lead: writing (original draft preparation), visualization; equal: conceptualization, data curation, formal analysis; supporting: writing (review and editing). HH: lead: writing (review and editing); equal: conceptualization, data curation; supporting: writing (original draft preparation), DR: equal: conceptualization, data curation. XZ: equal: data curation, formal analysis. LW: equal: data curation; supporting: writing (review and editing). ZZ: equal: data curation; supporting: visualization. LL: equal: conceptualization; supporting: data curation. CZ: equal: data curation; supporting: writing (review and editing. XH: equal: writing (review and editing); supporting: visualization. JH: equal: writing (review and editing); supporting: visualization. WW: equal: formal analysis; supporting: writing (review and editing). JM: equal: conceptualization; supporting: data curation. MW: equal: data curation; supporting: visualization. DQ: equal: data curation; supporting: methodology. YS: lead: conceptualization, methodology; equal: funding acquisition, supervision, writing (original draft preparation), writing (review and editing. XZ: lead: funding acquisition; equal: conceptualization, methodology, supervision, writing (original draft preparation), writing (review and editing). All authors read and approved the final manuscript. A total of 8,150 psoriasis patients (patients with psoriatic arthritis were not included) and 9,906 matched control individuals of Han Chinese ethnicity were recruited through collaboration with multiple medical institutions in China. All selected patients had disease diagnosed by at least two dermatologists, and all control individuals were healthy, without psoriasis, autoimmune diseases, systemic diseases, or a family history of genetic diseases (including first-, second- and third-degree relatives). We classified participants into northern Han Chinese ancestry (3,657 patients and 5,257 control individuals) and southern Han Chinese ancestry (4,493 patients and 4,649 control individuals), using the Huai River–Qin Mountains as the dividing line. All participating individuals provided written informed consent and were recruited according to the protocols approved by the institutional ethics committee of the First Affiliated Hospital of Anhui Medical University. The samples enrolled for analysis including 8,150 patient samples and 9,906 control samples for which HLA sequences were generated and stringent quality control were conducted in our previous MHC targeted sequencing research (Zhou et al., 2016Zhou F. Cao H. Zuo X. Zhang T. Zhang X. Liu X. et al.Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease.Nat Genet. 2016; 48: 740-746Crossref PubMed Scopus (67) Google Scholar). In addition to this, we excluded variants with MAF less than 1%, Hardy-Weinberg equilibrium of P less than 10–4, or with a call rate of less than 99%. Finally, we obtained a total of 26,775 markers including single nucleotide polymorphisms, classic alleles, and amino acid residues. The statistical power of our study was evaluated by Power and Sample Size software, version 3.1.2 (www.power-analysis.com/). Type I error was set to 1.87 × 10–6 (after Bonferroni correction, P = 0.05/26,775). We had a power of greater than 80% for detecting MAF greater than 2.5%, and odds ratio of the effect allele was equal to 2, but we did not have sufficient power (<80%) to detect low-frequency variants (MAF <0.5%) even if the odds ratio was greater than 3.5. Using genome-wide complex trait analysis (Yang et al., 2011Yang J. Lee S.H. Goddard M.E. Visscher P.M. GCTA: a tool for genome-wide complex trait analysis.Am J Hum Genet. 2011; 88: 76-82Abstract Full Text Full Text PDF PubMed Scopus (3814) Google Scholar, Yang et al., 2015Yang J. Bakshi A. Zhu Z. Hemani G. Vinkhuyzen A.A. Lee S.H. et al.Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index.Nat Genet. 2015; 47: 1114-1120Crossref PubMed Scopus (467) Google Scholar), we calculated the heritability of PsS and PsN conferred by MHC region with the assumption of a disease prevalence of 0.47%. Analysis of all variants were performed with logistic regression, assuming an additive effect of the allele dosages in the log-odds scale. We included sex, family history, and age at onset as covariates to correct for population stratification. HLA associations with three binary psoriasis susceptibility characteristics defined by region were analyzed as follows: i) associations between HLA alleles with psoriasis susceptibility in northern China (PsN vs control individuals from northern China), ii) associations between HLA alleles with psoriasis susceptibility in southern China (PsS vs control individuals from southern China), iii) directly comparing psoriasis patients in northern China with patients in southern China (PsN vs PsS). For stepwise conditional analysis, we continuously included the identified marker as a covariant to find additional associated variants with independent effect. This was repeated in a forward stepwise approach until no variation exceeded the predefined significance threshold after Bonferroni correction (P < 1.87 × 10–6). Logistic regression and stepwise conditional analysis were conducted by PLINK, version 1.9 (Chang et al., 2015Chang C.C. Chow C.C. Tellier L.C. Vattikuti S. Purcell S.M. Lee J.J. Second-generation PLINK: rising to the challenge of larger and richer datasets.Gigascience. 2015; 4: 7Crossref PubMed Scopus (4498) Google Scholar).Supplementary Figure S2A Manhattan plot of the stepwise conditional association for PsS in the MHC region. Stepwise analysis of HLA-C, HLA-B, HLA-A, and HLA-DPB1 in PsS versus control individuals. For each plot, the horizontal axis shows genomic position and the vertical axis shows negative log10-transformed P-values for association. The dashed horizontal line corresponds to the significance threshold of P = 1.87 × 10−6. MHC, major histocompatibility complex; PsS, psoriatic patients from southern China.View Large Image Figure ViewerDownload Hi-res image Download (PPT)Supplementary Table S1Association of HLA variants with psoriasis in northern ChinaHLA VariantFrequencyPsN vs NCFreq_PsNFreq_NCOR (95% CI)PHLA-C AllelesHLA-C*06:020.4540.12712.31 (11.11–13.63)1.00 × 10–501HLA-C*07:040.0190.0124.383 (3.32–5.80)3.12 × 10–25HLA-B AllelesB-Y67C0.0900.0862.04(1.77–2.34)1.21 × 10–23B-Y116S0.3010.2951.45 (1.33–1.58)7.71 × 10–18HLA-DPB1 AllelesHLA-DPB1*05:010.3430.3451.31 (1.20–1.42)1.21 × 10–10Abbreviations: CI, confidence interval; Freq_PsN, allele frequency in psoriatic patients from northern China; NC, control individuals from northern China; OR, odds ratio. Open table in a new tab Supplementary Table S2Association of HLA variants with psoriasis in southern ChinaHLA variantFrequencyPsS vs SCFreq_PsSFreq_SCOR (95% CI)PHLA-C AllelesHLA-C*06:020.4230.09216.52 (14.88–18.33)1.00 × 10–604HLA-C*07:040.0130.0084.24 (3.08–5.84)9.55 × 10–19HLA-C*01:020.1270.1681.32 (1.19–1.46)2.45 × 10–07HLA-B AllelesB-Y67C0.0950.0892.43 (2.13–2.77)3.00 × 10–39B-Q32L0.1340.2340.69 (0.62–0.76)1.62 × 10–12HLA-A AllelesHLA-A*02:071Equivalent to A-Y99C in the Han Chinese population.0.0830.0871.67 (1.44–1.94)2.49 × 10–11HLA-DPB1 AllelesDPB1-F35L0.4370.471.231 (1.14–1.33)2.30 × 10–07Abbreviations: CI, confidence interval; Freq_PsN, allele frequency in psoriatic patients from northern China; SC, control individuals from southern China; OR, odds ratio.1 Equivalent to A-Y99C in the Han Chinese population. Open table in a new tab Supplementary Table S3Allele frequencies between the northern Han Chinese and southern Han Chinese populationsHLA AllelesAllele Frequency in Northern ChineseAllele Frequency in Southern ChineseFreq_PsNFreq_NCFreq_PsSFreq_SCHLA-C*01:020.081620.11250.12710.1683HLA-C*06:020.45430.12710.42330.09163HLA-C*07:040.018460.01170.013240.008496A-F33T0.074240.12450.088250.1347A-Y123C0.052780.056970.082680.08658A-L180X0.22310.32450.2740.4099B-Y140S0.30120.29520.33080.3612B-Y91C0.089960.08560.09470.08948B-Q56L0.13850.22690.13350.2337BTNL2:T165I0.059890.064590.056660.06528DPB1-F640.4260.45560.48350.5228DPB1*05:010.3430.34470.37280.3864DPB1-F64L0.38980.4080.4370.47Bold type indicates alleles that showed region-specific risk for psoriasis in southern China.Abbreviations: Freq_NC, allele frequency in control individuals from northern China; Freq_PsN, allele frequency in psoriatic patients from northern China; Freq_PsS, allele frequency in psoriatic patients from southern China; Freq_SC, allele frequency in control individuals from southern China. Open table in a new tab Abbreviations: CI, confidence interval; Freq_PsN, allele frequency in psoriatic patients from northern China; NC, control individuals from northern China; OR, odds ratio. Abbreviations: CI, confidence interval; Freq_PsN, allele frequency in psoriatic patients from northern China; SC, control individuals from southern China; OR, odds ratio. Bold type indicates alleles that showed region-specific risk for psoriasis in southern China. Abbreviations: Freq_NC, allele frequency in control individuals from northern China; Freq_PsN, allele frequency in psoriatic patients from northern China; Freq_PsS, allele frequency in psoriatic patients from southern China; Freq_SC, allele frequency in control individuals from southern China.

Abstract:Dermatitis herpetiformis (DH) is a rare autoimmune blistering disease linked to gluten sensitivity; DH follows a chronically relapsing and remitting course, characterized by pruritic polymorphic lesions and granular deposition of IgA in the dermal papillary tips (Bolotin and Petronic-Rosic, 2011Bolotin D. Petronic-Rosic V. Dermatitis herpetiformis: part I. Epidemiology, pathogenesis, and clinical presentation.J Am Acad Dermatol. 2011; 64: 1017-1024Abstract Full Text Full Text PDF PubMed Scopus (120) Google Scholar). Currently, DH is considered as the specific cutaneous manifestation of celiac disease (CD) because of the presence of circulating IgA antitissue, epidermal transglutaminase antibodies, and the same typical histological features of small bowel villous atrophy (Reunala, 1998Reunala T. Dermatitis herpetiformis: coeliac disease of the skin.Ann Med. 1998; 30: 416-418Crossref PubMed Scopus (74) Google Scholar). The incidence and prevalence of DH varies across ethnic groups and is more common in Caucasians (Salmi et al., 2011Salmi T.T. Hervonen K. Kautiainen H. Collin P. Reunala T. Prevalence and incidence of dermatitis herpetiformis: a 40-year prospective study from Finland.Br J Dermatol. 2011; 165: 354-359Crossref PubMed Scopus (107) Google Scholar) but rare in Asians (Zhang et al., 2012Zhang F. Yang B. Lin Y. Chen S. Zhou G. Wang G. et al.Dermatitis herpetiformis in China: a report of 22 cases.J Eur Acad Dermatol Venereol. 2012; 26: 903-907Crossref PubMed Scopus (23) Google Scholar). For the original Chinese population, only a few sporadic cases have been reported in the population of mainland China (Zhang et al., 2012Zhang F. Yang B. Lin Y. Chen S. Zhou G. Wang G. et al.Dermatitis herpetiformis in China: a report of 22 cases.J Eur Acad Dermatol Venereol. 2012; 26: 903-907Crossref PubMed Scopus (23) Google Scholar), Hong Kong, and Singapore (Ratnam, 1988Ratnam K.V. IgA dermatosis in an adult Chinese population. A 10-year study of linear IgA and dermatitis herpetiformis in Singapore.Int J Dermatol. 1988; 27: 21-24Crossref PubMed Scopus (18) Google Scholar). The heritability of DH has not been well understood, but a genetic predisposition has been suggested by the reported associations of both HLA class I and class II molecules in various populations. Case-control studies in early stage patients with DH showed that DH is strongly associated with HLA-B8, DR3, and DQw2 haplotypes (Hall et al., 1989Hall R.P. Sanders M.E. Duquesnoy R.J. Katz S.I. Shaw S. Alterations in HLA-DP and HLA-DQ antigen frequency in patients with dermatitis herpetiformis.J Invest Dermatol. 1989; 93: 501-505Abstract Full Text PDF PubMed Scopus (41) Google Scholar, Hall et al., 1996Hall M.A. Lanchbury J.S. Ciclitira P.J. HLA class II region genes and susceptibility to dermatitis herpetiformis: DPB1 and TAP2 associations are secondary to those of the DQ subregion.Eur J Immunogenet. 1996; 23: 285-296Crossref PubMed Scopus (11) Google Scholar, Katz, 1977Katz S.I. HLA-B8 and dermatitis herpetiformis in patients with IgA deposits in skin.Arch Dermatol. 1977; 113: 155-156Crossref PubMed Scopus (39) Google Scholar). As in CD in Caucasians, more than 90% of patients share the HLA-DQ2 haplotype, with the majority of the remainder expressing HLA-DQ8 (Kaukinen et al., 2002Kaukinen K. Partanen J. Mäki M. Collin P. HLA-DQ typing in the diagnosis of celiac disease.Am J Gastroenterol. 2002; 97: 695-699Crossref PubMed Google Scholar). NOD DQ8+ murine models have subsequently confirmed these associations (Marietta et al., 2004Marietta E. Black K. Camilleri M. Krause P. Rogers R.S. David C. et al.A new model for dermatitis herpetiformis that uses HLA-DQ8 transgenic NOD mice.J Clin Invest. 2004; 114: 1090-1097Crossref PubMed Scopus (124) Google Scholar). By contrast, in Asians, literature reviews of previous studies (Ohata et al., 2012Ohata C. Ishii N. Hamada T. Shimomura Y. Niizeki H. Dainichi T. et al.Distinct characteristics in Japanese dermatitis herpetiformis: a review of all 91 Japanese patients over the last 35 years.Clin Dev Immunol. 2012; 2012: 1-9Crossref Scopus (43) Google Scholar) have revealed distinct characteristics of DH in Asian from Caucasian, which may be due to the absence of the HLA-B8/DR3/DQ2 haplotype in Japan (Shibahara et al., 2002Shibahara M. Nanko H. Shimizu M. Kanda N. Kubo M. Ikeda M. et al.Dermatitis herpetiformis in Japan: an update.Dermatology. 2002; 204: 37-42Crossref PubMed Scopus (37) Google Scholar). To our knowledge, there has been no genetic investigation of Chinese patients with DH published, nor any relative HLA association reported. To establish genetic risk factors for DH in the Chinese population, we conducted a comprehensive genetic study consisting of a two-stage genome-wide variant screen and a next-generation sequence-based HLA typing of DH in the Chinese population. The discovery analysis included 24 patients with DH and 486 controls who were all recruited from the Shandong Provincial Institute of Dermatology and Venereology in China from 2006 (Supplementary Table S1 online). An additional 12 patients with DH and 244 controls were recruited in the same way for validation analysis. The clinical diagnostic criteria and the serological test of patients with DH were described in the Supplementary Methods online. The controls for both the discovery and validation analyses were clinically confirmed to lack any autoimmune diseases. The study was approved by the institutional review board at the Shandong Provincial Institute of Dermatology and Venereology. All the subjects gave written informed consent. In our patients, men seemed to be more frequently affected than women with a male-to-female ratio of 2.27:1. The mean age of DH onset was in the fourth decade of the patients and the mean duration of symptoms was approximately 44.9 months. Twenty-five percent of patients experienced gastrointestinal symptoms, mostly comprising chronic diarrhea, steatorrhea, weight loss, and nonspecific abdominal discomfort. Among these patients, only two patients received a duodenal biopsy and both specimens showed mild-to-moderate villous atrophy on histological examination (Supplementary Figure S1 online). For the serological data (Supplementary Table S1), 11 of 25 patients (44%) were found to be anti-endomysial antibody positive, with titers ranging from 1:10 to 1:320; and 10 of 17 patients (58.8%) were found to be anti-tissue transglutaminase antibody positive (cutoff value > 20 RU/ml), of whom 8 patients also showed anti-endomysial antibody positivity. For each participant, genomic DNA was extracted for genotyping analysis. We genotyped all 24 cases and 486 controls of the discovery analysis using Illumina Omni Zhonghua chips with 900,015 single nucleotide polymorphisms (SNPs). After standard sample and SNP QC filtering (Zhang et al., 2011Zhang F. Liu H. Chen S. Low H. Sun L. Cui Y. et al.Identification of two new loci at IL23R and RAB32 that influence susceptibility to leprosy.Nat Genet. 2011; 43: 1247-1251Crossref PubMed Scopus (136) Google Scholar), a total of 24 cases and 480 controls with 806,342 SNPs were included for genome-wide association analyses. Rs2187668, located at the intronic site of HLA-DQA1, was identified as the most significant association. Furthermore, a group of SNPs within the major histocompatibility complex region also showed significant associations (P < 1.0 × 10−5). To clarify the extensive associations observed within the major histocompatibility complex region, we imputed the genotyping SNPs to classical HLA alleles and corresponding amino acid polymorphisms using the reference panel published previously (Okada et al., 2014Okada Y. Kim K. Han B. Pillai N.E. Ong R.T. Saw W.Y. et al.Risk for ACPA-positive rheumatoid arthritis is driven by shared HLA amino acid polymorphisms in Asian and European populations.Hum Mol Genet. 2014; 23: 6916-6926Crossref PubMed Scopus (111) Google Scholar), which contains genotype data for 6,173 SNPs and classical alleles of HLA-A, B, C, DRB1, DQA1, DQB1, DPA1, and DPB1 at 4-digit resolution in 530 unrelated individuals of Asian descent. Two strong associations were found for HLA-B*0801 and HLA-DRB1*0301, whereas conditioning for either could not eliminate the effect of the other one. In addition, no further associations (P < 1.0 × 10−5) were discovered after conditioning for the effects of HLA-DRB1*0301 and HLA-B*0801. To validate the results of the imputation analysis, we sequenced HLA-DRB1 and HLA-B in all discovery samples as well as an additional 12 DH cases and 244 controls by using long-range PCR with previously published primers (Bentley et al., 2009Bentley G. Higuchi R. Hoglund B. Goodridge D. Sayer D. Trachtenberg E.A. et al.High-resolution, high-throughput HLA genotyping by next-generation sequencing.Tissue Antigens. 2009; 74: 393-403Crossref PubMed Scopus (157) Google Scholar). Individual libraries of each amplicon were prepared by enzymatic fragmentation, end repair, adenylation, and ligation of indexed adaptors. The 12 indexed libraries were pooled before paired-end 2 × 250 bp sequencing on an illumina Miseq platform. Sequence reads were aligned and analyzed to the HLA region using Omixon target software (www.omixon.com) with default settings. A mean coverage of ×114 was identified for exon 2 of HLA-DRB1 and exon 2+3 of HLA-B per sample across all the samples. A direct comparison between the results of imputation and sequencing analyses revealed a high concordance of 87.31% at HLA-DRB1 and 92.57% at HLA-B for the 4-digit resolution. Analysis of the validation sample and the joint analysis of sequenced samples confirmed the strong association of HLA-B*0801 (P = 6.87 × 10−15, OR = 28.67) and HLA-DRB1*0301 (P = 2.1 × 10−13, OR = 9.7) with DH (Table 1). In the patients with DH, there was a significant increase in the HLA-B*0801 allele frequency (Freq = 18.1%) and in HLA-DRB1*0301 allele frequency (25%) compared with the healthy population. Conditioning on HLA-B*0801 could not eliminate the effect of HLA-DRB1*0301 (PconB*0801 < 0.05), and vice versa.Table 1Results of association analyses for significant classical HLA allelesClassical HLA alleleDiscovery setValidation setJoint analysesF_AF_UPORF_AF_UPORF_AF_UPORHLA-B*08010.1670.0077.33E-0827.230.2080.0081.77E-0531.320.1810.0086.87E-1528.67HLA-DRB1*03010.2500.0353.92E-079.080.2500.0291.48E-0411.140.2500.0332.10E-139.70F_A represents the frequency of present alleles in cases, whereas F_U represents the frequency of present alleles in controls. P and OR, P value and corresponding odds ratio per present allele, respectively. The P values and odds ratios for the joint analysis were adjusted by the study variable. Open table in a new tab F_A represents the frequency of present alleles in cases, whereas F_U represents the frequency of present alleles in controls. P and OR, P value and corresponding odds ratio per present allele, respectively. The P values and odds ratios for the joint analysis were adjusted by the study variable. As 25% of the patients with DH suffered from gastrointestinal symptoms, we performed subgroup analysis of the patients with and without gastrointestinal symptoms on the two identified risk alleles. Interestingly, we found that HLA-DRB1*0301 seemed to be a protective factor for patients with DH against gastrointestinal symptoms with nominal significance (P = 0.03, OR = 0.128, Supplementary Table S2 online). Finally, we evaluated the effectiveness of these two risk variants (B*0801 and DRB1*0301) for risk prediction using all the sequenced samples (Figure 1). B*0801 was present in 36.1% of the patients with DH (13 of 36), but only 1.5% of the healthy controls (11 of 720). Meanwhile, the presence of DRB1*0301 has a sensitivity of 50% and specificity of 93.3% (with an area under the curve value of 0.717) for risk prediction. Inclusion of these two risk variants (presence of either) could improve the performance of risk predictions but in a moderate way, increasing the sensitivity to 61.1% and the area under the curve value to 0.776 and decreasing the specificity to 93.1%. Based on the estimated prevalence of 1/100,000, the negative predictive value and positive predictive value were almost the same when employing the two risk variants (negative predictive value = 99.9%, positive predictive value = 0.009%) compared with DRB1*14 alone (negative predictive value = 99.9%, positive predictive value = 0.008%) or HLA-B*0801 alone (negative predictive value = 99.9%, positive predictive value = 0.024%) for predicting DH. Human HLA genes are extremely polymorphic, with specific alleles determined by complex gene sequence polymorphisms with distinct affinity for antigenic peptides for an effective immune response. In our study, we confirmed two independent classical HLA alleles associated with DH, of which HLA-DRB1*0301 was associated with the avoidance of gastrointestinal symptoms in patients with DH. Although DH was currently considered as the specific cutaneous manifestation of CD, only a very small number of patients with DH were considered to have CD because only 25% of patients with DH had gastrointestinal symptoms in our data and the rarity of CD incidence in China. Epidemiologically, HLA-B*0801 and HLA-DRB1*0301 are rare in Chinese populations (Qin et al., 2011Qin P.Q. Su F. Yan W.X. Xing Z. Meng P. Chengya W. et al.Distribution of human leucocyte antigen-A, -B and -DR alleles and haplotypes at high resolution in the population from Jiangsu province of China.Int J Immunogenet. 2011; 38: 475-481Crossref PubMed Scopus (32) Google Scholar) (B*0801 = 0.88%, DRB1*0301 = 4.9%), whereas the frequency of these alleles is much higher in the Caucasians (B*0801 = 11.44%, DRB1*0301 = 12.16%). Previously reported risk alleles in the Caucasians, such as HLA-DQ2 and DQ8, showed highly genetic heterogeneity between Caucasian and Chinese populations (Supplementary Table S3 online). In our data, the most significant signals encoding the DQ2 group are DQB1*0201 (29.2% vs. 11%, P = 7.64 × 10−4) and DQA1*0501 (27.1% vs. 14.1%, P = 0.02), whereas genes encoding DQ8 showed no significance (DQA1*0301 [25% vs. 25.8%, P = 1], DQB1*0302 [6.25% vs. 7.6%, P = 1]). This might be explained by the population heterogeneity between Caucasian and Chinese populations. Although many studies have emphasized that genetic factors other than HLA play an important role in the pathogenesis of DH, we have not discovered any significant associations outside of the major histocompatibility complex region. This finding strongly supports that the pathogenic mechanism of DH is highly driven by an HLA-mediated autoimmune response. The authors state no conflicts of interest. We thank the study participants and research staffs (Paul I.W. de Bakker at University Medical Center Utrecht, Yukinori Okada at Tokyo Medical and Dental University, Soumya Raychaudhuri at Brigham and Women’s Hospital) for their contributions on the HLA imputation analyses. We also thank the Shandong Computer Science Center (National Supercomputer Center in Jinan) for their contribution to the data processing. This work was funded by grants from the Shandong Provincial Advanced Taishan Scholar Construction Project, the Science and Technology Plan of Shandong Academy of Medical Science (2014-09, 201201), and the Health and Family Planning Commission of Shandong Province (2013WS0368, 2014WS0064). Download .pdf (2.46 MB) Help with pdf files Supplementary Data

Association Analysis of Allergic Sensitization Susceptibility Loci with Atopic Dermatitis in Chinese Population.

FOUR NOVEL SUSCEPTIBILITY VARIANTS IDENTIFIED IN CHINESE ARE ASSOCIATED WITH SPORADIC ALZHEIMER’S DISEASE

Association Analyses Identify Two Susceptibility Loci 5q31 and 5q22.1 for Atopic Dermatitis in Chinese Han Population.

Association Study on ADAM33 Polymorphisms in Mite-Sensitized Persistent Allergic Rhinitis in a Chinese Population.

Identification of Candidate Genes in Atopic Dermatitis Based on Bioinformatic Methods

A Common Variant in 11q23.3 Is Associated with Susceptibility to Atopic Dermatitis in the Han Chinese Population

Association Analysis of Single Nucleotide Polymorphisms at Five Loci: Comparison Between Atopic Dermatitis and Asthma in the Chinese Han Population

GWAS Follow-up Study Discovers a Novel Genetic Signal on 10q21.2 for Atopic Dermatitis in Chinese Han Population

A Variant on Chromosome 2p13.3 is Associated with Atopic Dermatitis in Chinese Han Population

Association of the Chromosome 11p13.5 Variant and Atopic Dermatitis with a Family History of Atopy in the Chinese Han Population.

Association Between the IL-10-1082G/A, IL-10-592A/C, and IL-10-819G/A Polymorphisms and Atopic Dermatitis Susceptibility: A Meta-Analysis.

Associations between interleukin-13, interleukin-4 and their receptor gene polymorphisms and susceptibility to atopic dermatitis in a Chinese Han population

HLA-C01:02 and HLA-A02:07 Confer Risk Specific for Psoriatic Patients in Southern China

Association between IL-6 polymorphisms and Atopic Dermatitis in Chinese Han children

Validating Gwas-Identified Risk Loci for Alzheimer'S Disease in Han Chinese Populations

Meta-analysis of Genome-Wide Linkage Studies of Atopic Dermatitis.

Association of IL-4 and IL-18 genetic polymorphisms with atopic dermatitis in Chinese children

Scanning Indels in the 5q22.1 Region and Identification of the TMEM232 Susceptibility Gene That is Associated with Atopic Dermatitis in the Chinese Han Population

The HLA Alleles B0801 and DRB10301 Are Associated with Dermatitis Herpetiformis in a Chinese Population

The Association Between Sixteen Genome-Wide Association Studies-Related Allergic Diseases Loci and Childhood Allergic Rhinitis in a Chinese Han Population.

Association of UBASH3A Gene Polymorphism and Atopic Dermatitis in the Chinese Han Population