Potential Role of Rhizobacteria Isolated from Citrus Rhizosphere for Biological Control of Citrus Dry Root Rot

Citrus trees face threats from several diseases that affect its production, in particular dry root rot (DRR). DRR is a multifactorial disease mainly attributed to Neocosmospora (Fusarium) solani and other several species of Neocosmospora and Fusarium spp. Nowadays, biological control holds a promising control strategy that showed its great potential as a reliable eco-friendly method for managing DRR disease. In the present study, antagonist rhizobacteria isolates were screened based on in vitro dual culture bioassay with N. solani. Out of 210 bacterial isolates collected from citrus rhizosphere, twenty isolates were selected and identified to the species level based on the 16S rRNA gene. Molecular identification based on 16S rRNA gene revealed nine species belonging to Bacillus, Stenotrophomonas, and Sphingobacterium genus. In addition, their possible mechanisms involved in biocontrol and plant growth promoting traits were also investigated. Results showed that pectinase, cellulose, and chitinase were produced by eighteen, sixteen, and eight bacterial isolates, respectively. All twenty isolates were able to produce amylase and protease, only four isolates produced hydrogen cyanide, fourteen isolates have solubilized tricalcium phosphate, and ten had the ability to produce indole-3-acetic acid (IAA). Surprisingly, antagonist bacteria differed substantially in their ability to produce antimicrobial substances such as bacillomycin (five isolates), iturin (ten isolates), fengycin (six isolates), surfactin (fourteen isolates), and bacteriocin (subtilosin A (six isolates)). Regarding the PGPR capabilities, an increase in the growth of the bacterial treated canola plants, used as a model plant, was observed. Interestingly, both bacterial isolates Bacillus subtilis K4-4 and GH3-8 appear to be more promising as biocontrol agents, since they completely suppressed the disease in greenhouse trials. Moreover, these antagonist bacteria could be used as bio-fertilizer for sustainable agriculture.


Introduction
Citrus is an important economic crop in Morocco, covering an area of 126.600 ha with a production of 2.6 million tons of citrus fruits a year. The annual exportation of citrus fruits to Europe and other western countries has reached 755.000 tons [1,2]. However, citrus plantations are threatened by several pathogens of which the Fusarium species, mainly F. solani, is causing citrus dry root rot (DRR) disease, which is responsible for significant economic loss [1,3,4]. In addition, Fusarium species can cause other serious diseases on citrus plantation such as twig rot, decline dieback, twig blight, and vascular wilt, thereby are a major threat to citrus production worldwide [1,5,6]. Fusarium sensu latu was recently Studies on the biological control of citrus DRR disease are almost absent worldwide, which might urge the exploration of the core microbiome of citrus rhizosphere for this purpose. The rhizosphere is an environment with high microbial diversity, and it is known to be a reservoir of microorganism that possess some important biological activities [37]. In addition, microorganisms that are present in the rhizosphere of healthy citrus trees may have the possibility to compete for nutrients, improve plant growth, and improve their abilities to resist to different biotic and abiotic stress. Therefore, the present study aimed at exploring the potential biocontrol of bacterial isolates, from the citrus rhizosphere of different Moroccan regions, against Neocosmospora (Fusarium) solani, the most dominant species associated with DRR disease. Furthermore, the most effective bacterial isolates selected from dual culture bioassay will be characterized to the species level using 16S rDNA gene. These bacteria will be screened for (i) biochemical traits and genes involved in production of antifungal compounds, (ii) their PGPR capabilities with canola (Brassica napus), and (iii) in planta greenhouse trials to control DRR disease.

Antagonistic Activity
A total of 210 morphologically different bacterial colonies isolated from different citrus rhizosphere areas across the country were tested for their antagonistic activity against N. solani using dual culture bioassay ( Figure 1). The zone of inhibition offered by each bacterium was assessed 7 days post-incubation period (Table 1). Results underlined that twenty isolates displayed an important antifungal activity against N. solani (>60%) with very highly significant difference. The selected antagonistic bacteria were also tested for their antifungal activity against Fusarium oxysporum, F. brachygibbosum, and F. equiseti ( Table 1). The most promising bacterial isolates were k4-4, TD1, and GH3-8 with inhibition rates of mycelial growth of 72.22, 71.50, and 68.93%, respectively. main families of lipopeptides [35]. The bacteriocins, including subtilosin A, are ribosomally produced, with potent antimicrobial peptides, by diverse groups of micro-organisms able to form pores in the membrane [35,36]. Studies on the biological control of citrus DRR disease are almost absent worldwide, which might urge the exploration of the core microbiome of citrus rhizosphere for this purpose. The rhizosphere is an environment with high microbial diversity, and it is known to be a reservoir of microorganism that possess some important biological activities [37]. In addition, microorganisms that are present in the rhizosphere of healthy citrus trees may have the possibility to compete for nutrients, improve plant growth, and improve their abilities to resist to different biotic and abiotic stress. Therefore, the present study aimed at exploring the potential biocontrol of bacterial isolates, from the citrus rhizosphere of different Moroccan regions, against Neocosmospora (Fusarium) solani, the most dominant species associated with DRR disease. Furthermore, the most effective bacterial isolates selected from dual culture bioassay will be characterized to the species level using 16S rDNA gene. These bacteria will be screened for (i) biochemical traits and genes involved in production of antifungal compounds, (ii) their PGPR capabilities with canola (Brassica napus), and (iii) in planta greenhouse trials to control DRR disease.

