# A comprehensive functional portrait of two heat shock factor-type transcriptional regulators involved in Candida albicans morphogenesis and virulence.

Sadri Znaidi, Audrey Nesseir, Murielle Chauvel, Tristan Rossignol, Christophe d'Enfert

► To cite this version:

Sadri Znaidi, Audrey Nesseir, Murielle Chauvel, Tristan Rossignol, Christophe d'Enfert. A comprehensive functional portrait of two heat shock factor-type transcriptional regulators involved in Candida albicans morphogenesis and virulence. PLoS Pathogens, 2013, 9 (8), pp.e1003519. 10.1371/journal.ppat.1003519. pasteur-01522813

HAL Id: pasteur-01522813

https://pasteur.hal.science/pasteur-01522813

Submitted on 15 May 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Comprehensive Functional Portrait of Two Heat Shock Factor-Type Transcriptional Regulators Involved in *Candida albicans* Morphogenesis and Virulence

Sadri Znaidi¹,², Audrey Nesseir¹,²,³, Murielle Chauvel¹,², Tristan Rossignol¹,²*<sup>¶</sup>, Christophe d'Enfert¹,²*

¹ Institut Pasteur, Unité Biologie et Pathogénicité Fongiques, Département Génomes et Génétique, Paris, France, ²INRA, USC2019, Paris, France, ³ Université Paris Diderot, Sorbonne Paris Cité, Cellule Pasteur, Paris, France

## Abstract

Sfl1p and Sfl2p are two homologous heat shock factor-type transcriptional regulators that antagonistically control morphogenesis in *Candida albicans*, while being required for full pathogenesis and virulence. To understand how Sfl1p and Sfl2p exert their function, we combined genome-wide location and expression analyses to reveal their transcriptional targets *in vivo* together with the associated changes of the *C. albicans* transcriptome. We show that Sfl1p and Sfl2p bind to the promoter of at least 113 common targets through divergent binding motifs and modulate directly the expression of key transcriptional regulators of *C. albicans* morphogenesis and/or virulence. Surprisingly, we found that Sfl2p additionally binds to the promoter of 75 specific targets, including a high proportion of hyphal-specific genes (HSGs; HWP1, HYR1, ECE1, others), revealing a direct link between Sfl2p and hyphal development. Data mining pointed to a regulatory network in which Sfl1p and Sfl2p act as both transcriptional activators and repressors. Sfl1p directly represses the expression of positive regulators of hyphal growth (*BRG1*, *UME6*, *TEC1*, *SFL2*), while upregulating both yeast form-associated genes (*RME1*, *RHD1*, *YWP1*) and repressors of morphogenesis (*SSN6*, *NRG1*). On the other hand, Sfl2p directly upregulates HSGs and activators of hyphal growth (*UME6*, *TEC1*), while downregulating yeast form-associated genes and repressors of morphogenesis (*NRG1*, *RFG1*, *SFL1*). Using genetic interaction analyses, we provide further evidences that Sfl1p and Sfl2p antagonistically control *C. albicans* morphogenesis through direct modulation of the expression of important regulators of hyphal growth. Bioinformatic analyses suggest that binding of Sfl1p and Sfl2p to their targets occurs with the co-binding of Efg1p and/or Ndt80p. We show, indeed, that Sfl1p and Sfl2p targets are bound by Efg1p and that both Sfl1p and Sfl2p associate *in vivo* with Efg1p. Taken together, our data suggest that Sfl1p and Sfl2p act as central "switch on/off" proteins to coordinate the regulation of *C. albicans* morphogenesis.

**Citation:** Znaidi S, Nesseir A, Chauvel M, Rossignol T, d'Enfert C (2013) A Comprehensive Functional Portrait of Two Heat Shock Factor-Type Transcriptional Regulators Involved in *Candida albicans* Morphogenesis and Virulence. PLoS Pathog 9(8): e1003519. doi:10.1371/journal.ppat.1003519

**Editor:** Robert A. Cramer, Geisel School of Medicine at Dartmouth, United States of America

**Received** January 30, 2013; **Accepted** June 8, 2013; **Published** August 15, 2013

**Copyright:** © 2013 Znaidi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by research grants from the European Commission (FinSysB PITN-GA-2008-214004) and Agence Nationale de la Recherche (KANJI, ANR-08-MIE-033-01) to CdE. SZ was supported by postdoctoral fellowships from the European Commission (FinSysB PITN-GA-2008-214004) and the Agence Nationale de la Recherche (KANJI, ANR-08-MIE-033-01). AN was supported by a doctoral studentship from the DIM-Malinf Région Ile-de-France. This study has received funding from the French Government's Investissement d'Avenir program, Laboratoire d'Excellence "Integrative Biology of Emerging Infectious Diseases" (ANR-10-LABX-62-IBEID). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Introduction

*Candida albicans* is the most frequent causative agent of superficial as well as disseminated, life-threatening fungal infections [1]. The success of *C. albicans* as a major fungal pathogen of humans relies on a number of pathogenic traits, among which its capacity to grow and switch between at least three distinctive morphological forms: budding yeast, pseudohyphae and hyphae [2–5]. The morphogenetic transition has been commonly described as a critical trait for survival and virulence in the host, even though the analysis of a wide array of *C. albicans* knock-out mutants suggests that pathogenesis can be dissociated to some extent from morphological switching [6–8].

The yeast-to-hyphae transition is triggered by a variety of environmental stimuli including nutrient availability, temperature, pH, CO₂ and serum [9–13]. This process correlates with the

coordinated expression of a set of hyphal-specific genes (HSGs) with roles in orchestrating hyphal development. Consequently, the transition is highly regulated and involves multiple interconnected signalling pathways, including the cyclic AMP-dependent Protein Kinase A (cAMP-PKA, regarded as playing a central role in the control of morphogenesis), the Cph1p-mediated Mitogen-Activated Protein Kinase (MAPK) and the Rim101p-mediated pH cascade pathways, all of which positively regulate hyphal development through the modulation of the activity of transcription factors to control the expression of HSGs (see [13] for a recent review). These transcription factors include (among others) Efg1p/ Flo8p, acting downstream of cAMP-PKA [14–20], Tec1p [21] and Ume6p [22,23]. Hyphal morphogenesis is also subject to negative regulation mostly by the general corepressor Tup1p through interaction with the transcriptional repressors Nrg1p and Rfg1p [4,12,24–27].

<footnote>* E-mail: christophe.denfert@pasteur.fr</footnote>

<footnote>¶ Current address: Micalis - BimLip, Thiverval-Grignon, France.</footnote>
## Author Summary

Candida albicans can switch from a harmless colonizer of body organs to a life-threatening invasive pathogen. This switch is linked to the ability of C. albicans to undergo a yeast-to-filament shift induced by various cues, including temperature. Sfl1p and Sfl2p are two transcription factors required for C. albicans virulence, but antagonistically regulate morphogenesis: Sfl1p represses it, whereas Sfl2p activates it in response to temperature. We show here that Sfl1p and Sfl2p bind in vivo, via divergent motifs, to the regulatory region of a common set of targets encoding key determinants of morphogenesis and virulence and exert both activating and repressing effects on gene expression. Additionally, Sfl2p binds to specific targets, including genes essential for hyphal development. Bioinformatic analyses suggest that Sfl1p and Sfl2p control C. albicans morphogenesis by cooperating with two important regulators of filamentous growth, Efg1p and Ndt80p, a premise that was confirmed by the observation of concomitant binding of Sfl1p, Sfl2p and Efg1p to the promoter of target genes and the demonstration of direct or indirect physical association of Sfl1p and Sfl2p with Efg1p, in vivo. Our data suggest that Sfl1p and Sfl2p act as central "switch on/off" proteins to coordinate the regulation of C. albicans morphogenesis.

In the yeast *Saccharomyces cerevisiae*, which has been used as a model for studying the transcriptional control of the morphological transition [28,29], Sfl1p (ScSfl1p, for suppressor gene for flocculation 1) is a target of the cAMP-PKA pathway [30]. ScSFL1 encodes a negative regulator of pseudohyphal growth and invasion [31] and was isolated based on its ability to suppress flocculation defects in yeast [32]. ScSfl1p carries a putative heat shock factor (HSF)-type DNA binding domain and binds *in vitro* to a GAA triplet motif [33] characteristic of heat shock elements (HSEs) [34], while exerting its negative regulation through the recruitment of the Ssn6p-Tup1p corepressor complex [35]. ScSfl1p has dual activator/repressor functions, acting as a transcriptional repressor of flocculation-related genes and as an activator of stress-responsive genes [35,36]. Interestingly, the C. albicans genome encodes two structural homologs of ScSfl1p, namely Sfl1p and Sfl2p [37–40]. Either SFL1 or SFL2 functionally complement an *S. cerevisiae* sfl1 mutation [38,39] and encode important regulators of morphogenesis and virulence in C. albicans [37–40]. Intriguingly, although sharing structural homologies, Sfl1p and Sfl2p have antagonistic functions: while Sfl1p acts as a negative regulator of hyphal development, Sfl2p acts as a positive regulator of this process [37–40]. Functional analyses of C. albicans Sfl1p showed that deletion of SFL1 promoted filamentous growth and cell flocculation and correlated with induction of HSGs (*ECE1*, *HWP1*) and genes involved in cell adhesion (*ALS1*, *ALS3*), whereas its overexpression inhibited hyphal formation [37,38]. Consistent with a transcriptional regulatory function, an Sfl1p-GFP fusion localized to the nucleus, while one hybrid *lacZ* reporter analyses in C. albicans correlated with a repressor function [37]. Importantly, either deletion or overexpression of SFL1 attenuated C. albicans virulence in a mouse model of systemic infection [38]. On the other hand, we and others have shown that deletion of SFL2 impaired filamentation in response to different cues, whereas SFL2 overexpression promoted hyphal growth, even under non hyphae-stimulating conditions [39–41]. Noteworthy, an *sfl2Δ/sfl2Δ* strain exhibited reduced damage in a reconstituted human oral epithelium model and displayed attenuated virulence in a mouse model of gastrointestinal colonization and dissemination

model [39,40], indicating that Sfl2p also plays an important role in C. albicans pathogenesis. Similar to Sfl1p, an Sfl2p-GFP fusion localized to the nucleus, in line with a role in transcriptional regulation [39].

It is still unknown how Sfl1p and Sfl2p exert their antagonistic functions. Both SFL1 and SFL2 were shown to genetically interact with at least transcription factor *FLO8*. Hyphal development in *sfl1Δ/sfl1Δ* was abolished upon deletion of *FLO8* but enhanced upon *FLO8* overexpression [38] while overexpression of *SFL2* triggered filamentation in a *FLO8*- and *EFG1*-dependent manner [39], suggesting the implication of the cAMP-PKA pathway. It was also shown that *SFL2* is required for hyphal maintenance at high temperature and that a temperature increase from 25°C to 37°C leads to upregulation of both the RNA and protein levels of Sfl2p, indicating that Sfl2p is a temperature-responsive regulator [39]. In contrast, no clear association was determined between temperature and Sfl1p function. Interestingly, Song *et al.* showed that the putative HSF domains of Sfl1p and Sfl2p were required for their functional divergence by testing HSF domain-swapped hybrids for their ability to retain their effect on filamentation [39]. This suggests that the two putative HSF domains in Sfl1p and Sfl2p mediate the specific recognition of divergent target sites that determine the activation or repression roles of Sfl1p and Sfl2p [39]. To shed more light on Sfl1p and Sfl2p functions, we provide a comprehensive functional portrait of these two regulators using a combination of genome-wide location, genome-wide expression and genetic interaction analyses. We provide evidences that Sfl1p and Sfl2p act as central "switch on-off" proteins to coordinate the regulation of C. albicans morphogenesis and, potentially, pathogenesis and virulence.

## Results

### Epitope-tagging of Sfl1p and Sfl2p

To better characterize the function of Sfl1p and Sfl2p, we sought to identify their DNA-binding locations, *in vivo*, by chromatin immunoprecipitation. To this end, we generated triple-hemagglutinin epitope (HA₃)-tagged versions of SFL1 and SFL2 and used the pCaEXP system [42] to drive *MET3* promoter-dependent expression of the tagged alleles in *sfl1Δ/sfl1Δ* (Table 1; strain *sfl1-CaEXP-SFL1-HA₃*) and *sfl2Δ/sfl2Δ* (Table 1, strain *sfl2-CaEXP-SFL2-HA₃*) mutant strains, respectively (Figure 1A, see Materials and Methods for specific details). We also generated *sfl1Δ/sfl1Δ* and *sfl2Δ/sfl2Δ* mutants carrying the empty pCaEXP vector (*sfl1-CaEXP* and *sfl2-CaEXP*, respectively, see Table 1) to serve as negative controls for immunoprecipitation. Western blot analyses of strains grown under *P*.*MET3*-inducing conditions showed that both Sfl1p-HA₃ and Sfl2p-HA₃ fusion proteins were expressed (Figure 1B, lanes 4 and 6). As an additional control for signal specificity, immunoblotting of total extracts from a C. albicans strain expressing the Cap1p-HA₃ fusion (Figure 1, lane 2) or the corresponding empty-vector negative control (Figure 1, lane 1) was used [43].

