Cotton fiber growth in planta and in vitro. Models for plant cell elongation and cell wall biogenesis

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Genome resources for three modern cotton lines guide future breeding efforts

Cotton (Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars ('UGA230', 'UA48' and 'CSX8308') and updating the 'TM-1' cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from 'Pima' cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

Co-localization and analysis of miR477b with fiber length quantitative trait loci in cotton

Cotton is an important cash crop for the textile industry. However, the understanding of natural genetic variation of fiber elongation in relation to miRNA is lacking. A miRNA gene (miR477b) was found to co-localize with a previously mapped fiber length (FL) quantitative trait locus (QTL). The miR477b was differentially expressed during fiber elongation between two backcross inbred lines (BILs) differing in FL and its precursor sequences. Bioinformatics and qRT-PCR analysis were further used to analyse the miRNA genes, which could produce mature miR477b. Cotton plants with virus-induced gene silencing (VIGS) constructs to over-express the allele of miR477b from the BIL with longer fibers had significantly longer fibers as compared with negative control plants, while the VIGS plants with suppressed miRNA expression had significantly shorter fibers. The expression level of the target gene (DELLA) and related genes (RDL1 and EXPA1 for DELLA through HOX3 protein) in the two BILs and/or the VIGS plants were generally congruent, as expected. This report represents one of the first comprehensive studies to integrate QTL linkage mapping and physical mapping of small RNAs with both small and mRNA transcriptome analysis, followed by VIGS, to identify candidate small RNA genes affecting the natural variation of fiber elongation in cotton.

GhMYB52 Like: A Key Factor That Enhances Lint Yield by Negatively Regulating the Lignin Biosynthesis Pathway in Fibers of Upland Cotton (Gossypium hirsutum L.)

In the context of sustainable agriculture and biomaterial development, understanding and enhancing plant secondary cell wall formation are crucial for improving crop fiber quality and biomass conversion efficiency. This is especially critical for economically important crops like upland cotton (Gossypium hirsutum L.), for which fiber quality and its processing properties are essential. Through comprehensive genome-wide screening and analysis of expression patterns, we identified a particularly high expression of an R2R3 MYB transcription factor, GhMYB52 Like, in the development of the secondary cell wall in cotton fiber cells. Utilizing gene-editing technology to generate a loss-of-function mutant to clarify the role of GhMYB52 Like, we revealed that GhMYB52 Like does not directly contribute to cellulose synthesis in cotton fibers but instead represses a subset of lignin biosynthesis genes, establishing it as a lignin biosynthesis inhibitor. Concurrently, a substantial decrease in the lint index, a critical measure of cotton yield, was noted in parallel with an elevation in lignin levels. This study not only deepens our understanding of the molecular mechanisms underlying cotton fiber development but also offers new perspectives for the molecular improvement of other economically important crops and the enhancement of biomass energy utilization.

Fourier Transform Infrared (FT-IR) Spectroscopy and Simple Algorithm Analysis for Rapid and Non-Destructive Assessment of Cotton Fiber Maturity and Crystallinity for Plant Mapping

Information on boll distribution within a cotton plant is critical to evaluate the adaptation and response of cotton plants to environmental and biotic stress in cotton production. Cotton researchers have applied available conventional fiber measurements, such as the high volume instrument (HVI) and advanced fiber information system (AFIS), to map the location and the timing of boll development and distribution within plants and further to determine within-plant variability of cotton fiber properties. Both HVI and AFIS require numerous cotton bolls combined for the measurement. As an alternative approach, attenuated total reflection Fourier transform infrared (ATR FT-IR) spectroscopy was proposed to measure fiber maturity (MIR) and crystallinity (CIIR) of a sample as little as 0.5 mg lint. Extending fiber maturity and crystallinity measurement into a single boll for node-by-node mapping, FT-IR method might be advantageous due to less sampling amount compared with HVI and AFIS methods. Results showed that FT-IR technique enabled the evaluation of fiber MIR and CIIR at a boll level, which resulted in average MIR and CIIR values highly correlated with HVI micronaire (MIC) and AFIS maturity ratio (M). Hence, FT-IR technique possesses a good potential for a rapid and non-destructive node-by-node mapping of cotton boll maturity and crystallinity distribution.

Phosphatidic acid interacts with an HD-ZIP transcription factor GhHOX4 to influence its function in fiber elongation of cotton (Gossypium hirsutum)

Upland cotton, the mainly cultivated cotton species in the world, provides over 90% of natural raw materials (fibers) for the textile industry. The development of cotton fibers that are unicellular and highly elongated trichomes on seeds is a delicate and complex process. However, the regulatory mechanism of fiber development is still largely unclear in detail. In this study, we report that a homeodomain-leucine zipper (HD-ZIP) IV transcription factor, GhHOX4, plays an important role in fiber elongation. Overexpression of GhHOX4 in cotton resulted in longer fibers, while GhHOX4-silenced transgenic cotton displayed a "shorter fiber" phenotype compared with wild type. GhHOX4 directly activates two target genes, GhEXLB1D and GhXTH2D, for promoting fiber elongation. On the other hand, phosphatidic acid (PA), which is associated with cell signaling and metabolism, interacts with GhHOX4 to hinder fiber elongation. The basic amino acids KR-R-R in START domain of GhHOX4 protein are essential for its binding to PA that could alter the nuclear localization of GhHOX4 protein, thereby suppressing the transcriptional regulation of GhHOX4 to downstream genes in the transition from fiber elongation to secondary cell wall (SCW) thickening during fiber development. Thus, our data revealed that GhHOX4 positively regulates fiber elongation, while PA may function in the phase transition from fiber elongation to SCW formation by negatively modulating GhHOX4 in cotton.

Synergistic interplay of redox homeostasis and polysaccharide synthesis promotes cotton fiber elongation

Cell polarity is the foundation of cell development and tissue morphogenesis. The investigation of polarized growth provides opportunities to gain profound insights into morphogenesis and tissue functionality in organisms. Currently, there are still many mysteries surrounding the mechanisms that regulate polarized cell growth. Cotton fiber cells serve as an excellent model for studying polarized growth, and provide important clues for unraveling the molecular mechanisms, signaling pathways, and regulatory networks of polarized growth. In this study, we characterized two functional genes, GhMDHAR1AT/DT and GhDHAR2AT/DT with predominant expression during fiber elongation. Loss of function of both genes contributed to a significant increase in fiber length. Transcriptomic data revealed up-regulated expression of antioxidant genes in CRISPR mutant lines, along with delayed expression of secondary wall-related genes and temporally prolonged expression of primary wall-related genes. Experimental evidence demonstrated that the increase in GSH content and glutathione peroxidase (GPX) enzyme activity led to enhanced total antioxidant capacity (T-AOC), resulting in reduced H2O2 levels, which contributed to the extension of fiber elongation stage in CRISPR mutant lines. Moreover, the increased polysaccharide synthesis in CRISPR mutant lines was found to provide an abundant supply of raw materials for fiber cell wall elongation, suggesting that synergistic interplay between redox homeostasis and polysaccharide synthesis in fiber cells may facilitate cell wall remodeling and fiber elongation. This study provides valuable insights for deciphering the mechanisms of cell polarized growth and improving cotton fiber quality.

Cotton BLH1 and KNOX6 antagonistically modulate fiber elongation via regulation of linolenic acid biosynthesis

BEL1-LIKE HOMEODOMAIN (BLH) proteins are known to function in various plant developmental processes. However, the role of BLHs in regulating plant cell elongation is still unknown. Here, we identify a BLH gene, GhBLH1, that positively regulates fiber cell elongation. Combined transcriptomic and biochemical analyses reveal that GhBLH1 enhances linolenic acid accumulation to promote cotton fiber cell elongation by activating the transcription of GhFAD7A-1 via binding of the POX domain of GhBLH1 to the TGGA cis-element in the GhFAD7A-1 promoter. Knockout of GhFAD7A-1 in cotton significantly reduces fiber length, whereas overexpression of GhFAD7A-1 results in longer fibers. The K2 domain of GhKNOX6 directly interacts with the POX domain of GhBLH1 to form a functional heterodimer, which interferes with the transcriptional activation of GhFAD7A-1 via the POX domain of GhBLH1. Overexpression of GhKNOX6 leads to a significant reduction in cotton fiber length, whereas knockout of GhKNOX6 results in longer cotton fibers. An examination of the hybrid progeny of GhBLH1 and GhKNOX6 transgenic cotton lines provides evidence that GhKNOX6 negatively regulates GhBLH1-mediated cotton fiber elongation. Our results show that the interplay between GhBLH1 and GhKNOX6 modulates regulation of linolenic acid synthesis and thus contributes to plant cell elongation.

Genome-wide identification of the key kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L.)

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (G. hirsutum L.). Results showed that 159 kinesin genes, including 15 genes of the kinesin-13 gene subfamily, were identified in upland cotton; of which 157 kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 kinesin genes in upland cotton, including 10 kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 kinesin genes were significantly associated with three fiber traits, among which a kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

