Msn2 coordinates a stoichiometric gene expression program

N88S seipin-related seipinopathy is a lipidopathy associated with loss of iron homeostasis

BACKGROUND

Seipin is a protein encoded by the BSCL2 gene in humans and SEI1 gene in yeast, forming an Endoplasmic Reticulum (ER)-bound homo-oligomer. This oligomer is crucial in targeting ER-lipid droplet (LD) contact sites, facilitating the delivery of triacylglycerol (TG) to nascent LDs. Mutations in BSCL2, particularly N88S and S90L, lead to seipinopathies, which correspond to a cohort of motor neuron diseases (MNDs) characterized by the accumulation of misfolded N88S seipin into inclusion bodies (IBs) and cellular dysfunctions.

METHODS

Quantitative untargeted mass spectrometric proteomic and lipidomic analyses were conducted to examine changes in protein and lipid abundance in wild-type (WT) versus N88S seipin-expressing mutant cells. Differentially expressed proteins were categorized into functional networks to highlight altered protein functions and signaling pathways. Statistical comparisons were made using unpaired Student's t-tests or two-way ANOVA followed by Tukey´s / Šídák's multiple comparisons tests. P-values < 0.05 are considered significant.

RESULTS

In a well-established yeast model of N88S seipinopathy, misfolded N88S seipin forms IBs and exhibits higher levels of ER stress, leading to decreased cell viability due to increased reactive oxygen species (ROS), oxidative damage, lipid peroxidation, and reduced antioxidant activity. Proteomic and lipidomic analyses revealed alterations in phosphatidic acid (PA) levels, associated with disrupted inositol metabolism and decreased flux towards phospholipid biosynthesis. Importantly, deregulation of lipid metabolism contributed to ER stress beyond N88S seipin misfolding and IB formation. Additionally, the model exhibited deregulated iron (Fe) homeostasis during lifespan. N88S seipin-expressing cells showed impaired ability to cope with iron deficiency. This was linked to changes in the expression of Aft1p-controlled iron regulon genes, including the mRNA-binding protein CTH2 and the high-affinity iron transport system member FET3, in a p38/Hog1p- and Msn2p/Msn4p-dependent manner. Importantly, we unraveled a novel link between inositol metabolism and activation of the iron regulon in cells expressing the N88S seipin mutation. Despite iron accumulation, this was not associated with oxidative stress.

CONCLUSIONS

The study highlights that the effects of N88S seipin mutation extend beyond protein misfolding, with significant disruptions in lipid metabolism and iron homeostasis. This research marks a substantial advance in understanding and defining the roles of proteins and signaling pathways that contribute to human seipinopathy. Altered cellular processes, as well as potential therapeutic targets and biomarkers, were identified and can be explored in translational studies using human cell models.

Regulated resource reallocation is transcriptionally hard wired into the yeast stress response

Many organisms maintain generalized stress responses activated by adverse conditions. Although details vary, a common theme is the redirection of transcriptional and translational capacity away from growth-promoting genes and toward defense genes. Yet the precise roles of these coupled programs are difficult to dissect. Here we investigated Saccharomyces cerevisiae responding to salt as a model stressor. We used molecular, genomic, and single-cell microfluidic methods to examine the interplay between transcription factors Msn2 and Msn4 that induce stress-defense genes and Dot6 and Tod6 that transiently repress growth-promoting genes during stress. Surprisingly, loss of Dot6/Tod6 led to slower acclimation to salt, whereas loss of Msn2/4 produced faster growth during stress. This supports a model where transient repression of growth-promoting genes accelerates the Msn2/4 response, which is essential for acquisition of subsequent peroxide tolerance. Remarkably, we find that Msn2/4 regulate DOT6 mRNA production, influence Dot6 activation dynamics, and are required for full repression of growth-promoting genes. Thus, Msn2/4 directly regulate resource reallocation needed to mount their own response. We discuss broader implications for common stress responses across organisms.

The acyl-CoA-binding protein VdAcb1 is essential for carbon starvation response and contributes to virulence in Verticillium dahliae

In the face of carbon, nitrogen, and phosphorus starvation, microorganisms have evolved adaptive mechanisms to maintain growth. In a previous study, we identified a protein predicted to contain acyl-CoA-binding domains in the plant pathogenic fungus Verticillium dahliae. The predicted protein, designated VdAcb1, possesses an atypical signal peptide. However, the functions of this acyl-CoA-binding protein in V. dahliae are not clear. In this research, in vivo or in vitro assays confirmed that VdAcb1 is secreted extracellularly from V. dahliae, although it does not have the typical signal peptide. Furthermore, the unconventional secretion of VdAcb1 was dependent on VdGRASP, a member of the compartment for unconventional protein secretion (CUPS). The deletion mutant strain of VdAcb1 (ΔVdAcb1) exhibited significant sensitivity to carbon starvation. RNA-seq revealed that the expression of genes related to filamentous growth (MSB2 pathway) and sugar transport were regulated by VdAcb1 under conditions of carbon starvation. Yeast one-hybrid experiments further showed that the expression of VdAcb1 was positively regulated by the transcription factor VdMsn4. The ΔVdAcb1 strain showed significantly reduced virulence on Gossypium hirsutum and Nicotiana benthamiana. We hypothesize that under conditions of carbon starvation, the expression of VdAcb1 is activated by VdMsn4 and VdAcb1 is secreted into the extracellular space. In turn, this activates the downstream MAPK pathway to enhance filamentous growth and virulence of V. dahliae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-024-00175-3.

Optimization of the Irf8 +32-kb enhancer disrupts dendritic cell lineage segregation

