A network of multiple regulatory layers shapes gene expression in fission yeast

The impact of multiple abiotic stresses on ns-LTP2.8 gene transcript and ns-LTP2.8 protein accumulation in germinating barley (Hordeum vulgare L.) embryos

Abiotic stresses occur more often in combination than alone under regular field conditions limiting in more severe way crop production. Stress recognition in plants primarily occurs in the plasma membrane, modification of which is necessary to maintain homeostasis in response to it. It is known that lipid transport proteins (ns-LTPs) participate in modification of the lipidome of cell membranes. Representative of this group, ns-LTP2.8, may be involved in the reaction to abiotic stress of germinating barley plants by mediating the intracellular transport of hydrophobic particles, such as lipids, helping to maintain homeostasis. The ns-LTP2.8 protein was selected for analysis due to its ability to transport not only linear hydrophobic molecules but also compounds with a more complex spatial structure. Moreover, ns-LTP2.8 has been qualified as a member of pathogenesis-related proteins, which makes it particularly important in relation to its high allergenic potential. This paper demonstrates for the first time the influence of various abiotic stresses acting separately as well as in their combinations on the change in the ns-LTP2.8 transcript, ns-LTP2.8 protein and total soluble protein content in the embryonal axes of germinating spring barley genotypes with different ns-LTP2.8 allelic forms and stress tolerance. Tissue localization of ns-LTP2.8 transcript as well as ns-LTP2.8 protein were also examined. Although the impact of abiotic stresses on the regulation of gene transcription and translation processes remains not fully recognized, in this work we managed to demonstrate different impact on applied stresses on the fundamental cellular processes in very little studied tissue of the embryonal axis of barley.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Broad functional profiling of fission yeast proteins using phenomics and machine learning

Many proteins remain poorly characterized even in well-studied organisms, presenting a bottleneck for research. We applied phenomics and machine-learning approaches with Schizosaccharomyces pombe for broad cues on protein functions. We assayed colony-growth phenotypes to measure the fitness of deletion mutants for 3509 non-essential genes in 131 conditions with different nutrients, drugs, and stresses. These analyses exposed phenotypes for 3492 mutants, including 124 mutants of 'priority unstudied' proteins conserved in humans, providing varied functional clues. For example, over 900 proteins were newly implicated in the resistance to oxidative stress. Phenotype-correlation networks suggested roles for poorly characterized proteins through 'guilt by association' with known proteins. For complementary functional insights, we predicted Gene Ontology (GO) terms using machine learning methods exploiting protein-network and protein-homology data (NET-FF). We obtained 56,594 high-scoring GO predictions, of which 22,060 also featured high information content. Our phenotype-correlation data and NET-FF predictions showed a strong concordance with existing PomBase GO annotations and protein networks, with integrated analyses revealing 1675 novel GO predictions for 783 genes, including 47 predictions for 23 priority unstudied proteins. Experimental validation identified new proteins involved in cellular aging, showing that these predictions and phenomics data provide a rich resource to uncover new protein functions.

Genomic alterations involved in fluoroquinolone resistance development in Staphylococcus aureus

AIM

Fluoroquinolone (FQ) is a potent antibiotic class. However, resistance to this class emerges quickly which hinders its application. In this study, mechanisms leading to the emergence of multidrug-resistant (MDR) Staphylococcus aureus (S. aureus) strains under FQ exposure were investigated.

METHODOLOGY

S. aureus ATCC 29213 was serially exposed to ciprofloxacin (CIP), ofloxacin (OFL), or levofloxacin (LEV) at sub-minimum inhibitory concentrations (sub-MICs) for 12 days to obtain S. aureus -1 strains and antibiotic-free cultured for another 10 days to obtain S. aureus-2 strains. The whole genome (WGS) and target sequencing were applied to analyze genomic alterations; and RT-qPCR was used to access the expressions of efflux-related genes, alternative sigma factors, and genes involved in FQ resistance.

RESULTS

A strong and irreversible increase of MICs was observed in all applied FQs (32 to 128 times) in all S. aureus-1 and remained 16 to 32 times in all S. aureus-2. WGS indicated 10 noticeable mutations occurring in all FQ-exposed S. aureus including 2 insdel mutations in SACOL0573 and rimI; a synonymous mutation in hslO; and 7 missense mutations located in an untranslated region. GrlA, was found mutated (R570H) in all S. aureus-1 and -2. Genes encoding for efflux pumps and their regulator (norA, norB, norC, and mgrA); alternative sigma factors (sigB and sigS); acetyltransferase (rimI); methicillin resistance (fmtB); and hypothetical protein BJI72_0645 were overexpressed in FQ-exposed strains.

CONCLUSION

The emergence of MDR S. aureus was associated with the mutations in the FQ-target sequences and the overexpression of efflux pump systems and their regulators.

Translation in the cell under fierce competition for shared resources: a mathematical model

During translation, mRNAs 'compete' for shared resources. Under stress conditions, during viral infection and also in high-throughput heterologous gene expression, these resources may become scarce, e.g. the pool of free ribosomes is starved, and then the competition may have a dramatic effect on the global dynamics of translation in the cell. We model this scenario using a network that includes m ribosome flow models (RFMs) interconnected via a pool of free ribosomes. Each RFM models ribosome flow along an mRNA molecule, and the pool models the shared resource. We assume that the number of mRNAs is large, so many ribosomes are attached to the mRNAs, and the pool is starved. Our analysis shows that adding an mRNA has an intricate effect on the total protein production. The new mRNA produces new proteins, but the other mRNAs produce less proteins, as the pool that feeds these mRNAs now has a smaller abundance of ribosomes. As the number of mRNAs increases, the marginal utility of adding another mRNA diminishes, and the total protein production rate saturates to a limiting value. We demonstrate our approach using an example of insulin protein production in a cell-free system.

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high-throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing Translation > Regulation RNA Evolution and Genomics > Computational Analyses of RNA.

Three-layered control of mRNA poly(A) tail synthesis in Saccharomyces cerevisiae

Biogenesis of most eukaryotic mRNAs involves the addition of an untemplated polyadenosine (pA) tail by the cleavage and polyadenylation machinery. The pA tail, and its exact length, impacts mRNA stability, nuclear export, and translation. To define how polyadenylation is controlled in S. cerevisiae, we have used an in vivo assay capable of assessing nuclear pA tail synthesis, analyzed tail length distributions by direct RNA sequencing, and reconstituted polyadenylation reactions with purified components. This revealed three control mechanisms for pA tail length. First, we found that the pA binding protein (PABP) Nab2p is the primary regulator of pA tail length. Second, when Nab2p is limiting, the nuclear pool of Pab1p, the second major PABP in yeast, controls the process. Third, when both PABPs are absent, the cleavage and polyadenylation factor (CPF) limits pA tail synthesis. Thus, Pab1p and CPF provide fail-safe mechanisms to a primary Nab2p-dependent pathway, thereby preventing uncontrolled polyadenylation and allowing mRNA export and translation.

Factors affecting the rapid changes of protein under short-term heat stress

BACKGROUND

Protein content determines the state of cells. The variation in protein abundance is crucial when organisms are in the early stages of heat stress, but the reasons affecting their changes are largely unknown.

