Updates on molecular mechanisms in the development of branched trichome in Arabidopsis and nonbranched in cotton

GhHDZ76, a cotton HD-Zip transcription factor, involved in regulating the initiation and early elongation of cotton fiber development in G. hirsutum

In this study, the whole HD-Zip family members of G. hirsutum were identified, and GhHDZ76 was classified into the HD-Zip IV subgroup. GhHDZ76 was predominantly expressed in the 0-5 DPA of fiber development stage and localized in the nucleus. Overexpression of GhHDZ76 significantly increased the length and density of trichomes in Arabidopsis thaliana. The fiber length of GhHDZ76 knockout lines by CRISPR/Cas9 was significantly shorter than WT at the early elongation and mature stage, indicating that GhHDZ76 positively regulate the fiber elongation. Scanning electron microscopy showed that the number of ovule surface protrusion of 0 DPA of GhHDZ76 knockout lines was significantly lower than WT, suggesting that GhHDZ76 can also promote the initiation of fiber development. The transcript level of GhWRKY16, GhRDL1, GhEXPA1 and GhMYB25 genes related to fiber initiation and elongation in GhHDZ76 knockout lines were significantly decreased. Yeast two-hybrid and Luciferase complementation imaging (LCI) assays showed that GhHDZ76 can interact with GhWRKY16 directly. As a transcription factor, GhHDZ76 has transcriptional activation activity, which could bind to L1-box elements of the promoters of GhRDL1 and GhEXPA1. Double luciferase reporter assay showed that the GhWRKY16 could enhance the transcriptional activity of GhHDZ76 to pGhRDL1, but it did not promote the transcriptional activity of GhHDZ76 to pGhEXPA1. GhHDZ76 protein may also promote the transcriptional activity of GhWRKY16 to the downstream target gene GhMYB25. Our results provided a new gene resource for fiber development and a theoretical basis for the genetic improvement of cotton fiber quality.

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

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

Genome-wide identification of the key kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L.)

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (G. hirsutum L.). Results showed that 159 kinesin genes, including 15 genes of the kinesin-13 gene subfamily, were identified in upland cotton; of which 157 kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 kinesin genes in upland cotton, including 10 kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 kinesin genes were significantly associated with three fiber traits, among which a kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

Graphene enhances artemisinin production in the traditional medicinal plant Artemisia annua via dynamic physiological processes and miRNA regulation

We investigated the effects of graphene on the model herb Artemisia annua, which is renowned for producing artemisinin, a widely used pharmacological compound. Seedling growth and biomass were promoted when A. annua was cultivated with low concentrations of graphene, an effect which was attributed to a 1.4-fold increase in nitrogen uptake, a 15%-22% increase in chlorophyll fluorescence, and greater abundance of carbon cycling-related bacteria. Exposure to 10 or 20 mg/L graphene resulted in a ∼60% increase in H2O2, and graphene could act as a catalyst accelerator, leading to a 9-fold increase in catalase (CAT) activity in vitro and thereby maintaining reactive oxygen species (ROS) homeostasis. Importantly, graphene exposure led to an 80% increase in the density of glandular secreting trichomes (GSTs), in which artemisinin is biosynthesized and stored. This contributed to a 5% increase in artemisinin content in mature leaves. Interestingly, expression of miR828 was reduced by both graphene and H2O2 treatments, resulting in induction of its target gene AaMYB17, a positive regulator of GST initiation. Subsequent molecular and genetic assays showed that graphene-induced H2O2 inhibits micro-RNA (miRNA) biogenesis through Dicers and regulates the miR828-AaMYB17 module, thus affecting GST density. Our results suggest that graphene may contribute to yield improvement in A. annua via dynamic physiological processes together with miRNA regulation, and it may thus represent a new cultivation strategy for increasing yield capacity through nanobiotechnology.

The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum

KEY MESSAGE

GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

Identification of bZIP Transcription Factors That Regulate the Development of Leaf Epidermal Cells in Arabidopsis thaliana by Single-Cell RNA Sequencing

Epidermal cells are the main avenue for signal and material exchange between plants and the environment. Leaf epidermal cells primarily include pavement cells, guard cells, and trichome cells. The development and distribution of different epidermal cells are tightly regulated by a complex transcriptional regulatory network mediated by phytohormones, including jasmonic acid, and transcription factors. How the fate of leaf epidermal cells is determined, however, is still largely unknown due to the diversity of cell types and the complexity of their regulation. Here, we characterized the transcriptional profiles of epidermal cells in 3-day-old true leaves of Arabidopsis thaliana using single-cell RNA sequencing. We identified two genes encoding BASIC LEUCINE-ZIPPER (bZIP) transcription factors, namely bZIP25 and bZIP53, which are highly expressed in pavement cells and early-stage meristemoid cells. Densities of pavement cells and trichome cells were found to increase and decrease, respectively, in bzip25 and bzip53 mutants, compared with wild-type plants. This trend was more pronounced in the presence of jasmonic acid, suggesting that these transcription factors regulate the development of trichome cells and pavement cells in response to jasmonic acid.

GhMYB30-GhMUR3 affects fiber elongation and secondary wall thickening in cotton

Xyloglucan, an important hemicellulose, plays a crucial role in maintaining cell wall structure and cell elongation. However, the effects of xyloglucan on cotton fiber development are not well understood. GhMUR3 encodes a xyloglucan galactosyltransferase that is essential for xyloglucan synthesis and is highly expressed during fiber elongation. In this study, we report that GhMUR3 participates in cotton fiber development under the regulation of GhMYB30. Overexpression GhMUR3 affects the fiber elongation and cell wall thickening. Transcriptome showed that the expression of genes involved in secondary cell wall synthesis was prematurely activated in OE-MUR3 lines. In addition, GhMYB30 was identified as a key regulator of GhMUR3 by Y1H, Dual-Luc, and electrophoretic mobility shift assay (EMSA) assays. GhMYB30 directly bound the GhMUR3 promoter and activated GhMUR3 expression. Furthermore, DAP-seq of GhMYB30 was performed to identify its target genes in the whole genome. The results showed that many target genes were associated with fiber development, including cell wall synthesis-related genes, BR-related genes, reactive oxygen species pathway genes, and VLCFA synthesis genes. It was demonstrated that GhMYB30 may regulate fiber development through multiple pathways. Additionally, GhMYB46 was confirmed to be a target gene of GhMYB30 by EMSA, and GhMYB46 was significantly increased in GhMYB30-silenced lines, indicating that GhMYB30 inhibited GhMYB46 expression. Overall, these results revealed that GhMUR3 under the regulation of GhMYB30 and plays an essential role in cotton fiber elongation and secondary wall thickening. Additionally, GhMYB30 plays an important role in the regulation of fiber development and regulates fiber secondary wall synthesis by inhibiting the expression of GhMYB46.

Genome-wide identification of the key Kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L)

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of Kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of Kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key Kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (Gossypium hirsutum L.). Results showed that 159 Kinesin genes, including 15 genes of the Kinesin-13 gene subfamily, were identified in upland cotton; of which 157 Kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 Kinesin genes in upland cotton, including 10 Kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven Kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 Kinesin genes were significantly associated with three fiber traits, among which a Kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one Kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the Kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited Kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

