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

Transmembrane formins as active cargoes of membrane trafficking

Formins are a large, evolutionarily old family of cytoskeletal regulators whose roles include actin capping and nucleation, as well as modulation of microtubule dynamics. The plant class I formin clade is characterized by a unique domain organization, as most of its members are transmembrane proteins with possible cell wall-binding motifs exposed to the extracytoplasmic space-a structure that appears to be a synapomorphy of the plant kingdom. While such transmembrane formins are traditionally considered mainly as plasmalemma-localized proteins contributing to the organization of the cell cortex, we review, from a cell biology perspective, the growing evidence that they can also, at least temporarily, reside (and in some cases also function) in endomembranes including secretory and endocytotic pathway compartments, the endoplasmic reticulum, the nuclear envelope, and the tonoplast. Based on this evidence, we propose that class I formins may thus serve as 'active cargoes' of membrane trafficking-membrane-embedded proteins that modulate the fate of endo- or exocytotic compartments while being transported by them.

Advances in Single-Cell Transcriptome Sequencing and Spatial Transcriptome Sequencing in Plants

"Omics" typically involves exploration of the structure and function of the entire composition of a biological system at a specific level using high-throughput analytical methods to probe and analyze large amounts of data, including genomics, transcriptomics, proteomics, and metabolomics, among other types. Genomics characterizes and quantifies all genes of an organism collectively, studying their interrelationships and their impacts on the organism. However, conventional transcriptomic sequencing techniques target population cells, and their results only reflect the average expression levels of genes in population cells, as they are unable to reveal the gene expression heterogeneity and spatial heterogeneity among individual cells, thus masking the expression specificity between different cells. Single-cell transcriptomic sequencing and spatial transcriptomic sequencing techniques analyze the transcriptome of individual cells in plant or animal tissues, enabling the understanding of each cell's metabolites and expressed genes. Consequently, statistical analysis of the corresponding tissues can be performed, with the purpose of achieving cell classification, evolutionary growth, and physiological and pathological analyses. This article provides an overview of the research progress in plant single-cell and spatial transcriptomics, as well as their applications and challenges in plants. Furthermore, prospects for the development of single-cell and spatial transcriptomics are proposed.

Genome-wide analyses of member identification, expression pattern, and protein-protein interaction of EPF/EPFL gene family in Gossypium

BACKGROUND

Epidermal patterning factor / -like (EPF/EPFL) gene family encodes a class of cysteine-rich secretory peptides, which are widelyfound in terrestrial plants.Multiple studies has indicated that EPF/EPFLs might play significant roles in coordinating plant development and growth, especially as the morphogenesis processes of stoma, awn, stamen, and fruit skin. However, few research on EPF/EPFL gene family was reported in Gossypium.

RESULTS

We separately identified 20 G. raimondii, 24 G. arboreum, 44 G. hirsutum, and 44 G. barbadense EPF/EPFL genes in the 4 representative cotton species, which were divided into four clades together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 17 Selaginella moellendorffii ones based on their evolutionary relationships. The similar gene structure and common motifs indicated the high conservation among the EPF/EPFL members, while the uneven distribution in chromosomes implied the variability during the long-term evolutionary process. Hundreds of collinearity relationships were identified from the pairwise comparisons of intraspecifc and interspecific genomes, which illustrated gene duplication might contribute to the expansion of cotton EPF/EPFL gene family. A total of 15 kinds of cis-regulatory elements were predicted in the promoter regions, and divided into three major categories relevant to the biological processes of development and growth, plant hormone response, and abiotic stress response. Having performing the expression pattern analyses with the basic of the published RNA-seq data, we found most of GhEPF/EPFL and GbEPF/EPFL genes presented the relatively low expression levels among the 9 tissues or organs, while showed more dramatically different responses to high/low temperature and salt or drought stresses. Combined with transcriptome data of developing ovules and fibers and quantitative Real-time PCR results (qRT-PCR) of 15 highly expressed GhEPF/EPFL genes, it could be deduced that the cotton EPF/EPFL genes were closely related with fiber development. Additionally, the networks of protein-protein interacting among EPF/EPFLs concentrated on the cores of GhEPF1 and GhEPF7, and thosefunctional enrichment analyses indicated that most of EPF/EPFLs participate in the GO (Gene Ontology) terms of stomatal development and plant epidermis development, and the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways of DNA or base excision repair.