To test the functionality of the Sfl1p-HA₃ and Sfl2p-HA₃ fusions, both tagged and empty-vector control strains were grown overnight at 30°C in YPD then transferred to Lee's medium (hyphae-inducing medium) lacking methionine (*P*.*MET3*-inducing condition) at 37°C and allowed to resume growth during 4 h prior to microscopic examination (Figure 1C). It was previously shown that *P*.*MET3*-driven expression of wild-type SFL1 in a homozygous *sfl1* mutant strain under hyphae-inducing conditions abolished hyphal formation [37]. As expected, hyphal formation was induced in either the control strain SC5314 or the *sfl1Δ/sfl1Δ* mutant carrying the empty vector (Figure 1C, top left and middle
<table><caption>**Table 1.** Strains used in this study.</caption><thead><tr><th>Strain name</th><th>Lab identifier</th><th>Parental strain</th><th>Relevant genotype</th><th>Reference</th></tr></thead><tbody><tr><td>SC5314</td><td>CEC1462</td><td></td><td>Prototrophic</td><td>[84]</td></tr><tr><td>CAI4</td><td>CEC2095</td><td>SC5314</td><td>ura3Δ::λimm434/ura3Δ::λimm434</td><td>[85]</td></tr><tr><td>BWP17H</td><td>CEC157</td><td>BWP17</td><td>ura3Δ::λimm434/ura3Δ::λimm434, his1Δ::hisG/HIS1, arg4Δ::hisG/arg4Δ::hisG</td><td>Lab collection</td></tr><tr><td>BWP17AH</td><td>CEC161</td><td>BWP17</td><td>ura3Δ::λimm434/ura3Δ::λimm434, his1Δ::hisG/HIS1, arg4Δ::hisG/ARG4</td><td>[86]</td></tr><tr><td>SN76</td><td>CEC805</td><td></td><td>arg4Δ/arg4Δ, his1Δ/his1Δ, ura3Δ::λimm434/ura3Δ::λimm434, iro1Δ::λimm434/iro1Δ::λimm434</td><td>[87]</td></tr><tr><td>HLC52</td><td>CEC150</td><td></td><td>ura3Δ::λimm434/ura3Δ::λimm434, efg1Δ::hisG/efg1Δ::hisG-URA3-hisG</td><td>[17]</td></tr><tr><td>HLCEEG1</td><td>CEC3891</td><td>CAI4</td><td>ura3Δ::λimm434/ura3Δ::λimm434, efg1Δ::hisG/efg1Δ::EFG1-HA-URA3</td><td>[18]</td></tr><tr><td>AVL12</td><td>CEC3894</td><td>BWP17</td><td>ura3Δ::λimm434/ura3Δ::λimm434, arg4Δ::hisG/arg4Δ::hisG, his1Δ::hisG/his1Δ::hisG, efg1Δ::hisG/efg1Δ::EFG1-HA-URA3</td><td>[18]</td></tr><tr><td>AVL12-SFL1-TAP</td><td>CEC3923</td><td>AVL12</td><td>ura3Δ::λimm434/ura3Δ::λimm434, arg4Δ::hisG/arg4Δ::hisG, his1Δ::hisG/his1Δ::hisG, efg1Δ::hisG/efg1Δ::EFG1-HA-URA3, SFL1/SFL1-TAP-HIS1</td><td>This study</td></tr><tr><td>AVL12-SFL2-TAP</td><td>CEC3916</td><td>AVL12</td><td>ura3Δ::λimm434/ura3Δ::λimm434, arg4Δ::hisG/arg4Δ::hisG, his1Δ::hisG/his1Δ::hisG, efg1Δ::hisG/efg1Δ::EFG1-HA-URA3, SFL2/SFL2-TAP-HIS1</td><td>This study</td></tr><tr><td>AVL12-pHIS</td><td>CEC3913</td><td>AVL12</td><td>ura3Δ::λimm434/ura3Δ::λimm434, arg4Δ::hisG/arg4Δ::hisG, his1Δ::HIS1/his1Δ::hisG, efg1Δ::hisG/efg1Δ::EFG1-HA-URA3</td><td>This study</td></tr><tr><td>SGY243-CaEXP-B</td><td>CEC2894</td><td>SGY243</td><td>RPS1::(pCaEXP) URA3 P<sub>MET3</sub></td><td>[43]</td></tr><tr><td>SGY243-CaEXP-CAP1-HA</td><td>CEC2895</td><td>SGY243</td><td>RPS1::(pCaEXP) URA3 P<sub>MET3</sub>-CAP1-HA<sub>3</sub></td><td>[43]</td></tr><tr><td>CEC1561</td><td>CEC1561</td><td>SN76</td><td>sfl1Δ::ARG4/SFL1</td><td>This study</td></tr><tr><td>SFL1-TAP</td><td>CEC1922</td><td>CEC1561</td><td>sfl1Δ::ARG4/SFL1-TAP-HIS1 RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC1422</td><td>CEC1422</td><td>SN76</td><td>sfl2Δ::ARG4/SFL2</td><td>This study</td></tr><tr><td>SFL2-TAP</td><td>CEC1918</td><td>CEC1422</td><td>sfl1Δ::ARG4/SFL2-TAP-HIS1 RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC3075</td><td>CEC3075</td><td>CEC1561</td><td>sfl1Δ::ARG4/sfl1Δ::SFL1-HA<sub>3</sub>-URA3-HA<sub>3</sub></td><td>This study</td></tr><tr><td>CEC3076</td><td>CEC3076</td><td>CEC1422</td><td>sfl2Δ::ARG4/sfl2Δ::SFL2-HA<sub>3</sub>-URA3-HA<sub>3</sub></td><td>This study</td></tr><tr><td>sfl1Δ/sfl1Δ</td><td>CEC1910</td><td>CEC1561</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1</td><td>This study</td></tr><tr><td>CEC1997</td><td>CEC1997</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, RPS1/rps1::(Clp10) URA3 P<sub>PCKI</sub>-SFL1-TAP</td><td>This study</td></tr><tr><td>sfl1-CaEXP</td><td>CEC3283</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, RPS1/RPS1::(pCaEXP) URA3 P<sub>MET3</sub></td><td>This study</td></tr><tr><td>sfl1-CaEXP-SFL1-HA<sub>3</sub></td><td>CEC3284</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, RPS1/RPS1::(pCaEXP) URA3 P<sub>MET3</sub>-SFL1-HA<sub>3</sub></td><td>This study</td></tr><tr><td>sfl2Δ/sfl2Δ</td><td>CEC1503</td><td>CEC1422</td><td>sfl2Δ::ARG4/sfl2Δ::HIS1</td><td>This study</td></tr><tr><td>sfl2-CaEXP</td><td>CEC3253</td><td>CEC1503</td><td>sfl2Δ::ARG4/sfl2Δ::HIS1, RPS1/RPS1::(pCaEXP) URA3 P<sub>MET3</sub></td><td>This study</td></tr><tr><td>sfl2-CaEXP-SFL2-HA<sub>3</sub></td><td>CEC3254</td><td>CEC1503</td><td>sfl2Δ::ARG4/sfl2Δ::HIS1, RPS1/RPS1::(pCaEXP) URA3 P<sub>MET3</sub>-SFL2-HA<sub>3</sub></td><td>This study</td></tr><tr><td>sfl1ΔΔ</td><td>CEC2001</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>sfl1ΔΔ sfl2ΔΔ</td><td>CEC2658</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, sfl2Δ::URA3/sfl2Δ::SAT1</td><td>This study</td></tr><tr><td>sfl2ΔΔ</td><td>CEC1535</td><td>CEC1503</td><td>sfl2Δ::ARG4/sfl2Δ::HIS1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC1509</td><td>CEC1509</td><td>CEC1503</td><td>sfl2Δ::ARG4/sfl2Δ::HIS1, RPS1/rps1::(Clp10) URA3 P<sub>PCKI</sub>-SFL2-TAP</td><td>This study</td></tr><tr><td>sfl1ΔΔ brg1ΔΔ</td><td>CEC2840</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, brg1Δ::URA3/brg1Δ::SAT1</td><td>This study</td></tr><tr><td>brg1Δ/brg1Δ</td><td>CEC2009</td><td>SN76</td><td>brg1Δ::ARG4/brg1Δ::HIS1</td><td>This study</td></tr><tr><td>brg1ΔΔ</td><td>CEC2058</td><td>CEC2009</td><td>brg1Δ::ARG4/brg1Δ::HIS1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC3485</td><td>CEC3485</td><td>BWP17AH</td><td>ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC2988</td><td>CEC2988</td><td>BWP17AH</td><td>ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3 P<sub>TET</sub>-SFL2</td><td>This study</td></tr><tr><td>CEC3431</td><td>CEC3431</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC3484</td><td>CEC3484</td><td>CEC1910</td><td>sfl1Δ::ARG4/sfl1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3 P<sub>TET</sub>-SFL2</td><td>This study</td></tr><tr><td>CEC3435</td><td>CEC3435</td><td>CEC1503</td><td>sfl2Δ::ARG4/sfl2Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC3437</td><td>CEC3437</td><td>CEC1503</td><td>sfl2Δ::ARG4/sfl2Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3 P<sub>TET</sub>-SFL2</td><td>This study</td></tr><tr><td>ume6Δ/ume6Δ</td><td>CEC2656</td><td>SN76</td><td>ume6Δ::ARG4/ume6Δ::HIS1</td><td>This study</td></tr><tr><td>CEC3583</td><td>CEC3583</td><td>CEC2656</td><td>ume6Δ::ARG4/ume6Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC3585</td><td>CEC3585</td><td>CEC2656</td><td>ume6Δ::ARG4/ume6Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3 P<sub>TET</sub>-SFL2</td><td>This study</td></tr><tr><td>tec1Δ/tec1Δ</td><td>CEC2335</td><td>SN76</td><td>tec1Δ::ARG4/tec1Δ::HIS1</td><td>This study</td></tr><tr><td>CEC3589</td><td>CEC3589</td><td>CEC2335</td><td>tec1Δ::ARG4/tec1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC3591</td><td>CEC3591</td><td>CEC2335</td><td>tec1Δ::ARG4/tec1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3 P<sub>TET</sub>-SFL2</td><td>This study</td></tr><tr><td>CEC3581</td><td>CEC3581</td><td>CEC2009</td><td>brg1Δ::ARG4/brg1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr></tbody></table>
<table><caption>**Table 1.** Cont.</caption><thead><tr><th>Strain name</th><th>Lab identifier</th><th>Parental strain</th><th>Relevant genotype</th><th>Reference</th></tr></thead><tbody><tr><td>CEC3642</td><td>CEC3642</td><td>CEC2009</td><td>brg1Δ::ARG4/brg1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3 P<sub>TET</sub>-SFL2</td><td>This study</td></tr><tr><td>efg1Δ/efg1Δ</td><td>CEC1439</td><td>HLC52</td><td>ura3Δ::λimm434/ura3Δ::λimm434, efg1Δ::hisG/efg1Δ::hisG</td><td>This study</td></tr><tr><td>CEC3581</td><td>CEC3581</td><td>CEC1439</td><td>efg1Δ::ARG4/efg1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3</td><td>This study</td></tr><tr><td>CEC3156</td><td>CEC3156</td><td>CEC1439</td><td>efg1Δ::ARG4/efg1Δ::HIS1, ADH1/adh1::P<sub>TDH3</sub>-carTA::SAT1, RPS1/RPS1::(Clp10) URA3 P<sub>TET</sub>-SFL2</td><td>This study</td></tr></tbody></table>

doi:10.1371/journal.ppat.1003519.t001

panels, respectively). Conversely, hyphal formation was strongly impaired in the strain expressing *SFL1-HA3* (Figure 1C, top right panel), therefore phenocopying the effect of *P<sub>MET3</sub>*-driven wild-type *SFL1* expression as observed in Bauer *et al.* [37]. Under the same growth conditions the *sfl2Δ/sfl2Δ* strain carrying the empty vector was unable to form hyphae (Figure 1C, bottom middle panel), whereas expression of the *SFL2-HA3* allele allowed induction of hyphal formation as observed in strain SC5314 (Figure 1, compare bottom left and right panels). Taken together, these results show that epitope-tagging of Sfl1p and Sfl2p at their C-termini using the pCaEXP system allowed the production of fully functional proteins.

## Genome-wide location map of Sfl1p and Sfl2p at a single nucleotide resolution

We performed genome-wide location of Sfl1p or Sfl2p under hyphae-inducing conditions by chromatin immunoprecipitation coupled to massively parallel high-throughput sequencing (ChIP-Seq, see Materials and Methods), which allows to detect binding events at a single nucleotide resolution. The resulting reads were mapped to the *C. albicans* Assembly 21 genome and alignments were visualized using the Integrative Genomics Viewer (IGV) software [44,45] (see Materials and Methods for details). Using the Model-Based Analysis for ChIP-Seq (MACS) peak-finding algorithm [46], we identified 163 and 213 binding peaks for Sfl1p and

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Strategy for tagging Sfl1p and Sfl2p with a triple hemagglutinin (3×HA) epitope tag and characterization of the tagged strains. (A) Schematic representation of the *SFL1-HA3* or *SFL2-HA3* tagging cassette allowing expression of the Sfl1p-HA3 or Sfl2p-HA3 fusion proteins following a Stul digestion (Stul) and integration at the RPS1 locus (RPS1, black rectangles) [42]. A triple HA tag (dark grey box) was inserted in frame with the *SFL1* or *SFL2* coding sequences (*SFL1* or *SFL2*; black arrowed rectangle) in plasmid pCaEXP [42]. The tagged alleles are placed under the control of the *MET3* promoter (*MET3p*; ligh grey rectangle), which is induced in the absence of methionine and cysteine, and are followed by the *C. albicans* URA3 marker (open rectangle). (B) Western blot analysis of homozygous *sfl1* or *sfl2* mutants (*sfl1Δ/sfl1Δ* or *sfl2Δ/sfl2Δ*) expressing HA3-tagged versions of the *SFL1* or *SFL2* genes, respectively (*SFL1-HA3* or *SFL2-HA3*) together with the corresponding empty vector controls (Vector). The SGY243 strain expressing the *CAP1-HA3* (*CAP1-HA3*) or carrying the empty vector (Vector) were used as a positive control [43]. Strains were grown overnight in SD medium (*P<sub>MET3</sub>*-inducing conditions) and total protein extracts were prepared then subjected to SDS-PAGE. Western blotting was performed using an anti-HA antibody. Positions of the molecular mass standards are indicated on the left (kDa). Immunopositive signals from the Sfl1p-HA3 and Sfl2p-HA3 fusions are indicated with black arrows (C) Phenotypic analysis of the strains expressing the HA3-tagged *SFL1* or *SFL2* alleles. Strain SC5314 (control) together with the homozygous *sfl1* or *sfl2* mutants expressing the *SFL1-HA3* or *SFL2-HA3* alleles (*SFL1-HA3*, *SFL2-HA3*), respectively, or carrying the empty vector (Vector) were grown overnight in YPD at 30°C then transferred to Lee's medium lacking methionine and cysteine and allowed to grow during 4 h at 37°C before being examined microscopically (40× magnification). doi:10.1371/journal.ppat.1003519.g001</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Genome-wide location of *Candida albicans* Sfl1p and Sfl2p, *in vivo*, at a single-nucleotide resolution. (A) Venn diagram of the overlap between Sfl1p and Sfl2p binding targets. All 113 Sfl1p targets are also bound by Sfl2p, while 75 target promoters are Sfl2p-</figcaption></figure>

specific. The total number of Sfl1p or Sfl2p target promoters are indicated between parentheses. Target promoters include those that are clearly associated with given ORFs as well as those that are shared by two ORFs in opposite orientations. (B) A single-nucleotide resolution of Sfl1p and Sfl2p binding at selected *C. albicans* genomic regions *in vivo*. Plotted are read-count signal intensities of HA₃-tagged *SFL1*-*sfl1-CaEXP-SFL1-HA₃* or *SFL2*-*sfl2-CaEXP-SFL2-HA₃* coimmunoprecipitated DNA and the corresponding empty-vector control signals (*sfl1-CaEXP*, *sfl2-CaEXP*, respectively) from merged BAM files of two independent biological replicates. Some read-count signals extend beyond the maximum graduation (not shown) that ranges between 0–500 reads for Sfl1 data (*sfl1-CaEXP* and *sfl1-CaEXP-SFL1-HA₃*) and 0–1000 reads for Sfl2 data (*sfl2-CaEXP* and *sfl2-CaEXP-SFL2-HA₃*). The position of each signal in selected *C. albicans* genomic regions from assembly 21 is shown on the x-axis. The location of each selected region from the corresponding chromosome (Chr) is indicated at the top of each panel (limits are shown between parentheses in base pairs). The orientation of each ORF is depicted by the arrowed black rectangle. (C) Enrichment scores of the Gene Ontology (GO) terms to which are assigned Sfl1p and Sfl2p common (shaded area) or Sfl2p-specific (unshaded area) binding targets. GO term enrichment scores are calculated as the negative value of the log₁₀-transformed *P*-value. The number of genes of each category is shown at the right of each horizontal bar.
doi:10.1371/journal.ppat.1003519.g002

Sfl2p, respectively (see Tables S1–S6 in Text S1, Legends to Supplementary Tables S1–S8 in Text S1 and Materials and Methods for details).

As expected, most of Sfl1p or Sfl2p binding peaks were located at 'intergenic' regions (Tables S1–S6 in Text S1), consistent with a transcriptional regulatory function. Among the 163 Sfl1p binding peaks, 76 clearly associated with individual ORFs, while 34 were located at promoter regions shared by two ORFs in opposite orientations and the remaining 53 peaks were not clearly associated with ORFs. In particular, spurious binding overlapping with highly transcribed regions [47], mostly tRNA-encoding genes, or regions with repeated DNA sequence (Table S3 in Text S1), was observed. Among the 213 Sfl2p binding peaks, 140 clearly associated with unique ORFs, while 54 were located in promoter regions shared by two ORFs in opposite orientations and the remaining 19 peaks were not clearly linked to defined ORFs (Table S6 in Text S1). Additional bona fide Sfl1p (14 peaks) and Sfl2p (28 peaks) binding peaks were not detected by the peak-finding algorithm and were added to our target lists (Tables S3 and S6 in Text S1, see column entitled "comments" and Legends to Supplementary Tables S1–S8 in Text S1). Overall, examination of Sfl1p and Sfl2p binding peaks allowed to identify 113 and 188 target promoters (Figure 1A) including 39 and 56 promoter regions shared by two ORFs, respectively. Interestingly, all 113 Sfl1p targets were also bound by Sfl2p, suggesting functional interactions between the two regulators, while 75 additional targets were specific to Sfl2p (Figure 2A). In many occurrences, Sfl2p binding at promoter regions strongly overlapped with that of Sfl1p (Figure 2B, top panel as an example). In other cases, Sfl2p binding showed partial (Figure 2B, middle panel as an example) or no overlap (Figure 2B, bottom panel as an example) with Sfl1p binding. Noteworthy, Sfl2p and Sfl1p binding peaks were often lying across relatively long regions, particularly in the vicinity of transcription factor-encoding genes such as *EGF1* (Figure 2B, top panel), *UME6*, *NRG1* or *TEC1*, suggesting the presence of more than one binding site or the existence of functional interactions with other regulatory proteins at these sites.

We used the GO Term Finder tool from the CGD [48] to identify functional enrichment among Sfl1p and Sfl2p targets relative to the annotated *C. albicans* genome (Table 2; see Materials and Methods). Strikingly, we found that the most significantly
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Binding of Sfl1p-HA₃ and Sfl2p-HA₃ to selected target promoters. Strains *sfl1-CaEXP-SFL1-HA₃* (Sfl1p-HA₃) and *sfl2-CaEXP-SFL2-HA₃* (Sfl2p-HA₃) together with their respective untagged control strains (Vector) were grown under the same conditions as those for the ChIP-Seq experiment prior to ChIP followed by PCR to detect specific Sfl1p and Sfl2p binding enrichment at selected target promoters (See Materials and Methods for details). PCR was performed using primers corresponding to the promoter region of the indicated genes. The *URA3* and *YAK1* genes were used as a negative control for ChIP enrichment. Primer efficiency (shown on the right panel) was tested by the ability of the corresponding primers to quantify 10-fold serially diluted whole cell extract DNA (WCE, ChIP input samples, dilution factors are indicated at the top of the right panel). doi:10.1371/journal.ppat.1003519.g003</figcaption></figure>

enriched functional category among Sfl1p and Sfl2p common
targets was "Sequence-specific DNA-binding transcription activ-
ity" (21 genes, P=1.42×10⁻⁸; Figure 2C, grey shading), including
a large number of genes encoding major transcription factors
involved in C. albicans morphogenesis and virulence such as
UME6, TEC1, NRG1, RFG1, BRG1, FLO8, others (Figure 2C and
Table 2). In line with this finding, the functional grouping
"Filamentous growth" (30 genes, P=1.83×10⁻⁶) was also among
the most overrepresented categories of the identified GO terms
and included the above-mentioned transcription factors in
addition to HMS1, encoding a transcription factor that controls
C. albicans morphogenesis mediated by HSP90 compromise or high
temperature [49], as well as many genes encoding effectors or
signal transducers of this process such as MSB2, CHT2, GAP1,
ALS1, RAS2, others (Figure 2C). As expected, "Pathogenesis" (16
genes; P=2.40×10⁻⁴) was also among the most significantly
enriched functional categories among Sfl1p and Sfl2p common
targets and is consistent with the known roles of Sfl1p and Sfl2p in
C. albicans virulence [38,39]. Interestingly, Sfl1p and Sfl2p bound

to genes encoding transcription factors involved in white/opaque
switching, including WOR2, FLO8, EFG1 and AHR1 (“Regulation
of phenotypic switching”; 4 genes; P = 4.34×10⁻²), as well as
genes involved in biofilm formation (“Biofilm formation”; 12
genes; P = 6.40×10⁻⁴), suggesting wider functions for these two
regulators in C. albicans. These functions may include the ability to
respond to a variety of stimuli, such as drug treatment (“Cellular
response to drug”; 17 genes; P = 2.48×10⁻³), nutrient availability
(“Cellular response to nutrient levels”; 18 genes; P = 3.00×10⁻³
and “Galactose catabolic process via UDP-Galactose”; 3 genes;
P = 2.23×10⁻³) and pH levels (“Cellular response to pH; 9 genes;
P = 3.62×10⁻³).

