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.

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.

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

Background: DNA methylation is an important part of epigenetic regulation and plays an important role in the response of plants to adverse stress. Methods: In this study, whole-genome bisulfite sequencing (WGBS) was performed on the high-temperature-resistant material Xinluzao 36 and the high-temperature-sensitive material Che 61–72 at 0 h and 12 h under high-temperature stress conditions. Results: The results revealed that the Gossypium hirsutum methylation levels of CG and CHG (H = A, C, or T) decreased after the high-temperature stress treatment, and the methylation level of the A subgenome was significantly greater than that of the D subgenome. The methylation level of CHH increased, and the methylation level of CHH in the D subgenome was significantly greater than that in the A subgenome after high-temperature stress treatment. The methylation density of CG is lower than that of CHG and CHH, and the methylation density of the middle region of chromosomes is greater than that of both ends, which is opposite to the distribution density of genes. There were 124 common differentially methylated genes in the CG, CHG, and CHH groups, and 5130 common DEGs and differentially methylated genes were found via joint analysis with RNA-seq; these genes were significantly enriched in the biosynthesis of plant hormones, thiamine metabolism, glutathione metabolism, and tyrosine metabolism pathways. DNA methylation did not affect the expression of many genes (accounting for 85.68% of the differentially methylated genes), DNA methylation-promoted gene expression was located mainly in the downstream region of the gene or gene body, and the expression of inhibitory genes was located mainly in the upstream region of the gene. Conclusions: This study provides a theoretical basis for further exploration of the gene expression and functional regulatory mechanism of G. hirsutum DNA methylation under high-temperature stress conditions.

The occurrence and spread of Verticillium dahliae (V. dahliae) in cotton depends on the combined effects of pathogens, host plants, and the environment, among which temperature is one of the most important environmental factors. Studying how temperature impacts the occurrence of V. dahliae in cotton and the mechanisms governing host defense responses is crucial for disease prevention and control. Understanding the dual effects of temperature on both pathogens and hosts can provide valuable insights for developing effective strategies to manage this destructive fungal infection in cotton. This study was based on the deciduous V. dahliae Vd-3. Through cultivation at different temperatures, Vd-3 formed the most microsclerotia and had the largest colony diameter at 25 °C. Endospore toxins were extracted, and 48 h was determined to be the best pathogenic time point for endotoxins to infect cotton leaves through a chlorophyll fluorescence imaging system and phenotypic evaluation. Transcriptome sequencing was performed on cotton leaves infected with Vd-3 endotoxins for 48 h at different culture temperatures. A total of 34,955 differentially expressed genes (DEGs) were identified between each temperature and CK (no pathogen inoculation), including 17,422 common DEGs. The results of the enrichment analysis revealed that all the DEGs were involved mainly in photosynthesis and sugar metabolism. Among the 34,955 DEGs, genes in the biosynthesis and signal transduction pathways of jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) were identified, and their expression patterns were determined. A total of 5652 unique DEGs were clustered into six clusters using the k-means clustering algorithm, and the functions and main transcription factors (TFs) of each cluster were subsequently annotated. In addition, we constructed a gene regulatory network via weighted correlation network analysis (WGCNA) and identified twelve key genes related to cotton defense against V. dahliae at different temperatures, including four genes encoding transcription factors. These findings provide a theoretical foundation for investigating temperature regulation in V. dahliae infecting cotton and introduce novel genetic resources for enhancing resistance to this disease in cotton plants.

Background: In recent years, changes in climate conditions and long-term continuous cropping have led to the increased occurrence of Verticillium wilt in various cotton-growing regions, causing significant economic losses in cotton production. Research has shown that volatile substances are closely linked to plant disease resistance; however, studies on their roles in the response of cotton to Verticillium wilt, including their relationship with gene regulation, are limited. Methods: In this study, the transcriptomes and metabolomes of Xinluzao 57 (a highly susceptible Verticillium wilt variety) and 192,868 (a highly resistant Verticillium wilt variety) were sequenced at different time points after inoculation with Verticillium wilt. Results: A total of 21,911 commonly differentially expressed genes (DEGs) were identified within and between the materials, and they were clustered into eight groups. Significant annotations were made in pathways related to amino acids and anthocyanins. Metabolomics identified and annotated 26,200 volatile metabolites across nine categories. A total of 158 differentially accumulated metabolites (DAMs) were found within and between the materials; three clusters were identified, and the 10 metabolites with the most significant fold changes were highlighted. Weighted gene coexpression network analysis (WGCNA) revealed that 13 genes were significantly correlated with guanosine, 6 genes were correlated with 2-deoxyerythritol, and 32 genes were correlated with raffinose. Conclusions: Our results provide a foundation for understanding the role of volatile substances in the response of cotton to Verticillium wilt and offer new gene resources for future research on Verticillium wilt resistance.