Genomic insights into the taxonomic status of the three subspecies of Bacillus subtilis

Isolation and Characterization of the Zearalenone-Degrading Strain, Bacillus spizizenii B73, Inspired by Esterase Activity

Zearalenone (ZEN) is a widespread mycotoxin found in grain and feed, presenting a serious threat to animal and human health. This study investigated the ability of the novel strain B73, isolated from petroleum-contaminated soil, to detoxify ZEN. B73 was identified as Bacillus spizizenii through physiological and biochemical tests, and further confirmed based on the 16S rRNA gene sequence and the complete genome sequence. B. spizizenii B73 was capable of degrading up to 99.3% of ZEN at a concentration of 10 μg/mL in a minimal medium (pH = 7.0) within 8 h at 37 °C via HPLC-UV. In addition, B. spizizenii B73 was used to treat ZEN-contaminated wheat bran, dried distillers grains (DDGS), and corn meal, whereby the respective degradation rates reached 96.32%, 98.73%, and 80.31% after 36 h of treatment. HPLC-Q-Exactive-MS/MS analysis revealed one of the degradation products to have the formula C17H24O4. B. spizizenii B73 is a novel strain isolated from petroleum-contaminated soil, and the extracellular enzymes secreted by this strain show a remarkable ability to degrade ZEN.

Broad-Spectrum Antifungal, Biosurfactants and Bioemulsifier Activity of Bacillus subtilis subsp. spizizenii-A Potential Biocontrol and Bioremediation Agent in Agriculture

In this study, the antifungal, biosurfactant and bioemulsifying activity of the lipopeptides produced by the marine bacterium Bacillus subtilis subsp. spizizenii MC6B-22 is presented. The kinetics showed that at 84 h, the highest yield of lipopeptides (556 mg/mL) with antifungal, biosurfactant, bioemulsifying and hemolytic activity was detected, finding a relationship with the sporulation of the bacteria. Based on the hemolytic activity, bio-guided purification methods were used to obtain the lipopeptide. By TLC, HPLC and MALDI-TOF, the mycosubtilin was identified as the main lipopeptide, and it was further confirmed by NRPS gene clusters prediction based on the strain's genome sequence, in addition to other genes related to antimicrobial activity. The lipopeptide showed a broad-spectrum activity against ten phytopathogens of tropical crops at a minimum inhibitory concentration of 400 to 25 μg/mL and with a fungicidal mode of action. In addition, it exhibited that biosurfactant and bioemulsifying activities remain stable over a wide range of salinity and pH and it can emulsify different hydrophobic substrates. These results demonstrate the potential of the MC6B-22 strain as a biocontrol agent for agriculture and its application in bioremediation and other biotechnological fields.

Perception of Biocontrol Potential of Bacillus inaquosorum KR2-7 against Tomato Fusarium Wilt through Merging Genome Mining with Chemical Analysis

Simple Summary

Bacillus is a bacterial genus that is widely used as a promising alternative to chemical pesticides due to its protective activity toward economically important plant pathogens. Fusarium wilt of tomato is a serious fungal disease limiting tomato production worldwide. Recently, the newly isolated B. inaquosorum strain KR2-7 considerably suppressed Fusarium wilt of tomato plants. The present study was performed to perceive potential direct and indirect biocontrol mechanisms implemented by KR2-7 against this disease through genome and chemical analysis. The potential direct biocontrol mechanisms of KR2-7 were determined through the identification of genes involved in the synthesis of antibiotically active compounds suppressing tomato Fusarium wilt. Furthermore, the indirect mechanisms of this bacterium were perceived through recognizing genes that contributed to the resource acquisition or modulation of plant hormone levels. This is the first study that aimed at the modes of actions of B. inaquosorum against Fusarium wilt of tomatoes and the results strongly indicate that strain KR2-7 could be a good candidate for microbial biopesticide formulations to be used for biological control of plant diseases and plant growth promotion.

Abstract

