Genome sequencing of biocontrol strain Bacillus amyloliquefaciens Bam1 and further analysis of its heavy metal resistance mechanism

Abstract

Plant growth-promoting rhizobacteria (PGPR) or Biocontrol strains inevitably encounter heavy metal excess stress during the product’s processing and application. Bacillus amyloliquefaciens Bam1 was a potential biocontrol strain with strong heavy metal resistant ability. To understand its heavy metal resistance mechanism, the complete genome of Bam1 had been sequenced, and the comparative genomic analysis of Bam1 and FZB42, an industrialized PGPR and biocontrol strain with relatively lower heavy metal tolerance, was conducted. The comparative genomic analysis of Bam1 and the other nine B. amyloliquefaciens strains as well as one Bacillus velezensis (genetically and physiologically very close to B. amyloliquefaciens) was also performed. Our results showed that the complete genome size of Bam1 was 3.95 Mb, 4219 coding sequences were predicted, and it possessed the highest number of unique genes among the eleven analyzed strains. Nine genes related to heavy metal resistance were detected within the twelve DNA islands of Bam1, while only two of them were detected within the seventeen DNA islands of FZB42. When compared with B. amyloliquefaciens type strain DSM7, Bam1 lacked contig L, whereas FZB42 lacked contig D and I, as well as just possessed contig B with a very small size. Our results could also deduce that Bam1 promoted its essential heavy metal resistance mainly by decreasing the import and increasing the export of heavy metals with the corresponding homeostasis systems, which are regulated by different metalloregulators. While Bam1 promoted its non-essential heavy metal resistance mainly by the activation of some specific or non-specific exporters responding to different heavy metals. The variation of the genes related to heavy metal resistance and the other differences of the genomes, including the different number and arrangement of contigs, as well as the number of the heavy metal resistant genes in Prophages and Genomic islands, led to the significant different resistance of Bam1 and FZB42 to heavy metals.

Keywords

Bacillus amyloliquefaciens / Heavy metal resistance mechanism / Comparative genomic analysis / Resistance to essential heavy metals / Resistance to non-essential heavy metals

References

[1] Aguilar-Barajas E, Jacobo-Arreola S, Jiménez-Mejıa R, Riveros-Rosas H, Verduzco-Rosas LA, Ramırez-Dıaz MI, Julian-Sanchez A, Cervantes C. An Lrp-type transcriptional regulator controls expression of the Bacillus subtilis chromate transporter. Antonie Van Leeuwenhoek, 2013, 104: 941-948.

[2] Alloway BJ. Micronutrient deficiencies in global crop production: an introduction, 2008, Dordrecht: Springer

[3] Arndt D, Marcu A, Liang Y, Wishart DS. PHAST, PHASTER and PHASTEST: tools for finding prophage in bacterial genomes. Brief Bioinform, 2019, 20(4): 1560-1567.

[4] Banci L, Rosato A. Structural genomics of proteins involved in copper homeostasis. Acc Chem Res, 2003, 36: 215-221.

[5] Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol, 2012, 19(5): 455-477.

[6] Bertelli C, Laird MR, Williams KP, Simon F University Research Computing Group4 Lau BY, Hoad G, Winsor GL, Brinkman FSL. IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res, 2017, 45: W30-W35.

[7] Borriss R, Chen XH, Rueckert C, Blom J, Becker A, Baumgarth BG, Fan B, Pukall R, Schumann P, Spröer C, Junge H, Vater J, Pühler A, Klenk HP. Relationship of Bacillus amyloliquefaciens clades associated with strains DSM 7T and FZB42T: a proposal for Bacillus amyloliquefaciens subsp. amyloliquefaciens subsp. nov. and Bacillus amyloliquefaciens subsp. plantarum subsp. nov. based on complete genome sequence comparisons. Int J Syst Evol Microbiol, 2011, 61(Pt 8): 1786-1801.

[8] Cai S, Deng K, Li J, Xu R. Effects of metal ions on phosphatase activity of periphytic biofilm in paddy fields. Soils, 2020, 52(3): 525-531. Chinese Journal

[9] Carrondo MA. Ferritins, iron uptake and storage from the bacterioferritin viewpoint. EMBO J, 2003, 22(9): 1959-1968.

[10] Chan PP, Lowe TM. tRNAscan-SE: searching for tRNA genes in genomic sequences. Methods Mol Biol, 2019, 1962: 1-14.

