The MYB transcription factor GhMYB25 regulates early fibre and trichome development

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.

A high-resolution model of gene expression during Gossypium hirsutum (cotton) fiber development

BACKGROUND

Cotton fiber development relies on complex and intricate biological processes to transform newly differentiated fiber initials into the mature, extravagantly elongated cellulosic cells that are the foundation of this economically important cash crop. Here we extend previous research into cotton fiber development by employing controlled conditions to minimize variability and utilizing time-series sampling and analyses to capture daily transcriptomic changes from early elongation through the early stages of secondary wall synthesis (6 to 24 days post anthesis; DPA).

RESULTS

A majority of genes are expressed in fiber, largely partitioned into two major coexpression modules that represent genes whose expression generally increases or decreases during development. Differential gene expression reveals a massive transcriptomic shift between 16 and 17 DPA, corresponding to the onset of the transition phase that leads to secondary wall synthesis. Subtle gene expression changes are captured by the daily sampling, which are discussed in the context of fiber development. Coexpression and gene regulatory networks are constructed and associated with phenotypic aspects of fiber development, including turgor and cellulose production. Key genes are considered in the broader context of plant secondary wall synthesis, noting their known and putative roles in cotton fiber development.

CONCLUSIONS

The analyses presented here highlight the importance of fine-scale temporal sampling on understanding developmental processes and offer insight into genes and regulatory networks that may be important in conferring the unique fiber phenotype.

Spatial transcriptome and single-cell RNA sequencing reveal the molecular basis of cotton fiber initiation development

Recent advances in single-cell transcriptomics have greatly expanded our knowledge of plant development and cellular responses. However, analyzing fiber cell differentiation in plants, particularly in cotton, remains a complex challenge. A spatial transcriptomic map of ovule from -1 DPA, 0 DPA, and 1 DPA in cotton was successfully constructed, which helps to explain the important role of sucrose synthesis and lipid metabolism during early fiber development. Additionally, single-cell RNA sequencing (scRNA-seq) further highlighted the cellular heterogeneity and identified clusters of fiber developmental marker genes. Integration of spatial and scRNA-seq data unveiled key genes SVB and SVBL involved in fiber initiation, suggesting functional redundancy between them. These findings provide a detailed molecular landscape of cotton fiber development, offering valuable insights for enhancing lint yield.

Diverse roles of MYB transcription factors in plants

MYB transcription factors (TFs), one of the largest TF families in plants, are involved in various plant-specific processes as the central regulators, such as in phenylpropanoid metabolism, cell cycle, formation of root hair and trichome, phytohormones responses, reproductive growth and abiotic or biotic stress responses. Here we summarized multiple roles and explained the molecular mechanisms of MYB TFs in plant development and stress adaptation. The exploration of MYB TFs contributes to a better comprehension of molecular regulation in plant development and environmental adaptability.

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.

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.

GhGRF4/GhARF2-GhGASA24 module regulates fiber cell wall thickness by modulating cellulose biosynthesis in upland cotton (Gossypium hirsutum)

Fiber elongation rate is an essential characteristic of cotton fiber in the textile industry, yet it has been largely overlooked in genetic studies. Gibberellins (GAs) and auxin (IAA) are recognized for their role in directing numerous developmental processes in plants by influencing cell differentiation and elongation. However, the degree to which GA-IAA interaction governs cellular elongation in cotton fiber cells remains to be fully understood. In this study, we identified a causal gene, Gibberellic Acid-Stimulated in Arabidopsis 24 (GhGASA24), that appears to be responsible for fiber elongation rate via regulating fiber cell wall thickness. Subsequent experiments revealed that GhGASA24 influences cell wall formation by promoting the expression of GhCesA8 and GhCesA10. Our findings suggest that Auxin Response Factor 2 (GhARF2) regulates fiber elongation rate by directly binding to the AuxRE elements in GhGASA24 promoter. In addition, we identified Growth Regulation Factor 4 (GhGRF4) as a transcription factor that interacts with GhARF2 to form a heterodimer complex, which also transcriptionally activates GhGASA24. Intriguingly, GhGRF4 regulates GhARF2 expression by directly binding to its promoter, thereby acting as a cascade regulator to enhance the transcriptional levels of GhGASA24. We propose that the GhGRF4/GhARF2-GhGASA24-GhCesAs module may contribute to fiber cell wall thickness by modulating cellulose biosynthesis, and provide a theoretical basis for improvement of fiber quality.

The 4Fs of cotton: genome editing of cotton for fiber, food, feed, and fuel to achieve zero hunger

