YeTFaSCo: a database of evaluated yeast transcription factor sequence specificities

Ordered and disordered regions of the Origin Recognition Complex direct differential in vivo binding at distinct motif sequences

The Origin Recognition Complex (ORC) seeds replication-fork formation by binding to DNA replication origins, which in budding yeast contain a 17bp DNA motif. High resolution structure of the ORC-DNA complex revealed two base-interacting elements: a disordered basic patch (Orc1-BP4) and an insertion helix (Orc4-IH). To define the ORC elements guiding its DNA binding in vivo, we mapped genomic locations of 38 designed ORC mutants, revealing that different ORC elements guide binding at different sites. At silencing-associated sites lacking the motif, ORC binding and activity were fully explained by a BAH domain. Within replication origins, we reveal two dominating motif variants showing differential binding modes and symmetry: a non-repetitive motif whose binding requires Orc1-BP4 and Orc4-IH, and a repetitive one where another basic patch, Orc1-BP3, can replace Orc4-IH. Disordered basic patches are therefore key for ORC-motif binding in vivo, and we discuss how these conserved, minor-groove interacting elements can guide specific ORC-DNA recognition.

Species-aware DNA language models capture regulatory elements and their evolution

BACKGROUND

The rise of large-scale multi-species genome sequencing projects promises to shed new light on how genomes encode gene regulatory instructions. To this end, new algorithms are needed that can leverage conservation to capture regulatory elements while accounting for their evolution.

RESULTS

Here, we introduce species-aware DNA language models, which we trained on more than 800 species spanning over 500 million years of evolution. Investigating their ability to predict masked nucleotides from context, we show that DNA language models distinguish transcription factor and RNA-binding protein motifs from background non-coding sequence. Owing to their flexibility, DNA language models capture conserved regulatory elements over much further evolutionary distances than sequence alignment would allow. Remarkably, DNA language models reconstruct motif instances bound in vivo better than unbound ones and account for the evolution of motif sequences and their positional constraints, showing that these models capture functional high-order sequence and evolutionary context. We further show that species-aware training yields improved sequence representations for endogenous and MPRA-based gene expression prediction, as well as motif discovery.

CONCLUSIONS

Collectively, these results demonstrate that species-aware DNA language models are a powerful, flexible, and scalable tool to integrate information from large compendia of highly diverged genomes.

Hungry for control: metabolite signaling to chromatin in Plasmodium falciparum

The human malaria parasite Plasmodium falciparum undergoes a complex life cycle in two hosts, mammalian and mosquito, where it is constantly subjected to environmental changes in nutrients. Epigenetic mechanisms govern transcriptional switches and are essential for parasite persistence and proliferation. Parasites infecting red blood cells are auxotrophic for several nutrients, and mounting evidence suggests that various metabolites act as direct substrates for epigenetic modifications, with their abundance directly relating to changes in parasite gene expression. Here, we review the latest understanding of metabolic changes that alter the histone code resulting in changes to transcriptional programmes in malaria parasites.

Diverse mechanisms control amino acid-dependent environmental alkalization by Candida albicans

Candida albicans has the capacity to neutralize acidic growth environments by releasing ammonia derived from the catabolism of amino acids. The molecular components underlying alkalization and its physiological significance remain poorly understood. Here, we present an integrative model with the cytosolic NAD+-dependent glutamate dehydrogenase (Gdh2) as the principal ammonia-generating component. We show that alkalization is dependent on the SPS-sensor-regulated transcription factor STP2 and the proline-responsive activator Put3. These factors function in parallel to derepress GDH2 and the two proline catabolic enzymes PUT1 and PUT2. Consistently, a double mutant lacking STP2 and PUT3 exhibits a severe alkalization defect that nearly phenocopies that of a gdh2-/- strain. Alkalization is dependent on mitochondrial activity and in wild-type cells occurs as long as the conditions permit respiratory growth. Strikingly, Gdh2 levels decrease and cells transiently extrude glutamate as the environment becomes more alkaline. Together, these processes constitute a rudimentary regulatory system that counters and limits the negative effects associated with ammonia generation. These findings align with Gdh2 being dispensable for virulence, and based on a whole human blood virulence assay, the same is true for C. glabrata and C. auris. Using a transwell co-culture system, we observed that the growth and proliferation of Lactobacillus crispatus, a common component of the acidic vaginal microenvironment and a potent antagonist of C. albicans, is unaffected by fungal-induced alkalization. Consequently, although Candida spp. can alkalinize their growth environments, other fungal-associated processes are more critical in promoting dysbiosis and virulent fungal growth.

A universal system for boosting gene expression in eukaryotic cell-lines

We demonstrate a transcriptional regulatory design algorithm that can boost expression in yeast and mammalian cell lines. The system consists of a simplified transcriptional architecture composed of a minimal core promoter and a synthetic upstream regulatory region (sURS) composed of up to three motifs selected from a list of 41 motifs conserved in the eukaryotic lineage. The sURS system was first characterized using an oligo-library containing 189,990 variants. We validate the resultant expression model using a set of 43 unseen sURS designs. The validation sURS experiments indicate that a generic set of grammar rules for boosting and attenuation may exist in yeast cells. Finally, we demonstrate that this generic set of grammar rules functions similarly in mammalian CHO-K1 and HeLa cells. Consequently, our work provides a design algorithm for boosting the expression of promoters used for expressing industrially relevant proteins in yeast and mammalian cell lines.

In vitro reconstitution of chromatin domains shows a role for nucleosome positioning in 3D genome organization

Eukaryotic genomes are organized into chromatin domains. The molecular mechanisms driving the formation of these domains are difficult to dissect in vivo and remain poorly understood. Here we reconstitute Saccharomyces cerevisiae chromatin in vitro and determine its 3D organization at subnucleosome resolution by micrococcal nuclease-based chromosome conformation capture and molecular dynamics simulations. We show that regularly spaced and phased nucleosome arrays form chromatin domains in vitro that resemble domains in vivo. This demonstrates that neither loop extrusion nor transcription is required for basic domain formation in yeast. In addition, we find that the boundaries of reconstituted domains correspond to nucleosome-free regions and that insulation strength scales with their width. Finally, we show that domain compaction depends on nucleosome linker length, with longer linkers forming more compact structures. Together, our results demonstrate that regular nucleosome positioning is important for the formation of chromatin domains and provide a proof-of-principle for bottom-up 3D genome studies.

Regulatory activity is the default DNA state in eukaryotes

Genomes encode for genes and non-coding DNA, both capable of transcriptional activity. However, unlike canonical genes, many transcripts from non-coding DNA have limited evidence of conservation or function. Here, to determine how much biological noise is expected from non-genic sequences, we quantify the regulatory activity of evolutionarily naive DNA using RNA-seq in yeast and computational predictions in humans. In yeast, more than 99% of naive DNA bases were transcribed. Unlike the evolved transcriptome, naive transcripts frequently overlapped with opposite sense transcripts, suggesting selection favored coherent gene structures in the yeast genome. In humans, regulation-associated chromatin activity is predicted to be common in naive dinucleotide-content-matched randomized DNA. Here, naive and evolved DNA have similar co-occurrence and cell-type specificity of chromatin marks, challenging these as indicators of selection. However, in both yeast and humans, extreme high activities were rare in naive DNA, suggesting they result from selection. Overall, basal regulatory activity seems to be the default, which selection can hone to evolve a function or, if detrimental, repress.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

