Involvement of mitogen-activated protein kinases in the symbiosis Bradyrhizobium-Lupinus

Changes in the m6A RNA methylome accompany the promotion of soybean root growth by rhizobia under cadmium stress

Cadmium (Cd) is the most widely distributed heavy metal pollutant in soil and has significant negative effects on crop yields and human health. Rhizobia can enhance soybean growth in the presence of heavy metals, and the legume-rhizobia symbiosis has been used to promote heavy-metal phytoremediation, but much remains to be learned about the molecular networks that underlie these effects. Here, we demonstrated that soybean root growth was strongly suppressed after seven days of Cd exposure but that the presence of rhizobia largely eliminated this effect, even prior to nodule development. Moreover, rhizobia did not appear to promote root growth by limiting plant Cd uptake: seedlings with and without rhizobia had similar root Cd concentrations. Previous studies have demonstrated a role for m6A RNA methylation in the response of rice and barley to Cd stress. We therefore performed transcriptome-wide m6A methylation profiling to investigate changes in the soybean RNA methylome in response to Cd with and without rhizobia. Here, we provide some of the first data on transcriptome-wide m6a RNA methylation patterns in soybean; m6A modifications were concentrated at the 3' UTR of transcripts and showed a positive relationship with transcript abundance. Transcriptome-wide m6A RNA methylation peaks increased in the presence of Cd, and the integration of m6A methylome and transcriptome results enabled us to identify 154 genes whose transcripts were both differentially methylated and differentially expressed in response to Cd stress. Annotation results suggested that these genes were associated with Ca²⁺ homeostasis, ROS pathways, polyamine metabolism, MAPK signaling, hormones, and biotic stress responses. There were 176 differentially methylated and expressed transcripts under Cd stress in the presence of rhizobia. In contrast to the Cd-only gene set, they were also enriched in genes related to auxin, jasmonic acid, and brassinosteroids, as well as abiotic stress tolerance. They contained fewer genes related to Ca²⁺ homeostasis and also included candidates with known functions in the legume-rhizobia symbiosis. These findings offer new insights into how rhizobia promote soybean root growth under Cd stress; they provide candidate genes for research on plant heavy metal responses and for the use of legumes in phytoremediation.

Glycine max Homologs of DOESN'T MAKE INFECTIONS 1, 2, and 3 Function to Impair Heterodera glycines Parasitism While Also Regulating Mitogen Activated Protein Kinase Expression

Glycine max root cells developing into syncytia through the parasitic activities of the pathogenic nematode Heterodera glycines underwent isolation by laser microdissection (LM). Microarray analyses have identified the expression of a G. max DOESN'T MAKE INFECTIONS3 (DMI3) homolog in syncytia undergoing parasitism but during a defense response. DMI3 encodes part of the common symbiosis pathway (CSP) involving DMI1, DMI2, and other CSP genes. The identified DMI gene expression, and symbiosis role, suggests the possible existence of commonalities between symbiosis and defense. G. max has 3 DMI1, 12 DMI2, and 2 DMI3 paralogs. LM-assisted gene expression experiments of isolated syncytia under further examination here show G. max DMI1-3, DMI2-7, and DMI3-2 expression occurring during the defense response in the H. glycines-resistant genotypes G.max _{[Peking/PI548402]} and G.max _[PI88788] indicating a broad and consistent level of expression of the genes. Transgenic overexpression (OE) of G. max DMI1-3, DMI2-7, and DMI3-2 impairs H. glycines parasitism. RNA interference (RNAi) of G. max DMI1-3, DMI2-7, and DMI3-2 increases H. glycines parasitism. The combined opposite outcomes reveal a defense function for these genes. Prior functional transgenic analyses of the 32-member G. max mitogen activated protein kinase (MAPK) gene family has determined that 9 of them act in the defense response to H. glycines parasitism, referred to as defense MAPKs. RNA-seq analyses of root RNA isolated from the 9 G. max defense MAPKs undergoing OE or RNAi reveal they alter the relative transcript abundances (RTAs) of specific DMI1, DMI2, and DMI3 paralogs. In contrast, transgenically-manipulated DMI1-3, DMI2-7, and DMI3-2 expression influences MAPK3-1 and MAPK3-2 RTAs under certain circumstances. The results show G. max homologs of the CSP, and defense pathway are linked, apparently involving co-regulated gene expression.

