Basic leucine zipper transcription factor Hac1 binds DNA in two distinct modes as revealed by microfluidic analyses

Evolutionary-Scale Enzymology Enables Biochemical Constant Prediction Across a Multi-Peaked Catalytic Landscape

Quantitatively mapping enzyme sequence-catalysis landscapes remains a critical challenge in understanding enzyme function, evolution, and design. Here, we expand an emerging microfluidic platform to measure catalytic constants- k cat and K M -for hundreds of diverse naturally occurring sequences and mutants of the model enzyme Adenylate Kinase (ADK). This enables us to dissect the sequence-catalysis landscape's topology, navigability, and mechanistic underpinnings, revealing distinct catalytic peaks organized by structural motifs. These results challenge long-standing hypotheses in enzyme adaptation, demonstrating that thermophilic enzymes are not slower than their mesophilic counterparts. Combining the rich representations of protein sequences provided by deep-learning models with our custom high-throughput kinetic data yields semi-supervised models that significantly outperform existing models at predicting catalytic parameters of naturally occurring ADK sequences. Our work demonstrates a promising strategy for dissecting sequence-catalysis landscapes across enzymatic evolution and building family-specific models capable of accurately predicting catalytic constants, opening new avenues for enzyme engineering and functional prediction.

HacA, a key transcription factor for the unfolded protein response, is required for fungal development, aflatoxin biosynthesis and pathogenicity of Aspergillus flavus

Aspergillus flavus is a fungus notorious for contaminating food and feed with aflatoxins. As a saprophytic fungus, it secretes large amounts of enzymes to access nutrients, making endoplasmic reticulum (ER) homeostasis important for protein folding and secretion. The role of HacA, a key transcription factor in the unfolded protein response pathway, remains poorly understood in A. flavus. In this study, the hacA gene in A. flavus was knockout. Results showed that the absence of hacA led to a decreased pathogenicity of the strain, as it failed to colonize intact maize kernels. This may be due to retarded vegetable growth, especially the abnormal development of swollen tips and shorter hyphal septa. Deletion of hacA also hindered conidiogenesis and sclerotial development. Notably, the mutant strain failed to produce aflatoxin B1. Moreover, compared to the wild type, the mutant strain showed increased sensitivity to ER stress inducer such as Dithiothreitol (DTT), and heat stress. It also displayed heightened sensitivity to other environmental stresses, including cell wall, osmotic, and pH stresses. Further transcriptomic analysis revealed the involvement of the hacA in numerous biological processes, including filamentous growth, asexual reproduction, mycotoxin biosynthetic process, signal transduction, budding cell apical bud growth, invasive filamentous growth, response to stimulus, and so on. Taken together, HacA plays a vital role in fungal development, pathogenicity and aflatoxins biosynthesis. This highlights the potential of targeting hacA as a novel approach for early prevention of A. flavus contamination.

Tailored UPRE2 variants for dynamic gene regulation in yeast

Genetic elements are foundational in synthetic biology serving as vital building blocks. They enable programming host cells for efficient production of valuable chemicals and recombinant proteins. The unfolded protein response (UPR) is a stress pathway in which the transcription factor Hac1 interacts with the upstream unfolded protein response element (UPRE) of the promoter to restore endoplasmic reticulum (ER) homeostasis. Here, we created a UPRE2 mutant (UPRE2m) library. Several rounds of screening identified many elements with enhanced responsiveness and a wider dynamic range. The most active element m84 displayed a response activity 3.72 times higher than the native UPRE2. These potent elements are versatile and compatible with various promoters. Overexpression of HAC1 enhanced stress signal transduction, expanding the signal output range of UPRE2m. Through molecular modeling and site-directed mutagenesis, we pinpointed the DNA-binding residue Lys60 in Hac1(Hac1-K60). We also confirmed that UPRE2m exhibited a higher binding affinity to Hac1. This shed light on the mechanism underlying the Hac1-UPRE2m interaction. Importantly, applying UPRE2m for target gene regulation effectively increased both recombinant protein production and natural product synthesis. These genetic elements provide valuable tools for dynamically regulating gene expression in yeast cell factories.

Amoeba predation of Cryptococcus: A quantitative and population genomic evaluation of the accidental pathogen hypothesis

The "Amoeboid Predator-Fungal Animal Virulence Hypothesis" posits that interactions with environmental phagocytes shape the evolution of virulence traits in fungal pathogens. In this hypothesis, selection to avoid predation by amoeba inadvertently selects for traits that contribute to fungal escape from phagocytic immune cells. Here, we investigate this hypothesis in the human fungal pathogens Cryptococcus neoformans and Cryptococcus deneoformans. Applying quantitative trait locus (QTL) mapping and comparative genomics, we discovered a cross-species QTL region that is responsible for variation in resistance to amoeba predation. In C. neoformans, this same QTL was found to have pleiotropic effects on melanization, an established virulence factor. Through fine mapping and population genomic comparisons, we identified the gene encoding the transcription factor Bzp4 that underlies this pleiotropic QTL and we show that decreased expression of this gene reduces melanization and increases susceptibility to amoeba predation. Despite the joint effects of BZP4 on amoeba resistance and melanin production, we find no relationship between BZP4 genotype and escape from macrophages or virulence in murine models of disease. Our findings provide new perspectives on how microbial ecology shapes the genetic architecture of fungal virulence, and suggests the need for more nuanced models for the evolution of pathogenesis that account for the complexities of both microbe-microbe and microbe-host interactions.

Fundamental and Applicative Aspects of the Unfolded Protein Response in Yeasts

Upon the dysfunction or functional shortage of the endoplasmic reticulum (ER), namely, ER stress, eukaryotic cells commonly provoke a protective gene expression program called the unfolded protein response (UPR). The molecular mechanism of UPR has been uncovered through frontier genetic studies using Saccharomyces cerevisiae as a model organism. Ire1 is an ER-located transmembrane protein that directly senses ER stress and is activated as an RNase. During ER stress, Ire1 promotes the splicing of HAC1 mRNA, which is then translated into a transcription factor that induces the expression of various genes, including those encoding ER-located molecular chaperones and protein modification enzymes. While this mainstream intracellular UPR signaling pathway was elucidated in the 1990s, new intriguing insights have been gained up to now. For instance, various additional factors allow UPR evocation strictly in response to ER stress. The UPR machineries in other yeasts and fungi, including pathogenic species, are another important research topic. Moreover, industrially beneficial yeast strains carrying an enforced and enlarged ER have been produced through the artificial and constitutive induction of the UPR. In this article, we review canonical and up-to-date insights concerning the yeast UPR, mainly from the viewpoint of the functions and regulation of Ire1 and HAC1.

Overexpression of genes by stress-responsive promoters increases protein secretion in Saccharomyces cerevisiae

Recombinant proteins produced by cell factories are now widely used in various fields. Many efforts have been made to improve the secretion capacity of cell factories to meet the increasing demand for recombinant proteins. Recombinant protein production usually causes cell stress in the endoplasmic reticulum (ER). The overexpression of key genes possibly removes limitations in protein secretion. However, inappropriate gene expression may have negative effects. There is a need for dynamic control of genes adapted to cellular status. In this study, we constructed and characterized synthetic promoters that were inducible under ER stress conditions in Saccharomyces cerevisiae. The unfolded protein response element UPRE2, responding to stress with a wide dynamic range, was assembled with various promoter core regions, resulting in UPR-responsive promoters. Synthetic responsive promoters regulated gene expression by responding to stress level, which reflected the cellular status. The engineered strain using synthetic responsive promoters P_{4UPRE2 - TDH3} and P_{4UPRE2 - TEF1} for co-expression of ERO1 and SLY1 had 95% higher α-amylase production compared with the strain using the native promoters P_TDH3 and P_TEF1. This work showed that UPR-responsive promoters were useful in the metabolic engineering of yeast strains for tuning genes to support efficient protein production.