Phytochemical analysis and antifungal potential of two Launaea mucronata (Forssk.) Muschl and Launaea nudicaulis (L.) Hook.fil. wildly growing in Anbar province, Iraq

ETHNOPHARMACOLOGICAL RELEVANCE

Plant fungi are a serious problem in agriculture. Even though synthetic fungicides are an efficient control method, several negative side effects emerge from their extensive use, such as health problems, environmental pollution, and the emergence of resistant microorganisms. Thus, it is becoming more and more urgent to search for alternative control methods.

AIM OF THE STUDY

The aim of our study was to analyze phytochemical composition and antifungal potential of Launaea mucronata (Forssk.) Muschl. and Launaea nudicaulis (L.) Hook. fil. wildly growing in Anbar province, Iraq. In addition, molecular analysis was used to identify the plants species and molecular docking analysis was investigated between the major phytochemicals present in these plants and three selected fungal proteins, in order to assess the antifungal activity of the selected biochemicals against these proteins.

MATERIALS AND METHODS

Molecular analysis was performed using ITS sequencing protocol. The phytochemical analysis was done using GC-MS technique, while molecular docking analysis was performed by FRED application between selected compounds from each plant and three enzymes: 17β-hydroxysteroid dehydrogenase, endochitinase, and 14-α-demethylase. Finally, the antifungal activity was assessed by measuring inhibition percentage of Fusarium solani and Macrophomina phaseolina growth treated with ethanomethanolic extract of each plant.

RESULTS

Molecular analysis identified the selected plants as L. mucronata and L. nudicaulis, with an ITS region of 600 bp. Phytochemical analysis of Launaea spp. reported the presence of 35 compounds in each ethanomethanolic extract, belonging to different classes. L. mucronata was mainly formed of lupeol (9.33%), whereas L. nudicaulis extract was dominated by the heterocyclic compound 4-(3-methoxyphenoxy)-1,2,5-oxadiazol-3-amine (20.2%). Furthermore, molecular docking analysis showed that 4H-pyran-4-one,2,3-dihydro-3,5-dihydroxy-6-methyl from L. mucronata and gulonic acid Ƴ-lactone from L. nudicaulis possessed the highest affinity score to 17-β-hydroxysteroid dehydrogenase (-4.584 and -7.811 kcal/mol, respectively). Sucrose from L. mucronata and glutaric acid, di(3,4-difluorobenzyl) ester from L. nudicaulis gave the highest affinity to endochitinase (-7.979 and - 8.446 kcal/mol, respectively). In addition, sterol 14-α-demethylase was affinitive to sucrose from L. mucronata and glutaric acid, di(3,4-difluorobenzyl) ester from L. nudicaulis via energetic score of -10.002 and -9.582 kcal/mol, respectively. Both extracts exhibited antifungal activity against F. solani and M. phaseolina in a dose-dependent manner.

CONCLUSIONS

This study confirms the antifungal potential of both Launaea spp. explained by the presence of phytochemicals with antimicrobial properties. These compounds have potential to be used as new drugs to treat infectious diseases caused by pathogens. Consequently, plants from Launaea genus could be a raw material for many studies such as therapeutic, taxonomical, drug modelling, and antifungal agent.

AsHSP26.2, a creeping bentgrass chloroplast small heat shock protein positively regulates plant development

KEY MESSAGE

The creeping bentgrass small heat shock protein AsHSP26.2 positively regulates plant growth and is a novel candidate for use in crop genetic engineering for enhanced biomass production and grain yield. Small heat shock proteins (sHSPs), a family of proteins with high level of diversity, significantly influence plant stress tolerance and plant development. We have cloned a creeping bentgrass chloroplast-localized sHSP gene, AsHSP26.2 responsive to IAA, GA and 6-BA stimulation. Transgenic creeping bentgrass overexpressing AsHSP26.2 exhibited significantly enhanced plant growth with increased stolon number and length as well as enlarged leaf blade width and leaf sheath diameters, but inhibited leaf trichomes initiation and development in the abaxial epidermis. These phenotypes are completely opposite to those displayed in the transgenic plants overexpressing AsHSP26.8, another chloroplast sHSP26 isoform that contains additional seven amino acids (AEGQGDG) between the consensus regions III and IV (Sun et al., Plant Cell Environ 44:1769-1787, 2021). Furthermore, AsHSP26.2 overexpression altered phytohormone biosynthesis and signaling transduction, resulting in elevated auxin and gibberellins (GA) accumulation. The results obtained provide novel insights implicating the sHSPs in plant growth and development regulation, and strongly suggest AsHSP26.2 to be a novel candidate for use in crop genetic engineering for enhanced plant biomass production and grain yield.

Beyond skin-deep: targeting the plant surface for crop improvement

The above-ground plant surface is a well-adapted tissue layer that acts as an interface between the plant and its surrounding environment. As such, its primary role is to protect against desiccation and maintain the gaseous exchange required for photosynthesis. Further, this surface layer provides a barrier against pathogens and herbivory, while attracting pollinators and agents of seed dispersal. In the context of agriculture, the plant surface is strongly linked to post-harvest crop quality and yield. The epidermal layer contains several unique cell types adapted for these functions, while the non-lignified above-ground plant organs are covered by a hydrophobic cuticular membrane. This review aims to provide an overview of the latest understanding of the molecular mechanisms underlying crop cuticle and epidermal cell formation, with focus placed on genetic elements contributing towards quality, yield, drought tolerance, herbivory defence, pathogen resistance, pollinator attraction, and sterility, while highlighting the inter-relatedness of plant surface development and traits. Potential crop improvement strategies utilizing this knowledge are outlined in the context of the recent development of new breeding techniques.

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.

AXR1 modulates trichome morphogenesis through mediating ROP2 stability in Arabidopsis

Cell differentiation and morphogenesis are crucial for the establishment of diverse cell types and organs in multicellular organisms. Trichome cells offer an excellent paradigm for dissecting the regulatory mechanisms of plant cell differentiation and morphogenesis due to their unique growth characteristics. Here, we report the isolation of an Arabidopsis mutant, aberrantly branched trichome 3-1 (abt3-1), with a reduced trichome branching phenotype. Positional cloning and molecular complementation experiments confirmed that abt3-1 is a new mutant allele of Auxin resistant 1 (AXR1), which encodes the N-terminal half of ubiquitin-activating enzyme E1 and functions in auxin signaling pathway. Meanwhile, we found that transgenic plants expressing constitutively active version of ROP2 (CA-ROP2) caused a reduction of trichome branches, resembling that of abt3-1. ROP2 is a member of Rho GTPase of plants (ROP) family, serving as versatile signaling switches involved in a range of cellular and developmental processes. Our genetic and biochemical analyses showed AXR1 genetically interacted with ROP2 and mediated ROP2 protein stability. The loss of AXR1 aggravated the trichome defects of CA-ROP2 and induced the accumulation of steady-state ROP2. Consistently, elevated AXR1 expression levels suppressed ROP2 expression and partially rescued trichome branching defects in CA-ROP2 plants. Together, our results presented a new mutant allele of AXR1, uncovered the effects of AXR1 and ROP2 during trichome development, and revealed a pathway of ROP2-mediated regulation of plant cell morphogenesis in Arabidopsis.

