MicroRNAs (miRNAs) are small, noncoding RNAs that can contribute to cancer development and progression by acting as oncogenes or tumor suppressor genes. Recent studies have also linked different sets of miRNAs to metastasis through either the promotion or suppression of this malignant process. Interestingly, epigenetic silencing of miRNAs with tumor suppressor features by CpG island hypermethylation is also emerging as a common hallmark of human tumors. Thus, we wondered whether there was a miRNA hypermethylation profile characteristic of human metastasis. We used a pharmacological and genomic approach to reveal this aberrant epigenetic silencing program by treating lymph node metastatic cancer cells with a DNA demethylating agent followed by hybridization to an expression microarray. Among the miRNAs that were reactivated upon drug treatment, miR-148a, miR-34b/c, and miR-9 were found to undergo specific hypermethylation-associated silencing in cancer cells compared with normal tissues. The reintroduction of miR-148a and miR-34b/c in cancer cells with epigenetic inactivation inhibited their motility, reduced tumor growth, and inhibited metastasis formation in xenograft models, with an associated down-regulation of the miRNA oncogenic target genes, such as C-MYC, E2F3, CDK6, and TGIF2. Most important, the involvement of miR-148a, miR-34b/c, and miR-9 hypermethylation in metastasis formation was also suggested in human primary malignancies ( n = 207) because it was significantly associated with the appearance of lymph node metastasis. Our findings indicate that DNA methylation-associated silencing of tumor suppressor miRNAs contributes to the development of human cancer metastasis.
Weighted gene coexpression network analysis (WGCNA) is a powerful ‘guilt-by-association’-based method to extract coexpressed groups of genes from large heterogeneous messenger RNA expression data sets. We have utilized WGCNA to identify 11 coregulated gene clusters across 2342 breast cancer samples from 13 microarray-based gene expression studies. A number of these transcriptional modules were found to be correlated to clinicopathological variables (e.g. tumor grade), survival endpoints for breast cancer as a whole (disease-free survival, distant disease-free survival and overall survival) and also its molecular subtypes (luminal A, luminal B, HER2+ and basal-like). Examples of findings arising from this work include the identification of a cluster of proliferation-related genes that when upregulated correlated to increased tumor grade and were associated with poor survival in general. The prognostic potential of novel genes, for example, ubiquitin-conjugating enzyme E2S (UBE2S) within this group was confirmed in an independent data set. In addition, gene clusters were also associated with survival for breast cancer molecular subtypes including a cluster of genes that was found to correlate with prognosis exclusively for basal-like breast cancer. The upregulation of several single genes within this coexpression cluster, for example, the potassium channel, subfamily K, member 5 (KCNK5) was associated with poor outcome for the basal-like molecular subtype. We have developed an online database to allow user-friendly access to the coexpression patterns and the survival analysis outputs uncovered in this study (available at http://glados.ucd.ie/Coexpression/ ).