Impact of Functional Polymorphisms on Drug Survival of Biological Therapies in Patients with Moderate-to-Severe Psoriasis

Biological therapies (BTs) indicated for psoriasis are highly effective; however, not all patients obtain good results, and loss of effectiveness is the main reason for switching. Genetic factors may be involved. The objective of this study was to evaluate the influence of single-nucleotide polymorphisms (SNPs) on the drug survival of tumor necrosis factor inhibitors (anti-TNF) medications and ustekinumab (UTK) in patients diagnosed with moderate-to-severe psoriasis. We conducted an ambispective observational cohort study that included 379 lines of treatment with anti-TNF (n = 247) and UTK (132) in 206 white patients from southern Spain and Italy. The genotyping of the 29 functional SNPs was carried out using real-time polymerase chain reaction (PCR) with TaqMan probes. Drug survival was evaluated with Cox regression and Kaplan–Meier curves. The multivariate analysis showed that the HLA-C rs12191877-T (hazard ratio [HR] = 0.560; 95% CI = 0.40–0.78; p = 0.0006) and TNF-1031 (rs1799964-C) (HR = 0.707; 95% CI = 0.50–0.99; p = 0.048) polymorphisms are associated with anti-TNF drug survival, while TLR5 rs5744174-G (HR = 0.589; 95% CI = 0.37–0.92; p = 0.02), CD84 rs6427528-GG (HR = 0.557; 95% CI = 0.35–0.88; p = 0.013) and PDE3A rs11045392-T together with SLCO1C1 rs3794271-T (HR = 0.508; 95% CI = 0.32–0.79; p = 0.002) are related to UTK survival. The limitations are the sample size and the clustering of anti-TNF drugs; we used a homogeneous cohort of patients from 2 hospitals only. In conclusion, SNPs in the HLA-C, TNF, TLR5, CD84, PDE3A, and SLCO1C1 genes may be useful as biomarkers of drug survival of BTs indicated for psoriasis, making it possible to implement personalized medicine that will reduce financial healthcare costs, facilitate medical decision-making and improve patient quality of life. However, further pharmacogenetic studies need to be conducted to confirm these associations.


Introduction
Psoriasis is a chronic autoimmune skin disease with a prevalence of 2-4% in Europe and the United States [1]. It is characterized by the development of lesions in the form of erythematous plaques with silvery-white scales, mainly on the scalp, elbows, knees and back [2,3]. Moreover, it is associated with other potentially disabling diseases [4,5].
positions (TNF-238, TNF-308, TNF-857, TNF-1031) or its receptor (TNFRSF1B), or even associated proteins (TNFAIP3, TIRAP). Additionally, SNPs in cytokine genes (IL-1B, IL-6, IL-12B and IL- 17) or for receptors that regulate or trigger this cytokine cascade (CD84,  BCL2,CDKAL1,PGLYRP4-24, FCGR2A, FCGR3A, IL-17RA, IL-23R, SLCO1C1-PDE3A,TLR2,  TLR5, TLR9 and LY96), which leads to hyperstimulation of the immune system and, consequently, excessive proliferation of keratinocytes in the epidermis, typical of psoriatic lesions, have proved to have an impact on the pharmacological response of BTs (Table 1) [28][29][30]. Although we know the possible impact of the above-mentioned biomarkers, it is currently not possible to transfer this information to everyday clinical practice because of the low level of evidence (small sample size studies, where the drugs studied are grouped and the results have not been confirmed with other studies). Therefore, the effect of these SNPs needs to be confirmed in further studies, so that they can be established as predictive biomarkers of response to the BTs indicated for psoriasis [31]. On the basis of all the foregoing, the object of this pharmacogenetic study was to evaluate the influence of polymorphisms in functional genes on the survival of anti-TNF and UTK drugs in patients diagnosed with moderate-to-severe psoriasis.

