Splicing factor ESRP1 controls ER-positive breast cancer by altering metabolic pathways.

The epithelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2) control the epithelial-to-mesenchymal transition (EMT) splicing program in cancer. However, their role in breast cancer recurrence is unclear. In this study, we report that high levels of ESRP1, but not ESRP2, are associated with poor prognosis in estrogen receptor positive (ER+) breast tumors. Knockdown of ESRP1 in endocrine-resistant breast cancer models decreases growth significantly and alters the EMT splicing signature, which we confirm using TCGA SpliceSeq data of ER+ BRCA tumors. However, these changes are not accompanied by the development of a mesenchymal phenotype or a change in key EMT-transcription factors. In tamoxifen-resistant cells, knockdown of ESRP1 affects lipid metabolism and oxidoreductase processes, resulting in the decreased expression of fatty acid synthase (FASN), stearoyl-CoA desaturase 1 (SCD1), and phosphoglycerate dehydrogenase (PHGDH) at both the mRNA and protein levels. Furthermore, ESRP1 knockdown increases the basal respiration and spare respiration capacity. This study reports a novel role for ESRP1 that could form the basis for the prevention of tamoxifen resistance in ER+ breast cancer.

Expression levels of SF3B3 correlate with prognosis and endocrine resistance in estrogen receptor-positive breast cancer.

De novo or acquired resistance to endocrine therapy limits its utility in a significant number of estrogen receptor-positive (ER-positive) breast cancers. It is crucial to identify novel targets for therapeutic intervention and improve the success of endocrine therapies. Splicing factor 3b, subunit 1 (SF3B1) mutations are described in luminal breast cancer albeit in low frequency. In this study, we evaluated the role of SF3B1 and SF3B3, critical parts of the SF3b splicing complex, in ER-positive endocrine resistance. To ascertain the role of SF3B1/SF3B3 in endocrine resistance, their expression levels were evaluated in ER-positive/endocrine-resistant cell lines (MCF-7/LCC2 and MCF-7/LCC9) using a real-time quantitative reverse transcription PCR (qRT-PCR). To further determine their clinical relevance, expression analysis was performed in a cohort of 60 paraffin-embedded ER-positive, node-negative breast carcinomas with low, intermediate, and high Oncotype DX recurrence scores. Expression levels of SF3B1 and SF3B3 and their prognostic value were validated in large cohorts using publicly available gene expression data sets including The Cancer Genome Atlas. SF3B1 and SF3B3 levels were significantly increased in ERα-positive cells with acquired tamoxifen (MCF-7/LCC2; both P<0.0002) and fulvestrant/tamoxifen resistance (MCF-7/LCC9; P=0.008 for SF3B1 and P=0.0006 for SF3B3). Expression levels of both MCF-7/LCC2 and MCF-7/LCC9 were not affected by additional treatments with E2 and/or tamoxifen. Furthermore, qRT-PCR analysis confirmed that SF3B3 expression is significantly upregulated in Oncotype DX high-risk groups when compared with low risk (P=0.019). Similarly, in publicly available breast cancer gene expression data sets, overexpression of SF3B3, but not SF3B1, was significantly correlated with overall survival. Furthermore, the correlation was significant in ER-positive, but not in ER-negative tumors.This is the first study to document the role of SF3B3 in endocrine resistance and prognosis in ER-positive breast cancer. Potential strategies for therapeutic targeting of the splicing mechanism(s) need to be evaluated.

Establishment and characterization of a novel cell line derived from human thymoma AB tumor.

