Identification of microRNA-mRNA functional interactions in UVB-induced senescence of human diploid fibroblasts
© Greussing et al.; licensee BioMed Central Ltd. 2013
Received: 29 March 2012
Accepted: 8 March 2013
Published: 4 April 2013
Cellular senescence can be induced by a variety of extrinsic stimuli, and sustained exposure to sunlight is a key factor in photoaging of the skin. Accordingly, irradiation of skin fibroblasts by UVB light triggers cellular senescence, which is thought to contribute to extrinsic skin aging, although molecular mechanisms are incompletely understood. Here, we addressed molecular mechanisms underlying UVB induced senescence of human diploid fibroblasts.
We observed a parallel activation of the p53/p21WAF1 and p16INK4a/pRb pathways. Using genome-wide transcriptome analysis, we identified a transcriptional signature of UVB-induced senescence that was conserved in three independent strains of human diploid fibroblasts (HDF) from skin. In parallel, a comprehensive screen for microRNAs regulated during UVB-induced senescence was performed which identified five microRNAs that are significantly regulated during the process. Bioinformatic analysis of miRNA-mRNA networks was performed to identify new functional mRNA targets with high confidence for miR-15a, miR-20a, miR-20b, miR-93, and miR-101. Already known targets of these miRNAs were identified in each case, validating the approach. Several new targets were identified for all of these miRNAs, with the potential to provide new insight in the process of UVB-induced senescence at a genome-wide level. Subsequent analysis was focused on miR-101 and its putative target gene Ezh2. We confirmed that Ezh2 is regulated by miR-101 in human fibroblasts, and found that both overexpression of miR-101 and downregulation of Ezh2 independently induce senescence in the absence of UVB irradiation. However, the downregulation of miR-101 was not sufficient to block the phenotype of UVB-induced senescence, suggesting that other UVB-induced processes induce the senescence response in a pathway redundant with upregulation of miR-101.
We performed a comprehensive screen for UVB-regulated microRNAs in human diploid fibroblasts, and identified a network of miRNA-mRNA interactions mediating UVB-induced senescence. In addition, miR-101 and Ezh2 were identified as key players in UVB-induced senescence of HDF.
KeywordsAging Stress-induced senescence miR-101 Ezh2 Epigenetics UVB
Cellular senescence is controlled through a variety of regulatory mechanisms with particular contribution of the p53  and Rb  tumor suppressor pathways. Briefly, p53 is activated in response to telomere damage or other kinds of DNA damage and orchestrates the transcriptional activation of a variety of downstream genes contributing to growth arrest, such as the CDK inhibitor p21WAF1. In a parallel pathway, upregulation of the CDK inhibitor p16INK4a leads to the inhibition of phosphorylation of the retinoblastoma protein (pRb), thereby enforcing prolonged cell cycle arrest , which is also a useful marker for senescence in human tissues . Activation of p16INK4a triggers the appearance of senescence-associated heterochromatin foci , at least in some cell types. Activation of the p53/p21WAF1 and p16INK4a/pRb pathways represent two alternative scenarios for senescence initiation, and there is evidence that parallel activation of both pathways leads to enforced senescence , also referred to as “deep senescence” . The relative importance of each of these pathways for cellular senescence seems to depend on the cell type. It should be noted that, besides these central pathways, a variety of signals have been shown to initiate a senescence response .
Stress-induced premature senescence (SIPS) has been identified as a model for telomere-independent senescence that can be initiated by various forms of stress, including oxidative stress, irradiation, replicative stress and oncogene activation . In the case of skin derived human diploid fibroblasts (HDF), repeated mild treatment with UVB has been established as an experimental model for extrinsic skin aging , which depends on the accumulation of senescent cells, in particular fibroblasts in the dermis . UVB irradiation is known to induce aspects of a DNA damage response, and it was reported that knocking down p53 alleviates but does not abrogate the senescence response of HDF to repeated UVB stress , suggesting that both p53-dependent and -independent pathways may cooperate to enforce UVB-induced senescence. However, UVB irradiation is known to affect a variety of intracellular signal transduction pathways [12, 13] and the precise sequence of events during establishment of UVB-induced premature senescence remained elusive.
Recent studies have emphasized an important role of non-coding RNAs, also referred to as microRNAs, as regulators of gene expression . MicroRNAs are generated from larger precursor RNAs and were shown to interfere with the expression of protein coding genes by several mechanisms, including i) destabilization of specific mRNAs and ii) prevention of translation of specific mRNAs. Messenger RNAs targeted by microRNAs usually contain short sequences of homology . It is known that changes in microRNA expression are contributing to cellular senescence  and organismic aging [17, 18]; however, the role of microRNAs, if any, in UVB-induced senescence of human fibroblasts is only poorly understood. In the present communication, we have addressed molecular mechanisms underlying the establishment of premature senescence in human fibroblasts exposed to repeated series of mild UVB irradiations, as a model system to monitor molecular processes involved in extrinsic skin aging. We identified a small set of microRNAs, which are differentially regulated during the process and performed transcriptome analysis, in combination with advanced bioinformatics, to identify potential targets for these microRNAs.
UVB-induced changes in gene expression
Human diploid fibroblasts (HDF, strain HFF-2) derived from newborn foreskin were subjected to eight consecutive UVB treatments of 4000 J/m2 during four days. Under these conditions no overt cell death was observed (data not shown). Cell proliferation was strongly inhibited, and UVB-treated cells performed less than 2 population doublings (PDL) over 15 days of the experiment, whereas untreated cells underwent 12 PDL in the same time period. The growth arrest phenotype reached by repeated mild UVB stress resembled cellular senescence, as judged by the changes of cell morphology and the accumulation of a large percentage of cells that stained positive for senescence-associated β-galactosidase (SA-ß-gal) (Additional file 1: Figure S1). As could be expected, UVB treatment led to a strong but transient phosphorylation of p53 on serine 15, most notable at day 4, i.e. after application of the last stress, indicative of p53 activation, consistent with previous observations . The overall levels of p53 protein were also increased, in line with previous findings . Activation of p53 also resulted in the upregulation of its downstream effector p21WAF1, which is known to enforce cellular senescence in response to DNA damage . We also addressed effects of UVB treatment on the p16INK4A/pRb pathway, representing the other important pathway relevant for cellular senescence . As could be expected from the observed upregulation of p21WAF1, UVB treatment led to the rapid disappearance of phosphorylated species of the retinoblastoma protein (pRb), which was accompanied by a delayed but significant upregulation of the CDK inhibitor p16INK4A (Additional file 1: Figure S1). Together the data suggest that mild UVB stress induces premature cellular senescence and that both the p53/p21WAF1 axis and the p16INK4A/pRb pathway are involved in the senescence response, consistent with previous reports [23–25]. The same protocol of repeated mild UVB stress induced premature senescence also in two other strains of human diploid fibroblasts, HFF-1 and PFF, with a very similar kinetics (Additional file 2: Figure S2), along with the activation of the p53/p21WAF1 and pRb/p16INK4A pathways (Additional file 3: Figure S3).
microRNAs regulated in UVB-induced senescence
miRNA-mRNA regulatory networks in UVB-induced senescence
miR-15a, along with miR-16, is commonly deleted in human chronic lymphocytic leukemia  and known to target multiple oncogenes, including BCL2, MCL1, CCND1, and WNT3A . Our analysis confirmed a high confidence interaction between miR-15a and vEGF-A in UVB-induced senescence (Figure 5B), thereby validating the analytical procedure. Results of the bioinformatic analysis suggested PAPPA, APOD, RRAS2, Runx1, RARB, BTG2, Notch3 and SFRP1 as potential targets for miR-15 (Figure 5B).
To further validate the bioinformatics-based target selection, regulation of several newly identified candidate genes by specific microRNAs was addressed in HDF overexpressing miR-15a, miR-20a, and miR-93, respectively. In all cases, microRNA levels were significantly ( > 20 fold) increased by transfection (data not shown). For the genes RARB, RUNX1 and CDKN2B, we found that overexpression of the appropriate microRNA species (miR-15a and miR-20a, respectively) reduced protein levels of the respective gene products (Additional file 6: Figure S4). Whereas these findings do not prove binding of microRNAs 15a and 20a to the 3'-UTR of the target genes, downregulation of protein expression by overexpression of selected miRNAs provides a functional validation of the bioinformatics approach. The mRNAs for several other genes, including VEGFA, RHOJ, and NOTCH3 were significantly down-regulated by the appropriate microRNAs (data not shown); however, the limited availability of high-quality antibodies precluded determination of protein expression levels in these cases.
A role of miR-101/Ezh2 in UVB-induced senescence?
In this communication, we addressed molecular mechanisms underlying UVB-induced senescence of human diploid fibroblasts. Using genome-wide transcriptome analysis, we identified a transcriptional signature of UVB-induced senescence. In parallel, a comprehensive screen for microRNAs regulated during UVB-induced senescence was performed and five microRNAs were identified that are significantly regulated during the process. Subsequent analysis revealed several well established miRNA-mRNA regulatory interactions including miR-101/Ezh2, thereby validating the assay. In addition, several new miRNA-mRNA regulatory interactions were identified to occur in UVB-induced senescence. Overexpression of miR-101 and downregulation of Ezh2 independently induced senescence in the absence of UVB irradiation. However, the downregulation of miR-101 was not sufficient to block the phenotype of UVB-induced senescence, suggesting that other UVB-induced processes contribute to the senescence response, partially redundant with upregulation of miR-101. These findings extend our knowledge of miR-mRNA regulatory interactions, and at the same time provide a rich resource for data mining to identify new players in UVB-induced senescence and potential targets for interventions in extrinsic skin aging.
Transcriptional profiling of UVB-induced premature senescence
Whereas a low-density DNA array was used previously to study changes in the level of 240 senescence-related genes in UVB-induced senescence of HDF , the current study provides the first genome-wide transcriptional analysis of UVB-induced senescence. As could be expected , activation of p53 signaling and suppression of the G1/S transition were observed along with a pronounced DNA damage response, as revealed by pathway analysis tools. We also observed concerted regulation of genes in the aryl hydrocarbon receptor signaling and hepatic fibrosis/heaptic stellate activation pathways, which conceivably reflect specific responses to UVB irradiation, used here to trigger the senescence response. Of note, there were also distinct changes in the regulation of genes coding for enzymes in the pyrimidine and nicotinamide metabolism; however, more work will be required to fully understand the implications of these alterations for UVB-induced senescence. Conceivably, the transcriptional signature of UVB-induced senescence, as shown here, will pinpoint new targets for intervention in extrinsic photoaging, which will be investigated in follow-up studies.
miRNA-mRNA regulatory networks in UVB-induced senescence
miR-34, a known transcriptional target of p53 , was strongly upregulated at day 7, as expected from the observed activation of p53. Similarly, downregulation of miR-20a and miR-20b is consistent with previous observations suggesting that members of the miR-17-92 cluster are commonly downregulated in various senescence models as well as in organismal aging . These observations validated the miRNA screening approach. It is known that miRNAs regulate both stability and translation of mRNAs, and in most cases upregulation of miRNAs leads to the inhibition of gene function . The miR-17-92 cluster, containing the microRNA-17-18-19-20-92 polycistron, cooperates with Myc in tumorigenesis , probably by silencing of p21WAF1 expression at the post-transcriptional level . On the other hand, miR-17-92 inhibited proliferation and metastasis of pancreatic carcinoma cells , blocked E2F-dependent steps in the regulation of angiogenesis , and repressed endothelial cell migration , suggesting cell type specific effects. Our finding that E2F1 and VEGFA are relevant targets for miR-20a/b in the context of UVB-induced senescence is consistent with these data.
Although frequently overexpressed in human malignancies, miR-93 may actually function as a tumor suppressor gene. Accordingly, miR-93 abrogated VEGF protein secretion, suggesting that miR-93 interferes with angiogenesis , blocked tumor development in mammary fat pads , and suppressed proliferation of human colon cancer stem cells . However, miR-93 was also shown to i) promote tumor growth and angiogenesis by targeting integrin-β8 , ii) to contribute to silencing of p21WAF1 gene expression after DNA damage , and iii) to promote cell proliferation and clonogenicity of HepG2 Cells , suggesting that effects of miR-93 depend on cellular context. Upregulation of p21WAF1 gene expression in UVB-treated cells was correlated with the downregulation of miRNAs 20a/b and 93, known to target p21WAF1[27, 30, 44–46], and the identification of E2F1 and VEGFA as relevant targets for miR-93 in the context of UVB-induced senescence is consistent with the available data.
miR-15a, along with miR-16, was the first microRNA linked to cancer because both genes are commonly deleted in human chronic lymphocytic leukemia . Expression of miRNAs encoded by the miR-15/16 cluster inhibits cell proliferation, promotes apoptosis of cancer cells, and suppresses tumorigenicity both in vitro and in vivo. miR-15a and miR-16-1 function by targeting multiple oncogenes, including BCL2, MCL1, CCND1, and WNT3A [33, 47]. In contrast to their function as tumor suppressors, miR-15 can also promote tumor growth and progression, when expressed in cancer-associated fibroblasts . The identification of VEGFA as relevant target for miR-15 in the context of UVB-induced senescence is consistent with these data. On the other hand, mRNA levels of E2F1 and Bcl2 were downregulated in UVB-treated cells (Additional file 4: Table S1 and data not shown), although three miRNAs known to target E2F1 (miR 20a, 20b and 93) and one miRNA known to target Bcl2 (miR-15a) were all significantly downregulated in UVB-treated cells (Figure 3). It is conceivable that, in these cases, additional regulatory processes take place, which remain to be identified. A particular complication in the interpretation of the current results lies in pleiotropic effects of UVB irradiation, which are known to affect the expression of selected miRNAs and mRNAs (also shown here) but also may influence protein translation  and/or protein stability . Accordingly, the assignment of miRNAs and gene expression changes reflects plausibility rather than a strict mechanistic dependence. Bioinformatics analysis and preliminary validation analysis performed in the current study highlighted several miRNA targets as potential modulators of UVB-induced fibroblast senescence that were not previously described, including the cdk inhibitor CDKN2B and the transcription factors RUNX1 and RARB. Whereas these findings establish CDKN2B, RUNX1 and RARB as functional target genes for miR-15a and miR-20A, respectively, the binding of these microRNAs to the 3'-UTR of the target genes remains to be confirmed by additional experiments.
Regulation of UVB-induced senescence: the role of miR-101 and Ezh2
The here reported data demonstrate for the first time the implication of miR-101/Ezh2 signaling in UVB-induced senescence of human dermal fibroblasts. In control experiments, upregulation of miR-101 and the concomitant downregulation of Ezh2 was also observed in two independent HDF strains, HFF-1 and PFF (Additional file 3: Figure S3), indicating that the regulation of these components is a conserved feature in the senescence response of human dermal fibroblasts to UVB irradiation. The inverse relationship between miR-101 and Ezh2 expression levels was noticed before [51, 52] and Ezh2 has been shown by others to play a role in cellular senescence [53, 54]. When this communication was under revision, it was reported by others that miR-101 controlled Ezh2 function in cellular senescence of mouse embryonic fibroblasts . In the present communication, functional interactions between miR-101 and Ezh2 in UVB-induced senescence of HDF were analyzed in more detail. For example, we tried by knockdown of miR-101 to rescue the expression levels of Ezh2 in UVB irradiated cells. Although knockdown of miR-101 was very efficient also in UVB treated cells, no corresponding increase in the level of Ezh2 mRNA or protein was observed, suggesting that expression of the Ezh2 gene is regulated by additional signals, which remain to be established. It is conceivable that p53 is responsible for the effect, since p53 was shown to suppress the Ezh2 gene promoter [56, 57]. Conversely, we also tried to rescue the phenotype of UVB-induced senescence by overexpression of Ezh2 from a lentiviral overexpression vector. Whereas Ezh2 mRNA was significantly upregulated in infected cells, this did not lead to any detectable increase of Ezh2 protein levels, suggesting that overexpression of Ezh2 protein is not well tolerated, at least in HDF. The mechanisms underlying downregulation of Ezh2 in this cell type remain to be understood. After UVB irradiation, Ezh2 protein levels were rapidly reduced, even prior to the induction of miR-101 (Figure 3), strongly suggesting that several different UVB-dependent pathways converge to downregulate Ezh2, the upregulation of miR-101 being just one of several triggers for this process.
We report here a comprehensive screen for microRNAs and mRNAs regulated during UVB-induced senescence in human diploid fibroblasts. Using advanced bioinformatics solutions, we identified a network of miRNA-mRNA interactions mediating UVB-induced senescence in this cell type, providing a rich resource for future data mining. The data reported in this communication illustrate the regulation of five distinct miRNAs during UVB-induced cellular senescence. Together the results obtained in this study suggest important roles for microRNAs miR-15, miR-20a/b, miR-93 and miR-101, and their mRNA targets, during UVB-induced senescence of human diploid fibroblasts.
All chemicals were purchased from Sigma, unless indicated otherwise.
Human diploid foreskin fibroblasts (HDF) were either purchased from ATCC (Manassas, VA) (HFF-1 #SCRC-1041; HFF-2 #SCRC-1042) or isolated from newborn foreskin (PFF) in our laboratory, as described . Cells were used at passage 6 (HFF-2, PFF) and passage 10 (HFF-1). HDFs and human osteosarcoma cells (U2-OS; obtained from ATCC, Manassas, VA) were cultured in the same way in DMEM (Sigma) as described . The cumulative population doublings (cPDL) were calculated using the following equation: cPDL = (log(A) – log(B))/0.301 (A: number of cells at the end of one passage; B: number of cells that were seeded at the beginning of one passage).
Senescence-associated-β-galactosidase (SA-ß-gal) staining was used to determine the senescent status of the cells. To stain for SA-β-gal, cells were grown on 6-well plates and washed three times with PBS. Afterwards, the cells were fixed with 2% formaldehyde and 0.4% glutaraldehyde in PBS for 5 minutes at room temperature. Cells were washed three times with PBS and prepared for staining as described previously . Therefore, cells were covered with staining solution (150 mM NaCl, 2 mM MgCl, 5 mM potassium ferricyanide, 5 mM potassium ferrocyanide, 40 mM citric acid, 12 mM sodium phosphate, pH 6.0, adding 1 mg/mL 5-bromo-4-chloro-3-indolyl-b-D-galactoside [X-gal] directly before use) and incubated for 24 h at 37°C without light exposition. The reaction was stopped by washing off the staining solution with PBS. Cells were covered with PBS and blue staining indicating the presence of SA-β-gal can be detected under the microscope. To calculate the percentage of SA-β-gal positive cells, stained cells were counted and related to the total cell number.
For UVB treatment fibroblasts were seeded out in 10 cm dishes (Greiner Bio One, Austria) at a density of 6 × 105 (UVB) and 3 × 105 (control). Cells were washed with HBSS (Sigma) and covered with 2 ml HBSS. To calculate the irradiation time, power per area [W/m2] was measured by a UVX radiometer (Thermo Fisher) and the following equation was used: irradiation time [s] = energy per area [J/m2] divided by power per area [W/m2]. The output of a Philips TL20W/01 lamp (Philips, The Netherlands) for 10 cm dishes was determined as 14.2 ± 0.5 W/m2. To test the sublethal dosage, cells were irradiated with 3000, 3500, 4000 and 5000 J/m2. For experiments, cells were irradiated twice a day with a dose of 4000 J/m2 for 4 consecutive days, where irradiation time was 282 ± 10 s. UVB treatment of miRNA knockdown cells were performed in 6-well plates at a density of 1.3 × 105 cells (UVB) and 5 × 104 cells (control). The output of a Philips TL20W/01 lamp for 6-well plates was determined as 17.3 ± 0.5 W/m2 and irradiation time was 231 ± 6 s. After the first (day 1) and the last (day 4) irradiation, as well as after day 7, day 9 and day 15 after the first irradiation, cells were lysed and RNA and protein were isolated.
For the preparation of whole cell lysates, HDFs were washed twice with cold PBS and scraped off on ice in lysis buffer (50 mM Tris–HCl, 150 mM NaCl, 1% NP-40, 0.25% Na-deoxycholate, 1 mM EDTA, 100 nM Na3VO4, 1 mM NaF, 10 mM β-glycerophosphate, pH 7.4) from the 6-well plate or 10 cm dish with a rubber policeman. Cells were three times deep-frozen in liquid nitrogen and thawed and further kept on ice for 30 minutes. After centrifugation at 20,000 × g for 10 minutes at 4°C, supernatant was used to determine the protein concentration by DC Protein Assay Kit (Biorad, Austria).
mRNA and miRNA isolation
Total RNA was isolated using either TRIzol® Reagent (Invitrogen; for miRNA) or the RNeasy® Mini Kit (Qiagen; RNA) according the manufacturer’s protocol. For the isolation of RNA with TRIzol® Reagent cells were trysinized, centrifuged, washed with PBS and re-suspended in 1 ml TRIzol® Reagent. After 5 minutes incubation at room temperature 200 μl of Chloroform was added and vigorously vortexed. After 3 minutes incubation at room temperature a centrifugation step was performed for 10 minutes at 14 000 × g at 4°C. The aqueous phase was transferred to a micro-centrifuge tube and mixed with 500 μl isopropanol and incubated at room temperature for 10 minutes to precipitate the RNA. After 10 minutes centrifugation at 14 000 × g (4°C) the supernatant was discarded and the pellet was washed with ethanol (70%) and centrifuged again for 5 minutes at 6 000 × g (4°C). The supernatant was discarded, the pellet air dried and re-suspended in 30 μl of RNAse free water. Using the RNeasy® Mini Kit up to 5 × 106 cells were trypsinized and lysed in 350 μl RLT buffer (including 10 μl β-mercaptoethanol per ml). The lysate was homogenously mixed by pipetting up and down a few times. One volume of 70% ethanol was added, mixed by pipetting and transferred to RNeasy spin column. After a centrifugation step at 12 900 × g for 30 seconds at room temperature, three washing steps were performed. First 700 μl RW1 buffer, two times 500 μl RPE (one volume RPE added to four volumes ethanol) whereas the last washing step was carried out for 2 minutes at maximum speed to dry the RNA-binding membrane. RNA was eluted using 30 μl of RNAse-free water. RNA concentration was quantified by photometric measurement at 260 nm and 280 nm.
Equal amounts of protein were subjected to SDS gel electrophoresis (10–12.5% SDS/polyacrylamide gel) and transferred to PVDF membrane by wet electro-blotting (300 mA, 1 h) using the standard Western blot protocol. The membranes were blocked with 5% skim milk in PBS-T (Phosphate buffered saline + 0.1% Tween 20) or in 5% BSA in PBS-T for the detection of phosphorylated proteins for 1 h at room temperature. Primary antibodies were incubated for 1 h at room temperature or overnight at 4°C. After two times of washing with PBS-T, the second antibody was incubated for 45 min at room temperature. After a few washing steps with PBS-T, immune-reactive proteins were detected using an enhanced chemiluminescence system (ECL+, Amersham Life Science, Germany). The following antibodies were used: mouse monoclonal anti-p21WAF1 (Pharmingen, #556430) mouse monoclonal anti-pRb (Pharmingen, #554136), mouse monoclonal anti-p53 (Santa Cruz, #sc-126), rabbit polyclonal anti-GAPDH (Santa Cruz, #sc-25778), rabbit polyclonal phospho-p53 (Ser15; Cell Signaling, #9284), mouse monoclonal α-Tubulin (Sigma, #t-5168), mouse monoclonal ant-Cyclin D1 (Neomarkers, #MS-210-P) and mouse monoclonal anti-Ezh2 (BD Biosciences, #612666), mouse monoclonal anti-p16INK4A (BD Biosciences, #511325GR), rabbit polyclonal anti-p15INK4b (Abcam, #ab53034), mouse monoclonal anti-RUNX1/AML1 (Abcam, #ab54869), rabbit monoclonal anti-retinoic acid receptor beta – RARB (abcam, #ab124701). As secondary antibodies, polyclonal antibodies from Dako were used. For positive control, U2-OS cells were transiently transfected with a pcDNA3.1 vector carrying the sequence of either Ezh2 or p16INK4a.
Real-time q-PCR analysis of mRNAs
RNA was isolated using the RNeasy® Mini Kit (Qiagen) and quantified as described above. For cDNA synthesis 0.4-1 μg of total RNA was reverse transcribed with the Transcriptor First Strand cDNA Synthesis Kit purchased from Roche Applied Science and diluted 1:4. Amplification of Ezh2 mRNA was carried out with the following primers: (forward) 5’-CAT TCG GTA AAT CCA AAC TGC-3’ and (reverse) 5’-CGA CAT ACT TCA GGG CAT CA-3’. RUNX1 (Runt-related transcription factor 1) (forward) 5’-CCC TCG TGC CTC CCT GAA CCA-3’ ; (reverse) 5’-GGC TGG GGA GAG GGA TGG ACA-3’. VEGFA (Vascular endothelial growth factor A) (forward) 5’-TTT GCT TGC CAT TCC CCA C-3’; (reverse) 5’-GCT CTT GCT ACC TCT TTC CTC-3’. CDKN2B (Cyclin-dependent kinase 4 inhibitor B) (forward) 5’-GAT GAG GAC AAT GAG GCA AAG-3’ ; (reverse) 5’-TGG GAA GAA AAG CAA GAC AAC-3’. NOTCH3 (Neurogenic locus notch homolog protein 3) (forward) 5’-TCA TCC TCT TCT CTT TCC ACC-3’; (reverse) 5’-TCC CAG ACT CTT CAC AAG AC-3’. As a housekeeper B2M (β-2 microglobulin) was amplified with the following primers: (forward) 5’-GAA TTC ACC CCC ACT GAA AA-3’ and (reverse) 5’-CTC CAT GAT GCT GCT TAC A-3’ and GAPDH (Glyceraldehyde 3-phosphate dehydrogenase) (forward) 5’-GAG TCA ACG GAT TTG GTC GT-3’ ; (reverse) 5’-GAT CTC GCT CCT GGA AGA TG-3’. Real-time q-PCR was performed in duplicates using the SYBR Green I Master Mix in the LightCycler® 480 Instrument (Roche Applied Science).
Real-time q-PCR analysis of microRNAs
RNA was extracted with TRIzol® Reagent (Invitrogen, as described above). Using the Universal cDNA Synthesis Kit (Exiqon, Denmark) cDNA synthesis was performed according the manufacturer’s protocol. To perform real-time PCR SYBR Green master mix and specific primer for miR-101 (hsa-miR-101, LNA PCR primer set, UniRT; Exiqon, Denmark) were used. Expression levels of miR-101 were normalized to endogenous 5S rRNA (Exiqon, Denmark). In addition, 7 other microRNAs (hsa-miR-155, hsa-miR-15a, hsa-miR-17, hsa-miR-20a, hsa-miR-20b, hsa-miR-34a, hsa-miR-93) were chosen for qPCR confirmation of array data using the Taqman qPCR platform (Life Technologies). In brief, specific reverse transcription (RT) reactions were performed with each microRNA primer and RNU44 as endogenous control using 10 ng of total RNA as input material. Following RT, qPCRs were run in 4 replicates using the Taqman Universal Mastermix (Life Technologies) on a Rotor-Gene Q (Qiagen) according to the manufacturer’s protocol. Data analysis was performed using the ddCt method.
Genome-wide RNA profiling (microarray)
Total RNA was isolated using the RNeasy® Mini Kit (Qiagen, as described above). After quantification RNA was sent to The Microarray Facility (Tübingen, Germany). The obtained data set was analyzed using CARMAweb 1.5 (https://carmaweb.genome.tugraz.at/carma/). A cut-off of 1.5 fold change was used to determine the total number of up- and downregulated genes.
Low density array (TaqMan q-PCR)
RT-q-PCR was performed using Taqman® Low Density Array technology (Applied Biosystems). A total of 93 candidate and three housekeeper genes were used to design Taqman® custom array. According to the manufacturer’s protocol q-PCR was performed, as described .
Locked nucleic acid (LNA)-miRNA microarray
For miRNA expression profiling, LNA-miRNA microarrays (Exiqon, Denmark) consisting of 559 human and 170 mouse LNA-modified probes  against miRNAs annotated in Sanger miRBase v9.2  as well as 77 not yet annotated probes (miRPlus, Exiqon, Denmark) were used as previously described . In brief, total RNA extracts were end-labeled using Cy3 dye and hybridized against a common reference RNA-pool end-labeled with Cy5 on a Tecan HS 400 hybridization station (Tecan, Austria). Arrays were scanned at 10 μM resolution (Axon Genepix 4000B, Axon Instruments) and raw intensities were extracted using GenePixPro 4.1 software (Axon Instruments). The acquired array signal data was further processed under R/Bioconductor using linear models for microarray analysis  and differential expression between control and UVB treated cells was calculated using t-statistics and p-value adjustment to multiple testing according to Benjamini-Hochberg.
Stable overexpression and knockdown of Ezh2
For the overexpression of Ezh2 the lentiviral pLenti6/V5-DEST Gateway vector (Invitrogen) was used. Cloning included the TOPO cloning of Ezh2 into pENTR/D-TOPO. This vector was used to introduce the Ezh2 coding sequence into pLenti6/V5-DEST by recombination to generate the transfer vector pLenti6-Ezh2 (for further details see Invitrogen’s ViraPower Lentiviral Expression System manual). The GFP-control vector was generated in the same way. As a transfer vector for knockdown of Ezh2 lentiviral pLKO.1-TRC short-hairpin vector were purchased from Addgene/Open Biosystems (United Kindom). The following sequences were chosen: pLKO-73 within the 3’UTR of Ezh2 (5’-TAT TGC CTT CTC ACC AGC TGC-3’) and pLKO-75 within the coding region (5’-CCA ACA CAA GTC ATC CCA TTA-3’). As a control the empty vector was used. For packaging of the lentivirus in HEK293FT cells (Invitrogen), 3 μg of the corresponding overexpression or knockdown vector was combined with 7.5 μg vector backbone psPAX2, 2.5 μg envelope encoding plasmid pMD2.G. HEK293FT cells were cultivated in T75 flasks to 90% confluence and transfected with the mixture using Lipofectamine 2000 (Invitrogen). The next day medium was replaced by 10 ml of growth medium without antibiotics. After 48 hours (day 5) the supernatant was harvested, centrifuged at 300 × g for 5 minutes at room temperature and filtered through a 0.45 μm PVDF filter (Millipore, Vienna, Austria). Afterwards, the virus containing supernatant was concentrated using Polyethylene glycol. Therefore, to one volume Polyethylene glycol solution (50 mM Polyethylene glycol, 41 mM NaCl, autoclave, pH =7.2; PEG) four volumes of supernatant was added and incubated for 2 hours at 4°C, carefully mixing every 20–30 minutes by inverting. The solution was centrifuged at 1500 × g for 30 minutes at 4°C and a white pellet should be visible. After aspirating the supernatant an additional centrifugation step was carried out at 1500 × g for 5 minutes at 4°C to collect remaining PEG which is aspirated carefully. The white pellet was re-suspended in medium by pipetting up and down and vigorously vortexing for 20 to 30 seconds. As a guideline 500 μl for one T75 flask was used, aliquoted to 100 μl and stored at −80°C. The titer of the concentrated lentiviral supernatant was determined by seeding 5 × 104 U2-OS in 6-well plates with 8 mg/mL hexadimethrine bromide (polybrene; Sigma-Aldrich) as transduction enhancer. A stock of 4 mg/mL of polybrene was prepared by resolving polybrene in sterile water and filtered through a 0.22 μm sterile filter. The following day the medium was replaced by 2 ml of DMEM. After day 3 the selection was started. For overexpression selection 10 μg/ml Blasticidin, for knockdown 500 ng/ml puromycin was added to each well. The antibiotic containing medium was replaced every second day. Approximately 6–7 days after starting the selection the untransduced cells are dead. For the staining the cells are washed three times with PBS, covered with crystal violet and incubated for 5–10 minutes at room temperature. Wells are washed twice with ddH2O and air dried. The colonies are counted and multiplied with 1 000 and the corresponding dilution. This procedure allows to calculate transfection units (TU) of the virus/ml.
Overexpression and knockdown of miRNAs
To achieve the overexpression of microRNA in HDFs, cells were reverse transfected with Pre-miR™ miRNA Precursor for miR-15a, miR-20a, miR-93, miR-101, and Pre-miR™ miRNA Precursor Molecules-Negative Control #2 for negative control (Applied Biosystems, Austria) using siPORT™ NeoFX™ Transfection Agent (Ambion, Austria) according the manufacturer’s protocol. Cells were trysinized, counted and 6 × 104 cells were centrifuged to be re-suspended in 1.8 ml medium. The transfection reagent siPORT™ NeoFX™ and Opti-MEM medium (Invitrogen) were adjusted to room temperature. Opti-MEM medium was mixed with Precursors and with the transfection reagent and incubated for 10 minutes incubation at room temperature. Both reagents are mixed and additionally incubated at room temperature for 10 minutes. Afterwards, the mixture is added together with the cell suspension to 6-well plates . Cells were UVB treated as described above. Experiments were done in 6-well plates with 1.3 × 105 (UVB) and 5 × 104. After the last exposure, cells were trypsinized and prepared for transfection. For miR-101 knockdown, cells were reverse transfected with miRCURY LNA™ microRNA Power Inhibitor and miRCURY LNA™ microRNA Power Inhibitor Negative Control A for negative control (Exiqon, Denmark) using siPORT™ NeoFX™ Transfection Agent (Ambion, Austria) according the manufacturer’s protocol. UVB treated cells were trypsinized, counted, centrifuged and re-suspended in 1.8 ml medium. Opti-MEM medium was mixed with the miRNA inhibitor or control and with the transfection reagent and incubated for 10 minutes at room temperature. The protocol was continued as described above. For both, Power Inhibitor and Negative Control A 15 nM (day 4 and day 7) were used.
Transduction of cells by lentiviral particles
For transduction 5 × 104 human diploid fibroblasts (HDFs) are seeded to 6-well plates the day before (day 1). A multiplicity of infection of two together with 8 μg/ml polybrene as transduction enhancer was used in a total of one milliliter. The medium was changed the next day and selection was started at day 3. For the selection of Ezh2 overexpressing cells 10 μg/ml blasticidin, for Ezh2 knockdown cells 500 ng/ml puromycin was used.
Bioinformatics analysis for predicting miRNA target genes
For miRNA target prediction the following tools were used: TargetScan , PicTar , miRanda , PITA , ElMMo , RNA22 , DIANA-microT  and GenMiR++. The predictions of these tools reflect miRNA:mRNA pairing, site location, conservation, site accessibility, multiple sites and expression profile . The results were visualized using in house developed tools (Genesis  and ClueGO ) as well as publicly available tools: Cytoscape .
We used limma package  for miRNA microarray analysis. T-test and adjusted p-value were used to identify genes differentially expressed between control and UVB treated cells. Furthermore, we also investigated the linear relationship between miRNAs and their high confidence target genes by using Pearson correlation coefficients. All analyses were done with the use of statistical software programs R/Bioconductor.
Pathway analysis and network visualization
The data set obtained from the genome-wide RNA profiling was uploaded to Ingenuity Systems (http://www.ingenuity.com/) and pathway analysis was done. The visualization of the network was done using the open-source Cytoscape software platform  for visualizing biomolecular interaction networks. The color and shape of nodes are based on different node attributes available for the analyzed dataset. The red triangles, purple circles and orange diamonds in the network are indicating miRNAs, target genes, and transcription factors, respectively. Edges represent correlation between miRNAs and mRNAs, the color of the edges designate the type of interaction. Red is for positive and green is for negative correlation. The biological functions triggered by the miRNA through its target genes were visualized using ClueGO  a Cytoscape plugin.
Enhancer of zeste homolog 2
Hanks’ balanced salt slution
Human diploid fibroblasts
Mitogen-activated protein kinase
Real-time quantitative PCR
Small hairpin RNA
Ultraviolet light band B.
We acknowledge excellent technical support by Michael Neuhaus and Hans-Peter Viertler. The work was supported by grants from the Austrian Science Funds (NFN S93 and SFB021). Genome Research Austria GEN-AU (Project 820982 “Non-coding RNAs” and “Bioinformatics Integration Network”); and grants by the Herzfelder’sche Familienstiftung and CE.R.I.E.S to JG, and a BOKU Doc scholarship to MH, and Lise-Meitner Scholarship to RM.
- Vigneron A, Vousden KH: p53, ROS and senescence in the control of aging. Aging (Albany NY). 2010, 2 (8): 471-474.Google Scholar
- Thomas DM, Yang HS, Alexander K, Hinds PW: Role of the retinoblastoma protein in differentiation and senescence. Cancer Biol Ther. 2003, 2 (2): 124-130.View ArticlePubMedGoogle Scholar
- Lukas J, Parry D, Aagaard L, Mann DJ, Bartkova J, Strauss M, Peters G, Bartek J: Retinoblastoma-protein-dependent cell-cycle inhibition by the tumour suppressor p16. Nature. 1995, 375 (6531): 503-506. 10.1038/375503a0.View ArticlePubMedGoogle Scholar
- Ressler S, Bartkova J, Niederegger H, Bartek J, Scharffetter-Kochanek K, Jansen-Durr P, Wlaschek M: p16INK4A is a robust in vivo biomarker of cellular aging in human skin. Aging Cell. 2006, 5 (5): 379-389. 10.1111/j.1474-9726.2006.00231.x.View ArticlePubMedGoogle Scholar
- Narita M, Nunez S, Heard E, Lin AW, Hearn SA, Spector DL, Hannon GJ, Lowe SW: Rb-mediated heterochromatin formation and silencing of E2F target genes during cellular senescence. Cell. 2003, 113 (6): 703-716. 10.1016/S0092-8674(03)00401-X.View ArticlePubMedGoogle Scholar
- Beausejour CM, Krtolica A, Galimi F, Narita M, Lowe SW, Yaswen P, Campisi J: Reversal of human cellular senescence: roles of the p53 and p16 pathways. EMBO J. 2003, 22 (16): 4212-4222. 10.1093/emboj/cdg417.PubMed CentralView ArticlePubMedGoogle Scholar
- Passos JF, Nelson G, Wang C, Richter T, Simillion C, Proctor CJ, Miwa S, Olijslagers S, Hallinan J, Wipat A, et al: Feedback between p21 and reactive oxygen production is necessary for cell senescence. Mol Syst Biol. 2010, 6: 347-PubMed CentralView ArticlePubMedGoogle Scholar
- Campisi J: Cellular senescence: putting the paradoxes in perspective. Curr Opin Genet Dev. 2011, 21 (1): 107-112. 10.1016/j.gde.2010.10.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Toussaint O, Remacle J, Dierick JF, Pascal T, Frippiat C, Zdanov S, Magalhaes JP, Royer V, Chainiaux F: From the hayflick mosaic to the mosaics of ageing. Role of stress-induced premature senescence in human ageing. Int J Biochem Cell Biol. 2002, 34 (11): 1415-1429. 10.1016/S1357-2725(02)00034-1.View ArticlePubMedGoogle Scholar
- Chainiaux F, Magalhaes JP, Eliaers F, Remacle J, Toussaint O: UVB-induced premature senescence of human diploid skin fibroblasts. Int J Biochem Cell Biol. 2002, 34 (11): 1331-1339. 10.1016/S1357-2725(02)00022-5.View ArticlePubMedGoogle Scholar
- Borlon C, Vankoningsloo S, Godard P, Debacq-Chainiaux F, Toussaint O: Identification of p53-dependent genes potentially involved in UVB-mediated premature senescence of human skin fibroblasts using siRNA technology. Mech Ageing Dev. 2008, 129 (3): 109-119. 10.1016/j.mad.2007.10.014.View ArticlePubMedGoogle Scholar
- Takashima A: UVB-dependent modulation of epidermal cytokine network: roles in UVB-induced depletion of langerhans cells and dendritic epidermal T cells. J Dermatol. 1995, 22 (11): 876-887.View ArticlePubMedGoogle Scholar
- el-Ghorr AA, Norval M: Biological effects of narrow-band (311 nm TL01) UVB irradiation: a review. J Photochem Photobiol B. 1997, 38 (2–3): 99-106.View ArticlePubMedGoogle Scholar
- Almeida MI, Reis RM, Calin GA: MicroRNA history: discovery, recent applications, and next frontiers. Mutat Res. 2011, 717 (1–2): 1-8.View ArticlePubMedGoogle Scholar
- Thomas M, Lieberman J, Lal A: Desperately seeking microRNA targets. Nat Struct Mol Biol. 2010, 17 (10): 1169-1174. 10.1038/nsmb.1921.View ArticlePubMedGoogle Scholar
- Hackl M, Brunner S, Fortschegger K, Schreiner C, Micutkova L, Muck C, Laschober GT, Lepperdinger G, Sampson N, Berger P, et al: miR-17, miR-19b, miR-20a, and miR-106a are down-regulated in human aging. Aging Cell. 2010, 9 (2): 291-296. 10.1111/j.1474-9726.2010.00549.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Liang R, Bates DJ, Wang E: Epigenetic control of MicroRNA expression and aging. Curr Genomics. 2009, 10 (3): 184-193. 10.2174/138920209788185225.PubMed CentralView ArticlePubMedGoogle Scholar
- Grillari J, Grillari-Voglauer R: Novel modulators of senescence, aging, and longevity: Small non-coding RNAs enter the stage. Exp Gerontol. 2010, 45 (4): 302-311. 10.1016/j.exger.2010.01.007.View ArticlePubMedGoogle Scholar
- Chen W, Kang J, Xia J, Li Y, Yang B, Chen B, Sun W, Song X, Xiang W, Wang X, et al: p53-related apoptosis resistance and tumor suppression activity in UVB-induced premature senescent human skin fibroblasts. Int J Mol Med. 2008, 21 (5): 645-653.PubMedGoogle Scholar
- Dumaz N, Drougard C, Quilliet X, Mezzina M, Sarasin A, Daya-Grosjean L: Recovery of the normal p53 response after UV treatment in DNA repair-deficient fibroblasts by retroviral-mediated correction with the XPD gene. Carcinogenesis. 1998, 19 (9): 1701-1704. 10.1093/carcin/19.9.1701.View ArticlePubMedGoogle Scholar
- Di Leonardo A, Linke SP, Clarkin K, Wahl GM: DNA damage triggers a prolonged p53-dependent G1 arrest and long-term induction of Cip1 in normal human fibroblasts. Genes Dev. 1994, 8 (21): 2540-2551. 10.1101/gad.8.21.2540.View ArticlePubMedGoogle Scholar
- England NL, Cuthbert AP, Trott DA, Jezzard S, Nobori T, Carson DA, Newbold RF: Identification of human tumour suppressor genes by monochromosome transfer: rapid growth-arrest response mapped to 9p21 is mediated solely by the cyclin-D-dependent kinase inhibitor gene, CDKN2A (p16INK4A). Carcinogenesis. 1996, 17 (8): 1567-1575. 10.1093/carcin/17.8.1567.View ArticlePubMedGoogle Scholar
- Ahmed NU, Ueda M, Ichihashi M: Induced expression of p16 and p21 proteins in UVB-irradiated human epidermis and cultured keratinocytes. J Dermatol Sci. 1999, 19 (3): 175-181. 10.1016/S0923-1811(98)00068-1.View ArticlePubMedGoogle Scholar
- Medrano EE, Im S, Yang F, Abdel-Malek ZA: Ultraviolet B light induces G1 arrest in human melanocytes by prolonged inhibition of retinoblastoma protein phosphorylation associated with long-term expression of the p21Waf-1/SDI-1/Cip-1 protein. Cancer Res. 1995, 55 (18): 4047-4052.PubMedGoogle Scholar
- Fang Y, Hu XH, Jia ZG, Xu MH, Guo ZY, Gao FH: Tiron protects against UVB-induced senescence-like characteristics in human dermal fibroblasts by the inhibition of superoxide anion production and glutathione depletion. Australas J Dermatol. 2012, 53 (3): 172-180. 10.1111/j.1440-0960.2012.00912.x.View ArticlePubMedGoogle Scholar
- Tagawa H, Karube K, Tsuzuki S, Ohshima K, Seto M: Synergistic action of the microRNA-17 polycistron and Myc in aggressive cancer development. Cancer Sci. 2007, 98 (9): 1482-1490. 10.1111/j.1349-7006.2007.00531.x.View ArticlePubMedGoogle Scholar
- Wang Z, Liu M, Zhu H, Zhang W, He S, Hu C, Quan L, Bai J, Xu N: Suppression of p21 by c-Myc through members of miR-17 family at the post-transcriptional level. Int J Oncol. 2010, 37 (5): 1315-1321.View ArticlePubMedGoogle Scholar
- Biyashev D, Qin G: E2F and microRNA regulation of angiogenesis. Am J Cardiovasc Dis. 2011, 1 (2): 110-118.PubMed CentralPubMedGoogle Scholar
- Long J, Wang Y, Wang W, Chang BH, Danesh FR: Identification of microRNA-93 as a novel regulator of vascular endothelial growth factor in hyperglycemic conditions. J Biol Chem. 2010, 285 (30): 23457-23465. 10.1074/jbc.M110.136168.PubMed CentralView ArticlePubMedGoogle Scholar
- Choi OR, Lim IK: Loss of p21(Sdi1) expression in senescent cells after DNA damage accompanied with increase of miR-93 expression and reduced p53 interaction with p21(Sdi1) gene promoter. Biochem Biophys Res Commun. 2011, 407 (2): 406-411. 10.1016/j.bbrc.2011.03.038.View ArticlePubMedGoogle Scholar
- Fu X, Tian J, Zhang L, Chen Y, Hao Q: Involvement of microRNA-93, a new regulator of PTEN/Akt signaling pathway, in regulation of chemotherapeutic drug cisplatin chemosensitivity in ovarian cancer cells. FEBS Lett. 2012, 586 (9): 1279-1286. 10.1016/j.febslet.2012.03.006.View ArticlePubMedGoogle Scholar
- Klein U, Lia M, Crespo M, Siegel R, Shen Q, Mo T, Ambesi-Impiombato A, Califano A, Migliazza A, Bhagat G, et al: The DLEU2/miR-15a/16-1 cluster controls B cell proliferation and its deletion leads to chronic lymphocytic leukemia. Cancer Cell. 2010, 17 (1): 28-40. 10.1016/j.ccr.2009.11.019.View ArticlePubMedGoogle Scholar
- Aqeilan RI, Calin GA: Croce CM: miR-15a and miR-16-1 in cancer: discovery, function and future perspectives. Cell Death Differ. 2010, 17 (2): 215-220. 10.1038/cdd.2009.69.View ArticlePubMedGoogle Scholar
- Cao R, Zhang Y: SUZ12 is required for both the histone methyltransferase activity and the silencing function of the EED-EZH2 complex. Mol Cell. 2004, 15 (1): 57-67. 10.1016/j.molcel.2004.06.020.View ArticlePubMedGoogle Scholar
- Carvalho J, van Grieken NC, Pereira PM, Sousa S, Tijssen M, Buffart TE, Diosdado B, Grabsch H, Santos MA, Meijer G, et al: Lack of microRNA-101 causes E-cadherin functional deregulation through EZH2 up-regulation in intestinal gastric cancer. J Pathol. 2012, 228 (1): 31-44.PubMedGoogle Scholar
- He L, He X, Lim LP, de Stanchina E, Xuan Z, Liang Y, Xue W, Zender L, Magnus J, Ridzon D, et al: A microRNA component of the p53 tumour suppressor network. Nature. 2007, 447 (7148): 1130-1134. 10.1038/nature05939.PubMed CentralView ArticlePubMedGoogle Scholar
- Chekulaeva M, Filipowicz W: Mechanisms of miRNA-mediated post-transcriptional regulation in animal cells. Curr Opin Cell Biol. 2009, 21 (3): 452-460. 10.1016/j.ceb.2009.04.009.View ArticlePubMedGoogle Scholar
- Yan H, Wu J, Liu W, Zuo Y, Chen S, Zhang S, Zeng M, Huang W: MicroRNA-20a overexpression inhibited proliferation and metastasis of pancreatic carcinoma cells. Hum Gene Ther. 2010, 21 (12): 1723-1734. 10.1089/hum.2010.061.View ArticlePubMedGoogle Scholar
- Pin AL, Houle F, Guillonneau M, Paquet ER, Simard MJ: Huot J: miR-20a represses endothelial cell migration by targeting MKK3 and inhibiting p38 MAP kinase activation in response to VEGF. Angiogenesis. 2012, 15 (4): 593-608. 10.1007/s10456-012-9283-z.View ArticlePubMedGoogle Scholar
- Liu S, Patel SH, Ginestier C, Ibarra I, Martin-Trevino R, Bai S, McDermott SP, Shang L, Ke J, Ou SJ, et al: MicroRNA93 regulates proliferation and differentiation of normal and malignant breast stem cells. PLoS Genet. 2012, 8 (6): e1002751-10.1371/journal.pgen.1002751.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu XF, Zou J, Bao ZJ: Dong J: miR-93 suppresses proliferation and colony formation of human colon cancer stem cells. World J Gastroenterol. 2011, 17 (42): 4711-4717. 10.3748/wjg.v17.i42.4711.PubMed CentralView ArticlePubMedGoogle Scholar
- Fang L, Deng Z, Shatseva T, Yang J, Peng C, Du WW, Yee AJ, Ang LC, He C, Shan SW, et al: MicroRNA miR-93 promotes tumor growth and angiogenesis by targeting integrin-beta8. Oncogene. 2011, 30 (7): 806-821. 10.1038/onc.2010.465.View ArticlePubMedGoogle Scholar
- Xu D, He XX, Chang Y, Sun SZ, Xu CR, Lin JS: Downregulation of MiR-93 expression reduces cell proliferation and clonogenicity of HepG2 cells. Hepatogastroenterology. 2012, 59 (120): 2367-2373.PubMedGoogle Scholar
- Ivanovska I, Ball AS, Diaz RL, Magnus JF, Kibukawa M, Schelter JM, Kobayashi SV, Lim L, Burchard J, Jackson AL, et al: MicroRNAs in the miR-106b family regulate p21/CDKN1A and promote cell cycle progression. Mol Cell Biol. 2008, 28 (7): 2167-2174. 10.1128/MCB.01977-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong P, Iwasaki M, Somervaille TC, Ficara F, Carico C, Arnold C, Chen CZ, Cleary ML: The miR-17-92 microRNA polycistron regulates MLL leukemia stem cell potential by modulating p21 expression. Cancer Res. 2010, 70 (9): 3833-3842. 10.1158/0008-5472.CAN-09-3268.PubMed CentralView ArticlePubMedGoogle Scholar
- Ciccia A, Elledge SJ: The DNA damage response: making it safe to play with knives. Mol Cell. 2010, 40 (2): 179-204. 10.1016/j.molcel.2010.09.019.PubMed CentralView ArticlePubMedGoogle Scholar
- Finnerty JR, Wang WX, Hebert SS, Wilfred BR, Mao G, Nelson PT: The miR-15/107 group of microRNA genes: evolutionary biology, cellular functions, and roles in human diseases. J Mol Biol. 2010, 402 (3): 491-509. 10.1016/j.jmb.2010.07.051.PubMed CentralView ArticlePubMedGoogle Scholar
- Musumeci M, Coppola V, Addario A, Patrizii M, Maugeri-Sacca M, Memeo L, Colarossi C, Francescangeli F, Biffoni M, Collura D, et al: Control of tumor and microenvironment cross-talk by miR-15a and miR-16 in prostate cancer. Oncogene. 2011, 30 (41): 4231-4242. 10.1038/onc.2011.140.View ArticlePubMedGoogle Scholar
- Pedley J, Ablett EM, Pettit A, Meyer J, Dunn IS, Sturm RA, Parsons PG: Inhibition of retinoblastoma protein translation by UVB in human melanocytic cells and reduced cell cycle arrest following repeated irradiation. Oncogene. 1996, 13 (6): 1335-1342.PubMedGoogle Scholar
- Tong X, Van Dross RT, Abu-Yousif A, Morrison AR, Pelling JC: Apigenin prevents UVB-induced cyclooxygenase 2 expression: coupled mRNA stabilization and translational inhibition. Mol Cell Biol. 2007, 27 (1): 283-296. 10.1128/MCB.01282-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Varambally S, Cao Q, Mani RS, Shankar S, Wang X, Ateeq B, Laxman B, Cao X, Jing X, Ramnarayanan K, et al: Genomic loss of microRNA-101 leads to overexpression of histone methyltransferase EZH2 in cancer. Science. 2008, 322 (5908): 1695-1699. 10.1126/science.1165395.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao P, Deng Z, Wan M, Huang W, Cramer SD, Xu J, Lei M, Sui G: MicroRNA-101 negatively regulates Ezh2 and its expression is modulated by androgen receptor and HIF-1alpha/HIF-1beta. Mol Cancer. 2010, 9: 108-10.1186/1476-4598-9-108.PubMed CentralView ArticlePubMedGoogle Scholar
- Bracken AP, Kleine-Kohlbrecher D, Dietrich N, Pasini D, Gargiulo G, Beekman C, Theilgaard-Monch K, Minucci S, Porse BT, Marine JC, et al: The polycomb group proteins bind throughout the INK4A-ARF locus and are disassociated in senescent cells. Genes Dev. 2007, 21 (5): 525-530. 10.1101/gad.415507.PubMed CentralView ArticlePubMedGoogle Scholar
- Fan T, Jiang S, Chung N, Alikhan A, Ni C, Lee CC, Hornyak TJ: EZH2-dependent suppression of a cellular senescence phenotype in melanoma cells by inhibition of p21/CDKN1A expression. Mol Cancer Res. 2011, 9 (4): 418-429. 10.1158/1541-7786.MCR-10-0511.PubMed CentralView ArticlePubMedGoogle Scholar
- Tzatsos A, Paskaleva P, Lymperi S, Contino G, Stoykova S, Chen Z, Wong KK, Bardeesy N: A lysine (K)-specific demethylase 2B (KDM2B)-let-7-enhancer of zester homolog 2 (EZH2) pathway regulates cell cycle progression and senescence in primary cells. J Biol Chem. 2011, 286 (38): 33061-33069. 10.1074/jbc.M111.257667.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang X, Milyavsky M, Shats I, Erez N, Goldfinger N, Rotter V: Activated p53 suppresses the histone methyltransferase EZH2 gene. Oncogene. 2004, 23 (34): 5759-5769. 10.1038/sj.onc.1207706.View ArticlePubMedGoogle Scholar
- Tarasov V, Jung P, Verdoodt B, Lodygin D, Epanchintsev A, Menssen A, Meister G, Hermeking H: Differential regulation of microRNAs by p53 revealed by massively parallel sequencing: miR-34a is a p53 target that induces apoptosis and G1-arrest. Cell Cycle. 2007, 6 (13): 1586-1593. 10.4161/cc.6.13.4436.View ArticlePubMedGoogle Scholar
- Hutter E, Renner K, Pfister G, Stockl P, Jansen-Durr P, Gnaiger E: Senescence-associated changes in respiration and oxidative phosphorylation in primary human fibroblasts. Biochem J. 2004, 380 (Pt 3): 919-928.PubMed CentralView ArticlePubMedGoogle Scholar
- Dimri GP, Lee X, Basile G, Acosta M, Scott G, Roskelley C, Medrano EE, Linskens M, Rubelj I, Pereira-Smith O, et al: A biomarker that identifies senescent human cells in culture and in aging skin in vivo. Proc Natl Acad Sci U S A. 1995, 92 (20): 9363-9367. 10.1073/pnas.92.20.9363.PubMed CentralView ArticlePubMedGoogle Scholar
- Laschober GT, Ruli D, Hofer E, Muck C, Carmona-Gutierrez D, Ring J, Hutter E, Ruckenstuhl C, Micutkova L, Brunauer R, et al: Identification of evolutionarily conserved genetic regulators of cellular aging. Aging Cell. 2010, 9 (6): 1084-1097. 10.1111/j.1474-9726.2010.00637.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Castoldi M, Schmidt S, Benes V, Noerholm M, Kulozik AE, Hentze MW, Muckenthaler MU: A sensitive array for microRNA expression profiling (miChip) based on locked nucleic acids (LNA). RNA. 2006, 12 (5): 913-920. 10.1261/rna.2332406.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36 (Database issue): D154-D158.PubMed CentralPubMedGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMedGoogle Scholar
- Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27 (1): 91-105. 10.1016/j.molcel.2007.06.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, et al: Combinatorial microRNA target predictions. Nat Genet. 2005, 37 (5): 495-500. 10.1038/ng1536.View ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34 (Database issue): D140-D144.PubMed CentralView ArticlePubMedGoogle Scholar
- Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genet. 2007, 39 (10): 1278-1284. 10.1038/ng2135.View ArticlePubMedGoogle Scholar
- Gaidatzis D, van Nimwegen E, Hausser J, Zavolan M: Inference of miRNA targets using evolutionary conservation and pathway analysis. BMC Bioinforma. 2007, 8: 69-10.1186/1471-2105-8-69.View ArticleGoogle Scholar
- Miranda KC, Huynh T, Tay Y, Ang YS, Tam WL, Thomson AM, Lim B, Rigoutsos I: A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell. 2006, 126 (6): 1203-1217. 10.1016/j.cell.2006.07.031.View ArticlePubMedGoogle Scholar
- Maragkakis M, Reczko M, Simossis VA, Alexiou P, Papadopoulos GL, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, et al: DIANA-microT web server: elucidating microRNA functions through target prediction. Nucleic Acids Res. 2009, 37 (Web Server issue): W273-W276.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang JC, Babak T, Corson TW, Chua G, Khan S, Gallie BL, Hughes TR, Blencowe BJ, Frey BJ, Morris QD: Using expression profiling data to identify human microRNA targets. Nat Methods. 2007, 4 (12): 1045-1049. 10.1038/nmeth1130.View ArticlePubMedGoogle Scholar
- Saito T, Saetrom P: MicroRNAs-targeting and target prediction. N Biotechnol. 2010, 27 (3): 243-249. 10.1016/j.nbt.2010.02.016.View ArticlePubMedGoogle Scholar
- Sturn A, Quackenbush J, Trajanoski Z: Genesis: cluster analysis of microarray data. Bioinformatics. 2002, 18 (1): 207-208. 10.1093/bioinformatics/18.1.207.View ArticlePubMedGoogle Scholar
- Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman WH, Pages F, Trajanoski Z, Galon J: ClueGO: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009, 25 (8): 1091-1093. 10.1093/bioinformatics/btp101.PubMed CentralView ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.