ORIGINAL RESEARCH article
Distinct and Specific Role of NlpC/P60 Endopeptidases LytA and LytB in Cell Elongation and Division of Lactobacillus plantarum
- 1Louvain Institute of Biomolecular Science and Technology, Université Catholique de Louvain, Louvain-La-Neuve, Belgium
- 2Micalis Institute, INRA, AgroParisTech, Université Paris-Saclay, Jouy-en-Josas, France
Peptidoglycan (PG) is an essential lattice of the bacterial cell wall that needs to be continuously remodeled to allow growth. This task is ensured by the concerted action of PG synthases that insert new material in the pre-existing structure and PG hydrolases (PGHs) that cleave the PG meshwork at critical sites for its processing. Contrasting with Bacillus subtilis that contains more than 35 PGHs, Lactobacillus plantarum is a non-sporulating rod-shaped bacterium that is predicted to possess a minimal set of 12 PGHs. Their role in morphogenesis and cell cycle remains mostly unexplored, except for the involvement of the glucosaminidase Acm2 in cell separation and the NlpC/P60 D, L-endopeptidase LytA in cell shape maintenance. Besides LytA, L. plantarum encodes three additional NlpC/P60 endopeptidases (i.e., LytB, LytC and LytD). The in silico analysis of these four endopeptidases suggests that they could have redundant functions based on their modular organization, forming two pairs of paralogous enzymes. In this work, we investigate the role of each Lyt endopeptidase in cell morphogenesis in order to evaluate their distinct or redundant functions, and eventually their synthetic lethality. We show that the paralogous LytC and LytD enzymes are not required for cell shape maintenance, which may indicate an accessory role such as in PG recycling. In contrast, LytA and LytB appear to be key players of the cell cycle. We show here that LytA is required for cell elongation while LytB is involved in the spatio-temporal regulation of cell division. In addition, both PGHs are involved in the proper positioning of the division site. The absence of LytA activity is responsible for the asymmetrical positioning of septa in round cells while the lack of LytB results in a lateral misplacement of division planes in rod-shaped cells. Finally, we show that the co-inactivation of LytA and LytB is synthetically affecting cell growth, which confirms the key roles played by both enzymes in PG remodeling during the cell cycle of L. plantarum. Based on the large distribution of NlpC/P60 endopeptidases in low-GC Gram-positive bacteria, these enzymes are attractive targets for the discovery of novel antimicrobial compounds.
The cell wall is a rigid structure that protects bacteria against external and internal pressures while giving them a proper shape (Delcour et al., 1999). A major component of the cell wall is the peptidoglycan (PG). The PG is a polymer composed of glycan strands linked together by peptide side chains to form a meshwork. The composition of the glycan strands is the same in all bacterial species: N-acetyl-muramic acid (MurNAc) alternating with N-acetyl-glucosamine (GlcNAc), which are connected by β-1,4 linkage (Vollmer, 2008). The composition of the peptide side chain is more flexible and can vary between bacterial species (Vollmer et al., 2008a). In Lactobacillus plantarum, the composition of the stem peptide is L-Ala, D-Glu, meso-diaminopimelate (meso-DAP), D-Ala, and D-lactate (Deghorain et al., 2007; Kleerebezem et al., 2010). The PG of L. plantarum is decorated with additional elements such as wall teichoic acids (WTA), O-acetylation of MurNAc (39%) and GlcNAc (9%) (Bernard et al., 2011a), and amidation of D-Glu (100%) and meso-DAP (94%) (Bernard et al., 2011b).
Bacteria have to continuously adapt the composition and organization of the PG in order to grow and divide properly. This task is performed by the concerted action between PG synthases and PG hydrolases (PGHs). While PG synthases mediate two main enzymatic reactions, i.e., transglycosylation and transpeptidation, PGHs have a broad range of enzymatic activities with respect to their role in PG processing (Vollmer et al., 2008b). They are separated into different families according to their PG cleavage sites (Vollmer et al., 2008b). Three families of enzymes cleave the glycan strand: glucosaminidases hydrolyze the sugar bond between GlcNAc and MurNAc (Litzinger et al., 2010); muramidases and lytic transglycosylases cleave between MurNAc and GlcNAc (Höltje et al., 1975; Barrett et al., 1984). Amidases cleave the amide bond between the glycan strand and the peptide side chain (Vollmer et al., 2008b). The last family is composed of endopeptidases (Vollmer et al., 2008b) and carboxypeptidases (Sauvage et al., 2008) that hydrolyze specific links in the peptide side chain or the interpeptide bridge.
In the last family, the NlpC/P60 D,L-endopeptidases deserves a special interest in Gram-positive bacteria since they were shown to have a major role in morphogenesis and cell cycle of Bacillus subtilis. In this species, seven NlpC/P60 endopeptidases were identified (Figure 1A). With the exception of PgdS that is strictly dedicated to the hydrolysis of poly-γ-glutamate, all the other D,L-endopeptidases were shown to hydrolyze PG (i.e., LytE, LytF, CwlS, CwlO, and CwlT) or PG fragments for their potential recycling (i.e., YkfC) (Schmidt et al., 2001; Fukushima et al., 2018). All of them were shown to display a γ-D-Glu-meso-DAP muropeptidase activity (Ishikawa et al., 1998; Margot et al., 1999; Schmidt et al., 2001; Yamaguchi et al., 2004; Fukushima et al., 2006, 2007). With the exception of the bifunctional CwlT enzyme that is involved in the conjugation of ICEBs1 (Fukushima et al., 2007), LytE, LytF, CwlS, and CwlO are modular enzymes implicated in morphogenesis (Hashimoto et al., 2012). CwlO and LytE, whose co-inactivation is synthetically lethal, are required for cell elongation (Hashimoto et al., 2012). However, they perform specific roles and they are differentially controlled by players of the elongation machinery (Domínguez-Cuevas et al., 2013; Meisner et al., 2013). Inactivation of CwlO leads to slightly bent and wider cells than the wild type while inactivation of LytE leads to slightly longer and thinner cells (Domínguez-Cuevas et al., 2013; Meisner et al., 2013). Besides its role in cell elongation, LytE was also reported to play a role in cell separation (Carballido-López et al., 2006). In addition, CwlO, which contains two coiled-coil domains, is activated by the membrane protein complex FtsEX (Domínguez-Cuevas et al., 2013; Meisner et al., 2013), while LytE, which contains three LysM PG-binding domains, was proposed to be guided by the actin-like cytoskeleton protein MreBH (Carballido-López et al., 2006; Domínguez-Cuevas et al., 2013; Meisner et al., 2013). Concerning the two last D,L-endopeptidases, LytF and CwlS, which contain five and four LysM domains, respectively, they were shown to be strictly implicated in the cell separation process (Yamamoto et al., 2003; Fukushima et al., 2006).
Figure 1. In silico analysis of NlpC/P60 endopeptidases of B. subtilis and L. plantarum. Schematic representation of LytE, LytF, CwlS, CwlO, CwlT and PgdS from B. subtilis 168 (A) and LytA, LytB, LytC, and LytD from L. plantarum WCFS1 (B) Protein size (aa) corresponds to the precursor with its exportation signal-peptide. SH3b-like domain structures were predicted by Phyre 2.0 (www.sbg.bio.ic.ac.uk/phyre2/). (C) Alignment of conserved regions of NlpC/P60 domains that contain the essential residues for catalysis (Cys, His, Asp; boxed). NlpC/P60 endopeptidases of B. subtilis and L. plantarum are indicated in black and red, respectively. (D) Conserved regions between L. plantarum LytC, B. subtilis YkfC (YkfCBsu), and B. cereus YkfC (YkfCBce) that are involved in the selection of PG stem peptides ending with L-Ala. Key residues involved in substrate selection (Tyr and Asp) are boxed. Alignments in (C,D) were obtained with PRALINE (http://www.ibi.vu.nl/programs/pralinewww/).
Lactobacillus plantarum is a commensal Gram-positive bacterium of the human gastrointestinal tract, which displays a B. subtilis–like rod shape (Kleerebezem et al., 2010). As a non-sporulating bacterium, L. plantarum only harbors 12 predicted PGHs, while Bacillus subtilis genome encodes at least 35 PGHs (Smith et al., 2000; Rolain et al., 2012). The predicted autolysome of L. plantarum contains 2 glucosaminidases (Acm1 and Acm2), 2 muramidases (Lys1 and Lys2), 3 lytic transglycosidases (MltA, MltB, MltC), 1 amidase (LytH), and 4 D,L-endopeptidases of the NlpC/P60 family (LytA, LytB, LytC, and LytD) (Rolain et al., 2012). In a previous work, we have investigated the functional role of these PGHs through systematic gene inactivation (Rolain et al., 2012). Knockout mutants of 9 out of the 12 predicted PGHs were successfully obtained (Rolain et al., 2012). Two PGHs, Acm2 and LytA, were shown to be differentially involved in the cell cycle. The glucosaminidase Acm2 is responsible of cell separation at a post-divisional stage, while the predicted D,L-endopeptidase LytA was identified as a morphogenic PGH given that its loss resulted in round cells (Rolain et al., 2012). Besides LytA, the function of the three other putative D,L-endopeptidases remain mostly unexplored (Figure 1B; Rolain et al., 2012). By comparison with the modular organization of B. subtilis D,L-endopeptidases, LytC and LytD display a similar organization as YkfC with SH3 PG-binding domains, while LytA and LytB possess interesting features such as the absence of coiled-coil domains, a limited number of LysM domains, and the presence of a glycosylated domain rich in alanine, serine, threonine (AST) (Figure 1B). Compared to B. subtilis, these specificities suggest that they could have different functions and/or be differently regulated.
The aim of this work is to dissect the role of the four Lyt endopeptidases of L. plantarum. To this end, cell morphogenesis was carefully examined in simple mutants of each individual lyt gene and double mutants of paralogous genes (i.e., lytA lytB and lytC lytD). While LytC and LytD did not seem to contribute to cell morphogenesis, LytA and LytB are key morphogenic PGHs. LytA is required for cell elongation, while LytB plays a role in the timing of cell division. Both PGHs are needed for the correct positioning of the division site and display a synthetic growth defect when co-inactivated.
Materials and Methods
Strains, Plasmids, and Growth Conditions
Strains and plasmids used in this work are listed in Table 1. Plasmid constructions were performed in Escherichia coli (strains AbleK and DH5α). Functional study of Lyt enzymes was performed in L. plantarum NZ7100. E. coli strains were grown at 37°C with shaking in LB (Lysogeny Broth) medium and L. plantarum strains were grown at 30°C in MRS broth (Difco). When appropriate, antibiotics were added to the medium. Chloramphenicol was used at the concentration of 10 μg ml−1 for E. coli and L. plantarum; and erythromycin was used at the concentration of 250 μg ml−1 for E. coli and 12.5 μg ml−1 for L. plantarum.
Nisin and ComS Induction
Nisin A (Sigma Aldrich) was used to induce PnisA in the conditional mutants. Conditional mutant strains were cultured overnight in presence of nisin 25 ng ml−1. Then, after dilution at an OD600 of 0.1, strains were either not induced or induced with nisin (1 or 25 ng ml−1).
The ComS-inducible system was used for complementation assays and localization of the FtsZ ring (fusion FtsZ-GFP+). ComS from S. thermophilus (LPYFAGCL) was synthesized by Peptide 2.0 (Chantilly, VA, USA). The dehydrated peptide was suspended in sterile milliQ water at a stock concentration of 100 μM. Complementation strains were diluted at an OD600 of 0.1, induced 2 h later with 8 μM of ComS, and observed between 1 h 30 min and 6 h after induction. Strains producing FtsZ-GFP+ were diluted at an OD600 of 0.1, induced 2 h later as reported above, and shaken for improving GFP+ maturation until their observation 2 h later.
Bacteria were cultured overnight in MRS supplemented with antibiotics and inducers when needed, diluted at an OD600 of 0.05 in the same medium, and separated in 96-wells plates. The growth was monitored in a multi-plate reader Infinity Pro-200 (Tecan) every 10 min at 600 nm during 12 h.
Cells were collected in exponential phase from MRS cultures (with antibiotics and inducers when needed) and resuspended in PBS buffer. Bacteria were observed on agarose pads composed of 1% agarose PBS buffer for static observations and of 1% agarose MRS for time-lapse observations. Cellular membranes were stained with FM4-64 (Life Technologies) as reported before (Andre et al., 2011). Images were obtained using an Axio I inverted microscope (Zeiss) equipped with an α Plan-Apochromat objective (100 × /1.46 Oil DIC M27) (Zeiss), a HXP 120 C lighting unit (Zeiss) and C10600 ORCA-R2 camera (Hamamatsu). The fluorescence of FM4-64 and FtsZ-GFP+ was respectively detected with filter sets Cy3 (43 HE) and GFP (38 HE), displaying bandpass excitation (nm): 550/25 (Cy3) or 470/40 (GFP) and bandpass emission: 605/70 (Cy3) or 525/50 (GFP) (Zeiss). Images were analyzed using Axiovision 4.8 (Zeiss), MicrobeTracker (Sliusarenko et al., 2011), or MicrobeJ (Ducret et al., 2016).
Profiling of Muropeptides
PG from L. plantarum strains was prepared by treating a bacterial pellet with SDS, nucleases, and proteases solutions in order to eliminate all the cellular components except PG, according to a protocol previously described (Courtin et al., 2006). This protocol was slightly modified by applying DNase (50 μg ml−1) and RNase (50 μg ml−1) treatments before hydrofluoric acid extraction. PG was digested with mutanolysin from Streptomyces globisporus (Sigma-Aldrich). The resulting muropeptides were analyzed by RP-HPLC as previously reported (Courtin et al., 2006). Muropeptides were identified according to their retention times by comparison to the previously published reference chromatogram for L. plantarum PG (Bernard et al., 2011a). In addition, disaccharide-dipeptide (Di) purified from Lactobacillus casei (Regulski et al., 2012) was used as standard. The relative abundance (in %) of Di (with or without O-acetylation) was calculated as the ratio of the areas of the two peaks over the sum of the areas of all the identified peaks on the chromatogram.
DNA Manipulations and Transformation
Classical methods of molecular biology were used as previously described (Sambrook et al., 1989). Preparation of electro-competent cells and transformation of E. coli and L. plantarum were performed as previously reported (Dower et al., 1988; Josson et al., 1989). The Phusion High Fidelity polymerase (NEB) was used for amplification by PCR of inserts used for plasmid constructions. PCR amplifications for validation were performed with the GoTaq polymerase (Promega) in a GeneAmp PCR system (Applied Biosystem). Primers were synthesized by Eurogentec (Belgium) and are listed in Supplementary Table 1.
General Strategy for the Construction of Plasmids and Mutant Strains
All expression and disruption plasmids were constructed in E. coli before their transfer by electroporation in L. plantarum. The correct construction of the vectors was validated by PCR (validation primers in Supplementary Table 1) and enzymatic digestion. DNA sequencing was performed to validate the final product (validation primers in Supplementary Table 1). Suicide vectors for gene inactivation were recombined by simple homologous recombination into the chromosome of L. plantarum as previously described (Palumbo et al., 2004). The conditional mutants contained a 3′-end truncated copy of the targeted gene (inactive, truncation of the NlpC/P60 domain for LytA and LytC) and a second intact copy under the control of the nisin-inducible promoter. The validation of mutant strains was performed by PCR amplification of regions flanking the site of vector integration and sequencing of the PCR product (validation primers in Supplementary Table 1).
Construction details of plasmids for conditional/disruption mutants, for complementation, for ComS-inducible expression, and for expression of truncated/hybrid proteins and FtsZ-GFP+ fusion are presented in Supplementary Text.
In silico Analysis of NlpC/P60 Endopeptidases From L. plantarum
The four D,L-endopeptidases of L. plantarum can be separated in two groups based on their modular organization (Figure 1B). The first group contains LytA (Lp_3421, 370 aa) and LytB (Lp_2162, 496 aa). In their mature form, these two enzymes contain LysM domain(s) at the N-terminus; a low complexity AST central domain; and a catalytic domain of the NlpC/P60 family (pfam00877) at the C-terminus. The second group is composed of LytC (Lp_2520, 297 aa) and LytD (Lp_1242, 243 aa). These two proteins contain putative SH3b domains in their N-terminal part that were identified by a structural prediction (Phyre 2.0) and an NlpC/P60 catalytic domain at their C-terminus.
All members of the NlpC/P60 family contain three conserved residues involved in catalysis: Cys, His and a polar residue (His, Asp or Gln) (Anantharaman and Aravind, 2003). The alignment of the NlpC/P60 domains of the four L. plantarum Lyt enzymes and the 7 B. subtilis D,L-endopeptidases showed that Cys and His residues are fully conserved, as well as the third polar residue that is an Asp in this case (Figure 1C). The 3D structure prediction of the four Lyt catalytic domains suggests that the three conserved catalytic residues are in close proximity (Supplementary Figure 1A), consistent with their contribution to catalysis. Interestingly, LytC and LytD displayed a modular organization, which is similar to YkfC of B. subtilis and Bacillus cereus (Xu et al., 2010, 2015). Although YkfC of B. subtilis has been proposed to be intracellular due to the absence of a predicted signal-sequence, the B. cereus enzyme is predicted to be extracellular (Xu et al., 2010), as for LytC and LytD. YkfC is a γ-D-Glu-meso-DAP endopeptidase specific for PG peptides with free N-terminal L-Ala, suggesting a role in PG peptide recycling (Schmidt et al., 2001; Xu et al., 2010). Some NlpC/P60 endopeptidases were identified as recycling enzymes based on the presence of two conserved residues: an aspartate in their active site and a tyrosine at the junction between catalytic and SH3 PG-binding domains, which both restrict the access of the catalytic site to short PG peptides (Xu et al., 2015). Interestingly, LytC is the sole Lyt enzyme that contains these two conserved residues (Figure 1C). We also examined the phylogenetic proximity of NlpC/P60 domains between the D,L-endopeptidases of L. plantarum and B. subtilis (Supplementary Figure 1B). Interestingly, LytC clustered with YkfC while LytD clustered with the conjugation-related CwlT autolysin. In addition, the NlpC/P60 domains of LytA and LytB form a separate group, while LytE, LytF, CwlS, and CwlO involved in morphogenesis cluster altogether (Supplementary Figure 1B).
Finally, the accessory domains of LytA and LytB display interesting features compared to the morphogenic D,L-endopeptidases of B. subtilis. Both enzymes contain a central AST domain of a different size (186 and 250 aa, respectively), which is predicted to be unstructured. Both domains were shown to be O-glycosylated, suggesting a regulation mechanism of the enzyme activity or stability (Fredriksen et al., 2013; Rolain et al., 2013). LytA and LytB also contain a lower number of LysM domains (i.e., 1 or 2) compared to the morphogenic endopeptidases of B. subtilis (i.e., 3–5) (Figures 1A,B). The LysM domain is a well-established PG-binding domain in PGHs that recognizes the GlcNAc – X – GlcNAc motif of PG glycan strands (Visweswaran et al., 2011; Mesnage et al., 2014). Interestingly, the LysM domains from LytA and LytB have different isoelectrical points (IP). The single LysM of LytA is basic (IP of 8.34) as found for LysM domains of B. subtilis D,L-endopeptidases (IP comprised between 9.22 and 10.14) while the two LysM of LytB are in the acidic range (IPs of 5.43 and 5.18). As the ability of LysM domains to bind PG is linked to their IP (Visweswaran et al., 2011), this could indicate that LytA and LytB bind PG under different pH conditions.
LytC and LytD Are Not Involved in Morphogenesis
In the first part of this work, we investigated the role of the putative D,L-endopeptidases LytC and LytD. A conditional mutant containing a chromosomal copy of the lytC gene under the control of the nisin-inducible promoter (PnisA-lytC fusion) was constructed. Growth and cell morphology of the mutant strain was monitored in the absence of the nisin inducer (N0). The mutant strain did not show any obvious difference compared to the WT (strain NZ7100) (Supplementary Figure 2). Similar results were obtained for a mutant in which the chromosomal copy of lytC was disrupted by single crossover with a suicide plasmid carrying an internal fragment of the gene (data not shown). The stable deletion of lytD (ΔlytD) has previously been obtained (Rolain et al., 2012). The phenotype of the mutant strain was reinvestigated regarding growth and cell morphology. As previously reported, it was indistinguishable from the WT (Supplementary Figure 2; Rolain et al., 2012). Since LytC and LytD may have redundant functions due to their similar organization (Figure 1B), a double lytC lytD mutant was also constructed. Again, this mutant grown in absence of nisin behaved as the WT (Supplementary Figure 2).
Altogether, these data indicate that LytC and LytD are not involved in L. plantarum morphogenesis under the tested conditions.
LytA Is a Major Peptidoglycan Hydrolase in Cell Morphogenesis
We previously reported the construction of a stable lytA null mutant (ΔlytA) (Rolain et al., 2012). The cell morphology of this mutant was strongly affected, with round and aggregated cells instead of rod-shaped and well-separated cells observed for wild type (WT) cells (Figure 2A). In addition, this mutant displayed a growth defect and cannot be transformed to perform further analyses such as complementation studies (Rolain et al., 2012). For these reasons, we decided to construct a conditional lytA mutant by placing the chromosomal copy of lytA under the control of PnisA (PnisA-lytA fusion). In absence of nisin (N0), cell morphology and growth of this conditional mutant strain was similar to the stable ΔlytA mutant (Figures 2A,B). Under nisin induction (25 ng ml−1, N25), the strain recovered a WT phenotype regarding growth and morphology, except for a slightly larger diameter (Figures 2A,B).
Figure 2. Effect of LytA deficiency on cell morphology, growth, and PG composition. (A) Images of L. plantarum cells obtained by phase contrast (PC) microscopy and epifluorescence microscopy for membrane labeling with FM4-64. Left panel, WT and ΔlytA mutant; middle panel, conditional PnisA-lytA mutant without (N0) or with nisin 25 ng ml−1 (N25); right panel, complementation of PnisA-lytA mutant with LytA (Pshp0064-lytA; + LytA) and a catalytic mutant of LytA (Pshp0064-lytA*; + LytA*), grown without nisin (N0) in presence of ComS (8 μM, C8). Cells were collected in exponential phase from MRS cultures (with chloramphenicol when needed) and observed on agarose pads after suspension in PBS. Similar observations were obtained from at least 3 independent experiments. The scale bar is 2 μm. (B) Growth curves of WT, ΔlytA mutant, and PnisA-lytA mutant (N0 and N25) in MRS medium. Curves were generated from triplicates (mean values + standard deviations). (C) Percentage of disaccharides-dipeptides without and with O-acetylation (Di and Di+Ac, respectively) in the PG of WT and PnisA-lytA mutant (N0 and N25) after mutanolysin digestion. Mean values of three independent extractions ± standard deviations. Significance with respect to the WT (R, reference) is based on Student's t-test. *P < 0.05 and ***P < 0.001, respectively.
The conditional mutant strain was complemented by a copy of lytA expressed under the control of the ComS-inducible system (Pshp0064-lytA fusion) carried by a low copy number plasmid (see Supplementary Text for details on the ComS-inducible system). In the presence of the inducer peptide ComS (8 μM; C8), the LytA-depleted strain recovered the WT rod-shaped morphology (Figure 2A). In addition, we showed that LytA activity is crucial since complementation with a catalytic inactive mutant of LytA [Cys284 to Ala (Figure 1C); named LytA*] was unable to restore the WT rod-shaped morphology (Figure 2A).
Previous analyses of PG extracted from the ΔlytA mutant suggested that LytA was cleaving the bond between D-Glu and meso-DAP in the peptide side chains of the PG network (Rolain et al., 2012). To further correlate the observed morphological phenotypes with LytA activity, analyses of PG composition were performed in triplicates for the conditional mutant strain (PnisA-lytA) in presence or absence of nisin (Figure 2C, Supplementary Figure 3, and Supplementary Table 2). PG analysis of the WT strain showed that the pool of disaccharides with two amino acids [GlcNAc-MurNAc-L-Ala-D-Glu with and without O-acetylation, (Di and Di+Ac, respectively)], which are potential products of D,L-endopeptidase activity on PG stem peptides, represents 1.83% of total soluble muropeptides. For the nisin-depleted strain (N0), the combined amount of Di and Di+Ac was 5-fold lower (0.37%), consistent with a lower production of LytA. When the conditional mutant strain was cultured in the presence of nisin (25 ng ml−1, N25), the combined amount of Di and Di+Ac was 3-fold higher (1.13%) than in non-induced conditions, reaching a level close to the WT but remaining slightly lower. These data strongly suggest that LytA is a γ-D-Glu-meso-DAP endopeptidase that cleaves the bond between D-Glu and meso-DAP in the PG meshwork, and that its absence or decreased activity significantly alter the PG structure.
Altogether, these results demonstrate that LytA is a major PGH of L. plantarum morphogenesis and that its endopeptidase activity, more than its physical presence in a PG biosynthetic complex, is required for proper growth and cell cycle progression.
LytA Is Required for Cell Elongation and Septum Positioning
We previously reported that the morphology of the ΔlytA mutant was strongly affected with the presence of aggregated round cells of variable diameters and alterations of PG thickness at division sites (Rolain et al., 2012). Here, we used the conditional lytA mutant to understand the impact of a progressive depletion of LytA on cell morphogenesis. Cells were first stained with FM4-64 to observe membranes and septum position (Figure 3A) or labeled with an FtsZ-GFP+ fusion to visualize the divisional Z rings (Figure 3B). Figure 3 shows the degenerative process of the conditional lytA mutant taken at different time points after nisin removal. In early steps (0–2 h), cells remained similar to WT, forming chains of rod-shaped bacteria with division planes and Z rings at mid-cell. After 4 h of depletion, cells formed short chains of 2–3 deformed cells. Dividing cells started to bend and septa seemed to be placed in an abnormal orientation (Figure 3, Supplementary Figure 4A). After prolonged depletion (6 h), a large proportion of cells were round and aggregated. Membrane labeling appeared brighter and less homogenous, which may indicate that membrane biogenesis is altered [Figure 3 (6 h), Supplementary Figure 5A]. In round cells, FtsZ does not seem to be correctly localized, forming dots or arcs at the cell periphery (Figure 3, 6 h). These observations show that LytA is required for cell elongation. In addition, population analysis showed that more than 70% of dividing cells (n > 300) have an altered positioning of the division plane after 6 h of LytA depletion (deviation >10% from the median position), while only ~10% of WT cells have slightly misplaced septa (a sample of LytA-depleted cells with misplaced septa is shown in Supplementary Figures 4B, 5A). These asymmetrical division events likely account for the round cells of different diameters (Supplementary Figure 5B) that were previously reported for the stable ΔlytA mutant (Rolain et al., 2012). Mispositioning of division sites might be due to an indirect effect resulting from the lack of LytA-mediated PG remodeling activity, or from the absence of elongation.
Figure 3. Effect of progressive LytA depletion on cell morphology and division site positioning. (A) PnisA-lytA mutant cells observed during nisin depletion (0, 2, 4, 6 h) by phase contrast (PC) microscopy and FM4-64 staining. The scale bar is 2 μm. (B) PnisA-lytA mutant cells expressing an FtsZ-GFP+ fusion (Pshp0064-ftsZ-fgp+; + FtsZ-GFP+) observed during nisin depletion (0, 2, 4, 6 h) by phase contrast (PC) microscopy and epifluorescence microscopy for the localization of FtsZ-GFP+. White arrows indicate Z-rings in normal rod-shaped cells. Bacteria were cultured in MRS with erythromycin and chloramphenicol, and induced with ComS (8 μM, C8). Cultures were moderately shaken after ComS induction. The scale bar is 2 μm. For (A,B), similar observations were obtained from at least two independent experiments.
Morphological alterations during the cell cycle of lytA-depleted cells were also followed by time-lapse microscopy. The conditional lytA mutant was cultured overnight in the presence of nisin and then diluted on an MRS-containing agarose pad in the absence of the inducer. The cell cycle of the LytA-depleted strain was compared to the cell cycle of the WT strain (Figure 4, Supplementary Movies 1, 2). WT cells were shown to elongate until they reached twice their original length and then divided into two daughter cells of equal length (Figure 4A). For the LytA-depleted strain, cells started as rod-shaped bacteria and progressively reached the round and aggregated morphology after 2–3 generations (Figure 4B, left panel). Before the first division, a short elongation phase was observed but cells never doubled their length. Then, cells started to bend and division occurred like a break in the middle of the cell (V-shaped dividing cells). At the next generation, elongation is inhibited and new daughter cells were generated by asymmetrical division. This resulted in the formation of round cells of different diameters, including mini cells. At the end of the depletion process, only round and aggregated cells were observed due to absence of elongation, misplacement of division sites, and the lack of cell separation.
Figure 4. Comparison of the cell cycle of WT and LytA-depleted strain. (A) Time-lapse microscopy of the WT strain. (B) Time-lapse microscopy of the PnisA-lytA mutant strain without nisin (N0). The cell cycle of both strains is schematically represented on the right. Bacteria were grown on MRS-containing agarose pads at 30°C. The scale bar is 2 μm. For (A,B) one representative time lapse from at least two independent experiments.
To get further insight into the role of LytA during the cell cycle, we chose to compare the LytA-depleted phenotype with that obtained from the conditional inactivation of the mreB1CD locus whose homolog from B. subtilis has been shown to be specifically involved in cell elongation (Formstone and Errington, 2005). To this end, the chromosomal copy of L. plantarum mreB1CD operon was placed under the control of the nisin-inducible promoter (PnisA-mreB1CD fusion). The phenotype of the resulting MerB1CD-depleted cells was observed in time lapse (Supplementary Figure 6 and Supplementary Movie 3). From these observations, it appeared that cells became rapidly unable to elongate, giving rise to small spherical cells after one or two division cycles. This phenotype is not exactly the same as that observed for the LytA-depleted mutant in which impairment of elongation was accompanied with a defect in septum positioning. Thus, the data suggest that the lack of elongation for the lytA mutant is not, by itself, responsible for the mispositioning of the septum.
Altogether, the results demonstrate that LytA is a crucial player of L. plantarum cell cycle required for the elongation process. In addition, its absence has probably an indirect effect on the correct placement of the division site.
LytB Plays a Key Role in Septum Maturation and Timing of Division
We previously constructed a stable deletion mutant of lytB (ΔlytB) (Rolain et al., 2012). Growth and overall cell morphology of the mutant were initially reported to be similar to the WT (Figure 5A) (Rolain et al., 2012). However, a more thorough analysis of the cell size revealed that the mutant has a mean cell length significantly longer than the WT [3.5 ± 0.85 vs. 2.8 ± 0.67 μm, respectively; mean values ± standard deviations, n > 600 cells from triplicates, Kolmogorov–Smirnov (KS) test, P < 0.001], with a larger proportion of cells exceeding 4 μm (6-fold) and a lower proportion of cells shorter than 2 μm (16-fold) (Figure 5B). For the cell diameter, it remained unchanged in comparison to the WT (0.85 ± 0.12 vs. 0.82 ± 0.13 μm). The ΔlytB mutant was also complemented by a lytB copy under the control of the ComS-inducible system (Pshp0064-lytB fusion). In the presence of ComS, the mean cell length of the complemented strain was significantly reduced compared to growth conditions without ComS (3.0 ± 0.80 vs. 3.6 ± 0.70, respectively; n > 900 cells from triplicates, KS test, P < 0.001). However, the distribution of length frequencies of the complemented strain compared to the WT showed a partial complementation (Figure 5B), suggesting that the amount of LytB is inadequate for a full reversion of the phenotype.
Figure 5. Effect of LytB deficiency and its complementation on cell morphology and division site positioning. (A) Images of cells of WT, ΔlytB mutant, and ΔlytB mutant complemented with LytB (Pshp0064-lytB; + LytB) obtained by phase contrast (PC) microscopy. Yellow arrowheads show asymmetrical divisions in ΔlytB mutant cells. Bacteria were grown in MRS with erythromycin and ComS (8 μM, C8) when needed. Scale bar is 2 μm. (B) Cell length (μm) of WT, ΔlytB mutant and ΔlytB mutant + LytB without (C0) and with ComS (C8) measured in exponential growth phase. Cells were cultured in MRS with erythromycin and ComS when needed. (C) Number of septa (0 to 3; bars) related to cell length (lanes; mean values ± standard deviations) in WT and ΔlytB mutant cells. (D) Relative septum deviation (% of deviation from the median position of cells stained with FM4-64) in mono-septal cells of WT and ΔlytB mutant. For (B–D) measures were obtained from triplicates by using MicrobeJ with n > 500 cells.
The observed cell length increase in the LytB-deficient strain was correlated with septum positioning using FM4-64 labeling (Figure 5C, Supplementary Figure 7). LytB-deficient cells that are longer than 4 μm contained 2 or 3 septa (Figure 5C). The presence of two septa instead of three in more than 50% of long cells suggests that the formation of septa in daughter cells was not synchronized, while synchronization was observed in the WT (Supplementary Figure 7). In addition, lateral septa in WT daughter cells appeared at the end of the constriction of the median septum, while the maturation of the mid-cell septum seemed delayed in LytB-deficient cells (Supplementary Figure 7). This delayed maturation in dividing cells was confirmed by examining FtsZ positioning (FtsZ-GFP+). In WT, Z rings localized in daughter cells at the end of the constriction of the septum of the mother cell while the migration of Z rings in LytB-deficient cells took place at an earlier stage when maturation of the median septum was largely incomplete (Supplementary Figure 8). We also examined the lateral positioning of septa in dividing cells, which was affected in LytB-deficient cells compared to WT cells (Figure 5D, Supplementary Figure 7). In mono-septal cells, lateral mispositioning of the septum (>10% deviation from the median position) was observed in ~25% of mutant cells, including 5% of cells with major misplacement (>30% of deviation) (Figure 5D). In multi-septal cells, a similar situation was observed with numerous lateral misplacements of septa, leading in extreme cases to the production of minicells (Supplementary Figure 7). All these observations suggest that LytB is involved in septum maturation and that its absence affects the timing of division and the lateral positioning of the septum.
To further document the role of LytB in cell-cycle dynamics, the ΔlytB mutant was observed by time-lapse microscopy. In Figure 6, two behaviors of ΔlytB mutant cells are shown. In the first case (Figure 6A, Supplementary Movie 4), which is representative of the behavior of many cells, we observed that the cell elongated until it reached 6 μm prior to divide, which is 2 μm longer than the WT at the same pre-divisional state (Figure 4A). At the next generation, the daughter cells did not elongate significantly before starting division (Figure 6A). In an alternative rarer scenario, the left cell elongated until 5 μm and then a first asymmetrical division occurred, which gave rise to a mini cell (Figure 6B, Supplementary Movie 5). Then, a second asymmetrical division took place at the opposite pole in the long daughter cell (Figure 6B). From various time-lapse experiments, ~25% of division events were asymmetrical (n > 100). In addition, we observed that some small cells stopped to grow and finally lysed, which may explain the lower proportion of short cells (< 2 μm) at the whole population level of the ΔlytB mutant (Figure 5B, Supplementary Figure 9).
Figure 6. Cell cycle of the LytB-deficient strain. Time-lapse microscopy of ΔlytB mutant cells showing a deregulation between elongation and division (A) and asymmetrical divisions (B). In (B) the selected cell for the scheme is indicated with a white arrow. Scheme of the cell cycle of two selected cells are displayed on the right. Bacteria were grown on MRS-containing agarose pads at 30°C. The scale bar is 2 μm.
While elongation and division are well-separated processes in the WT (Figure 4A), the ΔlytB mutant is not able to regulate the timing of division during elongation. In addition, LytB deficiency affects the lateral placement of the division site, which could be asymmetrically localized in short or long cells (Figure 6B).
The Combined Inactivation of LytA and LytB Severely Affects Cell Growth
The above results suggest that LytA and LytB play different roles in the cell cycle of L. plantarum. This contrasts with the situation of B. subtilis where elongation is controlled by two D,L-endopeptidases (i.e., CwlO and LytE) (Domínguez-Cuevas et al., 2013; Meisner et al., 2013). Since LytB does not appear to be a rescuer of LytA, we decided to construct a double lytA lytB mutant to determine to what extent their activities interact with each other. Previous attempts to construct a stable mutant deleted for both lytA and lytB (ΔlytB ΔlytA) failed, which was an indirect indication that the combined presence of LytA and lytB is essential for growth (Rolain et al., 2012). To solve this issue, we transferred the PnisA-lytA conditional fusion into the ΔlytB mutant. The resulting double mutant (PnisA-lytA ΔlytB) was obtained in the presence of nisin. Compared to the singly depleted PnisA-lytA mutant, the double PnisA-lytA ΔlytB mutant displayed a severe growth defect (Figure 7A). No growth of the double mutant was observed during the first 6 h while the simple mutant strain started to grow after 2 h post-inoculation (Figure 7A). The retarded growth of the double mutant is likely due to suppressor mutations since the corresponding cells displayed a typical rod-shape morphology (Figure 7A and data not shown). These suppressors are probably resulting from a reversion of the nisin-controlled expression system as previously observed for another conditional mutant of an essential gene (alanine racemase) in L. plantarum (unpublished data). So, this result indicates that the combined presence of LytA and LytB is (nearly) essential for growth and reinforces the key role played by both enzymes in PG remodeling during the cell cycle.
Figure 7. Effect of the double LytA LytB deficiency on growth and cell morphology. (A) Growth curves of WT, PnisA-lytA mutant, ΔlytB mutant and double PnisA-lytA ΔlytB mutant. N0, non-induced; N1 and N25, nisin 1 and 25 ng ml−1, respectively. Bacteria were grown in MRS with chloramphenicol and nisin when appropriate. Curves were generated from triplicates (mean values + standard deviations). (B) Cells of the double PnisA-lytA ΔlytB mutant observed by phase contrast (PC) microscopy and FM4-64 staining under nisin depletion (2 h). Cells were collected from MRS cultures (with chloramphenicol) and observed on agarose pads after suspension in PBS. The scale bar is 2 μm. (C) PnisA-lytA ΔlytB mutant cells expressing FtsZ-GFP+ (Pshp0064-ftsZ-fgp+; + FtsZ-GFP+) observed by epifluorescence microscopy. Bacteria were cultured in MRS with erythromycin and chloramphenicol, and induced with ComS (8 μM, C8) in absence of nisin (N0). Cultures were moderately shaken after induction. The scale bar is 2 μm. (D) Cells of PnisA-lytA ΔlytB mutant complemented with LytA (Pshp0064-lytA; + LytA, left panels) and LytB (Pshp0064-lytB; + LytB, right panels) observed by phase contrast microscopy and FM4-64 staining under nisin depletion (2 h). Bacteria were cultured in MRS with erythromycin and chloramphenicol, and induced with ComS (8 μM, C8) in absence of nisin. Scale bar is 2 μm. For (B–D) similar observations were obtained from at least three independent experiments.
The cell morphology of the double mutant was then examined after 2 h of nisin depletion. Intriguingly, a mixed situation was observed with round aggregated cells and long cells forming chains, which is reminiscent of either LytA or LytB deficiency, respectively (Figure 7B). This mixed phenotype is likely due to cell-to-cell variations in the extent of LytA depletion. We also examined the positioning of septa and FtsZ-GFP+ in the double mutant under nisin depletion (Figures 7B,C). In the long cells, septa and Z ring-like structures formed perpendicularly to the long axis without being localized at mid-cell. This mispositioning is reminiscent of the abnormalities observed with the simple ΔlytB mutant (see Figure 6B). Akin to simple LytA-deficient cells, no clear positioning of FtsZ-GFP+ was observed in round cells. To investigate this heterogeneous cell morphology, we performed time-lapse experiments (Supplementary Figure 10 and Supplementary Movies 6, 7). At the beginning of the depletion, cells were long and formed chains as found with the ΔlytB mutant. Then, they became twisted before their separation in round or unshaped cells resembling to LytA-deficient cells (Supplementary Figure 10). So, the initial observation of a mixed cell morphology under static conditions likely corresponds to different levels of LytA at the early stages of the depletion process.
In order to validate that LytA and LytB have no overlapping function, the double PnisA-lytA ΔlytB mutant was complemented by either lytA or lytB. When the double mutant was grown in the presence of ComS to express the complementing partner and under nisin depletion condition, lytA-complemented cells recovered a rod shape but most cells were longer than WT cells (Figure 7D, left panels); whereas lytB-complemented cells remained round and aggregated (Figure 7D, right panels). In addition, complementation assays of the simple PnisA-lytA mutant (N0) by lytB or the simple ΔlytB mutant by lytA did not restore the morphology of WT cells (data not shown). Finally, complementation of the simple PnisA-lytA mutant by chimeric proteins where the catalytic domains (NlpC/P60 domains, see Figure 1B) were swapped between both proteins (i.e., LytA-NlpC/P60LyB and LytB-NlpC/P60LytA) revealed that the functionality of the LytA protein mainly relies on its accessory domains (Supplementary Figure 11). Indeed, the LytA-NlpC/P60LyB hybrid restored cell elongation while the LytB-NlpC/P60LytA fusion protein did not (Supplementary Figure 11).
Taken together, these results show that the combined PG hydrolytic activity of LytA and LytB is of major importance for L. plantarum growth. They also demonstrate that LytA and LytB play separate and non-redundant functions during the cell cycle.
This work shed a new light on our understanding of the roles played by the four NlpC/P60 endopeptidases identified in L. plantarum. While LytC and LytD do not appear to be involved in morphogenesis, the paralogous LytA and lytB enzymes play well-separated roles in the cell cycle (Figure 8). We reveal that LytA is a key PGH mainly dedicated to cell elongation while LytB is involved in the division process. Notably, the synthetic co-inhibition of these two processes showed that they are both important for L. plantarum growth.
Figure 8. Roles of LytA and LytB in the cell cycle. LytA is mainly involved in cell elongation and in the control of cell diameter. Its absence (probably indirectly) also affects septum positioning. LytA potentially interacts with the elongasome and its basophilic LysM domain suggests that LytA could act in the external PG layers where the pH is higher. LytB is implicated in division (timing, septum maturation, and lateral positioning) and potentially interacts with the divisome. The presence of two acidophilic LysM domains suggests that LytB could be localized close to the cell membrane where the pH is lower. LytC is proposed to be involved in the recycling/catabolism of PG stem peptides while LytD could be an accessory endopeptidase unrelated to the cell cycle. Labeling of accessory domains of Lyt enzymes are as in Figure 1. Dotted arrows indicate potential interactions. PG, peptidoglycan; M, cytoplasmic membrane.
LytC and LytD Are Not Implicated in Morphogenesis During Standard Growth Conditions
Mutant strains deficient for LytC and/or LytD did not show any observable defect in growth and cell morphology under the tested conditions (Supplementary Figure 1). LytC and LytD display a similar structural organization with accessory SH3b PG-binding domains (Figure 1B). Interestingly, LytC is predicted to be a member of the YkfC subfamily of NlpC/P60 D,L-endopeptidases that only hydrolyze the γ-D-Glu-meso-DAP bond of free stem peptides ending with L-Ala (tri, tetra and pentapeptides) (Xu et al., 2015). This activity suggests that members of this subfamily are either involved in the recycling of stem peptide fragments for PG biosynthesis or for their catabolism under nutrient deprivation (Schmidt et al., 2001; Xu et al., 2015). While the physiological role of YkfC remains unknown in B. subtilis, the expression of its encoded gene seemed to be regulated and specifically induced in late vegetative growth phase (Smith et al., 2000). Similarly, the lytC gene was reported to be upregulated by stress conditions (i.e., T° and NaCl) while the other lyt genes were stress-unresponsive (Rolain et al., 2012). In addition, in standard growth conditions, proteomic studies of surface proteins of L. plantarum WCFS1 did not reveal the presence of LytC while LytA, LytB and LytD were shown to be produced (Fredriksen et al., 2013). These observations suggest that LytC may play a physiological role in very specific conditions that remain to be determined and could explain the absence of phenotypes regarding growth and morphogenesis in our study.
Intriguingly, the NlpC/P60 domains of LytD and CwlT are closer compared to any other NlpC/P60 domain of B. subtilis D,L-endopeptidases (Supplementary Figure 2B). CwlT is not involved in the cell cycle as it is involved in the biogenesis of the DNA transfer machinery required for the conjugation of ICEBs1 (Fukushima et al., 2007). The similarity of LytD with CwlT may indicate its involvement in an accessory process such as predation (e.g., fratricins) or the trans-envelope assembly of a secretion/uptake apparatus as previously reported for some members of the NlpC/P60 family (Xu et al., 2015). It is noteworthy that a range of streptococcal fratricins that are specifically induced during competence development are modular enzymes including SH3 domains and closely-related NlpC/P60 domains (Berg et al., 2012). Based on these in silico analyses, additional investigations are thus needed to determine the exact function of LytC and LytD.
LytA Is a Major NlpC/P60 Endopeptidase Involved in the Control of Cell Length
Conditional inactivation of lytA and time-lapse experiments revealed its major role in the cell elongation pathway of L. plantarum (Figure 4). The absence of LytA resulted in elongation arrest together with misplaced septa, leading to unseparated round cells of various sizes (Figures 3, 4, Supplementary Figures 4, 5). Interestingly, complementation studies of the conditional lytA mutant with LytA variants carrying the LytB NlpC/P60 domain, or deleted of either the LysM domain or the AST domain, showed a recovery of elongation without a restoration of cell width, which indicates an additional implication of LytA in the control of the cell diameter (Supplementary Figures 11, 12). Moreover, these dramatic alterations of cell shape were more probably due to a lack of the LytA endopeptidase activity, as shown by complementation studies with a LytA catalytic inactive mutant (Figure 2). In addition, alterations of division planes and cell aggregation do not seem to be the consequence of a lack of elongation as the inactivation of MreBCD, a well-known player of elongation in other rod-shaped Gram-positive bacteria only leads to the formation of well-defined round cells in L. plantarum (Supplementary Figure 6). In B. subtilis, which displays a similar composition of the stem peptide and the same cross-linking (Vollmer et al., 2008a), CwlO and LytE were shown to, respectively control the cell longitudinal axis and cell diameter, being altogether required for the elongation process (Hashimoto et al., 2012; Domínguez-Cuevas et al., 2013; Meisner et al., 2013). Thus, the LytA enzyme of L. plantarum seems to recapitulate the combined activity of CwlO and LytE, which might be consistent with a lower level of PGH redundancy in L. plantarum (Smith et al., 2000; Rolain et al., 2012). However, an intriguing consequence of the absence of LytA activity is the misplacement of septa (Supplementary Figures 4, 5) and a lack of a clear localization of FtsZ in round cells (Figure 3). Apart from the Min system (Bernard et al., 2012), no positive or negative players regulating the positioning of the divisome have been identified in L. plantarum. Our analysis of PG composition suggests that LytA is responsible of the cleavage of most γ-D-Glu-meso-DAP bonds. As products of this hydrolysis represent around 10% of total disaccharides-peptides in the WT profile of muropeptides (Supplementary Table 2), this could indicate that LytA plays an important role in PG remodeling that could indirectly affect either the Min system or unknown players involved in the positioning of the division plane.
LytB Is Involved in Septum Maturation and Timing of Division
Cell measurements, positioning of septa, and time-lapse experiments showed that LytB is a key player in the regulation of cell division in L. plantarum (Figure 8). Indeed, LytB-deficient cells displayed heterogeneous lengths with misplaced division sites along the longitudinal cell axis (Figures 5, 6, Supplementary Figures 7, 8). Our results also suggest that the maturation of the septum is delayed as septa/Z-rings are present in daughter cells while the maturation of the mid-cell septum is still ongoing (Supplementary Figures 7, 8). This delayed septum maturation could contribute to the altered timing between elongation and division such as observed in time-lapse experiments (Figure 6). Moreover, the positioning of septa in daughter cells is laterally affected and seems desynchronized as bi-septal long cells were observed (see Supplementary Figures 7, 8). In B. subtilis, the individual inactivation of LytF, CwlS or LytE, yields filaments of unseparated cells (Fukushima et al., 2006), a phenotype that is reminiscent of that reported for the LytB inactivation in L. plantarum. Similarly, the inactivation of NlpC/P60 endopeptidases of L. casei BL23 (Lc-P75) and L. rhamnosus GG (P75/Msp1) lead to cell separation defects and cell chaining (Claes et al., 2012; Regulski et al., 2012). However, the cell-cycle alterations observed for LytB inactivation seem more complex than a mere cell separation defect. By comparison with LytA inactivation, the absence of LytB did not reveal any alteration in the muropeptide profile compared to the WT (data not shown), which is compatible with a specific role of LytB in PG remodeling (e.g., septum maturation). In addition, we cannot exclude that LytB is intimately associated to the division machinery with a structural role. This could indirectly affect the synchronization and placement of division planes in daughter cells. Alternatively, a local remodeling of PG architecture by LytB could provide a specific signature important for the anchoring of early players and the correct placement of the division machinery. Intriguingly, cell-cycle alterations due to LytB depletion are reminiscent of those previously reported for the MurNAc O-acetyl-transferase mutant of L. plantarum (OatA−) (Bernard et al., 2011a, 2012). As for LytB, the absence of OatA led to an altered timing between division and elongation while OatA overproduction resulted in misplaced septa (Bernard et al., 2012). This suggests that OatA and LytB might contribute to the same process. It is tempting to speculate that PG O-acetylation might locally disturb PG hydrolysis, for example, by affecting LysM-mediated binding of LytB to its substrate. Local PG O-acetylation and PG hydrolysis at mid-cell could in turn act as a cue for the initiation of cytokinesis.
Co-inactivation of LytA and LytB Leads to a Synthetic Growth Defect
Interestingly, co-inactivation of LytA and LytB severely affects growth, which might indicate a synthetic lethality (Figure 7). Time-lapse experiments of LytA depletion in the LytB-deficient strain clearly showed that LytA and LytB are involved in two separate processes (Supplementary Figure 10) since the cell morphology defect results in an evolving mixed phenotype with a dominancy of LytA over LytB regarding cell shape alteration. In addition to a lack of cell elongation due to LytA deficiency, LytA is also involved in the correct positioning of the septum in round cells. Consequently, the combination of a second defect in the division process due to LytB inactivation (i.e., delayed septum maturation and altered selection of the division site) may thus explain the synthetic effect of a double mutation. It is well-established that rod-shaped cells can tolerate a defect in cell elongation but are much more sensitive to alterations in cell division (Szwedziak and Löwe, 2013). This severe growth defect is also reminiscent of the co-lethality observed for the co-inactivation of LytE and CwlO in B. subtilis (Hashimoto et al., 2012). In this case, while both LytE and CwlO are required for cell elongation, LytE is also involved in cell separation (Domínguez-Cuevas et al., 2013). Thus, the synthetic effect of LytA-LytB deficiency on growth in L. plantarum reinforces the essential role played by D,L-endopeptidases of the NlpC/P60 family during the cell cycle of rod-shaped Gram-positive bacteria.
How Are LytA and LytB Controlled?
Bacillus subtilis contains three actin-like cytoskeleton proteins (i.e., MreB, MreBH, and Mbl) that are involved in the control of cell elongation. The activity of the LysM-containing enzyme LytE was proposed to be guided by MreBH (and possibly MreB) while the activity of the coiled-coil domain-containing enzyme CwlO was proposed to be associated to Mbl (Domínguez-Cuevas et al., 2013; Meisner et al., 2013). In addition, CwlO was shown to be under the control of the ABC transporter FtsEX, the latter being required for cell elongation in B. subtilis as opposed to cell division in E. coli (Yang et al., 2011; Domínguez-Cuevas et al., 2013; Meisner et al., 2013). PGHs (e.g., PcsB of Streptococcus pneumoniae) or PGH adaptor proteins (e.g., EnvC in E. coli) that are controlled by the FtsEX complex contain at least one coiled-coil domain that directly interacts with the surface-exposed FtsX protein (Sham et al., 2011; Yang et al., 2011). A similar interaction has been proposed between CwlO and FtsX (Domínguez-Cuevas et al., 2013). Intriguingly, none of the 12 putative PGHs of L. plantarum contains predicted coiled-coil domain (COILS prediction tool, https://toolkit.tuebingen.mpg.de/#/tools/pcoils). Moreover, no canonical FtsEX could be identified in the L. plantarum genome based on a search of conserved domains COG2884 and COG2177 for FtsE and FtsX, respectively. This suggests that LytA and/or LytB might be controlled by another mechanism or different players.
In this work, we also performed complementation experiments of the lytA mutant with chimeric proteins (i.e., LytA-NlpC/P60LyB and LytB-NlpC/P60LytA). The results showed that the accessory domains predominates over the catalytic domain to confer their specific function to the Lyt enzymes (Supplementary Figure 11). In various LysM-containing PGHs, it has been shown that deletion of LysM domain(s) leads to less active or even inactive enzymes (Steen et al., 2005; Layec et al., 2009; Frankel and Schneewind, 2012). In addition, it was observed that the number of LysM domains per PGH (cooperative binding) and their IPs play a major role in their PG-binding capacity (Visweswaran et al., 2011). The presence of one LysM domain with a basic IP in LytA compared to two LysM domains with acidic IPs in LytB may indicate that the LysM domain of LytA have a lower PG-binding capacity in vicinity of the cell membrane, which is more acidic compared to external PG layers (Jolliffe et al., 1981). In this case, LytA would be able to bind PG and be more active in external layers of the cell wall, where PG is more stretched and its hydrolysis required for cell elongation (Figure 8; Lee and Huang, 2013). Conversely, the presence of acidophilic LysM domains in LytB suggests that it could be active in close proximity with the cell membrane, which would be compatible with its role in the maturation and the positioning of the septum.
The roles of the glycosylated AST domains of LytA and LytB remains unexplored. In the Acm2 glucosaminidase of L. plantarum, the AST domain and more specifically its glycosyl residues were shown to negatively impact on the enzymatic activity (Rolain et al., 2013). The purification of glycosylated and non-glycosylated variants of LytA and LytB was attempted in different hosts but were unsuccessful so far (data not shown). Besides a direct impact on enzymatic activity, AST-related PGH domains were hypothesized to be involved in subcellular targeting (Huard et al., 2003; Eckert et al., 2006; Claes et al., 2012; Lebeer et al., 2012; Regulski et al., 2012). For instance, the AST domains of LytA or LytB and their glycolytic decorations might interact with other cell-cycle proteins in order to associate LytA and/or LytB to specific complexes of cell wall biosynthesis.
PH, BH, M-PC-C, and YD conceived and designed the study. M-CD, TR, PC, and AK carried out the laboratory work. M-CD, PC, PH, BH, M-PC-C, and YD analyzed the data. M-CD, PH, and BH wrote the manuscript. All authors read and approved the final manuscript.
Conflict of Interest Statement
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.
Work in the teams of PH, BH, and YD was supported by the National Foundation for Scientific Research (FNRS) and the Research Department of the Communauté française de Belgique (Concerted Research Action). Work in the team of M-PC-C was supported by INRA and Région Ile de France. M-CD and TR held a doctoral fellowship from FRIA. YD and PH are Research Director and Senior Research Associate of the FNRS.
We are grateful to Dr. Johann Mignolet for his critical reading of the manuscript. We warmly thank Hervé Degand and Sylvie Derclaye for their technical assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00713/full#supplementary-material
Supplementary Text. Construction of strains and plasmids.
Supplementary Figure 1. Structural prediction and phylogenic relationships of NlpC/P60 catalytic domains.
Supplementary Figure 2. Impact of LytC and LytD depletion on cell morphology and growth.
Supplementary Figure 3. RP-HPLC separation of muropeptides from L. plantarum WT and PnisA-lytA mutant (without, N0; and with nisin, N25).
Supplementary Figure 4. Effect of LytA depletion on septum misplacement in dividing cells.
Supplementary Figure 5. Effect of LytA depletion on septum misplacement in cell aggregates.
Supplementary Figure 6. Cell cycle of the MreB1CD-deficient strain.
Supplementary Figure 7. Effect of LytB inactivation on septum misplacement in dividing long cells.
Supplementary Figure 8. Effect of LytB inactivation on Z rings positioning.
Supplementary Figure 9. Time-lapse microscopy of ΔlytB mutant cells showing lysis of small non-growing cells.
Supplementary Figure 10. Cell cycle of the double LytA LytB deficient strain.
Supplementary Figure 11. Importance of accessory LytA domains for its morphogenic function.
Supplementary Figure 12. Effect of the individual deletion of LytA domains on cell morphology.
Supplementary Table 1. Primers used in this study.
Supplementary Table 2. Disaccharide (Ds)-peptide composition of PG from L. plantarum WT and PnisA-lytA mutant (without, N0; and with nisin, N25).
Supplementary Movie 1. Cell cycle of the WT (NZ7100).
Supplementary Movie 2. Cell cycle of the conditional lytA mutant (MCD202, N0).
Supplementary Movie 3. Cell cycle of the conditional mreB1CD mutant (MCD203, N0).
Supplementary Movie 4. Cell cycle of the ΔlytB mutant (TR0015), Cell 1.
Supplementary Movie 5. Cell cycle of the ΔlytB mutant (TR0015), Cell 2.
Supplementary Movie 6. Cell cycle of the double LytA-LytB deficient strain (MCD20215, N0), Cell 1.
Supplementary Movie 7. Cell cycle of the double LytA-LytB deficient strain (MCD20215, N0), Cell 2.
Andre, G., Deghorain, M., Bron, P. A., van Swam, I. I., Kleerebezem, M., Hols, P., et al. (2011). Fluorescence and atomic force microscopy imaging of wall teichoic acids in Lactobacillus plantarum. ACS Chem. Biol. 6, 366–376. doi: 10.1021/cb1003509
Barrett, J. F., Dolinger, D. L., Schramm, V. L., and Shockman, G. D. (1984). The mechanism of soluble peptidoglycan hydrolysis by an autolytic muramidase. A processive exodisaccharidase. J. Biol. Chem. 259, 11818–11827.
Berg, K. H., Biørnstad, T. J., Johnsborg, O., and Håvarstein, L. S. (2012). Properties and biological role of streptococcal fratricins. Appl. Environ. Microbiol. 78, 3515–3522. doi: 10.1128/AEM.00098-12
Bernard, E., Rolain, T., Courtin, P., Guillot, A., Langella, P., Hols, P., et al. (2011a). Characterization of O-acetylation of N-acetylglucosamine: a novel structural variation of bacterial peptidoglycan. J. Biol. Chem. 286, 23950–23958. doi: 10.1074/jbc.M111.241414
Bernard, E., Rolain, T., Courtin, P., Hols, P., and Chapot-Chartier, M. P. (2011b). Identification of the amidotransferase AsnB1 as being responsible for meso-diaminopimelic acid amidation in Lactobacillus plantarum peptidoglycan. J. Bacteriol. 193, 6323–6330. doi: 10.1128/JB.05060-11
Bernard, E., Rolain, T., David, B., André, G., Dupres, V., Dufrêne, Y. F., et al. (2012). Dual role for the O-acetyltransferase OatA in peptidoglycan modification and control of cell septation in Lactobacillus plantarum. PLoS ONE 7:e47893. doi: 10.1371/journal.pone.0047893
Carballido-López, R., Formstone, A., Li, Y., Ehrlich, S. D., Noirot, P., and Errington, J. (2006). Actin homolog MreBH governs cell morphogenesis by localization of the cell wall hydrolase LytE. Dev. Cell 11, 399–409. doi: 10.1016/j.devcel.2006.07.017
Claes, I. J., Schoofs, G., Regulski, K., Courtin, P., Chapot-Chartier, M. P., Rolain, T., et al. (2012). Genetic and biochemical characterization of the cell wall hydrolase activity of the major secreted protein of Lactobacillus rhamnosus GG. PLoS ONE 7:e31588. doi: 10.1371/journal.pone.0031588
Courtin, P., Miranda, G., Guillot, A., Wessner, F., Mézange, C., Domakova, E., et al. (2006). Peptidoglycan structure analysis of Lactococcus lactis reveals the presence of an L,D-carboxypeptidase involved in peptidoglycan maturation. J. Bacteriol. 188, 5293–5298. doi: 10.1128/JB.00285-06
Deghorain, M., Goffin, P., Fontaine, L., Mainardi, J. L., Daniel, R., Errington, J., et al. (2007). Selectivity for D-lactate incorporation into the peptidoglycan precursors of Lactobacillus plantarum: role of Aad, a VanX-like D-alanyl-D-alanine dipeptidase. J. Bacteriol. 189, 4332–4337. doi: 10.1128/JB.01829-06
Delcour, J., Ferain, T., Deghorain, M., Palumbo, E., and Hols, P. (1999). The biosynthesis and functionality of the cell-wall of lactic acid bacteria. Antonie Van Leeuwenhoek 76, 159–184. doi: 10.1023/A:1002089722581
Desguin, B., Goffin, P., Bakouche, N., Diman, A., Viaene, E., Dandoy, D., et al. (2015). Enantioselective regulation of lactate racemization by LarR in Lactobacillus plantarum. J. Bacteriol. 197, 219–230. doi: 10.1128/JB.02192-14
Domínguez-Cuevas, P., Porcelli, I., Daniel, R. A., and Errington, J. (2013). Differentiated roles for MreB-actin isologues and autolytic enzymes in Bacillus subtilis morphogenesis. Mol. Microbiol. 89, 1084–1098. doi: 10.1111/mmi.12335
Eckert, C., Lecerf, M., Dubost, L., Arthur, M., and Mesnage, S. (2006). Functional analysis of AtlA, the major N-acetylglucosaminidase of Enterococcus faecalis. J. Bacteriol. 188, 8513–8519. doi: 10.1128/JB.01145-06
Formstone, A., and Errington, J. (2005). A magnesium-dependent mreB null mutant: implications for the role of mreB in Bacillus subtilis. Mol. Microbiol. 55, 1646–1657. doi: 10.1111/j.1365-2958.2005.04506.x
Frankel, M. B., and Schneewind, O. (2012). Determinants of murein hydrolase targeting to cross-wall of Staphylococcus aureus peptidoglycan. J. Biol. Chem. 287, 10460–10471. doi: 10.1074/jbc.M111.336404
Fredriksen, L., Moen, A., Adzhubei, A. A., Mathiesen, G., Eijsink, V. G., and Egge-Jacobsen, W. (2013). Lactobacillus plantarum WCFS1 O-linked protein glycosylation: an extended spectrum of target proteins and modification sites detected by mass spectrometry. Glycobiology 23, 1439–1451. doi: 10.1093/glycob/cwt071
Fukushima, T., Afkham, A., Kurosawa, S., Tanabe, T., Yamamoto, H., and Sekiguchi, J. (2006). A new D,L-endopeptidase gene product, YojL (renamed CwlS), plays a role in cell separation with LytE and LytF in Bacillus subtilis. J. Bacteriol. 188, 5541–5550. doi: 10.1128/JB.00188-06
Fukushima, T., Uchida, N., Ide, M., Kodama, T., and Sekiguchi, J. (2018). DL-endopeptidases function as both cell wall hydrolases and poly-gamma-glutamic acid hydrolases. Microbiology 164, 277–286. doi: 10.1099/mic.0.000609
Fukushima, T., Yao, Y., Kitajima, T., Yamamoto, H., and Sekiguchi, J. (2007). Characterization of new L,D-endopeptidase gene product CwlK (previous YcdD) that hydrolyzes peptidoglycan in Bacillus subtilis. Mol. Genet. Genomics 278, 371–383. doi: 10.1007/s00438-007-0255-8
Hashimoto, M., Ooiwa, S., and Sekiguchi, J. (2012). Synthetic lethality of the lytE cwlO genotype in Bacillus subtilis is caused by lack of D,L-endopeptidase activity at the lateral cell wall. J. Bacteriol. 194, 796–803. doi: 10.1128/JB.05569-11
Huard, C., Miranda, G., Wessner, F., Bolotin, A., Hansen, J., Foster, S. J., et al. (2003). Characterization of AcmB, an N-acetylglucosaminidase autolysin from Lactococcus lactis. Microbiology 149, 695–705. doi: 10.1099/mic.0.25875-0
Josson, K., Scheirlinck, T., Michiels, F., Platteeuw, C., Stanssens, P., Joos, H., et al. (1989). Characterization of a gram-positive broad-host-range plasmid isolated from Lactobacillus hilgardii. Plasmid 21, 9–20. doi: 10.1016/0147-619X(89)90082-6
Kleerebezem, M., Boekhorst, J., van Kranenburg, R., Molenaar, D., Kuipers, O. P., Leer, R., et al. (2003). Complete genome sequence of Lactobacillus plantarum WCFS1. Proc. Natl. Acad. Sci. U. S. A 100, 1990–1995. doi: 10.1073/pnas.0337704100
Kleerebezem, M., Hols, P., Bernard, E., Rolain, T., Zhou, M., Siezen, R. J., et al. (2010). The extracellular biology of the lactobacilli. FEMS Microbiol. Rev. 34, 199–230. doi: 10.1111/j.1574-6976.2009.00208.x
Kuipers, O. P., de Ruyter, P. G., Kleerebezem, M., and de Vos, W. M. (1997). Controlled overproduction of proteins by lactic acid bacteria. Trends Biotechnol. 15, 135–140. doi: 10.1016/S0167-7799(97)01029-9
Layec, S., Gérard, J., Legué, V., Chapot-Chartier, M. P., Courtin, P., Borges, F., et al. (2009). The CHAP domain of Cse functions as an endopeptidase that acts at mature septa to promote Streptococcus thermophilus cell separation. Mol. Microbiol. 71, 1205–1217. doi: 10.1111/j.1365-2958.2009.06595.x
Lebeer, S., Claes, I. J., Balog, C. I., Schoofs, G., Verhoeven, T. L., Nys, K., et al. (2012). The major secreted protein Msp1/p75 is O-glycosylated in Lactobacillus rhamnosus GG. Microb. Cell Fact. 11:15. doi: 10.1186/1475-2859-11-15
Litzinger, S., Fischer, S., Polzer, P., Diederichs, K., Welte, W., and Mayer, C. (2010). Structural and kinetic analysis of Bacillus subtilis N-acetylglucosaminidase reveals a unique Asp-His dyad mechanism. J. Biol. Chem. 285, 35675–35684. doi: 10.1074/jbc.M110.131037
Margot, P., Pagni, M., and Karamata, D. (1999). Bacillus subtilis 168 gene lytF encodes a gamma-D-glutamate-meso-diaminopimelate muropeptidase expressed by the alternative vegetative sigma factor, sigmaD. Microbiology 145(Pt 1), 57–65. doi: 10.1099/13500872-145-1-57
Meisner, J., Montero Llopis, P., Sham, L. T., Garner, E., Bernhardt, T. G., and Rudner, D. Z. (2013). FtsEX is required for CwlO peptidoglycan hydrolase activity during cell wall elongation in Bacillus subtilis. Mol. Microbiol. 89, 1069–1083. doi: 10.1111/mmi.12330
Mesnage, S., Dellarole, M., Baxter, N. J., Rouget, J. B., Dimitrov, J. D., Wang, N., et al. (2014). Molecular basis for bacterial peptidoglycan recognition by LysM domains. Nat. Commun. 5:4269. doi: 10.1038/ncomms5269
Palumbo, E., Deghorain, M., Cocconcelli, P. S., Kleerebezem, M., Geyer, A., Hartung, T., et al. (2006). D-alanyl ester depletion of teichoic acids in Lactobacillus plantarum results in a major modification of lipoteichoic acid composition and cell wall perforations at the septum mediated by the Acm2 autolysin. J. Bacteriol. 188, 3709–3715. doi: 10.1128/JB.188.10.3709-3715.2006
Palumbo, E., Favier, C. F., Deghorain, M., Cocconcelli, P. S., Grangette, C., Mercenier, A., et al. (2004). Knockout of the alanine racemase gene in Lactobacillus plantarum results in septation defects and cell wall perforation. FEMS Microbiol. Lett. 233, 131–138. doi: 10.1016/j.femsle.2004.02.001
Regulski, K., Courtin, P., Meyrand, M., Claes, I. J., Lebeer, S., Vanderleyden, J., et al. (2012). Analysis of the peptidoglycan hydrolase complement of Lactobacillus casei and characterization of the major gamma-D-glutamyl-L-lysyl-endopeptidase. PLoS ONE 7:e32301. doi: 10.1371/journal.pone.0032301
Rolain, T., Bernard, E., Beaussart, A., Degand, H., Courtin, P., Egge-Jacobsen, W., et al. (2013). O-glycosylation as a novel control mechanism of peptidoglycan hydrolase activity. J. Biol. Chem. 288, 22233–22247. doi: 10.1074/jbc.M113.470716
Rolain, T., Bernard, E., Courtin, P., Bron, P. A., Kleerebezem, M., Chapot-Chartier, M. P., et al. (2012). Identification of key peptidoglycan hydrolases for morphogenesis, autolysis, and peptidoglycan composition of Lactobacillus plantarum WCFS1. Microb. Cell Fact. 11:137. doi: 10.1186/1475-2859-11-137
Sauvage, E., Kerff, F., Terrak, M., Ayala, J. A., and Charlier, P. (2008). The penicillin-binding proteins: structure and role in peptidoglycan biosynthesis. FEMS Microbiol. Rev. 32, 234–258. doi: 10.1111/j.1574-6976.2008.00105.x
Schmidt, D. M., Hubbard, B. K., and Gerlt, J. A. (2001). Evolution of enzymatic activities in the enolase superfamily: functional assignment of unknown proteins in Bacillus subtilis and Escherichia coli as L-Ala-D/L-Glu epimerases. Biochemistry 40, 15707–15715. doi: 10.1021/bi011640x
Serrano, L. M., Molenaar, D., Wels, M., Teusink, B., Bron, P. A., de Vos, W. M., et al. (2007). Thioredoxin reductase is a key factor in the oxidative stress response of Lactobacillus plantarum WCFS1. Microb. Cell Fact. 6:29. doi: 10.1186/1475-2859-6-29
Sham, L. T., Barendt, S. M., Kopecky, K. E., and Winkler, M. E. (2011). Essential PcsB putative peptidoglycan hydrolase interacts with the essential FtsXSpn cell division protein in Streptococcus pneumoniae D39. Proc. Natl. Acad. Sci. U. S. A 108, E1061–E1069. doi: 10.1073/pnas.1108323108
Sliusarenko, O., Heinritz, J., Emonet, T., and Jacobs-Wagner, C. (2011). High-throughput, subpixel precision analysis of bacterial morphogenesis and intracellular spatio-temporal dynamics. Mol. Microbiol. 80, 612–627. doi: 10.1111/j.1365-2958.2011.07579.x
Sørvig, E., Mathiesen, G., Naterstad, K., Eijsink, V. G., and Axelsson, L. (2005). High-level, inducible gene expression in Lactobacillus sakei and Lactobacillus plantarum using versatile expression vectors. Microbiology 151, 2439–2449. doi: 10.1099/mic.0.28084-0
Steen, A., Buist, G., Horsburgh, G. J., Venema, G., Kuipers, O. P., Foster, S. J., et al. (2005). AcmA of Lactococcus lactis is an N-acetylglucosaminidase with an optimal number of LysM domains for proper functioning. FEBS J. 272, 2854–2868. doi: 10.1111/j.1742-4658.2005.04706.x
Visweswaran, G. R., Dijkstra, B. W., and Kok, J. (2011). Murein and pseudomurein cell wall binding domains of bacteria and archaea–a comparative view. Appl. Microbiol. Biotechnol. 92, 921–928. doi: 10.1007/s00253-011-3637-0
Xu, Q., Abdubek, P., Astakhova, T., Axelrod, H. L., Bakolitsa, C., Cai, X., et al. (2010). Structure of the gamma-D-glutamyl-L-diamino acid endopeptidase YkfC from Bacillus cereus in complex with L-Ala-gamma-D-Glu: insights into substrate recognition by NlpC/P60 cysteine peptidases. Acta Crystallogr. Sect. F. Struct. Biol. Cryst. Commun. 66, 1354–1364. doi: 10.1107/S1744309110021214
Xu, Q., Mengin-Lecreulx, D., Liu, X. W., Patin, D., Farr, C. L., Grant, J. C., et al. (2015). Insights into substrate specificity of NlpC/P60 cell wall hydrolases containing bacterial SH3 domains. MBio. 6:e02327-14. doi: 10.1128/mBio.02327-14
Yamaguchi, H., Furuhata, K., Fukushima, T., Yamamoto, H., and Sekiguchi, J. (2004). Characterization of a new Bacillus subtilis peptidoglycan hydrolase gene, yvcE (named cwlO), and the enzymatic properties of its encoded protein. J. Biosci. Bioeng. 98, 174–181. doi: 10.1016/S1389-1723(04)00262-2
Yamamoto, H., Kurosawa, S., and Sekiguchi, J. (2003). Localization of the vegetative cell wall hydrolases LytC, LytE, and LytF on the Bacillus subtilis cell surface and stability of these enzymes to cell wall-bound or extracellular proteases. J. Bacteriol. 185, 6666–6677. doi: 10.1128/JB.185.22.6666-6677.2003
Yang, D. C., Peters, N. T., Parzych, K. R., Uehara, T., Markovski, M., and Bernhardt, T. G. (2011). An ATP-binding cassette transporter-like complex governs cell-wall hydrolysis at the bacterial cytokinetic ring. Proc. Natl. Acad. Sci. U.S.A 108, E1052–E1060. doi: 10.1073/pnas.1107780108
Keywords: Lactobacillus, cell cycle, cell wall, muropeptidase, NlpC/P60 endopeptidase, peptidoglycan hydrolase, morphogenesis
Citation: Duchêne M-C, Rolain T, Knoops A, Courtin P, Chapot-Chartier M-P, Dufrêne YF, Hallet BF and Hols P (2019) Distinct and Specific Role of NlpC/P60 Endopeptidases LytA and LytB in Cell Elongation and Division of Lactobacillus plantarum. Front. Microbiol. 10:713. doi: 10.3389/fmicb.2019.00713
Received: 17 September 2018; Accepted: 21 March 2019;
Published: 12 April 2019.
Edited by:Tobias Dörr, Cornell University, United States
Reviewed by:Christoph Mayer, University of Tübingen, Germany
Cara C. Boutte, University of Texas at Arlington, United States
Felipe Cava, Umeå University, Sweden
Copyright © 2019 Duchêne, Rolain, Knoops, Courtin, Chapot-Chartier, Dufrêne, Hallet and Hols. 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) and the copyright owner(s) 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: Pascal Hols, firstname.lastname@example.org