Dissecting the major genetic components underlying cotton lint development

Numerous genetic loci and several functionally characterized genes have been linked to determination of lint percentage (lint%), one of the most important cotton yield components, but we still know little about the major genetic components underlying lint%. Here, we first linked the genetic loci containing MYB25-like_At and HD1_At to the fiberless seed trait of 'SL1-7-1' and found that MYB25-like_At and HD1_At were very lowly expressed in 'SL1-7-1' ovules during fiber initiation. We then dissected the genetic components involved in determination of lint% using segregating populations derived from crosses of fuzzless mutants and intermediate segregants with different lint%, which not only confirmed the HD1_At locus but identified the HD1_Dt locus as being the major genetic components contributing to fiber initiation and lint%. The segregating populations also allowed us to evaluate the relative contributions of MYB25-like_At, MYB25-like_Dt, HD1_At, and HD1_Dt to lint%. Haplotype analysis of an Upland cotton (Gossypium hirsutum) population with 723 accessions (including 81 fuzzless seed accessions) showed that lint% of the accessions with the LP allele (higher lint%) at MYB25-like_At, MYB25-like_Dt, or HD1_At was significantly higher than that with the lp allele (lower lint%). The lint% of the Upland cotton accessions with 3 or 4 LP alleles at MYB25-like and HD1 was significantly higher than that with 2 LP alleles. The results prompted us to propose a strategy for breeding high-yielding cotton varieties, i.e. pyramiding the LP alleles of MYB25-like and HD1 with new lint% LP alleles without negative impact on seed size and fiber quality.

Genome-wide identification of the key Kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L)

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of Kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of Kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key Kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (Gossypium hirsutum L.). Results showed that 159 Kinesin genes, including 15 genes of the Kinesin-13 gene subfamily, were identified in upland cotton; of which 157 Kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 Kinesin genes in upland cotton, including 10 Kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven Kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 Kinesin genes were significantly associated with three fiber traits, among which a Kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one Kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the Kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited Kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

EB1C forms dimer and interacts with protein phosphatase 2A (PP2A) to regulate fiber elongation in upland cotton (Gossypium hirsutum)

Cotton is the most economically important natural fiber crop grown in more than sixty-five countries of the world. Fiber length is the main factor affecting fiber quality, but the existing main varieties are short in length and cannot suit the higher demands of the textile industry. It is necessary to discover functional genes that enable fiber length improvement in cotton through molecular breeding. In this study, overexpression of GhEB1C in Arabidopsis thaliana significantly promotes trichomes, tap roots, and root hairs elongation. The molecular regulation of GhEB1C involves its interactions with itself and GhB'ETA, and the function of GhEB1C regulation mainly depends on the two cysteine residues located at the C-terminal. In particular, the function activity of GhEB1C protein triggered with the regulation of protein phosphatase 2A, while silencing of GhEB1C in cotton significantly influenced the fiber protrusions and elongation mechanisms., Further, influenced the expression of MYB-bHLH-WD40 complex, brassinosteroids, and jasmonic acid-related genes, which showed that transcriptional regulation of GhEB1C is indispensable for cotton fiber formation and elongation processes. Our study analyzed the brief molecular mechanism of GhEB1C regulation. Further elucidated that GhEB1C can be a potential target gene to improve cotton fiber length through transgenic breeding.

A cell wall-localized β-1,3-glucanase promotes fiber cell elongation and secondary cell wall deposition

β-1,3-glucanase functions in plant physiological and developmental processes. However, how β-1,3-glucanase participates in cell wall development remains largely unknown. Here, we answered this question by examining the role of GhGLU18, a β-1,3-glucanase, in cotton (Gossypium hirsutum) fibers, in which the content of β-1,3-glucan changes dynamically from 10% of the cell wall mass at the onset of secondary wall deposition to <1% at maturation. GhGLU18 was specifically expressed in cotton fiber with higher expression in late fiber elongation and secondary cell wall (SCW) synthesis stages. GhGLU18 largely localized to the cell wall and was able to hydrolyze β-1,3-glucan in vitro. Overexpression of GhGLU18 promoted polysaccharide accumulation, cell wall reconstruction, and cellulose synthesis, which led to increased fiber length and strength with thicker cell walls and shorter pitch of the fiber helix. However, GhGLU18-suppressed cotton resulted in opposite phenotypes. Additionally, GhGLU18 was directly activated by GhFSN1 (fiber SCW-related NAC1), a NAC transcription factor reported previously as the master regulator in SCW formation during fiber development. Our results demonstrate that cell wall-localized GhGLU18 promotes fiber elongation and SCW thickening by degrading callose and enhancing polysaccharide metabolism and cell wall synthesis.

Regulatory controls of duplicated gene expression during fiber development in allotetraploid cotton

