Bacillus velezensis K-9 antagonizes S. scabies
The antagonistic activity to S.scabies was screened by the double culture method. The maximum antibacterial rate of strain K-9 was detected for fermentation liquid, and the minimum antibacterial activity was in cell-free filtrate. The antibacterial rate of fermentation broth was 53.92%. The antibacterial rate of strain K-9 was 44.90% and of cell-free filtrate was 27.25% (Fig. 1).
Genomic identification of Bacillus velezensis K-9
The length of 16S rRNA of strain K-9 was 1447 bp. The 16S rRNA gene sequence of strain K-9 was related to Bacillus velezensis, Bacillus amyloliquefaciens and Bacillus subtilis and other strains, sharing more than 97% of the genetic sequences. The results of phylogenetic tree showed that K-9 and Bacillus velezensis SA3 strain branched closely, with 63% similarity rate(Fig. 2).
In the CGView circle of Bacillus velezensis K-9 genome (Fig. 3), corresponding information from outer circle to inner circle is as follows: the first circle (outermost circle), GC content; circle 2, sequencing depth; circle 3, gene category; the fourth circle, CDS information on the positive chain; circle 5, CDS information on the negative chain; the sixth circle, negative chain COG functional classification annotation information; circle 7 (innermost circle), plus link COG functional classification annotation information.
The genome size was 3 891 530 bp. The GC content was 46.45%, and there were 11 contigs. The number of tRNA and rRNA was 79 and 10, respectively. There were 3915 protein-coding genes in Bacillus velezensis K-9, and the proportion of protein -coding genes was 89.87%(Table 1).
Table 1
Complete genome features of Bacillus velezensisK-9
Property
|
Value
|
Sequencing platform
|
PacBio RS II
|
Genome size (bp)
|
3891530
|
Contig numbers
|
11
|
DNA G + C (%)
|
46.45
|
Number of tRNA genes
|
79
|
Number of rRNA genes
|
10
|
Number of protein coding genes
|
3915
|
Ratio of protein coding genes
|
89.87
|
Annotation of the Clusters of Orthologous Groups (COG) of proteins has been developed based on protein sequence alignment based on gene sequences, with each COG representing a gene family15–16. The COG classification number of Bacillus velezensis K-9 genome was 4, and the number of COG types was 20. The number of genes with COG annotation was 1238, accounting for 33.05% of all genes. As can be seen from COG annotation information (Fig. 4), apart from the predicted gene function represented by "R" and the unknown function represented by "S", "K" represents the transcription gene family with the largest number of members (239), followed by "E" representing the amino acid transfer and metabolism family with 238 members.
Through genomic GO (Gene Ontology) annotations, biological terms can be standardized for easy communication. Statistical results of GO annotation of Bacillus velezensis K-9 genome are shown in Fig. 5. In Bacillus velezensis K-9 genome, the number of GO classifications was 3, and the number of genes involved in annotation was 3004, accounting for 80.19% of all genes.
In the KEGG (Kyoto Encyclopedia of Genes and Genomes) annotations of Bacillus velezensis K-9 genome, the gene functions were divided into six major categories. Among them, 82 genes were annotated to cellular processes, 374 genes to environmental information processing, 240 to genetic information processing, 2066 to metabolism, 68 to human diseases, and 28 genes to organismal systems. Detailed breakdown of each category is shown in Fig. 6.
As shown in Fig. 7, among the genes coding for fermentation activity, the glycoside hydrolase genes were the most abundant in Bacillus velezensis K-9, accounting for 29.65%, followed by 48 genes of glycosyl transferases, accounting for 27.90%. The gene numbers for carbohydrate esterases, carbohydrate-binding modules, auxiliary oxidoreductases, and polysaccharide lyases were 39 (22.67%), 18 (10.47%), 11 (6.40%), and 5 (2.91%), respectively. The main components of plant diseases and residues in soil are sugars and proteins, and the strain K-9 contains β-glucosidase (EC3.2.1.21), chitin (β-1,6-glucanosyltransferase, EC 2.4.1.-), α-amylase (EC 3.2.1.1), cellulase (EC 3.2.1.4) and other enzymes. The biocontrol mechanism of strain K-9 was further analyzed from the perspective of gene.
AntiSMASH6.0.1 software was used to predict the secondary metabolites of Bacillus velezensis K-9. The results showed that there were 12 secondary metabolite synthesis gene clusters in the genome of Bacillus velezensis K-9. Except for the unknown metabolites of the third polyketone (T3PKS), terpene and lanthipeptide-class-II, the other eight gene clusters are associated with biosynthesis of secondary metabolites that have antibacterial activity and are common in Bacillus: Clusters 1 and 5, polyketone synthesis genes; Cluster 4, synthesis gene cluster of Camelina; Cluster 6, macrolide synthesis; Cluster 9 butamectin synthesis; Cluster 10, lysobacterin synthesis; Cluster 11, siderophores synthesis; and Cluster 12, surfactant synthesis. Gene clusters 4, 11 and 12 are associated with synthesis of secondary metabolites through non-ribosomal pathways (Table 2, Fig. 8-A).
Table 2
Prediction of secondary metabolites of Bacillus velezensisK-9 strain
Gene clustera
|
Typeག
|
Locationགྷ
|
Secondary metabolitesང
|
Functione
|
Similarityf
|
1
|
TransAT-PKS
|
546640–638999
|
Difficidin
|
Anti-bacterial
|
100%
|
2
|
T3PKS
|
767634–808360
|
—
|
—
|
—
|
3
|
Terpene
|
873700–893826
|
—
|
—
|
—
|
4
|
NRPS
|
922563–1056873
|
Fengycin
|
Anti-fungal
|
100%
|
5
|
TransAT-PKS
|
1130615–1231180
|
Bacillaene
|
Anti-fungal, anti-bacterial
|
100%
|
6
|
TransAT-PKS
|
1450368–1538601
|
Macrolactin
|
Anti-fungal, anti-bacterial, anti-virus
|
100%
|
7
|
Lanthipeptide-class-ii
|
1705164–1734052
|
—
|
—
|
—
|
8
|
Terpene
|
1854545–1875285
|
—
|
—
|
—
|
9
|
PKS-like
|
1957329–1998573
|
Butirosin
|
Anti-bacterial
|
7%
|
10
|
Other
|
83195–124613
|
Bacilysin
|
Anti-bacterial
|
100%
|
11
|
NRPS
|
660923–712714
|
Bacillibactin
|
Accumulate and take up iron ions
|
100%
|
12
|
NRPS
|
199691–265098
|
Surfactin
|
Anti-virus, anti-mycoplasma,anti-tumour
|
82%
|
a : the secondary metabolite synthesis gene cluster annotated with antiSMASH6.0.1, b : the gene cluster type,c : the location of gene clusters in the strain genome, d : the secondary metabolite that may be produced based on the gene cluster, e : the bioactive function of metabolites synthesized by gene cluster, f : similarity to known gene clusters,“—”: indicates unknown.
According to the bacteriostatic activity test results of crude protein extract of Bacillus velezensis K-9 against Scab and protein gene analysis (Figure 8-B), It is presumed to be the main antibacterial substance of the K-9 against Scab.
Subsequently, the crude protein extract was detected by SDS-PAGE. As can be seen from Fig. 9, the protein subunits were distributed below 100kDa, and there were many bands. The gel electrophoresis pattern is clear and the bands are rich, which can be used for subsequent LC-MS/MS metabolic test.
Each protein was enzymatically digested to produce peptides of different lengths, which were then identified by LC-MS/MS(Table 3). Based on Uniprot-Bacillus Velezensis database, 4001 peptides and 500 proteins were found. The protein information related to antimicrobial activity is shown in Table 3. Polyketide synthetase (gene name: PksN、PksL), plastin synthetase, subtilisin, flagellar biosynthetic protein (gene name: FliP) and so on were identified. Polyketide synthetase, coptigenin synthetase and non-ribosomal peptide synthetase have been identified from crude proteins to catalyze the synthesis of antibacterial substances. The genes of these three enzymes are often used as indicator genes to screen for the production of novel structural natural active compounds17–18. Subtilisin, an antimicrobial substance synthesized by the ribosome pathway, has the characteristics of good stability and strong acid-base adaptability, and has antibacterial effect on most Gram-positive bacteria19. FliP genes play a role in the flagellum-specific transport system. Flagellin is a typical virulence factor, which helps bacterial infection and colonization and plays an important role in bacteriostatic activity20. Chitin binding protein can help improve plant defense system, thereby enhancing plant disease resistance21. Chitosanase is a hydrolytic enzyme that exclusively degrades chitosan. However, it is also known as Pathogenesis related protein, which can increase the disease resistance of plants. Bacillus Velezensis K−9 inhibited the activity of potato scab pathogen by producing these five antibacterial substances.
Table 3
Bacillus velezensis K-9 results of qualitative analysis of protein substance by LS-MS
Login number
|
Protein information
|
The name of the protein
|
A0A7W4QCN4
|
Polyketide synthase PksN OS = Bacillus velezensis OX = 492670 GN = pksN_2 PE = 4 SV = 1
|
Polyketide synthase PksN
|
A0A7W4LRG1
|
Polyketide synthase PksL OS = Bacillus velezensis OX = 492670 GN = pksL_5 PE = 4 SV = 1
|
Polyketide synthase PksL
|
A0A1D9PPV3
|
Bacilysin biosynthesis protein BacB OS = Bacillus velezensis OX = 492670 GN = bacB PE = 4 SV = 1
|
Bacilysin biosynthesis protein BacB
|
L7WLS3
|
Subtilisin OS = Bacillus velezensis OX = 492670 GN = kerBPN PE = 3 SV = 1
|
Subtilisin
|
A0A8G1PD45
|
Flagellar biosynthetic protein FliP OS = Bacillus velezensis OX = 492670 GN = fliP PE = 3 SV = 1
|
Flagellar biosynthetic protein FliP
|
A0A7S8F0Q0
|
Fengycin synthetase B (Fragment) OS = Bacillus velezensis OX = 492670 GN = fenB PE = 4 SV = 1
|
Fengycin synthetase B
|
A0A6H3B0Z0
|
Non-ribosomal peptide synthetase OS = Bacillus velezensis OX = 492670 GN = CHR37_17550 PE = 3 SV = 1
|
Non-ribosomal peptide synthetase NRPS
|
A0A8G1LH67
|
Flagellin OS = Bacillus velezensis OX = 492670 GN = K4L72_17430 PE = 3 SV = 1
|
Flagellin
|
A0A8F0EXT8
|
Chitosanase OS = Bacillus velezensis OX = 492670 GN = csn PE = 3 SV = 1
|
Chitosanase
|
A0A1D9PQS1
|
Chitin-binding protein OS = Bacillus velezensis OX = 492670 GN = gbpA PE = 4 SV = 1
|
Chitin-binding protein
|
Note: OS is the Latin name of the species, GN is the gene name, PE is the protein existence certificate (PE = 1: experimental evidence at protein level, PE = 2: experimental evidence at tranlevel, PE = 3: protein inferred from homology, PE = 4: protein predicted, PE = 5: protein uncertain), SV is the protein sequence version. |