Identification on mitogen-activated protein kinase signaling cascades by integrating protein interaction with transcriptional profiling analysis in cotton

Genome-wide identification and analysis of abiotic stress responsiveness of the mitogen-activated protein kinase gene family in Medicago sativa L

BACKGROUND

The mitogen-activated protein kinase (MAPK) cascade is crucial cell signal transduction mechanism that plays an important role in plant growth and development, metabolism, and stress responses. The MAPK cascade includes three protein kinases, MAPK, MAPKK, and MAPKKK. The three protein kinases mediate signaling to downstream response molecules by sequential phosphorylation. The MAPK gene family has been identified and analyzed in many plants, however it has not been investigated in alfalfa.

RESULTS

In this study, Medicago sativa MAPK genes (referred to as MsMAPKs) were identified in the tetraploid alfalfa genome. Eighty MsMAPKs were divided into four groups, with eight in group A, 21 in group B, 21 in group C and 30 in group D. Analysis of the basic structures of the MsMAPKs revealed presence of a conserved TXY motif. Groups A, B and C contained a TEY motif, while group D contained a TDY motif. RNA-seq analysis revealed tissue-specificity of two MsMAPKs and tissue-wide expression of 35 MsMAPKs. Further analysis identified MsMAPK members responsive to drought, salt, and cold stress conditions. Two MsMAPKs (MsMAPK70 and MsMAPK75) responds to salt and cold stresses; two MsMAPKs (MsMAPK60 and MsMAPK73) responds to cold and drought stresses; four MsMAPKs (MsMAPK1, MsMAPK33, MsMAPK64 and MsMAPK71) responds to salt and drought stresses; and two MsMAPKs (MsMAPK5 and MsMAPK7) responded to all three stresses.

CONCLUSION

This study comprehensively identified and analysed the alfalfa MAPK gene family. Candidate genes related to abiotic stresses were screened by analysing the RNA-seq data. The results provide key information for further analysis of alfalfa MAPK gene functions and improvement of stress tolerance.

GhMPK9-GhRAF39_1-GhWRKY40a Regulates the GhERF1b- and GhABF2-Mediated Pathways to Increase Cotton Disease Resistance

Mitogen-activated protein kinase (MAPK) cascade is the center of plant signal transduction system that amplify immune signals into cellular responses by phosphorylating diverse substrates. The MAPK cascade consisting of MAPK kinase kinases (MAPKKKs), MAPK kinases (MAPKKs), and MAPKs is well characterized in plants, in which Raf-like kinases are generally regarded as MAPKKKs. However, it is rarely reported that Raf-like MAPKKKs function as middle regulators to link MAPK and its downstream transcription factors in plant immunity. Verticillium wilt, caused by the soil-borne vascular fungus Verticillium dahliae, is a serious disease in many plants, including cotton. The previous studies showed that GhMPK9 (a MAPK) is involved in the response to Verticillium wilt. Here, the Raf-like kinase GhRAF39_1 is reported as helper regulates the phosphorylation of WRKY transcription factor GhWRKY40a by GhMPK9. The phosphorylated GhWRKY40a can further activate the transcription of GhERF1b to up-regulate defense-related genes while inhibit the transcription of GhABF2 to regulate the stomatal opening, thus improving the resistance to Verticillium wilt in cotton. This study reveals a new signaling module of GhMPK9-GhRAF39_1-GhWRKY40a to regulate GhERF1b- and GhABF2-mediated defense responses, which triggers plant defense against Verticillium wilt.

Identification and expression analysis of MAPK cascade gene family in foxtail millet (Setaria italica)

