Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 29 August 2025

Sec. Microbe and Virus Interactions with Plants

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1643204

Genomic and phenotypic profiling of an Artemisia endophyte: dual biostimulant and biocontrol activities in pea under Fusarium stress

  • 1Microbiology Laboratory, Institute of Agriculture, Lithuanian Research Centre for Agriculture and Forestry, Kėdainiai, Lithuania
  • 2Department of Plant Science, McGill University, Montreal, QC, Canada

Objective: To evaluate the plant growth-promoting and disease-suppressing potential of Serratia sp. AR11, an endophytic bacterium isolated from Artemisia absinthium, through phenotypic assessment and whole-genome analysis in pea (Pisum sativum) under normal and Fusarium-stress conditions.

Materials and methods: Greenhouse experiments were conducted to assess the effects of AR11 inoculation on shoot and root growth, biomass, chlorophyll content, and Fusarium-induced stunting. Whole-genome sequencing was performed using the PacBio SMRT platform, followed by functional annotation to identify genes related to nutrient metabolism, secondary metabolite biosynthesis, and stress adaptation. Biosafety assessment included screening for virulence and antibiotic resistance genes.

Results: AR11 inoculation significantly increased shoot and root length and biomass, while reducing Fusarium-induced stunting by over 70%. Under pathogen stress, treated plants showed a 67% increase in SPAD index compared to controls. Genome analysis revealed a complete 5.49 Mb circular genome with 5,175 protein-coding genes, including those for nitrogen metabolism, phosphate solubilization, siderophore production, and antifungal secondary metabolite biosynthesis. Additional features included stress-responsive regulatory systems, CRISPR-Cas loci, and DNA methylation patterns. No high-risk virulence or acquired antibiotic resistance genes were detected.

Conclusion: Serratia sp. AR11 is a safe, multifunctional endophytic bacterium with strong genomic and phenotypic evidence for enhancing pea growth and suppressing Fusarium-related diseases, making it a promising candidate for sustainable low-input agriculture.

Introduction

Promoting sustainable agriculture requires innovative solutions that minimize chemical inputs while enhancing crop productivity and resilience (Sudheer et al., 2020). Endophytic bacteria are a key part of the plant microbiome. Most live inside plant tissues without causing harm. These bacteria can improve plant growth, increase stress tolerance, and help protect against disease (Anand et al., 2023). These bacteria can fix atmospheric nitrogen, solubilize phosphate, produce plant hormones, and synthesize antimicrobial metabolites traits essential for biostimulant and bioprotective applications (Timofeeva et al., 2023). However, identifying highly effective strains and understanding the molecular mechanisms behind their beneficial effects remain major challenges.

Medicinal plants are an underexplored source of endophytic diversity with high biosynthetic potential (Aswani et al., 2020). Among them, Artemisia species have drawn attention to their rich secondary metabolite profiles and ethnobotanical importance (Lantzouraki et al., 2023). Artemisia species, such as A. absinthium, A. alba, A. annua, A. verlotiorum, and A. vulgaris, are recognized for their secondary metabolic. Caffeoyl and feruloyl quinic acid esters are identified as their potential bioactive compounds, which notably characterize the most antioxidant-rich species (Negri et al., 2024). Host plants and their endophytes often develop specialized relationships, a result of long-term co-evolution. These interactions can influence both the plant’s metabolic profile and the functional traits of the endophytes (Jia et al., 2016). The chemically rich environment of medicinal plants such as Artemisia may shape endophyte communities capable of producing diverse secondary metabolites (Tripathi et al., 2025). In some cases, the therapeutic properties of the plant may result, at least in part from metabolites synthesized by its endophytes. For example, Artemisia nilagirica and A. absenithium host endophytic bacteria known to produce antibacterial compounds (Ashitha et al., 2019; Damavandi et al., 2023). Endophytes, particularly those found in medicinal plants with ethnobotanical significance, act as promising sources of novel bioactive compounds (Nongkhlaw and Joshi, 2015). They produce a range of compounds that support plant growth and enhance protection against pathogens (Basit et al., 2021). Endophytes from Artemisia spp. not only contribute to host metabolism but may also enhance plant health in other crops (Hadian et al., 2024). Despite increasing recognition of their bioactive potential, few studies have investigated the functional genomics of Artemisia derived endophytes, particularly in relation to agriculturally important crops like pea (Pisum sativum L.).

Pea is a globally important legume that contributes to soil fertility through biological nitrogen fixation (Dhillon et al., 2022). However, it is highly susceptible to root rot caused by Fusarium spp., leading to yield losses exceeding 90% in severe cases (El-Saadony et al., 2021). While chemical fungicides remain a common control strategy, they pose environmental risks and are often ineffective under field variability. Harnessing endophytic microbes as dual function bioinoculants enhances growth while suppressing pathogens represents a promising but underdeveloped strategy in legume systems.

Endophytic microorganisms can suppress invading pathogens through multiple mechanisms, including competition for nutrients, and ecological niches within the rhizosphere, rhizoplane, and internal plant tissues. Certain beneficial endophytic bacteria produce antibiotics and lytic enzymes that actively inhibit pathogen growth. Additionally, introducing endophytic microbes into plants through artificial inoculation can significantly reduce infections caused by fungi, bacteria, viruses, insects, and nematodes, serving as an effective biological control strategy (Fontana et al., 2021; Hallmann et al., 1997). Numerous studies highlight the potential of endophytic bacteria in promoting plant growth and protecting against Fusarium infections. For example, it has been reported that Pseudomonas aeruginosa and Aneurinibacillus aneurinilyticus isolated from Ocimum sanctum significantly enhances pea growth and provides strong protection against Fusarium oxysporum f.sp. pisi, a major root rot pathogen (Gupta et al., 2022).

In our previous work, we isolated 61 endophytic bacterial strains from the roots, stems, and leaves of four Artemisia species (A. absinthium, A. campestris, A. dubia, and A. vulgaris) collected across diverse agroecological zones. Initial in vitro screening for key plant growth-promoting traits, such as phosphate solubilization, nitrogen fixation, and seed germination enhancement in pea (Pisum sativum) led to the selection of seven promising strains for further study (Hadian et al., 2024). Among these, strains AR11 and AR32, initially classified as Bacillus spp. via 16S rRNA gene sequencing, showed potential ecological relevance by modulating microbial diversity in the pea rhizosphere (Hadian et al., 2025). Despite growing recognition of the potential of endophytic bacteria in sustainable agriculture, the ecological functions and genomic characteristics of Artemisia-associated endophytes remain largely underexplored, particularly within legume cropping systems. While traits such as nitrogen fixation and phosphate solubilization have been individually reported, integrated studies that link genomic insights with plant performance under biotic stress are still limited. This is especially relevant for managing Fusarium wilt, a major challenge in pea (Pisum sativum) cultivation worldwide.

While many studies stop at phenotypic evaluation, our research integrates functional screening with whole-genome analysis to determine the molecular basis of beneficial traits. This phenotypic-to-genomic approach remains underexplored for endophytes from medicinal plants like Artemisia, particularly in pea systems. To address this, the present study combines phenotypic evaluation of seven Artemisia-derived endophytic bacterial strains with genome-based analysis under controlled greenhouse conditions, both with and without Fusarium inoculation. This approach enabled the identification of strains with consistent growth-promoting and disease-suppressing effects. The most effective strain was further investigated at the genomic level, revealing key traits related to plant growth support, stress response, and antifungal potential. The results provide valuable insights into the functional capacity of Artemisia-associated endophytes and support their potential use as bioinoculants in pea production.

We hypothesize that endophytic bacteria from Artemisia promote pea growth and reduce Fusarium infection through specific molecular traits. These include genes for nitrogen metabolism, phosphate solubilization, and antifungal compound production. By linking plant performance with genomic data, we aim to identify the key mechanisms behind this dual functionality.

Materials and methods

Isolation and selection of endophytic bacteria

Endophytic bacterial strains used in this study were previously isolated from root, stem, and leaf tissues of four Artemisia species (A. absinthium, A. campestris, A. dubia, and A. vulgaris), collected from three different regions in Lithuania: Kaunas, Kėdainiai, and Šiauliai. A total of 61 isolates were obtained, and seven strains were selected based on beneficial traits, including phosphate solubilization, nitrogen fixation, enhancement of pea seed germination, and antagonistic activity against Fusarium spp., as assessed by dual-culture assays. These results have been previously published (Hadian et al., 2024). The selected strains were: AR11, AR32, and AR35 (from A. absinthium root); VL13 (from A. vulgaris leaf); VR24 (from A. vulgaris root); CR25 (from A. campestris root); and DR31 (from A. dubia root).

Assessment of plant growth-promoting activity and biocontrol potential in pea

Seed inoculation with endophytic bacteria

Pea seeds (Pisum sativum L.) were surface sterilized following a standard protocol. Sterilized seeds were immersed in bacterial suspensions (OD₆₀₀ = 1.0, approximately 108 CFU mL−1) and incubated for 3 h at room temperature with gentle shaking to ensure uniform coating. Control seeds underwent the same procedure but were treated with sterile distilled water.