CONCLUSION

Totally, 132 EPF/EPFL genes were identified for the first time in cotton, whose bioinformatic analyses of cis-regulatory elements and expression patterns combined with qRT-PCR experiments to prove the potential functions in the biological processes of plant growth and responding to abiotic stresses, specifically in the fiber development. These results not only provide comprehensive and valuable information for cotton EPF/EPFL gene family, but also lay solid foundation for screening candidate EPF/EPFL genes in further cotton breeding.

Spatiotemporal transcriptomic atlas of rhizome formation in Oryza longistaminata

Rhizomes are modified stems that grow underground and produce new individuals genetically identical to the mother plant. Recently, a breakthrough has been made in efforts to convert annual grains into perennial ones by utilizing wild rhizomatous species as donors, yet the developmental biology of this organ is rarely studied. Oryza longistaminata, a wild rice species featuring strong rhizomes, provides a valuable model for exploration of rhizome development. Here, we first assembled a double-haplotype genome of O. longistaminata, which displays a 48-fold improvement in contiguity compared to the previously published assembly. Furthermore, spatiotemporal transcriptomics was performed to obtain the expression profiles of different tissues in O. longistaminata rhizomes and tillers. Two spatially reciprocal cell clusters, the vascular bundle 2 cluster and the parenchyma 2 cluster, were determined to be the primary distinctions between the rhizomes and tillers. We also captured meristem initiation cells in the sunken area of parenchyma located at the base of internodes, which is the starting point for rhizome initiation. Trajectory analysis further indicated that the rhizome is regenerated through de novo generation. Collectively, these analyses revealed a spatiotemporal transcriptional transition underlying the rhizome initiation, providing a valuable resource for future perennial crop breeding.

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.

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.

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Co-localization and analysis of miR477b with fiber length quantitative trait loci in cotton

Cotton is an important cash crop for the textile industry. However, the understanding of natural genetic variation of fiber elongation in relation to miRNA is lacking. A miRNA gene (miR477b) was found to co-localize with a previously mapped fiber length (FL) quantitative trait locus (QTL). The miR477b was differentially expressed during fiber elongation between two backcross inbred lines (BILs) differing in FL and its precursor sequences. Bioinformatics and qRT-PCR analysis were further used to analyse the miRNA genes, which could produce mature miR477b. Cotton plants with virus-induced gene silencing (VIGS) constructs to over-express the allele of miR477b from the BIL with longer fibers had significantly longer fibers as compared with negative control plants, while the VIGS plants with suppressed miRNA expression had significantly shorter fibers. The expression level of the target gene (DELLA) and related genes (RDL1 and EXPA1 for DELLA through HOX3 protein) in the two BILs and/or the VIGS plants were generally congruent, as expected. This report represents one of the first comprehensive studies to integrate QTL linkage mapping and physical mapping of small RNAs with both small and mRNA transcriptome analysis, followed by VIGS, to identify candidate small RNA genes affecting the natural variation of fiber elongation in cotton.

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

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

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Characterization and expression analysis of GLABRA3 (GL3) genes in cotton: insights into trichome development and hormonal regulation

BACKGROUND

GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) genes encode a typical helix-loop-helix (bHLH) transcription factors that primarily regulate trichome branching and root hair development, DNA endoreduplication, trichoblast size, and stomatal formation. The functions of GL3 genes in cotton crop have been poorly characterized. In this study, we performed comprehensive genome-wide scans for GL3 and EGL3 homologs to enhance our comprehension of their potential roles in trichome and fiber development in cotton crop.

METHODS AND RESULTS

