PriME-PGx: La Princesa University Hospital Multidisciplinary Initiative for the Implementation of Pharmacogenetics

SLCO1B1 and ABCG2 genotype-informed phenotypes are related to variation in ramipril exposure

Ramipril is an angiotensin-converting enzyme inhibitor used for hypertension and heart failure management. To date, scarce literature is available on pharmacogenetic associations affecting ramipril. The goal of this study was to investigate the effect of 120 genetic variants in 34 pharmacogenes (i.e., genes encoding for enzymes like CYPs or UGTs and transporters like ABC or SLC) on ramipril pharmacokinetic variability and adverse drug reaction (ADR) incidence. Twenty-nine healthy volunteers who had participated in a single-dose bioequivalence clinical trial of two formulations of ramipril were recruited. A univariate and multivariate analysis searching for associations between genetic variants and ramipril pharmacokinetics was performed. SLCO1B1 and ABCG2 genotype-informed phenotypes strongly predicted ramipril exposure. Volunteers with the SLCO1B1 decreased function (DF) phenotype presented around 1.7-fold higher dose/weight-corrected area under the curve (AUC/DW) than volunteers with the normal function (NF) phenotype (univariate p-value [puv] < 0.001, multivariate p-value [pmv] < 0.001, β = 0.533, R2 = 0.648). Similarly, volunteers with ABCG2 DF + poor function (PF) phenotypes presented around 1.6-fold higher AUC/DW than those with the NF phenotype (puv = 0.011, pmv < 0.001, β = 0.259, R2 = 0.648). Our results suggest that SLCO1B1 and ABCG2 are important transporters to ramipril pharmacokinetics, and their genetic variation strongly alters its pharmacokinetics. Further studies are required to confirm these associations and their clinical relevance.

Impact of STAT6 Variants on the Response to Proton Pump Inhibitors and Comorbidities in Patients with Eosinophilic Esophagitis

Proton pump inhibitors (PPIs) are the first-line drug for eosinophilic esophagitis (EoE), although it is estimated that there is a lack of histological remission in 50% of patients. This research aimed to identify pharmacogenetic biomarkers predictive of PPI effectiveness and to study their association with disease features. Peak eosinophil count (PEC) and the endoscopic reference score (EREFS) were determined before and after an eight-week PPI course in 28 EoE patients. The impact of the signal transducer and activator of transcription 6 (STAT6), CYP2C19, CYP3A4, CYP3A5, and ABCB1 genetic variations on baseline PEC and EREFS, their reduction and histological response, and on EoE symptoms and comorbidities was analyzed. PEC reduction was higher in omeprazole-treated patients (92.5%) compared to other PPIs (57.9%, p = 0.003). STAT6 rs12368672 (g.18453G>C) G/G genotype showed higher baseline PEC values compared to G/C and C/C genotypes (83.2 vs. 52.9, p = 0.027). EREFS reduction in STAT6 rs12368672 G/G and G/C genotypes was higher than in the C/C genotype (36.7% vs. -75.0% p = 0.011). However, significance was lost after Bonferroni correction. Heartburn incidence was higher in STAT6 rs167769 (g.27148G>A) G/G patients compared to G/A (54.55% vs. 11.77%, p = 0.030). STAT6 rs12368672G>C and rs167769G>A variants might have a relevant impact on EoE status and PPI response. Further research is warranted to clarify the clinical relevance of these variants.

Developments in pharmacogenetics, pharmacogenomics, and personalized medicine

The development of Pharmacogenetics and Pharmacogenomics in Western Europe is highly relevant in the worldwide scenario. Despite the usually low institutional support, many research groups, composed of basic and clinical researchers, have been actively working for decades in this field. Their contributions made an international impact and paved the way for further studies and pharmacogenomics implementation in clinical practice. In this manuscript, that makes part of the Special Issue entitled Spanish Pharmacology, we present an analysis of the state of the art of Pharmacogenetics and Pharmacogenomics research in Europe, we compare it with the developments in Spain, and we summarize the most salient contributions since 1988 to the present, as well as recent developments in the clinical application of pharmacogenomics knowledge. Finally, we present some considerations on how we could improve translation to clinical practice in this specific scenario.

