Genetic diversity of indigenous rhizobial symbionts of the Lupinus mariae-josephae endemism from alkaline-limed soils within its area of distribution in Eastern Spain

The genomic diversity of a collection of 103 indigenous rhizobia isolates from Lupinus mariae-josephae (Lmj), a recently described Lupinus species endemic to alkaline-limed soils from a restricted habitat in Eastern Spain, was investigated by molecular methods. Isolates were obtained from soils of four geographic locations in the Valencia province that harbored the known Lmj plant populations. Using an M13 RAPD fingerprinting technique, 19 distinct RAPD profiles were identified. Phylogenetic analysis based on 16S rDNA and the housekeeping genes glnII, recA and atpD showed a high diversity of native Bradyrhizobium strains that were able to establish symbiosis with Lmj. All the strains grouped in a clade unrelated to strains of the B. canariense and B. japonicum lineages that establish symbioses with lupines in acid soils of the Mediterranean area. The phylogenetic tree based on concatenated glnII, recA and atpD gene sequences grouped the Lmj isolates in six different operational taxonomic units (OTUs) at the 93% similarity level. These OTUs were not associated to any specific geographical location, and their observed divergence predicted the existence of different Bradyrhizobium genomic species. In contrast, phylogenetic analysis of symbiotic genes based on nodC and nodA gene sequences, defined only two distinct clusters among the Lmj strains. These two Lmj nod gene types were largely distinct from nod genes of bradyrhizobia nodulating other Old World lupine species. The singularity and large diversity of these strains in such a small geographical area makes this an attractive system for studying the evolution and adaptation of the rhizobial symbiont to the plant host.