Messenger RNA in HeLa cells: kinetics of formation and decay

Time-resolved profiling of RNA binding proteins throughout the mRNA life cycle

mRNAs continually change their protein partners throughout their lifetimes, yet our understanding of mRNA-protein complex (mRNP) remodeling is limited by a lack of temporal data. Here, we present time-resolved mRNA interactome data by performing pulse metabolic labeling with photoactivatable ribonucleoside in human cells, UVA crosslinking, poly(A)+ RNA isolation, and mass spectrometry. This longitudinal approach allowed the quantification of over 700 RNA binding proteins (RBPs) across ten time points. Overall, the sequential order of mRNA binding aligns well with known functions, subcellular locations, and molecular interactions. However, we also observed RBPs with unexpected dynamics: the transcription-export (TREX) complex recruited posttranscriptionally after nuclear export factor 1 (NXF1) binding, challenging the current view of transcription-coupled mRNA export, and stress granule proteins prevalent in aged mRNPs, indicating roles in late stages of the mRNA life cycle. To systematically identify mRBPs with unknown functions, we employed machine learning to compare mRNA binding dynamics with Gene Ontology (GO) annotations. Our data can be explored at chronology.rna.snu.ac.kr.

High-throughput translational profiling with riboPLATE-seq

Protein synthesis is dysregulated in many diseases, but we lack a systems-level picture of how signaling molecules and RNA binding proteins interact with the translational machinery, largely due to technological limitations. Here we present riboPLATE-seq, a scalable method for generating paired libraries of ribosome-associated and total mRNA. As an extension of the PLATE-seq protocol, riboPLATE-seq utilizes barcoded primers for pooled library preparation, but additionally leverages anti-rRNA ribosome immunoprecipitation on whole polysomes to measure ribosome association (RA). We compare RA to its analogue in ribosome profiling and RNA sequencing, translation efficiency, and demonstrate both the performance of riboPLATE-seq and its utility in detecting translational alterations induced by specific inhibitors of protein kinases.

m6A mRNA modifications are deposited in nascent pre-mRNA and are not required for splicing but do specify cytoplasmic turnover

Understanding the biologic role of N⁶-methyladenosine (m⁶A) RNA modifications in mRNA requires an understanding of when and where in the life of a pre-mRNA transcript the modifications are made. We found that HeLa cell chromatin-associated nascent pre-mRNA (CA-RNA) contains many unspliced introns and m⁶A in exons but very rarely in introns. The m⁶A methylation is essentially completed upon the release of mRNA into the nucleoplasm. Furthermore, the content and location of each m⁶A modification in steady-state cytoplasmic mRNA are largely indistinguishable from those in the newly synthesized CA-RNA or nucleoplasmic mRNA. This result suggests that quantitatively little methylation or demethylation occurs in cytoplasmic mRNA. In addition, only ∼10% of m⁶As in CA-RNA are within 50 nucleotides of 5' or 3' splice sites, and the vast majority of exons harboring m⁶A in wild-type mouse stem cells is spliced the same in cells lacking the major m⁶A methyltransferase Mettl3. Both HeLa and mouse embryonic stem cell mRNAs harboring m⁶As have shorter half-lives, and thousands of these mRNAs have increased half-lives (twofold or more) in Mettl3 knockout cells compared with wild type. In summary, m⁶A is added to exons before or soon after exon definition in nascent pre-mRNA, and while m⁶A is not required for most splicing, its addition in the nascent transcript is a determinant of cytoplasmic mRNA stability.

Deciphering the Structure and Function of Nuclear Pores Using Single-Molecule Fluorescence Approaches

Due to its central role in macromolecular trafficking and nucleocytoplasmic information transfer, the nuclear pore complex (NPC) has been studied in great detail using a wide spectrum of methods. Consequently, many aspects of its architecture, general function, and role in the life cycle of a cell are well understood. Over the last decade, fluorescence microscopy methods have enabled the real-time visualization of single molecules interacting with and transiting through the NPC, allowing novel questions to be examined with nanometer precision. While initial single-molecule studies focused primarily on import pathways using permeabilized cells, it has recently proven feasible to investigate the export of mRNAs in living cells. Single-molecule assays can address questions that are difficult or impossible to answer by other means, yet the complexity of nucleocytoplasmic transport requires that interpretation be based on a firm genetic, biochemical, and structural foundation. Moreover, conceptually simple single-molecule experiments remain technically challenging, particularly with regard to signal intensity, signal-to-noise ratio, and the analysis of noise, stochasticity, and precision. We discuss nuclear transport issues recently addressed by single-molecule microscopy, evaluate the limits of existing assays and data, and identify open questions for future studies. We expect that single-molecule fluorescence approaches will continue to be applied to outstanding nucleocytoplasmic transport questions, and that the approaches developed for NPC studies are extendable to additional complex systems and pathways within cells.

Apoglobin Stability Is the Major Factor Governing both Cell-free and in Vivo Expression of Holomyoglobin

Expression levels in animal muscle tissues and in Escherichia coli vary widely for naturally occurring mammalian myoglobins (Mb). To explore this variation, we developed an in vitro transcription and wheat germ extract-based translation assay to examine quantitatively the factors that govern expression of holoMb. We constructed a library of naturally occurring Mbs from two terrestrial and four deep-diving aquatic mammals and three distal histidine mutants designed to enhance apoglobin stability but decrease hemin affinity. A strong linear correlation is observed between cell-free expression levels of holo-metMb variants and their corresponding apoglobin stabilities, which were measured independently by guanidine HCl-induced unfolding titrations using purified proteins. In contrast, there is little dependence of expression on hemin affinity. Our results confirm quantitatively that deep diving mammals have highly stable Mbs that express to higher levels in animal myocytes, E. coli, and the wheat germ cell-free system than Mbs from terrestrial mammals. Our theoretical analyses show that the rate of aggregation of unfolded apoMb is very large, and as a result, the key factor for high level expression of holoMb, and presumably other heme proteins, is an ultra high fraction of folded, native apoglobin that is capable of rapidly binding hemin. This fraction is determined by the overall equilibrium folding constant and not hemin affinity. These results also demonstrate that the cell-free transcription/translation system can be used as a high throughput platform to screen for apoglobin stability without the need to generate large amounts of protein for in vitro unfolding measurements.

Temporal and spatial characterization of nonsense-mediated mRNA decay

