ORIGINAL RESEARCH article

Front. Plant Sci., 03 September 2025

Sec. Plant Biotechnology

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1657785

Ecological and genetic diversity of plant growth-promoting genes in rhizobacteria isolated from the rhizosphere of wild flora on Mount Erciyes, Türkiye

  • 1. Department of Crop and Animal Production, Safiye Cikrikcioglu Vocational College, Kayseri University, Kayseri, Türkiye

  • 2. Genome and Stem Cell Research Center, Erciyes University, Kayseri, Türkiye

Abstract

Introduction:

This study is the first to investigate the genetic diversity of plant growth-promoting genes in rhizobacteria isolated from the wild ecology of Mount Erciyes, Türkiye. It has a flora rich in flowering plants, with 1170 plant taxa, 194 of which are endemic to the area.

Methods:

A total of 165 bacterial isolates, including Azotobacter sp., Azospirillum sp., and Bacillus sp., were screened for genes associated with plant growth promotion: nitrogen fixation (nif), indole pyruvate decarboxylase (ipdC), 1-aminocyclopropane-1-carboxylate deaminase (accd), phosphate-solubilizing genes (Acpho, Alpho and phy), and siderophore biosynthesis (sd).

Results:

The analysis revealed significant genetic variability across isolates, particularly for nif and sd genes, with distinct band patterns indicating genetic diversity among Azospirillum and Bacillus isolates.

Discussion:

The findings emphasize the role of these rhizobacteria in nutrient cycling and stress resilience, potentially enhancing plant growth in nutrient-limited soils. In the current study, it contributes to understanding microbial biodiversity in Mount Erciyes and suggests a promising potential for sustainable agriculture through plant-microbe interactions.

1 Introduction

Recent advancements in agricultural research have increasingly focused on sustainable practices to enhance crop productivity while minimizing environmental impacts. In semi-arid regions like the southeastern part of the Central Anatolian Plateau, including Kayseri, Türkiye, cereal and legume cultivation faces challenges due to nutrient-poor aridisols and limited water availability. Plant growth-promoting rhizobacteria (PGPR) have emerged as a promising solution to these challenges, offering eco-friendly alternatives to chemical fertilizers. Studies have demonstrated that PGPR strains, such as Azospirillum and Bacillus, can significantly enhance crop yields. For instance, a study by reported that Azospirillum brasilense and Bacillus subtilis increased wheat grain yields by up to 30% in semi-arid conditions through mechanisms like nitrogen fixation and phosphate solubilization. Similarly, found that co-inoculation of Bacillus and Rhizobium improved chickpea nodulation and yield by 25% in nutrient-deficient soils, highlighting the potential of PGPR for sustainable agriculture in regions like Kayseri.

The application of PGPR extends beyond yield enhancement to include stress tolerance and soil health improvement. Research by emphasized that PGPR, including Azospirillum and Bacillus, promote legume growth by producing phytohormones like indole-3-acetic acid (IAA) and solubilizing phosphorus, reducing reliance on synthetic fertilizers. Additionally, PGPR contribute to biocontrol by producing siderophores and antimicrobial compounds, as noted in a meta-analysis by , which reported an 11.4% increase in soybean nodulation with Bradyrhizobium and Bacillus co-inoculation. These findings underscore the multifaceted roles of PGPR in enhancing nutrient cycling and plant resilience, critical for sustainable farming in semi-arid ecosystems.

Türkiye’s diverse ecological landscape, particularly around Mount Erciyes in Kayseri, supports a rich biodiversity that remains underexplored for microbial resources. The southeastern part of the Central Anatolian Plateau, including the Kayseri region, is dominated by cereal and legume cultivation on semi-arid, nutrient-poor aridisols. Plant growth-promoting rhizobacteria (PGPR) strains such as Azospirillum and Bacillus have been shown to increase wheat and barley grain yields by up to 40%, providing a sustainable, low-input alternative to chemical fertilizers that could greatly benefit local farmers (). Türkiye’s diverse ecological structure fosters high levels of endemism and genetic diversity. Mount Erciyes, a prominent stratovolcano situated approximately 25 km southwest of Kayseri, rises from the plains of Sultansazlığı and hosts a flora rich in flowering plant species. Of the 1170 plant taxa identified in this region, 194 are endemic (). Although many studies have explored the biodiversity of plant and animal life on Mount Erciyes, a significant knowledge gap remains regarding the rhizobacteria inhabiting this unique environment. To date, there has been no research on PGPR or their associated functional genes in this area. PGPR are beneficial bacteria commonly isolated from soil and are known to promote plant growth through various mechanisms. These bacteria belong to diverse genera, including Alcaligenes, Agrobacterium, Azospirillum, Azotobacter, Arthrobacter, Bacillus, Bradyrhizobium, Burkholderia, Caulobacter, Chromobacterium, Enterobacter, Erwinia, Flavobacterium, Herbaspirillum, Klebsiella, Mesorhizobium, Micrococcus, Pseudomonas, Rhizobium, Rhodococcus, and Serratia (). Among these genera, Azospirillum and Azotobacter are noteworthy for their nitrogen-fixing abilities, while Bacillus is notable for producing plant hormones, siderophores, and phosphate-solubilizing enzymes.

In this study, we focused on characterizing the genetic diversity of plant growth-promoting genes within rhizobacterial populations, specifically targeting Azospirillum sp., Azotobacter sp., and Bacillus sp., isolated from the unique flora of Mount Erciyes. To assess their plant growth-promoting potential, we screened for genes associated with key PGPR functions. These included nitrogen fixation genes (nif), indole pyruvate decarboxylase (ipdC) involved in auxin (IAA) biosynthesis, 1-aminocyclopropane-1-carboxylate deaminase (accd) associated with ethylene modulation, and genes involved in phosphate solubilization, such as acid phosphatase (Acpho), alkaline phosphatase (Alpho), and phytase (phy). Additionally, we screened for siderophore biosynthesis genes (sd) essential for iron acquisition and plant growth enhancement.

