Genome-Wide Identification and Analysis of Expression Profiles of Maize Mitogen-Activated Protein Kinase Kinase Kinase

Abstract

Mitogen-activated protein kinase (MAPK) cascades are highly conserved signal transduction model in animals, yeast and plants. Plant MAPK cascades have been implicated in development and stress responses. Although MAPKKKs have been investigated in several plant species including Arabidopsis and rice, no systematic analysis has been conducted in maize. In this study, we performed a bioinformatics analysis of the entire maize genome and identified 74 MAPKKK genes. Phylogenetic analyses of MAPKKKs from maize, rice and Arabidopsis have classified them into three subgroups, which included Raf, ZIK and MEKK. Evolutionary relationships within subfamilies were also supported by exon-intron organizations and the conserved protein motifs. Further expression analysis of the MAPKKKs in microarray databases revealed that MAPKKKs were involved in important signaling pathways in maize different organs and developmental stages. Our genomics analysis of maize MAPKKK genes provides important information for evolutionary and functional characterization of this family in maize.

About the Authors:

Xiangpei Kong

Affiliation: State Key Laboratory of Crop Biology, Shandong Key Laboratory of Crop Biology, College of Life Sciences, Shandong Agricultural University, Tai’an, Shandong, China

Wei Lv

Affiliation: State Key Laboratory of Crop Biology, Shandong Key Laboratory of Crop Biology, College of Life Sciences, Shandong Agricultural University, Tai’an, Shandong, China

Dan Zhang

Affiliation: State Key Laboratory of Crop Biology, Shandong Key Laboratory of Crop Biology, College of Life Sciences, Shandong Agricultural University, Tai’an, Shandong, China

Shanshan Jiang

Affiliation: State Key Laboratory of Crop Biology, Shandong Key Laboratory of Crop Biology, College of Life Sciences, Shandong Agricultural University, Tai’an, Shandong, China

Shizhong Zhang

Affiliation: State Key Laboratory of Crop Biology, Shandong Key Laboratory of Crop Biology, College of Horticulture Science and Engineering, Shandong Agricultural University, Tai’an, Shandong, China

Dequan Li

* E-mail: [email protected]

Affiliation: State Key Laboratory of Crop Biology, Shandong Key Laboratory of Crop Biology, College of Life Sciences, Shandong Agricultural University, Tai’an, Shandong, China

Introduction

Mitogen-activated protein kinase (MAPK) cascades are conserved signal transduction pathways to translate external stimuli into cellular responses in all eukaryotes [1]. A typical MAPK cascade consists of three sequentially activated kinases. Upstream signals activate MAPK kinase kinases (MAPKKK or MEKK) which in turn phosphorylate a MAPK kinase (MKK or MEK); Subsequently, MKKs activate a specific MAPK. The downstream targets of MAPKs can be transcription factors or cytoskeletal proteins [2]–[5]. In plants, different MAPK pathways can be activated in various signaling pathways, including development, cell division, hormone responses, plant innate immunity, and abiotic stress [6]–[23].

The genome of Arabidopsis contains approximate 80 MAPKKKs, 10 MKKs, and 20 MAPKs [2]. Completion of the rice genome project has identified the existence of 75 MAPKKKs, 8 MKKs and 17 MAPKs [12], [24]. MAPKKKs act at the top of MAPK cascades and show great sequence diversity. In plants, MAPKKKs are divided into three groups according to their sequence alignment: the MEKK-like family, Raf-like family and ZIK-like family [2]. Compared with MAPKs and MKKs, plant MAPKKKs contain long N- or C- terminal regions that might function in scaffolding to recruit MKKs and MAPKs [20]. In Arabidopsis, several MAPKKKs have been involved in defense, ethylene signaling and plant development.

MEKK1-MKK4/5-MPK3/6-WRKY22/WRKY29 plays an important role in plant innate immunity [25]. MEKK1-MKK1/MKK2–MPK4 has previously been shown to play important roles in oxidative stress signaling, salt and cold stresses, whereas negatively regulates plant innate immunity [26]–[29]. More recently, Kong et al. (2012) and Zhang et al. (2012) identified that SUMM1 (MEKK2) functions as a positive regulator of the R protein SUMM2, and its activity is negatively regulated by the MEKK1-MKK1/2-MPK4 cascade [30], [31]. Recent genetic evidence has indicated that YDA–MKK4/MKK5–MPK3/MPK6 negatively regulates stomatal development through phosphorylation SPEECHLESS (SPCH) [17], [32]. CTR1 is able to inhibit MKK9–MPK3/MPK6 activation in ethylene signaling and probably acts as an unconventional MAPKKK [33], [34]. ANP2/3-MKK6-MPK4/11/13 plays roles in the regulation of cytokinesis [35]–[39]. In addition, Gao and Xiang (2008) reported that At1g73660 (AtRaf5) mutant exhibited an enhanced tolerance to salt in Arabidopsis [40]. In rice, overexpression DSM1 (OsMAPKKK6) increased the tolerance to dehydration stress due to ROS scavenging [41]. Another Raf-like MAPKKK ILA1 (OsMAPKKK43) was identified involved in mechanical tissue formation in the leaf lamina joint in rice [42].

