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.

The Regulation of Pea (Pisum sativum L.) Symbiotic Nodule Infection and Defense Responses by Glutathione, Homoglutathione, and Their Ratio.

In this study, the roles of glutathione (GSH), homoglutathione (hGSH), and their ratio in symbiotic nodule development and functioning, as well as in defense responses accompanying ineffective nodulation in pea (Pisum sativum) were investigated. The expression of genes involved in (h)GSH biosynthesis, thiol content, and localization of the reduced form of GSH were analyzed in nodules of wild-type pea plants and mutants sym33-3 (weak allele, "locked" infection threads, occasional bacterial release, and defense reactions) and sym33-2 (strong allele, "locked" infection threads, defense reactions), and sym40-1 (abnormal bacteroids, oxidative stress, early senescence, and defense reactions). The effects of (h)GSH depletion and GSH treatment on nodule number and development were also examined. The GSH:hGSH ratio was found to be higher in nodules than in uninoculated roots in all genotypes analyzed, with the highest value being detected in wild-type nodules. Moreover, it was demonstrated, that a hGSHS-to-GSHS switch in gene expression in nodule tissue occurs only after bacterial release and leads to an increase in the GSH:hGSH ratio. Ineffective nodules showed variable GSH:hGSH ratios that correlated with the stage of nodule development. Changes in the levels of both thiols led to the activation of defense responses in nodules. The application of a (h)GSH biosynthesis inhibitor disrupted the nitrogen fixation zone in wild-type nodules, affected symbiosome formation in sym40-1 mutant nodules, and meristem functioning and infection thread growth in sym33-3 mutant nodules. An increase in the levels of both thiols following GSH treatment promoted both infection and extension of defense responses in sym33-3 nodules, whereas a similar increase in sym40-1 nodules led to the formation of infected cells resembling wild-type nitrogen-fixing cells and the disappearance of an early senescence zone in the base of the nodule. Meanwhile, an increase in hGSH levels in sym40-1 nodules resulting from GSH treatment manifested as a restriction of infection similar to that seen in untreated sym33-3 nodules. These findings indicated that a certain level of thiols is required for proper symbiotic nitrogen fixation and that changes in thiol content or the GSH:hGSH ratio are associated with different abnormalities and defense responses.

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.

Efficacy of a Plant-Microbe System: Pisum sativum (L.) Cadmium-Tolerant Mutant and Rhizobium leguminosarum Strains, Expressing Pea Metallothionein Genes PsMT1 and PsMT2, for Cadmium Phytoremediation.

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 SGECdt, which is characterized by increased tolerance to and accumulation of cadmium (Cd) in plants, was analyzed. Without Cd treatment, the wild type and mutant SGECdt 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 CdCl2 were characterized by an enlarged senescence zone. It was in stark contrast to Cd-treated nodules of the mutant SGECdt that maintained their proper organization. Cadmium treatment of either wild-type SGE or mutant SGECdt 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 SGECdt 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.

New sources of Sym2A allele in the pea (Pisum sativum L.) carry the unique variant of candidate LysM-RLK gene LykX.

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 Sym2 A 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 Sym2 E 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 Sym2 A-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 Sym2 A, as has been confirmed by an allelism test with Sym2 A 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.

Mapping-by-sequencing using NGS-based 3'-MACE-Seq reveals a new mutant allele of the essential nodulation gene Sym33 (IPD3) in pea (Pisum sativum L.).

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 F 2-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 F 2-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.

Mechanism of FACT removal from transcribed genes by anticancer drugs curaxins.

Human FACT (facilitates chromatin transcription) is a multifunctional protein complex that has histone chaperone activity and facilitates nucleosome survival and transcription through chromatin. Anticancer drugs curaxins induce FACT trapping on chromatin of cancer cells (c-trapping), but the mechanism of c-trapping is not fully understood. Here, we show that in cancer cells, FACT is highly enriched within the bodies of actively transcribed genes. Curaxin-dependent c-trapping results in redistribution of FACT from the transcribed chromatin regions to other genomic loci. Using a combination of biochemical and biophysical approaches, we have demonstrated that FACT is bound to and unfolds nucleosomes in the presence of curaxins. This tight binding to the nucleosome results in inhibition of FACT-dependent transcription in vitro in the presence of both curaxins and competitor chromatin, suggesting a mechanism of FACT trapping on bulk nucleosomes (n-trapping).

