Single-cell transcriptome atlas identified novel regulators for pigment gland morphogenesis in cotton

GhHAM regulates GoPGF-dependent gland development and contributes to broad-spectrum pest resistance in cotton

Cotton is a globally cultivated crop, producing 87% of the natural fiber used in the global textile industry. The pigment glands, unique to cotton and its relatives, serve as a defense structure against pests and pathogens. However, the molecular mechanism underlying gland formation and the specific role of pigment glands in cotton's pest defense are still not well understood. In this study, we cloned a gland-related transcription factor GhHAM and generated the GhHAM knockout mutant using CRISPR/Cas9. Phenotypic observations, transcriptome analysis, and promoter-binding experiments revealed that GhHAM binds to the promoter of GoPGF, regulating pigment gland formation in cotton's multiple organs via the GoPGF-GhJUB1 module. The knockout of GhHAM significantly reduced gossypol production and increased cotton's susceptibility to pests in the field. Feeding assays demonstrated that more than 80% of the cotton bollworm larvae preferred ghham over the wild type. Furthermore, the ghham mutants displayed shorter cell length and decreased gibberellins (GA) production in the stem. Exogenous application of GA3 restored stem cell elongation but not gland formation, thereby indicating that GhHAM controls gland morphogenesis independently of GA. Our study sheds light on the functional differentiation of HAM proteins among plant species, highlights the significant role of pigment glands in influencing pest feeding preference, and provides a theoretical basis for breeding pest-resistant cotton varieties to address the challenges posed by frequent outbreaks of pests.

Global dynamics and cytokinin participation of salt gland development trajectory in recretohalophyte Limonium bicolor

Salt gland is an epidermal Na+ secretory structure that enhances salt resistance in the recretohalophyte sea lavender (Limonium bicolor). To elucidate the salt gland development trajectory and related molecular mechanisms, we performed single-cell RNA sequencing of L. bicolor protoplasts from young leaves at salt gland initiation and differentiation stages. Dimensionality reduction analyses defined 19 transcriptionally distinct cell clusters, which were assigned into 4 broad populations-promeristem, epidermis, mesophyll, and vascular tissue-verified by in situ hybridization. Cytokinin was further proposed to participate in salt gland development by the expression patterns of related genes and cytological evidence. By comparison analyses of Single-cell RNA sequencing with exogenous application of 6-benzylaminopurine, we delineated 5 salt gland development-associated subclusters and defined salt gland-specific differentiation trajectories from Subclusters 8, 4, and 6 to Subcluster 3 and 1. Additionally, we validated the participation of TRIPTYCHON and the interacting protein Lb7G34824 in salt gland development, which regulated the expression of cytokinin metabolism and signaling-related genes such as GLABROUS INFLORESCENCE STEMS 2 to maintain cytokinin homeostasis during salt gland development. Our results generated a gene expression map of young leaves at single-cell resolution for the comprehensive investigation of salt gland determinants and cytokinin participation that helps elucidate cell fate determination during epidermis formation and evolution in recretohalophytes.

Regulation of Glandular Size and Phytoalexin Biosynthesis by a Negative Feedback Loop in Cotton

Plants have evolved diverse defense mechanisms encompassing physical and chemical barriers. Cotton pigment glands are known for containing various defense metabolites, but the precise regulation of gland size to modulate defense compound levels remains enigmatic. Here, it is discovered that the VQ domain-containing protein JAVL negatively regulates pigment gland size and the biosynthesis of defense compounds, while the MYC2-like transcription factor GoPGF has the opposite effect. Notably, GoPGF directly activates the expression of JAVL, whereas JAVL suppresses GoPGF transcription, establishing a negative feedback loop that maintains the expression homeostasis between GoPGF and JAVL. Furthermore, it is observed that JAVL negatively regulates jasmonate levels by inhibiting the expression of jasmonate biosynthetic genes and interacting with GoPGF to attenuate its activation effects, thereby maintaining homeostatic regulation of jasmonate levels. The increased expression ratio of GoPGF to JAVL leads to enlarged pigment glands and elevated jasmonates and defense compounds, enhancing insect and pathogen resistance in cotton. These findings unveil a new mechanism for regulating gland size and secondary metabolites biosynthesis, providing innovative strategies for strengthening plant defense.

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

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

In vivo maternal haploid induction system in cotton

The ghdmp mutant of cotton, generated through the CRISPR system, exhibits a haploid induction rate of 1.06% in F1 progeny as the haploid inducer line.

The dual role of GoPGF reveals that the pigment glands are synthetic sites of gossypol in aerial parts of cotton

Gossypol and the related terpenoids are stored in the pigment gland to protect cotton plants from biotic stresses, but little is known about the synthetic sites of these metabolites. Here, we showed that GoPGF, a key gene regulating gland formation, was expressed in gland cells and roots. The chromatin immunoprecipitation sequencing (ChIP-seq) analysis demonstrated that GoPGF targets GhJUB1 to regulate gland morphogenesis. RNA-sequencing (RNA-seq) showed high accumulation of gossypol biosynthetic genes in gland cells. Moreover, integrated analysis of the ChIP-seq and RNA-seq data revealed that GoPGF binds to the promoter of several gossypol biosynthetic genes. The cotton callus overexpressing GoPGF had dramatically increased the gossypol levels, indicating that GoPGF can directly activate the biosynthesis of gossypol. In addition, the gopgf mutant analysis revealed the existence of both GoPGF-dependent and -independent regulation of gossypol production in cotton roots. Our study revealed that the pigment glands are synthetic sites of gossypol in aerial parts of cotton and that GoPGF plays a dual role in regulating gland morphogenesis and gossypol biosynthesis. The study provides new insights for exploring the complex relationship between glands and the metabolites they store in cotton and other plant species.

Cell type-specific cytonuclear coevolution in three allopolyploid plant species

Cytonuclear disruption may accompany allopolyploid evolution as a consequence of the merger of different nuclear genomes in a cellular environment having only one set of progenitor organellar genomes. One path to reconcile potential cytonuclear mismatch is biased expression for maternal gene duplicates (homoeologs) encoding proteins that target to plastids and/or mitochondria. Assessment of this transcriptional form of cytonuclear coevolution at the level of individual cells or cell types remains unexplored. Using single-cell (sc-) and single-nucleus (sn-) RNAseq data from eight tissues in three allopolyploid species, we characterized cell type-specific variations of cytonuclear coevolutionary homoeologous expression and demonstrated the temporal dynamics of expression patterns across development stages during cotton fiber development. Our results provide unique insights into transcriptional cytonuclear coevolution in plant allopolyploids at the single-cell level.