Single-cell transcriptome atlas reveals developmental trajectories and a novel metabolic pathway of catechin esters in tea leaves

Tea Aphid-Induced β-Elemene Biosynthesis by CsELE Enhances JA-Dependent Herbivore Resistance in Tea Plants

Herbivore-induced plant volatiles play a crucial role in enhancing plant resistance and adaptability. However, the key volatiles triggering this response and their activation mechanisms remain largely unknown. Although tea aphids (Toxoptera aurantii) are common tea plantation pests, the underlying biochemical pathways by which aphid-induced volatiles in Camellia sinensis contribute to aphid resistance are not well understood. In this study, we measured jasmonic acid (JA) accumulation induced by aphid feeding using LC-MS and DESI-IMS analysis. We also identified a diverse array of aphid-induced volatiles through GC-MS. Notably, tea aphids stimulated β-Elemene accumulation, enhancing 'Shuchazao' resistance to tea aphids. We also evaluated the insecticidal activity of β-Elemene through contact toxicity and repellent bioassays. For the first time, we identified aphid-induced CsELE as a terpene synthase enzyme catalyzing the final step in β-Elemene biosynthesis from farnesyl pyrophosphate. Transient overexpression and gene silencing of CsELE affected β-Elemene accumulation, thereby influencing tea plant resistance to aphids. Furthermore, Dual-Luciferase and Y1H assays revealed that MYC2a, a core transcription factor in the JA signaling pathway, positively regulates CsELE expression. Overall, this study advances our understanding of the transcriptional regulation of tea volatiles in response to pest attacks, providing new insights into the ecological significance of terpene volatiles under biotic stress.

Single-Cell Transcriptome Reveals the Cellular Response to PEG-Induced Stress in Wheat Leaves

Drought is a major factor limiting the production and yield of wheat bread (Triticum aestivum). Therefore, investigating the wheat drought-related response mechanism is an urgent priority. Here, the single-cell transcriptome of drought-nonsusceptible and susceptible wheat seedlings subjected to PEG-induced stress was systematically analyzed to study the drought-related response at the cellular level. We identified five major cell types using known marker genes and constructed a wheat leaf cell atlas. On this foundation, several potential specific marker genes for each cell type were identified, which provide a reference for further cell type annotation. Moreover, we identified cellular heterogeneity in the drought-related response mechanisms and regulatory networks among cell types. Specifically, the drought response of mesophyll cells was correlated with the photosynthetic pathway. Pseudotime trajectory analysis revealed the transition of epidermal cells from their normal function to a defense response under stress. Moreover, we also characterized the genes associated with the drought response. Notably, we identified two transcription factors (TraesCS1D02G223600 and TraesCS1D02G072900) as master regulators in most cell types. Our study provides detailed insights into the response heterogeneity of cells under PEG-induced stress. The gene resources obtained in our study could be applied to breed better crop plants with improved drought tolerance.

Identification of novel drug targets and small molecule discovery for MRSA infections

Introduction

The topmost deadliest microorganism, namely, methicillin-resistant Staphylococcus aureus (MRSA), causes dreadful infections like bacteremia, pneumonia, endocarditis, and systemic inflammations. The virulence factors associated with MRSA exhibit multidrug-resistant characteristics, complicating treatment choices. So, the primary objective of this study is to identify the MRSA virulence factors and inhibiting its activity utilizing bioinformatic approaches.

Methods

The screening of novel therapeutic MRSA targets was conducted based on the predictions retrieved from non-homologous, physicochemical analysis, subcellular localization, druggability, and virulence factor examinations. Following that, flavonoid compounds were docked against specific MRSA targets using AutoDock Vina. Further, molecular dynamic simulations and binding free energy calculations were performed using simulation software.

Results

After examining 2,640 virulence factors that presented in MRSA, the heme response regulator R (HssR) was found to be a novel protein that greatly controls the levels of heme in MRSA infections. Subsequently, the binding energy calculations for flavonoid compounds and HssR revealed that the catechin provided -7.9 kcal/mol, which surpassed the standard drug, namely, vancomycin (-5.9 kcal/mol). Further, the results were validated by evaluating molecular dynamic simulation parameters like RMSD, RMSF, ROG, SASA, and PCA. Through analyzing these parameters, catechin provided a more stable, compact nature and less solvent exposure with HssR than vancomycin. Moreover, the predicted binding free energy for HssR-catechin was found to be -23.0 kcal/mol, which was less compared to the HssR-vancomycin (-16.91 kcal/mol) complex. The results suggested that the catechin was able to modulate the activity of the HssR protein effectively.

Conclusion

These potential findings revealed that heme response regulator R as a promising therapeutic target while the flavonoid compound catechin could act as alternative therapeutic inhibitor that target MRSA infections.

Not just flavor: Insights into the metabolism of tea plants

Tea, one of the world's most popular beverages, boasts a rich cultural history and distinctive flavor profiles. With advances in genomics and plant metabolism research, significant progress has been made in understanding the biosynthetic pathways and the underlying regulatory mechanisms of tea plants (Camellia sinensis). Tea metabolites play a pivotal role in determining tea flavor, and functional properties, while also being closely tied to the plant's stress resistance, environmental adaptability, and other newly discovered biological functions. In recent years, research has expanded beyond the well-characterized metabolites, such as catechins, l-theanine, and caffeine, to include volatile compounds, hormones, photosynthetic pigments, lignin, and other recently discovered metabolites, shedding new light on the intricate tea plant metabolism. This review highlights the biosynthetic pathways and regulatory mechanisms of key metabolites in tea plants, with a focus on the critical enzyme genes and regulatory factors. Additionally, emerging technologies and methodologies applied in tea plant metabolism research are briefly introduced. By further exploring the biological functions of tea metabolites and their upstream regulatory networks, future studies may offer theoretical insights and technological support for tea plant cultivation, tea quality improvement, and the sustainable development of the tea industry.