Impact of CYP2D6 and CYP2B6 phenotypes on the response to tramadol in patients with acute post-surgical pain

Tramadol is an important minor opioid prescribed for pain management. In this study, we analyzed the well-known impact of CYP2D6 genetic variation and 60 additional variants in eight candidate genes (i.e., ABCG2, SLCO1B1, CYP2D6, CYP2B6, CYP2C19, CYP2C9, CYP3A5, and CYP3A4) on tramadol efficacy and safety. Some 108 patients with pain after surgery admitted to a post-anesthesia care unit (PACU) and prescribed tramadol were recruited. They were genotyped, and tramadol M1/M2 metabolite concentrations were determined by a newly validated HPLC-MS/MS method. CYP2D6 intermediate (IM) and poor (PM) metabolizers showed lower M1 concentrations adjusted for dose/weight at 30 and 120 min compared to ultrarapid (UM) and normal (NM) metabolizers (univariate p < 0.001 and 0.020, multivariate p < 0.001 and 0.001, unstandardized β coefficients = 0.386 and 0.346, R2  = 0.146 and 0.120, respectively). CYP2B6 PMs (n = 10) were significantly related to a higher reduction in pain 30 min after tramadol intake (univariate p = 0.038, multivariate p = 0.016, unstandardized β coefficient = 0.224, R2  = 0.178), to lower PACU admission time (p = 0.007), and to lower incidence of adverse drug reactions (p = 0.038) compared to the other phenotypes. CYP3A4 IMs and PMs showed a higher prevalence of drowsiness and dizziness (p = 0.028 and 0.005, respectively). Our results suggest that the interaction of CYP2B6 and CYP2D6 phenotypes may be clinically relevant, pending validation of these results in large, independent cohorts. Additional research is required to clarify the impact of CYP3A4 genetic variation on tramadol response.

Clinical pharmacology facing the real-world setting: Pharmacovigilance, pharmacoepidemiology and the economic evaluation of drugs

Traditionally, clinical pharmacology has focused its activities on drug-organism interaction, from an individual or collective perspective. Drug efficacy assessment by performing randomized clinical trials and analysis of drug use in clinical practice by carrying out drug utilization studies have also been other areas of interest. From now on, Clinical pharmacology should move from the analysis of the drug-individual interaction to the analysis of the drug-individual-society interaction. It should also analyze the clinical and economic consequences of the use of drugs in the conditions of normal clinical practice, beyond clinical trials. The current exponential technological development that facilitates the analysis of real-life data offers us a golden opportunity to move to all these other areas of interest. This review describes the role that clinical pharmacology has played at the beginning and during the evolution of pharmacovigilance, pharmacoepidemiology and economic drug evaluations in Spain. In addition, the challenges that clinical pharmacology is going to face in the following years in these three areas are going to be outlined too.

Impact of Sex and Genetic Variation in Relevant Pharmacogenes on the Pharmacokinetics and Safety of Valsartan, Olmesartan and Hydrochlorothiazide

Drug combination therapy is the most common pharmacological strategy for hypertension management. No pharmacogenetic biomarkers for guiding hypertension pharmacotherapy are available to date. The study population were 64 volunteers from seven bioequivalence trials investigating formulations with valsartan, olmesartan and/or hydrochlorothiazide. Every volunteer was genotyped for 10 genetic variants in different transporters' genes. Additionally, valsartan-treated volunteers were genotyped for 29 genetic variants in genes encoding for different metabolizing enzymes. Variability in pharmacokinetic parameters such as maximum concentration (Cmax) and time to reach it (tmax), the incidence of adverse drug reactions (ADRs) and blood pressure measurements were analyzed as a function of pharmacogenetic and demographic parameters. Individuals with the ABCB1 rs1045642 T/T genotype were associated with a higher valsartan tmax compared to those with T/G and G/G genotypes (p < 0.001, β = 0.821, R2 = 0.459) and with a tendency toward a higher postural dizziness incidence (11.8% vs. 0%, p = 0.070). A higher hydrochlorothiazide dose/weight (DW)-corrected area under the curve (AUC∞/DW) was observed in SLC22A1 rs34059508 G/A volunteers compared to G/G volunteers (p = 0.050, β = 1047.35, R2 = 0.051), and a tendency toward a higher postural dizziness incidence (50% vs. 1.6%, p = 0.063). Sex impacted valsartan and hydrochlorothiazide pharmacokinetics, showing a lower exposure in women, whereas no significant differences were found for olmesartan pharmacokinetics.

