# Identification of DNA binding motifs of the Mycobacterium tuberculosis PhoP/PhoR two-component signal transduction system.

Mena Cimino, Christophe Thomas, Amine Namouchi, Sarah Dubrac, Brigitte Gicquel, Deshmukh N Gopaul

► To cite this version:

Mena Cimino, Christophe Thomas, Amine Namouchi, Sarah Dubrac, Brigitte Gicquel, et al.. Identification of DNA binding motifs of the Mycobacterium tuberculosis PhoP/PhoR two-component signal transduction system.. PLoS ONE, 2012, 7 (8), pp.e42876. 10.1371/journal.pone.0042876. pasteur-01234357

HAL Id: pasteur-01234357

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

Submitted on 26 Nov 2015

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>
# Identification of DNA Binding Motifs of the *Mycobacterium tuberculosis* PhoP/PhoR Two-Component Signal Transduction System

Mena Cimino¹,³, Christophe Thomas²,⁵, Amine Namouchi¹,⁵, Sarah Dubrac³, Brigitte Gicquel¹,*,
Deshmukh N. Gopaul²¶

¹Unité de Génétique Mycobacterienne, Institut Pasteur, Paris, France, ²Unité de Plasticité du Génome bactérien, Institut Pasteur, CNRS UMR 3525, Paris, France, ³Unité de Recherche Biologie des Bactéries pathogènes à Gram-positif, Institut Pasteur, Paris, France

## Abstract

**Background:** The *Mycobacterium tuberculosis* PhoP/PhoR two-component signal transduction system controls the expression of about 2% of the genome and plays a major role in pathogenicity. However, its regulon has not been well characterized.

**Methodology/Principal Findings:** The binding site of PhoP transcription regulator was identified in the upstream regions of msl3, pks2, lipF and fadD21 genes, by using gene fusions, electrophoretic mobility shift assays and DNase I footprinting experiments. A consensus sequence for PhoP binding was deduced. It consists of two direct repeats, DR1/DR2, associated with a third repeat, DR3, important in some cases for PhoP binding to DR1/DR2 but located at a variable distance from these direct repeats. DR1/DR2 and DR3 consensus sequences were used to screen the whole-genome sequence for other putative binding sites potentially corresponding to genes directly regulated by PhoP. The identified 87 genes, encoding transcription regulators, and proteins involved in secondary metabolites biosynthesis, transport and catabolism are proposed to belong to the PhoP regulon.

**Conclusions/Significance:** A consensus sequence derived from the analysis of PhoP binding to four gene promoter regions is proposed. We show for the first time the involvement of a third direct repeat motif in this binding reaction. The consensus sequence was instrumented to study the global regulation mediated by PhoP in *M. tuberculosis*. This analysis leads to the identification of several genes that are potentially regulated by this key player.

**Citation:** Cimino M, Thomas C, Namouchi A, Dubrac S, Gicquel B, et al. (2012) Identification of DNA Binding Motifs of the *Mycobacterium tuberculosis* PhoP/PhoR Two-Component Signal Transduction System. PLoS ONE 7(8): e42876. doi:10.1371/journal.pone.0042876

**Editor:** Petros C. Karakousis, Johns Hopkins University School of Medicine, United States of America

**Received:** March 15, 2012; **Accepted:** July 13, 2012; **Published:** August 7, 2012

**Copyright:** © 2012 Cimino 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 study received funding from the 7th Framework Program NEWTBVAC contract, HEALTH-F3-2009-241745, the FP7-REGPOT EUMEDNET vs TB project, the ANR-09-MIEN-024-02 project and the Institut Pasteur Starter grant/CNRS URA 2185. 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

*Mycobacterium tuberculosis* (MTB), the causal agent of tuberculosis in humans, is one of the leading causes of mortality due to a single infectious agent. MTB is a successful intracellular pathogen that can adapt to changing environmental conditions within the host. It can infect different cell types, in which it may replicate or remain dormant for years. Differential expression of MTB genes is observed during the infection of various cells of the immune system, such as dendritic cells, macrophages and alveolar epithelial cells, during disease progression [1,2,3,4]. We have previously shown that bacterial stress responses are more strongly induced in dendritic cells, whereas genes encoding ribosomal proteins are overexpressed, indicating bacterial multiplication, in macrophages [5].

As in other prokaryotes, two-component signal transduction systems (TCS) are key elements of the adaptive response to various stimuli in the tubercle bacillus. TCS contain an environmentally sensitive histidine kinase (HK) and a response regulator (RR) that is activated by the cognate HK. These systems play a major role in bacterial responses to changing growth contexts. The MTB genome encodes 11 TCS [6]. This relatively small number of TCS probably reflects the intracellular lifestyle of MTB, as the cell environment is less variable than that confronted by soil bacteria, or a certain degree of overlap in signal processing.

The MTB *phoP* gene has been shown to encode one of the components of a TCS (PhoP/R) playing a major role in virulence [7]. Inactivation of *phoP* results in an attenuated mutant unable to replicate in animal models and in cells cultured *in vitro*, but able to persist, unlike auxotrophs. This led to the construction of live candidate vaccines against tuberculosis based on the inactivation

<footnote>* E-mail: brigitte.gicquel@pasteur.fr</footnote>

<footnote>¶ These authors contributed equally to this work.</footnote>

<footnote>‡‡ These authors also contributed equally to this work.</footnote>
of *phoP* [8]. Analyses of the genomes of known avirulent strains, such as the BCG vaccine or the H37Ra attenuated clinical variant, showed mutations in *phoP* or *phoR* resulting in a loss of function, thus confirming the role of *phoP/phoR* in virulence [9,10,11].

The MTB PhoP regulator belongs to the PhoB/OmpR subfamily, the largest one among response regulators. Members of this subfamily have two domains, an N-terminal regulatory domain and a C-terminal DNA-binding domain (also called the effector domain). Many members of this subfamily have been studied in detail, and most reports have indicated that these response regulators bind DNA as a dimer recognizing tandem repeat binding sites at the -35 region of the regulated promoter [12]. The *Bacillus subtilis* TCS PhoP/PhoR, which belongs to this subfamily, senses phosphate. It is activated in conditions in which phosphate is limited. It regulates the expression of 31 genes involved in phosphate utilization [13]. However, in *Salmonella enterica typhimurium*, PhoP/PhoQ responds to other stimuli, such as $Mg^{2+}$ concentration [14], low pH [15,16] and antibacterial peptides [17,18], thereby regulating genes involved in $Mg^{2+}$ homeostasis and virulence [19]. The stimuli sensed by the MTB PhoP/PhoR TCS are unknown.

Efforts have been made to identify genes regulated by this PhoP/PhoR TCS, by comparing the transcriptomes of *phoP* null mutants and wild-type parental strains. Two studies identified a large number of genes positively or negatively regulated by PhoP, 70 of which were upregulated [20,21]. However, these two studies identified different, though overlapping, sets of regulated genes. Differences could also be due to the use of two different strains in the two studies (H37Rv for one study and MT103 for the other one), or to the lack of knowledge of the metabolites acting as stimuli for this TCS, resulting in non optimal experimental conditions.

Many of the genes upregulated by PhoP are involved in general or lipid metabolism, substrate transport across the plasma membrane and the synthesis of regulators, such as DosR, controlling dormancy. *phoP* null mutants display deficiencies in the synthesis of sulfatides and diacyl and polyacyl trehaloses. *msl3*, a polyketide beta-ketoacyl synthase gene, is involved in the synthesis of polyacyl trehaloses, and the *pks2* and *mmp18* genes are involved in sulfatide synthesis [21,22]. The disruption of *pks2* generated a sulfolipid-deficient mutant that was unable to synthesize hydroxyphthioceranic and phthioceranic acids [23]. The *msl3*-disrupted mutant was unable to produce the mycolipanoic and mycolipenic acids required for the synthesis of a major class of polyacylated trehaloses. In the absence of these classes of polyacylated trehaloses, which anchor trehaloses to the cell surface, the mutants grew in bead-like aggregates, with no discernable decrease in either growth rate or in virulence [24]. *phoP* null mutants and H37Ra lack these complex lipids. PhoP/PhoR TCS also controls secretion of the early-secreted 6 kDa antigen (ESAT-6), an important virulence factor and antigenic component of *M. tuberculosis* [9].