RESULTS

We quantified 47,535 mRNAs and 3742 proteins in the filling grains of wheat in two different thermal environments. The impact of mRNA abundance and sequence features involved in protein translation and degradation on protein expression was evaluated by regression analysis. Transcription, codon usage and amino acid frequency were the main drivers of changes in protein expression under heat stress, and their combined contribution explains 58.2 and 66.4% of the protein variation at 30 and 40 °C (20 °C as control), respectively. Transcription contributes more to alterations in protein content at 40 °C (31%) than at 30 °C (6%). Furthermore, the usage of codon AAG may be closely related to the rapid alteration of proteins under heat stress. The contributions of AAG were 24 and 13% at 30 and 40 °C, respectively.

CONCLUSION

In this study, we analyzed the factors affecting the changes in protein expression in the early stage of heat stress and evaluated their influence.

Stochastic microswimming model of ribosome motion on the polysome

In this work we assume that the ribosome propels itself during the translocation step of the translation process of protein synthesis by running a cycle of stochastically generated conformational changes involving its two subunits. This cycle includes only two experimentally found ribosome shape changes. The main result is an analytic expression for ribosome's average swimming speed on a polysome, where the ribosome is in the presence of other ribosomes. Relevant geometric parameters of ribosome deformations are calculated first by solving a deterministic problem where the ribosome runs a cycle of prescribed conformational changes. The method of reflections and pairwise additivity are used to obtain the stresses and forces needed to apply the multiparticle reciprocal theorem. Ribosome's average velocity when it runs the corresponding stochastic cycle of deformations is calculated assuming independence among the conformational cycles of different ribosomes on the polysome. The results obtained show that swimming in tandem on the polysome allows the ribosome to reach any typical subcellular speed with deformations whose amplitude is of a smaller size than when it swims alone in the fluid. Also, the flow organized by its swimming stroke becomes more determinant for its motion than random diffusion, compared to the solitary ribosome.

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.

Ribosome profiling reveals ribosome stalling on tryptophan codons and ribosome queuing upon oxidative stress in fission yeast

Translational control is essential in response to stress. We investigated the translational programmes launched by the fission yeast Schizosaccharomyces pombe upon five environmental stresses. We also explored the contribution of defence pathways to these programmes: The Integrated Stress Response (ISR), which regulates translation initiation, and the stress-response MAPK pathway. We performed ribosome profiling of cells subjected to each stress, in wild type cells and in cells with the defence pathways inactivated. The transcription factor Fil1, a functional homologue of the yeast Gcn4 and the mammalian Atf4 proteins, was translationally upregulated and required for the response to most stresses. Moreover, many mRNAs encoding proteins required for ribosome biogenesis were translationally downregulated. Thus, several stresses trigger a universal translational response, including reduced ribosome production and a Fil1-mediated transcriptional programme. Surprisingly, ribosomes stalled on tryptophan codons upon oxidative stress, likely due to a decrease in charged tRNA-Tryptophan. Stalling caused ribosome accumulation upstream of tryptophan codons (ribosome queuing/collisions), demonstrating that stalled ribosomes affect translation elongation by other ribosomes. Consistently, tryptophan codon stalling led to reduced translation elongation and contributed to the ISR-mediated inhibition of initiation. We show that different stresses elicit common and specific translational responses, revealing a novel role in Tryptophan-tRNA availability.

eIF3k Domain-Containing Protein Regulates Conidiogenesis, Appressorium Turgor, Virulence, Stress Tolerance, and Physiological and Pathogenic Development of Magnaporthe oryzae Oryzae

The eukaryotic translation initiation factor 3 (eIF3) complex consists of essential and non-essential sub-complexes. Non-essential eIF3 complex subunits, such as eIF3e, eIF3j, eIF3k, and eIF3l, modulate stress tolerance and enhance the lifespan of Neurospora crassa and Caenorhabditis elegans. However, there is limited knowledge of the role of the non-essential eIF3 sub-complex in the pathophysiological development of plant fungal pathogens. Here, we deployed genetic and biochemical techniques to explore the influence of a hypothetical protein containing eIF3k domain in Magnaporthe oryzae Oryzae (MoOeIF3k) on reproduction, hyphae morphogenesis, stress tolerance, and pathogenesis. Also, the targeted disruption of MoOeIF3k suppressed vegetative growth and asexual sporulation in ΔMoOeif3k strains significantly. We demonstrated that MoOeIF3k promotes the initiation and development of the rice blast disease by positively regulating the mobilization and degradation of glycogen, appressorium integrity, host penetration, and colonization during host-pathogen interaction. For the first time, we demonstrated that the eIF3k subunit supports the survival of the blast fungus by suppressing vegetative growth and possibly regulating the conversions and utilization of stored cellular energy reserves under starvation conditions. We also observed that the deletion of MoOeIF3k accelerated ribosomal RNA (rRNA) generation in the ΔMoOeif3k strains with a corresponding increase in total protein output. In summary, this study unravels the pathophysiological significance of eIF3k filamentous fungi. The findings also underscored the need to systematically evaluate the individual subunits of the non-essential eIF3 sub-complex during host-pathogen interaction. Further studies are required to unravel the influence of synergetic coordination between translation and transcriptional regulatory machinery on the pathogenesis of filamentous fungi pathogens.

Gcn2 eIF2α kinase mediates combinatorial translational regulation through nucleotide motifs and uORFs in target mRNAs

The protein kinase Gcn2 is a central transducer of nutritional stress signaling important for stress adaptation by normal cells and the survival of cancer cells. In response to nutrient deprivation, Gcn2 phosphorylates eIF2α, thereby repressing general translation while enhancing translation of specific mRNAs with upstream ORFs (uORFs) situated in their 5'-leader regions. Here we performed genome-wide measurements of mRNA translation during histidine starvation in fission yeast Schizosaccharomyces pombe. Polysome analyses were combined with microarray measurements to identify gene transcripts whose translation was up-regulated in response to the stress in a Gcn2-dependent manner. We determined that translation is reprogrammed to enhance RNA metabolism and chromatin regulation and repress ribosome synthesis. Interestingly, translation of intron-containing mRNAs was up-regulated. The products of the regulated genes include additional eIF2α kinase Hri2 amplifying the stress signaling and Gcn5 histone acetyl transferase and transcription factors, together altering genome-wide transcription. Unique dipeptide-coding uORFs and nucleotide motifs, such as '5'-UGA(C/G)GG-3', are found in 5' leader regions of regulated genes and shown to be responsible for translational control.

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

PAT-Seq: A Method for Simultaneous Quantitation of Gene Expression, Poly(A)-Site Selection and Poly(A)-Length Distribution in Yeast Transcriptomes