The mitogen-activated protein kinase (MAPK) cascade pathway is a highly conserved plant cell signaling pathway that plays an important role in plant growth and development and stress response. Currently, MAPK cascade genes have been identified and reported in a variety of plants including Arabidopsis thaliana, Oryza sativa, and Triticum aestivum, but have not been identified in foxtail millet (Setaria italica). In this study, a total of 93 MAPK cascade genes, including 15 SiMAPKs, 10 SiMAPKKs and 68 SiMAPKKKs genes, were identified by genome-wide analysis of foxtail millet, and these genes were distributed on nine chromosomes of foxtail millet. Using phylogenetic analysis, we divided the SiMAPKs and SiMAPKKs into four subgroups, respectively, and the SiMAPKKKs into three subgroups (Raf, ZIK, and MEKK). Whole-genome duplication analysis revealed that there are 14 duplication pairs in the MAPK cascade family in foxtail millet, and they are expanded by segmental replication events. Results from quantitative real-time PCR (qRT-PCR) revealed that the expression levels of most SiMAPKs and SiMAPKKs were changed under both exogenous hormone and abiotic stress treatments, with SiMAPK3 and SiMAPKK4-2 being induced under almost all treatments, while the expression of SiMAPKK5 was repressed. In a nutshell, this study will shed some light on the evolution of MAPK cascade genes and the functional mechanisms underlying MAPK cascade genes in response to hormonal and abiotic stress signaling pathways in foxtail millet (Setaria italica).

Comprehensive genome-wide identification and functional characterization of MAPK cascade gene families in Nelumbo

Mitogen-activated protein kinase (MAPK) cascade signaling pathway plays pivotal roles in various plant biological processes. However, systematic study of MAPK cascade gene families is yet to be conducted in lotus. Herein, 198 putative MAPK genes, including 152 MAP3Ks, 15 MKKs, and 31 MPKs genes were identified in Nelumbo. Segmental duplication was identified as the predominant factor driving MAPK cascade gene family expansion in lotus. MAPK cascade genes in N. nucifera and N. lutea shared high degree of sequence homologies, with 84, 9, and 19 homologous MAP3K, MKK, and MPK gene pairs being detected between the two species, respectively, with most genes predominantly undergoing purifying selection. Gene expression profiling indicated that NnMAPK cascade genes were extensively involved in plant development and submergence stress response. Co-expression analysis revealed potential interaction between transcription factors (TFs) and NnMAPK cascade genes in various biological processes. NnMKK showed predicted interactions with multiple NnMAP3K or NnMPK proteins, which suggested that functional diversity of MAPK cascade genes could be as a result of their complex protein interaction mechanisms. This first systematic analysis of MAPK cascade families in lotus provides deeper insights into their evolutionary dynamics and functional properties, which potentially could be crucial for lotus genetic improvement.

Plant Disease Resistance-Related Signaling Pathways: Recent Progress and Future Prospects

Plant-pathogen interactions induce a signal transmission series that stimulates the plant's host defense system against pathogens and this, in turn, leads to disease resistance responses. Plant innate immunity mainly includes two lines of the defense system, called pathogen-associated molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). There is extensive signal exchange and recognition in the process of triggering the plant immune signaling network. Plant messenger signaling molecules, such as calcium ions, reactive oxygen species, and nitric oxide, and plant hormone signaling molecules, such as salicylic acid, jasmonic acid, and ethylene, play key roles in inducing plant defense responses. In addition, heterotrimeric G proteins, the mitogen-activated protein kinase cascade, and non-coding RNAs (ncRNAs) play important roles in regulating disease resistance and the defense signal transduction network. This paper summarizes the status and progress in plant disease resistance and disease resistance signal transduction pathway research in recent years; discusses the complexities of, and interactions among, defense signal pathways; and forecasts future research prospects to provide new ideas for the prevention and control of plant diseases.

Comparative Transcriptome Analysis of Deep-Rooting and Shallow-Rooting Potato (Solanum tuberosum L.) Genotypes under Drought Stress

