Re-utilization of pyrimidine nucleotides during rat liver regeneration

Global SLAM-seq for accurate mRNA decay determination and identification of NMD targets

Gene expression analysis requires accurate measurements of global RNA degradation rates, earlier problematic with methods disruptive to cell physiology. Recently, metabolic RNA labeling emerged as an efficient and minimally invasive technique applied in mammalian cells. Here, we have adapted SH-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) for a global mRNA stability study in yeast using 4-thiouracil pulse-chase labeling. We assign high-confidence half-life estimates for 67.5% of expressed ORFs, and measure a median half-life of 9.4 min. For mRNAs where half-life estimates exist in the literature, their ranking order was in good agreement with previous data, indicating that SLAM-seq efficiently classifies stable and unstable transcripts. We then leveraged our yeast protocol to identify targets of the nonsense-mediated decay (NMD) pathway by measuring the change in RNA half-lives, instead of steady-state RNA level changes. With SLAM-seq, we assign 580 transcripts as putative NMD targets, based on their measured half-lives in wild-type and upf3Δ mutants. We find 225 novel targets, and observe a strong agreement with previous reports of NMD targets, 61.2% of our candidates being identified in previous studies. This indicates that SLAM-seq is a simpler and more economic method for global quantification of mRNA half-lives. Our adaptation for yeast yielded global quantitative measures of the NMD effect on transcript half-lives, high correlation with RNA half-lives measured previously with more technically challenging protocols, and identification of novel NMD regulated transcripts that escaped prior detection.

Metabolic turnover and dynamics of modified ribonucleosides by 13C labeling

Tandem mass spectrometry (MS/MS) is an accurate tool to assess modified ribonucleosides and their dynamics in mammalian cells. However, MS/MS quantification of lowly abundant modifications in non-ribosomal RNAs is unreliable, and the dynamic features of various modifications are poorly understood. Here, we developed a ¹³C labeling approach, called ¹³C-dynamods, to quantify the turnover of base modifications in newly transcribed RNA. This turnover-based approach helped to resolve mRNA from ncRNA modifications in purified RNA or free ribonucleoside samples and showed the distinct kinetics of the N6-methyladenosine (m⁶A) versus 7-methylguanosine (m⁷G) modification in polyA+-purified RNA. We uncovered that N6,N6-dimethyladenosine (m⁶₂A) exhibits distinct turnover in small RNAs and free ribonucleosides when compared to known m⁶₂A-modified large rRNAs. Finally, combined measurements of turnover and abundance of these modifications informed on the transcriptional versus posttranscriptional sensitivity of modified ncRNAs and mRNAs, respectively, to stress conditions. Thus, ¹³C-dynamods enables studies of the origin of modified RNAs at steady-state and subsequent dynamics under nonstationary conditions. These results open new directions to probe the presence and biological regulation of modifications in particular RNAs.

Thioguanosine Conversion Enables mRNA-Lifetime Evaluation by RNA Sequencing Using Double Metabolic Labeling (TUC-seq DUAL)

Temporal information about cellular RNA populations is essential to understand the functional roles of RNA. We have developed the hydrazine/NH₄ Cl/OsO₄ -based conversion of 6-thioguanosine (6sG) into A', where A' constitutes a 6-hydrazino purine derivative. A' retains the Watson-Crick base-pair mode and is efficiently decoded as adenosine in primer extension assays and in RNA sequencing. Because 6sG is applicable to metabolic labeling of freshly synthesized RNA and because the conversion chemistry is fully compatible with the conversion of the frequently used metabolic label 4-thiouridine (4sU) into C, the combination of both modified nucleosides in dual-labeling setups enables high accuracy measurements of RNA decay. This approach, termed TUC-seq DUAL, uses the two modified nucleosides in subsequent pulses and their simultaneous detection, enabling mRNA-lifetime evaluation with unprecedented precision.

Methods for the analysis of transcriptome dynamics

