Transcriptomic and functional analysis provides molecular insights into multicellular trichome development

Integrated multi-omics analysis reveals genes involved in flavonoid biosynthesis and trichome development of Artemisia argyi

Artemisia argyi is an herbaceous plant of the genus Artemisia. Its young and mature leaves are used as food and medicine, respectively. Glandular trichomes (GTs) are distributed on the leaf surface in A. argyi and are generally considered the location of flavonoid biosynthesis and accumulation. However, the mechanism of flavonoid biosynthesis and accumulation in A. argyi remains unclear. In this study, the coregulatory genes involved in flavonoid biosynthesis and trichome development in this species were screened and evaluated, and the biosynthetic pathways for key flavonoids in A. argyi were uncovered. AaMYB1 and AaYABBY1 were screened using weighted gene co-expression network analysis, and both genes were then genetically transformed into Nicotiana tabacum L. cv. K326 (tobacco). Simultaneously, AaYABBY1 was also genetically transformed into Arabidopsis thaliana. The total flavonoid and rutin contents were increased in tobacco plants overexpressing AaMYB1 and AaYABBY1, and the expression levels of genes participating in the flavonoid synthesis pathway, such as PAL, FLS, and F3H, were significantly up-regulated in plants overexpressing these genes. These results indicated that AaMYB1 and AaYABBY1 promote flavonoid biosynthesis in tobacco. Furthermore, compared to that in the wild-type, the trichome density was significantly increased in tobacco and A. thaliana plants overexpressing AaYABBY1. These results confirm that AaYABBY1 might be involved in regulating trichome formation in A. argyi. This indicates the potential genes involved in and provides new insights into the development of trichome cellular factories based on the "development-metabolism" interaction network and the cultivation of high-quality A. argyi.

An HD-Zip III transcription factor, BjPHVa, negatively regulates non-glandular trichome formation in Brassica juncea

Trichomes are specialized structures derived from epidermal cells. Non-glandular trichomes primarily protect plants from herbivores and intense ultraviolet radiation by acting as a physical barrier. Recent research has highlighted the importance of homeodomain leucine zipper (HD-Zip) IV transcription factors (TFs) in promoting trichome development. In this study, an HD-Zip III TF called PHAVOLUTA (BjPHVa) was identified as a negative regulator of non-glandular trichome initiation in Brassica juncea. Genome editing of BjPHVa resulted in a significant increase in trichome number in B. juncea. Co-expression networks revealed a strong association between trichome development and the HD-Zip family, which was supported by transcriptomic analysis findings. An R2R3-MYB TF, BjGL1a, a key regulator of trichome development, was found to be associated with BjPHVa-regulated trichome development. Knockdown of BjGL1a expression resulted in reduced trichome number in B. juncea. BjPHVa was observed to interact directly with BjGL1a while binding to the BjGL1a promoter, resulting in the inhibition of BjGL1a transcription. These results provide new insights into the identification of regulators involved in trichome development and offer new opportunities to enhance resistance to predicted stresses through genome editing targeting PHVa within Brassicaceae.

Deep learning-based quantification and transcriptomic profiling reveal a methyl jasmonate-mediated glandular trichome formation pathway in Cannabis sativa

Cannabis glandular trichomes (GTs) are economically and biotechnologically important structures that have a remarkable morphology and capacity to produce, store, and secrete diverse classes of secondary metabolites. However, our understanding of the developmental changes and the underlying molecular processes involved in cannabis GT development is limited. In this study, we developed Cannabis Glandular Trichome Detection Model (CGTDM), a deep learning-based model capable of differentiating and quantifying three types of cannabis GTs with a high degree of efficiency and accuracy. By profiling at eight different time points, we captured dynamic changes in gene expression, phenotypes, and metabolic processes associated with GT development. By integrating weighted gene co-expression network analysis with CGTDM measurements, we established correlations between phenotypic variations in GT traits and the global transcriptome profiles across the developmental gradient. Notably, we identified a module containing methyl jasmonate (MeJA)-responsive genes that significantly correlated with stalked GT density and cannabinoid content during development, suggesting the existence of a MeJA-mediated GT formation pathway. Our findings were further supported by the successful promotion of GT development in cannabis through exogenous MeJA treatment. Importantly, we have identified CsMYC4 as a key transcription factor that positively regulates GT formation via MeJA signaling in cannabis. These findings provide novel tools for GT detection and counting, as well as valuable information for understanding the molecular regulatory mechanism of GT formation, which has the potential to facilitate the molecular breeding, targeted engineering, informed harvest timing, and manipulation of cannabinoid production.

