Functional Analysis and Genome Mining Reveal High Potential of Biocontrol and Plant Growth Promotion in Nodule-Inhabiting Bacteria Within Paenibacillus polymyxa Complex

Bacteria belonging to the genus Paenibacillus were frequently isolated from legume nodules. The nodule-inhabiting Paenibacillus as a resource of biocontrol and plant growth-promoting endophytes has rarely been explored. This study explored the nodule-inhabiting Paenibacillus’ antifungal activities and biocontrol potentials against broad-spectrum important phytopathogenic fungi. We collected strains which were isolated from nodules of Robinia pseudoacacia, Dendrolobium triangulare, Ormosia semicastrata, Cicer arietinum, Acacia crassicarpa, or Acacia implexa and belong to P. peoriae, P. kribbensis, P. endophyticus, P. enshidis, P. puldeungensis, P. taichungensis, or closely related to P. kribbensis, or P. anseongense. These nodule-inhabiting Paenibacillus showed diverse antagonistic activities against five phytopathogenic fungi (Fusarium graminearum, Magnaporthe oryzae, Rhizoctonia solani, Sclerotinia sclerotiorum, and Botrytis cinerea). Six strains within the P. polymyxa complex showed broad-spectrum and potent activities against all the five pathogens, and produced multiple hydrolytic enzymes, siderophores, and lipopeptide fusaricidins. Fusaricidins are likely the key antimicrobials responsible for the broad-spectrum antifungal activities. The nodule-inhabiting strains within the P. polymyxa complex were able to epiphytically and endophytically colonize the non-host wheat plants, produce indole acetic acids (IAA), and dissolve calcium phosphate and calcium phytate. P. peoriae strains RP20, RP51, and RP62 could fix N2. P. peoriae RP51 and Paenibacillus sp. RP31, which showed potent plant colonization and plant growth-promotion competence, effectively control fungal infection in planta. Genome mining revealed that all strains (n = 76) within the P. polymyxa complex contain ipdC gene encoding indole-3-pyruvate decarboxylase for biosynthesis of IAA, 96% (n = 73) contain the fus cluster for biosynthesis of fusaricidins, and 43% (n = 33) contain the nif cluster for nitrogen fixation. Together, our study highlights that endophytic strains within the P. polymyxa complex have a high probability to be effective biocontrol agents and biofertilizers and we propose an effective approach to screen strains within the P. polymyxa complex.


INTRODUCTION
The intense use of agrochemicals including chemical pesticides, synthetic fertilizers, and plant growth regulators to increase crop yields to meet the increasing food demand for increasing population has intensified the side effects of agrochemicals on agriculture, human health, and ecosystems (Lamichhane et al., 2016;Carvalho, 2017). Some microbes naturally in association with plants are able to produce antimicrobials and plant growth regulators, provide nutrients to plants by nitrogen fixation and phosphate solubilization, or induce plant systemic resistance to biotic and abiotic stresses, and thus can protect crops and promote crop growth (Grady et al., 2016;Alori et al., 2017;Cassán et al., 2020). The microbe-based biological control agents and biofertilizers are ecofriendly alternatives to control plant diseases and promote crop growth and are growingly developed to reduce the use of agrochemicals and support the sustainable agriculture (Grady et al., 2016;Alori et al., 2017;Cassán et al., 2020).
Endophytes are microbes that reside inside plants for part or full of their life cycle and are apparently not harmful to the host plants (Wilson, 1995;Hallman et al., 1997). Endophytes are apparently competent plant colonizers and adaptive to niches inside plants and may get sufficient nutrition and protection from plants; they may adapt to plant immune response and form close association or mutualistic relationship with plants and thus can be potent biocontrol agents and plant growth promoters (Rosenblueth and Martínez-Romero, 2006;Hardoim et al., 2015;Rybakova et al., 2016).
Nodule-inhabiting Paenibacillus as a resource of biocontrol and plant growth-promoting endophytes has rarely been explored. The aim of this study was to screen nodule-inhabiting Paenibacillus strains having antagonistic activities against broadspectrum phytopathogenic fungi and explore their antifungal mechanisms, plant colonization and plant growth-promotion competences, and biocontrol potentials. Five important fungal pathogens Fusarium graminearum, Magnaporthe oryzae, Botrytis cinerea, Rhizoctonia solani, and Sclerotinia sclerotiorum were selected as target pathogens because they have broad plant hosts, cause devastating damage to major staple food crops and economic crops, currently require high doses of chemical fungicides for control, and thus demand effective biocontrol agents to reduce the use of chemical fungicides.