Overexpression of the Transcription Factor Azf1 Reveals Novel Regulatory Functions and Impacts β-Glucosidase Production in Trichoderma reesei

The fungus Trichoderma reesei is an essential producer of enzymes that degrade lignocellulosic biomass to produce value-added bioproducts. The cellulolytic system of T. reesei is controlled by several transcription factors (TFs) that efficiently regulate the production of these enzymes. Recently, a new TF named Azf1 was identified as a positive regulator of cellulase expression. Here, we investigated novel regulatory functions of Azf1 by its overexpression. In the mutant strain OEazf1, overexpression of azf1 was achieved under both repression and induction conditions. Although azf1 was more abundant in transcript and protein, overexpression of this TF did not activate transcription of the cellulase gene in the presence of the repressor glucose, suggesting that Azf1 may be subject to posttranslational regulation. In cellulose, the expression of swo, encoding the accessory protein swollenin, and the β-glucosidases cel1a, cel1b, cel3b, and cel3g increases in the early stages of cultivation. The increased production of these β-glucosidases increases the hydrolysis rate of cellobiose and sophorose, which activates carbon catabolite repression (CCR) and causes repression of cellulase genes and the regulator Xyr1 in the later stages of cultivation. Moreover, overexpression of azf1 led to increased cellulase activity in T. reesei during long-term cultivation in cellulose and sugarcane bagasse. Our results provide new insights into the mechanisms regulating Azf1 and novel genes that are important targets of this TF. This work contributes to a better understanding of the complex mechanisms regulating cellulase expression in T. reesei. It will contribute to the development of strains with higher production of these essential enzymes.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

Metronidazole Potentiation by Panax Ginseng and Symphytum officinale: A New Strategy for P. gingivalis Infection Control

The important periodontal disease pathogen Porphyromonas gingivalis produces thick biofilms that increase its pathogenicity. Finding natural antimicrobial agents is crucial because of the rise in antibiotic resistance. The purpose of this study was to determine if plant extracts such as Symphytum officinale (S) and Panax Ginseng (G) were effective against P. gingivalis separately and in combination with a common antibiotic, metronidazole (F). Six different dilutions were produced using the plant extracts in different concentrations and antibiotics separately and in combination with F, G, and S using the two-fold serial dilution technique. To evaluate the effects of these substances, biofilm inhibition experiments were conducted. Plaque samples were collected from periodontitis patients to isolate P. gingivalis, and a standard strain of P. gingivalis (ATCC 33277) was purchased. Additionally, Acylated Homoserine Lactones (AHLs) detection was carried out to look for any activity that would interfere with quorum sensing. GraphPad Prism was used for statistical analysis with a p-value < 0.05. The combinations of Symphytum officinale and metronidazole (S+F) showed the maximum effectiveness in biofilm inhibition (98.7%), which was slightly better than G+F (98.2%), with substantial variations in biofilm inhibition levels in different treatment regimes. Notably, the patient isolate was more active than the standard strain. Additionally, the plant extracts and their combinations at particular dilutions had notable inhibitory effects on the generation of AHL (p < 0.05). The study highlights the possibility of Symphytum officinale and Panax Ginseng as effective treatments for P. gingivalis biofilm and AHLs, both on their own and in combination with metronidazole. These organic substances may open the door to cutting-edge methods of treating periodontal disorders.

Adaptive responses of yeast strains tolerant to acidic pH, acetate, and supraoptimal temperature

Ethanol fermentations can be prematurely halted as Saccharomyces cerevisiae faces adverse conditions, such as acidic pH, presence of acetic acid, and supraoptimal temperatures. The knowledge on yeast responses to these conditions is essential to endowing a tolerant phenotype to another strain by targeted genetic manipulation. In this study, physiological and whole-genome analyses were conducted to obtain insights on molecular responses which potentially render yeast tolerant towards thermoacidic conditions. To this end, we used thermotolerant TTY23, acid tolerant AT22, and thermo-acid tolerant TAT12 strains previously generated by adaptive laboratory evolution (ALE) experiments. The results showed an increase in thermoacidic profiles in the tolerant strains. The whole-genome sequence revealed the importance of genes related to: H⁺, iron, and glycerol transport (i.e., PMA1, FRE1/2, JEN1, VMA2, VCX1, KHA1, AQY3, and ATO2); transcriptional regulation of stress responses to drugs, reactive oxygen species and heat-shock (i.e., HSF1, SKN7, BAS1, HFI1, and WAR1); and adjustments of fermentative growth and stress responses by glucose signaling pathways (i.e., ACS1, GPA1/2, RAS2, IRA2, and REG1). At 30 °C and pH 5.5, more than a thousand differentially expressed genes (DEGs) were identified in each strain. The integration of results revealed that evolved strains adjust their intracellular pH by H⁺ and acetic acid transport, modify their metabolism and stress responses via glucose signaling pathways, control of cellular ATP pools by regulating translation and de novo synthesis of nucleotides, and direct the synthesis, folding and rescue of proteins throughout the heat-shock stress response. Moreover, the motifs analysis in mutated transcription factors suggested a significant association of SFP1, YRR1, BAS1, HFI1, HSF1, and SKN7 TFs with DEGs found in thermoacidic tolerant yeast strains. KEY POINTS: • All the evolved strains overexpressed the plasma membrane H⁺ -ATPase PMA1 at optimal conditions • Tolerant strain TAT12 mutated genes encoding weak acid and heat response TFs HSF1, SKN7, and WAR1 • TFs HSF1 and SKN7 likely controlled the transcription of metabolic genes associated to heat and acid tolerance.

Substrate-specific effects of natural genetic variation on proteasome activity

Protein degradation is an essential biological process that regulates protein abundance and removes misfolded and damaged proteins from cells. In eukaryotes, most protein degradation occurs through the stepwise actions of two functionally distinct entities, the ubiquitin system and the proteasome. Ubiquitin system enzymes attach ubiquitin to cellular proteins, targeting them for degradation. The proteasome then selectively binds and degrades ubiquitinated substrate proteins. Genetic variation in ubiquitin system genes creates heritable differences in the degradation of their substrates. However, the challenges of measuring the degradative activity of the proteasome independently of the ubiquitin system in large samples have limited our understanding of genetic influences on the proteasome. Here, using the yeast Saccharomyces cerevisiae, we built and characterized reporters that provide high-throughput, ubiquitin system-independent measurements of proteasome activity. Using single-cell measurements of proteasome activity from millions of genetically diverse yeast cells, we mapped 15 loci across the genome that influence proteasomal protein degradation. Twelve of these 15 loci exerted specific effects on the degradation of two distinct proteasome substrates, revealing a high degree of substrate-specificity in the genetics of proteasome activity. Using CRISPR-Cas9-based allelic engineering, we resolved a locus to a causal variant in the promoter of RPT6, a gene that encodes a subunit of the proteasome's 19S regulatory particle. The variant increases RPT6 expression, which we show results in increased proteasome activity. Our results reveal the complex genetic architecture of proteasome activity and suggest that genetic influences on the proteasome may be an important source of variation in the many cellular and organismal traits shaped by protein degradation.

Predicting nucleosome positioning using statistical equilibrium models in budding yeast

We present a protocol using thermodynamic models to predict nucleosome positioning with transcription factors (TFs) and chromatin remodelers. We describe step-by-step approaches to annotate genome-wide nucleosome-depleted regions (NDRs), compute nucleosome and TF occupancy, optimize parameters, and evaluate model performance. These models identify nucleosome-displacing TFs in the budding yeast genome and predict the locations and sizes of NDRs solely based on DNA sequence and TF motifs. The protocol can be applied to all organisms with prior knowledge of TF motifs. For complete details on the use and execution of this protocol, please refer to Kharerin and Bai (2021).1.

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