The two transcriptomic analyses suggested that the *lipF* and *fadD21* genes were positively regulated by PhoP [20,21]. The *lipF* gene of MTB encodes an esterase that has been shown to be important in pathogenesis. The insertion of a transposon between the *lipF* promoter region and the transcription start site significantly decreased the ability of the bacterium to grow in mouse lungs [25].

PhoP is autoregulated, as it binds to two adjacent 9-bp direct repeat motifs binding to PhoP in a sequence-specific manner. These motifs are located downstream from the transcription +1 site, contrasting with the models proposed for other members of the PhoB/OmpR TCS family [26]. We investigated PhoP binding

sites and their molecular regulation in more detail, by looking for such sites in other promoters thought to be regulated by PhoP on the basis of previous transcriptome analyses.

We first used gene fusions to identify the *msl3*, *pks2*, *lipF* and *fadD21* gene promoters. DNA-affinity electrophoretic mobility shift studies and DNase I protection assays have been used to identify PhoP binding sites. These sites were found to contain tandem repeat sequences displaying similarities that could be used to define a consensus-binding motif. We then carried out whole-genome analysis with these consensus sequences, to identify genes regulated directly by PhoP.

## Results

### Identification of promoter regions

We used *lacZ* as a reporter gene, to measure the activity of *M. tuberculosis* promoters thought to be regulated by PhoP in *M. smegmatis*, a non-pathogenic and fast-growing mycobacterial species.

Several fragments corresponding to upstream and structural gene regions of *pks2* (-86,+60), *msl3* (-282,+28) and *fadD21* (-198, +10) were generated (Table 1). These fragments were inserted upstream from the *lacZ* reporter gene in the pJEM15 *E. coli*-mycobacterial reporter shuttle plasmid. For *lipF*, a 611 bp region (-596, +14) was used, because this region has previously been reported to be required for the upregulation of transcriptional activity in response to exposure to an acid stress. This upregulation has been demonstrated in both pathogenic *M. tuberculosis* and non-pathogenic *M. smegmatis* strains [27]. Cultures of *M. smegmatis* transformed with the empty pJEM15 vector, a positive control, the pJEM31 (PAN) vector containing a promoter sequence, PAN, isolated from *M. paratuberculosis* [28] or recombinant pJEM15 plasmids carrying *lipF*, *pks2*, *msl3* or *fadD21* gene fragments were set up. Aliquots of cultures were collected during the exponential growth phase and beta-galactosidase activities were assessed to determine transcription levels.

As indicated in Fig. 1, expression levels were higher for *M. smegmatis* transformed with the recombinant plasmids carrying *pks2*, *msl3*, *lipF* or *fadD21* gene fusions than for *M. smegmatis* transformed with the empty pJEM15 vector. No significant expression was detected for *M. smegmatis* transformed with the empty pJEM15 vector. With the plasmid pJEM31 (PAN) used as a positive control, we observed expression of similar magnitude as observed for *pks2*, *msl3* and *fadD21*. A much lower expression was observed for *M. smegmatis* carrying the *lipF* fusion, under the conditions tested. These results demonstrate that there was a promoter present in the various cloned regions in fusion with *lacZ* gene.

### Identification of PhoP binding sites

We investigated whether PhoP bound directly to the DNA sequences upstream from *lipF*, *pks2*, *msl3* and *fadD21* identified in gene fusion studies, by carrying out electrophoretic mobility shift assays (EMSA) with phosphorylated PhoP protein (PhoP-P). In fact, previous studies have shown that PhoP must be phosphorylated for efficient binding to specific sites. We also carried out EMSA with the regulatory binding region of the PhoP gene itself, Protect40, as a positive control [29,30]. Only positive bindings were shown in Fig. 2. We started by using PCR fragments of about 200 bp in size for EMSA (Table 1, Table 2, and Fig. 2A). For the *lipF* promoter region, we generated three 200 bp subfragments: *lipFa* (-596 to -375), *lipFb* (-394 to -168) and *lipFc* (-167 to +14). Only *lipFa* displayed a clear shift in electrophoretic mobility in the presence of PhoP-P (Fig. 2A and Fig.S1). For *msl3*, we
<table><caption>**Table 1.** List of primers used in this study.</caption><thead><tr><th>Fragments</th><th>Primers</th><th>Sequence (5'-3' position relative to the translation start site)</th><th>Fragments size (bp)</th></tr></thead><tbody><tr><td rowspan="2">lipF<sup>a</sup></td><td>lipFF</td><td>(-596) GGCGAATTCGCACTTACCA (-578)</td><td rowspan="2">611</td></tr><tr><td>lipFR</td><td>(+14) CCAGGCGCACGCACCAAC (-3)</td></tr><tr><td rowspan="2">lipFa<sup>b</sup></td><td>lipFaF</td><td>(-596) GGCGAATTCGCACTTACCA (-578)</td><td rowspan="2">222</td></tr><tr><td>lipFaR</td><td>(-375) ATCAATCGCAAATCTCGCAC (-394)</td></tr><tr><td rowspan="2">lipFb<sup>b</sup></td><td>lipFbF</td><td>(-394) GTGCGAGATTTGCGATTGAT (-375)</td><td rowspan="2">227</td></tr><tr><td>lipFbR</td><td>(-168) CAACCGACCTGCAACACTA (-186)</td></tr><tr><td rowspan="2">lipFc<sup>b</sup></td><td>lipFcF</td><td>(-167) GCGGCACTGGCATCGCGCAT (-148)</td><td rowspan="2">182</td></tr><tr><td>lipFcR</td><td>(+14) CCAGGCGCACGCACCAAC (-3)</td></tr><tr><td rowspan="2">lipFa2<sup>b</sup></td><td>lipFa2F</td><td>(-465) CACCGAATGATTGATGG (-449)</td><td rowspan="2">89</td></tr><tr><td>lipFa2R</td><td>(-377) CAATCGCAAATCTCGCA (-393)</td></tr><tr><td rowspan="2">lipF<sup>c</sup></td><td>lipFDF</td><td>(-639) GCCGATATCGGCTTGCGATTC (-619)</td><td rowspan="2">189</td></tr><tr><td>lipFDR</td><td>(-470) CATCAATCATTCGGTGGCGCG (-450)</td></tr><tr><td rowspan="2">pks2<sup>ab</sup></td><td>pks2F</td><td>(-86) TAGCACAGCCGCTTAGAAC (-68)</td><td rowspan="2">147</td></tr><tr><td>pks2R</td><td>(+60) TCAGCCAACGTCCATGC (+44)</td></tr><tr><td rowspan="2">pks2c<sup>b</sup></td><td>pks2cF</td><td>(+10) GTCGGCGGCGTCAGGCACT (+28)</td><td rowspan="2">77</td></tr><tr><td>pks2cR</td><td>(+86) TAACAGCAACCGGAGTCA (+69)</td></tr><tr><td rowspan="2">pks2d<sup>b</sup></td><td>pks2dF</td><td>(-86) TAGCACAGCCGCTTAGAAC (-68)</td><td rowspan="2">102</td></tr><tr><td>pks2dR</td><td>(+15) GCCGACCCCAAGCCCAAT (-2)</td></tr><tr><td rowspan="2">pks2<sup>ab</sup></td><td>pks2DF</td><td>(-137) ACCGCGGGCTGCCCATGCCA (-118)</td><td rowspan="2">198</td></tr><tr><td>pks2DR</td><td>(+60) TCAGCCAACGTCCATGCA (+43)</td></tr><tr><td rowspan="2">msl3<sup>a</sup></td><td>msl3F</td><td>(-282) GTCGAGCAGTGTTGCTACC (-264)</td><td rowspan="2">311</td></tr><tr><td>msl3R</td><td>(+28) CGACCGATGTGGCGGTTGCGG (+8)</td></tr><tr><td rowspan="2">msl3a<sup>b</sup></td><td>msl3aF</td><td>(-282) GTCGAGCAGTGTTGCTACC (-264)</td><td rowspan="2">234</td></tr><tr><td>msl3aR</td><td>(-49) TACCAACACATTCGGGCTC (-67)</td></tr><tr><td rowspan="2">msl3b<sup>b</sup></td><td>msl3bF</td><td>(-67) GAGCCCGAATGTGTTGGTA (-49)</td><td rowspan="2">98</td></tr><tr><td>msl3bR</td><td>(+30) CGACCGATGTGGCGGTTGCGG (+10)</td></tr><tr><td rowspan="2">msl3a2<sup>b</sup></td><td>msl3a2F</td><td>(-269) GCTACCTTAACTTTCCCAG (-251)</td><td rowspan="2">82</td></tr><tr><td>msl3a2R</td><td>(-188) CGAAGAACACATTATCCAG (-206)</td></tr><tr><td rowspan="2">msl3<sup>c</sup></td><td>msl3DF</td><td>(-311) TCACACGCGCGGTGCATGC (-293)</td><td rowspan="2">222</td></tr><tr><td>msl3DR</td><td>(-90) TCATGGTAGTGGCTCAAGA (-108)</td></tr><tr><td rowspan="2">fadD21<sup>ab</sup></td><td>fadD21F</td><td>(-198) GGCCATGGAACGCACGGCTTA (-178)</td><td rowspan="2">209</td></tr><tr><td>fadD21R</td><td>(+10) GAGTCGGACATTGGTGCTA (-8)</td></tr><tr><td rowspan="2">fadD21b<sup>b</sup></td><td>fadD21bF</td><td>(-145) CTGCGTACACCGACTTCGA (-127)</td><td rowspan="2">85</td></tr><tr><td>fadD21bR</td><td>(-61) GCGATAGCGATGCTGATC (-78)</td></tr><tr><td rowspan="2">fadD21<sup>c</sup></td><td>fadD21DF</td><td>(-194) ATGGAACGCACGGCTTAGCAACT (-172)</td><td rowspan="2">201</td></tr><tr><td>fadD21DR</td><td>(-16) CGGACATTGGTGCTACATTACCG (+7)</td></tr></tbody></table>