A novel transcription factor CsSNACA2 plays a pivotal role within nitrogen assimilation in tea plants

Tea (Camellia sinensis) is a globally renowned economic crop, with organs such as leaves and buds utilized for consumption. As a perennial foliage crop, tea plants have high-nitrogen consumption and demand but exhibit relatively low nitrogen use efficiency. Exploring the genetic factors involved in nitrogen assimilation in tea plants could lead to improvements in both tea yield and quality. Here, we first conducted transcriptome sequencing on two tissues (roots and young leaves) under two different nitrate levels (0.2 and 2.5 mm KNO3) and at six time points (0, 15, and 45 min; 2 and 6 h and 2 days). Differential gene expression patterns were observed for several genes that exhibited altered expression at 2 h. Clustering and enrichment analyses, along with co-expression network construction, provided evidence for the crucial involvement of CsSNACA2 in nitrogen assimilation. CsSNACA2 overexpression elicited pronounced phenotypic changes in nitrogen-deficient plants. Furthermore, CsSNACA2 suppressed the expression of CsNR (encoding nitrate reductase) and CsCLCa (encoding aNO 3 - /H+ exchanger). Moreover, CsSNACA2 served as a downstream target of CsSPL6.1. In addition, we characterized Csi-miR156e and Csi-miR156k, which directly cleave CsSPL6.1. This study identified a transcription factor module participating in nitrogen assimilation in tea plants, providing a genetic foundation for future innovations in tea cultivar improvement. These results broaden our understanding of the genetic mechanisms governing nitrogen assimilation in dicotyledonous plants.

Single-cell atlases reveal leaf cell-type-specific regulation of metal transporters in the hyperaccumulator Sedum alfredii under cadmium stress

Hyperaccumulation in plants is a complex and dynamic biological process. Sedum alfredii, the most studied Cd hyperaccumulator, can accumulate up to 9000 mg kg-1 Cd in its leaves without suffering toxicity. Although several studies have reported the molecular mechanisms of Cd hyperaccumulation, our understanding of the cell-type-specific transcriptional regulation induced by Cd remains limited. In this study, the first full-length transcriptome of S. alfredii was generated using the PacBio Iso-Seq technology. A total of 18,718,513 subreads (39.90 Gb) were obtained, with an average length of 2133 bp. The single-cell RNA sequencing was employed on leaves of S. alfredii grown under Cd stress. A total of 12,616 high-quality single cells were derived from the control and Cd-treatment samples of S. alfredii leaves. Based on cell heterogeneity and the expression profiles of previously reported marker genes, seven cell types with 12 transcriptionally distinct cell clusters were identified, thereby constructing the first single-cell atlas for S. alfredii leaves. Metal transporters such as CAX5, COPT5, ZIP5, YSL7, and MTP1 were up-regulated in different cell types of S. alfredii leaves under Cd stress. The distinctive gene expression patterns of metal transporters indicate special gene regulatory networks underlying Cd tolerance and hyperaccumulation in S. alfredii. Collectively, our findings are the first observation of the cellular and molecular responses of S. alfredii leaves under Cd stress and lay the cornerstone for future hyperaccumulator scRNA-seq investigations.

A Systematic Review: Quercetin-Secondary Metabolite of the Flavonol Class, with Multiple Health Benefits and Low Bioavailability

The main goal of this systematic review on the flavonol class secondary metabolite quercetin is to evaluate and summarize the existing research on quercetin's potential health benefits, therapeutic properties, and effectiveness in disease prevention and treatment. In addition to evaluating quercetin's potential for drug development with fewer side effects and lower toxicity, this type of review attempts to collect scientific evidence addressing quercetin's roles as an antioxidant, anti-inflammatory, antibacterial, and anticancer agent. In the first part, we analyze various flavonoid compounds, focusing on their chemical structure, classification, and natural sources. We highlight their most recent biological activities as reported in the literature. Among these compounds, we pay special attention to quercetin, detailing its chemical structure, physicochemical properties, and process of biosynthesis in plants. We also present natural sources of quercetin and emphasize its health benefits, such as its antioxidant and anti-inflammatory effects. Additionally, we discuss methods to enhance its bioavailability, analyzing the latest and most effective delivery systems based on quercetin.

The Importance of Argan Oil in Medicine and Cosmetology

Argan oil, rich in unsaturated fatty acids and polyphenols, exerts beneficial effects on both the intestinal and skin microbiotas. In the gut, it promotes the growth of beneficial bacteria, such as lactobacilli, while reducing pathogenic bacteria, due to its anti-inflammatory properties that help maintain microbial balance. Additionally, it improves the integrity of the intestinal mucosa, reducing the risk of dysbiosis. On the skin, argan oil hydrates and balances the lipid environment, creating a favorable setting for beneficial microorganisms, while also possessing antimicrobial and anti-inflammatory properties that soothe conditions like eczema and acne. Thus, argan oil is valuable for overall health, supporting digestion and skin health. The objective of this review is to provide a summary of the benefits of argan oil for alternative and complementary medicine. An exhaustive search of the literature was carried out using targeted keywords. A set of 83 articles were selected and analyzed. As the mechanisms of action of argan oil are not completely understood, this work highlighted the benefits of this oil by analyzing its nutritional properties and its beneficial effects on the intestinal and skin microbiotas. Indeed, argan oil is valuable for overall health.