2 Materials and methods

2.1 Collection of soil samples

Soil samples were collected from various soil types (clayey, sandy, and loamy soils) within the vegetative regions of Mount Erciyes, with each sample consisting of a 1 kg mixture obtained by combining 100 g of soil from 10 different spots within an area. The sampling sites were strategically distributed to capture a range of ecological niches, including areas near endemic plant species and agricultural fields, to ensure a comprehensive representation of the rhizobacterial populations associated with the region’s unique flora. Samples were taken from the plant root-soil interface area (5-10 cm depth) in locations free from commercial microbial fertilizers. In total, 22 soil samples were collected, with GPS coordinates recorded for each sampling location (Table 1).

Table 1

Location nameLatitudeLongitudeAltitude (m)
1. KıranardıN 38 36’ 51,90640”E 35 32’ 10,47160”1616
2. HisarcıkN 38 37’ 32,01250”E 35 30’ 25,428301542
3. EndürlükN 38 38’ 9,70420”E 35 33’ 4,56050”1376
4. KepezN 38 36’ 59,83560”E 35 39’ 55,80450”1697
5. YazyurduN 38 34’ 48,22060”E 35 40’ 40,70250”1641
6. CebirN 38 34’ 5,21330”E 35 37’ 0,670201942
7. KoçcağızN 38 30’ 17,01790”E 35 42’ 14,21860”1603
8. ÇaylıcaN 38 24’ 57,44160”E 35 37’ 11,02440”1475
9. SoysallıN 38 23’ 24,73720”E 35 21’ 43,81860”1074
10. ÇayırözüN 38 24’ 58,48490”E 35 17’ 20,58580”1073
11. KulpakN 38 27’ 27,01850”E 35 19’ 40,23290”1287
12. Ahmet pınarıN 38 27’ 49,09050”E 35 24’ 7,49670”1863
13. TekirN 38 30’ 45,68570”E 35 31’ 15,93220”2209
14. HacılarN 38 35’ 26,66210”E 35 27’ 54,21670”2074
15. Sakar ÇiftliğiN 38 37’ 29,16490”E 35 23’ 52,13680”1461
16. KızılörenN 38 36’ 6,01990”E 35 20’ 21,60570”1542
17. SarıkürklüN 38 39’ 5,77110”E 35 18’ 8,17520”1133
18. SubaşıN 38 33’ 6,09570”E 35 14’ 13,52140”1141
19. KızıliniN 38 32’ 7,63710”E 35 17’ 28,36790”1623
20. ŞehşabanN 38 28’ 52,61930”E 35 14’ 48,15300”1156
21. KarpuzsekisiN 38 41’ 4,99640”E 35 18’ 50,90630”1031
22. HürmetçiN 38 42’ 0,09030”E 35 18’ 46,35280”1031

Soil samples and their GPS locations.

2.2 Physicochemical characterization of soil samples

A total of 22 different soil samples were analyzed for texture, pH (potential of hydrogen), lime (calcium carbonate equivalent, CaCO3), organic matter (OM), electrical conductivity (EC), phosphorus (P), total nitrogen (N), iron (Fe), copper (Cu), zinc (Zn), and manganese (Mn) properties. Soil texture was determined using the Bouyoucos hydrometer method (). Soil pH was measured potentiometrically with a glass electrode pH meter (Cole-Parmer Jenway 3510 Standard Digital pH Meter) using a 1:2.5 soil-to-water ratio (). Lime content was determined volumetrically with a Scheibler calcimeter. Organic matter content was measured using the Walkley–Black wet oxidation method (). Exchangeable cations (Na+, Ca²+, Mg²+, and K+) were extracted by shaking the soils with 1 N ammonium acetate (pH 7.0) and quantified by inductively coupled plasma optical emission spectrometry (ICP-OES) (). Available phosphorus was determined using the molybdenum blue colorimetric method by measuring the absorbance of the resulting blue solution at 660 nm with a spectrophotometer (). Total nitrogen content was calculated using the micro-Kjeldahl method after wet digestion with a mixture of salicylic acid, sulfuric acid, and salt. Micronutrients, including Fe, Cu, Zn, and Mn, were extracted using the DTPA (diethylenetriaminepentaacetic acid) method and quantified by ICP-OES ().

2.3 Bacteria isolation and culture conditions

2.3.1 Isolation and identification of Azotobacter sp.

To isolate Azotobacter sp. from soil samples, 10 grams of each soil sample were mixed with 90 ml of dH2O. After thorough vortex mixing, the samples were left to settle at room temperature for approximately 20 min. Following this, 1 ml of the supernatant was transferred to 50 ml flasks containing Ashby liquid medium (pH 7.3) (per liter: 20 g mannitol, 0.2 g K2HPO4, 0.2 g MgSO4•7H2O, 0.2 g NaCl, 0.1 g K2SO4, 5 g CaCO3, and 1 ml of a microelement solution containing: 1 g H3BO3, 2 g FeSO47H2O, 1 g MnCl24H2O, 1 g Na2MoO4•2H2O, 0.5 g NaBr, and 0.2 g ZnSO4•7H2O) and incubated at 28°C for 5-7 days at 200 rpm in a shaking incubator (). From the resulting culture, 1 ml was taken and serially diluted (10-4 to 10-7) before being plated on Ashby agar (Azotobacter Agar with Mannitol). Plates were incubated at 28°C for 5-7 days. Azotobacter colonies grown on the selective medium were stored at +4°C until DNA isolation ().

2.3.2 Isolation and identification of Azospirillum sp.

