Chad J. Creighton

Genomics data from The Cancer Genome Atlas (TCGA) project has led to the comprehensive molecular characterization of multiple cancer types. The large sample numbers in TCGA offer an excellent opportunity to address questions associated with tumo heterogeneity. Exploration of the data by cancer researchers and clinicians is imperative to unearth novel therapeutic/diagnostic biomarkers. Various computational tools have been developed to aid researchers in carrying out specific TCGA data analyses; however there is need for resources to facilitate the study of gene expression variations and survival associations across tumors. Here, we report UALCAN, an easy to use, interactive web-portal to perform to in-depth analyses of TCGA gene expression data. UALCAN uses TCGA level 3 RNA-seq and clinical data from 31 cancer types. The portal's user-friendly features allow to perform: 1) analyze relative expression of a query gene(s) across tumor and normal samples, as well as in various tumor sub-groups based on individual cancer stages, tumor grade, race, body weight or other clinicopathologic features, 2) estimate the effect of gene expression level and clinicopathologic features on patient survival; and 3) identify the top over- and under-expressed (up and down-regulated) genes in individual cancer types. This resource serves as a platform for in silico validation of target genes and for identifying tumor sub-group specific candidate biomarkers. Thus, UALCAN web-portal could be extremely helpful in accelerating cancer research. UALCAN is publicly available at http://ualcan.path.uab.edu.

Liver cancer has the second highest worldwide cancer mortality rate and has limited therapeutic options. We analyzed 363 hepatocellular carcinoma (HCC) cases by whole-exome sequencing and DNA copy number analyses, and we analyzed 196 HCC cases by DNA methylation, RNA, miRNA, and proteomic expression also. DNA sequencing and mutation analysis identified significantly mutated genes, including LZTR1, EEF1A1, SF3B1, and SMARCA4. Significant alterations by mutation or downregulation by hypermethylation in genes likely to result in HCC metabolic reprogramming (ALB, APOB, and CPS1) were observed. Integrative molecular HCC subtyping incorporating unsupervised clustering of five data platforms identified three subtypes, one of which was associated with poorer prognosis in three HCC cohorts. Integrated analyses enabled development of a p53 target gene expression signature correlating with poor survival. Potential therapeutic targets for which inhibitors exist include WNT signaling, MDM4, MET, VEGFA, MCL1, IDH1, TERT, and immune checkpoint proteins CTLA-4, PD-1, and PD-L1.

Some breast cancers have been shown to contain a small fraction of cells characterized by CD44(+)/CD24(-/low) cell-surface antigen profile that have high tumor-initiating potential. In addition, breast cancer cells propagated in vitro as mammospheres (MSs) have also been shown to be enriched for cells capable of self-renewal. In this study, we have defined a gene expression signature common to both CD44(+)/CD24(-/low) and MS-forming cells. To examine its clinical significance, we determined whether tumor cells surviving after conventional treatments were enriched for cells bearing this CD44(+)/CD24(-/low)-MS signature. The CD44(+)/CD24(-/low)-MS signature was found mainly in human breast tumors of the recently identified "claudin-low" molecular subtype, which is characterized by expression of many epithelial-mesenchymal-transition (EMT)-associated genes. Both CD44(+)/CD24(-/low)-MS and claudin-low signatures were more pronounced in tumor tissue remaining after either endocrine therapy (letrozole) or chemotherapy (docetaxel), consistent with the selective survival of tumor-initiating cells posttreatment. We confirmed an increased expression of mesenchymal markers, including vimentin (VIM) in cytokeratin-positive epithelial cells metalloproteinase 2 (MMP2), in two separate sets of postletrozole vs. pretreatment specimens. Taken together, these data provide supporting evidence that the residual breast tumor cell populations surviving after conventional treatment may be enriched for subpopulations of cells with both tumor-initiating and mesenchymal features. Targeting proteins involved in EMT may provide a therapeutic strategy for eliminating surviving cells to prevent recurrence and improve long-term survival in breast cancer patients.

Oncogenic mutations in the small GTPase Ras are highly prevalent in cancer, but an understanding of the vulnerabilities of these cancers is lacking. We undertook a genome-wide RNAi screen to identify synthetic lethal interactions with the KRAS oncogene. We discovered a diverse set of proteins whose depletion selectively impaired the viability of Ras mutant cells. Among these we observed a strong enrichment for genes with mitotic functions. We describe a pathway involving the mitotic kinase PLK1, the anaphase-promoting complex/cyclosome, and the proteasome that, when inhibited, results in prometaphase accumulation and the subsequent death of Ras mutant cells. Gene expression analysis indicates that reduced expression of genes in this pathway correlates with increased survival of patients bearing tumors with a Ras transcriptional signature. Our results suggest a previously underappreciated role for Ras in mitotic progression and demonstrate a pharmacologically tractable pathway for the potential treatment of cancers harboring Ras mutations.

Cancer genomic, transcriptomic, and proteomic profiling has generated extensive data that necessitate the development of tools for its analysis and dissemination. We developed UALCAN to provide a portal for easy exploring, analyzing, and visualizing these data, allowing users to integrate the data to better understand the gene, proteins, and pathways perturbed in cancer and make discoveries. UALCAN web portal enables analyzing and delivering cancer transcriptome, proteomics, and patient survival data to the cancer research community. With data obtained from The Cancer Genome Atlas (TCGA) project, UALCAN has enabled users to evaluate protein-coding gene expression and its impact on patient survival across 33 types of cancers. The web portal has been used extensively since its release and received immense popularity, underlined by its usage from cancer researchers in more than 100 countries. The present manuscript highlights the task we have undertaken and updates that we have made to UALCAN since its release in 2017. Extensive user feedback motivated us to expand the resource by including data on a) microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and promoter DNA methylation from TCGA and b) mass spectrometry-based proteomics from the Clinical Proteomic Tumor Analysis Consortium (CPTAC). UALCAN provides easy access to pre-computed, tumor subgroup-based gene/protein expression, promoter DNA methylation status, and Kaplan-Meier survival analyses. It also provides new visualization features to comprehend and integrate observations and aids in generating hypotheses for testing. UALCAN is accessible at http://ualcan.path.uab.edu.

We describe the landscape of somatic genomic alterations of 66 chromophobe renal cell carcinomas (ChRCCs) on the basis of multidimensional and comprehensive characterization, including mtDNA and whole-genome sequencing. The result is consistent that ChRCC originates from the distal nephron compared with other kidney cancers with more proximal origins. Combined mtDNA and gene expression analysis implicates changes in mitochondrial function as a component of the disease biology, while suggesting alternative roles for mtDNA mutations in cancers relying on oxidative phosphorylation. Genomic rearrangements lead to recurrent structural breakpoints within TERT promoter region, which correlates with highly elevated TERT expression and manifestation of kataegis, representing a mechanism of TERT upregulation in cancer distinct from previously observed amplifications and point mutations.

MicroRNAs (miRNAs) are small regulatory RNAs that can regulate gene expression by binding to mRNA sequences and repressing target-gene expression post-transcriptionally, either by inhibiting translation or promoting RNA degradation. We have analysed expression of 328 known and 152 novel human miRNAs in 10 benign peripheral zone tissues and 16 prostate cancer tissues using microarrays and found widespread, but not universal, downregulation of miRNAs in clinically localized prostate cancer relative to benign peripheral zone tissue. These findings have been verified by real-time RT-PCR assays on select miRNAs, including miR-125b, miR-145 and let-7c. The downregulated miRNAs include several with proven target mRNAs whose proteins have been previously shown to be increased in prostate cancer by immunohistochemistry, including RAS, E2F3, BCL-2 and MCL-1. Using a bioinformatics approach, we have identified additional potential mRNA targets of one of the miRNAs, (miR-125b) that are upregulated in prostate cancer and confirmed increased expression of one of these targets, EIF4EBP1, in prostate cancer tissues. Our findings indicate that changes in miRNA expression may have an important role in the biology of human prostate cancer.

Adoptive transfer of T lymphocytes expressing a CD19-specific chimeric antigen receptor (CAR.CD19) induces complete tumor regression in patients with lymphoid malignancies. Although in vivo persistence of CAR-T cells correlates with clinical responses, it remains unknown whether specific cell subsets within the CAR-T-cell product correlate with their subsequent in vivo expansion and persistence. We analyzed 14 patients with B-cell malignancies infused with autologous CAR.CD19-redirected T cells expanded ex vivo using IL-2, and found that their in vivo expansion only correlated with the frequency within the infused product of a CD8(+)CD45RA(+)CCR7(+) subset, whose phenotype is closest to "T-memory stem cells." Preclinical models showed that increasing the frequency of CD8(+)CD45RA(+)CCR7(+) CAR-T cells in the infused line by culturing the cells with IL-7 and IL-15 produced greater antitumor activity of CAR-T cells mediated by increased resistance to cell death, following repetitive encounters with the antigen, while preserving their migration to secondary lymphoid organs. This trial was registered at www.clinicaltrials.gov as #NCT00586391 and #NCT00709033.

Dysregulated metabolism is a hallmark of cancer, manifested through alterations in metabolites. We performed metabolomic profiling on 138 matched clear cell renal cell carcinoma (ccRCC)/normal tissue pairs and found that ccRCC is characterized by broad shifts in central carbon metabolism, one-carbon metabolism, and antioxidant response. Tumor progression and metastasis were associated with metabolite increases in glutathione and cysteine/methionine metabolism pathways. We develop an analytic pipeline and visualization tool (metabolograms) to bridge the gap between TCGA transcriptomic profiling and our metabolomic data, which enables us to assemble an integrated pathway-level metabolic atlas and to demonstrate discordance between transcriptome and metabolome. Lastly, expression profiling was performed on a high-glutathione cluster, which corresponds to a poor-survival subgroup in the ccRCC TCGA cohort.

Molecular alterations involving the PI3K/AKT/mTOR pathway (including mutation, copy number, protein, or RNA) were examined across 11,219 human cancers representing 32 major types. Within specific mutated genes, frequency, mutation hotspot residues, in silico predictions, and functional assays were all informative in distinguishing the subset of genetic variants more likely to have functional relevance. Multiple oncogenic pathways including PI3K/AKT/mTOR converged on similar sets of downstream transcriptional targets. In addition to mutation, structural variations and partial copy losses involving PTEN and STK11 showed evidence for having functional relevance. A substantial fraction of cancers showed high mTOR pathway activity without an associated canonical genetic or genomic alteration, including cancers harboring IDH1 or VHL mutations, suggesting multiple mechanisms for pathway activation.

Metastatic disease is a primary cause of cancer-related death, and factors governing tumor cell metastasis have not been fully elucidated. Here, we address this question by using tumor cell lines derived from mice that develop metastatic lung adenocarcinoma owing to expression of mutant K-ras and p53. Despite having widespread somatic genetic alterations, the metastasis-prone tumor cells retained a marked plasticity. They transited reversibly between epithelial and mesenchymal states, forming highly polarized epithelial spheres in three-dimensional culture that underwent epithelial-to-mesenchymal transition (EMT) following treatment with transforming growth factor-beta or injection into syngeneic mice. This transition was entirely dependent on the microRNA (miR)-200 family, which decreased during EMT. Forced expression of miR-200 abrogated the capacity of these tumor cells to undergo EMT, invade, and metastasize, and conferred transcriptional features of metastasis-incompetent tumor cells. We conclude that tumor cell metastasis is regulated by miR-200 expression, which changes in response to contextual extracellular cues.

Abstract Not all breast cancers respond to tamoxifen, and many develop resistance despite initial benefit. We used an in vivo model of estrogen receptor (ER)–positive breast cancer (MCF-7 xenografts) to investigate mechanisms of this resistance and develop strategies to circumvent it. Epidermal growth factor receptor (EGFR) and HER2, which were barely detected in control estrogen-treated tumors, increased slightly with tamoxifen and were markedly increased when tumors became resistant. Gefitinib, which inhibits EGFR/HER2, improved the antitumor effect of tamoxifen and delayed acquired resistance, but had no effect on estrogen-stimulated growth. Phosphorylated levels of p42/44 and p38 mitogen-activated protein kinases (both downstream of EGFR/HER2) were increased in the tamoxifen-resistant tumors and were suppressed by gefitinib. There was no apparent increase in phosphorylated AKT (also downstream of EGFR/HER2) in resistant tumors, but it was nonetheless suppressed by gefitinib. Phosphorylated insulin-like growth factor-IR (IGF-IR), which can interact with both EGFR and membrane ER, was elevated in the tamoxifen-resistant tumors compared with the sensitive group. However, ER-regulated gene products, including total IGF-IR itself and progesterone receptor, remained suppressed even at the time of acquired resistance. Tamoxifen's antagonism of classic ER genomic function was retained in these resistant tumors and even in tumors that overexpress HER2 (MCF-7 HER2/18) and are de novo tamoxifen-resistant. In conclusion, EGFR/HER2 may mediate tamoxifen resistance in ER-positive breast cancer despite continued suppression of ER genomic function by tamoxifen. IGF-IR expression remains dependent on ER but is activated in the tamoxifen-resistant tumors. This study provides a rationale to combine HER inhibitors with tamoxifen in clinical studies, even in tumors that do not initially overexpress EGFR/HER2. [Cancer Res 2008;68(3):826–33]

Astrocytes are the most abundant cell type in the brain, where they perform a wide array of functions, yet the nature of their cellular heterogeneity and how it oversees these diverse roles remains shrouded in mystery. Using an intersectional fluorescence-activated cell sorting-based strategy, we identified five distinct astrocyte subpopulations present across three brain regions that show extensive molecular diversity. Application of this molecular insight toward function revealed that these populations differentially support synaptogenesis between neurons. We identified correlative populations in mouse and human glioma and found that the emergence of specific subpopulations during tumor progression corresponded with the onset of seizures and tumor invasion. In sum, we have identified subpopulations of astrocytes in the adult brain and their correlates in glioma that are endowed with diverse cellular, molecular and functional properties. These populations selectively contribute to synaptogenesis and tumor pathophysiology, providing a blueprint for understanding diverse astrocyte contributions to neurological disease.

Aberrant expression of microRNAs (miRNAs) and the enzymes that control their processing have been reported in multiple biological processes including primary and metastatic tumours, but the mechanisms governing this are not clearly understood. Here we show that TAp63, a p53 family member, suppresses tumorigenesis and metastasis, and coordinately regulates Dicer and miR-130b to suppress metastasis. Metastatic mouse and human tumours deficient in TAp63 express Dicer at very low levels, and we found that modulation of expression of Dicer and miR-130b markedly affected the metastatic potential of cells lacking TAp63. TAp63 binds to and transactivates the Dicer promoter, demonstrating direct transcriptional regulation of Dicer by TAp63. These data provide a novel understanding of the roles of TAp63 in tumour and metastasis suppression through the coordinate transcriptional regulation of Dicer and miR-130b and may have implications for the many processes regulated by miRNAs.

