Integrative epigenetic and genetic pan-cancer somatic alteration portraits

Genetic and epigenetic alterations are required for carcinogenesis and the mutation burden across tumor types has been investigated. Here, we investigate epigenetic alterations with a novel measure of global DNA methylation dysregulation, the methylation dysregulation index (MDI), across 14 cancer types in The Cancer Genome Atlas (TCGA) database. DNA methylation data—obtained using Illumina HumanMethylation450 BeadChip—was accessed from TCGA. We calculated the MDI in 14 tumor types (n = 5,592 tumors), using adjacent normal tissues (n = 701) from each tumor site. Copy number alteration, and mutation burden were retrieved from cBioportal (n = 5,152). We tested the relation of subject MDI across tumors and with age, gender, tumor stage, estimated tumor purity, and copy number alterations for both overall MDI and genomic-context-specific MDI. We also investigated the top most dysregulated loci shared across tumor types. There was a broad range of extent in methylation dysregulation across tumor types (P < 2.2E-16). However, a consistent pattern of methylation dysregulation stratified by genomic context was observed across tumor types where the highest dysregulation occurred at non-CpG island regions. Considering other summary measures of somatic alteration, MDI was correlated with copy number alterations but not with mutation burden. Using the top dysregulated CpG sites in common across tumors, 4 classes of cancer types were observed, and the functional consequences of these alterations to gene expression were confirmed. This work identified the global DNA methylation dysregulation patterns across 14 cancer types showing a higher impact for the non-CpG island areas. The most dysregulated loci across cancer types identified common clusters across cancer types that may have implications for future treatment and prevention measures.     

Age-related DNA methylation in normal breast tissue and its relationship with invasive breast tumor methylation

Age is a key risk factor for breast cancer and epigenetic alterations may contribute to age-related increases in breast cancer risk, though the relation of age-related methylation in normal breast tissues with altered methylation in breast tumors is unclear. We investigated the relation of age with DNA methylation in normal breast tissues genome-wide using two data sets from the Gene Expression Omnibus (GEO) database (GSE32393 and GSE31979). We validated our observations in an independent set of normal breast tissues, examined age-related methylation in normal breast for enrichment of genomic features, and compared age-related methylation in normal tissue with methylation alterations in breast tumors. Between the two array-based methylation data sets, there were 204 CpG loci with significant (P < 0.05) and consistent age-related methylation, 97% of which were increases in methylation. Our validation sets confirmed the direction of age-related DNA methylation changes in all measured regions. Among the 204 age-related CpG loci, we observed a significant enrichment for CpG islands (P = 8.7E-6) and polycomb group protein target genes (P = 0.03). In addition, 24 of the 204 CpGs with age-related methylation in normal breast were significantly differentially methylated between normal and breast tumor tissues. We identified consistent age-related methylation changes in normal breast tissue that are further altered in breast tumors and may represent early events contributing to breast carcinogenesis. This work identifies age-related methylation in normal breast tissue and begins to deconstruct the contribution of aging to epigenetic alterations present in breast tumors.     

Breast cancer family history and allele-specific DNA methylation in the legacy girls study

Family history, a well-established risk factor for breast cancer, can have both genetic and environmental contributions. Shared environment in families as well as epigenetic changes that also may be influenced by shared genetics and environment may also explain familial clustering of cancers. Epigenetic regulation, such as DNA methylation, can change the activity of a DNA segment without a change in the sequence; environmental exposures experienced across the life course can induce such changes. However, genetic-epigenetic interactions, detected as methylation quantitative trait loci (mQTLs; a.k.a. meQTLs) and haplotype-dependent allele-specific methylation (hap-ASM), can also contribute to inter-individual differences in DNA methylation patterns. To identify differentially methylated regions (DMRs) associated with breast cancer susceptibility, we examined differences in white blood cell DNA methylation in 29 candidate genes in 426 girls (ages 6–13 years) from the LEGACY Girls Study, 239 with and 187 without a breast cancer family history (BCFH). We measured methylation by targeted massively parallel bisulfite sequencing (bis-seq) and observed BCFH DMRs in two genes: ESR1 (Δ4.9%, P = 0.003) and SEC16B (Δ3.6%, P = 0.026), each of which has been previously implicated in breast cancer susceptibility and pubertal development. These DMRs showed high inter-individual variability in methylation, suggesting the presence of mQTLs/hap-ASM. Using single nucleotide polymorphisms data in the bis-seq amplicon, we found strong hap-ASM in SEC16B (with allele specific-differences ranging from 42% to 74%). These findings suggest that differential methylation in genes relevant to breast cancer susceptibility may be present early in life, and that inherited genetic factors underlie some of these epigenetic differences.     