Meiotic resetting of the cellular Sod1 pool is driven by protein aggregation, degradation, and transient LUTI-mediated repression

Gametogenesis requires packaging of the cellular components needed for the next generation. In budding yeast, this process includes degradation of many mitotically stable proteins, followed by their resynthesis. Here, we show that one such case-Superoxide dismutase 1 (Sod1), a protein that commonly aggregates in human ALS patients-is regulated by an integrated set of events, beginning with the formation of pre-meiotic Sod1 aggregates. This is followed by degradation of a subset of the prior Sod1 pool and clearance of Sod1 aggregates. As degradation progresses, Sod1 protein production is transiently blocked during mid-meiotic stages by transcription of an extended and poorly translated SOD1 mRNA isoform, SOD1LUTI. Expression of SOD1LUTI is induced by the Unfolded Protein Response, and it acts to repress canonical SOD1 mRNA expression. SOD1LUTI is no longer expressed following the meiotic divisions, enabling a resurgence of canonical mRNA and synthesis of new Sod1 protein such that gametes inherit a full complement of Sod1 protein. Failure to aggregate and degrade Sod1 results in reduced gamete fitness in the presence of oxidants, highlighting the importance of this regulation. Investigation of Sod1 during yeast gametogenesis, an unusual cellular context in which Sod1 levels are tightly regulated, could shed light on conserved aspects of its aggregation and degradation, with relevance to understanding Sod1's role in human disease.

Fast-Growing Saccharomyces cerevisiae Cells with a Constitutive Unfolded Protein Response and Their Potential for Lipidic Molecule Production

In Saccharomyces cerevisiae cells, dysfunction of the endoplasmic reticulum (ER), so-called ER stress, leads to conversion of HAC1 mRNA to the spliced form (HAC1i), which is translated into a transcription factor that drastically changes the gene expression profile. This cellular response ultimately enhances ER functions and is named the unfolded protein response (UPR). Artificial evocation of the UPR is therefore anticipated to increase productivity of beneficial materials on and in the ER. However, as demonstrated here, cells constitutively expressing HAC1i mRNA (HAC1i cells), which exhibited a strong UPR even under nonstress conditions, grew considerably slowly and frequently yielded fast-growing and low-UPR progeny. Intriguingly, growth of HAC1i cells was faster in the presence of weak ER stress that was induced by low concentrations of the ER stressor tunicamycin or by cellular expression of the ER-located version of green fluorescent protein (GFP). HAC1i cells producing ER-localized GFP stably exhibited a strong UPR activity, carried a highly expanded ER, and abundantly produced triglycerides and heterogenous carotenoids. We therefore propose that our findings provide a basis for metabolic engineering to generate cells producing valuable lipidic molecules. IMPORTANCE The UPR is thought to be a cellular response to cope with the accumulation of unfolded proteins in the ER. In S. cerevisiae cells, the UPR is severely repressed under nonstress conditions. The findings of this study shed light on the physiological significance of the tight regulation of the UPR. Constitutive UPR induction caused considerable growth retardation, which was partly rescued by the induction of weak ER stress. Therefore, we speculate that when the UPR is inappropriately induced in unstressed cells lacking aberrant ER client proteins, the UPR improperly impairs normal cellular functions. Another important point of this study was the generation of S. cerevisiae strains stably exhibiting a strong UPR activity and abundantly producing triglycerides and heterogenous carotenoids. We anticipate that our findings may be applied to produce valuable lipidic molecules using yeast cells as a potential next-generation technique.

Characterization and modulation of endoplasmic reticulum stress response target genes in Kluyveromyces marxianus to improve secretory expressions of heterologous proteins

BACKGROUND

Kluyveromyces marxianus is a promising cell factory for producing bioethanol and that raised a demand for a high yield of heterologous proteins in this species. Expressions of heterologous proteins usually lead to the accumulation of misfolded or unfolded proteins in the lumen of the endoplasmic reticulum (ER) and then cause ER stress. To cope with this problem, a group of ER stress response target genes (ESRTs) are induced, mainly through a signaling network called unfolded protein response (UPR). Characterization and modulation of ESRTs direct the optimization of heterologous expressions. However, ESRTs in K. marxianus have not been identified so far.

RESULTS

In this study, we characterized the ER stress response in K. marxianus for the first time, by using two ER stress-inducing reagents, dithiothreitol (DTT) and tunicamycin (TM). Results showed that the Kar2-Ire1-Hac1 pathway of UPR is well conserved in K. marxianus. About 15% and 6% of genes were upregulated during treatment of DTT and TM, respectively. A total of 115 upregulated genes were characterized as ESRTs, among which 97 genes were identified as UPR target genes and 37 UPR target genes contained UPR elements in their promoters. Genes related to carbohydrate metabolic process and actin filament organization were identified as new types of UPR target genes. A total of 102 ESRTs were overexpressed separately in plasmids and their effects on productions of two different lignocellulolytic enzymes were systematically evaluated. Overexpressing genes involved in carbohydrate metabolism, including PDC1, PGK and VID28, overexpressing a chaperone gene CAJ1 or overexpressing a reductase gene MET13 substantially improved secretion expressions of heterologous proteins. Meanwhile, overexpressing a novel gene, KLMA_50479 (named ESR1), as well as overexpressing genes involved in ER-associated protein degradation (ERAD), including HRD3, USA1 andYET3, reduced the secretory expressions. ESR1 and the aforementioned ERAD genes were deleted from the genome. Resultant mutants, except the yet3Δ mutant, substantially improved secretions of three different heterologous proteins. During the fed-batch fermentation, extracellular activities of an endoxylanase and a glucanase in hrd3Δ cells improved by 43% and 28%, respectively, compared to those in wild-type cells.

CONCLUSIONS

Our results unveil the transcriptional scope of the ER stress response in K. marxianus and suggest efficient ways to improve productions of heterologous proteins by manipulating expressions of ESRTs.

The unfolded protein response in Pichia pastoris without external stressing stimuli

Dysfunction or capacity shortage of the endoplasmic reticulum (ER) is cumulatively called ER stress and provokes the unfolded protein response (UPR). In various yeast species, the ER-located transmembrane protein Ire1 is activated upon ER stress and performs the splicing reaction of HAC1 mRNA, the mature form of which is translated into a transcription factor protein that is responsible for the transcriptome change on the UPR. Here we carefully assessed the splicing of HAC1 mRNA in Pichia pastoris (Komagataella phaffii) cells. We found that, inconsistent with previous reports by others, the HAC1 mRNA was substantially, but partially, spliced even without ER-stressing stimuli. Unlike Saccharomyces cerevisiae, growth of P. pastoris was significantly retarded by the IRE1-gene knockout mutation. Moreover, P. pastoris cells seemed to push more abundant proteins into the secretory pathway than S. cerevisiae cells. We also suggest that P. pastoris Ire1 has the ability to control its activity stringently in an ER stress-dependent manner. We thus propose that P. pastoris cells are highly ER-stressed possibly because of the high load of endogenous proteins into the ER.

Revealing enzyme functional architecture via high-throughput microfluidic enzyme kinetics

Systematic and extensive investigation of enzymes is needed to understand their extraordinary efficiency and meet current challenges in medicine and engineering. We present HT-MEK (High-Throughput Microfluidic Enzyme Kinetics), a microfluidic platform for high-throughput expression, purification, and characterization of more than 1500 enzyme variants per experiment. For 1036 mutants of the alkaline phosphatase PafA (phosphate-irrepressible alkaline phosphatase of Flavobacterium), we performed more than 670,000 reactions and determined more than 5000 kinetic and physical constants for multiple substrates and inhibitors. We uncovered extensive kinetic partitioning to a misfolded state and isolated catalytic effects, revealing spatially contiguous regions of residues linked to particular aspects of function. Regions included active-site proximal residues but extended to the enzyme surface, providing a map of underlying architecture not possible to derive from existing approaches. HT-MEK has applications that range from understanding molecular mechanisms to medicine, engineering, and design.