Tomato Fusarium wilt, caused by Fusarium oxysporum f. sp. lycopersici (Fol), is a destructive disease that threatens the agricultural production of tomatoes. In the present study, the biocontrol potential of strain KR2-7 against Fol was investigated through integrated genome mining and chemical analysis. Strain KR2-7 was identified as B. inaquosorum based on phylogenetic analysis. Through the genome mining of strain KR2-7, we identified nine antifungal and antibacterial compound biosynthetic gene clusters (BGCs) including fengycin, surfactin and Bacillomycin F, bacillaene, macrolactin, sporulation killing factor (skf), subtilosin A, bacilysin, and bacillibactin. The corresponding compounds were confirmed through MALDI-TOF-MS chemical analysis. The gene/gene clusters involved in plant colonization, plant growth promotion, and induced systemic resistance were also identified in the KR2-7 genome, and their related secondary metabolites were detected. In light of these results, the biocontrol potential of strain KR2-7 against tomato Fusarium wilt was identified. This study highlights the potential to use strain KR2-7 as a plant-growth promotion agent.

Bacillus velezensis: A Treasure House of Bioactive Compounds of Medicinal, Biocontrol and Environmental Importance

Bacillus velezensis gram-positive bacterium, is frequently isolated from diverse niches mainly soil, water, plant roots, and fermented foods. B. velezensis is ubiquitous, non-pathogenic and endospore forming. Being frequently isolated from diverse plant holobionts it is considered host adapted microorganism and recognized of high economic importance given its ability to promote plant growth under diverse biotic and abiotic stress conditions. Additionally, the species suppress many plant diseases, including bacterial, oomycete, and fungal diseases. It is also able after plant host root colonization to induce unique physiological situation of host plant called primed state. Primed host plants are able to respond more rapidly and/or effectively to biotic or abiotic stress. Moreover, B. velezenis have the ability to resist diverse environmental stresses and help host plants to cope with, including metal and xenobiotic stresses. Within species B. velezensis strains have unique abilities allowing them to adopt different life styles. Strain level abilities knowledge is warranted and could be inferred using the ever-expanding new genomes list available in genomes databases. Pangenome analysis and subsequent identification of core, accessory and unique genomes is actually of paramount importance to decipher species full metabolic capacities and fitness across diverse environmental conditions shaping its life style. Despite the crucial importance of the pan genome, its assessment among large number of strains remains sparse and systematic studies still needed. Extensive knowledge of the pan genome is needed to translate genome sequencing efforts into developing more efficient biocontrol agents and bio-fertilizers. In this study, a genome survey of B. velezensis allowed us to (a) highlight B. velezensis species boundaries and show that Bacillus suffers taxonomic imprecision that blurs the debate over species pangenome; (b) identify drivers of their successful acquisition of specific life styles and colonization of new niches; (c) describe strategies they use to promote plant growth and development; (d) reveal the unlocked strain specific orphan secondary metabolite gene clusters (biosynthetic clusters with corresponding metabolites unknown) that product identification is still awaiting to amend our knowledge of their putative role in suppression of pathogens and plant growth promotion, and (e) to describe a dynamic pangenome with a secondary metabolite rich accessory genome.

6S-2 RNA deletion in the undomesticated B. subtilis strain NCIB 3610 causes a biofilm derepression phenotype

Bacterial 6S RNA regulates transcription via binding to the active site of RNA polymerase holoenzymes. 6S RNA has been identified in the majority of bacteria, in most cases encoded by a single gene. Firmicutes including Bacillus subtilis encode two 6S RNA paralogs, 6S-1 and 6S-2 RNA. Hypothesizing that the regulatory role of 6S RNAs may be particularly important under natural, constantly changing environmental conditions, we constructed 6S RNA deletion mutants of the undomesticated B. subtilis wild-type strain NCIB 3610. We observed a strong phenotype for the ∆6S-2 RNA strain that showed increased biofilm formation on solid media and the ability to form surface-attached biofilms in liquid culture. This phenotype remained undetected in derived laboratory strains (168, PY79) that are defective in biofilm formation. Quantitative RT-PCR data revealed transcriptional upregulation of biofilm marker genes such as tasA, epsA and bslA in the ∆6S-2 RNA strain, particularly during transition from exponential to stationary growth phase. Salt stress, which blocks sporulation at a very early stage, was found to override the derepressed biofilm phenotype of the ∆6S-2 RNA strain. Furthermore, the ∆6S-2 RNA strain showed retarded swarming activity and earlier spore formation. Finally, the ∆6S-1&2 RNA double deletion strain showed a prolonged lag phase of growth under oxidative, high salt and alkaline stress conditions, suggesting that the interplay of both 6S RNAs in B. subtilis optimizes and fine-tunes transcriptomic adaptations, thereby contributing to the fitness of B. subtilis under the unsteady and temporarily harsh conditions encountered in natural habitats.

