Saturated very-long-chain fatty acids promote cotton fiber and Arabidopsis cell elongation by activating ethylene biosynthesis

Two duplicated GhMML3 genes coordinately control development of lint and fuzz fibers in cotton

Cotton produces two types of fibers: fuzz and lint. Cotton yield is determined by the number of epidermal cells that develop into lint fibers. Despite numerous studies, the genetic and molecular mechanisms that control lint and fuzz fiber development remain unclear. Here, using the recessive naked-seed or fuzzless-linted mutant (n2NSM) in combination with gene editing and complementation, we found that the recessive fuzzless gene n2 encodes the MYBMIXTA-like (MML) transcription factor GhMML3_D12. Overexpression of GhMML3_D12 in n2NSM restored fuzz fiber development, whereas CRISPR-Cas9 knockout of GhMML3_D12 in wild-type cotton (J668) resulted in a fuzzless-linted phenotype. Interestingly, simultaneous edits to GhMML3_D12 and its duplicate GhMML3_A12 resulted in plants with a fiberless (fuzzless-lintless) phenotype. Detailed investigation of the seed fiber phenotypes of segregating progeny derived from a cross between J668 and a fiberless gene-edited mutant of GhMML3 (#mml3s) not only identified progeny that mimicked natural fuzzless and fiberless mutants but also revealed that the duplicated GhMML3_A12 and GhMML3_D12 regulate the development of fuzz and lint fibers in a dose-dependent manner. Comparative transcriptome analysis and single-cell RNA sequencing identified GhMML3 as the central hub of the gene network that regulates fiber initiation and early-stage elongation. The gene regulatory network revealed potential candidate genes and key regulators that may contribute to fiber initiation and development, and a model for the control of lint and fuzz fiber development by GhMML3 was proposed. We also found that the GhMML3_D12 protein can bind directly to the promoters of GhHD-1 and GhMYB25, two key genes involved in fiber initiation, thereby activating their expression. This study provides new insights into the fundamental mechanisms that underlie cotton fiber development.

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

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

Cytoplasmic G6PDs modulate callus formation in Arabidopsis root explants through regulation of very-long-chain fatty acids accumulation

Glucose-6-phosphate dehydrogenase (G6PD), the rate-limiting enzyme in the pentose phosphate pathway, impacts cancer cell proliferation and plant stress responses. However, its role in plant cell dedifferentiation and callus formation is not well understood. This study explores the function of cytoplasmic G6PD isoforms in Arabidopsis pericycle cell reprogramming into callus by employing a suite of mutant analyses, qRT-PCR, and GC-MS. Our findings demonstrate that g6pd5/6 double mutants exhibit enhanced callus formation compared to wild-type and single mutants, implicating cytoplasmic G6PDs as negative regulators of callus development. The double mutant showed reduced NADPH levels and increased expression of very-long-chain fatty acid (VLCFA) biosynthesis genes and the VLCFA-downstream gene Aberrant Lateral Root Formation 4 (ALF4) on callus-inducing medium (CIM). Notably, VLCFA concentrations were decreased in g6pd5/6 mutants, and supplementation of VLCFA reduced callus area. Additionally, callus formation in the alf4/g6pd5/6 triple mutant aligned with wild-type, suggesting a redundant inhibitory function of G6PD5 and G6PD6 in the regulation of VLCFA accumulation and related signaling. Contrasting with their roles in cancer cell proliferation, our study unveils novel insights into the G6PD signaling pathway, highlighting its unique function in negatively regulating plant callus formation.

Polyploidization leads to salt stress resilience via ethylene signaling in citrus plants

Polyploidization is a common occurrence in the evolutionary history of flowering plants, significantly contributing to their adaptability and diversity. However, the molecular mechanisms behind these adaptive advantages are not well understood. Through comprehensive phenotyping of diploid and tetraploid clones from Citrus and Poncirus genera, we discovered that genome doubling significantly enhances salt stress resilience. Epigenetic and transcriptomic analyses revealed that increased ethylene production in the roots of tetraploid plants was associated with hypomethylation and enhanced chromatin accessibility of the ACO1 gene. This increased ethylene production activates the transcription of reactive oxygen species scavenging genes and stress-related hormone biosynthesis genes. Consequently, tetraploid plants exhibited superior root functionality under salt stress, maintaining improved cytosolic K+/Na+ homeostasis. To genetically validate the link between salt stress resilience and ACO1 expression, we generated overexpression and knockout lines, confirming the central role of ACO1 expression regulation following genome doubling in salt stress resilience. Our work elucidates the molecular mechanisms underlying the role of genome doubling in stress resilience. We also highlight the importance of chromatin dynamics in fine-tuning ethylene gene expression and activating salt stress resilience pathways, offering valuable insights into plant adaptation and crop genome evolution.

NAC transcription factor GbNTL9 modifies the accumulation and organization of cellulose microfibrils to enhance cotton fiber strength

INTRODUCTION

Fiber strength is a critical determinant of fiber quality, with stronger fibers being highly preferred in the cotton textile industry. However, the genetic basis and the specific regulatory mechanism underlying the formation of cotton fiber strength remain largely unknown.

OBJECTIVES

To explore fiber strength-related genes, QTL mapping, map-based cloning, and gene function verification were conducted in a backcross inbred line BS41 derived from interspecific hybridization between upland cotton and sea-island cotton.

METHODS

Upland cotton Emian22 (E22) and an interspecific backcross inbred line (BIL) BS41 were used as parents to construct secondary segregation populations for BSA and QTL mapping of fiber strength. The candidate gene GbNTL9 was identified through map-based cloning and expression analysis. The function of NTL9 was determined through transgenic experiments and cytological observations. The regulatory mechanisms of NTL9 were explored using RNA-seq, RT-qPCR, yeast two-hybrid, bimolecular fluorescence complementation, and yeast one-hybrid.

RESULTS

A major QTL for fiber strength, qFS-A11-1, was mapped to a 14.6-kb genomic region using segregating populations from E22 × BS41. GbNTL9, which encodes a NAC transcription factor, was identified as the candidate gene. Overexpression of both upland cotton genotype NTL9E22 and sea-island genotype NTL9BS41 in upland cotton enhanced fiber strength by facilitating the dense accumulation and orderly organization of cellulose microfibrils within the cell wall. Transcriptomic analysis revealed that NTL9 inhibited the expression of genes involved in secondary wall synthesis, such as CESA4, CESA7, and CESA8, thereby delaying cell wall cellulose deposition and altering the microfibril deposition pattern. NTL9 interacted with MYB6 and functioned as a downstream gene in the ethylene signaling pathway. Additionally, an effective gene marker NTL9-24 was developed to distinguish haplotypes from G. barbadense and G. hirsutum for fiber quality breeding program.

CONCLUSION

Our findings demonstrate that GbNTL9 positively regulates fiber strength through altering the microfibril deposition pattern, and provide a new insight into the molecular mechanism underlying fiber strength.

The Disruptions of Sphingolipid and Sterol Metabolism in the Short Fiber of Ligon-Lintless-1 Mutant Revealed Obesity Impeded Cotton Fiber Elongation and Secondary Cell Wall Deposition

Boosting evidence indicated lipids play important roles in plants. To explore lipid function in cotton fiber development, the lipid composition and content were detected by untargeted and targeted lipidomics. Compared with rapid elongation fibers, the lipid intensity of 16 sub-classes and 56 molecular species decreased, while only 7 sub-classes and 26 molecular species increased in the fibers at the stage of secondary cell wall deposition. Unexpectedly, at the rapid elongation stage, 20 sub-classes and 60 molecular species increased significantly, while only 5 sub-classes and 8 molecular species decreased in the ligon lintless-1 (li-1) mutant compared with its wild-type Texas Maker-1 (TM-1). Particularly, campesteryl, sitosteryl, and total steryl ester increased by 21.8-, 48.7-, and 45.5-fold in the li-1 fibers, respectively. All the molecular species of sphingosine-1-P, phytoceramide-OHFA, and glucosylceramide increased while all sphingosine, phytosphingosine, and glycosyl inositol phospho ceramides decreased in the li-1 fibers. Similarly, the different expression genes between the mutant and wild type were enriched in many pathways involved in the lipid metabolism. Furthermore, the number of lipid droplets also increased in the li-1 leaf and fiber cells when compared with the wild type. These results illuminated that fiber cell elongation being blocked in the li-1 mutant was not due to a lack of lipids, but rather lipid over-accumulation (obesity), which may result from the disruption of sphingolipid and sterol metabolism. This study provides a new perspective for further studying the regulatory mechanisms of fiber development.

The transcription factor GhMYB4 represses lipid transfer and sucrose transporter genes and inhibits fiber cell elongation in cotton

Cotton (Gossypium hirsutum) fiber is a highly elongated single cell with a thickened cell wall. MYB transcription factors are important regulators of plant cell elongation; however, the molecular mechanism involved in regulating fiber elongation remains to be explored. Here, we present evidence that the R2R3-MYB transcription factor GhMYB4 negatively regulates cotton fiber cell elongation by suppressing the expression of 2 crucial genes previously reported to affect fiber development: lipid transfer protein 4 (GhLTP4) and sucrose transporter 12 (GhSWEET12). GhMYB4 is preferentially expressed in elongating fiber cells. Knockdown of GhMYB4 in cotton results in longer fiber cells, whereas overexpression of GhMYB4 in Arabidopsis leads to reduced plant height and root length. Transcriptomic and lipidomic analyses revealed that GhMYB4 is involved in coordinating 3 interconnected biological processes, namely lipid content regulation, auxin signaling, and sugar metabolism. Additionally, we showed that GhMYB4 inhibits the expression of GhLTP4 and GhSWEET12 by binding to the MYB cis-element (TTTAGTG) in their respective promoters. Interestingly, basic helix-loop-helix transcription factor 105 (GhbHLH105) and MYB transcription factor 212 (GhMYB212) counteract the inhibitory effects of GhMYB4 on the expression of GhLTP4 and GhSWEET12, respectively. These findings provide insights into the complex molecular mechanisms regulating plant cell elongation.

Genome-Wide In Silico Analysis of 1-Aminocyclopropane-1-carboxylate oxidase (ACO) Gene Family in Rice (Oryza sativa L.)

The plant hormone ethylene elicits crucial regulatory effects on plant growth, development, and stress resistance. As the enzyme that catalyzes the final step of ethylene biosynthesis, 1-Aminocyclopropane-1-carboxylic acid oxidase (ACO) plays a key role in precisely controlling ethylene production. However, the functional characterization of the ACO gene family in rice remains largely unexplored. In this study, we performed a phylogenetic analysis of seven OsACO genes (OsACO1-OsACO7), which were classified into three subfamilies (Types I, II, and III). The members within the same clades exhibited similar tertiary structures and conserved protein motifs. We conducted inter/intraspecies covariance assays of OsACOs to elucidate their evolutionary and duplication events. Numerous cis-acting elements identified in OsACO promoter regions are associated with development, hormonal stimuli, and environmental responses. The expression assay by RT-qPCR revealed that OsACO genes exhibited tissue-specific expression and were significantly altered under various abiotic stresses, indicating their potential involvement in these processes regulated at the transcriptional level. Additionally, we predicted candidate-targeting miRNAs and identified putative cysteine sites of S-nitrosylation (SNO) and S-sulfhydration (SSH) in OsACOs, providing insights into their post-transcriptional and post-translational regulatory mechanisms. These findings pave the way for the further investigation of OsACO functions and their potential applications in improving rice growth and stress resilience by modulating ethylene biosynthesis.

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.

Brassinosteroids in cotton: orchestrating fiber development

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

Research progress on the mechanisms of fruit glossiness in cucumber

Cucumber (Cucumis sativus L.) is an important horticultural crop in China. Consumer requirements for aesthetically pleasing appearances of horticultural crops are gradually increasing, and cucumbers having a good visual appearance, as well as flavor, are important for breeding and industry development. The gloss of cucumber fruit epidermis is an important component of its appeal, and the wax layer on the fruit surface plays important roles in plant growth and forms a powerful barrier against external biotic and abiotic stresses. The wax of the cucumber epidermis is mainly composed of alkanes, and the luster of cucumber fruit is mainly determined by the alkane and silicon contents of the epidermis. Several genes, transcription factors, and transporters affect the synthesis of ultra-long-chain fatty acids and change the silicon content, further altering the gloss of the epidermis. However, the specific regulatory mechanisms are not clear. Here, progress in research on the luster of cucumber fruit epidermis from physiological, biochemical, and molecular regulatory perspectives are reviewed. Additionally, future research avenues in the field are discussed.

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.

Evolutionary comparison of lncRNAs in four cotton species and functional identification of LncR4682-PAS2-KCS19 module in fiber elongation

Long non-coding RNAs (lncRNAs) play an important role in various biological processes in plants. However, there have been few reports on the evolutionary signatures of lncRNAs in closely related cotton species. The lncRNA transcription patterns in two tetraploid cotton species and their putative diploid ancestors were compared in this paper. By performing deep RNA sequencing, we identified 280 429 lncRNAs from 21 tissues in four cotton species. lncRNA transcription evolves more rapidly than mRNAs, and exhibits more severe turnover phenomenon in diploid species compared to that in tetraploid species. Evolutionarily conserved lncRNAs exhibit higher expression levels, and lower tissue specificity compared with species-specific lncRNAs. Remarkably, tissue expression of homologous lncRNAs in Gossypium hirsutum and G. barbadense exhibited similar patterns, suggesting that these lncRNAs may be functionally conserved and selectively maintained during domestication. An orthologous lncRNA, lncR4682, was identified and validated in fibers of G. hirsutum and G. barbadense with the highest conservatism and expression abundance. Through virus-induced gene silencing in upland cotton, we found that lncR4682 and its target genes GHPAS2 and GHKCS19 positively regulated fiber elongation. In summary, the present study provides a systematic analysis of lncRNAs in four closely related cotton species, extending the understanding of transcriptional conservation of lncRNAs across cotton species. In addition, LncR4682-PAS2-KCS19 contributes to cotton fiber elongation by participating in the biosynthesis of very long-chain fatty acids.

Comparative Lipidomics Analysis Provides New Insights into the Metabolic Basis of Color Formation in Green Cotton Fiber

Green fiber (GF) is a naturally colored fiber. A limited understanding of its color formation mechanism restricts the improvement of colored cotton quality. This experiment used upland cotton green fiber germplasm 1-4560 and genetic inbred line TM-1; the lipid profiles of green fibers at 30 (white stage) and 35 days post-anthesis (DPA) (early greening stage), as well as those of TM-1 at the same stages, were revealed. Among the 109 differential types of lipids (DTLs) unique to GF, the content of phosphatidylserine PS (16:0_18:3) was significantly different at 30 and 35 DPA. It is speculated that this lipid is crucial for the pigment accumulation and color formation process of green fibers. The 197 DTLs unique to TM-1 may be involved in white fiber (WF) development. Among the shared DTLs in GF35 vs. GF30 and WF35 vs. WF30, sulfoquinovosyldiacyl-glycerol SQDG (18:1_18:1) displays a significant difference in the content change between green fibers and white fibers, potentially affecting color formation through changes in content. The enriched metabolic pathways in both comparison groups are relatively conserved. In the most significantly enriched glycerophospholipid metabolic pathway, 1-acyl-sn-glycero-3-phosphocholine (C04230) only appears in white cotton. This indicates differences in the metabolic pathways between white and green fibers, potentially related to different mechanisms of color formation and fiber development. These findings provide a new theoretical basis for studying cotton fiber development and offer important insights into the specific mechanism of green fiber color formation.

GhATL68b regulates cotton fiber cell development by ubiquitinating the enzyme required for β-oxidation of polyunsaturated fatty acids

E3 ligases are key enzymes required for protein degradation. Here, we identified a C3H2C3 RING domain-containing E3 ubiquitin ligase gene named GhATL68b. It is preferentially and highly expressed in developing cotton fiber cells and shows greater conservation in plants than in animals or archaea. The four orthologous copies of this gene in various diploid cottons and eight in the allotetraploid G. hirsutum were found to have originated from a single common ancestor that can be traced back to Chlamydomonas reinhardtii at about 992 million years ago. Structural variations in the GhATL68b promoter regions of G. hirsutum, G. herbaceum, G. arboreum, and G. raimondii are correlated with significantly different methylation patterns. Homozygous CRISPR-Cas9 knockout cotton lines exhibit significant reductions in fiber quality traits, including upper-half mean length, elongation at break, uniformity, and mature fiber weight. In vitro ubiquitination and cell-free protein degradation assays revealed that GhATL68b modulates the homeostasis of 2,4-dienoyl-CoA reductase, a rate-limiting enzyme for the β-oxidation of polyunsaturated fatty acids (PUFAs), via the ubiquitin proteasome pathway. Fiber cells harvested from these knockout mutants contain significantly lower levels of PUFAs important for production of glycerophospholipids and regulation of plasma membrane fluidity. The fiber growth defects of the mutant can be fully rescued by the addition of linolenic acid (C18:3), the most abundant type of PUFA, to the ovule culture medium. This experimentally characterized C3H2C3 type E3 ubiquitin ligase involved in regulating fiber cell elongation may provide us with a new genetic target for improved cotton lint production.

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.

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.

Transcriptome and metabolome atlas reveals contributions of sphingosine and chlorogenic acid to cold tolerance in Citrus

Citrus is one of the most important fruit crop genera in the world, but many Citrus species are vulnerable to cold stress. Ichang papeda (Citrus ichangensis), a cold-hardy citrus species, holds great potential for identifying valuable metabolites that are critical for cold tolerance in Citrus. However, the metabolic changes and underlying mechanisms that regulate Ichang papeda cold tolerance remain largely unknown. In this study, we compared the metabolomes and transcriptomes of Ichang papeda and HB pummelo (Citrus grandis "Hirado Buntan", a cold-sensitive species) to explore the critical metabolites and genes responsible for cold tolerance. Metabolomic analyses led to the identification of common and genotype-specific metabolites, consistent with transcriptomic alterations. Compared to HB pummelo under cold stress, Ichang papeda accumulated more sugars, flavonoids, and unsaturated fatty acids, which are well-characterized metabolites involved in stress responses. Interestingly, sphingosine and chlorogenic acid substantially accumulated only in Ichang papeda. Knockdown of CiSPT (C. ichangensis serine palmitoyltransferase) and CiHCT2 (C. ichangensis hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyltransferase2), two genes involved in sphingosine and chlorogenic acid biosynthesis, dramatically decreased endogenous sphingosine and chlorogenic acid levels, respectively. This reduction in sphingosine and chlorogenic acid notably compromised the cold tolerance of Ichang papeda, whereas exogenous application of these metabolites increased plant cold tolerance. Taken together, our findings indicate that greater accumulation of a spectrum of metabolites, particularly sphingosine and chlorogenic acid, promotes cold tolerance in cold-tolerant citrus species. These findings broaden our understanding of plant metabolic alterations in response to cold stress and provide valuable targets that can be manipulated to improve Citrus cold tolerance.

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.

Dual transcriptomic analysis reveals early induced Castanea defense-related genes and Phytophthora cinnamomi effectors

Phytophthora cinnamomi Rands devastates forest species worldwide, causing significant ecological and economic impacts. The European chestnut (Castanea sativa) is susceptible to this hemibiotrophic oomycete, whereas the Asian chestnuts (Castanea crenata and Castanea mollissima) are resistant and have been successfully used as resistance donors in breeding programs. The molecular mechanisms underlying the different disease outcomes among chestnut species are a key foundation for developing science-based control strategies. However, these are still poorly understood. Dual RNA sequencing was performed in C. sativa and C. crenata roots inoculated with P. cinnamomi. The studied time points represent the pathogen's hemibiotrophic lifestyle previously described at the cellular level. Phytophthora cinnamomi expressed several genes related to pathogenicity in both chestnut species, such as cell wall-degrading enzymes, host nutrient uptake transporters, and effectors. However, the expression of effectors related to the modulation of host programmed cell death (elicitins and NLPs) and sporulation-related genes was higher in the susceptible chestnut. After pathogen inoculation, 1,556 and 488 genes were differentially expressed by C. crenata and C. sativa, respectively. The most significant transcriptional changes occur at 2 h after inoculation (hai) in C. sativa and 48 hai in C. crenata. Nevertheless, C. crenata induced more defense-related genes, indicating that the resistant response to P. cinnamomi is controlled by multiple loci, including several pattern recognition receptors, genes involved in the phenylpropanoid, salicylic acid and ethylene/jasmonic acid pathways, and antifungal genes. Importantly, these results validate previously observed cellular responses for C. crenata. Collectively, this study provides a comprehensive time-resolved description of the chestnut-P. cinnamomi dynamic, revealing new insights into susceptible and resistant host responses and important pathogen strategies involved in disease development.

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.

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.

Combined genome and transcriptome analysis of elite fiber quality in Gossypium barbadense

Gossypium barbadense, which is one of several species of cotton, is well known for its superior fiber quality. However, the genetic basis of its high-quality fiber remains largely unexplored. Here, we resequenced 269 G. barbadense accessions. Phylogenetic structure analysis showed that the set of accessions was clustered into 3 groups: G1 and G2 mainly included modern cultivars from Xinjiang, China, and G3 was related to widely introduced accessions in different regions worldwide. A genome-wide association study of 5 fiber quality traits across multiple field environments identified a total of 512 qtls (main-effect QTLs) and 94 qtlEs (QTL-by-environment interactions) related to fiber quality, of which 292 qtls and 57 qtlEs colocated with previous studies. We extracted the genes located in these loci and performed expression comparison, local association analysis, and introgression segment identification. The results showed that high expression of hormone-related genes during fiber development, introgressions from Gossypium hirsutum, and the recombination of domesticated elite allelic variation were 3 major contributors to improve the fiber quality of G. barbadense. In total, 839 candidate genes with encoding region variations associated with elite fiber quality were mined. We confirmed that haplotype GB_D03G0092H traced to G. hirsutum introgression, with a 1-bp deletion leading to a frameshift mutation compared with GB_D03G0092B, significantly improved fiber quality. GB_D03G0092H is localized in the plasma membrane, while GB_D03G0092B is in both the nucleus and plasma membrane. Overexpression of GB_D03G0092H in Arabidopsis (Arabidopsis thaliana) significantly improved the elongation of longitudinal cells. Our study systematically reveals the genetic basis of the superior fiber quality of G. barbadense and provides elite segments and gene resources for breeding high-quality cotton cultivars.

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.

Comprehensive analysis of Eleutherococcus senticosus (Rupr. & Maxim.) Maxim. fruits based on UPLC-MS/MS and GC-MS: A rapid qualitative analysis

Eleutherococcus senticosus (Rupr. & Maxim.) Maxim. fruits (ESF), as a natural edible fruit, has long been popularized. However, few studies have conducted comprehensive chemical analyses of it. This study aimed to assess nonvolatile, volatile, and fatty oil components of ESF and to preliminarily explore the antioxidant activities. The qualitative and quantitative analyses of volatile and fatty oil components of ESF from 15 different regions were performed by the gas chromatography-mass spectrometry (GC-MS). Totally, 37 and 28 compounds were identified from volatile oil and fatty oil, respectively. The ultra-high-performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UPLC-QTOF-MS/MS) was used to accurately detect 43 compounds of nonvolatile components. The volatile and fatty oil components and nonvolatile components of ESF were used as samples to determine the antioxidant activity of 2,2-diphenyl-1-picrylhydrazyl (DPPH) in vitro. The components of ESF had antioxidant activity, and the nonvolatile components had stronger antioxidant activity. The results revealed that the proposed method, which is of great significance for the screening of new active ingredients, is valuable for the identification of pharmaceutical component and further development of food industry.

