Comparative phosphoproteomic analysis of BR-defective mutant reveals a key role of GhSK13 in regulating cotton fiber development

Cotton metabolism regulatory network: Unraveling key genes and pathways in fiber development and growth regulation

Cotton (Gossypium hirsutum L.) is one of the world's most important commercial crops. However, the dynamics of metabolite abundance and potential regulatory networks throughout its life cycle remain poorly understood. In this study, we developed a cotton metabolism regulatory network (CMRN) that spans various developmental stages and encompasses 2138 metabolites and 90 309 expressed genesin upland cotton. By integrating high-resolution spatiotemporal metabolome and transcriptome data, we identified 1958 differentially accumulated metabolites and 13 597 co-expressed differentially expressed genes between the dwarf mutant pagoda1 and its wild-type counterpart Zhongmiansuo 24. These metabolites and genes were categorized into seven clusters based on tissue-specific accumulation patterns and gene expression profiles across different developmental stages. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed significant differential enrichment in the fatty acid elongation pathway, particularly in fibers. The differential involvement of genes and metabolites in very-long-chain fatty acid (VLCFA) synthesis led to the identification of GhKCS1b_Dt as a key gene. Overexpression of GhKCS1b_Dt significantly promoted fiber elongation, while its silencing markedly inhibited cotton fiber growth, affirming its positive regulatory role in fiber elongation. This dataset provides a valuable resource for further research into metabolic pathways and gene regulatory networks, offering novel insights for advancing cotton breeding strategies.

Brassinosteroids in cotton: orchestrating fiber development

Cotton cultivation spans over 30 million hectares across 85 countries and regions, with more than half participating in the global cotton textile trade. The elongated cotton fiber cell is an ideal model for studying cell elongation and understanding plant growth and development. Brassinosteroids (BRs), recognized for their role in cell elongation, offer the potential for improving cotton fiber quality and yield. Despite extensive research highlighting BR's positive impact on fiber development, a comprehensive review on this topic has been lacking. This review addresses this gap, providing a detailed analysis of the latest advancements in BR signaling and its effects on cotton fiber development. We explore the complex network of BR biosynthesis components, signaling molecules, and regulators, including crosstalk with other pathways and transcriptional control mechanisms. Additionally, we propose molecular strategies and highlight key genetic elements for optimizing BR-related genes to enhance fiber quality and yield. The review emphasizes the importance of BR homeostasis and the hormonal landscape during cotton fiber development, offering insights into targeted manipulation opportunities and challenges. This consolidation offers a comprehensive understanding of BR's multifaceted roles in fiber development, outlining a strategic approach for BR optimization in cotton fiber quality and yield.

Cotton BOP1 mediates SUMOylation of GhBES1 to regulate fibre development and plant architecture

The Arabidopsis BLADE-ON-PETIOLE (BOP) genes are primarily known for their roles in regulating leaf and floral patterning. However, the broader functions of BOPs in regulating plant traits remain largely unexplored. In this study, we investigated the role of the Gossypium hirsutum BOP1 gene in the regulation of fibre length and plant height through the brassinosteroid (BR) signalling pathway. Transgenic cotton plants overexpressing GhBOP1 display shorter fibre lengths and reduced plant height compared to the wild type. Conversely, GhBOP1 knockdown led to increased plant height and longer fibre, indicating a connection with phenotypes influenced by the BR pathway. Our genetic evidence supports the notion that GhBOP1 regulates fibre length and plant height in a GhBES1-dependent manner, with GhBES1 being a major transcription factor in the BR signalling pathway. Yeast two-hybrid, luciferase complementation assay and pull-down assay results demonstrated a direct interaction between GhBOP1 and GhSUMO1, potentially forming protein complexes with GhBES1. In vitro and in vivo SUMOylation analyses revealed that GhBOP1 functions in an E3 ligase-like manner to mediate GhBES1 SUMOylation and subsequent degradation. Therefore, our study not only uncovers a novel mechanism of GhBES1 SUMOylation but also provides significant insights into how GhBOP1 regulates fibre length and plant height by controlling GhBES1 accumulation.

Co-Expression Network Analysis and Introgressive Gene Identification for Fiber Length and Strength Reveal Transcriptional Differences in 15 Cotton Chromosome Substitution Segment Lines and Their Upland and Sea Island Parents