Patient Characteristics
A total of 206 patients diagnosed with moderate-to-severe psoriasis who had been under treatment with anti-TNF or anti-IL12/23 (UTK) agents for at least 4 months were included in this study. Overall, data on 379 lines of treatment were collected, 247 with anti-TNF (ADA, ETN, INF, CTL) and 132 with UTK. The demographic, clinical and treatment variables are summarized in Table 2. The patients' mean follow-up time was 7.8 years.  Table S1).
The patients treated with UTK were younger (51.53 ± 13.13 years), had been diagnosed earlier (26 years), had a higher body mass index (BMI) (29.26 ± 6.13), and more of them had comorbidities than the other patients (94/132; 71.21%). In addition, the patients treated with UTK had a higher percentage of lesions in difficult-to-treat areas: nails (50%), flexures (37.12%) and genital area (16.28%). However, the bio-naive patients were generally treated with anti-TNF drugs (55.87%) as first-line therapy. Finally, there was more switching with anti-TNF drugs (78.14%), although the patients showed slightly greater adherence than those treated with UTK (80.58% vs. 74.81%) ( Table 2).
Next, we studied the reasons for discontinuation or switching of BTs (Supplementary Table S3). Loss of effectiveness was the main reason for discontinuation (188; 67.14%), followed by adverse events (55; 19.64%). The patients treated with UTK were 1.9 times more likely to interrupt treatment because of lack of effectiveness (OR = 1.895; 95% CI = 1.05-3.52; p = 0.024), whereas those treated with anti-TNF drugs interrupted the therapy because of remission of the disease (OR = 3.835; 95% CI = 1.05-21.07; p = 0.027) ( Table S3).

Anti-TNF
The bivariate analysis of survival of anti-TNF drugs showed an association between the HLA-C and TNF SNPs, specifically TNF-238, TNF-308 and TNF-1031 (Table S9) (Table S9). Furthermore, an influence of the various polymorphisms of the TNF gene on increased survival of anti-TNF drugs was observed (Table S9). Subsequently, the Cox multivariate regression confirmed the association between the TNF-1031 SNP and anti-TNF survival; however, in the TNF-238 and TNF-308 SNPs statistical significance was not main-    (Table S9).

Anti-IL12/23
The bivariate survival analysis showed an association between the polymorphisms in the IL1B, PDE3A, SLCO1C1, CD84 and TLR5 genes, as well as in HLA-C, and longer UTK survival (Table S10). The multivariate confirmed that the PDE3A rs11045392 and SLCO1C1 rs3794271 polymorphisms were associated with longer UTK survival; however, statistical significance was not maintained in the IL1B rs1143627 polymorphism. Finally, the Cox multivariate regression, adjusted for the patient's age and bio-naivety, showed that the T allele of the HLA-C rs12191877 SNP (HR = 0.560; 95% CI = 0.40-0.78; p = 0.0006) and the C allele of the TNF-1031 (rs1799964) SNP (HR = 0.707; 95% CI = 0.50-0.99; p = 0.048) were associated with longer anti-TNF drug survival (Table 3).

Anti-IL12/23
The bivariate survival analysis showed an association between the polymorphisms in the IL1B, PDE3A, SLCO1C1, CD84 and TLR5 genes, as well as in HLA-C, and longer UTK survival (Table S10). The multivariate confirmed that the PDE3A rs11045392 and SLCO1C1 rs3794271 polymorphisms were associated with longer UTK survival; however, statistical significance was not maintained in the IL1B rs1143627 polymorphism.

