Genome-wide identification and analysis of MAPK and MAPKK gene families in Brachypodium distachyon

A close-up of regulatory networks and signaling pathways of MKK5 in biotic and abiotic stresses

Mitogen-activated protein Kinase Kinase 5 (MKK5) is a central hub in the complex phosphorylation chain reaction of the Mitogen-activated protein kinases (MAPK) cascade, regulating plant responses to biotic and abiotic stresses. This review manuscript aims to provide a comprehensive analysis of the regulatory mechanism of the MKK5 involved in stress adaptation. This review will delve into the intricate post-transcriptional and post-translational modifications of the MKK5, discussing how they affect its expression, activity, and subcellular localization in response to stress signals. We also discuss the integration of the MKK5 into complex signaling pathways, orchestrating plant immunity against pathogens and its modulating role in regulating abiotic stresses, such as: drought, cold, heat, and salinity, through the phytohormonal signaling pathways. Furthermore, we highlight potential applications of the MKK5 for engineering stress-resilient crops and provide future perspectives that may pave the way for future studies. This review manuscript aims to provide valuable insights into the mechanisms underlying MKK5 regulation, bridge the gap from numerous previous findings, and offer a firm base in the knowledge of MKK5, its regulating roles, and its involvement in environmental stress regulation.

Genome-wide identification, classification and expression analysis of the heat shock transcription factor family in Garlic (Allium sativum L.)

BACKGROUND

The heat shock transcription factor (HSF) plays a crucial role in the regulatory network by coordinating responses to heat stress as well as other stress signaling pathways. Despite extensive studies on HSF functions in various plant species, our understanding of this gene family in garlic, an important crop with nutritional and medicinal value, remains limited. In this study, we conducted a comprehensive investigation of the entire garlic genome to elucidate the characteristics of the AsHSF gene family.

RESULTS

In this study, we identified a total of 17 AsHSF transcription factors. Phylogenetic analysis classified these transcription factors into three subfamilies: Class A (9 members), Class B (6 members), and Class C (2 members). Each subfamily was characterized by shared gene structures and conserved motifs. The evolutionary features of the AsHSF genes were investigated through a comprehensive analysis of chromosome location, conserved protein motifs, and gene duplication events. These findings suggested that the evolution of AsHSF genes is likely driven by both tandem and segmental duplication events. Moreover, the nucleotide diversity of the AsHSF genes decreased by only 0.0002% from wild garlic to local garlic, indicating a slight genetic bottleneck experienced by this gene family during domestication. Furthermore, the analysis of cis-acting elements in the promoters of AsHSF genes indicated their crucial roles in plant growth, development, and stress responses. qRT-PCR analysis, co-expression analysis, and protein interaction prediction collectively highlighted the significance of Asa6G04911. Subsequent experimental investigations using yeast two-hybridization and yeast induction experiments confirmed its interaction with HSP70/90, reinforcing its significance in heat stress.

CONCLUSIONS

This study is the first to unravel and analyze the AsHSF genes in garlic, thereby opening up new avenues for understanding their functions. The insights gained from this research provide a valuable resource for future investigations, particularly in the functional analysis of AsHSF genes.

Identification and analysis of MAPK cascade gene families of Camellia oleifera and their roles in response to cold stress

BACKGROUND

Low-temperature severely limits the growth and development of Camellia oleifera (C. oleifera). The mitogen-activated protein kinase (MAPK) cascade plays a key role in the response to cold stress.

METHODS AND RESULTS

Our study aims to identify MAPK cascade genes in C. oleifera and reveal their roles in response to cold stress. In our study, we systematically identified and analyzed the MAPK cascade gene families of C. oleifera, including their physical and chemical properties, conserved motifs, and multiple sequence alignments. In addition, we characterized the interacting networks of MAPKK kinase (MAPKKK)-MAPK kinase (MAPKK)-MAPK in C. oleifera. The molecular mechanism of cold stress resistance of MAPK cascade genes in wild C. oleifera was analyzed by differential gene expression and real-time quantitative reverse transcription-PCR (qRT-PCR).

CONCLUSION

In this study, 21 MAPKs, 4 MAPKKs and 55 MAPKKKs genes were identified in the leaf transcriptome of C. oleifera. According to the phylogenetic results, MAPKs were divided into 4 groups (A, B, C and D), MAPKKs were divided into 3 groups (A, B and D), and MAPKKKs were divided into 2 groups (MEKK and Raf). Motif analysis showed that the motifs in each subfamily were conserved, and most of the motifs in the same subfamily were basically the same. The protein interaction network based on Arabidopsis thaliana (A. thaliana) homologs revealed that MAPK, MAPKK, and MAPKKK genes were widely involved in C. oleifera growth and development and in responses to biotic and abiotic stresses. Gene expression analysis revealed that the CoMAPKKK5/CoMAPKKK43/CoMAPKKK49-CoMAPKK4-CoMAPK8 module may play a key role in the cold stress resistance of wild C. oleifera at a high-elevation site in Lu Mountain (LSG). This study can facilitate the mining and utilization of genetic resources of C. oleifera with low-temperature tolerance.

Genome-wide identification of MAPK family in papaya (Carica papaya) and their involvement in fruit postharvest ripening

BACKGROUND

Papaya (Carica papaya) is an economically important fruit cultivated in the tropical and subtropical regions of China. However, the rapid softening rate after postharvest leads to a short shelf-life and considerable economic losses. Accordingly, understanding the mechanisms underlying fruit postharvest softening will be a reasonable way to maintain fruit quality and extend its shelf-life.

RESULTS

Mitogen-activated protein kinases (MAPKs) are conserved and play essential roles in response to biotic and abiotic stresses. However, the MAPK family remain poorly studied in papaya. Here, a total of nine putative CpMAPK members were identified within papaya genome, and a comprehensive genome-wide characterization of the CpMAPKs was performed, including evolutionary relationships, conserved domains, gene structures, chromosomal locations, cis-regulatory elements and expression profiles in response to phytohormone and antioxidant organic compound treatments during fruit postharvest ripening. Our findings showed that nearly all CpMAPKs harbored the conserved P-loop, C-loop and activation loop domains. Phylogenetic analysis showed that CpMAPK members could be categorized into four groups (A-D), with the members within the same groups displaying high similarity in protein domains and intron-exon organizations. Moreover, a number of cis-acting elements related to hormone signaling, circadian rhythm, or low-temperature stresses were identified in the promoters of CpMAPKs. Notably, gene expression profiles demonstrated that CpMAPKs exhibited various responses to 2-chloroethylphosphonic acid (ethephon), 1-methylcyclopropene (1-MCP) and the combined ascorbic acid (AsA) and chitosan (CTS) treatments during papaya postharvest ripening. Among them, both CpMAPK9 and CpMAPK20 displayed significant induction in papaya flesh by ethephon treatment, and were pronounced inhibition after AsA and CTS treatments at 16 d compared to those of natural ripening control, suggesting that they potentially involve in fruit postharvest ripening through ethylene signaling pathway or modulating cell wall metabolism.

CONCLUSION

This study will provide some valuable insights into future functional characterization of CpMAPKs, and hold great potential for further understanding the molecular mechanisms underlying papaya fruit postharvest ripening.

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).

Impact of high atmospheric carbon dioxide on the biotic stress response of the model cereal species Brachypodium distachyon

Losses due to disease and climate change are among the most important issues currently facing crop production. It is therefore important to establish the impact of climate change, and particularly of high carbon dioxide (hCO2), on plant immunity in cereals, which provide 60% of human calories. The aim of this study was to determine if hCO2 impacts Brachypodium distachyon immunity, a model plant for temperate cereals. Plants were grown in air (430 ppm CO2) and at two high CO2 conditions, one that is relevant to projections within the coming century (1000 ppm) and a concentration sufficient to saturate photosynthesis (3000 ppm). The following measurements were performed: phenotyping and growth, salicylic acid contents, pathogen resistance tests, and RNAseq analysis of the transcriptome. Improved shoot development was observed at both 1000 and 3000 ppm. A transcriptomic analysis pointed to an increase in primary metabolism capacity under hCO2. Alongside this effect, up-regulation of genes associated with secondary metabolism was also observed. This effect was especially evident for the terpenoid and phenylpropanoid pathways, and was accompanied by enhanced expression of immunity-related genes and accumulation of salicylic acid. Pathogen tests using the fungus Magnaporthe oryzae revealed that hCO2 had a complex effect, with enhanced susceptibility to infection but no increase in fungal development. The study reveals that immunity in B. distachyon is modulated by growth at hCO2 and allows identification of pathways that might play a role in this effect.

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.

The effect of a-Lipoic acid (ALA) on oxidative stress, inflammation, and apoptosis in high glucose-induced human corneal epithelial cells

PURPOSE

