Cotton fiber as a model for understanding shifts in cell development under domestication

Bifunctional transcription factor GhMYB4 orchestrates transition from elongation to secondary cell wall synthesis trade-off in cotton fiber

The integration of cell wall expansion and reinforcement is vital for plant cell wall development. Cotton fiber is characterized by the synchronized development of fiber growth and cell wall formation, providing an excellent system for investigating plant cell walls. Here, we found that the often-overlooked transition stage coordinates fiber elongation and secondary cell wall (SCW) synthesis through antagonistic effects of the transcription factor GhMYB4. Knockout of GhMYB4 produced longer and finer fibers, contrasting with overexpression of GhMYB4. We show that GhMYB4 represses fiber elongation through fatty acid and brassinosteroid pathways involving a negative feedback loop with GhRAP2 and GhTCP15, while activating SCW synthesis by inducing cellulose biosynthesis in succession with the GhMYB52-GhFSN1-GhILR3 module. We identify that the bifunctionality of GhMYB4 depends on its interaction with different cis-elements and is executed through distinct transcriptional regulation motifs. Our findings propose a strategy to improve fiber quality by orchestrating wall expansion and stiffness.

A GDSL esterase/lipase gene GbGELP identified from a fiber micronaire QTL qMIC-A11 modulates cell elongation and fiber development

KEY MESSAGE

A fiber micronaire QTL qMIC-A11 was fine-mapped, and the GDSL esterase/lipase gene GbGELP was identified as the causal gene of the QTL. GbGELP modulates cell elongation and cotton fiber development. Fine mapping and map-based cloning of fiber micronaire (MIC)-related quantitative trait loci (QTL) have not been reported to date. Here, we utilized a G. hirsutum (Gh) acc. TM-1-G. barbadense (Gb) acc. Hai7124 introgression line CSSL47, which exhibits a significant decrease in MIC compared to TM-1, to cross with TM-1 and develop the F2 and F2:3 secondary segregating populations. Further, a stable MIC QTL qMIC-A11 was simultaneously detected in the F2 and F2:3 populations and anchored within a 407 kb region. Among them, GB_A11G1593 encoding a GDSL esterase/lipase, exhibited substantially higher expression levels at fiber elongation period in CSSL47 compared to TM-1, which was temporally identified as the causal gene for qMIC-A11 and named as GbGELP. The heterologous expression of GbGELP in Arabidopsis showed increased root length, root cell length, rosette leaf growth, and trichome density. However, knockdown of GbGELP homologs in CSSL47 significantly decreased the fiber length. Further investigation found that there was an A/T single-nucleotide polymorphism variation (SNP) in the promoter of GELP orthologs between CSSL47 and TM-1, which results in a differential CATTAAATT/CATTTAATT HAHR1-box cis-acting element, a binding site for the homeodomain-leucine zipper IV (HD-ZIP IV). GbGELP was regulated by a HD-ZIP IV transcription factor GhHDG2 via binding to the CATTAAATT element in the GbGELP promoter, while GhGELP could be activated due to GhHDG2 unable to bind the CATTTAATT element in the GhGELP promoter. The fine-mapped MIC QTL qMIC-A11, along with the causal gene GbGELP, will be utilized to improve the fiber quality in cotton breeding.

Genome-Wide Identification and Classification of Arabinogalactan Proteins Gene Family in Gossypium Species and GhAGP50 Increases Numbers of Epidermal Hairs in Arabidopsis

