Introduction
Psoriasis (Ps) is a common inflammatory dermatologic condition that affects from 1% to 12% of population, depending on ethnicity and geographical location [1]. Aetiology of the disease is complex and it is generally considered a skin-specific T-cell mediated autoimmunological disease with risk factors involving genetic factors and environmental triggers [2]. The best characterized genetic risk factor for psoriasis is HLA-C*06 allele, encoding a variant of a major histocompatibility complex I (MHC I) antigen. However, HLA-C*06 does explain only a part of hereditability of psoriasis. Therefore, further investigations were undertaken in the last years to identify other genetic variants associated with the disease. Many single nucleotide polymorphisms (SNPs) were revealed by genome-wide association studies (GWAS), mostly localized in gene regions related to adaptive immunity, but including those connected with innate immune system pathways [3]. GWAS together with target genotyping performed in different ethnic populations resulted in identification of approximately 50 genetic markers associated with psoriasis [4]. Since reproducibility is very important in genetic investigations, in the current study we tried to verify an association between nine selected SNPs, previously linked to altered disease risk in GWAS. The investigated variants are located in genes encoding for proteins mainly involved in innate immunity. IFNLR1 (interferon lambda receptor 1, also known as IL28RA) may have an essential role in anti-viral activity and immune regulation, and it has been recently reported that it inhibits human epidermal keratinocyte proliferation in vitro [5]. STAT2, encoding one of ubiquitous transcription factors, was associated with psoriasis vulgaris in one of the GWAS, and subsequently it was demonstrated that it can be involved in the pathogenesis of psoriasis [6, 7]. The NFKBIA gene product is an inhibitor of the nuclear factor κB (NF-κB) pathway that may play a role in psoriasis pathogenesis, as reported increased expression of NF-κB in psoriatic skin may drive the proliferation of immune cells [8]. Similarly, TNFAIP3 (tumor necrosis factor, alpha-induced protein 3), TNIP1 (TNFAIP3-interacting protein 1), TRAF3IP2 (TRAF3 interacting protein 2) and REL (REL proto-oncogene, NF-kB subunit) are also involved in regulation of NF-κB associated inflammatory signalling pathways [3]. NOS2 encodes for inducible nitric oxide synthase playing an important role in metabolic and inflammatory processes, which is significantly upregulated in psoriatic lesional skin [9]. Finally, the IL12B gene product is a p40 subunit, shared by IL-23 and IL-12 cytokines. IL-23 induces differentiation of naïve T cells into Th17 and Th22 cells that in turn mediate the epidermal hyperplasia, keratinocyte immune activation and tissue inflammation [10].
Despite the strong evidence for an association between genetic factors with psoriasis, not many studies were focused on their impact on treatment efficacy up-to-date, and those published mostly refer to biopharmaceutics and psoriatic arthritis. Although the availability of monoclonal antibodies for psoriasis treatment is increasing, topical therapy combined with narrow-band ultraviolet B (NB-UVB) phototherapy is still one of the most effective treatments, especially in the case of mild or moderate disease severity, and an impact on genetic factors on the therapy outcome is unknown.
Aim
In the current paper we have investigated nine genetic variants in relation to their association with the disease, and subsequently in relation to response to topical and NB-UVB therapy in a cohort of Polish Caucasian patients with mild to moderate plaque psoriasis.
Material and methods
Study subjects
The current study involved 507 Polish psoriatic patients (361 type I and 146 type II psoriasis) and 396 healthy controls from the same geographical region, who were enrolled in the study after signing informed consent. The protocol of the study was approved by the local ethics committee. The diagnosis of psoriasis vulgaris was based on clinical examination and histopathological findings in selected cases. The main enrolment criteria included more than 1 year duration of psoriasis and presence of the active form of the disease. The patients were initially recruited for the purpose of our previous study [11]. More detailed characteristics of the subjects are presented in Table 1.
Table 1
Evaluation of genetic factors in relation to treatment outcome was performed in a subset of 306 patients with complete clinical data (49.7% females, mean disease onset at 30.6 ±13.7 years). The patients were classified as mild (n = 154) or moderate (n = 152) psoriasis based on the Psoriasis Area and Severity Index at the start of therapy (PASI ≤ 12 or 12–18, respectively). All patients were subjected to topical therapy, which consisted of vitamin D analogues (calcipotriol), steroids (betamethasone dipropionate glycol 0.05% or clobetasol 17-propionate 0.05%), salicylic acid 5%, urea 5–10%, and dithranol, administered according to the recommendations [12]. Additionally, 155 from 306 (50.7%) patients received NB-UVB phototherapy (Cosmedico Medizintechnik GmbH, Germany, or Series 3 Daavlin Philips AG) 3 times a week throughout the study period, based on skin phototype according to the manufacturer’s recommendations. Patients medicated with systemic immunosuppressive regimens and PUVA (psoralen + ultraviolet A) treatment were excluded from the study. Patients were evaluated at baseline, and then after 8 weeks of therapy. The percent reduction of PASI score was evaluated as a variable in the analysis of the treatment outcome. Mean PASI value (± SD) at the start of therapy was 11.52 ±4.52, range: 3.2–18.0, and after 8 weeks of therapy: 3.78 ±4.05, range: 0.0–17.2 (mean PASI reduction: 68.64 ±29.56%, 7.74 ±4.53 points, range: 2.52–100.0%, 0.2–19.6 points).
Genotyping
Genomic DNA was extracted from buccal swab samples using Genomic Micro AX SWAB Gravity kit (A&A Biotechnology, Gdynia, Poland) and from peripheral blood samples using a GeneMATRIX Quick Blood DNA Purification Kit (EURx, Poland) following the manufacturer’s protocol. DNA was subsequently standardized to equal concentrations of 20 ng/µl, based on spectrophotometric absorbance measurement (260/280 nm). Genotyping for the following single nucleotide polymorphisms (SNPs): IFNLR1 rs7552167, IL12B rs12188300, STAT2 rs2066819, NFKBIA rs8016947, NOS2 rs28998802, REL rs702873, TRAF3IP2 rs33980500, TNIP1 rs2233278 and TNFAIP3 rs582757 was performed using a pre-validated allelic discrimination TaqMan real-time PCR assay (Life Technologies, USA), and TaqMan GTXpress Master Mix (Life Technologies, USA). Fluorescence data were captured using the ViiA7 Real-Time PCR System (Applied Biosystems, USA) after 40 reaction cycles. Specific genotypes were assigned to individual samples after analysis with TaqMan Genotyper software (Thermo Fisher Scientific, USA). The information about SNP location and TaqMan assay used is given in Table 2.
Table 2
Statistical analysis
The χ2 and Fisher’s exact test was used to compare qualitative variables between genotype groups, alleles as well as clinical data in genotype groups, and accordance of genotypes with Hardy-Weinberg equilibrium was calculated. Analysis of PASI reduction was performed by means of the Kruskal-Wallis test. P < 0.05 was considered statistically significant. Statistical analysis was performed using Statistica ver. 13.1 software (TIBCO Software Inc., Palo Alto, CA, USA).
Results and Discussion
Distribution of genotypes was in accordance with Hardy-Weinberg equilibrium for all studied SNPs, both in psoriatic patients and healthy controls (p > 0.05). Among the analysed variants, IL12B rs12188300, TRAF3IP2 rs33980500 and TNFAIP3 rs582757 minor alleles have been significantly more frequent among psoriatic patients (p < 0.05) (Table 3). When analysing genotypes we observed an association of the studied IL12B, TRAF3IP2, TNIP1 polymorphisms in a dominant model (significantly higher frequency of carriers of at least one risk allele among psoriatic patients), while for NOS2 the odds for disease were significantly increased only in homozygotes with two minor alleles (recessive model). A marginally significant difference was also observed in case of TRAF3IP2 in the recessive model. Polymorphic variants of IFNLR1, STAT2, NFKBIA and REL genes did not show any association with psoriasis in the analysed population. More detailed results of genotype analysis is given in Table 4. The variants investigated in the current study were initially identified as psoriasis risk modifiers in GWAS studies [13–19] and the associations were further investigated in many confirmation studies, where both positive and negative results were reported. Kisiel et al. confirmed association of IFNLR1 rs7552167 (p = 0.017), NFKBIA rs8016947 (p = 7.5 × 10–3) and REL rs702873 (p = 2.5 × 10–4) in the Polish population [20]. However, we have not observed any difference in the frequency of genotypes or alleles between psoriatic patients and controls (p = 0.276, p = 0.474 and p = 0.307, respectively), it might be due to a less numerous control group in our study or due to exclusion of patients with psoriatic arthritis from the patient group. In the German population IL12B rs12188300 and TRAF3IP2 rs33980500 SNPs were strongly associated with psoriasis vulgaris (p = 0.027, p = 1.5 × 10–9), as confirmed in the present study on the Polish population (p = 0.036, p = 0.021) [21]. Large scale meta-analysis of the psoriasis-associated TNFAIP3 SNPs pointed rs582757 as top Ps variant (p = 6.07 × 10-12, OR = 1.23) [22]. Our results confirm the TNFAIP3 rs582757 association with psoriasis vulgaris (p = 0.046). STAT2 rs2066819, one of the novel Ps susceptibility SNP loci did not reach statistical significance in our study, possibly due to a limited sample size [4, 19]. Similar to our study, replicated GWAS analysis in the Japanese population showed TNIP1 rs2233278 association with psoriasis (p = 7.0 × 10–6). TNFAIP3-interacting protein 1, TNIP1 gene coded protein, is an inflammation regulator in autoimmune diseases, involved in internal signalling pathways, like regulation of an NF-κB pathway and microRNA [23, 24].
Table 3
Gene | SNP id | MAF psoriasis | MAF controls | P-value* | Genotypes – dominant model** | Genotypes – recessive model*** | ||
---|---|---|---|---|---|---|---|---|
OR (95% CI) | P-value | OR (95% CI) | P-value | |||||
IFNLR1 | rs7552167:A>G | 0.133 | 0.152 | 0.276 | 0.86 (0.64–1.17) | 0.357 | 1.41 (0.62–3.23) | 0.524 |
IL12B | rs12188300:A>T | 0.118 | 0.087 | 0.036 | 1.53 (1.09–2.16) | 0.014 | 1.29 (0.41–4.01) | 0.773 |
STAT2 | rs2066819:C>T | 0.096 | 0.090 | 0.684 | 1.07 (0.75–1.51) | 0.790 | 0.80 (0.26–2.46) | 0.784 |
NFKBIA | rs8016947:T>G | 0.437 | 0.455 | 0.474 | 1.07 (0.81–1.42) | 0.664 | 1.37 (0.99–1.89) | 0.066 |
NOS2 | rs28998802:G>A | 0.150 | 0.122 | 0.099 | 1.19 (0.88–1.62) | 0.281 | 0.31 (0.10–0.94) | 0.038 |
REL | rs702873:C>T | 0.376 | 0.400 | 0.307 | 0.81 (0.61–1.06) | 0.126 | 0.99 (0.68–1.45) | 1.000 |
TRAF3IP2 | rs33980500:C>T | 0.130 | 0.095 | 0.021 | 1.42 (1.02–1.97) | 0.040 | 0.42 (0.14–1.32) | 0.137 |
TNIP1 | rs2233278:G>C | 0.083 | 0.059 | 0.067 | 1.49 (1.01–2.19) | 0.045 | 1.28 (0.08–20.55) | 1.000 |
TNFAIP3 | rs582757:C>T | 0.304 | 0.260 | 0.046 | 1.22 (0.94–1.59) | 0.141 | 0.60 (0.37–0.99) | 0.053 |
Table 4
[i] P-values calculated by means of Fisher exact test, p1 – heterozygotes vs. major homozygotes, p2 – minor homozygotes vs. major homozygotes, p3 – heterozygotes + minor homozygotes vs. major homozygotes (dominant model), p4 – minor homozygotes vs. major homozygotes + heterozygotes (recessive model).
The vast majority of psoriasis risk SNP variants are situated near the genes encoding molecules involved in the innate immune system [4]. As many of them do not affect directly the protein sequence or function (intronic, intragenic and synonymous SNPs), it is probable that interethnic differences can affect the outcome of association studies. Hence, only extensive meta-analyses taking into account the results of many heterogeneous studies performed in different ethnic groups could finally verify the impact of individual SNPs on the psoriasis risk. For that reason the results of that study not only provide data for psoriasis genetics in the Polish population, but also could be of some value for future meta-analyses.
Subsequently, the response to therapy was analysed in relation to the patients’ genotype. However, no significant impact on efficacy of topical therapy combined with NB-UVB was noted in case of any investigated SNPs (Table 5). Additional analysis was performed for possible association of genetic factors with the age of disease onset, but also in this case the results were negative (data not shown). Contrary to the current negative results in relation to topical therapy, some evidence was provided that the genetic polymorphism may be responsible for different patient response to treatment with biologic drugs that are targeted to inflammatory cytokines involved in developing psoriasis [25]. Among the genes investigated in the current study, IL12B, IFNLR1, TNFAIP3 polymorphisms were reported to affect treatment with anti-TNF (adalimumab, etanercept) or anti-IL-12/IL-23 (ustekinumab) agents. However, biologic therapy is much more specific, and the polymorphism within genes involved in pathogenesis can directly affect potency of therapeutic antibodies and their interaction with targets. On the contrary, topical treatment is generally less specific, and the relation between drugs and genes involved in psoriasis pathogenesis is not so clear. Possibly, other genetic variants may play a role, like it was described earlier in case of vitamin D receptor (VDR) and calcipotriol therapy [26].
Table 5
Conclusions
The results of the current study confirmed an association of IL12B, TNFAIP3, TRAF3IP2, TNIP1 and NOS2 variants with psoriasis vulgaris in the Polish population. However, the analysed SNPs do not seem to have any impact on efficacy of topical/NB-UVB therapy, and are not potentially useful in personalized therapy.