Pea Marker Database (PMD) - A new online database combining known pea (Pisum sativum L.) gene-based markers.

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.

Nucleosome-free DNA regions differentially affect distant communication in chromatin.

Communication between distantly spaced genomic regions is one of the key features of gene regulation in eukaryotes. Chromatin per se can stimulate efficient enhancer-promoter communication (EPC); however, the role of chromatin structure and dynamics in this process remains poorly understood. Here we show that nucleosome spacing and the presence of nucleosome-free DNA regions can modulate chromatin structure/dynamics and, in turn, affect the rate of EPC in vitro and in silico. Increasing the length of internucleosomal linker DNA from 25 to 60 bp results in more efficient EPC. The presence of longer nucleosome-free DNA regions can positively or negatively affect the rate of EPC, depending upon the length and location of the DNA region within the chromatin fiber. Thus the presence of histone-free DNA regions can differentially affect the efficiency of EPC, suggesting that gene regulation over a distance could be modulated by changes in the length of internucleosomal DNA spacers.

Overcoming a nucleosomal barrier to replication.

Efficient overcoming and accurate maintenance of chromatin structure and associated histone marks during DNA replication are essential for normal functioning of the daughter cells. However, the molecular mechanisms of replication through chromatin are unknown. We have studied traversal of uniquely positioned mononucleosomes by T7 replisome in vitro. Nucleosomes present a strong, sequence-dependent barrier for replication, with particularly strong pausing of DNA polymerase at the +(31-40) and +(41-65) regions of the nucleosomal DNA. The exonuclease activity of T7 DNA polymerase increases the overall rate of progression of the replisome through a nucleosome, likely by resolving nonproductive complexes. The presence of nucleosome-free DNA upstream of the replication fork facilitates the progression of DNA polymerase through the nucleosome. After replication, at least 50% of the nucleosomes assume an alternative conformation, maintaining their original positions on the DNA. Our data suggest a previously unpublished mechanism for nucleosome maintenance during replication, likely involving transient formation of an intranucleosomal DNA loop.

Large-scale ATP-independent nucleosome unfolding by a histone chaperone.

DNA accessibility to regulatory proteins is substantially influenced by nucleosome structure and dynamics. The facilitates chromatin transcription (FACT) complex increases the accessibility of nucleosomal DNA, but the mechanism and extent of its nucleosome reorganization activity are unknown. Here we determined the effects of FACT from the yeast Saccharomyces cerevisiae on single nucleosomes by using single-particle Förster resonance energy transfer (spFRET) microscopy. FACT binding results in dramatic ATP-independent, symmetrical and reversible DNA uncoiling that affects at least 70% of the DNA within a nucleosome, occurs without apparent loss of histones and proceeds via an 'all-or-none' mechanism. A mutated version of FACT is defective in uncoiling, and a histone mutation that suppresses phenotypes caused by this FACT mutation in vivo restores the uncoiling activity in vitro. Thus, FACT-dependent nucleosome unfolding modulates the accessibility of nucleosomal DNA, and this activity is an important function of FACT in vivo.

De Novo Assembly of the Pea (Pisum sativum L.) Nodule Transcriptome.

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.

Structure of transcribed chromatin is a sensor of DNA damage.

Early detection and repair of damaged DNA is essential for cell functioning and survival. Although multiple cellular systems are involved in the repair of single-strand DNA breaks (SSBs), it remains unknown how SSBs present in the nontemplate strand (NT-SSBs) of DNA organized in chromatin are detected. The effect of NT-SSBs on transcription through chromatin by RNA polymerase II was studied. NT-SSBs localized in the promoter-proximal region of nucleosomal DNA and hidden in the nucleosome structure can induce a nearly quantitative arrest of RNA polymerase downstream of the break, whereas more promoter-distal SSBs moderately facilitate transcription. The location of the arrest sites on nucleosomal DNA suggests that formation of small intranucleosomal DNA loops causes the arrest. This mechanism likely involves relief of unconstrained DNA supercoiling accumulated during transcription through chromatin by NT-SSBs. These data suggest the existence of a novel chromatin-specific mechanism that allows the detection of NT-SSBs by the transcribing enzyme.

