Genome-wide identification and analysis of the evolution and expression patterns of the cellulose synthase gene superfamily in Gossypium species

The cellulose synthase gene superfamily, which includes the cellulose synthase (Ces) and cellulose synthase-like (Csl) families, is involved in the synthesis of cellulose and hemicellulose. This superfamily is critical for cotton fiber development in Gossypium species. Applying a series of bioinformatic methods, we identified 228 Ces/Csl genes from four Gossypium species (G. hirsutum, G. barbadense, G. arboreum, and G. raimondii). These genes were then grouped into 11 subfamilies based on phylogenetic relationships. A subsequent analysis of gene evolution revealed sites in CSLG and CSLJ genes that were under long-term positive selection pressure, with a posterior probability >0.95. Moreover, the dN:dS value for the CSLJ clade was 1.305, suggesting this subfamily was under positive selection pressure. Our data indicated that the dN:dS value ranged from 0.0084 to 0.9693 among the homologous Ces/Csl genes, implying they were under purifying selection pressure. Our transcriptome and qRT-PCR analyses revealed that CesA genes were more highly expressed in tetraploids than in diploids. However, the Csl expression levels exhibited the opposite trend. Furthermore, changes to promoter sequences may have influenced the expression of homologous Ces/Csl genes. Our findings may provide novel insights into the evolutionary relationships and expression patterns of the Ces/Csl genes in Gossypium species.