Root-specific theanine metabolism and regulation at the single-cell level in tea plants (Camellia sinensis)

Root-synthesized secondary metabolites are critical quality-conferring compounds of foods, plant-derived medicines, and beverages. However, information at a single-cell level on root-specific secondary metabolism remains largely unexplored. L-Theanine, an important quality component of tea, is primarily synthesized in roots, from which it is then transported to new shoots of tea plant. In this study, we present a single-cell RNA sequencing (scRNA-seq)-derived map for the tea plant root, which enabled cell-type-specific analysis of glutamate and ethylamine (two precursors of theanine biosynthesis) metabolism, and theanine biosynthesis, storage, and transport. Our findings support a model in which the theanine biosynthesis pathway occurs via multicellular compartmentation and does not require high co-expression levels of transcription factors and their target genes within the same cell cluster. This study provides novel insights into theanine metabolism and regulation, at the single-cell level, and offers an example for studying root-specific secondary metabolism in other plant systems.

The multicellular compartmentation of plant specialized metabolism

The phenomenon of multicellular compartmentation in biosynthetic pathways has been documented for only a limited subset of specialized metabolites, despite its hypothesized significance in facilitating plant survival and adaptation to environmental stress. Transporters that shuttle metabolic intermediates between cells are hypothesized to be integral components enabling compartmentalized biosynthesis. Nevertheless, our understanding of the multicellular compartmentation of plant specialized metabolism and the associated intermediate transporters remains incomplete. The emergence of single-cell and spatial multiomics techniques holds promise for shedding light on unresolved questions in this field, such as the prevalence of multicellular compartmentation across the plant kingdom and the specific types of specialized metabolites whose biosynthetic pathways are prone to compartmentation. Advancing our understanding of the mechanisms underlying multicellular compartmentation will contribute to improving the production of specialized target metabolites through metabolic engineering or synthetic biology.

Single-cell RNA sequencing reveals a high-resolution cell atlas of petals in Prunus mume at different flowering development stages

Prunus mume (mei), a traditional ornamental plant in China, is renowned for its fragrant flowers, primarily emitted by its petals. However, the cell types of mei petals and where floral volatile synthesis occurs are rarely reported. The study used single-cell RNA sequencing to characterize the gene expression landscape in petals of P. mume 'Fenhong Zhusha' at budding stage (BS) and full-blooming stage (FS). Six major cell types of petals were identified: epidermal cells (ECs), parenchyma cells (PCs), xylem parenchyma cells, phloem parenchyma cells, xylem vessels and fibers, and sieve elements and companion cells complex. Cell-specific marker genes in each cell type were provided. Floral volatiles from mei petals were measured at four flowering development stages, and their emissions increased from BS to FS, and decreased at the withering stage. Fifty-eight differentially expressed genes (DEGs) in benzenoid/phenylpropanoid pathway were screened using bulk RNA-seq data. Twenty-eight DEGs expression increased from BS to FS, indicating that they might play roles in floral volatile synthesis in P. mume, among which PmBAHD3 would participate in benzyl acetate synthesis. ScRNA-seq data showed that 27 DEGs mentioned above were expressed variously in different cell types. In situ hybridization confirmed that PmPAL2, PmCAD1, PmBAHD3,5, and PmEGS1 involved in floral volatile synthesis in mei petals are mainly expressed in EC, PC, and most vascular tissues, consistent with scRNA-seq data. The result indicates that benzyl acetate and eugenol, the characteristic volatiles in mei, are mostly synthesized in these cell types. The first petal single-cell atlas was constructed, offering new insights into the molecular mechanism of floral volatile synthesis.

ScRNA-seq reveals the spatiotemporal distribution of camptothecin pathway and transposon activity in Camptotheca acuminata shoot apexes and leaves

Camptotheca acuminata Decne., a significant natural source of the anticancer drug camptothecin (CPT), synthesizes CPT through the monoterpene indole alkaloid (MIA) pathway. In this study, we used single-cell RNA sequencing (scRNA-seq) to generate datasets encompassing over 60,000 cells from C. acuminata shoot apexes and leaves. After cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Analysis of MIA pathway gene expression revealed that most of them exhibited heightened expression in proliferating cells (PCs) and vascular cells (VCs). In contrast to MIA biosynthesis in Catharanthus roseus, CPT biosynthesis in C. acuminata did not exhibit multicellular compartmentalization. Some putative genes encoding enzymes and transcription factors (TFs) related to the biosynthesis of CPT and its derivatives were identified through co-expression analysis. These include 19 cytochrome P450 genes, 8 O-methyltransferase (OMT) genes, and 62 TFs. Additionally, these pathway genes exhibited dynamic expression patterns during VC and EC development. Furthermore, by integrating gene and transposable element (TE) expression data, we constructed novel single-cell transcriptome atlases for C. acuminata. This approach significantly facilitated the identification of rare cell types, including peripheral zone cells (PZs). Some TE families displayed cell type specific, tissue specific, or developmental stage-specific expression patterns, suggesting crucial roles for these TEs in cell differentiation and development. Overall, this study not only provides novel insights into CPT biosynthesis and spatial-temporal TE expression characteristics in C. acuminata, but also serves as a valuable resource for further comprehensive investigations into the development and physiology of this species.

Single-cell RNA-sequencing provides new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass leaf blades

