Skip to main content

A Lotus japonicus E3 ligase interacts with the Nod Factor Receptor 5 and positively regulates nodulation



Post-translational modification of receptor proteins is involved in activation and de-activation of signalling systems in plants. Both ubiquitination and deubiquitination have been implicated in plant interactions with pathogens and symbionts.


Here we present LjPUB13, a PUB-ARMADILLO repeat E3 ligase that specifically ubiquitinates the kinase domain of the Nod Factor receptor NFR5 and has a direct role in nodule organogenesis events in Lotus japonicus. Phenotypic analyses of three LORE1 retroelement insertion plant lines revealed that pub13 plants display delayed and reduced nodulation capacity and retarded growth. LjPUB13 expression is spatially regulated during symbiosis with Mesorhizobium loti, with increased levels in young developing nodules.


LjPUB13 is an E3 ligase with a positive regulatory role during the initial stages of nodulation in L. japonicus.


The legume-rhizobia symbiosis leads to the formation of novel organs on the plant root, termed nodules. Rhizobia within nodule cells differentiate into bacteroids that fix atmospheric dinitrogen in exchange for plant carbohydrates. The symbiotic signalling process is initiated when rhizobia secrete nodulation (Nod) factors upon sensing flavonoids produced by compatible legumes. Lotus japonicus Nod factor receptors NFR1 and NFR5 and the corresponding proteins LYK3 and NFP in Medicago truncatula and, GmNFR1 and GmNFR5 in soybean [1,2,3,4,5,6,7,8] are crucial for perception of rhizobial Nod factors. Rhizobia enter into roots through infection threads (ITs) that, in most cases, initiate in epidermal root hair cells and progress to inner root tissues (reviewed in [9]). Formation of functional nodules requires two separate but tightly coordinated developmental processes: bacterial infection and nodule organogenesis. Protein ubiquitination has been identified as being critical for these two signalling pathways.

Ubiquitination of proteins, by which proteins are tagged by ubiquitin and subsequently destined to be degraded by the proteasome, is a regulatory process essential for eukaryotic growth, development and response to interacting microbes [10]. In some cases, ubiquitination is important for regulating the activity or trafficking of the target protein [11]. Post-translational modification by ubiquitination is accomplished by a three-step process that involves ATP-dependent activation of ubiquitin by an E1 enzyme, followed by conjugation by an E2 enzyme and specific ubiquitin ligation to substrate proteins by direct interaction by an E3 ligase [12, 13].

E3 ligases are divided into families based on their mechanism of action and on the presence of specific E2 interacting domains, such as HECT, RING and U-box (reviewed in [14, 15]). The largest class of plant U-box (PUB) proteins is the ARMADILLO (ARM) domain-containing PUB proteins that contain tandemly-repeated ARM motifs located at the C-terminal [16, 17]. ARM repeat proteins are known to be involved in a number of different cellular processes including signal transduction, cytoskeletal regulation, nuclear import, transcriptional regulation, and ubiquitination. PUB-ARM proteins have been implicated in plant receptor-like kinase (RLK) signalling, with the ARM repeat domain mediating the binding of PUBs to the kinase domain [17, 18].

E3 ligases have been shown to play roles in the establishment of legume-rhizobium symbiosis (reviewed in [19]). M. truncatula LIN [20] and the orthologous gene from L. japonicus, CERBERUS [21], encode E3 ligases containing U-box, ARM and WD-40 repeats, and have been reported to control rhizobial infection inside root hairs. A second PUB-ARM E3 ubiquitin ligase in M. truncatula, MtPUB1, was identified as a negative regulator of nodulation by direct interaction with the receptor-like kinase LYK3 [22]. M. truncatula PUB1 is required for both rhizobial and arbuscular mycorrhiza (AM) endosymbiosis as it also directly interacts with the receptor kinase DMI2, a key component of the common symbiosis signalling pathway [23]. A member of the SEVEN IN ABSENTIA (SINA) family of E3 ligases, SINA4, was shown to interact with SYMRK receptor-like kinase in L. japonicus and be a negative regulator of rhizobial infection [24]. LjnsRING, a RING-H2 E3 ubiquitin ligase from L. japonicus, was also reported to be required for both rhizobial infection and nodule function [25]. Nevertheless, the mechanistic mode of action of the E3 ligases in the symbiotic interactions has not been elucidated, and the targets of all the symbiotic E3 ligases remain to be identified.

Here, we report the involvement of a PUB-ARM protein, LjPUB13, in the symbiotic interaction of L. japonicus with rhizobia. LjPUB13 is phylogenetically a close relative of Arabidopsis PUB13, which has been implicated in plant responses to bacterial flagellin (flg22) [26] and, recently, fungal long-chain chitooligosaccharides (chitooctaose) [27]. In L. japonicus, LjPUB13 is involved in the establishment of a successful symbiosis, through the interaction of its ARM domain with the Nod factor receptor NFR5 and direct ubiquitination of the NFR5 kinase domain.


LjPUB13 encodes an active E3 ligase

