Genome-wide identification, expression profiles and regulatory network of MAPK cascade gene family in barley

Functional characterization of geranyl/farnesyl diphosphate synthase in Wurfbainia villosa and Wurfbainia longiligularis

Wurfbainia villosa and Wurfbainia longiligularis are the two primary plant sources of Fructus Amomi, a traditional Chinese medicine. Both plants are rich in volatile terpenoids, including monoterpenes and sesquiterpenes, which are the primary medicinal components of Fructus Amomi. The trans-isopentenyl diphosphate synthase (TIDS) gene family plays a key part in determining terpenoid diversity and accumulation. However, the TIDS gene family have not been identified in W. villosa and W. longiligularis. This study identified thirteen TIDS genes in W. villosa and eleven TIDS genes in W. longiligularis, which may have expanded through segmental replication events. Based on phylogenetic analysis and expression levels, eight candidate WvTIDSs and five WlTIDSs were selected for cloning. Functional characterization in vitro demonstrated that four homologous geranyl diphosphate synthases (GPPSs) (WvGPPS1, WvGPPS2, WlGPPS1, WlGPPS2) and two geranylgeranyl diphosphate synthases (GGPPSs) (WvGGPPS and WlGGPPS) were responsible for catalyzing the biosynthesis of geranyl diphosphate (GPP), whereas two farnesyl diphosphate synthases (FPPSs) (WvFPPS and WlFPPS) catalysed the biosynthesis of the farnesyl diphosphate (FPP). A comparison of six proteins with identified GPPS functions showed that WvGGPPS and WlGGPPS exhibited the highest activity levels. These findings indicate that homologous GPPS and GGPPS together promote the biosynthesis of GPP in W. villosa and W. longiligularis, thus providing sufficient precursors for the synthesis of monoterpenes and providing key genetic elements for Fructus Amomi variety improvement and molecular breeding.

A close-up of regulatory networks and signaling pathways of MKK5 in biotic and abiotic stresses

Mitogen-activated protein Kinase Kinase 5 (MKK5) is a central hub in the complex phosphorylation chain reaction of the Mitogen-activated protein kinases (MAPK) cascade, regulating plant responses to biotic and abiotic stresses. This review manuscript aims to provide a comprehensive analysis of the regulatory mechanism of the MKK5 involved in stress adaptation. This review will delve into the intricate post-transcriptional and post-translational modifications of the MKK5, discussing how they affect its expression, activity, and subcellular localization in response to stress signals. We also discuss the integration of the MKK5 into complex signaling pathways, orchestrating plant immunity against pathogens and its modulating role in regulating abiotic stresses, such as: drought, cold, heat, and salinity, through the phytohormonal signaling pathways. Furthermore, we highlight potential applications of the MKK5 for engineering stress-resilient crops and provide future perspectives that may pave the way for future studies. This review manuscript aims to provide valuable insights into the mechanisms underlying MKK5 regulation, bridge the gap from numerous previous findings, and offer a firm base in the knowledge of MKK5, its regulating roles, and its involvement in environmental stress regulation.

Characterization and expression profiles of the ZmLBD gene family in Zea mays

BACKGROUND

The Lateral Organ Boundaries Domain (LBD) gene family is a family of plant-specific transcription factors (TFs) that are widely involved in processes such as lateral organ formation, stress response, and nutrient metabolism. However, the function of LBD genes in maize remains poorly understood.

METHODS AND RESULTS

In this study, a total of 49 ZmLBD genes were identified at the genome-wide level of maize, they were classified into nine branches based on phylogenetic relationships, and all of them were predicted to be nuclear localized. The 49 ZmLBD genes formed eight pairs of segmental duplicates, and members of the same branches' members had similar gene structure and conserved motif composition. The promoters of ZmLBD genes contain multiple types of cis-acting elements. In addition, by constructing the regulatory network of ZmLBD and other genes and miRNAs, 12 and 22 ZmLBDs were found to be involved in the gene regulatory network and miRNA regulatory network, respectively. The expression pattern analysis suggests that ZmLBD genes may be involved in different biological pathways, and drought stress induced the expressions of two inbred lines.

CONCLUSIONS

The findings enhance our comprehension of the potential roles of the ZmLBD gene family in maize growth and development, which is pivotal for genetic enhancement and breeding efforts pertaining to this significant crop.

Bioinformatic Identification and Expression Analyses of the MAPK-MAP4K Gene Family Reveal a Putative Functional MAP4K10-MAP3K7/8-MAP2K1/11-MAPK3/6 Cascade in Wheat (Triticum aestivum L.)

The mitogen-activated protein kinase (MAPK) cascades act as crucial signaling modules that regulate plant growth and development, response to biotic/abiotic stresses, and plant immunity. MAP3Ks can be activated through MAP4K phosphorylation in non-plant systems, but this has not been reported in plants to date. Here, we identified a total of 234 putative TaMAPK family members in wheat (Triticum aestivum L.). They included 48 MAPKs, 17 MAP2Ks, 144 MAP3Ks, and 25 MAP4Ks. We conducted systematic analyses of the evolution, domain conservation, interaction networks, and expression profiles of these TaMAPK-TaMAP4K (representing TaMAPK, TaMAP2K, TaMAP3K, and TaMAP4K) kinase family members. The 234 TaMAPK-TaMAP4Ks are distributed on 21 chromosomes and one unknown linkage group (Un). Notably, 25 of these TaMAP4K family members possessed the conserved motifs of MAP4K genes, including glycine-rich motif, invariant lysine (K) motif, HRD motif, DFG motif, and signature motif. TaMAPK3 and 6, and TaMAP4K10/24 were shown to be strongly expressed not only throughout the growth and development stages but also in response to drought or heat stress. The bioinformatics analyses and qRT-PCR results suggested that wheat may activate the MAP4K10-MEKK7-MAP2K11-MAPK6 pathway to increase drought resistance in wheat, and the MAP4K10-MAP3K8-MAP2K1/11-MAPK3 pathway may be involved in plant growth. In general, our work identified members of the MAPK-MAP4K cascade in wheat and profiled their potential roles during their response to abiotic stresses and plant growth based on their expression pattern. The characterized cascades might be good candidates for future crop improvement and molecular breeding.