Thymomas are low-grade epithelial tumors of the anterior mediastinum. The complexity of the disease and the lack of in vitro and in vivo models hamper the development of better therapeutics. In this study, we report a novel cell line, designated as IU-TAB-1, which was established from a patient with stage II thymoma (World Health Organization-type AB). The IU-TAB-1 cell line was established in vitro and characterized using histological and immunohistochemical staining, fluorescence-activated cell sorting, cytogenetic analyses and functional assays including in vitro and a NOD/SCID xenograft model. A whole-genome gene expression analysis (Illumina) was performed on the IU-TAB-1 cell line and 34 thymomas to determine the clinical relevance of the cell line. The IU-TAB-1 cell line was positive for epithelial markers (pan-cytokeratin and EpCAM/CD326) including thymic epithelial (TE) surface markers (such as CD29, CD9, CD54/ICAM-1, CD58 and CD24) and p63, and negative for B- and T-cell lineage markers. Gene expression profiling demonstrated overlapping and distinct genes between IU-TAB-1 and primary thymomas including the primary tumor (from which the cell line was derived). IU-TAB-1 cells are tumorigenic when implanted in immunodeficient mice with tumors reaching a volume of 1000 mm³ at around 130 days. The established cell line represents a biologically relevant new tool to investigate the molecular pathology of thymic malignancies and to evaluate the efficacy of novel therapeutics both in vitro and in vivo.

Persistent upregulation of U6:SNORD44 small RNA ratio in the serum of breast cancer patients.

INTRODUCTION: Serum microRNAs have the potential to be valuable biomarkers of cancer. This investigation addresses two issues that impact their utility: a) appropriate normalization controls and b) whether their altered levels persist in patients who are clinically free of the disease. METHODS: Sera from 40 age-matched healthy women and 39 breast cancer patients without clinical disease at the time of serum collection were analyzed for microRNAs let-7f, miR-16, miR-21 and miR-155 using quantitative real-time PCR. U6 and 5S, which are transcribed by RNA polymerase III (RNAP-III) and the small nucleolar RNU44 (SNORD44), were also analyzed for normalization. Significant results from the initial study were verified using a second set of sera from 15 healthy patients, 15 breast cancer patients without clinical disease and 15 with metastatic disease, and a third set of 12 healthy and 18 patients with metastatic disease. U6 was further verified in the extended second cohort of 75 healthy and 68 breast cancer patients without clinical disease. RESULTS: U6:SNORD44 ratio was consistently higher in breast cancer patients with or without active disease (fold change range 1.5-6.6, p value range 0.0003 to 0.05). This increase in U6:SNORD44 ratio was observed in the sera of both estrogen receptor-positive (ER+) and ER-negative breast cancer patients. MiR-16 and 5S, which are often used as normalization controls for microRNAs, showed remarkable experimental variability and thus are not ideal for normalization. CONCLUSIONS: Elevated serum U6 levels in breast cancer patients irrespective of disease activity at the time of serum collection suggest a new paradigm in cancer; persistent systemic changes during cancer progression, which result in elevated activity of RNAP-III and/or the stability/release pathways of U6 in non-cancer tissues. Additionally, these results highlight the need for developing standards for normalization between samples in microRNA-related studies for healthy versus cancer and for inter-laboratory reproducibility. Our studies rule out the utility of miR-16, U6 and 5S RNAs for this purpose.

Prostate cancer sheds the αvß3 integrin in vivo through exosomes.

The αvß3 integrin has been shown to promote aggressive phenotypes in many types of cancers, including prostate cancer. We show that GFP-labeled αvß3 derived from cancer cells circulates in the blood and is detected in distant lesions in NOD scid gamma (NSG) mice. We, therefore, hypothesized that αvß3 travels through exosomes and tested its levels in pools of vesicles, which we designate extracellular vesicles highly enriched in exosomes (ExVs), and in exosomes isolated from the plasma of prostate cancer patients. Here, we show that the αvß3 integrin is found in patient blood exosomes purified by sucrose or iodixanol density gradients. In addition, we provide evidence that the αvß3 integrin is transferred through ExVs isolated from prostate cancer patient plasma to ß3-negative recipient cells. We also demonstrate the intracellular localization of ß3-GFP transferred via cancer cell-derived ExVs. We show that the ExVs present in plasma from prostate cancer patients contain higher levels of αvß3 and CD9 as compared to plasma ExVs from age-matched subjects who are not affected by cancer. Furthermore, using PSMA antibody-bead mediated immunocapture, we show that the αvß3 integrin is expressed in a subset of exosomes characterized by PSMA, CD9, CD63, and an epithelial-specific marker, Trop-2. Finally, we present evidence that the levels of αvß3, CD63, and CD9 remain unaltered in ExVs isolated from the blood of prostate cancer patients treated with enzalutamide. Our results suggest that detecting exosomal αvß3 integrin in prostate cancer patients could be a clinically useful and non-invasive biomarker to follow prostate cancer progression. Moreover, the ability of αvß3 integrin to be transferred from ExVs to recipient cells provides a strong rationale for further investigating the role of αvß3 integrin in the pathogenesis of prostate cancer and as a potential therapeutic target.

