Endophytic Bacterial Isolates From Halophytes Demonstrate Phytopathogen Biocontrol and Plant Growth Promotion Under High Salinity

Halophytic endophytes potentially contribute to the host’s adaptation to adverse environments, improving its tolerance against various biotic and abiotic stresses. Here, we identified the culturable endophytic bacteria of three crop wild relative (CWR) halophytes: Cakile maritima, Matthiola tricuspidata, and Crithmum maritimum. In the present study, the potential of these isolates to improve crop adaptations to various stresses was investigated, using both in vitro and in-planta approaches. Endophytic isolates were identified by their 16S rRNA gene sequence and evaluated for their ability to: grow in vitro in high levels of NaCl; inhibit the growth of the economically important phytopathogens Verticillium dahliae, Ralstonia solanacearum, and Clavibacter michiganensis and the human pathogen Aspergillus fumigatus; provide salt tolerance in-planta; and provide growth promoting effect in-planta. Genomes of selected isolates were sequenced. In total, 115 endophytic isolates were identified. At least 16 isolates demonstrated growth under increased salinity, plant growth promotion and phytopathogen antagonistic activity. Three showed in-planta suppression of Verticillium growth. Furthermore, representatives of three novel species were identified: two Pseudomonas species and one Arthrobacter. This study provides proof-of-concept that the endophytes from CWR halophytes can be used as “bio-inoculants,” for the enhancement of growth and stress tolerance in crops, including the high-salinity stress.

Introduction

Bacterial endophytes are widespread among plants and colonize intercellular and intracellular spaces of all host compartments. Each individual plant is a host to bacterial and fungal endophytes that colonize its tissues for all or part of its life cycle without causing any apparent pathogenesis (Ryan et al., 2008). Various studies have shown how microbial communities contribute to plant defense and the substantial beneficial effects they have on host plants, including improved nutrient acquisition, accelerated growth, resilience against pathogens and improved resistance against abiotic stress such as heat, drought, and salinity (Rodriguez et al., 2019).

The diversity and structure of endophytic microbiomes are dynamic and directly affected by ecological characteristics of the host plant and soil such as geographic location, environmental factors and interactions within the host plant (Edwards et al., 2015). Most characterized members of bacterial endophytic communities belong to the Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria (Bulgarelli et al., 2013, 2015; Edwards et al., 2015). However, endophytic microbiome structure can be affected by the host-plant species’ genotype, plant organ or tissue type, developmental stage, growing season, geographic, and field conditions, soil type, nutrient status of the host species and cultivation practices (Rodríguez-Blanco et al., 2015; Liu et al., 2017; Rodriguez et al., 2019).

Endophytic microbes hold enormous potential to increase plant health. Interestingly, endophytic bacteria can be used to overcome the effect of salinity stress, promote plant growth and nutrient uptake; these approaches can provide beneficial and environmentally friendly solutions for a sustainable global food security (Glick, 2014; Tkacz and Poole, 2015; Vaishnav et al., 2019). For successful exploitation of endophytes, we need a deeper understanding of endophytic community composition and the mechanisms that underlie their plant growth promotion, in order to successfully select the most efficient bacterial isolates.

Members of endophytic bacterial communities influence each other with antagonistic, competitive, and mutualistic interactions (Toju et al., 2018). This results from nutritional competition, exchange and even metabolic interdependence. This, in turn can influence microbiome composition and its effect on the host-plant (Rodriguez et al., 2019). Host-plant genotype can also have a dramatic impact on microbial members; individual cultivars can influence the microbial community structure and the beneficial effects of endophytic bacteria (Haney et al., 2015; Marques et al., 2015; Pérez-Jaramillo et al., 2016; Rodriguez et al., 2019). Thus, for the utilization of endophytic bacterial isolates, an optimum approach is to isolate key bacterial strains from crop wild relatives (CWRs) (Mendes et al., 2013).

Halophytes could be valuable sources of novel endophytic isolates that can be used to overcome various biotic and abiotic stresses (Ruppel et al., 2013; Shabala, 2013; Yuan et al., 2016; Etesami and Beattie, 2018). High salinity in plants results in ionic and osmotic stress due increased extracellular hypertonic conditions and accumulation of Na⁺ and Cl^– intracellularly (Vaishnav et al., 2019). The resulting stress affects intracellular water balance, rate of cell division, hormonal imbalance, changes in photosynthesis, nutrient translocation, processes that decrease plant growth (Munns, 2002). Plant-associated microorganism can contribute to plant health impeded by salinity stress, by influencing phytohormonal levels and signaling, contributing to homeostasis maintenance of toxic ions under salinity stress, enhancing photosynthesis, and contributing to biomass production and allocation (Dodd and Pérez-Alfocea, 2012). Since soil salinity disrupts the physiological and morphological plant processes and increases pathogen susceptibility (Etesami and Beattie, 2018), the use of plant growth promoting endophytes in crops can be a more eco-friendly approach than agricultural chemicals.

Here, we tested the hypothesis that cultivated endophytic bacteria isolated from three CWR halophytic plant species have properties of salinity stress tolerance, plant growth promotion and phytopathogen growth inhibition. These species included two members of the Brassicaceae family (Matthiola tricuspidata and Cakile maritima), and one of the Apiaceae family (Crithmum maritimum). To test this hypothesis, we cultured and identified 115 different bacterial isolates and functionally characterized them in in vitro and in-planta assays. The bacterial isolates were tested in vitro for their ability to grow on salinity levels up to 17.5%, their biocontrol of the economically important plant phytopathogens Verticillium dahliae, Ralstonia solanacearum, and Clavibacter michiganensis ssp. michiganensis. Subsequently, isolates with demonstrated in vitro salt tolerance, were tested in-planta to demonstrate whether they promoted plant growth under no stress conditions and under high salinity. Furthermore, bacterial isolates were tested in-planta to check their biocontrol properties against Verticillium dahliae. This is the first study of bacterial endophytes obtained from M. tricuspidata, Cr. maritimum, and Ca. maritima, and identifies their potential as bacterial bio-inoculants in commercial crops to overcome salinity stress and plant diseases caused by the economically important pathogens.

Materials and Methods

Site Description and Plant Sample Collection

Samples were collected during summer 2018 in three distinct sites in Crete, Greece: site 1 (S1: N35°25′. E24°41′), site 2 (S2: N35°06, E25°48), and site 3 (S3: N35°00, E25°44) (Supplementary Figure 1). At S1, a natural beach area favoring salt-marsh vegetation, three Matthiola tricuspidata individuals were collected. At S2, a beach area close to Pachia Ammos village, three Crithmum maritimum individuals were collected. At S3, a popular beach area located in the town of Ierapetra, three individuals of Cakile maritima were collected. Each sample was collected with sterile gloves, forceps and gloves, placed in separate plastic bags to avoid cross contamination and immediately transported to the laboratory for processing.

Plant Surface Sterilization and Endophytic Cell Isolation

Leaf and root materials from each species were cut and processed individually. Plant material was gently washed with sterile distilled water repeatedly to remove soil and dust particles. For surface sterilization, plant roots and leaves were placed into sterile Erlenmeyer flasks containing ethanol 75% v/v for 60 s with shaking and then in sterile Erlenmeyer flasks containing sodium hypochlorite solution 3% w/v^–1 (NaClO) for 10 min. The plant materials were then placed again in ethanol 75% v/v for 60 s. To remove any remaining NaClO, plant materials were rinsed 10 times with sterile distilled water (dH₂O). The sterilization and transfer procedures were carried out in a type II laminar flow hood. About 100 μL of the last rinse (for each analyzed sample) was plated on Nutrient Agar (NA) medium and monitored for microbial growth to evaluated surface sterilization efficiency. Only successfully sterilized root material was used further. Approximately 500 mg of leaves and roots per each species were weighed and slashed to small parts for further processing using a sterile scalpel and further grounded into a slurry with an autoclaved pestle and mortar. The slurry was transferred into sterile petri dishes and 30 mL of autoclaved dH₂O was added. The petri dishes were sealed and placed onto a rotary shaker (150 rpm) at 25°C for 2 h. After shaking, 100 μl of the material in triplicate were inoculated on NA plates and incubated at 28°C. Colony forming units (cfus) were chosen from each plate based on their color, texture and morphology. Pure bacterial colonies were grown in Nutrient Broth (NB) and cells stocks were stored in 50% v/v glycerol at −80°C. A total of 115 isolates were identified.