As an important warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) flourishes in warm areas around the world due to the existence of the C4 photosynthetic pathway. However, how C4 photosynthesis operates in bermudagrass leaves is still poorly understood. In this study, we performed single-cell RNA-sequencing on 5296 cells from bermudagrass leaf blades. Eight cell clusters corresponding to mesophyll, bundle sheath, epidermis and vascular bundle cells were successfully identified using known cell marker genes. Expression profiling indicated that genes encoding NADP-dependent malic enzymes (NADP-MEs) were highly expressed in bundle sheath cells, whereas NAD-ME genes were weakly expressed in all cell types, suggesting C4 photosynthesis of bermudagrass leaf blades might be NADP-ME type rather than NAD-ME type. The results also indicated that starch synthesis-related genes showed preferential expression in bundle sheath cells, whereas starch degradation-related genes were highly expressed in mesophyll cells, which agrees with the observed accumulation of starch-filled chloroplasts in bundle sheath cells. Gene co-expression analysis further revealed that different families of transcription factors were co-expressed with multiple C4 photosynthesis-related genes, suggesting a complex transcription regulatory network of C4 photosynthesis might exist in bermudagrass leaf blades. These findings collectively provided new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass.

ScRNA-seq reveals dark- and light-induced differentially expressed gene atlases of seedling leaves in Arachis hypogaea L

Although the regulatory mechanisms of dark and light-induced plant morphogenesis have been broadly investigated, the biological process in peanuts has not been systematically explored on single-cell resolution. Herein, 10 cell clusters were characterized using scRNA-seq-identified marker genes, based on 13 409 and 11 296 single cells from 1-week-old peanut seedling leaves grown under dark and light conditions. 6104 genes and 50 transcription factors (TFs) displayed significant expression patterns in distinct cell clusters, which provided gene resources for profiling dark/light-induced candidate genes. Further pseudo-time trajectory and cell cycle evidence supported that dark repressed the cell division and perturbed normal cell cycle, especially the PORA abundances correlated with 11 TFs highly enriched in mesophyll to restrict the chlorophyllide synthesis. Additionally, light repressed the epidermis cell developmental trajectory extending by inhibiting the growth hormone pathway, and 21 TFs probably contributed to the different genes transcriptional dynamic. Eventually, peanut AHL17 was identified from the profile of differentially expressed TFs, which encoded protein located in the nucleus promoted leaf epidermal cell enlargement when ectopically overexpressed in Arabidopsis through the regulatory phytohormone pathway. Overall, our study presents the different gene atlases in peanut etiolated and green seedlings, providing novel biological insights to elucidate light-induced leaf cell development at the single-cell level.

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.

Single-cell transcriptome profiling reveals the spatiotemporal distribution of triterpenoid saponin biosynthesis and transposable element activity in Gynostemma pentaphyllum shoot apexes and leaves

Gynostemma pentaphyllum (Thunb.) Makino is an important producer of dammarene-type triterpenoid saponins. These saponins (gypenosides) exhibit diverse pharmacological benefits such as anticancer, antidiabetic, and immunomodulatory effects, and have major potential in the pharmaceutical and health care industries. Here, we employed single-cell RNA sequencing (scRNA-seq) to profile the transcriptomes of more than 50,000 cells derived from G. pentaphyllum shoot apexes and leaves. Following cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Each cell type displayed substantial transcriptomic heterogeneity both within and between tissues. Examining gene expression patterns across various cell types revealed that gypenoside biosynthesis predominantly occurred in mesophyll cells, with heightened activity observed in shoot apexes compared to leaves. Furthermore, we explored the impact of transposable elements (TEs) on G. pentaphyllum transcriptomic landscapes. Our findings the highlighted the unbalanced expression of certain TE families across different cell types in shoot apexes and leaves, marking the first investigation of TE expression at the single-cell level in plants. Additionally, we observed dynamic expression of genes involved in gypenoside biosynthesis and specific TE families during epidermal and vascular cell development. The involvement of TE expression in regulating cell differentiation and gypenoside biosynthesis warrant further exploration. Overall, this study not only provides new insights into the spatiotemporal organization of gypenoside biosynthesis and TE activity in G. pentaphyllum shoot apexes and leaves but also offers valuable cellular and genetic resources for a deeper understanding of developmental and physiological processes at single-cell resolution in this species.

Single-cell RNA-seq reveals a link of ovule abortion and sugar transport in Camellia oleifera

Camellia oleifera is the most important woody oil crop in China. Seed number per fruit is an important yield trait in C. oleifera. Ovule abortion is generally observed in C. oleifera and significantly decreases the seed number per fruit. However, the mechanisms of ovule abortion remain poorly understood at present. Single-cell RNA sequencing (scRNA-seq) was performed using mature ovaries of two C. oleifera varieties with different ovule abortion rates (OARs). In total, 20,526 high-quality cells were obtained, and 18 putative cell clusters were identified. Six cell types including female gametophyte, protoxylem, protophloem, procambium, epidermis, and parenchyma cells were identified from three main tissue types of ovule, placenta, and pericarp inner layer. A comparative analysis on scRNA-seq data between high- and low-OAR varieties demonstrated that the overall expression of CoSWEET and CoCWINV in procambium cells, and CoSTP in the integument was significantly upregulated in the low-OAR variety. Both the infertile ovule before pollination and the abortion ovule producing after compatible pollination might be attributed to selective abortion caused by low sugar levels in the apoplast around procambium cells and a low capability of hexose uptake in the integument. Here, the first single-cell transcriptional landscape is reported in woody crop ovaries. Our investigation demonstrates that ovule abortion may be related to sugar transport in placenta and ovules and sheds light on further deciphering the mechanism of regulating sugar transport and the improvement of seed yield in C. oleifera.

Single-cell transcriptomic analysis of pea shoot development and cell-type-specific responses to boron deficiency