Maize (Zea mays L.) is one of the oldest and most important world-wide crops that are relied upon for human food, animal feed and for starch ethanol production. So far, seven MAPKs and 4 MKKs have been characterized in maize [43]–[51]. However, to our knowledge, the maize MAPKKK gene family has not been characterized in detail. In this study, we performed a bioinformatics analysis of the entire maize genome and identified 74 MAPKKK genes. In addition, we provide detailed information on the genomic structures, chromosomal locations and phylogenetic tree of maize MAPKKK genes. Subsequently, we investigated their transcript profiles in different organs and developmental stages using microarray data, which will help future studies for elucidating the precise roles of MAPKKKs in maize growth and development.

Materials and Methods

Identification of MAPKKK Gene Family in Maize

The completed genome sequence of Zea mays was downloaded from the maize sequence database (http://www.maizesequence.org/index.html). For the identification of maize MAPKKK gene family, Arabidopsis and rice MAPKKK protein sequences were firstly used as query sequences to search against the maize genome database and NCBI using BLASTP program. And self BLAST of the sequences was carried out to remove the redundancy. The Pfam (http://pfam.sanger.ac.uk/search) and SMART (http://smart.embl-heidelberg.de/) databases were used to confirm each predicted maize MAPKKK protein sequence.

Gene Structure Analysis of Maize MAPKKK Genes

The information of maize MAPKKK genes, including accession number, chromosomal location, ORF length, exon-intron structure, were retrieved from the B73 maize sequencing database (http://www.maizesequence.org/index.html).

Phylogenetic Analysis of Maize MAPKKK Proteins

Multiple alignments of MAPKKK proteins were carried out using the Clustal X v1.83 program. The protein sequences of Arabidopsis and rice MAPKKK were obtained from the TIGR database and phylogenetic analysis was performed with MEGA5.0 program by neighbor-joining method and the bootstrap test was carried out with 1000 replicates.

Chromosomal Locations and Gene Duplication of MAPKKK Genes

Genes were mapped on chromosomes by identifying their chromosomal position provided in the maize sequence database. Gene duplication events of MAPKKK genes in maize B73 were also investigated. We defined the gene duplication in accordance with the criteria: 1) the alignment length covered >80% of the longer gene; 2) the aligned region had an identity >80%; 3) only one duplication event was counted for tightly linked genes. All of the relevant genes identified in the maize genomes were aligned using Clustal X v1.83 and calculated using MEGA v5.0.

Expression Analyses of the MAPKKK Genes