Polyploidy complicates transcriptional regulation and increases phenotypic diversity in organisms. The dynamics of genetic regulation of gene expression between coresident subgenomes in polyploids remains to be understood. Here we document the genetic regulation of fiber development in allotetraploid cotton Gossypium hirsutum by sequencing 376 genomes and 2,215 time-series transcriptomes. We characterize 1,258 genes comprising 36 genetic modules that control staged fiber development and uncover genetic components governing their partitioned expression relative to subgenomic duplicated genes (homoeologs). Only about 30% of fiber quality-related homoeologs show phenotypically favorable allele aggregation in cultivars, highlighting the potential for subgenome additivity in fiber improvement. We envision a genome-enabled breeding strategy, with particular attention to 48 favorable alleles related to fiber phenotypes that have been subjected to purifying selection during domestication. Our work delineates the dynamics of gene regulation during fiber development and highlights the potential of subgenomic coordination underpinning phenotypes in polyploid plants.

LIPID TRANSFER PROTEIN4 regulates cotton ceramide content and activates fiber cell elongation

Cell elongation is a fundamental process for plant growth and development. Studies have shown lipid metabolism plays important role in cell elongation; however, the related functional mechanisms remain largely unknown. Here, we report that cotton (Gossypium hirsutum) LIPID TRANSFER PROTEIN4 (GhLTP4) promotes fiber cell elongation via elevating ceramides (Cers) content and activating auxin-responsive pathways. GhLTP4 was preferentially expressed in elongating fibers. Over-expression and down-regulation of GhLTP4 led to longer and shorter fiber cells, respectively. Cers were greatly enriched in GhLTP4-overexpressing lines and decreased dramatically in GhLTP4 down-regulating lines. Moreover, auxin content and transcript levels of indole-3-acetic acid (IAA)-responsive genes were significantly increased in GhLTP4-overexpressing cotton fibers. Exogenous application of Cers promoted fiber elongation, while NPA (N-1-naphthalic acid, a polar auxin transport inhibitor) counteracted the promoting effect, suggesting that IAA functions downstream of Cers in regulating fiber elongation. Furthermore, we identified a basic helix-loop-helix transcription factor, GhbHLH105, that binds to the E-box element in the GhLTP4 promoter region and promotes the expression of GhLTP4. Suppression of GhbHLH105 in cotton reduced the transcripts level of GhLTP4, resulting in smaller cotton bolls and decreased fiber length. These results provide insights into the complex interactions between lipids and auxin-signaling pathways to promote plant cell elongation.

Over-expression of GhACTIN1 under the control of GhSCFP promoter improves cotton fiber and yield

Actin dynamics is pivotal in controlling cotton fiber elongation and the onset of secondary wall biosynthesis. We report that overexpression of GhACTIN1 under fiber fiber-specific promoter, GhSCFP, improves cotton fiber length, strength, and micronaire value. However, the effect of transgene has a more positive effect on fiber strength and micronaire value than fiber length. F-actin quantification and cellulose contents measurement in transgenic developing cotton fiber during the elongation phase showed an increase of up to 8.7% and 4.7% respectively. Additionally, physiological factors such as water use efficiency showed no significant change in transgenic cotton lines, while stomatal conductance and photosynthetic rate were significantly increased. Moreover, agronomical data determined that lint percentage (GOT) and seed cotton yield also increased up to 4.6% and 29.5% respectively, in transgenic cotton lines compared to the control lines. Our data demonstrate that the GhACTIN1 gene is a strong candidate gene for cotton fiber and yield improvement.

A Simplified Microscopy Technique to Rapidly Characterize Individual Fiber Traits in Cotton

Recent advances in phenotyping techniques have substantially improved the ability to mitigate type-II errors typically associated with high variance in phenotyping data sets. In particular, the implementation of automated techniques such as the High-Volume Instrument (HVI) and the Advanced Fiber Information System (AFIS) have significantly enhanced the reproducibility and standardization of various fiber quality measurements in cotton. However, micronaire is not a direct measure of either maturity or fineness, lending to limitations. AFIS only provides a calculated form of fiber diameter, not a direct measure, justifying the need for a visual-based reference method. Obtaining direct measurements of individual fibers through cross-sectional analysis and electron microscopy is a widely accepted standard but is time-consuming and requires the use of hazardous chemicals and specialized equipment. In this study, we present a simplified fiber histology and image acquisition technique that is both rapid and reproducible. We also introduce an automated image analysis program that utilizes machine learning to differentiate good fibers from bad and to subsequently collect critical phenotypic measurements. These methods have the potential to improve the efficiency of cotton fiber phenotyping, allowing for greater precision in unravelling the genetic architecture of critical traits such as fiber diameter, shape, areas of the secondary cell wall/lumen, and others, ultimately leading to larger genetic gains in fiber quality and improvements in cotton.

GoSTR, a negative modulator of stem trichome formation in cotton

Trichomes, the outward projection of plant epidermal tissue, provide an effective defense against stress and insect pests. Although numerous genes have been identified to be involved in trichome development, the molecular mechanism for trichome cell fate determination is not well enunciated. Here, we reported GoSTR functions as a master repressor for stem trichome formation, which was isolated by map-based cloning based on a large F2 segregating population derived from a cross between TM-1 (pubescent stem) and J220 (smooth stem). Sequence alignment revealed a critical G-to-T point mutation in GoSTR's coding region that converted codon 2 from GCA (Alanine) to TCA (Serine). This mutation occurred between the majority of Gossypium hirsutum with pubescent stem (GG-haplotype) and G. barbadense with glabrous stem (TT-haplotype). Silencing of GoSTR in J220 and Hai7124 via virus-induced gene silencing resulted in the pubescent stems but no visible change in leaf trichomes, suggesting stem trichomes and leaf trichomes are genetically distinct. Yeast two-hybrid assay and luciferase complementation imaging assay showed GoSTR interacts with GoHD1 and GoHOX3, two key regulators of trichome development. Comparative transcriptomic analysis further indicated that many transcription factors such as GhMYB109, GhTTG1, and GhMYC1/GhDEL65 which function as positive regulators of trichomes were significantly upregulated in the stem from the GoSTR-silencing plant. Taken together, these results indicate that GoSTR functions as an essential negative modulator of stem trichomes and its transcripts will greatly repress trichome cell differentiation and growth. This study provided valuable insights for plant epidermal hair initiation and differentiation research.

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.

Integration of eQTL and machine learning to dissect causal genes with pleiotropic effects in genetic regulation networks of seed cotton yield

