Microarray and deep sequencing cross-platform analysis of the mirRNome and isomiR variation in response to epidermal growth factor
- Franc Llorens1, 2, 3, 4, 5,
- Manuela Hummel1, 2,
- Lorena Pantano1, 2, 6, 7,
- Xavier Pastor1, 2, 7,
- Ana Vivancos1, 2, 8,
- Ester Castillo1, 2,
- Heidi Mattlin1, 2, 9,
- Anna Ferrer1, 2,
- Matthew Ingham1, 2, 9,
- Marc Noguera7, 10,
- Robert Kofler1, 2, 11, 12,
- Juliane C Dohm1, 2, 11,
- Raquel Pluvinet7,
- Mònica Bayés1, 2, 9,
- Heinz Himmelbauer1, 2,
- José Antonio del Rio3, 4, 5,
- Eulàlia Martí1, 2, 6 and
- Lauro Sumoy1, 2, 7Email author
© Llorens et al.; licensee BioMed Central Ltd. 2013
Received: 8 October 2012
Accepted: 20 May 2013
Published: 1 June 2013
Epidermal Growth Factor (EGF) plays an important function in the regulation of cell growth, proliferation, and differentiation by binding to its receptor (EGFR) and providing cancer cells with increased survival responsiveness. Signal transduction carried out by EGF has been extensively studied at both transcriptional and post-transcriptional levels. Little is known about the involvement of microRNAs (miRNAs) in the EGF signaling pathway. miRNAs have emerged as major players in the complex networks of gene regulation, and cancer miRNA expression studies have evidenced a direct involvement of miRNAs in cancer progression.
In this study, we have used an integrative high content analysis approach to identify the specific miRNAs implicated in EGF signaling in HeLa cells as potential mediators of cancer mediated functions. We have used microarray and deep-sequencing technologies in order to obtain a global view of the EGF miRNA transcriptome with a robust experimental cross-validation. By applying a procedure based on Rankprod tests, we have delimited a solid set of EGF-regulated miRNAs. After validating regulated miRNAs by reverse transcription quantitative PCR, we have derived protein networks and biological functions from the predicted targets of the regulated miRNAs to gain insight into the potential role of miRNAs in EGF-treated cells. In addition, we have analyzed sequence heterogeneity due to editing relative to the reference sequence (isomiRs) among regulated miRNAs.
We propose that the use of global genomic miRNA cross-validation derived from high throughput technologies can be used to generate more reliable datasets inferring more robust networks of co-regulated predicted miRNA target genes.
EGF signaling activates an extensive network of signal transduction pathways leading to: 1) Activation or inhibition of genes regulating DNA synthesis, cell proliferation and pro- or anti- apoptotic pathways 2) rises in intracellular calcium levels, 3) increases in glycolysis and protein synthesis [1, 2]. EGF signaling provides cancer cells with increased pro-survival responses and deregulated activity of this network leads to diverse types of tumors [2, 3]. By targeting EGF survival pathways, cancer cells can be forced to undergo apoptosis or become sensitive to radiation or chemotherapy. Hence, targeted therapies that block EGF signaling have been successfully applied towards anti-cancer treatment .
Although the EGFR mediated signaling response is one of the best understood, questions such as how the specificity and dynamics of response to EGF is achieved or what are the regulatory feed-back mechanisms have been the subject of intense investigation to reach a global and integrative view of EGF-dependent cellular regulation [5–7]. Indeed, since EGF signaling pathways can be deregulated at many levels during cancer progression, a multi-faceted approach to target these pathways and their regulatory mechanisms for cancer treatment is required.
Over recent years, miRNAs have emerged as major players in the complex networks controlling gene regulation. These small, endogenous non-coding RNA molecules that regulate post-transcriptionally the expression of protein coding genes  have been implicated in various spectra of human pathobiology, such as cancer, metabolic disorders and infectious diseases ([9, 10] for review). Cancer miRNA expression studies have proven that miRNA expression profiles can classify tumors in a very robust fashion, and mutational and functional analyses have evidenced a direct role of miRNAs in cancer progression [11, 12]. Deregulation of miRNAs promotes carcinogenesis; hence these molecules can act as oncogenes or tumor suppressor genes . Moreover, both miRNA inhibition and activation show great promise in the treatment of various types of cancer, as well as viral and metabolic diseases. Aberrant gene expression is the main reason for miRNA dysfunction in cancer, which results in abnormal miRNA levels in tumor samples. miRNA germline and somatic mutations or polymorphisms in the protein-coding mRNAs targeted by miRNAs also contribute to cancer predisposition, initiation or progression . These findings have revealed that, besides considering the role of oncogenes and tumor suppressor protein coding genes, it will be essential to understand how miRNAs affect the responsiveness of cells to signaling molecules involved in such processes in order to advance in the knowledge of the mechanisms underlying malignant transformation . While the EGF pathway has been investigated in depth at the phosphorylation , and gene expression levels [6, 16] much less is known about how EGF may regulate miRNA expression or the role of such miRNAs in the regulation of the EGF-related cellular functions such as cell growth, proliferation and differentiation. Hayashi and collaborators  profiled the EGF-dependent microRNAs of the fetal mouse submandibular gland at embryonic day 13. Avraham and colleagues,  described a coordinated transcriptional program between miRNA and transcription factors revealed through an EGF treatment time course in HeLa cells. Similarly, a very recent study in lung cancer cell lines has addressed the miRNA response to EGFR inhibition by shRNA .
Optimal exploitation of genomics and bioinformatics technologies have provided many tools that have been successfully applied to the study of EGF signaling (establishing an elaborate model of an EGF-dependent transcriptional modulator network)  and to many types of cancers . In addition, ongoing progress in ultra-sequencing technologies has eased and expanded the possibilities to perform measurements of millions of molecules in a single assay, allowing for even more precise modeling.
In this study, we present the complete view of the EGF-induced miRNA transcription in HeLa cells which has been cross-validated using different array and ultra-sequencing technologies, including the analysis of the miRNA variability in the EGF-induced miRNome, a phenomenon that may specifically influence the mechanisms of gene silencing or gene targeting under both physiological and pathological conditions. In parallel, we have developed specific bioinformatics and statistical tools to assist in the handling and analysis of the vast amounts of data generated.
Reagents and antibodies
EGF and anti-Tubulin (1:10000) were purchased from Sigma. Anti phosho-ERK1/2 (1:2000), anti phosho-AKT (1:1000) were from Cell Signaling. U0126, AG1478 and wortmannin were from Calbiochem.
Cell culture and sample preparation
HeLa cells were cultured at 37°C in a 95/5 Air/CO2 water saturated atmosphere in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% heat inactivated fetal bovine serum (FBS), 2 mM L-glutamine and 100U/ml Penicillin/streptomycin. For treatments, the cells were transferred to 60 mm dishes and after 48 hours starved for 24 hours in DMEM containing 0.5% FBS. The cells were incubated (if indicated) with the protein kinase inhibitors U0126 (10 μM) or AG1478 (10 μM) for 30 min, and then stimulated with EGF (150 ng/ml) for the indicated times. Cells were harvested, washed twice with cold phosphate-buffered saline and lysed with either 2x Laemmli sample buffer (Sigma), for protein extraction, or miRVANA lysis/binding buffer (Ambion), for total RNA extraction following manufacturer’s instructions.
For transfections, HeLa cells were cultured as indicated above and transfected with either an empty plasmid (pcDNA), or with the plasmid encoding a constitutively active form of Ras (pcDNA-RasV12) with Lipofectamine 2000 (Invitrogen) following the manufacturer’s instructions. After transfection, cells were cultured for 24 h in DMEM containing 0.5% FBS.
Total RNA was quantified with NanoDrop ND-1000 followed by quality assessment with 2100 Bioanalyzer (Agilent Technologies) Nano 6000 assay according to the manufacturer’s protocol. Acceptable A260/A280 ratios were in the range 1.8-2.2. Acceptable rRNA ratio (28S/18S) had to be >0.9 and RIN (RNA Integrity Number) value has to be >8.0.
For Western blotting 50 μg of cell extracts from HeLa cells were subjected to 8-10% SDS-PAGE. Gels were transferred onto PVDF membranes and processed for specific immunodetection by ECL (Pharmacia) following manufacturer’s instructions using antibodies as indicated above.
One μg of total RNA from sample (individual EGF treated or untreated control) and reference (pool of three control samples) were labeled with Hy3™ and Hy5™ fluorescent label, respectively, using the miRCURY™ LNA Array power labeling kit (Exiqon, Denmark) following the procedure described by the manufacturer. The Hy3™-labeled samples and a Hy5™-labeled reference RNA sample were mixed pair-wise and hybridized to the miRCURY™ LNA array version 9.2 (Exiqon, Denmark), which contains capture probes targeting all miRNAs for all species registered in the miRBase version 9.2 at the Sanger Institute. The hybridization was performed according to the miRCURY™ LNA array manual using a Tecan HS4800 hybridization station (Tecan, Austria). After hybridization the microarray slides were scanned and stored in an ozone free environment (ozone level below 2.0 ppb) in order to prevent potential bleaching of the fluorescent dyes. The miRCURY™ LNA array microarray slides were scanned using the Agilent G2565BA Microarray Scanner System (Agilent Technologies, Inc., USA) and the image analysis was carried out using the ImaGene 7.0 software (BioDiscovery, Inc., USA).
Total RNA (0.5 μg) was dephosphorylated with CIP at 37°C for 30 min, samples were denatured and ligation was carried out for 2 hours at 16°C, where a molecule of Cyanine 3-pCp is incorporated to the 3′ end of RNA molecules. Labeled RNA was dried, resuspended with Hybridization Buffer and Blocking Agent, incubated 10 min at 100°C and transferred to an ice water bath for 5 min. Samples were hybridized in a volume of 45 μl to the Human miRNA V2 Oligo Microarray (Agilent) for 20 hours at 55°C with 20 rpm rotation. Post-hybridization washes were in GE Wash Buffer 1 (Agilent) at RT to remove the cover slip, followed by one wash with GE Wash Buffer 1 (Agilent) at RT for 5 min and one wash with GE Wash Buffer 2 (Agilent) at 37°C for 5 min.
Arrays were scanned on an Agilent G2565BA microarray scanner under default settings recommended by Agilent Technologies for miRNA microarrays at 100% PMT setting and 5 μm resolution. Image derived raw intensity data was extracted using Agilent Feature Extraction Software (Agilent).
Illumina miRNA sequencing
Starting from 1 μg of total RNA, small RNAs in the range of 18–30nt were separated by 15% Novex TBE-urea PAGE, excised from the gel, and eluted out of the gel slice. 5′ RNA adapters (5′-GUUCAGAGUUCUACAGUCCGACGAUC-3′) were ligated using T4 RNA ligase, and ligated fragments in the range of 40-60nt were separated and recovered as before. Thereafter, 3′ RNA adapters (5′-UCGUAUGCCGUCUUCUGCUUGUidT-3′) were ligated to the RNA. Ligation products were isolated by 10% Novex TBE-urea PAGE, recovering fragments of 70-90nt. SuperScript II reverse transcriptase was used to generate cDNA constructs with the SRA RT primer (5′-CAAGCAGAAGACGGCATACGA-3′), from the small RNA ligated with 5′ and 3′ adapters. Single stranded cDNA with adapters at both ends were selectively amplified by 15-cycle PCR reaction using Phusion DNA polymerase and primers GX1 (5′-CAAGCAGAAGACGGCATACGA-3′) and GX2 (5′-AATGATACGGCGACCACCGACAGGTTCAGAGTTCTACAGTCCGA-3′). Amplified cDNA was resolved by 6% Novex TBE PAGE and amplicon fragments of 92 nt were recovered as before. Library quality was assessed on the Agilent Technologies 2100 Bioanalyzer. DNA was loaded into a lane of a single-read flow cell at a concentration of 3–3.5 pM for cluster generation using a single-read cluster generation kit (Illumina). The sequencing primer (5′-CGACAGGTTCAGAGTTCTACAGTCCGACGATC-3′) was annealed to the clusters and the flow cell was then mounted on a Genome Analyzer (GA) I or GA II instrument for sequencing, and 36–41 sequencing cycles were performed. A PhiX control lane loaded at a concentration of 2 pM was used to monitor run quality. Image processing and base calling was performed using Illumina sequencing analysis pipelines v0.3.0 or v1.3.2. Comment: replicate 1 was sequenced in May 2008; replicates 2 and 3 were sequenced in December 2008. Some of the technicalities had changed (upgrade of GA I to GA II, different pipeline versions, different number of sequencing cycles). Raw sequencing data were further processed to specifically analyze miRNA sequences using the MIRO software (http://seq.crg.es/download/software/Miro/) (Koffler et al., unpublished). 3′ adapters were recognized and trimmed using a perl script that penalizes mismatches to a lesser extent at read ends, following the distribution of sequencing errors along Solexa reads . Where no adaptor could be recognized, sequence was cut down to 22 nt. The compiled collection of reads with removed adapters was aligned against the reduced-complexity set of miRBase version 12.0 entries (including miRBase primary miRNA, mature miRNA and star miRNAs). The mapping was performed by applying Eland iteratively in order to include all possible product sizes. Reads mapping unambiguously were counted for each unique miRBase entry within the reduced-complexity miRBase reference set. The per-miRNA count data from Illumina sequencing was normalized according to estimated effective library sizes .
Quantitative real time PCR was performed using the miRCURY LNA™ microRNA PCR System (Exiqon) on total RNA extracted from HeLa cells treated at different times with EGF (with or without protein kinase inhibitors) with miRVana’s isolation kit (Ambion) following the manufacturer’s instructions. PCR amplification and detection were performed with the ROCHE LightCycler 480 detector, using 2x SYBR GREEN Master Mix. The reaction profile was: Polymerase Activation/Denaturation cycle (95°C for 10 min) followed by 40 amplification cycles (95°C-10 sec, 60°C-20 sec). miRNA levels were calculated using the LightCycler 480 software. The data analysis was carried out using the ∆∆Ct method that provides the target gene expression values as fold changes in the problem sample compared to a calibrator sample. Both problem and calibrator samples were normalized by the relative expression of housekeeping genes (SNORD44 and SNORD48).
Analysis of miRNA variability (isomiRs)
Deep characterization of miRNA was performed using a stand-alone version of the Seqbuster software (http://estivill_lab.crg.es/seqbuster; ). Several filters were applied for miRNA variant analysis. First, the sequences considered in the analysis presented a frequency above 3. Second, 10 was chosen as the ‘Contribution Cut-Off’ parameter, meaning that every isomiR considered in the analysis contributes by more than 10% to the total number of variants annotated in the same miRNA locus. The contribution cutoff was 10% to highlight the most abundant species annotating onto a specific miRNA locus, as previously shown . However, the results did not significantly differ when using a decreased contribution cutoff. Third, we applied the Z-score option to exclude sequencing errors as the possible cause of the nucleotide changes observed in some variants .
Statistical analysis of differential miRNA expression
Log2ratio values were computed for all pairs of control and EGF stimulated samples. Analysis for differential expression on a miRNA-by-miRNA basis was performed using SAM and limma (SAM: , limma: ), including correction for multiple testing using the False Discovery Rate (FDR) method.
We also aimed at defining a consensus list of regulated genes using information from all platforms simultaneously. Since expression measures are not directly comparable between different platforms we used the Rankprod approach  that is based on log2ratio ranks. Only genes present in all the platforms under consideration can be included in each analysis. Therefore we applied the Rankprod analysis for all combinations of platforms as given by the complete merge data matrix. p-value adjustment according to  (FDR) was then applied to the union of all genes.
Network and pathway analysis
Ingenuity pathway analysis 3.1 software (IPA; Ingenuity Systems) was used for evaluating the functional significance of regulated gene profiles. Protein coding genes identified by at least 3 of 4 miRNA target prediction methods (TargetScan 4.0, PicTar, miRBase, and miRanda) as hypothetical targets of the EGF responsive miRNAs were used for network inference and pathway analyses implemented in IPA tools. The list of corresponding HUGO official gene symbols was uploaded into the IPA web tool, and each gene was mapped to the Ingenuity Pathway Knowledge Base. Significant interaction networks were generated by IPA for genes found together with higher likelihood than by random chance. Using a 99% confidence level, IPA network scores ≥2 (reflecting the negative logarithm of P values <0.01) were considered significant. Significances for biological functions were then assigned to each network by determining a P value for the enrichment of the genes in the network for such functions compared with the whole Ingenuity Pathway Knowledge Base as a reference set.
Permutation testing for enrichment in experimentally verified miRNA targets
We used the miRTarBase database  to call miRNA targets. This database collects all validated miRNA targets that have been published up to November 2012.
We used the list of de-regulated genes after EGF treatment published in our previous work . All mRNA targets from the 6 h EGF-regulated miRNAs were called from the database and crossed with the de-regulated gene list, dividing between up-regulated and down-regulated genes.
where e is the number of miRNA target genes in the de-regulated (up and down) gene list coming from the simulated data in each permutation, o is the number of miRNA target genes in the de-regulated (up and down) gene list coming from the real data, and N is the number of permutations, in this case 400.
Microarray cross-platform comparison
The fact that Agilent platform only showed 6 regulated genes suggests this platform may be less sensitive than the Exiqon one. Indeed the Exiqon platform allowed us to detect double the regulated genes than the Agilent platform (2.8% vs. 1.27%) from the total of mRNAs represented on each platform.
Illumina deep sequencing
The main limitation of array technologies is that the set of regulated genes that can be found is reduced due to the constrained number of probes included in each platform. In our experimental design, miRNAs not covered by both microarray platforms had a higher likelihood of being false positives or false negatives, with a very small regulated set of miRNAs found shared among both platforms. The only way to extend the validation without being limited by the probe content of each platform was to use of deep sequencing, which could fix the lower resolution associated to microarray measurements and should allow a cross-platform validation using different technologies. For this purpose, we used the small RNA-seq methodology developed by Illumina. We re-analyzed aliquots of total RNA from the exact same three replicate experiments that had been tested on microarrays, for the control (serum-starved) and the EGF treated at the 6 h time point.
On average, 5.5 million raw sequences were obtained per sample, which after running the MIRO analysis pipeline allowed us to monitor the expression of an average of 3.7 million unambiguous tags matching to mature miRNA with less than 3 mismatches, corresponding to 490 different miRNA genes (see Additional file 4). We found a good correlation between replicates with coefficients ranging from 0.85 to 0.96 (see Additional file 5).
Extended comparison of microarray and small RNA-seq results
Top concordant regulated miRNAs across different platforms
Upon EGF activation, multiple pathways are engaged including Ras/ERK MAPK, PI3K/AKT [1, 2]. To evaluate the role of these signaling pathways on EGF-dependent miRNA regulation, selective pharmacological inhibitors were tested for their effect on the expression changes of the four most up-regulated miRNAs in response to EGF treatment. In order to evaluate the effects of protein kinase inhibitors used in this study, as well as the effect of RasV12 transfection on ERK and AKT activation, western-blotting assays against phospho-ERK and phospho-AKT were performed (see Additional file 8). To assess the involvement of transactivation of EGFR in EGF-induced miRNAs expression, we used the specific EGFR tyrosine kinase inhibitor, AG1478. As expected, pre-treatment of HeLa cells with AG1478 abolished EGF-induced tyrosine phosphorylation of EGFR (data not shown) and totally blocked EGF-induced miR-29a, miR-29b, miR-221 and miR-222 expression (Figure 4B). The ERK inhibitor U0126 also effectively blocked overexpression of all 4 miRNAs bringing them to levels equivalent to AG1478 inhibitor treatment. In contrast, the PI3K-AKT pathway inhibitor, wortmannin, only partially blocked miRNA up-regulation provoked by EGF. Furthermore, transfection of cells with RasV12, a constitutively active Ras mutant, which activates the ERK pathway, is sufficient to promote miR-29a, miR-29b, miR-221 and miR-222 overexpression without the addition of EGF.
Inferred biological function of EGF-regulated miRNAs
Ingenuity pathways analysis
Diseases and disorders
4,26E-09 - 1,42E-02
2,72E-05 - 1,42E-02
Connective Tissue Disorders
4,26E-09 - 1,42E-02
2,22E-06 - 1,42E-02
Dermatological Diseases and Conditions
3,49E-08 - 1,42E-02
Molecular and cellular functions
Cellular Growth and Proliferation
5,25E-06 - 1,42E-02
1,28E-08 - 1,42E-02
6,24E-09 - 1,42E-02
2,80E-06 - 1,42E-02
DNA Replication, Recombination, and Repair
2,71E-06 - 1,42E-02
Sequence heterogeneity in the EGF-dependent miRNome
IsomiR expression profiles between the serum-starved and the EGF-treated cells present a high correlation for all the regulated miRNAs suggesting that EGF is not involved in the molecular mechanisms of biosynthesis leading to sequence heterogeneity of the miRNAs analyzed.
The results of profiling miRNA using genomic methods have proven to be very platform dependent. This can be explained by differences the probe design methods which at times address all possible isoforms or sometimes just target one very specific isoform.
Some of the differences observed between Agilent and Exiqon may be related to this fact. Agilent arrays use hairpin probes that require a perfect match of the 3′ end of the probe, therefore being very specific for the major reference mature isoform and less likely to detect isomiRs different from the reference isoform. On the other hand, Exiqon probe design in principle could cross-hybridize to pre-miRNA and different isomiRs from a specific class. This could explain why more miRNAs show differential expression in Exiqon. For example, miR-21 was not found to be differentially expressed on the Agilent platform, and this could be explained by the fact that the miRBase reference sequence interrogated by the hairpin loop probes is the minor isoform in HeLa cells.
Small RNA-seq is a powerful methodology that provides a qualitative improvement on the microarray technology. We have shown that it is able to detect changes similar to those detected by microarrays, as well as revealing the underlying isoform complexity that cannot be captured using microarray based assays.
In our hands, both platform approaches, high-throughput sequencing and microarrays, are complementary techniques and not mutually exclusive. We observed high agreement between both strategies for the highest ranked targets using Rankprod, which were further cross-validated using RT-qPCR.
We compared our results with those reported in a recent study describing the temporal response to EGF treatment of HeLa cells and found little overlap (Avraham et al., 2010). We believe this is due to the fact that in our study we compare the treated and untreated transcriptome at the time of interest whereas the experimental design of this other work compares to a baseline at the onset of treatment. Whereas our design would focus on the specific EGF dependent response, the comparison to the baseline may be biased toward temporal rather than treatment specific responses. In agreement with this, subtraction of control 15 minute time point in our dataset showed similar temporal profiles for most genes with or without exposure to EGF, indicating that exposure to reduced serum used in EGF exposure experiments may account for most of the responses reported previously, and therefore the response attributed to EGF in the work of Avraham may be confounded. This would explain why in our hands the number of EGF responsive miRNAs may be reduced relative to theirs.
In agreement with our observations, recent results derived from lung cancer cell lines where EGFR transcripts were targeted by shRNA expression resulted in miRNA targets highly overlapping with our own targets (including miR-221, miR-222, miR-29a,b,c) also found to be targets of MET tyrosine kinase receptor . At the functional level this would identify such targets as core to EGFR (and MET) tyrosine kinase signaling in other cell types.
EGF-dependent miRNAs in cancer
EGF treatment leads to the activation of a network of signaling pathways which leads to the up-regulation of anti-apoptotic proteins or to the inactivation of pro-apoptotic proteins. The post-transcriptional and expression events implicated in these pathways have been a matter of intense study during the last decade [5–7, 16]. A critical question about signal transduction is how a weak or transient activation achieves a robust and long-term switch in gene expression. The discovery of miRNAs opened a new door in the study of the regulation of the downstream events related to these pathways that may help disclose some of these open questions, such us the inhibitory mechanisms leading to attenuation of gene expression . This hypothesis is supported by the fact that Gene Ontology analysis indicates that several of the genes targeted by the up-regulated miRNAs after 6 hours of EGF treatment in our study are involved in the mechanisms that down-regulate the MAP kinase pathway after a mitogenic stimulus.
Recent studies have emphasized causative links between microRNA deregulation and cancer development. In fact, many oncogenes and tumor suppressor genes are regulated by miRNAs ([8, 30–33] for review) and significant down-regulation of most of the miRNAs has been reported in various tumors as compared to normal tissues .
Few studies so far have focused on the potential link between EGF signaling and miRNA [7, 17] although a direct link between EGF signaling and miRNA biosynthesis has been described. Argonaute-2 expression is regulated by epidermal growth factor receptor and mitogen-activated protein kinase signaling and correlates with a transformed phenotype in breast cancer cells  and miR-7 has been found to mediate EGF receptor signaling and promote photoreceptor differentiation in the Drosophila eye . In addition, miR-7 down-regulates EGFR mRNA and protein expression in cancer cell lines (lung, breast, and glioblastoma) and the AKT pathway [37, 38] inhibiting schwannoma cell growth  and inducing cell cycle arrest and cell death, which suggests a negative regulatory mechanism acting on the MAPK pathway . Furthermore, increasing evidence shows that the links between microRNA deregulation and cancer development may be often mediated through the MAPK pathway [41–44]. Interestingly, the 8 miRNAs described to be up-regulated in the present study had already been described to be involved in several processes such as cell growth and proliferation, cell cycle and cellular development among others, besides their implication in several types of cancers.
In vitro studies indicate that miR-21 promotes proliferation and invasion in cultured cells . miR-21 has also been found to be overexpressed under many pathophysiological conditions, especially in cancer ( for review). Indeed, miR-21 has been implicated in various aspects of carcinogenesis, including cellular proliferation, apoptosis, and migration . miR-21 is up-regulated in breast cancer  and this overexpression is correlated with specific breast cancer biopathological features such as advanced tumor stage, lymph node metastasis and poor survival of the patients. miR-21 is also up-regulated in pancreatic cancer  and in hepatocellular carcinoma (HCC) tumors . Furthermore, an aberrant expression of miR-21 can contribute to HCC growth and spreading by modulating Phosphatase and TENsin homolog (PTEN) expression which mediates phenotypic characteristics of cancer cells such as cell growth, migration and invasion .
miR-21 has been also found to target tumor suppressor Programmed Cell Death 4 (PDCD4) and stimulates invasion, intravasation and metastasis in colorectal cancer  and its inhibition in breast cancer cells causes reduced cell growth . Interestingly PDCD4, suppresses tumor progression in human colon carcinoma cells by down-regulating the MEK kinase kinase 1 (MAP4K1) gene transcription , with the consequent inhibition of c-Jun activation and Activator Protein-1 (AP-1)-dependent transcription. Given the role of PDCD4 as a negative regulator of AP-1, the miR-21-mediated down-regulation of PDCD4 seems to be essential for the maximal induction of AP-1 activity in response to RAS transformation . Interestingly we have also found that c-jun expression was down-regulated after 6 h of EGF treatment . Recent studies have implicated the down-regulation of the tumor suppressor gene acidic (leucine-rich) nuclear phosphoprotein 32 family, member A (ANP32A)  and Spry2  of miR-21 to its oncogenic function.
miR-29 regulates the anti-apoptotic Bcl-2 family member, Mcl-1  and DNA methyltransferases (DNMT) 3A and 3B  and is down-regulated in lung cancer where DNMT is frequently up-regulated. In addition, several studies have implicated miR-29 in several types cancers such as acute myeloid leukemia , chronic lymphocytic leukemia , hepatocellular carcinoma  and cutaneous melanoma .
miR-132 is elevated in pancreatic cancer  and in the monocytic leukemia cell line THP-1 by Lipopolysaccharide (LPS) treatment. LPS triggers activation of NF-κB and AP-1 transcription factors and results in up-regulation of immune-responsive genes. REST (Re1 silencing transcription factor) and its cofactor complex are targets of miR-132 . Given the role of both REST and miRNA as repressors, this suggests a double-negative feedback loop between REST and the miR-132 for stabilizing and maintaining neuronal gene expression. miR-132 is also induced by photic entrainment cues via a MAPK/CREB-dependent mechanism, modulates clock-gene expression, and attenuates the entraining effects of light .
miR-31 is a pleiotropic miRNA up-regulated in colorectal tumors and its expression levels are correlated with the stage of the tumor . On the contrary, it has been found to be down-regulated greater than two-fold in breast cancer . miR-31 also inhibits breast cancer metastasis  and over-expression of miR-31 significantly enhanced proliferation and tumorigenicity of lung cancer cells . An interesting recently published study indicates that restoring miR-31 function in already-established metastases induces metastatic regression and opens the possibility of intervention strategies centered in restoring miR-31 function in order to combat metastasis .
Distinct up-regulations of miR-222, miR-221, and miR-31 have been observed in HCC . miR-221 and miR-222 are potent regulators of proliferation, mediated by regulation of p27Kip, a cell cycle inhibitor and tumor suppressor . High levels of miR-221 and miR-222 appear in human thyroid papillary carcinomas , prostate cancer  and in glioblastomas [36, 73] and correlate with low levels of p27kip. Indeed, sustained activation of ERK1/2 by NGF induces miR-221 and miR-222 in PC12 cells .
Finally, miR-215 has been described to be a p53-responsive miRNA  able to induce cell cycle arrest. miR-215 suppresses Denticleless protein homolog expression, inducing a decreased cell proliferation by causing G2-arrest  and contributing to enhanced CDKN1A/p21 levels, colony suppression, cell cycle arrest, and cell detachment from a solid support .
Blocking of miR-29a, miR-29b, miR-221 and miR-222 expression by the selective MEK inhibitor U0126 in an analogous manner to AG1478 and the overexpression of these miRNAs in RasV12 transfected cells suggest that the MAPK pathway is a critical effector downstream of EGFR for the transmission of the information leading to miRNAs regulation. The same effect but to a lesser extent is observed after PI3K/AKT pathway inhibition with wortmannin. Indeed, both MAPK and PI3K/AKT pathways are activated by EGF-treatment controlling cell fates such as cell proliferation, survival and apoptosis [78, 79]. Expression of most of the miRNAs regulated by EGF found in the present work induces a decrease (miR-7, miR-215) or increase in cell proliferation (miR-21, miR-221, miR-222, miR-31).
Further studies are required to identify their in vivo targets and to elucidate the specific regulatory networks where they are involved.
The results obtained from isomiR analysis show the presence of sequence heterogeneity for all the miRNA analyzed, supporting the idea of a novel degree of complexity for the miRNA transcriptome [67, 80] as a result of inexact Dicer processing. Furthermore, some of the miRNA species studied have a majority of mature miRNA sequences that differ from that listed on databases. Lee et al. catalogued miRNA end heterogeneity from different biological samples and found that from 34% to 51% of the detectable miRNA species have a predominant sequence differing from that listed in databases . However, the equal proportion of the different types of variants in serum-starved versus EGF-treated indicates that the molecular mechanisms involved in the generation of these variants are not regulated by EGF. Several studies already confirmed the presence of these variants in multiple and diverse kinds of samples, tissues and species, however, up to this date, there is no data regarding the functional role of isomiR diversity as well as their generation mechanisms besides a recent study suggesting that 5′ end variations may produce specificity changes in the type of AGOs with which these miRNA isoforms preferentially associate . In addition a biological role for isomiRs is supported by the finding of regulation of variability in isomiR expression .
In the present study we have used microarrays and deep sequencing analysis to present the EGF-induced miRNome and its integration into functional cellular networks.
Knowing that there are biases in genomic studies that are platform dependent, in our study we have attempted to get around this limitation to increase the confidence in the detected miRNome changes. By applying a procedure for cross-platform data tests we are able to present a dataset, also validated by RT-qPCR, which grants more reliable analyses at the functional genomics level and more robust predictions of networks of co-regulated microRNA target genes.
Specifically, we find that miR-29a, miR-29b, miR-221, miR-222, miR-21, miR-132, miR-215 and miR-7 are all up-regulated 6 hours after EGF treatment in HeLa cells, with an inferred effect on genes with a function in the regulation of cell proliferation, cell cycle, and DNA repair as well as gene expression. This supports a prominent role for these miRNA in modulating the effects of EGF on cancer-related cellular processes.
In addition, we make a comparison of performance between microarray and high throughput sequencing profiling that suggests the two methodologies combined may survey the miRNome in a better way than each on its own.
Finally, we observed the presence of sequence heterogeneity for all the EGF-dependent microRNAs, supporting the idea of a novel degree of complexity for the miRNA transcriptome. There is no apparent effect of EGF on isomiR regulation but we can confirm cell type specific differences in the prevalent edited isoforms relative to the miRBase reference miRNAs.
We thank other CRG microarray laboratory members for advice and discussions. This work was supported by grants from Ministerio de Ciencia e Innovación [MICINN; BFU2012-32617] to JADR and FL, Generalitat de Catalunya [SGR2009-366] to JADR and FL, and the EU-FP7 PRIORITY to JADR. Start up funds from the Institute for Predictive and Personalized Medicine of Cancer and the Center for Genomic Regulation to LS; by the Spanish Ministry of Science and Technology [SAF2004-06976] to LS, Juan de la Cierva researcher contract to FL, and by excellence in research team recognitions by the Catalan government, Departament de Innovació, Universitats i Ensenyament, Generalitat de Catalunya [SGR2005-404] to LS, and by the Instituto de Salud Carlos III Fondo de Investigaciones Sanitarias to JADR and to LS [PI10/01154].
- Burgess AW, Cho HS, Eigenbrot C, Ferguson KM, Garrett TP, Leahy DJ, Lemmon MA, Sliwkowski MX, Ward CW, Yokoyama S: An open-and-shut case? Recent insights into the activation of EGF/ErbB receptors. Mol Cell. 2003, 12 (3): 541-552. 10.1016/S1097-2765(03)00350-2.View ArticlePubMedGoogle Scholar
- Henson ES, Gibson SB: Surviving cell death through epidermal growth factor (EGF) signal transduction pathways: implications for cancer therapy. Cell Signal. 2006, 18 (12): 2089-2097. 10.1016/j.cellsig.2006.05.015.View ArticlePubMedGoogle Scholar
- Olayioye MA, Neve RM, Lane HA, Hynes NE: The ErbB signaling network: receptor heterodimerization in development and cancer. EMBO J. 2000, 19 (13): 3159-3167. 10.1093/emboj/19.13.3159.PubMed CentralView ArticlePubMedGoogle Scholar
- Bianco R, Gelardi T, Damiano V, Ciardiello F, Tortora G: Rational bases for the development of EGFR inhibitors for cancer treatment. Int J Biochem Cell Biol. 2007, 39 (7–8): 1416-1431.View ArticlePubMedGoogle Scholar
- Olsen JV, Blagoev B, Gnad F, Macek B, Kumar C, Mortensen P, Mann M: Global, in vivo, and site-specific phosphorylation dynamics in signaling networks. Cell. 2006, 127 (3): 635-648. 10.1016/j.cell.2006.09.026.View ArticlePubMedGoogle Scholar
- Amit I, Citri A, Shay T, Lu Y, Katz M, Zhang F, Tarcic G, Siwak D, Lahad J, Jacob-Hirsch J, Amariglio N, Vaisman N, Segal E, Rechavi G, Alon U, Mills GB, Domany E, Yarden Y: A module of negative feedback regulators defines growth factor signaling. Nat Genet. 2007, 39 (4): 503-512. 10.1038/ng1987.View ArticlePubMedGoogle Scholar
- Avraham R, Sas-Chen A, Manor O, Steinfeld I, Shalgi R, Tarcic G, Bossel N, Zeisel A, Amit I, Zwang Y, Enerly E, Russnes HG, Biagioni F, Mottolese M, Strano S, Blandino G, Borresen-Dale AL, Pilpel Y, Yakhini Z, Segal E, Yarden Y: EGF decreases the abundance of microRNAs that restrain oncogenic transcription factors. Sci Signal. 2010, 3 (124): ra43-10.1126/scisignal.2000876.View ArticlePubMedGoogle Scholar
- Singh SK, Pal Bhadra M, Girschick HJ, Bhadra U: MicroRNAs–micro in size but macro in function. FEBS J. 2008, 275 (20): 4929-4944. 10.1111/j.1742-4658.2008.06624.x.View ArticlePubMedGoogle Scholar
- Chang TC, Mendell JT: microRNAs in vertebrate physiology and human disease. Annu Rev Genomics Hum Genet. 2007, 8: 215-239. 10.1146/annurev.genom.8.080706.092351.View ArticlePubMedGoogle Scholar
- Hammell M, Long D, Zhang L, Lee A, Carmack CS, Han M, Ding Y, Ambros V: mirWIP: microRNA target prediction based on microRNA-containing ribonucleoprotein-enriched transcripts. Nat Methods. 2008, 5 (9): 813-819. 10.1038/nmeth.1247.PubMed CentralView ArticlePubMedGoogle Scholar
- Blenkiron C, Goldstein LD, Thorne NP, Spiteri I, Chin SF, Dunning MJ, Barbosa-Morais NL, Teschendorff AE, Green AR, Ellis IO, Tavare S, Caldas C, Miska EA: MicroRNA expression profiling of human breast cancer identifies new markers of tumor subtype. Genome Biol. 2007, 8 (10): R214-10.1186/gb-2007-8-10-r214.PubMed CentralView ArticlePubMedGoogle Scholar
- Davalos V, Esteller M: MicroRNAs and cancer epigenetics: a macrorevolution. Curr Opin Oncol. 2010, 22 (1): 35-45. 10.1097/CCO.0b013e328333dcbb.View ArticlePubMedGoogle Scholar
- Cherqui D, Laurent A, Mocellin N, Tayar C, Luciani A, Van Nhieu JT, Decaens T, Hurtova M, Memeo R, Mallat A, Duvoux C: Liver resection for transplantable hepatocellular carcinoma: long-term survival and role of secondary liver transplantation. Ann Surg. 2009, 250 (5): 738-746. 10.1097/SLA.0b013e3181bd582b.View ArticlePubMedGoogle Scholar
- Deng S, Calin GA, Croce CM, Coukos G, Zhang L: Mechanisms of microRNA deregulation in human cancer. Cell Cycle. 2008, 7 (17): 2643-2646. 10.4161/cc.7.17.6597.View ArticlePubMedGoogle Scholar
- Inui M, Martello G, Piccolo S: MicroRNA control of signal transduction. Nat Rev Mol Cell Biol. 2010, 11 (4): 252-263.View ArticlePubMedGoogle Scholar
- Llorens F, Hummel M, Pastor X, Ferrer A, Pluvinet R, Vivancos A, Castillo E, Iraola S, Mosquera AM, Gonzalez E, Lozano J, Ingham M, Dohm JC, Noguera M, Kofler R, del Rio JA, Bayes M, Himmelbauer H, Sumoy L: Multiple platform assessment of the EGF dependent transcriptome by microarray and deep tag sequencing analysis. BMC Genomics. 2011, 12: 326-10.1186/1471-2164-12-326.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayashi T, Koyama N, Gresik EW, Kashimata M: Detection of EGF-dependent microRNAs of the fetal mouse submandibular gland at embryonic day 13. J Med Invest. 2009, 56 (Suppl): 250-252.View ArticlePubMedGoogle Scholar
- Garofalo M, Romano G, Di Leva G, Nuovo G, Jeon Y-J, Ngankeu A, Sun J, Lovat F, Alder H, Condorelli G, Engelman JA, Ono M, Rho JK, Cascione L, Volinia S, Nephew KP, Croce CM: EGFR and MET receptor tyrosine kinase-altered microRNA expression induces tumorigenesis and gefitinib resistance in lung cancers. Nat Med. 2012, 18 (1): 74-82.Google Scholar
- Dohm JC, Lottaz C, Borodina T, Himmelbauer H: Substantial biases in ultra-short read data sets from high-throughput DNA sequencing. Nucleic Acids Res. 2008, 36 (16): e105-10.1093/nar/gkn425.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson MD, Oshlack A: A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010, 11 (3): R25-10.1186/gb-2010-11-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Pantano L, Estivill X, Marti E: SeqBuster, a bioinformatic tool for the processing and analysis of small RNAs datasets, reveals ubiquitous miRNA modifications in human embryonic cells. Nucleic Acids Res. 2010, 38 (5): e34-10.1093/nar/gkp1127.PubMed CentralView ArticlePubMedGoogle Scholar
- Marti E, Pantano L, Banez-Coronel M, Llorens F, Minones-Moyano E, Porta S, Sumoy L, Ferrer I, Estivill X: A myriad of miRNA variants in control and Huntington's disease brain regions detected by massively parallel sequencing. Nucleic Acids Res. 2010, 38 (20): 7219-7235. 10.1093/nar/gkq575.PubMed CentralView ArticlePubMedGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.PubMed CentralView ArticlePubMedGoogle Scholar
- Wettenhall JM, Smyth GK: limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics. 2004, 20 (18): 3705-3706. 10.1093/bioinformatics/bth449.View ArticlePubMedGoogle Scholar
- Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J: RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics. 2006, 22 (22): 2825-2827. 10.1093/bioinformatics/btl476.View ArticlePubMedGoogle Scholar
- Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125 (1–2): 279-284.View ArticlePubMedGoogle Scholar
- Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, Tsai WT, Chen GZ, Lee CJ, Chiu CM, Chien CH, Wu MC, Huang CY, Tsou AP, Huang HD: miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011, 39 (Database issue): D163-D169.PubMed CentralView ArticlePubMedGoogle Scholar
- Raymond CK, Roberts BS, Garrett-Engele P, Lim LP, Johnson JM: Simple, quantitative primer-extension PCR assay for direct monitoring of microRNAs and short-interfering RNAs. RNA. 2005, 11 (11): 1737-1744. 10.1261/rna.2148705.PubMed CentralView ArticlePubMedGoogle Scholar
- Wieduwilt MJ, Moasser MM: The epidermal growth factor receptor family: biology driving targeted therapeutics. Cell Mol Life Sci. 2008, 65 (10): 1566-1584. 10.1007/s00018-008-7440-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Skaftnesmo KO, Prestegarden L, Micklem DR, Lorens JB: MicroRNAs in tumorigenesis. Curr Pharm Biotechnol. 2007, 8 (6): 320-325. 10.2174/138920107783018390.View ArticlePubMedGoogle Scholar
- Zimmerman AL, Wu S: MicroRNAs, cancer and cancer stem cells. Cancer Lett. 2011, 300 (1): 10-19. 10.1016/j.canlet.2010.09.019.PubMed CentralView ArticlePubMedGoogle Scholar
- Fabbri M, Calin GA: Epigenetics and miRNAs in human cancer. Adv Genet. 2010, 70: 87-99.View ArticlePubMedGoogle Scholar
- Garzon R, Marcucci G, Croce CM: Targeting microRNAs in cancer: rationale, strategies and challenges. Nat Rev Drug Discov. 2010, 9 (10): 775-789. 10.1038/nrd3179.PubMed CentralView ArticlePubMedGoogle Scholar
- Hwang HW, Mendell JT: MicroRNAs in cell proliferation, cell death, and tumorigenesis. Br J Cancer. 2006, 94 (6): 776-780. 10.1038/sj.bjc.6603023.PubMed CentralView ArticlePubMedGoogle Scholar
- Adams BD, Claffey KP, White BA: Argonaute-2 expression is regulated by epidermal growth factor receptor and mitogen-activated protein kinase signaling and correlates with a transformed phenotype in breast cancer cells. Endocrinology. 2009, 150 (1): 14-23.PubMed CentralView ArticlePubMedGoogle Scholar
- Gillies JK, Lorimer IA: Regulation of p27Kip1 by miRNA 221/222 in glioblastoma. Cell Cycle. 2007, 6 (16): 2005-2009. 10.4161/cc.6.16.4526.View ArticlePubMedGoogle Scholar
- Kefas B, Godlewski J, Comeau L, Li Y, Abounader R, Hawkinson M, Lee J, Fine H, Chiocca EA, Lawler S, Purow B: microRNA-7 inhibits the epidermal growth factor receptor and the Akt pathway and is down-regulated in glioblastoma. Cancer Res. 2008, 68 (10): 3566-3572. 10.1158/0008-5472.CAN-07-6639.View ArticlePubMedGoogle Scholar
- Giles KM, Barker A, Zhang PM, Epis MR, Leedman PJ: MicroRNA regulation of growth factor receptor signaling in human cancer cells. Methods Mol Biol. 2011, 676: 147-163. 10.1007/978-1-60761-863-8_11.View ArticlePubMedGoogle Scholar
- Saydam O, Senol O, Wurdinger T, Mizrak A, Ozdener GB, Stemmer-Rachamimov AO, Yi M, Stephens RM, Krichevsky AM, Saydam N, Brenner GJ, Breakefield XO: miRNA-7 attenuation in Schwannoma tumors stimulates growth by upregulating three oncogenic signaling pathways. Cancer Res. 2011, 71 (3): 852-861. 10.1158/0008-5472.CAN-10-1219.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (9): 5731-5741.View ArticlePubMedGoogle Scholar
- Roccaro AM, Sacco A, Chen C, Runnels J, Leleu X, Azab F, Azab AK, Jia X, Ngo HT, Melhem MR, Burwick N, Varticovski L, Novina CD, Rollins BJ, Anderson KC, Ghobrial IM: microRNA expression in the biology, prognosis, and therapy of Waldenstrom macroglobulinemia. Blood. 2009, 113 (18): 4391-4402. 10.1182/blood-2008-09-178228.PubMed CentralView ArticlePubMedGoogle Scholar
- Dangi-Garimella S, Yun J, Eves EM, Newman M, Erkeland SJ, Hammond SM, Minn AJ, Rosner MR: Raf kinase inhibitory protein suppresses a metastasis signalling cascade involving LIN28 and let-7. EMBO J. 2009, 28 (4): 347-358. 10.1038/emboj.2008.294.PubMed CentralView ArticlePubMedGoogle Scholar
- Rahadiani N, Takakuwa T, Tresnasari K, Morii E, Aozasa K: Latent membrane protein-1 of Epstein-Barr virus induces the expression of B-cell integration cluster, a precursor form of microRNA-155, in B lymphoma cell lines. Biochem Biophys Res Commun. 2008, 377 (2): 579-583. 10.1016/j.bbrc.2008.10.007.View ArticlePubMedGoogle Scholar
- Chou YT, Lin HH, Lien YC, Wang YH, Hong CF, Kao YR, Lin SC, Chang YC, Lin SY, Chen SJ, Chen HC, Yeh SD, Wu CW: EGFR promotes lung tumorigenesis by activating miR-7 through a Ras/ERK/Myc pathway that targets the Ets2 transcriptional repressor ERF. Cancer Res. 2010, 70 (21): 8822-8831. 10.1158/0008-5472.CAN-10-0638.View ArticlePubMedGoogle Scholar
- Zheng J, Xue H, Wang T, Jiang Y, Liu B, Li J, Liu Y, Wang W, Zhang B, Sun M: miR-21 downregulates the tumor suppressor P12(CDK2AP1) and stimulates cell proliferation and invasion. J Cell Biochem. 2010, 112 (3): 872-80.View ArticleGoogle Scholar
- Jazbutyte V, Thum T: MicroRNA-21: from cancer to cardiovascular disease. Curr Drug Targets. 2010, 11 (8): 926-935. 10.2174/138945010791591403.View ArticlePubMedGoogle Scholar
- Papagiannakopoulos T, Shapiro A, Kosik KS: MicroRNA-21 targets a network of key tumor-suppressive pathways in glioblastoma cells. Cancer Res. 2008, 68 (19): 8164-8172. 10.1158/0008-5472.CAN-08-1305.View ArticlePubMedGoogle Scholar
- Yan LX, Huang XF, Shao Q, Huang MY, Deng L, Wu QL, Zeng YX, Shao JY: MicroRNA miR-21 overexpression in human breast cancer is associated with advanced clinical stage, lymph node metastasis and patient poor prognosis. RNA. 2008, 14 (11): 2348-2360. 10.1261/rna.1034808.PubMed CentralView ArticlePubMedGoogle Scholar
- Dillhoff M, Liu J, Frankel W, Croce C, Bloomston M: MicroRNA-21 is overexpressed in pancreatic cancer and a potential predictor of survival. J Gastrointest Surg. 2008, 12 (12): 2171-2176. 10.1007/s11605-008-0584-x.PubMed CentralView ArticlePubMedGoogle Scholar
- Connolly E, Melegari M, Landgraf P, Tchaikovskaya T, Tennant BC, Slagle BL, Rogler LE, Zavolan M, Tuschl T, Rogler CE: Elevated expression of the miR-17-92 polycistron and miR-21 in hepadnavirus-associated hepatocellular carcinoma contributes to the malignant phenotype. Am J Pathol. 2008, 173 (3): 856-864. 10.2353/ajpath.2008.080096.PubMed CentralView ArticlePubMedGoogle Scholar
- Meng F, Henson R, Wehbe-Janek H, Ghoshal K, Jacob ST, Patel T: MicroRNA-21 regulates expression of the PTEN tumor suppressor gene in human hepatocellular cancer. Gastroenterology. 2007, 133 (2): 647-658. 10.1053/j.gastro.2007.05.022.PubMed CentralView ArticlePubMedGoogle Scholar
- Asangani IA, Rasheed SA, Nikolova DA, Leupold JH, Colburn NH, Post S, Allgayer H: MicroRNA-21 (miR-21) post-transcriptionally downregulates tumor suppressor Pdcd4 and stimulates invasion, intravasation and metastasis in colorectal cancer. Oncogene. 2008, 27 (15): 2128-2136. 10.1038/sj.onc.1210856.View ArticlePubMedGoogle Scholar
- Tang Z, Bauer JA, Morrison B, Lindner DJ: Nitrosylcobalamin promotes cell death via S nitrosylation of Apo2L/TRAIL receptor DR4. Mol Cell Biol. 2006, 26 (15): 5588-5594. 10.1128/MCB.00199-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Talotta F, Cimmino A, Matarazzo MR, Casalino L, De Vita G, D'Esposito M, Di Lauro R, Verde P: An autoregulatory loop mediated by miR-21 and PDCD4 controls the AP-1 activity in RAS transformation. Oncogene. 2009, 28 (1): 73-84. 10.1038/onc.2008.370.View ArticlePubMedGoogle Scholar
- Schramedei K, Morbt N, Pfeifer G, Lauter J, Rosolowski M, Tomm JM, von Bergen M, Horn F, Brocke-Heidrich K: MicroRNA-21 targets tumor suppressor genes ANP32A and SMARCA4. Oncogene. 2011, 30 (26): 2975-85. 10.1038/onc.2011.15.PubMed CentralView ArticlePubMedGoogle Scholar
- Kwak HJ, Kim YJ, Chun KR, Woo YM, Park SJ, Jeong JA, Jo SH, Kim TH, Min HS, Chae JS, Choi EJ, Kim G, Shin SH, Gwak HS, Kim SK, Hong EK, Lee GK, Choi KH, Kim JH, Yoo H, Park JB, Lee SH: Downregulation of Spry2 by miR-21 triggers malignancy in human gliomas. Oncogene. 2011, 30 (21): 2433-42. 10.1038/onc.2010.620.View ArticlePubMedGoogle Scholar
- Mott JL, Kobayashi S, Bronk SF, Gores GJ: mir-29 regulates Mcl-1 protein expression and apoptosis. Oncogene. 2007, 26 (42): 6133-6140. 10.1038/sj.onc.1210436.PubMed CentralView ArticlePubMedGoogle Scholar
- Fabbri M, Garzon R, Cimmino A, Liu Z, Zanesi N, Callegari E, Liu S, Alder H, Costinean S, Fernandez-Cymering C, Volinia S, Guler G, Morrison CD, Chan KK, Marcucci G, Calin GA, Huebner K, Croce CM: MicroRNA-29 family reverts aberrant methylation in lung cancer by targeting DNA methyltransferases 3A and 3B. Proc Natl Acad Sci U S A. 2007, 104 (40): 15805-15810. 10.1073/pnas.0707628104.PubMed CentralView ArticlePubMedGoogle Scholar
- Garzon R, Heaphy CE, Havelange V, Fabbri M, Volinia S, Tsao T, Zanesi N, Kornblau SM, Marcucci G, Calin GA, Andreeff M, Croce CM: MicroRNA 29b functions in acute myeloid leukemia. Blood. 2009, 114 (26): 5331-5341. 10.1182/blood-2009-03-211938.PubMed CentralView ArticlePubMedGoogle Scholar
- Santanam U, Zanesi N, Efanov A, Costinean S, Palamarchuk A, Hagan JP, Volinia S, Alder H, Rassenti L, Kipps T, Croce CM, Pekarsky Y: Chronic lymphocytic leukemia modeled in mouse by targeted miR-29 expression. Proc Natl Acad Sci U S A. 2010, 107 (27): 12210-12215. 10.1073/pnas.1007186107.PubMed CentralView ArticlePubMedGoogle Scholar
- Xiong Y, Fang JH, Yun JP, Yang J, Zhang Y, Jia WH, Zhuang SM: Effects of microRNA-29 on apoptosis, tumorigenicity, and prognosis of hepatocellular carcinoma. Hepatology. 2010, 51 (3): 836-845.PubMedGoogle Scholar
- Nguyen T, Kuo C, Nicholl MB, Sim MS, Turner RR, Morton DL, Hoon DS: Downregulation of microRNA-29c is associated with hypermethylation of tumor-related genes and disease outcome in cutaneous melanoma. Epigenetics. 2011, 6 (3): 388-394. 10.4161/epi.6.3.14056.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu J, Xie X: Comparative sequence analysis reveals an intricate network among REST, CREB and miRNA in mediating neuronal gene expression. Genome Biol. 2006, 7 (9): R85-10.1186/gb-2006-7-9-r85.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheng HY, Papp JW, Varlamova O, Dziema H, Russell B, Curfman JP, Nakazawa T, Shimizu K, Okamura H, Impey S, Obrietan K: microRNA modulation of circadian-clock period and entrainment. Neuron. 2007, 54 (5): 813-829. 10.1016/j.neuron.2007.05.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Bandres E, Cubedo E, Agirre X, Malumbres R, Zarate R, Ramirez N, Abajo A, Navarro A, Moreno I, Monzo M, Garcia-Foncillas J: Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues. Mol Cancer. 2006, 5: 29-PubMed CentralView ArticlePubMedGoogle Scholar
- Valastyan S, Reinhardt F, Benaich N, Calogrias D, Szasz AM, Wang ZC, Brock JE, Richardson AL, Weinberg RA: A pleiotropically acting microRNA, miR-31, inhibits breast cancer metastasis. Cell. 2009, 137 (6): 1032-1046. 10.1016/j.cell.2009.03.047.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee LW, Zhang S, Etheridge A, Ma L, Martin D, Galas D, Wang K: Complexity of the microRNA repertoire revealed by next-generation sequencing. RNA. 2010, 16 (11): 2170-2180. 10.1261/rna.2225110.PubMed CentralView ArticlePubMedGoogle Scholar
- Valastyan S, Chang A, Benaich N, Reinhardt F, Weinberg RA: Activation of miR-31 function in already-established metastases elicits metastatic regression. Genes Dev. 2011, 25 (6): 646-659. 10.1101/gad.2004211.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong QW, Lung RW, Law PT, Lai PB, Chan KY, To KF, Wong N: MicroRNA-223 is commonly repressed in hepatocellular carcinoma and potentiates expression of Stathmin1. Gastroenterology. 2008, 135 (1): 257-269. 10.1053/j.gastro.2008.04.003.View ArticlePubMedGoogle Scholar
- Miller TE, Ghoshal K, Ramaswamy B, Roy S, Datta J, Shapiro CL, Jacob S, Majumder S: MicroRNA-221/222 confers tamoxifen resistance in breast cancer by targeting p27Kip1. J Biol Chem. 2008, 283 (44): 29897-29903. 10.1074/jbc.M804612200.PubMed CentralView ArticlePubMedGoogle Scholar
- Visone R, Russo L, Pallante P, De Martino I, Ferraro A, Leone V, Borbone E, Petrocca F, Alder H, Croce CM, Fusco A: MicroRNAs (miR)-221 and miR-222, both overexpressed in human thyroid papillary carcinomas, regulate p27Kip1 protein levels and cell cycle. Endocr Relat Cancer. 2007, 14 (3): 791-798. 10.1677/ERC-07-0129.View ArticlePubMedGoogle Scholar
- Mercatelli N, Coppola V, Bonci D, Miele F, Costantini A, Guadagnoli M, Bonanno E, Muto G, Frajese GV, De Maria R, Spagnoli LG, Farace MG, Ciafre SA: The inhibition of the highly expressed miR-221 and miR-222 impairs the growth of prostate carcinoma xenografts in mice. PLoS One. 2008, 3 (12): e4029-10.1371/journal.pone.0004029.PubMed CentralView ArticlePubMedGoogle Scholar
- le Sage C, Nagel R, Egan DA, Schrier M, Mesman E, Mangiola A, Anile C, Maira G, Mercatelli N, Ciafre SA, Farace MG, Agami R: Regulation of the p27(Kip1) tumor suppressor by miR-221 and miR-222 promotes cancer cell proliferation. EMBO J. 2007, 26 (15): 3699-3708. 10.1038/sj.emboj.7601790.PubMed CentralView ArticlePubMedGoogle Scholar
- Terasawa K, Ichimura A, Sato F, Shimizu K, Tsujimoto G: Sustained activation of ERK1/2 by NGF induces microRNA-221 and 222 in PC12 cells. FEBS J. 2009, 276 (12): 3269-3276. 10.1111/j.1742-4658.2009.07041.x.View ArticlePubMedGoogle Scholar
- Braun CJ, Zhang X, Savelyeva I, Wolff S, Moll UM, Schepeler T, Orntoft TF, Andersen CL, Dobbelstein M: p53-Responsive micrornas 192 and 215 are capable of inducing cell cycle arrest. Cancer Res. 2008, 68 (24): 10094-10104. 10.1158/0008-5472.CAN-08-1569.PubMed CentralView ArticlePubMedGoogle Scholar
- Song B, Wang Y, Titmus MA, Botchkina G, Formentini A, Kornmann M, Ju J: Molecular mechanism of chemoresistance by miR-215 in osteosarcoma and colon cancer cells. Mol Cancer. 2010, 9: 96-10.1186/1476-4598-9-96.PubMed CentralView ArticlePubMedGoogle Scholar
- Georges SA, Biery MC, Kim SY, Schelter JM, Guo J, Chang AN, Jackson AL, Carleton MO, Linsley PS, Cleary MA, Chau BN: Coordinated regulation of cell cycle transcripts by p53-Inducible microRNAs, miR-192 and miR-215. Cancer Res. 2008, 68 (24): 10105-10112. 10.1158/0008-5472.CAN-08-1846.View ArticlePubMedGoogle Scholar
- Takeuchi K, Ito F: EGF receptor in relation to tumor development: molecular basis of responsiveness of cancer cells to EGFR-targeting tyrosine kinase inhibitors. FEBS J. 2010, 277 (2): 316-326. 10.1111/j.1742-4658.2009.07450.x.View ArticlePubMedGoogle Scholar
- Sanui A, Yotsumoto F, Tsujioka H, Fukami T, Horiuchi S, Shirota K, Yoshizato T, Kawarabayashi T, Kuroki M, Miyamoto S: HB-EGF inhibition in combination with various anticancer agents enhances its antitumor effects in gastric cancer. Anticancer Res. 2010, 30 (8): 3143-3149.PubMedGoogle Scholar
- Ebhardt HA, Tsang HH, Dai DC, Liu Y, Bostan B, Fahlman RP: Meta-analysis of small RNA-sequencing errors reveals ubiquitous post-transcriptional RNA modifications. Nucleic Acids Res. 2009, 37 (8): 2461-2470. 10.1093/nar/gkp093.PubMed CentralView ArticlePubMedGoogle Scholar
- Ebhardt HA, Fedynak A, Fahlman RP: Naturally occurring variations in sequence length creates microRNA isoforms that differ in argonaute effector complex specificity. Silence. 2010, 1 (1): 12-10.1186/1758-907X-1-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Fernandez-Valverde SL, Taft RJ, Mattick JS: Dynamic isomiR regulation in Drosophila development. RNA. 2010, 16 (10): 1881-1888. 10.1261/rna.2379610.PubMed CentralView ArticlePubMedGoogle Scholar