Quantification of the yeast transcriptome by single-molecule sequencing

YeastFab Cloning of Toxic Genes and Protein Expression Optimization in Yeast

YeastFab is a Golden Gate-based cloning standard and parts repository. It is designed for modular, hierarchical assembly of transcription units and multi-gene assemblies for expression in Saccharomyces cerevisiae. This makes it a suitable toolbox to optimize the expression strength of heterologous genes in yeast. When cloning heterologous coding sequences into YeastFab vectors, in several cases we have observed toxicity to the cloning host Escherichia coli. The provided protocol details how to clone such toxic genes from multiple synthetic DNA fragments while adhering to the YeastFab standard. The presented cloning strategy includes a C-terminal FLAG tag that allows screening for constructs with a desired protein expression in yeast by western blot. The design allows scarlessly removing the tag through a Golden Gate reaction to facilitate cloning of expression constructs with the native, untagged transgene.

Polarized exocytosis and anionic phospholipid species implicated in the initiation of clathrin-mediated endocytosis

Understanding of the mechanisms that initiate clathrin-mediated endocytosis (CME) is incomplete. Recent studies in budding yeast identified the endocytic adaptor protein Yap1801/Yap1802 (budding yeast AP180) as a key CME factor that promotes CME initiation in daughter cells during polarized growth, but how Yap1801/2 is recruited preferentially to the plasma membrane of daughter cells is not clear. The only known cargos for Yap1801/2 in yeast are the synaptobrevins Snc1 and Snc2, which act as v-SNARES for exocytic vesicles arriving at the plasma membrane and are essential for polarized cell growth. In this study, we analyze the spatiotemporal dynamics of functional, fluorescently-tagged Snc1/2 expressed from their endogenous loci and provide evidence that, in concert with anionic phospholipids, Snc1/2 recruit Yap1801/2 preferentially to growing daughter cells. These findings suggest that the coincidence of anionic phospholipids and Snc1/2 facilitates CME initiation in growing daughter cells and directly links polarized CME to polarized secretion.

Exploiting phenotypic heterogeneity to improve production of glutathione by yeast

BACKGROUND

Gene expression noise (variation in gene expression among individual cells of a genetically uniform cell population) can result in heterogenous metabolite production by industrial microorganisms, with cultures containing both low- and high-producing cells. The presence of low-producing individuals may be a factor limiting the potential for high yields. This study tested the hypothesis that low-producing variants in yeast cell populations can be continuously counter-selected, to increase net production of glutathione (GSH) as an exemplar product.

RESULTS

A counter-selection system was engineered in Saccharomyces cerevisiae based on the known feedback inhibition of gamma-glutamylcysteine synthetase (GSH1) gene expression, which is rate limiting for GSH synthesis: the GSH1 ORF and the counter-selectable marker GAP1 were expressed under control of the TEF1 and GSH-regulated GSH1 promoters, respectively. An 18% increase in the mean cellular GSH level was achieved in cultures of the engineered strain supplemented with D-histidine to counter-select cells with high GAP1 expression (i.e. low GSH-producing cells). The phenotype was non-heritable and did not arise from a generic response to D-histidine, unlike that with certain other test-constructs prepared with alternative markers.

CONCLUSIONS

The results corroborate that the system developed here improves GSH production by targeting low-producing cells. This supports the potential for exploiting end-product/promoter interactions to enrich high-producing cells in phenotypically heterogeneous populations, in order to improve metabolite production by yeast.

Enrichment of rare codons at 5' ends of genes is a spandrel caused by evolutionary sequence turnover and does not improve translation