Cotton is globally known for its high-priority cellulose-rich natural fiber. In addition to providing fiber for the textile industry, it is an important source material for edible oil, livestock feed, and fuel products. Global warming and the growing population are the major challenges to the world's agriculture and the potential risks to food security. In this context, improving output traits in cotton is necessary to achieve sustainable cotton production. During the last few years, high throughput omics techniques have aided in identifying crucial genes associated with traits of cotton fiber, seed, and plant architecture which could be targeted with more precision and efficiency through the CIRPSR/Cas-mediated genome editing technique. The various CRISPR/Cas systems such as CRISPR/Cas9, CRISPR/nCas9, and CRISPR/Cas12a have been employed to edit cotton genes associated with a wide range of traits including fiber length, flowering, leaf colour, rooting, seed oil, plant architecture, gossypol content, somatic embryogenesis, and biotic and abiotic stresses tolerance, highlighting its effectiveness in editing the cotton genome. Thus, CRISPR/Cas-mediated genome editing has emerged as a technique of choice to tailor crop phenotypes for better yield potential and environmental resilience. The review covers a comprehensive analysis of cotton phenotypic traits and their improvement with the help of the latest genome editing tools to improve fiber, food, feed, and fuel-associated genes of cotton to ensure food security.

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.

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

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

A dominant negative mutation of GhMYB25-like alters cotton fiber initiation, reducing lint and fuzz

Cotton (Gossypium hirsutum) fibers, vital natural textile materials, are single-cell trichomes that differentiate from the ovule epidermis. These fibers are categorized as lint (longer fibers useful for spinning) or fuzz (shorter, less useful fibers). Currently, developing cotton varieties with high lint yield but without fuzz remains challenging due to our limited knowledge of the molecular mechanisms underlying fiber initiation. This study presents the identification and characterization of a naturally occurring dominant negative mutation GhMYB25-like_AthapT, which results in a reduced lint and fuzzless phenotype. The GhMYB25-like_AthapT protein exerts its dominant negative effect by suppressing the activity of GhMYB25-like during lint and fuzz initiation. Intriguingly, the negative effect of GhMYB25-like_AthapT could be alleviated by high expression levels of GhMYB25-like. We also uncovered the role of GhMYB25-like in regulating the expression of key genes such as GhPDF2 (PROTODERMAL FACTOR 2), CYCD3; 1 (CYCLIN D3; 1), and PLD (Phospholipase D), establishing its significance as a pivotal transcription factor in fiber initiation. We identified other genes within this regulatory network, expanding our understanding of the determinants of fiber cell fate. These findings offer valuable insights for cotton breeding and contribute to our fundamental understanding of fiber development.

Sphingosine Promotes Fiber Early Elongation in Upland Cotton

Sphingolipids play an important role in cotton fiber development, but the regulatory mechanism is largely unclear. We found that serine palmitoyltransferase (SPT) enzyme inhibitors, myriocin and sphingosine (dihydrosphingosine (DHS) and phytosphingosine (PHS)), affected early fiber elongation in cotton, and we performed a sphingolipidomic and transcriptomic analysis of control and PHS-treated fibers. Myriocin inhibited fiber elongation, while DHS and PHS promoted it in a dose-effect manner. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS), we found that contents of 22 sphingolipids in the PHS-treated fibers for 10 days were changed, of which the contents of 4 sphingolipids increased and 18 sphingolipids decreased. The transcriptome analysis identified 432 differentially expressed genes (238 up-regulated and 194 down-regulated) in the PHS-treated fibers. Among them, the phenylpropanoid biosynthesis pathway is the most significant enrichment. The expression levels of transcription factors such as MYB, ERF, LBD, and bHLH in the fibers also changed, and most of MYB and ERF were up-regulated. Auxin-related genes IAA, GH3 and BIG GRAIN 1 were up-regulated, while ABPs were down-regulated, and the contents of 3 auxin metabolites were decreased. Our results provide important sphingolipid metabolites and regulatory pathways that influence fiber elongation.

GhWER controls fiber initiation and early elongation by regulating ethylene signaling pathway in cotton (Gossypium hirsutum)

Cotton fibers are specialized single-cell trichomes derived from epidermal cells, similar to root hairs and trichomes in Arabidopsis. While the MYB-bHLH-WD40 (MBW) complex has been shown to regulate initiation of both root hairs and trichomes in Arabidopsis, the role of their homologous gene in cotton fiber initiation remains unknown. In this study, we identified a R2R3 MYB transcription factor (TF), GhWER, which exhibited a significant increase in expression within the outer integument of ovule at -1.5 DPA (days post anthesis). Its expression peaked at -1 DPA and then gradually decreased. Knockout of GhWER using CRISPR technology inhibited the initiation and early elongation of fiber initials, resulting in the shorter mature fiber length. Additionally, GhWER interacted with two bHLH TF, GhDEL65 and GhbHLH121, suggesting a potential regulatory complex for fiber development. RNA-seq analysis of the outer integument of the ovule at -1.5 DPA revealed that the signal transduction pathways of ethylene, auxin and gibberellin were affected in the GhWER knockout lines. Further examination demonstrated that GhWER directly activated ethylene signaling genes, including ACS1 and ETR2. These findings highlighted the biological function of GhWER in regulating cotton fiber initiation and early elongation, which has practical significance for improving fiber quality and yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01477-6.

