Inferring direct regulatory targets from expression and genome location analyses: a comparison of transcription factor deletion and overexpression

Branched-chain amino acid biosynthesis in fungi

Branched-chain amino acids (BCAAs)-isoleucine, leucine, and valine-are synthesized by fungi. These amino acids are important components of proteins and secondary metabolites. The biochemical pathway for BCAA biosynthesis is well-characterized in the yeast Saccharomyces cerevisiae. The biosynthesis of these three amino acids is interconnected. Different precursors are metabolized in multiple steps through shared enzymes to produce isoleucine and valine, and the valine biosynthesis pathway branches before the penultimate step to a series of leucine biosynthesis-specific steps to produce leucine. Recent efforts have made advances toward characterization of the BCAA biosynthesis pathway in several fungi, revealing diversity in gene duplication and functional divergence in the genes for these enzymatic steps in different fungi. The BCAA biosynthesis pathway is regulated by the transcription factor LEU3 in S. cerevisiae, and LeuB in Aspergillus nidulans and Aspergillus fumigatus, and the activity of these transcription factors is modulated by the leucine biosynthesis pathway intermediate α-isopropylmalate. Herein, we discuss recent advances in our understanding of the BCAA pathway and its regulation, focusing on filamentous ascomycete fungi and comparison with the well-established process in yeast.

Novel function of transcription factor Uga3 as an activator of branched-chain amino acid permease BAP2 gene expression

Gene regulation in yeast occurs at the transcription level, i.e. the basal level of expression is very low and increased transcription requires gene-specific transcription factors allowing the recruitment of basal transcriptional machinery. Saccharomyces cerevisiae BAP2 gene encodes the permease responsible for most uptake of leucine, valine and isoleucine, amino acids that this yeast can use as nitrogen sources. Moreover, BAP2 expression is known to be induced by the presence of amino acids such as leucine. In this context, the results presented in this paper show that BAP2 is an inducible gene in the presence of nitrogen-non-preferred source proline but exhibits high constitutive non-inducible expression in nitrogen-preferred source ammonium. BAP2 expression is regulated by the SPS sensor system and transcription factors Leu3, Gcn4 and Dal81. This can be achieved or not through a direct binding to the promoter depending on the quality of the nitrogen source. We further demonstrate here that an interaction occurs in vivo between Uga3 ‒ the transcriptional activator responsible for γ-aminobutyric acid (GABA)-dependent induction of the GABA genes ‒ and the regulatory region of the BAP2 gene, which leads to an increase in BAP2 transcription.

Yeast response and tolerance to benzoic acid involves the Gcn4- and Stp1-regulated multidrug/multixenobiotic resistance transporter Tpo1

The action of benzoic acid in the food and beverage industries is compromised by the ability of spoilage yeasts to cope with this food preservative. Benzoic acid occurs naturally in many plants and is an intermediate compound in the biosynthesis of many secondary metabolites. The understanding of the mechanisms underlying the response and resistance to benzoic acid stress in the eukaryotic model yeast is thus crucial to design more suitable strategies to deal with this toxic lipophilic weak acid. In this study, the Saccharomyces cerevisiae multidrug transporter Tpo1 was demonstrated to confer resistance to benzoic acid. TPO1 transcript levels were shown to be up-regulated in yeast cells suddenly exposed to this stress agent. This up-regulation is under the control of the Gcn4 and Stp1 transcription factors, involved in the response to amino acid availability, but not under the regulation of the multidrug resistance transcription factors Pdr1 and Pdr3 that have binding sites in TPO1 promoter region. Benzoic acid stress was further shown to affect the intracellular pool of amino acids and polyamines. The observed decrease in the concentration of these nitrogenous compounds, registered upon benzoic acid stress exposure, was not found to be dependent on Tpo1, although the limitation of yeast cells on nitrogenous compounds was found to activate Tpo1 expression. Altogether, the results described in this study suggest that Tpo1 is one of the key players standing in the crossroad between benzoic acid stress response and tolerance and the control of the intracellular concentration of nitrogenous compounds. Also, results can be useful to guide the design of more efficient preservation strategies and the biotechnological synthesis of benzoic acid or benzoic acid-derived compounds.

A synthetic gene circuit for measuring autoregulatory feedback control