ELONGATED HYPOTCOTYL5 and SPINE BASE SIZE1 together mediate light-regulated spine expansion in cucumber

Plant trichome development is influenced by diverse developmental and environmental signals, but the molecular mechanisms involved are not well understood in most plant species. Fruit spines (trichomes) are an important trait in cucumber (Cucumis sativus L.), as they affect both fruit smoothness and commercial quality. Spine Base Size1 (CsSBS1) has been identified as essential for regulating fruit spine size in cucumber. Here, we discovered that CsSBS1 controls a season-dependent phenotype of spine base size in wild-type plants. Decreased light intensity led to reduced expression of CsSBS1 and smaller spine base size in wild-type plants, but not in the mutants with CsSBS1 deletion. Additionally, knockout of CsSBS1 resulted in smaller fruit spine base size and eliminated the light-induced expansion of spines. Overexpression of CsSBS1 increased spine base size and rescued the decrease in spine base size under low light conditions. Further analysis revealed that ELONGATED HYPOTCOTYL5 (HY5), a major transcription factor involved in light signaling pathways, directly binds to the promoter of CsSBS1 and activates its expression. Knockout of CsHY5 led to smaller fruit spine base size and abolished the light-induced expansion of spines. Taken together, our study findings have clarified a CsHY5-CsSBS1 regulatory module that mediates light-regulated spine expansion in cucumber. This finding offers a strategy for cucumber breeders to develop fruit with stable appearance quality under changing light conditions.

The CsTM alters multicellular trichome morphology and enhances resistance against aphid by interacting with CsTIP1;1 in cucumber

The multicellular trichomes of cucumber (Cucumis sativus L.) serve as the primary defense barrier against external factors, whose impact extends beyond plant growth and development to include commercial characteristics of fruits. The aphid (Aphis gossypii Glover) is one of prominent pests in cucumber cultivation. However, the relationship between physical properties of trichomes and the aphid resistance at molecular level remains largely unexplored. Here, a spontaneous mutant trichome morphology (tm) was characterized by increased susceptibility towards aphid. Further observations showed the tm exhibited a higher and narrower trichome base, which was significantly distinguishable from that in wild-type (WT). We conducted map-based cloning and identified the candidate, CsTM, encoding a C-lectin receptor-like kinase. The knockout mutant demonstrated the role of CsTM in trichome morphogenesis. The presence of SNP does not regulate the relative expression of CsTM, but diminishes the CsTM abundance of membrane proteins in tm. Interestingly, CsTM was found to interact with CsTIP1;1, which encodes an aquaporin with extensive reports in plant resistance and growth development. The subsequent aphid resistance experiments revealed that both CsTM and CsTIP1;1 regulated the development of trichomes and conferred resistance against aphid by affecting cytoplasmic H2O2 contents. Transcriptome analysis revealed a significant enrichment of genes associated with pathogenesis, calcium binding and cellulose synthase. Overall, our study elucidates an unidentified mechanism that CsTM-CsTIP1;1 alters multicellular trichome morphology and enhances resistance against aphid, thus providing a wholly new perspective for trichome morphogenesis in cucumber.

Advances in understanding epigenetic regulation of plant trichome development: a comprehensive review

Plant growth and development are controlled by a complex gene regulatory network, which is currently a focal point of research. It has been established that epigenetic factors play a crucial role in plant growth. Trichomes, specialized appendages that arise from epidermal cells, are of great significance in plant growth and development. As a model system for studying plant development, trichomes possess both commercial and research value. Epigenetic regulation has only recently been implicated in the development of trichomes in a limited number of studies, and microRNA-mediated post-transcriptional regulation appears to dominate in this context. In light of this, we have conducted a review that explores the interplay between epigenetic regulations and the formation of plant trichomes, building upon existing knowledge of hormones and transcription factors in trichome development. Through this review, we aim to deepen our understanding of the regulatory mechanisms underlying trichome formation and shed light on future avenues of research in the field of epigenetics as it pertains to epidermal hair growth.

QTL mapping and stability analysis of trichome density in zucchini (Cucurbita pepo L.)