Microarray expression data from various datasets were obtained making use of Genevestigator (https://www.genevestigator.com/gv/) with the Maize Gene Chip platform. The maize MAPKKK expression data was obtained through searching the Maize Gene Chip using identified MAPKKK ID (Table 1).

[Figure omitted. See PDF.]

Table 1. Characteristics of MAPKKKs from maize.

https://doi.org/10.1371/journal.pone.0057714.t001

Plant Materials and Growth Conditions

For maize inbred line Qi 319 (from Shandong Academy of Agricultural Sciences), embryo of 25 days after pollination was harvested from greenhouse-grown plants in sand under 16 h of light (25°C) and 8 h of dark (20°C), and eight-week-old seedling tissues and organs were harvested for expression analysis. Samples were collected and were immediately frozen in liquid N2 for further use. Two biological replicates were performed for each sample.

RNA Isolation and Real-time Quantitative RT-PCR Expression Analysis

Total RNAs were extracted according to the instructions of Trizol reagent (Invitrogen, Carlsbad, CA, USA) from leaves of maize seedlings with different treatments. The first strand cDNAs were synthesized using First Strand cDNA Synthesis kit (Fermentas, USA).

Real-time quantification RT-PCR reactions were performed in Bio-RAD MyiQ™ Real-time PCR Detection System (Bio-Rad, USA) using the TransStart Top Green qPCR SuperMix (TransGen, China) according to the manufacturer’s instructions. Each PCR reaction (20 µl) contained 10 µl 2×real-time PCR Mix (containing SYBR Green I), 0.5 µl of each primer, and appropriately diluted cDNA. The thermal cycling conditions were 95°C for 30 s followed by 45 cycles of 95°C for 15 s, 55°C −60°C for 30 s, and 72°C for 15 s. The Zmactin gene was used as internal reference for all the qRT–PCR analysis. Each treatment was repeated three times independently. Relative gene expression was calculated according to the delta-delta Ct method of the system. The primers used are described in Table S1 in File S1.

Results and Discussion

Genome-wide Identification of MAPKKK Family in Maize

Availability of complete maize genome sequences has made it possible for the first time to identify all the MAPKKK family members in this plant species. BLAST searches of the maize sequences database and NCBI database were performed using 80 Arabidopsis and 75 rice MAPKKK sequences as query and this analysis has identified 74 putative MAPKKK gene family members in the complete maize genome, designated as ZmMAPKKK1-ZmMAPKKK74 according to their group, since there was no standard nomenclature followed for MAPKKKs neither in Arabidopsis nor in rice. All the 74 MAPKKKs had conserved protein kinase domains. Because there were alternative splice variants in some genes of the family, the following analysis was restricted to only a single variant for further analysis. The detailed information of maize MAPKKK genes identified in the present study, including accession numbers, number of amino acids, molecular weight, and isoelectric point (pI), was listed in Table 1. ZmMAPKKK ORF lengths ranged from 1062 bp (ZmMAPKKK57) to 4014 bp (ZmMAPKKK14) and the molecular weights ranged from 39.8 kDa (ZmMAPKKK57) to 148.1 kDa (ZmMAPKKK14). Since the size of maize genome (∼2300 Mb) is much larger than the genomes of Arabidopsis (125 Mb) and rice (389 Mb), MAPKKK genes in maize would be larger than that in Arabidopsis and rice. However, according to the present study, the number of maize MAPKKK genes was even smaller than that of Arabidopsis and rice (Figure 1).

[Figure omitted. See PDF.]

Figure 1. Phylogenetic tree of MAPKKKs from maize, rice and Arabidopsis.

Neighbor-joining tree was created using MEGA5.0 program with 1,000 bootstrap using full length sequences of 74 maize, 75 rice, and 80 Arabidopsis MAPKKK proteins.

https://doi.org/10.1371/journal.pone.0057714.g001

Comparative Phylogenetic Analysis of MAPKKK Gene in Maize, Arabidopsis and Rice

To examine the evolutionary relationships between different MAPKKK members in maize, Arabidopsis and rice, an unrooted tree was constructed from alignments of the full MAPKKK amino acid sequences using Neighbor-Joining (NJ) method by MEGA5.0 and phylogenetic analysis indicated that ZmMAPKKKs can be divided into three major groups: MEKK, Raf and ZIK. There were 46 MAPKKKs from maize, 43 from rice and 48 from Arabidopsis in Raf group. MEKK group contained 22 maize MAPKKKs, 22 rice MAPKKKs and 21 Arabidopsis MAPKKKs. Only 6 MAPKKKs from maize, 10 from rice and 11 from Arabidopsis were grouped into ZIK group (Figure 1).

The inspection of the phylogenetic tree indicated 19 ZmMAPKKK paralogous gene pairs and these gene pairs represented 52% of the maize MAPKKK genes family members (Figure S1 in File S1), suggesting maize MAPKKK gene family may have undergone multiple duplications during the evolution history. Phylogenetic analysis also showed that there were 16 pairs of maize/rice MAPKKK proteins in the same clade of the phylogenetic tree (Figure 1).

Gene Structural Organization and Analysis of Conserved Domain in MAPKKK Genes

Based on the predicted sequences, the maize MAPKKK gene structures were determined. As shown in Figure 2, there were 8–17 exons in most maize MEKK group genes, whereas six genes (ZmMAPKKK17, ZmMAPKKK18, ZmMAPKKK19, ZmMAPKKK20, ZmMAPKKK21 and ZmMAPKKK22) only had one exon, and one gene (ZmMAPKKK14) had 24 exons, which were consistent with the exon numbers of their orthologs in Arabidopsis and rice. All members from Raf and ZIK possessed 2–17 exons and 7–9 exons respectively. This conserved exon numbers in each subgroup among all three species supported their close evolutionary relationship and the introduced classification of subgroups.

[Figure omitted. See PDF.]

Figure 2. Exon–intron structures of maize MAPKKK genes.

Boxes, exons; green boxes, open reading frames; lines, introns. A, MEKK; B, ZIK; C, Raf.

https://doi.org/10.1371/journal.pone.0057714.g002

Using Clustal X to analyze the full protein sequences of all MAPKKKs, we found that the most of the Raf group proteins had a C-terminal kinase domain and extended N-terminal domains. However, most of the ZIK group members had N-terminal kinase domain whereas kinase domain of MEKK family protein were located either at N- or C-terminal or central part of the protein, which were consistent with their orthologs in Arabidopsis and rice (data not shown) [24]. In addition, we also investigated the conserved motif in their kinase domains. Among the three families MEKK family is relatively well characterized. Most MEKK-like proteins seem to participate in canonical MAP kinase cascades that activate downstream MKKs. AtMEKK1 and AtMEKK2 were shown to play important roles in plant innate immunity [28], [30], [52]. More recently, Hashimoto et al. (2012) reported that NbMAPKKKα, NbMAPKKKβ and NbMAPKKKγ functioned as positive regulators of PCD [53]. All the members of maize MEKK family shared conserved motif G (T/S) Px (W/F) MAPEV, which confirmed their association with MEKK family [24] (Figure 3A). ZIK-like kinases also known as WNK (With No lysine (K)), which have not been shown to phosphorylate MKKs in plants, are involved in internal rhythm. AtWNK1 phosphorylated the putative circadian clock component APRR3 in vitro and might be involved in a signal transduction cascade regulating its biological activity [54]. AtWNK2/5/8 regulated flowering time by modulating the photoperiod pathway [55]. Recently, OsWNK1 was found to respond differentially under various abiotic stresses and also showed rhythmic expression profile under diurnal and circadian conditions at the transcription level [56]. The conserved motif of ZIK family proteins in maize were investigated using Clustal X and as shown in Figure 3B, a conserved signature motif GTPEFMAPE (L/V) (Y/F) was found in all members [24]. Compared with ZIK and MEKK like families, Raf family has many more members. Two of the best-studied Arabidopsis Raf-like MAPKKKs, CTR1 and EDR1 are known to participate in ethylene-mediated signaling and defense responses. However, neither CTR1 nor EDR1 have been confirmed to participate in a classic MAPK cascade. As shown in Figure 3C, all the members of Raf family have the conserved motif GTXX (W/Y) MAPE except ZmMAPKKK47, which strongly supported their identity as members of Raf subfamily [24].

[Figure omitted. See PDF.]

Figure 3. Alignment of MAPKKK family from maize.

The highlighted part shows the conserved motif. A, MEKK; B, ZIK; C, Raf.

https://doi.org/10.1371/journal.pone.0057714.g003

Genomic Distribution and Gene Duplication

The physical locations of the MAPKKK genes on maize chromosomes were depicted in Figure 4. It was found that 73 ZmMAPKKKs were mapped on all 10 chromosomes of maize and 1 MAPKKK (ZmMAPKKK29) was situated on unanchored contigs (chromosome unknown). Ten were present on chromosomes 3 and 5; nine on chromosomes 1, 2, 4; four on chromosomes 6, 7, 10; In addition, chromosome 8 had 8 MAPKKK members, whereas chromosome 9 encoded 6 MAPKKKs members.

[Figure omitted. See PDF.]

Figure 4. Chromosomal distributions of MAPKKK genes in the maize genome.

https://doi.org/10.1371/journal.pone.0057714.g004

Gene duplication events play a significant role in the amplification of gene family members in the genome. Several rounds of genome duplication events have been found in maize genome [57]. The expansion mechanism of the maize MAPKKK gene family was analyzed to understand gene duplication events. As shown in Figure 4, nineteen paralogs of the 74 maize MAPKKKs were identified, including 17 segmental duplication events between chromosomes and the other 2 duplication events within the same chromosome (ZmMAPKKK30 and ZmMAPKKK31, ZmMAPKKK50 and ZmMAPKKK51). Furthermore, these gene pairs shared similar exon-intron structures. This result suggested the duplication events play vital roles in MAPKKK genes expansion in maize genome.

Expression Pattern of the Maize MAPKKK Genes in Different Tissues and Developmental Stages

To observe expression profiles of the MAPKKK in maize development, we analyzed the expression of the MAPKKK genes under normal growth conditions by a Genevestigator analysis (https://www.genevestigator.ethz.ch/) in 18 different tissues, including the seedlings, coleoptiles, radicles, tassel, anther, ear, silk, caryopsis, embryo, endosperm, pericarp, culm, internode, foliar leaf, juvenile leaf, adult leaf, blade and primary root. Fifty seven genes correspond to probes and there were 17 MAPKKK genes whose corresponding probes were not found. Heatmap representation of expression profile of 57 MAPKKK genes during maize development was shown in Figure 5. Eight MAPKKKs (ZmMAPKKK25, ZmMAPKKK28, ZmMAPKKK36, ZmMAPKKK43, ZmMAPKKK52, ZmMAPKKK53, ZmMAPKKK67 and ZmMAPKKK72) had higher expression in anther than that of in other organs. Eight MAPKKKs (ZmMAPKKK4, ZmMAPKKK9, ZmMAPKKK32, ZmMAPKKK33, ZmMAPKKK37, ZmMAPKKK46, ZmMAPKKK69 and ZmMAPKKK73) had higher expression in embryo than that of in endosperm, whereas ZmMAPKKK22, ZmMAPKKK29, ZmMAPKKK39 and ZmMAPKKK49 had the opposite expression profiles in embryo and endosperm. In addition, five MAPKKKs (ZmMAPKKK18, ZmMAPKKK22, ZmMAPKKK55, ZmMAPKKK63 and ZmMAPKKK62) were expressed with high abundance in primary roots which was consistent with their expression in radicle. Specifically, ZmMAPKKK10 and ZmMAPKKK11 demonstrated a unique expression pattern in silk. Furthermore, MAPKKK duplicated gene pair expression patterns were also investigated, only seven pairs (ZmMAPKKK33 and ZmMAPKKK32, ZmMAPKKK44 and ZmMAPKKK45, ZmMAPKKK52 and ZmMAPKKK53, ZmMAPKKK64 and ZmMAPKKK65, ZmMAPKKK63 and ZmMAPKKK62, ZmMAPKKK67 and ZmMAPKKK68, ZmMAPKKK70 and ZmMAPKKK71) shared the similar expression patterns in nearly all the organs, whereas other paralogs were not the case. These results showed that although the duplicated genes had higher similarities in amino acid, they may not have similar function or are involved in the same signaling pathway.

[Figure omitted. See PDF.]

Figure 5. The expression profile of 57 MAPKKK genes in maize different tissues.

The deep and light red shading represents the relative high or low expression levels, respectively.

https://doi.org/10.1371/journal.pone.0057714.g005

In addition, we also identified the expression profiles of MAPKKK family genes under different developmental stages through analysis of publicly available microarray data sets. All the 57 genes were expressed in at least one of developmental stages (Figure 6). Nine MAPKKK genes (ZmMAPKKK10, ZmMAPKKK11, ZmMAPKKK16, ZmMAPKKK25, ZmMAPKKK28, ZmMAPKKK30, ZmMAPKKK56, ZmMAPKKK70, ZmMAPKKK71) were expressed in all developmental stages mentioned in the Figure 6 except for inflorescence formation stage, whereas another nine MAPKKK genes (ZmMAPKKK1, ZmMAPKKK9, ZmMAPKKK26, ZmMAPKKK33, ZmMAPKKK45, ZmMAPKKK49, ZmMAPKKK65, ZmMAPKKK68, ZmMAPKKK73) had higher expression in inflorescence formation than that of in other developmental stages. In addition, ZmMAPKKK37, ZmMAPKKK46 and ZmMAPKKK58 had higher expression in germination stage than other genes, whereas ZmMAPKKK43, ZmMAPKKK69 and ZmMAPKKK71 had highest expression in anthesis stage. Specifically, ZmMAPKKK52 was expressed with low abundance in all stages. Moreover, several paralogs (ZmMAPKKK15 and ZmMAPKKK16, ZmMAPKKK71 and ZmMAPKKK70, ZmMAPKKK52 and ZmMAPKKK53, ZmMAPKKK62 and ZmMAPKKK63, ZmMAPKKK64 and ZmMAPKKK65) showed highly similar expression profiles, which may indicate subfunctionalization in the course of evolution. However, other gene pairs showed quit different under the maize developmental stages.

[Figure omitted. See PDF.]

Figure 6. The expression profile of 57 MAPKKK genes in maize developmental stages.

The deep and light red shading represents the relative high or low expression levels, respectively.

https://doi.org/10.1371/journal.pone.0057714.g006

Next, we used quantitative real-time RT-PCR to validate the expression patterns in different tissues resulting from microarray database. Nine genes (ZmMAPKKK10, ZmMAPKKK11, ZmMAPKKK16, ZmMAPKKK18, ZmMAPKKK27, ZmMAPKKK47, ZmMAPKKK51, ZmMAPKKK55, and ZmMAPKKK63) were selected to confirm their expression in primary root, pericarp, internode, adult leaf, silk, culm, seedling, endosperm, embryo and tassel. Surprisingly, most our qRT-PCR data did not correspond with microarray data (Figures 5 and 7). For example, our qRT-PCR results showed that ZmMAPKKK11, ZmMAPKKK18, ZmMAPKKK47 and ZmMAPKKK51 exhibited a highest expression level in embryo (Figure 7), and the microarray data showed that these four genes had higher expression in silk, root and coleoptiles than that of in embryo (Figure 5). However, ZmMAPKKK11, ZmMAPKKK18, ZmMAPKKK47 and ZmMAPKKK51 showed higher expression in dough stage (Figure 6), suggesting they may play important roles in seed development and which was consistent with our qRT-PCR data (Figure 7). The conflicting results between our qRT-PCR and microarray database may be due to the different plant materials and growth conditions, and different experimental conditions. From these results, it is speculated that most of MAPKKK genes with different expression levels in all the maize detected organs might play key roles in plant development and several MAPKKK genes may uniquely function in maize developmental stages. However, more researches are needed to determine the functions of the MAPKKK family by additional biological experiments.

[Figure omitted. See PDF.]

Figure 7. Expression patterns of the 9 MAPKKK genes in various tissues by quantitative real-time RT-PCR analysis.

Error bars indicate standard deviations (n = 3). 1, primary root; 2, pericarp; 3, internode; 4, adult leaf; 5, silk; 6, culm; 7, seedling; 8, endosperm; 9, embryo; 10, tassel.

https://doi.org/10.1371/journal.pone.0057714.g007

Conclusion

An increasing body of evidence has shown that the mitogen-activated protein kinase (MAPK) cascades are involved in plant development and stress responses. So far, MAPKKKs have been investigated in several plant species including Arabidopsis and rice, no systematic analysis has been conducted in maize. In this present study, we performed a genome-wide survey and identified 74 MAPKKK genes from maize. Phylogenetic analysis of MAPKKKs from maize, rice and Arabidopsis has classified them into three subgroups. Members within each subgroup may have recent common evolutionary origins since they shared conserved protein motifs and exon-intron structures. Furthermore, microarray analysis showed that a number of maize MAPKKK genes differentially expressed across different tissues and developmental stages. In addition, quantitative real-time RT-PCR was performed on nine selected MAPKKK genes to confirm their expression patterns in different tissues. Our observations may lay the foundation for future functional analysis of maize MAPKKK genes to unravel their biological roles.

Supporting Information

[Figure omitted. See PDF.]

File S1.

Supporting Information file contains Figure S1 and Table S1.

https://doi.org/10.1371/journal.pone.0057714.s001

(DOC)

Author Contributions

Conceived and designed the experiments: XK DL. Performed the experiments: XK DZ SJ. Analyzed the data: XK WL SZ. Contributed reagents/materials/analysis tools: XK. Wrote the paper: XK DL.

Citation: Kong X, Lv W, Zhang D, Jiang S, Zhang S, Li D (2013) Genome-Wide Identification and Analysis of Expression Profiles of Maize Mitogen-Activated Protein Kinase Kinase Kinase. PLoS ONE 8(2): e57714. https://doi.org/10.1371/journal.pone.0057714

References

1. Nakagami H, Pitzschke A, Hirt H (2005) Emerging MAP kinase pathways in plant stress signalling. Trends Plant Sci 10(7): 339–346.

2. MAPK Group (2002) Mitogen-activated protein kinase cascades in plants: a new nomenclature. Trends Plant Sci 7(7): 301–308.

3. Bethke G, Unthan T, Uhrig JF, Poschl Y, Gust AA, et al. (2009) Flg22 regulates the release of an ethylene response factor substrate from MAP kinase 6 in Arabidopsis thaliana via ethylene signaling. Proc Natl Acad Sci USA 106(19): 8067–8072.

4. Fiil BK, Petersen K, Petersen M, Mundy J (2009) Gene regulation by MAP kinase cascades. Curr Opin Plant Biol 12(5): 615–621.

5. Popescu SC, Popescu GV, Bachan S, Zhang Z, Gerstein M, et al. (2009) MAPK target networks in Arabidopsis thaliana revealed using functional protein microarrays. Gene Dev 23(1): 80–92.

6. Huang Y, Li H, Gupta R, Morris PC, Luan S, Kieber JJ (2000) ATMPK4, an Arabidopsis homolog of mitogen-activated protein kinase, is activated in vitro by AtMEK1 through threonine phosphorylation. Plant Physiol 122(4): 1301–1310.

7. Ichimura K, Mizoguchi T, Yoshida R, Yuasa T, Shinozaki K (2000) Various abiotic stresses rapidly activate Arabidopsis MAP kinases ATMPK4 and ATMPK6. Plant J 24(5): 655–665.

8. Kiegerl S, Cardinale F, Siligan C, Gross A, Baudouin E, et al. (2000) SIMKK, a mitogen-activated protein kinase (MAPK) kinase, is a specific activator of the salt stress-induced MAPK, SIMK. Plant Cell 12(11): 2247–2258.

9. Munnik T, Meijer HJ (2001) Osmotic stress activates distinct lipid and MAPK signalling pathways in plants. FEBS Lett 498(2–3): 172–178.

10. Jin H, Axtell MJ, Dahlbeck D, Ekwenna O, Zhang S, et al. (2002) NPK1, an MEKK1-like mitogen-activated protein kinase kinase kinase, regulates innate immunity and development in plants. Dev Cell 3(2): 291–297.

11. Ren D, Yang H, Zhang S (2002) Cell death mediated by MAPK is associated with hydrogen peroxide production in Arabidopsis. J Biol Chem 277(1): 559–565.

12. Hamel LP, Nicole MC, Sritubtim S, Morency MJ, et al. (2006) Ancient signals: comparative genomics of plant MAPK and MAPKK gene families. Trends Plant Sci 11(4): 192–198.

13. Xing Y, Jia W, Zhang J (2008) AtMKK1 mediates ABA-induced CAT1 expression and H2O2 production via AtMPK6-coupled signaling in Arabidopsis. Plant J 54(3): 440–451.

14. Jammes F, Song C, Shin D, Munemasa S, Takeda K, et al. (2009) MAP kinases MPK9 and MPK12 are preferentially expressed in guard cells and positively regulate ROS-mediated ABA signaling. Proc Natl Acad Sci USA 106(48): 20520–20525.

15. Zhang L, Xi D, Luo L, Meng F, Li Y, et al. (2011) Cotton GhMPK2 is involved in multiple signaling pathways and mediates defense responses to pathogen infection and oxidative stress. FEBS J 278(8): 1367–1378.

16. Kong F, Wang J, Cheng L, Liu S, Wu J, et al. (2012) Genome-wide analysis of the mitogen-activated protein kinase gene family in Solanum lycopersicum. Gene 499(1): 108–120.

17. Takahashi Y, Soyano T, Kosetsu K, Sasabe M, Machida Y (2010) HINKEL kinesin, ANP MAPKKKs and MKK6/ANQ MAPKK, which phosphorylates and activates MPK4 MAPK, constitute a pathway that is required for cytokinesis in Arabidopsis thaliana. Plant Cell Physiol 51(10): 1766–1776.

18. Colcombet J, Hirt H (2008) Arabidopsis MAPKs: a complex signalling network involved in multiple biological processes. Biochem J 413(2): 217–226.

19. Lee JS, Huh KW, Bhargava A, Ellis BE (2008) Comprehensive analysis of protein-protein interactions between Arabidopsis MAPKs and MAPK kinases helps define potential MAPK signalling modules. Plant Signal Behav 3(12): 1037–1041.

20. Rodriguez MC, Petersen M, Mundy J (2010) Mitogen-activated protein kinase signaling in plants. Annu Rev Plant Biol 61: 621–649.

21. Song D, Chen J, Song F, Zheng Z (2006) A novel rice MAPK gene, OsBIMK2, is involved in disease-resistance responses. Plant Biol 8(5): 587–596.

22. Melech-Bonfil S, Sessa G (2010) Tomato MAPKKKepsilon is a positive regulator of cell-death signaling networks associated with plant immunity. Plant J 64(3): 379–391.

23. Chen L, Hu W, Tan S, Wang M, Ma Z, et al. (2012) Genome-Wide Identification and Analysis of MAPK and MAPKK Gene Families in Brachypodium distachyon. PloS ONE 7(10): e46744.

24. Rao KP, Richa T, Kumar K, Raghuram B, Sinha AK (2010) In silico analysis reveals 75 members of mitogen-activated protein kinase kinase kinase gene family in rice. DNA Res 17(3): 139–153.

25. Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, et al. (2002) MAP kinase signalling cascade in Arabidopsis innate immunity. Nature 415(6875): 977–983.

26. Teige M, Scheikl E, Eulgem T, Doczi R, Ichimura K, et al. (2004) The MKK2 pathway mediates cold and salt stress signaling in Arabidopsis. Mol Cell 15(1): 141–152.

27. Ichimura K, Casais C, Peck SC, Shinozaki K, Shirasu K (2006) MEKK1 is required for MPK4 activation and regulates tissue-specific and temperature-dependent cell death in Arabidopsis. J Biol Chem 281(48): 36969–36976.

28. Gao M, Liu J, Bi D, Zhang Z, Cheng F, et al. (2008) MEKK1, MKK1/MKK2 and MPK4 function together in a mitogen-activated protein kinase cascade to regulate innate immunity in plants. Cell Res 18(12): 1190–1198.

29. Meszaros T, Helfer A, Hatzimasoura E, Magyar Z, Serazetdinova L, et al. (2006) The Arabidopsis MAP kinase kinase MKK1 participates in defence responses to the bacterial elicitor flagellin. Plant J 48(4): 485–498.

30. Kong Q, Qu N, Gao M, Zhang Z, Ding X, et al. (2012) The MEKK1-MKK1/MKK2-MPK4 kinase cascade negatively regulates immunity mediated by a mitogen-activated protein kinase kinase kinase in Arabidopsis. Plant Cell 24(5): 2225–2236.

31. Zhang Z, Wu Y, Gao M, Zhang J, Kong Q, 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(3): 253–263.

32. Lampard GR, Lukowitz W, Ellis BE, Bergmann DC (2009) Novel and expanded roles for MAPK signaling in Arabidopsis stomatal cell fate revealed by cell type-specific manipulations. Plant Cell 21(11): 3506–3517.

33. Xu J, Li Y, Wang Y, Liu H, Lei L, et al. (2008) Activation of MAPK kinase 9 induces ethylene and camalexin biosynthesis and enhances sensitivity to salt stress in Arabidopsis. J Biol Chem 283(40): 26996–27006.

34. Yoo SD, Cho YH, Tena G, Xiong Y, Sheen J (2008) Dual control of nuclear EIN3 by bifurcate MAPK cascades in C2H4 signalling. Nature 451(7180): 789–795.

35. Sasabe M, Soyano T, Takahashi Y, Sonobe S, Igarashi H, et al. (2006) Phosphorylation of NtMAP65–1 by a MAP kinase down-regulates its activity of microtubule bundling and stimulates progression of cytokinesis of tobacco cells. Gene Dev 20(8): 1004–1014.

36. Beck M, Komis G, Muller J, Menzel D, Samaj J (2010) Arabidopsis homologs of nucleus- and phragmoplast-localized kinase 2 and 3 and mitogen-activated protein kinase 4 are essential for microtubule organization. Plant Cell 22(3): 755–771.

37. Beck M, Komis G, Ziemann A, Menzel D, Samaj J (2011) Mitogen-activated protein kinase 4 is involved in the regulation of mitotic and cytokinetic microtubule transitions in Arabidopsis thaliana. New Phytol 189(4): 1069–1083.

38. Soyano T, Nishihama R, Morikiyo K, Ishikawa M, Machida Y (2003) NQK1/NtMEK1 is a MAPKK that acts in the NPK1 MAPKKK-mediated MAPK cascade and is required for plant cytokinesis. Gene Dev 17(8): 1055–1067.

39. Zeng Q, Chen JG, Ellis BE (2011) AtMPK4 is required for male-specific meiotic cytokinesis in Arabidopsis. Plant J 67(5): 895–906.

40. Gao L, Xiang CB (2008) The genetic locus At1g73660 encodes a putative MAPKKK and negatively regulates salt tolerance in Arabidopsis. Plant Mol Biol 67(1–2): 125–134.

41. Ning J, Li X, Hicks LM, Xiong L (2010) A Raf-like MAPKKK gene DSM1 mediates drought resistance through reactive oxygen species scavenging in rice. Plant Physiol 152(2): 876–890.

42. Ning J, Zhang B, Wang N, Zhou Y, Xiong L (2012) Increased leaf angle1, a Raf-like MAPKKK that interacts with a nuclear protein family, regulates mechanical tissue formation in the Lamina joint of rice. Plant Cell 23(12): 4334–4347.

43. Berberich T, Sano H, Kusano T (1999) Involvement of a MAP kinase, ZmMPK5, in senescence and recovery from low-temperature stress in maize. Mol Gen Genet 262(3): 534–542.

44. Lin F, Ding H, Wang J, Zhang H, Zhang A, et al. (2009) Positive feedback regulation of maize NADPH oxidase by mitogen-activated protein kinase cascade in abscisic acid signalling. J Exp Bot 60(11): 3221–3238.

45. Kong X, Pan J, Zhang M, Xing X, Zhou Y, et al. (2011) ZmMKK4, a novel group C mitogen-activated protein kinase kinase in maize (Zea mays), confers salt and cold tolerance in transgenic Arabidopsis. Plant Cell Environ 34(8): 1291–1303.

46. Kong X, Sun L, Zhou Y, Zhang M, Liu Y, et al. (2011) ZmMKK4 regulates osmotic stress through reactive oxygen species scavenging in transgenic tobacco. Plant Cell Rep 30(11): 2097–2104.

47. Wu T, Kong XP, Zong XJ, Li DP, Li DQ (2011) Expression analysis of five maize MAP kinase genes in response to various abiotic stresses and signal molecules. Mol Biol Rep 38(6): 3967–3975.

48. Liu Y, Zhou Y, Liu L, Sun L, Zhang M, et al. (2012) Maize ZmMEK1 is a single-copy gene. Mol Biol Rep 39(3): 2957–2966.

49. Pan J, Zhang M, Kong X, Xing X, Liu Y, et al. (2012) ZmMPK17, a novel maize group D MAP kinase gene, is involved in multiple stress responses. Planta 235(4): 661–676.

50. Zhang M, Pan J, Kong X, Zhou Y, Liu Y, et al. (2012) ZmMKK3, a novel maize group B mitogen-activated protein kinase kinase gene, mediates osmotic stress and ABA signal responses. J Plant Physiol 169(15): 1501–1510.

51. Zhou Y, Zhang D, Pan J, Kong X, Liu Y, et al. (2012) Overexpression of a multiple stress-responsive gene, ZmMPK4, enhances tolerance to low temperature in transgenic tobacco. Plant Physiol Biochem 58: 174–181.

52. Pitzschke A, Djamei A, Bitton F, Hirt H (2009) A major role of the MEKK1-MKK1/2-MPK4 pathway in ROS signalling. Mol Plant 2(1): 120–137.

53. Hashimoto M, Komatsu K, Maejima K, Okano Y, Shiraishi T, et al. (2012) Identification of three MAPKKKs forming a linear signaling pathway leading to programmed cell death in. BMC Plant Biol 12(1): 103.

54. Murakami-Kojima M, Nakamichi N, Yamashino T, Mizuno T (2002) The APRR3 component of the clock-associated APRR1/TOC1 quintet is phosphorylated by a novel protein kinase belonging to the WNK family, the gene for which is also transcribed rhythmically in Arabidopsis thaliana. Plant Cell Physiol 43(6): 675–683.

55. Wang Y, Liu K, Liao H, Zhuang C, Ma H, Yan X (2008) The plant WNK gene family and regulation of flowering time in Arabidopsis. Plant Biol 10(5): 548–562.

56. Kumar K, Rao KP, Biswas DK, Sinha AK (2011) Rice WNK1 is regulated by abiotic stress and involved in internal circadian rhythm. Plant Signal Behav 6(3): 316–320.

57. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, et al. (2009) The B73 maize genome: complexity diversity and dynamics. Science 326(5956): 1112–1115.