Search

Tumors caused by loss-of-function mutations in genes encoding TCA cycle enzymes have been recently discovered and are now of great interest. Mutations in succinate dehydrogenase (SDH) subunits cause pheochromocytoma/paraganglioma (PCPG) and syndromically associated tumors, which differ phenotypically and clinically from more common SDH-intact tumors of the same types. Consequences of SDH deficiency include rewired metabolism, pseudohypoxic signaling and altered redox balance. PCPG with SDHB mutations are particularly aggressive, and development of treatments has been hampered by lack of valid experimental models. Attempts to develop mouse models have been unsuccessful. Using a new strategy, we developed a xenograft and cell line model of SDH-deficient pheochromocytoma from rats with a heterozygous germline Sdhb mutation. The genome, transcriptome and metabolome of this model, called RS0, closely resemble those of SDHB-mutated human PCPGs, making it the most valid model now available. Strategies employed to develop RS0 may be broadly applicable to other SDH-deficient tumors.

Simple Summary Triple-negative breast cancer (TNBC) has a poor outcome compared to the other major breast cancer subtypes and new therapies are needed. We sought to clarify the functions of a ubiquitous protein, Annexin A1, in the development and progression of TNBC. We found that Annexin A1 expression correlated with poor patient prognosis in basal-like breast tumors and also in the basal like-2 subset of TNBCs. Stable knockdown of Annexin A1 attenuated the growth of SUM149 xenografts, which model basal-like 2 tumors. In a polyoma middle T antigen-driven allograft model of breast cancer, Annexin A1 depletion markedly delayed tumor formation, induced epithelial to mesenchymal transition and upregulated basal markers. Finally, loss of Annexin A1 resulted in the loss of a discrete CD24+/Sca1− population containing putative tumor-initiating cells. Collectively, our data demonstrate a novel cell-autonomous role for Annexin A1 in the promotion of tumor-forming capacity in certain TNBC tumors. Abstract Triple-negative breast cancer (TNBC) has a poor outcome compared to other breast cancer subtypes, and new therapies that target the molecular alterations driving tumor progression are needed. Annexin A1 is an abundant multi-functional Ca2+ binding and membrane-associated protein. Reported roles of Annexin A1 in breast cancer progression and metastasis are contradictory. Here, we sought to clarify the functions of Annexin A1 in the development and progression of TNBC. The association of Annexin A1 expression with patient prognosis in subtypes of TNBC was examined. Annexin A1 was stably knocked down in a panel of human and murine TNBC cell lines with high endogenous Annexin A1 expression that were then evaluated for orthotopic growth and spontaneous metastasis in vivo and for alterations in cell morphology in vitro. The impact of Annexin A1 knockdown on the expression of genes involved in mammary epithelial cell differentia tion and epithelial to mesenchymal transition was also determined. Annexin A1 mRNA levels correlated with poor patient prognosis in basal-like breast tumors and also in the basal-like 2 subset of TNBCs. Unexpectedly, loss of Annexin A1 expression had no effect on either primary tumor growth or spontaneous metastasis of MDA-MB-231_HM xenografts, but abrogated the growth rate of SUM149 orthotopic tumors. In an MMTV-PyMT driven allograft model of breast cancer, Annexin A1 depletion markedly delayed tumor formation in both immuno-competent and immuno-deficient mice and induced epithelial to mesenchymal transition and upregulation of basal markers. Finally, loss of Annexin A1 resulted in the loss of a discrete CD24+/Sca1− population containing putative tumor initiating cells. Collectively, our data demonstrate a novel cell-autonomous role for Annexin A1 in the promotion of tumor-forming capacity in a model of human breast cancer and suggest that some basal-like TNBCs may require high endogenous tumor cell Annexin A1 expression for continued growth.