Amplification and Phylogenetic Analysis of 16S rRNA Gene Sequences of Paenibacillus Strains
Bacterial 16S rRNA gene sequences were amplified with the primers 27F (5 -AGAGTTTGATCCTGGCTCAG-3 ) and 1492R (5 -GGTTACCTTGTTACGACTT-3 ) from colonies grown on the LB agar as previously described (Ali et al., 2020). Nontype Paenibacillus strains were identified based on the identity between their16S rRNA gene sequences and those of the type strains at the EzBioCloud 1 and the phylogenetic status of their 16S rRNA gene sequences. Nucleotide sequences were aligned using the MUSCLE program; positions containing gaps and missing data were eliminated; final 1319 positions were constructed to a phylogenetic tree with the Neighbor-Joining method and the Kimura 2-parameter model integrated in the MEGA5 software (Tamura et al., 2011).

Screening Antagonistic Paenibacillus Against Fungal Pathogens
Paenibacillus antagonistic activities against fungal pathogens were examined by the bacteria-fungi confrontation assay on PDA. A 5-mm mycelial plug from the edge of a fresh 7d fungal colony on PDA was transferred to the center of a fresh PDA in a 90-mm Petri dish. Paenibacillus strains were cultured in the LB broth at 30 • C and 200 rpm for 48 h. Bacterial suspension was adjusted to 1 × 10 8 CFU ml −1 and 5 µl of the bacterial suspension was inoculated at 25 mm away from the fungal plug on the PDA plate. A fungal mycelial plug alone on the PDA plate was used as control. The PDA plates were kept at 28 • C for 7 days. Paenibacillus antagonistic activities were measured by inhibition of mycelial growth according to Riungu et al. (2008). Four replications were tested for each Paenibacillus strain and the experiment was repeated three times. 1 www.ezbiocloud.net

Assay of Hydrolytic Enzyme Activities From Paenibacillus Against Broad-Spectrum Fungi
Hydrolytic enzyme activities of the Paenibacillus strains inhibiting all the five tested fungal pathogens were assessed using minimal salt agar media supplemented with the enzyme substrate. The minimal salt medium per liter contains (NH 4 ) 2 SO 4 1 g, KH 2 PO 4 0.2 g, K 2 HPO 4 1.6 g, MgSO 4 ·7H 2 O 0.2 g, NaCl 0.1 g, FeSO 4 ·7H 2 O 0.01 g, and CaCl 2 ·2H 2 O 0.02 g. The medium for assessing chitinase activity was supplemented with 1% (w/v) colloidal chitin and 1.5% (w/v) agar and adjusted to the final pH 7.0 and sterilized by autoclaving at 121 • C for 15 min. The medium for assessing β-1,3-glucanase activity was supplemented with 0.5 g of glucose, 6.7 g of yeast extract, 60 mg of aniline blue and 4 g of pachyman powder (Megazyme Ltd., Wicklow, Ireland) per liter. The pH was adjusted to 6.8 and 1.2 g of agar was added before autoclaving at 121 • C for 10 min (Mahasneh and Stewart, 1980). The medium for assessing protease was supplemented with 1% (w/v) 1% (w/v) casein. Bacterial suspension was adjusted to 1 × 10 8 CFU ml −1 and 5 µl of the bacterial suspension was inoculated onto the agar plates and incubated at 28 • C for 2-3 days. The chitin plates were stained with 0.1% (w/v) Congo red and washed with 1% NaCl. Bacterial hydrolytic enzyme activities were indicated by hydrolysis zones around the bacterial colonies.
Colloidal chitin was prepared according to Souza et al. (2009). Five grams of chitin powder (J&K Scientific Ltd., Beijing, China) was added slowly to 60 ml of concentrated hydrochloric acid in an Erlenmeyer flask and kept with shaking at 180 rpm at 37 • C for 1 h. The mixture was filtered through glass wool and the filtrate was added to a 200 ml of 50% (v/v) ethanol with vigorous stirring. The precipitate was transfer to a glass funnel with filter paper and washed with sterile distilled water until the colloidal chitin reached a neutral pH 7. The colloidal chitin retained on the filter paper was collected, weighed and stored in dark at 4 • C before use.