Genetic and Epigenetic Somatic Alterations in Head and Neck Squamous Cell Carcinomas Are Globally Coordinated but Not Locally Targeted

Background

Solid tumors, including head and neck squamous cell carcinomas (HNSCC), arise as a result of genetic and epigenetic alterations in a sustained stress environment. Little work has been done that simultaneously examines the spectrum of both types of changes in human tumors on a genome-wide scale and results so far have been limited and mixed. Since it has been hypothesized that epigenetic alterations may act by providing the second carcinogenic hit in gene silencing, we sought to identify genome-wide DNA copy number alterations and CpG dinucleotide methylation events and examine the global/local relationships between these types of alterations in HNSCC.

Methodology/Principal Findings

We have extended a prior analysis of 1,413 cancer-associated loci for epigenetic changes in HNSCC by integrating DNA copy number alterations, measured at 500,000 polymorphic loci, in a case series of 19 primary HNSCC tumors. We have previously demonstrated that local copy number does not bias methylation measurements in this array platform. Importantly, we found that the global pattern of copy number alterations in these tumors was significantly associated with tumor methylation profiles (p<0.002). However at the local level, gene promoter regions did not exhibit a correlation between copy number and methylation (lowest q = 0.3), and the spectrum of genes affected by each type of alteration was unique.

Conclusion/Significance

This work, using a novel and robust statistical approach demonstrates that, although a “second hit” mechanism is not likely the predominant mode of action for epigenetic dysregulation in cancer, the patterns of methylation events are associated with the patterns of allele loss. Our work further highlights the utility of integrative genomics approaches in exploring the driving somatic alterations in solid tumors.     

Age-related DNA methylation in normal breast tissue and its relationship with invasive breast tumor methylation

Age is a key risk factor for breast cancer and epigenetic alterations may contribute to age-related increases in breast cancer risk, though the relation of age-related methylation in normal breast tissues with altered methylation in breast tumors is unclear. We investigated the relation of age with DNA methylation in normal breast tissues genome-wide using two data sets from the Gene Expression Omnibus (GEO) database ({"type":"entrez-geo","attrs":{"text":"GSE32393","term_id":"32393"}}GSE32393 and {"type":"entrez-geo","attrs":{"text":"GSE31979","term_id":"31979"}}GSE31979). We validated our observations in an independent set of normal breast tissues, examined age-related methylation in normal breast for enrichment of genomic features, and compared age-related methylation in normal tissue with methylation alterations in breast tumors. Between the two array-based methylation data sets, there were 204 CpG loci with significant (P < 0.05) and consistent age-related methylation, 97% of which were increases in methylation. Our validation sets confirmed the direction of age-related DNA methylation changes in all measured regions. Among the 204 age-related CpG loci, we observed a significant enrichment for CpG islands (P = 8.7E-6) and polycomb group protein target genes (P = 0.03). In addition, 24 of the 204 CpGs with age-related methylation in normal breast were significantly differentially methylated between normal and breast tumor tissues. We identified consistent age-related methylation changes in normal breast tissue that are further altered in breast tumors and may represent early events contributing to breast carcinogenesis. This work identifies age-related methylation in normal breast tissue and begins to deconstruct the contribution of aging to epigenetic alterations present in breast tumors.     

Global hypomethylation and promoter methylation in small intestinal neuroendocrine tumors