Discussion
The biological therapies indicated for the treatment of moderate-to-severe psoriasis are highly effective and safe; however, certain patients experience a short-term lack of response and/or a gradual long-term loss of response [15]. This variability in response may be due to genetic factors, because abnormalities in the genes involved in the pathological environment of the disease, metabolism or mechanism of action may influence the effectiveness of BTs. Specifically, polymorphisms in the genes that encode certain human leukocyte antigens, cytokines, receptors and transporters related to the immune system have proved to play a crucial role in interindividual variability in the response to anti-TNF drugs and UTK [25,26,28,31]. In this study, 379 lines of treatment (anti-TNF: 247 and UTK: 132) in 198 patients diagnosed with moderate-to-severe psoriasis were evaluated to determine the influence of SNPs in functional genes on drug survival.
Firstly, we evaluated the survival of the drugs studied and observed that UTK showed longer survival than the anti-TNF drugs (3 vs. 2 years; p Log rank = 0.07; Figure 1; Table S1), the main reason for discontinuation being loss of effectiveness of UTK (OR = 1.895; 95% CI = 1.05-3.52; p = 0.024; Table S2). Similarly, the various systematic reviews and meta-analyses carried out to date have observed longer survival of UTK compared to anti-TNF [20,22].
On the other hand, we evaluated the impact of the SNPs in functional genes selected for this study on the survival of the BTs used in psoriasis, finding that patients carrying the T allele of the HLA-C rs12191877 polymorphism and those with the C allele of the TNF-1031 (rs1799964) polymorphism, adjusted for age and bio-naivety, were associated with longer pharmacological survival of anti-TNF drugs, and that patients who had not developed psoriatic arthritis, and those carrying the G allele of the TLR5 rs5744174 polymorphism, the A allele of the CD84 rs6427528 polymorphism or the T allele of the PDE3A rs11045392 polymorphism, which was in linkage disequilibrium with SLCO1C1 rs3794271, were associated with longer survival of UTK.
The association observed in the patients in this study carrying the T allele of the HLA-C rs12191877 SNP with longer survival of anti-TNF drugs, compared to patients carrying the wild-type homozygous (CC) genotype (T vs. CC; HR = 0.56; 95% CI = 0.40-0.78; p = 0.0006) is in line with previous studies [25]. Human leukocyte antigens (HLAs) form part of the major histocompatibility complex (MHC) and collaborate in identifying exogenous proteins that can trigger an immune response [33,34]. Specifically, the HLA-C glycoprotein promotes the production of cytokines and activates monocytes. High HLA-C levels have been related to the risk of suffering from inflammatory bowel disease or spondyloarthropathy [35]. Consequently, multiple pharmacogenetic studies have been conducted to investigate the impact of HLA genetic variants on the development of psoriasis and its response to the various BTs [25,27,36]. In particular, the HLA-C rs12191877 polymorphism was associated with the response to anti-TNF drugs (a reduction of 75% in psoriasis area and severity index score [PASI 75] at 3 months) in 144 white patients (from Spain) (OR = 0.30, 95% CI = 0.11-0.88, p = 0.027) [37].
Moreover, this SNP is in linkage disequilibrium with the HLA-C*0602 allele (also known as the HLA-Cw*06 allele) [38]. The presence of this allele has proved to confer a high risk of suffering from psoriasis, but its association with the response to anti-TNF drugs has yet to be confirmed [25,26,28,39,40]. In two studies with white patients it was observed that patients with the HLA-Cw*06 allele showed a higher risk of not responding to anti-TNF therapy (109 Spanish patients treated with ADA, ETN or INF, PASI75 in week 24, OR = 58.1, 95% CI = 71.7-93.8, p = 0.049; 1326 British patients treated with ADA, PASI90 at 6 months, OR = 2.95, p < 0.0001) [27,41]. However, there is disagreement with other studies published by Ryan et al., Talamonti et al. and Caldarola et al., because no statistically significant association was found for this allele (p > 0.05) [39,40,42].
Our study observed an association between the C allele of the TNF-1031 (rs1799964) polymorphism of the tumor necrosis factor gene promoter region and survival of anti-TNF (HR = 0.707; 95% CI = 0.50-0.99; p = 0.048, C vs. TT) [28]. The importance of this SNP lies in the proinflammatory cytokine encoded by the gene in question (TNF), which is directly associated with the pathogenesis of psoriatic lesions, as well as with the target of four biological drugs (ADA, ETN, INF and CTL). TNF promotes the production of T cells, which trigger local proliferation of epidermal keratinocytes, and therefore genetic abnormalities of this cytokine can directly influence the functioning of anti-TNF drugs [43][44][45]. However, there is a previous study that evaluated the impact of this polymorphism and that associates its TT genotype with a greater short-term response to anti-TNF ( [41]. With regard to treatment with UTK, it was found that patients carrying the G variant allele of the TLR5 rs5744174 SNP showed longer survival of UTK than those carrying the wild-type AA genotype (G allele vs. AA: HR = 0.589, 95% CI = 0.37-0.92, p = 0.02). These results are in line with the study conducted by Loft and his team on 230 white patients (from Denmark) with psoriasis treated with UTK. It was observed that the patients carrying the TLR5 rs5744174-C allele (complementary to the G allele) showed a better response (∆PASI after 3 months) to UTK treatment (OR = 5.26, 95% CI = 1.93-14.38, p = 0.0012, q = 0.19) [30]. The TLR5 gene, belonging to the toll-like receptor (TLR) family, codes for a transmembrane protein in which the toll interleukin receptor (TIR)-type intracellular domain binds to IL-1, generating an inflammatory cascade when the exterior domain recognizes bacterial flagellin [46,47]. Therefore, when the variant allele is present in the TLR5 rs5744174 polymorphism, IL-1 levels decrease and CCL20 and IFN-γ levels increase, reducing the inflammatory response [48][49][50].
The effect of the CD84 rs6427528 receptor gene polymorphism on the response of autoimmune diseases to BTs has been studied. In particular, a meta-analysis of genomewide association studies (GWAS) with white patients (13 studies/2706 patients) diagnosed with rheumatoid arthritis and treated with anti-TNF showed that the CD84 rs6427528-AG genotype was associated with greater effectiveness of treatment with ETN (733) (p = 0.004) [51]. In psoriasis, a pharmacogenetic study has been conducted in a white population (from the Netherlands) in which the CD84 rs6427528-AG genotype was associated with a better response (∆PASI at 12 weeks) to treatment with ETN (25 patients with 161 ETN episodes; beta = −2.028, 95% CI = −3.794-0.261, p = 0.025) [29]. Our study was unable to confirm this association, but it was observed that the GG genotype of the CD84 rs6427528 polymorphism was associated with UTK survival (GG vs. A: HR = 0.557; 95% CI = 0.35-0.88; p = 0.013). These results are interesting, since genetic abnormalities in the CD84 gene, belonging to the CD2 superfamily of cell surface receptors, even in the non-coding region (3 UTR), mean greater expression of this gene in mononuclear peripheral blood cells, and consequently deregulation of T and B cell signaling, as well as in the adhesion and activation of immune cells [29,52].
Finally, the PDE3A rs11045392 SNP, which is in linkage disequilibrium with SLCO1C1 rs3794271, showed a statistically significant association between the T allele or the TT and CT genotypes, compared to the CC genotype, and UTK survival (T allele vs. CC: HR = 0.508; 95% CI = 0.32-0.79; p = 0.002). In psoriasis, a study evaluated the influence of the SLCO1C1 rs3794271 and PDE3A rs11045392 polymorphisms on anti-TNF response in 130 white patients (from Spain) and showed that patients carrying the C allele obtained a better response (∆PASI after 3 months) (p = 0.00057) [53]. This may be due to the fact that both genes are related to the physiopathogenesis of psoriasis and therefore to the response to various BTs. Specifically, the PDE3A gene is expressed mainly in cardiac tissue and encodes a phosphodiesterase, which is responsible for internal control of nucleotide signaling, while the SLCO1C1 gene encodes a sodium-independent transporter with high affinity for the thyroid hormone in brain tissue.
Among the main limitations of our study that should be highlighted is the limited sample size, which may be responsible for the lack of statistical association found between the SNPs and pharmacological survival. Nevertheless, we obtained important and representative results which show the influence that certain genetic polymorphisms have on the probability of changing treatment with anti-TNF drugs or UTK in patients with moderate-to-severe psoriasis. The fact that the four anti-TNF drugs (ADA, ETN, INF and CTL) were grouped together in the statistical analysis may seem to be a limitation; however, grouping them in this way increases the statistical power, and they are usually studied together. In addition, it should be pointed out that this study was carried out with patients from only two hospitals, all diagnosed with moderate-to-severe psoriasis, following the same therapeutic protocols, by the same team of dermatologists. This makes it possible to achieve a high degree of homogeneity in the cohort, as well as in the process of collecting the variables.
To sum up, the results suggest that the HLA-C rs12191877 and TNF-1031 (rs1799964) SNPs are associated with anti-TNF drug survival, while TLR5 rs5744174, CD84 rs6427528 and PDE3A rs11045392 together with SLCO1C1 rs3794271 are associated with UTK survival. Consequently, these polymorphisms could act as survival biomarkers for the various BTs in patients with moderate-to-severe psoriasis and could be useful for direct selection of the appropriate biological therapy, thereby achieving more individualized medicine with lower pharmaceutical cost and better results for the patient. Although more studies in larger cohorts will be needed to confirm the prognostic value of the biomarkers, the information on the persistence, safety and tolerability of treatment in the short and long term may help patients and healthcare professionals in making decisions to begin treatment with a BT [18].

Study Design
We conducted a prospective and retrospective observational cohort study.

Ethics Statements
This study was carried out with the approval of the Ethics and Research Committee of the Sistema Andaluz de Salud (SAS: Andalusian Health Service) and in accordance with the Declaration of Helsinki (code: FG-DERM-SNP-0702-N-21, 27 October 2015). The patients signed a written informed consent form for collection of saliva samples and their donation to the biobank. The samples were coded and treated confidentially.

Study Population
The study included 260 white patients diagnosed with modern-to-severe psoriasis treated with anti-TNF and UTK drugs, recruited at the Hospital Universitario Virgen de las Nieves, Granada, Spain, and at the Dermatology Clinic, Department of Clinical, Internal, Anesthesiological and Cardiovascular Sciences, La Sapienza Università di Roma, Italy, during the period from January 2019 to November 2021 ( Figure 6). The inclusion criteria for the group of patients were age 18 years or more, diagnosis of moderate-to-severe psoriasis (body surface area [BSA] and/or PASI > 10) in treatment with BT (ADA, INF, ETN, CTL or UTK) and clinical data available. The patients were treated in accordance with the EuroGuiDerm Guideline [54].
Italy, during the period from January 2019 to November 2021 ( Figure 6). The inclusion criteria for the group of patients were age 18 years or more, diagnosis of moderate-tosevere psoriasis (body surface area [BSA] and/or PASI > 10) in treatment with BT (ADA, INF, ETN, CTL or UTK) and clinical data available. The patients were treated in accordance with the EuroGuiDerm Guideline [54].

Socio-Demographic and Clinical Variables
The socio-demographic information, including sex, age and BMI, was collected from the patients' medical records at the beginning of the study. We collected the clinical features of the psoriasis, including age at diagnosis of psoriasis, family history of psoriasis, baseline PASI, type of psoriasis (vulgaris or plaque, pustular, inverse, guttate or a combination of 2, 3 or 4 types) and location of the lesions, according to the consensus document on evaluation and treatment of moderate-to-severe psoriasis: Psoriasis Group of the Spanish Academy of Dermatology and Venereology [3]. We also studied the patients' comorbidities: psoriatic arthritis, hypertension, dyslipidaemia, and others, such as diabetes mellitus, other autoimmune diseases, history of cancer, respiratory diseases (asthma, chronic obstructive pulmonary disease [COPD]), mental health disorders and metabolic syndrome. In addition, we studied variables related to the lines of treatment, such as bionaivety, the line of treatment in general (systemic and biological treatments), the reason for the interruption of the biological therapy, concomitant treatment, therapeutic adherence and duration of biological therapy in months and in days. Therapeutic adherence

Socio-Demographic and Clinical Variables
The socio-demographic information, including sex, age and BMI, was collected from the patients' medical records at the beginning of the study. We collected the clinical features of the psoriasis, including age at diagnosis of psoriasis, family history of psoriasis, baseline PASI, type of psoriasis (vulgaris or plaque, pustular, inverse, guttate or a combination of 2, 3 or 4 types) and location of the lesions, according to the consensus document on evaluation and treatment of moderate-to-severe psoriasis: Psoriasis Group of the Spanish Academy of Dermatology and Venereology [3]. We also studied the patients' comorbidities: psoriatic arthritis, hypertension, dyslipidaemia, and others, such as diabetes mellitus, other autoimmune diseases, history of cancer, respiratory diseases (asthma, chronic obstructive pulmonary disease [COPD]), mental health disorders and metabolic syndrome. In addition, we studied variables related to the lines of treatment, such as bio-naivety, the line of treatment in general (systemic and biological treatments), the reason for the interruption of the biological therapy, concomitant treatment, therapeutic adherence and duration of biological therapy in months and in days. Therapeutic adherence was evaluated with objective measures (dispensing percentage in the biological therapy record) and subjective measures (Morisky-Green and Haynes-Sackett tests).

DNA Isolation
After the patients' inclusion and signing of the informed consent, saliva samples were collected with buccal swabs (OCR-100 Kit from BD, Plymouth, UK). The DNA was extracted using the QIAamp DNA Mini extraction kit (Qiagen GmbH, Hilden, Germany), following the specifications provided by the manufacture for purifying DNA from saliva, and stored at −20 • C. The DNA concentration and purity were measured using a NanoDrop 2000 UV spectrophotometer with the absorbance ratio at 280/260 and 280/230. Table 5 shows the 29 SNPs that were genotyped by real-time polymerase chain reaction (PCR) allelic discrimination assay using TaqMan probes (ABI Applied Biosystems, QuantStudio 3 Real-Time PCR System, 96 wells), following the manufacturer's instructions. In the analysis of HLA-B/MICA and BCL2, customized probes by ThermoFisher Scientific (Waltham, MA, USA) were used. Table 5. Information on the TaqMan SNP genotyping assay used in our study.

Detection of Gene Polymorphisms and Quality Control
rs11465996 C__30755344_10 dbSNP ID: identifier of Single Nucleotide Polymorphism Database. SNPs that were not included in the statistical analysis, because they did not pass quality control, are shaded in grey. * These SNPs were analyzed using custom assays by ThermoFisher Scientific (Waltham, MA, USA).
To assess the internal consistency of the genotyping, 10% of the results were confirmed by Sanger sequencing analysis. The real-time PCR and the Sanger sequencing were performed at the Pharmacogenetics Unit of the Hospital Universitario Virgen de las Nieves, Granada, Spain. Subsequently, quality control to eliminate rare variants and samples with missing data of the genotypes was carried out. The criteria for SNP quality control were (1) missing genotype rate per SNP < 0.05; (2) minor allele frequency > 0.01; (3) p-value > 0.05 in the Hardy-Weinberg equilibrium test; (4) missing genotype rate between cases and control < 0.05.

Survival Variables
The data on treatment with anti-TNF and UTK biological therapies were recorded, including the dates of starting and finishing the therapy and the reasons for discontinuation.
When treatment with a drug is discontinued, a new line of treatment with another drug is initiated. Therefore, the survival of the drug is defined as the time from initiation to discontinuation of treatment with the particular BT. The data were obtained from the patients' medical records.

Statistical Analysis
First, the study sample size was calculated to obtain a precision of 1 unit in the estimate of the logarithm of an odds ratio using a two-sided 95% normal asymptotic confidence interval, based on previous studies.
The patient characteristics were collected and presented as frequencies with percentages for the qualitative variables. The quantitative variables were expressed as mean (plus or minus standard deviation) for the variables with normal distribution or medians and percentiles (25 and 75) for those with non-normal distribution. The Lilliefors (Kolmogorov-Smirnov) test was used to evaluate normality.
We used the Kaplan-Meier survival curve with the log-rank test to compare survival for patients treated with different BTs and the associations of demographic, clinical and genetic variables. The association with SNPs was evaluated in multiple models (genotypic, dominant and recessive), which were defined as follows: genotypic (DD vs. Dd vs. dd), dominant ((DD, Dd) vs. dd) and recessive (DD vs. (Dd, dd)), where D is the major (wildtype) allele and d the minor (variant) allele. Moreover, a univariate and multivariate Cox regression analysis was performed by calculating the hazard ratio (HR) and the 95% confidence interval (95% CI) to predict the risk factors of discontinuing treatment with BTs. For the multivariate analysis, we used the Cox proportional risks regression model, taking into account the variables associated with statistically significant results in the bivariate analysis of drug survival (backward stepwise selection method).
All the tests were two-sided, with a significance level of p < 0.05. The data analysis was performed with the R 4.0.2 program [55]. The Hardy-Weinberg equilibrium and the linkage disequilibrium were determined with the D' and r2 coefficients and estimated using the PLINK, Haploview 4.2 and SNPStat programs [56][57][58].

Conclusions
Our results suggest that patients with moderate-to-severe psoriasis under treatment with TNF inhibitor drugs or UTK show greater pharmacological survival of UTK. Furthermore, various genetic polymorphisms associated with BT survival have been found. Patients carrying the T allele of the HLA-C rs12191877 SNP or the C allele of the TNF-1031 (rs1799964) SNP are associated with greater pharmacological survival of anti-TNF drugs, while patients who have not developed psoriatic arthritis or are carriers of the G allele of the TLR5 rs5744174 polymorphism, the A allele of CD84 rs6427528 or the T allele of PDE3A rs11045392, which is in linkage disequilibrium with SLCO1C1 rs3794271, are associated with longer UTK survival. Consequently, these results confirm that these genetic polymorphisms could be useful as pharmacological survival biomarkers for biological therapies indicated for psoriasis, increasing the current scientific evidence. This information will make it possible to implement personalized medicine that would reduce economic healthcare costs, facilitate medical decision-making and improve patient quality of life. However, further pharmacogenetic studies need to be conducted to confirm these associations.