QTL Mapping and Data Mining to Identify Genes Associated With the Sinorhizobium fredii HH103 T3SS Effector NopD in Soybean

In some legume-rhizobium symbioses, host specificity is influenced by rhizobial type III effectors-nodulation outer proteins (Nops). However, the genes encoding host proteins that interact with Nops remain unknown. In this study, we aimed to identify candidate soybean genes associated with NopD, one of the type III effectors of Sinorhizobium fredii HH103. The results showed that the expression pattern of NopD was analyzed in rhizobia induced by genistein. We also found NopD can be induced by TtsI, and NopD as a toxic effector can induce tobacco leaf death. In 10 soybean germplasms, NopD played a positively effect on nodule number (NN) and nodule dry weight (NDW) in nine germplasms, but not in Kenjian28. Significant phenotype of NN and NDW were identified between Dongnong594 and Charleston, Suinong14 and ZYD00006, respectively. To map the quantitative trait locus (QTL) associated with NopD, a recombinant inbred line (RIL) population derived from the cross between Dongnong594 and Charleston, and chromosome segment substitution lines (CSSLs) derived from Suinong14 and ZYD00006 were used. Two overlapping conditional QTL associated with NopD on chromosome 19 were identified. Two candidate genes were identified in the confident region of QTL, we found that NopD could influence the expression of Glyma.19g068600 (FBD/LRR) and expression of Glyma.19g069200 (PP2C) after HH103 infection. Haplotype analysis showed that different types of Glyma.19g069200 haplotypes could cause significant nodule phenotypic differences, but Glyma.19g068600 (FBD/LRR) was not. These results suggest that NopD promotes S. fredii HH103 infection via directly or indirectly regulating Glyma.19g068600 and Glyma.19g069200 expression during the establishment of symbiosis between rhizobia and soybean plants.

Connection the Rhizomicrobiome and Plant MAPK Gene Expression Response to Pathogenic Fusarium oxysporum in Wild and Cultivated Soybean

Little known the connections between soybeans mitogen-activated protein kinase (MAPK) gene expression and the rhizomicrobiome upon invasion of the root pathogen Fusarium oxysporum. To address this lack of knowledge, we assessed the rhizomicrobiome and root transcriptome sequencing of wild and cultivated soybean during the invasion of F. oxysporum. Results indicated F. oxysporum infection enriched Bradyrhizobium spp. and Glomus spp. and induced the expression of more MAPKs in the wild soybean than cultivated soybean. MAPK gene expression was positively correlated with Pseudomonadaceae but negatively correlated with Sphingomonadaceae and Glomeraceae in both cultivated and wild soybean. Specifically, correlation profiles revealed that Pseudomonadaceae was especially correlated with the induced expression of GmMAKKK13-2 (Glyma.14G195300) and GmMAPK3-2 (Glyma.12G073000) in wild and cultivated soybean during F. oxysporum invasion. Main fungal group Glomeraceae was positively correlated with GmMAPKKK14-1 (Glyma.18G060900) and negatively correlated with GmRaf6-4 (Glyma.02G215300) in the wild soybean response to pathogen infection; while there were positive correlations between Hypocreaceae and GmMAPK3-2 (Glyma.12G073000) and between Glomeraceae and GmRaf49-3 (Glyma.06G055300) in the wild soybean response, these correlations were strongly negative in the response of cultivated soybean to F. oxysporum. Taken together, MAPKs correlated with different rhizomicrobiomes indicating the host plant modulated by the host self-immune systems in response to F. oxysporum.

