The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis.

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

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

BACKGROUND: Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors. RESULTS: In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms. CONCLUSIONS: These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

A predictive model for ethylene-mediated auxin and cytokinin patterning in the Arabidopsis root.

The interaction between auxin and cytokinin is important in many aspects of plant development. Experimental measurements of both auxin and cytokinin concentration and reporter gene expression clearly show the coexistence of auxin and cytokinin concentration patterning in Arabidopsis root development. However, in the context of crosstalk among auxin, cytokinin, and ethylene, little is known about how auxin and cytokinin concentration patterns simultaneously emerge and how they regulate each other in the Arabidopsis root. This work utilizes a wide range of experimental observations to propose a mechanism for simultaneous patterning of auxin and cytokinin concentrations. In addition to revealing the regulatory relationships between auxin and cytokinin, this mechanism shows that ethylene signaling is an important factor in achieving simultaneous auxin and cytokinin patterning, while also predicting other experimental observations. Combining the mechanism with a realistic in silico root model reproduces experimental observations of both auxin and cytokinin patterning. Predictions made by the mechanism can be compared with a variety of experimental observations, including those obtained by our group and other independent experiments reported by other groups. Examples of these predictions include patterning of auxin biosynthesis rate, changes in PIN1 and PIN2 patterns in pin3,4,7 mutants, changes in cytokinin patterning in the pls mutant, PLS patterning, and various trends in different mutants. This research reveals a plausible mechanism for simultaneous patterning of auxin and cytokinin concentrations in Arabidopsis root development and suggests a key role for ethylene pattern integration.

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

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

Single-Cell Transcriptome Atlas and Regulatory Dynamics in Developing Cotton Anthers.

Plant anthers are composed of different specialized cell types with distinct roles in plant reproduction. High temperature (HT) stress causes male sterility, resulting in crop yield reduction. However, the spatial expression atlas and regulatory dynamics during anther development and in response to HT remain largely unknown. Here, the first single-cell transcriptome atlas and chromatin accessibility survey in cotton anther are established, depicting the specific expression and epigenetic landscape of each type of cell in anthers. The reconstruction of meiotic cells, tapetal cells, and middle layer cell developmental trajectories not only identifies novel expressed genes, but also elucidates the precise degradation period of middle layer and reveals a rapid function transition of tapetal cells during the tetrad stage. By applying HT, heterogeneity in HT response is shown among cells of anthers, with tapetal cells responsible for pollen wall synthesis are most sensitive to HT. Specifically, HT shuts down the chromatin accessibility of genes specifically expressed in the tapetal cells responsible for pollen wall synthesis, such as QUARTET 3 (QRT3) and CYTOCHROME P450 703A2 (CYP703A2), resulting in a silent expression of these genes, ultimately leading to abnormal pollen wall and male sterility. Collectively, this study provides substantial information on anthers and provides clues for heat-tolerant crop creation.

The elicitor VP2 from Verticillium dahliae triggers defence response in cotton.

Verticillium dahliae is a widespread and destructive soilborne vascular pathogenic fungus that causes serious diseases in dicot plants. Here, comparative transcriptome analysis showed that the number of genes upregulated in defoliating pathotype V991 was significantly higher than in the non-defoliating pathotype 1cd3-2 during the early response of cotton. Combined with analysis of the secretome during the V991-cotton interaction, an elicitor VP2 was identified, which was highly upregulated at the early stage of V991 invasion, but was barely expressed during the 1cd3-2-cotton interaction. Full-length VP2 could induce cell death in several plant species, and which was dependent on NbBAK1 but not on NbSOBIR1 in N. benthamiana. Knock-out of VP2 attenuated the pathogenicity of V991. Furthermore, overexpression of VP2 in cotton enhanced resistance to V. dahliae without causing abnormal plant growth and development. Several genes involved in JA, SA and lignin synthesis were significantly upregulated in VP2-overexpressing cotton. The contents of JA, SA, and lignin were also significantly higher than in the wild-type control. In summary, the identified elicitor VP2, recognized by the receptor in the plant membrane, triggers the cotton immune response and enhances disease resistance.

Construction of Host Plant Insect-Resistance Mutant Library by High-Throughput CRISPR/Cas9 System and Identification of A Broad-Spectrum Insect Resistance Gene.