Natural Variation in Diauxic Shift between Patagonian Saccharomyces eubayanus Strains

The study of natural variation can untap novel alleles with immense value for biotechnological applications. Saccharomyces eubayanus Patagonian isolates exhibit differences in the diauxic shift between glucose and maltose, representing a suitable model to study their natural genetic variation for novel strains for brewing. However, little is known about the genetic variants and chromatin regulators responsible for these differences. Here, we show how genome-wide chromatin accessibility and gene expression differences underlie distinct diauxic shift profiles in S. eubayanus. We identified two strains with a rapid diauxic shift between glucose and maltose (CL467.1 and CBS12357) and one strain with a remarkably low fermentation efficiency and longer lag phase during diauxic shift (QC18). This is associated in the QC18 strain with lower transcriptional activity and chromatin accessibility of specific genes of maltose metabolism and higher expression levels of glucose transporters. These differences are governed by the HAP complex, which differentially regulates gene expression depending on the genetic background. We found in the QC18 strain a contrasting phenotype to those phenotypes described in S. cerevisiae, where hap4Δ, hap5Δ, and cin5Δ knockouts significantly improved the QC18 growth rate in the glucose-maltose shift. The most profound effects were found between CIN5 allelic variants, suggesting that Cin5p could strongly activate a repressor of the diauxic shift in the QC18 strain but not necessarily in the other strains. The differences between strains could originate from the tree host from which the strains were obtained, which might determine the sugar source preference and the brewing potential of the strain. IMPORTANCE The diauxic shift has been studied in budding yeast under laboratory conditions; however, few studies have addressed the diauxic shift between carbon sources under fermentative conditions. Here, we study the transcriptional and chromatin structure differences that explain the natural variation in fermentative capacity and efficiency during diauxic shift of natural isolates of S. eubayanus. Our results show how natural genetic variants in transcription factors impact sugar consumption preferences between strains. These variants have different effects depending on the genetic background, with a contrasting phenotype to those phenotypes previously described in S. cerevisiae. Our study shows how relatively simple genetic/molecular modifications/editing in the lab can facilitate the study of natural variations of microorganisms for the brewing industry.

Transcription factor binding process is the primary driver of noise in gene expression

Noise in expression of individual genes gives rise to variations in activity of cellular pathways and generates heterogeneity in cellular phenotypes. Phenotypic heterogeneity has important implications for antibiotic persistence, mutation penetrance, cancer growth and therapy resistance. Specific molecular features such as the presence of the TATA box sequence and the promoter nucleosome occupancy have been associated with noise. However, the relative importance of these features in noise regulation is unclear and how well these features can predict noise has not yet been assessed. Here through an integrated statistical model of gene expression noise in yeast we found that the number of regulating transcription factors (TFs) of a gene was a key predictor of noise, whereas presence of the TATA box and the promoter nucleosome occupancy had poor predictive power. With an increase in the number of regulatory TFs, there was a rise in the number of cooperatively binding TFs. In addition, an increased number of regulatory TFs meant more overlaps in TF binding sites, resulting in competition between TFs for binding to the same region of the promoter. Through modeling of TF binding to promoter and application of stochastic simulations, we demonstrated that competition and cooperation among TFs could increase noise. Thus, our work uncovers a process of noise regulation that arises out of the dynamics of gene regulation and is not dependent on any specific transcription factor or specific promoter sequence.

N6-methyladenosine (m6A) reader Pho92 is recruited co-transcriptionally and couples translation to mRNA decay to promote meiotic fitness in yeast

N6- methyladenosine (m6A) RNA modification impacts mRNA fate primarily via reader proteins, which dictate processes in development, stress, and disease. Yet little is known about m6A function in Saccharomyces cerevisiae, which occurs solely during early meiosis. Here, we perform a multifaceted analysis of the m6A reader protein Pho92/Mrb1. Cross-linking immunoprecipitation analysis reveals that Pho92 associates with the 3'end of meiotic mRNAs in both an m6A-dependent and independent manner. Within cells, Pho92 transitions from the nucleus to the cytoplasm, and associates with translating ribosomes. In the nucleus Pho92 associates with target loci through its interaction with transcriptional elongator Paf1C. Functionally, we show that Pho92 promotes and links protein synthesis to mRNA decay. As such, the Pho92-mediated m6A-mRNA decay is contingent on active translation and the CCR4-NOT complex. We propose that the m6A reader Pho92 is loaded co-transcriptionally to facilitate protein synthesis and subsequent decay of m6A modified transcripts, and thereby promotes meiosis.

Motif and conserved module analysis in DNA (promoters, enhancers) and RNA (lncRNA, mRNA) using AlModules

Nucleic acid motifs consist of conserved and variable nucleotide regions. For functional action, several motifs are combined to modules. The tool AIModules allows identification of such motifs including combinations of them and conservation in several nucleic acid stretches. AIModules recognizes conserved motifs and combinations of motifs (modules) allowing a number of interesting biological applications such as analysis of promoter and transcription factor binding sites (TFBS), identification of conserved modules shared between several gene families, e.g. promoter regions, but also analysis of shared and conserved other DNA motifs such as enhancers and silencers, in mRNA (motifs or regulatory elements e.g. for polyadenylation) and lncRNAs. The tool AIModules presented here is an integrated solution for motif analysis, offered as a Web service as well as downloadable software. Several nucleotide sequences are queried for TFBSs using predefined matrices from the JASPAR DB or by using one's own matrices for diverse types of DNA or RNA motif discovery. Furthermore, AIModules can find TFBSs common to two or more sequences. Demanding high or low conservation, AIModules outperforms other solutions in speed and finds more modules (specific combinations of TFBS) than alternative available software. The application also searches RNA motifs such as polyadenylation site or RNA-protein binding motifs as well as DNA motifs such as enhancers as well as user-specified motif combinations ( https://bioinfo-wuerz.de/aimodules/ ; alternative entry pages: https://aimodules.heinzelab.de or https://www.biozentrum.uni-wuerzburg.de/bioinfo/computing/aimodules ). The application is free and open source whether used online, on-site, or locally.

RSC and GRFs confer promoter directionality by restricting divergent noncoding transcription

The directionality of gene promoters-the ratio of protein-coding over divergent noncoding transcription-is highly variable. How promoter directionality is controlled remains poorly understood. Here, we show that the chromatin remodelling complex RSC and general regulatory factors (GRFs) dictate promoter directionality by attenuating divergent transcription relative to protein-coding transcription. At gene promoters that are highly directional, depletion of RSC leads to a relative increase in divergent noncoding transcription and thus to a decrease in promoter directionality. We find that RSC has a modest effect on nucleosome positioning upstream in promoters at the sites of divergent transcription. These promoters are also enriched for the binding of GRFs such as Reb1 and Abf1. Ectopic targeting of divergent transcription initiation sites with GRFs or the dCas9 DNA-binding protein suppresses divergent transcription. Our data suggest that RSC and GRFs play a pervasive role in limiting divergent transcription relative to coding direction transcription. We propose that any DNA-binding factor, when stably associated with cryptic transcription start sites, forms a barrier which represses divergent transcription, thereby promoting promoter directionality.

TFLink: an integrated gateway to access transcription factor-target gene interactions for multiple species