Arabinogalactan proteins (AGPs) constitute a diverse class of hydroxyproline-rich glycoproteins implicated in various aspects of plant growth and development. However, their functional characterization in cotton (Gossypium spp.) remains limited. As a globally significant economic crop, cotton serves as the primary source of natural fiber, making it essential to understand the genetic mechanisms underlying its growth and development. This study aims to perform a comprehensive genome-wide identification and characterization of the AGP gene family in Gossypium spp., with a particular focus on elucidating their structural features, evolutionary relationships, and functional roles. A genome-wide analysis was conducted to identify AGP genes in Gossypium spp., followed by classification into distinct subfamilies based on sequence characteristics. Protein motif composition, gene structure, and phylogenetic relationships were examined to infer potential functional diversification. Subcellular localization of a key candidate gene, GhAGP50, was determined using fluorescent protein tagging, while gene expression patterns were assessed through β-glucuronidase (GUS) reporter assays. Additionally, hormonal regulation of GhAGP50 was investigated via treatments with methyl jasmonate (MeJA), abscisic acid (ABA), indole-3-acetic acid (IAA), and gibberellin (GA). A total of 220 AGP genes were identified in Gossypium spp., comprising 19 classical AGPs, 28 lysine-rich AGPs, 55 AG peptides, and 118 fasciclin-like AGPs (FLAs). Structural and functional analyses revealed significant variation in gene organization and conserved motifs across subfamilies. Functional characterization of GhAGP50, an ortholog of AGP18 in Arabidopsis thaliana, demonstrated its role in promoting epidermal hair formation in leaves and stalks. Subcellular localization studies indicated that GhAGP50 is targeted to the nucleus and plasma membrane. GUS staining assays revealed broad expression across multiple tissues, including leaves, inflorescences, roots, and stems. Furthermore, hormonal treatment experiments showed that GhAGP50 expression is modulated by MeJA, ABA, IAA, and GA, suggesting its involvement in hormone-mediated developmental processes. This study presents a comprehensive genome-wide analysis of the AGP gene family in cotton, providing new insights into their structural diversity and functional significance. The identification and characterization of GhAGP50 highlight its potential role in epidermal hair formation and hormonal regulation, contributing to a deeper understanding of AGP functions in cotton development. These findings offer a valuable genetic resource for future research aimed at improving cotton growth and fiber quality through targeted genetic manipulation.

Allele and transcriptome mining in Gossypium hirsutum reveals variation in candidate genes at genetic loci affecting cotton fiber quality and textile flammability

BACKGROUND

Breeding valuable traits in crop plants requires identifying diverse alleles in the germplasm that are likely to affect desirable characteristics. The genetic diversity of historic cultivars of cotton is a reservoir of potentially important genes for crop improvement and genetic research. Diversity in the characteristics of harvested cotton fibers affects their suitability for end-use applications. Candidate loci and genes have been identified that affect the length, strength, and maturity of cotton fibers which affect the quality and value of the yarn, thread and textile. Natural genetic mechanisms in the plant may also affect the flammability of the produced textiles.

RESULTS

Here we show that a combination of allele mining and transcriptome analysis can identify candidate genes for cotton fiber traits including strength and perhaps flammability. We found novel DNA variants in fiber-expressed gene families in 132 newly sequenced cotton varieties and identified genes with genotype-specific RNA expression.

CONCLUSIONS

Among these, we identified novel variation in DNA sequence and RNA expression in genes at major QTL qD04-ELO-WLIM (JGI-Gohir.D04G160000), qA13-MIC (Gohir.A13G157500), qA07-STR (Gohir.A07G191600), supported the candidacy of qD11-UHML-KRP6 (Gohir.D11G197900) and qD13-STR (Gohir.D13G17450), and identified an additional A03-WLIM transcription factor gene (Gohir.A03G182100) and several RNA expression variant candidates of potential flammability genes that may be useful for plant biologists and cotton breeders. Candidate genes for traits like flame resistance that are likely due to the combination of many small effect QTL can benefit from this multi-mining approach. We provide an annotated variant call format (vcf) file with variations at 24,996 loci that are predicted to affect 10,418 cotton fiber genes in the historic breeding germplasm.

A high-resolution model of gene expression during Gossypium hirsutum (cotton) fiber development

BACKGROUND

Cotton fiber development relies on complex and intricate biological processes to transform newly differentiated fiber initials into the mature, extravagantly elongated cellulosic cells that are the foundation of this economically important cash crop. Here we extend previous research into cotton fiber development by employing controlled conditions to minimize variability and utilizing time-series sampling and analyses to capture daily transcriptomic changes from early elongation through the early stages of secondary wall synthesis (6 to 24 days post anthesis; DPA).

RESULTS

A majority of genes are expressed in fiber, largely partitioned into two major coexpression modules that represent genes whose expression generally increases or decreases during development. Differential gene expression reveals a massive transcriptomic shift between 16 and 17 DPA, corresponding to the onset of the transition phase that leads to secondary wall synthesis. Subtle gene expression changes are captured by the daily sampling, which are discussed in the context of fiber development. Coexpression and gene regulatory networks are constructed and associated with phenotypic aspects of fiber development, including turgor and cellulose production. Key genes are considered in the broader context of plant secondary wall synthesis, noting their known and putative roles in cotton fiber development.

CONCLUSIONS

The analyses presented here highlight the importance of fine-scale temporal sampling on understanding developmental processes and offer insight into genes and regulatory networks that may be important in conferring the unique fiber phenotype.

The GhCEWT1-GhCEWT2-GhCes4D/GhCOBL4D module orchestrates plant cell elongation and cell wall thickness