Nonsense-mediated mRNA decay (NMD) is a quality control mechanism responsible for "surveying" mRNAs during translation and degrading those that harbor a premature termination codon (PTC). Currently the intracellular spatial location of NMD and the kinetics of its decay step in mammalian cells are under debate. To address these issues, we used single-RNA fluorescent in situ hybridization (FISH) and measured the NMD of PTC-containing β-globin mRNA in intact single cells after the induction of β-globin gene transcription. This approach preserves temporal and spatial information of the NMD process, both of which would be lost in an ensemble study. We determined that decay of the majority of PTC-containing β-globin mRNA occurs soon after its export into the cytoplasm, with a half-life of <1 min; the remainder is degraded with a half-life of >12 h, similar to the half-life of normal PTC-free β-globin mRNA, indicating that it had evaded NMD. Importantly, NMD does not occur within the nucleoplasm, thus countering the long-debated idea of nuclear degradation of PTC-containing transcripts. We provide a spatial and temporal model for the biphasic decay of NMD targets.

Brazilin induces apoptosis and G2/M arrest via inactivation of histone deacetylase in multiple myeloma U266 cells

Although brazilin [7,11b-dihydrobenz(b)indeno[1,2-d]pyran-3,6a,9,10(6H)-tetrol] isolated from Caesalpinia sappan was known to have various biological activities, including anti-inflammation, antibacteria, and antiplatelet aggregation, there is no report yet on its anticancer activity. In the present study, the anticancer mechanism of brazilin was elucidated in human multiple myeloma U266 cells. We found that brazilin significantly inhibited the activity of histone deacetylases (HDACs), transcription factors involved in the regulation of apoptosis and cell cycle arrest in U266 cells. Consistently, brazilin enhanced acetylation of histone H3 at Lys 23, indicating activation of histone acetyltransferase (HAT), and also suppressed the expressions of HDAC1 and HDAC2 at both protein and mRNA levels. Additionally, brazilin significantly increased the number of sub-G1 cell population and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)-positive cells undergoing apoptosis and also activated caspase-3 and regulated the expression of Bcl-2 family proteins, including Bax, Bcl-x(L), and Bcl-2 in U266 cells, indicating that brazilin induces apoptosis through the mitochondria-dependent pathway. Interestingly, cell cycle analysis revealed that brazilin induced G2/M phase arrest along with apoptosis induction. Consistently, brazilin attenuated the expression of cyclin-dependent kinases (CDKs), such as cyclin D1, cyclin B1, and cyclin E, and also activated p21 and p27 in U266 cells. Furthermore, HAT inhibitor anacardic acid reversed activation of acetyl-histone H3 and cleavage of PARP induced by brazilin, while pan-caspase inhibitor Z-VAD-FMK001 did not affect the expression of HDAC induced by brazilin that brazilin mediates apoptosis via inactivation of HDAC in U266 cells. Notably, brazilin significantly potentiated the cytotoxic effect of standard chemotherapeutic agents, such as bortezomib or doxorubicin, in U266 cells. When our findings are taken together, they suggest that brazilin has potential as a chemotherapeutic agent alone or in combination with an anticancer agent for multiple myeloma treatment.

Eigen-genomic system dynamic-pattern analysis (ESDA): modeling mRNA degradation and self-regulation

High-throughput methods systematically measure the internal state of the entire cell, but powerful computational tools are needed to infer dynamics from their raw data. Therefore, we have developed a new computational method, Eigen-genomic System Dynamic-pattern Analysis (ESDA), which uses systems theory to infer dynamic parameters from a time series of gene expression measurements. As many genes are measured at a modest number of time points, estimation of the system matrix is underdetermined and traditional approaches for estimating dynamic parameters are ineffective; thus, ESDA uses the principle of dimensionality reduction to overcome the data imbalance. Since degradation rates are naturally confounded by self-regulation, our model estimates an effective degradation rate that is the difference between self-regulation and degradation. We demonstrate that ESDA is able to recover effective degradation rates with reasonable accuracy in simulation. We also apply ESDA to a budding yeast dataset, and find that effective degradation rates are normally slower than experimentally measured degradation rates. Our results suggest that either self-regulation is widespread in budding yeast and that self-promotion dominates self-inhibition, or that self-regulation may be rare and that experimental methods for measuring degradation rates based on transcription arrest may severely overestimate true degradation rates in healthy cells.

Actinomycin D Decreases Mcl-1 Expression and Acts Synergistically with ABT-737 against Small Cell Lung Cancer Cell Lines

Purpose: ABT-737, which blocks the function of Bcl-2 and Bcl-XL but not Mcl-1, has shown single-agent activity in preclinical models of small cell lung cancer (SCLC). Elevated expression of Mcl-1 induces resistance to ABT-737 in SCLC. Based on the short half-life of Mcl-1 mRNA and protein, we hypothesized that the actinomycin D could reverse Mcl-1–induced resistance to ABT-737.

Experimental Design: The dose-response of multiple SCLC cell lines to actinomycin D in the absence and presence of ABT-737 was followed by the assessment of Bcl-2 family expression and poly ADP ribose polymerase cleavage by Western blot, viability by tetrazolium dye reduction and clonogenic assay, and cell cycle kinetics by flow cytometry.

Results: Actinomycin D decreased Mcl-1 expression and resulted in a cell line–dependent increase in Noxa expression. Clinically relevant concentrations of actinomycin D from 0.4 to 4 ng/mL showed single-agent activity across a panel of SCLC cell lines. When combined with low micromolar doses of ABT-737, near complete loss of viability was seen with synergistic combination indices of 0.5 to 0.7. Exposure to 4 ng/mL actinomycin was only required for the first 24 hours of the combined incubation, mimicking a clinically achievable area under the curve, but the presence of ABT-737 was required for an additional 48 hours to obtain maximal effect.

Conclusions: Clinically relevant concentrations of actinomycin D act synergistically with ABT-737 to induce SCLC apoptosis, which can be at least partially attributed to the actinomycin D–induced decrease in Mcl-1 and increase in Noxa expression. Taken together, these data suggest the feasibility of combining actinomycin D with BH3-mimetic drugs in the clinical setting. Clin Cancer Res; 16(17); 4392–400. ©2010 AACR.

Kinetic analysis of RNA interference for lamin A/C in HeLa cells