[11] Chandrangsu P, Helmann JD. Intracellular Zn(II) intoxication leads to dysregulation of the PerR regulon resulting in heme toxicity in Bacillus subtilis. PLoS Genet, 2016, 12(12

[12] Chandrangsu P, Rensing C, Helmann DJ. Metal homeostasis and resistance in bacteria. Nat Rev Microbiol, 2017, 15(6): 338-350.

[13] Chen F, Mackey AJ, Stoeckert CJ, Roos DS. OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res, 2006, 34: D363-D368.

[14] Chen XH, Koumoutsi A, Scholz R, Eisenreich A, Schneider K, Heinemeyer I, Morgenstern B, Voss B, Hess WR, Reva O, Junge H, Voigt B, Jungblut PR, Vater J, Sussmuth R, Liesegang H, Strittmatter A, Gottschalk G, Borriss R. Comparative analysis of the complete genome sequence of the plant growth-promoting bacterium Bacillus amyloliquefaciens FZB42. Nat Biotechnol, 2007, 25: 1007-1014.

[15] Chillappagari S, Miethke M, Trip H, Kuipers OP, Marahiel MA. Copper acquisition is mediated by YcnJ and regulated by YcnK and CsoR in Bacillus subtilis. J Bacteriol, 2009, 191(7): 2362-2370.

[16] Choi J, Nam J, Seo MH. Complete genome sequence of Bacillus velezensis NST6 and comparison with the species belonging to operational group B. amyloliquefaciens. Genomics, 2020, 113(1 Pt 1): 380-386.

[17] Chun BH, Kim KH, Jeong SE, Jeon CO. Genomic and metabolic features of the Bacillus amyloliquefaciens group-B. amyloliquefaciens, B. velezensis, and B. siamensis-revealed by pan-genome analysis. Food Microbiol, 2019, 77: 146-157.

[18] Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res, 1999, 27: 4636-4641.

[19] Didovyk A, Borek B, Tsimring L, Hasty J. Transcriptional regulation with CRISPR-Cas9: principles, advances, and applications. Curr Opin Biotechnol, 2016, 40: 177-184.

[20] Fan B, Wang C, Song XF, Ding XL, Wu LM, Wu HJ, Gao XW, Borriss R. Bacillus velezensis FZB42 in 2018: the Gram-positive model strain for plant growth promotion and biocontrol. Front Microbiol, 2018, 9: 2491.

[21] Fouts DE. Phage_Finder: automated identification and classification of prophage regions in complete bacterial genome sequences. Nucleic Acids Res, 2006, 34(20): 5839-5851.

[22] Gaballa A, Helmann JD. Bacillus subtilis CPx-type ATPases: characterization of Cd, Zn, Co and Cu efflux systems. Biometals, 2003, 16: 497-505.

[23] Gaballa A, MacLellan S, Helmann JD. Transcription activation by the siderophore sensor Btr is mediated by ligand-dependent stimulation of promoter clearance. Nucleic Acids Res, 2012, 40: 3585-3595.

[24] Gautam S, Chauhan A, Sharma R, Sehgal R, Shirkot CK. Potential of Bacillus amyloliquefaciens for biocontrol of bacterial canker of tomato incited by Clavibacter michiganensis ssp. michiganensis. Microb Pathog, 2019, 130: 196-203.

[25] Gendy S, Chauhan A, Agarwal MA, Rathore PRS, Rajneesh J. Is long-term heavy metal exposure driving carriage of antibiotic resistance in environmental opportunistic pathogens: a comprehensive phenomic and genomic assessment using Serratia sp. SRS-8-S-2018. Front Microbiol, 2020, 11: 1923.

[26] Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res, 2007, 35: W52-W57.

[27] Guan G, Pinochet-Barros A, Gaballa A, Patel SJ, Argüello JM, Helmann JD. PfeT, a P1B4 -type ATPase, effluxes ferrous iron and protects Bacillus subtilis against iron intoxication. Mol Microbiol, 2015, 98: 787-803.

[28] Hamlett NV, Landale EC, Davids BH, Summer AO. Roles of the Tn 21 merT, merP, and merC gene produce in mercury resistance and mercury binding. J Bacterial, 1992, 174: 6377-6385.

[29] Harvie DR, Andreini C, Cavallaro G, Meng W, Connolly BA, Yoshida K, Fujita Y, Harwood CR, Radford DS, Tottey S, Cavet JS, Robinson NJ. Predicting metals sensed by ArsR-SmtB repressors: allosteric interference by a non-effector metal. Mol Microbiol, 2006, 59(4): 1341-1356.

[30] He Y, Dong L, Zhou S, Jia Y, Gu R, Bai Q, Gao J, Li Y, Xiao H. Chromium resistance characteristics of Cr(VI) resistance genes ChrA and ChrB in Serratia sp. S2. Ecotox Environ Safe, 2018, 157: 417-423.

[31] Hooda PS. Trace elements in soil, 2010, Hoboken: Blackwell Publishing Ltd

[32] Huang XJ, Shin JH, Azul PB, Su TT, Helmann JD. Bacillus subtilis MntR coordinates the transcriptional regulation of manganese uptake and efflux systems. Mol Microbiol, 2017, 103(2): 253-268.

[33] Huang YY, He CT, Chuang S, Guo JJ, Mubeen S, Yuan JG, Yang ZY. Toxicity of cadmium and its health risks from leafy vegetable consumption. Food Funct, 2017, 8(4): 1373-1401.

[34] Islam MA, Romić D, Akber MA, Romic M. Trace metals accumulation in soil irrigated with polluted water and assessment of human health risk from vegetable consumption in Bangladesh. Environ Geochem Health, 2018, 40: 59-85.

[35] Jarosławiecka A, Zofia PS. Lead resistance in micro-organisms. Microbiology, 2014, 160: 12-25.

[36] Kim EJ, Chung HJ, Suh B, Hah YC, Roe JH. Expression and regulation of the sodFgene encoding iron and zinc-containing superoxide dismutase in Streptomyces coelicolor Müller. J Bacteriol, 1998, 180(8): 2014-2020.

[37] Kröber M, Verwaaijen B, Wibberg D, Winkler A, Pühler A, Schlüter A. Bacillus amyloliquefaciens FZB42 as response to biofilm formation analyzed by RNA sequencing. J Biotechnol, 2016, 231: 212-223.

[38] Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res, 2009, 19: 1639-1645.

[39] Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL. Versatile and open software for comparing large genomes. Genome Biol, 2004, 5(2): R12.

[40] Li Z, Liu BF, Cui H, Ding J, Li HH, Xie GJ, Ren NQ, Xing DF. The complete genome sequence of Ethanoligenens harbinense reveals the metabolic pathway of acetate–ethanol fermentation: a novel understanding of the principles of anaerobic biotechnology. Environ Int, 2019, 131.

[41] Liu N, Zhu L, Zhang Z, Huang H, Jiang L. Draft genome sequence of a multidrug-resistant blaOXA-69-producing Acinetobacter baumannii L13 isolated from Tarim River sample in China. J Glob Antimicrob Resist, 2019, 18: 145-147.

[42] Lu BX, Leong HW. Computational methods for predicting genomic islands in microbial genomes. Comput Struct Biotechnol J, 2016, 14: 200-206.

[43] Luo YC, Cheng YJ, Yi JC, Zhang ZJ, Luo Q, Zhang DJ, Li YG. Complete genome sequence of industrial biocontrol strain Paenibacillus polymyxa HY96-2 and further analysis of its biocontrol mechanism. Front Microbiol, 2018, 9: 1520.

[44] Ma Z, Faulkner MJ, Helmann JD. Origins of specificity and cross-talk in metal ion sensing by Bacillus subtilis Fur. Mol Microbiol, 2012, 86: 1144-1155.

[45] Matsui K, Yoshinami S, Narita M, Chien MF, Phung LT, Silver S, Endo G. Mercury resistance transposons in Bacilli strains from different geographical regions. FEMS Microbiol Lett, 2016, 363(5

[46] Mitra B, Sharma R. The cysteine-rich amino-terminal domain of ZntA, a Pb(II)/Zn(II)/Cd(II)-translocating ATPase from Escherichia coli, is not essential for its function. Biochemistry, 2011, 40(25): 7694-7699.

[47] Moore CM, Helmann JD. Metal ion homeostasis in Bacillus subtilis. Curr Opin Microbiol, 2005, 8(2): 188-195.

[48] Moore CM, Gaballa A, Hui M, Ye RW, Helmann JD. Genetic and physiological responses of Bacillus subtilis to metal ion stress. Mol Microbiol, 2005, 57(1): 27-40.

[49] Mu C, Zhu W. Antibiotic effects on gut microbiota, metabolism, and beyond. Appl Microbiol Biotechnol, 2019, 103(23–24): 9277-9285.

[50] Nguyen LT, Heiko A, Schmidt AH, Bui QM. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol, 2015, 32(1): 268-274.

[51] Nicolaou SA, Fast AG, Nakamaru-Ogiso E, Papoutsakis ET. Overexpression of fetA (ybbL) and fetB (ybbM), encoding an iron exporter, enhances resistance to oxidative stress in Escherichia coli. Appl Environ Microbiol, 2013, 79(23): 7210-7219.

[52] Pi H, Helmann DJ. Sequential induction of Fur-regulated genes in response to iron limitation in Bacillus subtilis. Proc Natl Acad Sci USA, 2017, 114(48): 12785-12790.

[53] Pi H, Wendel BM, Helmann JD. Dysregulation of magnesium transport protects Bacillus subtilis against manganese and cobalt intoxication. J Bacteriol, 2020, 202(7): e00711-19.

[54] Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol, 2008, 25(7): 1253-1256.

[55] Ravi A, Rajan S, Khalid NK, Jose MS, Aravindakumar CT, Krishnankutty RE. Impact of supplements on enhanced activity of Bacillus amyloliquefaciens Bmb1 against Pythium aphanidermatum through lipopeptide modulation. Probiotics Antimicrob Proteins, 2021, 13(2): 367-374.

[56] Sachla AJ, Luo YC, Helmann JD. Manganese impairs the QoxABCD terminal oxidase leading to respiration-associated toxicity. Mol Microbiol, 2021

[57] Silva FJ, Ferreira LC, Campos VP, Cruz-Magalhães V, Barros AF, Andrade JP, Roberts DP, de Souza JT. Complete genome sequence of the biocontrol agent Bacillus velezensis UFLA258 and its comparison with related species: diversity within the commons. Genome Biol Evol, 2019, 11(10): 2818-2823.

[58] Smaldone GT, Revelles O, Gaballa A, Sauer U, Antelmann H, Helmann JD. A global investigation of the Bacillus subtilis iron-sparing response identifies major changes in metabolism. J Bacteriol, 2012, 194(10): 2594-2605.

[59] Sullivan MJ, Goh KGK, Gosling D, Katupitiya L, Ulett GC. Copper intoxication in group B Streptococcus triggers transcriptional activation of the cop operon that contributes to enhanced virulence during acute infection. J Bacteriol, 2021, 203(19

[60] Tsai KJ, Lin YF, Wong MD, Yang HHC, Fu HL, Rosen BP. Membrane topology of the pl258 CadA Cd(II)/Pb(II)/Zn(II)-translocating P-Type ATPase. J Bioenerg Biomembr, 2002, 34: 147-156.

[61] VanderWal AR, Makthal N, Pinochet-Barros A, Helmann JD, Olsen RJ, Kumaraswami M. Iron efflux by PmtA is critical for oxidative stress resistance and contributes significantly to group A Streptococcus virulence. Infect Immun, 2017, 85(6): e00091-17.

[62] Vij S, Tyag AK. A20/AN1 zinc-finger domain-containing proteins in plants and animals represent common elements in stress response. Funct Integr Genom, 2008, 8: 301-307.

[63] Wang M, Chen SB, Zheng H, Li SS, Chen L, Wang D. The responses of cadmium phytotoxicity in rice and the microbial community in contaminated paddy soils for the application of different long-term N fertilizers. Chemosphere, 2020, 238.

[64] Watanabe M, Lee K, Goto K, Kumagai S, Sugita-Konishi Y, HaraKudo Y. Rapid and effective DNA extraction method with bead grinding for a large amount of fungal DNA. J Food Prot, 2010, 73: 1077-1084.

[65] Welker NE, Campbell LL. Unrelatedness of Bacillus amyloliquefaciens and Bacillus subtilis. J Bacteriol, 1967, 94(4): 1124-1130.

[66] Xu J, Wang XY, Guo WZ. The cytochrome P450 superfamily: key players in plant development and defense. J Integr Agric, 2015, 14(9): 1673-1686.

[67] Yazdia M, Kolahib M, Kazemic EM, Barnaby AG. Study of the contamination rate and change in growth features of lettuce (Lactuca sativa Linn.) in response to cadmium and a survey of its phytochelatin synthase gene. Ecotoxicol Environ Saf, 2019, 180: 295-308.

[68] Yi JC, Zhang DJ, Cheng YJ, Tan JJ, Luo YC. The impact of Paenibacillus polymyxa HY96-2 luxS on biofilm formation and control of tomato bacterial wilt. Appl Microbiol Biotechnol, 2019, 103: 9643-9657.

[69] Yuan SP, Sun XH, Ji AM, Deng XF, Zhou JF. Analysts on heavy metal pollution in vegetables of big and middle city in China. J Anhui Agn Sci, 2013, 41(1): 247-248 Chinese Journal

[70] Zeng F, Wu X, Qiu B, Wu F, Jiang L, Zhang G. Physiological and proteomic alterations in rice (Oryza sativa L.) seedlings under hexavalent chromium stress. Planta, 2014, 240(2): 291-308.

[71] Zhai LX, Ren RX, Meng D, Tian QP, Guan ZB, Cai YJ, Liao XG. Comparison of aminotransferases of three Bacillus strains Bacillus altitudinis W3, Bacillus velezensis SYBC H47, and Bacillus amyloliquefaciens YP6 via genome analysis and bioinformatics. J Appl Genet, 2019, 60(3–4): 427-430.

Funding

Shanghai Agriculture Applied Technology Development Program(2021-02-08-00-12-F00771)

Key Research and Development Program of Ningxia Hui Autonomous Region(2020BBF03004)

the Fundamental Research Funds for the Central Universities(22221818014)

the 111 Project(B18022)

Open Funding Project of the State Key Laboratory of Bioreactor Engineering