Bacterial Isolation and Identification of Isolates

To identify the 115 bacterial isolates, 16S rRNA gene Sanger sequencing method was employed. To extract crude genomic DNA, 1 mL of liquid bacterial culture in NB was placed in liquid nitrogen for 15 s. After room temperature incubation, the lysate was centrifuged at 10,000 × g for 1 min. Two microliter of the lysate were used to amplify the 16S rRNA gene using primers 27F: 5′-AGAGTTTGATCCTGGCTCAG-3′ (White et al., 1990) and 1492R: 5′-GTTTACCTTGTTACGACTT-3′ (Lane, 1991). PCR reactions of 20 μL were amplified in a BioRad T-100 Thermocycler with initial denaturation at 94°C for 2 min, followed by 35 cycles of 5 s at 94°C, 30 s annealing at 47°C, 2 min primer extension at 72°C, and a final extension at 72°C for 5 min. Apart from the lysate, each tube contained, Bac-Free PCR Buffer, 250 nM of each primer, 0.2 mM of each deoxy- ribonucleotide triphosphate and 0.1 U BAC-Free HotStart Taq polymerase (Nippon Genetics, Europe). PCR products were purified using Nucleo Spin Gel and PCR Clean up (Macherey-Nagel, Germany). Cleaned-up PCR products were sent to Macrogen (Europe) for sequencing with primer 27F.

The resulting chromatograms were quality inspected using MEGA 5 (Kumar et al., 2018) and the start/end regions of low quality were manually trimmed off. Cleaned-up fasta files were aligned in SILVA (Quast et al., 2013). The resulting sequences of the 16S rRNA gene were queried against ezBioCloud (Yoon et al., 2017) reference database for identification and documentation of the described bacterial isolate with the closest sequence similarity.

Bacterial Salt Tolerance Assay

