CDC25B partners with PP2A to induce AMPK activation and tumor suppression in triple negative breast cancer

Cell division cycle 25 (CDC25) dual specificity phosphatases positively regulate the cell cycle by activating cyclin-dependent kinase/cyclin complexes. Here, we demonstrate that in addition to its role in cell cycle regulation, CDC25B functions as a regulator of protein phosphatase 2A (PP2A), a major cellular Ser/Thr phosphatase, through its direct interaction with PP2A catalytic subunit. Importantly, CDC25B alters the regulation of AMP-activated protein kinase signaling (AMPK) by PP2A, increasing AMPK activity by inhibiting PP2A to dephosphorylate AMPK. CDC25B depletion leads to metformin resistance by inhibiting metformin-induced AMPK activation. Furthermore, dual inhibition of CDC25B and PP2A further inhibits growth of 3D organoids isolated from patient derived xenograft model of breast cancer compared to CDC25B inhibition alone. Our study identifies CDC25B as a regulator of PP2A, and uncovers a mechanism of controlling the activity of a key energy metabolism marker, AMPK

Deep Sequence Analysis of Non-Small Cell Lung Cancer: Integrated Analysis of Gene Expression, Alternative Splicing, and Single Nucleotide Variations in Lung Adenocarcinomas with and without Oncogenic KRAS Mutations

KRAS mutations are highly prevalent in non-small cell lung cancer (NSCLC), and tumors harboring these mutations tend to be aggressive and resistant to chemotherapy. We used next-generation sequencing technology to identify pathways that are specifically altered in lung tumors harboring a KRAS mutation. Paired-end RNA-sequencing of 15 primary lung adenocarcinoma tumors (8 harboring mutant KRAS and 7 with wild-type KRAS) were performed. Sequences were mapped to the human genome, and genomic features, including differentially expressed genes, alternate splicing isoforms and single nucleotide variants, were determined for tumors with and without KRAS mutation using a variety of computational methods. Network analysis was carried out on genes showing differential expression (374 genes), alternate splicing (259 genes), and SNV-related changes (65 genes) in NSCLC tumors harboring a KRAS mutation. Genes exhibiting two or more connections from the lung adenocarcinoma network were used to carry out integrated pathway analysis. The most significant signaling pathways identified through this analysis were the NFκB, ERK1/2, and AKT pathways. A 27 gene mutant KRAS-specific sub network was extracted based on gene–gene connections from the integrated network, and interrogated for druggable targets. Our results confirm previous evidence that mutant KRAS tumors exhibit activated NFκB, ERK1/2, and AKT pathways and may be preferentially sensitive to target therapeutics toward these pathways. In addition, our analysis indicates novel, previously unappreciated links between mutant KRAS and the TNFR and PPARγ signaling pathways, suggesting that targeted PPARγ antagonists and TNFR inhibitors may be useful therapeutic strategies for treatment of mutant KRAS lung tumors. Our study is the first to integrate genomic features from RNA-Seq data from NSCLC and to define a first draft genomic landscape model that is unique to tumors with oncogenic KRAS mutations

Pathway-Based Analysis of Genome-Wide Association Data Identified SNPs in HMMR as Biomarker for Chemotherapy-Induced Neutropenia in Breast Cancer Patients

Neutropenia secondary to chemotherapy in breast cancer patients can be life-threatening and there are no biomarkers available to predict the risk of drug-induced neutropenia in those patients. We previously performed a genome-wide association study (GWAS) for neutropenia events in women with breast cancer who were treated with 5-fluorouracil, epirubicin and cyclophosphamide and recruited to the SUCCESS A trial. A genome-wide significant single-nucleotide polymorphism (SNP) signal in the tumor necrosis factor superfamily member 13B (TNFSF13B) gene, encoding the cytokine B-cell activating factor (BAFF), was identified in that GWAS. Taking advantage of these existing GWAS data, in the present study we utilized a pathway-based analysis approach by leveraging knowledge of the pharmacokinetics and pharmacodynamics of drugs and breast cancer pathophysiology to identify additional SNPs/genes associated with the underlying etiology of chemotherapy-induced neutropenia. We identified three SNPs in the hyaluronan mediated motility receptor (HMMR) gene that were significantly associated with neutropenia (p < 1.0E-04). Those three SNPs were trans-expression quantitative trait loci for the expression of TNFSF13B (p < 1.0E-04). The minor allele of these HMMR SNPs was associated with a decreased TNFSF13B mRNA level. Additional functional studies performed with lymphoblastoid cell lines (LCLs) demonstrated that LCLs possessing the minor allele for the HMMR SNPs were more sensitive to drug treatment. Knock-down of TNFSF13B in LCLs and HL-60 promyelocytic cells and treatment of those cells with BAFF modulated the cell sensitivity to chemotherapy treatment. These results demonstrate that HMMR SNP-dependent cytotoxicity of these chemotherapeutic agents might be related to TNFSF13B expression level. In summary, utilizing a pathway-based approach for the analysis of GWAS data, we identified additional SNPs in the HMMR gene that were associated with neutropenia and also were correlated with TNFSF13B expression

TREAT: a bioinformatics tool for variant annotations and visualizations in targeted and exome sequencing data

