Genome-wide survey and evolutionary history of the pectin methylesterase (PME) gene family in the Dothideomycetes class of fungi

Abstract

Once deposited in the plant cell wall, pectin undergoes demethylesterification by endogenous pectin methylesterases (PMEs), which play various roles in growth and development, including defense against pathogen attacks. Pathogen PMEs can alter pectin's methylesterification pattern, increasing its susceptibility to degradation by other fungal pectinases and thus playing a critical role as virulence factors during early infection stages. To investigate the evolutionary history of PMEs in the Dothideomycetes class of fungi, we obtained genomic data from 15 orders (79 species) and added genomic data from 61 isolates of Corynespora cassiicola. Our analyses involved maximum likelihood phylogenies, gene genealogies, and selection analyses. Additionally, we measured PME gene expression levels of C. cassiicola using soybean as a host through RT-qPCR assays. We recovered 145 putative effector PMEs and 57 putative non-effector PMEs from across the Dothideomycetes. The PME gene family exhibits a small size (up to 5 members per genome) and comprises three major clades. The evolutionary patterns of the PME1 and PME2 clades were largely shaped by duplications and recurring gene retention events, while biased gene loss characterized the small-sized PME3 clade. The presence of five members in the PME gene family of C. cassiicola suggests that the family may play a key role in the evolutionary success of C. cassiicola as a polyphagous plant pathogen. The haplogroups Cc_PME1.1 and Cc_PME1.2 exhibited an accelerated rate of evolution, whereas Cc_PME2.1, Cc_PME2.2, and Cc_PME2.3 seem to be under strong purifying selective constraints. All five PME genes were expressed during infection of soybean leaves, with the highest levels during from six to eight days post-inoculation. The highest relative expression level was measured for CC_29_g7533, a member of the Cc_PME2.3 clade, while the remaining four genes had relatively lower levels of expression.