Genomic and Transcriptomic Analysis of Pea (Pisum sativum L.) Breeding Line 'Triumph' with High Symbiotic Responsivity.

Pea (Pisum sativum L.), like most legumes, forms mutualistic symbioses with nodule bacteria and arbuscular mycorrhizal (AM) fungi. The positive effect of inoculation is partially determined by the plant genotype; thus, pea varieties with high and low symbiotic responsivity have been described, but the molecular genetic basis of this trait remains unknown. Here, we compare the symbiotically responsive breeding line 'Triumph' of grain pea with its parental cultivars 'Vendevil' (a donor of high symbiotic responsivity) and 'Classic' (a donor of agriculturally valuable traits) using genome and transcriptome sequencing. We show that 'Triumph' inherited one-fourth of its genome from 'Vendevil', including the genes related to AM and nodule formation, and reveal that under combined inoculation with nodule bacteria and AM fungi, 'Triumph' and 'Vendevil', in contrast to 'Classic', demonstrate similar up-regulation of the genes related to solute transport, hormonal regulation and flavonoid biosynthesis in their roots. We also identify the gene PsGLP2, whose expression pattern distinguishing 'Triumph' and 'Vendevil' from 'Classic' correlates with difference within the promoter region sequence, making it a promising marker for the symbiotic responsivity trait. The results of this study may be helpful for future molecular breeding programs aimed at creation of symbiotically responsive cultivars of pea.

Identification and Expression Analysis of the C-TERMINALLY ENCODED PEPTIDE Family in Pisum sativum L.

The C-TERMINALLY ENCODED PEPTIDE(CEP) peptides play crucial roles in plant growth and response to environmental factors. These peptides were characterized as positive regulators of symbiotic nodule development in legume plants. However, little is known about the CEP peptide family in pea. Here, we discovered in pea genome 21 CEP genes (PsCEPs), among which three genes contained additional conserved motifs corresponding to the PIP (PAMP-induced secreted peptides) consensus sequences. We characterized the expression patterns of pea PsCEP genes based on transcriptomic data, and for six PsCEP genes with high expression levels in the root and symbiotic nodules the detailed expression analysis at different stages of symbiosis and in response to nitrate treatment was performed. We suggest that at least three PsCEP genes, PsCEP1, PsCEP7 and PsCEP2, could play a role in symbiotic nodule development, whereas the PsCEP1 and PsCEP13 genes, downregulated by nitrate addition, could be involved in regulation of nitrate-dependent processes in pea. Further functional studies are required to elucidate the functions of these PsCEP genes.

A variable gene family encoding nodule-specific cysteine-rich peptides in pea (Pisum sativum L.).

Various legume plants form root nodules in which symbiotic bacteria (rhizobia) fix atmospheric nitrogen after differentiation into a symbiotic form named bacteroids. In some legume species, bacteroid differentiation is promoted by defensin-like nodule-specific cysteine-rich (NCR) peptides. NCR peptides have best been studied in the model legume Medicago truncatula Gaertn., while in many other legumes relevant information is still fragmentary. Here, we characterize the NCR gene family in pea (Pisum sativum L.) using genomic and transcriptomic data. We found 360 genes encoding NCR peptides that are expressed in nodules. The sequences of pea NCR genes and putative peptides are highly variable and differ significantly from NCR sequences of M. truncatula. Indeed, only one pair of orthologs (PsNCR47-MtNCR312) has been identified. The NCR genes in the pea genome are located in clusters, and the expression patterns of NCR genes from one cluster tend to be similar. These data support the idea of independent evolution of NCR genes by duplication and diversification in related legume species. We also described spatiotemporal expression profiles of NCRs and identified specific transcription factor (TF) binding sites in promoters of "early" and "late" NCR genes. Further, we studied the expression of NCR genes in nodules of Fix- mutants and predicted potential regulators of NCR gene expression, one among them being the TF ERN1 involved in the early steps of nodule organogenesis. In general, this study contributes to understanding the functions of NCRs in legume nodules and contributes to understanding the diversity and potential antibiotic properties of pea nodule-specific antimicrobial molecules.