Analysis of transcriptional regulatory interactions and their comparisons across multiple species are crucial for progress in various fields in biology, from functional genomics to the evolution of signal transduction pathways. However, despite the rapidly growing body of data on regulatory interactions in several eukaryotes, no databases exist to provide curated high-quality information on transcription factor-target gene interactions for multiple species. Here, we address this gap by introducing the TFLink gateway, which uniquely provides experimentally explored and highly accurate information on transcription factor-target gene interactions (∼12 million), nucleotide sequences and genomic locations of transcription factor binding sites (∼9 million) for human and six model organisms: mouse, rat, zebrafish, fruit fly, worm and yeast by integrating 10 resources. TFLink provides user-friendly access to data on transcription factor-target gene interactions, interactive network visualizations and transcription factor binding sites, with cross-links to several other databases. Besides containing accurate information on transcription factors, with a clear labelling of the type/volume of the experiments (small-scale or high-throughput), the source database and the original publications, TFLink also provides a wealth of standardized regulatory data available for download in multiple formats. The database offers easy access to high-quality data for wet-lab researchers, supplies data for gene set enrichment analyses and facilitates systems biology and comparative gene regulation studies. Database URL https://tflink.net/.

Partitioned usage of chromatin remodelers by nucleosome-displacing factors

Nucleosome-displacing-factors (NDFs) in yeast, similar to pioneer factors in higher eukaryotes, can open closed chromatin and generate nucleosome-depleted regions (NDRs). NDRs in yeast are also affected by ATP-dependent chromatin remodelers (CRs). However, how NDFs and CRs coordinate in nucleosome invasion and NDR formation is still unclear. Here, we design a high-throughput method to systematically study the interplay between NDFs and CRs. By combining an integrated synthetic oligonucleotide library with DNA methyltransferase-based, single-molecule nucleosome mapping, we measure the impact of CRs on NDRs generated by individual NDFs. We find that CRs are dispensable for nucleosome invasion by NDFs, and they function downstream of NDF binding to modulate the NDR length. A few CRs show high specificity toward certain NDFs; however, in most cases, CRs are recruited in a factor-nonspecific and NDR length-dependent manner. Overall, our study provides a framework to investigate how NDFs and CRs cooperate to regulate chromatin opening.

Chromatin accessibility in yeast is regulated by nucleosome-displacing-factors (NDFs) and chromatin remodelers (CRs). Chen et al. show that NDFs first invade into nucleosomes and then recruit CRs to modulate the NDR length. NDF-specific and NDR length-dependent recruitment of CRs allow partitioned usage of CRs by NDFs.

The S. cerevisiae m6A-reader Pho92 promotes timely meiotic recombination by controlling key methylated transcripts

N6-Methyladenosine (m6A), one of the most abundant internal modification of eukaryotic mRNAs, participates in the post-transcriptional control of gene expression through recruitment of specific m6A readers. In Saccharomyces cerevisiae, the m6A methyltransferase Ime4 is expressed only during meiosis and its deletion impairs this process. To elucidate how m6A control gene expression, we investigated the function of the budding yeast m6A reader Pho92. We show that Pho92 is an early meiotic factor that promotes timely meiotic progression. High-throughput RNA sequencing and mapping of Pho92-binding sites following UV-crosslinking reveal that Pho92 is recruited to specific mRNAs in an m6A-dependent manner during the meiotic prophase, preceding their down-regulation. Strikingly, point mutations altering m6A sites in mRNAs targeted by Pho92 are sufficient to delay their down-regulation and, in one case, to slow down meiotic progression. Altogether, our results indicate that Pho92 facilitate the meiotic progression by accelerating the down-regulation of timely-regulated mRNAs during meiotic recombination.

Evolution of transcription factor binding through sequence variations and turnover of binding sites

Variations in noncoding regulatory sequences play a central role in evolution. Interpreting such variations, however, remains difficult even in the context of defined attributes such as transcription factor (TF) binding sites. Here, we systematically link variations in cis-regulatory sequences to TF binding by profiling the allele-specific binding of 27 TFs expressed in a yeast hybrid, in which two related genomes are present within the same nucleus. TFs localize preferentially to sites containing their known consensus motifs but occupy only a small fraction of the motif-containing sites available within the genomes. Differential binding of TFs to the orthologous alleles was well explained by variations that alter motif sequence, whereas differences in chromatin accessibility between alleles were of little apparent effect. Motif variations that abolished binding when present in only one allele were still bound when present in both alleles, suggesting evolutionary compensation, with a potential role for sequence conservation at the motif's vicinity. At the level of the full promoter, we identify cases of binding-site turnover, in which binding sites are reciprocally gained and lost, yet most interspecific differences remained uncompensated. Our results show the flexibility of TFs to bind imprecise motifs and the fast evolution of TF binding sites between related species.

Molecular and experimental tools to design synthetic enhancers

Understanding the grammar of enhancers and how they regulate gene expression is key for both basic research and for the pharma and biotech industries. The design and characterization of synthetic enhancers can expand the known regulatory space. This is achieved by the utilization of DNA Oligo Libraries (OLs), which facilitates screening of as many as millions of synthetic enhancer variants simultaneously. This review includes the latest commercial DNA OL synthesis technology and its capabilities, and a general 'know-how' guide for the design, construction, and analysis of OL-based synthetic enhancer characterization experiments. Specifically, we focus on synthetic-enhancer-based massively parallel reporter assay, Sort-seq methodologies (e.g. flow cytometry, deep sequencing), and a brief description of machine learning-based attempts for OL-analysis and follow-up validation experiments.

The microprotein Nrs1 rewires the G1/S transcriptional machinery during nitrogen limitation in budding yeast

Commitment to cell division at the end of G1 phase, termed Start in the budding yeast Saccharomyces cerevisiae, is strongly influenced by nutrient availability. To identify new dominant activators of Start that might operate under different nutrient conditions, we screened a genome-wide ORF overexpression library for genes that bypass a Start arrest caused by absence of the G1 cyclin Cln3 and the transcriptional activator Bck2. We recovered a hypothetical gene YLR053c, renamed NRS1 for Nitrogen-Responsive Start regulator 1, which encodes a poorly characterized 108 amino acid microprotein. Endogenous Nrs1 was nuclear-localized, restricted to poor nitrogen conditions, induced upon TORC1 inhibition, and cell cycle-regulated with a peak at Start. NRS1 interacted genetically with SWI4 and SWI6, which encode subunits of the main G1/S transcription factor complex SBF. Correspondingly, Nrs1 physically interacted with Swi4 and Swi6 and was localized to G1/S promoter DNA. Nrs1 exhibited inherent transactivation activity, and fusion of Nrs1 to the SBF inhibitor Whi5 was sufficient to suppress other Start defects. Nrs1 appears to be a recently evolved microprotein that rewires the G1/S transcriptional machinery under poor nitrogen conditions.

Commitment to cell division at the end of G1 phase in the budding yeast Saccharomyces cerevisiae is strongly influenced by nutrient availability. This study identifies a micro-protein that promotes G1/S transcription activation and cell cycle entry in yeast under nitrogen-limited conditions.

The evolution, evolvability and engineering of gene regulatory DNA