The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum

KEY MESSAGE

GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

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.

Genome-wide in-silico analysis of ethylene biosynthesis gene family in Musa acuminata L. and their response under nutrient stress

Ethylene is a gaseous phytohormone involved in plants' growth and developmental processes, including seed germination, root initiation, fruit ripening, flower and leaf senescence, abscission, and stress responses. Ethylene biosynthesis (EB) gene analysis in response to nitrogen (N) and potassium (K) stress has not yet been conducted in Musa acuminata (banana) roots. The genome mining of banana (Musa acuminata L.) revealed 14 putative 1-aminocyclopropane-1-carboxylate synthase (ACS), 10 1-aminocyclopropane-1-carboxylate oxidase (ACO), and 3 Ethylene overproducer 1 (ETO1) genes. ACS, ACO, and ETO1 proteins possessed amino acid residues ranging from 422-684, 636-2670, and 893-969, respectively, with molecular weight (Mw) ranging from 4.93-7.55 kD, 10.1-8.3 kD and 10.1-10.78 kD. The number of introns present in ACS, ACO, and ETO1 gene sequences ranges from 0-14, 1-6, and 0-6, respectively. The cis-regulatory element analysis revealed the presence of light-responsive, abscisic acid, seed regulation, auxin-responsive, gibberellin element, endosperm-specific, anoxic inducibility, low-temperature responsiveness, salicylic acid responsiveness, meristem-specific and stress-responsive elements. Comprehensive phylogenetic analyses ACS, ACO, and ETO1 genes of Banana with Arabidopsis thaliana revealed several orthologs and paralogs assisting in understanding the putative functions of these genes. The expression profile of Musa acuminata genes in root under normal and low levels of nitrogen and potassium shows that MaACS14 and MaACO6 expressed highly at normal nitrogen supply. MaACS1 expression was significantly upregulated at low potassium levels, whereas, MaACO6 gene expression was significantly downregulated. The functional divergence and site-specific selective pressures on specific gene sequences of banana have been investigated. The bioinformatics-based genome-wide assessment of the family of banana attempted in the present study could be a significant step for deciphering novel ACS, ACO, and ETO1 genes based on genome-wide expression profiling.

