Prediction models of persistent taxane-induced peripheral neuropathy among breast cancer survivors using whole-exome sequencing.

Persistent taxane-induced peripheral neuropathy (TIPN) is highly prevalent among early-stage breast cancer survivors (ESBCS) and has detrimental effect on quality of life. We leveraged logistic regression models to develop and validate polygenic prediction models to estimate the risk of persistent PN symptoms in a training cohort and validation cohort taking clinical risk factors into account. Based on 337 whole-exome sequenced ESBCS two of five prediction models for individual PN symptoms obtained AUC results above 60% when validated. Using the model for numbness in feet (35 SNVs) in the test cohort, 73% survivors were correctly predicted. For tingling in feet (55 SNVs) 70% were correctly predicted. Both models included SNVs from the ADAMTS20, APT6V0A2, CCDC88C, CYP2C8, EPHA5, NR1H3, PSKH2/APTV0D2, and SCN10A genes. For cramps in feet, difficulty climbing stairs and difficulty opening a jar the validation was unsuccessful. Polygenic prediction models including clinical risk factors can estimate the risk of persistent taxane-induced numbness in feet and tingling in feet in ESBCS.

Enhancer mutations modulate the severity of chemotherapy-induced myelosuppression.

Non-small cell lung cancer is often diagnosed at advanced stages, and many patients are still treated with classical chemotherapy. The unselective nature of chemotherapy often results in severe myelosuppression. Previous studies showed that protein-coding mutations could not fully explain the predisposition to myelosuppression. Here, we investigate the possible role of enhancer mutations in myelosuppression susceptibility. We produced transcriptome and promoter-interaction maps (using HiCap) of three blood stem-like cell lines treated with carboplatin or gemcitabine. Taking advantage of publicly available enhancer datasets, we validated HiCap results in silico and in living cells using epigenetic CRISPR technology. We also developed a network approach for interactome analysis and detection of differentially interacting genes. Differential interaction analysis provided additional information on relevant genes and pathways for myelosuppression compared with differential gene expression analysis at the bulk level. Moreover, we showed that enhancers of differentially interacting genes are highly enriched for variants associated with differing levels of myelosuppression. Altogether, our work represents a prominent example of integrative transcriptome and gene regulatory datasets analysis for the functional annotation of noncoding mutations.

Early changes in gene expression profiles in AML patients during induction chemotherapy.

BACKGROUND: Elucidation of the genetic mechanisms underlying treatment response to standard induction chemotherapy in AML patients is warranted, in order to aid in risk-adapted treatment decisions as novel treatments are emerging. In this pilot study, we explored the treatment-induced expression patterns in a small cohort of AML patients by analyzing differential gene expression (DGE) over the first 2 days of induction chemotherapy. METHODS: Blood samples were collected from ten AML patients at baseline (before treatment initiation) and during the first 2 days of treatment (Day 1; approximately 24 h, and Day 2; approximately 48 h after treatment initiation, respectively) and RNA was extracted for subsequent RNA sequencing. DGE between time points were assessed by pairwise analysis using the R package edgeR version 3.18.1 in all patients as well as in relation to treatment response (complete remission, CR, vs non-complete remission, nCR). Ingenuity Pathway Analysis (Qiagen) software was used for pathway analysis and visualization. RESULTS: After initial data quality control, two patients were excluded from further analysis, resulting in a final cohort of eight patients with data from all three timepoints. DGE analysis demonstrated activation of pathways with genes directly or indirectly associated with NF-κB signaling. Significant activation of the NF-κB pathway was seen in 50% of the patients 2 days after treatment start, while iNOS pathway effects could be identified already after 1 day. nCR patients displayed activation of pathways associated with cell cycle progression, oncogenesis and anti-apoptotic behavior, including the STAT3 pathway and Salvage pathways of pyrimidine ribonucleotides. Notably, a significant induction of cytidine deaminase, an enzyme responsible for the deamination of Ara-C, could be observed between baseline and Day 2 in the nCR patients but not in patients achieving CR. CONCLUSIONS: In conclusion, we show that time-course analysis of gene expression represents a feasible approach to identify relevant pathways affected by standard induction chemotherapy in AML patients. This poses as a potential method for elucidating new drug targets and biomarkers for categorizing disease aggressiveness and evaluating treatment response. However, more studies on larger cohorts are warranted to elucidate the transcriptional basis for drug response.