Mutations in non-coding regulatory DNA sequences can alter gene expression, organismal phenotype and fitness^1-3. Constructing complete fitness landscapes, in which DNA sequences are mapped to fitness, is a long-standing goal in biology, but has remained elusive because it is challenging to generalize reliably to vast sequence spaces^4-6. Here we build sequence-to-expression models that capture fitness landscapes and use them to decipher principles of regulatory evolution. Using millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Saccharomyces cerevisiae, we learn deep neural network models that generalize with excellent prediction performance, and enable sequence design for expression engineering. Using our models, we study expression divergence under genetic drift and strong-selection weak-mutation regimes to find that regulatory evolution is rapid and subject to diminishing returns epistasis; that conflicting expression objectives in different environments constrain expression adaptation; and that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for using such models to detect signatures of selection on expression from natural variation in regulatory sequences and use it to discover an instance of convergent regulatory evolution. We assess mutational robustness, finding that regulatory mutation effect sizes follow a power law, characterize regulatory evolvability, visualize promoter fitness landscapes, discover evolvability archetypes and illustrate the mutational robustness of natural regulatory sequence populations. Our work provides a general framework for designing regulatory sequences and addressing fundamental questions in regulatory evolution.

Loss of major nutrient sensing and signaling pathways suppresses starvation lethality in electron transport chain mutants

The electron transport chain (ETC) is a well-studied and highly conserved metabolic pathway that produces ATP through generation of a proton gradient across the inner mitochondrial membrane coupled to oxidative phosphorylation. ETC mutations are associated with a wide array of human disease conditions and to aging-related phenotypes in a number of different organisms. In this study, we sought to better understand the role of the ETC in aging using a yeast model. A panel of ETC mutant strains that fail to survive starvation was used to isolate suppressor mutants that survive. These suppressors tend to fall into major nutrient sensing and signaling pathways, suggesting that the ETC is involved in proper starvation signaling to these pathways in yeast. These suppressors also partially restore ETC-associated gene expression and pH homeostasis defects, though it remains unclear whether these phenotypes directly cause the suppression or are simply effects. This work further highlights the complex cellular network connections between metabolic pathways and signaling events in the cell and their potential roles in aging and age-related diseases.

Yeast Plasma Membrane Fungal Oligopeptide Transporters Display Distinct Substrate Preferences despite Their High Sequence Identity

Fungal Oligopeptide Transporters (Fot) Fot1, Fot2 and Fot3 have been found in Saccharomyces cerevisiae wine strains, but not in strains from other environments. In the S. cerevisiae wine strain EC1118, Fot1 and Fot2 are responsible for a broader range of oligopeptide utilization in comparison with strains not containing any Fot. This leads to better fermentation efficiency and an increased production of desirable organoleptic compounds in wine. Despite the benefits associated with Fot activity in S. cerevisiae within the wine environment, little is known about this family of transporters in yeast. The presence of Fot1, Fot2 and Fot3 in S. cerevisiae wine strains is due to horizontal gene transfer from the yeast Torulaspora microellipsoides, which harbors Fot2Tm, FotX and FotY proteins. Sequence analyses revealed that Fot family members have a high sequence identity in these yeast species. In this work, we aimed to further characterize the different Fot family members in terms of subcellular localization, gene expression in enological fermentation and substrate specificity. Using CRISPR/Cas9, we constructed S. cerevisiae wine strains containing each different Fot as the sole oligopeptide transporter to analyze their oligopeptide preferences by phenotype microarrays. The results of oligopeptide consumption show that Fot counterparts have different di-/tripeptide specificities, suggesting that punctual sequence divergence between FOT genes can be crucial for substrate recognition, binding and transport activity. FOT gene expression levels in different S. cerevisiae wine strains during enological fermentation, together with predicted binding motifs for transcriptional regulators in nitrogen metabolism, indicate that these transporters may be under the control of the Nitrogen Catabolite Repression (NCR) system. Finally, we demonstrated that Fot1 is located in the yeast plasma membrane. This work contributes to a better understanding of this family of oligopeptide transporters, which have demonstrated a key role in the utilization of oligopeptides by S. cerevisiae in enological fermentation.

Development of an Haa1-based biosensor for acetic acid sensing in Saccharomyces cerevisiae

Acetic acid is one of the main inhibitors of lignocellulosic hydrolysates and acetic acid tolerance is crucial for the development of robust cell factories for conversion of biomass. As a precursor of acetyl-coenzyme A, it also plays an important role in central carbon metabolism. Thus, monitoring acetic acid levels is a crucial aspect when cultivating yeast. Transcription factor-based biosensors represent useful tools to follow metabolite concentrations. Here, we present the development of an acetic acid biosensor based on the Saccharomyces cerevisiae transcription factor Haa1 that upon binding to acetic acid relocates to the nucleus. In the biosensor, a synthetic transcription factor consisting of Haa1 and BM3R1 from Bacillus megaterium was used to control expression of a reporter gene under a promoter containing BM3R1 binding sites. The biosensor did not drive expression under a promoter containing Haa1 binding sites and responded to acetic acid over a linear range spanning from 10 to 60 mM. To validate its applicability, the biosensor was integrated into acetic acid-producing strains. A direct correlation between biosensor output and acetic acid production was detected. The developed biosensor enables high-throughput screening of strains producing acetic acid and could also be used to investigate acetic acid-tolerant strain libraries.

A unique GCN5 histone acetyltransferase complex controls erythrocyte invasion and virulence in the malaria parasite Plasmodium falciparum

The histone acetyltransferase GCN5-associated SAGA complex is evolutionarily conserved from yeast to human and functions as a general transcription co-activator in global gene regulation. In this study, we identified a divergent GCN5 complex in Plasmodium falciparum, which contains two plant homeodomain (PHD) proteins (PfPHD1 and PfPHD2) and a plant apetela2 (AP2)-domain transcription factor (PfAP2-LT). To dissect the functions of the PfGCN5 complex, we generated parasite lines with either the bromodomain in PfGCN5 or the PHD domain in PfPHD1 deleted. The two deletion mutants closely phenocopied each other, exhibiting significantly reduced merozoite invasion of erythrocytes and elevated sexual conversion. These domain deletions caused dramatic decreases not only in histone H3K9 acetylation but also in H3K4 trimethylation, indicating synergistic crosstalk between the two euchromatin marks. Domain deletion in either PfGCN5 or PfPHD1 profoundly disturbed the global transcription pattern, causing altered expression of more than 60% of the genes. At the schizont stage, these domain deletions were linked to specific down-regulation of merozoite genes involved in erythrocyte invasion, many of which contain the AP2-LT binding motif and are also regulated by AP2-I and BDP1, suggesting targeted recruitment of the PfGCN5 complex to the invasion genes by these specific factors. Conversely, at the ring stage, PfGCN5 or PfPHD1 domain deletions disrupted the mutually exclusive expression pattern of the entire var gene family, which encodes the virulent factor PfEMP1. Correlation analysis between the chromatin state and alteration of gene expression demonstrated that up- and down-regulated genes in these mutants are highly correlated with the silent and active chromatin states in the wild-type parasite, respectively. Collectively, the PfGCN5 complex represents a novel HAT complex with a unique subunit composition including an AP2 transcription factor, which signifies a new paradigm for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Epigenetic regulation of gene expression plays essential roles in orchestrating the general and parasite-specific cellular pathways in the malaria parasite Plasmodium falciparum. To better understand the epigenetic mechanisms in this parasite, we characterized the histone acetyltransferase GCN5-mediated transcription regulation during intraerythrocytic development of the parasite. Using tandem affinity purification and proteomic characterization, we identified that the PfGCN5-associated complex contains nine core components, including two PHD domain proteins (PfPHD1 and PfPHD2) and an AP2-domain transcription factor, which is divergent from the canonical GCN5 complexes evolutionarily conserved from yeast to human. To understand the functions of the PfGCN5 complex, we performed domain deletions in two subunits of this complex, PfGCN5 and PfPHD1. We found that the two deletion mutants displayed very similar growth phenotypes, including significantly reduced merozoite invasion rates and elevated sexual conversion. These two mutants were associated with dramatic decreases in histone H3K9 acetylation and H3K4 trimethylation, which led to global changes in chromatin states and gene expression. Consistent with the phenotypes, genes significantly affected by the PfGCN5 and PfPHD1 gene disruption include those participating in parasite-specific pathways such as invasion, virulence, and sexual development. In conclusion, this study presents a new model of the PfGCN5 complex for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors

Activator proteins 1 (AP-1) comprise one of the largest families of eukaryotic basic leucine zipper transcription factors. Despite advances in the characterization of AP-1 DNA-binding sites, our ability to predict new binding sites and explain how the proteins achieve different gene expression levels remains limited. Here we address the role of sequence-specific DNA flexibility for stability and specific binding of AP-1 factors, using microsecond-long molecular dynamics simulations. As a model system, we employ yeast AP-1 factor Yap1 binding to three different response elements from two genetic environments. Our data show that Yap1 actively exploits the sequence-specific flexibility of DNA within the response element to form stable protein–DNA complexes. The stability also depends on the four to six flanking nucleotides, adjacent to the response elements. The flanking sequences modulate the conformational adaptability of the response element, making it more shape-efficient to form specific contacts with the protein. Bioinformatics analysis of differential expression of the studied genes supports our conclusions: the stability of Yap1–DNA complexes, modulated by the flanking environment, influences the gene expression levels. Our results provide new insights into mechanisms of protein–DNA recognition and the biological regulation of gene expression levels in eukaryotes.

HSF1 induces RNA polymerase II synthesis of ribosomal RNA in S. cerevisiae during nitrogen deprivation

The resource intensive process of accurate ribosome synthesis is essential for cell viability in all organisms. Ribosome synthesis regulation centers on RNA polymerase I (pol I) transcription of a 35S rRNA precursor that is processed into the mature 18S, 5.8S and 25S rRNAs. During nutrient deprivation or stress, pol I synthesis of rRNA is dramatically reduced. Conversely, chronic stress such as mitochondrial dysfunction induces RNA polymerase II (pol II) to transcribe functional rRNA using an evolutionarily conserved cryptic pol II rDNA promoter suggesting a universal phenomenon. However, this polymerase switches and its role in regulation of rRNA synthesis remain unclear. In this paper, we demonstrate that extended nitrogen deprivation induces the polymerase switch via components of the environmental stress response. We further show that the switch is repressed by Sch9 and activated by the stress kinase Rim15. Like stress-induced genes, the switch requires not only pol II transcription machinery, including the mediator, but also requires the HDAC, Rpd3 and stress transcription factor Hsf1. The current work shows that the constitutive allele, Hsf1^PO4* displays elevated levels of induction in non-stress conditions while binding to a conserved site in the pol II rDNA promoter upstream of the pol I promoter. Whether the polymerase switch serves to provide rRNA when pol I transcription is inhibited or fine-tunes pol I initiation via RNA interactions is yet to be determined. Identifying the underlying mechanism for this evolutionary conserved phenomenon will help understand the mechanism of pol II rRNA synthesis and its role in stress adaptation.

Advances in promoter engineering: novel applications and predefined transcriptional control

Synthetic biology continues to progress by relying on more robust tools for transcriptional control, of which promoters are the most fundamental component. Numerous studies have sought to characterize promoter function, determine principles to guide their engineering, and create promoters with stronger expression or tailored inducible control. In this review, we will summarize promoter architecture and highlight recent advances in the field, focusing on the novel applications of inducible promoter design and engineering towards metabolic engineering and cellular therapeutic development. Additionally, we will highlight how the expansion of new, machine learning techniques for modeling and engineering promoter sequences are enabling more accurate prediction of promoter characteristics. This article is protected by copyright. All rights reserved.

Hybrid-architectured promoter design to engineer expression in yeast

Engineered promoters are key components that allow engineered expression of genes in the cell-factory design. Promoters having exceptional strength are attractive candidates for designing metabolic engineering strategies for tailoring de novo production strategies that require directed evolution methods by engineering with de novo synthetic biology tools. Engineered promoter variants (EPVs) of naturally occurring promoters (NOPs) can be designed using metabolic engineering strategies and synthetic biology tools if the genes encoding the activating transcription factors (TFs) exist in the genome and are expressed and synthesized at non-limiting concentrations within the cell. The hybrid-architectured EPV design method targets an essential and predetermined part of the general transcription machinery. That is cis-acting DNA site(s) in coordination with the trans-acting factor(s) that must bind for the regulated transcription machinery activation. The method needs genomic and functional information that can lead to the discovery of the master TF(s) and synthetic cis-acting DNA elements, enabling the engineering of binding of master regulator TF(s). The method aims to generate EPVs that combine the advantages of being an exceptional stronger EPV(s) than the NOPs and permit "green-and-clean production" on a non-toxic carbon source, such as ethanol or glucose. By introducing our recent work on the engineering of ADH2 hybrid-promoter architectures to enhance recombinant protein expression on ethanol, we provide the method and protocol for the design of ADH2 hybrid-promoter architectures that can be adapted to other promoters in different substrate utilization pathways in Pichia pastoris (syn. Komagataella phaffii), as well as in other yeasts.

Transcription at a Distance in the Budding Yeast Saccharomyces cerevisiae

Proper transcriptional regulation depends on the collaboration of multiple layers of control simultaneously. Cells tightly balance cellular resources and integrate various signaling inputs to maintain homeostasis during growth, development and stressors, among other signals. Many eukaryotes, including the budding yeast Saccharomyces cerevisiae, exhibit a non-random distribution of functionally related genes throughout their genomes. This arrangement coordinates the transcription of genes that are found in clusters, and can occur over long distances. In this work, we review the current literature pertaining to gene regulation at a distance in budding yeast.

Yeast cell fate control by temporal redundancy modulation of transcription factor paralogs

Recent single-cell studies have revealed that yeast stress response involves transcription factors that are activated in pulses. However, it remains unclear whether and how these dynamic transcription factors temporally interact to regulate stress survival. Here we show that budding yeast cells can exploit the temporal relationship between paralogous general stress regulators, Msn2 and Msn4, during stress response. We find that individual pulses of Msn2 and Msn4 are largely redundant, and cells can enhance the expression of their shared targets by increasing their temporal divergence. Thus, functional redundancy between these two paralogs is modulated in a dynamic manner to confer fitness advantages for yeast cells, which might feed back to promote the preservation of their redundancy. This evolutionary implication is supported by evidence from Msn2/Msn4 orthologs and analyses of other transcription factor paralogs. Together, we show a cell fate control mechanism through temporal redundancy modulation in yeast, which may represent an evolutionarily important strategy for maintaining functional redundancy between gene duplicates.

A case of convergent-gene interference in the budding yeast knockout library causing chromosome instability