A Lotus japonicus sequence encoding for a novel ARMADILLO repeat-containing protein was initially identified as the TC63883 in the DFCI Gene Index Database (, while searching for ARMADILLO repeat proteins in L. japonicus. Initial search at transcript data in Lotus Base ( showed that this gene is expressed in nodules. The predicted protein has a typical Plant U-box (PUB) E3 ubiquitin ligase sequence with a U-box motif followed by a C-terminal ARMADILLO (ARM) repeat domain (Fig. 1a). It encodes a 671 aminoacid long protein with an estimated molecular mass of 73 kDa. Phylogenetic analysis showed that this protein is distinct from other previously characterised PUBs in legumes, but it is the closest homolog to Arabidopsis thaliana PUB13 (Additional file 1: Figure S1 and S2). We therefore designated the protein as LjPUB13. The encoded polypeptide shares 69% and 65% amino acid identity with AtPUB13 (AT3G46510) and AtPUB12 (AT2G28830), respectively. At the genomic level, this gene displays an exon-intron structure similar to AtPUB13 [28], but different from AtPUB12 (, with four exons spanning a 4547 bp region on chromosome 3 (Fig. 1a).

Fig. 1

LjPUB13 is an active E3 ubiquitin ligase. a L. japonicus PUB13 gene structure. PUB13 protein contains a U-box motif and an ARM-repeat domain. The pub13 mutants carry LORE1 insertions in the coding region of the PUB13 gene. The positions of the LORE1 insertions in the pub13 mutants are marked by arrows. Numbers indicate nucleotides. b LjPUB13 auto-ubiquitination. PUB13 was purified as a GST fusion protein. The ubiquitination was detected by both anti-GST and anti-Ub antibodies. UbcH6 was used as the E2 enzyme in the reactions. The experiment was repeated four times with similar results

The predicted E3 ligase enzymatic activity of LjPUB13 was analysed by an in vitro ubiquitination assay, performed with LjPUB13 purified as a glutathione S-transferase (GST) fusion protein from Escherichia coli. Incubation of LjPUB13 with E1 ubiquitin-activating and E2 ubiquitin-conjugated enzymes, ubiquitin and ATP resulted in LjPUB13 polyubiquitination (Fig. 1b). Bands with a larger molecular weight than LjPUB13-GST and a protein ladder were detected by anti-GST and anti-ubiquitin antibodies, respectively, when LjPUB13 was added to the reaction (Fig. 1b, lane 4). The observed E3 ligase activity was abolished in the absence of E1 or E2 enzymes (Fig. 1b, lanes 1 & 2). These results show that L. japonicus PUB13 is a functional E3 ligase, possessing auto-ubiquitination activity.

The expression of LjPUB13 is symbiotically regulated

To test if LjPUB13 is regulated during symbiosis, we monitored its transcript levels in roots and nodules at 7, 14, 21 and 28 days post inoculation (dpi) with Mesorhizobium loti. No significant increase was detected in inoculated plants compared to uninoculated ones. Functional nodules, in general, had approximately 2 times less LjPUB13 transcript compared to roots (Fig. 2a).

Fig. 2

LjPUB13 is expressed in symbiotically active root and nodules. a LjPUB13 expression in uninoculated L. japonicus roots vs roots inoculated by M. loti at 7, 14, 21 and 28 days post inoculation and 14, 21 and 28-day-old nodules. Transcript levels were normalized to those of UBQ. Bars represent means (+SE) of three biological replications (n = 8). Significant differences (P < 0.05) are indicated by asterisk. b, c, d, f Expression of LjPUB13 in L. japonicus transgenic hairy roots transformed with a ProPUB13::GUS construct, detected after GUS staining. In both uninoculated (b, c) and inoculated roots (d, f), PUB13 promoter activity was observed in the vascular bundle. At uninoculated plants strong LjPUB13 expression is detected in lateral root initiation sites (c, arrows). After rhizobial inoculation (d, f), LjPUB13 promoter was active in nodule promordia and young developing nodules at 7dpi (d, arrows) but not in mature nodules at 21dpi (f). e and g show the tissues colonised by DsRed-expressing rhizobia in (d) and (f), respectively. Bars, 500 μm

To complement the results from gene transcript analysis, we constructed a LjPUB13pro::GUS fusion and analysed the spatial and temporal regulation of LjPUB13 promoter activity in Agrobacterium rhizogenes transformed roots expressing the transgene. In both uninoculated and inoculated roots, the LjPUB13 promoter was active throughout the whole root, but strongest in the vascular bundle (Fig. 2b) and at the lateral root initiation sites (Fig. 2c). In M. loti inoculated plants, an increase in the promoter activity in the root zones, where nodule primordia formation occurs, was observed (Fig. 2d). Although strong LjPUB13Pro::GUS expression was detected in young developing nodules at 7 dpi (Fig. 2d), in the fully developed mature nodules (28 dpi), LjPUB13Pro::GUS expression was restrained to the nodule-root connection zone (Fig. 2f). Collectively, these results illustrate that LjPUB13 transcript is present in symbiotically active root and young nodule tissues.

Ljpub13 mutants display growth defects and a delayed and reduced nodulation capacity

To investigate whether LjPUB13 plays a role in the interaction between L. japonicus and the M. loti symbiont, we identified three mutant lines that possess LORE1 retroelement insertions [29, 30] in LjPUB13 (Fig. 1a) and homozygous insertion lines were obtained for phenotypic analyses. The expression levels of LjPUB13 in uninoculated plants of pub13.1, pub13.2 and pub13.3 mutant lines are 4-, 4.5- and 19-fold reduced, respectively, compared to wild-type plants of the same age (Additional file 1: Figure S3).

We observed that in the absence of the M. loti symbiont, L. japonicus pub13 mutants displayed a reduced growth phenotype; the root and shoot length are significantly shorter in all pub13 mutants when compared to wild type plants of the same age (Fig. 3). In the presence of M. loti, pub13 mutants formed a significantly lower number of nodules compared to the wild type Gifu, both in absolute number (Fig. 4a) and when normalized to root length (Fig. 4b). The latter indicates that the reduced nodulation phenotype does not correlate to the shorter root phenotype and is further supported by the direct comparison of pub13.1 with the har1 hypernodulation mutant [31, 32], where in har1 very short roots (shorter than pub13.1) a high number of nodules is formed (Additional file 1: Figure S4). Notably, nodulation was not only significantly reduced but also delayed in pub13 mutants (Fig. 4c, Additional file 1: Figure S5). At 14 dpi only 15% of pub13–1 and pub13–2 plants had initiated symbiosis and had nodules compared to approximately 60% of wild type plants. This frequency in nodule appearance in the different plants of each plant line was higher at later stages of nodulation (Fig. 4c). Sections of mature nodules showed that they were fully infected and morphologically normal (Additional file 1: Figure S6). In addition, the root infection process was not evidently perturbed in the L. japonicus pub13 mutants since both wild type plants and pub13 mutants formed a similar number of infection threads in root hairs (Fig. 4d). Together, these results show that LjPUB13 is involved in plant growth and nodule organogenesis.

Fig. 3

L. japonicus pub13 mutants have a reduced growth phenotype. a L. japonicus pub13.3 plants have shorter roots and shoots compared to wild type. b Shoot length and c root length of pub13 mutants is significantly shorter compared to wild type plants. Graphs show means (+SE) of three biological replications (n = 10)

Fig. 4

Mutants of LjPUB13 have reduced nodulation. a Number of nodules per plant and b number of nodules per root cm in pub13.1 and pub13.2 mutants compared to wild type plants at 14 and 21 dpi, c frequency of plants carrying nodules, d number of infection threads per root cm in pub13.1 and pub13.2 mutants compared to wild type plants at 10 and 14 dpi. Graphs show means (+SE) of at least three biological replications (n = 10). Significant differences (P < 0.05) are indicated by asterisk. Comparisons are between wild-type and pub13 mutants

L. japonicus pub13 mutants display normal root PTI responses to flg22 treatment

Based on its gene structure, amino acid identity and phylogenetic relationship, LjPUB13 could be considered a putative ortholog of AtPUB13 in L. japonicus. In Arabidopsis, a number of phenotypes have been reported for Atpub13 mutants, which primarily indicates for a role in PTI (pathogen-associated molecular patterns (PAMP)/pattern-triggered immunity) response to flg22 treatment [26, 33]. Although the studies in Arabidopsis were focused in leaves, we searched for similar functions in L. japonicus roots as our aim was to investigate the function of LjPUB13 in root responses to interacting microbes or microbial signals.

Firstly, we investigated the involvement of LjPUB13 in the induction of reactive oxygen species (ROS) after flg22 treatment. Our analysis of wild type and pub13 mutant roots revealed that similar levels of ROS were produced by wild type and mutant L. japonicus plants after treatment with 0.5 μM flg22 (Additional file 1: Figure S7).

Secondly, we analysed the transcriptional activation of the defence marker genes LjMPK3, LjPEROXIDASE and LjPR1, together with LjFLS2 and LjBAK1 after 1-h treatment with flg22. All marker genes were induced in the treated samples (Additional file 1: Figure S8); LjMPK3, LjPEROXIDASE and LjPR1 had a 4- to 6-fold increase in transcript levels while a 2- to 3-fold increase was observed for LjFLS2 and LjBAK1. Nevertheless, the expression levels of these marker genes were comparable in wild type and pub13 mutants in both treated and non-treated samples (Additional file 1: Figure S8), suggesting that L. japonicus pub13 mutants neither express immunity-related genes constitutively nor overexpress them in the presence of flg22 as seen in Arabidopsis mutants [26].

Collectively, these results show that in L. japonicus roots the PUB13 gene is not directly involved in PTI responses induced by flg22.

Lotus PUB13 interacts with LjBAK1 but fails to ubiquitinate LjFLS2

It is known that in Arabidopsis, PUB13, phosphorylated by BAK1 in the presence of flg22, polyubiquitinates the FLS2 receptor. This results in FLS2 degradation and regulation of the PTI signalling downstream of FLS2-BAK1 [26]. The apparent absence of PUB13-dependent PTI responses in L. japonicus roots prompted us to investigate the molecular basis of this differential response in L. japonicus.

We identified the corresponding L. japonicus BAK1, and the LjBAK1cyt-GST or LjFLS2cyt-GST fusion proteins were produced in E. coli. Next, we investigated the ability of the LjPUB13-HIS protein to interact with LjBAK1 by performing an in vitro binding assay. Similar to Arabidopsis [26], we observed that LjPUB13 interacts with LjBAK1, while the LjPUB13-LjFLS2 interaction is barely detectable (Fig. 5a). LjPUB13 appears to not be able to ubiquitinate LjFLS2cyt, although we tested potential capacity in the presence of three different E2 enzymes (Additional file 1: Figure S9). These results indicate that PUB13 has different molecular capacities in L. japonicus and Arabidopsis.

Fig. 5

LjPUB13 interacts with LjBAK1, LjFLS2 and NFR5 in vitro. a In vitro binding assay of the HIS-tagged LjPUB13 with the GST-fused cytoplasmic regions of LjBAK1 and LjFLS2. LjPUB13 interacts with LjBAK1cyt and LjFLS2cyt. b In vitro binding assays of the GST-fused LjPUB13 and LjBAK1cyt with the HIS-tagged cytoplasmic regions of the Nod Factor receptors NFR1 and NFR5. LjPUB13 interacts strongly with NFR5cyt, while the LjPUB13-NFR1cyt interaction is undetectable. Sometimes a weak interaction of LjPUB13-NFR1cyt was observed but this result was not always reproducible. LjBAK1cyt also interacts with NFR5cyt, but no interaction was ever observed between LjBAK1cyt and NFR1cyt. After GST pulldown, bead-bound proteins were analyzed by immunoblotting using anti-HIS and anti-GST antibodies. These experiments were repeated five times

L. japonicus PUB13 interacts with and ubiquitinates NFR5

Since LjPUB13 is transcriptionally regulated during symbiosis and the gene is implicated in sustaining nodule organogenesis, we explored the involvement of LjPUB13 in Nod factor signalling. Thus, we tested the possibility of LjPUB13 to interact with the Nod factor receptors NFR1 and NFR5.

First, we investigated these supposed interactions in vitro. The cytoplasmic regions of NFR1 or NFR5 were expressed in E. coli as fusions to a HIS tag and were used in an in vitro binding assay together with the LjPUB13-GST fusion protein. We found that LjPUB13 interacted strongly with NFR5cyt, while the LjPUB13-NFR1cyt interaction was usually undetectable (Fig. 5b). Sometimes a weak interaction of LjPUB13-NFR1cyt was observed but this result was not always reproducible, and therefore cannot be fully considered as possible. Moreover, we tested whether L. japonicus BAK1, can interact with the receptors NFR1 and NFR5. LjPUB13 strongly interacts with LjBAK1 and it is plausible that LjPUB13 acts together with LjBAK1 as a complex. In addition, BAK1 is a well-known membrane co-receptor for many membrane receptor kinases [34]. Thus, in in vitro binding assays, GST-LjBAK1cyt was indeed found to interact with HIS-NFR5cyt but no interaction was observed between GST-LjBAK1cyt and HIS-NFR1cyt (Fig. 5b).

We verified the interaction of LjPUB13 with the Nod factor receptors, using a BiFC (Bimolecular Fluorescence Complementation) assay in Nicotiana benthamiana. Since the ARM domain of PUB13 is responsible for the specificity in protein-protein interactions of E3 ligases, we created a truncated fusion protein, where only the ARM domain of LjPUB13 (Additional file 1: Figure S2) was linked to N- or C-terminal half of YFP. We observed a strong reconstituted YFP signal on N. benthamiana leaf cells when LjPUB13ARM was co-expressed with LjBAK1, NFR1 or NFR5 (Fig. 6). Two other receptor-like kinases, LjCLAVATA2 [35] and LjLYS11 [36], were used as negative control interactions. Indeed, no YFP signal was detected when LjPUB13ARM was co-expressed with either of the two control proteins. These results show that the ARM domain of LjPUB13 recognizes and interacts with L. japonicus receptor kinases in a specific manner. In addition, the LjBAK1-NFR5 interaction was also confirmed by BiFC in planta (Fig. 6).

Fig. 6

In planta interactions of LjPUB13ARM and LjBAK1 with the Nod Factor receptors NFR1 and NFR5. YFP split into N- and C-terminal halves were fused to the ARM domain of LjPUB13 or to full length LjBAK1, NFR1 and NFR5. N. benthamiana leaves were cotransformed with different construct combinations and leaf epidermal cells were observed via confocal microscopy. The experiment was repeated four times with similar results. Bars, 50 μm

Finally, we examined whether LjPUB13 can ubiquitinate the Nod factor receptors NFR1 and NFR5 since LjPUB13 was shown to be a functional E3 ligase (Fig. 1b). The E. coli-produced HIS-NFR1cyt or HIS-NFR5cyt fusion proteins were used together with GST-LjPUB13 in an in vitro ubiquitination assay. Interestingly, LjPUB13 polyubiquitinated the cytosolic region of NFR5 and we detected high-molecular-weight proteins above the HIS-NFR5cyt (Fig. 7). The NFR5 ubiquitination by LjPUB13 was demonstrated by using two different E2 enzymes, UbcH6 and UbcH5b. On the contrary, when HIS-NFR1cyt was tested as a substrate of LjPUB13, no ubiquitination activity was observed (Fig. 7). Thus, L. japonicus PUB13 specifically ubiquitinates the Nod factor receptor NFR5.

Fig. 7

LjPUB13 specifically ubiquitinates the Nod Factor receptor NFR5. LjPUB13 was purified as GST fusion, while the cytoplasmic regions of NFR1 and NFR5 were purified as HIS fusion proteins. a, b The NFR5 ubiquitination was detected by an anti-HIS antibody, using two different E2 ligases: UbcH6 (a) and UbcH5b (b). No ubiquitination was observed when NFR1 was used as a substrate. c The presence of LjPUB13 at the ubiquitination reactions was detected by an anti-GST antibody (input). d Purified NFR1-HIS and NFR5-HIS used at the ubiquitination reactions. The experiments were repeated three times with similar results


E3 ligases were identified as important signalling components in nitrogen-fixing symbiosis in model legumes [20,21,22,23,24,25]. However, the targets of all the symbiotic E3 ligases remain to be identified and the signalling events downstream of their action remain unexplored. Here, we report that a PUB-ARM E3 ligase that possesses ubiquitination activity (Fig. 1) plays a direct role in Nod factor signalling in Lotus japonicus, through its interaction and ubiquitination of the Nod factor receptor NFR5. LjPUB13 is, thus, involved in the successful establishment of L. japonicus-rhizobium symbiosis, likely modulating the NFR5 protein levels or activity.

We found that L. japonicus PUB13 is involved in the organogenesis of nitrogen-fixing nodules. pub13 mutant alleles displayed reduced and delayed nodule organogenetic capacity (Fig. 4). On the other hand, bacterial infection inside root hairs was normally sustained (Fig. 4d). This suggests that LjPUB13 is required for the signalling events that lead to nodule organogenesis in the cortex, rather than the infection events. Furthermore, considering the low accumulation of LjPUB13 gene transcripts in mature nodules (Fig. 2) and the successful colonization of nodules by rhizobium in the pub13 mutants (Additional file 1: Figure S6), we envision that PUB13 E3 ligase plays a role in the initial stages of L. japonicus-rhizobium symbiosis establishment, rather than the later stages of nodule development and colonization.

The reduced nodulation phenotype could be attributed to the general reduced plant growth exhibited by pub13 mutants (Fig. 3). However, the number of nodules was significantly different in pub13 mutants from that of the wild-type plants also when expressed per unit of root length (Fig. 4b). This strongly suggests that a reduced nodulation is not directly related to a shorter root. The fact that the nodulation capacity is independent of the shoot/root length is further supported by studies in hypernodulation mutants, like har1 [31, 32], where an excessive number of nodules are produced by a very short plant (Additional file 1: Figure S4). Moreover, RNAi knockdown of LjnsRING E3 ligase, a nodule-specific gene involved in the early infection and mature nodule function [25], and mutation of Amsh1, a deubiquitination enzyme [37], affect plant growth in L. japonicus; this indicates that some proteins involved in (de)ubiquitination may act as nodes for plant growth and nodulation pathways. This is further supported by the expression of LjPUB13 observed in developing lateral roots at non-inoculated plants (Fig. 2c). The growth defective phenotype of L. japonicus pub13 mutants has also been reported for the Arabidopsis pub13 mutants [38].

The expression of LjPUB13 observed in nodule primordia (Fig. 2), together with the findings that LjPUB13 interacts with (Figs. 5, 6) and directly ubiquitinates the NFR5 receptor (Fig. 7), strongly supports the requirement for LjPUB13 during the early stages of nodule organogenesis in L. japonicus-rhizobia symbiosis. The defective nodulation phenotype observed in pub13 mutants suggests a positive regulatory role of LjPUB13 in rhizobial symbiosis. On the contrary, the symbiotic E3s MtPUB1 and SINA4 [22, 24] exhibit a clear negative role on the early steps of the infection process. This indicates that different roles may be performed by various E3 ligases during symbiosis.

LjPUB13 is the closest putative ortholog of Arabidopsis PUB13 in L. japonicus. In Arabidopsis, PUB13 has been implicated in FLS2-mediated flg22 signalling [26, 33] and, recently, in LYK5-mediated chitooctaose responses [27]. We show here that the defence responses downstream of flg22 perception are independent of PUB13 in the roots of this model legume. In contrast to what has been reported for Arabidopsis seedlings [26, 33], our results from ROS production, and defence gene regulation (Additional file 1: Figure S7 and S8) show that LjPUB13 does not appear to be involved in flg22-dependent defence responses in L. japonicus roots. In line with this, our in vitro assays show that LjPUB13 may not directly ubiquitinate LjFLS2 (Additional file 1: Figure S9), as has been shown for its Arabidopsis counterpart [26]. However, we cannot rule out that the involvement of LjPUB13 in plant immunity may be manifested in other parts of the plant or under conditions that have not been addressed in this study.

The differences observed in plant immunity responses in L. japonicus pub13 mutant lines compared to Arabidopsis mutants, could be attributed to a host-dependent specialisation of function for these PUB proteins. Functional differentiation of orthologous PUB genes was also found for Arabidopsis PUB17 and Brassica napus ARC1 [39]. Alternatively, a redundancy in the roles of PUB proteins in L. japonicus could mean that paralogs may be responsible for the roles that have been assigned to PUB13 in Arabidopsis. Blast analyses against L. japonicus PUB13 in Lotus base v.3.0 revealed the presence of, yet uncharacterized, proteins with similarity to PUB13. Proteins with highest similarity (appr. 50% identity) are presented in Additional file 1: Figure S1. Future studies are needed to examine possible implication of these proteins in defence responses in L. japonicus.

Based on our results, we propose a plausible mechanism where LjPUB13 acts on NFR5 post-translationally. Ubiquitination of NFR5 may lead to degradation, modulation of activity or re-localization. In any case, it is expected that NFR5 turnover is essential for efficient nodule organogenesis, and the recruitment of LjPUB13 ensures the onset and/or continuation of NFR5-mediated signalling and the successful initiation of nodule formation.

Along this line, and although NFR5 internalization from the plasma membrane has not been shown directly as yet, an association of NFR5 with the clathrin-mediated endocytosis has been proposed [40]. A clathrin protein (CHC1) was shown to interact with the Rho-like GTPase ROP6 [41], an interaction partner of NFR5 [42] in L. japonicus. Interestingly, ROP6 silencing in roots by RNAi did not affect the rhizobium entry in root hairs, but inhibited the IT growth through the root cortex, which resulted in the development of fewer nodules per plant [42].

We also show that LjPUB13 physically and specifically interacts with the cytosolic regions of LjBAK1 and NFR5 in vitro and in vivo. Interestingly, the cytosolic region of LjBAK1 was also shown to interact with the cytosolic region of NFR5 (Figs. 5 and 6), suggesting that PUB13/BAK1 may act as a complex in the establishment of L. japonicus-rhizobium symbiosis. The recruitment of LjBAK1 in this case is anticipated, considering the ubiquitous role of BAK1 in interactions of multiple membrane receptors kinases [26, 34].


Ubiquitination and deubiquitination of proteins play a major role in root nodule symbiosis [20,21,22,23,24,25, 37]. Although many symbiotic E3 ligases have been identified to date, the mechanistic mode of action in the symbiotic interactions has not been elucidated, and their targets remain to be identified.

In the present study, the L. japonicus PUB-ARM E3 ligase PUB13 was found to directly ubiquitinate the Nod factor receptor NFR5, and the delayed nodulation phenotype of the pub13 mutants reveal the implication of this gene in nodule development. Based on our results, we suggest that LjPUB13 has a role in the establishment of the L. japonicus-rhizobium symbiosis and acts as a positive regulator in nodule formation through the post-transcriptional control of NFR5.


Biological material, growth conditions and inoculation

L. japonicus ecotype Gifu B-129 was used as a wild type control. Wild type, pub13 and har1–3 [31] seeds were surface-sterilized as previously described [43] and were grown on wet filter paper for 3 to 4 days. Plants were then grown in petri dishes with solid quarter-strength B&D medium [44] on filter paper. The plants were grown in a vertical position in growth boxes, keeping the roots in the dark. Growth chamber conditions were 16-h day and 8-h night cycles at 21 °C.

For nodulation kinetics and IT counting, each petri dish was inoculated with 500 μl of a 0.02 OD600 culture of Mesorhizobium loti cv. R7A DsRed [45]. Nodule numbers were scored on at least 30 plants and ITs were counted on at least 30 cm of root.

For gene expression analyses the plants were grown in pots with 2:1 mix of sand and vermiculite. Half of the plants were inoculated with a 0.1 OD600 culture of M. loti cv. R7A. The plants were watered periodically with Hoagland nutrient solution.

Protein purification

For GST-LjPUB13, GST-LjBAK1cyt, and GST-LjFLS2cyt, the corresponding ORFs were PCR amplified and the resulting fragments were cloned between the EcoRI and XhoI sites of pGEX-6P-1 (GE Healthcare). For HIS-LjPUB13, the ORF was PCR amplified and the resulting fragment was cloned between the EcoRI and XhoI sites of pET21a (Novagen). Primers are listed in Additional file 1: Table S1. GST- and HIS- tagged proteins were expressed in E. coli Rosetta (DE3) (Merck Chemicals). GST-tagged proteins were then purified with Glutathione-Sepharose 4B beads (GE Healthcare) and eluted with 0.2 M reduced glutathione, whilst HIS-tagged proteins were purified with Ni Sepharose High Performance beads (GE Healthcare) and eluted with 250 mM imidazole.

Plasmids with the NFR1cyt and NFR5cyt fragments cloned in pProEX-1 vector were kindly provided by Mickael Blaise. HIS-NFR1cyt and HIS-NFR5cyt proteins were expressed in Rosetta 2 E. coli (DE3) competent cells (Novagen). IMAC purification was performed using Ni-NTA columns (Qiagen). The proteins were eluted with an elution buffer containing 50 mM Tris-HCl pH 8, 500 mM NaCl, 500 mM imidazole pH 8, 1 mM benzamidine, 5 mM β-mercaptoethanol and 10% glycerol. The eluted proteins were then injected onto a Superdex 200 increase 10/300 GL (GE Healthcare) column connected to an ÄKTA PURIFIER system (GE Healthcare) and eluted with a buffer containing 50 mM Tris-HCl pH 8, 500 mM NaCl, 5 mM β-mercaptoethanol and 10% glycerol.

In vitro ubiquitination assay

The in vitro ubiquitination assays were performed with an Ubiquitinylation kit (BML-UW9920, Enzo Life Sciences), using either UbcH6 or UbcH5b E2 enzymes and following the manufacturer’s protocol. For LjFLS2 ubiquitination tests UbcH5c E2 enzyme was also used.

All the reactions were incubated at 30 °C for 3 h, and then stopped by adding SDS sample buffer and boiled at 98 °C for 5 min. The samples were then separated by SDS–PAGE and analyzed by Western blotting, using the a-GST antibody (GE Healthcare) and the a-HIS antibody (Roche) to detect the tagged proteins or the anti-ubiquitin antibody (Santa Cruz Biotechnology P4D1) to detect the ubiquitinated fraction.

In vitro binding assay

For the in vitro binding assays, 2 μg of a GST tagged protein was incubated with 2 μg of a HIS tagged protein, in different protein combinations, in 150 μl of cold buffer A (50 mM Tris-HCl pH 7.5, 100 mM NaCl, 10% glycerol) with 0.1% Triton X-100, for 1.5 h at 4 °C. Glutathione-Sepharose 4B beads (GE Healthcare) were added to the mixtures and incubated with the proteins for 2 more hrs at 4 °C. The beads were then washed 3 times with cold buffer A. The bead-bound proteins were analysed by immunoblotting, according to standard protocols, using anti-GST (GE Healthcare) and anti-HIS (Roche) antibodies.

Expression analysis by qRT-PCR

To test the temporal and spatial expression of LjPUB13 in non-inoculated and M. loti inoculated L. japonicus plants and the expression of defence genes after treatment with flg22, analysis by qRT-PCR was performed as previously described [45, 46]. Gene primers are listed in Additional file 1: Table S2.

LjPUB13 promoter activity in L. japonicus

For the promoter-GUS-terminator construction, a 1475 bp promoter with a 5′ untranslated region (UTR) and a 310 bp terminator region was amplified from L. japonicus genomic DNA (the primers are listed in Additional file 1: Table S3) and cloned into GoldenGate vectors [47]. The PCR fragments were firstly cloned into GoldenGate Level 0 vectors before being assembled as a construct (promoter:GUS:terminator) in a modified pIV10 vector [48]. The construct was transferred into A. rhizogenes AGL1.

Hairy root induction using A. rhizogenes was performed as described previously [49]. Chimeric plants were transferred into magenta growth boxes containing a sterilized 4:1 mix of clay granules and vermiculite as well as ¼ strength B&D medium supplemented with 1 mM KNO3. For inoculation, liquid cultures of M. loti cv. R7A expressing DsRed were grown to an optical density of 0.02 and applied directly to the root systems (0.7 ml per plant). Plants were grown at 21 °C (16 h light, 8 h dark) and harvested at indicated times post inoculation. GUS staining was performed as described previously [50]. Whole roots were visualized on a Leica M165 FC stereomicroscope.

Infection thread formation

For inspection of infection-thread formation, roots were harvested 10 and 14 days after inoculation with M. loti strain cv. R7A DsRed. Sections (1 cm) of at least 30 plants were examined under a Zeiss Axioplan 2 fluorescent microscope.

ROS accumulation

Seven-day-old roots were cut into 0.5 cm pieces and incubated overnight (with shaking) in 200 μl water in 96-well plates (Grenier Bio-one). Before the measurements, the water was exchanged with 200 μl of buffer (20 mM luminol, Sigma; 5 μg/ml horseradish peroxidase, Sigma), supplemented with either H2O or 0.5 μM flg22 peptide (FLS22-P-1, Alpha Diagnostic). Luminescence was recorded with a Varioskan™ Flash Multimode Reader (Thermo).

Protein-protein interaction studies by bimolecular fluorescence complementation (BiFC)

The ARM domain of LjPUB13 (Additional file 1: Figure S2) and the full length LjBAK1 were cloned, using Gateway technology (Invitrogen), into the pGREEN029:35S:GW:nYFP/cYFP vectors creating N- and C-terminal fusions to YFP. The primers are listed in Additional file 1: Table S3. The NFR1 and NFR5 fused to nYFP/cYFP constructs used were the same as those described in [51].

Agrobacterium tumefaciens AGL1 cells transformed with the protein expression plasmids were grown in 5 ml LB medium supplemented with appropriate antibiotics at 28 °C. Bacteria were pelleted by centrifugation at 4000 g for 20 min at room temperature and re-suspended in agroinfiltration medium (10 mM MgCl2, 10 mM MES and 450 μM acetosyringone), incubated for 2–3 h in the dark and finally re-suspended to an OD600 of 0.2. A 1:1 mixture of cultures was prepared for each construct combination together with the P19 construct (to an OD600 of 0.02). A 1-ml syringe was used for infiltration of the bacterial mixture to the abaxial side of the Nicotiana benthamiana leaves. Fluorescence was detected using a Zeiss LSM510 confocal microscope.

Statistical analyses

Differences in the tested biological parameters between mutant and wild type plants were analyzed by Student’s t-test. A significant level of 5% was applied.





Bimolecular Fluorescence Complementation




Glutathione S-transferase




Pathogen-associated molecular patterns (PAMP)/pattern-triggered immunity


Plant U-box


Reactive oxygen species


Yellow Fluorescent Protein


  1. 1.

    Limpens E, Franken C, Smit P, Willemse J, Bisseling T, Geurts R. LysM domain receptor kinases regulating rhizobial nod factor-induced infection. Science. 2003;302:630–3.

    CAS  Article  Google Scholar 

  2. 2.

    Madsen EB, Madsen LH, Radutoiu S, Olbryt M, Rakwalska M, Szczyglowski K, Sato S, Kaneko T, Tabata S, Sandal N, Stougaard J. A receptor kinase gene of the LysM type is involved in legume perception of rhizobial signals. Nature. 2003;425:637–40.

    CAS  Article  Google Scholar 

  3. 3.

    Radutoiu S, Madsen LH, Madsen EB, Felle HH, Umehara Y, Grønlund M, Sato S, Nakamura Y, Tabata S, Sandal N, Stougaard J. Plant recognition of symbiotic bacteria requires two LysM receptor-like kinases. Nature. 2003;425:585–92.

    CAS  Article  Google Scholar 

  4. 4.

    Arrighi JF, Barre A, Ben Amor B, Bersoult A, Soriano LC, Mirabella R, de Carvalho-Niebel F, Journet EP, Ghérardi M, Huguet T, Geurts R, Dénarié J, Rougé P, Gough C. The Medicago truncatula lysine motif receptor-like kinase gene family includes NFP and new nodule-expressed genes. Plant Physiol 2006;142:265–279.

    CAS  Article  Google Scholar 

  5. 5.

    Mulder L, Lefebvre B, Cullimore J, Imberty A. LysM domains of Medicago truncatula NFP protein involved in nod factor perception. Glycosylation state, molecular modeling and docking of chitooligosaccharides and nod factors. Glycobiology. 2006;16:801–9.

    CAS  Article  Google Scholar 

  6. 6.

    Smit P, Limpens E, Geurts R, Fedorova E, Dolgikh E, Gough C, Bisseling T. Medicago LYK3, an entry receptor in rhizobial nodulation factor signalling. Plant Physiol. 2007;145:183–91.

    CAS  Article  Google Scholar 

  7. 7.

    Indrasumunar A, Kereszt A, Searle I, Miyagi M, Li D, Nguyen CD, Men A, Carroll BJ, Gresshoff PM. Inactivation of duplicated nod factor receptor 5 (NFR5) genes in recessive loss-of-function non-nodulation mutants of allotetraploid soybean (Glycine max L. Merr.). Plant Cell Physiol. 2010;51:201–14.

    CAS  Article  Google Scholar 

  8. 8.

    Indrasumunar A, Searle I, Lin MH, Kereszt A, Men A, Carroll BJ, Gresshoff PM. Nodulation factor receptor kinase 1α controls nodule organ number in soybean (Glycine max L. Merr). Plant J. 2011;65:39–50.

    CAS  Article  Google Scholar 

  9. 9.

    Oldroyd GE, Murray JD, Poole PS, Downie JA. The rules of engagement in the legume-rhizobial symbiosis. Annu Rev Genet. 2011;45:119–44.

    CAS  Article  Google Scholar 

  10. 10.

    Vierstra RD. The ubiquitin–26S proteasome system at the nexus of plant biology. Nat Rev Mol Cell Biol. 2009;10:385–97.

    CAS  Article  Google Scholar 

  11. 11.

    Komander D. The emerging complexity of protein ubiquitination. Biochem Soc Trans. 2009;37:937–53.

    CAS  Article  Google Scholar 

  12. 12.

    Hershko A, Ciechanover A. The ubiquitin system. Annu Rev Biochem. 1998;67:425–79.

    CAS  Article  Google Scholar 

  13. 13.

    Pickart CM, Eddins MJ. Ubiquitin, structures, functions, mechanisms. Biochim Biophys Acta. 2004;1695:55–72.

    CAS  Article  Google Scholar 

  14. 14.

    Smalle J, Vierstra RD. The ubiquitin 26S proteasome proteolytic pathway. Annu Rev Plant Biol. 2004;55:555–90.

    CAS  Article  Google Scholar 

  15. 15.

    Stone SL, Callis J. Ubiquitin ligases mediate growth and development by promoting protein death. Curr Opin Plant Biol. 2007;10:624–32.

    CAS  Article  Google Scholar 

  16. 16.

    Mudgil Y, Shiu SH, Stone SL, Salt JN, Goring DR. A large complement of the predicted Arabidopsis ARM repeat proteins are members of the U-box E3 ubiquitin ligase family. Plant Physiol. 2004;134:59–66.

    CAS  Article  Google Scholar 

  17. 17.

    Samuel MA, Salt JN, Shiu SH, Goring DR. Multifunctional arm repeat domains in plants. Int Rev Cell Mol Biol. 2006;253:1–26.

    CAS  Google Scholar 

  18. 18.

    Gu T, Mazzurco M, Sulaman W, Matias DD, Goring DR. Binding of an arm repeat protein to the kinase domain of the S-locus receptor kinase. Proc Natl Acad Sci U S A. 1998;95:382–7.

    CAS  Article  Google Scholar 

  19. 19.

    Hervé C, Lefebvre B, Cullimore J. How many E3 ubiquitin ligase are involved in the regulation of nodulation? Plant Signal Behav. 2011;6:660–4.

    Article  Google Scholar 

  20. 20.

    Kiss E, Oláh B, Kaló P, Morales M, Heckmann AB, Borbola A, Lózsa A, Kontár K, Middleton P, Downie JA, et al. LIN, a novel type of U-box/WD40 protein, controls early infection by rhizobia in legumes. Plant Physiol. 2009;151:1239–49.

    CAS  Article  Google Scholar 

  21. 21.

    Yano K, Shibata S, Chen WL, Sato S, Kaneko T, Jurkiewicz A, Sandal N, Banba M, Imaizumi-Anraku H, Kojima T, et al. CERBERUS, a novel U-box protein containing WD-40 repeats, is required for formation of the infection thread and nodule development in the legume-rhizobium symbiosis. Plant J. 2009;60:168–80.

    CAS  Article  Google Scholar 

  22. 22.

    Mbengue M, Camut S, de Carvalho-Niebel F, Deslandes L, Froidure S, Klaus-Heisen D, Moreau S, Rivas S, Timmers T, Hervé C, et al. The Medicago truncatula E3 ubiquitin ligase PUB1 interacts with the LYK3 symbiotic receptor and negatively regulates infection and nodulation. Plant Cell. 2010;22:3474–88.

    CAS  Article  Google Scholar 

  23. 23.

    Vernié T, Camut S, Camps C, Rembliere C, de Carvalho-Niebel F, Mbengue M, Timmers T, Gasciolli V, Thompson R, le Signor C, Lefebvre B, Cullimore J, Hervé C. PUB1 interacts with the receptor kinase DMI2 and negatively regulates rhizobial and arbuscular mycorrhizal symbioses through its ubiquitination activity in Medicago truncatula. Plant Physiol. 2016;170:2312–24.

    Article  Google Scholar 

  24. 24.

    Den Herder G, Yoshida S, Antolín-Llovera M, Ried MK, Parniske M. Lotus japonicus E3 ligase SEVEN IN ABSENTIA4 destabilizes the symbiosis receptor-like kinase SYMRK and negatively regulates rhizobial infection. Plant Cell. 2012;24:1691–707.

    CAS  Article  Google Scholar 

  25. 25.

    Shimomura K, Nomura M, Tajima S, Kouchi H. LjnsRING, a novel RING finger protein, is required for symbiotic interactions between Mesorhizobium loti and Lotus japonicus. Plant Cell Physiol. 2006;47:1572–81.

    CAS  Article  Google Scholar 

  26. 26.

    Lu D, Lin W, Gao X, Wu S, Cheng C, Avila J, Heese A, Devarenne TP, He P, Shan L. Direct ubiquitination of pattern recognition receptor FLS2 attenuates plant innate immunity. Science. 2011;332:1439–42.

    CAS  Article  Google Scholar 

  27. 27.

    Liao D, Cao Y, Sun X, Espinoza C, Nguyen C, Liang Y, Stacey G. Arabidopsis E3 ubiquitin ligase PLANT U-BOX13 (PUB13) regulates chitin receptor LYSIN MOTIF RECEPTOR KINASE5 (LYK5) protein abundance. New Phytol. 2017;214:1646–56.

    CAS  Article  Google Scholar 

  28. 28.

    Li W, Ahn IP, Ning Y, Park CH, Zeng L, Whitehill JGA, Lu H, Zhao Q, Ding B, Xie Q, et al. The U-box/ARM E3 ligase PUB13 regulates cell death, defense, and flowering time in Arabidopsis. Plant Physiol. 2012;159:239–50.

    CAS  Article  Google Scholar 

  29. 29.

    Fukai E, Soyano T, Umehara Y, Nakayama S, Hirakawa H, Tabata S, Sato S, Hayashi M. Establishment of a Lotus japonicus gene tagging population using the exon-targeting endogenous retrotransposon LORE1. Plant J. 2012;69:720–30.

    CAS  Article  Google Scholar 

  30. 30.

    Urbański DF, Małolepszy A, Stougaard J, Andersen SU. Genome-wide LORE1 retrotransposon mutagenesis and high-throughput insertion detection in Lotus japonicus. Plant J. 2012;69:731–41.

    Article  Google Scholar 

  31. 31.

    Krusell L, Madsen LH, Sato S, Aubert G, Genua A, Szczyglowski K, Duc G, Kaneko T, Tabata S, de Bruijn F, Pajuelo E, Sandal N, Stougaard J. Shoot control of root development and nodulation is mediated by a receptor-like kinase. Nature. 2002;420:422–6.

    CAS  Article  Google Scholar 

  32. 32.

    Nishimura R, Hayashi M, Wu GJ, Kouchi H, Imaizumi-Anraku H, Murakami Y, Kawasaki S, Akao S, Ohmori M, Nagasawa M, Harada K, Kawaguchi M. HAR1 mediates systemic regulation of symbiotic organ development. Nature. 2002;420:426–9.

    CAS  Article  Google Scholar 

  33. 33.

    Zhou J, Lu D, Xu G, Finlayson SA, He P, Shan L. The dominant negative ARM domain uncovers multiple functions of PUB13 in Arabidopsis immunity, flowering, and senescence. J Exp Bot. 2015;66:3353–66.

    CAS  Article  Google Scholar 

  34. 34.

    Chinchilla D, Shan L, He P, de Vries S, Kemmerling B. One for all, the receptor-associated kinase BAK1. Trends Plant Sci. 2009;14:535–41.

    CAS  Article  Google Scholar 

  35. 35.

    Krusell L, Sato N, Fukuhara I, Koch BE, Grossmann C, Okamoto S, Oka-Kira E, Otsubo Y, Aubert G, Nakagawa T, Sato S, Tabata S, Duc G, Parniske M, Wang TL, Kawaguchi M, Stougaard J. The Clavata2 genes of pea and Lotus japonicus affect autoregulation of nodulation. Plant J. 2011;65:861–71.

    CAS  Article  Google Scholar 

  36. 36.

    Rasmussen SR, Füchtbauer W, Novero M, Volpe V, Malkov N, Genre A, Bonfante P, Stougaard J, Radutoiu S. Intraradical colonization by arbuscular mycorrhizal fungi triggers induction of a lipochitooligosaccharide receptor. Sci Rep. 2016;6:29733.

    CAS  Article  Google Scholar 

  37. 37.

    Małolepszy A, Urbański DF, James EK, Sandal N, Isono E, Stougaard J, Andersen SU. The deubiquitinating enzyme AMSH1 is required for rhizobial infection and nodule organogenetis in Lotus japonicus. Plant J. 2015;83:719–31.

    Article  Google Scholar 

  38. 38.

    Antignani V, Klocko AL, Bak G, Chandrasekaran SD, Dunivin T, Nielsen E. Recruitment of PLANT U-BOX13 and the PI4Kb1/b2 phosphatidylinositol-4 kinases by the small GTPase RabA4B plays important roles during salicylic acid-mediated plant defense signalling in Arabidopsis. Plant Cell. 2015;27:243–61.

    CAS  Article  Google Scholar 

  39. 39.

    Yang CW, González-Lamothe R, Ewan RA, Rowland O, Yoshioka H, Shenton M, Ye H, O'Donnell E, Jones JD, Sadanandom A. The E3 ubiquitin ligase activity of Arabidopsis PLANT U-BOX17 and its functional tobacco homolog ACRE276 are required for cell death and defense. Plant Cell. 2006;18:1084–98.

    CAS  Article  Google Scholar 

  40. 40.

    Wang C, Xu X, Hong Z, Feng Y, Zhang Z. Involvement of ROP6 and clathrin in nodulation factor signalling. Plant Signal Behav. 2015b;10:e1033127.

    PubMed  PubMed Central  Google Scholar 

  41. 41.

    Wang C, Zhu M, Duan L, Yu H, Chang X, Li L, Kang H, Feng Y, Zhu H, Hong Z, Zhang Z. Lotus japonicus clathrin heavy chain 1 is associated with rho-like GTPase ROP6 and involved in nodule formation. Plant Physiol. 2015a;167:1497–510.

    CAS  Article  Google Scholar 

  42. 42.

    Ke D, Fang Q, Chen C, Zhu H, Chen T, Chang X, Yuan S, Kang H, Ma L, Hong Z, Zhang Z. The small GTPase ROP6 interacts with NFR5 and is involved in nodule formation in Lotus japonicus. Plant Physiol. 2012;159:131–43.

    CAS  Article  Google Scholar 

  43. 43.

    Handberg K, Stougaard J. Lotus japonicus, an autogamous, diploid legume species for classical and molecular genetics. Plant J. 1992;2:487–96.

    Article  Google Scholar 

  44. 44.

    Broughton WJ, Dilworth MJ. Control of leghaemoglobin synthesis in snake beans. Biochem J. 1971;125:1075–80.

    CAS  Article  Google Scholar 

  45. 45.

    Kawaharada Y, Kelly S, Nielsen MW, Hjuler CT, Gysel K, Muszyński A, Carlson RW, Thygesen MB, Sandal N, Asmussen MH, et al. Receptor-mediated exopolysaccharide perception controls bacterial infection. Nature. 2015;523:308–12.

    CAS  Article  Google Scholar 

  46. 46.

    Tanou G, Minas IS, Karagiannis E, Tsikou D, Audebert S, Papadopoulou KK, Molassiotis A. The impact of sodium nitroprusside and ozone in kiwifruit ripening physiology, a combined gene and protein expression profiling approach. Ann Bot. 2015;116:649–62.

    CAS  Article  Google Scholar 

  47. 47.

    Weber E, Engler C, Gruetzner R, Werner S, Marillonnet S. A modular cloning system for standardized assembly of multigene constructs. PLoS One. 2011;6:e16765.

    CAS  Article  Google Scholar 

  48. 48.

    Stougaard J, Abildsten D, Marcker KA. The Agrobacterium rhizogenes pRi TL-DNA segment as a gene vector system for transformation of plants. Mol Gen Genomics. 1987;207:251–5.

    CAS  Article  Google Scholar 

  49. 49.

    Hansen J, Jørgensen JE, Stougaard J, Marcker K. Hairy roots - a short cut to transgenic root nodules. Plant Cell Rep. 1989;8:12–5.

    CAS  Article  Google Scholar 

  50. 50.

    Vitha S, Benes K, Phillips JP, Gartland KM. Histochemical GUS analysis. Methods Mol Biol. 1995;44:185–93.

    CAS  PubMed  Google Scholar 

  51. 51.

    Madsen EB, Antolín-Llovera M, Grossmann C, Ye J, Vieweg S, Broghammer A, Krusell L, Radutoiu S, Jensen ON, Stougaard J, Parniske M. Autophosphorylation is essential for the in vivo function of the Lotus japonicus nod factor receptor 1 and receptor-mediated signalling in cooperation with nod factor receptor 5. Plant J. 2011;65:404–17.

    CAS  Article  Google Scholar 

Download references


The authors wish to thank Prof. Jens Stougaard (Aarhus University) for scientific support and for making corrections and comments on the manuscript; and Prof. Claus Schwechheimer (Technische Universität München) for hosting DT in his laboratory by providing financial support through the SFB924 and for making comments on the manuscript. The authors thank Mickael Blaise for providing the plasmids for HIS-NFR1cyt and HIS-NFR5cyt protein expression, Zoltán Bozsóki for helping with the ROS assay, and Finn Pedersen for taking care of the plants in the greenhouse.


This work was partially supported by the Postgraduate Programs 3817 & 3439 of the Department of Biochemistry and Biotechnology, University of Thessaly (to KKP). DT and SR were supported by the Danish National Research Foundation grant no. DNRF79. DT was supported by a STSM (290713–031130) from Cost Action FA1103.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.

The LjPUB13 sequence is available in GenBank under the accession number KY131979 and in Lotus base v.3.0 ( as Lj3g3v3189730.1. Other sequence data from this article can be found in the GenBank under the following accession numbers: LjBAK1 (KY131980), LjFLS2 (JN099749), NFR1 (AJ575249), NFR5 (AJ575255).

Author information




KKP, SR and DT conceived the research plan; DT, EER, ISP, JEW and DBJ conducted experiments; KKP, SR and EI supervised the experiments; EI provided training to DT in biochemistry; DT, SR and KKP analysed data; DT wrote the article with contributions of all the authors; SR and KKP complemented and revised the writing. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Simona Radutoiu or Kalliope K. Papadopoulou.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. Primers used for cloning into expression vectors. Table S2. Primers used in qRT-PCR. Table S3. Primers used for cloning into GoldenGate and Gateway vectors. Figure S1. Phylogenetic tree of amino acid sequences of LjPUB13 with previously characterized E3s of other species and L. japonicus uncharacterized PUBs. Figure S2. Amino acid sequence alignment of LjPUB13 with AtPUB13. Figure S3. Expression levels of LjPUB13 in pub13 LORE1 mutants. Figure S4. Formation of nodules in pub13.1 vs. har1 mutants. Figure S5. Inoculated 28-day-old wild type and homozygous pub13.3 mutants. Figure S6. Nodule sections of wild type, pub13.1 and pub13.2 plants. Figure S7. ROS accumulation in roots of pub13 mutants vs wt plants. Figure S8. Expression of defence genes in L. japonicus wild type and pub13 mutants. Figure S9. LjFLS2 ubiquitination tests. (DOCX 39900 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tsikou, D., Ramirez, E.E., Psarrakou, I.S. et al. A Lotus japonicus E3 ligase interacts with the Nod Factor Receptor 5 and positively regulates nodulation. BMC Plant Biol 18, 217 (2018).

Download citation


  • E3 ligase
  • Nodulation
  • Lotus japonicus
  • PUB13
  • Ubiquitination
  • Symbiosis