Yeast Synthetic Minimal Biosensors for Evaluating Protein Production

The unfolded protein response (UPR) is a highly conserved cellular response in eukaryotic cells to counteract endoplasmic reticulum (ER) stress, typically triggered by unfolded protein accumulation. In addition to its relevance to human diseases like cancer, the induction of the UPR has a significant impact on the recombinant protein production in eukaryotic cell factories, including the industrial workhorseSaccharomyces cerevisiae. Being able to accurately detect and measure this ER stress response in single cells, enables the rapid optimization of protein production conditions and high-throughput strain selection strategies. Current methodologies to monitor the UPR in S. cerevisiae are often temporally and spatially removed from the cultivation stage or lack updated systematic evaluation. To this end, we constructed and systematically evaluated a series of high-throughput UPR sensors by different designs, incorporating either yeast native UPR promoters or novel synthetic minimal UPR promoters. The native promoters of DER1 and ERO1 were identified to have suitable UPR biosensor properties and served as an expression level guide for orthogonal sensor benchmarking. Our best synthetic minimal sensor is only 98 bp in length, has minimal homology to other native yeast sequences and displayed superior sensor characteristics. The synthetic minimal UPR sensor was able to accurately distinguish between cells expressing different heterologous proteins and between the different secretion levels of the same protein. This work demonstrated the potential of synthetic UPR biosensors as high-throughput tools to predict the protein production capacity of strains, interrogate protein properties hampering their secretion, and guide rational engineering strategies for optimal heterologous protein production.

The Yeast eIF2 Kinase Gcn2 Facilitates H2O2-Mediated Feedback Inhibition of Both Protein Synthesis and Endoplasmic Reticulum Oxidative Folding during Recombinant Protein Production

Recombinant protein production is a known source of oxidative stress. However, knowledge of which reactive oxygen species are involved or the specific growth phase in which stress occurs remains lacking. Using modern, hypersensitive genetic H₂O₂-specific probes, microcultivation, and continuous measurements in batch culture, we observed H₂O₂ accumulation during and following the diauxic shift in engineered Saccharomyces cerevisiae, correlating with peak α-amylase production. In agreement with previous studies supporting a role of the translation initiation factor kinase Gcn2 in the response to H₂O₂, we find that Gcn2-dependent phosphorylation of eIF2α increases alongside translational attenuation in strains engineered to produce large amounts of α-amylase. Gcn2 removal significantly improved α-amylase production in two previously optimized high-producing strains but not in the wild type. Gcn2 deficiency furthermore reduced intracellular H₂O₂ levels and the Hac1 splicing ratio, while expression of antioxidants and the endoplasmic reticulum (ER) disulfide isomerase PDI1 increased. These results suggest protein synthesis and ER oxidative folding are coupled and subject to feedback inhibition by H₂O₂.

IMPORTANCE Recombinant protein production is a multibillion dollar industry. Optimizing the productivity of host cells is, therefore, of great interest. In several hosts, oxidants are produced as an unwanted side product of recombinant protein production. The buildup of oxidants can result in intracellular stress responses that could compromise the productivity of the host cell. Here, we document a novel protein synthesis inhibitory mechanism that is activated by the buildup of a specific oxidant (H₂O₂) in the cytosol of yeast cells upon the production of recombinant proteins. At the center of this inhibitory mechanism lies the protein kinase Gcn2. By removing Gcn2, we observed a doubling of recombinant protein productivity in addition to reduced H₂O₂ levels in the cytosol. In this study, we want to raise awareness of this inhibitory mechanism in eukaryotic cells to further improve protein production and contribute to the development of novel protein-based therapeutic strategies.

Bayesian Markov models improve the prediction of binding motifs beyond first order

Transcription factors (TFs) regulate gene expression by binding to specific DNA motifs. Accurate models for predicting binding affinities are crucial for quantitatively understanding of transcriptional regulation. Motifs are commonly described by position weight matrices, which assume that each position contributes independently to the binding energy. Models that can learn dependencies between positions, for instance, induced by DNA structure preferences, have yielded markedly improved predictions for most TFs on in vivo data. However, they are more prone to overfit the data and to learn patterns merely correlated with rather than directly involved in TF binding. We present an improved, faster version of our Bayesian Markov model software, BaMMmotif2. We tested it with state-of-the-art motif discovery tools on a large collection of ChIP-seq and HT-SELEX datasets. BaMMmotif2 models of fifth-order achieved a median false-discovery-rate-averaged recall 13.6% and 12.2% higher than the next best tool on 427 ChIP-seq datasets and 164 HT-SELEX datasets, respectively, while being 8 to 1000 times faster. BaMMmotif2 models showed no signs of overtraining in cross-cell line and cross-platform tests, with similar improvements on the next-best tool. These results demonstrate that dependencies beyond first order clearly improve binding models for most TFs.

Does Saccharomyces cerevisiae Require Specific Post-Translational Silencing against Leaky Translation of Hac1up?

HAC1 encodes a key transcription factor that transmits the unfolded protein response (UPR) from the endoplasmic reticulum (ER) to the nucleus and regulates downstream UPR genes in Saccharomyces cerevisiae. In response to the accumulation of unfolded proteins in the ER, Ire1p oligomers splice HAC1 pre-mRNA (HAC1^u) via a non-conventional process and allow the spliced HAC1 (HAC1ⁱ) to be translated efficiently. However, leaky splicing and translation of HAC1^u may occur in non-UPR cells to induce undesirable UPR. To control accidental UPR activation, multiple fail-safe mechanisms have been proposed to prevent leaky HAC1 splicing and translation and to facilitate rapid degradation of translated Hac1^up and Hac1ⁱp. Among proposed regulatory mechanisms is a degron sequence encoded at the 5′ end of the HAC1 intron that silences Hac1^up expression. To investigate the necessity of an intron-encoded degron sequence that specifically targets Hac1^up for degradation, we employed publicly available transcriptomic data to quantify leaky HAC1 splicing and translation in UPR-induced and non-UPR cells. As expected, we found that HAC1^u is only efficiently spliced into HAC1ⁱ and efficiently translated into Hac1ⁱp in UPR-induced cells. However, our analysis of ribosome profiling data confirmed frequent occurrence of leaky translation of HAC1^u regardless of UPR induction, demonstrating the inability of translation fail-safe to completely inhibit Hac1^up production. Additionally, among 32 yeast HAC1 surveyed, the degron sequence is highly conserved by Saccharomyces yeast but is poorly conserved by all other yeast species. Nevertheless, the degron sequence is the most conserved HAC1 intron segment in yeasts. These results suggest that the degron sequence may indeed play an important role in mitigating the accumulation of Hac1^up to prevent accidental UPR activation in the Saccharomyces yeast.

High-Throughput Affinity Measurements of Transcription Factor and DNA Mutations Reveal Affinity and Specificity Determinants

Transcription factors (TFs) bind regulatory DNA to control gene expression, and mutations to either TFs or DNA can alter binding affinities to rewire regulatory networks and drive phenotypic variation. While studies have profiled energetic effects of DNA mutations extensively, we lack similar information for TF variants. Here, we present STAMMP (simultaneous transcription factor affinity measurements via microfluidic protein arrays), a high-throughput microfluidic platform enabling quantitative characterization of hundreds of TF variants simultaneously. Measured affinities for ∼210 mutants of a model yeast TF (Pho4) interacting with 9 oligonucleotides (>1,800 K_ds) reveal that many combinations of mutations to poorly conserved TF residues and nucleotides flanking the core binding site alter but preserve physiological binding, providing a mechanism by which combinations of mutations in cis and trans could modulate TF binding to tune occupancies during evolution. Moreover, biochemical double-mutant cycles across the TF-DNA interface reveal molecular mechanisms driving recognition, linking sequence to function. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

