Previous Article | Next Article ![]()
Molecular and Cellular Biology, May 2009, p. 2636-2643, Vol. 29, No. 10
0270-7306/09/$08.00+0 doi:10.1128/MCB.01506-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
CBD, UMR 5547, CNRS-UPS, 31062 Toulouse Cedex 4, France,1 Biomolecular Research Program, King Faisal Specialist Hospital and Research Center, Riyadh 11211, Saudi Arabia2
Received 26 September 2008/ Returned for modification 11 November 2008/ Accepted 23 February 2009
|
|
|---|
|
|
|---|
ARE-mediated decay (AMD) is regulated in response to intra- or extracellular cues. Molecular mechanisms underlying this regulation involve changes in the subcellular localization or posttranslational modifications of ARE-binding proteins (ARE-BPs), leading to the alteration of their ARE-binding activity and/or interaction with the degradation machinery (14, 18). While most ARE-BPs, such as AUF1, TTP, or KSRP, attract the degradation machinery, others, like the Hu/ELAV family of proteins, protect ARE mRNAs from rapid destruction (5, 6). An interesting finding is that ARE-BPs are evolutionarily conserved from yeast to mammals, as exemplified by Cth2, the yeast orthologue of TTP, which regulates the global expression of ARE genes in response to nutritional cues (26). Notably, recent studies point to the conservation of ARE-mediated gene regulation in invertebrates. In the leech Helobdella robusta, notch mRNA stability might be regulated by conserved AUUUA pentamers in its 3'UTR (20). In addition, mammalian AREs, such as those from tumor necrosis factor
and interleukin-6, induce reporter mRNA decay in Drosophila melanogaster cells (22), in which orthologues to both mammalian ARE-BPs and components of AMD pathways (e.g., 5' to 3' and 3' to 5') are conserved (16, 22). However, the genome-wide occurrence of AREs has not yet been investigated in invertebrates, so the potential contribution of AMD to global gene expression is unknown.
Using a computation-based approach, we have searched for AREs in the 3'UTRs of D. melanogaster genes and generated the first repertoire of Drosophila ARE genes. Combining both ex vivo and in vivo approaches, we have validated experimentally that our analysis has uncovered genuine AMD targets. Our analysis shows that AMD regulates temporal gene expression in D. melanogaster and predicts the widespread contribution of AREs to posttranscriptional regulation in this organism.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. The ARE motifs
|
(iii) ARE conservation between D. melanogaster and Drosophila pseudoobscura. D. pseudoobscura 3'UTR computationally predicted sequences assembled previously (34) were obtained from Alexander Stark (personal communication). Since the data set contains a complete set of one-to-one matched and aligned 3'UTR sequences for both D. melanogaster and D. pseudoobscura, we ran both 3'UTR sets through the ARED pipeline described above, producing an ARE transcript list for each species. The D. melanogaster list of ARE transcripts (data not shown) and the D-ARED (see Table S1 in the supplemental material) do not exactly match due to their assembling from different releases of Flybase. To assess the conservation of ARE motifs between the two species, we tabulated the joint distribution of non-ARE, class I, and class II ARE transcripts in D. melanogaster and D. pseudoobscura in a three-by-three contingency table (Table 2). The reported P value for a chi-squared test of independence (–log P > 15; all logs to base 10) indicates that D. melanogaster and D. pseudoobscura orthologues are extremely unlikely to contain an ARE by random chance alone. ARE conservation in orthologous branchless-RC (bnl-RC) and cecropinA1 (cecA1) 3'UTRs was assessed by performing sequence alignments using the UCSC Genome Bioinformatics Browser (http://genome.ucsc.edu/).
|
View this table: [in a new window] |
TABLE 2. D. melanogaster versus D. pseudoobscura cross-species ARE conservation analysisa
|
RNA analysis.
The expression of the GFP reporter mRNAs was induced in S2 or mbn-2 cells for 2 h with 0.7 mM of CuSO4. For the determination of mRNA half-life (t1/2), cells were incubated with actinomycin D (5 µg/ml) for the times indicated in the figures. Total RNA then was prepared using an RNeasy RNA isolation kit (Qiagen). For Northern blot analysis,
8 µg of total RNA was heat denatured and separated by electrophoresis on a nondenaturing 1% agarose gel in Tris-acetate-EDTA and electrotransferred to a nylon membrane (Hybond-N+; Amersham Pharmacia Biotech). Hybridization was carried out using in vitro-transcribed probes directed against GFP and rp49 mRNAs synthesized using 32P-labeled CTP. Total RNAs from whole flies or fat bodies were prepared with TRIzol reagent (Invitrogen). For quantitative reverse transcription-PCR (RT-PCR) analysis, 1 µg of total DNase-free RNA was retrotranscribed using Superscript II RT (Invitrogen). Real-time PCR was performed on cDNA obtained from 0.1 µg of total RNA using platinum SYBR green (Invitrogen) according to the manufacturer's instructions using an iCycler IQ real-time detection system (Bio-Rad). All samples were analyzed in triplicate. The levels of mRNAs were determined by cycling threshold analysis and were normalized using rp49 as a control. Normalized data then were used to quantify the relative levels of mRNA using the 
CT method (see Table S4 in the supplemental material for primers used for PCRs).
Transgenic flies. The cecA1 promoter region from –760 to +30 bp was amplified from the pA10 plasmid (28) and was used to drive the firefly luciferase coding sequence fused to the wild-type (WT) or mutant cecA1 3'UTR. The cecA1-luc-3'UTRcecA1 constructs were cloned into a pattB-white+ derivative plasmid (V. Daburon, unpublished data) to allow for phiC31-targeted site-specific recombination into the attP landing site 3R-86F on the Drosophila melanogaster genome (7). luciferase (luc) reporter plasmids (250 ng/µl) were injected into embryos, and transgenic lines were established. The UAS-Tis11 RNA interference (RNAi) stock is from the Vienna Drosophila RNAi Center.
Septic injury and in vivo luciferase monitoring. Five-day-old females were fed for 24 h with beetle luciferin (1 mM final concentration; Promega). Flies were collected under CO2 and individually infected by being pricked with a concentrated culture of E. coli, as described previously (8). Infected flies were placed in 96-well microtiter assay plates prefilled with luciferin-containing food as described previously (9). Plates were placed in a LucyI luminometer (Anthos) programmed for the automatic monitoring of relative luminescence units emanating from individual wells for 17 s each for 15 min.
|
|
|---|
16%) contain this motif in at least one transcript, and we compiled the results in a Drosophila ARE gene database (D-ARED) (see Table S1 in the supplemental material). The occurrence of this motif within 3'UTRs is significant, as 3'UTRs are 4- and 20-fold enriched in this motif compared to 5'UTRs and coding sequences, respectively. Interestingly, in analyzing the transcripts of 1,104 genes whose 3'UTR ends contained a poly(A) signal motif, we found that the mean length of ARE 3'UTRs is 2.5-fold longer (813 nucleotide [nt]; n = 1,104) than that of non-ARE 3'UTRs (323 nt; n = 5,745) (–log P > 6 by t test). These data suggest that shorter 3'UTRs could have been selected to escape ARE-mediated control and support the hypothesis that 3'UTR-mediated gene regulation impacts 3'UTR length (35). A closer analysis of the distribution of the AUUUA pentamer in 3'UTRs reveals that 352 genes (
20% of all ARE genes) contain at least two contiguous repeats of the AUUUA pentamer, a feature of mammalian class II ARE genes (13). Those 352 genes were further categorized into four groups (I to IV) (Table 1) depending on their number of AUUUA repeats (see Table S1 in the supplemental material). A number of conserved sequence motifs have previously been characterized by the alignment of 12 Drosophila genomes from species up to 45 million years distant (36). To test the functional significance of the mammalian consensus ARE sequence used to establish the D-ARED, we first examined whether the ARE sites overlap with positions of the evolutionarily conserved motifs in 3'UTRs. We found that 63.4% of the class II AREs in the D-ARED (or 53.9% of both class I and II AREs) overlap with conserved motifs, some of them containing the AUUUA core sequence (see Table S2 in the supplemental material). This result most probably is underestimated, since only AREs with strictly conserved positions are retained in this analysis, whereas no positional constraint in 3'UTRs has been reported for ARE activity in mammals or in yeast.
We then performed a direct cross-species analysis of ARE conservation between D. melanogaster and D. pseudoobscura using a previously assembled data set of computationally aligned 3'UTRs (10). We first found that AREs are similarly enriched in 3'UTRs of these two distantly related species (
25 million years) (29.8 and 24.9% of transcripts in D. pseudoobscura and D. melanogaster, respectively) (Table 2). In addition, as found in D. melanogaster, the mean length of non-ARE 3'UTRs in D. pseudoobscura (1,191 nt; n = 7,583) is significantly shorter than that of ARE-containing 3'UTRs (1,549 nt; n = 3,224). Finally and most remarkably, 52.1 and 66.8% of D. melanogaster class I and II ARE transcripts, respectively, exhibit a conserved ARE (class I or II) in orthologous 3'UTRs in D. pseudoobscura (Table 2). This result is highly significant (P
0 by a binomial test) compared to the expected random values of 30%, assuming factor independence (Table 2). Taken together, our data establish that AREs are widespread in D. melanogaster 3'UTRs and evolutionarily conserved across Drosophila species.
To assess the functional diversity of the ARE genes in our database, we determined gene ontology categories and analyzed the overrepresented molecular functions, biological processes, and pathways using the Panther database tool (http://www.pantherdb.org) (39). We found that D. melanogaster ARE genes are involved in a large variety of cellular and developmental processes, ranging from cell homeostasis to neurogenesis (see Table S3 in the supplemental material). We observed a marked enrichment of genes encoding transcription factors, as previously reported for ARE genes in human cells (4), where transcription factor-coding mRNAs are largely short lived (40). We also observed a significant enrichment of ARE genes in several receptor tyrosine kinase pathways (platelet-derived growth factor, epidermal growth factor, and fibroblast growth factor [FGF]) sharing the common RAS/mitogen-activated protein kinase signaling cascade and in the Wnt signaling pathway. Thus, ARE-mediated gene regulation could modulate the magnitude and/or duration of signaling pathway activation in Drosophila through the concerted regulation of several transduction components. It is noteworthy that a major effector in the Wnt/β-catenin signaling cascade is encoded by the β-catenin mRNA, whose stability is regulated by AMD in mammals (19). Our finding that D. melanogaster β-catenin mRNA contains a conserved ARE in its 3'UTR (see Tables S1 and S2 in the supplemental material) suggests that the investigation of the ARE-mediated regulation of signaling pathways in D. melanogaster will be relevant to related signaling pathways in mammals.
Drosophila AREs mediate mRNA decay in cultured cells. As contiguous repeats of the AUUUA pentamer are efficient destabilizing cis elements in mammals (17, 41), we have focused on class II ARE genes as promising AMD targets in Drosophila. RNA expression data were scarce and available for only 6 out of the 16 group 1 and 2 genes of this class. We selected branchless (bnl) for further analysis, as it exhibits a remarkable dynamic mRNA expression pattern in vivo that could reflect a well-orchestrated combination of transcription and stability controls (38). bnl codes for an FGF homologue playing a crucial role in the morphogenesis of the tracheal (respiratory) system (38). bnl produces two transcripts with distinct 3'UTRs, bnl-RA and bnl-RC (see Fig. S1 in the supplemental material). bnl-RC contains eight AUUUA pentamers in its 3'UTR, with four in tandem (see Table S1 in the supplemental material; also see Fig. 2A). bnl-RA is absent from the D-ARED even though it contains four dispersed AUUUA (see Fig. 2A), a pattern more related to the one present in class I ARE genes in mammals (13). D. melanogaster class I ARE genes (see Table S1 in the supplemental material) thus are likely to be underestimated in our database due to the high stringency of the ARE motif used for screening. In addition to bnl, we also selected cecA1 from group III for the further examination of ARE activity in living flies, since cecA1 mRNA decay is likely to be rapid (23) and can easily be assessed in vivo (see below). cecA1 codes for a small antimicrobial peptide (AMP) produced in response to bacterial and/or fungal infection (31). Its 3'UTR contains five AUUUA motifs, with two in tandem (see Table S1 in the supplemental material).
![]() View larger version (52K): [in a new window] |
FIG. 2. Drosophila AREs accelerate mRNA decay in cultured cells. (A) Schemes of the GFP reporters bearing the WT or the ARE mutant 3'UTRs of bnl-RC and bnl-RA and cecA1. AUUUA pentamers are indicated in boldface. AUUUA-containing regions (indicated above the constructs) have been entirely deleted ( ARE) or mutagenized at four AUUUA repeats into ACCCA (indicated by asterisks in AREmut). MT, metallothionein promoter. (B and C) Stability of the GFP reporters was assessed by Northern blotting (B), using probes specific for GFP and rp49 mRNAs, and by quantitative RT-PCR (qRT-PCR) (C). Reporters harboring bnl or cecA1 3'UTRs were transiently transfected and expressed in S2 or mbn-2 cells, respectively. Total RNA was isolated at different points following the addition of actinomycin D (5 µg/ml). Levels of reporter mRNA shown in panel C were normalized to those of rp49 mRNA, and averages and standard deviations from three independent experiments are plotted.
|
![]() View larger version (44K): [in a new window] |
FIG. 1. bnl-RC and cecA1 class II ARE conservation across Drosophila species. (A) Phylogenic tree relating the seven Drosophila species used for the evaluation of ARE conservation. (B) Sequences (in boldface) and groups of AREs (as defined in Table 1) are shown. The asterisk indicates the loss of class II ARE not converted into a class I ARE (group V).
|
ARE and 3'bnl RC
ARE) and cecA1 3'UTRs (3'cecA1AREmut), respectively (Fig. 2A). WT (control) and mutant 3'UTRs were fused downstream of the GFP coding sequence, and the mRNA decay of resulting reporters was analyzed in D. melanogaster cultured cells. We found that all ARE mutant reporters are stabilized compared to controls (Fig. 2B and C). This result shows that Drosophila AREs are bona fide cis regulators of mRNA decay. Notably, bnl-RA and bnl-RC 3'UTRs induce similar rates of reporter decay (Fig. 2C), suggesting that tandem and dispersed AUUUA repeats are equivalent destabilizing sequences under these conditions. However, it remains possible that different spacing and/or numbers of AUUUA pentamers elicit a range of regulation in vivo, depending on the tissue-specific expression and regulation of ARE-BPs. An independent work recently has proposed that as little as a single AUUUA pentamer can regulate gcm mRNA decay in D. melanogaster (33). Our present work establishes D-ARED as a novel and valuable predictive tool to identify genuine AMD target genes in D. melanogaster, together with available data on gene expression and ARE sequence conservation. Acute expression profile of cecA1 in cultured cells and in vivo relies on AMD. To validate that Drosophila AREs regulate gene expression in a physiological context, we examined the cecA1 ARE activity in LPS-stimulated cultured hemocytes (mbn-2), a cellular system allowing for the expression of immune-induced genes, including cecA1 (30). The cecA1 mRNA level increases and peaks 3 h after LPS treatment (inductive phase) and then rapidly declines (recovery phase) (Fig. 3A). We assessed cecA1 mRNA decay and found that its t1/2 drops from 136 to 54 min between the inductive and recovery phases, showing that transcription and mRNA turnover coordinately regulate cecA1 mRNA in response to LPS (Fig. 3B). To determine whether cecA1 AREs accelerates mRNA turnover during the recovery phase, the stability of both GFP-3'cecA1WT and GFP-3'cecA1AREmut reporters was analyzed 4 h following LPS treatment. The ARE mutant reporter is significantly stabilized (t1/2 = 289 min) compared to the control (t1/2 = 43 min) (Fig. 3C), strongly suggesting that AMD is responsible for the rapid decrease of cecA1 mRNA levels after an immune stimulation.
![]() View larger version (19K): [in a new window] |
FIG. 3. Rapid downregulation of endogenous cecA1 expression involves AMD. (A) Time course expression of cecA1 mRNA in mbn-2 cells exposed to LPS (10 µg/ml), as monitored by quantitative RT-PCR (qRT-PCR). The highest level of cecA1 mRNA was set as 100%. (B) Decay of cecA1 mRNA in mbn-2 cells following 1.5 or 4 h of LPS treatment. RNA was isolated at different points following the addition of actinomycin D and analyzed by qRT-PCR. (C) Decay of GFP-3'cecA1WT and GFP-3'cecA1AREmut mRNAs, as analyzed by qRT-PCR in mbn-2 cells following 4 h of LPS treatment. RNA was analyzed as described in the legend to Fig. 2B. In all panels, cecA1 mRNA levels were normalized to those of rp49 mRNA. Mean values ± standard deviations are given.
|
2-fold slower mRNA clearance of luc-3'cecA1AREmut (t1/2 = 118 min) compared to that of luc-3'cecA1WT (t1/2 = 60 min). We conclude that the ARE of cecA1 accelerates mRNA turnover in vivo and likely contributes to the appropriate temporal cecA1 turnoff during the immune response in D. melanogaster. It is noteworthy that the four cec genes of D. melanogaster (cecA1, cecA2, cecB, and cecC) (31) are found in the D-ARED (see Table S1 in the supplemental material) and that cecB mRNA also exhibits short-lived expression (30). Interestingly, a recent report by Ryu et al. (29) revealed that constitutive and high levels of AMP expression, including that of CecA1, are detrimental for the commensal microbial homeostasis in the gut of adult flies, resulting in pathological consequences. Taken together, these data indicate that AMD is essential to regulate AMP expression temporally during the Drosophila immune response, thus preventing tissue damage.
![]() View larger version (23K): [in a new window] |
FIG. 4. cecA1 ARE accelerates gene expression turnoff in vivo. (A) Schemes of the transgenic reporters inserted in the genome of D. melanogaster. Transgenes contain the luc cDNA fused to either the WT or mutagenized cecA1 3'UTR (as indicated in the legend to Fig. 2A). Inducible expression of luc is driven by cecA1 regulatory sequences in the fat body and hemocytes upon septic injury. (B) Time course analysis of bioluminescence in living 5-day-old females expressing luc transgenes after septic injury with E. coli. Shown is the mean luminescence (± standard deviations) from 10 individuals plotted as a function of time. For a given fly, the highest level of luminescence was set as 100%, and other luminescence levels were expressed relative to the peak expression value. As the time required for reaching the maximal level of luminescence (Tmax) varies between flies (between 5 and 7 h after infection), curves were arbitrarily aligned on their Tmax. (C) Quantification of luc mRNA decay in whole adult flies by quantitative RT-PCR. Total RNA was prepared from pools of 10 infected flies starting 6 h after infection. luc mRNA levels were normalized to those of rp49 mRNA. Mean values ± standard deviations from 3 independent experiments are shown.
|
![]() View larger version (23K): [in a new window] |
FIG. 5. Tis11-dependent AMD in vivo. (A) Semiquantitative (left) and quantitative (right) RT-PCR analysis of Tis11 mRNA silencing in mbn-2 cells treated with a specific Tis11 dsRNA. (B) Decay of cecA1 mRNA in untreated or Tis11-depleted mbn-2 cells analyzed by quantitative RT-PCR. Total RNA was prepared at the indicated time following 4 h of LPS exposure and the addition of actinomycin D. cecA1 mRNA levels were normalized to those of rp49 mRNA. Mean values ± standard deviations from three experiments are shown. (C) qRT-PCR analysis of Tis11 mRNA silencing in the dissected fat body of third-instar larvae carrying the Cg-Gal4 driver alone or together with the UAS-Tis11 RNAi transgene. Larvae were raised at 25°C 2 days before RNA isolation. Tis11 mRNA levels were normalized to those of rp49 mRNA and set arbitrarily as 100% in flies carrying the Gal4 driver alone. Mean values ± standard deviations are shown. (D) Prolonged cecA1 mRNA accumulation in Tis11 RNAi adult flies as analyzed by quantitative RT-PCR. Cg-Gal4 and Cg-Gal4/Tis11 RNAi transgenic females were raised for 5 days at 25°C before septic injury. Total RNA was prepared at the indicated time from pools of 10 infected flies starting 6 h after infection. cecA1 mRNA levels were normalized to those of rp49 mRNA. Mean values ± standard deviations are shown. (E) Time course analysis of bioluminescence in luc-3'cecA1WT reporter flies carrying only Cg-Gal4 drivers or together with UAS-Tis11 RNAi transgenes. Adults were raised for 5 days at 25°C before infection. The relative mean bioluminescence (± standard deviations) is represented as indicated in the legend to Fig. 4B.
|
3-fold slower decline of cecA1 mRNA after infection (Fig. 5D). To validate that Tis11 regulates cecA1 mRNA decay through its ARE-containing 3'UTR, the bioluminescence of WT and Tis11 RNAi flies expressing the luc-3'cecA1WT reporter was analyzed. We observed a significant prolonged bioluminescence of the reporter in Tis11 RNAi flies during the recovery phase (starting after the time to maximum level of bioluminescence [Tmax]) (Fig. 5E). Taken together, these data strongly suggest that Tis11 regulates cecA1 mRNA turnover in D. melanogaster through its ARE. This function is very reminiscent of the TTP-mediated downregulation of immune genes in mammals (1), and it will be of interest to further investigate whether Tis11 is a specific or general posttranscriptional regulator of immune genes in Drosophila.
AREs confer transient expression to early response genes during the inflammatory response in mammals (1). In D. melanogaster, many early immune-induced genes, including cec family genes, have a transient expression pattern (8). Remarkably, we found an
3-fold enrichment (–log P > 5 by a hypergeometric test) of genes containing the ARE motif WWWU[AUUUA]UWWW (W = A/U) among the Drosophila early immune-induced genes found previously (8). Thus, we propose that AMD contributes to terminate the expression of a wide variety of genes during D. melanogaster immune response, as described for mammals. Also tempting is the hypothesis that the conserved AUUUA pentamers in the 3'UTRs of D. melanogaster bnl and its mammalian orthologue FGF-10 (data not shown) contribute to the remarkable dynamic expression of those genes in vivo, which is essential to pattern the respiratory systems of both organisms (25). Thus, AMD might be conserved from Drosophila to humans to regulate physiological processes requiring the fine-tuning of gene expression. Due to its wide repertoire of genetic and molecular tools, D. melanogaster constitutes a powerful model system to investigate in vivo the relevance of ARE-mediated regulation.
F.C. is supported by the Fundação para a Ciência e a Tecnologia (SFRH/BPD/21460/2005) and the Fundação Calouste Gulbenkian.
Published ahead of print on 9 March 2009. ![]()
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»