Trichomes provide an excellent model for studying cell differentiation and proliferation. The aboveground tissues of plants with long dense trichomes (LDTs) can cause skin itching in people working in a zucchini field, in which management, pollination, and fruit harvesting are difficult. In this study, an F2 population was constructed with the LDT inbred line "16" and the sparse micro trichome (SMT) inbred line "63" for QTL analysis of type I and II trichome density. Two QTLs were identified on chromosomes 3 and 15 using the QTL-seq method. Additionally, 191 InDel markers were developed on 20 chromosomes, a genetic map was constructed for QTL mapping, and three QTLs were identified on chromosomes 3, 6, and 15. Two QTLs, CpTD3.1 and CpTD15.1, were identified in both QTL-seq and genetic map-based QTL analyses, and CpTD15.1 was the major-effect QTL. The stability of CpTD3.1 and CpTD15.1 was confirmed using data from F2 plants under different environmental conditions. The major-effect QTL CpTD15.1 was located between markers chr15-4991349 and chr15-5766791, with a physical distance of 775.44 kb, and explained 12.71%-29.37% of the phenotypic variation observed in the three environments. CpTD3.1 was located between markers chr3-218350 and chr3-2891236, in a region with a physical distance of 2,672.89 kb, and explained 5.00%-10.64% of the phenotypic variation observed in the three environments. The functional annotations of the genes within the CpTD15.1 region were predicted, and five genes encoding transcription factors regulating trichome development were selected. Cp4.1LG15g04400 encoded zinc finger protein (ZFP) and harbored nonsynonymous SNPs in the conserved ring finger domain between the two parental lines. There were significant differences in Cp4.1LG15g04400 expression between "16" and "63", and a similar pattern was found between germplasm resources of LDT lines and SMT lines. It was presumed that Cp4.1LG15g04400 might regulate trichome density in zucchini. These results lay a foundation for better understanding the density of multicellular nonglandular trichomes and the regulatory mechanism of trichome density in zucchini.

Identification of candidate genes that regulate the trade-off between seedling cold tolerance and fruit quality in melon (Cucumis melo L.)

Trade-offs between survival and growth are widely observed in plants. Melon is an annual, trailing herb that produces economically valuable fruits that are traditionally cultivated in early spring in China. Melon seedlings are sensitive to low temperatures, and thus usually suffer from cold stress during the early growth period. However, little is known about the mechanism behind the trade-offs between seedling cold tolerance and fruit quality in melon. In this study, a total of 31 primary metabolites were detected from the mature fruits of eight melon lines that differ with respect to seedling cold tolerance; these included 12 amino acids, 10 organic acids, and 9 soluble sugars. Our results showed that concentrations of most of the primary metabolites in the cold-resistant melons were generally lower than in the cold-sensitive melons; the greatest difference in metabolite levels was observed between the cold-resistant line H581 and the moderately cold-resistant line HH09. The metabolite and transcriptome data for these two lines were then subjected to weighted correlation network analysis, resulting in the identification of five key candidate genes underlying the balancing between seedling cold tolerance and fruit quality. Among these genes, CmEAF7 might play multiple roles in regulating chloroplast development, photosynthesis, and the ABA pathway. Furthermore, multi-method functional analysis showed that CmEAF7 can certainly improve both seedling cold tolerance and fruit quality in melon. Our study identified an agriculturally important gene, CmEAF7, and provides a new insight into breeding methods to develop melon cultivars with seedling cold tolerance and high fruit quality.

Novel players in organogenesis and flavonoid biosynthesis in cucumber glandular trichomes

Glandular trichomes (GTs) are outgrowths of plant epidermal cells that secrete and store specialized secondary metabolites that protect plants against biotic and abiotic stresses and have economic importance for human use. While extensive work has been done to understand the molecular mechanisms of trichome organogenesis in Arabidopsis (Arabidopsis thaliana), which forms unicellular, non-glandular trichomes (NGTs), little is known about the mechanisms of GT development or regulation of secondary metabolites in plants with multicellular GTs. Here, we identified and functionally characterized genes associated with GT organogenesis and secondary metabolism in GTs of cucumber (Cucumis sativus). We developed a method for effective separation and isolation of cucumber GTs and NGTs. Transcriptomic and metabolomic analyses showed that flavonoid accumulation in cucumber GTs is positively associated with increased expression of related biosynthesis genes. We identified 67 GT development-related genes, the functions of 7 of which were validated by virus-induced gene silencing. We further validated the role of cucumber ECERIFERUM1 (CsCER1) in GT organogenesis by overexpression and RNA interference transgenic approaches. We further show that the transcription factor TINY BRANCHED HAIR (CsTBH) serves as a central regulator of flavonoid biosynthesis in cucumber glandular trichomes. Work from this study provides insight into the development of secondary metabolite biosynthesis in multi-cellular glandular trichomes.

