Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

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.

Introduction: Cotton, one of the most important economic crops worldwide, has long been bred mainly for improvements in yield and quality, with relatively little focus on salt–alkali resistance.Methods: In this study, transcriptomic and metabolomic sequencing were performed on Gossypium hirsutum exposed to alkaline stress for different durations.Results: The results of sample clustering, principal component analysis (PCA), and the number of differentially expressed genes (DEGs) revealed that 12 hours and 24 hours were the periods during which upland cotton presented the strongest response to salt stress, with flavonoid biosynthesis and alpha-linolenic acid metabolism playing significant roles during this time. A total of 6,610 DEGs were identified via comparison to the 0 h time point, including 579 transcription factors (TFs) that were significantly enriched in pathways such as flavonoid biosynthesis, the cell cycle, the cytochrome P450 pathway, phenylalanine metabolism, phototransduction, and alpha-linolenic acid metabolism. Through ultrahigh-performance liquid chromatography–MS (UPLC-MS), 4,225 metabolites were identified, and 1,684 differentially accumulated metabolites (DAMs) were identified by comparison to the levels at 0 h. A joint analysis of RNA-seq and metabolomic data revealed that the flavonoid biosynthesis and alpha-linolenic acid metabolism pathways play key roles in the response of G. hirsutum to alkaline stress, and the key genes in these pathways were identified. The weighted gene correlation network analysis (WGCNA) revealed 15 candidate genes associated with alkali tolerance in cotton, including 4 TFs and 4 genes related to flavonoid and anthocyanin biosynthesis.Conclusion: In conclusion, our study provides a theoretical foundation for understanding the molecular mechanisms underlying alkali tolerance in cotton and offers new gene resources for future research.

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.

Drought stress significantly affects the growth, development, and yield of cotton, triggering the response of multiple genes. Among them, ascorbate peroxidase (APX) is one of the important antioxidant enzymes in the metabolism of reactive oxygen species in plants, and APX enhances the ability of plants to resist oxidation, thus increasing plant stress tolerance. Therefore, enhancing the activity of APX in cells is crucial to improving plant stress resistance. Previous studies have isolated differentially expressed proteins under drought stress (GhAPX7) in drought-resistant (KK1543) and drought-sensitive (XLZ26) plants. Thus, this study analyzed the expression patterns of GhAPX7 in different cotton tissues to verify the drought resistance function of GhAPX7 and explore its regulatory pathways. GhAPX7 had the highest expression in cotton leaves, which significantly increased under drought stress, suggesting that GhAPX7 is essential for improving antioxidant capacity and enzyme activities in cotton. GhAPX7 silencing indirectly affects pronounced leaf yellowing and wilting in drought-resistant and drought-sensitive plants under drought stress. Malondialdehyde (MDA) content was significantly increased and chlorophyll and proline content and APX enzyme activity were generally decreased in silenced plants compared to the control. This result indicates that GhAPX7 may improve drought resistance by influencing the contents of MDA, chlorophyll, proline, and APX enzyme activity through increased expression levels. Transcriptome analysis revealed that the drought-related differentially expressed genes between the control and treated groups enriched plant hormone signal transduction, MAPK signaling, and plant–pathogen interaction pathways. Therefore, the decreased expression of GhAPX7 significantly affects the expression levels of genes in these three pathways, reducing drought resistance in plants. This study provides insights into the molecular mechanisms of GhAPX7 and its role in drought resistance and lays a foundation for further research on the molecular mechanisms of response to drought stress in cotton.