Abstract

Psoriasis is a multifactorial immune-mediated skin disease with a strong genetic background. Previous studies reported that psoriasis with a family history (PFH) and sporadic psoriasis (SP) have a distinct manifestation and genetic predisposition. However, the genetic heterogeneity of PFH and SP in the major histocompatibility complex (MHC) region has not been fully elucidated. To explore genetic variants in the MHC region that drive family aggregation of psoriasis, we included a total of 8,127 psoriasis cases and 9,906 healthy controls from Han Chinese and divided psoriasis into two subtypes, PFH () and SP (). Then, we calculated the heritability of PFH and SP and performed a large-scale stratified association analysis. We confirmed that variants in the MHC region collectively explained a higher heritability of PFH (16.8%) than SP (13.3%). Further stratified association analysis illustrated that HLA-C06:02 and NOTCH4:G511S contribute to the family aggregation of psoriasis, and BTNL2:R281K specifically confers risk for SP. HLA-C06:02 and NOTCH4:G511S could partially explain why patients with PFH have a stronger genetic predisposition, more complex phenotypes, and more frequent other autoimmune diseases. The identification of the SP-specific variant BTNL2:R281K revealed that the genetic architecture of SP is not just a subset of PFH.

1. Introduction

Psoriasis is a chronic and multifactorial immune-mediated skin disease characterized by epidermal hyperplasia. It affects 0.09%–11.4% of worldwide population [1, 2]. Psoriasis has a strong genetic predisposition. Evidence shows monozygotic twins have higher concordance rates (0.2–0.73) than dizygotic twins (0.09–0.2) [38]. Moreover, 23.1%–31.9% of psoriasis patients have a family history [9, 10]. To date, linkage analysis and genome-wide association studies (GWAS) have identified more than 80 disease-associated loci, with the most prominent risk loci observed in the human leukocyte antigen (HLA) region [1114].

Although patients with psoriasis tend to cluster in families, most of them are sporadic. Evidence has shown psoriasis with a family history (PFH) is different from sporadic psoriasis (SP) in terms of both genetic background and phenotype. Genes outside of the major histocompatibility complex (MHC) region such as ERAP1 and NFKB1 were found to be associated with the risk of familial aggregation of psoriasis [15, 16]. A clinical investigation reported that PFH patients have an earlier onset of psoriasis and more frequently nail disease, enthesitis, and other autoimmune diseases [10, 17]. It seems that PFH has a stronger genetic background and more complicated phenotypes than SP.

However, the effects of the genetic heterogeneity of the HLA alleles on PFH and SP had been largely underexplored in spite of the strong associations between HLA alleles and psoriasis. One of the greatest challenges is that HLA genes are highly polymorphic and can be in strong linkage disequilibrium, which has complicated the determination of the independent risk signals. In this study, we devised an analytical approach to identify the heterogeneous effects of HLA alleles on PFH and SP by deep sequencing of the MHC region. We first calculated the heritability of PFH and SP and confirmed that variants in the MHC region collectively explained a higher heritability of PFH than SP. Then, we considered two hypothetical models of the genetic structure in PFH versus SP: (1) We considered that shared variants might be associated with both PFH and SP; such variants would confer more risk of PFH than of SP. (2) We considered that specific variants might specifically confer risk for PFH; such variants would not be associated with SP. Based on two hypothetical models, we conducted a step-wise association analysis of HLA with PFH and SP to identify those shared and specific variants for psoriasis subtypes. Finally, we performed a case-case analysis (PFH vs. SP) to determine whether those shared and specific variants contribute to the family aggregation of psoriasis.

2. Materials and Methods

2.1. Participants

We included a total of 8,127 patients with psoriasis (patients with psoriasis arthritis were not included) and 9,906 geographically matched healthy controls of Chinese ancestry, as previously described [18]. All cases were diagnosed by at least two dermatologists, and all controls were healthy individuals without psoriasis, any autoimmune diseases, or a family history of genetic diseases (including first-, second- and third-degree relatives). Clinical information for all individuals was gathered using a previously developed structured questionnaire. All participants signed informed consent. The study was approved by the institutional ethical committee of the First Affiliated Hospital of Anhui Medical University (No.20131349) and was conducted in accordance with the Declaration of Helsinki principles.