HvMPK4 phosphorylates HvWRKY1 to enhance its suppression of barley immunity to powdery mildew fungus

Mitogen-activated protein kinase (MAPK) cascades play important roles in disease resistance in model plant species. However, the functions of MAPK signaling pathways in crop disease resistance are largely unknown. Here we report the function of HvMKK1-HvMPK4-HvWRKY1 module in barley immune system. HvMPK4 is identified to play a negative role in barley immune response against Bgh, as virus-induced gene silencing of HvMPK4 results in enhanced disease resistance whilst stably overexpressing HvMPK4 leads to super-susceptibility to Bgh infection. Furthermore, the barley MAPK kinase HvMKK1 is found to specifically interact with HvMPK4, and the activated HvMKK1DD variant specifically phosphorylates HvMPK4 in vitro. Moreover, the transcription factor HvWRKY1 is identified to be a downstream target of HvMPK4 and phosphorylated by HvMPK4 in vitro in the presence of HvMKK1DD. Phosphorylation assay coupled with mutagenesis analyses identifies S122, T284, and S347 in HvWRKY1 as the major residues phosphorylated by HvMPK4. HvWRKY1 is phosphorylated in barley at the early stages of Bgh infection, which enhances its suppression on barley immunity likely due to enhanced DNA-binding and transcriptional repression activity. Our data suggest that the HvMKK1-HvMPK4 kinase pair acts upstream of HvWRKY1 to negatively regulate barley immunity against powdery mildew.

Identification and expression analysis of MAPK cascade gene family in foxtail millet (Setaria italica)

The mitogen-activated protein kinase (MAPK) cascade pathway is a highly conserved plant cell signaling pathway that plays an important role in plant growth and development and stress response. Currently, MAPK cascade genes have been identified and reported in a variety of plants including Arabidopsis thaliana, Oryza sativa, and Triticum aestivum, but have not been identified in foxtail millet (Setaria italica). In this study, a total of 93 MAPK cascade genes, including 15 SiMAPKs, 10 SiMAPKKs and 68 SiMAPKKKs genes, were identified by genome-wide analysis of foxtail millet, and these genes were distributed on nine chromosomes of foxtail millet. Using phylogenetic analysis, we divided the SiMAPKs and SiMAPKKs into four subgroups, respectively, and the SiMAPKKKs into three subgroups (Raf, ZIK, and MEKK). Whole-genome duplication analysis revealed that there are 14 duplication pairs in the MAPK cascade family in foxtail millet, and they are expanded by segmental replication events. Results from quantitative real-time PCR (qRT-PCR) revealed that the expression levels of most SiMAPKs and SiMAPKKs were changed under both exogenous hormone and abiotic stress treatments, with SiMAPK3 and SiMAPKK4-2 being induced under almost all treatments, while the expression of SiMAPKK5 was repressed. In a nutshell, this study will shed some light on the evolution of MAPK cascade genes and the functional mechanisms underlying MAPK cascade genes in response to hormonal and abiotic stress signaling pathways in foxtail millet (Setaria italica).

Genome-wide identification of the mitogen-activated kinase gene family from Limonium bicolor and functional characterization of LbMAPK2 under salt stress

BACKGROUND

Mitogen-activated protein kinases (MAPKs) are ubiquitous signal transduction components in eukaryotes. In plants, MAPKs play an essential role in growth and development, phytohormone regulation, and abiotic stress responses. The typical recretohalophyte Limonium bicolor (Bunge) Kuntze has multicellular salt glands on its stems and leaves; these glands secrete excess salt ions from its cells to mitigate salt damage. The number, type, and biological function of L. bicolor MAPK genes are unknown.

RESULTS

We identified 20 candidate L. bicolor MAPK genes, which can be divided into four groups. Of these 20 genes, 17 were anchored to 7 chromosomes, while LbMAPK18, LbMAPK19, and LbMAPK20 mapped to distinct scaffolds. Structure analysis showed that the predicted protein LbMAPK19 contains the special structural motif TNY in its activation loop, whereas the other LbMAPK members harbor the conserved TEY or TDY motif. The promoters of most LbMAPK genes carry cis-acting elements related to growth and development, phytohormones, and abiotic stress. LbMAPK1, LbMAPK2, LbMAPK16, and LbMAPK20 are highly expressed in the early stages of salt gland development, whereas LbMAPK4, LbMAPK5, LbMAPK6, LbMAPK7, LbMAPK11, LbMAPK14, and LbMAPK15 are highly expressed during the late stages. These 20 LbMAPK genes all responded to salt, drought and ABA stress. We explored the function of LbMAPK2 via virus-induced gene silencing: knocking down LbMAPK2 transcript levels in L. bicolor resulted in fewer salt glands, lower salt secretion ability from leaves, and decreased salt tolerance. The expression of several genes [LbTTG1 (TRANSPARENT TESTA OF GL1), LbCPC (CAPRICE), and LbGL2 (GLABRA2)] related to salt gland development was significantly upregulated in LbMAPK2 knockdown lines, while the expression of LbEGL3 (ENHANCER OF GL3) was significantly downregulated.

CONCLUSION