Beyond Trees: Regulons and Regulatory Motif Characterization

Trees and their seeds regulate their germination, growth, and reproduction in response to environmental stimuli. These stimuli, through signal transduction, trigger transcription factors that alter the expression of various genes leading to the unfolding of the genetic program. A regulon is conceptually defined as a set of target genes regulated by a transcription factor by physically binding to regulatory motifs to accomplish a specific biological function, such as the CO-FT regulon for flowering timing and fall growth cessation in trees. Only with a clear characterization of regulatory motifs, can candidate target genes be experimentally validated, but motif characterization represents the weakest feature of regulon research, especially in tree genetics. I review here relevant experimental and bioinformatics approaches in characterizing transcription factors and their binding sites, outline problems in tree regulon research, and demonstrate how transcription factor databases can be effectively used to aid the characterization of tree regulons.

The ancillary N-terminal region of the yeast AP-1 transcription factor Yap8 contributes to its DNA binding specificity

Activator protein 1 (AP-1) is one of the largest families of basic leucine zipper (bZIP) transcription factors in eukaryotic cells. How AP-1 proteins achieve target DNA binding specificity remains elusive. In Saccharomyces cerevisiae, the AP-1-like protein (Yap) family comprises eight members (Yap1 to Yap8) that display distinct genomic target sites despite high sequence homology of their DNA binding bZIP domains. In contrast to the other members of the Yap family, which preferentially bind to short (7–8 bp) DNA motifs, Yap8 binds to an unusually long DNA motif (13 bp). It has been unclear what determines this unique specificity of Yap8. In this work, we use molecular and biochemical analyses combined with computer-based structural design and molecular dynamics simulations of Yap8–DNA interactions to better understand the structural basis of DNA binding specificity determinants. We identify specific residues in the N-terminal tail preceding the basic region, which define stable association of Yap8 with its target promoter. We propose that the N-terminal tail directly interacts with DNA and stabilizes Yap8 binding to the 13 bp motif. Thus, beside the core basic region, the adjacent N-terminal region contributes to alternative DNA binding selectivity within the AP-1 family.

Stress sensor Ire1 deploys a divergent transcriptional program in response to lipid bilayer stress

Membrane integrity at the endoplasmic reticulum (ER) is tightly regulated, and its disturbance is implicated in metabolic diseases. Using an engineered sensor that activates the unfolded protein response (UPR) exclusively when normal ER membrane lipid composition is compromised, we identified pathways beyond lipid metabolism that are necessary to maintain ER integrity in yeast and in C. elegans. To systematically validate yeast mutants that disrupt ER membrane homeostasis, we identified a lipid bilayer stress (LBS) sensor in the UPR transducer protein Ire1, located at the interface of the amphipathic and transmembrane helices. Furthermore, transcriptome and chromatin immunoprecipitation analyses pinpoint the UPR as a broad-spectrum compensatory response wherein LBS and proteotoxic stress deploy divergent transcriptional UPR programs. Together, these findings reveal the UPR program as the sum of two independent stress responses, an insight that could be exploited for future therapeutic intervention.

MODER2: first-order Markov modeling and discovery of monomeric and dimeric binding motifs

MOTIVATION

Position-specific probability matrices (PPMs, also called position-specific weight matrices) have been the dominating model for transcription factor (TF)-binding motifs in DNA. There is, however, increasing recent evidence of better performance of higher order models such as Markov models of order one, also called adjacent dinucleotide matrices (ADMs). ADMs can model dependencies between adjacent nucleotides, unlike PPMs. A modeling technique and software tool that would estimate such models simultaneously both for monomers and their dimers have been missing.

RESULTS

We present an ADM-based mixture model for monomeric and dimeric TF-binding motifs and an expectation maximization algorithm MODER2 for learning such models from training data and seeds. The model is a mixture that includes monomers and dimers, built from the monomers, with a description of the dimeric structure (spacing, orientation). The technique is modular, meaning that the co-operative effect of dimerization is made explicit by evaluating the difference between expected and observed models. The model is validated using HT-SELEX and generated datasets, and by comparing to some earlier PPM and ADM techniques. The ADM models explain data slightly better than PPM models for 314 tested TFs (or their DNA-binding domains) from four families (bHLH, bZIP, ETS and Homeodomain), the ADM mixture models by MODER2 being the best on average.

AVAILABILITY AND IMPLEMENTATION

Software implementation is available from https://github.com/jttoivon/moder2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

HacA Governs Virulence Traits and Adaptive Stress Responses in Trichophyton rubrum

The ability of fungi to sense environmental stressors and appropriately respond is linked to secretory system functions. The dermatophyte infection process depends on an orchestrated signaling regulation that triggers the transcription of genes responsible for adherence and penetration of the pathogen into host-tissue. A high secretion system is activated to support the host-pathogen interaction and assures maintenance of the dermatophyte infection. The gateway of secretion machinery is the endoplasmic reticulum (ER), which is the primary site for protein folding and transport. Current studies have shown that ER stress that affects adaptive responses is primarily regulated by UPR and supports fungal pathogenicity; this has been assessed for yeasts and Aspergillus fumigatus, in regard to how these fungi cope with host environmental stressors. Fungal UPR consists of a transmembrane kinase sensor (Ire1/IreA) and a downstream target Hac1/HacA. The active form of Hac is achieved via non-spliceosomal intron removal promoted by endonuclease activity of Ire1/IreA. Here, we assessed features of HacA and its involvement in virulence and susceptibility in Trichophyton rubrum. Our results showed that exposure to antifungals and ER-stressing agents initiated the activation of HacA from T. rubrum. Interestingly, the activation occurs when a 20 nt fragment is removed from part of the exon-2 and part of intron-2, which in turn promotes the arisen of the DNA binding site motif and a dimer interface domain. Further, we found changes in the cell wall and cellular membrane composition in the ΔhacA mutant as well as an increase in susceptibility toward azole and cell wall disturbing agents. Moreover, the ΔhacA mutant presented significant defects in important virulence traits like thermotolerance and growth on keratin substrates. For instance, the development of the ΔhacA mutant was impaired in co-culture with keratinocytes or human nail fragments. Changes in the pro-inflammatory cytokine release were verified for the ΔhacA mutant during the co-culture assay, which might be related to differences in pathogen-associated molecular patterns (PAMPs) in the cell wall. Together, these results suggested that HacA is an integral part of T. rubrum physiology and virulence, implying that it is an important molecular target for antidermatophytic therapy.

micrIO: an open-source autosampler and fraction collector for automated microfluidic input-output

Microfluidic devices are an enabling technology for many labs, facilitating a wide range of applications spanning high-throughput encapsulation, molecular separations, and long-term cell culture. In many cases, however, their utility is limited by a 'world-to-chip' barrier that makes it difficult to serially interface samples with these devices. As a result, many researchers are forced to rely on low-throughput, manual approaches for managing device input and output (IO) of samples, reagents, and effluent. Here, we present a hardware-software platform for automated microfluidic IO (micrIO). The platform, which is uniquely compatible with positive-pressure microfluidics, comprises an 'AutoSipper' for input and a 'Fraction Collector' for output. To facilitate widespread adoption, both are open-source builds constructed from components that are readily purchased online or fabricated from included design files. The software control library, written in Python, allows the platform to be integrated with existing experimental setups and to coordinate IO with other functions such as valve actuation and assay imaging. We demonstrate these capabilities by coupling both the AutoSipper and Fraction Collector to two microfluidic devices: a simple, valved inlet manifold and a microfluidic droplet generator that produces beads with distinct spectral codes. Analysis of the collected materials in each case establishes the ability of the platform to draw from and output to specific wells of multiwell plates with negligible cross-contamination between samples.