2.2. Sequencing and Quality Control

The samples enrolled for our analysis underwent targeted sequencing and stringent quality control as described in our previous study [18]. The raw sequencing data from samples have been deposited in the Sequence Read Archive (SRA) under the accession number SRA205317. Furthermore, we removed variants with a call , , or Hardy–Weinberg equilibrium of . A total of 26,775 variants remained for the following analysis.

2.3. Phenotype Classification

We classified psoriasis () into two subtypes, PFH () and SP (). To rule out the influence of the confounder of the age of onset, we removed patients with age of onset after 40 years, retaining 1,538 PFH and 5,262 SP for the following analysis.

2.4. Statistical Analyses

The association of age of onset with a family history of psoriasis was evaluated by the Pearson chi-square test and Student’s test. The association analysis of HLA with psoriasis subtypes was performed using logistic regression, assuming an additive effect of the allele dosages in the log-odds scale. We included sex and region as covariates to adjust for population stratification. Association analysis was conducted as follows: (1) association analysis of variants with PFH and (2) association analysis of variants with SP. For conditional association analysis, we consecutively included the detected variants as covariates to search for additional variants with independent association. This was repeated in a forward stepwise approach until no marker satisfied the significance threshold (Bonferroni correction: ). For those variants that were independently associated with PFH or SP, we further conducted a case-case association analysis by using logistic regression after adjusting for sex and region. For epistasis analysis, we defined a logistic regression model that simultaneously included interaction terms between significant variants. Association analysis and conditional association analysis were conducted by PLINK (v.1.9) [19]. Plots of HLA associations were drawn by R (v.3.6.1).

2.5. Annotation of Variants

For variants with amino acid changes, we used SIFT [20], PolyPhen-2 [21], CADD [22], and DANN [23] to measure whether an amino acid is deleterious based on changes in the structure/function of proteins, allelic diversity, pathogenicity, and conservation. The criteria used are as follows: SIFT: deleterious (), tolerated (); PolyPhen-2: probably damaging (), possibly damaging (); benign (); CADD: top 10% most deleterious (); top 1% most deleterious (); DANN: top 10% most deleterious (); top 1% most deleterious ().

3. Results

3.1. Sample Selection and Phenotype Classification

A total of 8,127 psoriasis patients were included. Clinically, psoriasis can be classified into type 1 psoriasis (defined as age of onset before 40) and type 2 psoriasis (defined as age of onset after 40) [24]. It is well known that genetic factors are more prevalent in type 1 psoriasis, and type 1 psoriasis is more likely to be familial [24]. In our study, we observed that type 1 psoriasis was strongly associated with familial aggregation (, ; Table 1). In order to rule out the influence of age of onset, we excluded patients with age of onset after 40, retaining 1,538 PFH (mean age of onset: ) and 5,262 SP (mean age of onset: ) cases for the following analysis. There was no substantial correlation between age of onset and family history () after removing patients with type 2 psoriasis.

3.2. Variants Identification

Our target sequencing of the entire 5-Mb MHC region, from upstream of HLA-A to downstream of HLA-DPB1, detected a total of 36,300 indels and 224,872 SNPs [18]. Next, we generated accurate genotyping results (i.e., HLA alleles, amino acid polymorphisms, etc.) for genes including the HLA classic genes (i.e., HLA-A, HLA-B, HLA-C, etc.) and the non-HLA genes (i.e., NOTCH4, TAP2, MICA, etc.). After stringent quality control as previously described [18], we further removed variants with a call , minor allele frequency , or Hardy–Weinberg equilibrium of . A total of 26,775 variants, including HLA alleles, amino acid polymorphisms, SNPs, and short indels remained for the following analysis, with most of them (63.6%) being common variants ().

3.3. Heritability of PFH and SP

To determine whether the genetic background of PFH is stronger than that of SP, GCTA (v1.93.1) [25] was used to calculate the heritability of PFH and SP under the assumption of a disease prevalence of 0.47%. We observed that variants in the MHC region collectively explained a higher heritability of PFH (, standard ) than SP (, standard ). This confirmed that PFH has a stronger genetic predisposition than SP in terms of the MHC region.

