Catechol-O-methyltransferase gene (COMT) is associated with neurocognitive functioning in patients with sickle cell disease

PURPOSE

Neurocognitive impairment is a common and debilitating complication of sickle cell disease (SCD) resulting from a combination of biological and environmental factors. The catechol-O-methyltransferase (COMT) gene modulates levels of dopamine availability in the prefrontal cortex. COMT has repeatedly been implicated in the perception of pain stimuli and frequency of pain crises in patients with SCD and is known to be associated with neurocognitive functioning in the general population. The current study aimed to examine the associations of genetic variants in COMT and neurocognitive functioning in patients with SCD.

PATIENTS AND METHODS

The Sickle Cell Clinical Research and Intervention Program (SCCRIP) longitudinal cohort was used as a discovery cohort (n = 166). The genotypes for 5 SNPs (rs6269, rs4633, rs4818, rs4680, and rs165599) in COMT were extracted from whole genome sequencing data and analyzed using a dominant model. A polygenic score for COMT (PGSCOMT) integrating these 5 SNPs was analyzed as a continuous variable. The Cooperative Study of Sickle Cell Disease (CSSCD, n = 156) and the Silent Cerebral Infarction Transfusion (SIT, n = 114) Trial were used as 2 independent replication cohorts. Due to previously reported sex differences, all analyses were conducted separately in males and females. The Benjamini and Hochberg approach was used to calculate false discovery rate adjusted p-value (q-value).

RESULTS

In SCCRIP, 1 out of 5 SNPs (rs165599) was associated with IQ at q<0.05 in males but not females, and 2 other SNPs (rs4633 and rs4680) were marginally associated with sustained attention at p<0.05 in males only but did not maintain at q<0.05. PGSCOMT was negatively associated with IQ and sustained attention at p<0.05 in males only. Using 3 cohorts' data, 4 out of 5 SNPs (rs6269, rs4633, rs4680, rs165599) were associated with IQ (minimum q-value = 0.0036) at q<0.05 among male participants but not female participants. The PGSCOMT was negatively associated with IQ performance among males but not females across all cohorts.

CONCLUSION

Select COMT SNPs are associated with neurocognitive abilities in males with SCD. By identifying genetic predictors of neurocognitive performance in SCD, it may be possible to risk-stratify patients from a young age to guide implementation of early interventions.

Genetic Variation and Sickle Cell Disease Severity: A Systematic Review and Meta-Analysis

Importance

Sickle cell disease (SCD) is a monogenic disorder, yet clinical outcomes are influenced by additional genetic factors. Despite decades of research, the genetics of SCD remain poorly understood.

Objective

To assess all reported genetic modifiers of SCD, evaluate the design of associated studies, and provide guidelines for future analyses according to modern genetic study recommendations.

Data Sources

PubMed, Web of Science, and Scopus were searched through May 16, 2023, identifying 5290 publications.

Study Selection

At least 2 reviewers identified 571 original, peer-reviewed English-language publications reporting genetic modifiers of human SCD phenotypes, wherein the outcome was not treatment response, and the comparison was not between SCD subtypes or including healthy controls.

Data Extraction and Synthesis

Data relevant to all genetic modifiers of SCD were extracted, evaluated, and presented following STREGA and PRISMA guidelines. Weighted z score meta-analyses and pathway analyses were conducted.

Main Outcomes and Measures

Outcomes were aggregated into 25 categories, grouped as acute complications, chronic conditions, hematologic parameters or biomarkers, and general or mixed measures of SCD severity.

Results

The 571 included studies reported on 29 670 unique individuals (50% ≤ 18 years of age) from 43 countries. Of the 17 757 extracted results (4890 significant) in 1552 genes, 3675 results met the study criteria for meta-analysis: reported phenotype and genotype, association size and direction, variability measure, sample size, and statistical test. Only 173 results for 62 associations could be cross-study combined. The remaining associations could not be aggregated because they were only reported once or methods (eg, study design, reporting practice) and genotype or phenotype definitions were insufficiently harmonized. Gene variants regulating fetal hemoglobin and α-thalassemia (important markers for SCD severity) were frequently identified: 19 single-nucleotide variants in BCL11A, HBS1L-MYB, and HBG2 were significantly associated with fetal hemoglobin (absolute value of Z = 4.00 to 20.66; P = 8.63 × 10-95 to 6.19 × 10-5), and α-thalassemia deletions were significantly associated with increased hemoglobin level and reduced risk of albuminuria, abnormal transcranial Doppler velocity, and stroke (absolute value of Z = 3.43 to 5.16; P = 2.42 × 10-7 to 6.00 × 10-4). However, other associations remain unconfirmed. Pathway analyses of significant genes highlighted the importance of cellular adhesion, inflammation, oxidative and toxic stress, and blood vessel regulation in SCD (23 of the top 25 Gene Ontology pathways involve these processes) and suggested future research areas.

