Complexity of bacterial phosphorylation interaction network
- 1Institut National de la Recherche Agronomique, UMR1319 Micalis, Jouy-en-Josas, France
- 2Systems and Synthetic Biology, Department of Chemical and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
Signal transduction in eukaryotes is generally transmitted through phosphorylation cascades that involve a complex interplay of transmembrane receptors, protein kinases, phosphatases and their targets. Our previous work indicated that bacterial protein-tyrosine kinases and phosphatases may exhibit similar properties, since they act on many different substrates. To capture the complexity of this phosphorylation-based network, we performed a comprehensive interactome study focused on the protein-tyrosine kinases and phosphatases in the model bacterium Bacillus subtilis. The resulting network identified many potential new substrates of kinases and phosphatases, some of which were experimentally validated. Our study highlighted the role of tyrosine and serine/threonine kinases and phosphatases in DNA metabolism, transcriptional control and cell division. This interaction network reveals significant crosstalk among different classes of kinases. We found that tyrosine kinases can bind to several modulators, transmembrane or cytosolic, consistent with a branching of signaling pathways. Most particularly, we found that the division site regulator MinD can form a complex with the tyrosine kinase PtkA and modulate its activity in vitro. In vivo, it acts as a scaffold protein which anchors the kinase at the cell pole. This network highlighted a role of tyrosine phosphorylation in the spatial regulation of the Z-ring during cytokinesis.
Introduction
In eukarya, protein phosphorylation on serine/threonine and tyrosine residues is catalyzed by the Hanks-type family of protein kinases (Hanks et al., 1988). In bacteria, Hanks-type kinases catalyze phosphorylation of proteins only on serine and threonine residues (Pereira et al., 2011), and a distinct family of kinases, termed bacterial tyrosine kinases (BY-kinases), phosphorylate proteins on tyrosine (Grangeasse et al., 2007). BY-kinases have been initially described as autophosphorylating enzymes involved in exopolysaccharide synthesis and export (Whitfield, 2006). More recently, it became apparent that they phosphorylate many protein substrates and regulate their activity (Shi et al., 2010). Arguably the best characterized bacterial system with respect to physiological substrates of BY-kinases is the Firmicute model organism Bacillus subtilis (Mijakovic et al., 2005b). B. subtilis possesses a well characterized BY-kinase PtkA (Mijakovic et al., 2003) and a putative BY-kinase PtkB (EpsB). In order to phosphorylate its substrates, PtkA requires a transmembrane activator TkmA (Mijakovic et al., 2003), which is encoded by the same operon as the kinase. PtkB-encoding operon also codes for a transmembrane activator termed TkmB (Mijakovic et al., 2005b). While structurally the TkmB/PtkB pair constitutes a bona fide BY-kinase, its activity toward substrates has never been experimentally demonstrated. Recent study suggested that TkmB interaction with PtkB could modulate its function in biofilm formation in B. subtilis (Gerwig et al., 2014). Homologs of TkmA and TkmB exist also in proteobacteria, but there they are “fused” with the BY-kinases in a single polypeptide chain, encoded by a single gene (Jadeau et al., 2012). The rationale behind these different architectures has puzzled researchers for years. One hypothesis suggests that dissociation of BY-kinases from transmembrane activators in Firmicutes may enable these kinases to interact with alternative (possibly cytosolic) modulator proteins (Shi et al., 2010).
B. subtilis PtkA is known to phosphorylate a broad spectrum of substrates, including UDP-glucose dehydrogenases (Mijakovic et al., 2003), single-stranded DNA-binding proteins (Mijakovic et al., 2006), transcription factors (Derouiche et al., 2013) and other enzymes (Jers et al., 2010). Cells devoid of PtkA exhibit a pleiotropic phenotype with defects in cell cycle (Petranovic et al., 2007) and biofilm formation (Kiley and Stanley-Wall, 2010). In addition to regulating the enzyme activity of certain substrates, PtkA also seems to affect sub-cellular localization of others (Jers et al., 2010). Bacterial Hanks-type serine/threonine (Ser/Thr) kinases are equally promiscuous toward their substrates (Pereira et al., 2011). This is exemplified by the Hanks-type kinase PrkC from B. subtilis, which phosphorylates the essential translation factor EF-G (Gaidenko et al., 2002) and the enzymes involved in carbohydrate metabolism: the transaldolase YwjH, the glutamine synthetase GlnA, the isocitrate dehydrogenase Icd and the acetolactate-decarboxylase AlsD (Pietack et al., 2010). While recent phosphoproteomic studies in B. subtilis have revealed about a dozen proteins phosphorylated on tyrosine, and many more on serine and threonine (Macek et al., 2007; Soufi et al., 2010), identification of kinases responsible for these phosphorylation events in vivo remains a major challenge. The full complement of protein kinases in B. subtilis includes BY-kinases PtkA and PtkB (Mijakovic et al., 2003), the serine/threonine kinases of the Hanks type PrkC, PrkD, and YabT (Madec et al., 2002; Bidnenko et al., 2013), the unique serine/threonine kinase HprK/P (Galinier et al., 1998), the two-component-like serine/threonine kinases RsbT, RsbW, and SpoIIAB (Garsin et al., 1998; Pane-Farre et al., 2005), and 37 standard two component system kinases (Fabret et al., 1999; Kobayashi et al., 2001). The complement of corresponding phosphatases is much smaller, including the phosphotyrosine-protein phosphatases PtpZ (YwqE), YfkJ, and YwlE (Mijakovic et al., 2005a; Musumeci et al., 2005), and phosphoserine/threonine-protein phosphatases PrpC and SpoIIE (Duncan et al., 1995; Obuchowski et al., 2000). In addition, an arginine phosphorylation system has been described recently in B. subtilis, involving the kinase McsB and the phosphatase YwlE (Elsholz et al., 2012).
To address these crucial questions on the physiological relationships of the kinases, their substrates and alternative kinase modulators, we set out to construct a high confidence protein-protein interaction network centered on known BY-kinases, activators and phosphatases in B. subtilis. Our yeast two-hybrid interactome approach identified numerous candidates of high biological relevance for interactions with protein kinases and cognate phosphatases. In addition, it detected clustering of shared substrates at the interface between kinase and phosphatase nodes. It confirmed several known kinase-substrate relationships, and revealed a number of new ones, notably substrates involved in DNA replication, transcriptional regulation and cell division. Most importantly, we identified the division site selection factor MinD as an anchoring protein which modulates PtkA activity. Our approach also pointed to the existence of cross-talk between the BY-kinase and the Hanks-type Ser/Thr kinase interaction networks, revealing that the two are tightly interconnected. This suggests that bacterial protein kinases (BY- and Ser/Thr Hanks-type) could be regulated by complex signaling cascades, similar to the ones found so far only in eukarya. This study constitutes the most comprehensive view of a bacterial signaling network based on protein phosphorylation to date.
Material and Methods
Construction of Vectors for Protein Expression
All PCR amplifications were performed using B. subtilis 168 genomic DNA as template and primers as listed in Table S4. To construct the knock-out minD mutant, a PCR-amplified internal fragment of the minD gene was inserted into pMUTIN2 (Vagner et al., 1998). To obtain the CFP-N terminally fused proteins, minD and ptkA gene coding sequences were first inserted into a pSG1911 plasmid derivative (Feucht and Lewis, 2001) prior to insertion at the amyE platform locus. For protein expression and purification, PCR fragments were inserted into plasmid pQE-30 Xa (Qiagen) to get the 6xHis-tag fusion proteins. Strep-tagged versions of proteins were obtained from a pQE-30 vector with His6-tag replaced by the strep-tag (Jers et al., 2010). For yeast two hybrid the PCR-amplified gene coding sequences (listed in Table S4-A) were cloned into the bait vector pGBDU-C1, in frame with the binding domain of Gal4 (James et al., 1996), and transformed in PJ69-4A yeast haploid strain (James et al., 1996).
Strain Constructions
Stains used in this study as well as primers employed for constructions are precised in Table S4-A,B. The zapA::SpcR disruption mutation was transferred into ΔptkA background by transformation of the competent ΔptkA cells with chromosomal DNA from LH165 (L. Hamoen; Kawai and Ogasawara, 2006) to spectinomycin resistance. The minD gene was inactivated using the vector pMUTIN2 (Vagner et al., 1998). MinD inactivation mutation was then transferred into the ΔptkA background by transformation with pMUTIN::minD to erythromycin resistance.
Yeast Two Hybrid Screenings of A B. subtilis Genomic Library
The bait vectors were used to screen the pGAD-expressed B. subtilis genomic library in yeast using a previously described mating strategy (Noirot-Gros et al., 2002; Marchadier et al., 2011). Colonies were then tested for their ability to grow on SC –LUH and -LUA plates (lacking leucine, uracil, histidine or adenine). The prey candidates were identified by PCR amplification and sequencing of the DNA inserts from the pGAD plasmids. The potential false-positive interactions were eliminated experimentally using a specificity assays as previously described (Noirot-Gros et al., 2002). About 3–5 × 107 yeast haploid cells expressing the BD-bait protein were mated against 3 x108 haploid cells of complementary mating type, expressing the B. subtilis genomic library of AD-fusion prey protein fragments (Noirot-Gros et al., 2002). Mating efficiencies (15–30%) allowed to test up to 108 possible interactions in one screen. For each screen, the mixture of mated cells was plated on rich medium and incubated for 4–5 h at 30°C, then collected, washed and spread on selective SC-LUH (synthetic complete medium lacking leucine, uracil, and histidine) plates supplemented with 0.5 mM 3-aminotriazole. After 7–10 days at 30°C, the His+ colonies were transferred on SC-LUA (synthetic complete medium lacking leucine, uracil, and adenine) plates for 3 days at 30°C. All the His+ Ade+ colonies were subjected to PCR amplification targeting AD prey-inserts were sequenced and compared with the B. subtilis genome.
Interaction Specificity Assay of Identified Bait-Preys Protein Pairs
False-positive interactions were experimentally removed by a specificity assay. Main sources of false positives in yeast two hybrid are (i) auto-activation of the reporter genes independently to the ability of the targeted protein to bind to a protein partner, (ii) proteins with stickiness properties and (iii) modification of the basal reporter gene expression in the yeast diploid cell. This last category mostly gives non reproducible growth phenotypes. False positive were eliminated experimentally by rescuing the protein-encoding prey plasmids from the His+ Ade+ colonies and reintroducing in PJ69–4α strain in order to perform binary interaction test assays. The rescue of the prey plasmid from diploid cells was carried out using a “gap-repair” cloning procedure as described elsewhere (Weir and Keeney, 2014). Using oligonucleotides priming upstream and downstream of the DNA insert, the prey gene fusion flanked by 150 pb homologous to the vector was amplified by PCR and co-transformed with pGAD linearized vector (digested with appropriated restriction sites within the MCS), in fresh yeast haploid cells. The yeast cell expressing the recombination-based reassembled prey plasmid were then mated against compatible haploid cells carrying (i) an empty pGBDU vector expressing and infused Gal4-BD domain, (ii) a vector expressing the targeted BD-bait protein and (iii) expressing more than three unrelated BD-protein fusions. An interaction is qualified as specific when reproduced twice with the original bait protein, without exhibiting auto-activation nor any unspecific stickiness phenotypes. As already observed (Noirot-Gros et al., 2002) 100% of the prey candidates identified as overlapping fragments of a given gene (class C1, Table S1) were found to be specific while only 20% of the prey proteins found as an unique fragment (class C2, Table S1) have passed the specificity test. In total, a given protein pair qualified as “specific” would have trigger interacting phenotypes up to three time (once during the genomic screen and at least twice during the matricial specificity test assay). The resulting network of specific bait-prey interactions is thus strongly enriched in complex with biological relevance. The protein interactions from this publication have been submitted to the IMEx consortium through IntAct (http://www.imexconsortium.org) (Orchard et al., 2014) and assigned the identifier IM-22270.
PPI Network Construction
In the network the specific protein pairs were represented as nodes connected with edges using the open source software tool Cytoscape 3.0. The network was built from data on connected bait and preys as listed in Table S2.
Proteins Expression and Purification
E. coli K12 NM522 and M15 (expressing chaperonin GroEL/GroES) were used for vector construction and protein purification, respectively. Cells were routinely grown in LB medium supplemented with appropriated antibiotics. Protein synthesis and purifications were carried out as described previously (Mijakovic et al., 2003). Induction of expression was initiated at OD600 = 0.6 by the addition of IPTG (1 mM). Cells were harvested 3 h later and sonicated. From crude extracts, the 6xHis- or Strep-tagged proteins were purified using Ni-NTA (Qiagen), or Strep-Tactin affinity chromatography (Novagen), prior to desalting by PD-10 columns (GE Healthcare). Protein purity was estimated by scanning densitometry of coomassie stained SDS-Page gels (Table S5).
Far-Western Blotting
Proteins were separated by 12% SDS-PAGE and transferred onto nitrocellulose membranes. Proteins were then renatured (Wu et al., 2007) and incubated with MinD labeled with Strep-Tactin (conjugated to horse radish peroxidase, IBA) (Machida and Mayer, 2009) in TBS-T buffer (50 mM Tris, pH 7.5, 150 mM NaCl, 0.05% Tween 20) containing 1% BSA. For labeling, 50 μl of 0.1 μg/μl Strep-tagged MinD was incubated with 5 μ g Strep-Tactin on ice for 1 h. Signal was detected by AEC staining Kit (Sigma).
In vitro Phosphorylation Assays
Proteins were incubated at 37°C in a reacting buffer (50 mM TrisHCl pH 7.5, 100 mM NaCl, 5 mM MgCl2, 5% glycerol) containing 50 μM ATP and 20 μCi/mmol [γ−32P]-ATP, for 1 h, prior to migration on an 8–12% SDS-polyacrylamide gel. The radioactive signal was captured by a FUJI phosphoimager (Mijakovic et al., 2003). In vitro phosphorylation of RecA by PtkA was performed by mixing PtkA (2.5 μM) and RecA (5 μM) in the absence or presence of TkmA (2.5 μM). Phosphorylation and dephosphorylation of RecA by YabT and SpoIIE, was performed by incubating RecA (1 μM) with YabT (0.5 μM). Then SpoIIE (5 μM) was added to the reaction. Phosphorylation and dephosphorylation of RacA by YabT and SpoIIE was performed by mixing RacA (0.8 μM) with YabT (0.15 μM) at 37°C for 1 h. Centricon was used to remove ATP from the sample, prior to addition of SpoIIE (2 μM) and DnaK (0.8 μM). Phosphorylation of DnaC by PrkD was performed by mixing 1 μM of DnaC/I stable co-purified complex with 0.1 μM PrkD. Only DnaC (50.6 kD vs. 30.1 kD for DnaI) was found to be phosphorylated by PrkD. MinD potential kinase activity was assayed by mixing equimolar concentrations of TkmA, PtkA K59M, and MinD (2 μM). In vitro phosphorylation of DivIVA was assayed by mixing equimolar concentrations of PtkA and TkmA (2 μM) with DivIVA (3.5 μM). MinD-mediated activation of PtkA was tested by incubating PtkA (2 μM) with 0.5, 1.0, 2.0, 4.0, and 8.0 μM of MinD. All the in vitro phosphorylation assays were reproduced at least three time.
Fluorescent Microscopy
Strains were grown in liquid LB medium containing appropriated antibiotic overnight at 37°C prior to 1/1000 dilution in LB. Cell were then grown up to OD600 0.3. Samples were mounted on 1.2% agarose pads for examination of living cells by epifluorescent microscopy. The plasma membrane of B. subtilis cells were labeled using the vital stain FM 4–64 and nucleoids were stained with DAPI. CFP/GFP fluorescence was observed using appropriated dichroic filter sets.
Results
A Protein-Protein Interaction Network Centered on BY-Kinases
A network of high biological significance was constructed using a yeast two-hybrid interactome walking approach, comprising a strategy for systematic elimination of false positives (Noirot-Gros et al., 2002; Marchadier et al., 2011). We first performed genome-wide yeast two-hybrid screens of a B. subtilis library to identify the proteins that physically associated with the BY-kinases PtkA and PtkB, their cognate modulators TkmA and TkmB, and the three identified phosphotyrosine-protein phosphatases PtpZ (Mijakovic et al., 2005a), YwlE and YfkJ (Musumeci et al., 2005). To increase the likelihood of identifying potential kinase substrates, we also used PtkA mutant derivatives as baits, carrying point mutations either within the catalytic pocket (D81A or K59M) or within the C-terminal cluster with autophosphorylation sites (Y225F, Y227F, and Y228F) (Mijakovic et al., 2003).
In addition to the expected interactions between the BY-kinases and their respective modulators PtkA-TkmA and PtkB-TkmB, we observed a cross-connection between PtkB and TkmA as well as between TkmA and TkmB (Figure S1). The phosphotyrosine-protein phosphatase PtpZ appeared to interact with TkmA, whereas YwlE and YfkJ, the two low molecular weight phosphotyrosine-protein phosphatases (LMPTPs), lay separately from the BY-kinase and modulators (Figure S1). These findings are in agreement with previous studies indicating that PtpZ is the phosphatase dedicated to dephosphorylating substrates of PtkA/TkmA (Mijakovic et al., 2003, 2005a, 2006). B. subtilis LMPTPs presumably have a different role in stress resistance (Musumeci et al., 2005). Furthermore, we noticed that PtpZ shares common interactants with PtkB, suggesting that PtpZ might also act on putative substrates of PtkB (Figure S1). Common prey proteins were also found between the phosphatases YwlE and YfkJ, indicating they may have redundant activities. The BY-kinase PtkA also interacted with SalA and MinD, two MRP-like/MinD-family ATPases that bear structural resemblance with BY-kinases (Mijakovic et al., 2005b). Additionally, MinD was found to interact with PtkB and PtpZ. Last but not least, this first round of screens also revealed an interaction between the Ser/Thr kinase of the Hanks-type, YabT, and TkmA, the activator of the BY-kinase PtkA.
We then performed a second round of genomic screenings using MinD, SalA, and YabT as baits (Figure S1). In addition, we completed this set with the YabT-cognate phosphoserine-protein phosphatase SpoIIE. Interestingly, SalA and MinD were found to connect with Hanks-type Ser/Thr kinases. Indeed, save for its known interaction with the cell division regulator MinC, MinD also made contact with YabT while SalA interacted with PrkD. Both YabT and SpoIIE were highly connected, thus emerging like hubs in this network. Particularly, YabT made contact with TkmA and PtkB, suggesting the possibility of interference between the serine/threonine and the tyrosine phosphorylation networks. This assumption was strengthened by two outcomes from our third and last rounds of screenings. We observed that YvcJ, a P-loop containing GTPase with identified kinase activity (Pompeo et al., 2011), interacted with TkmB. Furthermore, we identified several connections between the stressosome controlling Ser/Thr protein kinase RsbT (Kang et al., 1998) and kinase and phosphatase proteins from both the Tyr (PtkA, PtpZ) and Ser/Thr (PrkD) phosphorylation pathways (Figure S1).
After addition of connecting edges from the B. subtilis interactome (Marchadier et al., 2011), the resulting interaction map of proteins involved in phosphorylation pathways comprises 137 specific interactions linking 82 proteins (Figure 1A, Tables S1, S2). This network is of high confidence and contains many proteins with already documented phosphosites. These positive controls include recently characterized phosphoproteins RecA, phophorylated by YabT (Bidnenko et al., 2013), and FatR, phosphorylated by PtkA (Derouiche et al., 2013). SrfAA and YqbO, which were detected in a previous phosphoproteome study (Macek et al., 2007), also appear in the network connected with YabT and SpoIIE, respectively. DegS, which was previously characterized as a substrate for both YabT and PrkD, was added to the network (Jers et al., 2011). About 44% of proteins were detected as overlapping fragments, delineating discrete minimal interacting domains (Table S2). These high-confidence interactions provide an important validation of our approach. A significant fraction (20%) of the prey proteins were found to interact with both a kinase and a cognate phosphatase (Table 1), suggesting they could be regulated by reversible phosphorylation. These proteins have been grouped into functional categories (Figure 1B, Table S3). Two thirds of them belong to four main functional classes: post translational modifications, transport of metabolites (ions or peptides), transcription and DNA metabolism.
Figure 1. Protein tyrosine kinases and phosphatases interaction network in B. subtilis. (A) Binary interactions: The network of protein-protein interaction was illustrated using the software Cytoscape 3.0.1. Protein interacting pairs are connected by full lines. Interactions determined in other studies are illustrated by blue lines (Marchadier et al., 2011) or dashed lines (Wu and Errington, 2003; Jers et al., 2011). Reciprocal interactions have been omitted for more clarity (see Table S1). Members of the tyrosine phosphorylation pathway in B. subtilis are indicated as red nodes. Phosphoproteins characterized in the literature are labeled in red. (B) Pie chart diagram illustrating the functional classification of interacting proteins according to the biological processes as shown in Table S4.
New Substrates of BY-Kinases Involved in Transcriptional Control and DNA Metabolism
In the first set of screens, the network revealed the interaction of PtkA-activator TkmA with a transcription regulator FatR. In a separate study we reported that FatR is phosphorylated by PtkA in the presence of TkmA, at residue Y45 located in the DNA binding domain (Derouiche et al., 2013). The consequence of FatR phosphorylation is its dissociation from its DNA operator sequence with derepression of the fatR-cyp102A3 operon. The MRP-like gene regulator SalA which was found to interact with PtkA was also proved to be a substrate for phosphorylation (Derouiche A. and Mijakovic I. unpublished result). In this case PtkA-dependent phosphorylation stimulated DNA-binding of SalA, which in turn repressed the transition state regulator scoC (Derouiche A. and Mijakovic I. unpublished result). Furthermore, we found that TkmA interacted with the N-terminal domain of the general DNA recombinase RecA (Figure 2A). Supporting this observation, we show here that RecA can be phosphorylated by PtkA in the presence of TkmA in vitro (Figure 2B). To date, no direct substrates of the putative BY-kinase PtkB have been identified in B. subtilis, since solubility issues precluded in vitro studies with PtkB (Mijakovic et al., 2003). Our network revealed several interaction partners of PtkB involved in various processes involving the DNA metabolism (Figures 2D, 3C). These include the RNA polymerase β-subunit RpoB, the DNA helicase PcrA, the DNA polymerase PolA, the ATP-dependent DNA helicase/nuclease AddB and the DNA mismatch repair protein MutL. In fact, with the exception of AddB, these proteins also interact with the phosphotyrosine-protein phosphatase PtpZ, highlighting that they indeed could be targeted for modifications through phosphorylation/dephosphorylation on tyrosine residues. RpoB was previously identified as phosphorylated on tyrosine (Ge et al., 2011) and serine (Sun et al., 2010) residues in Helicobacter pylori and Streptococcus pneumoniae, respectively. Both RpoB and the UvrD-like helicase PcrA were also reported to be phosphorylated on tyrosine in Klebsiella pneumoniae (Lin et al., 2009). In our network, PtkB was also found to interact not only with its cognate modulator TkmB but also with TkmA, the cognate modulator of PtkA, suggesting a cross-talk between the two tyrosine kinase systems. In conclusion, our PPI network revealed a number of new substrates of PtkA which we have experimentally validated. In addition it pointed to potential substrates of PtkB.
Figure 2. Phosphorylation of the recombinase protein RecA. (A) Protein-protein interaction between TkmA and the N-terminal domain (NTD) of RecA assayed by yeast two-hybrid (see experimental procedures). AD and BD refer to activating and DNA binding domain of Gal4 fused to TkmA and RecA, respectively. Interacting domains of TkmA and RecA are indicated. (B) In vitro phosphorylation of RecA by PtkA in the presence of TkmA. Experiments were performed as described in the experimental procedures. The presence or absence of RecA is indicated as ± above each lane. Bands corresponding to autophosphorylation of kinases and phosphorylated RecA are indicated by arrows. (C) In vitro phosphorylation and dephosphorylation of RecA mediated by YabT and SpoIIE. Phosphorylation was performed by mixing RecA and YabT as described in the experimental procedures. Note that SpoIIE-mediated dephosphorylation of both RecA (YabT-mediated phosphorylation) and YabT (autophosphorylation) is observed. (D) Protein-protein interactions between tyrosine sphosphorylation and DNA-binding proteins. Proteins are indicated by nodes and interactions by edges. Red nodes indicate the components of protein tyrosine phosphorylation/dephosphorylation machinery. Gray nodes indicate DNA-binding proteins. Identified phosphoproteins are labeled in red.
Figure 3. YabT/SpoIIE centered interacting network highlight membrane protein substrates. (A) Membrane spanning and cytosolic proteins are symbolized by octagonal and round shapes, respectively. Phosphoproteins are labeled in red. (B) Pie-Chart illustrating the repartition of protein substrates in the membrane and cytosol. HTH-containing proteins among the cytosolic proteins are represented. (C) Cross talk between Ser/Thr/Tyr kinases and pathways of DNA metabolism. Interactions between kinases and partners are indicated by nodes connected with edges. Nodes color refers to Tyr-family (red) and Ser/Thr (yellow) family of kinase. PPI are symbolized by black edges; red edges correspond to in vitro phosphorylation validations; blue edges indicate in vitro/in vivo interactions; green edges symbolize high correlation of expressions (as described in Nicolas et al., 2012). Red node labels indicate phosphoproteins.
New Substrates of the Hanks-Type Ser/Thr Kinases and Phosphatases
The interaction of the activator TkmA and the Ser/Thr-kinase YabT established in the first round of screens suggested that the tyrosine and serine/threonine phosphorylation pathways could be connected. This connection was further substantiated in the second round, via the Ser/Thr kinase PrkD, the GTPase YvcJ, the RNA polymease β-subunit RpoB, the recombinase RecA, and the transcriptional factors YesS, and YdeL which were found to connect with both BY-kinase and Ser/Thr-kinase nodes. We explored the Ser/Thr-kinase nodes with the same approach to examine whether PPI network allows identification of new substrates. YabT and its cognate phosphatase SpoIIE are encoded by the same operon and are are strongly up-regulated during sporulation. SpoIIE is known to modulate the phosphorylation state of the anti-anti sigma F factor SpoIIAA (Arigoni et al., 1996). YabT is an unusual Hanks-type kinase possessing a transmembrane helix, but no extracellular signal-binding domain. In a recent study we showed that YabT localizes at the septum during spore formation and is activated by binding ssDNA (Bidnenko et al., 2013). Activated YabT phosphorylates the recombinase RecA at the onset of sporulation and regulates its role in maintaining chromosome integrity. Interaction profiles of YabT and SpoIIE revealed numerous protein partners. Indeed, YabT appears as a hub protein that connects 32 proteins, about half of which (including RecA) also interact with SpoIIE (Figure 3A, Table 2). We speculated that proteins interacting with both a kinase and a phosphatase are very likely to undergo phosphorylation. Supporting this hypothesis, we found that RecA, a known substrate of YabT, can be dephosphorylated by SpoIIE (Figure 2C). The DNA-binding developmental protein RacA was also investigated (Figure 4). During sporulation RacA anchors the segregating chromosome at the cell pole in a DivIVA-dependent manner (Ben-Yehuda et al., 2003). Examination of the condition-dependent transcription profiles of racA and spoIIE shows they are members of the same synexpression group mainly up-regulated during sporulation, indicating that they are involved the same biological process (Figure 4C, Nicolas et al., 2012). The phosphorylation and dephosphorylation of RacA by YabT and SpoIIE, respectively, were confirmed in vitro (Figure 4B). This observation highlights a potential role of the kinase YabT in the regulation of RacA activity, and further validates our assumption that interaction partners of a kinase and a phosphatase are presumably substrates for both enzymes.
Figure 4. Reversible Ser/Thr phosphorylation of RacA. (A) RacA interaction with YabT and SpoIIE by yeast 2HB. Gal4 BD- (line) and AD-fusions (column) expressed in yeast haploid cells of complementary mating type were mated and assayed for expression of interaction phenotype as described in the experimental procedures. Diploids co-expressing an AD-RacA fusion (aa 81-194/194) and BD-YabT or BD-SpoIIE are able to grow onto –LUH selective media. (B) In vitro phosphorylation/dephosphorylation of RacA by YabT, and SpoIIE. Phosphorylation was performed as described in the experimental procedures. (C) Condition-dependent transcriptomes of yabT (black), spoIIE (gray) and racA (light red), as obtained in Nicolas et al. (2012).
Classification of all proteins of the YabT hub revealed that 69% of partners are transmembrane proteins, with a high proportion of metabolite- and ion-transporters. Among cytosolic interacting proteins, 67% exhibited HTH motifs suggesting DNA binding properties (Figure 3B). Of note, this parallels the properties of YabT itself, a transmembrane kinase activated by DNA binding (Bidnenko et al., 2013). In addition to RecA and RacA, YabT/SpoIIE were found to interact with SbcE, a structural maintenance of chromosome (SMC)-like protein involved in DNA double-strand break repair (Krishnamurthy et al., 2010) and CisA, a site-specific DNA recombinase necessary for reconstitution of the sigma factor K during sporulation (Kunkel et al., 1990) (Figure 3C). Given the ability of YabT to be stimulated by ssDNA, our identification of potential substrates with DNA-binding properties and/or proteins involved in DNA metabolism substantiates its observed involvement in maintenance of DNA integrity (Bidnenko et al., 2013).
Further, our PPI network revealed a direct interaction between the YabT paralogue PrkD and the replicative DNA helicase DnaC (Figure 5, Tables S1, S2). Interaction of the kinase with DnaC in yeast two-hybrid requires the DnaC C-terminal domain comprising the helicase function. We confirmed the ability of PrkD to phosphorylate DnaC in vitro (Figure 5B), suggesting that PrkD could be involved in regulating DnaC activity. DnaC loading onto a melted single strand DNA fork is assisted by DnaB and DnaI (56). Analysis of condition-dependent transcriptome of B. subtilis shows that prkD (old name ybdM) and dnaB exhibit a high pairwise correlation throughout ~100 different physiological conditions (Nicolas et al., 2012; Figure 5C), further strengthening the notion that PrkD could participate to the regulation of chromosomal replication. In conclusion, our PPI network revealed that both YabT and PrkD interact with, and phosphorylate, proteins with crucial roles in DNA metabolism.
Figure 5. Phosphorylation of the DNA helicase DnaC. (A) Protein-protein interaction between PrkD and DnaC assayed by yeast two hybrid. AD and BD refer to activating and DNA binding domain of Gal4 fused to PrkD and DnaC, respectively. DnaC-Nter refers to the N-terminal domain of dnaC, spanning residues aa 1–170. (B) In vitro phosphorylation of DnaC by PrkD. Phosphorylation was performed as described in experimental procedures. Bands corresponding to autophosphorylation of PkrD and phosphorylated DnaC are indicated by arrows. (C) Correlations between transcriptional levels of prkD and the components of the DNA helicase loader complex. Pairwise correlations between DnaC, DnaB, and DnaI values were calculated over the 104 physiological conditions (Nicolas et al., 2012).
MinD, a New Modulator of the BY-Kinase PtkA
The MinD ATPase, required for the correct positioning of the division site, was identified as interacting partner of the two BY-kinases, PtkA and PtkB, as well as the phosphotyrosine-protein phosphatase PtpZ (Figure 6A). This pointed to MinD as a potential target for tyrosine phosphorylation. However, we found that MinD was not phosphorylated by PtkA in vitro (data not shown). We first validated the interaction between MinD and PtkA by far Western blotting (Figure 6B). In this assay we detected the expected MinD-MinC complex and confirmed the physical interaction between MinD and PtkA, but not TkmA, in agreement with the two-hybrid. MinD is an ATPase from the P-loop NTPase superfamily previously reported to be a close structural homolog of BY-kinases, but lacking the C-terminal tyrosine cluster necessary for kinase autophosphorylation (Mijakovic et al., 2005b). We therefore explored potential kinase activity of MinD in vitro (Figure 6C). MinD was found unable to autophosphorylate or to phosphorylate a PtkA derivative carrying a substituted catalytic residue (K59M), corroborating that MinD does not act as a kinase. During cytokinesis, MinD is known to functionally interact with the cell division initiator protein DivIVA in the presence of the bridging protein MinJ, by preventing the formation of a septum near the cell pole (Levin et al., 1992; Lee and Price, 1993; Marston et al., 1998). Interestingly, we found that DivIVA could undergo phosphorylation by PtkA in the presence of TkmA (Figure 6D).
Figure 6. Functional analysis of MinD/PtkA interaction. (A) MinD-centered interaction networks. Interactions are illustrated by edges connecting nodes (proteins). Interactions were identified by yeast 2HB (black), genetic (blue) and in vitro (red). (B) Interaction between PtkA and MinD detected by far-western blotting. PtkA, MinC, and TkmA were separated by SDS-PAGE and transferred onto nitrocellulose membrane. Interaction was determined by strep tagged MinD which was revealed by Strep-Tactin. MinC acts as a positive control and TkmA acts as a negative control. (C) MinD is not acting as a kinase. All proteins were mixed in equimolar concentration. The presence or absence of TkmA is indicated as ± above each lane. Bands corresponding to phosphorylated PtkA are indicated by arrows. (D) In vitro phosphorylation of DivIVA by PtkA in the presence of TkmA. The presence or absence of proteins is indicated as ± above each lane. Time dependence of the reaction was followed and the incubation times (15, 30, or60 mins) are given above the lanes. Bands corresponding to autophosphorylated PtkA and phosphorylated DivIVA are indicated by arrows. (E) MinD acts as an activator of PtkA. Autoradiography of SDS-polyacrylamide gels showing the influence of MinD on autophosphorylation of PtkA. (a) PtkA was incubated with increasing concentrations of MinD. The ratio of MinD to PtkA was indicated above each lane. (b) MinD activates PtkA less efficiently than TkmA. PtkA was incubated with equivalent concentration of TkmA or MinD. (c) MinD can compete with TkmA for PtkA activation. TkmA/PtkA were incubated in the presence of increasing amount of MinD. The ratio of MinD to TkmA/PtkA was indicated above each lane. Bands corresponding to autophosphorylated PtkA are indicated by arrows. (F) Yeast two-hybrid interaction between MinD and PtkA. C-terminal domain of TkmA (aa187–248) and C-terminal domains of MinD (I: aa135–268; II: aa89–268) were fused to the activating domain of Gal4 (AD, labeled in black). PtkA wt and mutant derivatives (as indicated) are fused with the binding domain of Gal4 (BD, labeled in green). PGAD and pGBDU are control empty plasmids expressing the Gal4-AD and BD domains.
We then tested the ability of MinD to promote the activation of PtkA in the absence of TkmA in vitro (Figure 6E). Incubation of PtkA with increasing amounts of MinD activated its autophosphorylation, indicating that MinD is able to act as modulator of PtkA activity similar to TkmA (Figure 6Ea). However, incubation of PtkA with TkmA or MinD revealed that MinD-mediated PtkA activation was much less efficient compared to TkmA (Figure 6Eb). Finally, we showed that in the presence of excess of MinD, the activation of PtkA in the presence of TkmA was inhibited, indicating that MinD was able to compete with TkmA for binding to PtkA (Figure 6Ec). We investigated the functional elements of PtkA required for interaction with MinD and TkmA by looking for their ability to interact with several PtkA mutant derivatives in yeast two-hybrid (Figure 6F). The mutants included substitutions of key active site residues (catalytic K59M and Mg-coordinating D81A) and the tyrosine autophosphorylation sites (Y225, Y227, and Y228 changed into F). Mutations within the C-terminal tyrosine cluster did not affect the PtkA interaction with either TkmA or MinD. Interestingly, the substitution of either active site residues specifically impaired PtkA interaction with MinD but not with TkmA. This result points to different modes of binding of TkmA and MinD to PtkA.
PtkA Plays a Role in Cell Division
Transcription profiles of minD, divIVA, ptkA, and tkmA exhibit a high degree of correlation throughout all physiological conditions, suggesting they function in the same biological process (Figure 7A). A close examination of ptkA single mutant revealed a mild cell division defect with a noticeable proportion of cells containing double septa (7%), compare to the wild type cells (<1%) (Figures 8A,B). We investigated for PtkA-associated phenotypes in cells deficient in the division factor zapA. B. subtilis zapA mutant exhibits no remarkable division phenotype (Gueiros-Filho and Losick, 2002). However, a zapA null mutation was described to exacerbate cytokinesis block-induced phenotypes when combined with other actors of cell division (Gueiros-Filho and Losick, 2002; Dempwolff et al., 2012; Surdova et al., 2013). In rich media, zapA ptkA double mutant exhibited a higher proportion of cell with double septa (12%) as well as cells containing asymmetric septa (3%), sometimes leading to anucleated cells (Figures 8A,B). These observations suggest a role of PtkA in the regulation of cytokinesis.
Figure 7. PtkA/MinD cellular expression and localization. (A) Condition-dependent transcription profiles of minD, divIVA, ptkA, and tkmA (Nicolas et al., 2012). PtkA and TkmA expression profiles are in black and gray, respectively. Expressions of DivIVA and MinD are in red (dark and light, respectively). (B) MinD-dependent localization of PtkA: GFP fluorescence (upper panel) and bright field (lower panel) images are shown. Localization of MinD in wild type cells (left panel) and localization of PtkA in wild type (middle panel) and in ΔminD B. subtilis cells (right panel) as indicated. Arrows point to minicells.
Figure 8. PtkA is involed in regulation of correct septal positioning. (A) ΔptkA cells in wild type or zapA deficient background were observed by fluorescent microscopy. Cells were visualized by membrane staining (FM4-64, or left panel), nucleoid staining (Dapi, or right panel), or bright field (BF). Septum defects are indicated by white arrows (double septum) and white triangles (asymmetric septum or anucleated cell). (B) Defaults in septal formation in each strain are estimated as % over a total number of cells >500. AS: asymmetric septum, red; DS: double septum, green.
To further understand the role of MinD in PtkA activity, we examined their subcellular localization. MinD is known to exhibit a very specific localization pattern at the cell poles as well as at mid cell. When fused to CFP, PtkA was found to localize as a single focus, preferentially at the old cell poles of elongating cells (Figure 7B). The absence of MinD typically leads to the formation of anucleated minicells, resulting from the mislocalisation of septum at the cell poles (Burmann et al., 2011). Here we found that the inactivation of minD also abolished the polar localization of PtkA foci, indicating it is MinD-dependant. This result provides in vivo validation of the MinD/PtkA interaction and hints at a role of MinD in holding PtkA at the cell pole during cytokinesis.
Discussion
A hallmark of BY-kinases from Gram-positive bacteria is that the cytosolic kinase and its transmembrane activator are individual proteins encoded by separate adjacent genes. This differs from their Gram-negative counterparts where the two activities are encoded by the same gene. A functional consequence of this particular architectural feature could be that it enables the kinases from Gram-positive bacteria to interact with alternative activators, which may provide different substrate specificities. We presented here a protein-protein interaction landscape of the tyrosine phosphorylation network of B. subtilis. Combined with a systematic and experimental elimination of false positives, our approach has proved effective for identifying numerous substrates of BY-kinases and phosphatases, as well as new proteins that modulate their activity. This PPi network not only provided a kinase context for already identified phosphoproteins, but also identified new protein partners for kinases, including other kinases and phosphatases. This network can be considered of high biological relevance, with numerous proteins being contacted by both a kinase and cognate phosphatase, strongly suggesting they are true substrates or they belong to the same regulatory assembly. Our examination of kinase/substrates and kinase/modulators from the tyrosine but also the serine/threonine phosphorylation pathways corroborated these assumptions. This study highlights the role of protein phosphorylation in the regulation of various aspects of the B. subtilis cell cycle as DNA replication, transcription, cell division and sporulation.
Common signaling mechanisms often involve multisite protein phosphorylation to regulate the various activities of a protein whether carried out by the same kinase or by the sequential action of several kinases. In B. subtilis, the best characterized example is the two-component histidine kinase DegS, which was shown to be phosphorylated by the two Ser/Thr kinases PrkD and YabT in vitro (Jers et al., 2011). PrkD-dependent phosphorylation was specific to serine 76 located in the signal sensing domain and led to increased efficiency of phosphotransfer to DegU. In vivo, DegS phosphorylation at serine 76 was shown to affect various aspects of cellular behavior such as competence, swarming and motility. Our B. subtilis interaction network also revealed that other proteins can be phosphorylated by more than one kinase (Table S3-B). A notable example is the multifunctional DNA recombinase protein RecA recently found regulated by phosphorylation on a serine (Lusetti and Cox, 2002; Butala et al., 2009; Bidnenko et al., 2013). Here, our data suggest that tyrosine phosphorylation could also play a role in the regulation of RecA activity. RecA binding to TkmA involves its N-terminal moiety, which is required for the formation of a presynaptic nucleoprotein filament (Lee and Wang, 2009). Whether this activity could be modulated by tyrosine phosphorylation in the cell remains to be determined. Dual phosphorylation involving modification on tyrosine, serine or threonine could also be proposed for the transcriptional regulators SalA and YesS as well as for the RNA polymerase subunit RpoB, which was found to interact with more than one component of the Ser/Thr and Tyr phosphorylation pathways. Together, these observations suggest the existence of dynamic regulations of protein activities in B. subtilis cells, where multiple and sequential phosphorylation events could mediate integration of several input signals to better adjust the cellular response to environmental changes, similar to what is described in eukarya.
Remarkably, our network reveals important crosstalk among different classes of kinases. Interactions were observed between components of the Tyrosine and Ser/Thr phosphorylation pathways, as well as among Ser/Thr-kinases. These connections strongly suggest that kinases might phosphorylate each other, generating signaling cascades similar to what was observed in eukarya. This hypothesis has been explored in a separate study, in which the ability of all soluble B. subtilis BY- and Ser/Thr kinases to phosphorylate each other in vitro was explored. This assay revealed many trans-phosphorylations at functionally critical residues indicative of kinase regulation by phosphorylation (Shi et al., 2014). A similar study carried out recently in M. tuberculosis also provided strong evidence for a hierarchical and spatially organized crosstalk among Ser/Thr kinases in this bacteria (Baer et al., 2014). These results support the existence of kinase cross-phosphorylation to form complex eukaryotic-like signaling networks in bacteria. Furthermore, promiscuous binding of BY-kinases PtkA and PtkB with different modulators was observed throughout the network suggesting a functional interplay. Whether PtkB performs in vivo as a kinase remains to be determined. An alternative hypothesis could be that PtkB acts as a pseudokinase to mediate interaction with protein substrates for subsequent tyrosine modification by PtkA. Our data also hints at a potential a role of modulators in recruiting some substrates. In agreement with this notion is the finding that the transcriptional regulator FatR, which interacts with TkmA, is regulated by PtkA (Derouiche et al., 2013). On the whole, our PPi network highlights elaborate mechanisms for phosphorylation of tyrosine in the signal transduction pathways that regulate various aspects of DNA metabolism.
A remarkable feature of this PPi network is the central position of the cell division regulator protein MinD which appears as a hub protein densely connected with many members of the tyrosine phosphorylation pathways. MinD shares sequence and structural homology with BY-kinases, but lacks crucial tyrosine residues involved in phosphoryltransfer activity (Mijakovic et al., 2005b). As such, MinD could be considered as a pseudokinase with a non enzymatic-like function. Our resultes indicate a role of MinD in the activation of PtkA in vitro, as well as in its in vivo localization as a focus at the cell poles. It has been previously reported that PtkA co-localizes with its transmembrane adapter TkmA at the cell membrane in exponentially growing cells (Jers et al., 2010). However, in the previous study ptkA and tkmA genes were simultaneously overexpressed which may account for this discrepancy. In our present study the localization was detected with wild type expression levels from the ptkA natural promoter, without overexpression of tkmA. Based on these results, we propose that MinD could act as a platform protein that would tether PtkA at a specific cellular area. Additionally MinD would also act as an activator of PtkA for signal transmission to DivIVA. The biological significance of DivIVA phosphorylation by PtkA during cystokinesis in B. subtilis remains to be investigated, but our data clearly points to a role of PtkA in the regulation of cytokinesis. This is in agreement with the cell cycle phenotype of the ptkA mutant which was reported previously (Petranovic et al., 2007). A parallel can be drawn between the MinD/PtkA interplay and the A parallel can be drawn with scaffold-directed assembly of signaling components in eukarya. Eukaryotic scaffold proteins are often described as catalytically inactive kinases that can bring together multiple components of signaling cascades, regulating and promoting their interactions at particular localizations in the cell. An example of this process is the Ste5 protein scaffold in the mitogen-activated protein kinase (MAPK) pathway, which directs mating signaling by promoting an activation cascade mediated by three insulated protein kinases MAPKKK Ste11, MAPKK ste7 and MAPK Fus3 (Good et al., 2009). Here, MinD was found connected with other kinases besides PtkA (PktB, YabT) and a phosphatase (PtpZ) thus bringing various components of tyrosine, but also serine/threonine phosphorylation pathways in close proximity. A hypothesis worth exploring would be whether MinD could orchestrate the intracellular location of several kinases and phosphatases, similar to eukaryal anchoring proteins (Pawson and Scott, 1997).
In several bacteria, the phosphorylation of DivIVA was shown to be a key step in the regulation of cytokinesis. So far, this regulation is known to be mediated by Ser/Thr Hank-type kinases. In Streptomyces coelicolor, the polarization machinery is regulated by AfsK, a Ser/Thr protein kinase which localizes at hyphal tips and phosphorylates DivIVA on its C-terminus (Hempel et al., 2012; Saalbach et al., 2013). In S. pneumoniae, StkP localizes at the division site and specifically phosphorylates DivIVA on threonine 201 (Fleurie et al., 2012). In Mycoplasma pneumoniae, phosphorylation of the DivIVA homolog Wag31 by S/T kinases PknA/PknB affects cell shape control (Kang et al., 2005). Our results indicate that in B. subtilis, this regulation may involve the phosphorylation of DivIVA on a tyrosine residue.
In conclusion, this study shows that bacterial regulatory network based on protein phosphorylation is considerably more complex and inter-connected than previously believed. The BY-kinase sub-network is inextricably connected to the Ser/Thr-kinase sub-network. Together with the separate study (Shi et al., 2014) providing evidence that kinases can cross-phosphorylate each other, our data highlighted that BY-kinases are activated by alternative modulators, and that substrates are phosphorylated by different kinases, leading to a better integration of a multitude of signals. Significant additional efforts will be required to decipher the functioning of this complex and dynamic regulatory network.
Conflict of Interest Statement
The Associate Editor, Christophe Grangeasse, declares that, despite having collaborated with authors Ivan Mijakovic and Lei Shi, the review process was handled objectively and no conflict of interest exists. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Stephen McGovern for providing purified DnaC protein. This work was supported by the Agence Nationale de la Recherche (2010-BLAN-1303-01 to Ivan Mijakovic and Marie-Françoise Noirot-Gros) and the Federation pour la Recherche Medicale (SPF20111223359 to Ivan Mijakovic and Lei Shi). The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Protein-protein interaction data reported in this paper has been submitted to the IMEx consortium through IntAct (http://www.imexconsortium.org) and assigned the identifier IM-22270.
Supplementary Material
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb.2014.00538/abstract
References
Arigoni, F., Duncan, L., Alper, S., Losick, R., and Stragier, P. (1996). SpoIIE governs the phosphorylation state of a protein regulating transcription factor sigma F during sporulation in Bacillus subtilis. Proc. Natl. Acad. Sci. U.S.A. 93, 3238–3242. doi: 10.1073/pnas.93.8.3238
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Baer, C. E., Iavarone, A. T., Alber, T., and Sassetti, C. M. (2014). Biochemical and spatial coincidence in the provisional Ser/Thr protein kinase interaction network of mycobacterium tuberculosis. J. Biol. Chem. 289, 20422–20433. doi: 10.1074/jbc.M114.559054
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Ben-Yehuda, S., Rudner, D. Z., and Losick, R. (2003). RacA, a bacterial protein that anchors chromosomes to the cell poles. Science 299, 532–536. doi: 10.1126/science.1079914
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Bidnenko, V., Shi, L., Kobir, A., Ventroux, M., Pigeonneau, N., Henry, C., et al. (2013). Bacillus subtilis serine/threonine protein kinase YabT is involved in spore development via phosphorylation of a bacterial recombinase. Mol. Microbiol. 88, 921–935. doi: 10.1111/mmi.12233
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Burmann, F., Ebert, N., Van Baarle, S., and Bramkamp, M. (2011). A bacterial dynamin-like protein mediating nucleotide-independent membrane fusion. Mol. Microbiol. 79, 1294–1304. doi: 10.1111/j.1365-2958.2011.07523.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Butala, M., Zgur-Bertok, D., and Busby, S. J. (2009). The bacterial LexA transcriptional repressor. Cell Mol. Life Sci. 66, 82–93. doi: 10.1007/s00018-008-8378-6
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Dempwolff, F., Wischhusen, H. M., Specht, M., and Graumann, P. L. (2012). The deletion of bacterial dynamin and flotillin genes results in pleiotrophic effects on cell division, cell growth and in cell shape maintenance. BMC Microbiol. 12:298. doi: 10.1186/1471-2180-12-298
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Derouiche, A., Bidnenko, V., Grenha, R., Pigonneau, N., Ventroux, M., Franz-Wachtel, M., et al. (2013). Interaction of bacterial fatty-acid-displaced regulators with DNA is interrupted by tyrosine phosphorylation in the helix-turn-helix domain. Nucleic Acids Res. 41, 9371–9381. doi: 10.1093/nar/gkt709
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Duncan, L., Alper, S., Arigoni, F., Losick, R., and Stragier, P. (1995). Activation of cell-specific transcription by a serine phosphatase at the site of asymmetric division. Science 270, 641–644. doi: 10.1126/science.270.5236.641
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Elsholz, A. K., Turgay, K., Michalik, S., Hessling, B., Gronau, K., Oertel, D., et al. (2012). Global impact of protein arginine phosphorylation on the physiology of Bacillus subtilis. Proc. Natl. Acad. Sci. U.S.A. 109, 7451–7456. doi: 10.1073/pnas.1117483109
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Fabret, C., Feher, V. A., and Hoch, J. A. (1999). Two-component signal transduction in Bacillus subtilis: how one organism sees its world. J. Bacteriol. 181, 1975–1983.
Feucht, A., and Lewis, P. J. (2001). Improved plasmid vectors for the production of multiple fluorescent protein fusions in Bacillus subtilis. Gene 264, 289–297. doi: 10.1016/S0378-1119(01)00338-9
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Fleurie, A., Cluzel, C., Guiral, S., Freton, C., Galisson, F., Zanella-Cleon, I., et al. (2012). Mutational dissection of the S/T-kinase StkP reveals crucial roles in cell division of Streptococcus pneumoniae. Mol. Microbiol. 83, 746–758. doi: 10.1111/j.1365-2958.2011.07962.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Gaidenko, T. A., Kim, T. J., and Price, C. W. (2002). The PrpC serine-threonine phosphatase and PrkC kinase have opposing physiological roles in stationary-phase Bacillus subtilis cells. J. Bacteriol. 184, 6109–6114. doi: 10.1128/JB.184.22.6109-6114.2002
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Galinier, A., Kravanja, M., Engelmann, R., Hengstenberg, W., Kilhoffer, M. C., Deutscher, J., et al. (1998). New protein kinase and protein phosphatase families mediate signal transduction in bacterial catabolite repression. Proc. Natl. Acad. Sci. U.S.A. 95, 1823–1828. doi: 10.1073/pnas.95.4.1823
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Garsin, D. A., Duncan, L., Paskowitz, D. M., and Losick, R. (1998). The kinase activity of the antisigma factor SpoIIAB is required for activation as well as inhibition of transcription factor sigmaF during sporulation in Bacillus subtilis. J. Mol. Biol. 284, 569–578. doi: 10.1006/jmbi.1998.2202
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Ge, R., Sun, X., Xiao, C., Yin, X., Shan, W., Chen, Z., et al. (2011). Phosphoproteome analysis of the pathogenic bacterium Helicobacter pylori reveals over-representation of tyrosine phosphorylation and multiply phosphorylated proteins. Proteomics 11, 1449–1461. doi: 10.1002/pmic.201000649
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Gerwig, J., Kiley, T. B., Gunka, K., Stanley-Wall, N., and Stulke, J. (2014). The protein tyrosine kinases EpsB and PtkA differentially affect biofilm formation in Bacillus subtilis. Microbiology 160, 682–691. doi: 10.1099/mic.0.074971-0
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Good, M., Tang, G., Singleton, J., Remenyi, A., and Lim, W. A. (2009). The Ste5 scaffold directs mating signaling by catalytically unlocking the Fus3 MAP kinase for activation. Cell 136, 1085–1097. doi: 10.1016/j.cell.2009.01.049
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Grangeasse, C., Cozzone, A. J., Deutscher, J., and Mijakovic, I. (2007). Tyrosine phosphorylation: an emerging regulatory device of bacterial physiology. Trends Biochem. Sci. 32, 86–94. doi: 10.1016/j.tibs.2006.12.004
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Gueiros-Filho, F. J., and Losick, R. (2002). A widely conserved bacterial cell division protein that promotes assembly of the tubulin-like protein FtsZ. Genes Dev. 16, 2544–2556. doi: 10.1101/gad.1014102
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Hanks, S. K., Quinn, A. M., and Hunter, T. (1988). The protein kinase family: conserved features and deduced phylogeny of the catalytic domains. Science 241, 42–52. doi: 10.1126/science.3291115
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Hempel, A. M., Cantlay, S., Molle, V., Wang, S. B., Naldrett, M. J., Parker, J. L., et al. (2012). The Ser/Thr protein kinase AfsK regulates polar growth and hyphal branching in the filamentous bacteria Streptomyces. Proc. Natl. Acad. Sci. U.S.A. 109, E2371–E2379. doi: 10.1073/pnas.1207409109
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Jadeau, F., Grangeasse, C., Shi, L., Mijakovic, I., Deleage, G., and Combet, C. (2012). BYKdb: the Bacterial protein tYrosine Kinase database. Nucleic Acids Res. 40, D321–D324. doi: 10.1093/nar/gkr915
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
James, P., Halladay, J., and Craig, E. A. (1996). Genomic libraries and a host strain designed for highly efficient two-hybrid selection in yeast. Genetics 144, 1425–1436.
Jers, C., Kobir, A., Sondergaard, E. O., Jensen, P. R., and Mijakovic, I. (2011). Bacillus subtilis two-component system sensory kinase DegS is regulated by serine phosphorylation in its input domain. PLoS ONE 6:e14653. doi: 10.1371/journal.pone.0014653
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Jers, C., Pedersen, M. M., Paspaliari, D. K., Schutz, W., Johnsson, C., Soufi, B., et al. (2010). Bacillus subtilis BY-kinase PtkA controls enzyme activity and localization of its protein substrates. Mol. Microbiol. 77, 287–299. doi: 10.1111/j.1365-2958.2010.07227.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Kang, C. M., Abbott, D. W., Park, S. T., Dascher, C. C., Cantley, L. C., and Husson, R. N. (2005). The Mycobacterium tuberculosis serine/threonine kinases PknA and PknB: substrate identification and regulation of cell shape. Genes Dev. 19, 1692–1704. doi: 10.1101/gad.1311105
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Kang, C. M., Vijay, K., and Price, C. W. (1998). Serine kinase activity of a Bacillus subtilis switch protein is required to transduce environmental stress signals but not to activate its target PP2C phosphatase. Mol. Microbiol. 30, 189–196. doi: 10.1046/j.1365-2958.1998.01052.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Kawai, Y., and Ogasawara, N. (2006). Bacillus subtilis EzrA and FtsL synergistically regulate FtsZ ring dynamics during cell division. Microbiology 152, 1129–1141. doi: 10.1099/mic.0.28497-0
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Kiley, T. B., and Stanley-Wall, N. R. (2010). Post-translational control of Bacillus subtilis biofilm formation mediated by tyrosine phosphorylation. Mol. Microbiol. 78, 947–963. doi: 10.1111/j.1365-2958.2010.07382.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Kobayashi, K., Ogura, M., Yamaguchi, H., Yoshida, K., Ogasawara, N., Tanaka, T., et al. (2001). Comprehensive DNA microarray analysis of Bacillus subtilis two-component regulatory systems. J. Bacteriol. 183, 7365–7370. doi: 10.1128/JB.183.24.7365-7370.2001
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Krishnamurthy, M., Tadesse, S., Rothmaier, K., and Graumann, P. L. (2010). A novel SMC-like protein, SbcE (YhaN), is involved in DNA double-strand break repair and competence in Bacillus subtilis. Nucleic Acids Res. 38, 455–466. doi: 10.1093/nar/gkp909
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Kunkel, B., Losick, R., and Stragier, P. (1990). The Bacillus subtilis gene for the development transcription factor sigma K is generated by excision of a dispensable DNA element containing a sporulation recombinase gene. Genes Dev. 4, 525–535. doi: 10.1101/gad.4.4.525
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Lee, C. D., and Wang, T. F. (2009). The N-terminal domain of Escherichia coli RecA have multiple functions in promoting homologous recombination. J. Biomed. Sci. 16:37. doi: 10.1186/1423-0127-16-37
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Lee, S., and Price, C. W. (1993). The minCD locus of Bacillus subtilis lacks the minE determinant that provides topological specificity to cell division. Mol. Microbiol. 7, 601–610. doi: 10.1111/j.1365-2958.1993.tb01151.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Levin, P. A., Margolis, P. S., Setlow, P., Losick, R., and Sun, D. (1992). Identification of Bacillus subtilis genes for septum placement and shape determination. J. Bacteriol. 174, 6717–6728.
Lin, M. H., Hsu, T. L., Lin, S. Y., Pan, Y. J., Jan, J. T., Wang, J. T., et al. (2009). Phosphoproteomics of Klebsiella pneumoniae NTUH-K2044 reveals a tight link between tyrosine phosphorylation and virulence. Mol. Cell. Proteomics 8, 2613–2623. doi: 10.1074/mcp.M900276-MCP200
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Lusetti, S. L., and Cox, M. M. (2002). The bacterial RecA protein and the recombinational DNA repair of stalled replication forks. Annu. Rev. Biochem. 71, 71–100. doi: 10.1146/annurev.biochem.71.083101.133940
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Macek, B., Mijakovic, I., Olsen, J. V., Gnad, F., Kumar, C., Jensen, P. R., et al. (2007). The serine/threonine/tyrosine phosphoproteome of the model bacterium Bacillus subtilis. Mol. Cell. Proteomics 6, 697–707. doi: 10.1074/mcp.M600464-MCP200
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Machida, K., and Mayer, B. J. (2009). Detection of protein-protein interactions by far-western blotting. Methods Mol. Biol. 536, 313–329. doi: 10.1007/978-1-59745-542-8_34
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Madec, E., Laszkiewicz, A., Iwanicki, A., Obuchowski, M., and Seror, S. (2002). Characterization of a membrane-linked Ser/Thr protein kinase in Bacillus subtilis, implicated in developmental processes. Mol. Microbiol. 46, 571–586. doi: 10.1046/j.1365-2958.2002.03178.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Marchadier, E., Carballido-Lopez, R., Brinster, S., Fabret, C., Mervelet, P., Bessieres, P., et al. (2011). An expanded protein-protein interaction network in Bacillus subtilis reveals a group of hubs: exploration by an integrative approach. Proteomics 11, 2981–2991. doi: 10.1002/pmic.201000791
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Marston, A. L., Thomaides, H. B., Edwards, D. H., Sharpe, M. E., and Errington, J. (1998). Polar localization of the MinD protein of Bacillus subtilis and its role in selection of the mid-cell division site. Genes Dev. 12, 3419–3430. doi: 10.1101/gad.12.21.3419
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Mijakovic, I., Musumeci, L., Tautz, L., Petranovic, D., Edwards, R. A., Jensen, P. R., et al. (2005a). In vitro characterization of the Bacillus subtilis protein tyrosine phosphatase YwqE. J. Bacteriol. 187, 3384–3390. doi: 10.1128/JB.187.10.3384-3390.2005
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Mijakovic, I., Petranovic, D., Bottini, N., Deutscher, J., and Ruhdal Jensen, P. (2005b). Protein-tyrosine phosphorylation in Bacillus subtilis. J. Mol. Microbiol. Biotechnol. 9, 189–197. doi: 10.1159/000089647
Mijakovic, I., Petranovic, D., Macek, B., Cepo, T., Mann, M., Davies, J., et al. (2006). Bacterial single-stranded DNA-binding proteins are phosphorylated on tyrosine. Nucleic Acids Res. 34, 1588–1596. doi: 10.1093/nar/gkj514
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Mijakovic, I., Poncet, S., Boel, G., Maze, A., Gillet, S., Jamet, E., et al. (2003). Transmembrane modulator-dependent bacterial tyrosine kinase activates UDP-glucose dehydrogenases. EMBO J. 22, 4709–4718. doi: 10.1093/emboj/cdg458
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Musumeci, L., Bongiorni, C., Tautz, L., Edwards, R. A., Osterman, A., Perego, M., et al. (2005). Low-molecular-weight protein tyrosine phosphatases of Bacillus subtilis. J. Bacteriol. 187, 4945–4956. doi: 10.1128/JB.187.14.4945-4956.2005
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Nicolas, P., Mader, U., Dervyn, E., Rochat, T., Leduc, A., Pigeonneau, N., et al. (2012). Condition-dependent transcriptome reveals high-level regulatory architecture in Bacillus subtilis. Science 335, 1103–1106. doi: 10.1126/science.1206848
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Noirot-Gros, M. F., Dervyn, E., Wu, L. J., Mervelet, P., Errington, J., Ehrlich, S. D., et al. (2002). An expanded view of bacterial DNA replication. Proc. Natl. Acad. Sci. U.S.A. 99, 8342–8347. doi: 10.1073/pnas.122040799
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Obuchowski, M., Madec, E., Delattre, D., Boel, G., Iwanicki, A., Foulger, D., et al. (2000). Characterization of PrpC from Bacillus subtilis, a member of the PPM phosphatase family. J. Bacteriol. 182, 5634–5638. doi: 10.1128/JB.182.19.5634-5638.2000
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Orchard, S., Ammari, M., Aranda, B., Breuza, L., Briganti, L., Broackes-Carter, F., et al. (2014). The MIntAct project–IntAct as a common curation platform for 11 molecular interaction databases. Nucleic Acids Res. 42, D358–D363. doi: 10.1093/nar/gkt1115
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Pane-Farre, J., Lewis, R. J., and Stulke, J. (2005). The RsbRST stress module in bacteria: a signalling system that may interact with different output modules. J. Mol. Microbiol. Biotechnol. 9, 65–76. doi: 10.1159/000088837
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Pawson, T., and Scott, J. D. (1997). Signaling through scaffold, anchoring, and adaptor proteins. Science 278, 2075–2080. doi: 10.1126/science.278.5346.2075
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Pereira, S. F., Goss, L., and Dworkin, J. (2011). Eukaryote-like serine/threonine kinases and phosphatases in bacteria. Microbiol. Mol. Biol. Rev. 75, 192–212. doi: 10.1128/MMBR.00042-10
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Petranovic, D., Michelsen, O., Zahradka, K., Silva, C., Petranovic, M., Jensen, P. R., et al. (2007). Bacillus subtilis strain deficient for the protein-tyrosine kinase PtkA exhibits impaired DNA replication. Mol. Microbiol. 63, 1797–1805. doi: 10.1111/j.1365-2958.2007.05625.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Pietack, N., Becher, D., Schmidl, S. R., Saier, M. H., Hecker, M., Commichau, F. M., et al. (2010). In vitro phosphorylation of key metabolic enzymes from Bacillus subtilis: PrkC phosphorylates enzymes from different branches of basic metabolism. J. Mol. Microbiol. Biotechnol. 18, 129–140. doi: 10.1159/000308512
Pompeo, F., Luciano, J., Brochier-Armanet, C., and Galinier, A. (2011). The GTPase function of YvcJ and its subcellular relocalization are dependent on growth conditions in Bacillus subtilis. J. Mol. Microbiol. Biotechnol. 20, 156–167. doi: 10.1159/000329298
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Saalbach, G., Hempel, A. M., Vigouroux, M., Flardh, K., Buttner, M. J., and Naldrett, M. J. (2013). Determination of phosphorylation sites in the DivIVA cytoskeletal protein of Streptomyces coelicolor by targeted LC-MS/MS. J. Proteome Res. 12, 4187–4192. doi: 10.1021/pr400524d
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Shi, L., Kobir, A., Jers, C., and Mijakovic, I. (2010). Bacterial Protein-Tyrosine Kinases. Curr. Proteomics 7, 188–194. doi: 10.2174/157016410792928198
Shi, L., Pigeonneau, N., Ravikumar, V., Dobrinic, P., Macek, B., Franjevic, D., et al. (2014). Cross-phosphorylation of bacterial serine/threonine and tyrosine protein kinases on key regulatory residues. Front. Microbiol. 5:495. doi: 10.3389/fmicb.2014.00495
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Soufi, B., Kumar, C., Gnad, F., Mann, M., Mijakovic, I., and Macek, B. (2010). Stable isotope labeling by amino acids in cell culture (SILAC) applied to quantitative proteomics of Bacillus subtilis. J. Proteome Res. 9, 3638–3646. doi: 10.1021/pr100150w
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Sun, X., Ge, F., Xiao, C. L., Yin, X. F., Ge, R., Zhang, L. H., et al. (2010). Phosphoproteomic analysis reveals the multiple roles of phosphorylation in pathogenic bacterium Streptococcus pneumoniae. J. Proteome Res. 9, 275–282. doi: 10.1021/pr900612v
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Surdova, K., Gamba, P., Claessen, D., Siersma, T., Jonker, M. J., Errington, J., et al. (2013). The conserved DNA-binding protein WhiA is involved in cell division in Bacillus subtilis. J. Bacteriol. 195, 5450–5460. doi: 10.1128/JB.00507-13
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Vagner, V., Dervyn, E., and Ehrlich, S. D. (1998). A vector for systematic gene inactivation in Bacillus subtilis. Microbiology 144 (Pt 11), 3097–3104. doi: 10.1099/00221287-144-11-3097
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Weir, M., and Keeney, J. B. (2014). PCR mutagenesis and gap repair in yeast. Methods Mol. Biol. 1205, 29–35. doi: 10.1007/978-1-4939-1363-3_3
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Whitfield, C. (2006). Biosynthesis and assembly of capsular polysaccharides in Escherichia coli. Annu. Rev. Biochem. 75, 39–68. doi: 10.1146/annurev.biochem.75.103004.142545
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Wu, L. J., and Errington, J. (2003). RacA and the Soj-Spo0J system combine to effect polar chromosome segregation in sporulating Bacillus subtilis. Mol. Microbiol. 49, 1463–1475. doi: 10.1046/j.1365-2958.2003.03643.x
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Wu, Y., Li, Q., and Chen, X. Z. (2007). Detecting protein-protein interactions by Far western blotting. Nat Protoc. 2, 3278–3284. doi: 10.1038/nprot.2007.459
Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar
Keywords: protein phosphorylation, protein-protein network, bacterial protein kinase, bacterial cell division, bacterial signaling
Citation: Shi L, Pigeonneau N, Ventroux M, Derouiche A, Bidnenko V, Mijakovic I and Noirot-Gros M-F (2014) Protein-tyrosine phosphorylation interaction network in Bacillus subtilis reveals new substrates, kinase activators and kinase cross-talk. Front. Microbiol. 5:538. doi: 10.3389/fmicb.2014.00538
Received: 19 August 2014; Accepted: 26 September 2014;
Published online: 22 October 2014.
Edited by:
Christophe Grangeasse, Centre National de la Recherche Scientifique, FranceReviewed by:
Bettina Siebers, University of Duisburg-Essen, GermanyImrich Barak, Slovak Academy of Sciences, Slovakia
Copyright © 2014 Shi, Pigeonneau, Ventroux, Derouiche, Bidnenko, Mijakovic and Noirot-Gros. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marie-Françoise Noirot-Gros, Institut National de la Recherche Agronomique, UMR1319 Micalis, Domaine de Vilvert, 78350 Jouy-en-Josas, France e-mail: marie-francoise.gros@jouy.inra.fr
†These authors have contributed equally to this work.