3.4. Associations of HLA with PFH

We performed a stepwise association analysis of variants with PFH susceptibility (1,538 PFH vs. 9,906 control). The most significantly associated variant mapped to the classical HLA-C06:02 allele (, ; Figure 1, Table 2), followed by HLA-C07:04 (, ). After conditioned on HLA-C06:02 and HLA-C07:04, there was no variant in HLA-C that reached the predefined threshold (). To look for additional independent variants, we conditioned on the HLA-C region and observed a novel significant risk variant for NOTCH4:G511S (, ). When we conditioned on HLA-C and NOTCH4, we observed an independent association at HLA-B amino acid position 67 (B-Y67C, , ), which was previously reported to be associated with psoriasis in both European and Han Chinese populations [18, 26]. After conditioning on HLA-C, NOTCH4, and HLA-B, we discovered an independent association for HLA-DPB1 amino acid position 35 (DPB1-F35L, , ). DPB1-F35L was in strong linkage disequilibrium with a previously reported psoriasis-associated allele, HLA-DPB105:01 () [18]. Further stepwise analysis identified another significant allele for TAP2:01:03 (, ). When we conditioned on all the above genes, no additional variant remained significantly associated with PFH susceptibility.

3.5. Associations of HLA with SP

As expected, in the fine-mapping analysis of SP (5,262 SP cases vs. 9,906 controls), we observed HLA-C06:02 (, , Figure 2, Table 3) as the most significantly SP-associated allele, followed by HLA-C07:04 (, ). After the condition on HLA-C06:02 and HLA-C07:04, we found an independent association for HLA-C01:02 (, ), which was recently reported to be associated with psoriasis in southern China [27]. We then investigated additional risk variants independent of HLA-C. When we included HLA-C as covariant, we observed significant independent associations at HLA-B amino acids 67 (B-Y67C, , ), 116 (B-Y116S, , ), and 70 (B-Q70K, , ). When we conditioned on HLA-C and HLA-B, we discovered a significant allele for HLA-DPB105:01 (, ). Then, we conditioned on HLA-C, HLA-B, and HLA-DPB1. We observed that HLA-A amino acid 95 (A-I95V, , ) exceeded the significance threshold. Further stepwise conditional analysis found that BTNL2:R281K (, ) was an independent risk variant. After conditioning on all the above genes, we observed no other significant associations.

3.6. Case-Case Association Analysis

Fine-mapping analysis of PFH and SP has identified a total of 11 independent susceptibility variants. To further determine if these 11 variants contribute to the family aggregation of psoriasis, we conducted a case-case analysis (PFH vs. SP), using a value of (Bonferroni correction for 11 variants, 0.05/11) as the cutoff for statistical significance. We found that HLA-C06:02 (, ) and NOTCH4:G511S (, ) confer risk for PFH, while BTNL2:R281K (, ) confers risk for SP (Table 4).

3.7. Annotation of NOTCH4:G511S and BTNL2:R281K

For variants NOTCH4:G511S and BTNL2:R281K, we used four annotation programs (see Materials and Methods) to measure the deleteriousness of variants. SIFT predicted that NOTCH4:G511S () and BTNL2:R281K () are deleterious. PolyPhen-2 predicted that BTNL2:R281K is probably damaging (), and NOTCH4:G511S is possible damaging (), respectively. Both CADD and DANN predicted that BTNL2:R281K (scaled and ) and NOTCH4:G511S (scaled and ) are the top 1% most deleterious variants.

3.8. No Epistasis between HLA-C06:02 and NOTCH4:G511S

