|Home | About | Journals | Submit | Contact Us | Français|
RNA levels in a cell are regulated by the relative rates of RNA synthesis and decay. We recently developed a new approach for measuring both RNA synthesis and decay in a single experimental setting by biosynthetic labeling of newly transcribed RNA. Here, we show that this provides measurements of RNA half-lives from microarray data with a so far unreached accuracy. Based on such measurements of RNA half-lives for human B-cells and mouse fibroblasts, we identified conserved regulatory principles for a large number of biological processes. We show that different regulatory patterns between functionally similar proteins are characterized by differences in the half-life of the corresponding transcripts and can be identified by measuring RNA half-life. We identify more than 100 protein families which show such differential regulatory patterns in both species. Additionally, we provide strong evidence that the activity of protein complexes consisting of subunits with overall long transcript half-lives can be regulated by transcriptional regulation of individual key subunits with short-lived transcripts. Based on this observation, we predict more than 100 key regulatory subunits for human complexes of which 28% could be confirmed in mice (P < 10−9). Therefore, this atlas of transcript half-lives provides new fundamental insights into many cellular processes.
mRNA levels in a cell are determined by the relative rates of RNA synthesis by polymerases and degradation by nucleases. Constant transcript levels reflect an equilibrium of RNA synthesis and decay while changes in transcript levels may be caused by alterations in either of them (1). State-of-the-art gene expression profiling allows precise measurements of total transcript abundance on whole transcriptome level but cannot distinguish whether changes in total mRNA are due to alterations in de novo transcription or in decay. RNA decay rates have previously been determined by blocking transcription, e.g. using actinomycin D (act-D), and subsequently monitoring ongoing RNA decay over time (2–8). If RNA decay continues at the same rate after inhibition of transcription, decay rates for thousands of transcripts can be obtained. However, transcriptional arrest induces a major stress response in the cell. This influences key regulatory mechanisms governing RNA decay which leads to substantial stabilization of individual transcripts (9–12).
De novo transcription can be measured in a non-disruptive way by introducing 4-thiouridine (4sU) into newly transcribed RNA utilizing nucleoside salvage pathways (13) followed by thiol-mediated isolation of newly transcribed RNA from total RNA (13–18). By combining this technique with standard microarray techniques, newly transcribed RNA can be directly measured for thousands of genes at the same time (14,15,18). Furthermore, with the integrative approach we developed recently total cellular RNA can be separated into both newly transcribed, labeled RNA and pre-existing, unlabeled RNA with high specificity (14). RNA half-lives can then be determined based on both newly transcribed RNA/total RNA ratios as well as pre-existing RNA/total RNA ratios.
In this article, we demonstrate that half-life measurements based on RNA decay, e.g. after blocking transcription, are inherently imprecise for medium- to long-lived transcripts. In contrast, RNA half-lives determined from newly transcribed/total RNA ratios are precise independent of transcript half-life. In this study, we present the first atlas of RNA half-lives in human B-cells and murine fibroblasts determined with this superior precision. This atlas was used to identify patterns in transcript half-life conserved in mammals across species and cell types. A transcriptome-wide comparison between species revealed that transcript half-lives are conserved and specifically correlated to gene function. This enabled us to better characterize the regulation of important genes and identify conserved regulatory principles for a broad range of biological processes.
We show that differences in regulation between functionally similar proteins are reflected in differences in the corresponding transcript half-lives. As a consequence, such differences in transcript half-lives can be used to detect differential regulatory patterns for these genes. Most importantly, we provide strong evidence that the activity of protein complexes consisting of multiple subunits with overall long transcript half-lives can be both fast and efficiently regulated by transcriptional regulation of individual key subunits characterized by very short-lived transcripts. Based on this concept, we could identify more than 100 potential key regulators for protein complexes in human for which at least 28% were confirmed in mice. Accordingly, analysis of RNA half-life provides fundamental new insights into conserved regulatory mechanisms of many biological processes. Therefore, the atlas of transcript half-lives we provide in this study will be valuable for further studies on a large variety of biological processes.
Newly transcribed RNA was labeled in human B-cells (BL41) and murine NIH-3T3 fibroblasts by culturing cells in the presence of 4-thiouridine (4sU) for 1 h. Total cellular RNA was isolated, thiol-specifically biotinylated and separated into labeled, newly transcribed RNA and unlabeled, pre-existing RNA using streptavidin coated magnetic beads as described (14) (see also Supplementary Data for details). In addition, RNA decay rates were obtained for NIH-3T3 cells by blocking RNA synthesis for 1, 2 and 3 h using actinomycin D at a final concentration of 5 µg/ml. Three biological replicates per condition were analyzed using Affymetrix HG U133 Plus 2.0 arrays (human) and MG 430 2.0 arrays (mouse). For murine fibroblasts, an additional three replicates were performed for newly transcribed and total RNA to assess the reproducibility of our approach.
Microarray data were pre-processed with R and Bioconductor (19,20). A first normalization incorporated background correction, normalization and probe-level summarization by GCRMA. As these standard methods assume equal overall intensities for all arrays, a second normalization step is required to compensate for the different amounts of template mRNA present in newly transcribed RNA, pre-existing RNA and total RNA samples. In previous studies based on blocking transcription (e.g. using actinomycin D) this was either done by using reference genes (8) or by fitting the exponential decay model to time series measurements (5). We performed this normalization based on the combined analysis of total, newly transcribed and unlabeled pre-existing RNA from a single RNA sample (14). Since total RNA (N) is quantitatively separated into labeled, newly transcribed (L) and unlabeled, pre-existing RNA (U), normalized ratios of newly transcribed/total RNA and pre-existing/total RNA should add up to 100% and normalization factors can be obtained by a simple linear regression analysis (see Supplementary Data).
RNA decay has been shown to follow first-order kinetics (21) with
where λ is the decay rate for a given transcript. The transcript half-life then is t1/2 = ln 2/λ. At the beginning of labeling, U(0) = N(0). We assume that total RNA at time t is a multiple or fraction of the original concentration, i.e. N(t) = α (t)N(0) with α (t) a function of time. The RNA half-life of a specific probe set is then calculated as (see Supplementary Data for details)
where L(t)/N(t) and U(t)/N(t) are the normalized ratios of newly transcribed/total RNA and pre-existing/total RNA. The proportionality factor α (t) can be defined in different ways to model different scenarios: α (t) = 1 for the steady state and α (t) = 2t/CCL to model cell growth and division where CCL is the cell-cycle length of the cell. In this study, we used steady-state assumptions as reproducibility between replicates was higher than in the cell division model which amplifies measurement errors for long-lived transcripts (Supplementary Figure S1).
To calculate the number of uracils in the spliced transcript for each gene, cDNA sequences for human and mouse were downloaded from the Ensembl site (release 54, May 2009) (22). The uracil number was then calculated as the number of thymines in the cDNA sequence for each gene. For genes with alternatively spliced transcripts uracil numbers were averaged.
As genes may be represented by more than one probe set on the array, we defined a probe set quality score PQS based on the difference between 1 and the sum of normalized pre-existing/total RNA ratios and newly transcribed/total RNA ratios (for details see Supplementary Data):
Gene half-life for a gene g is determined using the probe set with the maximum PQS for gene g.
Differences in transcript half-life between two genes gi and gj were calculated as the RNA half-life ratio
where t1/2 (gi) is the RNA half-life for gene gi.
To identify functional groups significantly over-represented among short- or long-lived transcripts, we compared the overall distribution of RNA half-lives against the distribution of RNA half-lives for specific functional categories. For this purpose the functional categories of the Gene Ontology (GO) (23) were used and GO annotations were taken from the GO website. Only GO terms were analyzed with at least 10 annotated genes. Significance of differences in the distributions was determined with the Kolmogorov–Smirnov test (K–S test) in R (19). P-values were corrected for multiple testing using the method by Benjamini and Yekutieli (24), a more conservative version of the method of Benjamini and Hochberg (25), which controls the false discovery rate (FDR) and does not require the tests to be independent. Correction of P-values was performed for all ontologies taken together and statistically significant results were determined at a significance level of 0.001.
Protein family annotations were taken from the Pfam database (26) (downloaded 1 August 2008). In total, we obtained 3170 families for human and 3031 families for mouse. Using orthology mappings from the mouse genome database (MGD) (27), protein family members were mapped between species. After removing redundant families, we obtained a final list of 738 families consisting of at least two members with half-lives in both human and mouse. Average RNA half-life ratios for members of the same family or subunits of the same protein complex (see below) were compared against results for randomized half-lives or families/complexes (10 000 randomizations each). P-values were calculated as the fraction of randomizations with lower or equal average half-life ratios than the observed average. Transcripts were defined as fast-decaying (slow-decaying) if the corresponding half-life was among the 20% shortest (longest) half-lives in at least one cell line and the 40% shortest (longest) ones in the other cell line. Functional similarity between members of the same family were calculated using the relevance similarity measure defined by Schlicker et al. (28) which determines the similarity between GO annotations of two proteins. For our purposes, similarity was calculated separately for the molecular function and biological process ontologies of the GO.
Protein complexes for human and mouse (1185 and 285 protein complexes, respectively) were taken from the CORUM database (29) (downloaded 4 June 2008). Complexes were mapped between species and redundant complexes identical to another complex for either species were eliminated. This resulted in a large set containing 1434 non-redundant but partially overlapping protein complexes on which all analysis for both human and mouse were based.
Protein complex subunits with significantly shorter transcript half-life than the rest of the complex were identified by calculating for each subunit p in each complex C the difference ratio:
Here, hlr(p, C - p) is the average RNA half-life ratio between subunit p and all the other subunits of C and hlr(C - p, C - p) the average RNA half-life ratio between these other subunits. A subunit p was predicted as a regulatory subunit for a complex C if the following conditions were fulfilled: (i) the average RNA half-life ratio for p to the remaining subunits of C was at least 40% higher than the average half-life ratio between these subunits [dr(p, C) > 1.4]. (ii) RNA half-life was shorter than the median half-life in human B-cells and murine fibroblasts (~5 h) and shorter than the median transcript half-life in complex C.
To investigate RNA turnover rates in mouse and human, we analyzed RNA half-lives in human B-cells (BL41) and murine NIH-3T3 fibroblasts (see ‘Materials and Methods’ section). The murine measurements have recently been published (14). Human B-cells were chosen to compare the data from murine fibroblasts with a cell line of both a different species and cell type to identify conserved regulatory patterns. RNA half-lives were obtained based on both newly transcribed/total RNA and pre-existing/total RNA ratios (Figure 1A and B, and Supplementary Figure S2). In addition, we performed microarray measurements on total RNA following 1, 2 and 3 h of transcriptional arrest by actinomycin D (act-D) in murine fibroblasts (Supplementary Figure S3) in order to compare our new approach with this standard method used to determine RNA decay rates (2–8). While all three approaches provided highly reproducible data for short-lived transcripts (t½ < 1–2 h) only newly transcribed/total RNA ratios yielded reliable data for medium- to long-lived transcripts. Although reproducibility of half-lives increased with longer act-D treatment, differences between 2 and 3 h act-D treatment were still considerable (Supplementary Figure S3) indicating that transcriptional arrest would have to be prolonged by another several hours to obtain accurate transcript half-lives. By simulating the effect of noise on RNA half-life determination (Supplementary Figure S4), we confirmed that this is not a problem of the individual measurements but an inherent feature of RNA half-life measurements based on monitoring the decay of transcripts.
Efficient capture of nascent transcripts is dependent on the incorporation of sufficient amounts of 4sU. This is particularly important for short transcripts with low uracil content. Previously, efficiency of capture of nascent transcripts by the streptavidin-coated magnetic beads was shown to be very high as 80–90% of radioactively labeled nascent RNA could be recovered. Although so far this was not specifically evaluated for short transcripts, reduced capture rates of small transcripts may create a substantial bias when calculating RNA half-lives based on newly transcribed/total RNA ratios (30). To check for this kind of bias, we compared the uracil number for transcripts with the measured ratios of newly transcribed/total RNA for our microarray data. No significant correlation for human B-cells (Figure 1C) and only a weak correlation for murine fibroblasts (Supplementary Figure S2) was observed. Thus, the employed concentrations of 4sU were sufficient for highly efficient capture even of transcripts with low uracil content.
Accuracy and reproducibility of half-lives determined from newly transcribed/total RNA ratios can be further increased by assessing probe set quality based on the relationship between newly transcribed, pre-existing and total RNA. On both the Affymetrix MG 430 2.0 arrays (mouse) and the HG U133 Plus 2.0 arrays (human) many genes are represented by multiple probe sets (Supplementary Figure S5). Due to quality differences between probe sets and experimental noise, half-life measurements of different probe sets for a single gene often result in dramatically different results. We solved this problem by calculating a probe set quality score (PQS) for each probe set based on the difference between the sum of measurements for newly transcribed and pre-existing RNA and total RNA levels for each probe set (see ‘Materials and Methods’ section and Supplementary Data). To evaluate the performance of this procedure, we determined probe set quality scores independently for each replicate of total, newly transcribed and pre-existing RNA of human B-cells and murine fibroblasts and identified the optimal probe set for each gene for the corresponding replicate. In our study, the same probe sets were independently identified as optimal for the corresponding genes in all three replicates significantly more often than expected by chance (binomial test, FDR corrected P < 0.05). In addition, when considering only the optimal probe set for each gene in each replicate, we observed decreased variations between replicates and increased reproducibility of results compared to the standard averaging approach (Supplementary Figure S5). This indicates that our approach can identify individual, incorrect measurements but also distinct quality differences between probe sets in an experimental setting.
For both murine fibroblasts and human B-cells, we determined RNA half-lives for more than 8000 genes based on newly transcribed/total RNA ratios using only the optimal probe set for each gene (Supplementary Tables S1 and S2). Median RNA half-life t1/2m was determined at 315 min (95% confidence interval: 240–382 min) in human B-cells and 274 min in murine fibroblasts (225–323 min) (see ‘Materials and Methods’ section and Supplementary Data). Although median half-lives differ between the two cell lines, this difference is not significant (t-test for unequal variance, P ~0.38). Furthermore, the distribution of half-lives is similar (Figure 2A) and follows approximately a log-normal distribution in both species. Using orthology tables between human and mouse genes from the MGD (27), transcript half-lives were compared for 4825 genes with RNA half-lives for both species (see Figure 2B). The average half-life ratio (hlr, see ‘Materials and Methods’ section) between human and mouse was 1.8 and ~67% of genes were within the two-fold range. Although variation between species was significantly larger than between different replicates for the same species and experiment (hlr = 1.2–1.3), it was significantly lower than the variation between different replicates after 1 h Act-D treatment in murine fibroblasts (hlr = 1.97). In addition, differences in transcript half-life between human and mouse were reduced by filtering genes based on the variation between different replicates for the same species. The lower the half-life ratios between different replicates in one species, the lower were the half-life ratios between the two species (Figure 2C). Furthermore, for only 18 genes (~0.37%) the deviation between species was significant (t-test for unequal variance, FDR corrected P < 0.01). A list of these genes is provided in Supplementary Table S3.
Previous studies of RNA half-lives based on inhibition of transcription have shown that mRNAs of transcriptional genes are preferentially short-lived while transcripts involved in the metabolism of the cell are quite stable (5,8,31). Based on the RNA half-lives determined in this study, we performed a comprehensive analysis of GO terms (for detailed results see Supplementary Table S4 and S5) to associate functional categories with differences in half-life distribution (see ‘Materials and Methods’ section). In both species, short-lived transcripts were characteristic for genes involved in the regulation of transcription (P < 10−16) and signal transduction (P < 0.007) (Figure 3A and Supplementary Figure S6). Signal transduction has so far only been associated with fast transcript decay in Arabidopsis (5). Interestingly, short half-life was specific only for regulators of transcription and signal transduction as a comparison of the half-life distribution for regulator genes neither involved in transcription nor signal transduction against the distribution for all genes showed no significant difference. Furthermore, previous observations of fast transcript decay for apoptosis and cell cycle transcripts (8) could not be confirmed as no significant difference in the half-life distribution was observed (Supplementary Figure S6).
Significant enrichment for very long RNA half-lives was found for genes involved in cellular respiration and energy metabolism as well as translation and protein decay by the proteasome. Enzymes and protein complexes involved in all parts of energy metabolism consistently had half-lives of more than 5 h both in human and mouse (Supplementary Figure S6 and Table S6). Interestingly, our method revealed a behavior of translational genes not previously described. While the frequency of transcription and signal transduction regulators steadily decreased with increasing RNA half-life (Figure 3A) and the frequency of genes involved in metabolic processes like energy metabolism and the proteasome steadily increased (Figure 3B), translational genes, encoding e.g. for ribosomal subunits or translation initiation factors, clustered in the medium-to-long half-life range and showed a lower frequency on either side (Figure 3C).
Since protein family members generally have similar functions, we investigated whether this translates to similar transcript half-lives. Indeed, we found that similarity of transcript half-lives was significantly increased in protein families (P < 10−4). However, we also identified a large number of protein families with substantial variations in transcript half-life. 111 of 738 protein families (15%) with at least two members for which we obtained half-lives in both human and mouse contained family members with both fast- and slow-decaying transcripts (see ‘Materials and Methods’ section). A detailed list of these families containing both the genes with the shortest and longest transcript half-life conserved between mice and men is provided in Supplementary Table S7. The median size of these families  was significantly larger than for families without such large differences in transcript half-lives (5, Wilcoxon rank sum test P < 10−16). The most likely explanation for this finding is that these larger families simply contain more diverse family members with a wider range of functions. Indeed, we found that both similarity of molecular functions and biological processes was significantly lower than for the other families (Wilcoxon rank sum test P < 0.05, Supplementary Figure S7). Surprisingly however, we observed that the proteins with the largest difference in half-life did not necessarily show the largest differences in their molecular functions or biological processes. For 44 (39.6%) of the 111 families, functional similarity was actually higher when comparing proteins with short-lived transcripts to proteins with long-lived transcripts than when the other proteins in the family were compared. This indicates that these proteins have a similar function but differential types of regulation. A good example for this finding is the hexokinase gene family. While hexokinase I (HK-I) transcripts decayed slowly (t1/2 ~ 9 h) as do those of most other genes involved in cellular respiration, transcripts of hexokinase II (HK-II) decayed fast (t1/2 ~ 1–3.6 h) in both human and mouse. Phosphorylation of glucose by hexokinases is the first step of the glycolysis and hexokinase I is considered a ‘housekeeping gene’ whose mRNA levels remain stable despite alterations in glucose or insulin levels (32) or feeding conditions (33). Contrary to that, expression of HK-II is induced by a variety of stimuli (32,34–36) thereby accelerating hexose catabolism and regulating blood sugar levels. Fast changes in the expression of HK-II are supported by the fast turnover of HK-II mRNA. In contrast, slow decay of HK-I transcripts prevents any rapid changes in gene expression. Thus, transcript half-lives for these enzymes are specifically adjusted according to their functional role. We noted the same phenomenon for the family of cytosine triphosphate (CTP) synthases. Here, CTP synthase 1 (CTPS) has a short transcript half-life (t1/2 ~ 3 h) in both mice and men while CTPS 2 has a long transcript half-life (t1/2 > 11 h). Although CTP synthases are essential enzymes (37), little is known about functional differences and transcriptional regulation. From our results we predict a similar regulatory pattern as for the two hexokinases. One isoform (CTPS) can be rapidly and transiently induced by transcriptional regulation, e.g. during cell cycle, while RNA concentrations of the other isoform (CTPS2) are stable and provide basal enzyme activity levels.
Short- and long-lived transcripts were also observed for the BCL-2 family which contains both pro- and anti-apoptotic proteins (38) (Table 1). Interestingly, long RNA half-lives were only observed for pro-apoptotic family members—although not for all of them. This provides additional evidence that an arrest in transcription following severe stress conditions could lead to a selective decline of the short-lived anti-apoptotic genes but not the pro-apoptotic ones thereby promoting apoptosis (21). Interestingly, the pro-apoptotic BCL-2 family members BAX and BAK, which share a common domain structure and are generally assumed to substitute for each other (39) stand out with respect to their different transcript half-lives (>12 h versus 2–3 h, respectively) in both mice and men. There is evidence that BAX and BAK play non-redundant roles and are regulated in different ways (40–43). Our results support this non-redundant role indicating that it is facilitated by the differences in transcript turnover. This implies that transcriptional regulation is important for BAK-mediated regulation of apoptosis while BAX activity may be preferentially regulated by post-transcriptional means.
These examples show that analysis of transcript half-life can reveal different types of regulation of proteins which otherwise appear to have very similar function.
For yeast, it has been previously reported that decay rates of transcripts encoding subunits of the same protein complex are similar (31). To investigate this for human and mouse, we analyzed transcript half-lives for the 1434 non-redundant but partially overlapping human and mouse protein complexes taken from the CORUM database (29) (see ‘Materials and Methods’ section). In both human B-cells and murine fibroblasts, transcripts with long half-lives were significantly over-represented among transcripts encoding subunits of these protein complexes (P < 0.0003). As expected, this was particularly prominent for subunits of the large protein complexes involved in energy and protein metabolism (see Supplementary Table S6). In addition, we found that transcript half-lives for subunits of the same protein complex were significantly more similar to each other compared to random expectation (hlr = 1.8–2.0, P < 10−4) even if we accounted for the overall high transcript half-life in complexes by randomizing complex memberships instead of half-lives.
Interestingly, despite this relatively high similarity of RNA half-lives in complexes, individual subunits of some complexes deviated substantially in transcript half-life from the remaining subunits. In principle, a single protein may be involved in several different protein complexes which may be regulated in different ways and, thus, may be characterized by different median transcript half-lives. Therefore, for some subunits their deviation in transcript half-life to one complex might be explained by similarity to another. Nevertheless, even when excluding these proteins, we still identified 155 complexes in human and 164 in mouse (out of 698 and 650 complexes, respectively, for which we had half-lives for at least three subunits) which contained protein subunits with a considerably shorter half-life than the remaining subunits (Figure 4, see ‘Materials and Methods’ section for details). 61 (~37%) of the complexes with deviating subunits were identified in both species which are significantly more than expected by chance (hypergeometric test, P < 10−5). Notably, we found that not only these complexes but also their deviating subunits were conserved. In total, we identified 102 and 108 proteins in human and mouse, respectively, which showed significantly shorter half-lives than the other members of the complexes they were part of (Figure 4). A complete list of these proteins and the complexes they are involved in is provided in Supplementary Tables S8 and S9. For 29 (28%) of these proteins identified in human B-cells the significantly shorter RNA half-life in the corresponding complex could be confirmed in murine fibroblasts (hypergeometric test, P < 10−9). For an additional 26 (25%) proteins, RNA half-life in the murine fibroblasts was shorter than the median half-life of the complex but the difference was not sufficiently pronounced.
A short transcript half-life allows both fast up- and down-regulation of gene activity on the transcriptional level (see Supplementary Figure S8) (44). Our results suggest that transcriptional regulation of individual key subunits is an evolutionary conserved mechanism to regulate the activity of protein complexes despite overall long RNA half-lives. One example for this type of regulation is the PBAF (Polybromo- and BAF containing complex) chromatin remodeling complex. Here, only the ARID2 (AT rich interactive domain 2) protein is characterized by short transcript half-life (Figure 5A). This subunit has previously been shown to be essential for the stability of the PBAF complex (45). Due to its short transcript half-life, the activity of the PBAF complex can be efficiently regulated by transcriptional regulation of the ARID2 protein alone. Another example is the E3 ubiquitin ligase complex (Figure 5B) consisting of CUL5 (Cullin 5) and RNF7 (Ring finger protein 7) linked by the heterodimeric Elongin BC complex (TCEB1 and TCEB2) to an Ankyrin repeat and SOCS box (ASB) protein which serves as substrate-recognition component (46,47). Here, different ASB proteins are responsible for the recognition of different substrates. In our study, we found that two such ASB proteins, ASB6 and ASB7, showed significantly shorter transcript half-lives than the remaining subunits. As these proteins are responsible for substrate recognition, transcriptional regulation of only these subunits suffices to regulate the activity of the complex with regard to specific substrates. As all the other ligase subunits show long transcript half-lives, our results indicate that they are continuously available for binding with various ASB subunits targeting different substrates.
Regulation of biological processes occurs at the transcriptional, translational and post-translational level. Optimal control is only achieved by coordinated regulation at all levels. Yet, important information on functional characteristics on many biological processes can already be obtained by analyzing transcriptional regulation. While measurements of differential gene expression indicate which genes are regulated on the transcriptional level in a specific condition, we show in this study that analysis of RNA decay and turnover can provide insights on transcriptional regulation on a more general level.
RNA decay has been studied in a wide range of species: E. coli (3), yeast (31), Arabidopsis (4,5) and human (8). These studies were based on measurements of RNA decay after transcriptional arrest. In this article, we showed that these RNA half-lives, although quite accurate for short half-lives, are unreliable for medium to long half-lives. Contrary to that, measurements of RNA de novo transcription provide reliable and precise results on the whole range of RNA half-lives. Probe set quality scores determined for every probe set based on the combined analysis of newly transcribed, unlabeled pre-existing RNA and total cellular RNA further improved data quality.
One potential bias, which might affect half-life measurements based on newly transcribed RNA, is insufficient capture of short transcripts due to their low number of uracil residues. This would result in underestimation of newly transcribed/total RNA ratios and overestimation of corresponding half-lives. Such a bias was noted in a recent study by Miller et al. (30) which used 4-thiouracil (4tU) instead of 4-thiouridine (4sU) which we used in our study. By correlating uracil number of transcripts with newly transcribed/total RNA ratios we demonstrated that 4sU labeling for both human B-cells and murine fibroblasts resulted in sufficient 4sU incorporation to ensure efficient capture of transcript even with rather few (<100) uracil residues. Note that 4sU incorporation into nascent RNAs can be easily enhanced by increasing the applied 4sU concentration in the cell culture medium and, thus, transcript size bias can be experimentally controlled. In contrast, 4tU labeling requires the co-expression of uracil phosphoribosyltransferase (UPRT) of the protozoa Toxoplasma gondii (18). We found 4tU/UPRT based labeling to be strongly dependent on UPRT expression levels as well as the cell type under study but not on the concentration of 4tU (unpublished data). Therefore, labeling efficiency can not be significantly increased by simply adding more 4tU but transcript length bias needs to be controlled for by bioinformatic means (30).
Using our new approach, we determined precise RNA half-lives for more than 8000 genes in both human B-cells and mouse fibroblasts. By choosing two completely unrelated cell types, we focused on regulatory mechanisms not specific for only individual cell types. For about 5000 orthologous genes, we obtained RNA half-lives in both species and cell types. For the large majority of these orthologous genes, transcript half-lives are conserved across species and cell types. Only 18 out of the 4825 genes compared, i.e. only ~0.37%, actually showed a significant difference in transcript half-life between the two species. Furthermore, variation between species was correlated to the variation observed in the individual experiments for each species. This suggests that to a large degree the observed variations between species were due to variations within the individual experiments and do not constitute important inter-species differences. This does not imply that RNA decay is a static process and that no significant differences in transcript half-life exist in between these two cell lines or species. Our results only show that for conserved genes expressed in both cell lines, transcript half-life is also conserved.
Fast transcript decay allows rapid alterations of steady-state RNA concentrations due to transcriptional changes. At the same time, these changes can also be rapidly reversed. Thus, a short transcript half-life is important for efficient regulation at transcriptional level. Assuming that protein levels and transcript levels are correlated, protein levels of these genes can be efficiently regulated by alteration in transcription rates alone. In contrast, up- or down-regulation of stable transcripts takes a very long time to result in altered total RNA levels which, once established, also persist much longer. Consistent with previous reports, we confirmed a shift towards short half-lives for genes involved in the regulation of transcription and observed this also for regulators of signal transduction. However, a similar shift for genes involved in the regulation of cell cycle or apoptosis as proposed earlier (8) or regulating genes in general (apart from transcriptional and signal transduction regulators) could not be confirmed. Thus, a short transcript half-life is not characteristic for regulators as such but only for regulators of transcription and signal transduction. The most stable transcripts were found for genes encoding for energy metabolism and protein translation and degradation. Interestingly, we observed that RNA half-lives of genes involved in translation cluster in the medium- to long-lived range and decrease in frequency on either side. This indicates that a greater degree of transcriptional control may be required for constituents of the translational machinery than for transcripts coding for proteins involved in protein degradation and energy metabolism. As measurements of RNA decay rates based on transcriptional arrest are inherently imprecise for medium- to long-lived transcripts, it is not surprising that this has been missed by previous studies.
So far, the biological processes and sequence features determining RNA decay are only poorly understood. Previous studies have suggested that certain RNA motifs in untranslated regions (8) or miRNA binding and the presence of introns (5) may play a role. Our approach now allows the analysis of RNA sequence features and motifs which determine fast and slow but also intermediate fast RNA decay. Therefore, the method and data we provide in this study will be valuable for more systematic studies on the mechanisms governing RNA decay.
We identified many biological processes in which closely related members of the same protein family with overlapping function differ significantly in RNA half-life. Here, differences in transcript half-life likely correspond to differences in regulation and, accordingly, functional roles of the corresponding genes. This is best exemplified by hexokinase I and II as well as the pro-apoptotic proteins BAX and BAK. These examples show that transcript half-lives are fine-tuned to support the regulation of cooperative but non-redundant roles of closely related family members. Based on these results, we predict similar regulatory patterns and provide a database for a large number of functionally less characterized genes and processes.
Most proteins function by interacting with other proteins in protein complexes. In this study, we confirmed previous observations in yeast that transcript half-lives for subunits of protein complexes are very similar (31). Furthermore, decay of transcripts for these subunits was found to be generally slow. This implies that most protein complexes are pre-dominantly regulated at the post-translational level. Nevertheless, for more than 150 complexes with overall long transcript half-lives in both human and mouse we identified individual key subunits with a short transcript half-life which deviate significantly from the remaining subunits in all complexes they are part of. For almost a third of these proteins, we found this pattern to be conserved across species. The probability of finding the same complexes and subunits in both species by chance is negligibly small. Therefore, we propose a generalized mechanism employed by cells to facilitate regulation of protein complexes in an efficient and targeted way. For complexes depending on the availability of specific essential components, regulation of complex activity is accomplished by regulating the abundance of only one or few of these key subunits. Thereby, complex activity can be regulated both faster—as most subunits are available and may have already assembled—and more energy efficient than by regulating all complex members. With the two examples of the PBAF complex and the E3 ubiquitin ligase complex, we demonstrated how transcriptional regulation of individual key subunits which are e.g. critical for either complex formation and stability (for the PBAF complex) or specificity (for the E3 ubiquitin ligase complex) can support efficient regulation of complex activity at the transcriptional level. A similar observation was made by de Lichtenberg et al. (48) for protein complexes of the yeast cell cycle. As most of these complexes contained both periodically and constitutively expressed subunits, they suggested a ‘just-in-time assembly’ (instead of ‘just-in-time synthesis’) in which the timing of the complex assembly is regulated by transcriptional regulation of only some subunits. Our results indicate that this may not be specific for the cell cycle but a general mechanism by which the function of large protein complex is regulated.
Based on this concept, we predict altogether about 100 key regulatory subunits in more than 150 complexes for both human and mouse. This list can be further extended by 85 and 67 proteins in human and mouse, respectively, which show a significantly shorter RNA half-life in at least one complex but not all complexes they are contained in. In these cases, their short transcript half-life may be explained by the fact that they are part of a complex for which all subunits have to be regulated strongly on the transcriptional level and, accordingly, have short transcript half-lives. Although these subunits were not included in the predictions for key regulatory subunits, they probably also have an important regulative function within the other complexes they are part of. Further studies on the regulatory subunits we predict in this study are required for a better understanding of the regulation of the involved protein complexes and the biological processes they govern.
Supplementary Data are available at NAR Online.
German Federal Ministry of Education and Research (BMBF NGFNplus 01GS0801 to L.D, U.K, C.C.F, R.Z.); the Friedrich-Baur Stiftung (to L.D.). Funding for open access charge: Ludwig-Maximilians-Universität München and BMBF.
Conflict of interest statement. None declared.
We would like to thank Bernd Rädle for his excellent technical assistance.