{'fraction_non_alphanumeric': 0.044094103600792335, 'fraction_numerical': 0.032670338914063984, 'mean_word_length': 4.866602992677492, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
cel-mir-237 and its homologue, hsa-miR-125b, modulate the cellular response to ionizing radiation

Elucidating the mechanisms involved in sensitizing radioresistant tumors to ionizing radiation (IR) treatments while minimizing injury to surrounding normal tissue is an important clinical goal. Due to their sequence-derived specificity and properties as gene regulators in IR-affected pathways, microRNAs (miRNAs) could serve as adjuvant therapeutic agents that alter cellular sensitivity to radiation treatment. To identify radiosensitizing miRNAs, we initially utilized the C. elegans vulval cell model, an in vivo system developed to study IR-dependent radiosensitivity as a measure of clonogenic cell death. We tested several candidate miRNA deletion mutants post γ-irradiation and identified cel-mir-237 as a miRNA which when deleted caused animals to be more resistant to IR, while cel-mir-237 overexpressing strains were IR-sensitive. Additionally, wild-type animals downregulated cel-mir-237 levels post IR in a time-dependent manner. We identified jun-1 (JUN transcription factor homolog) as a novel target of cel-mir-237. Specifically, jun-1 transcript levels increased in wild-type animals post-γ-irradiation, and loss of cel-mir-237 also resulted in higher jun-1 expression. As expected, loss of jun-1 resulted in IR sensitivity, similar to the phenotype of cel-mir-237 overexpressors. Since miR-237 is the homologue of human miR-125, we validated our findings in MCF-7 and MDA-MB-231 breast cancer cell lines, which harbor lower hsa-miR-125b levels than normal HMECs. Forced expression of hsa-miR-125b in these cells resulted in radiosensitivity, as seen by reduced clonogenic survival, enhanced apoptotic activity, and enhanced senescence post IR. Finally, re-expression of c-JUN in MDA-MB-231 cells promoted radio-resistance and abrogated miR-125-mediated radio-sensitization. Our findings suggest that overexpression of cel-mir-237 and its homologue, hsa-miR-125b, functions as sensitizers to γ-irradiation in both a nematode in vivo model and breast cancer cells, and could potentially be utilized as an adjuvant therapeutic to enhance radiation sensitivity.


Introduction
Ionizing radiation (IR) affects cells by generating DNA double-strand breaks (DSB) resulting in a cellular response involving a complex network of pathways that includes DNA replication, DNA-repair, and programmed cell death 1 . These DSBs are repaired through homologous recombination (HR) or non-homologous end-joining (NHEJ) DNA-repair mechanisms 2,3 . This process of normal DNA-repair is disrupted in cancer cells, and is one of the reasons why cytotoxic radiation-based therapeutic approaches are still a primary treatment modality for many cancer patients. Unfortunately, despite a number of advancements in the field of radiobiology, a significant number of patients suffer from locoregional recurrence due to the outgrowth of radioresistant cells 4 . Dose-escalation approaches to eradicating residual tumor cells are met with significant acute and chronic toxicities that arise from the effects of radiotherapy on surrounding normal tissue, limiting the use and subsequent effectiveness of this treatment. Therefore it is essential to develop novel strategies to sensitize cancer cells to radiation-based therapies.
MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression posttranscriptionally, either through mRNA degradation or translational inhibition [5][6][7] . MiRNAs regulate gene expression by base pairing with complementary sequences in the 3′ untranslated region (3′ UTR) of protein-coding transcripts, and function as master regulators of gene expression 7,8 . Dysregulation of a particular miRNA family can promote disease states such as cancer by modulating cellular pathways that control differentiation, apoptosis, and survival [9][10][11][12] . Restoration of these cancer-associated miRNAs either by inhibition (i.e., antagomiRs) or re-introduction (i.e., miRNA mimics) strategies has been shown to reduce tumorigenic properties such as cell growth and invasion, while promoting apoptosis and sensitivity to chemotherapeutic agents [13][14][15][16][17] . Importantly, miRNAs have previously been shown to be critical in the cellular response to radiation 18 , and to act as modifiers of the radiation response both in vivo in a Caenorhabditis elegans model of radiosensitivity, and in cancer cells [19][20][21] .
In this study, we utilized a C. elegans model of radiosensitivity to assay for radiationinduced reproductive cell death 22 . Tissue multipotent precursor cells have an ability to selfrenew and divide to generate invariant cell lineages in a manner analogous to the cancer stem cell system [23][24][25] , and this model has been shown to represent non-apoptotic, clonogenic cell death, analogous to the most common form of cell death in tumors. Therefore, use of this model to further confirm the genetic components of radiosensitivity of these cells is a useful step toward altering radio-therapeutic outcome.
We sought to identify miRNAs that could function as sensitizers to γ-irradiation, which we envision could be used as a neoadjuvant during IR treatment. We tested several candidate miRNA deletion mutant nematodes after γ-irradiation for a radiosensitive phenotype, and identified cel-mir-237 as a miRNA able to alter radiation sensitivity. Upon further investigation the cel-miR-237 homologue, hsa-miR-125b, was also found to promote γirradiation sensitivity in MCF-7 and MDA-MB-231 breast cancer cell lines. Both cel-miR-237 and hsa-miR-125b induce this radiation-sensitivity in part by targeting and downregulating JUN. Importantly, altering c-JUN expression in either C. elegans or breast cancer cells lines resulted in altered sensitivity to radiation treatments. This data suggests that miR-125b could be used as a potential adjuvant to radiation therapies to enhance radiosensitivity.

cel-mir-237(tm2238) deletion results in radioresistance in C. elegans cells after γirradiation
We previously identified miRNAs with altered expression in response to γ-irradiation 20 . To confirm which of these responsive miRNAs could function as radio-sensitizers in vivo, we identified the C. elegans homologues of these miRNAs (Table 1). We then performed experiments where various C. elegans strains harboring loss-of-function mutations in these miRNA genes (herein termed deletion mutants) were irradiated and scored for radiosensitivity based on the presence of previously defined γ-irradiation-dependent morphological vulval defects 22 . In this model, vulval defects are the metric of radiosensitivity, and represent reproductive cell death. While a majority of the deletion mutant strains tested resulted in no change or in sensitization to γ-irradiation (Fig. S1A &  Table 1), only the mir-237(tm2238) deletion mutant displayed a significant radioresistant phenotype (Fig. 1A). Furthermore, cel-mir-237 levels significantly decreased in wild-type N2 animals 6-9hrs post IR, and remained low at later time points (Fig. 1B). This is in contrast to other miRNAs tested, where the respective levels change ∼16-24hrs post γ-irradiation (Fig.  S1B). The only other miRNA whose levels decreased post-IR (Fig. 1B) and whose loss of function promotes an IR-resistant phenotype in C. elegans is let-7 20 . This finding suggests that loss of cel-mir-237 predicts radioresistance.
We then set out to verify whether jun-1 and cel-mir-237 operated within the same genetic pathway. jun-1 is predicted to be a target of mir-237/lin-4 family of miRNAs (Fig. S2A), and the mir-237-complemetary region within the C. elegans jun-1 3′ UTR is conserved in other nematode species (Fig. S2B). Furthermore, the murine miR-237 homologue, mmu-miR-125b, has been previously shown to directly target Jun 33 , and in situ analyses by Feng et al. and Esquela-Kerscher et al. confirmed that both jun-1 (Fig. S2C) 34 and mir-237 35 are expressed in a number of tissues including vulval muscle cells. We therefore examined jun-1 levels in N2 animals in comparison to the mir-237 deletion mutant and overexpressor strains. The levels of jun-1 increased in animals with the mir-237 locus deleted, while jun-1 levels decreased in the cel-mir-237 o/e strains, as compared to N2 animals (Fig. 2B). Furthermore, in N2 animals that had undergone γ-irradiation, the levels of jun-1 significantly increased 9-24hrs post γ-irradiation (Fig. 2C). In the same samples, this expression profile was inversely correlated with cel-mir-237. Notably, jun-1 levels became upregulated only at time-points after the observed changes in cel-mir-237 levels, consistent with the notion of jun-1 being controlled by cel-mir-237.
To test whether jun-1 functioned as an integral part of the γ-irradiation response, jun-1(gk551) and jun-1(gk557) mutant strains were utilized to determine whether these lossof-function lines would phenocopy the cel-mir-237-induced γ-irradiation phenotype, and show radiosensitivity. jun-1 mutants do not have any obvious vulval defects prior to γirradiation. However, when subjected to increasing doses of γ-irradiation, both jun-1 mutant lines were more radiosensitive than N2 animals (Fig. 2D). To assess whether jun-1 and celmir-237 operated within the same genetic pathway, epistatic experiments were performed. Here, jun-1 mutants were each crossed with the radioresistant mir-237(tm2238) mutant strain. The resultant mir-237;jun-1 double mutants still exhibited radiosensitivity (Fig. 2E). These results indicate that jun-1 is both epistatic to and genetically downstream of celmir-237, and that loss of jun-1 due to regulation by cel-mir-237 sensitizes cells to DNA damaging agents such as γ-irradiation.

The cel-mir-237 homologue, hsa-miR-125b, sensitizes breast cancer cells to γ-irradiation
Given the radiosensitizing role of cel-mir-237 in our in vivo model of radiosensitivity, we investigated whether the homologous human miRNA, hsa-miR-125 (Fig. 3A), could promote a similar phenotype in cancer cells. While it is unclear which miR-125 isoform is most closely related to cel-mir-237 we focused on hsa-miR-125b, since meta-data from miRBase indicated that it was expressed at higher levels than hsa-miR-125a. Furthermore, in METABRIC, the mean log2 expression of hsa-miR-125b-5p across all 1262 breast cancer samples was 13.38, as compared to 7.2, 8.59, or not detectable for hsa-miR-125a-3p, hsa-miR-125a-5p, or hsa-miR-125b-3p, respectively. Since proper development of C. elegans vulval precursor cells depend upon EGFR and NOTCH1 signaling, two aberrantly activated pathways in human breast cancer, we profiled both normal and breast cancer cell lines. In a vast majority of breast cancer lines miR-125b was downregulated compared to HMECs, MCF-10As, and a panel of conditionally reprogrammed normal breast cells (CRCs) 36 (Fig.  3B & Fig. S3A), supporting the notion that miR-125b broadly functions as a tumor suppressor in breast cancer cell lines [37][38][39] .
To investigate the role of miR-125b during γ-irradiation, HMECs, MCF-7, and MDA-MB-231 cell lines were profiled after treatment. Normal HMECs downregulated miR-125b levels as early as 4hrs post-IR, which rebounded by 16hrs post-IR( Fig. 3E & 3H). Interestingly, re-introduction of miR-125b alone promoted a pro-differentiation phenotype where cells became large, formed vacuolar structures, and had a heightened granular appearance (Fig. S3B). The persistent endogenous miR-125b levels observed post-IR coupled with the pro-differentiation phenotype could, in part, explain why ectopic expression of miR-125b did not sensitize these cells to γ-irradiation, but rather promoted a mild protective effect( Fig. S3C-E & 3J, right panel). Furthermore, addition of a miR-125b inhibitor did not protect HMECs from IR-induced cell death, as compared a control miRNA inhibitor (Fig. S3H). Overall, these findings suggest that while miR-125b itself functions by regulating cellular proliferation and differentiation pathways in normal mammary epithelial cells, it does not alter the sensitivity of these cells to DNA-damaging agents such as γirradiation.
Contrary to the effects in HMECs, breast cancer cell lines have variable responses to γirradiation. In MDA-MB-231 cells, miR-125b levels rapidly decreased in a dose-dependent manner and remained low after exposure to γ-irradiation( Fig. 3D & 3G), while in MCF-7 cells miR-125b levels fluctuated post γ-irradiation, which is most likely reflective of the low baseline levels of miR-125b( Fig. 3C & 3F). Importantly, in contrast to HMECs, both MCF-7 and MDA-MB-231 cell lines ectopically expressing miR-125b were sensitized to γirradiation( Fig. 3I), as compared to lines expressing a scrambled RNA control (miR-Scr). When MCF-7 cells were transfected with a miR-125b inhibitor γ-irradiation sensitivity was not altered (Fig. S3F), as expected given the low basal levels of miR-125b in these cells, while the miR-125b inhibitor promoted a mild protective effect in MDA-MB-231 cells (Fig.  S3G).

hsa-miR-125b promotes γ-irradiation induced senescence
To investigate the anti-growth phenotypes associated with miR-125b re-introduction post γirradiation, we functionally assessed breast cancer cell lines ectopically expressing miR-125b. We first tested cell proliferation by PI staining and found that while γ-irradiation promoted G2/M arrest 40,41 , miR-125b did not significantly alter cell cycle distribution either alone or in combination with γ-irradiation (Fig. S4A). Instead, miR-125b transfected MDA-MB-231 cells underwent striking morphological changes post-γ-irradiation indicative of senescence( Fig. 4A-B), and SA-β-gal assays performed on miR-125b transfected cells confirmed this senescence-promoting effect( Fig. 4C-D). We also observed significantly enhanced SA-β-gal levels in MCF-7 cells co-treated with miR-125b and 4-6-Gy γirradiation, as compared to miR-Scr and γ-irradiation treatments (Fig. 4E). These miR-125bdependent effects were not observed in HMECs after γ-irradiation (Fig. S4C). Additionally, Annexin V staining assays indicated miR-125b transfected MDA-MB-231 cells were less likely to undergo γ-irradiation-induced apoptosis and DNA damage( Fig. S4B & S4D), further supporting the notion that miR-125b promotes a cytostatic phenotype when coupled with γ-irradiation. Together, these results suggest that while miR-125b alone does not induce cellular senescence, it does play an important role in priming breast cancer cells for subsequent radiation-induced senescence.

Aberrant expression of c-JUN confers γ-irradiation sensitivity
Given mir-237/miR-125b modulation of jun-1 was important in γ-irradiation sensitivity in C. elegans, we next asked whether the mammalian homologue of jun-1, JUN, was important during this process. We first screened for JUN levels in breast cancer lines, and found MCF-7 cells, which harbored the lowest miR-125b levels, had the highest levels of JUN( To further assess the role of c-JUN, we generated c-JUN-ORF-IRES-GFP-MDA-MB-231 cell lines (Fig. S5A). We chose this line because it expressed low endogenous levels of JUN, and we observed a robust increase in JUN levels after γ-irradiation treatment, suggesting that c-JUN is important for cell survival during this process. This line expressed >14-fold more JUN transcript than an Empty-GFP-control line, and was not responsive to miR-125b mimic treatments (Fig. 5C). Endogenous JUN did not change in c-JUN-ORF-IRES-GFP versus parental lines, but decreased in abundance in the presence of miR-125b, as determined by a JUN 5′ UTR primer set, which is expected given miR-125b targets JUN 42 .
Also as expected, c-JUN-ORF-IRES-GFP lines were more resistant to γ-irradiation as compared to GFP-control lines (Fig. 5D), and did not demonstrate sensitivity to miR-125binduced γ-irradiation as compared to an Empty-GFP control line (Fig. 5E). We subsequently performed loss of function experiments in JUN high expressing MCF-7 cells, and found that two independent siRNAs to JUN could weakly sensitize these cells to γ-irradiation( Fig. 5F & S5B). Since there are multiple c-JUN family members that can generate an AP-1 transcription factor complex with c-FOS 43 , MCF-7 cells were instead treated with Tanshinone IIA (Fig. S5C), an inhibitor of AP-1 DNA binding activity, which resulted in enhanced γ-irradiation sensitivity (Fig. S5D-E).

hsa-miR-125b regulates MAPK signaling and activity of c-JUN
While investigating the role of the miR-125b-c-JUN regulatory axis during γ-irradiationinduced stress, we found that miR-125b inhibited the activity of c-JUN. Specifically, after 4Gy treatment, transfection with miR-125b mimic reduced c-JUN protein expression ∼1.5to 1.8-fold in MCF-7 and MDA-MB-231 cells, respectively( Fig. 6A-B). This is expected given previous reports of miR-125b targeting JUN 42 . However, these reductions were not observed in the absence of γ-irradiation( Fig. 6C-D). This indicated that miR-125b might also regulate a process that alters the activity and subsequently the stability of c-JUN. Indeed pSer73-c-JUN levels were downregulated 1.5-and 3.0-fold in miR-125b transfected MCF-7 and MDA-MB-231 cells, respectively (Fig. 6C-D), indicating the activity of upstream Jun Nterminal kinases (JNKs) was attenuated by miR-125b. To identify which signaling pathways were most affected by miR-125b, we performed qPCR array experiments on cells treated with γ-irradiation or miR-125b. We then overlapped genes that were 2-fold downregulated by miR-125b and 2-fold upregulated after γ-irradiation treatment, with all putative miR-125b targets as predicted by TargetScan (Fig. 6E). Using this analysis we identified several upstream kinases such as MAP3K1 and MAP2K7, that were strongly elevated after γ-irradiation in both cell lines, yet were downregulated by miR-125b (Table S3).
Interestingly, MAP2K7 (JNK-activating kinase), MAP3K1, and MEKK1, have already been shown to be targets of miR-125b 44,45 . This analysis also identified several genes that were downregulated by miR-125b independent of γ-irradiation such as MOS and HSPB1 in MCF-7 cells, CDKN1A and SFN in MDA-MB-231 cells, and EGR1 in both cell lines.
Given many of these genes are considered pro-tumorigenic, these findings support the notion that miR-125b has tumor-suppressive capabilities in breast cancer independent of promoting γ-irradiation sensitivity. While further studies are needed to determine which genes from our analysis are also targeted and regulated during miR-125b-induced γ-irradiation sensitivity, overall these data indicate that miR-125b can promote sensitivity to γ-irradiation by modulating the MAPK-c-JUN axis resulting in reduced expression of genes important for cell survival during DNA-damage induced cell stress (Fig. 6F).

Discussion
Studying γ-irradiation response in vivo utilizing gene knockout models offers opportunities to determine how one genetic factor affects an entire organism after γ-irradiation-induced stress. Here we report the parallel roles of the two homologous miRNAs, cel-mir-237 and hsa-miR-125b, in radiosensitization of a C. elegans model of clonogenic cell death and human breast cancer cells to γ-irradiation. In N2 nematodes, cel-mir-237 was downregulated upon γ-irradiation, and functionally the loss of cel-mir-237 protected cells from γirradiation, while overexpression of cel-mir-237 promoted a radiosensitive phenotype. These results suggest that cells may actively downregulate cel-mir-237 as a mechanism to promote survival during and after γ-irradiation. A similar mechanism may be in place in breast cancer cells, as we identified miR-125b to be downregulated after γ-irradiation treatment, and the forced expression of miR-125b to impart γ-irradiation sensitivity in MDA-MB-231 and MCF-7 cells, but not in normal HMEC cells. We also demonstrated that the JUN family proteins: jun-1 in C. elegans and c-JUN in humans, were regulated by cel-mir-237 and hsa-miR-125b, respectively, and that altered expression of this family of proteins could modulate γ-irradiation sensitivity.
cel-mir-237 is part of a miRNA family that includes lin-4 and mir-237; is upregulated during the L2 to L3 transition in C. elegans 35,46,47 , and may share certain functions with the heterochronic miRNA gene lin-4 48 . However, no phenotype has yet been reported for loss of this particular miRNA. Initially, we did not observe any gross morphologic or phenotypic defects in cel-mir-237(tm2238) deletion mutants maintained in an unstressed state, including those phenotypes tested by Miska et al 49 . However, once cel-mir-237(tm2238) deletion mutants were stressed with γ-irradiation, we observed striking vulval morphological defects due radiosensitivity and reproductive cell death. This is interesting, given most miRNArelated phenotypes found in vivo require compound deletions of related miRNAs within a particular miRNA family 50 . However, once coupled with a stress condition, such as γirradiation, dysregulation of mir-237 can confer in vivo-related morphological defects, in addition to those already identified for miR-34 and let-7 19,20 .
Our results also indicate that cel-mir-237 may function by regulating jun-1. Using sequence and conservation analysis, jun-1 was predicted to be a target of cel-mir-237 by two different algorithms 27,28 . In addition, it was previously reported that jun-1 levels were upregulated when the miRNA biogenesis gene dcr-1 was mutated, but did not change in rde-1 or rde-4 RNAi pathway-specific gene mutants 51 . Here, we have demonstrated an inverse relationship between jun-1 and cel-mir-237 levels not only in overexpressor and deletion mutants, but also in N2 animals following γ-irradiation treatment. Notably, following γ-irradiation, jun-1 levels increased hours after cel-mir-237 levels decreased. This result hints at a possible temporal correlation between cel-mir-237 and jun-1 post-γ-irradiation. Functionally, celmir-237 overexpressor strains and jun-1 loss-of-function mutant strains similarly exhibited higher sensitivities to γ-irradiation than N2 animals. Finally, in an epistatic study, jun-1 lossof-function and mir-237(tm2238) double mutants were more sensitive to γ-irradiation than N2 animals, which was similar to the phenotype of jun-1 loss-of-function mutants, but opposite of that of mir-237(tm2238) mutants, indicating jun1 and cel-mir-237 are within the same genetic pathway. Taken together, these results are consistent with the hypothesis that cel-mir-237 downregulates jun-1 levels and that jun1 expression is critically important in modulating cellular survival cues post-γ-irradiation.
Development of C. elegans vulval precursor cells depends upon signaling pathways that are also aberrantly activated in human breast cancer; therefore, we hypothesized that the mammalian cel-mir-237 homologue, hsa-miR-125b, could promote a similar radiosensitizing phenotype. Indeed, MCF-7 and MDA-MB-231 breast cancer cells transfected with miR-125b mimic were more sensitive to γ-irradiation as compared to miR-Scr control treated cells. Our results are consistent with a study in oral squamous cell carcinoma cell lines in which forced expression of miR-125b, acting through the target gene, ICAM2, sensitized cells to X-ray-irradiation 52 . On the other hand, Tan et al. observed that overexpression of miR-125b promoted cell survival after UV-based irradiation through the downregulation of p38α 53 . These conflicting reports might be due to the fact that Tan et al.
investigated the effects of UV-based irradiation, which promotes a different form of DNA damage, activating nucleotide excision repair rather than HR or NHEJ. Furthermore, Tan et al. studied the acute effects of irradiation in miR-125b transfected cells, which was mainly a phenotype driven by apoptosis. However, our study utilized clonogenic assays to assess the long-term consequences of γ-irradiation in miR-125b expressing cells, which would capture both IR-induced apoptosis as well as delayed cell death.
The role of miR-125b in modulating irradiation-sensitivity in breast cancer is most likely a cell-context specific phenotype, given miR-125b itself has been described as a tumor suppressor in breast cancer 37,54-56 , melanoma 57 , and bladder cancer 58 , but also as an oncogene in prostate cancer 59 and leukemia 60,61 . miR-125b has also been shown to promote stem cell-like properties in cholangiocarcinoma 62 and expand progenitor cell populations in acute megakaryoblastic leukemia 63 . Furthermore, miR-125b is significantly enriched in hair follicle stem cells and promotes initiation and malignant progression of skin tumors 45,64 . The diversity of miR-125b-induced phenotypes can largely be explained by the target mRNA genes that play an essential role in each of the aforementioned processes. As an example, miR-125b has been demonstrated to target genes in pathways involving cell proliferation, cell differentiation, and apoptosis including, TP53, TDG, ERBB2/3, EPO, proapoptotic genes such as BAK1 and PUMA, and cell cycle regulators such as CDC25C 54,59,[65][66][67] .
In this study, we observed that the interaction between miR-125b and JUN controlled γirradiation sensitivity in MCF-7 and MDA-MB-231 breast cancer cells. The consequence of disrupting this regulatory network, either by forced miR-125b expression or inhibition of JUN activity during IR-induced stress, was increased senescence. We did not observe this phenomenon in normal HMEC cells. While it has been previously reported that miR-125b targets and reduces JUN expression during granulocyte differentiation 33,42 and promotes cellular senescence in melanoma 57 , this is the first report identifying these two genes interacting within the same regulatory network to confer γ-irradiation sensitivity in breast cancer. We also identified miR-125b regulated genes involved in stress-response MAPK signaling, which are induced in cancer cells after γ-irradiation, presumably to promote cell survival. Certain breast cancer disease subtypes such as triple-negative breast cancer (TNBC), which are aggressive and have shorter intervals for locoregional recurrence, distant metastasis and disease-free survival [68][69][70] , also harbor a heightened MAPK-gene signature. Future work will be required to determine if use of miR-125b as an adjuvant during IR treatment paradigms could enhance the radiosensitivity of tumor cells in these patient populations a well.
Primers for cloning and characterization are shown in Table S4. Animals were assessed for stable integration, copy number, and levels of integrated gene expression. jun-1(gk551);mir-237(tm2238), and jun-1(gk557);mir-237(tm2238) strains were generated by genetic crossing of each of the two indicated strains. Male nematodes were only used for genetic crossing; all animals examined in the radiation treatment experiments were hermaphrodites.

Determination of radiation treatment dose response
To minimize cell-cycle dependent variation in irradiation sensitivity experiments, synchronized C. elegans strains were examined at a fixed radiation dosage using a Cesium 137 source to determine peak radioresistance post-hatching as discussed in 22 . The time-point showing greatest radioresistance was the point at which the greatest proportion of animals had wild-type vulvae 22 . All experiments were performed twice and radioresistance peak times are shown in Table S5.
Animals were synchronized and placed on OP50-NGM plates for the duration of time as determined by the radioresistance peak assay. Animals were placed in a 15-ml conical tube with OP50-seeded NGM and irradiated at 0-400Gy. Animals were scored for vulval malformations as adults 22 . Normal vulval structures were scored as a wild-type (WT) phenotype; conversely, a protruding vulval (Pvl) or the absence of a vulva structure (Vul) was categorized as a mutant (Mut) phenotype. Dose curves are presented as ratios of animals with WT:Mut vulva. For each dose-response, one data-point represents 3 biological replicates each containing ∼150-300 animals (averages of 245.35 animals per replicate for Fig. 1A; 215.80 for Fig. 1C; 174.33 for Fig. 2D; and 190.27 for Fig. 2E). Statistical analysis was done by 2-tailed paired t-test. All strains were normalized to 0Gy data-points.

Cell lines and transfections
MCF-7 and MDA-MB-231 cells were obtained from the American Type Culture Collection (ATCC), and cultured at 37 o C with 5% CO 2 in RPMI supplemented with 10% FBS and 1% penicillin/streptomycin. HMECs were maintained as per manufacture's protocol (Life Technologies). All experiments were performed within 1-2 years upon receipt of cells, and characterized by STR analysis via ATCC. Cell lines were tested for mycoplasma via the Lonza MycoAlert TM Detection Kit. For RNA-based functional experiments, mirVana mimics and validated Silencer Select siRNAs (Life Technologies) were transfected into 5×10 4 cells using DharmaFECT 1 (GE Healthcare Dharmacon), as per manufacture's protocol.

Plasmids and generation of stable cell lines
To generate the stable c-JUN-ORF-IRES-GFP MDA-MB-231 line, HEK-293T cells were co-transfected with pCMV-VSVG and pCMV-dR8.2 (gag/pol) helper constructs (Addgene) along with a pMIEG3-c-Jun expression plasmid (Addgene) using Fugene 6 (Promega). A pMIEG3-empty control construct was also generated by excising the c-Jun sequence via NotI / XhoI endonuclease digestion. Retrovirus was produced using The RNAi Consortium protocols (www.broadinstitute.org/rnai/public). Stable transductants were double sorted for GFP on a BD FACS Aria flow cytometer.

Clonogenic assays
Unless otherwise indicated, 1×10 5 cells were seeded into 12-well plates and transfected with 30nM miR-125b mimic or negative control mimic, 15nM miR-125b inhibitor or negative control inhibitor, and 15nM si-Jun or scrambled siRNA (Life Technologies) in antibioticfree RPMI media. 24hrs post-transfection cells were re-seeded into 6-well plates for expansion. In some experiments cells were pre-treated with Tanshinone IIA. 72hrs later cells were irradiated at 0-6Gy. Cells were immediately seeded into 6-well plates at ∼100-900 cells per well, depending upon γ-irradiation condition from the lowest number of cells (100-200) at 0Gy, (275) at 2Gy, (400) at 4Gy, to the highest (900) at 6 Gy in order to produce colonies with clearly defined borders and minimal overlapping. After 10 days of undisturbed growth, colonies were fixed, stained with crystal violet, and counted on a dissection scope (Zeiss). Plating efficiency and surviving fraction were calculated from technical quadruplicates of each γ-irradiation dose and miRNA condition. Experiments were performed a minimum of three times. Statistical significance was evaluated by two-way ANOVA (two-tailed) of the data for each transfection or drug treatment across all doses.

RNA-based functional experiments
For cytotoxicity assessment and development of IC 50 curves, Sulforhodamine Blue (SRB) assays were performed as previously described 36 . Briefly, transfected or drug treated cells were seeded into 24-well plates, HMECs were additionally treated with γ-irradiation, and at indicated times cells were fixed, stained, and dye intensity quantified. Additionally, transfected cells underwent γ-irradiation treatment and were analyzed for presence of apoptosis by the Annexin V Apoptosis Detection Kit (BD Pharmigen), presence of senescence by the SA-β-Galactosidase Staining Kit (Cell Signaling Technology) 36,76 , and for cell cycle analysis as previously described 36 . Unless otherwise stated, data was made relative to the initial time point for each treatment condition and presented as the average of at least two independent experiments ±/-SD. p-values were calculated by two-way ANOVA.
For DNA damage assays, transfected cells were treated with either 0 or 4Gy γ-irradiation, embedded in low-melt agarose, and spread onto slides for alkaline lysis and electrophoresis as described previously 77 . Individual cells were scored from 0 to 4, where 0=no DNA damage and 4=extensive tailing indicative of severe DNA damage. Cells were scored from 10 fields per slide with 3 replicates for each condition. p-values were calculated by Fisher's exact test.

RNA isolation and qRT-PCR
C. elegans were synchronized and sample collections occurred at the indicated time postirradiation. 3600 animals were lysed for RNA extraction using the mirVana miRNA Isolation Kit. Cell lines were lysed in TRIzol. Total RNA was treated with Turbo DNase. For mRNA analysis, cDNA synthesis was performed using SuperScript III with oligo dT primers. 2.5ng of input cDNA was used for SYBR Green-based real-time PCR. Gene expression was normalized to RPL19 for mammalian cells, and pmp-3 78 for C. elegans.
For C. elegans miRNA analysis, qRT-PCR was performed using the TaqMan® miRNA system as per manufacturer's instructions, and gene expression was normalized to U18. For mammalian miRNA analysis, qRT-PCR was performed using the miScript II system (Qiagen). MiRNA expression was normalized to RNU6B. qPCR reactions were analyzed on a Roche LightCycler 480 instrument. Primers used in this study are listed in Table S6. All experiments include 3 biological replicates unless otherwise stated in the figure legends.
For qPCR arrays, MCF-7 and MDA-MB-231 breast cancer cells were transfected with 15nM miR-125b or miR-Scr control in duplicate and harvested for RNA 72 hours posttransfection. Lysates were converted to cDNA using the RT 2 First Strand Kit, and run on RT 2 Profiler Human MAPK Signaling PCR Arrays (Qiagen) measuring expression of 84 genes in this pathway. Gene expression was normalized using the geometric mean of GAPDH, HPRT1, and RPLP0 expression, and was then further processed for miRNA targetenrichment analysis. All genes that were 2-fold downregulated by miR-125b and 2-fold upregulated after γ-irradiation treatment were considered candidates for a miR-125b-γirradiation gene signature. Genes that were strongly elevated after γ-irradiation in both cells lines, yet were downregulated by miR-125b are shown in Table S3.

Bioinformatic analysis and miRNA target prediction
Predicted miRNA gene targets were acquired from three published databases: miRanda 26,27 ,TargetScanWorm 28,29 , and mirWIP 30 . Bioinformatic analysis of putative miR-125b targets were performed using TargetScan 7 . All predicted targets were further selected based on, 1) being predicted by more than one database (Table S1) and 2) having known or predicted functions that are involved in DNA repair and DNA damage response 31 , cell survival and cell proliferation, AP-1 complex formation, or were genes from a recent RNAi screen identifying radio-sensitizers in C. elegans 32 . Overlapping miRNA-targets with irradiationresponse genes, or to genes changed after miR-125b transfection, or γ-irradiation treatment was performed using GeneVenn.

