Histone deacetylases (HDACs) play a significant role in a plant's development and response to various environmental stimuli by regulating the gene transcription. However, HDACs remain unidentified in cotton. In this study, a total of 29 HDACs were identified in allotetraploid Gossypium hirsutum, while 15 and 13 HDACs were identified in Gossypium arboretum and Gossypium raimondii, respectively. Gossypium HDACs were classified into three groups (reduced potassium dependency 3 (RPD3)/HDA1, HD2-like, and Sir2-like (SRT) based on their sequences, and Gossypium HDACs within each subgroup shared a similar gene structure, conserved catalytic domains and motifs. Further analysis revealed that Gossypium HDACs were under a strong purifying selection and were unevenly distributed on their chromosomes. Gene expression data revealed that G. hirsutum HDACs were differentially expressed in various vegetative and reproductive tissues, as well as at different developmental stages of cotton fiber. Furthermore, some G. hirsutum HDACs were co-localized with quantitative trait loci (QTLs) and single-nucleotide polymorphism (SNPs) of fiber-related traits, indicating their function in fiber-related traits. We also showed that G. hirsutum HDACs were differentially regulated in response to plant hormones (abscisic acid (ABA) and auxin), DNA damage agent (methyl methanesulfonate (MMS)), and abiotic stresses (cold, salt, heavy metals and drought), indicating the functional diversity and specification of HDACs in response to developmental and environmental cues. In brief, our results provide fundamental information regarding G. hirsutum HDACs and highlight their potential functions in cotton growth, fiber development and stress adaptations, which will be helpful for devising innovative strategies for the improvement of cotton fiber and stress tolerance.
Keywords: abiotic stress; cotton; fiber; genome-wide analysis; histone deacetylases.