The p38-like MAP kinase modulated H2O2 accumulation in wounding signaling pathways of sweet potato

In sweet potato (Ipomoea batatas cv Tainung 57), MAPK cascades are involved in the regulation of Ipomoelin (IPO) expression upon wounding. p38 MAPK plays an important role in plant's responses to various environmental stresses. However, the role of p38-like MAPK in wounding response is still unknown. In this study, the levels of phosphorylated-p38-like MAPK (pp38-like MAPK) in sweet potato were noticeably reduced after wounding. In addition, SB203580 (SB), a specific inhibitor blocking p38 MAPK phosphorylation, considerably decreased the accumulation of pp38-like MAPK. Expression of a wound-inducible gene IPO was elevated by SB. Moreover, it stimulated hydrogen peroxide (H₂O₂) production rather than cytosolic Ca²⁺ elevation in sweet potato leaves. However, NADPH oxidase (NOX) inhibitor diphenyleneiodonium could not inhibit IPO induction stimulated by SB. These results indicated a p38-like MAPK mechanism was involved in the regulation of IPO expression through NOX-independent H₂O₂ generation. In addition, the presence of the protein phosphatase inhibitor okadaic acid or the MEK1/ERK inhibitor PD98059 repressed the H₂O₂- or SB-induced IPO expression, demonstrating phosphatase(s) and MEK1/ERK functioning in the downstream of H₂O₂ and pp38-like MAPK in the signal transduction pathway stimulating IPO. Conclusively, wounding decreased the amount of pp38-like MAPK, stimulated H₂O₂ production, and then induced IPO expression.

Identification of Soybean Genes Whose Expression is Affected by the Ensifer fredii HH103 Effector Protein NopP

In some legume⁻rhizobium symbioses, host specificity is influenced by rhizobial nodulation outer proteins (Nops). However, the genes encoding host proteins that interact with Nops remain unknown. We generated an Ensifer fredii HH103 NopP mutant (HH103ΩNopP), and analyzed the nodule number (NN) and nodule dry weight (NDW) of 10 soybean germplasms inoculated with the wild-type E. fredii HH103 or the mutant strain. An analysis of recombinant inbred lines (RILs) revealed the quantitative trait loci (QTLs) associated with NopP interactions. A soybean genomic region containing two overlapping QTLs was analyzed in greater detail. A transcriptome analysis and qRT-PCR assay were used to identify candidate genes encoding proteins that interact with NopP. In some germplasms, NopP positively and negatively affected the NN and NDW, while NopP had different effects on NN and NDW in other germplasms. The QTL region in chromosome 12 was further analyzed. The expression patterns of candidate genes Glyma.12g031200 and Glyma.12g073000 were determined by qRT-PCR, and were confirmed to be influenced by NopP.

Cell and Developmental Biology of Plant Mitogen-Activated Protein Kinases

Plant mitogen-activated protein kinases (MAPKs) constitute a network of signaling cascades responsible for transducing extracellular stimuli and decoding them to dedicated cellular and developmental responses that shape the plant body. Over the last decade, we have accumulated information about how MAPK modules control the development of reproductive tissues and gametes and the embryogenic and postembryonic development of vegetative organs such as roots, root nodules, shoots, and leaves. Of key importance to understanding how MAPKs participate in developmental and environmental signaling is the characterization of their subcellular localization, their interactions with upstream signal perception mechanisms, and the means by which they target their substrates. In this review, we summarize the roles of MAPK signaling in the regulation of key plant developmental processes, and we survey what is known about the mechanisms guiding the subcellular compartmentalization of MAPK modules.

Differential Signaling and Sugar Exchanges in Response to Avirulent Pathogen- and Symbiont-Derived Molecules in Tobacco Cells