Fiber length (FL) and strength (FS) are the core indicators for evaluating cotton fiber quality. The corresponding stages of fiber elongation and secondary wall thickening are of great significance in determining FL and FS formation, respectively. QTL mapping and high-throughput sequencing technology have been applied to dissect the molecular mechanism of fiber development. In this study, 15 cotton chromosome segment substitution lines (CSSLs) with significant differences in FL and FS, together with their recurrent parental Gossypium hirsutum line CCRI45 and donor parent G. barbadense line Hai1, were chosen to conduct RNA-seq on developing fiber samples at 10 days post anthesis (DPA) and 20 DPA. Differentially expressed genes (DEGs) were obtained via pairwise comparisons among all 24 samples (each one with three biological repeats). A total of 969 DEGs related to FL-high, 1285 DEGs to FS-high, and 997 DEGs to FQ-high were identified. The functional enrichment analyses of them indicated that the GO terms of cell wall structure and ROS, carbohydrate, and phenylpropanoid metabolism were significantly enriched, while the GO terms of glucose and polysaccharide biosynthesis, and brassinosteroid and glycosylphosphatidylinositol metabolism could make great contributions to FL and FS formation, respectively. Weighted gene co-expressed network analyses (WGCNA) were separately conducted for analyzing FL and FS traits, and their corresponding hub DEGs were screened in significantly correlated expression modules, such as EXPA8, XTH, and HMA in the fiber elongation and WRKY, TDT, and RAC-like 2 during secondary wall thickening. An integrated analysis of these hub DEGs with previous QTL identification results successfully identified a total of 33 candidate introgressive DEGs with non-synonymous mutations between the Gh and Gb species. A common DEG encoding receptor-like protein kinase 1 was reported to likely participate in fiber secondary cell thickening regulation by brassionsteroid signaling. Such valuable information was conducive to enlightening the developing mechanism of cotton fiber and also provided an abundant gene pool for further molecular breeding.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

A conserved brassinosteroid-mediated BES1-CERP-EXPA3 signaling cascade controls plant cell elongation

Continuous plant growth is achieved by cell division and cell elongation. Brassinosteroids control cell elongation and differentiation throughout plant life. However, signaling cascades underlying BR-mediated cell elongation are unknown. In this study, we introduce cotton fiber, one of the most representative single-celled tissues, to decipher cell-specific BR signaling. We find that gain of function of GhBES1, a key transcriptional activator in BR signaling, enhances fiber elongation. The chromatin immunoprecipitation sequencing analysis identifies a cell-elongation-related protein, GhCERP, whose transcription is directly activated by GhBES1. GhCERP, a downstream target of GhBES1, transmits the GhBES1-mediated BR signaling to its target gene, GhEXPA3-1. Ultimately, GhEXPA3-1 promotes fiber cell elongation. In addition, inter-species functional analysis of the BR-mediated BES1-CERP-EXPA3 signaling cascade also promotes Arabidopsis root and hypocotyl growth. We propose that the BES1-CERP-EXPA3 module may be a broad-spectrum pathway that is universally exploited by diverse plant species to regulate BR-promoted cell elongation.

Regulatory network of GSK3-like kinases and their role in plant stress response

Glycogen synthase kinase 3 (GSK3) family members are evolutionally conserved Ser/Thr protein kinases in mammals and plants. In plants, the GSK3s function as signaling hubs to integrate the perception and transduction of diverse signals required for plant development. Despite their role in the regulation of plant growth and development, emerging research has shed light on their multilayer function in plant stress responses. Here we review recent advances in the regulatory network of GSK3s and the involvement of GSK3s in plant adaptation to various abiotic and biotic stresses. We also discuss the molecular mechanisms underlying how plants cope with environmental stresses through GSK3s-hormones crosstalk, a pivotal biochemical pathway in plant stress responses. We believe that our overview of the versatile physiological functions of GSK3s and underlined molecular mechanism of GSK3s in plant stress response will not only opens further research on this important topic but also provide opportunities for developing stress-resilient crops through the use of genetic engineering technology.

Metabolomics-centered mining of plant metabolic diversity and function: Past decade and future perspectives

