Comparative genomic analysis of mitogen activated protein kinase gene family in grapevine

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

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

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

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

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Expression Profiling of Mitogen-Activated Protein Kinase Genes Reveals Their Evolutionary and Functional Diversity in Different Rubber Tree (Hevea brasiliensis) Cultivars

Rubber tree (Heveabrasiliensis) is the only commercially cultivated plant for producing natural rubber, one of the most essential industrial raw materials. Knowledge of the evolutionary and functional characteristics of kinases in H. brasiliensis is limited because of the long growth period and lack of well annotated genome information. Here, we reported mitogen-activated protein kinases in H. brasiliensis (HbMPKs) by manually checking and correcting the rubber tree genome. Of the 20 identified HbMPKs, four members were validated by proteomic data. Protein motif and phylogenetic analyses classified these members into four known groups comprising Thr-Glu-Tyr (TEY) and Thr-Asp-Tyr (TDY) domains, respectively. Evolutionary and syntenic analyses suggested four duplication events: HbMPK3/HbMPK6, HbMPK8/HbMPK9/HbMPK15, HbMPK10/HbMPK12 and HbMPK11/HbMPK16/HbMPK19. Expression profiling of the identified HbMPKs in roots, stems, leaves and latex obtained from three cultivars with different latex yield ability revealed tissue- and variety-expression specificity of HbMPK paralogues. Gene expression patterns under osmotic, oxidative, salt and cold stresses, combined with cis-element distribution analyses, indicated different regulation patterns of HbMPK paralogues. Further, Ka/Ks and Tajima analyses suggested an accelerated evolutionary rate in paralogues HbMPK10/12. These results revealed HbMPKs have diverse functions in natural rubber biosynthesis, and highlighted the potential possibility of using MPKs to improve stress tolerance in future rubber tree breeding.

Plant Genome Duplication Database

Genome duplication, widespread in flowering plants, is a driving force in evolution. Genome alignments between/within genomes facilitate identification of homologous regions and individual genes to investigate evolutionary consequences of genome duplication. PGDD (the Plant Genome Duplication Database), a public web service database, provides intra- or interplant genome alignment information. At present, PGDD contains information for 47 plants whose genome sequences have been released. Here, we describe methods for identification and estimation of dates of genome duplication and speciation by functions of PGDD.The database is freely available at http://chibba.agtec.uga.edu/duplication/.

Characterization, expression patterns and functional analysis of the MAPK and MAPKK genes in watermelon (Citrullus lanatus)

BACKGROUND

Mitogen-activated protein kinase (MAPK) cascades, which consist of three functionally associated protein kinases, namely MEKKs, MKKs and MPKs, are universal signaling modules in all eukaryotes and have been shown to play critical roles in many physiological and biochemical processes in plants. However, little or nothing is known about the MPK and MKK families in watermelon.

RESULTS

In the present study, we performed a systematic characterization of the ClMPK and ClMKK families including the identification and nomenclature, chromosomal localization, phylogenetic relationships, ClMPK-ClMKK interactions, expression patterns in different tissues and in response to abiotic and biotic stress and transient expression-based functional analysis for their roles in disease resistance. Genome-wide survey identified fifteen ClMPK and six ClMKK genes in watermelon genome and phylogenetic analysis revealed that both of the ClMPK and ClMKK families can be classified into four distinct groups. Yeast two-hybrid assays demonstrated significant interactions between members of the ClMPK and ClMKK families, defining putative ClMKK2-1/ClMKK6-ClMPK4-1/ClMPK4-2/ClMPK13 and ClMKK5-ClMPK6 cascades. Most of the members in the ClMPK and ClMKK families showed differential expression patterns in different tissues and in response to abiotic (e.g. drought, salt, cold and heat treatments) and biotic (e.g. infection of Fusarium oxysporum f. sp. niveum) stresses. Transient expression of ClMPK1, ClMPK4-2 and ClMPK7 in Nicotiana benthamiana resulted in enhanced resistance to Botrytis cinerea and upregulated expression of defense genes while transient expression of ClMPK6 and ClMKK2-2 led to increased susceptibility to B. cinerea. Furthermore, transient expression of ClMPK7 also led to hypersensitive response (HR)-like cell death and significant accumulation of H2O2 in N. benthamiana.

CONCLUSION

We identified fifteen ClMPK and six ClMKK genes from watermelon and analyzed their phylogenetic relationships, expression patterns and protein-protein interactions and functions in disease resistance. Our results demonstrate that ClMPK1, ClMPK4-2 and ClMPK7 positively but ClMPK6 and ClMKK2-2 negatively regulate the resistance to B. cinerea when transiently expressed in N. benthamiana and that ClMPK7 functions as a regulator of HR-like cell death through modulating the generation of H2O2.