Our findings paraded that Gossypium hirsutum and G. barbadense have 6 GL3s each, unevenly distributed on 4 chromosomes whereas, G. arboreum, and G. raimondii have 3 GL3s each, unevenly distributed on 2 chromosomes. Gh_A08G2088 and Gb_A09G2187, despite having the same bHLH domain as the other GL3 genes, were excluded due to remarkable short sequences and limited number of motifs, indicating a lack of potential functional activity. The phylogenetic analysis categorized remaining 16 GL3s into three subfamilies (Group I-III) closely related to A. thaliana. The 16 GL3s have complete bHLH domain, encompassing 590-631 amino acids, with molecular weights (MWs) ranging from 65.92 to 71.36 kDa. Within each subfamily GL3s depicted shared similar gene structures and motifs, indicating conserved characteristics within respective groups. Promoter region analysis revealed 27 cis-acting elements, these elements were responsive to salicylic acid, abscisic acid (ABA), methyl jasmonate (MeJA), and gibberellin. The expression of GL3 genes was analyzed across 12 tissues in both G. barbadense and G. hirsutum using the publicly available RNA-seq data. Among GL3s, Gb_D11G0219, Gb_D11G0214, and Gb_D08G2182, were identified as relatively highly expressed across different tissues, consequently selected for hormone treatment and expression validation in G. barbadense. RT-qPCR results demonstrated significant alterations in the expression levels of Gb_D11G0219 and Gb_D11G0214 following MeJA, GA, and ABA treatment. Subcellular localization prediction revealed that most GL3 proteins were predominantly expressed in the nucleus, while a few were localized in the cytoplasm and chloroplasts.

CONCLUSIONS

In summary, this study lays the foundation for subsequent functional validation of GL3 genes by identifying hormonal regulation patterns and probable sites of action in cotton trichome formation and fiber development. The results stipulate a rationale to elucidate the roles and regulatory mechanisms of GL3 genes in the intricate process of cotton fibre and trichome development.

Advances in the Application of Single-Cell Transcriptomics in Plant Systems and Synthetic Biology

Plants are complex systems hierarchically organized and composed of various cell types. To understand the molecular underpinnings of complex plant systems, single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for revealing high resolution of gene expression patterns at the cellular level and investigating the cell-type heterogeneity. Furthermore, scRNA-seq analysis of plant biosystems has great potential for generating new knowledge to inform plant biosystems design and synthetic biology, which aims to modify plants genetically/epigenetically through genome editing, engineering, or re-writing based on rational design for increasing crop yield and quality, promoting the bioeconomy and enhancing environmental sustainability. In particular, data from scRNA-seq studies can be utilized to facilitate the development of high-precision Build-Design-Test-Learn capabilities for maximizing the targeted performance of engineered plant biosystems while minimizing unintended side effects. To date, scRNA-seq has been demonstrated in a limited number of plant species, including model plants (e.g., Arabidopsis thaliana), agricultural crops (e.g., Oryza sativa), and bioenergy crops (e.g., Populus spp.). It is expected that future technical advancements will reduce the cost of scRNA-seq and consequently accelerate the application of this emerging technology in plants. In this review, we summarize current technical advancements in plant scRNA-seq, including sample preparation, sequencing, and data analysis, to provide guidance on how to choose the appropriate scRNA-seq methods for different types of plant samples. We then highlight various applications of scRNA-seq in both plant systems biology and plant synthetic biology research. Finally, we discuss the challenges and opportunities for the application of scRNA-seq in plants.

A comprehensive review on Gossypium hirsutum resistance against cotton leaf curl virus