Identification and Functional Characterization of CsMYCs in Cucumber Glandular Trichome Development

Glandular trichomes (GTs), specialized structures formed by the differentiation of plant epidermal cells, are known to play important roles in the resistance of plants to external biotic and abiotic stresses. These structures are capable of storing and secreting secondary metabolites, which often have important agricultural and medicinal values. In order to better understand the molecular developmental mechanisms of GTs, studies have been conducted in a variety of crops, including tomato (Solanum lycopersicum), sweetworm (Artemisia annua), and cotton (Gossypium hirsutum). The MYC transcription factor of the basic helix-loop-helix (bHLH) transcription factor family has been found to play an important role in GT development. In this study, a total of 13 cucumber MYC transcription factors were identified in the cucumber (Cucumis sativus L.) genome. After performing phylogenetic analyses and conserved motifs on the 13 CsMYCs in comparison to previously reported MYC transcription factors that regulate trichome development, seven candidate MYC transcription factors were selected. Through virus-induced gene silencing (VIGS), CsMYC2 is found to negatively regulate GT formation while CsMYC4, CsMYC5, CsMYC6, CsMYC7, and CsMYC8 are found to positively regulate GT formation. Furthermore, the two master effector genes, CsMYC2 and CsMYC7, are observed to have similar expression patterns indicating that they co-regulate the balance of GT development in an antagonistic way.

A HD-ZIP transcription factor specifies fates of multicellular trichomes via dosage-dependent mechanisms in tomato

Hair-like structures are shared by most living organisms. The hairs on plant surfaces, commonly referred to as trichomes, form diverse types to sense and protect against various stresses. However, it is unclear how trichomes differentiate into highly variable forms. Here, we show that a homeodomain leucine zipper (HD-ZIP) transcription factor named Woolly controls the fates of distinct trichomes in tomato via a dosage-dependent mechanism. The autocatalytic reinforcement of Woolly is counteracted by an autoregulatory negative feedback loop, creating a circuit with a high or low Woolly level. This biases the transcriptional activation of separate antagonistic cascades that lead to different trichome types. Our results identify the developmental switch of trichome formation and provide mechanistic insights into the progressive fate specification in plants, as well as a path to enhancing plant stress resistance and the production of beneficial chemicals.

Comparative Transcriptomic Analysis Revealed the Suppression and Alternative Splicing of Kiwifruit (Actinidia latifolia) NAP1 Gene Mediating Trichome Development

Kiwifruit (Actinidia chinensis) is commonly covered by fruit hairs (trichomes) that affect kiwifruit popularity in the commercial market. However, it remains largely unknown which gene mediates trichome development in kiwifruit. In this study, we analyzed two kiwifruit species, A. eriantha (Ae) with long, straight, and bushy trichomes and A. latifolia (Al) with short, distorted, and spare trichomes, by second- and third-generation RNA sequencing. Transcriptomic analysis indicated that the expression of the NAP1 gene, a positive regulator of trichome development, was suppressed in Al compared with that in Ae. Additionally, the alternative splicing of AlNAP1 produced two short transcripts (AlNAP1-AS1 and AlNAP1-AS2) lacking multiple exons, in addition to a full-length transcript of AlNAP1-FL. The defects of trichome development (short and distorted trichome) in Arabidopsis nap1 mutant were rescued by AlNAP1-FL but not by AlNAP1-AS1. AlNAP1-FL gene does not affect trichome density in nap1 mutant. The qRT-PCR analysis indicated that the alternative splicing further reduces the level of functional transcripts. These results indicated that the short and distorted trichomes in Al might be caused by the suppression and alternative splicing of AlNAP1. Together, we revealed that AlNAP1 mediates trichome development and is a good candidate target for genetic modification of trichome length in kiwifruit.

A chitinase CsChi23 promoter polymorphism underlies cucumber resistance against Fusarium oxysporum f. sp. cucumerinum