Background Cancer of unknown primary (CUP), representing approximately 3-5% of all malignancies, is defined as metastatic cancer where a primary site of origin cannot be found despite a standard diagnostic workup. Because knowledge of a patient's primary cancer remains fundamental to their treatment, CUP patients are significantly disadvantaged and most have a poor survival outcome. Developing robust and accessible diagnostic methods for resolving cancer tissue of origin, therefore, has significant value for CUP patients. Methods We developed an RNA-based classifier called CUP-AI-Dx that utilizes a 1D Inception convolutional neural network (1D-Inception) model to infer a tumor's primary tissue of origin. CUP-AI-Dx was trained using the transcriptional profiles of 18,217 primary tumours representing 32 cancer types from The Cancer Genome Atlas project (TCGA) and International Cancer Genome Consortium (ICGC). Gene expression data was ordered by gene chromosomal coordinates as input to the 1D-CNN model, and the model utilizes multiple convolutional kernels with different configurations simultaneously to improve generality. The model was optimized through extensive hyperparameter tuning, including different max-pooling layers and dropout settings. For 11 tumour types, we also developed a random forest model that can classify the tumour's molecular subtype according to prior TCGA studies. The optimised CUP-AI-Dx tissue of origin classifier was tested on 394 metastatic samples from 11 tumour types from TCGA and 92 formalin-fixed paraffin-embedded (FFPE) samples representing 18 cancer types from two clinical laboratories. The CUP-AI-Dx molecular subtype was also independently tested on independent ovarian and breast cancer microarray datasets Findings CUP-AI-Dx identifies the primary site with an overall top-1-accuracy of 98.54% in cross-validation and 96.70% on a test dataset. When applied to two independent clinical-grade RNA-seq datasets generated from two different institutes from the US and Australia, our model predicted the primary site with a top-1-accuracy of 86.96% and 72.46% respectively. Interpretation The CUP-AI-Dx predicts tumour primary site and molecular subtype with high accuracy and therefore can be used to assist the diagnostic work-up of cancers of unknown primary or uncertain origin using a common and accessible genomics platform. Funding NIH R35 GM133562, NCI P30 CA034196, Victorian Cancer Agency Australia.

Merkel cell carcinomas (MCC) are immunogenic skin cancers associated with viral infection or UV mutagenesis. To study T-cell infiltrates in MCC, we analyzed 58 MCC lesions from 39 patients using multiplex-IHC/immunofluorescence (m-IHC/IF). CD4+ or CD8+ T cells comprised the majority of infiltrating T lymphocytes in most tumors. However, almost half of the tumors harbored prominent CD4/CD8 double-negative (DN) T-cell infiltrates (>20% DN T cells), and in 12% of cases, DN T cells represented the majority of T cells. Flow cytometric analysis of single-cell suspensions from fresh tumors identified DN T cells as predominantly Vδ2− γδ T cells. In the context of γδ T–cell inflammation, these cells expressed PD-1 and LAG3, which is consistent with a suppressed or exhausted phenotype, and CD103, which indicates tissue residency. Furthermore, single-cell RNA sequencing (scRNA-seq) identified a transcriptional profile of γδ T cells suggestive of proinflammatory potential. T-cell receptor (TCR) analysis confirmed clonal expansion of Vδ1 and Vδ3 clonotypes, and functional studies using cloned γδ TCRs demonstrated restriction of these for CD1c and MR1 antigen-presenting molecules. On the basis of a 13-gene γδ T–cell signature derived from scRNA-seq analysis, gene-set enrichment on bulk RNA-seq data showed a positive correlation between enrichment scores and DN T-cell infiltrates. An improved disease-specific survival was evident for patients with high enrichment scores, and complete responses to anti–PD-1/PD-L1 treatment were observed in three of four cases with high enrichment scores. Thus, γδ T–cell infiltration may serve as a prognostic biomarker and should be explored for therapeutic interventions. See related Spotlight on p. 600

