Compensation for differences in gene copy number among yeast ribosomal proteins is encoded within their promoters

Dedicated chaperones coordinate co-translational regulation of ribosomal protein production with ribosome assembly to preserve proteostasis

The biogenesis of eukaryotic ribosomes involves the ordered assembly of around 80 ribosomal proteins. Supplying equimolar amounts of assembly-competent ribosomal proteins is complicated by their aggregation propensity and the spatial separation of their location of synthesis and pre-ribosome incorporation. Recent evidence has highlighted that dedicated chaperones protect individual, unassembled ribosomal proteins on their path to the pre-ribosomal assembly site. Here, we show that the co-translational recognition of Rpl3 and Rpl4 by their respective dedicated chaperone, Rrb1 or Acl4, reduces the degradation of the encoding RPL3 and RPL4 mRNAs in the yeast Saccharomyces cerevisiae. In both cases, negative regulation of mRNA levels occurs when the availability of the dedicated chaperone is limited and the nascent ribosomal protein is instead accessible to a regulatory machinery consisting of the nascent-polypeptide-associated complex and the Caf130-associated Ccr4-Not complex. Notably, deregulated expression of Rpl3 and Rpl4 leads to their massive aggregation and a perturbation of overall proteostasis in cells lacking the E3 ubiquitin ligase Tom1. Taken together, we have uncovered an unprecedented regulatory mechanism that adjusts the de novo synthesis of Rpl3 and Rpl4 to their actual consumption during ribosome assembly and, thereby, protects cells from the potentially detrimental effects of their surplus production.

Living cells are packed full of molecules known as proteins, which perform many vital tasks the cells need to survive and grow. Machines called ribosomes inside the cells use template molecules called messenger RNAs (or mRNAs for short) to produce proteins. The newly-made proteins then have to travel to a specific location in the cell to perform their tasks. Some newly-made proteins are prone to forming clumps, so cells have other proteins known as chaperones that ensure these clumps do not form.

The ribosomes themselves are made up of several proteins, some of which are also prone to clumping as they are being produced. To prevent this from happening, cells control how many ribosomal proteins they make, so there are just enough to form the ribosomes the cell needs at any given time. Previous studies found that, in yeast, two ribosomal proteins called Rpl3 and Rpl4 each have their own dedicated chaperone to prevent them from clumping. However, it remained unclear whether these chaperones are also involved in regulating the levels of Rpl3 and Rpl4.

To address this question, Pillet et al. studied both of these dedicated chaperones in yeast cells. The experiments showed that the chaperones bound to their target proteins (either units of Rpl3 or Rpl4) as they were being produced on the ribosomes. This protected the template mRNAs the ribosomes were using to produce these proteins from being destroyed, thus allowing further units of Rpl3 and Rpl4 to be produced. When enough Rpl3 and Rpl4 units were made, there were not enough of the chaperones to bind them all, leaving the mRNA templates unprotected. This led to the destruction of the mRNA templates, which decreased the numbers of Rpl3 and Rpl4 units being produced.

The work of Pillet et al. reveals a feedback mechanism that allows yeast to tightly control the levels of Rpl3 and Rpl4. In the future, these findings may help us understand diseases caused by defects in ribosomal proteins, such as Diamond-Blackfan anemia, and possibly also neurodegenerative diseases caused by clumps of proteins forming in cells. The next step will be to find out whether the mechanism uncovered by Pillet et al. also exists in human and other mammalian cells.

Yeast Crf1p: An activator in need is an activator indeed

Ribosome biogenesis is an energetically costly process, and tight regulation is required for stoichiometric balance between components. This requires coordination of RNA polymerases I, II, and III. Lack of nutrients or the presence of stress leads to downregulation of ribosome biogenesis, a process for which mechanistic target of rapamycin complex I (mTORC1) is key. mTORC1 activity is communicated by means of specific transcription factors, and in yeast, which is a primary model system in which transcriptional coordination has been delineated, transcription factors involved in regulation of ribosomal protein genes include Fhl1p and its cofactors, Ifh1p and Crf1p. Ifh1p is an activator, whereas Crf1p has been implicated in maintaining the repressed state upon mTORC1 inhibition. Computational analyses of evolutionary relationships have indicated that Ifh1p and Crf1p descend from a common ancestor. Here, we discuss recent evidence, which suggests that Crf1p also functions as an activator. We propose a model that consolidates available experimental evidence, which posits that Crf1p functions as an alternate activator to prevent the stronger activator Ifh1p from re-binding gene promoters upon mTORC1 inhibition. The correlation between retention of Crf1p in related yeast strains and duplication of ribosomal protein genes suggests that this backup activation may be important to ensure gene expression when Ifh1p is limiting. With ribosome biogenesis as a hallmark of cell growth, failure to control assembly of ribosomal components leads to several human pathologies. A comprehensive understanding of mechanisms underlying this process is therefore of the essence.

Mediator dynamics during heat shock in budding yeast

The Mediator complex is central to transcription by RNA polymerase II (Pol II) in eukaryotes. In budding yeast (Saccharomyces cerevisiae), Mediator is recruited by activators and associates with core promoter regions, where it facilitates preinitiation complex (PIC) assembly, only transiently before Pol II escape. Interruption of the transcription cycle by inactivation or depletion of Kin28 inhibits Pol II escape and stabilizes this association. However, Mediator occupancy and dynamics have not been examined on a genome-wide scale in yeast grown in nonstandard conditions. Here we investigate Mediator occupancy following heat shock or CdCl₂ exposure, with and without depletion of Kin28. We find that Pol II occupancy shows similar dependence on Mediator under normal and heat shock conditions. However, although Mediator association increases at many genes upon Kin28 depletion under standard growth conditions, little or no increase is observed at most genes upon heat shock, indicating a more stable association of Mediator after heat shock. Unexpectedly, Mediator remains associated upstream of the core promoter at genes repressed by heat shock or CdCl₂ exposure whether or not Kin28 is depleted, suggesting that Mediator is recruited by activators but is unable to engage PIC components at these repressed targets. This persistent association is strongest at promoters that bind the HMGB family member Hmo1, and is reduced but not eliminated in hmo1Δ yeast. Finally, we show a reduced dependence on PIC components for Mediator occupancy at promoters after heat shock, further supporting altered dynamics or stronger engagement with activators under these conditions.

