Skip to main content

Influence of genotype, nodule position, and edaphic factors on microbial diversity and assembly of pigeonpea (Cajanus cajan) root nodules in Indian soils

Abstract

Background

Pigeonpea (Cajanus cajan) is an important legume crop in semi-arid regions with multiple uses. The microbial diversity within its root nodules in Indian soils remains poorly explored. We investigated the bacterial diversity of pigeonpea root nodules across different genotypes and soil types to identify the factors driving their assembly. Using a metagenomic approach and high-throughput sequencing of the 16S rRNA gene, we analyzed the nodule microbiomes of three pigeonpea genotypes (Asha, Durga, and Mannem Konda Kandi) grown in three different soil types (Alfisol, Vertisol, and Inceptisol) and wild pigeonpea (C. scarabaeoides) in its native soil.

Results

Our results indicated that pigeonpea nodules harbor diverse rhizobial and non-rhizobial endophytes and that host genotype, nodule position, soil type, and other edaphic factors influence significant variation in the microbial community structure. The core nodule microbiome was dominated by Proteobacteria and Bacteroidetes. Bradyrhizobium and Ensifer were predominant among the rhizobial taxa, and non-rhizobial genera such as Pseudomonas, Chitinophaga, and Limnobacter were also abundant. Edaphic factors, particularly soil type, pH, and nutrient availability, had a stronger influence on the nodule bacterial community composition than the host genotype. Although bulk soil exhibited higher bacterial diversity, nodule microbiomes were less diverse but more specialized, indicating host-mediated selection. A comparison of the nodule microbiomes of wild and cultivated pigeonpea revealed distinct differences, with the core nodule microbiome of wild pigeonpea dominated by Bradyrhizobium, while that of cultivated pigeonpea exhibited a diverse bacterial community.

Conclusions

These findings demonstrate that soil properties play a more critical role than host genetics in shaping the pigeonpea nodule microbiome, emphasizing the importance of environmental conditions in symbiotic interactions. The differences between wild and cultivated genotypes suggest that domestication has altered microbial recruitment strategies. This study provides foundational insights into the factors driving microbial assembly in pigeonpea nodules, with implications for improving crop productivity through targeted microbial management. Future research should explore the functional roles of these microbial communities to optimize their use in sustainable agriculture.

Background

Pigeonpea [Cajanus cajan (L.) Millspaugh], a vital legume food crop, has diverse applications as food, feed, fodder, and fuel. Like other legumes, it also enriches soil through biological nitrogen fixation. Globally, pigeonpea cultivation spans approximately 6.0 million hectares [1], primarily as a rain-fed crop in the semi-arid tropical and subtropical regions of South Asia, East Africa, Latin America, and the Caribbean [2, 3]. Over a billion people in developing countries rely on pigeonpea as their primary dietary protein source. Millions of resource-poor smallholder farmers cultivate this multipurpose crop to sustain their livelihoods, often with limited inputs. Domestication of its wild progenitor species, Cajanus cajanifolius (endemic to the Indian subcontinent), gave rise to cultivated pigeonpea in central India more than 3,500 years ago, subsequently leading to its global expansion [2, 3].

India is the world’s leading producer of pigeonpea, accounting for 72% of the global supply [1]. Among cultivated legumes in India, it ranks second after chickpea, contributing 15% of the area under cultivation and 17% of the total pulse production [4]. The major pigeonpea growing zones in India can be divided into three distinct regions: the southern zone (comprising Andhra Pradesh, Karnataka, and Telangana), the central zone (encompassing Madhya Pradesh, Gujarat, and Maharashtra), and the northern plain zone (primarily Uttar Pradesh) [5]. Notably, pigeonpea yields in Andhra Pradesh, Madhya Pradesh, and Uttar Pradesh consistently surpass the national average [4]. The predominant soil types in these states are red soil (Alfisol) in Andhra Pradesh, black soil (Vertisol) in Madhya Pradesh, and alluvial soil (Inceptisol) in Uttar Pradesh [6].

Rhizobia were once considered to be the sole inhabitants of legume nodules. Several rhizobia can colonize nodules of the same plant and can even co-occupy the same nodule [7, 8]. Recent culture-dependent and culture-independent approaches have confirmed the existence of a diverse nodule microbiome, suggesting that rhizobia coexist with non-rhizobial nodule endophytes [9, 10], challenging the widely accepted notion that individual nodules typically harbor only one rhizobial strain. Instead, nodules may harbor multiple species, including non-rhizobial endophytes [11], also known as ‘non-rhizobial root nodule endophytes’ (NREs) [12].

Pigeonpea can establish symbiotic associations with diverse rhizobial genera [13]. Culture-dependent studies suggest that pigeonpea can be nodulated by Bradyrhizobium spp. [14], Ensifer (formerly Sinorhizobium) spp. [15, 16], Rhizobium spp. [17, 18], Mesorhizobium spp. [15, 17], or even Burkholderia (reclassified as Paraburkholderia) spp. [17], indicating that pigeonpea may be promiscuous in recruiting rhizobia for nodulation in Indian soils. However, in certain geographical regions such as the Dominican Republic [19, 20], Côte d’Ivoire [21], Brazil [22], and Ethiopia [23], pigeonpea plants form nodules exclusively with Bradyrhizobium spp. or Ensifer spp. [13, 24], suggesting that additional species reported in nodulation studies may require re-evaluation. In addition to rhizobial symbionts, pigeonpea roots harbor a diverse community of non-rhizobial colonizers, including bacteria from the genera Agrobacterium, Azotobacter, Azospirillum, Bacillus, Brevibacillus, Chryseobacterium, Enterobacter, Klebsiella, Lactobacillus, Paenibacillus, Pseudomonas, Serratia, Streptomyces, and others [21, 25, 26].

Whole-genome sequencing has enabled the characterization of the genetic and genomic diversity of pigeonpea, which has a genome size of 833.07 Mb [3, 27, 28], while the nodule microbiome remains largely uncharacterized. The nodule microbiomes of legumes such as common bean [29, 30], Medicago truncatula [31], clover [32], soybean [33], lentil [34], groundnut [35], Sophora davidii [36], Alnus spp. [37, 38], Prosopis cineraria [39], Desmodium spp. [40], chickpea [41], and sea buckthorn [42] have been reported in different countries but poorly explored in India. A comprehensive study of the nodule-associated bacterial community in Indian soils or high-throughput screening of common symbionts of pigeonpea is required for a crop that has evolutionary roots in India, the largest pigeonpea-producing nation. The present study analyzed the microbial diversity of pigeonpea root nodules in major Indian soils and evaluated the factors driving nodule microbiome assembly.

We assessed the bacterial diversity associated with pigeonpea in Indian soils to identify whether there is a core nodule microbiome. In addition, the influences of edaphic factors, geoclimatic conditions, host characteristics, and agricultural practices on the nodule microbiota and its spatial distribution across different pigeonpea genotypes were examined. We also compared the nodule microbiomes among multiple pigeonpea genotypes (considering three genotypes and three cultivated soil types) and their wild relative Cajanus scarabaeoides from their native soils. We hypothesized that the composition and diversity of the nodule bacterial community, especially the relative abundance of rhizobial and non-rhizobial taxa, may vary between the different genotypes and soil types and be influenced by environmental and host factors. Further, we presumed wild pigeonpea to host a distinct nodule microbiome due to its adaptation to native soil environments.

Methods

Seed material

Three popular pigeonpea cultivars (genotypes) with unique physiological and agronomical traits viz. Mannem Konda Kandi (MKK; ICPH-2740), Asha (ICPL-87119), and Durga (ICPL-84031), were selected for this study. The seeds were procured from the International Crop Research Institute for Semi-Arid Tropics (ICRISAT), Patancheru, Hyderabad, Telangana, India. Before sowing, seeds were surface-sterilized using 0.1% mercuric chloride (HgCl2) and 70% ethanol, then germinated on Murashige and Skoog agar medium.

Soil sampling

Soil samples were collected from pigeonpea fields across multiple locations to assess the impact of soil type and genotype on the pigeonpea nodule microbiome. Three distinct soil types—Alfisols (Rompicharla, Andhra Pradesh; 16.213900 N, 79.921386 E), Vertisols (Athner, Madhya Pradesh; 21.6406552 N, 77.91300 E), and Inceptisols (Sitamarhi, Uttar Pradesh; 25.2782289 N, 82.28691 E)—were sampled from farmers’ fields during the pre-sowing season in June 2017. Additional soil samples were collected from 92 locations across major pigeonpea-growing states in India (Andhra Pradesh, Telangana, Madhya Pradesh, and Uttar Pradesh). These regions include diverse soil types, geoclimatic conditions, and agricultural practices (e.g., cropping patterns, irrigation, and fertigation). Physicochemical characterization of all soils was performed at the Charles Renard Analytical Laboratory, ICRISAT, using standard analytical methods (Additional File 1).

Pigeonpea nodule microbiome in multiple genotypes

Three pigeonpea genotypes (Asha, Durga, and MKK) were grown in three different soil types (Alfisols, Vertisols, and Inceptisols) to assess the relative contributions of soil type and genotype to the nodule microbiome. Three plants of each genotype were transplanted into pots (7.5 kg capacity) filled with respective soil types. The plants were grown in six biological replicates in a glasshouse under identical light, temperature, and humidity conditions until the peak vegetative stage (four weeks). Six pots of soil for each soil type (without growing any plants) were used as bulk soil controls. The plants and control pots were watered as required with sterilized distilled water every other day without any further fertilization. Nodules were harvested at the vegetative stage (28 days after seedling emergence) and stored at -80 °C until DNA isolation.