DNA Methylation Patterns Differ between Free-Living Rhizobium leguminosarum RCAM1026 and Bacteroids Formed in Symbiosis with Pea (Pisum sativum L.).

Rhizobium leguminosarum (Rl) is a common name for several genospecies of rhizobia able to form nitrogen-fixing nodules on the roots of pea (Pisum sativum L.) while undergoing terminal differentiation into a symbiotic form called bacteroids. In this work, we used Oxford Nanopore sequencing to analyze the genome methylation states of the free-living and differentiated forms of the Rl strain RCAM1026. The complete genome was assembled; no significant genome rearrangements between the cell forms were observed, but the relative abundances of replicons were different. GANTC, GGCGCC, and GATC methylated motifs were found in the genome, along with genes encoding methyltransferases with matching predicted target motifs. The GGCGCC motif was completely methylated in both states, with two restriction-modification clusters on different replicons enforcing this specific pattern of methylation. Methylation patterns for the GANTC and GATC motifs differed significantly depending on the cell state, which indicates their possible connection to the regulation of symbiotic differentiation. Further investigation into the differences of methylation patterns in the bacterial genomes coupled with gene expression analysis is needed to elucidate the function of bacterial epigenetic regulation in nitrogen-fixing symbiosis.

Transposons-Based Clonal Diversity in Trematode Involves Parts of CR1 (LINE) in Eu- and Heterochromatin.

Trematode parthenitae have long been believed to form clonal populations, but clonal diversity has been discovered in this asexual stage of the lifecycle. Clonal polymorphism in the model species Himasthla elongata has been previously described, but the source of this phenomenon remains unknown. In this work, we traced cercarial clonal diversity using a simplified amplified fragment length polymorphism (SAFLP) method and characterised the nature of fragments in diverse electrophoretic bands. The repetitive elements were identified in both the primary sequence of the H. elongata genome and in the transcriptome data. Long-interspersed nuclear elements (LINEs) and long terminal repeat retrotransposons (LTRs) were found to represent an overwhelming majority of the genome and the transposon transcripts. Most sequenced fragments from SAFLP pattern contained the reverse transcriptase (RT, ORF2) domains of LINEs, and only a few sequences belonged to ORFs of LTRs and ORF1 of LINEs. A fragment corresponding to a CR1-like (LINE) spacer region was discovered and named CR1-renegade (CR1-rng). In addition to RT-containing CR1 transcripts, we found short CR1-rng transcripts in the redia transcriptome and short contigs in the mobilome. Probes against CR1-RT and CR1-rng presented strikingly different pictures in FISH mapping, despite both being fragments of CR1. In silico data and Southern blotting indicated that CR1-rng is not tandemly organised. CR1 involvement in clonal diversity is discussed.

Transcriptome Analysis of Alternative Splicing Events Induced by Arbuscular Mycorrhizal Fungi (Rhizophagus irregularis) in Pea (Pisum sativum L.) Roots.

Alternative splicing (AS), a process that enables formation of different mRNA isoforms due to alternative ways of pre-mRNA processing, is one of the mechanisms for fine-tuning gene expression. Currently, the role of AS in symbioses formed by plants with soil microorganisms is not fully understood. In this work, a comprehensive analysis of the transcriptome of garden pea (Pisum sativum L.) roots in symbiosis with arbuscular mycorrhiza was performed using RNAseq and following bioinformatic analysis. AS profiles of mycorrhizal and control roots were highly similar, intron retention accounting for a large proportion of the observed AS types (67%). Using three different tools (SUPPA2, DRIMSeq and IsoformSwitchAnalyzeR), eight genes with AS events specific for mycorrhizal roots of pea were identified, among which four were annotated as encoding an apoptosis inhibitor protein, a serine/threonine-protein kinase, a dehydrodolichyl diphosphate synthase, and a pre-mRNA-splicing factor ATP-dependent RNA helicase DEAH1. In pea mycorrhizal roots, the isoforms of these four genes with preliminary stop codons leading to a truncated ORFs were up-regulated. Interestingly, two of these four genes demonstrating mycorrhiza-specific AS are related to the process of splicing, thus forming parts of the feedback loops involved in fine-tuning of gene expression during mycorrhization.