Advances in systems biology modeling: 10 years of crowdsourcing DREAM challenges

Computational and mathematical models are key to obtain a system-level understanding of biological processes, but their limitations have to be clearly defined to allow their proper application and interpretation. Crowdsourced benchmarks in the form of challenges provide an unbiased assessment of methods, and for the past decade, the Dialogue for Reverse Engineering Assessment and Methods (DREAM) organized more than 15 systems biology challenges. From transcription factor binding to dynamical network models, from signaling networks to gene regulation, from whole-cell models to cell-lineage reconstruction, and from single-cell positioning in a tissue to drug combinations and cell survival, the breadth is broad. To celebrate the 5-year anniversary of Cell Systems, we review the genesis of these systems biology challenges and discuss how interlocking the forward- and reverse-modeling paradigms allows to push the rim of systems biology. This approach will persist for systems levels approaches in biology and medicine.

Mechanisms coordinating ribosomal protein gene transcription in response to stress

While expression of ribosomal protein genes (RPGs) in the budding yeast has been extensively studied, a longstanding enigma persists regarding their co-regulation under fluctuating growth conditions. Most RPG promoters display one of two distinct arrangements of a core set of transcription factors (TFs) and are further differentiated by the presence or absence of the HMGB protein Hmo1. However, a third group of promoters appears not to be bound by any of these proteins, raising the question of how the whole suite of genes is co-regulated. We demonstrate here that all RPGs are regulated by two distinct, but complementary mechanisms driven by the TFs Ifh1 and Sfp1, both of which are required for maximal expression in optimal conditions and coordinated downregulation upon stress. At the majority of RPG promoters, Ifh1-dependent regulation predominates, whereas Sfp1 plays the major role at all other genes. We also uncovered an unexpected protein homeostasis-dependent binding property of Hmo1 at RPG promoters. Finally, we show that the Ifh1 paralog Crf1, previously described as a transcriptional repressor, can act as a constitutive RPG activator. Our study provides a more complete picture of RPG regulation and may serve as a paradigm for unravelling RPG regulation in multicellular eukaryotes.

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

Modulating transcription through development of semi-synthetic yeast core promoters

Altering gene expression regulation by promoter engineering is a very effective way to fine-tune heterologous pathways in eukaryotic hosts. Typically, pathway building approaches in yeast still use a limited set of long, native promoters. With the today’s introduction of longer and more complex pathways, an expansion of this synthetic biology toolbox is necessary. In this study we elucidated the core promoter structure of the well-characterized yeast TEF1 promoter and determined the minimal length needed for sufficient protein expression. Furthermore, this minimal core promoter sequence was used for the creation of a promoter library covering different expression strengths. This resulted in a group of short, 69 bp promoters with an 8.0-fold expression range. One exemplar had a two and four times higher expression compared to the native CYC1 and ADH1 promoter, respectively. Additionally, as it was described that the protein expression range could be broadened by upstream activating sequences (UASs), we integrated earlier described single and multiple short, synthetic UASs in front of the strongest yeast core promoter. This approach resulted to further variation in protein expression and an overall promoter library spanning a 20-fold activity range and covering a length from 69 bp to maximally 129 bp. Furthermore, the robustness of this library was assessed on three alternative carbon sources besides glucose. As such, the suitability of short yeast core promoters for metabolic engineering applications on different media, either in an individual context or combined with UAS elements, was demonstrated.

Synthetic 5' UTRs Can Either Up- or Downregulate Expression upon RNA-Binding Protein Binding

The construction of complex gene-regulatory networks requires both inhibitory and upregulatory modules. However, the vast majority of RNA-based regulatory "parts" are inhibitory. Using a synthetic biology approach combined with SHAPE-seq, we explored the regulatory effect of RNA-binding protein (RBP)-RNA interactions in bacterial 5' UTRs. By positioning a library of RNA hairpins upstream of a reporter gene and co-expressing them with the matching RBP, we observed a set of regulatory responses, including translational stimulation, translational repression, and cooperative behavior. Our combined approach revealed three distinct states in vivo: in the absence of RBPs, the RNA molecules can be found in either a molten state that is amenable to translation or a structured phase that inhibits translation. In the presence of RBPs, the RNA molecules are in a semi-structured phase with partial translational capacity. Our work provides new insight into RBP-based regulation and a blueprint for designing complete gene-regulatory circuits at the post-transcriptional level.

An in Vivo Binding Assay for RNA-Binding Proteins Based on Repression of a Reporter Gene

We study translation repression in bacteria by engineering a regulatory circuit that functions as a binding assay for RNA binding proteins (RBP) in vivo. We do so by inducing expression of a fluorescent protein-RBP chimera, together with encoding its binding site at various positions within the ribosomal initiation region (+11-13 nt from the AUG) of a reporter module. We show that when bound by their cognate RBPs, the phage coat proteins for PP7 (PCP) and Qβ (QCP), strong repression is observed for all hairpin positions within the initiation region. Yet, a sharp transition to no-effect is observed when positioned in the elongation region, at a single-nucleotide resolution. Employing in vivo Selective 2'-hydroxyl acylation analyzed by primer extension followed by sequencing (SHAPE-seq) for a representative construct, established that in the translationally active state the mRNA molecule is nonstructured, while in the repressed state a structured signature was detected. We then utilize this regulatory phenomena to quantify the binding affinity of the coat proteins of phages MS2, PP7, GA, and Qβ to 14 cognate and noncognate binding sites in vivo. Using our circuit, we demonstrate qualitative differences between in vitro to in vivo binding characteristics for various variants when comparing to past studies. Furthermore, by introducing a simple mutation to the loop region for the Qβ-wt site, MCP binding is abolished, creating the first high-affinity QCP site that is completely orthogonal to MCP. Consequently, we demonstrate that our hybrid transcriptional-post-transcriptional circuit can be utilized as a binding assay to quantify RNA-RBP interactions in vivo.

The Genetic Makeup and Expression of the Glycolytic and Fermentative Pathways Are Highly Conserved Within the Saccharomyces Genus