Aberrant DNA methylation is a feature of human cancer affecting gene expression and tumor phenotype. Here, we quantified promoter methylation of candidate genes and global methylation in 44 small intestinal-neuroendocrine tumors (SI-NETs) from 33 patients by pyrosequencing. Findings were compared with gene expression, patient outcome and known tumor copy number alterations. Promoter methylation was observed for WIF1, RASSF1A, CTNNB1, CXCL14, NKX2–3, P16, LAMA1, and CDH1. By contrast APC, CDH3, HIC1, P14, SMAD2, and SMAD4 only had low levels of methylation. WIF1 methylation was significantly increased (P = 0.001) and WIF1 expression was reduced in SI-NETs vs. normal references (P = 0.003). WIF1, NKX2–3, and CXCL14 expression was reduced in metastases vs. primary tumors (P < 0.02). Low expression of RASSF1A and P16 were associated with poor overall survival (P = 0.045 and P = 0.011, respectively). Global methylation determined by pyrosequencing of LINE1 repeats was reduced in tumors vs. normal references, and was associated with loss in chromosome 18. The tumors fell into three clusters with enrichment of WIF1 methylation and LINE1 hypomethylation in Cluster I and RASSF1A and CTNNB1 methylation and loss in 16q in Cluster II. In Cluster III, these alterations were low-abundant and NKX2-3 methylation was low. Similar analyses in the SI-NET cell lines HC45 and CNDT2 showed methylation for CDH1 and WIF1 and/or P16, CXCL14, NKX2-3, LAMA1, and CTNNB1. Treatment with the demethylating agent 5-azacytidine reduced DNA methylation and increased expression of these genes in vitro. In conclusion, promoter methylation of tumor suppressor genes is associated with suppressed gene expression and DNA copy number alterations in SI-NETs, and may be restored in vitro.     

Identification of tumor suppressors and oncogenes from genomic and epigenetic features in ovarian cancer

The identification of genetic and epigenetic alterations from primary tumor cells has become a common method to identify genes critical to the development and progression of cancer. We seek to identify those genetic and epigenetic aberrations that have the most impact on gene function within the tumor. First, we perform a bioinformatic analysis of copy number variation (CNV) and DNA methylation covering the genetic landscape of ovarian cancer tumor cells. We separately examined CNV and DNA methylation for 42 primary serous ovarian cancer samples using MOMA-ROMA assays and 379 tumor samples analyzed by The Cancer Genome Atlas. We have identified 346 genes with significant deletions or amplifications among the tumor samples. Utilizing associated gene expression data we predict 156 genes with altered copy number and correlated changes in expression. Among these genes CCNE1, POP4, UQCRB, PHF20L1 and C19orf2 were identified within both data sets. We were specifically interested in copy number variation as our base genomic property in the prediction of tumor suppressors and oncogenes in the altered ovarian tumor. We therefore identify changes in DNA methylation and expression for all amplified and deleted genes. We statistically define tumor suppressor and oncogenic features for these modalities and perform a correlation analysis with expression. We predicted 611 potential oncogenes and tumor suppressors candidates by integrating these data types. Genes with a strong correlation for methylation dependent expression changes exhibited at varying copy number aberrations include CDCA8, ATAD2, CDKN2A, RAB25, AURKA, BOP1 and EIF2C3. We provide copy number variation and DNA methylation analysis for over 11,500 individual genes covering the genetic landscape of ovarian cancer tumors. We show the extent of genomic and epigenetic alterations for known tumor suppressors and oncogenes and also use these defined features to identify potential ovarian cancer gene candidates.     

Hierarchical Clustering of Breast Cancer Methylomes Revealed Differentially Methylated and Expressed Breast Cancer Genes

Oncogenic transformation of normal cells often involves epigenetic alterations, including histone modification and DNA methylation. We conducted whole-genome bisulfite sequencing to determine the DNA methylomes of normal breast, fibroadenoma, invasive ductal carcinomas and MCF7. The emergence, disappearance, expansion and contraction of kilobase-sized hypomethylated regions (HMRs) and the hypomethylation of the megabase-sized partially methylated domains (PMDs) are the major forms of methylation changes observed in breast tumor samples. Hierarchical clustering of HMR revealed tumor-specific hypermethylated clusters and differential methylated enhancers specific to normal or breast cancer cell lines. Joint analysis of gene expression and DNA methylation data of normal breast and breast cancer cells identified differentially methylated and expressed genes associated with breast and/or ovarian cancers in cancer-specific HMR clusters. Furthermore, aberrant patterns of X-chromosome inactivation (XCI) was found in breast cancer cell lines as well as breast tumor samples in the TCGA BRCA (breast invasive carcinoma) dataset. They were characterized with differentially hypermethylated XIST promoter, reduced expression of XIST, and over-expression of hypomethylated X-linked genes. High expressions of these genes were significantly associated with lower survival rates in breast cancer patients. Comprehensive analysis of the normal and breast tumor methylomes suggests selective targeting of DNA methylation changes during breast cancer progression. The weak causal relationship between DNA methylation and gene expression observed in this study is evident of more complex role of DNA methylation in the regulation of gene expression in human epigenetics that deserves further investigation.     