Physiological and transcriptional analyses reveal formation of memory under recurring drought stresses in seedlings of cotton (Gossypium hirsutum)

Plants are frequently subjected to a range of environmental stresses, including drought, salinity, cold, pathogens, and herbivore attacks. To survive in such conditions, plants have evolved a novel adaptive mechanism known as 'stress memory'. The formation of stress memories necessitates coordinated responses at the cellular, genetic/genomic, and epigenetic levels, involving altered physiological responses, gene activation, hyper-induction and chromatin modification. Cotton (Gossypium spp.) is an important economic crop with numerous applications and high economic value. In this study, we establish G. hirsutum drought memory following cycles of mild drought and re-watering treatments and analyzed memory gene expression patterns. Our findings reveal the physiological, biochemical, and molecular mechanisms underlying drought stress memory formation in G. hirsutum. Specifically, H3K4me3, a histone modification, plays a crucial role in regulating [+ /+ ] transcriptional memory. Moreover, we investigated the intergenerational inheritance of drought stress memory in G. hirsutum. Collectively, our data provides theoretical guidance for cotton breeding.

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.

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.

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.

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

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

Physiological and transcriptional regulation in Taxodium hybrid 'Zhongshanshan' leaves in acclimation to prolonged partial submergence

MAIN CONCLUSION