Soil preparation

Soil was mixed with commercial compost and sieved through a 2 mm mesh to remove debris. The final soil mix had an organic carbon content of 11.99%, total nitrogen of 1.084%, and a slightly alkaline pH of 7.79, which is favorable for nutrient availability and microbial activity. Macronutrient analysis showed 460 mg kg−1 of potassium oxide (K₂O) and 205 mg kg−1 of phosphorus pentoxide (P₂O₅), indicating high fertility suitable for plant growth and microbial interaction studies.

Preparation of pathogenic inoculum

The pathogenic Fusarium sp. used in this study was previously isolated from infected pea roots and identified at the Microbiology Department of LAMMC (Hadian et al., 2024). Molecular analysis showed high similarity to Fusarium oxysporum, though species-level resolution was not the primary focus. Pathogenicity tests confirming its virulence on pea were thoroughly performed as described in Hadian et al. (2024), ensuring its suitability for the present experiments. Mycelial plugs (~5 mm) from PDA stock cultures were transferred to fresh potato dextrose agar (PDA) plates and incubated at 25 °C for 4 days. After 1 week of growth, conidia were harvested by flooding the plates with sterile distilled water and gently scraping the surface. The resulting suspension was filtered through sterile muslin cloth, and spore concentration was adjusted to 1 × 104 conidia mL−1 using a hemocytometer. The inoculum was applied in two stages to ensure effective root-zone colonization: the first application was conducted at the first true leaf stage, and the second 1 week later by injecting 5 mL of the spore suspension into the soil near the root zone. Control plants received an equal volume of sterile distilled water (Raza et al., 2024).

Experimental design and pot trials

The experiment was conducted under controlled greenhouse conditions (25 ± 2 °C, 16 h photoperiod, light intensity of ~250 μmol m−2 s−1, 60% relative humidity), where pea seeds were subjected to various treatments and sown in pots containing sterilized soil. The pots were maintained under optimal environmental conditions with regular irrigation to ensure adequate soil moisture. The experimental design included several treatment groups. The control group consisted of water-primed seeds without any bacterial or fungal inoculation. Another group received bacterial inoculation only, in which seeds were treated with one of the selected bacterial strains (AR11, AR32, AR35, VL13, VR24, CR25, and DR31) without exposure to Fusarium. A separate group, designated as Fusarium inoculation only (Control + Fu), included seeds that were challenged solely with the Fusarium pathogen, without bacterial treatment. The final group involved co-inoculation, where seeds were treated with a bacterial strain and simultaneously challenged with Fusarium. Each treatment was replicated 4 times, and the plants were observed for a period of 14 days following pathogen inoculation. Fourteen days after Fusarium inoculation, the plants were carefully uprooted and washed with tap water to remove soil residues. The root systems were then visually examined and scored for disease severity using a standardized 1–5 scale, where a score of 1 indicated healthy roots with no visible symptoms, 2 represented up to 25% root tissue damage, 3 indicated approximately 50% damage, 4 corresponded to around 75% root tissue damage, and 5 denoted complete damage or decay of root tissues.

Measurement of plant growth parameters

Morphological parameters including root length, shoot length, fresh and dry weight of both roots and shoots were measured for each plant at the time of harvest. These values were used to assess the overall impact of treatments on plant growth and development.

Chlorophyll content estimation

Chlorophyll content was determined using a SPAD-502 Plus chlorophyll meter (Konica Minolta, Japan). SPAD readings were taken from three fully expanded, healthy leaves, one each from the upper, middle, and lower canopy, to minimize variability. The meaning of the three readings was calculated to obtain the final SPAD value for each replicate.

DNA extraction

Genomic DNA was extracted from the endophytic bacterial strain Serratia sp. AR11 using the GeneJET Genomic DNA Purification Kit (Thermo Scientific, United States) following the manufacturer’s protocol. Briefly, single colonies were inoculated into Luria-Bertani (LB) broth and incubated overnight at 30 °C with shaking at 200 rpm. Cells were harvested by centrifugation at 10,000 × g for 5 min. Lysis and protein precipitation steps were followed by ethanol-based DNA purification. The quantity and quality of the extracted genomic DNA were assessed using a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific, United States), and DNA was then submitted for high-throughput whole-genome sequencing.

Whole-genome sequencing and bioinformatics analysis

High-molecular-weight genomic DNA from strain AR11 was submitted to Novogene (Germany) for long-read sequencing on the PacBio Sequel II platform using Circular Consensus Sequencing (CCS) mode to generate high-fidelity (HiFi) reads. Raw data underwent quality control and filtering prior to downstream analysis. De novo genome assembly was performed using Canu v1.5, optimized for long-read sequencing data. Circularization of contigs was completed using Circlator v1.5.5, and genome completeness was evaluated using BUSCO v5.4.3 with the bacterial single copy ortholog database. Structural genome annotation included gene prediction with Prodigal v2.6.3 (Hyatt et al., 2010), rRNA identification with Infernal v1.1.3 (Nawrocki and Eddy, 2013) using Rfam covariance models, and tRNA detection using tRNAscan-SE v2.0 (Chan and Lowe, 2019). Interspersed and tandem repetitive elements were annotated with RepeatMasker v4.0.5 (Tarailo-Graovac and Chen, 2009). Mobile genetic elements were identified using multiple tools: PhiSpy v2.3 (Akhter et al., 2012) for prophage prediction, IslandPath-DIMOB v0.2 (Bertelli and Brinkman, 2018) for genomic island detection, and CRT v1.2 (Bland et al., 2007) for CRISPR array identification. Secondary metabolite biosynthetic gene clusters (BGCs) were predicted with antiSMASH v5.0.0 (Blin et al., 2019). Functional annotation of predicted coding sequences was conducted using BLASTp (Altschul et al., 1990) (E-value ≤1 × 10−5) against the NCBI NR, Swiss-Prot, KEGG, and eggNOG databases. Gene Ontology (GO) terms were assigned via Blast2GO (Conesa et al., 2005), and functional classification into Clusters of Orthologous Groups (COGs) was performed using eggNOG-mapper. Carbohydrate-active enzymes (CAZymes) were identified using the CAZy (Cantarel et al., 2009) database, and membrane transporters were classified via the Transporter Classification Database (TCDB) (Saier, 2006). SignalP (Petersen et al., 2011) and TMHMM (Krogh et al., 2001) were used to predict signal peptides and transmembrane domains, respectively, to infer subcellular localization and secretion potential. Finally, DNA methylation profiling was carried out using the raw PacBio sequencing signal via SMRT Link v8.0.0, enabling genome-wide detection of modified bases including N6-methyladenine (m6A) and N4-methylcytosine (m4C), as well as associated methylation motifs and potential regulatory functions.

Statistical and data analysis

