Yunqing Ren1,2
Ling Wang3Huatuo Dai1,2Guiying Qiu1,2Jipeng Liu1Dianhe Yu1,2
Jianjun Liu3,4Cheng-Zhi Lyu5
Lunfei Liu2,6*
Min Zheng2*- 1Department of Dermatology, The Children’s Hospital, Zhejiang University School of Medicine, National Clinical Research Center for Child Health, Hangzhou, Zhejiang, China
- 2Department of Dermatology, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- 3Genome Institute of Singapore, Agency for Science, Technology and Research, Singapore, Singapore
- 4Department of Medicine, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore
- 5Department of Dermatology, Dalian Dermatosis Hospital, Dalian, China
- 6Department of Dermatology, The Fourth Affiliated Hospital, Zhejiang University School of Medicine, Yiwu, Zhejiang, China
Background: TNF-α inhibitors are effective biological agents for treating psoriasis, but the treatment responses differ across patients. This study aimed to identify genetic biomarkers of anti-TNF-α response in Chinese psoriasis patients using a genome-wide association approach.
Methods: We recruited two independent cohorts of Chinese psoriasis patients administered etanercept biosimilar (with or without methotrexate). We identified 61 and 87 good responders (PASI improvement ≥75%), 19 and 10 poor responders (PASI improvement <50%) after 24 weeks treatment in the two cohorts, respectively. Then we performed genome-wide association studies (GWAS) on anti-TNF-α response in each cohort independently, followed by a fixed-effects inverse-variance meta-analysis in the 148 good and 29 poor responders.
Results: We tested genetic associations with >3 million genetic variants in either cohort. Meta-analysis identified significant associations within seven loci at p < 10−5, which also showed consistent association evidence in the two cohorts. These seven loci include rs2431355 (OR = 6.65, p = 4.46 × 10−7, IQGAP2-F2RL2 on 5q13.3), rs11801616 (OR = 0.11, p = 1.75 × 10−6, SDC3 on 1p35.2), rs3754679 (OR = 0.17, p = 7.71 × 10−6, CNOT11 on 2q11.2), rs13166823 (OR = 0.09, p = 3.71 × 10−6, IRF1-AS1 on 5q31.1), rs10220768 (OR = 5.49, p = 1.48 × 10−6, NPAP1 on 15q11.2), rs4796752 (OR = 5.56, p = 1.49 × 10−6, KRT31 on 17q21.2), and rs13045590 (OR = 0.08, p = 9.67 × 10−7, CTSZ on 20q13.3). Of the seven SNPs, six SNPs showed significant eQTL effect (p < 1 × 10−6) for several genes in multiple tissues.
Conclusion: These results suggest novel biological mechanisms and potential biomarkers for the response to anti-TNF therapies. These findings warrant further validation.
Introduction
Psoriasis is an immune-mediated, chronic, inflammatory skin disorder influenced by genetic and environmental factors (Boehncke and Schon 2015; Greb et al., 2016). It is characterized by infiltration of leucocytes into the dermis and epidermis, excessive proliferation and changes in keratinocytes differentiation, leading to red, scaly skin, abnormal desquamation and defects in epidermal barrier. The incurable disease and the need for long-term medication bring a heavy life and psychological burden on patients. In the past decade, several biological agents have been exploited for the treatment of moderate to severe psoriasis greatly improved the situation of the patients. Among them, TNF-α inhibitors such as etanercept, adalimumab and infliximab are the most commonly used biologic drugs. However, there is considerable heterogeneity in the response to this drugs whose efficacy is range from 50% to 80% (Boehncke and Schon 2015). Accompanied with the expensive treatment expense and the potential side effects associated with the treatment, therefore, screening patients who are suitable for TNF-α inhibitors, in other words, patients with better efficacy of treatment may have a certain clinical significance.
Previous studies have identified several clinical factors such as dermoscopic hemorrhagic dots and BMI at baseline that influence the response to anti-TNF-α therapy in psoriasis (Lallas et al., 2016; Zweegers et al., 2017). Consequently, effort has been also put into the identification of genetic variation predicting anti-TNF treatment outcome. Most of these pharmacogenetics studies are candidate gene based, focusing on polymorphisms in genes known to be involved in psoriasis susceptibility or immune system such as TNF-α signaling pathways, including TNF, TNF receptor superfamily member 1A (TNFRSF1A), TNFRSF1B, TNF-α-induced protein 3 (TNFAIP3), IL12B/IL23R, CARD14, LCE, HLA-Cw6 and so on (Vasilopoulos et al., 2012; Gallo et al., 2013; Batalla et al., 2015; Coto-Segura et al., 2016; Liu et al., 2016; Ovejero-Benito et al., 2017; Dand et al., 2019; van den Reek et al., 2017; Song et al., 2015). However, the majority of these studies have generated inconsistent results with small sample size and without independent validation.
GWAS has been recently used to identify new genetic biomarkers associated with the differential response to anti-TNF agents (Liu et al., 2008; Umicevic Mirkov et al., 2013; Honne et al., 2016; Julia et al., 2016), and only two studies were performed in psoriasis. The first study included 65 psoriasis patients from Japan (Nishikawa et al., 2016), and the other one included 243 psoriasis patients from Spain (Ovejero-Benito et al., 2020). However, no genome wide significant SNPs have been identified. The majority of the SNPs reported by GWAS are not known psoriasis susceptibility genes. In addition, none of the SNPs reported by candidate gene studies were confirmed by the GWAS analyses.
Here, we report a meta-GWAS analysis of two independent cohorts with a total of 212 patients treated with etanercept biosimilar. We have identified seven associated genetic biomarkers of anti-TNF drug response in Chinese Han population, which provided novel biological insight for clinical response of TNF-α inhibitor treatment.
Materials and methods
Patient cohorts and clinical evaluation
Two cohorts of psoriasis patients were included in this study. Cohort 1 included 103 patients and cohort 2 comprised 109 patients. All the patients 1) were Chinese descent; 2) had a psoriasis area and severity index (PASI) score >10 and an affected body surface area (BSA) > 10% at baseline; 3) were treated with etanercept biosimilar 50 mg per week. While the patients of the cohort 1 were treated only with etanercept biosimilar, the cohort 2 was from a clinic trial of combination of etanercept with methotrexate, where the patients were treated with etanercept biosimilar 50 mg per week and a stable dose of methotrexate 15 mg per week (definite dose depends on the patient’s tolerance). Methotrexate is one of the most common and conventional disease-modifying anti-rheumatic drugs used for psoriasis and psoriatic arthritis. Detailed clinical parameters collected includes present age, sex, disease duration, smoking and alcohol drinking history, family history, psoriatic arthritis, previous systemic therapies, and PASI score at first screening and during subsequent visits. The study was approved by the local institutional ethics review board, the second affiliated hospital of Zhejiang University school of medicine, and written informed consent was obtained from all patients.
Efficacy measurements
We used PASI score to evaluate the clinical severity of disease for all the patients. We calculated response as the percentage improvement in the PASI score at 12 and 24 weeks compared with that at baseline. In this study, the patients had a better therapeutic effect at week 24 than that at week 12. Patients were categorized as poor responders, if their PASI improvement rates were <50%, or good responders, if their PASI improvement rates were ≥75%. Similar to previous studies, moderate responders (whose PASI improvement rates were between 50% and 75%) were excluded for association analysis (Masouri et al., 2016; Mendrinou et al., 2016).
Genotyping and quality control
Genomic DNA of each sample was extracted from peripheral blood using standard protocols. Genotyping analysis of the two cohorts was conducted using the Illumina Infinium Asian Screening Array-24 v1.0 BeadChip. We performed quality control (QC) on the two cohorts respectively. We excluded SNPs from X, Y, and mitochondrial chromosomes from GWAS analysis. All samples had good quality with high call rate (>95%) and no extreme heterozygosity (<3SD). Cryptic relatedness between individuals was identified through identity by descent analysis using PLINK (see URL) and two first degree relative pairs were identified in cohort 2. One sample of each pair with lower SNP call rate were removed. Principle component analysis was performed using PLINK, one outlier was identified and removed in cohort 1. Totally, 102 samples for cohort 1 and 107 samples for cohort 2 passed QC and remained in further analysis.
Genotype imputation and quality control
To increase overall coverage of the genome, we performed imputation using IMPUTE v2 based on the 1,000 Genomes Project reference panel from all populations (dated October 2014). We imputed on basis of samples with high call rate and the genotyped SNPs that passed quality control (call rate≥95%, HWE p ≥ 1 × 10−6, MAF≥0.01). SNP quality control after imputation were as follows: imputation quality score ≥0.8, HWE p ≥ 1 × 10−6, MAF≥10% in both good responders and poor responders.
Statistical analyses
The difference of clinical characteristics of the patients between the two cohorts was calculated by Chisq test or t test using R 4.0.2. The genome-wide case-control association analysis was conducted using SNPTEST v2 (see URLs) and imputed dosage data in each cohort. The association test was performed by assuming an additive model for allelic risk, with the age, gender, BMI at baseline and the first three principal components included as covariates to adjust for study effect and population stratification. To combine the association statistics from the two independent GWASs, we performed a fixed-effects inverse-variance meta-analysis of 3478293 overlapped SNPs using META (see URLs). Cochran’s Q test and I2 were used to assess between-study heterogeneity.
Power calculation was performed for the top SNPs detected in the current study using Quanto (see URL) under following assumption: disease prevalence of 0.01% and a log-additive risk model. For the calculation, we used effect size estimates observed in the combined population and the observed allele frequencies in all samples. The current sample size has a moderate power (ranging from 45% to 93%) to detect the effect at a p-value of 1 × 10−5.
We first determined the SNPs and indels that are in high LD (r2 ≥0.8) with seven lead SNPs using the LD information provided by HaploReg V3 on the 1,000 Genomes phase 1 Asian (ASN) population, which identified 86 SNPs and small indels. Then the mRNA expression level (eQTL) effect was evaluated through Genotype-Tissue Expression (GTEx) portal (see URL).
Results
Study subjects
The baseline characteristics of the two cohorts were summarized in Table 1. For the Cohort 1, the average age of the patients at baseline was 42.5 years, with 72.6% being male. The mean PASI score at baseline was 24.4 ± 12.3, which gradually decreased during the treatment and achieved the best curative effect at 24 weeks. For the Cohort 2, the baseline average age, proportion of male and PASI score was 41.8 ± 13.2, 78.5% and 22.3 ± 9.3, respectively. There was no significant difference in baseline age, gender ratio, BMI, personal history between the two cohorts as well as the PASI score at baseline. However, the overall response to treatment in term of PASI reduction was better for the Cohort 2 than the Cohort one due to the combined using of methotrexate.
TABLE 1. Basic clinical characteristics of the patients in this study.
Genome-wide association analysis
To minimize the adverse impact of the misclassification of treatment response, we just analyzed the patients of good responders (PASI improvement ≥75%) and poor responders (PASI improvement <50%) after 24 weeks treatment. There were 61 good responders and 19 poor responders for the Cohort 1, and 87 good responders and 10 poor responders for the Cohort 2. Because the number of poor responders was moderate for both cohorts (19 for the cohort 1 and 10 for the cohort 2), only the SNPs whose MAFs were above 10% in the samples of poor responders were included in the association analysis. We further performed PCA analysis of all the patients and found that good and poor responders were well matched genetically without evidence of population stratification (Supplementary Figures S1–S3).
We first performed the genome-wide association analysis in the two cohorts independently: 3815035 autosomal SNPs for the Cohort 1, and 3754575 autosomal SNPs for the Cohort 2. To assess and control the genetic heterogeneity of association effect between the two independent cohorts, we performed a Meta-GWAS analysis of the two cohorts where the association results of 3478293 overlapping SNPs from the two cohorts were combined by a meta-analysis under a fixed-effects model. The λgc of Meta-GWAS analysis was 1.065, indicating a minimal inflation of the genome-wide association results due to population stratification, which is consistent with the result of PCA analysis. Consistently, the quantile-quantile plot of the observed p values for association (Figure 1) showed that the distribution of the logarithms of the observed p values largely fits the null distribution.
FIGURE 1. Results of meta-GWAS of the two cohorts. (A) The −log10 genome-wide p values using a fixed-effects inverse-variance meta-analysis calculated on the basis of data from overlapped 3478293 SNPs in the 148 good responders and 29 poor responders. The blue horizontal line presents a liberal threshold of 10−5 for suggestive significance. The red horizontal line indicates genome-wide significance (p < 5 × 10−8). (B) Quantile-quantile plots of the observed p values versus the expected values form p value of association. λgc = 1.065.
The Meta-GWAS analysis revealed significant association at seven loci with a p-value less than 10−5 (Table 2). The top association was observed at rs2431355 (OR = 6.65, p = 4.46 × 10−7, IQGAP2-F2RL2 on 5q13.3), and the other significant associations were found at rs11801616 (OR = 0.11, p = 1.75 × 10−6, SDC3 on 1p35.2), rs3754679 (OR = 0.17, p = 7.71 × 10−6, CNOT11 on 2q11.2), rs13166823 (OR = 0.09, p = 3.71 × 10−6, IRF1-AS1 on 5q31.1), rs10220768 (OR = 5.49, p = 1.48 × 10−6, NPAP1 on 15q11.2), rs4796752 (OR = 5.56, p = 1.49 × 10−6, KRT31 on 17q21.2), and rs13045590 (OR = 0.08, p = 9.67 × 10−7, CTSZ on 20q13.3). For all the seven loci, the minor allele of the top SNP was associated with poor response. For each of the seven loci, the association was detected in both cohorts, and the estimated genetic effect of association (OR beta value) was consistent between the two cohorts without any evidence of genetic heterogeneity (Table 2). Furthermore, for each of the seven loci, the evidence of association was supported by both genotyped and imputed SNPs. The regional plots of the seven loci are provided in Supplementary Figures S4–S10.
TABLE 2. Meta-analysis results of the seven loci for anti-TNF response in psoriasis.
To further investigate the association of these seven SNPs at 12 weeks, a fixed-effects inverse-variance meta-analysis was also performed in the patients of good responders (PASI improvement ≥75%) and poor responders (PASI improvement <50%) after 12 weeks treatment. There were 46 good responders and 31 poor responders for the Cohort 1, and 65 good responders and 15 poor responders for the Cohort 2. These seven top SNPs association analysis results at 12 weeks are shown in Supplementary Table S1, which showed consistent trends of association as the analysis at 24 weeks.
In addition, we checked the associations for the 19 SNPs that were reported with p < 4 × 10−5 in the two published GWAS on anti-TNF-α response (Nishikawa et al., 2016; Ovejero-Benito et al., 2020); Two of these SNPs (rs77497886 and rs1487419 on CDH12) achieved nominal significance in our study (p = 0.031 and 0.032).
EQTL analysis and functional annotation
To better understand the biological function, we have investigated the regulatory functions of these seven top significant SNPs by analyzing the information from the HaploReg (v3) and the GTEx database. Of the seven SNPs, four SNPs were found in the database, rs2431355, rs13045590, rs3754679 and rs4796752, and showed significant eQTL effect (p < 1 × 10−6) for several genes in multiple tissues. The top SNP rs2431355 conferred significant eQTL effect in multiple tissues for F2RL2 gene (p = 3.70 × 10−19). rs13045590 was associated with significant eQTL effect (p = 9.50 × 10−23) on the expression level of CTSZ gene in muscle-skeletal and heart. rs4796752 was showed most significant eQTL effect (p = 1.80 × 10−11) with KRT31 gene in skin, and rs3754679 was showed significant eQTL with CNOT11 gene in breast tissue with p = 1.00 × 10−6. rs11801616 and rs10220768 were absent in the database, but 13 SNPs in high LD (r2 > 0.8) were also analyzed, which showed conferred eQTL effect for SDC3 (p = 6.40 × 10−9) and NPAP1 gene (p = 7.80 × 10−6) respectively. rs13166823 and SNPs in high LD had no eQTL results in GTEx portal. The eQTL analysis results were summarized in Supplementary Table S2.
Discussion
In this study, we performed the first genome-wide association study to evaluate pharmacogenetic effect on the response to etanercept biosimilar in Chinese population. By investigating two independent cohorts of patients of psoriasis treated with etanercept biosimilar (with or without methotrexate), we were able to evaluate the robustness of the association by comparing the consistence of statistical evidences from independent samples using genetic heterogeneity test. By interrogating 3478293 SNPs in a total of 148 good responders and 29 poor responders, the current study is the largest effort in Asian population and has discovered seven loci that were significantly associated with treatment outcome. Although the statistical evidences of association for these seven loci did not reach genome-wide significance (p < 5E-08), the robustness of the association evidence is reflected by the consistent evidence from the two independent cohorts as well as supported by multiple SNPs within each locus. Furthermore, the associations at these seven loci were also supported by the treatment outcome data of these two cohorts at 12 W. As previous studies, we have excluded the subjects with moderate response (PASI improvement between 50% and 75%), to minimize the misclassification of treatment response and thus enhance the robustness of association results.
The strongest association with treatment response was found in the IQGAP2- F2RL2 locus on 5q13.3. IQGAP2 gene encodes the IQ motif containing GTPase activating protein 2, which regulates a wide variety of cellular signaling pathways. IQGAP2 gene has been shown to function in thrombin-induced platelet cytoskeletal reorganization (Schmidt et al., 2003). In addition, the top SNP rs2431355 conferred significant eQTL effect in multiple tissues for F2RL2, which also plays an essential role in hemostasis and thrombosis. Variants in IQGAP2 have previously been associated with platelet count and mean platelet volume (MPV) (Eicher et al., 2016). Interestingly, accumulating evidence has shown that platelet dysfunction plays an important role in the pathogenesis of psoriasis (Fan et al., 2020). Besides as the principal mediator of hemostasis and thrombosis, platelets are also immune cells that involve in immune and inflammatory processes (Zamora, Canto, and Vidal 2021). Platelets aggregation has been detected in patients with psoriasis. Several studies have suggested that MPV values were significantly higher in patients with psoriasis and correlated with disease severity, also decreased after effective biological therapy (Capo et al., 2014; Kilic et al., 2017; Fan et al., 2020). Moreover, the levels of platelet-lymphocyte complexes and platelet to lymphocyte ratio were also indicated to be prognostic biomarkers of response to anti-TNF-α treatment (Najar Nobari et al., 2020; Sanz-Martinez et al., 2020). Taking together, our current finding and previous evidences suggest that genetic variation influencing platelets activity could influence the resistance or efficacy to anti-TNF agents. Further studies confirming the functional implication of IQGAP2/F2RL2 in psoriasis and the response to anti-TNF agents are warranted.
In addition, 12 SNPs on chromosome 20q13 were found to be associated with treatment outcome. The most significant SNP rs13045590 was in an intron of the cathepsin Z (CTSZ) gene. Using the GTex data, we also found that rs13045590 was associated with significant eQTL effect on the expression level of CTSZ gene in muscle-skeletal and heart. CTSZ is a member of the cysteine cathepsins, which regulates various cellular physiological functions, including adhesion, migration, invasion, and maturation of immune cells. Recently this SNP was found to be associated in Platelet distribution width (PDW) (Vuckovic et al., 2020). It is also interesting to note that several variants within this locus (CTSZ, TUBB1 and PRELID3B) had been previously reported to be associated with MPV and platelet count (Gieger et al., 2011; Qayyum et al., 2012). Therefore, the association identified in 20q13 has further implicated the link between platelets activity and anti-TNF-α treatment response in psoriasis. Validation studies will be required to explore this further.
Another suggestive locus is identified at 5q31.1, which is immediately adjacent to the locus that was previously reported to be associated with psoriasis and Crohn’s disease (Rioux et al., 2001; Chang et al., 2008). This locus contains a cluster of cytokine and immune-related genes including interferon regulatory factor 1 (IRF1), IL5 and IL13 as well as SLC22A4 and SLC22A5. Rs13166823 is located within IRF1 antisense RNA 1(IRF1-AS1). The protein encoded by IRF1 gene is a transcriptional regulator and tumor suppressor, playing a role in cell proliferation, apoptosis, DNA damage response as well as innate and acquired immune responses. Previous study revealed that IRF1 expression decreased in lesional psoriatic skin compared with nonlesional or normal skin. And the induction of IRF1 activity to IFN-γ was also greatly reduced in psoriatic keratinocytes vs. normal keratinocytes (Jackson et al., 1999). Recent studies indicated that IRF1 could mediate M1 polarization of macrophages which may contribute to psoriasis. However, Lin et al. found that TNF-α blockage inhibits M1 polarization through STAT1-and IRF1-independent pathways. Furthermore, rs2070729 in IRF1 has been previously associated with platelet count (Gieger et al., 2011). These findings suggest that IRF1 may be involved in the pathogenesis of psoriasis whose genetic variation may affect the response during psoriasis treatment.
Rs10220768 is located in Prader-Willi region non-protein coding RNA 4 (PWRN4). Using the published eQTL data, we have shown that rs10220768 and six linked SNP were associated with significant eQTL effect on the expression level of nuclear pore associated protein 1 (NPAP1) gene in thyroid. NPAP1 gene encodes nuclear pore associated protein one which also located chromosome 15 and defects in this gene may be associated with Prader-Willi syndrome. A recent study has identified that NPAP1 is one of the potential genes which associated with the response to immune checkpoint inhibitors therapy in patients with non-small cell lung cancer (Yu et al., 2019). Considering the high efficacy in psoriasis therapy of TNF-α inhibitors, it is likely that these SNP-associated genes could be involved in the response of immunotherapy in psoriasis.
The SNP rs4796752 on 17q21.2 is located within the region of the keratins genes clusters. Using the GTex data, eQTL analysis showed that rs4796752 was most significantly associated with KRT31 gene in skin. KRT31 is a type I acidic keratin and played an essential role in maintaining the health state of hair follicles. Application of K31 improved the damaged hair in smoothness, diameter and mechanical strength (Basit et al., 2018). Previous evidence showed K31 peptides could contribute to alopecia areata induction by supporting T cell activation when presented by DCs to syngeneic naive T cells in mice (Erb, Freyschmidt-Paul, and Zoller 2013). More and more evidence showed that connection between psoriasis and hair follicles on shared molecular regulatory mechanism and therapeutic targets (Suzuki et al., 2022). Patients with psoriasis have an OR of 2.5 for developing alopecia areata. It has also been reported that alopecia may be related to psoriasis itself or systemic therapies used to treat it especially anti-TNF-α agents. Of relevance, treatment with anti-TNF-α inhibitors can also cause alopecia. Therefore, a potentially interesting connection exists between KRT31 and psoriasis as well as the efficacy of TNF blockade.
Twelve SNPs in the locus on chromosome 1p35.2 had p value < 10–5, with the most significant SNP rs11801616 located in the 5 kb upstream of 5′ of RN7SK pseudogene 91 (RN7SKP91). Interestingly, seven SNPs in high LD with rs11801616 were associated with significant eQTL effect on the expression level of SDC3 gene in Artery. Syndecan-3, encoded by SDC3 gene, belongs to the family of heparan sulfate proteoglycans (HSPGs) and mainly involves in brain development and feeding behaviors. Recent studies suggested that SDC3 has important roles in inflammation and angiogenesis (Arokiasamy et al., 2019). Increased expression of SDC3 has been shown on endothelia in rheumatoid arthritis and psoriasis arthritis synovium comparing to normal synovium, which is thought to have roles in angiogenesis, chemokine presentation and leukocytes extravasation (Patterson et al., 2008). Furthermore, we recently identified several pathways and functional processes including HSPG and VEGF production which are correlated with anti-TNF-a treatment (Liu et al., 2019). Therefore, SDC3 has the possibility to be a crucial part in the pathogenesis of psoriasis and influence the therapeutic response in psoriasis.
Rs3754679 on 2q11.2 is located within Neuronal PAS domain protein 2 (NPAS2) gene. The eQTL analysis showed that rs3754679 was associated with significant eQTL effect on the expression level of CNOT11 gene in mammary breast. CNOT11 was one of subunits of the CCR4 (carbon catabolite repressor 4)-NOT (Negative on TATA) complex, which plays a crucial role in post-transcriptional mRNA regulation in eukaryotes. Previous evidence showed CCR4-NOT complex contributes to repression of Major Histocompatibility Complex class II transcription (Rodriguez-Gil et al., 2017). Although the role of CNOT11 subunit and how it is integrated into the complex are still unknown, it was found to be significantly upregulated in psoriatic skin lesions. These findings suggested a potential association between CCR4-NOT and psoriasis as well as the efficacy of TNF blockade.
A major limitation of the current study is that the seven loci did not reach genome-wide significance, the same as the two previous GWAS study carried out in 65 Japanese patients and 243 Caucasian patients. Although we reduced the adverse impact of misclassification on the results by excluding the partial responders in the association analysis, the reduced sample size would decrease the statistical power of this study and thus weaken the statistical significance of the results. In addition, the treatment regime was not identical between the two cohorts. The patients of the Cohort 1 were treated with etanercept biosimilar alone, while the patients of the Cohort 2 were treated with the combination therapy of etanercept biosimilar and methotrexate. Methotrexate, a folic acid antagonist, is a conventional systemic agent used widely for psoriasis. Adding methotrexate to etanercept has been proved to increase the efficacy compared with etanercept alone (Gottlieb et al., 2012). Although the meta-analysis could control the genetic heterogeneity of association effect, it is still difficult to avoid the effect of increased efficacy on results of Cohort 2 due to combined using of methotrexate. Therefore, our current findings require further validation in independent samples and/or at a functional level.
In the present study, we have identified seven genetic loci that show evidence of influencing anti-TNF treatment response based on a meta-GWAS in 177 Chinese psoriasis patients, which might serve as predictors for treatment and guide personalized therapeutic decisions. These genetic associations also provided biological insight for how TNF-α inhibitor treatment, suggesting the existence of different biological mechanisms through which TNF-α inhibitors exert their immunomodulatory role in psoriasis.
URLs:
PLINK: http://www.cog-genomics.org/plink/1.9/
IMPUTE: https://mathgen.stats.ox.ac.uk/impute/impute_v2.html
SNPTEST: https://mathgen.stats.ox.ac.uk/genetics_software/snptest/snptest.html#download
META: http://mathgen.stats.ox.ac.uk/genetics_software/meta/meta.html#
Quanto: https://bio.tools/QUANTO
GTEx: https://gtexportal.org/home/
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://ngdc.cncb.ac.cn/gvm/getProjectDetail?project=GVP000002, GVP000002.
Ethics statement
The studies involving human participants were reviewed and approved by the local institutional ethics review board, the second affiliated hospital of Zhejiang University school of medicine. The patients/participants provided their written informed consent to participate in this study.
Author contributions
YR and MZ conceived and designed the experiments. LL, C-ZL, and YR conducted sample recruitment and quality control. YR, HD, and GQ performed the experiments. HD, LW, JL, and DY undertook data processing and statistical analysis. YR, LL, MZ, and JL interpreted the results. YR, LW, HD, and LL wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the National Nature Science Foundation of China (81872520, 81573040), Basic Public Welfare Research Project of Zhejiang Province (LGF20H110002), Science and Technology program of Yiwu Science and Technology Bureau (Grant No.20-3-243).
Acknowledgments
The authors are most grateful to all the individuals who participated in this study. They would also like to thank Prof. Gang Wang, Jinhua Xu, Jun Gu, Xiaoming Liu, Xiaowen Pang, Hongfu Xie for their assistance in recruited samples.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.968935/full#supplementary-material
References
Arokiasamy, S., Balderstone, M. J. M., De Rossi, G., and Whiteford, J. R. (2019). Syndecan-3 in inflammation and angiogenesis. Front. Immunol. 10, 3031. doi:10.3389/fimmu.2019.03031
Basit, A., Asghar, F., Sadaf, S., and Akhtar, M. W. (2018). Health improvement of human hair and their reshaping using recombinant keratin K31. Biotechnol. Rep. 20, e00288. doi:10.1016/j.btre.2018.e00288
Batalla, A., Coto, E., Gonzalez-Fernandez, D., Gonzalez-Lara, L., Gomez, J., Santos-Juanes, J., et al. (2015). The Cw6 and late-cornified envelope genotype plays a significant role in anti-tumor necrosis factor response among psoriatic patients. Pharmacogenet. Genomics 25 (6), 313–316. doi:10.1097/FPC.0000000000000136
Boehncke, W. H., and Schon, M. P. (2015). Psoriasis. Lancet 386 (9997), 983–994. doi:10.1016/S0140-6736(14)61909-7
Capo, A., Di Nicola, M., Auriemma, M., Piaserico, S., Cuccurullo, C., Santilli, F., et al. (2014). Mean platelet volume variation after biologic therapy in psoriasis and psoriatic arthritis. Eur. J. Dermatol. 24 (1), 133–135. doi:10.1684/ejd.2014.2269
Chang, M., Li, Y., Yan, C., Callis-Duffin, K. P., Matsunami, N., Garcia, V. E., et al. (2008). Variants in the 5q31 cytokine gene cluster are associated with psoriasis. Genes Immun. 9 (2), 176–181. doi:10.1038/sj.gene.6364451
Coto-Segura, P., Gonzalez-Fernandez, D., Batalla, A., Gomez, J., Gonzalez-Lara, L., Queiro, R., et al. (2016). Common and rare CARD14 gene variants affect the antitumour necrosis factor response among patients with psoriasis. Br. J. Dermatol. 175 (1), 134–141. doi:10.1111/bjd.14461
Dand, N., Duckworth, M., Baudry, D., Russell, A., Curtis, C. J., Lee, S. H., et al. (2019). HLA-C*06:02 genotype is a predictive biomarker of biologic treatment response in psoriasis. J. Allergy Clin. Immunol. 143 (6), 2120–2130. doi:10.1016/j.jaci.2018.11.038
Eicher, J. D., Chami, N., Kacprowski, T., Nomura, A., Chen, M. H., Yanek, L. R., et al. (2016). Platelet-related variants identified by exomechip meta-analysis in 157, 293 individuals. Am. J. Hum. Genet. 99 (1), 40–55. doi:10.1016/j.ajhg.2016.05.005
Erb, U., Freyschmidt-Paul, P., and Zoller, M. (2013). Tolerance induction by hair-specific keratins in murine alopecia areata. J. Leukoc. Biol. 94 (4), 845–857. doi:10.1189/jlb.0413196
Fan, Z., Wang, L., Jiang, H., Lin, Y., and Wang, Z. (2020). Platelet dysfunction and its role in the pathogenesis of psoriasis. Dermatology 237, 56–65. doi:10.1159/000505536
Gallo, E., Cabaleiro, T., Roman, M., Solano-Lopez, G., Abad-Santos, F., Garcia-Diez, A., et al. (2013). The relationship between tumour necrosis factor (TNF)-alpha promoter and IL12B/IL-23R genes polymorphisms and the efficacy of anti-TNF-alpha therapy in psoriasis: A case-control study. Br. J. Dermatol. 169 (4), 819–829. doi:10.1111/bjd.12425
Gieger, C., Radhakrishnan, A., Cvejic, A., Tang, W., Porcu, E., Pistis, G., et al. (2011). New gene functions in megakaryopoiesis and platelet formation. Nature 480 (7376), 201–208. doi:10.1038/nature10659
Gottlieb, A. B., Langley, R. G., Strober, B. E., Papp, K. A., Klekotka, P., Creamer, K., et al. (2012). A randomized, double-blind, placebo-controlled study to evaluate the addition of methotrexate to etanercept in patients with moderate to severe plaque psoriasis. Br. J. Dermatol. 167 (3), 649–657. doi:10.1111/j.1365-2133.2012.11015.x
Greb, J. E., Goldminz, A. M., Elder, J. T., Lebwohl, M. G., Gladman, D. D., Wu, J. J., et al. (2016). Psoriasis. Nat. Rev. Dis. Prim. 2, 16082. doi:10.1038/nrdp.2016.82
Honne, K., Hallgrimsdottir, I., Wu, C., Sebro, R., Jewell, N. P., Sakurai, T., et al. (2016). A longitudinal genome-wide association study of anti-tumor necrosis factor response among Japanese patients with rheumatoid arthritis. Arthritis Res. Ther. 18, 12. doi:10.1186/s13075-016-0920-6
Jackson, M., Howie, S. E., Weller, R., Sabin, E., Hunter, J. A., and McKenzie, R. C. (1999). Psoriatic keratinocytes show reduced IRF-1 and STAT-1alpha activation in response to gamma-IFN. FASEB J. 13 (3), 495–502. doi:10.1096/fasebj.13.3.495
Julia, A., Fernandez-Nebro, A., Blanco, F., Ortiz, A., Canete, J. D., Maymo, J., et al. (2016). A genome-wide association study identifies a new locus associated with the response to anti-TNF therapy in rheumatoid arthritis. Pharmacogenomics J. 16 (2), 147–150. doi:10.1038/tpj.2015.31
Kilic, S., Resorlu, H., Isik, S., Oymak, S., Akbal, A., Hiz, M. M., et al. (2017). Association between mean platelet volume and disease severity in patients with psoriasis and psoriatic arthritis. Postepy Dermatol. Alergol. 34 (2), 126–130. doi:10.5114/ada.2017.67076
Lallas, A., Argenziano, G., Zalaudek, I., Apalla, Z., Ardigo, M., Chellini, P., et al. (2016). Dermoscopic hemorrhagic dots: An early predictor of response of psoriasis to biologic agents. Dermatol. Pract. Concept. 6 (4), 7–12. doi:10.5826/dpc.0604a02
Liu, C., Batliwalla, F., Li, W., Lee, A., Roubenoff, R., Beckman, E., et al. (2008). Genome-wide association scan identifies candidate polymorphisms associated with differential response to anti-TNF treatment in rheumatoid arthritis. Mol. Med. 14 (9-10), 575–581. doi:10.2119/2008-00056.Liu
Liu, J., Dong, Z., Zhu, Q., He, D., Ma, Y., Du, A., et al. (2016). TNF-Alpha promoter polymorphisms predict the response to etanercept more powerfully than that to infliximab/adalimumab in spondyloarthritis. Sci. Rep. 6, 32202. doi:10.1038/srep32202
Liu, L., Liu, W., Zheng, Y., Chen, J., Zhou, J., Dai, H., et al. (2019). Anti-TNF- αtreatment-related pathways and biomarkers revealed by transcriptome analysis in Chinese psoriasis patients. BMC Syst. Biol. 13 (2), 29. doi:10.1186/s12918-019-0698-7
Masouri, S., Stefanaki, I., Ntritsos, G., Kypreou, K. P., Drakaki, E., Evangelou, E., et al. (2016). A pharmacogenetic study of psoriasis risk variants in a Greek population and prediction of responses to anti-TNF-alpha and anti-IL-12/23 agents. Mol. Diagn. Ther. 20 (3), 221–225. doi:10.1007/s40291-016-0198-z
Mendrinou, E., Patsatsi, A., Zafiriou, E., Papadopoulou, D., Aggelou, L., Sarri, C., et al. (2016). FCGR3A-V158F polymorphism is a disease-specific pharmacogenetic marker for the treatment of psoriasis with Fc-containing TNFα inhibitors. Pharmacogenomics J. 17, 237–241. doi:10.1038/tpj.2016.16
Najar Nobari, N., Shahidi Dadras, M., Nasiri, S., Abdollahimajd, F., and Gheisari, M. (2020). Neutrophil/platelet to lymphocyte ratio in monitoring of response to TNF-alpha inhibitors in psoriatic patients. Dermatol. Ther. 33, e13457. doi:10.1111/dth.13457
Nishikawa, R., Nagai, H., Bito, T., Ikeda, T., Horikawa, T., Adachi, A., et al. (2016). Genetic prediction of the effectiveness of biologics for psoriasis treatment. J. Dermatol. 43 (11), 1273–1277. doi:10.1111/1346-8138.13412
Ovejero-Benito, M. C., Muñoz-Aceituno, E., Sabador, D., Almoguera, B., Prieto-Pérez, R., Hakonarson, H., et al. (2020). Genome-wide association analysis of psoriasis patients treated with anti-TNF drugs. Exp. Dermatol. 29 (12), 1225–1232. doi:10.1111/exd.14215
Ovejero-Benito, M. C., Prieto-Perez, R., Llamas-Velasco, M., Belmonte, C., Cabaleiro, T., Roman, M., et al. (2017). Polymorphisms associated with etanercept response in moderate-to-severe plaque psoriasis. Pharmacogenomics 18 (7), 631–638. doi:10.2217/pgs-2017-0014
Patterson, A. M., Cartwright, A., David, G., Fitzgerald, O., Bresnihan, B., Ashton, B. A., et al. (2008). Differential expression of syndecans and glypicans in chronically inflamed synovium. Ann. Rheum. Dis. 67 (5), 592–601. doi:10.1136/ard.2006.063875
Qayyum, R., Snively, B. M., Ziv, E., Nalls, M. A., Liu, Y., Tang, W., et al. (2012). A meta-analysis and genome-wide association study of platelet count and mean platelet volume in african americans. PLoS Genet. 8 (3), e1002491. doi:10.1371/journal.pgen.1002491
Rioux, J. D., Daly, M. J., Silverberg, M. S., Lindblad, K., Steinhart, H., Cohen, Z., et al. (2001). Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease. Nat. Genet. 29 (2), 223–228. doi:10.1038/ng1001-223
Rodriguez-Gil, A., Ritter, O., Saul, V. V., Wilhelm, J., Yang, C. Y., Grosschedl, R., et al. (2017). The CCR4-NOT complex contributes to repression of Major Histocompatibility Complex class II transcription. Sci. Rep. 7 (1), 3547. doi:10.1038/s41598-017-03708-7
Sanz-Martinez, M. T., Moga, E., Sanchez Martinez, M. A., Zamora Atenza, C., Vidal, S., Juarez, C., et al. (2020). High levels of platelet-lymphocyte complexes in patients with psoriasis are associated with a better response to anti-TNF-alpha therapy. J. Invest. Dermatol. 140 (6), 1176–1183. doi:10.1016/j.jid.2019.08.457
Schmidt, V. A., Scudder, L., Devoe, C. E., Bernards, A., Cupit, L. D., and Bahou, W. F. (2003). IQGAP2 functions as a GTP-dependent effector protein in thrombin-induced platelet cytoskeletal reorganization. Blood 101 (8), 3021–3028. doi:10.1182/blood-2002-09-2807
Song, G. G., Seo, Y. H., Kim, J. H., Choi, S. J., Ji, J. D., and Lee, Y. H. (2015). Association between TNF-alpha (-308 A/G, -238 A/G, -857 C/T) polymorphisms and responsiveness to TNF-alpha blockers in spondyloarthropathy, psoriasis and crohn's disease: A meta-analysis. Pharmacogenomics 16 (12), 1427–1437. doi:10.2217/pgs.15.90
Suzuki, T., Ito, T., Gilhar, A., Tokura, Y., Reich, K., and Paus, R. (2022). The hair follicle-psoriasis axis: Shared regulatory mechanisms and therapeutic targets. Exp. Dermatol. 31 (3), 266–279. doi:10.1111/exd.14462
Umicevic Mirkov, M., Cui, J., Vermeulen, S. H., Stahl, E. A., Toonen, E. J., Makkinje, R. R., et al. (2013). Genome-wide association analysis of anti-TNF drug response in patients with rheumatoid arthritis. Ann. Rheum. Dis. 72 (8), 1375–1381. doi:10.1136/annrheumdis-2012-202405
van den Reek, Jmpa, Coenen, M. J. H., van de L'Isle Arias, M., Zweegers, J., Rodijk-Olthuis, D., Schalkwijk, J., et al. (2017). Polymorphisms in CD84, IL12B and TNFAIP3 are associated with response to biologics in patients with psoriasis. Br. J. Dermatol. 176 (5), 1288–1296. doi:10.1111/bjd.15005
Vasilopoulos, Y., Manolika, M., Zafiriou, E., Sarafidou, T., Bagiatis, V., Kruger-Krasagaki, S., et al. (2012). Pharmacogenetic analysis of TNF, TNFRSF1A, and TNFRSF1B gene polymorphisms and prediction of response to anti-TNF therapy in psoriasis patients in the Greek population. Mol. Diagn. Ther. 16 (1), 29–34. doi:10.1007/BF03256427
Vuckovic, D., Bao, E. L., Akbari, P., Lareau, C. A., Mousas, A., Jiang, T., et al. (2020). The polygenic and monogenic basis of blood traits and diseases. Cell 182 (5), 1214–1231. e1211. doi:10.1016/j.cell.2020.08.008
Yu, Y., Zeng, D., Ou, Q., Liu, S., Li, A., Chen, Y., et al. (2019). Association of survival and immune-related biomarkers with immunotherapy in patients with non-small cell lung cancer: A meta-analysis and individual patient-level analysis. JAMA Netw. Open 2 (7), e196879. doi:10.1001/jamanetworkopen.2019.6879
Zamora, C., Canto, E., and Vidal, S. (2021). The dual role of platelets in the cardiovascular risk of chronic inflammation. Front. Immunol. 12, 625181. doi:10.3389/fimmu.2021.625181
Zweegers, J., Roosenboom, B., van de Kerkhof, P. C., van den Reek, J. M., Otero, M. E., Atalay, S., et al. (2017). Frequency and predictors of a high clinical response in patients with psoriasis on biological therapy in daily practice: Results from the prospective, multicenter BioCAPTURE cohort. Br. J. Dermatol. 176 (3), 786–793. doi:10.1111/bjd.14888
Keywords: psoriasis, genome-wide association analysis (GWAS), anti-TNF-α agents, treatment response, SNP
Citation: Ren Y, Wang L, Dai H, Qiu G, Liu J, Yu D, Liu J, Lyu C-Z, Liu L and Zheng M (2022) Genome-wide association analysis of anti-TNF-α treatment response in Chinese patients with psoriasis. Front. Pharmacol. 13:968935. doi: 10.3389/fphar.2022.968935
Received: 14 June 2022; Accepted: 28 July 2022;
Published: 19 August 2022.
Edited by:
Yuanting Zheng, Fudan University, ChinaReviewed by:
Xianbo Zuo, China-Japan Friendship Hospital, ChinaHui Shen, Tulane University, United States
Ming Li, Shanghai Jiao Tong University, China
Copyright © 2022 Ren, Wang, Dai, Qiu, Liu, Yu, Liu, Lyu, Liu and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lunfei Liu, 2197055@zju.edu.cn; Min Zheng, minz@zju.edu.cn