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

The phosphorylation of a WD40-repeat protein negatively regulates flavonoid biosynthesis in Camellia sinensis under drought stress

Flavonoids constitute the main nutraceuticals in the leaves of tea plants (Camellia sinensis). To date, although it is known that drought stress can negatively impact the biosynthesis of flavonoids in tea leaves, the mechanism behind this phenomenon is unclear. Herein, we report a protein phosphorylation mechanism that negatively regulates the biosynthesis of flavonoids in tea leaves in drought conditions. Transcriptional analysis revealed the downregulation of gene expression of flavonoid biosynthesis and the upregulation of CsMPK4a encoding a mitogen-activated protein kinase in leaves. Luciferase complementation and yeast two-hybrid assays disclosed that CsMPK4a interacted with CsWD40. Phosphorylation assay in vitro, specific protein immunity, and analysis of protein mass spectrometry indicated that Ser-216, Thr-221, and Ser-253 of CsWD40 were potential phosphorylation sites of CsMPK4a. Besides, the protein immunity analysis uncovered an increased phosphorylation level of CsWD40 in tea leaves under drought conditions. Mutation of the three phosphorylation sites generated dephosphorylated CsWD403A and phosphorylated CsWD403D variants, which were introduced into the Arabidopsis ttg1 mutant. Metabolic analysis showed that the anthocyanin and proanthocyanidin content was lower in ttg1:CsWD403D transgenic plants than ttg1::CsWD403A transgenic and wild type plants. The transient overexpression of CsWD403D downregulated the anthocyanidin biosynthesis in tea leaves. The dual-fluorescein protein complementation experiment showed that CsWD403D did not interact with CsMYB5a and CsAN2, two key transcription factors of procyanidins and anthocyanidins biosynthesis in tea plant. These findings indicate that the phosphorylation of CsWD40 by CsMPK4a downregulates the flavonoid biosynthesis in tea plants in drought stresses.

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

BACKGROUND

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

METHODS AND RESULTS

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

CONCLUSION

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

Genome-wide identification of B-box zinc finger (BBX) gene family in Medicago sativa and their roles in abiotic stress responses

BACKGROUND

B-box (BBX) family is a class of zinc finger transcription factors (TFs) that play essential roles in regulating plant growth, development, as well as abiotic stress. However, no systematic analysis of BBX genes has yet been conducted in alfalfa (Medica go sativa L.), and their functions have not been elucidated up to now.

RESULTS

In this study, 28 MsBBX genes were identified from the alfalfa genome, which were clustered into 4 subfamilies according to an evolutionary tree of BBX proteins. Exon-intron structure and conserved motif analysis reflected the evolutionary conservation of MsBBXs in alfalfa. Collinearity analysis showed that segmental duplication promoted the expansion of the MsBBX family. Analysis of cis-regulatory elements suggested that the MsBBX genes possessed many growth/development-, light-, phytohormone-, and abiotic stress-related elements. MsBBX genes were differentially expressed in leaves, flowers, pre-elongated stems, elongated stems, roots and nodules, and most MsBBXs were remarkably induced by drought, salt and various plant growth regulators (ABA, JA, and SA). Further functional verification demonstrated that overexpressing of the MsBBX11 gene clearly promoted salt tolerance in transgenic Arabidopsis by regulating growth and physiological processes of seedlings.

CONCLUSIONS

This research provides insights into further functional research and regulatory mechanisms of MsBBX family genes under abiotic stress of alfalfa.

Multilayer omics landscape analyses reveal the regulatory responses of tea plants to drought stress

Adverse environments, especially drought conditions, deeply influence plant development and growth in all aspects, and the yield and quality of tea plants are largely dependent on favorable growth conditions. Although tea plant responses to drought stress (DS) have been studied, a comprehensive multilayer epigenetic, transcriptomic, and proteomic investigation of how tea responds to DS is lacking. In this study, we generated DNA methylome, transcriptome, proteome, and phosphoproteome data to explore multiple regulatory landscapes in the tea plant response to DS. An integrated multiomics analysis revealed the response of tea plants to DS at multiple regulatory levels. Furthermore, a set of DS-responsive genes involved in photosynthesis, transmembrane transportation, phytohormone metabolism and signaling, secondary metabolite pathways, transcription factors, protein kinases, posttranslational and epigenetic modification, and other key stress-responsive genes were identified for further functional investigation. These results reveal the multilayer regulatory landscape of the tea plant response to DS and provide insight into the mechanisms of these DS responses.