We also performed functional category enrichment analyses of
the 75 Sfl2p-specific targets (Figure 2C, unshaded area). Interest-
ingly, these targets were grouped into functional categories
pertaining to interaction with the host, including "Multi-organism
process" (19 genes; P = 2.06×10⁻⁵), "Symbiosis, encompassing
mutualism through parasitism" (9 genes; P = 2.18×10⁻³), "Adhe-
sion to host" (6 genes; P = 2.69×10⁻³) and "Fungal-type cell wall"
<table><caption>**Table 2.** Overrepresented functional categories in Sfl1p and Sfl2p ChIP-Seq data.</caption><thead><tr><th>GO term<sup>a</sup></th><th>CGD accession # (ontology classification)<sup>b</sup></th><th>% Frequency<sup>c</sup> (# of genes)</th><th>% Genome frequency<sup>d</sup> (# of genes)</th><th>P value<sup>e</sup></th><th>Genes<sup>f</sup></th></tr></thead><tbody><tr><td colspan="6"><strong>Sfl1p and Sfl2p common targets</strong></td></tr><tr><td>Sequence-specific DNA binding transcription factor activity</td><td>GO:0003700</td><td>19.1 (21)</td><td>3.5 (230)</td><td>0.0000</td><td>FLO8, UME6, FGR15, CRZ2, RFG1, SEF1, SFL2, BRG1, MIG1, RME1, STP2, TEC1, ZCF31, WOR2, EFG1, CUP9, FCR1, NRG1, BCR1, CTA4, AHR1</td></tr><tr><td>Filamentous growth</td><td>GO:0030447</td><td>27.3 (30)</td><td>8.4 (550)</td><td>0.0000</td><td>FLO8, MSB2, UME6, FGR15, RFG1, GAL10, SEF1, CHT2, SFL2, BRG1, GAP1, orf19.4459, STP2, ALS1, RAS2, TEC1, WOR2, RHB1, EFG1, CUP9, TCC1, SSN6, FCR1, orf19.6874, NRG1, BCR1, CTA4, AHR1, AAF1, HMS1</td></tr><tr><td>Regulation of transcription, DNA-dependent</td><td>GO:0006355</td><td>26.4 (29)</td><td>9.2 (601)</td><td>0.0001</td><td>FLO8, UME6, FGR15, CRZ2, RFG1, GAL1, SEF1, SFL2, CTA24, BRG1, MIG1, RME1, STP2, TEC1, ZCF31, WOR2, EFG1, CUP9, TCC1, SSN6, FCR1, orf19.6874, NRG1, BCR1, CTA4, AHR1, HAP41, AAF1, HMS1</td></tr><tr><td>Response to stimulus</td><td>GO:0050896</td><td>40.9 (45)</td><td>19.8 (1290)</td><td>0.0001</td><td>FLO8, MSB2, UME6, HNM1, REG1, FGR15, SIT1, CRZ2, orf19.2726, orf19.2822, RFG1, GSC1, DIP5, GAL1, GAL10, GAL102, SEF1, CHT2, SFL2, BRG1, FET3, FET34, MIG1, orf19.4459, SWE1, orf19.4883, STP2, MDR1, ALS1, RAS2, TEC1, ZCF31, RHB1, EFG1, HSP104, TCC1, SSN6, FCR1, GAC1, NRG1, BCR1, CTA4, AHR1, GPX2, HMS1</td></tr><tr><td>Pathogenesis</td><td>GO:0009405</td><td>14.5 (16)</td><td>3.3 (215)</td><td>0.0002</td><td>FLO8, UME6, FRG1, GSC1, SFL2, BRG1, SWE1, MDR1, ALS1, TEC1, EFG1, HSP104, TCC1, SSN6, NRG1, AHR1</td></tr><tr><td>Biofilm formation</td><td>GO:0042710</td><td>10.9 (12)</td><td>2.0 (128)</td><td>0.0006</td><td>FLO8, CRZ2, YWP1, BRG1, ALS1, TEC1, ZCF31, EFG1, HSP104, NRG1, BCR1, AHR1</td></tr><tr><td>Galactose catabolic process via UDP-galactose</td><td>GO:0033499</td><td>2.7 (3)</td><td>0.0 (3)</td><td>0.0022</td><td>GAL1, GAL10, GAL7</td></tr><tr><td>Cellular response to drug</td><td>GO:0035690</td><td>15.5 (17)</td><td>4.4 (287)</td><td>0.0025</td><td>HNM1, SIT1, GSC1, DIP5, GAL102, FET3, MIG1, SWE1, STP2, MDR1, ZCF31, RHB1, EFG1, SSN6, FCR1, NRG1, AHR1</td></tr><tr><td>Cellular response to nutrient levels</td><td>GO:0031669</td><td>16.4 (18)</td><td>5.0 (323)</td><td>0.0030</td><td>UME6, REG1, FGR15, orf19.2822, RFG1, GAL1, GAL10, CHT2, BRG1, MIG1, orf19.4459, RAS2, RHB1, EFG1, FCR1, NRG1, BCR1, AHR1</td></tr><tr><td>Cellular response to pH</td><td>GO:0071467</td><td>8.2 (9)</td><td>1.2 (81)</td><td>0.0036</td><td>UME6, CRZ2, SEF1, SFL2, STP2, ALS1, EFG1, TCC1, NRG1</td></tr><tr><td>Regulation of flocculation</td><td>GO:0060256</td><td>3.6 (4)</td><td>0.2 (10)</td><td>0.0071</td><td>FLO8, GAL10, SFL2, ALS1</td></tr><tr><td>Regulation of phenotypic switching</td><td>GO:1900239</td><td>3.6 (4)</td><td>0.2 (15)</td><td>0.0434</td><td>FLO8, WOR2, EFG1, AHR1</td></tr><tr><td colspan="6"><strong>Sfl2p-specific targets</strong></td></tr><tr><td>Fungal-type cell wall</td><td>GO:0009277</td><td>15.1 (11)</td><td>2.2 (142)</td><td>0.0000</td><td>HWP1, EAP1, ALS3, PIR1, HYR1, SIM1, RBR3, PGA31, RHD3, WSC1, ALS6</td></tr><tr><td>Multi-organism process</td><td>GO:0051704</td><td>26.0 (19)</td><td>6.4 (418)</td><td>0.0000</td><td>CPH2, HWP1, EAP1, ALS3, CZF1, FCR3, ECE1, SFL1, RFX2, HYR1, ROB1, RHD3, SAP4, SRR1, ADE2, HGC1, RBT4, ALS6, SAP7</td></tr><tr><td>Symbiosis, encompassing mutualism through parasitism</td><td>GO:0044403</td><td>12.3 (9)</td><td>1.9 (126)</td><td>0.0022</td><td>CPH2, HWP1, EAP1, ALS3, RFX2, HYR1, SAP4, HGC1, ALS6</td></tr><tr><td>Biofilm formation</td><td>GO:0042710</td><td>12.3 (9)</td><td>2.0 (128)</td><td>0.0025</td><td>HWP1, EAP1, ALS3, CZF1, FCR3, ECE1, HYR1, ROB1, ALS6</td></tr><tr><td>Sequence-specific DNA binding transcription factor activity</td><td>GO:0003700</td><td>15.1 (11)</td><td>3.5 (230)</td><td>0.0027</td><td>CPH2, orf19.1604, ECM22, CZF1, FCR3, orf19.3328, GRF10, orf19.4342, SFL1, RFX2, ROB1</td></tr><tr><td>Adhesion to host</td><td>GO:0044406</td><td>8.2 (6)</td><td>0.7 (47)</td><td>0.0027</td><td>HWP1, EAP1, ALS3, RFX2, HYR1, SAP4</td></tr></tbody></table>

<footnote><sup>a</sup>Grouping of the Sfl1p and/or Sfl2p targets identified in ChIP-Seq data according to GO terminology determined by using the online CGD GO Term Finder tool (http://www.candidagenome.org/cgi-bin/GO/goTermFinder). Analysis conducted in October 2012 (See Materials and Methods).</footnote>

<footnote><sup>b</sup>Ontology classification according to the three GO terminologies (biological process, cellular component and molecular function).</footnote>

<footnote><sup>c</sup>Percentages were calculated based on the number of genes in each GO category divided by the total number (110 genes for Sfl1p and Sfl2p common targets, 73 genes for Sfl2p specific targets, see Materials and Methods for details).</footnote>

<footnote><sup>d</sup>Percentages were calculated based on the number of genes in each category divided by the total number of annotated genes of the C. albicans genome, according to CGD (6,513 genes).</footnote>

<footnote><sup>e</sup>P values for the overrepresented categories were calculated using a hypergeometric distribution with multiple hypothesis correction (i.e., Bonferroni's correction) as described in the GO Term Finder tool website (http://www.candidagenome.org/help/goTermFinder.shtml). The P value cutoff used was ≈0.05.</footnote>

<footnote><sup>f</sup>Gene name or orf19 nomenclature according to CGD. Some genes were attributed to more than one GO term.
doi:10.1371/journal.ppat.1003519.t002</footnote>
A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /><figcaption>**Figure 4. Sfl1p and Sfl2p transcriptomics.** (A) GeneSpring expression profile plots of each of the three biological replicates from the *sfl1-CaEXP-SFL1-HA<sub>3</sub>* versus *sfl1-CaEXP* (*sfl1-CaEXP-SFL1-HA<sub>3</sub>* vs. *sfl1-CaEXP*) and the *sfl2-CaEXP-SFL2-HA<sub>3</sub>* versus *sfl2-CaEXP* (*sfl2-CaEXP-SFL2-HA<sub>3</sub>* vs. *sfl2-CaEXP*) transcriptomics data. The log<sub>2</sub>-transformed relative expression level of each gene from averaged signal intensities of two nonoverlapping gene-specific microarray probes (See Materials and Methods for details), is shown on the y-axis and the corresponding biological replicate sample for each condition (1, 2 and 3) is shown on the x-axis. The profile plot is coloured according to the ratio observed for replicate 1 in the *sfl1-CaEXP-SFL1-HA<sub>3</sub>* vs. *sfl1-CaEXP* condition. (B) Heat maps of the 30 highest log<sub>2</sub>-transformed relative gene expression levels in the *sfl1-CaEXP-SFL1-HA<sub>3</sub>* versus *sfl1-CaEXP* (*sfl1-CaEXP-SFL1-HA<sub>3</sub>* vs *sfl1-CaEXP*, left panels, UP and DWN) and the *sfl2-CaEXP-SFL2-HA<sub>3</sub>* versus *sfl2-CaEXP* (*sfl2-CaEXP-SFL2-HA<sub>3</sub>* vs *sfl2-CaEXP*, right panels, UP and DWN) transcriptomics data (combination of the 3 biological replicates in each condition). The most upregulated (UP, descending signal intensity) or downregulated (DWN, ascending signal intensity) genes in *sfl1-CaEXP-SFL1-HA<sub>3</sub>* vs. *sfl1-CaEXP* (left panels, *SFL1* column) or *sfl2-CaEXP-SFL2-HA<sub>3</sub>* vs. *sfl2-CaEXP* (*SFL2*, right panels) transcriptomics data and their matching probe intensities from the *sfl2-CaEXP-SFL2-HA<sub>3</sub>* vs. *sfl2-CaEXP* condition (left panels, *SFL2* column) or the *sfl1-CaEXP-SFL1-HA<sub>3</sub>* vs. *sfl1-CaEXP* (right panels, *SFL1* column), respectively, are indicated with their corresponding name or *orf19* nomenclature. Heat maps were constructed using Genesis version 1.7.6 [83]. doi:10.1371/journal.ppat.1003519.g004</figcaption></figure>

(11 genes; *P* = 1.92×10<sup>-5</sup>). Sfl2p also bound specifically to 11 genes encoding transcription factors such as *CPH2*, *ECM22*, *CZF1*, *FCR3*, *RFX2* and *ROB1* (Table 2). We also found that Sfl2p bound specifically to the *SFL1* promoter, while both Sfl1p and Sfl2p bound to the promoter of *SFL2*, suggesting an autoregulatory loop controlling *SFL2* expression.

To validate our ChIP-Seq data, we performed additional independent ChIP experiments and measured Sfl1p and Sfl2p binding by PCR (ChIP-PCR) on selected targets (Figure 3). The *URA3* and *YAK1* genes were used as negative controls for ChIP enrichment. As expected, Sfl1p and Sfl2p binding was detected at the promoter of their targets, including *BRG1*, *EFG1*, *SFL2*, *UME6* and *TEC1* (Figure 3). The promoter region of Sfl2p-specific targets was also enriched by Sfl2p-HA<sub>3</sub> immunoprecipitation, including *SFL1*, *RBT1* and *FAV2*, but not by the immunoprecipitation of Sfl1p-HA<sub>3</sub> (Figure 3).

Taken together, our results suggest that Sfl1p and Sfl2p regulate *C. albicans* morphogenesis and potentially confer virulence through direct binding to the promoter of genes encoding key regulators of these processes. They also revealed that, while both transcription factors bind to common targets, Sfl2p specifically binds to additional target genes that appear to be involved in processes pertaining to interaction with the host.

## Global gene expression profiling reflects the antagonistic functions of *SFL1* and *SFL2* in regulating *C. albicans* morphogenesis

To determine whether Sfl1p and Sfl2p binding targets were also transcriptionally modulated, we performed global gene expression analyses of strains *sfl1-CaEXP-SFL1-HA<sub>3</sub>* versus *sfl1-CaEXP* and *sfl2-CaEXP-SFL2-HA<sub>3</sub>* versus *sfl2-CaEXP* grown 3 times independently under the same conditions than those in the ChIP-Seq experiments (see Materials and Methods for details). We found 643 upregulated and 579 downregulated genes (expression fold-change ≥1.5; *P*≤0.05) in the *sfl1-CaEXP-SFL1-HA<sub>3</sub>* strain as compared to strain *sfl1-CaEXP* (Table S7 in Text S1). On the other hand, 354 genes were upregulated and 478 genes were downregulated (expression fold-change ≥1.5; *P*≤0.05) in strain *sfl2-CaEXP-SFL2-HA<sub>3</sub>* relative to *sfl2-CaEXP* (Table S8 in Text S1). Data were visualized using an expression profile plot (GeneSpring version 12, Agilent Technologies), which allows to get a global view of gene expression variation and thus to compare the expression patterns in *SFL1* and *SFL2* data sets (Figure 4A). Interestingly, most of the highly upregulated genes in *pCaEXP-SFL1-HA<sub>3</sub>* vs. *pCaEXP* data were strongly downregulated in *pCaEXP-SFL2-HA<sub>3</sub>* vs. *pCaEXP* data (Figures 4A and 4B left panel). Many of these genes are
A

<figure><img src="image_7.png" /></figure>

B

<figure><img src="image_8.png" /></figure>

**Figure 5.** Quantitative real-time RT-PCR analysis of selected genes from *SFL1* and *SFL2* transcriptomics data. (A) Expression of the *SFL1*, *SFL2* and *BRG1* genes was quantified by qRT-PCR experiments in *SFL1*-deficient strains carrying or not a functional *SFL1-TAP* fusion [41] and grown during 0, 2 and 4 hours under gluconeogenic conditions (CEC2001 and CEC1997, respectively, Table 1). Expression of the *SFL2* and *BRG1* genes is repressed upon *SFL1* expression. Logarithmic scales are shown in each panel. (B) Expression of the *SFL2*, *SFL1*, *UME6* and *ALS3* genes was quantified by qRT-PCR experiments in an *SFL2*-deficient strain carrying or not a functional *SFL2-TAP* fusion [41] and grown during 0, 2 and 4 hours under gluconeogenic conditions (CEC1509 and CEC1535, respectively). Expression of the *SFL1* gene is repressed at time point 2 h, whereas those of *UME6* and *ALS3* are induced. Logarithmic scales are shown in each panel. Bars in each graph indicate log-transformed relative changes in RNA expression of the indicated samples as compared to the *CEF3* calibrator control (see Materials and Methods). Error bars denote standard deviations.
doi:10.1371/journal.ppat.1003519.g005

markers of the yeast form growth phase, such as *RME1*, *YWP1*, *RHD1* and *orf19.557*. On the other hand, most of the strongly downregulated genes in pCaEXP-*SFL1-HA*<sub>3</sub> vs. pCaEXP data were actually upregulated in pCaEXP-*SFL2-HA*<sub>3</sub> vs. pCaEXP data (Figure 4A), including the HSGs *ECE1*, *ALS3*, *IHD1*, *HWP1*, *HYR1* and *SAP5* (Figure 4B). Examination of the genes that were strongly modulated in pCaEXP-*SFL2-HA*<sub>3</sub> vs. pCaEXP data also revealed similar gene expression dynamics: many of the upregulated genes were found to be downregulated in pCaEXP-*SFL1-HA*<sub>3</sub> vs. pCaEXP data sets, and *vice versa* (Figure 4B, right panel).

We independently confirmed the microarray data by qRT-PCR analyses of selected genes using homozygous *sfl1* or *sfl2* mutant strains expressing (or not) functional TAP (tandem affinity purification)-tagged *SFL1* or *SFL2* alleles [41], respectively, under the control of the *PCK1* promoter (Figure 5, Table 1). Strains were grown under gluconeogenic (*PCK1*-inducing) conditions during 0, 2 and 4 hours and total RNA was isolated followed by qRT-PCR (See Materials and Methods for details). As expected, expression of *SFL1-TAP* gradually increased from time points 0 h to 4 h (Figure 5A, left panel). This increased *SFL1* expression correlated with decreased *SFL2* and *BRG1* expression (Figure 5A, middle and right panels), consistent with a negative regulation of *SFL2* and

*BRG1* expression. On the other hand, *P*<sub>*PCK1*</sub>-induced *SFL2-TAP* expression (Figure 5B, left panel) correlated with decreased expression of *SFL1* (Figure 5B, *SFL1* panel) and increased expression of *UME6* and *ALS3* (Figure 5B, *UME6* and *ALS3* panels), consistent with our microarray data (Figure 4).

Taken together, our transcriptomics data reflect the antagonistic functions of Sfl1p and Sfl2p in regulating *C. albicans* morphogenesis, with *SFL1* promoting the yeast-form growth which correlates with upregulation of yeast-form-specific genes and downregulation of HSGs, and *SFL2* promoting hyphal growth which correlates with upregulation of HSGs and downregulation of yeast-form-specific genes.

## The Sfl1p and Sfl2p regulatory network

We combined the transcriptomics and the ChIP-Seq data in order to get a genome-wide view of the transcriptional modules associated with Sfl1p and Sfl2p regulatory functions (Figure 6). We were expecting to find a substantial amount of genes that are bound by Sfl1p and downregulated in pCaEXP-*SFL1-HA*<sub>3</sub> vs. pCaEXP microarray data, as Sfl1p is thought to act as a repressor. In line with the function of Sfl2p as an activator of hyphal growth, we were also hypothesizing that binding of Sfl2p to its targets would correlate with increased expression of these target genes.
A

<figure><img src="image_9.png" /></figure>

B

<figure><img src="image_10.png" /><figcaption>**Figure 6. Sfl1p and Sfl2p transcriptional modules.** Venn diagrams of the overlap between the genes that are modulated in (A) SFL1 or SFL2 transcriptomics (light red circles, upregulated; light green circles, downregulated; gene expression fold-change cut-off ≥1.5; P-value cut-off ≤0.05) and commonly bound by Sfl1p and Sfl2p (light blue circle) or (B) SFL2 transcriptomics (light red circle, upregulated; light green circle, downregulated;</figcaption></figure>
gene expression fold-change cut-off ≥1.5; P-value cut-off ≤0.05) and specifically bound by Sfl2p (light grey circle). Numbers in the Venn diagrams indicate the number of genes. Circled numbers indicate the number of genes that are (A) both modulated in *SFL1* or *SFL2* transcriptomics data and commonly bound by Sfl1p and Sfl2p or (B) both modulated in *SFL2* transcriptomics data and specifically bound by Sfl2p. The name of these genes (or their orf19 nomenclature) and the functional categories to which they belong are shown in the linked boxes. *, *DCK1* is required for hyphal formation; orf19.3475 is a hyphal induced gene.
doi:10.1371/journal.ppat.1003519.g006

Surprisingly, among the 113 targets commonly bound by Sfl1p and Sfl2p, 40 genes were upregulated and only 22 genes were downregulated in pCaEXP-*SFL1-HA*<sub>3</sub> vs. pCaEXP data (Figure 6A). Conversely, 39 genes were downregulated in pCaEXP-*SFL2-HA*<sub>3</sub> vs. pCaEXP data and only 15 genes were upregulated (Figure 6A), indicating that Sfl1p and Sfl2p have dual transcriptional regulatory functions; acting as both transcriptional activators and transcriptional repressors.

As Sfl1p and Sfl2p respectively act as a repressor and an activator of hyphal growth, we examined the set of genes that were commonly bound by these two regulators and whose expression was both downregulated by *SFL1* and upregulated by *SFL2*. We found 9 genes matching these criteria (Figure 6A, middle right box), among which the key regulators of hyphal growth *UME6* and *TEC1*.

We also examined the set of genes that were both bound by Sfl1p and Sfl2p and upregulated in pCaEXP-*SFL1-HA*<sub>3</sub> vs. pCaEXP and/or downregulated in pCaEXP-*SFL2-HA*<sub>3</sub> vs. pCaEXP microarray data (Figure 6A, left boxes). This is consistent with Sfl1p acting as a transcriptional activator for these genes and/or Sfl2p functioning as their transcriptional repressor. Interestingly, we found that many of these genes encode (or are predicted to encode, e.g. orf19.6874) negative regulators of hyphal growth, including *SSN6*, orf19.6874 [50], *NRG1* and *RFG1* (Figure 6A, left boxes). Of particular interest, *EFG1*, the major regulator of *C. albicans* morphogenesis that functions as both a transcriptional activator and a repressor depending on the growth condition [51] was found to be upregulated by Sfl1p but not modulated in *SFL2* microarray data.

Sfl1p and Sfl2p also bound to the promoter of *BRG1*, *AHRI*, *HMS1* and *SFL2* (Figure 6A), all encoding transcriptional activators of hyphal growth. The expression of *BRG1* and *AHRI* was downregulated by Sfl1p but not modulated by Sfl2p (Figure 6A, bottom right box), whereas the expression of *HMS1* was downregulated by Sfl2p but not modulated by Sfl1p (Figure 6A, bottom left box). Interestingly, Sfl1p binding to the *SFL2* promoter correlates with decreased expression of *SFL2*, indicating a direct negative regulation of *SFL2* expression by Sfl1p (Figures 5A and 6A).

Sfl2p binding to its 75 specific target genes correlated with increased and decreased expression of 24 and 25 genes, respectively (Figure 6B). Strikingly, a significant subset of the genes that are both bound and transcriptionally induced by Sfl2p were the HSGs *ALS3*, *HGC1*, *HWP1*, *HTR1*, *ECE1*, *SAP4*, *IHD1*, *FAV2* and *RBT4* in addition to *DCK1* encoding a putative guanine nucleotide exchange factor required for filamentous growth and the hyphal induced gene orf19.3475 (Figure 6B, upper box). Moreover, Sfl2p directly upregulated genes encoding (or predicted to encode) transcription factors, including *FCR3*, encoding a positive regulator of *C. albicans* adherence [52], orf19.217, encoding a positive regulator of hyphal growth [41] and *RFX2*, encoding a regulator of DNA damage response, adhesion and virulence [53]. On the other hand, Sfl2p directly downregulated the expression of transcription factors *SFL1*, *ECM22*, *ROB1*, encoding a regulator of biofilm formation [54], and many genes involved or predicted to be involved in cell wall integrity (*EAP1*, *FUN31*, *SIM1*, *PIR1* and *RHD3*) as well as genes encoding or

predicted to encode permeases or transporters (*PHO86*, putative inorganic phosphate transporter; *HGT1*, high-affinity glucose transporter; *FLC3*, putative heme transporter; *HIP1* and orf19.7566, putative amino acid transporters).

Taken together, combination of the ChIP-Seq and the transcriptomics data i) indicate that Sfl1p and Sfl2p have dual transcriptional regulatory functions, acting as both activators and repressors, ii) suggest that Sfl1p and Sfl2p antagonistic functions in regulating hyphal morphogenesis is mediated through direct transcriptional modulation of genes encoding key regulators of *C. albicans* morphogenesis, iii) show that Sfl2p additionally specifically controls the expression of HSGs and iv) reveal a direct *SFL1-SFL2* cross-factor negative control.

