The RNA Degradation Pathway Regulates the Function of GAS5 a Non-Coding RNA in Mammalian Cells

Abstract

Studies of various mRNAs have revealed that changes in the abundance of transcripts, through mRNA degradation, act as a critical step in the control of various biological pathways. Similarly, the regulation of non-coding RNA (ncRNA) levels is also considered to be important for their biological functions; however, far less is known about the mechanisms and biological importance of ncRNA turnover for the regulation of ncRNA functions. The growth arrest-specific 5 (GAS5) ncRNA accumulates during growth arrest induced by serum starvation and its transcript is degraded by the well characterized nonsense-mediated RNA decay (NMD) pathway. Historically, NMD was discovered as a RNA quality control system to eliminate aberrant transcripts; however, accumulating evidence shows that NMD also regulates the abundance of physiological transcripts. Interestingly, the GAS5 transcript has the ability to bind the glucocorticoid receptor (GR), resulting in the inhibition of its ligand-dependent association with DNA. The GR binds the promoters of various glucocorticoid-responsive genes, including apoptosis-related genes. In this study, we examined whether the RNA degradation pathway can regulate this function of GAS5. We measured the steady-state abundance and the decay rate of GAS5 in UPF1-depleted human cells using the 5′-bromo-uridine immunoprecipitation chase (BRIC) method, an inhibitor-free method for directly measuring RNA stability. We found that levels of the GAS5 transcript were elevated owing to prolonged decay rates in response to UPF1 depletion, and consequently the apoptosis-related genes, cIAP2 and SGK1, were down-regulated. In addition, serum starvation also increased the transcript levels of GAS5 because of prolonged decay rates, and conversely decreased levels of cIAP2 and SGK1 mRNA. Taken together, we found that the RNA degradation pathway can regulate the function of the GAS5 ncRNA in mammalian cells.

About the Authors:

Hidenori Tani

* E-mail: [email protected] (HT); [email protected] (NA)

Affiliation: Research Institute for Environmental Management Technology, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan

Masaki Torimura

Affiliation: Research Institute for Environmental Management Technology, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan

Nobuyoshi Akimitsu

* E-mail: [email protected] (HT); [email protected] (NA)

Affiliation: Radioisotope Center, The University of Tokyo, Bunkyo, Tokyo, Japan

Introduction

The regulation of mRNA degradation and stability in the cytoplasm has been well characterized. This regulation is a critical step in changing the abundance of certain transcripts controlling various biological pathways [1]–[3]. It has been estimated that 5%–10% of human genes are controlled at the level of RNA stability [4]. Indeed, there are different ways that mRNA stability is regulated. Each individual mRNA has its intrinsic stability under a given condition [5]–[7], and the stability of an individual mRNA may change in response to a variety of extracellular stimuli [8]–[10]. Similarly, the regulation of non-coding RNA (ncRNA) levels is considered important for their biological functions. ncRNAs can be roughly classified into two groups: small transcripts, such as microRNAs (miRNAs) [11], [12] and long transcripts, such as long intergenic non-coding RNAs (lincRNAs) [13], [14]. miRNAs are endogenous, small ncRNAs that play important regulatory roles in gene expression. Recent discoveries indicate that some miRNAs have a long half-life, but miRNA dynamics are important in specific cellular contexts, as well as the requirement of sequence elements for miRNA decay [15], [16]. In agreement with the varied life spans of mRNAs, long ncRNAs also have half-lives that vary. It was found that long ncRNAs with short half-lives included known regulatory ncRNAs, such as HOTAIR, CDKN2B-AS1/ANRIL, and NEAT1/MENε/β [7], [17]. Thus, it is expected there will be regulatory mechanisms for RNA degradation of long ncRNAs, as is the case for mRNA. However, little is known about the regulation of long ncRNAs.

Growth arrest-specific 5 (GAS5) is a long (∼650 bases in humans) ncRNA that was originally isolated from a screen for potential tumor suppressor genes expressed at high levels during growth arrest [18]. The human GAS5 gene is a multiple small nucleolar RNA (snoRNA) host gene that encodes 10 box C/D snoRNAs within 11 introns, and has been classified as a member of the 5′-terminal oligopyrimidine tract (5′ TOP) gene family, characterized by an upstream oligopyrimidine tract sequence [19]. The GAS5 transcript is increased during growth arrest induced by either serum starvation or treatment with translation inhibitors [19]–[22]. GAS5 has little protein-coding potential but its RNA is spliced, polyadenylated, and associates with ribosomes [19]. The GAS5 transcript is ubiquitously expressed, but it is unstable in human and murine proliferating cells [7], [17], [23]. The instability of GAS5 transcripts in proliferating cells can be explained by the phenomenon of nonsense-mediated RNA decay (NMD), whereby ribosome-associated GAS5 transcripts become targeted for degradation because of their premature stop codons [19]. Historically, NMD was discovered as a RNA quality control system to eliminate aberrant transcripts, growing evidences have shown that NMD also regulates the abundance of physiological transcripts [24]–[26]. Elimination of UPF1, an essential NMD factor, remarkably stabilizes the GAS5 transcript [25], [27]. Moreover, the overexpression of GAS5 causes both an increase in apoptosis induced by UV or cisplatin and a reduction in the rate of progression through the cell cycle. Consistently, the downregulation of endogenous GAS5 inhibits apoptosis and maintains a more rapid cell cycle, indicating that GAS5 expression is both necessary and sufficient for normal growth arrest in human cells, and it suggests that GAS5 can act as a tumor suppressor [28], [29]. Recently, Kino et al. reported that GAS5 functions as a starvation- or growth arrest-linked riborepressor for the glucocorticoid receptor (GR) by binding to the DNA-binding domain of the GR, acting as a decoy glucocorticoid response element (GRE), thus competing with DNA GREs for binding to the GR [30]. This results in the inhibition of the ligand-dependent association of the GR with GREs of various glucocorticoid-responsive genes including apoptosis-related genes.

