Transcriptional buffering and 3'UTR lengthening are shaped during human neurodevelopment by shifts in mRNA stability and microRNA load

The contribution of mRNA half-life is commonly overlooked when examining changes in mRNA abundance during development. mRNA levels of some genes are regulated by transcription rate only, but others may be regulated by mRNA half-life only shifts. Furthermore, transcriptional buffering is predicted when changes in transcription rates have compensating shifts in mRNA half-life resulting in no change to steady-state levels. Likewise, transcriptional boosting should result when changes in transcription rate are accompanied by amplifying half-life shifts. During neurodevelopment there is widespread 3'UTR lengthening that could be shaped by differential shifts in the stability of existing short or long 3'UTR transcript isoforms. We measured transcription rate and mRNA half-life changes during induced human Pluripotent Stem Cell (iPSC)-derived neuronal development using RATE-seq. During transitions to progenitor and neuron stages, transcriptional buffering occurred in up to 50%, and transcriptional boosting in up to 15%, of genes with changed transcription rates. The remaining changes occurred by transcription rate only or mRNA half-life only shifts. Average mRNA half-life decreased two-fold in neurons relative to iPSCs. Short gene isoforms were more destabilized in neurons and thereby increased the average 3'UTR length. Small RNA sequencing captured an increase in microRNA copy number per cell during neurodevelopment. We propose that mRNA destabilization and 3'UTR lengthening are driven in part by an increase in microRNA load in neurons. Our findings identify mRNA stability mechanisms in human neurodevelopment that regulate gene and isoform level abundance and provide a precedent for similar post-transcriptional regulatory events as other tissues develop.

RNAPII driven post-translational modifications of nucleosomal histones

The current understanding of how specific distributions of histone post-translational modifications (PTMs) are achieved throughout the chromatin remains incomplete. This review focuses on the role of RNA polymerase II (RNAPII) in establishing H2BK120/K123 ubiquitination and H3K4/K36 methylation distribution. The rate of RNAPII transcription is mainly a function of the RNAPII elongation and recruitment rates. Two major mechanisms link RNAPII's transcription rate to the distribution of PTMs. First, the phosphorylation patterns of Ser2P/Ser5P in the C-terminal domain of RNAPII change as a function of time, since the start of elongation, linking them to the elongation rate. Ser2P/Ser5P recruits specific histone PTM enzymes/activators to the nucleosome. Second, multiple rounds of binding and catalysis by the enzymes are required to establish higher methylations (H3K4/36me3). Thus, methylation states are determined by the transcription rate. In summary, the first mechanism determines the location of methylations in the gene, while the second mechanism determines the methylation state.

A Trans-Omics Comparison Reveals Common Gene Expression Strategies in Four Model Organisms and Exposes Similarities and Differences between Them

The ultimate goal of gene expression regulation is on the protein level. However, because the amounts of mRNAs and proteins are controlled by their synthesis and degradation rates, the cellular amount of a given protein can be attained by following different strategies. By studying omics data for six expression variables (mRNA and protein amounts, plus their synthesis and decay rates), we previously demonstrated the existence of common expression strategies (CESs) for functionally related genes in the yeast Saccharomyces cerevisiae. Here we extend that study to two other eukaryotes: the yeast Schizosaccharomyces pombe and cultured human HeLa cells. We also use genomic data from the model prokaryote Escherichia coli as an external reference. We show that six-variable profiles (6VPs) can be constructed for every gene and that these 6VPs are similar for genes with similar functions in all the studied organisms. The differences in 6VPs between organisms can be used to establish their phylogenetic relationships. The analysis of the correlations among the six variables supports the hypothesis that most gene expression control occurs in actively growing organisms at the transcription rate level, and that translation plays a minor role. We propose that living organisms use CESs for the genes acting on the same physiological pathways, especially for those belonging to stable macromolecular complexes, but CESs have been modeled by evolution to adapt to the specific life circumstances of each organism.

A Protocol for Transcriptome-Wide Inference of RNA Metabolic Rates in Mouse Embryonic Stem Cells

