Tatiana A. Serova
Anna V. Tsyganova
Igor A. Tikhonovich
Viktor E. Tsyganov- 1Laboratory of Molecular and Cellular Biology, Department of Biotechnology, All-Russia Research Institute for Agricultural Microbiology, Russian Academy of Agricultural Sciences, Saint Petersburg, Russia
- 2Department of Genetics and Biotechnology, Saint Petersburg State University, Saint Petersburg, Russia
The development of nitrogen-fixing nodules formed during Rhizobium–legume symbiosis is strongly controlled by phytohormones. In this study, we investigated the effect of gibberellins (GAs) on senescence of pea (Pisum sativum) symbiotic nodules. Pea wild-type line SGE, as well as corresponding mutant lines SGEFix--1 (sym40), SGEFix--2 (sym33), SGEFix--3 (sym26), and SGEFix--7 (sym27), blocked at different stages of nodule development, were used in the study. An increase in expression of the GA2ox1 gene, encoding an enzyme involved in GA deactivation (GA 2-oxidase), and a decrease in the transcript abundance of the GA20ox1 gene, encoding one of the enzymes involved in GA biosynthesis (GA 20-oxidase), were observed in analyzed genotypes during nodule aging. A reduction in the amount of bioactive GA3 was demonstrated by immunolocalization in the early senescent mutant and wild-type lines during aging of symbiotic nodules. Down-regulated expression of senescence-associated genes encoding cysteine proteases 1 and 15a, thiol protease, bZIP transcription factor, 1-aminocyclopropane-1-carboxylate (ACC) synthase, ACC oxidase, and aldehyde oxidase was observed in the nodules of wild-type plants treated with exogenous GA3 relative to the untreated plants. GA3-treated plants also showed increases in nodule size and the nitrogen fixation zone, and decreases in the number of nodules and the senescence zone. Immunogold localization revealed higher levels of GA3 in the peribacteroid spaces in symbiosomes than in the matrix of infection threads. Furthermore, a decrease in GA3 label in mature and senescent symbiosomes in comparison with juvenile symbiosomes was observed. These results suggest a negative effect of GAs on the senescence of the pea symbiotic nodule and possible involvement of GAs in functioning of the mature nodule. Simultaneously, GA3 treatment led to nodule meristem bifurcation, indicating a possible role of GAs in nodule meristem functioning.
Introduction
Legume–Rhizobium interactions culminate in the formation of nitrogen-fixing nodules. Rhizobia growing near the plant adhere to the root hair and eventually penetrate the root, initiating the formation of an infection thread (Brewin, 2004; Tsyganova and Tsyganov, 2017). When the infection thread reaches reactivated cortical root cells, which form a nodule primordium, rhizobia are released into the host cell cytoplasm from unwalled infection droplets (Brewin, 2004). After release, rhizobia differentiate into bacteroids and become surrounded by a peribacteroid membrane; these form symbiosomes, organelle-like structures in which bacteroids fix nitrogen (Tsyganova et al., 2018).
If the nodule meristem functions for a long time, nodules of an indeterminate type are formed with different histological zones, which include the meristem (zone I), the infection zone (zone II), the nitrogen fixation zone (zone III), and the senescence zone (zone IV) (Guinel, 2009).
Senescence completes symbiotic nodule development and is accompanied by the destruction of symbiotic partners, large-scale protein degradation, and remobilization of nutrients to other plant organs (Puppo et al., 2005; Serova and Tsyganov, 2014). In particular, the catabolism of leghemoglobin, which is one of the most abundant proteins in the nodule, is observed during nodule senescence resulting in a color change of aged nodules from pink to green. Senescence in the indeterminate nodule is associated with the senescence zone formed at the base of the nodule and spreads toward its apical part and periphery (Pérez Guerra et al., 2010; Dupont et al., 2012). Hormonal regulation has a major impact on symbiotic nodule development (Ferguson and Mathesius, 2014; Tsyganova and Tsyganov, 2015, 2018). Current data suggest that both ethylene and abscisic acid (ABA) contribute to the aging of the symbiotic nodule (Puppo et al., 2005; Van de Velde et al., 2006; Karmarkar, 2014; Serova et al., 2017). In contrast, based on expression analysis of the nodules of Medicago truncatula (Van de Velde et al., 2006) and pea (Serova et al., 2017), it has been suggested that gibberellins (GAs) may have a negative impact on nodule senescence.
Gibberellins are a large group of diterpenoid carboxylic acids in higher plants. GAs stimulate organ growth, causing the enhancement of cell elongation and cell division (Hedden and Thomas, 2012). GA biosynthesis includes several steps catalyzed by terpene cyclases (Hedden and Thomas, 2012). The first steps involve the production of GA12, the common precursor of all types of GAs in plants (Hedden and Phillips, 2000). GA12 can be converted to another GA precursor, GA53. The final stages of GA biosynthesis are catalyzed by GA 20-oxidase and GA 3-oxidase. Their activity contributes to the content of bioactive forms of GA in the plant. In pea, GA 20-oxidases encoded by PsGA20ox genes (PsGA20ox1, PsGA20ox2) are involved in different stages of GA biosynthesis. They mainly catalyze the conversion of GA12 to GA9 and GA53 to GA20, where GA9 and GA20 are precursors of bioactive GAs (García-Martínez et al., 1997). The conversion of GA9 and GA20 to bioactive forms of GAs is catalyzed by GA 3-oxidases (García-Martínez et al., 1997; Hedden and Thomas, 2012).
GA1, GA3, GA4, GA5, and GA7 are the most common biologically active forms in higher plants (Hayashi et al., 2014). Along with bioactive forms, plants also contain inactive forms of GAs, including the precursors and metabolites of active GAs. Inactive GAs are present at higher concentrations and may perform yet unknown functions (Hedden and Thomas, 2012; Hayashi et al., 2014). There are several mechanisms of GA inactivation, the most prevalent of which is 2β-hydroxylation, catalyzed by GA 2-oxidases (GA2oxs) (Thomas et al., 1999; Hedden and Thomas, 2012). In pea, conversion of bioactive GAs, GA1 and GA4, and of its precursors GA9 and GA20, to inactive catabolites occurred by C19-GA 2-oxidases encoded by the PsGA2ox1 and PsGA2ox2 genes (Lester et al., 1999; Martin et al., 1999; Hedden and Thomas, 2012). Inactivation of the precursors GA12 and GA53 is catalyzed by C20-GA 2-oxidases (Hedden and Thomas, 2012).
Optimal GA levels differ during various stages of plant development and are maintained through feed-back and feed-forward regulation of GA metabolism (Weston et al., 2008; Hedden and Thomas, 2012). Bioactive GAs reduce GA biosynthesis and enhance GA deactivation (Weston et al., 2008). A GA signal transduction pathway is triggered by the binding of GAs to the soluble receptor GID1 (GIBBERELLIN-INSENSITIVE DWARF 1) (Ueguchi-Tanaka et al., 2005). Downstream signal transduction pathways involve DELLA proteins, which are key repressors of GA responses (Davière and Achard, 2013). After GA binding, the formation of a GA-GID1-DELLA complex occurs with its subsequent degradation by the 26S proteasome in the nucleus (Sun, 2011).
During nodulation, GAs are involved in the negative control of rhizobial infection and the positive regulation of nodule development (Lievens et al., 2005; Maekawa et al., 2009; Ferguson et al., 2011; Hayashi et al., 2014; McAdam et al., 2018). Up-regulation of GA biosynthetic genes, GA20ox and GA3ox, was observed in the early stages of nodulation in Sesbania rostrata (Lievens et al., 2005), Glycine max (Hayashi et al., 2012), and Lotus japonicus (Kouchi et al., 2004). The negative effect of GAs on rhizobial infection in L. japonicus and S. rostrata was accompanied by suppressed expression of the transcription factors NIN (Nodule Inception) and NSP (Nodulation Signaling Pathway) (Kouchi et al., 2004; Lievens et al., 2005; Maekawa et al., 2009). Analysis of the M. truncatula root transcriptome revealed that Nod factor perception led to spatial-temporal activation of genes involved in GA biosynthesis and catabolism (Larrainzar et al., 2015). Within the first hours of infection, an increase in the expression of GA deactivation genes was shown, which may promote the Nod factor signaling pathway (Larrainzar et al., 2015). However, up-regulation of GA biosynthesis genes was later observed, which likely causes an increase in GA content and may limit further infection (Larrainzar et al., 2015). The negative impact of GAs on rhizobial infection is mediated through destruction of DELLA proteins. It was previously shown that DELLA proteins interact with and activate a wide set of transcription factors involved in Nod factor signaling [NSP1, NSP2, IPD3 (interacting protein with DMI3), NF-YA1 (nuclear factor-YA1), and ERN1 (ERF required for nodulation 1)] (Fonouni-Farde et al., 2016; Jin et al., 2016). Furthermore, the positive effect of DELLA proteins was confirmed by poor nodulation of della mutants in M. truncatula (Jin et al., 2016).
There is evidence to suggest that different optimal levels of bioactive GAs are required at different stages of nodulation (Ferguson et al., 2005, 2011). Moreover, optimal levels seem to be species-specific and to depend on growth conditions. Inhibition of GA biosynthesis reduced nodulation in S. rostrata (Lievens et al., 2005). In pea, reduced GA levels inhibited nodulation, as demonstrated using a series of GA biosynthesis mutants, including the severely GA-deficient na-1 pea mutant (Ferguson et al., 2005). In L. japonicus, treatment of wild-type plants with exogenous GA3 induced the initiation of divisions in the pericycle leading to formation of pseudo-nodules (Kawaguchi et al., 1996). However, in the snf1 and snf2 mutants of L. japonicus, application of exogenous GA3 suppressed spontaneous nodulation by inhibiting the NSP (Maekawa et al., 2009). High concentrations of exogenous GA3 inhibited root hair infection in S. rostrata (10-5 M) (Lievens et al., 2005) and L. japonicus (10-6, 10-7 M) (Maekawa et al., 2009), and also reduced nodulation in pea (10-3 M) (Ferguson et al., 2005) and M. truncatula (10-7–10-4 M) (Fonouni-Farde et al., 2016; Jin et al., 2016). In contrast, low concentrations (10-9 M) of exogenous GA3 promoted nodulation in pea plants (Ferguson et al., 2005). Further overexpression of a GA signaling component, SLEEPY1, resulted in fewer nodules in transgenic roots of L. japonicus than in wild-type roots (Maekawa et al., 2009). However, the GA1-overproducing sln mutant of pea, which is strongly blocked in GA deactivation, formed nodules only on lateral roots and had the same number of nodules as the wild-type (Ross et al., 1995; Lester et al., 1999; Ferguson et al., 2005, 2011).
An ambivalent role for GAs in nodulation was clarified in a recent study (McAdam et al., 2018). The authors confirmed the well-known negative effect of GAs on development of infection threads and demonstrated the positive effect on nodule tissue development and functioning of nitrogen-fixing nodules. The action of GAs was shown to be mediated by DELLA proteins, which, on the contrary, promote infection thread growth and inhibit the initiation of cortex cell division and nodule development (McAdam et al., 2018).
Little is known about the interactions between GAs and other hormones during nodulation. In an earlier study using the strongly GA-deficient na-1 pea mutant, it was shown that low GA content resulted in elevated ethylene levels and decreased nodulation (Ferguson et al., 2011). This was confirmed by partial recovery of nodulation in the mutant na-1 treated with an inhibitor of ethylene biosynthesis (Ferguson et al., 2011), while an increase in the level of GAs rescued the number and structure of mutant nodules (Ferguson et al., 2005). Recently, a detailed phenotypic characterization of the na-1 ein2 double mutant and an analysis of epistatic interactions of the mutations revealed that, to some extent, GAs inhibit infection thread formation independently of ethylene, and that GAs facilitate the formation of nodules by partial inhibition of ethylene (McAdam et al., 2018).
Thus, a number of previous studies confirm involvement of GAs in the initial stages of Rhizobium–legume symbiosis. However, at present, the effect of GAs on the later stages of nodule development is not well understood. Therefore, the aim of this study was to identify a role for GAs in senescence of the symbiotic nodule of pea. It is known that in wild-type pea nodules, senescence is initiated at 4 weeks after inoculation (WAI) and actively developed in 6-week-old nodules (Kardailsky and Brewin, 1996; Serova et al., 2017, 2018). The study was carried out using pea wild-type line SGE and a number of mutant lines, SGEFix--1 (sym40), SGEFix--2 (sym33), SGEFix--3 (sym26), and SGEFix--7 (sym27), which are blocked at different stages of nodule development. Recently, we demonstrated the activation of nodule senescence in all analyzed mutants, which starts in 2-week-old nodules and is especially active in 4-week-old nodules (Serova et al., 2018). Based on data obtained from expression, immunolocalization, and pharmacological analyses, we suggest a negative effect of GAs on the senescence of the pea symbiotic nodule and possible involvement of GAs in the functioning of nodule meristem and cells from the nitrogen fixation zone.
Materials and Methods
Plant Material, Bacterial Strain, and Plant Growth Conditions
The pea (Pisum sativum L.) laboratory line SGE (Kosterin and Rozov, 1993) and corresponding mutant lines SGEFix--1 (sym40), SGEFix--2 (sym33), SGEFix--3 (sym26), and SGEFix--7 (sym27) were used in this study (Table 1). Rhizobium leguminosarum bv. viciae strain 3841 (Wang et al., 1982) was used as an inoculant. Methods for sterilization of seeds, plant inoculation, and growth conditions were described previously (Serova et al., 2017). For microscopic and expression analyses and quantitative measurement of nodulation, nodules were collected at 2, 4, and 6 WAI.
Table 1. Plant material used in the study.
Expression Analysis of Genes Associated With GA Metabolism in SGE and Mutant Lines
Primers for one GA biosynthesis gene, PsGA20ox1, and one GA deactivation gene, PsGA2ox1, were designed previously (Weston et al., 2008; Serova et al., 2017). For RNA extraction, 2, 4, and 6 WAI nodules of SGE, mutants SGEFix--1 (sym40), SGEFix--2 (sym33), SGEFix--3 (sym26), and SGEFix--7 (sym27) were ground in liquid nitrogen. Total RNA from each sample was isolated using PureZol reagent (Bio-Rad, Hercules, CA, United States) according to the manufacturer’s recommendations. RNA quantity and purity were checked using a MultiNA electrophoresis system on microchips for analysis of nucleic acids (Shimadzu Corporation, Kyoto, Japan). Reverse transcription was performed on 1.5 μg total RNA treated with DNAse I (MBI Fermentas, Vilnius, Lithuania) using 200 U RevertAid Reverse Transcriptase and 0.5 μg Oligo(dT)18 (MBI Fermentas) for cDNA synthesis under manufacturer-recommended conditions. cDNA synthesis was carried out in 20 μl of reaction mix, and the resulting cDNAs were diluted five times for following use. The reaction was carried out in an automated C1000TM Thermal Cycler (Bio-Rad).
For gene expression quantification, relative real-time PCR was performed in a C1000TM Thermal Cycler combined with the optical module CFX96TM Real-Time System (Bio-Rad), using iQ SYBR Green Supermix (Bio-Rad) according to manufacturer’s instructions. Results of reactions were processed using Bio-Rad CFX Manager software (Bio-Rad). Relative expression was calculated with the 2-ΔΔCT method using the reference gene PsGapC1 (accession number L07500.1). Two-WAI nodules of SGE were used as a calibrator for calculation of relative transcript abundance. For PsGA20ox1 transcript abundance, 6-WAI nodules were used as a calibrator. Reactions were carried out in three technical replicates and averaged. Statistical treatment of experimental results was carried out with Microsoft Excel software. Statistically significant differences were calculated using one-way ANOVA at P-value ≤ 0.05. Experiments were performed in three replicates with six to eight plants per variant.
GA3 Immunolabeling and Confocal Microscopy
Nodules of SGE and of the mutants SGEFix--1 (sym40), SGEFix--2 (sym33), SGEFix--3 (sym26), and SGEFix--7 (sym27) at 2, 4, and 6 WAI were fixed in freshly prepared 4% paraformaldehyde buffered in PBS (136 mM NaCl, 2.68 mM KCl, 10 mM Na2HPO4, 1.7 mM KH2PO4, pH 7.4) with addition of 3% N-ethyl-N′-(3-dimethylaminopropyl) carbodiimide hydrochloride and 0.1% Triton X-100 (Sigma-Aldrich, Dorset, United Kingdom). Samples were fixed under vacuum (-0.9 bar), for 7 min three times with 15 min intervals using a VacuuBrand ME 1C vacuum pump (Vacuubrand, Wertheim, Germany), and incubated overnight at 4°C. Nodules were rinsed with PBS three times, with 15 min intervals, and stained with 0.5% toluidine blue solution in PBS for 1 h. Washing of the residual dye was carried out with PBS two times, with 15 min intervals. Samples were subsequently molded in 3% agarose gel.
Sections (50 μm) were prepared at room temperature with a HM650V microtome (Microm, Walldorf, Germany). After this, GA3 immunolabeling and nuclei and bacteria staining were carried out in accordance with Serova et al. (2018). Anti-GA3 rat antibodies (Agrisera, Vännäs, Sweden) and goat anti-rat IgG Alexa Fluor 488 (Thermo Fisher Scientific, Waltham, MA, United States) were used as primary and secondary antibodies for GA3 immunolabeling, respectively. Sections were mounted in ProLong Gold antifade reagent (Thermo Fisher Scientific).
Control of GA3-specific signal was carried out using GA3-BSA conjugate (Agrisera) (Supplementary Figure S1). Anti-GA3 antibodies were incubated with GA3-BSA conjugate in a 1:40 ratio in PBS for 24 h at 4°C in a total volume of 100 μl. Then, the mix was used as primary antibodies for immunolabeling of pea nodules (Supplementary Figures S1A–C). Anti-GA3 antibodies were omitted as a control for specific binding of secondary antibodies in the absence of primary antibodies (Supplementary Figures S1D–F). Also, the control of specificity of GA3-antibodies in nuclei is shown (Supplementary Figures S1G–L).
Sections were analyzed using the laser scanning confocal system LSM 510 META (Carl Zeiss, Oberkochen, Germany) and ZEN2009 software (Carl Zeiss).
Immunogold Labeling, Transmission Electron Microscopy, and Quantitative Analysis
Nodules of SGE and mutant SGEFix--3 (sym26) at 2 WAI were used. Preparation of samples for transmission electron microscopy (TEM) and immunogold labeling of GA3 was done with ultrathin sections on gold grids as described by Tsyganova et al. (2009). A total of 15–20 nodules from at least five different plants were fixed in 2.5% glutaraldehyde in 0.06 M phosphate buffer (pH 7.2) at 4°C overnight. Samples were then rinsed in buffer and dehydrated in increasing concentrations of ethanol (30% at room temperature; 50%, 70%, 90%, and twice with 100% at 35°C) for 20 min at each step. Subsequently, specimens were gradually infiltrated with increasing concentrations of LR-White resin (Polysciences Europe, Eppelheim, Germany) at a ratio of 1:1, 1:2, and 1:3 mixed with ethanol (100%) at -20°C and finally embedded in LR-White resin and polymerized at -20°C for 48 h in small plastic containers using UV polymerization in a Leica EM AFS2 (Leica Microsystems, Wetzlar, Germany).
Ultrathin sections (90 nm) of the samples were obtained with a Leica EM UC7 ultramicrotome (Leica Microsystems) and blocked with a blocking solution (5% BSA, 0.5% goat serum, 0.05% cold water fish skin) and they were then washed in 0.1% acetylated BSA (BSA-C) in PBS. Sections were treated with the primary anti-GA3 rat antibodies (Agrisera) diluted 1:25 in 0.1% BSA-C in PBS at 4°C overnight. After four rinses in 0.1% BSA-C in PBS, samples were incubated with a 10 nm gold-conjugated secondary antibody goat anti-rat IgG (Amersham International, Little Chalfont, United Kingdom), diluted 1:50 in 0.1% BSA-C in PBS, for 4 h at 37°C. After short washes in PBS and distilled water, labeled grids were post-stained with uranyl acetate for 15 s.
The specificity of the immunogold labeling procedures was tested by several negative controls. Negative controls were treated either with pre-immune serum instead of the primary antibody and with non-specific secondary antibody (goat anti-mouse IgG). Negative controls for gibberellic acid revealed that no labeling occurred on the section when they were treated with pre-immune serum instead of the primary antibody (Supplementary Figure S2A) and with non-specific secondary antibody (Supplementary Figure S2B).
Nodule tissues were analyzed in a JEM–1400 EM transmission electron microscope (JEOL Ltd., Tokyo, Japan) at 80 kV. Electron micrographs were obtained by Veleta CCD camera (Olympus, Münster, Germany). Micrographs of randomly photographed immunogold-labeled sections were digitized, and gold particles were counted in visually identified cell structures. For statistical analysis, at least 10 different samples of root nodules and at least 50 sectioned symbiosomes were examined for SGE or mutant SGEFix--3 (sym26). The SGEFix--3 (sym26) mutant was chosen as an example of mutants with premature degradation of symbiotic structures (early senescence phenotype) for comparison with wild-type nodules. Previously, we demonstrated that early senescence is more pronounced in this mutant (Serova et al., 2018). Morphometrical data were obtained as described previously (Ivanova et al., 2015). Briefly, the number of gold particles per unit area was calculated. The areas and the number of gold particles were measured using software Zen 2 Core version 2.5 (Carl Zeiss). The data are presented as the number of gold particles/μm2. Data were analyzed by one-way ANOVA using the software SigmaPlot for Windows version 12.5 (Systat Software, Inc., San Jose, CA, United States). Means were separated by the Tukey multiple range test (P-value ≤ 0.001).
Pharmacological Treatment
GA3 treatment was carried out for evaluation of GA action on pea nodule senescence. After inoculation with rhizobia, 100 ml of a 10-6 M GA3 (Sigma-Aldrich) aqueous solution was applied into substrate to SGE plants every 3 days until plants were harvested. Control plants were watered without GA3. Nodules were harvested for expression and light microscopy analyses. Photographs of pea nodules were made with a SteREO Lumar.V12 stereomicroscope equipped with an Axiocam ICc 1 video camera (Carl Zeiss).
Quantitative Measurements of Nodulation
GA3-treated and untreated pea plants were harvested for quantitative measurement. The shoots and roots were separated, and cotyledons were removed. Nodules on the primary and secondary roots were removed with a blade and counted. For weight measurements, separated shoots, roots, and nodules were dried at 42°C for 3 days. Experiments were performed with 9–16 plants per GA3-treated and untreated variants. Statistically significant differences were calculated using one-way ANOVA at P-value ≤ 0.01.
The sizes of nodules and senescence zones were measured with AxioVision Rel. 4.8 software (Carl Zeiss). Pictures of separate nodules of GA3-treated and untreated plants were taken as described above. Supplementary Figure S3 illustrates the selection and measurement of nodule projection areas using the example of analysis for a nodule of untreated plant at 6 WAI. Projection areas were measured as average projection area of nodules from the main root of a single plant. Experiments were performed with 8–12 plants per GA3-treated and untreated variants. Statistically significant differences were calculated using one-way ANOVA at P-value ≤ 0.01.
Light Microscopy Analysis of Nodules of GA3-Treated and Untreated Plants
Fixation of pea nodules for light microscopy was the same as described above. Samples were subsequently dehydrated with a series of ethanol solutions in water and embedded in Steedman’s wax as previously described (Serova et al., 2017). Sections of 10 μm were obtained with a HM360 microtome (Microm) and placed on slides in a few drops of water. After drying at 28°C for 30 min, slices were de-waxed in ethanol solutions and placed in PBS (Serova et al., 2017). Slices were then stained with toluidine blue (0.1% solution in PBS) for 10 min and washed in PBS (two times for 10 min). Sections were mounted in PBS.
Light microscopy analysis of nodule sections of GA3-treated and untreated plants was carried out with Axio Imager.Z1 (Carl Zeiss). Photographs were taken with a microscope camera Axiocam 506 color (Carl Zeiss) and analyzed using ZEN 2 core SP1 software (Carl Zeiss).
Expression Analysis of Senescence-Associated Genes
Primers for selected senescence-associated marker genes encoding cysteine protease 1 and 15a (PsCyp1, PsCyp15a), thiol protease (PsTPP), bZIP transcription factor (PsATB2), enzymes of ethylene (PsACS2, PsACO1) and ABA (PsAO3) biosynthesis, and an enzyme of bioactive GAs deactivation (PsGA2ox1) were designed previously (Serova et al., 2017). In addition, primers for a GA biosynthesis gene (PsGA20ox1) designed by Weston et al. (2008) were used.
RNA was extracted from 2, 4, and 6 WAI nodules of GA3-treated and untreated SGE as described above (see section “Expression Analysis of Genes Associated With GA Metabolism in SGE and Mutant Lines”). Reverse transcription and relative real-time PCR were performed in accordance with the previous description (see section “Expression Analysis of Genes Associated With GA Metabolism in SGE and Mutant Lines”).
Results
Expression Analysis of Genes Encoding Enzymes Involved in Biosynthesis and Deactivation of GAs
To study the involvement of GAs in senescence of pea symbiotic nodules, an expression analysis of one GA biosynthesis gene (PsGA20ox1) and one GA deactivation gene (PsGA2ox1) was carried out in SGE and the mutant lines (Figure 1).
Figure 1. Relative transcript levels of PsGA20ox1 (A) and PsGA2ox1 (B) genes in nodules of wild-type SGE and corresponding pea mutants SGEFix--1 (sym40), SGEFix--2 (sym33), SGEFix--3 (sym26), and SGEFix--7 (sym27) at 2, 4, and 6 weeks after inoculation (WAI). Letters indicate significant differences (one-way ANOVA, P-value ≤ 0.05) a, from wild-type SGE at one time point; b, within genotype compared with 2 WAI; c, within genotype compared with 4 WAI.
The transcript level of PsGA20ox1 was down-regulated due to the aging of the wild-type and the mutant nodules with the exception of the mutant SGEFix--2 (sym33), which showed statistically insignificant differences (Figure 1A). In the wild-type nodules, a significant (6.8-fold) decrease in transcript abundance was observed from 2 to 6 WAI. In 2- and 4-week-old nodules of the mutant SGEFix--2 (sym33), the PsGA20ox1 mRNA level was significantly reduced in comparison with the wild-type. In the mutant SGEFix--3 (sym26), a slight down-regulation was detected from 4 to 6 WAI only. Furthermore, in 4-week-old nodules of the mutant SGEFix--3 (sym26), the PsGA20ox1 mRNA level was 3.8-fold higher than that in the wild-type. A significant (7.6-fold) down-regulation of the transcript level was observed in the mutant SGEFix--7 (sym27) already at 4 WAI. In addition, PsGA20ox1 transcript abundance was 2.6-fold lower in 4-week-old nodules of the mutant SGEFix--7 (sym27) than in those of the wild-type. The difference in transcript levels in the nodules of the mutant SGEFix--1 (sym40) from the wild-type was insignificant.
During aging of wild-type nodules, PsGA2ox1 transcript abundance was significantly (48.4-fold) up-regulated at 6 WAI only. In contrast, the mutants SGEFix--3 (sym26) and SGEFix--7 (sym27) showed 23.5- and 49.4-fold elevation of the expression level, respectively, already at 4 WAI. However, in the case of the mutant SGEFix--7 (sym27), a decrease in PsGA2ox1 transcript abundance was detected at 6 WAI. In 4-week-old nodules of the mutant SGEFix--1 (sym40), an increase in PsGA2ox1 mRNA was less pronounced. PsGA2ox1 transcript abundance was significantly higher in the mutants SGEFix--1 (sym40), SGEFix--3 (sym26), and SGEFix--7 (sym27) than in the wild-type nodules at 4 WAI. The expression level of PsGA2ox1 gene was slightly higher than in the wild-type in 4-week-old nodules of the mutant SGEFix--2 (sym33). However, changes in the transcript level of the mutant SGEFix--2 (sym33) were statistically insignificant.
Overall, the expression of the GA biosynthesis gene PsGA20ox1 was decreased, while the expression of the GA deactivation gene PsGA2ox1 was increased during aging of the pea nodules of wild-type and early senescent mutants.
GA3 Immunolocalization in Wild-Type and Mutant Nodules
To complement the data obtained by expression analysis of GA metabolism genes, GA3 immunolocalization was carried out in wild-type and mutant nodules of different ages (Figures 2–4 and Supplementary Figures S4–S6).
Figure 2. Immunolocalization of gibberellin (GA3) in nodules of wild-type SGE (A–C) and corresponding pea mutants SGEFix--1 (sym40) (D–F), SGEFix--2 (sym33) (G–I), SGEFix--3 (sym26) (J–L), and SGEFix--7 (sym27) (M–O) at 2 weeks after inoculation. Zones of nodules are designated by Roman numerals: I – meristem, II – infection zone, III – fixation zone, III’ – zone corresponding to nitrogen fixation zone in wild-type. A differential interference contrast microscopy image merged with laser scanning confocal microscopy image in red channel (A,D,G,J,M). Merged images of laser scanning confocal microscopy in green and red channels (B,E,H,K,N). Heat map provides a color code of fluorescence signal intensities (C,F,I,L,O). Visualization of GA3 by the Alexa Fluor 488 conjugated secondary antibody (green), nuclei, and bacteria stained with propidium iodide (red). Scale bar = 100 μm.
Figure 3. Immunolocalization of gibberellin (GA3) in nodules of wild-type SGE (A–C) and corresponding pea mutants SGEFix--1 (sym40) (D–F), SGEFix--2 (sym33) (G–I), SGEFix--3 (sym26) (J–L), and SGEFix--7 (sym27) (M–O) at 4 weeks after inoculation. Zones of nodule are designated by Roman numerals: I – meristem, II – infection zone, III – fixation zone, III’ – zone corresponding to nitrogen fixation zone in wild-type, IV – senescence zone. A differential interference contrast microscopy image merged with laser scanning confocal microscopy image in red channel (A,D,G,J,M). Merged images of laser scanning confocal microscopy in green and red channels (B,E,H,K,N). Heat map provides a color code of fluorescence signal intensities (C,F,I,L,O). Visualization of GA3 by the Alexa Fluor 488 conjugated secondary antibody (green), nuclei, and bacteria stained with propidium iodide (red). Scale bar = 100 μm.
Figure 4. Immunolocalization of gibberellin (GA3) in cells in central part of nodules of wild-type SGE (A–C) and corresponding pea mutants SGEFix--1 (sym40) (D–F), SGEFix--2 (sym33) (G–I), SGEFix--3 (sym26) (J–L), and SGEFix--7 (sym27) (M–O) at 2 weeks after inoculation. ic, infected cell; uic, uninfected cell; n, nucleus. Arrow indicates infection thread, arrowhead indicates infection droplet. A differential interference contrast microscopy image merged with laser scanning confocal microscopy image in red channel (A,D,G,J,M). Merged images of laser scanning confocal microscopy in green and red channels (B,E,H,K,N). Heat map provides color code of fluorescence signal intensities (C,F,I,L,O). Visualization of GA3 by the Alexa Fluor 488 conjugated secondary antibody (green), nuclei, and bacteria stained with propidium iodide (red). Scale bar = 10 μm.
In 2- and 4-week-old wild-type nodules, GA3 labeling was detected in cells from the meristem, and in the infection and nitrogen fixation zones (Figures 2A–C, 3A–C and Supplementary Figure S4). The GA3 label was observed in both infected and uninfected cells; in the latter, it was associated with cytoplasm surrounding starch granules (Figures 4A–C). In wild-type nodules, a high signal of GA3 labeling was also detected in nuclei, mostly in cells from the nitrogen fixation zone (Figures 4A–C and Supplementary Figures S1G–L). The intensity of labeling was much lower in 6-week-old wild-type nodules (Supplementary Figures S5A–C,G–I), especially in the senescence zone (Supplementary Figures S5D–F). Traces of GA3 label were detected in senescent uninfected cells (Supplementary Figures S5J–L). In 2-week-old nodules of the mutant SGEFix--1 (sym40), a high intensity of labeling was observed in nuclei and cytoplasm, but was absent in infection threads and infection droplets (Figures 4D–F). Whole nodule intensity of labeling was lower in the mutant SGEFix--1 (sym40) than in the wild-type (Figures 2D–F). In 4-week-old nodules of the mutant SGEFix--1 (sym40), the intensity of GA3 labeling was significantly decreased (Figures 3D–F); traces of the GA3 label were observed in nuclei and cytoplasm, particularly at the cell periphery (Supplementary Figures S6A–C). In nodules of the mutant SGEFix--2 (sym33), the lowest level of GA3 was observed (Figures 2G–I). The GA3 label was detected mainly in nuclei and cytoplasm, but was absent in infection threads (Figures 4G–I). In 4-week-old mutant nodules, a high level of fluorescence was associated with nuclei (Figures 3G–I and Supplementary Figures S6D–F). A high intensity of fluorescence was detected in the mutants SGEFix--3 (sym26) (Figures 2J–L) and SGEFix--7 (sym27) (Figures 2M–O) in 2-week-old nodules. In SGEFix--3 (sym26), it was even higher than that of the wild-type nodules of the same age. In the mutants SGEFix--3 (sym26) and SGEFix--7 (sym27), the maximum GA3 signal was observed in infected and uninfected cells of regions corresponding to the nitrogen fixation zone of wild-type nodules (Figures 4J–O). The amount of GA3 was significantly decreased in the senescence zone occupying the dominant part of 4-week-old mutant nodules (Figures 3J–O). Trace levels of GA3 label were observed in uninfected cells (Supplementary Figures S6G–L).
Thus, a decline in the amount of detectable GA3 during aging of pea symbiotic nodules of wild-type and mutant lines and its distribution in different histological nodule zones were demonstrated.
Immunogold Localization of GA3 in Nodules of Wild-Type SGE and Mutant SGEFix--3 (sym26)
To compare the localization of GA3 in the infection structures in 2-week-old nodules, immunogold localization of GA3 was carried out in the wild-type and SGEFix--3 (sym26) nodules (Figure 5).
Figure 5. Immunogold localization of gibberellin (GA3) in nodules of wild-type SGE (A–D) and mutant SGEFix--3 (sym26) (E–H) at 2 weeks after inoculation. Secondary goat anti-rat IgG MAb conjugated to 10 nm diameter colloidal gold was used. it, infection thread; itw, infection thread wall; b, bacterium; ba, bacteroid; jba, juvenile bacteroid; dba, degrading bacteroid; asterisk indicates vesicle with GA3 label. (A,E) Infection threads, (B,F) juvenile symbiosomes, (C,G) mature symbiosomes, (D,H) senescent symbiosomes. Scale bar (A,E) = 500 nm, (B–D,F–H) = 200 nm.
Gold particles were observed in nuclei, vacuoles, and plastids, and were very abundant in cytoplasm (data not shown). Rare gold particles were observed in the matrix of infection threads in wild-type (Figure 5A) and SGEFix--3 (sym26) (Figure 5E) as well as in the infection droplets (Table 2). Gold particles were mainly found in the infection thread walls and in bacteria embedded in the infection thread matrix (Figures 5A,E). With respect to symbiosomes, numerous gold particles were observed mainly in bacteroids of both SGE and mutant SGEFix--3 (sym26) (Figures 5B–D,F–H), although some labeling was also found in the symbiosome membrane especially in the juvenile symbiosomes in nodules of wild-type SGE (Figure 5B). In nodules of SGE, some vesicles with GA3 label were found in the peribacteroid spaces (Figure 5C). The highest amount of gold particles was observed in juvenile symbiosomes in both analyzed genotypes (Figures 5B,F and Table 2). In mature symbiosomes and especially in senescent symbiosomes in nodules of SGE and mutant SGEFix--3 (sym26), the amount of gold particles of the GA3 label was lower than that in juvenile symbiosomes (Figures 5C,D,G,H and Table 2). It is necessary to note that the amount of gold particles of the GA3 label was lower in symbiosomes of the mutant SGEFix--3 (sym26) than in wild-type nodules (Table 2).
Table 2. Distribution of gold particles in 2-week-old nodules of pea wild-type SGE and mutant SGEFix--3 (sym26) (immunogold localization).
Thus, it was shown that with an increase in the age of wild-type and mutant nodules, the amount of GA3 decreased in symbiosomes.
Analysis of Nodules of SGE Treated With Exogenous GA3 Relative to Untreated Plants
The involvement of bioactive GAs in senescence of pea symbiotic nodules was also measured via GA3 treatment of wild-type plants (Figure 6). A quantitative analysis of 2-, 4-, and 6-week-old nodules of GA3-treated and untreated plants was carried out (Tables 3, 4). Additionally, an analysis of the effect of GA3-treatment on the histological structure of the symbiotic nodule was performed (Figure 7). To identify the involvement of bioactive GAs in senescence of pea nodules at the transcriptional level, the analysis of senescence-associated genes was carried out in nodules of GA3-treated and untreated plants at 2, 4, and 6 WAI (Figures 8, 9).
Figure 6. Nodulated main and lateral roots of wild-type SGE pea plants untreated (A,C,E) and treated with exogenous gibberellin (GA3) (B,D,F) at 2, 4, and 6 weeks after inoculation. Arrow indicates main roots; arrowheads indicate lateral roots. Scale bar = 2 mm.
Table 3. Shoot, root, and nodule dry weight (DW), nodule number per plant, percentage of nodules with meristem bifurcation and nodules without senescence zone of GA3-treated and untreated wild-type SGE at 2, 4, and 6 weeks after inoculation (WAI).
Table 4. Projection areas of whole nodule and senescence zone of GA3-treated and untreated wild-type SGE pea plants at 2, 4, and 6 weeks after inoculation (WAI).
Figure 7. Histological organization of nodules in wild-type SGE pea plants untreated (A) and treated with exogenous gibberellin (GA3) (B) at 6 weeks after inoculation. Zones in nodules are designated by Roman numerals: II – infection zone, III – fixation zone, IV – senescence zone. Sections were stained with toluidine blue. Scale bar = 500 μm.
Figure 8. Relative expression of (A–D) PsCyp1, PsCyp15a, PsTPP, and PsATB2 genes in nodules of GA3-treated and untreated wild-type SGE at 2, 4, and 6 weeks after inoculation (WAI). Letters indicate significant differences (one-way ANOVA, P-value ≤ 0.05, n = 3) a, from untreated plants at one time point; b, within variant (GA3-treated or untreated) compared with 2 WAI; c, within variant (GA3-treated or untreated) compared with 4 WAI.
Figure 9. Relative expression of (A–E) PsACS2, PsACO1, PsGA20ox1, PsGA2ox1, and PsAO3 genes in nodules of GA3-treated and untreated wild-type SGE at 2, 4, and 6 weeks after inoculation (WAI). Letters indicate significant differences (one-way ANOVA, P-value ≤ 0.05, n = 3) a, from untreated plants at one time point; b, within variant (GA3-treated or untreated) compared with 2 WAI; c, within variant (GA3-treated or untreated) compared with 4 WAI. Breaks of histogram bars indicate a change in scale.
Quantitative Measurements of Nodulation of GA3-Treated and Untreated Plants
During nodule aging, an increase in dry weight (DW) and size of nodules was observed in both GA3-treated and untreated plants (Tables 3, 4). In the untreated plants, the average nodule DW and the projection area of the whole nodule were increased 2.8 and 4 times, respectively, during aging from 2 to 6 WAI (Figures 6A,C,E and Tables 3, 4). In 6-week-old untreated plants, 10- and 3.7-fold increases in the total nodule DW and the number of nodules, respectively, were observed relative to 2-week-old plants (Table 3). A low degree (9.13 and 3.63%) of meristem bifurcation was detected in 4- and 6-week-old nodules of untreated plants, respectively (Figures 6C,E and Table 3). At 6 WAI, the senescence zone dominated in the nodules of untreated plants and occupied approximately 62% of the entire nodule (Figure 6E and Table 4).
In the case of GA3-treatment, 3- and 3.6-fold increases in the average nodule DW and projection area of the whole nodule, respectively, were observed at 6 WAI, in comparison with 2 WAI (Figures 6B,D,F and Tables 3, 4). In GA3-treated plants, the total nodule DW and the number of nodules were 5.4- and 2.7-times higher, respectively, during aging from 2 to 6 WAI (Table 3). Pronounced meristem bifurcation (32.52 and 23.01%) was observed in 4- and 6-week-old nodules of GA3-treated plants, respectively (Figures 6D,F and Table 3). Despite the presence of the senescence zone occupying approximately 34% of the entire nodule, the biggest part of the nodules in GA3-treated plants was represented by the nitrogen fixation zone (Figure 6F and Table 4).
The average DW and projection area of whole 2-week-old nodules of GA3-treated plants were 2.1 and 1.8 times, respectively, larger than those of nodules of untreated plants (Tables 3, 4). At 2 WAI, visual manifestation of the senescence zone (green color) was absent in nodules of GA3-treated and untreated plants (Figures 6A,B). At 4 WAI, the average nodule DW and projection area of nodules of plants treated with exogenous GA3 were 2.5 and 2 times larger, respectively, than the nodules of the untreated plants (Tables 3, 4). In addition to the increase in size, the percentage of nodules exhibiting meristem bifurcation was 3.6 times greater in GA3-treated plants than in untreated plants at 4 WAI (Figures 6C,D and Table 3). At 6 WAI, GA3-treated plants also had a high percentage of nodules with branching due to meristem bifurcation being 6.3 times higher than that in untreated plants (Figure 6F and Table 3). At 6 WAI, nodules of plants treated with exogenous GA3 were 2.3 and 1.6 times (average nodule DW and nodule projection area, respectively) larger than the nodules of the untreated plants (Figure 6F and Tables 3, 4). In addition, in 6-week-old GA3-treated plants, the number of nodules without visible signs of senescence was 2.8 times higher than that of untreated plants (Table 3). Also, it was noted that the number of nodules of plants treated with exogenous GA3 was 1.5, 2.4, and 2.1 times smaller than that of the nodules of untreated plants at 2, 4, and 6 WAI, respectively (Table 3). In GA3-treated plants, a 1.6-fold increase in shoot DW was observed at 4 WAI relative to the control (Table 3). At 2 and 6 WAI, these differences were insignificant.
Microscopic Analysis of Nodule Structure of GA3-Treated and Untreated Plants
Longitudinal toluidine blue-stained sections of nodules of untreated and GA3-treated plants were compared at 6 WAI. The nodules of untreated and GA3-treated plants had the typical zonation of indeterminate nodules (Figure 7). In 6-week-old nodules of untreated plants, the senescence zone, which was characterized by flattened cells, loss of cell rigidity and cellular integrity, and plant cell wall breakage, occupied most of the nodule and contained a cavity due to complete decay of nodule cells (Figure 7A). While nodules of GA3-treated plants had a vast nitrogen fixation zone, the senescence zone contained both infected and uninfected cells and occupied a smaller part than that in the nodules of untreated plants (Figure 7B).
Expression Analysis of Senescence-Associated Genes in Wild-Type Nodules Treated With Exogenous GA3
In this study, an analysis of the expression of senescence-associated genes (Serova et al., 2017, 2018) and genes encoding enzymes of GA biosynthesis (PsGA20ox1) and deactivation (PsGA2ox1) was carried out. Changes in the mRNA levels of these genes during GA3 treatment allowed us to evaluate the possible effect of GAs on nodule senescence at the level of transcription.
A significant increase in the expression levels of selected cysteine protease genes was observed during nodule aging of the untreated plants. However, elevation of transcript abundance in nodules of GA3-treated plants was less pronounced (Figures 8A–C). In 6-week-old nodules of the untreated plants, expression levels of PsCyp1, PsCyp15a, and PsTPP genes were significantly increased relative to 2-week-old nodules, while in nodules of GA3-treated plants, no statistically significant differences in PsCyp1 and PsTPP transcript abundances were detected during nodule aging (Figures 8A,C). At 6 WAI, an up-regulation in PsCyp15a mRNA was detected in nodules of GA3-treated plants; however, the expression level was two times lower than in the untreated control (Figure 8B). In the case of PsCyp1 and PsTPP genes, expression levels in mature nodules of GA3-treated plants were 3- and 2.4-times lower, respectively, than in the untreated plants (Figures 8A,C).
A significant up-regulation of PsATB2 transcript abundance was observed in the nodules of untreated plants from 4 to 6 WAI. A less pronounced elevation of expression was detected during aging of nodules of GA3-treated plants, in comparison with the untreated control. Also, it was noted that at 6 WAI, PsATB2 transcript abundance was two times lower in nodules of GA3-treated plants than in the untreated plants (Figure 8D).
During nodule aging of the untreated plants, the expression levels of genes encoding key enzymes of ethylene biosynthesis, PsACS2 and PsACO1, were significantly increased from 2 to 6 WAI. A less pronounced up-regulation of PsACS2 and PsACO1 mRNA levels was observed in 6-week-old nodules of GA3-treated plants in comparison with 2-week-old nodules. At 6 WAI, PsACS2 and PsACO1 expression levels were 7.6- and 1.7-times lower, respectively, in nodules of GA3-treated plants than in those of untreated plants (Figures 9A,B).
A significant decrease in the expression level of the PsGA20ox1 gene was observed during nodule aging in both untreated and GA3-treated plants. It is worth noting that in 6-week-old nodules of GA3-treated plants, PsGA20ox1 transcript abundance was 2.15-times lower than that in the nodules of untreated plants (Figure 9C). At 6 WAI, a significant elevation in the expression of PsGA2ox1 was observed in the nodules of untreated and GA3-treated plants, relative to 2 WAI. However, the difference between nodules of untreated and GA3-treated plants at 6 WAI was not statistically significant. In 2-week-old nodules of GA3-treated plants, however, PsGA2ox1 transcript abundance was slightly higher than in the nodules of untreated plants (Figure 9D).
Finally, up-regulated expression of PsAO3, the gene encoding an enzyme of the final stage of ABA biosynthesis, was detected in the nodules of untreated and GA3-treated plants from 2 to 6 WAI. In 6-week-old nodules of GA3-treated plants, PsAO3 transcript abundance was slightly reduced relative to that in untreated plants (Figure 9E).
Thus, a down-regulation of senescence-associated genes, a decrease of senescence zone, an extension of nitrogen-fixation zone, an increase of nodule size, and the pronounced meristem bifurcation were observed in wild-type nodules treated with exogenous GA3 in comparison with the untreated plants.
Discussion
In this study, we investigated the involvement of GAs in senescence of the symbiotic nodule using SGE and corresponding pea mutants blocked at different stages of nodule development. The mutant SGEFix--1 (sym40) forms numerous white nodules with hypertrophied infection threads and droplets, and abnormal bacteroids (Tsyganov et al., 1998). White nodules colonized with “locked” suberinized infection threads are typical for SGEFix--2 (sym33) (Tsyganov et al., 1998; Ivanova et al., 2015). SGEFix--2 carries a weak allele of the gene sym33 and manifests a leaky phenotype. Usually, bacteria are not released from these infection threads; however, in some cells or some nodules, bacterial release occurs (Voroshilova et al., 2001). The mutants SGEFix--3 (sym26) and SGEFix--7 (sym27) are characterized by premature degradation of symbiotic structures, with degradation more pronounced in the former than in the latter (Serova et al., 2018). Recently, we demonstrated activation of nodule senescence in nodules of all analyzed mutants (Serova et al., 2018).
Previously, an increase in the expression level of GA biosynthesis genes (GA20ox, GA3ox) during early stages of the development of both determinate (Kouchi et al., 2004; Lievens et al., 2005; Maekawa et al., 2009; Hayashi et al., 2012) and indeterminate (Larrainzar et al., 2015) nodules was shown. A decrease in the transcript abundance of the GA biosynthesis genes was observed in mature nodules (Lievens et al., 2005; Hayashi et al., 2012). In our study, a significant decline in the transcript level of the PsGA20ox1 gene was detected in wild-type nodules at 4 and 6 WAI (Figure 10A). One of the mechanisms to maintain optimal concentrations of bioactive GAs is 2β-hydroxylation, which leads to a decrease in content of bioactive GAs (Thomas et al., 1999; Hedden and Thomas, 2012). GA2ox transcript levels were found to be elevated during the aging of symbiotic nodules of M. truncatula (Van de Velde et al., 2006) and pea (Serova et al., 2017) plants. Furthermore, GA2ox transcripts were detected in maturing nodules of L. japonicus (Kouchi et al., 2004). In this study, the PsGA2ox1 mRNA level was elevated at 6 WAI (Figure 10A). Thus, we suggest that during nodule aging, bioactive GA levels decrease; maintaining the optimal concentration of bioactive GAs appears to occur via the down-regulation of its biosynthetic genes and the up-regulation of its deactivation genes. The results of expression analysis were complemented by those obtained from immunolocalization of bioactive GA3. In the wild-type plants, a high intensity of GA3 label was detected in 2-week-old nodules and was maintained in mature 4-week-old nodules, but it was significantly decreased in old 6-week-old nodules (Figure 10A). The high GA3 content in meristematic cells may indicate the involvement of GAs in cell cycle activation, cell division, and persistence of the nodule meristem. Previously, it was shown that GAs modulate cell cycle activity in Arabidopsis roots (Achard et al., 2009). In addition, pea mutants impaired in GA biosynthetic genes displayed decreased nodulation and aberrant nodule meristem formation (Ferguson et al., 2005, 2011).
Figure 10. Scheme of endogenous gibberellin influence on nodule senescence in wild-type (A) and mutant SGEFix--3 (sym26) nodules (B) and influence of exogenous GA3 on nodule senescence and meristem bifurcation in wild-type (C). Zones of nodule are designated by Roman numerals: I – meristem (in gray), II – infection zone (in gray), III – fixation zone, (in red) III’ – zone corresponding to nitrogen fixation zone in wild-type (in red), IV – senescence zone (in green), WAI – weeks after inoculation. An increase in the color intensity of the arrows indicates elevated levels of gene expression, a decrease in color intensity indicates a reduction in expression levels. A decrease in orange intensity indicates a reduction in GA levels. The scheme is based on the results of this study and Serova et al. (2018).
A high intensity of GA3 labeling was detected in the infected cells in the nitrogen fixation zone, where it was associated mainly with cytoplasm and symbiosomes. In juvenile symbiosomes, the amount of label was higher than in mature and senescent symbiosomes. It is worth noting that McAdam et al. (2018) recently revealed that GAs promote the functioning of nitrogen-fixing nodules. However, the exact function of GAs in infected cell and symbiosome development merits elucidation. Previously, a negative impact of GAs on infection thread development was revealed (McAdam et al., 2018). In our study, we observed that the amount of GA3 in the infection threads was lower than in symbiosomes. However, GA3 was detected in infection thread walls. The function of GAs in infection thread development is currently unknown. It is interesting to note that GA involvement in facilitating bacterial release from the infection threads was recently suggested (Tatsukami and Ueda, 2016).
In uninfected cells, labeling of GA3 was associated with cytoplasm around starch granules. Ferguson and Mathesius (2003) suggested involvement of GAs in the hydrolysis of nodule starch through enhancement of α-amylase production. Substrates formed during this reaction might support the energetic requirements of rhizobia. In general, GAs may benefit both symbiotic partners through involvement in cell division and elongation, and by providing energy to support plant growth. GAs may enhance photosynthetic capacity, leading to an increase in photoassimilate content for growth and functioning of the symbiotic nodule (Ferguson and Mathesius, 2003). Six-week-old wild-type nodules contained a low amount of bioactive GA3 (Figure 10A). Traces of GA3 label were detected in the uninfected cells of the senescence zone; these cells might be degraded later that the infected cells of that same zone (Pladys et al., 1991).
Expression analysis of GA metabolism genes revealed down-regulation of a GA biosynthesis gene, PsGA20ox1, and significant up-regulation of a GA deactivation gene, PsGA2ox1, during nodule aging of the early senescent mutants SGEFix--7 (sym27) and SGEFix--3 (sym26), which confirms involvement of GAs in senescence of symbiotic nodules (Figure 10B). However, for mutant SGEFix--3 (sym26), an increase in PsGA20ox1 expression in 4-week-old nodules in comparison with 2-week-old nodules was observed. This increase is difficult to explain, but future elucidation of the gene sym26 molecular product may be informative. Data obtained with GA3 immunolocalization coincided with the observed transcriptional patterns. In 2-week-old nodules of the mutants SGEFix--3 (sym26) and SGEFix--7 (sym27), the maximum GA3 signal was observed in cells of the region corresponding to the nitrogen fixation zone of the wild-type nodules (Figure 10B). In the vast senescence zone of 4-week-old nodules, the intensity of GA3 label was low and predominantly found in uninfected cells. Both mutants form nodules with morphologically differentiated bacteroids, which characterized with early senescence (Serova et al., 2018). Recently, it was shown that in defective nodules of GA-deficient mutant na-1, bacteroids undergo premature degradation (McAdam et al., 2018). Together, these data strongly indicate that aging of symbiotic nodules is accompanied by a decrease in the level of GAs, which suggests involvement of GAs in the delay of nodule senescence (Figures 10A,B).
No significant differences in PsGA20ox1 and PsGA2ox1 transcript abundance were seen during nodule aging of the mutant SGEFix--2 (sym33), blocked at the earliest stage of symbiosis. In addition, a low level of GA3 labeling in nodules of the mutant SGEFix--2 (sym33) was observed. Recently, we have shown early activation of senescence-associated genes in this mutant at 2 WAI (Serova et al., 2018). This activation may explain the low level of GAs in the mutant. The low level of bioactive GAs in the mutant SGEFix--2 (sym33) is accompanied by the absence of infected cells due to the absence of bacterial release (Tsyganov et al., 1998). These results may reflect the necessity of GAs for bacterial release (Tatsukami and Ueda, 2016) and infected cell development.
Analysis of the mutant SGEFix--1 (sym40), blocked at a later stage of symbiosis than mutant SGEFix--2 (sym33) (Tsyganov et al., 2011) and manifesting signs of premature nodule senescence (Tsyganov et al., 1998; Serova et al., 2018), demonstrated a similar transcriptional pattern of PsGA20ox1 to that observed in wild-type nodules and a significant increase in PsGA2ox1 mRNA as well as low GA3 label content in 4-week-old nodules. These data are suggestive of a negative role of GAs during nodule senescence.
Aging of the symbiotic nodule is regulated by changes in gene expression. In particular, the expression levels of genes encoding cysteine proteases, transcription factors, enzymes of GA deactivation, and enzymes of ethylene and ABA biosynthesis were shown to increase during senescence of M. truncatula and pea nodules (Kardailsky and Brewin, 1996; Van de Velde et al., 2006; D’haeseleer et al., 2010; de Zélicourt et al., 2012; Karmarkar, 2014; Serova et al., 2017, 2018). To assess the effect of GAs on the senescence of the pea symbiotic nodules, we analyzed the mRNA levels of senescence-associated genes in 2-, 4-, and 6-week-old nodules of wild-type plants treated with exogenous GA3 relative to untreated plants. A less pronounced increase in transcript levels of genes PsCyp1, PsCyp15, and PsTPP was shown during aging of the nodules of GA3-treated plants, in contrast to the nodules of the untreated plants. In addition, the transcript abundance of all analyzed genes was significantly higher in the nodules of untreated plants than in nodules of plants treated with GA3 at 6 WAI, which may indicate a delay of nodule senescence upon treatment with exogenous GA3 (Figure 10C). It is known that cysteine proteases carry out large-scale protein degradation during nodule senescence (Pladys et al., 1991; Granell et al., 1992; Kardailsky and Brewin, 1996; Van de Velde et al., 2006; Pérez Guerra et al., 2010). Thus, our data suggest that a decrease in GAs is required to induce the degradation processes during nodule aging. Previously, it was shown that the expression level of the ATB2 gene was up-regulated during senescence of M. truncatula (D’haeseleer et al., 2010) and pea (Serova et al., 2017) nodules. The transcript level of PsATB2 was reduced in the nodules of GA3-treated plants in comparison with the nodules of untreated plants (Figure 10C). This may indicate possible regulation of nodule senescence by GAs through an effect on the bZIP transcription factor. Similar expression patterns were observed for genes encoding key enzymes of ethylene biosynthesis, ACC synthase (PsACS2), ACC oxidase (PsACO1), and aldehyde oxidase (PsAO3), an enzyme catalyzing the last step of ABA biosynthesis (Figure 10C). It is known that ethylene and ABA promote nodule senescence (Swamy and Smith, 1999; Guinel, 2015; Serova et al., 2017, 2018). In addition, elevation of ethylene levels in the na-1 mutant was previously demonstrated (Ferguson et al., 2011). Thus, a lower expression level of PsACS2, PsACO1, and PsAO3 genes in the nodules of GA3-treated plants may also suggest integration of hormonal signaling during nodule senescence. During the nodule aging of both GA3-treated and untreated plants, down-regulation of one GA biosynthetic gene, PsGA20ox1, and up-regulation of one GA deactivation gene, PsGA2ox1, were demonstrated (Figure 10C). However, only PsGA2ox1 transcript abundance was significantly lower in the nodules of GA3-treated plants than in those of untreated plants at 6 WAI. Thus, it can be assumed that regulation of bioactive GAs in the nodules of GA3-treated plants occurred mainly by regulation of GA biosynthesis.
Two-week-old nodules of GA3-treated plants were about two times larger than the nodules of untreated plants, probably due to involvement of GAs in cell division and cell elongation (Hedden and Thomas, 2012). The number of nodules formed by GA3-treated plants was nearly two times smaller than that in untreated plants. Previously, a similar concentration (10-6 M) of GA3 had no effect on nodule number in the wild type, but a higher concentration (10-3 M) inhibited nodule formation (Ferguson et al., 2005). This contradiction may be due to differences in plant genotypes and experimental conditions. However, GA3 treatment of the severely inhibited GA mutant na-1, characterized as having a low number of nodules, restored nodule number, which suggests a direct role for GAs in nodule development (Ferguson et al., 2005). Most of the GA mutants were characterized by a reduced number of nodules, but increased nodule dry weight. It has been suggested by Ferguson et al. (2005) that there is a compensation mechanism to regulate the size of individual nodules, depending on the number of nodules per plant.
Pronounced meristem bifurcation was observed in 4- and especially 6-week-old nodules of GA3-treated plants (Figure 10C). This may indicate the involvement of GAs in functioning of the nodule meristem. Previously, it was shown that meristem bifurcation during root branching is under phytohormonal control, including auxin, cytokinins, and ethylene (Gola, 2014). GA3 immunolocalization in pea nodules performed in this study and previous studies carried out on Arabidopsis and pea plants (Ferguson et al., 2005, 2011; Achard et al., 2009) confirms the involvement of GAs in nodule meristem functioning. It should be noted that 6-week-old nodules of GA3-treated plants also had a senescence zone, but it occupied a smaller part (about 34%) of the nodule than in nodules of the untreated plants (about 62%). This may be due to the meristem bifurcation observed in the nodules of GA3-treated plants. It is known that senescence of the nodule is associated with arrested division of meristem cells and the discontinuance of rhizobial release from infection droplets (Guinel, 2015).
The main part of mature nodules of GA3-treated plants was represented by the nitrogen fixation zone (Figure 10C). Also, the high content of GA3 label in cells of the nitrogen fixation zone in young and mature pea nodules indicates a possible involvement of GAs in the functioning of nitrogen fixation nodules and, consequently, delay of nodule senescence. In A. thaliana, it was shown that DELLAs repress the inhibitory effect of JASMONATE ZIM-domain proteins (JAZ), which are transcriptional regulators, on the expression of jasmonate (JA) responsive genes, such as lipoxygenase and defense genes (Hou et al., 2010). On the contrary, GAs suppress the expression of JA-responsive genes via DELLA degradation. Thus, the involvement of GAs in suppressing the response to JA-induced signaling may be one of the mechanisms by which GAs contribute to the delay of nodule senescence, which is regarded as a delayed response of the plant to rhizobia as a potential pathogen (Mellor, 1989). Furthermore, an up-regulation of transcripts of JA biosynthesis genes was observed during M. truncatula nodule senescence (Van de Velde et al., 2006).
Conclusion
In this study, the involvement of bioactive GAs in nodule senescence of pea wild-type and nodule development mutants was studied by assessing transcriptional patterns of GA metabolism genes, and through GA3 immunolocalization and pharmacological analyses. A decrease in GA content during nodule aging was demonstrated at the transcriptional level via a down-regulation of the GA biosynthesis gene, PsGA20ox1, and an up-regulation of the GA deactivation gene, PsGA2ox1, and also by the immunolocalization of bioactive GA3 in the mutant and wild-type nodules. These results indicate a role of GAs in a delay of nodule senescence. A down-regulation of senescence-associated genes, a decrease of the senescence zone, and an increase of the nitrogen fixation zone in nodules of wild-type plants treated with exogenous GA3 confirm a negative regulation of nodule senescence by GAs and involvement of GAs in the functioning of the mature nodule.
Data Availability
All datasets generated for this study are included in the manuscript and/or the Supplementary Files.
Author Contributions
IT designed the experiments. TS and AT performed the experiments. TS and VT analyzed the data and wrote the manuscript.
Funding
This work was supported by Russian Science Foundation (grant 17-76-30016).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The research was performed using equipment of the Core Centrum “Genomic Technologies, Proteomics and Cell Biology” in ARRIAM and the “Molecular and Cell Technologies” Research Resource Centre at Saint Petersburg State University. We thank Charlesworth Author Services for English editing of a draft of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00285/full#supplementary-material
References
Achard, P., Gusti, A., Cheminant, S., Alioua, M., Dhondt, S., Coppens, F., et al. (2009). Gibberellin signaling controls cell proliferation rate in Arabidopsis. Curr. Biol. 19, 1188–1193. doi: 10.1016/j.cub.2009.05.059
Brewin, N. J. (2004). Plant cell wall remodelling in the Rhizobium–legume symbiosis. Crit. Rev. Plant Sci. 23, 293–316. doi: 10.1080/07352680490480734
Davière, J.-M., and Achard, P. (2013). Gibberellin signaling in plants. Development 140, 1147–1151. doi: 10.1242/dev.087650
de Zélicourt, A., Diet, A., Marion, J., Laffont, C., Ariel, F., Moison, M., et al. (2012). Dual involvement of a Medicago truncatula NAC transcription factor in root abiotic stress response and symbiotic nodule senescence. Plant J. 70, 220–230. doi: 10.1111/j.1365-313X.2011.04859.x
D’haeseleer, K., De Keyser, A., Goormachtig, S., and Holsters, M. (2010). Transcription factor MtATB2: about nodulation, sucrose and senescence. Plant Cell Physiol. 51, 1416–1424. doi: 10.1093/pcp/pcq104
Dupont, L., Alloing, G., Pierre, O., El Msehli, S., Hopkins, J., Hérouart, D., et al. (2012). “The legume root nodule: from symbiotic nitrogen fixation to senescence,” in Senescence, ed. T. Nagata (Rijeka: IntechOpen), 137–168. doi: 10.5772/34438
Ferguson, B. J., Foo, E., Ross, J. J., and Reid, J. B. (2011). Relationship between gibberellin, ethylene and nodulation in Pisum sativum. New Phytol. 189,829–842. doi: 10.1111/j.1469-8137.2010.03542.x
Ferguson, B. J., and Mathesius, U. (2003). Signaling interactions during nodule development. J. Plant Growth Regul. 22, 47–72. doi: 10.1007/s00344-003-0032-9
Ferguson, B. J., and Mathesius, U. (2014). Phytohormone regulation of legume-rhizobia interactions. J. Chem. Ecol. 40, 770–790. doi: 10.1007/s10886-014-0472-7
Ferguson, B. J., Ross, J. J., and Reid, J. B. (2005). Nodulation phenotypes of gibberellin and brassinosteroid mutants of pea. Plant Physiol. 138, 2396–2405. doi: 10.1104/pp.105.062414
Fonouni-Farde, C., Tan, S., Baudin, M., Brault, M., Wen, J., Mysore, K. S., et al. (2016). DELLA-mediated gibberellin signalling regulates Nod factor signalling and rhizobial infection. Nat. Commun. 7:12636. doi: 10.1038/ncomms12636
García-Martínez, J. L., López-Diaz, I., Sánchez-Beltrán, M. J., Phillips, A. L., Ward, D. A., Gaskin, P., et al. (1997). Isolation and transcript analysis of gibberellin 20-oxidase genes in pea and bean in relation to fruit development. Plant Mol. Biol. 33, 1073–1084. doi: 10.1023/A:1005715722193
Gola, E. M. (2014). Dichotomous branching: the plant form and integrity upon the apical meristem bifurcation. Front. Plant Sci. 5:263. doi: 10.3389/fpls.2014.00263
Granell, A., Harris, N., Pisabarro, A. G., and Carbonell, J. (1992). Temporal and spatial expression of a thiolprotease gene during pea ovary senescence, and its regulation by gibberellin. Plant J. 2, 907–915. doi: 10.1046/j.1365-313X.1992.t01-5-00999.x
Guinel, F. C. (2009). Getting around the legume nodule: I. The structure of the peripheral zone in four nodule types. Botany 87, 1117–1138. doi: 10.1139/B09-074
Guinel, F. C. (2015). Ethylene, a hormone at the center-stage of nodulation. Front. Plant Sci. 6:1121. doi: 10.3389/fpls.2015.01121
Hayashi, S., Gresshoff, P. M., and Ferguson, B. J. (2014). Mechanistic action of gibberellins in legume nodulation. J. Integr. Plant Biol. 56, 971–978. doi: 10.1111/jipb.12201
Hayashi, S., Reid, D. E., Lorenc, M. T., Stiller, J., Edwards, D., Gresshoff, P. M., et al. (2012). Transient Nod factor-dependent gene expression in the nodulation-competent zone of soybean (Glycine max [L.] Merr.) roots. Plant Biotechnol. J. 10, 995–1010. doi: 10.1111/j.1467-7652.2012.00729.x
Hedden, P., and Phillips, A. L. (2000). Gibberellin metabolism: new insights revealed by the genes. Trends Plant Sci. 5, 523–530. doi: 10.1016/S1360-1385(00)01790-8
Hedden, P., and Thomas, S. G. (2012). Gibberellin biosynthesis and its regulation. Biochem. J. 444, 11–25. doi: 10.1042/BJ20120245
Hou, X., Lee, L. Y. C., Xia, K., Yan, Y., and Yu, H. (2010). DELLAs modulate jasmonate signaling via competitive binding to JAZs. Dev. Cell 19, 884–894. doi: 10.1016/j.devcel.2010.10.024
Ivanova, K. A., Tsyganova, A. V., Brewin, N. J., Tikhonovich, I. A., and Tsyganov, V. E. (2015). Induction of host defences by Rhizobium during ineffective nodulation of pea (Pisum sativum L.) carrying symbiotically defective mutations sym40 (PsEFD), sym33 (PsIPD3/PsCYCLOPS) and sym42. Protoplasma 252, 1505–1517. doi: 10.1007/s00709-015-0780-y
Jin, Y., Liu, H., Luo, D., Yu, N., Dong, W., Wang, C., et al. (2016). DELLA proteins are common components of symbiotic rhizobial and mycorrhizal signalling pathways. Nat. Commun. 7:12433. doi: 10.1038/ncomms12433
Kardailsky, I. V., and Brewin, N. J. (1996). Expression of cysteine protease genes in pea nodule development and senescence. Mol. Plant Microbe Interact. 9, 689–695. doi: 10.1094/MPMI-9-0689
Karmarkar, V. M. (2014). Transcriptional Regulation of Nodule Development and Senescence in Medicago Truncatula. Doctoral thesis, Wageningen University, Wageningen, D.C.
Kawaguchi, M., Imaizumi-Anraku, H., Fukai, S., and Syono, K. (1996). Unusual branching in the seedlings of Lotus japonicus—gibberellins reveal the nitrogen-sensitive cell divisions within the pericycle on roots. Plant Cell Physiol. 37, 461–470. doi: 10.1093/oxfordjournals.pcp.a028968
Kosterin, O. E., and Rozov, S. M. (1993). Mapping of the new mutation blb and the problem of integrity of linkage group I. Pisum Genet. 25, 27–31.
Kouchi, H., Shimomura, K., Hata, S., Hirota, A., Wu, G.-J., Kumagai, H., et al. (2004). Large-scale analysis of gene expression profiles during early stages of root nodule formation in a model legume, Lotus japonicus. DNA Res. 11, 263–274. doi: 10.1093/dnares/11.4.263
Larrainzar, E., Riely, B. K., Kim, S. C., Carrasquilla-Garcia, N., Yu, H.-J., Hwang, H.-J., et al. (2015). Deep sequencing of the Medicago truncatula root transcriptome reveals a massive and early interaction between nodulation factor and ethylene signals. Plant Physiol. 169, 233–265. doi: 10.1104/pp.15.00350
Lester, D. R., Ross, J. J., Smith, J. J., Elliott, R. C., and Reid, J. B. (1999). Gibberellin 2-oxidation and the SLN gene of Pisum sativum. Plant J. 19, 65–73. doi: 10.1046/j.1365-313X.1999.00501.x
Lievens, S., Goormachtig, S., Den Herder, J., Capoen, W., Mathis, R., Hedden, P., et al. (2005). Gibberellins are involved in nodulation of Sesbania rostrata. Plant Physiol. 139, 1366–1379. doi: 10.1104/pp.105.066944
Maekawa, T., Maekawa-Yoshikawa, M., Takeda, N., Imaizumi-Anraku, H., Murooka, Y., and Hayashi, M. (2009). Gibberellin controls the nodulation signaling pathway in Lotus japonicus. Plant J. 58, 183–194. doi: 10.1111/j.1365-313X.2008.03774.x
Martin, D. N., Proebsting, W. M., and Hedden, P. (1999). The SLENDER gene of pea encodes a gibberellin 2-oxidase. Plant Physiol. 121, 775–781. doi: 10.1104/pp.121.3.775
McAdam, E. L., Reid, J. B., and Foo, E. (2018). Gibberellins promote nodule organogenesis but inhibit the infection stages of nodulation. J. Exp. Bot. 69, 2117–2130. doi: 10.1093/jxb/ery046
Mellor, R. B. (1989). Bacteroids in the Rhizobium-legume symbiosis inhabit a plant internal lytic compartment: implications for other microbial endosymbioses. J. Exp. Bot. 40, 831–839. doi: 10.1093/jxb/40.8.831
Nemankin, N. (2011). Analysis of Pea (Pisum sativum L.) Genetic System, Controlling Development of Arbuscular Mycorrhiza and Nitrogen-Fixing Symbiosis. Ph. D thesis, Saint-Petersburg State University, Saint-Petersburg, Russia.
Ovchinnikova, E., Journet, E.-P., Chabaud, M., Cosson, V., Ratet, P., Duc, G., et al. (2011). IPD3 controls the formation of nitrogen-fixing symbiosomes in pea and Medicago Spp. Mol. Plant Microbe Interact. 24, 1333–1344. doi: 10.1094/MPMI-01-11-0013
Pérez Guerra, J. C., Coussens, G., De Keyser, A., De Rycke, R., De Bodt, S., Van De Velde, W., et al. (2010). Comparison of developmental and stress-induced nodule senescence in Medicago truncatula. Plant Physiol. 152, 1574–1584. doi: 10.1104/pp.109.151399
Pladys, D., Dimitrijevic, L., and Rigaud, J. (1991). Localization of a protease in protoplast preparations in infected cells of French bean nodules. Plant Physiol. 97, 1174–1180. doi: 10.1104/pp.97.3.1174
Puppo, A., Groten, K., Bastian, F., Carzaniga, R., Soussi, M., Lucas, M. M., et al. (2005). Legume nodule senescence: roles for redox and hormone signalling in the orchestration of the natural aging process. New Phytol. 165, 683–701. doi: 10.1111/j.1469-8137.2004.01285.x
Ross, J. J., Reid, J. B., Swain, S. M., Hasan, O., Poole, A. T., Hedden, P., et al. (1995). Genetic regulation of gibberellin deactivation in Pisum. Plant J. 7, 513–523. doi: 10.1046/j.1365-313X.1995.7030513.x
Serova, T. A., Tikhonovich, I. A., and Tsyganov, V. E. (2017). Analysis of nodule senescence in pea (Pisum sativum L.) using laser microdissection, real-time PCR, and ACC immunolocalization. J. Plant Physiol. 212, 29–44. doi: 10.1016/j.jplph.2017.01.012
Serova, T. A., and Tsyganov, V. E. (2014). Symbiotic nodule senescence in legumes: molecular-genetic and cellular aspects (review). Agric. Biol. 5, 3–15. doi: 10.15389/agrobiology.2014.5.3eng
Serova, T. A., Tsyganova, A. V., and Tsyganov, V. E. (2018). Early nodule senescence is activated in symbiotic mutants of pea (Pisum sativum L.) forming ineffective nodules blocked at different nodule developmental stages. Protoplasma 255, 1443–1459. doi: 10.1007/s00709-018-1246-9
Sun, T. P. (2011). The molecular mechanism and evolution of the GA–GID1–DELLA signaling module in plants. Curr. Biol. 21, R338–R345. doi: 10.1016/j.cub.2011.02.036
Swamy, P., and Smith, B. N. (1999). Role of abscisic acid in plant stress tolerance. Curr. Sci. 76, 1220–1227.
Tatsukami, Y., and Ueda, M. (2016). Rhizobial gibberellin negatively regulates host nodule number. Sci. Rep. 6:27998. doi: 10.1038/srep27998
Thomas, S. G., Phillips, A. L., and Hedden, P. (1999). Molecular cloning and functional expression of gibberellin 2-oxidases, multifunctional enzymes involved in gibberellin deactivation. Proc. Natl. Acad. Sci. U.S.A. 96, 4698–4703. doi: 10.1073/pnas.96.8.4698
Tsyganov, V. E., Borisov, A. Y., Rozov, S. M., and Tikhonovich, I. A. (1994). New symbiotic mutants of pea obtained after mutagenesis of laboratory line SGE. Pisum Genet. 26, 36–37.
Tsyganov, V. E., Morzhina, E. V., Stefanov, S. Y., Borisov, A. Y., Lebsky, V. K., and Tikhonovich, I. A. (1998). The pea (Pisum sativum L.) genes sym33 and sym40 control infection thread formation and root nodule function. Mol. Gen. Genet. 259, 491–503. doi: 10.1007/s004380050840
Tsyganov, V. E., Seliverstova, E. V., Voroshilova, V. A., Tsyganova, A. V., Pavlova, Z. B., Lebskii, V. K., et al. (2011). Double mutant analysis of sequential functioning of pea (Pisum sativum L.) genes Sym13, Sym33, and Sym40 during symbiotic nodule development. Russ. J. Genet. Appl. Res. 1, 343–348. doi: 10.1134/S2079059711050145
Tsyganov, V. E., Voroshilova, V. A., Borisov, A. Y., Tikhonovich, I. A., and Rozov, S. M. (2000). Four more symbiotic mutants obtained using EMS mutagenesis of line SGE. Pisum Genet. 32:63.
Tsyganov, V. E., Voroshilova, V. A., Rozov, S. M., Borisov, A. Y., and Tikhonovich, I. A. (2013). A new series of pea symbiotic mutants induced in the line SGE. Russ. J. Genet. Appl. Res. 3, 156–162. doi: 10.1134/S2079059713020093
Tsyganova, A. V., Kitaeva, A. B., and Tsyganov, V. E. (2018). Cell differentiation in nitrogen-fixing nodules hosting symbiosomes. Funct. Plant Biol. 45, 47–57. doi: 10.1071/FP16377
Tsyganova, A. V., and Tsyganov, V. E. (2015). Negative hormonal regulation of symbiotic nodule development. I. Ethylene (review). Agric. Biol. 50, 267–277. doi: 10.15389/agrobiology.2015.3.267eng
Tsyganova, A. V., and Tsyganov, V. E. (2018). Negative hormonal regulation of symbiotic nodule development. II. Salicilic, jasmonic and abscisic acids (review). Agric. Biol. 53, 3–14. doi: 10.15389/agrobiology.2018.1.3eng
Tsyganova, A. V., and Tsyganov, V. E. (2017). “Plant genetic control over infection thread development during legume-Rhizobium symbiosis,” in Symbiosis, ed. E. C. Rigobelo (London: IntechOpen), 23–52. doi: 10.5772/intechopen.70689
Tsyganova, A. V., Tsyganov, V. E., Findlay, K.C., Borisov, A. Y., Tikhonovich, I. A., and Brewin, N. G. (2009). Distribution of legume arabinogalactanprotein-extensin (AGPE) glycoproteins in symbiotically defective pea mutants with abnormal infection threads. Cell Tissue Biol. 51, 53–62. doi: 10.1134/S1990519X09010131
Ueguchi-Tanaka, M., Ashikari, M., Nakajima, M., Itoh, H., Katoh, E., Kobayashi, M., et al. (2005). GIBBERELLIN INSENSITIVE DWARF1 encodes a soluble receptor for gibberellin. Nature 437, 693–698. doi: 10.1038/nature04028
Van de Velde, W., Guerra, J. C. P., Keyser, A. D., De Rycke, R., Rombauts, S., Maunoury, N., et al. (2006). Aging in legume symbiosis. A molecular view on nodule senescence in Medicago truncatula. Plant Physiol. 141, 711–720. doi: 10.1104/pp.106.078691
Voroshilova, V. A., Boesten, B., Tsyganov, V. E., Borisov, A. Y., Tikhonovich, I. A., and Priefer, U. B. (2001). Effect of mutations in Pisum sativum L. Genes blocking different stages of nodule development on the expression of late symbiotic genes in Rhizobium leguminosarum bv. viciae. Mol. Plant Microbe Interact. 14, 471–476. doi: 10.1094/mpmi.2001.14.4.471
Wang, T. L., Wood, E. A., and Brewin, N. J. (1982). Growth regulators, Rhizobium and nodulation in peas. Planta 155, 345–349. doi: 10.1007/bf00429463
Keywords: Rhizobium–legume symbiosis, nodule senescence, cysteine protease, ethylene, abscisic acid, gibberellins, meristem
Citation: Serova TA, Tsyganova AV, Tikhonovich IA and Tsyganov VE (2019) Gibberellins Inhibit Nodule Senescence and Stimulate Nodule Meristem Bifurcation in Pea (Pisum sativum L.). Front. Plant Sci. 10:285. doi: 10.3389/fpls.2019.00285
Received: 26 October 2018; Accepted: 20 February 2019;
Published: 15 March 2019.
Edited by:
Eloise Foo, University of Tasmania, AustraliaReviewed by:
Frederique Catherine Guinel, Wilfrid Laurier University, CanadaErik Limpens, Wageningen University & Research, Netherlands
Copyright © 2019 Serova, Tsyganova, Tikhonovich and Tsyganov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Viktor E. Tsyganov, tsyganov@arriam.spb.ru