Bacillus velezensis FZB42 in 2018: The Gram-Positive Model Strain for Plant Growth Promotion and Biocontrol

1Bacteria that are associated with plant roots and exert beneficial effects on plant development are referred to as plant-growth-promoting rhizobacteria (PGPR; ). It is well accepted today, that numerous PGPR are also enabled to control plant diseases.

2Main subject of present and past research about microbial inoculants with beneficial action on plant health and growth are plant-associated representatives of the bacterial genus Pseudomonas, known as strong and persistent colonizer of plant roots. However, its commercial use is limited by difficulties in preparing stable and long-living bioformulations. As early as at the end of the 19th century a bacterial soil-fertilizing preparation Alinit® consisting of spores of the soil bacterium Bacillus ellenbachensis, later reclassified as Bacillus subtilis, was introduced by the German landowner Albert Caron (1853–1933) on his estate in Ellenbach. Alinit was marketed as “bacteriological fertilizer for the inoculation of cereals” by “Farbenfabriken former Friedrich Bayer,” the later Bayer AG, in Elberfeld, Germany. The history of these early attempts in using bacterial inoculants is comprehensively described by. After a long period of silence, the plant-growth-promoting effect of Bacillus spp. was rediscovered in. Today, formulations based on plant-beneficial endospore-forming Bacilli are by far the most widely used agents on the biopesticide market. Especially, members of the B. subtilis species complex (rRNA group 1) which includes at present more than 20 closely related species, and, to a minor extent, of the genus Paenibacillus spp., are able to suppress efficiently plant pathogens, such as viruses, bacteria, fungi and nematodes in vicinity of plant roots. This review describes the current ‘state of the art’ of the model strain for PGPR – and biocontrol, Bacillus velezensis FZB42, and the integrative data bank ‘AmyloWiki,’ recently established for this bacterium.

3FZB42 (=BGSC 10A6, DSM23117), the prototype of gram-positive bacteria with phytostimulatory and biocontrol action, has been genome sequenced in and is subject of intensive research. Since its isolation from beet rhizosphere more than 140 articles about FZB42 have been published1. FZB42 and its closely related ‘cousin’ FZB24, are successfully used as biofertilizer and biocontrol bacteria in agriculture being especially efficient against fungal and bacterial pathogens2. Beneficial effects of FZB42/FZB24 on plant growth and disease suppression in field trials were reported for potato, cotton, strawberry, wheat, lettuce, and tomato, for example.

4In past, FZB42 and related phytostimulatory Bacilli were subjects of intensive efforts to clarify their taxonomic position. The group of plant-associated, endo-spore forming rhizobacteria is known as member of the B. subtilis species complex, which included originally B. subtilis, B. licheniformis, and B. pumilus. In 1987, the species B. amyloliquefaciens was added, and FZB42 and some other biocontrol bacteria were found as belong to this species. By taking advantage of availability of an increasing number of genome sequences, we distinguished two subspecies: B. amyloliquefaciens subsp. amyloliquefaciens (type strain DSM7T) and B. amyloliquefaciens subsp. plantarum (type strain FZB42T). According to extended phylogenomic analysis B. amyloliquefaciens subsp. plantarum was shown as a later heterotypic synonym of B. velezensis, Recently, we proposed to establish an “operational group B. amyloliquefaciens,” which includes B. amyloliquefaciens, known for its ability to produce industrial enzymes (amylases, glucanases and proteases), B. siamensis, mainly occurring in Asian food, and PGPR B. velezensis, the main source for bioformulations increasingly used in agriculture for protecting plant health and to stimulate plant growth (, Figure 1).

5Biocontrol effects exerted by B. velezensis FZB42 and other antagonistic acting Bacilli are due to different mechanisms: besides direct antibiosis and competition by secretion of a spectrum of secondary metabolites in the rhizosphere, the beneficial action on the host-plant microbiome, and stimulation of plant induced systemic resistance (ISR,; ) are of similar importance.

6Remarkably, in contrast to Gram-negative biocontrol bacteria and fungal plant pathogens, application of FZB42 did not lead to durable changes in composition of rhizosphere microbial community. Moreover, application of FZB42 was shown to compensate negative changes within composition of the root microbiome caused by plant pathogens.

