Complete Genome Sequence of Paenibacillus polymyxa SQR-21, a Plant Growth-Promoting Rhizobacterium with Antifungal Activity and Rhizosphere Colonization Ability

Genome sequencing and analysis of Bacillus velezensis VJH504 reveal biocontrol mechanism against cucumber Fusarium wilt

One major issue in reducing cucumber yield is the destructive disease Cucumber (Cucumis sativus L.) wilt disease caused by Fusarium oxysporum f. sp. cucumerinum (Foc). When using the isolate VJH504 isolated from cucumber rhizosphere soil and identified as Bacillus velezensis, the growth of Foc in the double culture experiment was effectively inhibited. Phenotypic, phylogenetic, and genomic analyses were conducted to enhance understanding of its biocontrol mechanism. According to the result of the phenotype analysis, B. velezensis VJH504 could inhibit cucumber fusarium wilt disease both in vitro and in vivo, and significantly promote cucumber seed germination and seedling growth. Additionally, the tests of growth-promoting and biocontrol characteristics revealed the secretion of proteases, amylases, β-1,3-glucanases, cellulases, as well as siderophores and indole-3-acetic acid by B. velezensis VJH504. Using the PacBio Sequel II system, we applied the complete genome sequencing for B. velezensis VJH504 and obtained a single circular chromosome with a size of 3.79 Mb. A phylogenetic tree was constructed based on the 16S rRNA gene sequences of B. velezensis VJH504 and 13 other Bacillus species, and Average Nucleotide Identity (ANI) analysis was performed using their whole-genome sequences, confirming isolateVJH504 as B. velezensis. Following this, based on the complete genome sequence od B. velezensis VJH504, specific functional analysis, Carbohydrate-Active Enzymes (CAZymes) analysis, and secondary metabolite analysis were carried out, predicting organism's abilities for biofilm formation, production of antifungal CAZymes, and synthesis of antagonistic secondary metabolites against pathogens. Afterwards, a comparative genomic analysis was performed between B. velezensis VJH504 and three other B. velezensis strains, revealing subtle differences in their genomic sequences and suggesting the potential for the discovery of novel antimicrobial substances in B. velezensis VJH504. In conclusion, the mechanism of B. velezensis VJH504 in controlling cucumber fusarium wilt was predicted to appear that B. velezensis VJH504is a promising biocontrol agent, showcasing excellent application potential in agricultural production.

Biodegradation of PAHs by Bacillus marsiflavi, genome analysis and its plant growth promoting potential

The biodegradation of hazardous petroleum hydrocarbons has recently received a lot of attention because of its many possible applications. Bacillus marsiflavi strain was isolated from oil contaminated soil of Rawalpindi, Pakistan. Initial sequencing was done by 16s rRNA sequencing technique. Bac 144 had shown 78% emulsification index and 72% hydrophobicity content. Further, the strain displayed production of 15.5 mg/L phosphate sloubilization and 30.25 μg/ml indole acetic acid (IAA) in vitro assay. The strain showed 65% biodegradation of crude oil within 5 days by using Gas Chromatography-Mass Spectrometry (GC-MS) analysis. Whole Genome analysis of Bac 144 was performed by PacBio sequencing and results indicated that Bacillus marsiflavi Bac144 strain consisted of size of 4,417,505bp with closest neighbor Bacillus cereus ATCC 14579. The number of the coding sequence was 4662 and number of RNAs was 141. The GC content comprised 48.1%. Various genes were detected in genome responsible for hydrocarbon degradation and plant defense mechanism. The toxic effect of petroleum hydrocarbons in soil and its mitigation with Bac 144 was tested by soil experiment with three levels of oil contamination (5%, 10% and 15%). Soil enzymatic activity such as dehydrogenase and fluorescein diacetate (FDA) increased up to 49% and 40% with inoculation of Bac 144, which was considered to be correlated with hydrocarbon degradation recorded as 46%. An increase of 20%, 14% and 9% in shoot length of plant at 5%, 10% and 15% level of oil was recorded treated with Bac 144 as compared to untreated plants. A percent increase of 14.89%, 16.85%, and 13.87% in chlorophyll, carotenoid, and proline content of plant was observed by inoculation with Bac 144 under oil stress. Significant reduction of 14% and 18%, 21% was recorded in the malondialdehyde content of plant due to inoculation of Bac 144. A considerable increase of 21.33%, 19.5%, and 24.5% in super oxide dismutase, catalase, and peroxidase dismutase activity was also observed in plants inoculated with strain Bac 144. These findings suggested that Bac-144 can be considered as efficient candidate for bioremediation of hydrocarbons.

