Using Whole-Exome Sequencing to Identify Genetic Markers for Carboplatin and Gemcitabine-Induced Toxicities

Recent Advances in Drug Discovery Toxicology

Major advances in scientific discovery and insights that stem from the development and use of new techniques and models can bring remarkable progress to conventional toxicology. Although animal testing is still considered as the "gold standard" in traditional toxicity testing, there is a necessity for shift from animal testing to alternative methods regarding the drug safety testing owing to the emerging state-of-art techniques and the proposal of 3Rs (replace, reduce, and refine) towards animal welfare. This review describes some recent research methods in drug discovery toxicology, including in vitro cell and organ-on-a-chip, imaging systems, model organisms (C. elegans, Danio rerio, and Drosophila melanogaster), and toxicogenomics in modern toxicology testing.

Whole-exome sequencing reveals novel candidate single nucleotide variations for preventing adverse effects of levonorgestrel implantation

Aim: To identify novel genes associated with adverse effects of levonorgestrel (LNG) implants based on comparative whole-exome sequencing. Materials & methods: A cohort comprising 104 participants, including 52 controls and 52 women with LNG-related adverse effects, was recruited. Seven cases and eight controls were selected for whole-exome sequencing. We verified 13 single nucleotide variations (SNVs) related with integrin-mediated signaling pathway and cell proliferation using the MassARRAY platform. Results: Finally, we screened 49 cases and 52 controls for analyses. Two SNVs (rs7255721 and rs1042522) were located in ADAMTS10 and TP53, respectively, and significantly different between two groups. These two SNVs lead to changes in protein structure and physicochemical parameters. Conclusion: Here, we defined two pathogenic mutations related to adverse LNG effects.

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.

Genetic variation in the Mauritian cynomolgus macaque population reflects variation in the human population

The cynomolgus macaque is an important species for preclinical research, however the extent of genetic variation in this population and its similarity to the human population is not well understood. Exome sequencing was conducted for 101 cynomolgus macaques to characterize genetic variation. The variant distribution frequency was 7.81 variants per kilobase across the sequenced regions, with a total of 2,770,009 single nucleotide variants identified from 2,996,041 loci. A large portion (85.6%) had minor allele frequencies greater than 5%. Enriched pathways for genes with high genetic diversity (≥10 variants per kilobase) were those involving signaling peptides and immune response. Compared to human, the variant distribution frequency and nucleotide diversity in the macaque exome was approximately 4 times greater; however the ratio of non-synonymous to synonymous variants was similar (0.735 and 0.831, respectively). Understanding genetic variability in cynomolgus macaques will enable better interpretation and human translation of phenotypic variability in this species.

The influence of genetic polymorphisms on the toxicity of platinum-based chemotherapy in the treatment of non-small cell lung cancer

Lung cancer remains one of the most frequent and the deadliest of malignant diseases throughout the world. Target and immune therapy have revolutionalized the treatment of this disease, but platinum-based chemotherapy still has a place in the treatment algorithm. The toxicity profile of cisplatin is well known and can be a limiting factor in the adequate treatment delivery of the drug. There are important inter-individual differences in the efficacy and the toxicity of all chemotherapy drugs, which cannot be explained solely by the characteristics of the tumor. In order to define predictive factors for the occurrence of toxic effects, numerous genetic alterations have been investigated - especially single nucleotide polymorphisms (SNPs). The investigated genes are those involved in DNA repair mechanisms, signal pathways of apoptosis, DNA synthesis, transport mechanisms, but often with inconclusive and opposing results. It is clear that the effect of SNPs on the occurrence of cisplatin toxicity cannot be explained by investigating just one or several genes alone, but epigenetic interactions must be investigated, as well as interactions with outside factors. The study of SNPs is, however, a relatively simple and inexpensive method and, as such, can be used as one of the prognostic tools for everyday practice.

Current management of chemotherapy-induced neutropenia in adults: key points and new challenges: Committee of Neoplastic Supportive-Care (CONS), China Anti-Cancer Association Committee of Clinical Chemotherapy, China Anti-Cancer Association

