Population structure of the clover rhizobia Rhizobium leguminosarum bv. trifolii upon transition from soil into the nodular niche

High-throughput sequencing of the amplicon gene library revealed variations in the population structure of clover rhizobia (Rhizobium leguminosarum bv. trifolii) upon transition from soil into the root nodules of the host plant (Trifolium hybridum). Analysis of rhizobial diversity using the nodA gene revealed 3258 and 1449 nucleotide sequences (allelic variants) for the soil and root nodule population, respectively. They were combined into 29 operational taxonomic units (OTU) according to the 97% identity level; 24 OTU were found in the soil population, 12 were present in the root nodule population, and 7 were common. The predominant OTE13 (77.4 and 91.5% of the soil and root nodule populations, respectively) contained 155 and 200 variants of the soil and root nodule populations, respectively, with the nucleotide diversity increasing significantly upon the “soil → root” transition. The “moving window” approach was used to reveal the sites of the nodA gene in which polymorphism, including that associated with increased frequency of non-synonymous substitution frequency, increased sharply upon transition from soil into root nodules. PCR analysis of the IGS genotypes of individual strains revealed insignificant changes in rhizobial diversity upon transition from soil into root nodules. These results indicate that acceleration of rhizobial evolution in the course of symbiosis may be associated with development of highly polymorphic virulent subpopulations subjected to directional selection in the “plant-soil” system.