Genome-wide analysis of the P450 gene family in tea plant (Camellia sinensis) reveals functional diversity in abiotic stress

BACKGROUND

Cytochrome P450 (Cytochrome P450s) genes are involved in the catalysis of various reactions, including growth, development, and secondary metabolite biosynthetic pathways. However, little is known about the characteristics and functions of the P450 gene family in Camellia sinensis (C. sinensis).

RESULTS

To reveal the mechanisms of tea plant P450s coping with abiotic stresses, analyses of the tea plant P450 gene family were conducted using bioinformatics-based methods. In total, 273 putative P450 genes were identified from the genome database of C. sinensis. The results showed that P450s were well-balanced across the chromosomes I to XV of entire genome, with amino acid lengths of 268-612 aa, molecular weights of 30.95-68.5 kDa, and isoelectric points of 4.93-10.17. Phylogenetic analysis divided CsP450s into 34 subfamilies, of which CYP71 was the most abundant. The predicted subcellular localization results showed that P450 was distributed in a variety of organelles, with chloroplasts, plasma membrane,,and cytoplasm localized more frequently. The promoter region of CsP450s contained various cis-acting elements related to phytohormones and stress responses. In addition, ten conserved motifs (Motif1-Motif10) were identified in the CsP450 family proteins, with 27 genes lacking introns and only one exon. The results of genome large segment duplication showed that there were 37 pairs of genes with tandem duplication. Interaction network analysis showed that CsP450 could interact with multiple types of target genes, and there are protein interactions within the family. Tissue expression analysis showed that P450 was highly expressed in roots and stems. Moreover, qPCR analysis of the relative expression level of the gene under drought and cold stress correlated with the sequencing results.

CONCLUSIONS

This study lays the foundation for resolving the classification and functional study of P450 family genes and provides a reference for the molecular breeding of C. sinensis.

Integrating Full-Length Transcriptome and RNA Sequencing of Siberian Wildrye (Elymus sibiricus) to Reveal Molecular Mechanisms in Response to Drought Stress

Drought is one of the most significant limiting factors affecting plant growth and development on the Qinghai-Tibet Plateau (QTP). Mining the drought-tolerant genes of the endemic perennial grass of the QTP, Siberian wildrye (Elymus sibiricus), is of great significance to creating new drought-resistant varieties which can be used in the development of grassland livestock and restoring natural grassland projects in the QTP. To investigate the transcriptomic responsiveness of E. sibiricus to drought stress, PEG-induced short- and long-term drought stress was applied to two Siberian wildrye genotypes (drought-tolerant and drought-sensitive accessions), followed by third- and second-generation transcriptome sequencing analysis. A total of 40,708 isoforms were detected, of which 10,659 differentially expressed genes (DEGs) were common to both genotypes. There were 2107 and 2498 unique DEGs in the drought-tolerant and drought-sensitive genotypes, respectively. Additionally, 2798 and 1850 DEGs were identified in the drought-tolerant genotype only under short- and long-term conditions, respectively. DEGs numbering 1641 and 1330 were identified in the drought-sensitive genotype only under short- and long-term conditions, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that all the DEGs responding to drought stress in E. sibiricus were mainly associated with the mitogen-activated protein kinase (MAKP) signaling pathway, plant hormone signal transduction, the linoleic acid metabolism pathway, the ribosome pathway, and plant circadian rhythms. In addition, Nitrate transporter 1/Peptide transporter family protein 3.1 (NPF3.1) and Auxin/Indole-3-Acetic Acid (Aux/IAA) family protein 31(IAA31) also played an important role in helping E. sibiricus resist drought. This study used transcriptomics to investigate how E. sibiricus responds to drought stress, and may provide genetic resources and references for research into the molecular mechanisms of drought resistance in native perennial grasses and for breeding drought-tolerant varieties.

Genome-Wide Identification of the MAPK and MAPKK Gene Families in Response to Cold Stress in Prunus mume