Genome‐wide promoter methylation analysis identifies epigenetic silencing of MAPK13 in primary cutaneous melanoma

The involvement of epigenetic alterations in the pathogenesis of melanoma is increasingly recognized. Here, we performed genome‐wide DNA methylation analysis of primary cutaneous melanoma and benign melanocytic nevus interrogating 14 495 genes using BeadChip technology. This genome‐wide view of promoter methylation in primary cutaneous melanoma revealed an array of recurrent DNA methylation alterations with potential diagnostic applications. Among 106 frequently hypermethylated genes, there were many novel methylation targets and tumor suppressor genes. Highly recurrent methylation of the HOXA9, MAPK13, CDH11, PLEKHG6, PPP1R3C, and CLDN11 genes was established. Promoter methylation of MAPK13, encoding p38δ, was present in 67% of primary and 85% of metastatic melanomas. Restoration of MAPK13 expression in melanoma cells exhibiting epigenetic silencing of this gene reduced proliferation, indicative of tumor suppressive functions. This study demonstrates that DNA methylation alterations are widespread in melanoma and suggests that epigenetic silencing of MAPK13 contributes to melanoma progression.     

Human liver epigenetic alterations in non-alcoholic steatohepatitis are related to insulin action

Both genetic and lifestyle factors contribute to the risk of non-alcoholic steatohepatitis (NASH). Additionally, epigenetic modifications may also play a key role in the pathogenesis of NASH. We therefore investigated liver DNA methylation, as a marker for epigenetic alterations, in individuals with simple steatosis and NASH, and further tested if these alterations were associated with clinical phenotypes. Liver biopsies obtained from 95 obese individuals (age: 49.5 ± 7.7 years, BMI: 43 ± 5.7 kg/m², type 2 diabetes [T2D]: 35) as a wedge biopsy during a Roux-en-Y gastric bypass operation were investigated. Thirty-four individuals had a normal liver phenotype, 35 had simple steatosis, and 26 had NASH. Genome-wide DNA methylation pattern was analyzed using the Infinium HumanMethylation450 BeadChip. mRNA expression was analyzed from 42 individuals using the HumanHT-12 Expression BeadChip. We identified 1,292 CpG sites representing 677 unique genes differentially methylated in liver of individuals with NASH (q < 0.001), independently of T2D, age, sex, and BMI. Focusing on the top-ranking 30 and another 37 CpG sites mapped to genes enriched in pathways of metabolism (q = 0.0036) and cancer (q = 0.0001) all together, 59 NASH-associated CpG sites correlated with fasting insulin levels independently of age, fasting glucose, or T2D. From these, we identified 30 correlations between DNA methylation and mRNA expression, for example LDHB (r = ?0.45, P = 0.003). We demonstrated that NASH, more than simple steatosis, associates with differential DNA methylation in the human liver. These epigenetic alterations in NASH are linked with insulin metabolism.     

The Effects of Nucleoside Analogues on Promoter Methylation of Selected Tumor Suppressor Genes in MCF-7 and MDA-MB-231 Breast Cancer Cell Lines

The effects of 2-chloro-2′-deoxyadenosine, 9-β-D-arabinofuranosyl-2-fluoroadenine, and 5-aza-2′-deoxycytidine on promoter methylation of the selected tumor suppressor genes (i.e., ERα, BRCA1, RARβ2, E-cadherin, PTEN, and APC) were estimated using methylation-sensitive restriction analysis. The studies were carried out in hormone-responsive, low-invasive cell line MCF-7 and hormone-insensitive, highly invasive cell line MDA-MB-231. The results demonstrate an implication of the tested adenosine analogues and 5-aza-dCyd in regulation of DNA methylation process. Moreover, the effects of nucleoside analogues on PTEN promoter methylation suggest distinct mechanism of regulation of the epigenetic DNA modification in low-invasive compared to highly invasive breast cancer cells.     

加勒比松种源苗期DNA甲基化多样性分析