Oxidative stress and inflammation had been proved to play important role in the progression of diabetic keratopathy (DK). The excessive accumulation of AGEs and their bond to AGE receptor (RAGE) in corneas that cause the formation of oxygen radicals and the release of inflammatory cytokines, induce cell apoptosis. Our current study was aimed to evaluate the effect of ALA on AGEs accumulation as well as to study the molecular mechanism of ALA against AGE-RAGE axis mediated oxidative stress, apoptosis, and inflammation in HG-induced HCECs, so as to provide cytological basis for the treatment of DK.

METHODS

HCECs were cultured in a variety concentration of glucose medium (5.5, 10, 25, 30, 40, and 50 mM) for 48 h. The cell proliferation was evaluated by CCK-8 assay. Apoptosis was investigated with the Annexin V- fluorescein isothiocyanate (V-FITC)/PI kit, while, the apoptotic cells were determined by flow cytometer and TUNEL cells apoptosis Kit. According to the results of cell proliferation and cell apoptosis, 25 mM glucose medium was used in the following HG experiment. The effect of ALA on HG-induced HCECs was evaluated. The HCECs were treated with 5.5 mM glucose (normal glucose group, NG group), 5.5 mM glucose + 22.5 mM mannitol (osmotic pressure control group, OP group), 25 mM glucose (high glucose group, HG group) and 25 mM glucose + ALA (HG + ALA group) for 24 and 48 h. The accumulation of intracellular AGEs was detected by ELISA kit. The RAGE, catalase (CAT), superoxide dismutase 2 (SOD2), cleaved cysteine-aspartic acid protease-3 (Cleaved caspase-3), Toll-like receptors 4 (TLR4), Nod-like receptor protein 3 (NLRP3) inflammasome, interleukin 1 beta (IL-1 ß), and interleukin 18 (IL-18) were quantified by RT-PCR, Western blotting, and Immunofluorescence, respectively. Reactive oxygen species (ROS) production was evaluated by fluorescence microscope and fluorescence microplate reader.

RESULTS

When the glucose medium was higher than 25 mM, cell proliferation was significantly inhibited and apoptosis ratio was increased (P < 0.001). In HG environment, ALA treatment alleviated the inhibition of HCECs in a dose-dependent manner, 25 μM ALA was the minimum effective dose. ALA could significantly reduce the intracellular accumulation of AGEs (P < 0.001), activate protein and genes expression of CAT and SOD2 (P < 0.001), and therefore inhibited ROS-induced oxidative stress and cells apoptosis. Besides, ALA could effectively down-regulate the protein and gene level of RAGE, TLR4, NLRP3, IL-1B, IL-18 (P < 0.05), and therefore alleviated AGEs-RAGE-TLR4-NLRP3 pathway-induced inflammation in HG-induced HCECs.

CONCLUSION

Our study indicated that ALA could be a desired treatment for DK due to its potential capacity of reducing accumulation of advanced glycation end products (AGEs) and down-regulating AGE-RAGE axis-mediated oxidative stress, cell apoptosis, and inflammation in high glucose (HG)-induced human corneal epithelial cells (HCECs), which may provide cytological basis for therapeutic targets that are ultimately of clinical benefit.

Genome-wide identification evolution and expression of vestigial-like gene family in chicken

Vestigial-like (Vgll) genes are widespread in vertebrates and play an important role in muscle development. In this study, we used bioinformatics methods to systematically identify the chicken VGLL family in the whole genome and investigated its evolutionary history and gene structure features. Tissue expression spectra combined with real-time PCR data were used to analyze the organizational expression pattern of the genes. Based on the maximum likelihood method, a phylogenetic tree of the VGLL family was constructed, and 94 VGLL genes were identified in 24 breeds, among which four VGLL family genes were identified in the chicken genome. Ten motifs were detected in the VGLL genes, and the analysis of introns combined with gene structure revealed that the family was conserved during evolution. Tissue expression analysis suggested that the expression profiles of the VGLL family genes in 16 tissues differed between LU Shi and AA broilers. In addition, a single gene (VGLL2) showed increased expression in chickens at embryonic days 10-16 and was involved in the growth and development of skeletal muscle in chickens in the embryonic stage. In summary, VGLL genes are involved in chicken muscle growth and development, which provides useful information for subsequent functional studies of VGLL genes.

Genome-wide identification of the mitogen-activated protein kinase kinase kinase (MAPKKK) in pear (Pyrus bretschneideri) and their functional analysis in response to black spot

Identification of MAPKKK genes in pear and functional characterization of PbrMAPKKK82 in response to pear black spot. Mitogen-activated protein kinase kinase kinase (MAPKKK) is located upstream of the MAPK cascade pathway. This region senses extracellular stimuli via the signaling molecule or by themselves and is activated by phosphorylation. In this study, we identified 108 PbrMAPKKK genes from the pear genome. The genes were divided into three subfamilies and contained the conserved domain. Except for chromosome 7, there were 93 PbrMAPKKK genes randomly distributed on 16 out of the 17 chromosomes, while 15 PbrMAPKKK genes were detected on unknown chromosomes. They largely originated from whole-genome duplication (WGD) and dispersed events. In the expression analysis of PbrMAPKKK genes in seven pear tissue types by using a database, 20 PbrMAPKKK genes were selected to verify the expression associated with different resistance in two varieties by quantitative real-time PCR (qRT-PCR). The results showed that PbrMAPKKK12, PbrMAPKKK13, PbrMAPKKK53, PbrMAPKKK60, PbrMAPKKK65, PbrMAPKKK82, PbrMAPKKK83, and PbrMAPKKK96 were correlated with black spot resistance. PbrMAPKKK3, PbrMAPKKK9, PbrMAPKKK11, PbrMAPKKK34, PbrMAPKKK80, PbrMAPKKK81, PbrMAPKKK99, and PbrMAPKKK100 were correlated with black spot susceptibility, while the PbrMAPKKK gene positively responded to the life process of pear resistance to black spot. Furthermore, virus-induced gene silencing (VIGS) indicated that the PbrMAPKKK82 gene enhanced resistance to pear black spot disease.

Arbuscular mycorrhizal fungi enhance photosynthesis and drought tolerance by regulating MAPK genes expressions of Populus simonii × P. nigra

Arbuscular mycorrhizal fungi (AMF) promote plants to absorb more water and nutrients and improve their stress resistance. As the main signal transducer, the mitogen-activated protein kinase (MAPK) cascade plays a vital role in drought stress. However, how the MAPK family genes of mycorrhizal plants respond to stress is still not clear. Our study analyzed physiological indexes and expression profiles of MAPK family genes of Populus simonii × P. nigra under two inoculation treatments (inoculated with or without Rhizophagus irregularis) and two water conditions (well-watered or drought stress). The results showed that the stronger photosynthesis of mycorrhizal plants may be mediated by MAPK genes induced by AMF. Mycorrhizal plants showed lower oxidative damage and drought sensitivity. Mycorrhiza downregulated the expression of PsnMAPK7-2, PsnMAPK16-1, PsnMAPK19-2, and PsnMAPK20-2 which negatively regulate drought tolerance and induced specific PsnMAPKs in roots which activate transcription factors to regulate downstream gene expressions, enhancing drought tolerance. This is the first time to identify part of the MAPK gene family of P. simonii × P. nigra at the genome level and study MAPK genes in mycorrhizal forest trees. This is helpful to understand the function of the MAPK gene family in response to drought of mycorrhizal plants and lays a foundation for afforestation by using mycorrhizal saplings.

Genome-wide identification and characterization of the MAPKKK, MKK, and MPK families in Chinese elite maize inbred line Huangzaosi

Mitogen-activated protein kinase (MAPK or MPK) cascades consist of three protein kinase components, MAPK kinase kinases (MAPKKKs), MAPK kinases (MKKs and MPKs), which are indispensable for various plant physiological processes. The functions of MAPK families have been extensively studied in maize (Zea mays L.) and other plant species, but little is known about MAPK families in the elite Chinese maize line Huangzaosi (hzs). In this study, we observed that overall performance of Huangzaosi was substantially better than that of B73 under drought conditions at the seedling and V16 stages with a favorable root/canopy ratio. In silico analyses identified 72, 10, and 24 MAPKKKs, MKKs, and MPKs, respectively, in Huangzaosi. Examinations of phylogenetic relationships among Arabidopsis thaliana (L.) Heynh., rice (Oryza sativa L.), and maize (lines B73 and hzs), gene structures, conserved protein motifs, and chromosomal locations revealed their evolutionary relationships. The basal gene expression levels and tissue specificities of all three MAPK families in hzs reflected the diversity in the MAPK functions related to growth and development. The quantitative real-time polymerase chain reaction (qPCR) assay indicated that certain MAPK genes with high basal expression levels in the primary and crown roots responded differentially to drought between B73 and hzs, suggesting that these genes may contribute to their distinct drought tolerance at different developmental stages. The important information regarding the evolution and expression of hzs MAPK family members generated in this study provides a new avenue for the better understanding on the regulatory mechanism of MAPK cascade in the core inbred line hzs, which may be useful to guide the development of new maize cultivars with desirable traits (e.g., drought resistance).