Cotton (Gossypium hirsutum L.) is a significant fiber crop. Being a major contributor to the textile industry requires continuous care and attention. Cotton is subjected to various biotic and abiotic constraints. Among these, biotic factors including cotton leaf curl virus (CLCuV) are dominant. CLCuV is a notorious disease of cotton and is acquired, carried, and transmitted by the whitefly (Bemisia tabaci). A cotton plant affected with CLCuV may show a wide range of symptoms such as yellowing of leaves, thickening of veins, upward or downward curling, formation of enations, and stunted growth. Though there are many efforts to protect the crop from CLCuV, long-term results are not yet obtained as CLCuV strains are capable of mutating and overcoming plant resistance. However, systemic-induced resistance using a gene-based approach remained effective until new virulent strains of CLCuV (like Cotton Leaf Curl Burewala Virus and others) came into existence. Disease control by biological means and the development of CLCuV-resistant cotton varieties are in progress. In this review, we first discussed in detail the evolution of cotton and CLCuV strains, the transmission mechanism of CLCuV, the genetic architecture of CLCuV vectors, and the use of pathogen and nonpathogen-based approaches to control CLCuD. Next, we delineate the uses of cutting-edge technologies like genome editing (with a special focus on CRISPR-Cas), next-generation technologies, and their application in cotton genomics and speed breeding to develop CLCuD resistant cotton germplasm in a short time. Finally, we delve into the current obstacles related to cotton genome editing and explore forthcoming pathways for enhancing precision in genome editing through the utilization of advanced genome editing technologies. These endeavors aim to enhance cotton's resilience against CLCuD.

Dissecting the major genetic components underlying cotton lint development

Numerous genetic loci and several functionally characterized genes have been linked to determination of lint percentage (lint%), one of the most important cotton yield components, but we still know little about the major genetic components underlying lint%. Here, we first linked the genetic loci containing MYB25-like_At and HD1_At to the fiberless seed trait of 'SL1-7-1' and found that MYB25-like_At and HD1_At were very lowly expressed in 'SL1-7-1' ovules during fiber initiation. We then dissected the genetic components involved in determination of lint% using segregating populations derived from crosses of fuzzless mutants and intermediate segregants with different lint%, which not only confirmed the HD1_At locus but identified the HD1_Dt locus as being the major genetic components contributing to fiber initiation and lint%. The segregating populations also allowed us to evaluate the relative contributions of MYB25-like_At, MYB25-like_Dt, HD1_At, and HD1_Dt to lint%. Haplotype analysis of an Upland cotton (Gossypium hirsutum) population with 723 accessions (including 81 fuzzless seed accessions) showed that lint% of the accessions with the LP allele (higher lint%) at MYB25-like_At, MYB25-like_Dt, or HD1_At was significantly higher than that with the lp allele (lower lint%). The lint% of the Upland cotton accessions with 3 or 4 LP alleles at MYB25-like and HD1 was significantly higher than that with 2 LP alleles. The results prompted us to propose a strategy for breeding high-yielding cotton varieties, i.e. pyramiding the LP alleles of MYB25-like and HD1 with new lint% LP alleles without negative impact on seed size and fiber quality.

scPlantDB: a comprehensive database for exploring cell types and markers of plant cell atlases

Recent advancements in single-cell RNA sequencing (scRNA-seq) technology have enabled the comprehensive profiling of gene expression patterns at the single-cell level, offering unprecedented insights into cellular diversity and heterogeneity within plant tissues. In this study, we present a systematic approach to construct a plant single-cell database, scPlantDB, which is publicly available at https://biobigdata.nju.edu.cn/scplantdb. We integrated single-cell transcriptomic profiles from 67 high-quality datasets across 17 plant species, comprising approximately 2.5 million cells. The data underwent rigorous collection, manual curation, strict quality control and standardized processing from public databases. scPlantDB offers interactive visualization of gene expression at the single-cell level, facilitating the exploration of both single-dataset and multiple-dataset analyses. It enables systematic comparison and functional annotation of markers across diverse cell types and species while providing tools to identify and compare cell types based on these markers. In summary, scPlantDB serves as a comprehensive database for investigating cell types and markers within plant cell atlases. It is a valuable resource for the plant research community.

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.

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.

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.

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.

Application of single-cell multi-omics approaches in horticulture research

Cell heterogeneity shapes the morphology and function of various tissues and organs in multicellular organisms. Elucidation of the differences among cells and the mechanism of intercellular regulation is essential for an in-depth understanding of the developmental process. In recent years, the rapid development of high-throughput single-cell transcriptome sequencing technologies has influenced the study of plant developmental biology. Additionally, the accuracy and sensitivity of tools used to study the epigenome and metabolome have significantly increased, thus enabling multi-omics analysis at single-cell resolution. Here, we summarize the currently available single-cell multi-omics approaches and their recent applications in plant research, review the single-cell based studies in fruit, vegetable, and ornamental crops, and discuss the potential of such approaches in future horticulture research.

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.