Pigeonpea nodule microbiome in multiple soil types

A single pigeonpea genotype (MKK) was grown in 92 different soil types to study the influence of agroclimatic conditions on nodule microbiomes. The MKK genotype (ICPH-2740) was chosen based on its widespread cultivation and well-documented agronomic performance in diverse agroclimatic conditions. Plants were grown in three biological replicates in different soils in a glasshouse under identical light, temperature, and humidity conditions for up to 28 days until the plants formed mature nodules. The plants and bulk soil control (pre-sowing soils without plants) pots were supplied with sterilized distilled water, as needed, without fertilizer. The plants were harvested during the mature nodule formation stage (28th day). Nodules from three biological replicates were pooled together, and four bulk soil samples (control) were preserved at -80 °C until processing.

Spatial distribution of the nodule microbiome

To study the spatial distribution of nodule microbiomes and assess whether nodule microbiome composition varied with nodule position and host genotype, three pigeonpea genotypes (MKK, Asha, and Durga) were grown in a single soil type (red soil) in a glasshouse under identical conditions of light, temperature, and humidity for up to 28 days. Individual nodules were carefully harvested and stored at -80 °C until metagenomic DNA isolation. Nodules were classified into three groups based on their position on the root system: primary roots (taproots), secondary roots, and tertiary roots (Supplementary Fig. S1, Additional File 2). Three root nodules per plant were randomly selected across the entire root system, with one nodule each from primary, secondary, and tertiary roots. The microbiomes of 27 individual pigeonpea root nodules from nine plants (three plants of each genotype) and three bulk soil (pre-sowing) samples were analyzed via high-throughput DNA sequencing.

Wild pigeonpea (C. scarabaeoides) nodule microbiome

Wild pigeonpea (C. scarabaeoides) plants were sampled in situ from their natural habitat at the University of Hyderabad campus (17.457462 N, 78.314313 E) in Rangareddy district, Telangana, India, during the peak vegetative stage (60 days after seedling emergence). University of Hyderabad has a tropical environment with a dry deciduous biome and scrubby vegetation, with an average annual temperature of 24 °C (16.1–38.3 °C) and an average annual precipitation (rainfall) of 956.55 mm. Nodules were harvested manually from the wild pigeonpea roots and processed for DNA isolation.

Metagenomic DNA extraction, PCR, and sequencing

Metagenomic DNA was extracted from the bulk soil and nodules (0.3–0.5 g for each) using a NucleoSpin® Soil Kit (Machery Nagel, Germany) according to the manufacturer’s instructions. The V4 hypervariable region of the bacterial 16S rRNA gene was amplified using the 515F/806R primer pair [43]. The PCR mixture consisted of 0.2 µL of Phusion high fidelity (0.2 µL), 4 µL of high fidelity (HF) buffer (F520l; Thermo Scientific, Waltham, MA, USA), 0.4 µL of dinucleotide triphosphates (dNTPs), 1 µL of primers, 1.5 µL of template DNA (5 ng/µL), and 20 µL of H2O. For the nodule fraction, 1 µM peptide nucleic acid (PNA) was used to target plastid (pPNA, 5′-GGCTCAACCCTGGACAG-3′) and mitochondrial (mPNA, 5′-GGCAAGTGTTCTTCGGA-3′) DNA (PNA Bio, Newbury Park, CA, USA) as PCR clamps [44]. The PCR conditions were as follows: 98 °C for 1 min; 35 cycles of 98 °C for 30 s, 57 °C for 30 s, and 72 °C for 45 s; and a final elongation step of 72 °C for 7 min. Each DNA sample was amplified in triplicate, pooled, and then purified using a PCR clean-up kit (D4014, ZymoResearch). For each amplification run in 96-well plates, PCR-grade water was used as a negative control (no-DNA control). We proceeded with 16S rRNA gene diversity only because less than 1% of the samples were positive for ITS gene amplification for fungal community detection (data not shown). The samples were pooled and sequenced on the Illumina MiSeq platform via the V3 chemistry of 300PE run at M/s. Molecular Research DNA Laboratory in Texas, USA.

Processing of sequencing data

Initial quality filtering and read alignment were performed via USEARCH 10 fastq_mergepairs with fastq_maxee with an EE score of 1 [45]. After barcode removal, only reads above the desired length of 265 bp were used for further analysis. Reads were filtered from plant chloroplasts and mitochondria (approximately 2% of the initial reads were of plant origin) using a custom-made Bash script [46] (Additional File 3). The reads were binned into zero-radius operational taxonomic units (zOTUs), and chimeras were removed according to the Usearch10 pipeline with Unoise3 [47]. Bacterial zOTUs were annotated using the SILVA SSU132 16S rRNA database [48]. Additional details on the sequencing data processing, including quality control measures, read counts, and zOTU clustering procedures for all phases of this study, are provided in Additional File 2.

Bioinformatic and statistical analyses of microbiome