Genome-Wide Identification of MAPK, MAPKK, and MAPKKK Gene Families in Fagopyrum tataricum and Analysis of Their Expression Patterns Under Abiotic Stress

The mitogen-activated protein kinase (MAPK) cascade is a highly conserved signal transduction pathway, ubiquitous in eukaryotes, such as animals and plants. The MAPK cascade has a dominant role in regulating plant adaptation to the environment, such as through stress responses, osmotic adjustment, and processes that modulate pathogenicity. In the present study, the MAPK cascade gene family was identified in Fagopyrum tataricum (Tartary buckwheat), based on complete genome sequence data. Using phylogenetic tree, conservative motif, and chromosome location analyses, a total of 65 FtMAPK cascade genes, distributed on five chromosomes, were classified into three families: MAPK (n = 8), MAPKK (n = 1), and MAPKKK (n = 56). Transcriptome data from Tartary buckwheat seedlings grown under different light conditions demonstrated that, under blue and red light, the expression levels of 18 and 36 FtMAPK cascade genes were up-regulated and down-regulated, respectively. Through qRT-PCR experiments, it was observed that FtMAPK5, FtMAPKK1, FtMAPKKK8, FtMAPKKK10, and FtMAPKKK24 gene expression levels in the Tartary buckwheat seedlings increased under three types of abiotic stress: drought, salt, and high temperature. A co-expression network of FtMAPK cascade genes was constructed, based on gene expression levels under different light conditions, and co-expressed genes annotated by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses, which identified numerous transcription factors related to plant abiotic stress. The authors conclude that FtMAPK cascade genes have important roles in the growth and development of Tartary buckwheat, as well as its responses to abiotic stress.

Ectopic Expression of OsJAZs Alters Plant Defense and Development

A key step in jasmonic acid (JA) signaling is the ligand-dependent assembly of a coreceptor complex comprising the F-box protein COI1 and JAZ transcriptional repressors. The assembly of this receptor complex results in proteasome-mediated degradation of JAZ repressors, which in turn bind and repress MYC transcription factors. Many studies on JAZs have been performed in Arabidopsis thaliana, but the function of JAZs in rice is largely unknown. To systematically reveal the function of OsJAZs, in this study, we compared the various phenotypes resulting from 13 OsJAZs via ectopic expression in Arabidopsis thaliana and the phenotypes of 12 AtJAZs overexpression (OE) lines. Phylogenetic analysis showed that the 25 proteins could be divided into three major groups. Yeast two-hybrid (Y2H) assays revealed that most OsJAZ proteins could form homodimers or heterodimers. The statistical results showed that the phenotypes of the OsJAZ OE plants were quite different from those of AtJAZ OE plants in terms of plant growth, development, and immunity. As an example, compared with other JAZ OE plants, OsJAZ11 OE plants exhibited a JA-insensitive phenotype and enhanced resistance to Pst DC3000. The protein stability after JA treatment of OsJAZ11 emphasized the specific function of the protein. This study aimed to explore the commonalities and characteristics of different JAZ proteins functions from a genetic perspective, and to screen genes with disease resistance value. Overall, the results of this study provide insights for further functional analysis of rice JAZ family proteins.

Genome-wide identification and expression profiling of thes MAPK, MAPKK, and MAPKKK gene families in Ophiocordyceps sinensis

Mitogen-activated protein kinase (MAPK) cascades have a universal cell signaling mechanism in eukaryotes. A typical MAPK signal transduction module comprises three kinds of sequentially phosphorylated protein kinases: MAPK, Mitogen-activated protein kinase kinase (MAPKK), and Mitogen-activated protein kinase kinase kinase (MAPKKK). However, little is known regarding the genes involved in MAPK cascades in Ophiocordyceps sinensis. Nine genes (three MAPK, three MAPKK, and three MAPKKK) were identified in this study. The MAPK, MAPKK, and MAPKKK genes were divided into three subfamilies, according to the phylogenetic analysis. TEY and TGY represented the activation domains of the MAPKs; the corresponding domains in MAPKKs were SDIWS and SDVWS, and those in the MAPKKs were GSVFYWMAPEV and GTPMYMSPEV. Transcription data analysis and quantitative real-time polymerase chain reaction showed that the MAPK cascade was related to the growth of the fruiting body. This is the first study to report a genome-wide identification of the MAPK, MAPKK, and MAPKK gene families in O. sinensis.

Mitogen-Activated Protein Kinase Expression Profiling Revealed Its Role in Regulating Stress Responses in Potato (Solanum tuberosum)

Mitogen-activated protein kinase (MAPK) cascades are the universal signal transduction networks that regulate cell growth and development, hormone signaling, and other environmental stresses. However, their essential contribution to plant tolerance is very little known in the potato (Solanum tuberosum) plant. The current study carried out a genome-wide study of StMAPK and provided a deep insight using bioinformatics tools. In addition, the relative expression of StMAPKs was also assessed in different plant tissues. The similarity search results identified a total of 22 StMAPK genes in the potato genome. The sequence alignment also showed conserved motif TEY/TDY in most StMAPKs with conserved docking LHDXXEP sites. The phylogenetic analysis divided all 22 StMAPK genes into five groups, i.e., A, B, C, D, and E, showing some common structural motifs. In addition, most of the StMAPKs were found in a cluster form at the terminal of chromosomes. The promoter analysis predicted several stress-responsive Cis-acting regulatory elements in StMAPK genes. Gene duplication under selection pressure also indicated several purifying and positive selections in StMAPK genes. In potato, StMAPK2, StMAPK6, and StMAPK19 showed a high expression in response to heat stress. Under ABA and IAA treatment, the expression of the total 20 StMAPK genes revealed that ABA and IAA played an essential role in this defense process. The expression profiling and real-time qPCR (RT-qPCR) exhibited their high expression in roots and stems compared to leaves. These results deliver primary data for functional analysis and provide reference data for other important crops.

A phylogenetic study of the members of the MAPK and MEK families across Viridiplantae

Protein phosphorylation is regulated by the activity of enzymes generically known as kinases. One of those kinases is Mitogen-Activated Protein Kinases (MAPK), which operate through a phosphorylation cascade conformed by members from three related protein kinase families namely MAPK kinase kinase (MEKK), MAPK kinase (MEK), and MAPK; these three acts hierarchically. Establishing the evolution of these proteins in the plant kingdom is an interesting but complicated task because the current MAPK, MAPKK, and MAPKKK subfamilies arose from duplications and subsequent sub-functionalization during the early stage of the emergence of Viridiplantae. Here, an in silico genomic analysis was performed on 18 different plant species, which resulted in the identification of 96 genes not previously annotated as components of the MAPK (70) and MEK (26) families. Interestingly, a deeper analysis of the sequences encoded by such genes revealed the existence of putative domains not previously described as signatures of MAPK and MEK kinases. Additionally, our analysis also suggests the presence of conserved activation motifs besides the canonical TEY and TDY domains, which characterize the MAPK family.

Comparative analysis of MAPK and MKK gene families reveals differential evolutionary patterns in Brachypodium distachyon inbred lines

BACKGROUND

Mitogen-activated protein kinase (MAPK) cascades are involved with signal transduction in almost every aspect of plant growth and development, as well as biotic and abiotic stress responses. The evolutionary analysis of MAPKs and MKKs in individual or entire plant species has been reported, but the evolutionary patterns in the diverse inbred lines of Brachypodium distachyon are still unclear.

RESULTS

We conducted the systematical molecular evolutionary analysis of B. distachyon. A total of 799 MAPKs and 618 MKKs were identified from 53 B. distachyon inbred lines. Remarkably, only three inbred lines had 16 MPKs and most of those inbred lines lacked MPK7-2 members, whereas 12 MKKs existed in almost all B. distachyon inbred lines. Phylogenetic analysis indicated that MAPKs and MKKs were divided into four groups as previously reported, grouping them in the same branch as corresponding members. MPK21-2 was the exception and fell into two groups, which may be due to their exon-intron patterns, especially the untranslated regions (UTRs). We also found that differential evolution patterns of MKK10 paralogues from ancient tandem duplicates may have undergone functional divergence. Expression analyses suggested that MAPKs and MKKs likely played different roles in different genetic contexts within various tissues and with abiotic stresses.

CONCLUSION

Our study revealed that UTRs affected the structure and evolution of MPK21-2 genes and the differential evolution of MKK10 paralogues with ancient tandem duplication might have functional divergences. Our findings provide new insights into the functional evolution of genes in closely inbred lines.

Genome-wide Identification, Evolution and Expression Analysis of Basic Helix-loop-helix (bHLH) Gene Family in Barley (Hordeum vulgare L.)