Autoactivation of lineage-determining transcription factors mediates bistable expression, generating distinct cell phenotypes essential for complex body plans. Classical type 1 dendritic cell (cDC1) and type 2 dendritic cell (cDC2) subsets provide nonredundant functions for defense against distinct immune challenges. Interferon regulatory factor 8 (IRF8), the cDC1 lineage-determining transcription factor, undergoes autoactivation in cDC1 progenitors to establish cDC1 identity, yet its expression is downregulated during cDC2 differentiation by an unknown mechanism. This study reveals that the Irf8 +32-kb enhancer, responsible for IRF8 autoactivation, is naturally suboptimized with low-affinity IRF8 binding sites. Introducing multiple high-affinity IRF8 sites into the Irf8 +32-kb enhancer causes a gain-of-function effect, leading to erroneous IRF8 autoactivation in specified cDC2 progenitors, redirecting them toward cDC1 and a novel hybrid DC subset with mixed-lineage phenotypes. Further, this also causes a loss-of-function effect, reducing Irf8 expression in cDC1s. These developmental alterations critically impair both cDC1-dependent and cDC2-dependent arms of immunity. Collectively, our findings underscore the significance of enhancer suboptimization in the developmental segregation of cDCs required for normal immune function.

Anatomy of a superenhancer

Interferon regulatory factor-8 (IRF8) is the lineage determining transcription factor for the type one classical dendritic cell (cDC1) subset, a terminal selector for plasmacytoid dendritic cells and important for the function of monocytes. Studies of Irf8 gene regulation have identified several enhancers controlling its activity during development of progenitors in the bone marrow that precisely regulate expression at distinct developmental stages. Each enhancer responds to distinct transcription factors that are expressed at each stage. IRF8 is first expressed in early progenitors that form the monocyte dendritic cell progenitor (MDP) in response to induction of the transcription factor CCAAT/enhancer-binding protein alpha (C/EBPα) acting at the Irf8 +56 kb enhancer. IRF8 levels increase further as the MDP transits into the common dendritic cell progenitor (CDP) in response to E protein activity at the Irf8 +41 kb enhancer. Upon Nfil3-induction in CDPs leading to specification of the cDC1 progenitor, abrupt induction of BATF3 forms the JUN/BATF3/IRF8 heterotrimer that activates the Irf8 +32 kb enhancer that sustains Irf8 autoactivation throughout the cDC1 lifetime. Deletions of each of these enhancers has revealed their stage dependent activation. Surprisingly, studies of compound heterozygotes for each combination of enhancer deletions revealed that activation of each subsequent enhancer requires the successful activation of the previous enhancer in strictly cis-dependent mechanism. Successful progression of enhancer activation is finely tuned to alter the functional accessibility of subsequent enhancers to factors active in the next stage of development. The molecular basis for these phenomenon is still obscure but could have implications for genomic regulation in a broader developmental context.

Mutations of ribosomal protein genes induce overexpression of catalase in Saccharomyces cerevisiae

Ribosome assembly defects result in ribosomopathies, primarily caused by inadequate protein synthesis and induced oxidative stress. This study aimed to investigate the link between deleting one ribosomal protein gene (RPG) paralog and oxidative stress response. Our results indicated that RPG mutants exhibited higher oxidant sensitivity than the wild type (WT). The concentrations of H2O2 were increased in the RPG mutants. Catalase and superoxide dismutase (SOD) activities were generally higher at the stationary phase, with catalase showing particularly elevated activity in the RPG mutants. While both catalase genes, CTT1 and CTA1, consistently exhibited higher transcription in RPG mutants, Ctt1 primarily contributed to the increased catalase activity. Stress-response transcription factors Msn2, Msn4, and Hog1 played a role in regulating these processes. Previous studies have demonstrated that H2O2 can cleave 25S rRNA via the Fenton reaction, enhancing ribosomes' ability to translate mRNAs associated with oxidative stress-related genes. The cleavage of 25S rRNA was consistently more pronounced, and the translation efficiency of CTT1 and CTA1 mRNAs was altered in RPG mutants. Our results provide evidence that the mutations in RPGs increase H2O2 levels in vivo and elevate catalase expression through both transcriptional and translational controls.

Integrating analog and digital modes of gene expression at Arabidopsis FLC

Quantitative gene regulation at the cell population level can be achieved by two fundamentally different modes of regulation at individual gene copies. A 'digital' mode involves binary ON/OFF expression states, with population-level variation arising from the proportion of gene copies in each state, while an 'analog' mode involves graded expression levels at each gene copy. At the Arabidopsis floral repressor FLOWERING LOCUS C (FLC), 'digital' Polycomb silencing is known to facilitate quantitative epigenetic memory in response to cold. However, whether FLC regulation before cold involves analog or digital modes is unknown. Using quantitative fluorescent imaging of FLC mRNA and protein, together with mathematical modeling, we find that FLC expression before cold is regulated by both analog and digital modes. We observe a temporal separation between the two modes, with analog preceding digital. The analog mode can maintain intermediate expression levels at individual FLC gene copies, before subsequent digital silencing, consistent with the copies switching OFF stochastically and heritably without cold. This switch leads to a slow reduction in FLC expression at the cell population level. These data present a new paradigm for gradual repression, elucidating how analog transcriptional and digital epigenetic memory pathways can be integrated.

Molecular encoding of stimulus features in a single sensory neuron type enables neuronal and behavioral plasticity

Neurons modify their transcriptomes in response to an animal's experience. How specific experiences are transduced to modulate gene expression and precisely tune neuronal functions are not fully defined. Here, we describe the molecular profile of a thermosensory neuron pair in C. elegans experiencing different temperature stimuli. We find that distinct salient features of the temperature stimulus, including its duration, magnitude of change, and absolute value, are encoded in the gene expression program in this single neuron type, and we identify a novel transmembrane protein and a transcription factor whose specific transcriptional dynamics are essential to drive neuronal, behavioral, and developmental plasticity. Expression changes are driven by broadly expressed activity-dependent transcription factors and corresponding cis-regulatory elements that nevertheless direct neuron- and stimulus-specific gene expression programs. Our results indicate that coupling of defined stimulus characteristics to the gene regulatory logic in individual specialized neuron types can customize neuronal properties to drive precise behavioral adaptation.