Next-generation sequencing (NGS) and its application to RNA (RNA-seq) has opened up multiple aspects of RNA processing to deep transcriptome-wide analysis at nucleotide resolution. This has been useful in delineating the transcribed areas of the genome, and in quantitation of RNA isoforms. Such isoforms can diversify the regulatory repertoire of mRNAs. For example, the 3'-end of mRNA can vary in two important ways, in the position chosen for cleavage and polyadenylation, and in the length of the poly(A)-tail. Accordingly, the step-up in resolution made possible by NGS has revealed an unexpectedly high level of alternative polyadenylation (APA). Moreover, it has massively simplified the transcriptome-wide detection of poly(A)-tail length changes. Here we present our approach to the study of 3'-end dynamics using a 3'-focused RNA-seq method called PAT-seq (for poly(A)-test sequencing). The approach records gene expression, APA, and poly(A)-tail changes between transcriptomes to reveal complex interplay between transcriptional and posttranscriptional control mechanisms.

Establishment of 5'-3' interactions in mRNA independent of a continuous ribose-phosphate backbone

Functions of eukaryotic mRNAs are characterized by intramolecular interactions between their ends. We have addressed the question whether 5' and 3' ends meet by diffusion-controlled encounter "through solution" or by a mechanism involving the RNA backbone. For this purpose, we used a translation system derived from Drosophila embryos that displays two types of 5'-3' interactions: Cap-dependent translation initiation is stimulated by the poly(A) tail and inhibited by Smaug recognition elements (SREs) in the 3' UTR. Chimeric RNAs were made consisting of one RNA molecule carrying a luciferase coding sequence and a second molecule containing SREs and a poly(A) tail; the two were connected via a protein linker. The poly(A) tail stimulated translation of such chimeras even when disruption of the RNA backbone was combined with an inversion of the 5'-3' polarity between the open reading frame and poly(A) segment. Stimulation by the poly(A) tail also decreased with increasing RNA length. Both observations suggest that contacts between the poly(A) tail and the 5' end are established through solution, independently of the RNA backbone. In the same chimeric constructs, SRE-dependent inhibition of translation was also insensitive to disruption of the RNA backbone. Thus, tracking of the backbone is not involved in the repression of cap-dependent initiation. However, SRE-dependent repression was insensitive to mRNA length, suggesting that the contact between the SREs in the 3' UTR and the 5' end of the RNA might be established in a manner that differs from the contact between the poly(A) tail and the cap.

Systematic Comparisons of Positively Selected Genes between Gossypium arboreum and Gossypium raimondii Genomes

Background:

Studies of Positively Selected Genes (PSGs) in microorganisms and mammals have provided insights into the dynamics of genome evolution and the genetic basis of differences between species by using whole genome-wide scans. Systematic investigations and comparisons of PSGs in plants, however, are still limited.

Objective:

A systematic comparison of PSGs between the genomes of two cotton species, Gossypium arboreum (G. arboreum) and G. raimondii, will give the key answer for revealing molecular evolutionary differences in plants.

Methods:

Genome sequences of G. arboreum and G. raimondii were compared, including Whole Genome Duplication (WGD) events and genomic features such as gene number, gene length, codon bias index, evolutionary rate, number of expressed genes, and retention of duplicated copies.

Results:

Unlike the PSGs in G. raimondii, G. arboreum comprised more PSGs, smaller gene size and fewer expressed gene. In addition, the PSGs evolved at a higher rate of synonymous substitutions, but were subjected to lower selection pressure. The PSGs in G. arboreum were also retained with a lower number of duplicate gene copies than G. raimondii after a single WGD event involving Gossypium.

Conclusion:

These data indicate that PSGs in G. arboreum and G. raimondii differ not only in Ka/Ks, but also in their evolutionary, structural, and expression properties, indicating that divergence of G. arboreum and G. raimondii was associated with differences in PSGs in terms of evolutionary rates, gene length, expression patterns, and WGD retention in Gossypium.

CFIm25 and alternative polyadenylation: Conflicting roles in cancer

Alternative polyadenylation (APA) is now widely recognized to regulate gene expression. APA is an RNA-processing mechanism that generates distinct 3' termini on mRNAs, producing mRNA isoforms. Different factors influence the initiation and development of this process. CFIm25 (among others) is a cleavage and polyadenylation factor that plays a key role in the regulation of APA. Shortening of the 3'UTRs on mRNAs leads to enhanced cellular proliferation and tumorigenicity. One reason may be the up-regulation of growth promoting factors, such as Cyclin D1. Different studies have reported a dual role of CFIm25 in cancer (both oncogenic and tumor suppressor). microRNAs (miRNAs) may be involved in CFIm25 function as well as competing endogenous RNAs (ceRNAs). The present review focuses on the role of CFIm25 in cancer, cancer treatment, and possible involvement in other human diseases. We highlight the involvement of miRNAs and ceRNAs in the function of CFIm25 to affect gene expression. The lack of understanding of the mechanisms and regulation of CFIm25 and APA has underscored the need for further research regarding their role in cancer and other diseases.

Does functional specialization of ribosomes really exist?

It has recently become clear that ribosomes are much more heterogeneous than previously thought, with diversity arising from rRNA sequence and modifications, ribosomal protein (RP) content and posttranslational modifications (PTMs), as well as bound nonribosomal proteins. In some cases, the existence of these diverse ribosome populations has been verified by biochemical or structural methods. Furthermore, knockout or knockdown of RPs can diversify ribosome populations, while also affecting the translation of some mRNAs (but not others) with biological consequences. However, the effects on translation arising from depletion of diverse proteins can be highly similar, suggesting that there may be a more general defect in ribosome function or stability, perhaps arising from reduced ribosome numbers. Consistently, overall reduced ribosome numbers can differentially affect subclasses of mRNAs, necessitating controls for specificity. Moreover, in order to study the functional consequences of ribosome diversity, perturbations including affinity tags and knockouts are introduced, which can also affect the outcome of the experiment. Here we review the available literature to carefully evaluate whether the published data support functional diversification, defined as diverse ribosome populations differentially affecting translation of distinct mRNA (classes). Based on these observations and the commonly observed cellular responses to perturbations in the system, we suggest a set of important controls to validate functional diversity, which should include gain-of-function assays and the demonstration of inducibility under physiological conditions.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

Layers of regulation of cell-cycle gene expression in the budding yeast Saccharomyces cerevisiae

In the budding yeast Saccharomyces cerevisiae, transcription factors (TFs) regulate the periodic expression of many genes during the cell cycle, including gene products required for progression through cell-cycle events. Experimental evidence coupled with quantitative models suggests that a network of interconnected TFs is capable of regulating periodic genes over the cell cycle. Importantly, these dynamical models were built on transcriptomics data and assumed that TF protein levels and activity are directly correlated with mRNA abundance. To ask whether TF transcripts match protein expression levels as cells progress through the cell cycle, we applied a multiplexed targeted mass spectrometry approach (parallel reaction monitoring) to synchronized populations of cells. We found that protein expression of many TFs and cell-cycle regulators closely followed their respective mRNA transcript dynamics in cycling wild-type cells. Discordant mRNA/protein expression dynamics was also observed for a subset of cell-cycle TFs and for proteins targeted for degradation by E3 ubiquitin ligase complexes such as SCF (Skp1/Cul1/F-box) and APC/C (anaphase-promoting complex/cyclosome). We further profiled mutant cells lacking B-type cyclin/CDK activity ( clb1-6) where oscillations in ubiquitin ligase activity, cyclin/CDKs, and cell-cycle progression are halted. We found that a number of proteins were no longer periodically degraded in clb1-6 mutants compared with wild type, highlighting the importance of posttranscriptional regulation. Finally, the TF complexes responsible for activating G1/S transcription (SBF and MBF) were more constitutively expressed at the protein level than at periodic mRNA expression levels in both wild-type and mutant cells. This comprehensive investigation of cell-cycle regulators reveals that multiple layers of regulation (transcription, protein stability, and proteasome targeting) affect protein expression dynamics during the cell cycle.