HECT-type ubiquitin ligase KAKTUS mediates the proteasome-dependent degradation of cyclin-dependent kinase inhibitor KRP2 during trichome morphogenesis in Arabidopsis

Trichome development is a fascinating model to elaborate the plant cell differentiation and growth processes. A wealth of information has pointed to the contributions of the components associated with cell cycle control and ubiquitin/26S proteasome system (UPS) to trichome morphogenesis, but how these two pathways are connected remains obscure. Here, we report that HECT‐type ubiquitin ligase KAKTUS (KAK) targets the cyclin‐dependent kinase (CDK) inhibitor KRP2 (for kip‐related protein 2) for proteasome‐dependent degradation during trichome branching in Arabidopsis. We show that over‐expression of KRP2 promotes trichome branching and endoreduplication which is similar to kak loss of function mutants. KAK directly interacts with KRP2 and mediates KRP2 degradation. Mutation of KAK results in the accumulation of steady‐state KRP2. Consistently, in kak pKRP2:KRP2‐GFP plants, the trichome branching is further induced compared with the single mutant. Taken together, our studies bridge the cell cycle control and UPS pathways during trichome development and underscore the importance of post‐translational control in epidermal differentiation.

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.

Cannabis sativa: origin and history, glandular trichome development, and cannabinoid biosynthesis

Is Cannabis a boon or bane? Cannabis sativa has long been a versatile crop for fiber extraction (industrial hemp), traditional Chinese medicine (hemp seeds), and recreational drugs (marijuana). Cannabis faced global prohibition in the twentieth century because of the psychoactive properties of ∆9-tetrahydrocannabinol; however, recently, the perspective has changed with the recognition of additional therapeutic values, particularly the pharmacological potential of cannabidiol. A comprehensive understanding of the underlying mechanism of cannabinoid biosynthesis is necessary to cultivate and promote globally the medicinal application of Cannabis resources. Here, we comprehensively review the historical usage of Cannabis, biosynthesis of trichome-specific cannabinoids, regulatory network of trichome development, and synthetic biology of cannabinoids. This review provides valuable insights into the efficient biosynthesis and green production of cannabinoids, and the development and utilization of novel Cannabis varieties.

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.

CRISPR/Cas9-mediated mutagenesis of the mediator complex subunits MED5a and MED5b genes impaired secondary metabolite accumulation in hop (Humulus lupulus)

Hop (Humulus lupulus L.) is an important commercial crop known for the biosynthesis of valuable specialized secondary metabolites in glandular trichomes (lupulin glands), which are used for the brewing industry. To achieve burgeoning market demands is the essentiality of comprehensive understanding of the mechanisms of biosynthesis of secondary metabolites in hop. Over the past year, several studies using structural biology and functional genomics approaches have shown that Mediator (MED) serves as an integrative hub for RNAP II-mediated transcriptional regulation of various physiological and cellular processes, including involvement of MED5a and MED5b in hyperaccumulation of phenylpropanoid in A. thaliana. In the present work, an unprecedented attempt was made to generate Hlmed5a/med5b double loci mutant lines in hop using a CRISPR/Cas9-based genome editing system. The Hlmed5a/med5b double loci mutant lines showed reduced expression of structural genes of the flavonoid, humulone, and terpenoid biosynthetic pathways, which was more pronounced in the lupulin gland compared to leaf tissue and was consistent with their reduced accumulation. Phenotypic and anatomical observations revealed that Hlmed5a/med5b double loci mutant line exhibited robust growth, earlier flowering, earlier cone maturity, reduced cone size, variations in floral structure patterns, and distorted lupulin glands without any remarkable changes in leaf morphology, intensity of leaf color, and chlorophyll content. Comparative transcriptome analysis of leaf and lupulin gland tissues indicates that the expression of enzymatic genes related to secondary metabolite biosynthesis, phytohormone biosynthesis, floral organs, flowering time, and trichome development, including other genes related to starch and sucrose metabolism and defense mechanisms, were differentially modulated in the Hlmed5a/med5b lines. The combined results from functional and transcriptomic analyses illuminates the pivotal function of HlMED5a and HlMED5b in homeostasis of secondary meatbolites accumulation in hop.

GhSINA1, a SEVEN in ABSENTIA ubiquitin ligase, negatively regulates fiber development in Upland cotton

Protein ubiquitination is essential for plant growth and responses to the environment. The SEVEN IN ABSENTIA (SINA) ubiquitin ligases have been extensively studied in plants, but information on their roles in fiber development is limited. Here, we identified GhSINA1 in Upland cotton (Gossypium hirsutum), which has a conserved RING finger domain and SINA domain. Quantitative real-time PCR (qRT-PCR) analysis showed that GhSINA1 was preferentially expressed during fiber initiation and elongation, especially during initiation in the fuzzless-lintless cotton mutant. Subcellular localization experiments indicated that GhSINA1 localized to the nucleus. In vitro ubiquitination analysis revealed that GhSINA1 has E3 ubiquitin ligase activity. Ectopic overexpression of GhSINA1 in Arabidopsis thaliana reduced the number and length of root hairs and trichomes. Yeast two-hybrid (Y2H), firefly luciferase complementation imaging (LCI), and bimolecular fluorescence complementation (BiFC) assays demonstrated that the GhSINA1 proteins could interact with each other to form homodimers and heterodimers. Overall, these results suggest that GhSINA1 may act as a negative regulator in cotton fiber development through homodimerization and heterodimerization.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

Genome-wide analysis of the serine carboxypeptidase-like protein family reveals Ga09G1039 is involved in fiber elongation in cotton

The Gossypium is a model genus for understanding polyploidy and the evolutionary pattern of inheritance. This study aimed to investigate the characteristics of SCPLs in different cotton species and their role in fiber development. A total of 891 genes from one typical monocot and ten dicot species were naturally divided into three classes based on phylogenetic analysis. The SCPL gene family in cotton has undergone intense purifying selection with some functional variation. Segmental duplication and whole genome duplication were shown to be the two main reasons for the increase in the number of genes during cotton evolution. The identification of Gh_SCPL genes exhibiting differential expression in particular tissues or response to environmental stimuli provides a new measure for the in-depth characterization of selected genes of importance. Ga09G1039 was involved in the developmental process of fibers and ovules, and it is significantly different from proteins from other cotton species in terms of phylogenetic, gene structure, conserved protein motifs and tertiary structure. Overexpression of Ga09G1039 significantly increased the length of stem trichomes. Ga09G1039 may be a serine carboxypeptidase protein with hydrolase activity, according to functional region, prokaryotic expression, and western blotting analysis. The results provide a comprehensive overview of the genetic basis of SCPLs in Gossypium and further our knowledge in understanding the key aspects of SCPLs in cotton with their potential role in fiber development and stress resistance.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton

BACKGROUND

The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development.

RESULTS

We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content.

CONCLUSIONS

Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

Function and molecular mechanism of a poplar placenta limited MIXTA gene in regulating differentiation of plant epidermal cells

The placenta in fruits of most plants either desiccate and shrink as the fruits mature or develop further to form the fleshy tissues. In poplars, placental epidermal cells protrude collectively to produce catkin fibers. In this study, three carpel limited MIXTA genes, PdeMIXTA02, PdeMIXTA03, PdeMIXTA04, were find to specifically expressed in carpel immediately after pollination. Heterologous expression of the three genes in Arabidopsis demonstrated that PdeMIXTA04 significantly promoted trichomes density and could restore trichomes in the trichomeless mutant. By contrast, such functions were not observed with PdeMIXTA02, PdeMIXTA03. In situ hybridization revealed that PdeMIXTA04 was explicitly expressed in poplar placental epidermal cells. We also confirmed trichome-specific expression of the PdeMIXTA04 promoter. Multiple experimental proofs have confirmed the interaction between PdeMIXTA04, PdeMYC and PdeWD40, indicating PdeMIXTA04 functioned through the MYB-bHLH-WD40 ternary complex. Our work provided distinctive understanding of the molecular mechanism triggering differentiation of poplar catkins.

AaSEPALLATA1 integrates jasmonate and light-regulated glandular secretory trichome initiation in Artemisia annua

Glandular secretory trichomes (GSTs) can secrete and store a variety of specific metabolites. By increasing GST density, valuable metabolites can be enhanced in terms of productivity. However, the comprehensive and detailed regulatory network of GST initiation still needs further investigation. By screening a complementary DNA library derived from young leaves of Artemisia annua, we identified a MADS-box transcription factor, AaSEPALLATA1 (AaSEP1), that positively regulates GST initiation. Overexpression of AaSEP1 in A. annua substantially increased GST density and artemisinin content. The HOMEODOMAIN PROTEIN 1 (AaHD1)-AaMYB16 regulatory network regulates GST initiation via the jasmonate (JA) signaling pathway. In this study, AaSEP1 enhanced the function of AaHD1 activation on downstream GST initiation gene GLANDULAR TRICHOME-SPECIFIC WRKY 2 (AaGSW2) through interaction with AaMYB16. Moreover, AaSEP1 interacted with the JA ZIM-domain 8 (AaJAZ8) and served as an important factor in JA-mediated GST initiation. We also found that AaSEP1 interacted with CONSTITUTIVE PHOTOMORPHOGENIC 1 (AaCOP1), a major repressor of light signaling. In this study, we identified a MADS-box transcription factor that is induced by JA and light signaling and that promotes the initiation of GST in A. annua.

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.

Functional characterization of TBL genes revealed the role of GhTBL7 and GhTBL58 in cotton fiber elongation

TBL (Trichome Birefringence Like) gene family members are involved in trichome initiation and xylan acetylation in several plant species. In our research, we identified 102 TBLs from G. hirsutum. The phylogenetic tree classified TBL genes into five groups. Collinearity analysis of TBL genes indicated 136 paralogous gene pairs in G. hirsutum. Gene duplication indicated that WGD or segmental duplication contributed to the GhTBL gene family expansion. Promoter cis-elements of GhTBLs were related to growth and development, seed-specific regulation, light, and stress responses. GhTBL genes (GhTBL7, GhTBL15, GhTBL21, GhTBL25, GhTBL45, GhTBL54, GhTBL67, GhTBL72, and GhTBL77) exhibited upregulated response under exposure to cold, heat, NaCl, and PEG. GhTBL genes exhibited high expression during fiber development stages. Two GhTBL genes (GhTBL7 and GhTBL58) showed differential expression at 10 DPA fiber, as 10 DPA is a fast fiber elongation stage and fiber elongation is a very important stage of cotton fiber development. Subcellular localization of GhTBL7 and GhTBL58 revealed that these genes reside inside the cell membrane. Promoter GUS activity of GhTBL7 and GhTBL58 exhibited deep staining in roots. To further validate the role of these genes in cotton fiber elongation, we silenced these genes and observed a significant reduction in the fiber length at 10 DPA. In conclusion, the functional study of cell membrane-associated genes (GhTBL7 and GhTBL58) showed deep staining in root tissues and potential function during cotton fiber elongation at 10 DPA fiber.

Comparative transcriptome analysis reveals hormone, transcriptional and epigenetic regulation involved in prickle formation in Zanthoxylum armatum

Zanthoxylum armatum is an evergreen plant with high economical and medicinal values. The presence of prickles on stems and leaves is undesirable for them make picking difficult. To date, little is known of prickle formation in Z. armatum. Herein, the morphological and molecular features of prickle initiation in prickless (WC) and three types of prickly Z. armatum were characterized. Compared to WC, the levels of cytokinin and auxin were increased, while GA and JA declined in prickly Z. armatum. Transcriptome analysis identified 6258 differentially expressed genes (DEGs)between prickless and prickly Z. armatum. Among them, several DEGs related to hormone metabolism and signaling, including LOG7, CKX3, AHK1, three DELLAs, six JAZs and TIR1, were candidate genes involved in prickle formation. Transcription factors associated with prickle formation was screened, including MYB6-1/MYB6-2, WER, GL3-2, SPL4/5, SOC1, and SCL32. Of them, MYB6-1 and WER might negatively regulate prickles initiation via interacting with GL3-2. Additionally, the histone acetylation and DNA methylation levels, the transcripts of histone acetyltransferase/deacetylase and DNA methyltransferases showed significant differences between prickless and prickly plants, indicating their involvements in prickle initiation. These findings illustrate the regulation of prickle formation might be mediated by phytohormones (especially cytokinin), transcription factors and epigenetic modifications in Z. armatum.

Comparative Transcriptome Analysis of Gleditsia sinensis Thorns at Different Stages of Development

G. sinensis thorn (called "zào jiǎo cì", ZJC) has important medicinal and economic value, however, little is known about the molecular mechanisms behind the development of ZJC. In this study, we measured the content of soluble sugar and starch during the growth and development of the thorn, and performed transcriptome sequencing of the thorn segment, non-thorn segment, apex, and root tip at five distinct stages of thorn formation. The results showed that, with the growth of ZJC, the soluble sugar content of the roots, hypocotyls, thorn stems, thornless stems, leaves, and the starch content of the roots and leaves all firstly increased and then decreased after the basic structure of thorns was formed; the soluble sugar content and starch content of ZJC showed an overall downward trend (decreased by 59.26% and 84.56%, respectively). Myb-like, YABBY2, Growth-regulating factor 3, TCP2, Zinc transporter 8, and another 25 genes may be related to the maintenance and growth of thorns. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) between stems with thorn and thorn-free stems found that a significant number of DEGs were annotated with terms related to the positive regulation of development, heterochronic (GO:0045962), the positive regulation of photomorphogenesis (GO:2000306), and other biological process (BP) terms. The developmental initiation regulation of ZJC may be regulated by TCP transcription factors (TFs). Eight genes were selected randomly to validate the RNA-seq results using real-time quantitative PCR (RT-qPCR) and they indicated that the transcriptome data were reliable. Our work provided a comprehensive review of the thorn development of G. sinensis.

