Oncotarget

Research Papers:

miR-509-3p is clinically significant and strongly attenuates cellular migration and multi-cellular spheroids in ovarian cancer

PDF |  HTML  |  Supplementary Files  |  Order a Reprint  |  This article has been corrected. Oncotarget. 2017; 8(10):17406.

Oncotarget. 2016; 7:25930-25948. https://doi.org/10.18632/oncotarget.8412

Metrics: HTML 1042 views  |   ?  

Yinghong Pan, Gordon Robertson, Lykke Pedersen, Emilia Lim, Anadulce Hernandez-Herrera, Amy C. Rowat, Sagar L. Patil, Clara K. Chan, Yunfei Wen, Xinna Zhang, Upal Basu-Roy, Alka Mansukhani, Andy Chu, Payal Sipahimalani, Reanne Bowlby, Denise Brooks, Nina Thiessen, Cristian Coarfa, Yussanne Ma, Richard A. Moore, Jacquie E. Schein, Andrew J. Mungall, Jinsong Liu, Chad V. Pecot, Anil K. Sood, Steven J.M. Jones, Marco A. Marra, Preethi H. Gunaratne _

Abstract

Yinghong Pan1,*, Gordon Robertson2,*, Lykke Pedersen3, Emilia Lim2, Anadulce Hernandez-Herrera1, Amy C. Rowat4, Sagar L. Patil1, Clara K. Chan4, Yunfei Wen5, Xinna Zhang5,6, Upal Basu-Roy9, Alka Mansukhani9, Andy Chu2, Payal Sipahimalani2, Reanne Bowlby2, Denise Brooks2, Nina Thiessen2, Cristian Coarfa14, Yussanne Ma2, Richard A. Moore2, Jacquie E. Schein2, Andrew J. Mungall2, Jinsong Liu8, Chad V. Pecot10, Anil K. Sood5,6,7,8, Steven J.M. Jones2,12,13, Marco A. Marra2,13, Preethi H. Gunaratne1,11

1Department of Biochemistry and Biology, University of Houston, Houston, TX, USA

2Canada’s Michael Smith Genome Sciences Centre, BC Cancer Agency, Vancouver, BC, Canada

3Department of Biology, University of Copenhagen, Copenhagen, Denmark

4Department of Integrative Biology and Physiology, University of California, Los Angeles, CA, USA

5Department of Gynecologic Oncology and Reproductive Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX, USA

6Center for RNA Interference and Non-Coding RNA, The University of Texas MD Anderson Cancer Center, Houston, TX, USA

7Department of Cancer Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA

8Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA

9Department of Microbiology, New York University School of Medicine, New York, NY, USA

10UNC Lineberger Comprehensive Cancer Center, Thoracic Medical Oncology, University of North Carolina, Chapel Hill, NC, USA

11Department of Pathology and Immunology, Baylor College of Medicine, Houston, TX, USA

12Department of Molecular Biology and Biochemistry, Simon Fraser University, Burnaby, BC, Canada

13Department of Medical Genetics, University of British Columbia, Vancouver, BC, Canada

14Department of Molecular and Cellular Biology, Baylor College of Medicine, Houston, TX, USA

*These authors contributed equally to this work

Correspondence to:

Preethi H. Gunaratne, email: phgunaratne@uh.edu

Keywords: microRNA 509-3p, ovarian cancer, extracellular matrix (ECM), YAP1, spheroid formation

Received: August 26, 2015     Accepted: March 13, 2016     Published: March 27, 2016

ABSTRACT

Ovarian cancer presents as an aggressive, advanced stage cancer with widespread metastases that depend primarily on multicellular spheroids in the peritoneal fluid. To identify new druggable pathways related to metastatic progression and spheroid formation, we integrated microRNA and mRNA sequencing data from 293 tumors from The Cancer Genome Atlas (TCGA) ovarian cancer cohort. We identified miR-509-3p as a clinically significant microRNA that is more abundant in patients with favorable survival in both the TCGA cohort (P = 2.3E–3), and, by in situ hybridization (ISH), in an independent cohort of 157 tumors (P < 1.0E–3). We found that miR-509-3p attenuated migration and disrupted multi-cellular spheroids in HEYA8, OVCAR8, SKOV3, OVCAR3, OVCAR4 and OVCAR5 cell lines. Consistent with disrupted spheroid formation, in TCGA data miR-509-3p’s most strongly anti-correlated predicted targets were enriched in components of the extracellular matrix (ECM). We validated the Hippo pathway effector YAP1 as a direct miR-509-3p target. We showed that siRNA to YAP1 replicated 90% of miR-509-3p-mediated migration attenuation in OVCAR8, which contained high levels of YAP1 protein, but not in the other cell lines, in which levels of this protein were moderate to low. Our data suggest that the miR-509-3p/YAP1 axis may be a new druggable target in cancers with high YAP1, and we propose that therapeutically targeting the miR-509-3p/YAP1/ECM axis may disrupt early steps in multi-cellular spheroid formation, and so inhibit metastasis in epithelial ovarian cancer and potentially in other cancers.


miR-509-3p is clinically significant and strongly attenuates cellular migration and multi-cellular spheroids in ovarian cancer | Pan | Oncotarget

INTRODUCTION

Epithelial ovarian cancer (EOC) is the most lethal gynecological malignancy, typically presenting as an aggressive, advanced stage cancer with widespread metastases in the peritoneum [1]. It metastasizes primarily by cells exfoliating into the peritoneal fluid, and aggregating into multicellular spheroids that invade the peritoneal membrane [2, 3] and can resist chemotherapy [46]. While patients typically respond initially to cytoreductive surgery and to platinum- and taxane-based chemotherapy, most develop relapsed disease that becomes refractory to all current therapies. There is a pressing need for new therapeutic strategies that sensitize these cancers to current therapies.

Recently, The Cancer Genome Atlas (TCGA) Research Network characterized 489 high-grade serous ovarian cancer (HGSOC) tumors [7], using microarrays for gene and microRNA expression. Here, we report results from 293 TCGA ovarian tumors, using messenger RNA (mRNA) and microRNA (miRNA/miR) sequencing data. miRNAs are ~22-nucleotide non-coding RNAs that regulate transcript stability and translation, primarily by binding sites in 3′-UTRs [8]. We anticipated that the sequencing data could yield new insights by offering a wider dynamic range than microarrays, higher spatial resolution and sensitivity, isoform-specific mRNA expression, and better discrimination between miRNA stem-loops and mature miRNA strands [9].

We identified the X-linked miR-509-3p as a clinically significant miRNA. The miR was more abundant in patients with favorable survival in both the TCGA and an independent cohort. In attenuated migration, invasion, and aggregation into three-dimensional (3D) spheroids in the six ovarian cancer cell lines tested: HEYA8, OVCAR8, SKOV3, OVCAR3, OVCAR4 and OVCAR5. We showed that YAP1, a key oncogene in a subset of ovarian and other cancers, is a direct downstream target of miR-509-3p. We also found that siRNA to YAP1 was necessary and sufficient to replicate ~90% of miR-509-3p-mediated attenuation of migration in OVCAR8, which contained high levels of YAP1 protein. However, siYAP1 had no impact on migration attenuation or multi-cellular spheroids in the other cell lines, suggesting that the miR-509-3p/YAP1 axis is relevant only in cells containing high YAP1.

Since miR-509-3p’s strongest anti-correlated predicted targets in TCGA ovarian tumors were enriched for components of the extracellular matrix (ECM), and effects on ECM components are consistent with effects on 3D spheroid formation, we propose that miR-509-3p may act in ovarian cancers by targeting both ECM components and YAP1. We anticipate that the miR-509-3p/YAP1/ECM axis offers a novel therapeutic focus for controlling metastatic progression in HGSOC.

RESULTS

An independent component that discriminates overall survival is enriched in Xq27.3 miR cluster members, and in ECM and cell adhesion genes

We identified mRNA transcript isoforms from the mRNA-seq data using reference-based assembly [10], and calculated TargetScan 6.0 [11] miRNA binding sites on the full length reconstructed transcripts. For the 293 tumor samples with both mRNA-seq and miRNA-seq data (Table S1), we identified anti-correlated miRNA:mRNA pairs by calculating Spearman correlation coefficients (r) between the abundance profiles for mRNA isoforms and for 5p and 3p miRNA strands, thresholding at a q-value (FDR) < 0.05 (Table S2). Consistent with previous work [12], the miR-9 family and miR-29a/b were important for interactions with correlations below –0.3. Stronger interactions were also enriched in members of a ~100-kb cluster in Xq27.3 that contains hsa-mir-506, 507, 508, 509, 510, 513 and 514 (Figures 1A1C, and S1A), particularly for miR-506, miR-508, miR-509, miR-513 and miR-514. For the 250 target genes with the strongest anti-correlations to members of this cluster, enriched Gene Ontology terms highlighted the ECM, and enriched KEGG pathways included ECM-receptor interactions (FDR = 6E–8) and focal adhesion (FDR = 0.02) (Table S3).

Significant miRNAs, mRNAs and pathways.

