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.

The large size and complexity of the garden pea (Pisum sativum L.) genome hamper its sequencing and the discovery of pea gene resources. Although transcriptome sequencing provides extensive information about expressed genes, some tissue-specific transcripts can only be identified from particular organs under appropriate conditions. In this study, we performed RNA sequencing of polyadenylated transcripts from young pea nodules and root tips on an Illumina GAIIx system, followed by de novo transcriptome assembly using the Trinity program. We obtained more than 58,000 and 37,000 contigs from “Nodules” and “Root Tips” assemblies, respectively. The quality of the assemblies was assessed by comparison with pea expressed sequence tags and transcriptome sequencing project data available from NCBI website. The “Nodules” assembly was compared with the “Root Tips” assembly and with pea transcriptome sequencing data from projects indicating tissue specificity. As a result, approximately 13,000 nodule-specific contigs were found and annotated by alignment to known plant protein-coding sequences and by Gene Ontology searching. Of these, 581 sequences were found to possess full CDSs and could thus be considered as novel nodule-specific transcripts of pea. The information about pea nodule-specific gene sequences can be applied for gene-based markers creation, polymorphism studies, and real-time PCR.

At the onset of legume-rhizobial symbiosis, the mutual recognition of partners occurs based on a complicated interaction between signal molecules and receptors. Bacterial signal molecules named Nod factors (“nodulation factors”) are perceived by the plant LysM-containing receptor-like kinases (LysM-RLKs) that recognize details of its structure (i.e., unique substitutions), thus providing the conditions particular to symbiosis. In the garden pea (Pisum sativum L.), the allelic state of Sym2 gene has long been reported to regulate the symbiotic specificity: for infection to be successful, plants with the Sym2A allele (for “Sym2 Afghan”, as these genotypes originate mostly from Afghanistan) require an additional acetylation of the Nod factor which is irrelevant for genotypes with the Sym2E allele (for “Sym2 European”). Despite being described about 90 years ago, Sym2 has not yet been cloned, though phenotypic analysis suggests it probably encodes a receptor for the Nod factor. Recently, we described a novel pea gene LykX (PsLykX) from the LysM-RLK gene family that demonstrates a perfect correlation between its allelic state and the symbiotic specificity of the Sym2A-type. Here we report on a series of Middle-Eastern pea genotypes exhibiting the phenotype of narrow symbiotic specificity discovered in the VIR plant genetic resources gene bank (Saint-Petersburg, Russia). These genotypes are new sources of Sym2A, as has been confirmed by an allelism test with Sym2A pea cv. Afghanistan. Within these genotypes, LykX is present either in the allelic state characteristic for cv. Afghanistan, or in another, minor allelic state found in two genotypes from Tajikistan and Turkmenistan. Plants carrying the second allele demonstrate the same block of rhizobial infection as cv. Afghanistan when inoculated with an incompatible strain. Intriguingly, this “Tajik” allele of LykX differs from the “European” one by a single nucleotide polymorphism leading to an R75P change in the receptor part of the putative protein. Thus, our new data are in agreement with the hypothesis concerning the identity of LykX and the elusive Sym2 gene.

Two transgenic strains of Rhizobium leguminosarum bv. viciae, 3841-PsMT1 and 3841-PsMT2, were obtained. These strains contain the genetic constructions nifH-PsMT1 and nifH-PsMT2 coding for two pea (Pisum sativum L.) metallothionein genes, PsMT1 and PsMT2, fused with the promoter region of the nifH gene. The ability of both transgenic strains to form nodules on roots of the pea wild-type SGE and the mutant SGECd t , which is characterized by increased tolerance to and accumulation of cadmium (Cd) in plants, was analyzed. Without Cd treatment, the wild type and mutant SGECd t inoculated with R. leguminosarum strains 3841, 3841-PsMT1, or 3841-PsMT2 were similar histologically and in their ultrastructural organization of nodules. Nodules of wild-type SGE inoculated with strain 3841 and exposed to 0.5 μM CdCl 2 were characterized by an enlarged senescence zone. It was in stark contrast to Cd-treated nodules of the mutant SGECd t that maintained their proper organization. Cadmium treatment of either wild-type SGE or mutant SGECd t did not cause significant alterations in histological organization of nodules formed by strains 3841-PsMT1 and 3841-PsMT2. Although some abnormalities were observed at the ultrastructural level, they were less pronounced in the nodules of strain 3841-PsMT1 than in those formed by 3841-PsMT2. Both transgenic strains also differed in their effects on pea plant growth and the Cd and nutrient contents in shoots. In our opinion, combination of Cd-tolerant mutant SGECd t and the strains 3841-PsMT1 or 3841-PsMT2 may be used as an original model for study of Cd tolerance mechanisms in legume-rhizobial symbiosis and possibilities for its application in phytoremediation or phytostabilization technologies.

Large collections of pea symbiotic mutants were accumulated in the 1990s, but the causal genes for a large portion of the mutations are still not identified due to the complexity of the task. We applied a Mapping-by-Sequencing approach including Bulk Segregant Analysis and Massive Analysis of cDNA Ends (MACE-Seq) sequencing technology for genetic mapping the Sym11 gene of pea which controls the formation of symbioses with both nodule bacteria and arbuscular-mycorrhizal fungi. For mapping we developed an F2-population from the cross between pea line N24 carrying the mutant allele of sym11 and the wild type NGB1238 (=JI0073) line. Sequencing libraries were prepared from bulks of 20 plants with mutant and 12 with wild-type phenotype. MACE-Seq differential gene expression analysis between mutant-phenotype and wild-type-phenotype bulks revealed 2,235 genes, of which 514 (23%) were up-regulated and 1,721 (77%) were down-regulated in plant roots inoculated with rhizobia as a consequence of sym11 mutation. MACE-Seq also detected single nucleotide variants between bulks in 217 pea genes. Using a novel mathematical model we calculated the recombination frequency (RF) between the Sym11 gene and these 217 polymorphic genes. Six genes with the lowest RF were converted into CAPS or dCAPS markers and genetically mapped on the complete mapping population of 108 F2-plants which confirmed their tight linkage to Sym11 and to each other. The Medicago truncatula Gaertn. (Mt) homologs of these genes are located in a distinct region of Mt chromosome 5, which corresponds to linkage group I of pea. Among 94 candidate genes from this region only one was down-regulated—the pea Sym33 homolog of the Mt IPD3 gene which is essential for nodulation. Sequencing of the Sym33 allele of the N24 (sym11) mutant revealed a single nucleotide deletion (c.C319del) in its third exon resulting in a codon shift in the open reading frame and premature translation termination. Thus, we identified a novel mutant allele sym33-4 most probably responsible for the mutant phenotype of the N24 (sym11) line, thereby demonstrating that mapping by MACE-Seq can be successfully used for genetic mapping of mutations and identification of candidate genes in pea.