BackgroundMetastatic Merkel cell carcinoma (mMCC) is an aggressive neuroendocrine malignancy of the skin with a poor prognosis. Immune checkpoint inhibitors (ICIs) have shown substantial efficacy and favorable safety in clinical trials.MethodsMedical records of patients (pts) with mMCC treated with ICIs from August 2015 to December 2018 at Peter MacCallum Cancer Centre in Australia were analyzed. Response was assessed with serial imaging, the majority with FDG-PET/CT scans. RNA sequencing and immunohistochemistry for PD-L1, CD3 and Merkel cell polyomavirus (MCPyV) on tumor samples was performed.Results23 pts with mMCC were treated with ICIs. A median of 8 cycles (range 1 to 47) were administered, with treatment ongoing in 6 pts. Objective responses (OR) were observed in 14 pts (61%): 10 (44%) complete responses (CR) and 4 (17%) partial responses (PR). Median time to response was 8 weeks (range 6 to 12) and 12-month progression-free survival rate was 39%. Increased OR were seen in pts aged less than 75 (OR 80% vs 46%), no prior history of chemotherapy (OR 64% vs 50%), patients with an immune-related adverse event (OR 100% vs 43%) and in MCPyV-negative tumors (OR 69% vs 43%). Pts with a CR had lower mean metabolic tumor volume on baseline FDG-PET/CT scan (CR: 35.7 mL, no CR: 187.8 mL, p=0.05). There was no correlation between PD-L1 positivity and MCPyV status (p=0.764) or OR (p=0.245). 10 pts received radiation therapy (RT) during ICI: 4 pts started RT concurrently (OR 75%, CR 50%), 3 pts had isolated ICI-resistant lesions successfully treated with RT and 3 pts with multisite progression continued to progress despite RT. Overall, 6 pts (26%) had grade 1–2 immune-related adverse events.ConclusionICIs showed efficacy and safety in mMCC consistent with trial data. Clinical and imaging predictors of response were identified.

Transcriptional control shapes a cell's transcriptome composition, but it is RNA processing that refines its expression. The untranslated regions (UTRs) of mRNA are hotspots for regulatory control. Features in these can impact mRNA stability, localisation and translation. Here we describe how alternative cleavage and polyadenylation can change mRNA fate by changing the length of its 3'UTR.

Triple-negative breast cancer (TNBC) represents 10-20% of all human ductal adenocarcinomas and has a poor prognosis relative to other subtypes, due to the high propensity to develop distant metastases. Hence, new molecular targets for therapeutic intervention are needed for TNBC. We recently conducted a rigorous phenotypic and genomic characterization of four isogenic populations of MDA-MB-231 human triple-negative breast cancer cells that possess a range of intrinsic spontaneous metastatic capacities in vivo, ranging from nonmetastatic (MDA-MB-231_ATCC) to highly metastatic to lung, liver, spleen and spine (MDA-MB-231_HM). Gene expression profiling of primary tumours by RNA-Seq identified the fibroblast growth factor homologous factor, FGF13, as highly upregulated in aggressively metastatic MDA-MB-231_HM tumours. Clinically, higher FGF13 mRNA expression was associated with significantly worse relapse free survival in both luminal A and basal-like human breast cancers but was not associated with other clinical variables and was not upregulated in primary tumours relative to normal mammary gland. Stable FGF13 depletion restricted in vitro colony forming ability in MDA-MB-231_HM TNBC cells but not in oestrogen receptor (ER)-positive MCF-7 or MDA-MB-361 cells. However, despite augmenting MDA-MB-231_HM cell migration and invasion in vitro, FGF13 suppression almost completely blocked the spontaneous metastasis of MDA-MB-231_HM orthotopic xenografts to both lung and liver while having negligible impact on primary tumour growth. Together, these data indicate that FGF13 may represent a therapeutic target for blocking metastatic outgrowth of certain TNBCs. Further evaluation of the roles of individual FGF13 protein isoforms in progression of the different subtypes of breast cancer is warranted.