Transcription factor localization dynamics and DNA binding drive distinct promoter interpretations

Environmental information may be encoded in the temporal dynamics of transcription factor (TF) activation and subsequently decoded by gene promoters to enact stimulus-specific gene expression programs. Previous studies of this behavior focused on the encoding and decoding of information in TF nuclear localization dynamics, yet cells control the activity of TFs in myriad ways, including by regulating their ability to bind DNA. Here, we use light-controlled mutants of the yeast TF Msn2 as a model system to investigate how promoter decoding of TF localization dynamics is affected by changes in the ability of the TF to bind DNA. We find that yeast promoters directly decode the light-controlled localization dynamics of Msn2 and that the effects of changing Msn2 affinity on that decoding behavior are highly promoter dependent, illustrating how cells could regulate TF localization dynamics and DNA binding in concert for improved control of gene expression.

Molecular encoding of stimulus features in a single sensory neuron type enables neuronal and behavioral plasticity

Neurons modify their transcriptomes in response to an animal’s experience. How specific experiences are transduced to modulate gene expression and precisely tune neuronal functions are not fully defined. Here, we describe the molecular profile of a thermosensory neuron pair in C. elegans experiencing different temperature stimuli. We find that distinct salient features of the temperature stimulus including its duration, magnitude of change, and absolute value are encoded in the gene expression program in this single neuron, and identify a novel transmembrane protein and a transcription factor whose specific transcriptional dynamics are essential to drive neuronal, behavioral, and developmental plasticity. Expression changes are driven by broadly expressed activity-dependent transcription factors and corresponding cis -regulatory elements that nevertheless direct neuron- and stimulus-specific gene expression programs. Our results indicate that coupling of defined stimulus characteristics to the gene regulatory logic in individual specialized neuron types can customize neuronal properties to drive precise behavioral adaptation.

Response and regulatory mechanisms of heat resistance in pathogenic fungi

Abstract

Both the increasing environmental temperature in nature and the defensive body temperature response to pathogenic fungi during mammalian infection cause heat stress during the fungal existence, reproduction, and pathogenic infection. To adapt and respond to the changing environment, fungi initiate a series of actions through a perfect thermal response system, conservative signaling pathways, corresponding transcriptional regulatory system, corresponding physiological and biochemical processes, and phenotypic changes. However, until now, accurate response and regulatory mechanisms have remained a challenge. Additionally, at present, the latest research progress on the heat resistance mechanism of pathogenic fungi has not been summarized. In this review, recent research investigating temperature sensing, transcriptional regulation, and physiological, biochemical, and morphological responses of fungi in response to heat stress is discussed. Moreover, the specificity thermal adaptation mechanism of pathogenic fungi in vivo is highlighted. These data will provide valuable knowledge to further understand the fungal heat adaptation and response mechanism, especially in pathogenic heat-resistant fungi.

Key points

• Mechanisms of fungal perception of heat pressure are reviewed.

• The regulatory mechanism of fungal resistance to heat stress is discussed.

• The thermal adaptation mechanism of pathogenic fungi in the human body is highlighted.

Signal amplification in the KEAP1-NRF2-ARE antioxidant response pathway

The KEAP1-NRF2-ARE signaling pathway plays a central role in mediating the adaptive cellular stress response to oxidative and electrophilic chemicals. This canonical pathway has been extensively studied and reviewed in the past two decades, but rarely was it looked at from a quantitative signaling perspective. Signal amplification, i.e., ultrasensitivity, is crucially important for robust induction of antioxidant genes to appropriate levels that can adequately counteract the stresses. In this review article, we examined a number of well-known molecular events in the KEAP1-NRF2-ARE pathway from a quantitative perspective with a focus on how signal amplification can be achieved. We illustrated, by using a series of mathematical models, that redox-regulated protein sequestration, stabilization, translation, nuclear trafficking, DNA promoter binding, and transcriptional induction - which are embedded in the molecular network comprising KEAP1, NRF2, sMaf, p62, and BACH1 - may generate highly ultrasensitive NRF2 activation and antioxidant gene induction. The emergence and degree of ultrasensitivity depend on the strengths of protein-protein and protein-DNA interaction and protein abundances. A unique, quantitative understanding of signal amplification in the KEAP1-NRF2-ARE pathway will help to identify sensitive targets for the prevention and therapeutics of oxidative stress-related diseases and develop quantitative adverse outcome pathway models to facilitate the health risk assessment of oxidative chemicals.

Primordial super-enhancers: heat shock-induced chromatin organization in yeast

Specialized mechanisms ensure proper expression of critically important genes such as those specifying cell identity or conferring protection from environmental stress. Investigations of the heat shock response have been critical in elucidating basic concepts of transcriptional control. Recent studies demonstrate that in response to thermal stress, heat shock-responsive genes associate with high levels of transcriptional activators and coactivators and those in yeast intensely interact across and between chromosomes, coalescing into condensates. In mammalian cells, cell identity genes that are regulated by super-enhancers (SEs) are also densely occupied by transcriptional machinery that form phase-separated condensates. We suggest that the stress-remodeled yeast nucleome bears functional and structural resemblance to mammalian SEs, and will reveal fundamental mechanisms of gene control by transcriptional condensates.

Paralogs in the PKA Regulon Traveled Different Evolutionary Routes to Divergent Expression in Budding Yeast

