Search

Supplementary Figure Legends

FIGURES. Fig S1, Bioinformatics reveal potential CXCL14-receptor candidates; Fig. S2, Expression of CXCL14 correlates with EMT in breast cancer patients; Fig. S3, Expression of CXCL14 correlates with EMT in breast, prostate and ovarian cancer patients; Fig. S4, Enhanced expression of CXCL14 and ACKR2 correlates with EMT in ovarian and prostate cancer patients; Fig. S5, CXCL14-fibroblasts stimulate the EMT program in xenograft tumors; Fig. S6, Soluble factors secreted by CXCL14-fibroblasts induce EMT of breast cancer cells in vitro; Fig. S7, CXCL14-fibroblasts enhance invasion of breast cancer cells; Fig. S8, CXCL14-responsive and non-responsive breast cancer cell lines; Fig. S9, CXCL14 signals through a Gαi-coupled receptor; Fig. S10, CXCL14 responses depend on ACKR2; Fig. S11, Downregulation of ACKR2 abolishes CXCL14-induced ERK signaling; Fig. S12, CXCL14 is not a direct ligand of ACKR2; Fig. S13, Correlation of CXCL14 expression levels and EMT does not reflect stroma abundance; Fig. S14, Enhanced expression of CXCL14 and ACKR2 predicts worse outcome in breast cancer patients; Fig. S15, Identification of potential mediators of CXCL14-fibroblast- stimulated breast cancer EMT; Fig. S16, A CXCL14-ACKR2-NOS1 axis in fibroblasts stimulates EMT and invasion of breast cancer cells // TABLES. Table S1, Primers for the amplification of human and murine genes by qRT-PCR; Table S2, Scored levels of E-Cadherin and CytK8/18 in MCF7 breast xenograft tumors, and CytK8/18 in LnCaP prostate xenograft tumors; Table S3, CXCL14 candidate receptor expression in CXCL14-responsive and non-responsive cell lines; Table S4, Multivariable analysis for CXCL14high/ACKR2high subgroup of breast cancer patients from the TCGA gene expression dataset; Table S5, Univariable cox-regression analyses for the CXCL14 (high)/ACKR2(high) subgroup (90-percentile) in publicly available gene expression datasets (TCGA) from 13 different tumor types.

Purpose:Fibroblasts expressing the orphan chemokine CXCL14 have been previously shown to associate with poor breast cancer prognosis and promote cancer growth. This study explores the mechanism underlying the poor survival associations of stromal CXCL14.Experimental Design:Tumor cell epithelial-to-mesenchymal transition (EMT), invasion, and metastasis were studied in in vitro and in vivo models together with fibroblasts overexpressing CXCL14. An approach for CXCL14 receptor identification included loss-of-function studies followed by molecular and functional endpoints. The clinical relevance was further explored in publicly available gene expression datasets.Results:CXCL14 fibroblasts stimulated breast cancer EMT, migration, and invasion in breast cancer cells and in a xenograft model. Furthermore, tumor cells primed by CXCL14 fibroblasts displayed enhanced lung colonization after tail-vein injection. By loss-of function experiments, the atypical G-protein–coupled receptor ACKR2 was identified to mediate CXCL14-stimulated responses. Downregulation of ACKR2, or CXCL14-induced NOS1, attenuated the pro-EMT and migratory capacity. CXCL14/ACKR2 expression correlated with EMT and survival in gene expression datasets.Conclusions:Collectively, the findings imply an autocrine fibroblast CXCL14/ACKR2 pathway as a clinically relevant stimulator of EMT, tumor cell invasion, and metastasis. The study also identifies ACKR2 as a novel mediator for CXCL14 function and thereby defines a pathway with drug target potential.See related commentary by Zhang et al., p. 3476

PDF file - 135KB, Recruitment of blood vessels into fibroblast-stimulated prostate and breast tumors.

PDF file - 122KB, Lymphatic vessel recruitment into fibroblast-stimulated prostate and breast tumors.