Because of the high risk of recurrence in high-grade serous ovarian carcinoma (HGS-OvCa), the development of outcome predictors could be valuable for patient stratification.Using the catalog of The Cancer Genome Atlas (TCGA), we developed subtype and survival gene expression signatures, which, when combined, provide a prognostic model of HGS-OvCa classification, named "Classification of Ovarian Cancer" (CLOVAR).We validated CLOVAR on an independent dataset consisting of 879 HGS-OvCa expression profiles.The worst outcome group, accounting for 23% of all cases, was associated with a median survival of 23 months and a platinum resistance rate of 63%, versus a median survival of 46 months and platinum resistance rate of 23% in other cases.Associating the outcome prediction model with BRCA1/BRCA2 mutation status, residual disease after surgery, and disease stage further optimized outcome classification.Ovarian cancer is a disease in urgent need of more effective therapies.The spectrum of outcomes observed here and their association with CLOVAR signatures suggests variations in underlying tumor biology.Prospective validation of the CLOVAR model in the context of additional prognostic variables may provide a rationale for optimal combination of patient and treatment regimens.

Taking the Myc Despite nearly 30 years of research into the mechanisms by which Myc oncogene dysregulation contributes to tumorigenesis, there are still no effective therapies that inhibit Myc activity. Kessler et al. (p. 348 , published online 8 December; see the Perspective by Evan ) searched for gene products that support Myc-driven tumorigenesis. One pharmacologically tractable target that emerged from the screen was the SUMO-activating enzyme complex SAE1/2, which catalyzes a posttranslational modification (SUMOylation) that alters protein behavior and function. SUMOylation was found to control the Myc transcriptional response, and its inhibition caused mitotic defects and apoptosis in Myc-dependent breast cancer cells.

Abstract Motivation: Global gene expression profiling, both at the transcript level and at the protein level, can be a valuable tool in the understanding of genes, biological networks, and cellular states. As larger and larger gene expression data sets become available, data mining techniques can be applied to identify patterns of interest in the data. Association rules, used widely in the area of market basket analysis, can be applied to the analysis of expression data as well. Association rules can reveal biologically relevant associations between different genes or between environmental effects and gene expression. An association rule has the form LHS ⇒ RHS, where LHS and RHS are disjoint sets of items, the RHS set being likely to occur whenever the LHS set occurs. Items in gene expression data can include genes that are highly expressed or repressed, as well as relevant facts describing the cellular environment of the genes (e.g. the diagnosis of a tumor sample from which a profile was obtained). Results: We demonstrate an algorithm for efficiently mining association rules from gene expression data, using the data set from Hughes et al. (2000, Cell, 102, 109–126) of 300 expression profiles for yeast. Using the algorithm, we find numerous rules in the data. A cursory analysis of some of these rules reveals numerous associations between certain genes, many of which make sense biologically, others suggesting new hypotheses that may warrant further investigation. In a data set derived from the yeast data set, but with the expression values for each transcript randomly shifted with respect to the experiments, no rules were found, indicating that most all of the rules mined from the actual data set are not likely to have occurred by chance. Availability: An implementation of the algorithm using Microsoft SQL Server with Access 2000 is available at http://dot.ped.med.umich.edu:2000/pub/assoc_rules/assoc_rules.zip. Our results from mining the yeast data set are available at http://dot.ped.med.umich.edu:2000/pub/assoc_rules/yeast_results.zip. Contact: ccreight@umich.edu * To whom correspondence should be addressed.

Abstract Mass-spectrometry-based proteomic profiling of human cancers has the potential for pan-cancer analyses to identify molecular subtypes and associated pathway features that might be otherwise missed using transcriptomics. Here, we classify 532 cancers, representing six tissue-based types (breast, colon, ovarian, renal, uterine), into ten proteome-based, pan-cancer subtypes that cut across tumor lineages. The proteome-based subtypes are observable in external cancer proteomic datasets surveyed. Gene signatures of oncogenic or metabolic pathways can further distinguish between the subtypes. Two distinct subtypes both involve the immune system, one associated with the adaptive immune response and T-cell activation, and the other associated with the humoral immune response. Two additional subtypes each involve the tumor stroma, one of these including the collagen VI interacting network. Three additional proteome-based subtypes—respectively involving proteins related to Golgi apparatus, hemoglobin complex, and endoplasmic reticulum—were not reflected in previous transcriptomics analyses. A data portal is available at UALCAN website.

Although ovarian cancer is the most lethal gynecologic malignancy in women, little is known about how the cancer initiates and metastasizes. In the last decade, new evidence has challenged the dogma that the ovary is the main source of this cancer. The fallopian tube has been proposed instead as the primary origin of high-grade serous ovarian cancer, the subtype causing 70% of ovarian cancer deaths. By conditionally deleting Dicer, an essential gene for microRNA synthesis, and Pten, a key negative regulator of the PI3K pathway, we show that high-grade serous carcinomas arise from the fallopian tube in mice. In these Dicer-Pten double-knockout mice, primary fallopian tube tumors spread to engulf the ovary and then aggressively metastasize throughout the abdominal cavity, causing ascites and killing 100% of the mice by 13 mo. Besides the clinical resemblance to human serous cancers, these fallopian tube cancers highly express genes that are known to be up-regulated in human serous ovarian cancers, also demonstrating molecular similarities. Although ovariectomized mice continue to develop high-grade serous cancers, removal of the fallopian tube at an early age prevents cancer formation--confirming the fallopian tube origin of the cancer. Intriguingly, the primary carcinomas are first observed in the stroma of the fallopian tube, suggesting that these epithelial cancers have a mesenchymal origin. Thus, this mouse model demonstrates a paradigm for the origin and initiation of high-grade serous ovarian carcinomas, the most common and deadliest ovarian cancer.

Hippo signaling has been recognized as a key tumor suppressor pathway. Here, we perform a comprehensive molecular characterization of 19 Hippo core genes in 9,125 tumor samples across 33 cancer types using multidimensional "omic" data from The Cancer Genome Atlas. We identify somatic drivers among Hippo genes and the related microRNA (miRNA) regulators, and using functional genomic approaches, we experimentally characterize YAP and TAZ mutation effects and miR-590 and miR-200a regulation for TAZ. Hippo pathway activity is best characterized by a YAP/TAZ transcriptional target signature of 22 genes, which shows robust prognostic power across cancer types. Our elastic-net integrated modeling further reveals cancer-type-specific pathway regulators and associated cancer drivers. Our results highlight the importance of Hippo signaling in squamous cell cancers, characterized by frequent amplification of YAP/TAZ, high expression heterogeneity, and significant prognostic patterns. This study represents a systems-biology approach to characterizing key cancer signaling pathways in the post-genomic era.

On the basis of multidimensional and comprehensive molecular characterization (including DNA methalylation and copy number, RNA, and protein expression), we classified 894 renal cell carcinomas (RCCs) of various histologic types into nine major genomic subtypes. Site of origin within the nephron was one major determinant in the classification, reflecting differences among clear cell, chromophobe, and papillary RCC. Widespread molecular changes associated with TFE3 gene fusion or chromatin modifier genes were present within a specific subtype and spanned multiple subtypes. Differences in patient survival and in alteration of specific pathways (including hypoxia, metabolism, MAP kinase, NRF2-ARE, Hippo, immune checkpoint, and PI3K/AKT/mTOR) could further distinguish the subtypes. Immune checkpoint markers and molecular signatures of T cell infiltrates were both highest in the subtype associated with aggressive clear cell RCC. Differences between the genomic subtypes suggest that therapeutic strategies could be tailored to each RCC disease subset.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

MicroRNAs are short non-coding RNAs that regulate the stability and translation of mRNAs.Profiling experiments, using microarray or deep sequencing technology, have identified microRNAs that are preferentially expressed in certain tissues, specific stages of development, or disease states such as cancer.Deep sequencing utilizes massively parallel sequencing, generating millions of small RNA sequence reads from a given sample.Profiling of microRNAs by deep sequencing measures absolute abundance and allows for the discovery of novel microRNAs that have eluded previous cloning and standard sequencing efforts.Public databases provide in silico predictions of microRNA gene targets by various algorithms.To better determine which of these predictions represent true positives, microRNA expression data can be integrated with gene expression data to identify putative microRNA:mRNA functional pairs.Here we discuss tools and methodologies for the analysis of microRNA expression data from deep sequencing.

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smoking and/or human papillomavirus (HPV). SCCs harbor 3q, 5p, and other recurrent chromosomal copy-number alterations (CNAs), DNA mutations, and/or aberrant methylation of genes and microRNAs, which are correlated with the expression of multi-gene programs linked to squamous cell stemness, epithelial-to-mesenchymal differentiation, growth, genomic integrity, oxidative damage, death, and inflammation. Low-CNA SCCs tended to be HPV(+) and display hypermethylation with repression of TET1 demethylase and FANCF, previously linked to predisposition to SCC, or harbor mutations affecting CASP8, RAS-MAPK pathways, chromatin modifiers, and immunoregulatory molecules. We uncovered hypomethylation of the alternative promoter that drives expression of the ΔNp63 oncogene and embedded miR944. Co-expression of immune checkpoint, T-regulatory, and Myeloid suppressor cells signatures may explain reduced efficacy of immune therapy. These findings support possibilities for molecular classification and therapeutic approaches.

Abstract Mitochondria are essential cellular organelles that play critical roles in cancer. Here, as part of the International Cancer Genome Consortium/The Cancer Genome Atlas Pan-Cancer Analysis of Whole Genomes Consortium, which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumor types, we performed a multidimensional, integrated characterization of mitochondrial genomes and related RNA sequencing data. Our analysis presents the most definitive mutational landscape of mitochondrial genomes and identifies several hypermutated cases. Truncating mutations are markedly enriched in kidney, colorectal and thyroid cancers, suggesting oncogenic effects with the activation of signaling pathways. We find frequent somatic nuclear transfers of mitochondrial DNA, some of which disrupt therapeutic target genes. Mitochondrial copy number varies greatly within and across cancers and correlates with clinical variables. Co-expression analysis highlights the function of mitochondrial genes in oxidative phosphorylation, DNA repair and the cell cycle, and shows their connections with clinically actionable genes. Our study lays a foundation for translating mitochondrial biology into clinical applications.

Purpose To determine the pathologic complete response (pCR) rate in estrogen receptor (ER) –positive primary breast cancer triaged to chemotherapy when the protein encoded by the MKI67 gene (Ki67) level was > 10% after 2 to 4 weeks of neoadjuvant aromatase inhibitor (AI) therapy. A second objective was to examine risk of relapse using the Ki67-based Preoperative Endocrine Prognostic Index (PEPI). Methods The American College of Surgeons Oncology Group (ACOSOG) Z1031A trial enrolled postmenopausal women with stage II or III ER-positive (Allred score, 6 to 8) breast cancer whose treatment was randomly assigned to neoadjuvant AI therapy with anastrozole, exemestane, or letrozole. For the trial ACOSOG Z1031B, the protocol was amended to include a tumor Ki67 determination after 2 to 4 weeks of AI. If the Ki67 was > 10%, patients were switched to neoadjuvant chemotherapy. A pCR rate of > 20% was the predefined efficacy threshold. In patients who completed neoadjuvant AI, stratified Cox modeling was used to assess whether time to recurrence differed by PEPI = 0 score (T1 or T2, N0, Ki67 < 2.7%, ER Allred > 2) versus PEPI > 0 disease. Results Only two of the 35 patients in ACOSOG Z1031B who were switched to neoadjuvant chemotherapy experienced a pCR (5.7%; 95% CI, 0.7% to 19.1%). After 5.5 years of median follow-up, four (3.7%) of the 109 patients with a PEPI = 0 score relapsed versus 49 (14.4%) of 341 of patients with PEPI > 0 (recurrence hazard ratio [PEPI = 0 v PEPI > 0], 0.27; P = .014; 95% CI, 0.092 to 0.764). Conclusion Chemotherapy efficacy was lower than expected in ER-positive tumors exhibiting AI-resistant proliferation. The optimal therapy for these patients should be further investigated. For patients with PEPI = 0 disease, the relapse risk over 5 years was only 3.6% without chemotherapy, supporting the study of adjuvant endocrine monotherapy in this group. These Ki67 and PEPI triage approaches are being definitively studied in the ALTERNATE trial (Alternate Approaches for Clinical Stage II or III Estrogen Receptor Positive Breast Cancer Neoadjuvant Treatment in Postmenopausal Women: A Phase III Study; clinical trial information: NCT01953588).

Among breast cancers, triple-negative breast cancer (TNBC) is the most poorly understood and is refractory to current targeted therapies. Using a genetic screen, we identify the PTPN12 tyrosine phosphatase as a tumor suppressor in TNBC. PTPN12 potently suppresses mammary epithelial cell proliferation and transformation. PTPN12 is frequently compromised in human TNBCs, and we identify an upstream tumor-suppressor network that posttranscriptionally controls PTPN12. PTPN12 suppresses transformation by interacting with and inhibiting multiple oncogenic tyrosine kinases, including HER2 and EGFR. The tumorigenic and metastatic potential of PTPN12-deficient TNBC cells is severely impaired upon restoration of PTPN12 function or combined inhibition of PTPN12-regulated tyrosine kinases, suggesting that TNBCs are dependent on the proto-oncogenic tyrosine kinases constrained by PTPN12. Collectively, these data identify PTPN12 as a commonly inactivated tumor suppressor and provide a rationale for combinatorially targeting proto-oncogenic tyrosine kinases in TNBC and other cancers based on their profile of tyrosine-phosphatase activity.

Abstract MicroRNAs (miRNA) regulate complex patterns of gene expression, and the relevance of altered miRNA expression to ovarian cancer remains to be elucidated. By comprehensively profiling expression of miRNAs and mRNAs in serous ovarian tumors and cell lines and normal ovarian surface epithelium, we identified hundreds of potential miRNA-mRNA targeting associations underlying cancer. Functional overexpression of miR-31, the most underexpressed miRNA in serous ovarian cancer, repressed predicted miR-31 gene targets including the cell cycle regulator E2F2. MIR31 and CDKN2A, which encode p14ARF and p16INK4A, are located at 9p21.3, a genomic region commonly deleted in ovarian and other cancers. p14ARF promotes p53 activity, and E2F2 overexpression in p53 wild-type cells normally leads via p14ARF to an induction of p53-dependent apoptosis. In a number of serous cancer cell lines with a dysfunctional p53 pathway (i.e., OVCAR8, OVCA433, and SKOV3), miR-31 overexpression inhibited proliferation and induced apoptosis; however, in other lines (i.e., HEY and OVSAYO) with functional p53, miR-31 had no effect. Additionally, the osteosarcoma cell line U2OS and the prostate cancer cell line PC3 (p14ARF-deficient and p53-deficient, respectively) were also sensitive to miR-31. Furthermore, miR-31 overexpression induced a global gene expression pattern in OVCAR8 associated with better prognosis in tumors from patients with advanced stage serous ovarian cancer, potentially affecting many genes underlying disease progression. Our findings reveal that loss of miR-31 is associated with defects in the p53 pathway and functions in serous ovarian cancer and other cancers, suggesting that patients with cancers deficient in p53 activity might benefit from therapeutic delivery of miR-31. Cancer Res; 70(5); 1906–15