We kinetically analyzed RNA interference (RNAi) for lamin A/C in HeLa cells, assuming suppression and recovery phases of gene expression, and dilution of transfected small interfering RNA (siRNA) by cell divisions. We observed the inhibitory effect of RNAi over a period of 6 days using various siRNA concentrations, and the maximum gene silencing efficiency occurred at Day 2 or 3. The gene silencing efficiency as a function of time and siRNA concentration was further quantitatively evaluated using a kinetic analysis method, demonstrating that RNAi for lamin A/C can be understood as a conventional drug-response system. This work provides potentially important guidelines for future applications using nanomaterials as delivery vehicles of siRNA in RNAi for lamin A/C.

Aberrant herpesvirus-induced polyadenylation correlates with cellular messenger RNA destruction

Inhibition of host cell gene expression by the human herpesvirus KSHV occurs via a novel mechanism involving polyadenylation-linked RNA turnover.

Regulation of messenger RNA (mRNA) stability plays critical roles in controlling gene expression, ensuring transcript fidelity, and allowing cells to respond to environmental cues. Unregulated enhancement of mRNA turnover could therefore dampen cellular responses to such signals. Indeed, several herpesviruses instigate widespread destruction of cellular mRNAs to block host gene expression and evade immune detection. Kaposi's sarcoma-associated herpesvirus (KSHV) promotes this phenotype via the activity of its viral SOX protein, although the mechanism of SOX-induced mRNA turnover has remained unknown, given its apparent lack of intrinsic ribonuclease activity. Here, we report that KSHV SOX stimulates cellular transcriptome turnover via a unique mechanism involving aberrant polyadenylation. Transcripts in SOX-expressing cells exhibit extended poly(A) polymerase II-generated poly(A) tails and polyadenylation-linked mRNA turnover. SOX-induced polyadenylation changes correlate with its RNA turnover function, and inhibition of poly(A) tail formation blocks SOX activity. Both nuclear and cytoplasmic poly(A) binding proteins are critical cellular cofactors for SOX function, the latter of which undergoes striking nuclear relocalization by SOX. SOX-induced mRNA turnover therefore represents both a novel mechanism of host shutoff as well as a new model system to probe the regulation of poly(A) tail-stimulated mRNA turnover in mammalian cells.

During viral infection, many essential cellular functions are targets for viral manipulation, yet aside from RNA interference, surprisingly few examples of viruses disrupting RNA turnover have been documented. Kaposi's sarcoma-associated herpesvirus (KSHV) is an oncogenic virus that induces widespread cellular messenger RNA destabilization during lytic infection. The viral protein SOX is a critical effector of this phenotype, yet it lacks ribonuclease activity, so presumably it targets cellular factors governing RNA stability. Here, we show that SOX stimulates host mRNA destruction via a unique mechanism involving polyadenylation. During SOX expression, newly formed messages have longer than normal poly(A) tails, leading to their retention in the nucleus. Coincident with this hyperadenylation, poly(A) binding protein (PABPC) is relocalized from the cytoplasm to the nucleus. PABPC has prominent roles in translation, messenger RNA stabilization, and quality control in the cytoplasm; we find its nuclear relocalization by SOX correlates with enhanced mRNA turnover in the cytoplasm. Thus, KSHV appears to have evolved distinct polyadenylation-linked mechanisms to target both new messages in the nucleus and preexisting cytoplasmic messages for destruction, thereby effectively inhibiting cellular gene expression.

Reconstruction of Escherichia coli transcriptional regulatory networks via regulon-based associations

Background

Network reconstruction methods that rely on covariance of expression of transcription regulators and their targets ignore the fact that transcription of regulators and their targets can be controlled differently and/or independently. Such oversight would result in many erroneous predictions. However, accurate prediction of gene regulatory interactions can be made possible through modeling and estimation of transcriptional activity of groups of co-regulated genes.

Results

Incomplete regulatory connectivity and expression data are used here to construct a consensus network of transcriptional regulation in Escherichia coli (E. coli). The network is updated via a covariance model describing the activity of gene sets controlled by common regulators. The proposed model-selection algorithm was used to annotate the likeliest regulatory interactions in E. coli on the basis of two independent sets of expression data, each containing many microarray experiments under a variety of conditions. The key regulatory predictions have been verified by an experiment and literature survey. In addition, the estimated activity profiles of transcription factors were used to describe their responses to environmental and genetic perturbations as well as drug treatments.

Conclusion

Information about transcriptional activity of documented co-regulated genes (a core regulon) should be sufficient for discovering new target genes, whose transcriptional activities significantly co-vary with the activity of the core regulon members. Our ability to derive a highly significant consensus network by applying the regulon-based approach to two very different data sets demonstrated the efficiency of this strategy. We believe that this approach can be used to reconstruct gene regulatory networks of other organisms for which partial sets of known interactions are available.

Influence of mRNA decay rates on the computational prediction of transcription rate profiles from gene expression profiles

The abundance of an mRNA species depends not only on the transcription rate at which it is produced, but also on its decay rate,which determines how quickly it is degraded. Both transcription rate and decay rate are important factors in regulating gene expression. With the advance of the age of genomics, there are a considerable number of gene expression datasets, in which the expression profiles of tens of thousands of genes are often non-uniformly sampled. Recently,numerous studies have proposed to infer the regulatory networks from expression profiles. Nevertheless, how mRNA decay rates affect the computational prediction of transcription rate profiles from expression profiles has not been well studied. To understand the influences, we present a systematic method based on a gene dynamic regulation model by taking mRNA decay rates, expression profiles and transcription profiles into account. Generally speaking,an expression profile can be regarded as a representation of a biological condition. The rationale behind the concept is that the biological condition is reflected in the changing of gene expression profile. Basically,the biological condition is either associated to the cell cycle or associated to the environmental stresses. The expression profiles of genes that belong to the former, so-called cell cycle data, are characterized by periodicity, whereas the expression profiles of genes that belong to the latter, so-called condition-specific data, are characterized by a steep change after a specific time without periodicity. In this paper, we examine the systematic method on the simulated expression data as well as the real expression data including yeast cell cycle data and condition-specific data (glucose-limitation data). The results indicate that mRNA decay rates do not significantly influence the computational prediction of transcription-rate profiles for cell cycle data. On the contrary,the magnitudes and shapes of transcription-rate profiles for condition specific data are significantly affected by mRNA decay rates. This analysis provides an opportunity for researchers to conduct future research on inferring regulatory networks computationally with available expression profiles under different biological conditions.