Autoregulatory feedback loops occur in the regulation of molecules ranging from ATP to MAP kinases to zinc. Negative feedback loops can increase a system's robustness, while positive feedback loops can mediate transitions between cell states. Recent genome-wide experimental and computational studies predict hundreds of novel feedback loops. However, not all physical interactions are regulatory, and many experimental methods cannot detect self-interactions. Our understanding of regulatory feedback loops is therefore hampered by the lack of high-throughput methods to experimentally quantify the presence, strength and temporal dynamics of autoregulatory feedback loops. Here we present a mathematical and experimental framework for high-throughput quantification of feedback regulation and apply it to RNA binding proteins (RBPs) in yeast. Our method is able to determine the existence of both direct and indirect positive and negative feedback loops, and to quantify the strength of these loops. We experimentally validate our model using two RBPs which lack native feedback loops and by the introduction of synthetic feedback loops. We find that RBP Puf3 does not natively participate in any direct or indirect feedback regulation, but that replacing the native 3'UTR with that of COX17 generates an auto-regulatory negative feedback loop which reduces gene expression noise. Likewise, RBP Pub1 does not natively participate in any feedback loops, but a synthetic positive feedback loop involving Pub1 results in increased expression noise. Our results demonstrate a synthetic experimental system for quantifying the existence and strength of feedback loops using a combination of high-throughput experiments and mathematical modeling. This system will be of great use in measuring auto-regulatory feedback by RNA binding proteins, a regulatory motif that is difficult to quantify using existing high-throughput methods.

Parameter Estimation for Gene Regulatory Networks from Microarray Data: Cold Shock Response in Saccharomyces cerevisiae

We investigated the dynamics of a gene regulatory network controlling the cold shock response in budding yeast, Saccharomyces cerevisiae. The medium-scale network, derived from published genome-wide location data, consists of 21 transcription factors that regulate one another through 31 directed edges. The expression levels of the individual transcription factors were modeled using mass balance ordinary differential equations with a sigmoidal production function. Each equation includes a production rate, a degradation rate, weights that denote the magnitude and type of influence of the connected transcription factors (activation or repression), and a threshold of expression. The inverse problem of determining model parameters from observed data is our primary interest. We fit the differential equation model to published microarray data using a penalized nonlinear least squares approach. Model predictions fit the experimental data well, within the 95 % confidence interval. Tests of the model using randomized initial guesses and model-generated data also lend confidence to the fit. The results have revealed activation and repression relationships between the transcription factors. Sensitivity analysis indicates that the model is most sensitive to changes in the production rate parameters, weights, and thresholds of Yap1, Rox1, and Yap6, which form a densely connected core in the network. The modeling results newly suggest that Rap1, Fhl1, Msn4, Rph1, and Hsf1 play an important role in regulating the early response to cold shock in yeast. Our results demonstrate that estimation for a large number of parameters can be successfully performed for nonlinear dynamic gene regulatory networks using sparse, noisy microarray data.

Analysis of Candida albicans mutants defective in the Cdk8 module of mediator reveal links between metabolism and biofilm formation

Candida albicans biofilm formation is a key virulence trait that involves hyphal growth and adhesin expression. Pyocyanin (PYO), a phenazine secreted by Pseudomonas aeruginosa, inhibits both C. albicans biofilm formation and development of wrinkled colonies. Using a genetic screen, we identified two mutants, ssn3Δ/Δ and ssn8Δ/Δ, which continued to wrinkle in the presence of PYO. Ssn8 is a cyclin-like protein and Ssn3 is similar to cyclin-dependent kinases; both proteins are part of the heterotetrameric Cdk8 module that forms a complex with the transcriptional co-regulator, Mediator. Ssn3 kinase activity was also required for PYO sensitivity as a kinase dead mutant maintained a wrinkled colony morphology in the presence of PYO. Furthermore, similar phenotypes were observed in mutants lacking the other two components of the Cdk8 module-Srb8 and Srb9. Through metabolomics analyses and biochemical assays, we showed that a compromised Cdk8 module led to increases in glucose consumption, glycolysis-related transcripts, oxidative metabolism and ATP levels even in the presence of PYO. In the mutant, inhibition of respiration to levels comparable to the PYO-treated wild type inhibited wrinkled colony development. Several lines of evidence suggest that PYO does not act through Cdk8. Lastly, the ssn3 mutant was a hyperbiofilm former, and maintained higher biofilm formation in the presence of PYO than the wild type. Together these data provide novel insights into the role of the Cdk8 module of Mediator in regulation of C. albicans physiology and the links between respiratory activity and both wrinkled colony and biofilm development.

Amino acid homeostasis and chronological longevity in Saccharomyces cerevisiae

