Bacteria currently included in Rhizobium leguminosarum are too diverse to be considered a single species, so we can refer to this as a species complex (the Rlc). We have found 429 publicly available genome sequences that fall within the Rlc and these show that the Rlc is a distinct entity, well separated from other species in the genus. Its sister taxon is R. anhuiense. We constructed a phylogeny based on concatenated sequences of 120 universal (core) genes, and calculated pairwise average nucleotide identity (ANI) between all genomes. From these analyses, we concluded that the Rlc includes 18 distinct genospecies, plus 7 unique strains that are not placed in these genospecies. Each genospecies is separated by a distinct gap in ANI values, usually at approximately 96% ANI, implying that it is a ‘natural’ unit. Five of the genospecies include the type strains of named species: R. laguerreae, R. sophorae, R. ruizarguesonis, “R. indicum” and R. leguminosarum itself. The 16S ribosomal RNA sequence is remarkably diverse within the Rlc, but does not distinguish the genospecies. Partial sequences of housekeeping genes, which have frequently been used to characterize isolate collections, can mostly be assigned unambiguously to a genospecies, but alleles within a genospecies do not always form a clade, so single genes are not a reliable guide to the true phylogeny of the strains. We conclude that access to a large number of genome sequences is a powerful tool for characterizing the diversity of bacteria, and that taxonomic conclusions should be based on all available genome sequences, not just those of type strains.

Arbuscular mycorrhiza (AM) is known to be a mutually beneficial plant-fungal symbiosis; however, the effect of mycorrhization is heavily dependent on multiple biotic and abiotic factors. Therefore, for the proper employment of such plant-fungal symbiotic systems in agriculture, a detailed understanding of the molecular basis of the plant developmental response to mycorrhization is needed. The aim of this work was to uncover the physiological and metabolic alterations in pea (Pisum sativumL.) leaves associated with mycorrhization at key plant developmental stages. Plants of pea cv. Finale were grown in constant environmental conditions under phosphate deficiency. The plants were analyzed at six distinct time points, which corresponded to certain developmental stages of the pea: I: 7 days post inoculation (DPI) when the second leaf is fully unfolded with one pair of leaflets and a simple tendril; II: 21 DPI at first leaf with two pairs of leaflets and a complex tendril; III: 32 DPI when the floral bud is enclosed; IV: 42 DPI at the first open flower; V: 56 DPI when the pod is filled with green seeds; and VI: 90–110 DPI at the dry harvest stage. Inoculation withRhizophagus irregularishad no effect on the fresh or dry shoot weight, the leaf photochemical activity, accumulation of chlorophylla,bor carotenoids. However, at stage III (corresponding to the most active phase of mycorrhiza development), the number of internodes between cotyledons and the youngest completely developed leaf was lower in the inoculated plants than in those without inoculation. Moreover, inoculation extended the vegetation period of the host plants, and resulted in increase of the average dry weight per seed at stage VI. The leaf metabolome, as analyzed with GC-MS, included about three hundred distinct metabolites and showed a strong correlation with plant age, and, to a lesser extent, was influenced by mycorrhization. Metabolic shifts influenced the levels of sugars, amino acids and other intermediates of nitrogen and phosphorus metabolism. The use of unsupervised dimension reduction methods showed that (i) at stage II, the metabolite spectra of inoculated plants were similar to those of the control, and (ii) at stages IV and V, the leaf metabolic profiles of inoculated plants shifted towards the profiles of the control plants at earlier developmental stages. At stage IV the inoculated plants exhibited a higher level of metabolism of nitrogen, organic acids, and lipophilic compounds in comparison to control plants. Thus, mycorrhization led to the retardation of plant development, which was also associated with higher seed biomass accumulation in plants with an extended vegetation period. The symbiotic crosstalk between host plant and AM fungi leads to alterations in several biochemical pathways the details of which need to be elucidated in further studies.

Legume crops represent the major source of food protein and contribute to human nutrition and animal feeding. An essential improvement of their productivity can be achieved by symbiosis with beneficial soil microorganisms—rhizobia (Rh) and arbuscular mycorrhizal (AM) fungi. The efficiency of these interactions depends on plant genotype. Recently, we have shown that, after simultaneous inoculation with Rh and AM, the productivity gain of pea (Pisum sativum L) line K-8274, characterized by high efficiency of interaction with soil microorganisms (EIBSM), was higher in comparison to a low-EIBSM line K-3358. However, the molecular mechanisms behind this effect are still uncharacterized. Therefore, here, we address the alterations in pea seed proteome, underlying the symbiosis-related productivity gain, and identify 111 differentially expressed proteins in the two lines. The high-EIBSM line K-8274 responded to inoculation by prolongation of seed maturation, manifested by up-regulation of proteins involved in cellular respiration, protein biosynthesis, and down-regulation of late-embryogenesis abundant (LEA) proteins. In contrast, the low-EIBSM line K-3358 demonstrated lower levels of the proteins, related to cell metabolism. Thus, we propose that the EIBSM trait is linked to prolongation of seed filling that needs to be taken into account in pulse crop breeding programs. The raw data have been deposited to the ProteomeXchange with identifier PXD013479.