Genome-Wide Survey and Expression Profile Analysis of the Mitogen-Activated Protein Kinase (MAPK) Gene Family in Brassica rapa

Mitogen-activated protein kinase (MAPK) cascades are fundamental signal transduction modules in plants, controlling cell division, development, hormone signaling, and biotic and abiotic stress responses. Although MAPKs have been investigated in several plant species, a comprehensive analysis of the MAPK gene family has hitherto not been performed in Brassica rapa. In this study, we identified 32 MAPKs in the B. rapa genome by conducting BLASTP and syntenic block analyses, and screening for the essential signature motif (TDY or TEY) of plant MAPK proteins. Of the 32 BraMAPK genes retrieved from the Brassica Database, 13 exhibited exon splicing errors, excessive splicing of the 5' sequence, excessive retention of the 5' sequence, and sequencing errors of the 3' end. Phylogenetic trees of the 32 corrected MAPKs from B. rapa and of MAPKs from other plants generated by the neighbor-joining and maximum likelihood methods suggested that BraMAPKs could be divided into four groups (groups A, B, C, and D). Gene number expansion was observed for BraMAPK genes in groups A and D, which may have been caused by the tandem duplication and genome triplication of the ancestral genome of the Brassica progenitor. Except for five members of the BraMAPK10 subfamily, the identified BraMAPKs were expressed in most of the tissues examined, including callus, root, stem, leaf, flower, and silique. Quantitative real-time PCR demonstrated that at least six and five BraMAPKs were induced or repressed by various abiotic stresses and hormone treatments, respectively, suggesting their potential roles in the abiotic stress response and various hormone signal transduction pathways in B. rapa. This study provides valuable insight into the putative physiological and biochemical functions of MAPK genes in B. rapa.

Mitogen-activated protein kinase cascades in Vitis vinifera

Protein phosphorylation is one of the most important mechanisms to control cellular functions in response to external and endogenous signals. Mitogen-activated protein kinases (MAPK) are universal signaling molecules in eukaryotes that mediate the intracellular transmission of extracellular signals resulting in the induction of appropriate cellular responses. MAPK cascades are composed of four protein kinase modules: MAPKKK kinases (MAPKKKKs), MAPKK kinases (MAPKKKs), MAPK kinases (MAPKKs), and MAPKs. In plants, MAPKs are activated in response to abiotic stresses, wounding, and hormones, and during plant pathogen interactions and cell division. In this report, we performed a complete inventory of MAPK cascades genes in Vitis vinifera, the whole genome of which has been sequenced. By comparison with MAPK, MAPK kinases, MAPK kinase kinases and MAPK kinase kinase kinase kinase members of Arabidopsis thaliana, we revealed the existence of 14 MAPKs, 5 MAPKKs, 62 MAPKKKs, and 7 MAPKKKKs in Vitis vinifera. We identified orthologs of V. vinifera putative MAPKs in different species, and ESTs corresponding to members of MAPK cascades in various tissues. This work represents the first complete inventory of MAPK cascades in V. vinifera and could help elucidate the biological and physiological functions of these proteins in V. vinifera.

Identification, nomenclature, and evolutionary relationships of mitogen-activated protein kinase (MAPK) genes in soybean

Mitogen-activated protein kinase (MAPK) genes in eukaryotes regulate various developmental and physiological processes including those associated with biotic and abiotic stresses. Although MAPKs in some plant species including Arabidopsis have been identified, they are yet to be identified in soybean. Major objectives of this study were to identify GmMAPKs, assess their evolutionary relationships, and analyze their functional divergence. We identified a total of 38 MAPKs, eleven MAPKKs, and 150 MAPKKKs in soybean. Within the GmMAPK family, we also identified a new clade of six genes: four genes with TEY and two genes with TQY motifs requiring further investigation into possible legume-specific functions. The results indicated the expansion of the GmMAPK families attributable to the ancestral polyploidy events followed by chromosomal rearrangements. The GmMAPK and GmMAPKKK families were substantially larger than those in other plant species. The duplicated GmMAPK members presented complex evolutionary relationships and functional divergence when compared to their counterparts in Arabidopsis. We also highlighted existing nomenclatural issues, stressing the need for nomenclatural consistency. GmMAPK identification is vital to soybean crop improvement, and novel insights into the evolutionary relationships will enhance our understanding about plant genome evolution.