The salt tolerance of all bacterial isolates was estimated on the basis of the population density of these isolates at different concentrations of NaCl (ranging from 0.5, 5, 10, 15, and 17.5% (w/v) in NA. Ten microliter drops of freshly prepared NB cultures of each isolate were inoculated on sterilized petri plates, containing 25 mL NA with increasing NaCl concentrations and incubated at 28°C. For each NaCl concentration, an Escherichia coli laboratory isolate was inoculated as a negative control. After 24 h of incubation, the growth of each isolate was estimated compared to E. coli growth.

In vitro Growth Inhibition of Phytopathogens

Antibacterial activity of the bacterial isolates against the phytopathogenic bacteria Ralstonia solanacearum, and Clavibacter michiganensis was evaluated by co-culturing each of the bacterial isolates on NA lawn covered by R. solanacearum or C. michiganensis. The inhibition zone indicating inhibition by bacterial growth was recorded as the antibacterial effect. Antifungal activity of the isolates against Verticillium dahliae was investigated. Potato dextrose agar (PDA) was inoculated with each bacterial isolate for 24 h at 28°C and then V. dahliae was inoculated at room temperature for 3–4 weeks. Fungal growth inhibition was determined by measuring the inhibition zone of V. dahliae hyphae on the media.

In vitro Hemolysis Screening Assay

In order to assay the bacterial isolates for hemolytic activity, each isolate was grown on blood agar plates. The bacterial isolates were inoculated with the spot test method and were incubated at room temperature for 48 h. The known non-mammalian-pathogenic species Ensifer meliloti was employed as a negative control (Supplementary Figure 2B).

In vitro Growth Inhibition of Fungal Human Pathogen

Antifungal activity of specific isolates against anthropopathogenic fungus A spergillus fumigatus was evaluated by co-culturing 11 bacterial isolates on NA plate lawn covered by A. fumigatus for 72 h at room temperature under absence of light. The following isolates were tested: CML04, CMR11, CMR22, CMR25, CrR12, CrR25, MTR12, MTR17a, MTR17b, MTR17c, MTR17d. Fungal growth inhibition was determined by the growth inhibition zone of the A. fumigatus hyphae on the media.

In-planta Salt Tolerance Assays

Twelve of the bacterial isolates were selected, according to their ability to grow in high salinity conditions (up to 17.5% w/v NaCl), in order to test their plant growth promotion capacity of the model plant Arabidopsis thaliana. Firstly, the experiment was performed with no abiotic stress conditions. Bacterial isolates were cultured in NB for 46 h at 25°C with stirring. NB cultures were centrifuged at 224 × g for 15 min, the supernatant was discarded and the cells were resuspended in 50 mL sterilized dH₂O. Seeds of A. thaliana ecotype Columbia (Col-0) were grown in plastic pots (6 × 6 × 7 cm) filled with vermiculite: soil (1:1), at 25°C (16 h light/8 h dark). For each isolate and the corresponding control, 5 individual plants were grown in each pot. A. thaliana plants were watered with dH₂O for 10 days. Then, plants were watered with 10 mL suspensions of the 12 bacterial isolate liquid cultures, and were left for 7 days to let the bacterial isolates adapt. Subsequently, for a 30 day span, plants were watered every 2–3 days with dH₂O. At the end of the experiment, the fresh weight of the leaves from each plant was measured. The leaves were then dried at 65–70°C for 2 days and their dry weight was measured.

The same experiment was performed under salt treatment. Specifically, after the 7 day period of bacterial isolate inoculation, instead of dH₂O, the plants were watered with 10 mL of 250 mM NaCl. Fresh and dry weight of the leaves was measured.

For both experiments, mock samples were employed where no bacterial isolates were inoculated and control plants were inoculated with the isolate Escherichia coli (Control-E. coli), to check that the plants would not use the bacteria as a fertilizer.

Confrontation and Volatile Tests of Selected Bacterial Isolates Against Verticillium dahliae

A total of 16 isolates were selected for direct in vitro antagonism of V. dahliae. Fourteen of these isolates were selected due to their strong inhibition of V. dahliae in initial tests (Supplementary Table 1; CrR14, CrR18, MTR18, CMR01, CMR03, CML04, CMR25, MTR17a, MTR17d, MTR17f, MTR17g, MTR17h, and MTR17b, MTR17c). Two additional isolates were selected (CrR04 and MTR12) with medium inhibition in initial tests (Supplementary Table 1) for comparison. Direct in vitro antagonism of V. dahliae was evaluated by dual-culture assays (confrontation test) on PDA (Lahlali et al., 2007). In particular, a 6 mm diameter mycelial disc taken from the periphery of a 2 week-old PDA fungal culture was placed on a new PDA plate (90 mm in diameter) at approximately 25 mm-distance from the center of the plate. Then, a 30 mm-long line from each bacterial isolate (taken from a 48 h-old tripticasein soy broth (TSB) liquid culture with an inoculation loop) was streaked on the opposite site of the plate at equal distance from the center (one isolate per plate). Moreover, Trichoderma harzianum strain T22 was isolated from the commercial biofungicide TRIANUM-P (Koppert B.V. Hellas) and included in in vitro bioassays for comparison. Plates inoculated only with V. dahliae agar discs were served as controls. Plates (three per bacterial isolate plus controls) were incubated at 24°C in the dark. The radius of fungal colonies toward the direction of the test isolate and that of controls was measured 5, 7, 9, and 12 days post inoculation (d.p.i.) and radial growth rates were expressed in mm/day. At the end of the bioassays (12 d.p.i.) the underside of the plates was scanned using a Samsung Xpress SL-M2875ND Laser Multifunction Printer at 1200 dpi and microsclerotial (black) area on each plate image was determined manually using the image processing software ImageJ 1.46r (National Institutes of Health, United States). Then, the number of spores was estimated by transferring a 6 mm-diameter disc taken from the periphery of each culture into a 1.5 mL Eppendorf tube with 1 mL of water, and vortexed for 30 s. The number of spores was measured using a haematocytometer under a light microscope. Moreover, actively growing mycelia from cultures’ periphery (located closer to test isolate) were prepared and microscopic observations (30 readings per culture) were carried out to estimate hyphae width.

To evaluate the capacity of bacterial isolates to affect V. dahliae growth via the production of volatile compounds, dual-plate assays (Chaurasia et al., 2005) were conducted (volatile test). In brief, one 6 mm-diameter agar disc of actively growing mycelium of the fungus was placed in the center of a new PDA plate (90 mm in diameter), whilst each bacterial isolate (taken from a 48 h-old TSB liquid culture) was streaked on another PDA plate. The covers of the two plates were removed and resultant plates were adjusted together (bacterial culture was upturned) and sealed with cellophane membrane so the two microbes would share the same headspace without coming in contact with each other. Dual plates (upright and upturned) inoculated only with V. dahliae served as controls. Similarly, in dual-culture assays, dual-plates (three per bacterial isolate) were incubated at 24°C in the dark and the growth, microclerotial area, sporulation and hyphae width of fungal colonies were measured as described above.

Radial growth inhibition (RGI), microsclerotia formation inhibition (MFI), sporulation inhibition (SI) and hyphae thinning (HT) were calculated according to the formula: [(Vc-Vt)/Vc] × 100 where Vc = the microscopic value of V. dahliae in control plates and Vt = the respective value of V. dahliae against the antagonistic isolate in dual-culture or dual-plate assays.

Bacterial Isolates and Fungal Inoculum Preparation for in-planta Bioassays

The 16 selected (see above) bacterial isolates (CrR14, CrR18, CrR04, MTR12, MTR18, CMR01, CMR03, CML04, CMR25, MTR17a, MTR17d, MTR17f, MTR17g, MTR17h, and MTR17b, MTR17c) were used in in-planta bioassays. The isolates were grown in Erlenmeyer flasks with 200 mL liquid TSB, in an orbital incubator at 180 rpm and 28°C for 48 h in the dark. Bacterial suspensions were centrifuged at 3,000 × g for 10 min and cells were re-suspended in water reaching a final concentration of 10⁸ cfu mL^–1 (measured by dilution plating).

The highly virulent V. dahliae isolate 999-1 (Markakis et al., 2016), which originated from symptomatic eggplants (Solanum melongena L.), was used. V. dahliae conidial suspension for eggplant—V. dahliae bioassays was prepared as previously described (Markakis et al., 2016). In brief, conidia were produced by growing each V. dahliae strain in potato dextrose broth (PDB) at 160 rpm and 25°C in the dark for 5 days. Then, conidia were harvested by filtrating through three layers of cheesecloth and the suspensions centrifuged at 3,000 × g for 10 min. Spores were re-suspended in sterilized dH₂O and their concentration was adjusted to 5 × 10⁶ conidia mL^–1.

In-planta Verticillium Wilt Suppression Bioassays

Eggplant seedlings (cv. Black Beauty) were used in the in-planta bioassays. Plants at the one-true-leaf stage, grown in 100 mL-capacity pots containing soil substrate (HuminSubstrat, Klasmann-3 Deilmann GmbH, Germany) were root-drenched with bacterial suspension (20 mL of 10⁸ cfu mL^–1 of each isolate per plant), whereas plants that served as controls (negative = no bacterium/no fungus assigned as “C−” and positive = no bacterium/plus pathogen assigned as “V.D.”) were treated with 20 mL of water. One week later, eggplants (at the second-true-leaf stage) were inoculated with V. dahliae by drenching the soil substrate in each pot with conidial suspension (20 mL of 5 × 10⁶ conidia mL^–1 per pot). Negative control plants (C−) were treated with 20 mL of water. Eggplants were maintained under controlled conditions at 23 ± 2°C with a 12 h light and dark cycle.

Two independent experiments (experiments I and II ) were conducted to evaluate the suppressive effect of the aforementioned bacterial isolates against V. dahliae. In experiment I, 11 treatments (C−, V.d., V.d. + CrR14, V.d. + CrR18, V.d. + CrR04, V.d. + MTR12, V.d. + MTR18, V.d. + CMR01, V.d. + CMR03, V.d. + CML04 and V.d. + CMR25) were applied; whereas in experiment II, 10 treatments were conducted (C−, V.d., V.d. + MTR17a, V.d. + MTR17d, V.d. + BMTR17f, V.d. + MTR17g, V.d. + MTR17h, V.d. + MTR17b, V.d. + MTR17c, and V.d. + TRIANUM-P). The commercial biofungicide TRIANUM-P was included in experiment II (assigned as V.d. + TRIANUM-P) and applied according to manufacturer’s instruction (20 mL of 3 × 10⁷ cfu mL^–1 per plant). TRIANUM-P was served as a V. dahliae-suppressive reference treatment. Within each experiment, each treatment consisted of seven plants and experiments were replicated three times.

Disease Assessment

Verticillium wilt symptoms on eggplant were recorded at 2-, 3-, and 4- day intervals from 12 to 30 d.p.i with V. dahliae. Bioassays were evaluated by estimating disease severity, disease incidence, mortality and relative area under disease progress curve (RAUDPC). Disease parameters were recorded as previously described (Markakis et al., 2016). Briefly, disease severity at each observation was calculated from the number of wilting leaves, as a percentage of total number of leaves per each plant. Disease ratings were plotted over time to generate disease progress curves. Subsequently the area under disease progress curve (AUDPC) was calculated by the trapezoidal integration method (Campbell and Madden, 1990). Disease was expressed as a percentage of the maximum possible area with reference to the maximum value potential reached over the whole period of each experiment and is referred to as RAUDPC. Disease incidence was estimated as the percentage of infected plants. Only plants with a final disease severity of ≥20% were considered infected, to discriminate between V. dahliae-associated disease symptoms and other weak symptoms occasionally observed (Supplementary Table 2). Mortality was estimated as the percentage of dead plants.

Plant Growth

Growth parameters were evaluated at the end of bioassays (at 24 and 30 d.p.i. for experiments I and II, respectively). To estimate the effect of the aforementioned treatments on plant growth, all plants were clipped off at the soil surface level and their height, fresh weight and leaf number were measured.

Fungal Pathogen Re-isolation

To verify the presence of the applied V. dahliae strain in plant tissues, five plants per treatment in each experiment were randomly selected. Eggplant leaves which had been cut above soil level previously were removed and their stems were surface-disinfected by spraying with 95% ethyl alcohol and by quickly passing them through flame three times. For each plant, 3 xylem chips taken from different sites along the stem (base, middle and upper part of the stem) and aseptically placed onto acidified PDA after removing the phloem. Plates were then incubated at 24°C in the dark for 14 days. The emerging fungi that grew out of tissue excisions were examined visually and under a light microscope and identified according to their morphological characteristics (Pegg and Brady, 2002). Pathogen isolation ratio was expressed as the frequency of positive V. dahliae isolations of each plant.

Statistics

Analysis of variance (ANOVA) was employed to determine the effects of replication (1, 2, or 3), treatment (C−, V.d., V.d. + CrR14, V.d. + CrR18, V.d. + CrR04, V.d. + MTR12, V.d. + MTR18, V.d. + CMR01, V.d. + CMR03, V.d. + CML04, V.d. + CMR25 in Experiment I and C−, V.d., V.d. + MTR17a, V.d. + MTR17d, V.d. + BMTR17f, V.d. + MTR17g, V.d. + MTR17h, V.d. + MTR17b, V.d. + MTR17c, V.d. + TRIANUM-P in experiment II ) and their interaction on disease incidence (DI), final disease severity (FDS), mortality (M), RAUDPC and isolation ratio (IR), and on plant height, fresh weight and total number of leaves (Supplementary Tables 2, 3). Prior to ANOVA, normality of data and homogeneity of variance across treatments was evaluated and an arcsine transformation was applied to normalize variance. When a significant F test was obtained for treatments (P ≤ 0.05), the data were subjected to means separation by Tukey’s honestly significant difference test. Morphological and physiological characteristics of V. dahliae in dual-culture and dual-plate assays were also analyzed by Tukey’s test (P ≤ 0.05). Moreover, standard errors of means were calculated.

Bacterial Genome Sequencing and Annotation

Twelve bacterial isolates were selected for whole genome sequencing (CrR16, CMR16, CrR07, CMR13, CrR06, CrR18, CrR14, CMR27, CMR25, CML04, CMR29, CrR25). The isolates were selected when they met more than two of the following criteria: (a) the 16S rRNA gene sequence of the isolates being 99.6% similar to their closest relative or lower, (b) exhibiting salt tolerance higher than the 5% threshold, (c) exhibiting medium or strong inhibition against the growth of at least 2 of the 3 tested phytopathogens Verticillium dahliae, Ralstonia solanacearum, and Clavibacter michiganensis ssp. michiganensis (Supplementary Table 1). For each isolate a 250 bp paired-end library was produced for use with the Illumina MiSeq sequencing system (University of Exeter Sequencing Service, Exeter, United Kingdom). Reads were assembled using SPAdes 3.12.0 (Bankevich et al., 2012) and the assembled sequence was annotated using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP). Raw sequence reads and assembled genomes were uploaded to the Sequence Read Archive (Leinonen et al., 2011) and GenBank (Dennis Benson et al., 2017) and are available under BioProject accession number PRJNA634334. RAST (Rapid Annotation using Subsystem Technology) (Overbeek et al., 2014) was employed for genome analysis and annotation.