In this study, we examined whether the NMD RNA degradation pathway regulates the function of GAS5 as a model long ncRNA regulation. RNA decay has been assessed by blocking global transcription with traditional transcriptional inhibitors such as actinomycin D (ActD), 5,6-Dichloro-1-β-D-ribofuranosyl-benzimidazole (DRB), and α-amanitin (α-Am) [1], [31]. However, inhibitor-mediated global transcriptional arrest has a profound impact on cellular physiology, and affects RNA decay rates [7], [23], [32]. Thus, we used a novel inhibitor-free method of pulse-labeling nascent transcripts with 5′-bromouridine [7], [26], [33] or with 5′-ethynyluridine [33]–[36] to enable the measurement of RNA decay under non-disruptive conditions. We found that GAS5 levels were elevated owing to prolonged decay rates in response to UPF1 depletion, and that multiple apoptosis-related genes were downregulated. Similarly, serum starvation increased the expression level of GAS5 because of prolonged decay rates, and decreased the mRNA levels of apoptosis-related genes.

Materials and Methods

Cell culture

HEK293T and HEK293 cells were grown in Dulbecco's modified Eagle's medium (DMEM), supplemented with 10% fetal bovine serum (FBS) and antibiotics at 37°C in a humidified incubator with 5% CO2.

siRNA treatments

The sequences of the siRNAs were: control siRNA (sense: 5′-GUACCUGACUAGUCGCAGAAG-3′, antisense: 5′-GUACCUGACUAGUCGCAGAAG-3′), UPF1 siRNA (sense: 5′-GAUGCAGUUCCGCUCCAUUdTdT-3′, antisense: 5′-AAUGGAGCGGAACUGCAUCdTdT-3′), UPF1 siRNA (2) (sense: 5′-AAUUUCUGUAACUUGUUUCCU-3′, antisense: 5′-GAAACAAGUUACAGAAAUUAC-3′), and GAS5 siRNA (sense: 5′-CUUGCCUGGACCAGCUUAAUU-3′, antisense: 5′-UUAAGCUGGUCCAGGCAAGUU-3′). These siRNAs were transfected into cells using Lipofectamine RNAiMAX (Invitrogen), according to the manufacturer's instructions. siRNA duplexes were used at a final concentration of 10 nM. Cells were harvested 48 h after transfection and total RNAs were isolated using RNAiso Plus (Takara), according to the manufacturer's instructions. RT-qPCR analysis was used to determine whether RNA interference achieved significant depletion of each target sequence.

Over-expression of GAS5

A GAS5 (1−631) expression vector was constructed by subcloning the full-length GAS5 sequence lacking a poly A tail (based on the GAS5 sequence, NR_002578, in NCBI). The GAS5 cDNA was amplified from total RNA purified from HeLa cells, and then cloned into pcDNA3.1 (+) (Invitrogen). Additional information on vector construction can be provided upon request. The expression vector, at 1 µg/mL, was transfected into cells using Lipofectamine 2000 (Invitrogen), according to the manufacturer's instructions. Cells were harvested 48 h after transfection and total RNAs were isolated using RNAiso Plus (Takara), according to the manufacturer's instructions. RT-qPCR analysis was used to determine the levels of GAS5 over-expression.

Reverse Transcription - quantitative real-time Polymerase Chain Reaction (RT-qPCR)

The isolated RNA was reverse transcribed into cDNA using the PrimeScript RT Master Mix (Perfect Real Time) (Takara). The cDNA was amplified using the following primer sets: GAPDH (forward: 5′-GCACCGTCAAGGCTGAGAAC-3′, reverse: 5′-TGGTGAAGACGCCAGTGGA-3′), ACTB (forward: 5′-CCAACCGCGAGAAGATGA-3′, reverse: 5′-CCAGAGGCGTACAGGGATAG-3′), GAS5 (forward: 5′-CTTGCCTGGACCAGCTTAAT-3′, reverse: 5′-CAAGCCGACTCTCCATACCT-3′), UPF1 (forward: 5′-AGATCACGGCACAGCAGAT-3′, reverse: 5′-TGGCAGAAGGGTTTTCCTT-3′), MALAT1 (forward: 5′-GCTGTGGAGTTCTTAAATATCAACC-3′, reverse: 5′-TTCTCAATCCTGAAATCCCCTA-3′), cIAP2 (forward: 5′-TCTAGTGTTCTAGTTAATCC-3′, reverse: 5′-ACCACTTGGCATGTTGAACC-3′), SGK1 (forward: 5′-CTATGCTGCTGAAATAGC-3′, reverse: 5′-GTCCGAAGTCAGTAAGG-3′), ADRP (forward: 5′-CTGTGGCCAGCACGATCAC-3′, reverse: 5′-CTCACGAGCTGCATCATC-3′), GADD45A (forward: 5′-TTTGCAATATGACTTTGGAGGA-3′, reverse: 5′-CATCCCCCACCTTATCCAT-3′), SMG1 (forward: 5′-TCGAAGTCAAGAACACGTTGA-3′, reverse: 5′-GGGTGATGCAAAACTCACTAAA-3′), EIF4A3 (forward: 5′-GGATGAAGCTGATGAAATGTTG-3′, reverse: 5′-TGGTCATCTCCAGAATCTCGT-3′), and CASC3 (forward: 5′-GGGGTTCCAGTTAATACAAGTTTC-3′, reverse: 5′-GCCAGCTGTATTTCTCTTCTGAG-3′). GAPDH and ACTB were used for normalization. SYBR Premix Ex Taq II (Perfect Real Time) (TaKaRa) or THUNDERBIRD SYBR qPCR mix (TOYOBO) was used according to the manufacturer's instructions. RT-qPCR analysis was performed using a Thermal Cycler Dice Real Time System (Takara) or MyiQ2 (BIO-RAD).

BRIC