Antagonistic Activity
A total of 210 morphologically different bacterial colonies isolated from different citrus rhizosphere areas across the country were tested for their antagonistic activity against N. solani using dual culture bioassay ( Figure 1). The zone of inhibition offered by each bacterium was assessed 7 days post-incubation period (Table 1). Results underlined that twenty isolates displayed an important antifungal activity against N. solani (>60%) with very highly significant difference. The selected antagonistic bacteria were also tested for their antifungal activity against Fusarium oxysporum, F. brachygibbosum, and F. equiseti (Table 1). The most promising bacterial isolates were k4-4, TD1, and GH3-8 with inhibition rates of mycelial growth of 72.22, 71.50, and 68.93%, respectively.

Effect of Bacterial Isolates on the Mycelial Structure of Neocosmospora Solani
Microscopic observations of the mycelium of N. solani in co-culture with the antagonist bacteria displayed damaged morphology and distinct cytological alterations when compared to the untreated control ( Figure 2A). In most cases, the alterations correspond to hyphal swelling, deformation, and vacuolation of the mycelial structure ( Figure 2B-D), and massive conglobated along hyphae with uneven surface ( Figure 2E,F) and sometimes associated with a degradation of the mycelium and a release of the cytoplasmic contents ( Figure 2E).

Biochemical Traits
All isolates were tested for their ability to produce cell wall-degrading enzymes (Figures S1 and S2; Table 3).
Amylase Activity It was seen that all bacterial isolates were able to produce the amylase ( Table 3). The three bacterial isolates TD7, BM3-5, and GH1-5 showed a high amylolytic index of 2.10, 2.07, and 1.98, respectively.
Cellulase Degradation Results showed that among twenty bacterial isolates tested, sixteen had the ability to produce the cellulase ( Table 3). The three bacterial isolates TD7, BM4-3, and TG5 displayed higher cellulolytic index of 2.62, 2.62, and 2.43, respectively.

Phosphate Solubilization
The phosphate solubilization test showed that among twenty rhizobacterial isolates tested, fourteen had the ability to produce clear zones (clear halo > 1 mm) of phosphate solubilization on PVK agar medium around the colony on the plate after four days of  (Table 3). The highest solubilization index was obtained with three bacterial isolates GH3-8 (2.19), TD1 (2.15), and Bel3-4 (2.14).

Indole-3-acetic Acid Production
The results of IAA production were found positive for ten of twenty bacterial isolates and for which the supernatant culture changed to red color after the addition of the Salkowski's reagent.
Siderophores Production Results evinced that only three bacterial isolates k4-4, Bel3-4, and GH1-2 were able to produce siderophores. The change in the color of CAS from blue to orange or purple explained by the transfer of ferric ions from it to the siderophores.

Hypersensitivity Test
By comparing the tobacco leaves with the positive control (leaf injected with the bacterial suspension of Agrobacterium tumefaciens), our results showed that no bacterial isolates have induced a hypersensitivity reaction.

Effect of Rhizobacteria on Plant Growth of Brassica Napus
The selected bacteria were evaluated subsequently for their ability to form a beneficial association with canola by promoting its growth. Statistical analysis showed highly significant difference (p < 0.001) between treated and untreated canola seedlings with bacterial isolates for the plant length, root length, fresh root, and dry root weight. Interestingly, significant difference was also observed between antagonist bacteria. However, no significant difference was observed for the number of leaves. Our results showed that some rhizobacteria has a positive effect on the plant growth ( Figure S3). In fact, canola seedlings treated with BM3-2, K3-7, TG5, B2-1, and K4-4 exhibited the highest increase in plant length, while those with B2-1, TD7, BM3-4, and bel3-4 showed an increase in root length. However, plants treated with TD7, BM3-2, B2-1, and Bel3-4 displayed a significant increase in total fresh weight, while those with B2-1, TM10, bel3-4, TD7 expressed a substantial increase in total dry weight. For fresh and dry root weight, only plants treated with bacterial isolates BM4-1, TD7, GH1-2, and K4-4 have registered a significant increase compared to untreated control plants ( Figure S3; Table 4).

Hypersensitivity Test
By comparing the tobacco leaves with the positive control (leaf injected with the bacterial suspension of Agrobacterium tumefaciens), our results showed that no bacterial isolates have induced a hypersensitivity reaction.

Effect of Rhizobacteria on Plant Growth of Brassica Napus
The selected bacteria were evaluated subsequently for their ability to form a beneficial association with canola by promoting its growth. Statistical analysis showed highly significant difference (p < 0.001) between treated and untreated canola seedlings with bacterial isolates for the plant length, root length, fresh root, and dry root weight. Interestingly, significant difference was also observed between antagonist bacteria. However, no significant difference was observed for the number of leaves. Our results showed that some rhizobacteria has a positive effect on the plant growth ( Figure S3). In fact, canola seedlings treated with BM3-2, K3-7, TG5, B2-1, and K4-4 exhibited the highest increase in plant length, while those with B2-1, TD7, BM3-4, and bel3-4 showed an increase in root length. However, plants treated with TD7, BM3-2, B2-1, and Bel3-4 displayed a significant increase in total fresh weight, while those with B2-1, TM10, bel3-4, TD7 expressed a substantial increase in total dry weight. For fresh and dry root weight, only plants treated with bacterial isolates BM4-1, TD7, GH1-2, and K4-4 have registered a significant increase compared to untreated control plants ( Figure S3; Table 4).