Sequence Alignment and Phylogenetic Tree Construction

Selected gene sequences were aligned with ClustalX v2.0 (Larkin et al., 2007) and subsequently manually corrected. Sequence relationships were inferred using the maximum-likelihood (ML) method. ML phylogenies were constructed using MEGA 5.2 (Tamura et al., 2011). Phylogenetic trees were constructed using the concatenated recA and gyrB genes and assuming the bootstrap value derived from 1,500 replicates to represent the evolutionary history of the included taxa.

The evolutionary history of Arthrobacter recA-gyrB genes were inferred by using the Maximum Likelihood method based on the Tamura-Nei model (Tamura and Nei, 1993). The bootstrap consensus tree inferred from 500 replicates (Felsenstein, 1985) was taken to represent the evolutionary history of the taxa analyzed (Felsenstein, 1985). Branches corresponding to partitions reproduced in less than 50% bootstrap replicates are collapsed. Initial tree(s) for the heuristic search were obtained automatically as follows. When the number of common sites was <100 or less than one fourth of the total number of sites, the maximum parsimony method was used; otherwise BIONJ method with MCL distance matrix was used. Trees were drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 25 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd + Non-coding. All positions with less than 95% site coverage were eliminated. That is, fewer than 5% alignment gaps, missing data, and ambiguous bases were allowed at any position. Evolutionary analyses were conducted in MEGA5 (Tamura et al., 2011).

The evolutionary history Arthrobacter recA-gyrB genes was inferred with the Maximum Likelihood method as above. The tree with the highest log likelihood (−7831.6808) was selected. Initial tree(s) for the heuristic search were obtained automatically as follows. When the number of common sites was <100 or less than one fourth of the total number of sites, the maximum parsimony method was used; otherwise BIONJ method with MCL distance matrix was used. The bootstrap consensus tree inferred from 500 replicates (Felsenstein, 1985) is taken to represent the evolutionary history of the taxa analyzed (Felsenstein, 1985). The tree was drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 22 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd + Non-coding. All positions with less than 95% site coverage were eliminated as described above. Evolutionary analyses conducted in MEGA5 as described above.

Results

Identification and Abundance of Culturable Endophytic Bacteria

Endophytic bacteria were cultivated from different surface-sterilized tissue samples from all three halophytes. A total of 115 pure bacterial cultures showing different colony morphology (from root or leaf) were obtained; 91 were retrieved from roots and 24 from leaves. In detail, 45, 31, and 39 isolates were obtained from M. tricuspidate, Ca. maritima, and Cr. Maritimum, respectively (Supplementary Table 1).

For all 115 isolates, total 16S rRNA gene sequencing allowed for taxonomic analysis (Figure 1 and Supplementary Table 1). Bacterial isolates were assigned to 5 different classes (Supplementary Table 1) and 24 genera (Figure 1 and Supplementary Table 1). The most prevalent genus was Bacillus, accounting for 24% of the isolates, followed by Enterobacter (19%) and Pseudomonas (12%). The highest number of bacteria were isolated from roots of M. tricuspidata (38), followed by Ca. maritima roots (28) and Cr. maritimum roots (25), in contrast to the number of bacteria isolated from leaf samples (M. tricuspidata: 7, Ca. maritima: 4 and Cr. maritimum: 13). Isolates of genus Bacillus were isolated from all plants and tissues except the roots of Cr. maritimum, while Pseudomonas isolates were only isolated from root samples.

FIGURE 1

Figure 1. Abundance of genera of bacterial isolates obtained from leaves and roots of Matthiola tricuspidata, Crithmum maritimuma, and Cakile maritima plants.

Apart from Bacillus, which was isolated from both roots and leaves, the rest of the genera were isolated only from roots or leaves (Supplementary Table 1). Similarly, all isolates from the genera Enterobacter, Pseudomonas, Microbacterium, Paenarthrobacter, Pantoea, Arthrobacter, Brachybacterium, Arenivirga, and Glutamicibacter were isolated from root samples, whereas the isolates Oceanobacillus, Planomicrobium, and Staphylococcus were isolated only from leaf samples (Supplementary Table 1).

The most frequently isolated genera were hosted in at least two of the three halophyte species. M. tricuspidata hosted the largest number of genera (18). Ca. maritima hosted the largest number of Enterobacter (13 out of 22) and the smallest number of Bacillus (2 out of 28). Members of the genera Brachybacterium and Arenivirga were isolated only from M. tricuspidata plants, whilst Glutamicibacter were isolated only from Cr. maritimum and the Planomicrobium from Ca. maritima.

Bacterial Growth Under Salinity Stress

Bacterial isolates were tested for their ability to grow in elevated NaCl concentrations (5, 7.5, 10, 15, and 17.5%) (Supplementary Figure 2A). Most isolates showed growth at 5% NaCl (96 isolates). From the 28 isolates from Ca. maritima roots, 26 isolates (92.9% of the total) showed ability to grow at 5% salinity, and 14 of these (50%) showed growth at 7.5% salinity; all four isolates from the leaves of the same plant showed growth at 10% salinity and two of these could grow at 17.5%. 21 out of 25 isolates (84%) of the roots of Cr. maritimum could grow at 5% salt and 16 (64%) could grow at 10% salinity. 12 out of 13 (92.3%) isolates from the leaves of Cr. maritimum could grow at the 5% level and eight (61.5%) could grow at 10% salt. From the 38 isolates obtained from the root of M. triscupidata, 27 (71%) could grow at the 5% salt threshold and three could grow in 10% salt.

Of the six isolates that managed to grow at 17.5% salinity, four were isolated from leaf tissues (Supplementary Table 1): Staphylococcus saprophyticus (CML12) and Oceanobacillus picturae (CML15) isolated from Ca. maritima leaves, Oceanobacillus picturae (CrL11) from Cr. maritimum leaves, and Micrococcus aloeverae (MTL04) from M. triscupidata leaves (Supplementary Table 1 and Supplementary Figure 2A). The two isolates from root tissues that could grow on 17.5% are Enterobacter hormaechei subsp. hoffmannii (CMR13) isolated from Ca. maritima and Bacillus hwajinpoensis (CrR23) isolated from Cr. maritimum (Supplementary Table 1 and Supplementary Figure 2A). Another three Bacilli isolates (CrR16: Bacillus haikouensis, CrR22: Bacillus haikouensis, MTR05: Terribacillus saccharophilus) showed growth on 15% salinity (Supplementary Table 1 and Supplementary Figure 2A).

Phytopathogens Growth Inhibition Ability

All bacterial isolates were subjected to in vitro inhibition assays against three known phytopathogens: the bacteria Ralstonia solanacearum and Clavibacter michiganensis subsp. michiganensis and the fungus Verticillium dahliae. In the assays against the phytopathogenic bacteria, the bacterial isolates that showed any kind of inhibition were characterized as having “weak,” “medium” or “strong” inhibitory activity based on the size of the inhibition zone around the bacterial colony (Figure 2). In the in vitro assay against Verticillium the inhibitory activity was similarly judged as “weak,” “medium” or “strong” based on the linear distance between the bacterial and the fungal colonies (Figure 2).

FIGURE 2

Figure 2. In vitro growth inhibition experiments: against phytopathogenic bacteria Ralstonia solanacearum or Clavibacter michiganensis (A), against phytopathogenic fungus Verticillium dahliae (B) and against human pathogenic fungus Aspergillus fumigatus (C). Representative experiments are shown to display their “weak” (A1–C1), “medium” (A2–C2), or “strong” (A3–C3) inhibitory activity.

Twenty-five (21.7%) out of 115 bacterial isolates demonstrated inhibition of the Ralstonia solanacearum growth (Supplementary Table 1). These 25 isolates belong to six genera: Bacillus, Enterobacter, Erwinia, Glutamicibacter, Paenarthrobacter, and Pseudomonas. Isolate CML04 (Bacillus altitudinis), obtained from leaf tissues of Ca. maritima, was the only leaf-derived isolate that showed antagonistic activity against all 3 tested phytopathogens (Supplementary Table 1). Of the 45 isolates isolated from M. triscupidata, three isolates showed a weak inhibitory zone against Ralstonia. A total of 10 isolates belonging to the genera Enterobacter and Pseudomonas (Supplementary Table 1) showed a strong inhibition (Figure 2).