7Induced systemic resistance is triggered by a range of secondary metabolites, which are called ‘elicitors.’ Different signaling pathways, such as jasmonic acid (JA), ethylene (ET), and salicylic acid (SA) are activated to induce plant resistance. Mutant strains of FZB42, devoid in synthesis of surfactin (srf), were found impaired in triggering of JA/ET dependent ISR in lettuce plants, when challenged with plant pathogen Rhizoctonia solani. The lower expression of the JA/ET-inducible plant defensin factor (PDF1.2) in a sfp mutant strain, completely devoid in non-ribosomal synthesis of lipopeptides and polyketides, compared to the srf mutant strain, only impaired in surfactin synthesis, suggests that secondary metabolites other than surfactin might also trigger plant response.

8Gray leaf spot disease caused by Magnaporthe oryzae is a serious disease in perennial ryegrass (Lolium perenne). A mutant strain of FZB42 (AK3) only able to produce surfactin but no other lipopeptides such as bacillomycin D, and fengycin was shown to induce systemic resistance (ISR). Similarly, treatment with crude surfactin suppressed the disease in perennial ryegrass. ISR defense response was found connected with enhanced hydrogen peroxide (H2O2) development, elevated cell wall/apoplastic peroxidase activity, and deposition of callose and phenolic/polyphenolic compounds. Moreover, a hypersensitive response reaction and enhanced expression of different defense factors, such as peroxidase, oxalate oxidase, phenylalanine ammonia lyase, lipoxygenase, and defensins were caused by surfactin and also the surfactin producing mutant strain.

9Recent studies performed with mutant strains of B. velezensis SQR9, which is closely related with FZB42, revealed that non-ribosomal synthesized lipopeptides fengycin and bacillomycinD, the non-ribosomal synthesized polyketides macrolactin, difficidin, and bacillaene, the dipeptide bacilysin, exopolysaccharides, and volatile organic compounds (VOCs) contribute to ISR response in Arabidopsis plantlets after infection with plant pathogens Pseudomonas syringae pv. tomato and Botrytis cinerea.

10Volatile organic compounds produced by B. velezensis GB03 have been reported to trigger synthesis of ET/JA-responsive plant defense gene PDF1.2. Thirteen VOCs produced by FZB42 were identified using gas chromatography-mass spectrometry analysis. A direct effect against plant pathogens was registered: benzaldehyde, 1,2-benzisothiazol-3(2 H)-one and 1,3-butadiene significantly inhibited the colony size, cell viability, and motility of Ralstonia solanacearum, the causative agent of bacterial wilt in a wide variety of potential host plants. Furthermore, transcription of type III (T3SS) and type IV secretion (T4SS) systems were down regulated. In addition, synthesis of other genes contributing to pathogenicity, such as eps-genes responsible for extracellular polysaccharides, and genes involved in chemotaxis (motA, fliT) were found repressed. Simultaneously, the VOCs significantly up-regulated the expression of plant genes related to wilt resistance and pathogen defense. Over-expression of plant defense genes EDS1 and NPR1 suggested that the SA pathway is involved in the ISR response elicited by surfactin.

11A recent analysis performed with FZB42 VOCs confirmed that signal pathways involved in plant systemic resistance were positively affected. JA response (VSP1 and PDF1.2) and SA response genes (PR1 and FMO1) were triggered in Arabidopsis plantlets after incubation with the volatiles. Noteworthy, defense against nematodes were elicited by volatiles in Arabidopsis roots.

12An interesting mechanism of FZB42 to avoid leaf pathogen infection has been recently described. The foliar pathogen Phytophthora nicotianae is able to penetrate inside of plant tissues by using natural entry sites, such as stomata. Recently it was shown that colonizing of plant roots by FZB42 restricted entry of the pathogen into leave tissues of Nicotiana benthamiana. It was found that FZB52 turned on the abscisic acid (ABA) and SA-regulated pathways to induce stomatal closure after pathogen infection. In addition, it was shown, that several SA- and JA/ET-responsive genes in the leaves became activated in presence of FZB42, suggesting that these signaling pathways are also contributing to plant defenses against P. nicotianae.