Insects pose significant challenges in cotton-producing regions. Here, they describe a high-throughput CRISPR/Cas9-mediated large-scale mutagenesis library targeting endogenous insect-resistance-related genes in cotton. This library targeted 502 previously identified genes using 968 sgRNAs, generated ≈2000 T0 plants and achieved 97.29% genome editing with efficient heredity, reaching upto 84.78%. Several potential resistance-related mutants (10% of 200 lines) their identified that may contribute to cotton-insect molecular interaction. Among these, they selected 139 and 144 lines showing decreased resistance to pest infestation and targeting major latex-like protein 423 (GhMLP423) for in-depth study. Overexpression of GhMLP423 enhanced insect resistance by activating the plant systemic acquired resistance (SAR) of salicylic acid (SA) and pathogenesis-related (PR) genes. This activation is induced by an elevation of cytosolic calcium [Ca2+ ]cyt flux eliciting reactive oxygen species (ROS), which their demoted in GhMLP423 knockout (CR) plants. Protein-protein interaction assays revealed that GhMLP423 interacted with a human epidermal growth factor receptor substrate15 (EPS15) protein at the cell membrane. Together, they regulated the systemically propagating waves of Ca2+ and ROS, which in turn induced SAR. Collectively, this large-scale mutagenesis library provides an efficient strategy for functional genomics research of polyploid plant species and serves as a solid platform for genetic engineering of insect resistance.

Polyethyleneimine-coated MXene quantum dots improve cotton tolerance to Verticillium dahliae by maintaining ROS homeostasis.

Verticillium dahliae is a soil-borne hemibiotrophic fungal pathogen that threatens cotton production worldwide. In this study, we assemble the genomes of two V. dahliae isolates: the more virulence and defoliating isolate V991 and nondefoliating isolate 1cd3-2. Transcriptome and comparative genomics analyses show that genes associated with pathogen virulence are mostly induced at the late stage of infection (Stage II), accompanied by a burst of reactive oxygen species (ROS), with upregulation of more genes involved in defense response in cotton. We identify the V991-specific virulence gene SP3 that is highly expressed during the infection Stage II. V. dahliae SP3 knock-out strain shows attenuated virulence and triggers less ROS production in cotton plants. To control the disease, we employ polyethyleneimine-coated MXene quantum dots (PEI-MQDs) that possess the ability to remove ROS. Cotton seedlings treated with PEI-MQDs are capable of maintaining ROS homeostasis with enhanced peroxidase, catalase, and glutathione peroxidase activities and exhibit improved tolerance to V. dahliae. These results suggest that V. dahliae trigger ROS production to promote infection and scavenging ROS is an effective way to manage this disease. This study reveals a virulence mechanism of V. dahliae and provides a means for V. dahliae resistance that benefits cotton production.

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

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

The chromosome-scale reference genome of mirid bugs (Adelphocoris suturalis) genome provides insights into omnivory, insecticide resistance, and survival adaptation.

BACKGROUND: Adelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides. RESULTS: Here, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering. CONCLUSIONS: The high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies.

Single-cell resolution analysis reveals the preparation for reprogramming the fate of stem cell niche in cotton lateral meristem.

BACKGROUND: Somatic embryogenesis is a major process for plant regeneration. However, cell communication and the gene regulatory network responsible for cell reprogramming during somatic embryogenesis are still largely unclear. Recent advances in single-cell technologies enable us to explore the mechanism of plant regeneration at single-cell resolution. RESULTS: We generate a high-resolution single-cell transcriptomic landscape of hypocotyl tissue from the highly regenerable cotton genotype Jin668 and the recalcitrant TM-1. We identify nine putative cell clusters and 23 cluster-specific marker genes for both cultivars. We find that the primary vascular cell is the major cell type that undergoes cell fate transition in response to external stimulation. Further developmental trajectory and gene regulatory network analysis of these cell clusters reveals that a total of 41 hormone response-related genes, including LAX2, LAX1, and LOX3, exhibit different expression patterns in the primary xylem and cambium region of Jin668 and TM-1. We also identify novel genes, including CSEF, PIS1, AFB2, ATHB2, PLC2, and PLT3, that are involved in regeneration. We demonstrate that LAX2, LAX1 and LOX3 play important roles in callus proliferation and plant regeneration by CRISPR/Cas9 editing and overexpression assay. CONCLUSIONS: This study provides novel insights on the role of the regulatory network in cell fate transition and reprogramming during plant regeneration driven by somatic embryogenesis.

N6-methyladenosine RNA modification regulates cotton drought response in a Ca2+ and ABA-dependent manner.

