Cyclophilins (CYPs) are a member of the immunophilin superfamily (in addition to FKBPs and parvulins) and play a significant role in peptidyl-prolyl cis-trans isomerase (PPIase) activity. Previous studies have shown that CYPs have important functions in plants, but no genome-wide analysis of the cotton CYP gene family has been reported, and the specific biological function of this gene is still elusive. Based on the release of the cotton genome sequence, we identified 75, 78, 40 and 38 CYP gene sequences from G. barbadense, G. hirsutum, G. arboreum, and G. raimondii, respectively; 221 CYP genes were unequally located on chromosomes. Phylogenetic analysis showed that 231 CYP genes clustered into three major groups and eight subgroups. Collinearity analysis showed that segmental duplications played a significant role in the expansion of CYP members in cotton. There were light-responsiveness, abiotic-stress and hormone-response elements upstream of most of the CYPs. In addition, the motif composition analysis revealed that 49 cyclophilin proteins had extra domains, including TPR (tetratricopeptide repeat), coiled coil, U-box, RRM (RNA recognition motif), WD40 (RNA recognition motif) and zinc finger domains, along with the cyclophilin-like domain (CLD). The expression patterns based on qRT-PCR showed that six CYP expression levels showed greater differences between Xinhai21 (long fibres, G. barbadense) and Ashmon (short fibres, G. barbadense) at 10 and 20 days postanthesis (DPA). These results signified that CYP genes are involved in the elongation stage of cotton fibre development. This study provides a valuable resource for further investigations of CYP gene functions and molecular mechanisms in cotton.

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.

Phenylalanine (Phe), an aromatic amino acid, is a key precursor of flavonoids, which are crucial for plant growth and development. Arogenate dehydratase (ADT) catalyzes the final step in Phe biosynthesis. This study identified eleven ADT genes in G. hirsutum, twelve in G. barbadense, six in G. arboreum, and six in G. raimondii. Among them, GhADT5 exhibited the highest upregulation under alkali stress. Silencing GhADT5 using virus-induced gene silencing (VIGS) reduced cotton tolerance to alkali stress. GhADT5 silencing also led to decreased plant phenylalanine content, total flavonoid content, and activities of superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD). These reductions caused intracellular accumulation of Malondialdehyde (MDA) and reactive oxygen species (ROS). This oxidative damage ultimately reduced tolerance to alkali stress. In addition, silenced plants displayed reduced stomatal aperture, cellular deformation, and irregular intercellular breaks in the leaf epidermis. In summary, these findings suggest that GhADT5 may enhance resistance to alkali stress by regulating enzymatic and non-enzymatic antioxidant systems. This study highlights the role of GhADT5 under alkali stress and provides novel insights for breeding cotton varieties with improved stress tolerance.

Extreme high temperatures are threatening cotton production around the world due to the intensification of global warming. To cope with high-temperature stress, heat-tolerant cotton cultivars have been bred, but the heat-tolerant mechanism remains unclear. This study selected heat-tolerant (‘Xinluzao36′) and heat-sensitive (‘Che61-72′) cultivars of cotton treated with high-temperature stress as plant materials and performed comparative nanopore sequencing transcriptome analysis to reveal the potential heat-tolerant mechanism of cotton. Results showed that 120,605 nonredundant sequences were generated from the raw reads, and 78,601 genes were annotated. Differentially expressed gene (DEG) analysis showed that a total of 19,600 DEGs were screened; the DEGs involved in the ribosome, heat shock proteins, auxin and ethylene signaling transduction, and photosynthesis pathways may be attributed to the heat tolerance of the heat-tolerant cotton cultivar. This study also predicted a total of 5118 long non-coding RNAs (lncRNAs)and 24,462 corresponding target genes. Analysis of the target genes revealed that the expression of some ribosomal, heat shock, auxin and ethylene signaling transduction-related and photosynthetic proteins may be regulated by lncRNAs and further participate in the heat tolerance of cotton. This study deepens our understandings of the heat tolerance of cotton.

Gossypium barbadense possesses a superior fiber quality because of its fiber length and strength. An in-depth analysis of the underlying genetic mechanism could aid in filling the gap in research regarding fiber strength and could provide helpful information for Gossypium barbadense breeding. Three quantitative trait loci related to fiber strength were identified from a Gossypium barbadense recombinant inbred line (PimaS-7 × 5917) for further analysis. RNA sequencing was performed in the fiber tissues of PimaS-7 × 5917 0–35 days postanthesis. Four specific modules closely related to the secondary wall-thickening stage were obtained using the weighted gene coexpression network analysis. In total, 55 genes were identified as differentially expressed from 4 specific modules. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes were used for enrichment analysis, and Gbar_D11G032910, Gbar_D08G020540, Gbar_D08G013370, Gbar_D11G033670, and Gbar_D11G029020 were found to regulate fiber strength by playing a role in the composition of structural constituents of cytoskeleton and microtubules during fiber development. Quantitative real-time PCR results confirmed the accuracy of the transcriptome data. This study provides a quick strategy for exploring candidate genes and provides new insights for improving fiber strength in cotton.

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.