Rare versus common variants in pharmacogenetics: SLCO1B1 variation and methotrexate disposition

Identification of Genetic Variants Associated with Pravastatin and Pitavastatin Pharmacokinetics

A clinical trial was carried out to investigate the pharmacogenetics of single-dose pravastatin and pitavastatin pharmacokinetics in 173 and 164 healthy participants. Additionally, 96 participants were included from previous pharmacogenetic studies with pravastatin. In a genome-wide meta-analysis of pravastatin including all 269 participants, SLCO1B1 c.521T>C (rs4149056) was associated with increased AUC0-∞ (P = 9.8 × 10-12). Similarly, SLCO1B1 c.521T>C was genome-wide significantly associated with increased AUC0-∞ of pitavastatin (P = 9.7 × 10-15). Candidate gene analyses suggested that participants with increased function SLCO1B1 variants had decreased pravastatin exposure. Furthermore, decreased function CYP2C9 variants may increase pitavastatin and pitavastatin lactone exposure. Compared to participants with normal function SLCO1B1 genotype, the AUC0-∞ of pravastatin was 140% (90% confidence interval: 86-210%; P = 4.7 × 10-8) and 37% (20-56%; P = 1.1 × 10-4) greater in participants with poor and decreased function SLCO1B1 genotype, respectively, while participants with highly increased function SLCO1B1 genotype had a 60% (39-75%; P = 6.0 × 10-4) lower AUC0-∞. The AUC0-∞ of pitavastatin was 153% (100-222%; P = 1.6 × 10-9) and 35% (8-69%; P = 0.027) greater in participants with poor and decreased function SLCO1B1 genotype, respectively, than in those with normal function SLCO1B1 genotype. Participants with intermediate metabolizer CYP2C9 genotype had 18% (3-34%; P = 0.046) greater AUC0-∞ of pitavastatin than those with normal metabolizer CYP2C9 genotype. This study demonstrates the important role of SLCO1B1 in pravastatin and pitavastatin pharmacokinetics and suggests that CYP2C9 variants also affect the pharmacokinetics of pitavastatin.

Translating pharmacogenomic sequencing data into drug response predictions-How to interpret variants of unknown significance

The rapid development of sequencing technologies during the past 20 years has provided a variety of methods and tools to interrogate human genomic variations at the population level. Pharmacogenes are well known to be highly polymorphic and a plethora of pharmacogenomic variants has been identified in population sequencing data. However, so far only a small number of these variants have been functionally characterized regarding their impact on drug efficacy and toxicity and the significance of the vast majority remains unknown. It is therefore of high importance to develop tools and frameworks to accurately infer the effects of pharmacogenomic variants and, eventually, aggregate the effect of individual variations into personalized drug response predictions. To address this challenge, we here first describe the technological advances, including sequencing methods and accompanying bioinformatic processing pipelines that have enabled reliable variant identification. Subsequently, we highlight advances in computational algorithms for pharmacogenomic variant interpretation and discuss the added value of emerging strategies, such as machine learning and the integrative use of omics techniques that have the potential to further contribute to the refinement of personalized pharmacological response predictions. Lastly, we provide an overview of experimental and clinical approaches to validate in silico predictions. We conclude that the iterative feedback between computational predictions and experimental validations is likely to rapidly improve the accuracy of pharmacogenomic prediction models, which might soon allow for an incorporation of the entire pharmacogenetic profile into personalized response predictions.

Implementation of pre-emptive testing of a pharmacogenomic panel in clinical practice: Where do we stand?

Adverse drug reactions (ADRs) account for a large proportion of hospitalizations among adults and are more common in multimorbid patients, worsening clinical outcomes and burdening healthcare resources. Over the past decade, pharmacogenomics has been developed as a practical tool for optimizing treatment outcomes by mitigating the risk of ADRs. Some single-gene reactive tests are already used in clinical practice, including the DPYD test for fluoropyrimidines, which demonstrates how integrating pharmacogenomic data into routine care can improve patient safety in a cost-effective manner. The evolution from reactive single-gene testing to comprehensive pre-emptive genotyping panels holds great potential for refining drug prescribing practices. Several implementation projects have been conducted to test the feasibility of applying different genetic panels in clinical practice. Recently, the results of a large prospective randomized trial in Europe (the PREPARE study by Ubiquitous Pharmacogenomics consortium) have provided the first evidence that prospective application of a pre-emptive pharmacogenomic test panel in clinical practice, in seven European healthcare systems, is feasible and yielded a 30% reduction in the risk of developing clinically relevant toxicities. Nevertheless, some important questions remain unanswered and will hopefully be addressed by future dedicated studies. These issues include the cost-effectiveness of applying a pre-emptive genotyping panel, the role of multiple co-medications, the transferability of currently tested pharmacogenetic guidelines among patients of non-European origin and the impact of rare pharmacogenetic variants that are not detected by currently used genotyping approaches.

Evaluation of methotrexate Pharmacogenomic variation to predict acute neurotoxicity in children with acute lymphoblastic leukemia

BACKGROUND

Methotrexate is an important component of curative therapy in childhood acute lymphoblastic leukemia (ALL), but the role of genetic variation influencing methotrexate clearance and transport in toxicity susceptibility in children with ALL is not well established. Therefore, we evaluated the association between suspected methotrexate pharmacogenomic variants and methotrexate-related neurotoxicity.

METHODS

This study included children (aged 2-20 years) diagnosed with ALL (2005-2019) at six treatment centers in the southwest United States. Clinical information was abstracted from medical records. Suspected neurotoxic events occurring within 21 days of intravenous and/or intrathecal methotrexate delivered between the end of induction and start of maintenance therapy were independently reviewed by at least two pediatric oncologists. Germline DNA was genotyped and 97 methotrexate pharmacogenomic variants of interest with at least grade 3 evidence were identified using the Pharmacogenomics Knowledge Base. Associations between variants and neurotoxicity were assessed by logistic regression. Data were randomly split (80/20) and random forest was constructed to estimate the ability of the variants to correctly classify neurotoxicity.

RESULTS

Of the 763 patients included in the study, 8.2% (n = 63) developed methotrexate-associated neurotoxicity. In logistic models, none of the 97 available pharmacogenomic variants reached adjusted statistical significance. However, two variants, rs17222723 (odds ratio [OR] = 2.83 [ref. = T allele], 95% confidence interval [CI]: 1.20-6.15) in ABCC2 and rs1045642 (OR = 0.66 [ref. = minor A allele], 95% CI: 0.44-0.98) in ABCB1, were nominally associated (p-value < 0.05) with neurotoxicity susceptibility. The addition of pharmacogenomic variants did not improve the predictive performance of random forest model (AUC = 0.73) compared to clinical information alone (AUC = 0.74).

CONCLUSION

Overall, our results suggest that associations between neurotoxicity susceptibility and methotrexate pharmacogenomic variants are generally modest and these variants do not significantly improve neurotoxicity risk stratification among children with ALL.

Integration of genomics, clinical characteristics and baseline biological profiles to predict the risk of liver injury induced by high-dose methotrexate

Background

High-dose methotrexate (HD-MTX) is commonly employed in the treatment of malignant tumors in children and young adults due to its distinctive therapeutic efficacy. Nonetheless, the systemic exposure to MTX often results in liver injury (drug induced liver injury, DILI), thereby imposing limitations on the sustained administration of HD-MTX. Additionally, individual variations including genetic underpinnings attributable to disparities in therapeutic effects and clinical toxicity remain to be elucidated.

Methods

A total of 374 patients receiving initial HD-MTX treatment were selected for this study, which aimed to establish a predictive model using binary logistic regression and a visual nomogram for DILI risk assessment. Demographic and clinical characteristics were collected at baseline and post-HD-MTX to explore their correlations with the occurrence of DILI. Additionally, genotyping of 25 single nucleotide polymorphisms from drug transporters and enzymes in the folic acid cycle was performed.

Result