These findings increase our understanding of the LbMAPK gene family and will be useful for in-depth studies of the molecular mechanisms behind salt gland development and salt secretion in L. bicolor. In addition, our analysis lays the foundation for exploring the biological functions of MAPKs in an extreme halophyte.

Deciphering mitogen activated protein kinase pathway activated during insect attack in Nicotiana attenuata

Plant yields are compromised due to abiotic and biotic stresses. A crucial biotic stress instigated by insect attack, is a major concern that limits crop production. To overcome the deleterious effect of herbivory, pesticides are used but long-term usage of pesticides can be harmful to the environment and human health. Understanding the plants' inherent defense mechanism by interpreting the interaction pattern of defense-related proteins and signalling components and manipulating them to strengthen defense status, is one of the alternative approaches of green biotechnology. During insect attack, host plants initiate innumerable signalling pathways to activate defense response; Mitogen Activated Protein Kinase (MAPK) Pathway is a crucial component of signalling pathway that regulate the expression of downstream defense-related genes. MAPK pathway has three components: MAPKKK, MAPKK and MAPK. Earlier studies have shown participation of SIPK and WIPK (MAPKs) as well as MEK2 (MAPKK) during insect infestation and its association with plant defense. However, information on the third component and elucidation of the complete MAPK pathway are still elusive. Therefore, this study aims to identify the unknown component and decipher MAPK pathway in Nicotiana attenuata involved in plant defense against herbivory by identifying herbivory-inducible MAPKKKs and and their interaction with known partners of the MAPK pathway by docking and MD simulation. The possible pathway was predicted to be MAPKKK Na12134/Na04522-MEK2-SIPK/WIPK. Further, validation of the above interaction by in vitro and in vivo methods is highly recommended.Communicated by Ramaswamy H. Sarma.

Genome-wide identification and functional characterization of borneol dehydrogenases in Wurfbainia villosa

MAIN CONCLUSION

Genome-wide screening of short-chain dehydrogenases/reductases (SDR) family reveals functional diversification of borneol dehydrogenase (BDH) in Wurfbainia villosa. Wurfbainia villosa is an important medicinal plant, the fruits of which accumulate abundant terpenoids, especially bornane-type including borneol and camphor. The borneol dehydrogenase (BDH) responsible for the conversion of borneol to camphor in W. villosa remains unknown. BDH is one member of short-chain dehydrogenases/reductases (SDR) family. Here, a total of 115 classical WvSDR genes were identified through genome-wide screening. These WvSDRs were unevenly distributed on different chromosomes. Seven candidate WvBDHs based on phylogenetic analysis and expression levels were selected for cloning. Of them, four BDHs can catalyze different configurations of borneol and other monoterpene alcohol substrates to generate the corresponding oxidized products. WvBDH1 and WvBDH2, preferred (+)-borneol to (-)-borneol, producing the predominant ( +)-camphor. WvBDH3 yielded approximate equivalent amount of (+)-camphor and (-)-camphor, in contrast, WvBDH4 generated exclusively (+)-camphor. The metabolic profiles of the seeds showed that the borneol and camphor present were in the dextrorotatory configuration. Enzyme kinetics and expression pattern in different tissues suggested WvBDH2 might be involved in the biosynthesis of camphor in W. villosa. All results will increase the understanding of functional diversity of BDHs.

Genome-wide identification and characterization of lipoxygenase genes related to the English grain aphid infestation response in wheat

44 wheat LOX genes were identified by silico genome-wide search method. TaLOX5, 7, 10, 24, 29, 33 were specifically expressed post aphid infestation, indicating their participation in wheat-aphid interaction. In plants, LOX genes play important roles in various biological progresses including seed germination, tuber development, plant vegetative growth and most crucially in plant signal transduction, stress response and plant defense against plant diseases and insects. Although LOX genes have been characterized in many species, the importance of the LOX family in wheat has still not been well understood, hampering further improvement of wheat under stress conditions. Here, we identified 44 LOX genes (TaLOXs) in the whole wheat genome and classified into three subfamilies (9-LOXs, Type I 13-LOXs and Type II 13-LOXs) according to phylogenetic relationships. The TaLOXs belonging to the same subgroup shared similar gene structures and motif organizations. Synteny analysis demonstrated that segmental duplication events mainly contributed to the expansion of the LOX gene family in wheat. The results of protein-protein interaction network (PPI) and miRNA-TaLOXs predictions revealed that three TaLOXs (TaLOX20, 22 and 37) interacted mostly with proteins related to methyl jasmonate (MeJA) signaling pathway. The expression patterns of TaLOXs in different tissues (root, stem, leaf, spike and grain) under diverse abiotic stresses (heat, cold, drought, drought and heat combined treatment, and salt) as well as under diverse biotic stresses (powdery mildew pathogen, Fusarium graminearum and stripe rust pathogen) were systematically analyzed using RNA-seq data. We obtained aphid-responsive candidate genes by RNA-seq data of wheat after the English grain aphid infestation. Aphid-responsive candidate genes, including TaLOX5, 7, 10, 24, 29 and 33, were up-regulated in the wheat aphid-resistant genotype (Lunxuan144), while they were little expressed in the susceptible genotype (Jimai22) during late response (48 h and 72 h) to the English grain aphid infestation. Meanwhile, qRT-PCR analysis was used to validate these aphid-responsive candidate genes. The genetic divergence and diversity of all the TaLOXs in bread wheat and its relative species were investigated by available resequencing data. Finally, the 3D structure of the TaLOX proteins was predicted based on the homology modeling method. This study not only systematically investigated the characteristics and evolutionary relationships of TaLOXs, but also provided potential candidate genes in response to the English grain aphid infestation and laid the foundation to further study the regulatory roles in the English grain aphid infestation of LOX family in wheat and beyond.