The dissection of a gene regulatory network (GRN) that complements the genome-wide association study (GWAS) locus and the crosstalk underlying multiple agronomical traits remains a major challenge. In this study, we generate 558 transcriptional profiles of lint-bearing ovules at one day post-anthesis from a selective core cotton germplasm, from which 12,207 expression quantitative trait loci (eQTLs) are identified. Sixty-six known phenotypic GWAS loci are colocalized with 1,090 eQTLs, forming 38 functional GRNs associated predominantly with seed yield. Of the eGenes, 34 exhibit pleiotropic effects. Combining the eQTLs within the seed yield GRNs significantly increases the portion of narrow-sense heritability. The extreme gradient boosting (XGBoost) machine learning approach is applied to predict seed cotton yield phenotypes on the basis of gene expression. Top-ranking eGenes (NF-YB3, FLA2, and GRDP1) derived with pleiotropic effects on yield traits are validated, along with their potential roles by correlation analysis, domestication selection analysis, and transgenic plants.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development

Introduction

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels.

Methods

Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes.

Results

Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers.

Discussion

Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

Role of Actin Dynamics and GhACTIN1 Gene in Cotton Fiber Development: A Prototypical Cell for Study

Cotton crop is considered valuable for its fiber and seed oil. Cotton fiber is a single-celled outgrowth from the ovule epidermis, and it is a very dynamic cell for study. It has four distinct but overlapping developmental stages: initiation, elongation, secondary cell wall synthesis, and maturation. Among the various qualitative characteristics of cotton fiber, the important ones are the cotton fiber staple length, tensile strength, micronaire values, and fiber maturity. Actin dynamics are known to play an important role in fiber elongation and maturation. The current review gives an insight into the cotton fiber developmental stages, the qualitative traits associated with cotton fiber, and the set of genes involved in regulating these developmental stages and fiber traits. This review also highlights some prospects for how biotechnological approaches can improve cotton fiber quality.

GhSINA1, a SEVEN in ABSENTIA ubiquitin ligase, negatively regulates fiber development in Upland cotton

Protein ubiquitination is essential for plant growth and responses to the environment. The SEVEN IN ABSENTIA (SINA) ubiquitin ligases have been extensively studied in plants, but information on their roles in fiber development is limited. Here, we identified GhSINA1 in Upland cotton (Gossypium hirsutum), which has a conserved RING finger domain and SINA domain. Quantitative real-time PCR (qRT-PCR) analysis showed that GhSINA1 was preferentially expressed during fiber initiation and elongation, especially during initiation in the fuzzless-lintless cotton mutant. Subcellular localization experiments indicated that GhSINA1 localized to the nucleus. In vitro ubiquitination analysis revealed that GhSINA1 has E3 ubiquitin ligase activity. Ectopic overexpression of GhSINA1 in Arabidopsis thaliana reduced the number and length of root hairs and trichomes. Yeast two-hybrid (Y2H), firefly luciferase complementation imaging (LCI), and bimolecular fluorescence complementation (BiFC) assays demonstrated that the GhSINA1 proteins could interact with each other to form homodimers and heterodimers. Overall, these results suggest that GhSINA1 may act as a negative regulator in cotton fiber development through homodimerization and heterodimerization.

Deciphering the complex cotton genome for improving fiber traits and abiotic stress resilience in sustainable agriculture

BACKGROUND

Understanding the complex cotton genome is of paramount importance in devising a strategy for sustainable agriculture. Cotton is probably the most economically important cash crop known for its cellulose-rich fiber content. The cotton genome has become an ideal model for deciphering polyploidization due to its polyploidy, setting it apart from other major crops. However, the main challenge in understanding the functional and regulatory functions of many genes in cotton is still the complex cotton polyploidy genome, which is not limited to a single role. Cotton production is vulnerable to the sensitive effects of climate change, which can alter or aggravate soil, pests, and diseases. Thus, conventional plant breeding coupled with advanced technologies has led to substantial progress being made in cotton production.

GENOMICS APPROACHES IN COTTON

In the frontier areas of genomics research, cotton genomics has gained momentum accomplished by robust high-throughput sequencing platforms combined with novel computational tools to make the cotton genome more tractable. Advances in long-read sequencing have allowed for the generation of the complete set of cotton gene transcripts giving incisive scientific knowledge in cotton improvement. In contrast, the integration of the latest sequencing platforms has been used to generate multiple high-quality reference genomes in diploid and tetraploid cotton. While pan-genome and 3D genomic studies are still in the early stages in cotton, it is anticipated that rapid advances in sequencing, assembly algorithms, and analysis pipelines will have a greater impact on advanced cotton research.

CONCLUSIONS

This review article briefly compiles substantial contributions in different areas of the cotton genome, which include genome sequencing, genes, and their molecular regulatory networks in fiber development and stress tolerance mechanism. This will greatly help us in understanding the robust genomic organization which in turn will help unearth candidate genes for functionally important agronomic traits.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

Genome-wide analysis of the serine carboxypeptidase-like protein family reveals Ga09G1039 is involved in fiber elongation in cotton

The Gossypium is a model genus for understanding polyploidy and the evolutionary pattern of inheritance. This study aimed to investigate the characteristics of SCPLs in different cotton species and their role in fiber development. A total of 891 genes from one typical monocot and ten dicot species were naturally divided into three classes based on phylogenetic analysis. The SCPL gene family in cotton has undergone intense purifying selection with some functional variation. Segmental duplication and whole genome duplication were shown to be the two main reasons for the increase in the number of genes during cotton evolution. The identification of Gh_SCPL genes exhibiting differential expression in particular tissues or response to environmental stimuli provides a new measure for the in-depth characterization of selected genes of importance. Ga09G1039 was involved in the developmental process of fibers and ovules, and it is significantly different from proteins from other cotton species in terms of phylogenetic, gene structure, conserved protein motifs and tertiary structure. Overexpression of Ga09G1039 significantly increased the length of stem trichomes. Ga09G1039 may be a serine carboxypeptidase protein with hydrolase activity, according to functional region, prokaryotic expression, and western blotting analysis. The results provide a comprehensive overview of the genetic basis of SCPLs in Gossypium and further our knowledge in understanding the key aspects of SCPLs in cotton with their potential role in fiber development and stress resistance.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Characterization of chromatin accessibility and gene expression reveal the key genes involved in cotton fiber elongation

Cotton (Gossypium hirsutum L.) is an important economic crop, and cotton fiber is one of the longest plant cells, which provides an ideal model for the study of cell elongation and secondary cell wall synthesis. Cotton fiber length is regulated by a variety of transcription factors (TF) and their target genes; however, the mechanism of fiber elongation mediated by transcriptional regulatory networks is still unclear to a large extent. Here, we used a comparative assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) assay and RNA-seq analysis to identify fiber elongation transcription factors and genes using the short-fiber mutant ligon linless-2 (Li2 ) and wild type (WT). A total of 499 differential target genes were identified and GO analysis shows that differential genes are mainly involved in plant secondary wall synthesis and microtubule-binding processes. Analysis of the genomic regions preferentially accessible (Peak) has identified a number of overrepresented TF-binding motifs, highlighting sets of TFs that are important for cotton fiber development. Using ATAC-seq and RNA-seq data, we have constructed a functional regulatory network of each TF regulatory target gene and also the network pattern of TF regulating differential target genes. Further, to obtain the genes related to fiber length, the differential target genes were combined with FLGWAS data to identify the genes highly related to fiber length. Our work provides new insights into cotton fiber elongation.