Protein kinases of the MAPK cascade family (MAPKKK-MAPKK-MAPK) play an essential role in plant stress response and hormone signal transduction. However, their role in the cold hardiness of Prunus mume (Mei), a class of ornamental woody plant, remains unclear. In this study, we use bioinformatic approaches to assess and analyze two related protein kinase families, namely, MAP kinases (MPKs) and MAPK kinases (MKKs), in wild P. mume and its variety P. mume var. tortuosa. We identify 11 PmMPK and 7 PmMKK genes in the former species and 12 PmvMPK and 7 PmvMKK genes in the latter species, and we investigate whether and how these gene families contribute to cold stress responses. Members of the MPK and MKK gene families located on seven and four chromosomes of both species are free of tandem duplication. Four, three, and one segment duplication events are exhibited in PmMPK, PmvMPK, and PmMKK, respectively, suggesting that segment duplications play an essential role in the expansion and evolution of P. mume and its gene variety. Moreover, synteny analysis suggests that most MPK and MKK genes have similar origins and involved similar evolutionary processes in P. mume and its variety. A cis-acting regulatory element analysis shows that MPK and MKK genes may function in P. mume and its variety's development, modulating processes such as light response, anaerobic induction, and abscisic acid response as well as responses to a variety of stresses, such as low temperature and drought. Most PmMPKs and PmMKKs exhibited tissue-specifific expression patterns, as well as time-specific expression patterns that protect them through cold. In a low-temperature treatment experiment with the cold-tolerant cultivar P. mume 'Songchun' and the cold-sensitive cultivar 'Lve', we find that almost all PmMPK and PmMKK genes, especially PmMPK3/5/6/20 and PmMKK2/3/6, dramatically respond to cold stress as treatment duration increases. This study introduces the possibility that these family members contribute to P. mume's cold stress response. Further investigation is warranted to understand the mechanistic functions of MAPK and MAPKK proteins in P. mume development and response to cold stress.

Molecular regulation of immunity in tea plants

Tea, which is mainly produced using the young leaves and buds of tea plants (Camellia sinensis (L.) O. Kuntze), is one of the most common non-alcoholic beverages consumed in the world. The standard of tea mostly depends on the variety and quality of tea plants, which generally grow in subtropical areas, where the warm and humid conditions are also conducive to the occurrence of diseases. In fighting against pathogens, plants rely on their sophisticated innate immune systems which has been extensively studied in model plants. Many components involved in pathogen associated molecular patterns (PAMPs) triggered immunity (PTI) and effector triggered immunity (ETI) have been found. Nevertheless, the molecular regulating network against pathogens (e.g., Pseudopestalotiopsis sp., Colletotrichum sp. and Exobasidium vexans) causing widespread disease (such as grey blight disease, anthracnose, and blister blight) in tea plants is still unclear. With the recent release of the genome data of tea plants, numerous genes involved in tea plant immunity have been identified, and the molecular mechanisms behind tea plant immunity is being studied. Therefore, the recent achievements in identifying and cloning functional genes/gene families, in finding crucial components of tea immunity signaling pathways, and in understanding the role of secondary metabolites have been summarized and the opportunities and challenges in the future studies of tea immunity are highlighted in this review.

Comprehensive analysis of MAPK gene family in Populus trichocarpa and physiological characterization of PtMAPK3-1 in response to MeJA induction

Mitogen-activated protein kinases (MAPKs) play important roles in plant growth and development, as well as hormone and stress responses by signaling to eukaryotic cells, through MAPK cascade, the presence of various cues; thereby, regulating various responses. The MAPK cascade consists mainly of three gene families, MAPK, MAPKK, and MAPKKK, which activate downstream signaling pathways through sequential phosphorylation. Although the MAPK cascade gene family has been reported in several species, there is a lack of comprehensive analysis in poplar. We identified 21 MAPK genes, 11 MAPKK genes, and 104 MAPKKK genes in Populus trichocarpa. The phylogenetic classification was supported by conservative motif, gene structure and motif analysis. Whole genome duplication has an important role in the expansion of MAPK cascade genes. Analysis of promoter cis-elements and expression profiles indicates that MAPK cascade genes have important roles in plant growth and development, abiotic and biotic stresses, and phytohormone response. Expression profiling revealed a significant upregulation of PtMAPK3-1 expression in response to drought, salt and disease stresses. Poplar transiently overexpressing PtMAPK3-1 and treated with methyl jasmonic acid (MeJA) had higher catalase and peroxidase levels than non-overexpressing poplar. This work represents the first complete inventory of the MAPK cascade in P. trichocarpa, which reveals that PtMAPK3-1 is induced by the MeJA hormone and participates in the MeJA-induced enhancement of the antioxidant enzyme system.