A lower number of isolates showed any kind of inhibition against the phytopathogenic Clavibacter michiganensis subsp. michiganensis: 17 isolates out of 115 (14.8%). All 9 isolates from M. triscupidata (Bacillus licheniformis or Bacillus sonorensis isolates) demonstrated a strong inhibition zone (Figure 2). Similarly, two isolates from Cr. maritimum roots (both Pseudomonas glareae, CrR12 and CrR13) showed a strong inhibition zone (Supplementary Table 1). An additional two and four isolates showed weak and medium inhibition zone against Clavibacter, respectively (Supplementary Table 1).

The majority (76.5%) of the bacterial isolates demonstrated inhibition of the phytopathogenic fungus Verticillium dahliae (Supplementary Table 1). These isolates originate from both leaf and root tissues from all three halophytes. A strong inhibition zone was demonstrated by 34 of these 88 isolates, all of which except one, were isolated from halophytic plant roots (Supplementary Table 1). From these 34, 11 isolates belong to the genus Bacillus, another 11 to Pseudomonas and five to Enterobacter (Supplementary Table 1).

Eleven isolates were tested for inhibition against the human pathogenic fungus Aspergillus fumigatus. Interestingly, five isolates were able to inhibit the growth of A. fumigatus (Figure 2). Isolate CML04 (Bacillus altitudinis) isolated from Ca. maritima was the only isolate from leaf tissues able to show inhibitory effect, while the remaining isolates were isolated from M. triscupidata roots (Supplementary Table 1). Two isolates MTR17d (Bacillus sonorensis) and MTR17b (Bacillus licheniformis) showed a strong inhibition zone (Figure 2 and Supplementary Table 1).

In-planta Assay for Plant Growth Promotion and Salt Tolerance

Bacterial isolates with in vitro 10% and 17.5% NaCl salt tolerance were selected for the in-planta assays to demonstrate potential plant growth promotion under “no stress.” Arabidopsis thaliana plants were imbued with bacterial cultures and left for 7 days for the bacteria to adapt. Then, after watering the plants for a month, fresh and dry leaf weight were calculated (Table 1). The same experiment was repeated where after the 7 day mark, the plants were watered with 10 mL of 250 mM NaCl solution every 2–3 days for 30 days.

TABLE 1

Table 1. In-planta plant growth promotion and salt tolerance assays in Arabidopsis thaliana plants. For the salt tolerance assays plants were watered with and without NaCl solution every 2-3 days and fresh and dry plant weight was measured.

Under no stress conditions, plants inoculated with isolates CML12, CML15, CrL01, CrL11, CrR16, CrR23, MTR05 showed an increase in fresh leaf weight between 1.1 and 2.6 times to the non-inoculated plants and between 1.0 and 2.3 times increase in dry leaf weight (Table 1). Under salt stress, the growth promotion effect was less accentuated, since plants imbued with isolates CML15, CrR22, MTR05 had less increased fresh and dry leaf (Table 1).

Isolates CML15 and MTR05 conferred an increase in fresh and dry leaf weight both under no stress and under salt stress whereas isolate CrR22 had a positive affect only under salt stress condition (Table 1). On the other hand, isolates CML12, CrL01, CrL11, CrR16, and CrR23 had a positive effect on fresh and dry weight under no stress condition (Table 1).

Direct and Indirect in vitro Effects of Verticillium dahliae Growth

The selected 16 bacterial isolates with the exception of MTR17h inhibited significantly V. dahliae growth rate in dual-culture assays. However, only MTR17h, MTR17b, and MTR17c could suppress fungal growth by means of volatile compounds (Table 2). Likewise, nearly all isolates were capable of inhibiting fungal sporulation (except of MTR17g) in dual-culture assays. Most isolates significantly inhibited spore production in dual-plate assays. Interestingly, three isolates caused a significant induction of V. dahliae sporulation in such assays (MTR17h, MTR17b, and MTR17c), indicating that fungal growth suppression induces fungal sporulation (Table 2). Moreover, six out of 16 isolates could significantly reduce hyphae width in direct culture conditions, whereas seven out of 16 were capable of hyphae width reduction by the mean of volatiles. Additionally, nine isolates significantly inhibited microsclerotia formation in dual-culture assays; however, only three isolates significantly reduced microsclerotia formation in dual-plate assays (Table 2). MTR17h caused significant induction in microsclerotia formation both in dual-culture and in dual-plate assays.

TABLE 2

Table 2. Values of fungal parameters of Verticillium dahliae treated with 16 different bacterial isolates (CrR14, CrR18, CrR04, MTR12, MTR18, CM0R1, CMR03, CML04, CMR25, MTR17a, MTR17d, MTR17f, MTR17g, MTR17h, MTR17b, MTR17c) and Trichoderma harzianum strain T22 in dual-culture and dual-plate assays. Values were estimated as the percentage of inhibition compared to control (V.d.).

Suppression of Verticillium Wilt Symptoms in-planta

For the suppression of Verticillium dahliae wilt symptoms in-planta we used a well-established fungus/plant system, the Verticillium/eggplant system. We selected 16 bacterial isolates that showed promising in vitro growth inhibition effect to Verticillium.

Two distinct assays were performed (hereafter known as “experiment I” and “experiment II”). V. dahliae wilt symptoms on eggplant started 12 days after inoculation (d.p.i.), with V. dahliae conidial suspension and were recorded periodically for another 12 days in experiment I. Isolates CrR4, MTR12, MTR18, and CMR01 suppressed significantly disease severity at 18 and 21 d.p.i. whereas MTR18 and CM1 treatments caused significant reduction of disease severity at most observation time points (Table 2 and Figures 3, 4). Considering all disease parameters, CMR01 was the most effective isolate in terms of disease suppression (Table 2, Figure 3, and Supplementary Tables 2, 3). First disease symptoms in experiment II were also observed on 12 d.p.i. and recorded until 30 d.p.i. Disease severity progressed rapidly in the control (V.d.) and the non-suppressive treatments (MTR17d, MTR17f, and MTR17g), whereas MTR17a-, MTR17h-, MTR17b-, and MTR17c-treated plants showed less prominent symptoms and slower disease development (Table 2 and Supplementary Figure 3). Disease parameters indicated that isolate MTR17h, is comparable to the positive control (fungus Trichoderma harzianum isolate + TRIANUM-P), as the most effective in symptom suppression (Table 2 and Figure 4). While observed decrease in symptom severity in MTR17h-treated plants was associated with significantly lower V. dahliae re-isolation ratio compared to positive control (V.d.) plants, MTR17h isolate did not show strong growth inhibition effect on V. dahliae in in vitro assays (Figure 4), indicating less active growth of the pathogen into the xylem vessels. This finding could suggest that the plant innate immunity activation/reinforcement effect by MTR17h, needs to be further investigated in the future. Neither symptoms nor positive isolations were observed in negative control plants.

FIGURE 3

Figure 3. Verticillium wilt disease severity index on eggplant treated with various bacterial isolates at 12, 14, 18, 21, and 24 days post inoculation with Verticillium dahliae conidial suspension (20 mL of 5 × 10⁶ conidia mL^–1). Each column represents the mean of 21 plants after combining the results of 3 replicated experiments (experiment I). Columns at each observation time point followed by the same letter are not significantly different according to Tukey’s HSD test at P ≤ 0.05. Vertical bars indicate standard errors.

FIGURE 4

Figure 4. Verticillium wilt disease severity index on eggplant treated with various bacterial isolates and the commercial biofungicide TRIANUM-P (Koppert B.V. Hellas) at 12, 14, 16, 19, 23, 26, 28, and 30 days post inoculation with Verticillium dahliae conidial suspension (20 mL of 5 × 10⁶ conidia mL^–1). Each column represents the mean of 21 plants after combining the results of 3 replicated experiments (experiment II). Columns at each observation time point followed by the same letter are not significantly different according to Tukey’s HSD test at P ≤ 0.05. Vertical bars indicate standard errors.

Effects of Treatments in Plant Growth