The selection and breeding of deep rooting and drought-tolerant varieties has become a promising approach for improving the yield and adaptability of potato (Solanum tuberosum L.) in arid and semiarid areas. Therefore, the discovery of root-development-related genes and drought tolerance signaling pathways in potato is important. In this study, we used deep-rooting (C119) and shallow-rooting (C16) potato genotypes, with different levels of drought tolerance, to achieve this objective. Both genotypes were treated with 150 mM mannitol for 0 h (T0), 2 h (T2), 6 h (T6), 12 h (T12), and 24 h (T24), and their root tissues were subjected to comparative transcriptome analysis. A total of 531, 1571, 1247, and 3540 differentially expressed genes (DEGs) in C16 and 1531, 1108, 674, and 4850 DEGs in C119 were identified in T2 vs. T0, T6 vs. T2, T12 vs. T6, and T24 vs. T12 comparisons, respectively. Gene expression analysis indicated that a delay in the onset of drought-induced transcriptional changes in C16 compared with C119. Functional enrichment analysis revealed genotype-specific biological processes involved in drought stress tolerance. The metabolic pathways of plant hormone transduction and MAPK signaling were heavily involved in the resistance of C16 and C119 to drought, while abscisic acid (ABA), ethylene, and salicylic acid signal transduction pathways likely played more important roles in C119 stress responses. Furthermore, genes involved in root cell elongation and division showed differential expression between the two genotypes under drought stress. Overall, this study provides important information for the marker-assisted selection and breeding of drought-tolerant potato genotypes.

The Protein Phosphatase GhAP2C1 Interacts Together with GhMPK4 to Synergistically Regulate the Immune Response to Fusarium oxysporum in Cotton

The plant mitogen-activated protein kinase (MAPK) cascade plays an important role in mediating responses to biotic and abiotic stresses and is the main pathway through which extracellular stimuli are transduced intracellularly as signals. Our previous research showed that the GhMKK6-GhMPK4 cascade signaling pathway plays an important role in cotton immunity. To further analyze the role and regulatory mechanism of the GhMKK6-GhMPK4 cascade signaling pathway in cotton resistance to Fusarium wilt, we functionally analyzed GhMPK4. Our results show that silencing GhMPK4 reduces cotton tolerance to Fusarium wilt and reduces the expression of several resistance genes. Further experiments revealed that GhMPK4 is similar to GhMKK6, both of whose overexpression cause unfavorable cotton immune response characteristics. By using a yeast two-hybrid screening library and performing a bioinformatics analysis, we screened and identified a negative regulator of the MAPK kinase-protein phosphatase AP2C1. Through the functional analysis of AP2C1, it was found that, after being silenced, GhAP2C1 increased resistance to Fusarium wilt, but GhAP2C1 overexpression caused sensitivity to Fusarium wilt. These findings show that GhAP2C1 interacts together with GhMPK4 to regulate the immune response of cotton to Fusarium oxysporum, which provides important data for functionally analyzing and studying the feedback regulatory mechanism of the MAPK cascade and helps to clarify the regulatory mechanism through which the MAPK cascade acts in response to pathogens.

Boosting of Antioxidants and Alkaloids in Catharanthus roseus Suspension Cultures Using Silver Nanoparticles with Expression of CrMPK3 and STR Genes

Global agricultural systems are under unprecedented pressures due to climate change. Advanced nano-engineering can help increase crop yields while ensuring sustainability. Nanotechnology improves agricultural productivity by boosting input efficiency and reducing waste. Alkaloids as one of the numerous secondary metabolites that serve variety of cellular functions essential for physiological processes. This study tests the competence of silver nanoparticles (AgNPs) in boosting alkaloids accumulation in Catharanthus roseus suspension cultures in relation to the expression of C. roseus Mitogen Activated Protein Kinase 3 (CrMPK3) and Strictosidine Synthase (STR) genes. Five concentrations (5, 10, 15, 20 and 25 mg·L-1) of AgNPs were utilized in addition to deionized water as control. Results reflected binary positive correlations among AgNPs concentration, oxidative stress indicated with increase in hydrogen peroxide and malondialdehyde contents, activities of ascorbate peroxidase and superoxide dismutase, expression of the regulatory gene CrMPK3 and the alkaloid biosynthetic gene STR as well as alkaloids accumulation. These correlations add to the growing evidence that AgNPs can trigger the accumulation of alkaloids in plant cells through a signaling pathway that involves hydrogen peroxide and MAPKs, leading to up-regulation of the biosynthetic genes, including STR gene.