Chemotherapy-induced neutropenia (CIN) is a potentially fatal and common complication in myelosuppressive chemotherapy. The timing and grade of CIN may play prognostic and predictive roles in cancer therapy. CIN is associated with older age, poor functional and nutritional status, the presence of significant comorbidities, the type of cancer, previous chemotherapy cycles, the stage of the disease, specific chemotherapy regimens, and combined therapies. There are many key points and new challenges in the management of CIN in adults including: (1) Genetic risk factors to evaluate the patient's risk for CIN remain unclear. However, these risk factors urgently need to be identified. (2) Febrile neutropenia (FN) remains one of the most common reasons for oncological emergency. No consensus nomogram for FN risk assessment has been established. (3) Different assessment tools [e.g., Multinational Association for Supportive Care in Cancer (MASCC), the Clinical Index of Stable Febrile Neutropenia (CISNE) score model, and other tools] have been suggested to help stratify the risk of complications in patients with FN. However, current tools have limitations. The CISNE score model is useful to support decision-making, especially for patients with stable FN. (4) There are still some challenges, including the benefits of granulocyte colony stimulating factor treatment and the optimal antibiotic regimen in emergency management of FN. In view of the current reports, our group discusses the key points, new challenges, and management of CIN.

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⁻³ 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.

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.

Genetic Polymorphisms and Platinum-Based Chemotherapy-Induced Toxicities in Patients With Lung Cancer: A Systematic Review and Meta-Analysis

Background: Platinum-based agents, including cisplatin, carboplatin, and oxaliplatin, are indispensable for the treatment of lung cancer. The development of toxicity frequently necessitates dose reduction or discontinuation of therapy, despite the clinical response. Pharmacogenomics studies were reviewed to identify the possible genetic variants that underlie individual susceptibility to platinum-related toxicities. Method: We conducted a systematic search in PubMed and Embase for pharmacogenomics reports that focused on commonly reported platinum-induced toxicities, such as gastrointestinal (GI), hematological, neurological, and other toxicities, in patients diagnosed with lung cancer. Meta-analyses were conducted to determine the association between genetic polymorphisms and platinum-induced toxicity by checking the odds ratio (OR) and 95% confidence interval (CI) using random or fixed-effects models as appropriate. Results: Twenty eligible studies that met the inclusion criteria with sufficient data were extracted and presented comprehensively. A total of 16 polymorphisms from 11 genes were included in the meta-analysis. MTHFR rs1801131 and MDM2 rs1690924 were significantly correlated with platinum-induced GI toxicity (P = 0.04 and P = 0.02, respectively). Patients with the MTHFR rs1801131AA and MDM2 rs1690924TC/CC genotype tended to have a higher risk of GI toxicity than patients with other genotypes did (OR = 1.73, 95% CI = 0.86-2.18; and OR = 0.51, 95% CI = 0.29-0.88, respectively). Compared to carriers of the MTHFR rs1801133CC genotype, carriers of the CT/TT genotype had a significantly increased risk of hematological toxicity (P = 0.01, OR = 1.68, 95% CI = 1.12-2.52). Conclusion: In the future, physicians should pay careful attention to MTHFR and MDM2 for personalized chemotherapy treatment among patients with lung cancer.

Metastatic adrenocortical carcinoma displays higher mutation rate and tumor heterogeneity than primary tumors

Adrenocortical cancer (ACC) is a rare cancer with poor prognosis and high mortality due to metastatic disease. All reported genetic alterations have been in primary ACC, and it is unknown if there is molecular heterogeneity in ACC. Here, we report the genetic changes associated with metastatic ACC compared to primary ACCs and tumor heterogeneity. We performed whole-exome sequencing of 33 metastatic tumors. The overall mutation rate (per megabase) in metastatic tumors was 2.8-fold higher than primary ACC tumor samples. We found tumor heterogeneity among different metastatic sites in ACC and discovered recurrent mutations in several novel genes. We observed 37–57% overlap in genes that are mutated among different metastatic sites within the same patient. We also identified new therapeutic targets in recurrent and metastatic ACC not previously described in primary ACCs.