Genomics assisted functional characterization of Paenibacillus polymyxa HK4 as a biocontrol and plant growth promoting bacterium

The diseases caused by phytopathogens account for huge economic losses in the agricultural sector. Paenibacillus polymyxa is one of the agriculturally important biocontrol agents and plant growth promoting bacterium. This study describes the antifungal potential of P. polymyxa HK4 against an array of fungal phytopathogens and its ability to stimulate seed germination of cumin and groundnut under in vitro conditions. The cumin and groundnut seeds bacterized with HK4 exhibited enhanced germination efficiency in comparison to controls. The use of HK4 as a soil inoculant significantly promoted the shoot length and fresh weight of groundnut plants in pot studies. The draft genome analysis of HK4 revealed the genetic attributes for motility, root colonization, antagonism, phosphate solubilization, siderophore production and production of volatile organic compounds. The bacterium HK4 harnessed several hydrolytic enzymes that may assist its competence in the rhizosphere. The PCR amplification and sequence analysis of the conserved region of the fusA gene amplicon revealed the ability of HK4 to produce fusaricidin. Furthermore, the LC-ESI-MS/MS of crude cell pellet extract of HK4 confirmed the presence of fusaricidin as a major antifungal metabolite. This study demonstrated the potential of HK4 as a biocontrol agent and a plant growth promoter.

Comparative and Functional Analyses of Two Sequenced Paenibacillus polymyxa Genomes Provides Insights Into Their Potential Genes Related to Plant Growth-Promoting Features and Biocontrol Mechanisms

Many bacteria belonging to Paenibacillus polymyxa are plant growth-promoting rhizobacteria (PGPR) with the potential to promote plant growth and suppress phytopathogens and have been used as biological control agents (BCAs). However, the growth promotion and biocontrol mechanisms of P. polymyxa have not been thoroughly elucidated thus far. In this investigation, the genome sequences of two P. polymyxa strains, ZF129 and ZF197, with broad anti-pathogen activities and potential for growth promotion were comparatively studied. Comparative and functional analyses of the two sequenced P. polymyxa genomes showed that the ZF129 genome consists of one 5,703,931 bp circular chromosome and two 79,020 bp and 37,602 bp plasmids, designated pAP1 and pAP2, respectively. The complete genome sequence of ZF197 consists of one 5,507,169 bp circular chromosome and one 32,065 bp plasmid, designated pAP197. Phylogenetic analysis revealed that ZF129 is highly similar to two P. polymyxa strains, HY96-2 and SQR-21, while ZF197 is highly similar to P. polymyxa strain J. The genes responsible for secondary metabolite synthesis, plant growth-promoting traits, and systemic resistance inducer production were compared between strains ZF129 and ZF197 as well as other P. polymyxa strains. The results indicated that the variation of the corresponding genes or gene clusters between strains ZF129 and ZF197 may lead to different antagonistic activities of their volatiles or cell-free supernatants against Fusarium oxysporum. This work indicates that plant growth promotion by P. polymyxa is largely mediated by phytohormone production, increased nutrient availability and biocontrol mechanisms. This study provides an in-depth understanding of the genome architecture of P. polymyxa, revealing great potential for the application of this bacterium in the fields of agriculture and horticulture as a PGPR.

Plant growth promoting rhizobacteria isolated from halophytes and drought-tolerant plants: genomic characterisation and exploration of phyto-beneficial traits