Genome-wide characterization of trichome birefringence-like genes provides insights into fiber yield improvement

Cotton is an important fiber crop. The cotton fiber is an extremely long trichome that develops from the epidermis of an ovule. The trichome is a general and multi-function plant organ, and trichome birefringence-like (TBL) genes are related to trichome development. At the genome-wide scale, we identified TBLs in four cotton species, comprising two cultivated tetraploids (Gossypium hirsutum and G. barbadense) and two ancestral diploids (G. arboreum and G. raimondii). Phylogenetic analysis showed that the TBL genes clustered into six groups. We focused on GH_D02G1759 in group IV because it was located in a lint percentage-related quantitative trait locus. In addition, we used transcriptome profiling to characterize the role of TBLs in group IV in fiber development. The overexpression of GH_D02G1759 in Arabidopsis thaliana resulted in more trichomes on the stems, thereby confirming its function in fiber development. Moreover, the potential interaction network was constructed based on the co-expression network, and it was found that GH_D02G1759 may interact with several genes to regulate fiber development. These findings expand our knowledge of TBL family members and provide new insights for cotton molecular breeding.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum

The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC₄F_3:5). A single segment substitution line (CSSL-106) screened from BC₆F₂ was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC₇F₂ individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

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.

Chromosomal-level genome assembly of Melastoma candidum provides insights into trichome evolution

Melastoma, consisting of ~100 species diversified in tropical Asia and Oceania in the past 1-2 million years, represents an excellent example of rapid speciation in flowering plants. Trichomes on hypanthia, twigs and leaves vary markedly among species of this genus and are the most important diagnostic traits for species identification. These traits also play critical roles in contributing to differential adaptation of these species to their own habitats. Here we sequenced the genome of M. candidum, a common, erect-growing species from southern China, with the aim to provide genomic insights into trichome evolution in this genus. We generated a high-quality, chromosome-level genome assembly of M. candidum, with the genome size of 256.2 Mb and protein-coding gene number of 40,938. The gene families specific to, and significantly expanded in Melastoma are enriched for GO terms related to trichome initiation and differentiation. We provide evidence that Melastoma and its sister genus Osbeckia have undergone two whole genome duplications (WGDs) after the triplication event (γ) shared by all core eudicots. Preferential retention of trichome development-related transcription factor genes such as C2H2, bHLH, HD-ZIP, WRKY, and MYB after both WGDs might provide raw materials for trichome evolution and thus contribute to rapid species diversification in Melastoma. Our study provides candidate transcription factor genes related to trichome evolution in Melastoma, which can be used to evolutionary and functional studies of trichome diversification among species of this genus.

Phytoremediation potential of Solanum viarum Dunal and functional aspects of their capitate glandular trichomes in lead, cadmium, and zinc detoxification

In the present scenario, remediation of heavy metals (HMs) contaminated soil has become an important work to be done for the well-being of human and their environment. Phytoremediation can be regarded as an excellent method in environmental technologies. The present contemporary research explores the Solanum viarum Dunal function as a potential accumulator of hazardous HMs viz. lead (Pb), cadmium (Cd), zinc (Zn), and their combination (CHM). On toxic concentrations of Pb, Cd, Zn, and their synergistic exposure, seeds had better germination percentage and their 90d old aerial tissues accumulated Pb, Cd, and Zn concentrations ranging from 44.53, 84.06, and 147.29 mg kg^-1 DW, respectively. Pattern of accumulation in roots was as Zn 70.08 > Pb 48.55 > Cd 42.21 mg kg^-1DW. Under HMs treatment, positive modulation in physiological performances, antioxidant activities suggested an enhanced tolerance along with higher membrane stability due to increased levels of lignin, proline, and sugar. Phenotypic variations were recorded in prickles and roots of 120 d old HM stressed plants, which are directly correlated with better acclimation. Interestingly, trichomes of the plant also showed HM accumulation. Later, SEM-EDX microanalysis suggested involvement of S. viarum capitate glandular trichomes as excretory organs for Cd and Zn. Thus, the present study provides an understanding of the mechanism that makes S. viarum to function as potent accumulator and provides information to generate plants to be used for phytoremediation.

Transcriptomic analysis reveals the key role of histone deacetylation via mediating different phytohormone signalings in fiber initiation of cotton

Background

Histone deacetylation is one of the most important epigenetic modifications and plays diverse roles in plant development. However, the detailed functions and mechanisms of histone deacetylation in fiber development of cotton are still unclear. HDAC inhibitors (HDACi) have been commonly used to study the molecular mechanism underlying histone deacetylation or to facilitate disease therapy in humans through hindering the histone deacetylase catalytic activity. Trichostatin A (TSA)—the most widely used HDACi has been extensively employed to determine the role of histone deacetylation on different developmental stages of plants.

Results

Through in vitro culture of ovules, we observed that exogenous application of TSA was able to inhibit the fiber initiation development. Subsequently, we performed a transcriptomic analysis to reveal the underlying mechanisms. The data showed that TSA treatment resulted in 4209 differentially expressed genes, which were mostly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, and carbon metabolism pathways. The phytohormone signal transduction pathways harbor the most differentially expressed genes. Deeper studies showed that some genes promoting auxin, Gibberellic Acid (GA) signaling were down-regulated, while some genes facilitating Abscisic Acid (ABA) and inhibiting Jasmonic Acid (JA) signaling were up-regulated after the TSA treatments. Further analysis of plant hormone contents proved that TSA significantly promoted the accumulation of ABA, JA and GA₃.

Conclusions

Collectively, histone deacetylation can regulate some key genes involved in different phytohormone pathways, and consequently promoting the auxin, GA, and JA signaling, whereas repressing the ABA synthesis and signaling to improve the fiber cell initiation. Moreover, the genes associated with energy metabolism, phenylpropanoid, and glutathione metabolism were also regulated by histone deacetylation. The above results provided novel clues to illuminate the underlying mechanisms of epigenetic modifications as well as related different phytohormones in fiber cell differentiation, which is also very valuable for the molecular breeding of higher quality cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00840-4.

Characterization of anatomical features, developmental roadmaps, and key genes of bamboo leaf epidermis

The exact developmental roadmaps of bamboo leaf epidermis and the regulating genes are largely unknown. In this study, we comprehensively investigated the morphological features of the leaf epidermis of bamboo, Pseudosasa japonica. We also established the developmental roadmaps of the abaxial epidermis along the linearly growing leaf. A variant of P. japonica, P. japonica var. tsutsumiana, with smaller stomata and higher stomata density, was identified. Further analysis revealed that the higher stomata density of the variant was due to the abnormal increase in stomata columns within the single stomata band. This abnormal development of stomata bands was observed as early as the guard mother cell stage in the leaf division zone (DZ). Interestingly, the developmental pattern of the single stomata was similar in P. japonica and the variant. Molecular data showed that PjDLT (Dwarf and Low Tillering) was significantly downregulated in leaves DZ of the variant. Overexpression of PjDLT in Arabidopsis and rice results in smaller plants with lower stomata density, whereas downregulation or mutation of OsDLT results in increased stomata density. Our results highlight the morphological features and developmental schedule of the leaf epidermis of bamboo and provide evidence that DLT plays an important role in regulating stomata in bamboo and rice.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Differential environmental and genomic architectures shape the natural diversity for trichome patterning and morphology in different Arabidopsis organs