generated two subfragments: msl3a (-282 to -49), which was about 200 bp in length, and msl3b (-67 to +28), which was about 100 bp long. Only msl3a displayed a clear shift in electrophoretic mobility in the presence of PhoP-P (Fig. 2A).

The 147 bp (-86 to +60) pks2 and the 209 bp (-198 to +10) fadD21 fragments used in beta-galactosidase experiments also displayed a clear shift in electrophoretic mobility in the presence of PhoP-P (Fig. 2A). To address relative affinity differences between the DNA substrates, we varied the protein to DNA ratio to obtain a complete shift. The protein concentrations used ranged from 0 to 8 µM, for a constant DNA concentration of 40 nM. All five

regulatory regions studied (lipF, pks2, msl3, fadD21 and phoP) bound to PhoP-P resulted in the detection of higher molecular weight bands. For msl3a, fadD21 and pks2, a PhoP-P:DNA ratio of 100:1 was sufficient to obtain an almost complete shift, whereas a ratio of 200:1 was required for lipFa. The various DNA segments thus have different affinities for PhoP-P (Fig. 2). A faint second band of higher molecular weight was observed for Protect-40, but not for the other targets analyzed here. All reactions were tested for the requirement for phosphorylation (Fig. 2). Non-phosphorylated assays are shown (Fig.S2). No significant binding was observed in absence of phosphorylation.

<footnote><sup>a</sup>Fragments used for transcriptional fusions.</footnote>

<footnote><sup>b</sup>Fragments used for electrophoretic mobility shift assays.</footnote>

<footnote><sup>c</sup>Fragments used for DNase protection assay.</footnote>

<footnote>doi:10.1371/journal.pone.0042876.t001</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Beta-galactosidase activity in *M. smegmatis* expressing *lacZ* under the control of the promoter regions of the *M. tuberculosis H37Rv* *lipF*, *pks2*, *msl3* and *fadD21* genes. The *lipF*, *pks2*, *msl3* and *fadD21* promoter-lacZ fusions carried by the pJEM15 plasmid were introduced into *M. smegmatis*. The positive control is *M. smegmatis* transformed with pJEM31(PAN) containing *lacZ* under the control of a mobile genetic element normalized against *M. smegmatis* transformed with pJEM15 empty vector. The results shown are the mean of two independent experiments. doi:10.1371/journal.pone.0042876.g001</figcaption></figure>

Using the data obtained with the large fragments, we carried out a sequence alignment analysis of all the fragments shown to bind PhoP-P. The regions similar to DR1/DR2 and DR3 previously identified in the *phoP* promoter [26] were observed.

We applied trimming and walking techniques to the initial oligonucleotides to generate subsequent fragments of 50 to 100 bp in size (Fig. 2B, Fig. 3) and proceeded with PhoP-P binding.

The results from the EMSAs on the subfragments were analyzed in the context of the presence or absence of the DR repeats. In addition to lipFa, we also tested two shorter DNA regions for binding: lipFa1 (–575 to –508) and lipFa2 (–465 to –377). Binding was observed only with lipFa1, which contains DR1–3, at a PhoP-P concentration of 2 µM, but not for lipFa2. For *pks2*, we tested four additional shorter DNA regions: pks2a (–52 to +6), pks2b (–95 to –56), pks2c (+10 to +86) and pks2d (–86 to +15). Binding was observed for pks2d, which contains a partial DR1 and DR2–3 (at a PhoP-P concentration of 2 µM). Only very weak binding was observed for pks2a, which contained a DR3 sequence alone (at a PhoP-P concentration of 10 µM), and for pks2b, which contained only DR1 and DR2 sequences (at a PhoP-P concentration of 20 µM). pks2c which does not contain any of the DRs did not give any binding (Fig. 2B). We studied *msl3a* binding in more detail with two additional DNA fragments: *msl3a1* (–242 to –198) and *msl3a2* (–269 to –188). Binding was observed for both of these fragments; *msl3a2* which contains the entire DR1–3 (binding observed at a PhoP-P concentration of 2 µM), whereas *msl3a1* contains only the first 3 bp of DR3 (binding at a PhoP-P concentration of 2 µM). Finally, *fadD21* was studied by testing two additional DNA regions: *fadD21a* (–152 to –98), and *fadD21b* (–145 to –61). Binding was observed for *fadD21b*, which contained all three DR sequences (PhoP-P concentration of 2 µM). Only very weak binding was observed for *fadD21a*, which contained only DR1 and partial DR2 (PhoP-P concentration of 10 µM), as for pks2b, which also contained only DR1–2 (PhoP-P concentration of 20 µM). However, in assays of PhoP-P binding to Protect40, full binding was observed when we used a DR1–2-containing fragment (PhoP-P concentration of

4 µM). Interestingly, reducing the size of the DNA targets (45–102 bp) resulted in a decrease in the protein/DNA ratio (only 1:50) required to obtain a shift, (Fig. 2B), showing that we indeed were able to increase the affinity as we narrowed down the targets to match the DRs.

The transcription initiation sites were determined for *lipF*, *pks2* and *msl3* in previous studies [31,32]. This allowed us to locate DR1, DR2 and DR3 with regard to the transcription start.