Genome-wide identification of MAPK family genes and their response to abiotic stresses in tea plant (Camellia sinensis)

Mitogen-activated protein kinase (MAPK) cascades are conserved and universal signal transduction modules that play important roles in regulating stress responses in plants. Although MAP3K, MP2K, and MPK family in tea plant (Camellia sinensis) have been investigated, little is known about MPK family genes responding to various abiotic stresses in tea plant. In this study, we performed a comprehensive genome-wide analysis of the tea plant MAPKs (CsMPKs) family gene based on the genomic data of tea plants by bioinformatics-based methods. Here, 21 putative CsMPK genes were identified in the tea plant and divided into 4 subfamilies according to the homologous to Arabidopsis and their phylogenetic relationships. The gene structure and conserved motifs of these CsMPKs in the same group showed high similarity, suggesting that they were highly conserved and might have a similar function. The expression profiles of the CsMPK genes were further investigated by quantitative real-time reverse transcription PCR, indicating that many CsMPK genes were involved in response to cold, drought, heat, or heat combined with drought treatment, suggesting their potential roles in abiotic stress responses in tea plant. These results would provide valuable information for further exploring the functional characterization of CsMPK genes in tea plants.

Comprehensive analysis of MAPK cascade genes in sorghum (Sorghum bicolor L.) reveals SbMPK14 as a potential target for drought sensitivity regulation

The mitogen-activated protein kinase (MAPK) cascade plays a crucial role in regulating many important biological processes in plants. Here, we identified and characterized eight MAPKK and 49 MAPKKK genes in sorghum and analyzed their differential expression under drought treatment; we also characterized 16 sorghum MAPK genes. RNA-seq analysis revealed that 10 MAPK cascade genes were involved in drought stress response at the transcriptome level in sorghum. Overexpression of SbMPK14 in Arabidopsis and maize resulted in hypersensitivity to drought by promoting water loss, indicating that SbMPK14 functions as a negative regulator of the drought response. Subsequent transcriptome analysis and qRT-PCR verification of maize SbMPK14 overexpression lines revealed that SbMPK14 likely increases plant drought sensitivity by suppressing the activity of specific ERF and WRKY transcription factors. This comprehensive study provides valuable insight into the mechanistic basis of MAPK cascade gene function and their responses to drought in sorghum.

Genome-wide identification of MAPK gene family members in Fagopyrum tataricum and their expression during development and stress responses

BACKGROUND

Mitogen-activated protein kinases (MAPKs) plays essential roles in the development, hormone regulation and abiotic stress response of plants. Nevertheless, a comprehensive study on MAPK family members has thus far not been performed in Tartary buckwheat.

RESULTS

Here, we identified 16 FtMAPKs in the Fagopyrum tataricum genome. Phylogenetic analysis showed that the FtMAPK family members could be classified into Groups A, B, C and D, in which A, B and C members contain a Thr-Glu-Tyr (TEY) signature motif and Group D members contain a Thr-Asp-Tyr (TDY) signature motif. Promoter cis-acting elements showed that most ProFtMAPks contain light response elements, hormone response elements and abiotic stress response elements, and several ProFtMAPks have MYB-binding sites, which may be involved in the regulation of flavonoid biosynthesis-related enzyme gene expression. Synteny analysis indicated that FtMAPKs have a variety of biological functions. Protein interaction prediction suggested that MAPKs can interact with proteins involved in development and stress resistance. Correlation analysis further confirmed that most of the FtMAPK genes and transcription factors involved in the stress response have the same expression pattern. The transient transformation of FtMAPK1 significantly increased the antioxidant enzymes activity in Tartary buckwheat leaves. In addition, we also found that FtMAPK1 can respond to salt stress by up-regulating the transcription abundance of downstream genes.

CONCLUSIONS