Integrated mass spectrometry imaging and single-cell transcriptome atlas strategies provide novel insights into taxoid biosynthesis and transport in Taxus mairei stems

Taxol, which is a widely used important chemotherapeutic agent, was originally isolated from Taxus stem barks. However, little is known about the precise distribution of taxoids and the transcriptional regulation of taxoid biosynthesis across Taxus stems. Here, we used MALDI-IMS analysis to visualize the taxoid distribution across Taxus mairei stems and single-cell RNA sequencing to generate expression profiles. A single-cell T. mairei stem atlas was created, providing a spatial distribution pattern of Taxus stem cells. Cells were reordered using a main developmental pseudotime trajectory which provided temporal distribution patterns in Taxus stem cells. Most known taxol biosynthesis-related genes were primarily expressed in epidermal, endodermal, and xylem parenchyma cells, which caused an uneven taxoid distribution across T. mairei stems. We developed a single-cell strategy to screen novel transcription factors (TFs) involved in taxol biosynthesis regulation. Several TF genes, such as endodermal cell-specific MYB47 and xylem parenchyma cell-specific NAC2 and bHLH68, were implicated as potential regulators of taxol biosynthesis. Furthermore, an ATP-binding cassette family transporter gene, ABCG2, was proposed as a potential taxoid transporter candidate. In summary, we generated a single-cell Taxus stem metabolic atlas and identified molecular mechanisms underpinning the cell-specific transcriptional regulation of the taxol biosynthesis pathway.

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.

Heterogeneity of cancer-associated fibroblasts in head and neck squamous cell carcinoma

Cancer-associated fibroblasts (CAFs) consist of heterogeneous cellular populations that contribute critical roles in head and neck squamous cell carcinoma (HNSCC). A series of computer-aided analyses were performed to determine various aspects of CAFs in HNSCC, including their cellular heterogeneity, prognostic value, relationship with immune suppression and immunotherapeutic response, intercellular communication, and metabolic activity. The prognostic significance of CKS2+ CAFs was verified using immunohistochemistry. Our findings revealed that fibroblasts group demonstrated prognostic significance, with the CKS2+ subset of inflammatory CAFs (iCAFs) exhibiting a significant correlation with unfavorable prognosis and being localized in close proximity to cancer cells. Patients with a high infiltration of CKS2+ CAFs had a poor overall survival rate. There is a negative correlation between CKS2+ iCAFs and cytotoxic CD8+ T cells and natural killer (NK) cells, while a positive correlation was found with exhausted CD8+ T cells. Additionally, patients in Cluster 3, characterized by a high proportion of CKS2+ iCAFs, and patients in Cluster 2, characterized by a high proportion of CKS2- iCAFs and CENPF-/MYLPF- myofibroblastic CAFs (myCAFs), did not exhibit significant immunotherapeutic responses. Moreover, close interactions was confirmed to exist between cancer cells and CKS2+ iCAFs/ CENPF+ myCAFs. Furthermore, CKS2+ iCAFs demonstrated the highest level of metabolic activity. In summary, our study enhances the understanding of the heterogeneity of CAFs and provided insights into improving the efficacy of immunotherapies and prognostic accuracy for HNSCC patients.

Genetic regulation of nitrogen use efficiency in Gossypium spp