Metal tolerance gene family in barley: an in silico comprehensive analysis

Metal-tolerance proteins (MTPs) are divalent cation transporters that play critical roles in metal tolerance and ion homeostasis in plants. However, a comprehensive study of MTPs is still lacking in crop plants. The current study aimed to comprehensively identify and characterize the MTP gene family in barley (Hordeum vulgare, Hv), an important crop. In total, 12 HvMTPs were identified in the barley genome in this study. They were divided into three phylogenetic groups (Zn-cation diffusion facilitator proteins [CDFs], Fe/Zn-CDFs, and Mn-CDFs) and further subdivided into seven groups (G1, G5, G6, G7, G8, G9, and G12). The majority of MTPs were hydrophobic proteins found in the vacuolar membrane. Gene duplication analysis of HvMTPs revealed one pair of segmental-like duplications in the barley genome. Evolutionary analysis suggested that barley MTPs underwent purifying natural selection. Additionally, the HvMTPs were analyzed in the pan-genome sequences of barley (20 accessions), which suggests that HvMTPs are highly conserved in barley evolution. Cis-acting regulatory elements, microRNA target sites, and protein-protein interaction analysis indicated the role of HvMTPs in a variety of biological processes. Expression profiling suggests that HvMTPs play an active role in maintaining barley nutrient homeostasis throughout its life cycle, and their expression levels were not significantly altered by abiotic stresses like cold, drought, or heat. The expression of barley HvMTP genes in the presence of heavy metals such as Zn²⁺, Cu²⁺, As³⁺, and Cd²⁺ revealed that these MTPs were induced by at least one metal ion, implying their involvement in metal tolerance or transportation. The identification and comprehensive investigation of MTP gene family members will provide important gene resources for the genetic improvement of crops for metal tolerance, bioremediation, or biofortification of staple crops.

Knockout of MITOGEN-ACTIVATED PROTEIN KINASE 3 causes barley root resistance against Fusarium graminearum

The roles of mitogen-activated protein kinases (MAPKs) in plant-fungal pathogenic interactions are poorly understood in crops. Here, microscopic, phenotypic, proteomic, and biochemical analyses revealed that roots of independent transcription activator-like effector nuclease (TALEN)-based knockout lines of barley (Hordeum vulgare L.) MAPK 3 (HvMPK3 KO) were resistant against Fusarium graminearum infection. When co-cultured with roots of the HvMPK3 KO lines, F. graminearum hyphae were excluded to the extracellular space, the growth pattern of extracellular hyphae was considerably deregulated, mycelia development was less efficient, and number of appressoria-like structures and their penetration potential were substantially reduced. Intracellular penetration of hyphae was preceded by the massive production of reactive oxygen species (ROS) in attacked cells of the wild-type (WT), but ROS production was mitigated in the HvMPK3 KO lines. Suppression of ROS production in these lines coincided with elevated abundance of catalase (CAT) and ascorbate peroxidase (APX). Moreover, differential proteomic analysis revealed downregulation of several defense-related proteins in WT, and the upregulation of pathogenesis-related protein 1 (PR-1) and cysteine proteases in HvMPK3 KO lines. Proteins involved in suberin formation, such as peroxidases, lipid transfer proteins (LTPs), and the GDSL esterase/lipase (containing "GDSL" aminosequence motif) were differentially regulated in HvMPK3 KO lines after F. graminearum inoculation. Consistent with proteomic analysis, microscopic observations showed enhanced suberin accumulation in roots of HvMPK3 KO lines, most likely contributing to the arrested infection by F. graminearum. These results suggest that TALEN-based knockout of HvMPK3 leads to barley root resistance against Fusarium root rot.

Arbuscular mycorrhizal fungi enhance photosynthesis and drought tolerance by regulating MAPK genes expressions of Populus simonii × P. nigra

Arbuscular mycorrhizal fungi (AMF) promote plants to absorb more water and nutrients and improve their stress resistance. As the main signal transducer, the mitogen-activated protein kinase (MAPK) cascade plays a vital role in drought stress. However, how the MAPK family genes of mycorrhizal plants respond to stress is still not clear. Our study analyzed physiological indexes and expression profiles of MAPK family genes of Populus simonii × P. nigra under two inoculation treatments (inoculated with or without Rhizophagus irregularis) and two water conditions (well-watered or drought stress). The results showed that the stronger photosynthesis of mycorrhizal plants may be mediated by MAPK genes induced by AMF. Mycorrhizal plants showed lower oxidative damage and drought sensitivity. Mycorrhiza downregulated the expression of PsnMAPK7-2, PsnMAPK16-1, PsnMAPK19-2, and PsnMAPK20-2 which negatively regulate drought tolerance and induced specific PsnMAPKs in roots which activate transcription factors to regulate downstream gene expressions, enhancing drought tolerance. This is the first time to identify part of the MAPK gene family of P. simonii × P. nigra at the genome level and study MAPK genes in mycorrhizal forest trees. This is helpful to understand the function of the MAPK gene family in response to drought of mycorrhizal plants and lays a foundation for afforestation by using mycorrhizal saplings.

Unraveling the genetics of polyamine metabolism in barley for senescence-related crop improvement