Novel xylanase producing Bacillus strain X2: molecular phylogenetic analysis and its application for production of xylooligosaccharides

A Bacillus strain X2 that produced extracellular endo-xylanase (GH 11) (EC: 3.2.1.8) was isolated from the soil of the Northeast India region. This aerobic culture was Gram positive and endospore forming. Chemotaxonomic characterization showed variance with the fatty acid profile of related species in the Bacillus subtilis group. In Bacillus strain X2, distinct occurrence of iso-C14:0 lipids is absent in other related species. The 16S rRNA gene sequence homology showed 99% similarity with Bacillus subtilis subsp. inaquosorum. The phylogenetic analysis by the multilocus sequence analysis (MLSA) of the nucleotide sequence of six concatenated genes (16S rRNA, groEL, gyrA, polC, purH and rpoB) resolved the taxonomic position of the Bacillus strain X2 in the Bacillus subtilis subsp. group. The MLSA showed that it is a member of a clade that includes Bacillus subtilis subsp. stercoris. In in silico DNA-DNA hybridization (DDH), the highest matching score was obtained with Bacillus subtilis subsp. stercoris (87%). The in silico DDH of the genome (G + C 43.7 mol %) shared 48.5%, with Bacillus subtilis subsp. inaquosorum. The MLSA phylogenetic tree and the highest degree of DNA hybridization, indicating that it belongs to the Bacillus subtilis subspecies stercoris.

Genetic marker gene, recQ, differentiating Bacillus subtilis and the closely related Bacillus species

RecQ, which encodes a DNA helicase, was selected in searching for a marker gene of Bacillus subtilis and related species via genome mining. RecQ gene sequence similarity of type strains among Bacillus species used in this study ranged from 66.2% to 96.6%, whereas orthologous average nucleotide identity ranged from 72.6% to 95.8%. According to the phylogenetic tree based on recQ sequences, each type strain of all Bacillus species or subspecies used in this study was placed in a unique taxonomic position. Four B. subtilis subspecies, Bacillus tequilensis and Bacillus vallismortis were grouped in one cluster (cluster A). Strains of B. subtilis subsp. subtilis were classified into A1 cluster, and divided into subgroups. Isolates from Natto, Japanese fermented bean food, were classified into one subgroup, whereas those from Cheonggukjang, Korean fermented bean food, were divided into several subgroups within A1. Type strains of Bacillus halotolerans and Bacillus mojavensis were grouped into another cluster (cluster B), related to cluster A. Bacillus siamensis, Bacillus velezensis and Bacillus amyloliquefaciens were grouped into an independent cluster (cluster E). Sequencing of recQ was useful for the classification or differentiation of B. subtilis and closely related species. Therefore, recQ gene can be applied to the classification of these taxa.

Promotion of Bacillus subtilis subsp. inaquosorum, Bacillus subtilis subsp. spizizenii and Bacillus subtilis subsp. stercoris to species status