<h2>Summary</h2> Transmitochondrial cybrids and multiple OMICs approaches were used to understand mitochondrial reprogramming and mitochondria-regulated cancer pathways in triple-negative breast cancer (TNBC). Analysis of cybrids and established breast cancer (BC) cell lines showed that metastatic TNBC maintains high levels of ATP through fatty acid β oxidation (FAO) and activates Src oncoprotein through autophosphorylation at Y419. Manipulation of FAO including the knocking down of carnitine palmitoyltransferase-1A (CPT1) and 2 (CPT2), the rate-limiting proteins of FAO, and analysis of patient-derived xenograft models confirmed the role of mitochondrial FAO in Src activation and metastasis. Analysis of TCGA and other independent BC clinical data further reaffirmed the role of mitochondrial FAO and CPT genes in Src regulation and their significance in BC metastasis.

Accumulating evidence suggests that both levels and activity of the estrogen receptor (ER) and the progesterone receptor (PR) are dramatically influenced by growth-factor receptor (GFR) signaling pathways, and that this crosstalk is a major determinant of both breast cancer progression and response to therapy. The phosphatidylinositol 3-kinase (PI3K) pathway, a key mediator of GFR signaling, is one of the most altered pathways in breast cancer. We thus examined whether deregulated PI3K signaling in luminal ER+ breast tumors is associated with ER level and activity and intrinsic molecular subtype.We defined two independent molecular signatures of the PI3K pathway: a proteomic (reverse-phase proteomic array) PI3K signature, based on protein measurement for PI3K signaling intermediates, and a PI3K transcriptional (mRNA) signature based on the set of genes either induced or repressed by PI3K inhibitors. By using these signatures, we scored each ER+ breast tumor represented in multiple independent expression-profiling datasets (four mRNA, n = 915; one protein, n = 429) for activation of the PI3K pathway. Effects of PI3K inhibitor BEZ-235 on ER expression and activity levels and cell growth were tested by quantitative real-time PCR and cell proliferation assays.Within ER+ tumors, ER levels were negatively correlated with the PI3K activation scores, both at the proteomic and transcriptional levels, in all datasets examined. PI3K signature scores were also higher in ER+ tumors and cell lines of the more aggressive luminal B molecular subtype versus those of the less aggressive luminal A subtype. Notably, BEZ-235 treatment in four different ER+ cell lines increased expression of ER and ER target genes including PR, and treatment with IGF-I (which signals via PI3K) decreased expression of ER and target genes, thus further establishing an inverse functional relation between ER and PI3K. BEZ-235 had an additional effect on tamoxifen in inhibiting the growth of a number of ER+ cell lines.Our data suggest that luminal B tumors have hyperactive GFR/PI3K signaling associated with lower ER levels, which has been correlated with resistance to endocrine therapy. Targeting PI3K in these tumors might reverse loss of ER expression and signaling and restore hormonal sensitivity.

MicroRNAs (miRNAs) are short, non-coding RNAs that target and silence protein coding genes through 3'-UTR elements. Evidence increasingly assigns an immunosuppressive role for miRNAs in immunity, but relatively few miRNAs have been studied, and an overall understanding of the importance of these regulatory transcripts in complex in vivo systems is lacking. Here we have applied multiple technologies to globally analyze miRNA expression and function in allergic lung disease, an experimental model of asthma. Deep sequencing and microarray analyses of the mouse lung short RNAome revealed numerous extant and novel miRNAs and other transcript classes. Similar to mRNAs, lung miRNA expression changed dynamically during the transition from the naive to the allergic state, suggesting numerous functional relationships. A possible role for miRNA editing in altering the lung mRNA target repertoire was also identified. Multiple members of the highly conserved let-7 miRNA family were the most abundant lung miRNAs, and we confirmed in vitro that interleukin 13 (IL-13), a cytokine essential for expression for allergic lung disease, is regulated by mmu-let-7a. However, inhibition of let-7 miRNAs in vivo using a locked nucleic acid profoundly inhibited production of allergic cytokines and the disease phenotype. Our findings thus reveal unexpected complexity in the miRNAome underlying allergic lung disease and demonstrate a proinflammatory role for let-7 miRNAs.

Abstract MicroRNAs (miRNAs) are small noncoding RNAs that direct gene regulation through translational repression and degradation of complementary mRNA. Although miRNAs have been implicated as oncogenes and tumor suppressors in a variety of human cancers, functional roles for individual miRNAs have not been described in clear cell ovarian carcinoma, an aggressive and chemoresistant subtype of ovarian cancer. We performed deep sequencing to comprehensively profile miRNA expression in 10 human clear cell ovarian cancer cell lines compared with normal ovarian surface epithelial cultures and discovered 54 miRNAs that were aberrantly expressed. Because of the critical roles of the phosphatidylinositol 3-kinase/v-akt murine thymoma viral oncogene homolog 1/mammalian target of rapamycin (mTOR) pathway in clear cell ovarian cancer, we focused on mir-100, a putative tumor suppressor that was the most down-regulated miRNA in our cancer cell lines, and its up-regulated target, FRAP1/mTOR. Overexpression of mir-100 inhibited mTOR signaling and enhanced sensitivity to the rapamycin analog RAD001 (everolimus), confirming the key relationship between mir-100 and the mTOR pathway. Furthermore, overexpression of the putative tumor suppressor mir-22 repressed the EVI1 oncogene, which is known to suppress apoptosis by stimulating phosphatidylinositol 3-kinase/v-akt murine thymoma viral oncogene homolog 1 signaling. In addition to these specific effects, reversing the expression of mir-22 and the putative oncogene mir-182 had widespread effects on target and nontarget gene populations that ultimately caused a global shift in the cancer gene signature toward a more normal state. Our experiments have revealed strong candidate miRNAs and their target genes that may contribute to the pathogenesis of clear cell ovarian cancer, thereby highlighting alternative therapeutic strategies for the treatment of this deadly cancer.

Endometriosis is a common disease seen by gynecologists. Clinical features involve pelvic pain and unexplained infertility. Although endometriosis is pathologically characterized by endometrial tissue outside the normal uterine location, endometriosis is otherwise not easily explained. Endometriomas, endometriotic cysts of the ovary, typically cause pain and distortion of pelvic anatomy. To begin to understand the pathogenesis of endometriomas, we describe the first transcriptome-microRNAome analysis of endometriomas and eutopic endometrium using next-generation sequencing technology. Using this approach, we generated a total of more than 54 million independent small RNA reads from our 19 clinical samples. At the microRNA level, we found 10 microRNA that were up-regulated (miR-202, 193a-3p, 29c, 708, 509-3-5p, 574-3p, 193a-5p, 485-3p, 100, and 720) and 12 microRNA that were down-regulated (miR-504, 141, 429, 203, 10a, 200b, 873, 200c, 200a, 449b, 375, and 34c-5p) in endometriomas compared with endometrium. Using in silico prediction algorithms, we correlated these microRNA with their corresponding differentially expressed mRNA targets. To validate the functional roles of microRNA, we manipulated levels of miR-29c in an in vitro system of primary cultures of human endometrial stromal fibroblasts. Extracellular matrix genes that were potential targets of miR-29c in silico were significantly down-regulated using this biological in vitro system. In vitro functional studies using luciferase reporter constructs further confirmed that miR-29c directly affects specific extracellular matrix genes that are dysregulated in endometriomas. Thus, miR-29c and other abnormally regulated microRNA appear to play important roles in the pathophysiology of uterine function and dysfunction.

Vascular endothelial growth factor-A (VEGF), a potent angiogenic factor, is also implicated in self-renewal in several normal tissue types. VEGF has been shown to drive malignant stem cells but mechanisms thereof and tumor types affected are not fully characterized. Here, we show VEGF promotes breast and lung cancer stem cell (CSC) self-renewal via VEGF receptor-2 (VEGFR-2)/STAT3-mediated upregulation of Myc and Sox2. VEGF increased tumor spheres and aldehyde dehydrogenase activity, both proxies for stem cell function in vitro, in triple-negative breast cancer (TNBC) lines and dissociated primary cancers, and in lung cancer lines. VEGF exposure before injection increased breast cancer-initiating cell abundance in vivo yielding increased orthotopic tumors, and increased metastasis from orthotopic primaries and following tail vein injection without further VEGF treatment. VEGF rapidly stimulated VEGFR-2/JAK2/STAT3 binding and activated STAT3 to bind MYC and SOX2 promoters and induce their expression. VEGFR-2 knockdown or inhibition abrogated VEGF-mediated STAT3 activation, MYC and SOX2 induction and sphere formation. Notably, knockdown of either STAT3, MYC or SOX2 impaired VEGF-upregulation of pSTAT3, MYC and SOX2 expression and sphere formation. Each transcription factor, once upregulated, appears to promote sustained activation of the others, creating a feed-forward loop to drive self-renewal. Thus, in addition to angiogenic effects, VEGF promotes tumor-initiating cell self-renewal through VEGFR-2/STAT3 signaling. Analysis of primary breast and lung cancers (>1300 each) showed high VEGF expression, was prognostic of poor outcome and strongly associated with STAT3 and MYC expression, supporting the link between VEGF and CSC self-renewal. High-VEGF tumors may be most likely to escape anti-angiogenics by upregulating VEGF, driving CSC self-renewal to re-populate post-treatment. Our work highlights the need to better define VEGF-driven cancer subsets and supports further investigation of combined therapeutic blockade of VEGF or VEGFR-2 and JAK2/STAT3.

An elevated Neutrophil-to-lymphocyte ratio (NLR) is associated with worse outcomes in several malignancies. However, its role with contemporary immune checkpoint blockade (ICB) is unknown. We investigated the utility of NLR in metastatic renal cell carcinoma (mRCC) patients treated with PD-1/PD-L1 ICB.We examined NLR at baseline and 6 (±2) weeks later in 142 patients treated between 2009 and 2017 at Dana-Farber Cancer Institute (Boston, USA). Landmark analysis at 6 weeks was conducted to explore the prognostic value of relative NLR change on overall survival (OS), progression-free survival (PFS), and objective response rate (ORR). Cox and logistic regression models allowed for adjustment of line of therapy, number of IMDC risk factors, histology and baseline NLR.Median follow up was 16.6 months (range: 0.7-67.8). Median duration on therapy was 5.1 months (<1-61.4). IMDC risk groups were: 18% favorable, 60% intermediate, 23% poor-risk. Forty-four percent were on first-line ICB and 56% on 2nd line or more. Median NLR was 3.9 (1.3-42.4) at baseline and 4.1 (1.1-96.4) at week 6. Patients with a higher baseline NLR showed a trend toward lower ORR, shorter PFS, and shorter OS. Higher NLR at 6 weeks was a significantly stronger predictor of all three outcomes than baseline NLR. Relative NLR change by ≥25% from baseline to 6 weeks after ICB therapy was associated with reduced ORR and an independent prognostic factor for PFS (p < 0.001) and OS (p = 0.004), whereas a decrease in NLR by ≥25% was associated with improved outcomes.Early decline and NLR at 6 weeks are associated with significantly improved outcomes in mRCC patients treated with ICB. The prognostic value of the readily-available NLR warrants larger, prospective validation.

Allergy Induction Proteinases found in fungi and other allergens elicit allergic inflammation, but how they do so is far from clear. It is also unclear how pattern recognition receptors, which detect invading microbes, drive allergic inflammation. Millien et al. (p. 792 ) shed light on this puzzle by showing that, in mice, induction of allergic inflammation requires proteinase-dependent cleavage of the clotting factor fibrinogen, leading to generation of a ligand that activates the pattern-recognition receptor, Toll-like receptor 4 (TLR4). Cleaved fibrinogen signals through TLR4 to activate the innate immune system and recruit cells to the airway, which drives both allergic responses and antifungal immunity.

PBRM1 is the second most commonly mutated gene after VHL in clear cell renal cell carcinoma (ccRCC). However, the biological consequences of PBRM1 mutations for kidney tumorigenesis are unknown. Here, we find that kidney-specific deletion of Vhl and Pbrm1, but not either gene alone, results in bilateral, multifocal, transplantable clear cell kidney cancers. PBRM1 loss amplified the transcriptional outputs of HIF1 and STAT3 incurred by Vhl deficiency. Analysis of mouse and human ccRCC revealed convergence on mTOR activation, representing the third driver event after genetic inactivation of VHL and PBRM1. Our study reports a physiological preclinical ccRCC mouse model that recapitulates somatic mutations in human ccRCC and provides mechanistic and therapeutic insights into PBRM1 mutated subtypes of human ccRCC.

Lung cancer is the leading cause of cancer-related deaths, primarily due to distant metastatic disease. Metastatic lung cancer cells can undergo an epithelial-to-mesenchymal transition (EMT) regulated by various transcription factors, including a double-negative feedback loop between the microRNA-200 (miR-200) family and ZEB1, but the precise mechanisms by which ZEB1-dependent EMT promotes malignancy remain largely undefined. Although the cell-intrinsic effects of EMT are important for tumor progression, the reciprocal dynamic crosstalk between mesenchymal cancer cells and the extracellular matrix (ECM) is equally critical in regulating invasion and metastasis. Investigating the collaborative effect of EMT and ECM in the metastatic process reveals increased collagen deposition in metastatic tumor tissues as a direct consequence of amplified collagen gene expression in ZEB1-activated mesenchymal lung cancer cells. In addition, collagen fibers in metastatic lung tumors exhibit greater linearity and organization as a result of collagen crosslinking by the lysyl oxidase (LOX) family of enzymes. Expression of the LOX and LOXL2 isoforms is directly regulated by miR-200 and ZEB1, respectively, and their upregulation in metastatic tumors and mesenchymal cell lines is coordinated to that of collagen. Functionally, LOXL2, as opposed to LOX, is the principal isoform that crosslinks and stabilizes insoluble collagen deposition in tumor tissues. In turn, focal adhesion formation and FAK/SRC signaling is activated in mesenchymal tumor cells by crosslinked collagen in the ECM. Our study is the first to validate direct regulation of LOX and LOXL2 by the miR-200/ZEB1 axis, defines a novel mechanism driving tumor metastasis, delineates collagen as a prognostic marker, and identifies LOXL2 as a potential therapeutic target against tumor progression.