Characterisation of the R2R3 Myb subgroup 9 family of transcription factors in tomato

Tomato (Solanum lycopersicum) has many epidermal cell outgrowths including conical petal cells and multiple types of trichomes. These include the anther-specific trichome mesh which holds the anthers connate. The R2R3 Myb Subgroup 9 family of transcription factors is involved in development of epidermal cell outgrowths throughout the angiosperms. No previous study has examined all members of this transcription factor family in a single species. All 7 R2R3 Myb Subgroup 9 genes were isolated from tomato. They were ectopically expressed in tobacco to assess their ability to induce epidermal cell outgrowth. Endogenous expression patterns were examined by semi-quantitative RT-PCR at different stages of floral development relative to the development of anther trichomes. We report variation in the degree of epidermal cell outgrowth produced in transgenic tobacco by each ectopically expressed gene. Based on expression profile and ectopic activity, SlMIXTA-2 is likely involved in the production of leaf trichomes. SlMIXTA-2 is expressed most strongly in the leaves, and not expressed in the floral tissue. SlMYB17-2 is the best candidate for the regulation of the anther trichome mesh. SlMYB17-2 is expressed strongly in the floral tissue and produces a clear phenotype of epidermal cell outgrowths when ectopically expressed in tobacco. Analysis of the phenotypes of transgenic plants ectopically expressing all 7 genes has revealed the different extent to which members of the same transcription factor subfamily can induce cellular outgrowth.

Improvement of qualitative and quantitative traits in cotton under normal and stressed environments using genomics and biotechnological tools: A review

Due to the increasing demand for high-quality and high fiber-yielding cotton (Gossypium spp.), research into the development of stress-resilient cotton cultivars has acquired greater significance. Various biotic and abiotic stressors greatly affect cotton production and productivity, posing challenges to the future of the textile industry. Moreover, the content and quality of cottonseed oil can also potentially be influenced by future environmental conditions. Apart from conventional methods, genetic engineering has emerged as a potential tool to improve cotton fiber quality and productivity. Identification and modification of genome sequences and the expression levels of yield-related genes using genetic engineering approaches have enabled to increase both the quality and yields of cotton fiber and cottonseed oil. Herein, we evaluate the significance and molecular mechanisms associated with the regulation of cotton agronomic traits under both normal and stressful environmental conditions. In addition, the importance of gossypol, a toxic phenolic compound in cottonseed that can limit consumption by animals and humans, is reviewed and discussed.

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.

Genetic basis of nectar guide trichome variation between bumblebee- and self-pollinated monkeyflowers (Mimulus): role of the MIXTA-like gene GUIDELESS

Nectar guide trichomes play crucial ecological roles in bee-pollinated flowers, as they serve as footholds and guides for foraging bees to access the floral rewards. However, the genetic basis of natural variation in nectar guide trichomes among species remains poorly understood. In this study, we performed genetic analysis of nectar guide trichome variation between two closely related monkeyflower (Mimulus) species, the bumblebee-pollinated Mimulus lewisii and self-pollinated M. parishii. We demonstrate that a MIXTA-like R2R3-MYB gene, GUIDELESS, is a major contributor to the nectar guide trichome length variation between the two species. The short-haired M. parishii carries a recessive allele due to non-synonymous substitutions in a highly conserved motif among MIXTA-like MYB proteins. Furthermore, our results suggest that besides GUIDELESS, additional loci encoding repressors of trichome elongation also contribute to the transition from bumblebee-pollination to selfing. Taken together, these results suggest that during a pollination syndrome switch, changes in seemingly complex traits such as nectar guide trichomes could have a relatively simple genetic basis, involving just a few genes of large effects.

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

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

A eudicot MIXTA family ancestor likely functioned in both conical cells and trichomes