Molecular characterization, expression and interaction of MAPK, MAPKK and MAPKKK genes in upland cotton

Mitogen-activated protein kinase (MAPK) signaling cascades, consisting of three types of sequentially phosphorylated kinases (MAPKKK, MAPKK, and MAPK), play vital roles in various processes including plant development and stress response. In this study, 52 GhMAPKs, 23 GhMAPKKs, and 166 GhMAPKKKs were identified in upland cotton. Chromosomal locations, gene duplication and structure, motifs, cis-regulatory elements, and protein subcellular localization were further analyzed. With the identified MAPK cascade genes in G. arboretum and G. raimondii, a syntenic diagram of three cotton species was constructed. The interactions of seven GhMAPK cascade genes were investigated. Two complete signaling modules were defined: The GhMEKK24/GhMEKK31-GhMAPKK9-GhMAPK10 and GhMEKK3/GhMEKK24/GhMEKK31-GhMAPKK16-GhMAPK10/GhMAPK11 cascades. Moreover, interaction networks and the interaction pairs were combined with their expression patterns and demonstrated that the network mediated by the MAPK signaling cascade participates in abiotic stress signaling. Our research provides a foundation for studying the molecular mechanism of the MAPK signaling pathway under abiotic stress.

Genome-wide identification of MAPK cascade genes reveals the GhMAP3K14-GhMKK11-GhMPK31 pathway is involved in the drought response in cotton

The mitogen-activated protein kinase (MAPK) cascade pathway, which has three components, MAP3Ks, MKKs and MPKs, is involved in diverse biological processes in plants. In the current study, MAPK cascade genes were identified in three cotton species, based on gene homology with Arabidopsis. Selection pressure analysis of MAPK cascade genes revealed that purifying selection occurred among the cotton species. Expression pattern analysis showed that some MAPK cascade genes differentially expressed under abiotic stresses and phytohormones treatments, and especially under drought stress. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) experiments showed extensive interactions between different MAPK cascade proteins. Virus-induced gene silencing (VIGS) assays showed that some MAPK cascade modules play important roles in the drought stress response, and the GhMAP3K14-GhMKK11-GhMPK31 signal pathway was demonstrated to regulate drought stress tolerance in cotton. This study provides new information on the function of MAPK cascade genes in the drought response, and will help direct molecular breeding for improved drought stress tolerance in cotton.

Insights into Drought Stress Signaling in Plants and the Molecular Genetic Basis of Cotton Drought Tolerance

Drought stress restricts plant growth and development by altering metabolic activity and biological functions. However, plants have evolved several cellular and molecular mechanisms to overcome drought stress. Drought tolerance is a multiplex trait involving the activation of signaling mechanisms and differentially expressed molecular responses. Broadly, drought tolerance comprises two steps: stress sensing/signaling and activation of various parallel stress responses (including physiological, molecular, and biochemical mechanisms) in plants. At the cellular level, drought induces oxidative stress by overproduction of reactive oxygen species (ROS), ultimately causing the cell membrane to rupture and stimulating various stress signaling pathways (ROS, mitogen-activated-protein-kinase, Ca2+, and hormone-mediated signaling). Drought-induced transcription factors activation and abscisic acid concentration co-ordinate the stress signaling and responses in cotton. The key responses against drought stress, are root development, stomatal closure, photosynthesis, hormone production, and ROS scavenging. The genetic basis, quantitative trait loci and genes of cotton drought tolerance are presented as examples of genetic resources in plants. Sustainable genetic improvements could be achieved through functional genomic approaches and genome modification techniques such as the CRISPR/Cas9 system aid the characterization of genes, sorted out from stress-related candidate single nucleotide polymorphisms, quantitative trait loci, and genes. Exploration of the genetic basis for superior candidate genes linked to stress physiology can be facilitated by integrated functional genomic approaches. We propose a third-generation sequencing approach coupled with genome-wide studies and functional genomic tools, including a comparative sequenced data (transcriptomics, proteomics, and epigenomic) analysis, which offer a platform to identify and characterize novel genes. This will provide information for better understanding the complex stress cellular biology of plants.