The transcriptome is the complete set of transcripts in a cell or tissue and includes ribosomal RNA (rRNA), messenger RNA (mRNA), transfer RNA (tRNA), and regulatory noncoding RNA. At steady-state, the transcriptome results from a compensatory variation of the transcription and decay rate to maintain the RNA concentration constant. RNA transcription constitutes the first stage in gene expression, and thus is a major and primary mode of gene expression control. Nevertheless, regulation of RNA decay is also a key factor in gene expression control, involving either selective RNA stabilization or enhanced degradation. Transcriptome analysis allows the identification of gene expression alterations, providing new insights regarding the pathways and mechanisms involved in physiological and pathological processes. Upon perturbation of cell homeostasis, rapid changes in gene expression are required to adapt to new conditions. Thus, to better understand the regulatory mechanisms associated with gene expression alterations, it is vital to acknowledge the relative contribution of RNA synthesis and decay to the transcriptome. To the toxicology field, the study of gene expression regulation mechanisms can help identify the early and mechanistic relevant cellular events associated with a particular response. This review aims to provide a critical comparison of the available methods used to analyze the contribution of RNA transcription and decay to gene expression dynamics. Notwithstanding, an integration of the data obtained is necessary to understand the entire repercussions of gene transcription changes at a system-level. Thus, a brief overview of the methods available for the integration and analysis of the data obtained from transcriptome analysis will also be provided.

Global analysis of RNA metabolism using bio-orthogonal labeling coupled with next-generation RNA sequencing

Many open questions in RNA biology relate to the kinetics of gene expression and the impact of RNA binding regulatory factors on processing or decay rates of particular transcripts. Steady state measurements of RNA abundance obtained from RNA-seq approaches are not able to separate the effects of transcription from those of RNA decay in the overall abundance of any given transcript, instead only giving information on the (presumed steady-state) abundances of transcripts. Through the combination of metabolic labeling and high-throughput sequencing, several groups have been able to measure both transcription rates and decay rates of the entire transcriptome of an organism in a single experiment. This review focuses on the methodology used to specifically measure RNA decay at a global level. By comparing and contrasting approaches and describing the experimental protocols in a modular manner, we intend to provide both experienced and new researchers to the field the ability to combine aspects of various protocols to fit the unique needs of biological questions not addressed by current methods.

Determination of in vivo RNA kinetics using RATE-seq

The abundance of a transcript is determined by its rate of synthesis and its rate of degradation; however, global methods for quantifying RNA abundance cannot distinguish variation in these two processes. Here, we introduce RNA approach to equilibrium sequencing (RATE-seq), which uses in vivo metabolic labeling of RNA and approach to equilibrium kinetics, to determine absolute RNA degradation and synthesis rates. RATE-seq does not disturb cellular physiology, uses straightforward normalization with exogenous spike-ins, and can be readily adapted for studies in most organisms. We demonstrate the use of RATE-seq to estimate genome-wide kinetic parameters for coding and noncoding transcripts in Saccharomyces cerevisiae.

Heterogeneous nuclear ribonucleoprotein (hnRNP) E1 binds to hnRNP A2 and inhibits translation of A2 response element mRNAs

Heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a trans-acting RNA-binding protein that mediates trafficking of RNAs containing the cis-acting A2 response element (A2RE). Previous work has shown that A2RE RNAs are transported to myelin in oligodendrocytes and to dendrites in neurons. hnRNP E1 is an RNA-binding protein that regulates translation of specific mRNAs. Here, we show by yeast two-hybrid analysis, in vivo and in vitro coimmunoprecipitation, in vitro cross-linking, and fluorescence correlation spectroscopy that hnRNP E1 binds to hnRNP A2 and is recruited to A2RE RNA in an hnRNP A2-dependent manner. hnRNP E1 is colocalized with hnRNP A2 and A2RE mRNA in granules in dendrites of oligodendrocytes. Overexpression of hnRNP E1 or microinjection of exogenous hnRNP E1 in neural cells inhibits translation of A2RE mRNA, but not of non-A2RE RNA. Excess hnRNP E1 added to an in vitro translation system reduces translation efficiency of A2RE mRNA, but not of nonA2RE RNA, in an hnRNP A2-dependent manner. These results are consistent with a model where hnRNP E1 recruited to A2RE RNA granules by binding to hnRNP A2 inhibits translation of A2RE RNA during granule transport.

Pyrimidine nucleotide metabolism in rat hepatocytes: evidence for compartmentation of nucleotide pools