The ability of the yeast Saccharomyces cerevisiae to convert glucose, even in the presence of oxygen, via glycolysis and the fermentative pathway to ethanol has played an important role in its domestication. Despite the extensive knowledge on these pathways in S. cerevisiae, relatively little is known about their genetic makeup in other industrially relevant Saccharomyces yeast species. In this study we explore the diversity of the glycolytic and fermentative pathways within the Saccharomyces genus using S. cerevisiae, Saccharomyces kudriavzevii, and Saccharomyces eubayanus as paradigms. Sequencing data revealed a highly conserved genetic makeup of the glycolytic and fermentative pathways in the three species in terms of number of paralogous genes. Although promoter regions were less conserved between the three species as compared to coding sequences, binding sites for Rap1, Gcr1 and Abf1, main transcriptional regulators of glycolytic and fermentative genes, were highly conserved. Transcriptome profiling of these three strains grown in aerobic batch cultivation in chemically defined medium with glucose as carbon source, revealed a remarkably similar expression of the glycolytic and fermentative genes across species, and the conserved classification of genes into major and minor paralogs. Furthermore, transplantation of the promoters of major paralogs of S. kudriavzevii and S. eubayanus into S. cerevisiae demonstrated not only the transferability of these promoters, but also the similarity of their strength and response to various environmental stimuli. The relatively low homology of S. kudriavzevii and S. eubayanus promoters to their S. cerevisiae relatives makes them very attractive alternatives for strain construction in S. cerevisiae, thereby expanding the S. cerevisiae molecular toolbox.

Intrinsic cooperativity potentiates parallel cis-regulatory evolution

Convergent evolutionary events in independent lineages provide an opportunity to understand why evolution favors certain outcomes over others. We studied such a case where a large set of genes-those coding for the ribosomal proteins-gained cis-regulatory sequences for a particular transcription regulator (Mcm1) in independent fungal lineages. We present evidence that these gains occurred because Mcm1 shares a mechanism of transcriptional activation with an ancestral regulator of the ribosomal protein genes, Rap1. Specifically, we show that Mcm1 and Rap1 have the inherent ability to cooperatively activate transcription through contacts with the general transcription factor TFIID. Because the two regulatory proteins share a common interaction partner, the presence of one ancestral cis-regulatory sequence can 'channel' random mutations into functional sites for the second regulator. At a genomic scale, this type of intrinsic cooperativity can account for a pattern of parallel evolution involving the fixation of hundreds of substitutions.

Codon usage of highly expressed genes affects proteome-wide translation efficiency

Although the genetic code is redundant, synonymous codons for the same amino acid are not used with equal frequencies in genomes, a phenomenon termed "codon usage bias." Previous studies have demonstrated that synonymous changes in a coding sequence can exert significant cis effects on the gene's expression level. However, whether the codon composition of a gene can also affect the translation efficiency of other genes has not been thoroughly explored. To study how codon usage bias influences the cellular economy of translation, we massively converted abundant codons to their rare synonymous counterpart in several highly expressed genes in Escherichia coli This perturbation reduces both the cellular fitness and the translation efficiency of genes that have high initiation rates and are naturally enriched with the manipulated codon, in agreement with theoretical predictions. Interestingly, we could alleviate the observed phenotypes by increasing the supply of the tRNA for the highly demanded codon, thus demonstrating that the codon usage of highly expressed genes was selected in evolution to maintain the efficiency of global protein translation.

Inferring time derivatives including cell growth rates using Gaussian processes

Often the time derivative of a measured variable is of as much interest as the variable itself. For a growing population of biological cells, for example, the population's growth rate is typically more important than its size. Here we introduce a non-parametric method to infer first and second time derivatives as a function of time from time-series data. Our approach is based on Gaussian processes and applies to a wide range of data. In tests, the method is at least as accurate as others, but has several advantages: it estimates errors both in the inference and in any summary statistics, such as lag times, and allows interpolation with the corresponding error estimation. As illustrations, we infer growth rates of microbial cells, the rate of assembly of an amyloid fibril and both the speed and acceleration of two separating spindle pole bodies. Our algorithm should thus be broadly applicable.

Patterns of ribosomal protein expression specify normal and malignant human cells

Background

Ribosomes are highly conserved molecular machines whose core composition has traditionally been regarded as invariant. However, recent studies have reported intriguing differences in the expression of some ribosomal proteins (RPs) across tissues and highly specific effects on the translation of individual mRNAs.

Results

To determine whether RPs are more generally linked to cell identity, we analyze the heterogeneity of RP expression in a large set of human tissues, primary cells, and tumors. We find that about a quarter of human RPs exhibit tissue-specific expression and that primary hematopoietic cells display the most complex patterns of RP expression, likely shaped by context-restricted transcriptional regulators. Strikingly, we uncover patterns of dysregulated expression of individual RPs across cancer types that arise through copy number variations and are predictive for disease progression.

Conclusions

Our study reveals an unanticipated plasticity of RP expression across normal and malignant human cell types and provides a foundation for future characterization of cellular behaviors that are orchestrated by specific RPs.

Electronic supplementary material

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

Massively Parallel Interrogation of the Effects of Gene Expression Levels on Fitness

Data of gene expression levels across individuals, cell types, and disease states is expanding, yet our understanding of how expression levels impact phenotype is limited. Here, we present a massively parallel system for assaying the effect of gene expression levels on fitness in Saccharomyces cerevisiae by systematically altering the expression level of ∼100 genes at ∼100 distinct levels spanning a 500-fold range at high resolution. We show that the relationship between expression levels and growth is gene and environment specific and provides information on the function, stoichiometry, and interactions of genes. Wild-type expression levels in some conditions are not optimal for growth, and genes whose fitness is greatly affected by small changes in expression level tend to exhibit lower cell-to-cell variability in expression. Our study addresses a fundamental gap in understanding the functional significance of gene expression regulation and offers a framework for evaluating the phenotypic effects of expression variation.

Screening of an Escherichia coli promoter library for a phenylalanine biosensor