Taxodium 703 leaves activate fermentation, amino acids metabolism and ROS detoxification, and reduce TCA cycle and ABA biosynthesis in acclimation to prolonged partial submergence stress. Taxodium hybrid 'Zhongshanshan 703' (T. mucronatum × T. distichum; Taxodium 703) is a highly flooding-tolerant woody plant. To investigate the physiological and transcriptional regulatory mechanisms underlying its leaves in acclimation to long-term flooding, we exposed cuttings of Taxodium 703 to either non-flooding (control) or partial submergence for 2 months. The leaf tissues above (AL) and below (BL) flooding-water were separately harvested. Partial submergence decreased concentrations of chlorophyll (a + b) and dehydroascorbate (DHA) and lactate dehydrogenase (LDH) activity in AL, and reduced biomass, concentrations of succinic acid, fumaric acid and malic acid, and transcript levels of genes involved in tricarboxylic acid (TCA) cycle in BL. Under partial submergence, concentrations of starch, malondialdehyde and abscisic acid (ABA) decreased, and also mRNA levels of nine-cis-epoxycarotenoid dioxygenases that are involved in ABA biosynthesis in AL and BL of Taxodium 703. Partial submergence increased O2- content in AL, and improved concentrations of pyruvate and soluble sugars and activities of LDH and peroxidase in BL. In addition, partial submergence increased concentrations of ethanol, lactate, alanine, γ-aminobutyric acid, total amino acids and ascorbic acid (ASA), and ASA/DHA, activities of alcohol dehydrogenases (ADH) and ascorbate peroxidase, as well as transcript levels of ADH1A, ADH1B and genes involved in alanine biosynthesis and starch degradation in AL and BL of Taxodium 703. Overall, these results suggest that Taxodium 703 leaves activate fermentation, amino acids metabolism and reactive oxygen species detoxification, and maintain a steady supply of sugars, and reduce TCA cycle and ABA biosynthesis in acclimation to prolonged partial submergence stress.

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

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

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

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

