Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2012; 7(2): e31028.
Published online 2012 February 8. doi:  10.1371/journal.pone.0031028
PMCID: PMC3275573

Unraveling the Global microRNAome Responses to Ionizing Radiation in Human Embryonic Stem Cells

Majlinda Lako, Editor


MicroRNAs (miRNA) comprise a group of short ribonucleic acid molecules implicated in regulation of key biological processes and functions at the post-transcriptional level. Ionizing radiation (IR) causes DNA damage and generally triggers cellular stress response. However, the role of miRNAs in IR-induced response in human embryonic stem cells (hESC) has not been defined yet. Here, by using system biology approaches, we show for the first time, that miRNAome undergoes global alterations in hESC (H1 and H9 lines) after IR. Interrogation of expression levels of 1,090 miRNA species in irradiated hESC showed statistically significant changes in 54 genes following 1 Gy of X-ray exposures; global miRNAome alterations were found to be highly temporally and cell line - dependent in hESC. Time-course studies showed that the 16 hr miRNAome radiation response of hESC is much more robust compared to 2 hr-response signature (only eight genes), and may be involved in regulating the cell cycle. Quantitative real-time PCR performed on some miRNA species confirms the robustness of our miRNA microarray platform. Positive regulation of differentiation-, cell cycle-, ion transport- and endomembrane system-related processes were predicted to be negatively affected by miRNAome changes in irradiated hESC. Our findings reveal a fundamental role of miRNAome in modulating the radiation response, and identify novel molecular targets of radiation in hESC.


MicroRNAs (miRNAs) are considered to constitute a class of small noncoding RNAs vitally involved in regulation of gene expression and signal transduction [1]. Published data suggest miRNAs act as post-transcriptional regulators that may control the expression of about 60% of human genes [2] by means of messenger RNA (mRNA) decay and/or translational repression [3]. Deregulated expression of miRNAs was shown to underlie many diseased states; and miRNAome patterns prove to be highly specific, in many cases accurately reflecting the stage and prognosis of disease [4], [5]. Increasing body of evidence implies changes in miRNA expression profiles after genotoxic stress exposures, including ionizing radiation (IR) [6], [7], [8]. Previous studies analyzed miRNA expression following IR either in fully differentiated human normal somatic cells/artificial 3D tissues [9], [10], [11], [12], [13], [14] or in cancerous cells [15], [16], [17], [18], [19], [20]. Human embryonic stem cells (hESC) represent pluripotent cells with unique capabilities to differentiate into virtually all cell types of a human body. Human stem cells are believed to be endowed with mechanisms to ensure the superior genome fidelity that are just begun to be explored systematically [21], [22]. Exposures to medical, environmental or accidental sources of radiation can challenge the human well-being in the exposed populations, including pregnant women. We, and others, recently studied the IR response of hESC extensively [23], [24], [25], [26], [27], [28], [29]. However, none of the published papers focused on elucidating the miRNA signatures in hESC after IR exposures. Moreover, since the human genome is estimated to encode more than 1,000 of distinct miRNA species, it is imperative to employ high-throughput methodologies to fully analyze the role of miRNA species in biological processes. To the best of our knowledge, only one recent report analyzed miRNAome changes in irradiated human cells at a whole genome-wide level [20].