Cotton (Gossypium spp.) is the most important fibre crop, with desirable characteristics preferred for textile production. Cotton fibre output relies heavily on nitrate as the most important source of inorganic nitrogen (N). However, nitrogen dynamics in extreme environments limit plant growth and lead to yield loss and pollution. Therefore, nitrogen use efficiency (NUE), which involves the utilisation of the 'right rate', 'right source', 'right time', and 'right place' (4Rs), is key for efficient N management. Recent omics techniques have genetically improved NUE in crops. We herein highlight the mechanisms of N uptake and assimilation in the vegetative and reproductive branches of the cotton plant while considering the known and unknown regulatory factors. The phylogenetic relationships among N transporters in four Gossypium spp. have been reviewed. Further, the N regulatory genes that participate in xylem transport and phloem loading are also discussed. In addition, the functions of microRNAs and transcription factors in modulating the expression of target N regulatory genes are highlighted. Overall, this review provides a detailed perspective on the complex N regulatory mechanism in cotton, which would accelerate the research toward improving NUE in crops.

Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis

BACKGROUND

Homoeologs are defined as homologous genes resulting from allopolyploidy. Bread wheat, Triticum aestivum, is an allohexaploid species with many homoeologs. Homoeolog expression bias, referring to the relative contribution of homoeologs to the transcriptome, is critical for determining the traits that influence wheat growth and development. Asymmetric transcription of homoeologs has been so far investigated in a tissue or organ-specific manner, which could be misleading due to a mixture of cell types.

RESULTS

Here, we perform single nuclei RNA sequencing and ATAC sequencing of wheat root to study the asymmetric gene transcription, reconstruct cell differentiation trajectories and cell-type-specific gene regulatory networks. We identify 22 cell types. We then reconstruct cell differentiation trajectories that suggest different origins between epidermis/cortex and endodermis, distinguishing bread wheat from Arabidopsis. We show that the ratio of asymmetrically transcribed triads varies greatly when analyzing at the single-cell level. Hub transcription factors determining cell type identity are also identified. In particular, we demonstrate that TaSPL14 participates in vasculature development by regulating the expression of BAM1. Combining single-cell transcription and chromatin accessibility data, we construct the pseudo-time regulatory network driving root hair differentiation. We find MYB3R4, REF6, HDG1, and GATAs as key regulators in this process.

CONCLUSIONS

Our findings reveal the transcriptional landscape of root organization and asymmetric gene transcription at single-cell resolution in polyploid wheat.

Single-cell analysis opens a goldmine for plant functional studies

Functional studies in biology require the identification of genes, regulatory elements, and networks, followed by a deep understanding of how they orchestrate to specify cell types, mediate signaling, and respond to internal and external cues over evolutionary timescales. Advances in single-cell analysis have enabled biologists to tackle these questions at the resolution of the individual cell. Here, we highlight recent studies in plants that have embraced single-cell analyses to facilitate functional studies. This review will provide guidance and perspectives for incorporating these advanced approaches in plant research for the coming decades.

Recent Progress and Future Prospect of CRISPR/Cas-Derived Transcription Activation (CRISPRa) System in Plants

Genome editing technology has become one of the hottest research areas in recent years. Among diverse genome editing tools, the Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated proteins system (CRISPR/Cas system) has exhibited the obvious advantages of specificity, simplicity, and flexibility over any previous genome editing system. In addition, the emergence of Cas9 mutants, such as dCas9 (dead Cas9), which lost its endonuclease activity but maintains DNA recognition activity with the guide RNA, provides powerful genetic manipulation tools. In particular, combining the dCas9 protein and transcriptional activator to achieve specific regulation of gene expression has made important contributions to biotechnology in medical research as well as agriculture. CRISPR/dCas9 activation (CRISPRa) can increase the transcription of endogenous genes. Overexpression of foreign genes by traditional transgenic technology in plant cells is the routine method to verify gene function by elevating genes transcription. One of the main limitations of the overexpression is the vector capacity constraint that makes it difficult to express multiple genes using the typical Ti plasmid vectors from Agrobacterium. The CRISPRa system can overcome these limitations of the traditional gene overexpression method and achieve multiple gene activation by simply designating several guide RNAs in one vector. This review summarizes the latest progress based on the development of CRISPRa systems, including SunTag, dCas9-VPR, dCas9-TV, scRNA, SAM, and CRISPR-Act and their applications in plants. Furthermore, limitations, challenges of current CRISPRa systems and future prospective applications are also discussed.