Cell elongation defines cell size and shape, whereas the cell wall supports and protects it. However, the mechanism regulating cell elongation and cell wall thickness remains unknown. Here, taking advantage of a model for both cell elongation and cell wall biogenesis, cotton fiber, we identified a basic-helix-loop-helix (bHLH) factor, GhCEWT1, that contributes to both fiber cell elongation and cell wall thickness. Loss of function of GhCEWT1 reduced the fiber length and cell wall thickness. GhCEWT1 induced transcription of GhCEWT2. We also identified two target genes of GhCEWT2, cellulose synthase 4D (GhCes4D) and COBRA-LIKE 4D (GhCOBL4D). GhCEWT2 enhanced the transcription of GhCes4D and GhCOBL4D. GhCOBL4D overexpression significantly enhanced cotton fiber cell length and cell wall thickness. Our results revealed a GhCEWT1-GhCEWT2-GhCes4D/GhCOBL4D cascade functioning in both fiber cell elongation and cell wall thickness. These findings provide a comprehensive understanding of plant cell elongation and cell wall formation, as well as a theoretical basis for boosting the biomass on Earth.

COBRA-LIKE 9 modulates cotton cell wall development via regulating cellulose deposition

Plant cell walls are complex and dynamic cellular structures critical for plant growth, development, physiology, and adaptation. Cellulose is one of the most important components of the cell wall. However, how cellulose microfibrils deposit and assemble into crystalline cellulose remains elusive. The COBRA-LIKE plant-specific protein family plays a vital role in modulating the deposition and orientation of cellulose microfibril in plant cell walls. Here, we investigate the role of GhCOBL9 in cotton (Gossypium hirsutum) fiber development, an ideal model for studying cell elongation and cell wall thickening. The expression period of GhCOBL9 is consistent with the thickening stage of the secondary wall of cotton fibers. Overexpression of GhCOBL9 results in increased cellulose content in the cell wall and produces shorter, thicker, and stronger fibers, while RNA interference (RNAi)-mediated downregulation of GhCOBL9 leads to the opposite phenotypes, indicating its crucial role in cell wall development. Subcellular localization and binding activity assays reveal that GhCOBL9 targets the cell wall and binds to crystalline cellulose with high affinity. Transcriptomic analysis of GhCOBL9 transgenic lines uncovers expression alterations in genes related to cellulose and monosaccharide biosynthesis. Furthermore, we identify a fasciclin-like arabinogalactan protein 9 (GhFLA9) as an interacting partner of GhCOBL9 to modulate cell wall development. Additionally, the R2R3-MYB transcription factor GhMYB46-5 activates GhCOBL9 expression by binding to the MYB46-responsive cis-regulatory element in the GhCOBL9 promoter. These findings broaden our knowledge of COBL function in modulating plant cell wall development.

Developmentally controlled subcellular remodeling and VND-initiated vacuole-executed PCD module shape xylem-like cells in peat moss

Peat moss (Sphagnum) is a non-vascular higher plant with unique xylem-like hyaline (H) cells that are accompanied by photosynthetic chlorophyllous cells. These cellular structures play crucial roles in water storage and carbon sequestration. However, it is largely unknown how peat moss develops the H cells. This study systematically explored the Sphagnum Developmental Cell Atlas and Lineage and classified leaf cell development into two lineages with six stages (S0-S5) based on changes in key cellular traits, including the formation of spiral secondary cell walls (S4) and the presence of water pores (S5). Cell lineage-specific subcellular remodeling was transcriptionally regulated during leaf development, and vacuole-mediated clearance of organelles and cell death led to mature dead H cells. Interestingly, expression of land plant conserved Vascular-related NAC Domain (VND) genes correlated with H cell formation. Overall, these results suggest that the origination of xylem-like H cells is related to VND, likely through the neofunctionalization of vacuole-mediated cell death to attempt xylem formation in peat moss, suggesting potential uncoupling of xylem and phloem cell origins. This study positions peat moss as a potential model organism for studying integrative evolutionary cell biology.

Cotton HD-Zip I transcription factor GhHB4-like regulates the plant response to salt stress