Cloning and evolutionary analysis of tobacco MAPK gene family

The mitogen-activated protein (MAP) kinase cascade is an important signaling module which is involved in biotic and abiotic stress responses as well as plant growth and development. In this study, we identified 17 tobacco MAPKs including 11 novel tobacco MAPK genes that have not been identified before. Comparative analysis with MAPK gene families from other plants, such as Athaliana thaliana, rice and poplar, suggested that tobacco MAPKs (such as NtMPK1, NtMPK3 and NtMPK8) might play similar functions in response to abiotic and biotic stresses. QRT-PCR analysis revealed that a total of 14 NtMPKs were regulated by SA and/or MeJA, suggesting their potential roles involved in plant defense response. In addition, 6 NtMPKs were induced by drought treatment, implying their roles in response to drought stress. Our results indicated that most of tobacco MAPK might be involved in plant defense response, which provides the basis for further analysis on physiological functions of tobacco MAPKs.

PGDD: a database of gene and genome duplication in plants

Genome duplication (GD) has permanently shaped the architecture and function of many higher eukaryotic genomes. The angiosperms (flowering plants) are outstanding models in which to elucidate consequences of GD for higher eukaryotes, owing to their propensity for chromosomal duplication or even triplication in a few cases. Duplicated genome structures often require both intra- and inter-genome alignments to unravel their evolutionary history, also providing the means to deduce both obvious and otherwise-cryptic orthology, paralogy and other relationships among genes. The burgeoning sets of angiosperm genome sequences provide the foundation for a host of investigations into the functional and evolutionary consequences of gene and GD. To provide genome alignments from a single resource based on uniform standards that have been validated by empirical studies, we built the Plant Genome Duplication Database (PGDD; freely available at http://chibba.agtec.uga.edu/duplication/), a web service providing synteny information in terms of colinearity between chromosomes. At present, PGDD contains data for 26 plants including bryophytes and chlorophyta, as well as angiosperms with draft genome sequences. In addition to the inclusion of new genomes as they become available, we are preparing new functions to enhance PGDD.

MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity

MCScan is an algorithm able to scan multiple genomes or subgenomes in order to identify putative homologous chromosomal regions, and align these regions using genes as anchors. The MCScanX toolkit implements an adjusted MCScan algorithm for detection of synteny and collinearity that extends the original software by incorporating 14 utility programs for visualization of results and additional downstream analyses. Applications of MCScanX to several sequenced plant genomes and gene families are shown as examples. MCScanX can be used to effectively analyze chromosome structural changes, and reveal the history of gene family expansions that might contribute to the adaptation of lineages and taxa. An integrated view of various modes of gene duplication can supplement the traditional gene tree analysis in specific families. The source code and documentation of MCScanX are freely available at http://chibba.pgml.uga.edu/mcscan2/.

Toward a comprehensive phylogenetic reconstruction of the evolutionary history of mitogen-activated protein kinases in the plant kingdom

The mitogen-activated protein kinase (MAPK) pathway is a three-tier signaling cascade that transmits cellular information from the plasma membrane to the cytoplasm where it triggers downstream responses. The MAPKs represent the last step in this cascade and are activated when both tyrosine and threonine residues in a conserved TxY motif are phosphorylated by MAPK kinases, which in turn are themselves activated by phosphorylation by MAPK kinase kinases. To understand the molecular evolution of MAPKs in the plant kingdom, we systematically conducted a Hidden-Markov-Model based screen to identify MAPKs in 13 completely sequenced plant genomes. In this analysis, we included green algae, bryophytes, lycophytes, and several mono- and eudicotyledonous species covering >800 million years of evolution. The phylogenetic relationships of the 204 identified MAPKs based on Bayesian inference facilitated the retraction of the sequence of emergence of the four major clades that are characterized by the presence of a TDY or TEY-A/TEY-B/TEY-C type kinase activation loop. We present evidence that after the split of TDY- and TEY-type MAPKs, initially the TEY-C clade emerged. This was followed by the TEY-B clade in early land plants until the TEY-A clade finally emerged in flowering plants. In addition to these well characterized clades, we identified another highly conserved clade of 45 MAPK-likes, members of which were previously described as Mak-homologous kinases. In agreement with their essential functions, molecular population genetic analysis of MAPK genes in Arabidopsis thaliana accessions reveal that purifying selection drove the evolution of the MAPK family, implying strong functional constraints on MAPK genes. Closely related MAPKs most likely subfunctionalized, a process in which differential transcriptional regulation of duplicates may be involved.