Assay of Siderophore Production From Paenibacillus Against Broad-Spectrum Fungi
Siderophore production from the Paenibacillus strains inhibiting all the five tested fungal pathogens was assessed using the blue medium containing chrome azurol S (Schwyn and Neilands, 1987). Chelation of iron by siderophores was indicated by color changes from blue to purple (catechol type) or orange (hydroxamate type) in and around the bacterial colonies.

MALDI-TOF-MS Analysis for Lipopeptides From Paenibacillus Against Broad-Spectrum Fungi
Lipopeptides were detected by matrix-assisted laser desorption ionization-time of flight-mass spectrometry (MALDI-TOF-MS) (Masum et al., 2018). Bacteria were grown on LB agar at 30 • C for 48 h. A bacterial colony was suspended in a matrix solution [α-cyano-4-hydroxycinnamic acid (10 mg ml −1 ) in 30% (v/v) of acetonitrile and 70% (v/v) of 0.1% (w/v) TFA in water]. The suspension (1 µl) was spotted onto a target plate, air dried, and detected with a reflection-positive mode in the mass spectral range from 300 to 3000 Da using an Ultraflextreme MALDI-TOF mass spectrometer (Bruker, Bremen, Germany) equipped with a 355 nm nitrogen laser.

Scanning Electron Microscopy (SEM) and Transmission Electron Microscopy (TEM)
Blocks of 7-day-old fungal mycelia were taken from the edge of Fusarium graminearum mycelia grown on PDA alone or confrontation with P. peoriae RP51 and prepared for electron microscopy. The blocks were immersed in 2.5% (v/v) glutaraldehyde in 0.1 M phosphate buffer (pH 7.0) for 4 h, washed with the phosphate buffer for 15 min three times, then immersed in 1% (w/v) OsO 4 in the phosphate buffer for 2 h and washed three times with the phosphate buffer. The blocks were dehydrated in a graded series of ethanol for 15 min at each step and then in absolute acetone for 20 min twice. For SEM, the blocks were dehydrated in an HCP-2 critical point dryer (Hitachi, Tokyo, Japan) and then coated with gold-palladium in an E-1010 ion sputter (Hitachi, Tokyo, Japan) and observed with an Gemini SEM 300 scanning electron microscope (Carl Zeiss, Jena, Germany). For TEM, the dehydrated blocks were infiltrated with graded series of Spurr's resin mixed with absolute acetone and then the Spurr's resin overnight. The blocks embedded in the Spurr's resin were polymerized at 70 • C for 12 h. The specimen was sectioned with an EM UC7 ultratome (Leica Microsystems, Vienna, Austria). The ultrathin sections were stained with uranyl acetate and lead citrate and then observed with an H-7650 transmission electron microscope (Hitachi, Tokyo, Japan).

Assay of Plant Growth Promotion and Plant Colonization by Paenibacillus Against Broad-Spectrum Fungi
Wheat seeds (cv. Jimai22) were surface-sterilized by 70% ethanol for 1 min and 3% sodium hypochlorite for 6 min, and washed with sterile water for six times. Seeds were imbibed in sterile water overnight. The Paenibacillus strains against broad-spectrum fungi were grown in LB broth at 30 • C and 200 rpm for 48 h and then washed with sterile water and suspended to about 1 × 10 8 CFU ml −1 . The seeds (n = 100) were immersed in 10 ml of each bacterial suspension for 4 h or in sterile water as uninoculated control. The seeds were dried in a laminar flow hood and 15 seeds were placed on two layers of sterile 125mm Whatman filter paper moistened with 10 ml of sterile water in a 150-mm Petri dish for germination and growth at 25 • C under 16-h light and 8-h dark photoperiod. At 9 days after inoculation, root length, shoot length, fresh weight, and dry weight of the seedlings were measured. To determine epiphytic root colonization of the Paenibacillus strains, bacteria adhere to roots were washed off in 3 ml of sterile water and serially diluted; 100 µl of the bacterial suspensions were spread on LB agar. To determine endophytic colonization of the Paenibacillus strains in roots and stems, roots and stems were surface-sterilized by 70% ethanol for 1 min and 1% sodium hypochlorite for 2 min, washed with sterile water for six times. The sterilization efficacy was determined by no bacterial growth from the last washed water. The surface-sterilized roots and stems were ground in sterile water. The homogenate suspensions were serially diluted; 100 µl of the diluted suspensions were spread on LB agar. The LB agar plates were kept at 30 • C for 3 days and bacterial colonies showing the colony morphology of the inoculated strain were counted. The experiments were done with three replicates and repeated three times.