Plants are natural experts in organic synthesis, being able to generate large numbers of specific metabolites with widely varying structures that help them adapt to variable survival challenges. Metabolomics is a research discipline that integrates the capabilities of several types of research including analytical chemistry, statistics, and biochemistry. Its ongoing development provides strategies for gaining a systematic understanding of quantitative changes in the levels of metabolites. Metabolomics is usually performed by targeting either a specific cell, a specific tissue, or the entire organism. Considerable advances in science and technology over the last three decades have propelled us into the era of multi-omics, in which metabolomics, despite at an earlier developmental stage than genomics, transcriptomics, and proteomics, offers the distinct advantage of studying the cellular entities that have the greatest influence on end phenotype. Here, we summarize the state of the art of metabolite detection and identification, and illustrate these techniques with four case study applications: (i) comparing metabolite composition within and between species, (ii) assessing spatio-temporal metabolic changes during plant development, (iii) mining characteristic metabolites of plants in different ecological environments and upon exposure to various stresses, and (iv) assessing the performance of metabolomics as a means of functional gene identification , metabolic pathway elucidation, and metabolomics-assisted breeding through analyzing plant populations with diverse genetic variations. In addition, we highlight the prominent contributions of joint analyses of plant metabolomics and other omics datasets, including those from genomics, transcriptomics, proteomics, epigenomics, phenomics, microbiomes, and ion-omics studies. Finally, we discuss future directions and challenges exploiting metabolomics-centered approaches in understanding plant metabolic diversity.

Identification of Candidate Genes for Lint Percentage and Fiber Quality Through QTL Mapping and Transcriptome Analysis in an Allotetraploid Interspecific Cotton CSSLs Population

Upland cotton (Gossypium hirsutum) has long been an important fiber crop, but the narrow genetic diversity of modern G. hirsutum limits the potential for simultaneous improvement of yield and fiber quality. It is an effective approach to broaden the genetic base of G. hirsutum through introgression of novel alleles from G. barbadense with excellent fiber quality. In the present study, an interspecific chromosome segment substitution lines (CSSLs) population was established using G. barbadense cultivar Pima S-7 as the donor parent and G. hirsutum cultivar CCRI35 as the recipient parent. A total of 105 quantitative trait loci (QTL), including 85 QTL for fiber quality and 20 QTL for lint percentage (LP), were identified based on phenotypic data collected from four environments. Among these QTL, 25 stable QTL were detected in two or more environments, including four for LP, eleven for fiber length (FL), three for fiber strength (FS), six for fiber micronaire (FM), and one for fiber elongation (FE). Eleven QTL clusters were observed on nine chromosomes, of which seven QTL clusters harbored stable QTL. Moreover, eleven major QTL for fiber quality were verified through analysis of introgressed segments of the eight superior lines with the best comprehensive phenotypes. A total of 586 putative candidate genes were identified for 25 stable QTL associated with lint percentage and fiber quality through transcriptome analysis. Furthermore, three candidate genes for FL, GH_A08G1681 (GhSCPL40), GH_A12G2328 (GhPBL19), and GH_D02G0370 (GhHSP22.7), and one candidate gene for FM, GH_D05G1346 (GhAPG), were identified through RNA-Seq and qRT-PCR analysis. These results lay the foundation for understanding the molecular regulatory mechanism of fiber development and provide valuable information for marker-assisted selection (MAS) in cotton breeding.

Molecular studies of cellulose synthase supercomplex from cotton fiber reveal its unique biochemical properties

Cotton fiber is a highly elongated and thickened single cell that produces large quantities of cellulose, which is synthesized and assembled into cell wall microfibrils by the cellulose synthase complex (CSC). In this study, we report that in cotton (Gossypium hirsutum) fibers harvested during secondary cell wall (SCW) synthesis, GhCesA 4, 7, and 8 assembled into heteromers in a previously uncharacterized 36-mer-like cellulose synthase supercomplex (CSS). This super CSC was observed in samples prepared using cotton fiber cells harvested during the SCW synthesis period but not from cotton stem tissue or any samples obtained from Arabidopsis. Knock-out of any of GhCesA 4, 7, and 8 resulted in the disappearance of the CSS and the production of fiber cells with no SCW thickening. Cotton fiber CSS showed significantly higher enzyme activity than samples prepared from knock-out cotton lines. We found that the microfibrils from the SCW of wild-type cotton fibers may contain 72 glucan chains in a bundle, unlike other plant materials studied. GhCesA4, 7, and 8 restored both the dwarf and reduced vascular bundle phenotypes of their orthologous Arabidopsis mutants, potentially by reforming the CSC hexamers. Genetic complementation was not observed when non-orthologous CesA genes were used, indicating that each of the three subunits is indispensable for CSC formation and for full cellulose synthase function. Characterization of cotton CSS will increase our understanding of the regulation of SCW biosynthesis.

Genome-wide identification and expression analysis of GL2-interacting-repressor (GIR) genes during cotton fiber and fuzz development

