- Research article
- Open Access
Biochemical and functional characterization of SpdA, a 2′, 3′cyclic nucleotide phosphodiesterase from Sinorhizobium meliloti
BMC Microbiologyvolume 13, Article number: 268 (2013)
3′, 5′cAMP signaling in Sinorhizobium meliloti was recently shown to contribute to the autoregulation of legume infection. In planta, three adenylate cyclases CyaD1, CyaD2 and CyaK, synthesizing 3′, 5′cAMP, together with the Crp-like transcriptional regulator Clr and smc02178, a gene of unknown function, are involved in controlling plant infection.
Here we report on the characterization of a gene (smc02179, spdA) at the cyaD1 locus that we predicted to encode a class III cytoplasmic phosphodiesterase.
First, we have shown that spdA had a similar pattern of expression as smc02178 in planta but did not require clr nor 3′, 5′cAMP for expression.
Second, biochemical characterization of the purified SpdA protein showed that, contrary to expectation, it had no detectable activity against 3′, 5′cAMP and, instead, high activity against the positional isomers 2′, 3′cAMP and 2′, 3′cGMP.
Third, we provide direct experimental evidence that the purified Clr protein was able to bind both 2′, 3′cAMP and 3′, 5′cAMP in vitro at high concentration. We further showed that Clr is a 3′, 5′cAMP-dependent DNA-binding protein and identified a DNA-binding motif to which Clr binds. In contrast, 2′, 3′cAMP was unable to promote Clr specific-binding to DNA and activate smc02178 target gene expression ex planta.
Fourth, we have shown a negative impact of exogenous 2′, 3′cAMP on 3′, 5′cAMP-mediated signaling in vivo. A spdA null mutant was also partially affected in 3′, 5′cAMP signaling.
SpdA is a nodule-expressed 2′, 3′ specific phosphodiesterase whose biological function remains elusive. Circumstantial evidence suggests that SpdA may contribute insulating 3′, 5′cAMP-based signaling from 2′, 3′ cyclic nucleotides of metabolic origin.
Sinorhizobium meliloti is a soil-born α-proteobacterium that can enter a nitrogen-fixing symbiosis with Medicago sativa (alfalfa) and related legumes. The establishment of the symbiosis relies on a complex molecular dialogue between the two partners that triggers two essential and overlapping steps, nodulation and infection (see [1, 2] for reviews). During the infection process, bacteria colonize root hairs forming Infection Threads (ITs) that extend and proliferate towards the nodule primordium that is formed in the root cortex. Ultimately, rhizobia are released from ITs within nodule cells where they fix molecular dinitrogen. Nodulation and infection are tightly controlled processes and we have shown recently that bacterial adenylate cyclases (ACs) contribute to the negative autoregulation of infection .
ACs (EC 188.8.131.52) are enzymes that synthesize cAMP (3′, 5′-cyclic adenosine monophosphate) from ATP. There are 6 non-homologous classes of ACs as a typical example of convergent evolution [4, 5]. Class III is the universal class whose members can be found in both prokaryotes and eukaryotes although, to our knowledge, their presence in plants has not been established . The number of class III ACs strikingly varies in bacteria. E. coli has none whereas cyanobacteria, mycobacteria and rhizobia, a group of phylogenetically-diverse bacteria , have many, up to 32 in the soybean symbiont Bradyrhizobium japonicum. The biological function of class III ACs in bacteria remains poorly understood. Class III ACs synthesize cAMP in response to environmental cues such as light, oxygen, nitrogen and pH in Cyanobacteria  or high osmotic pressure in Myxococcus xanthus[9, 10]. Class III ACs are also involved in biotic interactions as they contribute to virulence in M. tuberculosis, P. aeruginosa and in some fungal pathogens [5, 11–13]. CO2 and Ca2+ are signals used by pathogens to sense their host environment through their AC–cAMP signaling systems. Candida albicans and mycobacteria express CO2-responsive ACs [5, 14] whereas CyaB from P. aeruginosa is Ca2+ sensitive. Another example of cAMP-associated signal being used by the human fungal pathogen C. albicans to sense the host environment is the bacterial peptidoglycan present in blood serum .
We have recently described the first instance of class III ACs contributing to a symbiotic (mutualistic) interaction, between Sinorhizobium meliloti and its host plant Medicago sativa. S. meliloti has 26 class III ACs of overall unknown biological functions with a variety of domain organization . In response to a plant signal present in nodules, three receptor-like adenylate cyclases CyaD1, CyaD2 and CyaK synthesize the secondary messenger molecule 3′, 5′cAMP. 3′, 5′cAMP together with the Crp-like transcriptional activator Clr in turn promote transcription of the target gene smc02178, of unknown biochemical function . We have recently found that this cascade contributes to the autoregulation of the symbiotic interaction. Specifically, activation of the cAMP cascade in nodules inhibits, by a mechanism that remains to be elucidated, secondary infection by rhizospheric bacteria. This control is lost in either a triple cyaD1cyaD2cyaK mutant, a clr or a smc02178 mutant resulting in a hyper-infection phenotype on plants–ie an abundance of abortive ITs on roots–as a consequence of a relaxed control of secondary infection .
The concentration of the second messenger 3′, 5′cAMP in cells is controlled at the level of its synthesis by ACs and/or by its degradation to 5′AMP by phosphodiesterases (PDEs). PDEs are a superfamily of enzymes divided in three, non-homologous, main classes. All mammalian PDEs as well as several enzymes identified in Drosophila, Caenorhabditis and Saccharomyces cerevisiae belong to class I, whose conserved carboxy-terminal catalytic domain contains two invariant motifs H(X)3H(X)25-35D/E . Class II PDEs are enzymes from Saccharomyces cerevisiae, Dictyostelium discoideum, Schizosaccharomyces pombe, C. albicans, and Vibrio fischeri. This class of enzymes shares the conserved motif HXHLDH. Class III PDEs belong to the superfamily of metallophosphoesterases . They share the conserved sequence motif D-(X)n-GD(X)n-GNH[E/D]-(X)n-H-(X)n-GHXH as well as a βαβαβ secondary structure signature .
Here we report on the characterization of a class III PDE from S. meliloti (SpdA, SMc02179) that we anticipated from the localization of the spdA gene at the cyaD1 locus to be involved in signal termination by turning-over the secondary messenger 3′, 5′cAMP. We have found that purified SpdA had actually no detectable activity against 3′, 5′cAMP and, instead, had high activity on the structural isomer 2′, 3′cAMP, which may occur in cells as a by-product of RNA degradation . We demonstrated that, contrary to 3′, 5′cAMP that promoted Clr binding to a cognate binding-site, 2′, 3′cAMP bound unproductively to Clr. Although SpdA biological function remains to be established, we present circumstantial evidence that SpdA may insulate 3′, 5′cAMP-mediated signaling from 2′, 3′-structural isomers.
SpdA, a putative PDE
Inspection of the cyaD1 locus (Figure 1A), that contains the clr gene as well as the clr–target gene smc02178, pointed to the smc02179 gene product as a potential PDE that we subsequently coined SpdA. SpdA belongs to a 15-member protein family sharing the IPR004843 domain characteristic of a wide range of metallophosphoesterases, among which phosphorine phosphatases, nucleotidases, and class III PDEs. We thus compared SpdA as well as the 14 other IPR004843-containing proteins to known PDEs from Mycobacterium tuberculosis (Rv0805), Haemophilus influenzae (Icc) and Escherichia coli (CpdA and CpdB) [20–22].
Overall analysis of the whole protein family indicated no clear phylogenetic relationship between the family members besides the fact that SMc04449 and SMc04018 behaved as an outgroup together with CpdB, a periplasmic 2′, 3′ cAMP-PDE from E. coli (see Additional file 1). SpdA closest homologue was M. tuberculosis Rv0805 and indeed closer sequence inspection indicated that SpdA contained the 5 sub-domains characteristic of Rv0805 and other class III PDEs  (Figure 1B) whereas all other S. meliloti proteins, except SMc02712, had fewer (see Additional file 1). SpdA had a predicted cytoplasmic location and missed the amino-terminal 200-aminoacid membrane anchoring domain of Rv0805 .
spdA is expressed in planta, independently of clrand 3’, 5’cAMP
We probed expression of a translational spdA-lacZ fusion (pGD2179, See Additional file 2) that contained the intergenic region between smc02178 and spdA (Figure 1A) as well as the first 12 codons of spdA. The spdA-lacZ fusion did not detectably express ex planta and instead expressed in Medicago sativa nodules with the same pattern as smc02178i.e. expression in young nodule primordia and in zones II and III of mature nodules (Figure 2A-F). However, spdA expression in planta was independent of clr, and ex planta expression could not be induced by exogenous 3′, 5′cAMP, in contrast to smc02178 expression (Figure 2G). None of the environmental conditions or compounds which we have tested was able to stimulate spdA expression ex planta, including 3′, 5′cGMP, 2′, 3′cAMP, 5′AMP, nodule extracts, root exudates or several growth and stress conditions (See Additional file 3).
Altogether these results indicated that spdA was expressed in planta from its own promoter and had the same expression pattern as smc02178 although the two genes were not co-regulated.
SpdA is a 2′, 3′cNMP PDE
We purified the SpdA protein as a carboxy-terminal His6-tagged fusion (Figure 3A). Under non-denaturing electrophoretic conditions the protein migrated as a monomer. Purified His6-SpdA protein displayed activity against the generic PDE substrate BispNPP in vitro (Figure 3B). SpdA had little or no activity against either 3′, 5′cAMP or 3′, 5′cGMP but significantly hydrolyzed the positional isomers 2′, 3′cAMP and 2′, 3′cGMP (Figure 3C) which are products of RNA degradation . The Km for 2′, 3′cAMP was 3.7 mM and kCat was 2 s-1 indicating a slow enzyme with low affinity for its substrate in vitro (See Additional file 4). We observed no inhibition of the enzyme by its substrate and found that 3′, 5′cAMP did not affect SpdA activity on 2′, 3′cAMP.
Despite IPR004843-containing proteins being documented metalloenzymes, the metal chelators EDTA, 1-10-Phenanthroline and Bipyridyl, or the addition of Fe2+ or Mn2+ metal ions, had no effect on SpdA activity (see Additional file 5). Mass spectrometry of isolated SpdA confirmed the absence of associated metal including Mg2+, Mn2+ and Co2+ together with the monomeric state of the protein. Indeed, a well resolved single mass peak corresponding to the monomer was observed after Max-Ent deconvolution of the spectra.
2′, 3′cAMP binds unproductively to Clr
In order to investigate a possible interference of 2′, 3′cyclic nucleotides with 3′, 5′ cAMP-signaling we assessed the capacity of 2′, 3′cAMP and 3′, 5′cAMP to bind Clr in vitro. For this purpose, we purified a GST-tagged version of Clr by affinity purification (Figure 4A). Purified Clr protein was loaded onto a 3′, 5′cAMP-agarose column. Bound Clr protein was then eluted with either the cognate 3′, 5′cAMP nucleotide or its 2′, 3′ isomer (30 mM). Both nucleotides displaced agarose-bound Clr thus suggesting that Clr could bind 3′, 5′cAMP and 2′, 3′cAMP at the same binding site (Figure 4B, C).
Clr is a predicted transcriptional activator of the Crp family . Inspection of the smc02178 promoter region pointed to a short palindromic sequence (TGTTCCGCGGGAAACA) centered ca. 68 bp upstream of the predicted start codon that was a potential binding site for Clr. Accordingly, deletion of this motif abolished activation of the smc02178 promoter by clr in the presence of exogenously provided 3′, 5′cAMP (Figure 5A). In order to directly assess whether this motif was a binding site for the Clr protein, we tested the ability of purified Clr-GST to bind DNA oligomers (28-mers) bracketing the putative Clr-binding motif (Figure 5B) or a mutated version (Figure 5C). We found that Clr induced a retard in oligomer migration that was strictly dependent on the presence of 3′, 5′cAMP, of an intact Clr-box and was Clr concentration-dependent. However, no clear shifted band was observed, irrespectively of the binding and gel electrophoresis conditions tested, which probably reflected dissociation of the Clr/cAMP/DNA complex. Nevertheless we interpreted this as evidence that Clr bound the predicted Clr-box in a 3′, 5′cAMP-dependent manner. 2′, 3′cAMP was unable to promote Clr binding to the Clr-box, at the same concentration as 3′, 5′cAMP. Mixed incubation of the two nucleotides (1/1) with Clr in vitro showed no detectable effect of 2′, 3′cAMP on DNA-binding by Clr (Figure 6A, B).
We tested the impact of exogenously provided 2′, 3′cAMP on smc02178 expression in vivo under different experimental conditions. Exogenous 2′, 3′cAMP alone was unable to promote activation of the smc02178-lacZ reporter fusion in vivo, even at high (7.5 mM) concentration (Figure 6C). In contrast 2′, 3′cAMP had a negative impact on 3′, 5′cAMP-driven smc02178 expression. Inhibition reached 50% (Figure 6C) when 3′, 5′cAMP was produced endogenously, as in normal physiological conditions, upon addition to the bacterial culture of a Medicago shoot extract containing the plant signal that triggers activity of the CyaD1CyaD2CyaK ACs . Inhibition was only 30% when 3′, 5′ cAMP was provided exogenously (See Additional file 6). Noteworthy, the negative impact of 2′, 3′cAMP was not observed on a constitutive hemA-lacZ reporter fusion (pXLGD4, see Additional file 2 and Additional file 6) suggesting a specific effect of 2′, 3′cAMP on 3′, 5′cAMP-mediated signaling.
Biological characterization of a S. meliloti spdAnull mutant
As to get an insight into SpdA biological function we inactivated the corresponding gene by cre-lox deletion . spdA inactivation decreased smc02178-lacZ expression by ca. 25% in the presence of plant shoot extracts, supposedly by increasing endogenous 2′, 3′cNMP concentration in vivo. Combining spdA inactivation together with exogenous 2′, 3′cAMP addition decreased smc02178 expression to 40% of wild-type (Figure 6C and See Additional file 6).
The spdA mutant had the same growth characteristics as wild-type both in rich complex medium (LBMC) and in synthetic Vincent medium with mannitol and glutamate (VGM) as carbon and nitrogen sources (see Additional file 7). We observed that exogenous 2′, 3′cAMP extended bacterial growth in VGM medium, suggesting that S. meliloti can grow by utilizing 2′, 3′cAMP, as Yersinia does . However the spdA mutant did not differ from wild-type in this respect. The spdA mutant also responded similarly to wild-type to various stress conditions including detergent (SDS) and heat shock (See Additional file 7).
spdA inactivation had no detectable effect on symbiotic performances, including nodulation, infection and nitrogen fixation (plant dry weight), on Medicago sativa nor on the level or pattern of smc02178 symbiotic expression in planta (See Additional file 8).
Hence we did not detect any phenotype associated with the spdA mutation besides its limited effect on 3’, 5’ cAMP-signaling.
Clr is a 3′, 5′cNMP-dependent DNA-binding transcriptional activator
The findings reported here give experimental support and extend the model proposed by , as we demonstrated that Clr binds to the smc02178 promoter region at a specific site in a 3′, 5′cAMP-dependent manner. The transcription start site (TSS) at the smc02178 promoter was not determined experimentally here. However a single smc02178 TSS was mapped in the closely related strain 2011 by RNA-sequencing of a pool of bacteria living in 16 different free-living and stress conditions . The TSS mapped 61.5 bp downstream of the center of the Clr-box which is the distance typically found in class I Crp(CAP)-dependent promoters. In Class I promoters, a single protein–protein interaction with CAP facilitates the binding of RNA-Polymerase to the promoter to yield the RNA-Polymerase–promoter closed complex .
One salient feature of Clr binding at the smc02178 promoter DNA was instability. In spite of the many binding and electrophoresis conditions tested, we consistently observed a smear instead of a clear-cut band shift upon binding of Clr to its target DNA. One feature that may account for this instability is that the Clr binding site is TGTTN8AACA, a shorter palindrome as compared to the consensus E. coli CRP(CAP)-binding site TGTGAN6TCACA. Identification of this binding motif, together with transcriptome analysis experiments, will help identification of new Clr targets in the S. meliloti genome.
The reason for which 2′, 3′cAMP did not promote DNA-binding of Clr is unclear. Although Clr bound 2′, 3′cAMP in vitro at high concentration (30 mM), it may not do so at the concentration of 2′, 3′cAMP that we used in EMSA assays (200 μM). Alternatively, 2′, 3′cAMP may not trigger the appropriate conformational change that allows Crp binding to DNA. Further experiments are needed to distinguish between these two possibilities.
SpdA encodes a 2′, 3′cNMP phosphodiesterase
Class III PDEs are metallophosphoesterases carrying the IPR004843 domain. IPR004843-containing proteins have a wide range of substrates, including cyclic nucleotides, and ensure a variety of biological functions . S. meliloti has 15 uncharacterized IPR004843-containing proteins (see Additional file 1) and we have demonstrated that purified SpdA has a PDE activity in vitro (Figure 3).
We have further found that SpdA had no or little activity against 3′, 5′cAMP or 3′, 5′cGMP and instead had high activity against 2′, 3′cAMP or 2′, 3′cGMP. Although this cannot be formally excluded it is unlikely that SpdA would have a predominant 3′, 5′cAMP PDE activity in vivo since a SpdA null mutant had lower, and not enhanced, smc02178 expression in vivo (Figure 6C).
Substrate specificity varies widely among class III PDEs. CpdA from E. coli and P. aeruginosa, Icc from Haemophilus influenzae are 3′, 5′cNMP PDEs [21, 22, 29] whereas E. coli CpdB was the first described 2′, 3′cNMP-specific PDE . Rv0805 from M. tuberculosis, although it was first reported as a 3′, 5′cNMP PDE , has a much stronger activity (150 times fold) against 2′, 3′cNMP than against 3′, 5′cNMP . Myxococcus xanthus PdeA and PdeB instead hydrolyse 2′, 3′cNMP and 3′, 5′cNMP with the same affinity . Hence class III PDEs substrate specificity cannot be predicted from simple primary sequence inspection. It is thus possible that several IPR004843 proteins of S. meliloti display a 2′, 3′cyclic phosphodiesterase activity, thus contributing a functional redundancy.
A surprising feature of SpdA was the absence of associated metal ion which is, to our knowledge, unique among IPR004843-containing proteins. Rv0805 activity for example was not inhibited by metal chelators but was boosted by Mn2+ addition . However, it has been already reported that the iPGM protein from castor bean that belongs to a superfamily of metalloenzymes  was actually metal-independent . Moreover, the Carboxy-terminal HD domain of the E. coli tRNA nucleotidyltransferase has a metal-independent phosphodiesterase activity toward 2′, 3′cAMP . Thus, the fact that SpdA displays metal-independent 2′, 3′cNMP-phosphodiesterase activity is not completely unprecedented. Mass spectrometric measurements performed under mild ionization conditions also pointed out that the well-defined monomeric form of the protein did not present any demetallation.
The 2′, 3′cNMP substrate specificity of SpdA leaves the question of 3′, 5′cAMP turnover intact. One option would be to identify a 3′, 5′cNMP PDE among the 14 other S. meliloti proteins containing the IPR004843 domain. Another, non-exclusive, possibility would be a regulation of 3′, 5′cAMP homeostasis by secretion rather than by degradation .
Possible biological functions for SpdA
Very little is known about the origin, role and fate of 2′, 3′cyclic nucleotides. One documented origin is RNA degradation and physiological or stressful conditions may indeed lead 2′, 3′cNMPs to accumulate in bacteria. We are not aware of any other origin such as, for example, isomerization of corresponding 3’, 5’ cyclic nucleotides. In this context, SpdA may serve at least three different, non-exclusive, functions: a metabolic function, a detoxifying function and a role in preventing cross talk with 3′, 5′cAMP signaling.
Although S. meliloti likely metabolized exogenous 2′, 3′cAMP (See Additional file 7), spdA was not critical for this since the mutant grew indistinctly from wild-type under these conditions.
2′, 3′cAMP was recently reported to be a toxic compound in kidney cells, that opens mitochondria permeability transition pores thus leading to a pre-apoptotic and necrotic stage . We thus considered whether SpdA may counteract a toxic effect of 2′, 3′cNMPs in S. meliloti. However the unaltered growth characteristics of the spdA mutant as compared to wild-type in various growth (including the presence of exogenous 2′, 3′cAMP) and stress conditions (see Additional file 7) did not give support to this possibility.
A third possibility would be SpdA preventing cross-talk between 2′, 3′cyclic nucleotides and 3′, 5′cAMP signaling. Several lines of evidence are in favor of this possibility: (i) the evolutionary-conserved physical location of spdA in close proximity to cyaD1, clr and the target gene smc02178 in all the sequenced strains of Sinorhizobium meliloti, Sinorhizobium saheli and Sinorhizobium fredii (https://www.genoscope.cns.fr/agc/microscope/mage/);  (ii) spdA expression in nodules at the same place where 3′, 5′cAMP signaling takes places  and where a massive RNA degradation occurs as part of the reorientation of the bacteroid transcriptome to the goal of nitrogen fixation  (iii) a significant and specific decrease in smc02178 expression upon providing exogenous 2′, 3′cAMP (iv) the spurious interaction of 2′, 3′cAMP with Clr.
Whatever SpdA function, the high Km value measured in vitro for the 2′, 3′cAMP substrate (3.7 mM) would imply that the cyclic nucleotide accumulates in high amounts in bacteroids, unless specific physiological or biochemical conditions lower Km value in vivo. Developing methods for direct measurements of 2′, 3′cNMP levels in bacteroids, where spdA preferentially expresses, is now needed to clarify this issue. A ribonucleic origin for 2′, 3′cAMP/cGMP would make sense physiologically given the extensive transcriptome reprofiling taking place in bacteroids  and the abundance of VapC-type ribonucleases in S. meliloti genome . Intriguingly, the human intracellular pathogen M. tuberculosis shares with S. meliloti, despite the large phylogenetic distance separating them, a wealth of ACs, a Clr-like transcriptional regulator as well as a close homolog of SpdA, Rv0805. Rv0805, like SpdA, has a preferential activity–and similar Km value-towards 2′, 3′ cyclic nucleotides  and contributes to overall bacterial virulence on macrophages, by a still obscure mechanism [11, 12, 24]. Interestingly, M. tuberculosis and S. meliloti have in commun a high number of VapC-type RNases of the VapC(B)-toxin (antitoxin) family [40, 41]. Altogether this suggests the intriguing possibility that SpdA, Rv0805 and other cytoplasmic PDEs may constitute a physiological adaptation in bacteria with a high RNA turnover, possibly in relationship with 3′, 5′cAMP-mediated signaling.
Signal transduction in bacteria is dominated by two-component regulatory systems . However, some bacteria, including important pathogens and symbionts, use cyclic or dicyclic nucleotide signaling for modulating interaction with their abiotic or biotic environment [43, 44]. Characterization of enzymes and mechanisms that synthesize and degrade secondary messenger molecules, restrict their diffusion within the cell and prevent cross-talking by diffusible isomers, is needed for fully understanding cyclic nucleotide signaling. In the context of characterizing 3′, 5′cAMP-mediated signaling in the S. meliloti-Medicago symbiosis, we have identified a plant-expressed 2′, 3′cAMP/cGMP specific phosphodiesterase whose biological function remains to be elucidated. Circumstantial evidence suggests that one SpdA function could be to insulate 3′, 5′cAMP-based signaling from 2′, 3′ cyclic nucleotides of metabolic origin.
Bacterial strains, plasmids, and growth conditions
Plasmids and bacterial strains used in this study are listed in Additional file 2 and Additional file 9 respectively. S. meliloti strains were grown at 28°C in rich LB medium supplemented with 2.5 mM CaCl2 and 2.5 mM MgSO4 (LBMC) or in modified Vincent synthetic medium with glutamate (0.1%) and mannitol (1%) as nitrogen and carbon sources, respectively (VGM) . E. coli strains were grown at 37°C in rich LB medium.
The concentrations of antibiotics used for S. meliloti cultures were 200 μg/ml for streptomycin, 100 μg/ml for neomycin, 10 μg/ml for tetracycline, and 30 μg/ml for gentamicin. The concentrations of antibiotic used for E. coli cultures were 50 μg/ml for ampicillin and 25 μg/ml for kanamycin.
Bacterial response to SDS and heat shock was evaluated by analysis of the growth curves of WT and ΔSpdA mutant in liquid LBMC. Strains were challenged with SDS (0.01% v/v) at OD600 0.1 and heat shock (50°C for 20 min) was applied to overnight cultures before dilution at OD600 0.1. Aliquots were collected at different time intervals, OD600 was measured and residual growth was determined .
Construction of plasmids and mutant strains
Primers used for DNA amplification are listed in Additional file 10. S. meliloti 1021 was used as template for DNA amplification. For deletion of the spdA gene, we used the cre-lox system . PCR fragments encompassing the upstream/amino-terminal coding region and the downstream/carboxyl-terminal coding region of spdA were amplified using CreLox 2179 up Left-CreLox 2179 up Right and 2179 Down NcoI-2179 Down HincII as primers (See Additional file 10), digested by SacI-SacII and NcoI-HincII, and cloned into the SacI-SacII and NcoI-HincII restriction sites of pCM351, respectively. The resulting plasmid was introduced into the S. meliloti 1021 strain by conjugation. Transconjugants sensitive to tetracycline and resistant to gentamicin were screened. A ΔspdA mutant was selected. The spdA-expressing construct pET::2179 was obtained after amplification of the spdA gene-coding region using S. meliloti 1021 genomic DNA as template and LNdeI2179 and RHindIII 2179 as primers. The PCR fragment was digested with NdeI and HindIII and cloned into the NdeI-HindIII digested pET-22b plasmid to yield pET::2179. The Clr-expressing construct pGEX::clr was obtained after amplification of the clr gene-coding region using S. meliloti 1021 genomic DNA as template and ClrBamHI and ClrEcoRI as primers. The PCR fragment was digested with BamHI and EcoRI and cloned into the BamHI-EcoRI digested pGEX-2T to yield pGEX::clr.
To construct pGD2179, that carries a spdA-lacZ translational fusion, a 177-bp PCR fragment encompassing the spdA promoter region was amplified using 2179left and 2179right primers, digested with HindIII and BamHI, and cloned in the in-frame orientation at the same sites of the lacZ translational fusion plasmid pGD926. The pAMG2178 plasmid was obtained after amplification of the smc02178 promoter-coding region using S. meliloti 1021 genomic DNA as template and BamHI 2178 and Hind BoxL as primers. For pAMG2178ΔClrbox, PCR fragments encompassing the upstream region Clr box and the downstream region Clr box of the smc02178 promoter were amplified using 2178 H-BoxLPstI and X 2178-BoxRPstI as primers. The two fragments obtained were digested by PstI and then ligated and amplified by PCR using BamHI 2178 and Hind BoxL as primers. The two fragments p2178WT (134 bp) and p2178ΔClrbox (128 bp) were then cloned into pGD926 vector.
All constructs were verified by PCR and Sanger sequencing in E. coli and by PCR in S. meliloti. Plasmids were transferred from E. coli to S. meliloti by triparental mating using pRK600 as the helper plasmid. pET::2179 and pGEX::clr were directly transferred into E. coli BL21(DE3) and SP850 respectively.
For His6-SpdA purification, an overnight culture of E. coli strain BL21(DE3) pET::2179 expressing wild-type S. meliloti spdA was diluted at OD600 0.1 in 250 ml of LB medium containing Ampicillin (Amp 50 μg/ml). Cultures were grown with shaking at 28°C. When the OD600 reached 0.8, 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) was added, and cultures were grown for 5 additional hours. Bacteria were collected by centrifugation (10,000x g for 30 min at 4°C), and pellets were washed with 60 ml Tris buffer (20 mM Tris–HCl [pH 8.0]). Bacteria were collected by centrifugation (10,000x g for 30 min at 4°C), and pellets were stored at−80°C. All of the subsequent procedures were performed at 4°C. Thawed bacteria were resuspended in 5 ml of buffer A (50 mM Tris–HCl [pH 8.0], 250 mM NaCl, 10% glycerol) and lysed by sonication. The lysates were centrifuged to remove the cell debris at 10,000x g for 30 min at 4°C. The supernatant was loaded to a Ni-NTA resin (Qiagen) equilibrated with buffer B (50 mM Tris–HCl [pH 8.0], 250 mM NaCl, 10% glycerol, 10 mM Imidazol, and 5 mM β-Mercaptoethanol). After washing with the buffer B containing 20 mM Imidazol, the bound protein was eluted using the buffer B containing 250 mM Imidazol. Protein was desalted into buffer A. Purified protein aliquots were stored at−80°C.
For Clr-GST purification, an overnight culture of E. coli strain SP850 pGEX::clr expressing wild-type S. meliloti clr was diluted at OD600 0.1 in 1 l of LB medium containing Ampicillin (Amp 50 μg/ml) and Kanamycin (Kan 25 μg/ml). Cultures were grown with shaking at 28°C. When the OD600 reached 0.8, 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) was added, and cultures were grown for 5 additional hours. Bacteria were collected by centrifugation (10,000x g for 30 min at 4°C), and pellets were washed with 60 ml PBS buffer (140 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4, [pH 7.3]). Bacteria were collected by centrifugation (10,000x g for 30 min at 4°C), and pellets were stored at−80°C. All of the subsequent procedures were performed at 4°C. Thawed bacteria were resuspended in 10 ml PBS buffer and lysed by sonication. The lysates were centrifuged to remove the cell debris at 10,000x g for 30 min at 4°C. The supernatant was loaded to a Glutathione sepharose 4B resin (GE Healthcare) equilibrated with PBS buffer. After washing with PBS buffer, the bound protein was eluted using 50 mM Tris–HCl buffer [pH 8.0] containing 10 mM reduced glutathione. Protein was desalted on Amicon CO 10,000 (Millipore) and buffer exchanged with 0.1 M Phosphate buffer [pH 7.0] containing 50 mM NaCl. Purified protein aliquots containing 10% glycerol were stored at−80°C.
Infusion ESI-Q-TOF experiment
ElectroSpray Ionization coupled with a quadrupole-time of flight tandem was used in the positive ion mode using a Q-TOF Ultima Instrument (Waters). The SpdA protein was dissolved in water with 0.05% formic acid and directly introduced into the source at a flow rate of 5 μl/min. Capillary entrance voltage was set to 2.7 kV, and dry gas temperature to 150°C. Voltages: Cone: 80 V, Rf lens: 40 V. MS profile [500 (10%), 1500 (60%), 2500 (20%), ramp 10%]. Scanning domain: m/z 1000-3000. Calibration was performed with orthophosphoric clusters. Continuum spectra exhibiting multicharged ions were transformed into molecular mass envelops using the MaxEnt 1 software.
Electromobility shift assay
A set of DNA probes covering the predicted Clr binding palindrome were obtained by annealing two complementary oligonucleotides. The annealing reactions were performed in water with 25 μM strand + (WTN8+ or MN8+ (see Additional file 10)) and 25 μM strand–(WTN8-or MN8-(see Additional file 10)) for each probe in a total reaction volume of 100 μl. Mixes were incubated at 95°C during 5 min following by slow cooling to 25°C. 175 nM double-strands probes were end labelled using 20 μCi of [ATPγ-32P] and 10 U of T4 polynucleotide kinase (Promega). Probes (1.75 nM each) were incubated in binding buffer (10 mM Tris [pH 8.0], 1 mM EDTA, 1 mM DTT, 10 μg/ml bovine serum albumin, 100 mM KCl) containing 50 μg/ml poly(2′-deoxyinosinic-2′-deoxycytidylic acid) (Sigma) and 10% glycerol for 30 min at room temperature with purified Clr and 3′, 5′cAMP or 2′, 3′cAMP added to the concentrations indicated in the figure legends in a final reaction volume of 15 μl. Samples were subjected to electrophoresis on a 10% polyacrylamide TBE 0.5 X gel containing 4% PEG-8000. Electrophoresis was conducted in TBE 0.5 X buffer at 80 V at room temperature. Gels were dried and analysed by autoradiography.
Plant assays and plant extracts preparation
Seeds of M. sativa cv. Europe were surface sterilized, germinated, and allowed to grow in 12-cm2 plates containing slanting nitrogen-free Fahraeus agar medium for 3 days at 22°C with day/night cycles of 16/8 h. The plants were inoculated with 2.103 bacteria per plant. Nodules were counted every day during 8 days then every 2 days until 35 days post-inoculation (dpi). At 35 dpi, shoots were collected and dried overnight at 65°C for weight measurements. Plant extracts were prepared as previously described .
S. meliloti strains carrying the pGD2178, pXLGD4 or pGD2179 plasmids were grown at 28°C in VGM. Overnight cultures were diluted to an OD600 of 0.1 in VGM and grown for an additional 2 h. 5 ml-cultures supplemented with 3′, 5′cAMP (2.5 mM or 5 mM), 2′, 3′cAMP (7.5 mM) or 5 mM 3′, 5′cGMP were grown for an additional 5 hours at 28°C. Overnight incubation was used for other potential inducers listed in Additional file 3. β-Galactosidase activities were measured (Miller units) using 1 ml of culture (or 0.1 ml for overnight cultures), as previously described .
Plants were inoculated with S. meliloti strains carrying the pGD2178 or the pGD2179 plasmid. Entire roots were collected 7 dpi or 14 dpi, fixed with 2% (vol/vol) glutaraldehyde solution for 1.5 h under vacuum, rinsed three times in Z buffer (0.1 M potassium phosphate buffer [pH 7.4], 1 mM MgSO4, and 10 mM KCl), and stained overnight at 28°C in Z buffer containing 0.08% 5-bromo-4-chloro-3-indolyl-D-galactoside (X-gal), 5 mM K3Fe(CN)6, and 5 mM K4Fe(CN)6. Nodules were harvested at 14 dpi, fixed with 2% (v/v) glutaraldehyde in Z buffer, and then sliced into 70 μm-thick longitudinal sections using a vibrating-blade microtome (VT1000S; Leica) before staining overnight at 28°C. Entire roots or nodule sections were observed under a light microscope.
Phosphodiesterase activity assays
Biochemical assays were performed in 50 mM Tris–HCl [pH 8], 5 mM β-Mercaptoethanol, 10 mM NaCl, 100 μM MnCl2, and 0 to 2.5 mM bis-P-nitrophenyl phosphate in a total volume of 50 μl. Reactions were initiated by the addition of 120 nM SpdA and the reaction was stopped after 10 min at 25°C by the addition of 10 μl of 200 mM NaOH. Release of p-nitrophenol was determined by measuring the absorbance at 405 nm. Cyclic NMP assays were performed in reaction mixtures containing 50 mM Tris–HCl [pH 8], 5 mM β-Mercaptoethanol, 10 mM NaCl, 10 mM cyclic nucleotides, 1 μM SpdA and 10 U calf intestine phosphatase (CIP) were incubated 10 min at 25°C, and were stopped by the addition of 1 ml Biomol Green Reagent (Enzo). Released of phosphate was determined by measuring the absorbance at 620 nm. The kinetic values were determined using the equation of v = Vmax [S]/(Km + [S]) where v, Vmax, Km and [S] represent the initial velocity, the maximum velocity, the Michaelis constant and the substrate concentration, respectively. The Kcat was calculated by dividing Vmax by the concentration of enzyme used in the reaction (Kcat = Vmax/[enzyme]).
3′, 5′cAMP affinity matrix was purchased from Sigma. 4.5 mM of purified Clr-GST was incubated in batch with 200 μl of 3′, 5′cAMP-agarose, previously equilibrated in buffer A (100 mM sodium phosphate buffer [pH 7], 50 mM NaCl, at 4°C during 30 min on a rotary mixer. After washing 7 times with 1 ml buffer A, bound protein was eluted by 30 min incubation in 1 ml buffer A supplemented with 30 mM 3′, 5′cAMP or 30 mM 2′, 3′cAMP at 4°C. Fractions were analysed by 12% SDS-PAGE.
Jones KM, Kobayashi H, Davies BW, Taga ME, Walker GC: How rhizobial symbionts invade plants: the Sinorhizobium-Medicago model. Nat Rev Microbiol. 2007, 5 (8): 619-633. 10.1038/nrmicro1705.
Murray JD: Invasion by invitation: rhizobial infection in legumes. Mol Plant Microbe Interact. 2011, 24 (6): 631-639. 10.1094/MPMI-08-10-0181.
Tian CF, Garnerone AM, Mathieu-Demazière C, Masson-Boivin C, Batut J: Plant-activated bacterial receptor adenylate cyclases modulate epidermal infection in the Sinorhizobium meliloti-Medicago symbiosis. Proc Natl Acad Sci USA. 2012, 109 (17): 6751-6756. 10.1073/pnas.1120260109.
He Y, Li N, Chen Y, Chen X, Bai J, Wu J, Xie J, Ying H: Cloning, expression, and characterization of an adenylate cyclase from Arthrobacter sp. CGMCC 3584. Appl Microbiol Biotechnol. 2012, 96 (4): 963-970. 10.1007/s00253-012-3890-x.
McDonough KA, Rodriguez A: The myriad roles of cyclic AMP in microbial pathogens: from signal to sword. Nat Rev Microbiol. 2012, 10 (1): 27-38.
Linder JU: Class III adenylyl cyclases: molecular mechanisms of catalysis and regulation. Cell Mol Life Sci. 2006, 63 (15): 1736-1751. 10.1007/s00018-006-6072-0.
Masson-Boivin C, Giraud E, Perret X, Batut J: Establishing nitrogen-fixing symbiosis with legumes: how many rhizobium recipes?. Trends Microbiol. 2009, 17 (10): 458-466. 10.1016/j.tim.2009.07.004.
Shenroy AR, Visweswariah SS: Class III nucleotide cyclases in bacteria and archaebacteria: lineage-specific expansion of adenylyl cyclases and a dearth of guanylyl cyclases. FEBS Lett. 2004, 561 (1–3): 11-21.
Kimura Y, Mishima Y, Nakano H, Takegawa K: An adenylyl cyclase, CyaA, of Myxococcus xanthus functions in signal transduction during osmotic stress. J Bacteriol. 2002, 184 (13): 3578-3585. 10.1128/JB.184.13.3578-3585.2002.
Kimura Y, Ohtani M, Takegawa K: An adenylyl cyclase, CyaB, acts as an osmosensor in Myxococcus xanthus. J Bacteriol. 2005, 187 (10): 3593-3598. 10.1128/JB.187.10.3593-3598.2005.
Agarwal N, Lamichhane G, Gupta R, Nolan S, Bishai WR: Cyclic AMP intoxication of macrophages by a Mycobacterium tuberculosis adenylate cyclase. Nature. 2009, 460 (7251): 98-102. 10.1038/nature08123.
Agarwal N, Bishai WR: cAMP signaling in Mycobacterium tuberculosis. Indian J Exp Biol. 2009, 47 (6): 393-400.
Topal H, Fulcher NB, Bitterman J, Salazar E, Buck J, Levin LR, Cann MJ, Wolfgang MC, Steegborn C: Crystal structure and regulation mechanisms of the CyaB adenylyl cyclase from the human pathogen Pseudomonas aeruginosa. J Mol Biol. 2012, 416 (2): 271-286. 10.1016/j.jmb.2011.12.045.
Hall RA, De Sordi L, Maccallum DM, Topal H, Eaton R, Bloor JW, Robinson GK, Levin LR, Buck J, Wang Y, et al: CO(2) acts as a signalling molecule in populations of the fungal pathogen Candida albicans. PLoS Pathog. 2010, 6 (11): e1001193-10.1371/journal.ppat.1001193.
Xu XL, Lee RT, Fang HM, Wang YM, Li R, Zou H, Zhu Y, Wang Y: Bacterial peptidoglycan triggers Candida albicans hyphal growth by directly activating the adenylyl cyclase Cyr1p. Cell Host Microbe. 2008, 4 (1): 28-39. 10.1016/j.chom.2008.05.014.
Capela D, Barloy-Hubler F, Gouzy J, Bothe G, Ampe F, Batut J, Boistard P, Becker A, Boutry M, Cadieu E, et al: Analysis of the chromosome sequence of the legume symbiont Sinorhizobium meliloti strain 1021. Proc Natl Acad Sci USA. 2001, 98 (17): 9877-9882. 10.1073/pnas.161294398.
Richter W: 3′,5′ Cyclic nucleotide phosphodiesterases class III: members, structure, and catalytic mechanism. Proteins. 2002, 46 (3): 278-286. 10.1002/prot.10049.
Shenoy AR, Visweswariah SS: New messages from old messengers: cAMP and mycobacteria. Trends Microbiol. 2006, 14 (12): 543-550. 10.1016/j.tim.2006.10.005.
Jackson EK: The 2′,3′-cAMP-adenosine pathway. Am J Physiol Renal Physiol. 2011, 301 (6): F1160-F1167. 10.1152/ajprenal.00450.2011.
Shenoy AR, Sreenath N, Podobnik M, Kovacevic M, Visweswariah SS: The Rv0805 gene from Mycobacterium tuberculosis encodes a 3′,5′-cyclic nucleotide phosphodiesterase: biochemical and mutational analysis. Biochemistry. 2005, 44 (48): 15695-15704. 10.1021/bi0512391.
Macfadyen LP, Ma C, Redfield RJ: A 3′,5′ cyclic AMP (cAMP) phosphodiesterase modulates cAMP levels and optimizes competence in Haemophilus influenzae Rd. J Bacteriol. 1998, 180 (17): 4401-4405.
Imamura R, Yamanaka K, Ogura T, Hiraga S, Fujita N, Ishihama A, Niki H: Identification of the cpdA gene encoding cyclic 3′,5′-adenosine monophosphate phosphodiesterase in Escherichia coli. J Biol Chem. 1996, 271 (41): 25423-25429. 10.1074/jbc.271.41.25423.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31 (13): 3497-3500. 10.1093/nar/gkg500.
Podobnik M, Tyagi R, Matange N, Dermol U, Gupta AK, Mattoo R, Seshadri K, Visweswariah SS: A mycobacterial cyclic AMP phosphodiesterase that moonlights as a modifier of cell wall permeability. J Biol Chem. 2009, 284 (47): 32846-32857. 10.1074/jbc.M109.049635.
Marx CJ, Lidstrom ME: Broad-host-range cre-lox system for antibiotic marker recycling in gram-negative bacteria. Biotechniques. 2002, 33 (5): 1062-1067.
Trülzsch K, Roggenkamp A, Pelludat C, Rakin A, Jacobi C, Heesemann J: Cloning and characterization of the gene encoding periplasmic 2′,3′-cyclic phosphodiesterase of Yersinia enterocolitica O:8. Microbiology. 2001, 147 (Pt 1): 203-213.
Sallet E, Roux B, Sauviac L, Jardinaud MF, Carrère S, Faraut T, De Carvalho-Niebel F, Gouzy J, Gamas P, Capela D, et al: Next-generation annotation of prokaryotic genomes with EuGene-P: application to Sinorhizobium meliloti 2011. DNA Res. 2013, 20 (4): 339-354. 10.1093/dnares/dst014.
Lawson CL, Swigon D, Murakami KS, Darst SA, Berman HM, Ebright RH: Catabolite activator protein: DNA binding and transcription activation. Curr Opin Struct Biol. 2004, 14 (1): 10-20. 10.1016/j.sbi.2004.01.012.
Fuchs EL, Brutinel ED, Klem ER, Fehr AR, Yahr TL, Wolfgang MC: In vitro and in vivo characterization of the Pseudomonas aeruginosa cyclic AMP (cAMP) phosphodiesterase CpdA, required for cAMP homeostasis and virulence factor regulation. J Bacteriol. 2010, 192 (11): 2779-2790. 10.1128/JB.00168-10.
Liu J, Burns DM, Beacham IR: Isolation and sequence analysis of the gene (cpdB) encoding periplasmic 2′,3′-cyclic phosphodiesterase. J Bacteriol. 1986, 165 (3): 1002-1010.
Keppetipola N, Shuman S: A phosphate-binding histidine of binuclear metallophosphodiesterase enzymes is a determinant of 2′,3′-cyclic nucleotide phosphodiesterase activity. J Biol Chem. 2008, 283 (45): 30942-30949. 10.1074/jbc.M805064200.
Kimura Y, Okazaki N, Takegawa K: Enzymatic characteristics of two novel Myxococcus xanthus enzymes, PdeA and PdeB, displaying 3′,5′- and 2′,3′-cAMP phosphodiesterase, and phosphatase activities. FEBS Lett. 2009, 583 (2): 443-448. 10.1016/j.febslet.2008.12.044.
Galperin MY, Bairoch A, Koonin EV: A superfamily of metalloenzymes unifies phosphopentomutase and cofactor-independent phosphoglycerate mutase with alkaline phosphatases and sulfatases. Protein Sci. 1998, 7 (8): 1829-1835. 10.1002/pro.5560070819.
Botha FC, Dennis DT: Isozymes of phosphoglyceromutase from the developing endosperm of Ricinus communis: isolation and kinetic properties. Arch Biochem Biophys. 1986, 245 (1): 96-103. 10.1016/0003-9861(86)90193-1.
Yakunin AF, Proudfoot M, Kuznetsova E, Savchenko A, Brown G, Arrowsmith CH, Edwards AM: The HD domain of the Escherichia coli tRNA nucleotidyltransferase has 2′,3′-cyclic phosphodiesterase, 2′-nucleotidase, and phosphatase activities. J Biol Chem. 2004, 279 (35): 36819-36827. 10.1074/jbc.M405120200.
Hantke K, Winkler K, Schultz JE: Escherichia coli exports cyclic AMP via TolC. J Bacteriol. 2011, 193 (5): 1086-1089. 10.1128/JB.01399-10.
Jackson EK, Ren J, Mi Z: Extracellular 2′,3′-cAMP is a source of adenosine. J Biol Chem. 2009, 284 (48): 33097-33106. 10.1074/jbc.M109.053876.
Vallenet D, Belda E, Calteau A, Cruveiller S, Engelen S, Lajus A, Le Fèvre F, Longin C, Mornico D, Roche D, et al: MicroScope–an integrated microbial resource for the curation and comparative analysis of genomic and metabolic data. Nucleic Acids Res. 2013, 41 (Database issue): D636-D647.
Capela D, Filipe C, Bobik C, Batut J, Bruand C: Sinorhizobium meliloti differentiation during symbiosis with alfalfa: a transcriptomic dissection. Mol Plant Microbe Interact. 2006, 19 (4): 363-372. 10.1094/MPMI-19-0363.
Arcus VL, McKenzie JL, Robson J, Cook GM: The PIN-domain ribonucleases and the prokaryotic VapBC toxin-antitoxin array. Protein Eng Des Sel. 2011, 24 (1-2): 33-40. 10.1093/protein/gzq081.
Min AB, Miallau L, Sawaya MR, Habel J, Cascio D, Eisenberg D: The crystal structure of the Rv0301-Rv0300 VapBC-3 toxin-antitoxin complex from M. tuberculosis reveals a Mg2+ ion in the active site and a putative RNA-binding site. Protein Sci. 2012, 21 (11): 1754-1767. 10.1002/pro.2161.
Jung K, Fried L, Behr S, Heermann R: Histidine kinases and response regulators in networks. Curr Opin Microbiol. 2012, 15 (2): 118-124. 10.1016/j.mib.2011.11.009.
Pesavento C, Hengge R: Bacterial nucleotide-based second messengers. Curr Opin Microbiol. 2009, 12 (2): 170-176. 10.1016/j.mib.2009.01.007.
Corrigan RM, Gründling A: Cyclic di-AMP: another second messenger enters the fray. Nat Rev Microbiol. 2013, 11 (8): 513-524. 10.1038/nrmicro3069.
Becker A, Bergès H, Krol E, Bruand C, Rüberg S, Capela D, Lauber E, Meilhoc E, Ampe F, De Bruijn FJ, et al: Global changes in gene expression in Sinorhizobium meliloti 1021 under microoxic and symbiotic conditions. Mol Plant Microbe Interact. 2004, 17 (3): 292-303. 10.1094/MPMI.2004.17.3.292.
Pandey SP, Minesinger BK, Kumar J, Walker GC: A highly conserved protein of unknown function in Sinorhizobium meliloti affects sRNA regulation similar to Hfq. Nucleic Acids Res. 2011, 39 (11): 4691-4708. 10.1093/nar/gkr060.
Miller JH: Molecular genetics experiments. 1972, Cold Spring Harbor, N.Y: Cold Spring Harbor Laboratory
We thank the Florimond-Desprez company (Cappelle en Perche, France) for generous gift of Medicago seeds. CMD was supported by a PhD fellowship from the French Ministère de l’Enseignement supérieur et de la Recherche. We thank Chang-Fu Tian (LIPM, State Key Laboratory of Agrobiotechnology, China Agricultural University) for his contribution to the early stages of this study and Ludovic Cottret (Plateforme Bioinformatique LIPM-SPE) for his help with phylogenetic analyses.
This work was funded in part by the ANR “RhizocAMP” (ANR-10-BLAN-1719) and the Pôle de Compétitivité “Agrimip Innovation Sud Ouest”. This work is part of the “Laboratoire d’Excellence” (LABEX) entitled TULIP (ANR-10-LABX-41).
The authors declared they have no competing interests.
CMD, VP and AMG collected and analysed data. AMG and JB directed the work. JB and CMB wrote and revised the manuscript. All authors have read and approved the final version of this manuscript.
Anne-Marie Garnerone and Jacques Batut contributed equally to this work.