The cotton protein GhIQD21 interacts with GhCaM7 and modulates organ morphogenesis in Arabidopsis by influencing microtubule stability

KEY MESSAGE

GhIQD21, a cotton IQ67-domain protein, interacts with GhCaM7 and alters organ shape in Arabidopsis by modulating microtubule stability. Calcium ion (Ca²⁺) and the calcium sensor calmodulin play crucial roles in the growth and development of plants. GhCaM7, a calmodulin in upland cotton (Gossypium hirsutum L.), is highly expressed in cotton fiber cells during the rapid elongation period and plays an important role in fiber cell development. In this study, we screened for GhCaM7-interacting proteins and identified GhIQD21, which contains a typical IQ67-domain. GhIQD21 was preferentially expressed at the fiber rapid elongation stage, and the protein localized to microtubules (MTs). Ectopic expression of GhIQD21 in Arabidopsis resulted in shorter leaves, petals, siliques, and plant height, thicker inflorescences, and more trichomes when compared with wild type (WT). Further investigation indicated that the morphogenesis of leaf epidermal cells and silique cells was altered. There was less consistency in the orientation of cortical microtubules in cotyledon and hypocotyl epidermal cells. Furthermore, compared with WT, transgenic seedling hypocotyls were more sensitive to oryzalin, a MT depolymerization drug. These results indicated that GhIQD21 is a GhCaM7-interacting protein located in MTs and that it plays a role in plant growth and potentially cotton fiber development. This study provides a foundation for further studies of the function and regulatory mechanism of GhIQD21 in fiber cell development.

Genetic regulation of nitrogen use efficiency in Gossypium spp

Cotton (Gossypium spp.) is the most important fibre crop, with desirable characteristics preferred for textile production. Cotton fibre output relies heavily on nitrate as the most important source of inorganic nitrogen (N). However, nitrogen dynamics in extreme environments limit plant growth and lead to yield loss and pollution. Therefore, nitrogen use efficiency (NUE), which involves the utilisation of the 'right rate', 'right source', 'right time', and 'right place' (4Rs), is key for efficient N management. Recent omics techniques have genetically improved NUE in crops. We herein highlight the mechanisms of N uptake and assimilation in the vegetative and reproductive branches of the cotton plant while considering the known and unknown regulatory factors. The phylogenetic relationships among N transporters in four Gossypium spp. have been reviewed. Further, the N regulatory genes that participate in xylem transport and phloem loading are also discussed. In addition, the functions of microRNAs and transcription factors in modulating the expression of target N regulatory genes are highlighted. Overall, this review provides a detailed perspective on the complex N regulatory mechanism in cotton, which would accelerate the research toward improving NUE in crops.

Brassinosteroids regulate cotton fiber elongation by modulating very-long-chain fatty acid biosynthesis

Brassinosteroid (BR), a growth-promoting phytohormone, regulates many plant growth processes including cell development. However, the mechanism by which BR regulates fiber growth is poorly understood. Cotton (Gossypium hirsutum) fibers are an ideal single-cell model in which to study cell elongation due to their length. Here we report that BR controls cotton fiber elongation by modulating very-long-chain fatty acid (VLCFA) biosynthesis. BR deficiency reduces the expression of 3-ketoacyl-CoA synthases (GhKCSs), the rate-limiting enzymes involved in VLCFA biosynthesis, leading to lower saturated VLCFA contents in pagoda1 (pag1) mutant fibers. In vitro ovule culture experiments show that BR acts upstream of VLCFAs. Silencing of BRI1-EMS-SUPPRESOR 1.4 (GhBES1.4), encoding a master transcription factor of the BR signaling pathway, significantly reduces fiber length, whereas GhBES1.4 overexpression produces longer fibers. GhBES1.4 regulates endogenous VLCFA contents and directly binds to BR RESPONSE ELEMENTS (BRREs) in the GhKCS10_At promoter region, which in turn regulates GhKCS10_At expression to increase endogenous VLCFA contents. GhKCS10_At overexpression promotes cotton fiber elongation, whereas GhKCS10_At silencing inhibits cotton fiber growth, supporting a positive regulatory role for GhKCS10_At in fiber elongation. Overall, these results uncover a mechanism of fiber elongation through crosstalk between BR and VLCFAs at the single-cell level.

Coordination of floral and fiber development in cotton (Gossypium) by hormone- and flavonoid-signalling associated regulatory miRNAs

Various plant development activities and stress responses are tightly regulated by various microRNAs (miRNA) and their target genes, or transcription factors in a spatiotemporal manner. Here, to exemplify how flowering-associated regulatory miRNAs synchronize their expression dynamics during floral and fiber development in cotton, constitutive expression diminution transgenic lines of auxin-signaling regulatory Gh-miR167 (35S-MIM167) were developed through target mimicry approach. 'Moderate' (58% to 80%)- and 'high' (> 80%)-Gh-miR167 diminution mimic lines showed dosage-dependent developmental deformities in anther development, pollen maturation, and fruit (= boll) formation. Cross pollination of 'moderate' 35S-MIM167 mimic lines with wild type (WT) plant partially restored boll formation and emergence of fiber initials on the ovule surface. Gh-miR167 diminution favored organ-specific transcription biases in miR159, miR166 as well as miR160, miR164, and miR172 along with their target genes during anther and petal development, respectively. Similarly, accumulative effect of percent Gh-miR167 diminution, cross regulation of its target ARF6/8 genes, and temporal mis-expression of hormone signaling- and flavonoid biosynthesis-associated regulatory miRNAs at early fiber initiation stage caused irregular fiber formation. Spatial and temporal transcription proportions of regulatory miRNAs were also found crucial for the execution of hormone- and flavonoid-dependent progression of floral and fiber development. These observations discover how assorted regulatory genetic circuits get organized in response to Gh-miR167 diminution and converge upon ensuing episodes of floral and fiber development in cotton.

Functional characterization of TBL genes revealed the role of GhTBL7 and GhTBL58 in cotton fiber elongation