Understanding how non-dividing cells remain viable over long periods of time, which may be decades in humans, is of central importance in understanding mechanisms of aging and longevity. The long-term viability of non-dividing cells, known as chronological longevity, relies on cellular processes that degrade old components and replace them with new ones. Key among these processes is amino acid homeostasis. Amino acid homeostasis requires three principal functions: amino acid uptake, de novo synthesis, and recycling. Autophagy plays a key role in recycling amino acids and other metabolic building blocks, while at the same time removing damaged cellular components such as mitochondria and other organelles. Regulation of amino acid homeostasis and autophagy is accomplished by a complex web of pathways that interact because of the functional overlap at the level of recycling. It is becoming increasingly clear that amino acid homeostasis and autophagy play important roles in chronological longevity in yeast and higher organisms. Our goal in this chapter is to focus on mechanisms and pathways that link amino acid homeostasis, autophagy, and chronological longevity in yeast, and explore their relevance to aging and longevity in higher eukaryotes.

Uga3 and Uga35/Dal81 transcription factors regulate UGA4 transcription in response to gamma-aminobutyric acid and leucine

The Saccharomyces cerevisiae UGA4 gene encodes a permease capable of importing gamma-aminobutyric acid (GABA) and delta-aminolevulinic acid (ALA) into the cell. GABA-dependent induction of this permease requires at least two positive-acting proteins, the specific factor Uga3 and the pleiotropic factor Uga35/Dal81. UGA4 is subjected to a very complex regulation, and its induction is affected by the presence of extracellular amino acids; this effect is mediated by the plasma membrane amino acid sensor SPS. Our results show that leucine affects UGA4 induction and that the SPS sensor and the downstream effectors Stp1 and Stp2 participate in this regulation. Moreover, we found that the Uga3 and Uga35/Dal81 transcription factors bind to the UGA4 promoter in a GABA-dependent manner and that this binding is impaired by the presence of leucine. We also found that the Leu3 transcription factor negatively regulates UGA4 transcription, although this seems to be through an indirect mechanism.

Characterizing regulatory path motifs in integrated networks using perturbational data

Pathicular – a Cytoscape plugin for analysing cellular responses to transcription factor perturbations is presented

We introduce Pathicular http://bioinformatics.psb.ugent.be/software/details/Pathicular, a Cytoscape plugin for studying the cellular response to perturbations of transcription factors by integrating perturbational expression data with transcriptional, protein-protein and phosphorylation networks. Pathicular searches for 'regulatory path motifs', short paths in the integrated physical networks which occur significantly more often than expected between transcription factors and their targets in the perturbational data. A case study in Saccharomyces cerevisiae identifies eight regulatory path motifs and demonstrates their biological significance.

Autophagy and amino acid homeostasis are required for chronological longevity in Saccharomyces cerevisiae

Following cessation of growth, yeast cells remain viable in a nondividing state for a period of time known as the chronological lifespan (CLS). Autophagy is a degradative process responsible for amino acid recycling in response to nitrogen starvation and amino acid limitation. We have investigated the role of autophagy during chronological aging of yeast grown in glucose minimal media containing different supplemental essential and nonessential amino acids. Deletion of ATG1 or ATG7, both of which are required for autophagy, reduced CLS, whereas deletion of ATG11, which is required for selective targeting of cellular components to the vacuole for degradation, did not reduce CLS. The nonessential amino acids isoleucine and valine, and the essential amino acid leucine, extended CLS in autophagy-deficient as well as autophagy-competent yeast. This extension was suppressed by constitutive expression of GCN4, which encodes a transcriptional regulator of general amino acid control (GAAC). Consistent with this, GCN4 expression was reduced by isoleucine and valine. Furthermore, elimination of the leucine requirement extended CLS and prevented the effects of constitutive expression of GCN4. Interestingly, deletion of LEU3, a GAAC target gene encoding a transcriptional regulator of branched side chain amino acid synthesis, dramatically increased CLS in the absence of amino acid supplements. In general, this indicates that activation of GAAC reduces CLS whereas suppression of GAAC extends CLS in minimal medium. These findings demonstrate important roles for autophagy and amino acid homeostasis in determining CLS in yeast.

Identification of the Candida albicans Cap1p regulon