Understanding how nutrient stress impacts plant growth is fundamentally important to the development of approaches to improve crop production under nutrient limitation. Here we applied single-cell RNA sequencing to shoot apices of Pisum sativum grown under boron (B) deficiency. We identified up to 15 cell clusters based on the clustering of gene expression profiles and verified cell identity with cell-type-specific marker gene expression. Different cell types responded differently to B deficiency. Specifically, the expression of photosynthetic genes in mesophyll cells (MCs) was down-regulated by B deficiency, consistent with impaired photosynthetic rate. Furthermore, the down-regulation of stomatal development genes in guard cells, including homologs of MUTE and TOO MANY MOUTHS, correlated with a decrease in stomatal density under B deficiency. We also constructed the developmental trajectory of the shoot apical meristem (SAM) cells and a transcription factor interaction network. The developmental progression of SAM to MC was characterized by up-regulation of genes encoding histones and chromatin assembly and remodeling proteins including homologs of FASCIATA1 (FAS1) and SWITCH DEFECTIVE/SUCROSE NON-FERMENTABLE (SWI/SNF) complex. However, B deficiency suppressed their expression, which helps to explain impaired SAM development under B deficiency. These results represent a major advance over bulk-tissue RNA-seq analysis in which cell-type-specific responses are lost and hence important physiological responses to B deficiency are missed. The reported findings reveal strategies by which plants adapt to B deficiency thus offering breeders a set of specific targets for genetic improvement. The reported approach and resources have potential applications well beyond P. sativum species and could be applied to various legumes to improve their adaptability to multiple nutrient or abiotic stresses.

A Single-Nucleus Resolution Atlas of Transcriptome and Chromatin Accessibility for Peanut (Arachis Hypogaea L.) Leaves

The peanut is an important worldwide cash-crop for edible oil and protein. However, the kinetic mechanisms that determine gene expression and chromatin accessibility during leaf development in peanut represented allotetraploid leguminous crops are poorly understood at single-cell resolution. Here, a single-nucleus atlas of peanut leaves is developed by simultaneously profiling the transcriptome and chromatin accessibility in the same individual-cell using fluorescence-activated sorted single-nuclei. In total, 5930 cells with 50 890 expressed genes are classified into 18 cell-clusters, and 5315 chromatin fragments are enriched with 26 083 target genes in the chromatin accessible landscape. The developmental trajectory analysis reveals the involvement of the ethylene-AP2 module in leaf cell differentiation, and cell-cycle analysis demonstrated that genome replication featured in distinct cell-types with circadian rhythms transcription factors (TFs). Furthermore, dual-omics illustrates that the fatty acid pathway modulates epidermal-guard cells differentiation and providescritical TFs interaction networks for understanding mesophyll development, and the cytokinin module (LHY/LOG) that regulates vascular growth. Additionally, an AT-hook protein AhAHL11 is identified that promotes leaf area expansion by modulating the auxin content increase. In summary, the simultaneous profiling of transcription and chromatin accessibility landscapes using snRNA/ATAC-seq provides novel biological insights into the dynamic processes of peanut leaf cell development at the cellular level.

A highly efficient protocol for isolation of protoplast from China, Assam and Cambod types of tea plants [Camellia sinensis (L.) O. Kuntze]

BACKGROUND

Tea is the most popular beverage worldwide second only to water. Its demand is tremendously rising due to increased awareness of its medicinal importance. The quality and uses of tea depend on the tea-types which are mainly three types including China, Assam and Cambod type having distinct compositions of secondary metabolites. Huge variation in secondary metabolites in different tea-types and cultivars limited the successful application of various approaches used for its trait improvement. The efficiency of a protocol for isolation of protoplast is specific to the types and cultivars of tea plants. The existing tea protoplast-isolation protocols [which were optimized for tea-types (China and Assam type) and Chinese cultivars grown in China] were found ineffective on types/cultivars grown in India due to type/cultivar variability. Therefore, optimization of protoplast-isolation protocol is essential for tea-types/cultivars grown in India, as it is the second largest producer of tea and the largest producer of black tea. Here, efforts were made to develop an efficient protoplast-isolation protocol from all major types of tea (China, Assam and Cambod types) grown in India and also from three types of tender leaves obtained from field-grown, hydroponically-grown and tissue culture-grown tea plants.

RESULTS

Developed protoplast-isolation protocol was effective for different types of leaf tissue obtained from the tender leaves of field-grown, hydroponically-grown and tissue culture-grown tea plants. Moreover, optimized protocol effectively worked on all three types of tea including China, Assam and Cambod types cultivated in India. The digestion of leaves with 3% cellulase R-10, 0.6% macerozyme, 1% hemicellulase and 4% polyvinylpyrrolidone for 12 h at 28ºC yielded approximately 3.8-4.6 × 107 protoplasts per gram fresh tissue and 80-95% viability in selected tea cultivars, and tissue culture plant material was found most appropriate for protoplast isolation.

CONCLUSIONS

In conclusion, we reported an efficient protocol for isolation of protoplasts from tender tea leaves of all major tea-types (China, Assam and Cambod) grown in India. Moreover, the protocol is also effective for tender-leaf tissue of field-grown, hydroponically-grown and tissue culture-grown tea plants. The findings are expected to contribute to the genetic improvement of tea traits widely.

Genome-wide identification and expression analysis of heat shock protein gene family in cassava