NAT2 phenotype alters pharmacokinetics of rivaroxaban in healthy volunteers

Rivaroxaban is a direct inhibitor of factor Xa, a member of direct oral anticoagulant group of drugs (DOACs). Despite being a widely extended alternative to vitamin K antagonists (i.e., acenocoumarol, warfarin) the interindividual variability of DOACs is significant, and may be related to adverse drug reaction occurrence or drug inefficacy, namely hemorrhagic or thromboembolic events. Since there is not a consistent analytic practice to monitor the anticoagulant activity of DOACs, previously reported polymorphisms in genes coding for proteins responsible for the activation, transport, or metabolism of DOACs were studied. The study population comprised 60 healthy volunteers, who completed two randomized, crossover bioequivalence clinical trials between two different rivaroxaban formulations. The effect of food, sex, biogeographical origin and 55 variants (8 phenotypes and 47 single nucleotide polymorphisms) in drug metabolizing enzyme genes (such as CYP2D6, CYP2C9, NAT2) and transporters (namely, ABCB1, ABCG2) on rivaroxaban pharmacokinetics was tested. Individuals dosed under fasting conditions presented lower tmax (2.21 h vs 2.88 h, β = 1.19, R2 =0.342, p = 0.012) compared to fed volunteers. NAT2 slow acetylators presented higher AUC∞ corrected by dose/weight (AUC∞/DW; 8243.90 vs 7698.20 and 7161.25 h*ng*mg /ml*kg, β = 0.154, R2 =0.250, p = 0.044), higher Cmax/DW (1070.99 vs 834.81 and 803.36 ng*mg /ml*kg, β = 0.245, R2 =0.320, p = 0.002), and lower tmax (2.63 vs 3.19 and 4.15 h, β = -0.346, R2 =0.282, p = 0.047) than NAT2 rapid and intermediate acetylators. No other association was statistically significant. Thus, slow NAT2 appear to have altered rivaroxaban pharmacokinetics, increasing AUC∞ and Cmax. Nonetheless, further research should be conducted to verify NAT2 involvement on rivaroxaban pharmacokinetics and to determine its clinical significance.

Preemptive TPMT Genotyping and Adherence to Genotype-Based Therapeutic Recommendations Reduces the Healthcare Cost in Patients Receiving Azathioprine or 6-Mercaptopurine for Autoimmune Diseases

A cost analysis of thiopurine treatment was carried out in 257 patients, with 153 preemptively genotyped for TPMT and 104 retrospectively genotyped in a Spanish setting. The healthcare cost was significantly higher in patients retrospectively genotyped compared to those who were preemptively genotyped (p < 0.001). TPMT intermediate metabolizers (IMs) (n = 23) showed a 3.3-fold higher healthcare cost when compared to normal metabolizers (NMs) (p < 0.001). The healthcare cost in patients with a TPMT IM phenotype whose physician adhered to the genotype-informed recommendation was similar than the cost in TPMT NMs and was significantly lower than IMs whose physician did not adhere to the therapeutic recommendation (3.8-fold, p = 0.016). Myelotoxicity occurrence was significantly lower in patients preemptively vs. retrospectively genotyped (2.0% and 21.2%, respectively, p < 0.001). Patients who developed myelotoxicity showed a significantly higher healthcare cost than those who did not (4.10-fold, p < 0.001). Overall, 87% of patients whose dose was not adjusted despite being TPMT IMs suffered myelotoxicity, while only one of the eight patients (13%) whose dose was adjusted suffered myelotoxicity (p < 0.001). In conclusion, TPMT preemptive genotyping and physician adherence to genotype-informed therapeutic recommendations prevents myelotoxicity and significantly reduces the healthcare cost, and it is therefore essential for the sustainability of the Spanish healthcare system.