Rapid, diffusional shuttling of poly(A) RNA between nuclear speckles and the nucleoplasm

Speckles are nuclear bodies that contain pre-mRNA splicing factors and polyadenylated RNA. Because nuclear poly(A) RNA consists of both mRNA transcripts and nucleus-restricted RNAs, we tested whether poly(A) RNA in speckles is dynamic or rather an immobile, perhaps structural, component. Fluorescein-labeled oligo(dT) was introduced into HeLa cells stably expressing a red fluorescent protein chimera of the splicing factor SC35 and allowed to hybridize. Fluorescence correlation spectroscopy (FCS) showed that the mobility of the tagged poly(A) RNA was virtually identical in both speckles and at random nucleoplasmic sites. This same result was observed in photoactivation-tracking studies in which caged fluorescein-labeled oligo(dT) was used as hybridization probe, and the rate of movement away from either a speckle or nucleoplasmic site was monitored using digital imaging microscopy after photoactivation. Furthermore, the tagged poly(A) RNA was observed to rapidly distribute throughout the entire nucleoplasm and other speckles, regardless of whether the tracking observations were initiated in a speckle or the nucleoplasm. Finally, in both FCS and photoactivation-tracking studies, a temperature reduction from 37 to 22 degrees C had no discernible effect on the behavior of poly(A) RNA in either speckles or the nucleoplasm, strongly suggesting that its movement in and out of speckles does not require metabolic energy.

Decay rates of human mRNAs: correlation with functional characteristics and sequence attributes

Although mRNA decay rates are a key determinant of the steady-state concentration for any given mRNA species, relatively little is known, on a population level, about what factors influence turnover rates and how these rates are integrated into cellular decisions. We decided to measure mRNA decay rates in two human cell lines with high-density oligonucleotide arrays that enable the measurement of decay rates simultaneously for thousands of mRNA species. Using existing annotation and the Gene Ontology hierarchy of biological processes, we assign mRNAs to functional classes at various levels of resolution and compare the decay rate statistics between these classes. The results show statistically significant organizational principles in the variation of decay rates among functional classes. In particular, transcription factor mRNAs have increased average decay rates compared with other transcripts and are enriched in "fast-decaying" mRNAs with half-lives <2 h. In contrast, we find that mRNAs for biosynthetic proteins have decreased average decay rates and are deficient in fast-decaying mRNAs. Our analysis of data from a previously published study of Saccharomyces cerevisiae mRNA decay shows the same functional organization of decay rates, implying that it is a general organizational scheme for eukaryotes. Additionally, we investigated the dependence of decay rates on sequence composition, that is, the presence or absence of short mRNA motifs in various regions of the mRNA transcript. Our analysis recovers the positive correlation of mRNA decay with known AU-rich mRNA motifs, but we also uncover further short mRNA motifs that show statistically significant correlation with decay. However, we also note that none of these motifs are strong predictors of mRNA decay rate, indicating that the regulation of mRNA decay is more complex and may involve the cooperative binding of several RNA-binding proteins at different sites.

Stability of alternatively spliced IGF-I mRNA in growth plate chondrocytes

Insulin-like growth factor-I (IGF-I) enhances the proliferation and hypertrophy of growth plate chondrocytes leading to increased growth plate width thus promoting bone elongation. Differing quantities of the multiple IGF-I transcripts within the growth plate suggest differential regulation of IGF-I. To assess this, the relative stabilities of the 1Ea, 1Eb, and 2Ea IGF-I mRNA classes in 2-6 week old rat growth plates were evaluated. The mean estimated half-life of class 1Ea was 3.7+/-0.05 h, while classes 1Eb and 2Ea decayed gradually over the course of the experiment. Estimated half-lives for each IGF-I mRNA species were unchanged at all ages examined. Incubation with Act D enhanced the transcription of class 1Ea, 1Eb, and 2Ea mRNAs to varying degrees. This implies that the differential stability of alternative IGF-I mRNA classes may be an inherent regulatory component that is not influenced by an animal's developmental state, and the differential abundance of alternative IGF-I mRNA species in the growth plate throughout development is due to transcriptional differences.

Mechanisms and control of mRNA decapping in Saccharomyces cerevisiae

The process of mRNA turnover is a critical component of the regulation of gene expression. In the past few years a discrete set of pathways for the degradation of polyadenylated mRNAs in eukaryotic cells have been described. A major pathway of mRNA degradation in yeast occurs by deadenylation of the mRNA, which leads to a decapping reaction, thereby exposing the mRNA to rapid 5' to 3' exonucleolytic degradation. A critical step in this pathway is decapping, since it effectively terminates the existence of the mRNA and is the site of numerous control inputs. In this review, we discuss the properties of the decapping enzyme and how its activity is regulated to give rise to differential mRNA turnover. A key point is that decapping appears to be controlled by access of the enzyme to the cap structure in a competition with the translation initiation complex. Strikingly, several proteins required for mRNA decapping show interactions with the translation machinery and suggest possible mechanisms for the triggering of mRNA decapping.

Degradation of mRNA in Escherichia coli: an old problem with some new twists

Metabolic instability is a hallmark property of mRNAs in most if not all organisms and plays an essential role in facilitating rapid responses to regulatory cues. This article provides a critical examination of recent progress in the enzymology of mRNA decay in Escherichia coli, focusing on six major enzymes: RNase III, RNase E, polynucleotide phosphorylase, RNase II, poly(A) polymerase(s), and RNA helicase(s). The first major advance in our thinking about mechanisms of RNA decay has been catalyzed by the possibility that mRNA decay is orchestrated by a multicomponent mRNA-protein complex (the "degradosome"). The ramifications of this discovery are discussed and developed into mRNA decay models that integrate the properties of the ribonucleases and their associated proteins, the role of RNA structure in determining the susceptibility of an RNA to decay, and some of the known kinetic features of mRNA decay. These models propose that mRNA decay is a vectorial process initiated primarily at or near the 5' terminus of susceptible mRNAs and propagated by successive endonucleolytic cleavages catalyzed by RNase E in the degradosome. It seems likely that the degradosome can be tethered to its substrate, either physically or kinetically through a preference for monphosphorylated RNAs, accounting for the usual "all or none" nature of mRNA decay. A second recent advance in our thinking about mRNA decay is the rediscovery of polyadenylated mRNA in bacteria. Models are provided to account for the role of polyadenylation in facilitating the 3' exonucleolytic degradation of structured RNAs. Finally, we have reviewed the documented properties of several well-studied paradigms for mRNA decay in E. coli. We interpret the published data in light of our models and the properties of the degradosome. It seems likely that the study of mRNA decay is about to enter a phase in which research will focus on the structural basis for recognition of cleavage sites, on catalytic mechanisms, and on regulation of mRNA decay.