Figure 1: Significant miRNAs, mRNAs and pathways. (A) Schematic transcriptome-wide view of the number of significant (FDR < 0.05) anti-correlations per miRNA (Spearman r < –0.3, for miRNAs with at least one significant anti-correlation). miRNA locations are not to scale. (B) The miRNAs in (A), ordered by the number of miR-gene interactions, and thresholded at 20 interactions. (C) miR-gene interactions with Spearman correlation coefficients less than –0.5. (D–G) Independent component analysis. (D) Kaplan-Meier plot for overall survival when samples are segregated by high vs. low component weights in IC 5 (see text). (E) Component weights of IC 5, for samples with favorable vs. unfavorable survival, for the mixing matrices for mRNAs (left) and miRNAs (right). Because samples with favorable survival have larger component weights than samples with favorable vs. unfavorable survival in IC 5, mRNAs with negative gene weights in this IC will be downregulated in samples with favorable survival. (F) Gene weights and q-values for the mRNAs and miRNAs that have significant gene weights in IC 5 (thresholds α = 6 for mRNAs, 1 for miRNAs) and are differentially expressed between samples with favorable and unfavorable overall survival. Green vs. orange: upregulated vs. downregulated in samples with favorable survival. (G) Enriched GO terms and pathways for IC 5 genes (q < 0.001, hypergeometric test, FDR-correction). R and K indicate REACTOME and KEGG, respectively.

We extended these results by using independent component analysis (ICA) [13] to retrieve significant associations between miRNAs, mRNAs, and pathways. ICA is robust to co-expressed genes that are not always co-regulated [14], and to anti-correlations between miRNAs and target mRNAs that are not always significant [15]. Independent components (ICs) have been shown to capture the differential regulation of biological processes and metabolic pathways in breast cancer [16], acute myeloid leukemia [17], and Alzheimer’s disease [18].

Briefly, ICA linearly decomposes an expression matrix (E) into ICs by solving the matrix equation E = CM. The columns of the component matrix C are the ICs in which each gene has a weight, and the weight of each IC in the samples is contained in the rows of the mixing matrix M. For miRNA-mRNA interactions, the gene weights of targeted mRNAs should be inversely correlated with the gene weights of targeting miRNAs in an IC that captures regulation by miRNAs. To identify such anti-correlations we modified the R function fastICA, which is widely used for expression data [19]. The resulting method, miICA, considers mRNA and miRNA data simultaneously.

ICA returned five ICs (Table S4). Using each IC’s component weights, i.e. the rows of Mm and Mmi, we divided the samples into groups with favorable vs. unfavorable survival or recurrence by identifying an optimal component weight threshold for a two-group Kaplan-Meier test (K-M test) against clinical outcomes [20]. IC 5 best discriminated samples by survival (P = 2.3E-3, Figures 1D, S1B and S1C). IC 5′s component weights were higher in samples with favorable survival (Figure 1E, K-W test, FDR correction), which means that mRNAs or miRNAs with positive gene weights in this component were more abundant in the group of samples with favorable survival (Figures 1C and S1D). IC 5 contained 65 mRNAs with very significant gene weights (Figure 1F, α = 6), all of which were also differentially expressed between samples with favorable and unfavorable survival (q < 0.05, K-W test, FDR correction). Since miRNAs generally had lower gene weights than mRNAs in IC 5, we used a lower threshold α cut-off (α = 1) to identify 62 miRNAs with significant gene weights in IC 5, 24 of which were also differentially expressed between samples with favorable and unfavorable survival (Figures 1F and S1D, q < 1.0E-2, K-W test, FDR correction). miRNAs from the Xq27.3 miR cluster were more abundant in samples with favorable survival and had positive gene weights in IC 5; thus, targets of these miRs that are relevant to HGSOC should be sought among mRNAs with negative gene weights in this IC.

Consistent with the correlation analysis, many of the most significantly enriched pathways for IC 5 (Figure 1G) involve genes related to the ECM and cell adhesion (Table S4, threshold α = 3, q < 1.0E-3, K-W test). In all of these pathways the majority of genes had negative gene weights that were anti-correlated with the IC 5 gene weights of members of the Xq27.3 miR cluster. Predicted targets of the cluster miRs were significantly overrepresented among genes with negative gene weights in IC 5 (Table S4), consistent with targets of the miR cluster being associated with the enriched pathways. Overall, these results indicate that members of the Xq27.3 miR cluster, and the significant genes from IC 5 that are annotated as members of significant pathways, may be associated with survival in HGSOC.

The Xq27.3 miRNA cluster is predicted to target structural and regulatory ECM components, and EMT regulators

To characterize in more detail ECM components that may be targeted by miR cluster members, we compared 1368 predicted target genes, and the 541 significant IC 5 genes, to six functional groups that have been defined for core matrisome and ECM-associated genes [21] (Table S5). For the core matrisome, predicted targets included 38 of 200 glycoproteins (P = 5.3E-10), 17 of 45 collagens (P = 6.2E-10) and 14 of 36 proteoglycans (P = 1.3E-8). For ECM-affiliated genes, predicted targets included 23 of 177 affiliated proteins (P = 6.7E-4), 51 of 254 regulators (P = 8.0E-14), and 50 of 353 secreted factors (P = 4.8E-8). Similarly, IC 5 genes were enriched in all ECM functional gene groups (P = 6.6E-24 to 6.6E-8) except ECM-affiliated genes.

We then assessed ECM-rich gene sets [22] that were prognostic for overall survival or recurrence in microarray-based gene expression data from the original TCGA study and two independent studies [7, 23, 24]. We compared six sets of highest-ranked 100 genes from: a) networks based on co-expression, functional linkage, and L2 penalty, and b) survival and recurrence outcomes, to the functional ECM groups, predicted miR cluster targets, and IC 5 genes (Table S5). The six gene sets were typically enriched (P < 0.001) for ECM glycoproteins, collagens, regulators and secreted factors. The 1368 predicted miR cluster targets were enriched in the gene sets for overall survival (P = 4.2E-52 to 1.3E-22) and for recurrence (P = 1.2E-33 to 1.4E-12). Finally, the 541 IC 5 genes were enriched in all six prognostic gene sets (P = 1E-25 to 3E-80).

Our sequencing data also identified epithelial-mesenchymal transition (EMT) regulators [25] as predicted targets of members of the miR cluster and of the miR-200 family: SNAI1 (miR-510-5p), SNAI2 (miR-507, miR-200b-3p), TWIST1 (miR-507, 508-3p, 509-3p), TWIST2 (miR-141-3p, 29b-3p, 29c-3p), ZEB1 (miR-141-3p, 200b-3p) and ZEB2 (miR-200b-3p, 507, 508-5p and 3p, 513c-5p, 514b-5p). Significant IC 5 genes included SNAI2, TWIST1, TWIST2 and ZEB1 (Table S2).

Given these results, we assessed how miR-509-3p mimics influenced levels of transcripts associated with ECM and EMT in TP53-wild type HEYA8 and TP53-mutated OVCAR8 ovarian cancer cell lines, in which we carried out time-series qPCRs on selected genes. TaqMan assays confirmed that miR-509-3p increased ~10,000 fold in abundance at 72, 96 and 120 h following transient transfection of miR-509-3p mimics (Figure S3A). Of the ECM and cytoskeletal genes, ACTA2, COL1A1 and COL3A1 were downregulated in response to miR-509-3p mimics in both cell lines at 72, 96 and 120 h after transfection, which correspond to 0, 24 and 48 h time points of the migration/invasion experiments (Figure S3B to S3D). Treatment with miR-509-3p mimics resulted in COL5A1 downregulation in HEYA8 but upregulation in OVCAR8, while COL5A2 was not significantly affected by the mimic in HEYA8, and was not detected in OVCAR8. ECM genes SPARC, FN1 and PRRX1 were upregulated in response to miR-509-3p mimics in both cell lines. miR-509-3p mimics downregulated TWIST at all three time points in both cell lines, and more significantly in HEYA8 than in OVCAR8 at each time point. SNAI2 was upregulated at all three time points in both cell lines, with fold changes generally more significant in OVCAR8 than in HEYA8.

miR-509-3p is positively associated with clinical outcomes and is localized to tumor cells

For members of the Xq27.3 miR cluster, and the miR-29 and miR-200-families, we identified compact subsets of miRNA 5p and 3p strands whose abundance discriminated overall survival or time to recurrence. After correcting for multiple testing, only miRs from the Xq27.3 cluster remained statistically significant. For example, miR-506, 507, 508 and 509 collectively stratified a group of 37 samples by overall survival (P < 0.049) (Figures 2A, 2B and S2; Table S6).

miR-509-3p and clinical outcomes.

Figure 2: miR-509-3p and clinical outcomes. (A, B) Kaplan-Meier plots for overall survival and time to recurrence for a discriminatory subset of Xq27.3 miRs, from TCGA sequencing data. (C, D) For an independent set of HGSOC samples, Kaplan-Meier plots for overall survival (N = 157) and progression-free time (N = 140) for sample groups separated by the median miR-509-3p in situ hybridization (ISH) score. (E, F) Representative ISH micrographs for miR-509-3p, showing (E) high (291) and (F) low (0) ISH staining intensities. Black and grey arrows indicate tumor and stroma cells respectively. Scale bars are 100 um.