## *SFL1* and *SFL2* genetically interact with transcriptional targets encoding major regulators of morphogenesis and virulence

Our finding that Sfl1p and Sfl2p directly control the expression of master regulators of *C. albicans* morphogenesis and virulence fostered us to assess the genetic interactions between *SFL1*, *SFL2* and these target genes. Data mining of our ChIP-Seq and transcriptomics results showed that Sfl1p directly negatively regulates *SFL2* expression (Figures 3, 5A and 6A). Moreover, Sfl1p directly negatively regulates the expression of *BRG1* (Figures 3, 5A and 6A), encoding a major regulator of hyphal growth. This suggests that *SFL1* represses filamentation through, at least, direct transcriptional repression of the *SFL2* and *BRG1* genes. To test this hypothesis, we constructed *sfl1Δ/sfl1Δ*, *sfl2Δ/sfl2Δ* and *sfl1Δ/sfl1Δ*, *brg1Δ/brg1Δ* double mutants and tested their ability to form hyphae (Figure 7A). All strains displayed yeast-form growth in SD medium at 30°C (Figure 7A, upper panels). In YP 10% FBS medium at 30°C (Figure 7A, middle and lower panels), which induces moderate filamentation, the homozygous *sfl1* mutant displayed highly dense cell aggregates of a mixture of hyphae and long pseudohyphae (Figure 7A, middle and lower panels), consistent with the function of *SFL1* as a transcriptional repressor of filamentous growth. Interestingly, deletion of *SFL2* or *BRG1* in the *sfl1* mutant strongly reduced filamentous growth as well as cell aggregation (Figure 7A, middle and lower panels), with the *sfl1 sfl2* double mutant cells growing as both yeast form and long to medium-size pseudohyphae and the *sfl1 brg1* double mutants growing as both yeast form and short pseudohyphae (Figure 7A, middle and lower panels). Single homozygous *sfl2* and *brg1* mutants showed phenotypes that were similar to those of the parental wild-type cells (Figure 7A, middle and lower panels).

We showed that Sfl2p directly upregulated *UME6* and *TEC1* expression (Figures 3, 5B and 6A), while specifically directly downregulating the expression of *SFL1* (Figures 3, 5B and 6B), suggesting that *SFL2* controls hyphal induction through at least *UME6*, *TEC1* and *SFL1*. We tested the effect of overexpressing *SFL2* on *C. albicans* morphogenesis in strains carrying the single homozygous deletions *sfl1*, *sfl2*, *ume6*, *tec1*, *brg1* and *efg1* (Figure 7B). We and others previously showed that *SFL2* overexpression in non-hypha-inducing conditions promotes hyphal growth [39,40]. We used the pNIMX system [41] to drive high levels of *SFL2* expression in the above-mentioned strain backgrounds grown in rich medium (Figure 7B). Overexpression of *SFL2* in the wild-type
A

<figure><img src="image_11.png" /></figure>

B

<figure><img src="image_12.png" /></figure>
**Figure 7. Genetic interactions of SFL1 and SFL2 with their transcriptional target genes encoding key regulators of hyphal development. (A)**
The wild-type SC5314 (WT) together with the homozygous *sfl1* (*sfl1ΔΔ*, CEC2001), *sfl2* (*sfl2ΔΔ*, CEC1535), *brg1* (*brg1ΔΔ*, CEC2058), the double homozygous *sfl1*, *sfl2* (*sfl1ΔΔ sfl2ΔΔ*, CEC2658) and *sfl1*, *brg1* (*sfl1ΔΔ brg1ΔΔ*, CEC2840) mutants were grown in yeast-promoting (SD at 30°C for 6 h 30 min) or sub-hypha-inducing (YP 10% FBS at 30°C for 6 h 30 min) conditions and observed microscopically. Scale bar = 10 µm. The detailed cell morphology of each strain grown in YP 10% FBS are shown (Morphological details, bottom panel) (B) The pNIMX expression system [41] was used to drive anhydrotetracycline-dependent overexpression of SFL2 (*P*<sub>TET</sub>-SFL2) in a wild-type (WT, BWP17AH complemented for uracil auxotrophy) or in different homozygous mutant backgrounds, including *sfl1Δ/sfl1Δ* (*sfl1ΔΔ*), *sfl2Δ/sfl2Δ* (*sfl2ΔΔ*), *ume6Δ/ume6Δ* (*ume6ΔΔ*), *tec1Δ/tec1Δ* (*tec1ΔΔ*), *brg1Δ/brg1Δ* (*brg1ΔΔ*) and *efg1Δ/efg1Δ* (*efg1ΔΔ*) (Table 1). All strains were grown in YPD medium at 30°C during 18 hours in the presence of 3 µg/ml of anhydrotetracycline before microscopic examination. As a control, the same growth conditions were also used with all strain backgrounds carrying the empty plasmid (Clp10, Control). Two different fields with detailed cell morphology of each strain overexpressing SFL2 are shown (Morphological details, right panels).
doi:10.1371/journal.ppat.1003519.g007

strain strongly induced filamentation, with cells displaying long pseudohyphae (Figure 7B, top panels). Interestingly, SFL2-driven filamentation was increased in the *sfl1Δ/sfl1Δ* mutant, as compared to that in the wild-type or the *sfl2Δ/sfl2Δ* strains (Figure 7B, compare the zoomed-out regions in lower left corners). Most of the *sfl1* mutant cells overexpressing SFL2 formed longer hyphae and pseudohyphae than those observed in the equivalent *sfl2* mutants (Figure 7B), suggesting that Sfl2p induces filamentous growth in part through repression of SFL1 expression. Conversely, filamentation was strongly reduced in the *ume6Δ/ume6Δ* strain, moderately reduced in either the *tec1Δ/tec1Δ* or *brg1Δ/brg1Δ* mutants and abolished in the *efg1Δ/efg1Δ* strain (Figure 7B). The *ume6* mutants overexpressing SFL2 formed significantly shorter pseudohyphae than those of the equivalent *tec1* and *brg1* mutants (Figure 7B).

Taken together, our results suggest that Sfl1p represses filamentation through at least direct negative regulation of SFL2 and BRG1 expression and indicate that Sfl2p regulates hyphal growth partly through *UME6*, *TEC1* and *BRG1* and totally through *EFG1*.

Motif discovery analyses suggest functional interactions
between Sfl1p, Sfl2p, Efg1p and Ndt80p

Many observations support the hypothesis that Sfl1p and Sfl2p
recognize different binding motifs. First, although sharing com-
mon transcriptional targets, Sfl1p and Sfl2p peak signals are
distributed differently along many of their common target
promoters (Figure 2B, middle panel as an example). Second,
Sfl2p binds specifically to the promoter of 75 targets (Figure 2B,
bottom panel as an example). Third, recent data by Song et al.
suggested that Sfl1p and Sfl2p mediate their functional divergence
through their HSF-type DNA binding domain [39], suggesting
divergent binding sites.