Soil salinity is a major environmental constraint to plant production. The homeodomain-leucine zipper I (HD-Zip I) transcription factors play a crucial role in growth, development and defence responses of plants. However, the function and underlying mechanism of HD-Zip I in cotton remain unexplored. This study investigated the role of GhHB4-like, a cotton HD-Zip I gene, in plant tolerance to salt stress. Ectopic expression of GhHB4-like gene enhanced, while its silencing impaired the salt tolerance in Arabidopsis. Y1H and effector-reporter assays revealed that GhHB4-like activated the expression of GhNAC007, which is essential for salt resistance. Knock-down of GhNAC007 also impaired salt resistance of cotton plants. In addition, GhHB4-like-GhNAC007 might have positively regulated the expression of GhMYB96 and ABA signalling-related genes, thereby leading to enhanced salt resistance. Interestingly, deleting motifs 3 and 5 near the 3'-end of GhHB4-like significantly enhanced GhNAC007 activation, indicating that both motifs acted as transcriptional activation inhibitory domains. The results suggest that GhHB4-like-GhNAC007 regulated plant response to salt stress, potentially by modulating GhMYB96 and ABA signalling-related genes.

Expression profile analysis of cotton fiber secondary cell wall thickening stage

To determine the genes associated with the fiber strength trait in cotton, three different cotton cultivars were selected: Sea Island cotton (Xinhai 32, with hyper-long fibers labeled as HL), and upland cotton (17-24, with long fibers labeled as L, and 62-33, with short fibers labeled as S). These cultivars were chosen to assess fiber samples with varying qualities. RNA-seq technology was used to analyze the expression profiles of cotton fibers at the secondary cell wall (SCW) thickening stage (20, 25, and 30 days post-anthesis (DPA)). The results showed that a large number of differentially expressed genes (DEGs) were obtained from the three assessed cotton cultivars at different stages of SCW development. For instance, at 20 DPA, Sea Island cotton (HL) had 6,215 and 5,364 DEGs compared to upland cotton 17-24 (L) and 62-33 (S), respectively. Meanwhile, there were 1,236 DEGs between two upland cotton cultivars, 17-24 (L) and 62-33 (S). Gene Ontology (GO) term enrichment identified 42 functions, including 20 biological processes, 11 cellular components, and 11 molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified several pathways involved in SCW synthesis and thickening, such as glycolysis/gluconeogenesis, galactose metabolism, propanoate metabolism, biosynthesis of unsaturated fatty acids pathway, valine, leucine and isoleucine degradation, fatty acid elongation pathways, and plant hormone signal transduction. Through the identification of shared DEGs, 46 DEGs were found to exhibit considerable expressional differences at different fiber stages from the three cotton cultivars. These shared DEGs have functions including REDOX enzymes, binding proteins, hydrolases (such as GDSL thioesterase), transferases, metalloproteins (cytochromatin-like genes), kinases, carbohydrates, and transcription factors (MYB and WRKY). Therefore, RT-qPCR was performed to verify the expression levels of nine of the 46 identified DEGs, an approach which demonstrated the reliability of RNA-seq data. Our results provided valuable molecular resources for clarifying the cell biology of SCW biosynthesis during fiber development in cotton.

Comparative Transcriptomic Analysis of Gossypium hirsutum Fiber Development in Mutant Materials (xin w 139) Provides New Insights into Cotton Fiber Development

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.

Transcriptome, Ectopic Expression and Genetic Population Analysis Identify Candidate Genes for Fiber Quality Improvement in Cotton

Comparative transcriptome analysis of fiber tissues between Gossypium barbadense and Gossypium hirsutum could reveal the molecular mechanisms underlying high-quality fiber formation and identify candidate genes for fiber quality improvement. In this study, 759 genes were found to be strongly upregulated at the elongation stage in G. barbadense, which showed four distinct expression patterns (I-IV). Among them, the 346 genes of group IV stood out in terms of the potential to promote fiber elongation, in which we finally identified 42 elongation-related candidate genes by comparative transcriptome analysis between G. barbadense and G. hirsutum. Subsequently, we overexpressed GbAAR3 and GbTWS1, two of the 42 candidate genes, in Arabidopsis plants and validated their roles in promoting cell elongation. At the secondary cell wall (SCW) biosynthesis stage, 2275 genes were upregulated and exhibited five different expression profiles (I-V) in G. barbadense. We highlighted the critical roles of the 647 genes of group IV in SCW biosynthesis and further picked out 48 SCW biosynthesis-related candidate genes by comparative transcriptome analysis. SNP molecular markers were then successfully developed to distinguish the SCW biosynthesis-related candidate genes from their G. hirsutum orthologs, and the genotyping and phenotyping of a BC3F5 population proved their potential in improving fiber strength and micronaire. Our results contribute to the better understanding of the fiber quality differences between G. barbadense and G. hirsutum and provide novel alternative genes for fiber quality improvement.