Background

The basic helix-loop-helix (bHLH) transcription factor is one of the most important gene families in plants, playing a key role in diverse metabolic, physiological, and developmental processes. Although it has been well characterized in many plants, the significance of the bHLH family in barley is not well understood at present.

Methods

Through a genome-wide search against the updated barley reference genome, the genomic organization, evolution and expression of the bHLH family in barley were systematically analyzed.

Results

We identified 141 bHLHs in the barley genome (HvbHLHs) and further classified them into 24 subfamilies based on phylogenetic analysis. It was found that HvbHLHs in the same subfamily shared a similar conserved motif composition and exon-intron structures. Chromosome distribution and gene duplication analysis revealed that segmental duplication mainly contributed to the expansion of HvbHLHs and the duplicated genes were subjected to strong purifying selection. Furthermore, expression analysis revealed that HvbHLHs were widely expressed in different tissues and also involved in response to diverse abiotic stresses. The co-expression network was further analyzed to underpin the regulatory function of HvbHLHs. Finally, 25 genes were selected for qRT-PCR validation, the expression profiles of HvbHLHs showed diverse patterns, demonstrating their potential roles in relation to stress tolerance regulation.

Conclusion

This study reported the genome organization, evolutionary characteristics and expression profile of the bHLH family in barley, which not only provide the targets for further functional analysis, but also facilitate better understanding of the regulatory network bHLH genes involved in stress tolerance in barley.

Genome-Wide Identification and Analysis of MKK and MAPK Gene Families in Brassica Species and Response to Stress in Brassica napus

Mitogen-activated protein kinase (MAPK) cascades are common and conserved signal transduction pathways and play important roles in various biotic and abiotic stress responses and growth and developmental processes in plants. With the advancement of sequencing technology, more systematic genetic information is being explored. The work presented here focuses on two protein families in Brassica species: MAPK kinases (MKKs) and their phosphorylation substrates MAPKs. Forty-seven MKKs and ninety-two MAPKs were identified and extensively analyzed from two tetraploid (B. juncea and B. napus) and three diploid (B. nigra, B. oleracea, and B. rapa) Brassica species. Phylogenetic relationships clearly distinguished both MKK and MAPK families into four groups, labeled A–D, which were also supported by gene structure and conserved protein motif analysis. Furthermore, their spatial and temporal expression patterns and response to stresses (cold, drought, heat, and shading) were analyzed, indicating that BnaMKK and BnaMAPK transcript levels were generally modulated by growth, development, and stress signals. In addition, several protein interaction pairs between BnaMKKs and C group BnaMAPKs were detected by yeast two-hybrid assays, in which BnaMKK3 and BnaMKK9 showed strong interactions with BnaMAPK1/2/7, suggesting that interaction between BnaMKKs and C group BnaMAPKs play key roles in the crosstalk between growth and development processes and abiotic stresses. Taken together, our data provide a deeper foundation for the evolutionary and functional characterization of MKK and MAPK gene families in Brassica species, paving the way for unraveling the biological roles of these important signaling molecules in plants.

Identification, Characterization, and Expression Profile Analysis of the mTERF Gene Family and Its Role in the Response to Abiotic Stress in Barley (Hordeum vulgare L.)

Plant mitochondrial transcription termination factor (mTERF) family regulates organellar gene expression (OGE) and is functionally characterized in diverse species. However, limited data are available about its functions in the agriculturally important cereal barley (Hordeum vulgare L.). In this study, we identified 60 mTERFs in the barley genome (HvmTERFs) through a comprehensive search against the most updated barley reference genome, Morex V2. Then, phylogenetic analysis categorized these genes into nine subfamilies, with approximately half of the HvmTERFs belonging to subfamily IX. Members within the same subfamily generally possessed conserved motif composition and exon-intron structure. Both segmental and tandem duplication contributed to the expansion of HvmTERFs, and the duplicated gene pairs were subjected to strong purifying selection. Expression analysis suggested that many HvmTERFs may play important roles in barley development (e.g., seedlings, leaves, and developing inflorescences) and abiotic stresses (e.g., cold, salt, and metal ion), and HvmTERF21 and HvmTERF23 were significant induced by various abiotic stresses and/or phytohormone treatment. Finally, the nucleotide diversity was decreased by only 4.5% for HvmTERFs during the process of barley domestication. Collectively, this is the first report to characterize HvmTERFs, which will not only provide important insights into further evolutionary studies but also contribute to a better understanding of the potential functions of HvmTERFs and ultimately will be useful in future gene functional studies.

MAPK cascade gene family in Camellia sinensis: In-silico identification, expression profiles and regulatory network analysis

BACKGROUND

Mitogen Activated Protein Kinase (MAPK) cascade is a fundamental pathway in organisms for signal transduction. Though it is well characterized in various plants, there is no systematic study of this cascade in tea.

RESULT

In this study, 5 genes of Mitogen Activated Protein Kinase Kinase (MKK) and 16 genes of Mitogen Activated Protein Kinase (MPK) in Camellia sinensis were found through a genome-wide search taking Arabidopsis thaliana as the reference genome. Also, phylogenetic relationships along with structural analysis which includes gene structure, location as well as protein conserved motifs and domains, were systematically examined and further, predictions were validated by the results. The plant species taken for comparative study clearly displayed segmental duplication, which was a significant candidate for MAPK cascade expansion. Also, functional interaction was carried out in C. sinensis based on the orthologous genes in Arabidopsis. The expression profiles linked to various stress treatments revealed wide involvement of MAPK and MAPKK genes from Tea in response to various abiotic factors. In addition, the expression of these genes was analysed in various tissues.

CONCLUSION

This study provides the targets for further comprehensive identification, functional study, and also contributed for a better understanding of the MAPK cascade regulatory network in C. sinensis.

Cloning, molecular and functional characterization by overexpression in Arabidopsis of MAPKK genes from grapevine (Vitis vinifera)

BACKGROUND

The mitogen-activated protein kinases (MAPKs), as a part of the MAPKKK-MAPKK-MAPK cascade, play crucial roles in plant development as an intracellular signal transduction pathway to respond various environmental signals. However, few MAPKK have been functionally characterized in grapevine.

RESULTS

In the study, five MAPKK (MKK) members were identified in grapevine (cultivar 'Pinot Noir'), cloned and designated as VvMKK1-VvMKK5. A phylogenetic analysis grouped them into four sub-families based on the similarity of their conserved motifs and gene structure to Arabidopsis MAPKK members. qRT-PCR results indicated that the expression of VvMKK1, VvMKK2, VvMKK4, and VvMKK5 were up-regulated in mature leaf and young blades, and roots, but exhibited low expression in leaf petioles. VvMKK2, VvMKK3, and VvMKK5 genes were differentially up-regulated when grapevine leaves were inoculated with spores of Erisyphe necator, or treated with salicylic acid (SA), ethylene (ETH), H2O2, or exposed to drought, indicating that these genes may be involved in a variety of signaling pathways. Over expression of VvMKK2 and VvMKK4 genes in transgenic Arabidopsis plants resulted in the production of seeds with a significantly higher germination and survival rate, and better seedling growth under stress conditions than wild-type plants. Overexpression of VvMKK2 in Arabidopsis improved salt and drought stress tolerance while overexpression of VvMKK4 only improved salt stress tolerance.

CONCLUSIONS

Results of the present investigation provide a better understanding of the interaction and function of MAPKKK-MAPKK-MAPK genes at the transcriptional level in grapevine and led to the identification of candidate genes for drought and salt stress in grapes.

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.

Genome-wide identification, expression profiles and regulatory network of MAPK cascade gene family in barley

BACKGROUND

Mitogen-activated protein kinase (MAPK) cascade is a conserved and universal signal transduction module in organisms. Although it has been well characterized in many plants, no systematic analysis has been conducted in barley.

RESULTS

Here, we identified 20 MAPKs, 6 MAPKKs and 156 MAPKKKs in barley through a genome-wide search against the updated reference genome. Then, phylogenetic relationship, gene structure and conserved protein motifs organization of them were systematically analyzed and results supported the predictions. Gene duplication analysis revealed that segmental and tandem duplication events contributed to the expansion of barley MAPK cascade genes and the duplicated gene pairs were found to undergone strong purifying selection. Expression profiles of them were further investigated in different organs and under diverse abiotic stresses using the available 173 RNA-seq datasets, and then the tissue-specific and stress-responsive candidates were found. Finally, co-expression regulatory network of MAPK cascade genes was constructed by WGCNA tool, resulting in a complicated network composed of a total of 72 branches containing 46 HvMAPK cascade genes and 46 miRNAs.

CONCLUSION

This study provides the targets for further functional study and also contribute to better understand the MAPK cascade regulatory network in barley and beyond.