Adrenocortical cancer (ACC) is a rarely diagnosed and aggressive cancer whose metastatic form has been scarcely studied. Here, the authors study primary and metastatic ACC to investigate genomic heterogeneity, discovering higher mutation rates in metastatic lesions and novel recurrent mutations.

Identification of a Genomic Region between SLC29A1 and HSP90AB1 Associated with Risk of Bevacizumab-Induced Hypertension: CALGB 80405 (Alliance)

Purpose: Bevacizumab is a VEGF-specific angiogenesis inhibitor indicated as an adjunct to chemotherapy for the treatment of multiple cancers. Hypertension is commonly observed during bevacizumab treatment, and high-grade toxicity can limit therapy or lead to cardiovascular complications. The factors that contribute to interindividual variability in blood pressure rise during bevacizumab treatment are not well understood.Experimental Design: To identify genomic regions associated with bevacizumab-induced hypertension risk, sequencing of candidate genes and flanking regulatory regions was performed on 61 patients treated with bevacizumab (19 cases developed early-onset grade 3 hypertension and 42 controls had no reported hypertension in the first six cycles of treatment). SNP-based tests for common variant associations and gene-based tests for rare variant associations were performed in 174 candidate genes.Results: Four common variants in independent linkage disequilibrium blocks between SLC29A1 and HSP90AB1 were among the top associations. Validation in larger bevacizumab-treated cohorts supported association between rs9381299 with early grade 3+ hypertension (P = 0.01; OR, 2.4) and systolic blood pressure >180 mm Hg (P = 0.02; OR, 2.1). rs834576 was associated with early grade 3+ hypertension in CALGB 40502 (P = 0.03; OR, 2.9). These SNP regions are enriched for regulatory elements that may potentially increase gene expression. In vitro overexpression of SLC29A1 in human endothelial cells disrupted adenosine signaling and reduced nitric oxide levels that were further lowered upon bevacizumab exposure.Conclusions: The genomic region between SLC29A1 and HSP90AB1 and its role in regulating adenosine signaling are key targets for further investigation into the pathogenesis of bevacizumab-induced hypertension. Clin Cancer Res; 24(19); 4734-44. ©2018 AACR.

Genetic polymorphism of SLC31A1 is associated with clinical outcomes of platinum-based chemotherapy in non-small-cell lung cancer patients through modulating microRNA-mediated regulation

SLC31A1 is the major transporter for platinum drug intake, its expression correlates with drug disposition and response. In 1004 Chinese NSCLC patients with platinum-based chemotherapy, we investigated the association between SLC31A1 polymorphisms and clinical outcomes. Heterozygotes of rs10759637 at 3′UTR was associated with severe thrombocytopenia (odds ratio [OR]: 2.69; P = 0.012) and shorter overall survival (hazard ratio [HR]: 1.24; P = 0.005). Variant homozygote of rs2233914 was correlated with longer overall survival (hazard ratio [HR]: 0.73; P = 0.008). Haplotype and diplotype of these linked SNPs were associated with hematologic toxicities. In stratification analyses, rs10759637 and rs2233914 consistently correlated with overall survival in specific subgroups such as men, smoker, patients older than 58 years, or with ECOG PS 0-1, or with squamous cell carcinoma. rs10759637 could change the local structure of 3′UTR harboring putative binding sites for hsa-miR-29, whose transfection into 16HBE cells resulted in remarkable suppression of gene expression. The rs10759637 variant significantly correlated with lowered luciferase activity in reporter assays and decreased expression of SLC31A1 transcript in tumorous tissues. The study thereby identified functional polymorphism of SLC31A1 that modulates miRNA-3′UTR interaction and gene expression as potential pharmacogenetic biomarker for clinical outcomes of platinum-based chemotherapy in NSCLC patients.

Use of exome sequencing to determine the full profile of genetic variants in the fluoropyrimidine pathway in colorectal cancer patients affected by severe toxicity

