Background: Serine carboxypeptidase-like protein (SCPL) plays an important role in response to stress in plant. However, our knowledge of the function of the SCPL gene family is limited. Results: In this study, a comprehensive and systematic analysis of SCPL gene family was conducted to explore the phylogeny and evolution of the SCPL gene in Gossypium hirsutum. The phenotype and molecular mechanism of silencing of the Gh_SCPL42 under Verticillium wilt stress was also studied. Our results showed that 96 SCPL genes were observed in genome of G. hirsutum, which distributed on 25 chromosomes and most of them were located in the nucleus. The phylogenetic tree analysis showed that members of SCPL gene family can be divided into three subgroups in G. hirsutum, which are relatively conservative in evolution. SCPL gene has a wide range of tissue expression types in G. hirsutum. Promoter analysis showed that the most cis-acting elements related to MeJA and ABA were contained. Through RNA-seq combined with genotyping, it was found that 11 GhSCPL genes not only had significant expression changes during Verticillium wilt stress but also had differential SNPs in the upstream, downstream, exonic or intronic regions. The expression of these 11 genes in the resistant (Zhongzhimian 2) and susceptible (Junmian 1) materials was further analyzed by qRT-PCR, it was found that 6 genes showed significant expression differences in the two materials. Among them, Gh_SCPL42 has the most obvious expression change. Furthermore, virus-induced gene silencing (VIGS) showed necrosis and yellowing of leaves and significantly higher disease severity index (DSI) and disease severity rate (DSR) values in VIGS plants than in control silenced Gh_SCPL42 plants. Moreover, the expression levels of genes related to the SA and JA pathways were significantly downregulated. These results show that Gh_SCPL42 might improve resistance to Verticillium wilt through the SA and JA pathways in G. hirsutum. Conclusion: In conclusion, our findings indicated that Gh_SCPL42 gene plays an important role in resistance to Verticillium wilt in cotton. It was provided an important theoretical basis for further research on the function of SCPL gene family and the molecular mechanism of resistance to Verticillium wilt in cotton.

The β-amylase (BAM) gene family encodes important enzymes that catalyze the conversion of starch to maltose in various biological processes of plants and play essential roles in regulating the growth and development of multiple plants. So far, BAMs have been extensively studied in Arabidopsis thaliana (A. thaliana). However, the characteristics of the BAM gene family in the crucial economic crop, cotton, have not been reported. In this study, 27 GhBAM genes in the genome of Gossypium hirsutum L (G. hirsutum) were identified by genome-wide identification, and they were divided into three groups according to sequence similarity and phylogenetic relationship. The gene structure, chromosome distribution, and collinearity of all GhBAM genes identified in the genome of G. hirsutum were analyzed. Further sequence alignment of the core domain of glucosyl hydrolase showed that all GhBAM family genes had the glycosyl hydrolase family 14 domain. We identified the BAM gene GhBAM7 and preliminarily investigated its function by transcriptional sequencing analysis, qRT-PCR, and subcellular localization. These results suggested that the GhBAM7 gene may influence fiber strength during fiber development. This systematic analysis provides new insight into the transcriptional characteristics of BAM genes in G. hirsutum. It may lay the foundation for further study of the function of these genes.

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.

Fuzzless Gossypium hirsutum mutants are ideal materials for investigating cotton fiber initiation and development. In this study, we used the fuzzless G. hirsutum mutant Xinluzao 50 FLM as the research material and combined it with other fuzzless materials for verification by RNA sequencing to explore the gene expression patterns and differences between genes in upland cotton during the fuzz period. A gene ontology (GO) enrichment analysis showed that differentially expressed genes (DEGs) were mainly enriched in the metabolic process, microtubule binding, and other pathways. A weighted gene co-expression network analysis (WGCNA) showed that two modules of Xinluzao 50 and Xinluzao 50 FLM and four modules of CSS386 and Sicala V-2 were highly correlated with fuzz. We selected the hub gene with the highest KME value among the six modules and constructed an interaction network. In addition, we selected some genes with high KME values from the six modules that were highly associated with fuzz in the four materials and found 19 common differential genes produced by the four materials. These 19 genes are likely involved in the formation of fuzz in upland cotton. Several hub genes belong to the arabinogalactan protein and GDSL lipase, which play important roles in fiber development. According to the differences in expression level, 4 genes were selected from the 19 genes and tested for their expression level in some fuzzless materials. The modules, hub genes, and common genes identified in this study can provide new insights into the formation of fiber and fuzz, and provide a reference for molecular design breeding for the genetic improvement of cotton fiber.

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT–PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Cotton fiber quality improvement remains a fundamental challenge in breeding programs due to the complex genetic architecture underlying fiber development. The narrow genetic base of upland cotton (Gossypium hirsutum L.) and the quantitative nature of fiber quality traits necessitate innovative approaches for identifying and incorporating superior alleles from related species. We developed a BC₆F₂ population by introgressing chromosome segments from the sea island cotton variety Xinhai 36 (G. barbadense) into the upland cotton variety Xinluzhong 60 (G. hirsutum). Based on fiber strength phenotyping, we constructed two DNA bulks representing extreme phenotypes (20 superior and 12 inferior individuals) for bulked segregant analysis sequencing (BSA-Seq). High-throughput sequencing generated 225.13 Gb of raw data with average depths of 20× for parents and 30× for bulks. SNP calling and annotation were performed using GATK and ANNOVAR against the upland cotton reference genome (TM-1). BSA-Seq analysis identified 13 QTLs primarily clustered within a 1.6 Mb region (20.6–22.2 Mb) on chromosome A10. Within this region, we detected nonsynonymous mutation genes involving a total of six genes. GO and KEGG enrichment analyses revealed significant enrichment for carbohydrate metabolic processes, protein modification, and secondary metabolite biosynthesis pathways. Integration with transcriptome data prioritized GH_A10G1043, encoding a β-amylase family protein, as the key candidate gene. Functional validation through overexpression and RNAi knockdown in Arabidopsis thaliana demonstrated that GH_A10G1043 significantly regulates starch content and β-amylase activity, though without visible morphological alterations. This study successfully identified potential genomic regions and candidate genes associated with cotton fiber strength using chromosome segment substitution lines combined with BSA-Seq. The key candidate gene GH_A10G1043 provides a valuable target for marker-assisted selection in cotton breeding programs. Our findings establish a foundation for understanding the molecular mechanisms of fiber quality formation and offer genetic resources for developing superior cotton varieties with enhanced fiber strength.