PDF file - 46KB, Nos1 knockdown in CXCL14-fibroblasts.

PDF file - 61KB, CXCL14-stimulation of fibroblasts induces Hmox1 and Nos1.

PDF file - 53KB, Enhanced activation of Hif1α-signaling in CXCL14-fibroblasts.

PDF file - 103KB, Nos1 is dispensable for the growth of CXCL14 fibroblasts under conditions of oxidative stress but important for their paracrine-mediated effects.

PDF file - 136KB

PDF file - 71KB, Expression of CAF marker by engineered fibroblasts.

PDF file - 104KB

PDF file - 42KB, Nos1 is induced upon oxidative stress in NIH3T3 fibroblasts.

PDF file - 59KB, NO and glutathione protect CXCL14-fibroblasts.

Increased expression of GLI1, the main Hedgehog signalling pathway effector, is related to unfavourable prognosis and progressive disease of certain breast cancer subtypes. We used conditional transgenic mice induced to overexpress GLI1 in the mammary epithelium either alone or in combination with deletion of one Trp53 allele to address the role of elevated GLI1 expression in breast tumour initiation and progression. Induced GLI1 expression facilitates mammary gland tumour formation and this was further increased upon heterozygous deletion of Trp53. The GLI1-induced primary tumours were of different murine molecular subtypes, including Normal-likeEx , Class8Ex , Claudin-LowEx and Erbb2-likeEx . The gene expression profiles of some of the tumours correlated well with the PAM50 subtypes for human breast cancer. Whole-exome sequencing revealed somatic mutation profiles with only little overlap between the primary tumours. Orthotopically serially transplanted GLI1-induced tumours maintained the main morphological characteristics of the primary tumours for ≥10 generations. Independent of Trp53 status and molecular subtype, the serially transplanted GLI1-induced tumours were able to grow both in the absence of transgenic GLI1 expression and in the presence of the GLI1 inhibitor GANT61. These data suggest that elevated GLI1 expression has a determinant role in tumour initiation; however, additional genetic events are required for tumour progression.

Cancer-associated fibroblasts (CAFs) regulate tumour growth, metastasis and response to treatment. Recent studies indicate the existence of functionally distinct CAF subsets. Suggested mechanisms whereby CAFs can impact on treatment response include paracrine signalling affecting cancer cell drug sensitivity and effects on tumour drug uptake. PDGFRβ is an important regulator of fibroblasts. Experimental studies have linked PDGFRβ-positive fibroblasts to metastasis and also to reduced tumour drug uptake. This study has investigated the potential role of PDGFRβ-positive fibroblasts in response to adjuvant tamoxifen treatment of breast cancer. Analyses of two breast cancer collections from randomised studies analysing adjuvant tamoxifen treatment in early breast cancer demonstrated significant benefit of tamoxifen in the group with low stromal PDGFRβ, which was not observed in the group with high stromal PDGFRβ. In general terms these findings provide novel evidence, derived from analyses of randomised clinical studies, of response-predictive capacity of a marker-defined subset of CAFs and, more specifically, identify stromal PDGFRβ as a marker related to tamoxifen benefit in early breast cancer.

microRNA (miRNA) dysregulation is a common feature of cancer cells, but the complex roles of miRNAs in cancer are not fully elucidated. Here, we used functional genomics to identify oncogenic miRNAs in non-small cell lung cancer and evaluate their impact on response to epidermal growth factor (EGFR)-targeting therapy. Our data demonstrate that miRNAs with an AAGUGC motif in their seed sequence increase both cancer cell proliferation and sensitivity to EGFR inhibitors. Global transcriptomics, proteomics and target prediction resulted in the identification of several tumor suppressors involved in the G1/S transition as AAGUGC-miRNA targets. The clinical implications of our findings were evaluated by analysis of AAGUGC-miRNA expression in multiple cancer types, supporting the link between this miRNA seed family, their tumor suppressor targets and cancer cell proliferation. In conclusion, we propose the AAGUGC seed motif as an oncomotif and that oncomotif-miRNAs promote cancer cell proliferation. These findings have potential therapeutic implications, especially in selecting patients for EGFR-targeting therapy.