Conclusions and Relevance

The findings of this comprehensive systematic review and meta-analysis of all published genetic modifiers of SCD indicated that implementation of standardized phenotypes, statistical methods, and reporting practices should accelerate discovery and validation of genetic modifiers and development of clinically actionable genetic profiles.

Discovery of novel predisposing coding and noncoding variants in familial Hodgkin lymphoma

Familial aggregation of Hodgkin lymphoma (HL) has been demonstrated in large population studies, pointing to genetic predisposition to this hematological malignancy. To understand the genetic variants associated with the development of HL, we performed whole genome sequencing on 234 individuals with and without HL from 36 pedigrees that had 2 or more first-degree relatives with HL. Our pedigree selection criteria also required at least 1 affected individual aged <21 years, with the median age at diagnosis of 21.98 years (3-55 years). Family-based segregation analysis was performed for the identification of coding and noncoding variants using linkage and filtering approaches. Using our tiered variant prioritization algorithm, we identified 44 HL-risk variants in 28 pedigrees, of which 33 are coding and 11 are noncoding. The top 4 recurrent risk variants are a coding variant in KDR (rs56302315), a 5' untranslated region variant in KLHDC8B (rs387906223), a noncoding variant in an intron of PAX5 (rs147081110), and another noncoding variant in an intron of GATA3 (rs3824666). A newly identified splice variant in KDR (c.3849-2A>C) was observed for 1 pedigree, and high-confidence stop-gain variants affecting IRF7 (p.W238∗) and EEF2KMT (p.K116∗) were also observed. Multiple truncating variants in POLR1E were found in 3 independent pedigrees as well. Whereas KDR and KLHDC8B have previously been reported, PAX5, GATA3, IRF7, EEF2KMT, and POLR1E represent novel observations. Although there may be environmental factors influencing lymphomagenesis, we observed segregation of candidate germline variants likely to predispose HL in most of the pedigrees studied.

Genetic Analysis of Lung Cancer and the Germline Impact on Somatic Mutation Burden

BACKGROUND

Germline genetic variation contributes to lung cancer (LC) susceptibility. Previous genome-wide association studies (GWAS) have implicated susceptibility loci involved in smoking behaviors and DNA repair genes, but further work is required to identify susceptibility variants.

METHODS

To identify LC susceptibility loci, a family history-based genome-wide association by proxy (GWAx) of LC (48 843 European proxy LC patients, 195 387 controls) was combined with a previous LC GWAS (29 266 patients, 56 450 controls) by meta-analysis. Colocalization was used to explore candidate genes and overlap with existing traits at discovered susceptibility loci. Polygenic risk scores (PRS) were tested within an independent validation cohort (1 666 LC patients vs 6 664 controls) using variants selected from the LC susceptibility loci and a novel selection approach using published GWAS summary statistics. Finally, the effects of the LC PRS on somatic mutational burden were explored in patients whose tumor resections have been profiled by exome (n = 685) and genome sequencing (n = 61). Statistical tests were 2-sided.

RESULTS

The GWAx-GWAS meta-analysis identified 8 novel LC loci. Colocalization implicated DNA repair genes (CHEK1), metabolic genes (CYP1A1), and smoking propensity genes (CHRNA4 and CHRNB2). PRS analysis demonstrated that these variants, as well as subgenome-wide significant variants related to expression quantitative trait loci and/or smoking propensity, assisted in LC genetic risk prediction (odds ratio = 1.37, 95% confidence interval = 1.29 to 1.45; P < .001). Patients with higher genetic PRS loads of smoking-related variants tended to have higher mutation burdens in their lung tumors.

CONCLUSIONS

This study has expanded the number of LC susceptibility loci and provided insights into the molecular mechanisms by which these susceptibility variants contribute to LC development.

A genome-wide association study of obstructive heart defects among participants in the National Birth Defects Prevention Study