From computational results, and the migration/invasion results described below, we focused on miR-509-3p. To confirm miR-509-3p’s association with overall survival, and to determine whether the miR was expressed by cells within the tumor and/or the stromal cells from the tumor microenvironment, we used in situ hybridization to stain tissue microarrays for miR-509-3p in an independent set of 157 HGSOC samples (Table 1). We scored staining intensity in high-powered image fields in a blinded fashion by quantitative image analysis with CellProfiler v2.0 [26]. Consistent with results from sequencing data, overall survival was favorable (P < 1.0E-3) for samples with higher miR-509-3p expression levels (Figure 2C); progression-free survival was also favorable (P < 1.0E-3, Figure 2D). When miR-509-3p was detected, we found that it was largely localized within the tumor cells and not in the stroma (Figure 2E and 2F).

Table 1: Summary of independent HGSOC cohort used for ISH

Variable

 

miR-509-3p

P-value

Low (N = 75)

High (N = 82)

Tumour stage

I or II

4

9

0.21

 

III or IV

70

73

 

Cytoreduction

Optimal

41

56

0.11

 

Suboptimal

30

24

 

Recurrence or progression

Yes

57

57

3E-2

 

No

8

21

 

Median PFS (yr)

 

7.7

21.3

< 1.0E-3

Median OS (yr)

 

40

81.2

< 1.0E-3

For miR-509-3p, low and high refer to ISH staining intensities that are below and above the median intensity.

However, the ovarian cancer cell lines we tested had relatively low levels of miR-509-3p expression by qPCR (Figure S4A). The microarray-based miRNA abundance data in the NCI-60 cell lines (GEO: GSE26375) [27] includes seven ovarian cell lines (IGROV1, NCI/ADR-RES, OVCAR-3, OVCAR-4, OVCAR-5, OVCAR-8, SK-OV-3), five of which are included in the work that we report here. In these data too, the expression of miR-509-3p is low compared to other microRNAs such as miR-29abc and members of the miR-200 family (Figure S4B and S4C).

miR-508-3p and miR-509-3p attenuate ovarian cancer cell migration and invasion

Because miR-508 and miR-509 were correlated with overall survival, and had predicted targets that were enriched for ECM components and for EMT regulators, we hypothesized that these miRNAs may be able to influence clinical outcomes in HGSOC by attenuating migration and invasion, and so metastatic progression. We assessed this in vitro, using a 96-well scratch-wound assay to compare HEYA8 and OVCAR8 cells that were untreated, versus those treated with either miRNA mimics from the Xq27.3 miR cluster or with a scrambled negative control (SCR) sequence.

Transient transfection with miR-509-3p strongly attenuated the ability of cells to migrate (Figure 3A and 3C) and invade a Matrigel matrix in both cell lines (Figure 3B and 3D). The relative wound density (RWD) of migration at 24 h was 39% lower in HEYA8 (p53-wild type) (P = 5.3E-13) and 64% lower in OVCAR8 (p53-mutant) cells (P = 5.8E-13) than in the same cells treated with the scrambled negative control sequence. Transient transfection with miR-508-3p had a smaller effect on migration for HEYA8 cells (6.3% reduction in RWD at 24 h, P = 0.20), and an intermediate effect in OVCAR8 cells (25% reduction in RWD at 24 h, P = 5.8E-13). In contrast, miR-509-5p only slightly reduced migration, as shown by the 11% decrease in RWD at 24 h for OVCAR8 cells (P = 6.6E-4), and showed no significant effect in HEYA8 (P = 1.00). The ability of cells to invade a Matrigel matrix was also most strongly attenuated by miR-509-3p, with a 70% reduction in RWD for HEYA8 (P = 5.4E-13), and a 29% reduction in OVCAR8 compared to SCR (P = 3.0E-3). miR-508-3p also reduced invasion, as indicated by the 6% decrease in RWD for HEYA8 (P = 1.9E-3) and 18% in OVCAR8 (P = 1.3E-2) and miR-508-5p did not significantly affect invasion for OVCAR8 (P = 0.89) or HEYA8 cells (P = 0.15). In contrast, miR-509-5p reduced invasion in OVCAR8 (14% reduction in RWD at 24 h, P = 3.2E-6), but did not significantly affect migration in HEYA8 cells (P = 1.0).

In vitro 96-well migration and invasion results after transfecting miR mimics into HEYA8 and OVCAR8 cells.

Figure 3: In vitro 96-well migration and invasion results after transfecting miR mimics into HEYA8 and OVCAR8 cells. (A) Migration and (B) invasion for HEYA8 cells (N = 3). (C) Migration and (D) invasion for OVCAR8 cells (N = 3). Error bars show standard deviations. The time courses start 72 h after transfection.

In independent experiments, we confirmed that miR-509-3p mimics attenuated migration strongly in OVCAR8 and SKOV3 and moderately in OVCAR3, OVCAR4, OVCAR5 and HEYA8 (Figures 5A, S5). Since the basal level of miR-509-3p in cells was very low in all six ovarian cancer cell lines tested (Figure S4A), the miR-509-3p inhibitor did not affect migration in these cell lines (Figure 6A).

YAP1 target validation, and correlation with survival and migration.

Figure 4: YAP1 target validation, and correlation with survival and migration. (A) Relative mRNA transcript levels by qPCR and (B) protein levels by Western blot after miR-509-3p transfection and siRNA for YAP1, for HEYA8 (H8) and OVCAR8 (O8) cell lines. Error bars represent standard deviations. *, **, and *** indicate P < 0.1, 0.01 and 0.001. (C) Predicted binding sites in the YAP1 3′ UTR. (D) YAP1 luciferase reporter activity for the two predicted binding sites, 48 h after transfecting a miR-509-3p mimic into HEYA8 and OVCAR8 cells. NT = no treatment. miR-NC = scrambled miR-509-3p mimic. Error bars represent standard deviations. *** indicates P < 0.001. (E) Yap1 immunofluorescence staining (middle panel) of non-transfected (NT) HEYA8 and OVCAR8 cells or cells transfected with miR-509-3p, Yap1 siRNA (siYap1) or a negative control scrambled miRNA (NC) for 72 hr. Cell nuclei were stained with DAPI (left panel). Scale bars represent 50 um. (F) Kaplan-Meier curves for overall patient survival for high vs. low YAP1 abundance, using TCGA sequencing data and a 17 RPKM threshold [20].

The effect of miR-509-3p and siYAP1 on cell migration and 3D aggregate formation for HEYA8 (left) and OVCAR8 (right) cells.

Figure 5: The effect of miR-509-3p and siYAP1 on cell migration and 3D aggregate formation for HEYA8 (left) and OVCAR8 (right) cells. (A) In vitro 96-well migration time-course results for miR-509-3p transfection and for siRNA for YAP1, in HEYA8 and OVCAR8 cells. (B) Time course evolution of spheroidal aggregation. (C, D) Phase contrast and (E, F) fluorescence micrographs of time-course experiments. For treated cases in (A–F), time courses start (0 hr) 72 h after transfection. NT: no treatment. NC: treatment with a scrambled miR-509-3p mimic.

The effect of miR-509-3p and siYAP1 on migration and 3D aggregate formation for SKOV3, OVCAR3, 4 and 5 cell lines.

Figure 6: The effect of miR-509-3p and siYAP1 on migration and 3D aggregate formation for SKOV3, OVCAR3, 4 and 5 cell lines. (A) In vitro 96-well migration time-course results for miR-509-3p transfection and siRNA for YAP1. The time courses start 72 h after transfection. (B) Fluorescence micrographs of time course evolution of spheroidal aggregation. NT: no treatment. NC: treatment with a scrambled miR-509-3p mimic.

miR-509-3p influences transcript levels of genes involved in migration and invasion, and directly targets YAP1

We identified which of 143 miR-509-3p target genes predicted by TargetScan v6.2 (www.targetscan.org) have been reported as involved in cellular migration and invasion, then used quantitative RT-PCR with gene-specific primers to evaluate transcript levels for the predicted targets BCAR1 [28], GPC6 [29], KCNMA1 [30], PBX3 [31], and YAP1 [32, 33]. Cells transfected with miR-509-3p mimics showed strong down-regulation of four of the five genes in OVCAR8 cells and three of the five genes in HEYA8 cells, relative to SCR treated and untreated controls (Figure S3E and S3F). Because YAP1 is an oncogene in ovarian cancer, and a high YAP1 level is associated with poor ovarian survival in TCGA data (Figure 4F), our independent cohort of 157 patients and in an independent study [34], we characterized YAP1 in more detail.