In planta Assay of Disease Control by Paenibacillus Against Fusarium
Fusarium root rot (FRR) and Fusarium foot rot (FFR) on wheat seedlings were used as the pathosystem (Colombo et al., 2019) to determine the biocontrol potential by selected Paenibacillus against broad-spectrum fungi in planta. Wheat seeds were surface-sterilized and inoculated by Paenibacillus, and prepared for germination as described above. Control seeds were treated with sterile water. Ten seeds were placed on the filter papers in a 150-mm Petri dish. When the primary roots of the seedlings reached about 30 mm, a 5-mm mycelial plug taken from the edge of an actively growing F. graminearum was put upside down on a primary root at a 10 mm distance from the seed. A sterile PDA plug was used as negative control. At 4 days after fungal inoculation, the extension (length) of the necrosis on the root (FRR) was measured (Supplementary Figure 1A). Inhibition of FRR by Paenibacillus was determined using the formula (Control necrosis -Paenibacillus-treated necrosis)/Control necrosis × 100 (Colombo et al., 2019). At 6 days after fungal inoculation, FFR was measured by scoring the symptoms at the crown level (0 = symptomless; 1 = slightly necrotic; 2 = moderately necrotic; 3 = severely necrotic; 4 = completely necrotic) (Supplementary Figure 1B)

Assays of Plant Growth-Promotion Traits of Paenibacillus Against Broad-Spectrum Fungi
Indole acetic acid (IAA) production by Paenibacillus was determined using the colorimetric assay developed by Sarwar and Kremer (1995). Paenibacillus strains were grown in LB broth supplemented with L-tryptophan (100 µg ml −1 ) (BBI Life Science, Shanghai, China) in dark at 30 • C for 3 days. After centrifugation, 150 µl of the culture supernatant was added into wells of 96-well microplate followed by addition of 100 µ1 of the Salkowski reagent and incubation in dark for 30 min. IAA (1 mg ml −1 ) (BBI Life Science, Shanghai, China) was diluted to 10 µg ml −1 to 50 µg ml −1 as standards and added into the wells. Each IAA standard and bacterial culture supernatant was tested in three replicate wells. After the 30-min reaction, the absorbance of the pink product was measured at 530 nm using a SpectraMax R Plus 384 Microplate Spectrophotometer (Molecular Devices Corporation, Sunnyvale, CA, United States). An IAA standard curve was generated and the IAA in the bacterial culture supernatant was quantified accordingly (Sarwar and Kremer, 1995). Bacterial solubilization of phosphate was determined by the assay on the NBRIP media (per liter contains glucose 10 g, MgCl 2 .6H 2 O 5 g, MgSO 4 .7H 2 O 0.25 g, KCl 0.2 g, (NH 4 ) 2 SO 4 0.1 g, agar 15 g; pH 7.0) containing 0.5% (w/v) calcium phosphate, ferric phosphate, or calcium phytate according to Nautiyal (1999) and Bashan et al. (2013).

Genome Mining
Gene clusters for biosynthesis of secondary metabolites (fusaricidins and siderophores) were mined from target whole genome sequences deposited in the NCBI database 2 as at 2020-07-30 using the antiSMASH 5.0 pipeline (Blin et al., 2019). Protein names "nitrogenase, " "ilvB, " and "1-aminocyclopropane-1-carboxylate deaminase" (ACC deaminase) were searched from the annotated protein coding genes (CDS) of the target whole genome sequence (Protein Table for the target organism) in the NCBI database (see text footnote 2). Gene ilvB encoding the biosynthetic-type acetolactate synthase large subunit is the gene ipdC encoding indole-3-pyruvate decarboxylase for the biosynthesis of indole acetic acid (IAA) via the IPyA pathway (Xie et al., 2016).

Statistical Analysis
All experiments were performed in completely randomized design. All the values were presented as mean ± standard error of at least three replications. One-way analysis of variance (ANOVA) was used to analyze the data following post hoc multiple comparisons using SPSS 16.0 software (SPSS Inc., Chicago, IL, United States). Least significant difference test was done to separate the treatments.