In Planta Bioassays
Based on the in vitro bioassay, ten antagonist bacteria were subjected to greenhouse trials to confirm their antifungal activity and their ability to suppress the citrus dry root rot disease. After eight weeks of post-inoculation periods, results indicated that bacterial isolate K4-4, GH3-8, and K4-4 were highly effective in suppressing the DRR disease (absence of disease symptoms). These bacterial isolates gave comparable results to that obtained with B. subtilis commercial product and negative control (C−) (without N. solani) (Figures 5 and 6). Interestingly, bacterial isolates K4-3, TD1, TG6, B2-1, and BM3-2 significantly reduced the disease severity without suppressing it completely. However, bothbacterial isolates TG5 and GH1-5 were shown to be less effective in controlling the disease ( Figure 6).

Discussion
In this study, we highlight that the plant rhizosphere is a valuable source of potent rhizobacteria that may serve as an eco-friendly solution for the control of soilborne diseases. In this context, a collection of 210 rhizobacteria were recovered from the citrus rhizospherre and screened for their inhibitory effect based on the in vitro dual culture bioassay between these bacterial isolates and N. solani. Twenty selected isolates have shown an important antifungal activity against N. solani. These bacterial isolates were also tested on other three other Fusarium species previously isolated and reported to be associated with the DRR disease [1]. Furthermore, microscopic examinations of N. solani mycelium in the presence of antagonist bacteria revealed the presence of structural changes including deformations, swelling, and vacuolation of the mycelium sometimes accompanied by degradation of the mycelium and release of cytoplasmic contents. The inhibition of the growth of the pathogenic fungus and damages observed on its hyphae are probably due to secretion of hydrolytic enzymes and lipopeptides by antagonist bacteria [38,39].
Previous studies have revealed the importance of antagonist bacteria in the control of Fusarium spp. [22,24,25,[40][41][42][43]. In our study, effective bacterial isolates were identified on the basis of the partial 16S rDNA genes. Three distinct genus were distinguished namely Bacillus, Stenotrophomonas and Sphingobacterium. The association of different species with citrus rhizosphere is linked to the diversity found in the citrus root exudates that may attract these species and favour their growth in order to colonise citrus root tissues [24,44]. Over time, several bacteria are emerging as novel PGPR and declared to have biological control potentials against different fungal plant diseases. Bacillus species were reported in several studies to have an important potential as BCAs [45][46][47]. As our study highlighted several species belonging to the same genus, Ali et al. [47] also reported the great ability of Bacillus tequilensis S5 to reduce the mycelial growth by 76.6%. Indeed, the cell-free filtrates of this bacterium lowered the growth of the pathogenic fungus by 82.2%. Chenniappan et al. [44] selected sixteen bacterial isolates capable of reducing the growth of several fungal pathogens of the Neocosmospora and Fusarium genus, in particular N. solani. These bacterial isolates were identified as Bacillus amyloliquefaciens, B. tequilensis, and B. subtilis. Authors have also identified the same genes involved in the mechanisms of biocontrol as we reported in our study. B. xiamenensis was recently shown to have antagonistic and PGP activities and could suppress red rot disease and enhance sugarcane growth [48]. Similarly, B. licheniformis was previously reported as a promising BCAs against plant pathogens [49][50][51]. B. velezensis was previously identified to the B. amyloliquefaciens group [52]. Damasceno et al. [53] underscored the high performance of B. velezensis to reduce the mycelial growth of Colletotrichum musae. This bacterium was able to mitigate the same performance as the fungicidal product thiabendazole. However, Rojas-Solís et al. [54] underlined a higher antagonist activity of Stenotrophomonas maltophilia CR71 against Botrytis cinerea, which was probably realted to its ability to emit volatile organic compounds (VOCs). Moreover, this bacterium was able to promote growth and achieve an effective biocontrol of B. cinerea through the production of potent volatiles such as DMDS. It was found that S. maltophilia lowered the mycelial growth of Colletotrichum nymphaeae by 60% [55]. For the antagonist Sphingobacterium multivorum, few studies have demonstrated its capabilities as BCAs [56]. Surprisingly, our study found that S. multivorum inhibited the mycelial growth of N. solani and reduced the DRR severity. This bacterium was previously reported to have important antifungal effect against Magnaporthe oryzae, the causal agent of rice blast [56]. This bacterial species was also found effective in degrading hexaconazole, thereby presenting a sustainable microbial bioremediation of persistent organic pollutants [57]. Bacillus spp. can produce several lytic enzymes and synthesise a wide range of metabolites including biosurfactants such as lipopeptides that have the ability to be effective in controlling different plant diseases [58]. The antibiosis via bacterial supernatant indicated that the twenty rhizobacterial isolates might be a pool of various secondary metabolites. Similar results were reported by Li et al. [38] who confirmed that cell-free filtrate from Bacillus megaterium inhibited the mycelial growth of Alternaria alternata. Bacterial antagonism by antibiosis appears to be the main mechanism involved in the biological control of plant pathogenic fungi [40], in addition to the production volatile antifungal compounds (VOCs) [59]. In this study, our findings indicated that effective bacterial isolates were able to produce VOCs and inhibit mycelial growth of N. solani. Similarly, Guevara-Avendaño et al. [60] highlighted that organic compounds (VOCs) emitted by B. amyloliquefaciens were found to inhibit N. solani mycelial growth and induced slight distortions of fungal hyphae. VOCs are chemical substances that can be easily evaporated into the air due to their low molecular weight, high vapor pressure, and low water solubility [61,62].
In this study, the effectiveness of selected bacterial isolates to control the disease in planta and their performance to promote plant growth were assessed. This feature is highly suitable for BCAs and helps to suppress the disease and promote plant growth and crops productivity. In this study, the above-mentioned mechanisms were considered as an important criterion for the selection of effective PGPR. Additionally, our results emphasized that selected antagonist rhizobacteria with high antifungal activity shared excellent attributes. Most of these bacterial isolates exerted different biocontrol mechanisms such as production of cell wall degrading enzymes, production of hydrocyanic acid (HCN), along with plant growth promotion traits. Production of lytic enzymes is among the major mechanisms employed by biocontrol agents to control fungal pathogens [40,44]. Our results showed that of twenty antagonist bacterial, eighteen displayed pectinase production, while cellulase and chitinase production was observed for sixteen and eight bacteria, respectively. Interestingly, all bacteria were able to produce amylase and protease. chitinases attack fungal cell wall and cause lysis by degrading chitin [40]. In addition to hydrolytic enzymes, PGPR produce chemical compounds with different benefits for the plant. Furthermore, HCN is an antimicrobial compound synthesized from glycine using HCN synthase, an enzyme encoded by a set of three genes (hcnA, hcnB, and hcnC) [63] and it is known to be deleterious to microorganisms [40]. According to Blume et al. [64], this metabolite acts on the cells of phytopathogenic fungi by blocking cytochrome oxidase in the respiratory chain. Our results indicated that four effective bacterial isolates produced HCN. Similar results were found for Pseudomonas fluorescens, which produce siderophores, HCN, and chitinase, and by which the antagonist more likely impedes the growth of Sclerotinia sclerotiorum on Brassica campestris [65]. Furthermore, promising bacterial isolates were tested in vitro for properties that are known to be important for promoting plant growth, such as solubilization of phosphate, the production of siderophores, and IAA production [23]. Siderophores were characterized by a strong affinity for iron and is one of the important mechanisms used by BCAs [66,67]. Among the twenty rhizobacterial antagonists tested in this study, three isolates were siderophores producers.
The phosphates solubilizing microorganisms (PSM) are considered as potential bioinoculants to increase crop productivity by transforming phosphorus of the soil (the monobasic (H 2 PO 4 − ) and diabetic (HPO 4 2− ) into soluble forms assimilable by plants, through the production of organic acids, the acidification process and chelation [68]. In this study, 14 isolates were found able to solubilize tricalcium phosphate. Sun et al. [69] found that out of nineteen bacteria, thirteen were able to solubilize phosphate and ten produced IAA. Phytohormones are known to play a key role in regulating the growth of all parts of the plant during its various development stages. PGPRs can produce phytohormones, similar to those produced by plants, which enhance its growth. Additionally, they are important in generating plant defence responses against invading pathogens [41]. IAA produced by PGPRs contribute in increasing root surface area and length, and thereby better soil nutrients uptake [42,68]. Interestingly, half of effective bacterial isolates produced different shades of pinkish-red, suggesting the production of this phytohormone. Results from this study demonstrated that rhizobacteria could be effective in controlling Neocosmospora and Fusarium spp. and promote plant growth, which is in agreement with other previous studies. Our finding corroborated those of Majeed et al. [70] who found that Stenotrophomonas and Bacillus produced IAA and solubilize phosphate. A Bacillus sp. strain was found to strongly inhibit the in vitro growth of several phytopathogens including N. solani and F. oxysporum [71]. This bacterial isolate possessed several biocontrol and PGP mechanisms such as production of lytic enzymes, IAA, and siderophores, as well as the ability to solubilize various sources of organic and inorganic phosphates [71]. The results appeared to be the same as the report of Paramanandham et al. [25] who underscored that Pseudomonas aeruginosa produced IAA, HCN, siderophores, chitinase, proteases, and cellulase which played a crucial role in the antifungal activity against F. oxysporum f. sp. lycopersici, and A. solani on tomatoes. This result was similar to the findings of Chenniappan et al. [44] who isolated 16 antagonist bacteria from turmeric rhizosphere with the ability to produce lytic enzymes and having PGPR attributes. Our findings corroborated those of Slama et al. [22] who confirmed that B. halotolerans produces lytic enzymes; amylase, protease, cellulase, and chitinase, auxin, siderophores production, and phosphate solubilization that can used as BCA against Bayoud disease caused by F. oxysporum f. sp. albedinis and an efficient bio-fertilizer of oasis ecosystems. In a similar study, Palmieri et al. [41] screened four antagonist bacteria Rhanella aquatilis, Pseudomonas fluorescens, Serratia marcescens, and B. amyloliquefaciens against Foc and Fsp of wilt and root rot causing chickpea and have proved their complementary capacities to solubilize tricalcium phosphate, chitinase production, and IAA production. Our result was also in accordance with those reported by Prajakta et al. [72] who reported that Bacillus mojavensis PB-35(R11) exhibited phosphate solubilization, hydrogen cyanide, chitinase, IAA, and siderophores production, which may be involved in antifungal activity against Rhizoctonia solani and in promoting plant growth of soybean. Lim et al. [73] recently reported a B. velezensis strain capable of producing antifungal volatile and agar-diffusible metabolites that inhibited mycelial growth of several plant pathogenic fungi. Similarly, Martínez-Raudales et al. [74] evidenced an antifungal activity of B. velezensis against N. solani, F. oxysporum, Phytophthora capsici, and R. solani, causative agents of chili pepper root rot.
In this study, the biosynthesis genes encoding for the production of lipopeptides (bacillomycin, iturin, fengycin, and surfactin) and bacteriocin (subtilosin) were investigated. Surprisingly, our selected antagonist bacteria harbour biocontrol genes responsible for the production of the lipopeptides and bacteriocin. These molecules are capable of decreasing pathogen growth [75]. Therefore, our study provided valid proof confirming that genes encoding for lipopeptides played an important role in fighting against N. solani. The ability of bacteria to produce lipopeptides is crucial in evaluating its potential as a BCAS against plant pathogens [76]. Lipopeptides from the fengycin, surfactin, and iturin families have been shown in several studies to have considerable potential for controlling plant pathogenic fungi [44,76]. Zhang et al. [43] reported the presence of genes encoding for lipopeptides biosynthesis in most of the Bacillus isolates tested. Fengycins have fungitoxic properties, particularly against filamentous fungi [77]. Subtilosin A was reported to be produced by B. amyloliquefaciens [36]. PGPR produce a variety of antibiotic compounds including lipopeptides, which are considered as major contributors to Bacillus antifungal activity [78]. Cao et al. [78] demonstrated that iturin and fengycin secreted by B. velezensis are responsible for its antimicrobial activity, while surfactin proposed to contribute in the formation of biofilms and cell motility, both of which are important for successful rhizosphere colonization. Bacillus velezensis was found capable to produce a variety of antibiotics compounds including surfactin, iturin, fengycin, ericin, and others [79,80]. Zalila-Kolsi et al. [81] denoted that B. amyloliquefaciens and B. subtilis that produce (iturin and surfactin), (surfactin and fengycin), respectively, have a wide range of action against several phytopathogenic fungi according to Gong et al. [39], both iturin A and plipastatin (fengycine) A have fungicidal activity, but iturin A is active at lower concentrations than plipastatin A. In addition, treatments with both molecules showed a deformed and damaged hyphae morphology [39]. Tora et al. [33] demonstrated that the biocontrol of Botrytis cinerea by Bacillus XT1 was facilitated by lipopeptides, suggesting that mycelial structure of B. cinerea was probably degenerated due to these compounds.
The present study highlighted that most effective bacterial isolates were able to exhibit more than PGP traits. Our results underlined an increase in the plant growth of B. napus seedlings treated with antagonist bacteria (bacterization) in comparison with untreated controls. Previous studies documented the existence of multiple biocontrol mechanisms among the studied bacteria that explains their potential as successful biocontrol [44,82]. Sun et al. [69] reported that bacteria isolated from the rhizosphere enhance the plant growth of B. napus, which was also confirmed by the findings of Syed-Ab-Rahman et al. [83] who found that the bacterial isolates tested have contributed to a significant increase in lettuce seedling length, with the largest increase being observed for the plant treated B. amyloliquefaciens.
The reduction of disease severity by antagonist bacteria is more likely to be linked to their mechanisms of biocontrol and to their adaptation in the host plant environment. It has been proven that B. subtilis has a suppressive effect on plant pathogenic fungi under both in vitro and in planta conditions [84]. Treatment of cassava with B. amyloliquefaciens decreased the disease incidence of Fusarium root rot by more than 50% [85]. Additionally, this bacterium was found significant in suppressing the incidence of Fusarium wilt in tomato caused by Fusarium oxysporum f. sp. lycopersici [85]. Parikh et al. [86] underscored that Bacillus spp. suppressed the mycelial growth of pathogenic Fusarium isolates and decreased Fusarium root rot disease in corn, soybean, and wheat. These results were explained by the ability of Bacillus spp. to produce volatile organic compounds (VOCs) resulting in mycelial growth reduction of Fusarium spp. and causing morphological alterations in fungal hyphae [59]. Abd-Elgawad et al. [87] highlighted that the application of B. subtilis and P. fluorescens against Fusarium dry root rot of citrus reduced the incidence of the disease in treated trees and increased the yield. Application of B. amyloliquefaciens and Microbacterium imperiale against Fusarium root rot in cassava resulted in a reduction of disease incidence of more than 50% in greenhouse trials [85]. This is more likely due to the ability of this bacterium to sporulate, in addition to its higher population stability, which facilitates its storage, encapsulation, and subsequent application in the field. Zhang et al. [43] emphasized that the antagonistic and plant-growth promotion activities of bacterial isolates; P. fluorescens, Pseudomonas sp., B. subtilis, and Paenibacillus polymyxa might be related to their production of several types of lytic enzymes (β-1,3-glucanases, chitinases, and proteases), antibiotics, phytohormones (IAA), secondary metabolites, HCN, siderophores, and VOCs, and their abilities to solubilize the phosphate as well.
Practical output of our study suggests that B. subtilis K4-4 and Gh3-8 provided best control of DRR and could be used as BCAs due to their antifungal activity and plant growth promoting traits. Therefore, our finding represents an added value for biological control of plant diseases using antagonist bacteria; the selected bacterial isolates have a great biocontrol potential due to a higher ability of inhibiting the fungal growth and reducing the disease. Moreover, Bacillus species were the basis of the BCAs framework with numerous reports of their potential. PGPRs are widely used since they represent promoting traits for plant growth and they did not represent any risks to human health [88]. Our results represent an important contribution for better understanding of the citrus rhizosphere biodiversity and highlight the importance of using microbiome rhizosphere to promote plant growth and fight plant pathogens.

Origin of Fungal Isolates
The plant pathogenic fungi used in the present study were Neocosmospora (Fusarium) solani (MH999444), F. oxysporum (MH999445), F. equiseti (MH999443), and F. brachygibbosum (MH999442). These fungal species were previously isolated from symptomatic roots of citrus trees during the growing season 2017 in Morocco [1]. Each fungal species was subcultured on Potato Dextrose Agar (PDA) supplemented with antibiotics (Chloramphenicol and streptomycin sulfate both at 50 µg/mL) medium and incubated 7 days in darkness at 25 • C prior to experiments.

Isolation of Bacteria from Citrus Rhizosphere
The soil samples were randomly collected from citricultural areas from five regions of Morocco (Taroudant, Meknes, Taounate, Sidi Kacem, and Berkane) (Figure 7). Samples of rhizosphere soil samples with citrus roots were carefully collected. Thereafter, serial dilutions were prepared; 10 g of the soil from each sample is suspended in an Erlenmeyer flask containing 100 mL of sterile distilled water (SDW). After stirring for 30 min, resulting dilutions were spread on petri dishes containing PDA medium and then incubated at 28 • C for 72-96 h [89]. Colonies with different morphologies were selected and re-streaked on Luria Bertani medium (LB) until pure cultures are obtained. A total of 210 colonies were collected and used in subsequent screenings.
of rhizosphere soil samples with citrus roots were carefully collected. Thereafter, serial dilutions were prepared; 10 g of the soil from each sample is suspended in an Erlenmeyer flask containing 100 mL of sterile distilled water (SDW). After stirring for 30 min, resulting dilutions were spread on petri dishes containing PDA medium and then incubated at 28 °C for 72-96 h [89]. Colonies with different morphologies were selected and re-streaked on Luria Bertani medium (LB) until pure cultures are obtained. A total of 210 colonies were collected and used in subsequent screenings.

Screening of Antagonist Bacteria
The bacterial collection of citrus rhizosphere soil was tested in this study. The antagonism bioassay using dual culture technique by direct confrontation between bacterial isolates and fungal pathogens to screen the suitable antagonist with substantial capacity

Screening of Antagonist Bacteria
The bacterial collection of citrus rhizosphere soil was tested in this study. The antagonism bioassay using dual culture technique by direct confrontation between bacterial isolates and fungal pathogens to screen the suitable antagonist with substantial capacity of restricting the growth of N. solani on PDA medium. Each bacterial isolate was streaked at four equidistant streaks along the perimeter of the Petri dish. Then, a 7 mm diameter mycelial plug was taken from the edge of a 7 day-old-colony of N. solani and deposited onto the center of the agar plate containing PDA medium between different streaks of each bacterium. Plates containing only the mycelial plug served as controls. The plates were sealed with parafilm and incubated at 28 • C for 7 days. The presence or absence of the inhibition zone was subsequently noted by calculating the inhibition rate of mycelial growth (IR) after one week of incubation. IR was calculated using the following formula: IR (%) = (C − T)/C × 100 With: IR: inhibition rate; C: Diameter of the fungal colony in the control plates; T: Diameter of the fungal colony in the presence of the antagonist [30].

Effect of Bacterial Isolates on the Mycelial Structure of Neocosmospora Solani
Interactions between the pathogenic fungus and bacterial isolates were investigated. After one week of co-culture isolates, a part of the mycelium was taken from the zone of inhibition and observed under light microscope (Ceti Microscope) to reveal the existing hyphal damages or cytological changes caused by the antagonistic bacteria compared to the control.

Bacterial Identification
The twenty bacteria showed to be effective against N. solani in the in vitro bioassay were identified by molecular tools. The genomic DNA of the bacteria was extracted using the protocol described by Llop et al. [90]. The partial 16S rDNA genes of the genomic DNA of antagonist isolates were amplified using universal primers: FD1: 5 AGAGTTT-GATCCTGGCT CAG 3 and RP2: 5 GGTTACCTTGTTACGACTT 3' [91]. The PCR reaction mixture was performed in a total volume of 25 µL containing 5µL of PCR buffer (5x), 1 µL (10 µM) of each primer, and 0.2 µL (5 U/µL) of Bioline taq DNA polymerase (Bioline, London, UK) and 2,5 µL of DNA template, the rest of the volume was completed with SDW. The following cycling conditions were used: Initial denaturation at 96 • C for 4 min, followed by 35 cycles of denaturation at 96 • C for 10 s, then annealing at 52 • C for 40 s and 72 • C for 2 min, and finally an extension at 72 • C for 4 min in Thermal Cycler. PCR products were sequenced in both directions using sanger sequencing method. The obtained sequences were assembled using DNAMAN software (version 6.0, Lynnon Biosoft, Quebec Canada), and compared with other bacterial DNA sequences in the National Centre for Biotechnology Information's (National Center for Biotechnology Information (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Partial sequences of 16S rDNA were deposited in the Genbank under accession numbers listed in Table 1.

Evaluation of Biocontrol Activities and Plant Growth Stimulating Attributes
The mechanisms of biocontrol activity and plant growth promoting traits exercised by the selected twenty rhizobacterial isolates were investigated by various in vitro experiments including production of lytic enzymes, antibiotic metabolites, HCN siderophores, IAA, and phosphate solubilization.

Indirect Antagonist Activity of Selected Bacteria
Volatile antifungal compounds (VOCs): The ability of bacterial isolates to produce volatile VOCs against N. solani was investigated. Each bacterium was first inoculated in three streaks on LB medium for 24 h. Then, the lid of each plate was replaced with the bottom of another plate containing the PDA medium, inoculated with a 7 mm fresh mycelial plug of the pathogen. Subsequently, both the bottoms were sealed with transparent adhesive tape (Parafilm ® ) to prevent any loss of volatile substances [24,92].
The control was prepared in the same way except that the bottom contained no bacterium. Incubation takes place at 25 • C and the inhibition rates were noted after 7 days of incubation according to the following formula: IR= (C − T)/C × 100.
Antibiosis via bacterial supernatant: Antibiosis via bacterial supernatant was carried out by incorporation of metabolites produced by bacteria in order to assess the involvement of diffusible substances in antifungal activity [38]. An aliquote (100 µL) of each bacterial suspension (1 × 10 8 CFU/mL) was inoculated into nutrient broth medium (NB), then incubated in a rotatory shaker at 30 • C for 3 days (130 r/min). The mixture was centrifuged for 25 min (5500 rpm) and the supernatant obtained from each isolate was filtered through a 0.22 µm diameter Millipore filter. The bacterial cell-free filtrate was incorporated into a PDA agar culture medium (45-50 • C) at a concentration of 10% (v/v). A 7 mm mycelial plug of the pathogen obtained from an actively growing culture of a 7 day-old colony was placed at the centre of the plates and further incubated at 25 • C and observations noted after 7 days of incubation. Control plates were prepared by adding a 10% concentration of liquid NB medium to the PDA instead of the bacterial supernatant. The inhibition rate was calculated as described above. There were two independent trials with 3 replicates.

Microbial Traits
Proteolytic activity: The ability of antagonist bacteria to produce protease was determined using a solid medium based on skim milk. The medium was inoculated with a 5 µL (1×10 8 CFU/mL) of each bacterial suspension. Plates were incubated at 28 • C for 48 h. Protease activity was revealed by the development of a clear halo around the colonies [83]. The proteolytic index was then calculated as the diameter of halo (mm) + diameter of a colony (mm)/diameter of a colony (mm) as described by Syed-Ab-Rahman et al. [83].
Amylase activity: The ability of bacterial isolates to produce amylase was assessed using a solid medium supplemented with soluble starch [93]. 5 µL of each bacterial culture (1×10 8 CFU/mL) was spotted in the Petri dish and incubated at 28 • C for 72 h. In order to reveal the hydrolysis of the starch, the agar surface was covered with 3 mL of the iodine solution. After 3 min, the appearance of a clear zone around the colony indicated a positive amylase activity. In the absence of amylase activity, starch turns a blackish blue color. Thus, the amylolytic index is calculated as previously described [83].
Cellulase degradation: The bacterial antagonists were tested for their ability to produce cellulase using a solid medium supplemented with carboxymethyl cellulose (CMC). The medium was inoculated with a 5 µL (1 × 10 8 CFU/mL) of each bacterial suspension. Plates were incubated at 28 • C for 72 h. Cellulase production was revealed by pouring Red Congo solution at a concentration of 0.1% on the surface of the Petri dishes for15 min and destaining with a solution of NaCl (1M) (by rinsing 3 times) [83]. Development of a clear halo around the colonies confirm the presence of cellulase activity, while the color remains red in the absence of cellulase activity. The cellulose index was calculated as previously described [83].
Pectinase activity: The ability of bacterial isolates to produce pectinases was determined using a solid medium containing pectin as previously described by Etesami et al. [94]. The medium was inoculated with a 5 µL of each bacterial suspension (1 × 10 8 CFU/mL). The Petri dishes were incubated at 28 • C for 72 h. The pectinolytic activity was revealed by the addition of Cetyl Trimethylammonium Bromide (CTAB) at 1%. The formation of a clear halo around the colonies indicated positive pectinolytic activity. Thus, the pectinolytic index was calculated as described above.
Phosphate Solubilization: Pikovskaya (PVK) medium amended with 5 g/l of tricalcium phosphate (Ca3 (PO4) 2), as the sole source of phosphate, was used to test the ability of bacterial isolates to solubilize inorganic phosphate, as previously described [22,83]. The medium was inoculated with a 5 µL (1 × 10 8 CFU/mL) of each bacterial suspension. Inoculated plates were incubated at 28 • C for 4 days. Bacterium capable of solubilizing phosphate will be surrounded by a clear halo, thus the phosphate solubilization index was calculated as described above.
Siderophores Production: Siderophores production was determined using chrome azurol S (CAS-shuttle) assay [40,97] with slight modifications. Each bacterial suspension (1×10 8 CFU/mL) was placed into 15 mL falcon tubes containing 10 mL of the liquid AT minimal medium devoid of iron [98]. After 96 h of incubation at 28 • C with stirring at 150 rpm, 0.5 mL of the supernatant was mixed with 0.5 mL of the CAS-HDTMA solution. The solution of CAS-HDTMA was prepared as follows: The ferric ion solution was freshly prepared by mixing 1.5 mL FeCl3. 6H2O (1 mM) with HCl (10 mM) and 7.5 mL of Chrome Azurol S (2 mM), stirred slowly with the addition of 6 mL hexadecyl trimethyl ammonium bromide (HDTMA), (10 mM) and shaken, the piperazine solution (4.307g of piperazine dissolved in 30 mL of water and 6.25 of HCl (12 M); pH = 5.6) was slowly added and made up to 100 mL. A tube containing 0.5 mL of the uninoculated minimal medium broth mixed with 0.5 mL of CAS-HDTMA served as a control. After 2 h of contact between the supernatant and the CAS-HDTMA solution, the blue coloring turns orange or purple when siderophores are produced.
IAA Production: The production of indole-3-acetic acid (IAA) was determined by the colorimetric method as described by Yuttavanichakul et al. [60]. An aliquote (100 µL) of each bacterial suspension were cultured in a liquid LB medium supplemented with L-tryptophan (1 g/L) and incubated at 28 • C on a rotatory shaker at 150 rpm for 4 days. Cultures were centrifuged at 5000 rpm for 20 min. Subsequently, 1 mL of the cell-free culture supernatant was mixed with 2 mL of Salkowski's reagent (12 g of FeCl3 per litre of 7.9 M H2SO4) and development of color was observed. The appearance of a pinkishred coloration indicates IAA production, while yellow coloration indicates a negative result [40,88].
Production of Hydrogen Cyanide (HCN): The ability of bacterial isolates to produce hydrocyanic acid (HCN) was examined following the protocol described by Lahlali et al. [30]. 100 µL of each bacterial suspension (1 × 10 8 CFU/mL) were spread on LPGA medium supplemented with 4.4 g glycine per litre (4.4 g/L). Subsequently, sterile Whatman filter paper (no.1) discs were saturated with a picrate solution (2.5% picric acid in 12.5% anhydrous sodium carbonate (Na 2 CO 3 ) solution) were placed on the lid of the Petri dish. The control plates were inoculated with SDW. The plates were sealed with parafilm and incubated at 28 • C for 4 days. The change of color from yellow to orange, red, or reddish brown indicated volatile HCN production [30].
Detection of Lipopeptides by the PCR Method: The total genomic DNA extracted from the twenty selected rhizobacteria were used for the detection of the biosynthesis lipopeptides (bacillomycin, fengycin, iturin, and surfactin) and bacteriocin (subtilosin A) genes. Each PCR amplifications were performed in a total volume of 25 µL of PCR mixture containing 5µL of PCR buffer (5×), 1 µL of each primer (10 µM), 0.25 µL of Taq DNA polymerase (5 U/µL) (Bioline, London, UK), 2.5 µL of genomic DNA, the rest of the volume was completed with SDW. Specific primers used for amplification of these gene was used (Table S1). PCR reactions were performed Thermal Cycler. PCR products were then visualized on a 1.5% agarose gel colored using cyber safe (Invitrogen, CA, USA) by electrophoresis and visualized with an ultraviolet illuminator and digitally recorded.

Hypersensitivity Test
The phytopathogenicity of the twenty antagonist bacteria was examined using the hypersensitivity test on tobacco (Nicotiana tabacum). An aliquote (100 µL) of each bacterial suspension (1 × 10 8 CFU/mL) was injected at the level of the midrib of the lower part of the leaf. For the positive and negative controls, a bacterial suspension of a plant pathogen (Agrobacterium tumefaciens) and SDW were used, respectively. The treated plants were kept in a growth chamber at a temperature of 27 • C and observed after 24 to 48 h. A phytopathogenic bacterium elicits a positive hypersensitivity response that results in leaf dryness and brown necrosis.

PGPR Test on Brassica Napus
The twenty effective bacterial isolates (IR > 60%) were tested for their PGPR capabilities in planta. Bacterial isolates were cultured in yeast extract peptone (YEP) medium at 28 • C with shaking. After 24 h of incubation, the cell cultures were centrifuged at 5000 rpm for 20 min, the supernatant was removed, and the pellet was suspended in phosphate buffered saline (PBS) at a concentration of 1 × 10 8 CFU/mL (OD = 0.1 to λ = 600 nm). Brassica napus was used due to its fast growth, large biomass production, and high germination rate [99,100]. Canola seeds, with uniform shape and size, used in this trial were surface sterilized with 95% ethanol for 30 s, washed with SDW 5 times, and air dried at room temperature. Subsequently, they were soaked in the bacterial culture diluted in PBS for 1 h at room temperature. Seeds soaked in the PBS solution were served as a control.
The soil was sterilized twice by autocalving at 121 • C for 90 min to destroy the microorganisms it contains. Subsequently, 140 g of soil was placed in pots (20 cm × 7 cm) previously disinfected with sodium hypochlorite and dried at ambient temperature. One seed was planted per pot and pots were arranged in complete randomized design with 5 replications. This trial was carried out in a controlled greenhouse at a temperature of 25 • C. After 30 days, three pots from each treatment were randomly selected and the plants/roots length and the plants/roots fresh/dry weight were measured and recorded [69,83].

Statistical Analysis
All in vitro and in planta experiments were carried out twice over time. All tests were carried out using a completely randomized design. All datasets were summarized as mean ± SD (standard deviation). The Arcsin transformation was applied to assess the disease severity. All statistical analysis were performed using SPSS statistical software (version 20, IBM SPSS Statistics 20, New York, NY, USA) and when the effect was revealed to be significant, Duncan test was performed for means separation at a significance level of p ≤ 0.05.

Conclusions
In the present study, research of effective biocontrol agents for the control of dry root rot of citrus, was started by in vitro screening of rhizobacteria. In light of the results obtained, twenty bacteria were selected to have an important antifungal activity against N. solani and tested also against three other Fusarium species associated with this disease. The selected bacteria were studied for their biochemical characteristic and evaluated for their abilities to enhance the growth of B. napus under greenhouse conditions. Ten isolates were also tested for their capacity to control the DRR disease under greenhouse conditions. Two bacteria appear promising since they were found to completely suppress the disease. These antagonist bacteria, isolated from the soil of healthy citrus, fulfil the hypothesis according to which microorganisms chosen as BCAs should be screened from their local niches, and used in the same environment in order to obtain the desired benefits. In addition, their plant growth traits were sought as possible additional mechanisms. Undoubtedly, the encouraging results of this study are important and open new alternatives toward the design of biocontrol strategies for managing DRR and limit losses in citrus crops. Therefore, two antagonists bacteria B. subtilis K4-4 and GH3-8 were proposed to control and prevent damages of DRR disease. However, to confirm their biocontrol potential at a large-scale, further experiments under natural conditions and during environmental stress would also be necessary through inoculations in citrus orchards and especially during stress conditions. Furthermore, the implementation of BCAs at a large-scale faces several challenges, which depends on the advances of interdisciplinary research, in particular mass production and formulation methods, which are crucial factors to have effective bacterial inoculums with higher reliability and competitiveness on the market.