Genome-wide identification and analysis of the evolution and expression patterns of the GATA transcription factors in three species of Gossypium genus

GATA transcription factors (TFs), which bind to DNA in regulatory regions, are involved in cell differentiation and possess a type-IV zinc finger and a DNA-binding domain. GATA genes have been characterized in plant species such as Arabidopsis thaliana, Oryza sativa, and Glycine max, and their functions have been elucidated in A. thaliana. Although many Gossypium quantitative trait loci for fiber quality harbor GATA TFs, GATA genes have not yet been characterized in cotton. In this study, we identified 179 GATA genes from the genomes of three Gossypium species. We analyzed the phylogenetic relationships, chromosomal distribution, gene structure, expression pattern, and predicted promoters of all 179 Gossypium GATA genes (46 in G. raimondii, 46 in G. arboreum, and 87 in G. hirsutum). Phylogenetic analysis grouped the 179 GATA genes into four subfamilies. Domain analysis revealed that GATA domains in subfamilies I, II, and III were located near the C-terminal, whereas those in subfamily IV were adjacent to the N-terminal. RNA-seq and (Real-time PCR) qRT-PCR revealed that 39.1% (34/87) of GATA genes were expressed in growing plant tissues in G. hirsutum, but only 12.6% (11/87) were expressed during fiber development. In addition, 45.7% (21/46) and 26.1% (12/46) of GATA genes were expressed in G. arboreum and G. raimondii, respectively. Our results may be useful for elucidating the evolution, expression patterns, and functional divergence of GATA genes in Gossypium.

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 d_N:d_S value for the CSLJ clade was 1.305, suggesting this subfamily was under positive selection pressure. Our data indicated that the d_N:d_S 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.