To isolate Azospirillum sp., 10 grams of each soil sample were mixed with 90 ml of sterile distilled water. After thorough vortex mixing, the samples were left to settle at room temperature for approximately 20 minutes. Then, 1 ml of the liquid portion was serially diluted (10-² to 10-7) (), and added to 10 ml of N-free medium (liquid NFb medium, 1 L, pH 6.5: 5 g malic acid, 0.05 g yeast extract, 0.5 g K2HPO4, 0.2 g MgSO47H2O, 0.1 g NaCl, 0.02 g CaCl22H2O, 0.04 g CuSO45H2O, 0.12 g ZnSO47H2O, 1.4 g H3BO3, 1 g Na2MoO42H2O, 1.175 g MnSO4H2O, 2 ml bromothymol blue [5 g/L in 0.2 N KOH], 2 ml Fe-EDTA solution [16.4 g/L], 4 ml vitamin solution [10 mg biotin; 20 mg pyridoxal-HCl, pH adjusted with KOH; for solid NFb medium, add 15 g agar]) (). Azospirillum sp. colonies grown on the selective medium were stored at +4°C until DNA isolation.

2.3.3 Isolation and identification of Bacillus sp.

To isolate Bacillus sp., 10 grams of each soil sample were mixed with 90 ml of sterile distilled water. After thorough vortexing, the samples were allowed to settle at room temperature for approximately 5 minutes. Then, 1 ml of the supernatant was transferred to 50 ml flasks containing 10 ml of Luria-Bertani (LB) Broth medium buffered with 0.25 M sodium acetate. The flasks were incubated in a shaking incubator at 30°C for 4 hours at 200 rpm. Following incubation, 1 ml of each sample was placed in sterile eppendorf tubes and heat-treated at 80°C for 5-10 minutes. Afterward, samples were plated on LB agar plates and incubated overnight at 30°C (; ). Bacillus bacterial colonies grown on the selective medium were stored at +4°C until DNA isolation.

2.4 DNA isolation and amplification of plant-growth-promoting genes

DNA isolation was performed using the Bio Basic Bacterial Genomic DNA Isolation Kit (BS624), following the manufacturer’s recommended protocol. From each isolated DNA sample, 5 μl were loaded onto a 1% (w/v) agarose gel prepared in 1x TAE buffer (Tris-Acetic Acid-EDTA, pH 8.0) containing 10 μl ethidium bromide, and electrophoresed at 80 V for 1 hour. DNA presence was verified using a Biorad gel imaging system.

To determine whether the rhizobacteria identified as Azotobacter, Azospirillum, and Bacillus carried genes promoting plant growth, PCR analysis was conducted using primers (nifH-1, nifH-2, nifH-3, nifH-4, nifH-5, nifH-6, nifH-7, nifH-8, nifH-9, Acpho-1, Acpho-2, Acpho-3, Acpho-4, Alpho-1, Alpho-2, Alpho-3, phy-1, phy-2, sd-1, sd-2, sd-3, accd, ipdC-1, ipdC-2, ipdC-3, ipdC-4, ipdC-5, ipdC-6, ipdC-7, ipdC-8 and ipdC-9) for the nif, ipdC, accd, Acpho, Alpho, phy, and sd genes (Table 2) (; ; ; ; ; ; ; ; ; ; ; ).

Table 2

PrimersPrimer sequence (5´→3´)Annealing (°C)Target gene
nifH-1FCAGACACGAAGAAGCCGGGC50Nitrogenase
RGACCAGCAGCTTGTTGTTGA
nifH-2FCGCCGGCGCAGTGTTTGCGG50
RCACTCGTTGCAGCTGTCGGC
nifH-3FGGTTGTGACCCGAAAGCTGA50
RGCGTACATGGCCATCATCTC
nifH-4FGGCTGCGATCCCGAAAGGCCGACTTCCGAACCCG55
RCTGGCAGCCTTGTTCTTCGCGGATCGGGCATGGC
nifH-5FGGCTGCGATCCAAGGCCGATCACCCG50
RCTGGCCTTGTTTCGCGGATGGCATGGC
nifH-6FGGCAAGGGCGGTATCGGCAAGTC55
RCCATCGTGATCGGGTCGGGATG
nifH-7FGACCCGCCTGATCCTGCACG55
RGTTCTCTTCCAGGAAGTTGATCGA
nifH-8FGACCCGCCTGATCCTGCACG55
RGTCGTAGGCGCCGTTC
nifH-9FGACCCGCCTGATCCTGCACG55
RCATGACGATGTAGATTTCCTGG
Acpho-1FAAGAGGGGCATTACCACTTTATTA55Acid phosphatase
RCGCCTTCCCAATCRCCATACAT
Acpho-2FCGGCTCCTATCCGTCCGG58
RCAACATCGCTTTGCCAGTG
Acpho-3FGTCGGCTTTGATATCGATGA55
RCGAGCGACCTCTTTTGGAAT
Acpho-4FTCTTCGGTGATAACTTAGGTGACTT55
RAATCGTATAAAGCGCCTTCC
Alpho-1FCAGTGGGACGACCACGAGGT57Alkaline phosphatase
RGAGGCCGATCGGCATGTCG
Alpho-2FGCCGGAGCTCATGAGTTTATTTAAACAGGTTCGAT50
RGCCGCTGCAGTTATTTGCCGCTTTTTAAGATG
Alpho-3FGGAATTCCATATGGGTTTCTTACGCAACAGAAT55
RATAGTTTAGCGGCCGCTCTGGCGTATTTTTTGAATAGCT
phy-1FTATGATTTTCCGTTGAAC48Phytase
RATTCCGTCTGTATCGCTTGT
phy-2FCTGTCTGATCCTTATCATTT52
RTCCGCTTCTGTCGGTCA
sd-1FGAGAATGGATTACAGAGGAT55Siderophore biosynthesis
RTTATGAACGAACAGCCACTT
sd-2FGGAGAATGGATTACAGAGG48
RGTCGTCATATAAATTTCCAG
sd-3FACGATTGCACAATATGAAAGA55
RCAATGGTTTGGAACTTCATG
accdFGTGAACCACCTGAATGTA50acc deaminase
RAAACGAGATGATTTACTTGG
ipdC-1FCACTTGAAAACGCAATATACTG50Indole pyruvate decarboxylase
RAAGAATTTGCTGGCCGAATCT
ipdC-2FCATTTGAAAACTCACTATACTG51
RAAGAATTTGTATGCCGAATCT
ipdC-3FCACTTGAAAACTCACTATACTG53
RAAGAATTTGCATGCCGAATCT
ipdC-4FCATTTGAAAACGCAATATACTG52
RAAGAATTTGTTGGCCGAATCT
ipdC-5FCACTTGAAAACTCAATATACTG50
RAAGAATTTGCTTGCCGAATCT
ipdC-6FCACTTGAAAACGCACTATACTG55
RAAGAATTTGCTGGCCGAATCT
ipdC-7FCACTTGAAAACGCAATATACTG54
RAAGAATTTGTTGGCCGAATCT
ipdC-8FCACTTGAAAACGCAATATACTG54
RAAGAATTTGTAGGCCGAATCT
ipdC-9FCACTTGAAAACGCAATATACTG50
RAAGAATTTGTATGCCGAATCT

