Original Research ARTICLE
Genome-wide identification and transcriptional expression analysis of mitogen-activated protein kinase and mitogen-activated protein kinase kinase genes in Capsicum annuum
- 1National Education Minster Key Laboratory of Plant Genetic Improvement and Comprehensive Utilization, Fujian Agriculture and Forestry University, Fuzhou, China
- 2College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou, China
- 3College of Crop Science, Fujian Agriculture and Forestry University, Fuzhou, China
- 4College of Life Science, Fujian Agriculture and Forestry University, Fuzhou, China
The tripartite mitogen-activated protein kinase (MAPK) signaling cascades have been implicated in plant growth, development, and environment adaptation, but a comprehensive understanding of MAPK signaling at genome-wide level is limited in Capsicum annuum. Herein, genome-wide identification and transcriptional expression analysis of MAPK and MAPK kinase (MAPKK) were performed in pepper. A total of 19 pepper MAPK (CaMAPKs) genes and five MAPKK (CaMAPKKs) genes were identified. Phylogenetic analysis indicated that CaMAPKs and CaMAPKKs could be classified into four groups and each group contains similar exon-intron structures. However, significant divergences were also found. Notably, five members of the pepper MAPKK family were much less conserved than those found in Arabidopsis, and 9 Arabidopsis MAPKs did not have orthologs in pepper. Additionally, 7 MAPKs in Arabidopsis had either two or three orthologs in the pepper genome, and six pepper MAPKs and one MAPKK differing in sequence were found in three pepper varieties. Quantitative real-time RT-PCR analysis showed that the majority of MAPK and MAPKK genes were ubiquitously expressed and transcriptionally modified in pepper leaves after treatments with heat, salt, and Ralstonia solanacearum inoculation as well as exogenously applied salicylic acid, methyl jasmonate, ethephon, and abscisic acid. The MAPKK-MAPK interactome was tested by yeast two-hybrid assay, the results showed that one MAPKK might interact with multiple MAPKs, one MAPK might also interact with more than one MAPKKs, constituting MAPK signaling networks which may collaborate in transmitting upstream signals into appropriate downstream cellular responses and processes. These results will facilitate future functional characterization of MAPK cascades in pepper.
Plant growth, development, and response to environmental cues are tightly and coordinately regulated by sophisticated signaling networks. Protein phosphorylation mediated by protein kinase cascades is one of the most important post-translational modifications in networks that efficiently transduce input signals into suitable outputs. One of the most important protein-kinase-based amplification cascades is the mitogen-activated protein kinase (MAPK) cascade, typically comprising the hierarchically organized MAPK kinase kinase (MAPKKK)-MAPK kinase (MAPKK)-MAPK, all encoded by multiple genes (Smékalová et al., 2014). In this cascade, MAPKs phosphorylate specific substrate proteins such as transcription factors and enzymes and subsequently trigger cellular responses. MAPKs are activated when both tyrosine and threonine residues in the TEY/TDY motif are phosphorylated by dual-specificity kinases (MAPKKs), which in turn are activated by MAPKKKs via phosphorylation of the conserved Thr/Ser motif (Group, 2002; Ishihama et al., 2011; Janitza et al., 2012). These cascades of activations rapidly transform upstream signals into appropriate intracellular responses (Group, 2002). Given the presence of multiple genes encoding MAPKKK (MAPKKKs), MAPKK (MAPKKs), and MAPK (MAPKs) family members in genomes of different plant species, one MAPK cascade may be associated with more than one upstream or downstream partner (Cardinale et al., 2002). Indeed, a myriad of MAPK cascade formations are possible, which may enable these modules to regulate a variety of different physiological and developmental processes. The signals perceived by MAPKKKs are integrated and transduced by MAPKKs, and the specificities of cellular responses regulated by MAPK modules are mainly determined by MAPKs. Extensive studies have demonstrated that multiple abiotic stimuli, such as wounding and exposure to drought, cold, salinity (Sinha et al., 2011; Zhang et al., 2014), pathogens, and microbe-associated molecular patterns (MAMPs); (Pitzschke et al., 2009), affect the expression of MAPKs, MAPKKs, and MAPKKKs. However, few of these genes have been associated with specific physiological functions. Limited numbers of MAPK cascades have been functionally identified and characterized. Some examples include drought, cold, and salt stress signal-associated cascades, such as AtMEKK1-AtMKK1/AtMKK2-AtMPK4 (Ichimura et al., 2000; Desikan et al., 2001; Hadiarto et al., 2006; Kong et al., 2012b; Zhang et al., 2012; Furuya et al., 2014), MEKK1-MKK4/5-MPK3/6-WRKY22/WRKY29 (a module involved in plant innate immunity); (Asai et al., 2002), AtMKK2-AtMPK10 MAPK (a module that regulates venation complexity by altering polar auxin transport efficiency); (Stanko et al., 2014), and OsMKK4-OsMPK1-OsWRKY53 (a module involved in the wounding signaling pathway); (Yoo et al., 2014), have been reported. Some MAPKs such as OsMPK20-4 and OsMPK3 might interact with each other without phosphorylation (Sheikh et al., 2013), and some MAPK cascades might crosstalk with CDPK signaling (Mehlmer et al., 2010). In some case, cascades in MAPK modules such as MEKK1 can take a shortcut to directly interact with the transcription factor WRKY53 (Miao et al., 2007). These crosstalks and shortcuts complicate the MAPK modules, which may facilitate adaptation and survival of plants in the face of diverse environmental stresses. The precise regulation mechanisms mediated by MAPK modules remain to be elucidated.
The completion of the genome sequencing of different plant species has made it possible to identify and characterize all members of the MAPK, MAPKK, and MAPKKK families. 20 MAPKs in Arabidopsis (Group, 2002), 17 MAPKs in rice (Hamel et al., 2005; Rao et al., 2010), 19 MAPKs in maize (Kong et al., 2013a), and 10 MAPKs in mulberry (Wei et al., 2014) have been identified. Depending on phylogeny, MAPKs from different plant species can be classified into four subgroups (groups A, B, C, and D); (Janitza et al., 2012; Asif et al., 2014). The TEY MAPKs constitute groups A–C; whereas, the TDY MAPKs are included in group D, and the estimated number in group D is diverse among different species with 8 in Arabidopsis (Cardinale et al., 2002), 10 in rice (Hamel et al., 2006; Rao et al., 2010), 14 in soybean (Neupane et al., 2013b), 11 in maize (Kong et al., 2013a; Liu et al., 2013), 9 in poplar (Hamel et al., 2005), and 9 in Brachypodium distachyon (Chen et al., 2012). Members in group A, B, C, and D (Sinha et al., 2011; Šamajová et al., 2013) have been implicated in responses of plants to abiotic and biotic stresses. Phylogenetic analyses of the MAPKK gene families showed that all MAPKK genes formed four groups (groups A, B, C, and D). The amino acid sequence of MAPKKs phosphorylation sites in plants is different from that of mammals. Plant MAPKKs contain conserved S/TxxxxxS/T activation motif, whereas mammalian enzymes have S/TxxxS/T motif. MAPK modules appear to be conserved throughout the evolution of higher plants; indeed, MAPK, MAPKK, and MAPKKK in different plant species shared sequence similarity (Janitza et al., 2012), and these evolutionary relationships within subfamilies are supported by exon-intron organizations and the architecture of conserved protein motifs. Sequence diversity also exists between members of MAPK and MAPKK family members in the same genome, as well as across genomes of different plant species (Chen et al., 2012; Wu et al., 2014; Wang et al., 2015). Genome-wide analysis of the MAPKs and MAPKKs might provide insights into the evolutionary relationships within subfamilies of MAPK and MAPKK and their potential roles in regulating growth, development, and responses to biotic and abiotic stresses in different plant species.
Pepper (Capsicum annuum), a vegetable crop of great agricultural and economic importance, is a typical Solanaceae species and is susceptible to soil-borne diseases. Heavy losses in pepper production are frequently caused by various soil-borne diseases as well as different abiotic stresses, including drought, extreme temperatures and salinity. The benefit of investigating the molecular mechanism of disease or stress resistance/tolerance is genetic improvement of the crop to better resist these diseases and stresses. Functional characterization of MAPK modules is a feasible approach for dissecting the molecular mechanism of crop disease or stress resistance/tolerance. So far, only two MAPKs have been characterized in pepper. Shin et al. cloned 2 mitogen-activated protein kinase (CaMK1 and CaMK2, renamed as CaMPK3 and CaMPK6-1 in our study) from pepper leaves and showed that CaMK1 and CaMK2 can be induced by wounding, UV-C and cold treatments (Shin et al., 2001). Furthermore, Huh et al. found that CaMK1 and CaMK2 can both be induced by Tobacco mosaic virus (TMV) infection and interact with CaWRKYa and phosphorylate CaWRKYa and proved that CaWRKYa is regulated by CaMK1 and CaMK2 at the posttranslational level in hot pepper (Huh et al., 2015), indicating that MAPK cascades may play important roles in the response of pepper to biotic and abiotic stresses. To get a more comprehensive understanding of the roles of MAPK cascades in pepper growth, development, and environment adaption, a genome-wide identification and expression analysis of MAPKs and MAPKKs was performed in the present study by searching the newly published pepper genome sequences (Kim et al., 2014). A total of 19 MAPK genes and five MAPKK were identified, annotated, and named according to their sequence similarity with Arabidopsis genes, and their chromosomal position and gene structure were also determined. Additionally, their transcript profiles in different organs were analyzed using published RNA-seq data, and their expression patterns in response to salt stress, heat shock, and Ralstonia solanacearum (RS) inoculation and exogenous application of phytohormones were also investigated using quantitative real time polymerase chain reaction (qPCR).
Our data showed that despite comparable structural conservation displayed by the conserved domains typically found in MAPKs and MAPKKs from other plant species, different degree of sequence divergences were also found in different MAPKs and MAPKKs between pepper and other plant species. Additionally, different transcriptional expressions of the members of MAPK and MAPKK families were found in different pepper organs and in the responses of pepper to the tested biotic and abiotic stresses, suggesting their important roles in the regulation of pepper growth, development, and responses to the tested biotic or abiotic stresses. To our knowledge, this is the first reported genome-wide and expression analysis of the pepper MAPK and MAPKK families, which would provide valuable information for further functional investigations of pepper MAPK and MAPKK families.
Materials and Methods
Identification of MAPK and MAPKK in Pepper
The completed genome sequences, CDS sequences, and peptide sequences of pepper version 1.55 were downloaded from the website (http://peppergenome.snu.ac.kr/) to construct a local genome database. To find all MAPK and MAPKK family members in pepper, NCBI BLASTP searches against a local database built using protein sequences were performed using sequences from all 35 known MAPKKs from Arabidopsis, tomato, rice, and soybean and 90 MAPKs from the above mentioned plant species (Table S1). In addition, a hidden Markov model (HMM) profile of these MAPK and MAPKK was generated using the Hmmer v 3.2 program (Finn et al., 2011) and was used as database to search similar proteins in the pepper protein database, sequences got from BLASTP search and analysis were further verified by HMM. Only the proteins with greater than 50% sequence identity were collected. Redundant sequences were removed manually, and the candidate sequences were further evaluated by identifying the putative functional domains using NCBI BLASTP (http://blast.ncbi.nlm.nih.gov/) and Arabidopsis information resource (http://www.arabidopsis.org/) and the domain analysis software SMART (Simple Modular Architecture Research Tool: http://smart.emblheidelberg.de/). Sequences without relevant domains or conserved motifs were removed. After multiple cycles of these analyses, MAPKK and MAPK family members in pepper were identified. All collected proteins were subjected to protein blast in the Arabidopsis information resource (http://arabidopsis.org/) for their nomenclature. The parameters of the deduced MAPK and MAPKK proteins, such as isoelectric points and molecular weights, were calculated on the ExPASy Proteomics Server (http://web.expasy.org/compute_pi/).
Analyses of Chromosomal Locations, Gene Structures, and Promoter Regions
Chromosomal locations and gene structures were obtained and analyzed from genome data downloaded from the pepper genome database (http://peppergenome.snu.ac.kr/ and http://peppersequence.genomics.cn/page/species/index.jsp). The promoter regions of the MAPK and MAPKK family members were identified from the whole-genome shotgun configurations available from the NCBI database. Cis-element analysis of the MAPK and MAPKK promoters were completed using the PlantCARE software (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/). The number of responsive cis-elements was calculated.
Sequence Homolog Alignment and Phylogenetic Analysis
To examine the evolutionary relationships between different MAPK and MAPKK members in pepper and Arabidopsis, the protein sequences of Arabidopsis and pepper MAPKs and MAPKKs were collected from the Arabidopsis information resource and pepper genome database (http://peppergenome.snu.ac.kr/), respectively. Phylogenetic analysis was performed using the Maximum likelihood method with bootstrapping analysis in the software MEGA version 6.0 (Kumar et al., 2008) and the bootstrap test was carried out with 1000 replicates.
Orthologs Sequence Analysis of MAPKS and MAPKKS among Pepper, Tomato, Arabidopsis, and Rice
To investigate the evolutionary relationships of MAPKs and MAPKKs among pepper, tomato, Arabidopsis and rice, the protein sequences of MAPKs and MAPKKs in pepper and their orthologs with highest sequence identities in tomato, Arabidopsis and rice were collected from the genome databases of pepper (http://peppergenome.snu.ac.kr/), tomato (http://solgenomics.net/tools/blast/), Arabidopsis (http://www.arabidopsis.org/), and rice (http://rice.plantbiology.msu.edu/analyses_search_blast.shtml), respectively. And the software DNAMAN was used to make multiple sequence alignment on MAPKs and MAPKKs of pepper with their orthologs in other plant species.
Plant Materials and Growth Conditions
The local cultivation pepper line 68-2 (from Quanzhou, Fujian, China) was used as our plant material. The pepper seeds were sowed on the soil for germination. After germination, the seedlings were transferred to the greenhouse and grown under standard conditions (16 h of light and 8 h of dark at 25°C), and 4-leaf-old pepper plants were used for treatment. For high-temperature treatment, the pepper plants were grown in a climate chamber at either high temperature (42°C) or normal temperature (25°C) for 24 h. For RS pathogen treatment, the plants were inoculated with a suspension of RS (cfu = 108) and maintained in a climate chamber at 28°C for 24 h. For plant hormones (ABA, SA, MeJA, and ETH) treatment, pepper plants were sprayed with either 100 μM abscisic acids (ABA), 1 mM salicylic acid (SA), 100 μM methyl jasmonate (MeJA), or 100 μM ethephon (ETH) dissolved in ddH2O. Meanwhile, control plants were sprayed with ddH2O. All plants were kept at 25°C for 24 h. Three repetitions were performed for each of the above groups. After treatment, the leaves were collected for RNA extraction. The pepper plant RNA-seq data of the two gene families (obtained from different tissues of pepper plants) were obtained from a previously published paper (Kim et al., 2014).
Quantitative Real-time PCR Analysis
The RNA-seq data of gene expression in different tissues in CM334 pepper plants was downloaded from Supplementary Table 22 of Kim's paper (Kim et al., 2014). Three individual plants were harvested for RNA extraction per treatment. The harvested leaves were frozen immediately by liquid nitrogen and transferred to the –80°C freezer. Frozen leaf tissues were disrupted in 2 mL RNAse-free microcentrifuge tubes (Axygen, USA) using three stainless steel beads and the Tissue Lyser II (Qiagen, Dusseldorf, Germany). Trizol (Invitrogen, Carlsbad, USA) was used to extract the total RNA from the disrupted tissues. Both RNA concentration and quality were confirmed by the NanoDrop 2000 (ThermoScientific, Massachusetts, USA). The RNA integrity was confirmed by agrarose electrophoresis. RNA (500 ng from each sample), 250 ng of oligo dT(15) primer, and 200 unit reverse transcriptase was used for a reverse transcription reaction using One Step PrimeScript™ cDNA Synthesis Kit (TaKaRa, Shigo, Japan), which include a genomic DNA digestion step, with the following programme: 42°C, 60 min; 85°C, 5 s; 4°C, forever. The synthesized cDNA products were diluted ten-fold for further qPCR analysis. qPCR was performed as a 20-μL reaction (10 μL 2 × SYBR Premix Ex Taq, 0.2 μM forward and reverse primers, and 100 ng cDNA) in a 200-μL PCR tube (Axygen, Hangzhou, China) using SYBR Premix Ex Taq™ (TaKaRa) and a CFX™ Real-Time System (Bio-rad, Hercules, USA) with the following program: 95°C, 30 s; 95°C, 5 s; 60°C, 34 s; 95°C, 15 s, 40 cycles. SYBR green fluorescence was used to detect the amplification of the target genes at each cycle. The Cq (threshold cycle), defined as the real-time PCR cycle at which a statistically significant increase of reporter fluorescence was first detected, was used as a measure for the starting copy numbers of the target gene. Three replicates of each experiment were performed, and three replicates were performed using independent plant individuals. Data were analyzed by the Livak method (Livak and Schmittgen, 2001; Bustin et al., 2009), and expressed as a normalized relative expression level (2−ΔΔCq) of the respective genes. CaActin (GQ339766) and 18S ribosomal RNA (EF564281), which have been confirmed by the published paper (Choi et al., 2012; Kim et al., 2015) and further confirmed by our own experiment, were used to normalize the transcript levels in pepper. Student's t-test was used to determine the statistical significance of the differential expression patterns between treatments and controls. Gene specific primers (GSPs) used in qPCR analyses are listed in Table S2 and the quality and specificity of the primers were confirmed by the software primer primer 5.0 and the genome database of pepper plants (http://peppergenome.snu.ac.kr/); (Table S1), respectively. The specificity of the PCR products was confirmed by the melting curve and the qPCR validation was also performed using the Bio-rad CFX manager 2.0 (Table S2).
Yeast Two-hybrid Assays
The interactions between MAPKs and MAPKKs proteins were identified in the Gal-4 reporter-based ProQuest™ two-hybrid system (Invitrogen, Darmstadt, Germany). CDS coding for all the MAPKs members of pepper were cloned into the prey vector (pDEST32) and CDS coding for MAPKKs were cloned into the bait vector (pDEST22). All the combinations of each bait and prey plasmid were PEG-transformed into the yeast strain MaV203 following the protocol of ProQuest™ two-hybrid system. Positive transformants were first selected in the synthetic dropout (SD) medium without leucine and tryptophan (SD/-Leu/-Trp) and the culture was transferred to the selection medium without Leucine, tryptophan, histidine and adenine. 3-amino-1,2,4-triazole (3-AT) was added to the selection plates to suppress the auto-activation of the prey vectors.
Identification of MAPK and MAPKK Genes in the Pepper Genome
The release of the complete pepper genome sequences makes it possible to identify pepper MAPK and MAPKK genes in a genome-wide manner (Kim et al., 2014). In the present study, the deduced amino acid sequences of MAPKs and MAPKKs from Arabidopsis, rice, tomato and soybean were used as queries to search for orthologs by BLASTP in the genome of pepper cultivar CM334, the redundant candidate sequences or sequences without relevant domains or conserved motifs were deleted manually. Sequences got from BLASTP search and analysis were further verified by HMM result. Finally, 19 CaMAPKs and five CaMAPKKs were acquired using this method. The protein sequences of the CaMAPKs and CaMAPKKs were further characterized by conserved domain analysis using the SMART software. The newly identified MAPK and MAPKK genes were designated following the nomenclature of Arabidopsis thaliana according to sequence similarities (Table 1). The number of MAPKs in pepper was similar to that of tomato (Kong et al., 2012a), Brachypodium distachyon (Chen et al., 2012), rice (Hamel et al., 2006; Rao et al., 2010), Arabidopsis (Asai et al., 2002; Colcombet and Hirt, 2008; Popescu et al., 2009), banana (Asif et al., 2014), and apple (Zhang et al., 2013); (Table 2). In contrast, the five MAPKKs in pepper were the same amount found in tomato (Wu et al., 2014), but differed significantly in number to the MAPKKs found in Arabidopsis, rice (Hamel et al., 2006; Rao et al., 2010), Brachypodium distachyon (Chen et al., 2012), and apple (Zhang et al., 2013). Sequence analysis revealed that several MAPKs in Arabidopsis, such as AtMPK2, 5, 8, 10, 11, 12, 14, 17, 18, had no orthologs in pepper; AtMPK4, 6, 9, 13, 15, 19, and 20 had more than one ortholog in the pepper genome. These orthologs in pepper were designated as CaMPK4-1, 4-2, 4-3, 6-1, 6-2, 9-1, 9-2, 13-1, 13-2, 17-1, 17-2, 19-1, 19-2, 20-1, and 20-2 based on their sequence similarities. Similarly, five MAPKKs had orthologs in the pepper genome, and these were designated as CaMKK2, –3, –5, –6 and –9 based on sequence similarities to their orthologs in Arabidopsis. The detailed information of two gene family members, including the number of introns, genome ID, genomic position, number of deduced amino acid residues, number of nucleotides, molecular weight, and isoelectric point (pI), is shown in Table 1. ORF (open reading frame) lengths of MAPK and MAPKK genes ranged from 984 to 1935 bp and 984 to 1140 bp, respectively, and the molecular weight of the resulting proteins ranged from 26.14 to 70.71 kDa and 36.76 to 41.72 kDa, respectively. The pIs of MAPKs and MAPKKs ranged from 5.06 to 9.54 and 5.68 to 8.89, respectively. The protein lengths of identified MAPKs in apple and maize were 370 amino acids (aa) to 865 aa and 369 aa to 615 aa, respectively and the corresponding pIs were 5.07–9.54 and 5.26–9.72 (Liu et al., 2013; Zhang et al., 2013). The lengths of identified MAPKKs in apple and maize were 314–554 aa and 227–613 aa and their pIs were 5.27–9.2 and 5.26–9.72, respectively (Liu et al., 2013; Zhang et al., 2013).
Multiple Alignments, Phylogenetic Analysis, and Domain Analysis of MAPK and MAPKK Families in the Pepper Genome
The alignments of the 19 MAPKs and five MAPKKs of pepper were performed using clustal X. The results of the alignment indicated that all 19 identified MAPKs contain 11 domains (I-XI) that are conserved in the serine/threonine protein kinases from other plant species (Figure 1). Several threonine and tyrosine residues located between domains VII and VIII of all 19 pepper MAPKs were conserved when compared to those of other plant species. These constitute the activation loop and are known to be phosphorylated during MAPK activation. Of the 19 MAPKs, nine contained the TDY motif, another nine contained the TEY motif, and only CaMPK4-3 contained the MEY motif. These results are consistent with that in rice and tomato, which each contains only one MAPK with a MEY domain (Zhang et al., 2013). Additionally, the majority of the TDY MAPKs had long C-terminal extensions compared to the TDY MAPKs in the pepper genome, which was consistent with the MAPKs from the other plant species (Zhang et al., 2013). The alignment of the five pepper MAPKKs showed the existence of 11 conserved motifs including S/T-X5-S/T (activation loop) and the DIK motif (active site); (Figure 2). These motifs were previously found in MAPKKs from other plant species (Hamel et al., 2006).
Figure 1. Multiple sequence alignment of the amino acid sequences of pepper MAPKs and MAPKs from other plant species (tomato, Arabidopsis, and rice). The alignment was performed using Clustal X. The 11 conserved domains (I–XI) present in pepper serine/threonine protein kinases are denoted by Roman numerals. The activation loop region is marked by a black line. The conserved threonine and tyrosine residues are indicated by asterisks.
Figure 2. Alignment of multiple pepper MAPKK amino acid sequences from pepper and other plant species (tomato, Arabidopsis, and rice). The alignment was performed using Clustal X. The 11 conserved domains (I–XI) present in pepper serine/threonine protein kinases are denoted by Roman numerals. The conserved S/T-X5-S/T motif and active site D (I/L/V) K motif are highlighted by a red box.
To analyze the evolutionary relationships between different MAPK and MAPKK family members in pepper and Arabidopsis, full-length amino acid sequences of 19 MAPKs and five MAPKKs from pepper, and 20 MAPKs and 10 MAPKKs from Arabidopsis, were together subjected to a multiple sequence alignment using the MEGA 6.0 program. An unrooted phylogenetic tree was constructed, and the 19 MAPKs were divided into four groups (A, B, C, and D) based on their predicted amino acid sequences (Figure 3). The size of groups A and B appeared to be similar; whereas, the size of group D differed significantly between pepper and Arabidopsis. The group D gene family is the largest subfamily in rice, poplar, tomato, and apple (Chen et al., 2012; Zhang et al., 2013; Wang et al., 2015). The five MAPKKs in the pepper genome were also classified into four groups with the MAPKKs of Arabidopsis. Like the MAPKKs in tomato (Wu et al., 2014), two MAPKKs were classified into group A, and groups B, C, and D each had one MAPKK in pepper (Figure 4). Similar to MAPKs, the number of MAPKKs in group D differed greatly between pepper and Arabidopsis.
Figure 3. Phylogenetic relationship among MAPK genes from pepper and Arabidopsis. The complete amino acid sequences of 19 pepper MAPKs and 20 Arabidopsis MAPKs were aligned using Clustal X, and the Maximum Likelihood tree was constructed using MEGA 6.0 with bootstrapping analysis (1000 replicates). The values above the branches are bootstrap support of 1000 replicates. Groups A–D indicate different gene clusters. To identify the plant species origin for each MAPK, a species acronym is included before the protein name: CaMAPK indicate MAPK from pepper and AtMAPK indicate MAPK from Arabidopsis. The green and red dots before the protein names indicate MAPKs from pepper and Arabidopsis, respectively.
Figure 4. Phylogenetic relationship among MAPKK genes from pepper and Arabidopsis. The complete amino acid sequences of five pepper MAPKKs combined with 10 Arabidopsis MAPKKs were aligned by Clustal X, and the Maximum Likelihood tree was constructed using MEGA 6.0 with bootstrapping analysis (1000 replicates). The values above the branches are bootstrap support of 1000 replicates. Groups A–D indicate different gene clusters. To identify the plant species origin for each MAPKK, a species acronym is included before the protein name: CaMKK indicate MAPKK from pepper and AtMKK indicate MAPKK from Arabidopsis. The green and red dots before the protein names indicate MAPKKs from pepper and Arabidopsis, respectively.
Genomic Distribution and Organization of MAPK and MAPKK Family Genes
The chromosomal locations of pepper MAPK and MAPKK genes were determined using BLASTN analysis of pepper chromosomes. The 19 CaMAPK genes were dispersed throughout 11 of the 12 chromosomes (Figure 5). Chromosomes 7 and 11 each contain 3 MAPK family genes, while chromosomes 1, 5, 6, and 8 each contain two MAPKs. Chromosomes 2, 3, 4, 10, and 12 each contain one MAPK. Among these MAPKs, two pairs of paralogous genes (CaMPK13-1/CaMPK13-2 and CaMPK20-1/CaMPK20-2) were located together on chromosomes 11 and 7, respectively. Other paralogs (CaMPK4-1/CaMPK4-2/CaMPK4-3, CaMPK6-1/CaMPK6-2, CaMPK17-1/CaMPK17-2, and CaMPK19-1/CaMPK19-2) were located on different chromosomes. Four of the five MAPKK family members, CaMKK3, CaMKK5, CaMKK6, and CaMKK9 were located on chromosome 3, while CaMKK2 was located on chromosome 2 (Figure 5). A similar observation was reported for tomato, where four MAPKKs were present on one chromosome and another MAPKK was located on a different chromosome (Wu et al., 2014).
Figure 5. Positions of MAPK and MAPKK gene family members on pepper chromosomes. The chromosome number is indicated at the top of each chromosome representation. Paralogs are indicated by green line. Scale represents a 50 Mb chromosomal distance. The paired MAPK genes are marked with a broken line.
To get insights into the gene structures of the MAPK and MAPKK gene families, the exon/intron organization of the two gene families was analyzed. MAPKs within the same groups exhibit similar exon-intron organizations. For example, MAPKs in group A contain 5–6 exons and 4–5 introns; MAPKs in group B contain 6 exons and 5 introns. MAPKs in group C contain 2 exons and one intron, and MAPKs in group D contain 7–10 exons and 6–9 introns (Figure 6). Similar to tomato, the MAPKKs in group A contain 8 exons and 7 introns. The group B MAPKK contains 6 exons and 5 introns, but MAPKK genes from groups C and D contain no introns (Figure 7).
Figure 6. The exon/intron structure of pepper MAPK genes. Introns and exons are represented by black lines and orange boxes, respectively. Introns and exons of the CaMAPKs are grouped according to the phylogenetic classification. The length in base pairs of each intron and exon is also indicated.
Figure 7. The exon/intron structure of pepper MAPKK genes. Introns and exons are represented by black lines and orange boxes, respectively. Introns and exons of the CaMAPKKs are grouped according to the phylogenetic classification. The length in base pairs of each intron and exon is also indicated.
Sequence Similarities of MAPKS and MAPKKS among Pepper, Tomato, Arabidopsis, and Rice
To investigate the evolutionary relationships of MAPKs and MAPKKs in pepper with their orthologs in other plant species, the sequences of MAPKs and MAPKKs in pepper genome were compared to their corresponding orthologs in tomato, Arabidopsis and rice (Table 3). Differences in the sequence similarities were found among different MAPK and MAPKK family members, for example, CaMPK4-3 and CaMKK9 both share more than 98.00% sequence identities with their orthologs in tomato, whereas, CaMPK16 and CaMKK3 both share only around 70% sequence identities with their orthologs in tomato. CaMPK7 and CaMKK6 share more than 86 and 73% sequence identities with their orthologs in rice, whereas, CaMPK17-2 and CaMKK3 share only around 50 and 16% with their orthologs in rice, respectively, suggesting that members in MAPK and MAPKK families in pepper exhibit different evolution rates. In addition, the sequence identities of orthologs in MAPK and MAPKK family between pepper and tomato are generally higher than those between pepper and Arabidopsis, and the sequence identities of orthologs of MAPK and MAPKK families between pepper and Arabidopsis are higher than those between pepper and rice. For example, the sequence identities of orthologs in MAPKs and MAPKKs between pepper and tomato range from 71.68 to 99.20% with an average of 89.70 and 69.32 to 98.25% with an average of 88.94%, respectively. The sequence identities of orthologs in MAPK and MAPKK families between pepper and Arabidopsis range from 54.56 to 87.12% (with an average of 73.57%) and 57.50 to 80.16% (with an average of 66.97%), respectively. While the sequence identities of orthologs in MAPK and MAPKK families between pepper and rice range from 50.00 to 78.17% (with an average of 68.81%) and 16.48 to 73.99% (with an average of 51.32%). But the sequence similarities of MAPK and MAPKK genes in pepper to their orthologs in Arabidopsis were not always higher than their sequence similarities to the corresponding orthologs in rice. For example, the sequence similarity of CaMPK1 and AtMPK1 is lower than that of CaMPK1 and OsMPK7. These results suggest that sequence similarities of orthologs in MAPK and MAPKK families in different plant species are dependent on the period of time that they have diverged.
Sequence Differences of MAPKS and MAPKKS among Different Pepper Cultivars
To study the evolution speed of MAPKs and MAPKKs in pepper genomes, the sequences of MAPKs and MAPKKs of CM334 released by Kim et al. (2014) were compared to those in two other cultivars, Chiltepin (Capsicum annuum var. glabriusculum) and Zunla-1 (Capsicum annuum L.), released by Qin et al. (2014); (Figure S1), as well as compared to the corresponding orthologs in Arabidopsis and tomato. The deduced amino acid sequences of the six MAPK genes, including CaMPK4-3, 17-1, 17-2, 19-1, 19-2, and 20-1 showed differences among the three pepper cultivars. For example, a fragment in the CDS of CaMPK4-3 encoding 20 amino acid residues was lost in CM334. In contrast, two fragments encoding 17 and 21 amino acid residues in the CDS of CaMPK19-1 were present in CM334 but were lost in Chiltepin and Zunla-1. Within the five MAPKKs, only CaMKK6 show differences among the three pepper cultivars.
Cis-Element Analysis in the Promoter Regions of MAPK and MAPKK Genes in Pepper
To investigate the possible roles of MAPKs and MAPKKs identified in the pepper genome, corresponding promoter regions (approximately 2 kb in length, the first nucleotide of the initiating ATG is designated +1) of the 19 MAPK and five MAPKK genes were subjected to cis-element analysis by PlantCARE online (Table 4). Cis-elements responsive to biotic [including W-box (Eulgem et al., 1999), EIRE (Fukuda, 1997), S box (pathogen-related cis-element Kunitomo et al., 2000)], and abiotic stresses [HSE (Mosser et al., 1988) and MBS (Nakagoshi et al., 1989)], as well as phytohormones such as SA [salicylic acid, TCA-element (Goldsbrough et al., 1993), JA jasmonic acid, CGTCA-motif (Beaudoin and Rothstein, 1997)], ET [ethylene, ERE (Itzhaki et al., 1994)], and ABA [abscisic acid, ABRE (Hamel et al., 2006)], were widely present in the promoters of the two gene families (Table 4). Several other cis-elements (listed in Table 4) were also found to be present only in the promoters of either MAPK genes or MAPKK genes, suggesting that transcriptional expression of MAPKK and MAPK genes may sometimes be co-regulated, but, at other times, regulated specifically.
MAPK and MAPKK Transcriptional Expression in Different Organs of Pepper Plants
The possible involvement of the MAPK modules in plant growth and development may be reflected in their expression patterns in different organs in pepper plants. To address this, the transcriptional expression profiles of the MAPKs and MAPKKs in different organs were analyzed to confirm their involvement in pepper development and growth using the published RNA-seq data of 6-week-old seedlings of a pepper variety (CM334), which were grown in a greenhouse on a 16-h light/8-h dark cycle (27°C/19°C). The majority of MAPK and MAPKK members expressed constitutively in leaves, roots, and stems (Figure 8). In leaves, the transcript level of CaMPK1 was the highest, followed by CaMPK19-2 and CaMPK6-2. In contrast, in roots, the transcript level of CaMPK19-2 was among the highest, followed by CaMPK4-1 and CaMPK20-1. In stems, CaMPK19-2 showed the highest transcript level, followed by CaMPK1 and CaMPK6-2. Among the five MAPKK genes, CaMKK5 exhibited the highest transcript levels in leaf and stem, whereas CaMKK9 exhibited the highest transcript level in roots. Taken together, these results indicate that MAPKs and MAPKKs are possibly involved in the growth or development in the tested pepper organs.
Figure 8. Expression levels of MAPKs and MAPKKs in different tissues in CM334 pepper plants (leave, stem, and root). The expression profiles of (A) MAPKs and (B) MAPKKs in different tissues of pepper plants. The pepper plant RNA-seq data of the two gene families were obtained from a previous study.
Expression Pattern Analysis of MAPKS and MAPKKS in Response to Phytohormones, Biotic, and Abiotic Stresses
MAPK modules have been implicated in plant responses to biotic and abiotic stresses, and the presence of stress-responsive cis-elements in the promoters of the MAPK and MAPKK genes suggest their involvement in pepper response to different stresses. To confirm this speculation, the transcriptional expression of the 19 MAPKs and five MAPKKs were analyzed in pepper plants by qPCR after application of salt stress, heat shock, and RS inoculation as well as exogenous application of SA, MeJA, ETH, and ABA (Figures 9, 10 and Figures S2–S8). Almost all of the MAPKK genes exhibited modified transcriptional expression upon these challenges, and were synergistically upregulated or downregulated by two of three stresses. Similarly, the expressions of almost all 19 MAPK genes were modified differently in level and timing by the three stresses and the four exogenously applied hormones, with the majority being upregulated and a few being downregulated. For the RS inoculation, the transcript levels of CaMPK1, CaMPK3, CaMPK4-1, 4-2, 4-3, 6-1, 6-2, 7, 9-1, 13-1, 13-2, and 17-1 were upregulated, however, the other MAPK members exhibited downregulated, except for CaMPK9-2 (no significant change). A total of 14 MAPK genes (CaMPK1, CaMPK3, CaMPK4-1, 4-3, 6-1, 6-2, 7, 9-1, 13-1, 13-2, 17-1, 17-2, 19-1, and 20-2) were significant upregulated and 4 MAPK genes were downregulated under the heat shock stress. Under the treatment of high salinity, mRNA level of 15 MAPK genes (CaMPK1, CaMPK3, CaMPK4-1, 4-3, 6-1, 6-2, 7, 9-1, 13-1, 13-2, 17-1, 17-2, 19-1, 20-1, and 20-2) exhibited significant upregulated and 4 MAPK genes were downregulated. The transcript levels of 11 MAPK genes (CaMPK1, CaMPK3, CaMPK4-1, 4-3, 6-1, 6-2, 7, 9-1, 13-1, 13-2, and 20-2) and 13 MAPK genes (CaMPK1, CaMPK3, CaMPK4-1, 4-3, 6-1, 6-2, 7, 9-1, 13-1, 13-2, 17-2, 19-2, and 20-2) can be significantly induced in response to exogenous application of SA and MeJA, respectively. Additionally, under the treatments of ET, the transcript levels of 16 MAPK genes (CaMPK1, CaMPK3, CaMPK4-1, 4-2, 4-3, 6-1, 6-2, 7, 9-1, 9-2, 13-1, 13-2, 17-2, 16, 19-2, and 20-2) can be significant induced, and the other 3 members showed significant downregulated. Only 2 MAPK genes (CaMPK9-1 and CaMPK17-1) can be significantly induced in exposure to ABA, but the other 17 members exhibited significant downregulated.
Figure 9. Heat maps showing the expression patterns of the MAPKs and MAPKKs in pepper in response to R.S (Ralstonia salanacearum) inoculation, high salinity (200 mM NaCl) and heat shock stresses (42°C). Quantitative RT-PCR (qPCR) was used to represent the transcript levels of MAPKs and MAPKKs in total RNA samples extracted from pepper leaves under the treatments of RS, high salinity and heat shock stresses. The relative expression (the mean fold changed between treated and control samples at each time points ± standard deviations) was log2 transformed. Each value represents the mean relative level of expression of three replicates. Blue represents low expression, red represents high expression. The heat map was generated using the software HemI.1.0.1.
Figure 10. Heat maps showing the expression patterns of the MAPKs and MAPKKs in pepper in response to phytohormones. The expression patterns of MAPKs and MAPKKs in response to ABA (100 μM), ETH (100 μM), MeJA (100 μM), and SA (1 mM). The relative expression (the mean fold changed between treated and control samples at each time points ± standard deviations) was log2 transformed. Each value represents the mean relative level of expression of three replicates. Blue represents low expression, red represents high expression. The heat map was generated using the software HemI.1.0.1.
The Interaction Pattern Identification of MAPKS and MAPKKS in Pepper Plant
It's well-known that the hierarchically organized MAP Kinase kinase kinase (MAPKKK)-MAP Kinase kinase (MAPKK)-MAP kinase, all encoded by multiple genes (Smékalová et al., 2014) constitute MAPK signaling cascades, in which MAPKs serve as substrates to MAPKK, and MAPKKs serve as substrates to MAPKKKs. The MAPK cascades play important roles in different plant biological processes ubiquitously. To determine MAPKK-MAPK patterns in pepper, we systematically assessed the interactions between all the MAPKs and MAPKKs using yeast two hybrid assay except for CaMPK17-2 and CaMPK19-2 that failed to be cloned possibly due to their low transcript levels (Figure 11). The positive interactions were first selected in the synthetic SD medium lacking leucine and tryptophan and the culture was transferred to the selection medium without Leucine, tryptophan, histidine and adenine with appropriate 3-amino-1,2,4-triazole. By this method, we got interacting MAPKs of only three of the five MAPKKs, (CaMKK3, CaMKK5, and CaMKK6) while no interacting MAPK was found using CaMKK2 and CaMKK9 as baits. Among three MAPKKs that we got interacting proteins, CaMKK3 interacted with CaMPK3; CaMKK5 interacted with CaMPK3, CaMPK6-1, CaMPK6-2, CaMPK4-1, CaMPK4-2, CaMPK4-3, CaMPK1, CaMPK7, CaMPK20-1, and CaMPK20-2; CaMKK6. To test if CaMKK5 can be autoactivated by itself, we performed another round of yeast two hybrid assay with CaMKK5 inserted into both of the prey and the bait vectors, the result showed that no positive clone grew on selection SD medium (Figure S9), suggesting that CaMKK5 can not be autoactivated by itself. The interaction of PcINF1 and SRC2-1 was used as a positive control (Liu et al., 2015).
Figure 11. Yeast two-hybrid (Y2H) analysis of MAPK-MAPKK interactions. Positive interactions are indicated by yeast growth on SD-His/-Ade selection media. The indicated concentrations of 3-amino-1,2,4-triazole (3-AT), were used to suppress the auto-activation background in the Y2H system. The pepper MAPKs are divided into 4 subgroups (A–D) on the basis of their phylogenetic tree (see Figure 3). The interaction of pDEST22 and pDEST32 was used as negative control.
MAPK modules have been implicated in plant growth, development, and adaption to environmental stresses. Recently, genome-wide analyses of gene families of MAPKs, MAPKKs, and MAPKKKs have been performed in several plant species; these analyses have provided insights into the mechanisms underlying the regulation of MAPK modules, laying an important foundation for further functional characterization (Liu et al., 2010; Yin et al., 2013; Zhang et al., 2013; Wu et al., 2014). The most noticeable features of identified MAPK and MAPKK families are conservation of their main motifs, classification into four conserved subgroups, and similar intron-exon organizations within the same families and across plant species (Chen et al., 2012; Zhang et al., 2013).
Evolution of Pepper MAPKS and MAPKKS Genes
In the present study, a total of 19 MAPKs and five MAPKKs were identified in the genome of the vegetable pepper, a typical Solanaceae with global agricultural and economic importance. These numbers were similar to the numbers of MAPKs and MAPKKs identified in other plant species, such as Arabidopsis (20/10), Oryza (17/8), Populus (21/11), tomato (16/5), apple (19/7), Medicago (18/4); (Ovečka et al., 2014), Phaseolus (15/9); (Neupane et al., 2013a), and soybean (38/11); (Neupane et al., 2013b). Consistent with previous results from other plant species, the 19 MAPKs in pepper possessed 11 conserved motifs, including the TEY, TDY, and MEY (only found in one MAPK) domains. The members of the MAPK and MAPKK gene families were both classified into four subfamilies, and the members within each family exhibited similar intron-exon structures. These data, together with previous studies, suggest similar origins and evolution patterns of the MAPK and MAPKK genes in different species. Possibly, the original family member duplicated prior to eudicot-monocot divergence (Zhang et al., 2013). Despite the sequence conservation, significant differences were found in the MAPK and MAPKK gene families in pepper compared to other plant species. Like tomato, pepper also possessed five MAPKKs, but this number is much less than that in Arabidopsis (10), rice (8), soybean (11), poplar (11), and Brachypodium distachyon (12). Although the 19 MAPK genes in pepper were similar to the 20 members in Arabidopsis, only a fraction of the MAPK genes in pepper have orthologs in Arabidopsis. In addition, several MAPK genes in Arabidopsis have more than one ortholog in the pepper genome. Perhaps, after the duplication of the origins of MAPK and MAPKK genes prior to eudicot-monocot divergence, some of the genes may have been lost during the evolution of pepper, while others might have duplicated along with the pepper genome expansion by rapid amplifications of retrotransposon elements (Kim et al., 2014). In fact, more than 81% of the pepper genome is composed of different transposable elements (TEs); (Kim et al., 2014). Additionally, differences in the sequences of MAPK and MAPKK orthologs among different plant species were also found in our study. Pepper shared much higher sequence similarities of MAPKs and MAPKKs with tomato than with Arabidopsis or rice, but the sequence similarities of MAPK and MAPKK genes in pepper to their orthologs in Arabidopsis were not always higher than to the corresponding orthologs in rice (For example, the sequence similarity between CaMPK1 and OsMPK7 is higher than that between CaMPK1 and AtMPK1), one explanation might be that some orthologs evolved slowly in the background of rice after eudicot-monocot divergence. Members of the MAPK and MAPKK families shared different amounts of sequence similarity to other plant species, indicating an imbalance in the evolutionary rates of different members of MAPK and MAPKK families, as some family members evolved more quickly than others. Furthermore, the sequences of six MAPKs were found to be different among the three pepper cultivars. Interestingly, these MAPK genes do not belong to the so-called quickly evolved genes, suggesting that mutation or evolution rates of these genes differ in different plant lineages. Further functional characterization of these divergent MAPKs in different pepper varieties might provide new insight into the mechanism underlying the regulation of MAPK modules in plant growth, development, and response to environmental stresses.
Expression Patterns of MAPKS and MAPKKS Genes in Pepper Tissues and in Response to Biotic and Abiotic Stresses and Exogenously Applied Hormones
Previous studies have established that plant's response to environment stress are largely regulated at the transcriptional level (Eulgem et al., 2004; Narsai et al., 2010), and the genes that transcriptionally modified against challenge of stress may play important roles in plant response to stress (Ramonell et al., 2005; Bartsch et al., 2006; Knoth et al., 2007). The RNA-seq data showed that the expression levels of the majority of pepper MAPKs and MAPKKs are relatively lower in leave, and higher in root and stem, similar observations that the majority of MAPK and MAPKK genes express preferentially in different tissues and developmental stages were also generally found (Zhang et al., 2013; Asif et al., 2014; Wei et al., 2014; Wu et al., 2014). These constitutively transcriptional expressions of MAPK and MAPKK genes suggest their possible involvement in the biological processes within the organs tested, but the precise roles and biological functions of their unequal expression levels remain to be elucidated in the future.
Under the challenges of heat shock, salt stress or RS inoculation, the majority of the MAPKs and MAPKKs in pepper genome were transcriptionally modified significantly with different timings and magnitudes in the present study. For example, within the 12 RS inoculation responsive MAPKs, CaMPK3, CaMPK4-1, CaMPK6-1, CaMPK6-2, CaMPK7, CaMPK9-1, and CaMPK13-1 were transciptionally upregulated by RS inoculation mainly at later time points, while CaMPK4-2, CaMPK17-1 and CaMPK19-1 were transcriptionally upregulated by RS inoculation at 3 and 6 h post inoculation. These results are consistent with that of the previous studies that MAPK and MAPKK genes are transcriptionally modified by various biotic or abiotic stresses (Liu et al., 2010; Chen et al., 2012; Yin et al., 2013; Zhang et al., 2013), for example, SpMPK3 (Li et al., 2014), MAP kinase 4 of barley (Abass and Morris, 2013), and GhMPK2 (Zhang et al., 2011a,b) were found to be involved in plant response to abiotic and biotic stresses. All these results indicate that MAPKs or MAPKKs in plants mainly act as positive regulators, only a small number act as negative regulators in response to biotic and abiotic stresses. Similar observations have been made in genomes of other plant species such as Brassica rapa (Lu et al., 2015), cucumber (Wang et al., 2015), Brachypodium distachyon (Jiang et al., 2015), maize (Sun et al., 2015), grapevine (Wang et al., 2014a), and tomato (Wu et al., 2014). By previous studies, a group of MAPK or MAPKK families' members may be required in the fine tune of plant response to a specific stress. For example, MPK9 and MPK12 (Jammes et al., 2011), MPK3, MPK4, and MPK6 (Pitzschke et al., 2009) in Arabidopsis function coordinately in defense reaction. However, the precise roles of these genes and how they coordinate in the response of plant to specific stress remain to be determined. In addition, despite the presence of genes specifically responding to a specific stress, members in MAPK and MAPKK families responding to two or three stresses were also found in the pepper genome in the present study, for example, CaMPK7 was synergistically induced by RS inoculation, heat stress and NaCl treatment, CaMPK4-1 was induced by heat stress, but was down-regulated by both RS inoculation and salinity, and CaMPK9-1 was upregulated by salt stress and RS inoculation, suggesting a diverse patterns of crosstalk between the tested three stresses mediated by different MAPK signaling modules. Consistently, the majority of mulberry MAPK family members were previously found to respond to different stresses (high/low temperature, salt, and drought); (Wei et al., 2014), and MAPKs or MAPKKs such as ZmMPK17 (Pan et al., 2012), GhMPK6a (Li et al., 2013), ZmMAPK1 (Wu et al., 2015), ZmSIMK1 (Wang et al., 2014b), ZmMKK1 (Cai et al., 2014), ZmMPK5 (Zhang et al., 2014) were also found to be involved in plant response to multiple stresses by individual studies, and some MAPKs were found to act as molecular linkers between abiotic stress signaling and plant reproduction(Sun et al., 2015). We speculate that MAPKs and MAPKKs might act as crucial nodes in crosstalk between plant response to biotic and abiotic stresses, which provide great potential in the fine tune of plant response to different stresses (Fujita et al., 2006; Schenke et al., 2011).
Since hormones such as SA, JA, ET, and ABA act as second messages ubiquitously in the response of plant to abiotic and biotic stresses (Seyfferth and Tsuda, 2014), and have been previously implicated in plant immunity, plant response to heat shock, salt stress, and pathogen infections, and crosstalk between plant response to biotic and abiotic stresses (Fujita et al., 2006; Clarke et al., 2009; Chan, 2012; van Zanten et al., 2012; Dang et al., 2014; Derksen et al., 2013; Hou et al., 2013; O'Brien and Benková, 2013; Singh and Jwa, 2013). In addition, MAPK cascades have been generally found to interact with signaling pathways mediated by hormones such as SA, JA, and ET (Smékalová et al., 2014), and hormone responsive-elements such as JA-responsive MeJA (Beaudoin and Rothstein, 1997), ethylene responsive-ERE (Itzhaki et al., 1994), and GCC-box (Ohme-Takagi and Shinshi, 1995), salicylic acid-responsive-TCA element (Goldsbrough et al., 1993), and abscisic acid-responsive ABRE (Guiltinan et al., 1990) were extensively found in the promoters of MAPK and MAPKK genes of pepper. To test the possible interactions between MAPK signaling in pepper and pathways mediated by hormones, the transcriptional modification of MAPKs and MAPKKs were determined by real-time PCR against the exogenous application of SA, MeJA, ETH, and ABA. Similarly to that challenged with the aforementioned abiotic stresses and RS inoculation, the majority of MAPKs and MAPKKs were found to be transcriptionally modified significantly by exogenously applied SA, MeJA, ETH, and ABA, and we got the similar results as the previous studies (Shin et al., 2001; Huh et al., 2015) that CaMPK3 (named as CaMK1 in the previous study) and CaMPK6-1 (named as CaMK2 in the previous study) can be induced in response to RS inoculation, high salinity and heat shock stresses, as well as exogenously applied SA, MeJA, ET, and ABA. We speculate that transcriptional modification of MAPK and MAPKK genes by heat shock, salt stress, and RS inoculation might correlate to signaling pathways mediated by SA, JA, ET, and ABA. Moreover, the MAPK modules may also act as components in hormone signal transduction (Singh and Jwa, 2013), and some of the MAPKs or MAPKKs responding to more than one exogenously applied hormones may act as convergent nodes in the integration of upstream signaling and in the orchestration of multiple cellular processes.
Patterns of Pepper MAPK and MAPKK Interactions and Their Possible Roles in Response of Pepper to Stresses
MAPK cascades, which consist of many distinct MAPKKK-MAPKK-MAPK modules linked to various upstream receptors and downstream targets through sequential phosphorylation and activation of the cascade components, are involved in different biological processes. Despite the occurrence of a large number of possible MAPK modules, no MAPK module has been experimentally identified in pepper so far. Yeast two hybrid system was generally employed to determine the MAPK interactome (Singh et al., 2012; Kong et al., 2013b). The MAPKK-MAPK interactome of pepper was assayed in this study, and the results showed that only CaMKK3, CaMKK5, and CaMKK6 got interacting MAPKs, one possibility may be that the other CaMKK2 and CaMKK9 might interacting with CaMPK17-2 and CaMPK19-2, which failed to be cloned possibly due to their low transcript level, another possibility can not be excluded that interacting MAPKs may not be completely found with the unphosphorylated MAPKKs as bait under the yeast two hybrid system, similarly to the results in maize (Kong et al., 2013b). Similarly to the previous studies in rice (Ding et al., 2009; Singh et al., 2012), our result in the present study showed that one MAPKKs may interact with multiple MAPKs. For example, the heat stress inducible CaMKK5 interacts significantly with CaMPK3, CaMPK6-1, CaMPK6-2, CaMPK4-1, CaMPK4-2, CaMPK4-3, CaMPK1, CaMPK7, CaMPK20-1, and CaMPK20-2. Similarly, CaMKK6 interacts significantly with CaMPK1, CaMPK20-1, and CaMPK20-2. Conversely, a single MAPK can interact with multiple upstream MAPKKs, for example, CaMPK1, CaMPK20-1 and CaMPK20-2 each interact with two MAPKKs, indicating that the MAPKK-MAPK modules constitute signaling networks. Our data also showed that the interaction of CaMKK5 to multiple MAPKs is not due to its autoactivation, to elucidate the underlying mechanism of the interactions of unphosphorylated MAPKKs, for example CaMKK5, with multiple MAPKs, further investigations are needed. Since both MAPKKs and MAPKs in pepper were upregulated or downregulated transcriptionally by heat, salt stress and RS inoculation as well as by exogenously applied SA, MeJA, ETH, and ABA with different timing, magnitude, these modules may collaborate in transmitting upstream signals into appropriate downstream cellular responses and processes.
Collectively, a total of 19 pepper MAPK genes and five MAPKK genes were identified in the present study. Compared to MAPK and MAPKK families in other plant species, comparable conservations and significant divergences among pepper and other plant species were found in the MAPK and MAPKK gene families. The majority of MAPKK and MAPK genes were expressed preferentially in all tested organs and were also transcriptionally modified in leaves by heat shock, salt stress, and RS inoculation, as well as the exogenous application of SA, MeJA, ETH, and ABA. By real time PCR, one specific stimulus simultaneously modified the expression of several MAPKs and MAPKKs, and one MAPK or MAPKK gene might respond synergistically to several stimuli. In addition, by interactions, MAPKKs and MAPKs constitute modules and signaling network, which might have important roles in orchestrating pepper growth, development, and response to specific biotic or abiotic stress. These results will facilitate the future functional characterization of MAPK cascades in pepper.
Conceived and designed the experiments: ZL and SH. Performed the experiments: ZL, LS, YL, LS, JC, QT, and JW. Analyzed the data: ZL, SY, XH, YL, YL, JH, CL, SM, and DG. Wrote the paper: ZL and SH. Provided guidance on the whole study: ZL and SH.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by grants from the National Natural Science Foundation of China (30971718, 31301254, and 31372061), and the Natural Science Foundation of Fujian Province, China (2014J01077). The funders had no role in study design, data collection, data analysis, decision to publish, or preparation of the manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.00780
Table S1. Nomenclature for MAPKs and MAPKKs in Arabidopsis, tomato, rice, and soybean.
Table S2. Primers used in this study and the qPCR validation.
Figure S1. Multiple sequence alignment of the amino acid sequences of pepper MAPKs in three different pepper cultivars.
Figure S2. The expression pattern of the MAPKs and MAPKKs in pepper in response to NaCl. (A) The expression pattern of MAPKs in response to NaCl. (B) The expression pattern of MAPKKs in response to NaCl. Asterisks indicate a significant difference as determined by Fisher's protected LSD test (P < 0.05).
Figure S3. The expression pattern of the MAPKs and MAPKKs in pepper in response to heat shock. (A) The expression pattern of MAPKs in response to heat shock. (B) The expression pattern of MAPKKs in response to heat shock. Asterisks indicate a significant difference as determined by Fisher's protected LSD test (P < 0.05).
Figure S4. The expression pattern of the MAPKs and MAPKKs in pepper in response to pathogen infection. (A) The expression pattern of MAPKs in response to pathogen infection. (B) The expression pattern of MAPKKs in response to pathogen infection. Asterisks indicate a significant difference as determined by Fisher's protected LSD test (P < 0.05).
Figure S5. The expression pattern of the MAPKs and MAPKKs in pepper in response to exogenous SA. (A) The expression pattern of MAPKs in response to exogenous SA. (B) The expression pattern of MAPKKs in response to exogenous SA. Asterisks indicate a significant difference as determined by Fisher's protected LSD test (P < 0.05).
Figure S6. The expression pattern of the MAPKs and MAPKKs in pepper in response to exogenous MeJA. (A) The expression pattern of MAPKs in response to exogenous MeJA. (B) The expression pattern of MAPKKs in response to exogenous MeJA. Asterisks indicate a significant difference as determined by Fisher's protected LSD test (P < 0.05).
Figure S7. The expression pattern of the MAPKs and MAPKKs in pepper in response to exogenous ET. (A) The expression pattern of MAPKs in response to exogenous ET. (B) The expression pattern of MAPKKs in response to exogenous ET. Asterisks indicate a significant difference as determined by Fisher's protected LSD test (P < 0.05).
Figure S8. The expression pattern of the MAPKs and MAPKKs in pepper in response to exogenous ABA. (A) The expression pattern of MAPKs in response to exogenously ABA. (B) The expression pattern of MAPKKs in response to exogenously ABA. Asterisks indicate a significant difference as determined by Fisher's protected LSD test (P < 0.05).
Figure S9. The autoactivation identification of CaMKK5 and CaMKK5 by yeast two hybird. Yeast cells were plated on SD/–Leu/–Trp/–His medium containing 100mM 3-amino-1,2,4-triazole. The interaction of PcINF1 and SRC2-1 was used as positive control.
Abass, M., and Morris, P. C. (2013). The Hordeum vulgare signalling protein MAP kinase 4 is a regulator of biotic and abiotic stress responses. J. Plant Physiol. 170, 1353–1359. doi: 10.1016/j.jplph.2013.04.009
Asai, T., Tena, G., Plotnikova, J., Willmann, M. R., Chiu, W. L., Gomez-Gomez, L., et al. (2002). MAP kinase signalling cascade in Arabidopsis innate immunity. Nature 415, 977–983. doi: 10.1038/415977a
Asif, M. H., Lakhwani, D., Pathak, S., Bhambhani, S., Bag, S. K., and Trivedi, P. K. (2014). Genome-wide identification and expression analysis of the mitogen-activated protein kinase gene family from banana suggest involvement of specific members in different stages of fruit ripening. Funct. Integr. Genomics 14, 161–175. doi: 10.1007/s10142-013-0349-9
Bartsch, M., Gobbato, E., Bednarek, P., Debey, S., Schultze, J. L., Bautor, J., et al. (2006). Salicylic acid-independent ENHANCED DISEASE SUSCEPTIBILITY1 signaling in Arabidopsis immunity and cell death is regulated by the monooxygenase FMO1 and the Nudix hydrolase NUDT7. Plant Cell 18, 1038–1051. doi: 10.1105/tpc.105.039982
Bustin, S. A., Benes, V., Garson, J. A., Hellemans, J., Huggett, J., Kubista, M., et al. (2009). The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 55, 611–622. doi: 10.1373/clinchem.2008.112797
Cai, G., Wang, G., Wang, L., Pan, J., Liu, Y., and Li, D. (2014). ZmMKK1, a novel group A mitogen-activated protein kinase kinase gene in maize, conferred chilling stress tolerance and was involved in pathogen defense in transgenic tobacco. Plant Sci. 214, 57–73. doi: 10.1016/j.plantsci.2013.09.014
Cardinale, F., Meskiene, I., Ouaked, F., and Hirt, H. (2002). Convergence and divergence of stress-induced mitogen-activated protein kinase signaling pathways at the level of two distinct mitogen-activated protein kinase kinases. Plant Cell 14, 703–711. doi: 10.1105/tpc.010256
Chan, Z. (2012). Expression profiling of ABA pathway transcripts indicates crosstalk between abiotic and biotic stress responses in Arabidopsis. Genomics 100, 110–115. doi: 10.1016/j.ygeno.2012.06.004
Chen, L., Hu, W., Tan, S., Wang, M., Ma, Z., Zhou, S., et al. (2012). Genome-wide identification and analysis of MAPK and MAPKK gene families in Brachypodium distachyon. PLoS ONE 7:e46744. doi: 10.1371/journal.pone.0046744
Choi, D. S., Hwang, I. S., and Hwang, B. K. (2012). Requirement of the cytosolic interaction between PATHOGENESIS-RELATED PROTEIN10 and LEUCINE-RICH REPEAT PROTEIN1 for cell death and defense signaling in pepper. Plant Cell 24, 1675–1690. doi: 10.1105/tpc.112.095869
Clarke, S. M., Cristescu, S. M., Miersch, O., Harren, F. J., Wasternack, C., and Mur, L. A. (2009). Jasmonates act with salicylic acid to confer basal thermotolerance in Arabidopsis thaliana. New Phytol. 182, 175–187. doi: 10.1111/j.1469-8137.2008.02735.x
Dang, F., Wang, Y., She, J., Lei, Y., Liu, Z., Eulgem, T., et al. (2014). Overexpression of CaWRKY27, a subgroup IIe WRKY transcription factor of Capsicum annuum, positively regulates tobacco resistance to Ralstonia solanacearum infection. Physiol. Plant. 150, 397–411. doi: 10.1111/ppl.12093
Desikan, R., Hancock, J. T., Ichimura, K., Shinozaki, K., and Neill, S. J. (2001). Harpin induces activation of the Arabidopsis mitogen-activated protein kinases AtMPK4 and AtMPK6. Plant Physiol. 126, 1579–1587. doi: 10.1104/pp.126.4.1579
Eulgem, T., Rushton, P. J., Schmelzer, E., Hahlbrock, K., and Somssich, I. E. (1999). Early nuclear events in plant defence signalling: rapid gene activation by WRKY transcription factors. EMBO J. 18, 4689–4699. doi: 10.1093/emboj/18.17.4689
Eulgem, T., Weigman, V. J., Chang, H. S., McDowell, J. M., Holub, E. B., Glazebrook, J., et al. (2004). Gene expression signatures from three genetically separable resistance gene signaling pathways for downy mildew resistance. Plant Physiol. 135, 1129–1144. doi: 10.1104/pp.104.040444
Fujita, M., Fujita, Y., Noutoshi, Y., Takahashi, F., Narusaka, Y., Yamaguchi-Shinozaki, K., et al. (2006). Crosstalk between abiotic and biotic stress responses: a current view from the points of convergence in the stress signaling networks. Curr. Opin. Plant Biol. 9, 436–442. doi: 10.1016/j.pbi.2006.05.014
Fukuda, Y. (1997). Interaction of tobacco nuclear protein with an elicitor-responsive element in the promoter of a basic class I chitinase gene. Plant Mol. Biol. 34, 81–87. doi: 10.1023/A:1005737128339
Furuya, T., Matsuoka, D., and Nanmori, T. (2014). Membrane rigidification functions upstream of the MEKK1-MKK2-MPK4 cascade during cold acclimation in Arabidopsis thaliana. FEBS Lett. 588, 2025–2030. doi: 10.1016/j.febslet.2014.04.032
Goldsbrough, A. P., Albrecht, H., and Stratford, R. (1993). Salicylic acid-inducible binding of a tobacco nuclear protein to a 10 bp sequence which is highly conserved amongst stress-inducible genes. Plant J. 3, 563–571. doi: 10.1046/j.1365-313X.1993.03040563.x
Hadiarto, T., Nanmori, T., Matsuoka, D., Iwasaki, T., Sato, K., Fukami, Y., et al. (2006). Activation of Arabidopsis MAPK kinase kinase (AtMEKK1) and induction of AtMEKK1-AtMEK1 pathway by wounding. Planta 223, 708–713. doi: 10.1007/s00425-005-0126-7
Hamel, L. P., Miles, G. P., Samuel, M. A., Ellis, B. E., Séguin, A., and Beaudoin, N. (2005). Activation of stress-responsive mitogen-activated protein kinase pathways in hybrid poplar (Populus trichocarpa x Populus deltoides). Tree Physiol. 25, 277–288. doi: 10.1093/treephys/25.3.277
Hamel, L. P., Nicole, M. C., Sritubtim, S., Morency, M. J., Ellis, M., Ehlting, J., et al. (2006). Ancient signals: comparative genomics of plant MAPK and MAPKK gene families. Trends Plant Sci. 11, 192–198. doi: 10.1016/j.tplants.2006.02.007
Huh, S. U., Lee, G. J., Jung, J. H., Kim, Y., Kim, Y. J., and Paek, K. H. (2015). Capsicum annuum transcription factor WRKYa positively regulates defense response upon TMV infection and is a substrate of CaMK1 and CaMK2. Sci. Rep. 5, 7981. doi: 10.1038/srep07981
Ichimura, K., Mizoguchi, T., Yoshida, R., Yuasa, T., and Shinozaki, K. (2000). Various abiotic stresses rapidly activate Arabidopsis MAP kinases ATMPK4 and ATMPK6. Plant J. 24, 655–665. doi: 10.1046/j.1365-313x.2000.00913.x
Ishihama, N., Yamada, R., Yoshioka, M., Katou, S., and Yoshioka, H. (2011). Phosphorylation of the Nicotiana benthamiana WRKY8 transcription factor by MAPK functions in the defense response. Plant Cell 23, 1153–1170. doi: 10.1105/tpc.110.081794
Itzhaki, H., Maxson, J. M., and Woodson, W. R. (1994). An ethylene-responsive enhancer element is involved in the senescence-related expression of the carnation glutathione-S-transferase (GST1) gene. Proc. Natl. Acad. Sci. U.S.A. 91, 8925–8929. doi: 10.1073/pnas.91.19.8925
Jammes, F., Yang, X., Xiao, S., and Kwak, J. M. (2011). Two Arabidopsis guard cell-preferential MAPK genes, MPK9 and MPK12, function in biotic stress response. Plant Signal. Behav. 6, 1875–1877. doi: 10.4161/psb.6.11.17933
Janitza, P., Ullrich, K. K., and Quint, M. (2012). Toward a comprehensive phylogenetic reconstruction of the evolutionary history of mitogen-activated protein kinases in the plant kingdom. Front. Plant Sci. 3:271. doi: 10.3389/fpls.2012.00271
Jiang, M., Wen, F., Cao, J., Li, P., She, J., and Chu, Z. (2015). Genome-wide exploration of the molecular evolution and regulatory network of mitogen-activated protein kinase cascades upon multiple stresses in Brachypodium distachyon. BMC Genomics 16, 228. doi: 10.1186/s12864-015-1452-1
Kim, D. S., Kim, N. H., and Hwang, B. K. (2015). The Capsicum annuum class IV chitinase ChitIV interacts with receptor-like cytoplasmic protein kinase PIK1 to accelerate PIK1-triggered cell death and defence responses. J. Exp. Bot. 66, 1987–1999. doi: 10.1093/jxb/erv001
Kim, S., Park, M., Yeom, S. I., Kim, Y. M., Lee, J. M., Lee, H. A., et al. (2014). Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat. Genet. 46, 270–278. doi: 10.1038/ng.2877
Knoth, C., Ringler, J., Dangl, J. L., and Eulgem, T. (2007). Arabidopsis WRKY70 is required for full RPP4-mediated disease resistance and basal defense against Hyaloperonospora parasitica. Mol. Plant Microbe Interact. 20, 120–128. doi: 10.1094/MPMI-20-2-0120
Kong, F., Wang, J., Cheng, L., Liu, S., Wu, J., Peng, Z., et al. (2012a). Genome-wide analysis of the mitogen-activated protein kinase gene family in Solanum lycopersicum. Gene 499, 108–120. doi: 10.1016/j.gene.2012.01.048
Kong, Q., Qu, N., Gao, M., Zhang, Z., Ding, X., Yang, F., et al. (2012b). The MEKK1-MKK1/MKK2-MPK4 kinase cascade negatively regulates immunity mediated by a mitogen-activated protein kinase kinase kinase in Arabidopsis. Plant Cell 24, 2225–2236. doi: 10.1105/tpc.112.097253
Kong, X., Lv, W., Zhang, D., Jiang, S., Zhang, S., and Li, D. (2013a). Genome-wide identification and analysis of expression profiles of maize mitogen-activated protein kinase kinase kinase. PLoS ONE 8:e57714. doi: 10.1371/journal.pone.0057714
Kong, X., Pan, J., Zhang, D., Jiang, S., Cai, G., Wang, L., et al. (2013b). Identification of mitogen-activated protein kinase kinase gene family and MKK-MAPK interaction network in maize. Biochem. Biophys. Res. Commun. 441, 964–969. doi: 10.1016/j.bbrc.2013.11.008
Kumar, S., Nei, M., Dudley, J., and Tamura, K. (2008). MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief. Bioinformatics 9, 299–306. doi: 10.1093/bib/bbn017
Kunitomo, H., Higuchi, T., Iino, Y., and Yamamoto, M. (2000). A zinc-finger protein, Rst2p, regulates transcription of the fission yeast ste11+ gene, which encodes a pivotal transcription factor for sexual development. Mol. Biol. Cell 11, 3205–3217. doi: 10.1091/mbc.11.9.3205
Li, C., Chang, P. P., Ghebremariam, K. M., Qin, L., and Liang, Y. (2014). Overexpression of tomato SpMPK3 gene in Arabidopsis enhances the osmotic tolerance. Biochem. Biophys. Res. Commun. 443, 357–362. doi: 10.1016/j.bbrc.2013.11.061
Li, Y., Zhang, L., Wang, X., Zhang, W., Hao, L., Chu, X., et al. (2013). Cotton GhMPK6a negatively regulates osmotic tolerance and bacterial infection in transgenic Nicotiana benthamiana, and plays a pivotal role in development. FEBS J. 280, 5128–5144. doi: 10.1111/febs.12488
Liu, N., Enkemann, S. A., Liang, P., Hersmus, R., Zanazzi, C., Huang, J., et al. (2010). Genome-wide gene expression profiling reveals aberrant MAPK and Wnt signaling pathways associated with early parthenogenesis. J. Mol. Cell Biol. 2, 333–344. doi: 10.1093/jmcb/mjq029
Liu, Y. K., Zhang, D., Wang, L., and Li, D. Q. (2013). Genome-wide analysis of mitogen-activated protein kinase gene family in maize. Plant Mol. Biol. Rep. 31, 1446–1460. doi: 10.1007/s11105-013-0623-y
Liu, Z. Q., Qiu, A. L., Shi, L. P., Cai, J. S., Huang, X. Y., Yang, S., et al. (2015). SRC2-1 is required in PcINF1-induced pepper immunity by acting as an interacting partner of PcINF1. J. Exp. Bot. 66, 3683–3698. doi: 10.1093/jxb/erv161
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2[-Delta Delta C(T)] Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lu, K., Guo, W., Lu, J., Yu, H., Qu, C., Tang, Z., et al. (2015). Genome-wide survey and expression profile analysis of the mitogen-activated protein kinase (MAPK) gene family in Brassica rapa. PLoS ONE 10:e0132051. doi: 10.1371/journal.pone.0132051
Mehlmer, N., Wurzinger, B., Stael, S., Hofmann-Rodrigues, D., Csaszar, E., Pfister, B., et al. (2010). The Ca(2+) -dependent protein kinase CPK3 is required for MAPK-independent salt-stress acclimation in Arabidopsis. Plant J. 63, 484–498. doi: 10.1111/j.1365-313X.2010.04257.x
Miao, Y., Laun, T. M., Smykowski, A., and Zentgraf, U. (2007). Arabidopsis MEKK1 can take a short cut: it can directly interact with senescence-related WRKY53 transcription factor on the protein level and can bind to its promoter. Plant Mol. Biol. 65, 63–76. doi: 10.1007/s11103-007-9198-z
Mosser, D. D., Theodorakis, N. G., and Morimoto, R. I. (1988). Coordinate changes in heat shock element-binding activity and HSP70 gene transcription rates in human cells. Mol. Cell. Biol. 8, 4736–4744.
Narsai, R., Castleden, I., and Whelan, J. (2010). Common and distinct organ and stress responsive transcriptomic patterns in Oryza sativa and Arabidopsis thaliana. BMC Plant Biol. 10:262. doi: 10.1186/1471-2229-10-262
Neupane, A., Nepal, M. P., Benson, B. V., Macarthur, K. J., and Piya, S. (2013a). Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus, Medicago, and Phaseolus. Plant Signal. Behav. 8:e27189. doi: 10.4161/psb.27189
Neupane, A., Nepal, M. P., Piya, S., Subramanian, S., Rohila, J. S., Reese, R. N., et al. (2013b). Identification, nomenclature, and evolutionary relationships of mitogen-activated protein kinase (MAPK) genes in soybean. Evol. Bioinform. Online 9, 363–386. doi: 10.4137/EBO.S12526
Ovečka, M., Takáč, T., Komis, G., Vadovič, P., Bekešová, S., Doskočilová, A., et al. (2014). Salt-induced subcellular kinase relocation and seedling susceptibility caused by overexpression of Medicago SIMKK in Arabidopsis. J. Exp. Bot. 65, 2335–2350. doi: 10.1093/jxb/eru115
Pan, J., Zhang, M., Kong, X., Xing, X., Liu, Y., Zhou, Y., et al. (2012). ZmMPK17, a novel maize group D MAP kinase gene, is involved in multiple stress responses. Planta 235, 661–676. doi: 10.1007/s00425-011-1510-0
Popescu, S. C., Popescu, G. V., Bachan, S., Zhang, Z., Gerstein, M., Snyder, M., et al. (2009). MAPK target networks in Arabidopsis thaliana revealed using functional protein microarrays. Genes Dev. 23, 80–92. doi: 10.1101/gad.1740009
Qin, C., Yu, C., Shen, Y., Fang, X., Chen, L., Min, J., et al. (2014). Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization. Proc. Natl. Acad. Sci. U.S.A. 111, 5135–5140. doi: 10.1073/pnas.1400975111
Ramonell, K., Berrocal-Lobo, M., Koh, S., Wan, J., Edwards, H., Stacey, G., et al. (2005). Loss-of-function mutations in chitin responsive genes show increased susceptibility to the powdery mildew pathogen Erysiphe cichoracearum. Plant Physiol. 138, 1027–1036. doi: 10.1104/pp.105.060947
Rao, K. P., Richa, T., Kumar, K., Raghuram, B., and Sinha, A. K. (2010). In silico analysis reveals 75 members of mitogen-activated protein kinase kinase kinase gene family in rice. DNA Res. 17, 139–153. doi: 10.1093/dnares/dsq011
Šamajová, O., Plíhal, O., Al-Yousif, M., Hirt, H., and Šamaj, J. (2013). Improvement of stress tolerance in plants by genetic manipulation of mitogen-activated protein kinases. Biotechnol. Adv. 31, 118–128. doi: 10.1016/j.biotechadv.2011.12.002
Schenke, D., Böttcher, C., and Scheel, D. (2011). Crosstalk between abiotic ultraviolet-B stress and biotic (flg22) stress signalling in Arabidopsis prevents flavonol accumulation in favor of pathogen defence compound production. Plant Cell Environ. 34, 1849–1864. doi: 10.1111/j.1365-3040.2011.02381.x
Seyfferth, C., and Tsuda, K. (2014). Salicylic acid signal transduction: the initiation of biosynthesis, perception and transcriptional reprogramming. Front. Plant Sci. 5:697. doi: 10.3389/fpls.2014.00697
Sheikh, A. H., Raghuram, B., Jalmi, S. K., Wankhede, D. P., Singh, P., and Sinha, A. K. (2013). Interaction between two rice mitogen activated protein kinases and its possible role in plant defense. BMC Plant Biol. 13:121. doi: 10.1186/1471-2229-13-121
Singh, R., and Jwa, N. S. (2013). The rice MAPKK-MAPK interactome: the biological significance of MAPK components in hormone signal transduction. Plant Cell Rep. 32, 923–931. doi: 10.1007/s00299-013-1437-y
Singh, R., Lee, M. O., Lee, J. E., Choi, J., Park, J. H., Kim, E. H., et al. (2012). Rice mitogen-activated protein kinase interactome analysis using the yeast two-hybrid system. Plant Physiol. 160, 477–487. doi: 10.1104/pp.112.200071
Smékalová, V., Doskočilová, A., Komis, G., and Samaj, J. (2014). Crosstalk between secondary messengers, hormones and MAPK modules during abiotic stress signalling in plants. Biotechnol. Adv. 32, 2–11. doi: 10.1016/j.biotechadv.2013.07.009
Stanko, V., Giuliani, C., Retzer, K., Djamei, A., Wahl, V., Wurzinger, B., et al. (2014). Timing is everything: highly specific and transient expression of a MAP kinase determines auxin-induced leaf venation patterns in Arabidopsis. Mol. Plant 7, 1637–1652. doi: 10.1093/mp/ssu080
Sun, W., Chen, H., Wang, J., Sun, H. W., Yang, S. K., Sang, Y. L., et al. (2015). Expression analysis of genes encoding mitogen-activated protein kinases in maize provides a key link between abiotic stress signaling and plant reproduction. Funct. Integr. Genomics 15, 107–120. doi: 10.1007/s10142-014-0410-3
van Zanten, M., Ritsema, T., Polko, J. K., Leon-Reyes, A., Voesenek, L. A., Millenaar, F. F., et al. (2012). Modulation of ethylene- and heat-controlled hyponastic leaf movement in Arabidopsis thaliana by the plant defence hormones jasmonate and salicylate. Planta 235, 677–685. doi: 10.1007/s00425-011-1528-3
Wang, G., Lovato, A., Polverari, A., Wang, M., Liang, Y. H., Ma, Y. C., et al. (2014a). Genome-wide identification and analysis of mitogen activated protein kinase kinase kinase gene family in grapevine (Vitis vinifera). BMC Plant Biol. 14:219. doi: 10.1186/s12870-014-0219-1
Wang, J., Pan, C., Wang, Y., Ye, L., Wu, J., Chen, L., et al. (2015). Genome-wide identification of MAPK, MAPKK, and MAPKKK gene families and transcriptional profiling analysis during development and stress response in cucumber. BMC Genomics 16, 386. doi: 10.1186/s12864-015-1621-2
Wang, L., Liu, Y., Cai, G., Jiang, S., Pan, J., and Li, D. (2014b). Ectopic expression of ZmSIMK1 leads to improved drought tolerance and activation of systematic acquired resistance in transgenic tobacco. J. Biotechnol. 172, 18–29. doi: 10.1016/j.jbiotec.2013.11.006
Wei, C., Liu, X., Long, D., Guo, Q., Fang, Y., Bian, C., et al. (2014). Molecular cloning and expression analysis of mulberry MAPK gene family. Plant Physiol. Biochem. 77, 108–116. doi: 10.1016/j.plaphy.2014.02.002
Wu, J., Wang, J., Pan, C., Guan, X., Wang, Y., Liu, S., et al. (2014). Genome-wide identification of MAPKK and MAPKKK gene families in tomato and transcriptional profiling analysis during development and stress response. PLoS ONE 9:e103032. doi: 10.1371/journal.pone.0103032
Wu, L., Zu, X., Zhang, H., Wu, L., Xi, Z., and Chen, Y. (2015). Overexpression of ZmMAPK1 enhances drought and heat stress in transgenic Arabidopsis thaliana. Plant Mol. Biol. 88, 429–443. doi: 10.1007/s11103-015-0333-y
Yin, Z., Wang, J., Wang, D., Fan, W., Wang, S., and Ye, W. (2013). The MAPKKK gene family in Gossypium raimondii: genome-wide identification, classification and expression analysis. Int. J. Mol. Sci. 14, 18740–18757. doi: 10.3390/ijms140918740
Yoo, S. J., Kim, S. H., Kim, M. J., Ryu, C. M., Kim, Y. C., Cho, B. H., et al. (2014). Involvement of the OsMKK4-OsMPK1 cascade and its downstream transcription factor OsWRKY53 in the wounding response in rice. Plant Pathol. J. 30, 168–177. doi: 10.5423/PPJ.OA.10.2013.0106
Zhang, D., Jiang, S., Pan, J., Kong, X., Zhou, Y., Liu, Y., et al. (2014). The overexpression of a maize mitogen-activated protein kinase gene (ZmMPK5) confers salt stress tolerance and induces defence responses in tobacco. Plant Biol. (Stuttg). 16, 558–570. doi: 10.1111/plb.12084
Zhang, L., Xi, D., Li, S., Gao, Z., Zhao, S., Shi, J., et al. (2011a). A cotton group C MAP kinase gene, GhMPK2, positively regulates salt and drought tolerance in tobacco. Plant Mol. Biol. 77, 17–31. doi: 10.1007/s11103-011-9788-7
Zhang, L., Xi, D., Luo, L., Meng, F., Li, Y., Wu, C. A., et al. (2011b). Cotton GhMPK2 is involved in multiple signaling pathways and mediates defense responses to pathogen infection and oxidative stress. FEBS J. 278, 1367–1378. doi: 10.1111/j.1742-4658.2011.08056.x
Zhang, S., Xu, R., Luo, X., Jiang, Z., and Shu, H. (2013). Genome-wide identification and expression analysis of MAPK and MAPKK gene family in Malus domestica. Gene 531, 377–387. doi: 10.1016/j.gene.2013.07.107
Zhang, Z., Wu, Y., Gao, M., Zhang, J., Kong, Q., Liu, Y., et al. (2012). Disruption of PAMP-induced MAP kinase cascade by a Pseudomonas syringae effector activates plant immunity mediated by the NB-LRR protein SUMM2. Cell Host Microbe 11, 253–263. doi: 10.1016/j.chom.2012.01.015
Keywords: pepper, MAPK, MAPKK, genome-wide analysis, gene expression
Citation: Liu Z, Shi L, Liu Y, Tang Q, Shen L, Yang S, Cai J, Yu H, Wang R, Wen J, Lin Y, Hu J, Liu C, Zhang Y, Mou S and He S (2015) Genome-wide identification and transcriptional expression analysis of mitogen-activated protein kinase and mitogen-activated protein kinase kinase genes in Capsicum annuum. Front. Plant Sci. 6:780. doi: 10.3389/fpls.2015.00780
Received: 06 June 2015; Accepted: 10 September 2015;
Published: 25 September 2015.
Edited by:Maria Carlota Vaz Patto, Instituto de Tecnologia Quimica e Biologica/Universidade Nova de Lisboa, Portugal
Reviewed by:Philipp Janitza, Leibniz Institute of Plant Biochemistry, Germany
Arsheed Sheikh, University of Warwick, UK
Copyright © 2015 Liu, Shi, Liu, Tang, Shen, Yang, Cai, Yu, Wang, Wen, Lin, Hu, Liu, Zhang, Mou and He. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shuilin He, College of Crop Science, Fujian Agriculture and Forestry University, No. 15 Shangxiadian Road, Cangshan District, Fuzhou, China, email@example.com
†These authors have contributed equally to this work.