In recent years, the application of transcription factor-based biosensors for the engineering of microbial production strains opened up new opportunities for industrial biotechnology. However, the design of synthetic regulatory circuits depends on the selection of suitable transcription factor-promoter pairs to convert the concentration of effector molecules into a measureable output. Here, we present an efficient strategy to screen promoter libraries for appropriate parts for biosensor design. To this end, we pooled the strains of the Alon library containing about 2000 different Escherichia coli promoter-gfpmut2 fusions, and enriched galactose- and L-phenylalanine-responsive promoters by toggled rounds of positive and negative selection using fluorescence-activated cell sorting (FACS). For both effectors, responsive promoters were isolated and verified by cultivation in microtiter plates. The promoter of mtr, encoding an L-tryptophan-specific transporter, was identified as suitable part for the construction of an L-phenylalanine biosensor. In the following, we performed a comparative analysis of different biosensor constructs based on the mtr promoter. The obtained data revealed a strong influence of the biosensor architecture on the performance characteristics. For proof-of-principle, the mtr sensor was applied in a FACS high-throughput screening of an E. coli MG1655 mutant library for the isolation of L-phenylalanine producers. These results emphasize the developed screening approach as a convenient strategy for the identification of effector-responsive promoters for the design of novel biosensors.

A Modeling Framework for Generation of Positional and Temporal Simulations of Transcriptional Regulation

We present a modeling framework aimed at capturing both the positional and temporal behavior of transcriptional regulatory proteins in eukaryotic cells. There is growing evidence that transcriptional regulation is the complex behavior that emerges not solely from the individual components, but rather from their collective behavior, including competition and cooperation. Our framework describes individual regulatory components using generic action oriented descriptions of their biochemical interactions with a DNA sequence. All the possible actions are based on the current state of factors bound to the DNA. We developed a rule builder to automatically generate the complete set of biochemical interaction rules for any given DNA sequence. Off-the-shelf stochastic simulation engines can model the behavior of a system of rules and the resulting changes in the configuration of bound factors can be visualized. We compared our model to experimental data at well-studied loci in yeast, confirming that our model captures both the positional and temporal behavior of transcriptional regulation.

Prediction of fine-tuned promoter activity from DNA sequence

The quantitative prediction of transcriptional activity of genes using promoter sequence is fundamental to the engineering of biological systems for industrial purposes and understanding the natural variation in gene expression. To catalyze the development of new algorithms for this purpose, the Dialogue on Reverse Engineering Assessment and Methods (DREAM) organized a community challenge seeking predictive models of promoter activity given normalized promoter activity data for 90 ribosomal protein promoters driving expression of a fluorescent reporter gene. By developing an unbiased modeling approach that performs an iterative search for predictive DNA sequence features using the frequencies of various k-mers, inferred DNA mechanical properties and spatial positions of promoter sequences, we achieved the best performer status in this challenge. The specific predictive features used in the model included the frequency of the nucleotide G, the length of polymeric tracts of T and TA, the frequencies of 6 distinct trinucleotides and 12 tetranucleotides, and the predicted protein deformability of the DNA sequence. Our method accurately predicted the activity of 20 natural variants of ribosomal protein promoters (Spearman correlation r = 0.73) as compared to 33 laboratory-mutated variants of the promoters (r = 0.57) in a test set that was hidden from participants. Notably, our model differed substantially from the rest in 2 main ways: i) it did not explicitly utilize transcription factor binding information implying that subtle DNA sequence features are highly associated with gene expression, and ii) it was entirely based on features extracted exclusively from the 100 bp region upstream from the translational start site demonstrating that this region encodes much of the overall promoter activity. The findings from this study have important implications for the engineering of predictable gene expression systems and the evolution of gene expression in naturally occurring biological systems.

Molecular mechanisms of ribosomal protein gene coregulation

The 137 ribosomal protein genes (RPG) of Saccharomyces provide a model for gene coregulation. Reja et al. examine the positional and functional organization of their regulators (Rap1, Fhl1, Ifh1, Sfp1, and Hmo1), the transcription machinery (TFIIB, TFIID, and RNA polymerase II), and chromatin at near-base-pair resolution using ChIP-exo.

The 137 ribosomal protein genes (RPGs) of Saccharomyces provide a model for gene coregulation. We examined the positional and functional organization of their regulators (Rap1 [repressor activator protein 1], Fhl1, Ifh1, Sfp1, and Hmo1), the transcription machinery (TFIIB, TFIID, and RNA polymerase II), and chromatin at near-base-pair resolution using ChIP-exo, as RPGs are coordinately reprogrammed. Where Hmo1 is enriched, Fhl1, Ifh1, Sfp1, and Hmo1 cross-linked broadly to promoter DNA in an RPG-specific manner and demarcated by general minor groove widening. Importantly, Hmo1 extended 20–50 base pairs (bp) downstream from Fhl1. Upon RPG repression, Fhl1 remained in place. Hmo1 dissociated, which was coupled to an upstream shift of the +1 nucleosome, as reflected by the Hmo1 extension and core promoter region. Fhl1 and Hmo1 may create two regulatable and positionally distinct barriers, against which chromatin remodelers position the +1 nucleosome into either an activating or a repressive state. Consistent with in vitro studies, we found that specific TFIID subunits, in addition to cross-linking at the core promoter, made precise cross-links at Rap1 sites, which we interpret to reflect native Rap1–TFIID interactions. Our findings suggest how sequence-specific DNA binding regulates nucleosome positioning and transcription complex assembly >300 bp away and how coregulation coevolved with coding sequences.

Nucleolar stress with and without p53

A veritable explosion of primary research papers within the past 10 years focuses on nucleolar and ribosomal stress, and for good reason: with ribosome biosynthesis consuming ~80% of a cell’s energy, nearly all metabolic and signaling pathways lead ultimately to or from the nucleolus. We begin by describing p53 activation upon nucleolar stress resulting in cell cycle arrest or apoptosis. The significance of this mechanism cannot be understated, as oncologists are now inducing nucleolar stress strategically in cancer cells as a potential anti-cancer therapy. We also summarize the human ribosomopathies, syndromes in which ribosome biogenesis or function are impaired leading to birth defects or bone narrow failures; the perplexing problem in the ribosomopathies is why only certain cells are affected despite the fact that the causative mutation is systemic. We then describe p53-independent nucleolar stress, first in yeast which lacks p53, and then in other model metazoans that lack MDM2, the critical E3 ubiquitin ligase that normally inactivates p53. Do these presumably ancient p53-independent nucleolar stress pathways remain latent in human cells? If they still exist, can we use them to target >50% of known human cancers that lack functional p53?

A Random Screen Using a Novel Reporter Assay System Reveals a Set of Sequences That Are Preferred as the TATA or TATA-Like Elements in the CYC1 Promoter of Saccharomyces cerevisiae