Nodule-Inhabiting Strains Within Paenibacillus polymyxa Complex Showed Broad-Spectrum Antifungal Activities
The phylogenetic tree of the 16S rRNA gene sequences of the nodule-inhabiting Paenibacillus strains and their relatives showed that P. polymyxa (the type species of the genus Paenibacillus), P. peoriae, P. kribbensis, P. ottowii, P. brasilensis, P. terrae, and "P. maysiensis" formed a monophyletic complex (Paenibacillus polymyxa complex) (Figure 1). Strains RP20, RP51, and RP62 isolated from nodules of Robinia pseudoacacia and strain CFCC 1854 isolated from Dendrolobium triangulare were classified to P. peoriae because their 16S rRNA gene sequences and that of P. peoriae type strain have identities above 99.6% and phylogenetically grouped together (Figure 1). Strain CFCC 1865 isolated from nodules of Ormosia semicastrata was classified to P. kribbensis because its 16S rRNA gene sequence and that of P. kribbensis type strain AM49 T have an identity of 99.7% and phylogenetically grouped together. Strain RP31 isolated from nodules of Robinia pseudo acacia and strain HKA-15 isolated from nodules of soybean (Annapurna et al., 2013) may belong to a same species because their 16S rRNA gene sequences showed identities of 99.2% and phylogenetically grouped together; they may belong to P. kribbensis or a novel species closely related to P. kribbensis because their 16S rRNA gene sequences showed identities about 98.8% to that of P. kribbensis AM49 T and phylogenetically grouped close to P. kribbensis AM49 T (Figure 1).
"Paenibacillus enshidis" CCTCC AB 2013275 T (=RP-207 T ) isolated from Robinia pseudoacacia (Yin et al., 2015) and P. medicaginis CC-Alfalfa-19 T isolated from nodules of alfalfa (Lai et al., 2015) are phylogenetic neighbors within a cluster closely related to the cluster containing the P. polymyxa complex (Figure 1). Strain CFCC 1991 isolated from nodules of Acacia implexa was classified to P. taichungensis because its 16S rRNA gene sequence and that of P. taichungensis type strain have an identity of 99.8% and phylogenetically grouped together within a standalone phylogenetic cluster (Figure 1). Strain CFCC 13938 isolated from nodules of Acacia crassicarpa was classified to P. puldeungensis because its partial 16S rRNA gene sequence (1395 bp) is identical to that of P. puldeungensis type strain CAU 9324 T within a standalone phylogenetic cluster (Figure 1).
Strain RP43 isolated from nodules of Robinia pseudoacacia may belong to P. anseongense or a novel species closely related to P. anseongense because its 16S rRNA gene sequence and that of P. anseongense type strain have an identity of 99.0% and phylogenetically grouped together (Figure 1). Strain RP43 is relatively close to P. periandrae PM10 T isolated from nodules of P. mediterranea (Menéndez et al., 2016).
The nodule-inhabiting strains within the P. polymyxa complex inhibited the growth of all the five tested fungal pathogens. P. peoriae strains RP20, RP51, and RP62, P. kribbensis CFCC 1865, and Paenibacillus sp. RP31 showed above 60% inhibition on the growth of the five fungal pathogens. P. peoriae CFCC 1854 showed relatively lower inhibition (below 60%) on F. graminearum, B. cinerea, R. solani, and S. sclerotiorum. RP20, RP51, RP62, RP31, and CFCC 1865 showed similar inhibition on M. oryzae. Strain RP51 showed the highest inhibition on F. graminearum, B. cinerea, and S. sclerotiorum. Strain RP31 showed the highest inhibition on B. cinerea and R. solani. In contrast, P. peoriae type strain CGMCC 1.3761 T , which was isolated from soil, inhibited only the growth of M. oryzae at a relatively lower extent than that by the six noduleinhabiting strains within the P. polymyxa complex ( Table 2 and Supplementary Figure 2). P. endophyticus CCTCC AB 2014195 T inhibited the growth of F. graminearum, M. oryzae, and R. solani but not B. cinerea and S. sclerotiorum. "P. enshidis" CCTCC AB 2013275 T and Paenibacillus sp. RP43 inhibited only the growth of F. graminearum. P. puldeungensis CFCC 13938 and P. taichungensis CFCC 1991 did not inhibit fugal growth ( Table 2).