A total of 16 MAPKs were identified in Tartary buckwheat, and the members of the MAPK family containing the TDY motif were found to have expanded. The same subfamily members have relatively conserved gene structures and similar protein motifs. Tissue-specific expression indicated that the expression of all FtMAPK genes varied widely in the roots, stems, leaves and flowers. Most FtMAPKs can regulate the expression of other transcription factors and participate in the abiotic stress response. Our findings comprehensively revealed the FtMAPK gene family and laid a theoretical foundation for the functional characterization of FtMAPKs.

The XTH Gene Family in Schima superba: Genome-Wide Identification, Expression Profiles, and Functional Interaction Network Analysis

Xyloglucan endotransglucosylase/hydrolase (XTH), belonging to glycoside hydrolase family 16, is one of the key enzymes in plant cell wall remodeling. Schima superba is an important timber and fireproof tree species in southern China. However, little is known about XTHs in S. superba. In the present study, a total of 34 SsuXTHs were obtained, which were classified into three subfamilies based on the phylogenetic relationship and unevenly distributed on 18 chromosomes. Furthermore, the intron-exon structure and conserved motif composition of them supported the classification and the members belonging to the same subfamily shared similar gene structures. Segmental and tandem duplication events did not lead to SsuXTH gene family expansion, and strong purifying selection pressures during evolution led to similar structure and function of SsuXTH gene family. The interaction network and cis-acting regulatory elements analysis revealed the SsuXTH expression might be regulated by multiple hormones, abiotic stresses and transcription factors. Finally, expression profiles and GO enrichment analysis showed most of the tandem repeat genes were mainly expressed in the phloem and xylem and they mainly participated in glycoside metabolic processes through the transfer and hydrolysis of xyloglucan in the cell wall and then regulated fiber elongation.

NHX Gene Family in Camellia sinensis: In-silico Genome-Wide Identification, Expression Profiles, and Regulatory Network Analysis

Salt stress affects the plant growth and productivity worldwide and NHX is one of those genes that are well known to improve salt tolerance in transgenic plants. It is well characterized in several plants, such as Arabidopsis thaliana and cotton; however, not much is known about NHXs in tea plant. In the present study, NHX genes of tea were obtained through a genome-wide search using A. thaliana as reference genome. Out of the 9 NHX genes in tea, 7 genes were localized in vacuole while the remaining 2 genes were localized in the endoplasmic reticulum (ER; CsNHX8) and plasma membrane (PM; CsNHX9), respectively. Furthermore, phylogenetic relationships along with structural analysis which includes gene structure, location, and protein-conserved motifs and domains were systematically examined and further, predictions were validated by the expression analysis. The dN/dS values show that the majority of tea NHX genes is subjected to strong purifying selection under the course of evolution. Also, functional interaction was carried out in Camellia sinensis based on the orthologous genes in A. thaliana. The expression profiles linked to various stress treatments revealed wide involvement of NHX genes from tea in response to various abiotic factors. This study provides the targets for further comprehensive identification, functional study, and also contributed for a better understanding of the NHX regulatory network in C. sinensis.

In-silico genome wide analysis of Mitogen activated protein kinase kinase kinase gene family in C. sinensis

Mitogen activated protein kinase kinase kinase (MAPKKK) form the upstream component of MAPK cascade. It is well characterized in several plants such as Arabidopsis and rice however the knowledge about MAPKKKs in tea plant is largely unknown. In the present study, MAPKKK genes of tea were obtained through a genome wide search using Arabidopsis thaliana as the reference genome. Among 59 candidate MAPKKK genes in tea, 17 genes were MEKK-like, 31 genes were Raf-like and 11 genes were ZIK- like. Additionally, phylogenetic relationships were established along with structural analysis, which includes gene structure, its location as well as conserved motifs, cis-acting regulatory elements and functional domain signatures that were systematically examined. Also, on the basis of one orthologous gene found between tea and Arabidopsis, functional interaction was carried out in C. sinensis based on an Arabidopsis association model. The expressional profiles indicated major involvement of MAPKKK genes from tea in response to various abiotic stress factors. Taken together, this study provides the targets for additional inclusive identification, functional study, and provides comprehensive knowledge for a better understanding of the MAPKKK cascade regulatory network in C. sinensis.