Obstructive heart defects (OHDs) share common structural lesions in arteries and cardiac valves, accounting for ~25% of all congenital heart defects. OHDs are highly heritable, resulting from interplay among maternal exposures, genetic susceptibilities, and epigenetic phenomena. A genome-wide association study was conducted in National Birth Defects Prevention Study participants (N_discovery  = 3978; N_replication  = 2507), investigating the genetic architecture of OHDs using transmission/disequilibrium tests (TDT) in complete case-parental trios (N_{discovery_TDT}  = 440; N_{replication_TDT}  = 275) and case-control analyses separately in infants (N_{discovery_CCI}  = 1635; N_{replication_CCI}  = 990) and mothers (case status defined by infant; N_{discovery_CCM}  = 1703; N_{replication_CCM}  = 1078). In the TDT analysis, the SLC44A2 single nucleotide polymorphism (SNP) rs2360743 was significantly associated with OHD (p_discovery  = 4.08 × 10^-9 ; p_replication  = 2.44 × 10^-4 ). A CAPN11 SNP (rs55877192) was suggestively associated with OHD (p_discovery  = 1.61 × 10^-7 ; p_replication  = 0.0016). Two other SNPs were suggestively associated (p < 1 × 10^-6 ) with OHD in only the discovery sample. In the case-control analyses, no SNPs were genome-wide significant, and, even with relaxed thresholds ( × _discovery  < 1 × 10^-5 and p_replication  < 0.05), only one SNP (rs188255766) in the infant analysis was associated with OHDs (p_discovery  = 1.42 × 10^-6 ; p_replication  = 0.04). Additional SNPs with p_discovery  < 1 × 10^-5 were in loci supporting previous findings but did not replicate. Overall, there was modest evidence of an association between rs2360743 and rs55877192 and OHD and some evidence validating previously published findings.

Fetal hemoglobin modulates neurocognitive performance in sickle cell anemia✰,✰✰

PURPOSE OF THE STUDY

Fetal hemoglobin (HbF) is a modifier of the clinical and hematologic phenotype of sickle cell anemia (SCA). Three quantitative trait loci (QTL) modulate HbF expression. The neurocognitive effects of variants in these QTL have yet to be explored. We evaluated the relation between 11 SNPs in the three HbF QTL: BCL11A, MYB, the HBB gene cluster, and full-scale intelligence (IQ) in SCA.

PATIENTS AND METHODS

The prospective longitudinal cohort study, Sickle Cell Clinical Research and Intervention Program, was used as a discovery cohort (n = 166). The genotypes for 11 SNPs were extracted through whole genome sequencing and were analyzed using an additive model. A polygenic score for HbF (PGS^HbF) integrating the numbers of low HbF alleles from 11 SNPs was analyzed as a continuous variable. The Cooperative Study of Sickle Cell Disease (n = 156) and the Silent Cerebral Infarction Transfusion (n = 114) Trial were used as two independent replication cohorts. Benjamini and Hochberg approach was used to calculate false discovery rate adjusted p-value (pFDR).

RESULTS

HbF was positively associated with IQ (minimum raw p = 0·0018) at pFDR<0·05. HbF mediated the relationship between two BCL11A SNPs, rs1427407 and rs7606173, HBS1L-MYB: rs9494142, and PGS^HbF with IQ (minimum raw p = 0·0035) at pFDR<0·05.

CONCLUSION

As the major modulator of the severity of SCA, HbF also influences neurocognition, which is done through mediation of its QTL. These findings have implications for early identification of neurocognitive risk and targeted intervention.

Generalization of a genetic risk score for time to first albuminuria in children with sickle cell anaemia: SCCRIP cohort study results

Albuminuria predicts kidney disease progression in individuals with sickle cell anaemia (SCA); however, earlier prediction of kidney disease with introduction of reno-protective therapies prior to the onset of albuminuria may attenuate disease progression. A genetic risk score (GRS) for SCA-related nephropathy may provide an improved one-time test for early identification of high-risk patients. We utilized a GRS from a recent, large, trans-ethnic meta-analysis to identify three single nucleotide polymorphisms that associate individually and in a GRS with time to first albuminuria episode in children with SCA.

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637.

Cross-cancer evaluation of polygenic risk scores for 16 cancer types in two large cohorts