YAP1 is a transcriptional co-activator and key downstream nuclear effector of the Hippo signaling pathway [32, 33]. YAP1 protein level was downregulated in response to miR-509-3p mimics in both HEYA8 and OVCAR8 cell lines at 72, 96 and 120 h after transfection (Figure S3G). To demonstrate that the gene is a critical downstream mediator of miR-509-3p, we first showed that the downregulation seen from qPCR, Western blots and immunofluorescence staining (Figure 4A, 4B and 4E) is due to direct repression through one or both of the two predicted miR-509-3p binding sites in the YAP1 3′ UTR (Figure 4C). In luciferase reporter assays, both predicted target sites responded independently to miR-509-3p mimics (Figure 4D), and a mutation at either site was sufficient to abrogate miR-509-3p-mediated repression of the reporter. Then, we confirmed that low YAP1 abundance was associated (P = 0.018) with better survival in the TCGA cohort (Figure 4F) [20]. Finally, we found that siRNA to YAP1 impaired migration approximately 90% in OVCAR8 as effectively as miR-509-3p (Figure 5A). However, siYAP1 had no observable effect on migration in HEYA8 (Figure 5A), SKOV3, OVCAR3, OVCAR4 and OVCAR5 (Figure 6A). Our qPCR and Western blots showed that miR-509-3p repression on YAP1 was largely at the protein rather than the transcript level (Figure 4A, 4B and Figure S6A, S6B). We also found that the level of YAP1 protein was over 2.4 fold higher in OVCAR8 than in HEYA8, SKOV3, OVCAR3, OVCAR4 and OVCAR5 (Figure S6C). Consistent with this, miR-509-3p had the largest influence on cellular migration in OVCAR8, and only marginal influence in the other cell lines, all of which expressed moderate to low levels of YAP1. We performed a rescue experiment to show that overexpression of exogenous YAP1 largely reversed the inhibitory effects of miR-509-3p on migration in OVCAR8 cells. OVCAR8 cells co-transfected with a miR-509-3p mimic and a YAP1 overexpression plasmid (pCDNA-FLAG-YAP1) had migrated 158% faster at 40 h than cells co-transfected with a miR-509-3p mimic and an empty vector C-FLAG-pCDNA3 (P = 0.040) (Figure S7). These data suggests that YAP1 is likely an important prognostic marker for EOC and that miR-509-3p/YAP1 axis may provide new druggable targets for ovarian tumors with high levels of YAP1.

miR-509-3p mimics and YAP1 siRNA impair spheroid formation in 3D culture models of ovarian cancer

In ovarian cancer, compact multicellular spheroids are important in metastasis, which occurs primarily through the peritoneal fluid or ascites [35]. We used a magnetic levitation-based 3D culturing system [36, 37] and image analysis of time-course phase contrast and fluorescence micrographs to evaluate the impact of miR-509-3p and YAP1 on spheroid formation in HEYA8 and OVCAR8 cells (Figure 5B5F). We compared untreated (NT) to cells treated with a negative control (NC), miR-509-3p or siRNA to YAP1, starting time courses 72 h after transfection. As we had seen for migration and invasion, the treatments had different effects on the two cell lines. For HEYA8 cells, miR-509-3p and siYAP1 treatments attenuated aggregation moderately more than NT and NC treatments (Figure 5B5D). For OVCAR8 cells miR-509-3p treatment strongly attenuated cellular aggregation into spheroids, with treated cells remaining largely dispersed, while siYAP1 treatment attenuated aggregation at a level intermediate between NC and miR-509-3p treatment (Figure 5B, 5E and 5F). We confirmed that spheroid formation was also attenuated in SKOV3, OVCAR3, OVCAR4 and OVCAR5 cell lines, with the magnitude of the attenuation varying between cell lines (Figure 6B). We also tested miR-508-3p, miR-508-5p, miR-509-3p and miR-509-5p effects on spheroid formation (Figure S8); of these, only miR-509-3p attenuated cellular aggregation into spheroids in OVCAR8 cells. This is consistent with our results showing that, of the Xq27.3 cluster members, miR-509-3p was the most clinically significant in relation to overall survival of ovarian cancer patients.

DISCUSSION

In the work reported here, we identified microRNA-509-3p as a clinically significant microRNA by integrating functional assays and clinical outcomes with microRNA and RNA sequencing data from over 250 tumors from The Cancer Genome Atlas (TCGA) ovarian cancer cohort. We used in situ hybridization (ISH) to confirm the association of miR-509-3p with survival in an independent cohort of 157 ovarian cancer tumors. miR-509-3p attenuated cellular migration and spheroid formation in all six ovarian cancer cell lines tested. In ovarian cancer cell lines, our results established that miR-509-3p can functionally alter YAP1 levels, and can attenuate migration, invasion, and aggregation into 3D spheroids.

miR-509-3p is expressed from a ~100-kb genomic cluster of miRNAs on Xq27.3. Members of this miR cluster have been associated with stage and survival in ovarian cancer, and are less abundant in omental metastases than in primary tumors [3846]. In renal cancer, cluster members are less abundant than in adjacent normal tissue, and overexpression can reduce proliferation and migration, while inducing apoptosis [47, 48]. In contrast, in melanoma, low abundance of cluster members is associated with impaired growth, invasiveness and colony formation, and with increased apoptosis [49]. Recently, miR-509-3p was confirmed to directly target XIAP and inhibit proliferation and increase sensitivity to cisplatin in chemoresistant ovarian cancer cells [50]. Members of the Xq27.3 miR cluster have been reported to be strongly anti-correlated to a multi-cancer ‘metastasis-associated fibroblast’ gene signature [51]. miR-509-3p has recently been shown to target CDK2 and to influence the cell cycle, colony formation and migration in human epithelial lung and cervical cancer cell lines [52]. Another family member, miR-506, and miR-200-family members, have been associated with EMT in a mesenchymal ovarian cancer subtype using the published microarray-based data [7, 43]. miR-509-3p and other Xq27.3 cluster miRs have recently been reported [53] as more abundant in a high-grade serous subtype (C5) that is defined by genes expressed in mesenchymal development [24], and less abundant in a subtype (C1) that has the poorest overall survival and is enriched in stromal genes.

Our results show that YAP1 is a direct downstream target of miR-509-3p, and is both a major effector of miR-509-3p-mediated attenuation of migration and invasion, and an important effector in spheroid formation in ovarian cancer cells that contain higher levels of YAP1. YAP1, on 11q22, is frequently subject to copy number gains, mutations, and inappropriate expression in diverse types of cancer, and is considered an oncogene in ovarian and other cancers [34, 5460]. Activation of YAP1 in ovarian tumors is negatively correlated with clinical outcomes and response to taxanes [61].

Curiously, while YAP1 is associated with the Hippo pathway, the mRNAs that were strongly inversely correlated to miR-509-3p in TCGA ovarian tumors were enriched not for Hippo pathway genes but for structural and regulatory components of the ECM. YAP1 is functionally associated with the ECM through its non-canonical, Hippo-independent role as a mechanotransducer, and, as a central effector in mechanotransduction, YAP1 influences how tumor cells sense and respond to the mechanical properties of the ECM and to their microenvironment [32, 33]. ECM changes have been reported to influence cell proliferation, differentiation, migration and invasion during normal development and differentiation [62], and mutations in mechanotransduction mediators, which are sensed by the cell as perturbations in ECM stiffness, have been associated with a number of diseases, including cancer [32, 33]. ECM rigidity and YAP1 activation can interact to maintain cancer-associated fibroblasts [63].

A tightly regulated balance between ECM deposition and turnover is critical for tissue homeostasis, and can be impaired in cancer cells, influencing migration, invasion and metastatic progression [35]. This is especially relevant to epithelial ovarian cancer, which primarily metastasizes by cells and multicellular spheroids exfoliating into the peritoneal fluid [35]. Multi-cellular spheroids can be drug resistant, due in part to the decreased penetrance of therapeutics to the protected interior cells [64]. Spheroid formation allows cancer cells to overcome anoikis, a form of apoptosis that can be triggered in individual cells that are unable to attach to other cells or to the ECM [65].

The model in Figure 7 summarizes our findings for miR-509-3p and the YAP1/ECM axis. Collectively, our results suggest that miR-509-3p-mediated changes in levels of YAP1 and ECM genes impair migration, invasion, and spheroid formation and so may attenuate metastatic progression in advanced stage ovarian cancer. Further, our results suggest that the direct downstream miR-509-3p target YAP1 is likely a critical driver of cellular migration and spheroid formation in ovarian cancers in which levels of YAP1 protein are high. We suggest that miR-509-3p and drugs that target the miR-509-3p/YAP1/ECM axis may offer novel therapeutic opportunities in ovarian cancer, and potentially in other cancers.

Model for the proposed mechanisms through which miR-509-3p alters the YAP1/ECM axis.

Figure 7: Model for the proposed mechanisms through which miR-509-3p alters the YAP1/ECM axis. Arrows and Tees indicate genes that were respectively up- and down-regulated in response to miR-509-3p (Figure S2). ECM functional gene groups are from [21], and names along the ECM/YAP1 axis were adapted from Figure 1 in [75].

MATERIALS AND METHODS

Library construction, sequencing, data processing

mRNA and miRNA sequencing data were generated using methods previously described [66]. For miRNA abundance we used a 5p and 3p mature strand data representation.

Analysis data sets

Analyses involving mRNA-seq and miRNA-seq used 289 datasets present in both data types. For miRNA-seq, unsupervised clustering, subset analysis and differential abundance calculations used data for 475 tumor samples. The hg19 (GRCh37) reference human genome was used for both data types. See Table S1.