TBL (Trichome Birefringence Like) gene family members are involved in trichome initiation and xylan acetylation in several plant species. In our research, we identified 102 TBLs from G. hirsutum. The phylogenetic tree classified TBL genes into five groups. Collinearity analysis of TBL genes indicated 136 paralogous gene pairs in G. hirsutum. Gene duplication indicated that WGD or segmental duplication contributed to the GhTBL gene family expansion. Promoter cis-elements of GhTBLs were related to growth and development, seed-specific regulation, light, and stress responses. GhTBL genes (GhTBL7, GhTBL15, GhTBL21, GhTBL25, GhTBL45, GhTBL54, GhTBL67, GhTBL72, and GhTBL77) exhibited upregulated response under exposure to cold, heat, NaCl, and PEG. GhTBL genes exhibited high expression during fiber development stages. Two GhTBL genes (GhTBL7 and GhTBL58) showed differential expression at 10 DPA fiber, as 10 DPA is a fast fiber elongation stage and fiber elongation is a very important stage of cotton fiber development. Subcellular localization of GhTBL7 and GhTBL58 revealed that these genes reside inside the cell membrane. Promoter GUS activity of GhTBL7 and GhTBL58 exhibited deep staining in roots. To further validate the role of these genes in cotton fiber elongation, we silenced these genes and observed a significant reduction in the fiber length at 10 DPA. In conclusion, the functional study of cell membrane-associated genes (GhTBL7 and GhTBL58) showed deep staining in root tissues and potential function during cotton fiber elongation at 10 DPA fiber.

Poplar coma morphogenesis and miRNA regulatory networks by combining ovary tissue sectioning and deep sequencing

Poplar coma, commonly referred to as "seed hairs", is a tuft of trichomes attached to the seed coat that helps seed dispersal. However, they can also trigger health impacts for humans, including sneezing, shortness of breath, and skin irritation. Despite efforts to study the regulatory mechanism of herbaceous trichome formation, poplar coma remains poorly understood. In this study, we showed that the epidermal cells of the funiculus and placenta are the origin of poplar coma based on observations of paraffin sections. Small RNA (sRNA) and degradome libraries were also constructed at three stages of poplar coma development, including initiation and elongation stages. Based on 7,904 miRNA-target pairs identified by small RNA and degradome sequencing, we constructed a miRNA-transcript factor and a stage-specific miRNA regulatory network. By combining paraffin section observation and deep sequencing, our research will provide greater insight into the molecular mechanisms of poplar coma development.

SMRT and Illumina sequencing provide insights into mechanisms of lignin and terpenoids biosynthesis in Pinus massoniana Lamb

Wood and oleoresin are important industrial raw materials with high economic value; however, their molecular formation and biosynthesis mechanisms in different tissues of Pinus massoniana remain unexplored. Therefore, we used single-molecule real-time sequencing technology (SMRT) and Illumina RNA sequencing to establish a transcriptome dataset and explore the expression pattern of genes related to secondary metabolites involved in wood formation and oleoresin biosynthesis in six different P. massoniana tissues. In total, 63.58 Gb of polymerase reads were obtained, including 41,407 isoforms with an average length of 1822 bp. We identified 3939 and 8785 isoforms and 161 and 481 transcription factors with tissue expression specificity and in the reproductive and vegetative organs, respectively. Eighty isoforms were annotated as cellulose synthases and 224 isoforms involved in lignin biosynthesis were enriched. Additionally, we identified 217 isoforms involved in the terpenoid biosynthesis pathway, with needles having the most tissue-specific genes for terpenoid biosynthesis. Some isoforms related to lignin biosynthesis were highly expressed in the xylem, according to the results of transcriptome sequencing and real-time quantitative reverse-transcription polymerase chain reaction. Our research confirmed the advantages of SMRT sequencing and provided valuable information for the transcriptional annotation of P. massoniana, which will be beneficial for producing better raw wood and oleoresin materials.

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

Cotton fiber provides the predominant plant textile in the world, and it is also a model for plant cell wall biosynthesis. The development of the single-celled cotton fiber takes place across several overlapping but discrete stages, including fiber initiation, elongation, the transition from elongation to secondary cell wall formation, cell wall thickening, and maturation and cell death. During each stage, the developing fiber undergoes a complex restructuring of genome-wide gene expression change and physiological/biosynthetic processes, which ultimately generate a strikingly elongated and nearly pure cellulose product that forms the basis of the global cotton industry. Here, we provide an overview of this developmental process focusing both on its temporal as well as evolutionary dimensions. We suggest potential avenues for further improvement of cotton as a crop plant.

Tissue specific expression of bacterial cellulose synthase (Bcs) genes improves cotton fiber length and strength

Fiber growing inside the cotton bolls is a highly demandable product and its quality is key to the success of the textile industry. Despite the various efforts to improve cotton fiber staple length Pakistan has to import millions of bales to sustain its industrial needs. To improve cotton fiber quality Bacterial cellulose synthase (Bcs) genes (acsA, acsB) were expressed in a local cotton variety CEMB-00. In silico studies revealed a number of conserved domains both in the cotton-derived and bacterial cellulose synthases which are essential for the cellulose synthesis. Transformation efficiency of 1.27% was achieved by using Agrobacterium shoot apex cut method of transformation. The quantitative mRNA expression analysis of the Bcs genes in transgenic cotton fiber was found to be many folds higher during secondary cell wall synthesis stage (35 DPA) than the expression during elongation phase (10 DPA). Average fiber length of the transgenic cotton plant lines S-00-07, S-00-11, S-00-16 and S-00-23 was calculated to be 13.02% higher than that of the non-transgenic control plants. Likewise, the average fiber strength was found to be 20.92% higher with an enhanced cellulose content of 22.45%. The mutated indigenous cellulose synthase genes of cotton generated through application of CRISPR/Cas9 resulted in 6.03% and 12.10% decrease in fiber length and strength respectively. Furthermore, mature cotton fibers of transgenic cotton plants were found to have increased number of twists with smooth surface as compared to non-transgenic control when analyzed under scanning electron microscope. XRD analysis of cotton fibers revealed less cellulose crystallinity index in transgenic cotton fibers as compared to control fibers due to deposition of more amorphous cellulose in transgenic fibers as a result of Bcs gene expression. This study paved the way towards unraveling the fact that Bcs genes influence cellulose synthase activity and this enzyme helps in determining the fate of cotton fiber length and strength.

Visualization of Cytoskeleton Organization and Dynamics in Elongating Cotton Fibers by Live-Cell Imaging

Cotton fibers are extremely elongated single cells and have long been regarded as an ideal model to investigate polarized plant cell elongation. Actin filaments (F-actin), as well as the cortical microtubules (CMTs), play vital roles in polarized cell growth and morphogenesis. We have generated stable transgenic cotton plants expressing fluorescent markers for the actin and microtubule cytoskeletons. Further live-cell imaging identified dynamic features of the F-actin and cortical microtubule (CMT) architectures and discovered that cotton fibers elongate in a unique tip-biased diffuse growth mode. Here, we describe methods for preparing growing cotton fiber samples, as well as the visualization of cytoskeletal organization and dynamics by live-cell imaging. Combined with comprehensive image analyses, these methods can be used to identify how cytoskeleton organization and dynamics determine cell morphogenesis in highly polarized cotton fibers.

Cotton microtubule-associated protein GhMAP20L5 mediates fiber elongation through the interaction with the tubulin GhTUB13

