Characterization and diverse evolution patterns of glycerol-3-phosphate dehydrogenase family genes in Dunaliella salina.

The glycerol-3-phosphate dehydrogenase (GPD) gene family plays a major role in glycerol synthesis and adaptation to abiotic stresses. Few studies on GPD family genes from the halotolerant algae Dunaliella salina are available. In this study, seven DsaGPD genes were identified by mining D. salina sequencing data. Among them, DsaGPD5 contained the canonical NAD⁺-GPD protein domain, called si-GPD. In comparison, DsaGPD1-4 not only contained the canonical NAD⁺-GPD domain but also a unique domain, the haloacid dehalogenase (HAD)-like superfamily domain, in their N-terminal region, called bi-GPD. DsaGPD6, 7 contained the FAD⁺-GPD domain. In the transient expression system, DsaGPD1, 3, 4 were found in the cytosol of Arabidopsis thaliana protoplast, DsaGPD2, 5 in the chloroplast, and DsaGPD6, 7 in the mitochondria. MEME analysis showed that six conserved motifs were present in both si-GPDs and bi-GPDs, whereas seven highly conserved motifs were only present in bi-GPDs. The quantitative real-time PCR results showed significant induction of the DsaGPD genes under abiotic stresses, indicating their tolerance-related role in D. salina. DsaGPD2 and DsaGPD5 may be the osmoregulator form and glyceride form in the chloroplast, respectively. The evolutionary forces acting on si-GPDs and bi-GPDs were different in the same organism: bi-GPDs were under purifying selection, while si-GPDs were mainly under positive selection. Furthermore, evolution of the N_HAD domain and C_GPD domain in bi-GPDs is highly correlated. In summary, this study characterizes DsaGPD gene family members and provides useful information for elucidating the salt tolerance mechanism in D. salina.