A brassinosteroid transcriptional regulatory network participates in regulating fiber elongation in cotton

Brassinosteroids (BRs) participate in the regulation of plant growth and development through BRI1-EMS-SUPPRESSOR1 (BES1)/BRASSINAZOLE-RESISTANT1 (BZR1) family transcription factors. Cotton (Gossypium hirsutum) fibers are highly elongated single cells, and BRs play a vital role in the regulation of fiber elongation. However, the mode of action on how BR is involved in the regulation of cotton fiber elongation remains unexplored. Here, we generated GhBES1.4 over expression lines and found that overexpression of GhBES1.4 promoted fiber elongation, whereas silencing of GhBES1.4 reduced fiber length. DNA affinity purification and sequencing (DAP-seq) identified 1,531 target genes of GhBES1.4, and five recognition motifs of GhBES1.4 were identified by enrichment analysis. Combined analysis of DAP-seq and RNA-seq data of GhBES1.4-OE/RNAi provided mechanistic insights into GhBES1.4-mediated regulation of cotton fiber development. Further, with the integrated approach of GWAS, RNA-seq, and DAP-seq, we identified seven genes related to fiber elongation that were directly regulated by GhBES1.4. Of them, we showed Cytochrome P450 84A1 (GhCYP84A1) and 3-hydroxy-3-methylglutaryl-coenzyme A reductase 1 (GhHMG1) promote cotton fiber elongation. Overall, the present study established the role of GhBES1.4-mediated gene regulation and laid the foundation for further understanding the mechanism of BR participation in regulating fiber development.