Plant growth promoting rhizobacteria (PGPR) are able to provide cross-protection against multiple stress factors and facilitate growth of their plant symbionts in many ways. The aim of this study was to isolate and characterize rhizobacterial strains under natural conditions, associated with naturally occurring representatives of wild plant species and a local tomato cultivar, growing in differently stressed Mediterranean ecosystems. A total of 85 morphologically different rhizospheric strains were isolated; twenty-five exhibited multiple in vitro PGP-associated traits, including phosphate solubilization, indole-3-acetic acid production, and 1-aminocyclopropane-1-carboxylate deaminase activity. Whole genome analysis was applied to eight selected strains for their PGP potential and assigned seven strains to Gammaproteobacteria, and one to Bacteroidetes. The genomes harboured numerous genes involved in plant growth promotion and stress regulation. They also support the notion that the presence of gene clusters with potential PGP functions is affirmative but not necessary for a strain to promote plant growth under abiotic stress conditions. The selected strains were further tested for their ability to stimulate growth under stress. This initial screening led to the identification of some strains as potential PGPR for increasing crop production in a sustainable manner.

Pan-genome analysis of Paenibacillus polymyxa strains reveals the mechanism of plant growth promotion and biocontrol

Rapid development of gene sequencing technologies has led to an exponential increase in microbial sequencing data. Genome research of a single organism does not capture the changes in the characteristics of genetic information within a species. Pan-genome analysis gives us a broader perspective to study the complete genetic information of a species. Paenibacillus polymyxa is a Gram-positive bacterium and an important plant growth-promoting rhizobacterium with the ability to produce multiple antibiotics, such as fusaricidin, lantibiotic, paenilan, and polymyxin. Our study explores the pan-genome of 14 representative P. polymyxa strains isolated from around the world. Heap's law model and curve fitting confirmed an open pan-genome of P. polymyxa. The phylogenetic and collinearity analyses reflected that the evolutionary classification of P. polymyxa strains are not associated with geographical area and ecological niches. Few genes related to phytohormone synthesis and phosphate solubilization were conserved; however, the nif cluster gene associated with nitrogen fixation exists only in some strains. This finding is indicative of nitrogen fixing ability is not stable in P. polymyxa. Analysis of antibiotic gene clusters in P. polymyxa revealed the presence of these genes in both core and accessory genomes. This observation indicates that the difference in living environment led to loss of ability to synthesize antibiotics in some strains. The current pan-genomic analysis of P. polymyxa will help us understand the mechanisms of biological control and plant growth promotion. It will also promote the use of P. polymyxa in agriculture.

Draft Genome Sequences of Bacillus and Paenibacillus Species Isolated from Seeds of Citrullus lanata (Watermelon), Cucurbita moschata (Butternut Squash), and Cucurbita pepo L. var. pepo L. (Pumpkin)

Here, we announce the draft genome sequences of four endophytic bacilli isolated from surface-sterilized seeds of three cucurbit species, Bacillus sp. strains EKM417B and EKM420B (from Citrullus lanata [watermelon]) and EKM501B (from Cucurbita moschata [butternut squash]) and Paenibacillus sp. strain EKM301P (from Cucurbita pepo L. var. pepo L. [pumpkin]). These strains previously demonstrated biostimulant and biocontrol activities.

Draft Genome Sequence of Bacillus sp. Strain EKM601B (Phylum Firmicutes), Living inside the Seeds of Luffa acutangula (Chinese Okra)

Presented here is the draft genome sequence of Bacillus sp. strain EKM601B, which contains 4,199,360 bp in 73 contigs. This candidate endophyte was isolated from surface-sterilized dry seeds of Luffa acutangula (Chinese okra) and demonstrated diverse plant-beneficial functions and antagonism against soilborne pathogens in vitro.

Presented here is the draft genome sequence of Bacillus sp. strain EKM601B, which contains 4,199,360 bp in 73 contigs. This candidate endophyte was isolated from surface-sterilized dry seeds of Luffa acutangula (Chinese okra) and demonstrated diverse plant-beneficial functions and antagonism against soilborne pathogens in vitro.

Complete Genome Sequence of Industrial Biocontrol Strain Paenibacillus polymyxa HY96-2 and Further Analysis of Its Biocontrol Mechanism

