Abstract
SpoIIE is a phosphatase involved in the activation of the first sigma factor of the forespore, σF, during sporulation. A ΔspoIIE mutant of Clostridium beijerinckii NCIMB 8052, previously generated by CRISPR-Cas9, did not sporulate but still produced granulose and solvents. Microscopy analysis also showed that the cells of the ΔspoIIE mutant are elongated with the presence of multiple septa. This observation suggests that in C. beijerinckii, SpoIIE is necessary for the completion of the sporulation process, as seen in Bacillus and Clostridium acetobutylicum. Moreover, when grown in reactors, the spoIIE mutant produced higher levels of solvents than the wild type strain. The impact of the spoIIE inactivation on gene transcription was assessed by comparative transcriptome analysis at three time points (4 h, 11 h and 23 h). Approximately 5% of the genes were differentially expressed in the mutant compared to the wild type strain at all time points. Out of those only 12% were known sporulation genes. As expected, the genes belonging to the regulon of the sporulation specific transcription factors (σF, σE, σG, σK) were strongly down-regulated in the mutant. Inactivation of spoIIE also caused differential expression of genes involved in various cell processes at each time point. Moreover, at 23 h, genes involved in butanol formation and tolerance, as well as in cell motility, were up-regulated in the mutant. In contrast, several genes involved in cell wall composition, oxidative stress and amino acid transport were down-regulated. These results indicate an intricate interdependence of sporulation and stationary phase cellular events in C. beijerinckii.
Introduction
Even though butanol is nowadays mainly produced through the petrochemical route, it used to be made industrially by a bioprocess called ABE fermentation in the first half of the 20th century. This process returned to the forefront at the end of the 1990s with the emerging interest for biobased chemicals. ABE fermentation relies on the ability of several bacteria from the Clostridium genus to convert carbohydrates to acetone, ethanol, butanol (ABE) and isopropanol. Clostridia are anaerobic bacteria that can form spores to protect themselves from unfavorable environmental conditions, including oxygen exposure. The main representatives of the solventogenic clostridia group are Clostridium acetobutylicum, C. beijerinckii, C. saccharobutylicum and C. saccharoperbutylacetonicum. These clostridia produce solvents while they form spores. Once the spores are mature, the solvent producing cells lyse, and the metabolically inactive spores are left behind (). That is why in industry, spores are seen as undesirable (Tracy et al., 2012; ), and many efforts were made to engineer asporogenous solvent producing strains (; ; ). However, the sporulation process and the associated regulatory network in these microorganisms are still poorly characterized ().
In solventogenic clostridia, several studies revealed a link between sporulation and solvent production (, ) but, none was able to explain the involved mechanism. The regulatory pathway controlling sporulation was first described and intensively studied in Bacillus subtilis, which is considered as a model organism for the sporulation process. Comparative studies between bacilli and clostridia show similarities in the sporulation process and its regulation, including the presence of the main actors such as Spo0A and σH as well as the four sporulation specific sigma factors, σF, σE, σG and σK (). However, differences in the sporulation regulatory networks are also observed between bacilli and clostridia and even among clostridia (). Important deviations from the B. subtilis paradigm exist in clostridial spore formers, especially concerning the communication between the forespore and the mother cell, a weaker connection between gene expression and morphogenesis, and modifications in the interplay between sigma factors (; ).
In B. subtilis, the SpoIIE protein is a phosphatase that plays a crucial role in the sporulation regulation mechanism (). SpoIIE acts in stage II of the sporulation process. SpoIIE plays a central role in the asymmetric septum formation separating the mother cell and the forespore. Studies in B. subtilis showed that SpoIIE interacts with cell division proteins and peptidoglycan synthesis proteins to enable a correct localization and thickness of the asymmetric septum (; ; ). Following the asymmetric division, SpoIIE enables the activation of σF, the first sigma factor of the forespore (Figure 1; ; ). σF is held inactive by the anti-sigma factor and kinase SpoIIAB. At the beginning of stage II of the sporulation cascade, SpoIIE dephosphorylates the anti-anti-sigma factor SpoIIAA to enable its interaction with SpoIIAB, which releases σF following asymmetric division (; Figure 1). In C. acetobutylicum, considered as the model solventogenic Clostridium, SpoIIE has also been reported to function as a phosphatase involved in the early stages of the sporulation regulation cascade. Previous studies have also shown that spoIIE mutants of C. acetobutylicum were asporogenous but still produced solvents (; ). However, no studies to date confirmed if this model can be applied to other ABE-producing strains.
FIGURE 1
Clostridium beijerinckii is the second most studied solventogenic specie and was used industrially for acetone production already at the beginning of the 20th century (
Materials and Methods
Bacterial Strains and Culture Conditions
Bacterial strains and plasmids are listed in Table 1. The C. beijerinckii wild-type (WT) strain was stored as spore suspension and the mutants as vegetative cells in 15% glycerol solution at – 20°C. Spore suspensions were heat-shocked 1 min at 98°C before inoculation in a liquid medium to kill any vegetative cell present and enable the germination of the spores. Except for fermentation assays, liquid cultures of the WT, mutant and complemented strains were grown in liquid modified CGM (mCGM) as described previously (
TABLE 1
| Strains or plasmids | Relevant characteristics | Reference |
| Strains | ||
| C. beijerinckii NCIMB 8052 (WT) | Wild type, sensitive to spectinomycin (650 μg/mL) and erythromycin (25–50 μg/mL) | NCIMB |
| C.beijerinckii ΔspoIIE | NCIMB 8052, ΔspoIIE | ( |
| E. coli XL1-blue | recA1 endA1 gyrA96 thi-1 hsdR17 supE44 relA1 lac [F′ proAB lacIq ZΔM15 Tn10 (Tetr)] | Agilent |
| Plasmids | ||
| pSpoIIE | pCB102, colE1, aad9, Cbei_0097 | ( |
| pRAN473 | repA, colE1, Ptet:mCherryOpt–MCS, catP | ( |
| pRAN473S | repA, colE1, Ptet:mCherryOpt–MCS, aad9 | This study |
| pRAN473S: cbei_0097 | repA, colE1, Ptet:mCherryOpt–cbei_0097, aad9 | This study |
Strains and plasmids used in this study; catP, chloramphenicol resistance gene; aad9, spectinomycin resistance gene.
Plasmid Construction
The primers used for plasmid construction are listed in the Supplementary Table 1 and were synthesized by Integrated DNA Technologies. To replace the catp gene by aad9 in pRAN473, pRAN473 was linearized by PCR using the primers M392 and M393. The aad9 gene was amplified from pSpoIIE with the primers M394 and M395. The fragments were fused using the Circular Polymerase Extension Cloning (
Fermentation
Fermentations were performed at 35°C in CM2 medium (van der Wal et al., 2013), which contains per liter: yeast extract, 5 g; KH2PO4, 1 g; K2HPO4, 0.76 g; ammonium acetate, 3 g; p-aminobenzoic acid, 0.10 g; MgSO4⋅7 H2O, 1 g; and FeSO4⋅7 H2O, 0.50 g, glucose, 60 g. Metabolites were determined in culture supernatants after removal of cells by centrifugation. Glucose, acetate, butyrate, lactate, acetone, butanol and ethanol concentration in the culture medium were determined by high-performance liquid chromatography (HPLC) as described previously using 4 methyl valeric acid (30 mM) as an internal standard (
DNA Extraction and Sequencing
Genomic DNA of C. beijerinckii NCIMB 8052 and C. beijerinckii ΔspoIIE mutants was purified using the GenElute bacterial genomic DNA kit (Sigma-Aldrich, United States). The concentration of genomic DNA was determined using a nanodrop spectrophotometer (Thermo Fisher Scientific, United States) and quality checked on 1% agarose gel. PCR reactions were carried out using the Q5 Master mix (New England Biolabs, United States). DNA sequencing of clones and genome assembly were performed by BaseClear (Leiden, Netherlands). The sequences of the WT and the ΔspoIIE clones were compared to the publicly available sequence of C. beijerinckii NCIMB 8052 on NCBI. SNPs between the genome of our WT and the ΔspoIIE mutant’s genome with a frequency above 98% were considered in our study.
Granulose Staining
Granulose accumulation was monitored by iodine staining. Each C. beijerinckii strain was grown on CM2 agar plates and incubated anaerobically at 37°C. After 24 h of incubation, the petri dish was opened and inverted over I2 crystals for approximately 1 min. The colonies of granulose-negative mutants were unstained by the sublimed I2 vapors, while the granulose-positive strains were labeled (
Spore Viability Assay
To verify the presence of viable spores, overnight cultures of each strain were prepared in 5 mL CM2 medium. The next day, two tubes containing 15 mL of fresh CM2 liquid medium were inoculated with 50 μL of overnight culture. Aliquots of 100 μL were collected at 24 h and 48 h, treated at 98oC for 1 min and used for 4 serial dilutions. 50 μL of each dilution was spread on CM2 agar plates in the anaerobic chamber and incubated for 24 h at 37°C. The colonies on each plate were then counted to determine the average number of viable spores per mL (spores/mL). This method was adapted from (Steiner et al., 2011).
Microscopy Analysis
Phase-contrast microscopy (Olympus BX51) was used to observe the morphology of WT and ΔspoIIE strains at ×400 and x1000 magnifications. Cells were cultivated for 72 h in liquid CM2 medium, samples were collected at 48 and 72 h, and centrifuged to stain them according to the Schaeffer-Fulton technique. A cell film was made on a glass slide and stained with malachite green and safranin to visualize spores (
Transcriptome Analysis
RNA Isolation and Sequencing Protocol
The same procedure was repeated three times in three different weeks to obtain three independent biological replicates. Each week, a fresh preculture was used to inoculate two identical bioreactors. The cultures were grown, as described in section “Fermentation.” Samples were taken over the early exponential, late exponential and stationary phases (samples at 4, 11 and 23 h). Following centrifugation of the samples, cell pellets were washed with chilled RNase free water and resuspended in RNase free water to obtain a suspension having an OD600nm of approximately 1. A 3-mL diluted sample was centrifuged, the supernatant was discarded, and the cell pellet was stored at −80°C for subsequent isolation. Frozen samples were thawed on ice, and RNA was isolated using High pure RNA isolation kit (Roche Diagnostics, Switzerland). Quality and concentration of RNA samples were checked using a nanodrop spectrophotometer (Thermo Fisher Scientific, United States). The absence of DNA in the RNA samples was evaluated by qPCR analysis performed with BioRad CFX 96 TouchTM (BioRad, Hercules, United States) and the PowerUP SYBr green reaction mix (Applied Biosystems, Thermo Fisher Scientific, United States). Reactions were performed in an overall volume of 10 μL with concentrations of components and reaction conditions, as described in the master mix protocol. RNA quality and integrity were determined using the Qsep 100 bioanalyzer (Bioptic Inc., Taiwan).
All the RNA samples collected were used for library construction and sequencing. The Ribo-zero kit (Illumina, United States) was used to enrich the samples in mRNA. The stranded library was prepared using the TruSeq Stranded mRNA Library Prep Kit (Illumina, United States) according to the manufacturer’s recommendation. AMPure XP beads (Beckman Coulter, United States) were used to clean up the cDNA fragments after each process. Library quality was checked using the bioanalyzer, and the library was then loaded onto HiSeq 2500 for high-throughput sequencing.
Bioinformatics Analysis
After Illumina sequencing, all the reads were mapped to the C. beijerinckii genome using Bowtie (
RT-qPCR analyses were performed with BioRad CFX 96 TouchTM and the PowerUP SYBr green reaction mix to confirm the RNA seq results. Primer3 website was used for oligonucleotide design (Supplementary Table S2). Relative expression at 23 h of 12 genes was monitored (Supplementary Table S2); gene gyrA was chosen as the reference gene from a selection of candidate genes (including gyrA, 16S rna, polIII and alaS) based on analysis by RefFinder algorithms to verify the stability of their expression. Reaction efficiency was determined for each assay using a ×5 serial dilution of cDNA samples. A sample collected after 4 h of cultivation in CM2 medium was chosen as a calibrator and a reference value for calculation of expression fold-change of each gene in other samples. All RT-qPCR analyses were performed in triplicate. The relative quantification was evaluated using the mathematical model described by Dr. Pfaffl (
Results
SpoIIE Disruption and Mutant Characterization During the Sporulation Cycle
The spoIIE gene (Cbei_0097) from C. beijerinckii NCIMB 8052 is 2412 bp long and encodes an 803 amino acid long transmembrane protein. The spoIIE gene is present in all the strains of C. beijerinckii with a genome available. In all C. beijerinckii strains, spoIIE shows a very high level of identity ranging from 96% to 100%. SpoIIE orthologs were found in all spore-forming Bacilli and Clostridia (
FIGURE 2

SpoIIE a sporulation protein conserved in spore-forming Firmicutes (A) Illustration of the predicted SpoIIE protein with its transmembrane regions (according to the TMHMM software) and its phosphatase domain terminus, a.a: amino acids (B) Synteny map of the spoIIE region from Cbei_0084 to Cbei_0107, with spoIIE (Cbei_0097), in C. beijerinckii NCIMB 8052 (CBEI) and its homologs in B. subtilis str.168 (BSU), C. acetobutylicum ATCC 824 (CAC), C. perfringens str.13 (CPE) and C. ljundahlii ATCC 49587 (CLUN), the homologous regions (HSP) are highlighted by color blocks for each organism. Image generated using the COGE platform https://genomevolution.org/r/1ceci, and the GEvo tool (Tang et al., 2015).
Using a xylose inducible CRISPR-Cas9 system, we recently constructed a spoIIE mutant by deleting a 2.379 kb fragment of the coding sequence of spoIIE (
FIGURE 3

Phenotypic comparison of ΔspoIIE mutant and WT. The ΔspoIIE mutant could not produce viable spores but still produced granulose and solvents. (A) Spore viability assay; (B) Granulose detection after 24 h of incubation; (C) Fermentation end products after 72 h of culture at 35°C in 50-mL serum bottles, the error bars indicate one standard deviation of the mean, which was determined based on the data from biological duplicates (n = 2).
We observed the morphology of the cells after 6 h, 24 h, 35 h, 48 h and 72 h of cultivation with phase-contrast microscopy (Figure 4) and after 20 h with fluorescent microscopy (Figure 5A, Supplementary Figures S2A,B). While mature spores were detected after 24 h of growth in WT cultures, no regular pre-spores nor spores were seen in ΔspoIIE mutant cultures, even after 72 h of incubation (Figure 4). At 48 h and 72 h, several mutant cells harbored phase-dark masses at the poles of the cells (Figure 4B black arrows) that were absent in WT cells. Even in the complemented ΔspoIIE mutant, in which sporulation was restored, some cells still harbored these phase-dark bodies at the poles (
FIGURE 4

Phase-contrast microscopy images of C. beijerinckii WT and ΔspoIIE cells during the fermentation. While mature spores were seen in the WT images after 24 h of cultivation, no spores were seen in ΔspoIIE cultures even after 72 h of cultivation. However phase-dark bodies, stained in blue by Shaeffer-Fulton staining, were observed in the ΔspoIIE cells after 48 h of cultivation (A) Pictures with x400 magnification at 6 h, 24 h and 35 h; (B) Pictures at 48 h and 72 h of culture with x1000 magnification with and without Shaeffer-Fulton stain. The short dark blue arrows indicate mature spores. The black arrows indicate the phase-dark masses at the poles of the mutant cells.
FIGURE 5

Fluorescence microscopy images of C. beijerinckii WT and ΔspoIIE cells (A) Fluorescence images of the wild type and the mutant cells after 20 h of cultivation and stained by the membrane staining MTG; (B) Fluorescence images of the ΔspoIIE mutant cells stained by MTG and harboring either the mCherry empty plasmid (pRAN73S) or the plasmid expressing mCherry fused to Cbei_0097 (pRAN73S:cbei_0097) at 18 h of cultivation (after 8 h of atc induction). The white arrows indicate septa observed in the mutant strain. The short white blue arrows indicate wild type cells with an asymmetric septum. The light blue arrows indicate the prespores observed in the wild type.
To localize the SpoIIE protein in the cells, SpoIIE was fused to the mCherry fluorescent protein (
Comparison of the Fermentation Profile of the Wild Type Strain and ΔspoIIE Mutant
The WT and ΔspoIIE strains were grown in a chemostat. Growth, substrate consumption, product formation and pH were monitored for 73 h (Figure 6 and Table 2). The disruption of spoIIE did not affect biomass creation as the mutant strain reached the stationary phase after 23 h, like the WT strain (Figure 6A). As observed in small scale fermentations (Figure 3C), the mutant consumed glucose to produce acids and solvents. However, the pH in the ΔspoIIE mutant culture dropped slightly earlier than in the WT culture (Figure 6A). The reassimilation of the acids started earlier in the mutant culture as the solvent titer at 11 h was twice higher than in the mutant culture compared to the WT. After the reassimilation of butyrate and acetate, the pH at 11 h rose only to 5.7 ± 0.1 in the mutant culture while reaching 6.1 ± 0.0 in the WT culture. From 11 h to 23 h, the pH decreased in both cultures abruptly to 5.0 ± 0.1 in the ΔspoIIE mutant culture and 5.3 ± 0.0 in the WT culture. After 23 h, the pH increased slightly to 5.6 ± 0.0 in the WT culture but still decreased in the mutant culture to 4.9 ± 0.2. At the end of the fermentation, higher amounts of butyrate and acetate accumulated in the mutant culture. Indeed at 49 h, the acetate concentration reached 0.8 g.L–1 ± 0.2 in the mutant culture and 0.3 g.L–1 ± 0.0 in the WT culture. The concentration of acids after 73 h of cultivation in the chemostat was higher in the ΔspoIIE mutant culture, as observed in the serum bottles (Figure 3C). This difference is probably due to the absence of reassimilation of these acids in the ΔspoIIE mutant. The mutant also produced more acetone and butyric acid than the WT (Table 2); but, we did not observe an elongation of the solvent production phase but rather a switch to acidogenesis after 35 h.
FIGURE 6

Fermentation profiles of C. beijerinckii strains. Fermentations were performed in biological and technical duplicates in CM2 medium in chemostats at 35°C during 73 h of cultivation. (A) Growth curve and pH variation in the cultures; (B) Acids and solvent titers during the fermentation. The error bars indicate one standard deviation of the mean, which was determined based on the data from biological and technical replicates (n = 4).
TABLE 2
| Wild type | ΔspoIIE | |
| Substrates consumed [g.L–1] | ||
| Glucose | 32.5 ± 0.5 | 37.0 ± 1.4 |
| Acetate | 1.1 ± 0.1 | 0.8 ± 0.1 |
| Products at the end [g.L–1] | ||
| Acetate* | 0.5 ± 0.0 | 0.7 ± 0.1 |
| Butyrate | 0.6 ± 0.0 | 0.7 ± 0.1 |
| Lactate | 0.5 ± 0.1 | 0.1 ± 0.0 |
| Acetone | 1.9 ± 0.1 | 2.7 ± 0.2 |
| Butanol | 5.8 ± 0.1 | 6.3 ± 0.3 |
| Ethanol | n.d | n.d |
Fermentation data on substrate consumption and product formation of cultures of WT and ΔspoIIE strains after 30 h of cultivation in bioreactors, one standard deviation of the mean was determined based on the data from independent biological and technical duplicates (n = 4) * acetate was also produced by the mutant strain, n.d = not detected.
Impact of the spoIIE Inactivation on the Transcriptome
Overview of the Transcription Data
To study the repercussion of spoIIE’s disruption on the transcriptome, samples for mRNA isolation were collected from three independent chemostat fermentations of the WT and the ΔspoIIE mutant at three time points (4, 11, and 23 h), corresponding to early exponential, mid-exponential and entry into the stationary phase. After RNA isolation, library construction and sequencing, the data were mapped against the published genome (NCBI). The differential expression was calculated using the SARtools pipeline on the Institut Pasteur Galaxy platform (Varet et al., 2016). Out of the 5026 coding genes annotated by NCBI Prokaryotic Genome Annotation Pipeline (PGAP), 5021 were detected in our transcriptomic data. We then compared the expression profile of the WT and ΔspoIIE mutant strains. The inactivation of spoIIE had a significant impact on the transcriptome, as about 40% of the total CDS was differentially expressed in the mutant at least at one time point (2005 out of the 5021 genes transcribed). The differentially expressed genes were clustered using the Clusters of Orthologous Groups (COGs) database (NCBI). Thirty-seven percent of these 2005 genes encode for proteins of unknown function. The rest encodes mainly proteins involved in sporulation, metabolism, signal transduction, and the membrane/cell wall biogenesis. At each time point, more than 70% of the genes were significantly down-regulated (log2 fold change < −1.5 and padj < 0.05) (Figure 7A). The difference between mutant and WT transcriptome was the largest at 23 h. Indeed, while the number of differentially expressed genes with a | log2 fold | > 1.5 and padj < 0.05 was equal to 589 and 578 at 4 h and 11 h, respectively, it reached 1312 genes at 23 h. This increase at 23 h was mainly due to a rise in the number of differentially expressed genes involved in metabolism (from 147 genes at 11 h to 432 at 23 h) and encoding proteins of unknown function (from 218 to 469).
FIGURE 7

Overall differential expression dynamics. (A) Volcano plots of each comparison. Red dots represent significantly differentially expressed features; (B) Venn diagram showing the number of differentially expressed genes in the mutant at each time point; (C) COG class repartition of the 243 genes differentially expressed at all time points in the ΔspoIIE strain: A, Sporulation; B, Cell wall/membrane/envelope biogenesis; C, Amino acid transport and metabolism; D, Acidogenesis and Solventogenesis; E, Stress response; F, Other posttranslational modification, protein turnover, chaperones; G, Secondary metabolites biosynthesis, transport and catabolism; H, Coenzyme transport and metabolism; I, Carbohydrate transport and metabolism; J, Signal transduction mechanisms; K, Replication, recombination and repair; L, Inorganic ion transport and metabolism; M, Intracellular trafficking, secretion, and vesicular transport; N, Lipid transport and metabolism; O, Energy production and conversion; P, Transcription; Q, Defense mechanisms. (D) Proportion up and down-regulated genes at each time point in the cluster of the genes differentially expressed at all time points in the ΔspoIIE strain. The RNAseq data analyzed contains three independent biological replicates per time point and per strain (n = 3).
Two hundred and forty-three genes were significantly up or down-regulated (|log2 fold | > 1.5 and padj < 0.05) in the mutant compared to the WT at the three time points (Figure 7B and Supplementary Table S4), 97% of them being down-regulated at each time point in the ΔspoIIE mutant. And 45% of these 243 genes encode proteins of unknown function while the remaining genes encode proteins involved in various cellular processes such as sporulation, solventogenesis, or chemotaxis (Figure 7C). Only four genes were up-regulated at the three time points. They form an operon of 4 genes annotated as a histidine kinase, two diguanylate cyclases, and a transcription regulator.
To analyze in more detail the data, the genes were clustered according to the putative function of the encoded proteins. In addition to the differential expression of the genes involved in sporulation, the differential expression of the genes linked to five additional cellular processes was analyzed: stress response, cell wall and membrane formation, signal transduction and motility, carbohydrate metabolism, and amino acid, ion and vitamin transport and metabolism. The COG clustering of the differential expressed genes revealed that these five functional clusters were the most affected by SpoIIE’s inactivity after the sporulation cluster (Figure 7C).
Impact of the spoIIE Disruption on the Expression of Sporulation Genes
The genes involved in the sporulation process were the most down-regulated in the ΔspoIIE mutant. Indeed 84% of them are differentially expressed at least at one time point, and they have the lowest log2 fold change at 11 h and 23 h (Supplementary Table S5). To better assess the impact of spoIIE’s disruption, the genes involved in sporulation regulation were clustered based on the genes present in C. beijerinckii’s genome and on previous studies on sporulation in C. acetobutylicum (
The expression levels of most genes involved in the initiation of the sporulation and the activation of Spo0A by phosphorylation did not change in the ΔspoIIE mutant (Supplementary Table S6). The expression of spo0A remains rather stable. Only the expression of two genes (Cbei_4885 and Cbei_3375) encoding AbrB-type regulators significantly differed in the mutant compared to the WT. In C. acetobutylicum, abrB genes were shown to have a crucial role in the transition from acidogenesis to solventogenesis (Xue et al., 2016). Out of the three abrB genes identified in C. acetobutylicum, two (Cac_3647 and Cac_0310) encode transcriptional regulators down-regulating genes necessary for solventogenesis: the sol operon as well as ald, bdh, adhE2 and adc genes. In C. beijerinckii, five abrB genes, including Cbei_3375 and Cbei_4885, are annotated. Cbei_4885 is highly similar to Cac_3647 and Cac_0310 with 88 and 89% identity, respectively. The expression of Cbei_4885 increased at 11 h and decreased at 23 h (log2 fold change of 2.09 and −2.75) while the expression of Cbei_3375, 54% identity with Cac_3647 and Cac_0310, increases slightly at 23 h (log2 fold change of 1.86).
In the ΔspoIIE strain, expression levels of the sporulation genes belonging to the Spo0A regulon were similar to those of the WT strain. As expected, no mRNA of the spoIIE gene was detected in the mutant cells. The expression level of spoIIGA and sigE was not differentially expressed in the mutant, at 4 h log2 fold change is low (−2.69 and −2.60) but was not statistically relevant, padj > 0.05. At 11 h and 23 h, their expression returned to WT levels. The spoIIAA-spoIIAB-sigF operon was also slightly less expressed at 4 h in the mutant compared to the WT, as observed for the spoIIGA-sigE operon. As seen in other spore formers (
By contrast, the expression of genes involved in the later stages of the sporulation was more severely impacted by spoIIE’s disruption. The genes belonging to the σF regulon were all down-regulated at the three time points except spoIIR, suggesting that spoIIR is not part of the σF regulon in C. beijerinckii, in contrary to its homolog in B. subtilis. Also, the fold-change of expression varied throughout the σF regulon. Indeed, while some genes like spoIIQ, spoIIP and spoIVB were strongly down-regulated at the three time points (Supplementary Table S6), other genes like gpr and dacF were mainly down-regulated at 4 h and 23 h. In the case of dacF, this could be due to the low coverage of the gene already in the WT at 11 h, but that was not the case for gprR. This result suggests that another transcription factor might be involved, enabling the expression of these genes during exponential growth in the spoIIE mutant. Most genes involved later in the sporulation process (Stage III to VII) belonging to the σE, σG, or σK regulon, were down-regulated. Only four genes (spoVD1, spoVD2, spoVE and yzbD), associated with the σE regulon in B. subtilis, were not differentially expressed in our strain.
However, out of the 50 most differentially expressed genes at 11 h and 23 h (Supplementary Table S5), only 14% and 24% of genes correspond to known proteins directly involved in the sporulation process, while 54% and 58% of the genes, respectively were of unknown function. This large number of genes with unknown function shows that other genes might be involved in the sporulation process and that spoIIE inactivation and the cells block at stage II of the sporulation cycle has an impact on other gene clusters.
The Disruption of spoIIE Affects Other Cellular Processes
As observed in the COG clustering of the genes three times differentially expressed, the inactivation of spoIIE has an impact on several cellular processes besides sporulation, including the central metabolism (Figure 7C). Indeed the results that were obtained in this study showed the differential expression of genes associated with five functional groups: stress response, cell wall and membrane formation, signal transduction mechanisms and motility, carbohydrate metabolism, and amino acid, ions, and vitamin transport and metabolism.
Stress response
Genes involved in stress response were also strongly differentially expressed in the ΔspoIIE mutant (Supplementary Table S7). The genes coding for rubrerythrins (Cbei_0569, Cbei_2325, and Cbei_3257) and superoxide dismutase (Cbei_1507, Cbei_1856) belonging to the σG regulon in C. difficile (
Cell wall and membrane formation
As shown in the microscopy pictures (Figures 4, 5 and Supplementary Figures S2A,B), the morphology of ΔspoIIE mutant cells is very different from WT cells’. Thus, a difference in the expression of genes encoding membrane- and cell wall-associated proteins (Supplementary Table S7) was expected. Among the 243 genes differentially expressed in the mutant at the three time points, 42 belonged to the cell wall/membrane cluster. These 42 genes were all down-regulated in the mutant (log2 fold change ranging from −12 to −1.9). Several of them encode glycosyltransferases, which might suggest a possible modification of the cell wall. An up-regulation of the genes involved in septum formation, fstZ, minD and refZ was also observed. These genes were proven to be crucial for the positioning of the septum formation and sporulation in B. subtilis (
Signal transduction mechanisms and motility
The initiation of the sporulation cycle is triggered by several environmental signals (
About 20% of the genes involved in signal transduction are differentially expressed at least one time point in the ΔspoIIE mutant. Several genes encoding serine/threonine kinases and PAS/PAC sensor hybrid histidine kinases were differentially expressed. In particular, six genes were differentially expressed at all three time points. Genes encoding two serine/threonine kinases and an unknown gene were consistently down-regulated (log2 fold change between −10 and −3.8). Three genes belonging to the same operon were up-regulated (log2 fold change between 2.7 and 4.2). This operon is unique to C. beijerinckii. It harbors a PAS/PAC sensor hybrid histidine kinase, a diguanylate cyclase/phosphodiesterase with PAS/PAC sensor(s) and a diguanylate cyclase with a transcription regulator.
In Clostridium, the loss of motility is usually coupled with the transition from exponential phase to stationary phase (
Carbohydrate metabolism
Granulose biosynthesis
Granulose, also called bacterial glycogen, is a starch-like polymer produced at the onset of the sporulation cycle (
FIGURE 8

Differential expression at 23 h in the mutant of genes involved in (A) granulose formation and (B) granulose degradation figures adapted from Microcyc. Gene names are in colored boxes according to the differential expression in the ΔspoIIE mutant compared to the WT; genes in gray boxes have no change in expression, genes in green boxes are up-regulated, and genes in red boxes are down-regulated.
Central metabolism
Most of the genes involved in glycolysis were not differentially expressed in the ΔspoIIE mutant (Supplementary Table S9). Only one of the two copies of the pyruvate kinase, Cbei_1412, was down-regulated at the three time points, but no change in expression was observed for the second copy Cbei_4851. However, half of the genes involved in acidogenesis and solventogenesis were differentially expressed in the mutant compared to the WT strain at 23 h (Figure 9). Indeed, while the genes putatively involved in lactate and ethanol production were down-regulated in the mutant strain, the genes involved in butyrate and butanol production were up-regulated. However, no change in the expression of the genes coding for the proteins involved in acetate or acetone production was detected. The up-regulation of the sol operon, as well as genes encoding several alcohol dehydrogenases, might be linked to the down-regulation of the abrB gene, Cbei_4885, as a similar control was observed in C. acetobutylicum when its homolog was disrupted. We saw no differential expression of the ctfA/B operon, which could be linked to the accumulation of butyrate observed in the fermentation broth after 23 h, despite the up-regulation of the other genes involved in the butanol production.
FIGURE 9

Differential expression in the mutant at 23 h of genes involved in acids and solvent formation, figure adapted from Microcyc. Gene names are in colored boxes according to the differential expression in the ΔspoIIE mutant compared to the WT; genes in gray boxes have no change in expression, genes in green boxes are up-regulated and genes in red are down-regulated.
Amino acid, ions and vitamin transport and metabolism
Several vitamins and ions are essential for solventogenesis and sporulation in Clostridium (
At 4 h genes involved in ascorbate, cobalt and iron transport were down-regulated while genes encoding riboflavin and vitamin D transporters were upregulated. Moreover, at 11 h and 23 h, genes coding for cysteine-, glutamine- and manganese-transporters and the metabolism of these compounds were strongly down-regulated in the mutant (log2 fold change between −11 and −4). By contrast, we observed an up-regulation of genes encoding a proline transporter (Cbei_2870 and Cbei_2871) that has been linked to a rise in proline in the cell, which may act as a stress protectant (Takagi, 2008;
At 23 h, the genes linked to the transport of phosphate were up-regulated. By contrast, the genes coding for methionine and glycine/betaine transporters (log2 fold change from −5 to −4) were strongly down-regulated as well as the genes encoding for ascorbate- and iron transporters. The up-regulation of ferritin, coupled with the down-regulation of iron transporters, indicates an accumulation of iron in the cell, requiring the action of ferritin and an interruption of the iron uptake. The role of iron in sporulation in Clostridium is not clear. In C. botulinum, iron is needed for the formation of heat resistant spores while in C. sporogenes, its addition impairs sporulation (
Genes of unknown function
Our transcriptomic analysis highlighted the presence of a large cluster of genes of unknown function impacted by spoIIE’s inactivation. Indeed, 45% of the genes differentially expressed at the three times points did not have a known function. Except for one gene coding for a histidine triad (HIT) protein up-regulated at 23 h, these genes were down-regulated at all time points. Moreover, when focusing only on the genes differentially expressed at the 23 h, 56% of the most differentially expressed genes (log2 fold change above 5) belonged to the unknown function cluster. Some of these strongly down-regulated genes might code for spore coat proteins. Indeed, some coat proteins are species-specific and of small size, and they remain to be identified in C. beijerinckii (
Discussion
Clostridium beijerinckii is a solventogenic bacterium that is well-studied because of its potential to produce biofuels and biochemicals (
Furthermore, the morphology of ΔspoIIE mutants in both clostridia is also different. In contrast to C. acetobutylicum ΔspoIIE cells, C. beijerinckii ΔspoIIE cells were elongated with several septa and showed phase-dark bodies at their polar ends. C. acetobutylicum mutant cells do not display any differentiation phenotype. Indeed C. acetobutylicumΔspoIIE cells looked like vegetative cells; no internal structures apart from symmetrical septal membranes or other changes in morphology were described, even after 72 h of cultivation. However, in B. subtilis, the described ΔspoIIE mutant has a morphology close to the one observed in our study (
FIGURE 10

Schematic representation of the differences in cell morphology of ΔspoIIE mutants in B. subtilis, C. beijerinckii and C. acetobutylicum. The figure was constructed from data reported in this study and in
To better understand the role of SpoIIE, we compared the expression profile between the spoIIE mutant and WT strains by RNA sequencing at early-, mid-exponential and stationary phases. Forty percent of the total CDS were differentially expressed between both strains, with the 23 h time point showing the most changes. The genes involved in the regulation of sporulation after stage II were down-regulated after 23 h of culture. In contrast to what was observed in C. acetobutylicum, sigF, sigE and sigG were still expressed at 11 h. In C. acetobutylicum, a significant decrease of sigF transcript and σF production in the ΔspoIIE mutant was observed (
The sigE and sigG transcripts were slightly less abundant in the mutant than in the WT. These results suggest that in C. beijerinckii, the expression of sigE and sigG is only partially regulated by σF. Despite being expressed, σE and σG might not be functional since the genes belonging to their respective regulons were weakly expressed in the mutant. This observation could be explained by (i) the need for σE and σG to reach a certain threshold level in the cell for them to allow the transcription of their regulons, (ii) a lack of active σE in the mother cell. Pro-σE might not be processed in the mutant and σE might then stay inactive. In B. subtilis, the SpoIIGA protease processes Pro-σE, after being activated by SpoIIR, which is expressed under the control of σF. It is interesting to note that in C. beijerinckii, we failed to detect a clear impact of spoIIE inactivation on spoIIGA and spoIIR expression, even though spoIIR belongs to the σF regulon in B. subtilis. Thus, another protein, belonging to the σF regulon might be required to process Pro-σE. It is worth noting that also in other clostridia, the expression of spoIIR is not strictly dependent on σF (
The transcriptomic data confirmed that next to sporulation, other stationary-phase phenomena such as stress response, signal transduction, motility, carbohydrate metabolism, and the transport of amino acids, vitamins and ions were modified in the ΔspoIIE mutant. The interruption of the sporulation cycle led to an extension of the metabolic activity of the cells. Thus, the asporogenous ΔspoIIE mutant displayed three exciting features that, as observed in other solventogenic clostridia, appear to be associated with an increase in solvent production. Firstly, several genes involved in butyrate and butanol production were up-regulated. While a definite rise in acid titer was observed in the medium at the end of the fermentation, no substantial increase in butanol was measured. This rise in acid titer might have been caused by a lack of CtfA/B proteins. The CtfA/B complex enables the reassimilation of acids; however, their genes were not up-regulated in the ΔspoIIE mutant. The introduction of these genes in the ΔspoIIE mutant, expressed under the control of a constitutive promotor, might enhance the solvent production, as previously observed in C. acetobutylicum (
In conclusion, this study has enabled us to identify genes potentially involved in the regulation of stationary phase phenomena. The expression of genes involved in secondary metabolism and signaling pathways, as well as a large number of genes of unknown functions, were strongly impacted in the spoIIE mutant. These genes should be investigated in more detail to obtain a better insight into the regulation of sporulation and other stationary-phase events. Indeed, the complete sporulation regulation pathway has not yet been elucidated; critical information such as activation of σG (
The impact of the interruption of the sporulation cycle at Stage II in C. beijerinckii was evaluated through fermentation, microscopy, genome and transcriptome analysis. The fermentation potential of the mutant strain could be evaluated. The transcriptomic analysis provided directions to better exploit this potential, by gene engineering or specific changes in media composition. Furthermore, this study reveals the complexity of the regulation of sporulation in Clostridium and its interconnection with other cellular regulatory networks. While some features are conserved, others seem to differ even between solventogenic clostridia. Nevertheless, this work constitutes a solid basis for further investigation of the molecular regulation of sporulation in C. beijerinckii and its potential for industrial application.
Statements
Data availability statement
The RNA seq data generated and analyzed for this study have been deposited in the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-7481. The genome sequence of Clostridium beijerinckii NCIMB 8052 ΔspoIIE strain described in this study is available on the European Nucleotide Archive (ENA) under the accession number PRJEB39199.
Author contributions
MD co-designed the study, performed the experiments, collected and analyzed the data, and wrote the manuscript. NK collected the transcriptomic data. MM contributed to the transcriptomic data analysis. FC supervised the experimental work. IM-V, JO and SK contributed to data interpretation, discussions and revised the manuscript. AMLC co-designed the study, supervised the work, contributed to data interpretation, discussions and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The work presented in this article was financed by the European Union Marie Sklodowska Curie Innovative Training Networks (ITN) – Contract number 642068.
Acknowledgments
The authors thank the Biomics Platform, C2RT, at the Institut Pasteur supported by France Génomique (ANR-10-INBS-09-09) and IBISA for sequencing the RNA-Seq library, Dr. François Wasels for his help on the genome analysis, Dr. Karel Sedlar for his advice on RNA seq data analysis, Arjen Bader for technical assistance with the fluorescence microscopy analysis and Dr. Eric M. Ransom for the gift of pRAN73 plasmid.
Conflict of interest
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.556064/full#supplementary-material
FIGURE S1Complementation of the ΔspoIIE mutant restores sporulation (A) Phase-contrast microscopy pictures (x400) after 35 h of cultivation of the ΔspoIIE mutant and the complemented strain; (B) Growth on plates after heat-shock treatment of 48 h old ΔspoIIE and complemented ΔspoIIE cultures. The short dark blue arrows indicate mature spores.
FIGURE S2Additional fluorescence microscopy images of wild type and ΔspoIIE mutant cells; (A) Images of several WT and ΔspoIIE mutant cells after 20 h of cultivation and stained by the membrane staining MTG; (B) Image of two ΔspoIIE mutant cells after 20 h of cultivation and stained by the membrane staining MTG; (C) Fluorescence images of the ΔspoIIE mutant cells stained by MTG and harboring either the mCherry empty plasmid (pRAN73S) or the plasmid expressing mCherry fused to Cbei_0097 (pRAN73S:cbei_0097) at 18 h of cultivation (after 8 hours of atc induction). The white arrows indicate septa observed in the mutant strain.
TABLE S1Oligonucleotides used in this study; the restriction sites are underlined.
TABLE S2Sequences of primers used in Q-RT-PCR experiments for the validation of the RNAseq analysis and Q-RT-PCR values.
TABLE S3SNPs detected in each ΔspoIIE mutant.
TABLE S4Genes significantly differentially expressed in the ΔspoIIE mutant at the three time points and their expression levels in the mutant relative to the wild type in log2 fold change.
TABLE S5Most differentially expressed genes at 11 h and 23 h and their expression levels in the mutant relative to the wild type in log2 fold change.
TABLE S6Differential expression of the genes belonging to the sporulation cluster in log2 fold change; the genes in italic were not significantly differentially expressed (| log2 fold| < 1.5 or padj > 0.05) during the fermentation.
TABLE S7Differential expression of the genes involved in stress response, cell wall/membrane composition, chemotaxis and motility in the ΔspoIIE mutant relative to the wild type, in log2 fold change; only the genes for which | log2 fold| > 1.5 and padj ≤ 0.05 at least once during the fermentation.
TABLE S8Differential expression of the PKS-NRPS gene cluster in the ΔspoIIE mutant relative to the wild type, in log2 fold change; the genes in italic were not significantly differentially expressed (| log2 fold| < 1.5 or padj > 0.05) during the fermentation.
TABLE S9Differential expression of the genes involved in the granulose formation, glycolysis, acidogenesis and solventogenesis in the spoIIE mutant relative to the wild type, in log2 fold change; the genes in italic were not significantly differentially expressed (| log2 fold| < 1.5 or padj > 0.05) during the fermentation.
TABLE S10Differential expression of the genes involved in the amino acid, ion and vitamin transport and metabolism in the spoIIE mutant relative to the wild type, in log2 fold change; only the genes for which | log2 fold| > 1.5 and padj ≤ 0.05 at least once during the fermentation.
TABLE S11Differential expression at 23 h of small species-specific genes of unknown function down-regulated in the ΔspoIIE mutant relative to the wild type, in log2 fold change; genes in bold are down-regulated at all time points.
TABLE S12Fermentations characteristics of WT, ΔspoIIE and ΔspoIIE pSpoIIE complemented mutant after 70 h of cultivation in CM2 media in serum bottles at 35°C, one standard deviation of the mean was determined based on the data from biological replicates (n = 2).
Footnotes
References
1
Al-HinaiM. A.JonesS. W.PapoutsakisE. T. (2014). σK of Clostridium acetobutylicum is the first known sporulation-specific sigma factor with two developmentally separated roles, one early and one late in sporulation.J. Bacteriol.196287–299. 10.1128/JB.01103-13
2
Al-HinaiM. A.JonesS. W.PapoutsakisE. T. (2015). The Clostridium sporulation programs: diversity and preservation of endospore differentiation.Microbiol. Mol. Biol. Rev.7919–37. 10.1128/MMBR.00025-14
3
BarákI.MuchováK. (2018). The positioning of the asymmetric septum during sporulation in Bacillus subtilis. Edited by Eric Cascales.PLoS One13:e0201979. 10.1371/journal.pone.0201979
4
BarákI.MuchováK.LabajováN. (2019). Asymmetric cell division during Bacillus subtilis sporulation.Future Microbiol.14353–363. 10.2217/fmb-2018-0338
5
BarákI.YoungmanP. (1996). SpoIIE mutants of Bacillus subtilis comprise two distinct phenotypic classes consistent with a dual functional role for the SpoIIE protein.J. Bacteriol.1784984–4989. 10.1128/JB.178.16.4984-4989.1996
6
BiC.JonesS. W.HessD. R.TracyB. P.PapoutsakisE. T. (2011). SpoIIE is necessary for asymmetric division, sporulation, and expression of σF, σE, and σG but does not control solvent production in Clostridium acetobutylicum ATCC 824.J. Bacteriol.1935130–5137. 10.1128/JB.05474-11
7
BrownE. E.MillerA. K.KriegerI. V.OttoR. M.SacchettiniJ. C.HermanJ. K. (2019). A DNA-Binding protein tunes septum placement during Bacillus subtilis sporulation.J. Bacteriol.201:e00287-19. 10.1128/jb.00287-19
8
CollasF.KuitW.ClementB.MarchalR.Lopez-ContrerasA. M.MonotF. (2012). Simultaneous production of isopropanol, butanol, ethanol and 2,3-butanediol by Clostridium acetobutylicum ATCC 824 engineered Strains.AMB Express2:45. 10.1186/2191-0855-2-45
9
de GérandoH. M.Fayolle-GuichardF.RudantL.MillahS. K.MonotF.Lopes FerreiraN.et al (2016). Improving isopropanol tolerance and production of Clostridium beijerinckii DSM 6423 by random mutagenesis and genome shuffling.Appl. Microbiol. Biotechnol.1005427–5436. 10.1007/s00253-016-7302-5
10
de HoonM. J. L.EichenbergerP.VitkupD. (2010). Hierarchical evolution of the bacterial sporulation network.Curr. Biol.20R735–R745. 10.1016/j.cub.2010.06.031
11
DialloM.HocqR.CollasF.ChartierG.WaselsF.WijayaH. S.et al (2020). Adaptation and application of a by two-Plasmid inducible CRISPR-Cas9 system in Clostridium beijerinckii.Methods17251–60. 10.1016/j.ymeth.2019.07.022
12
EdwardsA. N.NawrockiK. L.McBrideS. M. (2014). Conserved oligopeptide permeases modulate sporulation initiation in Clostridium difficile.Infect. Immun.824276–4291. 10.1128/IAI.02323-14
13
EydallinG.VialeA. M.Morán-ZorzanoM. T.MuñozF. J.MonteroM.Baroja-FernándezE.et al (2007). Genome-wide screening of genes affecting glycogen metabolism in Escherichia coli K-12.FEBS Lett.5812947–2953. 10.1016/j.febslet.2007.05.044
14
FeuchtA.AbbottsL.ErringtonJ. (2002). The cell differentiation protein SpoIIE contains a regulatory site that controls its phosphatase activity in response to asymmetric septation.Mol. Microbiol.451119–1130. 10.1046/j.1365-2958.2002.03082.x
15
FimlaidK. A.BondJ. P.SchutzK. C.PutnamE. E.LeungJ. M.LawleyT. D.et al (2013). Global analysis of the sporulation pathway of Clostridium difficile. Edited by Patrick H. Viollier.PLoS Genet.9:e1003660. 10.1371/journal.pgen.1003660
16
FimlaidK. A.ShenA. (2015). Diverse mechanisms regulate sporulation sigma factor activity in the Firmicutes.Curr. Opin. Microbiol.2488–95. 10.1016/j.mib.2015.01.006
17
FujitaM.González-PastorJ. E.LosickR. (2005). High- and low-threshold genes in the Spo0A regulon of Bacillus subtilis.J. Bacteriol.1871357–1368. 10.1128/JB.187.4.1357-1368.2005
18
GalperinM. Y.MekhedovS. L.PuigboP.SmirnovS.WolfY. I.RigdenD. J. (2012). Genomic determinants of sporulation in Bacilli and Clostridia: towards the minimal set of sporulation-specific genes.Environ. Microbiol.142870–2890. 10.1111/j.1462-2920.2012.02841.x
19
GrageK.McDermottP.RehmB. H. A. (2017). Engineering Bacillus megaterium for production of functional intracellular materials.Microb. Cell Fact.16:211. 10.1186/s12934-017-0823-5
20
GuY.JiangY.YangS.JiangW. (2014). Utilization of economical substrate-derived carbohydrates by solventogenic clostridia: pathway dissection, regulation and engineering.Curr. Opin. Biotechnol.29124–131. 10.1016/j.copbio.2014.04.004
21
GutierrezN. A.MaddoxI. S. (1987). Role of chemotaxis in solvent production by Clostridium acetobutylicum.Appl. Environ. Microbiol.531924–1927. 10.1128/aem.53.8.1924-1927.1987
22
HermanN. A.KimS. J.LiJ. S.CaiW.KoshinoH.ZhangW. (2017). The industrial anaerobe Clostridium acetobutylicum uses polyketides to regulate cellular differentiation.Nat. Commun.8:1514. 10.1038/s41467-017-01809-5
23
JamroskovicJ.ChromikovaZ.ListC.BartovaB.BarákI.Bernier-LatmaniR. (2016). Variability in DPA and calcium content in the spores of Clostridium species.Front. Microbiol.7:1791. 10.3389/fmicb.2016.01791
24
JonesD. T.WoodsD. R. (1986). Acetone-Butanol fermentation revisited.Microbiol. Rev.50484–524. 10.1128/mmbr.50.4.484-524.1986
25
JonesS. W.TracyB. P.GaidaS. M.PapoutsakisE. T. (2011). Inactivation of σF in Clostridium acetobutylicum ATCC 824 blocks sporulation prior to asymmetric division and abolishes σE and σG protein expression but does not block solvent formation.J. Bacteriol.1932429–2440. 10.1128/JB.00088-11
26
KempermanR.JonkerM.NautaA.KuipersO. P.KokJ. (2003). Functional analysis of the gene cluster involved in production of the bacteriocin circularin a by Clostridium beijerinckii ATCC 25752.Appl. Environ. Microbiol.695839–5848. 10.1128/aem.69.10.5839-5848.2003
27
KotteA. K.SevernO.BeanZ.SchwarzK.MintonN. P.WinzerK. (2017). RNPP-type quorum sensing regulates solvent formation and sporulation in Clostridium acetobutylicum.bioRxiv44:106666. 10.1101/106666
28
KroghA.LarssonB.Von HeijneG.SonnhammerE. L. L. (2001). Predicting transmembrane protein topology with a hidden markov model: application to complete genomes.J. Mol. Biol.305567–580. 10.1006/jmbi.2000.4315
29
KuitW.MintonN. P.López-ContrerasA. M.EgginkG. (2012). Disruption of the acetate kinase (Ack) gene of Clostridium acetobutylicum results in delayed acetate production.Appl. Microbiol. Biotechnol.94729–741. 10.1007/s00253-011-3848-4
30
LangmeadB.TrapnellC.PopM.SalzbergS. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome.Genome Biol.10:R25. 10.1186/gb-2009-10-3-r25
31
LeeS.LeeJ. H.MitchellR. J. (2015). Analysis of Clostridium beijerinckii NCIMB 8052’s transcriptional response to ferulic acid and its application to enhance the strain tolerance.Biotechnol.Biofuels8:68.
32
Lépiz-AguilarL.Rodríguez-RodríguezC. E.AriasM. L.LutzG.UlateW. (2011). Butanol production by Clostridium beijerinckii BA101 using cassava flour as fermentation substrate: enzymatic versus chemical pretreatments.World J. Microbiol. Biotechnol.271933–1939. 10.1007/s11274-010-0630-1
33
LetzelA. C.PidotS. J.HertweckC. (2013). A genomic approach to the cryptic secondary metabolome of the anaerobic world.Nat. Prod. Rep.30392–428. 10.1039/c2np20103h
34
LevinP. A.LosickR. (1994). Characterization of a cell division gene from Bacillus subtilis that is required for vegetative and sporulation septum formation.J. Bacteriol.1761451–1459. 10.1128/JB.176.5.1451-1459.1994
35
LiD.MengC.WuG.XieB.HanY.GuoY.et al (2018). Effects of zinc on the production of alcohol by Clostridium carboxidivorans P7 using model syngas.J. Ind. Microbiol. Biotechnol.4561–69. 10.1007/s10295-017-1992-2
36
LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN.et al (2009). The sequence alignment/map format and SAMtools.Bioinformatics252078–2079. 10.1093/bioinformatics/btp352
37
LiS.HuangL.KeC.PangZ.LiuL. (2020). Pathway dissection, regulation, engineering and application: lessons learned from biobutanol production by solventogenic clostridia.Biotechnol. Biofuels13:39. 10.1186/s13068-020-01674-3
38
LiT.YanY.HeJ. (2014). Reducing cofactors contribute to the increase of butanol production by a wild-type Clostridium Sp. strain BOH3.Bioresour. Technol.155220–228. 10.1016/j.biortech.2013.12.089
39
LiaoZ.ZhangY.LuoS.SuoY.ZhangS.WangJ. (2017). Improving cellular robustness and butanol titers of Clostridium acetobutylicum ATCC824 by introducing heat shock proteins from an extremophilic bacterium.J. Biotechnol.2521–10. 10.1016/j.jbiotec.2017.04.031
40
ListC.HosseiniZ.Lederballe MeibomK.HatzimanikatisV.Bernier-LatmaniR. (2019). Impact of iron reduction on the metabolism of Clostridium acetobutylicum.Environ. Microbiol.213548–3563. 10.1111/1462-2920.14640
41
LouieP.LeeA.StansmoreK.GrantR.GintherC.LeightonT. (1992). Roles of RpoD, SpoIIF, SpoIIJ, SpoIIN, and Sin in regulation of Bacillus subtilis Stage II sporulation-specific transcription.J. Bacteriol.1743570–3576. 10.1128/JB.174.11.3570-3576.1992
42
MahJ. H.KangD. H.TangJ. (2008). Effects of minerals on sporulation and heat resistance of Clostridium sporogenes.Int. J. Food Microbiol.128385–389. 10.1016/j.ijfoodmicro.2008.10.002
43
MaitiS.GallasteguiG.SarmaS. J.BrarS. K.Le BihanY.DroguiP.et al (2016). A re-look at the biochemical strategies to enhance butanol production.Biomass Bioenergy94187–200. 10.1016/j.biombioe.2016.09.001
44
MieleV.PenelS.DuretL. (2011). Ultra-Fast sequence clustering from similarity networks with SiLiX.BMC Bioinformatics12:116. 10.1186/1471-2105-12-116
45
MonotM.OrgeurM.CamiadeE.BrehierC.DupuyB. (2014). COV2HTML: a visualization and analysis tool of bacterial Next Generation Sequencing (NGS) data for postgenomics life scientists.OMICS18184–195. 10.1089/omi.2013.0119
46
MonteroM.EydallinG.VialeA. M.AlmagroG.MuñozF. J.RahimpourM.et al (2009). Escherichia coli glycogen metabolism is controlled by the PhoP-PhoQ regulatory system at submillimolar environmental Mg2+ concentrations, and is highly interconnected with a wide variety of cellular processes.Biochem. J.424129–141. 10.1042/BJ20090980
47
MuchováK.ChromikováZ.BarákI. (2020). Linking the peptidoglycan synthesis protein complex with asymmetric cell division during Bacillus subtilis sporulation.Int. J. Mol. Sci.21:4513. 10.3390/ijms21124513
48
MuchováK.ChromikováZ.BradshawN.WilkinsonA. J.BarákI. (2016). Morphogenic protein Rodz interacts with sporulation specific SpoIIE in Bacillus Subtilis. Edited by Claude Prigent.PLoS One11:e0159076. 10.1371/journal.pone.0159076
49
NimbalkarP. R.KhedkarM. A.ParulekarR. S.ChandgudeV. K.SonawaneK. D.ChavanP. V.et al (2018). Role of trace elements as cofactor: an efficient strategy toward enhanced biobutanol production.ACS Sustain. Chem. Eng.69304–9313. 10.1021/acssuschemeng.8b01611
50
ParedesC. J.AlsakerK. V.PapoutsakisE. T. (2005). A comparative genomic view of clostridial sporulation and physiology.Nat. Rev. Microbiol.3969–978. 10.1038/nrmicro1288
51
Paredes-SabjaD.ShenA.SorgJ. A. (2014). Clostridium difficile spore biology: sporulation, germination, and spore structural proteins.Trends Microbiol.22406–416. 10.1016/j.tim.2014.04.003
52
PatakovaP.LinhovaM.RychteraM.PaulovaL.MelzochK. (2013). Novel and neglected issues of Acetone-Butanol-Ethanol (ABE) fermentation by Clostridia: Clostridium metabolic diversity, tools for process mapping and continuous fermentation systems.Biotechnol. Adv.3158–67. 10.1016/j.biotechadv.2012.01.010
53
PfafflM. W. (2001). A new mathematical model for relative quantification in Real-Time RT-PCR.Nucleic Acids Res.29e45.
54
QuanJ.TianJ. (2014). Circular polymerase extension cloning.Methods Mol. Biol.1116103–117. 10.1007/978-1-62703-764-8_8
55
RansomE. M.EllermeierC. D.WeissD. S. (2015). Use of mCherry red fluorescent protein for studies of protein localization and gene expression in Clostridium difficile.Appl. Environ. Microbiol.811652–1660. 10.1128/AEM.03446-14
56
RansomE. M.WeissD. S.EllermeierC. D. (2016). Use of mCherryOpt fluorescent protein in Clostridium difficile.Methods Mol. Biol.147653–67. 10.1007/978-1-4939-6361-4_5
57
ReysenbachA. L.RavenscroftN.LongS.JonesD. T.WoodsD. R. (1986). Characterization, biosynthesis, and regulation of granulose in Clostridium acetobutylicum.Appl. Environ. Microbiol.52185–190. 10.1128/aem.52.1.185-190.1986
58
RobsonR. L.RobsonR. M.MorrisJ. G. (1974). The biosynthesis of granulose by Clostridium pasteurianum.Biochem. J.144503–511. 10.1042/bj1440503
59
SaujetL.PereiraF. C.HenriquesA. O.Martin-VerstraeteI. (2014). The regulatory network controlling spore formation in Clostridium difficile.FEMS Microbiol. Lett.3581–10. 10.1111/1574-6968.12540
60
SaujetL.PereiraF. C.SerranoM.SoutourinaO.MonotM.ShelyakinP. V.et al (2013). Genome-Wide analysis of cell type-specific gene transcription during spore formation in Clostridium difficile. Edited by Patrick H. Viollier.PLoS Genet.9:e1003756. 10.1371/journal.pgen.1003756
61
SchaefferA. B.FultonM. D. (1933). A simplified method of staining endospores.Science77:194. 10.1126/science.77.1990.194
62
ScotcherM. C.BennettG. N. (2005). SpoIIE regulates sporulation but does not directly affect solventogenesis in Clostridium acetobutylicum ATCC 824.J. Bacteriol.1871930–1936. 10.1128/JB.187.6.1930-1936.2005
63
SteinerE.DagoA.YoungD. I.HeapJ. T.MintonN. P.HochJ. A.et al (2011). Multiple orphan histidine kinases interact directly with Spo0A to control the initiation of endospore formation in Clostridium acetobutylicum.Mol. Microbiol.80641–654. 10.1111/j.1365-2958.2011.07608.x
64
SteinerE.ScottJ.MintonN. P.WinzerK. (2012). An Agr quorum sensing system that regulates granulose formation and sporulation in Clostridium acetobutylicum.Appl. Environ. Microbiol.781113–1122. 10.1128/AEM.06376-11
65
TakagiH. (2008). Proline as a stress protectant in yeast: physiological functions, metabolic regulations, and biotechnological applications.Appl. Microbiol. Biotechnol.81211–223. 10.1007/s00253-008-1698-5
66
TangH.BomhoffM. D.BrionesE.ZhangL.SchnableJ. C.LyonsE. (2015). SynFind: compiling syntenic regions across any set of genomes on demand.Genome Biol. Evol.73286–3298. 10.1093/gbe/evv219
67
TomasC. A.BeamishJ.PapoutsakisE. T. (2004). Transcriptional analysis of butanol stress and tolerance in Clostridium acetobutylicum.J. Bacteriol.1862006–2018. 10.1128/JB.186.7.2006-2018.2004
68
TomasC. A.WelkerN. E.PapoutsakisE. T. (2003). Overexpression of GroESL in Clostridium acetobutylicum results in increased solvent production and tolerance, prolonged metabolism, and changes in the cell’s transcriptional program.Appl. Environ. Microbiol.694951–4965. 10.1128/AEM.69.8.4951-4965.2003
69
TracyB. P.JonesS. W.FastA. G.IndurthiD. C.PapoutsakisE. T. (2012). Clostridia: the importance of their exceptional substrate and metabolite diversity for biofuel and biorefinery applications.Curr. Opin. Biotechnol.23364–381. 10.1016/j.copbio.2011.10.008
70
TracyB. P.JonesS. W.PapoutsakisE. T. (2011). Inactivation of σE and σG in Clostridium acetobutylicum illuminates their roles in clostridial-cell-form biogenesis, granulose synthesis, solventogenesis, and spore morphogenesis.J. Bacteriol.1931414–1426. 10.1128/JB.01380-10
71
van der WalH.SperberM.Houweling-TanB.BakkerR. R. C.BrandenburgW.López-ContrerasA. M. (2013). Production of Acetone, Butanol, and Ethanol from biomass of the green seaweed Ulva lactuca.Bioresour. Technol.128431–437. 10.1016/j.biortech.2012.10.094
72
VaretH.Brillet-GuéguenL.CoppéeJ. Y.DilliesM. A. (2016). SARTools: a DESeq2- and EdgeR-Based R Pipeline for Comprehensive Differential Analysis of RNA-Seq Data.PLoS One11:e0157022. 10.1371/journal.pone.0157022
73
VasylkivskaM.JureckovaK.BranskaB.SedlarK.KolekJ.ProvaznikI.et al (2019). Transcriptional analysis of amino acid, metal ion, vitamin and carbohydrate uptake in butanol-producing Clostridium beijerinckii NRRL B-598. Edited by Zhiqiang Wen.PLoS One14:e0224560. 10.1371/journal.pone.0224560
74
WangS. T.SetlowB.ConlonE. M.LyonJ. L.ImamuraD.SatoT.et al (2006). The forespore line of gene expression in Bacillus subtilis.J. Mol. Biol.35816–37. 10.1016/j.jmb.2006.01.059
75
WangY.LiX.MaoY.BlaschekH. P. (2012). Genome-Wide dynamic transcriptional profiling in Clostridium beijerinckii NCIMB 8052 using single-nucleotide resolution RNA-Seq.BMC Genomics13:102. 10.1186/1471-2164-13-102
76
WilsonW. A.RoachP. J.MonteroM.Baroja-FernándezE.MuñozF. J.EydallinG.et al (2010). Regulation of glycogen metabolism in yeast and bacteria.FEMS Microbiol. Rev.34952–985. 10.1111/j.1574-6976.2010.00220.x
77
XueQ.YangY.ChenJ.ChenL.YangS.JiangW.et al (2016). Roles of three AbrBs in regulating two-phase Clostridium acetobutylicum fermentation.Appl. Microbiol. Biotechnol.1009081–9089. 10.1007/s00253-016-7638-x
78
ZhuB.StülkeJ. (2018). SubtiWiki in 2018: from genes and proteins to functional network annotation of the model organism Bacillus subtilis.Nucleic Acids Res.46D743–D748. 10.1093/nar/gkx908
Summary
Keywords
Clostridium beijerinckii NCIMB 8052, sporulation, spoIIE, ABE production, CRISPR-Cas9, RNA seq, transcriptome analysis
Citation
Diallo M, Kint N, Monot M, Collas F, Martin-Verstraete I, van der Oost J, Kengen SWM and López-Contreras AM (2020) Transcriptomic and Phenotypic Analysis of a spoIIE Mutant in Clostridium beijerinckii. Front. Microbiol. 11:556064. doi: 10.3389/fmicb.2020.556064
Received
27 April 2020
Accepted
20 August 2020
Published
15 September 2020
Volume
11 - 2020
Edited by
Patrick Eichenberger, New York University, United States
Reviewed by
Imrich Barak, Institute of Molecular Biology (SAS), Slovakia; Charles P. Moran, Emory University, United States
Updates

Check for updates
Copyright
© 2020 Diallo, Kint, Monot, Collas, Martin-Verstraete, van der Oost, Kengen and López-Contreras.
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: Mamou Diallo, mamou.diallo@gmx.net
This article was submitted to Microbial Physiology and Metabolism, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.