Evolution, dynamic expression changes and regulatory characteristics of gene families involved in the glycerophosphate pathway of triglyceride synthesis in chicken (Gallus gallus)

It is well documented that four gene families, including the glycerophosphate acyltransferases (GPATs), acylglycerophosphate acyltransferases (AGPATs), lipid phosphate phosphohydrolases (LPINs) and diacylglycerol acyltransferases (DGATs), are involved in the glycerophosphate pathway of de novo triglyceride (TG) biosynthesis in mammals. However, no systematic analysis has been conducted to characterize the gene families in poultry. In this study, the sequences of gene family members in the glycerophosphate pathway were obtained by screening the public databases. The phylogenetic tree, gene structures and conserved motifs of the corresponding proteins were evaluated. Dynamic expression changes of the genes at different developmental stages were analyzed by qRT-PCR. The regulatory characteristics of the genes were analyzed by in vivo experiments. The results showed that the GPAT, AGPAT and LPIN gene families have 2, 7 and 2 members, respectively, and they were classified into 2, 4 and 2 cluster respectively based on phylogenetic analysis. All of the genes except AGPAT1 were extensively expressed in various tissues. Estrogen induction upregulated the expression of GPAM and AGPAT2, downregulated the expression of AGPAT3, AGPAT9, LPIN1 and LPIN2, and had no effect on the expression of the other genes. These findings provide a valuable resource for further investigation of lipid metabolism in liver of chicken.

Analysis of wheat gene expression related to the oxidative stress response and signal transduction under short-term osmotic stress

Water shortage is a major environmental stress that causes the generation of reactive oxygen species (ROS). The increase in ROS production induces molecular responses, which are key factors in determining the level of plant tolerance to stresses, including drought. The aim of this study was to determine the expression levels of genes encoding MAPKs (MAPK3 and MAPK6), antioxidant enzymes (CAT, APX and GPX) and enzymes involved in proline biosynthesis (P5CS and P5CR) in Triticum aestivum L. seedlings in response to short-term drought conditions. A series of wheat intervarietal substitution lines (ISCSLs) obtained by the substitution of single chromosomes from a drought-sensitive cultivar into the genetic background of a drought-tolerant cultivar was used. This source material allowed the chromosomal localization of the genetic elements involved in the response to the analyzed stress factor (drought). The results indicated that the initial plant response to drought stress resulted notably in changes in the expression of MAPK6 and CAT and both the P5CS and P5CR genes. Our results showed that the substitution of chromosomes 3B, 5A, 7B and 7D had the greatest impact on the expression level of all tested genes, which indicates that they contain genetic elements that have a significant function in controlling tolerance to water deficits in the wheat genome.

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.

Phylogenetic and expression analysis of histone acetyltransferases in Brachypodium distachyon

Histone acetylation is an important post-translational modification in eukaryotes and is regulated by two antagonistic enzymes, namely histone acetyltransferase (HAT) and histone deacetylase (HDAC). However, little has been done on the HAT superfamily in Brachypodium distachyon (B. distachyon), a new model plant of Poaceae. In this study, eight HATs were identified from B. distachyon and classified into four major families. Subcellular localization analysis showed that a majority of BdHATs were predominantly localized in the nucleus. Syntenic and phylogenetic analysis indicated there may be two common ancestral CREB-binding protein (p300/CBP, HAC) genes prior to the separation of monocots and dicots. Expression analysis revealed that the potential roles of BdHATs in B. distachyon development and responses to four abiotic stresses. Protein-protein network analysis identified some potential interactive genes with BdHATs. Thus, our results will provide solid basis for further study the function of HAT genes in B. distachyon and other monocot plants.

Genome-wide Identification of Jatropha curcas MAPK, MAPKK, and MAPKKK Gene Families and Their Expression Profile Under Cold Stress

Mitogen-activated protein kinase (MAPK) cascades are fundamental signal transduction modules in all eukaryotic organisms, controlling cell division, growth, development, and hormone signaling. Additionally, they can be activated in response to a variety of biotic and abiotic stressors. Although the evolution and expression patterns of MAPK cascade families have been systematically investigated in several model plants (e.g., Arabidopsis, rice, and poplar), we still know very little about MAPK, MAPKK, and MAPKKK families in Jatropha curcas, an economically important species. Therefore, this study performed genome-wide identification and transcriptional expression analysis of these three families in J. curcas. We identified 12 J. curcas MAPK (JcMAPKs), 5 JcMAPKKs, and 65 JcMAPKKKs. Phylogenetic analysis classified all JcMAPKs and JcMAPKKs into four subgroups, whereas JcMAPKKKs were grouped into three subfamilies (MEKK, RAF, and ZIK). Similarities in exon/intron structures supported the evolutionary relationships within subgroups and subfamilies. Conserved motif analysis indicated that all J. curcas MAPK cascades possessed typical, 200–300 amino-acid protein kinase domains. MAPK cascade genes were presented throughout all 11 chromosomes. Gene duplication analysis suggested that after JcMAPK and JcMAPKKK diverged, 3 and 19 tandem duplicates occurred under strong purifying selection. Furthermore, RNA-seq and qRT-PCR analyses revealed that some MAPK cascade genes are predominantly expressed in specific tissues. Moreover, their expression levels significantly increased under cold treatment. Our results should provide insight into the roles of MAPK cascade genes in regulating J. curcas stress responses and in hormonal signal transduction. Furthermore, these data have important applications in the genetic improvement of J. curcas.

Genome-Wide Bioinformatics Analysis of MAPK Gene Family in Kiwifruit (Actinidia Chinensis)

Mitogen activated protein kinase (MAPK) cascades are universal signal transduction modules that play crucial roles in various biotic and abiotic stresses, hormones, cell division, and developmental processes in plants. Mitogen activated protein kinase (MAPK/MPK), being a part of this cascade, performs an important function for further appropriate cellular responses. Although MAPKs have been investigated in several model plants, no systematic analysis has been conducted in kiwifruit (Actinidia chinensis). In the present study, we identified 18 putative MAPKs in the kiwifruit genome. This gene family was analyzed bioinformatically in terms of their chromosome locations, sequence alignment, gene structures, and phylogenetic and conserved motifs. All members possess fully canonical motif structures of MAPK. Phylogenetic analysis indicated that AcMAPKs could be classified into five subfamilies, and these gene motifs in the same group showed high similarity. Gene structure analysis demonstrated that the number of exons in AcMAPK genes ranged from 2 to 29, suggesting large variation among kiwifruit MAPK genes. The expression profiles of these AcMAPK genes were further investigated using quantitative real-time polymerase chain reaction (qRT-PCR), which demonstrated that AcMAPKs were induced or repressed by various biotic and abiotic stresses and hormone treatments, suggesting their potential roles in the biotic and abiotic stress response and various hormone signal transduction pathways in kiwifruit. The results of this study provide valuable insight into the putative physiological and biochemical functions of MAPK genes in kiwifruit.

Characterization of the APETALA2/Ethylene-responsive factor (AP2/ERF) transcription factor family in sunflower

One of the most prominent families of genes in plants is the AP2/ERF which play an important role in regulating plant growth and responses to various stresses. In this research, a genome-wide survey was conducted to recognize the AP2/ERF genes in sunflower (Helianthus annuus L.), and a total of 288 HaAP2/ERF was obtained. Phylogenetic analysis divided them into four sub-families, including 248 ERF, 4 RAV and 35 AP2, and one subgroup of the Soloist family. Localization of chromosome, gene structure, the conserved motif, gene ontology, interaction networks, homology modeling, the modeling of cis-regulatory elements and the analysis of events in the duplication of genes were carried out for HaAP2/ERF genes. Finally, 9AP2/ERF genes were chosen to confirm the gene expression of the selected genes in leaf and root tissues in various abiotic stress conditions by qPCR. The results confirmed that AP2/ERFs genes could effectively resist abiotic stress. Also, proline content was studied under drought, salinity, cold and heat stress. The results indicated that proline was increased under abiotic stress. This research has been done for the first time to determine the HaAP2/ERF family, which prepared valuable data for the evolutionary and practical research regarding AP2/ERF in sunflower.

Genome-Wide Analyses of the NAC Transcription Factor Gene Family in Pepper (Capsicum annuum L.): Chromosome Location, Phylogeny, Structure, Expression Patterns, Cis-Elements in the Promoter, and Interaction Network