BRIC was performed as previously described [7], [26], [33]. In brief, cells were incubated at 37°C in the presence of 150 µM 5′-bromo-uridine (BrU) (Wako) for 24 h in a humidified incubator with 5% CO2. At indicated time points after replacing BrU-containing medium with BrU-free medium, cells were harvested. Total RNA was isolated using RNAiso Plus (Takara). Two micrograms of BrU-labeled total RNA were denatured by heating at 80°C for 1 min and then added to 2 µg of anti-BrdU mAb-conjugated beads (clone 2B1, MBL). The mixture was incubated at room temperature for 1 h with rotation. Beads were washed three times with 0.1% BSA in PBS. ISOGEN LS (Nippon Gene) was added, followed by RNA isolation in accordance with the manufacturer's instructions. The isolated RNA was used for RT-qPCR.

5-Ethnyluridine (EU) pulse-labeling

Analysis of RNA half-lives was performed by EU pulse-labeling of RNA using the Click-iT Nascent RNA Capture Kit (Invitrogen) [33]–[36]. After serum withdrawal for 24 h, 200 µM of EU was added to medium and cells incubated for an additional 24 h. At indicated time points after replacing EU-containing medium with EU-free medium, cells were harvested. Total RNA was isolated using RNAiso Plus (Takara). Then, EU-labeled RNAs were biotinylated and captured using the Click-iT Nascent RNA Capture Kit (Invitrogen), according to the manufacturer's instructions. Isolated RNA s were used for RT-qPCR.

Results

Half-life of GAS5 RNA

We first determined the decay rate of GAS5 in HEK293T cells using the BRIC system [7], [26], [33]. In the BRIC system (Figure 1A), BrU-labeled RNAs were immunopurified from total RNA isolated from HEK293T cells at four time points after BrU pulse-labeling. Using this method we found the half-life of GAS5 in HEK293T cells to be 6.6 h (Figure 1B). A previous study using the BRIC system reported that the half-life of GAS5 in HeLa cells was 2.6 h [7]. The half-life of GAS5 in HEK293T cells is relatively long compared with that in HeLa cells.

[Figure omitted. See PDF.]

Figure 1. Decay rate of GAS5.

(A) Schematic representation of 5′-Bromo-uridine immunoprecipitation chase (BRIC). BrU is incorporated into newly transcribed RNA in place of uridine. Cells are collected and total RNA extracted at multiple time points after BrU washout. BrU-labeled RNAs are immunopurified using anti-BrdU beads. Non-labeled RNAs are discarded. Then, BrU-labeled RNAs are quantified by RT-qPCR. (B) The decay rate of GAS5 as determined by BRIC in HEK293T cells. Relative quantitative values at time 0 h were set to 100%. Values represent mean ± errors obtained from two independent experiments.

https://doi.org/10.1371/journal.pone.0055684.g001

GAS5 transcripts are degraded by the NMD pathway

To examine whether GAS5 transcripts are targeted for NMD, we determined the steady-state abundance and the decay rate of GAS5 with an abrogated NMD pathway by UPF1 depletion. Depletion of UPF1 by RNA interference (RNAi), using two different siRNAs to rule out off-target effects, reduced UPF1 mRNA levels in HEK293T cells to ∼28.6% of the control level (Figure 2A, S1). Knockdown of UPF1 increased the steady-state level of GAS5 by about ∼7.5 fold (Figure 2B, S1). This result is consistent with previous reports [25], [27]. Moreover, its half-life (t1/2) was prolonged from 5.6 h to 15.4 h in UPF-depleted cells (Figure 2C). The half-life of GAS5 in the control siRNA-treated cells (t1/2 = 5.6 h in Figure 2C) differed from that in normal cells (t1/2 = 6.6 h in Figure 1B). This difference is probably due to differences in cell culture conditions, e.g. addition of siRNAs and transfection reagents to siRNA-treated cells but not to normal cells. In contrast, knockdown of UPF1 did not affect the steady-state level or decay rate of MALAT1 [37]–[41], another nuclear long ncRNA (Figure 2B, 2C, S1). Moreover, the abundance of NMD-sensitive transcripts, GADD45A and SMG1, were significantly upregulated in response to UPF1 depletion, while those of non-NMD-sensitive transcripts, EIF4A3 and CASC3, were unchanged (Figure S1) [26]. These results validate UPF1 depletion for effective inhibition of the NMD pathway. Thus, we concluded that GAS5 is specifically degraded by the NMD pathway, and that NMD can regulate the abundance of a GAS5.

[Figure omitted. See PDF.]

Figure 2. Knockdown of UPF1 increases the expression level and prolongs the decay rate of GAS5, and decreases the expression of glucocorticoid-responsive genes.