miR-to-gene correlations

As described recently [67], we calculated correlations for RPM data for each miRNA-seq 5p or 3p strand and each mRNA-seq transcript from reference-based assembly [10]. The P-values obtained from the Spearman correlations were corrected for multiple testing using the Benjamini Hochberg method. Corrections for multiple testing were done per miRNA, and a q-value < 0.05 significance threshold was applied. TargetScan v6.0 binding site predictions [11] were calculated on the full-length sequences of the predicted transcripts. Enriched GO process, function and cellular locations, and enriched KEGG pathways, were calculated using DAVID [68] and mgsa [69].

Independent component analysis, miICA

In RPKM (mRNAs) and RPM (miRNA 5p or 3p strand) expression matrices, genes and miRs with a maximum RPKM or RPM < 10 respectively were removed, and only samples with survival data and both mRNA and miRNA expression profiles were retained in the expression matrices (denoted Em and Emi respectively). For genes for which the expression matrix generated from reference-based assembly contains results for more than one splice variant, we sum the expression for the individual spice-variants across the samples, generating a matrix of total mRNA expression that contains 12664 mRNAs across the 283 samples. We use the R-function rm.outlier() to remove outlier mRNAs and miRNAs with a coefficient of variation g (i.e. standard deviation divided by mean expression) greater than 5. As is done when ICA is applied to expression matrices from microarray experiments, we log2-transformed the miRNA and mRNA expression matrices, substituting small values for zeros.

Independent component analysis (ICA) [13] decomposes an expression matrix (E, GxS) into a component matrix (C) and a mixing matrix (M), i.e.

where C is the GxI component matrix containing the independent components (ICs) in the columns, and M is the IxS mixing matrix giving the linear mixing of each IC in each sample. The convention is to scale the ICs so that in each IC the entries Cgi have a mean of zero and a variance of one.

Component weights Mis measure the activity level of transcriptional programs, e.g. biological processes, in an IC across the samples. For example, when an IC maps to a cancer-related pathway, then this IC’s weights will be different in cancerous samples vs. non-cancerous samples. Component weights can be used to classify samples according to physiological, histological, or other features [70, 71]. Regarding gene weights Cgi, if one gene is inhibited by a transcriptional program while another is activated, then their gene weights will be anti-correlated, e.g. one being negative and the other positive.

Before applying ICA, the number of independent components (ICs), Imin (N,S), needs to be specified. To determine the optimal number of ICs to be estimated we used random matrix theory, as implemented in the EstDimRMT() R-function from the isva (v1.3) R package. The maximum of EstDimRMT() applied to the miRNA and the mRNA expression matrices, separately, was used as the number of ICs (I) to be estimated.

We used the default fastICA algorithm (v1.1–16) to first pre-whiten the expression matrix (E) by projecting it onto the first I principal components in the matrix (K). Subsequently we estimated an un-mixing matrix W such that C = EKW contains columns that are maximally independent, i.e. the independent components.

We extended the fastICA package to address genes and miRNAs simultaneously, i.e. to enable estimating an un-mixing matrix that simultaneously fulfills

and maximizes the independence of the ICs contained in both Cm and Cmi. Subscripts m and mi denote mRNA and miRNA, respectively. The mixing matrices are calculated from the estimated un-mixing matrix by solving

for Mm and Mmi.

For the ICs to be independent they must have a maximum negentropy [13, 72], i.e. the distribution of the entries Cgi in each IC is non-Gaussian. fastICA, and our modified fastICA, iteratively maximize the negentropy of the ICs using random seeding. Depending on the initial random guess for the un-mixing matrix, the ICs that do not represent structures in the data may converge to local, rather than to global maxima of negentropy. However, ICA is expected to always extract an IC that identifies structures in the data, i.e. to always converge to approximately the same negentropy maximum. To quantify this we run the modified fastICA 30 times. For each run, we calculate Spearman’s correlation coefficients between each of the original ICs and the ICs in the 29 other runs (algorithm adapted from [73]), and store the highest absolute correlation coefficient calculated for each IC in run r. We average the stored absolute correlation coefficients component-wise, take the average correlation coefficients as a measure for the stability of an IC, and order ICs according to their stability. We applied our modified fastICA() to the expression matrices containing the 283 samples that had survival and progression data, as well as miRNA and mRNA expression profiles

Following ICA, one way to identify significant genes in an IC is by finding genes that satisfy the inequality

where SD|C.i| is the standard deviation of the absolute gene weights in IC i, and α is often set to a value between 2 and 3 [16]. Here we use α = 3. Because we are interested in the biological processes that are related to the survival and recurrence times, for each IC we determine an optimal cutoff for the component weights by using Cutoff Finder [20] v2.0 to divide the patients into good and poor survival or progression groups.

We performed a pathway analysis on each IC by testing for overrepresentation of pathway gene-sets from the MsigDB (v3.0) (broadinstitute.org/gsea/msigdb/). Denoting the total number of genes annotated for a specific pathway (P) by Gp, and the number of significant genes in IC i that are annotated by P as Kp,I, we apply a hypergeometric test to test whether the probability P(X > Kp,i) is significant when there are N genes in total and Ni significant genes in IC i.

All p-values are FDR-corrected using the R function p.adjust.

For miR target enrichment calculations, a miR binding site was taken to be a 7mer seed match in the longest annotated Ensembl vGRCh37.p11 3′UTR. Statistical significance was calculated using q-values for a hypergeometric test.

miR subsets that discriminate clinical outcomes

To identify a concise set of miRNA mature strands whose expression was strongly associated with overall survival, we extracted combinatorial subsets of all mature strands from members of the Xq27.3 miR cluster, and from the miR-29 and miR-200-families, separately, and tested each subset for its log-rank P-value using the R survival package. For each of the three sets of 5p and 3p strands, we generated all possible subsets, and for each subset we separated samples into two groups using k-means clustering on abundance (RPM) data that we normalized between samples, and then again between miRNAs, such that the range of abundance was not a factor in the clustering. We tested Kaplan-Meier (KM) curves for each pair of sample groups for statistical significance, correcting log-rank P-values for multiple testing using the Benjamini Hochberg method. Representative subsets for which we displayed KM plots had a both high-ranked significant KM log-rank P-value and relatively large numbers of samples in the smaller of the two sample groups.

Cell cultures and transient transfections

Ovarian cancer cell lines OVCAR8 and HEYA8 were maintained in RPMI media (Life Technologies) supplemented with 10% heat inactivated fetal bovine serum (Atlanta Biologicals) and 1% of penicillin/streptomycin. Cell lines were incubated at 37°C with 90% humidity and 5% CO2. miRNA mimics/YAP1 siRNA (Life Technologies) and negative controls (Life Technologies) were transiently transfected at a 20 nM concentration using Lipofectamine 3000 (Life Technologies) following the manufacturer’s recommendations. For rescue experiments, pCDNA Flag-Yap1/empty vector C-FLAG-pCDNA3 were transfected at 3.6 ug/well in a 6-well plate using FuGene HD (Promega) 24 hrs after miRNA transfection, following the manufacturer’s recommendations.

Scratch-wound migration assay

Scratch-wound assays were performed using the Essen BioScience Cell Player Migration System. Seventy-two h after transfection, cells were seeded in 96-well ImageLock microplates (Essen BioScience) at 35,000 cells per well. Artificial scratch wounds were generated using the 96-well WoundMaker™. Each well was washed twice with PBS, and then overlaid with 200 μL of RPMI 1640 medium supplemented with 10% FBS. Images were acquired using the IncuCyte-FLR™ live cell imaging platform. Cell migration was calculated with IncuCyte’s automated image analysis algorithm using the Relative Wound Density (RWD) metric, which is defined by

where w(t) is the density of the wound region and c(t) is the density of the cell region at time t. Density was plotted as a function of time for three independent experiments, with error bars showing standard deviations. Statistical significance was determined using one-way ANOVA with the Tukey posthoc test (IBM SPSS Statistics 21). Results reported in Figures 3 and 5 were from independent experiments in two laboratories.

For the YAP1 overexpression rescue experiment, scratches were made using 200 ul pipette tips. Cells were washed two times by PBS and fresh media was added. The cell culture plates were transferred to a Biostation CT (Nikon) instrument programmed to take images every 4 h for up to 40 h. The images at different time points were aligned using the NIS-Element AR software. Scratch edges were detected and scratch area was masked with red color for all the images using NIS-Element AR software. Binary area fraction (fraction of masked area with red color) for each time point was measured using NIS-Element AR software. Area covered due to migration of cells after each time point was obtained by subtracting binary area fraction for each time point from that of initial time point. Average migration area for each time point was calculated for three biological replicates.

Matrigel invasion assay

Cell invasion assays were performed using Matrigel-coated 96-well microplates. Scratch-wounds were generated as described for the scratch-wound migration assay, and 50 μL of Matrigel Matrix (BD Biosciences), diluted to 1 mg/mL in cell culture medium, was added to each well. Microplates were incubated at 37°C to set the Matrigel, and was overlaid with 200 μL of culture medium 24 h later. Density was plotted as a function of time for three independent experiments, with error bars showing standard error deviations.