Functional divergence of duplicate genes, or paralogs, is an important driver of novelty in evolution. In the model yeast Saccharomyces cerevisiae, there are 547 paralog gene pairs that survive from an interspecies Whole Genome Hybridization (WGH) that occurred ~100MYA. In this work, we report that ~1/6th (110) of these WGH paralogs pairs (or ohnologs) are differentially expressed with a striking pattern upon Protein Kinase A (PKA) inhibition. One member of each pair in this group has low basal expression that increases upon PKA inhibition, while the other has moderate and unchanging expression. For these genes, expression of orthologs upon PKA inhibition in the non-WGH species Kluyveromyces lactis and for PKA-related stresses in other budding yeasts shows unchanging expression, suggesting that lack of responsiveness to PKA was likely the typical ancestral phenotype prior to duplication. Promoter sequence analysis across related budding yeast species further revealed that the subsequent emergence of PKA-dependence took different evolutionary routes. In some examples, regulation by PKA and differential expression appears to have arisen following the WGH, while in others, regulation by PKA appears to have arisen in one of the two parental lineages prior to the WGH. More broadly, our results illustrate the unique opportunities presented by a WGH event for generating functional divergence by bringing together two parental lineages with separately evolved regulation into one species. We propose that functional divergence of two ohnologs can be facilitated through such regulatory divergence.

An Improved CRISPR Interference Tool to Engineer Rhodococcus opacus

Rhodococcus opacus is a nonmodel bacterium that is well suited for valorizing lignin. Despite recent advances in our systems-level understanding of its versatile metabolism, studies of its gene functions at a single gene level are still lagging. Elucidating gene functions in nonmodel organisms is challenging due to limited genetic engineering tools that are convenient to use. To address this issue, we developed a simple gene repression system based on CRISPR interference (CRISPRi). This gene repression system uses a T7 RNA polymerase system to express a small guide RNA, demonstrating improved repression compared to the previously demonstrated CRISPRi system (i.e., the maximum repression efficiency improved from 58% to 85%). Additionally, our cloning strategy allows for building multiple CRISPRi plasmids in parallel without any PCR step, facilitating the engineering of this GC-rich organism. Using the improved CRISPRi system, we confirmed the annotated roles of four metabolic pathway genes, which had been identified by our previous transcriptomic analysis to be related to the consumption of benzoate, vanillate, catechol, and acetate. Furthermore, we showed our tool's utility by demonstrating the inducible accumulation of muconate that is a precursor of adipic acid, an important monomer for nylon production. While the maximum muconate yield obtained using our tool was 30% of the yield obtained using gene knockout, our tool showed its inducibility and partial repressibility. Our CRISPRi tool will be useful to facilitate functional studies of this nonmodel organism and engineer this promising microbial chassis for lignin valorization.

Of numbers and movement - understanding transcription factor pathogenesis by advanced microscopy

Transcription factors (TFs) are life-sustaining and, therefore, the subject of intensive research. By regulating gene expression, TFs control a plethora of developmental and physiological processes, and their abnormal function commonly leads to various developmental defects and diseases in humans. Normal TF function often depends on gene dosage, which can be altered by copy-number variation or loss-of-function mutations. This explains why TF haploinsufficiency (HI) can lead to disease. Since aberrant TF numbers frequently result in pathogenic abnormalities of gene expression, quantitative analyses of TFs are a priority in the field. In vitro single-molecule methodologies have significantly aided the identification of links between TF gene dosage and transcriptional outcomes. Additionally, advances in quantitative microscopy have contributed mechanistic insights into normal and aberrant TF function. However, to understand TF biology, TF-chromatin interactions must be characterised in vivo, in a tissue-specific manner and in the context of both normal and altered TF numbers. Here, we summarise the advanced microscopy methodologies most frequently used to link TF abundance to function and dissect the molecular mechanisms underlying TF HIs. Increased application of advanced single-molecule and super-resolution microscopy modalities will improve our understanding of how TF HIs drive disease.

A new protocol for single-cell RNA-seq reveals stochastic gene expression during lag phase in budding yeast

Current methods for single-cell RNA sequencing (scRNA-seq) of yeast cells do not match the throughput and relative simplicity of the state-of-the-art techniques that are available for mammalian cells. In this study, we report how 10x Genomics' droplet-based single-cell RNA sequencing technology can be modified to allow analysis of yeast cells. The protocol, which is based on in-droplet spheroplasting of the cells, yields an order-of-magnitude higher throughput in comparison to existing methods. After extensive validation of the method, we demonstrate its use by studying the dynamics of the response of isogenic yeast populations to a shift in carbon source, revealing the heterogeneity and underlying molecular processes during this shift. The method we describe opens new avenues for studies focusing on yeast cells, as well as other cells with a degradable cell wall.

Two roles for the yeast transcription coactivator SAGA and a set of genes redundantly regulated by TFIID and SAGA

Deletions within genes coding for subunits of the transcription coactivator SAGA caused strong genome-wide defects in transcription and SAGA-mediated chromatin modifications. In contrast, rapid SAGA depletion produced only modest transcription defects at 13% of protein-coding genes - genes that are generally more sensitive to rapid TFIID depletion. However, transcription of these 'coactivator-redundant' genes is strongly affected by rapid depletion of both factors, showing the overlapping functions of TFIID and SAGA at this gene set. We suggest that this overlapping function is linked to TBP-DNA recruitment. The remaining 87% of expressed genes that we term 'TFIID-dependent' are highly sensitive to rapid TFIID depletion and insensitive to rapid SAGA depletion. Genome-wide mapping of SAGA and TFIID found binding of both factors at many genes independent of gene class. Promoter analysis suggests that the distinction between the gene classes is due to multiple components rather than any single regulatory factor or promoter sequence motif.