Intratumoral collagen cross-links heighten stromal stiffness and stimulate tumor cell invasion, but it is unclear how collagen cross-linking is regulated in epithelial tumors. To address this question, we used Kras(LA1) mice, which develop lung adenocarcinomas from somatic activation of a Kras(G12D) allele. The lung tumors in Kras(LA1) mice were highly fibrotic and contained cancer-associated fibroblasts (CAF) that produced collagen and generated stiffness in collagen gels. In xenograft tumors generated by injection of wild-type mice with lung adenocarcinoma cells alone or in combination with CAFs, the total concentration of collagen cross-links was the same in tumors generated with or without CAFs, but coinjected tumors had higher hydroxylysine aldehyde-derived collagen cross-links (HLCC) and lower lysine-aldehyde-derived collagen cross-links (LCCs). Therefore, we postulated that an LCC-to-HLCC switch induced by CAFs promotes the migratory and invasive properties of lung adenocarcinoma cells. To test this hypothesis, we created coculture models in which CAFs are positioned interstitially or peripherally in tumor cell aggregates, mimicking distinct spatial orientations of CAFs in human lung cancer. In both contexts, CAFs enhanced the invasive properties of tumor cells in three-dimensional (3D) collagen gels. Tumor cell aggregates that attached to CAF networks on a Matrigel surface dissociated and migrated on the networks. Lysyl hydroxylase 2 (PLOD2/LH2), which drives HLCC formation, was expressed in CAFs, and LH2 depletion abrogated the ability of CAFs to promote tumor cell invasion and migration.CAFs induce a collagen cross-link switch in tumor stroma to influence the invasive properties of tumor cells.

Abstract Despite the development of combined modality treatments against liposarcoma in recent years, a significant proportion of patients respond only modestly to such approaches, possibly contributing to local or distant recurrence. Early detection of recurrent or metastatic disease could improve patient prognosis by triggering earlier clinical intervention. However, useful biomarkers for such purposes are lacking. Using both patient plasma samples and cell lines, we demonstrate here that miR-25-3p and miR-92a-3p are secreted by liposarcoma cells through extracellular vesicles and may be useful as potential biomarkers of disease. Both miR-25-3p and miR-92a-3p stimulated secretion of proinflammatory cytokine IL6 from tumor-associated macrophages in a TLR7/8-dependent manner, which in turn promoted liposarcoma cell proliferation, invasion, and metastasis via this interaction with the surrounding microenvironment. Our findings provide novel and previously unreported insight into liposarcoma progression, identifying communication between liposarcoma cells and their microenvironment as a process critically involved in liposarcoma progression. This study establishes the possibility that the pattern of circulating miRNAs may identify recurrence prior to radiological detectability while providing insight into disease outcome and as a possible approach to monitor treatment efficacy. Cancer Res; 77(14); 3846–56. ©2017 AACR.

Glioblastoma is a universally lethal form of brain cancer that exhibits an array of pathophysiological phenotypes, many of which are mediated by interactions with the neuronal microenvironment1,2. Recent studies have shown that increases in neuronal activity have an important role in the proliferation and progression of glioblastoma3,4. Whether there is reciprocal crosstalk between glioblastoma and neurons remains poorly defined, as the mechanisms that underlie how these tumours remodel the neuronal milieu towards increased activity are unknown. Here, using a native mouse model of glioblastoma, we develop a high-throughput in vivo screening platform and discover several driver variants of PIK3CA. We show that tumours driven by these variants have divergent molecular properties that manifest in selective initiation of brain hyperexcitability and remodelling of the synaptic constituency. Furthermore, secreted members of the glypican (GPC) family are selectively expressed in these tumours, and GPC3 drives gliomagenesis and hyperexcitability. Together, our studies illustrate the importance of functionally interrogating diverse tumour phenotypes driven by individual, yet related, variants and reveal how glioblastoma alters the neuronal microenvironment.

Mass-spectrometry-based proteomic data on human tumors-combined with corresponding multi-omics data-present opportunities for systematic and pan-cancer proteogenomic analyses. Here, we assemble a compendium dataset of proteomics data of 2002 primary tumors from 14 cancer types and 17 studies. Protein expression of genes broadly correlates with corresponding mRNA levels or copy number alterations (CNAs) across tumors, but with notable exceptions. Based on unsupervised clustering, tumors separate into 11 distinct proteome-based subtypes spanning multiple tissue-based cancer types. Two subtypes are enriched for brain tumors, one subtype associating with MYC, Wnt, and Hippo pathways and high CNA burden, and another subtype associating with metabolic pathways and low CNA burden. Somatic alteration of genes in a pathway associates with higher pathway activity as inferred by proteome or transcriptome data. A substantial fraction of cancers shows high MYC pathway activity without MYC copy gain but with mutations in genes with noncanonical roles in MYC. Our proteogenomics survey reveals the interplay between genome and proteome across tumor lineages.

Substantial evidence implicates insulin-like growth factor-I (IGF-I) signaling in the development and progression of breast cancer. To more clearly elucidate the role of IGF in human breast cancer, we identified and then examined gene expression patterns of IGF-I-treated breast cancer cells.MCF-7 cells were stimulated with IGF-I for 3 or 24 hours and were profiled for greater than 22,000 RNA transcripts. We defined an IGF-I signature pattern of more than 800 genes that were up- or downregulated at both time points. The gene signature was examined in clinical breast tumors and in experimental models that represented other oncogenic pathways. The signature was correlated with clinical and pathologic variables and with patient outcome.IGF-I caused temporal changes in gene expression that were strongly associated with cell proliferation, metabolism, and DNA repair. Genes with early and sustained regulation by IGF-I were highly enriched for transcriptional targets of the estrogen receptor (ER), Ras/extracellular signal-related kinase 1/2, and phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin pathways. In three large, independent data sets of profiled human breast tumors, the IGF-I signature was manifested in the majority of ER-negative breast tumors and in a subset (approximately 25%) of ER-positive breast tumors. Patients who had tumors that manifested the IGF-I signature (including patients who did not receive adjuvant therapy) had a shorter time to a poor outcome event. The IGF gene signature was highly correlated with numerous poor prognostic factors and was one of the strongest indicators of disease outcome.Transcriptional targets of IGF-I represent pathways of increased aggressiveness and possibly hormone independence in clinical breast cancers.

Breast cancer presents as either estrogen receptor alpha (ERalpha) positive or negative, with ERalpha+ tumors responding to antiestrogen therapy and having a better prognosis. By themselves, mRNA expression signatures of estrogen regulation in ERalpha+ breast cancer cells do not account for the vast molecular differences observed between ERalpha+ and ERalpha- cancers. In ERalpha- tumors, overexpression of epidermal growth factor receptor (EGFR) or c-erbB-2, leading to increased growth factor signaling, is observed such that mitogen-activated protein (MAP) kinase (MAPK) is significantly hyperactivated compared with ERalpha+ breast cancer. In ERalpha+/progesterone receptor-positive, estrogen-dependent MCF-7 breast cancer cells, we stably overexpressed EGFR or constitutively active erbB-2, Raf, or MAP/extracellular signal-regulated kinase kinase, resulting in cell lines exhibiting hyperactivation of MAPK, estrogen-independent growth, and the reversible down-regulation of ERalpha expression. By global mRNA profiling, we found a "MAPK signature" of approximately 400 genes consistently up-regulated or down-regulated in each of the MAPK+ cell lines. In several independent profile data sets of human breast tumors, the in vitro MAPK signature was able to accurately distinguish ER+ from ER- tumors. In addition, our in vitro mRNA profile data revealed distinct mRNA signatures specific to either erbB-2 or EGFR activation. A subset of breast tumor profiles was found to share extensive similarities with either the erbB-2-specific or the EGFR-specific signatures. Our results confirm that increased MAPK activation causes loss of ERalpha expression and suggest that hyperactivation of MAPK plays a role in the generation of the ERalpha- phenotype in breast cancer. These MAPK+ cell lines are excellent models for investigating the underlying mechanisms behind the ERalpha- phenotype.

The Polycomb group (PcG) protein EZH2 possesses oncogenic properties for which the underlying mechanism is unclear. We integrated in vitro cell line, in vivo tumor profiling, and genome-wide location data to nominate key targets of EZH2. One of the candidates identified was ADRB2 (Adrenergic Receptor, Beta-2), a critical mediator of beta-adrenergic signaling. EZH2 is recruited to the ADRB2 promoter and represses ADRB2 expression. ADRB2 inhibition confers cell invasion and transforms benign prostate epithelial cells, whereas ADRB2 overexpression counteracts EZH2-mediated oncogenesis. ADRB2 is underexpressed in metastatic prostate cancer, and clinically localized tumors that express lower levels of ADRB2 exhibit a poor prognosis. Taken together, we demonstrate the power of integrating multiple diverse genomic data to decipher targets of disease-related genes.

Molecular subtypes of serous ovarian cancer have been recently described. Using data from independent datasets including over 900 primary tumour samples, we show that deregulation of the Let-7 pathway is specifically associated with the C5 molecular subtype of serous ovarian cancer. DNA copy number and gene expression of HMGA2, alleles of Let-7, LIN28, LIN28B, MYC, MYCN, DICER1, and RNASEN were measured using microarray and quantitative reverse transcriptase PCR. Immunohistochemistry was performed on 127 samples using tissue microarrays and anti-HMGA2 antibodies. Fluorescence in situ hybridisation of bacterial artificial chromosomes hybridized to 239 ovarian tumours was used to measure translocation at the LIN28B locus. Short interfering RNA knockdown in ovarian cell lines was used to test the functionality of associations observed. Four molecular subtypes (C1, C2, C4, C5) of high-grade serous ovarian cancers were robustly represented in each dataset and showed similar pattern of patient survival. We found highly specific activation of a pathway involving MYCN, LIN28B, Let-7 and HMGA2 in the C5 molecular subtype defined by MYCN amplification and over-expression, over-expression of MYCN targets including the Let-7 repressor LIN28B, loss of Let-7 expression and HMGA2 amplification and over-expression. DICER1, a known Let-7 target, and RNASEN were over-expressed in C5 tumours. We saw no evidence of translocation at the LIN28B locus in C5 tumours. The reported interaction between LIN28B and Let-7 was recapitulated by siRNA knockdown in ovarian cancer cell lines. Our results associate deregulation of MYCN and downstream targets, including Let-7 and oncofetal genes, with serous ovarian cancer. We define for the first time how elements of an oncogenic pathway, involving multiple genes that contribute to stem cell renewal, is specifically altered in a molecular subtype of serous ovarian cancer. By defining the drivers of a molecular subtype of serous ovarian cancers we provide a novel strategy for targeted therapeutic intervention.

Epithelial tumor cells transit to a mesenchymal state in response to extracellular cues, in a process known as epithelial-to-mesenchymal transition (EMT).The precise nature of these cues has not been fully defined, an important issue given that EMT is an early event in tumor metastasis.Here, we have found that a population of metastasis-prone mouse lung adenocarcinoma cells expresses Notch and Notch ligands and that the Notch ligand Jagged2 promotes metastasis.Mechanistically, Jagged2 was found to promote metastasis by increasing the expression of GATA-binding (Gata) factors, which suppressed expression of the microRNA-200 (miR-200) family of microRNAs that target the transcriptional repressors that drive EMT and thereby induced EMT.Reciprocally, miR-200 inhibited expression of Gata3, which reversed EMT and abrogated metastasis, suggesting that Gata3 and miR-200 are mutually inhibitory and have opposing effects on EMT and metastasis.Consistent with this, high levels of Gata3 expression correlated with EMT in primary tumors from 2 cohorts of lung adenocarcinoma patients.These findings reveal what we believe to be a novel Jagged2/miR-200-dependent pathway that mediates lung adenocarcinoma EMT and metastasis in mice and may have implications for the treatment of human epithelial tumors.

The microRNA-200 (miR-200) family is part of a gene expression signature that predicts poor prognosis in lung cancer patients. In a mouse model of K-ras/p53-mutant lung adenocarcinoma, miR-200 levels are suppressed in metastasis-prone tumor cells, and forced miR-200 expression inhibits tumor growth and metastasis, but the miR-200 target genes that drive lung tumorigenesis have not been fully elucidated. Here, we scanned the genome for putative miR-200 binding sites and found them in the 3'-untranslated region (3'-UTR) of 35 genes that are amplified in human cancer. Mining of a database of resected human lung adenocarcinomas revealed that the levels of one of these genes, Flt1/VEGFR1, correlate inversely with duration of survival. Forced miR-200 expression suppressed Flt1 levels in metastasis-prone lung adenocarcinoma cells derived from K-ras/p53-mutant mice, and negatively regulated the Flt1 3'-UTR in reporter assays. Cancer-associated fibroblasts (CAFs) isolated from murine lung adenocarcinomas secreted abundant VEGF and enhanced tumor cell invasion in coculture studies. CAF-induced tumor cell invasion was abrogated by VEGF neutralization or Flt1 knockdown in tumor cells. Flt1 knockdown decreased the growth and metastasis of tumor cells in syngeneic mice. We conclude that miR-200 suppresses lung tumorigenesis by targeting Flt1.

Tamoxifen has been a frontline treatment for estrogen receptor alpha (ERα)-positive breast tumors in premenopausal women. However, resistance to tamoxifen occurs in many patients. ER still plays a critical role in the growth of breast cancer cells with acquired tamoxifen resistance, suggesting that ERα remains a valid target for treatment of tamoxifen-resistant (Tam-R) breast cancer. In an effort to identify novel regulators of ERα signaling, through a small-scale siRNA screen against histone methyl modifiers, we found WHSC1, a histone H3K36 methyltransferase, as a positive regulator of ERα signaling in breast cancer cells. We demonstrated that WHSC1 is recruited to the ERα gene by the BET protein BRD3/4, and facilitates ERα gene expression. The small-molecule BET protein inhibitor JQ1 potently suppressed the classic ERα signaling pathway and the growth of Tam-R breast cancer cells in culture. Using a Tam-R breast cancer xenograft mouse model, we demonstrated in vivo anti-breast cancer activity by JQ1 and a strong long-lasting effect of combination therapy with JQ1 and the ER degrader fulvestrant. Taken together, we provide evidence that the epigenomic proteins BRD3/4 and WHSC1 are essential regulators of estrogen receptor signaling and are novel therapeutic targets for treatment of Tam-R breast cancer.