We investigated the importance of the repeat motifs in sequence-specific DNA binding by PhoP, using oligonucleotides with altered nucleotide sequences Pho1, Pho2, Pho3 (Table 2, Fig. 2C and 3E) in place of DR1, DR2 or combined DR1/DR2 sequences, respectively whereby we removed one or all of the DR1–2 repeats. The Protect-40 wild-type sequence was used as a control (Fig. 2C). The Pho1–Pho3 altered sequences severely impaired PhoP-P binding. These results confirm that the binding of PhoP to the *phoP* promoter region is dependent on the direct repeat motifs.

## Characterization of the PhoP binding site by DNase I footprinting

We defined the DNA regions binding PhoP more precisely, by carrying out DNase I protection assays on the PhoP-P/DNA target complex. We used the *lipF*, *pks2*, *fadD21* and *msl3* promoter regions for these protection experiments. PhoP-P binding enhanced the protection of these four promoters against DNase I digestion (Fig. 4). The protected regions for *pks2*, *msl3* and *fadD21* were approximately 50 to 70 bp long (Fig. 4). For *pks2* and *fadD21*, these regions include coverage of the DRs, the spacer region between them, with a complete DR3 motif. In the case of *msl3*, the DR3 is partially protected until position A7 known for being present in all the DRs (Table 3, [26]). The protected region of *lipF* is only 30–35 bp in length and the DR3 region is not protected by PhoP protein.
A

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

B

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

C

<figure><img src="image_5.png" /></figure>
**Figure 2.** Electrophoretic mobility assays with the *lipF*, *pks2*, *msl3* and *fadD21* promoter fragments and PhoP-P. 2A. The large DNA fragments initially selected — *lipFa* (222 bp), *pks2* (148 bp), *msl3a* (235 bp), and *fadD21* (210 bp) — were incubated with PhoP-P in the presence of poly dl-dC at 10 µg/ml and run on a native polyacrylamide gel (8%), in 0.5 × TBE buffer. Complete binding was assessed on the basis of the amount of shifted material for each fragment in the presence of 4 µM of PhoP-P (protein/DNA ratio: 100/1), except for *lipFa*, for which 8 µM PhoP-P was required (protein/DNA ratio: 200/1). 2B. Smaller DNA fragments — *lipFa1* (68 bp), *pks2a* (59 bp), *pks2b* (40 bp), *pks2d* (104 bp), *msl3a1* (45 bp), *msl3a2* (83 bp), *fadD21a* (55 bp) and *fadD21b* (79 bp) — were designed and incubated with PhoP-P in the presence of poly dl-dC at 10 µg/ml. Protect40 (40 bp) was tested in all mobility shift assays, as a positive control [26]. 2C. Protect40 (40 bp) fragments with an altered sequence were used: Pho1 (DR1 altered sequence), Pho2 (DR2 altered sequence) and Pho3 (DR1 and DR2 altered sequences). Wild-type Protect40 (40 bp) was tested as a positive control [26].
doi:10.1371/journal.pone.0042876.g002

## DR sequence analyses

Using two different approaches (see Materials and Methods section), we identified three motifs matching the original DR1, DR2 and DR3 repeats of the PhoP promoter. These repeats are also present in all target regions upstream the *lipF*, *pks2*, *msl3* and *fadD21* genes (Fig. 5). The degree of sequence identity to the consensus of the canonical DR1–3 motifs of the *phoP* promoter region [23] ranged from 56% to 89% (Table 3). DR1 was the most conserved motif, followed by DR2, and then DR3. The DR1 and DR2 motifs were separated by 1 to 5 bp. By contrast, DR2 and DR3 were separated by 13 to 35 bp (Fig. 5). Taking this into account, we analyzed the distribution of the consensus sequences for each motif throughout the MTB genome. Motif configurations with only DR1 and DR2, or with DR1, DR2 and DR3 were observed in the upstream regions of 87 genes, accounting for 83.9% and 16.1% of total cases, respectively. The identified genes were analyzed according to the classes of the Clusters of Orthologous Groups (COGs) classification [29,30] (Fig. 6A, Table S1). We further investigated, using the Tuberculist data, the function of genes that are not listed in COGs (38 genes, Fig. 6B). Some of the identified genes have already been reported to be regulated by PhoP on the basis of transcriptomic data [20,21]. Based on the COG classification, we found that PhoP play a key role as a putative transcription regulator mostly for genes involved in transcription (21%) and Secondary metabolites biosynthesis, transport and catabolism (14%), lipid transport and metabolism (9%), and coenzyme transport and metabolism (9%) (Fig. 6A).

## Discussion

All living organisms have developed regulatory networks allowing them to survive in different environments, through rapid adaptation to hostile conditions and the conservation of energy by blocking pointless biochemical biosynthesis. TCS regulators, which have an environmentally responsive component, play a

major role in this process. One of the advantages of this mode of control is that it allows the expression of complex traits or actions, such as membrane biosynthesis, host defences escape or adaptation, to be coordinated. This regulation results in changes to the patterns of expression of multiple genes [33,34].

The MTB PhoP/PhoR TCS regulates the expression of more than 100 genes, based on differential gene expression data [20]. This includes genes encoding proteins involved in polyketide synthesis, for which differences in expression have been observed between PhoP null mutants and wild-type strains.

We aimed to identify genes directly regulated by PhoP. Transcriptome analyses have identified groups of genes regulated by PhoP, including genes encoding regulators that might function as intermediates in regulatory cascades [20,21]. An analysis and comparison of these groups showed that they had only a few genes in common. This result reflects the problems involved in dealing with genes displaying low-level regulation and the importance of taking into account different experimental conditions and different strains genotypes. Indeed, as PhoP plays a key role in the regulation of several virulence factors, its activation varies with stressful conditions and the nature of the stimulus of the cognate sensor. We used a combination of genetic, biochemical and bioinformatics approaches to define PhoP binding sites as hallmarks of PhoP-regulated genes. We studied *lipF*, *msl3*, *pks2* and *fadD21*, which were shown to be upregulated by PhoP in genetic [29] and transcriptomic studies [20,21]. However, there has been little evidence for the mechanism of regulation by PhoP in terms of DNA recognition sites, binding to selected genes or search for genes distributed throughout the genome with sub threshold regulation levels.

*M. smegmatis* has been used as a model strain to test *lipF*, *pks2*, *msl3*, and *fadD21* lacZ gene fusions for the identification of promoters. The flanking regions of these genes were used for direct assessments of PhoP interaction by *in vitro* EMSA and DNase I protection assays. All experiments have been done using

<table><caption>**Table 2.** List of DNA fragments used in this study.</caption><thead><tr><th>Oligonucleotides</th><th>Sequence (5'-3' position relative to the translation start)</th></tr></thead><tbody><tr><td>lipFa1<sup>b</sup></td><td>(-575) AGACGTACAGCAAACTCCCAGTCATACACCGATGCATTTCAAGGTGCCTGTGATCCCCTAAATGCAGG (-508)</td></tr><tr><td>pks2a<sup>b</sup></td><td>(-52) GACACAGCTACATCGAAGGATTGCTAGTTAGGCTACGTATCGCCTCCGGCATTGGGCTT (+6)</td></tr><tr><td>pks2b<sup>b</sup></td><td>(-95) GTGCCCCAGTAGCACAGCCGCTTAGAACTAAAGAGCCACC (-56)</td></tr><tr><td>msl3a1<sup>b</sup></td><td>(-242) CGTCTGGTAGCGGCATGGCAACGGCCTGTGAGTTGGCTGGATAAT (-198)</td></tr><tr><td>fadD21a<sup>b</sup></td><td>(-152) GGCGCCCCTGCGTACACCGACTTCGACTCTGCGTCAACCTGTTTCAGCACATGCA (-98)</td></tr><tr><td>Protect40<sup>b</sup></td><td>(-70) GACTGTTAGCAGACTACTGGCAACGAGCTTTCAGGAATTA (-31)</td></tr><tr><td>Pho1<sup>b</sup></td><td>(-70) GCTCAGGTAAAGACTACTGGCAACGAGCTTTCAGGAATTA (-31)</td></tr><tr><td>Pho2<sup>b</sup></td><td>(-70) GACTGTTAGCAGACTTGAACGTCGGAGCTTTCAGGAATTA (-31)</td></tr><tr><td>Pho3<sup>b</sup></td><td>(-70) GCTCAGGTAAAGACTTGAACGTCGGAGCTTTCAGGAATTA (-31)</td></tr></tbody></table>