Chromosome-wide co-fluctuation of stochastic gene expression in mammalian cells

Gene expression is subject to stochastic noise, but to what extent and by which means such stochastic variations are coordinated among different genes are unclear. We hypothesize that neighboring genes on the same chromosome co-fluctuate in expression because of their common chromatin dynamics, and verify it at the genomic scale using allele-specific single-cell RNA-sequencing data of mouse cells. Unexpectedly, the co-fluctuation extends to genes that are over 60 million bases apart. We provide evidence that this long-range effect arises in part from chromatin co-accessibilities of linked loci attributable to three-dimensional proximity, which is much closer intra-chromosomally than inter-chromosomally. We further show that genes encoding components of the same protein complex tend to be chromosomally linked, likely resulting from natural selection for intracellular among-component dosage balance. These findings have implications for both the evolution of genome organization and optimal design of synthetic genomes in the face of gene expression noise.

Considering Abundance, Affinity, and Binding Site Availability in the NF-κB Target Selection Puzzle

The NF-κB transcription regulation system governs a diverse set of responses to various cytokine stimuli. With tools from in vitro biochemical characterizations, to omics-based whole genome investigations, great strides have been made in understanding how NF-κB transcription factors control the expression of specific sets of genes. Nonetheless, these efforts have also revealed a very large number of potential binding sites for NF-κB in the human genome, and a puzzle emerges when trying to explain how NF-κB selects from these many binding sites to direct cell-type- and stimulus-specific gene expression patterns. In this review, we surmise that target gene transcription can broadly be thought of as a function of the nuclear abundance of the various NF-κB dimers, the affinity of NF-κB dimers for the regulatory sequence and the availability of this regulatory site. We use this framework to place quantitative information that has been gathered about the NF-κB transcription regulation system into context and thus consider questions it answers, and questions it raises. We end with a brief discussion of some of the future prospects that new approaches could bring to our understanding of how NF-κB transcription factors orchestrate diverse responses in different biological contexts.

An enriched network motif family regulates multistep cell fate transitions with restricted reversibility

Multistep cell fate transitions with stepwise changes of transcriptional profiles are common to many developmental, regenerative and pathological processes. The multiple intermediate cell lineage states can serve as differentiation checkpoints or branching points for channeling cells to more than one lineages. However, mechanisms underlying these transitions remain elusive. Here, we explored gene regulatory circuits that can generate multiple intermediate cellular states with stepwise modulations of transcription factors. With unbiased searching in the network topology space, we found a motif family containing a large set of networks can give rise to four attractors with the stepwise regulations of transcription factors, which limit the reversibility of three consecutive steps of the lineage transition. We found that there is an enrichment of these motifs in a transcriptional network controlling the early T cell development, and a mathematical model based on this network recapitulates multistep transitions in the early T cell lineage commitment. By calculating the energy landscape and minimum action paths for the T cell model, we quantified the stochastic dynamics of the critical factors in response to the differentiation signal with fluctuations. These results are in good agreement with experimental observations and they suggest the stable characteristics of the intermediate states in the T cell differentiation. These dynamical features may help to direct the cells to correct lineages during development. Our findings provide general design principles for multistep cell linage transitions and new insights into the early T cell development. The network motifs containing a large family of topologies can be useful for analyzing diverse biological systems with multistep transitions.

The functions of cells are dynamically controlled in many biological processes including development, regeneration and disease progression. Cell fate transition, or the switch of cellular functions, often involves multiple steps. The intermediate stages of the transition provide the biological systems with the opportunities to regulate the transitions in a precise manner. These transitions are controlled by key regulatory genes of which the expression shows stepwise patterns, but how the interactions of these genes can determine the multistep processes was unclear. Here, we present a comprehensive analysis on the design principles of gene circuits that govern multistep cell fate transition. We found a large network family with common structural features that can generate systems with the ability to control three consecutive steps of the transition. We found that this type of networks is enriched in a gene circuit controlling the development of T lymphocyte, a crucial type of immune cells. We performed mathematical modeling using this gene circuit and we recapitulated the stepwise and irreversible loss of stem cell properties of the developing T lymphocytes. Our findings can be useful to analyze a wide range of gene regulatory networks controlling multistep cell fate transitions.

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.

The stress-regulatory transcription factors Msn2 and Msn4 regulate fatty acid oxidation in budding yeast

The transcription factors Msn2 and Msn4 (multicopy suppressor of SNF1 mutation proteins 2 and 4) bind the stress-response element in gene promoters in the yeast Saccharomyces cerevisiae However, the roles of Msn2/4 in primary metabolic pathways such as fatty acid β-oxidation are unclear. Here, in silico analysis revealed that the promoters of most genes involved in the biogenesis, function, and regulation of the peroxisome contain Msn2/4-binding sites. We also found that transcript levels of MSN2/MSN4 are increased in glucose-depletion conditions and that during growth in nonpreferred carbon sources, Msn2 is constantly localized to the nucleus in wild-type cells. Of note, the double mutant msn2Δmsn4Δ exhibited a severe growth defect when grown with oleic acid as the sole carbon source and had reduced transcript levels of major β-oxidation genes. ChIP indicated that Msn2 has increased occupancy on the promoters of β-oxidation genes in glucose-depleted conditions, and in vivo reporter gene analysis indicated reduced expression of these genes in msn2Δmsn4Δ cells. Moreover, mobility shift assays revealed that Msn4 binds β-oxidation gene promoters. Immunofluorescence microscopy with anti-peroxisome membrane protein antibodies disclosed that the msn2Δmsn4Δ strain had fewer peroxisomes than the wild type, and lipid analysis indicated that the msn2Δmsn4Δ strain had increased triacylglycerol and steryl ester levels. Collectively, our data suggest that Msn2/Msn4 transcription factors activate expression of the genes involved in fatty acid oxidation. Because glucose sensing, signaling, and fatty acid β-oxidation pathways are evolutionarily conserved throughout eukaryotes, the msn2Δmsn4Δ strain could therefore be a good model system for further study of these critical processes.