为了解加勒比松(Pinus caribaea)种源的遗传多样性,利用甲基化敏感扩增多态性技术对加勒比松3个变种17个种源的DNA甲基化多样性进行了研究。结果表明,56对引物组合共扩增出425条谱带,其中多态性谱带422条,多态性百分率为99.25%。加勒比松种源幼苗半甲基化比率比全甲基化比率稍高,洪都拉斯加勒比松、古巴加勒比松和巴哈马加勒比松的DNA甲基化率分别为22.39%、22.29%和22.35%,差异不显著。加勒比松的DNA序列遗传多样性(H=0.4376)高于DNA甲基化多样性(H=0.3274),Mantel检验表明,基因组遗传变异与表观遗传变异不存在相关性(r=-0.171,P=0.16)。表观聚类与遗传聚类间存在较大差异,两种聚类分析结果均未将3个加勒比松变种分开。这表明加勒比松变种间的表观遗传变异极为丰富,能为加勒比松遗传改良提供优良种质资源。     

Epigenetic modulation of BRCA-1 and MGMT genes,and histones H4 and H3 are associated with breast tumors

Aberrant patterns in promoter methylation of tumor-suppressor genes and posttranslational modifications of histone proteins are considered as major features of malignancy. In this study, we aimed to investigate promoter methylation of three tumor-suppressor genes (BRCA-1, MGMT, and P16) and three histone marks (H3K9ac, H3K18ac, and H4K20me3) in patients with breast tumors. This case-control study included 27 patients with malignant breast tumors (MBT) and 31 patients with benign breast tumors (BBT). The methylation-specific PCR was used for determining promoter methylation of BRCA-1, MGMT, and P16 genes. Western blot analysis was performed to detect histone lysine acetylation (H3K9ac and H3K18ac) and lysine methylation (H4K20me3). BRCA-1 promoter methylation was detected in 44.4% of the MBT whereas this alteration was found in 9.7% of BBT (P = 0.005). The Kaplan-Meier analysis indicated that hypermethylation in BRCA-1 promoter was significantly associated with poor overall survival of patients with breast cancer (P = 0.039). MGMT promoter methylation was identified in 18.5% of MBT and 0.0% of the BBT (P = 0.01). The frequency of P16 promoter methylation was 25.8% in BBT and 11.1% in MBT (P = 0.12). As compared with BBT, MBT samples displayed the aberrant patterns of histones marks with hypomethylation of H4K20 and hypoacetylation of H3K18 (P = 0.03 and P = 0.04, respectively). There was a negative significant correlation between H3K9ac levels and tumor size in MBT group (r = −0.672; P = 0.008). The present findings suggest that promoter hypermethylation of MGMT and BRCA-1 genes along with alterations in H3K18ac and H4K20me3 levels may have prognostic values in patients with breast cancer. Moreover, the detection of these epigenetic modifications in breast tumors could be helpful in finding new methods for breast cancer therapy.     

Genome-wide combination profiling of copy number and methylation offers an approach for deciphering misregulation and development in cancer cells

Copy number changes and DNA methylation alterations are crucial to gene regulation in mammals. Recently, a number of microarray studies have been based on copy number and DNA methylation alterations in order to find clinical biomarkers of carcinogenesis. In this study, we attempted to combine profiles of copy number and methylation patterns in four human cancer cell lines using BAC microarray-based approaches and we detected several clinically important genes which showed genetic and epigenetic relationships. Within the clones analyzed, many contained cancer-related genes involved in cell cycle regulation, cell division, signal transduction, tumor necrosis, cell differentiation, and cell proliferation. One clone included the FHIT gene, a well-known tumor suppressor gene involved in various human cancers. Our combined profiling techniques may provide a method by which to find new clinicopathologic cancer biomarkers, and support the idea that systematic characterization of the genetic and epigenetic events in cancers may rapidly become a reality.     

The influence of one-carbon metabolism on gene promoter methylation in a population-based breast cancer study