The MIXTA family of MYB transcription factors modulate the development of diverse epidermal features in land plants. This study investigates the evolutionary history and function of the MIXTA gene family in the early-diverging eudicot model lineage Thalictrum (Ranunculaceae), with R2R3 SBG9-A MYB transcription factors representative of the pre-core eudicot duplication and thus hereby referred to as "paleoMIXTA" (PMX). Cloning and phylogenetic analysis of Thalictrum paleoMIXTA (ThPMX) orthologs across 23 species reveal a genus-wide duplication coincident with a whole-genome duplication. Expression analysis by qPCR confirmed that the highest expression is found in carpels, while newly revealing high expression in leaves and nuanced differences between paralogs in representative polyploid species. The single-copy ortholog from the diploid species T. thalictroides (TthPMX, previously TtMYBML2), which has petaloid sepals with conical-papillate cells and trichomes on leaves, was functionally characterized by virus-induced gene silencing (VIGS), and its role in leaves was also assessed from heterologous overexpression in tobacco. Another ortholog from a species with conical-papillate cells on stamen filaments, TclPMX, was also targeted for silencing. Overexpression assays in tobacco provide further evidence that the paleoMIXTA lineage has the potential for leaf trichome function in a core eudicot. Transcriptome analysis by RNA-Seq on leaves of VIGS-treated plants suggests that TthPMX modulates leaf trichome development and morphogenesis through microtubule-associated mechanisms and that this may be a conserved pathway for eudicots. These experiments provide evidence for a combined role for paleoMIXTA orthologs in (leaf) trichomes and (floral) conical-papillate cells that, together with data from other systems, makes the functional reconstruction of a eudicot ancestor most likely as also having a combined function.

Beyond skin-deep: targeting the plant surface for crop improvement

The above-ground plant surface is a well-adapted tissue layer that acts as an interface between the plant and its surrounding environment. As such, its primary role is to protect against desiccation and maintain the gaseous exchange required for photosynthesis. Further, this surface layer provides a barrier against pathogens and herbivory, while attracting pollinators and agents of seed dispersal. In the context of agriculture, the plant surface is strongly linked to post-harvest crop quality and yield. The epidermal layer contains several unique cell types adapted for these functions, while the non-lignified above-ground plant organs are covered by a hydrophobic cuticular membrane. This review aims to provide an overview of the latest understanding of the molecular mechanisms underlying crop cuticle and epidermal cell formation, with focus placed on genetic elements contributing towards quality, yield, drought tolerance, herbivory defence, pathogen resistance, pollinator attraction, and sterility, while highlighting the inter-relatedness of plant surface development and traits. Potential crop improvement strategies utilizing this knowledge are outlined in the context of the recent development of new breeding techniques.

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

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

Activation of Gossypium hirsutum ACS6 Facilitates Fiber Development by Improving Sucrose Metabolism and Transport

Cotton fiber yield depends on the density of fiber cell initials that form on the ovule epidermis. Fiber initiation is triggered by MYB-MIXTA-like transcription factors (GhMMLs) and requires a sucrose supply. Ethylene or its precursor ACC (1-aminocyclopropane-1-carboxylic acid) is suggested to affect fiber yield. The Gossypium hirsutum (L.) genome contains 35 ACS genes (GhACS) encoding ACC synthases. Here, we explored the role of a GhACS family member in the regulation of fiber initiation. Expression analyses showed that the GhACS6.3 gene pair was specifically expressed in the ovules during fiber initiation (3 days before anthesis to 5 days post anthesis, -3 to 5 DPA), especially at -3 DPA, whereas other GhACS genes were expressed at very low or undetectable levels. The expression profile of GhACS6.3 during fiber initial development was confirmed by qRT-PCR analysis. Transgenic lines overexpressing GhACS6.3 (GhACS6.3-OE) showed increased ACC accumulation in ovules, which promoted the formation of fiber initials and fiber yield components. This was accompanied by increased transcript levels of GhMML3 and increased transcript levels of genes encoding sucrose transporters and sucrose synthase. These findings imply that GhACS6.3 activation is required for fiber initial development. Our results lay the foundation for further research on increasing cotton fiber production.

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.

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

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

Morphology, sucrose metabolism and gene network reveal the molecular mechanism of seed fiber development in poplar

Poplar is an important tree species for ecological protection, wood production, bioenergy and urban greening; it has been widely planted worldwide. However, the catkin fibers produced by female poplars can cause environmental pollution and safety hazards during spring. This study focused on Populus tomentosa, and revealed the sucrose metabolism regulatory mechanism of catkin fibers development from morphological, physiological and molecular aspects. Paraffin section suggested that poplar catkin fibers were not seed hairs and produced from the epidermal cells of funicle and placenta. Sucrose degradation via invertase and sucrose synthase played the dominant role during poplar catkin fibers development. The expression patterns revealed that sucrose metabolism-related genes played important roles during catkin fibers development. Y1H analysis indicated that there was a potential interaction between sucrose synthase 2 (PtoSUS2)/vacuolar invertase 3 (PtoVIN3) and trichome-regulating MYB transcription factors in poplar. Finally, the two key genes, PtoSUS2 and PtoVIN3, had roles in Arabidopsis trichome density, indicating that sucrose metabolism is important in poplar catkin fibers development. This study is not only helpful for clarifying the mechanism of sucrose regulation during trichome development in perennial woody plants, but also establishes a foundation to solve poplar catkin fibers pollution through genetic engineering methods.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-Wide Association Study of Lint Percentage in Gossypium hirsutum L. Races