Cell-specific clock-controlled gene expression program regulates rhythmic fiber cell growth in cotton

BACKGROUND

The epidermis of cotton ovule produces fibers, the most important natural cellulose source for the global textile industry. However, the molecular mechanism of fiber cell growth is still poorly understood.

RESULTS

Here, we develop an optimized protoplasting method, and integrate single-cell RNA sequencing (scRNA-seq) and single-cell ATAC sequencing (scATAC-seq) to systematically characterize the cells of the outer integument of ovules from wild type and fuzzless/lintless (fl) cotton (Gossypium hirsutum). By jointly analyzing the scRNA-seq data from wildtype and fl, we identify five cell populations including the fiber cell type and construct the development trajectory for fiber lineage cells. Interestingly, by time-course diurnal transcriptomic analysis, we demonstrate that the primary growth of fiber cells is a highly regulated circadian rhythmic process. Moreover, we identify a small peptide GhRALF1 that circadian rhythmically controls fiber growth possibly through oscillating auxin signaling and proton pump activity in the plasma membrane. Combining with scATAC-seq, we further identify two cardinal cis-regulatory elements (CREs, TCP motif, and TCP-like motif) which are bound by the trans factors GhTCP14s to modulate the circadian rhythmic metabolism of mitochondria and protein translation through regulating approximately one third of genes that are highly expressed in fiber cells.

CONCLUSIONS

We uncover a fiber-specific circadian clock-controlled gene expression program in regulating fiber growth. This study unprecedentedly reveals a new route to improve fiber traits by engineering the circadian clock of fiber cells.

A Quantitative Structure-Activity Relationship Approach to Determine Biotoxicity of Amide Herbicides for Ecotoxicological Risk Assessment