Primer pairs used to identify growth-promoting genes by PCR analysis.

Each reaction contained the reagents at a final concentration as 2.4 mM MgCl2, 1× taq buffer, 0.25 mM dNTPmix, 0.3 pmol primers (each), 0.5 U taq DNA polymerase, and 30–100 ng of template DNA. The PCR amplification was performed under the following conditions: Initial denaturation at 95°C for 5 min, followed by 30 cycles at 94°C for 1 min, Tm C for 1 min, 72°C for 2 min, and a final extension step at 72°C for 10 min (). Following PCR, 5 μl of each product was loaded onto a 1% (w/v) agarose gel prepared in 1x TAE buffer (Tris-Acetic Acid-EDTA, pH 8.0) containing 10 μl ethidium bromide, and electrophoresed at 80 V for 1 hour. DNA bands were visualized on a Biorad gel imaging system. GeneRuler 100 bp Plus DNA Ladder (Thermo Scientific), ranging from 100 to 3000 bp, was used as the molecular size marker in all gels to estimate the fragment sizes.

2.5 Statistical analysis

Statistical analyses were performed using the MSTAT-C statistical software (Michigan State University, USA). A one-way analysis of variance (ANOVA) was applied to assess differences in soil physicochemical properties among the 22 sampling sites on Mount Erciyes. Mean comparisons were conducted using the Tukey-Kramer Honestly Significant Difference (HSD) post hoc test at a 99% confidence level (p < 0.01). All parameters presented in Table 3, including pH, electrical conductivity (EC), organic matter (OM), nitrogen (N), lime content, phosphorus (P), calcium (Ca), potassium (K), magnesium (Mg), sodium (Na), copper (Cu), iron (Fe), manganese (Mn), and zinc (Zn), were subjected to this analysis. However, due to the absence of replicate measurements per location (n=1), within-group variance could not be estimated, rendering traditional significance testing inconclusive. Consequently, the data were interpreted descriptively, focusing on ranges and patterns across locations.

Table 3

LocationsTexturepHEC µs/cmOM %N %Lime %P kg/daCa mg/kgK mg/kgMg mg/kgNa mg/kgCu mg/kgFe mg/kgMn mg/kgZn mg/kg
1. KıranardıClay loam7.36206.02.910.151.531.17739.6135.069.561.80.520.450.720.65
2. HisarcıkClay loam7.76365.02.450.121.232.23788.5143.043.074.50.620.620.750.78
3. EndürlükClay loam7.42272.02.630.132.653.451256.3132.537.564.50.540.690.810.81
4. KepezClay loam7.38212.02.510.133.441.561365.2122.556.078.80.510.570.860.95
5. YazyurduClay loam7.30220.02.020.103.212.581452.2131.557.5101.20.450.690.910.91
6. CebirClay loam7.16174.02.450.125.423.451574.5139.043.0123.30.630.750.560.88
7. KoçcağızLoam6.99222.03.120.163.5510.251956.0120.539.086.50.580.780.851.02
8. ÇaylıcaLoam6.97486.02.750.146.858.541845.0128.029.595.60.510.590.751.11
9. SoysallıClay loam8.101709.01.850.097.419.651923.0127.532.51523.60.541.020.461.06
10. ÇayırözüClay loam8.271708.05.020.2518.4217.222027.0180.5107.01648.10.430.560.841.15
11. KulpakClay loam8.13622.04.880.2415.2515.241856.3162.087.5755.20.420.680.750.78
12.Ahmet pınarıClay loam7.64234.05.120.2612.3510.262063.0163.084.0302.50.960.950.560.89
13. TekirLoam6.54168.44.560.2311.0012.472024.0167.562.5123.60.580.840.911.25
14. HacılarClay loam7.14186.54.750.2412.368.652156.0156.068.0154.80.650.710.881.14
15.Sakar ÇiftliğiClay loam7.40236.03.860.198.459.361863.0144.072.5112.30.430.800.760.78
16. KızılörenClay loam7.22152.13.950.209.1211.252036.0122.542.5163.90.570.750.940.94
17. SarıkürklüClay loam7.51145.14.120.219.768.541765.0143.056.0102.50.610.620.580.91
18. SubaşıClay loam7.55337.04.440.225.456.321742.0156.034.5185.20.250.690.910.66
19. KızıliniClay loam7.33116.04.450.224.665.481556.3182.539.0223.60.360.880.780.78
20. ŞehşabanClay loam7.63295.03.800.193.533.58968.5132.031.5148.80.410.570.521.02
21. KarpuzsekisiClay loam7.78424.03.950.205.456.531623.5144.043.0123.60.320.690.631.15
22. HürmetçiClay loam8.084760.04.120.216.124.771552.3146.038.54565.20.250.720.591.03

