Genome-wide identification and expression analysis of PYL family genes and functional characterization of GhPYL8D2 under drought stress in Gossypium hirsutum

Cotton is a crucial economic crop, serving as a natural fiber source for the textile industry. However, drought stress poses a significant threat to cotton fiber quality and productivity worldwide. Pyrabactin Resistance 1-Like (PYL) proteins, as abscisic acid (ABA) receptors, play a crucial role in adverse stress responses, but knowledge about the PYLs in cotton remains limited. In our study, we identified 40 GhPYL genes in Gossypium hirsutum through a genome-wide analysis of the cotton genome database. Our analysis revealed that the PYL family formed three distinct subfamilies with typical family characteristics in G. hirsutum. Additionally, through quantitative expression analysis, including transcriptome dataset and qRT-PCR, we found that all GhPYLs were expressed in all tissues of G. hirsutum, and all GhPYLs were differentially expressed under drought stress. Among them, GhPYL4A1, GhPY5D1, GhPY8D2, and a member of the type 2C protein phosphatases clade A family in Gossypium hirsutum (GhPP2CA), GhHAI2D, showed significant differences in expression levels within 12 h after stress treatment. Our protein interaction analysis and BiFC demonstrated the complex regulatory network between GhPYL family proteins and GhPP2CA proteins. We also found that there is an interaction between GhPYL8D2 and GhHAI2D, and through drought treatment of transgenic cotton, we found that GhPYL8D2 played a vital role in the response of G. hirsutum to drought through stomatal control via co-regulation with GhHAI2D. Our findings provide useful insights into the regulation of GhPYL family genes that occur in response to abiotic stresses in cotton.