Li Zheng
Shilian Huang
Tom Hsiang
Guohui Yu1
Zide Jiang
Jianguang Li- 1Innovative Institute for Plant Health, College of Agriculture and Biology, Zhongkai University of Agriculture and Engineering, Guangzhou, China
- 2Key Laboratory of South Subtropical Fruit Biology and Genetic Resource Utilization (MOA), Guangdong Province Key Laboratary of Tropical and Subtropical Fruit Tree Research, Institute of Fruit Tree Research, Guangdong Academy of Agricultural Sciences, Guangzhou, China
- 3School of Environmental Sciences, University of Guelph, Guelph, ON, Canada
- 4Department of Plant Pathology, Guangdong Province Key Laboratory of Microbial Signals and Disease Control, South China Agricultural University, Guangzhou, China
- 5Chinese Academy of Tropical Agricultural Sciences Guangzhou Experimental Station, Guangzhou, China
Bacillus amyloliquefaciens has been widely used in the agriculture, food, and medicine industries. Isolate PP19 was obtained from the litchi fruit carposphere and showed biocontrol efficacy against litchi downy blight (LDB) whether applied preharvest or postharvest. To further understand the underlying regulatory mechanisms, the genome of PP19 was sequenced and analyzed. The genome comprised a 3,847,565 bp circular chromosome containing 3990 protein-coding genes and 121 RNA genes. It has the smallest genome among 36 sequenced strains of B. amyloliquefaciens except for RD7-7. In whole genome phylogenetic analysis, PP19 was clustered into a group with known industrial applications, indicating that it may also produce high-yield metabolites that have yet to be identified. A large chromosome structural variation and large numbers of single nucleotide polymorphisms (SNPs) between PP19 (industrial strain) and UMAF6639 (plant-associated strain) were detected through comparative analysis, which may shed light on their functional differences. Preharvest treatment with PP19 enhanced resistance to LDB, by decreasing the plant H2O2 content and increasing the SOD activity. This is the first report of an industrial strain of B. amyloliquefaciens showing a plant-associated function and with major potential for the biocontrol of LDB.
Introduction
Litchi (Litchi chinensis) belongs to Sapindaceae and is native to South China. Currently, it is mainly distributed in China, India, Thailand, Vietnam, South Africa, Australia, and other tropical and subtropical regions (Tindall, 1994). Litchi is known in China as the “King of fruit,” with its taste and its high nutritional value; the pulp is rich in nutrients such as sugars, organic acids, dietary fibers, proteins, amino acids, polyphenols, and minerals (Jiang et al., 2006). During postharvest storage of litchi, the fruit can easily turn brown and rot due to litchi anthracnose, sour rot, and downy blight. Among these diseases, litchi downy blight (LDB) caused by Peronophythora litchii is the most serious (Wang et al., 2013). LDB alone may destroy 20–30% of litchi fruit every year (Zheng et al., 2019).
At present, the disease is mainly controlled by chemical methods. Due to the shortcomings of chemical control, such as issues arising from fungicide residues, biological control has become a desirable alternative. Pterostilbene, a phenolic compound found in plants, was found to inhibit the germination of Peronophythora litchii sporangia and affect the growth of its mycelia by deforming and shrinking the mycelia and sporangia, and damaging the cell wall, plasma membrane and organelles (Xu et al., 2018). Similarly, isoliquiritin, which is one of the flavonoid compounds found in licorice (Glycyrrhiza uralensis), can also damage the cell membrane of P. litchii (Luo et al., 2016). Zeamines, produced by Dickeya zeae, can damage the endomembrane of P. litchii, and zeamine-treated postharvest litchi fruit showed lower levels of P. litchii infection (Liao et al., 2015). Volatile organic compounds (VOCs) from Streptomyces fimicarius BWL-H1 can damage the endomembrane system and cell wall of P. litchii (Xing et al., 2018). Postharvest treatment of litchi fruit with Bacillus amyloliquefaciens LY-1 resulted in a lower pericarp browning index and percentage of fruit decay and higher fruit quality (Wu et al., 2017).
Bacillus amyloliquefaciens is used for commercial production of BamH I restriction enzyme (Wilson and Young, 1975), the antibiotic barnase (Hartley and Rogerson, 1972) and other antibiotics (Yu et al., 2002), alpha amylase (Tanaka and Hoshino, 2003), and proteinase subtilisin (Peng et al., 2004). It can successfully colonize the epidermis of fruit and vegetables, compete with pathogens for nutrition, secrete antibacterial substances to inhibit the growth of pathogens, and induce plant defense systems to resist pathogen invasion for successful biological control (Hong et al., 2013). Increasing research has shown that B. amyloliquefaciens could be used in the biological control of postharvest diseases. B. amyloliquefaciens NCPSJ7 reduced the disease incidence of grape gray mold caused by Botrytis cinerea (Zhou et al., 2020). Iturin A produced by B. amyloliquefaciens BUZ-14 showed strong in vitro and in planta antifungal activity against B. cinerea, Monilinia fructicola, Monilinia laxa, Penicillium digitatum, Penicillium expansum and Penicillium italicum, which are major pathogens of postharvest fruit (Calvo et al., 2019). Several other strains of B. amyloliquefaciens are known to decrease postharvest diseases, such as 9001 against apple ring rot (Li et al., 2013), BA3 against pear gray mold (Qu et al., 2016), BGP20 against vegetable soft rot (Zhao et al., 2013), BUZ-14 against postharvest diseases of orange, apple, grape, and stone fruit (Calvo et al., 2017), CPA-8 against postharvest diseases of sweet cherry fruit (Gotor-Vila et al., 2017; Vilanova et al., 2018), and DH-4 against citrus green mold (Chen et al., 2018).
The genomes of several B. amyloliquefaciens strains have been sequenced, which allows for comparative genomic analysis and association of particular genes with particular functional uses. Based on marker gene and whole-genome sequence comparisons, B. amyloliquefaciens has been subdivided into two groups, called “industrial” and “plant-associated” (Borriss et al., 2011; Magno-Pérez-Bryan et al., 2015). The industrial group is commonly used as industrial producers of primary metabolites, including enzymes, purine nucleosides, riboflavin, among others (Zakataeva et al., 2007; Sheremet et al., 2011; Meng et al., 2019), while the plant-associated group is considered plant-associated due to the ability to promote plant growth and exhibit biocontrol activity against pathogens (Fan et al., 2012; Molinatto et al., 2016; Xie et al., 2020).
Bacillus amyloliquefaciens PP19 was isolated from the carposphere of litchi fruit pericarp and showed antifungal activity against Peronophythora litchii SC18, and LDB biocontrol with either preharvest or postharvest treatments (unpublished data). The VOCs produced by PP19, especially the compounds benzothiazole (BTH) and α-farnesene (AF), can inhibit the growth of P. litchii and reduce the severity of LDB (Zheng et al., 2019). The purpose of this study was to further investigate PP19, by sequencing its genome and characterizing its genes to shed light on its functional abilities.
Materials and Methods
Bacterial Culture and DNA Extraction
Bacillus amyloliquefaciens PP19 cells were cultivated at 28°C in 500 ml flasks, containing 100 ml of LB media for 24 h, and subsequently harvested through centrifugation at 10,000 × g for 5 min at 4°C. Genomic DNA was extracted immediately following Hoffman and Winston (1987). The purity and integrity of DNA were assessed by agarose gel electrophoresis and then quantified by QubitTM Fluorometer (Invitrogen, CA, United States).
Genome Sequencing, Assembly, and Annotation
Genome sequencing and assembly were performed at Novogene Bioinformatics Technology Co., Ltd. (Beijing, China), on the Nanopore PromethION. All reads were assessed by NanoPlot with the threshold value Q > 7, and then the clean data were assembled with Unicycler (Wick et al., 2017). The coding genes of PP19 were predicted using GeneMarkS (Version 4.17) (Besemer et al., 2001). Interspersed nuclear elements and tandem repeats (TRs) were identified by RepeatMasker (Version open-4.0.5) (Saha et al., 2008) and TRF (Tandem Repeats Finder, Version 4.07b) (Benson, 1999), respectively. Transfer RNA (tRNA), ribosomal RNA (rRNA), and small RNA (sRNA) genes were detected by tRNAscan-SE (Version 1.3.1) (Lowe and Eddy, 1997), rRNAmmer (Version 1.2) (Lagesen et al., 2007), and CMsearch (Version 1.1rc4) (Cui et al., 2016) with default parameters, respectively. The prediction of genomic islands (GIs) and prophages was done with IslandPath-DIOMB (Version 0.2) (Hsiao et al., 2003) and phiSpy (Version 2.3) (Zhou et al., 2011), respectively. The clustered regularly interspaced short palindromic repeat (CRISPR) sequences were predicted by the CRISPR digger (Ge et al., 2016). For functional annotation, the predicted protein sequences were compared against the Gene Ontology (GO) database (Ashburner et al., 2000), Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2006), Clusters of Orthologous Groups of Proteins (COG) database (Galperin et al., 2015), Non-Redundant Protein Database (NR) (Li et al., 2002), Transporter Classification Database (TCDB) (Saier et al., 2016), Swiss-Prot (Bairoch and Apweiler, 2000), and Carbohydrate-Active enZYmes Database (CAZy) (Cantarel et al., 2009).
Phylogenetic Analyses and Comparative Genome Analysis
Complete amino acid gene sets of 36 B. amyloliquefaciens strains were downloaded from the NCBI database.1 The whole genomes were aligned, and a genome-based phylogenetic tree was constructed using CVTree3 (Zuo and Hao, 2015). MUMmerV3.22 (Kurtz et al., 2004), and used for identification of homologous regions and for collinearity analysis with default parameters. LASTZ (Large-Scale Genome Alignment Tool) (Harris, 2007) was used for detecting InDels (insertions-deletions) and SV (chromosomal structural variation).
Evaluation of Biocontrol Efficacy of Preharvest Treatment of Litchi Fruit With PP19
The pathogen P. litchii SC18 and the biocontrol agent B. amyloliquefaciens PP19 were cultured following Zheng et al. (2019). During summer 2019 in Danzhou city, Hainan province, litchi fruit of cv. “Feizixiao” at approximately 60% ripening was sprayed with PP19 at 5 × 107 CFU/mL or 1/10 strength LB broth on the litchi trees until runoff. After 7 days, the fruit was harvested, and brought back to the laboratory, 30 litchi fruits were placed in each plastic container (32 × 22 × 10 cm), in which the bottom was covered with filter paper moistened with 15 mL sterile water. After 24 h, fruits were sprayed with 30–60 mL of P. litchii SC18 at 5 × 104 sporangia/mL per container, and maintained in a greenhouse at 25°C with a day/night cycle of 12/12 h. The disease index was calculated following Zheng et al. (2019).
Enzyme Activity Detection in the Pericarp
Fruit of litchi cv. “Feizixiao” was sampled in 2019. Nine fruits per treatment were sampled at each time point (0, 48, 84, 96, 108, and 120 hpi). Subsequently, the pericarps were collected and immediately packed with tin foil, frozen in liquid nitrogen, and stored at −80°C until the assays. The content of H2O2 and the activities of catalase (CAT, EC 1.11.1.6), superoxide dismutase (SOD, EC 1.15.1.1), and peroxidase (POD, EC 1.11.1.7) were determined with corresponding detection kits (Nanjing Jiancheng Biological Engineering Institute, Nanjing, China).
Results
Genome Sequencing of B. amyloliquefaciens PP19
The Nanopore PromethION produced 342,009 reads with average read length of 13,054 bp. The 4.46 billion base pairs of reads, were assembled using Unicycler to generate a circular chromosome 3,847,565 bp in length with 46.27% GC content with no plasmids and 1160-fold genome coverage (Figure 1B and Supplementary Table 1). GeneMarkS was used to predict 3,990 protein coding sequences (CDSs), accounting for 89.85% of the genome, with an average length of 866 bp per gene (Supplementary Table 1). The results of genome component analyses showed that there were 226 interspersed repeats, 193 TRs, 86 tRNA genes, 27 rRNA genes, 8 sRNA genes, 8 GIs, and 6 prophages with no CRISPR sequences (Supplementary Table 2).
Figure 1. Phylogenetic analysis and circular chromosome of strain PP19. (A) Phylogenetic tree of PP19 based on the whole genome sequences of 37 B. amyloliquefaciens strains with Bacillus subtilis used as an outgroup. (B) The circular genome map revealing the genetic basis of PP19. From outside to inside, the map shows the (1) position of the genome, (2) coding genes on the + and – strand, (3) COG annotation, (4) GO annotation, (5) non-coding RNA, and (6) GC content. The inner red indicates that the GC content in this region is lower than the genome-wide average GC content, while the outer green indicates the opposite. The higher the peak is, the greater the difference with the average GC content (7) GC skew value (G–C/G + C). When the value is positive, the CDS is more likely to be transcribed from the positive chain; otherwise, the CDS is more likely to be transcribed from the negative chain.
Functional annotation was based on diamond alignment of the protein sequences of predicted genes with several commonly used databases. In the COG analysis, 2,864 proteins were classified into 25 categories. Except for the prediction for the general function category (304 genes), the amino acid transport and metabolism (303 genes), and the transcription (279 genes) categories showed the highest percentages among 25 categories (Supplementary Figure 1).
GO analysis divided 2,684 proteins into three categories: biological processes, molecular functions, and cellular components. The metabolic process and cellular process, catalytic activity and binding, cell part and cell showed the highest gene abundance in the three categories, respectively (Supplementary Figure 2). In the KEGG pathway analysis, the metabolic pathway accounted for the highest proportion (Supplementary Figure 3). CAZy is a database of carbohydrate enzymes, including related enzyme families that catalyze the degradation, modification, and biosynthesis of carbohydrates. It contains five main categories: glycoside hydrolase (GH), glycosyl transferase (GT), polysaccharide lyase (PL), carbohydrate esterase (CE), and auxiliary activity (AA). There were 58, 52, 3, 15, 1, and 41 genes classified as GH, GT, PL, CE, AA, and CBM (carbohydrate-binding module) categories, respectively (Supplementary Figure 4). TCDB analysis indicated that the genes annotated as porters (169 genes, including uniporters, symporters, antiporters) and P-P-bond-hydrolysis-driven transporters (153 genes) accounted for 71.88% of all predicted transporters (448 genes) (Supplementary Figure 5).
Phylogenetic and Comparative Genomic Analyses
Thirty-six B. amyloliquefaciens genome sequences were used in phylogenetic analysis. The characteristics of these genomes were enumerated (Supplementary Table 3), however, no obvious differences between plant-associated and industrially relevant strains were found in terms of genome size, number of genes, proteins, rRNAs, or tRNAs, while plant-associated strains showed a higher level of GC content (Supplementary Table 3). Aside from strain RD7-7, PP19 showed the smallest genome size, smaller by 3.85 to 4.24 Mb, but this did not result in a decreased number of predicted genes or proteins (Supplementary Tables 1, 3).
The phylogenetic tree showed that PP19 clustered with the economically important strains named here as the Industrial Group (Figure 1A). However, PP19 was confirmed to have a biocontrol effect against LDB (unpublished data). To the best of our knowledge, this is the first report that a strain grouping with economically important industrial-use strains could be used in biocontrol of a plant disease. Comparative genomic analysis was applied to reveal the differences in the genome between PP19 and DSM7 (industrial strain characterized by its enormous potential to produce extracellular enzymes of industrials importance including amylases and proteases) and between PP19 and UMAF6639 (plant-associated strain). There were 174 indels detected between PP19 and DSM7, while 133 indels between PP19 and UMAF6639. A large number of SNPs was detected between PP19 and UMAF6639 (195,510), while only 43,045 were found between PP19 and DSM7 (Supplementary Table 4). Several long fragment deletions and insertions were detected between PP19 and DSM7, resulting in minor chromosome SV between them (Figure 2A). Chromosome SV analysis between PP19 and UMAF6639 indicated that almost half of their genomes were not matched (Figure 2B). Synteny and collinearity were examined between PP19 and DSM7 and between PP19 and UMAF6639. The results demonstrated that nearly all regions of sequence similarity fell along the diagonals of the forward strand, except for several gene sites, indicating a generally similar gene and sequence order between PP19 and DSM7 (Figures 2C,D). In a comparison of PP19 and UMAF6639 genomic sequences, several translocation sites were detected (Figures 2E,F).
Figure 2. Chromosome structural variation (SV) and synteny analysis between PP19 and DSM7, PP19 and UMAF6619. (A) SV analysis between PP19 and DSM7, (B) SV analysis between PP19 and UMAF6619, (C,E) synteny analysis between PP19 and DSM7, PP19 and UMAF6619 (parallel collinear graph), (D,F) synteny analysis between PP19 and DSM7, PP19 and UMAF6619 (two dimensional collinear graph). Collinear, collinear region; Translocation, translocation region; Inversion, inversion region; Tran + Inver, translocation and inversion region; Insertion, region where the inserted fragment is greater than or equal to 50 bp; Deletion, region where the deleted fragment is greater than or equal to 50 bp; ComplexInDel, region that is unmatched but located in the same region; Forward_chain, the forward chain of the genome sequence; Reverse_chain, the reverse chain of the genome sequence; Forward_CDS, CDS translating on the forward chain of genome sequence; Reverse_CDS, CDS translating on the reverse chain of genome sequence; Subjoin_Forward_CDS, supplementary CDS translating on the forward chain of genome sequence; Subjoin_Reverse_CDS, supplementary CDS translating on the reverse chain of genome sequence; when the matching sequence of the deleted gene or inserted gene reaches 100% consistency, the matching sequence is considered a supplementary CDS.
Preharvest Treatment With PP19 Suppressed LDB and Enhanced the Activity of Defense-Related Enzymes
Genome analysis further revealed the major traits of strain PP19 involved in the biocontrol activity against LDB. There were 12 gene clusters related to secondary metabolite biosynthesis in PP19 (Supplementary Table 5), including one lantipeptide, two non-ribosomal peptides (NRPS), one type III polyketide synthetase (PKS) and one other type of PKS, two terpenes, one TransAT-PKS, three TransAT-PKS-NRPS, and one other cluster (Supplementary Table 5). Therefore, it can be inferred that PP19 could suppress the growth of P. litchii by producing diverse secondary metabolites, and we observed in vitro antagonism (not published).
To validate its use in controlling LDB at the preharvest stage and the possible underlying mechanism, PP19 was tested on preharvest litchi fruit with clean LB solution as a control. After 7 days, the treated fruit was harvested and inoculated with spores of P. litchii SC18. The disease severity was recorded from 60 to 96 hours post-inoculation (hpi) at 12 h intervals, and the results showed a significantly lower disease index at 60, 72, 84, 96 hpi in fruit pretreated with PP19; however, the difference decreased with increasing time (Figures 3A,B).
Figure 3. Suppression of LDB and enhancement of defense-related enzyme activity by preharvest treatment with PP19. (A) Visual comparison between PP19 preharvest treatment and control groups at 72 hpi. (B) LDB disease index for fruit treated with PP19 solution at preharvest and those treated with the control LB media. (C–F) Comparison of the H2O2 contents and SOD, CAT and POD activities detected after P. litchii inoculation in the PP19 preharvest treatment and control groups.
To explore the effects of PP19 on plant resistance to LDB, the plant H2O2 contents and four related enzyme activities were assessed from the litchi pericarp. The H2O2 content in the litchi pericarp treated with PP19 was lower than that of the control after P. litchii inoculation. With increasing time, the H2O2 content of PP19-treated litchi pericarps showed a small increase, peaking at 108 hpi with a 2.1-fold increase, and then decreasing at 120 hpi. For the control, accelerated accumulation of H2O2 was detected before 84 hpi, followed by a slow increase (Figure 3C). The activity of SOD trended to decrease with increasing time in both the control and PP19-treated litchi pericarps, and the SOD activity in the PP19 treatment was higher than that in the control (Figure 3D). There was a slight increase in CAT activity before 84 and 96 hpi in the PP19 treatment and control groups, respectively, and the CAT activity decreased after peaking. The CAT activity in the PP19 treatment was higher than that in the control before 84 hpi but lower after 96 hpi (Figure 3E). The POD activity increased with increasing time in both the PP19 treatment and control groups, except at 120 hpi in the PP19 treatment, and the POD activity in the PP19 treatment was lower than that in the control, except at 108 hpi (Figure 3F).
Discussion
The industrial group of B. amyloliquefaciens is often used industrially, as these strains are usually isolated from particular environments or generated by traditional mutation selection. Strains XH7 and TA208 were obtained by conventional mutagenesis, and these showed high yields of purine nucleosides, guanosine and ribavirin (Yang et al., 2011; Zhang et al., 2011). Strain LL3 was isolated from fermented food (Korean bibimbap) and could synthesize medicinally used poly-γ-glutamic acid (Geng et al., 2011), while RD7-7 isolated from rice doenjang (Korean fermented soybean paste) showed antimicrobial activity against Bacillus cereus (Eom and Choi, 2016). Strain HK1 was isolated from used agricultural mulch film, and exhibited a greater biodegradation of agricultural plastic film (Zhang et al., 2018). Strain YP6 was isolated from the rhizosphere of Lolium perenne from a phosphorus mine, and the alkaline phosphatase produced by YP6 showed high potential for biodegradation of organophosphorus pesticides (Meng et al., 2019). In this study, we sequenced the whole genome of B. amyloliquefaciens strain PP19 which was isolated from the carposphere of litchi fruit pericarp. Surprisingly, PP19 showed biocontrol activity against LDB but its genome was more similar to industrially relevant strains rather than plant-associated strains based on whole genome sequence comparisons.
We speculated that there were differences in genome structure and gene function between plant-associated and industrially relevant strains. Plant-associated strains showed a higher GC content, while no significant differences in genome size and gene number were observed between them (Supplementary Table 3). Comparative genomic analysis confirmed major chromosome SV between PP19 (industrial strain) and UMAF6639 (plant-associated strain) (Figure 2), which may result in functional differences. Comparative genome analyses of 37 genomes of B. amyloliquefaciens revealed differences which might be reflected in gene function. Unique genes common to the plant-associated clades are mainly found into the following COG categories: metabolism of amino acid transport and metabolism (E), carbohydrate transport and metabolism (G), synthesis of secondary metabolites (Q), and general function prediction (R) (Magno-Pérez-Bryan et al., 2015). COG analysis of the PP19 genome showed a similar classification (Supplementary Figure 1), which may indicate functional similarity between this and other plant-associated strains and partly explain the plant-associated function.
A previous study indicated that postharvest treatment with PP19 showed significant biocontrol against fruit LDB (unpublished data), and pre-exposure of postharvest litchi fruit to the VOCs produced by PP19 significantly reduced the severity of LDB (Zheng et al., 2019). Here, we assessed the biocontrol efficacy of preharvest treatment with PP19 to further evaluate its field application prospects, and the result confirmed PP19 showed a great potential for postharvest preservation (Figures 3A,B). H2O2, a representative reactive oxygen species that can be overproduced when induced by various abiotic and biotic stresses, is highly reactive and can cause damage to plants (Mhamdi and Van Breusegem, 2018). SOD, CAT, and POD are important enzymatic antioxidants for scavenging of reactive oxygen species (Gill and Tuteja, 2010). PP19 pretreatment with litchi resulted in a greater ability to scavenge H2O2 (Figure 3C), mainly due to the higher activity of SOD instead of CAT and POD (Figures 3D–F), as the SOD always showed a higher activity than control. The activity of CAT showed a little fluctuation in both control and PP19 treatment, so it may play a little part in LDB resistance. The activity of POD increased substantially in both control and PP19 treatment, which may play an important role in LDB resistance while not in difference between control and PP19 treatment.
Industrial strains tend to produce higher levels of a metabolite: XH7 (purine nucleoside guanosine) (Yang et al., 2011), TA208 (guanosine and ribavirin) (Zhang et al., 2011), LL3 (poly-γ-glutamic acid) (Geng et al., 2011), and MT45 (surfactin) (Zhi et al., 2017). In addition, a large number of genes were classified into amino acid transport and metabolism (E), coenzyme transport and metabolism (H) and secondary metabolite biosynthesis, and transport and catabolism (Q) categories. We speculate that PP19 should also produce a high yield of certain metabolites contributing to its antagonism of P. litchii, which can enhance the antioxidant related enzymes activities, and that will be a focus of future research.
Conclusion
In conclusion, preharvest treatment with PP19 suppressed LDB, showing great potential for the exploitation and utilization of this strain. We sequenced the genome of B. amyloliquefaciens PP19, and whole genome phylogenetic analysis placed it with industrially relevant strains. This is the first report of a industrially relevant strain being used for plant disease control. Substantial chromosomal SV between PP19 and UMAF6639 was detected through comparative analysis, which may be one of the reasons for the functional differences between plant-associated and industrial strains. COG analysis indicated functional similarity between PP19 and plant-associated strains. From the above, it is speculated that a certain high-yield metabolite with inhibitory action against Peronophythora litchii was produced by PP19, and this metabolite should be further identified and studied.
Data Availability Statement
The accession number of Bacillus amyloliquefaciens PP19 genome sequence is CP062984. All other datasets for this study are included in the article/Supplementary Material.
Author Contributions
LZ, SH, GY, and DG conducted the experiments. LZ and SH analyzed the data and wrote the manuscript. TH revised the manuscript. ZJ and JL designed the experiments and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the earmarked fund for the Central Public-interest Scientific Institution Basal Research Fund for Chinese Academy of Tropical Agricultural Sciences (1630112019004), Guangzhou Science and Technology Project (202002030278), Research Start-up Funds of Zhongkai College of Agricultural Engineering (KA2005408C7), Guangdong Modern Agricultural Industrial Technology System Special Project (2018LM1027), and Special Funds for Rural Revitalization Strategy [403-2018-XMZX-0002-90 (TS-1-4-1)].
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.
Acknowledgments
We thank Jiang Lin at Haian University for his valuable suggestions on the trials and photography and Xi Zhang, Zheng Kang, Yun Luo, Maoqing Cai, and Yuedong Geng at Haian University for their assistance in field and laboratory trials. We also thank Xiaofu Guo in Danzhou City for his work in the field management of the litchis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.619423/full#supplementary-material
Footnotes
References
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Bairoch, A., and Apweiler, R. (2000). The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 28, 45–48. doi: 10.1093/nar/28.1.45
Benson, G. (1999). Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573–580. doi: 10.1093/nar/27.2.573
Besemer, J., Lomsadze, A., and Borodovsky, M. (2001). GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 29, 2607–2618. doi: 10.1093/nar/29.12.2607
Borriss, R., Chen, X. H., Rueckert, C., Blom, J., Becker, A., Baumgarth, B., et al. (2011). Relationship of Bacillus amyloliquefaciens clades associated with strains DSM 7T and FZB42T: a proposal for Bacillus amyloliquefaciens subsp. amyloliquefaciens subsp. nov. and Bacillus amyloliquefaciens subsp. plantarum subsp. nov. based on complete genome sequence comparisons. Int. J. Syst. Evol. Microbiol. 61(Pt 8), 1786–1801. doi: 10.1099/ijs.0.023267-0
Calvo, H., Marco, P., Blanco, D., Oria, R., and Venturini, M. E. (2017). Potential of a new strain of Bacillus amyloliquefaciens BUZ-14 as a biocontrol agent of postharvest fruit diseases. Food Microbiol. 63, 101–110. doi: 10.1016/j.fm.2016.11.004
Calvo, H., Mendiara, I., Arias, E., Blanco, D., and Venturini, M. E. (2019). The role of iturin A from B. amyloliquefaciens BUZ-14 in the inhibition of the most common postharvest fruit rots. Food Microbiol. 82, 62–69. doi: 10.1016/j.fm.2019.01.010
Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The carbohydrate-active enzymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238. doi: 10.1093/nar/gkn663
Chen, K., Tian, Z., Luo, Y., Cheng, Y., and Long, C. A. (2018). Antagonistic activity and the mechanism of Bacillus amyloliquefaciens DH-4 against citrus green mold. Phytopathology 108, 1253–1262. doi: 10.1094/PHYTO-01-17-0032-R
Cui, X., Lu, Z., Wang, S., Jing-Yan Wang, J., and Gao, X. (2016). CMsearch: simultaneous exploration of protein sequence space and structure space improves not only protein homology detection but also protein structure prediction. Bioinformatics 32, i332–i340. doi: 10.1093/bioinformatics/btw271
Eom, J. S., and Choi, H. S. (2016). Inhibition of Bacillus cereus growth and toxin production by Bacillus amyloliquefaciens RD7-7 in fermented soybean products. J. Microbiol. Biotechnol. 26, 44–55. doi: 10.4014/jmb.1509.09090
Fan, B., Carvalhais, L. C., Becker, A., Fedoseyenko, D., von Wirén, N., and Borriss, R. (2012). Transcriptomic profiling of Bacillus amyloliquefaciens FZB42 in response to maize root exudates. BMC Microbiol. 12:116. doi: 10.1186/1471-2180-12-116
Galperin, M. Y., Makarova, K. S., Wolf, Y. I., and Koonin, E. V. (2015). Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 43, D261–D269. doi: 10.1093/nar/gku1223
Ge, R., Mai, G., Wang, P., Zhou, M., Luo, Y., Cai, Y., et al. (2016). CRISPRdigger: detecting CRISPRs with better direct repeat annotations. Sci. Rep. 6:32942. doi: 10.1038/srep32942
Geng, W., Cao, M., Song, C., Xie, H., Liu, L., Yang, C., et al. (2011). Complete genome sequence of Bacillus amyloliquefaciens LL3, which exhibits glutamic acid-independent production of poly-γ-glutamic acid. J. Bacteriol. 193, 3393–3394. doi: 10.1128/JB.05058-11
Gill, S. S., and Tuteja, N. (2010). Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48, 909–930. doi: 10.1016/j.plaphy.2010.08.016
Gotor-Vila, A., Teixidó, N., Di Francesco, A., Usall, J., Ugolini, L., Torres, R., et al. (2017). Antifungal effect of volatile organic compounds produced by Bacillus amyloliquefaciens CPA-8 against fruit pathogen decays of cherry. Food Microbiol. 64, 219–225. doi: 10.1016/j.fm.2017.01.006
Harris, R. S. (2007). Improved Pairwise Alignment of Genomic DNA. Ph.D. thesis, The Pennsylvania State Univeristy, University Park, PA.
Hartley, R. W., and Rogerson, D. L. Jr. (1972). Production and purification of the extracellular ribonuclease of Bacillus amyloliquefaciens (barnase) and its intracellular inhibitor (barstar). I. Barnase. Prep. Biochem. 2, 229–242. doi: 10.1080/00327487208061473
Hoffman, C. S., and Winston, F. (1987). A ten-minute DNA preparation from yeast efficiently releases autonomous plasmids for transformation of Escherichia coli. Gene 57, 267–272. doi: 10.1016/0378-1119(87)90131-4
Hong, P., An, G., Hu, M., and Wang, H. (2013). Advance in research on biological control of postharvest diseases of fruits and vegetables by Bacillus amyloliquefaciens. Chin. Agr. Sci. Bull. 29, 168–173. (in Chinese),Google Scholar
Hsiao, W., Wan, I., Jones, S. J., and Brinkman, F. S. (2003). IslandPath: aiding detection of genomic islands in prokaryotes. Bioinformatics 19, 418–420. doi: 10.1093/bioinformatics/btg004
Jiang, Y. M., Wang, Y., Song, L., Liu, H., Lichter, A., Kerdchoechuen, O., et al. (2006). Postharvest characteristics and handling of litchi fruit—an overview. Aust. J. Exp. Agr. 46, 1541–1556. doi: 10.1071/ea05108
Kanehisa, M., Goto, S., Hattori, M., Aoki-Kinoshita, K. F., Itoh, M., Kawashima, S., et al. (2006). From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 34, D354–D357. doi: 10.1093/nar/gkj102
Kurtz, S., Phillippy, A., Delcher, A. L., Smoot, M., Shumway, M., Antonescu, C., et al. (2004). Versatile and open software for comparing large genomes. Genome Biol. 5:R12. doi: 10.1186/gb-2004-5-2-r12
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
Li, W., Jaroszewski, L., and Godzik, A. (2002). Tolerating some redundancy significantly speeds up clustering of large protein databases. Bioinformatics 18, 77–82. doi: 10.1093/bioinformatics/18.1.77
Li, Y., Han, L. R., Zhang, Y., Fu, X., Chen, X., Zhang, L., et al. (2013). Biological control of apple ring rot on fruit by Bacillus amyloliquefaciens 9001. Plant Pathol. J. 29, 168–173. doi: 10.5423/PPJ.SI.08.2012.0125
Liao, L., Zhou, J., Wang, H., He, F., Liu, S., Jiang, Z., et al. (2015). Control of litchi downy blight by zeamines produced by Dickeya zeae. Sci. Rep. 5:15719. doi: 10.1038/srep15719
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
Luo, J., Li, Z., Wang, J., Weng, Q., Chen, S., and Hu, M. (2016). Antifungal activity of isoliquiritin and its inhibitory effect against Peronophythora litchi Chen through a membrane damage mechanism. Molecules 21:237. doi: 10.3390/molecules21020237
Magno-Pérez-Bryan, M. C., Martínez-García, P. M., Hierrezuelo, J., Rodríguez-Palenzuela, P., Arrebola, E., Ramos, C., et al. (2015). Comparative genomics within the Bacillus genus reveal the singularities of two robust Bacillus amyloliquefaciens biocontrol strains. Mol. Plant Microbe Interact. 28, 1102–1116. doi: 10.1094/MPMI-02-15-0023-R
Meng, D., Jiang, W., Li, J., Huang, L., Zhai, L., Zhang, L., et al. (2019). An alkaline phosphatase from Bacillus amyloliquefaciens YP6 of new application in biodegradation of five broad-spectrum organophosphorus pesticides. J. Environ. Sci. Heal. B 54, 336–343. doi: 10.1080/03601234.2019.1571363
Mhamdi, A., and Van Breusegem, F. (2018). Reactive oxygen species in plant development. Development 145:dev164376. doi: 10.1242/dev.164376
Molinatto, G., Puopolo, G., Sonego, P., Moretto, M., Engelen, K., Viti, C., et al. (2016). Complete genome sequence of Bacillus amyloliquefaciens subsp. plantarum S499, a rhizobacterium that triggers plant defences and inhibits fungal phytopathogens. J. Biotechnol. 238, 56–59. doi: 10.1016/j.jbiotec.2016.09.013
Peng, Y., Yang, X. J., Xiao, L., and Zhang, Y. Z. (2004). Cloning and expression of a fibrinolytic enzyme (subtilisin DFE) gene from Bacillus amyloliquefaciens DC-4 in Bacillus subtilis. Res. Microbiol. 155, 167–173. doi: 10.1016/j.resmic.2003.10.004
Qu, H., Zhao, L., Zhao, F., Liu, Y., and Yang, Z. (2016). Biocontrol of gray mold decay in pear by Bacillus amyloliquefaciens Strain BA3 and its effect on postharvest quality parameters. Pol. J. Microbiol. 65, 171–176. doi: 10.5604/17331331.1204476
Saha, S., Bridges, S., Magbanua, Z. V., and Peterson, D. G. (2008). Empirical comparison of ab initio repeat finding programs. Nucleic Acids Res. 36, 2284–2294. doi: 10.1093/nar/gkn064
Saier, M. H. Jr., Reddy, V. S., Tsu, B. V., Ahmed, M. S., Li, C., and Moreno-Hagelsieb, G. (2016). The transporter classification database (TCDB): recent advances. Nucleic Acids Res. 44, D372–D379. doi: 10.1093/nar/gkv1103
Sheremet, A. S., Gronskiy, S. V., Akhmadyshin, R. A., Novikova, A. E., Livshits, V. A., Shakulov, R. S., et al. (2011). Enhancement of extracellular purine nucleoside accumulation by Bacillus strains through genetic modifications of genes involved in nucleoside export. J. Ind. Microbiol. Biotechnol. 38, 65–70. doi: 10.1007/s10295-010-0829-z
Tanaka, A., and Hoshino, E. (2003). Secondary calcium-binding parameter of Bacillus amyloliquefaciens alpha-amylase obtained from inhibition kinetics. J. Biosci. Bioeng. 96, 262–267. doi: 10.1016/s1389-1723(03)80191-3
Tindall, H. D. (1994). Sapindaceous fruits: botany and horticulture. Hortic. Rev. 16, 143–196. doi: 10.1002/9780470650561.ch5
Vilanova, L., Teixidó, N., Usall, J., Balsells-Llauradó, M., Gotor-Vila, A., and Torres, R. (2018). Environmental fate and behaviour of the biocontrol agent Bacillus amyloliquefaciens CPA-8 after preharvest application to stone fruit. Pest Manag. Sci. 74, 375–383. doi: 10.1002/ps.4716
Wang, H., Qian, Z., Ma, S., Zhou, Y., Patrick, J. W., Duan, X., et al. (2013). Energy status of ripening and postharvest senescent fruit of litchi (Litchi chinensis Sonn.). BMC Plant Biol. 13:55. doi: 10.1186/1471-2229-13-55
Wick, R. R., Judd, L. M., Gorrie, C. L., and Holt, K. E. (2017). Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 13:e1005595. doi: 10.1371/journal.pcbi.1005595
Wilson, G. A., and Young, F. E. (1975). Isolation of a sequence-specific endonuclease (BamI) from Bacillus amyloliquefaciens H. J. Mol. Biol. 97, 123–125. doi: 10.1016/s0022-2836(75)80028-3
Wu, Y., Lin, H., Lin, Y., Shi, J., Xue, S., Hung, Y., et al. (2017). Effects of biocontrol bacteria Bacillus amyloliquefaciens LY-1 culture broth on quality attributes and storability of harvested litchi fruit. Postharvest Biol. Technol. 132, 81–87. doi: 10.1016/j.postharvbio.2017.05.021
Xie, S., Vallet, M., Sun, C., Kunert, M., David, A., Zhang, X., et al. (2020). Biocontrol potential of a novel endophytic bacterium from mulberry (Morus) tree. Front. Bioeng. Biotechnol. 7:488. doi: 10.3389/fbioe.2019.00488
Xing, M., Zheng, L., Deng, Y., Xu, D., Xi, P., Li, M., et al. (2018). Antifungal activity of natural volatile organic compounds against litchi downy blight pathogen Peronophythora litchii. Molecules 23:358. doi: 10.3390/molecules23020358
Xu, D., Deng, Y., Xi, P., Zhu, Z., Kong, X., Wan, L., et al. (2018). Biological activity of pterostilbene against Peronophythora litchii, the litchi downy blight pathogen. Postharvest Biol. Technol. 144, 29–35. doi: 10.1016/j.postharvbio.2018.05.011
Yang, H., Liao, Y., Wang, B., Lin, Y., and Pan, L. (2011). Complete genome sequence of Bacillus amyloliquefaciens XH7, which exhibits production of purine nucleosides. J. Bacteriol. 193, 5593–5594. doi: 10.1128/JB.05880-11
Yu, G. Y., Sinclair, J. B., Hartman, G. L., and Bertagnolli, B. L. (2002). Production of iturin A by Bacillus amyloliquefaciens suppressing Rhizoctonia solani. Soil Biol. Biochem. 34, 955–963. doi: 10.1016/s0038-0717(02)00027-5
Zakataeva, N. P., Gronskiy, S. V., Sheremet, A. S., Kutukova, E. A., Novikova, A. E., and Livshits, V. A. (2007). A new function for the Bacillus PbuE purine base efflux pump: efflux of purine nucleosides. Res. Microbiol. 158, 659–665. doi: 10.1016/j.resmic.2007.08.003
Zhang, G., Deng, A., Xu, Q., Liang, Y., Chen, N., and Wen, T. (2011). Complete genome sequence of Bacillus amyloliquefaciens TA208, a strain for industrial production of guanosine and ribavirin. J. Bacteriol. 193, 3142–3143. doi: 10.1128/JB.00440-11
Zhang, J., Chen, J., Jia, R., Dun, Z., Wang, B., Hu, X., et al. (2018). Selection and evaluation of microorganisms for biodegradation of agricultural plastic film. 3 Biotech 8:308. doi: 10.1007/s13205-018-1329-5
Zhao, Y., Li, P., Huang, K., Wang, Y., Hu, H., and Sun, Y. (2013). Control of postharvest soft rot caused by Erwinia carotovora of vegetables by a strain of Bacillus amyloliquefaciens and its potential modes of action. World J. Microbiol. Biotechnol. 29, 411–420. doi: 10.1007/s11274-012-1193-0
Zheng, L., Situ, J., Zhu, Q., Xi, P., Zheng, Y., Liu, H., et al. (2019). Identification of volatile organic compounds for the biocontrol of postharvest litchi fruit pathogen Peronophythora litchii. Postharvest Biol. Technol. 155, 37–46. doi: 10.1016/j.postharvbio.2019.05.009
Zhi, Y., Wu, Q., and Xu, Y. (2017). Genome and transcriptome analysis of surfactin biosynthesis in Bacillus amyloliquefaciens MT45. Sci. Rep. 7:40976. doi: 10.1038/srep40976
Zhou, Q., Fu, M., Xu, M., Chen, X., Qiu, J., Wang, F., et al. (2020). Application of antagonist Bacillus amyloliquefaciens NCPSJ7 against Botrytis cinerea in postharvest red globe grapes. Food Sci. Nutr. 8, 1499–1508. doi: 10.1002/fsn3.1434
Zhou, Y., Liang, Y., Lynch, K. H., Dennis, J. J., and Wishart, D. S. (2011). PHAST: a fast phage search tool. Nucleic Acids Res. 39, W347–W352. doi: 10.1093/nar/gkr485
Keywords: Bacillus amyloliquefaciens PP19, litchi downy blight, genome sequencing, comparative genomic analysis, preharvest treatment, biocontrol efficacy
Citation: Zheng L, Huang S, Hsiang T, Yu G, Guo D, Jiang Z and Li J (2021) Biocontrol Using Bacillus amyloliquefaciens PP19 Against Litchi Downy Blight Caused by Peronophythora litchii. Front. Microbiol. 11:619423. doi: 10.3389/fmicb.2020.619423
Received: 20 October 2020; Accepted: 08 December 2020;
Published: 12 January 2021.
Edited by:
Jacob George Malone, John Innes Centre, United KingdomReviewed by:
Oleg Reva, University of Pretoria, South AfricaJoachim Vater, Robert Koch Institute (RKI), Germany
Copyright © 2021 Zheng, Huang, Hsiang, Yu, Guo, Jiang and Li. 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: Zide Jiang, zdjiang@scau.edu.cn; Jianguang Li, ljg.cj.jr@163.com
†These authors have contributed equally to this work