Integrated Characterization of Cassava (Manihot esculenta) Pectin Methylesterase (MePME) Genes to Filter Candidate Gene Responses to Multiple Abiotic Stresses

Plant pectin methylesterases (PMEs) play crucial roles in regulating cell wall modification and response to various stresses. Members of the PME family have been found in several crops, but there is a lack of research into their presence in cassava (Manihot esculent), which is an important crop for world food security. In this research, 89 MePME genes were identified in cassava that were separated into two types (type-Ⅰ and type-Ⅱ) according to the existence or absence of a pro-region (PMEI domain). The MePME gene members were unevenly located on 17 chromosomes, with 19 gene pairs being identified that most likely arose via duplication events. The MePMEs could be divided into ten sub-groups in type-Ⅰ and five sub-groups in type-Ⅱ. The motif analysis revealed 11 conserved motifs in type-Ⅰ and 8 in type-Ⅱ MePMEs. The number of introns in the CDS region of type-Ⅰ MePMEs ranged between one and two, and the number of introns in type-Ⅱ MePMEs ranged between one and nine. There were 21 type-Ⅰ and 31 type-Ⅱ MePMEs that contained signal peptides. Most of the type-Ⅰ MePMEs had two conserved “RK/RLL” and one “FPSWVS” domain between the pro-region and the PME domain. Multiple stress-, hormone- and tissue-specific-related cis-acting regulatory elements were identified in the promoter regions of MePME genes. A total of five co-expressed genes (MePME1, MePME2, MePME27, MePME65 and MePME82) were filtered from different abiotic stresses via the use of UpSet Venn diagrams. The gene expression pattern analysis revealed that the expression of MePME1 was positively correlated with the degree of cassava postharvest physiological deterioration (PPD). The expression of this gene was also significantly upregulated by 7% PEG and 14 °C low-temperature stress, but slightly downregulated by ABA treatment. The tissue-specific expression analysis revealed that MePME1 and MePME65 generally displayed higher expression levels in most tissues than the other co-expressed genes. In this study, we obtain an in-depth understanding of the cassava PME gene family, suggesting that MePME1 could be a candidate gene associated with multiple abiotic tolerance.