Paenibacillus polymyxa (formerly known as Bacillus polymyxa) has been extensively studied for agricultural applications as a plant-growth-promoting rhizobacterium and is also an important biocontrol agent. Our team has developed the P. polymyxa strain HY96-2 from the tomato rhizosphere as the first microbial biopesticide based on P. polymyxa for controlling plant diseases around the world, leading to the commercialization of this microbial biopesticide in China. However, further research is essential for understanding its precise biocontrol mechanisms. In this paper, we report the complete genome sequence of HY96-2 and the results of a comparative genomic analysis between different P. polymyxa strains. The complete genome size of HY96-2 was found to be 5.75 Mb and 5207 coding sequences were predicted. HY96-2 was compared with seven other P. polymyxa strains for which complete genome sequences have been published, using phylogenetic tree, pan-genome, and nucleic acid co-linearity analysis. In addition, the genes and gene clusters involved in biofilm formation, antibiotic synthesis, and systemic resistance inducer production were compared between strain HY96-2 and two other strains, namely, SC2 and E681. The results revealed that all three of the P. polymyxa strains have the ability to control plant diseases via the mechanisms of colonization (biofilm formation), antagonism (antibiotic production), and induced resistance (systemic resistance inducer production). However, the variation of the corresponding genes or gene clusters between the three strains may lead to different antimicrobial spectra and biocontrol efficacies. Two possible pathways of biofilm formation in P. polymyxa were reported for the first time after searching the KEGG database. This study provides a scientific basis for the further optimization of the field applications and quality standards of industrial microbial biopesticides based on HY96-2. It may also serve as a reference for studying the differences in antimicrobial spectra and biocontrol capability between different biocontrol agents.

Characterization and antifungal activity against Pestalotiopsis of a fusaricidin-type compound produced by Paenibacillus polymyxa Y-1

Dendrobium nobile (D. nobile) is a valuable Chinese herbal medicine. The discovery of microbial resources from has provided a wealth of raw materials. Stalk rot, which is caused by Pestalotiopsis, is one of the most serious diseases of D nobile and has resulted in serious losses in production. However, an effective method for the prevention and control of stalk rot remains lacking. In this study, we aimed to identify a biocontrol strain against Pestalotiopsis. We isolated Paenibacillus polymyxa Y-1, an endophytic bacterium, from the stem of D. nobile. Three pairs of active metabolites isolated from this bacterium were identified as fusaricidin compounds. We then investigated the mechanism of fusaricidin compounds on Pestalotiopsis via proteomics. Proteomics data showed that the compounds mainly inhibit energy generation in the respiratory chain and amino acid biosynthesis of Pestalotiopsis.

Comprehensive genomic analysis of a plant growth-promoting rhizobacterium Pantoea agglomerans strain P5

In this study, we provide a comparative genomic analysis of Pantoea agglomerans strain P5 and 10 closely related strains based on phylogenetic analyses. A next-generation shotgun strategy was implemented using the Illumina HiSeq 2500 technology followed by core- and pan-genome analysis. The genome of P. agglomerans strain P5 contains an assembly size of 5082485 bp with 55.4% G + C content. P. agglomerans consists of 2981 core and 3159 accessory genes for Coding DNA Sequences (CDSs) based on the pan-genome analysis. Strain P5 can be grouped closely with strains PG734 and 299 R using pan and core genes, respectively. All the predicted and annotated gene sequences were allocated to KEGG pathways. Accordingly,  genes involved in plant growth-promoting (PGP) ability, including phosphate solubilization, IAA and siderophore production, acetoin and 2,3-butanediol synthesis and bacterial secretion, were assigned. This study provides an in-depth view of the PGP characteristics of strain P5, highlighting its potential use in agriculture as a biofertilizer.

Whole genome analysis of halotolerant and alkalotolerant plant growth-promoting rhizobacterium Klebsiella sp. D5A

This research undertook the systematic analysis of the Klebsiella sp. D5A genome and identification of genes that contribute to plant growth-promoting (PGP) traits, especially genes related to salt tolerance and wide pH adaptability. The genome sequence of isolate D5A was obtained using an Illumina HiSeq 2000 sequencing system with average coverages of 174.7× and 200.1× using the paired-end and mate-pair sequencing, respectively. Predicted and annotated gene sequences were analyzed for similarity with the Kyoto Encyclopedia of Genes and Genomes (KEGG) enzyme database followed by assignment of each gene into the KEGG pathway charts. The results show that the Klebsiella sp. D5A genome has a total of 5,540,009 bp with 57.15% G + C content. PGP conferring genes such as indole-3-acetic acid (IAA) biosynthesis, phosphate solubilization, siderophore production, acetoin and 2,3-butanediol synthesis, and N2 fixation were determined. Moreover, genes putatively responsible for resistance to high salinity including glycine-betaine synthesis, trehalose synthesis and a number of osmoregulation receptors and transport systems were also observed in the D5A genome together with numerous genes that contribute to pH homeostasis. These genes reveal the genetic adaptation of D5A to versatile environmental conditions and the effectiveness of the isolate to serve as a plant growth stimulator.