Lint percentage is one of the most essential yield components and an important economic index for cotton planting. Improving lint percentage is an effective way to achieve high-yield in cotton breeding worldwide, especially upland cotton (Gossypium hirsutum L.). However, the genetic basis controlling lint percentage has not yet been systematically understood. Here, we performed a genome-wide association mapping for lint percentage using a natural population consisting of 189 G. hirsutum accessions (188 accessions of G. hirsutum races and one cultivar TM-1). The results showed that 274 single-nucleotide polymorphisms (SNPs) significantly associated with lint percentage were detected, and they were distributed on 24 chromosomes. Forty-five SNPs were detected at least by two models or at least in two environments, and their 5 Mb up- and downstream regions included 584 makers related to lint percentage identified in previous studies. In total, 11 out of 45 SNPs were detected at least in two environments, and their 550 Kb up- and downstream region contained 335 genes. Through RNA sequencing, gene annotation, qRT-PCR, protein-protein interaction analysis, the cis-elements of the promotor region, and related miRNA prediction, Gh_D12G0934 and Gh_A08G0526 were selected as key candidate genes for fiber initiation and elongation, respectively. These excavated SNPs and candidate genes could supplement marker and gene information for deciphering the genetic basis of lint percentage and facilitate high-yield breeding programs of G. hirsutum ultimately.

Development and Biomechanics of Grewia lasiocarpa E. Mey. Ex Harv. Trichomes Exudate

Grewia lasiocarpa E. Mey. Ex Harv., Malvaceae (forest raisin) is a tropical small tree or shrub valued for its ecological importance as well as its nutritional, antioxidant, antibacterial, and anti-cancer properties as well as its ecological and ornamental importance. Glandular and non-glandular trichomes are present on the fruits, stem bark and leaves of G. lasiocarpa and these trichomes are the first line of defense. They are important structures that plants use to combat biotic and abiotic stress. The development of G. lasiocarpa trichomes and the biomechanics of the exudates present in the glandular (capitate) trichome were investigated for the first time using advanced microscopy techniques [Scanning electron microscope (SEM) and Transmission electron microscope (TEM)]. The pressurized cuticular striations may play a role in the exudates' biomechanics, i.e., releasing secondary metabolites present in the capitate trichome, which was observed to be multidirectional. The presence of many glandular trichomes on a plant implies an increase in the amount of phytometabolites. A common precursor for the development of trichomes (non-glandular and glandular) was observed to be DNA synthesis associated with a periclinal cell division, thus the final fate of the cell is determined by cell cycle regulation, polarity, and expansion. The glandular trichomes of G. lasiocarpa are multicellular and polyglandular, while the non-glandular (glandless) trichomes are either single-celled or multicellular. Since, trichomes 'house' phytocompounds of medicinal, nutritional, and agronomical benefits; the molecular and genetic study of the glandular trichomes of Grewia lasiocarpa will be beneficial to humanity.

Regulation of capsule spine formation in castor

Castor (Ricinus communis L.) is a dicotyledonous oilseed crop that can have either spineless or spiny capsules. Spines are protuberant structures that differ from thorns or prickles. The developmental regulatory mechanisms governing spine formation in castor or other plants have remained largely unknown. Herein, using map-based cloning in 2 independent F2 populations, F2-LYY5/DL01 and F2-LYY9/DL01, we identified the RcMYB106 (myb domain protein 106) transcription factor as a key regulator of capsule spine development in castor. Haplotype analyses demonstrated that either a 4,353-bp deletion in the promoter or a single nucleotide polymorphism leading to a premature stop codon in the RcMYB106 gene could cause the spineless capsule phenotype in castor. Results of our experiments indicated that RcMYB106 might target the downstream gene RcWIN1 (WAX INDUCER1), which encodes an ethylene response factor known to be involved in trichome formation in Arabidopsis (Arabidopsis thaliana) to control capsule spine development in castor. This hypothesis, however, remains to be further tested. Nevertheless, our study reveals a potential molecular regulatory mechanism underlying the spine capsule trait in a nonmodel plant species.