The haploinsufficient tumor suppressor, CUX1, acts as an analog transcriptional regulator that controls target genes through distal enhancers that loop to target promoters

One third of tumor suppressors are haploinsufficient transcriptional regulators, yet it remains unknown how a 50% reduction of a transcription factor is translated at the cis-regulatory level into a malignant transcriptional program. We studied CUX1, a haploinsufficient transcription factor that is recurrently mutated in hematopoietic and solid tumors. We determined CUX1 DNA-binding and target gene regulation in the wildtype and haploinsufficient states. CUX1 binds with transcriptional activators and cohesin at distal enhancers across three different human cell types. Haploinsufficiency of CUX1 altered the expression of a large number of genes, including cell cycle regulators, with concomitant increased cellular proliferation. Surprisingly, CUX1 occupancy decreased genome-wide in the haploinsufficient state, and binding site affinity did not correlate with differential gene expression. Instead, differentially expressed genes had multiple, low-affinity CUX1 binding sites, features of analog gene regulation. A machine-learning algorithm determined that chromatin accessibility, enhancer activity, and distance to the transcription start site are features of dose-sensitive CUX1 transcriptional regulation. Moreover, CUX1 is enriched at sites of DNA looping, as determined by Hi-C analysis, and these loops connect CUX1 to the promoters of regulated genes. We propose an analog model for haploinsufficient transcriptional deregulation mediated by higher order genome architecture.

Synthetic circuits that process multiple light and chemical signal inputs

Background

Multi-signal processing circuits are essential for rational design of sophisticated synthetic systems with good controllability and modularity, therefore, enable construction of high-level networks. Moreover, light-inducible systems provide fast and reversible means for spatiotemporal control of gene expression.

Results

Here, in HEK 293 cells, we present combinatory genetic circuits responding to light and chemical signals, simultaneously. We first constructed a dual input circuit converting different light intensities into varying of the sensitivity of the promoter to a chemical inducer (doxycycline). Next, we generated a ternary input circuit, which responded to light, doxycycline and cumate. This circuit allowed us to use different combinations of blue light and the two chemical inducers to generate gradual output values over two orders of magnitude.

Conclusions

Overall, in this study, we devise genetic circuits sensing and processing light and chemical inducers. Our work may provide insights into bio-computation and fine-tuning expression of the transgene.

Electronic supplementary material

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

Model-guided optogenetic study of PKA signaling in budding yeast

Optogenetic activation of the adenylate cyclase enzyme in Saccharomyces cerevisiae, paired with computational modeling, enables study of the dynamic quantitative properties of the cAMP/PKA signaling network. The ability to deliver such precise perturbation reveals fundamental dynamical features of PKA signaling, including the time scales of feedback.

In eukaryotes, protein kinase A (PKA) is a master regulator of cell proliferation and survival. The activity of PKA is subject to elaborate control and exhibits complex time dynamics. To probe the quantitative attributes of PKA dynamics in the yeast Saccharomyces cerevisiae, we developed an optogenetic strategy that uses a photoactivatable adenylate cyclase to achieve real-time regulation of cAMP and the PKA pathway. We capitalize on the precise and rapid control afforded by this optogenetic tool, together with quantitative computational modeling, to study the properties of feedback in the PKA signaling network and dissect the nonintuitive dynamic effects that ensue from perturbing its components. Our analyses reveal that negative feedback channeled through the Ras1/2 GTPase is delayed, pinpointing its time scale and its contribution to the dynamic features of the cAMP/PKA signaling network.

Dynamic control of gene regulatory logic by seemingly redundant transcription factors

Many transcription factors co-express with their homologs to regulate identical target genes, however the advantages of such redundancies remain elusive. Using single-cell imaging and microfluidics, we study the yeast general stress response transcription factor Msn2 and its seemingly redundant homolog Msn4. We find that gene regulation by these two factors is analogous to logic gate systems. Target genes with fast activation kinetics can be fully induced by either factor, behaving as an 'OR' gate. In contrast, target genes with slow activation kinetics behave as an 'AND' gate, requiring distinct contributions from both factors, upon transient stimulation. Furthermore, such genes become an 'OR' gate when the input duration is prolonged, suggesting that the logic gate scheme is not static but rather dependent on the input dynamics. Therefore, Msn2 and Msn4 enable a time-based mode of combinatorial gene regulation that might be applicable to homologous transcription factors in other organisms.

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

Cells respond to environmental signals by activating proteins called transcription factors. These bind to the DNA that is stored in the cell nucleus and turn on specific genes to make gene products. Many of these transcription factors move in and out of the nucleus once activated. Different environmental signals affect the amount of transcription factor that appears in the nucleus in different ways, and this is important in determining which genes should be turned on and how many copies of gene products should be made.

Many transcription factors co-exist with a similar version of themselves in the same cell. These closely related proteins, called homologous transcription factors, respond to the same signals and bind to the same place on the DNA to turn on the same genes. It was not clear what advantages the cells gain from having two molecules that perform the same roles.

Two homologous transcription factors called Msn2 and Msn4 are found in baker's yeast. These transcription factors respond to a wide variety of environmental stresses by moving rapidly into the nucleus, where they remain for a short time to turn on hundreds of target genes that are needed for the cell to survive.