Optimized Sequence Library Design for Efficient In Vitro Interaction Mapping

Sequence libraries that cover all k-mers enable universal, unbiased measurements of binding to both oligonucleotides and peptides. While the number of k-mers grows exponentially in k, space on all experimental platforms is limited. Here, we shrink k-mer library sizes by using joker characters, which represent all characters in the alphabet simultaneously. We present the JokerCAKE (joker covering all k-mers) algorithm for generating a short sequence such that each k-mer appears at least p times with at most one joker character per k-mer. By running our algorithm on a range of parameters and alphabets, we show that JokerCAKE produces near-optimal sequences. Moreover, through comparison with data from hundreds of DNA-protein binding experiments and with new experimental results for both standard and JokerCAKE libraries, we establish that accurate binding scores can be inferred for high-affinity k-mers using JokerCAKE libraries. JokerCAKE libraries allow researchers to search a significantly larger sequence space using the same number of experimental measurements and at the same cost.

Signal peptide peptidase activity connects the unfolded protein response to plant defense suppression by Ustilago maydis

The corn smut fungus Ustilago maydis requires the unfolded protein response (UPR) to maintain homeostasis of the endoplasmic reticulum (ER) during the biotrophic interaction with its host plant Zea mays (maize). Crosstalk between the UPR and pathways controlling pathogenic development is mediated by protein-protein interactions between the UPR regulator Cib1 and the developmental regulator Clp1. Cib1/Clp1 complex formation results in mutual modification of the connected regulatory networks thereby aligning fungal proliferation in planta, efficient effector secretion with increased ER stress tolerance and long-term UPR activation in planta. Here we address UPR-dependent gene expression and its modulation by Clp1 using combinatorial RNAseq/ChIPseq analyses. We show that increased ER stress resistance is connected to Clp1-dependent alterations of Cib1 phosphorylation, protein stability and UPR gene expression. Importantly, we identify by deletion screening of UPR core genes the signal peptide peptidase Spp1 as a novel key factor that is required for establishing a compatible biotrophic interaction between U. maydis and its host plant maize. Spp1 is dispensable for ER stress resistance and vegetative growth but requires catalytic activity to interfere with the plant defense, revealing a novel virulence specific function for signal peptide peptidases in a biotrophic fungal/plant interaction.

Biotrophic pathogens establish compatible interactions with their host to cause disease.

A critical step in this process is the suppression of plant defense responses by secreted effector proteins. In the maize infecting fungus Ustilago maydis expression of effector encoding genes is coordinately upregulated at defined stages of pathogenic development in so-called effector waves. Efficient secretion of the multitude of effectors relies on the unfolded protein response (UPR) to maintain homeostasis of the endoplasmic reticulum. Activation of the UPR is connected to the control of fungal proliferation through direct protein-protein interactions between the UPR regulator Cib1 and the developmental regulator Clp1. Here, we show that this interaction leads to functional modification of Cib1 and modulation of UPR gene expression to adapt the UPR for long-term activity in the plant. Within a core set of UPR regulated genes we identify the signal peptide peptidase Spp1 as a key factor for fungal virulence. We show that Spp1 requires its conserved catalytic activity to suppress the plant defense and cause disease. The virulence specific function of Spp1 does not involve pathways previously known to be associated with Spp1-like proteins or plant defense suppression, suggesting a novel role for Spp1 substrates in biotrophic interactions.

Stress-sensing and regulatory mechanism of the endoplasmic-stress sensors Ire1 and PERK

Ire1 and its family protein PERK are endoplasmic reticulum (ER)-stress sensors that initiate cellular responses against ER accumulation of unfolded proteins. As reviewed in this article, many publications describe molecular mechanisms by which yeast Ire1 senses ER conditions and gets regulated. We also cover recent studies which reveal that mammalian Ire1 (IRE1α) and PERK are controlled in a similar but not exactly the same manner. ER-located molecular chaperone BiP captures these ER-stress sensors and suppresses their activity. Intriguingly, Ire1 is associated with BiP not as a chaperone substrate, but as a unique ligand. Unfolded proteins accumulated in the ER promote dissociation of the Ire1-BiP complex. Moreover, Ire1 is directly bound with unfolded proteins, leading to its cluster formation and potent activation. PERK also captures unfolded proteins and then forms self-oligomers. Meanwhile, membrane-lipid aberrancy is likely to activate these ER-stress sensors independently of ER accumulation of unfolded proteins. In addition, there exist a number of reports that touch on other factors that control activity of these ER-stress sensors. Such a multiplicity of regulatory mechanisms for these ER-stress sensors is likely to contribute to fine tuning of their activity.

The Unfolded Protein Response Pathway in the Yeast Kluyveromyces lactis. A Comparative View among Yeast Species

Eukaryotic cells have evolved signalling pathways that allow adaptation to harmful conditions that disrupt endoplasmic reticulum (ER) homeostasis. When the function of the ER is compromised in a condition known as ER stress, the cell triggers the unfolded protein response (UPR) in order to restore ER homeostasis. Accumulation of misfolded proteins due to stress conditions activates the UPR pathway. In mammalian cells, the UPR is composed of three branches, each containing an ER sensor (PERK, ATF6 and IRE1). However, in yeast species, the only sensor present is the inositol-requiring enzyme Ire1. To cope with unfolded protein accumulation, Ire1 triggers either a transcriptional response mediated by a transcriptional factor that belongs to the bZIP transcription factor family or an mRNA degradation process. In this review, we address the current knowledge of the UPR pathway in several yeast species: Saccharomyces cerevisiae, Schizosaccharomyces pombe, Candida glabrata, Cryptococcus neoformans, and Candida albicans. We also include unpublished data on the UPR pathway of the budding yeast Kluyveromyces lactis. We describe the basic components of the UPR pathway along with similarities and differences in the UPR mechanism that are present in these yeast species.

Comprehensive, high-resolution binding energy landscapes reveal context dependencies of transcription factor binding

Transcription factors (TFs) are primary regulators of gene expression in cells, where they bind specific genomic target sites to control transcription. Quantitative measurements of TF-DNA binding energies can improve the accuracy of predictions of TF occupancy and downstream gene expression in vivo and shed light on how transcriptional networks are rewired throughout evolution. Here, we present a sequencing-based TF binding assay and analysis pipeline (BET-seq, for Binding Energy Topography by sequencing) capable of providing quantitative estimates of binding energies for more than one million DNA sequences in parallel at high energetic resolution. Using this platform, we measured the binding energies associated with all possible combinations of 10 nucleotides flanking the known consensus DNA target interacting with two model yeast TFs, Pho4 and Cbf1. A large fraction of these flanking mutations change overall binding energies by an amount equal to or greater than consensus site mutations, suggesting that current definitions of TF binding sites may be too restrictive. By systematically comparing estimates of binding energies output by deep neural networks (NNs) and biophysical models trained on these data, we establish that dinucleotide (DN) specificities are sufficient to explain essentially all variance in observed binding behavior, with Cbf1 binding exhibiting significantly more nonadditivity than Pho4. NN-derived binding energies agree with orthogonal biochemical measurements and reveal that dynamically occupied sites in vivo are both energetically and mutationally distant from the highest affinity sites.

An amebic protein disulfide isomerase (PDI) complements the yeast PDI1 mutation but is unable to support cell viability under ER or thermal stress