Subcellular localization is a main determinant of protein function; however, a global view of cellular proteome organization remains relatively unexplored. We have developed a robust mass spectrometry-based analysis pipeline to generate a proteome-wide view of subcellular localization for proteins mapping to 12,418 individual genes across five cell lines. Based on more than 83,000 unique classifications and correlation profiling, we investigate the effect of alternative splicing and protein domains on localization, complex member co-localization, cell-type-specific localization, as well as protein relocalization after growth factor inhibition. Our analysis provides information about the cellular architecture and complexity of the spatial organization of the proteome; we show that the majority of proteins have a single main subcellular location, that alternative splicing rarely affects subcellular location, and that cell types are best distinguished by expression of proteins exposed to the surrounding environment. The resource is freely accessible via www.subcellbarcode.org.

Purpose: Fibroblasts expressing the orphan chemokine CXCL14 have been previously shown to associate with poor breast cancer prognosis and promote cancer growth. This study explores the mechanism underlying the poor survival associations of stromal CXCL14. Experimental Design: Tumor cell epithelial-to-mesenchymal transition (EMT), invasion, and metastasis were studied in in vitro and in vivo models together with fibroblasts overexpressing CXCL14. An approach for CXCL14 receptor identification included loss-of-function studies followed by molecular and functional endpoints. The clinical relevance was further explored in publicly available gene expression datasets. Results: CXCL14 fibroblasts stimulated breast cancer EMT, migration, and invasion in breast cancer cells and in a xenograft model. Furthermore, tumor cells primed by CXCL14 fibroblasts displayed enhanced lung colonization after tail-vein injection. By loss-of function experiments, the atypical G-protein–coupled receptor ACKR2 was identified to mediate CXCL14-stimulated responses. Downregulation of ACKR2, or CXCL14-induced NOS1, attenuated the pro-EMT and migratory capacity. CXCL14/ACKR2 expression correlated with EMT and survival in gene expression datasets. Conclusions: Collectively, the findings imply an autocrine fibroblast CXCL14/ACKR2 pathway as a clinically relevant stimulator of EMT, tumor cell invasion, and metastasis. The study also identifies ACKR2 as a novel mediator for CXCL14 function and thereby defines a pathway with drug target potential. See related commentary by Zhang et al., p. 3476

In this study, we describe a novel gene expression signature of platelet-derived growth factor (PDGF)-activated fibroblasts, which is able to identify breast cancers with a PDGF-stimulated fibroblast stroma and displays an independent and strong prognostic significance. Global gene expression was compared between PDGF-stimulated human fibroblasts and cultured resting fibroblasts. The most differentially expressed genes were reduced to a gene expression signature of 113 genes. The biological significance and prognostic capacity of this signature were investigated using four independent clinical breast cancer data sets. Concomitant high expression of PDGFβ receptor and its cognate ligands is associated with a high PDGF signature score. This supports the notion that the signature detects tumors with PDGF-activated stroma. Subsequent analyses indicated significant associations between high PDGF signature score and clinical characteristics, including human epidermal growth factor receptor 2 positivity, estrogen receptor negativity, high tumor grade, and large tumor size. A high PDGF signature score is associated with shorter survival in univariate analysis. Furthermore, the high PDGF signature score acts as a significant marker of poor prognosis in multivariate survival analyses, including classic prognostic markers, Ki-67 status, a proliferation gene signature, or other recently described stroma-derived gene expression signatures.