Glucocorticoids regulate pituitary growth hormone-releasing hormone receptor messenger ribonucleic acid expression

Glucocorticoids regulate GH synthesis and secretion by influencing both hypothalamic and pituitary function. With respect to GH-releasing hormone (GHRH), an important GH secretagogue, glucocorticoids are reported not only to suppress hypothalamic GHRH expression but also to augment pituitary responsiveness to GHRH. To investigate further this latter observation, we have determined the effects of this steroid on expression of the GHRH receptor (GHRH-R) gene in the rat pituitary in vivo and in pituitary cells in vitro. Adult male rats were adrenalectomized or sham operated and treated with s.c. implants of cholesterol or corticosterone. Adrenalectomized animals showed substantially reduced pituitary GHRH-R mRNA levels, when compared with untreated sham-operated animals. Conversely, administration of corticosterone increased pituitary GHRH-R mRNA levels in intact, as well as adrenalectomized rats. We also analyzed the effects of the synthetic glucocorticoid, dexamethasone, on GHRH-R mRNA expression in cultured rat anterior pituitary cells. GHRH-R mRNA was significantly increased by dexamethasone, with a maximal response observed in the presence of 100 nM hormone. This dose of dexamethasone substantially elevated GHRH-R mRNA after 6 h, 12 h, and 24 h of treatment. Dexamethasone did not increase GHRH-R mRNA in the presence of the transcriptional inhibitor actinomycin D, indicating that the predominant effect of the hormone is to increase transcription of the GHRH-R gene. These data demonstrate that GHRH-R mRNA levels are directly stimulated by glucocorticoids, both in the presence and absence of hypothalamic influences, providing a probable explanation for the ability of this steroid to alter pituitary responsiveness to GHRH.

mRNA stability in mammalian cells

This review concerns how cytoplasmic mRNA half-lives are regulated and how mRNA decay rates influence gene expression. mRNA stability influences gene expression in virtually all organisms, from bacteria to mammals, and the abundance of a particular mRNA can fluctuate manyfold following a change in the mRNA half-life, without any change in transcription. The processes that regulate mRNA half-lives can, in turn, affect how cells grow, differentiate, and respond to their environment. Three major questions are addressed. Which sequences in mRNAs determine their half-lives? Which enzymes degrade mRNAs? Which (trans-acting) factors regulate mRNA stability, and how do they function? The following specific topics are discussed: techniques for measuring eukaryotic mRNA stability and for calculating decay constants, mRNA decay pathways, mRNases, proteins that bind to sequences shared among many mRNAs [like poly(A)- and AU-rich-binding proteins] and proteins that bind to specific mRNAs (like the c-myc coding-region determinant-binding protein), how environmental factors like hormones and growth factors affect mRNA stability, and how translation and mRNA stability are linked. Some perspectives and predictions for future research directions are summarized at the end.

Rearrangements of intranuclear structures involved in RNA processing in response to adenovirus infection

We have studied in HeLa cells at the electron microscope level the response to adenovirus infection of the RNA processing machinery. Components of the spliceosomes were localized by in situ hybridization with biotinylated U1 and U2 DNA probes and by immunolabeling with Y12 anti-Sm monoclonal antibody, whereas poly(A)+ RNAs were localized by specific binding of biotinylated poly(dT) probe. At early stages of nuclear transformation, the distribution of small nuclear RNPs was similar to that previously described in non-infected nuclei (Visa, N., Puvion-Dutilleul, F., Bachellerie, J.P. and Puvion, E., Eur. J. Cell Biol. 60, 308–321, 1993; Visa, N., Puvion-Dutilleul, F., Harper, F., Bachellerie, J. P. and Puvion, E., Exp. Cell Res. 208, 19–34, 1993). As the infection progresses, the large virus-induced inclusion body consists of a central storage site of functionally inactive viral genomes surrounded by a peripheral shell formed by clusters of interchromatin granules, compact rings and a fibrillogranular network in which are embedded the viral single-stranded DNA accumulation sites. Spliceosome components and poly(A)+ RNAs were then exclusively detected over the clusters of interchromatin granules and the fibrillogranular network whereas the viral single-stranded DNA accumulation sites and compact rings remained unlabeled, thus appearing to not be directly involved in splicing. Our data, therefore, suggest that the fibrillogranular network, in addition to being the site of viral transcription, is also a major site of viral RNA splicing. Like the clusters of interchromatin granules, which had been already involved in spliceosome assembly, they could also have a role in the sorting of viral spliced polyadenylated mRNAs before export to the cytoplasm. The compact rings, which contain non-polyadenylated viral RNA, might accumulate the non-used portions of the viral transcripts resulting from differential poly(A)+ site selection.

Regional mRNA changes in brain stem motor neurons from patients with amyotrophic lateral sclerosis

An antisense oligonucleotide (54 mer) from the mRNA to the midsize neurofilament protein (NFM) was labeled with 35S on the 3' end and purified by polyacrylamide gel electrophoresis (PAGE). In situ hybridization was performed on sections from medulla oblongata of patients with amyotrophic lateral sclerosis (ALS). The slides were dipped in photographic emulsion, developed, and stained. Neurons from both nucleus hypoglossus and nucleus ambiguous showed a marked reduction of silver grains when compared to normal. This indicates a reduction of mRNA, which may precede the reduction of ribosomal RNA and the changes in neurofilament proteins that have been described by several investigators in ALS. It does not settle the question of whether the reduction of mRNA is owing to reduced transcription or increased decay of mRNA.

A comparison of apparent mRNA half-life using kinetic labeling techniques vs decay following administration of transcriptional inhibitors