Main physico-chemical characteristics of the soil from Mount Erciyes.

No significant differences were observed among the locations for the measured physicochemical parameters, as replicate measurements were unavailable for statistical testing (p-values undefined; one-way ANOVA and Tukey-Kramer HSD test at p < 0.01). Data are presented descriptively to highlight variations in ranges.

3 Results

3.1 Physicochemical characterization of soil samples

Table 3 summarizes the main physicochemical properties of the soil samples collected from the sampling sites. Soil samples collected from 22 distinct locations within the vegetative regions of Mount Erciyes were subjected to comprehensive physicochemical characterization to assess their properties and suitability for supporting rhizobacterial populations. These locations, spanning diverse microenvironments across the slopes of Mount Erciyes, were selected to capture variations in soil texture, nutrient content, and ecological niches. Table 3 summarizes the primary physicochemical properties of the soil samples, including texture, pH, electrical conductivity (EC), organic matter (OM) content, total nitrogen (N), lime content, and concentrations of key nutrients such as phosphorus (P), calcium (Ca), potassium (K), magnesium (Mg), sodium (Na), copper (Cu), iron (Fe), manganese (Mn), and zinc (Zn).

Soil texture analysis revealed that the majority of the samples (19 out of 22) were classified as clay loam, while three locations (Koçcağız, Çaylıca, and Tekir) exhibited a loamy texture. This predominance of clay loam soils suggests a high capacity for nutrient and water retention, which is critical for supporting microbial activity and plant growth in the semi-arid conditions of the Central Anatolian Plateau. The pH values ranged from slightly acidic to moderately alkaline (6.54 to 8.27), with most samples falling within a neutral to slightly alkaline range (7.14–8.13). Notably, the Tekir location exhibited the most acidic pH (6.54), while Çayırözü had the highest pH (8.27), potentially influencing the solubility and availability of nutrients and the activity of rhizobacteria.

Electrical conductivity (EC), an indicator of soil salinity, varied significantly across the sampling sites, ranging from 116.0 µS/cm (Kızılini) to 4760.0 µS/cm (Hürmetçi). The exceptionally high EC at Hürmetçi suggests potential salinity stress, which could impact microbial communities and plant growth. Organic matter (OM) content ranged from 1.85% (Soysallı) to 5.12% (Ahmet Pınarı), indicating moderate to high organic content that supports microbial activity. Total nitrogen (N) content followed a similar trend, ranging from 0.09% (Soysallı) to 0.26% (Ahmet Pınarı), reflecting the variability in soil fertility across the sites.

Lime content, which affects soil structure and nutrient availability, was notably high at Çayırözü (18.42%) and Kulpak (15.25%), potentially influencing pH and microbial processes.

Nutrient analysis revealed variation in phosphorus (P), calcium (Ca), potassium (K), magnesium (Mg), and sodium (Na) concentrations. Phosphorus levels ranged from 1.17 kg/da (Kıranardı) to 17.22 kg/da (Çayırözü), indicating a wide range of phosphate availability that could influence the activity of phosphate-solubilizing rhizobacteria. Calcium concentrations were consistently high, ranging from 739.6 mg/kg (Kıranardı) to 2156.0 mg/kg (Hacılar), reflecting the calcareous nature of the region’s soils. Potassium and magnesium levels showed moderate variation, with K ranging from 120.5 mg/kg (Koçcağız) to 182.5 mg/kg (Kızılini) and Mg from 29.5 mg/kg (Çaylıca) to 107.0 mg/kg (Çayırözü). Sodium levels were notably elevated at Hürmetçi (4565.2 mg/kg) and Soysallı (1523.6 mg/kg), further indicating potential salinity challenges at these sites. Trace elements, including copper (Cu), iron (Fe), manganese (Mn), and zinc (Zn), were present in low but variable concentrations, with Zn showing the highest values at Tekir (1.25 mg/kg) and Çayırözü (1.15 mg/kg), which may support specific microbial functions such as siderophore production. Due to the lack of replicate samples per location, one-way ANOVA and Tukey-Kramer HSD tests (p < 0.01) could not detect significant differences among locations for any parameter, as within-group variance was inestimable.

3.2 Bacterial diversity

A total of 165 microbial isolates were collected from 22 distinct soil samples taken from the wild rhizosphere region of Mount Erciyes. Based on colony morphology and Gram staining, the isolates were identified as follows: 50 Azotobacter sp. (30.3%), 42 Azospirillum sp. (25.5%), and 73 Bacillus sp. (44.2%). Genomic DNA was extracted from all bacterial isolates, and PCR analysis was performed to assess the presence and genetic diversity of several plant growth-promoting genes, including nif, ipdC, accd, Acpho, Alpho, phy, and sd. This approach enabled the identification of functional gene variants associated with plant growth promotion across the microbial community. Soil samples collected from Mount Erciyes displayed varying physicochemical characteristics based on location, which likely influenced microbial diversity and genetic functions observed. For example, the Kıranardı and Hisarcık soils, both with clay loam textures and slightly alkaline pH (7.3–7.7), had moderate organic matter (2.45–2.91%) and nitrogen content (0.12–0.15%). These conditions are suitable for nitrogen-fixing bacteria like Azotobacter and Azospirillum, which were prevalent in these regions due to the availability of nutrients and optimal pH levels.