<footnote><sup>b</sup>Fragments used for electrophoretic mobility shift assays.
doi:10.1371/journal.pone.0042876.t002</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Schematic diagram of the DNA sequences tested in electrophoretic mobility assays. The DR1, DR2 and DR3 sites are indicated in blue, orange and red, respectively. Binding is indicated by (+), with (−) indicating an absence of binding and (−/+) indicating weak binding. Positions on the DNA regions are shown with respect to the translation start site, in brackets.
doi:10.1371/journal.pone.0042876.g003</figcaption></figure>

phosphorylated PhoP protein as unphosphorylated MTB PhoP binds weakly to oligonucleotides derived from the *lipF*, *pks2*, *msl3*, and *fadD21*, promoter sequences as previously described [31,32,35]. These observations are similar to those for PhoB from *E. coli*, in which PhoR and PhoB constitute a two-component system controlling phosphate uptake [35].

The first study on PhoP protein interaction with its own promoter revealed the presence of three direct repeat motifs (DRs) [26]. The first two direct repeats have been shown later to be sufficient for this interaction [30,32]. More recently, the analyses of the promoter regions of *pks2* and *msl3* genes, revealed the presence of two DRs motifs that are needed for recognition by PhoP protein [31]. Nevertheless the motifs published by Goyal *et al.* (2011) are different from those that we identified. While the previously identified DRs regions are different from those reported in our study, the two works identified, overall, the same protected region. The main difference resides in the fact that we took in

consideration the previously confirmed DRs regions identified in the promoter region of the *phoP* gene. This led us to clearly identify the presence of DR1–3 in some of the promoters that bound PhoP-P, and also allowed us to witness in some cases loss of binding when DR3 was missing. This is confirmed by our footprinting assays and gel shift experiments.

In this study, we show that for all studied genes, we were able to identify in their promoter regions three DRs namely DR1, DR2 and DR3. In addition, the DR3 repeat is essential for the binding of PhoP to *pks2* and *fadD21* DR1–2 repeats. This result is consistent with a potential cooperative interaction between DR3 and DR1–2. By comparing all identified DRs motifs we could suggest that DR3 assisted the binding to all PhoP-P-regulated sites showing low levels of similarity to *phoP*'s DR1 and DR2 sequences (Table 3). Moreover, the ratio of PhoP to DNA required to achieve this binding is variable, thus resulting in additional fine-tuning of expression for different genes under the control of PhoP.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** DNase I footprinting assay of the *pks2* and *msl3* promoter regions with PhoP-P. Radiolabeled PCR fragments corresponding to (A) *pks2* (−136 to +62), (B) *msl3* (−312 to −79), (C) *fadD21* (−194 to +7) and (D) *lipF* (−639 to −450) were used as DNA targets. Various amounts of PhoP-P were used with *pks2* and *msl3* (lane 1: 0; lane 2: 1.8 pmol; lane 3: 3.6 pmol; lane 4: 7.2 pmol and lane 5: 14.4 pmol) and were incubated with 0.2 pmol of DNA before DNase I digestion. For *fadD21* and *lipF* different amount of PhoP-P were used (lane 1: 0; lane 2: 1.8 pmol; lane 3: 18 pmol; lane 4: 180 pmol and lane 5: 900 pmol) and were incubated with 1 pmol of DNA before DNase I digestion. Lane 6: A+G Maxam and Gilbert reaction. Protected regions are indicated by a line with colored regions. The blue, orange and red segments correspond to the DR1, DR2 and DR3 sites, respectively.</figcaption></figure>

doi:10.1371/journal.pone.0042876.g004

The *lipF* promoter region contains, at positions −574, −563, and −530, sequences similar to the previously described DR1, DR2 and DR3 of PhoP. Richter et al. (2007) reported that the position of the transcription start was 511 bp upstream from the translation initiation site. Our results indicate that the *lipF* −10 site identified by Richter et al. (2007) overlaps with the newly assigned DR3 box described here. Based on the distance between the *lipF* DR1/DR2 sites and the transcription start site, we suggest that a tandem head-to-tail PhoP dimer, close to the −35, binds to the DR1/DR2 repeats that occupy adjacent same-DNA-face strands. The same configuration is also observed for *pks2*. A similar situation has been described for *E. coli* PhoB, which activates transcription by interacting with the sigma 70 subunit of the RNA polymerase in promoters in which the −35 sigma recognition element is replaced by the Pho box.

PhoP dimerization would bring the two copies of the DNA-binding domains close to each other, facilitating binding to DR1 and DR2, which are separated by only two nucleotides in this case. The presence of DR1/DR2 PhoP binding sites that overlaps with the −35 transcription region of *lipF* and *pks2* might be indicative of

a class II activation mechanism according to which the regulator PhoP binds the region 4 of the sigma factor [36]. This configuration might result in the recruitment of the RNA polymerase to the promoter and subsequent induction of transcription start. The presence of DR1/DR2 PhoP binding site upstream the −35 transcription region for *msl3* might involve another mechanism of transcription activation that remains to be investigated.

With the use of various bioinformatics analyses, we fine-tuned the three consensus sequences for the individual DR1, DR2 and DR3 sites (Fig. 6). Only DR1 and DR2 were highly conserved in all analysed regions. We used these confirmed consensus sequences by taking into account the length of the spacer sequences between them to investigate the distribution of these motifs throughout the *Mycobacterium tuberculosis* genome, in order to identify all genes that could be potentially regulated directly by the PhoP/PhoR TCS. We found these motifs in the upstream regions of 87 genes, including genes having already been reported to be regulated by PhoP [21]. The classification of all identified genes into COG categories revealed a predominance of genes encoding proteins

<table><caption>**Table 3.** Identification of the DR1–3 sites on the DNA promoters recognized by PhoP.</caption><thead><tr><th></th><th>DR1</th><th>DR2</th><th>DR3</th></tr></thead><tbody><tr><td>msl3</td><td>TCTGGTAGC 8/9: 89%</td><td>CATGGCAAC 7/9: 78%</td><td>AATGTGTTTC 5/9: 56%</td></tr><tr><td>fadD21</td><td>TGTTTCAGC 7/9: 78%</td><td>ATGCACAGC 6/9: 67%</td><td>ATGATCAGC 7/9: 78%</td></tr><tr><td>lipF</td><td>ACGTACAGC 8/9: 89%</td><td>ACTCCCAGT 6/9: 67%</td><td>CCTGTGATC 6/9: 67%</td></tr><tr><td>pks2</td><td>CCCAGTAGC 6/9: 67%</td><td>CCGCTTAGA 6/9: 67%</td><td>AGACACAGC 5/9: 56%</td></tr><tr><td>phoP</td><td>ACTGTTAGC 9/9: 100%</td><td>ACTGGCAAC 9/9: 100%</td><td>-</td></tr><tr><td>consensus sequence of phoP</td><td colspan="3">A C T/G T/G T/G Py A Pu C</td></tr></tbody></table>

For each DNA sequence shown to bind PhoP-P, sequences with a complete or partial match to the DR1, DR2 and DR3 sites were identified (underlined). The degree of conservation of the consensus sequence was calculated and is shown as the ratio of matching nucleotides to the consensus: % identity.