In eukaryotic cells, protein disulfide isomerases (PDI) are oxidoreductases that catalyze the proper disulfide bond formation during protein folding. The pathobiology of the protozoan parasite Entamoeba histolytica, the causative agent of human amebiasis, depends on secretion of several virulence factors, such as pore‐forming peptides and cysteine proteinases. Although the native conformation of these factors is stabilized by disulfide bonds, there is little information regarding the molecular machinery involved in the oxidative folding of amebic proteins. Whereas testing gene function in their physiological background would be the most suitable approach, we have taken advantage of the cellular benefits offered by the yeast Saccharomyces cerevisiae (as a model of eukaryotic cell) to examine the functional role of an amebic PDI (EhPDI). As the yeast PDI homolog is essential for cell viability, a functional complementation assay was carried out to test the ability of EhPDI to circumvent the lethal phenotype of a yeast PDI1 mutant. Also, its proficiency under stressful conditions was explored by examining the survival outcome following endoplasmic reticulum (ER) stress induced by a reductant agent (DTT) or thermal stress promoted by a nonpermissive temperature (37 °C). Our results indicate that EhPDI is functionally active when physiological conditions are stable. Nonetheless, when conditions are stressful (e.g., by the accumulation of misfolded proteins in the ER compartment), its functionality is exceeded, suggesting an inability to prevent unfolding, suppress aggregation, or assist refolding of proteins. Despite the latter, our findings constitute the initial step toward determining the participation of EhPDI in cellular mechanisms related to protein homeostasis.

Bayesian Markov models consistently outperform PWMs at predicting motifs in nucleotide sequences

Position weight matrices (PWMs) are the standard model for DNA and RNA regulatory motifs. In PWMs nucleotide probabilities are independent of nucleotides at other positions. Models that account for dependencies need many parameters and are prone to overfitting. We have developed a Bayesian approach for motif discovery using Markov models in which conditional probabilities of order k - 1 act as priors for those of order k This Bayesian Markov model (BaMM) training automatically adapts model complexity to the amount of available data. We also derive an EM algorithm for de-novo discovery of enriched motifs. For transcription factor binding, BaMMs achieve significantly (P    =  1/16) higher cross-validated partial AUC than PWMs in 97% of 446 ChIP-seq ENCODE datasets and improve performance by 36% on average. BaMMs also learn complex multipartite motifs, improving predictions of transcription start sites, polyadenylation sites, bacterial pause sites, and RNA binding sites by 26-101%. BaMMs never performed worse than PWMs. These robust improvements argue in favour of generally replacing PWMs by BaMMs.

Quantification of Cooperativity in Heterodimer-DNA Binding Improves the Accuracy of Binding Specificity Models

Many transcription factors (TFs) have the ability to cooperate on DNA elements as heterodimers. Despite the significance of TF heterodimerization for gene regulation, a quantitative understanding of cooperativity between various TF dimer partners and its impact on heterodimer DNA binding specificity models is still lacking. Here, we used a novel integrative approach, combining microfluidics-steered measurements of dimer-DNA assembly with mechanistic modeling of the implicated protein-protein-DNA interactions to quantitatively interrogate the cooperative DNA binding behavior of the adipogenic peroxisome proliferator-activated receptor γ (PPARγ):retinoid X receptor α (RXRα) heterodimer. Using the high throughput MITOMI (mechanically induced trapping of molecular interactions) platform, we derived equilibrium DNA binding data for PPARγ, RXRα, as well as the PPARγ:RXRα heterodimer to more than 300 target DNA sites and variants thereof. We then quantified cooperativity underlying heterodimer-DNA binding and derived an integrative heterodimer DNA binding constant. Using this cooperativity-inclusive constant, we were able to build a heterodimer-DNA binding specificity model that has superior predictive power than the one based on a regular one-site equilibrium. Our data further revealed that individual nucleotide substitutions within the target site affect the extent of cooperativity in PPARγ:RXRα-DNA binding. Our study therefore emphasizes the importance of assessing cooperativity when generating DNA binding specificity models for heterodimers.

Unfolded Protein Response (UPR) Regulator Cib1 Controls Expression of Genes Encoding Secreted Virulence Factors in Ustilago maydis

The unfolded protein response (UPR), a conserved eukaryotic signaling pathway to ensure protein homeostasis in the endoplasmic reticulum (ER), coordinates biotrophic development in the corn smut fungus Ustilago maydis. Exact timing of UPR activation is required for virulence and presumably connected to the elevated expression of secreted effector proteins during infection of the host plant Zea mays. In the baker’s yeast Saccharomyces cerevisiae, expression of UPR target genes is induced upon binding of the central regulator Hac1 to unfolded protein response elements (UPREs) in their promoters. While a role of the UPR in effector secretion has been described previously, we investigated a potential UPR-dependent regulation of genes encoding secreted effector proteins. In silico prediction of UPREs in promoter regions identified the previously characterized effector genes pit2 and tin1-1, as bona fide UPR target genes. Furthermore, direct binding of the Hac1-homolog Cib1 to the UPRE containing promoter fragments of both genes was confirmed by quantitative chromatin immunoprecipitation (qChIP) analysis. Targeted deletion of the UPRE abolished Cib1-dependent expression of pit2 and significantly affected virulence. Furthermore, ER stress strongly increased Pit2 expression and secretion. This study expands the role of the UPR as a signal hub in fungal virulence and illustrates, how biotrophic fungi can coordinate cellular physiology, development and regulation of secreted virulence factors.

Microfluidics: reframing biological enquiry

The underlying physical properties of microfluidic tools have led to new biological insights through the development of microsystems that can manipulate, mimic and measure biology at a resolution that has not been possible with macroscale tools. Microsystems readily handle sub-microlitre volumes, precisely route predictable laminar fluid flows and match both perturbations and measurements to the length scales and timescales of biological systems. The advent of fabrication techniques that do not require highly specialized engineering facilities is fuelling the broad dissemination of microfluidic systems and their adaptation to specific biological questions. We describe how our understanding of molecular and cell biology is being and will continue to be advanced by precision microfluidic approaches and posit that microfluidic tools - in conjunction with advanced imaging, bioinformatics and molecular biology approaches - will transform biology into a precision science.

The bZIP Transcription Factor HAC-1 Is Involved in the Unfolded Protein Response and Is Necessary for Growth on Cellulose in Neurospora crassa

High protein secretion capacity in filamentous fungi requires an extremely efficient system for protein synthesis, folding and transport. When the folding capacity of the endoplasmic reticulum (ER) is exceeded, a pathway known as the unfolded protein response (UPR) is triggered, allowing cells to mitigate and cope with this stress. In yeast, this pathway relies on the transcription factor Hac1, which mediates the up-regulation of several genes required under these stressful conditions. In this work, we identified and characterized the ortholog of the yeast HAC1 gene in the filamentous fungus Neurospora crassa. We show that its mRNA undergoes an ER stress-dependent splicing reaction, which in N. crassa removes a 23 nt intron and leads to a change in the open reading frame. By disrupting the N. crassa hac-1 gene, we determined it to be crucial for activating UPR and for proper growth in the presence of ER stress-inducing chemical agents. Neurospora is naturally found growing on dead plant material, composed primarily by lignocellulose, and is a model organism for the study of plant cell wall deconstruction. Notably, we found that growth on cellulose, a substrate that requires secretion of numerous enzymes, imposes major demands on ER function and is dramatically impaired in the absence of hac-1, thus broadening the range of physiological functions of the UPR in filamentous fungi. Growth on hemicellulose however, another carbon source that necessitates the secretion of various enzymes for its deconstruction, is not impaired in the mutant nor is the amount of proteins secreted on this substrate, suggesting that secretion, as a whole, is unaltered in the absence of hac-1. The characterization of this signaling pathway in N. crassa will help in the study of plant cell wall deconstruction by fungi and its manipulation may result in important industrial biotechnological applications.