Bacillus subtilis currently encompasses four subspecies, Bacillus subtilis subsp. subtilis, Bacillus subtilis subsp. inaquosorum, Bacillus subtilis subsp. spizizenii and Bacillus subtilis subsp. stercoris. Several studies based on genomic comparisons have suggested these subspecies should be promoted to species status. Previously, one of the main reasons for leaving them as subspecies was the lack of distinguishing phenotypes. In this study, we used comparative genomics to determine the genes unique to each subspecies and used these to lead us to the unique phenotypes. The results show that one difference among the subspecies is they produce different bioactive secondary metabolites. B. subtilis subsp. spizizenii is shown conserve the genes to produce mycosubtilin, bacillaene and 3,3'-neotrehalosadiamine. B. subtilis subsp. inaquosorum is shown conserve the genes to produce bacillomycin F, fengycin and an unknown PKS/NRPS cluster. B. subtilis subsp. stercoris is shown conserve the genes to produce fengycin and an unknown PKS/NRPS cluster. While B. subtilis subsp. subtilis is shown to conserve the genes to produce 3,3'-neotrehalosadiamine. In addition, we update the chemotaxonomy and phenotyping to support their promotion to species status.

Iturinic Lipopeptide Diversity in the Bacillus subtilis Species Group - Important Antifungals for Plant Disease Biocontrol Applications

Iturins and closely related lipopeptides constitute a family of antifungal compounds known as iturinic lipopeptides that are produced by species in the Bacillus subtilis group. The compounds that comprise the family are: iturin, bacillomycin D, bacillomycin F, bacillomycin L, mycosubtilin, and mojavensin. These lipopeptides are prominent in many Bacillus strains that have been commercialized as biological control agents against fungal plant pathogens and as plant growth promoters. The compounds are cyclic heptapeptides with a variable length alkyl sidechain, which confers surface activity properties resulting in an affinity for fungal membranes. Above a certain concentration, enough molecules enter the fungal cell membrane to create a pore in the cell wall, which leads to loss of cell contents and cell death. This study identified 330 iturinic lipopeptide clusters in publicly available genomes from the B. subtilis species group. The clusters were subsequently assigned into distinguishable types on the basis of their unique amino acid sequences and then verified by HPLC MS/MS analysis. The results show some lipopeptides are only produced by one species, whereas certain others can produce up to three. In addition, four species previously not known to produce iturinic lipopeptides were identified. The distribution of these compounds among the B. subtilis group species suggests that they play an important role in their speciation and evolution.

Genomic and metabolic features of the Bacillus amyloliquefaciens group- B. amyloliquefaciens, B. velezensis, and B. siamensis- revealed by pan-genome analysis

The genomic and metabolic features of the Bacillus amyloliquefaciens group comprising B. amyloliquefaciens, B. velezensis, and B. siamensis were investigated through a pan-genome analysis combined with an experimental verification of some of the functions identified. All B. amyloliquefaciens group genomes were retrieved from GenBank and their phylogenetic relatedness was subsequently investigated. Genome comparisons of B. amyloliquefaciens, B. siamensis, and B. velezensis showed that their genomic and metabolic features were similar; however species-specific features were also identified. Energy metabolism-related genes are more enriched in B. amyloliquefaciens, whereas secondary metabolite biosynthesis-related genes are enriched in B. velezensis. Compared to B. amyloliquefaciens and B. siamensis, B. velezensis harbors more genes in its core-genome which are involved in the biosynthesis of antimicrobial compounds, as well as genes involved in d-galacturonate and d-fructuronate metabolism. B. amyloliquefaciens, B. siamensis, and B. velezensis all harbor a xanthine oxidase gene cluster (xoABCDE) in their core-genomes that is involved in metabolizing xanthine and uric acid to glycine and oxalureate. A reconstruction of B. amyloliquefaciens group metabolic pathways using their individual pan-genomes revealed that the B. amyloliquefaciens group strains have the ability to metabolize diverse carbon sources aerobically, or anaerobically, and can produce various metabolites such as lactate, ethanol, acetate, CO₂, xylitol, diacetyl, acetoin, and 2,3-butanediol. This study therefore provides insights into the genomic and metabolic features of the B. amyloliquefaciens group.

Pan-genomic and transcriptomic analyses of Leuconostoc mesenteroides provide insights into its genomic and metabolic features and roles in kimchi fermentation