MicroRNAs are short (∼22 nucleotides) noncoding RNAs that regulate the stability and translation of mRNA targets. A number of computational algorithms have been developed to help predict which microRNAs are likely to regulate which genes. Gene expression profiling of biological systems where microRNAs might be active can yield hundreds of differentially expressed genes. The commonly used public microRNA target prediction databases facilitate gene-by-gene searches. However, integration of microRNA–mRNA target predictions with gene expression data on a large scale using these databases is currently cumbersome and time consuming for many researchers. We have developed a desktop software application which, for a given target prediction database, retrieves all microRNA:mRNA functional pairs represented by an experimentally derived set of genes. Furthermore, for each microRNA, the software computes an enrichment statistic for overrepresentation of predicted targets within the gene set, which could help to implicate roles for specific microRNAs and microRNA-regulated genes in the system under study. Currently, the software supports searching of results from PicTar, TargetScan, and miRanda algorithms. In addition, the software can accept any user-defined set of gene-to-class associations for searching, which can include the results of other target prediction algorithms, as well as gene annotation or gene-to-pathway associations. A search (using our software) of genes transcriptionally regulated in vitro by estrogen in breast cancer uncovered numerous targeting associations for specific microRNAs—above what could be observed in randomly generated gene lists—suggesting a role for microRNAs in mediating the estrogen response. The software and Excel VBA source code are freely available at http://sigterms.sourceforge.net .

Mutations in phosphatase and tensin homologue (PTEN) or genomic alterations in the phosphatidylinositol-3-OH kinase-signalling pathway are the most common genetic alterations reported in human prostate cancer. However, the precise mechanism underlying how indolent tumours with PTEN alterations acquire metastatic potential remains poorly understood. Recent studies suggest that upregulation of transforming growth factor (TGF)-β signalling triggered by PTEN loss will form a growth barrier as a defence mechanism to constrain prostate cancer progression, underscoring that TGF-β signalling might represent a pre-invasive checkpoint to prevent PTEN-mediated prostate tumorigenesis. Here we show that COUP transcription factor II (COUP-TFII, also known as NR2F2), a member of the nuclear receptor superfamily, serves as a key regulator to inhibit SMAD4-dependent transcription, and consequently overrides the TGF-β-dependent checkpoint for PTEN-null indolent tumours. Overexpression of COUP-TFII in the mouse prostate epithelium cooperates with PTEN deletion to augment malignant progression and produce an aggressive metastasis-prone tumour. The functional counteraction between COUP-TFII and SMAD4 is reinforced by genetically engineered mouse models in which conditional loss of SMAD4 diminishes the inhibitory effects elicited by COUP-TFII ablation. The biological significance of COUP-TFII in prostate carcinogenesis is substantiated by patient sample analysis, in which COUP-TFII expression or activity is tightly correlated with tumour recurrence and disease progression, whereas it is inversely associated with TGF-β signalling. These findings reveal that the destruction of the TGF-β-dependent barrier by COUP-TFII is crucial for the progression of PTEN-mutant prostate cancer into a life-threatening disease, and supports COUP-TFII as a potential drug target for the intervention of metastatic human prostate cancer.

Abstract Purpose: Marked reactive stroma formation, designated as grade 3 reactive stroma, is associated with poor outcome in clinically localized prostate cancer. To understand the biological processes and signaling mechanisms underlying the formation of such reactive stroma, we carried out microarray gene expression analysis of laser-captured reactive stroma and matched normal stroma. Experimental Design: Seventeen cases of reactive stroma grade 3 cancer were used to laser-capture tumor and normal stroma. Expression analysis was carried out using Agilent 44K arrays. Up-regulation of selected genes was confirmed by quantitative reverse transcription-PCR. Expression data was analyzed to identify significantly up- and down-regulated genes, and gene ontology analysis was used to define pathways altered in reactive stroma. Results: A total of 544 unique genes were significantly higher in the reactive stroma and 606 unique genes were lower. Gene ontology analysis revealed significant alterations in a number of novel processes in prostate cancer reactive stroma, including neurogenesis, axonogenesis, and the DNA damage/repair pathways, as well as evidence of increases in stem cells in prostate cancer reactive stroma. Conclusions: Formation of reactive stroma in prostate cancer is a dynamic process characterized by significant alterations in growth factor and signal transduction pathways and formation of new structures, including nerves and axons.

Molecular profiling studies have found that estrogen receptor-positive (ER+) human breast cancers are comprised of at least two distinct diseases with differing biologies.With the advent of DNA microarrays, global gene expression patterns were used to define the luminal A and luminal B subtypes of ER+ breast cancer, with luminal B cancers showing a more aggressive phenotype including substantially worse outcomes in patients.The luminal B subtype designation could be considered a surrogate for those ER+ tumors having low progesterone receptors, high proliferation, high grade, and predicted poor response to hormone therapy.While they express estrogen receptors, luminal B cancers do not show a corresponding expression of estrogenregulated genes, and may therefore rely upon alternative pathways for growth.At the molecular level, luminal B cancers appear dramatically distinct from luminal A cancers, at the levels of gene expression, gene copy, somatic mutation, and DNA methylation; luminal B cancers are also genetically and genomically altered to a greater extent than luminal A cancers.While, in the clinical setting, luminal B is typically regarded as an ER+, hormone-sensitive disease, more research is needed into how to better treat it.Comprehensive profiling initiatives, such as The Cancer Genome Atlas, have recently provided us a catalog of mutated or copy altered genes, from which new therapeutic targets could potentially be mined.Candidate pathways that might be targeted in luminal B include those involving growth factor receptors, including HER2 and EGFR, as well as PI3K/Akt/mTor.

Metastatic cancer is extremely difficult to treat, and the presence of metastases greatly reduces a cancer patient's likelihood of long-term survival. The ZEB1 transcriptional repressor promotes metastasis through downregulation of microRNAs (miRs) that are strong inducers of epithelial differentiation and inhibitors of stem cell factors. Given that each miR can target multiple genes with diverse functions, we posited that the prometastatic network controlled by ZEB1 extends beyond these processes. We tested this hypothesis using a mouse model of human lung adenocarcinoma metastasis driven by ZEB1, human lung carcinoma cells, and human breast carcinoma cells. Transcriptional profiling studies revealed that ZEB1 controls the expression of numerous oncogenic and tumor-suppressive miRs, including miR-34a. Ectopic expression of miR-34a decreased tumor cell invasion and metastasis, inhibited the formation of promigratory cytoskeletal structures, suppressed activation of the RHO GTPase family, and regulated a gene expression signature enriched in cytoskeletal functions and predictive of outcome in human lung adenocarcinomas. We identified several miR-34a target genes, including Arhgap1, which encodes a RHO GTPase activating protein that was required for tumor cell invasion. These findings demonstrate that ZEB1 drives prometastatic actin cytoskeletal remodeling by downregulating miR-34a expression and provide a compelling rationale to develop miR-34a as a therapeutic agent in lung cancer patients.

Large-scale sequencing efforts are uncovering the complexity of cancer genomes, which are composed of causal "driver" mutations that promote tumor progression along with many more pathologically neutral "passenger" events. The majority of mutations, both in known cancer drivers and uncharacterized genes, are generally of low occurrence, highlighting the need to functionally annotate the long tail of infrequent mutations present in heterogeneous cancers. Here we describe a mutation assessment pipeline enabled by high-throughput engineering of molecularly barcoded gene variant expression clones identified by tumor sequencing. We first used this platform to functionally assess tail mutations observed in PIK3CA, which encodes the catalytic subunit alpha of the phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K) frequently mutated in cancer. Orthogonal screening for PIK3CA variant activity using in vitro and in vivo cell growth and transformation assays differentiated driver from passenger mutations, revealing that PIK3CA variant activity correlates imperfectly with its mutation frequency across breast cancer populations. Although PIK3CA mutations with frequencies above 5% were significantly more oncogenic than wild-type in all assays, mutations occurring at 0.07% to 5.0% included those with and without oncogenic activities that ranged from weak to strong in at least one assay. Proteomic profiling coupled with therapeutic sensitivity assays on PIK3CA variant-expressing cell models revealed variant-specific activation of PI3K signaling as well as other pathways that include the MEK1/2 module of mitogen-activated protein kinase pathway. Our data indicate that cancer treatments will need to increasingly consider the functional relevance of specific mutations in driver genes rather than considering all mutations in drivers as equivalent.

Epithelial tumor metastasis is preceded by an accumulation of collagen cross-links that heighten stromal stiffness and stimulate the invasive properties of tumor cells. However, the biochemical nature of collagen cross-links in cancer is still unclear. Here, we postulated that epithelial tumorigenesis is accompanied by changes in the biochemical type of collagen cross-links. Utilizing resected human lung cancer tissues and a p21CIP1/WAF1-deficient, K-rasG12D-expressing murine metastatic lung cancer model, we showed that, relative to normal lung tissues, tumor stroma contains higher levels of hydroxylysine aldehyde-derived collagen cross-links (HLCCs) and lower levels of lysine aldehyde-derived cross-links (LCCs), which are the predominant types of collagen cross-links in skeletal tissues and soft tissues, respectively. Gain- and loss-of-function studies in tumor cells showed that lysyl hydroxylase 2 (LH2), which hydroxylates telopeptidyl lysine residues on collagen, shifted the tumor stroma toward a high-HLCC, low-LCC state, increased tumor stiffness, and enhanced tumor cell invasion and metastasis. Together, our data indicate that LH2 enhances the metastatic properties of tumor cells and functions as a regulatory switch that controls the relative abundance of biochemically distinct types of collagen cross-links in the tumor stroma.

Signal transducer and activator of transcription (STAT) 3 regulates many cardinal features of cancer including cancer cell growth, apoptosis resistance, DNA damage response, metastasis, immune escape, tumor angiogenesis, the Warburg effect and oncogene addiction and has been validated as a drug target for cancer therapy. Several strategies have been used to identify agents that target Stat3 in breast cancer but none has yet entered into clinical use. We used a high-throughput fluorescence microscopy search strategy to identify compounds in a drug-repositioning library (Prestwick library) that block ligand-induced nuclear translocation of Stat3 and identified piperlongumine (PL), a natural product isolated from the fruit of the pepper Piper longum. PL inhibited Stat3 nuclear translocation, inhibited ligand-induced and constitutive Stat3 phosphorylation, and modulated expression of multiple Stat3-regulated genes. Surface plasmon resonance assay revealed that PL directly inhibited binding of Stat3 to its phosphotyrosyl peptide ligand. Phosphoprotein antibody array analysis revealed that PL does not modulate kinases known to activate Stat3 such as Janus kinases, Src kinase family members or receptor tyrosine kinases. PL inhibited anchorage-independent and anchorage-dependent growth of multiple breast cancer cell lines having increased pStat3 or total Stat3, and induced apoptosis. PL also inhibited mammosphere formation by tumor cells from patient-derived xenografts. PL’s antitumorigenic function was causally linked to its Stat3-inhibitory effect. PL was non-toxic in mice up to a dose of 30 mg/kg/day for 14 days and caused regression of breast cancer cell line xenografts in nude mice. Thus, PL represents a promising new agent for rapid entry into the clinic for use in treating breast cancer, as well as other cancers in which Stat3 has a role.

Abstract The microRNA-200 (miR-200) family restricts epithelial–mesenchymal transition (EMT) and metastasis in tumor cell lines derived from mice that develop metastatic lung adenocarcinoma. To determine the mechanisms responsible for EMT and metastasis regulated by this microRNA, we conducted a global liquid chromatography/tandem mass spectrometry analysis to compare metastatic and nonmetastatic murine lung adenocarcinoma cells which had undergone EMT because of loss of miR-200. An analysis of syngeneic tumors generated by these cells identified multiple novel proteins linked to metastasis. In particular, the analysis of conditioned media, cell surface proteins, and whole-cell lysates from metastatic and nonmetastatic cells revealed large-scale modifications in the tumor microenvironment. Specific increases were documented in extracellular matrix (ECM) proteins, peptidases, and changes in distribution of cell adhesion proteins in the metastatic cell lines. Integrating proteomic data from three subproteomes, we defined constituents of a multilayer protein network that both regulated and mediated the effects of TGFβ. Lastly, we identified ECM proteins and peptidases that were directly regulated by miR-200. Taken together, our results reveal how expression of miR-200 alters the tumor microenvironment to inhibit the processes of EMT and metastasis. Cancer Res; 71(24); 7670–82. ©2011 AACR.

Malignant peripheral nerve sheath tumors (MPNSTs) are soft tissue sarcomas that arise in connective tissue surrounding peripheral nerves. They occur sporadically in a subset of patients with neurofibromatosis type 1 (NF1). MPNSTs are highly aggressive, therapeutically resistant, and typically fatal. Using comparative transcriptome analysis, we identified CXCR4, a G-protein-coupled receptor, as highly expressed in mouse models of NF1-deficient MPNSTs, but not in nontransformed precursor cells. The chemokine receptor CXCR4 and its ligand, CXCL12, promote MPNST growth by stimulating cyclin D1 expression and cell-cycle progression through PI3-kinase (PI3K) and β-catenin signaling. Suppression of CXCR4 activity either by shRNA or pharmacological inhibition decreases MPNST cell growth in culture and inhibits tumorigenesis in allografts and in spontaneous genetic mouse models of MPNST. We further demonstrate conservation of these activated molecular pathways in human MPNSTs. Our findings indicate a role for CXCR4 in NF1-associated MPNST development and identify a therapeutic target.

Primary myelofibrosis (PMF) is a fatal neoplastic disease characterized by clonal myeloproliferation and progressive bone marrow (BM) fibrosis thought to be induced by mesenchymal stromal cells stimulated by overproduced growth factors. However, tissue fibrosis in other diseases is associated with monocyte-derived fibrocytes. Therefore, we sought to determine whether fibrocytes play a role in the induction of BM fibrosis in PMF. In this study, we show that BM from patients with PMF harbors an abundance of clonal, neoplastic collagen- and fibronectin-producing fibrocytes. Immunodeficient mice transplanted with myelofibrosis patients' BM cells developed a lethal myelofibrosis-like phenotype. Treatment of the xenograft mice with the fibrocyte inhibitor serum amyloid P (SAP; pentraxin-2) significantly prolonged survival and slowed the development of BM fibrosis. Collectively, our data suggest that neoplastic fibrocytes contribute to the induction of BM fibrosis in PMF, and inhibiting fibrocyte differentiation with SAP may interfere with this process.