We explored the polyamine (PA) metabolic pathway genes in barley (Hv) to understand plant development and stress adaptation in Gramineae crops with emphasis on leaf senescence. Bioinformatics and functional genomics tools were utilized for genome-wide identification, comprehensive gene features, evolution, development and stress effects on the expression of the polyamine metabolic pathway gene families (PMGs). Three S-adenosylmethionine decarboxylases (HvSAMDCs), two ornithine decarboxylases (HvODCs), one arginine decarboxylase (HvADC), one spermidine synthase (HvSPDS), two spermine synthases (HvSPMSs), five copper amine oxidases (HvCuAOs) and seven polyamine oxidases (HvPAOs) members of PMGs were identified and characterized in barley. All the HvPMG genes were found to be distributed on all chromosomes of barley. The phylogenetic and comparative assessment revealed that PA metabolic pathway is highly conserved in plants and the prediction of nine H. vulgare miRNAs (hvu-miR) target sites, 18 protein-protein interactions and 961 putative CREs in the promoter region were discerned. Gene expression of HvSAMDC3, HvCuAO7, HvPAO4 and HvSPMS1 was apparent at every developmental stage. SPDS/SPMS gene family was found to be the most responsive to induced leaf senescence. This study provides a reference for the functional investigation of the molecular mechanism(s) that regulate polyamine metabolism in plants as a tool for future breeding decision management systems.

MAPKs trigger antiviral immunity by directly phosphorylating a rhabdovirus nucleoprotein in plants and insect vectors

Signaling by the evolutionarily conserved mitogen-activated protein kinase or extracellular signal-regulated kinase (MAPK/ERK) plays critical roles in converting extracellular stimuli into immune responses. However, whether MAPK/ERK signaling induces virus immunity by directly phosphorylating viral effectors remains largely unknown. Barley yellow striate mosaic virus (BYSMV) is an economically important plant cytorhabdovirus that is transmitted by the small brown planthopper (SBPH, Laodelphax striatellus) in a propagative manner. Here, we found that the barley (Hordeum vulgare) MAPK MPK3 (HvMPK3) and the planthopper ERK (LsERK) proteins interact with the BYSMV nucleoprotein (N) and directly phosphorylate N protein primarily on serine 290. The overexpression of HvMPK3 inhibited BYSMV infection, whereas barley plants treated with the MAPK pathway inhibitor U0126 displayed greater susceptibility to BYSMV. Moreover, knockdown of LsERK promoted virus infection in SBPHs. A phosphomimetic mutant of the N Ser290 (S290D) completely abolished virus infection because of impaired self-interaction of BYSMV N and formation of unstable N-RNA complexes. Altogether, our results demonstrate that the conserved MAPK and ERK directly phosphorylate the viral nucleoprotein to trigger immunity against cross-kingdom infection of BYSMV in host plants and its insect vectors.

An interolog-based barley interactome as an integration framework for immune signaling

The barley MLA nucleotide-binding leucine-rich-repeat (NLR) receptor and its orthologs confer recognition specificity to many fungal diseases, including powdery mildew, stem-, and stripe rust. We used interolog inference to construct a barley protein interactome (Hordeum vulgare predicted interactome, HvInt) comprising 66,133 edges and 7,181 nodes, as a foundation to explore signaling networks associated with MLA. HvInt was compared with the experimentally validated Arabidopsis interactome of 11,253 proteins and 73,960 interactions, verifying that the 2 networks share scale-free properties, including a power-law distribution and small-world network. Then, by successive layering of defense-specific "omics" datasets, HvInt was customized to model cellular response to powdery mildew infection. Integration of HvInt with expression quantitative trait loci (eQTL) enabled us to infer disease modules and responses associated with fungal penetration and haustorial development. Next, using HvInt and infection-time-course RNA sequencing of immune signaling mutants, we assembled resistant and susceptible subnetworks. The resulting differentially coexpressed (resistant - susceptible) interactome is essential to barley immunity, facilitates the flow of signaling pathways and is linked to mildew resistance locus a (Mla) through trans eQTL associations. Lastly, we anchored HvInt with new and previously identified interactors of the MLA coiled coli + nucleotide-binding domains and extended these to additional MLA alleles, orthologs, and NLR outgroups to predict receptor localization and conservation of signaling response. These results link genomic, transcriptomic, and physical interactions during MLA-specified immunity.

Genome-Wide Identification and Expression of MAPK Gene Family in Cultivated Strawberry and Their Involvement in Fruit Developing and Ripening

Studies on many plants have shown that mitogen-activated protein kinases (MAPKs) are key proteins involved in regulating plant responses to biotic and abiotic stresses. However, their involvement in cultivated strawberry development and ripening remains unclear. In this study, 43 FaMAPK gene family members were identified in the genome of cultivated strawberry (Fragaria × ananassa), phylogenetic analysis indicated that FaMAPKs could be classified into four groups. Systematic analysis of the conserved motif, exon-intron structure showed that there were significant varieties between different groups in structure, but in the same group they were similar. Multiple cis-regulatory elements associated with phytohormone response, and abiotic and biotic stresses were predicted in the promoter regions of FaMAPK genes. Transcriptional analysis showed that all FaMAPK genes were expressed at all developmental stages. Meanwhile, the effect of exogenous ABA and sucrose on the expression profile of FaMAPKs was investigated. Exogenous ABA, sucrose, and ABA plus sucrose treatments upregulated the expression of FaMAPK genes and increased the content of endogenous ABA, sucrose, and anthocyanin in strawberry fruits, suggesting that ABA and sucrose might be involved in the FaMAPK-mediated regulation of strawberry fruit ripening. Based on the obtained results, MAPK genes closely related to the ripening of strawberries were screened to provide a theoretical basis and support for future research on strawberries.

Global Identification and Characterization of C2 Domain-Containing Proteins Associated with Abiotic Stress Response in Rice (Oryza sativa L.)