Plants interact with microbes whose ultimate aim is to exploit plant carbohydrates for their reproduction. Plant–microbe interactions (PMIs) are classified according to the nature of their trophic exchanges: while mutualistic microbes trade nutrients with plants, pathogens unilaterally divert carbohydrates. The early responses following microbe recognition and the subsequent control of plant sugar distribution are still poorly understood. To further decipher PMI functionality, we used tobacco cells treated with microbial molecules mimicking pathogenic or mutualistic PMIs, namely cryptogein, a defense elicitor, and chitotetrasaccharide (CO4), which is secreted by mycorrhizal fungi. CO4 was perceived by tobacco cells and triggered widespread transient signaling components such as a sharp cytosolic Ca²⁺ elevation, NtrbohD-dependent H₂O₂ production, and MAP kinase activation. These CO4-induced events differed from those induced by cryptogein, i.e., sustained events leading to cell death. Furthermore, cryptogein treatment inhibited glucose and sucrose uptake but not fructose uptake, and promoted the expression of NtSUT and NtSWEET sugar transporters, whereas CO4 had no effect on sugar uptake and only a slight effect on NtSWEET2B expression. Our results suggest that microbial molecules induce different signaling responses that reflect microbial lifestyle and the subsequent outcome of the interaction.

Phylogenomic analysis of MKKs and MAPKs from 16 legumes and detection of interacting pairs in chickpea divulge MAPK signalling modules

The mitogen-activated protein kinase (MAPK)-mediated phosphorylation cascade is a vital component of plant cellular signalling. Despite this, MAPK signalling cascade is less characterized in crop legumes. To fill this void, we present here a comprehensive phylogeny of MAPK kinases (MKKs) and MAPKs identified from 16 legume species belonging to genistoid (Lupinus angustifolius), dalbergioid (Arachis spp.), phaseoloid (Glycine max, Cajanus cajan, Phaseolus vulgaris, and Vigna spp.), and galegoid (Cicer arietinum, Lotus japonicus, Medicago truncatula, Pisum sativum, Trifolium spp., and Vicia faba) clades. Using the genes of the diploid crop chickpea (C. arietinum), an exhaustive interaction analysis was performed between MKKs and MAPKs by split-ubiquitin based yeast two-hybrid (Y2H). Twenty seven interactions of varying strengths were identified between chickpea MKKs and MAPKs. These interactions were verified in planta by bimolecular fluorescence complementation (BiFC). As a first report in plants, four intra-molecular interactions of weak strength were identified within chickpea MKKs. Additionally; two TEOSINTE-BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factors of class I were identified as novel down-stream interacting partners of seven MAPKs. We propose that this highly reliable MAPK interaction network, presented here for chickpea, can be utilized as a reference for legumes and thus will help in deciphering their role in legume-specific events.

Immunolocalization of dually phosphorylated MAPKs in dividing root meristem cells of Vicia faba, Pisum sativum, Lupinus luteus and Lycopersicon esculentum

KEY MESSAGE

In plants, phosphorylated MAPKs display constitutive nuclear localization; however, not all studied plant species show co-localization of activated MAPKs to mitotic microtubules. The mitogen-activated protein kinase (MAPK) signaling pathway is involved not only in the cellular response to biotic and abiotic stress but also in the regulation of cell cycle and plant development. The role of MAPKs in the formation of a mitotic spindle has been widely studied and the MAPK signaling pathway was found to be indispensable for the unperturbed course of cell division. Here we show cellular localization of activated MAPKs (dually phosphorylated at their TXY motifs) in both interphase and mitotic root meristem cells of Lupinus luteus, Pisum sativum, Vicia faba (Fabaceae) and Lycopersicon esculentum (Solanaceae). Nuclear localization of activated MAPKs has been found in all species. Co-localization of these kinases to mitotic microtubules was most evident in L. esculentum, while only about 50% of mitotic cells in the root meristems of P. sativum and V. faba displayed activated MAPKs localized to microtubules during mitosis. Unexpectedly, no evident immunofluorescence signals at spindle microtubules and phragmoplast were noted in L. luteus. Considering immunocytochemical analyses and studies on the impact of FR180204 (an inhibitor of animal ERK1/2) on mitotic cells, we hypothesize that MAPKs may not play prominent role in the regulation of microtubule dynamics in all plant species.

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus, Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.