doi:10.1371/journal.pone.0042876.t003
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Identification of the DR1, DR2 and DR3 repeats in the upstream regions of the *phoP*, *lipF*, *pks2*, *msl3* and *fadD21* genes. (A) Nucleotide sequences and mapping of DR1, DR2 and DR3 sites. (B) Illustration of the consensus sequences of DR1, DR2 and DR3 obtained with MEME and Weblogo. The +1 indicate the transcription start site according to Richter *et al.*, (2007) for *lipF* and Goyal *et al.*, (2011) for *pks2* and *msl3*. doi:10.1371/journal.pone.0042876.g005</figcaption></figure>

involved in transcriptional regulation (Fig. 6). This result reflects and confirms the role of *phoP* as a key player in the virulence of tubercle bacilli by controlling the activity of several genes. Indeed, PhoP has been shown to regulate the differential transcription of up to 600 genes in *E. coli* [37]. Genes encoding proteins involved in the biosynthesis of polyketide-derived lipids were also well represented, providing additional support for our approach. These genes have already been reported to be under the control of the phoP/phoR regulon [20,21]. Some of the other genes identified were classified as involved in lipid biosynthesis pathways. This finding is consistent with previous studies showing that synthesis of the lipid components of the cell surface is regulated by PhoB in *E. coli*.

Our analysis provides information that cannot be obtained easily by transcriptomic approaches, because the stimuli of the PhoP cognate sensor have yet to be identified and because transcriptional regulators are often only weakly expressed. The identification of genes directly regulated by PhoP paves the way for the characterization of genes involved in pathogenicity. The role of these genes should be further investigated by genetic approaches and may be used as a starting point for the design of new drugs and vaccines.

## Materials and Methods

### Bacterial strains, media and growth conditions

*Mycobacterium tuberculosis* H37Rv [6], *M. smegmatis* mc²¹⁵⁵ [38] and their derivatives were cultured at 37°C in Middlebrook 7H9 medium supplemented with ADC (0.5% bovine serum albumin, 0.2% dextrose, 0.085% NaCl, 0.0003% beef catalase; Difco) and 0.05% Tween 80 or on solid Middlebrook 7H10 medium

supplemented with OADC (0.05% oleic acid, 0.5% bovine serum albumin, 0.2% dextrose, 0.085% NaCl, 0.0003% beef catalase; Difco). *Escherichia coli* cultures for cloning were grown in Luria-Bertani (LB) broth or on LB agar plates at 37°C. When required, antibiotics were added to the medium at the following concentrations: kanamycin (20 µg/ml), ampicillin (100 µg/ml). We added isopropyl-β-D-thiogalactopyranoside (IPTG) at a concentration of 0.5 mM and X-Gal at a concentration of 50 µg/ml, when required.

### Plasmid construction

For construction of the *lipF*, *pks2*, *msl3* and *fadD21* promoter-lacZ fusions, we amplified the corresponding fragments from *M. tuberculosis* H37Rv genomic DNA by PCR with the primers shown in Table 1. The PCR products were digested with *Bam*HI and *Kpn*I and inserted into the corresponding sites in the *pJEM15* *E. coli*-mycobacterial shuttle plasmid [39]. We used *pJEM31-PAN* containing a promoter sequence, PAN, isolated from *M. paratuberculosis* as a positive control. This promoter lies adjacent to and does not overlap the 3' end of an IS900 mobile genetic element [28]. This plasmid was used to transform *M. smegmatis* mc²¹⁵⁵. The strains were grown at 37°C in 7H9 medium supplemented with kanamycin (20 µg/ml), until they reached an optical density at 600 nm (OD₆₀₀ nm) of 0.6–0.8. The cells were collected by centrifugation and β-galactosidase activity in the cell extract was evaluated to determine the levels of transcription from the various promoter regions studied (*lipF*, *pks2*, *msl3*, and *fadD21*).

### Quantification of beta-galactosidase activity

Mycobacterial strains were grown to the exponential growth phase (OD₆₀₀ nm = 0.8). The cultures were centrifuged, and the
A

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

B

<figure><img src="image_11.png" /><figcaption>**Figure 6.** Distribution of the identified consensus sequences of DR1, DR2 and DR3 in the upstream regions of all *Mycobacterium tuberculosis* H37Rv genes. (A) Distribution according to the Clusters of Orthologous Groups (COG) classification. (B) Genes not listed in the COGs classification were classified on the basis of the functions listed in Tuberculist.
doi:10.1371/journal.pone.0042876.g006</figcaption></figure>