Mechanically Induced Trapping of Molecular Interactions and Its Applications

Measuring binding affinities and association/dissociation rates of molecular interactions is important for a quantitative understanding of cellular mechanisms. Many low-throughput methods have been developed throughout the years to obtain these parameters. Acquiring data with higher accuracy and throughput is, however, necessary to characterize complex biological networks. Here, we provide an overview of a high-throughput microfluidic method based on mechanically induced trapping of molecular interactions (MITOMI). MITOMI can be used to obtain affinity constants and kinetic rates of hundreds of protein-ligand interactions in parallel. It has been used in dozens of studies to measure binding affinities of transcription factors, map protein interaction networks, identify pharmacological inhibitors, and perform high-throughput, low-cost molecular diagnostics. This article covers the technological aspects of MITOMI and its applications.

Unfolded protein response in filamentous fungi-implications in biotechnology

The unfolded protein response (UPR) represents a mechanism to preserve endoplasmic reticulum (ER) homeostasis that is conserved in eukaryotes. ER stress caused by the accumulation of potentially toxic un- or misfolded proteins in the ER triggers UPR activation and the induction of genes important for protein folding in the ER, ER expansion, and transport from and to the ER. Along with this adaptation, the overall capacity for protein secretion is markedly increased by the UPR. In filamentous fungi, various approaches to employ the UPR for improved production of homologous and heterologous proteins have been investigated. As the effects on protein production were strongly dependent on the expressed protein, generally applicable strategies have to be developed. A combination of transcriptomic approaches monitoring secretion stress and basic research on the UPR mechanism provided novel and important insight into the complex regulatory cross-connections between UPR signalling, cellular physiology, and developmental processes. It will be discussed how this increasing knowledge on the UPR might stimulate the development of novel strategies for using the UPR as a tool in biotechnology.

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.

How duplicated transcription regulators can diversify to govern the expression of nonoverlapping sets of genes

The duplication of transcription regulators can elicit major regulatory network rearrangements over evolutionary timescales. However, few examples of duplications resulting in gene network expansions are understood in molecular detail. Here we show that four Candida albicans transcription regulators that arose by successive duplications have differentiated from one another by acquiring different intrinsic DNA-binding specificities, different preferences for half-site spacing, and different associations with cofactors. The combination of these three mechanisms resulted in each of the four regulators controlling a distinct set of target genes, which likely contributed to the adaption of this fungus to its human host. Our results illustrate how successive duplications and diversification of an ancestral transcription regulator can underlie major changes in an organism's regulatory circuitry.

Conserved and plant-unique strategies for overcoming endoplasmic reticulum stress

Stress caused by environmental conditions or physiological growth can lead to an accumulation of unfolded proteins in the endoplasmic reticulum (ER) causing ER stress, which in turn triggers a cytoprotective mechanism termed the unfolded protein response (UPR). Under mild-short stress conditions the UPR can restore ER functioning and cell growth, such as reducing the load of unfolded proteins through the upregulation of genes involved in protein folding and in degrading mis-folded proteins, and through autophagy activation, but it can also lead to cell death under prolonged and severe stress conditions. A diversified suite of sensors has been evolved in the eukaryotic lineages to orchestrate the UPR most likely to suit the cell's necessity to respond to the different kinds of stress in a conserved as well as species-specific manner. In plants three UPR sensors cooperate with non-identical signaling pathways: the protein kinase inositol-requiring enzyme (IRE1), the ER-membrane-associated transcription factor bZIP28, and the GTP-binding protein β1 (AGB1). In this mini-review, we show how plants differ from the better characterized metazoans and fungi, providing an overview of the signaling pathways of the UPR, and highlighting the overlapping and the peculiar roles of the different UPR branches in light of evolutionary divergences in eukaryotic kingdoms.

Protein-DNA binding: complexities and multi-protein codes

Binding of proteins to particular DNA sites across the genome is a primary determinant of specificity in genome maintenance and gene regulation. DNA-binding specificity is encoded at multiple levels, from the detailed biophysical interactions between proteins and DNA, to the assembly of multi-protein complexes. At each level, variation in the mechanisms used to achieve specificity has led to difficulties in constructing and applying simple models of DNA binding. We review the complexities in protein–DNA binding found at multiple levels and discuss how they confound the idea of simple recognition codes. We discuss the impact of new high-throughput technologies for the characterization of protein–DNA binding, and how these technologies are uncovering new complexities in protein–DNA recognition. Finally, we review the concept of multi-protein recognition codes in which new DNA-binding specificities are achieved by the assembly of multi-protein complexes.

Msn2 coordinates a stoichiometric gene expression program

BACKGROUND

Many cellular processes operate in an "analog" regime in which the magnitude of the response is precisely tailored to the intensity of the stimulus. In order to maintain the coherence of such responses, the cell must provide for proportional expression of multiple target genes across a wide dynamic range of induction states. Our understanding of the strategies used to achieve graded gene regulation is limited.

RESULTS

In this work, we document a relationship between stress-responsive gene expression and the transcription factor Msn2 that is graded over a large range of Msn2 concentrations. We use computational modeling and in vivo and in vitro analyses to dissect the roots of this relationship. Our studies reveal a simple and general strategy based on noncooperative low-affinity interactions between Msn2 and its cognate binding sites as well as competition over a large number of Msn2 binding sites in the genome relative to the number of Msn2 molecules.

CONCLUSIONS

In addition to enabling precise tuning of gene expression to the state of the environment, this strategy ensures colinear activation of target genes, allowing for stoichiometric expression of large groups of genes without extensive promoter tuning. Furthermore, such a strategy enables precise modulation of the activity of any given promoter by addition of binding sites without altering the qualitative relationship between different genes in a regulon. This feature renders a given regulon highly "evolvable."

Crosstalk between the unfolded protein response and pathways that regulate pathogenic development in Ustilago maydis

The unfolded protein response (UPR) is a conserved eukaryotic signaling pathway regulating endoplasmic reticulum (ER) homeostasis during ER stress, which results, for example, from an increased demand for protein secretion. Here, we characterize the homologs of the central UPR regulatory proteins Hac1 (for Homologous to ATF/CREB1) and Inositol Requiring Enzyme1 in the plant pathogenic fungus Ustilago maydis and demonstrate that the UPR is tightly interlinked with the b mating-type-dependent signaling pathway that regulates pathogenic development. Exact timing of UPR is required for virulence, since premature activation interferes with the b-dependent switch from budding to filamentous growth. In addition, we found crosstalk between UPR and the b target Clampless1 (Clp1), which is essential for cell cycle release and proliferation in planta. The unusual C-terminal extension of the U. maydis Hac1 homolog, Cib1 (for Clp1 interacting bZIP1), mediates direct interaction with Clp1. The interaction between Clp1 and Cib1 promotes stabilization of Clp1, resulting in enhanced ER stress tolerance that prevents deleterious UPR hyperactivation. Thus, the interaction between Cib1 and Clp1 constitutes a checkpoint to time developmental progression and increased secretion of effector proteins at the onset of biotrophic development. Crosstalk between UPR and the b mating-type regulated developmental program adapts ER homeostasis to the changing demands during biotrophy.

The plant-specific transcription factor gene NAC103 is induced by bZIP60 through a new cis-regulatory element to modulate the unfolded protein response in Arabidopsis