Several different techniques were used to determine the apparent half-lives of immunoglobulin gamma 2b heavy chain and kappa light chain mRNA's in mouse myeloma 4T001 and a mutant derived from 4T001, i.e., mutant I17. The mutant I17 Ig heavy chain mRNA lacks CH1 and has fused CH2 and CH3 domains resulting in a truncated protein. By all four techniques the Ig heavy chain mRNA from mutant I17 displays a half-life that is approximately 70% the half-life of Ig mRNA in 4T001 cells. However, the absolute values of apparent half-life varied by greater than twofold for both lines among several of the techniques employed. The half-life of Ig gamma 2b mRNA in 4T001 cells was found to be 6.4 h by measuring decay following administration of the adenosine analog DRB to block new mRNA synthesis and 5.7 hr by measuring accumulation in an approach to steady-state labeling protocol. In contrast, the observed Ig mRNA half-lives determined by measuring decay following administration of actinomycin D to block new mRNA synthesis, or in a pulse-chase analysis were 2.9 and 3.8 h, respectively. The apparent half-life for Ig kappa light chain mRNA was the same in the 4T001 and I17 lines using any one technique but the value varied depending on the technique from a high value of 5.9 h following DRB to a low value of 2.4 h with actinomycin decay. Approach to steady-state is theoretically the most accurate method to measure mRNA half-life when that value is less than the doubling time of the cells. Pulse-chase analyses are accurate for measuring mRNA half-life when that value is longer than the effective chase period. Measuring preformed message decay following administration of drugs to block new mRNA synthesis is adaptable over a range of half-lives, but the cells must be shown to retain correct RNA metabolism over the time frame of the experiment. Determining a correct half-life for a particular mRNA may not be feasible using only one method and may, in fact, require several different approaches until a consensus value emerges.

Messenger RNA stabilization in chicken lens development: a reexamination

Previous studies, using actinomycin D, suggested that the rate of mRNA degradation in chicken embryo lens epithelial cells decreased sharply between Days 11 and 13 of development (Yoshida and Katoh, Exp. Eye Res. 11, 1971; Exp. Cell Res. 71, 1972). We repeated these studies and directly measured the effects of actinomycin on lens mRNAs. Although actinomycin decreased [3H]uridine incorporation greater than 90%, only modest decreases in methionine incorporation were detected. Treatment for 6 hr had no detectable effect on delta-crystallin mRNA levels or on the relative amounts of proteins synthesized in an in vitro translation system. Unlike the previous studies, we did not find significant changes in the stability of lens epithelial mRNAs during development.

Newly synthesized RNA: simultaneous measurement in intact cells of transcription rates and RNA stability of insulin-like growth factor I, actin, and albumin in growth hormone-stimulated hepatocytes

The levels of several RNA transcripts in cultured hepatocytes are regulated by transcriptional and post-transcriptional mechanisms and are affected by growth hormone and insulin. We assessed the effects of these hormones on transcription rates and the stability of insulin-like growth factor I, actin, and albumin transcripts in intact cells of primary cultures of rat hepatocytes by analyzing thiol-labeled, newly synthesized RNA isolated by mercurated agarose affinity chromatography. The application of this concept to the measurement of transcript stability is presented in detail. The data indicate that growth hormone stimulates the transcription rates of insulin-like growth factor I, actin, and albumin genes. The stability of all three transcripts, particularly albumin, appears to be lower in growth hormone-containing medium than it is in insulin-containing medium. The experiments indicate that the rates of transcription and/or degradation of albumin mRNA are influenced by hormonal treatment. However, the cells maintain roughly constant albumin transcript levels independent of hormone treatment by compensatory changes in the rates of transcription and degradation.

Increased half-life of mu immunoglobulin mRNA during mouse B cell development increases its abundancy

When B cells encounter antigen, the cells mature into terminally differentiated plasma cells and the amount of steady-state immunoglobulin (Ig) mu mRNA is increased 23-60-fold over the amount seen in earlier B cell stages. Most of this dramatic increase in Ig gene mRNA accumulation could be due to post-transcriptional regulation. We have treated a series of mouse cell lines fixed at different stages of B cell differentiation with an adenosine nucleotide analog 5,6-dichloro-1-beta-D-ribofuranosylbenzimidazole (DRB) which specifically blocks synthesis of new RNA polymerase II transcripts. The amount of mu heavy chain cytoplasmic RNA, measured by quantitative Northern blot analysis at various times post DRB treatment, is reflective of the transcript's stability. The mu mRNA half-life values observed from the earliest-stage lymphomas (70Z/3 and WEHI-231) are about 1.9-4 hr, whereas the t1/2 of mu mRNA in the hybridomas (Hyb54.3C2 and IdG11) is about 13-17 hr. There is, therefore, a nine-fold maximal increase in half-life of the mu mRNA in the Hyb54.3C2 over that observed in the earliest stage (70Z/3) cells.

Translation and a 42-nucleotide segment within the coding region of the mRNA encoded by the MAT alpha 1 gene are involved in promoting rapid mRNA decay in yeast

In yeast, the mRNA encoded by the MAT alpha 1 gene is unstable (t1/2 = 5 min) and the mRNAs encoded by the ACT1 gene (t1/2 = 30 min) and the PGK1 gene (t1/2 = 45 min) are stable. To understand the RNA structural features that dictate mRNA decay rates in yeast, we have constructed PGK1/MAT alpha 1 and ACT1/MAT alpha 1 gene fusions and analyzed the decay rates of the resultant chimeric transcripts. Fusion of a MAT alpha 1 segment containing 73% of the coding region and the 3' untranslated region to either of the stable genes is sufficient to cause rapid decay of the chimeric mRNAs (t1/2 = 6-7.5 min). Sequences required for this rapid decay are not found in the MAT alpha 1 3' untranslated region but are located within a 42-nucleotide segment of the coding region that has a high content (8 out of 14) of rare codons. Introduction of a translational stop codon upstream of this region stabilizes the hybrid mRNAs, indicating that the rapid decay promoted by these sequences is dependent on ribosomal translocation.

Characterization of an inducible promoter system to investigate decay of stable mRNA molecules

We have developed a system in which the decay of stable mRNAs can be studied without the use of inhibitors of transcription. The Drosophila hsp70 heat shock promoter linked to the bovine growth hormone (BGH) gene was used to establish stable cell lines in which the BGH gene is transcribed in a conditional manner. The BGH mRNA is synthesized only after induction at 43 degrees C. Following a brief period of re-equilibration at 37 degrees C during which transcription of the heat shock-driven gene ceases, the stable BGH mRNA decays with typical first-order kinetics. Hence, the decay of the mRNA can be studied without assumptions regarding radioactive labeling of precursor pools or transcriptional inhibitors. The system is applicable to any stable mRNA.