The NAM, ATAF1/2, and CUC2 (NAC) transcription factors form a large plant-specific gene family, which is involved in the regulation of tissue development in response to biotic and abiotic stress. To date, there have been no comprehensive studies investigating chromosomal location, gene structure, gene phylogeny, conserved motifs, or gene expression of NAC in pepper (Capsicum annuum L.). The recent release of the complete genome sequence of pepper allowed us to perform a genome-wide investigation of Capsicum annuum L. NAC (CaNAC) proteins. In the present study, a comprehensive analysis of the CaNAC gene family in pepper was performed, and a total of 104 CaNAC genes were identified. Genome mapping analysis revealed that CaNAC genes were enriched on four chromosomes (chromosomes 1, 2, 3, and 6). In addition, phylogenetic analysis of the NAC domains from pepper, potato, Arabidopsis, and rice showed that CaNAC genes could be clustered into three groups (I, II, and III). Group III, which contained 24 CaNAC genes, was exclusive to the Solanaceae plant family. Gene structure and protein motif analyses showed that these genes were relatively conserved within each subgroup. The number of introns in CaNAC genes varied from 0 to 8, with 83 (78.9%) of CaNAC genes containing two or less introns. Promoter analysis confirmed that CaNAC genes are involved in pepper growth, development, and biotic or abiotic stress responses. Further, the expression of 22 selected CaNAC genes in response to seven different biotic and abiotic stresses [salt, heat shock, drought, Phytophthora capsici, abscisic acid, salicylic acid (SA), and methyl jasmonate (MeJA)] was evaluated by quantitative RT-PCR to determine their stress-related expression patterns. Several putative stress-responsive CaNAC genes, including CaNAC72 and CaNAC27, which are orthologs of the known stress-responsive Arabidopsis gene ANAC055 and potato gene StNAC30, respectively, were highly regulated by treatment with different types of stress. Our results also showed that CaNAC36 plays an important role in the interaction network, interacting with 48 genes. Most of these genes are in the mitogen-activated protein kinase (MAPK) family. Taken together, our results provide a platform for further studies to identify the biological functions of CaNAC genes.

Analysis of MAPK and MAPKK gene families in wheat and related Triticeae species

Background

The mitogen-activated protein kinase (MAPK) family is involved in signal transduction networks that underpin many different biological processes in plants, ranging from development to biotic and abiotic stress responses. To date this class of enzymes has received little attention in Triticeae species, which include important cereal crops (wheat, barley, rye and triticale) that represent over 20% of the total protein food-source worldwide.

Results

The work presented here focuses on two subfamilies of Triticeae MAPKs, the MAP kinases (MPKs), and the MAPK kinases (MKKs) whose members phosphorylate the MPKs. In silico analysis of multiple Triticeae sequence databases led to the identification of 152 MAPKs belonging to these two sub-families. Some previously identified MAPKs were renamed to reflect the literature consensus on MAPK nomenclature. Two novel MPKs, MPK24 and MPK25, have been identified, including the first example of a plant MPK carrying the TGY activation loop sequence common to mammalian p38 MPKs. An EF-hand calcium-binding domain was found in members of the Triticeae MPK17 clade, a feature that appears to be specific to Triticeae species. New insights into the novel MEY activation loop identified in MPK11s are offered. When the exon-intron patterns for some MPKs and MKKs of wheat, barley and ancestors of wheat were assembled based on transcript data in GenBank, they showed deviations from the same sequence predicted in Ensembl. The functional relevance of MAPKs as derived from patterns of gene expression, MPK activation and MKK-MPK interaction is discussed.

Conclusions

A comprehensive resource of accurately annotated and curated Triticeae MPK and MKK sequences has been created for wheat, barley, rye, triticale, and two ancestral wheat species, goat grass and red wild einkorn. The work we present here offers a central information resource that will resolve existing confusion in the literature and sustain expansion of MAPK research in the crucial Triticeae grains.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4545-9) contains supplementary material, which is available to authorized users.

Mitogen-Activated Protein Kinase Cascades in Plant Hormone Signaling

Mitogen-activated protein kinase (MAPK) modules play key roles in the transduction of environmental and developmental signals through phosphorylation of downstream signaling targets, including other kinases, enzymes, cytoskeletal proteins or transcription factors, in all eukaryotic cells. A typical MAPK cascade consists of at least three sequentially acting serine/threonine kinases, a MAP kinase kinase kinase (MAPKKK), a MAP kinase kinase (MAPKK) and finally, the MAP kinase (MAPK) itself, with each phosphorylating, and hence activating, the next kinase in the cascade. Recent advances in our understanding of hormone signaling pathways have led to the discovery of new regulatory systems. In particular, this research has revealed the emerging role of crosstalk between the protein components of various signaling pathways and the involvement of this crosstalk in multiple cellular processes. Here we provide an overview of current models and mechanisms of hormone signaling with a special emphasis on the role of MAPKs in cell signaling networks. One-sentence summary: In this review we highlight the mechanisms of crosstalk between MAPK cascades and plant hormone signaling pathways and summarize recent findings on MAPK regulation and function in various cellular processes.

Genome-wide identification and analysis of MAPK and MAPKK gene family in Chinese jujube (Ziziphus jujuba Mill.)

Background

Chinese jujube (Ziziphus jujuba Mill.) is one of the most important members in the Rhamnaceae family. The whole genome sequence and more than 30,000 proteins of Chinese jujube have been obtained in 2014. Mitogen-activated protein kinase cascades are universal signal transduction modules in plants, which is rapidly activated under various biotic and abiotic stresses. To date, there has been no comprehensive analysis of the MAPK and MAPKK gene family in Chinese jujube at the whole genome level.

Results

By performing a series of bioinformatics analysis, ten MAPK and five MAPKK genes were identified from the genome database of Chinese jujube, and then compared with the homologous genes from Arabidopsis. Phylogenetic analysis showed that ZjMAPKs was classified into four known groups, including A, B, C and D. ZjMAPKs contains five members of the TEY phosphorylation site and five members with the TDY motif. The ZjMAPKK family was subsequently divided into three groups, A, B and D. The gene structure, conserved motifs, functional annotation and chromosome distribution of ZjMAPKs and ZjMAPKKs were also predicted. ZjMAPKs and ZjMAPKKs were distributed on nine pseudo-chromosomes of Chinese jujube. Subsequently, expression analysis of ZjMAPK and ZjMAPKK genes using reverse transcription PCR and quantitative real-time PCR was carried out. The majority of ZjMAPK and ZjMAPKK genes were expressed in all tested organs/tissues with considerable differences in transcript levels indicating that they might be constitutively expressed. Moreover, ZjMKK5 was specific expressed in early development stage of jujube flower bud, indicating it plays some roles in reproductive organs development. The transcript expression of most ZjMAPK and ZjMAPKK genes was down-regulated in response to plant growth regulators, darkness treatment and phytoplasma infection.

Conclusions

We identified ten ZjMAPK and five ZjMAPKK genes from the genome database of Chinese jujube, the research results shown that ZjMPKs and ZjMKKs have the different expression patterns, indicating that they might play different roles in response to various treatments. The results provide valuable information for the further elucidation of physiological functions and biological roles of jujube MAPKs and MAPKKs.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4259-4) contains supplementary material, which is available to authorized users.

Genome-Wide Identification and Analysis of MAPK and MAPKK Gene Families in Bread Wheat (Triticum aestivum L.)

The mitogen-activated protein kinase (MAPK) cascade is a universal signal transduction module that plays a vital role in regulating growth and development, as well as environmental stress responses in plants. Wheat is one of the most important crops worldwide. Although the MAPK kinase kinase (MAP3K) family in wheat has been investigated, the MAPK and MAPK kinase (MAP2K) gene families remain unknown at present. Here, 54 MAPK and 18 MAPKK genes were identified in wheat using recent genomic information. Phylogenetic analysis of Triticum aestivum L. MAPKs and MAPKKs (TaMAPKs and TaMAPKKs) together with homologous genes from other species classified them into four groups, and the clustering was consistent with the genomic exon/intron structures. Conserved motif analysis found that MAPK proteins contained a typical TXY phosphorylation site and MAPKK proteins contained an S/T-X5-S/T motif. RNA-seq data mapping analysis showed that MAPK and MAPKK genes in group IV had tissue-specific expression profiles, whereas each group I member showed relatively high expression in all organs. Expression patterns of TaMAPK and TaMAPKK genes under stress conditions were also investigated and stress-responsive candidates were identified. Co-expression network analysis identified 11 TaMAPK genes and 6 TaMAPKK genes involved in the interaction network pathway. Overall, this study provided useful information for evolutionary and functional surveys of MAPK and MAPKK gene families in wheat and beyond.

Genome-wide analysis of basic helix-loop-helix (bHLH) transcription factors in Brachypodium distachyon

BACKGROUND

As a superfamily of transcription factors (TFs), the basic helix-loop-helix (bHLH) proteins have been characterized functionally in many plants with a vital role in the regulation of diverse biological processes including growth, development, response to various stresses, and so on. However, no systemic analysis of the bHLH TFs has been reported in Brachypodium distachyon, an emerging model plant in Poaceae.

RESULTS