Tumor metabolism is reorganized to support proliferation in the face of growth-related stress. Unlike the widespread profiling of changes to metabolic enzyme levels in cancer, comparatively less attention has been paid to the substrates/products of enzyme-catalyzed reactions, small-molecule metabolites. We developed an informatic pipeline to concurrently analyze metabolomics data from over 900 tissue samples spanning seven cancer types, revealing extensive heterogeneity in metabolic changes relative to normal tissue across cancers of different tissues of origin. Despite this heterogeneity, a number of metabolites were recurrently differentially abundant across many cancers, such as lactate and acyl-carnitine species. Through joint analysis of metabolomic data alongside clinical features of patient samples, we also identified a small number of metabolites, including several polyamines and kynurenine, which were associated with aggressive tumors across several tumor types. Our findings offer a glimpse onto common patterns of metabolic reprogramming across cancers, and the work serves as a large-scale resource accessible via a web application (http://www.sanderlab.org/pancanmet).

Abstract: Epithelial–mesenchymal transition (EMT) is involved in normal developmental cellular processes, but it may also be co-opted by a subset of cancer cells, to enable them to invade and form metastases at distant sites. Several gene transcription factors regulate EMT, including Snail1, Snail2, Zeb1, Zeb2, and Twist; ongoing studies continue to identify and elucidate other drivers. Specific micro ribonucleic acids (RNAs) have also been found to regulate EMT, including the microRNA-200 (miR-200) family, which targets Zeb1/Zeb2. Cancer “stem cells” – with the ability to self-renew and to regenerate all the cell types within the tumor – have been found to express EMT markers, further implicating both cancer stem cells and EMT with metastasis. Microenvironmental cues, including transforming growth factor-β, can direct EMT tumor metastasis, such as by regulating miR-200 expression. In human tumors, EMT markers and regulators may be expressed in a subset of tumor cells, such as in cells at the invasive front or tumor–microenvironment interface, though certain subtypes of cancer can show widespread mesenchymal-like features. In terms of therapeutic targeting of EMT in patients, potential areas of exploration could include targeting the cancer stem cell subpopulation, as well as microRNA-based therapeutics that reintroduce miR-200. This review will examine evidence for a role of EMT in invasion and metastasis, with the focus being on studies in lung and breast cancers. We also carry out analyses of publicly-available gene expression profiling datasets in order to show how EMT-associated genes appear coordinately expressed across human tumor specimens. Keywords: EMT, epithelial–mesenchymal transition, tumor microenvironment, miR-200, cancer stem cells

Abstract Bone metastasis from prostate cancer can occur years after prostatectomy, due to reactivation of dormant disseminated tumor cells (DTC) in the bone, yet the mechanism by which DTCs are initially induced into a dormant state in the bone remains to be elucidated. We show here that the bone microenvironment confers dormancy to C4-2B4 prostate cancer cells, as they become dormant when injected into mouse femurs but not under the skin. Live-cell imaging of dormant cells at the single-cell level revealed that conditioned medium from differentiated, but not undifferentiated, osteoblasts induced C4-2B4 cellular quiescence, suggesting that differentiated osteoblasts present locally around the tumor cells in the bone conferred dormancy to prostate cancer cells. Gene array analyses identified GDF10 and TGFβ2 among osteoblast-secreted proteins that induced quiescence of C4-2B4, C4-2b, and PC3-mm2, but not 22RV1 or BPH-1 cells, indicating prostate cancer tumor cells differ in their dormancy response. TGFβ2 and GDF10 induced dormancy through TGFβRIII to activate phospho-p38MAPK, which phosphorylates retinoblastoma (RB) at the novel N-terminal S249/T252 sites to block prostate cancer cell proliferation. Consistently, expression of dominant-negative p38MAPK in C4-2b and C4-2B4 prostate cancer cell lines abolished tumor cell dormancy both in vitro and in vivo. Lower TGFβRIII expression in patients with prostate cancer correlated with increased metastatic potential and decreased survival rates. Together, our results identify a dormancy mechanism by which DTCs are induced into a dormant state through TGFβRIII–p38MAPK–pS249/pT252–RB signaling and offer a rationale for developing strategies to prevent prostate cancer recurrence in the bone. Significance: These findings provide mechanistic insights into the dormancy of metastatic prostate cancer in the bone and offer a rationale for developing strategies to prevent prostate cancer recurrence in the bone. Cancer Res; 78(11); 2911–24. ©2018 AACR.

Forkhead box protein A1 (FOXA1) is a pioneer factor of estrogen receptor α (ER)-chromatin binding and function, yet its aberration in endocrine-resistant (Endo-R) breast cancer is unknown. Here, we report preclinical evidence for a role of FOXA1 in Endo-R breast cancer as well as evidence for its clinical significance. FOXA1 is gene-amplified and/or overexpressed in Endo-R derivatives of several breast cancer cell line models. Induced FOXA1 triggers oncogenic gene signatures and proteomic profiles highly associated with endocrine resistance. Integrated omics data reveal IL8 as one of the most perturbed genes regulated by FOXA1 and ER transcriptional reprogramming in Endo-R cells. IL-8 knockdown inhibits tamoxifen-resistant cell growth and invasion and partially attenuates the effect of overexpressed FOXA1. Our study highlights a role of FOXA1 via IL-8 signaling as a potential therapeutic target in FOXA1-overexpressing ER-positive tumors.

Abstract Radiotherapy is one of the mainstays of anticancer treatment, but the relationship between the radiosensitivity of cancer cells and their genomic characteristics is still not well defined. Here, we report the development of a high-throughput platform for measuring radiation survival in vitro and its validation in comparison with conventional clonogenic radiation survival analysis. We combined results from this high-throughput assay with genomic parameters in cell lines from squamous cell lung carcinoma, which is standardly treated by radiotherapy, to identify parameters that predict radiation sensitivity. We showed that activation of NFE2L2, a frequent event in lung squamous cancers, confers radiation resistance. An expression-based, in silico screen nominated inhibitors of phosphoinositide 3-kinase (PI3K) as NFE2L2 antagonists. We showed that the selective PI3K inhibitor, NVP-BKM120, both decreased NRF2 protein levels and sensitized NFE2L2 or KEAP1-mutant cells to radiation. We then combined results from this high-throughput assay with single-sample gene set enrichment analysis of gene expression data. The resulting analysis identified pathways implicated in cell survival, genotoxic stress, detoxification, and innate and adaptive immunity as key correlates of radiation sensitivity. The integrative and high-throughput methods shown here for large-scale profiling of radiation survival and genomic features of solid-tumor–derived cell lines should facilitate tumor radiogenomics and the discovery of genotype-selective radiation sensitizers and protective agents. Cancer Res; 73(20); 6289–98. ©2013 AACR.

Metastatic lung cancer is one of the most lethal forms of cancer and molecular pathways driving metastasis are still not clearly elucidated. Metastatic cancer cells undergo an epithelial–mesenchymal transition (EMT) where they lose their epithelial properties and acquire a migratory and invasive phenotype. Here we identify that the expression of microRNAs from the miR-200 family and the miR-183~96~182 cluster are significantly co-repressed in non-small cell lung cancer cell lines and primary tumors from multiple TCGA dataset with high EMT scores. Ectopic expression of the miR-183~96~182 cluster inhibited cancer cell migration and invasion, whereas its expression was tightly modulated by miR-200. We identified Foxf2 as a common, novel and direct target of both these microRNA families. Foxf2 expression tightly correlates with the transcription factor Zeb1 and is elevated in mesenchymal-like metastatic lung cancer cells. Foxf2 expression induced robust EMT, migration, invasion and metastasis in lung cancer cells, whereas Foxf2 inhibition significantly repressed these phenotypes. We also demonstrated that Foxf2 transcriptionally represses E-cadherin and miR-200, independent of Zeb1, to form a double-negative feedback loop. We, therefore, identified a novel mechanism whereby the miR-200 family and the miR-183~96~182 cluster inhibit lung cancer invasion and metastasis by targeting Foxf2.

Abstract Purpose: We previously reported an insulin-like growth factor (IGF) gene expression signature, based on genes induced or repressed by IGF-I, which correlated with poor prognosis in breast cancer. We tested whether the IGF signature was affected by anti–IGF-I receptor (IGF-IR) inhibitors and whether the IGF signature correlated with response to a dual anti–IGF-IR/insulin receptor (InsR) inhibitor, BMS-754807. Experimental Design: An IGF gene expression signature was examined in human breast tumors and cell lines and changes were noted following treatment of cell lines or xenografts with anti–IGF-IR antibodies or tyrosine kinase inhibitors. Sensitivity of cells to BMS-754807 was correlated with levels of the IGF signature. Human primary tumorgrafts were analyzed for the IGF signature and IGF-IR levels and activity, and MC1 tumorgrafts were treated with BMS-754807 and chemotherapy. Results: The IGF gene expression signature was reversed in three different models (cancer cell lines or xenografts) treated with three different anti–IGF-IR therapies. The IGF signature was present in triple-negative breast cancers (TNBC) and TNBC cell lines, which were especially sensitive to BMS-754807, and sensitivity was significantly correlated to the expression of the IGF gene signature. The TNBC primary human tumorgraft MC1 showed high levels of both expression and activity of IGF-IR and IGF gene signature score. Treatment of MC1 with BMS-754807 showed growth inhibition and, in combination with docetaxel, tumor regression occurred until no tumor was palpable. Regression was associated with reduced proliferation, increased apoptosis, and mitotic catastrophe. Conclusions: These studies provide a clear biological rationale to test anti–IGF-IR/InsR therapy in combination with chemotherapy in patients with TNBC. Clin Cancer Res; 17(8); 2314–27. ©2011 AACR.

Invasive bladder cancer, for which there have been few therapeutic advances in the past 20 years, is a significant medical problem associated with metastatic disease and frequent mortality. Although previous studies had identified many genetic alterations in invasive bladder cancer, recent genome-wide studies have provided a more comprehensive view. Here, we review those recent findings and suggest therapeutic strategies. Bladder cancer has a high mutation rate, exceeded only by lung cancer and melanoma. About 65% of all mutations are due to APOBEC-mediated mutagenesis. There is a high frequency of mutations and/or genomic amplification or deletion events that affect many of the canonical signaling pathways involved in cancer development: cell cycle, receptor tyrosine kinase, RAS, and PI-3-kinase/mTOR. In addition, mutations in chromatin-modifying genes are unusually frequent in comparison with other cancers, and mutation or amplification of transcription factors is also common. Expression clustering analyses organize bladder cancers into four principal groups, which can be characterized as luminal, immune undifferentiated, luminal immune, and basal. The four groups show markedly different expression patterns for urothelial differentiation (keratins and uroplakins) and immunity genes (CD274 and CTLA4), among others. These observations suggest numerous therapeutic opportunities, including kinase inhibitors and antibody therapies for genes in the canonical signaling pathways, histone deacetylase inhibitors and novel molecules for chromatin gene mutations, and immune therapies, which should be targeted to specific patients based on genomic profiling of their cancers.

The Cancer Genome Atlas (TCGA) Network recently comprehensively catalogued the molecular aberrations in 487 high-grade serous ovarian cancers, with much remaining to be elucidated regarding the microRNAs (miRNAs). Here, using TCGA ovarian data, we surveyed the miRNAs, in the context of their predicted gene targets.Integration of miRNA and gene patterns yielded evidence that proximal pairs of miRNAs are processed from polycistronic primary transcripts, and that intronic miRNAs and their host gene mRNAs derive from common transcripts. Patterns of miRNA expression revealed multiple tumor subtypes and a set of 34 miRNAs predictive of overall patient survival. In a global analysis, miRNA:mRNA pairs anti-correlated in expression across tumors showed a higher frequency of in silico predicted target sites in the mRNA 3'-untranslated region (with less frequency observed for coding sequence and 5'-untranslated regions). The miR-29 family and predicted target genes were among the most strongly anti-correlated miRNA:mRNA pairs; over-expression of miR-29a in vitro repressed several anti-correlated genes (including DNMT3A and DNMT3B) and substantially decreased ovarian cancer cell viability.This study establishes miRNAs as having a widespread impact on gene expression programs in ovarian cancer, further strengthening our understanding of miRNA biology as it applies to human cancer. As with gene transcripts, miRNAs exhibit high diversity reflecting the genomic heterogeneity within a clinically homogeneous disease population. Putative miRNA:mRNA interactions, as identified using integrative analysis, can be validated. TCGA data are a valuable resource for the identification of novel tumor suppressive miRNAs in ovarian as well as other cancers.

The ampulla of Vater is a complex cellular environment from which adenocarcinomas arise to form a group of histopathologically heterogenous tumors. To evaluate the molecular features of these tumors, 98 ampullary adenocarcinomas were evaluated and compared to 44 distal bile duct and 18 duodenal adenocarcinomas. Genomic analyses revealed mutations in the WNT signaling pathway among half of the patients and in all three adenocarcinomas irrespective of their origin and histological morphology. These tumors were characterized by a high frequency of inactivating mutations of ELF3, a high rate of microsatellite instability, and common focal deletions and amplifications, suggesting common attributes in the molecular pathogenesis are at play in these tumors. The high frequency of WNT pathway activating mutation, coupled with small-molecule inhibitors of β-catenin in clinical trials, suggests future treatment decisions for these patients may be guided by genomic analysis.

Osteosarcoma (OS) is the most frequent pediatric malignant bone tumor that has a high propensity for metastases. Through osteoblast-specific alteration of p53 status, we developed a genetically engineered mouse model of localized and metastatic OS to gain an understanding into the molecular pathogenesis of OS. Microarray analysis of both localized tumors and metastatic tumors identified the downregulation of the naked cuticle homolog 2 (NKD2) gene, a negative regulator of Wnt signaling. Overexpression of NKD2 in metastatic human and mouse OS cells significantly decreases cell proliferation, migration and invasion ability in vitro and drastically diminishes OS tumor growth and metastasis in vivo, whereas downregulation enhances migratory and invasive potential. Evaluation of NKD2-overexpressing tumors revealed upregulation of tumor-suppressor genes and downregulation of molecules involved in blood vessel formation and cell migration. Furthermore, assessment of primary human OS revealed downregulation of NKD2 in metastatic and recurrent OS. Finally, we provide biological evidence that use of small-molecule inhibitors targeting the Wnt pathway can have therapeutic efficacy in decreasing metastatic properties in OS. Our studies provide compelling evidence that downregulation of NKD2 expression and alterations in associated regulated pathways have a significant role in driving OS tumor growth and metastasis.

Epithelial tumor cells that have undergone epithelial-to-mesenchymal transition (EMT) are typically prone to metastasis and drug resistance and contribute to a poor clinical outcome. The transcription factor ZEB1 is a known driver of EMT, and mediators of ZEB1 represent potential therapeutic targets for metastasis suppression. Here, we have shown that phosphatidylinositol 3-kinase-targeted (PI3K-targeted) therapy suppresses metastasis in a mouse model of Kras/Tp53-mutant lung adenocarcinoma that develops metastatic disease due to high expression of ZEB1. In lung adenocarcinoma cells from Kras/Tp53-mutant animals and human lung cancer cell lines, ZEB1 activated PI3K by derepressing miR-200 targets, including amphiregulin (AREG), betacellulin (BTC), and the transcription factor GATA6, which stimulated an EGFR/ERBB2 autocrine loop. Additionally, ZEB1-dependent derepression of the miR-200 and miR-183 target friend of GATA 2 (FOG2) enhanced GATA3-induced expression of the p110α catalytic subunit of PI3K. Knockdown of FOG2, p110α, and RHEB ameliorated invasive and metastatic propensities of tumor cells. Surprisingly, FOG2 was not required for mesenchymal differentiation, suggesting that mesenchymal differentiation and invasion are distinct and separable processes. Together, these results indicate that ZEB1 sensitizes lung adenocarcinoma cells to metastasis suppression by PI3K-targeted therapy and suggest that treatments to selectively modify the metastatic behavior of mesenchymal tumor cells are feasible and may be of clinical value.

