Recent Advances and Future Perspectives in Cotton Research

GhHDZ76, a cotton HD-Zip transcription factor, involved in regulating the initiation and early elongation of cotton fiber development in G. hirsutum

In this study, the whole HD-Zip family members of G. hirsutum were identified, and GhHDZ76 was classified into the HD-Zip IV subgroup. GhHDZ76 was predominantly expressed in the 0-5 DPA of fiber development stage and localized in the nucleus. Overexpression of GhHDZ76 significantly increased the length and density of trichomes in Arabidopsis thaliana. The fiber length of GhHDZ76 knockout lines by CRISPR/Cas9 was significantly shorter than WT at the early elongation and mature stage, indicating that GhHDZ76 positively regulate the fiber elongation. Scanning electron microscopy showed that the number of ovule surface protrusion of 0 DPA of GhHDZ76 knockout lines was significantly lower than WT, suggesting that GhHDZ76 can also promote the initiation of fiber development. The transcript level of GhWRKY16, GhRDL1, GhEXPA1 and GhMYB25 genes related to fiber initiation and elongation in GhHDZ76 knockout lines were significantly decreased. Yeast two-hybrid and Luciferase complementation imaging (LCI) assays showed that GhHDZ76 can interact with GhWRKY16 directly. As a transcription factor, GhHDZ76 has transcriptional activation activity, which could bind to L1-box elements of the promoters of GhRDL1 and GhEXPA1. Double luciferase reporter assay showed that the GhWRKY16 could enhance the transcriptional activity of GhHDZ76 to pGhRDL1, but it did not promote the transcriptional activity of GhHDZ76 to pGhEXPA1. GhHDZ76 protein may also promote the transcriptional activity of GhWRKY16 to the downstream target gene GhMYB25. Our results provided a new gene resource for fiber development and a theoretical basis for the genetic improvement of cotton fiber quality.

Sphingolipid inhibitor response gene GhMYB86 controls fiber elongation by regulating microtubule arrangement

Although the cell membrane and cytoskeleton play essential roles in cellular morphogenesis, the interaction between the membrane and cytoskeleton is poorly understood. Cotton fibers are extremely elongated single cells, which makes them an ideal model for studying cell development. Here, we used the sphingolipid biosynthesis inhibitor, fumonisin B1 (FB1), and found that it effectively suppressed the myeloblastosis (MYB) transcription factor GhMYB86, thereby negatively affecting fiber elongation. A direct target of GhMYB86 is GhTUB7, which encodes the tubulin protein, the major component of the microtubule cytoskeleton. Interestingly, both the overexpression of GhMYB86 and GhTUB7 caused an ectopic microtubule arrangement at the fiber tips, and then leading to shortened fibers. Moreover, we found that GhMBE2 interacted with GhMYB86 and that FB1 and reactive oxygen species induced its transport into the nucleus, thereby enhancing the promotion of GhTUB7 by GhMYB86. Overall, we established a GhMBE2-GhMYB86-GhTUB7 regulation module for fiber elongation and revealed that membrane sphingolipids affect fiber elongation by altering microtubule arrangement.

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.

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.

Screening cotton genotypes for their drought tolerance ability based on the expression level of dehydration-responsive element-binding protein and proline biosynthesis-related genes and morpho-physio-biochemical responses

Drought stress adversely affects growth, development, productivity, and fiber quality of cotton (Gossypium hirsutum L). Breeding strategies to enhance drought tolerance require an improved knowledge of plant drought responses necessitating proper identification of drought-tolerant genotypes of crops, including cotton. The objective of this study was to classify the selected cotton genotypes for their drought tolerance ability based on morpho-physio-biochemical traits using Hierarchical Ward's cluster analysis. Five genotypes of cotton (Takfa 3, Takfa 6, Takfa 7, Takfa 84-4, and Takfa 86-5) were selected as plant materials, and were grown under well-watered (WW; 98 ± 2% field capacity) and water-deficit (WD; 50 ± 2% field capacity) conditions for 16 days during the flower initiation stage. Data on morpho-physio-biochemical parameters and gene expression levels for these parameters were collected, and subsequently genotypes were classified either as a drought tolerant or drought susceptible one. Upregulation of GhPRP (proline-rich protein), GhP5CS (Δ1-pyrroline-5-carboxylate synthetase), and GhP5CR (Δ1-pyrroline-5-carboxylate reductase) in relation to free proline enrichment was observed in Takfa 3 genotype under WD condition. An accumulation of free proline, total soluble sugar, and potassium in plants under WD conditions was detected, which played a key role as major osmolytes controlling cellular osmotic potential. Magnesium and calcium concentrations were also enriched in leaves under WD conditions, functioning as essential elements and regulating photosynthetic abilities. Leaf greenness, net photosynthetic rate, stomatal conductance, and transpiration rate were also declined under WD conditions, leading to growth retardation, especially aboveground traits of Takfa 6, Takfa 7, Takfa 84-4, and Takfa 86-5 genotypes. An increase in leaf temperature (1.1 - 4.0 °C) and crop water stress index (CWSI > 0.75) in relation to stomatal closure and reduced transpiration rate was recorded in cotton genotypes under WD conditions compared with WW conditions. Based on the increase of free proline, soluble sugar, leaf temperature, and CWSI, as well as the decrease of aboveground growth traits and physiological attributes, five genotypes were categorized into two cluster groups: drought tolerant (Takfa 3) and drought susceptible (Takfa 6, Takfa 7, Takfa 84-4, and Takfa 86-5). The identified drought-tolerant cotton genotype, namely, Takfa 3, may be grown in areas experiencing drought conditions. It is recommended to further validate the yield traits of Takfa 3 under rainfed field conditions in drought-prone environments.

Auxin resistant 2 and short hypocotyl 2 regulate cotton fiber initiation and elongation

Auxin, a pivotal regulator of diverse plant growth processes, remains central to development. The auxin-responsive genes auxin/indole-3-acetic acids (AUX/IAAs) are indispensable for auxin signal transduction, which is achieved through intricate interactions with auxin response factors (ARFs). Despite this, the potential of AUX/IAAs to govern the development of the most fundamental biological unit, the single cell, remains unclear. In this study, we harnessed cotton (Gossypium hirsutum) fiber, a classic model for plant single-cell investigation, to determine the complexities of AUX/IAAs. Our research identified 2 pivotal AUX/IAAs, auxin resistant 2 (GhAXR2) and short hypocotyl 2 (GhSHY2), which exhibit opposite control over fiber development. Notably, suppressing GhAXR2 reduced fiber elongation, while silencing GhSHY2 fostered enhanced fiber elongation. Investigating the mechanistic intricacies, we identified specific interactions between GhAXR2 and GhSHY2 with distinct ARFs. GhAXR2's interaction with GhARF6-1 and GhARF23-2 promoted fiber cell development through direct binding to the AuxRE cis-element in the constitutive triple response 1 promoter, resulting in transcriptional inhibition. In contrast, the interaction of GhSHY2 with GhARF7-1 and GhARF19-1 exerted a negative regulatory effect, inhibiting fiber cell growth by activating the transcription of xyloglucan endotransglucosylase/hydrolase 9 and cinnamate-4-hydroxylase. Thus, our study reveals the intricate regulatory networks surrounding GhAXR2 and GhSHY2, elucidating the complex interplay of multiple ARFs in AUX/IAA-mediated fiber cell growth. This work enhances our understanding of single-cell development and has potential implications for advancing plant growth strategies and agricultural enhancements.

GhGASA14 regulates the flowering time of upland cotton in response to GA3

KEY MESSAGE

The silencing of GhGASA14 and the identification of superior allelic variation in its coding region indicate that GhGASA14 may positively regulate flowering and the response to GA3. Gibberellic acid-stimulated Arabidopsis (GASA), a member of the gibberellin-regulated short amino acid family, has been extensively investigated in several plant species and found to be critical for plant growth and development. However, research on this topic in cotton has been limited. In this study, we identified 38 GhGASAs that were dispersed across 18 chromosomes in upland cotton, and all of these genes had a GASA core domain. Transcriptome expression patterns and qRT-PCR results revealed that GhGASA9 and GhGASA14 exhibited upregulated expression not only in the floral organs but also in the leaves of early-maturing cultivars. The two genes were functionally characterized by virus-induced gene silencing (VIGS), and the budding and flowering times after silencing the target genes were later than those of the control (TRV:00). Compared with that in the water-treated group (MOCK), the flowering period of the different fruiting branches in the GA3-treated group was more concentrated. Interestingly, allelic variation was detected in the coding sequence of GhGASA14 between early-maturing and late-maturing accessions, and the frequency of this favorable allele was greater in high-latitude cotton cultivars than in low-latitude ones. Additionally, a significant linear relationship was observed between the expression level of GhGASA14 and flowering time among the 12 upland cotton accessions. Taken together, these results indicated that GhGASA14 may positively regulate flowering time and respond to GA3. These findings could lead to the use of valuable genetic resources for breeding early-maturing cotton cultivars in the future.

Chromosome-level genome assembly and annotation of a potential model organism Gossypium arboreum ZB-1

Recent advancements in plant regeneration and synthetic polyploid creation have been documented in Gossypium arboreum ZB-1. These developments make ZB-1 a potential model within the Gossypium genus for investigating gene function and polyploidy. This work generated the sequence and annotation of the ZB-1 genome. The contig-level genome was constructed using the PacBio high-fidelity reads, encompassing 81 contigs with an N50 length of 112.12 Mb. The Hi-C data assisted the construction of the chromosome-level genome, which consists of 13 pseudo-chromosomes and 39 un-anchored contigs, with a total length of about 1.67 Gb. Repetitive sequences accounted for about 69.7% of the genome in length. Based on ab initio and evidence-based prediction, we have identified 48,021 protein-coding genes in the ZB-1 genome. Comparative genomics analysis revealed conserved gene content and arrangement between ZB-1 and G. arboreum SXY1. The single nucleotide polymorphism occurrence rate between ZB-1 and SXY1 was about 0.54 per 1,000 nucleotides. This study enriched the genomic resources for further exploration into cotton regeneration and polyploidy mechanisms.

Lightweight cotton diseases real-time detection model for resource-constrained devices in natural environments

Cotton, a vital textile raw material, is intricately linked to people's livelihoods. Throughout the cotton cultivation process, various diseases threaten cotton crops, significantly impacting both cotton quality and yield. Deep learning has emerged as a crucial tool for detecting these diseases. However, deep learning models with high accuracy often come with redundant parameters, making them challenging to deploy on resource-constrained devices. Existing detection models struggle to strike the right balance between accuracy and speed, limiting their utility in this context. This study introduces the CDDLite-YOLO model, an innovation based on the YOLOv8 model, designed for detecting cotton diseases in natural field conditions. The C2f-Faster module replaces the Bottleneck structure in the C2f module within the backbone network, using partial convolution. The neck network adopts Slim-neck structure by replacing the C2f module with the GSConv and VoVGSCSP modules, based on GSConv. In the head, we introduce the MPDIoU loss function, addressing limitations in existing loss functions. Additionally, we designed the PCDetect detection head, integrating the PCD module and replacing some CBS modules with PCDetect. Our experimental results demonstrate the effectiveness of the CDDLite-YOLO model, achieving a remarkable mean average precision (mAP) of 90.6%. With a mere 1.8M parameters, 3.6G FLOPS, and a rapid detection speed of 222.22 FPS, it outperforms other models, showcasing its superiority. It successfully strikes a harmonious balance between detection speed, accuracy, and model size, positioning it as a promising candidate for deployment on an embedded GPU chip without sacrificing performance. Our model serves as a pivotal technical advancement, facilitating timely cotton disease detection and providing valuable insights for the design of detection models for agricultural inspection robots and other resource-constrained agricultural devices.

Wall-associated kinase GhWAK13 mediates arbuscular mycorrhizal symbiosis and Verticillium wilt resistance in cotton

The cell wall is the major interface for arbuscular mycorrhizal (AM) symbiosis. However, the roles of cell wall proteins and cell wall synthesis in AM symbiosis remain unclear. We reported that a novel wall-associated kinase 13 (GhWAK13) positively regulates AM symbiosis and negatively regulates Verticillium wilt resistance in cotton. GhWAK13 transcription was induced by AM symbiosis and Verticillium dahliae (VD) infection. GhWAK13 is located in the plasma membrane and expressed in the arbuscule-containing cortical cells of mycorrhizal cotton roots. GhWAK13 silencing inhibited AM colonization and repressed gene expression of the mycorrhizal pathway. Moreover, GhWAK13 silencing improved Verticillium wilt resistance and triggered the expression of immunity genes. Therefore, GhWAK13 is considered an immune suppressor required for AM symbiosis and disease resistance. GhWAK7A, a positive regulator of Verticillium wilt resistance, was upregulated in GhWAK13-silenced cotton plants. Silencing GhWAK7A improved AM symbiosis. Oligogalacturonides application also suppressed AM symbiosis. Finally, GhWAK13 negatively affected the cellulose content by regulating the transcription of cellulose synthase genes. The results of this study suggest that immunity suppresses AM symbiosis in cotton. GhWAK13 affects AM symbiosis by suppressing immune responses.

Upgraded durian genome reveals the role of chromosome reshuffling during ancestral karyotype evolution, lignin biosynthesis regulation, and stress tolerance

Durian (Durio zibethinus) is a tropical fruit that has a unique flavor and aroma. It occupies a significant phylogenetic position within the Malvaceae family. Extant core-eudicot plants are reported to share seven ancestral karyotypes that have undergone reshuffling, resulting in an abundant genomic diversity. However, the ancestral karyotypes of the Malvaceae family, as well as the evolution trajectory leading to the 28 chromosomes in durian, remain poorly understood. Here, we report the high-quality assembly of the durian genome with comprehensive comparative genomic analyses. By analyzing the collinear blocks between cacao and durian, we inferred 11 Malvaceae ancestral karyotypes. These blocks were present in a single-copy form in cacao and mainly in triplicates in durian, possibly resulting from a recent whole genome triplication (WGT) event that led to hexaploidization of the durian genome around 20 (17-24) million years ago. A large proportion of the duplicated genes in durian, such as those involved in the lignin biosynthesis module for phenylpropane biosynthesis, are derived directly from whole genome duplication, which makes it an important force in reshaping its genomic architecture. Transcriptome studies have revealed that genes involved in feruloyl-CoA formations were highly preferentially expressed in fruit peels, indicating that the thorns produced on durian fruit may comprise guaiacyl and syringyl lignins. Among all the analyzed transcription factors (TFs), members of the heat shock factor family (HSF) were the most significantly upregulated under heat stress. All subfamilies of genes encoding heat shock proteins (HSPs) in the durian genome appear to have undergone expansion. The potential interactions between HSF Dzi05.397 and HSPs were examined and experimentally verified. Our study provides a high-quality durian genome and reveals the reshuffling mechanism of ancestral Malvaceae chromosomes to produce the durian genome. We also provide insights into the mechanism underlying lignin biosynthesis and heat stress tolerance.

GhFAD3-4 Promotes Fiber Cell Elongation and Cell Wall Thickness by Increasing PI and IP3 Accumulation in Cotton

The omega-3 fatty acid desaturase enzyme gene FAD3 is responsible for converting linoleic acid to linolenic acid in plant fatty acid synthesis. Despite limited knowledge of its role in cotton growth, our study focused on GhFAD3-4, a gene within the FAD3 family, which was found to promote fiber elongation and cell wall thickness in cotton. GhFAD3-4 was predominantly expressed in elongating fibers, and its suppression led to shorter fibers with reduced cell wall thickness and phosphoinositide (PI) and inositol triphosphate (IP3) levels. Transcriptome analysis of GhFAD3-4 knock-out mutants revealed significant impacts on genes involved in the phosphoinositol signaling pathway. Experimental evidence demonstrated that GhFAD3-4 positively regulated the expression of the GhBoGH3B and GhPIS genes, influencing cotton fiber development through the inositol signaling pathway. The application of PI and IP6 externally increased fiber length in GhFAD3-4 knock-out plants, while inhibiting PI led to a reduced fiber length in GhFAD3-4 overexpressing plants. These findings suggest that GhFAD3-4 plays a crucial role in enhancing fiber development by promoting PI and IP3 biosynthesis, offering the potential for breeding cotton varieties with superior fiber quality.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton

BACKGROUND

Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors.

RESULTS

In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms.

CONCLUSIONS

These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

Insights into the biocontrol and plant growth promotion functions of Bacillus altitudinis strain KRS010 against Verticillium dahliae

BACKGROUND

Verticillium wilt, caused by the fungus Verticillium dahliae, is a soil-borne vascular fungal disease, which has caused great losses to cotton yield and quality worldwide. The strain KRS010 was isolated from the seed of Verticillium wilt-resistant Gossypium hirsutum cultivar "Zhongzhimian No. 2."

RESULTS

The strain KRS010 has a broad-spectrum antifungal activity to various pathogenic fungi as Verticillium dahliae, Botrytis cinerea, Fusarium spp., Colletotrichum spp., and Magnaporthe oryzae, of which the inhibition rate of V. dahliae mycelial growth was 73.97% and 84.39% respectively through confrontation test and volatile organic compounds (VOCs) treatments. The strain was identified as Bacillus altitudinis by phylogenetic analysis based on complete genome sequences, and the strain physio-biochemical characteristics were detected, including growth-promoting ability and active enzymes. Moreover, the control efficiency of KRS010 against Verticillium wilt of cotton was 93.59%. After treatment with KRS010 culture, the biomass of V. dahliae was reduced. The biomass of V. dahliae in the control group (Vd991 alone) was 30.76-folds higher than that in the treatment group (KRS010+Vd991). From a molecular biological aspect, KRS010 could trigger plant immunity by inducing systemic resistance (ISR) activated by salicylic acid (SA) and jasmonic acid (JA) signaling pathways. Its extracellular metabolites and VOCs inhibited the melanin biosynthesis of V. dahliae. In addition, KRS010 had been characterized as the ability to promote plant growth.

CONCLUSIONS

This study indicated that B. altitudinis KRS010 is a beneficial microbe with a potential for controlling Verticillium wilt of cotton, as well as promoting plant growth.

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

GhUBC10-2 mediates GhGSTU17 degradation to regulate salt tolerance in cotton (Gossypium hirsutum)

Ubiquitin-conjugating enzyme (UBC) is a crucial component of the ubiquitin-proteasome system, which contributes to plant growth and development. While some UBCs have been identified as potential regulators of abiotic stress responses, the underlying mechanisms of this regulation remain poorly understood. Here, we report a cotton (Gossypium hirsutum) UBC gene, GhUBC10-2, which negatively regulates the salt stress response. We found that the gain of function of GhUBC10-2 in both Arabidopsis (Arabidopsis thaliana) and cotton leads to reduced salinity tolerance. Additionally, GhUBC10-2 interacts with glutathione S-transferase (GST) U17 (GhGSTU17), forming a heterodimeric complex that promotes GhGSTU17 degradation. Intriguingly, GhUBC10-2 can be self-polyubiquitinated, suggesting that it possesses E3-independent activity. Our findings provide new insights into the PTM of plant GST-mediated salt response pathways. Furthermore, we found that the WRKY transcription factor GhWRKY13 binds to the GhUBC10-2 promoter and suppresses its expression under salt conditions. Collectively, our study unveils a regulatory module encompassing GhWRKY13-GhUBC10-2-GhGSTU17, which orchestrates the modulation of reactive oxygen species homeostasis to enhance salt tolerance.

Cotton sphingosine kinase GhLCBK1 participates in fiber cell elongation by affecting sphingosine-1-phophate and auxin synthesis

Sphingolipids serve as essential components of biomembrane and possess significant bioactive properties. Sphingosine-1-phophate (S1P) plays a key role in plant resistance to stress, but its specific impact on plant growth and development remains to be fully elucidated. Cotton fiber cells are an ideal material for investigating the growth and maturation of plant cells. In this study, we examined the content and composition of sphingosine (Sph) and S1P throughout the progression of fiber cell development. The content of S1P elevated gradually during fiber elongation but declined during the transition stage. Exogenous application of S1P promoted fiber elongation while using of FTY720 (an antagonist of S1P), and DMS (an inhibitor of LCBK) hindered fiber elongation. Cotton Long Chain Base Kinase 1 (GhLCBK1) was notably expressed during the fiber elongation stage, containing all conserved domains of LCBK protein and localized in the endoplasmic reticulum. Overexpression GhLCBK1 increased the S1P content and promoted fiber elongation while retarded secondary cell wall (SCW) deposition. Conversely, downregulation of GhLCBK1 reduced the S1P levels, and suppressed fiber elongation, and accelerated SCW deposition. Transcriptome analysis revealed that upregulating GhLCBK1 or applying S1P induced the expression of GhEXPANSIN and auxin related genes. Furthermore, the levels of IAA were elevated and reduced in the fibers when up-regulating or down-regulating GhLCBK1, respectively. Our investigation demonstrated that GhLCBK1 and its product S1P facilitated the elongation of fiber cells by affecting auxin biosynthesis. This study contributes novel insights into the intricate regulatory pathways involved in fiber cell elongation, identifying GhLCBK1 as a potential target gene and laying the groundwork for enhancing fiber quality via genetic manipulation.

Gibberellic acid promotes single-celled fiber elongation through the activation of two signaling cascades in cotton

The agricultural green revolution spectacularly enhanced crop yield through modification of gibberellin (GA) signaling. However, in cotton, the GA signaling cascades remain elusive, limiting our potential to cultivate new cotton varieties and improve yield and quality. Here, we identified that GA prominently stimulated fiber elongation through the degradation of DELLA protein GhSLR1, thereby disabling GhSLR1's physical interaction with two transcription factors, GhZFP8 and GhBLH1. Subsequently, the resultant free GhBLH1 binds to GhKCS12 promoter and activates its expression to enhance VLCFAs biosynthesis. With a similar mechanism, the free GhZFP8 binds to GhSDCP1 promoter and activates its expression. As a result, GhSDCP1 upregulates the expression of GhPIF3 gene associated with plant cell elongation. Ultimately, the two parallel signaling cascades synergistically promote cotton fiber elongation. Our findings outline the mechanistic framework that translates the GA signal into fiber cell elongation, thereby offering a roadmap to improve cotton fiber quality and yield.

A comprehensive review on Gossypium hirsutum resistance against cotton leaf curl virus

Cotton (Gossypium hirsutum L.) is a significant fiber crop. Being a major contributor to the textile industry requires continuous care and attention. Cotton is subjected to various biotic and abiotic constraints. Among these, biotic factors including cotton leaf curl virus (CLCuV) are dominant. CLCuV is a notorious disease of cotton and is acquired, carried, and transmitted by the whitefly (Bemisia tabaci). A cotton plant affected with CLCuV may show a wide range of symptoms such as yellowing of leaves, thickening of veins, upward or downward curling, formation of enations, and stunted growth. Though there are many efforts to protect the crop from CLCuV, long-term results are not yet obtained as CLCuV strains are capable of mutating and overcoming plant resistance. However, systemic-induced resistance using a gene-based approach remained effective until new virulent strains of CLCuV (like Cotton Leaf Curl Burewala Virus and others) came into existence. Disease control by biological means and the development of CLCuV-resistant cotton varieties are in progress. In this review, we first discussed in detail the evolution of cotton and CLCuV strains, the transmission mechanism of CLCuV, the genetic architecture of CLCuV vectors, and the use of pathogen and nonpathogen-based approaches to control CLCuD. Next, we delineate the uses of cutting-edge technologies like genome editing (with a special focus on CRISPR-Cas), next-generation technologies, and their application in cotton genomics and speed breeding to develop CLCuD resistant cotton germplasm in a short time. Finally, we delve into the current obstacles related to cotton genome editing and explore forthcoming pathways for enhancing precision in genome editing through the utilization of advanced genome editing technologies. These endeavors aim to enhance cotton's resilience against CLCuD.

GhMYB30-GhMUR3 affects fiber elongation and secondary wall thickening in cotton

Xyloglucan, an important hemicellulose, plays a crucial role in maintaining cell wall structure and cell elongation. However, the effects of xyloglucan on cotton fiber development are not well understood. GhMUR3 encodes a xyloglucan galactosyltransferase that is essential for xyloglucan synthesis and is highly expressed during fiber elongation. In this study, we report that GhMUR3 participates in cotton fiber development under the regulation of GhMYB30. Overexpression GhMUR3 affects the fiber elongation and cell wall thickening. Transcriptome showed that the expression of genes involved in secondary cell wall synthesis was prematurely activated in OE-MUR3 lines. In addition, GhMYB30 was identified as a key regulator of GhMUR3 by Y1H, Dual-Luc, and electrophoretic mobility shift assay (EMSA) assays. GhMYB30 directly bound the GhMUR3 promoter and activated GhMUR3 expression. Furthermore, DAP-seq of GhMYB30 was performed to identify its target genes in the whole genome. The results showed that many target genes were associated with fiber development, including cell wall synthesis-related genes, BR-related genes, reactive oxygen species pathway genes, and VLCFA synthesis genes. It was demonstrated that GhMYB30 may regulate fiber development through multiple pathways. Additionally, GhMYB46 was confirmed to be a target gene of GhMYB30 by EMSA, and GhMYB46 was significantly increased in GhMYB30-silenced lines, indicating that GhMYB30 inhibited GhMYB46 expression. Overall, these results revealed that GhMUR3 under the regulation of GhMYB30 and plays an essential role in cotton fiber elongation and secondary wall thickening. Additionally, GhMYB30 plays an important role in the regulation of fiber development and regulates fiber secondary wall synthesis by inhibiting the expression of GhMYB46.

Identification, evolution, and expression of GDSL-type Esterase/Lipase (GELP) gene family in three cotton species: a bioinformatic analysis

BACKGROUND

GDSL esterase/lipases (GELPs) play important roles in plant growth, development, and response to biotic and abiotic stresses. Presently, an extensive and in-depth analysis of GELP family genes in cotton is still not clear enough, which greatly limits the further understanding of cotton GELP function and regulatory mechanism.

RESULTS

A total of 389 GELP family genes were identified in three cotton species of Gossypium hirsutum (193), G. arboreum (97), and G. raimondii (99). These GELPs could be classified into three groups and eight subgroups, with the GELPs in same group to have similar gene structures and conserved motifs. Evolutionary event analysis showed that the GELP family genes tend to be diversified at the spatial dimension and certain conservative at the time dimension, with a trend of potential continuous expansion in the future. The orthologous or paralogous GELPs among different genomes/subgenomes indicated the inheritance from genome-wide duplication during polyploidization, and the paralogous GELPs were derived from chromosomal segment duplication or tandem replication. GELP genes in the A/D subgenome underwent at least three large-scale replication events in the evolutionary process during the period of 0.6-3.2 MYA, with two large-scale evolutionary events between 0.6-1.8 MYA that were associated with tetraploidization, and the large-scale duplication between 2.6-9.1 MYA that occurred during diploidization. The cotton GELPs indicated diverse expression patterns in tissue development, ovule and fiber growth, and in response to biotic and abiotic stresses, combining the existing cis-elements in the promoter regions, suggesting the GELPs involvements of functions to be diversification and of the mechanisms to be a hormone-mediated manner.

CONCLUSIONS

Our results provide a systematic and comprehensive understanding the function and regulatory mechanism of cotton GELP family, and offer an effective reference for in-depth genetic improvement utilization of cotton GELPs.

Single-cell RNA sequencing reveals a hierarchical transcriptional regulatory network of terpenoid biosynthesis in cotton secretory glandular cells

Plants can synthesize a wide range of terpenoids in response to various environmental cues. However, the specific regulatory mechanisms governing terpenoid biosynthesis at the cellular level remain largely elusive. In this study, we employed single-cell RNA sequencing to comprehensively characterize the transcriptome profile of cotton leaves and established a hierarchical transcriptional network regulating cell-specific terpenoid production. We observed substantial expression levels of genes associated with the biosynthesis of both volatile terpenes (such as β-caryophyllene and β-myrcene) and non-volatile gossypol-type terpenoids in secretory glandular cells. Moreover, two novel transcription factors, namely GoHSFA4a and GoNAC42, are identified to function downstream of the Gossypium PIGMENT GLAND FORMATION genes. Both transcription factors could directly regulate the expression of terpenoid biosynthetic genes in secretory glandular cells in response to developmental and environmental stimuli. For convenient retrieval of the single-cell RNA sequencing data generated in this study, we developed a user-friendly web server . Our findings not only offer valuable insights into the precise regulation of terpenoid biosynthesis genes in cotton leaves but also provide potential targets for cotton breeding endeavors.

Comprehensive identification and functional characterization of GhpPLA gene family in reproductive organ development

BACKGROUND

Phospholipases As (PLAs) are acyl hydrolases that catalyze the release of free fatty acids in phospholipids and play multiple functions in plant growth and development. The three families of PLAs are: PLA1, PLA2 (sPLA), and patatin-related PLA (pPLA). The diverse functions that pPLAs play in the growth and development of a broad range of plants have been demonstrated by prior studies.

METHODS

Genome-wide analysis of the pPLA gene family and screening of genes for expression verification and gene silencing verification were conducted. Additionally, pollen vitality testing, analysis of the pollen expression pattern, and the detection of POD, SOD, CAT, MDA, and H2O2 were performed.

RESULT

In this study, 294 pPLAs were identified from 13 plant species, including 46 GhpPLAs that were divided into three subfamilies (I-III). Expression patterns showed that the majority of GhpPLAs were preferentially expressed in the petal, pistil, anther, and ovule, among other reproductive organs. Particularly, GhpPLA23 and GhpPLA44, were found to be potentially important for the reproductive development of G. hirsutum. Functional validation was demonstrated by VIGS which showed that reduced expression levels of GhpPLA23 and GhpPLA44 in the silenced plants were associated with a decrease in pollen activity. Moreover, a substantial shift in ROS and ROS scavengers and a considerable increase in POD, CAT, SOD, and other physiological parameters was found out in these silenced plants. Our results provide plausibility to the hypothesis that GhpPLA23 and GhpPLA44 had a major developmental impact on cotton reproductive systems. These results also suggest that pPLAs are important for G. hirsutum's reproductive development and suggest that they could be employed as potential genes for haploid induction.

CONCLUSIONS

The findings of the present research indicate that pPLA genes are essential for the development of floral organs and sperm cells in cotton. Consequently, this family might be important for the reproductive development of cotton and possibly for inducing the plant develop haploid progeny.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1

BACKGROUND

Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear.

RESULTS

In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR).

CONCLUSIONS

Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

Systematically and Comprehensively Understanding the Regulation of Cotton Fiber Initiation: A Review

Cotton fibers provide an important source of raw materials for the textile industry worldwide. Cotton fiber is a kind of single cell that differentiates from the epidermis of the ovule and provides a perfect research model for the differentiation and elongation of plant cells. Cotton fiber initiation is the first stage throughout the entire developmental process. The number of fiber cell initials on the seed ovule epidermis decides the final fiber yield. Thus, it is of great significance to clarify the mechanism underlying cotton fiber initiation. Fiber cell initiation is controlled by complex and interrelated regulatory networks. Plant phytohormones, transcription factors, sugar signals, small signal molecules, functional genes, non-coding RNAs, and histone modification play important roles during this process. Here, we not only summarize the different kinds of factors involved in fiber cell initiation but also discuss the mechanisms of these factors that act together to regulate cotton fiber initiation. Our aim is to synthesize a systematic and comprehensive review of different factors during fiber initiation that will provide the basics for further illustrating these mechanisms and offer theoretical guidance for improving fiber yield in future molecular breeding work.

GhXB38D represses cotton fibre elongation through ubiquitination of ethylene biosynthesis enzymes GhACS4 and GhACO1

Ethylene plays an essential role in the development of cotton fibres. Ethylene biosynthesis in plants is elaborately regulated by the activities of key enzymes, 1-aminocyclopropane-1-carboxylate oxidase (ACO) and 1-aminocyclopropane-1-carboxylate synthase (ACS); however, the potential mechanism of post-translational modification of ACO and ACS to control ethylene synthesis in cotton fibres remains unclear. Here, we identify an E3 ubiquitin ligase, GhXB38D, that regulates ethylene biosynthesis during fibre elongation in cotton. GhXB38D gene is highly expressed in cotton fibres during the rapid elongation stage. Suppressing GhXB38D expression in cotton significantly enhanced fibre elongation and length, accompanied by the up-regulation of genes associated with ethylene signalling and fibre elongation. We demonstrated that GhXB38D interacts with the ethylene biosynthesis enzymes GhACS4 and GhACO1 in elongating fibres and specifically mediates their ubiquitination and degradation. The inhibition of GhXB38D gene expression increased the stability of GhACS4 and GhACO1 proteins in cotton fibres and ovules, resulting in an elevated concentration of ethylene. Our findings highlight the role of GhXB38D as a regulator of ethylene synthesis by ubiquitinating ACS4 and ACO1 proteins and modulating their stability. GhXB38D acts as a negative regulator of fibre elongation and serves as a potential target for enhancing cotton fibre yield and quality through gene editing strategy.

Genome-wide association study of fiber yield-related traits uncovers the novel genomic regions and candidate genes in Indian upland cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is a major fiber crop that is cultivated worldwide and has significant economic importance. India harbors the largest area for cotton cultivation, but its fiber yield is still compromised and ranks 22nd in terms of productivity. Genetic improvement of cotton fiber yield traits is one of the major goals of cotton breeding, but the understanding of the genetic architecture underlying cotton fiber yield traits remains limited and unclear. To better decipher the genetic variation associated with fiber yield traits, we conducted a comprehensive genome-wide association mapping study using 117 Indian cotton germplasm for six yield-related traits. To accomplish this, we generated 2,41,086 high-quality single nucleotide polymorphism (SNP) markers using genotyping-by-sequencing (GBS) methods. Population structure, PCA, kinship, and phylogenetic analyses divided the germplasm into two sub-populations, showing weak relatedness among the germplasms. Through association analysis, 205 SNPs and 134 QTLs were identified to be significantly associated with the six fiber yield traits. In total, 39 novel QTLs were identified in the current study, whereas 95 QTLs overlapped with existing public domain data in a comparative analysis. Eight QTLs, qGhBN_SCY_D6-1, qGhBN_SCY_D6-2, qGhBN_SCY_D6-3, qGhSI_LI_A5, qGhLI_SI_A13, qGhLI_SI_D9, qGhBW_SCY_A10, and qGhLP_BN_A8 were identified. Gene annotation of these fiber yield QTLs revealed 2,509 unique genes. These genes were predominantly enriched for different biological processes, such as plant cell wall synthesis, nutrient metabolism, and vegetative growth development in the gene ontology (GO) enrichment study. Furthermore, gene expression analysis using RNAseq data from 12 diverse cotton tissues identified 40 candidate genes (23 stable and 17 novel genes) to be transcriptionally active in different stages of fiber, ovule, and seed development. These findings have revealed a rich tapestry of genetic elements, including SNPs, QTLs, and candidate genes, and may have a high potential for improving fiber yield in future breeding programs for Indian 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.

Genetic and Morpho-Physiological Differences among Transgenic and No-Transgenic Cotton Cultivars

Three carbon-chain extension genes associated with fatty acid synthesis in upland cotton (Gossypium hirsutum), namely GhKAR, GhHAD, and GhENR, play important roles in oil accumulation in cotton seeds. In the present study, these three genes were cloned and characterized. The expression patterns of GhKAR, GhHAD, and GhENR in the high seed oil content cultivars 10H1014 and 10H1041 differed somewhat compared with those of 10H1007 and 2074B with low seed oil content at different stages of seed development. GhKAR showed all three cultivars showed higher transcript levels than that of 2074B at 10-, 40-, and 45-days post anthesis (DPA). The expression pattern of GhHAD showed a lower transcript level than that of 2074B at both 10 and 30 DPA but a higher transcript level than that of 2074B at 40 DPA. GhENR showed a lower transcript level than that of 2074B at both 15 and 30 DPA. The highest transcript levels of GhKAR and GhENR were detected at 15 DPA in 10H1007, 10H1014, and 10H1041 compared with 2074B. From 5 to 45 DPA cotton seed, the oil content accumulated continuously in the developing seed. Oil accumulation reached a peak between 40 DPA and 45 DPA and slightly decreased in mature seed. In addition, GhKAR and GhENR showed different expression patterns in fiber and ovule development processes, in which they showed high expression levels at 20 DPA during the fiber elongation stage, but their expression level peaked at 15 DPA during ovule development processes. These two genes showed the lowest expression levels at the late seed maturation stage, while GhHAD showed a peak of 10 DPA in fiber development. Compared to 2074B, the oil contents of GhKAR and GhENR overexpression lines increased 1.05~1.08 folds. These results indicated that GhHAD, GhENR, and GhKAR were involved in both seed oil synthesis and fiber elongation with dual biological functions in cotton.

Single-cell RNA landscape of the special fiber initiation process in Bombax ceiba

As a new source of natural fibers, the Bombax ceiba tree can provide thin, light, extremely soft and warm fiber material for the textile industry. Natural fibers are an ideal model system for studying cell growth and differentiation, but the molecular mechanisms that regulate fiber initiation are not fully understood. In B. ceiba, we found that fiber cells differentiate from the epidermis of the inner ovary wall. Each initiated cell then divides into a cluster of fiber cells that eventually develop into mature fibers, a process very different from the classical fiber initiation process of cotton. We used high-throughput single-cell RNA sequencing (scRNA-seq) to examine the special characteristics of fiber initiation in B. ceiba. A total of 15 567 high-quality cells were identified from the inner wall of the B. ceiba ovary, and 347 potential marker genes for fiber initiation cell types were identified. Two major cell types, initiated fiber cells and epidermal cells, were identified and verified by RNA in situ hybridization. A developmental trajectory analysis was used to reconstruct the process of fiber cell differentiation in B. ceiba. Comparative analysis of scRNA-seq data from B. ceiba and cotton (Gossypium hirsutum) confirmed that the additional cell division process in B. ceiba is a novel species-specific mechanism for fiber cell development. Candidate genes and key regulators that may contribute to fiber cell differentiation and division in B. ceiba were identified. This work reveals gene expression signatures during B. ceiba fiber initiation at a single-cell resolution, providing a new strategy and viewpoint for investigation of natural fiber cell differentiation and development.

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.

Large-scale long terminal repeat insertions produced a significant set of novel transcripts in cotton

Genomic analysis has revealed that the 1,637-Mb Gossypium arboreum genome contains approximately 81% transposable elements (TEs), while only 57% of the 735-Mb G. raimondii genome is occupied by TEs. In this study, we investigated whether there were unknown transcripts associated with TE or TE fragments and, if so, how these new transcripts were evolved and regulated. As sequence depths increased from 4 to 100 G, a total of 10,284 novel intergenic transcripts (intergenic genes) were discovered. On average, approximately 84% of these intergenic transcripts possibly overlapped with the long terminal repeat (LTR) insertions in the otherwise untranscribed intergenic regions and were expressed at relatively low levels. Most of these intergenic transcripts possessed no transcription activation markers, while the majority of the regular genic genes possessed at least one such marker. Genes without transcription activation markers formed their+1 and -1 nucleosomes more closely (only (117±1.4)bp apart), while twice as big spaces (approximately (403.5±46.0) bp apart) were detected for genes with the activation markers. The analysis of 183 previously assembled genomes across three different kingdoms demonstrated systematically that intergenic transcript numbers in a given genome correlated positively with its LTR content. Evolutionary analysis revealed that genic genes originated during one of the whole-genome duplication events around 137.7 million years ago (MYA) for all eudicot genomes or 13.7 MYA for the Gossypium family, respectively, while the intergenic transcripts evolved around 1.6 MYA, resultant of the last LTR insertion. The characterization of these low-transcribed intergenic transcripts can facilitate our understanding of the potential biological roles played by LTRs during speciation and diversifications.

Alternative Splicing during Fiber Development in G. hirsutum

Cotton is a valuable cash crop in many countries. Cotton fiber is a trichome that develops from a single epidermal cell and serves as an excellent model for understanding cell differentiation and other life processes. Alternative splicing (AS) of genes is a common post-transcriptional regulatory process in plants that is essential for plant growth and development. The process of AS during cotton fiber formation, on the other hand, is mainly unknown. A substantial number of multi-exon genes were discovered to be alternatively spliced during cotton fiber formation in this study, accounting for 23.31% of the total number of genes in Gossypium hirsutum. Retention intron (RI) is not necessarily the most common AS type, indicating that AS genes and processes during fiber development are very temporal and tissue-specific. When compared to fiber samples, AS is more prevalent at the fiber initiation stages and in the ovule, indicating that development stages and tissues use different AS strategies. Genes involved in fiber development have gone through stage-specific AS, demonstrating that AS regulates cotton fiber development. Furthermore, AS can be regulated by trans-regulation elements such as splicing factor and cis-regulation elements such as gene length, exon numbers, and GC content, particularly at exon-intron junction sites. Our findings also suggest that increased DNA methylation may aid in the efficiency of AS, and that gene body methylation is key in AS control. Finally, our research will provide useful information about the roles of AS during the cotton fiber development process.

Response of in situ root phenotypes to potassium stress in cotton

Potassium plays a significant role in the basic functions of plant growth and development. Potassium uptake is closely associated with morphological characteristics of the roots. However, the dynamic characteristics of phenotype and lifespan of cotton (Gossypium hirsutum L.) lateral roots and root hairs under low and high potassium stress remain unclear. In this study, potassium stress experiments (low and high potassium, medium potassium as control) were conducted using RhizoPot (an in situ root observation device) to determine the response characteristics of lateral roots and root hairs in cotton under potassium stress. The plant morphology, photosynthetic characteristics, root phenotypic changes, and lifespan of lateral roots and root hairs were measured. Potassium accumulation, aboveground phenotype, photosynthetic capacity, root length density, root dry weight, root diameter, lateral root lifespan, and root hair lifespan under low potassium stress were significantly decreased compared to medium potassium treatment. However, the root hair length of the former was significantly increased than that of the latter. Potassium accumulation and the lateral root lifespan were significantly increased under high potassium treatment, while root length density, root dry weight, root diameter, root hair length, and root hair lifespan were significantly decreased compared to the medium potassium treatment. Notably, there were no significant differences in aboveground morphology and photosynthetic characters. Principal component analysis revealed that lateral root lifespan, root hair lifespan of the first lateral root, and root hair length significantly correlated with potassium accumulation. The root had similar regularity responses to low and high potassium stress except for lifespan and root hair length. The findings of this study enhance the understanding of the phenotype and lifespan of cotton's lateral roots and root hairs under low and high potassium stress.

Genome-wide identification and expression analysis of terpene synthases in Gossypium species in response to gossypol biosynthesis

Cottonseed is an invaluable resource, providing protein, oil, and abundant minerals that significantly contribute to the well-being and nutritional needs of both humans and livestock. However, cottonseed also contains a toxic substance called gossypol, a secondary metabolite in Gossypium species that plays an important role in cotton plant development and self-protection. Herein, genome-wide analysis and characterization of the terpene synthase (TPS) gene family identified 304 TPS genes in Gossypium. Bioinformatics analysis revealed that the gene family was grouped into six subgroups TPS-a, TPS-b, TPS-c, TPS-e, TPS-f, and TPS-g. Whole-genome, segmental, and tandem duplication contributed to the evolution of TPS genes. According to the analysis of selection pressure, it was predicted that TPS genes experience predominantly negative selection, with positive selection occurring subsequently. RT-qPCR analysis in TM-1 and CRI-12 lines revealed GhTPS48 gene as the candidate gene for silencing experiments. To summarize, comprehensive genome-wide studies, RT-qPCR, and gene silencing experiments have collectively demonstrated the involvement of the TPS gene family in the biosynthesis of gossypol in cotton.

Comparative Transcriptome Analysis Revealed Key Genes Regulating Gossypol Synthesis in Tetraploid Cultivated Cotton

Tetraploid cultivated cotton (Gossypium spp.) produces cottonseeds rich in protein and oil. Gossypol and related terpenoids, stored in the pigment glands of cottonseeds, are toxic to human beings and monogastric animals. However, a comprehensive understanding of the genetic basis of gossypol and gland formation is still lacking. We performed a comprehensive transcriptome analysis of four glanded versus two glandless tetraploid cultivars distributed in Gossypium hirsutum and Gossypium barbadense. A weighted gene co-expression network analysis (WGCNA) based on 431 common differentially expressed genes (DEGs) uncovered a candidate module that was strongly associated with the reduction in or disappearance of gossypol and pigment glands. Further, the co-expression network helped us to focus on 29 hub genes, which played key roles in the regulation of related genes in the candidate module. The present study contributes to our understanding of the genetic basis of gossypol and gland formation and serves as a rich potential source for breeding cotton cultivars with gossypol-rich plants and gossypol-free cottonseed, which is beneficial for improving food safety, environmental protection, and economic gains of tetraploid cultivated cotton.

Brassinosteroid signaling regulates phosphate starvation-induced malate secretion in plants

Inorganic phosphate (Pi) is often limited in soils due to precipitation with iron (Fe) and aluminum (Al). To scavenge heterogeneously distributed phosphorus (P) resources, plants have evolved a local Pi signaling pathway that induces malate secretion to solubilize the occluded Fe-P or Al-P oxides. In this study, we show that Pi limitation impaired brassinosteroid signaling and downregulated BRASSINAZOLE-RESISTANT 1 (BZR1) expression in Arabidopsis thaliana. Exogenous 2,4-epibrassinolide treatment or constitutive activation of BZR1 (in the bzr1-D mutant) significantly reduced primary root growth inhibition under Pi-starvation conditions by downregulating ALUMINUM-ACTIVATED MALATE TRANSPORTER 1 (ALMT1) expression and malate secretion. Furthermore, AtBZR1 competitively suppressed the activator effect of SENSITIVITY TO PROTON RHIZOTOXICITY 1 (STOP1) on ALMT1 expression and malate secretion in Nicotiana benthamiana leaves and Arabidopsis. The ratio of nuclear-localized STOP1 and BZR1 determined ALMT1 expression and malate secretion in Arabidopsis. In addition, BZR1-inhibited malate secretion is conserved in rice (Oryza sativa). Our findings provide insight into plant mechanisms for optimizing the secretion of malate, an important carbon resource, to adapt to Pi-deficiency stress.

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.

Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of pigment glands in Gossypium bickii

Comprehensive utilization of cottonseeds is limited by the presence of pigment glands and its inclusion gossypol. The ideal cotton has glandless seeds but a glanded plant, a trait found in only a few Australian wild cotton species, including Gossypium bickii. Introgression of this trait into cultivated species has proved to be difficult. Understanding the biological processes toward pigment gland morphogenesis and the associated underlying molecular mechanisms will facilitate breeding of cultivated cotton varieties with the trait of glandless seeds and glanded plant. In this study, single-cell RNA sequencing (scRNA-seq) was performed on 12 222 protoplasts isolated from cotyledons of germinating G. bickii seeds 48 h after imbibition. Clustered into 14 distinct clusters unsupervisedly, these cells could be grouped into eight cell populations with the assistance of known cell marker genes. The pigment gland cells were well separated from others and could be separated into pigment gland parenchyma cells, secretory cells, and apoptotic cells. By integrating the pigment gland cell developmental trajectory, transcription factor regulatory networks, and core transcription factor functional validation, we established a model for pigment gland formation. In this model, light and gibberellin were verified to promote the formation of pigment glands. In addition, three novel genes, GbiERF114 (ETHYLENE RESPONSE FACTOR 114), GbiZAT11 (ZINC FINGER OF ARABIDOPSIS THALIANA 11), and GbiNTL9 (NAC TRANSCRIPTION FACTOR-LIKE 9), were found to affect pigment gland formation. Collectively, these findings provide new insights into pigment gland morphogenesis and lay the cornerstone for future cotton scRNA-seq investigations.

Dirigent gene editing of gossypol enantiomers for toxicity-depleted cotton seeds

Axial chirality of biaryls can generate varied bioactivities. Gossypol is a binaphthyl compound made by cotton plants. Of its two axially chiral isomers, (-)-gossypol is the bioactive form in mammals and has antispermatogenic activity, and its accumulation in cotton seeds poses health concerns. Here we identified two extracellular dirigent proteins (DIRs) from Gossypium hirsutum, GhDIR5 and GhDIR6, which impart the hemigossypol oxidative coupling into (-)- and (+)-gossypol, respectively. To reduce cotton seed toxicity, we disrupted GhDIR5 by genome editing, which eliminated (-)-gossypol but had no effects on other phytoalexins, including (+)-gossypol, that provide pest resistance. Reciprocal mutagenesis identified three residues responsible for enantioselectivity. The (-)-gossypol-forming DIRs emerged later than their enantiocomplementary counterparts, from tandem gene duplications that occurred shortly after the cotton genus diverged. Our study offers insight into how plants control enantiomeric ratios and how to selectively modify the chemical spectra of cotton plants and thereby improve crop quality.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum

The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC₄F_3:5). A single segment substitution line (CSSL-106) screened from BC₆F₂ was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC₇F₂ individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

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.

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.

Genome-wide identification and expression analysis of the cotton patatin-related phospholipase A genes and response to stress tolerance

Patatin-related phospholipase A genes were involved in the response of Gossypium hirsutum to drought and salt tolerance. pPLA (patatin-related phospholipase A) is a key enzyme that catalyzes the initial step of lipid hydrolysis, which is involved in biological processes, such as drought, salt stress, and freezing injury. However, a comprehensive analysis of the pPLA gene family in cotton, especially the role of pPLA in the response to drought and salt tolerance, has not been reported so far. A total of 33 pPLA genes were identified in this study using a genome-wide search approach, and phylogenetic analysis classified these genes into three groups. These genes are unevenly distributed on the 26 chromosomes of cotton, and most of them contain a few introns. The results of the collinear analysis showed that G. hirsutum contained 1-5 copies of each pPLA gene found in G. arboreum and G. raimondii. The subcellular localization analysis of Gh_D08G061200 showed that the protein was localized in the nucleus. In addition, analysis of published upland cotton transcriptome data revealed that six GhPLA genes were expressed in various tissues and organs. Two genes (Gh_A04G142100.1 and Gh_D04G181000.1) were highly expressed in all tissues under normal conditions, showing the expression characteristics of housekeeping genes. Under different drought and salt tolerance stresses, we detected four genes with different expression levels. This study helps to clarify the role of pPLA in the response to drought and salt tolerance.

Genomic identification of cotton SAC genes branded ovule and stress-related key genes in Gossypium hirsutum

SAC genes have been identified to play a variety of biological functions and responses to various stresses. Previously, SAC genes have been recognized in animals and Arabidopsis. For the very first time, we identified 157 SAC genes in eight cotton species including three diploids and five tetraploids with 23 SAC members in G. hirsutum. Evolutionary analysis classified all cotton SAC gene family members into five distinct groups. Cotton SAC genes showed conserved sequence logos and WGD or segmental duplication. Multiple synteny and collinearity analyses revealed gene family expansion and purifying selection pressure during evolution. G. hirsutum SAC genes showed uneven chromosomal distribution, multiple exons/introns, conserved protein motifs, and various growth and stress-related cis-elements. Expression pattern analysis revealed three GhSAC genes (GhSAC3, GhSAC14, and GhSAC20) preferentially expressed in flower, five genes (GhSAC1, GhSAC6, GhSAC9, GhSAC13, and GhSAC18) preferentially expressed in ovule and one gene (GhSAC5) preferentially expressed in fiber. Similarly, abiotic stress treatment verified that GhSAC5 was downregulated under all stresses, GhSAC6 and GhSAC9 were upregulated under NaCl treatment, and GhSAC9 and GhSAC18 were upregulated under PEG and heat treatment respectively. Overall, this study identified key genes related to flower, ovule, and fiber development and important genetic material for breeding cotton under abiotic stress conditions.

Biotechnological detoxification: an unchanging source-sink balance strategy for crop improvement

The wide occurrence of natural phytotoxins renders many crops unfit for human consumption. To overcome this problem and produce detoxified crop varieties, we propose the use of biotechnological strategies that can enhance the harvest index without the need to increase crop biomass or alter whole plant architecture.