Identification and Characterization of Mitogen-Activated Protein Kinase (MAPK) Genes in Sunflower (Helianthus annuus L.)

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that regulate biotic and abiotic stresses in plants through signaling cascades comprised of three major subfamilies: MAP Kinase (MPK), MAPK Kinase (MKK), and MAPKK Kinase (MKKK). The main objectives of this research were to conduct genome-wide identification of MAPK genes in Helianthus annuus and examine functional divergence of these genes in relation to those in nine other plant species (Amborella trichopoda, Aquilegia coerulea, Arabidopsis thaliana, Daucus carota, Glycine max, Oryza sativa, Solanum lycopersicum, Sphagnum fallax, and Vitis vinifera), representing diverse taxonomic groups of the Plant Kingdom. A Hidden Markov Model (HMM) profile of the MAPK genes utilized reference sequences from A. thaliana and G. max, yielding a total of 96 MPKs and 37 MKKs in the genomes of A. trichopoda, A. coerulea, C. reinhardtii, D. carota, H. annuus, S. lycopersicum, and S. fallax. Among them, 28 MPKs and eight MKKs were confirmed in H. annuus. Phylogenetic analyses revealed four clades within each subfamily. Transcriptomic analyses showed that at least 19 HaMPK and seven HaMKK genes were induced in response to salicylic acid (SA), sodium chloride (NaCl), and polyethylene glycol (Peg) in leaves and roots. Of the seven published sunflower microRNAs, five microRNA families are involved in targeting eight MPKs. Additionally, we discussed the need for using MAP Kinase nomenclature guidelines across plant species. Our identification and characterization of MAP Kinase genes would have implications in sunflower crop improvement, and in advancing our knowledge of the diversity and evolution of MAPK genes in the Plant Kingdom.

Differential Responses of OsMPKs in IR56 Rice to Two BPH Populations of Different Virulence Levels

The conserved mitogen-activated protein kinase (MAPK) cascades play vital roles in plant defense responses against pathogens and insects. In the current study, the expression profiles of 17 OsMPKs were determined in the TN1 and IR56 rice varieties under the infestation of brown planthopper (BPH), one of the most destructive hemimetabolous rice pests. The virulent IR56 BPH population (IR56-BPH) and the avirulent TN1 BPH population (TN-BPH) were used to reveal the roles of OsMPKs in the compatible (IR56-BPH infested on the TN1 and IR56 rice varieties, and TN1-BPH infested on the TN1 rice variety) and the incompatible (TN1-BPH infested on the IR56 rice variety) interaction. The statistical analysis revealed that rice variety, BPH population type, and infestation period have significant effects on the transcription of OsMPKs. Out of these genes, five OsMPKs (OsMPK1, OsMPK3, OsMPK7, OsMPK14, and OsMPK16) were found to exhibit upregulated expression only during incompatible interaction. Six OsMPKs (OsMPK4, OsMPK5, OsMPK8, OsMPK9, OsMPK12, and OsMPK13) were associated with both incompatible and compatible interactions. The transcription analysis of salicylic acid, jasmonic acid, and ethylene phytohormone signaling genes revealed their roles during the rice⁻BPH interactions. The upregulated expression of OsC4H, OsCHS, and OsCHI in the incompatible interaction implied the potential defense regulatory roles of phenylpropanoids. In both varieties, the elevated transcript accumulations of OsGST and OsSOD, and the increased enzyme activities of POD, SOD, and GST at 1 day post-infestation (dpi), but not at 3 dpi, indicated that reactive oxygen species (ROS) signaling might be an early event in rice⁻BPH interactions. Furthermore, upregulated transcription of OsLecRK3 and OsLecRK4 was found only during an incompatible interaction, suggesting their involvement in the BPH resistance response in the IR56 rice variety. Lastly, based on the findings of this study, we have proposed a model of interactions of IR56 rice with TN1-BPH and IR56-BPH that depicts the resistance and susceptibility reactions, respectively.