AkhavanAghdam, Sinha, Tabbaa et al. investigated the roles of Msn2 and Msn4 by tracking where the proteins localized to and which genes they switched on inside the same single cell. Genes that can be turned on quickly could be activated by either Msn2 or Msn4, and both factors activated the genes to a similar extent. By contrast, both Msn2 and Msn4 were required to activate those genes that take a long time to be turned on. In these cases, Msn2 served as a 'switch' that governed the 'on' and 'off' state of the genes, while Msn4 behaved as a 'rheostat' to tune how much gene product was made. This cooperation between the two transcription factors is equivalent to a design commonly found in electrical circuits and may help the cell to survive in rapidly changing environments.

Further studies are now needed to investigate the mechanisms that provide Msn2 and Msn4 with distinct roles in gene regulation. Technological advances that allow the full genetic material of a single cell to be analyzed could also determine whether other homologous transcription factors regulate their target genes in similar ways.

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

Adaptation through proportion

Adaptation is a ubiquitous feature in biological sensory and signaling networks. It has been suggested that adaptive systems may follow certain simple design principles across diverse organisms, cells and pathways. One class of networks that can achieve adaptation utilizes an incoherent feedforward control, in which two parallel signaling branches exert opposite but proportional effects on the output at steady state. In this paper, we generalize this adaptation mechanism by establishing a steady-state proportionality relationship among a subset of nodes in a network. Adaptation can be achieved by using any two nodes in the sub-network to respectively regulate the output node positively and negatively. We focus on enzyme networks and first identify basic regulation motifs consisting of two and three nodes that can be used to build small networks with proportional relationships. Larger proportional networks can then be constructed modularly similar to LEGOs. Our method provides a general framework to construct and analyze a class of proportional and/or adaptation networks with arbitrary size, flexibility and versatile functional features.

Computational Framework for Modeling Multiple Noncooperative Transcription Factor Binding and Its Application to the Analysis of Nuclear Factor Kappa B Oscillatory Response

Recent studies have shown that regulation of many important genes is achieved with multiple transcription factor (TF) binding sites with low or no cooperativity. Additionally, noncooperative binding sites are gaining more and more importance in the field of synthetic biology. Here, we introduce a computational framework that can be applied to dynamical modeling and analysis of gene regulatory networks with multiple noncooperative TF binding sites. We propose two computational methods to be used within the framework, that is, average promoter state approximation and expression profiles based modeling. We demonstrate the application of the proposed framework on the analysis of nuclear factor kappa B (NF-κB) oscillatory response. We show that different promoter expression hypotheses in a combination with the number of TF binding sites drastically affect the dynamics of the observed system and should not be ignored in the process of quantitative dynamical modeling, as is usually the case in existent state-of-the-art computational analyses.

Probing and rearranging the transcription factor network controlling the C. elegans endoderm

The ELT-2 GATA factor is the predominant transcription factor regulating gene expression in the C. elegans intestine, following endoderm specification. We comment on our previous study (Wiesenfahrt et al., 2016) that investigated how the elt-2 gene is controlled by END-1, END-3 and ELT-7, the 3 endoderm specific GATA factors that lie upstream in the regulatory hierarchy. We also discuss the unexpected result that ELT-2, if expressed sufficiently early and at sufficiently high levels, can specify the C. elegans endoderm, replacing the normal functions of END-1 and END-3.

Dynamics of epigenetic regulation at the single-cell level

Chromatin regulators play a major role in establishing and maintaining gene expression states. Yet how they control gene expression in single cells, quantitatively and over time, remains unclear. We used time-lapse microscopy to analyze the dynamic effects of four silencers associated with diverse modifications: DNA methylation, histone deacetylation, and histone methylation. For all regulators, silencing and reactivation occurred in all-or-none events, enabling the regulators to modulate the fraction of cells silenced rather than the amount of gene expression. These dynamics could be described by a three-state model involving stochastic transitions between active, reversibly silent, and irreversibly silent states. Through their individual transition rates, these regulators operate over different time scales and generate distinct types of epigenetic memory. Our results provide a framework for understanding and engineering mammalian chromatin regulation and epigenetic memory.

Beyond antioxidant genes in the ancient Nrf2 regulatory network

Nrf2, a basic leucine zipper transcription factor encoded by the gene NFE2L2, is a master regulator of the transcriptional response to oxidative stress. Nrf2 is structurally and functionally conserved from insects to humans, and it heterodimerizes with the small MAF transcription factors to bind a consensus DNA sequence (the antioxidant response element, or ARE) and regulate gene expression. We have used genome-wide chromatin immunoprecipitation and gene expression data to identify direct Nrf2 target genes in Drosophila and humans. These data have allowed us to construct the deeply conserved ancient Nrf2 regulatory network-target genes that are conserved from Drosophila to human. The ancient network consists of canonical antioxidant genes, as well as genes related to proteasomal pathways and metabolism and a number of less expected genes. We have also used enhancer reporter assays and electrophoretic mobility-shift assays to confirm Nrf2-mediated regulation of ARE activity at a number of these novel target genes. Interestingly, the ancient network also highlights a prominent negative feedback loop; this, combined with the finding that Nrf2-mediated regulatory output is tightly linked to the quality of the ARE it is targeting, suggests that precise regulation of nuclear Nrf2 concentration is necessary to achieve proper quantitative regulation of distinct gene sets. Together, these findings highlight the importance of balance in the Nrf2-ARE pathway and indicate that Nrf2-mediated regulation of xenobiotic metabolism, glucose metabolism, and proteostasis has been central to this pathway since its inception.

Different Mechanisms Confer Gradual Control and Memory at Nutrient- and Stress-Regulated Genes in Yeast