Previously, Tuller et al. found that the first 30-50 codons of the genes of yeast and other eukaryotes are slightly enriched for rare codons. They argued that this slowed translation, and was adaptive because it queued ribosomes to prevent collisions. Today, the translational speeds of different codons are known, and indeed rare codons are translated slowly. We re-examined this 5' slow translation 'ramp.' We confirm that 5' regions are slightly enriched for rare codons; in addition, they are depleted for downstream Start codons (which are fast), with both effects contributing to slow 5' translation. However, we also find that the 5' (and 3') ends of yeast genes are poorly conserved in evolution, suggesting that they are unstable and turnover relatively rapidly. When a new 5' end forms de novo, it is likely to include codons that would otherwise be rare. Because evolution has had a relatively short time to select against these codons, 5' ends are typically slightly enriched for rare, slow codons. Opposite to the expectation of Tuller et al., we show by direct experiment that genes with slowly translated codons at the 5' end are expressed relatively poorly, and that substituting faster synonymous codons improves expression. Direct experiment shows that slow codons do not prevent downstream ribosome collisions. Further informatic studies suggest that for natural genes, slow 5' ends are correlated with poor gene expression, opposite to the expectation of Tuller et al. Thus, we conclude that slow 5' translation is a 'spandrel'--a non-adaptive consequence of something else, in this case, the turnover of 5' ends in evolution, and it does not improve translation.

Meta-analysis of dispensable essential genes and their interactions with bypass suppressors

Genes have been historically classified as essential or non-essential based on their requirement for viability. However, genomic mutations can sometimes bypass the requirement for an essential gene, challenging the binary classification of gene essentiality. Such dispensable essential genes represent a valuable model for understanding the incomplete penetrance of loss-of-function mutations often observed in natural populations. Here, we compiled data from multiple studies on essential gene dispensability in Saccharomyces cerevisiae to comprehensively characterize these genes. In analyses spanning different evolutionary timescales, dispensable essential genes exhibited distinct phylogenetic properties compared with other essential and non-essential genes. Integration of interactions with suppressor genes that can bypass the gene essentiality revealed the high functional modularity of the bypass suppression network. Furthermore, dispensable essential and bypass suppressor gene pairs reflected simultaneous changes in the mutational landscape of S. cerevisiae strains. Importantly, species in which dispensable essential genes were non-essential tended to carry bypass suppressor mutations in their genomes. Overall, our study offers a comprehensive view of dispensable essential genes and illustrates how their interactions with bypass suppressors reflect evolutionary outcomes.

Differential gene expression response to acute and chronic Cytauzxoon felis infection in domestic cats (Felis catus)

Cytauxzoonosis is a severe tick transmitted protozoan disease of domestic cats, caused by Cytauxzoon felis. The disease is characterized by acute onset of high fever, depression, lethargy, inappentence, anorexia, icterus, dehydration, hemolytic anemia, and alteration of immune response. The aim of our study was to further detail the immune response of domestic cats to C. felis infection by comparing the differential expression of feline immune transcriptional elements during acute and chronic cytauxzoonosis. True single molecule sequencing (tSMS) was used to analyze the whole genome of acutely and chronically infected C. felis cats, focusing on the analysis of genes involved on the immune response. Two C. felis donor cats were infested with Amblyomma americanum nymphs, which after repletion were collected and kept in humidity chambers until they molted. The resulting A. americanum were randomly selected to infest three C. felis naïve principal cats. Infection of these cats was confirmed by nested PCR of the 18S rRNA C. felis gene and clinical signs. RNA was extracted from whole blood at different time points and used for tSMS analyses, the results revealed overexpression in transcripts involved in type I interferon signaling, cellular and cytokine responses during the acute stage of infection, while cell cycle, and metabolic processes were downregulated. Genes involved in cell adhesion increased their expression in the chronic infected cats, whereas inflammatory and apoptotic related genes were downregulated. This study provided information on the host immune response to C. felis in domestic cats, demonstrating that inflammatory, apoptotic, and cell adhesion are some of the pathways altered during acute and chronic infection.

Spontaneous deamination of cytosine to uracil is biased to the non-transcribed DNA strand in yeast

Transcription in Saccharomyces cerevisiae is associated with elevated mutation and this partially reflects enhanced damage of the corresponding DNA. Spontaneous deamination of cytosine to uracil leads to CG>TA mutations that provide a strand-specific read-out of damage in strains that lack the ability to remove uracil from DNA. Using the CAN1 forward mutation reporter, we found that C>T and G>A mutations, which reflect deamination of the non-transcribed and transcribed DNA strands, respectively, occurred at similar rates under low-transcription conditions. By contrast, the rate of C>T mutations was 3-fold higher than G>A mutations under high-transcription conditions, demonstrating biased deamination of the non-transcribed strand (NTS). The NTS is transiently single-stranded within the ∼15 bp transcription bubble, or a more extensive region of the NTS can be exposed as part of an R-loop that can form behind RNA polymerase. Neither the deletion of genes whose products restrain R-loop formation nor the over-expression of RNase H1, which degrades R-loops, reduced the biased deamination of the NTS, and no transcription-associated R-loop formation at CAN1 was detected. These results suggest that the NTS within the transcription bubble is a target for spontaneous deamination and likely other types of DNA damage.

Fungal biotechnology: From yesterday to tomorrow

Fungi have been used to better the lives of everyday people and unravel the mysteries of higher eukaryotic organisms for decades. However, comparing progress and development stemming from fungal research to that of human, plant, and bacterial research, fungi remain largely understudied and underutilized. Recent commercial ventures have begun to gain popularity in society, providing a new surge of interest in fungi, mycelia, and potential new applications of these organisms to various aspects of research. Biotechnological advancements in fungal research cannot occur without intensive amounts of time, investments, and research tool development. In this review, we highlight past breakthroughs in fungal biotechnology, discuss requirements to advance fungal biotechnology even further, and touch on the horizon of new breakthroughs with the highest potential to positively impact both research and society.

A dynamical stochastic model of yeast translation across the cell cycle

Translation is a central step in gene expression, however its quantitative and time-resolved regulation is poorly understood. We developed a discrete, stochastic model for protein translation in S. cerevisiae in a whole-transcriptome, single-cell context. A "base case" scenario representing an average cell highlights translation initiation rates as the main co-translational regulatory parameters. Codon usage bias emerges as a secondary regulatory mechanism through ribosome stalling. Demand for anticodons with low abundancy is shown to cause above-average ribosome dwelling times. Codon usage bias correlates strongly both with protein synthesis rates and elongation rates. Applying the model to a time-resolved transcriptome estimated by combining data from FISH and RNA-Seq experiments, it could be shown that increased total transcript abundance during the cell cycle decreases translation efficiency at single transcript level. Translation efficiency grouped by gene function shows highest values for ribosomal and glycolytic genes. Ribosomal proteins peak in S phase while glycolytic proteins rank highest in later cell cycle phases.

RNAseq analysis of treatment-dependent signaling changes during inflammation in a mouse cutaneous wound healing model

BACKGROUND

Despite proven therapeutic effects in inflammatory conditions, the specific mechanisms of phytochemical therapies are not well understood. The transcriptome effects of Traumeel (Tr14), a multicomponent natural product, and diclofenac, a non-selective cyclooxygenase (COX) inhibitor, were compared in a mouse cutaneous wound healing model to identify both known and novel pathways for the anti-inflammatory effect of plant-derived natural products.

METHODS

Skin samples from abraded mice were analyzed by single-molecule, amplification-free RNAseq transcript profiling at 7 points between 12 and 192 h after injury. Immediately after injury, the wounds were treated with either diclofenac, Tr14, or placebo control (n = 7 per group/time). RNAseq levels were compared between treatment and control at each time point using a systems biology approach.

RESULTS

At early time points (12-36 h), both control and Tr14-treated wounds showed marked increase in the inducible COX2 enzyme mRNA, while diclofenac-treated wounds did not. Tr14, in contrast, modulated lipoxygenase transcripts, especially ALOX12/15, and phospholipases involved in arachidonate metabolism. Notably, Tr14 modulated a group of cell-type specific markers, including the T cell receptor, that could be explained by an overarching effect on the type of cells that were recruited into the wound tissue.

CONCLUSIONS

Tr14 and diclofenac had very different effects on the COX/LOX synthetic pathway after cutaneous wounding. Tr14 allowed normal autoinduction of COX2 mRNA, but suppressed mRNA levels for key enzymes in the leukotriene synthetic pathway. Tr14 appeared to have a broad 'phytocellular' effect on the wound transcriptome by altering the balance of cell types present in the wound.

Investigation of RNA metabolism through large-scale genetic interaction profiling in yeast

Gene deletion and gene expression alteration can lead to growth defects that are amplified or reduced when a second mutation is present in the same cells. We performed 154 genetic interaction mapping (GIM) screens with query mutants related with RNA metabolism and estimated the growth rates of about 700 000 double mutant Saccharomyces cerevisiae strains. The tested targets included the gene deletion collection and 900 strains in which essential genes were affected by mRNA destabilization (DAmP). To analyze the results, we developed RECAP, a strategy that validates genetic interaction profiles by comparison with gene co-citation frequency, and identified links between 1471 genes and 117 biological processes. In addition to these large-scale results, we validated both enhancement and suppression of slow growth measured for specific RNA-related pathways. Thus, negative genetic interactions identified a role for the OCA inositol polyphosphate hydrolase complex in mRNA translation initiation. By analysis of suppressors, we found that Puf4, a Pumilio family RNA binding protein, inhibits ribosomal protein Rpl9 function, by acting on a conserved UGUAcauUA motif located downstream the stop codon of the RPL9B mRNA. Altogether, the results and their analysis should represent a useful resource for discovery of gene function in yeast.

Systematic analysis of bypass suppression of essential genes

Essential genes tend to be highly conserved across eukaryotes, but, in some cases, their critical roles can be bypassed through genetic rewiring. From a systematic analysis of 728 different essential yeast genes, we discovered that 124 (17%) were dispensable essential genes. Through whole-genome sequencing and detailed genetic analysis, we investigated the genetic interactions and genome alterations underlying bypass suppression. Dispensable essential genes often had paralogs, were enriched for genes encoding membrane-associated proteins, and were depleted for members of protein complexes. Functionally related genes frequently drove the bypass suppression interactions. These gene properties were predictive of essential gene dispensability and of specific suppressors among hundreds of genes on aneuploid chromosomes. Our findings identify yeast's core essential gene set and reveal that the properties of dispensable essential genes are conserved from yeast to human cells, correlating with human genes that display cell line-specific essentiality in the Cancer Dependency Map (DepMap) project.

A genome-scale yeast library with inducible expression of individual genes

The ability to switch a gene from off to on and monitor dynamic changes provides a powerful approach for probing gene function and elucidating causal regulatory relationships. Here, we developed and characterized YETI (Yeast Estradiol strains with Titratable Induction), a collection in which > 5,600 yeast genes are engineered for transcriptional inducibility with single-gene precision at their native loci and without plasmids. Each strain contains SGA screening markers and a unique barcode, enabling high-throughput genetics. We characterized YETI using growth phenotyping and BAR-seq screens, and we used a YETI allele to identify the regulon of Rof1, showing that it acts to repress transcription. We observed that strains with inducible essential genes that have low native expression can often grow without inducer. Analysis of data from eukaryotic and prokaryotic systems shows that native expression is a variable that can bias promoter-perturbing screens, including CRISPRi. We engineered a second expression system, Z3 EB42, that gives lower expression than Z3 EV, a feature enabling conditional activation and repression of lowly expressed essential genes that grow without inducer in the YETI library.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

Promoter Architecture and Promoter Engineering in Saccharomyces cerevisiae

Promoters play an essential role in the regulation of gene expression for fine-tuning genetic circuits and metabolic pathways in Saccharomyces cerevisiae (S. cerevisiae). However, native promoters in S. cerevisiae have several limitations which hinder their applications in metabolic engineering. These limitations include an inadequate number of well-characterized promoters, poor dynamic range, and insufficient orthogonality to endogenous regulations. Therefore, it is necessary to perform promoter engineering to create synthetic promoters with better properties. Here, we review recent advances related to promoter architecture, promoter engineering and synthetic promoter applications in S. cerevisiae. We also provide a perspective of future directions in this field with an emphasis on the recent advances of machine learning based promoter designs.

Systematic genetics and single-cell imaging reveal widespread morphological pleiotropy and cell-to-cell variability

Our ability to understand the genotype-to-phenotype relationship is hindered by the lack of detailed understanding of phenotypes at a single-cell level. To systematically assess cell-to-cell phenotypic variability, we combined automated yeast genetics, high-content screening and neural network-based image analysis of single cells, focussing on genes that influence the architecture of four subcellular compartments of the endocytic pathway as a model system. Our unbiased assessment of the morphology of these compartments-endocytic patch, actin patch, late endosome and vacuole-identified 17 distinct mutant phenotypes associated with ~1,600 genes (~30% of all yeast genes). Approximately half of these mutants exhibited multiple phenotypes, highlighting the extent of morphological pleiotropy. Quantitative analysis also revealed that incomplete penetrance was prevalent, with the majority of mutants exhibiting substantial variability in phenotype at the single-cell level. Our single-cell analysis enabled exploration of factors that contribute to incomplete penetrance and cellular heterogeneity, including replicative age, organelle inheritance and response to stress.

Novel approach reveals genomic landscapes of single-strand DNA breaks with nucleotide resolution in human cells

Single-strand breaks (SSBs) represent the major form of DNA damage, yet techniques to map these lesions genome-wide with nucleotide-level precision are limited. Here, we present a method, termed SSiNGLe, and demonstrate its utility to explore the distribution and dynamic changes in genome-wide SSBs in response to different biological and environmental stimuli. We validate SSiNGLe using two very distinct sequencing techniques and apply it to derive global profiles of SSBs in different biological states. Strikingly, we show that patterns of SSBs in the genome are non-random, specific to different biological states, enriched in regulatory elements, exons, introns, specific types of repeats and exhibit differential preference for the template strand between exons and introns. Furthermore, we show that breaks likely contribute to naturally occurring sequence variants. Finally, we demonstrate strong links between SSB patterns and age. Overall, SSiNGLe provides access to unexplored realms of cellular biology, not obtainable with current approaches.

Single-molecule real-time sequencing facilitates the analysis of transcripts and splice isoforms of anthers in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

Anther development has been extensively studied at the transcriptional level, but a systematic analysis of full-length transcripts on a genome-wide scale has not yet been published. Here, the Pacific Biosciences (PacBio) Sequel platform and next-generation sequencing (NGS) technology were combined to generate full-length sequences and completed structures of transcripts in anthers of Chinese cabbage.

RESULTS

Using single-molecule real-time sequencing (SMRT), a total of 1,098,119 circular consensus sequences (CCSs) were generated with a mean length of 2664 bp. More than 75% of the CCSs were considered full-length non-chimeric (FLNC) reads. After error correction, 725,731 high-quality FLNC reads were estimated to carry 51,501 isoforms from 19,503 loci, consisting of 38,992 novel isoforms from known genes and 3691 novel isoforms from novel genes. Of the novel isoforms, we identified 407 long non-coding RNAs (lncRNAs) and 37,549 open reading frames (ORFs). Furthermore, a total of 453,270 alternative splicing (AS) events were identified and the majority of AS models in anther were determined to be approximate exon skipping (XSKIP) events. Of the key genes regulated during anther development, AS events were mainly identified in the genes SERK1, CALS5, NEF1, and CESA1/3. Additionally, we identified 104 fusion transcripts and 5806 genes that had alternative polyadenylation (APA).

CONCLUSIONS

Our work demonstrated the transcriptome diversity and complexity of anther development in Chinese cabbage. The findings provide a basis for further genome annotation and transcriptome research in Chinese cabbage.

Transcriptome profiling of mouse samples using nanopore sequencing of cDNA and RNA molecules

Our vision of DNA transcription and splicing has changed dramatically with the introduction of short-read sequencing. These high-throughput sequencing technologies promised to unravel the complexity of any transcriptome. Generally gene expression levels are well-captured using these technologies, but there are still remaining caveats due to the limited read length and the fact that RNA molecules had to be reverse transcribed before sequencing. Oxford Nanopore Technologies has recently launched a portable sequencer which offers the possibility of sequencing long reads and most importantly RNA molecules. Here we generated a full mouse transcriptome from brain and liver using the Oxford Nanopore device. As a comparison, we sequenced RNA (RNA-Seq) and cDNA (cDNA-Seq) molecules using both long and short reads technologies and tested the TeloPrime preparation kit, dedicated to the enrichment of full-length transcripts. Using spike-in data, we confirmed that expression levels are efficiently captured by cDNA-Seq using short reads. More importantly, Oxford Nanopore RNA-Seq tends to be more efficient, while cDNA-Seq appears to be more biased. We further show that the cDNA library preparation of the Nanopore protocol induces read truncation for transcripts containing internal runs of T's. This bias is marked for runs of at least 15 T's, but is already detectable for runs of at least 9 T's and therefore concerns more than 20% of expressed transcripts in mouse brain and liver. Finally, we outline that bioinformatics challenges remain ahead for quantifying at the transcript level, especially when reads are not full-length. Accurate quantification of repeat-associated genes such as processed pseudogenes also remains difficult, and we show that current mapping protocols which map reads to the genome largely over-estimate their expression, at the expense of their parent gene.

Promoter-dependent nuclear RNA degradation ensures cell cycle-specific gene expression

Cell cycle progression depends on phase-specific gene expression. Here we show that the nuclear RNA degradation machinery plays a lead role in promoting cell cycle-dependent gene expression by triggering promoter-dependent co-transcriptional RNA degradation. Single molecule quantification of RNA abundance in different phases of the cell cycle indicates that relative curtailment of gene expression in certain phases is attained even when transcription is not completely inhibited. When nuclear ribonucleases are deleted, transcription of the Saccharomyces cerevisiae G1-specific axial budding gene AXL2 is detected throughout the cell cycle and its phase-specific expression is lost. Promoter replacement abolished cell cycle-dependent RNA degradation and rendered the RNA insensitive to the deletion of nuclear ribonucleases. Together the data reveal a model of gene regulation whereby RNA abundance is controlled by promoter-dependent induction of RNA degradation.

Assessment of Bioleaching Microbial Community Structure and Function Based on Next-Generation Sequencing Technologies

It is widely known that bioleaching microorganisms have to cope with the complex extreme environment in which microbial ecology relating to community structure and function varies across environmental types. However, analyses of microbial ecology of bioleaching bacteria is still a challenge. To address this challenge, numerous technologies have been developed. In recent years, high-throughput sequencing technologies enabling comprehensive sequencing analysis of cellular RNA and DNA within the reach of most laboratories have been added to the toolbox of microbial ecology. The next-generation sequencing technology allowing processing DNA sequences can produce available draft genomic sequences of more bioleaching bacteria, which provides the opportunity to predict models of genetic and metabolic potential of bioleaching bacteria and ultimately deepens our understanding of bioleaching microorganism. High-throughput sequencing that focuses on targeted phylogenetic marker 16S rRNA has been effectively applied to characterize the community diversity in an ore leaching environment. RNA-seq, another application of high-throughput sequencing to profile RNA, can be for both mapping and quantifying transcriptome and has demonstrated a high efficiency in quantifying the changing expression level of each transcript under different conditions. It has been demonstrated as a powerful tool for dissecting the relationship between genotype and phenotype, leading to interpreting functional elements of the genome and revealing molecular mechanisms of adaption. This review aims to describe the high-throughput sequencing approach for bioleaching environmental microorganisms, particularly focusing on its application associated with challenges.

Decoupling Yeast Cell Division and Stress Defense Implicates mRNA Repression in Translational Reallocation during Stress

Stress tolerance and rapid growth are often competing interests in cells. Upon severe environmental stress, many organisms activate defense systems concurrent with growth arrest. There has been debate as to whether aspects of the stress-activated transcriptome are regulated by stress or an indirect byproduct of reduced proliferation. For example, stressed Saccharomyces cerevisiae cells mount a common gene expression program called the environmental stress response (ESR) [1] comprised of ∼300 induced (iESR) transcripts involved in stress defense and ∼600 reduced (rESR) mRNAs encoding ribosomal proteins (RPs) and ribosome biogenesis factors (RiBi) important for division. Because ESR activation also correlates with reduced growth rate in nutrient-restricted chemostats and prolonged G1 in slow-growing mutants, an alternate proposal is that the ESR is simply a consequence of reduced division [2-5]. A major challenge is that past studies did not separate effects of division arrest and stress defense; thus, the true responsiveness of the ESR-and the purpose of stress-dependent rESR repression in particular-remains unclear. Here, we decoupled cell division from the stress response by following transcriptome, proteome, and polysome changes in arrested cells responding to acute stress. We show that the ESR cannot be explained by changes in growth rate or cell-cycle phase during stress acclimation. Instead, failure to repress rESR transcripts reduces polysome association of induced transcripts, delaying production of their proteins. Our results suggest that stressed cells alleviate competition for translation factors by removing mRNAs and ribosomes from the translating pool, directing translational capacity toward induced transcripts to accelerate protein production.

Venomics: A Mini-Review

Venomics is the integration of proteomic, genomic and transcriptomic approaches to study venoms. Advances in these approaches have enabled increasingly more comprehensive analyses of venoms to be carried out, overcoming to some extent the limitations imposed by the complexity of the venoms and the small quantities that are often available. Advances in bioinformatics and high-throughput functional assay screening approaches have also had a significant impact on venomics. A combination of all these techniques is critical for enhancing our knowledge on the complexity of venoms and their potential therapeutic and agricultural applications. Here we highlight recent advances in these fields and their impact on venom analyses.

On the design and prospects of direct RNA sequencing

Throughout the past nearly a decade, the application of high-throughput sequencing to RNA molecules in the form of RNA sequencing (RNA-seq) and its many variations has revolutionized transcriptomic studies by enabling researchers to take a simultaneously deep and truly global look into the transcriptome. However, there is still considerable scope for improvement on RNA-seq data in its current form, primarily because of the short-read nature of the dominant sequencing technologies, which prevents the completely reliable reconstruction and quantification of full-length transcripts, and the sequencing library building protocols used, which introduce various distortions in the final data sets. The ideal approach toward resolving these remaining issues would involve the direct amplification-free sequencing of full-length RNA molecules. This has recently become practical with the advent of nanopore sequencing, which raises the possibility of yet another revolution in transcriptomics. I discuss the design considerations to be taken into account, the technical challenges that need to be addressed and the biological questions these advances can be expected to resolve.

Highly parallel direct RNA sequencing on an array of nanopores

Sequencing the RNA in a biological sample can unlock a wealth of information, including the identity of bacteria and viruses, the nuances of alternative splicing or the transcriptional state of organisms. However, current methods have limitations due to short read lengths and reverse transcription or amplification biases. Here we demonstrate nanopore direct RNA-seq, a highly parallel, real-time, single-molecule method that circumvents reverse transcription or amplification steps. This method yields full-length, strand-specific RNA sequences and enables the direct detection of nucleotide analogs in RNA.

Unification of Protein Abundance Datasets Yields a Quantitative Saccharomyces cerevisiae Proteome

Protein activity is the ultimate arbiter of function in most cellular pathways, and protein concentration is fundamentally connected to protein action. While the proteome of yeast has been subjected to the most comprehensive analysis of any eukaryote, existing datasets are difficult to compare, and there is no consensus abundance value for each protein. We evaluated 21 quantitative analyses of the S. cerevisiae proteome, normalizing and converting all measurements of protein abundance into the intuitive measurement of absolute molecules per cell. We estimate the cellular abundance of 92% of the proteins in the yeast proteome and assess the variation in each abundance measurement. Using our protein abundance dataset, we find that a global response to diverse environmental stresses is not detected at the level of protein abundance, we find that protein tags have only a modest effect on protein abundance, and we identify proteins that are differentially regulated at the mRNA abundance, mRNA translation, and protein abundance levels.

Single-cell RNA sequencing reveals intrinsic and extrinsic regulatory heterogeneity in yeast responding to stress

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells.

Genetically identical cells growing in the same environment can vary in their cellular state and behavior. Such heterogeneity may explain why some cells in an isogenic population can survive sudden severe environmental stress whereas other cells succumb. Cell-to-cell variation in gene expression has been linked to variable stress survival, but how and why transcript levels vary across the transcriptome in single cells is only beginning to emerge. Here, we used single-cell RNA sequencing (scRNA-seq) to measure cell-to-cell heterogeneity in the transcriptome of budding yeast (Saccharomyces cerevisiae). We find surprising patterns of variation across known sets of transcription factor targets, indicating that cells vary in their transcriptome profile both before and after stress exposure. scRNA-seq analysis combined with live-cell imaging of transcription factor activation dynamics revealed some cells in which the stress response was coordinately activated and other cells in which the traditional response was decoupled, suggesting unrecognized regulatory nuances that expand our understanding of stress response and survival.

Quantitating translational control: mRNA abundance-dependent and independent contributions and the mRNA sequences that specify them

Translation rate per mRNA molecule correlates positively with mRNA abundance. As a result, protein levels do not scale linearly with mRNA levels, but instead scale with the abundance of mRNA raised to the power of an ‘amplification exponent’. Here we show that to quantitate translational control, the translation rate must be decomposed into two components. One, TR_mD, depends on the mRNA level and defines the amplification exponent. The other, TR_mIND, is independent of mRNA amount and impacts the correlation coefficient between protein and mRNA levels. We show that in Saccharomyces cerevisiae TR_mD represents ∼20% of the variance in translation and directs an amplification exponent of 1.20 with a 95% confidence interval [1.14, 1.26]. TR_mIND constitutes the remaining ∼80% of the variance in translation and explains ∼5% of the variance in protein expression. We also find that TR_mD and TR_mIND are preferentially determined by different mRNA sequence features: TR_mIND by the length of the open reading frame and TR_mD both by a ∼60 nucleotide element that spans the initiating AUG and by codon and amino acid frequency. Our work provides more appropriate estimates of translational control and implies that TR_mIND is under different evolutionary selective pressures than TR_mD.

Genome-wide host RNA signatures of infectious diseases: discovery and clinical translation

The use of whole blood gene expression to derive diagnostic biomarkers capable of distinguishing between phenotypically similar diseases holds great promise but remains a challenge. Differential gene expression analysis is used to identify the key genes that undergo changes in expression relative to healthy individuals, as well as to patients with other diseases. These key genes can act as diagnostic, prognostic and predictive markers of disease. Gene expression ‘signatures’ in the blood hold the potential to be used for the diagnosis of infectious diseases, where current diagnostics are unreliable, ineffective or of limited potential. For diagnostic tests based on RNA signatures to be useful clinically, the first step is to identify the minimum set of gene transcripts that accurately identify the disease in question. The second requirement is rapid and cost‐effective detection of the gene expression levels. Signatures have been described for a number of infectious diseases, but ‘clinic‐ready’ technologies for RNA detection from clinical samples are limited, though existing methods such as RT‐PCR are likely to be superseded by a number of emerging technologies, which may form the basis of the translation of gene expression signatures into routine diagnostic tests for a range of disease states.

Deep Sequencing Transcriptome Analysis of Murine Wound Healing: Effects of a Multicomponent, Multitarget Natural Product Therapy-Tr14

Wound healing involves an orchestrated response that engages multiple processes, such as hemostasis, cellular migration, extracellular matrix synthesis, and in particular, inflammation. Using a murine model of cutaneous wound repair, the transcriptome was mapped from 12 h to 8 days post-injury, and in response to a multicomponent, multi-target natural product, Tr14. Using single-molecule RNA sequencing (RNA-seq), there were clear temporal changes in known transcripts related to wound healing pathways, and additional novel transcripts of both coding and non-coding genes. Tr14 treatment modulated >100 transcripts related to key wound repair pathways, such as response to wounding, wound contraction, and cytokine response. The results provide the most precise and comprehensive characterization to date of the transcriptome's response to skin damage, repair, and multicomponent natural product therapy. By understanding the wound repair process, and the effects of natural products, it should be possible to intervene more effectively in diseases involving aberrant repair.

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.

A note on statistical repeatability and study design for high-throughput assays

Characterizing the technical precision of measurements is a necessary stage in the planning of experiments and in the formal sample size calculation for optimal design. Instruments that measure multiple analytes simultaneously, such as in high-throughput assays arising in biomedical research, pose particular challenges from a statistical perspective. The current most popular method for assessing precision of high-throughput assays is by scatterplotting data from technical replicates. Here, we question the statistical rationale of this approach from both an empirical and theoretical perspective, illustrating our discussion using four example data sets from different genomic platforms. We demonstrate that such scatterplots convey little statistical information of relevance and are potentially highly misleading. We present an alternative framework for assessing the precision of high-throughput assays and planning biomedical experiments. Our methods are based on repeatability-a long-established statistical quantity also known as the intraclass correlation coefficient. We provide guidance and software for estimation and visualization of repeatability of high-throughput assays, and for its incorporation into study design. © 2016 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.

Multiple Transcript Properties Related to Translation Affect mRNA Degradation Rates in Saccharomyces cerevisiae

Degradation of mRNA contributes to variation in transcript abundance. Studies of individual mRNAs have shown that both cis and trans factors affect mRNA degradation rates. However, the factors underlying transcriptome-wide variation in mRNA degradation rates are poorly understood. We investigated the contribution of different transcript properties to transcriptome-wide degradation rate variation in the budding yeast, Saccharomyces cerevisiae, using multiple regression analysis. We find that multiple transcript properties are significantly associated with variation in mRNA degradation rates, and that a model incorporating these properties explains ∼50% of the genome-wide variance. Predictors of mRNA degradation rates include transcript length, ribosome density, biased codon usage, and GC content of the third position in codons. To experimentally validate these factors, we studied individual transcripts expressed from identical promoters. We find that decreasing ribosome density by mutating the first translational start site of a transcript increases its degradation rate. Using coding sequence variants of green fluorescent protein (GFP) that differ only at synonymous sites, we show that increased GC content of the third position of codons results in decreased rates of mRNA degradation. Thus, in steady-state conditions, a large fraction of genome-wide variation in mRNA degradation rates is determined by inherent properties of transcripts, many of which are related to translation, rather than specific regulatory mechanisms.

Analysis of the Yeast Peptidome and Comparison with the Human Peptidome

Peptides function as signaling molecules in species as diverse as humans and yeast. Mass spectrometry-based peptidomics techniques provide a relatively unbiased method to assess the peptidome of biological samples. In the present study, we used a quantitative peptidomic technique to characterize the peptidome of the yeast Saccharomyces cerevisiae and compare it to the peptidomes of mammalian cell lines and tissues. Altogether, 297 yeast peptides derived from 75 proteins were identified. The yeast peptides are similar to those of the human peptidome in average size and amino acid composition. Inhibition of proteasome activity with either bortezomib or epoxomicin led to decreased levels of some yeast peptides, suggesting that these peptides are generated by the proteasome. Approximately 30% of the yeast peptides correspond to the N- or C-terminus of the protein; the human peptidome is also highly represented in N- or C-terminal protein fragments. Most yeast and humans peptides are derived from a subset of abundant proteins, many with functions involving cellular metabolism or protein synthesis and folding. Of the 75 yeast proteins that give rise to peptides, 24 have orthologs that give rise to human and/or mouse peptides and for some, the same region of the proteins are found in the human, mouse, and yeast peptidomes. Taken together, these results support the hypothesis that intracellular peptides may have specific and conserved biological functions.

Constitutive turnover of histone H2A.Z at yeast promoters requires the preinitiation complex

The assembly of the preinitiation complex (PIC) occurs upstream of the +1 nucleosome which, in yeast, obstructs the transcription start site and is frequently assembled with the histone variant H2A.Z. To understand the contribution of the transcription machinery in the disassembly of the +1 H2A.Z nucleosome, conditional mutants were used to block PIC assembly. A quantitative ChIP-seq approach, which allows detection of global occupancy change, was employed to measure H2A.Z occupancy. Blocking PIC assembly resulted in promoter-specific H2A.Z accumulation, indicating that the PIC is required to evict H2A.Z. By contrast, H2A.Z eviction was unaffected upon depletion of INO80, a remodeler previously reported to displace nucleosomal H2A.Z. Robust PIC-dependent H2A.Z eviction was observed at active and infrequently transcribed genes, indicating that constitutive H2A.Z turnover is a general phenomenon. Finally, sites with strong H2A.Z turnover precisely mark transcript starts, providing a new metric for identifying cryptic and alternative sites of initiation.

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

To fit the genetic information of an animal, yeast or other eukaryote into cells, DNA is tightly wound around proteins called histones to form repeating units known as nucleosomes. However, this tight winding prevents proteins from accessing the DNA, and so prevents gene transcription – the first stage of producing the molecules encoded by a gene. For transcription to take place, nucleosomes at DNA sequences called promoters must be reorganized and disassembled, thereby allowing proteins to bind to and engage these sequences and to turn nearby genes on.

H2A is a histone protein that is found in the majority of nucleosomes in yeast cells. A different form of this histone – called H2A.Z – is found in nucleosomes near the promoter of almost every gene. It is thought that nucleosomes that contain H2A.Z are recognized and disassembled as the gene turns on, but it is unclear how this happens.

To investigate how H2A.Z nucleosomes are disassembled, Tramantano et al. depleted yeast cells of various proteins thought to play a role in the disassembly process. This indicated that the proteins that transcribe genes play crucial roles in the process of disassembling the H2A.Z nucleosomes, because H2A.Z accumulated at promoters in cells that are depleted of these proteins.

Further investigation revealed that disassembled H2A.Z nucleosomes are reassembled with H2A histones, before being converted back to the H2A.Z form by an enzyme called SWR1. This turnover of H2A.Z was seen at active genes and those that are infrequently transcribed, suggesting that it is a general phenomenon.

Tramantano et al. also found that the turnover rate of H2A.Z can be used to accurately predict the sites in the DNA where transcription starts. This observation could therefore help to identify previously unknown transcription start sites. Future work could address further questions about how H2A.Z nucleosomes are disassembled. For example, what is the mechanical force that drives this process? And at what step of the transcription process does it occur?

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

Jellyfish Bioactive Compounds: Methods for Wet-Lab Work

The study of bioactive compounds from marine animals has provided, over time, an endless source of interesting molecules. Jellyfish are commonly targets of study due to their toxic proteins. However, there is a gap in reviewing successful wet-lab methods employed in these animals, which compromises the fast progress in the detection of related biomolecules. Here, we provide a compilation of the most effective wet-lab methodologies for jellyfish venom extraction prior to proteomic analysis-separation, identification and toxicity assays. This includes SDS-PAGE, 2DE, gel chromatography, HPLC, DEAE, LC-MS, MALDI, Western blot, hemolytic assay, antimicrobial assay and protease activity assay. For a more comprehensive approach, jellyfish toxicity studies should further consider transcriptome sequencing. We reviewed such methodologies and other genomic techniques used prior to the deep sequencing of transcripts, including RNA extraction, construction of cDNA libraries and RACE. Overall, we provide an overview of the most promising methods and their successful implementation for optimizing time and effort when studying jellyfish.

Introns regulate the production of ribosomal proteins by modulating splicing of duplicated ribosomal protein genes

Most budding yeast introns exist in the many duplicated ribosomal protein genes (RPGs) and it has been posited that they remain there to modulate the expression of RPGs and cell growth in response to stress. However, the mechanism by which introns regulate the expression of RPGs and their impact on the synthesis of ribosomal proteins remain unclear. In this study, we show that introns determine the ratio of ribosomal protein isoforms through asymmetric paralog-specific regulation of splicing. Exchanging the introns and 3′ untranslated regions of the duplicated RPS9 genes altered the splicing efficiency and changed the ratio of the ribosomal protein isoforms. Mutational analysis of the RPS9 genes indicated that splicing is regulated by variations in the intron structure and the 3′ untranslated region. Together these data suggest that preferential splicing of duplicated RPGs provides a means for adjusting the ratio of different ribosomal protein isoforms, while maintaining the overall expression level of each ribosomal protein.

Cross-stress resistance in Saccharomyces cerevisiae yeast--new insight into an old phenomenon

Acquired stress resistance is the result of mild stress causing the acquisition of resistance to severe stress of the same or a different type. The mechanism of "same-stress" resistance (resistance to a second, strong stress after mild primary stress of the same type) probably depends on the activation of defense and repair mechanisms specific for a particular type of stress, while cross-stress resistance (i.e., resistance to a second, strong stress after a different type of mild primary stress) is the effect of activation of both a specific and general stress response program, which in Saccharomyces cerevisiae yeast is known as the environmental stress response (ESR). Advancements in research techniques have made it possible to study the mechanism of cross-stress resistance at various levels of cellular organization: stress signal transduction pathways, regulation of gene expression, and transcription or translation processes. As a result of this type of research, views on the cross-stress protection mechanism have been reconsidered. It was originally thought that cross-stress resistance, irrespective of the nature of the two stresses, was determined by universal mechanisms, i.e., the same mechanisms within the general stress response. They are now believed to be more specific and strictly dependent on the features of the first stress.

In vivo single-particle imaging of nuclear mRNA export in budding yeast demonstrates an essential role for Mex67p

Many messenger RNA export proteins have been identified; yet the spatial and temporal activities of these proteins and how they determine directionality of messenger ribonucleoprotein (mRNP) complex export from the nucleus remain largely undefined. Here, the bacteriophage PP7 RNA-labeling system was used in Saccharomyces cerevisiae to follow single-particle mRNP export events with high spatial precision and temporal resolution. These data reveal that mRNP export, consisting of nuclear docking, transport, and cytoplasmic release from a nuclear pore complex (NPC), is fast (∼ 200 ms) and that upon arrival in the cytoplasm, mRNPs are frequently confined near the nuclear envelope. Mex67p functions as the principal mRNP export receptor in budding yeast. In a mex67-5 mutant, delayed cytoplasmic release from NPCs and retrograde transport of mRNPs was observed. This proves an essential role for Mex67p in cytoplasmic mRNP release and directionality of transport.

Pericentromeric satellite repeat expansions through RNA-derived DNA intermediates in cancer

Aberrant transcription of the pericentromeric human satellite II (HSATII) repeat is present in a wide variety of epithelial cancers. In deriving experimental systems to study its deregulation, we observed that HSATII expression is induced in colon cancer cells cultured as xenografts or under nonadherent conditions in vitro, but it is rapidly lost in standard 2D cultures. Unexpectedly, physiological induction of endogenous HSATII RNA, as well as introduction of synthetic HSATII transcripts, generated cDNA intermediates in the form of DNA/RNA hybrids. Single molecule sequencing of tumor xenografts showed that HSATII RNA-derived DNA (rdDNA) molecules are stably incorporated within pericentromeric loci. Suppression of RT activity using small molecule inhibitors reduced HSATII copy gain. Analysis of whole-genome sequencing data revealed that HSATII copy number gain is a common feature in primary human colon tumors and is associated with a lower overall survival. Together, our observations suggest that cancer-associated derepression of specific repetitive sequences can promote their RNA-driven genomic expansion, with potential implications on pericentromeric architecture.

Single-Molecule Sequencing Reveals Estrogen-Regulated Clinically Relevant lncRNAs in Breast Cancer

Estrogen receptor (ER)α-positive tumors are commonly treated with ERα antagonists or inhibitors of estrogen synthesis, but most tumors develop resistance, and we need to better understand the pathways that underlie the proliferative and tumorigenic role of this estrogen-activated transcription factor. We here present the first single-molecule sequencing of the estradiol-induced ERα transcriptome in the luminal A-type human breast cancer cell lines MCF7 and T47D. Sequencing libraries were prepared from the polyadenylated RNA fraction after 8 hours of estrogen or vehicle treatment. Single-molecule sequencing was carried out in biological and technical replicates and differentially expressed genes were defined and analyzed for enriched processes. Correlation analysis with clinical expression and survival were performed, and follow-up experiments carried out using time series, chromatin immunoprecipitation and quantitative real-time PCR. We uncovered that ERα in addition to regulating approximately 2000 protein-coding genes, also regulated up to 1000 long noncoding RNAs (lncRNAs). Most of these were up-regulated, and 178 lncRNAs were regulated in both cell lines. We demonstrate that Long Intergenic Non-protein Coding RNA 1016 (LINC01016) and LINC00160 are direct transcriptional targets of ERα, correlate with ERα expression in clinical samples, and show prognostic significance in relation to breast cancer survival. We show that silencing of LINC00160 results in reduced proliferation, demonstrating that lncRNA expression have functional consequences. Our findings suggest that ERα regulation of lncRNAs is clinically relevant and that their functions and potential use as biomarkers for endocrine response are important to explore.

Ribosomal proteins as novel players in tumorigenesis

Ribosome biogenesis is the most demanding energetic and metabolic expenditure of the cell. The nucleolus, a nuclear compartment, coordinates rRNA transcription, maturation, and assembly into ribosome subunits. The transcription process is highly coordinated with ribosome biogenesis. In this context, ribosomal proteins (RPs) play a crucial role. In the last decade, an increasing number of studies have associated RPs with extraribosomal functions related to proliferation. Importantly, the expression of RPs appears to be deregulated in several human disorders due, at least in part, to genetic mutations. Although the deregulation of RPs in human malignancies is commonly observed, a more complex mechanism is believed to be involved, favoring the tumorigenic process, its progression and metastasis. This review explores the roles of the most frequently mutated oncogenes and tumor suppressor genes in human cancer that modulate ribosome biogenesis, including their interaction with RPs. In this regard, we propose a new focus for novel therapies.

Accounting for experimental noise reveals that mRNA levels, amplified by post-transcriptional processes, largely determine steady-state protein levels in yeast

Cells respond to their environment by modulating protein levels through mRNA transcription and post-transcriptional control. Modest observed correlations between global steady-state mRNA and protein measurements have been interpreted as evidence that mRNA levels determine roughly 40% of the variation in protein levels, indicating dominant post-transcriptional effects. However, the techniques underlying these conclusions, such as correlation and regression, yield biased results when data are noisy, missing systematically, and collinear---properties of mRNA and protein measurements---which motivated us to revisit this subject. Noise-robust analyses of 24 studies of budding yeast reveal that mRNA levels explain more than 85% of the variation in steady-state protein levels. Protein levels are not proportional to mRNA levels, but rise much more rapidly. Regulation of translation suffices to explain this nonlinear effect, revealing post-transcriptional amplification of, rather than competition with, transcriptional signals. These results substantially revise widely credited models of protein-level regulation, and introduce multiple noise-aware approaches essential for proper analysis of many biological phenomena.

Accounting for experimental noise reveals that mRNA levels, amplified by post-transcriptional processes, largely determine steady-state protein levels in yeast

Cells respond to their environment by modulating protein levels through mRNA transcription and post-transcriptional control. Modest observed correlations between global steady-state mRNA and protein measurements have been interpreted as evidence that mRNA levels determine roughly 40% of the variation in protein levels, indicating dominant post-transcriptional effects. However, the techniques underlying these conclusions, such as correlation and regression, yield biased results when data are noisy, missing systematically, and collinear---properties of mRNA and protein measurements---which motivated us to revisit this subject. Noise-robust analyses of 24 studies of budding yeast reveal that mRNA levels explain more than 85% of the variation in steady-state protein levels. Protein levels are not proportional to mRNA levels, but rise much more rapidly. Regulation of translation suffices to explain this nonlinear effect, revealing post-transcriptional amplification of, rather than competition with, transcriptional signals. These results substantially revise widely credited models of protein-level regulation, and introduce multiple noise-aware approaches essential for proper analysis of many biological phenomena.

The WTX Tumor Suppressor Interacts with the Transcriptional Corepressor TRIM28

WTX encodes a tumor suppressor implicated in the pediatric kidney cancer Wilms tumor and in mesenchymal differentiation with potentially distinct functions in the cytoplasm, at the plasma membrane, and in the nucleus. Although modulating components of the WNT signaling pathway is a proposed function for cytoplasmic and membrane-bound WTX, its nuclear properties are not well understood. Here we report that the transcriptional corepressor TRIM28 is the major binding partner for nuclear WTX. WTX interacted with the coiled coil domain of TRIM28 required for its binding to Krüppel-associated box domains of transcription factors and for its chromatin recruitment through its own coiled coil and proline-rich domains. Knockdown of endogenous WTX reduced the recruitment of TRIM28 to a chromatinized reporter sequence and its ability to repress a target transcript. In mouse embryonic stem cells where TRIM28 plays a major role in repressing endogenous retroviruses and long interspersed elements, knockdown of either TRIM28 or WTX combined with single molecule RNA sequencing revealed a highly significant shared set of differentially regulated transcripts, including derepression of non-coding repetitive sequences and their neighboring protein encoding genes (p < 1e-20). In mesenchymal precursor cells, depletion of WTX and TRIM28 resulted in analogous β-catenin-independent defects in adipogenic and osteogenic differentiation, and knockdown of WTX reduced TRIM28 binding to Pparγ promoter. Together, the physical and functional interaction between WTX and TRIM28 suggests that the nuclear fraction of WTX plays a role in epigenetic silencing, an effect that may contribute to its function as a regulator of cellular differentiation and tumorigenesis.

Estimating a structured covariance matrix from multi-lab measurements in high-throughput biology

We consider the problem of quantifying the degree of coordination between transcription and translation, in yeast. Several studies have reported a surprising lack of coordination over the years, in organisms as different as yeast and human, using diverse technologies. However, a close look at this literature suggests that the lack of reported correlation may not reflect the biology of regulation. These reports do not control for between-study biases and structure in the measurement errors, ignore key aspects of how the data connect to the estimand, and systematically underestimate the correlation as a consequence. Here, we design a careful meta-analysis of 27 yeast data sets, supported by a multilevel model, full uncertainty quantification, a suite of sensitivity analyses and novel theory, to produce a more accurate estimate of the correlation between mRNA and protein levels-a proxy for coordination. From a statistical perspective, this problem motivates new theory on the impact of noise, model mis-specifications and non-ignorable missing data on estimates of the correlation between high dimensional responses. We find that the correlation between mRNA and protein levels is quite high under the studied conditions, in yeast, suggesting that post-transcriptional regulation plays a less prominent role than previously thought.

Htt CAG repeat expansion confers pleiotropic gains of mutant huntingtin function in chromatin regulation

The CAG repeat expansion in the Huntington's disease gene HTT extends a polyglutamine tract in mutant huntingtin that enhances its ability to facilitate polycomb repressive complex 2 (PRC2). To gain insight into this dominant gain of function, we mapped histone modifications genome-wide across an isogenic panel of mouse embryonic stem cell (ESC) and neuronal progenitor cell (NPC) lines, comparing the effects of Htt null and different size Htt CAG mutations. We found that Htt is required in ESC for the proper deposition of histone H3K27me3 at a subset of 'bivalent' loci but in NPC it is needed at 'bivalent' loci for both the proper maintenance and the appropriate removal of this mark. In contrast, Htt CAG size, though changing histone H3K27me3, is prominently associated with altered histone H3K4me3 at 'active' loci. The sets of ESC and NPC genes with altered histone marks delineated by the lack of huntingtin or the presence of mutant huntingtin, though distinct, are enriched in similar pathways with apoptosis specifically highlighted for the CAG mutation. Thus, the manner by which huntingtin function facilitates PRC2 may afford mutant huntingtin with multiple opportunities to impinge upon the broader machinery that orchestrates developmentally appropriate chromatin status.