Identification of MYB gene family and functional analysis of GhMYB4 in cotton (Gossypium spp.)

Myeloblastosis (MYB) transcription factors (TFs) form a large gene family involved in a variety of biological processes in plants. Little is known about their roles in the development of cotton pigment glands. In this study, 646 MYB members were identified in Gossypium hirsutum genome and phylogenetic classification was analyzed. Evolution analysis revealed assymetric evolution of GhMYBs during polyploidization and sequence divergence of MYBs in G. hirustum was preferentially happend in D sub-genome. WGCNA (weighted gene co-expression network analysis) showed that four modules had potential relationship with gland development or gossypol biosynthesis in cotton. Eight differentially expressed GhMYB genes were identified by screening transcriptome data of three pairs of glanded and glandless cotton lines. Of these, four were selected as candidate genes for cotton pigment gland formation or gossypol biosynthesis by qRT-PCR assay. Silencing of GH_A11G1361 (GhMYB4) downregulated expression of multiple genes in gossypol biosynthesis pathway, indicating it could be involved in gossypol biosynthesis. The potential protein interaction network suggests that several MYBs may have indirect interaction with GhMYC2-like, a key regulator of pigment gland formation. Our study was the systematic analysis of MYB genes in cotton pigment gland development, providing candidate genes for further study on the roles of cotton MYB genes in pigment gland formation, gossypol biosynthesis and future crop plant improvement.

Overexpression of cotton Trihelix transcription factor GhGT-3b_A04 enhances resistance to Verticillium dahliae and affects plant growth in Arabidopsis thaliana

Verticillium wilt is a soil-borne fungal disease that severely affects cotton fiber yield and quality. Herein, a cotton Trihelix family gene, GhGT-3b_A04, was strongly induced by the fungal pathogen Verticillium dahliae. Overexpression of the gene in Arabidopsis thaliana enhanced the plant's resistance to Verticillium wilt but inhibited the growth of rosette leaves. In addition, the primary root length, root hair number, and root hair length increased in GhGT-3b_A04-overexpressing plants. The density and length of trichomes on the rosette leaves also increased. GhGT-3b_A04 localized to the nucleus, and transcriptome analysis revealed that it induced gene expression for salicylic acid synthesis and signal transduction and activated gene expression for disease resistance. The gene expression for auxin signal transduction and trichome development was reduced in GhGT-3b_A04-overexpressing plants. Our results highlight important regulatory genes for Verticillium wilt resistance and cotton fiber quality improvement. The identification of GhGT-3b_A04 and other important regulatory genes can provide crucial reference information for future research on transgenic cotton breeding.

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

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

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Comparative analysis of the MYB gene family in seven Ipomoea species

The MYB transcription factors regulate plant growth, development, and defense responses. However, information about the MYB gene family in Ipomoea species is rare. Herein, we performed a comprehensive genome-wide comparative analysis of this gene family among seven Ipomoea species, sweet potato (I. batatas), I. trifida, I. triloba, I. nil, I. purpurea, I. cairica, and I. aquatic, and identified 296, 430, 411, 291, 226, 281, and 277 MYB genes, respectively. The identified MYB genes were classified into five types: 1R-MYB (MYB-related), 2R-MYB (R2R3-MYB), 3R-MYB (R1R2R3-MYB), 4R-MYB, and 5R-MYB, and the MYB-related or R2R3-MYB type was the most abundant MYB genes in the seven species. The Ipomoea MYB genes were classed into distinct subgroups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis. Analysis of gene structure and protein motifs revealed that members within the same phylogenetic group presented similar exon/intron and motif organization. The identified MYB genes were unevenly mapped on the chromosomes of each Ipomoea species. Duplication analysis indicated that segmental and tandem duplications contribute to expanding the Ipomoea MYB genes. Non-synonymous substitution (Ka) to synonymous substitution (Ks) [Ka/Ks] analysis showed that the duplicated Ipomoea MYB genes are mainly under purifying selection. Numerous cis-regulatory elements related to stress responses were detected in the MYB promoters. Six sweet potato transcriptome datasets referring to abiotic and biotic stresses were analyzed, and MYB different expression genes' (DEGs') responses to stress treatments were detected. Moreover, 10 sweet potato MYB DEGs were selected for qRT-PCR analysis. The results revealed that four responded to biotic stress (stem nematodes and Ceratocystis fimbriata pathogen infection) and six responded to the biotic stress (cold, drought, and salt). The results may provide new insights into the evolution of MYB genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

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.

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

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

Genome scale analysis of 1-aminocyclopropane-1-carboxylate oxidase gene family in G. barbadense and its functions in cotton fiber development