N6 -methyladenosine (m6 A) is the most prevalent internal modification present in mRNAs, and is considered to participate in a range of developmental and biological processes. Drought response is highly regulated at the genomic, transcriptional and post-transcriptional levels. However, the biological function and regulatory mechanism of m6 A modification in the drought stress response is still poorly understood. We generated a transcriptome-wide m6 A map using drought-resistant and drought-sensitive varieties of cotton under different water deficient conditions to uncover patterns of m6 A methylation in cotton response to drought stress. The results reveal that m6 A represents a common modification and exhibit dramatic changes in distribution during drought stress. More 5'UTR m6 A was deposited in the drought-resistant variety and was associated with a positive effect on drought resistance by regulating mRNA abundance. Interestingly, we observed that increased m6 A abundance was associated with increased mRNA abundance under drought, contributing to drought resistance, and vice versa. The demethylase GhALKBH10B was found to decrease m6 A levels, facilitating the mRNA decay of ABA signal-related genes (GhZEP, GhNCED4 and GhPP2CA) and Ca2+ signal-related genes (GhECA1, GhCNGC4, GhANN1 and GhCML13), and mutation of GhALKBH10B enhanced drought resistance at seedling stage in cotton. Virus-induced gene silencing (VIGS) of two Ca2+ -related genes, GhECA1 and GhCNGC4, reduced drought resistance with the decreased m6 A enrichment on silenced genes in cotton. Collectively, we reveal a novel mechanism of post-transcriptional modification involved in affecting drought response in cotton, by mediating m6 A methylation on targeted transcripts in the ABA and Ca2+ signalling transduction pathways.

MERISTEM-DEFECTIVE regulates the balance between stemness and differentiation in the root meristem through RNA splicing control.

Plants respond to environmental stresses through controlled stem cell maintenance and meristem activity. One level of gene regulation is RNA alternative splicing. However, the mechanistic link between stress, meristem function and RNA splicing is poorly understood. The MERISTEM-DEFECTIVE (MDF) Arabidopsis gene encodes an SR-related family protein, required for meristem function and leaf vascularization, and is the likely orthologue of the human SART1 and yeast Snu66 splicing factors. MDF is required for the correct splicing and expression of key transcripts associated with root meristem function. We identified RSZ33 and ACC1, both known to regulate cell patterning, as splicing targets required for MDF function in the meristem. MDF expression is modulated by osmotic and cold stress, associated with differential splicing and specific isoform accumulation and shuttling between nucleus and cytosol, and acts in part via a splicing target SR34. We propose a model in which MDF controls splicing in the root meristem to promote stemness and to repress stress response, cell differentiation and cell death pathways.

GhWRKY41 forms a positive feedback regulation loop and increases cotton defence response against Verticillium dahliae by regulating phenylpropanoid metabolism.

Despite the established significance of WRKY proteins and phenylpropanoid metabolism in plant immunity, how WRKY proteins modulate aspects of the phenylpropanoid pathway remains undetermined. To understand better the role of WRKY proteins in plant defence, we identified a cotton (Gossypium hirsutum) protein, GhWRKY41, that is, universally and rapidly induced in three disease-resistant cotton cultivars following inoculation with the plant pathogenic fungus, Verticillium dahliae. We show that overexpression of GhWRKY41 in transgenic cotton and Arabidopsis enhances resistance to V. dahliae, while knock-down increases cotton more susceptibility to the fungus. GhWRKY41 physically interacts with itself and directly activates its own transcription. A genome-wide chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq), in combination with RNA sequencing (RNA-seq) analyses, revealed that 43.1% of GhWRKY41-binding genes were up-regulated in cotton upon inoculation with V. dahliae, including several phenylpropanoid metabolism master switches, receptor kinases, and disease resistance-related proteins. We also show that GhWRKY41 homodimer directly activates the expression of GhC4H and Gh4CL, thereby modulating the accumulation of lignin and flavonoids. This finding expands our understanding of WRKY-WRKY protein interactions and provides important insights into the regulation of the phenylpropanoid pathway in plant immune responses by a WRKY protein.

High-resolution sequencing of nine elite upland cotton cultivars uncovers genic variations and breeding improvement targets.

Structural variations (SVs) are critical factors affecting genome evolution and important traits. However, identification results and functional analyses of SVs in upland cotton are rare. Here, based on the genetic relationships, breeding history and cumulative planting area of upland cotton in China, nine predominant cultivars from the past 60 years (1950s-2010s) were selected for long read sequencing to uncover genic variations and breeding improvement targets for this crop. Based on the ZM24 reference genome, 0.88-1.47 × 104 SVs per cultivar were identified, and an SV set was constructed. SVs affected the expression of a large number of genes during fiber elongation, and a transposable element insertion resulted in the glandless phenotype in upland cotton. Six widespread inversions were identified based on nine draft genomes and high-throughput chromosome conformation capture data. Multiple haplotype blocks that were always associated with aggregated SVs were demonstrated to play a pivotal role in the agronomic traits of upland cotton and drove its adaptation to the northern planting region. Exotic introgression was the source of these haplotype blocks and increased the genetic diversity of upland cotton. Our results enrich the genome resources of upland cotton, and the identified SVs will promote genetic and breeding research in cotton.

The GhMAP3K62-GhMKK16-GhMPK32 kinase cascade regulates drought tolerance by activating GhEDT1-mediated ABA accumulation in cotton.