Heat shock proteins are important molecular chaperones that are involved in plant growth and stress responses. However, members of the Hsp family have been poorly studied in cassava. In this study, 225 MeHsp genes were identified in the cassava genome, and their genetic structures exhibited relatively conserved features within each subfamily. The 225 MeHsp genes showed random chromosomal distribution, and at least 74 pairs of segmentally duplicated MeHsp genes. Eleven tandemly duplicated MeHsp genes were identified. Cis-element analysis revealed the importance of MeHsps in plant adaptations to the environment. The prediction of protein interactions suggested that MeHsp70-20 may play a critical regulatory role in the interactive network. Furthermore, the expression profiles of MeHsps in different tissues and cell subsets were analyzed using bulk transcriptomics and single-cell transcriptomic data. Several subfamily genes exhibited unique expression patterns in the transcriptome and were selected for detailed analysis of the single-cell transcriptome. Quantitative real-time polymerase chain reaction (qRT-PCR) revealed the expression patterns of these genes under temperature stress, further supporting the prediction of cis-acting elements. This study provides valuable information for understanding the functional characteristics of MeHsp genes and the evolutionary relationships between MeHsps.

Application and prospects of single-cell and spatial omics technologies in woody plants

Over the past decade, high-throughput sequencing and high-resolution single-cell transcriptome sequencing technologies have undergone rapid development, leading to significant breakthroughs. Traditional molecular biology methods are limited in their ability to unravel cellular-level heterogeneity within woody plant tissues. Consequently, techniques such as single-cell transcriptomics, single-cell epigenetics, and spatial transcriptomics are rapidly gaining popularity in the study of woody plants. In this review, we provide a comprehensive overview of the development of these technologies, with a focus on their applications and the challenges they present in single-cell transcriptome research in woody plants. In particular, we delve into the similarities and differences among the results of current studies and analyze the reasons behind these differences. Furthermore, we put forth potential solutions to overcome the challenges encountered in single-cell transcriptome applications in woody plants. Finally, we discuss the application directions of these techniques to address key challenges in woody plant research in the future.

Genome-Wide Association Studies of Biluochun Tea Plant Populations in Dongting Mountain and Comprehensive Identification of Candidate Genes Associated with Core Agronomic Traits by Four Analysis Models

The elite germplasm resources are key to the beautiful appearance and pleasant flavor of Biluochun tea. We collected and measured the agronomic traits of 95 tea plants to reveal the trait diversity and breeding value of Biluochun tea plant populations. The results revealed that the agronomic traits of Biluochun tea plant populations were diverse and had high breeding value. Additionally, we resequenced these tea plant populations to reveal genetic diversity, population structure, and selection pressure. The Biluochun tea plant populations contained two groups and were least affected by natural selection based on the results of population structure and selection pressure. More importantly, four non-synonymous single nucleotide polymorphisms (nsSNPs) and candidate genes associated with (-)-gallocatechin gallate (GCG), (-)-gallocatechin (GC), and caffeine (CAF) were detected using at least two GWAS models. The results will promote the development and application of molecular markers and the utilization of elite germplasm from Biluochun populations.

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.

scRNA-seq Reveals the Mechanism of Fatty Acid Desaturase 2 Mutation to Repress Leaf Growth in Peanut (Arachis hypogaea L.)

Fatty Acid Desaturase 2 (FAD2) controls the conversion of oleic acids into linoleic acids. Mutations in FAD2 not only increase the high-oleic content, but also repress the leaf growth. However, the mechanism by which FAD2 regulates the growth pathway has not been elucidated in peanut leaves with single-cell resolution. In this study, we isolated fad2 mutant leaf protoplast cells to perform single-cell RNA sequencing. Approximately 24,988 individual cells with 10,249 expressed genes were classified into five major cell types. A comparative analysis of 3495 differentially expressed genes (DEGs) in distinct cell types demonstrated that fad2 inhibited the expression of the cytokinin synthesis gene LOG in vascular cells, thereby repressing leaf growth. Further, pseudo-time trajectory analysis indicated that fad2 repressed leaf cell differentiation, and cell-cycle evidence displayed that fad2 perturbed the normal cell cycle to induce the majority of cells to drop into the S phase. Additionally, important transcription factors were filtered from the DEG profiles that connected the network involved in high-oleic acid accumulation (WRKY6), activated the hormone pathway (WRKY23, ERF109), and potentially regulated leaf growth (ERF6, MYB102, WRKY30). Collectively, our study describes different gene atlases in high-oleic and normal peanut seedling leaves, providing novel biological insights to elucidate the molecular mechanism of the high-oleic peanut-associated agronomic trait at the single-cell level.

Mass spectrometry imaging and single-cell transcriptional profiling reveal the tissue-specific regulation of bioactive ingredient biosynthesis in Taxus leaves

Taxus leaves provide the raw industrial materials for taxol, a natural antineoplastic drug widely used in the treatment of various cancers. However, the precise distribution, biosynthesis, and transcriptional regulation of taxoids and other active components in Taxus leaves remain unknown. Matrix-assisted laser desorption/ionization-mass spectrometry imaging analysis was used to visualize various secondary metabolites in leaf sections of Taxus mairei, confirming the tissue-specific accumulation of different active metabolites. Single-cell sequencing was used to produce expression profiles of 8846 cells, with a median of 2352 genes per cell. Based on a series of cluster-specific markers, cells were grouped into 15 clusters, suggesting a high degree of cell heterogeneity in T. mairei leaves. Our data were used to create the first Taxus leaf metabolic single-cell atlas and to reveal spatial and temporal expression patterns of several secondary metabolic pathways. According to the cell-type annotation, most taxol biosynthesis genes are expressed mainly in leaf mesophyll cells; phenolic acid and flavonoid biosynthesis genes are highly expressed in leaf epidermal cells (including the stomatal complex and guard cells); and terpenoid and steroid biosynthesis genes are expressed specifically in leaf mesophyll cells. A number of novel and cell-specific transcription factors involved in secondary metabolite biosynthesis were identified, including MYB17, WRKY12, WRKY31, ERF13, GT_2, and bHLH46. Our research establishes the transcriptional landscape of major cell types in T. mairei leaves at a single-cell resolution and provides valuable resources for studying the basic principles of cell-type-specific regulation of secondary metabolism.

