Bacillus subtilis DB104, an extracellular protease-deficient derivative of B. subtilis 168, is widely used for recombinant protein expression. An understanding of the changes in gene expression during growth is essential for the commercial use of bacterial strains. Transcriptome and proteome analyses are ideal methods to study the genomic response of microorganisms. In this study, transcriptome analysis was performed to monitor changes in the gene expression level of B. subtilis DB104 while growing on a complete medium. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, K-mean cluster analysis, gene ontology (GO) enrichment analysis, and the function of sigma factors were used to divide 2122 differentially expressed genes (DEGs) into 10 clusters and identified gene functions according to expression patterns. The results of KEGG pathway analysis indicated that ABC transporter is down-regulated during exponential growth and metabolic changes occur at the transition point where sporulation starts. At this point, several stress response genes were also turned on. The genes involved in the lipid catabolic process were up-regulated briefly at 15 h as an outcome of the programmed cell death that postpones sporulation. The results suggest that changes in the gene expression of B. subtilis DB104 were dependent on the initiation of sporulation. However, the expression timing of the spore coat gene was only affected by the relevant sigma factor. This study can help to understand gene expression and regulatory mechanisms in B. subtilis species by providing an overall view of transcriptional changes during the growth of B. subtilis DB104.
Keywords: Bacillus subtilis DB104; RNA-seq; gene expression; gene regulation; transcriptome analysis.