Molecular systematics of rhizobia based on maximum likelihood and Bayesian phylogenies inferred from rrs, atpD, recA and nifH sequences, and their use in the classification of Sesbania microsymbionts from Venezuelan wetlands

A well-resolved rhizobial species phylogeny with 51 haplotypes was inferred from a combined atpD+ recA data set using Bayesian inference with best-fit, gene-specific substitution models. Relatively dense taxon sampling for the genera Rhizobium and Mesorhizobium was achieved by generating atpD and recA sequences for six type and 24 reference strains not previously available in GenBank. This phylogeny was used to classify nine nodule isolates from Sesbania exasperata, S. punicea and S. sericea plants native to seasonally flooded areas of Venezuela, and compared with a PCR-RFLP analysis of rrs plus rrl genes and large maximum likelihood rrs and nifH phylogenies. We show that rrs phylogenies are particularly sensitive to strain choice due to the high levels of sequence mosaicism found at this locus. All analyses consistently identified the Sesbania isolates as Mesorhizobium plurifarium or Rhizobium huautlense. Host range experiments on ten legume species coupled with plasmid profiling uncovered potential novel biovarieties of both species. This study demonstrates the wide geographic and environmental distribution of M. plurifarium, that R. galegae and R. huautlense are sister lineages, and the synonymy of R. gallicum, R. mongolense and R. yanglingense. Complex and diverse phylogeographic, inheritance and host-association patterns were found for the symbiotic nifH locus. The results and the analytical approaches used herein are discussed in the context of rhizobial taxonomy and molecular systematics.