HEK293T cells were treated with a control siRNA or with a siRNA targeting UPF1. (A) The expression levels of UPF1 in control cells (black bar) and in UPF1-depleted cells (grey bar) were determined by RT-qPCR. (B) The expression levels of GAS5 and MALAT1 in control cells (black bar) and in UPF1-depleted cells (grey bar) were determined by RT-qPCR. (C) The decay rates of GAS5 and MALAT1 were determined by BRIC in control cells (solid circle and black bar) and in UPF1-depleted cells (open circle and grey bar). Relative quantitative values at 0 h were set to 100%. (D) The expression levels of glucocorticoid-responsive genes, cIAP2 and SGK1, and the PPARδ-responsive gene, ADRP, in control cells (black bar) and in UPF1-depleted cells (grey bar) were determined by RT-qPCR. (A–D) The levels of GAPDH were used for normalization. Values represent mean±SD obtained from three independent experiments (*P<0.05, **P<0.01, Student's t test).

https://doi.org/10.1371/journal.pone.0055684.g002

The abundance of GAS5 modulates the transcriptional activity of apoptosis-related genes

GAS5 modulates the transcriptional activity of the glucocorticoid receptor and acts as a sensitizer of apoptosis by suppressing glucocorticoid-mediated induction of several responsive genes, including the cellular inhibitor of apoptosis 2 (cIAP2) and serum- and glucocorticoid-regulated kinase 1 (SGK1) [30]. The cIAP2 and the SGK1 are apoptosis inhibitors [42], [43]. The cIAP2 protein was found to bind caspases 3, 7, and 9 and inhibit their activities, thereby suppressing apoptosis triggered by treatment activation of the Fas receptor or the tumor necrosis factor–α (TNF–α) [42]. The cIAP2 gene promoter has imperfect tandem GREs ∼500 bp upstream of its transcriptional initiation site and GR activates its transcription by binding to these GREs [44], [45]. SGK1 also contains functional GREs in its promoter region [46]. To examine whether the mRNA levels of cIAP2 and SGK1 are decreased by the accumulation of GAS5, their mRNA levels were quantified in the cells where GAS5 was stabilized by the depletion of UPF1. GAS5 levels increased up to 7.5 fold in UPF1 knocked-down cells (Figure 2B, S1) and expectedly, the mRNA levels of cIAP2 and SGK1 were significantly decreased in the UPF1 knocked-down cells (Figure 2D, S1). In contrast, the mRNA level of the peroxisome proliferator–activated receptor δ (PPARδ)-responsive gene, adipose differentiation-related protein (ADRP) [47], whose transcription is not regulated by GAS5, was not significantly changed (Figure 2D, S1). We have also examined whether over-expression of GAS5 produces opposite effects to those of UPF1 knockdown. Consistent with the expectation from the knockdown experiments, GAS5 over-expression specifically down-regulated the expression levels of cIAP2 and SGK1 (Figure 3). These results suggest that the NMD pathway can regulate the levels of GAS5, resulting in a change in the transcriptional activity of apoptosis-related genes.

[Figure omitted. See PDF.]

Figure 3. Over-expression of GAS5 decreased the expression levels of glucocorticoid-responsive genes.

HEK293 cells were treated with plasmid vectors as indicated. The expression levels of the indicated RNAs in control cells with mock vector (black bar) and in cells over-expressing GAS5 (grey bar) were determined by RT-qPCR. The levels of GAPDH and ACTB were used for normalization. Values represent mean±SD obtained from three independent experiments (*P<0.05, **P<0.01, Student's t test).

https://doi.org/10.1371/journal.pone.0055684.g003

Serum starvation increases the expression level of GAS5 and decreases the GR-responsive genes

GAS5 is a 5′ TOP class gene [19]. Growth arrest by serum starvation is associated with attenuated translation of the 5′ TOP transcripts and inhibition of their degradation, resulting in accumulation of the GAS5 transcript [19]–[22]. To examine whether the mRNA levels of the GR-responsive genes, cIAP2 and SGK1, are decreased by the accumulation of GAS5 upon serum starvation, we cultured HEK293 cells in serum-free medium or in medium with a small amount of serum (0.5%) for up to 48 h to promote the accumulation of endogenous GAS5. GAS5 levels increased up to 3.7 fold in cells cultured for 48 h in serum-free medium, and 2.3 fold when cultured in medium containing 0.5% serum compared with cells cultured in the presence of 10% serum (Figure 4A). We have also determined the half-lives of GAS5 with or without serum starvation. Unfortunately, the BRIC method cannot be used for determining the half-life of GAS5 with serum starvation, because the recovery ratio of BrU-RNAs to non-labeled RNAs was too low. We suspect that BrU is not effectively incorporated into serum-starved cells. Instead, we used another inhibitor-free method, 5-Ethynyluridine (EU) pulse-labeling, which is commercially available as the Click-iT Nascent RNA Capture Kit (Invitrogen) [33]–[36]. We have successfully determined the half-life of GAS5 with or without serum starvation by the EU method, and found that t1/2 of GAS5 was prolonged from 4.2 h to >12 h in serum starved cells (Figure 4B). Similar to the UPF1-depletion experiments, the mRNA levels of cIAP2 and SGK1 were significantly decreased during serum starvation, and those of cIAP2 and SGK1 in medium containing 0.5% serum were slightly higher than those in serum-free medium (Figure 4C). Moreover, the mRNA level of the PPARδ-responsive gene, ADRP, was not significantly changed (Figure 4C). We have also examined whether knockdown of GAS5 suppresses the effects of serum withdrawal on the expression of glucocorticoid-responsive genes. Depletion of GAS5 by RNAi reduced GAS5 RNA levels in HEK293 cells to ∼21.9% of the control level (Figure 5A). The mRNA levels of cIAP2 and SGK1 were not significantly changed during serum starvation in GAS5-depleted cells (Figure 5B). These results indicate that the suppression of apoptosis-related gene expression in response to serum starvation is dependent on GAS5 transcripts.

[Figure omitted. See PDF.]

Figure 4. Serum starvation increases the expression level and prolongs the decay rate of GAS5, and decreases the expression of glucocorticoid-responsive genes.

HEK293 cells were cultured with 0, 0.5, and 10% serum as indicated. (A) The relative expression levels of GAS5 in cells with 0 and 0.5% serum relative to control cells (cultured with 10% serum) were determined by RT-qPCR. (B) The decay rates of GAS5 were determined using the Click-iT Nascent RNA Capture Kit (Invitrogen) in control cells (solid circle and black bar) and in serum-starved cells (open circle and grey bar). Relative quantitative values at 0 h were set to 100%. (C) The relative expression levels of glucocorticoid-responsive genes, cIAP2 and SGK1, and the PPARδ-responsive gene, ADRP, in cells with 0 and 0.5% serum relative to control cells (cultured with 10% serum) were determined by RT-qPCR. (A–C) The levels of GAPDH and ACTB were used for normalization. Values represent mean±SD obtained from two (in B) or three (in A and C) independent experiments (*P<0.05, **P<0.01, Student's t test).

https://doi.org/10.1371/journal.pone.0055684.g004

[Figure omitted. See PDF.]

Figure 5. Knockdown of GAS5 suppresses the effects of serum withdrawal on the expression of glucocorticoid-responsive genes.

HEK293 cells were treated with a control siRNA or with a siRNA targeting GAS5, and cultured either with or without serum for 48 h. (A, B) The relative expression levels of the indicated RNAs in cells with serum (black bar) or without serum (grey bar) relative to control (cells cultured at 0 h) were determined by RT-qPCR. The levels of GAPDH and ACTB were used for normalization. Values represent mean±SD obtained from three independent experiments (*P<0.05, **P<0.01, Student's t test).

https://doi.org/10.1371/journal.pone.0055684.g005

Discussion

In this study, we report the importance of stability-mediated ncRNA regulation using GAS5 as a model for long ncRNAs. The UPF1 knockdown experiments showed that GAS5 is specifically degraded by the NMD pathway, and that NMD can regulate the expression levels and the decay rates of GAS5, resulting in a change in the transcriptional activity of several apoptosis-related genes, such as cIAP2 and SGK1 (Figure 2). These effects were rescued by GAS5 over-expression (Figure 3). Moreover, we found that growth arrest by serum starvation elevates the levels of GAS5 and prolonged decay rates, resulting in decrease the levels of apoptosis-related gene mRNAs (Figure 4). These findings are supported by GAS5 knockdown, which suppress the effects of serum starvation on the expression of the apoptosis-related genes (Figure 5). NMD is inhibited by serum starvation because UPF1 phosphorylation is reduced [48]. NMD is also inhibited by other environmental insults, including amino acid starvation, hypoxia, ultraviolet irradiation, other DNA damage-inducing agents, and conditions that activate the unfolded protein response [24], [49]–[51]. Taken together, these observations strongly suggest that certain environmental insults, including serum starvation, inhibit the NMD pathway. This results in an accumulation of GAS5, which causes a suppression of apoptosis-related gene expression, via the GR. To our knowledge, this is the first clear demonstration that the RNA degradation pathway can regulate the function of a long ncRNA in mammalian cells.

Similar to the GAS5 ncRNA, 7SK RNA [52], [53] and steroid receptor RNA activator (SRA) RNA [54], [55] are also long ncRNAs that are involved in transcriptional regulation by binding with transcription factors. 7SK RNA (∼300 bases in human) is transcribed by RNA polymerase III and is located in the nucleus. Together with other associated cellular proteins including hexamethylene bisacetamide-induced 1/2 (HEXIM1/2), 7SK RNA inhibits the activity of the positive transcription elongation factor b (P-TEFb) [56], [57]. P-TEFb is required for RNA polymerase II transcription elongation. The structure of 7SK RNA is conserved across many species [58]. The long 5′ stem–loop that extends from nt 24 to 87 in 7SK RNA binds to the HEXIM1 protein [59], [60] as well as to the 3′ stem–loop from nt 296 to 331 [61]. 7SK RNA is released from P-TEFb/HEXIM/7SK complexes upon an arrest in transcription and in response to physiological stimulations, such as cardiac hypertrophy [62], inhibition of transcription by drugs or UV irradiation [56], [57], [63], which leads to P-TEFb activation.

SRA RNA was first identified in 1999 as a functional ncRNA able to co-activate steroid nuclear receptors [64]. SRA RNA is a long ncRNA that has multiple splice variants. The predominant SRA transcripts in normal tissue are approximately 0.7–0.9 kb while less abundant, larger, transcripts of 1.3–1.5 kb have been identified [54]. SRA RNA is a co-regulator not only of steroid and non-steroid nuclear receptors [64]–[68], but also of several other transcription factors [69], [70]. The SRA RNA folds to form multiple secondary structures, which contribute to its activity [64]. Over-expression of SRA RNA enhances transcription transactivation by steroid receptors. Conversely, knockdown of SRA RNA by siRNA suppressed androgen receptor activity in prostate cancer cells [71]. Many SRA-associated RNA-binding proteins have been found to either positively or negatively regulate the coactivator activity of SRA RNA [68], [69], [72]. Moreover, the SRA1 gene not only produces functional RNAs, but could also lead to the production of a protein [64], [73].

Recently, Yoon et al. reported that lincRNA-p21 stability is changed by the abundance of the RNA-binding protein, HuR, resulting in altered translational activity of a subset of mRNAs [74]. LincRNA-p21 (∼3.1 kb in human) was first noted in 2010 as a repressor in p53-dependent transcriptional responses [75]. In the presence of HuR, lincRNA-p21 is destabilized through the recruitment of let-7/Ago2. HuR then promotes the translation of CTNNB1 and JUNB mRNAs by favoring their association with polysomes. In the absence of HuR, lincRNA-p21 is stable and accumulates, and Rck promotes the association of lincRNA-p21 with CTNNB1 and JUNB mRNAs, repressing their translation through a mechanism that includes a reduction in the size of polysomes. Thus, HuR-dependent translation activation requires rapid degradation of lincRNA-p21 in order to prevent the recruitment of translation repressors onto target mRNAs. This report is particularly important because it shows that an RNA-binding protein can regulate the stability and function of a long ncRNA in mammalian cells.

In addition to the long ncRNAs mentioned above, recent studies have also revealed the functional involvement of long nuclear ncRNAs in the regulation of gene expression (e.g. XIST [76], HOTAIR [77], [78], CDKN2B-AS1/ANRIL [79], [80]) through the modification of chromatin, control of transcriptional and post-transcriptional gene expression (e.g. MALAT1 [37]–[41]), and maintenance of subnuclear structures (e.g. NEAT1/MENε/β [81]–[84]). Interestingly, most of these functional long ncRNAs have short life-spans [7], [17]. We believe that further studies are needed to elucidate the relationships between the RNA degradation pathway and the functions of nuclear long ncRNAs.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Knockdown of UPF1 increases the expression level of GAS5, and decreases the glucocorticoid-responsive genes. HEK293 cells were treated with a control siRNA or with a siRNA targeting UPF1. The expression levels of indicated genes in control cells (black bar) and in UPF1-depleted cells using UPF1 siRNA (2) (grey bar) was determined by RT-qPCR. The levels of GAPDH and ACTB were used for normalization. Values represent mean±SD obtained from three independent experiments (**P<0.01, student's test).

https://doi.org/10.1371/journal.pone.0055684.s001

(TIFF)

Author Contributions

Conceived and designed the experiments: HT NA. Performed the experiments: HT. Analyzed the data: HT NA. Wrote the paper: HT MT NA.

Citation: Tani H, Torimura M, Akimitsu N (2013) The RNA Degradation Pathway Regulates the Function of GAS5 a Non-Coding RNA in Mammalian Cells. PLoS ONE 8(1): e55684. https://doi.org/10.1371/journal.pone.0055684

References

1. Raghavan A, Bohjanen PR (2004) Microarray-based analyses of mRNA decay in the regulation of mammalian gene expression. Brief Funct Genomic Proteomic 3: 112–124.

2. Keene JD (2010) Minireview: global regulation and dynamics of ribonucleic Acid. Endocrinology 151: 1391–1397.

3. Schoenberg DR, Maquat LE (2012) Regulation of cytoplasmic mRNA decay. Nat Rev Genet 13: 246–259.

4. Bolognani F, Perrone-Bizzozero NI (2008) RNA-protein interactions and control of mRNA stability in neurons. J Neurosci Res 86: 481–489.

5. Sharova LV, Sharov AA, Nedorezov T, Piao Y, Shaik N, et al. (2009) Database for mRNA half-life of 19 977 genes obtained by DNA microarray analysis of pluripotent and differentiating mouse embryonic stem cells. DNA Res 16: 45–58.

6. Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, et al. (2011) Global quantification of mammalian gene expression control. Nature 473: 337–342.

7. Tani H, Mizutani R, Salam KA, Tano K, Ijiri K, et al. (2012) Genome-wide determination of RNA stability reveals hundreds of short-lived non-coding transcripts in mammals. Genome Res 22: 947–956.

8. Elkon R, Zlotorynski E, Zeller KI, Agami R (2010) Major role for mRNA stability in shaping the kinetics of gene induction. BMC Genomics 11: 259.

9. Rabani M, Levin JZ, Fan L, Adiconis X, Raychowdhury R, et al. (2011) Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells. Nat Biotechnol 29: 436–442.

10. Alonso CR (2012) A complex ‘mRNA degradation code’ controls gene expression during animal development. Trends Genet 28: 78–88.

11. Chekulaeva M, Filipowicz W (2009) Mechanisms of miRNA-mediated post-transcriptional regulation in animal cells. Curr Opin Cell Biol 21: 452–460.

12. Krol J, Loedige I, Filipowicz W (2010) The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet 11: 597–610.

13. Guttman M, Amit I, Garber M, French C, Lin MF, et al. (2009) Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature 458: 223–227.

14. Guttman M, Rinn JL (2012) Modular regulatory principles of large non-coding RNAs. Nature 482: 339–346.

15. Zhang Z, Qin YW, Brewer G, Jing Q (2012) MicroRNA degradation and turnover: regulating the regulators. Wiley Interdiscip Rev RNA 3: 593–600.

16. Rüegger S, Großhans H (2012) MicroRNA turnover: when, how, and why. Trends Biochem Sci 37: 436–446.

17. Clark MB, Johnston RL, Inostroza-Ponta M, Fox AH, Fortini E, et al. (2012) Genome-wide analysis of long noncoding RNA stability. Genome Res 22: 885–898.

18. Schneider C, King RM, Philipson L (1988) Genes specifically expressed at growth arrest of mammalian cells. Cell 54: 787–793.

19. Smith CM, Steitz JA (1998) Classification of gas5 as a multi-small-nucleolar-RNA (snoRNA) host gene and a member of the 5′-terminal oligopyrimidine gene family reveals common features of snoRNA host genes. Mol Cell Biol 18: 6897–6909.

20. Coccia EM, Cicala C, Charlesworth A, Ciccarelli C, Rossi GB, et al. (1992) Regulation and expression of a growth arrest-specific gene (gas5) during growth, differentiation, and development. Mol Cell Biol 12: 3514–3521.

21. Fleming JV, Hay SM, Harries DN, Rees WD (1998) Effects of nutrient deprivation and differentiation on the expression of growth-arrest genes (gas and gadd) in F9 embryonal carcinoma cells. Biochem J 330: 573–579.

22. Fontanier-Razzaq N, Harries DN, Hay SM, Rees WD (2002) Amino acid deficiency up-regulates specific mRNAs in murine embryonic cells. J Nutr 132: 2137–2142.

23. Friedel CC, Dölken L, Ruzsics Z, Koszinowski UH, Zimmer R (2009) Conserved principles of mammalian transcriptional regulation revealed by RNA half-life. Nucleic Acids Res 37: e115.

24. Mendell JT, Sharifi NA, Meyers JL, Martinez-Murillo F, Dietz HC (2004) Nonsense surveillance regulates expression of diverse classes of mammalian transcripts and mutes genomic noise. Nat Genet 36: 1073–1078.

25. Imamachi N, Tani H, Akimitsu N (2012) Up-frameshift protein 1 (UPF1): Multitalented entertainer in RNA decay. Drug Discov Ther 6: 55–61.

26. Tani H, Imamachi N, Salam KA, Mizutani R, Ijiri K, et al. (2012) Identification of hundreds of novel UPF1 target transcripts by direct determination of whole transcriptome stability. RNA Biol 9: 1370–1379.

27. Ideue T, Sasaki YT, Hagiwara M, Hirose T (2007) Introns play an essential role in splicing-dependent formation of the exon junction complex. Genes Dev 21: 1993–1998.

28. Mourtada-Maarabouni M, Hedge VL, Kirkham L, Farzaneh F, Williams GT (2008) Growth arrest in human T-cells is controlled by the non-coding RNA growth-arrest-specific transcript 5 (GAS5). J Cell Sci 121: 939–946.

29. Mourtada-Maarabouni M, Pickard MR, Hedge VL, Farzaneh F, Williams GT (2009) GAS5, a non-protein-coding RNA, controls apoptosis and is downregulated in breast cancer. Oncogene 28: 195–208.

30. Kino T, Hurt DE, Ichijo T, Nader N, Chrousos GP (2010) Noncoding RNA gas5 is a growth arrest- and starvation-associated repressor of the glucocorticoid receptor. Sci Signal 3: ra8.

31. Chen CY, Ezzeddine N, Shyu AB (2008) Messenger RNA half-life measurements in mammalian cells. Methods Enzymol 448: 335–357.

32. Blattner C, Kannouche P, Litfin M, Bender K, Rahmsdorf HJ, et al. (2000) UV-Induced stabilization of c-fos and other short-lived mRNAs. Mol Cell Biol 20: 3616–3625.

33. Tani H, Akimitsu N (2012) Genome-wide technology for determining RNA stability in mammalian cells: Historical perspective and recent advantages based on modified nucleotide labeling. RNA Biol 9: 1233–1238.

34. Jao CY, Salic A (2008) Exploring RNA transcription and turnover in vivo by using click chemistry. Proc Natl Acad Sci USA 105: 15779–15784.

35. Ideue T, Adachi S, Naganuma T, Tanigawa A, Natsume T, et al. (2012) U7 small nuclear ribonucleoprotein represses histone gene transcription in cell cycle-arrested cells. Proc Natl Acad Sci USA 109: 5693–5698.

36. Abe K, Ishigami T, Shyu AB, Ohno S, Umemura S, et al. (2012) Analysis of interferon-beta mRNA stability control after poly(I:C) stimulation using RNA metabolic labeling by ethynyluridine. Biochem Biophys Res Commun 428: 44–49.

37. Tani H, Nakamura Y, Ijiri K, Akimitsu N (2010) Stability of MALAT-1, a nuclear long non-coding RNA in mammalian cells, varies in various cancer cell. Drug Discov Ther 4: 235–239.

38. Bernard D, Prasanth KV, Tripathi V, Colasse S, Nakamura T, et al. (2010) A long nuclear-retained non-coding RNA regulates synaptogenesis by modulating gene expression. EMBO J 29: 3082–3093.

39. Tripathi V, Ellis JD, Shen Z, Song DY, Pan Q, et al. (2010) The nuclear-retained noncoding RNA MALAT1 regulates alternative splicing by modulating SR splicing factor phosphorylation. Mol Cell 39: 925–938.

40. Tano K, Mizuno R, Okada T, Rakwal R, Shibato J, et al. (2010) MALAT-1 enhances cell motility of lung adenocarcinoma cells by influencing the expression of motility-related genes. FEBS Lett 584: 4575–4580.

41. Miyagawa R, Tano K, Mizuno R, Nakamura Y, Ijiri K, et al. (2012) Identification of cis- and trans-acting factors involved in the localization of MALAT-1 noncoding RNA to nuclear speckles. RNA 18: 738–751.

42. Deveraux QL, Reed JC (1999) IAP family proteins--suppressors of apoptosis. Genes Dev 13: 239–252.

43. Mikosz CA, Brickley DR, Sharkey MS, Moran TW, Conzen SD (2001) Glucocorticoid receptor-mediated protection from apoptosis is associated with induction of the serine/threonine survival kinase gene, sgk-1. J Biol Chem 276: 16649–16654.

44. Hong SY, Yoon WH, Park JH, Kang SG, Ahn JH, et al. (2000) Involvement of two NF-kappa B binding elements in tumor necrosis factor alpha -, CD40-, and epstein-barr virus latent membrane protein 1-mediated induction of the cellular inhibitor of apoptosis protein 2 gene. J Biol Chem 275: 18022–18028.

45. Webster JC, Huber RM, Hanson RL, Collier PM, Haws TF, et al. (2002) Dexamethasone and tumor necrosis factor-alpha act together to induce the cellular inhibitor of apoptosis-2 gene and prevent apoptosis in a variety of cell types. Endocrinology 143: 3866–3874.

46. Itani OA, Liu KZ, Cornish KL, Campbell JR, Thomas CP (2002) Glucocorticoids stimulate human sgk1 gene expression by activation of a GRE in its 5′-flanking region. Am J Physiol Endocrinol Metab 283: E971–979.

47. Tachibana K, Kobayashi Y, Tanaka T, Tagami M, Sugiyama A, et al. (2005) Gene expression profiling of potential peroxisome proliferator-activated receptor (PPAR) target genes in human hepatoblastoma cell lines inducibly expressing different PPAR isoforms. Nucl Recept 3: 3.

48. Pal M, Ishigaki Y, Nagy E, Maquat LE (2001) Evidence that phosphorylation of human Upfl protein varies with intracellular location and is mediated by a wortmannin-sensitive and rapamycin-sensitive PI 3-kinase-related kinase signaling pathway. RNA 7: 5–15.

49. Sharifi NA, Dietz HC (2006) Physiologic substrates and functions for mammalian NMD. Nonsense-Mediated mRNA Decay. 2 97–109: 2006.

50. Gardner LB (2008) Hypoxic inhibition of nonsense-mediated RNA decay regulates gene expression and the integrated stress response. Mol Cell Biol 28: 3729–3741.

51. Gardner LB, Corn PG (2008) Hypoxic regulation of mRNA expression. Cell Cycle 7: 1916–1924.

52. Diribarne G, Bensaude O (2009) 7SK RNA, a non-coding RNA regulating P-TEFb, a general transcription factor. RNA Biol 6: 122–128.

53. Peterlin BM, Brogie JE, Price DH (2012) 7SK snRNA: a noncoding RNA that plays a major role in regulating eukaryotic transcription. Wiley Interdiscip Rev RNA 3: 92–103.

54. Colley SM, Leedman PJ (2011) Steroid Receptor RNA Activator - A nuclear receptor coregulator with multiple partners: Insights and challenges. Biochimie 93: 1966–1972.

55. Cooper C, Vincett D, Yan Y, Hamedani MK, Myal Y, et al. (2011) Steroid Receptor RNA Activator bi-faceted genetic system: Heads or Tails? Biochimie 93: 1973–1980.

56. Nguyen VT, Kiss T, Michels AA, Bensaude O (2001) 7SK small nuclear RNA binds to and inhibits the activity of CDK9/cyclin T complexes. Nature 414: 322–325.

57. Yang Z, Zhu Q, Luo K, Zhou Q (2001) The 7SK small nuclear RNA inhibits the CDK9/cyclin T1 kinase to control transcription. Nature 414: 317–322.

58. Marz M, Donath A, Verstraete N, Nguyen VT, Stadler PF, et al. (2009) Evolution of 7SK RNA and its protein partners in metazoa. Mol Biol Evol 26: 2821–2830.

59. Michels AA, Fraldi A, Li Q, Adamson TE, Bonnet F, et al. (2004) Binding of the 7SK snRNA turns the HEXIM1 protein into a P-TEFb (CDK9/cyclin T) inhibitor. EMBO J 23: 2608–2619.

60. Egloff S, Van Herreweghe E, Kiss T (2006) Regulation of polymerase II transcription by 7SK snRNA: two distinct RNA elements direct P-TEFb and HEXIM1 binding. Mol Cell Biol 26: 630–642.

61. Sedore SC, Byers SA, Biglione S, Price JP, Maury WJ, et al. (2007) Manipulation of P-TEFb control machinery by HIV: recruitment of P-TEFb from the large form by Tat and binding of HEXIM1 to TAR. Nucleic Acids Res 35: 4347–4358.

62. Kulkarni PA, Sano M, Schneider MD (2004) Phosphorylation of RNA polymerase II in cardiac hypertrophy: cell enlargement signals converge on cyclin T/Cdk9. Recent Prog Horm Res 59: 125–139.

63. Biglione S, Byers SA, Price JP, Nguyen VT, Bensaude O, et al. (2007) Inhibition of HIV-1 replication by P-TEFb inhibitors DRB, seliciclib and flavopiridol correlates with release of free P-TEFb from the large, inactive form of the complex. Retrovirology 4: 47.

64. Lanz RB, Razani B, Goldberg AD, O'Malley BW (2002) Distinct RNA motifs are important for coactivation of steroid hormone receptors by steroid receptor RNA activator (SRA). Proc Natl Acad Sci USA 99: 16081–16086.

65. Deblois G, Giguère V (2003) Ligand-independent coactivation of ERalpha AF-1 by steroid receptor RNA activator (SRA) via MAPK activation. J Steroid Biochem Mol Biol 85: 123–131.

66. Kawashima H, Takano H, Sugita S, Takahara Y, Sugimura K, et al. (2003) A novel steroid receptor co-activator protein (SRAP) as an alternative form of steroid receptor RNA-activator gene: expression in prostate cancer cells and enhancement of androgen receptor activity. Biochem J 369: 163–171.

67. Zhao X, Patton JR, Davis SL, Florence B, Ames SJ, et al. (2004) Regulation of nuclear receptor activity by a pseudouridine synthase through posttranscriptional modification of steroid receptor RNA activator. Mol Cell 15: 549–558.

68. Hatchell EC, Colley SM, Beveridge DJ, Epis MR, Stuart LM, et al. (2006) SLIRP, a small SRA binding protein, is a nuclear receptor corepressor. Mol Cell 22: 657–668.

69. Caretti G, Schiltz RL, Dilworth FJ, Di Padova M, Zhao P, et al. (2006) The RNA helicases p68/p72 and the noncoding RNA SRA are coregulators of MyoD and skeletal muscle differentiation. Dev Cell 11: 547–560.

70. Kelly VR, Xu B, Kuick R, Koenig RJ, Hammer GD (2010) Dax1 up-regulates Oct4 expression in mouse embryonic stem cells via LRH-1 and SRA. Mol Endocrinol 24: 2281–2291.

71. Kurisu T, Tanaka T, Ishii J, Matsumura K, Sugimura K, et al. (2006) Expression and function of human steroid receptor RNA activator in prostate cancer cells: role of endogenous hSRA protein in androgen receptor-mediated transcription. Prostate Cancer Prostatic Dis 9: 173–178.

72. Shi Y, Downes M, Xie W, Kao HY, Ordentlich P, et al. (2001) Sharp, an inducible cofactor that integrates nuclear receptor repression and activation. Genes Dev 15: 1140–1151.

73. Emberley E, Huang GJ, Hamedani MK, Czosnek A, Ali D, et al. (2003) Identification of new human coding steroid receptor RNA activator isoforms. Biochem Biophys Res Commun 301: 509–515.

74. Yoon JH, Abdelmohsen K, Srikantan S, Yang X, Martindale JL, et al. (2012) LincRNA-p21 suppresses target mRNA translation. Mol Cell 47: 648–655.

75. Huarte M, Guttman M, Feldser D, Garber M, Koziol MJ, et al. (2010) A large intergenic noncoding RNA induced by p53 mediates global gene repression in the p53 response. Cell 142: 409–419.

76. Erwin JA, Lee JT (2008) New twists in X-chromosome inactivation. Curr Opin Cell Biol 20: 349–355.

77. Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, et al. (2007) Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell 129: 1311–1323.

78. Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, et al. (2010) Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature 464: 1071–1076.

79. Yap KL, Li S, Muñoz-Cabello AM, Raguz S, Zeng L, et al. (2010) Molecular interplay of the noncoding RNA ANRIL and methylated histone H3 lysine 27 by polycomb CBX7 in transcriptional silencing of INK4a. Mol Cell 38: 662–674.

80. Kotake Y, Nakagawa T, Kitagawa K, Suzuki S, Liu N, et al. (2011) Long non-coding RNA ANRIL is required for the PRC2 recruitment to and silencing of p15(INK4B) tumor suppressor gene. Oncogene 30: 1956–1962.

81. Sasaki YT, Ideue T, Sano M, Mituyama T, Hirose T (2009) MENepsilon/beta noncoding RNAs are essential for structural integrity of nuclear paraspeckles. Proc Natl Acad Sci USA 106: 2525–2530.

82. Clemson CM, Hutchinson JN, Sara SA, Ensminger AW, Fox AH, et al. (2009) An architectural role for a nuclear noncoding RNA: NEAT1 RNA is essential for the structure of paraspeckles. Mol Cell 33: 717–726.

83. Sunwoo H, Dinger ME, Wilusz JE, Amaral PP, Mattick JS, et al. (2009) MEN epsilon/beta nuclear-retained non-coding RNAs are up-regulated upon muscle differentiation and are essential components of paraspeckles. Genome Res 19: 347–359.

84. Nakagawa S, Naganuma T, Shioi G, Hirose T (2011) Paraspeckles are subpopulation-specific nuclear bodies that are not essential in mice. J Cell Biol 193: 31–39.