Journal of Personalized Medicine Article Clinical Value of IL6R Gene Variants as Predictive Biomarkers for Toxicity to Tocilizumab in Patients with Rheumatoid Arthritis Luis Sainz 1,2,3,† , Pau Riera 3,4,5, *,† , Patricia Moya 1,3 , Sara Bernal 3,5,6 , Jordi Casademont 2,3,7 , Cesar Díaz-Torné 1,2,3 , Ana Milena Millán 1,3 , Hye Sang Park 1,3 , Adriana Lasa 3,5,6 and Héctor Corominas 1,2,3, * 1 Rheumatology Department, Hospital de la Santa Creu i Sant Pau (HSCSP), 08041 Barcelona, Spain 2 Faculty of Medicine, Universitat Autònoma de Barcelona (UAB), 08041 Barcelona, Spain 3 Institut d’Investigació Biomèdica Sant Pau (IIB SANT PAU), 08041 Barcelona, Spain 4 Pharmacy Department, Hospital de la Santa Creu i Sant Pau (HSCSP), 08025 Barcelona, Spain 5 CIBER de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029 Madrid, Spain 6 Genetics Department, Hospital de la Santa Creu i Sant Pau (HSCSP), 08025 Barcelona, Spain 7 Medicine Department, Hospital de la Santa Creu i Sant Pau (HSCSP), 08025 Barcelona, Spain * Correspondence: priera@santpau.cat (P.R.); hcorominas@santpau.cat (H.C.) † These authors contributed equally to this work. Abstract: Tocilizumab is a first-line biologic disease-modifying anti-rheumatic drug (bDMARD) that inhibits the interleukin-6 (IL-6) pathway by antagonizing the IL-6 receptor (IL-6R). Tocilizumab is widely used to treat rheumatoid arthritis (RA), a prevalent autoimmune disease that can cause irreversible joint damage and disability. Although many bDMARDs have been developed for RA, there is a lack of validated biomarkers which could guide personalized medicine strategies. To evaluate whether single-nucleotide polymorphisms (SNPs) in the IL6R gene could predict tocilizumab toxicity in patients with RA, we conducted a retrospective cohort study of 88 patients treated with tocilizumab. Six SNPs previously described in the IL6R gene were genotyped (rs12083537, rs11265618, Citation: Sainz, L.; Riera, P.; Moya, P.; rs4329505, rs2228145, rs4537545, and rs4845625). Using parametric tests, we studied the association Bernal, S.; Casademont, J.; Díaz-Torné, between the SNPs and hepatotoxicity, infection, hypersensitivity, gastrointestinal, hematological, and C.; Millán, A.M.; Park, H.S.; Lasa, A.; dyslipidemia adverse events (AEs). We found associations between dyslipidemia and rs4845625 and Corominas, H. Clinical Value of IL6R between hematological AEs and rs11265618 and rs4329505. No further associations were found for Gene Variants as Predictive Biomarkers the remaining SNPs and other AEs. Our findings support the potential clinical value of SNPs in the for Toxicity to Tocilizumab in Patients IL6R gene as predictive biomarkers for toxicity to tocilizumab in patients with RA. with Rheumatoid Arthritis. J. Pers. Med. 2023, 13, 61. https://doi.org/ Keywords: IL6R; genetic variants; tocilizumab; toxicity; rheumatoid arthritis; predictive biomarkers 10.3390/jpm13010061 Academic Editors: Luis Andrés López-Fernández, Sara Salvador-Martín and Matea Zajc 1. Introduction Petranović Rheumatoid arthritis (RA) is the most frequent autoimmune inflammatory disease Received: 5 December 2022 worldwide, with an overall prevalence of 0.2–1% [1,2]. It is characterized by synovial Revised: 22 December 2022 inflammation, hyperplasia, and the progressive destruction of bone and cartilage in multiple Accepted: 23 December 2022 joints [3], which, without appropriate treatment, may lead to progressive disability, systemic Published: 28 December 2022 complications, and reduced life expectancy [4]. The current therapeutic strategy, consisting of early treatment based on a treat-to-target approach, aims to achieve remission or low disease activity and ultimately prevent structural damage [5–10]. At present, there is lack of evidence to recommend the treatment of one biologic Copyright: © 2022 by the authors. disease-modifying anti-rheumatic drug (bDMARD) over another, as there are no validated Licensee MDPI, Basel, Switzerland. biomarkers that predict response or toxicity to bDMARDs. Consequently, treatment selec- This article is an open access article tion is a trial-and-error process, and this can delay the identification of the most effective distributed under the terms and conditions of the Creative Commons and well-tolerated treatment for a patient, which, in turn, may lead to disease progression Attribution (CC BY) license (https:// and/or a poorer quality of life [11–13]. Although there is growing interest in the role creativecommons.org/licenses/by/ of genetic markers as predictors of response/toxicity to biological drugs, most studies 4.0/). J. Pers. Med. 2023, 13, 61. https://doi.org/10.3390/jpm13010061 https://www.mdpi.com/journal/jpm J. Pers. Med. 2023, 13, 61 2 of 9 published to date have focused on methotrexate [14–17] and anti-tumor necrosis factor (anti-TNF) treatments [18–24]. Tocilizumab is a first-line bDMARD that competitively binds to both the membrane- bound and soluble forms of human interleukin-6 receptor (IL-6R), thereafter inhibiting the interleukin-6 (IL-6) pathway. It is an effective treatment for patients with moderate to severe RA, with reported response rates of 62–80% [13,25–27]. Despite being generally well- tolerated, approximately 80% of tocilizumab-treated patients experience adverse events (AEs) [28], most frequently infections, followed by gastrointestinal disorders. Although AEs are mainly mild or moderate, in 4% of patients, they are severe. To date, no biomarkers or predictors of toxicity have been identified. Given that tocilizumab inhibits IL-6R, we hypothesized that functional variations in its encoding gene could affect tocilizumab toxicity. While previous studies have assessed the influence of IL6R gene variants on tocilizumab response, the results are inconsistent and not applicable to clinical practice [28–33]. As for toxicity, no studies to date have evaluated the influence of IL6R gene variants on tocilizumab-induced toxicities. The aim of this study was to assess whether genetic variants in the IL6R gene are associated with tocilizumab toxicity in patients with RA. 2. Materials and Methods 2.1. Study Population We conducted a single-center retrospective cohort study including RA patients re- cruited in a tertiary referral hospital. The potential patients were identified using pharmacy registries in which we filtered for the patients who underwent tocilizumab treatment between January 2016 and January 2021. We applied the following inclusion criteria: RA diagnosed according to the criteria of the American College of Rheumatology (ACR)/European League Against Rheumatism (EULAR) 2010 [34], prescription of tocilizumab for RA, and age over 18 years old. The exclusion criteria were other rheumatic diagnoses such as connective tissue diseases or vasculitis and loss of patient follow-up. Sociodemographic and clinical data collected from electronic medical records were as follows: age, gender, age at diagnosis, previous treatments, comorbidities, body mass index, tocilizumab starting date, dose, and administration route, concomitant treatments, rheumatoid factor (RF), anti-citrullinated protein antibody (ACPA), reason for tocilizumab discontinuation, and AEs. AEs were classified in the following groups: hepatotoxicity, infections, hypersensitivity, gastrointestinal, hematological, and dyslipidemia. For hepatotoxicity, hematological, and dyslipidemia AEs, periodic blood tests were revised retrospectively and, when an AE occurred, numeric values for bilirubinemia, transaminases, leukocytes and neutrophils, platelets, and lipids (total cholesterol, low-density lipoprotein (LDL), and triglycerides) were recorded. AEs were graded as mild (grade 1), moderate (grade 2), severe (grade 3), life- threatening (grade 4), and death (grade 5), in accordance with the Common Terminology Criteria for Adverse Events (CTCAEs) v6.0 [35]. Patients who discontinued tocilizumab without relation to AEs, death, or external causes were classified in the primary failure group, when discontinued in the first 6 months of treatment, and in the secondary failure group, when discontinued after a period of efficacy (>6 months). The study was approved by the corresponding institutional ethics committees and registered at clinicaltrials.gov (protocol code: IIBSP-IIL-2020-148). All the participants gave written informed consent for blood sampling and genetic analyses. 2.2. Genetic Studies Considering published data and bearing in mind single-nucleotide polymorphisms’ (SNPs’) functionality in the IL6R, we selected SNPs as follows: rs12083537, rs11265618, rs4329505, rs2228145, rs4537545, and rs4845625 (Table 1). The rationale for selecting these J. Pers. Med. 2023, 13, 61 3 of 9 SNPs has been described in previous studies [28–31,36]. All the SNPs presented a minor allele frequency (MAF) over 0.10 in the European population according to the Allele Frequency Aggregator (ALFA) [37]. Table 1. Selected functional polymorphisms in the IL6R gene. Genomic Position Genetic Variant MAF Alleles (GRCh38) rs12083537 chr1:154408627 0.21 A>G rs11265618 chr1:154457616 0.17 C>T rs4329505 chr1:154459944 0.17 T>C rs2228145 chr1:154454494 0.40 A>C rs4537545 chr1:154446403 0.41 C>T rs4845625 chr1:154449591 0.43 T>C Abbreviations: IL6R, interleukin-6 receptor. MAF, minor allele frequency reported in the European population. The minor allele is highlighted in bold. Genomic DNA was automatically extracted from peripheral whole-blood samples (Autopure, Qiagen, Hilden, Germany). The SNPs were analyzed via real-time polymerase chain reaction (PCR) using TaqMan SNP genotyping assays (Applied Biosystems, Foster City, CA, USA). All the cases were successfully genotyped. 2.3. Statistical Analysis Hardy–Weinberg equilibrium was assessed for each analyzed SNP using the chi-square test (χ2). Codominant, dominant, and recessive models of inheritance were considered to assess associations with outcome variables whenever appropriate. Quantitative data were expressed as mean (± SD) for normally distributed variables. Normality was assessed using the Shapiro–Wilks test. Depending on the number of groups compared, the Student’s t-test or ANOVA were used for normally distributed variables. Bivariate association between qualitative dichotomous variables was assessed using χ2 or Fisher’s exact test. The associations between SNPs and qualitative response variables were tested using the χ2 test. All the tests were two-sided, significance was set to 5% (α = 0.05), and statistical analyses were performed using IBM-SPSS (IBM Corp. Released 2019. IBM SPSS Statistics for Windows, Version 26.0. Armonk, NY, USA: IBM Corp). 3. Results 3.1. Patient Characteristics Table 2 summarizes the baseline demographic and clinical characteristics for the 88 patients that were included in the study. The mean disease duration on tocilizumab initiation was 16.5 (SD 11.9) years, and treatment survival was 45.91 (SD 40.1) months. Just over a third of patients (35.2%, n =31) were in active treatment at the end of data collection. The primary reason for tocilizumab discontinuation was AEs, followed closely by secondary and primary failure of treatment. 3.2. Adverse Events Table 3 summarizes the frequency and severity of detected AEs, showing that hemato- logical and dyslipidemia AEs were the most frequent. Most hematological AEs were mild and only required observation (n = 16) or dose reduction (n = 9). Neutropenia was the most common toxicity, with a median value of 1.17 (SD 0.03) ×109 /L; however, there were no severe cases, and tocilizumab was discontinued in only two patients. When thrombocytope- nia occurred, the median value of platelets was 108.5 (SD 11.5) ×109 /L. As for dyslipidemia, the mean values of lipid values were as follows: total cholesterol, 284.5 (SD 42.0) mg/dL; LDL, 180.2 (SD 55.7) mg/dL; and triglycerides, 215.3 (SD 144.5) mg/dL. Regarding dys- J. Pers. Med. 2023, 13, 61 4 of 9 lipidemia management, 15 patients were initiated on hypolipemiant treatment during tocilizumab treatment. Table 2. Patient baseline characteristics (n = 88). Variables n (%) Mean (SD) Female/male 76 (86.4) / 12 (13.6) Age on diagnosis (years) 46.6 (15.9) Disease duration (years) 16.5 (11.9) RF positive 58 (65.9) ACPA positive 57 (64.8) Body mass index 28.6 (6.1) Number of previous cDMARDs 2.4 (1.3) Number of previous bDMARDs 1.5 (1.4) Tocilizumab administration (intravenous) 47 (53.4) Treatment survival (months) 45.91 (40.1) Reasons for discontinuation (n =57) Primary failure 14 (15.9) Secondary failure 16 (18.2) AE 19 (21.6) Unrelated death 1 (1.1) Other 7 (8.0) Abbreviations: ACPA, anti-citrullinated protein antibodies; AE, adverse event; bDMARD, biological disease- modifying anti-rheumatic drug; cDMARD, conventional disease-modifying anti-rheumatic drug; RF, rheumatoid factor; SD, standard deviation. Table 3. Reported adverse events (n = 88). Variables n (%) Mild n (%) Moderate n (%) Severe n (%) Hepatotoxicity 20 (22.7) 13 (14.7) 7 (7.9) 0 Hyperbilirubinemia 7 (7.9) Hypertransaminasemia 13 (14.7) Treatment discontinuation 5 (5.7) Infections 18 (20.5) 9 (10.2) 6 (6.8) 3 (3.4) Tocilizumab discontinuation 1 Hypersensitivity 7 (7.9) 0 7 (7.9) 0 Tocilizumab discontinuation 7 Gastrointestinal 3 (3.4) 1 (1.1) 2 (2.3) 0 Tocilizumab discontinuation 0 Hematological 27 (30.7) 19 (21.6) 8 (9.1) 0 Neutropenia 23 (26.1) Thrombocytopenia 6 (6.8) Tocilizumab discontinuation 2 (2.3) Dyslipidemia 27 (30.7) Statin initiation 15 (17) The third most frequent AE was hepatotoxicity, with more cases of hypertransami- nasemia (mean aspartate aminotransferase (AST), 62.7 (SD 21.2) IU/L and alanine amino- transferase (ALT), 92.1 (SD 42.2) IU/L) than of hyperbilirubinemia (mean bilirubin, 24.9 (SD 14.4) µmol/L). Only seven hepatotoxicity AEs were graded as moderate; four of those cases developed during concomitant methotrexate treatment and were managed effectively with methotrexate dose reduction. The most severe AEs were infections, with three patients requiring hospitalization; all the cases were managed with temporary tocilizumab suspension. All seven patients showing hypersensibility to tocilizumab were managed with permanent discontinuation. J. Pers. Med. 2023, 13, 61 5 of 9 3.3. Genetic Determinants and Adverse Events The genotypic frequencies of the six studied SNPs did not significantly deviate from the expected Hardy–Weinberg equilibrium. The mutated allele (C) of the SNP rs4845625 was more frequent than the wild-type (T) allele, as expected from allele frequency population- based studies [37]. In the univariate analysis, three SNPs (rs4845625, rs1126518, and rs432955) showed statistically significant associations with hepatotoxicity, hematological, and dyslipidemia AEs (Table 4). The CC carriers for rs4845625 developed dyslipidemia less frequently (16.7%) than the patients carrying the T allele (TT + CT) (36.7%) (p = 0.04). Considering a codominant model for rs4845625 and hepatotoxicity, statistically significant differences were found between the three genotypes (p = 0.048). Patients carrying the CC genotype for rs11265618 showed a higher rate of hematological AEs compared to T carriers (CT + TT) (36.7% vs. 14.3%, p =0.032). Comparably, the TT genotype for rs4329505 also demonstrated more hematological AEs than C carriers (CT + CC) (36.1% vs. 14.8%, p = 0.044). Table 4. Associations between rs4845625, rs11265618, and rs4329505 and hepatotoxicity, hematologi- cal, and dyslipidemia AEs. Genotype SNPs Hepatotoxicity AEs Hematological AEs Dyslipidemia AEs (n) Genetic OR Genetic OR Genetic OR % p % p % p Model (95% CI) Model (95% CI) Model (95% CI) T/T (13) 7.7 Cod - 0.048 38.5 Cod - 0.712 23.1 Cod - 0.051 rs4845625 C/T (45) 33.3 Dom 4.07 (0.49–33.42) 0.161 26.7 Rec 1.03 (0.39–2.71) 0.946 42.2 Rec 3.05 (1.02 -9.17) 0.04 C/C (30) 13.3 30 16.7 C/C (60) 16 Cod - 0.603 36.7 Cod - 0.075 33.3 Cod - 0.447 rs11265618 C/T (25) 25 Dom 0.65 (0.21–2.02) 0.456 12 Dom 3.47 (1.07–11.36) 0.032 28 Dom 0.67 (0.24–1.83) 0.43 T/T (3) 33.3 33.3 0 T/T (61) 24.6 Cod - 0.665 36.1 Cod - 0.099 32.8 Cod 0.388 rs4329505 C/T (24) 16.7 Dom 0.69 (0.23–2.16) 0.531 12.5 Dom 3.25 (1–10.64) 0.044 29.2 Dom 0.72 (0.26–1.98) 0.52 C/C (3) 33 33.3 0 Abbreviations: Cod: codominant; Dom: dominant; OR: odds ratio; Rec: recessive; SNP: single-nucleotide polymorphism. Further statistically significant associations were not found, either for the remain- ing AEs detected during tocilizumab treatment or for the remaining SNPs (rs12083537, rs2228145, and rs4537545). 4. Discussion We assessed associations between six SNPs in the IL6R gene and AEs during tocilizumab treatment in 88 patients with RA. We found an association between rs4845625 and dyslipi- demia rates and associations between rs11265618 and rs4329505 and hematological AEs. Previous studies focused on UGT1A1 polymorphisms and hyperbilirubinemia have reported results for Japanese patients with RA [38,39]. However, this is the first study demonstrating that SNPs in the IL6R gene may influence the development of AEs during long-term treatment with tocilizumab in patients with RA. Our results regarding rs4845625 show that carriers of the T allele (TT + TC) had a higher incidence of dyslipidemia compared to CC carriers (36.7% vs. 16.7%) when treated with tocilizumab. A recent study of patients with RA by our group also described a greater response to tocilizumab treatment at six months in carriers of the TT or CT genotypes for rs4845625 [33]. Additionally, in line with our results, in 2017, Arguinano et al. described an association between rs4845625*T and high levels of LDL in a healthy cohort, suggesting that IL6R genetic variations may constitute a cardiovascular risk factor [36]. Although it seems that genotypes associated with better efficacy outcomes are also related to greater rates of dyslipidemia, this should not be problematic from a treatment biomarker perspective, given that dyslipidemia secondary to tocilizumab treatment has not been linked to cardiovascular events [40–45]. In fact, the association between the inflammatory burden in RA with changes in the lipid profile and cardiovascular risk has been described J. Pers. Med. 2023, 13, 61 6 of 9 as the “lipid paradox”, as lower cholesterol is associated with higher cardiovascular risk in RA untreated patients [46]. Although the mechanism underlying tocilizumab effects on the lipid profile is unclear, it has been demonstrated that, in a proinflammatory state with the activation of IL-6 pathways, there is a lowering of serum lipids mediated by the modification of apolipoprotein synthesis [47] as well as IL-6 local effects on fat tissue and muscle [44]. It has also been hypothesized that tocilizumab treatment raises lipid levels by reversing the lowering of lipid levels induced by IL-6-mediated inflammation [47,48]. Another suggested explanation is that dyslipidemia is mediated by the direct effect of IL-6 levels, raised through the blockage of IL-6R on administering tocilizumab [44]. We also observed an increase in compounding hematological AEs for the CC genotype for rs11265618 and for the TT genotype for rs4329505. In previous research, these genotypes have been described as potential biomarkers to guide tocilizumab treatment in RA cohorts. CC carriers of rs11265618 showed improved low disease activity rates at 12 months [28], while TT carriers of rs4329505 showed less joint swelling at three months [30]. Although we found similar trends relating to the same genotypes with efficacy outcomes at six months in our previous study, the results were not statistically significant [33]. We found no association between the rs12083537, rs2228145, and rs4537545 genetic variants and AEs detected during tocilizumab treatment. The rs12083537 variant has previously been studied, with conflicting results; while Enevold et al. [30] found that carriers of the AA genotype presented a poorer response to tocilizumab treatment, other studies found the opposite [29,31], and a recent study found no association [33]. As for the IL6R SNP rs2228145, this has only been proposed as a biomarker of response to tocilizumab by Enevold et al. [30]. On the other hand, Arguinano et al. [36] have reported rs4537545*C to be associated with higher LDL levels in a cohort of healthy patients, although only when considered as a haplotype with rs4845625*T. In our research, we did not consider haplotypes, and we did not find increased dyslipidemia for the rs4537545 genetic variant. The main limitation of this study is its retrospective design. The relatively high incidence of AEs reported is possibly explained by our exhaustive search for hepatic, hematologic, and lipid profile abnormalities in routine blood tests over long periods of tocilizumab treatment. This would also explain the predominance of mild AEs that did not require tocilizumab discontinuation. Although the homogeneity of our cohort and the relatively large sample size reinforce the validity of our results, further prospective studies are warranted before the implementation of our findings in clinical practice. 5. Conclusions Our IL6R polymorphism results suggest that, during the treatment of patients with active RA with tocilizumab, rs4845625 is associated with the development of dyslipidemia, and rs11255618 and rs4329505 are associated with the development of hematological effects. Author Contributions: Conceptualization, L.S., P.R., P.M., S.B., C.D.-T., A.L. and H.C.; methodology, L.S., P.R., P.M., S.B., H.S.P., A.L., and H.C.; validation, P.R., P.M., S.B., A.L., J.C., and H.C.; investigation, L.S. and P.R.; resources, S.B. and A.L.; writing—original draft preparation, L.S. and P.R.; writing— review and editing, P.M., S.B., J.C., C.D.-T., A.L., A.M.M., H.S.P., and H.C.; visualization, L.S. and P.R.; supervision, P.M., S.B., A.L., H.C. and J.C.; project administration, L.S. and P.R.; funding acquisition, L.S. All authors have read and agreed to the published version of the manuscript. Funding: This study has been funded by the Societat Catalana de Reumatologia (Beca STADA 2020–2021). Institutional Review Board Statement: The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Research Ethics Committee of the Hospital de la Santa Creu i Sant Pau (IIBSP-IIL-2020-148; 1 March 2021). Informed Consent Statement: Informed consent was obtained from all patients involved in the study. J. Pers. Med. 2023, 13, 61 7 of 9 Data Availability Statement: The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to the inclusion of clinical and personal data. Acknowledgments: We thank Ignasi Gich for statistical analysis, Meritxell Cros for assisting with the genotyping process, and Ailish Maher for language editing. Conflicts of Interest: The authors declare no conflict of interest. References 1. Seoane-Mato, D.; Sánchez-Piedra, C.; Silva-Fernández, L.; Sivera, F.; Blanco, F.J.; Ruiz, F.P.; Juan-Mas, A.; Pego-Reigosa, J.M.; Narváez, J.; Martí, N.Q.; et al. Prevalence of rheumatic diseases in adult population in Spain (EPISER 2016 study): Aims and methodology. Reumatol. Clin. 2015, 15, 90–96. [CrossRef] [PubMed] 2. Cross, M.; Smith, E.; Hoy, D.; Carmona, L.; Wolfe, F.; Vos, T.; Williams, B.; Gabriel, S.; Lassere, M.; Johns, N.; et al. The global burden of rheumatoid arthritis: Estimates from the Global Burden of Disease 2010 study. Ann. Rheum. Dis. 2014, 73, 1316–1322. [CrossRef] 3. McInnes, I.B.; Schett, G. The Pathogenesis of Rheumatoid Arthritis. N. Engl. J. Med. 2011, 365, 2205–2219. [CrossRef] 4. Rodríguez-Valverde, V.; Cáliz, R.C.; Álvaro, J.M.-G.; de la Fuente, J.L.M.; Mendoza, J.M.; Molina, J.T.; Sánchez, J.A.; García, F.B.; Gualda, E.B.; Crespillo, J.C.; et al. III Actualización del Consenso de la Sociedad Española de Reumatología sobre terapia biológica en la artritis reumatoide. Reumatol. Clin. 2006, 2 (Suppl. S2), S52–S59. [CrossRef] [PubMed] 5. Goekoop-Ruiterman, Y.P.; De Vries-Bouwstra, J.K.; Allaart, C.F.; Van Zeben, D.; Kerstens, P.J.; Hazes, J.M.W.; Zwinderman, A.H.; Peeters, A.J.; De Jonge-Bok, J.M.; Mallée, C.; et al. Comparison of Treatment Strategies in Early Rheumatoid Arthritis: A randomized trial. Ann. Intern. Med. 2007, 146, 406–415. [CrossRef] [PubMed] 6. Saunders, S.A.; Capell, H.A.; Stirling, A.; Vallance, R.; Kincaid, W.; McMahon, A.D.; Porter, D.R. Triple therapy in early active rheumatoid arthritis: A randomized, single-blind, controlled trial comparing step-up and parallel treatment strategies. Arthritis Rheum. 2008, 58, 1310–1317. [CrossRef] 7. Möttönen, T.; Hannonen, P.; Korpela, M.; Nissilä, M.; Kautiainen, H.; Ilonen, J.; Laasonen, L.; Kaipiainen-Seppänen, O.; Franzen, P.; Helve, T.; et al. Delay to institution of therapy and induction of remission using single-drug or combination-disease-modifying antirheumatic drug therapy in early rheumatoid arthritis. Arthritis Rheum. 2002, 46, 894–898. [CrossRef] 8. Verstappen, S.M.M.; Jacobs, J.W.G.; van der Veen, M.J.; Heurkens, A.H.M.; Schenk, Y.; ter Borg, E.J.; Blaauw, A.A.M.; Bijlsma, J.W.J.; on the behalf of the Utrecht Rheumatoid Arthritis Cohort study group Intensive treatment with methotrexate in early rheumatoid arthritis: Aiming for remission. Computer Assisted Management in Early Rheumatoid Arthritis (CAMERA, an open-label strategy trial). Ann. Rheum. Dis. 2007, 66, 1443–1449. [CrossRef] 9. Grigor, C.; Capell, H.; Stirling, A.; McMahon, A.D.; Lock, P.; Vallance, R.; Porter, D.; Kincaid, W. Effect of a treatment strategy of tight control for rheumatoid arthritis (the TICORA study): A single-blind randomised controlled trial. Lancet 2004, 364, 263–269. [CrossRef] 10. Smolen, J.S.; Landewé, R.B.M.; Bijlsma, J.W.J.; Burmester, G.R.; Dougados, M.; Kerschbaumer, A.; McInnes, I.B.; Sepriano, A.; van Vollenhoven, R.F.; de Wit, M.; et al. EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2019 update. Ann. Rheum. Dis. 2020, 79, 685–699. [CrossRef] 11. Pers, Y.-M.; Fortunet, C.; Constant, E.; Lambert, J.; Godfrin-Valnet, M.; De Jong, A.; Mercier, G.; Prades, B.P.; Wendling, D.; Gaudin, P.; et al. Predictors of response and remission in a large cohort of rheumatoid arthritis patients treated with tocilizumab in clinical practice. Rheumatology 2013, 53, 76–84. [CrossRef] [PubMed] 12. Tarnowski, M.; Paradowska-Gorycka, A.; Dabrowska-Zamojcin, ˛ E.; Czerewaty, M.; Słuczanowska-Głabowska, ˛ S.; Pawlik, A. The effect of gene polymorphisms on patient responses to rheumatoid arthritis therapy. Expert Opin. Drug Metab. Toxicol. 2016, 12, 41–55. [CrossRef] [PubMed] 13. Narváez, J.; Magallares, B.; Torné, C.D.; Hernández, M.V.; Reina, D.; Corominas, H.; Sanmartí, R.; Llobet, J.M.; de la Serna, A.R.; Nolla, J.M. Predictive factors for induction of remission in patients with active rheumatoid arthritis treated with tocilizumab in clinical practice. Semin. Arthritis Rheum. 2016, 45, 386–390. [CrossRef] [PubMed] 14. Cao, M.; Guo, M.; Wu, D.-Q.; Meng, L. Pharmacogenomics of Methotrexate: Current Status and Future Outlook. Curr. Drug Metab. 2018, 19, 1182–1187. [CrossRef] 15. Salazar, J.; Moya, P.; Altés, A.; Díaz-Torné, C.; Casademont, J.; Cerdà-Gabaroi, D.; Corominas, H.; Baiget, M. Polymorphisms in genes involved in the mechanism of action of methotrexate: Are they associated with outcome in rheumatoid arthritis patients? Pharmacogenomics 2014, 15, 1079–1090. [CrossRef] 16. Moya, P.; Salazar, J.; Arranz, M.J.; Díaz-Torné, C.; del Río, E.; Casademont, J.; Corominas, H.; Baiget, M. Methotrexate phar- macokinetic genetic variants are associated with outcome in rheumatoid arthritis patients. Pharmacogenomics 2016, 17, 25–29. [CrossRef] 17. Jekic, B.; Maksimovic, N.; Damnjanovic, T.; Biljana, J.; Nela, M.; Tatjana, D. Methotrexate pharmacogenetics in the treatment of rheumatoid arthritis. Pharmacogenomics 2019, 20, 1235–1245. [CrossRef] J. Pers. Med. 2023, 13, 61 8 of 9 18. Márquez, A.; Ferreiro-Iglesias, A.; Dávila-Fajardo, C.L.; Montes, A.; Pascual-Salcedo, D.; Perez-Pampin, E.; Moreno-Ramos, M.J.; García-Portales, R.; Navarro, F.; Moreira, V.; et al. Lack of validation of genetic variants associated with anti–tumor necrosis factor therapy response in rheumatoid arthritis: A genome-wide association study replication and meta-analysis. Arthritis Res. Ther. 2014, 16, R66. [CrossRef] 19. Sode, J.; Vogel, U.; Bank, S.; Andersen, P.S.; Hetland, M.L.; Locht, H.; Heegaard, N.H.H.; Andersen, V. Genetic Variations in Pattern Recognition Receptor Loci Are Associated with Anti-TNF Response in Patients with Rheumatoid Arthritis. PLoS ONE 2015, 10, e0139781. [CrossRef] 20. Bek, S.; Bojesen, A.B.; Nielsen, J.V.; Sode, J.; Bank, S.; Vogel, U.; Andersen, V. Systematic review and meta-analysis: Pharmacoge- netics of anti-TNF treatment response in rheumatoid arthritis. Pharmacogenomics J. 2017, 17, 403–411. [CrossRef] 21. Wysocki, T.; Paradowska-Gorycka, A. Pharmacogenomics of Anti-TNF Treatment Response Marks a New Era of Tailored Rheumatoid Arthritis Therapy. Int. J. Mol. Sci. 2022, 23, 2366. [CrossRef] [PubMed] 22. Montes, A.; Perez-Pampin, E.; Joven, B.; Carreira, P.; Fernández-Nebro, A.; Ordóñez, M.D.C.; Navarro-Sarabia, F.; Moreira, V.; Vasilopoulos, Y.; Sarafidou, T.; et al. FCGR polymorphisms in the treatment of rheumatoid arthritis with Fc-containing TNF inhibitors. Pharmacogenomics 2015, 16, 333–345. [CrossRef] [PubMed] 23. Chen, W.; Xu, H.; Wang, X.; Gu, J.; Xiong, H.; Shi, Y. The tumor necrosis factor receptor superfamily member 1B polymorphisms predict response to anti-TNF therapy in patients with autoimmune disease: A meta-analysis. Int. Immunopharmacol. 2015, 28, 146–153. [CrossRef] [PubMed] 24. Canet, L.M.; Filipescu, I.; Cáliz, R.; Lupiañez, C.B.; Canhão, H.; Escudero, A.; Segura-Catena, J.; Soto-Pino, M.J.; Ferrer, M.A.; García, A.; et al. Genetic variants within the TNFRSF1B gene and susceptibility to rheumatoid arthritis and response to anti-TNF drugs: A multicenter study. Pharmacogenetics Genom. 2015, 25, 323–333. [CrossRef] 25. Iking-Konert, C.; Von Hinüber, U.; Richter, C.; Schwenke, H.; Gürtler, I.; Kästner, P.; Klapperich, B.; Peters, M.; Burmester, G.-R. ROUTINE—A prospective, multicentre, non-interventional, observational study to evaluate the safety and effectiveness of intravenous tocilizumab for the treatment of active rheumatoid arthritis in daily practice in Germany Rheumatology key messages. Rheumatology 2016, 55, 624–635. [CrossRef] 26. Forsblad-D’Elia, H.; Bengtsson, K.; Kristensen, L.E.; Jacobsson, L.T.H. Drug adherence, response and predictors thereof for tocilizumab in patients with rheumatoid arthritis: Results from the Swedish biologics register. Rheumatology 2014, 54, 1186–1193. [CrossRef] 27. Backhaus, M.; Kaufmann, J.; Richter, C.; Wassenberg, S.; Roske, A.-E.; Hellmann, P.; Gaubitz, M. Comparison of tocilizumab and tumour necrosis factor inhibitors in rheumatoid arthritis: A retrospective analysis of 1603 patients managed in routine clinical practice. Clin. Rheumatol. 2015, 34, 673–681. [CrossRef] 28. Maldonado-Montoro, M.; Cañadas-Garre, M.; González-Utrilla, A.; Calleja-Hernández, M. Influence of IL6R gene polymorphisms in the effectiveness to treatment with tocilizumab in rheumatoid arthritis. Pharmacogenomics J. 2018, 18, 167–172. [CrossRef] 29. Maldonado-Montoro, M.; Cañadas-Garre, M.; González-Utrilla, A.; Plaza-Plaza, J.C.; Calleja-Hernández, M. Genetic and clinical biomarkers of tocilizumab response in patients with rheumatoid arthritis. Pharmacol. Res. 2016, 111, 264–271. [CrossRef] 30. Enevold, C.; Baslund, B.; Linde, L.; Josephsen, N.L.; Tarp, U.; Lindegaard, H.; Jacobsen, S.; Nielsen, C.H. Interleukin-6- receptor polymorphisms rs12083537, rs2228145, and rs4329505 as predictors of response to tocilizumab in rheumatoid arthritis. Pharmacogenetics Genom. 2014, 24, 401–405. [CrossRef] 31. Luxembourger, C.; Ruyssen-Witrand, A.; Ladhari, C.; Rittore, C.; Degboe, Y.; Maillefert, J.-F.; Gaudin, P.; Marotte, H.; Wendling, D.; Jorgensen, C.; et al. A single nucleotide polymorphism of IL6-receptor is associated with response to tocilizumab in rheumatoid arthritis patients. Pharmacogenomics J. 2019, 19, 368–374. [CrossRef] 32. Wang, J.; Bansal, A.T.; Martin, M.A.; Germer, S.; Benayed, R.; Essioux, L.; Lee, J.S.; Begovich, A.B.; Hemmings, A.; Kenwright, A.; et al. Genome-wide association analysis implicates the involvement of eight loci with response to tocilizumab for the treatment of rheumatoid arthritis. Pharmacogenomics J. 2013, 13, 235–241. [CrossRef] 33. Sainz, L.; Riera, P.; Moya, P.; Bernal, S.; Casademont, J.; Díaz-Torné, C.; Millán, A.M.; Park, H.S.; Lasa, A.; Corominas, H. Role of IL6R Genetic Variants in Predicting Response to Tocilizumab in Patients with Rheumatoid Arthritis. Pharmaceutics 2022, 14, 1942. [CrossRef] 34. Aletaha, D.; Neogi, T.; Silman, A.J.; Funovits, J.; Felson, D.T.; Bingham, C.O., 3rd; Birnbaum, N.S.; Burmester, G.R.; Bykerk, V.P.; Cohen, M.D.; et al. 2010 Rheumatoid arthritis classification criteria: An American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum 2010, 62, 2569–2581. [CrossRef] 35. Common Terminology Criteria for Adverse Events (CTCAE) Version 5; US Department of Health and Human Services, National Institutes of Health, National Cancer Institute: Bethesda, MD, USA, 2017. 36. Arguinano, A.; Naderi, E.; Ndiaye, N.C.; Stathopoulou, M.; Dadé, S.; Alizadeh, B.; Visvikis-Siest, S. IL6R haplotype rs4845625*T/rs4537545*C is a risk factor for simultaneously high CRP, LDL and ApoB levels. Genes Immun. 2017, 18, 163–169. [CrossRef] 37. Phan, L.; Jin, Y.; Zhang, H.; Qiang, W.; Shekhtman, E.; Shao, D.; Revoe, D.; Villamarin, R.; Ivanchenko, E.; Kimura, M.; et al. ALFA: Allele Frequency Aggregator. National Center for Biotechnology Information, U.S. National Library of Medicine. 2020. Available online: www.ncbi.nlm.nih.gov/snp/docs/gsr/alfa/ (accessed on 10 March 2020). J. Pers. Med. 2023, 13, 61 9 of 9 38. Lee, J.S.; Wang, J.; Martin, M.; Germer, S.; Kenwright, A.; Benayed, R.; Spleiss, O.; Platt, A.; Pilson, R.; Hemmings, A.; et al. Genetic variation in UGT1A1 typical of Gilbert syndrome is associated with unconjugated hyperbilirubinemia in patients receiving tocilizumab. Pharmacogenetics Genom. 2011, 21, 365–374. [CrossRef] 39. Mori, S.; Terada, K.; Ueki, Y. Tocilizumab-induced hyperbilirubinemia in Japanese patients with rheumatoid arthritis: Its association with UDP glucuronosyltransferase 1A1 gene polymorphisms. Mod. Rheumatol. 2011, 22, 515–523. [CrossRef] 40. Hoffman, E.; Rahat, M.A.; Feld, J.; Elias, M.; Rosner, I.; Kaly, L.; Lavie, I.; Gazitt, T.; Zisman, D. Effects of Tocilizumab, an Anti-Interleukin-6 Receptor Antibody, on Serum Lipid and Adipokine Levels in Patients with Rheumatoid Arthritis. Int. J. Mol. Sci. 2019, 20, 4633. [CrossRef] 41. Lukas, C.; Redondin, M.; Pane, I.; Soubrier, M.; Houvenagel, E.; Sibilia, J.; Combe, B.; Morel, J. Cardiovascular events and change in cholesterol levels in patients with rheumatoid arthritis treated with tocilizumab: Data from the REGATE Registry. Clin. Exp. Rheumatol. 2021, 39, 501–507. [CrossRef] 42. Castagné, B.; Viprey, M.; Martin, J.; Schott, A.-M.; Cucherat, M.; Soubrier, M. Cardiovascular safety of tocilizumab: A systematic review and network meta-analysis. PLoS ONE 2019, 14, e0220178. [CrossRef] 43. Kim, S.C.; Solomon, D.H.; Rogers, J.R.; Gale, S.; Klearman, M.; Sarsour, K.; Schneeweiss, S. Cardiovascular Safety of Tocilizumab Versus Tumor Necrosis Factor Inhibitors in Patients with Rheumatoid Arthritis: A Multi-Database Cohort Study. Arthritis Rheumatol. 2017, 69, 1154–1164. [CrossRef] 44. Alsulaim, T.; Alhassan, N.; Khalil, H.; Almutlaq, A. Tocilizumab Effect on Lipid Profile in Correlation to Cardiovascular Events: A Retrospective Cohort Study. Int. J. Rheumatol. 2021, 2021, 5535486. [CrossRef] 45. Giles, J.T.; Sattar, N.; Gabriel, S.; Ridker, P.M.; Gay, S.; Warne, C.; Musselman, D.; Brockwell, L.; Shittu, E.; Klearman, M.; et al. Cardiovascular Safety of Tocilizumab Versus Etanercept in Rheumatoid Arthritis: A Randomized Controlled Trial. Arthritis Rheumatol. 2020, 72, 31–40. [CrossRef] 46. Behl, T.; Kaur, I.; Sehgal, A.; Zengin, G.; Brisc, C.; Brisc, M.C.; Munteanu, M.A.; Nistor-Cseppento, D.C.; Bungau, S. The Lipid Paradox as a Metabolic Checkpoint and Its Therapeutic Significance in Ameliorating the Associated Cardiovascular Risks in Rheumatoid Arthritis Patients. Int. J. Mol. Sci. 2020, 21, 9505. [CrossRef] 47. Hashizume, M.; Yoshida, H.; Koike, N.; Suzuki, M.; Mihara, M. Overproduced interleukin 6 decreases blood lipid levels via upregulation of very-low-density lipoprotein receptor. Ann. Rheum. Dis. 2010, 69, 741–746. [CrossRef] 48. Hashizume, M.; Mihara, M. IL-6 and lipid metabolism. Inflamm. Regen. 2011, 31, 325–333. [CrossRef] Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.