Even distinct cancer types share biological hallmarks. Here, we investigate polygenic risk score (PRS)-specific pleiotropy across 16 cancers in European ancestry individuals from the Genetic Epidemiology Research on Adult Health and Aging cohort (16,012 cases, 50,552 controls) and UK Biobank (48,969 cases, 359,802 controls). Within cohorts, each PRS is evaluated in multivariable logistic regression models against all other cancer types. Results are then meta-analyzed across cohorts. Ten positive and one inverse cross-cancer associations are found after multiple testing correction. Two pairs show bidirectional associations; the melanoma PRS is positively associated with oral cavity/pharyngeal cancer and vice versa, whereas the lung cancer PRS is positively associated with oral cavity/pharyngeal cancer, and the oral cavity/pharyngeal cancer PRS is inversely associated with lung cancer. Overall, we validate known, and uncover previously unreported, patterns of pleiotropy that have the potential to inform investigations of risk prediction, shared etiology, and precision cancer prevention strategies.

Pan-cancer analysis demonstrates that integrating polygenic risk scores with modifiable risk factors improves risk prediction

Cancer risk is determined by a complex interplay of environmental and heritable factors. Polygenic risk scores (PRS) provide a personalized genetic susceptibility profile that may be leveraged for disease prediction. Using data from the UK Biobank (413,753 individuals; 22,755 incident cancer cases), we quantify the added predictive value of integrating cancer-specific PRS with family history and modifiable risk factors for 16 cancers. We show that incorporating PRS measurably improves prediction accuracy for most cancers, but the magnitude of this improvement varies substantially. We also demonstrate that stratifying on levels of PRS identifies significantly divergent 5-year risk trajectories after accounting for family history and modifiable risk factors. At the population level, the top 20% of the PRS distribution accounts for 4.0% to 30.3% of incident cancer cases, exceeding the impact of many lifestyle-related factors. In summary, this study illustrates the potential for improving cancer risk assessment by integrating genetic risk scores.

The landscape of host genetic factors involved in immune response to common viral infections

BACKGROUND

Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities.

METHODS

We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort.

RESULTS

Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRβ1 at positions 11, 13, 71, and 74 for Epstein-Barr virus (EBV), Varicella zoster virus (VZV), human herpesvirus 7, (HHV7), and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P < 5.0 × 10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, and CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P = 5.0 × 10-15 (MCV), NTN5: P = 1.1 × 10-9 (BKV), and P2RY13: P = 1.1 × 10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases, from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions.

CONCLUSIONS

Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

The landscape of host genetic factors involved in immune response to common viral infections

INTRODUCTION

Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities.

METHODS

We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort.

RESULTS

Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRβ1 at positions 11, 13, 71, and 74 for Epstein-Barr Virus (EBV), Varicella Zoster Virus (VZV), Human Herpes virus 7, (HHV7) and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P<5.×10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, as well as CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P=5.0×10-15 (MCV), NTN5: P=1.1×10-9 (BKV), and P2RY13: P=1.1×10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases; from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions.

CONCLUSIONS

Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

Pan-cancer study detects genetic risk variants and shared genetic basis in two large cohorts

Deciphering the shared genetic basis of distinct cancers has the potential to elucidate carcinogenic mechanisms and inform broadly applicable risk assessment efforts. Here, we undertake genome-wide association studies (GWAS) and comprehensive evaluations of heritability and pleiotropy across 18 cancer types in two large, population-based cohorts: the UK Biobank (408,786 European ancestry individuals; 48,961 cancer cases) and the Kaiser Permanente Genetic Epidemiology Research on Adult Health and Aging cohorts (66,526 European ancestry individuals; 16,001 cancer cases). The GWAS detect 21 genome-wide significant associations independent of previously reported results. Investigations of pleiotropy identify 12 cancer pairs exhibiting either positive or negative genetic correlations; 25 pleiotropic loci; and 100 independent pleiotropic variants, many of which are regulatory elements and/or influence cross-tissue gene expression. Our findings demonstrate widespread pleiotropy and offer further insight into the complex genetic architecture of cross-cancer susceptibility.

Genomewide Meta-Analysis Validates a Role for S1PR1 in Microtubule Targeting Agent-Induced Sensory Peripheral Neuropathy

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 S1PR₁ ; 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); P_META  = 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 S1PR₁ 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 S1PR₁ as a therapeutic strategy for prevention and/or treatment of MTA-induced neuropathy.

Immune-mediated genetic pathways resulting in pulmonary function impairment increase lung cancer susceptibility

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

A Pharmacogenetic Prediction Model of Progression-Free Survival in Breast Cancer using Genome-Wide Genotyping Data from CALGB 40502 (Alliance)

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.