To maintain genome stability, organisms depend on faithful chromosome segregation, a process affected by diverse genetic pathways, some of which are not directly linked to mitosis. In this study, we set out to explore one such pathway represented by an undercharacterized gene, SNO1, identified previously in screens of the yeast knockout (YKO) library for mitotic fidelity genes. We found that the causative factor increasing mitotic error rate in the sno1Δ mutant is not loss of the Sno1 protein, but rather perturbation to the mRNA of the neighboring convergent gene, CTF13, encoding an essential component for forming the yeast kinetochore. This is caused by a combination of the Kanamycin resistance gene and the transcriptional terminator used in the YKO library affecting the CTF13 mRNA level and quality . We further provide a list of gene pairs potentially subjected to this artifact, which may be useful for accurate phenotypic interpretation of YKO mutants.

Hexose transport in Torulaspora delbrueckii: identification of Igt1, a new dual-affinity transporter

Torulaspora delbrueckii is a yeast species receiving increasing attention from the biotechnology industry, with particular relevance in the wine, beer and baking sectors. However, little is known about its sugar transporters and sugar transport capacity, frequently a rate-limiting step of sugar metabolism and efficient fermentation. Actually, only one glucose transporter, Lgt1, has been characterized so far. Here we report the identification and characterization of a second glucose transporter gene, IGT1, located in a cluster, upstream of LGT1 and downstream of two other putative hexose transporters. Functional characterization of IGT1 in a Saccharomyces cerevisiae hxt-null strain revealed that it encodes a transporter able to mediate uptake of glucose, fructose and mannose and established that its affinity, as measured by Km, could be modulated by glucose concentration in the medium. In fact, IGT1-transformed S. cerevisiae hxt-null cells, grown in 0.1% glucose displayed biphasic glucose uptake kinetics with an intermediate- (Km = 6.5 ± 2.0 mM) and a high-affinity (Km = 0.10 ± 0.01 mM) component, whereas cells grown in 2% glucose displayed monophasic kinetics with an intermediate-affinity (Km of 11.5 ± 1.5 mM). This work contributes to a better characterization of glucose transport in T. delbrueckii, with relevant implications for its exploitation in the food industry.

Toward the discovery of biological functions associated with the mechanosensor Mtl1p of Saccharomyces cerevisiae via integrative multi-OMICs analysis

Functional analysis of the Mtl1 protein in Saccharomyces cerevisiae has revealed that this transmembrane sensor endows yeast cells with resistance to oxidative stress through a signaling mechanism called the cell wall integrity pathway (CWI). We observed upregulation of multiple heat shock proteins (HSPs), proteins associated with the formation of stress granules, and the phosphatase subunit of trehalose 6-phosphate synthase which suggests that mtl1Δ strains undergo intrinsic activation of a non-lethal heat stress response. Furthermore, quantitative global proteomic analysis conducted on TMT-labeled proteins combined with metabolome analysis revealed that mtl1Δ strains exhibit decreased levels of metabolites of carboxylic acid metabolism, decreased expression of anabolic enzymes and increased expression of catabolic enzymes involved in the metabolism of amino acids, with enhanced expression of mitochondrial respirasome proteins. These observations support the idea that Mtl1 protein controls the suppression of a non-lethal heat stress response under normal conditions while it plays an important role in metabolic regulatory mechanisms linked to TORC1 signaling that are required to maintain cellular homeostasis and optimal mitochondrial function.

A Reporter System for Cytosolic Protein Aggregates in Yeast

Protein misfolding and aggregation are linked to neurodegenerative diseases of mammals and suboptimal protein expression within biotechnology. Tools for monitoring protein aggregates are therefore useful for studying disease-related aggregation and for improving soluble protein expression in heterologous hosts for biotechnology purposes. In this work, we developed a promoter-reporter system for aggregated protein on the basis of the yeast native response to misfolded protein. To this end, we first studied the proteome of yeast in response to the expression of folded soluble and aggregation-prone protein baits and identified genes encoding proteins related to protein folding and the response to heat stress as well as the ubiquitin-proteasome system that are over-represented in cells expressing an aggregation-prone protein. From these data, we created and validated promoter-reporter constructs and further engineered the best performing promoters by increasing the copy number of upstream activating sequences and optimization of culture conditions. Our best promoter-reporter has an output dynamic range of approximately 12-fold upon expression of the aggregation-prone protein and responded to increasing levels of aggregated protein. Finally, we demonstrate that the system can discriminate between yeast cells expressing different prion precursor proteins and select the cells expressing folded soluble protein from mixed populations. Our reporter system is thus a simple tool for diagnosing protein aggregates in living cells and should be applicable for the health and biotechnology industries.

CRISPRi screens reveal genes modulating yeast growth in lignocellulose hydrolysate

BACKGROUND

Baker's yeast is a widely used eukaryotic cell factory, producing a diverse range of compounds including biofuels and fine chemicals. The use of lignocellulose as feedstock offers the opportunity to run these processes in an environmentally sustainable way. However, the required hydrolysis pretreatment of lignocellulosic material releases toxic compounds that hamper yeast growth and consequently productivity.

RESULTS

Here, we employ CRISPR interference in S. cerevisiae to identify genes modulating fermentative growth in plant hydrolysate and in presence of lignocellulosic toxins. We find that at least one-third of hydrolysate-associated gene functions are explained by effects of known toxic compounds, such as the decreased growth of YAP1 or HAA1, or increased growth of DOT6 knock-down strains in hydrolysate.

CONCLUSION

Our study confirms previously known genetic elements and uncovers new targets towards designing more robust yeast strains for the utilization of lignocellulose hydrolysate as sustainable feedstock, and, more broadly, paves the way for applying CRISPRi screens to improve industrial fermentation processes.

Accumulation of cis- and trans-regulatory variations is associated with phenotypic divergence of a complex trait between yeast species

Gene regulatory variations accumulate during evolution and alter gene expression. While the importance of expression variation in phenotypic evolution is well established, the molecular basis remains largely unknown. Here, we examine two closely related yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus, which show phenotypical differences in morphology and cell cycle progression when grown in the same environment. By profiling the cell cycle transcriptome and binding of key transcription factors (TFs) in the two species and their hybrid, we show that changes in expression levels and dynamics of oscillating genes are dominated by upstream trans-variations. We find that multiple cell cycle regulators show both cis- and trans-regulatory variations, which alters their expression in favor of the different cell cycle phenotypes. Moreover, we show that variations in the cell cycle TFs, Fkh1, and Fkh2 affect both the expression of target genes, and the binding specificity of an interacting TF, Ace2. Our study reveals how multiple variations accumulate and propagate through the gene regulatory network, alter TFs binding, contributing to phenotypic changes in cell cycle progression.

Learning causal networks using inducible transcription factors and transcriptome-wide time series

We present IDEA (the Induction Dynamics gene Expression Atlas), a dataset constructed by independently inducing hundreds of transcription factors (TFs) and measuring timecourses of the resulting gene expression responses in budding yeast. Each experiment captures a regulatory cascade connecting a single induced regulator to the genes it causally regulates. We discuss the regulatory cascade of a single TF, Aft1, in detail; however, IDEA contains > 200 TF induction experiments with 20 million individual observations and 100,000 signal‐containing dynamic responses. As an application of IDEA, we integrate all timecourses into a whole‐cell transcriptional model, which is used to predict and validate multiple new and underappreciated transcriptional regulators. We also find that the magnitudes of coefficients in this model are predictive of genetic interaction profile similarities. In addition to being a resource for exploring regulatory connectivity between TFs and their target genes, our modeling approach shows that combining rapid perturbations of individual genes with genome‐scale time‐series measurements is an effective strategy for elucidating gene regulatory networks.

Speed-Specificity Trade-Offs in the Transcription Factors Search for Their Genomic Binding Sites