Structural analysis of nucleosomal barrier to transcription.

Thousands of human and Drosophila genes are regulated at the level of transcript elongation and nucleosomes are likely targets for this regulation. However, the molecular mechanisms of formation of the nucleosomal barrier to transcribing RNA polymerase II (Pol II) and nucleosome survival during/after transcription remain unknown. Here we show that both DNA-histone interactions and Pol II backtracking contribute to formation of the barrier and that nucleosome survival during transcription likely occurs through allosterically stabilized histone-histone interactions. Structural analysis indicates that after Pol II encounters the barrier, the enzyme backtracks and nucleosomal DNA recoils on the octamer, locking Pol II in the arrested state. DNA is displaced from one of the H2A/H2B dimers that remains associated with the octamer. The data reveal the importance of intranucleosomal DNA-protein and protein-protein interactions during conformational changes in the nucleosome structure on transcription. Mechanisms of nucleosomal barrier formation and nucleosome survival during transcription are proposed.

Preparation and analysis of positioned mononucleosomes.

Short DNA fragments containing single nucleosomes have been extensively employed as simple model experimental systems for analysis of many intranuclear processes, including binding of proteins to nucleosomes, covalent histone modifications, transcription, DNA repair, and ATP-dependent chromatin remodeling. Here we describe several recently developed procedures for obtaining and analysis of mononucleosomes assembled on 200-350-bp DNA fragments.

Preparation of mononucleosomal templates for analysis of transcription with RNA polymerase using spFRET.

Single positioned nucleosomes have been extensively employed as simple model experimental systems for analysis of various intranuclear processes. Here we describe an experimental system containing positioned mononucleosomes allowing transcription by various RNA polymerases. Each DNA template contains a pair of fluorescent labels (Cy3 and Cy5) allowing measuring relative distances between the neighboring coils of nucleosomal DNA using Forster resonance energy transfer (FRET). The single-particle FRET (spFRET) approach for analysis of DNA uncoiling from the histone octamer during transcription through chromatin is described in detail.

Experimental analysis of hFACT action during Pol II transcription in vitro.

FACT (facilitates chromatin transcription) is a histone chaperone that facilitates transcription through chromatin and promotes histone recovery during transcription. Here, we describe a highly purified experimental system that recapitulates many important properties of transcribed chromatin and the key aspects of hFACT action during this process in vitro. We present the protocols describing how to prepare different forms of nucleosomes, including intact nucleosome, covalently conjugated nucleosome, nucleosome missing one of the two H2A/2B dimers (hexasome) and tetrasome (a nucleosome missing both H2A/2B dimers). These complexes allow analysis of various aspects of FACT's function. These approaches and other methods described below can also be applied to the study of other chromatin remodelers and chromatin-targeted factors.

Analysis of the mechanism of nucleosome survival during transcription.

Maintenance of nucleosomal structure in the cell nuclei is essential for cell viability, regulation of gene expression and normal aging. Our previous data identified a key intermediate (a small intranucleosomal DNA loop, Ø-loop) that is likely required for nucleosome survival during transcription by RNA polymerase II (Pol II) through chromatin, and suggested that strong nucleosomal pausing guarantees efficient nucleosome survival. To evaluate these predictions, we analysed transcription through a nucleosome by different, structurally related RNA polymerases and mutant yeast Pol II having different histone-interacting surfaces that presumably stabilize the Ø-loop. The height of the nucleosomal barrier to transcription and efficiency of nucleosome survival correlate with the net negative charges of the histone-interacting surfaces. Molecular modeling and analysis of Pol II-nucleosome intermediates by DNase I footprinting suggest that efficient Ø-loop formation and nucleosome survival are mediated by electrostatic interactions between the largest subunit of Pol II and core histones.