INTRODUCTION: Drought is the principal abiotic stress that severely impacts cotton (Gossypium hirsutum) growth and productivity. Upon sensing drought, plants activate stress-related signal transduction pathways, including ABA signal and mitogen-activated protein kinase (MAPK) cascade. However, as the key components with the fewest members in the MAPK cascade, the function and regulation of GhMKKs need to be elucidated. In addition, the relationship between MAPK module and the ABA core signaling pathway remains incompletely understood. OBJECTIVE: Here we aim to elucidate the molecular mechanism of cotton response to drought, with a focus on mitogen-activated protein kinase (MAPK) cascades activating ABA signaling. METHODS: Biochemical, molecular and genetic analysis were used to study the GhMAP3K62-GhMKK16-GhMPK32-GhEDT1 pathway genes. RESULTS: A nucleus- and membrane-localized MAPK cascade pathway GhMAP3K62-GhMKK16-GhMPK32, which targets and phosphorylates the nuclear-localized transcription factor GhEDT1, to activate downstream GhNCED3 to mediate ABA-induced stomatal closure and drought response was characterized in cotton. Overexpression of GhMKK16 promotes ABA accumulation, and enhances drought tolerance via regulating stomatal closure under drought stress. Conversely, RNAi-mediated knockdown of GhMKK16 expression inhibits ABA accumulation, and reduces drought tolerance. Virus-induced gene silencing (VIGS)-mediated knockdown of either GhMAP3K62, GhMPK32 or GhEDT1 expression represses ABA accumulation and reduces drought tolerance through inhibiting stomatal closure. Expression knockdown of GhMPK32 or GhEDT1 in GhMKK16-overexpressing cotton reinstates ABA content and stomatal opening-dependent drought sensitivity to wild type levels. GhEDT1 could bind to the HD boxes in the promoter of GhNCED3 to activate its expression, resulting in ABA accumulation. We propose that the MAPK cascade GhMAP3K62-GhMKK16-GhMPK32 pathway functions on drought response through ABA-dependent stomatal movement in cotton.

The complex genome and adaptive evolution of polyploid Chinese pepper (Zanthoxylum armatum and Zanthoxylum bungeanum).

Zanthoxylum armatum and Zanthoxylum bungeanum, known as 'Chinese pepper', are distinguished by their extraordinary complex genomes, phenotypic innovation of adaptive evolution and species-special metabolites. Here, we report reference-grade genomes of Z. armatum and Z. bungeanum. Using high coverage sequence data and comprehensive assembly strategies, we derived 66 pseudochromosomes comprising 33 homologous phased groups of two subgenomes, including autotetraploid Z. armatum. The genomic rearrangements and two whole-genome duplications created large (~4.5 Gb) complex genomes with a high ratio of repetitive sequences (>82%) and high chromosome number (2n = 4x = 132). Further analysis of the high-quality genomes shed lights on the genomic basis of involutional reproduction, allomones biosynthesis and adaptive evolution in Chinese pepper, revealing a high consistent relationship between genomic evolution, environmental factors and phenotypic innovation. Our study provides genomic resources and new insights for investigating diversification and phenotypic innovation in Chinese pepper, with broader implications for the protection of plants under severe environmental changes.

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.

GhTCE1-GhTCEE1 dimers regulate transcriptional reprogramming during wound-induced callus formation in cotton.

Wounded plant cells can form callus to seal the wound site. Alternatively, wounding can cause adventitious organogenesis or somatic embryogenesis. These distinct developmental pathways require specific cell fate decisions. Here, we identify GhTCE1, a basic helix-loop-helix family transcription factor, and its interacting partners as a central regulatory module of early cell fate transition during in vitro dedifferentiation of cotton (Gossypium hirsutum). RNAi- or CRISPR/Cas9-mediated loss of GhTCE1 function resulted in excessive accumulation of reactive oxygen species (ROS), arrested callus cell elongation, and increased adventitious organogenesis. In contrast, GhTCE1-overexpressing tissues underwent callus cell growth, but organogenesis was repressed. Transcriptome analysis revealed that several pathways depend on proper regulation of GhTCE1 expression, including lipid transfer pathway components, ROS homeostasis, and cell expansion. GhTCE1 bound to the promoters of the target genes GhLTP2 and GhLTP3, activating their expression synergistically, and the heterodimer TCE1-TCEE1 enhances this activity. GhLTP2- and GhLTP3-deficient tissues accumulated ROS and had arrested callus cell elongation, which was restored by ROS scavengers. These results reveal a unique regulatory network involving ROS and lipid transfer proteins, which act as potential ROS scavengers. This network acts as a switch between unorganized callus growth and organized development during in vitro dedifferentiation of cotton cells.