Genetic Variation in CYP2D6 and SLC22A1 Affects Amlodipine Pharmacokinetics and Safety

Amlodipine is an antihypertensive drug with unknown pharmacogenetic biomarkers. This research is a candidate gene study that looked for associations between amlodipine pharmacokinetics and safety and pharmacogenes. Pharmacokinetic and safety data were taken from 160 volunteers from eight bioequivalence trials. In the exploratory step, 70 volunteers were genotyped for 44 polymorphisms in different pharmacogenes. CYP2D6 poor metabolizers (PMs) showed higher half-life (t1/2) (univariate p-value (puv) = 0.039, multivariate p-value (pmv) = 0.013, β = -5.31, R2 = 0.176) compared to ultrarapid (UMs), normal (NMs) and intermediate metabolizers (IMs). SLC22A1 rs34059508 G/A genotype was associated with higher dose/weight-corrected area under the curve (AUC72/DW) (puv = 0.025; pmv = 0.026, β = 578.90, R2 = 0.060) compared to the G/G genotype. In the confirmatory step, the cohort was increased to 160 volunteers, who were genotyped for CYP2D6, SLC22A1 and CYP3A4. In addition to the previous associations, CYP2D6 UMs showed a lower AUC72/DW (puv = 0.046, pmv = 0.049, β = -68.80, R2 = 0.073) compared to NMs, IMs and PMs and the SLC22A1 rs34059508 G/A genotype was associated with thoracic pain (puv = 0.038) and dizziness (puv = 0.038, pmv = 0.014, log OR = 10.975). To our knowledge, this is the first work to report a strong relationship between amlodipine and CYP2D6 and SLC22A1. Further research is needed to gather more evidence before its application in clinical practice.

CYP2C8*3 and *4 define CYP2C8 phenotype: An approach with the substrate cinitapride

Cinitapride is a gastrointestinal prokinetic drug, prescribed for the treatment of functional dyspepsia, and as an adjuvant therapy for gastroesophageal reflux disease. In this study, we aimed to explore the impact of relevant variants in CYP3A4 and CYP2C8 and other pharmacogenes, along with demographic characteristics, on cinitapride pharmacokinetics and safety; and to evaluate the impact of CYP2C8 alleles on the enzyme's function. Twenty-five healthy volunteers participating in a bioequivalence clinical trial consented to participate in the study. Participants were genotyped for 56 variants in 19 genes, including cytochrome P450 (CYP) enzymes (e.g., CYP2C8 or CYP3A4) or transporters (e.g., SLC or ABC), among others. CYP2C8*3 carriers showed a reduction in AUC of 42% and C_max of 35% compared to *1/*1 subjects (p = 0.003 and p = 0.011, respectively). *4 allele carriers showed a 45% increase in AUC and 63% in C_max compared to *1/*1 subjects, although these differences did not reach statistical significance. CYP2C8*3 and *4 alleles may be used to infer the following pharmacogenetic phenotypes: ultrarapid (UM) (*3/*3), rapid (RM) (*1/*3), normal (NM) (*1/*1), intermediate (IM) (*1/*4), and poor (PM) metabolizers (*4/*4). In this study, we properly characterized RMs, NMs, and IMs; however, additional studies are required to properly characterize UMs and PMs. These findings should be relevant with respect to cinitapride, but also to numerous CYP2C8 substrates such as imatinib, loperamide, montelukast, ibuprofen, paclitaxel, pioglitazone, repaglinide, or rosiglitazone.

Polymorphism of Drug Transporters, Rather Than Metabolizing Enzymes, Conditions the Pharmacokinetics of Rasagiline