A class of proteins, 1-aminocyclopropane-1-carboxylate oxidase (ACO), is required in the final step of production of ethylene from its immediate precursor 1-aminocyclopropane-1-carboxylic acid (ACC). Despite the crucial and regulatory role of ACO gene family in the fiber development, it has not been thoroughly analyzed and annotated in G. barbadense genome. In the present study, we have identified and characterized all isoforms of ACO gene family from genomes of Gossypium arboreum, G. barbadense, G. hirsutum and G. raimondii. Phylogenetic analysis classified all ACO proteins into six distinct groups on the basis of maximum likelihood. Gene locus analysis and circos plots indicated the distribution and relationship of these genes in cotton genomes. Transcriptional profiling of ACO isoforms in G. arboreum, G. barbadense and G. hirsutum fiber development exhibited the highest expression in G. barbadense during early fiber elongation. Moreover, the accumulation of ACC was found highest in developing fibers of G. barbadense in comparison with other cotton species. ACO expression and ACC accumulation correlated with the fiber length in cotton species. Addition of ACC to the ovule cultures of G. barbadense significantly increased fiber elongation while ethylene inhibitors hindered fiber elongation. These findings will be helpful in dissecting the role of ACOs in cotton fiber development and pave a way towards genetic manipulations for fiber quality improvement.

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

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

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.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes for Fuzz Development in Gossypium hirsutum

Fuzzless Gossypium hirsutum mutants are ideal materials for investigating cotton fiber initiation and development. In this study, we used the fuzzless G. hirsutum mutant Xinluzao 50 FLM as the research material and combined it with other fuzzless materials for verification by RNA sequencing to explore the gene expression patterns and differences between genes in upland cotton during the fuzz period. A gene ontology (GO) enrichment analysis showed that differentially expressed genes (DEGs) were mainly enriched in the metabolic process, microtubule binding, and other pathways. A weighted gene co-expression network analysis (WGCNA) showed that two modules of Xinluzao 50 and Xinluzao 50 FLM and four modules of CSS386 and Sicala V-2 were highly correlated with fuzz. We selected the hub gene with the highest KME value among the six modules and constructed an interaction network. In addition, we selected some genes with high KME values from the six modules that were highly associated with fuzz in the four materials and found 19 common differential genes produced by the four materials. These 19 genes are likely involved in the formation of fuzz in upland cotton. Several hub genes belong to the arabinogalactan protein and GDSL lipase, which play important roles in fiber development. According to the differences in expression level, 4 genes were selected from the 19 genes and tested for their expression level in some fuzzless materials. The modules, hub genes, and common genes identified in this study can provide new insights into the formation of fiber and fuzz, and provide a reference for molecular design breeding for the genetic improvement of cotton fiber.

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

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

Downregulation of Brassica napus MYB69 (BnMYB69) increases biomass growth and disease susceptibility via remodeling phytohormone, chlorophyll, shikimate and lignin levels

MYB transcription factors are major actors regulating plant development and adaptability. Brassica napus is a staple oil crop and is hampered by lodging and diseases. Here, four B. napus MYB69 (BnMYB69s) genes were cloned and functionally characterized. They were dominantly expressed in stems during lignification. BnMYB69 RNA interference (BnMYB69i) plants showed considerable changes in morphology, anatomy, metabolism and gene expression. Stem diameter, leaves, roots and total biomass were distinctly larger, but plant height was significantly reduced. Contents of lignin, cellulose and protopectin in stems were significantly reduced, accompanied with decrease in bending resistance and Sclerotinia sclerotiorum resistance. Anatomical detection observed perturbation in vascular and fiber differentiation in stems, but promotion in parenchyma growth, accompanied with changes in cell size and cell number. In shoots, contents of IAA, shikimates and proanthocyanidin were reduced, while contents of ABA, BL and leaf chlorophyll were increased. qRT-PCR revealed changes in multiple pathways of primary and secondary metabolisms. IAA treatment could recover many phenotypes and metabolisms of BnMYB69i plants. However, roots showed trends opposite to shoots in most cases, and BnMYB69i phenotypes were light-sensitive. Conclusively, BnMYB69s might be light-regulated positive regulators of shikimates-related metabolisms, and exert profound influences on various internal and external plant traits.

Fine-mapping and candidate gene analysis of qFL-c10-1 controlling fiber length in upland cotton (Gossypium hirsutum L.)