In Saccharomyces cerevisiae, the core promoters of class II genes contain either TATA or TATA-like elements to direct accurate transcriptional initiation. Genome-wide analyses show that the consensus sequence of the TATA element is TATAWAWR (8 bp), whereas TATA-like elements carry one or two mismatches to this consensus. The fact that several functionally distinct TATA sequences have been identified indicates that these elements may function, at least to some extent, in a gene-specific manner. The purpose of the present study was to identify functional TATA sequences enriched in one particular core promoter and compare them with the TATA or TATA-like elements that serve as the pre-initiation complex (PIC) assembly sites on the yeast genome. For this purpose, we conducted a randomized screen of the TATA element in the CYC1 promoter by using a novel reporter assay system and identified several hundreds of unique sequences that were tentatively classified into nine groups. The results indicated that the 7 bp TATA element (i.e., TATAWAD) and several sets of TATA-like sequences are preferred specifically by this promoter. Furthermore, we find that the most frequently isolated TATA-like sequence, i.e., TATTTAAA, is actually utilized as a functional core promoter element for the endogenous genes, e.g., ADE5,7 and ADE6. Collectively, these results indicate that the sequence requirements for the functional TATA or TATA-like elements in one particular core promoter are not as stringent. However, the variation of these sequences differs significantly from that of the PIC assembly sites on the genome, presumably depending on promoter structures and reflecting the gene-specific function of these sequences.

Control of relative timing and stoichiometry by a master regulator

Developmental processes in cells require a series of complex steps. Often only a single master regulator activates genes in these different steps. This poses several challenges: some targets need to be ordered temporally, while co-functional targets may need to be synchronized in both time and expression level. Here we study in single cells the dynamic activation patterns of early meiosis genes in budding yeast, targets of the meiosis master regulator Ime1. We quantify the individual roles of the promoter and protein levels in expression pattern control, as well as the roles of individual promoter elements. We find a consistent expression pattern difference between a non-cofunctional pair of genes, and a highly synchronized activation of a co-functional pair. We show that dynamic control leading to these patterns is distributed between promoter, gene and external regions. Through specific reciprocal changes to the promoters of pairs of genes, we show that different genes can use different promoter elements to reach near identical activation patterns.

Systematic dissection of the sequence determinants of gene 3' end mediated expression control

The 3'end genomic region encodes a wide range of regulatory process including mRNA stability, 3' end processing and translation. Here, we systematically investigate the sequence determinants of 3' end mediated expression control by measuring the effect of 13,000 designed 3' end sequence variants on constitutive expression levels in yeast. By including a high resolution scanning mutagenesis of more than 200 native 3' end sequences in this designed set, we found that most mutations had only a mild effect on expression, and that the vast majority (~90%) of strongly effecting mutations localized to a single positive TA-rich element, similar to a previously described 3' end processing efficiency element, and resulted in up to ten-fold decrease in expression. Measurements of 3' UTR lengths revealed that these mutations result in mRNAs with aberrantly long 3'UTRs, confirming the role for this element in 3' end processing. Interestingly, we found that other sequence elements that were previously described in the literature to be part of the polyadenylation signal had a minor effect on expression. We further characterize the sequence specificities of the TA-rich element using additional synthetic 3' end sequences and show that its activity is sensitive to single base pair mutations and strongly depends on the A/T content of the surrounding sequences. Finally, using a computational model, we show that the strength of this element in native 3' end sequences can explain some of their measured expression variability (R = 0.41). Together, our results emphasize the importance of efficient 3' end processing for endogenous protein levels and contribute to an improved understanding of the sequence elements involved in this process.

Gcn1 contacts the small ribosomal protein Rps10, which is required for full activation of the protein kinase Gcn2

General control non-derepressible 2 (Gcn2) is involved in many biological functions. We found that Gcn1 directly contacts small ribosomal protein (Rps)10 on the ribosome and this is required for maintaining basal Gcn2-activity and efficiently activating Gcn2 under amino acid starvation.

Molecular dissection of the genetic mechanisms that underlie expression conservation in orthologous yeast ribosomal promoters

Recent studies have shown a surprising phenomenon, whereby orthologous regulatory regions from different species drive similar expression levels despite being highly diverged in sequence. Here, we investigated this phenomenon by genomically integrating hundreds of ribosomal protein (RP) promoters from nine different yeast species into S. cerevisiae and accurately measuring their activity. We found that orthologous RP promoters have extreme expression conservation even across evolutionarily distinct yeast species. Notably, our measurements reveal two distinct mechanisms that underlie this conservation and which act in different regions of the promoter. In the core promoter region, we found compensatory changes, whereby effects of sequence variations in one part of the core promoter were reversed by variations in another part. In contrast, we observed robustness in Rap1 transcription factor binding sites, whereby significant sequence variations had little effect on promoter activity. Finally, cases in which orthologous promoter activities were not conserved could largely be explained by the sequence variation within the core promoter. Together, our results provide novel insights into the mechanisms by which expression is conserved throughout evolution across diverged promoter sequences.

Transcriptional response to the [ISP(+) ] prion of Saccharomyces cerevisiae differs from that induced by the deletion of its structural gene, SFP1

Currently, several protein-based genetic determinants, or prions, are described in yeast, and several hundred prion candidates have been predicted. Importantly, many known and potential prion proteins regulate transcription; therefore, prion induction should affect gene expression. While it is generally believed that the prion phenotype should mimic the deletion phenotype, this rule has exceptions. Formed by the transcription factor Sfp1p, [ISP(+) ] is one such exception as the [ISP(+) ] and sfp1Δ strains differ in many phenotypic traits. These data suggest that effects of prion formation by a transcription factor and its absence may affect gene expression in a different way. However, studies addressing this issue are practically absent. Here, we explore how [ISP(+) ] affects gene expression and how these changes correspond to the effect of SFP1 deletion. Our data indicate that the [ISP(+) ]-related expression changes cannot be explained by the inactivation of Sfp1p. Remarkably, most Sfp1p targets are not affected in the [ISP(+) ] strain; instead, the genes upregulated in the [ISP(+) ] strain are enriched in Gcn4p and Aft1p targets. We propose that Sfp1p serves as a part of a regulatory complex, and the activity of this complex may be modulated differently by the absence or prionization of Sfp1p.

Two distinct promoter architectures centered on dynamic nucleosomes control ribosomal protein gene transcription

