Expansion and diversification of the Glycine max (Gm) ERD15-like subfamily of the PAM2-like superfamily.

Main conclusion: Despite modulating senescence and drought responses, the GmERD15-like subfamily members are differentially induced by multiple stresses and diverge partially in stress signaling functions. The PAM2 motif represents a binding site for poly (A)-binding proteins (PABPs), often associated with RNA metabolism regulation. The PAM2-containing protein ERD15 stands out as a critical regulator of diverse stress responses in plants. Despite the relevance of the PAM2 motif, a comprehensive analysis of the PAM2 superfamily and ERD15-like subfamily in the plant kingdom is lacking. Here, we provide an extensive in silico analysis of the PAM2 superfamily and the ERD15-like subfamily in soybean, using Arabidopsis and rice sequences as prototypes. The Glycine max ERD15-like subfamily members were clustered in pairs, likely originating from DNA-based gene duplication, as the paralogs display high sequence conservation, similar exon/intron genome organization, and are undergoing purifying selection. Complementation analyses of an aterd15 mutant demonstrated that the plant ERD15-like subfamily members are functionally redundant in response to drought, osmotic stress, and dark-induced senescence. Nevertheless, the soybean members displayed differential expression profiles, biochemical activity, and subcellular localization, consistent with functional diversification. The expression profiles of Glyma04G138600 under salicylic acid (SA) and abscisic acid (ABA) treatments differed oppositely from those of the other GmERD15-like genes. Abiotic stress-induced coexpression analysis with soybean PABPs showed that Glyma04G138600 was clustered separately from other GmERD15s. In contrast to the AtERD15 stress-induced nuclear redistribution, Glyma04G138600 and Glyma02G260800 localized to the cytoplasm, while Glyma03G131900 fractionated between the cytoplasm and nucleus under normal and stress conditions. These data collectively indicate that despite modulating senescence and drought responses, the GmERD15-like subfamily members are differentially induced by multiple stresses and may diverge partially in stress signaling functions.