Despite the adaptive and taxonomic relevance of the natural diversity for trichome patterning and morphology, the molecular and evolutionary mechanisms underlying these traits remain mostly unknown, particularly in organs other than leaves. In this study, we address the ecological, genetic and molecular bases of the natural variation for trichome patterning and branching in multiple organs of Arabidopsis (Arabidopsis thaliana). To this end, we characterized a collection of 191 accessions and carried out environmental and genome-wide association (GWA) analyses. Trichome amount in different organs correlated negatively with precipitation in distinct seasons, thus suggesting a precise fit between trichome patterning and climate throughout the Arabidopsis life cycle. In addition, GWA analyses showed small overlapping between the genes associated with different organs, indicating partly independent genetic bases for vegetative and reproductive phases. These analyses identified a complex locus on chromosome 2, where two adjacent MYB genes (ETC2 and TCL1) displayed differential effects on trichome patterning in several organs. Furthermore, analyses of transgenic lines carrying different natural alleles demonstrated that TCL1 accounts for the variation for trichome patterning in all organs, and for stem trichome branching. By contrast, two other MYB genes (TRY and GL1), mainly showed effects on trichome patterning or branching, respectively.

Comparative transcriptome analysis reveals key pathways and genes involved in trichome development in tea plant (Camellia sinensis)

Trichomes, which develop from epidermal cells, are considered one of the important characteristics of the tea plant [Camellia sinensis (L.) O. Kuntze]. Many nutritional and metabolomic studies have indicated the important contributions of trichomes to tea products quality. However, understanding the regulation of trichome formation at the molecular level remains elusive in tea plants. Herein, we present a genome-wide comparative transcriptome analysis between the hairless Chuyeqi (CYQ) with fewer trichomes and the hairy Budiaomao (BDM) with more trichomes tea plant genotypes, toward the identification of biological processes and functional gene activities that occur during trichome development. In the present study, trichomes in both cultivars CYQ and BDM were unicellular, unbranched, straight, and soft-structured. The density of trichomes was the highest in the bud and tender leaf periods. Further, using the high-throughput sequencing method, we identified 48,856 unigenes, of which 31,574 were differentially expressed. In an analysis of 208 differentially expressed genes (DEGs) encoding transcription factors (TFs), five may involve in trichome development. In addition, on the basis of the Gene Ontology (GO) annotation and the weighted gene co-expression network analysis (WGCNA) results, we screened several DEGs that may contribute to trichome growth, including 66 DEGs related to plant resistance genes (PRGs), 172 DEGs related to cell wall biosynthesis pathway, 29 DEGs related to cell cycle pathway, and 45 DEGs related to cytoskeleton biosynthesis. Collectively, this study provided high-quality RNA-seq information to improve our understanding of the molecular regulatory mechanism of trichome development and lay a foundation for additional trichome studies in tea plants.

Overexpression of a ceramide synthase gene,GhCS1, inhibits fiber cell initiation and elongation by promoting the synthesis of ceramides containing dihydroxy LCB and VLCFA

Cotton is an important natural fiber crop worldwide. Cotton fiber cell is regarded as an ideal material for studying the growth and development of plant cells. Sphingolipids are important components of biomembrane and bioactive molecules which participate in many processes such as plant growth, development regulation, stimulus sensing, and stress response. However, the functions of sphingolipids in the cotton fiber development are still unclear. In the present study, we identified a cotton ceramide synthase gene, GhCS1, which is predominantly expressed in fiber cell. The GhCS1 is located in the endoplasmic reticulum and has the conserved domains of ceramide synthase. Overexpression of GhCS1 gene inhibited both vegetative and reproductive growth in cotton. Importantly, the fiber cell initiation and elongation were severely inhibited when compared with control. Comparison of the sphingolipid profile in the 0-DPA (days past anthesis) ovule (with fiber cell) between control and transgenic cotton plants showed that the content of sphingosines (Sph) decreased significantly in transgenic ovules, whereas the content of phyto-sphingosines (Phyto-Sph) had no change. Meanwhile, the content of ceramide containing Sph and very-long-chain fatty acid (VLCFA) increased significantly in transgenic ovules, while ceramide containing Phyto-Sph and long-chain fatty acids (LCFA)/VLCFA significantly decreased. These results indicated that GhCS1 was a functional ceramide synthase, which preferentially used Sph and VLCFA as substrates and was different from the Arabidopsis ceramide synthase AtLOH1/AtLOH3, which preferentially used Phyto-Sph and VLCFA as substrates, and also different from AtLOH2, which preferentially used Sph and LCFA as substrates. It is suggested that GhCS1 might be a new ceramide synthase gene in the plant, play some roles in the development of fiber cells and cotton plants.

RrTTG1 promotes fruit prickle development through an MBW complex in Rosa roxburghii

Fruit prickles are widely distributed on the pericarp and exhibit polymorphic traits at different developmental stages. Although they are multicellular appendages that are well-known for helping plants defend against biotic and abiotic stresses, their origination and molecular mechanism are still less known. Here, we studied the origination and molecular mechanism of fruit prickles in Rosa roxburghii. Using morphological and histological observations, we found that the fruit prickle primordium of R. roxburghii originated from the ground meristem that underwent cell division to form flagelliform prickles, continued to enlarge, and finally lignified to form mature fruit prickles. We amplified a homolog of candidate gene TRANSPARENT TESTA GLABRA1 (TTG1) from R. roxburghii, named RrTTG1. RrTTG1 harbored four conserved WD-repeat domains and was exclusively nuclear-localized. Using qRT-PCR and in situ hybridization, we found that RrTTG1 was constitutively expressed and highly expressed during the initiation and cell expansion phases of fruit prickles. Ectopic expression analysis in Arabidopsis proved that RrTTG1 substantially enhanced the number of trichome and pigmentation production and inhibited root hair formation. Besides, RrTTG1 complemented the phenotypes of the ttg1 mutant in Arabidopsis, thus indicating that RrTTG1 played pleiotropic roles akin to AtTTG1. We demonstrated that the RrTTG1 only interacted with RrEGL3, a homolog of ENHANCER OF GLABRA3 (EGL3), via yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays. Briefly, RrTTG1 might positively regulate the initiation of fruit prickle primordium and cell enlargement by forming the RrTTG1-RrEGL3-RrGL1 complex in R. roxburghii. Therefore, our results help characterize the RrTTG1 in R. roxburghii and also elucidate the establishment of the prickles regulatory system in the Rosaceae plants.

Tomato ARPC1 regulates trichome morphology and density and terpene biosynthesis