The genomic and metabolic features of Leuconostoc (Leu) mesenteroides were investigated through pan-genomic and transcriptomic analyses. Relatedness analysis of 17 Leu. mesenteroides strains available in GenBank based on 16S rRNA gene sequence, average nucleotide identity, in silico DNA-DNA hybridization, molecular phenotype, and core-genome indicated that Leu. mesenteroides has been separated into different phylogenetic lineages. Pan-genome of Leu. mesenteroides strains, consisting of 999 genes in core-genome, 1,432 genes in accessory-genome, and 754 genes in unique genome, and their COG and KEGG analyses showed that Leu. mesenteroides harbors strain-specifically diverse metabolisms, probably representing high evolutionary genome changes. The reconstruction of fermentative metabolic pathways for Leu. mesenteroides strains showed that Leu. mesenteroides produces various metabolites such as lactate, ethanol, acetate, CO₂, mannitol, diacetyl, acetoin, and 2,3-butanediol through an obligate heterolactic fermentation from various carbohydrates. Fermentative metabolic features of Leu. mesenteroides during kimchi fermentation were investigated through transcriptional analyses for the KEGG pathways and reconstructed metabolic pathways of Leu. mesenteroides using kimchi metatranscriptomic data. This was the first study to investigate the genomic and metabolic features of Leu. mesenteroides through pan-genomic and metatranscriptomic analyses, and may provide insights into its genomic and metabolic features and a better understanding of kimchi fermentations by Leu. mesenteroides.

Comparative Genomic Analysis of Bacillus amyloliquefaciens and Bacillus subtilis Reveals Evolutional Traits for Adaptation to Plant-Associated Habitats

Bacillus subtilis and its sister species B. amyloliquefaciens comprise an evolutionary compact but physiologically versatile group of bacteria that includes strains isolated from diverse habitats. Many of these strains are used as plant growth-promoting rhizobacteria (PGPR) in agriculture and a plant-specialized subspecies of B. amyloliquefaciens-B. amyloliquefaciens subsp. plantarum, has recently been recognized, here we used 31 whole genomes [including two newly sequenced PGPR strains: B. amyloliquefaciens NJN-6 isolated from Musa sp. (banana) and B. subtilis HJ5 from Gossypium sp. (cotton)] to perform comparative analysis and investigate the genomic characteristics and evolution traits of both species in different niches. Phylogenomic analysis indicated that strains isolated from plant-associated (PA) habitats could be distinguished from those from non-plant-associated (nPA) niches in both species. The core genomes of PA strains are more abundant in genes relevant to intermediary metabolism and secondary metabolites biosynthesis as compared with those of nPA strains, and they also possess additional specific genes involved in utilization of plant-derived substrates and synthesis of antibiotics. A further gene gain/loss analysis indicated that only a few of these specific genes (18/192 for B. amyloliquefaciens and 53/688 for B. subtilis) were acquired by PA strains at the initial divergence event, but most were obtained successively by different subgroups of PA stains during the evolutional process. This study demonstrated the genomic differences between PA and nPA B. amyloliquefaciens and B. subtilis from different niches and the involved evolutional traits, and has implications for screening of PGPR strains in agricultural production.

Propionibacterium namnetense sp. nov., isolated from a human bone infection

A polyphasic taxonomic study was performed on two Gram-positive-staining, anaerobic, pleomorphic, rod-shaped strains isolated from human bone and tissue samples. Sequencing of the 16S rRNA genes revealed that the strains belong to a novel species within the genus Propionibacterium, most closely related to Propionibacterium acnes subsp. acnes and Propionibacterium acnes subsp. elongatum with similarity values of 98.4 % and 98.1 %, respectively. In addition, protein-coding genes for rpoB, recA and gyrB clearly separated the novel organism from all species and subspecies of the genus Propionibacterium. However, a DNA-DNA hybridization analysis between the novel organism and the type strain P. acnes ATCC 6919T revealed a value of only 61.1 %. Furthermore, whole genome analysis using the program OrthoANI gave a value of 88.5 %, which is significantly below the cut-off value of 95 % for species delineation. The major fatty acids were iso-C15 : 0, anteiso-C15 : 0 and iso-C17 : 0. The DNA G+C content of the type strain was 59.7 mol%. When taken collectively, phenotypic, molecular genetic, chemotaxonomic and phylogenetic information demonstrate that the organism represents a distinct, albeit close relative of P. acnes On the basis of the results presented, the organism represents a novel member of the genus Propionibacterium for which the name Propionibacterium namnetense sp. nov. is proposed. The type strain is NTS 31307302T (=DSM 29427T=CCUG 66358T).