During the initial step of the symbiosis between legumes (Fabaceae) and nitrogen-fixing bacteria (rhizobia), the bacterial signal molecule known as the Nod factor (nodulation factor) is recognized by plant LysM motif-containing receptor-like kinases (LysM-RLKs). The fifth chromosome of barrel medic (Medicago truncatula Gaertn.) contains a cluster of paralogous LysM-RLK genes, one of which is known to participate in symbiosis. In the syntenic region of the pea (Pisum sativum L.) genome, three genes have been identified: PsK1 and PsSym37, two symbiosis-related LysM-RLK genes with known sequences, and the unsequenced PsSym2 gene which presumably encodes a LysM-RLK and is associated with increased selectivity to certain Nod factors. In this work, we identified a new gene encoding a LysM-RLK, designated as PsLykX, within the Sym2 genomic region. We sequenced the first exons (corresponding to the protein receptor domain) of PsSym37, PsK1, and PsLykX from a large set of pea genotypes of diverse origin. The nucleotide diversity of these fragments was estimated and groups of haplotypes for each gene were revealed. Footprints of selection pressure were detected via comparative analyses of SNP distribution across the first exons of these genes and their homologs MtLYK2, MtLYK3, and MtLYK4 from M. truncatula retrieved from the Medicago Hapmap project. Despite the remarkable similarity among all the studied genes, they exhibited contrasting selection signatures, possibly pointing to diversification of their functions. Signatures of balancing selection were found in LysM1-encoding parts of PsSym37 and PsK1, suggesting that the diversity of these parts may be important for pea LysM-RLKs. The first exons of PsSym37 and PsK1 displayed signatures of purifying selection, as well as MtLYK2 of M. truncatula. Evidence of positive selection affecting primarily LysM domains was found in all three investigated M. truncatula genes, as well as in the pea gene PsLykX. The data suggested that PsLykX is a promising candidate for PsSym2, which has remained elusive for more than 30 years.

Rhizobium leguminosarum bv. viciae RCAM1026 is a strain first isolated in 1964 from nodules of "Ramensky 77" cultivar of garden pea (Pisum sativum L.) now routinely used as a model strain in inoculation experiments on pea. Assembly with SPAdes yielded 133 contigs longer then 200 bp (N50 = 202,321, GC% = 60.84). Resulting annotated genome is 7,248,686 bp encoding 6792 genes.

Pea (Pisum sativum L.) is the oldest model object of plant genetics and one of the most agriculturally important legumes in the world. Since the pea genome has not been sequenced yet, identification of genes responsible for mutant phenotypes or desirable agricultural traits is usually performed via genetic mapping followed by candidate gene search. Such mapping is best carried out using gene-based molecular markers, as it opens the possibility for exploiting genome synteny between pea and its close relative Medicago truncatula Gaertn., possessing sequenced and annotated genome. In the last 5 years, a large number of pea gene-based molecular markers have been designed and mapped owing to the rapid evolution of “next-generation sequencing” technologies. However, the access to the complete set of markers designed worldwide is limited because the data are not uniformed and therefore hard to use. The Pea Marker Database was designed to combine the information about pea markers in a form of user-friendly and practical online tool. Version 1 (PMD1) comprises information about 2484 genic markers, including their locations in linkage groups, the sequences of corresponding pea transcripts and the names of related genes in M. truncatula. Version 2 (PMD2) is an updated version comprising 15944 pea markers in the same format with several advanced features. To test the performance of the PMD, fine mapping of pea symbiotic genes Sym13 and Sym27 in linkage groups VII and V, respectively, was carried out. The results of mapping allowed us to propose the Sen1 gene (a homologue of SEN1 gene of Lotus japonicus (Regel) K. Larsen) as the best candidate gene for Sym13, and to narrow the list of possible candidate genes for Sym27 to ten, thus proving PMD to be useful for pea gene mapping and cloning. All information contained in PMD1 and PMD2 is available at www.peamarker.arriam.ru.