- Research article
- Open access
- Published:
Differential miRNA expression profiles in proliferating or differentiated keratinocytes in response to gamma irradiation
BMC Genomics volume 14, Article number: 184 (2013)
Abstract
Background
MicroRNAs (miRNAs), a group of short non-coding RNAs that negatively regulate gene expression, have recently emerged as potential modulators of cellular response to ionizing radiations both in vitro and in vivo in various cell types and tissues. However, in epidermal cells, the involvement of the miRNA machinery in the cellular response to ionizing radiations remains to be clarified. Indeed, understanding the mechanisms of cutaneous radiosensitivity is an important issue since skin is the most exposed organ to ionizing radiations and among the most sensitive.
Results
We settled up an expression study of miRNAs in primary human skin keratinocytes using a microfluidic system of qPCR assay, which permits to assess the expression of almost 700 annotated miRNAs. The keratinocytes were cultured to a proliferative or a differentiated state mimicking basal or suprabasal layers of human epidermis. These cells were irradiated at 10 mGy or 6 Gy and RNA was extracted 3 hours after irradiation. We found that proliferative cells irradiated at 6 Gy display a global fall of miRNA expression whereas differentiated cells exposed to the same dose display a global increase of miRNAs expression. We identified twenty miRNAs weakly but significantly modulated after 6 Gy irradiation, whereas only 2 miRNAs were modulated after low-dose irradiation in proliferating cells. To go further into the biological meaning of this miRNA response, we over-expressed some of the responding miRNA in proliferating cells: we observed a significant decrease of cell viability 72 hours after irradiation. Functional annotation of their predicted targets revealed that G-protein related pathways might be regulated by these responding miRNAs.
Conclusions
Our results reveal that human primary keratinocytes exposed to ionizing irradiation expressed a miRNA pattern strongly related to the differentiation status of irradiated cells. We also demonstrate that some miRNAs play a role in the radiation response to ensure the short-term survival of irradiated keratinocytes.
Background
Ionizing radiations (IR) induce a large spectrum of damages that irradiated cells have to counteract by activating complex response pathways. The gene networks mobilized in the cellular response to radiations are subtly regulated. Several molecular actors have been shown to be involved in the regulation of gene expression after irradiation including chromatin-remodeling factors [1–3], transcription factors [4–6] or factors controlling protein synthesis [7], maturation [8] and degradation [9]. MicroRNAs (miRNAs) are short non-coding RNAs (21 to 23-nucleotide-long) which are predicted to control in mammals the activity of approximately 30% of all protein coding genes and have been shown to participate in the regulation of almost any cellular processes investigated so far [10]. By base pairing to the 3′ untranslated regions of target mRNAs, miRNAs mediate translational repression during the initiation or elongation step, and proteolysis of the nascent peptide or mRNA degradation after deadenylation [11].
miRNAs have been proposed to be important actors of the DNA damage response [12]. Several miRNAs expression studies have been conducted after irradiation in mice [13], human tumor samples [14], cancer cell lines [15, 16] or normal human cells [15–18]. By crossing these studies, it was not possible to identify a short-list of universal miRNAs responding to IR. The radiation effect on miRNA expression seems therefore to vary according to cell type, radiation dose and post-irradiation time point [19].
Skin is the first organ to be targeted by external radiations and among the most sensitive. The epidermis is a stratified squamous epithelium that forms the protective covering of the skin. It is composed of a basal layer of proliferating cells, a spinous postmitotic cell layer, a granular layer, and a stratum corneum of terminally differentiating keratinocytes. High doses of radiation can induce various deleterious effects in skin, which can appear a short time after exposure, such as erythema and desquamation, and several years later, such as carcinoma and late complications [20]. Radiation-induced damages in epidermal proliferating cells have been suspected to be involved in late effects such as skin cancer whereas defects in differentiated cells could lead to early effects such as hyperkeratosis or desquamation. It is therefore of great interest to investigate the role of miRNAs in the radiation response in proliferating versus differentiated keratinocytes. Although the biological effects of low doses are controversial, we have previously shown that a very low dose of IR is able to induce a strong and specific mRNA response in human keratinocytes [21]. By providing a global view of the miRNA response after strong and low doses of IR, our study should help answering in which condition a miRNA response is induced after irradiation.
In this study, miRNA expression profiling was investigated in human primary keratinocytes after acute exposure to a very low (10 mGy) or a strong (6 Gy) dose of gamma rays. We showed that radiation induced changes in miRNA levels varied according to the dose and the differentiation status of irradiated keratinocytes. Indeed, after a strong dose of 6 Gy, a global repression of miRNAs was observed in proliferating keratinocytes whereas a global induction was detected in differentiated keratinocytes. The very low dose of 10 mGy was also able to modulate some miRNAs in proliferating cells. Over-expression of the most stringently regulated miRNAs in proliferating cells modified their survival after 6 Gy irradiation. In differentiated keratinocytes, repression of the 3 most activated miRNAs had no detectable effect on survival. Our data suggest that the miRNA response is globally modulated in irradiated keratinocytes and could be involved in cellular pathways that are highly dependent to the differentiation state of the IR-exposed cells.
Results
Characterization of the cellular model
The cellular model used in this study consisted of primary cultures of human keratinocytes (HPK) isolated from infant foreskin. For modeling epidermal basal cells, second passage cultures were cultivated up to 50% of confluency where most of the cells are still proliferating. To induce epidermal differentiation, cells were cultured at low calcium concentration up to confluence, and then incubated 3 days post-confluence in a culture medium containing 1.8 mM of calcium. Differentiated cultures were subjected to expression analysis by real-time PCR of two well-known markers of epidermal suprabasal layers: keratin 1 (KRT1) and involucrin (IVL) as well as markers of proliferating basal cells: keratin 14 (KRT14) and PCNA. The expression was compared to that of proliferative cultures. As expected, a strong expression of KRT1 and IVL transcripts was observed in differentiated keratinocytes whereas PCNA and KRT14 were clearly repressed (Figure 1A). These results suggest that the majority of cells cultured 3 days post-confluence in a high-calcium medium exhibited a differentiation state close to that found in spinous and granular layers. On the contrary, semi-confluent keratinocytes express markers typical to the basal layers of the human epidermis. We next investigated the short-term radiosensitivity of proliferating and differentiated keratinocytes after the acute high dose of 6 Gy. The viability of irradiated cells was measured 24 to 72 hours after irradiation (Figures 1B and 1C). For proliferating cells (Figure 1B), we observed a significant decrease in cell viability starting 24 h post irradiation (Student t-test, P < 0.05) with a maximal effect intensity of 20% decrease 48 h and 72 h post-irradiation (Student t-test, P < 10-2). In differentiated cells (Figure 1C), a weaker but significant effect was observed from 24 h to 72 h post-irradiation with a maximal decrease in cell viability of 10% 3 days post-irradiation (Student t-test, P < 0.05). These data confirmed, as previously reported by our group [21], that differentiated keratinocytes are more radio-resistant than proliferative keratinocytes.
MiRNA expression profiling
We used TaqMan quantitative real-time PCR Low Density Array (TLDA) to profile miRNAs expression in control and irradiated keratinocytes. We investigated the miRNA response to a high dose of 6 Gy in proliferating and differentiated keratinocytes and to a very low dose of 10 mGy in proliferating cells, which are more radiosensitive. MiRNAs expression was assessed 3 hours post-irradiation since this time-point corresponds to the immediate gene response in irradiated keratinocytes [22] and miRNAs might be involved in the regulation of this gene response. Three to four human keratinocytes cultures established from independent donors were analyzed at each dose. Due to the high repeatability of the TLDA methodology [23, 24], we privileged in our experimental plan biological replicates from different donors than technical repeats. It is important to note that the experimental plan is constructed on paired data: for each donor, a reference condition response and a response following ionizing irradiation were quantified.
Out of 667 miRNAs studied, only those with a Cq (quantification Cycle) value below 32 cycles in all samples were considered as significantly expressed. This threshold was calculated using miRNA duplicate and snRNA (RNU6-1, SNORD44, SNORD48) tetraplicate repeatability, in agreement with literature [25]. 95 miRNAs for 0 Gy vs 6 Gy (0/6 Gy) irradiated proliferating keratinocytes, 100 miRNAs for 0/6 Gy-irradiated differentiated keratinocytes, and 64 miRNAs for 0 Gy vs 10 mGy (0/10 mGy) irradiated proliferating keratinocytes were therefore considered for further analysis (see Additional File 1 for a complete list of these miRNAs and their corresponding Cq in the various conditions). We first investigated if a bias exists in the average expression of miRNAs between HPK cell samples sorted by differentiation state and irradiation dose (Figure 2). We used all-paired Tukey-Kramer tests to compare Cq means within HPK samples. For 0/6 Gy-irradiated proliferating keratinocytes, there was no expression bias between the cell samples (ANOVA F-test, p-value = 0.578 and 0.813 respectively). For 0/6 Gy-irradiated differentiated keratinocytes, significant differences were observed: samples from HPK86 revealed higher Cq than the three others donors (ANOVA F-test, P < 10-4 and 6.10-4 respectively). A difference was also observed for 0/10 mGy-irradiated proliferating keratinocytes: samples from donor HPK65 showed in average lower Cq, than those from HPK55 and HPK98 (ANOVA F-test, P = 5.10-4 and < 10-4 respectively). However, these differences in the average expression of miRNAs were similar for reference and irradiated conditions (compare panels A to panels B), and could be due to inter-individual variability in primary cells. This observation underlined the importance to use paired data for statistical analyses. It further confirmed that there was no sample processing bias due to irradiation state and that raw data normalization was not necessary as previously proposed [26].
We then focused on miRNA profiles after irradiation in proliferating and differentiated keratinocytes. After the high dose of 6 Gy, a global trend of repression was observed in proliferating keratinocytes whereas an opposite effect was observed in differentiated cells where most of the expressed miRNAs were induced after irradiation (Figure 2C and Figure 3 top and medium panel). After the very low dose, a slight increase of the miRNA expression profile was observed in proliferating keratinocytes (Figure 2C and Figure 3 bottom panel). These data revealed that the miRNA response follows a global trend that is highly dependent on the radiation dose and the differentiation status of irradiated cells.
Individual miRNAs significantly modulated after irradiation among the 3 or 4 independent cell samples were identified (colored dots in Figure 3). These miRNAs covered a wide range of expression from highly (Cq around 22 cycles) to weakly expressed (Cq around 31 cycles). After 6 Gy, 8 and 12 miRNAs were found to be significantly modulated in proliferating and differentiated cells respectively (Figure 4A and B). As expected with regards to the global response, miRNAs modulated in proliferating cells were all repressed whereas 11 out 12 miRNAs modulated in differentiated cells were induced. The miRNAs responding to the high dose of 6 Gy in proliferating cells were all different from those induced in differentiated keratinocytes. After 10 mGy, only 2 induced miRNAs were identified so far in proliferating keratinocytes (Figure 4C).
To strengthen the validity of the data obtained from the large-scale TLDA analysis, we used individual QPCR to study the expression of 6 miRNAs in proliferating or differentiated keratinocytes among a time course between 1 and 24 hours post 6 Gy-irradiation. In proliferating keratinocytes (Figure 5A), we confirmed the down-regulation of the 3 tested miRNAs (miR-191-5p, miR-200b-3p and miR-331-3p) at 3 h whereas no modulation was detected at 1 or 24 h. In differentiated keratinocytes (Figure 5B), we also observed a significant response at 3 h for two of the 3 tested miRNAs: as observed after TLDA analysis (Figure 4B), let-7b-5p and miR-196b-5p were up-regulated 3 hours after irradiation. Globally, these data are in accordance with those obtained from the large-scale TLDA analysis. They also support the choice of the 3 hours time-point as a key moment in the miRNA response of irradiated keratinocytes.
Towards mechanisms of the miRNA response
To go further into the mechanisms of this miRNA response after IR, we analyzed the effect of radiations on two proteins known to involved in the processing of miRNAs: DICER1 and AGO2 [18]. Indeed, DICER1 is an RNAse III endonuclease that acts in the endonucleolytic processing of miRNAs by cleaving their loop to release a duplex RNA stem [27]. This duplex RNA is incorporated into the RISC complex, a multi-protein complex that also contains the Argonaute proteins (AGO1 to AGO4 in humans). Therefore, we analyzed the mRNA response of DICER1 and AGO2 genes in proliferating and differentiated keratinocytes 3 hours after 6 Gy irradiation. In proliferating cells, we observed a faint but significant decrease of mRNA expression level for both DICER1 and AGO2 (Student t-test, P < 0.05 and < 10-3 respectively) (Figure 6A), whereas protein amount was not significantly different (Figure 6C). In differentiated cells, no modulation was detected either at RNA or protein level for both DICER1 and AGO2 (Figure 6B and D). These data suggest that the late steps of the miRNA-processing pathway, that require AGO2 and DICER, are not directly impacted by ionizing radiations. To have a look to the first step of miRNA biogenesis, we analyzed the relative expression of the primary transcripts corresponding to 3 repressed miRNAs in proliferating keratinocytes. A significant down-regulation was observed for primiR-331 and primiR-200b in 6 Gy-irradiated keratinocytes (Figure 6E). These results indicate that some IR-responding miRNAs are directly regulated at the transcriptional level. As miR-191-5p is down-regulated whereas the level of its precursor remains constant, those results suggest that the global modulation of the miRNAs expression that we detected after irradiation (see Figure 3) is the result of multiple and complementary levels of regulation within the miRNA biogenesis pathway.
Functional relevance of IR-responding miRNAs
After identification of radiation-responsive miRNAs, our aim was to test their individual or global impact on the survival of irradiated cells. We decided to focus first on proliferating cells that are more radio-sensitive (see Figure 1B). Because of their higher level of repression in irradiated proliferating cells, miR-191-5p, miR-200b-3p and miR-331-3p were selected for further analysis. These 3 miRNAs were individually over-expressed before irradiation by transfection of the corresponding pre-miRNA. In parallel, a pool of the 3 miRNAs was also over-expressed in proliferating keratinocytes prior to irradiation. The experimental conditions were adapted to ensure a moderate over-expression (5 to 10 times) of the miRNAs in transfected cells (Figure 7A). Over-expression of individual miRNAs had no effect on cell survival (except for miR-331-3p at 48 h) whereas a combination of the 3 miRNAs leads to a significant reduction of cell viability especially at 72 h post-irradiation (Figure 7B and Table 1). We also observed that some combined over-expression of two of the 3 miRNAs were also able to significantly reduced cell survival at 72 h (Additional File 2). These data suggest that the miRNAs that are modulated in irradiated cells plays a role in the cellular radiation response through a combined mode of action. This action is probably exerted through the regulation of downstream targets genes. To go further into these mechanisms, we searched for putative gene targets of miR-191-5p, miR-200b-3p and miR-331-3p. By crossing the results of five prediction programs (miRanda, miRDB, miRWalk, RNA22 and TargetScan), we identified hundreds of genes targets for each miRNAs. Among these genes, the 3 miRNAs share 83 targets. 783 genes are potentially regulated by at least 2 of the 3 miRNAs (see Venn diagram in Additional File 3 and the complete list of these targets in Additional File 4). A functional annotation of these 783 putative targets using GOrilla (Gene Ontology enrichment analysis and vizualization tool) revealed an enrichment of genes involved in the modulation of G-protein coupled receptor-signaling pathway through adenylate cyclase (Enrichment p-value = 9.8 10-6, FDR q-value = 3.9 10-2).
We used the same strategy to explore the functional impact of the 6 Gy-responding miRNAs in differentiated cells. A combined silencing of the 3 more induced sequences (let7b-5p, miR-196b-5p and miR-19a-3p) by transfection of the corresponding antagomiR had no effect on cell survival at 72 h post-irradiation (data not shown).
Discussion
In this work, we identified several miRNAs responding to IR in human keratinocytes. In proliferating cells exposed to 6 Gy, we identified 8 miRNAs that are repressed 3 hours after exposure. Among these 8 miRNAs, only miR-17-5p [28] and miR-31-5p [29] have been previously found to be deregulated by IR in other primary or cancer cells. Moreover, miR-31-5p has been shown to repress the expression of several genes involved in DNA repair in adenocarcinoma cell lines [29]. Its down-regulation in proliferating keratinocytes might contribute to increase the DNA repair efficiency after irradiation. MiR-99b-5p was also described as playing a role in DNA repair. The miR-99 family, including miR-99b-5p, regulates the DNA-damage response in breast and prostate cancer cells by targeting the chromatin remodeling factor SNF2H [30]. When the miR-99 cluster is over-expressed in irradiated cells, the rate and the overall efficiency of repair by both NHEJ and homologous recombination are reduced [30]. In our study, the down-regulation of miR-99b-5p could contribute to an activation of double-strand breaks repair in irradiated keratinocytes. Irradiation of proliferating keratinocytes leads to a cell cycle arrest that is requested for DNA repair [31]. Some of the miRNAs repressed by IR in proliferating keratinocytes could be involved in this growth arrest: it is the case of miR-17-5p and miR-191-5p that are over-expressed in cancer cells where they contribute to activate cell proliferation [32, 33]. In differentiated keratinocytes exposed to 6 Gy, we identified 12 responding miRNAs, all except miR-125b-5p were induced in irradiated cells. Among these 12 miRNAs only let-7b-5p has been previously described as being modulated by ionizing radiations in primary or cancer cells [34]. This suggests that the miRNA response elicited by IR in human keratinocytes is very specific to this cell type and also dependent to the differentiation status of irradiated keratinocytes. In a recent study, Zhou et al. investigated the miRNA profile of UVB-irradiated normal human keratinocytes. Among the 44 miRNAs that they found modulated by UVB, only 4 (miR-31-5p, let-7b-5p, miR-125b-5p and miR-186-5p) also responded to IR in our study [35]. This could indicate that UVB and IR provoke specific effects in exposed keratinocytes that in turn elicit a stress-specific miRNA response. It has been previously observed that genotoxic stresses, such as UV or IR, induce a very specific response in differentiated keratinocytes including induction of terminal differentiation [36, 37]: among the 12 miRNAs that we found induced after irradiation in differentiated keratinocytes, 5 of them (miR-125b-5p, let-7b-5p, miR-181a-5p, miR-195-5p and miR-342-3p) have been previously observed as being also induced during the differentiation process of human keratinocytes [38]. These 5 radiation-responding miRNAs could be directly involved into the modification of the differentiation program that is elicited by gamma irradiation in differentiated keratinocytes.
The main question raised by our data is the biological significance of this miRNA response. We have re-expressed 3 down-regulated miRNAs in irradiated proliferating keratinocytes. We observed that a combined over-expression of these miRNAs reduced cell viability 72 hours post-IR. These results indicate that the down-regulation of these miRNAs is important for cell survival after IR. This could also evocate a coordinated function of these miRNAs in the cellular radiation response probably through the regulation of common genes targets. We identified more than 780 potential targets in common between at least two of these 3 responding miRNAs. Functional annotation of this list revealed an over-representation of genes involved in the adenylate cyclase G-protein coupled receptor-signaling pathway. The cAMP signaling pathway has been shown to modulate DNA-damaging agents induced apoptosis and DNA repair activity [39, 40]. This pathway could be one of the mechanisms downstream to ionizing-radiation-responding miRNAs and necessary for the survival of irradiated cells. In differentiated keratinocytes, the silencing of the 3 more responding miRNAs had no effects on immediate cell viability: this suggests that the miRNAs responding to IR in differentiated keratinocytes are not directly involved in immediate cell survival. Since these cells are already engaged in a differentiation process that will finally lead to cornification, a specialized form of programmed cell death, the activation of the differentiation process might be a way to eliminate damaged cells. As discussed above, this hypothesis was supported by the fact that, in differentiated keratinocytes, several miRNAs that we found induced by IR are also markers of epidermal differentiation.
Several studies have been conducted to study the effects of low doses ranging from 50 mGy to 100 mGy on miRNAs expression in various cellular models including fibroblasts [17, 41], blood cells [42] or thyroid cells [43]. However, to our knowledge, the present study is the first to address the miRNA response of human cells after a very low dose of 10 mGy. This dose is relevant in the context of medical imaging procedures where 20% of the patients can be potentially exposed to radiation between 3 and 20 mGy, especially after abdominal compute tomography [41]. Contrary to the high dose of 6 Gy where all the responding miRNAs are down-regulated in proliferating keratinocytes, the very low dose of 10 mGy leads to a particular response with only two up-regulated miRNAs: miR-342-3p and miR-708-5p. As previously observed for mRNA [21], a specific miRNA response is then detected in human keratinocytes after the very low dose of 10 mGy. This might reflect the emerging concept that the dose response is not linear and that very low doses of IR induce specific cellular mechanisms in irradiated cells [44]. Further experiments will be necessary to precise how miR-342-3p and miR-708-5p might be involved in this particular response.
Conclusion
Our results suggest that in human keratinocytes the expression of miRNAs in response to ionizing radiation is highly dependent to the differentiation status of irradiated cells. It is also clear from our study that the miRNAs response is dose dependent and that a very low dose of 10 mGy is able to modulate some miRNAs. Since deregulation of specific miRNAs has the consequence of changing cell survival after irradiation, this might represents a possible strategy to modulate cellular sensitivity to anti-cancer treatments.
Methods
Ethical considerations
Infant foreskins were collected according to the Declaration of Helsinki Principles. A written informed consent was obtained from infants’ parents according to french bioethical law of 2004 (loi 94–654 du 29 juillet 1994).
Cell culture
Foreskins from healthy boys were conserved at 4°C in DMEM supplemented by 10% FBS and penicillin/streptomycin till keratinocytes preparation. Keratinocyte were isolated from epidermis using overnight 4°C dispase/trypsin digestion. Cells were seeded on BioCoat™ Collagen I cellware (Becton-Dickinson) in KGM2 medium (PromoCell) supplemented with 100 μg.ml-1 antimicrobial Primocin™ (InvivoGen). Subconfluent cultures were passed using 0.05% trypsin/EDTA (GIBCO) and reseeded three-time up for studies in KGM2 medium. Keratinocytes differentiation was induced at confluence adding CaCl2 to KGM2 medium, raising calcium concentration from 0.06 to 1.8 mM.
Samples irradiation
Proliferating and three-days differentiated keratinocytes were 6 Gy-irradiated using a Cs137 source (IBL637 Cisbio international) with a dose-rate of 0.62 Gy.min-1 calculated on radioactive decay. Proliferating keratinocytes were 10 mGy-irradiated at the French army health service research facility (CRSSA) on IRDI 4000 (Alstom) with a dose-rate of 21.68 mGy.min-1 calculated on radioactive decay. Samples irradiation was performed in biological triplicate or tetraplicate.
RNA extraction and quality evaluation
Total RNA was prepared using mirVana™ miRNA Isolation Kit (Applied Biosystems) according to the manufacturer’s instructions and conserved at −80°C till use. RNA purity and concentration were assayed using a NanoDrop 2000 Spectrophotometer (Thermo Scientific). All samples presented A260/A280 ratios > 2 and A260/A230 ratios > 1.9 in DEPC-treated water. RNA integrity was checked using a 2100 Bioanalyser (Agilent). All samples harbored correct 18S/28S profile and a RIN factor between 8.9 and 10.
MiRNA profiling and analysis
RT-qPCR experiments were performed by profileXpert platform (Université Claude-Bernard Lyon 1). The quantification of the 667 human miRNAs was performed using TaqMan® Array Human MicroRNA Cards A and B v2.0 (4400239, Applied biosystems) according to manufacturer instructions. Briefly, miRNAs were retrotranscribed by Megaplex™ RT Primers kit. RT primers are based on miRBase v14. Complementary DNA (cDNA) were plated on cards and PCR were performed on 7900HT Fast Real-Time PCR System (Applied Biosystems). Raw Cq (quantification Cycle) values were calculated using the RQ Manager version 1.2 applying automatic baseline settings and a threshold of 0.15. Plant ath-miR159a served as negative control. Cq threshold of 32 was calculated using miRNA duplicate of card B and snRNA tetraplicate repeatability, confirmed by literature [25]. Only miRNA with Cq below the threshold in all condition by state and irradiation dose were sorted. Data analysis was performed using JMP 5.0.1.2 software (SAS Institute). For miRNA targets prediction, we used five different miRNA target prediction tools: miRanda, miRDB, miRWalk, RNA22 and TargetScan. From all predictions, we considered miRNA-transcript relationships that were predicted by a minimum of two different tools. In addition, we performed a Gene Ontology term enrichment analysis using the GOrilla application [45]. To correct for multiple testing, we used the Benjamini-Hochberg procedure and considered an association significant if the P value was less than 0.05.
Reverse transcription and quantitative PCR for mRNA
Contaminating DNA was removed from 500 ng of RNA using DNaseI (EN0521, Fermentas) at 0.1U.μl-1 in 10 μl at 37°C during 30 min. DNA-free RNA was reverse transcribed using SuperScript™ II Reverse Transcriptase (18064–014, Invitrogen) with 17 ng of oligo-d(T) (Amersham Pharmacia Biotech Inc.) and 5 ng of random hexamers d(N)6 (N8080127, Invitrogen). RT products were diluted to an equivalent RNA concentration of 10 ng per well. Samples and standards were manually plated in duplicate on 96 wells polypropylene plates (410088, Agilent Technologies). Non-template controls were performed for each assay in all experiment and displayed not detectable or significant amplification. Quantitative PCR (qPCR) were performed using ABsolute™ Blue QPCR SYBR® Green ROX Mix (AB-4163/A, ThermoScientific) with Mx3000 system (Stratagene). Primer design was made using PrimerBlast interface (http://www.ncbi.nlm.nih.gov/tools/primer-blast/). PCR reactions were initiated by the activation of Taq DNA polymerase at 95°C for 15 min, followed by 45 three-step amplification cycles consisting of 15 s denaturation at 95°C, 30 s annealing at 58-60°C and 30 s extension at 72°C. After the amplification, a dissociation stage was run to generate a melting curve for verification of amplification product specificity, which was also confirmed on agarose gels. ROX normalized and automatic baseline corrected quantification cycles (Cq) were collected using MxPro software (version 3.00b, Stratagene). Amplification efficiencies were calculated using linear regression of cDNA pool dilution and between 90% and 105%. Standard minus square correlations were above 0.995. Relative quantification was calculated using standard curve efficiency-corrected comparative quantification [46] using Microsoft® Excel® Mac 2008. 18S served as differentiation-related reference gene. Ubiquitin c (UBC) served as ionizing irradiation-related reference gene as one of the most invariant gene in this stress [21]. The nucleotide sequence and the amplicon length for of all the primers used in this study are given in the Additional file 5. The MIQE checklist [47] of our QPCR experiments is given is the Additional file 6.
Reverse transcription and quantitative PCR for miRNA and primiRNAs
PrimiRNA were reverse transcripted from 500 ng of total RNA using the PrimeScript™ RT reagent kit (RR037A, Takara) according to the manufacturer instructions. Mature miRNAs from 100 μg of total RNA were retrotranscribed using TaqMan® MiRNA Reverse Transcription Kit (4366596, Applied Biosystems™) according to the manufacturer protocol. QPCR were performed in duplicate. Non-template controls and No-RT controls were performed for each assay in all experiment and display not detectable amplification. Quantitative PCR were performed using TaqMan® primiRNA or miRNA Assays and TaqMan® Universal PCR Master Mix (Applied Biosystems). Relative quantification was calculated using 2ΔΔ Cq quantification method [48]. The UBC gene, that we have previously shown to be invariant after irradiation in human keratinocytes [21], was used as a reference to normalize QPCR for primiRNA quantification. MiR-125b-5p and miR-17-5p were used as reference to normalize individual QPCR after irradiation on proliferating and differentiated keratinocytes respectively. RNU6-1 served as reference for QPCR after pre-miR transfection assays.
Ionizing irradiation viability assay
Viability was performed with resazurin/resorufin-based alamarBlue® cell viability assay (Life Technologies™). Keratinocytes were seeded at 8000 cell per cm2. The day after, cells were transfected during 6 h with 50pM pre-miR-191-5p, 50pM pre-miR-331-3p or 500fM pre-miR-200b-3p, individually or pooled, in tetraplicate using HiPerfect Transfection Reagent (301705, Qiagen). Days after transfection, 10% v/v alamarBlue was added to culture medium and fluorescence was read in technical duplicate 4 h after addition using Xenius instrument (Safas).
Western blot
Protein extraction was performed using RadioImmunoPrecipitation Assay medium (RIPA, [49]). Protein concentrations were measured by Bradford method using Bio-Rad Protein Assay (500–0006, Bio-Rad). Electrophoresis was performed with 20 ng of protein. Western blot were performed as described previously [4]. Results were obtained from four independent donors. EIF2C2 (ab57113) and DICER1 (ab14601) antibodies were purchased from abcam® and used after 1/200 and 1/1000 dilution respectively. ACTB, clone C4 (MAB1501) antibody was purchased from Merck Millipore and used after 1/2000 dilution.
Abbreviations
- IR:
-
Ionizing radiations
- miRNA:
-
microRNA
- HPK:
-
Human Primary Keratinocytes
- RT-QPCR:
-
Reverse Transcriptase Quantitative Polymerase Chain reaction
- TLDA:
-
Taqman quantitative PCR Low Density Array
- Cq:
-
Quantification Cycle
- SD:
-
Standard Deviation
- PCNA:
-
Proliferating Cell Nuclear Antigen
References
Iyengar S, Farnham PJ: KAP1 protein: an enigmatic master regulator of the genome. J Biol Chem. 2011, 286 (30): 26267-26276. 10.1074/jbc.R111.252569.
Ziv Y, Bielopolski D, Galanty Y, Lukas C, Taya Y, Schultz DC, Lukas J, Bekker-Jensen S, Bartek J, Shiloh Y: Chromatin relaxation in response to DNA double-strand breaks is modulated by a novel ATM- and KAP-1 dependent pathway. Nat Cell Biol. 2006, 8 (8): 870-876. 10.1038/ncb1446.
Smeenk G, Wiegant WW, Vrolijk H, Solari AP, Pastink A, van Attikum H: The NuRD chromatin-remodeling complex regulates signaling and repair of DNA damage. J Cell Biol. 2010, 190 (5): 741-749. 10.1083/jcb.201001048.
Bonin F, Molina M, Malet C, Ginestet C, Berthier-Vergnes O, Martin MT, Lamartine J: GATA3 is a master regulator of the transcriptional response to low-dose ionizing radiation in human keratinocytes. BMC Genomics. 2009, 10: 417-10.1186/1471-2164-10-417.
Pawlik A, Alibert O, Baulande S, Vaigot P, Tronik-Le Roux D: Transcriptome characterization uncovers the molecular response of hematopoietic cells to ionizing radiation. Radiat Res. 2011, 175 (1): 66-82. 10.1667/RR2282.1.
Sashida G, Bae N, Di Giandomenico S, Asai T, Gurvich N, Bazzoli E, Liu Y, Huang G, Zhao X, Menendez S: The mef/elf4 transcription factor fine tunes the DNA damage response. Cancer Res. 2011, 71 (14): 4857-4865. 10.1158/0008-5472.CAN-11-0455.
Hayman TJ, Williams ES, Jamal M, Shankavaram UT, Camphausen K, Tofilon PJ: Translation initiation factor eIF4E is a target for tumor cell radiosensitization. Cancer Res. 2012, 72 (9): 2362-2372. 10.1158/0008-5472.CAN-12-0329.
Meerang M, Ritz D, Paliwal S, Garajova Z, Bosshard M, Mailand N, Janscak P, Hubscher U, Meyer H, Ramadan K: The ubiquitin-selective segregase VCP/p97 orchestrates the response to DNA double-strand breaks. Nat Cell Biol. 2011, 13 (11): 1376-1382. 10.1038/ncb2367.
Huang J, Yuan H, Lu C, Liu X, Cao X, Wan M: Jab1 mediates protein degradation of the Rad9-Rad1-Hus1 checkpoint complex. J Mol Biol. 2007, 371 (2): 514-527. 10.1016/j.jmb.2007.05.095.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nature reviews. 2008, 9 (2): 102-114.
Hu H, Gatti RA: MicroRNAs: new players in the DNA damage response. J Mol Cell Biol. 2011, 151-158. 3
Koturbash I, Zemp F, Kolb B, Kovalchuk O: Sex-specific radiation-induced microRNAome responses in the hippocampus, cerebellum and frontal cortex in a mouse model. Mutat Res. 2010, 722 (2): 114-118.
Wang XC, Du LQ, Tian LL, Wu HL, Jiang XY, Zhang H, Li DG, Wang YY, Wu HY, She Y: Expression and function of miRNA in postoperative radiotherapy sensitive and resistant patients of non-small cell lung cancer. Lung Cancer. 2011, 72 (1): 92-99. 10.1016/j.lungcan.2010.07.014.
Josson S, Sung SY, Lao K, Chung LW, Johnstone PA: Radiation modulation of microRNA in prostate cancer cell lines. Prostate. 2008, 68 (15): 1599-1606. 10.1002/pros.20827.
Marsit CJ, Eddy K, Kelsey KT: MicroRNA responses to cellular stress. Cancer Res. 2006, 66 (22): 10843-10848. 10.1158/0008-5472.CAN-06-1894.
Maes OC, An J, Sarojini H, Wu H, Wang E: Changes in MicroRNA expression patterns in human fibroblasts after low-LET radiation. J Cell Biochem. 2008, 105 (3): 824-834. 10.1002/jcb.21878.
Kraemer A, Anastasov N, Angermeier M, Winkler K, Atkinson MJ, Moertl S: MicroRNA-mediated processes are essential for the cellular radiation response. Radiat Res. 2011, 176 (5): 575-586. 10.1667/RR2638.1.
Joly-Tonetti N, Lamartine J: The role of MicroRNAs in the cellular response to ionizing radiations. Current topics in ionizing radiation research. 2012, 1: 149-174.
Ryan JL: Ionizing radiation: the good, the bad, and the ugly. J Invest Dermatol. 2012, 132 (3 Pt 2): 985-993.
Franco N, Lamartine J, Frouin V, Le Minter P, Petat C, Leplat JJ, Libert F, Gidrol X, Martin MT: Low-dose exposure to gamma rays induces specific gene regulations in normal human keratinocytes. Radiat Res. 2005, 163 (6): 623-635. 10.1667/RR3391.
Lamartine J, Franco N, Le Minter P, Soularue P, Alibert O, Leplat JJ, Gidrol X, Waksman G, Martin MT: Activation of an energy providing response in human keratinocytes after gamma irradiation. J Cell Biochem. 2005, 95 (3): 620-631. 10.1002/jcb.20394.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33 (20): e179-10.1093/nar/gni178.
Wang B, Howel P, Bruheim S, Ju J, Owen LB, Fodstad O, Xi Y: Systematic evaluation of three microRNA profiling platforms: microarray, beads array, and quantitative real-time PCR array. PLoS One. 2011, 6 (2): e17167-10.1371/journal.pone.0017167.
Gougelet A, Pissaloux D, Besse A, Perez J, Duc A, Dutour A, Blay JY, Alberti L: Micro-RNA profiles in osteosarcoma as a predictive tool for ifosfamide response. Int J Cancer. 2011, 129 (3): 680-690. 10.1002/ijc.25715.
Pradervand S, Weber J, Thomas J, Bueno M, Wirapati P, Lefort K, Dotto GP, Harshman K: Impact of normalization on miRNA microarray expression profiling. RNA. 2009, 15 (3): 493-501. 10.1261/rna.1295509.
Hutvagner G, McLachlan J, Pasquinelli AE, Balint E, Tuschl T, Zamore PD: A cellular function for the RNA-interference enzyme Dicer in the maturation of the let-7 small temporal RNA. Science. 2001, 293 (5531): 834-838. 10.1126/science.1062961.
Chaudhry MA, Omaruddin RA, Kreger B, de Toledo SM, Azzam EI: Micro RNA responses to chronic or acute exposures to low dose ionizing radiation. Mol Biol Rep. 2012, 39 (7): 7549-7558. 10.1007/s11033-012-1589-9.
Lynam-Lennon N, Reynolds JV, Marignol L, Sheils OM, Pidgeon GP, Maher SG: MicroRNA-31 modulates tumour sensitivity to radiation in oesophageal adenocarcinoma. J Mol Med (Berl). 2012, 90 (12): 1449-1458. 10.1007/s00109-012-0924-x.
Mueller AC, Sun D, Dutta A: The miR-99 family regulates the DNA damage response through its target SNF2H. Oncogene. 2012
Baghdoyan S, Lamartine J, Castel D, Pitaval A, Roupioz Y, Franco N, Duarte M, Martin MT, Gidrol X: Id2 reverses cell cycle arrest induced by {gamma}-irradiation in human HaCaT keratinocytes. J Biol Chem. 2005, 280 (16): 15836-15841. 10.1074/jbc.M414216200.
Luo H, Zou J, Dong Z, Zeng Q, Wu D, Liu L: Up-regulated miR-17 promotes cell proliferation, tumour growth and cell cycle progression by targeting the RND3 tumour suppressor gene in colorectal carcinoma. Biochem J. 2012, 442 (2): 311-321. 10.1042/BJ20111517.
Elyakim E, Sitbon E, Faerman A, Tabak S, Montia E, Belanis L, Dov A, Marcusson EG, Bennett CF, Chajut A: hsa-miR-191 is a candidate oncogene target for hepatocellular carcinoma therapy. Cancer Res. 2010, 70 (20): 8077-8087. 10.1158/0008-5472.CAN-10-1313.
Dickey JS, Zemp FJ, Martin OA, Kovalchuk O: The role of miRNA in the direct and indirect effects of ionizing radiation. Radiat Environ Biophys. 2011, 50 (4): 491-499. 10.1007/s00411-011-0386-5.
Zhou BR, Xu Y, Permatasari F, Liu WL, Li W, Guo XF, Huang QH, Guo Z, Luo D: Characterization of the miRNA profile in UVB-irradiated normal human keratinocytes. Exp Dermatol. 2012, 21 (4): 317-319. 10.1111/j.1600-0625.2012.01465.x.
Song HJ, Cho CK, Yoo SY, Park KS, Lee YS: Increased induction of Ca2+−mediated differentiation by gamma ray is mediated by endogenous activation of the protein kinase C signaling pathways in mouse epidermal cells. Int J Radiat Oncol Biol Phys. 1998, 41 (4): 897-904. 10.1016/S0360-3016(98)00143-6.
Yugawa T, Handa K, Narisawa-Saito M, Ohno S, Fujita M, Kiyono T: Regulation of Notch1 gene expression by p53 in epithelial cells. Mol Cell Biol. 2007, 27 (10): 3732-3742. 10.1128/MCB.02119-06.
Hildebrand J, Rutze M, Walz N, Gallinat S, Wenck H, Deppert W, Grundhoff A, Knott A: A comprehensive analysis of microRNA expression during human keratinocyte differentiation in vitro and in vivo. J Invest Dermatol. 2011, 131 (1): 20-29. 10.1038/jid.2010.268.
Choi YJ, Kim SY, Oh JM, Juhnn YS: Stimulatory heterotrimeric G protein augments gamma ray-induced apoptosis by up-regulation of Bak expression via CREB and AP-1 in H1299 human lung cancer cells. Exp Mol Med. 2009, 41 (8): 592-600. 10.3858/emm.2009.41.8.065.
Cho EA, Juhnn YS: The cAMP signaling system inhibits the repair of gamma-ray-induced DNA damage by promoting Epac1-mediated proteasomal degradation of XRCC1 protein in human lung cancer cells. Biochem Biophys Res Commun. 2012, 422 (2): 256-262. 10.1016/j.bbrc.2012.04.139.
Fazel R, Krumholz HM, Wang Y, Ross JS, Chen J, Ting HH, Shah ND, Nasir K, Einstein AJ, Nallamothu BK: Exposure to low-dose ionizing radiation from medical imaging procedures. N Engl J Med. 2009, 361 (9): 849-857. 10.1056/NEJMoa0901249.
Templin T, Amundson SA, Brenner DJ, Smilenov LB: Whole mouse blood microRNA as biomarkers for exposure to gamma-rays and (56)Fe ion. Int J Radiat Biol. 2011, 87 (7): 653-662. 10.3109/09553002.2010.549537.
Abou-El-Ardat K, Monsieurs P, Anastasov N, Atkinson M, Derradji H, De Meyer T, Bekaert S, Van Criekinge W, Baatout S: Low dose irradiation of thyroid cells reveals a unique transcriptomic and epigenetic signature in RET/PTC-positive cells. Mutat Res. 2012, 731 (1–2): 27-40.
Neumaier T, Swenson J, Pham C, Polyzos A, Lo AT, Yang P, Dyball J, Asaithamby A, Chen DJ, Bissell MJ: Evidence for formation of DNA repair centers and dose–response nonlinearity in human cells. Proc Natl Acad Sci USA. 2012, 109 (2): 443-448. 10.1073/pnas.1117849108.
Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinforma. 2009, 10: 48-10.1186/1471-2105-10-48.
Meijerink J, Mandigers C, van de Locht L, Tonnissen E, Goodsaid F, Raemaekers J: A novel method to compensate for different amplification efficiencies between patient DNA samples in quantitative real-time PCR. J Mol Diagn. 2001, 3 (2): 55-61. 10.1016/S1525-1578(10)60652-6.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55 (4): 611-622. 10.1373/clinchem.2008.112797.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta CT)) method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Alcaraz C, De Diego M, Pastor MJ, Escribano JM: Comparison of a radioimmunoprecipitation assay to immunoblotting and ELISA for detection of antibody to African swine fever virus. J Vet Dign Invest. 1990, 2 (3): 191-196. 10.1177/104063879000200307.
Acknowledgements
We thank Dr Clavel (Service d’Urologie, Clinique du Tonkin, Villeurbanne, France) for providing skin biopsies, Severine Croze and Dr Joël Lachuer (ProfileXpert-LCMT facility, Lyon, France) for their help in microfluidic QPCR, Dr Walid Rachidi (University of Grenoble) for his help in low-dose irradiation. We also thank Dr Ingrid Masse and Dr Odile Berthier-Vergnes for critical reading of the manuscript. This study was supported by grants from the French National Research Agency (ANR LODORA), Electricité de France (Commission Scientifique de Radioprotection) and by La Ligue contre le Cancer (Comité de Saône et Loire).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contribution
NJT performed experiments, analyzed data and wrote the manuscript. JV performed statistical analysis and wrote the manuscript. OG designed experiments and analyzed data. JL designed experiments, analyzed data and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
12864_2012_4858_MOESM1_ESM.xls
Additional file 1: Full data of miRNAs analyzed in our study. The Cq and relative expression for all analyzed miRNAs are given for the 3 conditions (Proliferating cells 6 Gy, Differentiated Cells 6 Gy, Proliferating cells 10 mGy). (XLS 105 KB)
12864_2012_4858_MOESM2_ESM.pdf
Additional file 2: Viability assays after combined up-regulation of two miRNAs in proliferating keratinocytes. The relative 6Gy/0Gy viability (reported to the pre-miR negative control transfection) is indicated for each combination of two pre-miRNAs transfection at 24 h, 48 h and 72 h post-irradiation (n=4 ; error bars show SD ; bilateral paired t-test : ***p≤0.001). (PDF 58 KB)
12864_2012_4858_MOESM3_ESM.pdf
Additional file 3: Venn diagram showing the predicted targets in common between the 3 miRNAs miR-191-5p, miR200b-3p and miR-331-3p. (PDF 16 KB)
12864_2012_4858_MOESM4_ESM.xls
Additional file 4: List of the predicted targets for the microRNAs miR-191-5p, miR200b-3p and miR-331-3p and list of the targets in common between these 3 miRNAs. (XLS 288 KB)
12864_2012_4858_MOESM5_ESM.docx
Additional file 5: Primers Table. Primers specificity was assayed by PrimerBLAST or Primer3. Primers were designed to amplify all mRNA variants if exist. Primers were purified by SePOP. (DOCX 64 KB)
12864_2012_4858_MOESM6_ESM.xls
Additional file 6: Minimum information for publication of quantitative real-time PCR experiments (MIQE) checklist. (XLS 40 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Joly-Tonetti, N., Viñuelas, J., Gandrillon, O. et al. Differential miRNA expression profiles in proliferating or differentiated keratinocytes in response to gamma irradiation. BMC Genomics 14, 184 (2013). https://doi.org/10.1186/1471-2164-14-184
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1471-2164-14-184