C2 domain-containing proteins (C2DPs) have been identified in different genomes that contain single or multiple C2 domains in their C- or N-terminal. It possesses higher functional activity in the transmembrane regions. The identification of C2 domains were reported in a previous study, such as multiple C2 domains and transmembrane-region proteins (MCTPs) and N-terminal-TM-C2 domain proteins (NTMC2s) of rice, Arabidopsis thaliana, and cotton, whereas the C2DP gene family in rice has not been comprehensively studied, and the role of the C2DP gene in rice in response to abiotic stress is not yet fully understood. In this study, we identified 82 C2DPs in the rice genome and divided them into seven groups through phylogenetic analysis. The synteny analysis revealed that duplication events were either exhibited within the genome of rice or between the genomes of rice and other species. Through the analysis of cis-acting elements in promoters, expression profiles, and qRT-PCR results, the functions of OsC2DPs were found to be widely distributed in diverse tissues and were extensively involved in phytohormones-related and abiotic stresses response in rice. The prediction of the microRNA (miRNA) targets of OsC2DPs revealed the possibility of regulation by consistent miRNAs. Notably, OsC2DP50/51/52 as a co-tandem duplication exhibited similar expression variations and involved the coincident miRNA-regulation pathway. Moreover, the results of the genotypic variation and haplotype analysis revealed that OsC2DP17, OsC2DP29, and OsC2DP49 were associated with cold stress responses. These findings provided comprehensive insights for characterizations of OsC2DPs in rice as well as for their roles for abiotic stress.

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.

MPK homolog GhNTF6 was involved in cotton against Verticillium wilt by interacted with VdEPG1

Mitogen-activated protein kinases (MPKs) are important in regulating plant development and stress response. Rapid activation of MPKs in plants usually depends on its phosphorylated. In view of this situation, a phosphorylated GhNTF6 belonged to MPKs family was screened in cotton roots under Verticillium dahliae challenge by phosphoproteomics analysis. Expression of GhNTF6 in cotton plants was did not induce by V. dahliae infection, while, silencing GhNTF6 results to enhance cotton plants susceptibility to V. dahliae, overexpression - GhNTF6 enhance Arabidopsis plants survivability to V. dahliae. Moreover, the mutation of GhNTF6 at site Thr195 and Thy197 with the phosphorylation decreased the plant resistance to V. dahliae. Therefore, GhNTF6 phosphorylation is important in plants against V. dahliae. Further analysis demonstrated that GhNTF6 interacted with a V. dahliae endopolygalacturonase (VdEPG1) on the cell nucleus. We propose that GhNTF6 is a potential molecular target for improving resistance to Verticillium wilt in cotton.

Network organization of the plant immune system: from pathogen perception to robust defense induction

The plant immune system has been explored essentially through the study of qualitative resistance, a simple form of immunity, and from a reductionist point of view. The recent identification of genes conferring quantitative disease resistance revealed a large array of functions, suggesting more complex mechanisms. In addition, thanks to the advent of high-throughput analyses and system approaches, our view of the immune system has become more integrative, revealing that plant immunity should rather be seen as a distributed and highly connected molecular network including diverse functions to optimize expression of plant defenses to pathogens. Here, we review the recent progress made to understand the network complexity of regulatory pathways leading to plant immunity, from pathogen perception, through signaling pathways and finally to immune responses. We also analyze the topological organization of these networks and their emergent properties, crucial to predict novel immune functions and test them experimentally. Finally, we report how these networks might be regulated by environmental clues. Although system approaches remain extremely scarce in this area of research, a growing body of evidence indicates that the plant response to combined biotic and abiotic stresses cannot be inferred from responses to individual stresses. A view of possible research avenues in this nascent biology domain is finally proposed.

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.

Genome-Wide Identification and Functional Characterization of the Trans-Isopentenyl Diphosphate Synthases Gene Family in Cinnamomum camphora

Trans-isopentenyl diphosphate synthases (TIDSs) genes are known to be important determinants for terpene diversity and the accumulation of terpenoids. The essential oil of Cinnamomum camphora, which is rich in monoterpenes, sesquiterpenes, and other aromatic compounds, has a wide range of pharmacological activities and has therefore attracted considerable interest. However, the TIDS gene family, and its relationship to the camphor tree (C. camphora L. Presl.), has not yet been characterized. In this study, we identified 10 TIDS genes in the genome of the C. camphora borneol chemotype that were unevenly distributed on chromosomes. Synteny analysis revealed that the TIDS gene family in this species likely expanded through segmental duplication events. Furthermore, cis-element analyses demonstrated that C. camphora TIDS (CcTIDS) genes can respond to multiple abiotic stresses. Finally, functional characterization of eight putative short-chain TIDS proteins revealed that CcTIDS3 and CcTIDS9 exhibit farnesyl diphosphate synthase (FPPS) activity, while CcTIDS1 and CcTIDS2 encode geranylgeranyl diphosphate synthases (GGPPS). Although, CcTIDS8 and CcTIDS10 were found to be catalytically inactive alone, they were able to bind to each other to form a heterodimeric functional geranyl diphosphate synthase (GPPS) in vitro, and this interaction was confirmed using a yeast two-hybrid assay. Furthermore, transcriptome analysis revealed that the CcTIDS3, CcTIDS8, CcTIDS9, and CcTIDS10 genes were found to be more active in C. camphora roots as compared to stems and leaves, which were verified by quantitative real-time PCR (qRT-PCR). These novel results provide a foundation for further exploration of the role of the TIDS gene family in camphor trees, and also provide a potential mechanism by which the production of camphor tree essential oil could be increased for pharmacological purposes through metabolic engineering.

Genome-wide identification and expression analysis of U-box gene family in wild emmer wheat (Triticum turgidum L. ssp. dicoccoides)