3D culture assay

To form three dimensional (3D) cell aggregates, cells, with or without miRNA/siRNA transfection, were incubated overnight in RPMI with 10% FBS, with a magnetic gold–polymer–iron oxide hydrogel (Nanoshuttle-PL, NS) (n3D Biosciences, Houston TX) to allow the attachment of NS to the cells [36, 37, 74]. 5 uM of CellTracker Green CMFDA Dye (Molecular Probes) was added to the cells and incubated for 30 minutes. The cells were washed with PBS, trypsinized, resuspended in RPMI with 10% FBS, and then magnetically levitated to the air-water interface for 4 hours, where they formed 3D structures with an ECM. Next, the 3D structures were mechanically disrupted by pipetting, and the resulting suspensions of aggregates were transferred to a 96-well plate, where they were formed or ‘printed’ into either a dot or a ring shape by placing a dot or ring magnet ‘drive’ underneath the plate for 15 minutes. Using either phase contrast or fluorescence, the 3D culture area was automatically imaged over time by an IncuCyte ZOOM live-cell imaging system (Essen BioSciences, Ann Arbor, Michigan), and the area of the dense 3D spheroidal regions was recorded for each well.

Quantitative polymerase chain reaction (qPCR)

Total RNA from cellular pellets was isolated using an miRNeasy Mini Kit (Qiagen) following the manufacturer’s instructions. Purity and RNA concentration were measured using a ND-100 Nanodrop spectrophotometer (Thermo Scientific). Reverse transcription of mRNAs was performed using the TaqMan Reverse Transcription Reagents (Invitrogen) using 300 ng of total RNA for each 20 μL reaction. qPCR was performed using the SYBR Green reagent (Applied Biosystems) in a 10 μL reaction with 1 μL of cDNA and 0.005 nmol of primers. Reactions were run on a StepOnePlus Real-Time PCR system (Life Technologies) and analyzed using the delta-delta CT method. The level of 18S ribosomal RNA was used to normalize the relative expression of genes.

In situ hybridization

A tissue microarray with high-grade serous ovarian cancer samples, as determined by a gynecologic pathologist, was used. Formalin-fixed paraffin-embedded tissue sections were dewaxed in xylene, and rehydrated through an ethanol dilution series. Tissue sections were digested with 15 μg/mL proteinase K for 20 min at room temperature, were then loaded onto a Discovery Ultra slide staining system (Ventana/Roche) for in situ hybridization. The tissue microarray slides were incubated with a double DIG-labeled miRCURY LNA microRNA probe (Exiqon) for 2 h. The digoxigenin was then detected with a polyclonal anti-DIG antibody and alkaline phosphatase-conjugated second antibody (Ventana) using NBT-BCIP as the substrate. We confirmed that ovarian tissue has very low background staining, using a double DIG-labeled LNA negative control probe (Exiqon) and the same procedure as above.

The expression levels for miR-509-3p were determined as previously described [26]. In brief, representative images from each tissue section were imaged at 10 × magnification. Blinded from all clinical data, CellProfiler 2.0 software [75] was trained with negative and positive controls to establish staining intensity threshold levels and cell diameter parameters for enumerating the number of positively stained cells at 100 × magnification. These parameters were used to score all images.

Luciferase reporter assay

The DNA sequence, flanked by ~ 600 nucleotides, with the two predicted miR-509-3p target sites at -329 to -336 nt in the 3483-nt YAP1 3′UTR (YAP1_1: -329 bp to -336 bp, YAP1_2: -3084 bp to -3091bp) were PCR-amplified from human genomic DNA and sub-cloned into the NotI restriction site of a psiCHECK-2 vector (Promega) using the Gibson Assembly Cloning Kit (New England Biolabs). We generated mutant clones by mutating the miR-509-3p seed interaction sites in the reporter plasmid construct, using the Q5 Site-Directed Mutagenesis kit (New England Biolabs) for YAP1_1 (ACCAATCA->ACgAtTCA) and the QuickChange II XL Site-Directed Mutagenesis kit (Agilent Technologies) for YAP1_2 (ACCAATCA->gCCAAgCA). Both clones were verified by Sanger sequencing. 100ng of each reporter plasmid was transfected alone, or co-transfected with either 20 nM of a miR-509-3p mimic, or a scrambled negative control miRNA mimic, into HEYA8 and OVCAR8 cells, using lipofectamine 2000 (Life Technologies) at 10000 cells/well in a 96-well plate, following the manufacturer’s recommended protocol. Luciferase assays were performed at 48 h post transfection using a Dual-Luciferase Reporter Assay System (Promega). Renilla luminescence signals were normalized by the firefly luciferase signals. Three replicates were performed for each assay.

Relationship of YAP1 abundance and survival

An RPKM abundance threshold giving the most significant high-vs-low log-rank P-value for overall survival was identified with CutoffFinder v2.0 [20].

ECM functional gene groups

Gene symbols for each ECM functional category [21] were taken from matrisome-table-Hs_CURRENT.xlsx, which was downloaded from web.mit.edu/hyneslab/matrisome.

Statistical significance of overlap between two gene lists

P-values were calculated using a custom Mathematica (Wolfram Research, Champaign, IL) script that implemented the method described at nemates.org/MA/progs/representation.stats.html. For example, given a population of npop = 21929 genes, d = 1368 predicted target genes for the Xq27.3 cluster members and n = 200 ECM glycoproteins from matrisome-table-Hs_CURRENT.xlsx from web.mit.edu/hyneslab/matrisome [21] had an overlap of x = 38 genes, for which P = 6.3E-10.

Immunofluorescence staining

Cells were seeded onto a Falcon culture slide (Corning) and cultured in complete RPMI medium under standard cell culture conditions. Cells were fixed in 4% paraformaldehyde at room temperature for 15 min, then were permeabilized in 1 × PBS containing 0.1% Tween-100 for 30 min at room temperature. The cells were blocked in a solution of 1 × PBS containing 10% normal goat serum and 0.1% Triton X-100 for 4 h. After being washed twice with 1 × PBS for 5 min each time, the cells were incubated overnight with a 1:200 dilution of monoclonal rabbit Yap1 antibody (Cell Signaling Technology) in the blocking solution at 4°C. After being washed twice with 1 × PBS for 5 min each time, a goat anti-rabbit IgG conjugated with Alexa Fluor 568 (Life Technologies) at 1:1000 in the blocking solution was incubated with the cells at room temperature for 1 hr. The cells were stained with DAPI for 5 min after being washed twice with 1 × PBS, and then mounted with ProLong Gold Antifade Mountant (Life Technologies). The fluorescence images were acquired using an Olympus imaging microscope with a 20 × objective.

Western blot

Total proteins were extracted from cellular pellets using protein extraction buffer (Triton X-100 1%, NaCl 150 mM, Tris 25 mM, pH 7.6) with Complete mini protease inhibitors (Roche). Thirty micrograms of protein were run on Nupage-NOVEX Bis-Tris 4–12% gels using a MOPS buffer. After electrophoresis proteins were transferred to a PVDF membrane (Invitrogen), and the membranes were then blocked for 1 h at room temperature with 5% non-fat milk dissolved on TBS-T (50 mM Tris-HCl, 150 mM NaCl, 0.1% Tween 20, pH 7.5). Primary antibodies Yap1 and GAPDH (Cell Signaling Technology) were added at a concentration of 1:1000 and were incubated overnight with the membrane at 4°C. HRP-linked anti-rabbit secondary antibodies (Cell Signaling Technology) were added at a concentration of 1:10,000 and incubated with membrane at room temperature for 1 hr. The proteins were visualized using a Pierce ECL Plus chemiluminescent substrate (Thermo Scientific).

ACKNOWLEDGMENTS

The authors gratefully acknowledge the data production teams at the Genome Sciences Centre.

GRANT SUPPORT

Research reported in this publication was supported by the National Cancer Institutes of Health under Award Number U24CA143866, awarded to MAM. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. PHG is supported by a grant from the Cancer Prevention and Research Institute of Texas (RP110355) and generous support from the McNair Foundation. AH-H is supported in part by a fellowship from The National Council on Science and Technology (CONACYT) in Mexico. ACR and CKC is supported by National Science Foundation (CAREER DBI-1254185 to ACR) and the University of California Cancer Research Coordinating Committee. This work is also supported in part by grants from the NIH (U54CA151668, UH2TR000943), and CPRIT (RP110595).

CONFLICTS OF INTEREST

The authors declare no conflicts of interests.

REFERENCES

1. Vaughan S, Coward JI, Bast RC Jr, Berchuck A, Berek JS, Brenton JD, Coukos G, Crum CC, Drapkin R, Etemadmoghadam D, Friedlander M, Gabra H, Kaye SB, et al. Rethinking ovarian cancer: recommendations for improving outcomes. Nat Rev Cancer. 2011; 11:719–25.

2. Shield K, Riley C, Quinn MA, Rice GE, Ackland ML, Ahmed N. Alpha2beta1 integrin affects metastatic potential of ovarian carcinoma spheroids by supporting disaggregation and proteolysis. J Carcinog. 2007; 6:11.