A 6-fold difference in the half-life of immunoglobulin mu heavy chain mRNA in cell lines representing two stages of B cell differentiation

When B cells differentiate into plasma cells, there is a large increase in the cellular content of mRNA coding for immunoglobulin. This increase cannot be fully accounted for by the increase in rate of transcription of the genes. We have investigated the possibility that the half-life of mu heavy chain mRNA increases during B cell differentiation, by measuring the rates of decay of the same endogenous mu gene in two cell lines representing the B cell and the plasma cell. Using a pulse-chase protocol, it was found that there was a significant increase in the half-life of mu mRNA between the B cell and the plasma cell, and no detectable difference in the average half-life of total poly(A)+ RNA in the two cell lines. The reduced rate of decay of mu mRNA in the more differentiated cell type is almost sufficient to account for the difference in steady state mu mRNA levels between the two cell lines.

Regulation of microM vs microS mRNA expression in an inducible B cell line

During the course of B lymphocyte differentiation into immunoglobulin secreting cells the abundance of mRNA for the heavy chain of secreted IgM (microS) increases dramatically. In order to understand the regulatory events responsible for the selective increase in micS mRNA we have looked for transcriptional alterations of VDJC mu gene segments as well as changes in the relative stability of microM and microS mRNA in BCL1 lymphoma cells which can be stimulated to increase the expression of microS mRNA. These experiments showed that although the transcriptional level of the mu gene locus is not preferentially augmented after stimulation, an alteration in the sites of polymerase termination is a significant factor contributing to the higher microS to microM ratio. This switch is dependent on new RNA synthesis. In addition, although the half-life of microS mRNA is not selectively increased, stimulation of the cells does result in a specific enhancement of the half-lives of both species of mu mRNA, which accounts for the higher steady state levels of total mu message.

High constitutive transcription of HSP86 gene in murine embryonal carcinoma cells

In order to investigate HSP86 heat-shock gene expression in embryonal carcinoma cell lines (EC), a partial mouse HSP86 cDNA clone was isolated and characterized. As observed for the corresponding protein, HSP86 RNA is shown to be constitutively more abundant in PCC4 and undifferentiated F9 EC cells than in fibroblasts, while its amount decreases upon F9 differentiation. Although mRNA stabilization is suggested to account in part of the high constitutive expression of the heat-shock-like protein HSC73 in F9 cells, HSP86 RNA appears as stable in fibroblasts as in F9 cells. Using run-on experiments we have established that high HSP86 expression in undifferentiated F9 cells in mainly due to enhanced transcription of the gene. Possible mechanisms responsible for this high level of transcription are discussed.

Highly localized tracks of specific transcripts within interphase nuclei visualized by in situ hybridization

Use of in situ hybridization optimized for fluorescent detection of nuclear RNA has revealed a striking localization of specific viral RNAs within nuclei of cells latently infected with EBV. Several hundred kb of specific transcripts is sharply restricted to a small region of the nucleus, frequently in a curvilinear "track". Detection of nuclear RNA was evidenced by hybridization without denaturation, sensitivity to RNAase, inhibition by actinomycin D, and specificity of transcribed sequences. Results indicate that RNA "tracks" extend from an internal genome into the nuclear periphery, and that RNA transport may be coupled to transcription. Localized nRNA is apparent for other viral sequences, different lymphoblastoid cell lines, nuclei prepared by two different methods, and an abundant, nonviral transfected sequence. Implications for understanding nuclear organization and the investigation of gene expression are discussed.

Determinants that contribute to cytoplasmic stability of human c-fos and beta-globin mRNAs are located at several sites in each mRNA

We have analyzed the contributions to cytoplasmic stability in an mRNA species with a very short half-life (human c-fos) and an mRNA species with a very long half-life (human beta-globin). When the human c-fos promoter was used to drive the expression of human c-fos, beta-globin, and chimeric DNAs between c-fos and beta-globin in transfected cells, a pulse of mRNA synthesis was obtained following induction of transcription by refeeding quiescent cells with medium containing 15% calf serum. The mRNA half-life was determined by using Northern (RNA) blot analysis of mRNAs prepared at various times following the pulse of transcription. Under these conditions human c-fos mRNA exhibited a half-life of 6.6 min and human beta-globin mRNA exhibited a half-life of 17.5 h. Replacement of the 3' end of the c-fos mRNA with the 3' end of the beta-globin mRNA increased the half-life of the resultant RNA from 6.6 to 34 min. The reciprocal chimera had a half-life of 34.6 min compared with the 17.5-h half-life of beta-globin mRNA. These results suggest that sequences which make a major contribution to mRNA stability reside in the 3' end of either or both molecules. A chimera in which the 5' untranslated region of globin was replaced by part of the 5' untranslated region of fos led to destabilization of the encoded mRNA. This construct produced an mRNA with a half-life of 6.8 h instead of the 17.5-h half-life of globin. This result suggests that additional determinants of stability reside in the 5' end of these mRNA molecules. Substitution of part of the 5' untranslated region of fos by the 5' untranslated region of beta-globin yielded an mRNA with stability similar to fos mRNA. These results suggest that interactions among sequences within each mRNA contribute to the stability of the respective molecules.

Multiple levels of regulation of tubulin gene expression during sea urchin embryogenesis

The expression of the tubulin genes during embryogenesis of the sea urchin Lytechinus pictus has been analyzed. Single strand tracer excess titrations of alpha- and beta-tubulin mRNA and RNA gel blot hybridizations indicate that tubulin mRNA remains at a constant 1.3 X 10(5) transcripts per embryo during cleavage stages, increases during ciliogenesis shortly before hatching (12 hr PF), declines until midgastrula (30-35 hr PF), and then gradually increases 3-fold to about 6 X 10(5) per pluteus larva (72 hr PF). Tubulin synthesis changes in concert with its mRNA, except that during cleavage the relative rate of tubulin synthesis increases without a corresponding increase in tubulin mRNA abundance. The relative rates of tubulin gene transcription were assayed by a run-on assay in isolated nuclei. The synthesis of alpha- or beta-tubulin RNA results in little supplementation of maternal tubulin RNA during cleavage stages, but the rate increases at least 18-fold during ciliogenesis and then gradually decreases thereafter. The accumulation of tubulin mRNA after gastrulation can be accounted for by an ontogenetic increase in tubulin RNA stability, assayed by actinomycin D chase and RNA gel blot hybridization. The rates of synthesis, stabilities, and abundances of alpha- and beta-tubulin mRNAs were similar, suggesting coordinate regulation. These observations indicate the importance of translational regulation during cleavage, transcriptional regulation during ciliogenesis, and regulation of mRNA stability by the level of unpolymerized tubulin during later development.