Nodule-Inhabiting Strains Within Paenibacillus polymyxa Complex Produced Hydrolytic Enzymes, Siderophores, and Fusaricidins
Six nodule-inhabiting Paenibacillus strains within the P. polymyxa complex, P. peoriae RP20, RP51, RP62, and CFCC 1854, P. kribbensis CFCC 1865, and Paenibacillus sp. RP31, showed antagonistic activities against all the five tested fungal pathogens and thus were noted as Paenibacillus against broadspectrum fungi. To know the mechanisms of the broad-spectrum antifungal activities of these nodule-inhabiting strains within the P. polymyxa complex, some known antifungal substances (hydrolytic enzymes, siderophores, and fusaricidins) were examined. P. peoriae strains CGMCC 1.3761 T within the P. polymyxa complex was also examined for comparison. All the seven strains within the P. polymyxa complex showed activities of chitinase, β-1,3-glucanase, and protease except that P. peoriae RP51 showing the highest inhibition on F. graminearum, B. cinerea, and S. sclerotiorum did not show protease activity (Figure 2). They show similar chitinase activities. P. peoriae CGMCC 1.3761 T and RP51, and P. kribbensis CFCC 1865 showed relatively stronger β-1,3-glucanase activities (Figure 2). The weak antifungal strain P. peoriae CGMCC 1.3761 T showed similar chitinase activity and β-1,3-glucanase activity as the potent antifungal strain RP51. Seemingly, the hydrolytic enzyme activities detected from these strains were not consistent with the extent of their antifungal activities. Therefore, the hydrolytic enzymes may not contribute to the broad-spectrum antifungal activities of the nodule-inhabiting strains within the P. polymyxa complex.
All the seven strains grown on the iron limited medium produced siderophores detected by the chrome azurol S assay (Figure 3). P. peoriae RP20, RP51, and RP62, and Paenibacillus sp. RP31 shown potent antifungal activities on PDA appeared to produce relatively more siderophores (Figure 3). However, the production of siderophores may not contribute to the broad-spectrum antifungal activities detected on the ironsufficient PDA.
The variety and quantity of the fusaricidins produced by the strains within the P. polymyxa complex are generally consistent with the extent of their antifungal activities. P. peoriae CGMCC 1.3761 T and CFCC 1854 produced fewer variety and amount of fusaricidins and showed weaker antifungal activities comparing with the other five strains. Therefore, fusaricidins are likely responsible for the antifungal activities against the broad-spectrum phytopathogenic fungi.

Most Strains Within Paenibacillus polymyxa Complex Contain the fus Gene Cluster
In the whole genome sequences of 396 strains assigned to 165 Paenibacillus species and deposited in the NCBI database, the gene cluster for the biosynthesis of fusaricidins (fus) was detected from 88 strains (22% of 396 strains) using antiSMASH. The fus cluster appeared to be aggregated in closely related species. In particular, the fus cluster was aggregated in 73 (96%) of 76 strains FIGURE 2 | Agar plate assays show hydrolysis of chitin, β-1,3-glucan, and casein by chitinase, β-1,3-glucanase, and protease from P. peoriae CGMCC 1.3761 T and nodule-inhabiting strains within Paenibacillus polymyxa complex.
FIGURE 3 | Chrome azurol S assay shows color changes from blue to light purple or orange indicating production of catechol type or hydroxamate type siderophores from P. peoriae CGMCC 1.3761 T and nodule-inhabiting strains within Paenibacillus polymyxa complex. within the P. polymyxa complex and in 8 of 10 strains within the P. elgii complex including P. elgii, P. tyrfis, and P. tianmuensis (Figure 1 and Supplementary Table 1). Gene clusters for biosynthesis of siderophores were detected from 99 strains (25% of 396 strains) and from 25 (33%) of 76 strains within the P. polymyxa complex (Supplementary Table 1). The genomes of nodule-inhabiting P. lupini type strain and P. prosopidis type strain out of the P. polymyxa complex (Figure 1)