Single-cell RNA sequencing analysis of the embryogenic callus clarifies the spatiotemporal developmental trajectories of the early somatic embryo in Dimocarpus longan

Plant embryogenic calli (ECs) can undergo somatic embryogenesis to regenerate plants. This process is mediated by regulatory factors, such as transcription factors and specifically expressed genes, but the precise molecular mechanisms underlying somatic embryogenesis at the single-cell level remain unclear. In this study, we performed high-resolution single-cell RNA sequencing analysis to determine the cellular changes in the EC of the woody plant species Dimocarpus longan (longan) and clarify the continuous cell differentiation trajectories at the transcriptome level. The highly heterogeneous cells in the EC were divided into 12 putative clusters (e.g., proliferating, meristematic, vascular, and epidermal cell clusters). We determined cluster-enriched expression marker genes and found that overexpression of the epidermal cell marker gene GDSL ESTERASE/LIPASE-1 inhibited the hydrolysis of triacylglycerol. In addition, the stability of autophagy was critical for the somatic embryogenesis of longan. The pseudo-timeline analysis elucidated the continuous cell differentiation trajectories from early embryonic cell division to vascular and epidermal cell differentiation during the somatic embryogenesis of longan. Moreover, key transcriptional regulators associated with cell fates were revealed. We found that ETHYLENE RESPONSIVE FACTOR 6 was characterized as a heat-sensitive factor that negatively regulates longan somatic embryogenesis under high-temperature stress conditions. The results of this study provide new spatiotemporal insights into cell division and differentiation during longan somatic embryogenesis at single-cell resolution.

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.

Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants

Single-cell and spatial transcriptomics have diverted researchers' attention from the multicellular level to the single-cell level and spatial information. Single-cell transcriptomes provide insights into the transcriptome at the single-cell level, whereas spatial transcriptomes help preserve spatial information. Although these two omics technologies are helpful and mature, further research is needed to ensure their widespread applicability in plant studies. Reviewing recent research on plant single-cell or spatial transcriptomics, we compared the different experimental methods used in various plants. The limitations and challenges are clear for both single-cell and spatial transcriptomic analyses, such as the lack of applicability, spatial information, or high resolution. Subsequently, we put forth further applications, such as cross-species analysis of roots at the single-cell level and the idea that single-cell transcriptome analysis needs to be combined with other omics analyses to achieve superiority over individual omics analyses. Overall, the results of this review suggest that combining single-cell transcriptomics, spatial transcriptomics, and spatial element distribution can provide a promising research direction, particularly for plant research.

Application of Multi-Perspectives in Tea Breeding and the Main Directions

Tea plants are an economically important crop and conducting research on tea breeding contributes to enhancing the yield and quality of tea leaves as well as breeding traits that satisfy the requirements of the public. This study reviews the current status of tea plants germplasm resources and their utilization, which has provided genetic material for the application of multi-omics, including genomics and transcriptomics in breeding. Various molecular markers for breeding were designed based on multi-omics, and available approaches in the direction of high yield, quality and resistance in tea plants breeding are proposed. Additionally, future breeding of tea plants based on single-cellomics, pangenomics, plant-microbe interactions and epigenetics are proposed and provided as references. This study aims to provide inspiration and guidance for advancing the development of genetic breeding in tea plants, as well as providing implications for breeding research in other crops.

5mC DNA methylation modification-mediated regulation in tissue functional differentiation and important flavor substance synthesis of tea plant (Camellia sinensis L.)

In plants, 5mC DNA methylation is an important and conserved epistatic mark involving genomic stability, gene transcriptional regulation, developmental regulation, abiotic stress response, metabolite synthesis, etc. However, the roles of 5mC DNA methylation modification (5mC methylation) in tea plant growth and development (in pre-harvest processing) and flavor substance synthesis in pre- and post-harvest processing are unknown. We therefore conducted a comprehensive methylation analysis of four key pre-harvest tissues (root, leaf, flower, and fruit) and two processed leaves during oolong tea post-harvest processing. We found that differential 5mC methylation among four key tissues is closely related to tissue functional differentiation and that genes expressed tissue-specifically, responsible for tissue-specific functions, maintain relatively low 5mC methylation levels relative to non-tissue-specifically expressed genes. Importantly, hypomethylation modifications of CsAlaDC and TS/GS genes in roots provided the molecular basis for the dominant synthesis of theanine in roots. In addition, integration of 5mC DNA methylationomics, metabolomics, and transcriptomics of post-harvest leaves revealed that content changes in flavor metabolites during oolong tea processing were closely associated with transcription level changes in corresponding metabolite synthesis genes, and changes in transcript levels of these important synthesis genes were strictly regulated by 5mC methylation. We further report that some key genes during processing are regulated by 5mC methylation, which can effectively explain the content changes of important aroma metabolites, including α-farnesene, nerolidol, lipids, and taste substances such as catechins. Our results not only highlight the key roles of 5mC methylation in important flavor substance synthesis in pre- and post-harvest processing, but also provide epimutation-related gene targets for future improvement of tea quality or breeding of whole-tissue high-theanine varieties.