Fusarium wilt disease, caused by Fusarium oxysporum f. sp. cucumerinum (Foc), leads to widespread yield loss and quality decline in cucumber. However, the molecular mechanisms underlying Foc resistance remain poorly understood. We report the mapping and functional characterisation of CsChi23, encoding a cucumber class I chitinase with antifungal properties. We assessed sequence variations at CsChi23 and the associated defence response against Foc. We functionally characterised CsChi23 using transgenic assay and expression analysis. The mechanism regulating CsChi23 expression was assessed by genetic and molecular approaches. CsChi23 was induced by Foc infection, which led to rapid upregulation in resistant cucumber lines. Overexpressing CsChi23 enhanced fusarium wilt resistance and reduced fungal biomass accumulation, whereas silencing CsChi23 causes loss of resistance. CsHB15, a homeodomain leucine zipper (HD-Zip) III transcription factor, was found to bind to the CsChi23 promoter region and activate its expression. Furthermore, silencing of CsHB15 reduces CsChi23 expression. A single-nucleotide polymorphism variation -400 bp upstream of CsChi23 abolished the HD-Zip III binding site in a susceptible cucumber line. Collectively, our study indicates that CsChi23 is sufficient to enhance fusarium wilt resistance and reveals a novel function of an HD-Zip III transcription factor in regulating chitinase expression in cucumber defence against fusarium wilt.

Single-cell transcriptome of Nepeta tenuifolia leaves reveal differentiation trajectories in glandular trichomes

The peltate glandular trichomes (PGTs) on Nepeta tenuifolia leaves can secrete and store bioactive essential oils. ScRNA-seq is a powerful tool for uncovering heterogeneous cells and exploring the development and differentiation of specific cells. Due to leaves rich in PGTs, the young leaves were used to isolated protoplasts and successfully captured 33,254 protoplasts for sequencing purposes. After cell type annotation, all the cells were partitioned into six broad populations with 19 clusters. Cells from PGTs were identified based on the expression patterns of trichome-specific genes, monoterpene biosynthetic genes, and metabolic analysis of PGT secretions. The developmental trajectories of PGTs were delineated by pseudotime analysis. Integrative analysis of scRNA-seq data from N. tenuifolia leaves and Arabidopsis thaliana shoot revealed that PGTs were specific to N. tenuifolia. Thus, our results provide a promising basis for exploring cell development and differentiation in plants, especially glandular trichome initiation and development.

Identification and characterization of a novel gene involved in glandular trichome development in Nepeta tenuifolia

Nepeta tenuifolia is a medicinal plant rich in terpenoids and flavonoids with antiviral, immunoregulatory, and anti-inflammatory activities. The peltate glandular trichome (PGT) is a multicellular structure considered to be the primary storage organ for monoterpenes; it may serve as an ideal model for studying cell differentiation and the development of glandular trichomes (GTs). The genes that regulate the development of GTs have not yet been well studied. In this study, we identified NtMIXTA1, a GT development-associated gene from the R2R3 MYB SBG9 family. NtMIXTA1 overexpression in tobacco resulted in the production of longer and denser GTs. Virus-induced gene silencing of NtMIXTA1 resulted in lower PGT density, a significant reduction in monoterpene concentration, and the decreased expression of genes related to monoterpene biosynthesis. Comparative transcriptome and widely targeted metabolic analyses revealed that silencing NtMIXTA1 significantly influenced the expression of genes, and the production of metabolites involved in the biosynthesis of terpenoids, flavonoids, and lipids. This study provides a solid foundation describing a mechanism underlying the regulation of GT development. In addition, this study further deepens our understanding of the regulatory networks involved in GT development and GT development-associated metabolite flux, as well as provides valuable reference data for studying plants with a high medicinal value without genetic transformation.

Molecular Mechanisms of Plant Trichome Development

Plant trichomes, protrusions formed from specialized aboveground epidermal cells, provide protection against various biotic and abiotic stresses. Trichomes can be unicellular, bicellular or multicellular, with multiple branches or no branches at all. Unicellular trichomes are generally not secretory, whereas multicellular trichomes include both secretory and non-secretory hairs. The secretory trichomes release secondary metabolites such as artemisinin, which is valuable as an antimalarial agent. Cotton trichomes, also known as cotton fibers, are an important natural product for the textile industry. In recent years, much progress has been made in unraveling the molecular mechanisms of trichome formation in Arabidopsis thaliana, Gossypium hirsutum, Oryza sativa, Cucumis sativus, Solanum lycopersicum, Nicotiana tabacum, and Artemisia annua. Here, we review current knowledge of the molecular mechanisms underlying fate determination and initiation, elongation, and maturation of unicellular, bicellular and multicellular trichomes in several representative plants. We emphasize the regulatory roles of plant hormones, transcription factors, the cell cycle and epigenetic modifications in different stages of trichome development. Finally, we identify the obstacles and key points for future research on plant trichome development, and speculated the development relationship between the salt glands of halophytes and the trichomes of non-halophytes, which provides a reference for future studying the development of plant epidermal cells.