Gene expression profiles among murine strains segregate with distinct differences in the progression of radiation-induced lung disease.

Molecular mechanisms underlying development of acute pneumonitis and/or late fibrosis following thoracic irradiation remain poorly understood. Here, we hypothesize that heterogeneity in disease progression and phenotypic expression of radiation-induced lung disease (RILD) across murine strains presents an opportunity to better elucidate mechanisms driving tissue response toward pneumonitis and/or fibrosis. Distinct differences in disease progression were observed in age- and sex-matched CBA/J, C57L/J and C57BL/6J mice over 1 year after graded doses of whole-thorax lung irradiation (WTLI). Separately, comparison of gene expression profiles in lung tissue 24 h post-exposure demonstrated >5000 genes to be differentially expressed (P<0.01; >twofold change) between strains with early versus late onset of disease. An immediate divergence in early tissue response between radiation-sensitive and -resistant strains was observed. In pneumonitis-prone C57L/J mice, differentially expressed genes were enriched in proinflammatory pathways, whereas in fibrosis-prone C57BL/6J mice, genes were enriched in pathways involved in purine and pyrimidine synthesis, DNA replication and cell division. At 24 h post-WTLI, different patterns of cellular damage were observed at the ultrastructural level among strains but microscopic damage was not yet evident under light microscopy. These data point toward a fundamental difference in patterns of early pulmonary tissue response to WTLI, consistent with the macroscopic expression of injury manifesting weeks to months after exposure. Understanding the mechanisms underlying development of RILD might lead to more rational selection of therapeutic interventions to mitigate healthy tissue damage.

Molecular Insights of Pathways Resulting from Two Common PIK3CA Mutations in Breast Cancer.

The PI3K pathway is activated in approximately 70% of breast cancers. PIK3CA gene mutations or amplifications that affect the PI3K p110α subunit account for activation of this pathway in 20% to 40% of cases, particularly in estrogen receptor alpha (ERα)-positive breast cancers. AKT family of kinases, AKT1-3, are the downstream targets of PI3K and these kinases activate ERα. Although several inhibitors of PI3K have been developed, none has proven effective in the clinic, partly due to an incomplete understanding of the selective routing of PI3K signaling to specific AKT isoforms. Accordingly, we investigated in this study the contribution of specific AKT isoforms in connecting PI3K activation to ERα signaling, and we also assessed the utility of using the components of PI3K-AKT isoform-ERα signaling axis as predictive biomarkers of response to PI3K inhibitors. Using a variety of physiologically relevant model systems with defined natural or knock-in PIK3CA mutations and/or PI3K hyperactivation, we show that PIK3CA-E545K mutations (found in ∼20% of PIK3CA-mutant breast cancers), but not PIK3CA-H1047R mutations (found in 55% of PIK3CA-mutant breast cancers), preferentially activate AKT1. Our findings argue that AKT1 signaling is needed to respond to estrogen and PI3K inhibitors in breast cancer cells with PIK3CA-E545K mutation, but not in breast cancer cells with other PIK3CA mutations. This study offers evidence that personalizing treatment of ER-positive breast cancers to PI3K inhibitor therapy may benefit from an analysis of PIK3CA-E545K-AKT1-estrogen signaling pathways. Cancer Res; 76(13); 3989-4001. ©2016 AACR.

HLA-genotyping of clinical specimens using Ion Torrent-based NGS.