Long noncoding RNA repertoire and targeting by nuclear exosome, cytoplasmic exonuclease, and RNAi in fission yeast

Long noncoding RNAs (lncRNAs), which are longer than 200 nucleotides but often unstable, contribute a substantial and diverse portion to pervasive noncoding transcriptomes. Most lncRNAs are poorly annotated and understood, although several play important roles in gene regulation and diseases. Here we systematically uncover and analyze lncRNAs in Schizosaccharomyces pombe. Based on RNA-seq data from twelve RNA-processing mutants and nine physiological conditions, we identify 5775 novel lncRNAs, nearly 4× the previously annotated lncRNAs. The expression of most lncRNAs becomes strongly induced under the genetic and physiological perturbations, most notably during late meiosis. Most lncRNAs are cryptic and suppressed by three RNA-processing pathways: the nuclear exosome, cytoplasmic exonuclease, and RNAi. Double-mutant analyses reveal substantial coordination and redundancy among these pathways. We classify lncRNAs by their dominant pathway into cryptic unstable transcripts (CUTs), Xrn1-sensitive unstable transcripts (XUTs), and Dicer-sensitive unstable transcripts (DUTs). XUTs and DUTs are enriched for antisense lncRNAs, while CUTs are often bidirectional and actively translated. The cytoplasmic exonuclease, along with RNAi, dampens the expression of thousands of lncRNAs and mRNAs that become induced during meiosis. Antisense lncRNA expression mostly negatively correlates with sense mRNA expression in the physiological, but not the genetic conditions. Intergenic and bidirectional lncRNAs emerge from nucleosome-depleted regions, upstream of positioned nucleosomes. Our results highlight both similarities and differences to lncRNA regulation in budding yeast. This broad survey of the lncRNA repertoire and characteristics in S. pombe, and the interwoven regulatory pathways that target lncRNAs, provides a rich framework for their further functional analyses.

Ribosomopathies: There's strength in numbers

Ribosomopathies are a group of human disorders most commonly caused by ribosomal protein haploinsufficiency or defects in ribosome biogenesis. These conditions manifest themselves as physiological defects in specific cell and tissue types. We review current molecular models to explain ribosomopathies and attempt to reconcile the tissue specificity of these disorders with the ubiquitous requirement for ribosomes in all cells. Ribosomopathies as a group are diverse in their origins and clinical manifestations; we use the well-described Diamond-Blackfan anemia (DBA) as a specific example to highlight some common features. We discuss ribosome homeostasis as an overarching principle that governs the sensitivity of specific cells and tissue types to ribosomal protein mutations. Mathematical models and experimental insights rationalize how even subtle shifts in the availability of ribosomes, such as those created by ribosome haploinsufficiency, can drive messenger RNA-specific effects on protein expression. We discuss recently identified roles played by ribosome rescue and recycling factors in regulating ribosome homeostasis.

Simplification of Markov chains with infinite state space and the mathematical theory of random gene expression bursts

Here we develop an effective approach to simplify two-time-scale Markov chains with infinite state spaces by removal of states with fast leaving rates, which improves the simplification method of finite Markov chains. We introduce the concept of fast transition paths and show that the effective transitions of the reduced chain can be represented as the superposition of the direct transitions and the indirect transitions via all the fast transition paths. Furthermore, we apply our simplification approach to the standard Markov model of single-cell stochastic gene expression and provide a mathematical theory of random gene expression bursts. We give the precise mathematical conditions for the bursting kinetics of both mRNAs and proteins. It turns out that random bursts exactly correspond to the fast transition paths of the Markov model. This helps us gain a better understanding of the physics behind the bursting kinetics as an emergent behavior from the fundamental multiscale biochemical reaction kinetics of stochastic gene expression.

An Integrated Polysome Profiling and Ribosome Profiling Method to Investigate In Vivo Translatome

Recent advances in global translatome analysis technologies enable us to understand how translational regulation of gene expression modulates cellular functions. In this chapter, we present an integrated method to measure various aspects of translatome by polysome profiling and ribosome profiling using purified B cells. We standardized our protocols to directly compare the results from these two approaches. Parallel assessment of translatome with these two approaches can generate a comprehensive picture on how translational regulation determines protein output.

Reconstitution of Targeted Deadenylation by the Ccr4-Not Complex and the YTH Domain Protein Mmi1

Ccr4-Not is a conserved protein complex that shortens the 3' poly(A) tails of eukaryotic mRNAs to regulate transcript stability and translation into proteins. RNA-binding proteins are thought to facilitate recruitment of Ccr4-Not to certain mRNAs, but lack of an in-vitro-reconstituted system has slowed progress in understanding the mechanistic details of this specificity. Here, we generate a fully recombinant Ccr4-Not complex that removes poly(A) tails from RNA substrates. The intact complex is more active than the exonucleases alone and has an intrinsic preference for certain RNAs. The RNA-binding protein Mmi1 is highly abundant in preparations of native Ccr4-Not. We demonstrate a high-affinity interaction between recombinant Ccr4-Not and Mmi1. Using in vitro assays, we show that Mmi1 accelerates deadenylation of target RNAs. Together, our results support a model whereby both RNA-binding proteins and the sequence context of mRNAs influence deadenylation rate to regulate gene expression.

Cis-regulatory elements explain most of the mRNA stability variation across genes in yeast

The stability of mRNA is one of the major determinants of gene expression. Although a wealth of sequence elements regulating mRNA stability has been described, their quantitative contributions to half-life are unknown. Here, we built a quantitative model for Saccharomyces cerevisiae based on functional mRNA sequence features that explains 59% of the half-life variation between genes and predicts half-life at a median relative error of 30%. The model revealed a new destabilizing 3' UTR motif, ATATTC, which we functionally validated. Codon usage proves to be the major determinant of mRNA stability. Nonetheless, single-nucleotide variations have the largest effect when occurring on 3' UTR motifs or upstream AUGs. Analyzing mRNA half-life data of 34 knockout strains showed that the effect of codon usage not only requires functional decapping and deadenylation, but also the 5'-to-3' exonuclease Xrn1, the nonsense-mediated decay genes, but not no-go decay. Altogether, this study quantitatively delineates the contributions of mRNA sequence features on stability in yeast, reveals their functional dependencies on degradation pathways, and allows accurate prediction of half-life from mRNA sequence.

Ribosome reinitiation can explain length-dependent translation of messenger RNA