A total of 146 bHLH TFs were identified in the Brachypodium distachyon genome and classified into 24 subfamilies. BdbHLHs in the same subfamily share similar protein motifs and gene structures. Gene duplication events showed a close relationship to rice, maize and sorghum, and segment duplications might play a key role in the expansion of this gene family. The amino acid sequence of the bHLH domains were quite conservative, especially Leu-27 and Leu-54. Based on the predicted binding activities, the BdbHLHs were divided into DNA binding and non-DNA binding types. According to the gene ontology (GO) analysis, BdbHLHs were speculated to function in homodimer or heterodimer manner. By integrating the available high throughput data in public database and results of quantitative RT-PCR, we found the expression profiles of BdbHLHs were different, implying their differentiated functions.

CONCLUSION

One hundred fourty-six BdbHLHs were identified and their conserved domains, sequence features, phylogenetic relationship, chromosomal distribution, GO annotations, gene structures, gene duplication and expression profiles were investigated. Our findings lay a foundation for further evolutionary and functional elucidation of BdbHLH genes.

Bioinformatics identification and transcript profile analysis of the mitogen-activated protein kinase gene family in the diploid woodland strawberry Fragaria vesca

Mitogen-activated protein kinases (MAPKs) play essential roles in mediating biotic and abiotic stress responses in plants. However, the MAPK gene family in strawberry has not been systematically characterized. Here, we performed a genome-wide survey and identified 12 MAPK genes in the Fragaria vesca genome. Protein domain analysis indicated that all FvMAPKs have typical protein kinase domains. Sequence alignments and phylogenetic analysis classified the FvMAPK genes into four different groups. Conserved motif and exon-intron organization supported the evolutionary relationships inferred from the phylogenetic analysis. Analysis of the stress-related cis-regulatory element in the promoters and subcellular localization predictions of FvMAPKs were also performed. Gene transcript profile analysis showed that the majority of the FvMAPK genes were ubiquitously transcribed in strawberry leaves after Podosphaera aphanis inoculation and after treatment with cold, heat, drought, salt and the exogenous hormones abscisic acid, ethephon, methyl jasmonate, and salicylic acid. RT-qPCR showed that six selected FvMAPK genes comprehensively responded to various stimuli. Additionally, interaction networks revealed that the crucial signaling transduction controlled by FvMAPKs may be involved in the biotic and abiotic stress responses. Our results may provide useful information for future research on the function of the MAPK gene family and the genetic improvement of strawberry resistance to environmental stresses.

CrMAPK3 regulates the expression of iron-deficiency-responsive genes in Chlamydomonas reinhardtii

BACKGROUND

Under iron-deficient conditions, Chlamydomonas exhibits high affinity for iron absorption. Nevertheless, the response, transmission, and regulation of downstream gene expression in algae cells have not to be investigated. Considering that the MAPK pathway is essential for abiotic stress responses, we determined whether this pathway is involved in iron deficiency signal transduction in Chlamydomonas.

RESULTS

Arabidopsis MAPK gene sequences were used as entry data to search for homologous genes in Chlamydomonas reinhardtii genome database to investigate the functions of mitogen-activated protein kinase (MAPK) gene family in C. reinhardtii under iron-free conditions. Results revealed 16 C. reinhardtii MAPK genes labeled CrMAPK2-CrMAPK17 with TXY conserved domains and low homology to MAPK in yeast, Arabidopsis, and humans. The expression levels of these genes were then analyzed through qRT-PCR and exposure to high salt (150 mM NaCl), low nitrogen, or iron-free conditions. The expression levels of these genes were also subjected to adverse stress conditions. The mRNA levels of CrMAPK2, CrMAPK3, CrMAPK4, CrMAPK5, CrMAPK6, CrMAPK8, CrMAPK9, and CrMAPK11 were remarkably upregulated under iron-deficient stress. The increase in CrMAPK3 expression was 43-fold greater than that in the control. An RNA interference vector was constructed and transformed into C. reinhardtii 2A38, an algal strain with an exogenous FOX1:ARS chimeric gene, to silence CrMAPK3. After this gene was silenced, the mRNA levels and ARS activities of FOX1:ARS chimeric gene and endogenous CrFOX1 were decreased. The mRNA levels of iron-responsive genes, such as CrNRAMP2, CrATX1, CrFTR1, and CrFEA1, were also remarkably reduced.

CONCLUSION

CrMAPK3 regulates the expression of iron-deficiency-responsive genes in C. reinhardtii.

StPBS2, a MAPK kinase gene, is involved in determining hyphal morphology, cell wall development, hypertonic stress reaction as well as the production of secondary metabolites in Northern Corn Leaf Blight pathogen Setosphaeria turcica

Mitogen activated protein kinase kinase (MAPKK) is a crucial component in the MAPK signaling pathway. However, the functions of MAPKKs in foliar pathogens remain poorly understood. In the current study, a MAPKK gene designated as StPBS2 was cloned from Setosphaeria turcica and the functions of this gene were investigated by RNAi technology. Four independent StPBS2 gene silence transformants with different efficiencies were confirmed by real time PCR. Compared to the wild type strain (WT), these transformants showed decreased colony growth, shortened hyphae cell length, broadened cell width and an obvious reduction in conidium yield. Moreover, the cell wall of the transformants was thicker and they were also more sensitive to substances that interfere with cell wall biosynthesis than WT. Additionally, the transformants displayed higher sensitivity to hypertonic stress than WT and the sensitivity was associated with the level of silencing of StPBS2. They were also resistant to the fungicides iprodione, procymidone and fludioxonil, to which WT almost completely sensitive. The transformants produced more red secondary metabolites than WT and the production was enhanced with increasing silencing level and increased glucose content in PDA medium. Our results suggest that StPBS2 is involved in morphogenesis, condiogenesis, cell wall development, hypertonic stress reaction and resistance to fungicides, as well as in the biosynthesis of secondary metabolites in S. turcica.

The MAPKKK and MAPKK gene families in banana: identification, phylogeny and expression during development, ripening and abiotic stress

The mitogen-activated protein kinase (MAPK) cascade, which is a major signal transduction pathway widely distributed in eukaryotes, has an important function in plant development and stress responses. However, less information is known regarding the MAPKKK and MAPKK gene families in the important fruit crop banana. In this study, 10 MAPKK and 77 MAPKKK genes were identified in the banana genome, and were classified into 4 and 3 subfamilies respectively based on phylogenetic analysis. Majority of MAPKKK and MAPKK genes in the same subfamily shared similar gene structures and conserved motifs. The comprehensive transcriptome analysis indicated that MAPKKK-MAPKK genes is involved in tissue development, fruit development and ripening, and response to abiotic stress of drought, cold and salt in two banana genotypes. Interaction networks and co-expression assays demonstrated that MAPK signaling cascade mediated network participates in multiple stress signaling, which was strongly activated in Fen Jiao (FJ). The findings of this study advance understanding of the intricately transcriptional control of MAPKKK-MAPKK genes and provide robust candidate genes for further genetic improvement of banana.

A Member of the 14-3-3 Gene Family in Brachypodium distachyon, BdGF14d, Confers Salt Tolerance in Transgenic Tobacco Plants

Plant 14-3-3 proteins are involved in diverse biological processes, but for the model monocotyledonous species, Brachypodium distachyon, their roles in abiotic stress tolerance are not well understood. In this study, a total of eight Bd14-3-3 genes were identified from B. distachyon and these were designated respectively as BdGF14a-BdGF14g. The qRT-PCR analyses of 3-month-old plants of B. distachyon showed that these genes were all expressed in the stems, leaves, and spikelets. By contrast, most of the plants had relatively lower transcriptional levels in their roots, except for the BdGF14g gene. The different expression profiles of the Bd14-3-3s under various stress treatments, and the diverse interaction patterns between Bd14-3-3s and BdAREB/ABFs, suggested that these gene products probably had a range of functions in the stress responses. The NaCl-induced Bd14-3-3 gene, BdGF14d, was selected for overexpression in tobacco. BdGF14d was found to be localized throughout the cell and it conferred enhanced tolerance to salt in the transgenic plants. Lowered contents of malondialdehyde, H₂O₂, and Na⁺, and lower relative electronic conductance (Rec%), yet greater activities of catalase and peroxidase, were observed in the overexpressing plants. Higher photosynthetic rate, transpiration rate, stomatal conductance, and water use efficiency were measured in the transgenic lines. Following abscisic acid (ABA) or NaCl treatment, stomatal aperture in leaves of the BdGF14d-overexpression plants was significantly lower than in leaves of the wild type (WT) controls. The stress-related marker genes involved in the ABA signaling pathway, the reactive oxygen species (ROS)-scavenging system, and the ion transporters were all up-regulated in the BdGF14d-overexpressing plants as compared with WT. Taken together, these results demonstrate that the Bd14-3-3 genes play important roles in abiotic stress tolerance. The ABA signaling pathway, the ROS-scavenging system, and ion transporters were all involved in enhancing the tolerance to salt stress in the BdGF14d-overexpression plants.