Single-cell transcriptome atlas reveals spatiotemporal developmental trajectories in the basal roots of moso bamboo (Phyllostachys edulis)

Roots are essential for plant growth and development. Bamboo is a large Poaceae perennial with 1642 species worldwide. However, little is known about the transcriptional atlas that underpins root cell-type differentiation. Here, we set up a modified protocol for protoplast preparation and report single-cell transcriptomes of 14 279 filtered single cells derived from the basal root tips of moso bamboo. We identified four cell types and defined new cell-type-specific marker genes for the basal root. We reconstructed the developmental trajectories of the root cap, epidermis, and ground tissues and elucidated critical factors regulating cell fate determination. According to in situ hybridization and pseudotime trajectory analysis, the root cap and epidermis originated from a common initial cell lineage, revealing the particularity of bamboo basal root development. We further identified key regulatory factors for the differentiation of these cells and indicated divergent root developmental pathways between moso bamboo and rice. Additionally, PheWOX13a and PheWOX13b ectopically expressed in Arabidopsis inhibited primary root and lateral root growth and regulated the growth and development of the root cap, which was different from WOX13 orthologs in Arabidopsis. Taken together, our results offer an important resource for investigating the mechanism of root cell differentiation and root system architecture in perennial woody species of Bambusoideae.

The glycoside hydrolase gene family profile and microbial function of Debaryomyces hansenii Y4 during South-road dark tea fermentation

Microbes are crucial to the quality formation of Sichuan South-road Dark Tea (SSDT) during pile-fermentation, but their mechanism of action has not yet been elucidated. Here, the glycoside hydrolase (GH) gene family and microbial function of Debaryomyces hansenii Y4 during solid-state fermentation were analyzed, and the results showed that many GH genes being distributed in comparatively abundant GH17, GH18, GH76, GH31, GH47, and GH2 were discovered in D. hansenii. They encoded beta-galactosidase, alpha-D-galactoside galactohydrolase, alpha-xylosidase, mannosidase, etc., and most of the GHs were located in the exocellular space and participated in the degradation of polysaccharides and oligosaccharides. D. hansenii Y4 could develop the mellow mouthfeel and "reddish brown" factors of SSDT via increasing the levels of water extracts, soluble sugars and amino acids but decreasing the tea polyphenols and caffeine levels, combined with altering the levels of thearubiins and brown index. It may facilitate the isomerization between epicatechin gallate and catechin gallate. Moreover, the expression levels of DEHA2G24860g (Beta-galactosidase gene) and DEHA2G08602g (Mannan endo-1,6-alpha-mannosidase DFG5 gene) were sharply up-regulated in fermentative anaphase, and they were significantly and negatively correlated with epicatechin content, especially, the expression of DEHA2G08602g was significantly and negatively correlated with catechin gallate level. It was hypothesized that D. hansenii Y4 is likely to be an important functional microbe targeting carbohydrate destruction and catechin transformation during SSDT pile-fermentation, with DEHA2G08602g as a key thermotolerant functional gene.

Single-cell transcriptomic analyses reveal cellular and molecular patterns of rubber tree response to early powdery mildew infection

As a perennial woody plant, the rubber tree (Hevea brasiliensis) must adapt to various environmental challenges through gene expression in multiple cell types. It is still unclear how genes in this species are expressed at the cellular level and the precise mechanisms by which cells respond transcriptionally to environmental stimuli, especially in the case of pathogen infection. Here, we characterized the transcriptomes in Hevea leaves during early powdery mildew infection using single-cell RNA sequencing. We identified 10 cell types and constructed the first single-cell atlas of Hevea leaves. Distinct gene expression patterns of the cell clusters were observed under powdery mildew infection, which was especially significant in the epidermal cells. Most of the genes involved in host-pathogen interactions in epidermal cells exhibited a pattern of dramatically increased expression with increasing pseudotime. Interestingly, we found that the HbCNL2 gene, encoding a nucleotide-binding leucine-rich repeat protein, positively modulated the defence of rubber leaves against powdery mildew. Overexpression of the HbCNL2 gene triggered a typical cell death phenotype in tobacco leaves and a higher level of reactive oxygen species in the protoplasts of Hevea leaves. The HbCNL2 protein was located in the cytomembrane and nucleus, and its leucine-rich repeat domain interacted with the histidine kinase-like ATPase domain of the molecular chaperone HbHSP90 in the nucleus. Collectively, our results provide the first observation of the cellular and molecular responses of Hevea leaves to biotrophic pathogen infection and can guide the identification of disease-resistance genes in this important tree species.

Genome-wide identification of UDP-glycosyltransferases in the tea plant (Camellia sinensis) and their biochemical and physiological functions

Tea (Camellia sinensis) has been an immensely important commercially grown crop for decades. This is due to the presence of essential nutrients and plant secondary metabolites that exhibit beneficial health effects. UDP-glycosyltransferases (UGTs) play an important role in the diversity of such secondary metabolites by catalysing the transfer of an activated sugar donor to acceptor molecules, and thereby creating a huge variety of glycoconjugates. Only in recent years, thanks to the sequencing of the tea plant genome, have there been increased efforts to characterise the UGTs in C. sinensis to gain an understanding of their physiological role and biotechnological potential. Based on the conserved plant secondary product glycosyltransferase (PSPG) motif and the catalytically active histidine in the active site, UGTs of family 1 in C. sinensis are identified here, and shown to cluster into 21 groups in a phylogenetic tree. Building on this, our current understanding of recently characterised C. sinensis UGTs (CsUGTs) is highlighted and a discussion on future perspectives made.

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.

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.