3. Zietarska M, Maugard CM, Filali-Mouhim A, Alam-Fahmy M, Tonin PN, Provencher DM, Mes-Masson AM. Molecular description of a 3D in vitro model for the study of epithelial ovarian cancer (EOC). Mol Carcinog. 2007; 46:872–85.

4. Burleson KM, Boente MP, Pambuccian SE, Skubitz AP. Disaggregation and invasion of ovarian carcinoma ascites spheroids. J Transl Med. 2006; 4:6.

5. L’Esperance S, Bachvarova M, Tetu B, Mes-Masson AM, Bachvarov D. Global gene expression analysis of early response to chemotherapy treatment in ovarian cancer spheroids. BMC Genomics. 2008; 9:99.

6. Liao J, Qian F, Tchabo N, Mhawech-Fauceglia P, Beck A, Qian Z, Wang X, Huss WJ, Lele SB, Morrison CD, Odunsi K. Ovarian cancer spheroid cells with stem cell-like properties contribute to tumor generation, metastasis and chemotherapy resistance through hypoxia-resistant metabolism. PLoS One. 2014; 9:e84941.

7. Cancer Genome Atlas Research N. Integrated genomic analyses of ovarian carcinoma. Nature. 2011; 474:609–15.

8. Farazi TA, Hoell JI, Morozov P, Tuschl T. MicroRNAs in human cancer. Adv Exp Med Biol. 2013; 774:1–20.

9. Git A, Dvinge H, Salmon-Divon M, Osborne M, Kutter C, Hadfield J, Bertone P, Caldas C. Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. 2010; 16:991–1006.

10. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28:511–5.

11. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP., MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007; 27:91–105.

12. Creighton CJ, Hernandez-Herrera A, Jacobsen A, Levine DA, Mankoo P, Schultz N, Du Y, Zhang Y, Larsson E, Sheridan R, Xiao W, Spellman PT, Getz G, et al. Integrated analyses of microRNAs demonstrate their widespread influence on gene expression in high-grade serous ovarian carcinoma. PLoS One. 2012; 7:e34546.

13. Comon P. Independent Component Analysis, a New Concept. Signal Processing. 1994; 36:287–314.

14. Kim RS, Ji H, Wong WH. An improved distance measure between the expression profiles linking co-expression and co-regulation in mouse. BMC Bioinformatics. 2006; 7:44.

15. Stanhope SA, Sengupta S, den Boon J, Ahlquist P, Newton MA. Statistical use of argonaute expression and RISC assembly in microRNA target identification. PLoS Comput Biol. 2009; 5:e1000516.

16. Teschendorff AE, Journee M, Absil PA, Sepulchre R, Caldas C. Elucidating the altered transcriptional programs in breast cancer using independent component analysis. PLoS Comput Biol. 2007; 3:e161.

17. Frigyesi A, Veerla S, Lindgren D, Hoglund M, Independent component analysis reveals new and biologically significant structures in micro array data. BMC Bioinformatics. 2006; 7:290.

18. Kong W, Mou X, Hu X. Exploring matrix factorization techniques for significant genes identification of Alzheimer’s disease microarray gene expression data. BMC Bioinformatics. 2011; 12:S7.

19. Zinovyev A, Kairov U, Karpenyuk T, Ramanculov E. Blind source separation methods for deconvolution of complex signals in cancer biology. Biochem Biophys Res Commun. 2013; 430:1182–7.

20. Budczies J, Klauschen F, Sinn BV, Gyorffy B, Schmitt WD, Darb-Esfahani S, Denkert C. Cutoff Finder: a comprehensive and straightforward Web application enabling rapid biomarker cutoff optimization. PLoS One. 2012; 7:e51862.

21. Hynes RO, Naba A Overview of the matrisome—an inventory of extracellular matrix constituents and functions. Cold Spring Harb Perspect Biol. 2012; 4:a004903.

22. Zhang W, Ota T, Shridhar V, Chien J, Wu B, Kuang R. Network-based survival analysis reveals subnetwork signatures for predicting outcomes of ovarian cancer treatment. PLoS Comput Biol. 2013; 9:e1002975.

23. Bonome T, Levine DA, Shih J, Randonovich M, Pise-Masison CA, Bogomolniy F, Ozbun L, Brady J, Barrett JC, Boyd J, Birrer MJ. A gene signature predicting for survival in suboptimally debulked patients with ovarian cancer. Cancer Res. 2008; 68:5478–86.

24. Tothill RW, Tinker AV, George J, Brown R, Fox SB, Lade S, Johnson DS, Trivett MK, Etemadmoghadam D, Locandro B, Traficante N, Fereday S, Hung JA, et al. Novel molecular subtypes of serous and endometrioid ovarian cancer linked to clinical outcome. Clin Cancer Res. 2008; 14:5198–208.

25. De Craene B, Berx G. Regulatory networks defining EMT during cancer initiation and progression. Nat Rev Cancer. 2013; 13:97–110.

26. Pecot CV, Rupaimoole R, Yang D, Akbani R, Ivan C, Lu C, Wu S, Han HD, Shah MY, Rodriguez-Aguayo C, Bottsford-Miller J, Liu Y, Kim SB, et al. Tumour angiogenesis regulation by the miR-200 family. Nat Commun. 2013; 4:2427.

27. Sokilde R, Kaczkowski B, Podolska A, Cirera S, Gorodkin J, Moller S, Litman T, Global microRNA analysis of the NCI-60 cancer cell panel. Mol Cancer Ther. 2011; 10:375–84.

28. Nick AM, Stone RL, Armaiz-Pena G, Ozpolat B, Tekedereli I, Graybill WS, Landen CN, Villares G, Vivas-Mejia P, Bottsford-Miller J, Kim HS, Lee JS, Kim SM, et al. Silencing of p130cas in ovarian carcinoma: a novel mechanism for tumor cell death. J Natl Cancer Inst. 2011; 103:1596–612.

29. Januchowski R, Zawierucha P, Rucinski M, Nowicki M, Zabel M. Extracellular matrix proteins expression profiling in chemoresistant variants of the A2780 ovarian cancer cell line. Biomed Res Int. 2014; 2014:365867.

30. Oeggerli M, Tian Y, Ruiz C, Wijker B, Sauter G, Obermann E, Guth U, Zlobec I, Sausbier M, Kunzelmann K, Bubendorf L. Role of KCNMA1 in breast cancer. PLoS One. 2012; 7:e41664.

31. Han HB, Gu J, Ji DB, Li ZW, Zhang Y, Zhao W, Wang LM, Zhang ZQ. PBX3 promotes migration and invasion of colorectal cancer cells via activation of MAPK/ERK signaling pathway. World J Gastroenterol. 2014; 20:18260–70.

32. Dupont S, Morsut L, Aragona M, Enzo E, Giulitti S, Cordenonsi M, Zanconato F, Le Digabel J, Forcato M, Bicciato S, Elvassore N, Piccolo S. Role of YAP/TAZ in mechanotransduction. Nature. 2011; 474:179–83.

33. Halder G, Dupont S, Piccolo S. Transduction of mechanical and cytoskeletal cues by YAP and TAZ. Nat Rev Mol Cell Biol. 2012; 13:591–600.

34. Xia Y, Chang T, Wang Y, Liu Y, Li W, Li M, Fan HY. YAP promotes ovarian cancer cell tumorigenesis and is indicative of a poor prognosis for ovarian cancer patients. PLoS One. 2014; 9:e91770.

35. Sodek KL, Murphy KJ, Brown TJ, Ringuette MJ. Cell-cell and cell-matrix dynamics in intraperitoneal cancer metastasis. Cancer Metastasis Rev. 2012; 31:397–414.

36. Souza GR, Molina JR, Raphael RM, Ozawa MG, Stark DJ, Levin CS, Bronk LF, Ananta JS, Mandelin J, Georgescu MM, Bankson JA, Gelovani JG, Killian TC. Three-dimensional tissue culture based on magnetic cell levitation. Nat Nanotechnol. 2010; 5:291–6.

37. Timm DM, Chen J, Sing D, Gage JA, Haisler WL, Neeley SK, Raphael RM, Dehghani M, Rosenblatt KP, Killian TC, Tseng H, Souza GR. A high-throughput three-dimensional cell migration assay for toxicity screening with mobile device-based macroscopic image analysis. Sci Rep. 2013; 3:3000.

38. Delfino KR, Rodriguez-Zas SL. Transcription factor-microRNA-target gene networks associated with ovarian cancer survival and recurrence. PLoS One. 2013; 8:e58608.

39. Eitan R, Kushnir M, Lithwick-Yanai G, David MB, Hoshen M, Glezerman M, Hod M, Sabah G, Rosenwald S, Levavi H. Tumor microRNA expression patterns associated with resistance to platinum based chemotherapy and survival in ovarian cancer patients. Gynecol Oncol. 2009; 114:253–9.

40. Lee CH, Subramanian S, Beck AH, Espinosa I, Senz J, Zhu SX, Huntsman D, van de Rijn M, Gilks CB. MicroRNA profiling of BRCA1/2 mutation-carrying and non-mutation-carrying high-grade serous carcinomas of ovary. PLoS One. 2009; 4:e7314.