Soybean MAPK, GMK1 is dually regulated by phosphatidic acid and hydrogen peroxide and translocated to nucleus during salt stress

Mitogen-activated protein kinase (MAPK) is activated by various biotic and abiotic stresses. Salt stress induces two well-characterized MAPK activating signaling molecules, phosphatidic acid (PA) via phospholipase D and phospholipase C, and reactive oxygen species (ROS) via nicotinamide adenine dinucleotide phosphate (NADPH)-oxidase. In our previous study, the activity of soybean MAPK, GMK1 was strongly induced within 5 min of 300 mM NaCl treatment and this early activity was regulated by PA. In this study, we focused on the regulation of GMK1 at the later stage of the salt stress, because its activity was strongly persistent for up to 30 min. H(2)O(2) activated GMK1 even in the presence of PA generation inhibitors, but GMK1 activity was greatly decreased in the presence of diphenyleneiodonium, an inhibitor of NADPH-oxidase after 5 min of the treatment. On the contrary, the n-butanol and neomycin reduced GMK1 activity within 5 min of the treatment. Thus, GMK1 activity may be sustained by H(2)O(2) 10 min after the treatment. Further, GMK1 was translocated into the nucleus 60 min after NaCl treatment. In the relationship between GMK1 and ROS generation, ROS generation was reduced by SB202190, a MAPK inhibitor, but was increased in protoplast overexpressing TESD-GMKK1. However, these effects were occurred at prolonged time of NaCl treatment. These data suggest that GMK1 indirectly regulates ROS generation. Taken together, we propose that soybean GMK1 is dually regulated by PA and H(2)O(2) at a time dependant manner and translocated to the nucleus by the salt stress signal.

Functional analysis of the type 3 effector nodulation outer protein L (NopL) from Rhizobium sp. NGR234: symbiotic effects, phosphorylation, and interference with mitogen-activated protein kinase signaling

Pathogenic bacteria use type 3 secretion systems to deliver virulence factors (type 3 effector proteins) directly into eukaryotic host cells. Similarly, type 3 effectors of certain nitrogen-fixing rhizobial strains affect nodule formation in the symbiosis with host legumes. Nodulation outer protein L (NopL) of Rhizobium sp. strain NGR234 is a Rhizobium-specific type 3 effector. Nodulation tests and microscopic analysis showed that distinct necrotic areas were rapidly formed in ineffective nodules of Phaseolus vulgaris (cv. Tendergreen) induced by strain NGRΩnopL (NGR234 mutated in nopL), indicating that NopL antagonized nodule senescence. Further experiments revealed that NopL interfered with mitogen-activated protein kinase (MAPK) signaling in yeast and plant cells (Nicotiana tabacum). Expression of nopL in yeast disrupted the mating pheromone (α-factor) response pathway, whereas nopL expression in N. tabacum suppressed cell death induced either by overexpression of the MAPK gene SIPK (salicylic acid-induced protein kinase) or by SIPK(DD) (mutation in the TXY motif resulting in constitutive MAPK activity). These data indicate that NopL impaired function of MAPK proteins or MAPK substrates. Furthermore, we demonstrate that NopL was multiply phosphorylated either in yeast or N. tabacum cells that expressed nopL. Four phosphorylated serines were confirmed by mass spectrometry. All four phosphorylation sites exhibit a Ser-Pro pattern, a typical motif in MAPK substrates. Taken together, data suggest that NopL mimics a MAPK substrate and that NopL suppresses premature nodule senescence by impairing MAPK signaling in host cells.

AM fungal exudates activate MAP kinases in plant cells in dependence from cytosolic Ca(2+) increase