BACKGROUNDThe program of embryonic development is launched by selective activation of a silent maternal transcriptome. InCaenorhabditis elegans, nuclei of the adult germline are responsible for the synthesis of at least two distinct mRNA populations; those required for housekeeping functions, and those that program the oocyte-to-embryo transition. We mapped this separation by changes to the length-distribution of poly(A)-tails that depend on GLD-2 mediated cytoplasmic polyadenylation and its regulators genome-wide.RESULTSMore than 1000 targets of cytoplasmic polyadenylation were identified by differential polyadenylation. Amongst mRNA with the greatest dependence on GLD-2 were those encoding RNA binding proteins with known roles in spatiotemporal patterning such asmex-5andpos-1. In General, the 3’ UTR of GLD-2 targets were longer, contained cytosine-patches, and were enriched for non-standard polyadenylation-motifs. To identify the deadenylase that initiated transcript silencing, we depleted the known deadenylases in thegld-2(0)mutant background. Only the loss of CCF-1 suppressed the short-tailed phenotype of GLD-2 targets suggesting that in addition to its general role in RNA turnover, this is the major deadenylase for regulatory silencing of maternal mRNA. Analysis of poly(A)-tail length-change in the embryo lacking specific RNA-binding proteins revealed new candidates for asymmetric expression in the first embryonic divisions.CONCLUSIONThe concerted action of RNA binding proteins exquisitely regulates GLD-2 activity in space and time. We present our data as interactive web resources for a model where GLD-2 mediated cytoplasmic polyadenylation regulates target mRNA at each stage of worm germline and early embryonic development.

Differential gene expression analysis may discover a set of genes too large to easily investigate, so a means of ranking genes by biological interest level is desired. p values are frequently abused for this purpose. As an alternative, we propose a method of ranking by confidence bounds on the log fold change, based on the previously developed TREAT test. These confidence bounds provide guaranteed false discovery rate and false coverage-statement rate control. When applied to a breast cancer dataset, the top-ranked genes by Topconfects emphasize markedly different biological processes compared to the top-ranked genes by p value. Electronic supplementary material The online version of this article (10.1186/s13059-019-1674-7) contains supplementary material, which is available to authorized users.

BackgroundA differential gene expression analysis may produce a set of significantly differentially expressed genes that is too large to easily investigate, so that a means of ranking genes by their biological interest level is desirable. The life-sciences have grappled with the abuse of p-values to rank genes for this purpose. As an alternative, a lower confidence bound on the magnitude of Log Fold Change (LFC) could be used to rank genes, but it has been unclear how to reconcile this with the need to perform False Discovery Rate (FDR) correction. The TREAT test of McCarthy and Smyth is a step in this direction, finding genes significantly exceeding a specified LFC threshold. Here we describe the use of test inversion on TREAT to present genes ranked by a confidence bound on the LFC, while still controlling FDR.ResultsTesting the Topconfects R package with simulated gene expression data shows the method outperforming current statistical approaches across a wide range of experiment sizes in the identification of genes with largest LFCs. Applying the method to a TCGA breast cancer data-set shows the method ranks some genes with large LFC higher than would traditional ranking by p-value. Importantly these two ranking methods lead to a different biological emphasis, in terms both of specific highly ranked genes and gene-set enrichment.ConclusionsThe choice of ranking method in differential expression analysis can affect the biological interpretation. The common default of ranking by p-value is implicitly by an effect size in which each gene is standardized to its own variability, rather than comparing genes on a common scale, which may not be appropriate. The Topconfects approach of presenting genes ranked by confident LFC effect size is a variation on the TREAT method with improved usability, removing the need to fine-tune a threshold parameter and removing the temptation to abuse p-values as a de-facto effect size.