A fiber length QTL, qFL-c10-1, was fine-mapped to a 96.5-kb region containing one gene that has not been characterized in plants. Fiber length is an important component of cotton fiber quality, which is associated with other quality properties such as fiber strength, fiber maturity, and fineness. In our previous studies, a stable QTL qFL-c10-1 controlling fiber length had been identified on chromosome A10 in an upland cotton recombinant inbred line (RIL) population from a cross between Jimian5 and DH962. To fine-map qFL-c10-1, an F₂ population with 1081 individual plants from a cross between a recombinant line DJ61 and Jimian5 was established. Using linkage analysis and progeny recombination experiment, qFL-c10-1 was mapped into a 96.5-kb genomic region that just contained one proper transcript Ghir_A10G022020 (described as GhFL10), an undescribed gene in plants. One 214-bp deletion was identified in the promoter region of DJ61 compared with Jimian5. Quantitative real-time PCR (qRT-PCR) and comparative analysis of parental sequences suggested that GhFL10 was the most promising candidate gene for qFL-c10-1. According to RNA-seq, yeast two-hybrid assay and bimolecular fluorescence complementation (BiFC), we speculate that GhFL10 interacts with NF-YA transcription factors to negatively regulate fiber elongation.

Single-cell RNA-seq reveals fate determination control of an individual fibre cell initiation in cotton (Gossypium hirsutum)

Cotton fibre is a unicellular seed trichome, and lint fibre initials per seed as a factor determines fibre yield. However, the mechanisms controlling fibre initiation from ovule epidermis are not understood well enough. Here, with single-cell RNA sequencing (scRNA-seq), a total of 14 535 cells were identified from cotton ovule outer integument of Xu142_LF line at four developmental stages (1.5, 1, 0.5 days before anthesis and the day of anthesis). Three major cell types, fibre, non-fibre epidermis and outer pigment layer were identified and then verified by RNA in situ hybridization. A comparative analysis on scRNA-seq data between Xu142 and its fibreless mutant Xu142 fl further confirmed fibre cluster definition. The developmental trajectory of fibre cell was reconstructed, and fibre cell was identified differentiated at 1 day before anthesis. Gene regulatory networks at four stages revealed the spatiotemporal pattern of core transcription factors, and MYB25-like and HOX3 were demonstrated played key roles as commanders in fibre differentiation and tip-biased diffuse growth respectively. A model for early development of a single fibre cell was proposed here, which sheds light on further deciphering mechanism of plant trichome and the improvement of cotton fibre yield.

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

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

The R2R3-MYB gene CgMYB4 is involved in the regulation of cell differentiation and fiber development in the stamens of Chelone glabra L

A Plantaginaceae flowering plant, Chelone glabra, is different from Arabidopsis thaliana and cotton (Gossypium hirsutum), as it produces fibers on the anther surface. However, the evolutionary molecular mechanism of how fiber development is controlled in the stamen is unclear. MYB genes are essential transcription factors for trichome and fiber development in plants. In this study, we isolated 29 MYB domain-containing sequences using early-stage anthers and several sets of degenerated primers conserved in the R2R3 domain of the MYB transcription factor. Among them, CgMYB4 is an R2R3-MYB gene encoding 281 amino acids. Phylogenetic analysis showed that CgMYB4 is closely related to GhMYB25L/AmMIXTA, which controls fiber initiation and development in cotton and epidermal cell differentiation in the petals of Antirrhinum. Semiquantitative RT-PCR analysis showed that CgMYB4 is strongly expressed at the stamens and carpels. Overexpression of CgMYB4 significantly enhanced root hair formation in transformed hairy roots, contrary to the root hair numbers, which were reduced in silenced CgMYB4 hairy roots. Moreover, overexpression of CgMYB4 also evidently promoted fiber development at filaments and conical cell-like epidermal cell increases at the anther wall. Our results showed that CgMYB4 is an R2R3-MYB gene and is positively involved in regulating cell division and fiber differentiation in the early stages of stamen development in C. glabra.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Cymbidium goeringii

The MYB gene family plays a vital regulatory role in plant metabolism, stress response, and floral color. The R2R3-MYB gene family of C. goeringii was identified, and its expression was analyzed using bioinformatics in this article. The R2R3-MYB genes of Arabidopsis thaliana were used as a reference to determine 104 CgMYB genes and categorize them into 22 subfamilies. Exon/intron organizations and conserved motif analysis revealed that the majority of CgMYB genes were highly conserved, and chromosome localization and collinearity analysis provided evidence of tandem duplication and segmental duplication events, indicating the phenomenon of gene family expansion and contraction. The function of CgMYB genes was analyzed by cis-acting element and gene ontology (GO) enrichment. In addition, we selected CgMYB91 and CgMYB32 for RT-qPCR, suggesting that CgMYB91 and CgMYB32 are associated with anthocyanin formation. In short, this study provides a comprehensive and specific function of the R2R3-MYB transcription factors (TFs) in orchids.