pellets were washed with phosphate-buffered saline. The cells were resuspended in 1 ml Z buffer (0.06 M Na₂HPO₄, 0.04 M NaH₂PO₄, 0.01 M KCl, 1 mM MgSO₄ and 0.05 M β-mercaptoethanol pH 7 at 25°C and sonicated to generate a cell extract. The chromogenic substrate o-nitrophenyl-D-galactoside (Sigma, France) was added to the cell extracts at a final concentration of 0.66 mg/ml. The mixtures were incubated at 28°C for 1 hour and the enzymatic reactions were stopped by adding 0.29 M Na₂CO₃. The A₄₂₀ nm of the supernatant was determined, and β-galactosidase activity was calculated in Miller units, according to the following formula: β-galactosidase activity = (1000×A₄₂₀ nm)/(time (min)×aliquot volume (ml)×OD₆₀₀ nm).

## Electrotransformation of mycobacteria

50 ml culture of bacteria was grown to an OD₆₀₀ nm of 0.6 to 0.8. The cells were collected by centrifugation, washed twice in 10% glycerol and resuspended in 2 ml of 10% glycerol. Aliquots (400 µl) were electroporated with vector DNA in 0.2 cm-path length cuvettes (Bio-Rad), with a single pulse (2.5 kV, 25 µF, 1000 Ω). Cells were transferred to 1 ml of 7H9-ADC-0.05% Tween 80 and incubated for 2 hours at 37°C.

## Production, purification and phosphorylation of PhoP

The PhoP protein was produced in *E. coli* Bl21(DE3) pLysS (Agilent Technologies), from the pET15b expression plasmid (Novagen, Merck Chemicals France) [30]. PhoP overproduction was induced by adding 1 mM IPTG (Euromedex, France) to LB

medium and incubating for 3 hours at 30°C. The cells were collected by centrifugation at 5,500 rpm (SLA-4000 rotor, RC5c (Sorvall) centrifuge) for 30 min and resuspended in 50 mM Tris-HCl pH 8.0, 1 M NaCl, 5 mM KCl, 10% glycerol, 5 mM imidazole, 100 µM phenyl methyl sulfonyl fluoride (PMSF) (Buffer A). The cells were lysed by passage through an Emulsiflex-C5 (Avestin Europe) at 4°C. PhoP was purified on a 5-ml HiTrap Chelating HP (GE Healthcare) affinity column loaded with 0.5 M Ni₂SO₄ and equilibrated with Buffer A, on an AKTA-prime FPLC machine (GE Healthcare). The histidine tag was cleaved overnight at 4°C, with a Thrombin CleanCleavage kit (Sigma Aldrich). The released PhoP was further purified with a second passage through the 5-ml HiTrap Chelating HP Ni-affinity column. Pure protein (>95%, based on SDS-PAGE) was obtained after several rounds of FPLC, and stored at -80°C in 50 mM Tris-HCl pH 8, 150 mM NaCl, 5 mM KCl, and 20% glycerol.

Aliquots of 200 nmoles of PhoP were incubated for 90 minutes at 30°C with 50 mM acetyl phosphate (AcP) (Sigma Aldrich) in 100 mM Tris-HCl pH 7.0, 10 mM MgCl₂, 150 mM KCl, to obtain phosphorylated PhoP (PhoP-P).

## DNA fragment synthesis and purification

Double-stranded DNA fragments (ds-DNAs) were labeled at the 5' end by PCR with commercially synthesized 5'-labeled IRD700 or DY682 (Eurofins MWG Operon, Germany) primer and the various templates described in Table 1, 2.
ds-DNAs were purified by electrophoresis in 8% polyacrylamide gels, electroeluted in 0.5× TBE (Biosolve, Netherlands), through 10 kDa MWCO (molecular weight cut off) dialysis tubing (Spectra/Por, Spectrum Labs, USA), and dried on a Speed Vac. The DNAs were resuspended in 60 mM sodium acetate pH 5.4, precipitated in 100% ethanol, washed in 70% ethanol and recovered in 10 mM Tris pH 7.8, 0.1 mM EDTA (Gibco). The DNA concentrations were deduced from A<sub>260 nm</sub> before use or storage.

## Electrophoretic mobility shift assay (EMSA)

Aliquots of 40 nM labeled ds-DNAs were separately mixed with various concentrations of PhoP or PhoP-P, in a final volume of 30 µl in the presence of poly dI-dC at 10 µg/ml. The mixture was dialyzed in several steps against 50 mM Tris-HCl, pH 8.0, 1 mM dithiothreitol (DTT), 1 mM EDTA, 100 µM PMSF, with decreasing concentrations of NaCl (750 mM, 250 mM, 150 mM respectively), at 4°C. The DNA/PhoP(-P) complex obtained at the end of dialysis was mixed with 3 µl of gel loading buffer (0.1% bromophenol blue and 40% sucrose), loaded onto a 20-cm×20-cm 8% native polyacrylamide gel and subjected to electrophoresis at 180 V/cm for 3 hours at 18°C, with 0.5× TBE used as the running buffer. The gels were then transferred onto blotting paper (Whatman 3MM CHR) and scanned on an Odyssey® Infrared Imaging System at 700 nm (Li-cor Inc, NE, USA).

## DNase I footprinting experiments

DNA probes were prepared by amplifying fragments of the *lipF*, *pks2*, *msl3* and *fadD21* upstream regions by PCR with the lipFDF/lipFDR, *msl3DF/msl3DR* *pks2DF/pks2DR* and *fadD21DF/fadD21DR* primers, respectively, using *M. tuberculosis* chromosome DNA as the template. In each case, the 5' end of the forward primer (*lipFDF*, *msl3DF*, *pks2DF* and *fadD21DF*) was labeled with [γ-³²P]ATP, with T4 polynucleotide kinase. Before the DNA binding reaction, PhoP-P was obtained as described above and dialyzed against a buffer containing 20 mM Tris pH 8, 1 mM EDTA, 200 mM NaCl, 50% glycerol. DNase I footprinting reactions were performed as previously described [25]. Briefly, 0 to 15 pmol of PhoP-P was mixed with 0.2 pmol of DNA and incubated with DNase I at room temperature (~24°C) for 1 min. The samples were analyzed by electrophoresis on a 6% polyacrylamide gel containing 7 M urea. Maxam and Gilbert sequencing ladders (G+A) were also loaded on the same gel.

## Sequence analysis

We narrowed down the DNA sequences potentially recognized by PhoP-P (to ~60 bp) by EMSA, with iterative sequence alignment to DR1 and DR2 (EMBOSS package, ClustalW).

## References

1. Bermudez LE, Goodman J (1996) *Mycobacterium tuberculosis* invades and replicates within type II alveolar cells. *Infect Immun* 64: 1400-1406.

2. Bodnar KA, Serbina NV, Flynn JL (2001) Fate of *Mycobacterium tuberculosis* within murine dendritic cells. *Infect Immun* 69: 800-809.

3. Hickman SP, Chan J, Salgame P (2002) *Mycobacterium tuberculosis* induces differential cytokine production from dendritic cells and macrophages with divergent effects on naive T cell polarization. *J Immunol* 168: 4636-4642.

4. Jiao X, Lo-Man R, Guermonprez P, Fiette L, Deriaud E, et al. (2002) Dendritic cells are host cells for mycobacteria in vivo that trigger innate and acquired immunity. *J Immunol* 168: 1294-1301.

5. Tailleux L, Waddell SJ, Pelizzola M, Mortellaro A, Withers M, et al. (2008) Probing host pathogen cross-talk by transcriptional profiling of both *Mycobacterium tuberculosis* and infected human dendritic cells and macrophages. *PLoS One* 3: e1403.

6. Cole ST, Barrell BG (1998) Analysis of the genome of *Mycobacterium tuberculosis* H37Rv. *Novartis Found Symp* 217: 160-172; discussion 172-167.

The results were confirmed by footprinting. The resulting sequence datasets were then analyzed further, by two different approaches. In the first, we used a Python script (this study) to search the identified regions for motifs identical and similar (no more than three nucleotides of difference) to those already identified in the upstream region of the PhoP gene [20,26,32]. We then analyzed these regions with the MEME (Multiple Em for Motif Elicitation) program [40], which uses the expectation-maximization (EM) algorithm for iterative improvement of a model of the motif. Finally, we constructed consensus motif sequences with WebLogo 3 [41,42]. The distribution of these consensus sequences in the *Mycobacterium tuberculosis* genome was analyzed with a Python script (this study) by scanning the 800 bp regions immediately upstream all *M. tuberculosis* H37Rv genes.

## Supporting Information

**Figure S1 Electrophoretic mobility assays with the lipFb, lipFc, msl3b, lipFa2 and pks2c fragments and PhoP-P.** lipFb (227 bp), lipFc (182 bp), msl3b (98 bp), lipFa2 (89 bp) and pks2c (77 bp) fragments were incubated with PhoP-P in the presence of poly dI-dC at 10 µg/ml and run on a native polyacrylamide gel (8%), in 0.5× TBE buffer. Each fragment was incubated in the presence of 8 µM PhoP-P, with the exception of msl3b, for which 10 µM PhoP-P. (EPS)

**Figure S2 Electrophoretic mobility assays for the lipFa, pks2, msl3 and fadD21 fragments with phosphorylated (P) and unphosphorylated (NP) PhoP.** The large DNA fragments initially selected — lipFa (222 bp), pks2 (148 bp), msl3a (235 bp), and fadD21 (210 bp) — were incubated with PhoP-P and unphosphorylated PhoP in the presence of poly dI-dC at 10 µg/ml and run on a native polyacrylamide gel (8%), in 0.5× TBE buffer. Each fragment was incubated in the presence of 4 µM PhoP-P and PhoP, with the exception of lipFa, for which 8 µM PhoP-P and PhoP was required. (TIFF)

**Table S1 List of genes with DR1, DR2 and DR3 in their upstream regions according to Figure 5.** First Excel spreadsheet: DR1+DR2; Second spreadsheet: DR1+DR2+DR3. (XLS)

## Author Contributions

Conceived and designed the experiments: MC BG DNG. Performed the experiments: MC CT SD. Analyzed the data: MC CT AN SD BG DNG. Contributed reagents/materials/analysis tools: MC CT AN SD BG DNG. Wrote the paper: MC CT AN SD BG DNG.

7. Perez E, Samper S, Bordas Y, Guilhot C, Gicquel B, et al. (2001) An essential role for phoP in *Mycobacterium tuberculosis* virulence. *Mol Microbiol* 41: 179-187.

8. Martin C, Williams A, Hernandez-Pando R, Cardona PJ, Gormley E, et al. (2006) The live *Mycobacterium tuberculosis* phoP mutant strain is more attenuated than BCG and confers protective immunity against tuberculosis in mice and guinea pigs. *Vaccine* 24: 3408-3419.

9. Frigui W, Bottai D, Majlessi L, Monot M, Josselin E, et al. (2008) Control of *M. tuberculosis* ESAT-6 secretion and specific T cell recognition by PhoP. *PLoS Pathog* 4: e33

10. Lee JS, Krause R, Schreiber J, Mollenkopf HJ, Kowall J, et al. (2008) Mutation in the transcriptional regulator PhoP contributes to virulence of *Mycobacterium tuberculosis* H37Ra strain. *Cell Host Microbe* 3: 97-103.

11. Zheng H, Lu L, Wang B, Pu S, Zhang X, et al. (2008) Genetic basis of virulence attenuation revealed by comparative genomic analysis of *Mycobacterium tuberculosis* strain H37Ra versus H37Rv. *PLoS One* 3: e2375.
12. Wang S, Engohang-Ndong J, Smith I (2007) Structure of the DNA-binding domain of the response regulator PhoP from *Mycobacterium tuberculosis*. *Biochemistry* 46: 14751-14761.

13. Birck C, Chen Y, Hulett FM, Samama JP (2003) The crystal structure of the phosphorylation domain in PhoP reveals a functional tandem association mediated by an asymmetric interface. *J Bacteriol* 185: 254-261.

14. Groisman EA, Cromie MJ, Shi Y, Latifi T (2006) A Mg2+-responding RNA that controls the expression of a Mg2+ transporter. Cold Spring Harb Symp Quant Biol 71: 251-258.

15. Alpuche Aranda CM SJ, Loomis WP, Miller SI (1992) *Salmonella typhimurium* activates virulence gene transcription within acidified macrophage phagosomes. Proc Natl Acad Sci U S A Nov 1; 89 (21): 10079-10083.

16. Martin-Orozco N, Touret N, Zaharik, ML Kopelman, R Miller, S, et al. (2006) Visualization of vacuolar acidification-induced transcription of genes of pathogens inside macrophages. Mol Biol Cell Jan; 17 (1): 498-510.

17. Bader MW, Navarre WW, Shiaw W, Nikaido H, Frye JG, et al. (2003) Regulation of *Salmonella typhimurium* virulence gene expression by cationic antimicrobial peptides. Mol Microbiol 50: 219-230.

18. Bader MW, Sanowar S, Daley ME, Schneider AR, Cho U, et al. (2005) Recognition of antimicrobial peptides by a bacterial sensor kinase. Cell 122: 461-472.

19. Miller SI, Pulkkinen WS, Selsted ME, Mekalanos JJ (1990) Characterization of defensin resistance phenotypes associated with mutations in the phoP virulence regulon of *Salmonella typhimurium*. *Infect Immun* 58: 3706-3710.

20. Gonzalo-Asensio J, Mostowy S, Harders-Westerveen J, Huygen K, Hernandez-Pando R, et al. (2008) PhoP: a missing piece in the intricate puzzle of *Mycobacterium tuberculosis* virulence. PLoS One 3: e3496.

21. Walters SB, Dubnau E, Kolesnikova I, Laval F, Daffe M, et al. (2006) The *Mycobacterium tuberculosis* PhoPR two-component system regulates genes essential for virulence and complex lipid biosynthesis. Mol Microbiol 60: 312-330.

22. Gonzalo Asensio J, Maia C, Ferrer NL, Barilone N, Laval F, et al. (2006) The virulence-associated two-component PhoP-PhoR system controls the biosynthesis of polyketide-derived lipids in *Mycobacterium tuberculosis*. *J Biol Chem* 281: 1313-1316.

23. Sirakova TD, Thirumala AK, Dubey VS, Sprecher H, Kolattukudy PE (2001) The *Mycobacterium tuberculosis* pks2 gene encodes the synthase for the hepta- and octamethyl-branched fatty acids required for sulfolipid synthesis. *J Biol Chem* 276: 16833-16839.

24. Dubey VS, Sirakova TD, Kolattukudy PE (2002) Disruption of *msl3* abolishes the synthesis of mycolipanoic and mycolipenic acids required for polyacyltylrehalose synthesis in *Mycobacterium tuberculosis* H37Rv and causes cell aggregation. Mol Microbiol 45: 1451-1459.

25. Camacho LR, Ensergueix D, Perez E, Gicquel B, Guilhot C (1999) Identification of a virulence gene cluster of *Mycobacterium tuberculosis* by signature-tagged transposon mutagenesis. Mol Microbiol 34: 257-267.

26. Gupta S, Sinha A, Sarkar D (2006) Transcriptional autoregulation by *Mycobacterium tuberculosis* PhoP involves recognition of novel direct repeat sequences in the regulatory region of the promoter. FEBS Lett 580: 5328-5338.

27. Saviola B, Woolwine SC, Bishai WR (2003) Isolation of acid-inducible genes of *Mycobacterium tuberculosis* with the use of recombinase-based in vivo expression technology. *Infect Immun* 71: 1379-1388.

28. Murray A, Winter N, Lagranderie M, Hill DF, Rauzier J, et al. (1992) Expression of Escherichia coli beta-galactosidase in *Mycobacterium bovis* BCG using an expression system isolated from *Mycobacterium paratuberculosis* which induced humoral and cellular immune responses. Mol Microbiol 6: 3331-3342.

29. Chesne-Seck ML, Barilone N, Boudou F, Gonzalo Asensio J, Kolattukudy PE, et al. (2008) A point mutation in the two-component regulator PhoP-PhoR accounts for the absence of polyketide-derived acyltrehaloses but not that of phthiocerol dimycocerosates in *Mycobacterium tuberculosis* H37Ra. *J Bacteriol* 190: 1329-1334.

30. Gonzalo-Asensio J, Soto CY, Arbues A, Sancho J, del Carmen Menendez M, et al. (2008) The *Mycobacterium tuberculosis* phoPR operon is positively autoregulated in the virulent strain H37Rv. *J Bacteriol* 190: 7068-7078.

31. Goyal R, Das AK, Singh R, Singh PK, Korpole S, et al. (2011) Phosphorylation of PhoP protein plays direct regulatory role in lipid biosynthesis of *Mycobacterium tuberculosis*. *J Biol Chem* 286: 45197-45208.

32. Gupta S, Pathak A, Sinha A, Sarkar D (2009) *Mycobacterium tuberculosis* PhoP recognizes two adjacent direct-repeat sequences to form head-to-head dimers. *J Bacteriol* 191: 7466-7476.

33. Chang C, Stewart RC (1998) The two-component system. Regulation of diverse signaling pathways in prokaryotes and eukaryotes. Plant Physiol 117: 723-731.

34. Parkinson JS, Kofoid EC (1992) Communication modules in bacterial signaling proteins. Annu Rev Genet 26: 71-112.

35. Makino K, Amemura M, Kawamoto T, Kimura S, Shinagawa H, et al. (1996) DNA binding of PhoB and its interaction with RNA polymerase. *J Mol Biol* 259: 15-26.

36. Dove SL, Darst SA, Hochschild A (2003) Region 4 of sigma as a target for transcription regulation. Mol Microbiol 48: 863-874.

37. Alteri CJ, Lindner JR, Reiss DJ, Smith SN, Mobley HL (2011) The broadly conserved regulator PhoP links pathogen virulence and membrane potential in *Escherichia coli*. *Mol Microbiol* 82: 145-163.

38. Snapper SB, Melton RE, Mustafa S, Kieser T, Jacobs WR, Jr. (1990) Isolation and characterization of efficient plasmid transformation mutants of *Mycobacterium smegmatis*. Mol Microbiol 4: 1911-1919.

39. Timm J, Lim EM, Gicquel B (1994) *Escherichia coli*-mycobacteria shuttle vectors for operon and gene fusions to lacZ: the pJEM series. *J Bacteriol* 176: 6749-6753.

40. Bailey TL, Boden M, Buske FA, Frith M, Grant CE, et al. (2009) MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res 37: W202-208.

41. Crooks GE, Hon G, Chandonia JM, Brenner SE (2004) WebLogo: a sequence logo generator. Genome Res 14: 1188-1190.

42. Schneider TD, Stephens RM (1990) Sequence logos: a new way to display consensus sequences. Nucleic Acids Res 18: 6097-6100.