The relative ease of mouse Embryonic Stem Cells (mESCs) culture and the potential of these cells to differentiate into any of the three primary germ layers: ectoderm, endoderm and mesoderm (pluripotency), makes them an ideal and frequently used ex vivo system to dissect how gene expression changes impact cell state and differentiation. These efforts are further supported by the large number of constitutive and inducible mESC mutants established with the aim of assessing the contributions of different pathways and genes to cell homeostasis and gene regulation. Gene product abundance is controlled by the modulation of the rates of RNA synthesis, processing, and degradation. The ability to determine the relative contribution of these different RNA metabolic rates to gene expression control using standard RNA-sequencing approaches, which only capture steady state abundance of transcripts, is limited. In contrast, metabolic labeling of RNA with 4-thiouridine (4sU) coupled with RNA-sequencing, allows simultaneous and reproducible inference of transcriptome wide synthesis, processing, and degradation rates. Here we describe, a detailed protocol for 4sU metabolic labeling in mESCs that requires short 4sU labeling times at low concentration and minimally impacts cellular homeostasis. This approach presents a versatile method for in-depth characterization of the gene regulatory strategies governing gene steady state abundance in mESC.

Mitigating Antagonism between Transcription and Proliferation Allows Near-Deterministic Cellular Reprogramming

Although cellular reprogramming enables the generation of new cell types for disease modeling and regenerative therapies, reprogramming remains a rare cellular event. By examining reprogramming of fibroblasts into motor neurons and multiple other somatic lineages, we find that epigenetic barriers to conversion can be overcome by endowing cells with the ability to mitigate an inherent antagonism between transcription and DNA replication. We show that transcription factor overexpression induces unusually high rates of transcription and that sustaining hypertranscription and transgene expression in hyperproliferative cells early in reprogramming is critical for successful lineage conversion. However, hypertranscription impedes DNA replication and cell proliferation, processes that facilitate reprogramming. We identify a chemical and genetic cocktail that dramatically increases the number of cells capable of simultaneous hypertranscription and hyperproliferation by activating topoisomerases. Further, we show that hypertranscribing, hyperproliferating cells reprogram at 100-fold higher, near-deterministic rates. Therefore, relaxing biophysical constraints overcomes molecular barriers to cellular reprogramming.

Single molecule approaches for studying gene regulation in metabolic tissues

Gene expression in metabolic tissues can be regulated at multiple levels, ranging from the control of promoter accessibilities, transcription rates, mRNA degradation rates and mRNA localization. Modulating these processes can differentially affect important performance criteria of cells. These include precision, cellular economy, rapid response and maintenance of DNA integrity. In this review we will describe how distinct strategies of gene regulation impact the trade-offs between the cells' performance criteria. We will highlight tools based on single molecule visualization of transcripts that can be used to measure promoter states, transcription rates and mRNA degradation rates in intact tissues. These approaches revealed surprising recurrent patterns in mammalian tissues, that include transcriptional bursting, nuclear retention of mRNA, and coordination of mRNA lifetimes to facilitate rapid adaptation to changing metabolic inputs. The ability to characterize gene expression at the single molecule level can uncover the design principles of gene regulation in metabolic tissues such as the liver and the pancreas.

Response to Nodal morphogen gradient is determined by the kinetics of target gene induction

Morphogen gradients expose cells to different signal concentrations and induce target genes with different ranges of expression. To determine how the Nodal morphogen gradient induces distinct gene expression patterns during zebrafish embryogenesis, we measured the activation dynamics of the signal transducer Smad2 and the expression kinetics of long- and short-range target genes. We found that threshold models based on ligand concentration are insufficient to predict the response of target genes. Instead, morphogen interpretation is shaped by the kinetics of target gene induction: the higher the rate of transcription and the earlier the onset of induction, the greater the spatial range of expression. Thus, the timing and magnitude of target gene expression can be used to modulate the range of expression and diversify the response to morphogen gradients.

DOI:http://dx.doi.org/10.7554/eLife.05042.001