13Besides their indirect action against pathogens via triggering of ISR, polyketides and lipopeptides act directly against bacterial and fungal plant pathogens. They comprise two families of secondary metabolites non-ribosomally synthesized by multimodular enzymes, polyketide synthases (PKSs) and Peptide synthetases (NRPS), acting in assembly line arrays. The monomeric building blocks are either organic acids (polyketides) or amino acids (lipopeptides), respectively. Their synthesis is depending on an enzyme (Sfp) that transfers 4′-phosphopantheine from coenzyme A to the carrier proteins of nascent peptide or polyketide chains. In Bacilli, e.g., FZB42, a special class of PKSs that lacks the cognate AT domain and require a discrete AT enzyme acting iteratively in trans (trans AT) was detected. The broadly conserved antiterminator protein LoaP (Nus G family) was identified as regulator of macrolactin and difficidin gene clusters in B. velezensis FZB42 on the level of transcription elongation. Unfortunately, structural instability of these polyketides excluded their use as antibacterial agents.

14Lipopeptides are another important class of secondary metabolites, also non-ribosomally synthesized by giant multifunctional enzymes (peptide synthetases, NRPS). Similar to PKS, three catalytic domains are involved in each elongation cycle: (1) The A-domain (adenylation domain) select its cognate amino acid; (2) The PCP domain (peptidyl-carrier domain) is equipped with a PPan prosthetic group to which the adenylated amino acid substrate is transferred and bound as thioester; (3) The condensation domain (C-domain) catalyzes formation of a new peptide bond. The lipopeptide bacillomycin D is an efficient antifungal compound produced by FZB42. Its 50% effective concentration against the fungal pathogen Fusarium graminearum was determined to be approximately 30 μg/ml. Bacillomycin D induced morphological changes in the plasma membranes and cell walls of F. graminearum hyphae and conidia. Furthermore, bacillomycin D induced the accumulation of reactive oxygen species and caused cell death in F. graminearum hyphae and conidia. Bacillomycin D suppresses F. graminearum on corn silks, wheat seedlings, and wheat heads.

15Today, B. subtilis is considered as being a plant-associated bacterium. A direct comparison between the genomes of B. subtilis 168 and B. velezensis FZB42 (Table 1) revealed that 534 FZB42 genes are not occurring in B. subtilis 168, but 3158 genes are shared between both species. By contrast, there are only 423 singletons defined for FZB42 vs. Bacillus subtilis 168. In this context one has to mention, that the singleton numbers don’t correspond to the numbers in the Venn diagram. The Venn diagram (Figure 2) shows the numbers of reciprocal best hits between subsets of genomes. However, a gene without reciprocal best hit to another genome does not necessarily have to be a singleton. A singleton is defined as a gene without any hit against any other genome than the own one.

16Many genes, essential for a plant-associated lifestyle, are shared between B. subtilis 168 and FZB42 as well. Relevant examples are YfmS, a chemotaxis sensory transducer, which is involved in plant root colonization, and BlrA (formerly YtvA) a blue light receptor related to plant phototropins. However, due to a century of ‘domestication’ under laboratory conditions, the type strain B. subtilis 168 has lost its ability to colonize roots and to control plant diseases. Its ability to form biofilms on solid surfaces (e.g., rhizoplane) is attenuated by several mutations detected in the genes sfp (necessary for production of lipopeptides and polyketides), epsC (required for extracellular polysaccharide synthesis), swrA (essential for swarming differentiation on solid surfaces), and degQ, which stimulates phosphorylation of DegU. By contrast, the closely related wild type B. subtilis 3610 forms robust biofilms and is able to produce antimicrobial compounds (Table 1). It was shown that by introducing wild type alleles of these four genes and the spo0F phosphatase encoding rapP gene, residing on a large plasmid occurring in B. subtilis 3610 but not in B. subtilis 168, the laboratory strain 168 forms biofilms which are essentially the same as in 3610. This demonstrates that domestication of B. subtilis 168 is only due to the four gene mutations mentioned above and loss of the plasmid occurring in strain 3610. Notably, FZB42 does not harbor a rapP containing plasmid, but is able to produce robust biofilm similar to B. subtilis 3610.

