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

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.