The unfolded protein response (UPR) plays important roles in plant development and plant-pathogen interactions, as well as in plant adaptation to adverse environmental stresses. Previously bZIP28 and bZIP60 have been identified as important UPR regulators for mitigating the endoplasmic reticulum (ER) stress in Arabidopsis thaliana. Here we report the biological function of NAC103 in a novel transcriptional regulatory cascade, connecting bZIP60 to the UPR downstream genes in Arabidopsis. Expression of NAC103 was induced by ER stress, and was completely abolished in the bZIP60 null mutant. A new ER stress-responsive cis-element UPRE-III (TCATCG) on the NAC103 promoter was identified, and trans-activation of UPRE-III by bZIP60 was confirmed in both yeast cells and Arabidopsis protoplasts. The direct binding of bZIP60 to UPRE-III-containing DNA was also demonstrated in an electrophoretic mobility shift assay. NAC103 formed homodimers in yeast two-hybrid and bimolecular fluorescence complementation assays. It had transcriptional activation activity and was localized in the nucleus. Over-expression of NAC103 had pleiotropic effects on plant growth, and induced expression of several UPR downstream genes in Arabidopsis under normal growth conditions. The activation of UPR gene promoters by NAC103 was also confirmed in effector/reporter protoplast assays. Thus, our study demonstrates a transcriptional regulatory cascade in which NAC103 relays ER stress signals from bZIP60 to UPR downstream genes through a newly identified ER stress cis-element (UPRE-III) and transcriptional activation activity of its encoded protein NAC103.

Structure of the transcriptional network controlling white-opaque switching in Candida albicans

The human fungal pathogen Candida albicans can switch between two phenotypic cell types, termed 'white' and 'opaque'. Both cell types are heritable for many generations, and the switch between the two types occurs epigenetically, that is, without a change in the primary DNA sequence of the genome. Previous work identified six key transcriptional regulators important for white-opaque switching: Wor1, Wor2, Wor3, Czf1, Efg1, and Ahr1. In this work, we describe the structure of the transcriptional network that specifies the white and opaque cell types and governs the ability to switch between them. In particular, we use a combination of genome-wide chromatin immunoprecipitation, gene expression profiling, and microfluidics-based DNA binding experiments to determine the direct and indirect regulatory interactions that form the switch network. The six regulators are arranged together in a complex, interlocking network with many seemingly redundant and overlapping connections. We propose that the structure (or topology) of this network is responsible for the epigenetic maintenance of the white and opaque states, the switching between them, and the specialized properties of each state.

LASAGNA-Search: an integrated web tool for transcription factor binding site search and visualization

The release of ChIP-seq data from the ENCyclopedia Of DNA Elements (ENCODE) and Model Organism ENCyclopedia Of DNA Elements (modENCODE) projects has significantly increased the amount of transcription factor (TF) binding affinity information available to researchers. However, scientists still routinely use TF binding site (TFBS) search tools to scan unannotated sequences for TFBSs, particularly when searching for lesser-known TFs or TFs in organisms for which ChIP-seq data are unavailable. The sequence analysis often involves multiple steps such as TF model collection, promoter sequence retrieval, and visualization; thus, several different tools are required. We have developed a novel integrated web tool named LASAGNA-Search that allows users to perform TFBS searches without leaving the web site. LASAGNA-Search uses the LASAGNA (Length-Aware Site Alignment Guided by Nucleotide Association) algorithm for TFBS alignment. Important features of LASAGNA-Search include (i) acceptance of unaligned variable-length TFBSs, (ii) a collection of 1726 TF models, (iii) automatic promoter sequence retrieval, (iv) visualization in the UCSC Genome Browser, and (v) gene regulatory network inference and visualization based on binding specificities. LASAGNA-Search is freely available at http://biogrid.engr.uconn.edu/lasagna_search/.

Mutual cross talk between the regulators Hac1 of the unfolded protein response and Gcn4 of the general amino acid control of Saccharomyces cerevisiae

Hac1 is the activator of the cellular response to the accumulation of unfolded proteins in the endoplasmic reticulum. Hac1 function requires the activity of Gcn4, which mainly acts as a regulator of the general amino acid control network providing Saccharomyces cerevisiae cells with amino acids. Here, we demonstrate novel functions of Hac1 and describe a mutual connection between Hac1 and Gcn4. Hac1 is required for induction of Gcn4-responsive promoter elements in haploid as well as diploid cells and therefore participates in the cellular amino acid supply. Furthermore, Hac1 and Gcn4 mutually influence their mRNA expression levels. Hac1 is also involved in FLO11 expression and adhesion upon amino acid starvation. Hac1 and Gcn4 act through the same promoter regions of the FLO11 flocculin. The results indicate an indirect effect of both transcription factors on FLO11 expression. Our data suggest a complex mutual cross talk between the Hac1- and Gcn4-controlled networks.

Microfluidic affinity and ChIP-seq analyses converge on a conserved FOXP2-binding motif in chimp and human, which enables the detection of evolutionarily novel targets

The transcription factor forkhead box P2 (FOXP2) is believed to be important in the evolution of human speech. A mutation in its DNA-binding domain causes severe speech impairment. Humans have acquired two coding changes relative to the conserved mammalian sequence. Despite intense interest in FOXP2, it has remained an open question whether the human protein’s DNA-binding specificity and chromatin localization are conserved. Previous in vitro and ChIP-chip studies have provided conflicting consensus sequences for the FOXP2-binding site. Using MITOMI 2.0 microfluidic affinity assays, we describe the binding site of FOXP2 and its affinity profile in base-specific detail for all substitutions of the strongest binding site. We find that human and chimp FOXP2 have similar binding sites that are distinct from previously suggested consensus binding sites. Additionally, through analysis of FOXP2 ChIP-seq data from cultured neurons, we find strong overrepresentation of a motif that matches our in vitro results and identifies a set of genes with FOXP2 binding sites. The FOXP2-binding sites tend to be conserved, yet we identified 38 instances of evolutionarily novel sites in humans. Combined, these data present a comprehensive portrait of FOXP2’s-binding properties and imply that although its sequence specificity has been conserved, some of its genomic binding sites are newly evolved.

Identification and characterization of a previously undescribed family of sequence-specific DNA-binding domains

Sequence-specific DNA-binding proteins are among the most important classes of gene regulatory proteins, controlling changes in transcription that underlie many aspects of biology. In this work, we identify a transcriptional regulator from the human fungal pathogen Candida albicans that binds DNA specifically but has no detectable homology with any previously described DNA- or RNA-binding protein. This protein, named White-Opaque Regulator 3 (Wor3), regulates white-opaque switching, the ability of C. albicans to switch between two heritable cell types. We demonstrate that ectopic overexpression of WOR3 results in mass conversion of white cells to opaque cells and that deletion of WOR3 affects the stability of opaque cells at physiological temperatures. Genome-wide chromatin immunoprecipitation of Wor3 and gene expression profiling of a wor3 deletion mutant strain indicate that Wor3 is highly integrated into the previously described circuit regulating white-opaque switching and that it controls a subset of the opaque transcriptional program. We show by biochemical, genetic, and microfluidic experiments that Wor3 binds directly to DNA in a sequence-specific manner, and we identify the set of cis-regulatory sequences recognized by Wor3. Bioinformatic analyses indicate that the Wor3 family arose more recently in evolutionary time than most previously described DNA-binding domains; it is restricted to a small number of fungi that include the major fungal pathogens of humans. These observations show that new families of sequence-specific DNA-binding proteins may be restricted to small clades and suggest that current annotations--which rely on deep conservation--underestimate the fraction of genes coding for transcriptional regulators.

Quantitative microfluidic biomolecular analysis for systems biology and medicine

In the postgenome era, biology and medicine are rapidly evolving towards quantitative and systems studies of complex biological systems. Emerging breakthroughs in microfluidic technologies and innovative applications are transforming systems biology by offering new capabilities to address the challenges in many areas, such as single-cell genomics, gene regulation networks, and pathology. In this review, we focus on recent progress in microfluidic technology from the perspective of its applications to promoting quantitative and systems biomolecular analysis in biology and medicine.