Fusaricidin-Producing Paenibacillus Damaged Fungal Cell Structures
Paenibacillus peoriae RP51 showing potent broad-spectrum antifungal activities and fusaricidin-producing activity was used   Kajimura and Kaneda (1997); Vater et al. (2015) to detect the antifungal action on fungal cell structures. Control F. graminearum hyphae grown on PDA were smooth and intact shown by SEM ( Figure 4A); the hyphal cells contained electrondense cell contents and electron-lucent lipid granules shown by TEM ( Figure 4C). In contrast, the hypha of F. graminearum confronting RP51 were shrunk and distorted ( Figure 4B); hyphal cell walls were broken while cell contents were leaked indicating by vacuolization (Figure 4D).

Nodule-Inhabiting Strains Within
Paenibacillus polymyxa Complex Can Produce IAA and Dissolve Phosphate or Fix N 2 The six nodule-inhabiting strains within the P. polymyxa complex and P. peoriae CGMCC 1.3761 T from soil were able to produce IAA about 5 µg from cultures at 1 × 10 9 CFU ml −1 . They were also able to dissolve calcium phosphate rich in alkaline soils, and calcium phytate in soils rich in organic phosphate (Figure 6) but not ferric phosphate rich in acidic soils. P. peoriae RP20, RP51, and RP62 were positive for nifH gene amplification (Figure 7) and could grow well on the nitrogen-free Jensen's agar medium. These nodule-inhabiting strains may stimulate the seed germination and seedling growth by producing IAA. P. peoriae RP20, RP51, and RP62 may fix N 2 in association with seedlings in the gnotobiotic culture without external nutrient supply but their phosphate solubilization activity was not involved in their FIGURE 5 | Colonization of wheat seedlings by Paenibacillus peoriae CGMCC 1.3761 T and nodule-inhabiting strains within P. polymyxa complex. Bacteria were recovered at 9 days after inoculation to imbibed wheat seeds. The different letters on the columns and standard error bars indicate significant difference between the treatments at p < 0.05. P. peoriae RP51 and Paenibacillus sp.