Models of mRNA translation usually presume that transcripts are linear; upon reaching the end of a transcript each terminating ribosome returns to the cytoplasmic pool before initiating anew on a different transcript. A consequence of linear models is that faster translation of a given mRNA is unlikely to generate more of the encoded protein, particularly at low ribosome availability. Recent evidence indicates that eukaryotic mRNAs are circularized, potentially allowing terminating ribosomes to preferentially reinitiate on the same transcript. Here we model the effect of ribosome reinitiation on translation and show that, at high levels of reinitiation, protein synthesis rates are dominated by the time required to translate a given transcript. Our model provides a simple mechanistic explanation for many previously enigmatic features of eukaryotic translation, including the negative correlation of both ribosome densities and protein abundance on transcript length, the importance of codon usage in determining protein synthesis rates, and the negative correlation between transcript length and both codon adaptation and 5' mRNA folding energies. In contrast to linear models where translation is largely limited by initiation rates, our model reveals that all three stages of translation-initiation, elongation, and termination/reinitiation-determine protein synthesis rates even at low ribosome availability.

Nuclear Compartmentalization Contributes to Stage-Specific Gene Expression Control in Trypanosoma cruzi

In the protozoan parasite Trypanosoma cruzi, as in other trypanosomatids, transcription of protein coding genes occurs in a constitutive fashion, producing large polycistronic transcription units. These units are composed of non-functionally related genes which are pervasively processed to yield each mRNA. Therefore, post-transcriptional processes are crucial to regulate gene expression. Considering that nuclear compartmentalization could contribute to gene expression regulation, we comparatively studied the nuclear, cytoplasmic and whole cell transcriptomes of the non-infective epimastigote stage of T. cruzi, using RNA-Seq. We found that the cytoplasmic transcriptome tightly correlates with the whole cell transcriptome and both equally correlate with the proteome. Nonetheless, 1,200 transcripts showed differential abundance between the nuclear and cytoplasmic fractions. For the genes with transcript content augmented in the nucleus, significant structural and compositional differences were found. The analysis of the reported epimastigote translatome and proteome, revealed scarce ribosome footprints and encoded proteins for them. Ontology analyses unveiled that many of these genes are distinctive of other parasite life-cycle stages. Finally, the relocalization of transcript abundance in the metacyclic trypomastigote infective stage was confirmed for specific genes. While gene expression is strongly dependent on transcript steady-state level, we here highlight the importance of the distribution of transcripts abundance between compartments in T. cruzi. Particularly, we show that nuclear compartmentation is playing an active role in the developmental stage determination preventing off-stage expression.

Unsupervised detection of regulatory gene expression information in different genomic regions enables gene expression ranking

Background

The regulation of all gene expression steps (e.g., Transcription, RNA processing, Translation, and mRNA Degradation) is known to be primarily encoded in different parts of genes and in genomic regions in proximity to genes (e.g., promoters, untranslated regions, coding regions, introns, etc.). However, the entire gene expression codes and the genomic regions where they are encoded are still unknown.

Results

Here, we employ an unsupervised approach to estimate the concentration of gene expression codes in different non-coding parts of genes and transcripts, such as introns and untranslated regions, focusing on three model organisms (Escherichia coli, Saccharomyces cerevisiae, and Schizosaccharomyces pombe). Our analyses support the conjecture that regions adjacent to the beginning and end of ORFs and the beginning and end of introns tend to include higher concentration of gene expression information relatively to regions further away. In addition, we report the exact regions with elevated concentration of gene expression codes. Furthermore, we demonstrate that the concentration of these codes in different genetic regions is correlated with the expression levels of the corresponding genes, and with splicing efficiency measurements and meiotic stage gene expression measurements in S. cerevisiae.

Conclusion

We suggest that these discoveries improve our understanding of gene expression regulation and evolution; they can also be used for developing improved models of genome/gene evolution and for engineering gene expression in various biotechnological and synthetic biology applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1497-z) contains supplementary material, which is available to authorized users.

Analysis of Ribosome-Associated mRNAs in Rice Reveals the Importance of Transcript Size and GC Content in Translation

Gene expression is controlled at transcriptional and post-transcriptional levels including decoding of messenger RNA (mRNA) into polypeptides via ribosome-mediated translation. Translational regulation has been intensively studied in the model dicot plant Arabidopsis thaliana, and in this study, we assessed the translational status [proportion of steady-state mRNA associated with ribosomes] of mRNAs by Translating Ribosome Affinity Purification followed by mRNA-sequencing (TRAP-seq) in rice (Oryza sativa), a model monocot plant and the most important food crop. A survey of three tissues found that most transcribed rice genes are translated whereas few transposable elements are associated with ribosomes. Genes with short and GC-rich coding regions are overrepresented in ribosome-associated mRNAs, suggesting that the GC-richness characteristic of coding sequences in grasses may be an adaptation that favors efficient translation. Transcripts with retained introns and extended 5′ untranslated regions are underrepresented on ribosomes, and rice genes belonging to different evolutionary lineages exhibited differential enrichment on the ribosomes that was associated with GC content. Genes involved in photosynthesis and stress responses are preferentially associated with ribosomes, whereas genes in epigenetic regulation pathways are the least enriched on ribosomes. Such variation is more dramatic in rice than that in Arabidopsis and is correlated with the wide variation of GC content of transcripts in rice. Taken together, variation in the translation status of individual transcripts reflects important mechanisms of gene regulation, which may have a role in evolution and diversification.

Analysis of the association between codon optimality and mRNA stability in Schizosaccharomyces pombe

Background

Recent experiments have shown that codon optimality is a major determinant of mRNA stability in Saccharomyces cerevisiae and that this phenomenon may be conserved in Escherichia coli and some metazoans, although work in Neurospora crassa is not consistent with this model.

Results

We examined the association between codon optimality and mRNA stability in the fission yeast Schizosaccharomyces pombe. Our analysis revealed the following points. First, we observe a genome-wide association between codon optimality and mRNA stability also in S. pombe, suggesting evolutionary conservation of the phenomenon. Second, in both S. pombe and S. cerevisiae, mRNA synthesis rates are also correlated at the genome-wide analysis with codon optimality, suggesting that the long-appreciated association between codon optimality and mRNA abundance is due to regulation of both mRNA synthesis and degradation. However, when we examined correlation of codon optimality and either mRNA half-lives or synthesis rates controlling for mRNA abundance, codon optimality was still positively correlated with mRNA half-lives in S. cerevisiae, but the association was no longer significant for mRNA half-lives in S. pombe or for synthesis rates in either organism. This illustrates how only the pairwise analysis of multiple correlating variables may limit these types of analyses. Finally, in S. pombe, codon optimality is associated with known DNA/RNA sequence motifs that are associated with mRNA production/stability, suggesting these two features have been under similar selective pressures for optimal gene expression.

Conclusions

Consistent with the emerging body of studies, this study suggests that the association between codon optimality and mRNA stability may be a broadly conserved phenomenon. It also suggests that the association can be explained at least in part by independent adaptations of codon optimality and other transcript features for elevated expression during evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3237-6) contains supplementary material, which is available to authorized users.

Label-free LC-MS/MS shotgun proteomics to investigate the anti-inflammatory effect of rCC16