All statistical analyses were conducted using R software (v4.2.3; https://cran.r-project.org/; accessed on 10 January 2024). Prior to hypothesis testing, data were evaluated for normality using the Shapiro–Wilk test and for homogeneity of variance using Levene’s test. For parametric variables (shoot height, root length, biomass, SPAD index), one-way analysis of variance (ANOVA) was performed, followed by Tukey’s honest significant difference (HSD) test for post hoc comparisons at a significance level of p < 0.05. For non-parametric data, such as disease severity scores, the Kruskal–Wallis test was applied, followed by Dunn’s post hoc test with Benjamini–Hochberg correction for multiple comparisons. To assess potential associations between physiological parameters, Pearson’s correlation analysis was performed, with significance set at p < 0.05. All values are reported as means ± standard error (SE). Each treatment group included four biological replicates (n = 4). Graphical visualizations were generated using the ggplot2 package in R.

Results

Impact of Artemisia derived endophytic bacteria on growth and biomass of pea plants under Fusarium stress

Seven bacterial strains were isolated from the endosphere of Artemisia spp. and selected based on previously demonstrated nitrogen fixation and phosphate solubilization activity (Hadian et al., 2024). Taxonomic analysis via 16S rRNA gene sequencing classified six isolates as Bacillus spp., while one strain (CR25) was identified as Pseudomonas.

Growth promotion under non-stress conditions

Under optimal (non-stress) conditions, inoculation with endophytic bacteria significantly enhanced pea plant growth compared to the uninoculated control (p < 0.05; Figure 1). Among all treatments, AR11 exhibited the most pronounced growth-promoting effect, with plants reaching 24.27 ± 0.5 cm root and 28.75 ± 0.7 cm shoot height. These values were assigned to the highest Tukey HSD groupings (“a”), confirming a statistically significant improvement (p < 0.05) over all other treatments. AR32 and AR35 also significantly increased root and shoot heights compared to the control (p < 0.05). Plants treated with AR32 reached 20.30 ± 0.6 cm root and 25.75 ± 0.6 cm shoot height, placing them in groups “b” and “b,” respectively. AR35 induced root and shoot heights of 19.75 ± 0.4 cm and 22.50 ± 0.6 cm, grouped as “bc” for both parameters. In contrast, control plants reached only 17.28 ± 0.4 cm root (group “cd”) and 22.30 ± 0.5 cm shoot (group “bc”), confirming a significant difference from AR11 and AR32 (p < 0.05). CR25 showed intermediate performance (19.40 ± 0.6 cm root, group “bc”; 21.90 ± 0.6 cm shoot, group “bcd”), which was not significantly different from AR35 but still improved over the control in root height (p < 0.05). In contrast, the remaining strains (VR24, VL13, and DR31) demonstrated only moderate or limited growth promotion. VR24-inoculated plants exhibited 15.05 ± 0.6 cm root (group “fg”) and 18.50 ± 0.6 cm shoot (group “de”), both significantly lower than AR11, AR32, and AR35 (p < 0.05). VL13 showed similar performance (14.50 ± 0.5 cm root, group “fgh”; 18.15 ± 0.6 cm shoot, group “de”), while DR31 achieved 12.00 ± 0.4 cm root (group “ghi”) and 17.25 ± 0.4 cm shoot (group “ef”), placing it among the least effective strains.

Figure 1
Bar chart titled

Figure 1. Effect of endophytic bacterial inoculation on pea root and shoot lengths under Fusarium stress. Bars show mean ± SE (n = 4). Significant differences (p < 0.05, one-way ANOVA, Tukey’s HSD). Different letters above bars indicate significant differences.

Growth performance under Fusarium stress

Fusarium infection significantly reduced plant growth, as evidenced by the non-inoculated infected control, which showed root and shoot lengths of 8.75 ± 0.3 cm and 13.00 ± 0.4 cm, respectively. These values fell into Tukey groups “hi” (root) and “ef” (shoot), confirming a strong negative effect of the pathogen on pea development (p < 0.05; Figure 1). In contrast, inoculation with AR11 significantly mitigated the effects of Fusarium stress, with treated plants reaching 15.50 ± 0.4 cm root height and 21.00 ± 0.5 cm shoot height, grouped as “de” and “cd,” respectively. These values were the highest among all Fusarium-inoculated treatments, indicating that AR11 provided the most effective protection against pathogen-induced growth suppression (p < 0.05). AR32 and AR35 also conferred notable, though comparatively less, protective effects. AR32-treated plants reached 13.00 ± 0.3 cm root (group “ef”) and 19.25 ± 0.5 cm shoot (group “de”), while AR35 recorded 12.00 ± 0.4 cm root (group “f”) and 18.75 ± 0.6 cm shoot (group “de”). Though both strains significantly improved growth compared to the infected control (p < 0.05), they remained statistically inferior to AR11, as confirmed by non-overlapping Tukey groupings.

Biomass accumulation

Biomass data reinforced the observed trends in plant height (Figure 2). Under non-stress conditions, AR11 led to the highest shoot (2.88 ± 0.3 g fresh, 0.92 ± 0.05 g dry) and root (0.75 ± 0.2 g fresh, 0.11 ± 0.04 g dry) biomass, with all values falling within the top statistical group (“a”). In contrast, AR32 and AR35 resulted in significantly lower values (e.g., AR32: 2.37 ± 0.3 g shoot fresh, group “b”; 0.76 ± 0.05 g dry, group “b”). Under Fusarium challenge, AR11 maintained robust biomass accumulation, with shoot fresh and dry weights of 2.46 ± 0.2 g (group “d”) and 0.78 ± 0.03 g (group “d”), and root fresh and dry weights of 0.49 ± 0.2 g and 0.09 ± 0.02 g (group “d”), respectively. These values were significantly higher than those observed in the pathogen-only control (shoot dry weight: 0.49 ± 0.02 g; group “f”). Compared to AR11, AR32 and AR35 retained moderate biomass under stress but were statistically inferior. AR32-treated plants under infection conditions recorded 0.55 ± 0.03 g (shoot dry, group “f”) and 0.065 ± 0.03 g (root dry, group “e”), whereas AR35 showed even lower values (shoot: 0.47 ± 0.02 g, group “g”; root: 0.056 ± 0.02 g, group “f”). These distinctions confirm that AR11 consistently sustained significantly higher biomass accumulation under pathogen pressure.

Figure 2
Bar graphs comparing weight parameters across different strains for dry root weight, dry shoot weight, fresh root weight, and fresh shoot weight. The charts show varying mean weights with error bars for each strain. Strains are labeled along the x-axis, with mean weight in grams on the y-axis. Color coding is used to differentiate between weight types: orange for dry root weight, red for dry shoot weight, blue for fresh root weight, and green for fresh shoot weight. Letters above bars indicate statistical groupings.

Figure 2. Effect of endophytic bacterial inoculation on pea root and shoot biomass (fresh and dry weight) under both non-stress and Fusarium-stressed conditions. Bars show mean ± SE (n = 4). Treatments include individual strains, co-inoculations with Fusarium (+Fu), and uninoculated controls. Different lowercase letters above bars indicate significant differences among means (p < 0.05).

Effect of endophytic bacteria on chlorophyll content (SPAD value) under Fusarium challenge

Chlorophyll content, measured by SPAD index, significantly varied among treatments under both normal and Fusarium-infected conditions (Figure 3). In the absence of the pathogen, inoculation with strains AR11, AR32, and CR25 resulted in significantly higher SPAD values compared to the non-inoculated control (p < 0.05), indicating improved chlorophyll accumulation. Among these, AR11 exhibited the highest median SPAD value, followed by AR32 and CR25, as shown in the boxplot. Under Fusarium stress, the treatments AR11 + Fusarium and AR32 + Fusarium maintained significantly higher SPAD values compared to the Fusarium-only control, demonstrating their ability to alleviate pathogen-induced chlorophyll reduction. These statistical differences are supported by the Dunnett test results shown in Figure 3. In contrast, treatments such as VL13 + Fusarium and DR31 + Fusarium did not significantly differ from the Fusarium control, suggesting limited effectiveness under pathogen pressure. Overall, the data highlight the superior performance of AR11 and AR32 in supporting chlorophyll retention both in the presence and absence of Fusarium infection.

Figure 3
Boxplot showing SPAD measurements by strain, using Dunnett's test. Strains include Control, AR11, AR11+FU, AR32, AR32+FU, AR35, AR35+FU, Control+FU, CR25, CR25+FU, DR31, DR31+FU, VL13, VL13+FU, VR24, and VR24+FU. SPAD values range approximately from 35 to 55, with colored boxes indicating variations and points showing outliers. Significance compared to the control is noted.

Figure 3. Effect of endophytic bacterial inoculation on chlorophyll content (SPAD index) in pea plants under non-stress (green boxplots) and Fusarium-stressed (yellow boxplots) conditions. Boxplots display the distribution of SPAD values across treatments, with individual data points representing biological replicates (n = 4).

Correlation between SPAD measurements and root/shoot development

Correlation analysis between SPAD index and plant growth traits revealed notable variation across bacterial treatments, with several strains showing strong associations between chlorophyll content and overall plant performance (Figure 4). AR32 + Fusarium exhibited the strongest and most consistent positive correlations across all measured traits, including fresh root weight (r = 0.99), dried shoot weight (r = 0.94), root height (r = 0.97), and shoot height (r = 0.95), indicating that SPAD values reliably reflected plant vigor under stress. AR11 + Fusarium also demonstrated strong positive correlations, particularly with fresh root weight (r = 0.79) and root height (r = 0.79), suggesting that AR11 not only preserved chlorophyll levels but also supported below-ground biomass development during Fusarium challenge. In the non-stressed condition, AR11 showed high correlations between SPAD and dried shoot weight (r = 0.75) and shoot height (r = 0.74), confirming its broader role in enhancing shoot growth. Additionally, VR24 + Fusarium displayed strong correlations between SPAD and most traits, including shoot weight and root-related parameters. In contrast, treatments such as Control + Fusarium and CR25 showed weak or negative correlations, indicating a disrupted or inconsistent relationship between chlorophyll content and growth under stress. These findings reinforce the role of AR32 and AR11 in maintaining a strong physiological link between photosynthetic capacity and biomass production under both normal and Fusarium-infected conditions.

Figure 4
Heatmap displaying the correlation between SPAD and traits per strain. Rows represent different strains and treatments, and columns represent traits: DRW, DSW, FRW, FSW, RH, and SH. Colors range from red (strong positive correlation, 1.0) to blue (strong negative correlation, -1.0), with the Pearson correlation values indicated in each cell.

Figure 4. Correlation of SPAD index with plant growth parameters across bacterial treatments. The heatmap displays Pearson correlation coefficients between SPAD values (chlorophyll content) and plant growth traits dried root weight (DRW), dried shoot weight (DSW), fresh root weight (FRW), fresh shoot weight (FSW), root height (RH), and shoot height (SH), for each bacterial strain under both Fusarium-infected (+Fu) and non-infected conditions. Color intensity represents the strength and direction of correlation, with red indicating strong positive and blue indicating strong negative correlations.

Disease severity

The impact of different bacterial treatments on Fusarium disease severity varied notably. Treatments with AR11 + Fu, AR32 + Fu, and AR35 + Fu resulted in lower disease severity scores compared to the Control + Fu, indicating moderate to strong biocontrol activity (Figure 5). In contrast, the CR25 + Fu treatment showed disease severity levels like or higher than the Control + Fu, indicating it did not reduce disease symptoms. Likewise, DR31 + Fu, VL13 + Fu, and VR24 + Fu did not show significant effect on disease severity.

Figure 5
Violin plot illustrating disease severity by treatment, showcasing variations in severity scores from one to five across different treatments including AR11+Fu, AR32+Fu, AR35+Fu, Control+Fu, CR25+Fu, DR31+Fu, VL13+Fu, and VR24+Fu. Color-coded legend aids in distinguishing treatments.

Figure 5. Disease severity by treatment group: violin plot showing the distribution of disease severity scores (y-axis is ranging from 1 to 5) for each treatment group. The boxplots within the violins represent the median and interquartile range.

Genome features of strain AR11

The whole genome sequencing strain AR11 using the PacBio SMRT platform produced a complete circular genome of 5.49 Mb (Table 1). Genome assembly and annotation predicted 5,175 coding DNA sequences (CDSs), along with 84 tRNA genes, seven copies each of the 16S and 23S rRNA genes, and 25 small RNAs (sRNAs). Additional structural features included 15 genomic islands, six prophage regions, three CRISPR loci, and six biosynthetic gene clusters (BGCs) related to secondary metabolism. Epigenomic profiling identified 61,116 N6-methyladenine (m6A) and 6,484 N4-methylcytosine (m4C) modifications, along with 251,055 unclassified modified bases.

Table 1
www.frontiersin.org

Table 1. Summary of key genomic features identified in the AR11 genome Serratia sp.

Functional annotation of the genome revealed extensive database coverage. A total of 5,023 protein-coding genes (97.1%) received at least one functional annotation (Table 2). Annotation breakdown across major databases is summarized below. The majority of annotated genes were ≥300 bp in length, indicating well-supported predictions.

Table 2
www.frontiersin.org

Table 2. Functional annotation coverage of Serratia sp. AR11 coding genes across reference databases.

Additional structural and regulatory features were identified across the AR11 genome. Repeat sequence analysis showed that approximately 0.14% of the genome consisted of interspersed and tandem repeats. DNA methylation profiling detected 61,116 N6-methyladenine (m6A) and 6,484 N4-methylcytosine (m4C) modifications, along with 251,055 unclassified modified bases; notably, GATC motifs were nearly fully methylated. Mobile genetic elements included 15 genomic islands (totaling 265.8 kb), six prophage regions (409 kb), and three CRISPR loci. Six secondary metabolite biosynthetic gene clusters ranging from 25 to 49 kb were also identified. The circular genome map of Serratia sp. AR11 (Figure 6) illustrates the genomic architecture, including coding sequence distribution, GC content and skew, and the functional categories of genes based on COG classification. The complete genome sequence has been deposited in the NCBI database under Bio Sample accession number SAMN49938280 (Submission ID: SUB15457842).

Figure 6
Circular genome map with colored tracks representing various gene functions. Each color corresponds to a category like RNA processing (red), energy production (brown), and carbohydrate metabolism (light green), among others. A legend on the right lists categories A to Z with corresponding colors and counts, detailing various biological functions and processes.

Figure 6. Circular genome map of strain AR11. The map displays genome coordinates, coding sequences (CDSs) color-coded by COG functional categories, and key categories such as amino acid transport, transcription, and carbohydrate metabolism. The innermost graph illustrates GC content and GC skew. Genes are color-coded based on their COG functional categories, with a comprehensive key detailing each category and its associated color provided immediately alongside the map.

Taxonomic affiliation based on NR annotation

Strain AR11 was initially selected as Bacillus based on 16S rRNA gene analysis. However, whole-genome annotation against the NCBI non-redundant (NR) protein database revealed a different taxonomic affiliation. The majority of predicted protein-coding genes showed highest similarity to sequences from the genus Serratia. Specifically, 2,365 genes matched unclassified Serratia strains, and 2,041 matched Serratia proteamaculans. Additional matches were observed for S. plymuthica (164 genes), S. liquefaciens (107), S. marcescens (79), S. grimesii (28), and S. fonticola (6). In contrast, only a few gene matches were found with other genera, such as Yersinia frederiksenii (9 genes), Klebsiella pneumoniae (8), and Escherichia coli (8), further supporting the affiliation of AR11 with the genus Serratia (Figure 7). Due to multiple heterogeneous 16S rRNA gene copies in the genome, we selected the most abundant and representative sequence, based on coverage and similarity to Serratia, to construct the phylogenetic tree. This tree placed AR11 within the Serratia clade, closely related to S. proteamaculans (Figure 8). Average Nucleotide Identity (ANI) analysis using OrthoANIu (Yoon et al., 2017) further confirmed this placement, showing 95.69% identity between AR11 and S. proteamaculans, consistent with the species-level classification. In contrast, ANI with a representative Bacillus genome was only 63.35%, with genome coverage below 1%. These findings confirm classification of AR11 as Serratia proteamaculans, primarily supported by NR-based genome annotation, and demonstrate the limitations of relying solely on 16S rRNA for taxonomic assignment.

Figure 7
Bar chart titled

Figure 7. Taxonomic distribution of top NR matches. The bar plot depicts the number of predicted genes in the AR11 genome that matched protein sequences from the top 20 reference species in the NCBI NR database.

Figure 8
Phylogenetic tree showing relationships between various bacterial strains, including Serratia, Yersinia, Klebsiella, Enterobacter, and Bacillus. Branch points are marked with bootstrap values indicating the confidence levels of each branch.

Figure 8. Neighbor-joining phylogenetic tree illustrating the evolutionary relationships of bacterial strains based on 16S rRNA gene sequences. Bootstrap support values (1,000 replicates) are shown at the nodes. Evolutionary distances were computed using the maximum composite likelihood method. The analysis included 11 nucleotide sequences, with a final dataset of 1,567 positions (pairwise deletion). All analyses were conducted in MEGA11. The position of isolate AR11 is indicated.

Functional annotation and KEGG pathway classification

KEGG pathway annotation of strain AR11 revealed a wide range of functional capabilities, highlighting its metabolic versatility and ecological adaptability. Among the annotated genes, 142 were associated with amino acid metabolism, 118 with carbohydrate metabolism, and 88 with energy metabolism (Figure 9). These dominant categories indicate strong potential for essential metabolic activities required for growth and survival. The genome also contained genes directly linked to plant growth-promoting traits. Specifically, 28 genes were involved in nitrogen metabolism, 37 in the phosphotransferase system (PTS), and 41 in oxidative phosphorylation indicating pathways critical for nutrient acquisition and energy production.

Figure 9
Bar graph depicting the distribution of annotated genes across various biological processes, including genetic information processing, metabolism, environmental information processing, and cellular processes. Each category is color-coded: red for genetic processing, green for metabolism, purple for environmental processing, and blue for cellular processes. The percentage of annotated genes is displayed along the horizontal axis, with numeric values indicating the gene count. Metabolism has the highest concentration, particularly in biosynthesis of amino acids and carbon metabolism. Environmental information processing also shows a significant number with ABC transporters.

Figure 9. The KEGG pathway classification of AR11’s genome reveals high functional diversity, with major representation in amino acid and carbohydrate metabolism, energy production, ABC transporters, and two-component systems, emphasizing the strain’s metabolic adaptability and ecological competence.

Strain AR11 has an enhanced ability to detect and respond to environmental changes. It contains 239 genes encoding ATP-binding cassette (ABC) transporters and 138 genes associated with two-component regulatory systems, which play key roles in environmental signal detection and cellular response. Additionally, genes related to bacterial chemotaxis (23), flagellar assembly (36), and multiple secretion systems were detected. Several genes associated with secondary metabolite biosynthesis and membrane transport were also identified, pointing toward potential antimicrobial activity and adaptation to plant-associated niches.

Screening of the AR11 genome against the Comprehensive Antibiotic Resistance Database (CARD) identified several genes associated with intrinsic resistance, including efflux pumps (adeF, msbA, emrR) and regulators (CRP), as well as sequence variants of elongation factor Tu and penicillin-binding protein PBP3.

Functional annotation based on eggNOG classification

The genome of strain AR11 was annotated using the eggNOG database, which assigned its predicted proteins to 25 functional categories based on the Clusters of Orthologous Groups (COG) system. The most common category was R: General function prediction only, which included 944 proteins (19.37%), indicating a high number of proteins with unknown or putative functions.

A large number of proteins were related to amino acid transport and metabolism (443 proteins, 9.09%), transcription (428 proteins, 8.78%), and carbohydrate transport and metabolism (383 proteins, 7.86%). These results suggest that strain AR11 has flexible metabolic pathways and effective systems for regulating gene expression. Among the functional categories, P: Inorganic ion transport and metabolism (330 proteins, 6.77%), C: Energy production and conversion (273 proteins), and H: Coenzyme transport and metabolism (142 proteins) were prominently represented. Additionally, proteins related to V: Defense mechanisms (55 proteins), U: Intracellular trafficking and secretion (84 proteins), and N: Cell motility (75 proteins) were identified, indicating potential roles in environmental interactions and motility (Figure 10). As expected for a prokaryotic organism, no proteins were assigned to categories W (extracellular structures), Y (nuclear structure), or Z (cytoskeleton).

Figure 10
Bar chart titled

Figure 10. EggNOG functional classification of predicted proteins. The bar chart displays the distribution of predicted proteins across 25 COG categories. The most abundant category is R: general function prediction only, followed by categories associated with metabolism, transcription, and nutrient transport. The y-axis represents the number of proteins in each category. Protein classification highlights metabolic and regulatory diversity.

Gene Ontology enrichment analysis

Gene Ontology (GO) enrichment analysis of the AR11 genome categorized coding genes into three principal GO domains: biological process, molecular function, and cellular component (Figure 11). In the biological process category, approximately 700 genes were associated with metabolic processes, 650 with cellular processes, and 550 with single-organism processes. Additionally, around 400 genes were linked to transport, 300 to the regulation of biological processes, and 250 to the response to stimulus. Within the Molecular Function domain, catalytic activity was the most represented category (~600 genes), followed by protein binding (~500), nucleic acid binding (~400), and transporter activity (~350). In the cellular component category, genes were predominantly assigned to the cell (~450), membrane (~300), and membrane part (~200), indicating a notable prevalence of membrane-associated systems.

Figure 11
Bar graph displaying gene ontology categories with percentage of genes on the left Y-axis and number of genes on the right Y-axis. X-axis categories include cellular component (blue), molecular function (red), and biological process (green). Each category shows varying bar heights, indicating different gene counts and percentages across functions like extracellular region, transcription factor activity, and cellular process.

Figure 11. Gene Ontology (GO) enrichment analysis of the AR11 genome. Bar plots represent the number of genes associated with enriched GO terms across the three domains: biological process, molecular function, and cellular component.

Distribution of carbohydrate-active enzyme families

The genome of strain AR11 encodes a diverse array of carbohydrate-active enzymes (CAZymes), classified into six major CAZy families (Figure 12). The glycoside hydrolases (GHs) family was the most predominant, comprising 37.96% of all CAZyme-encoding genes. These enzymes play a key role in the hydrolysis of complex polysaccharides into simpler sugars. The second most abundant family was glycosyltransferases (GTs), accounting for 28.34% of CAZyme genes, followed by carbohydrate esterases (CEs) at 13.9%, and carbohydrate-binding modules (CBMs) at 10.16%. Less abundant families included auxiliary activities (AAs), which make up 9.09% and are involved in oxidative reactions, and polysaccharide lyases (PLs), representing 0.53%, which cleave acidic polysaccharides through non-hydrolytic mechanisms.

Figure 12
Pie chart displaying enzyme distribution: Glycoside Hydrolases (GH) 37.96%, Glycosyl Transferases (GT) 28.34%, Carbohydrate Esterases (CE) 13.9%, Carbohydrate-Binding Modules (CBM) 10.16%, Auxiliary Activities (AA) 9.09%, and Polysaccharide Lyases (PL) 0.53%.

Figure 12. Distribution of CAZyme families in the AR11 genome. The pie chart illustrates the relative proportions of each CAZyme family identified in the genome. GHs and GTs were the most abundant, reflecting AR11’s capacity for complex carbohydrate metabolism and potential involvement in plant-microbe interactions.

CRISPR-Cas system identification

To assess the genomic defense potential and adaptive immunity of Serratia strain AR11, an analysis of the CRISPR-Cas system was performed. Three distinct CRISPR loci were identified on Contig00001, each containing varying numbers of repeats and spacers (Table 3). The repeat sequences ranged from 22 to 30 base pairs, while the spacer sequences ranged from 24 to 31 base pairs. CRISPR.1 consisted of three repeats and two spacers, CRISPR.2 contained five repeats and four spacers, and CRISPR.3 included four repeats and three spacers. The presence of multiple CRISPR loci suggests the potential for a functional CRISPR-Cas system, which could play a role in the strain’s ability to defend against foreign genetic elements.

Table 3
www.frontiersin.org

Table 3. CRISPR loci identification in Serratia strain AR11 genome.

Discussion

This study highlights Serratia sp. AR11 as a promising plant growth-promoting endophyte (PGPE) with dual biostimulant and biocontrol activity. Isolated from Artemisia absinthium, AR11 significantly enhanced pea plant growth and mitigated the deleterious effects of Fusarium infection. Through integrated phenotypic assays and high-throughput genomic analysis, it is indicated the molecular mechanisms are likely to underly these beneficial traits. These results are consistent with numerous studies reporting that endophytic and rhizosphere-associated Serratia species contribute to plant growth promotion and stress mitigation through diverse functional mechanisms (Kshetri et al., 2019; Kulkova et al., 2024; Patel et al., 2024).

Genomic foundations of biostimulant activity

Under both normal and Fusarium-challenged conditions, Serratia sp. AR11 significantly enhanced root and shoot development in pea plants. This suggests that AR11 improves nutrient acquisition and enhances stress tolerance through its interaction with the root system. These beneficial effects are supported by the genomic enrichment of genes related to primary metabolism, energy generation, and nutrient transport. Functional annotations based on KEGG pathways and COG categories revealed a substantial number of genes involved in amino acid metabolism (443), carbohydrate metabolism (383), and inorganic ion transport and metabolism (330) traits commonly found in nutrient cycling rhizobacteria. Notably, the presence of 28 genes associated with nitrogen metabolism and 37 genes encoding components of the phosphotransferase system (PTS) indicates that AR11 can facilitate nitrogen assimilation and sugar uptake, both critical for plant growth and vigor. Similar metabolic profiles have been reported in well-characterized PGPRs such as Streptomyces, Bacillus velezensis, Bacillus, and Pseudomonas fluorescens, which are known to improve biomass accumulation and nutrient-use efficiency in host plants (Olanrewaju et al., 2021; Sanow et al., 2023; Xu et al., 2025). Moreover, the identification of 41 genes related to oxidative phosphorylation and 273 genes involved in energy production and conversion (COG category C) highlights AR11’s strong metabolic activity. This energy generating capacity likely supports its persistence and activity in both the rhizosphere and endospheric. Comparable energy-related gene profiles have been observed in PGPRs adapted to oligotrophic or nutrient-poor soils (Colin et al., 2017), supporting AR11’s suitability for use in low-input agricultural systems. Although genes specifically annotated phytohormone biosynthesis were not identified, the presence of key genes in the tryptophan biosynthesis pathway suggests that Serratia sp. AR11 may possess the capacity for auxin (IAA) production, a well-known plant growth-promoting trait in PGPRs. This genomic indication aligns with previous reports of IAA synthesis in Serratia spp., particularly S. plymuthica, has been stimulate root elongation and enhance plant growth (Rico-Jiménez et al., 2023). Future studies involving hormone quantification assays could help validate the functional expression of these genomic traits in planta. This mechanism may also contribute to the improved root development observed in AR11-treated plants.

Genomic drivers of biocontrol against Fusarium

The ability of AR11 to reduce stunting in Fusarium-infected pea plants appears to be supported by several genomic traits associated with biocontrol functions. Genome analysis revealed six distinct biosynthetic gene clusters (BGCs), each ranging from 25 to 49 kb in size, that are involved in the synthesis of secondary metabolites. Genome analysis identified six biosynthetic gene clusters potentially involved in secondary metabolite production, which may include compounds with antifungal properties contributing to Fusarium inhibition. While direct chemical profiling was beyond the scope of this study, these genomic features provide a valuable basis for future metabolomic validation to confirm the bioactive compounds and their roles in biocontrol. Similar antifungal compounds such as prodigiosin (Su et al., 2016), serratomolides (Marques-Pereira et al., 2020), and siderophores (Chlebek et al., 2022) have previously been reported in Serratia strains with demonstrated antifungal activity. Additionally, AR11’s genome contains a high number of glycoside hydrolases (GHs), as identified in the CAZyme analysis. These enzymes are known to break down fungal cell wall components, representing a direct antagonistic mechanism against fungal pathogens. Although direct chemical profiling was not conducted in this study, genome mining identified six biosynthetic gene clusters (BGCs) related to secondary metabolite production, including antifungal compounds. These in silico predictions, supported by CAZyme annotations, provide strong genomic evidence for the potential synthesis of bioactive molecules with antifungal activity. Such genomic insights offer a focused basis for future metabolomic validation. This finding is consistent with earlier studies on Streptomyces sp. GD-4, where GH-rich genomic regions contributed to antifungal activity and rhizosphere colonization (Xu et al., 2025). AR11 also appears to compete effectively with Fusarium and other soil microbes for space and nutrients. Its genome encodes 239 ABC transporters, 138 two-component regulatory systems, and numerous genes involved in amino acid and carbohydrate metabolism. ABC transporters help bacteria survive by protecting them from harmful compounds, helping them move toward plant roots, control gene activity, release useful substances, and take in nutrients (Zeng and Charkowski, 2021). These traits are essential for successful root colonization, exclusion of competing pathogens, and resource acquisition. Comparable features have been observed in Bacillus subtilis strains used for biocontrol in cucumber, where competitive interactions help limit pathogen establishment (Samaras et al., 2021). An additional key mechanism contributing to AR11’s biocontrol activity is induced systemic resistance (ISR), in which beneficial microbes prepare the plant’s immune system to respond more quickly and effectively to pathogen attack. AR11 contains 142 genes related to signal transduction (COG category T), along with 23 genes for chemotaxis and 36 for flagellar assembly, indicating its ability to sense plant signals and efficiently colonize root tissues. While induced systemic resistance (ISR) markers were not experimentally quantified, the genome of AR11 encodes key traits commonly associated with ISR-inducing microorganisms, including signal transduction systems, chemotaxis, and motility genes. These features indicate a strong potential to activate plant immune responses through ISR-related pathways. These traits are important because ISR is often triggered by microbial activity in the rhizosphere or within plant tissues. Previous studies have shown that ISR is a key defense mechanism induced by certain plant growth-promoting bacteria and fungi, which strengthens the plant’s resistance to a wide range of pathogens (Pieterse et al., 2014). Moreover, endophytic bacteria typically have more genes involved in motility and host communication such as those for chemotaxis, flagellar assembly, and pilus formation supporting their enhanced colonization and interaction with host plants (Jia et al., 2023). Similar ISR-inducing abilities have been reported in Bacillus amyloliquefaciens, which enhances plant defenses without directly attacking the pathogen (Salwan et al., 2023).

Stress tolerance, genome plasticity, and epigenomic regulation

AR11 harbors several adaptive features for environmental resilience. The abundance of two-component systems and ABC transporters enables real-time sensing and response to rhizosphere changes (Borland et al., 2015). These systems are critical for PGPRs colonizing fluctuating environments such as sandy, degraded, or low-nutrient soils. Functional enrichment analyses of bacterial microbiota have shown that ABC transporters tend to be more abundant in rhizosphere communities compared to bulk soil, while two-component systems are often less represented. This shift may reflect the importance of ABC transporters in nutrient uptake and microbial competition in the rhizosphere, where microbes are exposed to intense biological and chemical interactions (Ruan et al., 2019). The high representation of both systems in AR11 suggests a robust ability to adapt to the dynamic conditions of the plant-root interface and maintain functional activity across diverse soil types. The genome of AR11 contains 15 genomic islands and six prophage regions, indicating past horizontal gene transfer events that may have contributed to its expanded metabolic functions and biocontrol capabilities. These regions often carry mobile genetic elements (MGEs), which are key drivers of bacterial evolution and adaptability. MGEs facilitate the acquisition and exchange of beneficial traits, such as stress tolerance, antimicrobial compound synthesis, and nutrient uptake. Their mobility allows bacteria to respond rapidly to environmental challenges, contributing to genetic diversity and ecological fitness (Weisberg and Chang, 2023). In beneficial microbes like AR11, this genomic flexibility enhances their ability to colonize the rhizosphere, interact with plant hosts, and adapt to changing soil conditions, all of which are critical for sustained plant growth promotion and pathogen suppression.

The CRISPR-Cas system serves as an adaptive immune mechanism in prokaryotes, enabling them to recognize and neutralize invading genetic elements such as bacteriophages and plasmids. This defense is achieved by incorporating short sequences from these invaders into the host genome, which are then used to guide the degradation of matching foreign DNA during subsequent attacks (Scrascia et al., 2023). In the genome of AR11, three CRISPR loci, with variable spacer-repeat architecture, were also identified. These provide a mechanism for defense against foreign DNA and enhance genome stability traits commonly reported in Serratia strains include Type I-E and I-F1 CRISPR systems (Scrascia et al., 2023). Epigenomic analysis of AR11 revealed widespread DNA methylation, particularly at GATC motifs targeted by Dam methyltransferases. This form of adenine methylation is a conserved epigenetic mechanism in bacteria, playing a crucial role in regulating gene expression, maintaining genomic stability, and facilitating adaptation to environmental stresses (Qiao et al., 2024). In Serratia and related enterobacteria, Dam-mediated methylation influences various cellular processes, including biofilm formation and stress responses, by modulating the transcription of specific genes (Bruneaux et al., 2022). The presence of these methylation patterns in AR11 suggests an additional layer of regulatory control that enhances its ability to thrive in diverse and challenging environments. The observed effects under controlled conditions highlight AR11 as a strong candidate for practical application. Given the genomic traits related to environmental adaptation, colonization, and functional versatility, field-level and semi-field trials are planned to evaluate performance under natural conditions, including interactions with native soil microbiota and environmental variability.

Accurate taxonomic placement and safety considerations

The whole genome phylogenomic analysis reclassified strain AR11 as a member of the Serratia proteamaculans species, correcting its initial misidentification as Bacillus based on 16S rRNA sequencing. This reclassification is supported by Average Nucleotide Identity (ANI) values above the 95% species threshold and shared genomic characteristics with other Serratia strains known for plant growth-promoting and antifungal properties (Campos et al., 2024; Kulkova et al., 2024). To assess biosafety comprehensively, additional aspects such as horizontal gene transfer potential, environmental persistence, and possible non-target effects were considered. Biosafety evaluation of AR11’s genome revealed intrinsic antibiotic resistance genes related mainly to efflux pumps and regulatory proteins common in environmental bacteria. Importantly, no acquired antibiotic resistance or virulent genes associated with clinical pathogens were detected. Moreover, no conjugative elements, plasmids, or mobile resistance islands were identified, suggesting a low likelihood of horizontal gene transfer under field conditions. This supports AR11’s potential as a safe bioinoculant for sustainable agricultural use. While potential non-target effects on native symbionts such as rhizobia warrant future investigation, such interactions were not directly assessed in the current study. Consistent with previous recommendations, rigorous screening for pathogenicity and resistance genes is essential before field application of plant growth-promoting bacteria (Keswani et al., 2019). Given its genomic traits and observed performance, strain AR11 could be further developed for formulation as a bioinoculant or used as a model for selecting similar strains from other medicinal plants.

Conclusion

This study presents Serratia sp. AR11, isolated from Artemisia absinthium, as a distinct endophytic strain with combined plant growth-promoting and disease-suppressing effects in pea under Fusarium stress. Its strong performance in controlled greenhouse conditions, including improved growth, biomass, chlorophyll content, and reduced disease severity, highlights its practical potential. Genome analysis supports these observations, revealing genes related to nutrient metabolism, stress response, and the predicted synthesis of antifungal compounds. While metabolite production was not experimentally confirmed, the genomic data provides a strong foundation for further functional validation. Unlike many previously reported Serratia strains, AR11 combines a unique plant origin with a clear genomic and phenotypic profile, making it a valuable candidate for developing bioinoculants tailored to pea cropping systems. The absence of high-risk virulence or resistance genes further supports its safe application. Follow-up studies focusing on metabolomic assays and field trials will help confirm and extend these findings. Future studies integrating metabolomic, microbiome, and formulation approaches may further enhance the application potential of strain AR11.

Data availability statement

The data presented in this study are available in the NCBI BioSample repository under accession numbers SAMN49938280 and PRJNA1290800.

Author contributions

SH: Formal analysis, Validation, Methodology, Data curation, Supervision, Project administration, Software, Conceptualization, Writing – review & editing, Resources, Writing – original draft, Visualization, Investigation, Funding acquisition. DS: Writing – review & editing, Conceptualization, Supervision, Methodology, Project administration. SS: Supervision, Methodology, Conceptualization, Project administration, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

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.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Akhter, S., Aziz, R. K., and Edwards, R. A. (2012). PhiSpy: a novel algorithm for finding prophages in bacterial genomes that combines similarity- and composition-based strategies. Nucleic Acids Res. 40:e126. doi: 10.1093/nar/gks406

PubMed Abstract | Crossref Full Text | Google Scholar

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2

PubMed Abstract | Crossref Full Text | Google Scholar

Anand, U., Pal, T., Yadav, N., Singh, V. K., Tripathi, V., Choudhary, K. K., et al. (2023). Current scenario and future prospects of endophytic microbes: promising candidates for abiotic and biotic stress management for agricultural and environmental sustainability. Microb. Ecol. 86, 1455–1486. doi: 10.1007/s00248-023-02190-1

PubMed Abstract | Crossref Full Text | Google Scholar

Ashitha, A., Midhun, S. J., Sunil, M. A., Nithin, T. U., Radhakrishnan, E. K., and Mathew, J. (2019). Bacterial endophytes from Artemisia nilagirica (Clarke) Pamp., with antibacterial efficacy against human pathogens. Microb. Pathog. 135:103624. doi: 10.1016/j.micpath.2019.103624

PubMed Abstract | Crossref Full Text | Google Scholar

Aswani, R., Jishma, P., and Radhakrishnan, E. K. (2020). “Endophytic bacteria from the medicinal plants and their potential applications” in Microbial endophytes (Amsterdam: Elsevier), 15–36.

Google Scholar

Basit, A., Shah, S. T., Ullah, I., Ullah, I., and Mohamed, H. I. (2021). “Microbial bioactive compounds produced by endophytes (bacteria and fungi) and their uses in plant health” in Plant growth-promoting microbes for sustainable biotic and abiotic stress management (Cham: Springer), 285–318.

Google Scholar

Bertelli, C., and Brinkman, F. S. L. (2018). Improved genomic island predictions with IslandPath-DIMOB. Bioinformatics 34, 2161–2167. doi: 10.1093/bioinformatics/bty095

PubMed Abstract | Crossref Full Text | Google Scholar

Bland, C., Ramsey, T. L., Sabree, F., Lowe, M., Brown, K., Kyrpides, N. C., et al. (2007). CRISPR Recognition Tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinformatics 8:209. doi: 10.1186/1471-2105-8-209

PubMed Abstract | Crossref Full Text | Google Scholar

Blin, K., Shaw, S., Steinke, K., Villebro, R., Ziemert, N., Lee, S. Y., et al. (2019). antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 47, W81–W87. doi: 10.1093/nar/gkz310

PubMed Abstract | Crossref Full Text | Google Scholar

Borland, S., Oudart, A., Prigent-Combaret, C., Brochier-Armanet, C., and Wisniewski-Dyé, F. (2015). Genome-wide survey of two-component signal transduction systems in the plant growth-promoting bacterium Azospirillum. BMC Genomics 16:833. doi: 10.1186/s12864-015-1962-x

PubMed Abstract | Crossref Full Text | Google Scholar

Bruneaux, M., Kronholm, I., Ashrafi, R., and Ketola, T. (2022). Roles of adenine methylation and genetic mutations in adaptation to different temperatures in Serratia marcescens. Epigenetics 17, 861–881. doi: 10.1080/15592294.2021.1966215

PubMed Abstract | Crossref Full Text | Google Scholar

Campos, D., Cottet, L., Santos, C., and Castillo, A. (2024). Antifungal activity of Serratia plymuthica against the phytopathogenic fungus Alternaria tenuissima. Microb. Pathog. 193:106750. doi: 10.1016/j.micpath.2024.106750

PubMed Abstract | Crossref Full Text | Google Scholar

Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The carbohydrate-active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238. doi: 10.1093/nar/gkn663

PubMed Abstract | Crossref Full Text | Google Scholar

Chan, P. P., and Lowe, T. M. (2019). tRNAscan-SE: searching for tRNA genes in genomic sequences. New York, NY: Humana, 1–14.

Google Scholar

Chlebek, D., Grebtsova, V., Piński, A., Żur-Pińska, J., and Hupert-Kocurek, K. (2022). Genetic determinants of antagonistic interactions and the response of new endophytic strain Serratia quinivorans KP32 to fungal phytopathogens. Int. J. Mol. Sci. 23:15561. doi: 10.3390/ijms232415561

PubMed Abstract | Crossref Full Text | Google Scholar

Colin, Y., Nicolitch, O., Van Nostrand, J. D., Zhou, J. Z., Turpault, M.-P., and Uroz, S. (2017). Taxonomic and functional shifts in the beech rhizosphere microbiome across a natural soil toposequence. Sci. Rep. 7:9604. doi: 10.1038/s41598-017-07639-1

PubMed Abstract | Crossref Full Text | Google Scholar

Conesa, A., Götz, S., García-Gómez, J. M., Terol, J., Talón, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610

PubMed Abstract | Crossref Full Text | Google Scholar

Damavandi, M. S., Shojaei, H., and Esfahani, B. N. (2023). The anticancer and antibacterial potential of bioactive secondary metabolites derived from bacterial endophytes in association with Artemisia absinthium. Sci. Rep. 13:18473. doi: 10.1038/s41598-023-45910-w

PubMed Abstract | Crossref Full Text | Google Scholar

Dhillon, L. K., Lindsay, D., Yang, T., Zakeri, H., Tar’an, B., Knight, J. D., et al. (2022). Biological nitrogen fixation potential of pea lines derived from crosses with nodulation mutants. Field Crop Res. 289:108731. doi: 10.1016/j.fcr.2022.108731

Crossref Full Text | Google Scholar

El-Saadony, M. T., Saad, A. M., Najjar, A. A., Alzahrani, S. O., Alkhatib, F. M., Shafi, M. E., et al. (2021). The use of biological selenium nanoparticles to suppress Triticum aestivum L. crown and root rot diseases induced by Fusarium species and improve yield under drought and heat stress. Saudi J. Biol. Sci. 28, 4461–4471. doi: 10.1016/j.sjbs.2021.04.043

PubMed Abstract | Crossref Full Text | Google Scholar

Fontana, D. C., de Paula, S., Torres, A. G., de Souza, V. H. M., Pascholati, S. F., Schmidt, D., et al. (2021). Endophytic fungi: biological control and induced resistance to phytopathogens and abiotic stresses. Pathogens 10:570. doi: 10.3390/pathogens10050570

PubMed Abstract | Crossref Full Text | Google Scholar

Gupta, S., Pandey, S., and Sharma, S. (2022). Decoding the plant growth promotion and antagonistic potential of bacterial endophytes from Ocimum sanctum Linn. against root rot pathogen Fusarium oxysporum in Pisum sativum. Front. Plant Sci. 13:813686. doi: 10.3389/fpls.2022.813686

PubMed Abstract | Crossref Full Text | Google Scholar

Hadian, S., Smith, D. L., Kopriva, S., Norkevičienė, E., and Supronienė, S. (2024). Exploring endophytic bacteria from Artemisia spp. and beneficial traits on pea plants. Plants 13:1684. doi: 10.3390/plants13121684

PubMed Abstract | Crossref Full Text | Google Scholar

Hadian, S., Smith, D. L., and Supronienė, S. (2025). Modulating the plant microbiome: effects of seed inoculation with endophytic bacteria on microbial diversity and growth enhancement in pea plants. Microorganisms 13:570. doi: 10.3390/microorganisms13030570

PubMed Abstract | Crossref Full Text | Google Scholar

Hallmann, J., Quadt-Hallmann, A., Mahaffee, W. F., and Kloepper, J. W. (1997). Bacterial endophytes in agricultural crops. Can. J. Microbiol. 43, 895–914. doi: 10.1139/m97-131

Crossref Full Text | Google Scholar

Hyatt, D., Chen, G.-L., LoCascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi: 10.1186/1471-2105-11-119

PubMed Abstract | Crossref Full Text | Google Scholar

Jia, M., Chen, L., Xin, H.-L., Zheng, C.-J., Rahman, K., Han, T., et al. (2016). A friendly relationship between endophytic fungi and medicinal plants: a systematic review. Front. Microbiol. 7:906. doi: 10.3389/fmicb.2016.00906

PubMed Abstract | Crossref Full Text | Google Scholar

Jia, H., Yan, Y., Ma, J., Xia, E., Ma, R., Li, Y., et al. (2023). Complete genome sequence of a plant growth-promoting endophytic bacterium V4 isolated from tea (Camellia sinensis) leaf. Beverage Plant Res. 3:24. doi: 10.48130/BPR-2023-0024

Crossref Full Text | Google Scholar

Keswani, C., Prakash, O., Bharti, N., Vílchez, J. I., Sansinenea, E., Lally, R. D., et al. (2019). Re-addressing the biosafety issues of plant growth promoting rhizobacteria. Sci. Total Environ. 690, 841–852. doi: 10.1016/j.scitotenv.2019.07.046

PubMed Abstract | Crossref Full Text | Google Scholar

Krogh, A., Larsson, B., von Heijne, G., and Sonnhammer, E. L. L. (2001). Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J. Mol. Biol. 305, 567–580. doi: 10.1006/jmbi.2000.4315

Crossref Full Text | Google Scholar

Kshetri, L., Naseem, F., and Pandey, P. (2019). Role of Serratia sp. as biocontrol agent and plant growth stimulator, with prospects of biotic stress management in plant. Singapore: Springer, 169–200.

Google Scholar

Kulkova, I., Wróbel, B., and Dobrzyński, J. (2024). Serratia spp. as plant growth-promoting bacteria alleviating salinity, drought, and nutrient imbalance stresses. Front. Microbiol. 15:1342331. doi: 10.3389/fmicb.2024.1342331

PubMed Abstract | Crossref Full Text | Google Scholar

Lantzouraki, D. Z., Amerikanou, C., Karavoltsos, S., Kafourou, V., Sakellari, A., Tagkouli, D., et al. (2023). Artemisia arborescens and Artemisia inculta from Crete; secondary metabolites, trace metals and in vitro antioxidant activities. Life 13:1416. doi: 10.3390/life13061416

PubMed Abstract | Crossref Full Text | Google Scholar

Marques-Pereira, C., Proença, D. N., and Morais, P. V. (2020). Genome sequences of Serratia strains revealed common genes in both Serratomolides gene clusters. Biology 9:482. doi: 10.3390/biology9120482

PubMed Abstract | Crossref Full Text | Google Scholar

Nawrocki, E. P., and Eddy, S. R. (2013). Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29, 2933–2935. doi: 10.1093/bioinformatics/btt509

PubMed Abstract | Crossref Full Text | Google Scholar

Negri, S., Pietrolucci, F., Andreatta, S., Chinyere Njoku, R., Antunes Silva Nogueira Ramos, C., Crimi, M., et al. (2024). Bioprospecting of Artemisia genus: from artemisinin to other potentially bioactive compounds. Sci. Rep. 14:4791. doi: 10.1038/s41598-024-55128-z

PubMed Abstract | Crossref Full Text | Google Scholar

Nongkhlaw, F. M., and Joshi, S. R. (2015). Investigation on the bioactivity of culturable endophytic and epiphytic bacteria associated with ethnomedicinal plants. J. Infect. Dev. Ctries. 9, 954–961. doi: 10.3855/jidc.4967

PubMed Abstract | Crossref Full Text | Google Scholar

Olanrewaju, O. S., Ayilara, M. S., Ayangbenro, A. S., and Babalola, O. O. (2021). Genome mining of three plant growth-promoting Bacillus species from maize rhizosphere. Appl. Biochem. Biotechnol. 193, 3949–3969. doi: 10.1007/s12010-021-03660-3

PubMed Abstract | Crossref Full Text | Google Scholar

Patel, S. K., Singh, S., Benjamin, J. C., Singh, V. R., Bisht, D., and Lal, R. K. (2024). Plant growth-promoting activities of Serratia marcescens and Pseudomonas fluorescens on Capsicum annuum L. plants. Ecol. Front. 44, 654–663. doi: 10.1016/j.ecofro.2024.01.002

Crossref Full Text | Google Scholar

Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat. Methods 8, 785–786. doi: 10.1038/nmeth.1701

PubMed Abstract | Crossref Full Text | Google Scholar

Pieterse, C. M. J., Zamioudis, C., Berendsen, R. L., Weller, D. M., Van Wees, S. C. M., and Bakker, P. A. H. M. (2014). Induced systemic resistance by beneficial microbes. Annu. Rev. Phytopathol. 52, 347–375. doi: 10.1146/annurev-phyto-082712-102340

PubMed Abstract | Crossref Full Text | Google Scholar

Qiao, S., Song, W., Hu, W., Wang, F., Liao, A., Tan, W., et al. (2024). The role of plant DNA methylation in development, stress response, and crop breeding. Agronomy 15:94. doi: 10.3390/agronomy15010094

Crossref Full Text | Google Scholar

Raza, A., Hassan, A., Akram, W., Anjum, T., Aftab, Z.-H., and Ali, B. (2024). Seed coating with the synthetic consortium of beneficial Bacillus microbes improves seedling growth and manages Fusarium wilt disease. Sci. Hortic. 325:112645. doi: 10.1016/j.scienta.2023.112645

Crossref Full Text | Google Scholar

Rico-Jiménez, M., Muñoz-Mira, S., Lomas-Martínez, C., Krell, T., and Matilla, M. A. (2023). Regulation of indole-3-acetic acid biosynthesis and consequences of auxin production deficiency in Serratia plymuthica. Microb. Biotechnol. 16, 1671–1689. doi: 10.1111/1751-7915.14296

PubMed Abstract | Crossref Full Text | Google Scholar

Ruan, R., Jiang, Z., Wu, Y., Xu, M., and Ni, J. (2019). High-throughput sequence analysis reveals variation in the relative abundance of components of the bacterial and fungal microbiota in the rhizosphere of Ginkgo biloba. Peer J. 7:e8051. doi: 10.7717/peerj.8051

PubMed Abstract | Crossref Full Text | Google Scholar

Saier, M. H. (2006). TCDB: the transporter classification database for membrane transport protein analyses and information. Nucleic Acids Res. 34, D181–D186. doi: 10.1093/nar/gkj001

Crossref Full Text | Google Scholar

Salwan, R., Sharma, M., Sharma, A., and Sharma, V. (2023). Insights into plant beneficial microorganism-triggered induced systemic resistance. Plant Stress 7:100140. doi: 10.1016/j.stress.2023.100140

Crossref Full Text | Google Scholar

Samaras, A., Nikolaidis, M., Antequera-Gómez, M. L., Cámara-Almirón, J., Romero, D., Moschakis, T., et al. (2021). Whole genome sequencing and root colonization studies reveal novel insights in the biocontrol potential and growth promotion by Bacillus subtilis MBI 600 on cucumber. Front. Microbiol. 11:600393. doi: 10.3389/fmicb.2020.600393

PubMed Abstract | Crossref Full Text | Google Scholar

Sanow, S., Kuang, W., Schaaf, G., Huesgen, P., Schurr, U., Roessner, U., et al. (2023). Molecular mechanisms of Pseudomonas-assisted plant nitrogen uptake: opportunities for modern agriculture. Mol. Plant Microbe Interact. 36, 536–548. doi: 10.1094/MPMI-10-22-0223-CR

PubMed Abstract | Crossref Full Text | Google Scholar

Scrascia, M., Roberto, R., D’Addabbo, P., Ahmed, Y., Porcelli, F., Oliva, M., et al. (2023). Bioinformatic survey of CRISPR loci across 15 Serratia species. MicrobiologyOpen 12:e1339. doi: 10.1002/mbo3.1339

Crossref Full Text | Google Scholar

Su, C., Xiang, Z., Liu, Y., Zhao, X., Sun, Y., Li, Z., et al. (2016). Analysis of the genomic sequences and metabolites of Serratia surfactantfaciens sp. nov. YD25T that simultaneously produces prodigiosin and serrawettin W2. BMC Genomics 17:865. doi: 10.1186/s12864-016-3171-7

PubMed Abstract | Crossref Full Text | Google Scholar

Sudheer, S., Bai, R. G., Usmani, Z., and Sharma, M. (2020). Insights on engineered microbes in sustainable agriculture: biotechnological developments and future prospects. Curr. Genomics 21, 321–333. doi: 10.2174/1389202921999200603165934

PubMed Abstract | Crossref Full Text | Google Scholar

Tarailo-Graovac, M., and Chen, N. (2009). Using RepeatMasker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinformatics 25, 4.10.1–4.10.14. doi: 10.1002/0471250953.bi0410s25

PubMed Abstract | Crossref Full Text | Google Scholar

Timofeeva, A. M., Galyamova, M. R., and Sedykh, S. E. (2023). Plant growth-promoting soil bacteria: nitrogen fixation, phosphate solubilization, siderophore production, and other biological activities. Plants 12:4074. doi: 10.3390/plants12244074

PubMed Abstract | Crossref Full Text | Google Scholar

Tripathi, A., Pandey, P., Tripathi, S. N., and Kalra, A. (2025). Plant growth-promoting endophytic consortium improved artemisinin biosynthesis via modulating antioxidants, gene expression, and transcriptional profile in Artemisia annua (L.) under stressed environments. Plant Stress 15:100757. doi: 10.1016/j.stress.2025.100757

Crossref Full Text | Google Scholar

Weisberg, A. J., and Chang, J. H. (2023). Mobile genetic element flexibility as an underlying principle to bacterial evolution. Ann. Rev. Microbiol. 77, 603–624. doi: 10.1146/annurev-micro-032521-022006

Crossref Full Text | Google Scholar

Xu, W., Liu, Y., Cheng, Y., and Zhang, J. (2025). Plant growth-promoting effect and complete genomic sequence analysis of the beneficial rhizosphere Streptomyces sp. GD-4 isolated from Leymus secalinus. Microorganisms 13:286. doi: 10.3390/microorganisms13020286

PubMed Abstract | Crossref Full Text | Google Scholar

Yoon, S.-H., Ha, S., Lim, J., Kwon, S., and Chun, J. (2017). A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie Van Leeuwenhoek 110, 1281–1286. doi: 10.1007/s10482-017-0844-4

PubMed Abstract | Crossref Full Text | Google Scholar

Zeng, Y., and Charkowski, A. O. (2021). The role of ATP-binding cassette transporters in bacterial phytopathogenesis. Phytopathology 111, 600–610. doi: 10.1094/PHYTO-06-20-0212-RVW

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Artemisia, biocontrol, biostimulants, endophytic bacteria, Serratia

Citation: Hadian S, Smith DL and Suproniene S (2025) Genomic and phenotypic profiling of an Artemisia endophyte: dual biostimulant and biocontrol activities in pea under Fusarium stress. Front. Microbiol. 16:1643204. doi: 10.3389/fmicb.2025.1643204

Received: 08 June 2025; Accepted: 07 August 2025;
Published: 29 August 2025.

Edited by:

Md. Motaher Hossain, Bangabandhu Sheikh Mujibur Rahman Agricultural University, Bangladesh

Reviewed by:

Abhinav Aeron, Chonbuk National University, Republic of Korea
Nagaraju Yalavarthi, Central Silk Board, India
Zhaogao Li, Zunyi Medical University, China
Chunli Zeng, Nanjing Agricultural University, China
Sukhminderjit Kaur, Chandigarh University, India

Copyright © 2025 Hadian, Smith and Suproniene. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Shervin Hadian, c2hlcnZpbi5oYWRpYW5AbGFtbWMubHQ=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.