Microtubule targeting agents (MTAs) are anticancer therapies commonly prescribed for breast cancer and other solid tumors. Sensory peripheral neuropathy (PN) is the major dose‐limiting toxicity for MTAs and can limit clinical efficacy. The current pharmacogenomic study aimed to identify genetic variations that explain patient susceptibility and drive mechanisms underlying development of MTA‐induced PN. A meta‐analysis of genomewide association studies (GWAS) from two clinical cohorts treated with MTAs (Cancer and Leukemia Group B (CALGB) 40502 and CALGB 40101) was conducted using a Cox regression model with cumulative dose to first instance of grade 2 or higher PN. Summary statistics from a GWAS of European subjects (n = 469) in CALGB 40502 that estimated cause‐specific risk of PN were meta‐analyzed with those from a previously published GWAS of European ancestry (n = 855) from CALGB 40101 that estimated the risk of PN. Novel single nucleotide polymorphisms in an enhancer region downstream of sphingosine‐1‐phosphate receptor 1 (S1PR1 encoding S1PR1; e.g., rs74497159, βCALGB 40101 per allele log hazard ratio (95% confidence interval (CI)) = 0.591 (0.254–0.928), βCALGB 40502 per allele log hazard ratio (95% CI) = 0.693 (0.334–1.053); PMETA = 3.62 × 10−7) were the most highly ranked associations based on P values with risk of developing grade 2 and higher PN. In silico functional analysis identified multiple regulatory elements and potential enhancer activity for S1PR1 within this genomic region. Inhibition of S1PR1 function in induced pluripotent stem cell‐derived human sensory neurons shows partial protection against paclitaxel‐induced neurite damage. These pharmacogenetic findings further support ongoing clinical evaluations to target S1PR1 as a therapeutic strategy for prevention and/or treatment of MTA‐induced neuropathy.

Genome-wide genotyping data are increasingly available for pharmacogenetic association studies, but application of these data for development of prediction models is limited. Prediction methods, such as elastic net regularization, have recently been applied to genetic studies but only limitedly to pharmacogenetic outcomes. An elastic net was applied to a pharmacogenetic study of progression-free survival (PFS) of 468 patients with advanced breast cancer in a clinical trial of paclitaxel, nab-paclitaxel, and ixabepilone. A final model included 13 single nucleotide polymorphisms (SNPs) in addition to clinical covariates (prior taxane status, hormone receptor status, disease-free interval, and presence of visceral metastases) with an area under the curve (AUC) integrated over time of 0.81, an increase compared to an AUC of 0.64 for a model with clinical covariates alone. This model may be of value in predicting PFS with microtubule targeting agents and may inform reverse translational studies to understand differential response to these drugs.Applications of genome-wide data beyond typical variant-byvariant association testing have been gaining popularity for pharmacogenetic phenotypes. Examples include gene-based or region-based tests, 1,2 heritability analyses, 3,4 transcriptome-based analyses, 5-7 and polygenic risk scores. 8-10 Another method is meta-analysis, combining summary statistics of multiple studies,

Microtubule targeting agents (MTAs) are anticancer therapies commonly prescribed for breast cancer and other solid tumors. Sensory peripheral neuropathy (PN) is the major dose-limiting toxicity for MTAs and can limit clinical efficacy. The current pharmacogenomic study aimed to identify genetic variations that explain patient susceptibility and drive mechanisms underlying development of MTA-induced PN. A meta-analysis of genome-wide association studies (GWAS) from two clinical cohorts treated with MTAs (CALGB 40502 and CALGB 40101) was conducted using a Cox regression model with cumulative dose to first instance of grade 2 or higher PN. Summary statistics from a GWAS of European subjects (n = 469) in CALGB 40502 that estimated cause-specific risk of PN were meta-analyzed with those from a previously published GWAS of European ancestry (n = 855) from CALGB 40101 that estimated the risk of PN. Novel single nucleotide polymorphisms in an enhancer region downstream of sphingosine-1-phosphate receptor 1 (S1PR1 encoding S1PR1; e.g., rs74497159, βCALGB 40101 per allele log hazard ratio (95% CI) = 0.591 (0.254 - 0.928), βCALGB 40502 per allele log hazard ratio (95% CI) = 0.693 (0.334 - 1.053); PMETA = 3.62x10-7) were the most highly ranked associations based on P-values with risk of developing grade 2 and higher PN. In silico functional analysis identified multiple regulatory elements and potential enhancer activity for S1PR1 within this genomic region. Inhibition of S1PR1 function in iPSC-derived human sensory neurons shows partial protection against paclitaxel-induced neurite damage. These pharmacogenetic findings further support ongoing clinical evaluations to target S1PR1 as a therapeutic strategy for prevention and/or treatment of MTA-induced neuropathy.