41. Liu Y, Devescovi V, Chen S, Nardini C. Multilevel omic data integration in cancer cell lines: advanced annotation and emergent properties. BMC Syst Biol. 2013; 7:14.

42. Vang S, Wu HT, Fischer A, Miller DH, MacLaughlan S, Douglass E, Steinhoff M, Collins C, Smith PJ, Brard L, Brodsky AS. Identification of ovarian cancer metastatic miRNAs. PLoS One. 2013; 8:e58226.

43. Yang D, Sun Y, Hu L, Zheng H, Ji P, Pecot CV, Zhao Y, Reynolds S, Cheng H, Rupaimoole R, Cogdell D, Nykter M, Broaddus R, et al. Integrated analyses identify a master microRNA regulatory network for the mesenchymal subtype in serous ovarian cancer. Cancer Cell. 2013; 23:186–99.

44. Yu X, Zhang X, Bi T, Ding Y, Zhao J, Wang C, Jia T, Han D, Guo G, Wang B, Jiang J, Cui S. MiRNA expression signature for potentially predicting the prognosis of ovarian serous carcinoma. Tumour Biol. 2013; 34:3501–8.

45. Zhang L, Volinia S, Bonome T, Calin GA, Greshock J, Yang N, Liu CG, Giannakakis A, Alexiou P, Hasegawa K, Johnstone CN, Megraw MS, Adams S, et al. Genomic and epigenetic alterations deregulate microRNA expression in human epithelial ovarian cancer. Proc Natl Acad Sci U S A, 2008; 105:7004–9.

46. Zhang WB, Pan ZQ, Yang QS, Zheng XM. Tumor suppressive miR-509-5p contributes to cell migration, proliferation and antiapoptosis in renal cell carcinoma. Ir J Med Sci. 2013; 182:621–7.

47. Hidaka H, Seki N, Yoshino H, Yamasaki T, Yamada Y, Nohata N, Fuse M, Nakagawa M, Enokida H. Tumor suppressive microRNA-1285 regulates novel molecular targets: aberrant expression and functional significance in renal cell carcinoma. Oncotarget. 2012; 3:44–57. doi: 10.18632/oncotarget.417.

48. Zhai Q, Zhou L, Zhao C, Wan J, Yu Z, Guo X, Qin J, Chen J, Lu R. Identification of miR-508-3p and miR-509-3p that are associated with cell invasion and migration and involved in the apoptosis of renal cell carcinoma. Biochem Biophys Res Commun. 2012; 419:621–6.

49. Streicher KL, Zhu W, Lehmann KP, Georgantas RW, Morehouse CA, Brohawn P, Carrasco RA, Xiao Z, Tice DA, Higgs BW, Richman L, Jallal B, Ranade K, et al. A novel oncogenic role for the miRNA-506-514 cluster in initiating melanocyte transformation and promoting melanoma growth. Oncogene. 2012; 31:1558–70.

50. Chen W, Zeng W, Li X, Xiong W, Zhang M, Huang Y, Zhou L, Jiang S. MicroRNA-509-3p increases the sensitivity of epithelial ovarian cancer cells to cisplatin-induced apoptosis. Pharmacogenomics. 2016; 17:187–97.

51. Kim H, Watkinson J, Varadan V, Anastassiou D. Multi-cancer computational analysis reveals invasion-associated variant of desmoplastic reaction involving INHBA, THBS2 and COL11A1. BMC Med Genomics. 2010; 3:51.

52. Yoon S, Han E, Choi YC, Kee H, Jeong Y, Yoon J, Baek K, Inhibition of cell proliferation and migration by miR-509-3p that targets CDK2, Rac1, and PIK3C2A. Mol Cells. 2014; 37:314–21.

53. Patch AM, Christie EL, Etemadmoghadam D, Garsed DW, George J, Fereday S, Nones K, Cowin P, Alsop K, Bailey PJ, Kassahn KS, Newell F, Quinn MC, et al. Whole-genome characterization of chemoresistant ovarian cancer. Nature. 2015; 521:489–94.

54. Cai H, Xu Y. The role of LPA and YAP signaling in long-term migration of human ovarian cancer cells. Cell Commun Signal. 2013; 11:31.

55. Fernandez LA, Squatrito M, Northcott P, Awan A, Holland EC, Taylor MD, Nahle Z, Kenney AM. Oncogenic YAP promotes radioresistance and genomic instability in medulloblastoma through IGF2-mediated Akt activation. Oncogene. 2012; 31:1923–37.

56. Hall CA, Wang R, Miao J, Oliva E, Shen X, Wheeler T, Hilsenbeck SG, Orsulic S, Goode S. Hippo pathway effector Yap is an ovarian cancer oncogene. Cancer Res. 2010; 70:8517–25.

57. Lorenzetto E, Brenca M, Boeri M, Verri C, Piccinin E, Gasparini P, Facchinetti F, Rossi S, Salvatore G, Massimino M, Sozzi G, Maestro R, Modena P. YAP1 acts as oncogenic target of 11q22 amplification in multiple cancer subtypes. Oncotarget. 2014; 5:2608–21. doi: 10.18632/oncotarget.1844.

58. Steinhardt AA, Gayyed MF, Klein AP, Dong J, Maitra A, Pan D, Montgomery EA, Anders RA. Expression of Yes-associated protein in common solid tumors. Hum Pathol. 2008; 39:1582–9.

59. Wang Z, Wu Y, Wang H, Zhang Y, Mei L, Fang X, Zhang X, Zhang F, Chen H, Liu Y, Jiang Y, Sun S, Zheng Y, et al. Interplay of mevalonate and Hippo pathways regulates RHAMM transcription via YAP to modulate breast cancer cell motility. Proc Natl Acad Sci U S A, 2014; 111:E89–98.

60. Zhang X, George J, Deb S, Degoutin JL, Takano EA, Fox SB, group AS, Bowtell DD, Harvey KF. The Hippo pathway transcriptional co-activator, YAP, is an ovarian cancer oncogene. Oncogene. 2011; 30:2810–22.

61. Jeong W, Kim SB, Sohn BH, Park YY, Park ES, Kim SC, Kim SS, Johnson RL, Birrer M, Bowtell DS, Mills GB, Sood A, Lee JS. Activation of YAP1 is associated with poor prognosis and response to taxanes in ovarian cancer. Anticancer Res. 2014; 34:811–7.

62. Eyckmans J, Boudou T, Yu X, Chen CS. A hitchhiker’s guide to mechanobiology. Dev Cell. 2011; 21:35–47.

63. Calvo F, Ege N, Grande-Garcia A, Hooper S, Jenkins RP, Chaudhry SI, Harrington K, Williamson P, Moeendarbary E, Charras G, Sahai E. Mechanotransduction and YAP-dependent matrix remodelling is required for the generation and maintenance of cancer-associated fibroblasts. Nat Cell Biol. 2013; 15:637–46.

64. Desoize B, J Jardillier. Multicellular resistance: a paradigm for clinical resistance? Crit Rev Oncol Hematol. 2000; 36:193–207.

65. Santini MT, G Rainaldi, PL Indovina. Apoptosis, cell adhesion and the extracellular matrix in the three-dimensional growth of multicellular tumor spheroids. Crit Rev Oncol Hematol. 2000; 36:75–87.

66. Cancer Genome Atlas Research N, Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013; 499:43–9.

67. Lim EL, Trinh DL, Scott DW, Chu A, Krzywinski M, Zhao Y, Robertson AG, Mungall AJ, Schein J, Boyle M, Mottok A, Ennishi D, Johnson NA, et al. Comprehensive miRNA sequence analysis reveals survival differences in diffuse large B-cell lymphoma patients. Genome Biol. 2015; 16:18.

68. Huang da W, BT Sherman, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009; 4:44–57.

69. Bauer S, J Gagneur, PN Robinson. GOing Bayesian: model-based gene set analysis of genome-scale data. Nucleic Acids Res. 2010; 38:3523–32.

70. Huang DS, CH Zheng. Independent component analysis-based penalized discriminant method for tumor classification using gene expression data. Bioinformatics. 2006; 22:1855–62.

71. Zhang XW, Yap YL, Wei D, Chen F, Danchin A. Molecular diagnosis of human cancer type by gene expression profiles and independent component analysis. Eur J Hum Genet. 2005; 13:1303–11.

72. Hyvarinen A, Oja E. Independent component analysis: algorithms and applications. Neural Netw. 2000; 13:411–30.

73. Chiappetta P, Roubaud MC, Torresani B. Blind source separation and the analysis of microarray data. J Comput Biol. 2004; 11:1090–109.

74. Haisler WL, Timm DM, Gage JA, Tseng H, Killian TC, Souza GR. Three-dimensional cell culturing by magnetic levitation. Nat Protoc. 2013; 8:1940–9.

75. Lamprecht MR, Sabatini DM, Carpenter AE. CellProfiler: free, versatile software for automated biological image analysis. Biotechniques. 2007; 42:71–75.

76. Schiller HB, Fassler R. Mechanosensitivity and compositional dynamics of cell-matrix adhesions. EMBO reports. 2013; 14:509–519.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.
PII: 8412