Clara cell protein (CC16) is an anti-inflammatory protein, which is expressed in the airway epithelium. It is involved in the development of airway inflammatory diseases, including chronic obstructive pulmonary disease and asthma. However, the exact molecular mechanism underlying its anti‑inflammatory action remains to be fully elucidated. The aim of the present study was to define the protein profiles of the anti‑inflammatory effect of CC16 in lipopolysaccharide (LPS)‑treated rat tracheal epithelial (RTE) cells using shotgun proteomics. Protein extracts were obtained from control RTE cells, RTE cells treated with LPS and RTE cells treated with LPS and recombinant CC16 (rCC16). Subsequent label‑free quantification and bioinformatics analyses identified 12 proteins that were differentially expressed in the three treatment groups as a cluster of five distinct groups according to their molecular functions. Five of the twelve proteins were revealed to be associated with the cytoskeleton: Matrix metalloproteinase‑9, myosin heavy chain 10, actin‑related protein‑3 homolog, elongation factor 1‑α‑1 (EF‑1‑α‑1), and acidic ribosomal phosphoprotein P0. Five of the twelve proteins were associated with cellular proliferation: DNA‑dependent protein kinase catalytic subunit, EF‑1‑α‑1, tyrosine 3‑monooxygenase, caspase recruitment domain (CARD) protein 12 and adenosylhomocysteinase (SAHH) 3. Three proteins were associated with gene regulation: EF‑1‑α‑1, SAHH 3 and acidic ribosomal phosphoprotein P0. Three proteins were associated with inflammation: Tyrosine 3‑monooxygenase, CARD protein 12 and statin‑related protein. ATPase (H+‑transporting, V1 subunit A, isoform 1) was revealed to be associated with energy metabolism, and uridine diphosphate glycosyltransferase 1 family polypeptide A8 with drug metabolism and detoxification. The identified proteins were further validated using reverse transcription‑quantitative polymerase chain reaction. These protein profiles, and their interacting protein network, may facilitate the elucidation of the molecular mechanisms underlying the anti‑inflammatory effects of CC16.

Predictive biophysical modeling and understanding of the dynamics of mRNA translation and its evolution

mRNA translation is the fundamental process of decoding the information encoded in mRNA molecules by the ribosome for the synthesis of proteins. The centrality of this process in various biomedical disciplines such as cell biology, evolution and biotechnology, encouraged the development of dozens of mathematical and computational models of translation in recent years. These models aimed at capturing various biophysical aspects of the process. The objective of this review is to survey these models, focusing on those based and/or validated on real large-scale genomic data. We consider aspects such as the complexity of the models, the biophysical aspects they regard and the predictions they may provide. Furthermore, we survey the central systems biology discoveries reported on their basis. This review demonstrates the fundamental advantages of employing computational biophysical translation models in general, and discusses the relative advantages of the different approaches and the challenges in the field.

Active Center Control of Termination by RNA Polymerase III and tRNA Gene Transcription Levels In Vivo

The ability of RNA polymerase (RNAP) III to efficiently recycle from termination to reinitiation is critical for abundant tRNA production during cellular proliferation, development and cancer. Yet understanding of the unique termination mechanisms used by RNAP III is incomplete, as is its link to high transcription output. We used two tRNA-mediated suppression systems to screen for Rpc1 mutants with gain- and loss- of termination phenotypes in S. pombe. 122 point mutation mutants were mapped to a recently solved 3.9 Å structure of yeast RNAP III elongation complex (EC); they cluster in the active center bridge helix and trigger loop, as well as the pore and funnel, the latter of which indicate involvement of the RNA cleavage domain of the C11 subunit in termination. Purified RNAP III from a readthrough (RT) mutant exhibits increased elongation rate. The data strongly support a kinetic coupling model in which elongation rate is inversely related to termination efficiency. The mutants exhibit good correlations of terminator RT in vitro and in vivo, and surprisingly, amounts of transcription in vivo. Because assessing in vivo transcription can be confounded by various parameters, we used a tRNA reporter with a processing defect and a strong terminator. By ruling out differences in RNA decay rates, the data indicate that mutants with the RT phenotype synthesize more RNA than wild type cells, and than can be accounted for by their increased elongation rate. Finally, increased activity by the mutants appears unrelated to the RNAP III repressor, Maf1. The results show that the mobile elements of the RNAP III active center, including C11, are key determinants of termination, and that some of the mutations activate RNAP III for overall transcription. Similar mutations in spontaneous cancer suggest this as an unforeseen mechanism of RNAP III activation in disease.

Complementary Post Transcriptional Regulatory Information is Detected by PUNCH-P and Ribosome Profiling

Two novel approaches were recently suggested for genome-wide identification of protein aspects synthesized at a given time. Ribo-Seq is based on sequencing all the ribosome protected mRNA fragments in a cell, while PUNCH-P is based on mass-spectrometric analysis of only newly synthesized proteins. Here we describe the first Ribo-Seq/PUNCH-P comparison via the analysis of mammalian cells during the cell-cycle for detecting relevant differentially expressed genes between G1 and M phase. Our analyses suggest that the two approaches significantly overlap with each other. However, we demonstrate that there are biologically meaningful proteins/genes that can be detected to be post-transcriptionally regulated during the mammalian cell cycle only by each of the approaches, or their consolidation. Such gene sets are enriched with proteins known to be related to intra-cellular signalling pathways such as central cell cycle processes, central gene expression regulation processes, processes related to chromosome segregation, DNA damage, and replication, that are post-transcriptionally regulated during the mammalian cell cycle. Moreover, we show that combining the approaches better predicts steady state changes in protein abundance. The results reported here support the conjecture that for gaining a full post-transcriptional regulation picture one should integrate the two approaches.

AnGeLi: A Tool for the Analysis of Gene Lists from Fission Yeast

Genome-wide assays and screens typically result in large lists of genes or proteins. Enrichments of functional or other biological properties within such lists can provide valuable insights and testable hypotheses. To systematically detect these enrichments can be challenging and time-consuming, because relevant data to compare against query gene lists are spread over many different sources. We have developed AnGeLi (Analysis of Gene Lists), an intuitive, integrated web-tool for comprehensive and customized interrogation of gene lists from the fission yeast, Schizosaccharomyces pombe. AnGeLi searches for significant enrichments among multiple qualitative and quantitative information sources, including gene and phenotype ontologies, genetic and protein interactions, numerous features of genes, transcripts, translation, and proteins such as copy numbers, chromosomal positions, genetic diversity, RNA polymerase II and ribosome occupancy, localization, conservation, half-lives, domains, and molecular weight among others, as well as diverse sets of genes that are co-regulated or lead to the same phenotypes when mutated. AnGeLi uses robust statistics which can be tailored to specific needs. It also provides the option to upload user-defined gene sets to compare against the query list. Through an integrated data submission form, AnGeLi encourages the community to contribute additional curated gene lists to further increase the usefulness of this resource and to get the most from the ever increasing large-scale experiments. AnGeLi offers a rigorous yet flexible statistical analysis platform for rich insights into functional enrichments and biological context for query gene lists, thus providing a powerful exploratory tool through which S. pombe researchers can uncover fresh perspectives and unexpected connections from genomic data. AnGeLi is freely available at: www.bahlerlab.info/AnGeLi.

Genome-Wide Analysis of PAPS1-Dependent Polyadenylation Identifies Novel Roles for Functionally Specialized Poly(A) Polymerases in Arabidopsis thaliana