Our stratified analysis suggested that NOTCH4:G511S and HLA-C06:02 contribute to the familial clustering of psoriasis. To further investigate whether there is an interaction between NOTCH4:G511S and HLA-C06:02, we performed an epistasis analysis by including interaction terms in the logistic regression model. In the case-case analysis, we included an interaction term between NOTCH4:G511S and HLA-C06:02 in our logistic regression model, and we observed no evidence of interaction between NOTCH4:G511S and HLA-C06:02 (). In the PFH versus control analysis, 15 interaction terms between all 6 significant variants were included in the logistic regression model; a value of 0.003 (0.05/15) was considered as statistical significance. We found strong interactions between HLA-C06:02 and TAP201:03 (, ) and HLA-C06:02 and DPB1-F35L (, ). But we still observed no evidence of interaction between HLA-C06:02 and NOTCH4:G511S after Bonferroni correction (, ). Our results suggested that HLA-C06:02 and NOTCH4:G511S may independently contribute to the familial clustering of psoriasis.

4. Discussion

In this study, we performed a large-scale stratified analysis to detect HLA and non-HLA variants in the MHC region that drive familial aggregation of psoriasis. In fine-mapping analysis of PFH and SP, we identified a total of 11 independent variants, of which eight variants (C06:02, C07:04, C01:02, HLA-B amino acid positions 67, 116, HLA-A amino acid position 95, DPB105:01, and BTNL2:R281K) were reported to be associated with psoriasis by previous studies [18, 2628]. Notably, HLA-C06:02, HLA-C07:04, HLA-B amino acid position 67, and HLA-DPB105:01 (strong linkage disequilibrium with DPB1-F35L) were shared susceptibility variants between PFH and SP, while other variants (NOTCH4:G511S and TAP201:03 for PFH, C01:02, B-Y116S, B-Q70K, A-I95V, and BTNL2:R281K for SP) showed specific associations with one of two psoriasis subtypes. Here, we have successfully decomposed the genetic architecture of PFH and SP into shared components and specific components. In order to further determine whether these shared or specific variants contribute to the family aggregation of psoriasis, we conducted a case-case analysis (PFH vs. SP) and eventually identified two variants (HLA-C06:02 and NOTCH4:G511S) in the MHC region which confer risk for family aggregation of psoriasis. Moreover, we also identified a variant (BTNL2:R281K) that specifically confers risk for SP.

For those shared variants between PFH and SP, only HLA-C06:02 was validated to be associated with family aggregation of psoriasis in the case-case analysis (, ). HLA-C06:02 is a well-known psoriasis-risk allele, and it was reported to be strongly associated with early age of onset of psoriasis [29]. This might explain why PFH have an earlier onset of psoriasis. Therefore, HLA-C06:02 conforms to our first hypothetical model that the shared variant contributes more risk to PFH than to SP. For those specific variants in PFH or SP, we found that NOTCH4:G511S (, ) contributes to the familial clustering of psoriasis. NOTCH4:G511S was proven to be an independent risk variant for PFH because the effect of NOTCH4:G511S with respect to PFH vs. control risk remained significant even after we conditioned on HLA-C, HLA-B, HLA-DPB1, HLA-A, TAP2, and BTNL2 (, Table S1A). In contrast, the impact of NOTCH4:G511S on SP vs. control risk disappeared after we conditioned on HLA-C and HLA-B (). Therefore, NOTCH4:G511S conforms to the second hypothetical model that the specific variant only confers risk for PFH and not for SP. NOTCH4 encodes members of the notch family of proteins that positively regulate notch signaling. The notch signaling pathway is a highly conserved molecular network that has a key role in angiogenesis [30], cellular proliferation and differentiation [31], and apoptosis [32]. Other studies have reported that the notch pathway regulates several critical aspects of epidermal renewal and development, including proliferation, differentiation, and cell-fate of keratinocytes [3336]. For example, a previous study showed that notch signaling induces differentiation of spinous cells into granulosa cells and, simultaneously, prevents premature differentiation of spinous cells, orchestrating the balance between differentiation and immature programs in suprabasal cells during epidermal development [36]. Moreover, the amino acid substitution glycine > serine at position 511 in NOTCH4 was predicted to be deleterious by all four annotation tools, suggesting this variant (NOTCH4:G511S) may cause aberrations in notch signaling and further lead to aggressive angiogenesis and the disruption of keratinocyte differentiation in psoriasis lesions. In addition, previous studies reported that NOTCH4 was associated with multiple autoimmune diseases, including Crohn’s disease, systemic scleroderma, rheumatoid arthritis, and ulcerative colitis [3740]. Of note, enthesitis is a key early manifestation of psoriasis arthritis. NOTCH4 was also reported to be associated with psoriasis arthritis [41], and notch signaling pathways were found to mediate synovial angiogenesis in psoriatic arthritis [42]. These findings suggested that notch signaling pathways were involved in the pathogenesis of PFH and NOTCH4:G511S may partially explain why PFH has more frequent enthesitis and other autoimmune diseases. We also conducted an epistasis analysis, finding no evidence of interaction between NOTCH4:G511S and HLA-C06:02. Therefore, NOTCH4:G511S and HLA-C06:02 may independently contribute to the family clustering of psoriasis based on our two hypothetical models.