IFNγ directly counteracts imatinib-induced apoptosis of primary human CD34+ CML stem/progenitor cells potentially through the upregulation of multiple key survival factors.

Tyrosine kinase inhibitors (TKIs) have dramatically improved the survival in chronic myeloid leukemia (CML), but residual disease typically persists even after prolonged treatment. Several lines of evidence suggest that TKIs administered to CML patients upregulate interferon γ (IFNγ) production, which may counteract the anti-tumorigenic effects of the therapy. We now show that activated T cell-conditioned medium (TCM) enhanced proliferation and counteracted imatinib-induced apoptosis of CML cells, and addition of a neutralizing anti-IFNγ antibody at least partially inhibited the anti-apoptotic effect. Likewise, recombinant IFNγ also reduced imatinib-induced apoptosis of CML cells. This anti-apoptotic effect of IFNγ was independent of alternative IFNγ signaling pathways, but could be notably diminished by STAT1-knockdown. Furthermore, IFNγ upregulated the expression of several anti-apoptotic proteins, including MCL1, PARP9, and PARP14, both in untreated and imatinib-treated primary human CD34+ CML stem/progenitor cells. Our results suggest that activated T cells in imatinib-treated CML patients can directly rescue CML cells from imatinib-induced apoptosis at least partially through the secretion of IFNγ, which exerts a rapid, STAT1-dependent anti-apoptotic effect potentially through the simultaneous upregulation of several key hematopoietic survival factors. These mechanisms may have a major clinical impact, when targeting residual leukemic stem/progenitor cells in CML.

The association of four genetic variants with myelosuppression in gemcitabine-treated Japanese is not evident in gemcitabine/carboplatin-treated Swedes.

Gemcitabine/carboplatin-induced myelosuppressive adverse drug reactions (ADRs) are clinical problems leading to patient suffering and dose alterations. There is a need for personalised medicine to improve treatment effects and patients' well-being. We tested four genetic variants, rs11141915, rs1901440, rs12046844 and rs11719165, previously suggested as potential biomarkers for gemcitabine-induced leukopenia/neutropenia in Japanese patients, in 213 Swedish gemcitabine/carboplatin-treated non-small cell lung cancer (NSCLC) patients. DNA was genotyped using TaqMan probes and real-time PCR. The relationships between the risk alleles and low toxicity (non-ADR: Common Terminology Criteria for Adverse Events [CTCAE] grades 0) or high toxicity (ADR: CTCAE grades 3-4) of platelets, leukocytes and neutrophils were evaluated using Fisher's exact test. The risk alleles did not correlate with myelosuppression, and the strongest borderline significance (not withstanding adjustment for multiple testing) was for rs1901440 (neutropenia, p = 0.043) and rs11719165 (leukopenia, p = 0.049) where the risk alleles trended towards lower toxicity, contrasting with previous study findings. Risk alleles and higher risk scores were more common among our patients. We conclude that the genetic variants do not apply to Swedish patients treated with gemcitabine/carboplatin. However, they can still be important in other populations and cohorts, especially in a gemcitabine monotherapy setting, where the causal genetic variation might influence myelosuppressive ADRs.

Combination of tyrosine kinase inhibitors and the MCL1 inhibitor S63845 exerts synergistic antitumorigenic effects on CML cells.

Tyrosine kinase inhibitor (TKI) treatment has dramatically improved the survival of chronic myeloid leukemia (CML) patients, but measurable residual disease typically persists. To more effectively eradicate leukemia cells, simultaneous targeting of BCR-ABL1 and additional CML-related survival proteins has been proposed. Notably, several highly specific myeloid cell leukemia 1 (MCL1) inhibitors have recently entered clinical trials for various hematologic malignancies, although not for CML, reflecting the insensitivity of CML cell lines to single MCL1 inhibition. Here, we show that combining TKI (imatinib, nilotinib, dasatinib, or asciminib) treatment with the small-molecule MCL1 inhibitor S63845 exerted strong synergistic antiviability and proapoptotic effects on CML lines and CD34+ stem/progenitor cells isolated from untreated CML patients in chronic phase. Using wild-type BCR-ABL1-harboring CML lines and their T315I-mutated sublines (generated by CRISPR/Cas9-mediated homologous recombination), we prove that the synergistic proapoptotic effect of the drug combination depended on TKI-mediated BCR-ABL1 inhibition, but not on TKI-related off-target mechanisms. Moreover, we demonstrate that colony formation of CML but not normal hematopoietic stem/progenitor cells became markedly reduced upon combination treatment compared to imatinib monotherapy. Our results suggest that dual targeting of MCL1 and BCR-ABL1 activity may efficiently eradicate residual CML cells without affecting normal hematopoietic stem/progenitors.