Based on transcriptomic analysis of wild-type and mutant tomato plants, ARPC1 was found to be important for trichome formation and development and it plays a key role in terpene synthesis. Trichomes are protruding epidermal cells in plant species. They function as the first defense layer against biotic and abiotic stresses. Despite the essential role of tomato trichomes in defense against herbivores, the understanding of their development is still incomplete. Therefore, the aim of this study was to identify genes involved in trichome formation and morphology and terpene synthesis, using transcriptomic techniques. To achieve this, we examined leaf morphology and compared the expression levels of some putative genes involved in trichome formation between wild-type (WT) and hairless-3 (hl-3) tomato mutant. The hl-3 plants displayed swollen and distorted trichomes and reduced trichome density (type I and IV) and terpene synthesis compared with that of the WT plants. Gene expression analysis showed that Actin-Related Protein Component1 (ARPC1) was expressed more highly in the WT than in the hl-3 mutant, indicating its critical role in trichome morphology and density. Additionally, the expression of MYC1 and several terpene synthase genes (TPS9, 12, 20), which are involved in type VI trichome initiation and terpene synthesis, was lower in the hl-3 mutant than in the WT plants. Moreover, transformation of the hl-3 mutant with WT ARPC1 restored normal trichome structure and density, and terpene synthesis. Structural and amino acid sequence analysis showed that there was a missplicing mutation in the hl-3 mutant, which was responsible for the abnormal trichome structure and density, and impaired terpene synthesis. Overall, the findings of this study demonstrated that ARPC1 is involved in regulating trichome structure and terpene synthesis in tomato.

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.

Transcriptome mining of genes in Zanthoxylum armatum revealed ZaMYB86 as a negative regulator of prickly development

Zanthoxylum armatum DC. is an important economic tree species. Prickle is a type of trichome with special morphology, and there are a lot of prickles on the leaves of Z. armatum, which seriously restricts the development of Z. armatum industry. In this study, the leaves of Z. armatum cv. Zhuye (ZY) and its budding variety 'Rongchangwuci' (WC) (A less prickly mutant variety) at different developmental stages were used as materials, and the transcriptome sequencing data were analyzed. A total of 96,931 differentially expressed genes (DEGs) were identified among the samples, among which 1560 were candidate DEGs that might be involved in hormone metabolism. The contents of JA, auxin and CK phytohormones in ZY leaves were significantly higher than those in WC leaves. Combined with weighted gene co-expression network analysis, eight genes (MYC, IAA, ARF, CRE/AHK, PP2C, ARR-A, AOS and LOX) were identified, including 25 transcripts, which might affect the metabolism of the three hormones and indirectly participate in the formation of prickles. Combining with the proteins successfully reported in other plants to regulate trichome formation, ZaMYB86, a transcription factor of R2R3 MYB family, was identified through local Blast and phylogenetic tree analysis, which might regulate prickle formation of Z. armatum. Overexpression of ZaMYB86 in mutant A. thaliana resulted in the reduction of trichomes in A. thaliana leaves, which further verified that ZaMYB86 was involved in the formation of pickles. Yeast two-hybrid results showed that ZaMYB86 interacted with ZaMYB5. Furthermore, ZaMYB5 was highly homologous to AtMYB5, a transcription factor that regulated trichomes development, in MYB DNA binding domain. Taken together, these results indicated that ZaMYB86 and ZaMYB5 act together to regulate the formation of prickles in Z. armatum. Our findings provided a new perspective for revealing the molecular mechanism of prickly formation.

Diverse roles of MYB transcription factors in regulating secondary metabolite biosynthesis, shoot development, and stress responses in tea plants (Camellia sinensis)

Tea (Camellia sinensis) is concocted from tea plant shoot tips that produce catechins, caffeine, theanine, and terpenoids, which collectively determine the rich flavors and health benefits of the infusion. However, little is known about the integrated regulation of shoot tip development and characteristic secondary metabolite biosynthesis in tea plants. Here, we demonstrate that MYB transcription factors (TFs) play key and yet diverse roles in regulating leaf and stem development, secondary metabolite biosynthesis, and environmental stress responses in tea plants. By integrating transcriptomic and metabolic profiling data in different tissues at a series of developmental stages or under various stress conditions, alongside biochemical and genetic analyses, we predicted the MYB TFs involved in regulating shoot development (CsMYB2, 98, 107, and 221), epidermal cell initiation (CsMYB184, 41, 139, and 219), stomatal initiation (CsMYB113 and 153), and the biosynthesis of flavonoids (including catechins, anthocyanins, and flavonols; CsMYB8 and 99), caffeine (CsMYB85 and 86), theanine (CsMYB9 and 49), carotenoids (CsMYB110), mono-/sesquiterpenoid volatiles (CsMYB68, 147, 148, and 193), lignin (CsMYB164 and 192), and indolic compounds (CsMYB139, 162, and 198), as well as the MYB TFs that are likely involved in hormone signaling-mediated environmental stress and defense responses. We characterized the functions of some key MYBs in regulating flavonoid and carotenoid biosynthesis for tea quality and flavor. This study provides a cross-family analysis of MYBs in tea alongside new insights into the coordinated regulation of tea plant shoot development and secondary metabolism, paving the way towards understanding of tea quality trait formation and genetic improvement of quality tea plants.

Characterisation of the Mentha canadensis R2R3-MYB transcription factor gene McMIXTA and its involvement in peltate glandular trichome development

BACKGROUND

Mentha canadensis L. has important economic value for the production of essential oils, which are synthesised, secreted and stored in peltate glandular trichomes. As a typical multicellular secretory trichome, glandular trichomes are important biological factories for the synthesis of some specialised metabolites. However, little is known about the molecular mechanism of glandular trichome development in M. canadensis.

RESULTS

In this study, the R2R3-MYB transcription factor gene McMIXTA was isolated to investigate its function in glandular trichome development. Bioinformatics analysis indicated that McMIXTA belonged to the subgroup 9 R2R3-MYB, with a R2R3 DNA-binding domain and conserved subgroup 9 motifs. A subcellular localisation assay indicated that McMIXTA was localised in the nucleus. Transactivation analysis indicated that McMIXTA was a positive regulator, with transactivation regions located between positions N253 and N307. Yeast two-hybrid and bimolecular fluorescence complementation assays showed that McMIXTA formed a complex with McHD-Zip3, a trichome development-related HD-ZIP IV transcription factor. Overexpression of McMIXTA in Mentha × piperita L. caused an increase in peltate glandular trichomes density of approximately 25% on the leaf abaxial surface.

CONCLUSIONS

Our results demonstrated that the subgroup 9 R2R3-MYB transcription factor McMIXTA has a positive effect on regulating peltate glandular trichome development and the MIXTA/HD-ZIP IV complexes might be conserved regulators for glandular trichome initiation. These results provide useful information for revealing the regulatory mechanism of multicellular glandular trichome development.

Multi-Dimensional Molecular Regulation of Trichome Development in Arabidopsis and Cotton