In yeast, ribosome production is controlled transcriptionally by tight coregulation of the 138 ribosomal protein genes (RPGs). RPG promoters display limited sequence homology, and the molecular basis for their coregulation remains largely unknown. Here we identify two prevalent RPG promoter types, both characterized by upstream binding of the general transcription factor (TF) Rap1 followed by the RPG-specific Fhl1/Ifh1 pair, with one type also binding the HMG-B protein Hmo1. We show that the regulatory properties of the two promoter types are remarkably similar, suggesting that they are determined to a large extent by Rap1 and the Fhl1/Ifh1 pair. Rapid depletion experiments allowed us to define a hierarchy of TF binding in which Rap1 acts as a pioneer factor required for binding of all other TFs. We also uncovered unexpected features underlying recruitment of Fhl1, whose forkhead DNA-binding domain is not required for binding at most promoters, and Hmo1, whose binding is supported by repeated motifs. Finally, we describe unusually micrococcal nuclease (MNase)-sensitive nucleosomes at all RPG promoters, located between the canonical +1 and -1 nucleosomes, which coincide with sites of Fhl1/Ifh1 and Hmo1 binding. We speculate that these "fragile" nucleosomes play an important role in regulating RPG transcriptional output.

Promoters maintain their relative activity levels under different growth conditions

Most genes change expression levels across conditions, but it is unclear which of these changes represents specific regulation and what determines their quantitative degree. Here, we accurately measured activities of ~900 S. cerevisiae and ~1800 E. coli promoters using fluorescent reporters. We show that in both organisms 60-90% of promoters change their expression between conditions by a constant global scaling factor that depends only on the conditions and not on the promoter's identity. Quantifying such global effects allows precise characterization of specific regulation-promoters deviating from the global scale line. These are organized into few functionally related groups that also adhere to scale lines and preserve their relative activities across conditions. Thus, only several scaling factors suffice to accurately describe genome-wide expression profiles across conditions. We present a parameter-free passive resource allocation model that quantitatively accounts for the global scaling factors. It suggests that many changes in expression across conditions result from global effects and not specific regulation, and provides means for quantitative interpretation of expression profiles.

Mapping the fine structure of a eukaryotic promoter input-output function

The precise tuning of gene expression levels is essential for the optimal performance of transcriptional regulatory networks. We created 209 variants of the Saccharomyces cerevisiae PHO5 promoter to quantify how different binding sites for the transcription factor Pho4 affect its output. We found that transcription-factor binding affinities determined in vitro could quantitatively predict the output of a complex yeast promoter. Promoter output was precisely tunable by subtle changes in binding-site affinity of less than 3 kcal mol(-1), which are accessible by modifying 1-2 bases. Our results provide insights into how transcription-factor binding sites regulate gene expression, their possible evolution and how they can be used to precisely tune gene expression. More generally, we show that in vitro binding-energy landscapes of transcription factors can precisely predict the output of a native yeast promoter, indicating that quantitative models of transcriptional regulatory networks are feasible.

Inferring gene expression from ribosomal promoter sequences, a crowdsourcing approach

The Gene Promoter Expression Prediction challenge consisted of predicting gene expression from promoter sequences in a previously unknown experimentally generated data set. The challenge was presented to the community in the framework of the sixth Dialogue for Reverse Engineering Assessments and Methods (DREAM6), a community effort to evaluate the status of systems biology modeling methodologies. Nucleotide-specific promoter activity was obtained by measuring fluorescence from promoter sequences fused upstream of a gene for yellow fluorescence protein and inserted in the same genomic site of yeast Saccharomyces cerevisiae. Twenty-one teams submitted results predicting the expression levels of 53 different promoters from yeast ribosomal protein genes. Analysis of participant predictions shows that accurate values for low-expressed and mutated promoters were difficult to obtain, although in the latter case, only when the mutation induced a large change in promoter activity compared to the wild-type sequence. As in previous DREAM challenges, we found that aggregation of participant predictions provided robust results, but did not fare better than the three best algorithms. Finally, this study not only provides a benchmark for the assessment of methods predicting activity of a specific set of promoters from their sequence, but it also shows that the top performing algorithm, which used machine-learning approaches, can be improved by the addition of biological features such as transcription factor binding sites.

Deciphering the rules by which 5'-UTR sequences affect protein expression in yeast

The 5'-untranslated region (5'-UTR) of mRNAs contains elements that affect expression, yet the rules by which these regions exert their effect are poorly understood. Here, we studied the impact of 5'-UTR sequences on protein levels in yeast, by constructing a large-scale library of mutants that differ only in the 10 bp preceding the translational start site of a fluorescent reporter. Using a high-throughput sequencing strategy, we obtained highly accurate measurements of protein abundance for over 2,000 unique sequence variants. The resulting pool spanned an approximately sevenfold range of protein levels, demonstrating the powerful consequences of sequence manipulations of even 1-10 nucleotides immediately upstream of the start codon. We devised computational models that predicted over 70% of the measured expression variability in held-out sequence variants. Notably, a combined model of the most prominent features successfully explained protein abundance in an additional, independently constructed library, whose nucleotide composition differed greatly from the library used to parameterize the model. Our analysis reveals the dominant contribution of the start codon context at positions -3 to -1, mRNA secondary structure, and out-of-frame upstream AUGs (uAUGs) to phenotypic diversity, thereby advancing our understanding of how protein levels are modulated by 5'-UTR sequences, and paving the way toward predictably tuning protein expression through manipulations of 5'-UTRs.

Determinants of nucleosome positioning

Nucleosome positioning is critical for gene expression and most DNA-related processes. Here we review the dominant patterns of nucleosome positioning that have been observed and summarize the current understanding of their underlying determinants. The genome-wide pattern of nucleosome positioning is determined by the combination of DNA sequence, ATP-dependent nucleosome remodeling enzymes and transcription factors that include activators, components of the preinitiation complex and elongating RNA polymerase II. These determinants influence each other such that the resulting nucleosome positioning patterns are likely to differ among genes and among cells in a population, with consequent effects on gene expression.

Promoter sequence determines the relationship between expression level and noise

A single transcription factor can activate or repress expression by three different mechanisms: one that increases cell-to-cell variability in target gene expression (noise) and two that decrease noise.