Growth parameters of eggplant inoculated with V. dahliae and treated with the 16 isolates and the T. harzianum isolate TRIANUM-P or not (C−), are shown on Table 3. V. dahliae-inoculated plants treated with MTR17c and T. harzianum TRIANUM-P developed significantly higher fresh weight compared with the V. dahliae-inoculated controls, whereas most of the plant growth parameters in non-inoculated plants were significantly higher than the inoculated ones.

TABLE 3

Table 3. Values (± standard errors) of disease parameters for eggplants inoculated with V. dahliae and treated with different bacterial isolates and TRIANUM-P (CrR14, CrR18, CrR04, MTR12, MTR18, CMR01, CMR03, CML04, CMR25 in experiment I, and MTR17a, MTR17d, MTR17f, MTR17g, MTR17h, MTR17b, MTR17c, TRIANUM-P in experiment II) or not (C−, V.d.).

Whole-Genome Sequencing and Analysis of Selected Endophytic Bacterial Isolates

Whole-genome sequencing (WGS) was performed on 12 selected isolates. Genomes were annotated using RAST (Supplementary Figure 5). All genes related to the virulence, disease and defense that were predicted are presented in Table 4 and Supplementary Table 2.

TABLE 4

Table 4. Number of genes related to Virulence, Disease and Defense features of the three new bacterial species identified in this study. The genome analysis and the annotation was performed using the RAST genome annotation software.

Genome-wide average nucleotide identity (ANI) calculations against all bacterial genome assemblies in GenBank pointed to the existence of three previously unknown bacterial species; their genomes showed less than 95% ANI with any previously sequenced genomes. Phylogenetic analysis of the recA and gyrB gene sequences extracted from the genomes (Figures 5, 6) had indicated that two isolates belong to Pseudomonadaceae, while the third was a member of Arthrobacter genus. Isolates CMR25 and CMR27 belonged to an unidentified species of the P. putida group (Figure 6) and isolate CrR25 was an undefined species of the P. mendocina group (Figure 6). Consistent with these results, 16S rRNA gene sequences of CMR25 and CMR27 were 99.8% identical to that of Pseudomonas plecoglossicida and CrR25 was 98.75% identical to that of Pseudomonas benzenivorans. WGS analysis places CMR16 as an unidentified Arthrobacter species (Figure 5), while its 16S rRNA gene sequence assigns the isolate to Paenarthrobacter nitroguajacolicus (Kotoučková et al., 2004) with 98.78% identity. This species has been previously isolated from leaves of maize (Pisarska and Pietr, 2012) and promoted growth of wheat under salt stress (Safdarian et al., 2019). Unfortunately, no genome sequence is available for the type strain of this species; however, the ANI between CMR16 genome and previously sequenced genomes (Yao et al., 2015) of Paenarthrobacter nitroguajacolicus (Kotoučková et al., 2004) range between 86.08 and 86.76%, well below the widely used threshold of 96% for species membership.

FIGURE 5

Figure 5. Molecular phylogenetic analysis of Arthrobacter recA-gyrB genes by Maximum Likelihood method. Phylogenetic tree is drawn to scale, with branch lengths measured in the number of substitutions per site.

FIGURE 6

Figure 6. Molecular Phylogenetic analysis of the recA-gyrB genes from Pseudomonads belonging to P. putida and P. mendocina groups by Maximum Likelihood method. The tree with the highest log likelihood (–7831.6808) is shown. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site.

Discussion

Utilization of endophytic microorganisms for the control of biotic/abiotic stresses is a relatively unexplored area of research. Endophytes have been studied for over two decades (Saikkonen et al., 1998; Hasegawa et al., 2006; Kaul et al., 2016), however, our understanding about their role in plant defense against biotic/abiotic stresses is still limited (Liu et al., 2020; Pascale et al., 2020). Isolation, identification and the study of endophytes from plants that undergo continued abiotic stress could be essential for the development of proper biocontrol strategy for sustainable agriculture and food security.

Here, we investigated the abundance of taxa of the culturable bacterial endophytes of three halophytic plants, endemic in Crete island, Greece, using culture-dependent techniques (Figure 7). We also investigated the proof-of-concept of using the halophytes as a valuable source of beneficial microbes that can potentially be used in agriculture, by testing our initial hypothesis that these endophytes have plant growth promotion and biocontrol properties.

FIGURE 7

Figure 7. A graphical abstract describing the steps of the procedure we followed to accomplish this work.

Taxonomically, 24 different genera were identified, the three most abundant ones were Bacillus, Enterobacter, and Pseudomonas, all of which have been previously observed in studies of the endophytic microbiome of halophytes (Shabala, 2013; Qin et al., 2014; Mora-Ruiz et al., 2016; Yuan et al., 2016).

In-planta testing of Oceanobacillus picturae (CML15), Terribacillus saccharophilus (MTR05), and Bacillus haikouensis (CrR22) demonstrated an increase in both dry and fresh leaf weight in Arabidopsis thaliana plant under salinity stress. These isolates are promising biofertilizers, since other isolates of the same species have also been shown to have plant growth promotion properties; Terribacillus saccharophilus, firstly reported at 2007, is a known halophilic bacterium able to grow on 0–16% NaCl (An et al., 2007; Liu et al., 2010). This species is a known endophytic bacterium (Han et al., 2011), shown to trigger an increase on monoterpenes, sesquiterpenes, tocopherols, and membrane sterols, compounds engaged in antioxidant capacity in leaf tissues of grape resulting in stress tolerance (Salomon et al., 2016). Bacillus haikouensis is halotolerant bacterium isolated from paddy soil, able to grow on up to 17% NaCl (Li et al., 2014). Oceanobacillus picturae is a halophilic phosphate-solubilizing species with demonstrated siderophore production potential, isolated from saline environments and shown to promote plant growth in mangroves and confer salinity stress tolerance in barley (El-Tarabily and Youssef, 2010; Mapelli et al., 2013; Orhan and Demirci, 2020). Many of our isolates were able to grow at high concentrations of salt (5–17% NaCl).

Isolates belonging to the species Bacillus licheniformis, Bacillus sonorensis, Pseudomonas glareae, Enterobacter hormaechei, Pseudomonas benzenivorans, Pseudomonas monteilii, Pseudomonas plecoglossicida were shown to have strong antagonistic activity against the phytopathogenic bacteria Ralstonia solanacearum and Clavibacter michiganensis subsp. michiganensis, two very important plant pathogens with high economic impact on agriculture (Gartemann et al., 2003; Peeters et al., 2013). Both are very important phytopathogens, since Ralstonia has a large host range able to infect more than 200 plant species easily adaptable in varying environmental conditions whereas C. michiganensis subsp. michiganensis is able to infect wheat, maize, potatoes, and red and green peppers, despite its main host being tomatoes (Eichenlaub and Gartemann, 2011; Peeters et al., 2013; Hwang et al., 2018). Moreover, specific isolates with in vitro growth inhibition effect against V. dahliae, were tested for their ability to inhibit V. dahliae in-planta. Several isolates demonstrated an in-planta suppression effect of the polyphagous pathogen V. dahliae. Interestingly, isolates with strong in vitro effect did not manage to inhibit V. dahliae in-planta, but other isolates with medium or low in vitro effect inhibited in-planta V. dahliae growth strongly. These data provide the proof of concept for our study but also indicate that in future studies all resulting isolates need to be investigated for their in-planta antifungal and/or antibacterial growth inhibition capacity.

Furthermore, the whole-genome sequencing (WGS) of selected isolates revealed three new previously unidentified bacterial species. The identification of three new species in a very small number isolates indicates the high potential of the wild halophytic endophytome in terms of identifying new microbial species with novel capabilities, that could be beneficial for both agriculture (stress tolerance, growth promotion, etc.) and potentially in clinical practice (identification of new antibiotics, antifungal compounds, etc.).

The results from the study of the microbial collection we generated, could be the basis for the future development of various synthetic “bio-inoculants,” as the isolates possess all of the following attributes for such usage: (a) they are not pathogenic and do not induce plant disease; (b) are able to colonize plants, and (c) are culturable, so they can be used in modern agriculture. Furthermore, these isolates can be the basis for future studies, including the investigation of the colonization strategies that these microbes use, as well as, the elucidation of the molecular dialogs that take place during host-root colonization; the growth promotion; the salt tolerance and the immunity activation, by unique beneficial endophytes or artificial endophytic communities.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

PS designed the research. CC, GD, AC, EM, AS, GM, and GR performed the research. CC, IL, DS, VC, EM, and PS analyzed the data. CC, EM, DS, VC, and PS wrote the manuscript. All authors have read and approved the manuscript.