We have evaluated and validated the NXType™ workflow (One Lambda, Inc.) and the accompanying TypeStream™ software on the Ion Torrent Next Generation Sequencing (NGS) platform using a comprehensive testing panel. The panel consisted of 285 genomic DNA (gDNA) samples derived from four major ethnic populations and contained 59 PT samples and 226 clinical specimens. The total number of alleles from the six loci interrogated by NGS was 3420. This validation panel provided a wide range of HLA sequence variations including many rare alleles, new variants and homozygous alleles. The NXType™ system (reagents and software) was able to correctly genotype the vast majority of these specimens. The concordance rate between SBT-derived genotypes and those generated by TypeStream™ auto-analysis ranged from 99.5% to 99.8% for the HLA-A, B, C, DRB1 and DQB1 loci, and was 98.9% for HLA-DPB1. A strategy for data review was developed that would allow correction of most of the few remaining typing errors. The entire NGS workflow from gDNA amplification to genotype assignment could be completed within 3 working days. Through this validation study, the limitations and shortcomings of the platform, specific assay system, and software algorithm were also revealed for further evaluation and improvement.

Validation of a next-generation-sequencing cancer panel for use in the clinical laboratory.

CONTEXT: Next-generation sequencing allows for high-throughput processing and sensitive variant detection in multiple genes from small samples. For many diseases, including cancer, a comprehensive mutational profile of a targeted list of genes can be used to simultaneously inform patient care, establish eligibility for ongoing clinical trials, and further research. OBJECTIVE: To validate a pan-cancer, next-generation-sequencing assay for use in the clinical laboratory. DESIGN: DNA was extracted from 68 clinical specimens (formalin-fixed, paraffin-embedded; fine-needle aspirates; peripheral blood; or bone marrow) and 5 normal controls. Sixty-four DNA samples (94%; 64 of 68) were successfully processed with the TruSeq Amplicon Cancer Panel (Illumina Inc, San Diego, California) and sequenced in 4 sequencing runs. The data were analyzed at 4 different filter settings for sequencing coverage and variant frequency cutoff. RESULTS: Libraries created from 40 specimens could be successfully sequenced in a single run and still yield sufficient coverage for robust data analysis of individual samples. Sensitivity for mutation detection down to 5% was demonstrated using dilutions of clinical specimens and control samples. The test was highly repeatable and reproducible and showed 100% concordance with clinically validated Sanger sequencing results. Comparison to an alternate next-generation sequencing technology was performed by also processing 9 of the specimens with the AmpliSeq Cancer Hotspot Panel (version 2; Life Technologies, Grand Island, New York). Thirty of the 31 (97%) TruSeq-detected variants covered by the designs of both panels were confirmed. CONCLUSIONS: A sensitive, high-throughput, pan-cancer mutation panel for sequencing of cancer hot-spot mutations in 42 genes was validated for routine use in clinical testing.

Organ-specific adaptive signaling pathway activation in metastatic breast cancer cells.

Breast cancer metastasizes to bone, visceral organs, and/or brain depending on the subtype, which may involve activation of a host organ-specific signaling network in metastatic cells. To test this possibility, we determined gene expression patterns in MDA-MB-231 cells and its mammary fat pad tumor (TMD-231), lung-metastasis (LMD-231), bone-metastasis (BMD-231), adrenal-metastasis (ADMD-231) and brain-metastasis (231-BR) variants. When gene expression between metastases was compared, 231-BR cells showed the highest gene expression difference followed by ADMD-231, LMD-231, and BMD-231 cells. Neuronal transmembrane proteins SLITRK2, TMEM47, and LYPD1 were specifically overexpressed in 231-BR cells. Pathway-analyses revealed activation of signaling networks that would enable cancer cells to adapt to organs of metastasis such as drug detoxification/oxidative stress response/semaphorin neuronal pathway in 231-BR, Notch/orphan nuclear receptor signals involved in steroidogenesis in ADMD-231, acute phase response in LMD-231, and cytokine/hematopoietic stem cell signaling in BMD-231 cells. Only NF-κB signaling pathway activation was common to all except BMD-231 cells. We confirmed NF-κB activation in 231-BR and in a brain metastatic variant of 4T1 cells (4T1-BR). Dimethylaminoparthenolide inhibited NF-κB activity, LYPD1 expression, and proliferation of 231-BR and 4T1-BR cells. Thus, transcriptome change enabling adaptation to host organs is likely one of the mechanisms associated with organ-specific metastasis and could potentially be targeted therapeutically.