Divergent evolutionary patterns of the MAPK cascade genes in Brassica rapa and plant phylogenetics

Mitogen-activated protein kinase (MAPK) cascade signal transduction modules play crucial roles in regulating many biological processes in plants. These cascades are composed of three classes of hierarchically organized protein kinases, MAPKKKs, MAPKKs and MAPKs. Here, we analyzed gene retention, phylogenetic, evolution and expression patterns of MAPK cascade genes in Brassica rapa. We further found that the MAPK branches, classes III and IV, appeared after the split of bryophytes and green algae after analyzing the MAPK cascade genes in 8 species, and their rapid expansion led to the great size of the families of MAPKs. In contrast, the ancestral class I subfamily of MAPKK gene families have been highly conserved from algae to angiosperms. For the MAPKKK family, the MEKK and Raf subfamily share a common evolutionary origin, and Raf plays a major role in the expansion of the MAPKKK gene family. The cis-elements and interaction network analyses showed the important function of MAPK cascade genes in development and stress responses in B. rapa. This study provides a solid foundation for molecular evolution analyses of MAPK cascade genes.

Genome-wide identification, phylogeny and expression analysis of AP2/ERF transcription factors family in Brachypodium distachyon

Background

The AP2/ERF transcription factor is one of the most important gene families in plants, which plays the vital role in regulating plant growth and development as well as in response to diverse stresses. Although AP2/ERFs have been thoroughly characterized in many plant species, little is known about this family in the model plant Brachypodium distachyon, especially those involved in the regulatory network of stress processes.

Results

In this study, a comprehensive genome-wide search was performed to identify AP2/ERF gene family in Brachypodium and a total of 141 BdAP2/ERFs were obtained. Phylogenetic analysis classified them into four subfamilies, of which 112 belonged to ERF, four to RAV and 24 to AP2 as well as one to soloist subfamily respectively, which was in accordance with the number of AP2 domains and gene structure analysis. Chromosomal localization, gene structure, conserved protein motif and cis-regulatory elements as well as gene duplication events analysis were further performed to systematically investigate the evolutionary features of these BdAP2/ERF genes. Furthermore, the regulatory network between BdAP2/ERF and other genes were constructed using the orthology-based method, and 39 BdAP2/ERFs were found to be involved in the regulatory network and 517 network branches were identified. The expression profiles of BdAP2/ERF during development and under diverse stresses were investigated using the available RNA-seq and microarray data and ten tissue-specific and several stress-responsive BdAP2/ERF genes were identified. Finally, 11 AP2/ERF genes were selected to validate their expressions in different tissues and under different stress treatments using RT-PCR method and results verified that these AP2/ERFs were involved in various developmental and physiological processes.

Conclusions

This study for the first time reported the characteristics of the BdAP2/ERF family, which will provide the invaluable information for further evolutionary and functional studies of AP2/ERF in Brachypodium, and also contribute to better understanding the molecular basis for development and stresses tolerance in this model species and beyond.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2968-8) contains supplementary material, which is available to authorized users.

Bioinformatics Analysis of MAPKKK Family Genes in Medicago truncatula

Mitogen-activated protein kinase kinase kinase (MAPKKK) is a component of the MAPK cascade pathway that plays an important role in plant growth, development, and response to abiotic stress, the functions of which have been well characterized in several plant species, such as Arabidopsis, rice, and maize. In this study, we performed genome-wide and systemic bioinformatics analysis of MAPKKK family genes in Medicago truncatula. In total, there were 73 MAPKKK family members identified by search of homologs, and they were classified into three subfamilies, MEKK, ZIK, and RAF. Based on the genomic duplication function, 72 MtMAPKKK genes were located throughout all chromosomes, but they cluster in different chromosomes. Using microarray data and high-throughput sequencing-data, we assessed their expression profiles in growth and development processes; these results provided evidence for exploring their important functions in developmental regulation, especially in the nodulation process. Furthermore, we investigated their expression in abiotic stresses by RNA-seq, which confirmed their critical roles in signal transduction and regulation processes under stress. In summary, our genome-wide, systemic characterization and expressional analysis of MtMAPKKK genes will provide insights that will be useful for characterizing the molecular functions of these genes in M. truncatula.

Expansion and stress responses of AP2/EREBP superfamily in Brachypodium distachyon

APETALA2/ethylene-responsive element binding protein (AP2/EREBP) transcription factors constitute one of the largest and most conserved gene families in plant, and play essential roles in growth, development and stress response. Except a few members, the AP2/EREBP family has not been characterized in Brachypodium distachyon, a model plant of Poaceae. We performed a genome-wide study of this family in B. distachyon by phylogenetic analyses, transactivation assays and transcript profiling. A total of 149 AP2/EREBP genes were identified and divided into four subfamilies, i.e., ERF (ethylene responsive factor), DREB (dehydration responsive element binding gene), RAV (related to ABI3/VP) and AP2. Tandem duplication was a major force in expanding B. distachyon AP2/EREBP (BdAP2/EREBP) family. Despite a significant expansion, genomic organizations of BdAP2/EREBPs were monotonous as the majority of them, except those of AP2 subfamily, had no intron. An analysis of transcription activities of several closely related and duplicated BdDREB genes showed their functional divergence and redundancy in evolution. The expression of BdAP2/EREBPs in different tissues and the expression of DREB/ERF subfamilies in B. distachyon, wheat and rice under abiotic stresses were investigated by next-generation sequencing and microarray profiling. Our results are valuable for further function analysis of stress tolerant AP2/EREBP genes in B. distachyon.

Genome-Wide Identification, Evolution, and Co-expression Network Analysis of Mitogen-Activated Protein Kinase Kinase Kinases in Brachypodium distachyon

Mitogen-activated protein kinase (MAPK) cascades are the conserved and universal signal transduction modules in all eukaryotes, which play the vital roles in plant growth, development, and in response to multiple stresses. In this study, we used bioinformatics methods to identify 86 MAPKKK protein encoded by 73 MAPKKK genes in Brachypodium. Phylogenetic analysis of MAPKKK family from Arabidopsis, rice, and Brachypodium has classified them into three subfamilies, of which 28 belonged to MEKK, 52 to Raf, and 6 to ZIK subfamily, respectively. Conserved protein motif, exon-intron organization, and splicing intron phase in kinase domains supported the evolutionary relationships inferred from the phylogenetic analysis. And gene duplication analysis suggested the chromosomal segment duplication happened before the divergence of the rice and Brachypodium, while all of three tandem duplicated gene pairs happened after their divergence. We further demonstrated that the MAPKKKs have evolved under strong purifying selection, implying the conservation of them. The splicing transcripts expression analysis showed that the splicesome translating longest protein tended to be adopted. Furthermore, the expression analysis of BdMAPKKKs in different organs and development stages as well as heat, virus and drought stresses revealed that the MAPKKK genes were involved in various signaling pathways. And the circadian analysis suggested there were 41 MAPKKK genes in Brachypodium showing cycled expression in at least one condition, of which seven MAPKKK genes expressed in all conditions and the promoter analysis indicated these genes possessed many cis-acting regulatory elements involved in circadian and light response. Finally, the co-expression network of MAPK, MAPKK, and MAPKKK in Brachypodium was constructed using 144 microarray and RNA-seq datasets, and ten potential MAPK cascades pathway were predicted. To conclude, our study provided the important information for evolutionary and functional characterization of MAPKKK family in Brachypodium, which will facilitate the functional analysis of BdMAPKKK genes, and also will facilitate better understanding the MAPK signal pathway in Brachypodium and beyond.

Genome-Wide Identification and Expression Analysis of the Mitogen-Activated Protein Kinase Gene Family in Cassava

Mitogen-activated protein kinases (MAPKs) play central roles in plant developmental processes, hormone signaling transduction, and responses to abiotic stress. However, no data are currently available about the MAPK family in cassava, an important tropical crop. Herein, 21 MeMAPK genes were identified from cassava. Phylogenetic analysis indicated that MeMAPKs could be classified into four subfamilies. Gene structure analysis demonstrated that the number of introns in MeMAPK genes ranged from 1 to 10, suggesting large variation among cassava MAPK genes. Conserved motif analysis indicated that all MeMAPKs had typical protein kinase domains. Transcriptomic analysis suggested that MeMAPK genes showed differential expression patterns in distinct tissues and in response to drought stress between wild subspecies and cultivated varieties. Interaction networks and co-expression analyses revealed that crucial pathways controlled by MeMAPK networks may be involved in the differential response to drought stress in different accessions of cassava. Expression of nine selected MAPK genes showed that these genes could comprehensively respond to osmotic, salt, cold, oxidative stressors, and abscisic acid (ABA) signaling. These findings yield new insights into the transcriptional control of MAPK gene expression, provide an improved understanding of abiotic stress responses and signaling transduction in cassava, and lead to potential applications in the genetic improvement of cassava cultivars.