Chalcone isomerase (CHI) is a key component of phenylalanine metabolism that can produce a variety of flavonoids. However, little information and no systematic analysis of CHI genes is available for cotton. Here, we identified 33 CHI genes in the complete genome sequences of four cotton species (Gossypium arboretum L., Gossypium raimondii L., Gossypium hirsutum L., and Gossypium barbadense L.). Cotton CHI proteins were classified into two main groups, and whole-genome/segmental and dispersed duplication events were important in CHI gene family expansion. qRT-PCR and semiquantitative RT-PCR results suggest that CHI genes exhibit temporal and spatial variation and respond to infection with Fusarium wilt race 7. A preliminary model of CHI gene involvement in cotton evolution was established. Pairwise comparison revealed that seven CHI genes showed higher expression in cultivar 06-146 than in cultivar Xinhai 14. Overall, this whole-genome identification unlocks a new approach to the comprehensive functional analysis of the CHI gene family, which may be involved in adaptation to plant pathogen stress.

Anthocyanidin reductase (ANR) is an important regulator of flavonoid metabolism, and proanthocyanidins, the secondary metabolites of flavonoids, play an important role in the response of plants to pathogenic stress. Therefore, in this study, the expression analysis of the ANR gene family of Gossypium barbadense after inoculation with Fusarium oxysporum f. sp. vasinfectum (FOV) was performed at different time points. It was found that Gb_ANR-47 showed significant differences in the disease-resistant cultivar 06-146 and the susceptible cultivar Xinhai 14, as well as in the highest root expression. It was found that the expression of Gb_ANR-47 in the resistant cultivar was significantly higher than that in the susceptible cultivar by MeJA and SA, and different amounts of methyl jasmonate (MeJA) and salicylic acid (SA) response elements were found in the promoter region of Gb_ANR-47. After silencing GbANR-47 in 06-146 material by VIGS technology, its resistance to FOV decreased significantly. The disease severity index (DSI) was significantly increased, and the anthocyanin content was significantly decreased in silenced plants, compared to controls. Our findings suggest that GbANR-47 is a positive regulator of FOV resistance in Gossypium barbadense. The research results provide an important theoretical basis for in-depth analysis of the molecular mechanism of GbANR-47 and improving the anti-FOV of Gossypium barbadense.

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

In this study, GH19 chitinase (Chi) gene family was systematically identified and characterized using genomic assemblies from four cotton species: Gossypium barbadense, G. hirsutum, G. arboreum, and G. raimondii. A suite of analyses was performed, including genome-wide gene identification, physicochemical property characterization of the encoded proteins, subcellular localization prediction, phylogenetic reconstruction, chromosomal mapping, promoter cis-element analysis, and comprehensive expression profiling using transcriptomic data and qRT-PCR (including tissue-specific expression, hormone treatments, and Fusarium oxysporum infection assays). A total of 107 GH19 genes were identified across the four species (35 in G. barbadense, 37 in G. hirsutum, 19 in G. arboreum, and 16 in G. raimondii). The molecular weights of GH19 proteins ranged from 9.9 to 97.3 kDa, and they were predominantly predicted to localize to the extracellular space. Phylogenetic analysis revealed three well-conserved clades within this family. In tetraploid cotton, GH19 genes were unevenly distributed across 12 chromosomes, often clustering in certain regions, whereas in diploid species, they were confined to five chromosomes. Promoter analysis indicated that GH19 gene promoters contain numerous stress- and hormone-responsive motifs, including those for abscisic acid (ABA), ethylene (ET), and gibberellin (GA), as well as abundant light-responsive elements. The expression patterns of GH19 genes were largely tissue-specific; for instance, GbChi23 was predominantly expressed in the calyx, whereas GbChi19/21/22 were primarily expressed in the roots and stems. Overall, this study provides the first comprehensive genomic and functional characterization of the GH19 family in G. barbadense, laying a foundation for understanding its role in disease resistance mechanisms and aiding in the identification of candidate genes to enhance plant defense against biotic stress.

Cottonseed is an important resource for edible oil and protein. Here, we evaluated cottonseed oil, protein, and gossypol contents using traditional chemical methods and near-infrared reflectance spectroscopy (NIRS) in diverse upland cotton (n = 456) and sea island cotton (n = 205) germplasm. In upland cotton, oil content averaged 21.23 ± 3.98% (12.74–43.56%), protein averaged 23.63 ± 4.63% (15.53–49.79%), and gossypol averaged 1.47 ± 0.21 mg/g (0.06–2.07). Correlation analysis showed a significant negative association between oil and protein (r = −0.125, p < 0.01; confirmed by NIRS: r = −0.171, p < 0.01), a significant negative association between protein and gossypol (r = −0.375, p < 0.01), and a significant positive association between oil and gossypol (r = 0.409, p < 0.01). In sea island cotton, oil, protein, and gossypol contents averaged 24.82 ± 6.15% (14.64–41.43%), 25.75 ± 2.04% (18.84–39.00%), and 1.60 ± 0.15 mg/g (1.22–2.08), respectively. The oil–protein association was strongly negative by NIRS (r = −0.744, p < 0.01), whereas correlations with gossypol were weak and not significant by the traditional method. After screening and evaluation, high oil and protein varieties were identified in upland cotton (n = 15) and sea island cotton (n = 6). Fourteen extreme-oil upland materials were further used to examine flowering-date effects on oil accumulation and physiological indicators, showing rapid oil accumulation and a flowering-date-dependent maximum. Finally, qRT-PCR analysis of lipid-metabolism-related candidate genes showed that seven genes were expressed significantly higher in high-oil than in low-oil materials (p < 0.05), peaking at the late stage of oil accumulation. GhDGAT1 and GhDGAT2 showed positive regulatory effects on oil accumulation, whereas GhFAD3 and GhKCR2 showed negative regulatory effects. Collectively, these findings provide compositional benchmarks, clarify trait relationships, and identify candidate genes useful for breeding cotton cultivars with improved seed oil/protein traits.