Triple-negative breast cancer (TNBC) represents 10-20% of all human ductal adenocarcinomas and has a poor prognosis relative to other subtypes. Hence, new molecular targets for therapeutic intervention are necessary. Analyses of panels of human or mouse cancer lines derived from the same individual that differ in their cellular phenotypes but not in genetic background have been instrumental in defining the molecular players that drive the various hallmarks of cancer. To determine the molecular regulators of metastasis in TNBC, we completed a rigorous in vitro and in vivo characterisation of four populations of the MDA-MB-231 human breast cancer line ranging in aggressiveness from non-metastatic to spontaneously metastatic to lung, liver, spleen and lymph node. Single nucleotide polymorphism (SNP) array analyses and genome-wide mRNA expression profiles of tumour cells isolated from orthotopic mammary xenografts were compared between the four lines to define both cell autonomous pathways and genes associated with metastatic proclivity. Gene set enrichment analysis (GSEA) demonstrated an unexpected association between both ribosome biogenesis and mRNA metabolism and metastatic capacity. Differentially expressed genes or families of related genes were allocated to one of four categories, associated with either metastatic initiation (e.g. CTSC, ENG, BMP2), metastatic virulence (e.g. ADAMTS1, TIE1), metastatic suppression (e.g. CST1, CST2, CST4, CST6, SCNNA1, BMP4) or metastatic avirulence (e.g. CD74). Collectively, this model system based on MDA-MB-231 cells should be useful for the assessment of gene function in the metastatic cascade and also for the testing of novel experimental therapeutics for the treatment of TNBC. This article has an associated First Person interview with the first author of the paper., Summary: Four MDA-MB-231 cell line variants were evaluated for spontaneous metastatic capacity. DNA SNP arrays and RNA-Seq analysis of xenografted tumour cells were conducted to identify genes associated with metastatic proclivity.

Background mMCC is a rare, highly aggressive neuroendocrine cancer of the skin with a poor prognosis. ICIs have favourable efficacy and safety in clinical trials. We outline single centre experience utilising ICIs in mMCC. Methods Medical records of patients (pts) with mMCC treated with ICIs from Aug 2015 to Dec 2018 at Peter MacCallum Cancer Centre in Australia were retrospectively analysed. RNA sequencing and immunohistochemistry for PD-L1, CD3 and Merkel cell polyomavirus (MCPyV) on tumor samples were performed. Baseline tumor volumes and responses were assessed with FDG-PET scans using the Hicks criteria. Results 23 pts with mMCC were treated with ICIs. Pt characteristics are summarised in the table. A median of 8 cycles (range 1 to 47) were administered, with treatment ongoing in 6 pts. Objective responses (OR) were observed in 14 pts (61%); 10 (44%) complete metabolic responses (CMR) and 4 (17%) partial metabolic responses (PMR). Median time to response was 8 weeks (range 6 to 12) and 12-month progression-free survival (PFS) rate was 39%. Increased OR were seen in pts aged less than 75 (OR 8/10, 80% vs 46%), pts with normal baseline LDH (OR 5/7, 71% vs 50%), pts with an immune-related adverse event (irAE) (OR 6/6, 100% vs 47%) and in MCPyV negative pts (OR 11/16, 69% vs 43%). Pts with a CMR had lower mean-tumor volume on baseline FDG-PET scan (CMR: 35.7mL, no CMR: 187.8mL, p = 0.05). No correlation was seen between tumor PD-L1 positivity and response to ICI (p = 0.25) or MCPyV status (p = 0.76). Similarly, no association was seen between OR and CD3 staining within the tumor (p = 0.07). 10 pts received radiation (RT) during ICI: 4 pts started RT concurrently (OR 75%, CMR 50%), 3 pts had isolated ICI-resistant lesions successfully treated with RT and 3 pts with multisite progression continued to progress despite RT. 6 pts (26%) had a Grade 1-2 irAE. Table . 341O Baseline characteristics Patient Characteristics n = 23 Age (years) Median (range) 75 (64 to 91) Gender Male Female 18 5 ECOG 0 12 3 8 12 3 ICI Avelumab Pembrolizumab Other (BGB-A317) 15 5 3 MCPyV status Positive Negative 7 16 Visceral disease 17 Prior chemotherapy 10 History of autoimmune disease 2 Previous immunosuppression 2 Conclusions ICIs showed efficacy and safety consistent with trial data. Younger age, negative MCPyV status, normal LDH, lower baseline FDG-PET tumor volume and irAEs are potentially associated with enhanced response. Legal entity responsible for the study The authors. Funding Has not received any funding. Disclosure R.J. Hicks: Shareholder / Stockholder / Stock options: Telix Pharmaceuticals; Travel / Accommodation / Expenses: Endocyte; Travel / Accommodation / Expenses: GE Medical Systems. R. Tothill: Honoraria (self), Travel / Accommodation / Expenses: Merck Serono. All other authors have declared no conflicts of interest.