Genome-wide identification and molecular evolution analysis of the heat shock transcription factor (HSF) gene family in four diploid and two allopolyploid Gossypium species

Heat shock transcription factors (HSFs) can regulate plant development and stress response. The comprehensive evolutionary history of the HSF family remains elusive in cotton. In this study, each cotton species had 78 members in Gossypium barbadense and Gossypium hirsutum. The diploid species had 39 GaHSFs in Gossypium arboreum, 31 GrHSFs in Gossypium raimondii, 34 GtHSFs in Gossypium turneri, and 34 GlHSFs in Gossypium longicalyx. The HSF family in cotton can be classified into three subfamilies, with seven groups in subfamily A and five groups in subfamily B. Different groups exhibited distinct gene proportions, conserved motifs, gene structures, expansion rates, gene loss rates, and cis-regulatory elements. The paleohexaploidization event led to the expansion of the HSF family in cotton, and the gene duplication events in six Gossypium species were inherited from their common ancestor. The HSF family in diploid species had a divergent evolutionary history, whereas two cultivated tetraploids presented a highly conserved evolution of the HSF family. The HSF members in At and Dt subgenomes of the cultivated tetraploids showed a different evolution from their corresponding diploid donors. Some HSF members were regarded as key candidates for regulating cotton development and stress response. This study provided the comprehensive information on the evolutionary history of the HSF family in cotton.