Cancer affects microRNA expression, release, and function in cardiac and skeletal muscle.

Circulating microRNAs (miRNA) are emerging as important biomarkers of various diseases, including cancer. Intriguingly, circulating levels of several miRNAs are lower in patients with cancer compared with healthy individuals. In this study, we tested the hypothesis that a circulating miRNA might serve as a surrogate of the effects of cancer on miRNA expression or release in distant organs. Here we report that circulating levels of the muscle-enriched miR486 is lower in patients with breast cancer compared with healthy individuals and that this difference is replicated faithfully in MMTV-PyMT and MMTV-Her2 transgenic mouse models of breast cancer. In tumor-bearing mice, levels of miR486 were relatively reduced in muscle, where there was elevated expression of the miR486 target genes PTEN and FOXO1A and dampened signaling through the PI3K/AKT pathway. Skeletal muscle expressed lower levels of the transcription factor MyoD, which controls miR486 expression. Conditioned media (CM) obtained from MMTV-PyMT and MMTV-Her2/Neu tumor cells cultured in vitro were sufficient to elicit reduced levels of miR486 and increased PTEN and FOXO1A expression in C2C12 murine myoblasts. Cytokine analysis implicated tumor necrosis factor α (TNFα) and four additional cytokines as mediators of miR486 expression in CM-treated cells. Because miR486 is a potent modulator of PI3K/AKT signaling and the muscle-enriched transcription factor network in cardiac/skeletal muscle, our findings implicated TNFα-dependent miRNA circuitry in muscle differentiation and survival pathways in cancer.

Identification and validation of genes with expression patterns inverse to multiple metastasis suppressor genes in breast cancer cell lines.

Metastasis suppressor genes (MSGs) have contributed to an understanding of regulatory pathways unique to the lethal metastatic process. When re-expressed in experimental models, MSGs block cancer spread to, and colonization of distant sites without affecting primary tumor formation. Genes have been identified with expression patterns inverse to a single MSG, and found to encode functional, druggable signaling pathways. We now hypothesize that common signaling pathways mediate the effects of multiple MSGs. By gene expression profiling of human MCF7 breast carcinoma cells expressing a scrambled siRNA, or siRNAs to each of 19 validated MSGs (NME1, BRMS1, CD82, CDH1, CDH2, CDH11, CASP8, MAP2K4, MAP2K6, MAP2K7, MAPK14, GSN, ARHGDIB, AKAP12, DRG1, CD44, PEBP1, RRM1, KISS1), we identified genes whose expression was significantly opposite to at least five MSGs. Five genes were selected for further analysis: PDE5A, UGT1A, IL11RA, DNM3 and OAS1. After stable downregulation of each candidate gene in the aggressive human breast cancer cell line MDA-MB-231T, in vitro motility was significantly inhibited. Two stable clones downregulating PDE5A (phosphodiesterase 5A), an enzyme involved in the regulation of cGMP-specific signaling, exhibited no difference in cell proliferation, but reduced motility by 47 and 66 % compared to the empty vector-expressing cells (p = 0.01 and p = 0.005). In an experimental metastasis assay, two shPDE5A-MDA-MB-231T clones produced 47-62 % fewer lung metastases than shRNA-scramble expressing cells (p = 0.045 and p = 0.009 respectively). This study demonstrates that previously unrecognized genes are inversely related to the expression of multiple MSGs, contribute to aspects of metastasis, and may stand as novel therapeutic targets.