In this study, 82 U-box genes were identified in wild emmer wheat (TdPUBs) through a genome-search method. Phylogenetic analysis classified them into seven groups and the genes belonging to the same group shared the similar exon-intron structure, motif organization and cis-element compositions. Synteny analysis of the U-box genes between different species revealed that segmental duplication and polyploidization mainly contributed to the expansion of TdPUBs. Furthermore, the genetic variations of U-box were investigated in wild emmer, domesticated emmer and durum wheat. Results showed that significant genetic bottleneck has occurred during domestication process of tetraploid emmer wheat. Meanwhile, 12 TdPUBs were co-located with known domestication related QTLs. Finally, the tissue-specific and stress-responsive TdPUB genes were identified through RNA-seq analysis. Combined with qPCR validation of 19 salt-responsive TdPUBs, the candidates involving in salt response were obtained. It lays the foundation to better understand the regulatory roles of U-box family in emmer wheat and beyond.

Genome-wide identification and characterization of caffeoyl-coenzyme A O-methyltransferase genes related to the Fusarium head blight response in wheat

Background

Lignin is one of the main components of the cell wall and is directly associated with plant development and defence mechanisms in plants, especially in response to Fusarium graminearum (Fg) infection. Caffeoyl-coenzyme A O-methyltransferase (CCoAOMT) is the main regulator determining the efficiency of lignin synthesis and composition. Although it has been characterized in many plants, to date, the importance of the CCoAOMT family in wheat is not well understood.

Results

Here, a total of 21 wheat CCoAOMT genes (TaCCoAOMT) were identified through an in silico genome search method and they were classified into four groups based on phylogenetic analysis, with the members of the same group sharing similar gene structures and conserved motif compositions. Furthermore, the expression patterns and co-expression network in which TaCCoAOMT is involved were comprehensively investigated using 48 RNA-seq samples from Fg infected and mock samples of 4 wheat genotypes. Combined with qRT-PCR validation of 11 Fg-responsive TaCCoAOMT genes, potential candidates involved in the FHB response and their regulation modules were preliminarily suggested. Additionally, we investigated the genetic diversity and main haplotypes of these CCoAOMT genes in bread wheat and its relative populations based on resequencing data.

Conclusions

This study identified and characterized the CCoAOMT family in wheat, which not only provided potential targets for further functional analysis, but also contributed to uncovering the mechanism of lignin biosynthesis and its role in FHB tolerance in wheat and beyond.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07849-y.

Genome-Wide Identification and Genetic Variations of the Starch Synthase Gene Family in Rice

Starch is a major ingredient in rice, and the amylose content of starch significantly impacts rice quality. OsSS (starch synthase) is a gene family related to the synthesis of amylose and amylopectin, and 10 members have been reported. In the present study, a synteny analysis of a novel family member belonging to the OsSSIV subfamily that contained a starch synthase catalytic domain showed that three segmental duplications and multiple duplications were identified in rice and other species. Expression data showed that the OsSS gene family is involved in diverse expression patterns. The prediction of miRNA targets suggested that OsSS are possibly widely regulated by miRNA functions, with miR156s targeted to OsSSII-3, especially. Haplotype analysis exhibited the relationship between amylose content and diverse genotypes. These results give new insight and a theoretical basis for the improved amylose content and eating quality of rice.

TALEN-Based HvMPK3 Knock-Out Attenuates Proteome and Root Hair Phenotypic Responses to flg22 in Barley

Mitogen activated protein kinases (MAPKs) integrate elicitor perception with both early and late responses associated with plant defense and innate immunity. Much of the existing knowledge on the role of plant MAPKs in defense mechanisms against microbes stems from extensive research in the model plant Arabidopsis thaliana. In the present study, we investigated the involvement of barley (Hordeum vulgare) MPK3 in response to flagellin peptide flg22, a well-known bacterial elicitor. Using differential proteomic analysis we show that TALEN-induced MPK3 knock-out lines of barley (HvMPK3 KO) exhibit constitutive downregulation of defense related proteins such as PR proteins belonging to thaumatin family and chitinases. Further analyses showed that the same protein families were less prone to flg22 elicitation in HvMPK3 KO plants compared to wild types. These results were supported and validated by chitinase activity analyses and immunoblotting for HSP70. In addition, differential proteomes correlated with root hair phenotypes and suggested tolerance of HvMPK3 KO lines to flg22. In conclusion, our study points to the specific role of HvMPK3 in molecular and root hair phenotypic responses of barley to flg22.

Genome-Wide Identification and Analysis of MKK and MAPK Gene Families in Brassica Species and Response to Stress in Brassica napus

Mitogen-activated protein kinase (MAPK) cascades are common and conserved signal transduction pathways and play important roles in various biotic and abiotic stress responses and growth and developmental processes in plants. With the advancement of sequencing technology, more systematic genetic information is being explored. The work presented here focuses on two protein families in Brassica species: MAPK kinases (MKKs) and their phosphorylation substrates MAPKs. Forty-seven MKKs and ninety-two MAPKs were identified and extensively analyzed from two tetraploid (B. juncea and B. napus) and three diploid (B. nigra, B. oleracea, and B. rapa) Brassica species. Phylogenetic relationships clearly distinguished both MKK and MAPK families into four groups, labeled A–D, which were also supported by gene structure and conserved protein motif analysis. Furthermore, their spatial and temporal expression patterns and response to stresses (cold, drought, heat, and shading) were analyzed, indicating that BnaMKK and BnaMAPK transcript levels were generally modulated by growth, development, and stress signals. In addition, several protein interaction pairs between BnaMKKs and C group BnaMAPKs were detected by yeast two-hybrid assays, in which BnaMKK3 and BnaMKK9 showed strong interactions with BnaMAPK1/2/7, suggesting that interaction between BnaMKKs and C group BnaMAPKs play key roles in the crosstalk between growth and development processes and abiotic stresses. Taken together, our data provide a deeper foundation for the evolutionary and functional characterization of MKK and MAPK gene families in Brassica species, paving the way for unraveling the biological roles of these important signaling molecules in plants.