17FZB42 releases several cellulases and hemicellulases degrading the external cellulosic and hemicellulosic substrates present in plant cell walls. Final products of enzymatic hydrolysis are free oligosaccharides, which act as elicitors of plant defense. Some genes encoding extracellular hydrolases, such as amyE (alpha-amylase), eglS (endo-1,4-β-glucanase), and xynA (xylanase) occurred only in the plant-associated representatives of the ‘B. amyloliquefaciens operational group’ but not in their soil-associated counterparts. Similarly, an operon involved in xylan degradation (xylA, xynP, xynB, xylR) is present in B. subtilis 168 and B. velezensis FZB42 but not in B. amyloliquefaciens DSM7T suggesting that both strains have in common some genes involved in plant macromolecule degradation.

18Bacillus velezensis harbored additional genes involved in hexuronate (galacturonate and fructuronate) degradation. Three genes were found unique for B. velezensis FZB42 and other members of this species: kdgK1, (2-dehydro-3-deoxygluconokinase), kdgA (2-dehydro-3-deoxyphosphogluconate aldolase), and the transcription regulator kdgR. They are part of the six-gene kdgKAR operon. In addition yjmD, a gene with putative galacticol-1-phosphate dehydrogenase function and two further genes: uxuA encoding mannonate dehydratase, and uxuB encoding mannonate oxidoreductase are part of the six-gene transcription unit. A second operon, containing the genes uxaC, uxaB, and uxaA encoding enzymes for metabolizing different hexuronates to D-altronate and D-fructuronate, occurs distantly from the kdgAR operon. In Escherichia coli K12 UxuA, KdgK, and KdgA are involved in a degradative pathway of aldohexuronates. Whilst the complete biochemical pathway from galacturonate to KDG is present, no gene encoding D-glucuronate isomerase was detected, suggesting that B. velezensis is not able to metabolize D-glucuronate.

19Nearly 10% of the FZB42 genome is involved in synthesizing antimicrobial compounds, such as the polyketides bacillaene, macrolactin and difficidin and the lipopeptides surfactin, bacillomycin D and fengycin. In total, the FZB42 genome harbors 13 gene clusters devoted to non-ribosomal and ribosomal synthesis of secondary metabolites with putative antimicrobial action. In two cases, in the nrs gene cluster and in the type III polyketide gene cluster their products are not identified till now (Table 1). Similar to B. subtilis 168T, the genome of the non-plant associated soil bacterium B. amyloliquefaciens DSM7T possesses a much lower number of gene clusters involved in synthesis of antimicrobial compounds than FZB42 (Table 1).

20Notably, the gene clusters involved in non-ribosomal synthesis of the antifungal lipopeptides bacillomycin D and fengycin, and the polyketides difficidin and macrolactin are missing or fragmentary in DSM7T and other representatives of B. amyloliquefaciens suggesting that synthesis of these secondary metabolites might be important for the plant associated life style. Five out of a total of 13 gene clusters are located within variable regions of the FZB42 chromosome, suggesting that they might be acquired via horizontal gene transfer. Most of them (bacillomycin D, macrolactin, difficidin, plantazolicin, and the orphan nrsA-F gene cluster) are without counterpart in DSM7T and B. subtilis 168T. Moreover, it has been shown experimentally that DSM7T, due to a deletion in the fengycin gene cluster, is unable to produce fengycin, notably the gene cluster for synthesis of iturinA is present in the DSM7T genome (Table 1).

21Besides type I PKS also genes encoding type III polyketide synthases are present in the genome of FZB42. By contrast to type I PKSs, type III PKSs catalyze priming, extension, and cyclization reactions iteratively to form a huge array of different polyketide products. In B. subtilis gene products of bspA-bspB operon were functionally characterized, and found to be involved in synthesis of triketide pyrones. The type III PKS BspA catalyzes synthesis of triketide pyrones and BspB (YpbQ) is a methyltransferase catalyzing its posttranslational modification to alkylpyrones ethers. However, their biological role needs further elucidation. Orthologs of bspA and bspB are present in FZB42 and DSM7T (Table 1).

22Another group of secondary metabolites are bacteriocins, which represent a class of post-translationally modified peptide antibiotics. Together with peptides without antibiotic activity, they are generally termed RiPPs (ribosomally synthesized and post-translationally modified peptides). RiPP precursor peptides are usually bipartite, being composed of an N-terminal leader and C-terminal core regions. RiPP precursor peptides can undergo extensive posttranslational modification, yielding structurally and functionally diverse products. In recent years, two RiPPs with antibacterial activity (bacteriocins) were identified in FZB42: plantazolicin and amylocyclicin.