The ability of cells to accurately control gene expression levels in response to extracellular cues is limited by the inherently stochastic nature of transcriptional regulation. A change in transcription factor (TF) activity results in changes in the expression of its targets, but the way in which cell-to-cell variability in expression (noise) changes as a function of TF activity, and whether targets of the same TF behave similarly, is not known. Here, we measure expression and noise as a function of TF activity for 16 native targets of the transcription factor Zap1 that are regulated by it through diverse mechanisms. For most activated and repressed Zap1 targets, noise decreases as expression increases. Kinetic modeling suggests that this is due to two distinct Zap1-mediated mechanisms that both change the frequency of transcriptional bursts. Notably, we found that another mechanism of repression by Zap1, which is encoded in the promoter DNA, likely decreases the size of transcriptional bursts, producing a unique transcriptional state characterized by low expression and low noise. In addition, we find that further reduction in noise is achieved when a single TF both activates and represses a single target gene. Our results suggest a global principle whereby at low TF concentrations, the dominant source of differences in expression between promoters stems from differences in burst frequency, whereas at high TF concentrations differences in burst size dominate. Taken together, we show that the precise amount by which noise changes with expression is specific to the regulatory mechanism of transcription and translation that acts at each gene.

In response to environmental changes, cells regulate the activity of transcription factors (TFs), which in turn change the expression of dozens of downstream target genes by binding to their promoters. The response of each target gene is determined by the interplay between TF concentration and the context in which TF binding sites occur in each target promoter. To examine the relationship between promoter sequence, mechanism of regulation, and response to TF activity, we measured expression of 16 target genes of a single TF in response to changes in TF concentration in single cells. We found that different native promoters that are all targets of the same TF exhibit diverse responses to changing TF levels in terms of both gene expression level and cell-to-cell variability (noise) in expression. Using computational modeling and mutations of specific promoter elements, we show that the molecular mechanisms of regulation can be inferred by measuring how noise changes with expression. These results show that a single TF can regulate transcription through multiple mechanisms, resulting in similar changes in mean expression but vastly different changes in cell-to-cell variability.

Measurements of the impact of 3' end sequences on gene expression reveal wide range and sequence dependent effects

A full understanding of gene regulation requires an understanding of the contributions that the various regulatory regions have on gene expression. Although it is well established that sequences downstream of the main promoter can affect expression, our understanding of the scale of this effect and how it is encoded in the DNA is limited. Here, to measure the effect of native S. cerevisiae 3′ end sequences on expression, we constructed a library of 85 fluorescent reporter strains that differ only in their 3′ end region. Notably, despite being driven by the same strong promoter, our library spans a continuous twelve-fold range of expression values. These measurements correlate with endogenous mRNA levels, suggesting that the 3′ end contributes to constitutive differences in mRNA levels. We used deep sequencing to map the 3′UTR ends of our strains and show that determination of polyadenylation sites is intrinsic to the local 3′ end sequence. Polyadenylation mapping was followed by sequence analysis, we found that increased A/T content upstream of the main polyadenylation site correlates with higher expression, both in the library and genome-wide, suggesting that native genes differ by the encoded efficiency of 3′ end processing. Finally, we use single cells fluorescence measurements, in different promoter activation levels, to show that 3′ end sequences modulate protein expression dynamics differently than promoters, by predominantly affecting the size of protein production bursts as opposed to the frequency at which these bursts occur. Altogether, our results lead to a more complete understanding of gene regulation by demonstrating that 3′ end regions have a unique and sequence dependent effect on gene expression.

A basic question in gene expression is the relative contribution of different regulatory layers and genomic regions to the differences in protein levels. In this work we concentrated on the effect of 3′ end sequences. For this, we constructed a library of yeast strains that differ only by a native 3′ end region integrated downstream to a reported gene driven by a constant inducible promoter. Thus we could attribute all differences in reporter expression between the strains to the different 3′ end sequences. Interestingly, we found that despite being driven by the same strong, inducible promoter, our library spanned a wide and continuous range of expression levels of more than twelve-fold. As these measurements represent the sole effect of the 3′ end region, we quantify the contribution of these sequences to the variance in mRNA levels by comparing our measurements to endogenous mRNA levels. We follow by sequence analysis to find a simple sequence signature that correlates with expression. In addition, single cell analysis reveals distinct noise dynamics of 3′ end mediated differences in expression compared to different levels of promoter activation leading to a more complete understanding of gene expression which also incorporates the effect of these regions.

Two DNA-encoded strategies for increasing expression with opposing effects on promoter dynamics and transcriptional noise

Individual cells from a genetically identical population exhibit substantial variation in gene expression. A significant part of this variation is due to noise in the process of transcription that is intrinsic to each gene, and is determined by factors such as the rate with which the promoter transitions between transcriptionally active and inactive states, and the number of transcripts produced during the active state. However, we have a limited understanding of how the DNA sequence affects such promoter dynamics. Here, we used single-cell time-lapse microscopy to compare the effect on transcriptional dynamics of two distinct types of sequence changes in the promoter that can each increase the mean expression of a cell population by similar amounts but through different mechanisms. We show that increasing expression by strengthening a transcription factor binding site results in slower promoter dynamics and higher noise as compared with increasing expression by adding nucleosome-disfavoring sequences. Our results suggest that when achieving the same mean expression, the strategy of using stronger binding sites results in a larger number of transcripts produced from the active state, whereas the strategy of adding nucleosome-disfavoring sequences results in a higher frequency of promoter transitions between active and inactive states. In the latter strategy, this increased sampling of the active state likely reduces the expression variability of the cell population. Our study thus demonstrates the effect of cis-regulatory elements on expression variability and points to concrete types of sequence changes that may allow partial decoupling of expression level and noise.

On the expansion of "dangerous" gene repertoires by whole-genome duplications in early vertebrates

The emergence and evolutionary expansion of gene families implicated in cancers and other severe genetic diseases is an evolutionary oddity from a natural selection perspective. Here, we show that gene families prone to deleterious mutations in the human genome have been preferentially expanded by the retention of "ohnolog" genes from two rounds of whole-genome duplication (WGD) dating back from the onset of jawed vertebrates. We further demonstrate that the retention of many ohnologs suspected to be dosage balanced is in fact indirectly mediated by their susceptibility to deleterious mutations. This enhanced retention of "dangerous" ohnologs, defined as prone to autosomal-dominant deleterious mutations, is shown to be a consequence of WGD-induced speciation and the ensuing purifying selection in post-WGD species. These findings highlight the importance of WGD-induced nonadaptive selection for the emergence of vertebrate complexity, while rationalizing, from an evolutionary perspective, the expansion of gene families frequently implicated in genetic disorders and cancers.