JMJD1C knockdown affects myeloid cell lines proliferation, viability, and gemcitabine/carboplatin-sensitivity.

OBJECTIVES: Chemotherapy-induced hematological toxicities are potentially life-threatening adverse drug reactions that vary between individuals. Recently, JMJD1C has been associated with gemcitabine/carboplatin-induced thrombocytopenia in non-small-cell lung cancer patients, making it a candidate marker for predicting the risk of toxicity. This study investigates if JMJD1C knockdown affects gemcitabine/carboplatin-sensitivity in cell lines. METHODS: Lentiviral transduction-mediated shRNA knockdown of JMJD1C in the cell lines K562 and MEG-01 were performed using shRNA#32 and shRNA#33. The knockdown was evaluated using qPCR. Cell proliferation, viability, and gemcitabine/carboplatin-sensitivity were subsequently determined using cell counts, trypan blue, and the MTT assay. RESULTS: ShRNA#33 resulted in JMJD1C downregulation by 56.24% in K562 and 68.10% in MEG-01. Despite incomplete knockdown, proliferation (reduction of cell numbers by 61-68%, day 7 post-transduction) and viability (reduction by 21-53%, day 7 post-transduction) were impaired in K562 and MEG-01 cells. Moreover, JMJD1C knockdown reduced the gemcitabine IC50-value for K562 cells (P < 0.01) and MEG-01 cells (P < 0.05) compared to scrambled shRNA control transduced cells. CONCLUSIONS: Our results suggest that JMJD1C is essential for proliferation, survival, and viability of K562 and MEG-01 cells. Further, JMJD1C also potentially affects the cells gemcitabine/carboplatin-sensitivity. Although further research is required, the findings show that JMJD1C could have an influential role for gemcitabine/carboplatin-sensitivity.

Activation of the µ-opioid receptor by alicyclic fentanyls: Changes from high potency full agonists to low potency partial agonists with increasing alicyclic substructure.

Fentanyl analogs represent an important group of new psychoactive substances and knowing their efficacy and potency might assist in interpreting observed concentrations. The potency of fentanyl analogs can be estimated from in vitro studies and can be used to establish structure-activity relationships. In this study, recombinant CHO-K1 cells (AequoScreen) expressing the human µ-opioid receptor were used to establish dose-response curves via luminescent analysis for cyclopropyl-, cyclobutyl-, cyclopentyl-, cyclohexyl-, and 2,2,3,3-tetramethylcyclopropylfentanyl (TMCPF), on three separate occasions, using eight different concentrations in an eight-fold serial dilution in triplicates starting at ~60 µM. Fentanyl was used as a full agonist reference while morphine and buprenorphine were included for comparison. Cyclopropylfentanyl (EC50 = 4.3 nM), cyclobutylfentanyl (EC50 = 6.2 nM), and cyclopentylfentanyl (EC50 = 13 nM) were full agonists slightly less potent than fentanyl (EC50 = 1.7 nM). Cyclohexylfentanyl (EC50 = 3.1 µM, efficacy 48%) and TMCPF (EC50 = 1.5 µM, efficacy 65%) were partial agonists less potent than morphine (EC50 = 430 nM). Based on the results, cyclopropyl-, cyclobutyl-, and cyclopentylfentanyl would be expected to induce intoxication or cause fatal poisonings at similar concentrations to fentanyl, while the toxic or fatal concentrations of cyclohexylfentanyl and TMCPF would be expected to be much higher.

Whole-genome sequencing and gene network modules predict gemcitabine/carboplatin-induced myelosuppression in non-small cell lung cancer patients.