The molecular dialogue occurring prior to direct contact between the fungal and plant partners of arbuscular-mycorrhizal (AM) symbioses begins with the release of fungal elicitors, so far only partially identified chemically, which can activate specific signaling pathways in the host plant. We show here that the activation of MAPK is also induced by exudates of germinating spores of Gigaspora margarita in cultured cells of the non-leguminous species tobacco (Nicotiana tabacum), as well as in those of the model legume Lotus japonicus. MAPK activity peaked about 15 min after the exposure of the host cells to the fungal exudates (FE). FE were also responsible for a rapid and transient increase in free cytosolic Ca(2+) in Nicotiana plumbaginifolia and tobacco cells, and pre-treatment with a Ca(2+)-channel blocker (La(3+)) showed that in these cells, MAPK activation was dependent on the cytosolic Ca(2+) increase. A partial dependence of MAPK activity on the common Sym pathway could be demonstrated for a cell line of L. japonicus defective for LjSym4 and hence unable to establish an AM symbiosis. Our results show that MAPK activation is triggered by an FE-induced cytosolic Ca(2+) transient, and that a Sym genetic determinant acts to modulate the intensity and duration of this activity.

MAPK machinery in plants: recognition and response to different stresses through multiple signal transduction pathways

The mitogen-activated protein kinase (MAPK) cascades play diverse roles in intra- and extra-cellular signaling in plants. MAP kinases are the component of kinase modules which transfer information from sensors to responses in eukaryotes including plants. They play a pivotal role in transduction of diverse extracellular stimuli such as biotic and abiotic stresses as well as a range of developmental responses including differentiation, proliferation and death. Several cascades are induced by different biotic and abiotic stress stimuli such as pathogen infections, heavy metal, wounding, high and low temperatures, high salinity, UV radiation, ozone, reactive oxygen species, drought and high or low osmolarity. MAPK signaling has been implicated in biotic stresses and has also been associated with hormonal responses. The cascade is regulated by various mechanisms, including not only transcriptional and translational regulation but through post-transcriptional regulation such as protein-protein interactions. Recent detailed analysis of certain specific MAP kinase pathways have revealed the specificity of the kinases in the cascade, signal transduction patterns, identity of pathway targets and the complexity of the cascade. The latest insights and finding are discussed in this paper in relation to the role of MAPK pathway modules in plant stress signaling.

A cell wall extract from the endophytic fungus Piriformospora indica promotes growth of Arabidopsis seedlings and induces intracellular calcium elevation in roots

Calcium (Ca2+), as a second messenger, is crucial for signal transduction processes during many biotic interactions. We demonstrate that cellular [Ca2+] elevations are early events in the interaction between the plant growth-promoting fungus Piriformospora indica and Arabidopsis thaliana. A cell wall extract (CWE) from the fungus promotes the growth of wild-type seedlings but not of seedlings from P. indica-insensitive mutants. The extract and the fungus also induce a similar set of genes in Arabidopsis roots, among them genes with Ca2+ signalling-related functions. The CWE induces a transient cytosolic Ca2+ ([Ca2+](cyt)) elevation in the roots of Arabidopsis and tobacco (Nicotiana tabacum) plants, as well as in BY-2 suspension cultures expressing the Ca2+ bioluminescent indicator aequorin. Nuclear Ca2+ transients were also observed in tobacco BY-2 cells. The Ca2+ response was more pronounced in roots than in shoots and involved Ca2+ uptake from the extracellular space as revealed by inhibitor studies. Inhibition of the Ca2+ response by staurosporine and the refractory nature of the Ca2+ elevation suggest that a receptor may be involved. The CWE does not stimulate H2O2 production and the activation of defence gene expression, although it led to phosphorylation of mitogen-activated protein kinases (MAPKs) in a Ca2+-dependent manner. The involvement of MAPK6 in the mutualistic interaction was shown for an mpk6 line, which did not respond to P. indica. Thus, Ca2+ is likely to be an early signalling component in the mutualistic interaction between P. indica and Arabidopsis or tobacco.