Genome-wide identification and analysis of mitogen activated protein kinase kinase kinase gene family in grapevine (Vitis vinifera)
© Wang et al.; licensee BioMed Central Ltd. 2014
Received: 29 April 2014
Accepted: 4 August 2014
Published: 27 August 2014
Mitogen-activated protein kinase kinase kinases (MAPKKKs; MAP3Ks) are important components of MAPK cascades, which are highly conserved signal transduction pathways in animals, yeast and plants, play important roles in plant growth and development. MAPKKKs have been investigated on their evolution and expression patterns in limited plants including Arabidopsis, rice and maize.
In this study, we performed a genome-wide survey and identified 45 MAPKKK genes in the grapevine genome. Chromosome location, phylogeny, gene structure and conserved protein motifs of MAPKKK family in grapevine have been analyzed to support the prediction of these genes. In the phylogenetic analysis, MAPKKK genes of grapevine have been classified into three subgroups as described for Arabidopsis, named MEKK, ZIK and RAF, also confirmed in grapevine by the analysis of conserved motifs and exon-intron organizations. By analyzing expression profiles of MAPKKK genes in grapevine microarray databases, we highlighted the modulation of different MAPKKKs in different organs and distinct developmental stages. Furthermore, we experimentally investigated the expression profiles of 45 grape MAPKKK genes in response to biotic (powdery mildew) and abiotic stress (drought), as well as to hormone (salicylic acid, ethylene) and hydrogen peroxide treatments, and identified several candidate MAPKKK genes that might play an important role in biotic and abiotic responses in grapevine, for further functional characterization.
This is the first comprehensive experimental survey of the grapevine MAPKKK gene family, which provides insights into their potential roles in regulating responses to biotic and abiotic stresses, and the evolutionary expansion of MAPKKKs is associated with the diverse requirement in transducing external and internal signals into intracellular actions in MAPK cascade in grapevine.
Plants are constantly confronted by various pathogenic and environmental stresses that challenge their survival. To deal with stresses, plants have evolved a variety of biochemical and physiological mechanisms. Stress-activated molecular pathways include multiple inter-linked regulatory networks such as protein kinase signaling cascades that can efficiently transduce input signals into suitable outputs . The best characterized protein-kinase-based amplification cascades rely on the mitogen activated protein kinases (MAPKs), which are conserved components of signal transduction in all eukaryotic organisms . The MAPK cascades rapidly transduce stress signals into various appropriate intracellular responses . The basic MAPK cascades are composed of three classes of protein kinases: MAPK kinase kinase (MAPKKK/MAP3K), MAP kinase (MAPKK/MKK) and MAPK (MAPK/MPK). MAPKKKs are the first component of the cascades that activate MAPKKs by phosphorylating two amino acids in the S/T-XXXXX-S/T (x represents any amino acid) motif of the MAPKK activation loop, and then MAPKKs become dual-specificity kinases that activate the downstream MAPK through double phosphorylation of the T-X-Y motif in the activation loop (T-loop) ,. The activated MAPK leads to the phosphorylation of transcription factors and other signaling components that regulate the expression of downstream target genes .
So far, different members in MAPK cascades have been identified and characterized by functional genomics approach in a variety of plant species, including Arabidopsis, tobacco, rice, alfalfa and poplar. Arabidopsis thaliana genome contains 80 MAPKKKs, 10 MAPKKs and 20 MAPKs ,, whereas the rice genome contains 75 MAPKKKs, 8 MAPKKs and 17 MAPKs ,. Compared with MAPKs and MAPKKs, MAPKKKs act at the top of MAPK cascades with much greater numbers and show more complexity and sequence diversity. According to characteristic sequence motifs, MAPKKKs are divided into three groups in higher plants: the MEKK-like subfamily, ZIK subfamily and Raf-like subfamily. Compared to ZIK subfamily and RAF-like subfamily, MEKK subfamily members have a less conserved protein structure . The RAF and ZIK subfamily proteins have a C-terminal kinase domain (KD) and a long N-terminal regulatory domain (RD) that might function in scaffolding to recruit MAPKKs and MAPKs ,.
In plants, MAPK cascades have been implicated in the signaling pathways related to various stresses, ethylene signaling, innate immunity and defense responses -. In Arabidopsis, the cascade MEKK1-MKK4/5-MPK3/6-WRKY22/WRKY29 plays an important role in plant innate immunity . Investigations in alfalfa (Medicago sativa) have indicated that OMTK1, a MAPKKK, was activated by hydrogen peroxide (H2O2) . Two well-studied MAPKKK are CTR1 (Constitutive Triple Response 1) and EDR1 (Enhanced Disease Resistance 1) of A. thaliana, both belonging to the RAF-like subfamily. The CTR1 multigene family encodes an essential negative regulator for ethylene-induced gene expression in Arabidopsis, while EDR1 was shown to be a negative regulator in salicylic acid-inducible defense responses  with edr1 mutants showing increased resistance to powdery mildew . In addition, it was reported that AtRaf5 mutant exhibited an enhanced tolerance to salt in Arabidopsis. Over-expression of Os-MAPKKK6 increased the tolerance to dehydration stress through ROS scavenging in rice . In contrast to several reports on MAPKKKs in Arabidopsis and rice, research on MAPKKKs in grapevine is still very limited.
Grapevine (Vitis vinifera L.) is one of the most economically valuable and most widely grown fruit crops in the world. Sequencing of the highly homozygous grapevine PN40024 genome  provides a great opportunity for analysis of the grapevine genome and gene family evolution. Previously we have validated 12 grapevine MAPK gene by gene isolation and expression . To further understand how the MAPK cascade operates in grapevine and their internal relationships, we surveyed the gene family of MAPKKKs, the top of the MAPK cascade in the grapevine genome. Fourty-five grapevine MAPKKK genes were identified by a detailed bioinformatics analysis, annotated and named according to their sequence similarity with Arabidopsis genes, as established by the grapevine scientific community ( personal communication) and their chromosomal position and gene structure were determined. In addition, we analyzed their transcript profiles in different organs and developmental stages using published microarray data. Finally, we examined their expression patterns in response to different stresses using quantitative real time polymerase chain reaction (qRT-PCR). These results indicate that the evolutionary expansion of MAPKKKs is associated with the diverse requirement in transducing external and internal signals into intracellular actions in MAPKKK-MAPKK-MAPK cascade in grapevine.
Results and discussion
Identification of MAPKKK family in grapevine and construction of a phylogenetic tree
Characteristics of MAPKKKs of grapevine
Gene length (bp)
Aminoacid length (AA)
When one or more genes in grapevine matched more than one gene in Arabidopsis, a new name was attributed consisting of the common MAPKKK term and an increasing numbering. The detailed information on the VviMAPKKK genes identified in the present study is listed in Table 1 and Additional file 1, including nomenclature, accession numbers, chromosomal localizations, gene length, number of amino acid in the protein, isoelectric point (PI) and molecular weight (MW). These genes were distributed over almost all chromosomes, except chromosome 9. The gene length ranged from 1,678 bp (VviMAPKKK34) to 52,549 bp (VviMAPKKK6). The open reading frames (ORFs) encoded polypeptides ranging from 211 AA (VviMAPKKK46) to 1425 AA (VviMAPKKK51). The predicted molecular masses ranged from 23.84 to 155.43 kD and isoelectric point value ranged 4.76-10.32 (Table 1). According to the present study, the number of grapevine MAPKKK genes is significantly smaller than those of Arabidopsis MAPKKKs (80)  and rice MAPKKKs (75) .
Phylogenetic analysis of VviMAPKKK genes
To investigate the evolutionary relationships between MAPKKK members in grapevine and Arabidopsis, and also to assign a name to grapevine MAPKKK genes (see below) a phylogenetic tree was constructed from alignments of the full coding sequences of all 125 MAPKKK genes (45 from grapevine and 80 from Arabidopsis, Additional file 2) with the procedure and parameters described in Materials and Methods (Figure 1). Based on the phylogenetic tree, grapevine MAPKKK were classified into the same corresponding categories in Arabidopsis, which include MEKK-like, RAF and ZIK subfamilies. There were 9 VviMAPKKKs and 21 AtMAPKKKs in the MEKK subfamily, only 9 VviMAPKKKs and 11 AtZIKs in ZIK subfamily, while 27 VviMAPKKKs and 48 AtRAFs grouped in the RAF subfamily (Figure 1).
In the three clades, there were many grapevine MAPKKKs clustering together, suggesting that these homologous genes may have derived from multiple duplications after the speciation of grape during the evolution. Moreover, many grapevine MAPKKK genes have their clear orthologues in the Arabidopsis genome, which suggests that these genes might be conserved for some specific functions in the two species. Interestingly, one grapevine gene, VviMAPKKK29, stands outside the main branches and was included in the ZIK family with a bootstrap values of 71%, just above the threshold of 70% established for the phylogenetic analysis (Figure 1).
Chromosomal location of VviMAPKKK genes
Gene structural analysis of VviMAPKKK genes
As shown in Figure 3 and Table S2 (Additional file 3), the number of introns in VviMAPKKK genes was highly variable, ranging from 4 to 16 introns, whereas two genes (VviMAPKKK47 and VviMAPKKK34) had only one intron. The large variation in structures of VviMAPKKK genes suggests that the grapevine genome has changed significantly during its long evolutionary history. However, a certain degree of similarity could be observed among subgroups, supporting evolutionary relationships among members of each clade. The majority of genes in the ZIK subfamily contain 6–7 introns, genes in the MEKK subfamily mostly ranged between 8 and 10 introns, the RAF subfamily showed a variable exon number but with a majority of genes ranging between 12 and 15 introns, often with very long introns. Within this frame, paralogous gene pairs generally shared highly similar exon-intron structures (Figure 3). Collectively, the divergent gene structures between the different phylogenetic subgroups suggest that duplication events of MAPKKK genes might have occurred in ancient times and that offspring genes evolved into diverse exon/intron structures, possibly to accomplish different functions in the grapevine genome.
Analysis of conserved domains among VviMAPKKKs
The pattern of amino acid residues found in many subdomains is conserved among the family members . All VviMAPKKK genes grouped under MEKK, ZIK and RAF subfamilies were further analyzed for the presence of specific signatures. Nine VviMAPKKKs and 21 AtMAPKKKs which belong to MEKK subfamily share the conserved signature motif G (T/S) Px (W/Y/F) MAPEV, as revealed by the amino acid sequence analysis of the protein kinase domain (Additional file 4: Figure S1A). Presence of this signature in 8 out of 9 VviMAPKKK further confirmed their grouping into the MEKK subfamily, while VviMAPKKK28 showed a substitution of the methionine residue with a threonine. The ZIK subfamily consists of 9 VviMAPKKKs and 11 AtZIKs. The characteristic feature of this subfamily consists of a conserved signature GTPEFMAPE (L/V) Y across all grapevine members (Additional file 4: Figure S1B). No additional kinase domains were identified in the grapevine MEKK or ZIK subfamilies, except for VviMAPKKK27 (Additional file 5: Table S3). One exception is VviMAPKKK29, which clusters together with ZIK-encoding genes at the nucleotide level (Figure 1) with a bootstrap values just above the threshold of 70%, but the alignment of the corresponding predicted protein with other grapevine MAPKKKs at the aminoacid level revealed the presence of a slightly modified RAF domain instead of the typical ZIK domain. For this reason it was simply named VviMAPKKK29, without any reference to subfamily.
The RAF subfamily is the largest of the 3 clades of MAPKKKs. Twenty-seven and 48 MAPKKKs were grouped in the RAF subfamily in grapevine and Arabidopsis, respectively. Multiple alignments of the kinase domains revealed the presence of the RAF specific signature GTxx (W/Y) MAPE in almost all grapevine MAPKKK proteins, with only slight variations in VviMAPKKK38 and 40 (Additional file 6: Figure S2). Moreover, the majority of proteins in the RAF subfamily contained additional protein domains (Additional file 2: Table S3), the most frequent one being the EDR1 domain (7 proteins) followed by the PB1 domain (5 proteins), and other additional domains with lower frequencies. Interestingly, two “stress/fungal response” domains were detected in the sequence of VviMPKKK29, which shows a relevant divergence from other members of both clades as already mentioned.
Among the components of the kinase cascade in plants, only a few MAPKKK genes have been characterized. It was shown that AtMAPKKK1 and AtMAPKKK2 played important roles in plant innate immunity ,. MAPKKK1 in Arabidopsis was found to be responsible for oxidative stress and to be involved in negative regulation of hormone signaling . It was reported that OMTK1, a MAPKKK from M. sativa, regulates oxidative stress signaling . Recently, the Arabidopsis AtZIK4 protein WNK1 (At3g04910) was demonstrated to phosphorylate the putative circadian clock component APRR3 in vitro and might be involved in the control of circadian rhythms by regulating its biological activity, suggesting a different function from that of other MAPKKKs . Two of the best-studied RAF-like MAPKKKs in Arabidopsis, CTR1 [AtRAF1] and EDR1 [AtRAF2], act as negative regulators in ethylene-induced gene expression , and in response to powdery mildew attack , respectively. However, neither CTR1 nor EDR1 have been confirmed to participate in a classic MAPK cascade . Among those genes, only EDR1 has a clear orthologue in grapevine (VviMAPKKK60 [VviRAF2]) and can be an interesting candidate to ascertain its possible analogous functions in this species, while other characterized Arabidopsis MAP3Ks show different degrees of similarity with several grapevine genes.
Expression profiles of VviMAPKKK genes in different developmental stages and tissues
The most peculiar expression of VviMAPKKK was in pollen samples where most genes showed a up- or down-regulation, in comparison to other organs. The clustering of VviMAPKKK s according to their expression profile (Figure 4) revealed that the expression of a group of genes was much higher in young tissues and organs than in ripening or senescing ones, suggesting that these VviMAPKKK are mostly related to signal transduction during development in metabolically active tissues. The decreased transcript levels of VviMAPKKK genes in Cluster A were especially evident in post-withering stages, in which berries are left to natural dehydration for about 3 months. On the opposite, VviMAPKKK transcripts in Cluster B showed a higher level in later stages of grape development and during withering, suggesting that this set of genes may be responsive to dehydration and putatively involved in the deep transcriptomic and metabolic changes controlling biosynthesis of secondary metabolites responsible for the typical aromas of wines. This information can be important for further dissection of the signal transduction pathways operating in the transition from vegetative to reproductive stages  and in the regulation of the biosynthesis of aromatic compounds in the berry, in which different groups of MAPKKK may be involved.
It should be noted that the clustering of expression profiles does not reflect phylogenetic similarities. We only found similar expression profiles for the couples of MAPKKKs 22/23, 4/26, 56/57, and 42/47. In general, genes within the 3 clades of MEKK, ZIK and RAF or even paralogous genes may have very different expression profiles and possibly serve different functions in each organ and stage. This could have resulted from post-duplication diversifications, including subfunctionalization, neofunctionalization, or sub-neofunctionalization . These results provide a basis for further investigations on the function of VviMAPKKK genes in grapevine developmental biology.
Expression profiles of VviMAPKKK genes in response to biotic and abiotic stresses
Only a limited number of genes in MAPKKK family have been functionally characterized in Arabidopsis and even less in other species . Among those characterized genes, some were shown to be involved in the response to biotic and abiotic stresses ,-. In particular, two members of Arabidopsis RAF-like MAPKKKs with a function in plant defense were characterized: CTR1 [AtRAF1], negatively regulating ethylene responses , and EDR1 [AtRAF2], acting as a negative regulator of disease resistance and ethylene-induced senescence in Arabidopsis. Gene expression patterns usually act as indicators of gene function. In the present study, we investigated the expression patterns of all VviMAPKKK genes by semi-quantitative real-time RT-PCR in response to biotic (powdery mildew) and abiotic (drought) stress conditions, as well as in response to hormones (SA, ETH) and H2O2 treatments. Powdery mildew caused by the biotrophic ascomycete Erysiphe necator Schw. adversely affects vine growth, berry quality and grape production worldwide . Salicylic acid (SA), ethylene (ETH) and hydrogen peroxide (H2O2) play central roles in biotic stress signaling upon pathogen infection. SA and ETH are signal molecules implicated in plant defense responses to pathogens ,. H2O2 is an important ROS and a critical signaling molecule in cascades leading to plant responses to pathogens and abiotic stress factors . Treated samples were collected in all cases at 6 time points, that is 4, 8, 12, 24, 48 and 72 hours post-treatment (hpt), except for samples subjected to drought stress, which were collected after 4, 8 and 12 days (dpt).
Only a few genes diverge from these 3 main groups: VviMAPKKK47, which is almost invariably repressed, and VviMAPKKK38, which is strongly induced by E. necator, by SA and by H2O2 at the same time points of 4 and 24 h post-treatment.
Examining each stress condition separately, it can be observed that E. necator caused a strong increase of transcripts of most genes in cluster A (VviMAPKKK46,50, 32, 39, 34) and additionally of VviMAPKKK31 and 38; in particular, VviMAPKKK50 showed the highest transcript abundance, between 6 and 27-fold the control (Figures 5 and 10). A few genes (VviMAPKKK4, 54 and 51) are significantly down regulated by powdery mildew infection, especially VviMAPKKK54, while other genes are variably but slightly modulated. It can be observed however that many VviMAPKKK transcripts showed a decreased abundance at the very early time point (4hpt), a trend to a more or less pronounced increase afterwards, and a new decrease at the last collection time (72 hpt). This observation might correlate with the full establishment of infection and a possible down-regulation of defense responses, as it was reported in barley that powdery mildew can induce susceptibility in infected cells . Although elucidating the exact roles of these VviMAPKKK genes in pathogen interactions requires further functional analysis, our findings provide the first gene-family-wide survey on the expression patterns of specific grapevine MAPKKK in pathological conditions, and these highly up- and down-regulated genes can be candidate genes for future investigations.
Salicylic acid and ethylene were chosen to investigate transcriptional responses of VviMAPKKKs to hormone treatments. Regarding the response to SA, Figures 6 and 10 shows a general picture of VviMAPKKKs down regulation for most genes at most time points, especially in Cluster A, in which VviMAPKKK34 and 46 show a decreased transcript abundance of more than 20 fold at 12 hpt, but increased afterwards, especially at 48 hpt. Clear increases could be detected for VviMAPKKK60 and 64 at 12 hpt and for VviMAPKKK34 and 50 at 48 hpt. Response to ethylene was striking for some VviMAPKKK genes (examples: VviMAPKKK39, with an increase above 10-fold at 8 hpt; VviMAPKKK60, induced by 9-fold at 12 hpt; VviMAPKKK34 and 46 with fold change values between 3 and 7 (Figure 7). The same VviMAPKKK60 and 64, responsive to SA at 12hpt, were also responsive to ethylene, especially at early time points, as well as VviMAPKKK52 and the paralogous couple VviMAPKKK 22/23. VviMAPKKK60 is the grape orthologue of Arabidopsis EDR1 gene. EDR1 exerts its negative control at a point of cross talk between ethylene and salicylic acid signaling . Therefore it seems interesting that treatments with both SA and ETH may induce an increase of VviMAPKKKK60 at early time points, possibly as a regulatory mechanism to keep a balance between the two pathways. SA and ETH are involved in different signal transduction pathways and their action is often considered antagonistic, but may also cooperate in regulating defense responses . Other MAPKKK genes were mostly down regulated by ethylene along the whole time course, such as VviMAPKKK30 and 36 (Figure 7). In some cases we observed a very similar expression profile in response to these two treatments, such as for VviMAPKKK30, 36 and 61, all down regulated and for VviMAPKKK27, 40 and 52 with variable but similar profiles along the time course of the experiments.
To further investigate the roles of VviMAPKKK genes, their expression levels were measured in response to the defense signaling compound H2O2. Several transcripts reacted strongly to this treatment along the whole experiment (Figures 8 and 10), especially VviMAPKKK39, which increased up to 50 fold at 24 hpt, and VviMAPKKK34, 46 and 50, with relevant increases at several collection times. A slight increase, of about 3-fold (24 hpt) was also observed for VviMAPKKK22, which shows a high degree of sequence similarity with Arabidopsis AtMAPKKK1 [MEKK1], involved in the regulation of reactive oxygen species (ROS) homeostasis . The level of this grapevine transcript was down regulated by E. necator infection and induced by ethylene. Since plant genomes generally contain a relatively large number of MAPKKK genes, this differential regulation of MAPKKK gene expression may be a mechanism by which stress responses are fine tuned, although additional work is needed to confirm this hypothesis.
To analyze the responsiveness of VviMAPKKK to abiotic stress, drought treatments were performed. Drought is a major environmental factor limiting productivity and distribution of plants . When grapevine plants were subjected to drought stress, almost all VviMAPKKK genes displayed significantly increased expression levels especially at 8 days post-treatment (Figures 9 and 10). Among all performed treatments, drought is the one which caused the strongest effect, with many VviMAPKKK transcripts increasing more than 10 fold, and 4 transcripts (VviMAPKKK22, 23, 51, and 54) showing a more than 20-fold increased expression (Figure 9). Down regulation, although rare, was also very strong in the case VviMAPKKK46, (up to 30-fold less than control). For several unregulated genes, the increase was sustained also in the last collection time of 12 dpt. The first collection time (4 days) did not reveal a significant modulation, suggesting that plants could maintain more or less stable mRNA levels in the first days after watering suspension, except for VviMAPKKK24, 34 and 46, strongly and early down regulated. Therefore, in comparison to other stresses, drought determines a peculiar response of strong transcriptional activation on most VviMAPKKKs, especially at 8 d after drought treatments, suggesting that grapevine MAPKKKs are very likely playing roles in response to drought stress. The involvement of MAPKKK in drought resistance was seldom investigated. Among the few examples there are a RAF-Like MAPKKK gene DSM1 with a function in drought and oxidative stresses signaling in rice  and a tobacco MAPKKK (NPK1) which enhanced drought tolerance in transgenic maize . Our data provide the first insight into the possible involvement of grapevine MAPKKKs in this type of stress. However, more research is needed to determine the specific functions of the MAPKKK family by additional biological experiments, particularly to investigate why such a high number of VviMAPKKKs genes are responsive to drought.
So far, MAPKKKs have only been investigated in some plant species including Arabidopsis, rice and maize, while no systematic analysis has been conducted in grapevine. In the present study we identified 45 MAPKKK coding genes in the V. vinifera genome, which were grouped into three subfamilies - MEKK, ZIK and RAF- and named according to their sequence similarity to Arabidopsis genes. The exon/intron structure, phylogeny and conserved domains strongly supported their identity as members of each subfamily. Furthermore, by exploring a previously published microarray analysis, we provided information about the expression profiles of all VviMAPKKK genes across different tissues and developmental stages. Finally, we experimentally determined the expression profiles of all grapevine MAPKKK genes in response to biotic and abiotic stress conditions, as well as hormone and H2O2 treatments. In conclusion, our work provides an inventory of VviMAPKKK genes potentially involved in environmental stresses, an initial insight into this important gene family and a number of possibly stress-related candidates for future functional analysis. This information provides a framework to unravel the biological roles of the VviMAPKKK genes family in grape and their regulatory mechanism, particularly their apparently wide implication in drought responses.
Identification of MAPKKK gene family in grapevine
MAPKKK protein sequences of Arabidopsis were used as query to search against the Proteome databases of V. vinifera (12X, V1) (http://genomes.cribi.unipd.it/grape) , Vitis-URGI (http://urgi.versailles.inra.fr/Species/Vitis) and NCBI databases (http://www.ncbi.nlm.nih.gov/) using the BLASTP program with e-values > 1E-5. The Protein family (Pfam) database (http://pfam.sanger.ac.uk/) was used to identify their protein domains using HMMER3.0. MAPKKK gene models were only accepted if they displayed the consensus sequences of dual-specificity protein kinases. Then, the online software SMART (http://smart.embl-heidelberg.de/) was used to further confirm the predicted MAPKKK protein sequences . To increase efficiency of the search, we also investigated the results of a recently published paper  reporting a comparison of gene predictions between the 8× and the 12X genome coverages. Three additional MAPKKK genes were identified in the 8X genome, which are not present in the 12X version, but we could not find evidence that these 3 genes are actually expresses in the ESTs databases, and therefore they were not included in the family.
Phylogenetic and motif analysis of the MAPKKK gene family in grapevine
The phylogenetic tree was constructed with the web service Phylogeny.fr (http://www.phylogeny.fr/) following the rules defined by the Grapevine Super Nomenclature Committee (; Grimplet J. personal communication); bootstrap values below 70% were collapsed. The Arabidopsis nomenclature used following the MAPKKK nomenclature reported in the TAIR database (https://www.arabidopsis.org/) in which MEKK genes were named MAPKKK while RAF and ZIK genes were annotated maintaining the subfamily name. Multiple-sequence alignments of computationally predicted MAPKKK proteins belonging to each group of both Arabidopsis and grapevine (including characteristic sequence motifs) were performed using ClustalX program (version 1.83)  and GeneDoc (http://www.nrbsc.org/gfx/genedoc/). The PFAM database (http://pfam.sanger.ac.uk/) was used for identification of additional conserved motifs outside the MAPKKK domain. The protein sequences of Arabidopsis MAPKKKs were obtained from the TAIR (http://www.arabidopsis.org/) database.
Gene structure, chromosomal location and gene duplication of grapevine MAPKKK genes
The information on MAPKKK genes in the grapevine genome, including accession number, chromosomal location, open reading frame (ORF) length, molecular masses, isoelectric point value (pI) and exon-intron structure were retrieved from the grapevine database (http://genomes.cribi.unipd.it/grape/). Gene duplication events of MAPKKK genes in grapevine were investigated based on three criteria: 1) The alignment length covered >90% of the longer gene; 2) The aligned region had an identity >90%; 3) Only one duplication event was counted for tightly linked genes . GSDS (Gene Structure Display Server, http://gsds.cbi.pku.edu.cn/) was exploited to illustrate exon-intron organization of MAPKKK genes .
Microarray data analysis of VviMAPKKK genes
To understand the spatial and temporal expression patterns of MAPKKK genes during the grapevine life cycle, the expression profiles of the MAPKKK genes was analyzed based on published high-throughput microarray data . In the data sets, a total of 54 grapevine samples were included (bud, inflorescence, carpel, petal, pollen, berry, withering berry, leaf, root, seed, seedling, rachis, stem, and tendril), covering most organs at several developmental stages. The expression data were transformed in log2 values. The heat map was made with software MeV4.8 (http://www.tm4.org/mev/).
Plant materials and stress treatments
PN40024 plants (V. vinifera inbred line of Pinot noir, sequenced genotype) were kindly provided by Dr. Anne-Françoise Adam-Blondon, INRA, France, and maintained in vitro on 1/2 MS medium supplied with 0.3 mg/L Indole 3-butyric acid (IBA, Sigma, USA), under a 16/8 h photoperiod (100 μmol m−2 s−1) at 25°C in the growing chamber. Five-week-old plants were used in all treatments.
For treatments with salicylic acid (SA, Sigma, USA), ethylene (ETH) (Ethephon, Sigma, USA) and hydrogen peroxide (H2O2, Sigma, USA), plants with fully expanded six to eight leaves per tissue-culture container (240 mL) were sprayed with 5 mM SA, 5 mM ethephon (as an ethylene donor) and 10 mM H2O2. All the chemicals were purchased from Sigma and dissolved in sterile distilled water. The samples (the second to fourth leaf counted from the top) were harvested at 4, 8, 12, 24, 48 and 72 h post-treatment. Each collected sample contained independent biological replicates (three independent treated plants) and three corresponding controls. For powdery mildew infection, a local strain of Erysiphe necator Schw. was maintained on PN40024 in a greenhouse. Young leaves of similar developmental stages were inoculated with E. necator by gently pressing and tapping conidia from infected leaves on healthy ones. The second to fourth leaves were sampled at 4, 8, 12, 24, 48 and 72 h post-inoculation. Each collected sample contained independent biological replicates (three independent treated plants) and three corresponding controls. For drought treatments, in vitro plants were acclimated to pots filled with a mixture of soil and sand (1:1) in the greenhouse until they grew to a length of about 40 cm with 14 leaves. The plants were watered thoroughly first and then not watered. The sixth leaves were collected at 4, 8 and 12 days after watering interruption, immediately frozen in liquid nitrogen and stored at −70°C until analysis. Every treated sample had a corresponding regularly-watered control. For each point, three independent biological replications (three independent plants) were sampled.
Total RNA isolation and qRT-PCR expression analysis
Total RNA was extracted from the collected samples according to described previously with some modifications . The concentration and purity of RNAs were examined by measuring optical density (OD) absorption ratio at 260 and 280 nm in a One Drop™ OD-1000 spectrophotometer (Thermo Fisher Scientific, USA). The RNA integrity was checked by electrophoresis on 1.0% agarose gels stained with ethidium bromide (EB). The first-strand cDNA templates were synthesized from1 μg total RNA using PrimeScritpt RT reagent Kit (TaKaRa, Japan) following the manufacturer’s instructions.
The expressions of VviMAPKKKs were examined by qRT-PCR using a SYBR Green method on an ABI 7300 Real-time PCR System (Applied Biosystems). Primers were designed by Beacon Designer 7.0 software (Premier Biosoft International, USA), based on the 3′-untranslated region and the 3′ terminal sequences of the coding region according to the predicted mRNA sequence. The amplification product of each reaction was about 200 bp. The reaction mix (total volume of 20 μL) contained: 10 μL SYBR® Premix Ex Taq™, 0.2 μL of each primer, 1 μL of template and 8.6 μL ddH2O. The PCR conditions were: pre-denaturation at 95°C for 30 s, followed by 40 cycles at 95°C for 20 s, 60°C for 20 s, and 72°C for 43 s. Grapevine actin gene (actin-101- like, VIT_12s0178g00200) was used as the internal normalize, which was previously shown to be a suitable internal standard . The relative gene expression level was calculated according to the 2-ΔΔCt method, where ΔΔCt = (Cttarget gene - Ctactin) treatment - (Cttarget gene - Ctactin) control,. To visualize the relative expression levels data are presented as the mean fold changes between treated and control samples at each time point ± standard deviations (SDs). Mean values and standard deviations (SDs) were obtained from three biological replicates, each with three technical replicates. The gene specific primers are listed in Table S1 (Additional file 9).
Statistical analyses were performed using the software SPSS version 13.0 (Chicago, IL) and Excel. All results were indicated as means ± standard deviations (SDs) based on Duncan’s multiple range test. P < 0.05 and P < 0.01 were taken as statistically significant or highly significant, respectively.
Availability of supporting data
Here we are with the supporting data (including sequence data, microarray data and expression data) as additional files.
This research was supported partly by the Ministry of Education-Nanjing Agricultural University Project # KYZ201311.
We are grateful to Dr. Jérôme Grimplet (Instituto de Ciencias de la Vid y del Vino, CSIC, Universidad de La Rioja, Gobierno de La Rioja, Logroño, 26006, Spain) for sharing unpublished information and kind support in nomenclature assignment.
- Hamel LP, Nicole MC, Sritubtim S, Morency MJ, Ellis M, Ehlting J, Beaudoin N, Barbazuk B, Klessig D, Lee J, Martin G, Mundy J, Ohashi Y, Scheel D, Sheen J, Xing T, Zhang S, Seguin A, Ellis BE: Ancient signals: comparative genomics of plant MAPK and MAPKK gene families. Trends Plant Sci. 2006, 11 (4): 192-198. 10.1016/j.tplants.2006.02.007.View ArticlePubMedGoogle Scholar
- Doczi R, Okresz L, Romero AE, Paccanaro A, Bogre L: Exploring the evolutionary path of plant MAPK networks. Trends Plant Sci. 2012, 17 (9): 518-525. 10.1016/j.tplants.2012.05.009.View ArticlePubMedGoogle Scholar
- Ichimura K, Shinozaki K, Tena G, Sheen J, Henry Y, Champion A, Kreis M, Zhang SQ, Hirt H, Wilson C, Heberle-Bors E, Ellis BE, Morris PC, Innes RW, Ecker JR, Scheel D, Klessig DF, Machida Y, Mundy J, Ohashi Y, Walker JC, Mapk G: Mitogen-activated protein kinase cascades in plants: a new nomenclature. Trends Plant Sci. 2002, 7 (7): 301-308. 10.1016/S1360-1385(02)02302-6.View ArticleGoogle Scholar
- Rodriguez MC, Petersen M, Mundy J: Mitogen-activated protein kinase signaling in plants. Annu Rev Plant Biol. 2010, 61: 621-649. 10.1146/annurev-arplant-042809-112252.View ArticlePubMedGoogle Scholar
- Fiil BK, Petersen K, Petersen M, Mundy J: Gene regulation by MAP kinase cascades. Curr Opin Plant Biol. 2009, 12 (5): 615-621. 10.1016/j.pbi.2009.07.017.View ArticlePubMedGoogle Scholar
- Colcombet J, Hirt H: Arabidopsis MAPKs: a complex signalling network involved in multiple biological processes. Biochem J. 2008, 413 (2): 217-226. 10.1042/BJ20080625.View ArticlePubMedGoogle Scholar
- Pitzschke A, Schikora A, Hirt H: MAPK cascade signalling networks in plant defence. Curr Opin Plant Biol. 2009, 12 (4): 421-426. 10.1016/j.pbi.2009.06.008.View ArticlePubMedGoogle Scholar
- Rao KP, Richa T, Kumar K, Raghuram B, Sinha AK: In silico analysis reveals 75 members of mitogen-activated protein kinase kinase kinase gene family in rice. DNA Res. 2010, 17 (3): 139-153. 10.1093/dnares/dsq011.PubMed CentralView ArticlePubMedGoogle Scholar
- Wankhede DP, Misra M, Singh P, Sinha AK: Rice mitogen activated protein kinase kinase and mitogen activated protein kinase interaction network revealed by in-silico docking and yeast two-hybrid approaches. PLoS One. 2013, 8 (5): e65011-10.1371/journal.pone.0065011.PubMed CentralView ArticlePubMedGoogle Scholar
- Gao L, Xiang CB: The genetic locus At1g73660 encodes a putative MAPKKK and negatively regulates salt tolerance in Arabidopsis. Plant Mol Biol. 2008, 67 (1–2): 125-134. 10.1007/s11103-008-9306-8.View ArticlePubMedGoogle Scholar
- Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, Gomez-Gomez L, Boller T, Ausubel FM, Sheen J: MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002, 415 (6875): 977-983. 10.1038/415977a.View ArticlePubMedGoogle Scholar
- Zhang S, Klessig DF: MAPK cascades in plant defense signaling. Trends Plant Sci. 2001, 6 (11): 520-527. 10.1016/S1360-1385(01)02103-3.View ArticlePubMedGoogle Scholar
- Nakagami H, Kiegerl S, Hirt H: OMTK1, a novel MAPKKK, channels oxidative stress signaling through direct MAPK interaction. J Biol Chem. 2004, 279 (26): 26959-26966. 10.1074/jbc.M312662200.View ArticlePubMedGoogle Scholar
- Adams-Phillips L, Barry C, Kannan P, Leclercq J, Bouzayen M, Giovannoni J: Evidence that CTR1-mediated ethylene signal transduction in tomato is encoded by a multigene family whose members display distinct regulatory features. Plant Mol Biol. 2004, 54 (3): 387-404. 10.1023/B:PLAN.0000036371.30528.26.View ArticlePubMedGoogle Scholar
- Frye CA, Innes RW: An Arabidopsis mutant with enhanced resistance to powdery mildew. Plant Cell. 1998, 10 (6): 947-956. 10.1105/tpc.10.6.947.PubMed CentralView ArticlePubMedGoogle Scholar
- Frye CA, Tang D, Innes RW: Negative regulation of defense responses in plants by a conserved MAPKK kinase. Proc Natl Acad Sci U S A. 2001, 98 (1): 373-378. 10.1073/pnas.98.1.373.PubMed CentralView ArticlePubMedGoogle Scholar
- Ning J, Li X, Hicks LM, Xiong L: A Raf-like MAPKKK gene DSM1 mediates drought resistance through reactive oxygen species scavenging in rice. Plant Physiol. 2010, 152 (2): 876-890. 10.1104/pp.109.149856.PubMed CentralView ArticlePubMedGoogle Scholar
- Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyere C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449 (7161): U463-U465. 10.1038/nature06148.View ArticleGoogle Scholar
- Wang G, Lovato A, Liang YH, Wang M, Chen F, Tornielli GB, Polverari A, Pezzotti M, Cheng ZM: Validation by isolation and expression analyses of MAPK gene family in grapevine (Vitis vinifera). Aust J Grape Wine Res. 2014, 20 (2): 255-262. 10.1111/ajgw.12081.View ArticleGoogle Scholar
- Adam-Blondon F, Grimplet J, Adam-Blondon AF, Bert PF, Bitz O, Cantu D, Cramer G, Pezzotti M, Rombauts S: Towards the Improvement of the Gene Centered Information in Grapevine Genomics. XXII° Plant and Animal Genome Conference, San Diego, CA, USA; 2014.Google Scholar
- Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, Claverie JM, Gascuel O: Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36 (Web Server issue): W465-W469. 10.1093/nar/gkn180.PubMed CentralView ArticlePubMedGoogle Scholar
- Shiu SH, Bleecker AB: Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis. Plant Physiol. 2003, 132 (2): 530-543. 10.1104/pp.103.021964.View ArticlePubMedGoogle Scholar
- Zhang Y, Gao M, Singer SD, Fei Z, Wang H, Wang X: Genome-wide identification and analysis of the TIFY gene family in grape. PLoS One. 2012, 7 (9): e44465-10.1371/journal.pone.0044465.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao J, Huang J, Yang Y, Hu X: Analyses of the oligopeptide transporter gene family in poplar and grape. BMC Genomics. 2011, 12: 465-10.1186/1471-2164-12-465.PubMed CentralView ArticlePubMedGoogle Scholar
- Kong Q, Qu N, Gao M, Zhang Z, Ding X, Yang F, Li Y, Dong OX, Chen S, Li X, Zhang Y: The MEKK1-MKK1/MKK2-MPK4 kinase cascade negatively regulates immunity mediated by a mitogen-activated protein kinase kinase kinase in Arabidopsis. Plant Cell. 2012, 24 (5): 2225-2236. 10.1105/tpc.112.097253.PubMed CentralView ArticlePubMedGoogle Scholar
- Pitzschke A, Djamei A, Bitton F, Hirt H: A major role of the MEKK1-MKK1/2-MPK4 pathway in ROS signalling. Mol Plant. 2009, 2 (1): 120-137. 10.1093/mp/ssn079.PubMed CentralView ArticlePubMedGoogle Scholar
- Nakagami H, Soukupova H, Schikora A, Zarsky V, Hirt H: A Mitogen-activated protein kinase kinase kinase mediates reactive oxygen species homeostasis in Arabidopsis. J Biol Chem. 2006, 281 (50): 38697-38704. 10.1074/jbc.M605293200.View ArticlePubMedGoogle Scholar
- Fung RW, Gonzalo M, Fekete C, Kovacs LG, He Y, Marsh E, McIntyre LM, Schachtman DP, Qiu W: Powdery mildew induces defense-oriented reprogramming of the transcriptome in a susceptible but not in a resistant grapevine. Plant Physiol. 2008, 146 (1): 236-249. 10.1104/pp.107.108712.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin Z, Alexander L, Hackett R, Grierson D: LeCTR2, a CTR1-like protein kinase from tomato, plays a role in ethylene signalling, development and defence. Plant J. 2008, 54 (6): 1083-1093. 10.1111/j.1365-313X.2008.03481.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Kong X, Lv W, Zhang D, Jiang S, Zhang S, Li D: Genome-wide identification and analysis of expression profiles of maize mitogen-activated protein kinase kinase kinase. PLoS One. 2013, 8 (2): e57714-10.1371/journal.pone.0057714.PubMed CentralView ArticlePubMedGoogle Scholar
- Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, Porceddu A, Venturini L, Bicego M, Murino V, Ferrarini A, Delledonne M, Pezzotti M: The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012, 24 (9): 3489-3505. 10.1105/tpc.112.100230.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang X, Tuskan GA, Cheng MZ: Divergence of the Dof gene families in poplar, Arabidopsis, and rice suggests multiple modes of gene evolution after duplication. Plant Physiol. 2006, 142 (3): 820-830. 10.1104/pp.106.083642.PubMed CentralView ArticlePubMedGoogle Scholar
- Ichimura K, Casais C, Peck SC, Shinozaki K, Shirasu K: MEKK1 is required for MPK4 activation and regulates tissue-specific and temperature dependent cell death in Arabidopsis. J Biol Chem. 2006, 281 (48): 36969-36976. 10.1074/jbc.M605319200.View ArticlePubMedGoogle Scholar
- Suarez-Rodriguez MC, Adams-Phillips L, Liu Y, Wang H, Su SH, Jester PJ, Zhang S, Bent AF, Krysan PJ: MEKK1 is required for flg22-induced MPK4 activation in Arabidopsis plants. Plant Physiol. 2007, 143 (2): 661-669. 10.1104/pp.106.091389.PubMed CentralView ArticlePubMedGoogle Scholar
- Clark KL, Larsen PB, Wang X, Chang C: Association of the Arabidopsis CTR1 Raf-like kinase with the ETR1 and ERS ethylene receptors. Proc Natl Acad Sci U S A. 1998, 95 (9): 5401-5406. 10.1073/pnas.95.9.5401.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao Y, Wei T, Yin KQ, Chen Z, Gu H, Qu LJ, Qin G: Arabidopsis RAP2.2 plays an important role in plant resistance to Botrytis cinerea and ethylene responses. New Phytol. 2012, 195 (2): 450-460. 10.1111/j.1469-8137.2012.04160.x.View ArticlePubMedGoogle Scholar
- Hamel LP, Miles GP, Samuel MA, Ellis BE, Seguin A, Beaudoin N: Activation of stress-responsive mitogen-activated protein kinase pathways in hybrid poplar (Populus trichocarpa x Populus deltoides). Tree Physiol. 2005, 25 (3): 277-288. 10.1093/treephys/25.3.277.View ArticlePubMedGoogle Scholar
- Overmyer K, Brosche M, Kangasjarvi J: Reactive oxygen species and hormonal control of cell death. Trends Plant Sci. 2003, 8 (7): 335-342. 10.1016/S1360-1385(03)00135-3.View ArticlePubMedGoogle Scholar
- Yahata E, Sugai K, Penjore K, Hanboonsong Y, Takada Y, Nishiguchi M, Yamaoka N: Susceptibility-inducing factor (suppressor) from Blumeria graminis f. sp hordei has no effect on the primary infection of the fungus. Physiol Mol Plant P. 2009, 74 (2): 161-166. 10.1016/j.pmpp.2009.11.002.View ArticleGoogle Scholar
- Wang KL, Li H, Ecker JR: Ethylene biosynthesis and signaling networks. Plant Cell. 2002, 14 (Suppl): S131-S151.PubMed CentralPubMedGoogle Scholar
- Tang D, Christiansen KM, Innes RW: Regulation of plant disease resistance, stress responses, cell death, and ethylene signaling in Arabidopsis by the EDR1 protein kinase. Plant Physiol. 2005, 138 (2): 1018-1026. 10.1104/pp.105.060400.PubMed CentralView ArticlePubMedGoogle Scholar
- Shi J, Zhang L, An H, Wu C, Guo X: GhMPK16, a novel stress-responsive group D MAPK gene from cotton, is involved in disease resistance and drought sensitivity. BMC Mol Biol. 2011, 12: 22-10.1186/1471-2199-12-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Shou H, Bordallo P, Wang K: Expression of the Nicotiana protein kinase (NPK1) enhanced drought tolerance in transgenic maize. J Exp Bot. 2004, 55 (399): 1013-1019. 10.1093/jxb/erh129.View ArticlePubMedGoogle Scholar
- Grimplet J, Van Hemert J, Carbonell-Bejerano P, Diaz-Riquelme J, Dickerson J, Fennell A, Pezzotti M, Martinez-Zapater JM: Comparative analysis of grapevine whole-genome gene predictions, functional annotation, categorization and integration of the predicted gene sequences. BMC Res Notes. 2012, 5: 213-10.1186/1756-0500-5-213.PubMed CentralView ArticlePubMedGoogle Scholar
- Letunic I, Doerks T, Bork P: SMART 7: recent updates to the protein domain annotation resource. Nucleic Acids Res. 2012, 40 (Database issue): D302-D305. 10.1093/nar/gkr931.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu Z, Cavalcanti A, Chen FC, Bouman P, Li WH: Extent of gene duplication in the genomes of Drosophila, nematode, and yeast. Mol Biol Evol. 2002, 19 (3): 256-262. 10.1093/oxfordjournals.molbev.a004079.View ArticlePubMedGoogle Scholar
- Guo AY, Zhu QH, Chen X, Luo JC: GSDS: a gene structure display server. Yi Chuan. 2007, 29 (8): 1023-1026. 10.1360/yc-007-1023.View ArticlePubMedGoogle Scholar
- Ye X, Kang BG, Osburn LD, Li Y, Zong-Ming C: Identification of the flavin-dependent monooxygenase-encoding YUCCA gene family in Populus trichocarpa and their expression in vegetative tissues and in response to hormone and environmental stresses. Plant Cell Tiss Org. 2009, 97 (3): 271-283. 10.1007/s11240-009-9526-x.View ArticleGoogle Scholar
- Udvardi MK, Czechowski T, Scheible WR: Eleven golden rules of quantitative RT-PCR. Plant Cell. 2008, 20 (7): 1736-1737. 10.1105/tpc.108.061143.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.