The poly(A) tail at 3’ ends of eukaryotic mRNAs promotes their nuclear export, stability and translational efficiency, and changes in its length can strongly impact gene expression. The Arabidopsis thaliana genome encodes three canonical nuclear poly(A) polymerases, PAPS1, PAPS2 and PAPS4. As shown by their different mutant phenotypes, these three isoforms are functionally specialized, with PAPS1 modifying organ growth and suppressing a constitutive immune response. However, the molecular basis of this specialization is largely unknown. Here, we have estimated poly(A)-tail lengths on a transcriptome-wide scale in wild-type and paps1 mutants. This identified categories of genes as particularly strongly affected in paps1 mutants, including genes encoding ribosomal proteins, cell-division factors and major carbohydrate-metabolic proteins. We experimentally verified two novel functions of PAPS1 in ribosome biogenesis and redox homoeostasis that were predicted based on the analysis of poly(A)-tail length changes in paps1 mutants. When overlaying the PAPS1-dependent effects observed here with coexpression analysis based on independent microarray data, the two clusters of transcripts that are most closely coexpressed with PAPS1 show the strongest change in poly(A)-tail length and transcript abundance in paps1 mutants in our analysis. This suggests that their coexpression reflects at least partly the preferential polyadenylation of these transcripts by PAPS1 versus the other two poly(A)-polymerase isoforms. Thus, transcriptome-wide analysis of poly(A)-tail lengths identifies novel biological functions and likely target transcripts for polyadenylation by PAPS1. Data integration with large-scale co-expression data suggests that changes in the relative activities of the isoforms are used as an endogenous mechanism to co-ordinately modulate plant gene expression.

The poly(A) tail of eukaryotic mRNAs promotes export from the nucleus, translation in the cytoplasm and stability of the mRNA, and changes in poly(A)-tail length can strongly impact on gene expression. The Arabidopsis thaliana genome encodes three nuclear canonical poly(A) polymerases (PAPS1, PAPS2, PAPS4) that fulfill different functions, presumably by preferentially polyadenylating certain subpopulations of pre-mRNAs. Here, we use a fractionation-based technique to assess the transcriptome-wide impact of reduced PAPS1 activity and identify functional classes of transcripts that are particularly sensitive to reduced PAPS1 activity. Analysis of these transcripts identifies two novel biological functions for PAPS1 in ribosome biogenesis and in redox homeostasis that we confirm experimentally. By overlaying our results with information about genome-wide co-expression, we demonstrate that genes co-expressed with PAPS1 are the most strongly affected in terms of poly(A)-tail length and total-abundance changes in the paps1 mutants. This provides strong evidence that the co-expression of these genes with PAPS1 that is seen across thousands of microarrays is at least partly caused by altered activity of the PAPS1 isoform, suggesting that the plant indeed uses modulation of the balance of isoform activities to coordinately regulate the expression of groups of genes.

PAT-seq: a method to study the integration of 3'-UTR dynamics with gene expression in the eukaryotic transcriptome

A major objective of systems biology is to quantitatively integrate multiple parameters from genome-wide measurements. To integrate gene expression with dynamics in poly(A) tail length and adenylation site, we developed a targeted next-generation sequencing approach, Poly(A)-Test RNA-sequencing. PAT-seq returns (i) digital gene expression, (ii) polyadenylation site/s, and (iii) the polyadenylation-state within and between eukaryotic transcriptomes. PAT-seq differs from previous 3' focused RNA-seq methods in that it depends strictly on 3' adenylation within total RNA samples and that the full-native poly(A) tail is included in the sequencing libraries. Here, total RNA samples from budding yeast cells were analyzed to identify the intersect between adenylation state and gene expression in response to loss of the major cytoplasmic deadenylase Ccr4. Furthermore, concordant changes to gene expression and adenylation-state were demonstrated in the classic Crabtree-Warburg metabolic shift. Because all polyadenylated RNA is interrogated by the approach, alternative adenylation sites, noncoding RNA and RNA-decay intermediates were also identified. Most important, the PAT-seq approach uses standard sequencing procedures, supports significant multiplexing, and thus replication and rigorous statistical analyses can for the first time be brought to the measure of 3'-UTR dynamics genome wide.

Condensin targets and reduces unwound DNA structures associated with transcription in mitotic chromosome condensation

Chromosome condensation is a hallmark of mitosis in eukaryotes and is a prerequisite for faithful segregation of genetic material to daughter cells. Here we show that condensin, which is essential for assembling condensed chromosomes, helps to preclude the detrimental effects of gene transcription on mitotic condensation. ChIP-seq profiling reveals that the fission yeast condensin preferentially binds to active protein-coding genes in a transcription-dependent manner during mitosis. Pharmacological and genetic attenuation of transcription largely rescue bulk chromosome segregation defects observed in condensin mutants. We also demonstrate that condensin is associated with and reduces unwound DNA segments generated by transcription, providing a direct link between an in vitro activity of condensin and its in vivo function. The human condensin isoform condensin I also binds to unwound DNA regions at the transcription start sites of active genes, implying that our findings uncover a fundamental feature of condensin complexes.

Chromosome condensation is a prerequisite for faithful segregation of chromosomes to daughter cells. Here, the authors show that the condensin complex binds to protein-coding genes in a transcription-dependent manner during condensation, and reduces unwound DNA segments generated by transcription.

Germination Potential of Dormant and Nondormant Arabidopsis Seeds Is Driven by Distinct Recruitment of Messenger RNAs to Polysomes

Dormancy is a complex evolutionary trait that temporally prevents seed germination, thus allowing seedling growth at a favorable season. High-throughput analyses of transcriptomes have led to significant progress in understanding the molecular regulation of this process, but the role of posttranscriptional mechanisms has received little attention. In this work, we have studied the dynamics of messenger RNA association with polysomes and compared the transcriptome with the translatome in dormant and nondormant seeds of Arabidopsis (Arabidopsis thaliana) during their imbibition at 25 °C in darkness, a temperature preventing germination of dormant seeds only. DNA microarray analysis revealed that 4,670 and 7,028 transcripts were differentially abundant in dormant and nondormant seeds in the transcriptome and the translatome, respectively. We show that there is no correlation between transcriptome and translatome and that germination regulation is also largely translational, implying a selective and dynamic recruitment of messenger RNAs to polysomes in both dormant and nondormant seeds. The study of 5' untranslated region features revealed that GC content and the number of upstream open reading frames could play a role in selective translation occurring during germination. Gene Ontology clustering showed that the functions of polysome-associated transcripts differed between dormant and nondormant seeds and revealed actors in seed dormancy and germination. In conclusion, our results demonstrate the essential role of selective polysome loading in this biological process.

AtCCR4a and AtCCR4b are Involved in Determining the Poly(A) Length of Granule-bound starch synthase 1 Transcript and Modulating Sucrose and Starch Metabolism in Arabidopsis thaliana