Besides, we identified a variant (BTNL2:R281K) that specifically confers risk for SP. The amino acid substitution arginine > lysine at position 281 in BTNL2 was predicted to be deleterious by all four annotation tools. Notably, this variant was not associated with either PFH () or SP () in the univariate analysis (Table S1A). However, it exhibited a significant association with SP after conditioning on HLA-C06:02 (). Subsequently, we classified SP and control into two subgroups based on the presence of HLA-C06:02, the C06:02-positive group and the C06:02-negative group. We found that although the raw frequency of BTNL2:R281K in SP (0.059) was lower than in controls (0.063), this was reversed in both subgroups (0.086 in SP and 0.070 in controls for the C06:02-positive group, and 0.054 in SP and 0.036 in controls for C06:02-negative group, Table S1B). Moreover, univariate analysis in these two subgroups showed a consistent effect ( for the C06:02-negative group, for the C06:02-positive group). This phenomenon can be explained by Simpson’s paradox [43], where two subgroups have the same association, but the overall population shows no association. In contrast, BTNL2:R281K did not reach the significance threshold in PFH vs. control after conditioning on other susceptibility variants. The case-case analysis also showed that the frequency of BTNL2:R281K was significantly higher in SP than in PFH (). Therefore, BTNL2:R281K is an SP-specific risk variant that does not conform to our hypothetical models, which suggests that the genetic architecture of SP is not just a subset of PFH despite the fact that PFH has a stronger genetic predisposition than SP. We speculate that BTNL2:R281K might be associated with certain phenotypes of SP. For example, a recent study found that patients with SP have a higher Body Mass Index (BMI) than those with PFH [44]. BTNL2 was reported to be associated with the waist–hip ratio [45] and might contribute to higher BMI in SP.

However, there are limitations to our study. We were unable to perform further genotype-phenotype analysis because the included samples lacked clinical information, such as BMI, enthesitis, and other phenotypes of psoriasis. Besides, the functional mechanisms by which these variants contribute to familial aggregation of psoriasis or disease phenotypes remain unclear. Further studies will be required to address these problems.

In conclusion, we have identified two variants that independently contribute to family clustering of psoriasis, and one variant specifically confers risk for SP. We found NOTCH4:G511S and HLA-C06:02 could partially explain why patients with a family history have an earlier onset of psoriasis and more frequent enthesitis and other autoimmune diseases. Our findings support the hypothetical model that PFH has a stronger genetic background than SP, but also indicate that the genetic architecture of SP is not just a subset of PFH. Three variants identified in our study enable to differentiate PFH from SP and might serve as markers to predict the risk of developing psoriasis for individuals at a young age who have a family history.

Data Availability

The raw sequencing data from samples have been deposited in the Sequence Read Archive (SRA) under the accession SRA205317 (http://www.ncbi.nlm.nih.gov/sra/).

Conflicts of Interest

The author(s) declare(s) that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

We are grateful to all participants, their families, and healthy controls who donated blood samples in this study. We also thank LetPub (http://www.letpub.com) for its linguistic assistance during the preparation of this manuscript. This research was funded by the National Natural Science Foundation of China (81872527, 81130031).

Supplementary Materials

Table S1A: Detailed stepwise association results of NOTCH4:G534S and BTNL2:R281K with psoriasis risk. Table S1B: Phenomenon of Simpson’s paradox in the association of BTNL2:R281K with SP. (Supplementary Materials)