Rapid synthesis of defined eukaryotic promoter libraries

Current gene synthesis methods allow the generation of long segments of dsDNA. We show that these techniques can be used to create synthetic regulatory elements and describe a method for the creation of completely defined, synthetic variants of the PHO5 promoter from the budding yeast Saccharomyces cerevisae. Overall, 128 promoters were assembled by high-temperature ligation, cloned into plasmids by isothermal assembly, maintained in E. coli, and consequently transformed into yeast by homologous recombination. Synthesis errors occurred at frequencies comparable to or lower than those achieved with current gene synthesis methods. The promoter synthesis method reported here is robust, fast, and readily accessible. Synthetically engineered promoter libraries will be useful tools for dissecting the intricacies of promoter input-output functions and may serve as tunable components for synthetic genetic networks.

Regulatory elements of Caenorhabditis elegans ribosomal protein genes

Background

Ribosomal protein genes (RPGs) are essential, tightly regulated, and highly expressed during embryonic development and cell growth. Even though their protein sequences are strongly conserved, their mechanism of regulation is not conserved across yeast, Drosophila, and vertebrates. A recent investigation of genomic sequences conserved across both nematode species and associated with different gene groups indicated the existence of several elements in the upstream regions of C. elegans RPGs, providing a new insight regarding the regulation of these genes in C. elegans.

Results

In this study, we performed an in-depth examination of C. elegans RPG regulation and found nine highly conserved motifs in the upstream regions of C. elegans RPGs using the motif discovery algorithm DME. Four motifs were partially similar to transcription factor binding sites from C. elegans, Drosophila, yeast, and human. One pair of these motifs was found to co-occur in the upstream regions of 250 transcripts including 22 RPGs. The distance between the two motifs displayed a complex frequency pattern that was related to their relative orientation.

We tested the impact of three of these motifs on the expression of rpl-2 using a series of reporter gene constructs and showed that all three motifs are necessary to maintain the high natural expression level of this gene. One of the motifs was similar to the binding site of an orthologue of POP-1, and we showed that RNAi knockdown of pop-1 impacts the expression of rpl-2. We further determined the transcription start site of rpl-2 by 5’ RACE and found that the motifs lie 40–90 bases upstream of the start site. We also found evidence that a noncoding RNA, contained within the outron of rpl-2, is co-transcribed with rpl-2 and cleaved during trans-splicing.

Conclusions

Our results indicate that C. elegans RPGs are regulated by a complex novel series of regulatory elements that is evolutionarily distinct from those of all other species examined up until now.

Large ribosomal protein 4 increases efficiency of viral recoding sequences

Expression of retroviral replication enzymes (Pol) requires a controlled translational recoding event to bypass the stop codon at the end of gag. This recoding event occurs either by direct suppression of termination via the insertion of an amino acid at the stop codon (readthrough) or by alteration of the mRNA reading frame (frameshift). Here we report the effects of a host protein, large ribosomal protein 4 (RPL4), on the efficiency of recoding. Using a dual luciferase reporter assay, we found that transfection of cells with a plasmid encoding RPL4 cDNA increases recoding efficiency in a dose-dependent manner, with a maximal enhancement of nearly twofold. Expression of RPL4 increases recoding of reporters containing retroviral readthrough and frameshift sequences, as well as the Sindbis virus leaky termination signal. RPL4-induced enhancement of recoding is cell line specific and appears to be specific to RPL4 among ribosomal proteins. Cotransfection of RPL4 cDNA with Moloney murine leukemia proviral DNA results in Gag processing defects and a reduction of viral particle formation, presumably caused by the RPL4-dependent alteration of the Gag-to-Gag-Pol ratio required for virion assembly and release.

Manipulating nucleosome disfavoring sequences allows fine-tune regulation of gene expression in yeast

Understanding how precise control of gene expression is specified within regulatory DNA sequences is a key challenge with far-reaching implications. Many studies have focused on the regulatory role of transcription factor-binding sites. Here, we explore the transcriptional effects of different elements, nucleosome-disfavoring sequences and, specifically, poly(dA:dT) tracts that are highly prevalent in eukaryotic promoters. By measuring promoter activity for a large-scale promoter library, designed with systematic manipulations to the properties and spatial arrangement of poly(dA:dT) tracts, we show that these tracts significantly and causally affect transcription. We show that manipulating these elements offers a general genetic mechanism, applicable to promoters regulated by different transcription factors, for tuning expression in a predictable manner, with resolution that can be even finer than that attained by altering transcription factor sites. Overall, our results advance the understanding of the regulatory code and suggest a potential mechanism by which promoters yielding prespecified expression patterns can be designed.

Inferring gene regulatory logic from high-throughput measurements of thousands of systematically designed promoters

Despite much research, our understanding of the rules by which cis-regulatory sequences are translated into expression levels is still lacking. We devised a method for obtaining parallel and highly accurate expression measurements of thousands of fully designed promoters, and applied it to measure the effect of systematic changes to location, number, orientation, affinity and organization of transcription factor (TF) binding sites and of nucleosome disfavoring sequences. Our analyses reveal a clear relationship between expression and binding site number, and TF-specific dependencies of expression on the distance between sites and gene starts including a striking ~10bp periodic relationship. We also demonstrate the utility of our approach for measuring TF sequence specificities and sensitivity of TF sites to surrounding sequence context, and for profiling the activity of most yeast transcription factors. Our method is readily applicable for studying both the cis and trans effects of genotype on transcriptional, post-transcriptional, and translational control.

Strategies to discover regulatory circuits of the mammalian immune system

Recent advances in technologies for genome- and proteome-scale measurements and perturbations promise to accelerate discovery in every aspect of biology and medicine. Although such rapid technological progress provides a tremendous opportunity, it also demands that we learn how to use these tools effectively. One application with great potential to enhance our understanding of biological systems is the unbiased reconstruction of genetic and molecular networks. Cells of the immune system provide a particularly useful model for developing and applying such approaches. Here, we review approaches for the reconstruction of signalling and transcriptional networks, with a focus on applications in the mammalian innate immune system.