Transcription factors (TFs) regulate gene expression by binding DNA sequences recognized by their DNA-binding domains (DBDs). DBD-recognized motifs are short and highly abundant in genomes. The ability of TFs to bind a specific subset of motif-containing sites, and to do so rapidly upon activation, is fundamental for gene expression in all eukaryotes. Despite extensive interest, our understanding of the TF-target search process is fragmented; although binding specificity and detection speed are two facets of this same process, trade-offs between them are rarely addressed. In this opinion article, we discuss potential speed-specificity trade-offs in the context of existing models. We further discuss the recently described 'distributed specificity' paradigm, suggesting that intrinsically disordered regions (IDRs) promote specificity while reducing the TF-target search time.

Genome-wide off-rates reveal how DNA binding dynamics shape transcription factor function

Protein-DNA interactions are dynamic, and these dynamics are an important aspect of chromatin-associated processes such as transcription or replication. Due to a lack of methods to study on- and off-rates across entire genomes, protein-DNA interaction dynamics have not been studied extensively. Here, we determine in vivo off-rates for the Saccharomyces cerevisiae chromatin organizing factor Abf1, at 191 sites simultaneously across the yeast genome. Average Abf1 residence times span a wide range, varying between 4.2 and 33 min. Sites with different off-rates are associated with different functional characteristics. This includes their transcriptional dependency on Abf1, nucleosome positioning and the size of the nucleosome-free region, as well as the ability to roadblock RNA polymerase II for termination. The results show how off-rates contribute to transcription factor function and that DIVORSEQ (Determining In Vivo Off-Rates by SEQuencing) is a meaningful way of investigating protein-DNA binding dynamics genome-wide.

A Step-by-Step Protocol for COMPASS, a Synthetic Biology Tool for Combinatorial Gene Assembly

For industry-scale production of high-value chemicals in microbial cell factories, the elimination of metabolic flux imbalances is a critical aspect. However, a priori knowledge about the genetic design of optimal production pathways is typically not available. COMPASS, COMbinatorial Pathway ASSembly, is a rapid cloning method for the balanced expression of multiple genes in biochemical pathways. The method generates thousands of individual DNA constructs in modular, parallel, and high-throughput manner. COMPASS employs inducible artificial transcription factors derived from plant (Arabidopsis thaliana) regulators to control the expression of pathway genes in yeast (Saccharomyces cerevisiae). It utilizes homologous recombination for parts assembly and employs a positive selection scheme to identify correctly assembled pathway variants after both in vivo and in vitro recombination. Finally, COMPASS is equipped with a CRISPR/Cas9 genome modification system allowing for the one-step multilocus integration of genes. Although COMPASS was initially developed for pathway engineering, it can equally be employed for balancing gene expression in other synthetic biology projects.

Transient intracellular acidification regulates the core transcriptional heat shock response

Heat shock induces a conserved transcriptional program regulated by heat shock factor 1 (Hsf1) in eukaryotic cells. Activation of this heat shock response is triggered by heat-induced misfolding of newly synthesized polypeptides, and so has been thought to depend on ongoing protein synthesis. Here, using the budding yeast Saccharomyces cerevisiae, we report the discovery that Hsf1 can be robustly activated when protein synthesis is inhibited, so long as cells undergo cytosolic acidification. Heat shock has long been known to cause transient intracellular acidification which, for reasons which have remained unclear, is associated with increased stress resistance in eukaryotes. We demonstrate that acidification is required for heat shock response induction in translationally inhibited cells, and specifically affects Hsf1 activation. Physiological heat-triggered acidification also increases population fitness and promotes cell cycle reentry following heat shock. Our results uncover a previously unknown adaptive dimension of the well-studied eukaryotic heat shock response.

An expanded cell wall damage signaling network is comprised of the transcription factors Rlm1 and Sko1 in Candida albicans

The human fungal pathogen Candida albicans is constantly exposed to environmental challenges impacting the cell wall. Signaling pathways coordinate stress adaptation and are essential for commensalism and virulence. The transcription factors Sko1, Cas5, and Rlm1 control the response to cell wall stress caused by the antifungal drug caspofungin. Here, we expand the Sko1 and Rlm1 transcriptional circuit and demonstrate that Rlm1 activates Sko1 cell wall stress signaling. Caspofungin-induced transcription of SKO1 and several Sko1-dependent cell wall integrity genes are attenuated in an rlm1Δ/Δ mutant strain when compared to the treated wild-type strain but not in a cas5Δ/Δ mutant strain. Genome-wide chromatin immunoprecipitation (ChIP-seq) results revealed numerous Sko1 and Rlm1 directly bound target genes in the presence of caspofungin that were undetected in previous gene expression studies. Notable targets include genes involved in cell wall integrity, osmolarity, and cellular aggregation, as well as several uncharacterized genes. Interestingly, we found that Rlm1 does not bind to the upstream intergenic region of SKO1 in the presence of caspofungin, indicating that Rlm1 indirectly controls caspofungin-induced SKO1 transcription. In addition, we discovered that caspofungin-induced SKO1 transcription occurs through self-activation. Based on our ChIP-seq data, we also discovered an Rlm1 consensus motif unique to C. albicans. For Sko1, we found a consensus motif similar to the known Sko1 motif for Saccharomyces cerevisiae. Growth assays showed that SKO1 overexpression suppressed caspofungin hypersensitivity in an rlm1Δ/Δ mutant strain. In addition, overexpression of the glycerol phosphatase, RHR2, suppressed caspofungin hypersensitivity specifically in a sko1Δ/Δ mutant strain. Our findings link the Sko1 and Rlm1 signaling pathways, identify new biological roles for Sko1 and Rlm1, and highlight the complex dynamics underlying cell wall signaling.

Candida albicans is the most common human fungal pathogen isolated in clinical settings. The echinocandin drug caspofungin is used to treat invasive candidiasis; however, the emergence of increasing echinocandin resistance underscores the need for new antifungal strategies. Elucidating the signaling mechanisms that govern caspofungin-induced tolerance has the potential to identify candidate proteins that could serve as novel therapeutic targets. Here, we expand the Rlm1 and Sko1 cell wall transcriptional network and find that Rlm1 indirectly regulates Sko1 signaling. Furthermore, we identify Sko1- and Rlm1-specific biological roles in caspofungin adaptation, such as osmoregulation and secretion. Lastly, we discover a protective role for glycerol in caspofungin tolerance. Overall, these findings provide mechanistic insight into the genetic and cellular bases underlying cell wall signaling in C. albicans.

Independent evolution of transcript abundance and gene regulatory dynamics

Changes in gene expression drive novel phenotypes, raising interest in how gene expression evolves. In contrast to the static genome, cells modulate gene expression in response to changing environments. Previous comparative studies focused on specific conditions, describing interspecies variation in expression levels, but providing limited information about variation across different conditions. To close this gap, we profiled mRNA levels of two related yeast species in hundreds of conditions and used coexpression analysis to distinguish variation in the dynamic pattern of gene expression from variation in expression levels. The majority of genes whose expression varied between the species maintained a conserved dynamic pattern. Cases of diverged dynamic pattern correspond to genes that were induced under distinct subsets of conditions in the two species. Profiling the interspecific hybrid allowed us to distinguish between genes with predominantly cis- or trans-regulatory variation. We find that trans-varying alleles are dominantly inherited, and that cis-variations are often complemented by variations in trans Based on these results, we suggest that gene expression diverges primarily through changes in expression levels, but does not alter the pattern by which these levels are dynamically regulated.