We performed motif-enrichment analyses using DNA sequences encompassing ±250 bp around peak summits in Sfl1p (Figure 8A) and Sfl2p (Figure 8B) binding data. Two independent motif discovery algorithms, the RSA-tools (RSAT) peak-motifs (http://rsat.ulb.ac.be/rsat/, [55]) and SCOPE (genie.dartmouth.edu/scope/, [56]) were used (See Materials and Methods for details). Strikingly, the highest scoring motifs in Sfl1p-enriched sequences included the Ndt80p (5'-ttACACAAA-3', mid-sporulation element, lowercase letters represent nucleotides with low-frequency occurrence) and the Efg1p (5'-taTGCAta-3') binding motifs [51,54,57] in addition to two high scoring motifs, 5'-TtCtaGaA-3' and 5'-TCGAACCC-3', carrying GAA triplets that are characteristic of HSEs (Figure 8A, shown are motifs found using the global overrepresentation of words relative to control sequences, significance index score (i.e. −log₁₀ *E*-value) >10 for RSAT analyses and >25 for SCOPE analyses). Ndt80p is a transcription factor that controls the expression of genes involved in many cellular processes, including drug resistance, cell separation, morphogenesis and virulence through the recognition of mid-sporulation elements on the promoter of its targets [57,58].

This suggests the presence of functional interactions between Sfl1p, Efg1p and Ndt80p and proposes that Sfl1p binds to two different motifs or that an additional factor binds either 5'-TCGAACCC-3' or 5'-TtCtaGaA-3'. We searched the YeT-FaSCo and the JASPAR databases for similarity with known transcription factor binding sites [59,60]. Interestingly, the 5'-TtCtaGaA-3' sequence was strongly similar to the *S. cerevisiae* Efg1p motif (P = 3.856×10⁻⁰⁴, using YeTFaSco), while database searches did not identify any known motif that closely resembled the 5'-TCGAACCC-3' sequence (data not shown). On the other hand, we found 3 high-scoring motifs in Sfl2p-enriched sequences, including the Efg1p and Ndt80p binding motifs as well as the GAA-containing sequence, 5'-aaNAATAGAA-3' (where N represents any nucleotide; shown are motifs found using the position-analysis program, significance index score >5) (Figure 8B). To confirm that the 5'-aaNAATAGAA-3' motif was specific to Sfl2p, we performed motif discovery analyses using DNA sequences encompassing ±250 bp around peak summits of the regions specifically bound by Sfl2p and found the similar high-scoring motif 5'-aANAATAGAA-3' (Figure 8C). The 5'-aANAATAGAA-3' motif shows moderate similarity with the *S. cerevisiae* Sfl1p and Mga1p motifs (scores = 17.75 and 17.36, respectively using the JASPAR database). All these identified motifs were distributed preferentially around the center of the sequences corresponding to peak locations (Figures 8A, 8B and 8C), suggesting that Sfl1p, Sfl2p, Efg1p and Ndt80p binding sites were very close to each other.

To determine if Efg1p and Ndt80p binding sites overlapped with the genome-wide occupancies of Sfl1p and Sfl2p, we compared Efg1p and Ndt80p binding profiles [51,57] to those of Sfl1p and Sfl2p (Figure 8D). Ndt80p binding was resolved by Sellam *et al.* under yeast-form growth conditions at 30°C [57], whereas Efg1p binding was analysed by Lassak *et al.* during both yeast-form growth (30°C) and hyphal induction (YP serum at 37°C) [51]. Strikingly, a high proportion of Sfl1p and Sfl2p binding sites overlapped with those of Ndt80p (Figure 8D), whereas Efg1p binding overlap was less frequent and depended on the morphological state of *C. albicans*, with rare or no overlap under hyphal induction and increased overlap under yeast-form growth (Figure 8D). Roughly, 90% of Sfl1p and Sfl2p common targets were bound by both Ndt80p and Efg1p (Figure 8D, upper panel as an example), whereas ~10% (10 out of 113 common targets) were bound by Ndt80p but not Efg1p. In at least two cases, Sfl1p and Sfl2p occupancy to common targets overlapped only with Efg1p binding: the promoter regions of *SIS1* and *PDE1*. On the other hand, ~47% of Sfl2p specific targets were bound by both Ndt80p and Efg1p, whereas ~42% overlapped only with Ndt80p binding (Figure 8D, middle panel as an example). On rare occasions (~11%), Sfl2p did not show significant overlap with the binding of any of the three regulators (Figure 8D, bottom panel as an example).

Taken together, our results indicate that Sfl1p and Sfl2p bind to
DNA via divergent motifs and suggest the co-binding of
A

<figure><img src="image_13.png" /></figure>

45 motifs in
36 sequences out of 112

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

47 motifs in
29 sequences out of 112

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

41 motifs in
23 sequences out of 112

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

24 motifs in
20 sequences out of 112

<figure><img src="image_20.png" /></figure>

B

<figure><img src="image_21.png" /></figure>

112 motifs in
81 sequences out of 193

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

98 motifs in
74 sequences out of 193

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

67 motifs in
39 sequences out of 193

<figure><img src="image_26.png" /></figure>

C

<figure><img src="image_27.png" /></figure>

54 motifs in
45 sequences out of 85

<figure><img src="image_28.png" /></figure>

D

<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /></figure>

<figure><img src="image_31.png" /></figure>
**Figure 8.** Sfl1p and Sfl2p binding locations overlap with those of Ndt80p and Efg1p. (A, B and C) Motif discovery analyses of Sfl1p and Sfl2p binding data. Motif logos of conserved sequences in (A) Sfl1p- and (B) Sfl2p-enriched DNA fragments as well as in (C) fragments overlapping with binding regions that are specific to Sfl2p. DNA sequences encompassing ±250 bp around peak summits in Sfl1p or Sfl2p binding data were used as input for motif discovery using two independent motif discovery algorithms, the RSA-tools (RSAT) peak-motifs (http://rsat.ulb.ac.be/rsat/, [55]) and SCOPE (genie.dartmouth.edu/scope/, [56]) (See Materials and Methods for details). High scoring motifs from either SCOPE or RSAT algorithms are shown. These include the Ndt80p and Efg1p binding motifs, suggesting a functional interaction between Sfl1p, Sfl2p, Ndt80p and Efg1p. The distribution of motif occurrences in the input sequences are shown at the right of each motif panel. Plotted are the number of occurrences of each motif (y-axis, motif occurrence) at a given position relative to peak center (distance to peak center in base pairs, x-axis). (D) Overlap of Ndt80p and Efg1p binding with Sfl1p and Sfl1p occupancies at selected locations from the *C. albicans* genome (selected genome interval shown above each panel). Genome-wide location data from Sellam *et al*. (Ndt80p, from 59-bp tiling array data, one of the two replicates of the study is shown [57]) and Lassak *et al*. (Efg1p, from 50-75-mer tiling array data for Efg1p binding in cells grown under yeast form and during hyphal induction [51], one of the three replicates in each condition is shown) are used to compare Ndt80p and Efg1p binding profiles to those of Sfl1p and Sfl2p (read counts in 10 bp windows from wiggle files of Sfl1p and Sfl2p binding data were used).
doi:10.1371/journal.ppat.1003519.g008

transcription factors Efg1p and Ndt80p to many Sfl1p and Sfl2p target promoters, either concomitantly or successively, depending on growth conditions.

## The Efg1p protein binds to the promoter of many Sfl1p and Sfl2p targets and co-immunoprecipitates with Sfl1p and Sfl2p, in vivo

Our bioinformatic analyses suggested the co-binding of Efg1p to many Sfl1p and Sfl2p target promoters. To test whether Sfl1p, Sfl2p and Efg1p concomitantly bind to common targets *in vivo*, strains individually expressing chromosomally TAP-tagged Sfl1p and Sfl2p (strains SFL1-TAP and SFL2-TAP, Table 1) and HA-tagged Efg1p (strain HLCEEFG1, [18], Table 1) under the control of their endogenous promoter were grown in SC medium at 30°C (yeast form-promoting condition) or in Lee's medium at 37°C (filamentous form-promoting condition) during 4 h before being subjected to ChIP-PCR analyses to detect differential binding of the three transcription factors to the promoter of selected Sfl1p and Sfl2p targets (*BRG1*, *EFG1*, *SFL2*, *UME6* and *TEC1*, Figure 9A, see Materials and Methods for details). All strains displayed similar hyphal growth phenotypes at 37°C in Lee's medium, whereas the yeast form growth phenotypes were similar for cells grown in SC medium at 30°C (Figure S1A). Immunoblotting confirmed the expression of the different fusion proteins under the corresponding growth conditions (Figure S1B). As expected, Sfl1p and Efg1p binding was detected at all tested promoters in SC medium at 30°C (Figure 9A, compare lanes 1 and 7 to lanes 2 and 8, respectively). Conversely, in Lee's medium at 37°C, Sfl1p and Efg1p binding was less efficient (Figure 9A, Sfl1p binding, compare lanes 1 and 2 to lanes 4 and 5; Efg1p binding, compare lanes 7 and 8 to lanes 9 and 10). Similarly, Sfl2p binding was detected at all tested promoters in Lee's medium at 37°C (Figure 9A, compare lane 4 to lane 6), whereas in SC medium at 30°C, Sfl2p binding was less efficient (Figure 9A, compare lanes 4 and 6 to lanes 1 and 3).

To further explore the functional interaction between Sfl1p, Sfl2p and Efg1p, we sought to verify if the Efg1p protein could be co-immunoprecipitated with Sfl1p or Sfl2p *in vivo*. To this end, we generated strains co-expressing C-terminally TAP-tagged Sfl1p or Sfl2p and HA-tagged Efg1p (AVL12-SFL1-TAP and AVL12-SFL2-TAP, respectively, Table 1) under the control of their chromosomal promoter together with control strains carrying individual Sfl1p-TAP, Sfl2p-TAP or Efg1p-HA fusions (strains SFL1-TAP, SFL2-TAP and AVL12-pHIS, Table 1, see Materials and Methods). Strains were grown during 4 h in SC medium at 30°C or in Lee's medium at 37°C, followed by crosslinking with formaldehyde to stabilize protein complexes and total extracts were incubated with IgG-coated beads for immunoprecipitation of the Sfl1p-TAP or Sfl2p-TAP proteins in the corresponding strain backgrounds. Immunoblotting with an anti-TAP antibody

(Figure 9B, IP, Anti-TAP panel) allowed to detect the Sfl1p-TAP signal in beads incubated with extracts from strains carrying the *SFL1-TAP* allele irrespective of the growth conditions (i.e. in both SC medium at 30°C and Lee's medium at 37°C) (Figure 9B, IP, Anti-TAP panel, lanes 2, 4, 7 and 9). On the other hand, very low amounts of the Sfl2p-TAP protein fusion were detected in beads incubated with extracts from strains carrying the *SFL2-TAP* allele and grown in SC medium at 30°C (Figure 9B, IP Anti-TAP panel, lanes 3 and 5), however, the Sfl2p-TAP signal strongly increased in Lee's medium at 37°C (Figure 9B, Anti-TAP panel, compare lanes 3 and 5 to lanes 8 and 10). Interestingly, immunoblotting of the bound fractions with an anti-HA antibody (Co-IP, Anti-HA panel) allowed to detect Efg1p-HA co-immunoprecipitation with Sfl1p-TAP under both growth conditions: in SC medium at 30°C and in Lee's medium at 37°C (Figure 9B, CoIP, Anti-HA panel, lanes 2 and 7). Efg1p-HA co-immunoprecipitation with Sfl2p-TAP was barely detectable in SC medium at 30°C but was significantly enhanced in Lee's medium at 37°C, a condition that triggers increased expression of Sfl2p (Figure 9B, CoIP, Anti-HA panel, compare lane 3 to lane 8). As expected, Efg1p-HA was undetectable from beads incubated with strains individually expressing *EFG1-HA*, *SFL1-TAP* or *SFL2-TAP* (Figure 9B, lanes 1, 4, 5, 6, 9 and 10).

Taken together, our results show that i) the Efg1p protein binds to many Sfl1p and Sfl2p targets, *in vivo* and ii) Both Sfl1p and Sfl2p proteins physically associate with Efg1p, *in vivo*.

## Discussion

The ChIP-Seq and transcriptomics technologies are powerful *in vivo* approaches that, when combined, allow to provide mechanistic insights into the function of transcriptional regulators. When associated with both genetic and physical interaction analyses, the overall generated data are cross-validated and provide a comprehensive view of the regulatory interactions within transcriptional networks. They also shed more light into the epistatic relationships to explain the phenotypes associated with transcription factor function. In the present report, we used such approaches to decipher the regulatory network of two HSF-type transcription factors, Sfl1p and Sfl2p, both required for *C. albicans* virulence but with antagonistic functions in regulating *C. albicans* morphogenesis. One limitation of our ChIP-Seq design was the use of ectopic promoter-driven expression of the *SFL1-HA3* and *SFL2-HA3* alleles (Figure 1). This may drive non physiological expression levels and some of the transcriptional changes and promoter occupancies may be altered from the situation where the genes are expressed from their endogenous promoters. Nevertheless, phenotypic analyses suggested that at least *P<sub>MET3</sub>*-driven expression of *SFL2-HA3* imparts filamentous growth in a manner similar to the wild-type SC5314 strain (Figure 1C). Furthermore, we generated strains expressing TAP-tagged *SFL1* and *SFL2* from their
A

<figure><img src="image_32.png" /></figure>

B

<figure><img src="image_33.png" /><figcaption>**Figure 9. Efg1p binds to the promoter of many Sfl1p and Sfl2p targets and co-immunoprecipitates with Sfl1p and Sfl2p, in vivo. (A)** ChIP-PCR assay of selected Sfl1p and Sfl2p target promoters. Strains SFL1-TAP (CEC1922), SFL2-TAP (CEC1918) and EFG1-HA (HLCEEFG1) were grown in SC medium at 30°C (30°C) or in Lee's medium at 37°C (37°C) together with the SC5314 control strain (Control) during 4 h before being subjected to chromatin immunoprecipitation (Anti-TAP, Anti-HA) followed by PCR using primers specific to the indicated promoter regions. The *URA3* and *YAK1* genes were used as negative controls for ChIP enrichment. (B) Co-Immunoprecipitation of Efg1p with Sfl1p and Sfl2p. Strains coexpressing *SFL1-TAP* and *EFG1-HA* (Lanes 2 and 3) or *SFL2-TAP* and *EFG1-HA* (Lanes 7 and 8) or controls (Lanes 1 and 6, *EFG1-HA* only; lanes 4 and 9, *SFL1-TAP* only; lanes 5 and 10, *SFL2-TAP* only) were cultivated in SC medium at 30°C or in Lee's medium at 37°C before crosslinking with formaldehyde. Total extracts were incubated with Dynal PanMouse IgG beads directed against TAP epitope tag prior to washing and Western blotting using anti-TAP (IP Anti-TAP, 10% of the beads/total extracts mixture) and anti-HA (Co-IP Anti-HA) antibodies. A portion of the total cell extracts (~2%) was included to verify the presence of the Efg1p-HA fusion (Total extracts Anti-HA).
doi:10.1371/journal.ppat.1003519.g009</figcaption></figure>
endogenous promoter and ChIP experiments using these strains confirmed some of our data that used the *P<sub>MET3</sub>* expression system (Figure 9A).

Our data allow to propose a model of Sfl1p and Sfl2p transcriptional network (Figure 10, for simplicity only binding associated with transcriptional modulation is shown) as well as a mechanism whereby Sfl1p and Sfl2p antagonistically regulate the yeast-to-hyphae transition (see below). Sfl2p, which responds to temperature increase, and Sfl1p bind to the promoter of common target genes (blue boxes in Figure 10) belonging to at least 3 functional groups involved in morphogenesis: transcriptional repressors of hyphal growth (*SSN6*, *NRG1*, *RFG1*, others), transcriptional activators of hyphal growth (*BRG1*, *UME6*, *TECI*, others) and yeast-form associated genes (*RME1*, *RHD1*, *TWP1*, others). While Sfl1p exerts direct negative and positive regulation on the expression of activators (*BRG1*, *UME6*, *TECI*) and repressors (*SSN6*, *NRG1*) of hyphal growth, respectively, Sfl2p directly upregulates and downregulates the expression of positive (*UME6*, *TECI*) and negative (*RFG1*, *NRG1*) regulators of hyphal growth, respectively (Figure 10). Additionally, Sfl1p directly upregulates the expression of yeast-form associated genes (*RME1*, *RHD1* and *TWP1*) whereas Sfl2p directly downregulates their expression (Figure 10). Moreover, Sfl1p and Sfl2p directly negatively regulate the expression of each other (Figure 10). As stated above, this model is consistent with the genetic interaction analyses performed between *SFL1* (genetically interacts with at least *BRG1* and *SFL2*), *SFL2* (genetically interacts with at least *UME6*, *TECI* and *BRG1*) and their target genes (Figure 7). Importantly, on the other hand Sfl2p exclusively binds to the promoter of specific target genes that belong to at least 2 functional groups involved in morphogenesis: HSGs (*ALS3*, *HGC1*, *HWP1*, *HYR1*, *ECE1*, *SAP4*, *IHD1*, *FAV2*, *RBT4*) and yeast-form specific genes (*PIR1*, *RHD3*) (Figure 10). We propose that binding of Sfl1p and Sfl2p to a high proportion of their transcriptional targets occurs with additional binding of transcription factors Ndt80p and/or Efg1p, depending on growth conditions (Figures 8, 9 and 10), presumably through direct or indirect physical interaction (Figures 8 and 9, see below). One could speculate that the requirement of a functional *EFG1* gene for Sfl1p and Sfl2p abilities to regulate morphogenesis under specific growth conditions (Figure 7 and [39]) could be explained by the need for Efg1p co-binding and/or physical interaction, as suggested by our study (Figures 7, 8 and 9). Indeed, we show here that Efg1p co-immunoprecipitates, *in vivo*, with Sfl1p and Sfl2p and binds to the promoter of many Sfl1p and Sfl2p target genes (Figure 9). On the other hand, our finding that Sfl2p binds exclusively to specific targets, including a high proportion of HSGs (Figure 6), provides additional insight into *SFL2* function. This might explain, for instance, why *SFL2* was able to bypass the need of *EFG1* and *FLO8* to induce hyphal growth in embedded conditions at 37°C [39]. We are currently testing whether Sfl1p and Sfl2p binding to their targets requires the presence of functional *EFG1* or *NDT80* genes. Overall, we propose that the execution of these single (including *SFL1*-*SFL2* cross-factor negative control) and multiple input motifs in Sfl1p or Sfl2p transcriptional network dictates the commitment of the *C. albicans* cells to form hyphae or yeast-form cells. This model is consistent with Sfl1p and Sfl2p acting as “switch on/off” proteins, with Sfl1p directly turning off the expression of positive regulators of hyphal growth while turning on the expression of both yeast-form associated genes and genes encoding repressors of hyphal development, whereas Sfl2p directly turns on the expression of HSGs and positive regulators of hyphal growth while turning off the expression of yeast-form associated genes as well as negative regulators of hyphal development (Figure 10).

The mechanisms whereby HSF-type transcription factors activate transcription involve homotrimerization, post-translation- al modifications (e.g. phosphorylation, others) as well as interaction with multiple protein partners, followed by recruitment of the co- activating mediator complex and initiation of the transcriptional process [61]. This mechanism may include or not nuclear translocation, as many HSFs were shown to reside in the nucleus under both activating and non-activating conditions or to be imported to the nucleus following activation [61]. It was shown that Sfl1p is constitutively localized to the nucleus under both yeast- and hyphae-promoting conditions and irrespective of temperature levels [37,38], whereas an Sfl2p-GFP fusion was undetectable at 25°C but displayed nuclear localization at 37°C [39]. Moreover, *SFL2* RNA levels were undetectable by Northern blotting at either 25°C or 30°C, but were greatly enhanced upon temperature increase [39] and this correlated with Sfl2p protein level variations [39]. Indeed, we show here that in SC medium at 30°C, Sfl2p protein levels are low, but are significantly enhanced upon temperature increase to 37°C in Lee's medium (Figure S1B). Moreover, we show that Sfl2p binding is more stable at 37°C in Lee's medium as compared to 30°C in SC medium, and *vice versa* for Sfl1p (Figure 9A). Based on these observations, we propose the following model of Sfl1p/Sfl2p activation: Sfl1p binds to its transcriptional targets to maintain the yeast form growth at low temperature by directly modulating the expression of genes involved in morphogenesis (Figure 10). A temperature increase to 37°C leads to an increase in both Sfl2p expression and binding to the promoter of Sfl1p targets in addition to specific targets (including HSGs) and induction of the hyphal development program (Figure 10). As we show here that Sfl1p and Sfl2p act as both activators and repressors of gene expression (Figures 6 and 10), it is likely that they alternatively recruit (directly or indirectly) co-repressors (e.g. Tup1p-Ssn6p) and co-activators (e.g. mediator- Swi/Snf complex) at different binding sites to regulate morpho- genesis. Our observation that Sfl2p binds to its own promoter, but not Sfl1p (Figures 3, 6A and 10) is consistent with this model as *SFL2* may undergo auto-induction which would lead to a rapid, amplified and sustained expression of *SFL2*, allowing an efficient response to temperature increase. On the other hand, *SFL1* expression, protein levels and nuclear localization remain constant under various conditions [38], which may dispense the need for autoregulation. The *SFL1*-*SFL2* cross-factor negative control is also consistent with this model. Under low temperature conditions, Sfl1p directly turns off *SFL2* expression to prevent activation of hyphal growth. Upon a temperature increase, *SFL2* expression is enhanced and Sfl2p binds to the *SFL1* promoter to turn off *SFL1* expression. This allows to relieve Sfl1p-mediated repression, thus contributing to activation of the hyphal development program.

Our motif discovery analyses suggested that Ndt80p co-binds together with Efg1p to the promoter of Sfl1p and Sfl2p targets (Figure 8). We also strikingly found that a high proportion of Sfl1p and Sfl2p binding sites overlapped with those of Ndt80p and/or Efg1p (Figure 8). However, since the Ndt80p ChIP-on-chip was performed on yeast-form grown cells at 30°C [57], one cannot exclude the possibility that Ndt80p binding is altered/lost upon hyphal induction, as is obviously the case for Efg1p ([51] and Figures 8D and 9A). Ndt80p occupies the promoter region of roughly a quarter of total *C. albicans* genes under yeast-form growth conditions, suggesting wide functions for Ndt80p [57]. Indeed, it was shown that Ndt80p regulates different processes including drug resistance, cell separation, hyphal differentiation, biofilm formation and virulence [54,57,58]. Importantly, the *C. albicans* *ndt80Δ/ndt80Δ* mutant is unable to form true hyphae under different filamentation-inducing conditions and, in the
<figure><img src="image_34.png" /><figcaption>**Figure 10. Model of Sfl1p and Sfl2p regulatory network.** Sfl2p (red oval), which induces hyphal growth in response to temperature increase or upon overexpression (red dashed arrow), and Sfl1p (orange oval) bind directly, together with Efg1p and Ndt80p depending on growth conditions (green and white ovals, respectively; dashed lines indicate hypothetical physical and/or functional interaction), to the promoter of common (blue boxes) target genes encoding major transcriptional activators (*UME6*, *TEC1* and *BRG1*) or repressors (*NRG1*, *RFG1*, *SSN6*) of hyphal growth as well as to the promoter of genes associated with yeast-form growth (*RME1*, *RHD1* and *YWP1*) and modulate the expression of many of them (for simplicity, only modulatory direct interactions are shown i.e. both binding at and transcriptional modulation of a given target; arrowed lines indicate direct upregulation whereas blunt lines indicate direct downregulation). On the other hand, Sfl2p directly upregulates the expression of specific targets (grey boxes), including a high proportion of hyphal-specific genes (HSGs), while exerting a direct negative regulation on the expression of yeast-form associated genes (*PIR1* and *RHD3*). Sfl1p and Sfl2p also exert a direct negative regulation on the expression of each other. The execution of Sfl1p or Sfl2p transcriptional control inputs allows to regulate the commitment (dashed line; blunt, inhibition; arrowed, activation) of *C. albicans* to form hyphae or yeast-form cells.</figcaption></figure>

doi:10.1371/journal.ppat.1003519.g010

presence of serum at 37°C, it fails to activate the expression of HSGs, including *HWP1*, *ECE1*, *RBT4*, *ALS3*, *HYR1* and *SAP4* [58], all directly regulated by Sfl2p (Figure 6), as well as the transcription factor-encoding genes *TEC1* and *UME6* which are both directly modulated by Sfl1p and Sfl2p (Figure 6). Additionally, under the same growth conditions, the homozygous *ndt80* mutant was unable to downregulate the yeast form-associated genes *YWP1*, *RHD3*, *RHD1* and the transcriptional repressor-encoding gene *NRG1* [58], which are also direct targets of Sfl1p or Sfl2p (Figure 6). These observations, together with our findings that i) Ndt80p binding motif was enriched among Sfl1p and Sfl2p bound sequences and that ii) a significant proportion of its genome-wide binding profile overlapped with Sfl1p and Sfl2p binding, suggest that Sfl1p, Sfl2p and Ndt80p cooperatively regulate *C. albicans* morphogenesis in response to temperature variation. Whether Sfl1p and Sfl2p regulate this process through physical interaction with Ndt80p and the associated sequence of molecular events occurring during the yeast-to-hyphal switch await further characterization. On the other hand, we found that Efg1p binding also overlapped with that of Sfl1p and Sfl2p, at a lesser extent, though, as compared to Ndt80p binding (Figure 8). It is intriguing that Efg1p binding undergoes alteration following the induction of hyphal development ([51] and Figures 8D and 9A). Our examination of Efg1p binding data by Lassak *et al.* [51] together with our ChIP experiments (Figure 9A) suggest that Efg1p binding to many targets is decreased/altered upon hyphal induction. We show here that during yeast-form growth, at low