Summary: TREAT (Targeted RE-sequencing Annotation Tool) is a tool for facile navigation and mining of the variants from both targeted resequencing and whole exome sequencing. It provides a rich integration of publicly available as well as in-house developed annotations and visualizations for variants, variant-hosting genes and host-gene pathways

A novel bioinformatics pipeline for identification and characterization of fusion transcripts in breast cancer and normal cell lines

SnowShoes-FTD, developed for fusion transcript detection in paired-end mRNA-Seq data, employs multiple steps of false positive filtering to nominate fusion transcripts with near 100% confidence. Unique features include: (i) identification of multiple fusion isoforms from two gene partners; (ii) prediction of genomic rearrangements; (iii) identification of exon fusion boundaries; (iv) generation of a 5′–3′ fusion spanning sequence for PCR validation; and (v) prediction of the protein sequences, including frame shift and amino acid insertions. We applied SnowShoes-FTD to identify 50 fusion candidates in 22 breast cancer and 9 non-transformed cell lines. Five additional fusion candidates with two isoforms were confirmed. In all, 30 of 55 fusion candidates had in-frame protein products. No fusion transcripts were detected in non-transformed cells. Consideration of the possible functions of a subset of predicted fusion proteins suggests several potentially important functions in transformation, including a possible new mechanism for overexpression of ERBB2 in a HER-positive cell line. The source code of SnowShoes-FTD is provided in two formats: one configured to run on the Sun Grid Engine for parallelization, and the other formatted to run on a single LINUX node. Executables in PERL are available for download from our web site: http://mayoresearch.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm

Integrated Analysis of Gene Expression, CpG Island Methylation, and Gene Copy Number in Breast Cancer Cells by Deep Sequencing

We used deep sequencing technology to profile the transcriptome, gene copy number, and CpG island methylation status simultaneously in eight commonly used breast cell lines to develop a model for how these genomic features are integrated in estrogen receptor positive (ER+) and negative breast cancer. Total mRNA sequence, gene copy number, and genomic CpG island methylation were carried out using the Illumina Genome Analyzer. Sequences were mapped to the human genome to obtain digitized gene expression data, DNA copy number in reference to the non-tumor cell line (MCF10A), and methylation status of 21,570 CpG islands to identify differentially expressed genes that were correlated with methylation or copy number changes. These were evaluated in a dataset from 129 primary breast tumors. Gene expression in cell lines was dominated by ER-associated genes. ER+ and ER− cell lines formed two distinct, stable clusters, and 1,873 genes were differentially expressed in the two groups. Part of chromosome 8 was deleted in all ER− cells and part of chromosome 17 amplified in all ER+ cells. These loci encoded 30 genes that were overexpressed in ER+ cells; 9 of these genes were overexpressed in ER+ tumors. We identified 149 differentially expressed genes that exhibited differential methylation of one or more CpG islands within 5 kb of the 5′ end of the gene and for which mRNA abundance was inversely correlated with CpG island methylation status. In primary tumors we identified 84 genes that appear to be robust components of the methylation signature that we identified in ER+ cell lines. Our analyses reveal a global pattern of differential CpG island methylation that contributes to the transcriptome landscape of ER+ and ER− breast cancer cells and tumors. The role of gene amplification/deletion appears to more modest, although several potentially significant genes appear to be regulated by copy number aberrations

Loss of Heterozygosity at the CYP2D6 Locus in Breast Cancer: Implications for Germline Pharmacogenetic Studies

Controversy exists regarding the impact of CYP2D6 genotype on tamoxifen responsiveness. We examined loss of heterozygosity (LOH) at the CYP2D6 locus and determined its impact on genotyping error when tumor tissue is used as a DNA source

Anastrozole has an association between degree of estrogen suppression and outcomes in early breast cancer and is a ligand for estrogen receptor α

Purpose: To determine if the degree of estrogen suppression with aromatase inhibitors (AI: anastrozole, exemestane, letrozole) is associated with efficacy in early-stage breast cancer, and to examine for differences in the mechanism of action between the three AIs. Experimental design: Matched case-control studies [247 matched sets from MA.27 (anastrozole vs. exemestane) and PreFace (letrozole) trials] were undertaken to assess whether estrone (E1) or estradiol (E2) concentrations after 6 months of adjuvant therapy were associated with risk of an early breast cancer event (EBCE). Preclinical laboratory studies included luciferase activity, cell proliferation, radio-labeled ligand estrogen receptor binding, surface plasmon resonance ligand receptor binding, and nuclear magnetic resonance assays. Results: Women with E1 ≥1.3 pg/mL and E2 ≥0.5 pg/mL after 6 months of AI treatment had a 2.2-fold increase in risk (P = 0.0005) of an EBCE, and in the anastrozole subgroup, the increase in risk of an EBCE was 3.0-fold (P = 0.001). Preclinical laboratory studies examined mechanisms of action in addition to aromatase inhibition and showed that only anastrozole could directly bind to estrogen receptor α (ERα), activate estrogen response element-dependent transcription, and stimulate growth of an aromatase-deficient CYP19A1-/- T47D breast cancer cell line. Conclusions: This matched case-control clinical study revealed that levels of estrone and estradiol above identified thresholds after 6 months of adjuvant anastrozole treatment were associated with increased risk of an EBCE. Preclinical laboratory studies revealed that anastrozole, but not exemestane or letrozole, is a ligand for ERα. These findings represent potential steps towards individualized anastrozole therapy