23An antibacterial substance still produced by a FZB42 mutant strain, unable to synthesize non-ribosomally any antimicrobial compound, was identified together with the gene cluster responsible for its biosynthesis. The pzn genes cluster encodes a small precursor peptide PznA that is post-translationally modified to contain thiazole and oxazole heterocycles. These rings are derived from Cys and Ser/Thr residues through the action of a modifying “BCD” synthetase complex, which consists of a cyclodehydratase (C), a dehydrogenase (B), and a docking protein (D). After modification and processing of the precursor peptide plantazolicin contains an unusual number of thiazoles and oxazoles. The structure variant plantazolicin A inhibits selectively Bacillus anthracis, and is efficient against plant pathogenic nematodes, whilst the precursor molecule PZNB is inactive.

24The head-to-tail cyclized bacteriocin amylocyclicin was firstly described in B. amyloliquefaciens FZB42. Circular bacteriocins are non-lanthionine containing bacteriocins with antimicrobial activity against Gram-positive food-borne pathogens. Amylocyclicin was highly efficient against Bacilli, especially against a sigW mutant of B. subtilis (Y2). An orthologous gene cluster was also detected in B. amyloliquefaciens DSM7T (Table 1).

25Lci was reported as an antimicrobial peptide synthesized by a B. subtilis strain with strong antimicrobial activity against plant pathogens, e.g., Xanthomonas campestris pv. oryzae and Pseudomonas solanacearum PE1. Its solution structure has a novel topology, containing a four-strand antiparallel β-sheet as the dominant secondary structure. The gene is not present in the B. subtilis 168 genome, but was detected in FZB42 and B. amyloliquefaciens DSM7T (Table 1). Lci was found highly expressed in FZB42 biofilms.

26Nowadays, global gene expression studies were increasingly performed to enlarge our knowledge base about effect of plants on gene expression in Gram-positive plant associated bacteria. The first combined transcriptome- and proteome analysis in Bacillus, using both, DNA-microarrays and 2-D protein gel electrophoresis, was conducted with B. subtilis 168. Plant-bacteria interactions were studied with B. subtilis OKB105 in presence of rice seedlings. Transcriptome analysis revealed that expression of 176 bacterial genes was affected by the host plant.

27In this context several studies were performed with FZB42, too. Transcription of many genes involved in carbon and amino acid metabolism was turned on, when maize root exudates were added to FZB42 cells growing in planktonic culture suggesting that nutrients present in root exudates are utilized by bacteria cells. Dependency of FZB42 from nutrient sources present in root exudates was corroborated in a second transcriptome study performed with DNA-microarrays. In this case root exudates with different composition obtained from maize plantlets growing under stress conditions (N, P, Fe, and K limitation) were used. In case of root exudates obtained from N-deprived maize plantlets containing decreased amounts of aspartate, valine and glutamate, FZB42 cells were found to be downregulated in transcription of genes involved in protein synthesis indicating a general stress response. By contrast, P-limited root exudates led to enhanced transcription of FZB42 genes involved in motility and chemotaxis, possibly suggesting a chemotactic response toward carbohydrates in root exudates. Transcriptional profiling via RNA-sequencing in the taxonomically related B. velezensis SQR9 revealed that maize root exudates stimulated at first expression of metabolism-relevant genes and then genes involved in production of the extracellular matrix.

28Response of FZB42 on maize root exudates during late exponential and stationary growth phase was also investigated on the level of protein synthesis applying 2-D gel electrophoresis and MALDI TOF MS for protein identification. Elicitors of plant innate immunity such as flagellins, elongation factor Tu, and cold shock proteins were detected in the extracellular fluid. Corresponding to the results obtained in our transcriptome studies, we found that the expression of genes involved in utilization of nutrients and transport was enhanced in presence of root exudates. The protein with the highest secretion in presence of maize root exudates was acetolactate synthase AlsS, an enzyme involved in post-exponential phase synthesis of acetoin and 2,3 butandiol.