Samples from Koçcağız and Çaylıca, featuring loamy soils with higher organic matter (3.12% and 2.75%, respectively), provided a nutrient-rich environment, which likely supports a diverse microbial community including Bacillus species. The variability in soil texture, organic matter, and pH across these sampling sites underscores how soil type influences the presence and diversity of plant growth-promoting genes, including those involved in nitrogen fixation and phosphate solubilization. Such diversity in microbial genes highlights the adaptive mechanisms and potential for these soils to support sustainable agriculture in response to nutrient demands.

3.3 Molecular screening of plant growth-promoting genes

3.3.1 Nitrogen fixation genes

Screening of 50 Azotobacter sp. isolates for nifH genes using various primers revealed diverse amplification patterns. The nifH-1 primer was expected to produce a 1102 bp band (Figure 1); however, 24 isolates exhibited bands of varying haplotypes with sizes deviating from this expectation (Table 4).

Figure 1

Table 4

PrimersAzotobacter sp.Azospirillum sp.Bacillus sp.Band size (bp)
nifH-1242331102
nifH-21816246
nifH-3102346
nifH-484323
nifH-5422 (different size)323
nifH-66360
nifH-7Variable haplotypes3205
nifH-8Variable haplotypes2205
nifH-9Variable haplotypesVariable haplotypes
ipdC-621Variable
ipdC-3 to 818 (500 bp larger than 200)200
accd391129800
acpho-11184734
alpho/phy
sd-2/31171685

Summary of bacterial species producing bands for plant growth-promoting genes.

Screening with the nifH-4 primer resulted in eight isolates producing the expected 323 bp band, while the nifH-5 primer yielded the anticipated 323 bp band in four isolates. The nifH-6 primer generated a 360 bp band in six isolates, consistent with expectations. In contrast, the nifH-7 primer produced gene patterns with a range of haplotypes significantly different from the expected 205 bp size. These findings indicate considerable genetic diversity in nifH genes among Azotobacter sp. isolates from different regions (Supplementary Figure 1).

Screening of 42 Azospirillum sp. isolates with nifH primers also showed varied results. The nifH-1 primer was expected to amplify a 1102 bp band, but 23 isolates displayed bands of different haplotypes with sizes deviating from this standard. The nifH-2 primer produced the expected 246 bp band in 18 isolates, while the nifH-3 primer yielded a 346 bp band in 10 isolates. For the nifH-4 primer, four isolates generated the expected 323 bp band, and the nifH-5 primer produced the anticipated 323 bp band in two isolates. The nifH-7 primer resulted in three isolates showing the expected 205 bp band. Additionally, screening with the nifH-8 primer revealed gene patterns with diverse haplotypes (Supplementary Figures 2A, B).

For 73 Bacillus sp. isolates, screening with nifH primers identified 37 isolates producing bands with various primers. The nifH-1 primer was expected to yield a 1102 bp band, and three isolates met this expectation. The nifH-2 primer produced the anticipated 246 bp band in 16 isolates, while the nifH-3 primer resulted in two isolates showing the expected 346 bp band. Screening with the nifH-5 primer revealed three isolates with bands deviating from the expected 323 bp size, though two isolates produced the anticipated 323 bp band. The nifH-8 primer identified two isolates with the expected 205 bp band. Finally, screening with the nifH-9 primer indicated gene patterns with diverse haplotypes (Supplementary Figure 3).

3.3.2 Indole pyruvate decarboxylase genes

Screening of 50 Azotobacter sp. isolates with nine ipdC primers yielded no amplification products, indicating the absence of detectable ipdC genes in these isolates. In contrast, screening of 42 Azospirillum sp. isolates revealed that 21 isolates produced bands of varying sizes with the ipdC-6 primer (Supplementary Figure 4), suggesting the presence of diverse ipdC gene haplotypes. Among 73 Bacillus sp. isolates, 18 were found to carry the ipdC gene. However, all positive isolates displayed band sizes approximately 500 bp larger than the expected 200 bp, indicating significant genetic variation from the anticipated profile (Table 4, Supplementary Figure 5).

3.3.3 1-aminocyclopropane-1-carboxylate deaminase genes

Screening of 50 Azotobacter sp. isolates for the accd gene revealed that 39 isolates produced bands of approximately 800 bp, consistent with the expected size (Supplementary Figure 6). Similarly, among 42 Azospirillum sp. isolates, 11 displayed bands of the anticipated size for the accd gene (Supplementary Figure 7). For Bacillus sp., 29 out of 73 isolates were found to carry the accd gene, all producing bands of the expected size (Table 4, Supplementary Figure 8).

3.3.4 Phosphate-solubilizing genes (acpho, alpho and phy)

Screening with the Acpho-1 primer identified 11 Azotobacter sp., 8 Azospirillum sp., and 4 Bacillus sp. isolates that produced a band of the expected 734 bp size, indicating the presence of the acid phosphatase gene (Supplementary Figures 9-11). In contrast, PCR amplification using Alpho and phy primers yielded no bands for any of the tested bacterial isolates. However, some Bacillus sp. isolates produced non-specific bands (Table 4), which do not confirm the presence of these genes due to their lack of specificity.

3.3.5 Siderophore biosynthesis genes

PCR screening for siderophore biosynthesis (sd) genes using sd-2 and sd-3 primers revealed that none of the 50 Azotobacter sp. isolates carried these genes. In contrast, 11 Azospirillum sp. isolates and 7 Bacillus sp. isolates produced the expected 1685 bp band, confirming the presence of sd genes (Figure 2). Additionally, some of these isolates exhibited extra bands of significantly larger sizes (Table 4, Supplementary Figures 12, 13), indicating potential genetic diversity in siderophore biosynthesis pathways.

Figure 2

4 Discussion