Gemcitabine/carboplatin chemotherapy commonly induces myelosuppression, including neutropenia, leukopenia, and thrombocytopenia. Predicting patients at risk of these adverse drug reactions (ADRs) and adjusting treatments accordingly is a long-term goal of personalized medicine. This study used whole-genome sequencing (WGS) of blood samples from 96 gemcitabine/carboplatin-treated non-small cell lung cancer (NSCLC) patients and gene network modules for predicting myelosuppression. Association of genetic variants in PLINK found 4594, 5019, and 5066 autosomal SNVs/INDELs with p ≤ 1 × 10-3 for neutropenia, leukopenia, and thrombocytopenia, respectively. Based on the SNVs/INDELs we identified the toxicity module, consisting of 215 unique overlapping genes inferred from MCODE-generated gene network modules of 350, 345, and 313 genes, respectively. These module genes showed enrichment for differentially expressed genes in rat bone marrow, human bone marrow, and human cell lines exposed to carboplatin and gemcitabine (p < 0.05). Then using 80% of the patients as training data, random LASSO reduced the number of SNVs/INDELs in the toxicity module into a feasible prediction model consisting of 62 SNVs/INDELs that accurately predict both the training and the test (remaining 20%) data with high (CTCAE 3-4) and low (CTCAE 0-1) maximal myelosuppressive toxicity completely, with the receiver-operating characteristic (ROC) area under the curve (AUC) of 100%. The present study shows how WGS, gene network modules, and random LASSO can be used to develop a feasible and tested model for predicting myelosuppressive toxicity. Although the proposed model predicts myelosuppression in this study, further evaluation in other studies is required to determine its reproducibility, usability, and clinical effect.

Genetic association of gemcitabine/carboplatin-induced leukopenia and neutropenia in non-small cell lung cancer patients using whole-exome sequencing.

OBJECTIVES: Gemcitabine/carboplatin treatment is known to cause severe adverse drug reactions which can lead to the need for reduction or cessation of chemotherapy. It would be beneficial to identify patients at risk of severe hematological toxicity in advance before treatment start. This study aims to identify genetic markers for gemcitabine/carboplatin-induced leukopenia and neutropenia in non-small cell lung cancer patients. MATERIAL AND METHODS: Whole-exome sequencing was performed on 215 patients. Association analysis was performed on single-nucleotide variants (SNVs) and genes, and the validation was based on an independent genome-wide association study (GWAS). Based on the association and validation analyses the genetic variants were then selected for and used in weighted genetic risk score (wGRS) prediction models for leukopenia and neutropenia. RESULTS: Association analysis identified 50 and 111 SNVs, and 12 and 20 genes, for leukopenia and neutropenia, respectively. Of these SNVS 20 and 19 were partially validated for leukopenia and neutropenia, respectively. The genes SVIL (p = 2.48E-06) and EFCAB2 (p = 4.63E-06) were significantly associated with leukopenia contain the partially validated SNVs rs3740003, rs10160013, rs1547169, rs10927386 and rs10927387. The wGRS prediction models showed significantly different risk scores for high and low toxicity patients. CONCLUSION: We have identified and partially validated genetic biomarkers in SNVs and genes correlated to gemcitabine/carboplatin-induced leukopenia and neutropenia and created wGRS models for predicting the risk of chemotherapy-induced hematological toxicity. These results provide a strong foundation for further studies of chemotherapy-induced toxicity.

Single-Cell RNA Sequencing of Hematopoietic Stem and Progenitor Cells Treated with Gemcitabine and Carboplatin.

Treatments that include gemcitabine and carboplatin induce dose-limiting myelosuppression. The understanding of how human bone marrow is affected on a transcriptional level leading to the development of myelosuppression is required for the implementation of personalized treatments in the future. In this study, we treated human hematopoietic stem and progenitor cells (HSPCs) harvested from a patient with chronic myelogenous leukemia (CML) with gemcitabine/carboplatin. Thereafter, scRNA-seq was performed to distinguish transcriptional effects induced by gemcitabine/carboplatin. Gene expression was calculated and evaluated among cells within and between samples compared to untreated cells. Cell cycle analysis showed that the treatments effectively decrease cell proliferation, indicated by the proportion of cells in the G2M-phase dropping from 35% in untreated cells to 14.3% in treated cells. Clustering and t-SNE showed that cells within samples and between treated and untreated samples were affected differently. Enrichment analysis of differentially expressed genes showed that the treatments influence KEGG pathways and Gene Ontologies related to myeloid cell proliferation/differentiation, immune response, cancer, and the cell cycle. The present study shows the feasibility of using scRNA-seq and chemotherapy-treated HSPCs to find genes, pathways, and biological processes affected among and between treated and untreated cells. This indicates the possible gains of using single-cell toxicity studies for personalized medicine.