Determinants that contribute to cytoplasmic stability of human c-fos and beta-globin mRNAs are located at several sites in each mRNA

We have analyzed the contributions to cytoplasmic stability in an mRNA species with a very short half-life (human c-fos) and an mRNA species with a very long half-life (human beta-globin). When the human c-fos promoter was used to drive the expression of human c-fos, beta-globin, and chimeric DNAs between c-fos and beta-globin in transfected cells, a pulse of mRNA synthesis was obtained following induction of transcription by refeeding quiescent cells with medium containing 15% calf serum. The mRNA half-life was determined by using Northern (RNA) blot analysis of mRNAs prepared at various times following the pulse of transcription. Under these conditions human c-fos mRNA exhibited a half-life of 6.6 min and human beta-globin mRNA exhibited a half-life of 17.5 h. Replacement of the 3' end of the c-fos mRNA with the 3' end of the beta-globin mRNA increased the half-life of the resultant RNA from 6.6 to 34 min. The reciprocal chimera had a half-life of 34.6 min compared with the 17.5-h half-life of beta-globin mRNA. These results suggest that sequences which make a major contribution to mRNA stability reside in the 3' end of either or both molecules. A chimera in which the 5' untranslated region of globin was replaced by part of the 5' untranslated region of fos led to destabilization of the encoded mRNA. This construct produced an mRNA with a half-life of 6.8 h instead of the 17.5-h half-life of globin. This result suggests that additional determinants of stability reside in the 5' end of these mRNA molecules. Substitution of part of the 5' untranslated region of fos by the 5' untranslated region of beta-globin yielded an mRNA with stability similar to fos mRNA. These results suggest that interactions among sequences within each mRNA contribute to the stability of the respective molecules.

Determinants of mRNA stability in Dictyostelium discoideum amoebae: differences in poly(A) tail length, ribosome loading, and mRNA size cannot account for the heterogeneity of mRNA decay rates

As an approach to understanding the structures and mechanisms which determine mRNA decay rates, we have cloned and begun to characterize cDNAs which encode mRNAs representative of the stability extremes in the poly(A)+ RNA population of Dictyostelium discoideum amoebae. The cDNA clones were identified in a screening procedure which was based on the occurrence of poly(A) shortening during mRNA aging. mRNA half-lives were determined by hybridization of poly(A)+ RNA, isolated from cells labeled in a 32PO4 pulse-chase, to dots of excess cloned DNA. Individual mRNAs decayed with unique first-order decay rates ranging from 0.9 to 9.6 h, indicating that the complex decay kinetics of total poly(A)+ RNA in D. discoideum amoebae reflect the sum of the decay rates of individual mRNAs. Using specific probes derived from these cDNA clones, we have compared the sizes, extents of ribosome loading, and poly(A) tail lengths of stable, moderately stable, and unstable mRNAs. We found (i) no correlation between mRNA size and decay rate; (ii) no significant difference in the number of ribosomes per unit length of stable versus unstable mRNAs, and (iii) a general inverse relationship between mRNA decay rates and poly(A) tail lengths. Collectively, these observations indicate that mRNA decay in D. discoideum amoebae cannot be explained in terms of random nucleolytic events. The possibility that specific 3'-structural determinants can confer mRNA instability is suggested by a comparison of the labeling and turnover kinetics of different actin mRNAs. A correlation was observed between the steady-state percentage of a given mRNA found in polysomes and its degree of instability; i.e., unstable mRNAs were more efficiently recruited into polysomes than stable mRNAs. Since stable mRNAs are, on average, "older" than unstable mRNAs, this correlation may reflect a translational role for mRNA modifications that change in a time-dependent manner. Our previous studies have demonstrated both a time-dependent shortening and a possible translational role for the 3' poly(A) tracts of mRNA. We suggest, therefore, that the observed differences in the translational efficiency of stable and unstable mRNAs may, in part, be attributable to differences in steady-state poly(A) tail lengths.

Determinants of mRNA stability in Dictyostelium discoideum amoebae: differences in poly(A) tail length, ribosome loading, and mRNA size cannot account for the heterogeneity of mRNA decay rates

As an approach to understanding the structures and mechanisms which determine mRNA decay rates, we have cloned and begun to characterize cDNAs which encode mRNAs representative of the stability extremes in the poly(A)+ RNA population of Dictyostelium discoideum amoebae. The cDNA clones were identified in a screening procedure which was based on the occurrence of poly(A) shortening during mRNA aging. mRNA half-lives were determined by hybridization of poly(A)+ RNA, isolated from cells labeled in a 32PO4 pulse-chase, to dots of excess cloned DNA. Individual mRNAs decayed with unique first-order decay rates ranging from 0.9 to 9.6 h, indicating that the complex decay kinetics of total poly(A)+ RNA in D. discoideum amoebae reflect the sum of the decay rates of individual mRNAs. Using specific probes derived from these cDNA clones, we have compared the sizes, extents of ribosome loading, and poly(A) tail lengths of stable, moderately stable, and unstable mRNAs. We found (i) no correlation between mRNA size and decay rate; (ii) no significant difference in the number of ribosomes per unit length of stable versus unstable mRNAs, and (iii) a general inverse relationship between mRNA decay rates and poly(A) tail lengths. Collectively, these observations indicate that mRNA decay in D. discoideum amoebae cannot be explained in terms of random nucleolytic events. The possibility that specific 3'-structural determinants can confer mRNA instability is suggested by a comparison of the labeling and turnover kinetics of different actin mRNAs. A correlation was observed between the steady-state percentage of a given mRNA found in polysomes and its degree of instability; i.e., unstable mRNAs were more efficiently recruited into polysomes than stable mRNAs. Since stable mRNAs are, on average, "older" than unstable mRNAs, this correlation may reflect a translational role for mRNA modifications that change in a time-dependent manner. Our previous studies have demonstrated both a time-dependent shortening and a possible translational role for the 3' poly(A) tracts of mRNA. We suggest, therefore, that the observed differences in the translational efficiency of stable and unstable mRNAs may, in part, be attributable to differences in steady-state poly(A) tail lengths.