Transcriptome profiling analysis reveals metabolic changes across various growth phases in Bacillus pumilus BA06
- Lin-Li Han†1, 2,
- Huan-Huan Shao†1, 2,
- Yong-Cheng Liu1, 2,
- Gang Liu1, 2,
- Chao-Ying Xie1, 2,
- Xiao-Jie Cheng1, 2,
- Hai-Yan Wang1, 2,
- Xue-Mei Tan1, 2 and
- Hong Feng1, 2Email authorView ORCID ID profile
© The Author(s). 2017
Received: 1 March 2017
Accepted: 4 July 2017
Published: 11 July 2017
Bacillus pumilus can secret abundant extracellular enzymes, and may be used as a potential host for the industrial production of enzymes. It is necessary to understand the metabolic processes during cellular growth. Here, an RNA-seq based transcriptome analysis was applied to examine B. pumilus BA06 across various growth stages to reveal metabolic changes under two conditions.
Based on the gene expression levels, changes to metabolism pathways that were specific to various growth phases were enriched by KEGG analysis. Upon entry into the transition from the exponential growth phase, striking changes were revealed that included down-regulation of the tricarboxylic acid cycle, oxidative phosphorylation, flagellar assembly, and chemotaxis signaling. In contrast, the expression of stress-responding genes was induced when entering the transition phase, suggesting that the cell may suffer from stress during this growth stage. As expected, up-regulation of sporulation-related genes was continuous during the stationary growth phase, which was consistent with the observed sporulation. However, the expression pattern of the various extracellular proteases was different, suggesting that the regulatory mechanism may be distinct for various proteases. In addition, two protein secretion pathways were enriched with genes responsive to the observed protein secretion in B. pumilus. However, the expression of some genes that encode sporulation-related proteins and extracellular proteases was delayed by the addition of gelatin to the minimal medium.
The transcriptome data depict global alterations in the genome-wide transcriptome across the various growth phases, which will enable an understanding of the physiology and phenotype of B. pumilus through gene expression.
Bacillus pumilus is an endospore-forming, gram-positive, rod-shaped bacterium. Due to its metabolic diversity and spore dispersal, B. pumilus is ubiquitous in various environments and commonly resistant to extreme environmental conditions [1, 2]. Similar to other Bacillus species, B. pumilus is able to secrete a large number of industrial enzymes, such as lipases , xylanases , and proteases [5–7]. Therefore, B. pumilus has attracted attention in biotechnology and was selected to engineer novel industrial production strain [8, 9]. In addition, some strains of B. pumilus were used to produce valuable small molecules [10, 11] and have served as biocontrol agents to manage plant diseases . However, studies pertaining to the physiological and metabolic processes of B. pumilus are extremely limited in comparison with the other Bacillus species.
Along with the great advances in genomics, many strains of B. pumilus have been selected for genome sequencing. To date, more than 30 genomes of B. pumilus are available at NCBI. Based on genomic alignments, B. pumilus is closer to B. subtilis, B. licheniformis, and B. amyloliquefaciens , suggesting that similar metabolic or physiological processes may exist between B. pumilus and the model organism B. subtilis. Although the genome is considered the blueprint of life, much information regarding the physiological or metabolic processes is not directly accessible from the genome . Therefore, various omics technologies such as transcriptomics, proteomics and metabolomics, have been employed as essential steps toward gaining insights into cell physiology from the genome. For instance, a combined omics-based approach has been applied to B. pumilus to try to understand cell physiology in response to oxidative stress and protein secretion [15, 16]. For this purpose, RNA-seq -based transcriptomics analysis is one of the most powerful tools that not only provides important insights into the functional elements of the genome, gene expression patterns and regulation , but also offers a simpler and more cost-effective approach . Therefore, the RNA-seq method has been widely applied to many Bacillus species. For example, B. subtilis [19, 20], B. licheniformis , and B. thuringiensis  were examined by RNA-seq analysis.
B. pumilus BA06 has been isolated from proteinaceous soil and is demonstrated to be able to secrete extracellular proteases that exhibited great potential in leather processing [5, 23, 24]. However, the production of extracellular proteases is considered to occur at the stationary growth phase and is extensively regulated in many species of Bacillus [25, 26]. For example, nitrogen and carbon sources have great impacts on the production of extracellular proteases . Thus, medium components have usually been optimized for fermentation of proteases . In addition, the other physiological processes of secondary metabolite synthesis, protein secretion, and sporulation occur during the stationary growth phase [29–31]. Based on studies in B. subtilis, a transition point occurs between the exponential growth and stationary growth phases . Over the various growth phases, the expression of many genes may be turned off and another set of genes may be turned on. For example, more than 100 genes, whose expression was induced at the onset of the stationary growth phase, have been assigned to the SigB regulon [33, 34].
By taking advantage of the RNA-seq technology, a time-resolved transcriptomic analysis to cover the various growth stages of B. pumilus BA06 was performed. The results clearly indicated that changes in gene expression or metabolic pathways were occurring during the various growth phases; this observation will be helpful for us toward understanding the physiology and phenotype of B. pumilus.
Cell growth, sporulation and extracellular protease activity
Meanwhile, the extracellular protease activity was also monitored during the growing period. Fig. 1b shows that the extracellular protease activity mainly occurred at the onset of the stationary phase with a peak at 48 h in MM, which was consistent with previous results . Moreover, the addition of gelatin led to a higher peak at the delayed time (72 h) in extracellular protease activity (Fig. 1b), indicating that the nitrogen source (gelatin) affected the production of extracellular proteases.
RNA sequencing and identification of differentially expressed genes (DEGs)
Summary of RNA-seq and the reads mapped to the genome of B. pumilus BA06
Left size (M)
Right size (M)
Number of reads
Mapped reads (%)
Paired reads (%)
Top 12 functional classification of the DEGs of B. pumilus BA06 during the entire growth curve by RAST analysis
Amino Acids and Derivatives
Motility and Chemotaxis
Dormancy and Sporulation
Cell Wall and Capsule
Iron Acquisition and Metabolism
Nucleosides and Nucleotides
Changes in the metabolic pathways at the transition point of growth
Continuous up-regulation of sporulation-related genes
A master regulator protein, Spo0A, governs sporulation initiation through phosphorylation by a series of kinases and phosphotransferase . All of these genes were identified in B. pumilus BA06 (Additional file 5) and were clustered into group G-2 (Fig. 5). The expression profile of these genes was generally different from the G-1 group. Many of these genes (G-2) were repressed at the later growth phases, perhaps indicating their leading role only in sporulation initiation.
During the procedure of endospore formation, four sigma factors (SigE, SigF, SigG, and SigK) play important roles in B. subtilis . Their counterparts in B. pumilus BA06 were identified (group G-3), and all of their expression levels were continuously up-regulated across all growth phases (Fig. 5 and Additional file 5). In contrast, expression of sigH (peg.2382) fluctuated over the growth course.
Extracellular proteases and the protein secretion systems
Relative expression level of the extracellular proteases of B. pumilus BA06 during the entirety of each growth phase
Expression level in FPKM value
serine alkaline protease (AprE)
alkaline serine protease (AprX)
subtilisin Carlsberg (SubE)
Alkaline serine proteinase (Epr)
extracellular protease (Vpr)
Wall-associated protease (WprA)
Sigma factors and regulator proteins
The other sigma factors exhibited different patterns of expression (Additional file 7). For example, the expression of sigV (peg.3701), rpoN (peg.1161) and rpoE (peg.399) was up-regulated mainly at the stationary growth phase (48 h). In contrast, the expression of sigV (peg.2600) and sig70 (peg.775) was repressed upon entry into the transition phase and thereafter (Fig. 7a).
In B. subtilis, several regulatory proteins were recognized as transition-state regulators, such as Hpr (ScoC), AbrB, and SinR. All their homologs were encoded by the BA06 genome. Their expression pattern is shown in Fig. 7b (Additional file 7). The gene for hpr was highly expressed at the exponential growth phase (4 h) and then declined. In contrast, expression of SinR and pai2 was increased at 48 h and later, especially in GM. Another transition-state regulator, AbrB, was strongly repressed at 48 h, which is consistent with a similar observation in B. subtilis .
Global regulators play important roles in the process of gene transcription by binding to the promoter elements. In Bacillus, CodY and Spo0A have been recognized as global regulators. Cody, a GTP-binding protein, regulates more than one hundred genes that are typically repressed during rapid (exponential) growth and induced when cells experience nutrient deprivation . Our data indicated that expression of codY was higher at the exponential growth phase and then declined. In contrast, the expression of Spo0A, a master regulator of sporulation initiation and secondary metabolism in B. subtilis , continuously increased up to 48 h. A similar expression pattern for Spo0A was also observed in B. subtilis .
Validation of the selected DEGs by real-time PCR
To confirm the accuracy and reproducibility of the transcriptome data, 8 genes were selected for qPCR validation. RNA samples from the cultures of MM and GM at different growth phases were used as templates. The data are shown in Additional file 8, indicating that the two sets of data between the RNA-seq and RT-PCR analyses were almost consistent.
In this study, we determined the transcriptome profiles of B. pumilus across the various growth phases in both MM and GM cultures. Through DEG analysis and KEGG enrichment, transcriptional changes to genes that are specific to certain metabolic pathways were unveiled at the various growth phases.
During the exponential growth phase, cellular growth is most active along with a quick depletion of nutrients. When the culture enters the transition phase, stress may occur at least through nutrient depletion . At this transition point, a critical alteration in TCA metabolism was revealed in B. pumilus BA06. Our results showed that the TCA cycle was most active during the exponential growth phase (Fig. 3). Similarly, TCA -cycle enzymes were frequently identified from the B. pumilus cells during the exponential growth phase in the proteomic analysis . The TCA cycle is a central metabolic pathway that not only unifies the carbohydrate, fat and protein metabolic processes, but also produces energy and reducing power. Therefore, the slowdown in the TCA cycle at the transition point and the later growth stage may be a signature of reduced cell activity. Previous studies showed that the expression levels of some genes involved in the TCA cycle decreased when B. subtilis was grown under anaerobic conditions . In addition, glucose may regulate the TCA cycle by repressing several genes for the TCA . In fact, at the transition phase and stationary growth phase, glucose in the culture is consumed. Therefore, regulation of the TCA cycle may include other factors. However, secondary metabolism and sporulation will be triggered in the stationary growth phase, which requires the TCA metabolism to maintain carbohydrate flux and energy. Based on our transcriptome data, the carbohydrate flux could be recycled via another branch pathway, from acetyl-CoA to phosphoenol-pyruvate via oxaloacetate, since the genes encoding pyruvate carboxyl transferase (peg.413) and phosphoenolpyruvate carboxykinase (peg.2100) were not repressed but were greatly induced in expression upon entry into the transition phase and later growth phases (Fig. 3).
Another significant change at the transition point is down-regulation of the genes involved in flagellar assembly and the chemotaxis signaling pathway. The flagellum is an organ for mobility, which enables bacterial cells to colonize an ecological niche . Mobility is also regulated by chemotaxis factors . There are several chemoreceptors known as methyl-accepting chemotaxis proteins (MCPs) that are encoded by B. pumilus, most of which were revealed as being down-regulated after the exponential growth phase (Additional file 4). However, for other factors, such as surfactin, extracellular proteolytic activity was required to support mobility [51, 52]. The operon srf (peg.2905–2910) for biosurfactin biosynthesis was also down-regulated in B. pumilus BA06 at the same growth phase (Additional file 4). In fact, B. pumilus BA06 could produce surfactant at 24 h of culture, although the yield was lower . Together, these data indicate the attenuation of the swarming mobility of B. pumilus upon entering the transition point and thereafter.
Sporulation is a complex cellular process and can be triggered by nutrient limitations . For B. pumilus BA06, the formation of endospores was quick during the stationary growth phase (Fig. 1). Meanwhile, our transcriptome data indicated that almost all the genes involved in sporulation continued to be up-regulated after exponential growth up to the stationary growth phase (Fig. 5 and Additional file 5). The addition of gelatin to the MM led to a delay in both sporulation and expression of the sporulation-related genes. However, sporulation is subject to highly hierarchical regulation. In general, the transcriptional factor Spo0A governs the decision to initiate sporulation by phosphorylation via a series of kinases and phosphotransferases . The expression level of Spo0A (peg.1484) increased slowly from 4 h to 48 h, which was consistent with a previous study in B. subtilis . Since the phosphorylation state is a key factor for Spo0A to initiate sporulation , the transcription level observed here may not be relevant to the explanation of its function in sporulation initiation. Five kinases were identified in B. pumilus, of which the expression of KinA (peg.311, peg.312), KinC (peg.373), KinD (peg.276), and KinE (peg.265) was up-regulated. In contrast, KinB (peg.909) was down-regulated at the transition point (Additional file 5). In addition, the sporulation initiation phosphotransferase (Spo0F, peg.3334) was also down-regulated at 12 h and then up-regulated at 48 h and thereafter. Overall, the expression levels of the conserved sporulation kinases and phosphotransferases and Spo0A displayed a pattern that was different from the other sporulation-related genes, whose expression was continuously up-regulated until the later stationary growth phase (Fig. 5). The interruption may be ascribed to their role as only being responsive toward initiating sporulation. Recently, KinD in B. subtilis was shown to delay the onset of sporulation . However, the expression level of KinD in B. pumilus was highest among the five sporulation kinases, which may directly phosphorylate Spo0A . Therefore, KinD may play a more important role in sporulation initiation in B. pumilus.
During the sporulation process, the roles of several specific sigma factors SigE, SigF, SigG, and SigK were well documented in regulation of sporulation. These four factors were sequentially activated after sporulation initiation: sigF in the forespore, sigE in the mother cell, sigG in the forespore, and sigK in the mother cell . Our transcriptome data indicated that all four sigma factors were expressed in a similar way with a peak at the later stationary growth phase (Fig. 5). However, another sigma factor, SigH (peg.2382), was found to participate in the regulation of sporulation in Clostridium difficile . In B. pumilus, sigH was expressed in a manner that was different from the above four sigma factors. SigH may regulate the other genes involved in mobility and cell division. Therefore, the sigH gene may not be specific to sporulation in B. pumilus.
Extracellular protease production and protein secretion are processes of interest to the field of biotechnology. Previously, AprE, Epr, WprA, and Vpr were identified to be extracellular protease in B. pumilus SCU11 . Our transcriptome data indicated that the expression patterns of various extracellular proteases were different in B. pumilus (Table 3). Three genes (aprE, aprX, and wprA) were expressed more highly during the stationary growth phase. Based on the expression levels, these proteases may contribute to the high extracellular proteolytic activity observed in the culture supernatant. In contrast, epr and subE were highly expressed in the exponential growth phase (Table 3), which could be explained by the fact that the expression of epr was controlled by SigD in B. subtilis . However, these results suggested that different regulatory mechanisms and secretion pathways may be employed for the various proteases. In B. subtilis, the regulatory mechanism of the aprE expression has been extensively studied, indicating that at least two positive regulatory proteins (DegU and Spo0A) are involved in regulation of the aprE expression [59, 60]. The expression of degU and spo0A in B. pumilus was up-regulated at 48 h in GM (Additional file 1 and Fig. 7), which may contribute to the observed increase of aprE expression in the medium supplemented with gelatin.
Two protein secretion pathways were enriched by KEGG analysis in B. pumilus: The Sec-dependent pathway and Tat system, which were also identified by the previous proteomics analysis . The Tat system responds by to secreting twin-arginine (RR/KRP) signal peptides, through which 44 proteins are predicated to be secreted in B. subtilis . WprA has been proposed to be secreted through the Tat system in B. pumilus . Regarding the expression pattern during the growth course, the Tat system (including TatA and TatC) may respond by secreting limited proteins and may play a major role in the transition phase. In contrast, the large number of genes involved in the Sec secretion pathway are down-regulated at the transition point, and then up-regulated at 48 h. In general, many proteins have been observed to be secreted extracellularly during the stationary growth phase in B. subtilis . Therefore, the Sec secretion system may respond by secreting many proteins during the stationary growth phase. For example, AprE, Vpr, and aprX may be secreted via this system . However, 65 proteins were detected from the late exponential growth cultures in B. pumilus; this number is much lower than the predicted protein species (513) . Therefore, the Sec-dependent protein secretion system may be the major pathway, especially during the stationary growth phase in B. pumilus.
In Bacillus, alternative sigma factors are involved in the regulation of certain genes or specific metabolic processes. For example, SigB is generally recognized in response to stress. In B. subtilis, approximately 150 general stress-associated genes have been identified as the SigB regulon . The expression of the sigB gene (peg.3035) is also up-regulated greatly at the transition phase in B. pumilus, which is similar with the observation in B. subtilis . Although the sigB regulon has not been identified for B. pumilus, the rsb operon (peg.3029–3033) is induced upon entry into the transition point (Additional file 9). Other sigma factors such as SigW and SigX have been reported to respond to stress. For example, SigW in B. subtilis responds by regulating detoxification and the production of antimicrobial compounds . However, the genes involved in the osmotic and oxidative responses, such as opuAA (peg. 2747, peg.3708), opuAB (peg.3707), opuCB (peg.1139), opuCA (peg.1140), perR (peg.820), and cat (peg.2205, peg.3703), have been shown to be up-regulated at the transition point in B. pumilus (Additional file 8: Table S8); some of these genes have also been revealed to be involved in the response to oxidative stress in another study using B. pumilus Jo2 and a microarray-based transcriptome analysis . Therefore, B. pumilus may suffer from stress when entering the transition phase and subsequently activates a set of stress-related genes to be expressed in a manner that is similar to B. subtilis.
In conclusion, an RNA-seq-based transcriptome analysis was first applied to B. pumilus BA06 to monitor the transcriptional profile over the entire growth course in a defined MM or GM medium. Changes to the specific metabolic pathways in terms of gene expression were shown to relate to the transition from the exponential growth phase and the onset of the stationary growth phase. Upon entry into the transition point, one of the striking changes was down-regulation of central TCA metabolism and oxidative phosphorylation. Furthermore, the cellular mobility ability was also reduced because of down-regulation of genes involved in flagellar assembly and the chemotaxis signaling pathway. In contrast, many stress-responding genes including the SigB regulon, were induced in expression after entering the transition point, indicating that the cells may suffer from stress. During the stationary growth phase, a significant change occurred in the genes involved in sporulation, and more than 100 sporulation-related genes were induced, which was consistent with the process of sporulation formation. However, the addition of gelatin to the MM medium did not cause a great impact on the transcriptome profile, except for the genes encoding sporulation-related proteins and extracellular proteases with delayed expression.
Bacterial strain and growth conditions
B. pumilus BA06 was routinely maintained on Luria-Bertani (LB, 10 g/l tryptone, 5 g/l yeast extract, 10 g/l NaCl, pH 7.5) agar plate. A single colony of B. pumilus BA06 was selected for transfer into 10 ml LB broth and incubated at 37 °C overnight with shaking at 140 rpm. Afterward, 500 μl of an overnight culture was transferred into 50 ml of minimal medium [MM, 1.0 g/l sodium citrate, 2.0 g/l (NH4)2SO4, 14.0 g/l K2HPO4, 6.0 g/l KH2PO4, 0.2 g/l MgSO4, 2.5 g/l yeast extract, 5.0 g/l D-glucose] and the gelatin-amended minimal medium (GM) (2.0 g/l gelatin) in 250-ml flasks. The cultures were incubated at 37 °C with shaking at 140 rpm for the indicated time points.
At various time points, the cell density was measured by reading the OD560 on a spectrometer. The total number of cells and endospores was also determined. To account for the total cell number, a 10-fold dilution of the fresh culture was achieved by serial dilution in sterile PBS buffer. Finally, 0.1 ml of the cell suspension was dispensed onto LB agar plates, and the colonies assigned as the total cells were counted after incubation at 37 °C. To count the number of endospores, an aliquot of fresh culture was sampled and incubated at 55 °C for 15 min to kill the vegetative cells; and the samples were diluted as above and dispensed onto LB agarose plates. The colonies formed on these plates were regarded as endospores. Meanwhile, extracellular protease activity was also assayed using casein as the substrate as described previously . All experiments were performed in triplicate.
RNA isolation, library construction and Illumina sequencing
The cell samples were pelleted by centrifugation at 8000 rpm at the indicated time points (4, 12, 24, 48, and 72 h) from the B. pumilus BA06 cultures in MM and GM, and the cells were then suspended in TE buffer supplemented with 1.5 mg/ml lysozyme and incubated at 37 °C for 10 min, to which the TRIzol reagent (Invitrogen, Invitrogen, Carlsbad, CA) was added. The cell suspension was mixed extensively while using the gauge to disrupt the cells completely. Finally, total RNA was isolated following the instructions provided with the TRIzol reagent. The genomic DNA was removed using the Genome DNA Eraser kit (Takara, Dalian, China).
Subsequently, the rRNA was removed using the Ribo-ZeroTM rRNA Removal kit (Epicentre Biotech, Madison, WI). The resulting mRNAs were fragmented and reverse transcribed using random hexamers as the primer. Second strand cDNA synthesis was performed using DNA Polymerase I and RNase H. The cDNA fragments were processed for end repair and ligated to paired end adaptors. Finally, the library was constructed and sequenced on an Illumina HiSeqTM2000 sequencing platform.
Mapping and identification of the differentially expressed genes
Clean data were obtained from the raw data by removing the sequences of the adapters and low-quality reads. The clean reads were aligned to the B. pumilus genome (GenBank accession number: AMDH00000000)  using Bowtie2 with default parameters, allowing up to one-base mismatches . The aligned read files were processed by Cufflinks v2.2.1 . The relative abundances of the transcripts were calculated as the fragments per kilobase of transcript per million fragments mapped (FPKM). Differentially expressed genes (DEGs) among the different samples were extracted by using edgeR in the Bioconductor package . The DEGs were defined with an FDR (false discovery rate) ≤ 0.05 and log2fold-change (log2FC) ≥ 1.
Functional classification and KEGG analysis
Functional classification of DEGs was performed online by RAST (Rapid Annotation using Subsystem Technology, http://www.nmpdr.org/FIG/wiki/view.cgi/FIG/ RapidAnnotationServer). RAST is a fully automated service that is especially designed for annotating bacterial and archaeal genomes . Once one gene is annotated, it can be classified into various subsystems.
Meanwhile, the DEGs were further assigned to KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways on the KEGG Automatic Annotation Server (http://www.genome.jp/kegg).
Quantitative real-time PCR analysis
The expression levels of the selected eight genes (degS, aprX, glnR, hpr, vpr, sinR, yqkD, and spo0A) of B. pumilus BA06 growing under the same conditions, were validated by real-time RT-PCR analysis. The real-time RT-PCR was performed using an iCyclerMyiQ Real-Time PCR System (Bio-Rad, Hercules, CA). The PCR conditions were set up as 95 °C for 2 mins, followed by 40 cycles of 95 °C for 10 s, 65 °C for 15 s and 72 °C for 20 s. A melting curve analysis of the amplification products was performed at the end of each PCR run to ensure that unique products were amplified. The specific primers used for the selected genes are listed in Additional file 10. The expression level was normalized to the internal control gene 16S rRNA, using the 2−ΔΔCt method .
We thank the American Journal Experts (www.aje.com) for language edition for our manuscript.
This study was financially supported by the National Natural Science Fund of China (31171204). The funding sponsor had no role in the design of the study, collections, analysis, interpretation of data, writing of the manuscript, or decision to publish the results.
Availability of data and materials
The raw Illumina sequencing dataset of the B. pumilus BA06 strain is available on the NCBI Sequence Read Archive (ARS) database under accession SRR5115698-707. The data supporting the conclusions of this article are included within the article and its additional files.
HF and HYW conceived and designed the experiments; LLH, HHS, GL, CYX and XJC performed the experiments; LLH, HHS, YCL and HF analyzed the data; HF and TXM contributed reagents and materials; HF, LLH and HHS wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Liu Y, Lai Q, Dong C, Sun F, Wang L, Li G, et al. Phylogenetic diversity of the Bacillus pumilus group and the marine ecotype revealed by multilocus sequence analysis. PLoS One. 2013;8:e80097.View ArticlePubMedPubMed CentralGoogle Scholar
- Connor N, Sikorski J, Rooney AP, Kopac S, Koeppel AF, Burger A, et al. Ecology of speciation in the genus Bacillus. Appl Environ Microbiol. 2010;76:1349–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim HK, Choi HJ, Kim MH, Sohn CB, Oh TK. Expression and characterization of Ca2+-independent lipase from Bacillus pumilus B26. Biochim Biophys Acta. 2002;1583:205–12.View ArticlePubMedGoogle Scholar
- Battan B, Sharma J, Dhiman SS, Kuhad RC. Enhanced production of cellulase-free thermostable xylanase by Bacillus pumilus ASH and its potential application in paper industry. Enzy Microbial Technol. 2007;41:733–9.View ArticleGoogle Scholar
- Huang Q, Peng Y, Li X, Wang H, Zhang Y. Purification and characterization of an extracellular alkaline serine protease with dehairing function from Bacillus pumilus. Curr Microbiol. 2003;46:169–73.View ArticlePubMedGoogle Scholar
- Kumar AG, Swarnalatha S, Gayathri S, Nagesh N, Sekaran G. Characterization of an alkaline active-thiol forming extracellular serine keratinase by the newly isolated Bacillus pumilus. J Appl Microbiol. 2011;104:411–9.Google Scholar
- Jaouadi B, Ellouz-Chaaboum S, Rhimi M, Bejar S. Biochemical and molecular characterization of a detergent-stable serine alkaline protease from Bacillus pumilus CBS with high catalytic efficiency. Biochimie. 2008;90:1291–305.View ArticlePubMedGoogle Scholar
- Küppers T, Steffen V, Hellmuth H, O’Connell T, Bongaerts J, Maurer KH, et al. Developing a new production host from a blueprint: Bacillus pumilus as an industrial enzyme producer. Microbial Cell Fact. 2014;13:46.View ArticleGoogle Scholar
- Wemhoff S, Meinhardt F. Generation of biologically contained, readily transformable, and genetically manageable mutants of the biotechnologically important Bacillus pumilus. Appl Microbiol Biotechnol. 2013;97:7805–19.View ArticlePubMedGoogle Scholar
- Hua D, Ma C, Lin S, Song L, Deng Z, Maomy Z, et al. Biotransformation of isoeugenol to vanillin by a newly isolated Bacillus pumilus strain: identification of major metabolites. J Biotechnol. 2007;130:463–70.View ArticlePubMedGoogle Scholar
- Srivastava RK, Jaiswal R, Panda D, Wangikar PP. Megacell phenotype and its relation to metabolic alterations in transketolase deficient strain of Bacillus pumilus. Biotechnol Bioeng. 2009;102:1387–97.View ArticlePubMedGoogle Scholar
- Bacillus pumilus strain GB 34 (006493) Fact Sheet. https://www3.epa.gov/pesticides/chem_search/reg_actions/registration/fs_PC-006493_13-Mar-03.pdf.
- Zhao CW, Wang HY, Zhang YZ, Feng H. Draft genome sequence of Bacillus pumilus BA06, a producer of alkaline serine protease with leather-dehairing function. J Bacteriol. 2012;194:6668–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Wiegand S, Dietrich S, Hertel R, Bongaerts J, Evers S, Volland S, et al. RNA-Seq of Bacillus licheniformis: active regulatory RNA features expressed within a productive fermentation. BMC Genomics. 2013;14:1.View ArticleGoogle Scholar
- Handtke S, Schroeter R, Jürgen B, Methling K, Schlüter R, Albercht D, et al. Bacillus pumilus reveals a remarkably high resistance to hydrogen peroxide provoked oxidative stress. PLoS One. 2014;9:e85625.View ArticlePubMedPubMed CentralGoogle Scholar
- Handtke S, Volland S, Methling K, Albrecht D, Becher D, Nehls J, et al. Cell physiology of the biotechnical relevant bacterium Bacillus pumilus—an omics-based approach. J Biotechnol. 2014;192:204–14.View ArticlePubMedGoogle Scholar
- Sorek R, Cossart P. Prokaryotic transcriptomics: a new view on regulation, physiology and pathogenicity. Nat Rev Genet. 2010;11:9–16.View ArticlePubMedGoogle Scholar
- Nagalakshmi U, Waern K, Snyder M. RNA-Seq: a method for comprehensive transcriptome analysis. Curr Protocol Mol Biol. 2010;4:4.11.11–14.11.13.Google Scholar
- Nicolas P, Mäder U, Dervyn E, Rochat T, Leduc A, Pigeonneau N, et al. Condition-dependent transcriptome reveals high-level regulatory architecture in Bacillus subtilis. Science. 2012;335:1103–6.View ArticlePubMedGoogle Scholar
- Brinsmade SR, Alexander EL, Livny J, Stettner AI, Segrè D, Rhee KY, et al. Hierarchical expression of genes controlled by the Bacillus subtilis global regulatory protein CodY. Proc Natl Acad Sci U S A. 2014;111:8227–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Guo J, Cheng G, Gou XY, Xing F, Li S, Han YC, et al. Comprehensive transcriptome and improved genome annotation of Bacillus licheniformis WX-02. FEBS Lett. 2015;589:2372–81.View ArticlePubMedGoogle Scholar
- Bassi D, Colla F, Gazzola S, Puglisi E, Delledonne M, Cocconcelli PS. Transcriptome analysis of Bacillus thuringiensis spore life, germination and cell outgrowth in a vegetable-based food model. Food Microbiol. 2016;55:73–85.View ArticlePubMedGoogle Scholar
- Wang HY, Liu DM, Liu Y, Cheng CF, Ma QY, Huang Q, et al. Screening and mutagenesis of a novel Bacillus pumilus strain producing alkaline protease for dehairing. Lett Appl Microbiol. 2006;44:1–6.View ArticleGoogle Scholar
- Wan MY, Wang HY, Zhang YZ, Feng H. Substrate specificity and thermostability of the dehairing alkaline protease from Bacillus pumilus. Appl Biochem Biotechnol. 2009;159:394–403.View ArticlePubMedGoogle Scholar
- Liu RF, Huang CL, Feng H. Salt stress represses production of extracellular proteases in Bacillus pumilus. Genet Mol Res. 2015;14:4339–948.Google Scholar
- Rao MB, Tanksale AM, Ghatge MS, Deshpande VV. Molecular and biotechnological aspects of microbial proteases. Microbiol Mol Biol Rev. 1998;62:597–635.PubMedPubMed CentralGoogle Scholar
- Gupta R, Beg QK, Lorenz P. Bacterial alkaline proteases: molecular approaches and industrial applications. Appl Microbiol Biotechnol. 2002;59:15–32.View ArticlePubMedGoogle Scholar
- Fakhfakh-Zouari N, Haddar A, Hmidet N, Frikha F, Nasri M. Application of statistical experimental design for optimization of keratinases production by Bacillus pumilus A1 grown on chicken feather and some biochemical properties. Process Biochem. 2010;45:617–26.View ArticleGoogle Scholar
- Schallmey M, Singh A, Ward OP. Developments in the use of Bacillus species for industrial production. Can J Microbiol. 2003;50:1–17.View ArticleGoogle Scholar
- Manabe K, Kageyama Y, Morimoto T, Shimizu E, Takahashi H, Kanaya S, et al. Improved production of secreted heterologous enzyme in Bacillus subtilis strain MGB874 via modification of glutamate metabolism and growth conditions. Microbial Cell Fact. 2013;12:18.View ArticleGoogle Scholar
- Setlow P. Spores of Bacillus subtilis: their resistance to and killing by radiation, heat and chemicals. J Appl Microbiol. 2006;101:514–25.View ArticlePubMedGoogle Scholar
- Strauch MA, Hoch JA. Transition-state regulators: sentinels of Bacillus subtilis post-exponential gene expression. Mol Microbiol. 1993;7:337–42.View ArticlePubMedGoogle Scholar
- Brigulla M, Hoffmann T, Krisp A, Völker A, Bremer E, Völker U. Chill induction of the SigB-dependent general stress response in Bacillus subtilis and its contribution to low-temperature adaptation. J Bacteriol. 2003;185:4305–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Akbar S, Gaidenko TA, Kang CM, O’Reilly M, Devine KM, Price CW. New family of regulators in the environmental signaling pathway which activates the general stress transcription factor sigma B of Bacillus subtilis. J Bacteriol. 2001;183:1329–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75.View ArticlePubMedPubMed CentralGoogle Scholar
- Blom E-J, Ridder AN, Lulko AT, Roerdink JB, Kuipers OP. Time-resolved transcriptomics and bioinformatic analyses reveal intrinsic stress responses during batch culture of Bacillus subtilis. PLoS One. 2011;6:e27160.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang CK, Tai PC, Lu CD. Time-related transcriptome analysis of B. subtilis 168 during growth with glucose. Curr Microbiol. 2014;68:12–20.View ArticlePubMedGoogle Scholar
- Zamboni N, Maaheimo H, Szyperski T, Hohmann HP, Sauer U. The phosphoenolpyruvate carboxykinase also catalyzes C3 carboxylation at the interface of glycolysis and the TCA cycle of Bacillus subtilis. Metabolic Eng. 2004;6:277–84.View ArticleGoogle Scholar
- Mukherjee S, Kearns DB. The structure and regulation of flagella in Bacillus subtilis. Ann Rev Genet. 2014;48:319–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Kearns DB, Losick R. Swarming mobility in undomesticated Bacillus subtilis. Mol Microbiol. 2003;49:581–90.View ArticlePubMedGoogle Scholar
- Higgins D, Dworkin J. Recent progress in Bacillus subtilis sporulation. FEMS Microbiol Rev. 2011;36:131–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu L, Liu Y, Shin H, Chen RR, Wang NS, Li J, et al. Developing Bacillus spp. as a cell factory for production of microbial enzymes and industrially important biochemicals in the context of systems and synthetic biology. Appl Microbiol Biotechnol. 2013;97:6113–27.View ArticlePubMedGoogle Scholar
- Voelker U, Voelker A, Maul B, Hecker M, Dufour A, Haldenwang WG. Separate mechanisms activate sigma B of Bacillus subtilis in response to environmental and metabolic stresses. J Bacteriol. 1995;177:3771–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang S, Haldenwang WG. Contributions of ATP, GTP, and redox state to nutritional stress activation of the Bacillus subtilis Sigma B transcription factor. J Bacteriol. 2005;187:7554–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Sonenshein AL. CodY, a global regulator of stationary phase and virulence in gram-positive bacteria. Curr Opin Microbiol. 2005;8:203–7.View ArticlePubMedGoogle Scholar
- Molle V, Fujita M, Jensen ST, Eichenberger P, González-Pastor JE, Liu JS, et al. The Spo0A regulon of Bacillus subtilis. Mol Microbiol. 2003;50:1683–701.View ArticlePubMedGoogle Scholar
- Sung H-M, Yasbin RE. Adaptive, or stationary-phase, mutagenesis, a component of bacterial differentiation in Bacillus subtilis. J Bacteriol. 2002;184:5641–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakano MM, Zuber P, Sonenshein AL. Anaerobic regulation of Bacillus subtilis Krebs cycle genes. J Bacteriol. 1998;180:3304–11.PubMedPubMed CentralGoogle Scholar
- Sonenshein AL. Control of key metabolic intersections in Bacillus subtilis. Nat Rev Microbiol. 2007;5:917–27.View ArticlePubMedGoogle Scholar
- Harshey RM, Matsuyama T. Dimorphic transition in Escherichia coli and Samonella typhimurium: surface-induced differentiation into hyperflagte swarmer cells. Proc Natl Acad Sci U S A. 1994;91:8631–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Connelly MB, Young GM, Sloma A. Extracellular proteolytic activity plays a central role in swarming mobility in Bacillus subtilis. J Bacteriol. 2004;186:4159–67.View ArticlePubMedGoogle Scholar
- Kerrns DB, Losick R. Swarming mobility in undomesticated Bacillus subtilis. Mol Microbiol. 2003;49:581–90.View ArticleGoogle Scholar
- Ji Y, Feng H. Optimal of fermentation medium for producing surfactin by Bacillus pumilus BA06. J Sichuan University (Nat Sci ed). 2016;53:925–30.Google Scholar
- Fimlaid KA, Shen A. Diverse mechanisms regulate sporulation sigma factor activity in the Firmicutes. Curr Opin Microbiol. 2015;24:88–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Aguilar C, Vlanakis H, Guzman A, Losick R, Kolter R. KinD is a checkpoint protein linking spore formation to extracellular-matrix production in Bacillus subtilis biofilms. mBio. 2010;1:e00035–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Saujet L, Monot M, Dupuy B, Soutourina O, Martin-Verstraete I. The key sigma factor of transition phase, sigH, control sporulation, metab0lism, and virulence factor expression in Clostridium difficile. J Bacteriol. 2011;193:3186–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang C, Yu S, Song T, He T, Shao H, Wang H. Extracellular proteome profiling of Bacillus pumilus SCU11 producing alkaline protease for dehairing. J Microbiol Biotechnol. 2016;26:1993–2005.View ArticlePubMedGoogle Scholar
- Dixit M, Murudkar CS, Rao KK. epr is transcribed from σD promoter and is involved in swarming of Bacillus subtilis. J Bacteriol. 2002;184:596–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Yasumura A, Abe S, Tanaka T. Involvement of nitrogen regulation in Bacillus subtilis degU expression. J Bacteriol. 2008;190:5162–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Kodama T, Endo K, Ara K, Ozaki K, Kakeshita H, Yamane K, et al. Effect of Bacillus Subtilis spo0A mutation on cell wall lytic enzymes and extracellular proteases, and prevention of cell lysis. J Biosci Bioeng. 2007;103:13–21.View ArticlePubMedGoogle Scholar
- Jongbloed JDH, Antelmann H, Hecker M, Nijland R, Pries F, Koski P, et al. Selective contribution of the twin-arginine translocation pathway to protein secretion in Bacillus subtilis. J Biol Chem. 2002;277:44068–78.View ArticlePubMedGoogle Scholar
- Tjalsma H, Antlmann H, Jongbloed JDH, Braun PG, Darmon E, Dorenbos R, et al. Proteomics of protein secretion by Bacillus subtilis: separating the “secrets” of the secrtome. Microbiol Mol Biol Rev. 2004;68:207–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Antelmann H, Tjalsma H, Voigt B, Ohlmeier S, Bron S, van Dijl JM, et al. A proteomic view on genome-based signal peptide predictions. Genome Res. 2001;11:14984–502.View ArticleGoogle Scholar
- Nannapaneni P, Hertwig F, Depke M, Hecker M, Mäder U, Völker U, et al. Defining the structure of the general stress regulon of Bacillus subtilis using targeted microarray analysis and random forest classification. Microbiol. 2012;158:696–707.View ArticleGoogle Scholar
- Turner MS, Helmann JD. Mutations in multidrug efflux homologs, sugar isomerases, and antimicrobial biosynthesis genes differentially elevate activity of the sigma X and sigma W factors in Bacillus subtilis. J Bacteriol. 2000;182:5202–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Method. 2012;9:357–9.View ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protocol. 2012;7:562–78.View ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Method. 2001;25:402–8.View ArticleGoogle Scholar