A molecular phylogenetic framework for Bacillus subtilis using genome sequences and its application to Bacillus subtilis subspecies stecoris strain D7XPN1, an isolate from a commercial food-waste degrading bioreactor

A thermophilic, heterotrophic and facultatively anaerobic bacterium designated strain D7XPN1 was isolated from Baku BakuKing™, a commercial food-waste degrading bioreactor (composter). The strain grew optimally at 45 °C (growth range between 24 and 50 °C) and pH 7 (growth pH range between pH 5 and 9) in Luria Broth supplemented with 0.3 % glucose. Strain D7XPN1 tolerated up to 7 % NaCl and showed amylolytic and xylanolytic activities. 16S rRNA gene analysis placed strain D7XPN1 in the cluster represented by Bacillus subtilis and the genome analysis of the 4.1 Mb genome sequence determined using RAST (Rapid Annotation using Subsystem Technology) indicated a total of 5116 genomic features were present of which 2320 features could be grouped into several subsystem categories. Of these, 615 features were related to carbohydrate metabolism which included a range of enzymes with potential in the biodegradation of food wastes, a property consistent with the ecological habitat of the isolate. ANIb (Average Nucleotide Identity based on BLAST) analysis with 49 Bacillus subtilis genomes indicated that it was distantly related to the three currently taxonomically validated B. subtilis subspecies namely B. subtilis subsp. subtilis (95.6 %), B. subtilis subsp. spizizenii (93 %) and B. subtilis subsp. inaquosorum (92 %) and based on our current knowledge warranted that it be included as a separate cluster together with strain JS which it was closely related (98.69 %). The close relationship of strains D7XPN1 and JS is also supported from our results from electronic DNA–DNA Hybridization (e-DDH) studies. Furthermore, our additional in-depth phylogenomic analyses using three different datasets unequivocally supported the creation of a fourth B. subtilis subspecies to include strains D7XPN1 and JS for which we propose strain D7XPN1^T (=KCTC 33554^T, JCM 30051^T) as the type strain, and designate it as B. subtilis subsp. stecoris.

Spore germination and germinant receptor genes in wild strains of Bacillus subtilis

AIMS

To compare the germination of laboratory and wild strains of Bacillus subtilis.

METHODS AND RESULTS

The spore germination of B. subtilis 168 (subsp. subtilis) was compared with that of the laboratory strain W23 (subsp. spizizenii) and desert-sourced isolates, including one member of subsp. subtilis (RO-NN-1), strains TU-B-10, RO-E-2, N10 and DV1-B-1, (all subsp. spizizenii), the B. mojavensis strain RO-H-1 and a B. subtilis natto strain. All germinated in L-alanine, although some were slower, and some 10-fold less sensitive to germinant. All germinated in calcium dipicolinate (CaDPA). Germination in asparagine, glucose, fructose + KCl was slow and incomplete in many of the strains, and decoating RO-NN-1 and W23 spores did not restore germination rates. Comparing the sequences of B. subtilis strains 168, RO-NN-1, W23, TU-B-10 and DV1-B-1, the operons encoding GerA, B and K germinant receptors were intact, although the two additional operons yndDEF and yfkQRST had suffered deletions or were absent in several spizizenii strains.

CONCLUSIONS

Wild strains possess an efficient germination machinery for L-alanine germination. AGFK germination is often less efficient, the gerB genes more diverged, and the two germinant receptor operons of unknown function have been lost from the genome in many subsp. spizizenii strains.

SIGNIFICANCE AND IMPACT OF THE STUDY

The two major subspecies of B. subtilis have conserved GerA receptor function, confirming its importance, at least in the natural environments of these strains.