29On the other hand, plants are also affected in their gene expression, when colonized by bacteria including representatives of the B. amyloliquefaciens operational group. Transcript analysis of rape seedlings confronted with a root-colonizing B. velezensis strain revealed that gene expression was more affected in leaves than in roots. Altogether the treatment caused a metabolic reprogramming in plant leaves. Similar effects on plant gene expression were reported for root-colonizing B. subtilis FB17. A microarray study performed with Arabidopsis plantlets exposed to FB17 showed that expression of auxin-regulated genes and genes involved in metabolism, stress response and plant defense were up-regulated. Some Arabidopsis mutants deficient in three of the up-regulated genes, were found less colonized by FB17. Further papers reporting about triggering of ISR response in plants by lipopeptides and VOCs from B. velezensis were already discussed in a previous section.

30Another study performed with FZB42 revealed that gene expression is dependent on life style. Ability to form biofilms is essential for colonizing plant root surfaces. Differential gene expression suggested that under biofilm-forming conditions transcription of 331 genes was increased and of 230 genes was decreased.

31The differential RNA-sequencing (dRNA-seq) technology was employed to unveil the structure of the FZB42 transcriptome. The unique feature of this technique is that two libraries split from the same RNA sample are compared. One library is subjected to terminator exonuclease that preferentially degrades processed RNAs with 5′-monophosphate group, thus primary transcripts with 5′-triphosphate group are enriched in relative terms. Applying this method, we obtained the first global transcription start sites (TSSs) map of a PGPR Bacillus species. We determined a comprehensive transcriptome profile for FZB42 by identifying 4,877 TSSs for protein-coding genes. This includes >2,000 primary TSSs, >700 secondary TSSs, and nearly 200 orphan TSSs. The primary TSSs have been identified for 60% of all FZB42 genes. In addition, >1,300 internal TSSs and >1,400 antisense TSSs were also identified. A lot of coding genes were shown to be transcribed from multiple TSSs and perhaps own different UTRs. Some mRNAs contained overlapped transcripts. The global charting of FZB42 TSSs can favor the identification of promoter regions, cis-acting regulatory elements, and cognate transcriptional regulators.

32By applying the dRNA-seq technique differentially expressed genes under different growth conditions were identified. For example, a large group of genes that are specifically regulated by root exudates during stationary growth were identified. The results obtained extended and corroborated our previous results obtained by using microarrays. Knowledge of the genes affected in their expression by plant root exudates contributes to our understanding of rhizobacterial physiology and its interaction with their host plants. They are listed as ‘Interaction with plants’ in AmyloWiki3. Moreover, this study allowed us to propose 46 previously unrecognized genes. 78 polycistronic transcripts covering 210 genes were identified and 10 previously mis-annotated genes were corrected.

33Over the last decade, a growing number of non-coding regulatory small RNAs (sRNAs) have been identified in bacteria, although the functions of most of them are still unknown. Most of sRNAs do not encode a protein, but function as an RNA regulator directly targeting multiple mRNAs. It is revealed that many sRNAs contribute to bacterial adaptation to changing environments and growth conditions, therefore it is feasible to expect that sRNAs may also coordinate mutual effects of rhizobacteria on plants.

34Besides graphing the profile of expressed protein-coding genes, dRNA-seq technology also offers a possibility to identify genome-wide sRNAs. We detected hundreds of non-coding RNAs in FZB42, including 136 antisense RNAs, 53 cis-encoded leader sequence or riboswitches, and 86 sRNA candidates. Among them 21 sRNAs were further validated by Northern blotting. According to their gene positions, the majority of the sRNAs perhaps act in-trans targeting the mRNAs encoded from a distant locus. Generally, sRNAs often binds to their target mRNAs, at 5′UTR in many cases, and thus modulate mRNA translation. Since the genome-wide TSS annotation of FZB42 informs about potential sRNA target sites of mRNAs, our study has provided a valuable basis for studying rhizobacterial sRNA regulation.

35The function of the identified sRNAs has not been characterized in detail. However, some of the sRNAs were found related to a specific growth phase or to respond to environmental cues (soil extract or maize root exudates). Furthermore, one sRNA was found to be involved in Bacillus sporulation and biofilm formation (data not shown). Since, sRNAs are more studied in Gram-negative than in Gram-positive bacteria, systematic detection of sRNAs in FZB42 extends our knowledge base about plant-associated Gram-positive bacteria, especially to rhizobacteria–plant interactions.