GL2-interacting-repressor (GIR) family members may contribute to fiber/fuzz formation via a newly discovered unique pathway in Gossypium arboreum. There are similarities between cotton fiber development and the formation of trichomes and root hairs. The GL2-interacting-repressors (GIRs) are crucial regulators of root hair and trichome formation. The GaFzl gene, annotated as GaGIR1, is negatively associated with trichome development and fuzz initiation. However, there is relatively little available information regarding the other GIR genes in cotton, especially regarding their effects on cotton fiber development. In this study, 21 GIR family genes were identified in the diploid cotton species Gossypium arboreum; these genes were divided into three groups. The GIR genes were characterized in terms of their phylogenetic relationships, structures, chromosomal distribution and evolutionary dynamics. These GIR genes were revealed to be unequally distributed on 12 chromosomes in the diploid cotton genome, with no GIR gene detected on Ga06. The cis-acting elements in the promoter regions were predicted to be responsive to light, phytohormones, defense activities and stress. The transcriptomic data and qRT-PCR results revealed that most GIR genes were not differentially expressed between the wild-type control and the fuzzless mutant line. Moreover, 14 of 21 family genes were expressed at high levels, indicating these genes may play important roles during fiber development and fuzz formation. Furthermore, Ga01G0231 was predominantly expressed in root samples, suggestive of a role in root hair formation rather than in fuzz initiation and development. The results of this study have enhanced our understanding of the GIR genes and their potential utility for improving cotton fiber through breeding.

GSK3s: nodes of multilayer regulation of plant development and stress responses

Glycogen synthase kinase 3 (GSK3) family members are highly conserved serine/threonine protein kinases in eukaryotes. Unlike animals, plants have evolved with multiple homologs of GSK3s involved in a diverse array of biological processes. Emerging evidence suggests that GSK3s act as signaling hubs for integrating perception and transduction of diverse signals required for plant development and responses to abiotic and biotic cues. Here we review recent advances in understanding the molecular interactions between GSK3s and an expanding spectrum of their upstream regulators and downstream substrates in plants. We further discuss how GSK3s act as key signaling nodes of multilayer regulation of plant development and stress response through either being regulated at the post-translational level or regulating their substrates via phosphorylation.

Recent Advances and Future Perspectives in Cotton Research

Cotton is not only the world's most important natural fiber crop, but it is also an ideal system in which to study genome evolution, polyploidization, and cell elongation. With the assembly of five different cotton genomes, a cotton-specific whole-genome duplication with an allopolyploidization process that combined the A- and D-genomes became evident. All existing A-genomes seemed to originate from the A₀-genome as a common ancestor, and several transposable element bursts contributed to A-genome size expansion and speciation. The ethylene production pathway is shown to regulate fiber elongation. A tip-biased diffuse growth mode and several regulatory mechanisms, including plant hormones, transcription factors, and epigenetic modifications, are involved in fiber development. Finally, we describe the involvement of the gossypol biosynthetic pathway in the manipulation of herbivorous insects, the role of GoPGF in gland formation, and host-induced gene silencing for pest and disease control. These new genes, modules, and pathways will accelerate the genetic improvement of cotton.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes Contributing to Fuzz Development in Gossypium arboreum

Fuzzless mutants are ideal materials to decipher the regulatory network and mechanism underlying fuzz initiation and formation. In this study, we utilized two Gossypium arboreum accessions differing in fuzz characteristics to explore expression pattern differences and discriminate genes involved in fuzz development using RNA sequencing. Gene ontology (GO) analysis was conducted and found that DEGs were mainly enriched in the regulation of transcription, metabolic processes and oxidation–reduction-related processes. Weighted gene co-expression network analysis discerned the MEmagenta module highly associated with a fuzz/fuzzless trait, which included a total of 50 hub genes differentially expressed between two materials. GaFZ, which negatively regulates trichome and fuzz formation, was found involved in MEmagenta cluster1. In addition, twenty-eight hub genes in MEmagenta cluster1 were significantly up-regulated and expressed in fuzzless mutant DPL972. It is noteworthy that Ga04G1219 and Ga04G1240, which, respectively, encode Fasciclin-like arabinogalactan protein 18(FLA18) and transport protein, showed remarkable differences of expression level and implied that they may be involved in protein glycosylation to regulate fuzz formation and development. This module and hub genes identified in this study will provide new insights on fiber and fuzz formation and be useful for the molecular design breeding of cotton genetic improvement.