Aim: To identify genetic variants associated with capecitabine toxicity in fluoropyrimidine pathway genes using exome sequencing. Patients & methods: Exomes from eight capecitabine-treated patients with severe adverse reactions (grade >2), among a population of 319, were sequenced (Ion Proton). SNPs in genes classified as potentially damaging (Sorting Intolerant from Tolerant and Polymorphism Phenotyping v2) were tested for association with toxicity in a validation cohort of 319 capecitabine-treated patients. Results: A total of 17 nonsynonymous genetic variants were identified. Of these, five putative damaging SNPs in DPYD, ABCC4 and MTHFR were genotyped in the validation cohort. DPYD rs1801160 was associated with the risk of toxicity (p = 0.029) and MTHFR rs1801133 with delayed administration of chemotherapy due to toxicity (p = 0.047). Conclusion: Exome sequencing revealed two specific biomarkers of the risk of toxicity to capecitabine.

Chemotherapy-induced neutropenia as a prognostic factor in patients with metastatic pancreatic cancer treated with gemcitabine

PURPOSE

Chemotherapy-induced neutropenia (CIN) is a common side effect of chemotherapy and an important dose-limiting factor. However, an association between CIN development and longer survival was recently reported in several solid cancers. In the present study, we aimed to assess whether CIN could be a prognostic factor and clarify other prognostic factors for patients with metastatic pancreatic cancer.

METHODS

We retrospectively analyzed the medical records of 84 patients who received gemcitabine monotherapy as first-line chemotherapy for metastatic pancreatic cancer to assess whether CIN could be a prognostic factor. Potential prognostic factors of survival were examined by univariate and multivariate analyses using the log-rank test and Cox proportional hazard model, respectively.

RESULTS

Median survival time was 170 days [95% confidence interval (CI), 147-193] in patients without CIN (grade 0), 301 days (95% CI, 152-450) in patients with grade 1-2 CIN, and 406 days (95% CI, 271-541) in patients with grade 3 CIN. The multivariate analysis revealed that a pretreatment C-reactive protein level of <0.50 mg/dL [hazard ratio (HR), 0.534; 95% CI, 0.323-0.758, P = 0.015] and grade 3 CIN (HR, 0.447; 95% CI, 0.228-0.875, P = 0.019) were independent favorable prognostic factors in patients with metastatic pancreatic cancer treated with gemcitabine.

CONCLUSIONS

Neutropenia during chemotherapy was associated with increased survival of patients with metastatic pancreatic cancer. Monitoring of CIN could be used to predict treatment responsiveness.

Gemcitabine reduces MDSCs, tregs and TGFβ-1 while restoring the teff/treg ratio in patients with pancreatic cancer

Background

Cancer immunotherapy can be potentiated by conditioning regimens such as cyclophosphamide, which reduces the level of regulatory T cells (tregs). However, myeloid suppressive cells are still remaining. Accordingly to previous reports, gemcitabine improves immune status of cancer patients. In this study, the role of gemcitabine was further explored to map its immunological target cells and molecules in patients with pancreatic cancer.

Methods

Patient blood was investigated by flow cytometry and cytokine arrays at different time points during gemcitabine treatment.

Results

The patients had elevated myeloid-derived suppressor cells (MDSCs), and Tregs at diagnosis. Myeloid cells were in general decreased by gemcitabine. The granulocytic MDSCs were significantly reduced while monocytic MDSCs were not affected. In vitro, monocytes responding to IL-6 by STAT3 phosphorylation were prevented to respond in gemcitabine medium. However, gemcitabine could not prevent STAT3 phosphorylation in IL-6-treated tumor cell lines. TGFβ-1 was significantly reduced after only one treatment and continued to decrease. At the same time, the effector T cell:Treg ratio was increased and the effector T cells had full proliferative capacity during the gemcitabine cycle. However, after a resting period, the level of suppressor cells and TGFβ-1 had been restored showing the importance of continuous conditioning.

Conclusions

Gemcitabine regulates the immune system in patients with pancreatic cancer including MDSCs, Tregs and molecules such as TGFβ-1 but does not hamper the ability of effector lymphocytes to expand to stimuli. Hence, it may be of high interest to use gemcitabine as a conditioning strategy together with immunotherapy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-016-1037-z) contains supplementary material, which is available to authorized users.