G allele mutation in ABCB1 rs1128503, *1b/*1b and *1b/*15 haplotypic mutation in SLCO1B1, female gender, and MTX dosage were identified as independent factors for moderate/severe DILI. Patients with GA or AA genotype in ABCB1 rs1128503 showed significant higher 24h MTX concentration than GG, and those with *1b/*1b haplotype group in SLCO1B1 exhibited lower dose adjusted concentration (C/D) than *1a/*1a group. Besides, patient administrated with HD-MTX were more prevalent to have higher C/D levels when using intravenous plus triple intrathecal injection route than those who were using intravenous injection alone. The composite predictive model (ROC curve: AUC = 0.805), comprising above four factors and 24h MTX concentration, exhibited high accuracy.

Conclusion

Female gender, recessive mutation in ABCB1 rs1128503, and a range of MTX concentration may be risk factors for increased susceptibility to DILI. Conversely, the *1b/*1b and *1b/*15 mutations in SLCO1B1 may have a protective effect against DILI. The proposed predictive model facilitates early individual risk assessment, enabling the implementation of proactive prevention strategies.

In silico and biological analyses of missense variants of the human biliary efflux transporter ABCC2: effects of novel rare missense variants

BACKGROUND AND PURPOSE

The ATP-dependent biliary efflux transporter ABCC2, also known as multidrug resistance protein 2 (MRP2), is essential for the cellular disposition and detoxification of various xenobiotics including drugs as well as endogenous metabolites. Common functionally relevant ABCC2 genetic variants significantly alter drug responses and contribute to side effects. The aim of this study was to determine functional consequences of rare variants identified in subjects with European ancestry using in silico tools and in vitro analyses.

EXPERIMENTAL APPROACH

Targeted next-generation sequencing of the ABCC2 gene was used to identify novel variants in European subjects (n = 143). Twenty-six in silico tools were used to predict functional consequences. For biological validation, transport assays were carried out with membrane vesicles prepared from cell lines overexpressing the newly identified ABCC2 variants and estradiol β-glucuronide and carboxydichlorofluorescein as the substrates.

KEY RESULTS

Three novel rare ABCC2 missense variants were identified (W227R, K402T, V489F). Twenty-five in silico tools predicted W227R as damaging and one as potentially damaging. Prediction of functional consequences was not possible for K402T and V489F and for the common linked variants V1188E/C1515Y. Characterisation in vitro showed increased function of W227R, V489F and V1188E/C1515Y for both substrates, whereas K402T function was only increased for carboxydichlorofluorescein.

CONCLUSION AND IMPLICATIONS

In silico tools were unable to accurately predict the substrate-dependent increase in function of ABCC2 missense variants. In vitro biological studies are required to accurately determine functional activity to avoid misleading consequences for drug therapy.

Clinicopathological characteristics of low-dose methotrexate-induced epidermal dysmaturation: A study of 22 patients

BACKGROUND AND OBJECTIVE

Erosions of the skin and mucous membranes with epidermal dysmaturation are a known side effect of cytostatic chemotherapy regimens and can also be observed during low-dose methotrexate (MTX) therapy. The study aimed to delineate the clinical and histopathological alterations.

PATIENTS AND METHODS

A database search of the archive for dermatopathology was conducted, identifying 22 patients who developed epidermal dysmaturation on low-dose MTX. Clinical and laboratory changes, along with an array of histologic parameters were analyzed and statistically evaluated using SPSS.

RESULTS

Patients were predominantly female with a mean age of 69.1 years. The main indications were psoriasis vulgaris and rheumatoid arthritis. Clinically, patients mostly presented erosive plaques at the injection site, on mucosal surfaces, and disseminated lesions. Most patients showed normal laboratory values. Histopathologically, key findings included enlarged keratinocytes with pale cytoplasm and enlarged nuclei with prominent nucleoli, along with the degeneration of the basal layer. Consistent observations in the dermal compartment included infiltration of neutrophilic granulocytes, lymphocytes, and histiocytes.

CONCLUSIONS

This study proposes clinicopathological criteria for the diagnosis of MTX-associated skin toxicity, aiming to increase awareness among clinicians and pathologists for early diagnosis. Early recognition can prevent potentially life-threatening progression.

Impact of liver diseases and pharmacological interactions on the transportome involved in hepatic drug disposition

The liver plays a pivotal role in drug disposition owing to the expression of transporters accounting for the uptake at the sinusoidal membrane and the efflux across the basolateral and canalicular membranes of hepatocytes of many different compounds. Moreover, intracellular mechanisms of phases I and II biotransformation generate, in general, inactive compounds that are more polar and easier to eliminate into bile or refluxed back toward the blood for their elimination by the kidneys, which becomes crucial when the biliary route is hampered. The set of transporters expressed at a given time, i.e., the so-called transportome, is encoded by genes belonging to two gene superfamilies named Solute Carriers (SLC) and ATP-Binding Cassette (ABC), which account mainly, but not exclusively, for the uptake and efflux of endogenous substances and xenobiotics, which include many different drugs. Besides the existence of genetic variants, which determines a marked interindividual heterogeneity regarding liver drug disposition among patients, prevalent diseases, such as cirrhosis, non-alcoholic steatohepatitis, primary sclerosing cholangitis, primary biliary cirrhosis, viral hepatitis, hepatocellular carcinoma, cholangiocarcinoma, and several cholestatic liver diseases, can alter the transportome and hence affect the pharmacokinetics of drugs used to treat these patients. Moreover, hepatic drug transporters are involved in many drug-drug interactions (DDI) that challenge the safety of using a combination of agents handled by these proteins. Updated information on these questions has been organized in this article by superfamilies and families of members of the transportome involved in hepatic drug disposition.

Systematic Evaluation of Tyrosine Kinase Inhibitors as OATP1B1 Substrates Using a Competitive Counterflow Screen

Although the primary elimination pathway for most tyrosine kinase inhibitors (TKI) involves CYP3A4-mediated metabolism, the mechanism by which these agents are brought into hepatocytes remains unclear. In this study, we optimized and validated a competitive counterflow (CCF) assay to examine TKIs as substrates of the hepatic uptake transporter OATP1B1. The CCF method was based on the stimulated efflux of radiolabeled estradiol-17β-glucuronide under steady-state conditions in HEK293 cells engineered to overexpress OATP1B1. Of the 62 approved TKIs examined, 13 agents were identified as putative substrates of OATP1B1, and pazopanib was selected as a representative hit for further validation studies. The transport of pazopanib by OATP1B1 was confirmed by decreased activity of its target VEGFR2 in OATP1B1-overexpressing cells, but not cells lacking OATP1B1, consistent with molecular docking analyses indicating an overlapping binding orientation on OATP1B1 with the known substrate estrone-3-sulfate. In addition, the liver-to-plasma ratio of pazopanib in vivo was decreased in mice with a deficiency of the orthologous transporters, and this was accompanied by diminished pazopanib-induced hepatotoxicity, as determined by changes in the levels of liver transaminases. Our study supports the utility of CCF assays to assess substrate affinity for OATP1B1 within a large set of agents in the class of TKIs and sheds light on the mechanism by which these agents are taken up into hepatocytes in advance of metabolism.

SIGNIFICANCE

Despite the established exposure-pharmacodynamic relationships for many TKIs, the mechanisms underlying the agents' unpredictable pharmacokinetic profiles remain poorly understood. We report here that the disposition of many TKIs depends on hepatic transport by OATP1B1, a process that has toxicologic ramifications for agents that are associated with hepatotoxicity.

Transporter-mediated drug-drug interactions: regulatory guidelines, in vitro and in vivo methodologies and translation, special populations, and the blood-brain barrier

This review, part of a special issue on drug-drug interactions (DDIs) spearheaded by the International Society for the Study of Xenobiotics (ISSX) New Investigators, explores the critical role of drug transporters in absorption, disposition, and clearance in the context of DDIs. Over the past two decades, significant advances have been made in understanding the clinical relevance of these transporters. Current knowledge on key uptake and efflux transporters that affect drug disposition and development is summarized. Regulatory guidelines from the FDA, EMA, and PMDA that inform the evaluation of potential transporter-mediated DDIs are discussed in detail. Methodologies for preclinical and clinical testing to assess potential DDIs are reviewed, with an emphasis on the utility of physiologically based pharmacokinetic (PBPK) modeling. This includes the application of relative abundance and expression factors to predict human pharmacokinetics (PK) using preclinical data, integrating the latest regulatory guidelines. Considerations for assessing transporter-mediated DDIs in special populations, including pediatric, hepatic, and renal impairment groups, are provided. Additionally, the impact of transporters at the blood-brain barrier (BBB) on the disposition of CNS-related drugs is explored. Enhancing the understanding of drug transporters and their role in drug disposition and toxicity can improve efficacy and reduce adverse effects. Continued research is essential to bridge remaining gaps in knowledge, particularly in comparison with cytochrome P450 (CYP) enzymes.

Clinical Implementation of Rare and Novel DPYD Variants for Personalizing Fluoropyrimidine Treatment: Challenges and Opportunities

Fluoropyrimidines (FLs) [5-Fluorouracil, Capecitabine] are used in the treatment of several solid tumors. Dihydropyrimidine dehydrogenase (DPD) is the rate-limiting enzyme for FL detoxification, and its deficiency could lead to severe, life-threatening or fatal toxicity after FL administration. Testing with a pharmacogenetic panel of four deleterious variants in the dihydropyrimidine dehydrogenase gene (DPYD) (DPYD*2A, DPYD*13, c.2846A > T, c.1129-5923C > G) prior to FL treatment, is recommended by scientific consortia (e.g., CPIC, DPWG) and drug regulatory agencies (e.g., EMA). However, this panel identifies < 20% of patients at risk of severe FL-related toxicity. Cumulative recent evidence highlights the potential clinical value of rare (minor allele frequency < 1%) and novel DPYD genetic variants for identifying an additional fraction of DPD-deficient patients at increased risk of severe FL-related toxicity. In this review, we aimed to comprehensively describe the available evidence regarding the potential clinical predictive role of novel and rare DPYD variants as toxicity markers in FL-treated patients, and to discuss the challenges and opportunities in tailoring FL treatment based upon clinical application of such markers. Although we must overcome existing barriers to the clinical implementation, the available data support that comprehensive assessment of the DPYD sequence, including rare and novel genetic variants, may significantly enhance the pre-emptive identification of at-risk patients, compared to the current targeted approach.

Clinical effects of tacrolimus blood concentrations early after allogeneic hematopoietic stem cell transplantation

BACKGROUND AIMS

Tacrolimus (TAC) plus short-term methotrexate (stMTX) is used for graft-versus-host disease (GVHD) prophylaxis after allogeneic hematopoietic stem cell transplantation (allo-HSCT). TAC blood concentrations are frequently adjusted to enhance the graft-versus-leukemia/lymphoma effect or attenuate severe GVHD. Limited information is available on the clinical impact of these adjustments and the optimal time to perform them in order to achieve good clinical outcomes.

METHODS

We retrospectively analyzed 211 patients who underwent allo-HSCT at our institutes.

RESULTS

Higher TAC concentrations in week 3 correlated with a significantly higher cumulative incidence of relapse (CIR) (P = 0.03) and lower nonrelapse mortality (P = 0.04). The clinical impact of high TAC concentrations in week 3 on CIR was detected in the refined disease risk index: low/intermediate (P = 0.04) and high (P < 0.01), and conditioning regimens other than cyclophosphamide/total body irradiation and busulfan/cyclophosphamide (P = 0.07). Higher TAC concentrations in week 1 correlated with a lower grade 2-4 acute GVHD rate (P = 0.01). Higher TAC concentrations in weeks 2 and 3 correlated with slightly lower (P = 0.05) and significantly lower (P = 0.02) grade 3-4 acute GVHD rates, respectively. Higher TAC concentrations in weeks 1 and 3 were beneficial for severe acute GVHD in patients with a human leukocyte antigen-matched donor (P = 0.03 and P < 0.01, respectively), not treated with anti-thymocyte globulin (P = 0.02 and P = 0.02, respectively), and receiving three stMTX doses (P = 0.03 and P = 0.02, respectively).

CONCLUSIONS

The clinical impact of TAC concentrations varied according to patient characteristics, including disease malignancy, conditioning regimens, donor sources, and GVHD prophylaxis. These results suggest that TAC management needs to be based on patient profiles.

SLCO1B1 Exome Sequencing and Statin Treatment Response in 64,000 UK Biobank Patients

The solute carrier organic anion transporter family member 1B1 (SLCO1B1) encodes the organic anion-transporting polypeptide 1B1 (OATP1B1 protein) that transports statins to liver cells. Common genetic variants in SLCO1B1, such as *5, cause altered systemic exposure to statins and therefore affect statin outcomes, with potential pharmacogenetic applications; yet, evidence is inconclusive. We studied common and rare SLCO1B1 variants in up to 64,000 patients from UK Biobank prescribed simvastatin or atorvastatin, combining whole-exome sequencing data with up to 25-year routine clinical records. We studied 51 predicted gain/loss-of-function variants affecting OATP1B1. Both SLCO1B1*5 alone and the SLCO1B1*15 haplotype increased LDL during treatment (beta*5 = 0.08 mmol/L, p = 6 × 10-8; beta*15 = 0.03 mmol/L, p = 3 × 10-4), as did the likelihood of discontinuing statin prescriptions (hazard ratio*5 = 1.12, p = 0.04; HR*15 = 1.05, p = 0.04). SLCO1B1*15 and SLCO1B1*20 increased the risk of General Practice (GP)-diagnosed muscle symptoms (HR*15 = 1.22, p = 0.003; HR*20 = 1.25, p = 0.01). We estimated that genotype-guided prescribing could potentially prevent 18% and 10% of GP-diagnosed muscle symptoms experienced by statin patients, with *15 and *20, respectively. The remaining common variants were not individually significant. Rare variants in SLCO1B1 increased LDL in statin users by up to 1.05 mmol/L, but replication is needed. We conclude that genotype-guided treatment could reduce GP-diagnosed muscle symptoms in statin patients; incorporating further SLCO1B1 variants into clinical prediction scores could improve LDL control and decrease adverse events, including discontinuation.

Risk factors associated with high-dose methotrexate induced toxicities

INTRODUCTION

High-dose methotrexate (HDMTX) therapy poses challenges in various neoplasms due to individualized pharmacokinetics and associated adverse effects. Our purpose is to identify early risk factors associated with HDMTX-induced toxicities, paving the way for personalized treatment.

AREAS COVERED

A systematic review of PubMed and Cochrane databases was conducted for articles from inception to July 2023. Eligible studies included reviews, clinical trials, and real-world analyses. Irrelevant studies were excluded, and manual searches and citation reviews were performed. Factors such as MTX exposure, drug interactions, demographics, serum albumin, urine pH, serum calcium, and genetic polymorphisms affecting MTX transport (e.g. SLCO1B1), intracellular folate metabolism (MTHFR), cell development (ARID5B), metabolic pathways (UGT1A1, PNPLA3), as well as epigenetics were identified.

EXPERT OPINION

This comprehensive review aids researchers and clinicians in early identification of HDMTX toxicity risk factors. By understanding the multifaceted risk factors associated with hematologic malignancies, personalized treatment approaches can be tailored to optimize therapeutic outcomes.

Large-scale next generation sequencing based analysis of SLCO1B1 pharmacogenetics variants in the Saudi population

BACKGROUND

SLCO1B1 plays an important role in mediating hepatic clearance of many different drugs including statins, angiotensin-converting enzyme inhibitors, chemotherapeutic agents and antibiotics. Several variants in SLCO1B1 have been shown to have a clinically significant impact, in relation to efficacy of these medications. This study provides a comprehensive overview of SLCO1B1 variation in Saudi individuals, one of the largest Arab populations in the Middle East.

METHODS

The dataset of 11,889 (9,961 exomes and 1,928 pharmacogenetic gene panel) Saudi nationals, was used to determine the presence and frequencies of SLCO1B1 variants, as described by the Clinical Pharmacogenetic Implementation Consortium (CPIC).

RESULTS

We identified 141 previously described SNPs, of which rs2306283 (50%) and rs4149056 (28%), were the most common. In addition, we observed six alleles [*15 (24.7%) followed by *20 (8.04%), *14 (5.86%), *5 (3.84%), *31 (0.21%) and *9 (0.03%)] predicted to be clinically actionable. Allele diplotype to phenotype conversion revealed 41 OATP1B1 diplotypes. We estimated the burden of rare, and novel predicted deleterious variants, resulting from 17 such alterations.

CONCLUSIONS

The data we present, from one of the largest Arab cohorts studied to date, provides the most comprehensive overview of SLCO1B1 variants, and the subsequent OATP1B1 activity of this ethnic group, which thus far remains relatively underrepresented in available international genomic databases. We believe that the presented data provides a basis for further clinical investigations and the application of personalized statin drug therapy guidance in Arabs.

Methotrexate level discrepancy post-glucarpidase: A pediatric case series and review of literature

Methotrexate is a common component of pediatric oncology treatment and delayed clearance increases risk of significant toxicities. Glucarpidase is indicated for patients with toxic plasma methotrexate concentrations with renal toxicity. Laboratory interference with immunoassay measurement post-glucarpidase administration is well established, with current product labeling indicating this persists for 48 h. However, recent experience in pediatric patients supports this discrepancy persists beyond 48 h. Three cases experienced delayed methotrexate clearance and received glucarpidase with subsequent measurement of methotrexate levels by liquid chromatography tandem mass spectrometry (LC-MS/MS) and/or immunoassay. Within this case series, discrepancies between LC-MS/MS and immunoassay levels persisted significantly longer than 48 h.

Non-targeted metabolomics for the identification of plasma metabolites associated with organic anion transporting polypeptide 1B1 function

Our aim was to evaluate biomarkers for organic anion transporting polypeptide 1B1 (OATP1B1) function using a hypothesis-free metabolomics approach. We analyzed fasting plasma samples from 356 healthy volunteers using non-targeted metabolite profiling by liquid chromatography high-resolution mass spectrometry. Based on SLCO1B1 genotypes, we stratified the volunteers to poor, decreased, normal, increased, and highly increased OATP1B1 function groups. Linear regression analysis, and random forest (RF) and gradient boosted decision tree (GBDT) regressors were used to investigate associations of plasma metabolite features with OATP1B1 function. Of the 9152 molecular features found, 39 associated with OATP1B1 function either in the linear regression analysis (p < 10-5) or the RF or GBDT regressors (Gini impurity decrease > 0.01). Linear regression analysis showed the strongest associations with two features identified as glycodeoxycholate 3-O-glucuronide (GDCA-3G; p = 1.2 × 10-20 for negative and p = 1.7 × 10-19 for positive electrospray ionization) and one identified as glycochenodeoxycholate 3-O-glucuronide (GCDCA-3G; p = 2.7 × 10-16). In both the RF and GBDT models, the GCDCA-3G feature showed the strongest association with OATP1B1 function, with Gini impurity decreases of 0.40 and 0.17. In RF, this was followed by one GDCA-3G feature, an unidentified feature with a molecular weight of 809.3521, and the second GDCA-3G feature. In GBDT, the second and third strongest associations were observed with the GDCA-3G features. Of the other associated features, we identified with confidence two representing lysophosphatidylethanolamine 22:5. In addition, one feature was putatively identified as pregnanolone sulfate and one as pregnenolone sulfate. These results confirm GCDCA-3G and GDCA-3G as robust OATP1B1 biomarkers in human plasma.

Next-generation sequencing in pharmacogenomics - fit for clinical decision support?

INTRODUCTION

The technological advances of sequencing methods during the past 20 years have fuelled the generation of large amounts of sequencing data that comprise common variations, as well as millions of rare and personal variants that would not be identified by conventional genotyping. While comprehensive sequencing is technically feasible, its clinical utility for guiding personalized treatment decisions remains controversial.

AREAS COVERED

We discuss the opportunities and challenges of comprehensive sequencing compared to targeted genotyping for pharmacogenomic applications. Current pharmacogenomic sequencing panels are heterogeneous and clinical actionability of the included genes is not a major focus. We provide a current overview and critical discussion of how current studies utilize sequencing data either retrospectively from biobanks, databases or repurposed diagnostic sequencing, or prospectively using pharmacogenomic sequencing.

EXPERT OPINION

While sequencing-based pharmacogenomics has provided important insights into genetic variations underlying the safety and efficacy of a multitude pharmacological treatments, important hurdles for the clinical implementation of pharmacogenomic sequencing remain. We identify gaps in the interpretation of pharmacogenetic variants, technical challenges pertaining to complex loci and variant phasing, as well as unclear cost-effectiveness and incomplete reimbursement. It is critical to address these challenges in order to realize the promising prospects of pharmacogenomic sequencing.

GENETIC PREDICTORS OF TOXIC EFFECTS OF METHOTREXATE IN CANCER PATIENTS

Today, methotrexate (MTX) is used in combination with other medicines to treat a wide range of malignancies. Despite its proven high efficacy, MTX often causes serious side effects, which may result in the need to reduce the dose of MTX or discontinue the drug altogether. This, in turn, can provoke the development of MTX resistance and cancer progression. Predicting the risk of MTX-induced toxicity is currently difficult due to the variability of pharmacokinetics and pharmacodynamics in different patients, so the scientific literature is intensively searching for potential biomarkers. Based on the data available in the current literature, we analyzed the relationship between variants in the genes encoding the key components of MTX intracellular metabolism and the MTX-induced side effects and drug response. According to the results of our work, the most studied variants are those of the SLC19A1 gene, which encodes the reduced folate carrier protein 1, and the MTHFR gene, which encodes the enzyme methylenetetrahydrofolate reductase. Studies of the effect of methylation of the promoter regions of genes on the therapeutic effect of MTX are also very promising. In conclusion, the study of molecular genetic markers of MTX toxicity is extremely relevant and necessary because it can help to avoid the effect of multidrug resistance and improve the quality of life and survival of patients.

Genetic risk factors for drug-induced long QT syndrome: findings from a large real-world case-control study

Aim: Drug-induced long QT syndrome (diLQTS), an adverse effect of many drugs, can lead to sudden cardiac death. Candidate genetic variants in cardiac ion channels have been associated with diLQTS, but several limitations of previous studies hamper clinical utility. Materials & methods: Thus, the purpose of this study was to assess the associations of KCNE1-D85N, KCNE2-I57T and SCN5A-G615E with diLQTS in a large observational case-control study (6,083 self-reported white patients treated with 27 different high-risk QT-prolonging medications; 12.0% with diLQTS). Results: KCNE1-D85N significantly associated with diLQTS (adjusted odds ratio: 2.24 [95% CI: 1.35-3.58]; p = 0.001). Given low minor allele frequencies, the study had insufficient power to analyze KCNE2-I57T and SCN5A-G615E. Conclusion: KCNE1-D85N is a risk factor for diLQTS that should be considered in future clinical practice guidelines.

Mining local exome and HLA data to characterize pharmacogenetic variants in Saudi Arabia

Pharmacogenomics (PGx) is a promising field of precision medicine where efficacy of drugs is maximized while side effects are minimized for individual patients. Knowledge of the frequency of PGx-relevant variants (pharmacovariants) in the local population is a pre-requisite to informed policy making. Unfortunately, such knowledge is largely lacking from the Middle East. Here, we describe the use of a large clinical exome database (n = 13,473) and HLA haplotypes (n = 64,737) from Saudi Arabia, one of the largest countries in the Middle East, along with previously published data from the local population to ascertain allele frequencies of known pharmacovariants. In addition, we queried another exome database (n = 816) of well-phenotyped research subjects from Saudi Arabia to discover novel candidate variants in known PGx genes (pharmacogenes). Although our results show that only 26% (63/242) of class 1A/1B PharmGKB variants were identified, we estimate that 99.57% of the local population have at least one such variant. This translates to a minimum estimated impact of 9% of medications dispensed by our medical center annually. We also highlight the contribution of rare variants where 71% of the pharmacogenes devoid of common pharmacovariants had at least one potentially deleterious rare variant. Thus, we show that approaches that go beyond the use of commercial PGx kits that have been optimized for other populations should be implemented to ensure universal and equitable access of all members of the local population to personalized prescription practices.

Pharmacogenetic aspects of efficacy and safety of methotrexate treatment in pediatric acute lymphoblastic leukemia

OBJECTIVES

To evaluate the role of ABCB1 (C3435T rs1045642, rs1128503, rs2032582, rs4148738), SLCO1B1 T521C rs4149056 genetic polymorphisms in the development of major types of methotrexate (MTX) toxicities and the occurrence of a terminal event (death, relapse) in pediatric АLL.

METHODS

The study included 124 patients diagnosed with pediatric ALL. All patients treated according to the protocols of the German BFM group (2002/2009) with high-dose (1,000, 2,000 and 5,000 mg/m2) methotrexate. MTX-related toxicities, including hematologic, hepatic and renal, were evaluated according to the common terminology criteria for adverse events version 5.0 (CTCAE v.5.0). Real-time PCR method was used to investigate polymorphisms of ABCB1 and SLCO1B1 genes. The study material was peripheral blood.

RESULTS

A competitive analysis demonstrated significant relationships between MTX ADRs. The results of the study support the existence of relationships between some ADRs and MTX kinetics. An associative analysis showed association with the development of AEs to methotrexate indicating their clinical significance from different genetic polymorphisms protein-transporters. The available results confirm the associations of the studied genes with the increased risk of high doses MTX toxic ADRs and terminal events.

CONCLUSIONS

Complementing the existing criteria for pediatric ALL risk groups with pharmacogenetic indicators will allow further individualization of therapy.

The burden of rare variants in DPYS gene is a novel predictor of the risk of developing severe fluoropyrimidine-related toxicity

BACKGROUND

Despite a growing number of publications highlighting the potential impact on the therapy outcome, rare genetic variants (minor allele frequency < 1%) in genes associated to drug adsorption, distribution, metabolism, and elimination are poorly studied. Previously, rare germline DPYD missense variants were shown to identify a subset of fluoropyrimidine-treated patients at high risk for severe toxicity. Here, we investigate the impact of rare genetic variants in a panel of 54 other fluoropyrimidine-related genes on the risk of severe toxicity.

METHODS

The coding sequence and untranslated regions of 54 genes related to fluoropyrimidine pharmacokinetics/pharmacodynamics were analyzed by next-generation sequencing in 120 patients developing grade 3-5 toxicity (NCI-CTC vs3.0) and 104 matched controls. Sequence Kernel Association Test (SKAT) analysis was used to select genes with a burden of genetic variants significantly associated with risk of severe toxicity. The statistical association of common and rare genetic variants in selected genes was further investigated. The functional impact of genetic variants was assessed using two different in silico prediction tools (Predict2SNP; ADME Prediction Framework).

RESULTS

SKAT analysis highlighted DPYS and PPARD as genes with a genetic mutational burden significantly associated with risk of severe fluoropyrimidine-related toxicity (Bonferroni adjusted P = 0.024 and P = 0.039, respectively). Looking more closely at allele frequency, the burden of rare DPYS variants was significantly higher in patients with toxicity compared with controls (P = 0.047, Mann-Whitney test). Carrying at least one rare DPYS variant was associated with an approximately fourfold higher risk of severe cumulative (OR = 4.08, P = 0.030) and acute (OR = 4.21, P = 0.082) toxicity. The burden of PPARD rare genetic variants was not significantly related to toxicity. Some common variants with predictive value in DPYS and PPARD were also identified: DPYS rs143004875-T and PPARD rs2016520-T variants predicted an increased risk of severe cumulative (P = 0.002 and P = 0.001, respectively) and acute (P = 0.005 and P = 0.0001, respectively) toxicity.

CONCLUSION

This work demonstrated that the rare mutational burden of DPYS, a gene strictly cooperating with DPYD in the catabolic pathway of fluoropyrimidines, is a promising pharmacogenetic marker for precision dosing of fluoropyrimidines. Additionally, some common genetic polymorphisms in DPYS and PPARD were identified as promising predictive markers that warrant further investigation.

Uptake Transporters at the Blood-Brain Barrier and Their Role in Brain Drug Disposition

Uptake drug transporters play a significant role in the pharmacokinetic of drugs within the brain, facilitating their entry into the central nervous system (CNS). Understanding brain drug disposition is always challenging, especially with respect to preclinical to clinical translation. These transporters are members of the solute carrier (SLC) superfamily, which includes organic anion transporter polypeptides (OATPs), organic anion transporters (OATs), organic cation transporters (OCTs), and amino acid transporters. In this systematic review, we provide an overview of the current knowledge of uptake drug transporters in the brain and their contribution to drug disposition. Here, we also assemble currently available proteomics-based expression levels of uptake transporters in the human brain and their application in translational drug development. Proteomics data suggest that in association with efflux transporters, uptake drug transporters present at the BBB play a significant role in brain drug disposition. It is noteworthy that a significant level of species differences in uptake drug transporters activity exists, and this may contribute toward a disconnect in inter-species scaling. Taken together, uptake drug transporters at the BBB could play a significant role in pharmacokinetics (PK) and pharmacodynamics (PD). Continuous research is crucial for advancing our understanding of active uptake across the BBB.

Association between SLCO1B1 genetic polymorphisms and bleeding risk in patients treated with edoxaban

Since SLCO1B1 encodes the uptake transporter OATP1B1, which can influence the pharmacokinetic and pharmacodynamic profiles of edoxaban, polymorphisms in SLCO1B1 may affect the edoxaban response. This study aimed to investigate the association between SLCO1B1 gene polymorphisms and the bleeding risk in patients receiving edoxaban. We genotyped 10 single-nucleotide polymorphisms (SNPs) from the SLCO1B1 gene in patients receiving edoxaban. We also analyzed rs3842 of ABCB1 as a confounder. The odds ratio (OR) and adjusted OR (AOR) were calculated from univariate and multivariable analysis, respectively. The area under the receiver operating characteristic curve (AUROC) was constructed for the discrimination of the model. A total of 159 patients receiving edoxaban were analyzed. Overdose and rs4149056 showed significant association with bleeding complications by around 11- and 5.5-fold, respectively. Additionally, patients with the rs4149057 variant allele (C) had a 3.9-fold increased bleeding risk compared with wild-type homozygote carriers (TT), whereas rs2306283 variant homozygote (GG) carriers had a 0.27-fold reduced bleeding risk compared with wild-type allele (A) carriers. Patients with the variant-type homozygote (CC) of ABCB1 rs3842 had a higher bleeding risk than T allele carriers (AOR = 5.3 and 5.9). The final models for multivariable analyses were acceptable based on the AUROC values (> 0.70). These findings may help predict bleeding risk in patients taking edoxaban and help personalize treatment.

Germline genetic variants and pediatric rhabdomyosarcoma outcomes: a report from the Children's Oncology Group

BACKGROUND

Relative to other pediatric cancers, survival for rhabdomyosarcoma (RMS) has not improved in recent decades, suggesting the need to enhance risk stratification. Therefore, we conducted a genome-wide association study for event-free survival (EFS) and overall survival (OS) to identify genetic variants associated with outcomes in individuals with RMS.

METHODS

The study included 920 individuals with newly diagnosed RMS who were enrolled in Children's Oncology Group protocols. To assess the association of each single nucleotide polymorphism (SNP) with EFS and OS, we estimated hazard ratios (HRs) and 95% confidence intervals (CIs) using multivariable Cox proportional hazards models, adjusted for clinical covariates. All statistical tests were two sided. We also performed stratified analyses by histological subtype (alveolar and embryonal RMS) and carried out sensitivity analyses of statistically significant SNPs by PAX3/7-FOXO1 fusion status and genetic ancestry group.

RESULTS

We identified that rs17321084 was associated with worse EFS (HR = 2.01, 95% CI = 1.59 to 2.53, P = 5.39 × 10-9) and rs10094840 was associated with worse OS (HR = 1.84, 95% CI = 1.48 to 2.27, P = 2.13 × 10-8). Using publicly available data, we found that rs17321084 lies in a binding region for transcription factors GATA2 and GATA3, and rs10094840 is associated with SPAG1 and RNF19A expression. We also identified that CTNNA3 rs2135732 (HR = 3.75, 95% CI = 2.34 to 5.99, P = 3.54 × 10-8) and MED31 rs74504320 (HR = 3.21, 95% CI = 2.12 to 4.86, P = 3.60 × 10-8) were associated with worse OS among individuals with alveolar RMS.

CONCLUSIONS

We demonstrated that common germline variants are associated with EFS and OS among individuals with RMS. Additional replication and investigation of these SNP effects may further support their consideration in risk stratification protocols.

External evaluation and systematic review of population pharmacokinetic models for high-dose methotrexate in cancer patients

Several population pharmacokinetic (PPK) models have been established to optimize the therapeutic regimen and reduce the toxicity of high-dose methotrexate (HDMTX) in patients with cancer. However, their predictive performance when extrapolated to different clinical centers was unknown. In this study, we aimed to externally evaluate the predictive ability of HDMTX PPK models and determine the potential influencing factors. We searched the literature and determined the predictive performance of the selected models using methotrexate concentrations in 721 samples from 60 patients in the First Affiliated Hospital of the Navy Medical University. Prediction-based diagnostics and simulation-based normalized prediction distribution errors (NPDE) were used to evaluate the predictive performance of the models. The influence of prior information was also assessed using Bayesian forecasting, and the potential factors affecting model predictability were investigated. Thirty models extracted from published PPK studies were assessed. Prediction-based diagnostics showed that the number of compartments potentially influenced model transferability, and simulation-based NPDE indicated model misspecification. Bayesian forecasting significantly improved the predictive performance of the models. Various factors, including bioassays, covariates, and population diagnosis, influence model extrapolation. The published models were unsatisfactory for all prediction-based diagnostics, except for the 24 h methotrexate concentration monitoring and simulation-based diagnostics, making them inappropriate for direct extrapolation. Moreover, Bayesian forecasting combined therapeutic drug monitoring could improve the predictive performance of the models.

PharmVar GeneFocus: SLCO1B1

The Pharmacogene Variation Consortium (PharmVar) is now providing star (*) allele nomenclature for the highly polymorphic human SLCO1B1 gene encoding the organic anion transporting polypeptide 1B1 (OATP1B1) drug transporter. Genetic variation within the SLCO1B1 gene locus impacts drug transport, which can lead to altered pharmacokinetic profiles of several commonly prescribed drugs. Variable OATP1B1 function is of particular importance regarding hepatic uptake of statins and the risk of statin-associated musculoskeletal symptoms. To introduce this important drug transporter gene into the PharmVar database and serve as a unified reference of haplotype variation moving forward, an international group of gene experts has performed an extensive review of all published SLCO1B1 star alleles. Previously published star alleles were self-assigned by authors and only loosely followed the star nomenclature system that was first developed for cytochrome P450 genes. This nomenclature system has been standardized by PharmVar and is now applied to other important pharmacogenes such as SLCO1B1. In addition, data from the 1000 Genomes Project and investigator-submitted data were utilized to confirm existing haplotypes, fill knowledge gaps, and/or define novel star alleles. The PharmVar-developed SLCO1B1 nomenclature has been incorporated by the Clinical Pharmacogenetics Implementation Consortium (CPIC) 2022 guideline on statin-associated musculoskeletal symptoms.

Functional Characterization of Six SLCO1B1 (OATP1B1) Variants Observed in Finnish Individuals with a Psychotic Disorder

Variants in the SLCO1B1 (solute carrier organic anion transporter family member 1B1) gene encoding the OATP1B1 (organic anion transporting polypeptide 1B1) protein are associated with altered transporter function that can predispose patients to adverse drug effects with statin treatment. We explored the effect of six rare SLCO1B1 single nucleotide variants (SNVs) occurring in Finnish individuals with a psychotic disorder on expression and functionality of the OATP1B1 protein. The SUPER-Finland study has performed exome sequencing on 9381 individuals with at least one psychotic episode during their lifetime. SLCO1B1 SNVs were annotated with PHRED-scaled combined annotation-dependent (CADD) scores and the Ensembl variant effect predictor. In vitro functionality studies were conducted for the SNVs with a PHRED-scaled CADD score of >10 and predicted to be missense. To estimate possible changes in transport activity caused by the variants, transport of 2',7'-dichlorofluorescein (DCF) in OATP1B1-expressing HEK293 cells was measured. According to the findings, additional tests with rosuvastatin and estrone sulfate were conducted. The amount of OATP1B1 in crude membrane fractions was quantified using a liquid chromatography tandem mass spectrometry-based quantitative targeted absolute proteomics analysis. Six rare missense variants of SLCO1B1 were identified in the study population, located in transmembrane helix 3: c.317T>C (p.106I>T), intracellular loop 2: c.629G>T (p.210G>V), c.633A>G (p.211I>M), c.639T>A (p.213N>L), transmembrane helix 6: 820A>G (p.274I>V), and the C-terminal end: 2005A>C (p.669N>H). Of these variants, SLCO1B1 c.629G>T (p.210G>V) resulted in the loss of in vitro function, abolishing the uptake of DCF, estrone sulfate, and rosuvastatin and reducing the membrane protein expression to 31% of reference OATP1B1. Of the six rare missense variants, SLCO1B1 c.629G>T (p.210G>V) causes a loss of function of OATP1B1 transport in vitro and severely decreases membrane protein abundance. Carriers of SLCO1B1 c.629G>T might be susceptible to altered pharmacokinetics of OATP1B1 substrate drugs and might have increased likelihood of adverse drug effects such as statin-associated musculoskeletal symptoms.

Novel and replicated clinical and genetic risk factors for toxicity from high-dose methotrexate in pediatric acute lymphoblastic leukemia

STUDY OBJECTIVE

Methotrexate (MTX) is a key component of treatment for high-risk pediatric acute lymphoblastic leukemia (ALL) but may cause acute kidney injury and prolonged hospitalization due to delayed clearance. The purpose of this study is to identify clinical and genetic factors that may predict which children are at risk for creatinine increase and prolonged MTX clearance.

DESIGN

We conducted a single-center, retrospective cohort study of pediatric patients with ALL who received 4000-5000 mg/m2 of MTX. Measurements We performed germline genotyping to determine genetic ancestry and allele status for 49 single nucleotide polymorphisms (SNPs) identified from the literature as related to MTX disposition. Bayesian hierarchical ordinal regression models for creatinine increase and for prolonged MTX clearance were developed.

MAIN RESULTS

Hispanic ethnicity, body mass index (BMI) < 3%, BMI between 85%-95%, and Native American genetic ancestry were found to be associated with an increased risk for creatinine elevation. Older age, Black race, and use of the intensive monitoring protocol were associated with a decreased risk for creatinine elevation. Older age, B- compared to T-ALL, and the minor alleles of rs2838958/SLC19A1 and rs7317112/ABCC4 were associated with an increased risk for delayed clearance. Black race, MTX dose reduction, and the minor allele of rs2306283/SLCO1B1 were found to be associated with a decreased risk for delayed clearance.

CONCLUSIONS

These predictors of MTX toxicities may allow for more precise individualized toxicity risk prediction.

Targeted next-generation sequencing of genes involved in Warfarin Pharmacodynamics and pharmacokinetics pathways using the Saudi Warfarin Pharmacogenetic study (SWAP)

BACKGROUND

Warfarin is an oral anticoagulant commonly used for treatment and prophylaxis against thromboembolic events. Warfarins's narrow therapeutic index window is one of the main challenges in clinical practice; thus, it requires frequent monitoring and dose adjustment to maintain patients' therapeutic range. Warfarin dose variation and response are attributed to several inter-and intra-individuals factors, including genetic variants in enzymes involved in warfarin pharmacokinetics (PK) and pharmacodynamics (PD) pathways. Thus, we aim to utilize the next-generation sequencing (NGS) approach to identify rare and common genetic variants that might be associated with warfarin responsiveness.

METHOD AND RESULTS

A predesigned NGS panel that included 16 genes involved in Warfarin PK/PD pathways was used to sequence 786 patients from the Saudi Warfarin Pharmacogenetic Cohort (SWAP). Identified variants were annotated using several annotation tools to identify the pathogenicity and allele frequencies of these variants. We conducted variants-level association tests with warfarin dose. We identified 710 variants within the sequenced genes; 19% were novel variants, with the vast majority being scarce variants. The genetic association tests showed that VKORC1 (rs9923231, and rs61742245), CYP2C9 (rs98332238, rs9332172, rs1057910, rs9332230, rs1799853, rs1057911, and rs9332119), CYP2C19 (rs28399511, and rs3758581), and CYP2C8 (rs11572080 and rs10509681) were significantly associated with warfarin weekly dose. Our model included genetics, and non-genetic factors explained 40.1% of warfarin dose variation.

CONCLUSION

The study identifies novel variants associated with warfarin dose in the Saudi population. These variants are more likely to be population-specific variants, suggesting that population-specific studies should be conducted before adopting a universal warfarin genotype-guided dosing algorithm.

Insights from a pharmacometric analysis of HDMTX in adults with cancer: Clinically relevant covariates for application in precision dosing

AIMS

High-dose methotrexate (HDMTX) is an essential part of the treatment of several adult and paediatric malignancies. Despite meticulous supportive care during HDMTX administration, severe toxicities, including acute kidney injury (AKI), may occur contributing to patient morbidity. Population pharmacokinetics provide a powerful tool to predict time to clear HDMTX and adjust subsequent doses. We sought to develop and validate pharmacokinetic models for HDMTX in adults with diverse malignancies and to relate systemic exposure with the occurrence of severe toxicity.

METHODS

Anonymized, de-identified data were provided from 101 US oncology practices that participate in the Guardian Research Network, a non-profit clinical research consortium. Modelled variables included clinical, laboratory, demographic and pharmacological data. Population pharmacokinetic analysis was performed by means of nonlinear mixed effects modelling using MonolixSuite.

RESULTS

A total of 693 HDMTX courses from 243 adults were analysed, of which 62 courses (8.8%) were associated with stage 2/3 acute kidney injury (43 stage 2, 19 stage 3). A three-compartment model adequately fitted the data. Time-dependent serum creatinine, baseline serum albumin and allometrically scaled bodyweight were clinically significant covariates related to methotrexate clearance. External evaluation confirmed a satisfactory predictive performance of the model in adults receiving HDMTX. Dose-normalized methotrexate concentration at 24 and 48 hours correlated with AKI incidence.

CONCLUSION

We developed a population pharmacometric model that considers weight, albumin and time-dependent creatinine that can be used to guide supportive care in adult patients with delayed HDMTX elimination.

Case report: Hepatotoxicity and nephrotoxicity induced by methotrexate in a paediatric patient, what is the role of precision medicine in 2023?

Methotrexate is an immunosuppressant and chemotherapeutic agent used in the treatment of a range of autoimmune disorders and cancers. Its main serious adverse effects, bone marrow suppression and gastrointestinal complications, arise from its antimetabolite effect. Nevertheless, hepatotoxicity and nephrotoxicity are two widely described adverse effects of methotrexate. Its hepatotoxicity has been studied mainly in the low-dose, chronic setting, where patients are at risk of fibrosis/cirrhosis. Studies of acute hepatoxicity of high dose methotrexate, such as during chemotherapy, are scarce. We present the case of a 14-year-old patient who received high-dose methotrexate and subsequently developed acute fulminant liver failure and acute kidney injury. Genotyping of MTHFR (Methylene tetrahydrofolate reductase gene), ABCB1 (codes for P-glycoprotein, intestinal transport and biliary excretion), ABCG2 (codes for BCRP, intestinal transporter and renal excretion) and SLCO1B1 (codes for OATP1B1, hepatic transporter) identified variants in all the genes analysed that predicted a reduced rate of methotrexate elimination and thus may have contributed to the clinical situation of the patient. Precision medicine involving pharmacogenomic testing could potentially avoid such adverse drug effects.

The genetic landscape of major drug metabolizing cytochrome P450 genes-an updated analysis of population-scale sequencing data

Genes encoding cytochrome P450 enzymes (CYPs) are extremely polymorphic and multiple CYP variants constitute clinically relevant biomarkers for the guidance of drug selection and dosing. We previously reported the distribution of the most relevant CYP alleles using population-scale sequencing data. Here, we update these findings by making use of the increasing wealth of data, incorporating whole exome and whole genome sequencing data from 141,614 unrelated individuals across 12 human populations. We furthermore extend our previous studies by systematically considering also uncharacterized rare alleles and reveal that they contribute between 1.5% and 17.5% to the overall genetically encoded functional variability. By using established guidelines, we aggregate and translate the available sequencing data into population-specific patterns of metabolizer phenotypes. Combined, the presented data refine the worldwide landscape of ethnogeographic variability in CYP genes and aspire to provide a relevant resource for the optimization of population-specific genotyping strategies and precision public health.

Association between high-dose methotrexate-induced toxicity and polymorphisms within methotrexate pathway genes in acute lymphoblastic leukemia

Methotrexate (MTX) is a folic acid antagonist, the mechanism of action is to inhibit DNA synthesis, repair and cell proliferation by decreasing the activities of several folate-dependent enzymes. It is widely used as a chemotherapy drug for children and adults with malignant tumors. High-dose methotrexate (HD-MTX) is an effective treatment for extramedullary infiltration and systemic consolidation in children with acute lymphoblastic leukemia (ALL). However, significant toxicity results in most patients treated with HD-MTX, which limits its use. HD-MTX-induced toxicity is heterogeneous, and this heterogeneity may be related to gene polymorphisms in related enzymes of the MTX intracellular metabolic pathway. To gain a deeper understanding of the differences in toxicity induced by HD-MTX in individuals, the present review examines the correlation between HD-MTX-induced toxicity and the gene polymorphisms of related enzymes in the MTX metabolic pathway in ALL. In this review, we conclude that only the association of SLCO1B1 and ARID5B gene polymorphisms with plasma levels of MTX and MTX-related toxicity is clearly described. These results suggest that SLCO1B1 and ARID5B gene polymorphisms should be evaluated before HD-MTX treatment. In addition, considering factors such as age and race, the other exact predictor of MTX induced toxicity in ALL needs to be further determined.

Effects of SLCO1B1 Genetic Variant on Metabolite Profile in Participants on Simvastatin Treatment

Organic-anion-transporting polypeptide 1B1 (OATP1B1), encoded by the solute carrier organic anion transporter family member 1B1 gene (SLCO1B1), is highly expressed in the liver and transports several endogenous metabolites into the liver, including statins. Previous studies have not investigated the association of SLCO1B1 rs4149056 variant with the risk of type 2 diabetes (T2D) or determined the metabolite signature of the C allele of SLCO1B1 rs4149056 (SLCO1B1 rs4149056-C allele) in a large randomly selected population. SLCO1B1 rs4149056-C inhibits OATP1B1 transporter and is associated with increased levels of blood simvastatin concentrations. Our study is to first to show that SLCO1B1 rs4149056 variant is not significantly associated with the risk of T2D, suggesting that simvastatin has a direct effect on the risk of T2D. Additionally, we investigated the effects of SLCO1B1 rs4149056-C on plasma metabolite concentrations in 1373 participants on simvastatin treatment and in 1368 age- and body-mass index (BMI)-matched participants without any statin treatment. We found 31 novel metabolites significantly associated with SLCO1B1 rs4149056-C in the participants on simvastatin treatment and in the participants without statin treatment. Simvastatin decreased concentrations of dicarboxylic acids, such as docosadioate and dodecanedioate, that may increase beta- and peroxisomal oxidation and increased the turnover of cholesterol into bile acids, resulting in a decrease in steroidogenesis due to limited availability of cholesterol for steroid synthesis. Our findings suggest that simvastatin exerts its effects on the lowering of low-density lipoprotein (LDL) cholesterol concentrations through several distinct pathways in the carriers of SLCO1B1 rs4149056-C, including dicarboxylic acids, bile acids, steroids, and glycerophospholipids.

Effects of genetic polymorphisms on methotrexate levels and toxicity in Chinese patients with acute lymphoblastic leukemia

Methotrexate (MTX) has an antitumor effect when used for the treatment of acute lymphoblastic leukemia (ALL). This study aims at evaluating the associations between 14 polymorphisms of six genes involved in MTX metabolism with serum MTX concentration and toxicity accompanying high-dose MTX. Polymorphisms in 183 Chinese patients with ALL were analyzed using TaqMan single nucleotide polymorphism genotyping assay. The serum MTX concentration was determined using homogeneous enzyme immunoassay. MTX-related toxicities were also evaluated. Renal toxicity was significantly associated with higher serum MTX concentrations at 24, 48, and 72 hours, and MTX elimination delay (P = 0.001, P < 0.001, P < 0.001, and P < 0.001, respectively), whereas SLCO1B1 rs4149056 was associated with serum MTX concentrations at 48 and 72 hours, and MTX elimination delay in candidate polymorphisms (P = 0.014, P = 0.019, and P = 0.007, respectively). SLC19A1 rs2838958 and rs3788200 were associated with serum MTX concentrations at 24 hours (P = 0.016, P = 0.043, respectively). MTRR rs1801394 was associated with serum MTX concentrations at 72 hours (P = 0.045). Neutropenia was related to SLC19A1 rs4149056 (odds ratio [OR]: 3.172, 95% confidence interval [CI]: 1.310-7.681, P = 0.011). Hepatotoxicity was associated with ABCC2 rs2273697 (OR: 3.494, 95% CI: 1.236-9.873, P = 0.018) and MTRR rs1801394 (OR: 0.231, 95% CI: 0.084-0.632, P = 0.004). Polymorphisms of SLCO1B1, SLC19A1, ABCC2, and MTRR genes help predict higher risk of increased MTX levels or MTX-related toxicities in adult ALL patients.

Organic Anion Transporter Gene Variants Associated With Plasma Exposure and Long-Term Response to Atrasentan in Patients With Diabetic Kidney Disease

Plasma exposure of the endothelin receptor antagonist atrasentan varies between individuals and is associated with nephroprotective effects and the risk of heart failure. We examined the influence of genetic polymorphisms on atrasentan plasma exposure and pharmacodynamic effects. We performed a substudy of the Study of Diabetic Nephropathy With Atrasentan (SONAR) trial which enrolled adults with type 2 diabetes and chronic kidney disease (estimated glomerular filtration rate: 25-75 mL/min/1.73 m² , and a urine albumin-to-creatinine ratio of 300-5,000 mg/g). Single nucleotide polymorphisms (SNPs) were determined for prespecified membrane transporters, metabolizing enzymes, and the endothelin-1 peptide. The associations among genotype, atrasentan plasma exposure, and the effect of atrasentan on the prespecified kidney and heart failure hospitalization (HHF) outcomes was assessed with Cox proportional hazards regression models. Of 3,668 patients randomized, 2,329 (63.5%) consented to genotype analysis. Two SNPs in the SLCO1B1 gene (rs4149056 and rs2306283), encoding the hepatic organic anion transporter 1B1 (OATP1B1), showed the strongest association with atrasentan plasma exposure. Based on their SLCO1B1 genotype, patients were classified into normal (atrasentan area under the plasma-concentration time curve from zero to infinity (AUC_0-inf ) 41.3 ng·h/mL) or slow (atrasentan AUC_0-inf 49.7 ng·h/mL, P < 0.001) OATP1B1 transporter phenotypes. Among patients with a normal OATP1B1 phenotype, the hazard ratio (HR) with atrasentan for the primary kidney and HHF outcomes were 0.61 (95% confidence interval (CI): 0.45-0.81) and 1.35 (95% CI: 0.84-2.13), respectively. In the slow transporter phenotype, HRs for kidney and HHF outcomes were 1.95 (95% CI: 0.95-4.03, P-interaction normal phenotype = 0.004), and 4.18 (95% CI: 1.37-12.7, P-interaction normal phenotype = 0.060), respectively. OATP1B1 gene polymorphisms are associated with significant between-patient variability in atrasentan plasma exposure and long-term efficacy and safety.

Pharmacometabolomics-Based Translational Biomarkers: How to Navigate the Data Ocean

Metabolome is the end point of the genome-environment interplay, and enables an important holistic overview of individual adaptability and host responses to environmental, ecological, as well as endogenous changes such as disease. Pharmacometabolomics is the application of metabolome knowledge to decipher the mechanisms of interindividual and intraindividual variations in drug efficacy and safety. Pharmacometabolomics also contributes to prediction of drug treatment outcomes on the basis of baseline (predose) and postdose metabotypes through mathematical modeling. Thus, pharmacometabolomics is a strong asset for a diverse community of stakeholders interested in theory and practice of evidence-based and precision/personalized medicine: academic researchers, public health scholars, health professionals, pharmaceutical, diagnostics, and biotechnology industries, among others. In this expert review, we discuss pharmacometabolomics in four contexts: (1) an interdisciplinary omics tool and field to map the mechanisms and scale of interindividual variability in drug effects, (2) discovery and development of translational biomarkers, (3) advance digital biomarkers, and (4) empower drug repurposing, a field that is increasingly proving useful in the current era of Covid-19. As the applications of pharmacometabolomics are growing rapidly in the current postgenome era, next-generation proteomics and metabolomics follow the example of next-generation sequencing analyses. Pharmacometabolomics can also empower data reliability and reproducibility through multiomics integration strategies, which use each data layer to correct, connect with, and inform each other. Finally, we underscore here that contextual data remain crucial for precision medicine and drug development that stand the test of time and clinical relevance.

The role of OATP1B1 and OATP1B3 transporter polymorphisms in drug disposition and response to anticancer drugs: a review of the recent literature

INTRODUCTION

Members of the solute carrier family of organic anion transporting polypeptides are responsible for the cellular uptake of a broad range of endogenous compounds and xenobiotics in multiple tissues. In particular, the polymorphic transporters OATP1B1 and OATP1B3 are highly expressed in the liver and have been identified as critical regulators of hepatic eliminaton. As these transporters are also expressed in cancer cells, the function alteration of these proteins have important consequences for an individual's susceptibility to certain drug-induced side effects, drug-drug interactions, and treatment efficacy.

AREAS COVERED

In this mini-review, we provide an update of this rapidly emerging field, with specific emphasis on the direct contribution of genetic variants in OATP1B1 and OATP1B3 to the transport of anticancer drugs, the role of these carriers in regulation of their disposition and toxicity profiles, and recent advances in attempts to integrate information on transport function in patients to derive individualized treatment strategies.

EXPERT OPINION

Based on currently available data, it appears imperative that different aspects of disease, physiology, and drugs of relevance should be evaluated along with an individual's genetic signature, and that tools such as biomarker levels can be implemented to achieve the most reliable prediction of clinically relevant pharmacodynamic endpoints.

Functional in vitro characterization of SLCO1B1 variants and simulation of the clinical pharmacokinetic impact of impaired OATP1B1 function

Organic Anion Transporting Polypeptide 1B1 is important to the hepatic elimination and distribution of many drugs. If OATP1B1 function is decreased, it can increase plasma exposure of e.g. several statins leading to increased risk of muscle toxicity. First, we examined the impact of three naturally occurring rare variants and the frequent SLCO1B1 c.388A>G variant on in vitro transport activity with cellular uptake assay using two substrates: 2', 7'-dichlorofluorescein (DCF) and rosuvastatin. Secondly, LC-MS/MS based quantitative targeted absolute proteomics measured the OATP1B1 protein abundance in crude membrane fractions of HEK293 cells over-expressing these SNVs. Additionally, we simulated the effect of impaired OATP1B1 function on rosuvastatin pharmacokinetics to estimate the need for genotype-guided dosing. R57Q impaired DCF and rosuvastatin transport significantly yet did not change protein expression considerably, while N130D and N151S did not alter activity but increased protein expression. R253Q did not change protein expression but reduced DCF uptake and increased rosuvastatin K_m. Based on pharmacokinetics simulations, doses of 30 mg (with 50% OATP1B1 function) and 20 mg (with 0% OATP1B1 function) result in plasma exposure similar to 40 mg dose (with 100% OATP1B1 function). Therefore dose reductions might be considered to avoid increased plasma exposure caused by function-impairing OATP1B1 genetic variants, such as R57Q.

Maintenance therapy for acute lymphoblastic leukemia: basic science and clinical translations

Maintenance therapy (MT) with oral methotrexate (MTX) and 6-mercaptopurine (6-MP) is essential for the cure of acute lymphoblastic leukemia (ALL). MTX and 6-MP interfere with nucleotide synthesis and salvage pathways. The primary cytotoxic mechanism involves the incorporation of thioguanine nucleotides (TGNs) into DNA (as DNA-TG), which may be enhanced by the inhibition of de novo purine synthesis by other MTX/6-MP metabolites. Co-medication during MT is common. Although Pneumocystis jirovecii prophylaxis appears safe, the benefit of glucocorticosteroid/vincristine pulses in improving survival and of allopurinol to moderate 6-MP pharmacokinetics remains uncertain. Numerous genetic polymorphisms influence the pharmacology, efficacy, and toxicity (mainly myelosuppression and hepatotoxicity) of MTX and thiopurines. Thiopurine S-methyltransferase (encoded by TPMT) decreases TGNs but increases methylated 6-MP metabolites (MeMPs); similarly, nudix hydrolase 15 (encoded by NUDT15) also decreases TGNs available for DNA incorporation. Loss-of-function variants in both genes are currently used to guide MT, but do not fully explain the inter-patient variability in thiopurine toxicity. Because of the large inter-individual variations in MTX/6-MP bioavailability and metabolism, dose adjustments are traditionally guided by the degree of myelosuppression, but this does not accurately reflect treatment intensity. DNA-TG is a common downstream metabolite of MTX/6-MP combination chemotherapy, and a higher level of DNA-TG has been associated with a lower relapse hazard, leading to the development of the Thiopurine Enhanced ALL Maintenance (TEAM) strategy-the addition of low-dose (2.5-12.5 mg/m²/day) 6-thioguanine to the 6-MP/MTX backbone-that is currently being tested in a randomized ALLTogether1 trial (EudraCT: 2018-001795-38). Mutations in the thiopurine and MTX metabolism pathways, and in the mismatch repair genes have been identified in early ALL relapses, providing valuable insights to assist the development of strategies to detect imminent relapse, to facilitate relapse salvage therapy, and even to bring about changes in frontline ALL therapy to mitigate this relapse risk.

Genomewide Association Study of Simvastatin Pharmacokinetics

We investigated genetic determinants of single-dose simvastatin pharmacokinetics in a prospective study of 170 subjects and a retrospective cohort of 59 healthy volunteers. In a microarray-based genomewide association study with the prospective data, the SLCO1B1 c.521T>C (p.Val174Ala, rs4149056) single nucleotide variation showed the strongest, genomewide significant association with the area under the plasma simvastatin acid concentration-time curve (AUC; P = 6.0 × 10^-10 ). Meta-analysis with the retrospective cohort strengthened the association (P = 1.6 × 10^-17 ). In a stepwise linear regression candidate gene analysis among all 229 participants, SLCO1B1 c.521T>C (P = 1.9 × 10^-13 ) and CYP3A4 c.664T>C (p.Ser222Pro, rs55785340, CYP3A4*2, P = 0.023) were associated with increased simvastatin acid AUC. Moreover, the SLCO1B1 c.463C>A (p.Pro155Thr, rs11045819, P = 7.2 × 10^-6 ) and c.1929A>C (p.Leu643Phe, rs34671512, P = 5.3 × 10^-4 ) variants associated with decreased simvastatin acid AUC. Based on these results and the literature, we classified the volunteers into genotype-predicted OATP1B1 and CYP3A4 phenotype groups. Compared with the normal OATP1B1 function group, simvastatin acid AUC was 273% larger in the poor (90% confidence interval (CI), 137%, 488%; P = 3.1 × 10^-6 ), 40% larger in the decreased (90% CI, 8%, 83%; P = 0.036), and 67% smaller in the highly increased function group (90% CI, 46%, 80%; P = 2.4 × 10^-4 ). Intermediate CYP3A4 metabolizers (i.e., heterozygous carriers of either CYP3A4*2 or CYP3A4*22 (rs35599367)), had 87% (90% CI, 39%, 152%, P = 6.4 × 10^-4 ) larger simvastatin acid AUC than normal metabolizers. These data suggest that in addition to no function SLCO1B1 variants, increased function SLCO1B1 variants and reduced function CYP3A4 variants may affect the pharmacokinetics, efficacy, and safety of simvastatin. Care is warranted if simvastatin is prescribed to patients carrying decreased function SLCO1B1 or CYP3A4 alleles.

A pharmacogenomic assessment of psychiatric adverse drug reactions to levetiracetam

OBJECTIVE

Levetiracetam (LEV) is an effective antiseizure medicine, but 10%-20% of people treated with LEV report psychiatric side-effects, and up to 1% may have psychotic episodes. Pharmacogenomic predictors of these adverse drug reactions (ADRs) have yet to be identified. We sought to determine the contribution of both common and rare genetic variation to psychiatric and behavioral ADRs associated with LEV.

METHODS

This case-control study compared cases of LEV-associated behavioral disorder (n = 149) or psychotic reaction (n = 37) to LEV-exposed people with no history of psychiatric ADRs (n = 920). All samples were of European ancestry. We performed genome-wide association study (GWAS) analysis comparing those with LEV ADRs to controls. We estimated the polygenic risk scores (PRS) for schizophrenia and compared cases with LEV-associated psychotic reaction to controls. Rare variant burden analysis was performed using exome sequence data of cases with psychotic reactions (n = 18) and controls (n = 122).

RESULTS

Univariate GWAS found no significant associations with either LEV-associated behavioural disorder or LEV-psychotic reaction. PRS analysis showed that cases of LEV-associated psychotic reaction had an increased PRS for schizophrenia relative to contr ols (p = .0097, estimate = .4886). The rare-variant analysis found no evidence of an increased burden of rare genetic variants in people who had experienced LEV-associated psychotic reaction relative to controls.

SIGNIFICANCE

The polygenic burden for schizophrenia is a risk factor for LEV-associated psychotic reaction. To assess the clinical utility of PRS as a predictor, it should be tested in an independent and ideally prospective cohort. Larger sample sizes are required for the identification of significant univariate common genetic signals or rare genetic signals associated with psychiatric LEV ADRs.

The Clinical Pharmacogenetics Implementation Consortium Guideline for SLCO1B1, ABCG2, and CYP2C9 genotypes and Statin-Associated Musculoskeletal Symptoms

Statins reduce cholesterol, prevent cardiovascular disease, and are among the most commonly prescribed medications in the world. Statin-associated musculoskeletal symptoms (SAMS) impact statin adherence and ultimately can impede the long-term effectiveness of statin therapy. There are several identified pharmacogenetic variants that impact statin disposition and adverse events during statin therapy. SLCO1B1 encodes a transporter (SLCO1B1; alternative names include OATP1B1 or OATP-C) that facilitates the hepatic uptake of all statins. ABCG2 encodes an efflux transporter (BCRP) that modulates the absorption and disposition of rosuvastatin. CYP2C9 encodes a phase I drug metabolizing enzyme responsible for the oxidation of some statins. Genetic variation in each of these genes alters systemic exposure to statins (i.e., simvastatin, rosuvastatin, pravastatin, pitavastatin, atorvastatin, fluvastatin, lovastatin), which can increase the risk for SAMS. We summarize the literature supporting these associations and provide therapeutic recommendations for statins based on SLCO1B1, ABCG2, and CYP2C9 genotype with the goal of improving the overall safety, adherence, and effectiveness of statin therapy. This document replaces the 2012 and 2014 Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines for SLCO1B1 and simvastatin-induced myopathy.

Combining familial hypercholesterolemia and statin genetic studies as a strategy for the implementation of pharmacogenomics. A multidisciplinary approach

The diagnostic process of familial hypercholesterolemia frequently involves the use of genetic studies. Patients are treated with lipid-lowering drugs, frequently statins. Although pharmacogenomic clinical practice guidelines focusing on genotype-based statin prescription have been published, their use in routine clinical practice remains very modest.We have implemented a new NGS strategy that combines a panel of genes related to familial hypercholesterolemia with genomic regions related to the pharmacogenomics of lipid-lowering drugs described in clinical practice guidelines and in EMA and FDA drug labels. A multidisciplinary team of doctors, biologists, and pharmacists creates a clinical report that provides diagnostic and therapeutic findings using a knowledge management and clinical decision support system, as well as an algorithm for treatment selection.For 12 months, a total of 483 genetic diagnostic studies for familial hypercholesterolemia were carried out, of which 221 (45.8%) requested a complementary pharmacogenomic test. Of these 221 patients, 66.5% were carriers of actionable variants in any of the studied pharmacogenomic pathways: 46.6% of patients in one pathway, 19.0% in two pathways, and 0.9% in three pathways. 45.7% of patients could have a response to atorvastatin different from that of the reference population, 45.7% for simvastatin and lovastatin, 29.0% for fluvastatin, and 6.7% patients for pitavastatin.This implementation approach facilitates the incorporation of pharmacogenomic studies in clinical care practice, it does not add complexity nor additional steps to laboratory processes, and improves the pharmacotherapeutic process of patients.

Interaction of Antifungal Drugs with CYP3A- and OATP1B-Mediated Venetoclax Elimination

Venetoclax, a BCL-2 inhibitor used to treat certain hematological cancers, exhibits low oral bioavailability and high interpatient pharmacokinetic variability. Venetoclax is commonly administered with prophylactic antifungal drugs that may result in drug interactions, of which the underlying mechanisms remain poorly understood. We hypothesized that antifungal drugs may increase venetoclax exposure through inhibition of both CYP3A-mediated metabolism and OATP1B-mediated transport. Pharmacokinetic studies were performed in wild-type mice and mice genetically engineered to lack all CYP3A isoforms, or OATP1B2 that received venetoclax alone or in combination with ketoconazole or micafungin. In mice lacking all CYP3A isoforms, venetoclax AUC was increased by 1.8-fold, and pretreatment with the antifungal ketoconazole further increased venetoclax exposure by 1.6-fold, despite the absence of CYP3A. Ensuing experiments demonstrated that the deficiency of OATP1B-type transporters is also associated with increases in venetoclax exposure, and that many antifungal drugs, including micafungin, posaconazole, and isavuconazole, are inhibitors of this transport mechanism both in vitro and in vivo. These studies have identified OATP1B-mediated transport as a previously unrecognized contributor to the elimination of venetoclax that is sensitive to inhibition by various clinically-relevant antifungal drugs. Additional consideration is warranted when venetoclax is administered together with agents that inhibit both CYP3A-mediated metabolism and OATP1B-mediated transport.

Unexpected severe hepatic and skin toxicities during high dose methotrexate course for osteosarcoma

INTRODUCTION

high dose methotrexate (HD-MTX) regimen is used in osteosarcoma, leukemia and lymphoma treatment. Osteosarcoma is mostly diagnosed in children and adolescents. Most frequent methotrexate toxicities are mucositis, myelosuppression, renal failure, hepatitis and necrotizing encephalopathy. Toxicities increase with renal impairment, denutrition, in older patients, with some pharmacogenetics factors or with drug interactions.

CASE REPORT

We report a 16th years old woman diagnosed with osteosarcoma and experienced an unexpected severe hepatic and skin toxicities as toxic epidermal necrolys, Steven Johnson syndrome.

MANAGEMENT AND OUTCOME

This toxicity occurred despite acid folinic rescue performed as good practice recommendation. Fourteen hours after methotrexate administration, renal failure was observed and after 72 h an erythematous rash and epidermal detachment with toxic epidermal necrolys. Seven days after methotrexate administration, hepatic failure began until grade IV cytolysis. High dose of folinic acid were administered during all severe toxicities. Methotrexate were not longer administered to this young patient and chemotherapy with ifosfamide (IFO), doxorubicine and cisplatin were performed in this patient and complete histologic response were observed in the surgical bone resection.

DISCUSSION

No classical toxicities risk factors were identified in this patient but a homozygote mutation of MTHFR gene and homozygote SLCO1B1 gene mutation were found. MTHFR and SLCO1B1 are both implicated in methotrexate metabolism.