Abnormal methylation in gene promoters is a hallmark of the cancer genome; however, factors that may influence promoter methylation have not been well elucidated. As the one-carbon metabolism pathway provides the universal methyl donor for methylation reactions, perturbation of this pathway might influence DNA methylation and, ultimately, affect gene functions. Utilizing approximately 800 breast cancer tumor tissues from a large population-based study, we investigated the relationships between dietary and genetic factors involved in the one-carbon metabolism pathway and promoter methylation of a panel of 13 breast cancer-related genes. We found that CCND2, HIN1 and CHD1 were the most “dietary sensitive” genes, as methylation of their promoters was associated with intakes of at least two out of the eight dietary methyl factors examined. On the other hand, some micronutrients (i.e., B₂ and B₆) were more “epigenetically active” as their intake levels correlated with promoter methylation status in 3 out of the 13 breast cancer genes evaluated. Both positive (hypermethylation) and inverse (hypomethylation) associations with high micronutrient intake were observed. Unlike what we saw for dietary factors, we did not observe any clear patterns between one-carbon genetic polymorphisms and the promoter methylation status of the genes examined. Our results provide preliminary evidence that one-carbon metabolism may have the capacity to influence the breast cancer epigenome. Given that epigenetic alterations are thought to occur early in cancer development and are potentially reversible, dietary modifications may offer promising venues for cancer intervention and prevention.     

Multi-Platform Whole-Genome Microarray Analyses Refine the Epigenetic Signature of Breast Cancer Metastasis with Gene Expression and Copy Number

Background

We have previously identified genome-wide DNA methylation changes in a cell line model of breast cancer metastasis. These complex epigenetic changes that we observed, along with concurrent karyotype analyses, have led us to hypothesize that complex genomic alterations in cancer cells (deletions, translocations and ploidy) are superimposed over promoter-specific methylation events that are responsible for gene-specific expression changes observed in breast cancer metastasis.

Methodology/Principal Findings

We undertook simultaneous high-resolution, whole-genome analyses of MDA-MB-468GFP and MDA-MB-468GFP-LN human breast cancer cell lines (an isogenic, paired lymphatic metastasis cell line model) using Affymetrix gene expression (U133), promoter (1.0R), and SNP/CNV (SNP 6.0) microarray platforms to correlate data from gene expression, epigenetic (DNA methylation), and combination copy number variant/single nucleotide polymorphism microarrays. Using Partek Software and Ingenuity Pathway Analysis we integrated datasets from these three platforms and detected multiple hypomethylation and hypermethylation events. Many of these epigenetic alterations correlated with gene expression changes. In addition, gene dosage events correlated with the karyotypic differences observed between the cell lines and were reflected in specific promoter methylation patterns. Gene subsets were identified that correlated hyper (and hypo) methylation with the loss (or gain) of gene expression and in parallel, with gene dosage losses and gains, respectively. Individual gene targets from these subsets were also validated for their methylation, expression and copy number status, and susceptible gene pathways were identified that may indicate how selective advantage drives the processes of tumourigenesis and metastasis.

Conclusions/Significance

Our approach allows more precisely profiling of functionally relevant epigenetic signatures that are associated with cancer progression and metastasis.     

Individual and combined effects of DNA methylation and copy number alterations on miRNA expression in breast tumors

Background

The global effect of copy number and epigenetic alterations on miRNA expression in cancer is poorly understood. In the present study, we integrate genome-wide DNA methylation, copy number and miRNA expression and identify genetic mechanisms underlying miRNA dysregulation in breast cancer.

Results

We identify 70 miRNAs whose expression was associated with alterations in copy number or methylation, or both. Among these, five miRNA families are represented. Interestingly, the members of these families are encoded on different chromosomes and are complementarily altered by gain or hypomethylation across the patients. In an independent breast cancer cohort of 123 patients, 41 of the 70 miRNAs were confirmed with respect to aberration pattern and association to expression. In vitro functional experiments were performed in breast cancer cell lines with miRNA mimics to evaluate the phenotype of the replicated miRNAs. let-7e-3p, which in tumors is found associated with hypermethylation, is shown to induce apoptosis and reduce cell viability, and low let-7e-3p expression is associated with poorer prognosis. The overexpression of three other miRNAs associated with copy number gain, miR-21-3p, miR-148b-3p and miR-151a-5p, increases proliferation of breast cancer cell lines. In addition, miR-151a-5p enhances the levels of phosphorylated AKT protein.

Conclusions

Our data provide novel evidence of the mechanisms behind miRNA dysregulation in breast cancer. The study contributes to the understanding of how methylation and copy number alterations influence miRNA expression, emphasizing miRNA functionality through redundant encoding, and suggests novel miRNAs important in breast cancer.