Bingyu Li1†Junyi Wan1†Jingjing Sha1†Mengyuan Tian1Mengchen Wang1Xinyuan Zhang1Wei Sun1Yanrong Mao1Jie Min1Yiyuan Qin1Yuqing Liu1
Wenhe Wang1
Xiangfeng He1,2,3*- 1Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing University of Agriculture, Beijing, China
- 2College of Landscape Architecture, Beijing University of Agriculture, Beijing, China
- 3Beijing Laboratory of Urban and Rural Ecological Environment, Beijing, China
Lily (Lilium spp.) is one of the most famous ornamental flowers globally. Lily basal rot (also known as root rot or stem rot) and lily gray mold have seriously affected the yield and quality of lily, resulting in huge economic losses. In this study, bacterial strain E was isolated from a continuous lily cropping field. Strain E displayed high control efficiency against lily basal rot and gray mold, caused by Fusarium oxysporum and Botrytis cinerea respectively, and promoted the occurrence of scale bulblets. Strain E displayed strong inhibitory effects against several other plant pathogenic fungi and two pathogenic bacteria in dual culture and disc diffusion assays, respectively. Whole genome sequencing revealed that strain E contained a 3,929,247 bp circular chromosome with 4,056 protein-coding genes and an average GC content of 47.32%. Strain E was classified as Bacillus velezensis using genome-based phylogenetic analysis and average nucleotide identity and digital DNA–DNA hybridization analyses. A total of 86 genes and 13 secondary metabolite biosynthetic gene clusters involved in antifungal and antibacterial activity, plant growth promotion, colonization, nutrient uptake and availability were identified in the genome of strain E. In vitro biochemical assays showed that strain E produced siderophores, proteases, cellulases, biofilms, antifungal and antibacterial substances, and exhibited organic phosphate solubilization and swimming and swarming motility, which were consistent with the results of the genome analysis. Colonization analysis showed that strain E could colonize the root of the lily, but not the leaf. Overall, these results demonstrate that B. velezensis strain E can be used as a potential biofertilizer and biocontrol agent for lily production.
Introduction
Lily (Lilium spp.) is a perennial bulbous flower belonging to the Liliaceae family. Owing to its large and colorful flowers, lily is one of the most famous ornamental flowers worldwide. China is one of the world’s largest producers of lily cut flowers. In recent years, diseases affecting lily have become increasingly serious because of continuous cropping. Lily basal rot (also known as root rot or stem rot) and lily gray mold have seriously affected the yield and quality of lily, and resulting in huge economic losses (Cui et al., 2018; He et al., 2019). Lily basal rot is a soil-borne disease that mainly affects the underground parts of plants and results in wilting, yellowing, and rotting of the bulbs. Fusarium oxysporum f. sp. lilii is the main pathogen that causes lily basal rot (Straathof and Löffler, 1994). Gray mold, caused by the fungus Botrytis cinerea, is a common fungal disease in flowers and one of the most serious diseases in Lilium (de Cáceres et al., 2015). Lily flowers infected by B. cinerea display white, brown, or light brown hygrophanous lesions, and even decay severely (Ma et al., 2018). Currently, chemical control methods are employed to control lily basal rot and lily gray mold. However, chemical control has shown dwindling effects because the pathogens are seed-borne or soil-borne, and pathogens develop resistance to the chemicals used (Harold et al., 1997; He et al., 2019). Simultaneously, the long-term, large-scale and repeated use of chemical pesticides has caused environmental pollution and safety problems in edible and medicinal lily. Therefore, pollution-free biological control is the preferred choice for controlling lily diseases.
Biocontrol has become a promising strategy for controlling various plant diseases (Zhang et al., 2017; Fan et al., 2018; Mao et al., 2020; Ding et al., 2021). To date, several biological control agents (BCAs) have been reported. However, only a limited number of strains can be commercialized, because storing and transporting BCAs as stable products is a difficult process (Fan et al., 2018; Ding et al., 2021). The Gram-positive Bacillus is one of the most important commercial BCAs because it can form long-lived, heat- and desiccation-resistant endospores (de Benito Armas et al., 2008; Fan et al., 2018). Some strains of Bacillus exhibit biological activity against plant pathogens, including bacteria, fungi, and nematodes, by secreting extracellular metabolites and enzymes (Rabbee et al., 2019; Miljakovi et al., 2020). Additionally, Bacillus strains can trigger induced systemic resistance (ISR) in plants to protect them from attacks by various plant pathogens (Fan et al., 2018; Miljakovi et al., 2020). In addition to being biocontrol agents, Bacillus strains can promote plant growth through diverse mechanisms, such as the production of 1-Aminocyclopropane-1-caroxylate deaminase (ACC), siderophores and phytohormones, and processes such as phosphate solubilization and nitrogen fixation (Bhattacharyya and Jha, 2012; Fan et al., 2018).
Bacillus velezensis is one of the most valuable BCAs in the genus Bacillus. B. velezensis was first isolated from environmental samples in 2005 (Ruiz-García et al., 2005). Based on the results of recent phylogenetic analyses, several Bacillus species, including all strains previously classified as B. velezensis, B. methylotrophicus, and B. amyloliquefaciens subsp. plantarum, have been reclassified as B. velezensis (Dunlap et al., 2016; Rabbee et al., 2019).
So far, many strains of B. velezensis have been reported to promote plant growth and control diseases in various plants, including shallot, soybean, lotus, walnut, Quercus, Japanese cypress, strawberry, sesame, apple, ginger, pepper, poplar, tomato, grape, tobacco, watermelon, wheat, maize, and cotton (Cai et al., 2017; Jiang et al., 2019; Aulia Rahma et al., 2020; Bayisa, 2020; Wang et al., 2020; Choub et al., 2021; Ding et al., 2021; Hamaoka et al., 2021; Han et al., 2021; Liang et al., 2021; Liu et al., 2021; Masmoudi et al., 2021; Moon et al., 2021; Sachin et al., 2021; Shin et al., 2021; Won et al., 2021; Adeniji and Babalola, 2022; Hong et al., 2022; Yuan et al., 2022). Although B. velezensis has attracted increasing attention as a beneficial microorganism, its application in the control of ornamental flower diseases has not yet been reported.
In this study, antagonistic bacteria against a wide range of fungal pathogens, including F. oxysporum and B. cinerea, were isolated from the soil. Among these, B. velezensis E had the widest spectrum of antifungal activity and exhibited the best control on lily basal rot and gray mold. B. velezensis E promoted plant growth and produced active substance, as was detected using plant growth and biochemical assays. Whole-genome analysis of strain E revealed the presence of genes involved in pathogens biocontrol and plant growth promotion. Overall, this research showed that B. velezensis E can potentially be used to control lily basal rot and gray mold and promote lily bulblet morphogenesis in scale propagation.
Materials and methods
Isolation and purification of bacteria from soil samples
Soil samples for bacterial isolation were obtained from approximately 5 cm below the soil surface in the continuous lily cropping field of Yanqing District, Beijing, China. Five grams of soil sample were suspended in 50 ml of sterile water and were incubated for 30 min at 28°C with shaking (150 rpm). After 10-fold serial dilutions, 0.5 ml of suspension was placed on Luria-Bertani (LB) plate and incubated at 28°C. Bacteria from the single colonies were further purified by streaking onto LB agar plates. Bacterial isolates were inoculated in LB broth medium and incubated overnight at 28°C and 180 rpm. After incubation, the isolates were mixed with an equal volume of 40% glycerol and stored at –80°C for further experimentation.
Antagonism assay against phytopathogenic fungi and bacteria
Bacterial isolates were assessed for antifungal activities against Fusarium oxysporum, Fusarium equiseti, Botrytis cinerea, Rizoctonia solani, Alternaria solani, and Phytophthora capsica using dual culture assay. A mycelium plug (0.5 cm diameter) of actively growing fungi was placed at the center of a Potato Dextrose Agar (PDA) plate. The suspension of bacterial isolates (OD600 = 1, 5 μl for each point) was then inoculated onto the PDA plate at three points at a distance of 2 cm from the mycelium plug, followed by growth at 28°C in the dark. Antifungal activity was evaluated using Singh’s method (Singh et al., 2013) until the mycelia of fungal pathogens grew to the edge of the Petri dish in the control. Three replicates were used for each treatment.
The antagonistic activities of the bacterial isolates against Xanthomonas arboricola and Pantoea agglomerans were determined using a disc diffusion assay with minor modifications (Durairaj et al., 2018). A total of 200 μl of overnight pathogenic bacterial culture grown in LB broth was evenly coated on the surface of the LB solid medium. Then, 5 μl of the fermentation suspension of bacterial isolates (OD600 = 1) was added to sterile filter paper with a diameter of 6 mm and placed on the culture medium with pathogenic bacteria. Sterile water was used as a negative control.
Inhibitory effect of bacterial culture filtrate (BCF) on spore germination, germ tube elongation, and mycelial growth of fungal pathogens
To prepare spore suspensions of fungal pathogens, 0.5 cm mycelial plugs of B. cinerea and F. oxysporum were placed on PDA medium and incubated at 28°C for 10 days. The fully sporulated fungal plate was washed with sterile distilled water containing 0.03% Tween 80 and gently scratched with sterile spatula. The spore suspension was filtered using sterile gauze and adjusted to 1 × 105 spores/ml by diluting with Potato Dextrose Broth (PDB).
To prepare BCF, a single colony of B. velezensis E was pre-inoculated in LB broth medium and incubated overnight at 28°C in a shaking incubator set at 180 rpm. Then, 100 μl of broth culture was inoculated into LB broth medium (100 ml) and incubated at 28°C and 180 rpm for 2 days. The bacterial suspension was then centrifuged at 12,000 rpm for 10 min. The bacterial supernatant was passed through 0.22 μm pore-size filter membranes to obtain a BCF without bacterial cells.
To estimate the inhibitory effect of BCF on fungal pathogens, the obtained BCF was added to the spore suspension of the fungal pathogens to obtain final concentration of 10% BCF. LB broth was added to the spore suspension of the fungal pathogen as a control. A mixture of fungal spores and BCF was incubated at 28°C and 180 rpm. Germination of 100 spores and germ tube lengths of 20 spores from each treatment were observed and recorded after 10 h. Morphological changes in the fungal hyphae were observed and photographed after incubation for 3 days. The fungal hyphae were harvested using Whatman No.1 filter paper after the mixture of fungal spores and BCF was incubated for 2 days and then dried at 60°C for 2 days, followed by the measurement of dry weight of hyphae.
Identification of bacterial isolate by 16S rDNA sequence analysis
Genomic DNA was extracted from the isolate using the Bacteria Genomic DNA Extraction Kit (Tiangen Corporation Ltd., Beijing, China) according to the manufacturer’s instructions. For 16S rRNA gene amplification, genomic DNA was amplified using the universal bacterial primers F27 (5′-AGAGTTTGATCCTGGCTGGCTCAG-3′) and R1492 (5′ -TACGGCTACCTTGTTACGACTT-3′). The PCR mix contained 2 μl of DNA template (100 ng/μl), 2 μl each of F27 and R1492 (10 μm), 25 μl of 2 × Taq PCR Green Mix (Ding-guo Changsheng Corporation Ltd., Beijing, China), and topped to a final volume of 50 μl using RNase-free water. The PCR reaction mix was incubated at 94°C for 5 min, and then 94°C for 1 min, 55°C for 1 min, 72°C for 2 min for 30 cycles, followed by a final extension at 72°C for 10 min. The PCR product was verified by agarose gel electrophoresis and purified using a TIANgel Midi Purification Kit (Tiangen Corporation Ltd., Beijing, China) according to the manufacturer’s instructions. The purified product was cloned into the pTOPO-TA simple vector (Aidlab Corporation Ltd., Beijing, China) and sequenced at BGI (Shenzhen, China). The obtained sequence (about 1,500 bp) was aligned to the GenBank database using NCBI BLAST.
Evaluation of biocontrol activities and plant growth-promoting traits of bacterial strain in vitro
Siderophore production: The Chrome Azurol S (CAS) assay was used to detect siderophore production (Pérez-Miranda et al., 2007). All glassware used for the detection of siderophores was washed with 6 M HCl (Esmaeel et al., 2016). Ten microliters of freshly grown bacterial culture was inoculated on CAS agar plates, followed by incubation at 28°C for 4 days. The formation of a yellow-orange halo around the bacterial colony indicated siderophore production. Fresh LB broth without bacteria was used as a negative control.
Phosphate solubilization: Inorganic phosphorus (IP medium) and organic phosphorus solid medium (OP medium) were used to determine inorganic and organic phosphate solubilization, respectively (Teng et al., 2019). Ten microliters of freshly grown bacterial culture was inoculated on IP or OP medium, followed by incubation at 28°C for 4 days. The formation of a halo zone around the bacterial colony in IP or OP medium indicated inorganic or organic phosphate solubilization.
Proteolytic activity: Skim milk agar plates were used to test the ability of antagonistic bacteria to produce protease. This experiment was performed using a previous described procedure with minor modifications (Adhikari et al., 2017). The bacteria were dipped from freshly grown single colonies with sterile toothpicks and inoculated onto skim milk agar plates, followed by incubation at 28°C for 2 days. The development of a clear halo around the bacterial colony indicated protease production.
Cellulase production: Carboxymethyl cellulose (CMC) agar was used to determine the ability of antagonistic bacteria to produce cellulase; this experiment was performed using a previously described method with minor modifications (Syed-Ab-Rahman et al., 2018). The bacteria were dipped from freshly grown single colonies with sterile toothpicks and inoculated onto CMC agar plates, followed by incubation at 28°C for 4 days. The Congo red solution (1 g/l) was poured onto the surface of the plate for 60 min. Finally, unbound Congo red was removed by rinsing the plate with NaCl solution (1 M). The formation of a clear halo around the bacterial colony confirmed cellulase production.
Motility assay: The swarming and swimming motility of antagonistic bacteria were analyzed using previously reported methods (Ding et al., 2021; Thomloudi et al., 2021). A total of 10 ml of bacteria were grown to mid-log phase at 28°C in LB broth, centrifuged at 4,000 rpm, and resuspended in 100 μl of PBS. Subsequently, a 5 μl drop of bacterial suspension was inoculated at the center of freshly prepared LB plates containing 0.7 or 0.3% agar for swarming and swimming assays. After incubation at 30°C for 8 h, the appearance of a bacterial zone was observed.
Biofilm production: The ability of antagonistic bacteria to form biofilms was determined according to Rondeau’s method (Rondeau et al., 2019) with slight modifications. The antagonistic bacteria were inoculated into 5 ml of LB broth and incubated overnight at 28°C. The bacterial culture was diluted to an OD 600 of 0.2 with LB broth and added to a sterile polypropylene tube (300 μl per tube). After incubation at 28°C for 1 day without agitation, the medium was removed and washed with distilled water thrice. Subsequently, 600 μl of 1% (w/v) crystal violet was added to each tube and incubated at room temperature for 30 min. The crystal violet was removed from the tube by rigorous washing with distilled water thrice. Following this, 1.2 ml of 100% methanol was added to each tube. After incubation for 30 min, the biofilm formation was quantified by measuring the absorbance at 540 nm. Five replicates were used for each treatment.
Genome sequencing, assembly, annotation, and analysis
The genomic DNA of strain E was extracted using the cetyltrimethylammonium bromide (CTAB) method (Wilson, 2001). The concentration of genomic DNA was determined using a Qubit fluorometer (Invitrogen, United States). Sample integrity and purity were detected using electrophoresis on 1% agarose gel and NanoDrop 2000 Spectrophotometer (Thermo Scientific, United States). The genome of strain E was sequenced using the PacBio sequel II and DNBSEQ platforms at the Beijing Genomics Institute (BGI, Shenzhen, China). After removing low-quality raw reads using SOAPnuke (v1.5.6; Chen et al., 2018), high-quality reads were assembled into a corrected circular consensus sequence using Canu (v1.5; Koren et al., 2017), Proovread (Hackl et al., 2014), and GATK.1 Gene prediction was performed on the assembled genome sequence of strain E using glimmer32 with Hidden Markov. tRNAscan-SE (Lowe and Eddy, 1997), RNAmmer (Lagesen et al., 2007), and the Rfam database3 were used to identify tRNA, rRNA, and sRNAs. Tandem repeat annotation was performed using the Tandem Repeat Finder,4 and minisatellite DNA and microsatellite DNA were selected based on the number and length of repeat units. Functional annotation was accomplished by analysis of protein sequences. To obtain their corresponding annotations, all genes from strain E were aligned with different databases, including the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), Non-Redundant Protein Database databases (NR), Swiss-Prot, Gene Ontology (GO), and InterPro using Diamond software.
The phylogenetic identity of strain E was established by calculating digital DNA–DNA hybridization (dDDH) using the genome-to-genome distance calculator website service (GGDC 2.1; Meier-Kolthoff et al., 2013)5 and orthologous average nucleotide identity (OrthoANI) using the web service EzBioCloud (Yoon et al., 2017).6 The OrthoANI values 95–96% and dDDH values 70% were considered as species delineation threshold values (Meier-Kolthoff et al., 2013; Yoon et al., 2017). Phylogenomic analysis of strain E was carried out using the Type (strain) Genome Server7 and the tree was generated using FastME 2.1.6.1 (Lefort et al., 2015) from GBDP distances calculated from genome sequences.
Biosynthetic gene clusters (BGCs) of potential antimicrobial compounds were identified by antiSMASH 6.08 and BAGEL 4 (van Heel et al., 2018; Blin et al., 2021).9 BGCs that shared less than 70% amino acid identity with known clusters were considered as novel (Petrillo et al., 2021).
Biocontrol evaluation of strain E against basal rot and gray mold of lily
The seedlings of Lilium ‘siberia’ from scale cutting propagation were used to evaluate the control effect of strain E on lily base rot. The seedlings were planted in plastic pots (15 cm diameter, 15 cm height) containing peat moss at 25°C/20°C over a 16/8 h light/dark period. The pathogen F. oxysporum (accession number: OP216527) was isolated from lily plants with typical symptoms of lily basal rot and was conserved in the Lily Breeding Laboratory, College of Landscape Architecture, Beijing University of Agriculture (Li et al., 2016). Fusarium oxysporum spores were prepared according to the aforementioned method, and the spore suspension was adjusted to 1 × 106 spores/ml by diluting with sterile water. The seedlings of Lilium ‘siberia’ with a bulb diameter of about 1.5 cm were used for F. oxysporum and isolated strains inoculation. For each treatment, nine pots with one seedling in each pot were used for inoculation. After adding 40 ml of bacteria suspension (OD600 = 1) of strain E to the plant root for 3 days, each seedling was inoculated with 40 ml of suspension of F. oxysporum spore. The seedlings were inoculated with sterile water as control group. The symptoms of lily basal rot were recorded 6 weeks after inoculation. Disease severity was scored on a scale of 0–5 according to Shahin’s method (Shahin et al., 2011) with slight modifications: 0, health; 1, slight decay; 2, moderate decay; 3, severe decay; 4, very severe decay; 5, completely rotten. Disease incidence was calculated according to Ding’s method (Ding et al., 2021). Disease index (DI) was calculated using the following formula:
where N is the number of plants at each scale, DS is the disease scale, and T is the total number of plants. The relative biocontrol efficacy of the bacterial strains was calculated according to the following equation (Purkayastha et al., 2018):
To evaluate the ability of strain E to control gray mold in lily, leaves were excised from Lilium ‘siberia’ at the mature flowering stage and disinfected for 30 s in 75% (v/v) alcohol, followed by rinsing with sterile water thrice. Sterilized leaf segments were placed in Petri dishes with sterile gauze for moisturization. Five microliters of suspension of strain E (OD600 = 1) containing 0.03% Tween 80 were dropped to near the center of the leaf segment. Sterile water containing 0.03% Tween 80 was used as a negative control. The pathogen Botrytis cinerea (accession number: OP216528) was isolated from lily plants with typical symptoms of gray mold and conserved in the Lily Breeding Laboratory. Mycelium plugs (0.5 cm in diameter) of Botrytis cinerea excised from the actively growing edge of the colonies were inoculated onto the center of the leaf segments. Each treatment consisted of three replicates, and each replicate consisted of three leaf segments per plate. Inoculated leaf segments were incubated at 25°C/20°C (16 h/8 h, light/dark). Lesion diameter (LD) at the site of inoculation was measured using the average X and Y lengths after 4 days of incubation. The relative biocontrol efficacy of strain E was calculated using the following equation:
Promoting effect of strain E on the occurrence of scale bulblets
Suspensions of strain E were prepared according to the above method and adjusted to 1 × 106 CFU/ml, 1 × 107 CFU/ml, and 1 × 108 CFU/ml using sterile water. Lily scales were disinfected in 2% sodium hypochlorite solution for 5 min, followed by rinsing with sterile water thrice. Sterilized lily scales were planted in sterilized peat, which was inoculated with 200 ml bacterial suspension per kg of peat. Sterile water was used as a negative control. After 7 weeks of growth at 22°C in the dark, the diameter, root length, and weight of scale bulblets were recorded. Each treatment consisted of three replicates, and each replicate consisted of three lily scales in one pot.
Construction of GFP-labeled strain E and colonization assay
The plasmid pGFP78, an Escherichia coli–Bacillus subtilis shuttle vector containing the GFP gene controlled by the GFP78 promoter, was electroporated into B. velezensis strain E, according to a previously reported method (Zhang et al., 2017; Ding et al., 2021). GFP-labeled B. velezensis strain E was selected on LB agar medium supplemented with tetracycline (10 μg/ml) and examined using a fluorescence microscope (Zeiss Axio Scope.A1, Carl Zeiss, Jena, Germany).
Suspensions of GFP-labeled strain E (1 × 10 9 CFU/ml) were prepared according to the aforementioned method. The Lilium ‘Siberia’ bulbs were soaked in the suspension of GFP-labeled strain E for 60 min. Sterile water was used as control. Each treatment consisted of three replicates, and each replicate consisted of three plantlets. The root colonization of strain E was detected using previously reported methods with minor modifications (Thomas, 2004; Valdameri et al., 2015; Zhang et al., 2022). To evaluate the epiphytic colonization ability of strain E, 1 g of the lily root from plantlets was collected, cleaned with sterilized water to remove the peat on the root surface, and then agitated in a vortex for 1 min in 5 ml of sterilized water containing 0.03% Tween 80 seven days after inoculation. To evaluate the endophytic colonization ability of strain E on lily root, 1 g of root was disinfected by washing for 30 s in 75% ethanol solution and 3 min in 5% NaClO solution, rinsed three times with sterilized water, and then ground in a sterile mortar and suspended in 9 ml of sterile water. The bacterial suspension obtained was serially diluted and plated on LB agar plates containing tetracycline (10 μg/ml). After incubation overnight at 28°C, the number of bacteria (CFUs) was counted and recorded.
Lilium ‘siberia’ seedling planted for 2 months were used to evaluate the ability of leaf colonization of strains E. The leaves of lily seedlings were sprayed with bacterial suspension (1 × 10 9 CFU/ml) of GFP-labeled strain E with 0.03% Tween 80. Sterile water was used as control. Each treatment consisted of three replicates, and each replicate consisted of three seedlings. The leaves were surface washed with sterilized water and dried on sterilized filter paper 7 days after inoculation. Then, 1 g of leaf from each seedling was ground in a sterile mortar and suspended in 9 ml of sterile water. The suspension was diluted tenfold with sterile water and plated on LB agar plates with tetracycline (10 μg/ml), followed by overnight incubation at 28°C for colony count (CFUs).
Statistical analysis
The data obtained in the study were analyzed using the statistical software SPSS 19.0 (SPSS Inc., Chicago, IL, United States). Duncan’s multiple range test was performed to determine significant differences among the mean values (p < 0.05).
Results
Isolation and identification of antagonistic bacteria against Fusarium oxysporum
Seven bacterial isolates with antagonistic activity against F. oxysporum were obtained. The mycelial growth inhibition rate of these bacterial isolates ranged from 42.7 to 86.0% using dual culture assay. Among these, strain E exhibited the highest antifungal activity against F. oxysporum (Table 1). The cell-free supernatants of strain E inhibited spore germination, germ tube length and hyphal dry weight of F. oxysporum (Figures 1A–C), and caused swelling of hyphae (Figure 2). Strain E was classified as Bacillus verezensis based on the BLASTn analysis of the 16S rRNA gene sequence. It shared more than 99% identity with several Bacillus verezensis strains.
Table 1. In vitro antagonistic activity of bacterial strains isolated from soil against Fusarium oxysporum.
Figure 1. Growth inhibition of Fusarium oxysporum and Botrytis cinerea by cell-free supernatants of B. verezensis strain E. (A–C) Fusarium oxysporum. (D–F) Botrytis cinerea. (A,D) Rate of spore germination. (C,E) Germ tube length. (D,F) Hyphae dry weight. LB broth was added to the spore suspension of the fungal pathogen as a control. Asterisks indicate significant difference according to Duncan’s multiple range test. (**p < 0.05).
Figure 2. Morphological changes observed in hyphae of fungal pathogens after treatment with strain E. LB broth was used as a control.
Antagonistic activity of isolated bacteria against other plant pathogenic fungi and bacteria
Using a dual culture assay, all seven bacterial strains tested against F. oxysporum were found to display antifungal activity against five economically important fungal pathogens of landscape plants. Among them, B. verezensis E had the widest spectrum of antifungal activity and exhibited a good inhibitory effect on Fusarium equiseti, Botrytis cinerea, Rizoctonia solani, Alternaria solani, and Phytophthora capsica, with percentage of growth inhibition ranging from 86.1 to 97.7% (Figure 3A). The phenomena of death, winding, curling, knotting or swelling of the hyphe of F. equiseti, B. cinerea and A. solani were observed after treatment with cell-free supernatants of strain E (Figure 2). The cell-free supernatants of strain E inhibited spore germination, germ tube length, and hyphal dry weight of B. cinerea (Figures 1D–F). The disc diffusion method was used to investigate the in vitro antagonistic effects of strain E on two economically important bacterial pathogens, Xanthomonas arboricola, and Pantoea agglomerans (Yang et al., 2011; Lamichhane, 2014; Fan et al., 2022). After 24 h of incubation, clear inhibition zones around the colonies of strain E were observed on LB plates coated with Xanthomonas arboricola and Pantoea agglomerans (Figure 3B).
Figure 3. Antagonistic activities of isolated bacterial strains against economically important fungal and bacterial pathogens. (A) Antagonistic activity of seven bacterial strains against five fungal pathogens using dual culture assay. Three replicates of each treatment were used. Fungal pathogens cultured on PDA plates was used as control. Bars labeled with different letters indicate significant difference (p < 0.05) according to Duncan’s multiple range test. (B) Antagonistic activity of strain E against two bacterial pathogens using disc diffusion assay. Sterile water was used as a control.
Effects of Bacillus verezensis E on controlling lily basal rot and gray mold
Three bacterial strains (E, 11, and Z1) with the highest antagonistic activity against F. oxysporum were selected to evaluate their ability to control lily basal rot in greenhouses. Strain E exhibited the highest disease control efficacy (80.2%). The disease control efficacies of strains 11 and Z1 were only 54.2 and 64.0%, respectively (Figures 4A–C). The biocontrol activities of B. verezensis E against lily gray mold were further evaluated using in vitro detached leaf assays. Lily leaves from the control group showed spot symptoms 24 h after inoculation with B. cinerea. However, the lily leaves treated with strain E did not show disease symptoms until 48 h after inoculation. The diameter of disease spots in leaves treated with strain E was significantly smaller than that of the control 4 days after inoculation with B. cinerea. Therefore, the disease control efficacy of B. verezensis strain E could reach 83.0% (Figures 4D,E).
Figure 4. Efficacy of isolated bacterial strains in controlling lily basal rot and gray mold caused by F. oxysporum and B. cinerea, respectively. (A) Lily basal rot. (B) Lily gray mold. (C) Disease index and (D) control efficacy of bacterial strains against lily basal rot. Bars labeled with different letters are significantly different (p < 0.05) according to Duncan’s multiple range test. (E) Lesion diameter of lily leaf after treatment with ddH2O and strain E. Asterisks indicate statistically significant difference according to Duncan’s multiple range test. (**p < 0.05).
Bacillus verezensis E promotes the occurrence of scale Bulblets
After the lily scales were treated with a bacterial solution of strain E for 6 weeks, the diameter, root length, and fresh weight of scale bulblets in the treatment groups were larger than those in the control group (Figures 5A–D). Among the scales treated with 1 × 107 CFU/ml bacterial suspension product, the scale bulblets had the heaviest fresh weight, largest diameter, and longest root (Figures 5B–D).
Figure 5. Strain E promotes the occurrence of scale bulblets. (A) The occurrence of scale bulblets after treatment with 1 × 106 CFU/ml, 1 × 107 CFU/ml, and 1 × 108 CFU/ml bacterial suspension. (B) Diameter of scale bulblets. (C) Root length of scale bulblets. (D) Fresh weight of scale bulblets. Sterile water was used as a control. Bars labeled with different letters are significantly different (p < 0.05) according to Duncan’s multiple range test.
Evaluation of biocontrol activities and plant growth stimulating traits of Bacillus verezensis E in vitro
Some traits related to biocontrol and plant growth were also detected in strain E. It produced siderophores, proteases and cellulases (Figures 6A,C,D). A very slight halo zone around the bacterial colony in the organic phosphorus solid medium showed that strain E could slightly solubilize organic phosphate (Figure 6B), but not inorganic phosphate (data not shown). The strain E displayed excellent swimming motility, and the plate was almost fully covered by the swimming cells after incubation at 30°C for 11 h (Figure 6E). It also exhibited swarming motility, showing a swarming zone of 3.5 cm diameter after incubation at 30°C for 11 h (Figure 6F), and could completely cover the plate in 33 h. In addition, the crystal violet staining assay showed that strain E formed biofilms in polypropylene tube (Figure 7A).
Figure 6. Plant growth promoting and environmental fitting traits of strain E. (A) Siderophore production. (B) Organic phosphate solubilization. (C) Protease production. (D) Cellulase production. (E) Swimming motility. (F) Swarming motility. All assays were performed in three or five repetitions in three independent experiments.
Figure 7. Assays of colonization and biofilm formation of strain E (A) Evaluation of the biofilm-forming capacity by crystal violet staining. LB broth was used as a control. Epiphytic (B) and endophytic (C) colonization of GFP-labeled strain E on the lily root. Sterile water was used as a control. Asterisks indicate statistically significant differences according to Duncan’s multiple range test. (**p < 0.05).
Colonization of Bacillus verezensis E on lily root and leaf
To investigate root and leaf colonization of B. velezensis strain E, a GFP-labeled strain E was constructed. GFP-labeled strain E exhibited strong colonization activity on the lily root. For epiphytic colonization, the GFP-labeled strain E (1.8 × 106 CFU per gram root) was isolated on the root surface at 7 days after inoculation (Figure 7B). For endophytic colonization, the density of GFP-labeled strain E was 0.5 × 106 CFU per gram of surface-disinfected root (Figure 7C). However, no colonies of GFP-labeled strain E were detected from the lily leaves using our method (data not shown). Therefore, the density of GFP-labeled strain E was less than 100 CFU per gram of leaf after inoculation for 7 days.
Sequencing and analysis of Bacillus verezensis E genome
To investigate the biocontrol mechanisms of B. velezensis strain E, its complete genome was sequenced and analyzed. The genome sequence of strain E consisted of a complete circular chromosome (Figure 8), which was 3,929,247 bp in size, and had a GC content of 47.32%. No plasmid was detected in strain E. The whole genome of strain E contained 4,056 protein-coding genes (88.76% of the genome length), 27 tRNA, 86 rRNA genes, and 272 repeat elements. The total length of protein-coding genes was 3,487,731 bp with an average length of 859.89 bp (Table 2).
Figure 8. Genomic completion map of strain E. The size and scale of genome is shown in the outermost circle. Outer to inner rings represent forward and reverse strand genes annotated according to COG, forward and reverse strand ncRNA, repeat, GC content and GC-SKEW.
Table 2. Statistics of strain E genome assembly.
Previously, strain E was identified as B. velezensis based on 16S rDNA gene sequences. To further verify this classification, the whole genome sequence of strain E was compared with those of other Bacillus strains using the Type (strain) Genome Server. Strain E was found to belong to B. velezensis (Figure 9). Strain E was further compared with other closely related Bacillus strains at the genome level using ANI and dDDH analyses. The OrthoANI (98.73, 98.22, 98.56, and 98.20%) and dDDH (89.0, 84.7, 84.6, and 84.3%) values between strain E and four Bacillus velezensis strains were higher than the published species thresholds (95–96% for OrthoANI, 70% for dDDH; Table 3; Meier-Kolthoff et al., 2013; Yoon et al., 2017). Because the analysis results from the ANI and dDDH methods were consistent with those from the genome-based phylogenetic analysis, strain E was classified as B. velezensis.
Figure 9. Sequence-based phylogenomic tree of strain E constructed on TYGS. Type strains are indicated by T. Reference sequences were obtained from the GenBank database. Bootstrapping was performed using 1,000 replicates.
Table 3. Pairwise dDDH and orthoANI values between strain E and other Bacillus strains.
Analysis of biosynthetic gene clusters of secondary metabolites in Bacillus verezensis E
To better understand this antagonistic activity, BGCs of secondary metabolites were predicted and located in the genome of strain E using antiSMASH. Twelve putative BGCs with potential antimicrobial activities were identified. Among them, seven BGCs encoding synthetases of surfactin, macrolactin H, bacillaene, fengycin, difficidin, bacillibactin and bacilysin were highly homologous to those of known compounds in the MiBIG database and had similar gene structures to other B. velezensis strains (Table 4; Supplementary Figure 1). The remaining five BGCs involved in the biosynthesis of terpenes, PKS-like, T3PKS, and lanthipeptide-class-ii did not match the clusters of known compounds in the MiBIG database (Table 4). However, they were common in other B. velezensis strains based on ClusterBlast analysis of antiSMASH, with more than 95% gene similarity (Supplementary Figure 1).
Table 4. Analysis of secondary metabolite biosynthetic gene clusters (BGCs) in strain E.
To detect potential novel bacteriocin BGCs, the genome sequence of B. velezensis strain E was mined using the BAGEL4. Three areas of interest (AOI) were predicted to produce potential bacteriocins (Supplementary Figure 2). Among them, bacteriocin BGCs from AOI_01 and AOI_03 corresponded to those of regions11 and region 4 from antiSMASH analysis, respectively. AOI_02, encoding LCI, was identified using BAGEL4 instead of antiSMASH. One product of AOI_03, which shared ≤ 70% of similarity with known sequences in BAGEL4 database, were classified as novel bacteriocin (Supplementary Figure 2).
Analysis of Bacillus verezensis E genes involved in colonization, plant growth, and biological control
The genes involved in biocontrol, plant growth, and colonization were identified in strain E (Supplementary Table 1). To exert beneficial effects on plants, bacteria must colonize the plants. Many studies have demonstrated that bacterial colonization depends on their motility, chemotaxis, and biofilm formation abilities (Beauregard et al., 2013; Allard-Massicotte et al., 2016; Xu et al., 2019). A set of genes involved in biofilm matrix formation (epsA-O, tapA-sipW-tasA, and bslA) and regulation (spo0A, sinI, sinR, ylbF, ymcA, and yaaT) were identified from the genome of strain E. Bacteria must adhere to plant surfaces before biofilm formation (Xu et al., 2019). The genome of strain E contained a complete set of genes implicated in chemotaxis and motility, including cheR, cheV, hfq, motA, motB, swrA, and swrB, six genes encoding methyl-accepting chemotaxis proteins, and 30 genes in the fla/che operon (EGL001674-EGL001705). Beneficial bacteria can promote plant growth and development by improving the solubilization of mineral nutrients and synthesizing specific compounds to promote plant growth (Bhattacharyya and Jha, 2012; Aeron et al., 2020; Samaras et al., 2021). Multiple genes associated with plant growth-promoting activity in strain E were predicted, including the production of indole-3-acetic acid (IAA)(trpABFCDE operon, dhaS, and ysnE), siderophores (dhbF, dhbB, dhbE, entC, and entA), organic phosphorus mineralization (phoA, phoB, phoD, phoR, phoP, and pstB), and nitrogen metabolism (nos, narGHJI operon, nasA, nirD, nirB, and nirK).
Discussion
Lily basal rot and lily gray mold are important diseases that endanger lily production (Cui et al., 2018; He et al., 2019). Although biological control has been widely used to control many plant diseases, the biological control of lily diseases has not been extensively studied. Previously, one endophytic fungus, Acremonium sp., and three endophytic bacteria, Bacillus velezensis Lle-9, Bacillus stratosphericus LW-03, and Paenibacillus polymyxa SK1, were isolated from Lilium davidii, Lilium leucanthum, Lilium wardii, and Lilium lancifolium, respectively, all of which showed antifungal effects against Fusarium oxysporum, Botrytis cinerea, Botryosphaeria dothidea, and Fusarium fujikuroi in vitro (Khan et al., 2020a,b,c, 2021). These four microorganisms promoted lily growth, including increases in plant height, leaf length, and root length, which were significantly higher than those of the non-inoculated control plants. However, the effects of these four antagonistic microorganisms on diseases affecting lily have not been studied. In the present study, a new biocontrol, B. velezensis strain E, with broad-spectrum antimicrobial activity against filamentous fungi and Gram-negative bacteria, was isolated from a continuous lily cropping field. Strain E could effectively control gray mold and basal rot of lily in greenhouse experiments and in vitro detached leaf assays, respectively. The scale cutting propagation of lily is a low-cost and easy operation; however, it has a low propagation coefficient and bulblet formation is slow (Sun et al., 2009). In this study, B. velezensis strain E could effectively promote in vitro regeneration of bulblets by scale cutting, which is an important way to produce lily bulbs. Hence, B. velezensis strain E could potentially be used not only as a BCA for lily to control basal rot and gray mold disease, but also as a microbial fertilizer in lily production.
Bacillus species, which are widely distributed in the environment, have been extensively studied, because they can produce biologically active secondary metabolites that inhibit pathogen growth and promote plant growth (Zahir et al., 2003; Rabbee et al., 2019). Bacillus velezensis, a latetype isoform of B. amyloliquefaciens, is a new type of biocontrol bacterium (Dunlap et al., 2016; Yan et al., 2021). Bacillus velezensis has wide-spectrum biological activities against various plant pathogens, such as Phytophthora infestans, Pectobacterium carotovorum, Glomerella glycines, Rhizoctonia solani, Alternaria alternate, and Blumeria graminis (Cai et al., 2017; Saravanan et al., 2021; Yan et al., 2021; Htwe Maung et al., 2022). However, there is no report of the antagonistic effect of B. velezensis on Fusarium equiseti, Xanthomonas arboricola, and Pantoea agglomerans. Among them, Pantoea agglomerans can cause bacterial shot-hole disease of plum, bacterial soft rot of Dioscorea opposita, blight disease on Solanum muricatum and oat, necrotic disease on Ziziphus jujuba, brown apical necrosis of walnut, and bacterial shot-hole disease of plum (Yang et al., 2011; She et al., 2019; Moloto et al., 2020; She et al., 2021; Fan et al., 2022; Wang et al., 2022). Xanthomonas arboricola can cause bacterial canker and spot and bacterial blight diseases of stone fruit, almond, and walnut trees, respectively (Lamichhane, 2014). Both of them are economically important bacterial pathogens and cause huge economic losses. In this study, B. velezensis strain E could effectively antagonize these pathogens. Our results provide a potentially microbial resource for the biocontrol of diseases caused by these pathogens.
To exert beneficial effects on plants, bacteria must colonize the plant. The GFP-labeled B. velezensis E could not be clearly distinguished from the roots and leaves of lily using confocal microscopy because of the autofluorescence of lily roots and leaves. Therefore, antibiotics were added to the plates to determine the colonization density of GFP-labeled strain E in the lily roots and leaves. In this study, lily roots could be colonized epiphytically and endophytically by strain E after 7 days of inoculation. This result was consistent with the previous report (Synek et al., 2021). In that study endophytic colonization of Arabidopsis root was found to occur by the passive and active entry of a beneficial endophyte Enterobacter sp. SA187. And the endophyte SA187 could be attracted to the meristematic zone of growing roots and colonized epiphytically in the elongation zone. Therefore, we believe that strain E is an endophytic bacterium. Although the fermentation suspension of strain E had a good control effect against Botrytis cinerea in in vitro detached leaf assays, it could not be detected on lily leaves in the colonization experiment. However, the BCF of strain E has a strong destructive effect on Botrytis cinerea; therefore, the practical application effect of strain E on controlling lily gray mold must be studied further.
Data availability statement
The data presented in the study are deposited in the NCBI database, accession number ON764297 (16S rDNA sequence of strain E, https://www.ncbi.nlm.nih.gov/nuccore/ON764297.1/) and CP099460 (whole genome sequence of strain E, https://www.ncbi.nlm.nih.gov/nuccore/CP099460).
Author contributions
XH planned the experiments and wrote the manuscript. XH, BL, JW, JS, MT, MW, XZ, WS, YM, JM, YQ, YL, and WW performed the experiments. XH, BL, JW, and JS analyzed the data. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the Beijing Municipal Education Commission (KM201810020009) and the Research Fund for Academic Degree and Graduate Education of Beijing University of Agriculture (2022YJS037).
Acknowledgments
We gratefully acknowledge Qi Wang (China Agricultural University, Beijing, China) for supplying the plasmid pGFP78 used in this study.
Conflict of interest
The 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.976918/full#supplementary-material
Footnotes
1. ^https://www.broadinstitute.org/gatk/
2. ^http://www.cbcb.umd.edu/software/glimmer/
4. ^http://tandem.bu.edu/trf/trf.html
5. ^http://ggdc.dsmz.de/home.php
References
Adeniji, A. A., and Babalola, O. O. (2022). Evaluation of Pseudomonas fulva PS9.1 and bacillus velezensis NWUMFkBS10.5 as candidate plant growth promoters during maize–fusarium interaction. Plants-Basel 11:324. doi: 10.3390/plants11030324
Adhikari, M., Yadav, D. R., Kim, S. W., Um, Y. H., Kim, H. S., Lee, S. C., et al. (2017). Biological control of bacterial fruit blotch of watermelon pathogen (Acidovorax citrulli) with rhizosphere associated bacteria. Plant Pathol. J. 33, 170–183. doi: 10.5423/PPJ.OA.09.2016.0187
Aeron, A., Khare, E., Jha, C. K., Meena, V. S., Aziz, S. M. A., Islam, M. T., et al. (2020). Revisiting the plant growth-promoting rhizobacteria: lessons from the past and objectives for the future. Arch. Microbiol. 202, 665–676. doi: 10.1007/s00203-019-01779-w
Allard-Massicotte, R., Tessier, L., Lécuyer, F., Lakshmanan, V., Lucier, J. F., Garneau, D., et al. (2016). Bacillus subtilis early colonization of Arabidopsis thaliana roots involves multiple chemotaxis receptors. MBio 7:e01664-16. doi: 10.1128/mBio.01664-16
Aulia Rahma, A., Suryanti, S., Somowiyarjo, S., and Joko, T. (2020). Induced disease resistance and promotion of shallot growth by bacillus velezensis B-27. Pak. J. Biol. Sci. 23, 1113–1121. doi: 10.3923/pjbs.2020.1113.1121
Bayisa, R. A. (2020). Enhancing resistance of Sesamum indicum against Alternaria sesami through bacillus velezensis AR1. Pest Manag. Sci. 76, 3577–3586. doi: 10.1002/ps.5890
Beauregard, P. B., Chai, Y., Vlamakis, H., Losick, R., and Kolter, R. (2013). Bacillus subtilis biofilm induction by plant polysaccharides. Proc. Natl. Acad. Sci. U. S. A. 110, E1621–E1630. doi: 10.1073/pnas.1218984110
Bhattacharyya, P. N., and Jha, D. K. (2012). Plant growth-promoting rhizobacteria (PGPR): emergence in agriculture. World J. Microbiol. Biotechnol. 28, 1327–1350. doi: 10.1007/s11274-011-0979-9
Blin, K., Shaw, S., Kloosterman, A. M., Charlop-Powers, Z., van Wezel, G. P., Medema, M. H., et al. (2021). AntiSMASH 6.0: improving cluster detection and comparison capabilities. Nucleic Acids Res. 49, W29–W35. doi: 10.1093/nar/gkab335
Cai, X. C., Liu, C. H., Wang, B. T., and Xue, Y. R. (2017). Genomic and metabolic traits endow bacillus velezensis CC09 with a potential biocontrol agent in control of wheat powdery mildew disease. Microbiol. Res. 196, 89–94. doi: 10.1016/j.micres.2016.12.007
Chen, Y. X., Chen, Y. S., Shi, C. M., Huang, Z. B., Zhang, Y., Li, S. K., et al. (2018). SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience 7, 1–6. doi: 10.1093/gigascience/gix120
Choub, V., Ajuna, H. B., Won, S. J., Moon, J. H., Choi, S. I., Maung, C. E. H., et al. (2021). Antifungal activity of bacillus velezensis CE 100 against anthracnose disease (Colletotrichum gloeosporioides) and growth promotion of walnut (Juglans regia L.) trees. Int. J. Mol. Sci. 22:10438. doi: 10.3390/ijms221910438
Cui, Q., Yan, X., Gao, X., Zhang, D. M., He, H. B., and Jia, G. X. (2018). Analysis of WRKY transcription factors and characterization of two Botrytis cinerea-responsive LrWRKY genes from Lilium regale. Plant Physiol. Biochem. 127, 525–536. doi: 10.1016/j.plaphy.2018.04.027
De Benito Armas, A., Padula, N. L., Setlow, B., and Setlow, P. (2008). Sensitization of Bacillus subtilis spores to dry heat and desiccation by pretreatment with oxidizing agents. Lett. Appl. Microbiol. 46, 492–497. doi: 10.1111/j.1472-765X.2008.02344.x
de Cáceres, N., González, F. F., Davey, M. R., Cancho Sanchez, E., and Wilson, Z. A. (2015). Conferred resistance to Botrytis cinerea in Lilium by overexpression of the RCH10 chitinase gene. Plant Cell Rep. 34, 1201–1209. doi: 10.1007/s00299-015-1778-9
Ding, H. X., Mo, W. D., Yu, S., Cheng, H. H., Peng, L. J., and Liu, Z. Y. (2021). Whole genome sequence of bacillus velezensis strain GUMT319: a potential biocontrol agent against tobacco black shank disease. Front. Microbiol. 12:658113. doi: 10.3389/fmicb.2021.658113
Dunlap, C. A., Kim, S. J., Kwon, S. W., and Rooney, A. P. (2016). Bacillus velezensis is not a later heterotypic synonym of Bacillus amyloliquefaciens; Bacillus methylotrophicus, Bacillus amyloliquefaciens subsp. plantarum and ‘Bacillus oryzicola’ are later heterotypic synonyms of bacillus velezensis based on phylogenomics. Int. J. Syst. Evol. Microbiol. 66, 1212–1217. doi: 10.1099/ijsem.0.000858
Durairaj, K., Velmurugan, P., Park, J. H., Chang, W. S., Park, Y. J., Senthilkumar, P., et al. (2018). Characterization and assessment of two biocontrol bacteria against Pseudomonas syringae wilt in Solanum lycopersicum and its genetic responses. Microbiol. Res. 206, 43–49. doi: 10.1016/j.micres.2017.09.003
Esmaeel, Q., Chevalier, M., Chataigné, G., Subashkumar, R., Jacques, P., and Leclère, V. (2016). Nonribosomal peptide synthetase with a unique iterative-alternative-optional mechanism catalyzes amonabactin synthesis in Aeromonas. Appl. Microbiol. Biotechnol. 100, 8453–8463. doi: 10.1007/s00253-016-7773-4
Fan, B., Wang, C., Song, X. F., Ding, X. L., Wu, L. M., Wu, H. J., et al. (2018). Bacillus velezensis FZB42 in 2018: the gram-positive model strain for plant growth promotion and biocontrol. Front. Microbiol. 9:2491. doi: 10.3389/fmicb.2018.02491
Fan, S. S., Zhou, F. Y., Xie, X. Y., Zhang, X. J., and Wang, T. L. (2022). First report of bacterial soft rot caused by Pantoea agglomerans on Chinese yam (Dioscorea opposita) in China. Plant Dis. 106:1977. doi: 10.1094/PDIS-10-21-2151-PDN
Hackl, T., Hedrich, R., Schultz, J., and Förster, F. (2014). Proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics 30, 3004–3011. doi: 10.1093/bioinformatics/btu392
Hamaoka, K., Aoki, Y., and Suzuki, S. (2021). Isolation and characterization of endophyte Bacillus velezensis KOF112 from grapevine shoot xylem as biological control agent for fungal diseases. Plants-Basel 10:1815. doi: 10.3390/plants10091815
Han, X. S., Shen, D. X., Xiong, Q., Bao, B. H., Zhang, W. B., Dai, T. T., et al. (2021). The plant-beneficial rhizobacterium bacillus velezensis FZB42 controls the soybean pathogen Phytophthora sojae due to bacilysin production. Appl. Environ. Microbiol. 87:e0160121. doi: 10.1128/AEM.01601-21
Harold, J. F. S., Fitt, B. D. L., and Landau, S. (1997). Temperature and effects of seed-borne Botrytis cinerea or Alternaria linicola on emergence of linseed (Linum usitatissimum) seedlings. J. Phytopathol. 145, 89–97. doi: 10.1111/j.1439-0434.1997.tb00369.x
He, X. F., Li, W. Y., Zhang, W. Z., Jin, X. T., Shenkute, A. G., Xu, S. F., et al. (2019). Transcriptome sequencing analysis provides insights into the response to Fusarium oxysporum in Lilium pumilum. Evol. Bioinforma. 15, 1–10. doi: 10.1177/1176934319838818
Hong, S., Kim, T. Y., Won, S. J., Moon, J. H., Ajuna, H. B., Kim, K. Y., et al. (2022). Control of fungal diseases and fruit yield improvement of strawberry using bacillus velezensis CE 100. Microorganisms 10:365. doi: 10.3390/microorganisms10020365
Htwe Maung, C. E., Choub, V., Cho, J. Y., and Kim, K. Y. (2022). Control of the bacterial soft rot pathogen, Pectobacterium carotovorum by Bacillus velezensis CE 100 in cucumber. Microb. Pathog. 173:105807. doi: 10.1016/j.micpath.2022.105807
Jiang, C. H., Yao, X. F., Mi, D. D., Li, Z. J., Yang, B. Y., Zheng, Y., et al. (2019). Comparative transcriptome analysis reveals the biocontrol mechanism of Bacillus velezensis F21 against fusarium wilt on watermelon. Front. Microbiol. 10:652. doi: 10.3389/fmicb.2019.00652
Khan, M. S., Gao, J., Chen, X., Zhang, M., Yang, F., Du, Y., et al. (2020a). The endophytic bacteria Bacillus velezensis Lle-9, isolated from Lilium leucanthum, harbors antifungal activity and plant growth-promoting effects. J. Microbiol. Biotechnol. 30, 668–680. doi: 10.4014/jmb.1910.10021
Khan, M. S., Gao, J., Chen, X., Zhang, M., Yang, F., Du, Y., et al. (2020b, 2020). Isolation and characterization of plant growth-promoting endophytic bacteria Paenibacillus polymyxa SK1 from Lilium lancifolium. Biomed. Res. Int. 2020:8650957. doi: 10.1155/2020/8650957
Khan, M. S., Gao, J., Munir, I., Zhang, M., Liu, Y., Moe, T. S., et al. (2021). Characterization of endophytic fungi, Acremonium sp., from Lilium davidii and analysis of its antifungal and plant growth-promoting effects. Biomed. Res. Int. 2021:9930210. doi: 10.1155/2021/9930210
Khan, M. S., Gao, J., Zhang, M., Chen, X., Moe, T. S., Du, Y., et al. (2020c). Isolation and characterization of plant growth-promoting endophytic bacteria Bacillus stratosphericus LW-03 from Lilium wardii. 3 Biotech. 10:305. doi: 10.1007/s13205-020-02294-2
Koren, S., Walenz, B. P., Berlin, K., Miller, J. R., Bergman, N. H., and Phillippy, A. M. (2017). Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736. doi: 10.1101/gr.215087.116
Lagesen, K., Hallin, P., Rødland, E. A., Staerfeldt, H. H., Rognes, T., and Ussery, D. W. (2007). RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–3108. doi: 10.1093/nar/gkm160
Lamichhane, J. R. (2014). Xanthomonas arboricola diseases of stone fruit, almond, and walnut trees: progress toward understanding and management. Plant Dis. 98, 1600–1610. doi: 10.1094/PDIS-08-14-0831-FE
Lefort, V., Desper, R., and Gascuel, O. (2015). FastME 2.0: a comprehensive, accurate, and fast distance-based phylogeny inference program. Mol. Biol. Evol. 32, 2798–2800. doi: 10.1093/molbev/msv150
Li, W. Y., Liu, Z. P., He, X. F., and Wang, W. H. (2016). The pathogen isolation, identification and pathogenicity of lily blight. J. Shenyang Agric. Univ. 47, 153–158. doi: 10.3969/j.issn.1000-1700.2016.02.005
Liang, L. Q., Fu, Y. J., Deng, S. S., Wu, Y., and Gao, M. Y. (2021). Genomic, antimicrobial, and aphicidal traits of bacillus velezensis ATR2, and its biocontrol potential against ginger rhizome rot disease caused by Bacillus pumilus. Microorganism 10:63. doi: 10.3390/microorganisms10010063
Liu, H. Y., Zeng, Q. C., Yalimaimaiti, N., Wang, W., Zhang, R. F., and Yao, J. (2021). Comprehensive genomic analysis of bacillus velezensis AL7 reveals its biocontrol potential against Verticillium wilt of cotton. Mol. Gen. Genomics. 296, 1287–1298. doi: 10.1007/s00438-021-01816-8
Lowe, T. M., and Eddy, S. R. (1997). tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955–964. doi: 10.1093/nar/25.5.955
Ma, S., Hu, Y., Liu, S., Sun, J., Irfan, M., Chen, L. J., et al. (2018). Isolation, identification and the biological characterization of Botrytis cinerea. Int. J. Agric. Biol. 20, 1033–1040. doi: 10.17957/IJAB/15.0600
Mao, T. T., Chen, X. J., Ding, H. X., Chen, X. Y., and Jiang, X. L. (2020). Pepper growth promotion and fusarium wilt biocontrol by Trichoderma hamatum MHT1134. Biocontrol Sci. Tech. 30, 1228–1243. doi: 10.1080/09583157.2020.1803212
Masmoudi, F., Tounsi, S., Dunlap, C. A., and Trigui, M. (2021). Endophytic halotolerant bacillus velezensis FMH2 alleviates salt stress on tomato plants by improving plant growth and altering physiological and antioxidant responses. Plant Physiol. Biochem. 165, 217–227. doi: 10.1016/j.plaphy.2021.05.025
Meier-Kolthoff, J. P., Auch, A. F., Klenk, H.-P., and Göker, M. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform. 14:60. doi: 10.1186/1471-2105-14-60
Miljakovi, D., Marinkovi, J., and Baleevi-Tubi, S. (2020). The significance of bacillus spp. in disease suppression and growth promotion of field and vegetable crops. Microorganisms 8:1037. doi: 10.3390/microorganisms8071037
Moloto, V. M., Goszczynska, T., Hassen, A. I., Pierneef, R., and Coutinho, T. (2020). Draft genome sequences of Pantoea agglomerans strains BD1274 and BD1212, isolated from onion seeds, reveal major differences in pathogenicity and functional genes. Microbiol. Resour. Announc. 9, 1–3. doi: 10.1128/MRA.01507-19
Moon, J. H., Won, S. J., Maung, C. E. H., Choi, J. H., Choi, S. I., Ajuna, H. B., et al. (2021). Bacillus velezensis CE 100 inhibits root rot diseases (Phytophthora spp.) and promotes growth of Japanese cypress (Chamaecyparis 19esame Endlicher) seedlings. Microorganisms 9:821. doi: 10.3390/microorganisms9040821
Pérez-Miranda, S., Cabirol, N., George-Téllez, R., Zamudio-Rivera, L. S., and Fernández, F. J. (2007). O-CAS, a fast and universal method for siderophore detection. J. Microbiol. Methods 70, 127–131. doi: 10.1016/j.mimet.2007.03.023
Petrillo, C., Castaldi, S., Lanzilli, M., Selci, M., Cordone, A., Giovannelli, D., et al. (2021). Genomic and physiological characterization of bacilli isolated from salt-pans with plant growth promoting features. Front. Microbiol. 12:715678. doi: 10.3389/fmicb.2021.715678
Purkayastha, G. D., Mangar, P., Saha, A., and Saha, D. (2018). Evaluation of the biocontrol efficacy of a serratia marcescens strain indigenous to tea rhizosphere for the management of root rot disease in tea. PLoS One 13:e0191761. doi: 10.1371/journal.pone.0191761
Rabbee, M. F., Ali, M. S., Choi, J., Hwang, B. S., Jeong, S. C., and Baek, K. H. (2019). Bacillus velezensis: a valuable member of bioactive molecules within plant microbiomes. Molecules 24:1046. doi: 10.3390/molecules24061046
Rondeau, M., Esmaeel, Q., Crouzet, J., Blin, P., Gosselin, I., Sarazin, C., et al. (2019). Biofilm-constructing variants of Paraburkholderia phytofirmans PsJN outcompete the wild-type form in free-living and static conditions but not in planta. Appl. Environ. Microbiol. 85, e02670–e02618. doi: 10.1128/AEM.02670-18
Ruiz-García, C., Béjar, V., Martínez-Checa, F., Llamas, I., and Quesada, E. (2005). Bacillus velezensis sp. Nov., a surfactant-producing bacterium isolated from the river Vélez in Málaga, southern Spain. Int. J. Syst. Evol. Microbiol. 55, 191–195. doi: 10.1099/ijs.0.63310-0
Sachin, N., Tsang, A., Shaanker, R. U., and Dayanandan, S. (2021). Genome sequence resource of bacillus velezensis EB14, a native endophytic bacterial strain with biocontrol potential against the poplar stem canker causative pathogen, Sphaerulina musiva. Phytopathology 111, 890–892. doi: 10.1094/PHYTO-09-20-0433-A
Samaras, A., Nikolaidis, M., Antequera-Gómez, M. L., Cámara-Almirón, J., Romero, D., Moschakis, T., et al. (2021). Whole genome sequencing and root colonization studies reveal novel insights in the biocontrol potential and growth promotion by Bacillus subtilis MBI 600 on cucumber. Front. Microbiol. 11:600393. doi: 10.3389/fmicb.2020.600393
Saravanan, R., Nakkeeran, S., Saranya, N., Senthilraja, C., Renukadevi, P., Krishnamoorthy, A. S., et al. (2021). Mining the genome of bacillus velezensis VB7 (CP047587) for MAMP genes and non-ribosomal peptide synthetase gene clusters conferring antiviral and antifungal activity. Microorganisms 9:2511. doi: 10.3390/microorganisms9122511
Shahin, A., Arens, P., Van Heusden, A. W., Van Der Linden, G., Van Kaauwen, M., Khan, N., et al. (2011). Genetic mapping in Lilium: mapping of major genes and quantitative trait loci for several ornamental traits and disease resistances. Plant Breed. 130, 372–382. doi: 10.1111/j.1439-0523.2010.01812.x
She, X. M., Yu, L., Lan, G. B., Tang, Y. F., Deng, M. G., Li, Z. G., et al. (2019). First report of necrotic disease caused by Pantoea agglomerans on Ziziphus jujuba in China. Plant Dis. 103:1405. doi: 10.1094/PDIS-01-19-0025-PDN
She, X. M., Yu, L., Lan, G. B., Tang, Y. F., Deng, M. G., Li, Z. G., et al. (2021). Pantoea agglomerans causing blight disease on Pepino melon (Solanum muricatum) in China. Crop Prot. 139:105385. doi: 10.1016/j.cropro.2020.105385
Shin, J. H., Park, B. S., Kim, H. Y., Lee, K. H., and Kim, K. S. (2021). Antagonistic and plant growth-promoting effects of bacillus velezensis BS1 isolated from rhizosphere soil in a pepper field. Plant Pathol. J. 37, 307–314. doi: 10.5423/PPJ.NT.03.2021.0053
Singh, R. K., Kumar, D. P., Solanki, M. K., Singh, P., Srivastava, S., Srivastva, A. K., et al. (2013). Multifarious plant growth promoting characteristics of chickpea rhizosphere associated bacilli help to suppress soil-borne pathogens. Plant Growth Regul. 73, 91–101. doi: 10.1007/s10725-013-9870-z
Straathof, T. P., and Löffler, H. J. M. (1994). Resistance to fusarium oxysporum at different developmental stages of Asiatic hybrid lilies. J. Am. Soc. Hortic. Sci. 119, 1068–1072. doi: 10.21273/JASHS.119.5.1068
Sun, H. M., Jia, Z. K., Yang, L. U., Tao, W. L., and Wang, C. X. (2009). Advances on the cutting propagation of scale in Lilium. Northern Hort. 2, 141–146. doi: 10.11937/bfyy.200902053
Syed-Ab-Rahman, S. F., Carvalhais, L. C., Chua, E., Xiao, Y., Wass, T. J., and Schenk, P. M. (2018). Identification of soil bacterial isolates suppressing different Phytophthora spp. and promoting plant growth. Front. Plant Sci. 9:1502. doi: 10.3389/fpls.2018.01502
Synek, L., Rawat, A., L'Haridon, F., Weisskopf, L., Saad, M. M., and Hirt, H. (2021). Multiple strategies of plant colonization by beneficial endophytic Enterobacter sp. SA187. Environ. Microbiol. 23, 6223–6240. doi: 10.1111/1462-2920.15747
Teng, Z. D., Chen, Z. P., Zhang, Q., Yao, Y., Song, M. Y., and Li, M. (2019). Isolation and characterization of phosphate solubilizing bacteria from rhizosphere soils of the Yeyahu wetland in Beijing. China. Environ. Sci. Pollut. Res. Int. 26, 33976–33987. doi: 10.1007/s11356-018-2955-5
Thomas, P. (2004). Isolation of Bacillus pumilus from in vitro grapes as a long-term alcohol-surviving and rhizogenesis inducing covert endophyte. J. Appl. Microbiol. 97, 114–123. doi: 10.1111/j.1365-2672.2004.02279.x
Thomloudi, E. E., Tsalgatidou, P. C., Baira, E., Papadimitriou, K., Venieraki, A., and Katinakis, P. (2021). Genomic and metabolomic insights into secondary metabolites of the novel bacillus halotolerans Hil4, an endophyte with promising antagonistic activity against gray mold and plant growth promoting potential. Microorganisms 9:2508. doi: 10.3390/microorganisms9122508
Valdameri, G., Kokot, T. B., Pedrosa, F. D. O., and De Souza, E. M. (2015). Rapid quantification of rice root-associated bacteria by flow cytometry. Lett. Appl. Microbiol. 60, 237–241. doi: 10.1111/lam.12351
Van Heel, A. J., De Jong, A., Song, C., Viel, J. H., Kok, J., and Kuipers, O. P. (2018). BAGEL4: a user-friendly web server to thoroughly mine RiPPs and bacteriocins. Nucleic Acids Res. 46, W278–W281. doi: 10.1093/nar/gky383
Wang, J. J., Chen, T. X., Xue, L. H., Wei, X. K., White, J. F., and Qin, Z. M.. (2022). A new bacterial leaf blight disease of oat (Avena sativa) caused by Pantoea agglomerans in China. Plant Pathol. 71, 470–478. doi: 10.1111/ppa.13479
Wang, G. F., Meng, J. F., Tian, T., Xiao, X. Q., Zhang, B., and Xiao, Y. N. (2020). Endophytic bacillus velezensis strain B-36 is a potential biocontrol agent against lotus rot caused by fusarium oxysporum. J. Appl. Microbiol. 128, 1153–1162. doi: 10.1111/jam.14542
Wilson, K. (2001). Preparation of genomic DNA from bacteria. Curr. Protoc. Mol. Biol. 56, 241–245. doi: 10.1002/0471142727.mb0204s56
Won, S. J., Moon, J. H., Ajuna, H. B., Choi, S. I., Maung, C. E. H., Lee, S., et al. (2021). Biological control of leaf blight disease caused by Pestalotiopsis maculans and growth promotion of Quercus acutissima Carruth container seedlings using bacillus velezensis CE 100. Int. J. Mol. Sci. 22:11296. doi: 10.3390/ijms222011296
Xu, Z. H., Zhang, H. H., Sun, X. L., Liu, Y., Yan, W. X., Xun, W. B., et al. (2019). Bacillus velezensis wall teichoic acids is required for biofilm formation and root colonization. Appl. Environ. Microbiol. 85, e02116–e02118. doi: 10.1128/AEM.02116-18
Yan, H. H., Qiu, Y., Yang, S. A., Wang, Y. Q., Wang, K. Y., Jiang, L. L., et al. (2021). Antagonistic activity of bacillus velezensis SDTB038 against Phytophthora infestans in potato. Plant Dis. 105, 1738–1747. doi: 10.1094/PDIS-08-20-1666-RE
Yang, K. Q., Qu, W. W., Liu, X., Liu, H. X., and Hou, L. Q. (2011). First report of Pantoea agglomerans causing brown apical necrosis of walnut in China. Plant Dis. 95:773. doi: 10.1094/PDIS-01-11-0060
Yoon, S. H., Ha, S. M., Kwon, S., Lim, J., Kim, Y., Seo, H., et al. (2017). Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67, 1613–1617. doi: 10.1099/ijsem.0.001755
Yuan, H. B., Shi, B. K., Wang, L., Huang, T. X., Zhou, Z. Q., Hou, H., et al. (2022). Isolation and characterization of bacillus velezensis strain P2-1 for biocontrol of apple postharvest decay caused by Botryosphaeria dothidea. Front. Microbiol. 12:808938. doi: 10.3389/fmicb.2021.808938
Zahir, Z. A., Arshad, M., and Frankenberger, W. T. (2003). Plant growth promoting rhizobacteria: applications and perspectives in agriculture. Adv. Agron. 81, 97–168. doi: 10.1016/S0065-2113(03)81003-9
Zhang, T., Xiong, J., Tian, R. C., Li, Y., Zhang, Q. Y., Li, K., et al. (2022). Effects of single-and mixed-bacterial inoculation on the colonization and assembly of endophytic communities in plant roots. Front. Plant Sci. 13:928367. doi: 10.3389/fpls.2022.928367
Keywords: lily basal rot, lily gray mold, Bacillus velezensis, genome analysis, plant growth promotion, biocontrol
Citation: Li B, Wan J, Sha J, Tian M, Wang M, Zhang X, Sun W, Mao Y, Min J, Qin Y, Liu Y, Wang W and He X (2022) Genomics assisted functional characterization of Bacillus velezensis E as a biocontrol and growth promoting bacterium for lily. Front. Microbiol. 13:976918. doi: 10.3389/fmicb.2022.976918
Edited by:
Yong Wang, Guizhou University, ChinaReviewed by:
Haixia Ding, Guizhou University, ChinaQassim Esmaeel, Université de Reims Champagne-Ardenne, France
Copyright © 2022 Li, Wan, Sha, Tian, Wang, Zhang, Sun, Mao, Min, Qin, Liu, Wang and He. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiangfeng He, aHhmNzkxMjMwQDE2My5jb20=
†These authors have contributed equally to this work