Plant trichomes are specialized epidermal cells that are widely distributed on plant aerial tissues. The initiation and progression of trichomes are controlled in a coordinated sequence of multiple molecular events. During the past decade, major breakthroughs in the molecular understanding of trichome development were achieved through the characterization of various trichomes defective mutants and trichome-associated genes, which revealed a highly complex molecular regulatory network underlying plant trichome development. This review focuses on the recent millstone in plant trichomes research obtained using genetic and molecular studies, as well as 'omics' analyses in model plant Arabidopsis and fiber crop cotton. In particular, we discuss the latest understanding and insights into the underlying molecular mechanisms of trichomes formation at multiple dimensions, including at the chromatin, transcriptional, post-transcriptional, and post-translational levels. We summarize that the integration of multi-dimensional trichome-associated genes will enable us to systematically understand the molecular regulation network that landscapes the development of the plant trichomes. These advances will enable us to address the unresolved questions regarding the molecular crosstalk that coordinate concurrent and ordered the changes in cotton fiber initiation and progression, together with their possible implications for genetic improvement of cotton fiber.

Dynamic roles and intricate mechanisms of ethylene in epidermal hair development in Arabidopsis and cotton

Ethylene affects many aspects of plant growth and development, including root hairs and trichomes growth in Arabidopsis, as well as fiber development in cotton, though the underlying mechanism is unclear. In this article, we update the research progress associated with the main genes in ethylene biosynthesis and signaling pathway, and we propose a clear ethylene pathway based on genome-wide identification of homologues in cotton. Expression pattern analysis using transcriptome data revealed that some candidate genes may contribute to cotton fiber development through the ethylene pathway. Moreover, we systematically summarized the effects of ethylene on the development of epidermal hair and the underlying regulatory mechanisms in Arabidopsis. Based on the knowledge of ethylene-promoted cell differentiation, elongation, and development in different tissues or plants, we advised a possible regulatory network for cotton fiber development with ethylene as the hub. Importantly, we emphasized the roles of ethylene as an important node in regulating cotton vegetative growth, and stress resistance, and suggested utilizing multiple methods to subtly modify ethylene synthesis or signaling in a tissue or spatiotemporal-specific manner to clarify its exact effect on architecture, adaptability of the plant, and fiber development, paving the way for basic research and genetic improvement of the cotton crop.

PaNAC089 is a membrane-tethered transcription factor (MTTF) that modulates flowering, chlorophyll breakdown and trichome initiation

Flowering and senescence are essential developmental stages of green plants, which are governed by complex molecular regulatory networks. However, the connection between flowering regulation and senescence regulation in London plane tree (Platanus acerifolia ) remains unknown. In this study, we identified a gene PaNAC089 from London plane tree, which encodes a membrane-tethered transcription factor (MTTF) belonging to the NAC (NAM, ATAF1/2, CUC2) transcription factor family. We investigated the functions of PaNAC089 in the regulation of flowering and senescence through the analysis of expression profiles and transgenic phenotypes. Heterologous overexpression of ΔPaNAC089 delayed flowering and inhibited chlorophyll breakdown to produce dark green rosette leaves in Arabidopsis . In addition, the trichome density of rosette leaves was decreased in transgenic lines. In ΔPaNAC089 overexpression plants, a series of functional genes with inhibited expression were identified by quantitative real-time polymerase chain reaction (qRT-PCR), including genes that regulate flowering, chlorophyll decomposition, and trichome initiation. Furthermore, Δ PaNAC089 directly binds to the promoter of CONSTANS (CO ) and NON-YELLOWING2 (NYE2 ) in the yeast one-hybrid assay. Consistent with this, luciferase (LUC) transient expression assays also showed that Δ PaNAC089 could inhibit the activity of NYE2 . To summarise, our data suggests that PaNAC089 is an MTTF that modulates flowering, chlorophyll breakdown and trichome initiation.

NbJAZ3 is required for jasmonate-meditated glandular trichome development in Nicotiana benthamiana

Exogenous methyl jasmonate (MeJA) treatment induces glandular trichome development in Nicotiana benthamiana, but the function of JAZ proteins, acting as core repressors, and their downstream genes have not been clearly shown in plants. Here, a bioinformatics analysis of 71 JAZ genes from tobacco, Arabidopsis thaliana, and tomato was carried out and shown to share highly conserved domains. Then, the expression profile of 17 NbJAZs in different tissues was analyzed, and NbJAZ3 was highly expressed in trichome. Through transgenic technology, we demonstrated that the glandular trichome density of NbJAZ3-overexpression lines significantly decreased with lower expression levels of NbWo, NbCycB2, and NbMIXTA. In contrast, the trichome density of NbJAZ3 RNAi lines slightly increased with higher expression level of NbWo. Given the negative protein feedback regulation relationship between NbCycB2 and NbWo, we verified that MeJA induced NbWo expression. NbWo was a direct target gene of NbJAZ3 and further demonstrated that NbJAZ3 inhibited the transcriptional activation of NbCycB2 by NbWo. Together, our findings outline a novel JA-meditated glandular trichome development model consisting of the NbJAZ3-NbWo-NbCycB2 axis.

RNA-Seq with a novel glabrous-ZM24fl reveals some key lncRNAs and the associated targets in fiber initiation of cotton

BACKGROUND

Cotton fiber is an important natural resource for textile industry and an excellent model for cell biology study. Application of glabrous mutant cotton and high-throughput sequencing facilitates the identification of key genes and pathways for fiber development and cell differentiation and elongation. LncRNA is a type of ncRNA with more than 200 nt in length and functions in the ways of chromatin modification, transcriptional and post-transcriptional modification, and so on. However, the detailed lncRNA and associated mechanisms for fiber initiation are still unclear in cotton.

RESULTS

In this study, we used a novel glabrous mutant ZM24fl, which is endowed with higher somatic embryogenesis, and functions as an ideal receptor for cotton genetic transformation. Combined with the high-throughput sequencing, fatty acid pathway and some transcription factors such as MYB, ERF and bHLH families were identified the important roles in fiber initiation; furthermore, 3,288 lncRNAs were identified, and some differentially expressed lncRNAs were also analyzed. From the comparisons of ZM24_0 DPA vs ZM24_-2 DPA and fl_0 DPA vs ZM24_0 DPA, one common lncRNA MSTRG 2723.1 was found that function upstream of fatty acid metabolism, MBY25-mediating pathway, and pectin metabolism to regulate fiber initiation. In addition, other lncRNAs MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1 were also showed potential important roles in fiber development; and the co-expression analysis between lncRNAs and targets showed the distinct models of different lncRNAs and complicated interaction between lncRNAs in fiber development of cotton.

CONCLUSIONS

From the above results, a key lncRNA MSTRG 2723.1 was identified that might mediate some key genes transcription of fatty acid metabolism, MYB25-mediating pathway, and pectin metabolism to regulate fiber initiation of ZM24 cultivar. Co-expression analysis implied that some other important lncRNAs (e.g., MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1) were also showed the different regulatory model and interaction between them, which proposes some valuable clues for the lncRNAs associated mechanisms in fiber development.