A role for Hedgehog (Hh) signalling in the development of colorectal cancer (CRC) has been proposed. In CRC and other solid tumours, Hh ligands are upregulated; however, a specific Hh antagonist provided no benefit in a clinical trial. Here we use Hh reporter mice to show that downstream Hh activity is unexpectedly diminished in a mouse model of colitis-associated colon cancer, and that downstream Hh signalling is restricted to the stroma. Functionally, stroma-specific Hh activation in mice markedly reduces the tumour load and blocks progression of advanced neoplasms, partly via the modulation of BMP signalling and restriction of the colonic stem cell signature. By contrast, attenuated Hh signalling accelerates colonic tumourigenesis. In human CRC, downstream Hh activity is similarly reduced and canonical Hh signalling remains predominantly paracrine. Our results suggest that diminished downstream Hh signalling enhances CRC development, and that stromal Hh activation can act as a colonic tumour suppressor., The Hedgehog signalling pathway can drive tumorigenesis. Here, the authors show that in a colitis-associated colon cancer model downstream Hedgehog signalling is restricted to the stroma and its over-activation can inhibit tumorigenesis, associated with activation of BMP signaling.

FunCoup (http://FunCoup.sbc.su.se) is a database that maintains and visualizes global gene/protein networks of functional coupling that have been constructed by Bayesian integration of diverse high-throughput data. FunCoup achieves high coverage by orthology-based integration of data sources from different model organisms and from different platforms. We here present release 2.0 in which the data sources have been updated and the methodology has been refined. It contains a new data type Genetic Interaction, and three new species: chicken, dog and zebra fish. As FunCoup extensively transfers functional coupling information between species, the new input datasets have considerably improved both coverage and quality of the networks. The number of high-confidence network links has increased dramatically. For instance, the human network has more than eight times as many links above confidence 0.5 as the previous release. FunCoup provides facilities for analysing the conservation of subnetworks in multiple species. We here explain how to do comparative interactomics on the FunCoup website. AuthorCount; 6

To identify new regulators of homologous recombination repair, we carried out a genome-wide short-interfering RNA screen combined with ionizing irradiation using RAD51 foci formation as readout. All candidates were confirmed by independent short-interfering RNAs and validated in secondary assays like recombination repair activity and RPA foci formation. Network analysis of the top modifiers identified gene clusters involved in recombination repair as well as components of the ribosome, the proteasome and the spliceosome, which are known to be required for effective DNA repair. We identified and characterized the RNA polymerase II-associated protein CDC73/Parafibromin as a new player in recombination repair and show that it is critical for genomic stability. CDC73 interacts with components of the SCF/Cullin and INO80/NuA4 chromatin-remodeling complexes to promote Histone ubiquitination. Our findings indicate that CDC73 is involved in local chromatin decondensation at sites of DNA damage to promote DNA repair. This function of CDC73 is related to but independent of its role in transcriptional elongation.

OncoTarget: open access impact journal 7(3), 2837-2854 (2015). doi:10.18632/oncotarget.6693, Published by Impact Journals LLC, [S.l.]

Network analysis is an important tool for functional annotation of genes and proteins. A common approach to discern structure in a global network is to infer network clusters, or modules, and assume a functional coherence within each module, which may represent a complex or a pathway. It is however not trivial to define optimal modules. Although many methods have been proposed, it is unclear which methods perform best in general. It seems that most methods produce far from optimal results but in different ways. MGclus is a new algorithm designed to detect modules with a strongly interconnected neighborhood in large scale biological interaction networks. In our benchmarks we found MGclus to outperform other methods when applied to random graphs with varying degree of noise, and to perform equally or better when applied to biological protein interaction networks. MGclus is implemented in Java and utilizes the JGraphT graph library. It has an easy to use command-line interface and is available for download from http://sonnhammer.sbc.su.se/download/software/MGclus/.

MOTIVATION: Analyzing groups of functionally coupled genes or proteins in the context of global interaction networks has become an important aspect of bioinformatic investigations. Assessing the statistical significance of crosstalk enrichment between or within groups of genes can be a valuable tool for functional annotation of experimental gene sets. RESULTS: Here we present CrossTalkZ, a statistical method and software to assess the significance of crosstalk enrichment between pairs of gene or protein groups in large biological networks. We demonstrate that the standard z-score is generally an appropriate and unbiased statistic. We further evaluate the ability of four different methods to reliably recover crosstalk within known biological pathways. We conclude that the methods preserving the second-order topological network properties perform best. Finally, we show how CrossTalkZ can be used to annotate experimental gene sets using known pathway annotations and that its performance at this task is superior to gene enrichment analysis (GEA). AVAILABILITY AND IMPLEMENTATION: CrossTalkZ (available at http://sonnhammer.sbc.su.se/download/software/CrossTalkZ/) is implemented in C++, easy to use, fast, accepts various input file formats, and produces a number of statistics. These include z-score, p-value, false discovery rate, and a test of normality for the null distributions.

We have investigated the distribution of PDGFR-a and -b in the stroma of normal breast tissue, ductal carcinoma in situ (DCIS) and invasive ductal carcinoma (IDC) by immunohistochemistry. Both receptors were expressed by fibroblasts of normal breast tissue. While the PDGFRb remained high during cancer progression, the PDGFRa was lost in most of the DCIS samples and not detectable in IDC. The patient group with PDGFRa(low)/PDGFRb(high) fibroblasts had the highest risk for local recurrences or metastasis. The characteristic loss of stromal PDGFRa expression was also observed in the PyMT mouse breast cancer model. In vitro co-culture experiments confirmed a tumor cell-induced differentiation of fibroblasts, resulting in a fibroblast population marked by PDGFRa(low)/PDGFRb(high) expression. We could identify the interaction of Jagged-1 expressed by the tumor cells with Notch-2 expressed on the fibroblasts, as an essential activator of TGFb signaling in fibroblasts. Gene expression analysis and reporter-gene assays of sorted cells from co-cultures revealed, that fibroblasts switch from a more anti-tumorigenic signaling (expression of BMPs and Wnt antagonists) towards a pro-tumorigenic signaling (expression of BMP inhibitors, TGFb and Wnt ligands). Co-cultured tumor cells show an increased migratory activity and the induction of a mesenchymal phenotype. Inhibition of Notch or TGF signaling in fibroblasts could block these phenotypes. To investigate the impact of fibroblasts on the progression of DCIS to IDC, a DCIS xenograft model was used. Co-injection of MCF10.DCIS cells with human breast fibroblasts, expressing both PDGFRs, into the mouse mammary fat pad resulted in an earlier tumor take with faster tumor growth. The corresponding co-injection experiments with fibroblast and tumor cells, that were depleted of Notch2 respectively Jagged-1 expression by CRISPR/CAS technology is currently ongoing. Collectively the study identifies a potentially druggable paracrine cross-talk between epithelial and stromal cells, which promotes DCIS progression and involves a Notch mediated up-regulation of TGFb signaling. The study also suggests that loss of PDGFRa marks a subset of fibroblasts with specific tumor-promoting activities. Citation Format: Carina Strell, Janna Paulsson, Artur Mezheyeuski, Shao-Bo Jin, Pernilla Roswall, Johanna Andrea, Oliver Frings, Fredrick Warnberg, Hemming Johansson, Karin Jirstrom, Kristian Pietras, Per Hall, Urban Lendahl, Christer Betsholtz, Derek C. Radisky, Arne Ostman. Cell-contact dependent epithelial-stromal cross talk, including a modulation of stromal PDGFR expression, drives the progression of early-stage breast cancer lesions. [abstract]. In: Proceedings of the AACR Special Conference: Function of Tumor Microenvironment in Cancer Progression; 2016 Jan 7–10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2016;76(15 Suppl):Abstract nr B41.

microRNA dysregulation is a common feature of cancer cells, but the complex roles of microRNAs in cancer are not fully elucidated. Here we used functional genomics to identify oncogenic microRNAs in non-small cell lung cancer (NSCLC) and to evaluate their impact on response to EGFR targeting therapy. Our data demonstrates that microRNAs with an AAGUGC-motif in their seed-sequence increase both cancer cell proliferation and sensitivity to EGFR inhibitors. Using miR-372 as a prototypic AAGUGC-miRNA we set out to identify target mRNAs and an explanation to the discovered phenotype. Global transcriptomics, proteomics and target prediction, resulted in the identification of more than 500 candidate target mRNAs, including several tumor suppressors involved in the G1/S transition. The clinical relevance of our findings were evaluated by analysis of public domain data of NSCLC patients, revealing that almost 200 of the candidate target mRNAs showed a negative correlation to AAGUGC-miRNAs. In addition, the analysis of NSCLC clinical data supported the link between this microRNA seed-family and cancer cell proliferation and indicated that high expression of this type of microRNAs is associated with shorter relapse free survival. Expanding the analysis to include additional types of cancer revealed large differences in expression of AAGUGC-miRNA both within and between cancer types. In general the results from the analysis of NSCLC was replicated in other cancer types, where expression of AAGUGC-miRNAs was associated with decreased mRNA levels of tumor suppressor targets and increased expression of genes involved in proliferation. In conclusion we propose that AAGUGC-microRNAs are an integral part of an oncogenic signaling network, and that these findings have potential therapeutic implications, especially in selecting patients for EGFR-targeting therapy. Citation Format: Yan Zhou, Oliver Frings, Erik Fredlund, Jorrit Boekel, Rui M. Branca, Lukas M. Orre. AAGUGC-microRNAs are an integral part of an oncogenic signaling network driving non-small cell lung cancer proliferation. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 1109.

Gene expression analysis is often used to investigate the molecular and functional underpinnings of a phenotype. However, differential expression of individual genes is limited in that it does not consider how the genes interact with each other in networks. To address this shortcoming we propose a number of network-based analyses that give additional functional insights into the studied process. These were applied to a dataset of sex-specific gene expression in the chicken gonad and brain at different developmental stages. We first constructed a global chicken interaction network. Combining the network with the expression data showed that most sex-biased genes tend to have lower network connectivity, that is, act within local network environments, although some interesting exceptions were found. Genes of the same sex bias were generally more strongly connected with each other than expected. We further studied the fates of duplicated sex-biased genes and found that there is a significant trend to keep the same pattern of sex bias after duplication. We also identified sex-biased modules in the network, which reveal pathways or complexes involved in sex-specific processes. Altogether, this work integrates evolutionary genomics with systems biology in a novel way, offering new insights into the modular nature of sex-biased genes.

The InParanoid project gathers proteomes of completely sequenced eukaryotic species plus Escherichia coli and calculates pairwise ortholog relationships among them. The new release 7.0 of the database has grown by an order of magnitude over the previous version and now includes 100 species and their collective 1.3 million proteins organized into 42.7 million pairwise ortholog groups. The InParanoid algorithm itself has been revised and is now both more specific and sensitive. Based on results from our recent benchmarking of low-complexity filters in homology assignment, a two-pass BLAST approach was developed that makes use of high-precision compositional score matrix adjustment, but avoids the alignment truncation that sometimes follows. We have also updated the InParanoid web site (http://InParanoid.sbc.su.se). Several features have been added, the response times have been improved and the site now sports a new, clearer look. As the number of ortholog databases has grown, it has become difficult to compare among these resources due to a lack of standardized source data and incompatible representations of ortholog relationships. To facilitate data exchange and comparisons among ortholog databases, we have developed and are making available two XML schemas: SeqXML for the input sequences and OrthoXML for the output ortholog clusters.

In the growing field of genomics, multiple alignment programs are confronted with ever increasing amounts of data. To address this growing issue we have dramatically improved the running time and memory requirement of Kalign, while maintaining its high alignment accuracy. Kalign version 2 also supports nucleotide alignment, and a newly introduced extension allows for external sequence annotation to be included into the alignment procedure. We demonstrate that Kalign2 is exceptionally fast and memory-efficient, permitting accurate alignment of very large numbers of sequences. The accuracy of Kalign2 compares well to the best methods in the case of protein alignments while its accuracy on nucleotide alignments is generally superior. In addition, we demonstrate the potential of using known or predicted sequence annotation to improve the alignment accuracy. Kalign2 is freely available for download from the Kalign web site (http://msa.sbc.su.se/).