Funding

This project utilized equipment funded by the Wellcome Trust Institutional Strategic Support Fund (WT097835MF), the Wellcome Trust Multi-User Equipment Award (WT101650MA) and BBSRC LOLA award (BB/K003240/1). This work benefitted from the support of the University of Exeter’s High-Performance Computing (HPC) facility. The project was partially supported by the Emblematic Action of the Greek General Secretariat for Research and Technology, “Agro4Crete,” Protocol Number: SAE 013, Operational Program: SAE 013.

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 acknowledge Dr. Karen Moore, Audrey Farbos, Georgina Morris of the Exeter Sequencing Service at University of Exeter, for DNA library preparation and genome sequencing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.681567/full#supplementary-material

References

An, S.-Y., Asahara, M., Goto, K., Kasai, H., and Yokota, A. (2007). Terribacillus saccharophilus gen. nov., sp. nov. and Terribacillus halophilus sp. nov., spore-forming bacteria isolated from field soil in Japan. Int. J. Syst. Evol. Microbiol. 57, 51–55. doi: 10.1099/ijs.0.64340-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulgarelli, D., Garrido-Oter, R., Münch, P. C., Weiman, A., Dröge, J., Pan, Y., et al. (2015). Structure and function of the bacterial root microbiota in wild and domesticated barley. Cell Host Microb. 17, 392–403. doi: 10.1016/j.chom.2015.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulgarelli, D., Spaepen, S., Schlaeppi, K., Schulze-Lefert, P., and van Themaat, E. V. L. (2013). Structure and functions of the bacterial microbiota of plants. Annu. Rev. Plant Biol. 64, 807–838. doi: 10.1146/annurev-arplant-050312-120106

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, C., and Madden, L. (1990). Introduction to Plant Disease Epidemiology. Hoboken, NJ: John Wiley & Sons.

Google Scholar

Chaurasia, B., Pandey, A., Palni, L. M. S., Trivedi, P., Kumar, B., and Colvin, N. (2005). Diffusible and volatile compounds produced by an antagonistic Bacillus subtilis strain cause structural deformations in pathogenic fungi in vitro. Microbiol. Res. 160, 75–81. doi: 10.1016/j.micres.2004.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Dennis Benson, G. A., Cavanaugh, M., Clark, K., Karsch-Mizrachi, I., Ostell, J., Pruitt, K. D., et al. (2017). GenBank. Nucleic Acids Res. 46, D41–D47. doi: 10.1093/nar/gkx1094

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodd, I. C., and Pérez-Alfocea, F. (2012). Microbial amelioration of crop salinity stress. J. Exp. Bot. 63, 3415–3428. doi: 10.1093/jxb/ers033

PubMed Abstract | CrossRef Full Text | Google Scholar

Edwards, J., Johnson, C., Santos-Medellín, C., Lurie, E., Podishetty, N. K., Bhatnagar, S., et al. (2015). Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. U.S.A. 112, E911–E920. doi: 10.1073/pnas.1414592112

PubMed Abstract | CrossRef Full Text | Google Scholar

Eichenlaub, R., and Gartemann, K.-H. (2011). The Clavibacter michiganensis subspecies: molecular investigation of gram-positive bacterial plant pathogens. Annu. Rev. Phytopathol. 49, 445–464. doi: 10.1146/annurev-phyto-072910-095258

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Tarabily, K. A., and Youssef, T. (2010). Enhancement of morphological, anatomical and physiological characteristics of seedlings of the mangrove Avicennia marina inoculated with a native phosphate-solubilizing isolate of Oceanobacillus picturae under greenhouse conditions. Plant Soil 332, 147–162. doi: 10.1007/s11104-010-0280-y

CrossRef Full Text | Google Scholar

Etesami, H., and Beattie, G. A. (2018). Mining halophytes for plant growth-promoting halotolerant bacteria to enhance the salinity tolerance of non-halophytic crops. Front. Microbiol. 9:148. doi: 10.3389/fmicb.2018.00148

PubMed Abstract | CrossRef Full Text | Google Scholar

Felsenstein, J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783. doi: 10.2307/2408678

CrossRef Full Text | Google Scholar

Gartemann, K.-H., Kirchner, O., Engemann, J., Gräfen, I., Eichenlaub, R., and Burger, A. (2003). Clavibacter michiganensis subsp. michiganensis: first steps in the understanding of virulence of a Gram-positive phytopathogenic bacterium. J. Biotechnol. 106, 179–191. doi: 10.1016/j.jbiotec.2003.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Glick, B. R. (2014). Bacteria with ACC deaminase can promote plant growth and help to feed the world. Microbiol. Res. 169, 30–39. doi: 10.1016/j.micres.2013.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, J., Song, Y., Liu, Z., and Hu, Y. (2011). Culturable bacterial community analysis in the root domains of two varieties of tree peony (Paeonia ostii). FEMS Microbiol. Lett. 322, 15–24. doi: 10.1111/j.1574-6968.2011.02319.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Haney, C. H., Samuel, B. S., Bush, J., and Ausubel, F. M. (2015). Associations with rhizosphere bacteria can confer an adaptive advantage to plants. Nat. Plants 1:15051. doi: 10.1038/nplants.2015.51

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasegawa, S., Meguro, A., Shimizu, M., Nishimura, T., and Kunoh, H. (2006). Endophytic actinomycetes and their interactions with host plants. Actinomycetologica 20, 72–81. doi: 10.3209/saj.20.72

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, I. S., Oh, E.-J., Kim, D., and Oh, C.-S. (2018). Multiple plasmid-borne virulence genes of Clavibacter michiganensis ssp. capsici critical for disease development in pepper. New Phytol. 217, 1177–1189. doi: 10.1111/nph.14896

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaul, S., Sharma, T., and Dhar, K. M. (2016). “Omics” tools for better understanding the plant-endophyte interactions. Front. Plant Sci. 7:955. doi: 10.3389/fpls.2016.00955

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotoučková, L., Schumann, P., Durnová, E., Spröer, C., Sedláček, I., Neča, J., et al. (2004). Arthrobacter nitroguajacolicus sp. nov., a novel 4-nitroguaiacol-degrading actinobacterium. Int. J. Syst. Evol. Microbiol. 54, 773–777. doi: 10.1099/ijs.0.02923-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096

PubMed Abstract | CrossRef Full Text | Google Scholar

Lahlali, R., Bajii, M., and Jijakli, M. H. (2007). Isolation and evaluation of bacteria and fungi as biological control agents against Rhizoctonia solani. Commun. Agric. Appl. Biol. Sci. 72, 973–982.

Google Scholar

Lane, J. D. (1991). “16S/23S rRNA sequencing,” in Nucleic Acid Techniques in Bacterial Systematics, eds E. Stackebrandt and M. Goodfellow (Hoboken, NJ: John Wiley & Sons Ltd), 115–175.

Google Scholar

Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., et al. (2007). Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948. doi: 10.1093/bioinformatics/btm404

PubMed Abstract | CrossRef Full Text | Google Scholar

Leinonen, R., Sugawara, H., and Shumway, M. (2011). The sequence read archive. Nucleic Acids Res. 39, D19–D21. doi: 10.1093/nar/gkq1019

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Yang, G., Lu, Q., Zhao, Y., and Zhou, S. (2014). Bacillus haikouensis sp. nov., a facultatively anaerobic halotolerant bacterium isolated from a paddy soil. Antonie van Leeuwenhoek Int. J. Gen. Mol. Microbiol. 106, 789–794. doi: 10.1007/s10482-014-0248-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Brettell, L. E., Qiu, Z., and Singh, B. K. (2020). Microbiome-mediated stress resistance in plants. Trends Plant Sci. 25, 733–743. doi: 10.1016/j.tplants.2020.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Carvalhais, L. C., Crawford, M., Singh, E., Dennis, P. G., Pieterse, C. M. J., et al. (2017). Inner plant values: diversity, colonization and benefits from endophytic bacteria. Front. Microbiol. 8:2552. doi: 10.3389/fmicb.2017.02552

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Jiang, L., Guo, C., and Yang, S. S. (2010). Terribacillus aidingensis sp. nov., a moderately Halophilic bacterium. Int. J. Syst. Evol. Microbiol. 60, 2940–2945. doi: 10.1099/ijs.0.017228-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mapelli, F., Marasco, R., Rolli, E., Barbato, M., Cherif, H., Guesmi, A., et al. (2013). Potential for plant growth promotion of Rhizobacteria associated with salicornia growing in tunisian hypersaline soils. Biomed Res. Int. 2013, 1–13. doi: 10.1155/2013/248078

