Comprehensive genome-wide identification and functional characterization of MAPK cascade gene families in Nelumbo

Mitogen-activated protein kinase (MAPK) cascade signaling pathway plays pivotal roles in various plant biological processes. However, systematic study of MAPK cascade gene families is yet to be conducted in lotus. Herein, 198 putative MAPK genes, including 152 MAP3Ks, 15 MKKs, and 31 MPKs genes were identified in Nelumbo. Segmental duplication was identified as the predominant factor driving MAPK cascade gene family expansion in lotus. MAPK cascade genes in N. nucifera and N. lutea shared high degree of sequence homologies, with 84, 9, and 19 homologous MAP3K, MKK, and MPK gene pairs being detected between the two species, respectively, with most genes predominantly undergoing purifying selection. Gene expression profiling indicated that NnMAPK cascade genes were extensively involved in plant development and submergence stress response. Co-expression analysis revealed potential interaction between transcription factors (TFs) and NnMAPK cascade genes in various biological processes. NnMKK showed predicted interactions with multiple NnMAP3K or NnMPK proteins, which suggested that functional diversity of MAPK cascade genes could be as a result of their complex protein interaction mechanisms. This first systematic analysis of MAPK cascade families in lotus provides deeper insights into their evolutionary dynamics and functional properties, which potentially could be crucial for lotus genetic improvement.