Non-small-cell lung cancer (NSCLC) demonstrates remarkable molecular diversity. With the completion of The Cancer Genome Atlas (TCGA), there is opportunity for systematic analyses of the entire TCGA NSCLC cohort, including comparisons and contrasts between different disease subsets. On the basis of multidimensional and comprehensive molecular characterization (including DNA methylation and copy, and RNA and protein expression), 1023 NSCLC cases—519 from TCGA adenocarcinoma (AD) project and 504 from TCGA squamous cell carcinoma (SQCC) project—were classified using a 'cluster-of-clusters' analytic approach. Patterns from TCGA NSCLC subsets were examined in independent external databases, including the PROSPECT (Profiling of Resistance patterns and Oncogenic Signaling Pathways in Evaluation of Cancers of the Thorax) NSCLC data set. Nine genomic subtypes of NSCLC were identified, three within SQCC and six within AD. SQCC subtypes were associated with transcriptional targets of SOX2 or p63. One predominately AD subtype (with a large proportion of SQCC) shared molecular features with neuroendocrine tumors. Two AD subtypes manifested a CpG island methylator phenotype. Three AD subtypes showed high p38 and mTOR pathway activation. AD subtypes associated with low differentiation showed relatively worse prognosis. SQCC subtypes and two of the AD subtypes expressed cancer testis antigen genes, whereas three AD subtypes expressed several immune checkpoint genes including PDL1 and PDL2, corresponding with patterns of greater immune cell infiltration. Subtype associations for several immune-related markers—including PD1, PDL1, CD3 and CD8—were confirmed in the PROSPECT cohort using immunohistochemistry. NSCLC molecular subtypes have therapeutic implications and lend support to a personalized approach to NSCLC management based on molecular characterization.

Abstract Liposarcoma can be an aggressive, debilitating, and fatal malignancy. In this study, we identifed miRNAs associated with the differentiation status of liposarcoma to gain insight into the basis for its progression. miRNA expression profiles determined in human tumors and normal fat specimens identified a dedifferentiated tumor expression signature consisting of 35 miRNAs. Deregulated miRNA expression was confirmed in a second independent sample cohort. The miR-155 was the most overexpressed miRNA and functional investigations assigned an important role in the growth of dedifferentiated liposarcoma cell lines. Transient or stable knockdown of miR-155 retarded tumor cell growth, decreased colony formation, and induced G1-S cell-cycle arrest in vitro and blocked tumor growth in murine xenografts in vivo. We identified casein kinase 1α (CK1α) as a direct target of miR-155 control which enhanced β-catenin signaling and cyclin D1 expression, promoting tumor cell growth. In summary, our results point to important functions for miR-155 and β-catenin signaling in progression of liposarcoma, revealing mechanistic vulnerabilities that might be exploited for both prognostic and therapeutic purposes. Cancer Res; 72(7); 1751–62. ©2012 AACR.

Abstract Over the past 20 years, classifications of kidney cancer have undergone major revisions based on morphological refinements and molecular characterizations. The 2016 WHO classification of renal tumors recognizes more than ten different renal cell carcinoma (RCC) subtypes. Furthermore, the marked inter‐ and intra‐tumor heterogeneity of RCC is now well appreciated. Nevertheless, contemporary multi‐omics studies of RCC, encompassing genomics, transcriptomics, proteomics, and metabolomics, not only highlight apparent diversity but also showcase and underline commonality. Here, we wish to provide an integrated perspective concerning the future ‘functional’ classification of renal cancer by bridging gaps among morphology, biology, multi‐omics, and therapeutics. This review focuses on recent progress and elaborates the potential value of contemporary pan‐omics approaches with a special emphasis on cancer genomics unveiled through next‐generation sequencing technology, and how an integrated multi‐omics approach might impact precision‐based personalized kidney cancer care in the near future. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley &amp; Sons, Ltd.

Myxoid liposarcoma (MLPS), a disease especially of young adults with potential for local recurrence and metastasis, currently lacks solid prognostic factors and therapeutic targets. The authors of this report evaluated the natural history and outcome of patients with MLPS and commonly deregulated protein biomarkers.Medical records were retrospectively reviewed for patients who presented to the authors' institution with localized (n = 207) or metastatic (n = 61) MLPS (1990 to 2010). A tissue microarray of MLPS patient specimens (n = 169) was constructed for immunohistochemical analysis of molecular markers.The 5-year and 10-year disease-specific survival rates among patients with localized disease were 93% and 87%, respectively; male gender, age >45 years, and recurrent tumor predicted poor outcome. The local recurrence rate was 7.4%, and the risk of local recurrence was associated with recurrent tumors and nonextremity disease location. Male gender was the main risk factor for metastatic disease, which occurred in 13% of patients. Forty percent of patients who had localized disease received chemotherapy, mostly in the neoadjuvant setting. Immunohistochemical analysis revealed significantly higher expression of C-X-C chemokine receptor type 4 (CXCR4) and platelet-derived growth factor beta (PDGFR-β) in metastatic lesions versus localized lesions. Tumors with a round cell phenotype expressed increased levels of CXCR4, p53, adipophilin, PDGFR-α, PDGFR-β, and vascular endothelial growth factor relative to myxoid phenotype. Only the receptor tyrosine kinase encoded by the AXL gene (AXL) was identified as a prognosticator of disease-specific survival in univariate analysis.In this study, the authors identified clinical and molecular outcome prognosticators for patients with MLPS as well as several potential therapeutic targets.

Background: While epithelial–mesenchymal transition (EMT) can be readily induced experimentally in cancer cells, the EMT process as manifested in human tumors needs to be better understood. Pan‐cancer genomic datasets from The Cancer Genome Atlas (TCGA), representing over 10,000 patients and 32 distinct cancer types, provide a rich resource for examining correlative patterns involving EMT mediators in the setting of human cancers. Results: Here, we surveyed a 16‐gene signature of canonical EMT markers in TCGA pan‐cancer cohort. The histology or cell‐of‐origin of a tumor sample may align more with mesenchymal or epithelial phenotype, and noncancer as well as cancer cells can contribute to the overall molecular patterns observed within a tumor sample; correlation models involving EMT markers can factor in both of the above variables. EMT‐associated genes appear coordinately expressed across all cancers and within each cancer type surveyed. Gene signatures of immune cells correlate highly with EMT marker expression in tumors. In pan‐cancer analysis, several EMT‐related genes can be significantly associated with worse patient outcome. Conclusions: Gene correlates of EMT phenotype in human tumors could include novel mediators of EMT that might be confirmed experimentally, by which TCGA datasets may serve as a platform for discovery in ongoing studies. Developmental Dynamics 247:555‐564, 2018 . © 2017 Wiley Periodicals, Inc.

Chromophobe renal cell carcinoma (chRCC) typically shows ~7 chromosome losses (1, 2, 6, 10, 13, 17, and 21) and ~31 exonic somatic mutations, yet carries ~5%-10% metastatic incidence. Since extensive chromosomal losses can generate proteotoxic stress and compromise cellular proliferation, it is intriguing how chRCC, a tumor with extensive chromosome losses and a low number of somatic mutations, can develop lethal metastases. Genomic features distinguishing metastatic from nonmetastatic chRCC are unknown. An integrated approach, including whole-genome sequencing (WGS), targeted ultradeep cancer gene sequencing, and chromosome analyses (FACETS, OncoScan, and FISH), was performed on 79 chRCC patients including 38 metastatic (M-chRCC) cases. We demonstrate that TP53 mutations (58%), PTEN mutations (24%), and imbalanced chromosome duplication (ICD, duplication of ≥ 3 chromosomes) (25%) were enriched in M-chRCC. Reconstruction of the subclonal composition of paired primary-metastatic chRCC tumors supports the role of TP53, PTEN, and ICD in metastatic evolution. Finally, the presence of these 3 genomic features in primary tumors of both The Cancer Genome Atlas kidney chromophobe (KICH) (n = 64) and M-chRCC (n = 35) cohorts was associated with worse survival. In summary, our study provides genomic insights into the metastatic progression of chRCC and identifies TP53 mutations, PTEN mutations, and ICD as high-risk features.

Metastatic lung cancer is the leading cause of cancer-associated mortality worldwide, therefore necessitating novel approaches to identify specific genetic drivers for lung cancer progression and metastasis. We recently performed an in vivo gain-of-function genetic screen to identify driver genes of lung cancer metastasis. In the study reported here, we identify TMEM106B as a primary robust driver of lung cancer metastasis. Ectopic expression of TMEM106B could significantly promote the synthesis of enlarged vesicular lysosomes that are laden with elevated levels of active cathepsins. In a TFEB-dependent manner, TMEM106B could modulate the expression of lysosomal genes of the coordinated lysosomal expression and regulation (CLEAR) pathway in lung cancer cells and patient samples. We also demonstrate that TMEM106B-induced lysosomes undergo calcium-dependent exocytosis, thereby releasing active lysosomal cathepsins necessary for TMEM106B-mediated cancer cell invasion and metastasis in vivo, which could be therapeutically prevented by pharmacological inhibition of cathepsins. Further, in TCGA LUAD data sets, 19% of patients show elevated expression of TMEM106B, which predicts for poor disease-free and overall-survival.

Implantation of a blastocyst in the uterus is a multistep process tightly controlled by an intricate regulatory network of interconnected ovarian, uterine, and embryonic factors. Bone morphogenetic protein (BMP) ligands and receptors are expressed in the uterus of pregnant mice, and BMP2 has been shown to be a key regulator of implantation. In this study, we investigated the roles of the BMP type 1 receptor, activin-like kinase 2 (ALK2), during mouse pregnancy by producing mice carrying a conditional ablation of Alk2 in the uterus (Alk2 cKO mice). In the absence of ALK2, embryos demonstrate delayed invasion into the uterine epithelium and stroma, and upon implantation, stromal cells fail to undergo uterine decidualization, resulting in sterility. Mechanistically, microarray analysis revealed that CCAAT/enhancer-binding protein β (Cebpb) expression is suppressed during decidualization in Alk2 cKO females. These findings and the similar phenotypes of Cebpb cKO and Alk2 cKO mice lead to the hypothesis that BMPs act upstream of CEBPB in the stroma to regulate decidualization. To test this hypothesis, we knocked down ALK2 in human uterine stromal cells (hESC) and discovered that ablation of ALK2 alters hESC decidualization and suppresses CEBPB mRNA and protein levels. Chromatin immunoprecipitation (ChIP) analysis of decidualizing hESC confirmed that BMP signaling proteins, SMAD1/5, directly regulate expression of CEBPB by binding a distinct regulatory sequence in the 3' UTR of this gene; CEBPB, in turn, regulates the expression of progesterone receptor (PGR). Our work clarifies the conserved mechanisms through which BMPs regulate peri-implantation in rodents and primates and, for the first time, uncovers a linear pathway of BMP signaling through ALK2 to regulate CEBPB and, subsequently, PGR during decidualization.

Infertility and adverse gynecological outcomes such as preeclampsia and miscarriage represent significant female reproductive health concerns. The spatiotemporal expression of growth factors indicates that they play an important role in pregnancy. The goal of this study is to define the role of the ERBB family of growth factor receptors in endometrial function. Using conditional ablation in mice and siRNA in primary human endometrial stromal cells, we identified the epidermal growth factor receptor (Egfr) to be critical for endometrial function during early pregnancy. While ablation of Her2 or Erbb3 led to only a modest reduction in litter size, mice lacking Egfr expression are severely subfertile. Pregnancy demise occurred shortly after blastocyst implantation due to defects in decidualization including decreased proliferation, cell survival, differentiation and target gene expression. To place Egfr in a genetic regulatory hierarchy, transcriptome analyses was used to compare the gene signatures from mice with conditional ablation of Egfr, wingless-related MMTV integration site 4 (Wnt4) or boneless morphogenic protein 2 (Bmp2); revealing that not only are Bmp2 and Wnt4 key downstream effectors of Egfr, but they also regulate distinct physiological functions. In primary human endometrial stromal cells, marker gene expression, a novel high content image-based approach and phosphokinase array analysis were used to demonstrate that EGFR is a critical regulator of human decidualization. Furthermore, inhibition of EGFR signaling intermediaries WNK1 and AKT1S1, members identified in the kinase array and previously unreported to play a role in the endometrium, also attenuate decidualization. These results demonstrate that EGFR plays an integral role in establishing the cellular context necessary for successful pregnancy via the activation of intricate signaling and transcriptional networks, thereby providing valuable insight into potential therapeutic targets.

Although early detection and treatment of prostate cancer (PCa) improves outcomes, many patients still die of metastatic PCa. Here, we report that metastatic PCa exhibits reduced levels of the microRNAsmiR-101 and miR-27a. These micro-RNAs (miRNAs) negatively regulate cell invasion and inhibit the expression of FOXM1 and CENPF, two master regulators of metastasis in PCa. Interestingly, the repression of FOXM1 and CENPF by these miRNAs occurs through COUP-TFII, a member of the orphan nuclear receptors family. Loss of miR-101 positively correlates with the increase of COUP-TFII-FOXM1-CENPF activity in clinical PCa data sets, implicating clinical relevance of such regulation. Further studies show that COUP-TFII is a critical factor controlling metastatic gene networks to promote PCa metastasis. Most importantly, this miRNA-COUP-TFII-FOXM1-CENPF regulatory axis is also involved in the development of enzalutaminde resistance. Taken together, our findings highlight the contribution of specific miRNAs through the regulation of the COUP-TFII-FOXM1-CENPF cascade in PCa metastasis and drug resistance.

KAP1 (TRIM28) is a transcriptional regulator in embryonic development that controls stem cell self-renewal, chromatin organization, and the DNA damage response, acting as an essential corepressor for KRAB family zinc finger proteins (KRAB-ZNF). To gain insight into the function of this large gene family, we developed an antibody that recognizes the conserved zinc fingers linker region (ZnFL) in multiple KRAB-ZNF. Here, we report that the expression of many KRAB-ZNF along with active SUMOlyated KAP1 is elevated widely in human breast cancers. KAP1 silencing in breast cancer cells reduced proliferation and inhibited the growth and metastasis of tumor xenografts. Conversely, KAP1 overexpression stimulated cell proliferation and tumor growth. In cells where KAP1 was silenced, we identified multiple downregulated genes linked to tumor progression and metastasis, including EREG/epiregulin, PTGS2/COX2, MMP1, MMP2, and CD44, along with downregulation of multiple KRAB-ZNF proteins. KAP1-dependent stabilization of KRAB-ZNF required direct interactions with KAP1. Together, our results show that KAP1-mediated stimulation of multiple KRAB-ZNF contributes to the growth and metastasis of breast cancer.