temperature, Efg1p co-immunoprecipitates with Sfl1p but not with Sfl2p, presumably due to the low levels of Sfl2p at low temperature (Figure 9B). One could speculate that, at low temperature, Sfl1p associates directly or indirectly with Efg1p on the promoter of its targets to repress hyphal development. Following a temperature increase, both Sfl2p levels and Sfl2p DNA binding are enhanced (Figures S1 and 9A), which in turn activates the hyphal development program. Although Efg1p binding is altered upon hyphal induction, Efg1p co-immunoprecipitated with Sfl2p (Figure 9B) at 37°C in Lee's medium, which may explain Sfl2p dependency on *EFG1* to regulate morphogenesis under certain conditions. Nobile *et al.* elegantly showed that an intricate transcriptional network involving Ndt80p, Efg1p, Brg1p, Bcr1p, Rob1p and Tec1p controls biofilm development in *C. albicans* [54]. Interestingly, with the exception of *BCR1*, all genes encoding these regulators are direct targets of Sfl1p or Sfl2p (Figure 6 and [54]). It is tempting to speculate that Sfl1p and Sfl2p may convey temperature regulation to the transcriptional network controlling biofilm formation.

*C. albicans* adaptation to temperature variation is one of the major critical traits of its ability to cause disease or to act as a commensal of warm-blooded species, as a temperature increase triggers hyphal development [2]. To date, three temperature-responsive transcription factors have been shown to play a role in *C. albicans* morphogenesis, Hsflp [62,63], Sfl2p [39,40] and Hms1p [49]. Importantly, all three transcription factors are required for full virulence in different host/tissue models [39,40,49,63], reinforcing
the link between temperature adaptation and pathogenesis in *C. albicans*. The *HMS1* gene, encoding a basic helix-loop-helix (bHLH) transcription factor, has been recently isolated in a screen aimed at identifying transcription factors whose function is required for the *HSP90*- or high temperature-mediated filamentous growth [49]. Hms1p acts downstream of the Pho85p-Pcl1p cyclin-dependent kinase pathway but its function was still dependent upon cAMP-PKA signalling [49]. Interestingly, both Sfl1p and Sfl2p bind to the promoter of the *HMS1* gene, while Sfl2p downregulates its expression (Figure 6A), suggesting that activation of Sfl2p turns off the *HSP90*-dependent filamentation response (at least under the conditions used in the present study). Similar to Sfl2p, Hsf1p is an HSF-type transcription factor that induces transcription following a temperature increase, but, unlike *SFL1* and *SFL2*, *HSF1* is essential for viability [62]. Hsf1p is required for the expression of essential chaperones, including *HSP104*, *HSP90*, *HSP70* as well as other classical heat-shock protein (HSP)-encoding genes such as *HSP60*, *HSP78*, others [62]. Although carrying HSF-type domains in their primary protein sequences and sharing relatively high sequence similarity levels with Hsf1p, speculating a role in the transcriptional regulation of HSP (or HSP-related) genes, the Sfl1p and Sfl2p binding targets did not show any significant enrichment of functional categories pertaining to the heat-shock response pathway (e.g. protein folding/refolding), including HSPs and chaperones (Figure 2C). This may have important evolutionary implications as it might reflect specific needs of *C. albicans* to efficiently act as an opportunistic yeast of warm-blooded animals through converting temperature-sensing inputs into a morphogenesis programming output using HSF-type regulators like Sfl1p and Sfl2p. Nevertheless, we detected Sfl1p and Sfl2p binding at the promoter of the *HSP104*, *HSP70* and *SIS1* genes (binding intensity below algorithm threshold used for *HSP70*), suggesting that a reminiscent classical heat-shock response may have been retained in Sfl1p and Sfl2p. It is intriguing that one of the two potential binding motifs of Sfl1p (Figure 8A), 5'-TtCtaGaA-3', is strikingly similar to the *S. cerevisiae* Hsf1p motif [64,65], in line with the hypothesis that transcriptional rewiring affected the regulation of the heat shock response and temperature adaptation between *S. cerevisiae* and *C. albicans*.

It is worth noting that the predicted protein sequences of Sfl1p and Sfl2p are highly similar to those of *S. cerevisiae* Sfl1p and Mga1p. The *MGA1* gene has been initially isolated as a multicopy suppressor of both the *snf2Δ* (component of the SWI/SNF remodelling complex, also known as *gam1*) [66] and the *mep1Δ/mep1Δ mep2Δ/mep2Δ* (encoding ammonium permeases) filamentous defect [67] mutations in *S. cerevisiae*. Interestingly, Mga1p was shown to act as a master regulator of *S. cerevisiae* pseudohyphal development through direct transcriptional control of key genes involved in morphogenesis [68]. Many intriguing functional similarities exist between Sfl2p and *S. cerevisiae* Mga1p, although either *SFL1* or *SFL2* could complement an *sfl1Δ* mutation and *SFL2* could not complement the pseudohyphal growth defect of an *mga1Δ* mutant [39]. First, both proteins recognize similar DNA binding motifs (5'-AtAGAAC-A-3' for Mga1p [33] and 5'-ANATAGAA-3' for Sfl2p (Figure 8)). Second, both transcription factors bind to the promoter of orthologous genes (*ScPHD1* and *ScSOK2/CaEFG1*, *HMS1*, *ScGAT2/CaBRG1*, *MSB2*, *ACH1*, *ScENA1/CaENA21*, *GCN4*, *CUP9*, *TPO4*, *ScSCW4/CaMP65*, others; binding to some genes is below peak-finding algorithm threshold). Third, the regulatory networks to which they belong are intriguingly similar: Mga1p establishes cross talks with major regulators of *S. cerevisiae* pseudohyphal growth including Phd1p, Sok2p (Efg1p orthologs), Flo8p and Tec1p, as in the case of Sfl2p (Figure 6) [39,68]. Fourth, overexpression of *MGA1* and *SFL2* is sufficient to induce morphogenesis in the respective species under

conditions that do not promote filamentation [39,68]. Fifth, Sfl2p requires *EFG1* and *FLO8* to induce filamentation under specific conditions (Figure 7B and [39]) and we show here that Efg1p co-immunoprecipitates with Sfl2p (Figure 9B). Similarly, Mga1p requires a functional *FLO8* gene for its ability to bind DNA and Mga1p and Flo8p interact with each other [68]. We suggest that transcriptional rewiring may have affected the functions of Sfl2p and Mga1p in their respective species: In diploid *S. cerevisiae* cells, Mga1p responds to nitrogen limitation to turn on pseudohyphal growth, whereas in *C. albicans* Sfl2p responds to temperature increase to induce hyphal development.

## Materials and Methods

### Strains and growth media

The *C. albicans* strains used in this study are listed in Table 1. Depending on experimental conditions, *C. albicans* strains were grown in YPD (1% yeast extract, 2% peptone, and 1% dextrose), YP (1% yeast extract, 2% peptone) supplemented with 10% Fetal Bovine Serum (FBS), SD (synthetic dextrose, 0.67% yeast nitrogen base (YNB; Difco) with 2% glucose) [69] supplemented if necessary with arginine, histidine or uridine (20 mg/l each and 2% agar for growth on solid medium), SC (synthetic complete) or Lee's medium supplemented or not with methionine [70]. Expression from the tetracycline-inducible promoter (*P*<sub>TET</sub>) was achieved through addition of 3 µg/ml anhydrotetracycline (ATC - Fisher Bioblock Scientific) in YPD at 30°C [41]. ATC-containing cultures were maintained in the dark as ATC is light sensitive. *Escherichia coli* strains TOP10 (Invitrogen) or DH5α were used for DNA cloning and maintenance of the plasmid constructs.

### Plasmid construction and generation of epitope-tagged or mutant strains

All *C. albicans* transformation experiments used the lithium-acetate transformation protocol of Walther and Wendland [71] and selection of transformants for uridine or histidine prototrophy (when using the *URA3* or the *HIS1* markers, respectively) or Nourseothricine resistance (when using the *SAT1* marker) [72]. Plasmid pCaMPY-3×HA and the SGY243 strains expressing the *CAP1-HA3* allele or carrying the empty vector (pCaEXP) were kindly provided by Dr Martine Raymond (Université de Montréal, Canada). Strains AVL12 and HLCEEFG1 (expressing EFG1-HA under the control of the endogenous promoter) were the kind gifts of Dr Joachim Ernst (Heinrich-Heine-Universität, Dusseldorf, Germany). We first attempted to generate epitope (HA3, triple hemagglutinin)-tagged strains expressing Sfl1-HA3 or Sfl2-HA3 under the control of their endogenous promoter at their chromosomal location. *SFL1*- or *SFL2*-tagging cassettes were PCR-amplified from plasmid pCaMPY-3×HA [73] using primers SFL1-HA-FWD (forward, Table S9 in Text S1, the lowercase sequence corresponds to positions +2316 to +2415 of the *SFL1* ORF) and SFL1-HA-REV (reverse, Table S9 in Text S1, the lowercase sequence corresponds to positions +2419 to +2518 of the *SFL1* ORF) or primers SFL2-HA-FWD (forward, Table S9 in Text S1, the lowercase sequence corresponds to positions +2043 to +2142 of the *SFL2* ORF) and SFL2-HA-REV (reverse, Table S9 in Text S1, the lowercase sequence corresponds to positions +2146 to +2245 of the *SFL2* ORF), which anneal specifically to the in-frame pCaMPY-3×HA vector sequences PET-up and PET-down (respective uppercase sequences in Table S9 in Text S1), as described previously [73]. The resulting fragments (1,853 bp), containing the *C. albicans* *URA3* marker flanked by direct repeats of the HA3-encoding sequences and 100 bp of sequences homologous to the 3' end of the *SFL1* or *SFL2* genes, were used to
respectively transform *ura3*-deficient *sfl1Δ/SFL1* and *sfl2Δ/SFL2* heterozygous mutants, yielding strains CEC3075 and CEC3076, respectively (Table 1). Expression of the Sfl1p-HA₃ and Sfl2p-HA₃ fusions in strains CEC3075 and CEC3076 was not detectable by Western blot analyses, suggesting that integration of the tagging cassette at the 3' untranslated regions of *SFL1* and *SFL2* had a knockdown effect. Despite many attempts, excision of the *URA3* marker through intramolecular recombination between the HA₃ sequences was not successful. We rather observed 100% loss of the entire tagging cassette at the *SFL1* and *SFL2* loci. We therefore used the pCaEXP system to drive expression of the tagged *SFL1* and *SFL2* alleles at the *RPS1* locus [42]. The *SFL1-HA3* or *SFL2-HA3* fusions were PCR amplified from CEC3075 or CEC3076 genomic DNA, respectively, using primers SFL1-HA-CaEXP-FWD (forward, Table S9 in Text S1, introduces a *Bg*I site [underlined]) or SFL2-HA-CaEXP-FWD (forward, Table S9 in Text S1, introduces a *Bg*I site [underlined]), respectively, and primer HA-CaEXP-REV (reverse, Table S9 in Text S1, introduces sequentially a *Bg*I site [underlined] and a TAA stop codon [in red lowercase letters]). The resulting fragments (*SFL1-HA3*, ~2,600 bp; *SFL2-HA3*, ~2,330 bp) were digested with *Bg*I and cloned into the compatible *Bam*HI site of plasmid pCaEXP, generating plasmids pCaEXP-*SFL1-HA3* and pCaEXP-*SFL2-HA3*. Plasmids pCaEXP (empty vector, control), pCaEXP-*SFL1-HA3* and pCaEXP-*SFL2-HA3* were digested with *Stu*I for integration at the *RSP1* locus [42] and the resulting fragments were used to transform strains CEC1910 and CEC1503 (Table 1), respectively, to generate strains *sfl1-CaEXP*, *sfl1-CaEXP-SFL1-HA3*, *sfl2-CaEXP* and *sfl2-CaEXP-SFL2-HA3* (Table 1).

Construction of *C. albicans* knock-out mutants (Table 1) used PCR-generated *ARG4*, *HIS1*, *URA3* and *SAT1* disruption cassettes flanked by 100 base pairs of target homology region (primer sequences are listed in Table S9 in Text S1) as described by Gola *et al.* [74] and Schaub *et al.* [75]. Independent transformants were produced and the gene replacements were verified by PCR on whole yeast cells as described previously [74,75]. If necessary, transformants were converted to uracil prototrophy using *Stu*I-linearized *CIP*10 [76]. Mutant strains carrying the pCIP-*P*<sub>TET</sub>-*SFL2* [41] plasmid (Table 1) were first transformed with the pNIMX construct as described in Chauvel *et al.* [41].