Amide herbicides have been widely applied in agriculture and found to be widespread and affect nontarget organisms in the environment. To better understand the biotoxicity mechanisms and determine the toxicity to the nontarget organisms for the hazard and risk assessment, five QSAR models were developed for the biotoxicity prediction of amide herbicides toward five aquatic and terrestrial organisms (including algae, daphnia, fish, earthworm and avian species), based on toxicity concentration and quantitative molecular descriptors. The results showed that the developed models complied with OECD principles for QSAR validation and presented excellent performances in predictive ability. In combination, the investigated QSAR relationship led to the toxicity mechanisms that eleven electrical descriptors (E_HOMO, E_LUMO, α_xx, α_yy, α_zz, μ, qN^-, Q_xx, Q_yy, qH⁺, and q^-), four thermodynamic descriptors (Cv, S^θ, H^θ, and ZPVE), and one steric descriptor (V_m) were strongly associated with the biotoxicity of amide herbicides. Electrical descriptors showed the greatest impacts on the toxicity of amide herbicides, followed by thermodynamic and steric descriptors.

Recent progression and future perspectives in cotton genomic breeding

Upland cotton is an important global cash crop for its long seed fibers and high edible oil and protein content. Progress in cotton genomics promotes the advancement of cotton genetics, evolutionary studies, functional genetics, and breeding, and has ushered cotton research and breeding into a new era. Here, we summarize high-impact genomics studies for cotton from the last 10 years. The diploid Gossypium arboreum and allotetraploid Gossypium hirsutum are the main focus of most genetic and genomic studies. We next review recent progress in cotton molecular biology and genetics, which builds on cotton genome sequencing efforts, population studies, and functional genomics, to provide insights into the mechanisms shaping abiotic and biotic stress tolerance, plant architecture, seed oil content, and fiber development. We also suggest the application of novel technologies and strategies to facilitate genome-based crop breeding. Explosive growth in the amount of novel genomic data, identified genes, gene modules, and pathways is now enabling researchers to utilize multidisciplinary genomics-enabled breeding strategies to cultivate "super cotton", synergistically improving multiple traits. These strategies must rise to meet urgent demands for a sustainable cotton industry.

De novo transcriptome sequencing and gene co-expression reveal a genomic basis for drought sensitivity and evidence of a rapid local adaptation on Atlas cedar (Cedrus atlantica)

Introduction

Understanding the adaptive capacity to current climate change of drought-sensitive tree species is mandatory, given their limited prospect of migration and adaptation as long-lived, sessile organisms. Knowledge about the molecular and eco-physiological mechanisms that control drought resilience is thus key, since water shortage appears as one of the main abiotic factors threatening forests ecosystems. However, our current background is scarce, especially in conifers, due to their huge and complex genomes.

Methods

Here we investigated the eco-physiological and transcriptomic basis of drought response of the climate change-threatened conifer Cedrus atlantica. We studied C. atlantica seedlings from two locations with contrasting drought conditions to investigate a local adaptation. Seedlings were subjected to experimental drought conditions, and were monitored at immediate (24 hours) and extended (20 days) times. In addition, post-drought recovery was investigated, depicting two contrasting responses in both locations (drought resilient and non-resilient). Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of drought resilience and investigate a rapid local adaptation of C. atlantica.

Results

De novo transcriptome assembly was performed for the first time in this species, providing differences in gene expression between the immediate and extended treatments, as well as among the post-drought recovery phenotypes. Weighted gene co-expression network analysis showed a regulation of stomatal closing and photosynthetic activity during the immediate drought, consistent with an isohydric dynamic. During the extended drought, growth and flavonoid biosynthesis inhibition mechanisms prevailed, probably to increase root-to-shoot ratio and to limit the energy-intensive biosynthesis of secondary metabolites. Drought sensitive individuals failed in metabolism and photosynthesis regulation under drought stress, and in limiting secondary metabolite production. Moreover, genomic differences (SNPs) were found between drought resilient and sensitive seedlings, and between the two studied locations, which were mostly related to transposable elements.

Discussion

This work provides novel insights into the transcriptomic basis of drought response of C. atlantica, a set of candidate genes mechanistically involved in its drought sensitivity and evidence of a rapid local adaptation. Our results may help guide conservation programs for this threatened conifer, contribute to advance drought-resilience research and shed light on trees' adaptive potential to current climate change.

Transcriptome Profiling Reveals Candidate Genes Related to Stipe Gradient Elongation of Flammulina filiformis

Stipe gradient elongation is an important and remarkable feature in the development of most mushroom fruiting bodies. However, its molecular mechanism has rarely been described. Here, the decreasing trend of stipe elongation and increasing trend of cell length in a gradient from the top to the base of the stipe were determined in a model basidiomycete mushroom: Flammulina filiformis. According to RNA-seq results, 1409 differentially expressed genes (DEGs) were identified among elongation region (ER), transition region (TR), and stable region (SR) samples, including 26 transcription factors (TFs). Based on Short Time-series Expression Miner (STEM) clustering of DEGs, clusters 1 and 3, with obvious expression trends that were consistent with or in contrast to the elongation rate, were screened. The cluster 1 DEGs were mainly involved in the GO cellular component category and KEGG genetic information processing class; however, the cluster 3 DEGs were mainly involved in metabolic processes. Furthermore, qRT-PCR confirmed that key genes of the long-chain fatty acid synthesis pathway were involved in stipe gradient elongation and regulated by NADPH oxidase-derived ROS signaling molecules. These findings provide an essential basis for understanding the molecular mechanism of stipe gradient elongation.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

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

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

The physiological and transcriptomic study of secondary growth in Neolamarckia cadamba stimulated by the ethylene precursor ACC

