Integrated Differential Expression Analysis and WGCNA Identify Hub Genes Underlying Cotton Plant Height Development

Plant height is a key agronomic trait that influences plant architecture and mechanical harvesting suitability in cotton; however, the molecular mechanisms underlying its dynamic development remain unclear. In this study, two recombinant inbred line (RIL) populations sharing CCRI127 as a common paternal parent (RIL-GH07, n = 150; RIL-2358B, n = 276) were developed. Based on stable plant-height performance across multiple environments, tall and short extreme lines were selected from the two RIL populations for transcriptome sequencing. By integrating differential expression analysis with weighted gene co-expression network analysis (WGCNA), we identified hub genes associated with cotton plant height development, characterized the molecular features and core pathways governing dynamic stem elongation at different growth stages, thereby providing insights into the transcriptional regulation of plant height development in cotton. The two RIL populations showed broadly similar plant-height growth patterns, with slow elongation at 15 DOS, rapid elongation during 30-60 DOS, and reduced growth after 70 DOS. Transcriptome differential expression analysis identified 15,052 non-redundant DEGs, which exhibited clear population- and stage-specific expression patterns. In the GH07 population, the largest number of DEGs was detected at 15 DOS (7193), whereas in the 2358B population relatively large numbers of DEGs were maintained at both 30 DOS (3839) and 70 DOS (3118). Analysis of DEGs shared by the two populations across four developmental stages showed that, in addition to genes with consistent expression trends, each stage also contained a substantial number of DEGs with opposite expression directions. WGCNA identified 25 gene expression modules, among which the green and yellow modules were significantly positively correlated with plant height. Functional enrichment analysis indicated that genes in these two modules were mainly enriched in hormone regulation and signal transduction, protein modification and degradation, and intracellular transport. Seven hub genes were identified by integrating intramodular connectivity and kME values. Functional prediction suggested that these genes may play important roles in cotton plant height development. This study provides genetic resources and a theoretical basis for subsequent functional validation of cotton plant height-related genes and the improvement of plant architecture in cotton.