Nitrogen (N2), a critical component of nucleotides, proteins, and chlorophyll in plants, is one of the most essential elements for life. Despite constituting approximately 78% of the atmosphere, nitrogen is not readily available in a form usable by plants (). Soil is prone to nitrogen loss through processes such as denitrification, volatilization, and leaching due to nitrogen’s high reactivity and mobility. The leached reactive forms can cause environmental issues and adverse health effects (). Nitrogen fixation, the process of converting atmospheric N2 into ammonia (NH3), varies widely among bacteria and is catalyzed by the nitrogenase enzyme complex, which consists of two main components: an iron (Fe) protein encoded by the nifH gene and a molybdenum-iron (MoFe) protein encoded by the nifDK genes. The nifH gene is evolutionarily conserved and serves as a marker for nitrogen fixation, with primers designed from its sequence widely used to detect the genetic potential for nitrogen fixation in bacteria. While nif genes are sometimes located on chromosomes, they are often found on plasmids alongside other nitrogen-fixing genes, such as nod genes (; ; ). In addition to encoding nitrogenase, nif genes regulate other enzymes involved in nitrogen fixation. Their expression is triggered by low nitrogen and oxygen levels, particularly around the root zone. Although the nif genes encoding nitrogenase are highly conserved across nitrogen-fixing bacteria, their regulation varies among diazotrophs depending on their evolutionary lineage. Approximately 20 nif genes encode the nitrogenase complex.

Transcription of nif genes is initiated under nitrogen stress, often activated by the nitrogen-sensitive nifA protein. When fixed nitrogen is scarce, nifC, an RNA polymerase, induces nifA expression, which in turn activates other transcription factors. In environments with sufficient reduced nitrogen or oxygen, the nifL protein suppresses nifA, halting nitrogenase production. Studies have identified nitrogenase activity in various Bacillus species isolated from rhizosphere soils, including B. megaterium, B. cereus, B. pumilus, B. circulans, B. licheniformis, B. subtilis, B. brevis, and B. firmus (; ; ; ; ). For instance, isolated 14 Bacillus strains capable of acetylene reduction from rice fields along the Yangtze River in China. reported high nitrogenase activity in B. fusiformis from Chungbuk Province, South Korea, while highlighted the nitrogenase activity of B. alkalidiazotrophicus, a low salt-tolerant alkaliphile from Mongolian soda soils, demonstrating the genetic diversity of nitrogenase activity in Bacillus species across different ecosystems. This nitrogen fixation capacity directly contributes to plant growth promotion by increasing the availability of nitrogen, a limiting nutrient for plant development. For example, nitrogen-fixing Bacillus species have been shown to enhance the growth of crops such as rice and wheat by improving nitrogen uptake, leading to increased biomass and grain yield ().

Phytohormones are vital for signaling and regulating plant growth and development. Auxins, particularly indole-3-acetic acid (IAA), are among the most studied plant growth regulators (; ). Bacteria associated with plants produce phytohormones that influence plant growth (), plant pathology (), and microbial interactions (). Tryptophan serves as a precursor for IAA biosynthesis (). Over 80% of rhizosphere bacteria can produce IAA, which enhances root branching, weight, size, and surface area, improving nutrient exchange in plants (). Bacillus species possess the ipdC gene, which encodes indole pyruvate decarboxylase, a key enzyme in IAA production via the indole-3-pyruvate (IPyA) pathway (). The production of IAA by Bacillus species promotes root system development, which enhances nutrient and water uptake, leading to improved plant vigor and yield. For instance, Bacillus subtilis inoculation has been shown to increase root biomass and grain yield in maize by up to 15% under field conditions, demonstrating its potential for agricultural applications ().

Azospirillum species, widely distributed in tropical, subtropical, and temperate soils, produce phytohormones such as IAA, gibberellins, cytokinins, and abscisic acid, promoting plant growth (; ). Of the 18 identified Azospirillum species (), few are known to enhance plant growth, with A. brasilense being the most commonly used to improve crop yields under field conditions. A. brasilense synthesizes IAA from tryptophan through three pathways: indole pyruvic acid (IPyA), tryptamine, and indole acetonitrile (; , ). The IPyA pathway, involving transamination of tryptophan to IPyA by aromatic amino acid aminotransferases, followed by decarboxylation to indole acetaldehyde and oxidation to IAA, is the best-characterized. The ipdC gene, encoding phenylpyruvate decarboxylase, is critical in this pathway, and its presence has been confirmed in many Azospirillum strains (). The application of A. brasilense has been associated with significant increases in crop productivity, including a 20–30% yield improvement in wheat and maize under field conditions, primarily due to enhanced root development and nutrient assimilation ().

1-Aminocyclopropane-1-carboxylic acid (ACC), a natural ethylene precursor, regulates seed germination, senescence, fruit ripening, wound healing, and plant development (). ACC deaminase, found in plant growth-promoting bacteria, converts ACC into α-ketobutyrate and ammonium, promoting root elongation during seed germination (). Bacillus species with the accd gene enhance root-shoot development and increase dry and wet weights in plants (; ). By lowering ethylene levels, ACC deaminase-producing bacteria mitigate stress responses in plants, leading to improved growth under adverse conditions such as drought or salinity. For example, Bacillus amyloliquefaciens has been shown to increase tomato plant biomass by 25% under salt stress, highlighting its role in enhancing plant resilience and productivity ().