RP31 Can Control Fungal Infection in planta
Paenibacillus peoriae RP51 and Paenibacillus sp. RP31 showed potent antifungal activities in vitro, potent plant colonization and plant growth-promoting competence. They were selected to evaluate their potential in biocontrol of fungal diseases in planta using the pathosystem of FRR and FFR on wheat seedlings (Colombo et al., 2019). After inoculation of a Fusarium plug on a primary root at 10 mm distance from the seed, necrosis of plant tissues occurred along the roots (root rot) and shoots (foot rot) of wheat seedlings (Supplementary Figure 1). At 4 days after fungal inoculation, strain RP51 and strain RP31, which were inoculated to wheat seeds, inhibited the FRR symptom about 66% and 60%, respectively ( Table 5). At 6 days after fungal inoculation, strain RP51 and strain RP31 inhibited the FFR symptom about 61% and 55%, respectively (Table 5).
Our MALDI-TOF-MS analysis identified multiple fusaricidins produced by the nodule-inhabiting strains within the P. polymyxa complex. The variety and quantity of the fusaricidins produced by these strains within the P. polymyxa complex are generally consistent with the extent of their antifungal activities. Likely, fusaricidins are responsible for the antifungal activities against broad-spectrum phytopathogenic fungi by the nodule-inhabiting strains within the P. polymyxa complex.
Fusaricidins and the single fus operon encoding a nonribosomal peptide synthetase (fusA) for the biosynthesis of a variety of fusaricidins (Han et al., 2012) have been demonstrated to play roles in P. polymyxa against broad-spectrum fungal pathogens (Choi et al., 2008;Li et al., 2013;. Fusaricidins produced by P. polymyxa have been demonstrated to inhibit Fusarium spore germination and disrupt hyphal membranes . Likely, fusaricidins are responsible for the damage of Fusarium hyphal cells by P. peoriae RP51 found in this study. Notably, our genome mining revealed that 83% (73 of 88) fus-containing strains in the genus Paenibacillus belong to the P. polymyxa complex and 96% (73 of 76) strains within the P. polymyxa complex contain the fus cluster. Consistently, Liu et al. (2019) showed that all isolated nitrogen-fixing strains within the P. polymyxa complex (n = 20) inhibited the growth of four or over four of six tested phytopathogenic fungi whereas other five strains outside of the P. polymyxa complex did not inhibit any of the six tested fungi. Fusaricidins produced by P. polymyxa have also been demonstrated to induce plant systemic resistance via salicylic acid pathway against Fusarium and Phytophthora pathogens (Lee et al., 2013;. Therefore, the P. polymyxa complex is a promising resource for fusaricidin-dependent control of broad-spectrum fungal pathogens and inducing plant systemic resistance. Our genome mining also revealed that all strains (n = 76) within the P. polymyxa complex contain the ilvB (=ipdC) gene and 43% of them (33 of 76) contain the nif cluster. Consistently, the six nodule-inhabiting strains within the P. polymyxa complex can produce IAA and half can fix N 2 . Likewise, Quyet-Tien et al.
(2010) found all tested P. polymyxa strains (n = 29) producing IAA and 38% of them (11 of 29) containing the nifH gene. Liu et al. (2019) showed that almost all nitrogen-fixing strains (19 of 20) within the P. polymyxa complex produced IAA. The nodule-inhabiting strain HKA-15 within the P. polymyxa complex also contains the nifH gene (Senthilkumar et al., 2009). Anand et al. (2013) and Padda et al. (2016) did 15 N dilution assay and revealed that P. polymyxa promoted plant growth via nitrogen fixation. In contrast, only five of 5216 strains assigned to the genus Bacillus contain the nif clusters, which are most closely related to those of Paenibacillus. Most likely, nitrogenfixing strains within the P. polymyxa complex have the advantage of nitrogen fixation over Bacillus strains as biocontrol agents or biofertilizers.
Together, production of fusaricidins and IAA and nitrogen fixation are likely conserved in the P. polymyxa complex and evolved in their adaptations to the plant-associated life (Xie et al., 2016;Wang et al., 2020;Zhou et al., 2020). Our genome mining and screening of Paenibacillus strains revealed that the strains within the P. polymyxa complex have a high probability to be plant growth promoters and biocontrol agents against broad-spectrum phytopathogenic fungi. However, there are exceptions for the strains within the P. polymyxa complex to have broad-spectrum antifungal activities (e.g., P. peoriae type strain) or promote plant growth (e.g., P. peoriae CFCC 1854). Therefore, screening of the strains within the P. polymyxa complex for biocontrol and biofertilization is necessary. We screened out P. peoriae RP51 and Paenibacillus sp. RP31 showing broad-spectrum antifungal activities, potent plant colonization competence and plant growth-promoting activities, and effective in planta control of fungal infection. Moreover, we propose an effective approach to screen Paenibacillus strains to be effective biocontrol agents and biofertilizers based on our study and other studies (e. g. Quyet-Tien et al., 2010;Colombo et al., 2019;Liu et al., 2019). First, phylogenetic analysis of nearly complete 16S rRNA gene sequences identifies strains within the P. polymyxa complex, which most likely produce fusaricidins and IAA. Second, confrontation culture of Paenibacillus and pathogens and amplification of nifH gene screen out Paenibacillus strains having broad-spectrum antifungal activity and nitrogen-fixing ability. Third, gnotobiotic cultures of plant seedlings, Paenibacillus, and pathogens screen out Paenibacillus strains having plant colonization competence and in planta biocontrol potential. Fourth, field tests on the promising Paenibacillus strains, such as P. peoriae RP51, determine effective biocontrol agents and biofertilizers, leading to reduced use of agrochemicals for sustainable agriculture.
Conclusively, our study highlights that endophytic strains within the P. polymyxa complex have a high probability to be effective biocontrol agents and biofertilizers.

DATA AVAILABILITY STATEMENT
The dataset (Supplementary Table 1) generated for this study can be found in the online repository. The nucleotide sequences generated in this study were deposited at: https://www.ncbi. nlm.nih.gov/genbank/ under the accession numbers MN715870, MN715872, MN715875, MT093458, MN715871, MT093459, MT093460, MT093461, and MT040706.

AUTHOR CONTRIBUTIONS
MA, LL, XL, YY, and QA contributed to the conceptualization. LL, YY, and QA contributed to the microbial resources. MA, YL, RH, AH, and TX contributed to the investigation. MA and QA contributed to the data analysis. MA and QA contributed to the manuscript preparation. XL, BL, JY, and QA contributed to the review and editing. BL, JY, and QA contributed to the supervision. LL, BL, JY, and QA contributed to the funding acquisition. All authors contributed to the article and approved the submitted version.