Rasagiline is a selective and irreversible inhibitor of monoamine oxidase type B with neuroprotective effect, indicated for the management of Parkinson’s disease. The aim of this work was to evaluate the impact of seven CYP1A2 alleles and of 120 additional variants located in other CYP enzymes (e.g., CYP2C19), UGT enzymes (e.g., UGT1A1) or other enzymes (e.g., NAT2), and transporters (e.g., SLCO1B1) on the pharmacokinetic variability and safety of rasagiline. A total of 118 healthy volunteers enrolled in four bioequivalence clinical trials consented to participate in this pharmacogenetic study. CYP1A2 alleles were not associated with the pharmacokinetic variability of rasagiline. Patients with ABCB1 rs1045642 G/A+A/A genotypes presented higher area under the curve adjusted by dose per weight (AUC_0-∞/DW) than those with the G/G genotype (p = 0.012) and lower volume of distribution (V_d/F) and clearance (Cl/F) (p = 0.001 and p = 0.012, respectively). Subjects with the ABCC2 rs2273697 A/A genotype presented lower t_max (i.e., the time to reach the maximum concentration, C_max) compared to those with G/G+G/A genotypes (p = 0.001). Volunteers with the SLC22A1 *1/*5 genotype exhibited lower C_max/DW and higher t_max (p = 0.003 and p = 0.018, respectively) than subjects with the *1/*1 diplotype. Only one adverse drug reaction was reported: headache. Our results suggest the genetic polymorphism of drug transporters, rather than metabolizing enzymes, conditions the pharmacokinetics of rasagiline.

Genotype-Guided Prescription of Azathioprine Reduces the Incidence of Adverse Drug Reactions in TPMT Intermediate Metabolizers to a Similar Incidence as Normal Metabolizers

Introduction

Thiopurine drugs are purine nucleoside analogues used for treatment of different immune-related conditions. To date, different studies highlighted the importance of thiopurine methyltransferase (TPMT) genotyping in patients who initiate treatment with thiopurines to make an adequate dose adjustment. We aimed to investigate the influence of TPMT phenotype, concomitant treatments, and demographic characteristics on the incidence of adverse reactions (ADRs) in patients who start treatment with azathioprine (AZA).

Methods

This was an observational and retrospective study. The study population comprised 109 patients who started treatment with AZA following routine TPMT genotyping before June 2019 and who were routinely followed up at Hospital Universitario de La Princesa. The incidence of ADRs and treatment duration were evaluated according to TPMT phenotype.

Results

Forty-five men and 64 women were recruited, with a mean age of 67.6 ± 18.5. The medical specialty with the most requests was dermatology (45.9%) and the most frequent disease for which genotyping was requested was bullous pemphigoid (27.5%). All patients were normal metabolizers (NM), except for eight intermediate metabolizers (IM) (7.3%); no poor metabolizers (PM) were found. The initial azathioprine dose was subtherapeutic in both groups (103.2 ± 45.4 mg in NMs and 75 ± 32.3 mg in IMs), increasing during the first months of treatment, especially in NMs (120.3 ± 41.3 vs. 78.6 ± 30.4 mg in IMs, p = 0.011). Most patients (73.4%) received corticosteroids to keep the disease under control; and for 41.2% of NMs, physicians were able to reduce the dose at 6 months post treatment. No IMs completed 6 months of treatment. Hepatotoxicity, gastric intolerance, and blood disorders were the most common ADRs. The incidence of ADRs in the sample was 28.4% (n = 31) with a similar trend between IMs (37.5%) and NMs (27.8%). Patients undergoing concomitant treatment with allopurinol were associated with a higher incidence of ADRs (n = 4, 100% vs. n = 105, 20%; p = 0.002).

Conclusion

TPMT genotyping before AZA prescription reduces ADR incidence in IMs to a similar level as NMs in the Spanish population. However, it is important to note no IMs completed 6 months of treatment, suggesting that there may be some differences in drug tolerability according to phenotype. In addition, most NMs are treated with subtherapeutic doses, are poorly followed up, and thus suffer avoidable ADRs. Finally, concomitant therapies that inhibit the xanthine oxidase enzyme (XDH), such as allopurinol, predispose to ADRs. Therefore, pharmacogenetic testing should be integrated as an additional clinical tool, in such a way that each patient receives personalized, precision treatment, where all factors influencing drug response are considered.

Patients’ Perceptions of Pharmacogenetic Testing and Access to Their Results: State of the Art in Spain and Systematic Review

The process of clinical pharmacogenetics implementation depends on patients’ and general population’s perceptions. To date, no study has been published addressing Spanish patients’ opinions on pharmacogenetic testing, the availability of the results, and the need for signing informed consent. In this work, we contacted 146 patients that had been previously genotyped at our laboratory and 46 healthy volunteers that had participated in a bioequivalence clinical trial at the Clinical Pharmacology Department of Hospital Universitario de La Princesa and consented to pharmacogenetic testing for research purposes. From the latter, 108 and 34, respectively, responded to the questionnaire (i.e., a response rate of 74%); Participants were scheduled for a face-to-face, telephone, or videoconference interview and were asked a total of 27 questions in Spanish. Great or almost complete acceptance of pharmacogenetic testing was observed (99.3%), age and university education level being the main predictors of acceptance rates and understanding (multivariate analysis, p = 0.004, R2 = 0.17, age being inversely proportional to acceptance rates and understanding and university level being related to higher acceptance rates and understanding compared to other education levels). Mixed perceptions were observed on the requirement of written informed consent (55.6% in favor); therefore, it seems recommendable to continue requesting it for the upcoming years until more perceptions are collected. The majority of participants (95.8%) preferred storing pharmacogenetic results in medical records rather than in electronic sources (55.6%) and highly agreed with the possibility of carrying their results on a portable card (91.5%). Patients agreed to broad genetic testing, including biomarkers unrelated to their disease (93.7%) or with little clinically relevant evidence (94.4%). Patients apparently rely on clinician’s or pharmacogeneticist’s interpretation and seem, therefore, open to the generation of ethically challenging information. Finally, although most patients (68.3%) agreed with universal population testing, some were reluctant, probably due to the related costs and sustainability of the Spanish Health System. This was especially evident in the group of patients who were older and with a likely higher proportion of pensioners.

Clinical Relevance of Novel Polymorphisms in the Dihydropyrimidine Dehydrogenase (DPYD) Gene in Patients with Severe Fluoropyrimidine Toxicity: A Spanish Case-Control Study

Among cancer patients treated with fluoropyrimidines, 10-40% develop severe toxicity. Polymorphism of the dihydropyrimidine dehydrogenase (DPYD) gene may reduce DPD function, the main enzyme responsible for the metabolism of fluoropyrimidines. This leads to drug accumulation and to an increased risk of toxicity. Routine genotyping of this gene, which usually includes DPYD *HapB3, *2A, *13 and c.2846A > T (D949V) variants, helps predict approximately 20-30% of toxicity cases. For DPD intermediate (IM) or poor (PM) metabolizers, a dose adjustment or drug switch is warranted to avoid toxicity, respectively. Societies such as the Spanish Society of Pharmacogenetics and Pharmacogenomics (SEFF), the Dutch Pharmacogenetics Working Group (DPWG) or the Clinical Pharmacogenetics Implementation Consortium (CPIC) and regulatory agencies (e.g., the Spanish Medicines Agency, AEMPS) already recommend DPYD routine genotyping. However, the predictive capacity of genotyping is currently still limited. This can be explained by the presence of unknown polymorphisms affecting the function of the enzyme. In this case-control work, 11 cases of severe fluoropyrimidine toxicity in patients who did not carry any of the four variants mentioned above were matched with 22 controls, who did not develop toxicity and did not carry any variant. The DPYD exome was sequenced (Sanger) in search of potentially pathogenic mutations. DPYD rs367619008 (c.187 A > G, p.Lys63Glu), rs200643089 (c.2324 T > G, p.Leu775Trp) and rs76387818 (c.1084G > A, p.Val362Ile) increased the percentage of explained toxicities to 38-48%. Moreover, there was an intronic variant considered potentially pathogenic: rs944174134 (c.322-63G > A). Further studies are needed to confirm its clinical relevance. The remaining variants were considered non-pathogenic.