Construction of chromosomally TAP-tagged *SFL1* and *SFL2* alleles (Table 1) used PCR-generated tagging cassettes from plasmid pFA-TAP-HIS, a derivative of the pFA-GFP-tagging plasmid series [74] (primers are listed in Table S9 in Text S1, oligos # 50-53) followed by targeted homologous recombination at the 3' untranslated regions of *SFL1* and *SFL2* to generate strains expressing C-terminally tagged Sfl1p (strains SFL1-TAP and AVL12-SFL1-TAP, Table 1) and Sfl2p (strains SFL2-TAP and AVL12-SFL2-TAP, Table 1) proteins.

## Total protein preparation and Western blotting

Total protein extracts were prepared from 24 OD₆₀₀ units of strains expressing (*sfl1*-CaEXP-SFL1-HA, *sfl2*-CaEXP-SFL2-HA) or not (empty vector; *sfl1*-CaEXP, *sfl2*-CaEXP) the *SFL1-HA3* or *SFL2-HA3* alleles (Table 1) grown overnight in SD medium (*P*<sub>MET3</sub>-inducing conditions). Cultured cells were centrifuged at 3,500 rpm during 5 min at room temperature and the pellets were resuspended in 150 µl of ice-cold TE buffer (10 mM Tris, [pH 7.5], 1.5 mM EDTA) supplemented with a protease inhibitor cocktail (Roche) and 1.5 mM phenylmethylsulfonyl fluoride (PMSF) then transferred to 1.5-ml tubes. The equivalent of 100 µl ice-cold glass beads was added to each tube and the suspensions were vortexed 5 times during 1 minute with 1-min incubations on ice in between. The extracts were clarified

by centrifugation at 5,000 rpm during 1 min, boiled for 1 min and separated (25 µl) by electrophoresis on a sodium dodecyl sulfate-8% polyacrylamide gel. Proteins were electrophoretically transferred to nitrocellulose membranes. The membranes were incubated with a mouse anti-HA monoclonal antibody (12CA5; Roche) for 1 h at a dilution of 1:1,000, followed by incubation with a horseradish peroxidase-conjugated secondary antibody (Sigma) during 30 min, washed, and developed with enhanced chemiluminescent detection reagents (ECL kit, GE Healthcare).

## Microscopy and image analyses

Cells were observed with a Leica DM RXA microscope (Leica Microsystems). Images were captured with a Hamamatsu ORCA II-ER cooled CCD camera, using the Openlab software version 3.5.1 (Improvision Inc.).

## ChIP-Seq, data preprocessing and analyses

Two independent cultures of strains *sfl1*-CaEXP or *sfl2*-CaEXP (untagged; control strains) and *sfl1*-CaEXP-*SFL1-HA3* or *sfl2*-CaEXP-*SFL2-HA3* (tagged strains) (Table 1) were grown overnight in 2 ml YPD at 30°C, diluted to an OD₆₀₀ of 0.3 in Lee's medium deprived of methionine and cysteine (to induce *P*<sub>MET3</sub>) and grown during 4 hours at 37°C (hyphae-inducing conditions). The subsequent steps of DNA cross-linking, DNA shearing and chromatin immunoprecipitation (ChIP) were conducted as described in Liu *et al.* [73], with some modifications. Briefly, cultures were treated with 1% formaldehyde (cross-linking) and snap-frozen in liquid nitrogen. Total cell extracts were prepared by bead beating using a FastPrep-24 instrument (MP Biomedicals) with 6 runs during 1 minute each at 6.0 m/sec and 1 minute on ice in between (these settings led to efficient breakage of hyphal cells). Preparation of soluble chromatin fragments was performed by sonicating the extracts 6 times during 20 sec at power 8 (knob position) for an output signal amplitude of 15 (Microns, Peak to Peak) using a probe sonicator (MSE), yielding ~200-bp DNA fragments on average. The extracts were then incubated at 4°C overnight with a mouse monoclonal anti-HA antibody (Santa Cruz Biotech) coupled to magnetic beads (pan-mouse immunoglobulin G Dynabeads; Dynal Biotech, Brown Deer, WI). The concentration of the purified immunoprecipitated DNA was ranging between 0.2 ng/µl and 1.5 ng/µl in 50 µl TE (10 mM Tris [pH 8.0], 1 mM EDTA). Library construction (10 ng of the immunoprecipitated DNA were used, adaptor-DNA fragments ranging from 150 to 350 bp) was performed using the TruSeq DNA sample preparation kit as recommended by the manufacturer (Illumina), followed by quality control analyses using a Bioanalyzer 2100 instrument (Agilent Technologies). DNA library samples were indexed and pools of the Sfl1p (4 samples, both tagged and control) or Sfl2p (4 samples, both tagged and control) ChIP samples were loaded onto two lanes of an Illumina HiSeq2000 sequencer flow cell for single-read (51 base pairs per read) high-throughput sequencing. The resulting 51-nucleotide sequence reads (FASTQ files) were imported into the Galaxy NGS data analysis software (https://main.g2 bx.psu.edu/) and the tools implemented in Galaxy were used for further processing via workflows [77,78]. Quality control analyses of the FASTQ files were performed using FastQC (version 0.10.0, Babraham Institute) and adaptor-contaminated sequences were trimmed. The reads were then mapped to the *C. albicans* assembly 21 genome using the Bowtie algorithm [79] and the files of mapped reads (BAM files) for the ChIP sample (2 biological replicates from samples *sfl1*-CaEXP-*SFL1-HA3* or *sfl2*-CaEXP-*SFL2-HA3*) and from the control (2 biological replicates from samples *sfl1*-CaEXP
or *sfl2-CaEXP* were processed using the command line version 1.4Orc2 of the Model-Based Analysis for ChIP-Seq (MACS) peak-finding algorithm [46] for peak finding with the following parameters: bandwidth = 250; mfold = 10,30; shiftsize = 100; P-value cutoff for Sfl1p peaks = 1e-14 and P-value cutoff for Sfl2p peaks = 1e-100. Replicates 1 and 2 from the two independently performed ChIP-Seq experiments were processed separately. Overlapping peak intervals (intersection) from replicates 1 and 2 of Sfl1p or Sfl2p binding data were generated using the Galaxy tool Intercept version 1.0.0 (https://main.g2.bx.psu.edu/). The complete Sfl1p and Sfl2p binding and expression datasets are provided in Tables S1–S8 in Text S1. The command line version of the PeakAnnotator (v 1.4) sub-package from the PeakAnalyzer suite of algorithms [80] was used to annotate the Sfl1p and Sfl2p binding peaks in Tables S1, S2, S4 and S5 in Text S1. The association of peaks to target genes was also conducted by human eye (Tables S3 and S6 in Text S1), based on the location of ORFs relative to binding peaks. We provide wiggle tracks with tag counts for every 10 bp segment (See Materials and Methods section entitled "Data accession numbers" below). Visualization of the ChIP-Seq results was conducted using the Integrated Genomics Viewer software [44,45].

## ChIP-PCR assays

Thirty cycles of PCR with 15 seconds at 95°C, 15 seconds at 50°C and 40 seconds at 70°C were performed on independently generated ChIP samples (Figures 3 and 9A) in a 50-µl reaction volume with 1 µl (5%) of immunoprecipitated material. Primers were designed to assay binding enrichment approximately around ChIP-Seq peak summits (primer sequences are provided in Table S9 in Text S1). The *URA3* and *YAK1* ORFs were used as negative controls.

## RNA isolation for microarray experiments

Strains *sfl1-CaEXP* or *sfl2-CaEXP* (control strains, for subsequent Cy3 labeling) and *sfl1-CaEXP-SFL1-HA3* or *sfl2-CaEXP-SFL2-HA3* (test strain, for subsequent Cy5-labeling) (Table 1) were grown overnight in 2 ml YPD at 30°C. The next day, an aliquot of the overnight culture was used to inoculate 50 ml of Lee's medium deprived of methionine and cysteine to a starting OD₆₀₀ of 0.3. This culture was grown for 4 hours at 37°C, cells were washed with diethyl pyrocarbonate (DEPC)-treated water, collected by centrifugation and pellets were immediately frozen and stored at -80°C until RNA isolation. Three independently obtained sets of cell cultures were used. RNA was isolated from frozen cell pellets using the hot-phenol method [81]. Briefly, cells were resuspended in 375 µl TES buffer (10 mM Tris [pH 7.5], 10 mM EDTA, 0.5% SDS) at room temperature, after which 375 µl acid Phenol:Chloroform (5:1, Amresco, Solon, OH) were added. Samples were then incubated for 1 hour at 65°C with vigorous vortexing during 20 sec each 10 min and subjected to centrifugation for 20 min at 14,000 rpm. The supernatants were transferred to new tubes containing 750 µl acid Phenol:Chloroform (5:1), mixed, and subjected to centrifugation at 14,000 rpm for 10 min. The aqueous phase was transferred to new tubes containing 750 µl Chloroform:Isoamyl alcohol (24:1, Interchim, Montluçon, France), mixed and centrifuged at 14,000 rpm during 10 min. RNA was precipitated from the resulting aqueous layer by mixing that portion in new tubes with 1 ml 99% ethanol (pre-cooled at -20°C) and 37 µl of 3 M sodium acetate [pH 5.0] and subjecting the mixture to centrifugation at 14,000 rpm for 40 min at 4°C. The supernatants were removed, the pellet was resuspended in 500 µl 70% ethanol, and the RNA was collected by centrifugation at 14,000 rpm for 20 min at 4°C. The

supernatants were again removed, and the RNA was resuspended in 150 to 300 µl DEPC-treated water. The RNA was stored at -80°C until needed.

## First-strand cDNA synthesis and microarray hybridization

Prior to first-strand cDNA synthesis, the purity and concentration of RNA samples were determined from A260/A280 readings (NanoVue Plus, GE Healthcare) and RNA integrity was determined by a Bioanalyzer 2100 instrument (Agilent Technologies) per the manufacturer's instructions (RNA concentration was ranging between 7.92 and 10.48 µg/µl). First-strand cDNA was synthesized from 20 µg total RNA, using the Superscript III indirect cDNA labeling system (Invitrogen) with the following minor modifications to the manufacturers' instructions. Briefly, the Qiagen PCR Purification kit was used to remove unincorporated aminoallyl-dUTP and free amines with substitution of the Qiagen-supplied buffers with phosphate wash (5 mM Phosphate buffer [K₂HPO₄/KH₂PO₄O₄] [pH 8.0], 80% ethanol) and elution (4 mM Phosphate buffer [K₂HPO₄/KH₂PO₄O₄] [pH 8.5]) buffers. The purified first-strand cDNAs were subsequently labelled with the monoreactive Cy dye N-hydroxysuccinimide esters Cy3 (control, cDNA from strains *sfl1-CaEXP* or *sfl2-CaEXP*) and Cy5 (cDNA from strains *sfl1-CaEXP-SFL1-HA3* or *sfl2-CaEXP-SFL2-HA3*) (GE Healthcare) and the uncoupled dye was removed using the standard Qiagen PCR purification kit protocol. The Cy3- and Cy5-labeled cDNA lyophilized pellets were resuspended in 10 µl of DNase-free water then 2.5 µl and 12.5 µl of 10X blocking agent and 2X hybridization buffer (Agilent Technologies), respectively, were added. The resulting samples were mixed, incubated at 95°C during 3 min and snap cooled on ice during 1 min then hybridized to a *Candida albicans* expression array (Agilent Technologies) designed such that two nonoverlapping probe sets are targeting each of 6,105 *C. albicans* ORFs for a total of 15,744 probes, thereby allowing two independent measurements of the mRNA level for a given gene (The EMBL-European Bioinformatics Institute ArrayExpress platform accession number: A-MEXP-2142, http://www.ebi.ac.uk/arrayexpress/arrays/A-MEXP-2142).

## Gene expression microarray data analysis

Images of Cy5 and Cy3 fluorescence intensities were generated by scanning the expression arrays using an Axon Autoloader 4200AL scanner (Molecular Devices, Downington, PA). Images were subsequently analyzed with the GenePix Pro 6.1.0.2 software (Molecular Devices, Downington, PA). GenePix Results (GPR) files were imported into the Arraypipe 2.0 [82] or the GeneSpring (Agilent Technologies) softwares. Following spot filtering and bad spot flagging, global signal intensities were normalized using Loess normalization and replicate slides (n = 3) were combined and the P-values calculated using a standard Student's t-test.

## Quantitative RT-PCR analyses

Total RNA was prepared from strains CEC2001 (*sfl1Δ/sfl1Δ*) and CEC1997 (*sfl1Δ/sfl1Δ P*<sub>*PCKT*</sub>-SFL1-TAP) or CEC1535 (*sfl2Δ/sfl2Δ*) and CEC1509 (*sfl2Δ/sfl2Δ P*<sub>*PCKT*</sub>-SFL2-TAP) (Table 1) during a kinetics experiment (0 h, 2 h and 4 h) in YNB plus 2% casaminoacids (*P*<sub>*PCKT*</sub>-inducing conditions). Cells from 100 mL cultures were mechanically disrupted with glass beads using a Fastprep (MP Biomedicals) and total RNA was extracted using RNAeasy (QIAGEN) according to the manufacturer's instructions. The quality and quantity of the isolated RNA were determined using an Agilent 2100 Bioanalyzer. Before cDNA synthesis, total RNA samples were DNase-treated using the Turbo DNA-free kit (Ambion). 2 µg of total RNA were used to perform cDNA
synthesis using Superscript II Reverse Transcriptase according to the manufacturer's instructions (Invitrogen). Quantitative PCR was carried out on a Mastercycler ep realplex (Eppendorf) with a 2X SYBR Green master mix (SYBR Green Power, Applied Biosystems). The oligonucleotide primers used are listed in Table S9 in Text S1 (oligos # 18–27). The reaction mixture contained 2.5 µM of each primer and 5 µL of cDNA at 1:10, 1:100 or 1:1000 dilutions. Each sample was processed in triplicate. Relative expression levels were calculated using the delta-delta Ct (ΔΔCt) method, with *C. albicans* translation elongation factor *CEF3* transcript as a calibrator. The relative expression was calculated as 2^(*Ct* target - *Ct* *CEF3*) / 2^(*Ct* target - *Ct* *CEF3*) where *Ct* target is the Ct value for the target gene and *Ct* *CEF3* is the Ct value for the *CEF3* gene. The relative expression was calculated as 2^(*Ct* target - *Ct* *CEF3*) / 2^(*Ct* target - *Ct* *CEF3*) where *Ct* target is the Ct value for the target gene and *Ct* *CEF3* is the Ct value for the *CEF3* gene.

## Co-immunoprecipitation experiments

Strains co-expressing Sfl1p-TAP and Efg1p-HA or Sfl2p-TAP and Efg1p-HA (AVL12-SFL1-TAP or AVL12-SFL2-TAP, respectively, Table 1) together with the control strains SFL1-TAP, SFL2-TAP and AVL12-pHIS (Table 1) were grown during 4 h in 50 ml SC medium at 30°C or Lee's medium at 37°C prior to crosslinking with formaldehyde. Cells were lysed with glass beads and total extracts were prepared in 700 µl lysis buffer (50 mM HEPES-KOH pH 7.5, 140 mM NaCl, 1 mM EDTA, 1% Triton X-100, 0.1% Na-deoxycholate) then sonicated as described for the ChIP-Seq experiment. Immunoprecipitation was performed with 500 µl of clarified sonicated extracts and 40 µl of IgG-coated magnetic beads (Dynabeads Pan mouse IgG, Invitrogen), previously preybed overnight with PBS-0.1% BSA. The beads were washed once with 1 ml lysis buffer and three times with lysis buffer supplemented with 150 mM NaCl. Reverse crosslinking was achieved by incubating beads at 100°C during 25 min in reverse-crosslinking buffer (2% SDS, 0.5 M 2-mercaptoethanol, 250 mM Tris, pH 8.8). The immunoprecipitates were resolved by electrophoresis on an 8% SDS-polyacrylamide gel. Proteins were electrophoretically transferred to nitrocellulose membranes. Blots were revealed with rat monoclonal anti-HA peroxidase conjugate - High Affinity (clone 3F10, Roche) for detection of co-immunoprecipitated Efg1p-HA or with Peroxydase-Anti-Peroxydase Soluble complex (Sigma Aldrich) for detection of immunoprecipitated Sfl1p-TAP and Sfl2p-TAP at a 1:2000 dilution.

## Bioinformatic analyses

Gene Ontology functional enrichment analyses were conducted using the CGD Gene Ontology (GO) Term Finder tool (http://www.candidagenome.org/cgi-bin/GO/goTermFinder). The orf19 list of the Sfl1p and Sfl2p common targets or the orf19 list of the Sfl2p-specific targets was used as input for functional grouping. To decide which of the two ORFs sharing the same bound promoter are included among the GO-term finder input list, we selected those ORFs showing differential expression in Sfl1p and Sfl2p transcriptomics data (expression level fold-change ≥1.5, P-value ≤0.05). This led to a list of 110 (Sfl1p and Sfl2p common targets) and 73 (Sfl2p specific targets) genes for GO term enrichment analyses (Table 2). If some GO terms contained overlapping gene lists, the GO term with the largest number of genes or with the best significance score was selected. The P-value cutoff for considering a functional grouping enrichment was P≤0.05. For motif discovery analyses, peak summit location files generated by the MACS algorithm [46] were imported into the Galaxy NGS analysis pipeline and DNA sequences encompassing ±250 bp around peak summits in Sfl1p or Sfl2p data sets were extracted using the Extract Genomic DNA tool version 2.2.2. The resulting sequences were used as input for motif discovery using

the SCOPE (Suite for Computational Identification of Promoter Elements, version 2.1.0) program (http://genie.dartmouth.edu/scope/) [56] or the Regulatory Sequence Analysis Tools ([RSAT] http://rsat.ulb.ac.be/rsat/) peak-motifs algorithm [55]. The parameters used in RSAT peak-motifs algorithm were as follows: oligo-analysis and position-analysis were selected; oligo length was 6 and 7; the Markov order (m) of the background model for oligo-analysis was set to automatically adapt to sequence length; the number of motifs per algorithm was 10 and both strands of the DNA sequence inputs were searched for motif discovery. For building a control set of sequences (that is sequences randomly chosen from the genome), we used the RSA tool "random genome fragments". The parameters used in SCOPE were as follows: species selected was *C. albicans* (genome sequence available at www.broad.mit.edu/annotation/genome/); "fixed" was selected for the upstream sequence control set and both strands of the DNA sequence inputs were searched for motif discovery.

## Data accession numbers

ChIP-Seq and microarray data can be found at the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/projects/geo/) or ArrayExpress (http://www.ebi.ac.uk/arrayexpress/) databases under series numbers GSE42886 or E-MEXP-3779, respectively.

## Supporting Information

**Figure S1 Characterization of strains carrying chromosomally tagged alleles of SFL1 and SFL2.** (A) Strains SFL1-TAP (CEC1922), SFL2-TAP (CEC1918) and EFG1-HA (HLCEEFG1), carrying chromosomally tagged *SFL1* (tandem affinity purification tag, TAP), *SFL2* (tandem affinity purification tag, TAP) and *EFG1* (haemagglutinin tag, HA) alleles were grown in SC medium at 30°C or Lee's medium at 37°C during 4 h together with the SC5314 strain as a control (CTRL) prior to microscopic examination (40× magnification). (B) Western blot (WB) analyses of strains SFL1-TAP, SFL2-TAP (upper panel) and EFG1-HA (lower panel) together with the SC5314 control strain (CTRL). Strains were grown in SC medium at 30°C (30°C) or in Lee's medium at 37°C (37°C) during 4 h and total protein extracts were prepared then subjected to SDS-PAGE. Western blotting was performed using an anti-TAP antibody (SFL1-TAP and SFL2-TAP, Peroxydase-Anti-Peroxydase Soluble complex, Roche) or an anti HA antibody (EFG1-HA, Monoclonal Anti-HA peroxidase conjugate - High Affinity (clone 3F10, Roche). Positions of the molecular mass standards are indicated on the left (kDa). Antibody cross-reacting signals were used as a loading control (Loading Control). (TIF)

**Text S1 Includes Tables S1–S9 and full description of Tables S1–S9.**

(XLSX)

## Acknowledgments

We are grateful to members of the Unité Biologie et Pathogénicité Fongiques for their constant support and numerous insights during the course of this project. We are indebted to Caroline Proux and Jean-Yves Coppée from the Institut Pasteur Transcriptomics and Epigenomics Platform (Génopole Institut Pasteur) for their assistance with the ChIP-Seq experiments; Dr Martine Raymond from the Institute for Research in Immunology and Cancer (IRIC), Montréal, Canada for the SGY243-strain derivatives and the pCaMPY-3×HA plasmid and Dr Joachim Ernst from the Heinrich-Heine-Universität, Dusseldorf, Germany for strains AVL12 and HLCEEFG1. We also thank the Candida Genome Database for
sequence information and Drs Gaelle Lelandais, Jawad Merhej, Frédéric Devaux and Emmanuelle Permal for stimulating discussions.

## References

1. Pfaller MA, Diekema DJ (2007) Epidemiology of invasive mycoses in North America. Crit Rev Microbiol 36: 1-53.

2. Biswas S, Van Dijck P, Datta A (2007) Environmental sensing and signal transduction pathways regulating morphopathogenic determinants of *Candida albicans*. Microbiol Mol Biol Rev 71: 348-376.

3. Sudbery P, Gow N, Berman J (2004) The distinct morphogenic states of *Candida albicans*. Trends Microbiol 12: 317-324.

4. Sudbery PE (2011) Growth of *Candida albicans* hyphae. Nat Rev Microbiol 9: 737-748.

5. Whiteway M, Bachewich C (2007) Morphogenesis in *Candida albicans*. AnnuRevMicrobiol 61: 529-553.

6. Gow NA, Brown AJ, Odds FC (2002) Fungal morphogenesis and host invasion. Curr Opin Microbiol 5: 366-371.

7. Noble SM, French S, Kohn LA, Chen V, Johnson AD (2010) Systematic screens of a *Candida albicans* homozygous deletion library decouple morphogenetic switching and pathogenicity. Nat Genet 42: 590-598.

8. Saville SP, Lazzell AL, Monteagudo C, Lopez-Ribot JL (2003) Engineered control of cell morphology *in vivo* reveals distinct roles for yeast and filamentous forms of *Candida albicans* during infection. Eukaryot Cell 2: 1053-1060.

9. Berman J (2006) Morphogenesis and cell cycle progression in *Candida albicans*. Curr Opin Microbiol 9: 595-601.

10. Cottier F, Muhlschlegel FA (2009) Sensing the environment: response of *Candida albicans* to the X factor. FEMS Microbiol Lett 295: 1-9.

11. Davis-Hanna A, Piispanen AE, Stateva LI, Hogan DA (2008) Farnesol and dodecanol effects on the *Candida albicans* Ras1-cAMP signalling pathway and the regulation of morphogenesis. Mol Microbiol 67: 47-62.

12. Enjalbert B, Whiteway M (2005) Release from quorum-sensing molecules triggers hyphal formation during *Candida albicans* resumption of growth. Eukaryot Cell 4: 1203-1210.

13. Shapiro RS, Robbins N, Cowen LE (2011) Regulatory circuitry governing fungal development, drug resistance, and disease. Microbiol Mol Biol Rev 75: 213-267.

14. Cao F, Lane S, Raniga PP, Lu Y, Zhou Z, et al. (2006) The Flo8 transcription factor is essential for hyphal development and virulence in *Candida albicans*. Mol Biol Cell 17: 295-307.

15. Doedt T, Krishnamurthy S, Bockmuhl DP, Tebarth B, Stempel C, et al. (2004) APSES proteins regulate morphogenesis and metabolism in *Candida albicans*. Mol Biol Cell 15: 3167-3180.

16. Giusani AD, Vinces M, Kumamoto CA (2002) Invasive filamentous growth of *Candida albicans* is promoted by Czflp-dependent relief of Efg1p-mediated repression. Genetics 160: 1749-1753.

17. Lo HJ, Kohler JR, DiDomenico B, Loebenberg D, Cacciapuoti A, et al. (1997) Nonfilamentous *C. albicans* mutants are avirulent. Cell 90: 939-949.

18. Noffz CS, Liedschulte V, Lengeler K, Ernst JF (2008) Functional mapping of the *Candida albicans* Efg1 regulator. Eukaryot Cell 7: 881-893.

19. Stoldt VR, Sonneborn A, Leuker CE, Ernst JF (1997) Efg1p, an essential regulator of morphogenesis of the human pathogen *Candida albicans*, is a member of a conserved class of bHLH proteins regulating morphogenetic processes in fungi. EMBO J 16: 1982-1991.

20. Wang Y (2009) CDKs and the yeast-hyphal decision. Curr Opin Microbiol 12: 644-649.

21. Schweizer A, Rupp S, Taylor BN, Rollinghoff M, Schroppel K (2000) The TEA/ATTS transcription factor CaTec1p regulates hyphal development and virulence in *Candida albicans*. Mol Microbiol 38: 435-445.

22. Banerjee M, Thompson DS, Lazzell A, Carlisle PL, Pierce C, et al. (2008) UME6, a novel filament-specific regulator of *Candida albicans* hyphal extension and virulence. MolBiolCell 19: 1354-1365.

23. Zeidler U, Lettner T, Lassnig C, Muller M, Lajko R, et al. (2009) UME6 is a crucial downstream target of other transcriptional regulators of true hyphal development in *Candida albicans*. FEMS Yeast Res 9: 126-142.

24. Braun BR, Johnson AD (1997) Control of filament formation in *Candida albicans* by the transcriptional repressor TUP1. Science 277: 105-109.

25. Cleary IA, Mulabagal P, Reinhard SM, Yadev NP, Murdoch C, et al. (2010) Pseudohyphal regulation by the transcription factor Rfg1p in *Candida albicans*. Eukaryot Cell 9: 1363-1373.

26. Kadosh D, Johnson AD (2005) Induction of the *Candida albicans* filamentous growth program by relief of transcriptional repression: a genome-wide analysis. Mol Biol Cell 16: 2903-2912.

27. Murad AM, d'Enfert C, Gaillardin C, Tournu H, Tekaia F, et al. (2001) Transcript profiling in *Candida albicans* reveals new cellular functions for the transcriptional repressors CaTup1, CaMig1 and CaNrg1. Mol Microbiol 42: 981-993.

28. Borneman AR, Leigh-Bell JA, Yu H, Bertone P, Gerstein M, et al. (2006) Target hub proteins serve as master regulators of development in yeast. Genes Dev 20: 435-448.

## Author Contributions

Conceived and designed the experiments: SZ AN TR CdE. Performed the experiments: SZ AN MC TR. Analyzed the data: SZ AN TR CdE. Wrote the paper: SZ AN CdE.

29. Ryan O, Shapiro RS, Kurat CF, Mayhew D, Baryshnikova A, et al. (2012) Global gene deletion analysis exploring yeast filamentous growth. Science 337: 1353-1356.

30. Pan X, Heitman J (2002) Protein kinase A operates a molecular switch that governs yeast pseudohyphal differentiation. Mol Cell Biol 22: 3981-3993.

31. Robertson LS, Fink GR (1998) The three yeast A kinases have specific signaling functions in pseudohyphal growth. Proc Natl Acad Sci U S A 95: 13783-13787.

32. Fujita A, Kikuchi Y, Kuhara S, Misumi Y, Matsumoto S, et al. (1989) Domains of the SFL1 protein of yeasts are homologous to Myc oncoproteins or yeast heat-shock transcription factor. Gene 85: 321-328.

33. Zhu C, Byers KJ, Cordell RP, Shi Z, Berger MF, et al. (2009) High-resolution DNA-binding specificity analysis of yeast transcription factors. Genome Res 19: 556-566.

34. Yamamoto A, Mizukami Y, Sakurai H (2005) Identification of a novel class of target genes and a novel type of binding sequence of heat shock transcription factor in *Saccharomyces cerevisiae*. J Biol Chem 280: 11911-11919.

35. Conlan RS, Tzamarias D (2001) Sfl1 functions via the co-repressor Ssn6-Tup1 and the cAMP-dependent protein kinase Tpk2. J Mol Biol 309: 1007-1015.

36. Galeote VA, Alexandre H, Bach B, Delobel P, Dequin S, et al. (2007) Sfl1 acts as an activator of the HSP30 gene in *Saccharomyces cerevisiae*. Curr Genet 52: 55-63.

37. Bauer J, Wendland J (2007) *Candida albicans* Sfl1 suppresses flocculation and filamentation. Eukaryot Cell 6: 1736-1744.

38. Li Y, Su C, Mao X, Cao F, Chen J (2007) Roles of *Candida albicans* Sfl1 in hyphal development. Eukaryot Cell 6: 2112-2121.

39. Song W, Wang H, Chen J (2011) *Candida albicans* Sfl2, a temperature-induced transcriptional regulator, is required for virulence in a murine gastrointestinal infection model. FEMS Yeast Res 11: 209-222.

40. Spiering MJ, Moran GP, Chauvel M, Maccallum DM, Higgins J, et al. (2010) Comparative transcript profiling of *Candida albicans* and *Candida dubliniensis* identifies SFL2, a *C. albicans* gene required for virulence in a reconstituted epithelial infection model. Eukaryot Cell 9: 251-265.

41. Chauvel M, Nesseir A, Cabral V, Znaidi S, Goyard S, et al. (2012) A Versatile Overexpression Strategy in the Pathogenic Yeast *Candida albicans*: Identification of Regulators of Morphogenesis and Fitness. PLoS One 7: e45912.

42. Care RS, Trevedickh J, Binley KM, Sudbery PE (1999) The MET3 promoter: a new tool for *Candida albicans* molecular genetics. Mol Microbiol 34: 792-798.

43. Znaidi S, Barker KS, Weber S, Alarco AM, Liu TT, et al. (2009) Identification of the *Candida albicans* Cap1p regulon. Eukaryot Cell 8: 806-820.

44. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, et al. (2011) Integrative genomics viewer. Nat Biotechnol 29: 24-26.

45. Thorvaldsdottir H, Robinson JT, Mesirov JP (2012) Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform 14:178-92.

46. Feng J, Liu T, Zhang Y (2011) Using MACS to identify peaks from ChIP-Seq data. Curr Protoc Bioinformatics Chapter 2: Unit 2 14.

47. Chen Y, Negre N, Li Q, Mieczkowska JO, Slattery M, et al. (2012) Systematic evaluation of factors influencing ChIP-seq fidelity. Nat Methods 9: 609-614.

48. Arnaud MB, Costanzo MC, Shah P, Skrzpek MS, Sherlock G (2009) Gene Ontology and the annotation of pathogen genomes: the case of *Candida albicans*. Trends Microbiol 17: 295-303.

49. Shapiro RS, Sellam A, Tebbji F, Whiteway M, Nantel A, et al. (2012) Pho85, Pc11, and Hms1 signaling governs *Candida albicans* morphogenesis induced by high temperature or Hsp90 compromise. Curr Biol 22: 461-470.

50. Homann OR, Dea J, Noble SM, Johnson AD (2009) A phenotypic profile of the *Candida albicans* regulatory network. PLoS Genet 5: e1000783.

51. Lassak T, Schneider E, Bussmann M, Kurtz D, Manak JR, et al. (2011) Target specificity of the *Candida albicans* Efg1 regulator. Mol Microbiol 82: 602-618.

52. Finkel JS, Xu W, Huang D, Hill EM, Desai JV, et al. (2012) Portrait of *Candida albicans* adherence regulators. PLoS Pathog 8: e1002525.

53. Hao B, Clancy CJ, Cheng S, Raman SB, Iczkowski KA, et al. (2009) *Candida albicans* RFX2 encodes a DNA binding protein involved in DNA damage responses, morphogenesis, and virulence. Eukaryot Cell 8: 627-639.

54. Nobile CJ, Fox EP, Nett JE, Sorrells TR, Mitrovich QM, et al. (2012) A recently evolved transcriptional network controls biofilm development in *Candida albicans*. Cell 148: 126-138.

55. Thomas-Chollier M, Herrmann C, Defrance M, Sand O, Thieffry D, et al. (2012) RSAT peak-motifs: motif analysis in full-size ChIP-seq datasets. Nucleic Acids Res 40: e31.

56. Carlson JM, Chakravarty A, DeZiel CE, Gross RH (2007) SCOPE: a web server for practical *de novo* motif discovery. Nucleic Acids Research 35: W259-W264.

57. Sellam A, Tebbji F, Nantel A (2009) Role of Ndt80p in sterol metabolism regulation and azole resistance in *Candida albicans*. Eukaryot Cell 8: 1174-1183.
58. Sellam A, Askew C, Epp E, Tebbji F, Mullick A, et al. (2010) Role of transcription factor CaNdt80p in cell separation, hyphal growth, and virulence in *Candida albicans*. *Eukaryot Cell* 9: 634–644.

59. Bryne JC, Valen E, Tang MH, Marstrand T, Winther O, et al. (2008) JASPAR, the open access database of transcription factor-binding profiles: new content and tools in the 2008 update. *Nucleic Acids Res* 36: D102–106.

60. de Boer CG, Hughes TR (2012) YeTFaSCo: a database of evaluated yeast transcription factor sequence specificities. *Nucleic Acids Res* 40: D169–179.

61. Akerfelt M, Morimoto RI, Sistonen L (2010) Heat shock factors: integrators of cell stress, development and lifespan. *Nat Rev Mol Cell Biol* 11: 545–555.

62. Nicholls S, Leach MD, Priest CL, Brown AJ (2009) Role of the heat shock transcription factor, Hsf1, in a major fungal pathogen that is obligately associated with warm-blooded animals. *Mol Microbiol* 74: 844–861.

63. Nicholls S, MacCallum DM, Kaffarnik FA, Selway L, Peck SC, et al. (2011) Activation of the heat shock transcription factor Hsf1 is essential for the full virulence of the fungal pathogen *Candida albicans*. *Fungal Genet Biol* 48: 297–305.

64. MacIsaac KD, Wang T, Gordon DB, Gifford DK, Stormo GD, et al. (2006) An improved map of conserved regulatory sites for *Saccharomyces cerevisiae*. *BMC Bioinformatics* 7: 113.

65. Morozov AV, Siggia ED (2007) Connecting protein structure with predictions of regulatory sites. *Proc Natl Acad Sci U S A* 104: 7068–7073.

66. Zhang S, Burkett TJ, Yamashita I, Garfinkel DJ (1997) Genetic redundancy between *SPT23* and *MGA2*: regulators of Ty-induced mutations and Ty1 transcription in *Saccharomyces cerevisiae*. *Mol Cell Biol* 17: 4718–4729.

67. Lorenz MC, Heitman J (1998) Regulators of pseudohyphal differentiation in *Saccharomyces cerevisiae* identified through multicopy suppressor analysis in ammonium permease mutant strains. *Genetics* 150: 1443–1457.

68. Borneman AR, Gianoulis TA, Zhang ZD, Yu H, Rozowsky J, et al. (2007) Divergence of transcription factor binding sites across related yeast species. *Science* 317: 815–819.

69. Sherman F (1991) Getting started with yeast. *Methods in Enzymology* 194: 3–21.

70. Lee KL, Buckley HR, Campbell CC (1975) An amino acid liquid synthetic medium for the development of mycelial and yeast forms of *Candida Albicans*. *Sabouraudia* 13: 148–153.

71. Walther A, Wendland J (2003) An improved transformation protocol for the human fungal pathogen *Candida albicans*. *Curr Genet* 42: 339–343.

72. Reuss O, Vik A, Kolter R, Morschhauser J (2004) The *SAT1* flipper, an optimized tool for gene disruption in *Candida albicans*. *Gene* 341: 119–127.

73. Liu TT, Znaidi S, Barker KS, Xu L, Homayouni R, et al. (2007) Genome-wide expression and location analyses of the *Candida albicans* Tac1p regulon. *Eukaryot Cell* 6: 2122–2138.

74. Gola S, Martin R, Walther A, Dunkler A, Wendland J (2003) New modules for PCR-based gene targeting in *Candida albicans*: rapid and efficient gene targeting using 100 bp of flanking homology region. *Yeast* 20: 1339–1347.

75. Schaub Y, Dunkler A, Walther A, Wendland J (2006) New pFA-cassettes for PCR-based gene manipulation in *Candida albicans*. *J Basic Microbiol* 46: 416–429.

76. Murad AM, Lee PR, Broadbent ID, Barelle CJ, Brown AJ (2000) CIP10, an efficient and convenient integrating vector for *Candida albicans*. *Yeast* 16: 325–327.

77. Blankenberg D, Von Kuster G, Coraor N, Ananda G, Lazarus R, et al. (2010) Galaxy: a web-based genome analysis tool for experimentalists. *Curr Protoc Mol Biol Chapter 19: Unit 19 10 11–21*.

78. Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, et al. (2005) Galaxy: a platform for interactive large-scale genome analysis. *Genome Res* 15: 1451–1455.

79. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. *Genome Biol* 10: R25.

80. Salmon-Divon M, Dvinge H, Tammoja K, Bertone P (2010) PeakAnalyzer: genome-wide annotation of chromatin binding and modification loci. *BMC Bioinformatics* 11: 415.

81. Wodicka L, Dong H, Mittmann M, Ho MH, Lockhart DJ (1997) Genome-wide expression monitoring in *Saccharomyces cerevisiae*. *NatBiotechnol* 15: 1359–1367.

82. Hokamp K, Roche FM, Acab M, Rousseau ME, Kuo B, et al. (2004) ArrayPipe: a flexible processing pipeline for microarray data. *Nucleic Acids Res* 32: W457–459.

83. Sturn A, Quackenbush J, Trajanoski Z (2002) Genesis: cluster analysis of microarray data. *Bioinformatics* 18: 207–208.

84. Gillum AM, Tsay EY, Kirsch DR (1984) Isolation of the *Candida albicans* gene for orotidine-5'-phosphate decarboxylase by complementation of *S. cerevisiae ura3* and *E. coli* pyrF mutations. *Molecular & General Genetics* 198: 179–182.

85. Fonzi WA, Irwin MY (1993) Isogenic strain construction and gene mapping in *Candida albicans*. *Genetics* 134: 717–728.

86. Firon A, Aubert S, Iraqui I, Guadagnini S, Goyard S, et al. (2007) The *SUN41* and *SUN42* genes are essential for cell separation in *Candida albicans*. *Mol Microbiol* 66: 1256–1275.

87. Noble SM, Johnson AD (2005) Strains and strategies for large-scale gene deletion studies of the diploid human fungal pathogen *Candida albicans*. *Eukaryot Cell* 4: 298–309.