Although phosphorus (P) is abundant in soil, its insoluble form limits plant uptake, making it a critical constraint for growth. Soil phosphorus exists as mineral phosphates (e.g., calcium phosphates, hydroxyapatite, and rock phosphate) and organic phosphates (e.g., phosphoesters, phytates, or inositol phosphates), but plants can only absorb it as soluble monobasic or dibasic ions (). Phosphate-solubilizing bacteria produce enzymes that convert insoluble phosphates into soluble forms (). The primary mechanism involves organic acid production, which acidifies the microbial environment, displacing Ca²+ ions to release ionic phosphate. Bacillus species are key phosphate-solubilizing bacteria, utilizing enzymes like acid phosphatases (Acpho), alkaline phosphatases (Alpho), and phytases (phy) for dephosphorylation (; ; ; ). Phosphatase genes in B. cereus and B. thuringiensis show 99% homology with those of B. cereus ATCC 14579T (). Alkaline phosphatases, encoded by phoA, phoD, and phoX gene families, and phytases, with varying catalytic mechanisms, exhibit high microbial diversity (; ; ; ). In marine systems, phoD and phoX predominate over phoA, while phoD is the most abundant alkaline phosphatase gene in terrestrial soils, correlating with potential phosphatase activity (; ; ; ). Phosphate solubilization by Bacillus species enhances phosphorus availability, leading to improved plant growth and yield. For instance, B. megaterium has been reported to increase phosphorus uptake in soybean, resulting in a 10–15% increase in seed yield under greenhouse conditions ().

Iron (Fe) deficiency, an abiotic stress, can reduce crop yields (). Plants release chelators and phytosiderophores to bind Fe³+, enhancing its solubility and reducing it to Fe²+ for uptake. Rhizospheric bacteria, including Bacillus species like B. cereus, B. anthracis, and B. thuringiensis, produce siderophores via non-ribosomal peptide synthesis (; ; ; ). These bacteria carry siderophore biosynthesis genes (sd) and produce siderophores like bacillibactin (a catecholate-type siderophore) in Bacillus and azotobactin in A. vinelandii (). Plants in non-sterile soils show no iron deficiency symptoms, indicating the role of microbial siderophores in enhancing iron uptake (). Siderophores also combat plant pathogens by competing for iron, thus supporting plant growth (; ; ). The application of siderophore-producing Bacillus species has been shown to improve iron nutrition in crops like maize, leading to a 10% increase in chlorophyll content and a 12% increase in biomass under iron-deficient conditions ().

In summary, the multifaceted contributions of nitrogen-fixing, phytohormone-producing, ACC deaminase-active, phosphate-solubilizing, and siderophore-producing bacteria, particularly Bacillus and Azospirillum species, significantly enhance plant growth and productivity. These mechanisms collectively improve nutrient availability, stimulate root development, mitigate stress, and enhance resistance to pathogens, leading to increased crop yields and sustainable agricultural practices. For example, field studies have demonstrated that inoculation with Bacillus and Azospirillum strains can increase crop yields by 10–30% across various crops, including maize, wheat, and soybean, under diverse environmental conditions (; ). These findings underscore the potential of these bacteria as biofertilizers to enhance plant production, offering a sustainable alternative to chemical fertilizers while addressing global food security challenges.

5 Conclusion

Plant growth-promoting rhizobacteria (PGPR) exhibit significant genetic diversity in key genes, including those for nitrogen fixation (nif), indole pyruvate decarboxylase (ipdC), 1-aminocyclopropane-1-carboxylate deaminase (accd), phosphate solubilization (Acpho, Alpho, phy), and siderophore biosynthesis (sd). These genes enable PGPR to fix atmospheric nitrogen, synthesize phytohormones, and enhance the availability of phosphate and iron, thereby promoting plant growth and resilience under diverse environmental conditions. This study investigates the genetic diversity of PGPR traits in bacterial isolates from Mount Erciyes, Türkiye, revealing the influence of local soil properties:

Hisarcık: Nutrient-rich, neutral to slightly alkaline soils harbor diverse nif genes, particularly in Azotobacter species, supporting robust nitrogen fixation in fertile environments.

Kıranardı: Soils with moderate pH and balanced micronutrients foster diversity in accd and phosphate-solubilizing genes, enhancing plant stress tolerance and nutrient uptake.

Kepez: Iron-deficient soils at higher altitudes show increased diversity in sd genes among Bacillus isolates, improving iron acquisition for both plants and microbes.

Endürlük: Acidic soils with variable mineral content promote diversity in phosphate-solubilizing genes, reflecting microbial adaptation to phosphorus scarcity.

The variability in PGPR traits closely corresponds to soil composition and altitude, highlighting the adaptability of rhizobacteria to specific ecological niches. These insights emphasize the potential of tailored microbial consortia to support sustainable agriculture.

Statements

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. Further inquiries can be directed to the corresponding author.

Author contributions

UA: Writing – original draft, Methodology, Formal Analysis, Funding acquisition, Supervision, Conceptualization, Writing – review & editing, Project administration, Data curation, Investigation, Validation.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This study was supported by Scientific and Technological Research Council of Türkiye (TÜBİTAK) under the Grant Number 123O376 and Kayseri University Scientific Project Unit the Grant Number FOA-2023-1115. The authors thank to TÜBİTAK and Kayseri University Scientific Project Unit for their supports.

Acknowledgments

The author thanks Dr. Adem Güneş for his valuable assistance about soil analysis and Dr. Salih Karabörklü for statistical data analysis.

Conflict of interest

The author declares 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 author(s) declare that no Generative 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.

Supplementary material

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

References

Summary

Keywords

plant hormones, gene analyses, biodiversity, rhizobacteria, PGPR

Citation

Azizoglu U (2025) Ecological and genetic diversity of plant growth-promoting genes in rhizobacteria isolated from the rhizosphere of wild flora on Mount Erciyes, Türkiye. Front. Plant Sci. 16:1657785. doi: 10.3389/fpls.2025.1657785

Received

01 July 2025

Accepted

18 August 2025

Published

03 September 2025

Volume

16 - 2025

Edited by

Debasis Mitra, Graphic Era University, India

Reviewed by

Chitra Bhattacharya, Atmiya University, India

Olivia Landau, USDA-ARS Wheat Health Genetics and Quality Research, United States

Updates

Copyright

*Correspondence: Ugur Azizoglu, ; ;

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics