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

Genetic linkage analysis of stable QTLs in Gossypium hirsutum RIL population revealed function of GhCesA4 in fiber development

INTRODUCTION

Upland cotton is an important allotetrapolyploid crop providing natural fibers for textile industry. Under the present high-level breeding and production conditions, further simultaneous improvement of fiber quality and yield is facing unprecedented challenges due to their complex negative correlations.

OBJECTIVES

The study was to adequately identify quantitative trait loci (QTLs) and dissect how they orchestrate the formation of fiber quality and yield.

METHODS

A high-density genetic map (HDGM) based on an intraspecific recombinant inbred line (RIL) population consisting of 231 individuals was used to identify QTLs and QTL clusters of fiber quality and yield traits. The weighted gene correlation network analysis (WGCNA) package in R software was utilized to identify WGCNA network and hub genes related to fiber development. Gene functions were verified via virus-induced gene silencing (VIGS) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 strategies.

RESULTS

An HDGM consisting of 8045 markers was constructed spanning 4943.01 cM of cotton genome. A total of 295 QTLs were identified based on multi-environmental phenotypes. Among 139 stable QTLs, including 35 newly identified ones, seventy five were of fiber quality and 64 yield traits. A total of 33 QTL clusters harboring 74 QTLs were identified. Eleven candidate hub genes were identified via WGCNA using genes in all stable QTLs and QTL clusters. The relative expression profiles of these hub genes revealed their correlations with fiber development. VIGS and CRISPR/Cas9 edition revealed that the hub gene cellulose synthase 4 (GhCesA4, GH_D07G2262) positively regulate fiber length and fiber strength formation and negatively lint percentage.

CONCLUSION

Multiple analyses demonstrate that the hub genes harbored in the QTLs orchestrate the fiber development. The hub gene GhCesA4 has opposite pleiotropic effects in regulating trait formation of fiber quality and yield. The results facilitate understanding the genetic basis of negative correlation between cotton fiber quality and yield.

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

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

Membrane Interactions of GET1 and GET2 Facilitate Fiber Cell Initiation through the Guided Entry of the TA Protein Pathway in Cotton

The guided entry of TA proteins (GET) pathway, which is responsible for the post-translational targeting and insertion of the tail-anchored (TA) protein into the endoplasmic reticulum (ER), plays an important role in physiological processes such as protein sorting, vesicle trafficking, cell apoptosis, and enzymatic reactions in which the GET1/2 complex is indispensable. However, a comprehensive study of the GET1 and GET2 genes and the GET pathway in cotton has not yet been carried out. Here, 12 GET1 and 21 GET2 genes were identified in nine representative plant species, and the phylogenetic relationships, gene structures, protein motifs, cis-regulatory elements (CREs), and temporal and spatial expression profiles were analyzed thoroughly. Our study indicated that GhGET1s and GhGET2s might be localized on ER membranes. According to expression profiling and CREs analysis, GhGET2-A02 was identified as a promising candidate for fiber cell development, interacting with two GhGET1s in the membrane, with a binding bias toward GhGET1-A06. Silencing of GhGET1-A06 or GhGET2-A02 reduced fiber initiation and elongation. In summary, our research provides important evidence for understanding the gene families and functions of GET1 and GET2 in cotton and provides clues for molecular breeding of high-quality cotton fiber varieties.

GhEXL3 participates in brassinosteroids regulation of fiber elongation in Gossypium hirsutum

Cotton fiber (Gossypium hirsutum) serves as an ideal model for investigating the molecular mechanisms of plant cell elongation at the single-cell level. Brassinosteroids (BRs) play a crucial role in regulating plant growth and development. However, the mechanism by which BR influences cotton fiber elongation remains incompletely understood. In this study, we identified EXORDIUM-like (GhEXL3) through transcriptome analysis of fibers from BR-deficient cotton mutant pagoda 1 (pag1) and BRI1-EMS-SUPPRESSOR 1 (GhBES1.4, encoding a central transcription factor of BR signaling) overexpression cotton lines. Knockout of GhEXL3 using CRISPR/Cas9 was found to impede cotton fiber elongation, while its overexpression promoted fiber elongation, suggesting a positive regulatory function for GhEXL3 in fiber elongation. Furthermore, in vitro ovule culture experiments revealed that the overexpression of GhEXL3 partially counteracted the inhibitory effects of brassinazole (BRZ) on cotton fiber elongation, providing additional evidence of GhEXL3 involvement in BR signaling pathways. Moreover, our findings demonstrate that GhBES1.4 directly binds to the E-box (CACGTG) motif in the GhEXL3 promoter region and enhances its transcription. RNA-seq analysis revealed that overexpression of GhEXL3 upregulated the expression of EXPs, XTHs, and other genes associated with fiber cell elongation. Overall, our study contributes to understanding the mechanism by which BR regulates the elongation of cotton fibers through the direct modulation of GhEXL3 expression by GhBES1.4.

Genome‑wide analysis of cotton SCAMP genes and functional characterization of GhSCAMP2 and GhSCAMP4 in salt tolerance

BACKGROUND

Secretory carrier membrane proteins (SCAMPs) form a family of integral membrane proteins and play a crucial role in mediating exocytosis in both animals and plants. While SCAMP genes have been studied in several plant species, their functions in cotton, particularly in response to abiotic stress, have not yet been reported.

RESULTS

In this study, a total of 53 SCAMP genes were identified in G. arboreum, G. raimondii, G. hirsutum, and G. barbadense. These genes were classified into five groups based on a phylogenetic analysis with SCAMPs from Arabidopsis thaliana. The main factor driving the expansion of the SCAMP gene family in G. hirsutum is tandem and segmental duplication events. Using MEME, in addition to the conserved SCAMP domain, we identified 3-13 other domains in each GhSCAMP. The cis-element analysis suggested that GhSCAMPs were widely involved in cotton growth and development, and responses to abiotic stresses. RNA sequencing (RNA-Seq) and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) results showed that most GhSCAMPs were expressed highly in many tissues and had differential expression responses to drought, cold, and heat stresses. Knock-down of GhSCAMP2 and GhSCAMP4 by virus-induced gene silencing (VIGS) lead to a salt-sensitive phenotype and had a lower content of CAT, POD, and SOD.

CONCLUSIONS

This study identified SCAMP genes in four cotton species, enhancing our understanding of the potential biological functions of SCAMPs. Additionally, we demonstrated that GhSCAMP2 and GhSCAMP4 positively regulate cotton tolerance to salt stress.

Regulatory networks of coresident subgenomes during rapid fiber cell elongation in upland cotton

Cotton, an intriguing plant species shaped by polyploidization, evolution, and domestication, holds particular interest due to the complex mechanisms governing fiber traits across its two subgenomes. However, the regulatory elements or transcriptional networks between subgenomes during fiber elongation remain to be fully clarified. Here, we analyzed 1462 cotton fiber samples to reconstruct the gene-expression regulatory networks that influence fiber cell elongation. Inter-subgenome expression quantitative trait loci (eQTLs) largely dictate gene transcription, with a notable tendency for the D subgenome to regulate A-subgenome eGenes. This regulation reveals synchronized homoeologous gene expression driven by co-localized eQTLs and divergent patterns that diminish genetic correlations, thus leading to preferential expression in the A and D subgenomes. Hotspot456 emerged as a key regulator of fiber initiation and elongation, and artificial selection of trans-eQTLs in hotspot456 that positively regulate KCS1 has facilitated cell elongation. Experiments designed to clarify the roles of trans-eQTLs in improved fiber breeding confirmed the inhibition of GhTOL9 by a specific trans-eQTL via GhWRKY28, which negatively affects fiber elongation. We propose a model in which the GhWRKY28-GhTOL9 module regulates this process through the ESCRT (endosomal sorting complex required for transport) pathway. This research significantly advances our understanding of cotton's evolutionary and domestication processes and the intricate regulatory mechanisms that underlie significant plant traits.

Strigolactone-induced degradation of SMXL7 and SMXL8 contributes to gibberellin- and auxin-mediated fiber cell elongation in cotton

Cotton (Gossypium) fiber length, a key trait determining fiber yield and quality, is highly regulated by a class of recently identified phytohormones, strigolactones (SLs). However, the underlying molecular mechanisms of SL signaling involved in fiber cell development are largely unknown. Here, we show that the SL signaling repressors MORE AXILLARY GROWTH2-LIKE7 (GhSMXL7) and GhSMXL8 negatively regulate cotton fiber elongation. Specifically, GhSMXL7 and GhSMXL8 inhibit the polyubiquitination and degradation of the gibberellin (GA)-triggered DELLA protein (GhSLR1). Biochemical analysis revealed that GhSMXL7 and GhSMXL8 physically interact with GhSLR1, which interferes with the association of GhSLR1 with the E3 ligase GA INSENSITIVE2 (GhGID2), leading to the repression of GA signal transduction. GhSMXL7 also interacts with the transcription factor GhHOX3, preventing its binding to the promoters of essential fiber elongation regulatory genes. Moreover, both GhSMXL7 and GhSMXL8 directly bind to the promoter regions of the AUXIN RESPONSE FACTOR (ARF) genes GhARF18-10A, GhARF18-10D, and GhARF19-7D to suppress their expression. Cotton plants in which GhARF18-10A, GhARF18-10D, and GhARF19-7D transcript levels had been reduced by virus-induced gene silencing (VIGS) displayed reduced fiber length compared with control plants. Collectively, our findings reveal a mechanism illustrating how SL integrates GA and auxin signaling to coordinately regulate plant cell elongation at the single-cell level.

R2R3 MYB Transcription Factor GhMYB201 Promotes Cotton Fiber Elongation via Cell Wall Loosening and Very-Long-Chain Fatty Acid Synthesis

Cotton fiber is the leading natural textile material, and fiber elongation plays an essential role in the formation of cotton yield and quality. Although a number of components in the molecular network controlling cotton fiber elongation have been reported, a lot of players still need to be functionally dissected to understand the regulatory mechanism of fiber elongation comprehensively. In the present study, an R2R3-MYB transcription factor gene, GhMYB201, was characterized and functionally verified via CRISPR/Cas9-mediated gene editing. GhMYB201 was homologous to Arabidopsis AtMYB60, and both coding genes (GhMYB201At and GhMYB201Dt) were preferentially expressed in elongating cotton fibers. Knocking-out of GhMYB201 significantly reduced the rate and duration of fiber elongation, resulting in shorter and coarser mature fibers. It was found that GhMYB201 could bind and activate the transcription of cell wall loosening genes (GhRDLs) and also β-ketoacyl-CoA synthase genes (GhKCSs) to enhance very-long-chain fatty acid (VLCFA) levels in elongating fibers. Taken together, our data demonstrated that the transcription factor GhMYB201s plays an essential role in promoting fiber elongation via activating genes related to cell wall loosening and VLCFA biosynthesis.

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.

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

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

Deciphering the dynamic expression network of fiber elongation and the functional role of the GhTUB5 gene for fiber length in cotton based on an introgression population of upland cotton

INTRODUCTION

Interspecific introgression between Gossypium hirsutum and G. barbadense allows breeding cotton with outstanding fiber length (FL). However, the dynamic gene regulatory network of FL-related genes has not been characterized, and the functional mechanism through which the hub gene GhTUB5 mediates fiber elongation has yet to be determined.

METHODS

Coexpression analyses of 277 developing fiber transcriptomes integrated with QTL mapping using 250 introgression lines of different FL phenotypes were conducted to identify genes related to fiber elongation. The function of GhTUB5 was determined by ectopic expression of two TUB5 alleles in Arabidopsis and knockout of GhTUB5 in upland cotton. Yeast two-hybrid, split-luciferase and pull-down assays were conducted to screen for interacting proteins, and upstream genes were identified by yeast one-hybrid, dual-LUC and electrophoretic mobility shift assays.

RESULTS

The 32,612, 30,837 and 30,277 genes expressed at 5, 10 and 15 days postanthesis (dpa) were grouped into 19 distinct coexpression modules, and 988 genes in the MEblack module were enriched in the cell wall process and exhibited significant associations with FL. A total of 20 FL-QTLs were identified, each explaining 3.34-16.04 % of the phenotypic variance in the FL. Furthermore, several FL-QTLs contained 15 genes that were differentially expressed in the MEblack module including the tubulin beta gene (TUB5). Compared with the wild type, the overexpression of GhTUB5 and GbTUB5 in Arabidopsis suppressed root cell length but promoted cellulose synthesis. Knockout of GhTUB5 resulted in longer fiber lines. Protein-based experiments revealed that GhTUB5 interacts with GhZFP6. Additionally, GhTUB5 was directly activated by GhHD-ZIP7, a homeobox-leucine zipper transcription factor, and its paralogous gene was previously reported to mediate fiber elongation.

CONCLUSION

This study opens a new avenue to dissect functional mechanism of cotton fiber elongation. Our findings provide some molecular details on how GhTUB5 mediates the FL phenotype in cotton.

A fatty acid elongase complex regulates cell membrane integrity and septin-dependent host infection by the rice blast fungus

Very-long-chain fatty acids (VLCFAs) regulate biophysical properties of cell membranes to determine growth and development of eukaryotes, such as the pathogenesis of the rice blast fungus Magnaporthe oryzae. The fatty acid elongase Elo1 regulates pathogenesis of M. oryzae by modulating VLCFA biosynthesis. However, it remains unknown whether and how Elo1 associates with other factors to regulate VLCFA biosynthesis in fungal pathogens. Here, we identified Ifa38, Phs1 and Tsc13 as interacting proteins of Elo1 by proximity labelling in M. oryzae. Elo1 associated with Ifa38, Phs1 and Tsc13 on the endoplasmic reticulum (ER) membrane to control VLCFA biosynthesis. Targeted gene deletion mutants Δifa38, Δphs1 and Δtsc13 were all similarly impaired as Δelo1 in vegetative growth, conidial morphology, stress responses in ER, cell wall and membrane. These deletion mutants also displayed severe damage in cell membrane integrity and failed to organize the septin ring that is essential for penetration peg formation and pathogenicity. Our study demonstrates that M. oryzae employs a fatty acid elongase complex to regulate VLCFAs for maintaining or remodelling cell membrane structure, which is important for septin-mediated host penetration.

Integrative Transcriptomics and Proteomics Analysis of a Cotton Mutant yl1 with a Chlorophyll-Reduced Leaf

Leaf color mutants serve as ideal materials for studying photosynthesis, chlorophyll metabolism, and other physiological processes. Here, we identified a spontaneous yellow-leaf mutant (yl1) with chlorophyll-reduced leaves from G. hirsutum L. cv ZM24. Compare to wild type ZM24 with green leaves, yl1 exhibited patchy yellow leaves and reduced chlorophyll content. To further explore the mechanisms of the patchy yellow phenotype of the mutant plant, the transcriptomics and proteomics profiles were conducted for the mutant and wild types. A total of 9247 differentially expressed genes (DEGs) and 1368 differentially accumulated proteins (DAPs) were identified. Following gene ontology (GO) annotation and KEGG enrichment, the DEGs/DAPs were found to be significantly involved in multiple important pathways, including the obsolete oxidation-reduction process, photosynthesis, light-harvesting, the microtubule-based process, cell redox homeostasis, and the carbohydrate metabolic process. In photosynthesis and the light-harvesting pathway, a total of 39 DAPs/DEGs were identified, including 9 genes in the PSI, 7 genes in the PS II, 9 genes in the light-harvesting chlorophyll protein complex (LHC), 10 genes in the PsbP family, and 4 genes in the cytochrome b6/f complex. To validate the reliability of the omics data, GhPPD1, a DAPs in the PsbP family, was knocked down in cotton using the TRV-based VIGS system, and it was observed that the GhPPD1-silenced plants exhibited patchy yellow color, accompanied by a significant decrease in chlorophyll content. In conclusion, this study integrated transcriptomic and proteomic approaches to gain a deeper understanding of the mechanisms underlying the chlorophyll-reduced leaf phenotype.

Single-cell transcriptome atlas reveals somatic cell embryogenic differentiation features during regeneration

Understanding somatic cell totipotency remains a challenge facing scientific inquiry today. Plants display remarkable cell totipotency expression, illustrated by single-cell differentiation during somatic embryogenesis (SE) for plant regeneration. Determining cell identity and exploring gene regulation in such complex heterogeneous somatic cell differentiation have been major challenges. Here, we performed high-throughput single-cell sequencing assays to define the precise cellular landscape and revealed the modulation mode of marker genes during embryogenic differentiation in cotton (Gossypium hirsutum L.) as the crop for biotechnology application. We demonstrated that nonembryogenic calli (NEC) and primary embryogenic calli (PEC) tissues were composed of heterogeneous cells that could be partitioned into four broad populations with six distinct cell clusters. Enriched cell clusters and cell states were identified in NEC and PEC samples, respectively. Moreover, a broad repertoire of new cluster-specific genes and associated expression modules were identified. The energy metabolism, signal transduction, environmental adaptation, membrane transport pathways, and a series of transcription factors were preferentially enriched in cell embryogenic totipotency expression. Notably, the SE-ASSOCIATED LIPID TRANSFER PROTEIN (SELTP) gene dose-dependently marked cell types with distinct embryogenic states and exhibited a parabolic curve pattern along the somatic cell embryogenic differentiation trajectory, suggesting that SELTP could serve as a favorable quantitative cellular marker for detecting embryogenic expression at the single-cell level. In addition, RNA velocity and Scissor analysis confirmed the pseudo-temporal model and validated the accuracy of the scRNA-seq data, respectively. This work provides valuable marker-genes resources and defines precise cellular taxonomy and trajectory atlases for somatic cell embryogenic differentiation in plant regeneration.

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.

Cotton transposon-related variome reveals roles of transposon-related variations in modern cotton cultivation

INTRODUCTION

Transposon plays a vital role in cotton genome evolution, contributing to the expansion and divergence of genomes within the Gossypium genus. However, knowledge of transposon activity in modern cotton cultivation is limited.

OBJECTIVES

In this study, we aimed to construct transposon-related variome within Gossypium genus and reveal role of transposon-related variations during cotton cultivation. In addition, we try to identify valuable transposon-related variations for cotton breeding.

METHODS

We utilized graphical genome construction to build up the graphical transposon-related variome. Based on the graphical variome, we integrated t-test, eQTL analysis and Mendelian Randomization (MR) to identify valuable transposon activities and elite genes. In addition, a convolutional neural network (CNN) model was constructed to evaluate epigenomic effects of transposon-related variations.

RESULTS

We identified 35,980 transposon activities among 10 cotton genomes, and the diversity of genomic and epigenomic features was observed among 21 transposon categories. The graphical cotton transposon-related variome was constructed, and 9,614 transposon-related variations with plasticity in the modern cotton cohort were used for eQTL, phenotypic t-test and Mendelian Randomization. 128 genes were identified as gene resources improving fiber length and strength simultaneously. 4 genes were selected from 128 genes to construct the elite gene panel whose utility has been validated in a natural cotton cohort and 2 accessions with phenotypic divergence. Based on the eQTL analysis results, we identified transposon-related variations involved in cotton's environmental adaption and human domestication, providing evidence of their role in cotton's adaption-domestication cooperation.

CONCLUSIONS

The cotton transposon-related variome revealed the role of transposon-related variations in modern cotton cultivation, providing genomic resources for cotton molecular breeding.

Genomic Analysis of Brassinosteroid Biosynthesis Gene Family Reveals Its Roles in Cotton Development across Gossypium Species

Cotton is a globally significant economic crop. Brassinosteroids (BRs) are crucial to cotton development. This study systematically analyzed the BR synthase gene family in four cotton species and identified 60 BR genes: 20 in Gossypium hirsutum (GhBRs), 20 in G. barbadense (GbBRs), 10 in G. arboreum (GaBRs), and 10 in G. raimondii (GrBRs). The analysis was extended to chromosomal localization, evolutionary relationships, domain features, and cis-regulatory elements in the promoter regions of BR synthase genes. The results showed that the BR synthase genes were evenly distributed across different subgenomes and chromosomes. Bioinformatics analyses revealed high conservation of amino acid sequences, secondary structures, and conserved domains among the subfamily members, which is closely linked to their pivotal roles in the BR biosynthesis pathway. Cis-element distribution analysis of the BR synthase genes further underscored the complexity of BR gene expression regulation, which is influenced by multiple factors, including plant hormones, abiotic stress, and transcription factors. Expression profiling of GhBRs genes in various cotton tissues and developmental stages highlighted the key roles of GhROT3-1 and GhDET2-1 in fiber elongation and initiation, respectively. Protein-protein interactions and transcription factor analyses further elucidated the regulatory mechanisms of GhROT3-1 and GhDET2-1 in cotton growth and development. This study lays a theoretical foundation for understanding the role of the BR signaling pathway in cotton development, facilitating molecular breeding.

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

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.

Organelle Ca2+/CAM1-SELTP confers somatic cell embryogenic competence acquisition and transformation in plant regeneration

Somatic cell totipotency in plant regeneration represents the forefront of the compelling scientific puzzles and one of the most challenging problems in biology. How somatic embryogenic competence is achieved in regeneration remains elusive. Here, we discover uncharacterized organelle-based embryogenic differentiation processes of intracellular acquisition and intercellular transformation, and demonstrate the underlying regulatory system of somatic embryogenesis-associated lipid transfer protein (SELTP) and its interactor calmodulin1 (CAM1) in cotton as the pioneer crop for biotechnology application. The synergistic CAM1 and SELTP exhibit consistent dynamical amyloplast-plasmodesmata (PD) localization patterns but show opposite functional effects. CAM1 inhibits the effect of SELTP to regulate embryogenic differentiation for plant regeneration. It is noteworthy that callus grafting assay reflects intercellular trafficking of CAM1 through PD for embryogenic transformation. This work originally provides insight into the mechanisms responsible for embryogenic competence acquisition and transformation mediated by the Ca2+/CAM1-SELTP regulatory pathway, suggesting a principle for plant regeneration and cell/genetic engineering.

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

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Transcriptional and translational landscape fine-tune genome annotation and explores translation control in cotton

INTRODUCTION

The unavailability of intergenic region annotation in whole genome sequencing and pan-genomics hinders efforts to enhance crop improvement.

OBJECTIVES

Despite advances in research, the impact of post-transcriptional regulation on fiber development and translatome profiling at different stages of fiber growth in cotton (G. hirsutum) remains unexplored.

METHODS

We utilized a combination of reference-guided de novo transcriptome assembly and ribosome profiling techniques to uncover the hidden mechanisms of translational control in eight distinct tissues of upland cotton.

RESULTS

Our study identified P-site distribution at three-nucleotide periodicity and dominant ribosome footprint at 27 nucleotides. Specifically, we have detected 1,589 small open reading frames (sORFs), including 1,376 upstream ORFs (uORFs) and 213 downstream ORFs (dORFs), as well as 552 long non-coding RNAs (lncRNAs) with potential coding functions, which fine-tune the annotation of the cotton genome. Further, we have identified novel genes and lncRNAs with strong translation efficiency (TE), while sORFs were found to affect mRNA transcription levels during fiber elongation. The reliability of these findings was confirmed by the high consistency in correlation and synergetic fold change between RNA-sequencing (RNA-seq) and Ribosome-sequencing (Ribo-seq) analyses. Additionally, integrated omics analysis of the normal fiber ZM24 and short fiber pag1 cotton mutant revealed several differentially expressed genes (DEGs), and fiber-specific expressed (high/low) genes associated with sORFs (uORFs and dORFs). These findings were further supported by the overexpression and knockdown of GhKCS6, a gene associated with sORFs in cotton, and demonstrated the potential regulation of the mechanism governing fiber elongation on both the transcriptional and post-transcriptional levels.

CONCLUSION

Reference-guided transcriptome assembly and the identification of novel transcripts fine-tune the annotation of the cotton genome and predicted the landscape of fiber development. Our approach provided a high-throughput method, based on multi-omics, for discovering unannotated ORFs, hidden translational control, and complex regulatory mechanisms in crop plants.

Zinc Finger Protein8 (GhZFP8) Regulates the Initiation of Trichomes in Arabidopsis and the Development of Fiber in Cotton

Cotton is one of the most important natural fibers used in the textile industry worldwide. It is important to identify the key factors involved in cotton fiber development. In this study, zinc finger protein8 (GhZFP8) encoding a C2H2 transcription factor (TF) was cloned from cotton. qPCR showed that the transcripts of GhZFP8 in cotton were detected in the leaves and fibers at 3, 6, and 30 days post-anthesis (DPA), but not in the roots, stems, or flowers. The overexpression of GhZFP8 increased the trichome number on the siliques, leaves, and inflorescence, but inhibited the growth. The expression of trichome development and cell-elongation-related genes decreased obviously in GhZFP8 overexpressor Arabidopsis. Indole-3-acetic acid (IAA) and 1-Aminocyclopropanecarboxylic acid (ACC) contents were much higher in GhZFP8 overexpressors than that found in the wild type, but the gibberellin (GA) content was lower. The interference of GhZFP8 in cotton caused smaller bolls and shorter fibers than that of the control. The results of DNA affinity purification (DAP)-seq showed that GhZFP8 could bind to the promoter, exon, intron, and intergenic region of the target genes, which are involved in photosynthesis, signal transduction, synthesis of biomass, etc. Our findings implied that GhZFP8 processed multiple biological functions and regulated the development of cotton fiber.

Isolation and Functional Characterization of a Constitutive Promoter in Upland Cotton (Gossypium hirsutum L.)

Multiple cis-acting elements are present in promoter sequences that play critical regulatory roles in gene transcription and expression. In this study, we isolated the cotton FDH (Fiddlehead) gene promoter (pGhFDH) using a real-time reverse transcription-PCR (qRT-PCR) expression analysis and performed a cis-acting elements prediction analysis. The plant expression vector pGhFDH::GUS was constructed using the Gateway approach and was used for the genetic transformation of Arabidopsis and upland cotton plants to obtain transgenic lines. Histochemical staining and a β-glucuronidase (GUS) activity assay showed that the GUS protein was detected in the roots, stems, leaves, inflorescences, and pods of transgenic Arabidopsis thaliana lines. Notably, high GUS activity was observed in different tissues. In the transgenic lines, high GUS activity was detected in different tissues such as leaves, stalks, buds, petals, androecium, endosperm, and fibers, where the pGhFDH-driven GUS expression levels were 3-10-fold higher compared to those under the CaMV 35S promoter at 10-30 days post-anthesis (DPA) during fiber development. The results indicate that pGhFDH can be used as an endogenous constitutive promoter to drive the expression of target genes in various cotton tissues to facilitate functional genomic studies and accelerate cotton molecular breeding.

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.

Cotton pedigree genome reveals restriction of cultivar-driven strategy in cotton breeding

BACKGROUND

Many elite genes have been identified from the available cotton genomic data, providing various genetic resources for gene-driven breeding. However, backbone cultivar-driven breeding is the most widely applied strategy. Revealing the genetic basis of cultivar-driven strategy's restriction is crucial for transition of cotton breeding strategy.

RESULT

CRI12 is a backbone cultivar in cultivar-driven breeding. Here we sequence the pedigree of CRI12 using Nanopore long-read sequencing. We construct a graphical pedigree genome using the high-quality CRI12 genome and 13,138 structural variations within 20 different pedigree members. We find that low hereditary stability of elite segments in backbone cultivars is a drawback of cultivar-driven strategy. We also identify 623 functional segments in CRI12 for multiple agronomic traits in presence and absence variation-based genome-wide association study on three cohorts. We demonstrate that 25 deleterious segments are responsible for the geographical divergence of cotton in pathogen resistance. We also characterize an elite pathogen-resistant gene (GhKHCP) utilized in modern cotton breeding. In addition, we identify 386 pedigree fingerprint segments by comparing the segments of the CRI12 pedigree with those of a large cotton population.

CONCLUSION

We characterize the genetic patterns of functional segments in the pedigree of CRI12 using graphical genome method, revealing restrictions of cultivar-driven strategies in cotton breeding. These findings provide theoretical support for transitioning from cultivar-driven to gene-driven strategy in cotton breeding.

Hub Genes in Stable QTLs Orchestrate the Accumulation of Cottonseed Oil in Upland Cotton via Catalyzing Key Steps of Lipid-Related Pathways

Upland cotton is the fifth-largest oil crop in the world, with an average supply of nearly 20% of vegetable oil production. Cottonseed oil is also an ideal alternative raw material to be efficiently converted into biodiesel. However, the improvement in kernel oil content (KOC) of cottonseed has not received sufficient attention from researchers for a long time, due to the fact that the main product of cotton planting is fiber. Previous studies have tagged QTLs and identified individual candidate genes that regulate KOC of cottonseed. The regulatory mechanism of oil metabolism and accumulation of cottonseed are still elusive. In the current study, two high-density genetic maps (HDGMs), which were constructed based on a recombinant inbred line (RIL) population consisting of 231 individuals, were used to identify KOC QTLs. A total of forty-three stable QTLs were detected via these two HDGM strategies. Bioinformatic analysis of all the genes harbored in the marker intervals of the stable QTLs revealed that a total of fifty-one genes were involved in the pathways related to lipid biosynthesis. Functional analysis via coexpression network and RNA-seq revealed that the hub genes in the co-expression network that also catalyze the key steps of fatty acid synthesis, lipid metabolism and oil body formation pathways (ACX4, LACS4, KCR1, and SQD1) could jointly orchestrate oil accumulation in cottonseed. This study will strengthen our understanding of oil metabolism and accumulation in cottonseed and contribute to KOC improvement in cottonseed in the future, enhancing the security and stability of worldwide food supply.

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.

Cell cycle-dependent kinase inhibitor GhKRP6, a direct target of GhBES1.4, participates in BR regulation of cell expansion in cotton

The steroidal hormone brassinosteroid (BR) has been shown to positively regulate cell expansion in plants. However, the specific mechanism by which BR controls this process has not been fully understood. In this study, RNA-seq and DAP-seq analysis of GhBES1.4 (a core transcription factor in BR signaling) were used to identify a cotton cell cycle-dependent kinase inhibitor called GhKRP6. The study found that GhKRP6 was significantly induced by the BR hormone and that GhBES1.4 directly promoted the expression of GhKRP6 by binding to the CACGTG motif in its promoter region. GhKRP6-silenced cotton plants had smaller leaves with more cells and reduced cell size. Furthermore, endoreduplication was inhibited, which affected cell expansion and ultimately decreased fiber length and seed size in GhKRP6-silenced plants compared with the control. The KEGG enrichment results of control and VIGS-GhKRP6 plants revealed differential expression of genes related to cell wall biosynthesis, MAPK, and plant hormone transduction pathways - all of which are related to cell expansion. Additionally, some cyclin-dependent kinase (CDK) genes were upregulated in the plants with silenced GhKRP6. Our study also found that GhKRP6 could interact directly with a cell cycle-dependent kinase called GhCDKG. Taken together, these results suggest that BR signaling influences cell expansion by directly modulating the expression of cell cycle-dependent kinase inhibitor GhKRP6 via GhBES1.4.

Low soil available phosphorus level reduces cotton fiber length via osmoregulation

Introduction

Phosphorus (P) deficiency hinders cotton (Gossypium hirustum L.) growth and development, seriously affecting lint yield and fiber quality. However, it is still unclear how P fertilizer affects fiber length.

Methods

Therefore, a two-year (2019-2020) pool-culture experiment was conducted using the split-plot design, with two cotton cultivars (CCRI-79; low-P tolerant and SCRC-28; low-P sensitive) as the main plot. Three soil available phosphorus (AP) contents (P0: 3 ± 0.5, P1: 6 ± 0.5, and P2 (control) with 15 ± 0.5 mg kg-1) were applied to the plots, as the subplot, to investigate the impact of soil AP content on cotton fiber elongation and length.

Results

Low soil AP (P0 and P1) decreased the contents of the osmotically active solutes in the cotton fibers, including potassium ions (K+), malate, soluble sugar, and sucrose, by 2.2-10.2%, 14.4-47.3%, 8.7-24.5%, and 10.1-23.4%, respectively, inhibiting the vacuoles from facilitating fiber elongation through osmoregulation. Moreover, soil AP deficiency also reduced the activities of enzymes participated in fiber elongation (plasma membrane H+-ATPase (PM-H+-ATPase), vacuole membrane H+-ATPase (V-H+-ATPase), vacuole membrane H+-translocating inorganic pyrophosphatase (V-H+-PPase), and phosphoenolpyruvate carboxylase (PEPC)). The PM-H+-ATPase, V-H+-ATPase, V-H+-PPase, and PEPC were reduced by 8.4-33.0%, 7.0-33.8%, 14.1-38.4%, and 16.9-40.2%, respectively, inhibiting the transmembrane transport of the osmotically active solutes and acidified conditions for fiber cell wall, thus limiting the fiber elongation. Similarly, soil AP deficiency reduced the fiber length by 0.6-3.0 mm, mainly due to the 3.8-16.3% reduction of the maximum velocity of fiber elongation (VLmax). Additionally, the upper fruiting branch positions (FB10-11) had higher VLmax and longer fiber lengths under low soil AP.

Discussion

Cotton fibers with higher malate content and V-H+-ATPase and V-H+-PPase activities yielded longer fibers. And the malate and soluble sugar contents and V-H+-ATPase and PEPC activities in the SCRC-28's fiber were more sensitive to soil AP deficiency in contrast to those of CCRI-79, possibly explaining the SCRC-28 fiber length sensitivity to low soil AP.