Cells respond to environmental stimuli by fine-tuned regulation of gene expression. Here we investigated the dose-dependent modulation of gene expression at high temporal resolution in response to nutrient and stress signals in yeast. The GAL1 activity in cell populations is modulated in a well-defined range of galactose concentrations, correlating with a dynamic change of histone remodeling and RNA polymerase II (RNAPII) association. This behavior is the result of a heterogeneous induction delay caused by decreasing inducer concentrations across the population. Chromatin remodeling appears to be the basis for the dynamic GAL1 expression, because mutants with impaired histone dynamics show severely truncated dose-response profiles. In contrast, the GRE2 promoter operates like a rapid off/on switch in response to increasing osmotic stress, with almost constant expression rates and exclusively temporal regulation of histone remodeling and RNAPII occupancy. The Gal3 inducer and the Hog1 mitogen-activated protein (MAP) kinase seem to determine the different dose-response strategies at the two promoters. Accordingly, GAL1 becomes highly sensitive and dose independent if previously stimulated because of residual Gal3 levels, whereas GRE2 expression diminishes upon repeated stimulation due to acquired stress resistance. Our analysis reveals important differences in the way dynamic signals create dose-sensitive gene expression outputs.

cis Determinants of Promoter Threshold and Activation Timescale

Although the relationship between DNA cis-regulatory sequences and gene expression has been extensively studied at steady state, how cis-regulatory sequences affect the dynamics of gene induction is not known. The dynamics of gene induction can be described by the promoter activation timescale (AcTime) and amplitude threshold (AmpThr). Combining high-throughput microfluidics with quantitative time-lapse microscopy, we control the activation dynamics of the budding yeast transcription factor, Msn2, and reveal how cis-regulatory motifs in 20 promoter variants of the Msn2-target-gene SIP18 affect AcTime and AmpThr. By modulating Msn2 binding sites, we can decouple AmpThr from AcTime and switch the SIP18 promoter class from high AmpThr and slow AcTime to low AmpThr and either fast or slow AcTime. We present a model that quantitatively explains gene-induction dynamics on the basis of the Msn2-binding-site number, TATA box location, and promoter nucleosome organization. Overall, we elucidate the cis-regulatory logic underlying promoter decoding of TF dynamics.

Photoaffinity labeling of transcription factors by DNA-templated crosslinking

Characterization of transcription factor-DNA interaction is of high importance in elucidating the molecular mechanisms of gene transcriptions. DNA-based affinity probes were developed to capture and identify transcription factors by covalent crosslinking; however, the requirement of a crosslinker on the affinity probe remains a disadvantage, as the crosslinker itself often interferes with the protein-DNA interactions. We report a dual-probe method able to capture DNA-binding transcription factors with unmodified protein-binding sites in scenarios where conventional probes have failed. We have also shown the method's converse application in selecting specific transcription factor-binding DNA sequences from a probe library and its extension to studying proteins recognizing epigenetic marks. This study may provide a new tool for exploring DNA-binding proteins in biology.

Low affinity binding site clusters confer hox specificity and regulatory robustness

In animals, Hox transcription factors define regional identity in distinct anatomical domains. How Hox genes encode this specificity is a paradox, because different Hox proteins bind with high affinity in vitro to similar DNA sequences. Here, we demonstrate that the Hox protein Ultrabithorax (Ubx) in complex with its cofactor Extradenticle (Exd) bound specifically to clusters of very low affinity sites in enhancers of the shavenbaby gene of Drosophila. These low affinity sites conferred specificity for Ubx binding in vivo, but multiple clustered sites were required for robust expression when embryos developed in variable environments. Although most individual Ubx binding sites are not evolutionarily conserved, the overall enhancer architecture-clusters of low affinity binding sites-is maintained and required for enhancer function. Natural selection therefore works at the level of the enhancer, requiring a particular density of low affinity Ubx sites to confer both specific and robust expression.

Pulsatile dynamics in the yeast proteome

The activation of transcription factors in response to environmental conditions is fundamental to cellular regulation. Recent work has revealed that some transcription factors are activated in stochastic pulses of nuclear localization, rather than at a constant level, even in a constant environment [1-12]. In such cases, signals control the mean activity of the transcription factor by modulating the frequency, duration, or amplitude of these pulses. Although specific pulsatile transcription factors have been identified in diverse cell types, it has remained unclear how prevalent pulsing is within the cell, how variable pulsing behaviors are between genes, and whether pulsing is specific to transcriptional regulators or is employed more broadly. To address these issues, we performed a proteome-wide movie-based screen to systematically identify localization-based pulsing behaviors in Saccharomyces cerevisiae. The screen examined all genes in a previously developed fluorescent protein fusion library of 4,159 strains [13] in multiple media conditions. This approach revealed stochastic pulsing in ten proteins, all transcription factors. In each case, pulse dynamics were heterogeneous and unsynchronized among cells in clonal populations. Pulsing is the only dynamic localization behavior that we observed, and it tends to occur in pairs of paralogous and redundant proteins. Taken together, these results suggest that pulsatile dynamics play a pervasive role in yeast and may be similarly prevalent in other eukaryotic species.

Gene regulation: when analog beats digital

Why do some genes seem to respond in a 'digital', on/off manner to a graded signal, while others produce an 'analog', graded response? A new study suggests that the DNA-binding properties of transcription factors can strongly influence the response patterns of gene networks.

Absence of a simple code: how transcription factors read the genome

Transcription factors (TFs) influence cell fate by interpreting the regulatory DNA within a genome. TFs recognize DNA in a specific manner; the mechanisms underlying this specificity have been identified for many TFs based on 3D structures of protein-DNA complexes. More recently, structural views have been complemented with data from high-throughput in vitro and in vivo explorations of the DNA-binding preferences of many TFs. Together, these approaches have greatly expanded our understanding of TF-DNA interactions. However, the mechanisms by which TFs select in vivo binding sites and alter gene expression remain unclear. Recent work has highlighted the many variables that influence TF-DNA binding, while demonstrating that a biophysical understanding of these many factors will be central to understanding TF function.