Comprehensive statistical, visual, and meta-analyses of the microbiome, including diversity analyses and comparisons along with graphical representation, were performed using MicrobiomeAnalyst [49, 50], an R-based online tool (https://www.microbiomeanalyst.ca/). The data were filtered for low count and low variance and normalized by cumulative sum scaling for marker gene (16S rRNA) analyses. The Shannon diversity index was used to measure the alpha diversity with the Mann–Whitney/Kruskal–Wallis (nonparametric tests) statistical method for significance testing. The Bray–Curtis index was calculated and visualized using principal coordinate analysis (PCoA) plots to assess beta diversity among the samples. The significance of the index was evaluated with permutational multivariate analysis of variance (PERMANOVA). The core taxa (phyla and genera) were visualized as heatmaps of compositional (relative) abundance. The key constituents of the core microbiome were identified at a relative abundance threshold of 0.01% and a sample prevalence of 20% [49, 50]. The unique and/or predictive features (biomarkers) were identified and classified using linear discriminant analysis (LDA) effect size (LEfSe) and random forest analysis. The classification errors (Out-of-bag [OOB] errors) were estimated to validate each random forest model.

The statistical significance of factors influencing the nodule microbial community was evaluated through permutations of residuals under a reduced model, the sum of squares type III (partial) with 9,999 permutations using unrestricted permutations of the raw data model of PERMANOVA. The pseudo-F values obtained from PERMANOVA served as proxies, indicating the relative importance of each factor in differentiating the samples. The pseudo-F values for each set of factors were plotted and visualized in Prism 9 (GraphPad, San Diego, USA).

Results

Influence of soil type and genotype on the pigeonpea nodule microbiome

The impact of soil type and genotype on the bacterial community associated with pigeonpea nodules was studied in three different pigeonpea genotypes grown in three distinct Indian soils. Nodules of different shapes, sizes, and quantities were obtained from pigeonpea plants (Supplementary Figs. S2-S3, Additional File 2).

The bulk soil fraction had a greater α-diversity (Shannon diversity index) than the pigeonpea nodules (Supplementary Fig. S4, Additional File 2). LEfSe analysis identified Bradyrhizobium in nodules as a biomarker associated with cultivated pigeonpea. Proteobacteria, Chloroflexi, and Planctomycetes were the predominant phyla in the nodules of pigeonpea. Bradyrhizobium, Unclassified_Roseiflexaceae, Unclassified_Planctomycetales, Unclassified_Rokubacteriales, and Unclassified_Comamonadaceae were the most abundant taxa at the genus level in the pigeonpea nodule microbiome. Among the genotypes, the hybrid genotype MKK exhibited the highest α-diversity (Shannon diversity index), followed by Durga and Asha (Supplementary Figs. S5-S7, Additional File 2).

LEfSe analysis identified Bradyrhizobium as a biomarker for the Asha genotype. Differential abundance at the genus level was observed across the selected genotypes. Bradyrhizobiumwas the most abundant genus among all the genotypes, but it was highly abundant in the Asha genotype. The soil type also influenced bacterial diversity. Vertisols presented the highest α-diversity (Shannon diversity index), followed by Inceptisols and Alfisols, although the α-diversity in Vertisols was lower than that in pre-sowing soils. In Alfisols and Vertisols, the dominant phyla were Proteobacteria, Chloroflexi, and Methylomirabilota, whereas in Inceptisols, Proteobacteria, Planctomycetes and Chloroflexi were the most abundant. All soils were dominated by Bradyrhizobium at the genus level, with a greater abundance detected in Alfisols (Supplementary Figs. S8-S9, Additional File 2).

The core nodule microbiome of multiple pigeonpea genotypes grown in three different soil types was dominated by Proteobacteria, Chloroflexi, and Planctomycetes at the phylum level. The key constituents of the core nodule microbiome across all the soils at detection thresholds (relative abundances) ranging from 0.6 to 1% included Bradyrhizobium, Unclassified_Roseiflexaceae, Unclassified_Rokubacteriales, Blastococcus, Unclassified_Planctomycetales, Unclassified_Comamonadaceae, Unclassified_Acidobacteria, Unclassified_Pirellulaceae, and Unclassified_Nitrosomonadaceae (Supplementary Fig. S10, Additional File 2). Analysis of similarity (ANOSIM) values for all ordered groups were calculated on the basis of the Bray‒Curtis similarity matrix and are depicted using PRIMER 7 (Supplementary Table S1, Additional File 2).

Edaphic factors shaping the pigeonpea nodule microbiome

Nodules from pigeonpea plants grown in laterite soil exhibited the highest α-diversity (Shannon diversity index), followed by those grown in black, red, alluvial, and mixed red and black soils (Supplementary Fig. S11, Additional File 2). LEfSe analysis revealed Limnobacter and Flavisolibacter as the differentially abundant bacterial genera in mixed black and red soils, and Novosphingobium in laterite soils. The bacterial genera Ideonella, Azospirillum, Flavobacterium, Ensifer, Priestia, and Curvibacter were identified as biomarkers for black soil; Novosphingobium, Caenimonas, and Aureimonas were identified for laterite soil; Limnobacter, Flavisolibacter, and Herbaspirillum were identified for mixed red and black soil; and Chitinophaga and Pseudomonas were identified for red soil. Amplicon sequencing revealed significant shifts in the relative abundance of these differentially abundant and biomarker bacterial genera identified by LEfSe and RF analyses (Supplementary Fig. S12, Additional File 2).

All the nodules across the soil types presented a greater relative abundance of several non-rhizobial endophytes (NREs), whereas the rhizobial taxa Ensifer, Bradyrhizobium, Microvirga, Mesorhizobium, and Shinella were less abundant (Supplementary Figs. S13-S14, Additional File 2). Nodules from plants grown in acidic soils presented the highest α-diversity (Shannon diversity index), followed by those from alkaline and neutral soils. LEfSe analysis identified Novosphingobium, Caenimonas, and Pseudomonas as differentially abundant in acidic soils and Limnobacter and Flavisolibacter as differentially abundant in alkaline soils. All the nodules, irrespective of the soil pH, presented greater relative abundances of several NRE taxa, whereas the rhizobial taxa Ensifer, Bradyrhizobium, Microvirga, Mesorhizobium, and Shinella were less abundant in the nodules. Notably, Ensifer had a greater relative abundance in alkaline soils than in acidic or neutral soils (Supplementary Figs. S15-S16, Additional File 2).

Nodules from plants grown in soils with moderate N levels (280–560 kg ha− 1) presented significantly greater α-diversity (Shannon diversity index) than those from soils with low N levels (< 280 kg ha− 1) did (Supplementary Figs. S17-S18, Additional File 2). Nodules from plants grown in soils with sufficient B concentrations (≥ 0.5 mg kg− 1 [0.5 ppm]) presented significantly greater α-diversity than those from soils deficient in B concentrations (< 0.5 mg kg− 1 [0.5 ppm]) (Supplementary Figs. S19-S20, Additional File 2). Among the nodules from plants grown in soils sampled from different geographical locations at varying altitudes (ranging from 50 m to 600 m above mean sea level), those from higher altitudes (≥ 300 m) exhibited slightly higher α-diversity than those from lower altitudes (< 300 m) (Supplementary Figs. S21-S22, Additional File 2). The α-diversity did not significantly differ between nodules from plants grown in soils with and without fertilizer application (Supplementary Figs. S23-S24, Additional File 2).

Among the nodules from plants grown in soils with different intercropping patterns, the pigeonpea–sunflower (Helianthus annuus) intercropping system exhibited the highest α-diversity, followed by the black gram (Vigna mungo), soybean (Glycine max), green gram (Vigna radiata), cotton (Gossypium herbaceum), horse gram (Macrotyloma uniflorum), and maize (Zea mays) intercropping systems. Solitary pigeonpea cropping resulted in greater α-diversity than intercropping with groundnut (Arachis hypogaea), guava (Psidium guajava), sorghum (Sorghum bicolor), cucumber (Cucumis sativus), or ridge gourd (Luffa acutangula). Ensifer was more abundant in the bulk soil, whereas Pseudomonas was relatively abundant in nodules across multiple intercropping systems (Supplementary Figs. S25-S26, Additional File 2).

Analysis of the core microbiome enabled the identification of core taxa that remained unchanged in composition across different sample groups on the basis of sample prevalence and relative abundance. Pigeonpea nodules were colonized primarily by Proteobacteria and Bacteroidetes (Fig. 1A). The main bacterial genera found in pigeonpea nodules included Limnobacter, Novosphingobium, Flavisolibacter, Caenimonas, Pseudomonas, Chitinophaga, Ensifer, Unclassified_Aurantimonadaceae, Bradyrhizobium, and Microvirga. These genera constituted the core nodule microbiome across all the studied soils, with detection thresholds (relative abundance %) ranging from 0.010 to 0.728 (Fig. 1B).

Fig. 1
figure 1

Heatmaps representing the core nodule microbiome of pigeonpea across diverse soils at the (A) phylum and (B) genus levels. The Y-axis represents the prevalence level of core bacterial taxa across the detection threshold (relative abundance) range on the X-axis. The variation in the prevalence of each phylum/genus is indicated by a color gradient from blue (decreased) to red (increased). The ‘Uncultured_’ taxa label in the figure represents unclassified bacterial taxa

Multifactor PERMANOVA analysis revealed that soil type, pH, macronutrient, and micronutrient levels were the most significant factors influencing the bacterial community composition in pigeonpea nodules (p < 0.05) (Fig. 2A–E and Supplementary Table S2, Additional File 2). Geoclimatic factors and agricultural practices also had notable but lesser impacts. The results regarding the influence of different factors, which were not included above, are provided in Supplementary Figs. S27-S56 (Additional File 2).

Fig. 2
figure 2

PERMANOVA output measuring the influence of different factors on the pigeonpea nodule microbiota using the pseudo-F values as proxies. (A) Edaphic factors, (B) Soil macronutrients, (C) Soil micronutrients, (D) Agricultural factors, (E) Geoclimatic factors, and (F) Host factors. * p < 0.05; ** p < 0.01; *** p < 0.001; ns p > 0.05

Single-nodule microbiomes of pigeonpea genotypes: spatial distribution across primary, secondary, and tertiary roots

The bulk soil exhibited significantly greater α-diversity (Shannon diversity index) than the nodules, regardless of their position on the roots. Among the nodule positions, primary and tertiary root nodules presented very similar α-diversity, with primary root nodules having a marginally higher Shannon index than did tertiary root nodules. In contrast, secondary root nodules had reduced α-diversity. The PCoA plots revealed distinct clustering patterns for each group, and PERMANOVA indicated significant differences in zOTU assemblage between nodule bacterial communities on the basis of their position on the roots. All the nodules presented a relatively high abundance of the rhizobial genera Bradyrhizobium and Rhizobium. Other rhizobial taxa, such as Shinella, Mesorhizobium, Microvirga, and Ensifer, were also abundant in the nodules, although to a lesser extent. The taxonomic composition of the nodule bacterial communities was influenced by both the host genotype and the position of the nodule on the roots (Supplementary Figs. S57-S60, Additional File 2).

Multifactor PERMANOVA analysis (Supplementary Table S3, Additional File 2) revealed that host genotype (pseudo-F = 4.1087) exerted the strongest influence on the nodule bacterial community composition, indicating greater variation between sample groups than within each group. This was followed by the influence of plant fraction (nodule and bulk soil; pseudo-F = 3.2245) and nodule position on the root (pseudo-F = 2.1938) (Fig. 2F). The interaction between plant genotype and nodule position had a significantly stronger influence than the individual factors did (pseudo-F = 4.2787).

ANOSIM revealed significant differences in the overall nodule bacterial community structures across pigeonpea genotypes and nodule positions. Among the genotypes, the most pronounced community variations were observed between MKK and Durga (ANOSIM R = 0.499, P < 0.001), followed by Durga and Asha (R = 0.306, P < 0.05), whereas there was less difference between Asha and MKK (R = 0.248, P < 0.01). Significant bacterial community variation was also detected within each genotype across different nodule positions (primary, secondary, and tertiary roots) (Supplementary Table S4, Additional File 2). A schematic representation of the variation in nodule microbial community composition between nodule positions across the genotypes (based on ANOSIM R values) is presented in Fig. 3.

Fig. 3
figure 3

Schematic representation of analysis of similarities (ANOSIM) between the pigeonpea nodule microbial communities across the three different genotypes. Numerals represent the percentage of variance in the nodule microbial community composition explained by the analyzed factors based on ANOSIM R values. The dotted lines in black represent the differences between nodule positions within each genotype, while the red dotted lines represent the overall variation in nodule microbial communities between the genotypes. Higher percentages indicate a stronger influence of the respective factor (nodule position or genotype) on the microbial community composition. The distances depicted by the dotted lines are representative and do not correspond to actual spatial distances in relation to the percentage variations

Nodule microbiome of wild pigeonpea (Cajanus scarabaeoides)

The core nodule microbiome of the wild pigeonpea (Supplementary Fig. S61, Additional File 2) was predominantly composed of Proteobacteria, Methylomirabilota, and Actinobacteria at the phylum level, with Bradyrhizobium being the most prevalent at the genus level (Fig. 4). A diverse range of rhizobia and NREs were observed in the core microbiome.

Fig. 4
figure 4

The core microbiome of wild pigeonpea nodules at the (A) phylum and (B) genus levels. The Y-axis represents the prevalence level of core bacterial taxa across the detection threshold (relative abundance) range on the X-axis. The variation in the prevalence of each phylum/genus is indicated by a color gradient from blue (increased) to yellow (decreased). The ‘Not_assigned/Uncultured_’ taxa label in the figure represents unclassified bacterial taxa

Comparative analysis of the nodule microbiomes of wild and cultivated pigeonpea genotypes using ANOSIM (Table 1) revealed significant differences in the bacterial community structure, with the wild genotype exhibiting a distinct microbiome compared to its cultivated counterpart. The most significant difference of 85.1% was observed between the wild and cultivated genotypes, with a total difference of 66.1% at p < 0.001. Additionally, a considerable difference of 35.3% was observed between the control and the wild type, which was still significant at p < 0.05. Principal coordinate analysis further supported the ANOSIM results, indicating that 28.2% of the variation in the bacterial communities was explained by the X-axis and 11.9% by the Y-axis (Fig. 5).

Table 1 ANOSIM differences of the nodule microbiome between pigeonpea genotypes. Comparative differences were calculated using a Bray‒Curtis similarity distance matrix by ANOSIM (PRIMER 7 software).
Fig. 5
figure 5

Principal coordinate analysis for the wild vs. cultivated pigeonpea nodule microbiome. The samples were transformed by the square root and distance matrix constructed with the Bray‒Curtis similarity metric. Total variations of 28.2% (X-axis) and 11.9% (Y-axis) were observed with principal coordinate analysis

Discussion

Nitrogen-fixing bacteria form symbiotic partnerships with leguminous plants (and non-legume Parasponia spp.) in nodules, where atmospheric nitrogen is reduced to ammonia, promoting plant growth and development. Plants exchange photosynthetic assimilates in the form of organic carbon to provide nutrients and support symbiotic bacteria, apart from the niche. They secrete flavonoids into the rhizosphere to recruit these rhizobia, leading to signal exchange, nodulation, and symbiosis [51, 52]. This symbiotic interaction occurs only between specific legume(s) and their corresponding symbionts(s).

Globally, pigeonpea and other important legume crops are rich sources of protein, nutrients, and other benefits to the rapidly increasing human population while being more sustainable than most other crops are [3]. Rhizobial endosymbionts associated with pigeonpea nodules supply the nitrogen required for plant growth [14]. Notably, legume root nodules host a diverse community of bacteria rather than a single species of rhizobia [10]. Pigeonpea nodules harbor diverse rhizobia and non-rhizobial endophyte (NRE) taxa [21, 25, 26].

The nodule bacterial community is shaped by the host genotype and nodule position on the root

Our results demonstrated considerable differences in the overall nodule bacterial community structure across the three genotypes and between nodule positions within each genotype. The nodule microbiota was similar at the phylum level but more distinct at the genus level in terms of the relative abundance of each taxon. The variation in the nodule bacterial community structure across various nodule positions within individual genotypes could be due to a random process since close root nodules show significant differences in their microbiome composition [53]. The host genotype can also influence the nodulation process [54].

Culture-dependent studies have shed light on the potential benefits of NREs on their host plants [53, 55,56,57]. However, their preferential selection by legume hosts is not known. Culture-independent (metagenomic) studies have demonstrated the vast diversity and abundance of NREs in legume nodules [53, 55, 58]. Our study found that many non-rhizobial taxa were relatively abundant in the nodules. It is unclear whether their presence is linked to potential biological and/or biogeochemical cues resulting from the tripartite interactions between the soil, host plants, and indigenous microbiota.

Edaphic and agro-climatic factors shape the Pigeonpea nodule bacterial community

Soil parameters and history are key drivers of microbial composition and diversity. Soil type, and to a lesser extent, host genotype, influence the abundance and structural and functional diversity of soil and root-associated microbial communities [31, 59,60,61]. We found that edaphic factors play a significant role and explain most of the variation in the bacterial community structure of the pigeonpea nodule microbiome. Variations in soil type, pH, and nutrient status were strongly associated with differences in nodule bacterial diversity and abundance, confirming previous studies that highlighted the influence of soil properties and host genotype on the composition of root-associated microbiomes [10, 62,63,64,65,66,67,68].

Our findings also indicate that agricultural practices significantly affect the nodule microbiome. Agricultural practices significantly influence the diversity, abundance, and richness of the bacterial community in the pigeonpea rhizosphere and drive the dynamics of dominant rhizobacterial taxa [69]. Although fertilizer amendments are known to influence the soil properties and microbial community composition [70,71,72,73,74], we did not observe a significant effect of fertilizer amendment (in the soil) on the pigeonpea nodule microbial community. However, significant variation in the relative abundance of the nodule microbiota was observed across the different cropping systems, particularly intercropping, likely due to increased root proximity and altered root exudate profiles that promote beneficial microbial recruitment [72, 75,76,77,78,79].

Additionally, climatic factors such as rainfall [62, 63] and geographical factors, including altitude and location [80,81,82], were correlated with variations in the nodule microbiome. Notably, location-driven biogeographical patterns in nodule bacterial communities, previously documented in soybean [83, 84], were also evident in pigeonpea, underscoring the complex interplay of environmental and agronomic factors in modulating these communities.

Pigeonpea nodules are colonized by Proteobacteria and Bacteroidetes

Bulk soil exhibited greater microbial diversity than pigeonpea nodules, as reflected by the α-diversity and bacterial community composition data. The soil serves as a reservoir or meta-community for microorganisms that comprise the nodule microbiomes. Our findings align with those of Miranda-Sánchez et al. [85] and Xiao et al. [55], who reported that legume hosts recruit nodule endophytes from bulk soil through a hierarchical filtering process, wherein a gating mechanism at the root possibly restricts the microbiota from entering the nodule endosphere. Root-associated microbiomes, especially nodule microbiomes, are typically less diverse than soil microbiomes because of plant selection and competition among microorganisms for plant-derived resources [51, 55]. Nevertheless, nodule microbiomes can still be extremely complex, containing hundreds of distinct bacterial species [86].

Pigeonpea nodules showed a specific selection for Proteobacteria and Bacteroidetes, regardless of geographical origin, soil nutrient status, soil type, or soil history. Members of these phyla are often abundant in plant tissues [87, 88]. Proteobacteria, especially Alphaproteobacteria, are common root colonizers across multiple soils and plant species, including barley, rice, Arabidopsis, and Lotus [89,90,91,92]. Alphaproteobacteria members likely metabolize plant-derived nutrients quickly [88] and share genomic similarities with plant symbionts [93].

The main nodule-inhabiting genera representing the core microbiome across diverse soils were Limnobacter, Novosphingobium, Flavisolibacter, Caenimonas, Pseudomonas, Chitinophaga, Ensifer, Bradyrhizobium, and Microvirga. This suggests that these genera are particularly attracted to pigeonpea roots, possibly through active plant secretion. Plants influence the microbiota around their roots, creating a gradient of impact that decreases with increasing distance from the roots [46]. While plant presence shapes the surrounding microbial community, other factors, such as soil type, determine its exact profile.

Diverse rhizobia are associated with pigeonpea nodules

We consistently detected at least three rhizobial genera across all pigeonpea nodules: Ensifer, Bradyrhizobium, and Microvirga. In contrast, Mesorhizobium and Shinella were less abundant. Among these, Ensifer was the most abundant nodule-colonizing rhizobial genus, as identified through our 16S rRNA gene amplicon assay. Ensifer spp. are common endosymbionts of native Indian legumes that grow in alkaline soils [94,95,96]. We observed a relatively high abundance of Ensifer in the nodules of pigeonpea grown in alkaline soils, suggesting that Ensifer might outcompete Bradyrhizobium, the most common pigeonpea endosymbiont [14, 21], in these native soils to colonize pigeonpea nodules. Moreover, the case of pigeonpea is not limited, as soybean roots, which typically host Bradyrhizobium, are colonized by a bacterial community where this symbiont constitutes merely about 1% of the total population [60, 97]. This relatively low symbiont abundance contrasts with that of other legumes, such as pea, where the symbiotic genus Rhizobium accounts for approximately 10–20% of the root microbiome; Medicago, with Ensifer accounting for approximately 10–60%; and Lotus, where Mesorhizobium accounts for approximately 10% of the root-associated bacteria [91, 98,99,100,101].

Bradyrhizobium, abundant in Indian soils [14], poorly colonizes pigeonpea roots under competition from other soil bacteria and likely includes many non-symbiotic strains. Nevertheless, pigeonpea plants can still nodulate [46]. Bradyrhizobium has evolved to be recognized by pigeonpea, infect its roots, and develop root nodules. Mesorhizobium association with pigeonpea nodules has been well documented in Indian soils [15, 17]. Nevertheless, the detection of zOTUs belonging to atypical rhizobia, such as Shinella and Microvirga in pigeonpea nodules, which are potential endosymbiotic bacteria associated with other legumes, indicates that these atypical taxa can be endophytes in non-host legumes such as pigeonpea [55, 102].

The core nodule microbiome of wild pigeonpea is dominated by Bradyrhizobium

Cajanus scarabaeoides (L.) Millspaugh is one of the few wild relatives of pigeonpea that can reproduce. It is highly tolerant to drought and salinity, is mostly resistant to insect pests, and has a high protein content in its grains [103]. Widely distributed in South Asia, this species has been largely overlooked. Unlike contemporary agricultural systems, which are artificially molded by human intervention and rely on fertigation and chemigation to maintain high yields, wild ecosystems offer greater genetic variety, soil heterogeneity, interspecies competition, and biodiversity [104]. Modern technology allows the creation of high-yielding varieties that incorporate microbiota from their wild relatives [105]. However, plant‒soil feedback in agricultural environments, particularly the effects of plant domestication on soil ecosystems and geochemical processes, remains underexplored.

We found that Bradyrhizobium dominates the core microbiome of wild pigeonpea nodules. Bradyrhizobium is associated with pigeonpea and other legumes such as cowpea, black gram, and green gram [106, 107]. In the Dominican Republic, Bradyrhizobium strains form a nitrogen-fixing symbiosis with pigeonpea and are effective biofertilizers to replace N fertigation [19]. Bradyrhizobium sp. is associated with high pigeonpea yields in Côte d’Ivoire [21]. Biological nitrogen fixation in groundnut was improved by B. yuanmingense isolated from Ghanaian native soils [108]. Similarly, pigeonpea performance was enhanced by inoculation with biochar-formulated Bradyrhizobium strains [109]. In our previous work, we observed poor competitiveness of Bradyrhizobium toward the cultivated pigeonpea root microbiome, as both cultivated and wild pigeonpea nodules were dominated by Bradyrhizobium, revealing specific recruitment by all genotypes [46]. This selective recruitment ensures that the plant prioritizes microbial partners that provide essential functions, such as nitrogen fixation, even if these partners are not the most competitive in the broader soil environment.

Domestication significantly alters plant physiology and microbiomes [29, 110, 111]. The microbiota associated with the nodules of wild legumes are phylogenetically diverse and often possess various plant growth promoting abilities, including osmotolerance [112]. They may also represent unique natural hotspots of antibiotic-resistance genes [113]. Interestingly, a study on nodule bacterial diversity conducted in continental Portugal with 16 wild legume species revealed that most isolated strains belong to non-rhizobial genera like Pseudomonas and Flavobacterium, rather than typical rhizobia [112]. Wild soybeans, which are ecologically more resilient than cultivated soybeans, may have evolved to recruit beneficial bacteria in their rhizosphere that may increase nutrient uptake, biostasis, and disease resistance [114]. These findings have important implications for understanding legume–rhizobia interactions and optimizing symbiotic nitrogen fixation in agriculture.

Limitations of the study

Our study primarily focused on spatial variation in the pigeonpea nodule microbiome. We did not extensively analyze temporal dynamics and seasonal variations in microbial community composition, which could influence symbiotic interactions, microbial succession, and stability. Additionally, the study’s focus on specific genotypes and soil types limits the broader applicability of the findings across diverse agroecosystems. The research was limited to four genotypes, including three cultivated varieties and one wild relative (Cajanus scarabaeoides), while many pigeonpea cultivars and wild accessions remain unexplored. Future studies should include these wild accessions, modern hybrids, and landraces, which may harbor unique microbial associations beneficial for crop resilience and symbiotic nitrogen fixation. Although edaphic factors were considered, other influences, such as plant physiological status and root exudate composition, were not examined. Addressing these gaps in future research could provide a more comprehensive understanding of the factors shaping the pigeonpea nodule microbiome.

Conclusions

In conclusion, this study elucidates the complex interactions between host genotype, nodule position, and edaphic factors in shaping the nodule microbiome of pigeonpea. Our findings reveal that soil properties, particularly pH, nutrient status, and soil type, exert a substantial influence on the nodule bacterial community than the host genotype, underscoring the critical role of environmental conditions in symbiotic relationships. The core nodule microbiome, dominated by Proteobacteria and Bacteroidetes, with Bradyrhizobium and Ensifer as the major rhizobial taxa and various non-rhizobial taxa, reflects the selective recruitment by the host plant. While bulk soil exhibited a higher microbial diversity, the nodule microbiome was less diverse but highly specialized, indicating a selective process driven by plant-microbe interactions. Comparative analysis between wild and cultivated pigeonpea revealed substantial differences in their nodule microbiomes, underscoring the influence of domestication on microbial community composition. These insights into the microbial diversity and community structure of pigeonpea nodules provide a foundation for developing strategies to boost crop productivity and resilience through targeted microbial management. Future research, with a focus on the functional characterization of these microbial communities and their interactions with the host plant, would help to exploit microbial potential in sustainable agriculture.

Data availability

The datasets generated and/or analysed during the current study are available in the SRA database of NCBI under BioProject accession number PRJNA1171966.

Abbreviations

ANOSIM:

Analysis of similarity

LDA:

Linear discriminant analysis

LefSe:

Linear discriminant analysis (LDA) effect size

MKK:

Mannem Konda Kandi (ICPH–2740)

NRE:

Non–rhizobial root nodule endophytes

OOB:

Out–of–bag

PCoA:

Principal coordinate analysis

PERMANOVA:

Permutational multivariate analysis of variance

PNA:

Peptide nucleic acid

zOTU:

Zero–radius operational taxonomic unit

References

  1. Online Agriculture FAOSTAT. Statistics 2021. https://www.fao.org/faostat/en/#data

  2. Fuller DQ, Murphy C, Kingwell-Banham E, Castillo CC, Naik S. Cajanus cajan (L.) Millsp. Origins and domestication: the South and Southeast Asian archaeobotanical evidence. Genet Resour Crop Evol. 2019;66:1175–88. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10722-019-00774-w.

    Article  CAS  Google Scholar 

  3. Varshney RK, Saxena RK, Upadhyaya HD, Khan AW, Yu Y, Kim C, et al. Whole-genome resequencing of 292 Pigeonpea accessions identifies genomic regions associated with domestication and agronomic traits. Nat Genet. 2017;49:1082–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ng.3872.

    Article  CAS  PubMed  Google Scholar 

  4. Tiwari AK, Shivhare AK. Pulses in India: retrospect and prospects– 2018. Directorate of pulses development. Ministry of Agriculture and Farmers Welfare, Government of India; 2019.

  5. Singh RP. Status paper on pulses. Directorate of pulses.Development. Ministry of Agriculture and Farmers Welfare, Government of India; 2013.

  6. Siddiqui SA, Fatima N. Indian soils: identification and classification. Earth Sci India. 2017;10:1–14.

    CAS  Google Scholar 

  7. Mendoza-Suárez M, Andersen SU, Poole PS, Sánchez-Cañizares C. Competition, nodule occupancy, and persistence of inoculant strains: key factors in the Rhizobium-legume symbioses. Front Plant Sci. 2021;12:1684. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2021.690567.

    Article  Google Scholar 

  8. Mendoza-Suárez MA, Geddes BA, Sánchez-Cañizares C, Ramírez-González RH, Kirchhelle C, Jorrin B et al. Optimizing Rhizobium- legume symbioses by simultaneous measurement of rhizobial competitiveness and N 2 fixation in nodules. Proceedings of the National Academy of Sciences. 2020;117:9822–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.1921225117

  9. Martínez-Hidalgo P, Hirsch AM. The nodule microbiome: N2-fixing rhizobia do not live alone. Phytobiomes J. 2017;1:70–82. https://doiorg.publicaciones.saludcastillayleon.es/10.1094/PBIOMES-12-16-0019-RVW.

    Article  Google Scholar 

  10. Custódio V, Gonin M, Stabl G, Bakhoum N, Oliveira MM, Gutjahr C, et al. Sculpting the soil microbiota. Plant J. 2022;109:508–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/tpj.15568.

    Article  CAS  PubMed  Google Scholar 

  11. Keet JH, Ellis AG, Hui C, Le Roux JJ. Legume-rhizobium symbiotic promiscuity and effectiveness do not affect plant invasiveness. Ann Bot. 2017;119:1319–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/aob/mcx028.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. De Meyer SE, De Beuf K, Vekeman B, Willems A. A large diversity of non-rhizobial endophytes found in legume root nodules in Flanders (Belgium). Soil Biol Biochem. 2015;83:1–11. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.soilbio.2015.01.002.

    Article  CAS  Google Scholar 

  13. Alaswad A, Oehrle N, Krishnan H. Classical soybean (Glycine max (L.) Merr) symbionts, Sinorhizobium fredii USDA191 and Bradyrhizobium diazoefficiens USDA110, reveal contrasting symbiotic phenotype on pigeon pea (Cajanus Cajan (L.) Millsp). Int J Mol Sci. 2019;20:1091. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms20051091.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Jorrin B, Maluk M, Atoliya N, Kumar SC, Chalasani D, Tkacz A, et al. Genomic diversity of pigeon pea (Cajanus Cajan L. Millsp.) endosymbionts in India and selection of potential strains for use as agricultural inoculants. Front Plant Sci. 2021;12:680981. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2021.680981.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Jain D, Kumari A, Saheewala H, Sanadhya S, Maheshwari D, Meena RH, et al. Biochemical, functional and molecular characterization of pigeon pea rhizobia isolated from semi-arid regions of India. Arch Microbiol. 2020;202:1809–16. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00203-020-01904-0.

    Article  CAS  PubMed  Google Scholar 

  16. Gosai J, Anandhan S, Bhattacharjee A, Archana G. Elucidation of quorum sensing components and their role in regulation of symbiotically important traits in Ensifer nodulating pigeon pea. Microbiol Res. 2020;231:126354. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.micres.2019.126354.

    Article  CAS  PubMed  Google Scholar 

  17. Singha B, Mazumder PB, Pandey P. Characterization of plant growth promoting rhizobia from root nodule of two legume species cultivated in Assam, India. Proceedings of the National Academy of Sciences, India Section B: Biological Sciences. 2018;88:1007–16. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s40011-016-0836-6

  18. Singh K, Dhull S, Gera R. Genetic diversity of diazotrophs nodulating pigeon pea in arid and semi-arid zones of Haryana, India. Int J Curr Microbiol Appl Sci. 2018;7:3447–67. https://doiorg.publicaciones.saludcastillayleon.es/10.20546/ijcmas.2018.703.397.

    Article  CAS  Google Scholar 

  19. Araujo J, Díaz-Alcántara CA, Velázquez E, Urbano B, González-Andrés F. Bradyrhizobium yuanmingense related strains form nitrogen-fixing symbiosis with Cajanus cajan L. in Dominican Republic and are efficient biofertilizers to replace N fertilization. Sci Hortic. 2015;192:421–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scienta.2015.06.009.

    Article  Google Scholar 

  20. Araújo J, Flores-Félix JD, Igual JM, Peix A, González-Andrés F, Díaz-Alcántara CA, et al. Bradyrhizobium cajani Sp. Nov. Isolated from nodules of Cajanus cajan. Int J Syst Evol Microbiol. 2017;67:2236–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1099/ijsem.0.001932.

    Article  CAS  PubMed  Google Scholar 

  21. Fossou RK, Ziegler D, Zézé A, Barja F, Perret X. Two major clades of bradyrhizobia dominate symbiotic interactions with Pigeonpea in fields of Côte D’Ivoire. Front Microbiol. 2016;7:1793. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2016.01793.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Rufini M, Oliveira DP, Trochmann A, Soares BL, de Andrade MJB, de Moreira FM. Bradyrhizobium spp. Strains in symbiosis with pigeon pea Cv. fava-larga under greenhouse and field conditions. Rev Bras Cienc Solo. 2016;40:e0160156. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/18069657rbcs20160156.

    Article  CAS  Google Scholar 

  23. Tulu D, Endalkachew W, Makka A, Asnake F, Tilahun A, Chris O. Morphophysiological diversity of rhizobia nodulating pigeon pea (Cajanus Cajan L. Millsp.) growing in Ethiopia. Afr J Biotechnol. 2018;17:167–77. https://doiorg.publicaciones.saludcastillayleon.es/10.5897/ajb2017.16338.

    Article  CAS  Google Scholar 

  24. López-Baena FJ, Ruiz-Sainz JE, Rodríguez-Carvajal MA, Vinardell JM. Bacterial molecular signals in the Sinorhizobium fredii-soybean symbiosis. Int J Mol Sci. 2016;17:755. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms17050755.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Bind M, Nema S. Isolation and molecular characterization of endophytic bacteria from pigeon pea along with antimicrobial evaluation against Fusarium Udum. Appl Microbiol Open Access. 2019;5:163. https://doiorg.publicaciones.saludcastillayleon.es/10.35248/2471-9315.19.5.163.

    Article  Google Scholar 

  26. Dubey RC, Maheshwari DK, Kumar H, Choure K. Assessment of diversity and plant growth promoting attributes of rhizobia isolated from Cajanus cajan L. Afr J Biotechnol. 2010;9:8619–29. https://doiorg.publicaciones.saludcastillayleon.es/10.4314/ajb.v9i50.

    Article  CAS  Google Scholar 

  27. Varshney RK, Chen W, Li Y, Bharti AK, Saxena RK, Schlueter JA, et al. Draft genome sequence of Pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers. Nat Biotechnol. 2012;30:83–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nbt.2022.

    Article  CAS  Google Scholar 

  28. Singh NK, Gupta DK, Jayaswal PK, Mahato AK, Dutta S, Singh S, et al. The first draft of the Pigeonpea genome sequence. J Plant Biochem Biotechnol. 2012;21:98–112. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s13562-011-0088-8.

    Article  PubMed  Google Scholar 

  29. Pérez-Jaramillo JE, de Hollander M, Ramírez CA, Mendes R, Raaijmakers JM, Carrión VJ. Deciphering rhizosphere Microbiome assembly of wild and modern common bean (Phaseolus vulgaris) in native and agricultural soils from Colombia. Microbiome. 2019;7:114. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-019-0727-1.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Stopnisek N, Shade A. Persistent Microbiome members in the common bean rhizosphere: an integrated analysis of space, time, and plant genotype. ISME J. 2021;15:2708–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41396-021-00955-5.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Lagunas B, Richards L, Sergaki C, Burgess J, Pardal AJ, Hussain RMF, et al. Rhizobial nitrogen fixation efficiency shapes endosphere bacterial communities and Medicago truncatula host growth. Microbiome. 2023;11:146. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-023-01592-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Fields B, Moeskjær S, Deakin WJ, Moffat EK, Roulund N, Andersen SU, et al. Rhizobium nodule diversity and composition are influenced by clover host selection and local growth conditions. Mol Ecol. 2023;32:4259–77. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/mec.17028.

    Article  CAS  PubMed  Google Scholar 

  33. Wang L, Jiao Y, Bi Y, Hu Y, Jiang Y, Wang S, et al. Nodulation number tempers the relative importance of stochastic processes in the assembly of soybean root-associated communities. ISME Commun. 2023;3. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s43705-023-00296-8.

  34. Debnath S, Chakraborty S, Langthasa M, Choure K, Agnihotri V, Srivastava A, et al. Non-rhizobial nodule endophytes improve nodulation, change root exudation pattern and promote the growth of lentil, for prospective application in fallow soil. Front Plant Sci. 2023;14. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2023.1152875.

  35. Hossain MS, DeLaune PB, Gentry TJ. Microbiome analysis revealed distinct microbial communities occupying different sized nodules in field-grown peanut. Front Microbiol. 2023;14. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2023.1075575.

  36. Ai J, Yu T, Liu X, Jiang Y, Wang E, Deng Z-S. Seed associated microbiota and vertical transmission of bacterial communities from seed to nodule in Sophora davidii. Plant Soil. 2023;491:285–302. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-023-06115-2.

    Article  CAS  Google Scholar 

  37. Yuan Y, Chen Z, Huang X, Wang F, Guo H, Huang Z, et al. Comparative analysis of nitrogen content and its influence on actinorhizal nodule and rhizospheric microorganism diversity in three Alnus species. Front Microbiol. 2023;14. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2023.1230170.

  38. Dove R, Wolfe ER, Stewart NU, Coleman A, Chavez SH, Ballhorn DJ. Root nodules of red alder (Alnus rubra) and Sitka alder (Alnus viridis Ssp. sinuata) are inhabited by taxonomically diverse cultivable microbial endophytes. Microbiologyopen 2024;13. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/mbo3.1422

  39. Ali R, Chaluvadi SR, Wang X, Hazzouri KM, Sudalaimuthuasari N, Rafi M, et al. Microbiome properties in the root nodules of Prosopis cineraria, a leguminous desert tree. Microbiol Spectr. 2024;12. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/spectrum.03617-23.

  40. Adan IH, Asudi GO, Niassy S, Jalloh AA, Mutua JM, Chidawanyika F, et al. Comparative Microbiome diversity in root-nodules of three Desmodium species used in push-pull cropping system. Front Microbiol. 2024;15. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2024.1395811.

  41. Matsumoto BL, Sieradzki ET, Greenlon A, Perilla-Henao LM, Ter Horst AM, Geonczy SE, et al. Viruses of nitrogen-fixing Mesorhizobium bacteria in globally distributed Chickpea root nodules. Phytobiomes J. 2024;8:216–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1094/PBIOMES-06-23-0042-R.

    Article  Google Scholar 

  42. Liu H, Ni B, Duan A, He C, Zhang J. High Frankia abundance and low diversity of microbial community are associated with nodulation specificity and stability of sea Buckthorn root nodule. Front Plant Sci. 2024;15. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2024.1301447.

  43. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci. 2011;108:4516–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.1000080107.

  44. Lundberg DS, Yourstone S, Mieczkowski P, Jones CD, Dangl JL. Practical innovations for high-throughput amplicon sequencing. Nat Methods. 2013;10:999–1002. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nmeth.2634.

    Article  CAS  PubMed  Google Scholar 

  45. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/bioinformatics/btq461.

    Article  CAS  PubMed  Google Scholar 

  46. Chalasani D, Basu A, Pullabhotla SVSRN, Jorrin B, Neal AL, Poole PS, et al. Poor competitiveness of Bradyrhizobium in pigeon pea root colonization in Indian soils. mBio. 2021;12:e00423–21. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/mBio.00423-21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Edgar R. UNOISE2: improved error-correction for illumina 16S and ITS amplicon sequencing. BioRxiv. 2016;081257. https://doiorg.publicaciones.saludcastillayleon.es/10.1101/081257.

  48. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gks1219.

    Article  CAS  PubMed  Google Scholar 

  49. Dhariwal A, Chong J, Habib S, King IL, Agellon LB, Xia J. MicrobiomeAnalyst: a web-based tool for comprehensive statistical, visual and meta-analysis of Microbiome data. Nucleic Acids Res. 2017;45:W180–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkx295.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Chong J, Liu P, Zhou G, Xia J. Using Microbiomeanalyst for comprehensive statistical, functional, and meta-analysis of Microbiome data. Nat Protoc. 2020;15:799–821. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41596-019-0264-1.

    Article  CAS  PubMed  Google Scholar 

  51. Clúa J, Roda C, Zanetti M, Blanco F. Compatibility between legumes and rhizobia for the establishment of a successful nitrogen-fixing symbiosis. Genes (Basel). 2018;9:125. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes9030125.

    Article  CAS  PubMed  Google Scholar 

  52. Poole P, Ramachandran V, Terpolilli J. Rhizobia: from saprophytes to endosymbionts. Nat Rev Microbiol. 2018;16:291–303. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nrmicro.2017.171.

    Article  CAS  PubMed  Google Scholar 

  53. Mayhood P, Mirza BS. Soybean root nodule and rhizosphere microbiome: distribution of rhizobial and nonrhizobial endophytes. Appl Environ Microbiol. 2021;87:e02884–20. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/AEM.02884-20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Plett KL, Bithell SL, Dando A, Plett JM. Chickpea shows genotype-specific nodulation responses across soil nitrogen environment and root disease resistance categories. BMC Plant Biol. 2021;21:310. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-021-03102-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Xiao X, Chen W, Zong L, Yang J, Jiao S, Lin Y, et al. Two cultivated legume plants reveal the enrichment process of the Microbiome in the rhizocompartments. Mol Ecol. 2017;26:1641–51. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/mec.14027.

    Article  CAS  PubMed  Google Scholar 

  56. Sharaf H, Rodrigues RR, Moon J, Zhang B, Mills K, Williams MA. Unprecedented bacterial community richness in soybean nodules vary with cultivar and water status. Microbiome. 2019;7:63. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-019-0676-8.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Etesami H. Root nodules of legumes: A suitable ecological niche for isolating non-rhizobial bacteria with biotechnological potential in agriculture. Curr Res Biotechnol. 2022;4:78–86. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.crbiot.2022.01.003.

    Article  CAS  Google Scholar 

  58. Wigley K, Moot D, Wakelin SA, Laugraud A, Blond C, Seth K, et al. Diverse bacterial taxa inhabit root nodules of Lucerne (Medicago sativa L.) in new Zealand pastoral soils. Plant Soil. 2017;420:253–62. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-017-3395-6.

    Article  CAS  Google Scholar 

  59. Sneha GR, Swarnalakshmi K, Sharma M, Reddy K, Bhoumik A, Suman A, et al. Soil type influence nutrient availability, microbial metabolic diversity, eubacterial and Diazotroph abundance in Chickpea rhizosphere. World J Microbiol Biotechnol. 2021;37:167. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11274-021-03132-0.

    Article  CAS  PubMed  Google Scholar 

  60. Liu F, Hewezi T, Lebeis SL, Pantalone V, Grewal PS, Staton ME. Soil Indigenous Microbiome and plant genotypes cooperatively modify soybean rhizosphere Microbiome assembly. BMC Microbiol. 2019;19:201. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12866-019-1572-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Leite J, Fischer D, Rouws LFM, Fernandes-Júnior PI, Hofmann A, Kublik S, et al. Cowpea nodules harbor non-rhizobial bacterial communities that are shaped by soil type rather than plant genotype. Front Plant Sci. 2017;7:2064. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2016.02064.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Dastogeer KMG, Tumpa FH, Sultana A, Akter MA, Chakraborty A. Plant microbiome–an account of the factors that shape community composition and diversity. Curr Plant Biol. 2020;23:100161. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.cpb.2020.100161.

    Article  Google Scholar 

  63. de Vries FT, Manning P, Tallowin JRB, Mortimer SR, Pilgrim ES, Harrison KA, et al. Abiotic drivers and plant traits explain landscape-scale patterns in soil microbial communities. Ecol Lett. 2012;15:1230–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/J.1461-0248.2012.01844.X.

    Article  PubMed  Google Scholar 

  64. Dong M, Kowalchuk GA, Liu H, Xiong W, Deng X, Zhang N, et al. Microbial community assembly in soil aggregates: A dynamic interplay of stochastic and deterministic processes. Appl Soil Ecol. 2021;163:103911. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.apsoil.2021.103911.

    Article  Google Scholar 

  65. Cao H, Chen R, Wang L, Jiang L, Yang F, Zheng S, et al. Soil pH, total phosphorus, climate and distance are the major factors influencing microbial activity at a regional Spatial scale. Sci Rep. 2016;6:25815. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/srep25815.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Ortúzar M, Riesco R, Criado M, Alonso M, del Trujillo P. Unraveling the dynamic interplay of microbial communities associated to Lupinus angustifolius in response to environmental and cultivation conditions. Sci Total Environ. 2024;946:174277. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/J.SCITOTENV.2024.174277.

    Article  PubMed  Google Scholar 

  67. Zhang Y, Yang Y, Osborne B, Zhou H, Wu J, Zhang W, et al. Edaphic factors control microbial biomass and elemental stoichiometry in alpine meadow soils of the Tibet plateau. Plant Soil. 2024. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-024-06578-x.

    Article  Google Scholar 

  68. Shah ASS, Wakelin SAA, Moot DJJ, Blond C, Laugraud A, Ridgway HJJ. Trifolium repens and T. subterraneum modify their nodule Microbiome in response to soil pH. J Appl Microbiol. 2021;131:1858–69. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/jam.15050.

    Article  CAS  PubMed  Google Scholar 

  69. Singh U, Choudhary AK, Sharma S. A 3-year field study reveals that agri-management practices drive the dynamics of dominant bacterial taxa in the rhizosphere of Cajanus cajan. Symbiosis. 2022;86:215–27. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s13199-022-00834-3.

    Article  Google Scholar 

  70. Vineela C, Wani SP, Srinivasarao C, Padmaja B, Vittal KPR. Microbial properties of soils as affected by cropping and nutrient management practices in several long-term manurial experiments in the semi-arid tropics of India. Appl Soil Ecol. 2008;40:165–73. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/J.APSOIL.2008.04.001.

    Article  Google Scholar 

  71. Sharaf H, Thompson AA, Williams MA, Peck GM. Compost applications increase bacterial community diversity in the Apple rhizosphere. Soil Sci Soc Am J. 2021;85:1105–21. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/saj2.20251.

    Article  CAS  Google Scholar 

  72. Liu B, Ahnemann H, Arlotti D, Huyghebaert B, Cuperus F, Tebbe CC. Impact of diversified cropping systems and fertilization strategies on soil microbial abundance and functional potentials for nitrogen cycling. Sci Total Environ. 2024;932:172954. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/J.SCITOTENV.2024.172954.

    Article  CAS  PubMed  Google Scholar 

  73. Su J-Y, Liu C-H, Tampus K, Lin Y-C, Huang C-H. Organic amendment types influence soil properties, the soil bacterial microbiome, and tomato growth. Agronomy. 2022;12:1236. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/agronomy12051236.

    Article  CAS  Google Scholar 

  74. Paudel D, Wang L, Poudel R, Acharya JP, Victores S, de Souza CHL, et al. Elucidating the effects of organic vs. conventional cropping practice and rhizobia inoculation on rhizosphere microbial diversity and yield of peanut. Environ Microbiome. 2023;18:60. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40793-023-00517-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Dang P, Lu C, Huang T, Zhang M, Yang N, Han X, et al. Enhancing intercropping sustainability: manipulating soybean rhizosphere Microbiome through cropping patterns. Sci Total Environ. 2024;931:172714. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/J.SCITOTENV.2024.172714.

    Article  CAS  PubMed  Google Scholar 

  76. Yu R-P, Yang H, Xing Y, Zhang W-P, Lambers H, Li L. Belowground processes and sustainability in agroecosystems with intercropping. Plant Soil 2022:1–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-022-05487-1

  77. Vora SM, Ankati S, Patole C, Podile AR, Archana G. Alterations of primary metabolites in root exudates of intercropped Cajanus cajan–Zea mays modulate the adaptation and proteome of Ensifer (Sinorhizobium) fredii NGR234. Microb Ecol. 2022;83:1008–25. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00248-021-01818-4

  78. Vora SM, Joshi P, Belwalkar M, Archana G. Root exudates influence chemotaxis and colonization of diverse plant growth promoting rhizobacteria in the pigeon pea– maize intercropping system. Rhizosphere. 2021;18:100331. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.rhisph.2021.100331.

    Article  Google Scholar 

  79. Liu Y, Yin X, Xiao J, Tang L, Zheng Y. Interactive influences of intercropping by nitrogen on flavonoid exudation and nodulation in faba bean. Sci Rep. 2019;9:4818. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-019-41146-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Bell-Dereske LP, Benucci GMN, da Costa PB, Bonito G, Friesen ML, Tiemann LK, et al. Regional biogeography versus intra-annual dynamics of the root and soil Microbiome. Environ Microbiome. 2023;18:50. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40793-023-00504-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Van Horn DJ, Van Horn ML, Barrett JE, Gooseff MN, Altrichter AE, Geyer KM, et al. Factors controlling soil microbial biomass and bacterial diversity and community composition in a cold desert ecosystem: role of geographic scale. PLoS ONE. 2013;8:e66103. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pone.0066103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Kumar S, Suyal DC, Yadav A, Shouche Y, Goel R. Microbial diversity and soil physiochemical characteristic of higher altitude. PLoS ONE. 2019;14:e0213844. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/JOURNAL.PONE.0213844.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Shiro S, Matsuura S, Saiki R, Sigua GC, Yamamoto A, Umehara Y, et al. Genetic diversity and geographical distribution of Indigenous soybean-nodulating bradyrhizobia in the united States. Appl Environ Microbiol. 2013;79:3610–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/AEM.00236-13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Zhang B, Du N, Li Y, Shi P, Wei G. Distinct biogeographic patterns of rhizobia and non-rhizobial endophytes associated with soybean nodules across China. Sci Total Environ. 2018;643:569–78. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2018.06.240.

    Article  CAS  PubMed  Google Scholar 

  85. Miranda-Sánchez F, Rivera J, Vinuesa P. Diversity patterns of Rhizobiaceae communities inhabiting soils, root surfaces and nodules reveal a strong selection of rhizobial partners by legumes. Environ Microbiol. 2016;18:2375–91. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/1462-2920.13061.

    Article  CAS  PubMed  Google Scholar 

  86. Sarao SK, Boothe V, Das BK, Gonzalez-Hernandez JL, Brözel VS. Bradyrhizobium and the soybean rhizosphere: species level bacterial population dynamics in established soybean fields, rhizosphere and nodules. Plant Soil 2024:1–16. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-024-06814-4

  87. Hirsch PR, Mauchline TH. Who’s who in the plant root microbiome? Nat Biotechnol. 2012;30:961–2. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nbt.2387.

    Article  CAS  PubMed  Google Scholar 

  88. Reinhold-Hurek B, Bünger W, Burbano CS, Sabale M, Hurek T. Roots shaping their microbiome: global hotspots for microbial activity. Annu Rev Phytopathol. 2015;53:403–24. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev-phyto-082712-102342.

    Article  CAS  PubMed  Google Scholar 

  89. Edwards J, Johnson C, Santos-Medellín C, Lurie E, Podishetty NK, Bhatnagar S, et al. Structure, variation, and assembly of the root-associated microbiomes of rice. Proc Natl Acad Sci. 2015;112:E911–20. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.1414592112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Bulgarelli D, Rott M, Schlaeppi K, Ver L, van Themaat E, Ahmadinejad N, Assenza F, et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature. 2012;488:91–5. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nature11336.

    Article  CAS  PubMed  Google Scholar 

  91. Zgadzaj R, Garrido-Oter R, Jensen DB, Koprivova A, Schulze-Lefert P, Radutoiu S. Root nodule symbiosis in Lotus japonicus drives the establishment of distinctive rhizosphere, root, and nodule bacterial communities. Proceedings of the National Academy of Sciences. 2016;113:E7996–8005. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.1616564113

  92. Garrido-Oter R, Nakano RT, Dombrowski N, Ma K-W, McHardy AC, Schulze-Lefert P. Modular traits of the rhizobiales root microbiota and their evolutionary relationship with symbiotic rhizobia. Cell Host Microbe. 2018;24:155–67. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.chom.2018.06.006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Pini F, Galardini M, Bazzicalupo M, Mengoni A. Plant-bacteria association and symbiosis: are there common genomic traits in alphaproteobacteria?? Genes (Basel). 2011;2:1017–32. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes2041017.

    Article  CAS  PubMed  Google Scholar 

  94. Tak N, Awasthi E, Bissa G, Meghwal RR, James EK, Sprent JS, et al. Multi locus sequence analysis and symbiotic characterization of novel Ensifer strains nodulating Tephrosia spp. Indian Thar Desert Syst Appl Microbiol. 2016;39:534–45. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.syapm.2016.08.002.

    Article  CAS  PubMed  Google Scholar 

  95. Rathi S, Tak N, Bissa G, Chouhan B, Ojha A, Adhikari D, et al. Selection of Bradyrhizobium or Ensifer symbionts by the native Indian caesalpinioid legume Chamaecrista pumila depends on soil pH and other edaphic and Climatic factors. FEMS Microbiol Ecol. 2018;94:fiy180. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/femsec/fiy180.

    Article  CAS  Google Scholar 

  96. Choudhary S, Tak N, Bissa G, Chouhan B, Choudhary P, Sprent JI, et al. The widely distributed legume tree Vachellia (Acacia) Nilotica subsp. Indica is nodulated by genetically diverse Ensifer strains in India. Symbiosis. 2020;80:15–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s13199-019-00658-8.

  97. Rascovan N, Carbonetto B, Perrig D, Díaz M, Canciani W, Abalo M, et al. Integrated analysis of root microbiomes of soybean and wheat from agricultural fields. Sci Rep. 2016;6:28084. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/srep28084.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Horner A, Browett SS, Antwis RE. Mixed-cropping between field pea varieties alters root bacterial and fungal communities. Sci Rep. 2019;9:16953. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-019-53342-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  99. Brown SP, Grillo MA, Podowski JC, Heath KD. Soil origin and plant genotype structure distinct Microbiome compartments in the model legume Medicago truncatula. Microbiome. 2020;8:139. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-020-00915-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Cordero J, de Freitas JR, Germida JJ. Bacterial Microbiome associated with the rhizosphere and root interior of crops in Saskatchewan, Canada. Can J Microbiol. 2020;66:71–85. https://doiorg.publicaciones.saludcastillayleon.es/10.1139/cjm-2019-0330.

    Article  CAS  PubMed  Google Scholar 

  101. Tkacz A, Bestion E, Bo Z, Hortala M, Poole PS. Influence of plant fraction, soil, and plant species on microbiota: a multikingdom comparison. mBio. 2020;11:e02785–19. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/mBio.02785-19.

    Article  PubMed  PubMed Central  Google Scholar 

  102. Peix A, Ramírez-Bahena MH, Velázquez E, Bedmar EJ. Bacterial associations with legumes. CRC Crit Rev Plant Sci. 2015;34:17–42. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/07352689.2014.897899.

    Article  Google Scholar 

  103. Saxena RK, Von Wettberg E, Upadhyaya HD, Sanchez V, Songok S, Saxena K, et al. Genetic diversity and demographic history of Cajanus spp. Illustrated from genome-wide SNPs. PLoS ONE. 2014;9:e88568. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pone.0088568.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  104. Pérez-Jaramillo JE, Mendes R, Raaijmakers JM. Impact of plant domestication on rhizosphere Microbiome assembly and functions. Plant Mol Biol. 2016;90:635–44. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11103-015-0337-7.

    Article  CAS  PubMed  Google Scholar 

  105. Pérez-Jaramillo JE, Carrión VJ, de Hollander M, Raaijmakers JM. The wild side of plant microbiomes. Microbiome. 2018;6:143. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-018-0519-z.

    Article  PubMed  PubMed Central  Google Scholar 

  106. Zhang YF, Wang ET, Tian CF, Wang FQ, Han LL, Chen WF, et al. Bradyrhizobium elkanii, Bradyrhizobium yuanmingense and Bradyrhizobium japonicum are the main rhizobia associated with Vigna unguiculata and Vigna radiata in the subtropical region of China. FEMS Microbiol Lett. 2008;285:146–54. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1574-6968.2008.01169.x.

    Article  CAS  PubMed  Google Scholar 

  107. Appunu C, N’Zoue A, Moulin L, Depret G, Laguerre G. Vigna mungo, V. radiata and V. unguiculata plants sampled in different agronomical–ecological–climatic regions of India are nodulated by Bradyrhizobium yuanmingense. Syst Appl Microbiol. 2009;32:460–70. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.syapm.2009.05.005.

    Article  CAS  PubMed  Google Scholar 

  108. Osei O, Abaidoo RC, Ahiabor BDK, Boddey RM, Rouws LFM. Bacteria related to Bradyrhizobium yuanmingense from Ghana are effective groundnut micro-symbionts. Appl Soil Ecol. 2018;127:41–50. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.apsoil.2018.03.003.

    Article  PubMed  PubMed Central  Google Scholar 

  109. Araujo J, Díaz-Alcántara C-A, Urbano B, González-Andrés F. Inoculation with native Bradyrhizobium strains formulated with Biochar as carrier improves the performance of Pigeonpea (Cajanus Cajan L). Eur J Agron. 2020;113:125985. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.eja.2019.125985.

    Article  CAS  Google Scholar 

  110. Alegria Terrazas R, Balbirnie-Cumming K, Morris J, Hedley PE, Russell J, Paterson E, et al. A footprint of plant eco-geographic adaptation on the composition of the barley rhizosphere bacterial microbiota. Sci Rep. 2020;10:12916. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-020-69672-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  111. Cordovez V, Dini-Andreote F, Carrión VJ, Raaijmakers JM. Ecology and evolution of plant microbiomes. Annu Rev Microbiol. 2019;73:69–88. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev-micro-090817-062524.

    Article  CAS  PubMed  Google Scholar 

  112. Cardoso P, Alves A, Silveira P, Sá C, Fidalgo C, Freitas R, et al. Bacteria from nodules of wild legume species: phylogenetic diversity, plant growth promotion abilities and osmotolerance. Sci Total Environ. 2018;645:1094–102. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2018.06.399.

    Article  CAS  PubMed  Google Scholar 

  113. Liu B, Zhang D, Pan X. Nodules of wild legumes as unique natural hotspots of antibiotic resistance genes. Sci Total Environ. 2022;839:156036. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/J.SCITOTENV.2022.156036.

    Article  CAS  PubMed  Google Scholar 

  114. Chang C, Chen W, Luo S, Ma L, Li X, Tian C. Rhizosphere microbiota assemblage associated with wild and cultivated soybeans grown in three types of soil suspensions. Arch Agron Soil Sci. 2019;65:74–87. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/03650340.2018.1485147.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We acknowledge Department of Science and Technology (DST), Government of India (GoI), for a JC Bose Fellowship to ARP (JCB/2017/000053), University Grants Commission (UGC) for research fellowships to AB and DC, and University of Hyderabad for IoE-PDRFs fellowship to DC. We acknowledge DST Funds for Infrastructure in Science and Technology (DST-FIST), Level II, and UGC Special Assistance Program (UGC-SAP) support to the Department of Plant Sciences at the University of Hyderabad. We thank Dr. R. Dileep Kumar (SHUATS, Prayagraj, India) for his help with the sample collection.

Funding

This work was financially supported by the Department of Science and Technology (DST), Government of India (GoI), under project number EMR/2017/002658, and by the University of Hyderabad– Institution of Eminence project, under grant number UoH-IoE-RC3-21-065, granted to ARP. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

ARP conceived and planned the study; AB, DC, and SPVSRN designed the study, conducted the experiments, and analyzed the data; SU, VRC, and PLC conducted the experiments; AB, DC, and SPVSRN drafted the manuscript; ARP reviewed and revised the manuscript. All authors have read, critically revised, and approved the final version of the manuscript.

Corresponding author

Correspondence to Appa Rao Podile.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Clinical Trial Number

Not applicable.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Additional File 1

: Soil characterization data.

Additional File 2

: Detailed processing of sequencing data across the various experimental phases; Supplementary Figures S1–S61; Supplementary Tables S1–S4.

Additional File 3

: Bash code, merging report, and final reads per sample for all samples.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Basu, A., Chalasani, D., Sarma, P. et al. Influence of genotype, nodule position, and edaphic factors on microbial diversity and assembly of pigeonpea (Cajanus cajan) root nodules in Indian soils. Environmental Microbiome 20, 41 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40793-025-00707-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40793-025-00707-4

Keywords