CRISPR/Cas9-Induced Loss-of-Function Mutation in the Barley Mitogen-Activated Protein Kinase 6 Gene Causes Abnormal Embryo Development Leading to Severely Reduced Grain Germination and Seedling Shootless Phenotype

The diverse roles of mitogen-activated protein kinases (MAPKs, MPKs) in plant development could be efficiently revealed by reverse genetic studies. In Arabidopsis, mpk6 knockout mutants complete the life cycle; however, ~40% of their embryos show defects in the development leading to abnormal phenotypes of seeds and seedlings' roots. Contrary to the Arabidopsis MPK6, the rice MPK6 (OsMPK6) is an essential gene as transfer DNA (T-DNA) insertion and CRISPR/Cas9 induced loss-of-function mutations in the OsMPK6 cause early embryo arrest. In this study, we successfully developed a viable transgenic barley line with the CRISPR/Cas9-induced heterozygous single base pair cytosine-guanine (CG) deletion [wild type (WT)/-1C] in the third exon of the HvMPK6 gene, a barley ortholog of the Arabidopsis and rice MPK6. There were no obvious macroscopic phenotype differences between the WT/-1C plants and WT plants. All the grains collected from the WT/-1C plants were of similar size and appearance. However, seedling emergence percentage (SEP) from these grains was substantially decreased in the soil in the T2 and T3 generation. The mutation analysis of the 248 emerged T2 and T3 generation plants showed that none of them was a biallelic mutant in the HvMPK6 gene, suggesting lethality of the -1C/-1C homozygous knockout mutation. In the soil, the majority of the -1C/-1C grains did not germinate and the minority of them developed into abnormal seedlings with a shootless phenotype and a reduced root system. Some of the -1C/-1C seedlings also developed one or more small chlorotic leaf blade-like structure/structures. The -1C/-1C grains contained the late-stage developed abnormal embryos with the morphologically obvious scutellum and root part of the embryonic axis but with the missing or substantially reduced shoot part of the embryonic axis. The observed embryonic abnormalities correlated well with the shootless phenotype of the seedlings and suggested that the later-stage defect is predetermined already during the embryo development. In conclusion, our results indicate that barley MPK6 is essential for the embryologically predetermined shoot formation, but not for the most aspects of the embryo and early seedling development.

Transcriptome-wide identification of MAPKKK genes in bermudagrass (Cynodon dactylon L.) and their potential roles in low temperature stress responses

As upstream components of MAPK cascades, mitogen-activated protein kinase kinase kinases (MAPKKKs) act as adaptors linking upstream signaling steps to the core MAPK cascades. MAPK cascades are universal modules of signal transduction in eukaryotic organisms and play crucial roles in plant development processes and in responses to biotic and abiotic stress and signal transduction. Members of the MAPKKK gene family have been identified in several plants,however, MAPKKKs have not been systematically studied in bermudagrass (Cynodon dactylon L.). In this study, 55 potential CdMAPKKKs were produced from bermudagrass transcriptome data, of which 13 belonged to the MEKK, 38 to the Raf, and 4 to the ZIK subfamily. Multiple alignment and conserved motif analysis of CdMAPKKKs supported the evolutionary relationships inferred from phylogenetic analyses. Moreover, the distribution pattern in Poaceae species indicated that members of the MAPKKK family were conserved among almost all diploid species, and species-specific polyploidy or higher duplication ratios resulted in an expansion of the MAPKKK family. In addition, 714 co-functional links which were significantly enriched in signal transduction, responses to temperature stimuli, and other important biological processes of 55 CdMAPKKKs were identified using co-functional gene networks analysis; 30 and 19 co-functional genes involved in response to cold or heat stress, respectively, were also identified. Results of promoter analyses, and interaction network investigation of all CdMAPKKKs based on the rice homologs suggested that CdMAPKKKs are commonly associated with regulation of numerous biological processes. Furthermore, 12 and 13 CdMAPKKKs were significantly up- and downregulated, respectively, in response to low temperature stress; among them, six CdMAPKKKs were significantly induced by low temperature stress, at least at one point in time. This is the first study to conduct identification and functional analysis of the MAPKKK gene family in bermudagrass, and our results provide a foundation for further research on the functions of CdMAPKKKs in response to low temperature stress.

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

BACKGROUND

Mitogen Activated Protein Kinase (MAPK) cascade is a fundamental pathway in organisms for signal transduction. Though it is well characterized in various plants, there is no systematic study of this cascade in tea.

RESULT

In this study, 5 genes of Mitogen Activated Protein Kinase Kinase (MKK) and 16 genes of Mitogen Activated Protein Kinase (MPK) in Camellia sinensis were found through a genome-wide search taking Arabidopsis thaliana as the reference genome. Also, phylogenetic relationships along with structural analysis which includes gene structure, location as well as protein conserved motifs and domains, were systematically examined and further, predictions were validated by the results. The plant species taken for comparative study clearly displayed segmental duplication, which was a significant candidate for MAPK cascade expansion. Also, functional interaction was carried out in C. sinensis based on the orthologous genes in Arabidopsis. The expression profiles linked to various stress treatments revealed wide involvement of MAPK and MAPKK genes from Tea in response to various abiotic factors. In addition, the expression of these genes was analysed in various tissues.

CONCLUSION

This study provides the targets for further comprehensive identification, functional study, and also contributed for a better understanding of the MAPK cascade regulatory network in C. sinensis.