Cardiac dysfunction is the second leading cause of death in myotonic dystrophy type 1 (DM1), primarily because of arrhythmias and cardiac conduction defects. A screen of more than 500 microRNAs (miRNAs) in a DM1 mouse model identified 54 miRNAs that were differentially expressed in heart. More than 80% exhibited downregulation toward the embryonic expression pattern and showed a DM1-specific response. A total of 20 of 22 miRNAs tested were also significantly downregulated in human DM1 heart tissue. We demonstrate that many of these miRNAs are direct MEF2 transcriptional targets, including miRNAs for which depletion is associated with arrhythmias or fibrosis. MEF2 protein is significantly reduced in both DM1 and mouse model heart samples, and exogenous MEF2C restores normal levels of MEF2 target miRNAs and mRNAs in a DM1 cardiac cell culture model. We conclude that loss of MEF2 in DM1 heart causes pathogenic features through aberrant expression of both miRNA and mRNA targets.

A major clinical hurdle for the management of advanced prostate cancer (PCa) in patients is the resistance of tumors to androgen deprivation therapy (ADT) and their subsequent development into castration-resistant prostate cancer (CRPC). While recent studies have identified potential pathways involved in CRPC development, the drivers of CRPC remain largely undefined. Here we determined that nuclear receptor coactivator 2 (NCoA2, also known as SRC-2), which is frequently amplified or overexpressed in patients with metastatic PCa, mediates development of CRPC. In a murine model, overexpression of NCoA2 in the prostate epithelium resulted in neoplasia and, in combination with Pten deletion, promoted the development of metastasis-prone cancer. Moreover, depletion of NCoA2 in PTEN-deficient mice prevented the development of CRPC. In human androgen-sensitive prostate cancer cells, androgen signaling suppressed NCoA2 expression, and NCoA2 overexpression in murine prostate tumors resulted in hyperactivation of PI3K/AKT and MAPK signaling, promoting tumor malignance. Analysis of PCa patient samples revealed a strong correlation among NCoA2-mediated signaling, disease progression, and PCa recurrence. Taken together, our findings indicate that androgen deprivation induces NCoA2, which in turn mediates activation of PI3K signaling and promotes PCa metastasis and CRPC development. Moreover, these results suggest that the inhibition of NCoA2 has potential for PCa therapy.

Abstract Purpose: Dedifferentiated liposarcoma (DDLPS) is an aggressive malignancy that can recur locally or disseminate even after multidisciplinary care. Genetically amplified and expressed MDM2, often referred to as a “hallmark” of DDLPS, mostly sustains a wild-type p53 genotype, substantiating the MDM2:p53 axis as a potential therapeutic target for DDLPS. Here, we report on the preclinical effects of SAR405838, a novel and highly selective MDM2 small-molecule inhibitor, in both in vitro and in vivo DDLPS models. Experimental Design: The therapeutic effectiveness of SAR405838 was compared with the known MDM2 antagonists Nutlin-3a and MI-219. The effects of MDM2 inhibition were assessed in both in vitro and in vivo. In vitro and in vivo microarray analyses were performed to assess differentially expressed genes induced by SAR405838, as well as the pathways that these modulated genes enriched. Results: SAR405838 effectively stabilized p53 and activated the p53 pathway, resulting in abrogated cellular proliferation, cell-cycle arrest, and apoptosis. Similar results were observed with Nutlin-3a and MI-219; however, significantly higher concentrations were required. In vitro effectiveness of SAR405838 activity was recapitulated in DDLPS xenograft models where significant decreases in tumorigenicity were observed. Microarray analyses revealed genes enriching the p53 signaling pathway as well as genomic stability and DNA damage following SAR405838 treatment. Conclusions: SAR405838 is currently in early-phase clinical trials for a number of malignancies, including sarcoma, and our in vitro and in vivo results support its use as a potential therapeutic strategy for the treatment of DDLPS.

<h2>Summary</h2> RNA sequencing (RNA-seq) detects estrogen receptor alpha gene (<i>ESR1</i>) fusion transcripts in estrogen receptor-positive (ER⁺) breast cancer, but their role in disease pathogenesis remains unclear. We examined multiple <i>ESR1</i> fusions and found that two, both identified in advanced endocrine treatment-resistant disease, encoded stable and functional fusion proteins. In both examples, ESR1-e6>YAP1 and ESR1-e6>PCDH11X, <i>ESR1</i> exons 1–6 were fused in frame to C-terminal sequences from the partner gene. Functional properties include estrogen-independent growth, constitutive expression of ER target genes, and anti-estrogen resistance. Both fusions activate a metastasis-associated transcriptional program, induce cellular motility, and promote the development of lung metastasis. ESR1-e6>YAP1- and ESR1-e6>PCDH11X-induced growth remained sensitive to a CDK4/6 inhibitor, and a patient-derived xenograft (PDX) naturally expressing the ESR1-e6>YAP1 fusion was also responsive. Transcriptionally active <i>ESR1</i> fusions therefore trigger both endocrine therapy resistance and metastatic progression, explaining the association with fatal disease progression, although CDK4/6 inhibitor treatment is predicted to be effective.

Renal medullary carcinoma (RMC) is a highly lethal malignancy that mainly afflicts young individuals of African descent and is resistant to all targeted agents used to treat other renal cell carcinomas. Comprehensive genomic and transcriptomic profiling of untreated primary RMC tissues was performed to elucidate the molecular landscape of these tumors. We found that RMC was characterized by high replication stress and an abundance of focal copy-number alterations associated with activation of the stimulator of the cyclic GMP-AMP synthase interferon genes (cGAS-STING) innate immune pathway. Replication stress conferred a therapeutic vulnerability to drugs targeting DNA-damage repair pathways. Elucidation of these previously unknown RMC hallmarks paves the way to new clinical trials for this rare but highly lethal malignancy.

Prostate inflammation has been suggested as an etiology for benign prostatic hyperplasia (BPH). We show that decreased expression of the androgen receptor (AR) in luminal cells of human BPH specimens correlates with a higher degree of regional prostatic inflammation. However, the cause-and-effect relationship between the two events remains unclear. We investigated specifically whether attenuating AR activity in prostate luminal cells induces inflammation. Disrupting luminal cell AR signaling in mouse models promotes cytokine production cell-autonomously, impairs epithelial barrier function, and induces immune cell infiltration, which further augments local production of cytokines and chemokines including Il-1 and Ccl2. This inflammatory microenvironment promotes AR-independent prostatic epithelial proliferation, which can be abolished by ablating IL-1 signaling or depleting its major cellular source, the macrophages. This study demonstrates that disrupting luminal AR signaling promotes prostate inflammation, which may serve as a mechanism for resistance to androgen-targeted therapy for prostate-related diseases.

Tumor cells gain metastatic capacity through a Golgi phosphoprotein 3–dependent (GOLPH3-dependent) Golgi membrane dispersal process that drives the budding and transport of secretory vesicles. Whether Golgi dispersal underlies the pro-metastatic vesicular trafficking that is associated with epithelial-to-mesenchymal transition (EMT) remains unclear. Here, we have shown that, rather than causing Golgi dispersal, EMT led to the formation of compact Golgi organelles with improved ribbon linking and cisternal stacking. Ectopic expression of the EMT-activating transcription factor ZEB1 stimulated Golgi compaction and relieved microRNA-mediated repression of the Golgi scaffolding protein PAQR11. Depletion of PAQR11 dispersed Golgi organelles and impaired anterograde vesicle transport to the plasma membrane as well as retrograde vesicle tethering to the Golgi. The N-terminal scaffolding domain of PAQR11 was associated with key regulators of Golgi compaction and vesicle transport in pull-down assays and was required to reconstitute Golgi compaction in PAQR11-deficient tumor cells. Finally, high PAQR11 levels were correlated with EMT and shorter survival in human cancers, and PAQR11 was found to be essential for tumor cell migration and metastasis in EMT-driven lung adenocarcinoma models. We conclude that EMT initiates a PAQR11-mediated Golgi compaction process that drives metastasis.

Abstract Non-small cell lung cancer (NSCLC) is the leading cause of cancer-related death worldwide, due in part to the propensity of lung cancer to metastasize. Aberrant epithelial-to-mesenchymal transition (EMT) is a proposed model for the initiation of metastasis. During EMT cell-cell adhesion is reduced allowing cells to dissociate and invade. Of the EMT-associated transcription factors, ZEB1 uniquely promotes NSCLC disease progression. Here we apply two independent screens, BioID and an Epigenome shRNA dropout screen, to define ZEB1 interactors that are critical to metastatic NSCLC. We identify the NuRD complex as a ZEB1 co-repressor and the Rab22 GTPase-activating protein TBC1D2b as a ZEB1/NuRD complex target. We find that TBC1D2b suppresses E-cadherin internalization, thus hindering cancer cell invasion and metastasis.

<h2>Summary</h2> A systematic cataloging of genes affected by genomic rearrangement, using multiple patient cohorts and cancer types, can provide insight into cancer-relevant alterations outside of exomes. By integrative analysis of whole-genome sequencing (predominantly low pass) and gene expression data from 1,448 cancers involving 18 histopathological types in The Cancer Genome Atlas, we identified hundreds of genes for which the nearby presence (within 100 kb) of a somatic structural variant (SV) breakpoint is associated with altered expression. While genomic rearrangements are associated with widespread copy-number alteration (CNA) patterns, approximately 1,100 genes—including overexpressed cancer driver genes (e.g., <i>TERT</i>, <i>ERBB2</i>, <i>CDK12</i>, <i>CDK4</i>) and underexpressed tumor suppressors (e.g., <i>TP53</i>, <i>RB1</i>, <i>PTEN</i>, <i>STK11</i>)—show SV-associated deregulation independent of CNA. SVs associated with the disruption of topologically associated domains, enhancer hijacking, or fusion transcripts are implicated in gene upregulation. For cancer-relevant pathways, SVs considerably expand our understanding of how genes are affected beyond point mutation or CNA.

MAPK4 is an atypical MAPK. Currently, little is known about its physiological function and involvement in diseases, including cancer. A comprehensive analysis of 8887 gene expression profiles in The Cancer Genome Atlas (TCGA) revealed that MAPK4 overexpression correlates with decreased overall survival, with particularly marked survival effects in patients with lung adenocarcinoma, bladder cancer, low-grade glioma, and thyroid carcinoma. Interestingly, human tumor MAPK4 overexpression also correlated with phosphorylation of AKT, 4E-BP1, and p70S6K, independent of the loss of PTEN or mutation of PIK3CA. This led us to examine whether MAPK4 activates the key metabolic, prosurvival, and proliferative kinase AKT and mTORC1 signaling, independent of the canonical PI3K pathway. We found that MAPK4 activated AKT via a novel, concerted mechanism independent of PI3K. Mechanistically, MAPK4 directly bound and activated AKT by phosphorylation of the activation loop at threonine 308. It also activated mTORC2 to phosphorylate AKT at serine 473 for full activation. MAPK4 overexpression induced oncogenic outcomes, including transforming prostate epithelial cells into anchorage-independent growth, and MAPK4 knockdown inhibited cancer cell proliferation, anchorage-independent growth, and xenograft growth. We concluded that MAPK4 can promote cancer by activating the AKT/mTOR signaling pathway and that targeting MAPK4 may provide a novel therapeutic approach for cancer.

Abstract Plasma and tumor caveolin-1 (Cav-1) are linked with disease progression in prostate cancer. Here we report that metabolomic profiling of longitudinal plasmas from a prospective cohort of 491 active surveillance (AS) participants indicates prominent elevations in plasma sphingolipids in AS progressors that, together with plasma Cav-1, yield a prognostic signature for disease progression. Mechanistic studies of the underlying tumor supportive onco-metabolism reveal coordinated activities through which Cav-1 enables rewiring of cancer cell lipid metabolism towards a program of 1) exogenous sphingolipid scavenging independent of cholesterol, 2) increased cancer cell catabolism of sphingomyelins to ceramide derivatives and 3) altered ceramide metabolism that results in increased glycosphingolipid synthesis and efflux of Cav-1-sphingolipid particles containing mitochondrial proteins and lipids. We also demonstrate, using a prostate cancer syngeneic RM-9 mouse model and established cell lines, that this Cav-1-sphingolipid program evidences a metabolic vulnerability that is targetable to induce lethal mitophagy as an anti-tumor therapy.

Metastatic castration-resistant prostate cancer (mCRPC) is a malignant and lethal disease caused by relapse after androgen-deprivation (ADT) therapy. Since enzalutamide is innovated and approved by US FDA as a new treatment option for mCRPC patients, drug resistance for enzalutamide is a critical issue during clinical usage. Although several underlying mechanisms causing enzalutamide resistance were previously identified, most of them revealed that drug resistant cells are still highly addicted to androgen and AR functions. Due to the numerous physical functions of AR in men, innovated AR-independent therapy might alleviate enzalutamide resistance and prevent production of adverse side effects. Here, we have identified that yes-associated protein 1 (YAP1) is overexpressed in enzalutamide-resistant (EnzaR) cells. Furthermore, enzalutamide-induced YAP1 expression is mediated through the function of chicken ovalbumin upstream promoter transcription factor 2 (COUP-TFII) at the transcriptional and the post-transcriptional levels. Functional analyses reveal that YAP1 positively regulates numerous genes related to cancer stemness and lipid metabolism and interacts with COUP-TFII to form a transcriptional complex. More importantly, YAP1 inhibitor attenuates the growth and cancer stemness of EnzaR cells in vitro and in vivo. Finally, YAP1, COUP-TFII, and miR-21 are detected in the extracellular vesicles (EVs) isolated from EnzaR cells and sera of patients. In addition, treatment with EnzaR-EVs induces the abilities of cancer stemness, lipid metabolism and enzalutamide resistance in its parental cells. Taken together, these results suggest that YAP1 might be a crucial factor involved in the development of enzalutamide resistance and can be an alternative therapeutic target in prostate cancer.

Abstract During early pregnancy in the mouse, nidatory estrogen (E2) stimulates endometrial receptivity by activating a network of signaling pathways that is not yet fully characterized. Here, we report that bone morphogenetic proteins (BMPs) control endometrial receptivity via a conserved activin receptor type 2 A (ACVR2A) and SMAD1/5 signaling pathway. Mice were generated to contain single or double conditional deletion of SMAD1/5 and ACVR2A/ACVR2B receptors using progesterone receptor (PR)-cre. Female mice with SMAD1/5 deletion display endometrial defects that result in the development of cystic endometrial glands, a hyperproliferative endometrial epithelium during the window of implantation, and impaired apicobasal transformation that prevents embryo implantation and leads to infertility. Analysis of Acvr2a- PRcre and Acvr2b -PRcre pregnant mice determined that BMP signaling occurs via ACVR2A and that ACVR2B is dispensable during embryo implantation. Therefore, BMPs signal through a conserved endometrial ACVR2A/SMAD1/5 pathway that promotes endometrial receptivity during embryo implantation.