Removing the poly(A) tail is the first and rate-limiting step of mRNA degradation and apparently an effective step not only for modulating mRNA stability but also for translation of many eukaryotic transcripts. Carbon catabolite repressor 4 (CCR4) has been identified as a major cytoplasmic deadenylase in Saccharomyces cerevisiae. The Arabidopsis thaliana homologs of the yeast CCR4, AtCCR4a and AtCCR4b, were identified by sequence-based analysis; however, their role and physiological significance in plants remain to be elucidated. In this study, we revealed that AtCCR4a and AtCCR4b are localized to cytoplasmic mRNA processing bodies, which are specific granules consisting of many enzymes involved in mRNA turnover. Double mutants of AtCCR4a and AtCCR4b exhibited tolerance to sucrose application but not to glucose. The levels of sucrose in the seedlings of the atccr4a/4b double mutants were reduced, whereas no difference was observed in glucose levels. Further, amylose levels were slightly but significantly increased in the atccr4a/4b double mutants. Consistent with this observation, we found that the transcript encoding granule-bound starch synthase 1 (GBSS1), which is responsible for amylose synthesis, is accumulated to a higher level in the atccr4a/4b double mutant plants than in the control plants. Moreover, we revealed that GBSS1 has a longer poly(A) tail in the double mutant than in the control plant, suggesting that AtCCR4a and AtCCR4b can influence the poly(A) length of transcripts related to starch metabolism. Our results collectively suggested that AtCCR4a and AtCCR4b are involved in sucrose and starch metabolism in A. thaliana.

Causal signals between codon bias, mRNA structure, and the efficiency of translation and elongation

Ribosome profiling data report on the distribution of translating ribosomes, at steady-state, with codon-level resolution. We present a robust method to extract codon translation rates and protein synthesis rates from these data, and identify causal features associated with elongation and translation efficiency in physiological conditions in yeast. We show that neither elongation rate nor translational efficiency is improved by experimental manipulation of the abundance or body sequence of the rare AGG tRNA. Deletion of three of the four copies of the heavily used ACA tRNA shows a modest efficiency decrease that could be explained by other rate-reducing signals at gene start. This suggests that correlation between codon bias and efficiency arises as selection for codons to utilize translation machinery efficiently in highly translated genes. We also show a correlation between efficiency and RNA structure calculated both computationally and from recent structure probing data, as well as the Kozak initiation motif, which may comprise a mechanism to regulate initiation.

Translatome profiling: methods for genome-scale analysis of mRNA translation

During the past decade, there has been a rapidly increased appreciation of the role of translation as a key regulatory node in gene expression. Thereby, the development of methods to infer the translatome, which refers to the entirety of mRNAs associated with ribosomes for protein synthesis, has facilitated the discovery of new principles and mechanisms of translation and expanded our view of the underlying logic of protein synthesis. Here, we review the three main methodologies for translatome analysis, and we highlight some of the recent discoveries made using each technique. We first discuss polysomal profiling, a classical technique that involves the separation of mRNAs depending on the number of bound ribosomes using a sucrose gradient, and which has been combined with global analysis tools such as DNA microarrays or high-throughput RNA sequencing to identify the RNAs in polysomal fractions. We then introduce ribosomal profiling, a recently established technique that enables the mapping of ribosomes along mRNAs at near-nucleotide resolution on a global scale. We finally refer to ribosome affinity purification techniques that are based on the cell-type-specific expression of tagged ribosomal proteins, allowing the capture of translatomes from specialized cells in organisms. We discuss the advantages and disadvantages of these three main techniques in the pursuit of defining the translatome, and we speculate about future developments.

In vitro analysis of RNA degradation catalyzed by deadenylase enzymes

In this chapter, we describe a method for purification and analysis of the enzymatic activity of deadenylase enzymes. Nearly all eukaryotic messenger RNAs are modified at the 3' end by the addition of an adenosine polymer: the poly-adenosine tail. The poly(A) tail plays a central role in protein expression and mRNA fate. The poly(A) tail promotes translation of the mRNA. Shortening of the poly(A) tail, referred to as deadenylation, reduces protein synthesis and initiates destruction of the mRNA. A specialized class of exoribonucleases, called deadenylase enzymes, carries out this process. Deadenylases are found throughout eukarya, but their functions remain largely unexplored. We present a detailed protocol to analyze deadenylase activity in vitro. First, recombinant deadenylase enzyme is over-expressed and purified from bacteria. Next, labeled RNA substrate is prepared. Deadenylation reactions are performed, and reaction products are analyzed by denaturing gel electrophoresis. Reaction rates are then determined quantitatively. Crucial controls and experimental parameters are described along with practical tips that promote success.

Using Klenow-mediated extension to measure poly(A)-tail length and position in the transcriptome

The poly(A)-tail that terminates most mRNA and many noncoding RNA is a convenient "hook" to isolate mRNA. However the length of this tail and its position within the primary RNA transcript can also hold diagnostic value for RNA metabolism. In general, mRNA with a long poly(A)-tail is well translated, whereas a short poly(A)-tail can indicate translational silencing. A short poly(A)-tail is also appended to RNA-decay intermediates via the TRAMP complex. A number of approaches have been developed to measure the length and position of the poly(A)-tail. Here, we describe a simple method to "tag" adenylated RNA using the native function of DNA polymerase I to extend an RNA primer on a DNA template in second-strand DNA synthesis. This function can be harnessed as a means to purify, visualize, and quantitate poly(A)-dynamics of individual RNA and the transcriptome en masse.

Genetic influences on translation in yeast

Heritable differences in gene expression between individuals are an important source of phenotypic variation. The question of how closely the effects of genetic variation on protein levels mirror those on mRNA levels remains open. Here, we addressed this question by using ribosome profiling to examine how genetic differences between two strains of the yeast S. cerevisiae affect translation. Strain differences in translation were observed for hundreds of genes. Allele specific measurements in the diploid hybrid between the two strains revealed roughly half as many cis-acting effects on translation as were observed for mRNA levels. In both the parents and the hybrid, most effects on translation were of small magnitude, such that the direction of an mRNA difference was typically reflected in a concordant footprint difference. The relative importance of cis and trans acting variation on footprint levels was similar to that for mRNA levels. There was a tendency for translation to cause larger footprint differences than expected given the respective mRNA differences. This is in contrast to translational differences between yeast species that have been reported to more often oppose than reinforce mRNA differences. Finally, we catalogued instances of premature translation termination in the two yeast strains and also found several instances where erroneous reference gene annotations lead to apparent nonsense mutations that in fact reside outside of the translated gene body. Overall, genetic influences on translation subtly modulate gene expression differences, and translation does not create strong discrepancies between genetic influences on mRNA and protein levels.

Individuals in a species differ from each other in many ways. For many traits, a fraction of this variation is genetic—it is caused by DNA sequence variants in the genome of each individual. Some of these variants influence traits by altering how much certain genes are expressed, i.e. how many mRNA and protein molecules are made in different individuals. Surprisingly, earlier work has found that the effects of genetic variants on mRNA and protein levels for the same genes appear to be very different. Many variants appeared to influence only mRNA (but not protein) levels, and vice versa. In this paper, we studied this question by using a technique called “ribosome profiling” to measure translation (the cellular process of reading mRNA molecules and synthesizing protein molecules) in two yeast strains. We found that the genetic differences between these two strains influence translation for hundreds of genes. Because most of these effects were small in magnitude, they explain at most a small fraction of the discrepancies between the effects of genetic variants on mRNA and protein levels.