36In recent years post-translational modifications (PTM) of proteins, such as protein phosphorylation, acetylation, methylation, and succinylation, attracted increasing attention due to their important physiological significance in organisms. Whereas most studies of PTM were performed in eukaryotic cells, nowadays the role of PTM in prokaryotes is increasingly investigated. Acetylation of lysine residues in FZB42 was studied using a combination of immune-affinity purification and high- resolution LC-MS/MS. A total of 3,268 acetylated lysine residues were detected in 1254 proteins, accounting for 32.9% of the entire proteins of FZB42. Remarkably, a high proportion (71.1 and 78.6%) of the proteins related to the synthesis of polyketides and lipopeptides were found acetylated. The finding implies an important role of lysine acetylation in the regulation of FZB42 antibiotic biosynthesis.

37Using a similar technique, we profiled lysine malonylation of proteins in FZB42. In total, we identified 809 malonyl-lysine sites in 382 proteins (Figure 3). Lysine malonylation targets the proteins implicated in a wide range of biological functions, such as fatty acid biosynthesis and metabolism, central carbon metabolism, translation processes, and NAD(P) binding. A group of proteins involved in bacterium-plant interaction was also malonylated. Moreover, malonylation seems to occur on proteins with higher surface accessibility, although the significance of the site preference remains unclear. Similar to lysine acetylation, 33 polyketide synthases (PKS) and polypeptide synthetases (NRPS) involved in non-ribosomal synthesis of bacillaene, difficidin, macrolactin, and bacillomycinD, fengycin and surfactin, were found highly malonylated. They account for 8.6% of all malonylated proteins. The PKSs and NRPSs possessed 128 malonylation sites, averagely 3.8 sites per protein, which is significantly higher than the mean of 2.1 malonylation sites per protein. The polyketide synthases, BmyA, BaeM, BaeN, and BaeR contain more than 10 malonylation sites. BaeR is the most highly malonylated protein carrying 17 malonylation sites.

38Together with the data obtained for acetylation, the high malonylation rate of PKSs and NRPSs indicates a potential effect of protein modification on biosynthesis of antibiotics in FZB42. Better understanding of the underlying mechanism of how PTM affects PKSs and NRPSs may facilitate the development of FZB42 antibiotic production and application.

39With the increasing reception of FZB42 as a model organism for Gram-positive PGPR, and in order to celebrate its whole genome sequencing around 10 years ago, we have established an integrated database ‘AmyloWiki’4 for collecting and gathering all the information known to date about this bacterium (Figure 4). More than 140 articles about FZB42 can be found in AmyloWiki5 and are in part assigned to the corresponding genes/proteins. AmyloWiki centers the achievement of FZB42 studies till now including diverse information such as its 3979 genes, its transcriptome structure, protein regulators and their targets. 595 genes of FZB42 involved in plant-bacteria interactions were listed6. It informs also about recently identified sRNA genes and post-translational modification sites (see previous sections). A growing list of FZB42-site directed mutant strains, available for scientific community, is also presented. AmyloWiki shares some features with SubtiWiki, the popular database for B. subtilis 168; however, specific features of FZB42 such as genes not occurring in B. subtilis 168, and genes involved in antagonism against plant pathogens and plant-microbe interaction, are highlighted in AmyloWiki. To facilitate communication and information exchange, a growing list of groups studying FZB42 is available, and many possibilities for interactive data exchange and feedback with the users are given.

40AmyloWiki is configured to be a comprehensive and user-friendly database, built upon typical XAMPP (X-Windows, Linux or Mac OS + Apache + MySQL + PHP + Perl) environment. Apache 2.4.23 was used to construct a webserver. All data sets were processed and stored in MySQL (5.0.11). PHP language (version 5.6.28) was used to built database management system and interface. Webpages were designed with HTML5, CSS3 and JavaScript techniques. AmyloWiki provides a series of functions such as data submission, resource downloading, searching, advanced retrieval, and feedback.