PubMed Abstract | CrossRef Full Text | Google Scholar

Markakis, E. A., Fountoulakis, M. S., Daskalakis, G. C., Kokkinis, M., and Ligoxigakis, E. K. (2016). The suppressive effect of compost amendments on Fusarium oxysporum f.sp. radicis-cucumerinum in cucumber and Verticillium dahliae in eggplant. Crop Prot. 79, 70–79. doi: 10.1016/j.cropro.2015.10.015

CrossRef Full Text | Google Scholar

Marques, J. M., da Silva, T. F., Vollú, R. E., de Lacerda, J. R. M., Blank, A. F., Smalla, K., et al. (2015). Bacterial endophytes of sweet potato tuberous roots affected by the plant genotype and growth stage. Appl. Soil Ecol. 96, 273–281. doi: 10.1016/j.apsoil.2015.08.020

CrossRef Full Text | Google Scholar

Mendes, R., Garbeva, P., and Raaijmakers, J. M. (2013). The rhizosphere microbiome: Significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663. doi: 10.1111/1574-6976.12028

PubMed Abstract | CrossRef Full Text | Google Scholar

Mora-Ruiz, M. D. R., Font-Verdera, F., Orfila, A., Rita, J., and Rosselló-Móra, R. (2016). Endophytic microbial diversity of the halophyte Arthrocnemum macrostachyum across plant compartments. FEMS Microbiol. Ecol. 92:fiw145. doi: 10.1093/femsec/fiw145

PubMed Abstract | CrossRef Full Text | Google Scholar

Munns, R. (2002). Comparative physiology of salt and water stress. Plant. Cell Environ. 25, 239–250. doi: 10.1046/j.0016-8025.2001.00808.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Orhan, F., and Demirci, A. (2020). Salt stress mitigating potential of Halotolerant/Halophilic plant growth promoting. Geomicrobiol. J. 37, 663–669. doi: 10.1080/01490451.2020.1761911

CrossRef Full Text | Google Scholar

Overbeek, R., Olson, R., Pusch, G. D., Olsen, G. J., Davis, J. J., Disz, T., et al. (2014). The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 42, D206–D214. doi: 10.1093/nar/gkt1226

PubMed Abstract | CrossRef Full Text | Google Scholar

Pascale, A., Proietti, S., Pantelides, I. S., and Stringlis, I. A. (2020). Modulation of the root microbiome by plant molecules: the basis for targeted disease suppression and plant growth promotion. Front. Plant Sci. 10:1741. doi: 10.3389/fpls.2019.01741

PubMed Abstract | CrossRef Full Text | Google Scholar

Peeters, N., Guidot, A., Vailleau, F., and Valls, M. (2013). Ralstonia solanacearum, a widespread bacterial plant pathogen in the post-genomic era. Mol. Plant Pathol. 14, 651–662. doi: 10.1111/mpp.12038

PubMed Abstract | CrossRef Full Text | Google Scholar

Pegg, G. G. F., and Brady, B. L. B. (2002). Verticillium Wilts. New York, NY: CABI.

Google Scholar

Pérez-Jaramillo, J. E., Mendes, R., and Raaijmakers, J. M. (2016). Impact of plant domestication on rhizosphere microbiome assembly and functions. Plant Mol. Biol. 90, 635–644. doi: 10.1007/s11103-015-0337-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Pisarska, K., and Pietr, S. J. (2012). Isolation and partial characterization of culturable endophytic Arthrobacter spp. from leaves of maize (Zea mays L.). Commun. Agric. Appl. Biol. Sci. 77, 225–233.

Google Scholar

Qin, S., Zhang, Y. J., Yuan, B., Xu, P. Y., Xing, K., Wang, J., et al. (2014). Isolation of ACC deaminase-producing habitat-adapted symbiotic bacteria associated with halophyte Limonium sinense (Girard) Kuntze and evaluating their plant growth-promoting activity under salt stress. Plant Soil 374, 753–766. doi: 10.1007/s11104-013-1918-3

CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez, P. A., Rothballer, M., Chowdhury, S. P., Nussbaumer, T., Gutjahr, C., and Falter-Braun, P. (2019). Systems biology of plant-microbiome interactions. Mol. Plant 12, 804–821. doi: 10.1016/j.molp.2019.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodríguez-Blanco, A., Sicardi, M., and Frioni, L. (2015). Plant genotype and nitrogen fertilization effects on abundance and diversity of diazotrophic bacteria associated with maize (Zea mays L.). Biol. Fertil. Soils 51, 391–402. doi: 10.1007/s00374-014-0986-8

CrossRef Full Text | Google Scholar

Ruppel, S., Franken, P., and Witzel, K. (2013). Properties of the halophyte microbiome and their implications for plant salt tolerance. Funct. Plant Biol. 40:940. doi: 10.1071/FP12355

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryan, R. P., Germaine, K., Franks, A., Ryan, D. J., and Dowling, D. N. (2008). Bacterial endophytes: recent developments and applications. FEMS Microbiol. Lett. 278, 1–9. doi: 10.1111/j.1574-6968.2007.00918.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Safdarian, M., Askari, H., Shariati, J. V., and Nematzadeh, G. (2019). Transcriptional responses of wheat roots inoculated with Arthrobacter nitroguajacolicus to salt stress. Sci. Rep. 9:1792. doi: 10.1038/s41598-018-38398-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Saikkonen, K., Faeth, S. H., Helander, M., and Sullivan, T. J. (1998). FUNGAL ENDOPHYTES: a continuum of interactions with host plants. Annu. Rev. Ecol. Syst. 29, 319–343. doi: 10.1146/annurev.ecolsys.29.1.319

CrossRef Full Text | Google Scholar

Salomon, M. V., Purpora, R., Bottini, R., and Piccoli, P. (2016). Rhizosphere associated bacteria trigger accumulation of terpenes in leaves of Vitis vinifera L. cv. Malbec that protect cells against reactive oxygen species. Plant Physiol. Biochem. 106, 295–304. doi: 10.1016/j.plaphy.2016.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Shabala, S. (2013). Learning from halophytes: physiological basis and strategies to improve abiotic stress tolerance in crops. Ann. Bot. 112, 1209–1221. doi: 10.1093/aob/mct205

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., and Nei, M. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 10, 512–526. doi: 10.1093/oxfordjournals.molbev.a040023

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121

PubMed Abstract | CrossRef Full Text | Google Scholar

Tkacz, A., and Poole, P. (2015). Role of root microbiota in plant productivity. J. Exp. Bot. 66, 2167–2175. doi: 10.1093/jxb/erv157

PubMed Abstract | CrossRef Full Text | Google Scholar

Toju, H., Peay, K. G., Yamamichi, M., Narisawa, K., Hiruma, K., Naito, K., et al. (2018). Core microbiomes for sustainable agroecosystems. Nat. Plants 4, 247–257. doi: 10.1038/s41477-018-0139-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaishnav, A., Shukla, A. K., Sharma, A., Kumar, R., and Choudhary, D. K. (2019). Endophytic bacteria in plant salt stress tolerance: current and future prospects. J. Plant Growth Regul. 38, 650–668. doi: 10.1007/s00344-018-9880-1

CrossRef Full Text | Google Scholar

White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in PCR Protocols, eds M. A. Innis, D. H. Gelfand, J. J. Sninsky, and T. J. White (San Diego: Academic Press), 315–322. doi: 10.1016/B978-0-12-372180-8.50042-1

CrossRef Full Text | Google Scholar

Yao, Y., Tang, H., Su, F., and Xu, P. (2015). Comparative genome analysis reveals the molecular basis of nicotine degradation and survival capacities of Arthrobacter. Sci. Rep. 5:8642. doi: 10.1038/srep08642

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, S.-H., Ha, S.-M., Kwon, S., Lim, J., Kim, Y., Seo, H., et al. (2017). Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67, 1613–1617. doi: 10.1099/ijsem.0.001755

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, Z., Druzhinina, I. S., Labbé, J., Redman, R., Qin, Y., Rodriguez, R., et al. (2016). Specialized microbiome of a halophyte and its role in helping non-host plants to withstand salinity. Sci. Rep. 6:32467. doi: 10.1038/srep32467

PubMed Abstract | CrossRef Full Text | Google Scholar