Gene Expression Analysis Reveals Distinct Pathways of Resistance to Bevacizumab in Xenograft Models of Human ER-Positive Breast Cancer.

Bevacizumab, the recombinant antibody targeting vascular endothelial growth factor (VEGF), improves progression-free but not overall survival in metastatic breast cancer. To seek further insights in resistance mechanisms to bevacizumab at the molecular level, we developed VEGF and non-VEGF-driven ER-positive MCF7-derived xenograft models allowing comparison of tumor response at different timepoints. VEGF gene (MV165) overexpressing xenografts were initially sensitive to bevacizumab, but eventually acquired resistance. In contrast, parental MCF7 cells derived tumors were de novo insensitive to bevacizumab. Microarray analysis with qRT-PCR validation revealed that Follistatin (FST) and NOTCH were the top signaling pathways associated with resistance in VEGF-driven tumors (P<0.05). Based on the presence of VEGF, treatment with bevacizumab resulted in altered patterns of metagenes and PAM50 gene expression. In VEGF-driven model after short and long-term bevacizumab treatments, a change in the intrinsic subtype (luminal to myoepithelial/basal-like) was observed in association with increased expression of genes implicated with cancer stem cell phenotype (P<0.05). Our results show that the presence or absence of VEGF expression affects the response to bevacizumab therapy and gene pathways. In particular, long-term bevacizumab treatment shifts the cancer cells to a more aggressive myoepithelial/basal subtype in VEGF-expressing model, but not in non-VEGF model. These findings could shed light on variable results to anti-VEGF therapy in patients and emphasize the importance of patient stratification based on the VEGF expression. Our data strongly suggest consideration of patient subgroups for treatment and designing novel combinatory therapies in the clinical setting.

Identification of FDA-approved drugs targeting breast cancer stem cells along with biomarkers of sensitivity.

Recently developed genomics-based tools are allowing repositioning of Food and Drug Administration (FDA)-approved drugs as cancer treatments, which were employed to identify drugs that target cancer stem cells (CSCs) of breast cancer. Gene expression datasets of CSCs from six studies were subjected to connectivity map to identify drugs that may ameliorate gene expression patterns unique to CSCs. All-trans retinoic acid (ATRA) was negatively connected with gene expression in CSCs. ATRA reduced mammosphere-forming ability of a subset of breast cancer cells, which correlated with induction of apoptosis, reduced expression of SOX2 but elevated expression of its antagonist CDX2. SOX2/CDX2 ratio had prognostic relevance in CSC-enriched breast cancers. K-ras mutant breast cancer cell line enriched for CSCs was resistant to ATRA, which was reversed by MAP kinase inhibitors. Thus, ATRA alone or in combination can be tested for efficacy using SOX2, CDX2, and K-ras mutation/MAPK activation status as biomarkers of response.

A molecular signature of normal breast epithelial and stromal cells from Li-Fraumeni syndrome mutation carriers.

Specific changes in gene expression during cancer initiation should enable discovery of biomarkers for risk assessment, early detection and targets for chemoprevention. It has been previously demonstrated that altered mRNA and proteome signatures of morphologically normal cells bearing a single inherited "hit" in a tumor suppressor gene parallel many changes observed in the corresponding sporadic cancer. Here, we report on the global gene expression profile of morphologically normal, cultured primary breast epithelial and stromal cells from Li-Fraumeni syndrome (LFS) TP53 mutation carriers. Our analyses identified multiple changes in gene expression in both morphologically normal breast epithelial and stromal cells associated with TP53 haploinsufficiency, as well as interlocking pathways. Notably, a dysregulated p53 signaling pathway was readily detectable. Pharmacological intervention with the p53 rescue compounds CP-31398 and PRIMA-1 provided further evidence in support of the central role of p53 in affecting these changes in LFS cells and treatment for this cancer. Because loss of signaling mediated by TP53 is associated with the development and survival of many human tumors, identification of gene expression profiles in morphologically normal cells that carry "one-hit" p53 mutations may reveal novel biomarkers, enabling the discovery of potential targets for chemoprevention of sporadic tumors as well.