Targeting proteins for Xklp2 (TPX2s) comprise a class of MAPs that are essential for plant growth and development by regulating the dynamic changes of microtubules (MTs) and proper formation of cytoskeleton. However, the function of TPX2 proteins in cotton fiber development remains poorly understood. Here, we identified the function of a fiber elongation-specific TPX2 protein, GhMAP20L5, in cotton. Suppressed GhMAP20L5 gene expression in cotton (GhMAP20L5i) significantly reduced fiber elongation rate, fiber length and lint percentage. GhMAP20L5i fibers had thinner and looser secondary cell walls (SCW), and incompact helix twists. GhMAP20L5 specifically interacted with the tubulin GhTUB13 on the cytoskeleton. Gene coexpression analysis showed that GhMAP20L5 involved in multiple pathways related to cytoskeleton establishment and fiber cell wall formation and affected cellulase genes expressions. In summary, our results revealed that GhMAP20L5 is important for fiber development by regulating cytoskeleton establishment and the cellulose deposition in cotton.

Comparative proteomic analysis identified proteins and the phenylpropanoid biosynthesis pathway involved in the response to ABA treatment in cotton fiber development

Abscisic acid (ABA) is a plant hormone that plays an important role in cotton fiber development. In this study, the physiological changes and proteomic profiles of cotton (Gossypium hirsutum) ovules were analyzed after 20 days of ABA or ABA inhibitor (ABAI) treatment. The results showed that compared to the control (CK), the fiber length was significantly decreased under ABA treatment and increased under ABAI treatment. Using a tandem mass tags-based quantitative technique, the proteomes of cotton ovules were comprehensively analyzed. A total of 7321 proteins were identified, of which 365 and 69 differentially accumulated proteins (DAPs) were identified in ABA versus CK and ABAI versus CK, respectively. Specifically, 345 and 20 DAPs were up- and down-regulated in the ABA group, and 65 and 4 DAPs were up- and down-regulated in the ABAI group, respectively. The DAPs in the ABA group were mainly enriched in the biosynthesis of secondary metabolites, phenylpropanoid biosynthesis and flavonoid secondary metabolism, whereas the DAPs in the ABAI group were mainly enriched in the indole alkaloid biosynthesis and phenylpropanoid biosynthesis pathways. Moreover, 9 proteins involved in phenylpropanoid biosynthesis were upregulated after ABA treatment, suggesting that this pathway might play important roles in the response to ABA, and 3 auxin-related proteins were upregulated, indicating that auxin might participate in the regulation of fiber development under ABAI treatment.

GhROP6 GTPase modulates auxin accumulation in cotton fibers by regulating cell-specific GhPIN3a localization

PIN-FORMED- (PIN) mediated polar auxin transport plays a predominant role in most auxin-triggered organogenesis in plants. Global control of PIN polarity at the plasma membrane contributes to the essential establishment of auxin maxima in most multicellular tissues. However, establishment of auxin maxima in single cells is poorly understood. Cotton fibers, derived from ovule epidermal cells by auxin-triggered cell protrusion, provide an ideal model to explore the underlying mechanism. Here, we report that cell-specific degradation of GhPIN3a, which guides the establishment of the auxin gradient in cotton ovule epidermal cells, is associated with the preferential expression of GhROP6 GTPase in fiber cells. In turn, GhROP6 reduces GhPIN3a abundance at the plasma membrane and facilitates intracellular proteolysis of GhPIN3a. Overexpression and activation of GhROP6 promote cell elongation, resulting in a substantial improvement in cotton fiber length.

Genome-wide analysis elucidates the roles of GhHMA genes in different abiotic stresses and fiber development in upland cotton

The heavy metal-binding domain is involved in heavy metal transporting and plays a significant role in plant detoxification. However, the functions of HMAs are less well known in cotton. In this study, a total of 143 GhHMAs (heavy metal-binding domain) were detected by genome-wide identification in G. hirsutum L. All the GhHMAs were classified into four groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GhHMA genes was highly conserved. 212 paralogous GhHMA gene pairs were identified, and the segmental duplications were the main role to the expansion of GhHMAs. The Ka/Ks values suggested that the GhHMA gene family has undergone purifying selection during the long-term evolutionary process. GhHMA3 and GhHMA75 were located in the plasma membrane, while GhHMA26, GhHMA117 and GhHMA121 were located in the nucleus, respectively. Transcriptomic data and qRT-PCR showed that GhHMA26 exhibited different expression patterns in each tissue and during fiber development or under different abiotic stresses. Overexpressing GhHMA26 significantly promoted the elongation of leaf trichomes and also improved the tolerance to salt stress. Therefore, GhHMA26 may positively regulate fiber elongation and abiotic stress. Yeast two-hybrid assays indicated that GhHMA26 and GhHMA75 participated in multiple biological functions. Our results suggest some genes in the GhHMAs might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GhHMA genes in cotton.

Phenomics and transcriptomics analyses reveal deposition of suberin and lignin in the short fiber cell walls produced from a wild cotton species and two mutants

Fiber length is one of the major properties determining the quality and commercial value of cotton. To understand the mechanisms regulating fiber length, genetic variations of cotton species and mutants producing short fibers have been compared with cultivated cottons generating long and normal fibers. However, their phenomic variation other than fiber length has not been well characterized. Therefore, we compared physical and chemical properties of the short fibers with the long fibers. Fiber characteristics were compared in two sets: 1) wild diploid Gossypium raimondii Ulbrich (short fibers) with cultivated diploid G. arboreum L and tetraploid G. hirsutum L. (long fibers); 2) G. hirsutum short fiber mutants, Ligon-lintless 1 (Li1) and 2 (Li2) with their near isogenic line (NIL), DP-5690 (long fibers). Chemical analyses showed that the short fibers commonly consisted of greater non-cellulosic components, including lignin and suberin, than the long fibers. Transcriptomic analyses also identified up-regulation of the genes related to suberin and lignin biosynthesis in the short fibers. Our results may provide insight on how high levels of suberin and lignin in cell walls can affect cotton fiber length. The approaches combining phenomic and transcriptomic analyses of multiple sets of cotton fibers sharing a common phenotype would facilitate identifying genes and common pathways that significantly influence cotton fiber properties.

Transcriptomic analysis reveals the key role of histone deacetylation via mediating different phytohormone signalings in fiber initiation of cotton

Background

Histone deacetylation is one of the most important epigenetic modifications and plays diverse roles in plant development. However, the detailed functions and mechanisms of histone deacetylation in fiber development of cotton are still unclear. HDAC inhibitors (HDACi) have been commonly used to study the molecular mechanism underlying histone deacetylation or to facilitate disease therapy in humans through hindering the histone deacetylase catalytic activity. Trichostatin A (TSA)—the most widely used HDACi has been extensively employed to determine the role of histone deacetylation on different developmental stages of plants.

Results