Statistical analysis
Statistical analysis was achieved by 2-tailed paired t-test for vulval radiation sensitivity assay, Fisher's exact test for DNA damage comet assay, 2-tailed Student's t-test for gene expression analysis by qRT-PCR, and 2-way ANOVA for clonogenic assays, Annexin V staining assay, and SA-β-gal assay. All graphs show means ± S.D. Additionally for qRT-PCR data, means and S.D. were calculated at the ddCt level before being converted to fold changes as presented in the graphs 79 .

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.   (I) Analysis of clonogenic assays on MCF-7 and MDA-MB-231 cells 10 days following γirradiation. Cell lines transfected with miR-125b mimic 72hrs prior to γ-irradiation harbored lower colony forming potential than those transfected with miR-Scr control. For MCF-7 and MDA-MB-231 cells, p-values were <0.001 and <0.01, respectively, as determined by 2-way ANOVA. (J) For HMEC cells, surviving fraction was determined by Trypan Blue counts 6 days post-γ-irradiation. No significant difference was observed between miR-125b vs. miR-Scr transfected cells.   Quantification of (C); ratios of p-Ser73, total c-Jun, and β-actin expression in miR-125b vs. miR-Scr transfected cells. (E) Venn diagrams depicting results of qPCR array experiments on MCF-7 and MDA-MB-231 cells 72hrs after miR-125b transfection, or after 4Gy treatment. Genes that changed 2-fold or more are reported and are overlapped with all putative miR-125b targets as determined by TargetScan. (F) A model for cel-mir-237/hsa-miR-125b regulation of cell survival after γirradiation. In C. elegans, vulval cells actively downregulate cel-miR-237 (green arrow) to promote jun-1 expression, which mediates vulval cell survival. In breast cancer cells miR-125b operates in a similar manner, but may also act through the MAPK pathway to promote senescence.