How a cell can tell where it is in a developing embryo has fascinated scientists for decades. The pioneering computer scientist and mathematical biologist Alan Turing was the first person to coin the term ‘morphogen’ to describe a protein that provides information about locations in the body. A morphogen is released from a group of cells (called the ‘source’) and as it moves away its activity (called the ‘signal’) declines gradually. Cells sense this signal gradient and use it to detect their position with respect to the source. Nodal is an important morphogen and is required to establish the correct identity of cells in the embryo; for example, it helps determine which cells should become a brain or heart or gut cell and so on.

The zebrafish is a widely used model to study animal development, in part because its embryos are transparent; this allows cells and proteins to be easily observed under a microscope. When Nodal acts on cells, another protein called Smad2 becomes activated, moves into the cell's nucleus, and then binds to specific genes. This triggers the expression of these genes, which are first copied into mRNA molecules via a process known as transcription and are then translated into proteins. The protein products of these targeted genes control cell identity and movement.

Several models have been proposed to explain how different concentrations of Nodal switch on the expression of different target genes; that is to say, to explain how a cell interprets the Nodal gradient. Dubrulle et al. have now measured factors that underlie how this gradient is interpreted. Individual cells in zebrafish embryos were tracked under a microscope, and Smad2 activation and gene expression were assessed. Dubrulle et al. found that, in contradiction to previous models, the amount of Nodal present on its own was insufficient to predict the target gene response. Instead, their analysis suggests that the size of each target gene's response depends on its rate of transcription and how quickly it is first expressed in response to Nodal.

These findings of Dubrulle et al. suggest that timing and transcription rate are important in determining the appropriate response to Nodal. Further work will be now needed to find out whether similar mechanisms regulate other processes that rely on the activity of morphogens.

DOI:http://dx.doi.org/10.7554/eLife.05042.002

Cytoplasmic 5'-3' exonuclease Xrn1p is also a genome-wide transcription factor in yeast

The 5′ to 3′ exoribonuclease Xrn1 is a large protein involved in cytoplasmatic mRNA degradation as a critical component of the major decaysome. Its deletion in the yeast Saccharomyces cerevisiae is not lethal, but it has multiple physiological effects. In a previous study, our group showed that deletion of all tested components of the yeast major decaysome, including XRN1, results in a decrease in the synthetic rate and an increase in half-life of most mRNAs in a compensatory manner. Furthermore, the same study showed that the all tested decaysome components are also nuclear proteins that bind to the 5′ region of a number of genes. In the present work, we show that disruption of Xrn1 activity preferentially affects both the synthesis and decay of a distinct subpopulation of mRNAs. The most affected mRNAs are the transcripts of the highly transcribed genes, mainly those encoding ribosome biogenesis and translation factors. Previously, we proposed that synthegradases play a key role in regulating both mRNA synthesis and degradation. Evidently, Xrn1 functions as a synthegradase, whose selectivity might help coordinating the expression of the protein synthetic machinery. We propose to name the most affected genes “Xrn1 synthegradon.”

External conditions inversely change the RNA polymerase II elongation rate and density in yeast

Elongation speed is a key parameter in RNA polymerase II (RNA pol II) activity. It affects the transcription rate, while it is conditioned by the physicochemical environment it works in at the same time. For instance, it is well-known that temperature affects the biochemical reactions rates. Therefore in free-living organisms that are able to grow at various environmental temperatures, such as the yeast Saccharomyces cerevisiae, evolution should have not only shaped the structural and functional properties of this key enzyme, but should have also provided mechanisms and pathways to adapt its activity to the optimal performance required. We studied the changes in RNA pol II elongation speed caused by alternations in growth temperature in yeast to find that they strictly follow the Arrhenius equation, and that they also provoke an almost inverse proportional change in RNA pol II density within the optimal growth temperature range (26-37 °C). Moreover, we discovered that yeast cells control the transcription initiation rate by changing the total amount of available RNA pol II.

What do you mean by transcription rate?: the conceptual difference between nascent transcription rate and mRNA synthesis rate is essential for the proper understanding of transcriptomic analyses

