Genome-wide identification and analysis of abiotic stress responsiveness of the mitogen-activated protein kinase gene family in Medicago sativa L

The mitogen-activated protein kinase (MAPK) cascade is crucial cell signal transduction mechanism that plays an important role in plant growth and development, metabolism, and stress responses. The MAPK cascade includes three protein kinases, MAPK, MAPKK, and MAPKKK. The three protein kinases mediate signaling to downstream response molecules by sequential phosphorylation. The MAPK gene family has been identified and analyzed in many plants, however it has not been investigated in alfalfa. In this study, Medicago sativa MAPK genes (referred to as MsMAPKs) were identified in the tetraploid alfalfa genome. Eighty MsMAPKs were divided into four groups, with eight in group A, 21 in group B, 21 in group C and 30 in group D. Analysis of the basic structures of the MsMAPKs revealed presence of a conserved TXY motif. Groups A, B and C contained a TEY motif, while group D contained a TDY motif. RNA-seq analysis revealed tissue-specificity of two MsMAPKs and tissue-wide expression of 35 MsMAPKs. Further analysis identified MsMAPK members responsive to drought, salt, and cold stress conditions. Two MsMAPKs (MsMAPK70 and MsMAPK75) responds to salt and cold stresses; two MsMAPKs (MsMAPK60 and MsMAPK73) responds to cold and drought stresses; four MsMAPKs (MsMAPK1, MsMAPK33, MsMAPK64 and MsMAPK71) responds to salt and drought stresses; and two MsMAPKs (MsMAPK5 and MsMAPK7) responded to all three stresses. This study comprehensively identified and analysed the alfalfa MAPK gene family. Candidate genes related to abiotic stresses were screened by analysing the RNA-seq data. The results provide key information for further analysis of alfalfa MAPK gene functions and improvement of stress tolerance.