In the present work, we aimed to study the effects of IR on global miRNAome in hESC, to perform prediction of the biological processes/themes affected by miRNAome changes, and to develop a miRNA-based gene expression signature specific for irradiated, but not for non-irradiated, hESC. For the first time, we found that miRNAome undergoes genome-wide alterations in hESC after IR. All miRNAs published in the Sanger miRBase release version 15.0 ( were interrogated. The expression levels of 1,090 miRNA species in irradiated hESC showed statistically significant changes in 54 genes following 1 Gy of X-ray exposures (p<0.05). We found that many miRNA species were modulated in a cell line-specific manner in hESC after IR, with H1 cell line being more radiation responsive than H9. Hierarchical clustering, class prediction and Gene Ontology analysis were performed to characterize hESC miRNAome response to radiation in more detail. Positive regulation of differentiation-, cell cycle-, ion transport- and endomembrane system-related processes were predicted as being negatively affected by alterations in miRNAome in irradiated hESC.

Materials and Methods

1. Cell culture and treatments

Human ESCs (H1 and H9 cell lines, WiCell, Madison, WI, passage 35–40; BG01V line, ATCC, Manassas, VA) were routinely cultured in mTeSR-1 medium (Stemcell Technologies, Vancouver, Canada) using cell culture vessels coated with BD Matrigel hESC-qualified Matrix (BD Biosciences, San Jose, CA) at 37°C and 5% CO2. Cells were grown following supplier's protocol and as described in [23], [24]. Cell cultures were exposed to 1 Gy of X-ray irradiation using X-RAD 320 Biological Irradiator unit (Precision X-Ray, Inc., North Branford, CT; dose rate about 1 Gy/min; 320 kV, 12.5 mA); then the cells were returned to CO2 incubator and harvested at 2 h and 16 h post-irradiation for analysis. The mock-irradiated cells for each time-point were used as a control. For functional analysis of hsa-mir-575, both gain- and loss-of-function studies were carried out. For overexpression studies, specific hsa-mir-575 mirVana miRNA mimic transfection experiments were performed with H1 hESC cultures according to manufacturer's protocol (Applied Biosystems, Carlsbad, CA). MirVana miRNA inhibitor (Applied Biosystems, Carlsbad, CA) was used to knock-down hsa-mir-575 expression in a separate set of experiments. Transfection studies, in parallel with negative control miRNA inhibitor experiments, were carried out with Stemfect RNA Transfection kit per vendor's protocol (Stemgent, Cambridge, MA). To examine the viability of hESC in colonies upon hsa-mir-575 overexpression, cells were incubated at 37°C for 1 h with Hoechst 33342 (8 µg/ml; Molecular Probes, Eugene, OR) and propidium iodide (PI, 20 µg/ml; Sigma, St. Louis, MO) as in [25]. Cell colonies were visualized using a fluorescence microscope (Axioplan 2, Zeiss, Thornwood, NY) equipped with a fluorescent light source.

2. RNA sample preparation, probe labeling and DNA microarray procedure

The extraction of total RNA was performed with miRNeasy kit (Qiagen, Valencia, CA) per manufacturers' instructions. The amount and quality of RNA samples were assessed on the Agilent 2100 Bioanalyzer with RNA 6000 Nano Reagents and Supplies (Agilent, Santa Clara, CA). Subsequently, 0.5 µg of total RNA was used in each reaction to generate labeled samples with miRCURY LNA microRNA Hy5 Power labeling kit (Exiqon, Woburn, MA). MiRNA spike control was added to the RNA samples prior to the labeling reactions following the manufacturer's protocol. The labeled targets corresponding either to experimental or control samples were separately hybridized to “3D-Gene” oligo microarrays provided by the manufacturer's (Toray Industries Inc., Tokyo, Japan; miRBase release 15.0) containing 1,090-elements spotted in duplicate using Takara Hybridization chambers (Takara Bio, Inc., Japan). Protocols for microarray hybridization and washing were as provided by manufacturer. Hybridized DNA microarrays were scanned on an Axon GenePix DNA microarray scanner (Molecular Devices, Inc., Sunnyvale, CA), and TIFF images were subsequently generated for further analysis. The sample labeling, hybridization, “3D-Gene” array washing and scanning were conducted by Toray Industries Inc. (Tokyo, Japan).

3. Microarray data analysis

Analyses were performed using BRB-ArrayTools Version 4.2.0 developed by Dr. Richard Simon and BRB-ArrayTools Development Team (Biometric Research Branch, National Cancer Institute, NIH). Time-matched irradiated versus sham-irradiated samples were used to determine the radiation–responsive miRNA species for each data point. The data from the microarray was collected and analyzed in accordance to the Minimum Information About a Microarray Experiment (MIAME) guidelines. MIAME-compliant raw data for this series of experiments have been deposited in the ArrayExpress database maintained by the European Bioinformatics Institute (accession no. E-MEXP-3366). Differentially expressed miRNA genes were identified using random-variance t-test [30] and as in [31], [32]. Changes in gene expression were considered statistically significant if the p values for corresponding genes were less than 0.05.

4. Prediction of miRNA targets

Candidate miRNA species with a microarray intensity signal ≥100 and p-value≤0.05, identified as being differentially expressed after IR exposures with BRB-ArrayTools, were chosen for target prediction analysis. MiRanda database was used to predict miRNA targets [33]. Targets were input to the Database for Annotation, Visualization and Integrated Discovery (DAVID) [34], [35] themes analysis. Biological themes/processes were determined using the Functional Annotation Clustering feature using p-value less than 0.05.

5. Quantitative real-time PCR

Total RNA samples used for the microarray analysis were reverse transcribed to cDNA template using RT-specific primers (Applied Biosystems, Carlsbad, CA). Then, TaqMan miRNA assays were performed in triplicate using miRNA Cells-to-Ct kit reagents (Applied Biosystems) per manufacturer's instructions. The U6 snRNA was used as an internal control. RT-PCR was performed on iCycler iQ (Bio-Rad, Inc., Hercules, CA) in 20-µl reactions by using TaqMan Assay-on-Demand primers/probe sets (Applied Biosystems) for the following miRNA genes: hsa-miR-302b, hsa-miR-575, hsa-miR-1274b, hsa-miR-1915 and hsa-miR-1973. In a separate set of experiments, H1 hESC cultures transfected either with hsa-mir-575 mirVana miRNA mimic or miRNA inhibitor were subjected to lysis 5 days post-transfection. TaqMan gene expression studies for POU5F1 and SOX2 were carried out in triplicate with Cells-to-Ct kit reagents (Applied Biosystems), and 18S RNA expression was used as a reference. Quantitative RT-PCR data were analyzed as in [32].

Results and Discussion

To identify radiation inducible miRNA species and characterize changes in miRNAome after IR in hESC, we exposed both H1 and H9 cell lines to 1 Gy of X-rays. These cell lines are most extensively studied among hESC. Our miRNA microarray data analysis revealed that 53 and six miRNAs, were differentially expressed in irradiated H1 and H9 cells respectively, compared to corresponding sham-irradiated hESC cultures (Figure S1). These results indicate that miRNA expression profiles of hESC after IR exposures are largely cell line-dependent.

Further bioinformatics analyses revealed that only seven miRNA genes were differentially expressed in H1 at the “early” (2 hr) time point post-IR (p<0.05) (Figure 1); among them, only four miRNA species showed more than 1.5 – fold induction compared to sham-irradiated cells (hsa-miR-15b, hsa-miR-1274b, hsa-miR-302b and hsa-miR-1973) (Table 1). Gene expression studies showed that IR-induced alterations in miRNAome in H9 cell line (2 hr post-IR) involve the up-regulation of only two miRNA genes, namely, hsa-miR-1973 and hsa-miR-92a (Table 2). In both cell lines, the level of upregulation of miRNA species at this time point was less or equal to 2-fold over corresponding non-irradiated baseline values. No significant down-regulation of miRNA was observed at 2 hr following IR exposures. Interestingly, lack of statistically proven repression of any miRNA species at relatively early 4 hr post 1.25 Gy IR exposure of human blood was reported recently [13].

Figure 1
Venn diagrams of significantly upregulated miRNA species in human embryonic stem cells after IR exposures (p<0.05).
Table 1
Differentially expressed miRNA species in H1 hESC at 2 hr post 1 Gy of X-ray exposures (p<0.05).
Table 2
Differentially expressed miRNA species in H9 hESC at 2 hr post 1 Gy of X-ray exposures (p<0.05).

On the contrary, at 16 hr timepoint after IR, we observed both statistically significant up-regulation (Table 3) and repression (Table 4) of select miRNA species in H1 hESC. Irradiated H1 cells showed up-regulation of 42 miRNA genes, and down-regulation of only eight miRNAs at this time post 1 Gy IR exposure (p<0.05). A much higher magnitude of changes in the level of expression of IR-modulated miRNAs was observed as part of the “late” hESC response to IR (Table 3); 33 miRNA species were induced more than 2 – fold in irradiated H1 compared to control sham-exposed cell cultures. In contrast, we found that only five miRNAs were significantly overexpressed in H9 cell line, with a single miRNA species, hsa-miR-575, being overexpressed more than 2 fold (p<0.05) (Table 5).

Table 3
Up-regulated miRNA species in H1 hESC at 16 hr post 1 Gy of IR exposures (p<0.05).
Table 4
Down-regulated miRNA species in H1 hESC at 16 hr post 1 Gy of IR exposures (p<0.05).
Table 5
Differentially expressed miRNA species in H9 hESC at 16 hr post 1 Gy of X-ray exposures (p<0.05).

To identify the factors that affect the patterns of microRNAome changes after IR exposures, all experimental samples were subjected to hierarchical clustering (Figure 2). This analysis grouped hESC samples in accordance to IR exposure conditions, with the most dramatic changes in miRNAome occurring as a result of dynamic of cellular response to radiation. The timing of post-IR exposure response, and not a cell type – specific changes, was the most prominent determinant of the clustering of samples analyzed.

Figure 2
Heatmap representation of differentially expressed miRNA species in human embryonic stem cells after IR exposures (p<0.05).

Our recent study demonstrated no statistically significant downregulation of mRNA gene transcripts after 1 Gy of IR exposure in H9 cells [24]. Intriguingly, in the present study we observed a few miRNAs that were upregulated in H9 under these same conditions (Table 2 and Table 5). Since miRNAs are considered to be one of the powerful epigenetic means to silence gene expression, it may seem puzzling at first glance. However, it is thought that in mammals the primary mode of miRNA-mediated gene repression is through inhibition of translation, not direct mRNA target degradation [36]. Therefore, it could explain the lack of downregulation of any protein-coding transcripts observed in our previous work [24]. The repression of genes involved in hESC differentiation after IR exposures may explain the maintenance of pluripotentiality in surviving irradiated hESC that we reported before [25]. Changes in cell cycle profile of IR-exposed hESC [24] may also be the consequence of miRNAome alterations. For example, one of the key genes responsible for cell cycle progression and about 26% of all phosphorylation events in hESC is CDK1/2 [37], which we found to be the major predicted target of the most highly upregulated miRNA (hsa-miR-575) in our present study. Interestingly, hsa-miR-575 is known to be overexpressed in some aggressive human cancers [38], and is downregulated upon human leukemic HL-60 cell differentiation induced by 4-hydroxynonenal, a product of lipid peroxidation [39]. We sought to investigate the role of hsa-miR-575 in hESC in more detail. Our studies aimed to mimic the overexpression of this miRNA species showed no evidence for an increase in cell killing in H1 hESC cultures (Figure S2). Moreover, we did not observe signs of differentiation in these cultures since continuous expression of markers of pluripotency, such as POU5F1 and SOX2, was evident at least 5 days post-transfection (Table S1). Interestingly, upon treatment of irradiated H1 cells with hsa-miR-575 inhibitor as part of our loss-of-function studies, we observed substantial down-regulation of POU5F1 and upregulation of SOX2 (Table S1). Both genes are known to play an important role in maintenance of pluripotency in undifferentiated hESC, and major alterations in their expression levels may result in differentiation [40], [41]. Our results may support the assumption that upregulation of hsa-miR-575 in irradiated hESC serves to prevent differentiation of these cells, but additional studies into the exact function of hsa-miR-575 are needed to be performed in a future.

In general, the majority of differentially expressed miRNAs that we identified in our global screen of irradiated hESC are not well-characterized yet. However, the involvement of some IR-responsive miRNAs in biological processes was already reported. For example, the most highly induced miRNA in H1 cells at 2 hr post-1 Gy of IR exposures, hsa-miR-15b, was shown to sensitize human cancer cells to apoptosis [42]. This might be relevant for cell fate choice in hESC as well, since massive apoptosis was observed in irradiated hESC after 1 Gy of X-ray exposures [25]. In addition, the up-regulation of this miRNA species could influence many other aspects of cellular homeostasis: miR-15b is overexpressed in bronchioalveolar lung stem cells [43], is a direct transcriptional target of E2F1 modulating cell cycle progression [44], and is known to affect cell metabolism through decreasing ATP levels [45].

Hsa-mir-302b was shown to be highly expressed in hESC [46] and is directly responsible for “stemness” characteristics in human cells being part of stem cell-enriched mir-302-367 cluster [47], [48]. Mir-302 species were determined to be highly integrated with Oct4/Sox2 transcriptional machinery in hESC being partly responsible for unique abbreviated cell cycle profile in hESCs [49]. The modest upregulation of this miRNA species following IR exposures may prevent hESC from undergoing spontaneous differentiation. Other IR-modulated miRNA species, such as hsa-mir-20a, is a component of yet another stem cell-enriched family, namely, mir-17-92, which is implicated in modulating E2F activity on cell cycle progression [50] and in the repression of execution of cellular senescence program [51], [52]. Hsa-mir-149* also regulates E2F1 activity and exerts its function through induction of apoptosis [53]. One more hESC highly expressed miRNA, hsa-mir-1909, was implicated in hESC physiology by targeting Notch [54].

A large number of overexpressed IR-responsive miRNAs that we identified in our work were found to be deregulated in human cancers, such as hsa-mir-513 [55], hsa-mir-744 [56], hsa-mir-92a [57], [58], hsa-mir-1228* [59], hsa-mir-671-5p [60], hsa-mir-638 [38], hsa-mir-370 [61], and hsa-mir-675 [62]. Hsa-mir-663 is involved in inflammatory conditions and suppression of cell proliferation [63], [64] which can underlie hESC response to IR exposures. Interestingly, two of significantly upregulated miRNAs after IR exposures in hESC, such as hsa-mir-1275 and hsa-mir-494, were recently found to localize to mitochondria [65]. Given the importance of mitochondrial compartment in the maintenance of stem cell homeostasis and responses to different stress factors [66], the functional validation of targets of these miRNA species might prove to be instrumental in unraveling the underpinnings of hESC unique biological capabilities.

We observed several miRNA species that were downregulated following irradiation (Table 4). Repression of hsa-mir-93 could block cell differentiation after IR exposures in part by relieving STAT3 [67]. Other biological effect of decreased expression of hsa-mir-93 at 16 hr post-IR in H1 could be a decreased cell survival [68], probably as a result of modulation of its targets such as E2F1 and CDKN1A [69]. Interestingly, CDKN1A was found to be a target for other down-regulated miRNA species we identified in our present study, namely, hsa-mir-302a, hsa-mir-302b, and hsa-mir-302d [70]. With CDKN1A being one of the key radiation-responsive genes in human cells, the possible involvement of miRNAome alterations in its regulation after IR definitely merits further investigations.

We attempted to compare our data on global miRNAome alterations in irradiated hESC with published datasets. However, since the knowledge about changes in miRNA levels in human cells following irradiation is very limited to date, we were able to identify only a few miRNA species in our subset of IR-modulated miRNAs which were reported before as responsive to IR exposures. Among them, we found hsa-mir-20a, shown to be IR-responsive in fully differentiated human endothelial cells [10] and implicated in regulation of IR-induced premature senescence [52]; hsa-mir-24, a member of the miR-23b cluster that interferes with TGF-β expression and shown to be up-regulated after IR [20]. Mir-24 is implicated in halting cell cycle progression and inhibition of apoptosis [71], and impairment of DDR and senescence program execution [72], [73]. In addition, hsa-mir-744 and hsa-mir-17 were also found to be deregulated by IR exposures by others [20]. Our studies with BG01V line of hESC, considered being karyotypically abnormal, identified hsa-mir-1915 and hsa-mir-1274b as differentially expressed at 2 hr and 16 hr post-IR, respectively (Table S2). We believe that further experiments will examine the relative contribution of cell type specificity to a repertoire of IR-modulated miRNAs at a whole-genome level in various human cells.

We analyzed the possible targets of differentially expressed miRNAs by running MiRanda algorithm search. The results are presented in Tables S3, S4, S5, S6, S7. At least some of the highly scored targets may represent an integrated network with miRNAome and act in concert in irradiated hESCs; for example, LIN9, a putative target of hsa-mir-1973, regulates MYB which is by itself is predicted to be a target of hsa-mir-15b (Table S3), and crucially affects the cell cycle machinery through the cyclins and CDK1 [74]. Hence, modulation of regulators of key biological processes by alterations in miRNAome may represent a powerful strategy used by IR-exposed hESC to cope with genotoxic stress.

Despite numerous recent advances, the human miRNAome still remains mainly unexplored regarding the physiological function of specific miRNA species within cells. It should be noted that deciphering the function of individual miRNAs is challenging. There are many families comprising microRNAs differing only in one-two nucleotides which make a functional assignment to these redundant genes a complicated task. Each miRNA is thought to target numerous putative transcripts that have non-related functions, thereby the experimental validation of specific targets for miRNAs is still in its infancy. Therefore, determining which biological processes might be prime candidates for miRNA-mediated regulation of gene expression after IR exposures might prove to be more informative than querying individual miRNA-mRNA pairs. To this end, we performed Gene Ontology (GO) analysis of biological themes/processes. The aim of the GO analysis was to predict the effect of overexpressed miRNAs on cellular functions. The downregulated genes predicted to be targeted by the radiation-induced miRNAs, were uploaded to the DAVID GO database, and gene-set enrichment analysis was carried out in the biologic process/molecular function/cellular component categories. The results of the GO analysis of the biological themes overrepresented in the pools of the genes predicted to be targeted by miRNAs are shown in Table 6 and Table 7. Our data indicate that IR exposures of H1 and H9 hESCs could result in an increased miRNA control of genes involved in positive regulation of cell differentiation, cell projection, transcription activation, alternative splicing, cell death and cell cycle regulation.

Table 6
Gene Ontology analysis of affected biological themes based on predicted miRNA targets in H1 cell line.
Table 7
Gene Ontology analysis of affected biological themes based on predicted miRNA targets in H9 cell line.

To partially validate miRNAome changes we obtained in our study, we compared data on miRNAs expression with two independent techniques, that is, miRNA microarray platform and quantitative RT-PCR (Table 8). In general, we observed a good concordance between these datasets, thus confirming the robustness of our approach.

Table 8
Verification of miRNA microarray data with Taqman qRT-PCR. Shown are means and standard errors for corresponding values obtained by two techniques.

In summary, we found that microRNAome of human embryonic stem cells undergoes global alterations following ionizing radiation exposures. We showed that the gene expression signature for global miRNAome alterations at a “late”, 16 hr time point is significantly different in comparison to that of “early” response. The changes that we observed in the levels of miRNA species constituting the whole human microRNAome are substantially cell line-dependent, with H1 cell line of hESC being more responsive to IR than H9 line. Surprisingly, there were just a few H9-specific IR-modulated miRNA species compared to H1. Since these two cell lines are most thoroughly analyzed to date among all published hESC lines, it might be of interest to probe miRNAome of these lines under more various experimental settings in a future. We observed consistent patterns of miRNA response to IR exposures, including several miRNA species that were systematically differentially expressed at both time intervals. But the majority of IR-modulated miRNA showed a time-dependent response in both H1 and H9 hESCs, with a more robust response occurring at a late time point. It should be noted that the degree of target repression imposed by miRNAs is probably quantitatively modest since most of specific miRNA species endogenous targets is usually downregulated by less than 50% [75]. Given these considerations, most proteins comprising the human proteome presumably remain effective over this degree of inhibition. Therefore, future studies will focus on how alterations in miRNAome affect the proteome of hESC following IR-exposures, and how it translates to the ultimate cell fate of these unique human cells that show great promise in both disease modeling and cell regenerative therapy approaches.

The strength of our study is based upon the comprehensive coverage of the whole human miRNAome in an experimental platform that we employed, in profiling the changes in miRNAs in two different hESC lines, and in analyzing the dynamics of response of hESC to IR exposures at the level of miRNAome. We identified a number of radiation-responsive miRNA species with as yet unknown functions that provide a broad avenue to future research in this area. We attempted to outline some of hypothesis that can foster the investigation of functional relevance of miRNAome alterations in hESC. Our study reveals new insights into how hESC respond to genotoxic stress, in particular, to IR exposures resulting in global alterations in microRNAome in these cells. The findings may contribute to improved understanding of the biology of hESC, and the mechanisms underlying enhanced genomic maintenance in hESC.

Supporting Information

Figure S1

Venn diagram of a total number of differentially expressed miRNA species in human embryonic stem cells after IR exposures (p<0.05).


Figure S2

H1 hESC culture staining for viability upon hsa-mir-575 changes in gene expression studies. Cell cultures were stained with Hoechst 33342 (shown in blue) and propidium iodide (in red). A – control, 0 Gy; B – mock transfection; C – hsa-mir-575 mimic transfection; D – 1 Gy, 24 hrs post- IR.


Table S1

Taqman qRT-PCR on cultured H1 hESCs. Shown are means and standard errors for gene expression changes obtained following hsa-mir-575 expression modulation in comparison to mock-treated cell cultures.


Table S2

Taqman qRT-PCR on cultured BG01V hESCs. Shown are means and standard errors for miRNA expression changes obtained for indicated timepoints after irradiation in comparison to mock-treated cell cultures.


Table S3

Up-regulated (>1.5 - fold) miRNA genes (1 Gy, 2 hr) in H1 as determined by microarray analysis (p<0.05).


Table S4

Selection of top 10 up-regulated miRNA genes (1 Gy, 16 hr) in H1 as determined by microarray analysis (p<0.05).


Table S5

Down-regulated (>1.5 - fold) miRNA genes (1 Gy, 16 hr) in H1 (p<0.05).


Table S6

Up-regulated (>1.5 - fold) miRNA genes (1 Gy, 2 hr) in H9 (p<0.05).


Table S7

List of up-regulated (>1.5 - fold) miRNA genes (1 Gy, 16 hr) in H9 (p<0.05).



We wish to thank William DeGraff, Radiation Biology Branch, CCR, National Cancer Institute, for his invaluable help with cell culture X-ray irradiations.


Competing Interests: The authors have declared that no competing interests exist.

Funding: This research was supported in part by the Intramural Research Program of the National Institutes of Health, Clinical Center. No additional external funding received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Inui M, Martello G, Piccolo S. MicroRNA control of signal transduction. Nat Rev Mol Cell Biol. 2010;11:252–263. [PubMed]
2. Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009;19:92–105. [PubMed]
3. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–233. [PubMed]
4. Croce CM. Causes and consequences of microRNA dysregulation in cancer. Nat Rev Genet. 2009;10:704–714. [PMC free article] [PubMed]
5. Latronico MV, Condorelli G. MicroRNAs and cardiac pathology. Nat Rev Cardiol. 2009;6:419–429. [PubMed]
6. Josson S, Sung SY, Lao K, Chung LW, Johnstone PA. Radiation modulation of microRNA in prostate cancer cell lines. Prostate. 2008;68:1599–1606. [PMC free article] [PubMed]
7. 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:824–834. [PubMed]
8. Chaudhry MA. Real-time PCR analysis of micro-RNA expression in ionizing radiation-treated cells. Cancer Biother Radiopharm. 2009;24:49–56. [PubMed]
9. Simone NL, Soule BP, Ly D, Saleh AD, Savage JE, et al. Ionizing radiation-induced oxidative stress alters miRNA expression. PLoS One. 2009;4:e6377. [PMC free article] [PubMed]
10. Wagner-Ecker M, Schwager C, Wirkner U, Abdollahi A, Huber PE. MicroRNA expression after ionizing radiation in human endothelial cells. Radiat Oncol. 2010;5:25. [PMC free article] [PubMed]
11. Kovalchuk O, Zemp FJ, Filkowski JN, Altamirano AM, Dickey JS, et al. microRNAome changes in bystander three-dimensional human tissue models suggest priming of apoptotic pathways. Carcinogenesis. 2010;31:1882–1888. [PMC free article] [PubMed]
12. Nikiforova MN, Gandi M, Kelly L, Nikiforov YE. MicroRNA dysregulation in human thyroid cells following exposure to ionizing radiation. Thyroid. 2011;21:261–266. [PMC free article] [PubMed]
13. Templin T, Paul S, Amundson SA, Young EF, Barker CA, et al. Radiation-Induced Micro-RNA Expression Changes in Peripheral Blood Cells of Radiotherapy Patients. Int J Radiat Oncol Biol Phys 2011 [PMC free article] [PubMed]
14. Vincenti S, Brillante N, Lanza V, Bozzoni I, Presutti C, et al. HUVEC Respond to Radiation by Inducing the Expression of Pro-angiogenic MicroRNAs. Radiat Res 2011 [PubMed]
15. Cha HJ, Shin S, Yoo H, Lee EM, Bae S, et al. Identification of ionizing radiation-responsive microRNAs in the IM9 human B lymphoblastic cell line. Int J Oncol. 2009;34:1661–1668. [PubMed]
16. Shin S, Cha HJ, Lee EM, Jung JH, Lee SJ, et al. MicroRNAs are significantly influenced by p53 and radiation in HCT116 human colon carcinoma cells. Int J Oncol. 2009;34:1645–1652. [PubMed]
17. Shin S, Cha HJ, Lee EM, Lee SJ, Seo SK, et al. Alteration of miRNA profiles by ionizing radiation in A549 human non-small cell lung cancer cells. Int J Oncol. 2009;35:81–86. [PubMed]
18. Cha HJ, Seong KM, Bae S, Jung JH, Kim CS, et al. Identification of specific microRNAs responding to low and high dose gamma-irradiation in the human lymphoblast line IM9. Oncol Rep. 2009;22:863–868. [PubMed]
19. Chaudhry MA, Sachdeva H, Omaruddin RA. Radiation-induced micro-RNA modulation in glioblastoma cells differing in DNA-repair pathways. DNA Cell Biol. 2010;29:553–561. [PubMed]
20. Niemoeller OM, Niyazi M, Corradini S, Zehentmayr F, Li M, et al. MicroRNA expression profiles in human cancer cells after ionizing radiation. Radiat Oncol. 2011;6:29. [PMC free article] [PubMed]
21. Maynard S, Swistowska AM, Lee JW, Liu Y, Liu ST, et al. Human embryonic stem cells have enhanced repair of multiple forms of DNA damage. Stem Cells. 2008;26:2266–2274. [PMC free article] [PubMed]
22. Adams BR, Hawkins AJ, Povirk LF, Valerie K. ATM-independent, high-fidelity nonhomologous end joining predominates in human embryonic stem cells. Aging (Albany NY) 2010;2:582–596. [PMC free article] [PubMed]
23. Sokolov MV, Neumann RD. Radiation-induced bystander effects in cultured human stem cells. PLoS One. 2010;5:e14195. [PMC free article] [PubMed]
24. Sokolov MV, Panyutin IV, Panyutin IG, Neumann RD. Dynamics of the transcriptome response of cultured human embryonic stem cells to ionizing radiation exposure. Mutat Res. 2011;709–710:40–48. [PMC free article] [PubMed]
25. Sokolov MV, Panyutin IV, Onyshchenko MI, Panyutin IG, Neumann RD. Expression of pluripotency-associated genes in the surviving fraction of cultured human embryonic stem cells is not significantly affected by ionizing radiation. Gene. 2010;455:8–15. [PMC free article] [PubMed]
26. Banuelos CA, Banath JP, MacPhail SH, Zhao J, Eaves CA, et al. Mouse but not human embryonic stem cells are deficient in rejoining of ionizing radiation-induced DNA double-strand breaks. DNA Repair (Amst) 2008;7:1471–1483. [PubMed]
27. Adams BR, Golding SE, Rao RR, Valerie K. Dynamic dependence on ATR and ATM for double-strand break repair in human embryonic stem cells and neural descendants. PLoS One. 2010;5:e10001. [PMC free article] [PubMed]
28. Momcilovic O, Choi S, Varum S, Bakkenist C, Schatten G, et al. Ionizing radiation induces ataxia telangiectasia mutated-dependent checkpoint signaling and G(2) but not G(1) cell cycle arrest in pluripotent human embryonic stem cells. Stem Cells. 2009;27:1822–1835. [PMC free article] [PubMed]
29. Momcilovic O, Knobloch L, Fornsaglio J, Varum S, Easley C, et al. DNA damage responses in human induced pluripotent stem cells and embryonic stem cells. PLoS One. 2010;5:e13410. [PMC free article] [PubMed]
30. Wright GW, Simon RM. A random variance model for detection of differential gene expression in small microarray experiments. Bioinformatics. 2003;19:2448–2455. [PubMed]
31. Sokolov MV, Neumann RD, Panyutin IG. Effects of DNA-targeted ionizing radiation produced by 5-[125I]iodo-2′-deoxyuridine on global gene expression in primary human cells. BMC Genomics. 2007;8:192. [PMC free article] [PubMed]
32. Sokolov MV, Smirnova NA, Camerini-Otero RD, Neumann RD, Panyutin IG. Microarray analysis of differentially expressed genes after exposure of normal human fibroblasts to ionizing radiation from an external source and from DNA-incorporated iodine-125 radionuclide. Gene. 2006;382:47–56. [PubMed]
33. Betel D, Wilson M, Gabow A, Marks DS, Sander C. The resource: targets and expression. Nucleic Acids Res. 2008;36:D149–153. [PMC free article] [PubMed]
34. Dennis G, Jr, Sherman BT, Hosack DA, Yang J, Gao W, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:P3. [PubMed]
35. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. [PubMed]
36. Williams AE. Functional aspects of animal microRNAs. Cell Mol Life Sci. 2008;65:545–562. [PubMed]
37. Van Hoof D, Munoz J, Braam SR, Pinkse MW, Linding R, et al. Phosphorylation dynamics during early differentiation of human embryonic stem cells. Cell Stem Cell. 2009;5:214–226. [PubMed]
38. Yao Y, Suo AL, Li ZF, Liu LY, Tian T, et al. MicroRNA profiling of human gastric cancer. Mol Med Report. 2009;2:963–970. [PubMed]
39. Pizzimenti S, Ferracin M, Sabbioni S, Toaldo C, Pettazzoni P, et al. MicroRNA expression changes during human leukemic HL-60 cell differentiation induced by 4-hydroxynonenal, a product of lipid peroxidation. Free Radic Biol Med. 2009;46:282–288. [PubMed]
40. Adachi K, Suemori H, Yasuda SY, Nakatsuji N, Kawase E. Role of SOX2 in maintaining pluripotency of human embryonic stem cells. Genes Cells. 2010;15:455–470. [PubMed]
41. Ji AR, Ku SY, Cho MS, Kim YY, Kim YJ, et al. Reactive oxygen species enhance differentiation of human embryonic stem cells into mesendodermal lineage. Exp Mol Med. 2010;42:175–186. [PMC free article] [PubMed]
42. Xia L, Zhang D, Du R, Pan Y, Zhao L, et al. miR-15b and miR-16 modulate multidrug resistance by targeting BCL2 in human gastric cancer cells. Int J Cancer. 2008;123:372–379. [PubMed]
43. Qian S, Ding JY, Xie R, An JH, Ao XJ, et al. MicroRNA expression profile of bronchioalveolar stem cells from mouse lung. Biochem Biophys Res Commun. 2008;377:668–673. [PubMed]
44. Ofir M, Hacohen D, Ginsberg D. miR-15 and miR-16 Are direct transcriptional targets of E2F1 that limit E2F-induced proliferation by targeting cyclin E. Mol Cancer Res. 2011;9:440–447. [PubMed]
45. Nishi H, Ono K, Iwanaga Y, Horie T, Nagao K, et al. MicroRNA-15b modulates cellular ATP levels and degenerates mitochondria via Arl2 in neonatal rat cardiac myocytes. J Biol Chem. 2010;285:4920–4930. [PMC free article] [PubMed]
46. Cao H, Yang CS, Rana TM. Evolutionary emergence of microRNAs in human embryonic stem cells. PLoS One. 2008;3:e2820. [PMC free article] [PubMed]
47. Subramanyam D, Lamouille S, Judson RL, Liu JY, Bucay N, et al. Multiple targets of miR-302 and miR-372 promote reprogramming of human fibroblasts to induced pluripotent stem cells. Nat Biotechnol. 2011;29:443–448. [PubMed]
48. Barroso-del Jesus A, Lucena-Aguilar G, Menendez P. The miR-302-367 cluster as a potential stemness regulator in ESCs. Cell Cycle. 2009;8:394–398. [PubMed]
49. Card DA, Hebbar PB, Li L, Trotter KW, Komatsu Y, et al. Oct4/Sox2-regulated miR-302 targets cyclin D1 in human embryonic stem cells. Mol Cell Biol. 2008;28:6426–6438. [PMC free article] [PubMed]
50. Trompeter HI, Abbad H, Iwaniuk KM, Hafner M, Renwick N, et al. MicroRNAs MiR-17, MiR-20a, and MiR-106b act in concert to modulate E2F activity on cell cycle arrest during neuronal lineage differentiation of USSC. PLoS One. 2011;6:e16138. [PMC free article] [PubMed]
51. Rizzo M, Mariani L, Pitto L, Rainaldi G, Simili M. miR-20a and miR-290, multi-faceted players with a role in tumourigenesis and senescence. J Cell Mol Med. 2010;14:2633–2640. [PubMed]
52. Wang Y, Scheiber MN, Neumann C, Calin GA, Zhou D. MicroRNA Regulation of Ionizing Radiation-Induced Premature Senescence. Int J Radiat Oncol Biol Phys 2010 [PMC free article] [PubMed]
53. Lin RJ, Lin YC, Yu AL. miR-149* induces apoptosis by inhibiting Akt1 and E2F1 in human cancer cells. Mol Carcinog. 2010;49:719–727. [PubMed]
54. Li Y, Mine T, Ioannides CG. Short GC-rich RNA similar to miR 1909 and 1915 folds in silico with the 5′-UTR and ORF of Notch and responders: potential for the elimination of cancer stem cells. Oncol Rep. 2010;24:1443–1453. [PubMed]
55. Zhao JJ, Yang J, Lin J, Yao N, Zhu Y, et al. Identification of miRNAs associated with tumorigenesis of retinoblastoma by miRNA microarray analysis. Childs Nerv Syst. 2009;25:13–20. [PubMed]
56. Nurul-Syakima AM, Yoke-Kqueen C, Sabariah AR, Shiran MS, Singh A, et al. Differential microRNA expression and identification of putative miRNA targets and pathways in head and neck cancers. Int J Mol Med 2011 [PubMed]
57. Shigoka M, Tsuchida A, Matsudo T, Nagakawa Y, Saito H, et al. Deregulation of miR-92a expression is implicated in hepatocellular carcinoma development. Pathol Int. 2010;60:351–357. [PubMed]
58. Chen ZL, Zhao XH, Wang JW, Li BZ, Wang Z, et al. microRNA-92a promotes lymph node metastasis of human esophageal squamous cell carcinoma via E-cadherin. J Biol Chem. 2011;286:10725–10734. [PMC free article] [PubMed]
59. Guled M, Lahti L, Lindholm PM, Salmenkivi K, Bagwan I, et al. CDKN2A, NF2, and JUN are dysregulated among other genes by miRNAs in malignant mesothelioma -A miRNA microarray analysis. Genes Chromosomes Cancer. 2009;48:615–623. [PubMed]
60. Nymark P, Guled M, Borze I, Faisal A, Lahti L, et al. Integrative analysis of microRNA, mRNA and aCGH data reveals asbestos- and histology-related changes in lung cancer. Genes Chromosomes Cancer. 2011;50:585–597. [PubMed]
61. Lo SS, Hung PS, Chen JH, Tu HF, Fang WL, et al. Overexpression of miR-370 and downregulation of its novel target TGFbeta-RII contribute to the progression of gastric carcinoma. Oncogene 2011 [PubMed]
62. Schmitz KJ, Helwig J, Bertram S, Sheu SY, Suttorp AC, et al. Differential expression of microRNA-675, microRNA-139-3p and microRNA-335 in benign and malignant adrenocortical tumours. J Clin Pathol. 2011;64:529–535. [PMC free article] [PubMed]
63. Ni CW, Qiu H, Jo H. MicroRNA-663 upregulated by oscillatory shear stress plays a role in inflammatory response of endothelial cells. Am J Physiol Heart Circ Physiol. 2011;300:H1762–1769. [PubMed]
64. Pan J, Hu H, Zhou Z, Sun L, Peng L, et al. Tumor-suppressive mir-663 gene induces mitotic catastrophe growth arrest in human gastric cancer cells. Oncol Rep. 2010;24:105–112. [PubMed]
65. Bandiera S, Ruberg S, Girard M, Cagnard N, Hanein S, et al. Nuclear outsourcing of RNA interference components to human mitochondria. PLoS One. 2011;6:e20746. [PMC free article] [PubMed]
66. Saretzki G, Walter T, Atkinson S, Passos JF, Bareth B, et al. Downregulation of multiple stress defense mechanisms during differentiation of human embryonic stem cells. Stem Cells. 2008;26:455–464. [PubMed]
67. Foshay KM, Gallicano GI. miR-17 family miRNAs are expressed during early mammalian development and regulate stem cell differentiation. Dev Biol. 2009;326:431–443. [PubMed]
68. Fang L, Deng Z, Shatseva T, Yang J, Peng C, et al. MicroRNA miR-93 promotes tumor growth and angiogenesis by targeting integrin-beta8. Oncogene. 2011;30:806–821. [PubMed]
69. Petrocca F, Visone R, Onelli MR, Shah MH, Nicoloso MS, et al. E2F1-regulated microRNAs impair TGFbeta-dependent cell-cycle arrest and apoptosis in gastric cancer. Cancer Cell. 2008;13:272–286. [PubMed]
70. Borgdorff V, Lleonart ME, Bishop CL, Fessart D, Bergin AH, et al. Multiple microRNAs rescue from Ras-induced senescence by inhibiting p21(Waf1/Cip1). Oncogene. 2010;29:2262–2271. [PubMed]
71. Lal A, Navarro F, Maher CA, Maliszewski LE, Yan N, et al. miR-24 Inhibits cell proliferation by targeting E2F2, MYC, and other cell-cycle genes via binding to “seedless” 3′UTR microRNA recognition elements. Mol Cell. 2009;35:610–625. [PMC free article] [PubMed]
72. Lal A, Kim HH, Abdelmohsen K, Kuwano Y, Pullmann R, Jr, et al. p16(INK4a) translation suppressed by miR-24. PLoS One. 2008;3:e1864. [PMC free article] [PubMed]
73. Lal A, Pan Y, Navarro F, Dykxhoorn DM, Moreau L, et al. miR-24-mediated downregulation of H2AX suppresses DNA repair in terminally differentiated blood cells. Nat Struct Mol Biol. 2009;16:492–498. [PMC free article] [PubMed]
74. Pilkinton M, Sandoval R, Song J, Ness SA, Colamonici OR. Mip/LIN-9 regulates the expression of B-Myb and the induction of cyclin A, cyclin B, and CDK1. J Biol Chem. 2007;282:168–175. [PubMed]
75. Baek D, Villen J, Shin C, Camargo FD, Gygi SP, et al. The impact of microRNAs on protein output. Nature. 2008;455:64–71. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science