- Open Access
Integrating genetics and epigenetics in breast cancer: biological insights, experimental, computational methods and therapeutic potential
BMC Systems Biology volume 9, Article number: 62 (2015)
Development of human cancer can proceed through the accumulation of different genetic changes affecting the structure and function of the genome. Combined analyses of molecular data at multiple levels, such as DNA copy-number alteration, mRNA and miRNA expression, can clarify biological functions and pathways deregulated in cancer. The integrative methods that are used to investigate these data involve different fields, including biology, bioinformatics, and statistics.
These methodologies are presented in this review, and their implementation in breast cancer is discussed with a focus on integration strategies. We report current applications, recent studies and interesting results leading to the identification of candidate biomarkers for diagnosis, prognosis, and therapy in breast cancer by using both individual and combined analyses.
This review presents a state of art of the role of different technologies in breast cancer based on the integration of genetics and epigenetics, and shares some issues related to the new opportunities and challenges offered by the application of such integrative approaches.
Breast Cancer (BC) is the most common cancer in women and the second most common cause of cancer mortality among females . Classification of BC is currently based on histological types and molecular subtypes in order to reflect the hormone-responsiveness of the tumour. The three most common histological types include invasive ductal carcinoma, ductal carcinoma in situ and invasive lobular carcinoma. The molecular subtypes of BC, which are based on the presence or absence of estrogen receptors (ER), progesterone receptors (PR), and human epidermal growth factor receptor-2 (HER2), include luminal A (ER+ and/or PR+; HER2–), luminal B (ER+ and/or PR+; HER2+), basal-like (ER–, PR–, and HER2–), and HER2-enriched (ER–, PR–, and HER2+) subtypes [2, 3]. This classification reflects the BC heterogeneity and the complexity of diagnosis, prognosis, and treatment of BC.
High-throughput approaches allow today a tumour to be investigated at multiple levels: (i) DNA with copy number alteration (CNA), ii) epigenetic alterations, specifically, DNA methylation, histone modifications and microRNA (miRNA) expression level alterations, and (iii) mRNA, with gene expression (GE) de-regulation. These high-throughput approaches redefined the different types of BC in terms of classification, showing the presence of only two BC profiles with different prognosis [4–6].
Development of human cancer can proceed through the accumulation of genetic and epigenetic changes affecting the structure and function of the genome. Several studies have reported that the epigenetic silencing of one allele may act in concert with an inactivating genetic alteration in the opposite allele, thus resulting in total allelic loss of the gene [7, 8]. Birgisdottir et al.  have reported hypermethylation and deletion of the BRCA1 promoter and suggested Knudson's two 'hits' in sporadic BC . Li et al.  were focused on the expression of beclin 1 mRNA and they demonstrated that loss of heterozygosity and aberrant DNA methylation might be the possible reasons of the decreased expression of beclin 1 in the BC. In BC, a biallelic inactivation of the FHIT gene could be a consequence of epigenetic inactivation of both parental alleles, or epigenetic modification of one allele and deletion of the remaining allele .
In 2006, Feinberg et al. suggested that epigenetics and genetics should be combined or integrated in order to achieve better understanding of cancer . A systems biology approach has been employed to explore the functional relationships among multidimensional “omics” technologies. This approach has been demonstrated to be important for addressing a patient to the optimal treatment in a personalized way, in order to improve the efficacy of the treatment for that patient .
This review refers to current studies of genetic and epigenetic changes associated with BC, focusing in particular on the processes controlled by CNA, epigenetic alterations (DNA methylation, histone modifications and miRNAs), and GE. Several approaches combining genetic and epigenetic data, in particular regarding CNA and miRNA deregulation, have been considered with the final purpose to identify new biomarkers for BC diagnosis and prognosis suitable to be translated into a clinical environment. Furthermore, experimental and computation methods used for the study and the analysis of these biomarkers are presented. We also discuss the biological insights and clinical impact from such analyses as well as the future challenges of these combination approaches.
Copy number alterations in BC
CNAs are alterations of the DNA of a genome that result in a cell having an abnormal number of copies of one or more sections of the DNA. They have been identified as causes of cancer diseases and developmental abnormalities (e.g. ). Changes in DNA copy number (CN) can occur in specific genes or involve whole chromosomes, usually genomic regions between 1kbp and 1Mbp in length .
Figure 1 shows an example of a wild type (WT) cell with two copies of DNA segments that suffer of alterations in tumour cells bringing deletions (CN = 0; CN = 1) or amplifications (CN = 3; CN = 4) of the DNA section.
The ability of cancer cells to accumulate genetic alterations is crucial for the development of cancer in order to inactivate tumour suppressor genes (TSGs) and activate oncogenes (OGs).
In BC, several genetic alterations have been found.
Frequent CN deletions between axillary lymph node metastasis and BC primary tumours were revealed, including aberrations at 6q15-16, containing the gene PNRC1 (a putative tumour suppressor) . Amplification and overexpression of the HER2 (HER2/neu, ERBB2) oncogene on chromosome 17q12 has been observed in 15–25 % of invasive BC . HER2-amplified (HER2+) has been associated with poor prognosis in BC , amplification of the HER2 gene leading to HER2 protein levels 10–100 times greater than normal levels .
EGFR amplification has been frequently associated with indices of poor prognosis in BC patients, such as large tumour size, high histological grade, high proliferative index, HER2 negative, upregulation of PR , and negative ER status .
In the same region of HER2 (17q12–21) other genes have been found co-amplified or deleted, e.g. topoisomerase (TOP2A) . Different studies observed the possibility of guiding therapy based on TOP2A status [22, 23].
A recent study has shown alterations of PIK3CA and MET in BC . High CN of PIK3CA and MET was associated to a poor prognosis, and these alterations occur often in triple receptor negative BC . Alterations were also found at 9q31.3-33.1, where the genes DBC1 and DEC1 (regulators of apoptosis) are located .
OGs activation by genomic amplification occurs in the members of different oncogene families, e.g. MYC and CCND. MYC is a key regulator of cell growth, proliferation, metabolism, differentiation, and apoptosis . This oncogene is located on chromosome 8q24, and several mechanisms are implicated in its deregulation in BC, including gene amplification and traslocations. MYC amplification plays a role in BC progression because it has been detected in the more aggressive phenotype of ductal carcinoma in situ  or in invasive processes [27–29].
Gene amplification of CCND1 has been observed in a subgroup of BCs with poor prognosis and associated with resistance to tamoxifen . Region of amplification is 11q13, and CCND1 acts as a cell cycle regulator, promoting progression through the G1-S phase .
Higher ESR1 gene amplification is found in BC with CCND1 gene amplification in comparison with tumours without CCND1 gene amplification . Amplification of ESR1 has been associated with negative ER . The gene TSPAN1 (on 1p34.1) has been found deleted in metastasizing BC and might represent an important TSG . Another gene, EMSY was found involved in sporadic BC. EMSY amplification has been shown to be associated with a poor prognosis .
Compared to non-metastatic invasive ductal carcinoma, metastatic invasive ductal carcinoma showed a unique pattern of CNAs, including gains at 2p24-13, 2q22-33, 9q21-31, 12q21-23, 17 q23-25 and loses at 11q23-ter, 14q23-31, 20p11-q12, 2q36-ter, 8q24-ter, 9q33-ter, 2p11-q11, and 12q13 [35, 36].
Table 1 reports a synthesis of the considered mutated genes in BC, with their genetic alterations due to CNs.
Current experimental methods for the identification of CNA include cytogenetic techniques, microarrays, and sequencing-based computational approaches.
Karyotyping is a cytogenetic technique performing a standardized and effective single cell screening in order to identify significant genomic aberrations in pathological and in normal samples.
In a standard karyotyping, a dye like Giesma or Quinacrine is used to stain bands on the chromosomes. Each chromosome presents banding pattern for detecting CNAs. Thus, any alteration in banding pattern represents a CNA .
Spectral karyotyping (SKY technique) is a novel technique for chromosome analysis , based on the approach of the fluorescence in situ hybridization technique (FISH). Sky refers to the multicolour-FISH technique where each chromosome is represented with different colours (a dye with different fluorophores). This technique is used to identify CNAs in cancer cells and in other disease conditions when other techniques are not enough accurate .
Resolution is the main limitation of both techniques, the chromosome profile obtained by karyotyping being not enough sensitive to notice short and relevant abnormalities .
Hybridization-based microarray approaches, including array comparative genomic hybridization (array CGH) and Single Nucleotide Polymorphism (SNP) microarrays, have been used as an alternative technology to conventional cytogenetic approaches . They are able to infer CNAs (amplifications and deletions) compared to a reference sample. Array CGH platforms compare quickly and efficiently two labelled samples (different fluorophores - test and reference). Denaturation of the DNA in single stranded allows the hybridization of the two samples to microarrays containing DNA sequence probes of known genome position (e.g. bacterial artificial chromosomes, cDNAs, or more recently, oligonucleotides). By using a fluorescence microscope and a dedicated computer software, the signal ratio of different coloured fluorescents is measured in order to identify chromosomal differences between the two sources. An important consideration is the consequence of the reference sample on the CN profile. A comprehensive-characterized reference is the key for the correct interpretation of array CGH data .
SNP-arrays have a higher resolution than CGH-arrays, and can be used to identify allele-specific information. SNP microarray has few key differences from CGH technologies. Probe designs are specific to single-nucleotide differences between DNA sequences.
Ultimately, next generation sequencing (NGS) have replaced microarrays as the platform for discovery and genotyping, and present considerable computational and bioinformatics challenges.
We can summarize CNA analysis from microarray in three steps: 1) normalization, 2) probe-level modelling, and 3) CN estimation .
The target of normalization is to remove non relevant effects, such as the GC content of the fragment amplified by PCR, technical variations between arrays occurring from differences in sample preparation or labelling, and array production or scanning differences .
Probe-level modelling is usually performed at two levels: single locus and multilocus. Single locus modelling measures the CN of a specific target fragment or DNA probe locus in order to produce a raw fragment CN. Multilocus modelling combines the raw CNs of neighbouring fragments or DNA probe loci into a “meta-probe set” which determines the CN of the whole region [41, 42].
Several methods are suitable for analysing CNA on microarray data.
The first CNA analysis method has been developed by Affymetrix: Chromosome Copy Number Analysis Tool . Normalization is performed by quantile normalization. Modelling uses robust multichip average. CN estimation can be done subsequently with an arbitrary algorithm.
DNA-Chip Analyzer (dChip)  normalizes using an invariant set method which corresponds to a normalization of the arrays based on the identification of a common baseline array and on adjustment of all the other arrays relative to the baseline array. Modelling is based on a model-based expression index (MBEI) for single-locus. This output is then used by a Hidden Markov Model (HMM) to infer CNs .
Copy Number Analyser for GeneChip arrays (CNAG)  normalizes the arrays in order to have the same mean signal intensity for all autosomal probes. This allows fragment probes comparable between arrays to be obtained. The signal intensity ratios is corrected for the differences in PCR product length and GC content. An HMM algorithm is applied to infer CNs along each chromosome.
Birdsuite's Birdseye  normalizes using quantile normalization. Modelling and segmentation are performed together at the multi-loci level. HMM estimates CNs.
Copy-number estimation using Robust Multichip Analysis (CRMA)  has been developed as an extension of the RMA model. Normalization is obtained by allelic cross-hybridization correction (ACC). Modelling uses robust multichip average (RMA). CNA analysis can be done using an arbitrary segmentation algorithm.
Given the different existing computational methods for CNA detection using SNP arrays, researchers have the problem to choose the optimal tool for their analyses.
With the aim of offering a support to bioinformatics researches and to answer to their emerging needs to choose among different CNA detection algorithms, the CNV Workshop was developed . It represents the first cohesive and convenient platform for detection, annotation, and assessment of the biological and clinical significance of structural variants . The purpose of the platform is to process data from a wide variety of SNP arrays, and to implement different normalization and CN estimation algorithms.
Since one of the main problem in the choice of the tool is the detection of discrepancies among different platforms , some studies have compared the different analysis using the same data set. Although limited to few methods, due to the high computational cost, several studies allowed the assessment of advantages and disadvantages of some techniques [49–51].
Baross et al. found that CNAG, dChip, CNAT and GLAD are suitable for high-throughput processing of Affymetrix 100 K SNP array data for CN analysis. However, the tools revealed considerable variations in the numbers of putative CNA. dChip found more CNA than the other tested tools. The highest rate of false positive candidate deletion calls was produced by CNAG. In general, the performance of all tools in the detection of single copy deletions was better than that of single copy duplications. The authors recommend also the use of reference data set for accurate analysis, processed in the same laboratory and ideally from samples with an ethnic composition similar to the sample set.
Eckel-Passow et al.  provided a description of four freely-available software packages (PennCNV, Aroma. Affymetrix, Affymetrix Power Tools (APT), and Corrected Robust Linear Model with Maximum Likelihood Distance (CRLMM)) that are commonly used for CNA analysis of data generated from Affymetrix Genome-Wide Human SNP Array 6.0 platform. APT obtained the best performance with respect to bias. However, PennCNV and Aroma.Affymetrix had the smallest variability associated with the median locus-level CN.
Zhang et al.  assessed four software programs currently used for CNA detection: Birdsuite (version 1.5.2), PennCNV-Affy (a trial version), HelixTree (Version 6.4.2), and Partek (Version, 6.09.0129). They evaluated the accuracy in detecting both rare and common CNVs in the Affymetrix 6.0 platform. They found considerable variations among the programs in the number of CNAs. Birdsuite obtained the highest percentages of known HapMap CNAs containing more than twenty markers in two reference CNA datasets. In the tested rare CNA data, Birdsuite and Partek had higher positive predictive values than the other tools.
Other methods exist for analysing CNA on NGS and they are not described in this review. However, most of the more recent algorithms for CNA discovery are modelled on computational methods which were first used to analyse capillary sequencing reads and fully sequenced large-insert clones .
A future challenging direction is the discovery of gene CN changes for the development of therapies. For example, duplication of one gene encoding a specific receptor can be associated with a particular pathology. Thus, compounds that down regulate receptor expression may lead benefit in patients.
Cancer is the prime case in which CNAs have been shown to drive disease  and therapies where overexpressed or amplified oncogenic drivers are targeted have been already considered. In particular, in BC, the gene encoding epidermal growth factor receptor (EGFR) results to be amplified, and small molecules such as gefitinib, erlotinib, lapatinib, and cetuximab have been applied to inhibit EGFR with benefits for patients [53, 54].
ERBB2, encoding HER2, is amplified in 30 % of BC [17, 55]. In the therapy of HER2-amplified BC, trastuzumab, an anti-HER2 antibody, has been used . Pertuzumab, a humanized monoclonal antibody, binds HER2, and like trastuzumab, it stimulates antibody-dependent, and cell mediated-cytotoxicity . Pertuzumab and trastuzumab binds to different HER2 epitopes acting in the same way. When given together, they operate reinforcing antitumor activity .
These proven benefits, although limited to few genes involved in BC, raise the exciting possibility that targeting amplified disease drivers may offer opportunities for therapy development in BC where effective treatments are still limited.
Epigenetic alterations in BC
DNA methylation and histone modifications
DNA methylation and histone modifications play a crucial role in the maintenance of cellular functions and identity. In particular, the main cellular networks affected by epigenetics are cell cycle, apoptosis, DNA repair, detoxification, inflammation, cell adhesion and invasion.
In cancer, the DNA methylation and histone modifications are perturbed, leading to significant changes in GE, which confer to the tumoral cells advantages in proliferation and maintenance of tumoral phenotype. For instance, the genomic inactivation of a tumor suppressor gene (p53, BRCA1,…) or the activation of an oncogene (i.e., Myc) contribute to the malignant transformation. Epigenetic changes differ from genetic changes mainly because they occur at a higher frequency than genetic changes, they are reversible upon treatment with pharmacological agents and occur at defined regions in a gene.
DNA methylation refers to the addition of a methyl group (−CH3) covalently to the base cytosine (C) in the dinucleotide 5′-CpG-3′. CpGs islands are in the promoter region of many genes [59, 60]. Most CpG dinucleotides in the human genome are methylated, and often leads to silencing of GE. The observation that CpGs islands of housekeeping genes are mainly unmethylated, and the methylation is associated with loss of GE led to the hypothesis that DNA methylation plays an important role in regulating GE [59, 60].
Figure 2 shows how DNA methylation affects GE. Methyl groups in the recognition elements of transcription factors inhibits the binding of transcription factors to DNA, thus resulting in reduced transcriptional activity.
Histones are considered DNA-packaging protein components of chromatin, able to regulate chromatin dynamics. In fact they are subjected to several post-translational modifications, occurring at the amino-terminal end of the histone tail protruding from the surface of the nucleosome . The modifications of histone tails, including lysine acetylation, lysine and arginine methylation, lysine ubiquitylation, phosphorylation, sumoylation, and ribosylation, can significantly affect the expression of genes in a dynamic manner . The most studied histone epigenetic alterations are acetylation/deacetylation, and methylation/demethylation. In BC, abnormal histone modification and DNA hypermethylation are frequently associated to epigenetic silencing of tumor suppressor genes and genomic instability [62, 63].
The distribution of methylated and unmethylated CpGs in the genome shows different patterns of methylation confirming tissue-specific manner .
DNA methylation biomarkers for early detection and prognosis of cancer have been studied in the last years. Table 2 shows genes differentially methylated in BC.
Fackler et al.  found that promoter methylation of 4 genes (RASSF1A, CCND2, TWIST, HIN1) was more frequently detected in tumor than in normal tissue. In another study , 4 genes CCND2, RASSF1A, APC and HIN1 were able to classify between invasive carcinomas, fibroadenomas, and normal tissue. 10 hypermethylated genes, APC, BIN1, BMP6, BRCA1, CST6, ESR-b, GSTP1, P16, P21 and TIMP3, were identified to distinguish between cancerous and normal tissues .
Several studies provide strong evidence of DNA methylation signatures with prognostic role. DNA methylation status of the PITX2 in BC cell lines is negatively associated with PITX2 mRNA expression and with poor prognosis .
Previous studies observed several candidate methylation sites that are associated with the hormone receptor status of BC. RASSF1A and CCND2 were significantly more methylated in the ER+ than ER− BC , whereas the inverse correlations were identified between hypermethylation of the PGR, TFF1, CDH13, TIMP3, HSD17B4, ESR1 and BCL2 genes and ER status . Hypermethylation of the ESR1, TGFBR2, PTGS2 and CDH13 genes was associated with PR status .
Li et al.  used 27 K arrays in a small sample of ER/PR+ and ER/PR BC samples, and identified and validated four genes: FAM124B and ST6GALNAC1 were significantly hypermethylated, and NAV1 and PER1 were significantly hypomethylated in ER+/PR+ BC.
Fang et al.  used genome wide analysis to characterize BCs based on their metastatic potential. The study found a coordinated methylation of a large number of genes discovering a “methylator” phenotype. The methylator phenotype was associated with low metastatic risk and survival.
Identification of promoter methylation of biomarker genes in the DNA of bodily fluids, like serum or plasma, is a rapidly growing research field in cancer detection.
The principle is based on evidence that solid malignant tumors release significant amounts of cell-free DNA into the bloodstream through cellular necrosis or apoptosis .
The analysis of the methylation patterns of cell-free DNA by a blood-based test could become a screening tool. In particular, DNA methylation in circulating free DNA from blood of BC was investigated. ITIH5, DKK3, and RASSF1A promoter methylation from serum were identified as candidate biomarkers for the early detection of BC .
CST6 has been identified by two independent dataset as being differentially methylated between BC and control plasma samples .
SOX17 promoter is highly methylated in primary BCs, in circulating tumor cells isolated from patients with BC, and in corresponding cell-free DNA samples .
Similar studies on plasma identified hypermethylation status of KIF1A , and HYLA2 locus  in BC suggesting methylation level in blood having a power to distinguish very early BC cases from controls.
Non-invasive technique such as blood-test screening is a more suitable and cost-efficient methodology compared to mammography and magnetic resonance imaging.
Actually, in clinical use no specific methylation biomarker has been yet validated, due to the reduced number of matched normal DNA samples in cohorts.
Characterizing more than 880 human BC, Elsheikh et al. have demonstrated that histone acetylation and methylation patterns represent an early sign of BC . Low levels of acethylated lysine and methylated lysine and arginine were described to have prognostic value, i.e. of triple-negative carcinomas and HER2-positive BC subtype [79, 80].
In recent years three major technologies have been employed in DNA methylation analysis: chemical treatment with bisulphite (BS), methylation-specific enzyme digestion, and affinity enrichment [81, 82].
The first category includes an assay to characterize methyl cytosine by treatment of genomic DNA with BS. BS treatment converts unmethylated cytosine residues to uracil, without recognizing methyl cytosine residues, which are protected against this treatment.
Methylated and unmethylated DNA can be distinguished by the employment of sequence analysis (e.g. NGS, microarray). PCR amplicons created after BS conversion can be hybridized to microarrays containing methylation-specific oligonucleotides (MSO; 19–23 nucleotides) to query DNA methylation status . BS-based methods cannot distinguish between methylcitosine and other variants (e.g. hydroxymethylcytosine) .
The second category includes methylation-sensitive restriction endonucleases, which distinguish sequences based on methylation status; furthermore methylcytosine could be identified by immunoprecipitation with antibodies or by affinity purification on methyl-binding protein beads.
Restriction endonucleases and microarray are also combined for high-throughput examination of the methylation status [85, 86]. A limitation in utilizing restriction endonucleases is that enzymes identify only a limited fraction of genome CpG sites [81, 82]. A methodology  with multiple enzyme-mediated restrictions was proposed, leads to a better coverage of all CpG dinucleotides in mammalian genomes.
A third category, enrichment techniques, include methylated DNA immunoprecipitation (MeDIP). Genomic DNA is immunoprecipitated with a monoclonal antibody that specifically identifies 5-methylcytidine. The immunoprecipitated fraction can be detected by PCR in order to identify the methylation state of individual regions .
A combination approach of MeDIP and methylation-sensitive restriction endonucleases was developed, promising to quickly compare methylomes at lower cost . Alternatively, MeDIP can be combined with large-scale analysis (e.g. microarrays) .
Many of the techniques proposed for DNA methylation profiling can be combined with NGS technologies .
Bioinformatics research has been focused on the prediction of DNA methylation information with a dual purpose: i) accurate DNA methylation predictions could replace experimental data, and ii) DNA methylation prediction algorithms from training data can give additional information of an epigenetic mechanism.
A large number of computational predictive models have been developed to identify CpG dinucleotides methylated or unmethylated [91, 92], CpG islands (or CpG-rich regions) methylated or unmethylated [93, 94], and CpG islands (or CpG-rich regions) differentially methylated in different tissue/cell types or phenotypes . Most of them use DNA sequence characteristics combined with a machine-learning algorithm.
Combination approaches of computational and experimental methods can speed up genome-wide DNA methylation profiling and detect crucial factors or pathways driving DNA methylation patterns. However DNA methylation prediction shows some difficulties: i) DNA methylation of the sampled cells need to be averages across cells, ii) there are differences across tissues, iii) DNA methylation can have unstable position, and iv) can be not well located in a genomic locus [96, 97].
A key step for accurate computational predictive models is a correct features selection.
The features can be grouped into two categories: genetic and epigenetic features. Given a region of interest (a CpG island or a genomic region around a particular CpG dinucleotide), the genetic features include: i) general features of the region of interest (e.g., length, and distribution of the CpG dinucleotides in the region), ii) DNA sequence composition of the region of interest, iii) patterns of conserved transcription factor binding sites or conserved elements within or near the region of interest, iv) structural and physicochemical properties of the region of interest, v) functional annotations of nearby genes, vi) single nucleotide polymorphisms of the region of interest, and vii) the conservation of the region of interest among species .
Epigenetic features are also crucial in order to fully characterize DNA methylation status.
DNA methylation, as an epigenetic phenomenon, is affected by some other epigenetic factors, such as histone methylation and histone acetylation.
Statistical methods related to differential DNA methylation data analysis cover a number of different approaches. In particular, these methods are accessible to the user by Bioconductor/R. Table 3 shows some methods and packages currently available for methylation differential analysis, such as Wilcoxon rank sum test (implemented in methyAnalysis package) , t-test (implemented in methyAnalysis, CpGAssoc, RnBeads, and IMA package [99–102]), Kolmogorov-Smirnov Tests (), permutation test (implemented in CpGAssoc package ), empirical Bayes method (implemented in RnBeads, IMA and minfi package [101–103]), and bump hunting method (implemented in bumphunter and minfi package [104, 105]).
Wilcoxon rank sum test detect statistically significant sites according to the absolute difference between the average methylation levels of the analysed groups [106, 107]. This method can have a limitation in case of low or unbalance number of samples groups .
Kolmogorov-Smirnov test is another commonly used test that quantifies distributional differences. However, the Kolmogorov-Smirnov test considers each CpG marker as a sampling unit and its naive application is not valid [110, 111].
Permutation test is a resampling-based nonparametric test which permutes data following the null hypothesis of equal data distributions between groups .
Different number of empirical Bayes models were proposed for differential methylation analysis, with different statistical distribution assumptions . Teng et al.  constructed five empirical Bayes models based on either a gamma distribution or a log-normal distribution, for the detection of differential methylated loci. They observed that log-normal, rather than gamma, could be a more accurate and precise method.
Bump hunting method used in bumphunter and minfi packages based the correlations of methylation levels between nearby CpG locus, and, for each locus, a linear model was used to estimate the coefficient of difference in methylation levels between the cancer group and the normal groups [105, 107].
A comparison study among these six statistical approaches was proposed . Finally, different approaches were recommended for different applications: the bump hunting method is better for small sample size; the empirical Bayes methods are suggested when DNA methylation levels are independent across CpG loci, while only the bump hunting method is suggested when DNA methylation levels are correlated across CpG loci. All methods are found suitable for medium or large sample sizes .
Cancer was the first group of diseases to be associated with DNA methylation. Numerous genes have been identified as being differentially methylated in BC, with a crucial roles in DNA repair, apoptosis, hormone receptor, and cell cycle. These TSGs may be good therapeutic targets through regulation of methylation activity by DNA methyltransferase inhibitors. Human DNA methylation is catalysed by enzymes of the DNA cytosine methyltransferases family including DNMT1, DNMT3A, DNMT3B and DNMT3L . A lower DNA methyltransferase activity increases expression of silenced genes such as TSGs reactivating expression of key genes.
Key targets for potential DNA demethylation agents are DNA methyltransferase inhibitor 5-aza-2′-deoxycytidine (decitabine), zebularine, and SGI-110 . The mechanism of action of these pro-drugs is similar since they need to be incorporated into DNA to act as inhibitors of DNMTs .
Decitabine shows activity against hematologic malignancy and low-dose correlates with changes in GE induced by a reduction in DNA methylation.
A phase I clinical and pharmacodynamic trial was proposed in order to assess the feasibility of delivering a dose of decitabine combined with carboplatin . Decitabine showed some limitations for treatment of advanced solid tumors (e.g. BC): i) weak stability, ii) lack of specificity for cancer cells, and iii) rapid inactivation by the action of cytidine deaminase .
Zebularine and SGI-110 are more selective for cancer cells and have higher resistance to deamination. In particular, Zebularine  showed an antitumor effect in a mouse model. In zebularine-treated mices, the oral treatment with zebularine showed a significant delay in tumor growth . In combination with decitabine, zebularine has proven a significant inhibitory effect on cell proliferation and colony formation in MDA-MB-231 BC cell line through induction of ER alpha and PR mRNA expression . Unfortunately, toxicity remain its main limitation .
SGI-110 , a 5′-AzapG-3′ dinucleotide, induces expression of the p16 tumor suppressor gene, and inhibit tumor cell growth. This short oligonucleotide is resistant to cytidine deaminase deamination which may potentially increase its resistance, enhance bioavailability, and make the drug more efficacious.
DNA methyltransferase inhibitors can have side effects as the concomitant activation of both TSGs and OGs. The combination of chemotherapeutic agents and of DNA methyltransferase inhibitors could be efficacious .
Although the benefits of DNA methyltransferase inhibitors were demonstrated, toxicity, lack of specificity and low stability are issues to be solved in order to improve BC treatment .
Histone acetylation process is controlled by the balanced activity of histone acetyltransferases and histone deacetylases (HDACs). The HDAC family is divided into zinc-dependent enzymes (classes I, IIa, IIb, and IV, of which there are 11 subtype enzymes) and zinc-independent enzymes (class III, also called sirtuins), requiring NAD+ for their catalytic activity. Over the past decade, a number of HDAC inhibitors have been designed and synthesized, based on HDAC chemical structures. Some of these HDAC inhibitors are able to modify the chromatin structure, causing re-expression of aberrantly silenced genes, which in turn is associated with growth inhibition and apoptosis in cancer cells . In ER-negative BC, the treatment with specific HDAC inhibitors reactivates ERα and progesterone receptor (PR) gene expression, which are known to be aberrantly silenced in BC. Preclinical studies of HDAC inhibitors combined with DNMT inhibitors or with anti-tumoral treatment (i.e., tamoxifen) have demonstrated a higher safety, tolerability and clinical effectiveness than single treatment [123, 124].
microRNA deregulation in BC
miRNAs are small noncoding RNAs (20–22 nucleotides long) that are excised from longer (60–110 nucleotides) RNA precursor [105, 106] and act in different biological functions including development, proliferation, differentiation and cell death [125, 126]. miRNAs are major regulators of GE. Many evidences indicate that their deregulation is associated to several steps of cancer initiation and progression. In comparison with other approaches targeting single genes, they are certainly more stable thanks to their small size , and are able to discriminate different BC subtypes.
Blenkiron et al. found deregulated miRNAs between basal and luminal BC . Iorio et al. , Lowery et al.  and Mattie et al.  identified miRNAs that were able to classify ER, PR and HER2/neu receptor status, respectively. Gregory et al.  found miR-200 associated with the BC luminal subtype. Reduced expression of miR-145 and miR-205 was found to play a role in basal like triple negative tumours (ER-/PR-/HER2-) while are normally expressed in normal myoepithelial cells .
miRNAs can be also prognostic and predictive biomarkers. Zhou et al.  found mir-125b as useful indicator for poorly response to taxol-based treatments in vivo. The overexpression of miR-181a has been correlated with lymph node metastasis . miR-106b-25 expression was proven significantly predictor of good relapse time , while miR-375 was found negatively regulate ER expression .
miRNAs with a role in metastasis in BC include miR-7 [138, 139], miR-17/20 [140, 141], miR-22 [142–144], miR-30 [145, 146], miR-31 [147–149], miR-126 , miR-145 , miR-146 , miR-193b , miR-205 , miR-206 , miR-335 , miR-448 , miR-661  and let-7 .
miRNAs can be easily extracted and detected from blood , circulating exomes , saliva [162, 163], and even sputum [164, 165]. Several studies demonstrated that circulating miRNAs reflect the pattern observed in the tumour tissues (e.g. ), thus opening the possibility to use circulating miRNAs as biomarkers for diagnosis and prognosis. Lodes et al.  provided an evidence on using serum miRNAs as biomarkers to discriminate between normal and patients in many cancer diseases including breast, prostate, colon, ovarian, and lung cancer. They showed that it is sufficient 1 mL of serum to detect miRNA expression patterns, without the need of amplification techniques. Recently, an analysis of circulating miRNAs have led to identify mir-21, miR-92a [167, 168], miR-10b, miR-125b, miR-155, miR-191, miR-382  and miR-30a  as candidate biomarkers for early detection of BC. Circulating miRNAs have also been associated with disease prognosis and response to treatment. Madhavan et al.  found circulating miRNAs as marker of disease free survival and overall survival. Plasma miR-10b and miR-373 were found associated with the development of metastases  while miR-125b  and miR-155  have been found correlated to chemotherapy response.
Table 4 reports a synthesis of the considered miRNAs deregulated in BC, with their principal biological effects.
Many technologies for detecting miRNAs have been developed, including RT-PCR, in situ hybridization, microarray, and NGS .
RT-PCR is a sensitive and precise technology but it is also an expensive and low-throughput method .
In situ hybridization is based on labelled complementary strands for the sequences of interest (e.g. miRNA) in a portion or section of tissue . The small size of the mature miRNA presents problems for conventional in situ hybridization methods and it is semi-quantitative.
Microarrays have several limitations as those due to background or cross-hybridization problems. Moreover, microarrays and other techniques can provide analyses only on known miRNAs .
Contrarily, sequence-based methods allow the identification of unknown miRNAs and early overcome other methods. Stark et al. , by using deep sequencing, discovered and quantified new miRNAs. Similarly, Farazi et al. , generated a miRNA signature able to differentiate ductal breast carcinoma in situ, invasive ductal breast carcinoma and normal tissue.
Also deep sequencing may be a powerful method to study circulating miRNAs. Several studies investigated correlations among miRNAs in the serum of BC with clinicopathological indices  and found miRNAs associations with overall survival .
Despite the high potential and promising results of these methods in clinical applications, there are still some problems that need to be addressed, e.g. the lack of inconsistency for some results between different studies. Standardization of procedures for sample conservation, preparation and/or processing , and the use of different quality controls for data normalization  could be effective in reducing these limitations.
There are two different approaches to examine both miRNA and mRNA expression profiles.
A first approach considers either a miRNA or an mRNA first, and then applies ad-hoc strategies, such as computational or experimental methods, in order to obtain miRNA-mRNA pair information [182, 183].
Computational methods play important roles in the identification of new miRNAs. These methods can be divided into three major categories: 1) sequence or structure conservation-based, 2) machine learning-based method, 3) and non-comparative methods.
Sequence or structure conservation-based methods are based on sequence/structure conservation as techniques to find miRNAs. The principle is the nature of conservation across different species for most of the known miRNAs. Comparative genomics filter out sequence/structure conservation that are not evolutionarily conserved in related species . Examples of such computational methods, focusing on the secondary structure of RNA and looking for conserved hairpin structures between related species, are Srnaloop , MiRscan , and miRseeker . One of the first study related to these methods was by Lee and Ambros . The authors, using bioinformatics techniques, searched for sequences conserved between the C. elegans and C. briggsae genomes. They focused on premiRNA sequences and secondary structures with similar characteristics to lin-4 and let-7, the first two miRNAs found on that time.
Several web based software tools have been developed to find new miRNA genes, based on sequence and secondary structure similarities with known miRNAs [189–191, 193]. However, the limit of these approaches was demonstrated by Bentwich et al., showing the possibility that large quantity of nonconserved miRNAs could be missed by the use of this tool .
Free energy (or Gibbs free energy) can be used as feature for miRNA target prediction. It shows how strong the binding of a miRNA with its target is by predicting how the miRNA and its candidate target will hybridize. The free energy of miRNA-mRNA binding is normally assigned by the RNAfold program-Vienna RNA Package .
Machine learning-based methods do not necessarily depend on sequence conservation. A classifier is constructed on a training dataset, that contains a set of known miRNA sequences (positive training dataset), and on a set consisting of mRNAs, tRNAs and rRNAs (negative training dataset). The information given to the classifier can be, for instance, the position of the mature sequence or the folding energy. The classifier, by describing a candidate miRNA with this set of features, is able to predict true and false miRNA sequences . The limit of this approach is the choice of negative set. As example, we do not know a priori if a particular sequence can generate functional miRNAs . Several studies have tried to overcome this kind of problems with the use of only positive models [198, 199]. However, the results were poorer than those found by approaches that consider both positive and negative training sets .
Different classification methods are currently available based on machine learning, e.g. SVM, neural networks, HMM, and Naive Bayes (NB), and several tools based on machine learning have been developed and released to the research community, e.g. RNAmicro , MiRFinder , ProMir , MiRRim , SSCprofiler , HHMMiR  and BayesMiRNAFind .
Non-comparative methods use intrinsic structural features of miRNA, and include algorithms like PalGrade , Triplet-SVM , miPred , miR-abela , and HHMMiR . These methods are able of detecting a large number of miRNAs that seem to be unique to primates.
Bentwich et al.  developed PalGrade by integrating bioinformatics predictions with microarray analysis and sequence-directed cloning. This approach allowed the detection of 89 human miRNAs, 53 of them being not conserved beyond primates.
Xue et al.  proposed an ab initio classification of real pre-miRNA from other hairpin sequences with similar stem-loop features. SVM was applied on these features to classify real vs pseudo pre-miRNAs achieving 90 % of accuracy.
Ng et al.  employed a Gaussian Radial Basis Function kernel (RBF) as a similarity measure for 29 global and intrinsic hairpin folding attributes. They tested the model on 123 human pre-miRs and 246 pseudo hairpins, reporting 84.55 %, 97.97 %, and 93.50 % in sensitivity, specificity and accuracy, respectively.
Sewer et al.  developed miR-abela to detect human miRNAs. They focused on particular properties of some genomic regions around already known miRNAs, and were able to predict between 50–100 novel pre-miRNAs, 30 % of them already found as new in other studies.
miRNAs may have a crucial role in guiding treatment decisions. miRNAs can be therapeutic agents in cancer for two major characteristics: (1) their expression is deregulated in cancer compared to normal tissues, and (2) cancer phenotype can be changed by targeting miRNA expression .
Compared to gene profiles, miRNA-based therapeutics have several advantages, as for example their ability to target multiple genes, frequently in the context of a network. miRNAs regulating the network of genes and cellular pathways play a crucial role in BC pathogenesis and therapy.
There are two strategies for developing miRNA-based therapies: i) by the introduction of miRNA-mimic oligonucleotides, which mimic miRNA expression, up-regulating miRNA, and ii) by the introduction of miRNA inhibitor oligonucleotides to inhibit the expression of the miRNA of interest. However, some major obstacles for the use of miRNA therapeutics exist, including the tissue-specific delivery [211, 212], and the fact that erroneous targeting of miRNAs may cause toxic phenotypes .
For an effective drug-design of miRNA-targeted therapies in BC, it could be useful to understand the interplay between miRNAs and mRNAs leading to BC, thus studying the networks of gene controlled by each miRNA of interest. miRNA and their targets can form complex regulatory networks, and the comprehension of miRNA-target relation will help the development of personalized and tailored therapies .
Gene expression deregulation in BC
GE profiling in BC has been widely demonstrated to generate different prognostic and diagnostic gene signatures. However, molecular tests have a potential not only for diagnosis but also for tailoring treatment plans, in particular with the aim of reducing resistance, non-response and toxicity . Most of the tests either focus on gene expression microarrays or quantitative reverse transcription (qRT)-PCR analyses.
van't Veer et al.  obtained one of the prognostic signature for BC currently available on the market: MammaPrint. Microarray analysis of 78 BC patients with no systemic therapy led to the identification of a list of 70 genes able to predict the prognosis of the disease. The test was independently validated in a cohort of 295 early stage invasive BC, and results proved that the signature was an independent prognostic marker in BC . A second independent validation study was performed by the TRANSBIG Consortium  in a cohort of 302 adjuvantly untreated patients, and was followed by additional validation studies [218–220]. MammaPrint was developed by Agendia, a laboratory in Amsterdam, approved in 2007 by the U.S. Food and Drug Administration (FDA) and then released commercially. This is a microarray-based test assessing the risk that a BC can metastasize to other parts of the body.
Paik et al.  developed Oncotype DX, a qRT PCR-based signature which measures the expression level of 21 genes (16 target + 5 reference genes). The test is able to predict chemotherapy benefits and the likelihood of distant BC recurrence. This is the first genomic biomarker assay which is commercially available for BC treatment as support of chemotherapy. Three separate studies containing 447 BC patients allowed to identify the 21-gene profile, which were divided into 16 target and 5 reference genes. The test was then validated using 668 node negative, ER positive, tamoxifen treated patients from NSABP B-14. An Oncotype DX Recurrence Score (ODRS) was defined and measured, as expression of a risk percentage for the development of distant metastases . Oncotype DX was subsequently evaluated in the NSA BP-B20 trial, a study that explored the benefit of chemotherapy plus tamoxifen, and proved the accuracy of the biomarker. Currently, the Oncotype DX assay is performed in the licensed Genomic Health laboratory, which is the laboratory where the assay was developed.
Prediction Analysis of Microarray (PAM50), by using qRT-PCR assay, measures the expression of 55 (50 target and 5 reference genes) to identify the intrinsic subtypes of BC: luminal A, luminal B, HER2-enriched, and basal-like . The gene signature was developed by analysing 189 BC samples, and was then validated on 761 BCs for prognosis and on 133 BCs for prediction of response to a taxane and anthracycline regimen . NanoString’s Prosigna™ received a CE-mark designation for selling BC PAM50 in 2012, and received FDA clearance in 2013.
Genomic Grade Index (GGI)  is a 97 gene which measures the histological tumour grade. This test is based on the assumption that histological grade is a strong prognostic factor in ER positive BC. Sotiriou et al.  found that GGI gene signature is able to classify BC as histological grades I and III. They used 64 samples of ER-positive BC tumours to select genes that were differentially expressed (DE) between histologic grade I and III tumours, and to generate the gene signature. Data from 597 independent tumours were then used to evaluate GGI and to also demonstrate that GGI can separate histological grade 2 BC into low or high categories with different clinical outcomes. The results of the BIG-1-98 study (55 endocrine-treated patients)  demonstrated that the GGI is also a potential predictor of relapse for endocrine-treated BC patients. Ipsogen launched the MapQuant Dx (TM) genomic grade test by incorporating GGI. The test is currently used, in particular when tumor grade information can be decisive for prescribing a chemotherapy.
Immunohistochemical (IHC) assay Mammostrat  uses 5 immunohistochemical markers (SLC7A5, HTF9C, P53, NDRG1, and CEACAM5) to stratify patients on tamoxifen therapy into different risk groups, in order to inform treatment decisions. In the validation study, an analysis was performed on two independent data sets of 299 and 344 BC samples . Clarient launched on the market the Insight® Dx Mammostrat® Breast Cancer Recurrence test in 2010.
Table 5 reports the considered commercially available tests, with their principal characteristics (e.g. number of genes, validation data sets).
Understanding GE and how it changes under normal and pathological conditions is necessary to provide information about the expressed genes. Large scale GE data provide the activity of thousands of genes at once.
Several techniques exist for studying and quantifying GE.
Traditional methods focus on measuring the expression of one gene at a time, as, for example, the Northern Blotting and the Real-Time Quantitative Reverse Transcription PCR (RT-PCR).
Northern blotting (called also RNA blot) was the first tool used to measure RNA levels, and, until the end of the 1990s, it was used extensively. It allows to quantify levels of mRNA by electrophoresis, which is able to separate RNA samples by size. The RNA of interest is revealed by a hybridization probe complementary to it. The first step of RNA blot is to denature the RNA into single strands. Hence, gel electrophoresis separates the RNA molecules according to their size. Subsequently, the RNA is transferred from the gel onto a blotting membrane, containing RNA bands originally on the gel. A probe complementary to the RNA of interest binds to a particular RNA sequence in the sample . The RNA-probe complexes can thus be detected using a variety of different chemistries or radionuclide labelling.
RT-PCR is a major development of PCR technology, overcoming Northern blot as the method for RNA detection and quantification . It enables to monitor and measure a targeted DNA molecule generated during each cycle of PCR process. In RT-PCR, the mRNA must be converted to a double-stranded molecule by using the enzyme reverse transcriptase. This phase is followed by quantitative PCR (qPCR) on the cDNA with the detection and quantification of amplified products . The quantity of each specific target is obtained by measuring the increase in fluorescence signal from DNA-binding dyes or probes, during successive rounds of enzyme-mediated amplification. The limitation of this technique is the quantification of few genes at a time .
Several technologies such as microarray, Serial Analysis of Gene (SAGE), Cap Analysis of Gene Expression (CAGE) and Massively Parallel Signature Sequencing (MPSS) allow the mRNA expression data for hundreds of genes to be obtained in one single experiment .
The most commonly used technology to profile the expression of thousands of transcripts simultaneously is microarray. DNA microarray is an array of oligonucleotide probes bound to a chip surface [231, 232]. Labelled cDNA from a sample is hybridized to complementary probe sequences on the chip, and strongly associated complexes are identified by detection of fluorophore-, silver-, or chemiluminescence-labelled targets [231, 232].
Many variables influence the outcome of the experiments in microarray analysis, thus contributing to experimental errors and biological variations (for more details see ).
In contrast to microarray methods, sequence-based approaches directly determine the cDNA sequence . SAGE , CAGE , and MPSS , all tag-based sequencing approaches, are based on Sanger sequencing technology.
The development of novel high-throughput DNA sequencing methods, such as RNA-Seq (RNA sequencing), has provided new approaches for both mapping and quantifying transcriptomes. It has clear advantage over existing approaches: RNA-Seq is not limited to the detection of transcripts that correspond to existing genomic sequence, and it is suitable to discovery genomic sequences that are still unknown .
In RNA-Seq analysis, RNA is converted to a library of cDNA fragments with adaptors attached to one or both ends. Each molecule is then sequenced in a high-throughput way in order to obtain short sequences (reads 30–400 bp). Following sequencing, the resulting reads are mapped to the genome in order to produce a genome-scale transcription map consisting of both the transcriptional structure and the level of expression for each gene . Although RNA-Seq has many advantages with respect to the other methods, other issues must be overcome to achieve best practices in the measurement of gene expression, for instance, the lack of accurate methods able to identify and track the expression changes of rare RNA isoforms from all genes .
Table 6 reports a synthesis of the considered experimental methods for studying and quantifying GE, with their principal advantages and limitations.
Microarray or RNA-sequencing technologies, as above reported, produce an overall design of all the transcriptional activity in a biological sample. However, these methods necessarily produce a large amount of data to be visualized, evaluated for their quality, normalized, filtered and interpreted.
Hence, the data originated by platforms (such as microarrays or RNA-seq) must be pre-processed. Pre-processing step is crucial to normalize the data and to clean biological signal values from experimental noise [238, 239].
Data must be also reduced prior to be used in advanced analysis, and this can be accomplished in two different ways: 1) by dimensionality reduction methods, that do not modify the original representation of data, and 2) by dimensionality reduction techniques which involve modification or loss of information from the original data. Among this second category, there are those methods based on projection (e.g. principal component analysis) or compression (e.g. using information theory) .
One of the most validated method of the first category is feature selection technique. It is often used to identify key genes able to separate the samples into different classes (e.g. cancerous and normal cells), and to remove irrelevant genes. Golub et al.  showed indeed that most genes are not significant in a problem of samples classification. However, feature selection is also important in order to obtain faster and efficient classification models, and to avoid over fitting.
Filter methods find subset of genes dependent on the class label, and do not consider the relevance of genes in combination with other genes . Usually they are simple and fast.
Filter methods include correlation-based feature selection (CFS) , t-test [243, 245], information gain [243, 245], mutual information , entropy-based methods , Euclidian distance , signal to noise ratio , and significant analysis of microarrays .
Wrapper methods try to achieve the best combination of genes that may offer high classification accuracy. They include hybrid genetic algorithms , particle swarm optimization , successive feature selection (SFS)  and GA-KDE-Bayes . However, this approach is less used, in particular in microarray analysis, due to its high computational costs .
Filter approach does not interact with the classifier, contrarily to wrapper and embedded techniques, usually resulting in lower performance.
An intermediate approach between the lowest results of the filter methods and the high computational cost of the wrapper methods is represented by the embedded method. With this method, the feature selection procedure is inbuilt to a classifier. Classification trees like ID3, random forest, and Support Vector Machine (SVM) based on Recursive Feature Elimination (RFE) are all examples of embedded methods [243, 245].
The principle of the feature selection and validation techniques is shown in Fig. 3. A pre-processing step is performed: i) the quality of data is evaluated, ii) outliers are removed, and iii) data are normalized. Feature selection is performed. Usually, original data are divided into two data sets: a training set, subjected to the feature selection, and a testing set, used to evaluate the feature selection of the model with different validation techniques. Feature selection finds a subset of genes of interest, (e.g. a gene signature), and the validation of genes is performed. The most used validation techniques are cross-validation or leave one out validation , even if several studies suggested the use of a 10-fold cross validation because they give a more biased but less variable estimate than the leave-one-out error (e.g. ). When the feature selection of the model satisfies the required validation performance, the genes are defined and can be interpreted.
Drug compounds that facilitate and control tightly therapeutic GE are a promising target. Transcriptional gene regulatory system has been encoded within several viral vectors (eg. Tetracycline-based systems can regulate GE of particular targets with the use of cell-type-specific promoters) .
The regulation of GE systems is an attractive target for gene therapy development, and potential applications have been assessed in a wide variety of preclinical laboratory models of disease. The first study was performed by Hallahan et al. , which described how TNF-a expression, under the control of the Egr-1 promoter, could be increased in response to ionizing X-ray radiation. This increase of TNF-a expression was associated with an improved control of tumour growth in comparison with X-ray radiation alone . Advantages from induction of GE by ionizing radiation include reduction of damage to adjacent healthy tissues .
Kan et al.  have constructed a novel retroviral vector (MetXia-P450) encoding CYP2B6. This vector was used to transfect the human tumour cell lines HT29 and T47D. CYP2B6 metabolizes the prodrug cyclophosphamide (CPA) to produce phosphoramide mustard that cross-links DNA, thus leading to cell death. In order to evaluate safety and clinical response, MetXia-P450 entered Phase I clinical trials for nine BC patients and three melanoma patients with cutaneous tumours, with encouraging results.
Although viral vectors are very efficient for gene transfer, their uses are still limited by safety concerns . As an alternative, non-viral BC gene therapy (e.g. naked DNA) is growing due to its safety profile, easy preparation procedures, and moderate costs. β-galactosidase (LacZ) expressing plasmid DNA has been successfully delivered in three patients by a needle-free jet injection to skin metastases from primary BC, and also to melanoma lesions in 14 patients. No side effects were observed. The transgene was detectable at messenger RNA (mRNA) and at protein levels in all patients.
Copy number alterations and gene expression in BC
Several studies demonstrated that changes in DNA CN are translated into corresponding changes in GE [258, 259]. Although it is possible that changes in specific DNA sequences (i.e. centromeres or telomeres) can have directly negative consequences , the main responsible for the malignant phenotype has been proven to be the gene dosage hypothesis: alterations of gene copies change the expression levels of the involved gene .
Figure 4 shows the principle consequences of an altered gene dosage. Specifically, figure 4.1) shows: i) WT condition where a correct number and expression of A and B gives a correct production of C; ii) how the amplification/over expression of gene copies (e.g. B) can cause an increased dosage of a single gene (e.g. C), and iii) how a deletion/under expression of gene copies (e.g. B) can cause a decreased dosage of a single gene (e.g. C) . Figure 4.2) shows how altered gene dosage can influence stoichiometry of protein complex DE that produces F. An amplification/over expression of protein D can inhibit the formation of protein complex DE, thus altering the pathway activity and the correct production of F. A deletion/under expression of protein D do not produce protein complex DE .
While useful information has been revealed by analysing GE profiles alone or CNA data alone, integrative analysis of CNA and GE data are necessary in order to have more information in gene characterization. Specifically, RNA data give information on genes that are up/down-regulated, but do not consider primary changes driving cancer from secondary modifications, such as proliferation and differentiation state. On the other hand, DNA data give information on amplifications and deletions that are drivers of cancer. Therefore, integrating DNA and RNA data can clarify genetic regulatory relationships in cancer cells . It is interesting that transcriptional changes for 10–63 % of genes occur in amplified regions, and, for 14–62 % of genes, in regions of loss .
Several studies showed that gains (or losses) in DNA genomics have consequences in the expression levels of genes in the implicated regions, which are increases or decreased, respectively [264–266]. If we consider individual genes, the situation is more complicated. For instances, 14 % of down-regulated genes can appear within regions of DNA gain, while 9 % of up-regulated genes can occur in regions of DNA loss . These findings suggest to take a particular attention in the integration of CNA and GE.
The Cancer Genome Atlas project  is generating multidimensional platforms including gene expression and CNA data for the same set of patients . Although it is possible to perform analysis with unpaired data [263, 268, 269], the analysis is much more accurate when both types of data are available from the same patient. In this condition, the paired data analysis allows better statistical power and a reduction of false positives [270, 271].
Some studies have shown that integrating CNA information with GE data can often provide a powerful tool for identifying functionally relevant genes in cancer [e.g. 275–282]. Chen et al.  found a list of eighteen genes for which a strong correlation between CNA and GE exists, using signal-to-noise ratio (SNR). They found one particular gene, RUNX3, which is involved in the control of the in vitro invasive potential of MDA-MB-231.
Zhang et al.  identified an 81-gene prognostic CN signature that was found highly correlated with GE levels (Cox regression P < 0.05). This signature identified a subgroup of patients with increased probability of distant metastasis in an independent validation set of 113 patients.
Andre et al.  reported the level of mRNA expression, significantly correlated to the CAN, for VEGF, EGFR, and PTEN, using Algorithm Array CGH Expression integration tool (ACE-it). These genes could be targeted in triple-negative BC in clinical trials, and one of them, E2F3, can have a major role in a subset of triple-negative BC.
Hyman et al.  studied CNAs in 14 BC cell lines, and identified 270 differently expressed genes using signal-to-noise statistics (α value <0.05). 91 of the 270 genes represented hypothetical proteins or genes with no functional annotation, whereas 179 genes had available functional information.
Orsetti et al.  presented a study on CNA on chromosome 1, the prevalent target of genetic anomalies in BC, and the CNA consequences at the RNA expression level in BC. They identified 30 genes showing significant over-expression. A discriminating score was applied by comparing the expression levels of the subgroup of samples presenting amplification and the expression levels of the subgroup of samples without amplification.
Chin et al.  associated CNA and GE profiles of genes linked to poor treatment response. They identified 66 genes in these regions whose expression levels were correlated with CN, using Pearson's correlation (FDR < 0.01, Wilcoxon rank-sum test). Gene Ontology analyses of these genes showed that they are involved in nucleic acid metabolism, protein modification, signalling, and in the cell cycle and/or protein transport.
Chin SF et al.  evaluated genome-wide correlations between GE and CN by following an approach based on the Wilcoxon test. They showed strong statistical associations between either CN gain and over-expression (196 genes) or CN loss and under-expression (63 genes). Many well-known and potentially novel oncogenes and tumour suppressors were included in their analysis.
Table 7 reports a synthesis of the considered genes based on the integration of CNA and GE.
No experimental methods actually exist giving, in one single analysis, results about the integration of CNA and GE.
Computational integrative methodologies between CNA and GE include a two-step approach, and joint analysis. Figure 5a) shows a two-step approach, combining the results from individual analysis of GE and CNA. Figure 5b shows a joint analysis obtaining directly the final result from the integration of GE and CNA.
There are different statistical measures to assess the CNA and GE relationship in order to quantify gene dosage effect. They include, in two-steps approaches, both regression and correlation-based analysis.
Regression approaches model the dependence of RNA levels from DNA CN, and consider RNA levels as responses and DNA CN as predictors . These methods can be divided into: 1) univariate linear regression models, proposed to model the associations between individual CN and GE probes , 2) multivariate linear regression models, integrating statistical power across multiple probes targeting adjacent genes or chromosomal positions , and 3) nonlinear regression models.
Most studies use linear regression models but regulatory mechanisms, contributing to gene expression changes (e.g. CNA, miRNA, DNA methylation), can give non linearity . Non linear relationship between CNA and GE have been investigated by Solvang et al. , which focused on the identification of nonlinear relationships to explain the regulatory mechanisms of alteration of mRNA expressions in the cancer process.
Correlation-based approaches have been used to study the relationship between CNAs and GE. For each pair of co-measured data, a correlation matrix was estimated reflecting the strength of association . Several studies have shown correlations between CNA and GE gene across samples e.g. . Other studies, like Tsafrir et al. , identified a correlation along the genome by using filtered CNA and filtered GE. DR-Correlate , a modified version of the Ortiz-Estevez algorithm,  was used in a correlation-based analysis to examine the genome and to detect genes with high associations between CNA and GE. In order to improved correlation results, Schäfer et al.  replaced the sample means with the reference medians in the correlation test, while Lipson et al.  used a quantile-based analysis to obtain improved correlation coefficients.
Table 8 reports a synthesis of the considered two-step analyses and types.
Joint analysis uses CNA and GE data as paired data entries and not as separate structures. The discrepancy between the sample size and the number of genes is a problem that can cause high noise. Techniques such as Singular value decomposition (SVD) or Principal Component Analysis (PCA) are the most popular ones for reducing the dimension of gene data [289, 290]. However, GE and CNA data are separately analysed using these methods.
The generalized singular value decomposition (GSVD) is a popular regression framework used in joint analysis. With the purpose to identify variation patterns between two biological inputs, Berger et al.  applied an iterative procedure based on the GSVD, projecting CNA/GSE data into different decomposition directions. GSVD was used in two BC cell lines and tumour datasets, thus obtaining gene subsets that were biologically validated.
Soneson et al.  applied PCA to reduce dimensions, and Canonical Correlation Analysis (CCA) to identify highly correlated CNA/GE pairs. Gonzalez et al.  implemented the regularized CCA to identify the correlation between paired datasets. iCluster is a method able to generate a single integrated cluster assignment based on a simultaneous inference analysis from multiple data types . In BC, iCluster has been used to align concordant DNA CNA and gene GE changes, showing encouraging results .
Table 9 reports some software for the CNA and GE analysis and their method of integration type.
Integrating genetics and epigenetics in BC
Sarkar et al. suggested that the epigenetic changes act as the initiating signal in the development of cancer progenitor cells and a combination of all genetic changes which are differentially expressed in the various cancer subtypes, could act on the cell vulnerable to epigenetic alterations .
Epigenetic mechanisms are tightly linked to one another and make the overall gene regulation system. The miR-29 family, for example, including miR-29a, miR-29b, and miR-29c, is a miRNA that collaborates with other epigenetic mechanisms. The expression of miR-29b is regulated by both histone modification  and DNA methylation . miR-7/miR-218 can regulate DNA methylation and histone modification status by decreasing homeobox B3 (HOXB3) expression .
However, while classical epigenetic mechanisms, such as histone modification and DNA methylation, regulate expression at the transcriptional level, miRNAs act at the posttranscriptional level.
Elucidating the basic mechanisms of post-transcriptional regulation of GE is essential to gain a full understanding of how GE is regulated at different levels, of the interplay between these mechanisms, and of the extensive contribution of post-transcriptional dysfunction in cancer.
An impressive number of papers have been published on miRNAs increasing the number of scientific challenges, and we focused on the studies and methods applied to the combination miRNAs-mRNA, CNAs-miRNAs, and GE-genetic alterations-miRNAs.
Integrated analysis of mRNA and miRNA in BC
The miRNA profile is more accurately associated with cell differentiation and cancer progression when compared with GE expression profile.
The aberrant expression of miRNAs in cancer can lead to the altered expression of target mRNAs. miRNAs can also modulate multiple genes regulating entire networks. The interaction of a miRNA with its target mRNAs can lead to the repression or incentive of GE.
Combination of miRNA and mRNA has still to be deeply explored in diagnostic and prognostic studies. Cascione et al.  proposed a large-scale analysis of miRNA and cancer-focused mRNA expression in normal, triple negative tumour, and associated metastatic tissues in BC. Two miRNA signatures were identified, predictive of overall survival (P = 0.05) and distant-disease free survival (P = 0.009), respectively. Volinia et al.  found 30 mRNAs and 7 miRNAs associated with overall survival, across different clinical and molecular subclasses of BCs. In addition, expression profiles from 8 BC datasets, different from those used for the miRNA extraction, were used for validation. Buffa et al.  matched mRNA and miRNA global expression profiling, and four miRNAs were found independently associated with DRFS in ER-positive BC (3 novel and 1 known miRNA- miR-128a) and six miRNAs in ER-negative BC (5 novel and 1 known miRNA; miR-210). Van der Auwera et al.  identified a set of 13 miRNAs whose expression differed between inflammatory BC (IBC) and non-IBC. Enerly et al.  demonstrated, from the joint analysis of miRNA and mRNA data, a central role for miRNAs in regulating particular pathways. Hannafon et al.  identified putative miRNAs by mRNA functional interactions in ductal carcinoma in situ: the three miRNAs miR-125b, miR-182 and miR-183, and six of their putative target genes, MEMO1, NRIP1, CBX7, DOK4, NMT2, and EGR1.
Luo et al.  performed an integrated analysis of miRNAs and mRNA expression profiles in 12 BC cell lines, identifying 35 functional target genes of three significantly down-regulated miRNAs in invasive cell lines (miR-200c, miR-205, and miR-375).
Several studies demonstrated the greater accuracy of miRNA expression levels compared with those of gene signatures. miRNA expression levels should directly represent the functional activity of the genes, while genes have to be translated to proteins to show their biological effects .
For a more detail review on the role of miRNAs and mRNA, see .
Table 10 reports a synthesis of the considered miRNAs biomarkers in BC as obtained by the integration of miRNA and mRNA.
Each miRNA can potentially regulate the expression of hundreds of genes, and a single gene can be targeted by multiple miRNAs .
Specific miRNAs has been identified as regulator of metastatic progression through miRNA regulatory networks. Yan et al.  found miR-21 as the most significantly up-regulated miRNA in BC when compared with normal adjacent tumour tissues (NAT). Its target prediction revealed the putative target genes by creating a small miRNA regulatory networks.
Figure 6 shows mTOR and STAT3 signalling acting on miR-21 up-regulation in cancer  and miR-21 promoting cancer cell invasion and metastasis through suppression of BCL-2, PTEN, PDCD4,TPM1, maspin . The introduction of anti-miR-21 to MCF-7 BC cells and in mouse model resulted in decreased cell growth (via increased apoptosis) and in reduced cell proliferation .
miR-10b was found highly expressed in BC metastatic cancer cells. In vivo studies demonstrated that miR-10b promotes cell migration and invasion [320, 321] and initiates tumour metastasis [320, 321]. miR-10b is induced by the transcription factor Twist. In turn, miR-10b inhibits HOXD10 and, through a cascade of cellular alterations, inhibits the expression of the prometastatic gene RHOC [320, 321].
let-7 has been found poorly expressed or deleted in many cancers. Known oncogenic targets of let-7 are H-RAS, HMGA2, and BACH1. These genes result down-regulated by let-7 over-expression . HMGA2, and BACH1 promote the transcription of pro-invasive genes, suppress cell invasion and metastasis to the bone . let-7 is regulated by LIN-28, MEK signalling, and RKIP .
miR-200 family is important in maintaining the tumour epithelial phenotype and in inhibiting the epithelial-to-mesenchymal transition (EMT). miR-200 family was found to inhibit cell migration by acting on the transcription factors ZEB1 and ZEB2, which suppress E-cadherin . Furthermore, miR-200 was found silencing Sec-23a and promoting metastases by inhibiting TINAGL1 and IGFBP4 .
Table 11 reports a synthesis of the considered mRNA-miRNA networks.
The most used experimental technique for determining miRNA targets is the transfection of mimic miRNAs or of miRNA inhibitors . The consequences of the modulation of miRNAs on the expression levels are measured by using different tools, including RT-PCR or microarrays. The most important disadvantage of these techniques is that they are not able to discriminate between indirect and direct interactions . Labelled miRNA pull-down (LAMP) assay system or luciferase report assays add reporters or labels to miRNAs on the 3'-UTR of transcripts of interest, allowing the identification and the analysis of direct interaction regions among miRNA and its target gene . The disadvantage of reporter assays is that they are laborious, sensitive upon the region chosen for cloning, and that they require hard and complex work for trasfection .
There are different approaches to examine both miRNA and mRNA expression profiles. In this paragraph we examine miRNA and mRNA regulatory pairs together [183, 185, 187, 188]. Several studies showed that the miRNA-mRNA interactions varies with the development of disease [331, 332].
The integrative methods employ a three-step procedure: 1) Identification of DE miRNAs and mRNAs in the biological condition of interest. It can be done as reported in section 2 c); 2) Selection of putative miRNA-mRNA pairs (for instance, a prediction algorithm can be used to obtain the DE miRNA from DE mRNA. It can be done as reported in section 3 c); and 3) Identification of statistically significant miRNA-mRNA pairs. This last step needs the selection of an appropriate association measure, and the determination of its significance. The common assumption is based on the idea that regulatory relationship between any miRNA and its target mRNAs is an inverse correlation .
The mathematic tools consider simple correlation analyses (Pearson, Spearman) [336, 337], mutual information , linear regression [338–340], regularized least squares [341, 342] and bayesian inference [343, 344]. These methods give a score for each interaction mRNA-miRNA.
van Iterson et al.  used the global test  to associate each miRNA with the expression levels of a set of predicted mRNA targets. They suggest global tests to be better suited for integrated analysis of miRNA and mRNA expression data, compared with either Pearson correlation or lasso-based approaches.
Pearson Correlation is a measure of linear-dependency, widely used to show miRNA-mRNA showing a statistically significant correlation . There are several web-tools that employ Pearson correlation for miRNA-mRNA target research (e.g. [349–353]).
Non-parametric (Spearman) correlation coefficient can be used as alternative measure of correlation. Usually it is chosen in case of outliers or with small number of measures. Contrary to Spearman correlations, Pearson coefficients require that both variables derive from a bi-variate normal distribution .
Mutual information is analogous to the Pearson Correlation but it is sensitive not just to linear dependencies, and can define whether two given variables are independent .
Multiple linear regression can evaluate the interactions between a set of miRNAs and a target mRNA, contrary to correlation measures which focuses on particular pairs interaction.
R-squared statistics is used for measuring the goodness of the fit of the data . When the number of samples with GE profiles is smaller than the number of covariates (e.g. miRNA), partial least squares can be applied . This model gives those miRNAs explaining the maximum variance in GE profiles by ensuring a good fit of the model.
Lasso-based approaches are used to deal with undetermined linear system .
Bayesian inference use a priori information to estimate parameters and predict values in a probability framework. Several studies use this method for scoring putative miRNA-mRNA targets based on miRNA and mRNA expression data [352–354].
Table 12 reports a synthesis of methods considered for the integration analysis mRNA-miRNA.
In the context of a network, miRNAs are able to regulate distinct biological cell processes like apoptosis, proliferation or receptor driven pathways, thus suggesting their possible use also as therapeutic targets or tools . The most important advantage, with respect to other approaches targeting single genes, is their ability to target multiple molecules.
There are two main approaches to target miRNA expression in cancer. Direct approaches involve the use of oligonucleotides or virus‐based vectors to either block the expression of an oncogenic miRNA or to reintroduce a TS miRNA lost in cancer. Indirect approaches involve the use of drugs to modulate miRNA expression by targeting their transcription and their processing.
We think that the miRNAs described in the following sections could be interesting for the development of possible therapies in BC.
Ma et al.  found miR-10b up-regulated in BC and explored a possible therapeutic application in an animal model of BC-bearing mice. The silencing of miR-10b with antagomiRs reduces miR-10b levels and increases miR-10b target, HOXD10. The therapy decreases metastases and was well tolerated by mice.
Multiple studies have also shown a significant association between miR-206 and ER in BC (e.g. ). In mouse models, the overexpression of miR-206 was found significantly decreasing metastatic activity for 2 BC cell lines: BOM1 (highly metastatic to bone) and LM2 (highly metastatic to lung) .
miR-125 was found to be significantly down regulated in BC patients . Experimentally, over-expression of miR-125 reduces ERBB2 and ERBB3 cell motility, and also reduces invasiveness of other numerous cancers [357, 358]
miR-34 is down regulated in BC cell lines and tissues, compared with normal cell lines and adjacent non-tumor tissues . Expression of miR-34 was found correlated with p53 status. In fact, silencing of p53 in human tumour cell lines decreases in miR-34 level . Moreover, as reported by Weidhaas et al.  miR-34 levels change levels significantly after irradiation. A potential use for miR-34 as radiosensitizing agent could be envisaged.
miR-155 is also linked to key cancer pathways as the gene is up-regulated by mutant p53 in BC, thus facilitating tumour cell invasion . miR-155 has also attracted considerable interest as a putative therapeutic target .
Table 13 reports a synthesis of the considered miRNAs, their potential target and function.
Integrated analysis of CNA and miRNA in BC
Many miRNAs are frequently located at fragile sites of the genome, which are usually either amplified or deleted in human cancer . The aberrant miRNA expression in BC, in part, is due to these genomic alternations.
Zhang and colleagues studied 283 known human miRNAs in BC and showed that 72.8 % of miRNAs are located in regions that reveal CNAs . In a recent study, miRNAs were shown to be up-regulated in gain regions compared to copy-neutral regions in BC, although the effect on miRNA expression was not incisive . Iorio et al.  compared BC CGH data with independent miRNA expression by miRNA microarrays, and demonstrated that 81.8 % of miRNAs increased expression level and showed high DNA CN, and that 60 % of miRNAs exhibit decreased expression level with loss of DNA CN.
Several miRNAs have been associated with cancers due to CNA, suggesting that miRNAs can act either as oncomiRs or oncosuppressor miRNA . Figure 7 shows amplification of chromosomal regions of miRNAs encoding oncomiRs and leading to their up-regulation. OncomiRs can act silencing TSG thus making possible the development of cancer.
The first miRNA found to act as a mammalian oncogene is polycistron miR-17-92, also known as OncomiR-1 because it was the first identified oncomiR . It is located in chromosome 13 and has been found amplified in human BC . It acts as an anti-apoptotic miR cluster by targeting intrinsic apoptotic protein Bim in B-cell lymphoma subtypes .
Other oncomiRs have been described since the first discovery. miR-21 is located in 3'UTR of VMP1 (vacuole membrane protein 1) gene at chromosome 17q23.2, a region amplified in BC and also in neuroblastomas, colon and lung cancers . miR-151a-5p is located on 8q24.3, a genomic site frequently associated with gain in BC . High expression of miR-151a-5p has been associated with gain, and functional experiments showed that over-expression induce cell proliferation and also increase the levels of p-AKT .
As for oncomiRs, also several miRNAs with oncosuppressor functions have been described. Figure 8 shows deletion of chromosome region of oncosuppressor miRNAs leading to their down-regulation. Down-regulation of oncosuppressor miRNAs results in up-expression of target oncogenes.
Chromosome 11 is frequently altered in BC and mirR-125b, that is located at 11q23-24, results one of the most frequently deleted regions . In a study of Muller et al. , mir-320 has been found to be located in regions with DNA CN loss in BC. The predicted target of miR-320 is MECP2 which is up-regulated in BC and serves as an oncogene promoting cell proliferation. Genetic deletion could contribute to miR-100 down-regulation  inducing epithelial-mesenchymal transition.
In several cancer types, including BC, genomic deletion or loss of heterozygosis of the region of the miR-34a have been described . miR-34a is highly expressed in normal tissues. Its expression level is under the control of the TS gene product p53 and it acts as a TS inducing cell cycle arrest in G1-phase, senescence and apoptosis .
Wang et al.  showed that CN deletion is an important mechanism leading to the down-regulation of expression of specific let-7 family members in BC. Also miR-33 expression was found to be strongly associated with the genomic alteration . Furthermore, the expression of the cluster miR-145/miR-143 family, miRNA located on a region involved in several types of translocations and deletions, has been found reduced or absent in various types of cancers, including BC [152, 368].
Table 14 shows the principal oncomiRs and oncogenes with their alterations considered in this section.
miRNAs that are silenced or amplified from CNA can have a cascade effect on the expression of different genes regulating entire pathways.
In the following paragraph, we give examples of important miRNAs that are altered in BC and of the consequences of their downregulation in the functional pathway.
Figure 9a shows miR-335 that suppresses BC metastasis by targeting SOX4 and Tenascin-C which promote cancer cell migration, invasion and ultimately metastasis [326–328]. miR-335 is silenced through CN deletions .
mir-320 is found to be located in regions with CN loss in BC. The predicted target of miR-320 is methyl CpG-binding protein 2 (MECP2), which is up-regulated in BC and is an oncogene promoting cell proliferation .
In a study of Volinia et al.  miR-21 was found as the only miRNA up-regulated in all six types of solid cancers (BC, colon, lung, prostate, stomach carcinomas and pancreas exocrine tumours). Figure 9b shows miR-21 network: it modulates gemcitabine-induced apoptosis by PTEN-dependent activation of PI 3-kinase and by activation of AKT/mTOR signalling . Inhibition of this miRNA should result in cell death .
Several studies showed that miRNA levels are influenced by CNAs.
No-experimental methods are usually used for their integration. Individual studies from miRNA and CNA are combined with statistically and/or computational analysis.
de rinaldis et al.  analysed association between miRNA expression and CNAs in a large triple-negative BC data set. This association was evaluated using Spearman correlation. In addition, for each miRNA-encoding DNA locus identified as altered in any of the samples, a separate non-parametric Wilcoxon rank sum test was applied to measure differences in expression between samples with deletions and amplifications, compared to samples with no CNAs. 64 miRNAs were found with statistically significant miRNA-CNA correlation, showing an overall influence of genetic alterations (amplifications and deletions) on the expression of the miRNAs.
Aure et al.  investigated individual and combined effects of CN and methylation on miRNA expression in BC. They identified 70 miRNAs whose expression was associated with CNAs or methylation, or with both conditions. 24 miRNAs were associated mainly with CNAs, 22 miRNAs with methylation aberrations and 24 miRNAs with a combination of CN and methylation aberrations. In order to identify miRNAs associated with hypomethylation or amplification, each miRNA in each patient was allocated to one of the two groups ‘altered’ or ‘non-altered’ based on CNA and DNA methylation. A Wilcoxon rank-sum test was used for each miRNA to underlie whether the miRNA expression was significantly different in the two groups.
Srivastava et al.  showed that H2AX was negatively correlated with miR-24-2 and not in accordance with the CNA status, both in cell lines and in sporadic BC tissues. The authors tried to explain the possible mechanisms of such non concordant relationship between expression and number of gene copies based on specific miR regulation of expression. They discussed a role of miR-24-2 in guiding H2AFX GE in the background of the differential status of CNA.
Combination of gene expression, genetic alterations and miRNAs in BC
Fearon and Vogelstein  proposed that accumulation of genetic alterations could determine a malignant phenotype and accompany cancer progression. However, this theory does not explain the great heterogeneity of observed genetic alterations, even within homogeneous histological groups .
Normal cells evolve progressively to a neoplastic state, based on a multistep process to acquire the traits that enable them to become tumorigenic and ultimately malignant. Tumors are not only masses of proliferating cancer cells, but complex tissues composed of multiple distinct molecular types that participate in an interaction with one another [383, 384].
The transitions in the malignant cancer progression are dynamic and reversible steps between multiple phenotypic states (e.g. epithelial and mesenchymal phenotype) . These reversible transitions are based on complex epigenetic regulatory mechanisms (e.g. the induction of changes in the modifications of chromatin-associated histones) during epithelial-mesenchymal transitions [385, 386].
Sarkar et al.  reported a review based on the role of epigenetic regulation in the steps from normal cell to cancer progenitor cells that, after growing, undergo an epithelial-mesenchymal transition. Epigenetic drugs could potentiate traditional therapeutics by inhibiting both the formation and growth of cancer progenitor cells .
We argue that tumour heterogeneity is due not only to a simple accumulation of genetic alterations but can be the cause of the combined effect of genetic and epigenetic alterations. Furthermore, Alfred Knudson  hypothesized that hereditary retinoblastoma involves two mutations, the first one in the germ line. Thus, non-hereditary retinoblastoma should be due to two somatic mutations, an hypothesis known as Knudson “two-hit” hypothesis. The two-hit hypothesis proposes that loss of a single functional allele, which may potentially results in expression of a truncated or mutated product, is insufficient to involve cellular functions.
Several studies support the validity of the "two hit theory" in BC. Meric-Bernstam et al.  applied this hypothesis in BC, and suggested that the second hit does not need to be a point of mutation or somatic loss, but it may be the epigenetic silencing of a gene.
Konishi et al.  showed that cell lines carrying one mutant and one normal copy of BRCA1 have a normal cell phenotype, and they are normal until the second allele is lost through somatic mutation or epigenetic silencing.
Genetic and epigenetic events are two complementary mechanisms that are involved in carcinogenesis. It is not clear at all how these mechanisms influence GE during tumorigenesis.
In BC, integration analysis of GE, genomic changes and miRNA expression was adopted in a limited number of studies (e.g. [397–399]). Eo et al.  proposed a pathway-based classification of BC which integrates data on DE genes, CNA and miRNA. Pathway information was incorporated in a condition-specific manner. A 215-gene signature was found from 327 tumours. By using an independent data set, this gene signature was validated.
Cancer Genome Atlas Network  analysed BC by genomic DNA CN arrays, DNA methylation, exome sequencing, messenger RNA arrays, microRNA sequencing and reverse-phase protein arrays. They found biomarkers for gene expression subtypes and the presence of four main BC classes.
Kristensen et al.  used an integrated approach to identify and classify BC according to the most deregulated pathways that provide the best predictive value with respect to prognosis, and identified key molecular and stromal signatures.
In a combined analysis of miRNA and mRNA expression data, Blenkiron et al.  found a number of miRNAs DE among molecular tumour subtypes. Furthermore, they found that changes in miRNA expression correlate with genomic loss or gain.
Cava et al.  assessed the potential of a new triple approach by integrating mRNA expression profile, CNAs, and miRNA expression levels to select a limited number of genomic BC biomarkers and to obtain a more accurate classification of BC grade.
CNAs have been demonstrated to be able also to identify genes DE between drug-sensitive and -resistant BC cells when integrated to GE and microRNA expression profiles.
Yamamoto et al.  focused on miRNAs and genes located on the genome-amplified and -deleted regions. These genes showed also an altered expression in GE profiles. The authors analysed MCF7 and a parental BC cell line drug-resistance MCF7-ADR. miR-505 was identified as a tumour suppressor, whose genomic region was found to be deleted in doxorubicin-resistant cells. Furthermore, miR-505 seems to be regulated by its predicted target Akt3 (an anti-apoptotic gene), by mRNA profiling coupled with downstream validation studies.
Despite promising initial results about the possible clinical implications of GE profiling, a more recent source of concern has been that gene signatures derived from the various studies show little overlap and poor reproducibility. This can be explained, from one side, by the complexity of the human genome which provides that different genes can be indices of the same message with identical outcomes. From the other side, one explanation can be the use of different types of arrays (of different sample quality) and the different parameters considered for the data analysis. However, GE analysis measures mRNA expression, which, by the central dogma of molecular biology, results from the transcription of DNA. Specifically, GE analysis give information on DE genes among different conditions, but do not consider primary alterations of DNA from secondary effects of disease, such as, in the case of cancer, proliferation and differentiation state. On the other hand, studies of DNA CNA allow important indices to be derived as drivers of cancer. Therefore, integrating DNA and RNA data has been proposed to clarify genetic some regulatory relationships in cancer cells.
Since 2001, a new term "microRNA" was introduced into the scientific literature, challenging the central dogma of molecular biology. miRNAs are segments of RNA that are transcribed from DNA in a way similar to mRNA but they are not translated into proteins. In short, instead of producing a protein, miRNA can block mRNA directly. Evidences demonstrated that their deregulation is associated to several steps of cancer initiation and progression. However, we think that the association of miRNAs and their mRNA targets is a more favourable approach to study cell differentiation and cancer progression when compared with GE expression or miRNA profile alone. It is therefore of great concern for researchers to investigate how miRNA expression is linked to known BC markers. Several advantages can be envisaged by miRNA analysis: i) miRNAs are certainly more stable due to their small size when compared to long mRNAs , ii) miRNA expression levels can characterize the functional activity of the target gene while genes have to be translated to proteins to be biologically functional iii) miRNA-based therapeutics have the ability to target multiple genes.
Misregulation of genes with consequence disruption of the gene function is often induced by epigenetic and genetics events. The epigenetic silencing of one allele may act in concert with an inactivating genetic alteration in the opposite allele, thus resulting in a total allelic loss of the gene [7, 8]. From this viewpoint a gene subjected to a different possible alterations (such as CNAs and target of miRNAs) and that presents DE levels between two conditions is a "weak" point of DNA and could be a key element for cancer development. In our opinion each cancer should have a signature with the description of a specific set of alterations. Based on these observations, targeting specifically and simultaneously multiple pathways subjected to different alterations may confer a greater therapeutic efficacy.
We argue that useful information has been revealed by analysing GE profiles alone, CNA data alone but or miRNAs, however, in order to have complimentary information in gene characterization, an integrative analysis of CNA and GE data and miRNA is necessary.
However, integrative analyses have some limitations: the most fundamental challenge is dimensionally, considering that more levels in the analysis increase the computational time and the dimension of unknown parameters . In addition, at every step, there are problems of compatibility of the data, such as normalization to the same scale, batch effects, and use of different platforms.
Large-scale integration is possible only for few projects worldwide, given the high cost for all analyses to be carried out simultaneously and on the entire data set.
In referring to current studies of genetic changes associated with BC, we focused in particular on the processes controlled by CNA. However, DNA changes include other genomic rearrangements, such as somatic point mutations.
The analysis of the genomes of 100 tumours revealed more than 7400 somatic point mutations in 21416 protein-coding genes . These mutations affect many of the well-established cancer related genes, such as BRCA1, RB1, TP53, PTEN, AKT1, CDH1, GATA3, PI3KCA. These genes control apoptosis, proliferation and cell cycle, and transcription. Other somatic mutations affect genes involved in signal transduction (APC, KRAS, MAPK2K4, SMAD4, CASP8, CDKN1B…). Somatic mutation in three main genes (TP53, PI3KCA, and GATA3) shows more than 10 % incidence across all BC . One of the most commonly mutated TSG in BC is P53 . It is localized to chromosome 17p13 and its inactivation is important also in other cancer diseases. Several studies have investigated the predictive power of P53 for response to treatments and outcome of BC patients [399–401]. Bertheau et al.  reported that P53 base-pair substitutions are highly linked to specific BC molecular subtypes, being found in 26 % of luminal tumours (17 % of luminal A, 41 % of luminal B), in 50 % of HER2 amplified tumours, and in 88 % of basal-like carcinomas. The type of mutations changes according to the tumour subtype. Basal-like tumours present higher frequency of deletions. Furthermore, the authors found that non inflammatory locally advanced BC with mutated P53 has a higher rate of response to dose-dense doxorubicin–cyclophosphamide chemotherapy than TP53-WT tumours. As recently reported , P53 is at the centre of the hallmarks of cancer, supporting genomic stability, exerting anti-angiogenic effects, controlling tumour inflammation and immune response, and repressing metastases. In BC, mutations in BRCA1 and BRCA2 result in protein truncations as consequence of small insertions, deletions or nonsense mutations. Although BRCA1 and BCRA2 mutations are hereditary, these genes would also be involved in the development of sporadic BC. Compared with normal breast epithelium, many BCs have shown low levels of the BRCA1 mRNA [403, 404], while BRCA2 has been found the target of frequent loss of heterozygosity (LOH) in BC [405, 406].
Other omics data could be further integrated for a more inclusive analysis. Considering that proteins translate effects of CNAs into the biological functions of the cell, further studies could integrate protein-protein interactions networks with gene-gene co-expression networks. For example, by dissecting the protein-protein interaction network into disjoint sub networks, van den Akkerb et al.  found sub-population of genes by using pair wise GE correlation measures. The obtained genes were consistently found across different studies.
Also the DNA methylation could be integrated in a pathway analysis and could be combined with other biological data. Andrews et al.  integrated results from CNAs, GE profiling and methylation to identify differentially regulated pathways between a highly metastatic BC cell line and low metastatic parental cell line. Validation experiments confirmed that hypermethylated genes correlated with decreased expression in the metastatic, compared to the parental cell line.
Results generated from whole-genome analyses have been submitted in The Cancer Genome Atlas (TCGA) database, which includes CNAs, DNA methylation and GE profiles [409, 410]. These data might be used for integrative analyses of results generated from a single technology platform .
Integrating genetics and epigenetics in BC may offer a powerful approach for the identification of biomarkers with diagnostic, prognostic and therapeutic potential. The experimental and computational methods presented in this review can be used to guide researchers for these integration studies.
Copy number alteration
Human epidermal growth factor receptor-2
Tumor suppressor genes
Loss of heterozygosity
Fluorescence in situ hybridization technique
Single nucleotide polymorphism
Next generation sequencing
Hidden Markov model
Copy number analyser for GeneChip arrays
Robust multichip average
- EGFR :
Epidermal growth factor receptor
Real-time quantitative reverse transcription PCR
Polyak K. Heterogeneity in breast cancer. J Clin Invest. 2011;10:3786–8.
Viale G. The current state of breast cancer classification. Ann Oncol. 2012;23 Suppl 10:207–10.
Hsiao YH, Chou MC, Fowler C, Mason JT, Man YG. Breast cancer heterogeneity: mechanisms, proofs, and implications. J Cancer. 2010;1(1):6–13.
Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, et al. Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst. 2006;98(4):262–72.
Cava C, Bertoli G, Ripamonti M, Mauri G, Zoppis I, Della Rosa PA, et al. Integration of mRNA Expression Profile, Copy Number Alterations, and microRNA Expression Levels in Breast Cancer to Improve Grade Definition. PLoS One. 2014;9(5):e97681.
Ivshina AV, George J, Senko O, Mow B, Putti TC, et al. Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006;66(21):10292–301.
Pinto R, De Summa S, Petriella D, Tudoran O, Danza K, Tommasi S. The value of new high-throughput technologies for diagnosis and prognosis in solid tumors. Cancer Biomark. 2014;14(2–3):103–17.
Esteller M, Fraga MF, Guo M, Garcia-Foncillas J, Hedenfalk I, Godwin AK, et al. DNA methylation patterns in hereditary human cancers mimic sporadic tumorigenesis. Hum Mol Genet. 2001;10(26):3001–7.
Birgisdottir V, Stefansson OA, Bodvarsdottir SK, Hilmarsdottir H, Jonasson JG, Eyfjord JE. Epigenetic silencing and deletion of the BRCA1 gene in sporadic breast cancer. Breast Cancer Res. 2006;4:R38.
Li Z, Chen B, Wu Y, Jin F, Xia Y, Liu X. Genetic and epigenetic silencing of the beclin 1 gene in sporadic breast tumors. BMC Cancer. 2010;10:98.
Yang Q, Nakamura M, Nakamura Y, Yoshimura G, Suzuma T, Umemura T, et al. Two-hit inactivation of FHIT by loss of heterozygosity and hypermethylation in breast cancer. Clin Cancer Res. 2002;9:2890–3.
Feinberg AP, Ohlsson R, Henikoff S. The epigenetic progenitor origin of human cancer. Nat Rev Genet. 2006;7(1):21–33.
Gonzalez-Angulo AM, Hennessy BT, Mills GB. Future of personalized medicine in oncology: a systems biology approach. J Clin Oncol. 2010;28(16):2777–83.
Beroukhim R, Mermel CH, Porter D, Wei G, Raychaudhuri S, Donovan J, et al. The landscape of somatic copy-number alteration across human cancers. Nature. 2010;463(7283):899–905.
Poplawski AB, Jankowski M, Erickson SW, Diaz de Stahl T, Partridge EC, Crasto C, et al. Frequent genetic differences between matched primary and metastatic breast cancer provide an approach to identification of biomarkers for disease progression. Eur J Human Genet. 2010;18:560–8.
Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL. Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science. 1987;235:177–82.
Staaf J, Jönsson G, Ringnér M, Vallon-Christersson J, Grabau D, Arason A, et al. Research article High-resolution genomic and expression analyses of copy number alterations in HER2-amplified breast cancer. Breast Cancer Res. 2010;12:R25.
Rizzolo P, Silvestri V, Falchetti M, Ottini L. Inherited and acquired alterations in development of breast cancer. Appl Clin Genet. 2011;4:145.
Faivre EJ, Lange CA. Progesterone receptors upregulate Wnt-1 to induceepidermal growth factor receptor transactivation and c-Src-dependent sustained activation of Erk1/2 mitogen-activated protein kinase in breast cancer cells. Mol Cell Biol. 2007;27(2):466–80.
Tsutsui S, Ohno S, Murakami S, Hachitanda Y, Oda S. Prognostic value of epidermal growth factor receptor (EGFR) and its relationship to the estrogen receptor status in 1029 patients with breast cancer. Breast Cancer Res Treat. 2002;71:67–75.
Knoop AS, Knudsen H, Balslev E, et al. Retrospective analysis of topoisomerase IIa amplifications and deletions as predictive markers in primary breast cancer patients randomly assigned to cyclophosphamide, methotrexate, and fluorouracil or cyclophosphamide, epirubicin, and fluorouracil: Danish Breast Cancer Cooperative Group. J Clin Oncol. 2005;23:7483–90.
O'Malley FP, Chia S, Tu D, et al. Topoisomerase II α and responsiveness of breast cancer to adjuvant chemotherapy. J Natl Cancer Inst. 2009;101:644–50.
Tanner M, Isola J, Wiklund T, et al. Topoisomerase IIα gene amplification predicts favorable treatment response to tailored and dose-escalated anthracycline-based adjuvant chemotherapy in HER-2/neu–amplified breast cancer: Scandinavian Breast Group Trial 9401. J Clin Oncol. 2006;24:2428–36.
Gonzalez‐Angulo AM, Chen H, Karuturi MS, Chavez‐MacGregor M, Tsavachidis S, Meric‐Bernstam F, et al. Frequency of mesenchymal‐epithelial transition factor gene (MET) and the catalytic subunit of phosphoinositide‐3‐kinase (PIK3CA) copy number elevation and correlation with outcome in patients with early stage breast cancer. Cancer. 2013;119(1):7–15.
Xu J, Chen Y, Olopade OI. MYC and breast cancer. Gene Canc. 2010;1(6):629–40.
Aulmann S, Bentz M, Sinn HP. C-myc oncogene amplification in ductal carcinoma in situ of the breast. Breast Cancer Res Treat. 2002;74:25–31.
Robanus-Maandag EC, Bosch CA, Kristel PM, et al. Association of C-MYC amplification with progression from the in situ to the invasive stage in C-MYC-amplified breast carcinomas. J Pathol. 2003;201:75–82.
Aulmann S, Adler N, Rom J, Helmchen B, Schirmacher P, Sinn HP. c-myc amplifications in primary breast carcinomas and their local recurrences. J Clin Pathol. 2006;59:424–8.
Corzo C, Corominas JM, Tusquets I, et al. The MYC oncogene in breast cancer progression: from benign epithelium to invasive carcinoma. Cancer Genet Cytogenet. 2006;165:151–6.
Lundgren K, Brown M, Pineda S, Cuzick J, Salter J, Zabaglo L, et al. Effects of cyclin D1 gene amplification and protein expression on time to recurrence in postmenopausal breast cancer patients treated with anastrozole or tamoxifen: a TransATAC study. Breast Cancer Res. 2012;14(2):R57.
Sherr CJ, Roberts JM. CDK inhibitors: positive and negative regulators of G1-phase progression. Genes Dev. 1999;13:1501–12.
Holst F, Stahl PR, Ruiz C, et al. Estrogen receptor alpha (ESR1) gene amplification is frequent in breast cancer. Nat Genet. 2007;39:655–60.
Desouki MM, Liao S, Huang H, Conroy J, Nowak NJ, Shepherd L, et al. Identification of metastasis-associated breast cancer genes using a high-resolution whole genome profiling approach. J Cancer Res Clin Oncol. 2011;137:795–809.
Rodriguez C, Hughes-Davies L, Vallès H, et al. Amplification of the BRCA2 pathway gene EMSY in sporadic breast cancer is related to negative outcome. Clin Cancer Res. 2004;10:5785–91.
Wang C, Iakovlev VV, Wong V, Leung S, Warren K, Iakovleva G, et al. Genomic alterations in primary breast cancers compared with their sentinel and more distal lymph node metastases: an aCGH study. Gene Chromosome Canc. 2009;48:1091–101.
Trapé AP, Gonzalez-Angulo AM. Breast cancer and metastasis: on the way toward individualized therapy. Cancer Genomics-Proteomics. 2012;9(5):297–310.
Imataka G, Arisaka O. Chromosome analysis using spectral karyotyping (SKY). Cell Biochem Biophys. 2012;62(1):13–7.
Salman M, Jhanwar SC, Ostrer H. Will the new cytogenetics replace the old cytogenetics? Clin Genet. 2004;66:265–75.
Alkan C, Coe BP, Eichler EE. Genome structural variation discovery and genotyping. Nat Rev Genet. 2011;12(5):363–76.
Park H et al. Discovery of common Asian copy number variants using integrated high-resolution array CGH and massively parallel DNA sequencing. Nat Genet. 2010;42:400–5.
Li W, Olivier M. Current analysis platforms and methods for detecting copy number variation. Physiol Genomics. 2013;45(1):1–16.
Clevert DA, Mitterecker A, Mayr A, Klambauer G, Tuefferd M, De Bondt A, et al. cn. FARMS: a latent variable model to detect copy number variations in microarray data with a low false discovery rate. Nucleic Acids Res. 2011;39(12):e79–9.
Huang J, Wei W, Zhang J, Liu G, Bignell GR, Stratton MR, et al. Whole genome DNA copy number changes identified by high density oligonucleotide arrays. Hum Genomics. 2004;1:287–99.
Zhao X, Li C, Paez JG, Chin K, Janne PA, Chen TH, et al. An integrated view of copy number and allelic alterations in the cancer genome using single nucleotide polymorphism arrays. Cancer Res. 2004;64(9):3060–71.
Nannya Y, Sanada M, Nakazaki K, Hosoya N, Wang L, Hangaishi A, et al. A robust algorithm for copy number detection using high-density oligonucleotide single nucleotide polymorphism genotyping arrays. Cancer Res. 2005;65(14):6071–9.
Korn JM, Kuruvilla FG, McCarroll SA, Wysoker A, Nemesh J, Cawley S, et al. Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs. Nat Genet. 2008;40:1253–60.
Bengtsson H, Irizarry R, Carvalho B, Speed TP. Estimation and assessment of raw copy numbers at the single locus level. Bioinformatics. 2008;24:759–67.
Gai X, Perin JC, Murphy K, O'Hara R, D'arcy M, Wenocur A, et al. CNV Workshop: an integrated platform for high-throughput copy number variation discovery and clinical diagnostics. BMC Bioinformatics. 2010;11(1):74.
Baross A, Delaney AD, Li HI, Nayar T, Flibotte S, Qian H, et al. Assessment of algorithms for high throughput detection of genomic copy number variation in oligonucleotide microarray data. BMC Bioinformatics. 2007;8:368.
Eckel-Passow JE, Atkinson EJ, Maharjan S, Kardia SL, de Andrade M. Software comparison for evaluating genomic copy number variation for Affymetrix 6.0 SNP array platform. BMC Bioinformatics. 2011;12:220.
Zhang D, Qian Y, Akula N, Alliey-Rodriguez N, Tang J, Gershon ES, et al. Accuracy of CNV detection from GWAS data. PLoS One. 2011;6:e14511.
Gordon DJ, Resio B, Pellman D. Causes and consequences of aneuploidy in cancer. Nat Rev Genet. 2012;13:189–203.
Carling D. The AMP-activated protein kinase cascade—a unifying system for energy control. Trends Biochem Sci. 2004;29:18–24.
Paez JG, Ja¨nne PA, Lee JC, Tracy S, Greulich H, Gabriel S, et al. EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004;304:1497–500.
Shih J, Bashir B, Gustafson KS, Andrake M, Dunbrack RL, Goldstein LJ, et al. Cancer Signature Investigation: ERBB2 (HER2)-Activating Mutation and Amplification-Positive Breast Carcinoma Mimicking Lung Primary. J Natl Compr Canc Netw. 2015;13(8):947–52.
Baselga J, Norton L, Albanell J, Kim YM, Mendelsohn J. Recombinant humanized anti-HER2 antibody (Herceptin) enhances the antitumor activity of paclitaxel and doxorubicin against HER2/neu overexpressing human breast cancer xenografts. Cancer Res. 1998;58:2825–31.
Baselga J, Cortés J, Kim SB, Im SA, Hegg R, Im YH, et al. Pertuzumab plus trastuzumab plus docetaxel for metastatic breast cancer. N Engl J Med. 2012;366(2):109–19.
Scheuer W, Friess T, Burtscher H, Bossenmaier B, Endl J, Hasmann M. Strongly enhanced antitumor activity of trastuzumab and pertuzumab combination treatment on HER2-positive human xenograft tumor models. Cancer Res. 2009;69:9330–6.
Bird AP. CpG-rich islands and the function of DNA methylation. Nature. 1986;321:209–13.
Ehrlich M, Gama-Sosa MA, Huang LH, Midgett RM, Kuo KC, McCune RA, et al. Amount and distribution of 5-methylcytosine in human DNA from different types of tissues of cells. Nucleic Acids Res. 1982;10:2709–21.
Fullgrabe J, Kavanagh E, Joseph B. Histone onco-modifications. Oncogene. 2011;30:3391–403.
Jones PA, Baylin SB. The epigenomics of cancer. Cell. 2007;128:683–92.
Stearns V, Zhou Q, Davidson NE. Epigenetic regulation as a new target for breast cancer therapy. Cancer Invest. 2007;25:659–65.
Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009;462:315–22.
Fackler MJ, McVeigh M, Mehrotra J, Blum MA, Lange J, Lapides A, et al. Quantitative multiplex methylation-specific PCR assay for the detection of promoter hypermethylation in multiple genes in breast cancer. Cancer Res. 2004;13:4442–52.
Jeronimo C, Monteiro P, Henrique R, Dinis-Ribeiro M, Costa I, Costa VL, et al. Quantitative hypermethylation of a small panel of genes augments the diagnostic accuracy in fine-needle aspirate washings of breast lesions. Breast Cancer Res Treat. 2008;1:27–34.
Radpour R, Kohler C, Haghighi MM, Fan AX, Holzgreve W, et al. Methylation profiles of 22 candidate genes in breast cancer using high-throughput MALDI-TOF mass array. Oncogene. 2009;28:2969–78.
Nimmrich I, Sieuwerts AM, Meijer-van Gelder ME, Schwope I, Bolt-de Vries J, Harbeck N, et al. DNA hypermethylation of PITX2 is a marker of poor prognosis in untreated lymph node-negative hormone receptor-positive breast cancer patients. Breast Cancer Res Treat. 2008;3:429–37.
Sunami E, Shinozaki M, Sim MS, Nguyen SL, Vu AT, Giuliano AE, et al. Estrogen receptor and HER2/neu status affect epigenetic differences of tumor-related genes in primary breast tumors. Breast Cancer Res. 2008;3:R46.
Widschwendter M, Siegmund KD, Muller HM, Fiegl H, Marth C, Muller-Holzner E, et al. Association of breast cancer DNA methylation profiles with hormone receptor status and response to tamoxifen. Cancer Res. 2004;64:3807–13.
Li L, Lee KM, Han W, Choi JY, Lee JY, Kang GH, et al. Estrogen and progesterone receptor status affect genome-wide DNA methylation profile in breast cancer. Hum Mol Genet. 2010;21:4273–7.
Fang F, Turcan S, Rimner A, Kaufman A, Giri D, Morris LG, et al. Breast cancer methylomes establish an epigenomic foundation for metastasis. Sci Transl Med. 2011;3:75ra25.
Leon SA, Shapiro B, Sklaroff DM, Yaros MJ. Free DNA in the serum of cancer patients and the effect of therapy. Cancer Res. 1977;37:646–50.
Kloten V, Becker B, Winner K, Schrauder MG, Fasching PA, Anzeneder T, et al. Promoter hypermethylation of the tumor-suppressor genes ITIH5, DKK3, and RASSF1A as novel biomarkers for blood-based breast cancer screening. Breast Cancer Res. 2013;1:R4.
Chimonidou M, Tzitzira A, Strati A, Sotiropoulou G, Sfikas C, Malamos N, et al. CST6 promoter methylation in circulating cell-free DNA of breast cancer patients. Clin Biochem. 2013;3:235–40.
Chimonidou M, Strati A, Malamos N, Georgoulias V, Lianidou ES. SOX17 promoter methylation in circulating tumor cells and matched cell-free DNA isolated from plasma of patients with breast cancer. Clin Chem. 2013;1:270–9.
Guerrero-Preston R, Hadar T, Ostrow KL, Soudry E, Echenique M, Ili-Gangas C, et al. Differential promoter methylation of kinesin family member 1a in plasma is associated with breast cancer and DNA repair capacity. Oncol Rep. 2014;32:505–12.
Yang R, Pfütze K, Zucknick M, Sutter C, Wappenschmidt B, Marme F, et al. DNA methylation array analyses identified breast cancer associated HYAL2 methylation in peripheral blood. Int J Cancer. 2015;136:1845–55.
Elsheikh SE, Green AR, Rakha EA, Powe DG, Ahmed RA, Collins HM, et al. Global histone modifications in breast cancer correlate with tumor phenotypes, prognostic factors, and patient outcome. Cancer Res. 2009;9:3802–9.
Yokoyama Y, Matsumoto A, Hieda M, Shinchi Y, Ogihara E, Hamada M, et al. Loss of histone H4K20 trimethylation predicts poor prognosis in breast cancer and is associated with invasive activity. Breast Cancer Res. 2014;3:R66.
Dhingra T, Mittal K, Sarma GS. Analytical Techniques for DNA Methylation–An Overview. Curr Pharm Anal. 2014;1:71–85.
Szyf M. DNA methylation signatures for breast cancer classification and prognosis. Genome Med. 2012;3:26.
Gitan RS, Shi H, Chen CM, Yan PS, Huang TH. Methylation-specific oligonucleotide microarray: a new potential for high-throughput methylation analysis. Genome Res. 2002;1:158–64.
Huang Y, Pastor WA, Shen Y, Tahiliani M, Liu DR, Rao A. The behaviour of 5-hydroxymethylcytosine in bisulfite sequencing. PLoS ONE. 2010;5:e8888.
Huang TH, Perry MR, Laux DE. Methylation profiling of CpG islands in human breast cancer cells. Hum Mol Genet. 1999;3:459–70.
Yan PS, Chen CM, Shi H, Rahmatpanah F, Wei SH, Huang TH. Applications of CpG island microarrays for high-throughput analysis of DNA methylation. J Nutr. 2002;132(8 Suppl):S2430–4.
Schumacher A, Kapranov P, Kaminsky Z, et al. Microarray-based DNA methylation profiling: technology and applications. Nucleic Acids Res. 2006;2:528–42.
Mohn F, Weber M, Schübeler D, Roloff TC. Methylated DNA immunoprecipitation (MeDIP). Methods Mol Biol. 2009;507:55–64.
Zhang B, Zhou Y, Lin N, Lowdon RF, Hong C, Nagarajan RP, et al. Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm. Genome Res. 2013;23:1522–40.
Zhang M, Smith A. Challenges in understanding genome-wide DNA methylation. J Comput Sci Technol. 2010;1:26–34.
Bhasin M, Zhang H, Reinherz E, Reche P. Prediction of methylated CpGs in DNA sequences using a support vector machine. FEBS Lett. 2005;579:4302–8.
Lu L, Lin K, Qian Z, Li H, Cai Y, Li Y. Predicting DNA methylation status using word composition. J Biomedical Science and Engineering. 2010;3:672–6.
Ali I, Seker H. Detailed methylation prediction of CpG islands on human chromosome 21. 10th WSEAS International Conference on Mathematics and Computers. In: Biology and Chemistry. 2009. p. 147–52.
Fan S, Zhang M, Zhang X. Histone methylation marks play important roles in predicting the methylation status of CpG islands. Biochem Biophys Res Commun. 2008;374:559–64.
Previti C, Harari O, Zwir I, del Val C. Profile analysis and prediction of tissue-specific CpG island methylation classes. BMC Bioinformatics. 2009;10:116.
Bell JT, Pai AA, Pickrell JK, Gaffney DJ, Pique-Regi R, Degner JF, et al. DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines. Genome Biol. 2011;12:R10.
Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, Burger M, et al. DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006;38:1378–85.
Zhang W, Spector TD, Deloukas P, Bell JT, Engelhardt BE. Predicting genome-wide DNA methylation using methylation marks, genomic position, and DNA regulatory elements. Genome Biol. 2015;16:14.
Du P, Bourgon R. methyAnalysis: DNA methylation data analysis and visualization. R package version 1.10.0. 2014.
Barfield RT, Kilaru V, Smith AK, Conneely KN. CpGassoc: an R function for analysis of DNA methylation microarray data. Bioinformatics. 2012;9:1280–1.
Assenov Y, Mueller F, Lutsik P, Walter J, Lengauer T, Bock C. Compehensive Analysis of DNA Methylation Data with RnBeads. Nat Methods. 2014;11:1138–40.
Wang D, Yan L, Hu Q, Sucheston LE, Higgins MJ, Ambrosone CB, et al. IMA: an R package for high-throughput analysis of Illumina's 450K Infinium methylation data. Bioinformatics. 2012;5:729–30.
Price EM, Cotton AM, Lam LL, Farré P, Emberly E, Brown CJ, et al. Additional annotation enhances potential for biologically-relevant analysis of the illumina infinium humanmethylation450 beadchip array. Epigenetics Chromatin. 2013;1:4.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP. Hansen KD and Irizarry RA Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA Methylation microarrays. Bioinformatics. 2014;10:1363–9.
Jaffe AE, Murakami P, Lee H, Leek JT, Fallin DM, Feinberg AP, et al. Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int J Epidemiol. 2012;1:200–9.
Kanduri M, Cahill N, Göransson H, Enström C, Ryan F, Isaksson A, et al. Differential genome-wide array-based methylation profiles in prognostic subsets of chronic lymphocytic leukemia. Blood. 2010;2:296–305.
Wessely F, Emes RD. Identification of DNA methylation biomarkers from Infinium arrays. Front Genet. 2012;3:161.
Wilhelm-Benartzi CS, Koestler DC, Karagas MR, Flanagan JM, Christensen BC, Kelsey KT, et al. Review of processing and analysis methods for DNA methylation array data. Br J Cancer. 2013;6:1394–402.
Phipson B, Oshlack A. DiffVar: a new method for detecting differential variability with application to methylation in cancer and aging. Genome Biol. 2014;9:465.
Goeman JJ, Bühlmann P. Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007;23(8):980–7.
Zhao N, Bell DA, Maity A, Staicu AM, Joubert BR, London SJ, et al. Global analysis of methylation profiles from high resolution CpG data. Genet Epidemiol. 2012;2:53–64.
Westfall PH, Stanley Young S. Resampling-Based Multiple Testing: Examples and Methods for p-Value Adjustment. New York: Wiley-Interscience; 1993.
Teng M, Wang Y, Kim S, Li L, Shen C, Wang G, et al. Empirical bayes model comparisons for differential methylation analysis. Comp Funct Genomics. 2012;2012:376706.
Li D, Xie Z, Pape ML, Dye T. An evaluation of statistical methods for DNA methylation microarray data analysis. BMC Bioinformatics. 2015;16:217.
Subramaniam D, Thombre R, Dhar A, Anant S. DNA methyltransferases: a novel target for prevention and therapy. Front Oncol. 2014;4:80.
Appleton K, Mackay HJ, Judson I, Plumb JA, McCormick C, Strathdee G, et al. Phase I and pharmacodynamic trial of the DNA methyltransferase inhibitor decitabine and carboplatin in solid tumors. J Clin Oncol. 2007;25:4603–9.
Pouliot MC, Labrie Y, Diorio C, Durocher F. The Role of Methylation in Breast Cancer Susceptibility and Treatment. Anticancer Res. 2015;9:4569–74.
Chen M, Shabashvili D, Nawab A, Yang SX, Dyer LM, Brown KD, et al. DNA methyltransferase inhibitor, zebularine, delays tumor growth and induces apoptosis in a genetically engineered mouse model of breast cancer. Mol Cancer Ther. 2012;11:370–82.
Billam M, Sobolewski MD, Davidson NE. Effects of a novel DNA methyltransferase inhibitor zebularine on human breast cancer cells. Breast Cancer Res Treat. 2010;120:581–92.
Yoo CB, Jeong S, Egger G, Liang G, Phiasivongsa P, Tang C, et al. Delivery of 5-aza-2’-deoxycytidine to cells using oligodeoxynucleotides. Cancer Res. 2007;67:6400–8.
Nie J, Liu L, Li X, Han W. Decitabine, a new star in epigenetic therapy: the clinical application and biological mechanism in solid tumors. Cancer Lett. 2014;354:12–20.
Marson CM. Histone deacetylase inhibitors: design, structure-activity relationships and therapeutic implications for cancer. Anticancer Agents Med Chem. 2009;9:661–92.
Munster PN, Thurn KT, Thomas S, Raha P, Lacevic M, Miller A, et al. A phase II study of the histone deacetylase inhibitor vorinostat combined with tamoxifen for the treatment of patients with hormone therapy-resistant breast cancer. Br J Cancer. 2011;104:1828–35.
Yardley DA, Ismail-Khan RR, Melichar B, Lichinitser M, Munster PN, Klein PM, et al. Randomized phase II, double-blind, placebo-controlled study of exemestane with or without entinostat in postmenopausal women with locally recurrent or metastatic estrogen receptor-positive breast cancer progressing on treatment with a nonsteroidal aromatase inhibitor. J Clin Oncol. 2013;17:2128–35.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Pasquinelli AE, Hunter S, Bracht J. MicroRNAs: a developing story. Curr Opin Genet Dev. 2005;15:200–5.
Iorio MV, Croce CM. MicroRNA dysregulation in cancer: diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol Med. 2012;4(3):143–59.
Blenkiron C, Goldstein LD, Thorne NP, Spiteri I, Chin SF, Dunning MJ, et al. MicroRNA expression profiling of human breast cancer identifies new markers of tumor subtype. Genome Biol. 2007;8:R214.
Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, et al. MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005;65:7065–70.
Lowery AJ, Miller N, Devaney A, McNeill RE, Davoren PA, Lemetre C, et al. MicroRNA signatures predict oestrogen receptor, progesterone receptor and HER2/neu receptor status in breast cancer. Breast Cancer Res. 2009;11:R27.
Mattie MD, Benz CC, Bowers J, Sensinger K, Wong L, Scott GK, et al. Optimized high-throughput microRNA expression profiling provides novel biomarker assessment of clinical prostate and breast cancer biopsies. Mol Cancer. 2006;5:24.
Gregory PA, Bracken CP, Bert AG, Goodall GJ. MicroRNAs as regulators of epithelial–mesenchymal transition. Cell Cycle. 2008;7:3112–8.
Sempere LF, Christensen M, Silahtaroglu A, Bak M, Heath CV, Schwartz G, et al. Altered MicroRNA expression confined to specific epithelial cell subpopulations in breast cancer. Cancer Res. 2007;67:11612–20.
Zhou M, Liu Z, Zhao Y, Ding Y, Liu H, Xi Y, et al. MicroRNA-125b confers the resistance of breast cancer cells to paclitaxel through suppression of pro-apoptotic Bcl-2 antagonist killer 1 (Bak1) expression. J Biol Chem. 2010;285:21496–2507.
Taylor MA, Sossey-Alaoui K, Thompson CL, Danielpour D, Schiemann WP. TGF-beta upregulates miR-181a expression to promote breast cancer metastasis. J Clin Invest. 2013;123:150–63.
Smith AL, Iwanaga R, Drasin DJ, et al. The miR-106b-25 cluster targets Smad7, activates TGF-beta signaling, and induces EMT and tumor initiating cell characteristics downstream of Six1 in human breast cancer. Oncogene. 2012;31:5162–71.
de Souza Rocha Simonini P, Breiling A, Gupta N, et al. Epigenetically deregulated microRNA-375 is involved in a positive feedback loop with estrogen receptor alpha in breast cancer cells. Cancer Res. 2010;70:9175–84.
Reddy SD, Ohshiro K, Rayala SK, Kumar R. MicroRNA-7, a homeobox D10 target, inhibits p21-activated kinase 1 and regulates its functions. Cancer Res. 2008;68:8195–200.
Webster RJ, Giles KM, Price KJ, Zhang PM, Mattick JS, Leedman PJ. Regulation of epidermal growth factor receptor signaling in human cancer cells by microRNA-7. J Biol Chem. 2009;284:5731–41.
Yu Z, Willmarth NE, Zhou J, Katiyar S, Wang M, Liu Y, et al. microRNA 17/20 inhibits cellular invasion and tumor metastasis in breast cancer by heterotypic signaling. Proc Natl Acad Sci U S A. 2010;107:8231–6.
Yu Z, Wang C, Wang M, Li Z, Casimiro MC, Liu M, et al. A cyclin D1/microRNA 17/20 regulatory feed-back loop in control of breast cancer cell proliferation. J Cell Biol. 2008;182:509–17.
Xu D, Takeshita F, Hino Y, Fukunaga S, Kudo Y, Tamaki A, et al. miR-22 represses cancer progres-sion by inducing cellular senescence. J Cell Biol. 2011;193:409–24.
Patel JB, Appaiah HN, Burnett RM, Bhat-Nakshatri P, Wang G, Mehta R, et al. Control of EVI-1 oncogene expression in metastatic breast cancer cells through microRNA miR-22. Oncogene. 2011;30:1290–301.
Pandey DP, Picard D. miR-22 inhibits estrogen signaling by directly targeting the estrogen receptor alpha mRNA. Mol Cell Biol. 2009;29:3783–90.
Wu F, Zhu S, Ding Y, Beck WT, Mo YY. MicroRNA- mediated regulation of Ubc9 expression in cancer cells. Clin Cancer Res. 2009;15:1550–7.
Yu F, Deng H, Yao H, Liu Q, Su F, Song E. Mir-30 reduction maintains self-renewal and inhibits apoptosis in breast tumor-initiating cells. Oncogene. 2010;29:4194–204.
Valastyan S, Reinhardt F, Benaich N, Calogrias D, Szász AM, Wang ZC, et al. A pleiotropically acting microRNA, miR-31, inhibits breast cancer metastasis. Cell. 2009;137:1032–46.
Valastyan S, Chang A, Benaich N, Reinhardt F, Weinberg RA. Concurrent suppression of integrin alpha5, radixin, and RhoA phenocopies the effects of miR-31 on metastasis. Cancer Res. 2010;70:5147–54.
Valastyan S, Benaich N, Chang A, Reinhardt F, Weinberg RA. Concomitant suppression of three target genes can explain the impact of a microRNA on metastasis. Genes Dev. 2009;23:2592–7.
Harris TA, Yamakuchi M, Ferlito M, Mendell JT, Lowenstein CJ. MicroRNA-126 regulates endothelial expression of vascular cell adhesion molecule 1. Proc Natl Acad Sci U S A. 2008;105:1516–21.
Sachdeva M, Zhu S, Wu F, Wu H, Walia V, Kumar S, et al. p53 represses c-Myc through induction of the tumor suppressor miR-145. Proc Natl Acad Sci U S A. 2009;106:3207–12.
Hurst DR, Edmonds MD, Scott GK, Benz CC, Vaidya KS, Welch DR. Breast cancer metastasis suppressor 1 up-regulates miR-146, which suppresses breast cancer metastasis. Cancer Res. 2009;69:1279–83.
Li XF, Yan PJ, Shao ZM. Downregulation of miR-193b contributes to enhance urokinase-type plasminogen activator (uPA) expression and tumor progression and invasion in human breast cancer. Oncogene. 2009;28:3937–48.
Wu H, Zhu S, Mo YY. Suppression of cell growth and invasion by miR-205 in breast cancer. Cell Res. 2009;19:439–48.
Song G, Zhang Y, Wang L. MicroRNA-206 targets notch3, activates apoptosis, and inhibits tumor cell migration and focus formation. J Biol Chem. 2009;284:31921–7.
Edmonds MD, Hurst DR, Vaidya KS, Stafford LJ, Chen D, Welch DR. Breast cancer metastasis suppressor 1 coordinately regulates metastasis-associated microRNA expression. Int J Cancer. 2009;125:1778–85.
Li QQ, Chen ZQ, Cao XX, Xu JD, Xu JW, Chen YY, et al. Involvement of NF-κB/miR-448 regulatory feedback loop in chemotherapy-induced epithelial-mesenchymal transition of breast cancer cells. Cell Death Differ. 2011;18:16–25.
Reddy SD, Pakala SB, Ohshiro K, Rayala SK, Kumar R. MicroRNA-661, a c/EBPalpha target, inhibits metastatic tumor antigen 1 and regulates its functions. Cancer Res. 2009;69:5639–42.
Yu F, Yao H, Zhu P, Zhang X, Pan Q, Gong C, et al. let-7 regulates self renewal and tumorigenicity of breast cancer cells. Cell. 2007;131:1109–23.
Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A. 2008;105:10513–8.
Taylor DD, Gercel-Taylor C. MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol Oncol. 2008;110:13–21.
Michael A, Bajracharya SD, Yuen PS, Zhou H, Star RA, Illei GG, et al. Exosomes from human saliva as a source of microRNA biomarkers. Oral Dis. 2010;16:34–8.
Park NJ, Zhou H, Elashoff D, Henson BS, Kastratovic DA, Abemayor E, et al. Salivary microRNA: discovery, characterization, and clinical utility for oral cancer detection. Clin Cancer Res. 2009;15:5473–7.
Xie Y, Todd NW, Liu Z, Zhan M, Fang H, Peng H, et al. Altered miRNA expression in sputum for diagnosis of nonsmall cell lung cancer. Lung Cancer. 2010;67:170–6.
Yu L, Todd NW, Xing L, Xie Y, Zhang H, Liu Z, et al. Early detection of lung adenocarcinoma in sputum by a panel of microRNA markers. Int J Cancer. 2010;127:2870–8.
Lodes MJ, Caraballo M, Suciu D, Munro S, Kumar A, Anderson B. Detection of cancer with serum miRNAs on an oligonucleotide microarray. PLoS One. 2009;4:e6229.
Madhavan D, Cuk K, Burwinkel B, Yang R. Cancer diagnosis and prognosis decoded by blood-based circulating microRNA signatures. Front Genet. 2013;4:116.
Si H, Sun X, Chen Y, Cao Y, Chen S, Wang H, et al. Circulating microRNA-92a and microRNA-21 as novel minimally invasive biomarkers for primary breast cancer. J Cancer Res Clin Oncol. 2013;139:223–9.
Mar-Aguilar F, Mendoza-Ramirez JA, Malagon-Santiago I, Espino-Silva PK, Santuario-Facio SK, Ruiz-Flores P, et al. Serum circulating microRNA profiling for identification of potential breast cancer biomarkers. Dis Markers. 2013;34:163–9.
Zeng RC, Zhang W, Yan XQ, Ye ZQ, Chen ED, Huang DP, et al. Down-regulation of miRNA-30a in human plasma is a novel marker for breast cancer. Med Oncol. 2013;30:477.
Madhavan D, Zucknick M, Wallwiener M, Cuk K, Modugno C, Scharpff M, et al. Circulating miRNAs as surrogate markers for circulating tumor cells and prognostic markers in metastatic breast cancer. Clin Cancer Res. 2012;18:5972–82.
Chen W, Cai F, Zhang B, Barekati Z, Zhong XY. The level of circulating miRNA-10b and miRNA-373 in detecting lymph node metastasis of breast cancer: potential biomarkers. Tumour Biol. 2012;34:455–62.
Wang H, Tan G, Dong L, Cheng L, Li K, Wang Z, et al. Circulating MiR-125b as a marker predicting chemoresistance in breast cancer. PLoS ONE. 2012;7:e34210.
Sun Y, Wang M, Lin G, Sun S, Li X, Qi J, et al. Serum microRNA-155 as a potential biomarker to track disease in breast cancer. PLoS ONE. 2012;7:e47003.
Thompson RC, Deo M, Turner DL. Analysis of microRNA expression by in situ hybridization with RNA oligonucleotide probes. Methods. 2007;43(2):153–61.
Stark MS, Tyagi DJ, Nancarrow GM, Boyle AL, Cook DC, Whiteman PG, et al. Characterization of the Melanoma miRNAome by Deep Sequencing. PLoS One. 2010;5(3):e9685.
Farazi TA, Horlings HM, Ten Hoeve JJ, Mihailovic A, Halfwerk H, Morozov P, et al. MicroRNA sequence and expression analysis in breast tumors by deep sequencing. Cancer Res. 2011;71:4443–53.
Wu Q, Lu Z, Li H, Lu J, Guo L, Ge Q. Next-generation sequencing of microRNAs for breast cancer detection. J Biomed Biotechnol. 2011;2011:597145.
Hu Z, Chen X, Zhao Y, Tian T, Jin G, Shu Y, et al. Serum microRNA signatures identified in a genome-wide serum microRNA expression profiling predict survival of non-small-cell lung cancer. J Clin Oncol Off J Am Soc Clin Oncol. 2010;28:1721–6.
Xu JZ, Wong CW. Hunting for robust gene signature from cancer profiling data: sources of variability, different interpretations, and recent methodological developments. Cancer Lett. 2010;296:9–16.
Peltier HJ, Latham GJ. Normalization of microRNA expression levels in quantitative RT-PCR assays: identification of suitable reference RNA targets in normal and cancerous human solid tissues. RNA. 2008;14:844–52.
Maire G, Martin JW, Yoshimoto M, Chilton-MacNeill S, Zielenska M, Squire JA. Analysis of miRNA-gene expression-genomic profiles reveals complex mechanisms of microRNA deregulation in osteosarcoma. Cancer Genetics. 2011;204(3):138–46.
Wang C, Su Z, Sanai N, et al. microRNA expression profile and differentially-expressed genes in prolactinomas following bromocriptine treatment. Oncol Rep. 2012;27(5):1312–20.
Lai EC, Wiel C, Rubin GM. Complementary miRNA pairs suggest a regulatory role for miRNA:miRNA duplexes. RNA. 2004;10(2):171–5.
Yu J, Liu F, Yin P, et al. Integrating miRNA and mRNA expression profiles in response to heat stress-induced injury in rat small intestine. Funct Integr Genomics. 2011;11(2):203–13.
Liu B, Liu L, Tsykin A, et al. Identifying functional miRNA mRNA regulatory modules with correspondence latent dirichlet allocation. Bioinformatics. 2010;26(24):3105–11.
Nielsen JA, Lau P, Maric D, Barker JL, Hudson LD. Integrating microRNA and mRNA expression profiles of neuronal progenitors to identify regulatory networks underlying the onset of cortical neurogenesis. BMC Neurosci. 2009;10:98.
Li L, Xu J, Yang D, Tan X, Wang H. Computational approaches for microRNA studies: a review. Mamm Genome. 2010;21(1–2):1–12.
Grad Y, Aach J, Hayes GD, Reinhart BJ, Church GM, Ruvkun G, et al. Computational and experimental identification of C. elegans microRNAs. Mol. Cell. 2003;11:1253–63.
Lim LP, Lau NC, Weinstein EG, Abdelhakim A, Yekta S, Rhoades MW, et al. The microRNAs of Caenorhabditis elegans. Genes Dev. 2003;17:991–1008.
Lai EC, Tomancak P, Williams RW, Rubin GM. Computational identification of Drosophila microRNA genes. Genome Biol. 2003;4:R42.
Lee RC, Ambros V. An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001;294:862–4.
Berezikov E, Guryev V, Van DE, Belt J, Wienholds E, Plasterk RH, et al. Phylogenetic shadowing and computational identification of human microRNA genes. Cell. 2005;120:21–4.
Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, et al. Identification of hundreds of conserved and nonconserved human microRNAs. Nat Genet. 2005;37:766–70.
Hofacker IL, Fontana W, Stadler PF, Bonhoeffer LS, Tacker M, Schus-Ter P. Fast folding and comparison of RNA secondary structures. Monatshefte fÄur Chemie/Chemial Monthly. 1994;125(2):167–88.
Lindow M, Gorodkin J. Principles and limitations of computational microRNA gene and target finding. DNA Cell Biol. 2007;26:339–51.
Allmer J, Yousef M. Computational methods for ab initio detection of microRNAs. Front Genet. 2012;3:209.
Wang C, Ding C, Meraz RF, Holbrook SR. PSoL: a positive sample only learning algorithm for finding non-coding RNA genes. Bioinformatics. 2006;22:2590–6.
Yousef M, Jung S, Showe LC, Showe MK. Learning from positive examples when the negative class is undetermined – microRNA gene identification. Algorithms Mol Biol. 2008;3:2.
Hertel J, Stadler PF. Hairpins in a Haystack: recognizing microRNA precursors in comparative genomics data. Bioinformatics. 2006;22:e197–202.
Huang TH, Fan B, Rothschild MF, Hu ZL, Li K, Zhao SH. MiRFinder: an improved approach and software implementation for genome-wide fast microRNA precursor scans. BMC Bioinformatics. 2007;8:341.
Nam JW, Shin KR, Han J, Lee Y, Kim VN, Zhang BT. Human microRNA prediction through a probabilistic co-learning model of sequence and structure. Nucleic Acids Res. 2005;33:3570–81.
Terai G, Komori T, Asai K, Kin T. miRRim: a novel system to find conserved miRNAs with high sensitivity and specificity. RNA. 2007;13:2081–90.
Oulas A, Boutla A, Gkirtzou K, Reczko M, Kalantidis K, Poirazi P. Prediction of novel microRNA genes in cancer-associated genomic regions – a combined computational and experimental approach. Nucleic Acids Res. 2009;37:3276–87.
Kadri S, Hinman V, Benos PV. HHMMiR: efficient de novo prediction of microRNAs using hierarchical hidden Markov models. BMC Bioinformatics. 2009;10 Suppl 1:S35.
Yousef M, Nebozhyn M, Shatkay H, Kanterakis S, Showe LC, Showe MK. Combining multi-species genomic data for microRNA identification using a Naive Bayes classifier. Bioinformatics. 2006;22:1325–34.
Bentwich I. Prediction and validation of microRNAs and their targets. FEBS Lett. 2005;579:5904–10.
Xue C, Li F, He T, Liu GP, Li Y, Zhang X. Classification of real and pseudo microRNA precursors using local structure-sequence features and support vector machine. BMC Bioinformatics. 2005;6:310.
Ng KL, Mishra SK. De novo SVM classification of precursor microRNAs from genomic pseudo hairpins using global and intrinsic folding measures. Bioinformatics. 2007;23:1321–30.
Sewer A, Paul N, Landgraf P, Aravin A, Pfeffer S, Brownstein MJ, et al. Identification of clustered microRNAs using an ab initio prediction method. BMC Bioinformatics. 2005;6:267.
Zhao X, Pan F, Holt CM, Lewis AL, Lu JR. Controlled delivery of antisense oligonucleotides: a brief review of current strategies. Expert Opin Drug Deliv. 2009;6:673–86.
Dias N, Stein CA. Antisense oligonucleotides: basic concepts and mechanisms. Mol Cancer Ther. 2002;1:347–55.
Samantarrai D, Dash S, Chhetri B, Mallick B. Genomic and epigenomic cross-talks in the regulatory landscape of miRNAs in breast cancer. Mol Cancer Res. 2013;11(4):315–28.
Zoon CK, Starker EQ, Wilson AM, Emmert-Buck MR, Libutti SK, Tangrea MA. Current molecular diagnostics of breast cancer and the potential incorporation of microRNA. Expert Rev Mol Diagn. 2009;9(5):455–67.
van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415(6871):530–6.
van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, et al. A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002;347(25):1999–2009.
Buyse M, Loi S, van't Veer L, Viale G, Delorenzi M, Glas AM, et al. TRANSBIG Consortium. Validation and clinical utility of a 70-gene prognostic signature for women with node-negative breast cancer. J Natl Cancer Inst. 2006;98(17):1183–92.
Bueno-de-Mesquita JM, van Harten WH, Retel VP, van't Veer LJ, van Dam FS, Karsenberg K, et al. Use of 70-gene signature to predict prognosis of patients with node-negative breast cancer: a prospective community-based feasibility study (RASTER). Lancet Oncol. 2007;8(12):1079–87.
Bueno-de-Mesquita JM, Linn SC, Keijzer R, Wesseling J, Nuyten DS, van Krimpen C, et al. Validation of 70-gene prognosis signature in node-negative breast cancer. Breast Cancer Res Treat. 2009;117(3):483–95.
Wittner BS, Sgroi DC, Ryan PD, Bruinsma TJ, Glas AM, Male A, et al. Analysis of the MammaPrint breast cancer assay in a predominantly postmenopausal cohort. Clin Cancer Res. 2008;14(10):2988–93.
Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, et al. A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004;351(27):2817–26.
Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, et al. Multi-gene RT-PCR assay for predicting recurrence in node negative breast cancer patients — NSABP studies B-20 and B-14. Breast Cancer Res Treat. 2003;82:A16.
Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27(8):1160–7.
Desmedt C, Giobbie-Hurder A, Neven P, Paridaens R, Christiaens MR, Smeets A, et al. The Gene expression Grade Index: a potential predictor of relapse for endocrine-treated breast cancer patients in the BIG 1–98 trial. BMC Med Genomics. 2009;2:40.
Bartlett JM, Thomas J, Ross DT, Seitz RS, Ring BZ, Beck RA, et al. Mammostrat as a tool to stratify breast cancer patients at risk of recurrence during endocrine therapy. Breast Cancer Res. 2010;12(4):R47.
Ring BZ, Seitz RS, Beck R, Shasteen WJ, Tarr SM, Cheang MC, et al. Novel prognostic immunohistochemical biomarker panel for estrogen receptor-positive breast cancer. J Clin Oncol. 2006;24(19):3039–47.
Streit S, Michalski CW, Erkan M, Kleeff J, Friess H. Northern blot analysis for detection and quantification of RNA in pancreatic cancer cells and tissues. Nat Protoc. 2009;4(1):37–43.
Bustin SA. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000;25(2):169–93.
Valasek MA, Repa JJ. The power of real-time PCR. Adv Physiol Educ. 2005;29(3):151–9.
Costa C, Giménez-Capitán A, Karachaliou N, Rosell R. Comprehensive molecular screening: from the RT-PCR to the RNA-seq. Trans Lung Cancer Res. 2013;2(2):87–91.
Taniguchi M, Miura K, Iwao H, Yamanaka S. Quantitative assessment of DNA microarrays—comparison with Northern blot analyses. Genomics. 2001;71(1):34–9.
Alwine JC, Kemp DJ, Stark GR. Method for detection of specific RNAs in agarose gels by transfer to diazobenzyloxymethyl-paper and hybridization with DNA probes. Proc Natl Acad Sci U S A. 1977;74(12):5350–4.
Pollock JD. Gene expression profiling: methodological challenges, results, and prospects for addiction research. Chem Phys Lipids. 2002;121(1–2):241–56.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.
Velculescu VE, Zhang L, Vogelstein B, Kinzler KW. Serial analysis of gene expression. Science. 1995;270(5235):484–7.
Kodzius R, Kojima M, Nishiyori H, Nakamura M, Fukuda S, Tagami M, et al. CAGE: cap analysis of gene expression. Nat Methods. 2006;3(3):211–22.
Brenner S, Johnson M, Bridgham J, Golda G, Lloyd DH, Johnson D, et al. Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nat Biotechnol. 2000;18(6):630–4.
Medvedev P, Stanciu M, Brudno M. Computational methods for discovering structural variation with next-generation sequencing. Nat Methods. 2009;6(11 Suppl):S13–20.
van de Wiel MA, Picard F, van Wieringen WN, Ylstra B. Preprocessing and downstream analysis of microarray DNA copy number profiles. Brief Bioinform. 2011;12(1):10–21.
Saeys Y, Inza I, Larrañaga P. A review of feature selection techniques in bioinformatics. Bioinformatics. 2007;23(19):2507–17.
Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999;286(5439):531–7.
Bolón-Canedo V, Sánchez-Maroño N, Alonso-Betanzos A, Benítez JM, Herrera F. A review of microarray datasets and applied feature selection methods. Inf Sci. 2014;282:111–35.
Kumar AP, Valsala P. Feature Selection for high Dimensional DNA Microarray data using hybrid approaches. Bioinformation. 2013;9(16):824–8.
Chuang LY, Yang CS, Wu KC, Yang CH. Correlation-based gene selection and classification using Taguchi-BPSO. Methods Inf Med. 2010;49(3):254