mRNA synthesis in all organisms is performed by RNA polymerases, which work as nanomachines on DNA templates. The rate at which their product is made is an important parameter in gene expression. Transcription rate encompasses two related, yet different, concepts: the nascent transcription rate, which measures the in situ mRNA production by RNA polymerase, and the rate of synthesis of mature mRNA, which measures the contribution of transcription to the mRNA concentration. Both parameters are useful for molecular biologists, but they are not interchangeable and they are expressed in different units. It is important to distinguish when and where each one should be used. We propose that for functional genomics the use of nascent transcription rates should be restricted to the evaluation of the transcriptional process itself, whereas mature mRNA synthesis rates should be employed to address the transcriptional input to mRNA concentration balance leading to variation of gene expression.

Eukaryotic mRNA decay: methodologies, pathways, and links to other stages of gene expression

mRNA concentration depends on the balance between transcription and degradation rates. On both sides of the equilibrium, synthesis and degradation show, however, interesting differences that have conditioned the evolution of gene regulatory mechanisms. Here, we discuss recent genome-wide methods for determining mRNA half-lives in eukaryotes. We also review pre- and posttranscriptional regulons that coordinate the fate of functionally related mRNAs by using protein- or RNA-based trans factors. Some of these factors can regulate both transcription and decay rates, thereby maintaining proper mRNA homeostasis during eukaryotic cell life.

Do DNA microarrays tell the story of gene expression?

Poor reproducibility of microarray measurements is a major obstacle to their application as an instrument for clinical diagnostics. In this paper, several aspects of poor reproducibility are analyzed. All of them belong to the category of interpretive weaknesses of DNA microarray technology. First, the attention is drawn to the fact that absence of the information regarding post-transcriptional mRNA stability makes it impossible to evaluate the level of gene activity from the relative mRNA abundances, the quantities available from microarray measurements. Second, irreducible intracellular variability with persistent patterns of stochasticity and burstiness put natural limits to reproducibility. Third, strong interactions within intracellular biomolecular networks make it highly problematic to build a bridge between transcription rates of individual genes and structural fidelity of their genetic codes. For these reasons, the microarray measurements of relative mRNA abundances are more appropriate in laboratory settings as a tool for scientific research, hypotheses generating and producing the leads for subsequent validation through more sophisticated technologies. As to clinical settings, where firm conclusive diagnoses, not the leads for further experimentation, are required, microarrays still have a long way to go until they become a reliable instrument in patient-related decision making.

Elevated FMR1 mRNA in premutation carriers is due to increased transcription

Carriers of premutation alleles (55-200 CGG repeats) of the fragile X mental retardation 1 (FMR1) gene have levels of FMR1 mRNA that are elevated by as much as 10-fold in peripheral blood leukocytes and CNS tissue. The excess expanded-repeat mRNA, per se, is now believed to result in forms of clinical involvement that are largely restricted to premutation carriers, including the neurodegenerative disorder, fragile X-associated tremor/ataxia syndrome (FXTAS). Although evidence to date suggests that the elevated mRNA is not due to increased stability, the basis for the increase is not known. In the current study, we have determined the relative transcriptional activities of premutation and normal FMR1 alleles using a highly sensitive nuclear run-on assay that involves immunocapture of digoxigenin-modified run-on transcripts followed by PCR amplification of the nascent transcripts. Using the nuclear run-on approach, we demonstrate that the rate of run-on synthesis of FMR1 transcripts is increased in premutation alleles. The current run-on assay should be broadly applicable to studies of other genes with promoters of weak to moderate strength. The fraction of capped FMR1 mRNA remains unaltered for premutation transcripts, indicating that elevated message levels are not due to premature escape from the cotranscriptional capping process. We also show that, in contrast to the situation with myotonic dystrophy, there is no net nuclear sequestration of premutation FMR1 mRNA. Finally, we have demonstrated that AGG interruptions within the CGG repeat element do not influence FMR1 mRNA levels.