Pyrimidine nucleotide metabolism in rat hepatocytes was studied by measurement of the labelling kinetics of the various intermediates after double labelling with [14C]orotic acid and [3H]cytidine, the precursors for the de novo and the salvage pathways respectively. For the uridine nucleotides, differences were found for the 14C/3H ratios in the UDP-sugars, in UMP (of RNA) and in their precursor UTP, suggesting the existence of separated flows of the radioactive precursors through the de novo and the salvage pathways. Higher ratios in the UDP-sugars, which are synthesized in the cytoplasm, and a lower ratio in UMP (of RNA) relative to the 14C/3H ratio in UTP indicated that UTP derived from orotic acid is preferentially used for the cytoplasmic biosynthesis of the UDP-sugars. Uridine, derived from cytidine, is preferentially used for the nuclear-localized synthesis of RNA. In contrast to these findings, the 14C/3H ratios in the cytidine derivatives CMP-NeuAc and CMP (of RNA), and in the liponucleotides CDP-choline and CDP-ethanolamine, were all lower than that in the precursor CTP. This indicates a preferential utilization of the salvage-derived CTP for the synthesis of the liponucleotides as well as for RNA and CMP-NeuAc. Similar conclusions could be drawn from experiments in which the intracellular amounts of several uridine- and cytidine-nucleotide-containing derivatives were increased by preincubating the hepatocytes with unlabelled pyrimidine nucleotides or ethanolamine. Based on these data, we propose a refined model for the intracellular compartmentation of pyrimidine nucleotide biosynthesis in which three pools of UTP are distinguished: a pool of de novo-derived molecules and a pool of salvage-derived molecules, both of which are channelled to the site of utilization; in addition an 'overflow' pool exists, consisting of molecules having escaped from channelling. An overflow pool could also be distinguished for CTP, but no discrimination between de novo and salvage-derived molecules could be made.

Half-life of cytoplasmic rRNA and tRNA, of plastid rRNA and of uridine nucleotides in heterotrophically and photoorganotrophically grown cells of Euglena gracilis and its apoplastic mutant W3BUL

1. For Euglena gracilis half-life data were calculated according to Greenberg's equation (Nature 240, 102-104, 1972) using steady-state specific radioactivities determined for cyt-rRNA, cyt-tRNA and pl-rRNA. 2. For all RNAs equal half-lives were found of 45 and 38 hr, respectively, in heterotrophically and photoorganotrophically grown cells. 3. Using the decay analysis equal half-lives were found for cyt-rRNA, cyt-tRNA and pl-rRNA being 79, 43 and 60 hr, respectively, in heterotrophically and photoorganotrophically grown wild-type cells and the mutant W3BUL. 4. As suggested by the specific radioactivity of intracellular [3H]UMP compared to that of [3H]uracil fed, the remarkable differences between RNA half-lives determined for heterotrophically and photoorganotrophically grown wild-type cells, seem to be caused by a different extent of the de novo synthesis of UMP. 5. Reutilization of RNA breakdown products suggested by increased half-lives of RNAs in the decay analysis compared to those determined by Greenberg's equation seems to occur mainly in heterotrophically grown cells.

Isotachophoretic and HPLC determination of nucleotides in rat liver cell nuclei isolated by non-aqueous technique

Rat liver whole cells and cell nuclei were prepared by a non-aqueous technique (glycerol). The nuclear preparations were of different purity as determined by RNA/DNA ratios (0.17-1.60) and accordingly were divided into 3 subgroups (mean values 0.29, 1.04 and 1.48). RNA nucleotides were separated by isotachophoresis and HPLC and calculated per mg DNA. Two of the nuclear subgroups (RNA/DNA = 1.04 and 1.48) had significantly elevated nucleotide values in relation to RNA/DNA. UDP-N-acetylhexosamine/DNA, on the contrary, was reduced in conformity with RNA in the preparations. Our findings may indicate different nucleotide concentrations in different parts of the cell.

Turnover of ribosomal proteins in regenerating rat liver after partial hepatectomy

The rates of synthesis and degradation of ribosomal proteins, prelabelled with [14C]bicarbonate, were determined as an index of the rate of ribosome turnover in regenerating rat liver. The half-life of ribosomes is about 178 and 75 hr in regenerating and normal liver, respectively. The comparison of turnover rates of ribosomal proteins with the corrected values of rRNA, based on re-utilization of nucleotides, suggests that ribosomes are metabolized as a unit in vivo. There is at least 70% overestimation for ribosome half-life when orotate-labelled RNA is used for turnover determinations. The absolute rate of synthesis is estimated as 3925 and 1081 ribosomes/min per cell in 24 hr regenerating and normal rat liver, respectively.