Through in vitro culture of ovules, we observed that exogenous application of TSA was able to inhibit the fiber initiation development. Subsequently, we performed a transcriptomic analysis to reveal the underlying mechanisms. The data showed that TSA treatment resulted in 4209 differentially expressed genes, which were mostly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, and carbon metabolism pathways. The phytohormone signal transduction pathways harbor the most differentially expressed genes. Deeper studies showed that some genes promoting auxin, Gibberellic Acid (GA) signaling were down-regulated, while some genes facilitating Abscisic Acid (ABA) and inhibiting Jasmonic Acid (JA) signaling were up-regulated after the TSA treatments. Further analysis of plant hormone contents proved that TSA significantly promoted the accumulation of ABA, JA and GA₃.

Conclusions

Collectively, histone deacetylation can regulate some key genes involved in different phytohormone pathways, and consequently promoting the auxin, GA, and JA signaling, whereas repressing the ABA synthesis and signaling to improve the fiber cell initiation. Moreover, the genes associated with energy metabolism, phenylpropanoid, and glutathione metabolism were also regulated by histone deacetylation. The above results provided novel clues to illuminate the underlying mechanisms of epigenetic modifications as well as related different phytohormones in fiber cell differentiation, which is also very valuable for the molecular breeding of higher quality cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00840-4.

Morphogenesis and cell wall composition of trichomes and their function in response to salt in halophyte Salsola ferganica

BACKGROUND

To survive harsh environmental conditions, desert plants show various adaptions, such as the evolution of trichomes, which are protective epidermal protrusions. Currently, the morphogenesis and function of trichomes in desert plants are not well understood. Salsola ferganica is an annual halophyte distributed in cold deserts; at the seedling stage, its rod-shaped true leaves are covered with long and thick trichomes and are affected by habitat conditions. Therefore, we evaluated the trichomes on morphogenesis and cell wall composition of S. ferganica compared to Arabidopsis thaliana and cotton, related gene expression, and preliminary function in salt accumulation of the leaves.

RESULTS

The trichomes of S. ferganica were initiated from the epidermal primordium, followed by two to three rounds of cell division to form a multicellular trichome, while some genes associated with them were positively involved. Cell wall composition analysis showed that different polysaccharides including heavily methyl-esterified and fully de-esterified pectins (before maturation, probably in the primary wall), xyloglucans (in the mid-early and middle stages, probably in the secondary wall), and extensin (during the whole developmental period) were detected, which were different from those found in trichomes of Arabidopsis and cotton. Moreover, trichome development was affected by abiotic stress, and might accumulate salt from the mesophyll cells and secrete outside.

CONCLUSIONS

S. ferganica has multicellular, non-branched trichomes that undergo two to three rounds of cell division and are affected by abiotic stress. They have a unique cell wall composition which is different from that of Arabidopsis and cotton. Furthermore, several genes positively or negatively regulate trichome development. Our findings should contribute to our further understanding of the biogenesis and adaptation of plant accessory structures in desert plant species.

Strigolactones act downstream of gibberellins to regulate fiber cell elongation and cell wall thickness in cotton (Gossypium hirsutum)

Strigolactones (SLs) are a class of phytohormones that regulate plant shoot branching and adventitious root development. However, little is known regarding the role of SLs in controlling the behavior of the smallest unit of the organism, the single cell. Here, taking advantage of a classic single-cell model offered by the cotton (Gossypium hirsutum) fiber cell, we show that SLs, whose biosynthesis is fine-tuned by gibberellins (GAs), positively regulate cell elongation and cell wall thickness by promoting the biosynthesis of very long-chain fatty acids (VLCFAs) and cellulose, respectively. Furthermore, we identified two layers of transcription factors (TFs) involved in the hierarchical regulation of this GA-SL crosstalk. The top-layer TF GROWTH-REGULATING FACTOR 4 (GhGRF4) directly activates expression of the SL biosynthetic gene DWARF27 (D27) to increase SL accumulation in fiber cells and GAs induce GhGRF4 expression. SLs induce the expression of four second-layer TF genes (GhNAC100-2, GhBLH51, GhGT2, and GhB9SHZ1), which transmit SL signals downstream to two ketoacyl-CoA synthase genes (KCS) and three cellulose synthase (CesA) genes by directly activating their transcription. Finally, the KCS and CesA enzymes catalyze the biosynthesis of VLCFAs and cellulose, respectively, to regulate development of high-grade cotton fibers. In addition to providing a theoretical basis for cotton fiber improvement, our results shed light on SL signaling in plant development at the single-cell level.

Carpel-specific down-regulation of GhCKXs in cotton significantly enhances seed and fiber yield

Cytokinin is considered to be an important driver of seed yield. To increase the yield of cotton while avoiding the negative consequences caused by constitutive overproduction of cytokinin, we down-regulated specifically the carpel genes for cytokinin oxidase/dehydrogenase (CKX), a key negative regulator of cytokinin levels, in transgenic cotton. The carpel-specific down-regulation of CKXs significantly enhanced cytokinin levels in the carpels. The elevated cytokinin promoted the expression of carpel- and ovule-development-associated genes, GhSTK2, GhAG1, and GhSHP, boosting ovule formation and thus producing more seeds in the ovary. Field experiments showed that the carpel-specific increase of cytokinin significantly increased both seed yield and fiber yield of cotton, without resulting in detrimental phenotypes. Our study details the regulatory mechanism of cytokinin signaling for seed development, and provides an effective and feasible strategy for yield improvement of seed crops.

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

A deletion/duplication in the Ligon lintless-2 locus induces siRNAs that inhibit cotton fiber cell elongation

Most cultivated cotton (Gossypium hirsutum L.) varieties have two types of seed fibers: short fuzz fiber strongly adhered to the seed coat, and long lint fiber used in the textile industry. The Ligon lintless-2 (Li2) cotton mutant has a normal vegetative phenotype but produces very short lint fiber on the seeds. The Li2 mutation is controlled by a single dominant gene. We discovered a large structural rearrangement at the end of chromosome D13 in the Li2 mutant based on whole-genome sequencing and genetic mapping of segregating populations. The rearrangement contains a 177-kb deletion and a 221-kb duplication positioned as a tandem inverted repeat. The gene Gh_D13G2437 is located at the junction of the inverted repeat in the duplicated region. During transcription such structure spontaneously forms self-complementary hairpin RNA of Gh_D13G2437 followed by production of small interfering RNA (siRNA). Gh_D13G2437 encodes a Ran-Binding Protein 1 (RanBP1) that preferentially expresses during cotton fiber elongation. The abundance of siRNA produced from Gh_D13G2437 reciprocally corresponds with the abundance of highly homologous (68%-98% amino acid sequence identity) RanBP1 family transcripts during fiber elongation, resulting in a shorter fiber phenotype in the Li2. Overexpression of Gh_D13G2437 in the Li2 mutant recovered the long lint fiber phenotype. Taken together, our findings revealed that siRNA-induced silencing of a family of RanBP1s inhibit elongation of cotton fiber cells in the Li2 mutant.