Genes and variants in hematopoiesis-related pathways are associated with gemcitabine/carboplatin-induced thrombocytopenia.

Chemotherapy-induced myelosuppression, including thrombocytopenia, is a recurrent problem during cancer treatments that may require dose alterations or cessations that could affect the antitumor effect of the treatment. To identify genetic markers associated with treatment-induced thrombocytopenia, we whole-exome sequenced 215 non-small cell lung cancer patients homogeneously treated with gemcitabine/carboplatin. The decrease in platelets (defined as nadir/baseline) was used to assess treatment-induced thrombocytopenia. Association between germline genetic variants and thrombocytopenia was analyzed at single-nucleotide variant (SNV) (based on the optimal false discovery rate, the severity of predicted consequence, and effect), gene, and pathway levels. These analyses identified 130 SNVs/INDELs and 25 genes associated with thrombocytopenia (P-value < 0.002). Twenty-three SNVs were validated in an independent genome-wide association study (GWAS). The top associations include rs34491125 in JMJD1C (P-value = 9.07 × 10-5), the validated variants rs10491684 in DOCK8 (P-value = 1.95 × 10-4), rs6118 in SERPINA5 (P-value = 5.83 × 10-4), and rs5877 in SERPINC1 (P-value = 1.07 × 10-3), and the genes CAPZA2 (P-value = 4.03 × 10-4) and SERPINC1 (P-value = 1.55 × 10-3). The SNVs in the top-scoring pathway "Factors involved in megakaryocyte development and platelet production" (P-value = 3.34 × 10-4) were used to construct weighted genetic risk score (wGRS) and logistic regression models that predict thrombocytopenia. The wGRS predict which patients are at high or low toxicity risk levels, for CTCAE (odds ratio (OR) = 22.35, P-value = 1.55 × 10-8), and decrease (OR = 66.82, P-value = 5.92 × 10-9). The logistic regression models predict CTCAE grades 3-4 (receiver operator characteristics (ROC) area under the curve (AUC) = 0.79), and large decrease (ROC AUC = 0.86). We identified and validated genetic variations within hematopoiesis-related pathways that provide a solid foundation for future studies using genetic markers for predicting chemotherapy-induced thrombocytopenia and personalizing treatments.

ABCB1 Variation Affects Myelosuppression, Progression-free Survival and Overall Survival in Paclitaxel/Carboplatin-treated Ovarian Cancer Patients.

The standard chemotherapy for ovarian cancer is paclitaxel/carboplatin. Patients often exhibit myelosuppressive toxicity, and the treatment response varies considerably. In this study, we investigated the previously reported SNPs 1199G>A (rs2229109), 1236C>T (rs1128503), 2677G>T/A (rs2032582), 3435C>T (rs1045642) in ABCB1, and 1196A>G (rs10509681) in CYP2C8 and their association with treatment-induced myelosuppression, progression-free survival (PFS) and overall survival (OS). From the phase III study, OAS-07OVA, 525 patients (All) treated with carboplatin and paclitaxel administered as Paclical (Arm A, n = 260) or Taxol® (Arm B, n = 265) were included and genotyped using pyrosequencing. Genotype associations with myelosuppression, PFS and OS were investigated using anova, Kaplan-Meier analysis and Cox proportional hazard models. The most prominent finding was for the ABCB1 variant 3435TT, which was significantly associated with increased PFS in All (hazard ratio (HR) = 0.623), in Arm A (HR = 0.590) and in Arm B (HR = 0.627), as well as increased OS in All (HR = 0.443) and in Arm A (HR = 0.372) compared to the wild-type, 3435CC. For toxicity, the most interesting finding concerned the haplotype, including 1236TT, 2677TT and 3435TT, which was associated with higher neutrophil values in Arm B (p = 0.039) and less neutrophil decrease in All (p = 0.048) and in Arm B (p = 0.021). It is noteworthy that the results varied depending on the treatment arm which indicates that the effects of ABCB1 variants vary with the treatment regimen. Our results reflect the contradictory results of previous studies, confirming that small variations in the composition of treatment regimens and patient populations may influence the interpretation of SNPs effects on treatment outcome.