Cap1p, a transcription factor of the basic region leucine zipper family, regulates the oxidative stress response (OSR) in Candida albicans. Alteration of its C-terminal cysteine-rich domain (CRD) results in Cap1p nuclear retention and transcriptional activation. To better understand the function of Cap1p in C. albicans, we used genome-wide location profiling (chromatin immunoprecipitation-on-chip) to identify its transcriptional targets in vivo. A triple-hemagglutinin (HA(3)) epitope was introduced at the C terminus of wild-type Cap1p (Cap1p-HA(3)) or hyperactive Cap1p with an altered CRD (Cap1p-CSE-HA(3)). Location profiling using whole-genome oligonucleotide tiling microarrays identified 89 targets bound by Cap1p-HA(3) or Cap1p-CSE-HA(3) (the binding ratio was at least twofold; P < or = 0.01). Strikingly, Cap1p binding was detected not only at the promoter region of its target genes but also at their 3' ends and within their open reading frames, suggesting that Cap1p may associate with the transcriptional or chromatin remodeling machinery to exert its activity. Overrepresented functional groups of the Cap1p targets (P < or = 0.02) included 11 genes involved in the OSR (CAP1, GLR1, TRX1, SOD1, CAT1, and others), 13 genes involved in response to drugs (PDR16, MDR1, FLU1, YCF1, FCR1, and others), 4 genes involved in phospholipid transport (PDR16, GIT1, RTA2, and orf19.932), and 3 genes involved in the regulation of nitrogen utilization (GST3, orf19.2693, and orf19.3121), suggesting that Cap1p has other cellular functions in addition to the OSR. Bioinformatic analyses of the bound sequences suggest that Cap1p recognizes the DNA motif 5'-MTKASTMA. Finally, transcriptome analyses showed that increased expression generally accompanies Cap1p binding at its targets, indicating that Cap1p functions as a transcriptional activator.

Inferring direct regulatory targets of a transcription factor in the DREAM2 challenge

In the DREAM2 community-wide experiment on regulatory network inference, one of the challenges was to identify which genes, in a list of 200, are direct regulatory targets of the transcription factor BCL6. The organizers of the challenge defined targets based on gene expression and chromatin immunoprecipitation experiments (ChIP-chip). The expression data were publicly available; the ChIP-chip data were not. In order to assess the likelihood that a gene is a BCL6 target, we used three classes of information: expression-level differences, over-representation of sequence motifs in promoter regions, and gene ontology annotations. A weight was attached to each analysis based on how well it identified BCL6-bound genes as defined by publicly available ChIP-chip data. By the organizers' criteria, our group, GenomeSingapore, performed best. However, our retrospective analysis indicates that this success was dominated by a gene expression analysis that was predicated on a regulatory model known to be favored by the organizers. We also noted that the 200-gene test set was enriched only in genes that are upregulated, while genes bound by BCL6 are enriched in both upregulated and downregulated genes. Together, these observations suggest possible model biases in the selection of the gold-standard gene set and imply that our success was attained in part by adhering to the same assumptions. We argue that model biases of this type are unavoidable in the inference of regulatory networks and, for that reason, we suggest that future community-wide experiments of this type should focus on the prediction of data, rather than models.

Genomewide location analysis of Candida albicans Upc2p, a regulator of sterol metabolism and azole drug resistance

Upc2p, a transcription factor of the zinc cluster family, is an important regulator of sterol biosynthesis and azole drug resistance in Candida albicans. To better understand Upc2p function in C. albicans, we used genomewide location profiling to identify the transcriptional targets of Upc2p in vivo. A triple hemagglutinin epitope, introduced at the C terminus of Upc2p, conferred a gain-of-function effect on the fusion protein. Location profiling identified 202 bound promoters (P < 0.05). Overrepresented functional groups of genes whose promoters were bound by Upc2p included 12 genes involved in ergosterol biosynthesis (NCP1, ERG11, ERG2, and others), 18 genes encoding ribosomal subunits (RPS30, RPL32, RPL12, and others), 3 genes encoding drug transporters (CDR1, MDR1, and YOR1), 4 genes encoding transcription factors (INO2, ACE2, SUT1, and UPC2), and 6 genes involved in sulfur amino acid metabolism (MET6, SAM2, SAH1, and others). Bioinformatic analyses suggested that Upc2p binds to the DNA motif 5'-VNCGBDTR that includes the previously characterized Upc2p binding site 5'-TCGTATA. Northern blot analysis showed that increased binding correlates with increased expression for the analyzed Upc2p targets (ERG11, MDR1, CDR1, YOR1, SUT1, SMF12, and CBP1). The analysis of ERG11, MDR1, and CDR1 transcripts in wild-type and upc2Delta/upc2Delta strains grown under Upc2p-activating conditions (lovastatin treatment and hypoxia) showed that Upc2p regulates its targets in a complex manner, acting as an activator or as a repressor depending upon the target and the activating condition. Taken together, our results indicate that Upc2p is a key regulator of ergosterol metabolism. They also suggest that Upc2p may contribute to azole resistance by regulating the expression of drug efflux pump-encoding genes in addition to ergosterol biosynthesis genes.