Upland cotton (Gossypium hirsutum) is grown for its elite fiber. Understanding differential gene expression patterns during fiber development will help to identify genes associated with fiber quality. In this study, we used two recombinant inbred lines (RILs) differing in fiber quality derived from an intra-hirsutum population to explore expression profiling differences and identify genes associated with high-quality fiber or specific fiber-development stages using RNA sequencing. Overall, 72/27, 1137/1584, 437/393, 1019/184, and 2555/1479 differentially expressed genes were up-/down-regulated in an elite fiber line (L1) relative to a poor-quality fiber line (L2) at 10, 15, 20, 25, and 30 days post-anthesis, respectively. Three-hundred sixty-three differentially expressed genes (DEGs) between two lines were colocalized in fiber strength (FS) quantitative trait loci (QTL). Short Time-series Expression Miner (STEM) analysis discriminated seven expression profiles; gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation were performed to identify difference in function between genes unique to L1 and L2. Co-expression network analysis detected five modules highly associated with specific fiber-development stages, especially for high-quality fiber tissues. The hub genes in each module were identified by weighted gene co-expression network analysis. Hub genes encoding actin 1, Rho GTPase-activating protein with PAK-box, TPX2 protein, bHLH transcription factor, and leucine-rich repeat receptor-like protein kinase were identified. Correlation networks revealed considerable interaction among the hub genes, transcription factors, and other genes.

Metacaspases (MCs) are cysteine proteases that are important for programmed cell death (PCD) in plants. In this study, we identified 89 MC genes in the genomes of four Gossypium species (Gossypium raimondii, Gossypium barbadense, Gossypium hirsutum, and Gossypium arboreum), and classified them as type-I or type-II genes. All of the type-I and type-II MC genes contain a sequence encoding the peptidase C14 domain. During developmentally regulated PCD, type-II MC genes may play an important role related to fiber elongation, while type-I genes may affect the thickening of the secondary wall. Additionally, 13 genes were observed to be differentially expressed between two cotton lines with differing fiber strengths, and four genes (GhMC02, GhMC04, GhMC07, and GhMC08) were predominantly expressed in cotton fibers at 5–30 days post-anthesis (DPA). During environmentally induced PCD, the expression levels of four genes were affected in the root, stem, and leaf tissues within 6 h of an abiotic stress treatment. In general, the MC gene family affects the development of cotton fibers, including fiber elongation and fiber thickening while four prominent fiber- expressed genes were identified. The effects of the abiotic stress and hormone treatments imply that the cotton MC gene family may be important for fiber development. The data presented herein may form the foundation for future investigations of the MC gene family in Gossypium species.

Microtubules (MTs) are of importance to fiber development. The Xklp2 (TPX2) proteins as a class of microtubule-associated proteins (MAPs) play a key role in plant growth and development by regulating the dynamic changes of microtubules (MTs). However, the mechanism underlying this is unknown. The interactions between TPX2 proteins and tubulin protein, which are the main structural components, have not been studied in fiber development of upland cotton. Therefore, a genome-wide analysis of the TPX2 family was firstly performed in Gossypium hirsutum L. This study identified 41 GhTPX2 sequences in the assembled G. hirsutum genome by a series of bioinformatic methods. Generally, this gene family is phylogenetically grouped into six subfamilies, and 41 G. hirsutum TPX2 genes (GhTPX2s) are distributed across 21 chromosomes. A heatmap of the TPX2 gene family showed that homologous GhTPX2 genes, GhWDLA2/7 and GhWDLA4/9, have large differences in expression levels between two upland cotton recombinant inbred lines (69307 and 69362) that are different in fiber quality at 15 and 20 days post anthesis. The relative data indicate that these four genes are down-regulated under oryzalin, which causes microtubule depolymerization, as determined via qRT-PCR. A subcellular localization experiment suggested that GhWDLA2 and GhWDLA7 are localized to the microtubule cytoskeleton, and GhWDLA4 and GhWDLA9 are only localized to the nucleus. However, only GhWDLA7 between GhWDLA2 and GhWDLA7 interacted with GhTUA2 in the yeast two-hybrid assay. These results lay the foundation for further function study of the TPX2 gene family.