41Briefly, most information of user’s interest can be returned by performing a searching. User can search with different of the query strings, such as gene name, gene locus, and PubMed ID. The items that matched the query string will be returned in the result page. This can be exemplified by searching a gene, as happens most often. The basic information of the gene such as its product, locus, synonyms, homolog in B. subtilis, position, length and others, will be provided on the top of the result page. The genomic context of the gene can be viewed in a visualized window with scrollable function to check its neighbor genes. The organization of the gene, if it is present in an operon, the functions the gene involved, and its functional categories/subcategories are offered next. Other associated information includes the phenotypes of the mutant, its transcriptional start sites, protein/non-coding RNA regulators, sigma factors, PTM sites and so on. The references concerning the gene are listed at the bottom of the retrieval page.

42For the convenience of the user, all datasets of AmyloWiki can be downloaded at the “Download” page. The data can be downloaded in an Excel-compatible format for their specific analysis. AmyloWiki will be maintained by us with a frequent update to improve its configuration and to keep the information comprehensive. For example, it is planned to add in future experimental protocols specifically worked out and used for FZB42, like transformation and bioassay. Here, support given by experienced groups dealing with FZB42 is highly welcomed. The pages for data submission and correction are designed for authorized users in order to update relevant information. Unauthorized users are encouraged to submit their latest data via E-mail to the authors of the website. Then their information will be verified and included in AmyloWiki.

43In order to improve consistency in performance of bioinoculants in a sustainable agriculture we have to integrate them as part of modern crop management programs allowing to decrease the amount of agrochemicals, including harmful chemopesticides. A full understanding of the complex relationship between plant, soil, climate, microbiota, and the microbial inoculant is a necessary precondition for application success of biologicals. Basic research which has been restricted in past to selected representatives of taxonomic groups (‘model organisms’) such as B. subtilis 168, E. coli K12, Saccharomyces cerevisiae, and Drosophila melanogaster has considerably deepen our understanding of those groups in general. We recommend using FZB42 as a model for research on Gram-positive rhizobacteria. This will greatly enhance scientific progress in the field and might contribute to a better consistency in application of environmental friendly beneficial Bacilli in modern agriculture. After 20 years of basic and applied research FZB42 has been proven as suitable for selecting to this task. The following features favor use of FZB42 as model organism:

44Apathogenicity: Concerning biosafety issues, no representatives of the B. subtilis species complex including B. velezensis have been listed as risk group in ‘The Approved List of biological agents’. However, B. cereus and B. anthracis were listed in human pathogen hazard group 3, excluding their use as biocontrol agents in agriculture.

45long term successful application of a commercialized product based on FZB42 in agriculture

46genetic amenable and a large collection of defined mutant is available for interested scientists

47large body of scientific knowledge (>140 articles) about FZB42 is already available

48an integrative data base ‘AmyloWiki’ has been established about FZB42 aimed to enhance collaboration of several groups dealing with Gram-positive PGPR and biocontrol bacteria

49Most of the biocontrol agents currently in use are based on living microbes. Representatives of the B. subtilis species complex, including B. velezensis, B. subtilis, and B. pumilus are increasingly used for commercial production of biofungicides. Most of them are stabilized liquid suspensions or dried formulations prepared from durable endo-spores. They are developed for seed coating, soil or leave application. Unfortunately, it is very unlikely that concentration of Bacillus synthesized cyclic lipopeptides in their natural environment is sufficient for antibiosis. A possibility for circumventing this problem are combined bioformulations consisting of both, Bacillus spores and antagonistic acting metabolites. However, only a few bioformulations currently on the market, such as SERENADE® prepared from B. subtilis QST713 and Double Nickel 55 prepared from B. amyloliquefaciens D747 (both strain names need to be corrected as B. velezensis, ), contain together with living spores antimicrobial compounds, such as cyclic lipopeptides (iturins, fengycin). Unfortunately, also in these products only the number of spores is declared as active ingredient of the biofungicide, but concentration of the metabolites is not indicated, excluding an exact treatment of pathogen infected plant parts. Labeling a fixed concentration of the active principle for suppressing the target would allow a better comparison of chemical and biological pesticides. To the best of our knowledge, no bioformulations containing exclusively antimicrobial metabolites are commercially available, although companies like ABiTEP performed extended large scale trials with concentrated and stabilized Bacillus supernatants in order to suppress plant pathogens.

50BF and RB outlined and wrote the manuscript. BF, CW, XLD, XFS, and RB developed the integrated data base AmyloWiki. LW, HW, and XG contributed essential scientific results reported in this review. All authors have approved and corrected the final version of the manuscript.

51The 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.