Though many biological roles of ethylene have been investigated intensively, the molecular mechanism of ethylene's action in woody plants remains unclear. In this study, we investigated the effects of exogenous 1-aminocyclopropane-1-carboxylic acid (ACC), the precursor of ethylene, on the growth of Neolamarckia cadamba seedlings, a fast-growing tropical tree. After 14 days of ACC treatment, the plants showed a reduced physiological morphology while stem diameter increased; however, this did not occur after the addition of 1-MCP. Meanwhile, the lignin content of N. cadamba also increased. Transcriptome analysis revealed that the expression of the ethylene biosynthesis and signaling genes ACC oxidase (ACO) and ethylene insensitive 3 (EIN3) were up-regulated mainly at the 6th hour and the 3rd day of the ACC treatment, respectively. The transcription levels of transcription factors, mainly in the basic helix-loop-helix (bHLH), ethylene response factor (ERF), WRKY and v-myb avian myeloblastosis viral oncogene homolog (MYB) families, involved in the ethylene signaling and secondary growth also increased significantly. Furthermore, in accordance to the increased lignification of the stem, the transcriptional level of key enzymes in the phenylalanine pathway were elevated after the ACC treatment. Our results revealed the physiological and molecular mechanisms underlying the secondary growth stimulated by exogenous ACC treatment on N. cadamba seedlings.

A very-long-chain fatty acid synthesis gene, SD38, influences plant height by activating ethylene biosynthesis in rice

As an important trait in crop breeding, plant height is associated with lodging resistance and yield. With the identification and cloning of several semi-dwarfing genes, increasing numbers of semi-dwarf cultivars have emerged, which has led to a 'green revolution' in rice (Oryza sativa) production. In this study, we identified a rice semi-dwarf mutant, semi-dwarf 38 (sd38), which showed significantly reduced cell length. SD38 encodes a fatty acid elongase, β-ketoacyl-CoA synthase, which is involved in the synthesis of very-long-chain fatty acids (VLCFAs). Expression analysis showed that SD38 was localized on the membrane of the endoplasmic reticulum, and was expressed in all analyzed tissues with differential abundance. The mutation of SD38 affected lipid metabolism in the sd38 mutant. A functional complementarity test in Saccharomyces cerevisiae indicated that SD38 was capable of complementing the deficiency of ELO3p activity in BY4741-elo3 knockout yeast cells by participating in the synthesis of C24:0 VLCFA. Significant changes were observed in the expression of genes involved in ethylene synthesis, which resulted in reduced content of the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) in the sd38 mutant. Exogenously supplied VLCFA (C24:0) increased the expression levels of OsACS3, OsACS4, and OsACO7 and the plant height of sd38 mutant seedlings, similar to the effect of exogenous application of ACC and ethephon. These results reveal a relationship among VLCFAs, ethylene biosynthesis, and plant height and improve our understanding of plant height development in crops.

Status and prospects of genome-wide association studies in cotton

Over the last two decades, the use of high-density SNP arrays and DNA sequencing have allowed scientists to uncover the majority of the genotypic space for various crops, including cotton. Genome-wide association study (GWAS) links the dots between a phenotype and its underlying genetics across the genomes of populations. It was first developed and applied in the field of human disease genetics. Many areas of crop research have incorporated GWAS in plants and considerable literature has been published in the recent decade. Here we will provide a comprehensive review of GWAS studies in cotton crop, which includes case studies on biotic resistance, abiotic tolerance, fiber yield and quality traits, current status, prospects, bottlenecks of GWAS and finally, thought-provoking question. This review will serve as a catalog of GWAS in cotton and suggest new frontiers of the cotton crop to be studied with this important tool.

Genome-wide identification and characterization of the KCS gene family in sorghum (Sorghum bicolor (L.) Moench)

The aboveground parts of plants are covered with cuticle, a hydrophobic layer composed of cutin polyester and cuticular wax that can protect plants from various environmental stresses. β-Ketoacyl-CoA synthase (KCS) is the key rate-limiting enzyme in plant wax synthesis. Although the properties of KCS family genes have been investigated in many plant species, the understanding of this gene family in sorghum is still limited. Here, a total of 25 SbKCS genes were identified in the sorghum genome, which were named from SbKCS1 to SbKCS25. Evolutionary analysis among different species divided the KCS family into five subfamilies and the SbKCSs were more closely related to maize, implying a closer evolutionary relationship between sorghum and maize. All SbKCS genes were located on chromosomes 1, 2, 3, 4, 5, 6, 9 and 10, respectively, while Chr 1 and Chr 10 contained more KCS genes than other chromosomes. The prediction results of subcellular localization showed that SbKCSs were mainly expressed in the plasma membrane and mitochondria. Gene structure analysis revealed that there was 0-1 intron in the sorghum KCS family and SbKCSs within the same subgroup were similar. Multiple cis-acting elements related to abiotic stress, light and hormone response were enriched in the promoters of SbKCS genes, which indicated the functional diversity among these genes. The three-dimensional structure analysis showed that a compact spherical space structure was formed by various secondary bonds to maintain the stability of SbKCS proteins, which was necessary for their biological activity. qRT-PCR results revealed that nine randomly selected SbKCS genes expressed differently under drought and salt treatments, among which SbKCS8 showed the greatest fold of expression difference at 12 h after drought and salt stresses, which suggested that the SbKCS genes played a potential role in abiotic stress responses. Taken together, these results provided an insight into investigating the functions of KCS family in sorghum and in response to abiotic stress.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.