Characterization of Novel Fusaricidins Produced by Paenibacillus polymyxa-M1 Using MALDI-TOF Mass Spectrometry

Paenibacillus polymyxa-M1 is a potent producer of bioactive compounds, such as lipopeptides, polyketides, and lantibiotics of biotechnological and medical interest. Genome sequencing revealed nine gene clusters for nonribosomal biosynthesis of such agents. Here we report on the investigation of the fusaricidins, a complex of cyclic lipopeptides containing 15-guanidino-3-hydroxypentadecanoic acid (GHPD) as fatty acid component by matrix-assisted laser-desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). More than 20 variants of these compounds were detected and characterized in detail. Mass spectrometric sequence analysis was performed by MALDI-LIFT-TOF/TOF fragment analysis. The obtained product ion spectra show a specific processing in the fatty acid part. GHPD is cleaved between the α- and ß-position yielding two fragments a and b, one bearing the end-standing guanidine group and another one comprising the residual two C-atoms of GHPD with the attached peptide moiety. The complete sequence of all fusaricidins was derived from sets of bn- and yn-ions. The fusaricidin complex can be divided into four lipopeptide families, three of them showing variations of the amino acid in position 3, Val or Ile for the first and Tyr or Phe for families 2 and 3, respectively. A collection of novel fusaricidins was detected differing from those of families 1-3 by an additional residue of 71 Da (family 4). LIFT-TOF/TOF fragment spectra of these species imply that in their peptide moiety, an Ala-residue is attached by an ester bond to the free hydroxyl group of Thr4. More than 10 novel fusaricidins were characterized mass spectrometrically.

Draft Genome Sequence of Paenibacillus polymyxa EBL06, a Plant Growth-Promoting Bacterium Isolated from Wheat Phyllosphere

Paenibacillus polymyxa strain EBL06 is a plant growth-promoting bacterium with high antifungal activity. The estimated genome of this strain is 5.68 Mb in size and harbors 4,792 coding sequences (CDSs).

Development and validation of an rDNA operon based primer walking strategy applicable to de novo bacterial genome finishing

Advances in sequencing technology have drastically increased the depth and feasibility of bacterial genome sequencing. However, little information is available that details the specific techniques and procedures employed during genome sequencing despite the large numbers of published genomes. Shotgun approaches employed by second-generation sequencing platforms has necessitated the development of robust bioinformatics tools for in silico assembly, and complete assembly is limited by the presence of repetitive DNA sequences and multi-copy operons. Typically, re-sequencing with multiple platforms and laborious, targeted Sanger sequencing are employed to finish a draft bacterial genome. Here we describe a novel strategy based on the identification and targeted sequencing of repetitive rDNA operons to expedite bacterial genome assembly and finishing. Our strategy was validated by finishing the genome of Paenibacillus polymyxa strain CR1, a bacterium with potential in sustainable agriculture and bio-based processes. An analysis of the 38 contigs contained in the P. polymyxa strain CR1 draft genome revealed 12 repetitive rDNA operons with varied intragenic and flanking regions of variable length, unanimously located at contig boundaries and within contig gaps. These highly similar but not identical rDNA operons were experimentally verified and sequenced simultaneously with multiple, specially designed primer sets. This approach also identified and corrected significant sequence rearrangement generated during the initial in silico assembly of sequencing reads. Our approach reduces the required effort associated with blind primer walking for contig assembly, increasing both the speed and feasibility of genome finishing. Our study further reinforces the notion that repetitive DNA elements are major limiting factors for genome finishing. Moreover, we provided a step-by-step workflow for genome finishing, which may guide future bacterial genome finishing projects.