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

Monocot-like leaf structure and trichome-water relations in early growth stages of the C3 plant sand rice (Agriophyllum squarrosum)

Sand rice (Agriophyllum squarrosum (L.) Moq.), a species of Amaranthaceae sensu lato, is an annual psammophyte found mainly on sand dunes in arid and semi-arid regions of Central Asia. Trichomes on the leaf surface are a key adaptive feature of sand rice, however, their ecophysiological functions remain to be elucidated. In this study, the anatomical analysis of a leaf reveals that the foliage is characterized by its isobilateral and amphistomatic nature. Together with an average δ13C value of -28.57 ‰, these findings indicate that sand rice is a C3 plant. Utilizing a trichome-defective mutant (Agriophyllum squarrosum trichomeless1, astcl1) identified from an EMS mutagenesis library by our group, alongside the wild type Shapotou (SPT), we observed that both epidermal surfaces on SPT leaves exhibit dense and dendritic trichomes, forming rough leaf boundary layers. Chlorophyll leaching analysis revealed that SPT has a higher epidermal permeability, with 86.32 % of chlorophylls leaching out within 30 minutes, compared to only 24.18 % for astcl1. The astcl1 leaf displayed a higher leaf wettability, lower water loss, and thicker cuticle compared to SPT leaves, which may be associated with trichomes serving as the primary pathway for nonstomatal water movement. This is supported by the observation that toluidine blue staining is more intense in trichomes than in pavement cells. Our findings reveal that the correlation between trichomes and the cuticle appears to be a critical factor in managing water loss and overall plant hydration, highlighting the significant role of trichomes in the adaptation of sand rice to desert environments.

A GDSL esterase/lipase gene GbGELP identified from a fiber micronaire QTL qMIC-A11 modulates cell elongation and fiber development

KEY MESSAGE

A fiber micronaire QTL qMIC-A11 was fine-mapped, and the GDSL esterase/lipase gene GbGELP was identified as the causal gene of the QTL. GbGELP modulates cell elongation and cotton fiber development. Fine mapping and map-based cloning of fiber micronaire (MIC)-related quantitative trait loci (QTL) have not been reported to date. Here, we utilized a G. hirsutum (Gh) acc. TM-1-G. barbadense (Gb) acc. Hai7124 introgression line CSSL47, which exhibits a significant decrease in MIC compared to TM-1, to cross with TM-1 and develop the F2 and F2:3 secondary segregating populations. Further, a stable MIC QTL qMIC-A11 was simultaneously detected in the F2 and F2:3 populations and anchored within a 407 kb region. Among them, GB_A11G1593 encoding a GDSL esterase/lipase, exhibited substantially higher expression levels at fiber elongation period in CSSL47 compared to TM-1, which was temporally identified as the causal gene for qMIC-A11 and named as GbGELP. The heterologous expression of GbGELP in Arabidopsis showed increased root length, root cell length, rosette leaf growth, and trichome density. However, knockdown of GbGELP homologs in CSSL47 significantly decreased the fiber length. Further investigation found that there was an A/T single-nucleotide polymorphism variation (SNP) in the promoter of GELP orthologs between CSSL47 and TM-1, which results in a differential CATTAAATT/CATTTAATT HAHR1-box cis-acting element, a binding site for the homeodomain-leucine zipper IV (HD-ZIP IV). GbGELP was regulated by a HD-ZIP IV transcription factor GhHDG2 via binding to the CATTAAATT element in the GbGELP promoter, while GhGELP could be activated due to GhHDG2 unable to bind the CATTTAATT element in the GhGELP promoter. The fine-mapped MIC QTL qMIC-A11, along with the causal gene GbGELP, will be utilized to improve the fiber quality in cotton breeding.

The DNA N6-adenine methylation target gene GhMAF1 promotes fiber initiation at the base of Gossypium ovules

DNA methylation consists of 5-methylcytosine and N6-methyl deoxyadenosine (6mA) and is crucial in plant development. However, its specific role and potential mechanism to initiate cotton fibers remain unclear. This study employed Oxford Nanopore Technologies (ONT) sequencing to analyze DNA methylation alterations in ZM24 and ZM24 fuzzless-lintless (ZM24fl) during fiber initiation. Our results indicated that DNA 6mA methylation exhibited the most remarkable difference among ovule samples at -2, 0, and 5 d post anthesis of ZM24 and ZM24fl. Subsequently, genes with significant changes in DNA 6mA methylation and transcription during fiber initiation were screened. We found that GhMAF1 displayed significant transcriptional upregulation and 6mA enrichment in its promoter, which could serve as a potential target for DNA 6mA in fiber initiation. Further, we knocked out GhMAF1 using CRISPR-Cas technology and demonstrated that GhMAF1 specifically promotes the initiation of fiber cells at the base of the ovule by mediating the downstream JAZ2/CPC-MML3/MML4 pathway. These findings unveil a novel spatial module of fiber cell initiation on the ovule surface that involves GhMAF1. Ultimately, this work provides significant knowledge for the regulatory network of DNA 6mA modification in fiber initiation to improve fiber yield and quality.

Genome-wide analysis of basic helix-loop-helix (bHLH) transcription factors in petunia and identification of the putative candidate member involved in floral volatile benzenoids/phenylpropanoids metabolism

The basic helix-loop-helix (bHLH) family, a prominent group of transcription factors, is involved in plant growth, development, and secondary metabolic processes. Petunia (Petunia hybrida), a beloved and widely cultivated garden flower, boasts a diverse array of varieties, some of which exude a captivating fragrance that has garnered immense popularity. The aromatic allure of petunias primarily stems from the presence of volatile benzenoids/phenylpropanoids, the principal floral scent compounds. But whether bHLH transcription factors regulate petunia floral scent compound synthesis is not clear. In this study, we sought to screen the putative candidate member of bHLH which can be involved in the biosynthesis of benzenoids/phenylpropanoids by examining 63 members of the petunia bHLH gene family. Phylogenetic analysis of the 63 petunia bHLH proteins them into 16 subgroups. Almost all bHLH members contained alkaline/helix-loop-helix domains. Based on the reported RNA sequencing data of P. hybrida 'Mitchell', 30 assembled sequences were mapped to the bHLH genes of P. axillaris. Further qRT-PCR assays suggested that PhbHLH19 might be the putative candidate member in the biosynthesis of benzenoids/phenylpropanoids. PhbHLH19 showed higher expression levels in the petal limb but was lowly expressed at the bud stage, with a rapid increase in the expression level when flowers opened. The expression of PhbHLH19 displayed a significant positive correlation with that of PhPAL2, and the yeast one-hybrid assay verified that PhbHLH19 can bind to the promoter of PhPAL2. Moreover, a dual-luciferase assay proved the transcriptional activation of PhbHLH19 on PhPAL2. These findings suggested that PhbHLH19 might be a putative candidate in the regulation of benzenoid/phenylpropanoid synthesis by activating PhPAL2 expression.

The E3 ubiquitin ligase CONSTITUTIVE PHOTOMORPHOGENIC 1 and transcription factors ELONGATED HYPOCOTYL 5 and ROOT HAIR DEFECTIVE6 integrate light signaling and root hair development

Light signaling plays a substantial role in regulating plant development, including the differentiation and elongation of single-celled tissue. However, the identity of the regulatory machine that affects light signaling on root hair cell (RHC) development remains unclear. Here, we investigated how darkness inhibits differentiation and elongation of RHC in Arabidopsis (Arabidopsis thaliana). We found that light promotes the growth and development of RHC. RNA-seq analysis showed that light signaling regulates the differentiation of RHC by promoting the expression of specific genes in the root epidermis associated with cell wall remodeling, jasmonic acid, auxin, and ethylene signaling pathways. Together, these genes integrate light and phytohormone signals with root hair (RH) development. Our investigation also revealed that the core light signal factor ELONGATED HYPOCOTYL 5 (HY5) directly interacts with the key RH development factor ROOT HAIR DEFECTIVE6 (RHD6), which promotes the transcription of RSL4. However, CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1) repressed the RHD6 function through the COP1-HY5 complex. Our genetic studies confirm associations between RHD6, HY5, and COP1, indicating that RHD6 largely depends on HY5 for RH development. Ultimately, our work suggests a central COP1-HY5-RHD6 regulatory module that integrates light signaling and RH development with several downstream pathways, offering perspectives to decipher single-celled RH development.

Genome Scan Analysis for Advancing Knowledge and Conservation Strategies of Primitivo Clones (Vitis vinifera L.)

The success of the Primitivo variety underscores the critical need for the managing of clone genetic conservation, utilization, and improvement. By combining genomic and environmental data, breeders can better predict the performance of varieties, thereby improving breeding efficiency and enabling more targeted development of high-quality grapevine cultivars. In this study, 35 Primitivo clones were analysed, including selected and certified clones that have been propagated over several years in Apulia. Genetic variability among the Primitivo clones was assessed through genotyping by sequencing. Using 38,387 filtered SNPs, pairwise identity-by-state (IBS) analysis demonstrated the uniqueness of the 35 clones (IBS < 0.75), indicating a high degree of variability among the samples. Genetic diversity analysis revealed three primary groups, which were differentiated based on geographic origin. The clones from Gioia del Colle were grouped into two distinct clusters, which aligns with the observed variation in grape-related traits. The fixation index (FST > 0.50) identified numerous loci putatively associated with stress responses and developmental traits, including genes involved in key plant biological processes, stress response regulation, and adaptation to environmental conditions such as glutamate receptors, auxin, and ethylene signalling.

AaMYC3 bridges the regulation of glandular trichome density and artemisinin biosynthesis in Artemisia annua

Artemisinin, the well-known natural product for treating malaria, is biosynthesised and stored in the glandular-secreting trichomes (GSTs) of Artemisia annua. While numerous efforts have clarified artemisinin metabolism and regulation, the molecular association between artemisinin biosynthesis and GST development remains elusive. Here, we identified AaMYC3, a bHLH transcription factor of A. annua, induced by jasmonic acid (JA), which simultaneously regulates GST density and artemisinin biosynthesis. Overexpressing AaMYC3 led to a substantial increase in GST density and artemisinin accumulation. Conversely, in the RNAi-AaMYC3 lines, both GST density and artemisinin content were markedly reduced. Through RNA-seq and analyses conducted both in vivo and in vitro, AaMYC3 not only directly activates AaHD1 transcription, initiating GST development, but also up-regulates the expression of artemisinin biosynthetic genes, including CYP71AV1 and ALDH1, thereby promoting artemisinin production. Furthermore, AaMYC3 acts as a co-activator, interacting with AabHLH1 and AabHLH113, to trigger the transcription of two crucial enzymes in the artemisinin biosynthesis pathway, ADS and DBR2, ultimately boosting yield. Our findings highlight a critical connection between GST initiation and artemisinin biosynthesis in A. annua, providing a new target for molecular design breeding of traditional Chinese medicine.

Association analyses reveal both anthropic and environmental selective events during eggplant domestication

Eggplant (Solanum melongena) is one of the four most important Solanaceous crops, widely cultivated and consumed in Asia, the Mediterranean basin, and Southeast Europe. We studied the genome-wide association of historical genebank phenotypic data on a genotyped worldwide collection of 3449 eggplant accessions. Overall, 334 significant associations for key agronomic traits were detected. Significant correlations were obtained between different types of phenotypic data, some of which were not obvious, such as between fruit size/yield and fruit color components, suggesting simultaneous anthropic selection for genetically unrelated traits. Anthropic selection of traits like leaf prickles, fruit color, and yield, acted on distinct genomic regions in the two domestication centers (India and Southeast Asia), further confirming the multiple domestication of eggplant. To discriminate anthropic from environmental selection in domestication centers, we conducted a genotype-environment association (GEA) on a subset of georeferenced accessions from the Indian subcontinent. The population structure in this area revealed four genetic clusters, corresponding to a latitudinal gradient, and environmental factors explained 31% of the population structure when the effect of spatial distances was removed. GEA and outlier association identified 305 candidate regions under environmental selection, containing genes for abiotic stress responses, plant development, and flowering transition. Finally, in the Indian domestication center anthropic and environmental selection acted largely independently, and on different genomic regions. These data allow a better understanding of the different effects of environmental and anthropic selection during domestication of a crop, and the different world regions where some traits were initially selected by humans.

PbMYB5 transcription factor plays a role in regulating anthocyanin biosynthesis in pear (Pyrus bretschneideri Rehd) skin

The phenylacetone pathway, which encompasses flavonoids, lignin, and other compounds, is of paramount importance in determining the quality of pear fruit. Nevertheless, the precise regulatory functions of R2R3-MYB transcription factors in the metabolic pathways that regulate pear color changes remain unclear. In this study, we isolated an R2R3-PbMYB5(PbMYB5) transcription factor from 'Red Zaosu' pears and demonstrated that it influenced the expression of several genes, including PbCAD1, PbF5H, PbLAR, PbANR, and PbUFGT. The overexpression of PbMYB5 resulted in a notable elevation in anthocyanin concentration within the pear epidermis. Further research has shown that PbMYB5 is able to bind to PbANS and also has interactions with PbbHLH3 and PbbHLH33.We proposed that PbMYB5 forms a complex with PbbHLH3, PbbHLH33, and PbWD40 to activate PbANS and promote anthocyanin accumulation. This study offers new insights into the regulation of various metabolic pathways that impact fruit coloration.

Regulatory mechanisms of trichome and root hair development in Arabidopsis

In plants, cell fate determination is regulated temporally and spatially via a complex of signals consisting of a large number of genetic interactions. Trichome and root hair formation are excellent models for studying cell fate determination in plants. Nowadays, the mysteries underlying the reprograming of trichome and root hair and how nature programs the development of trichome and root hair is an interesting topic in the scientific field. In this review, we discuss the spatial and temporal regulatory networks and cross-talk between phytohormones and epigenetic modifications in the regulation of trichome and root hair initiation in Arabidopsis. The discussion in this review provides a good model for understanding the regulatory mechanism of cell differentiation processes in plants. Moreover, we summarize recent advances in the modulation of trichome and root hair initiation in plants and compare different regulatory mechanisms to help illuminate key goals for future research.

GLABRA3-mediated trichome branching requires transcriptional repression of MICROTUBULE-DESTABILIZING PROTEIN25

Microtubules play pivotal roles in establishing trichome branching patterns, which is a model system for studying cell-shape control in Arabidopsis (Arabidopsis thaliana). However, the signaling pathway that regulates microtubule reorganization during trichome branching remains poorly understood. In this study, we report that MICROTUBULE-DESTABILIZING PROTEIN25 (MDP25) is involved in GLABRA3 (GL3)-mediated trichome branching by regulating microtubule stability. Loss of MDP25 function led to excessive trichome branching, and this phenotype in mdp25 could not be rescued by the MDP25 K7A or MDP25 K18A mutated variants. Pharmacological treatment and live-cell imaging revealed increased microtubule stability in the mdp25 mutant. Furthermore, the microtubule collar observed during trichome branching remained more intact in mdp25 compared to the WT under oryzalin treatment. Results of genetic assays further demonstrated that knocking out MDP25 rescued the reduced branching phenotype of gl3 trichomes. In gl3 trichomes, normal microtubule organization was disrupted, and microtubule stability was significantly compromised. Moreover, GL3 physically bound to the MDP25 promoter, thereby inhibiting its expression. Overexpression of GL3 negated the effects of PMDP25-driven MDP25 or its mutant proteins on trichome branching and microtubules in the mdp25 background. Overall, our study uncovers a mechanism by which GL3 inhibits MDP25 transcription, thereby influencing microtubule stability and regulating trichome branching. This mechanism provides a connection between early regulatory components and microtubules during trichome development.

The receptor-like cytoplasmic kinase OsBSK1-2 regulates immunity via an HLH/bHLH complex

Plants need to fine-tune defense responses to maintain a robust but flexible host barrier to various pathogens. Helix-loop-helix/basic helix-loop-helix (HLH/bHLH) complexes play important roles in fine-tuning plant development. However, the function of these genes in plant immunity and how they are regulated remain obscure. Here, we identified an atypical bHLH transcription factor, Oryza sativa (Os)HLH46, that interacts with rice receptor-like cytoplasmic kinase (RLCK) Os BRASSINOSTEROID-SIGNALING KINASE1-2 (OsBSK1-2), which plays a key role in rice blast resistance. OsBSK1-2 stabilized OsHLH46 both in vivo and in vitro. In addition, OsHLH46 positively regulates rice blast resistance, which depends on OsBSK1-2. OsHLH46 has no transcriptional activation activity and interacts with a typical bHLH protein, OsbHLH6, which negatively regulates rice blast resistance. OsbHLH6 binds to the promoter of OsWRKY45 and inhibits its expression, while OsHLH46 suppresses the function of OsbHLH6 by blocking its DNA binding and transcriptional inhibition of OsWRKY45. Consistent with these findings, OsWRKY45 was up-regulated in OsHLH46-overexpressing plants. In addition, the oshlh46 mutant overexpressing OsbHLH6 is more susceptible to Magnaporthe oryzae than is the wild type, suggesting that OsHLH46 suppresses OsbHLH6-mediated rice blast resistance. Our results not only demonstrated that OsBSK1-2 regulates rice blast resistance via the OsHLH46/OsbHLH6 complex, but also uncovered a new mechanism for plants to fine-tune plant immunity by regulating the HLH/bHLH complex via RLCKs.

Tracing the evolutionary and genetic footprints of atmospheric tillandsioids transition from land to air

Plant evolution is driven by key innovations of functional traits that enables their survivals in diverse ecological environments. However, plant adaptive evolution from land to atmospheric niches remains poorly understood. In this study, we use the epiphytic Tillandsioideae subfamily of Bromeliaceae as model plants to explore their origin, evolution and diversification. We provide a comprehensive phylogenetic tree based on nuclear transcriptomic sequences, indicating that core tillandsioids originated approximately 11.3 million years ago in the Andes. The geological uplift of the Andes drives the divergence of tillandsioids into tank-forming and atmospheric types. Our genomic and transcriptomic analyses reveal gene variations and losses associated with adaptive traits such as impounding tanks and absorptive trichomes. Furthermore, we uncover specific nitrogen-fixing bacterial communities in the phyllosphere of tillandsioids as potential source of nitrogen acquisition. Collectively, our study provides integrative multi-omics insights into the adaptive evolution of tillandsioids in response to elevated aerial habitats.

Membrane Interactions of GET1 and GET2 Facilitate Fiber Cell Initiation through the Guided Entry of the TA Protein Pathway in Cotton

The guided entry of TA proteins (GET) pathway, which is responsible for the post-translational targeting and insertion of the tail-anchored (TA) protein into the endoplasmic reticulum (ER), plays an important role in physiological processes such as protein sorting, vesicle trafficking, cell apoptosis, and enzymatic reactions in which the GET1/2 complex is indispensable. However, a comprehensive study of the GET1 and GET2 genes and the GET pathway in cotton has not yet been carried out. Here, 12 GET1 and 21 GET2 genes were identified in nine representative plant species, and the phylogenetic relationships, gene structures, protein motifs, cis-regulatory elements (CREs), and temporal and spatial expression profiles were analyzed thoroughly. Our study indicated that GhGET1s and GhGET2s might be localized on ER membranes. According to expression profiling and CREs analysis, GhGET2-A02 was identified as a promising candidate for fiber cell development, interacting with two GhGET1s in the membrane, with a binding bias toward GhGET1-A06. Silencing of GhGET1-A06 or GhGET2-A02 reduced fiber initiation and elongation. In summary, our research provides important evidence for understanding the gene families and functions of GET1 and GET2 in cotton and provides clues for molecular breeding of high-quality cotton fiber varieties.

Natural variation at the cotton HIC locus increases trichome density and enhances resistance to aphids

Plant trichomes are an excellent model for studying cell differentiation and development, providing crucial defenses against biotic and abiotic stresses. There is a well-established inverse relationship between trichome density and aphid prevalence, indicating that higher trichome density leads to reduced aphid infestations. Here we present the cloning and characterization of a dominant quantitative trait locus, HIC (hirsute cotton), which significantly enhances cotton trichome density. This enhancement leads to markedly improved resistance against cotton aphids. The HIC encodes an HD-ZIP IV transcriptional activator, crucial for trichome initiation. Overexpression of HIC leads to a substantial increase in trichome density, while knockdown of HIC results in a marked decrease in density, confirming its role in trichome regulation. We identified a variant in the HIC promoter (-810 bp A to C) that increases transcription of HIC and trichome density in hirsute cotton compared with Gossypium hirsutum cultivars with fewer or no trichomes. Interestingly, although the -810 variant in the HIC promoter is the same in G. barbadense and hirsute cotton, the presence of a copia-like retrotransposon insertion in the coding region of HIC in G. barbadense causes premature transcription termination. Further analysis revealed that HIC positively regulates trichome density by directly targeting the EXPANSIN A2 gene, which is involved in cell wall development. Taken together, our results underscore the pivotal function of HIC as a primary regulator during the initial phases of trichome formation, and its prospective utility in enhancing aphid resistance in superior cotton cultivars via selective breeding.

Single-cell landscape of long and short glandular trichomes in Nicotiana tabacum leaves

Glandular trichomes (GTs) play a crucial role in plant defenses and the synthesis of secondary metabolites. Understanding the developmental trajectory of GTs is essential for unraveling their functional significance and potential applications. Here we established a comprehensive single-cell atlas of Nicotiana tabacum leaves, a model plant for GT studies. The atlas included a total of 40,433 cells and successfully captured both long GTs (LGTs) and short GTs (SGTs) from Nicotiana leaves. The developmental trajectories of these trichomes were delineated, revealing potential disparities in epidermal development. Comparative analysis of Arabidopsis and Nicotiana trichome development indicated limited similarity between Arabidopsis epidermal non-glandular trichomes and Nicotiana LGTs and SGTs, implying the essentiality of studying the genes directly involved in the development of Nicotiana GTs for a proper and comprehensive understanding of GT biology. Overall, our results provide profound insights into the developmental intricacies of the specialized GTs.

Genome-wide identification of CaWD40 proteins reveals the involvement of a novel complex (CaAN1-CaDYT1-CaWD40-91) in anthocyanin biosynthesis and genic male sterility in Capsicum annuum

BACKGROUND

The WD40 domain, one of the most abundant in eukaryotic genomes, is widely involved in plant growth and development, secondary metabolic biosynthesis, and mediating responses to biotic and abiotic stresses. WD40 repeat (WD40) protein has been systematically studied in several model plants but has not been reported in the Capsicum annuum (pepper) genome.

RESULTS

Herein, 269, 237, and 257 CaWD40 genes were identified in the Zunla, CM334, and Zhangshugang genomes, respectively. CaWD40 sequences from the Zunla genome were selected for subsequent analysis, including chromosomal localization, phylogenetic relationships, sequence characteristics, motif compositions, and expression profiling. CaWD40 proteins were unevenly distributed on 12 chromosomes, encompassing 19 tandem duplicate gene pairs. The 269 CaWD40s were divided into six main branches (A to F) with 17 different types of domain distribution. The CaWD40 gene family exhibited diverse expression patterns, and several genes were specifically expressed in flowers and seeds. Yeast two-hybrid (Y2H) and dual-luciferase assay indicated that CaWD40-91 could interact with CaAN1 and CaDYT1, suggesting its involvement in anthocyanin biosynthesis and male sterility in pepper.

CONCLUSIONS

In summary, we systematically characterized the phylogeny, classification, structure, and expression of the CaWD40 gene family in pepper. Our findings provide a valuable foundation for further functional investigations on WD40 genes in pepper.

Deciphering the dynamic expression network of fiber elongation and the functional role of the GhTUB5 gene for fiber length in cotton based on an introgression population of upland cotton

INTRODUCTION

Interspecific introgression between Gossypium hirsutum and G. barbadense allows breeding cotton with outstanding fiber length (FL). However, the dynamic gene regulatory network of FL-related genes has not been characterized, and the functional mechanism through which the hub gene GhTUB5 mediates fiber elongation has yet to be determined.

METHODS

Coexpression analyses of 277 developing fiber transcriptomes integrated with QTL mapping using 250 introgression lines of different FL phenotypes were conducted to identify genes related to fiber elongation. The function of GhTUB5 was determined by ectopic expression of two TUB5 alleles in Arabidopsis and knockout of GhTUB5 in upland cotton. Yeast two-hybrid, split-luciferase and pull-down assays were conducted to screen for interacting proteins, and upstream genes were identified by yeast one-hybrid, dual-LUC and electrophoretic mobility shift assays.

RESULTS

The 32,612, 30,837 and 30,277 genes expressed at 5, 10 and 15 days postanthesis (dpa) were grouped into 19 distinct coexpression modules, and 988 genes in the MEblack module were enriched in the cell wall process and exhibited significant associations with FL. A total of 20 FL-QTLs were identified, each explaining 3.34-16.04 % of the phenotypic variance in the FL. Furthermore, several FL-QTLs contained 15 genes that were differentially expressed in the MEblack module including the tubulin beta gene (TUB5). Compared with the wild type, the overexpression of GhTUB5 and GbTUB5 in Arabidopsis suppressed root cell length but promoted cellulose synthesis. Knockout of GhTUB5 resulted in longer fiber lines. Protein-based experiments revealed that GhTUB5 interacts with GhZFP6. Additionally, GhTUB5 was directly activated by GhHD-ZIP7, a homeobox-leucine zipper transcription factor, and its paralogous gene was previously reported to mediate fiber elongation.

CONCLUSION

This study opens a new avenue to dissect functional mechanism of cotton fiber elongation. Our findings provide some molecular details on how GhTUB5 mediates the FL phenotype 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.

A dominant negative mutation of GhMYB25-like alters cotton fiber initiation, reducing lint and fuzz

Cotton (Gossypium hirsutum) fibers, vital natural textile materials, are single-cell trichomes that differentiate from the ovule epidermis. These fibers are categorized as lint (longer fibers useful for spinning) or fuzz (shorter, less useful fibers). Currently, developing cotton varieties with high lint yield but without fuzz remains challenging due to our limited knowledge of the molecular mechanisms underlying fiber initiation. This study presents the identification and characterization of a naturally occurring dominant negative mutation GhMYB25-like_AthapT, which results in a reduced lint and fuzzless phenotype. The GhMYB25-like_AthapT protein exerts its dominant negative effect by suppressing the activity of GhMYB25-like during lint and fuzz initiation. Intriguingly, the negative effect of GhMYB25-like_AthapT could be alleviated by high expression levels of GhMYB25-like. We also uncovered the role of GhMYB25-like in regulating the expression of key genes such as GhPDF2 (PROTODERMAL FACTOR 2), CYCD3; 1 (CYCLIN D3; 1), and PLD (Phospholipase D), establishing its significance as a pivotal transcription factor in fiber initiation. We identified other genes within this regulatory network, expanding our understanding of the determinants of fiber cell fate. These findings offer valuable insights for cotton breeding and contribute to our fundamental understanding of fiber development.

Sphingosine Promotes Fiber Early Elongation in Upland Cotton

Sphingolipids play an important role in cotton fiber development, but the regulatory mechanism is largely unclear. We found that serine palmitoyltransferase (SPT) enzyme inhibitors, myriocin and sphingosine (dihydrosphingosine (DHS) and phytosphingosine (PHS)), affected early fiber elongation in cotton, and we performed a sphingolipidomic and transcriptomic analysis of control and PHS-treated fibers. Myriocin inhibited fiber elongation, while DHS and PHS promoted it in a dose-effect manner. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS), we found that contents of 22 sphingolipids in the PHS-treated fibers for 10 days were changed, of which the contents of 4 sphingolipids increased and 18 sphingolipids decreased. The transcriptome analysis identified 432 differentially expressed genes (238 up-regulated and 194 down-regulated) in the PHS-treated fibers. Among them, the phenylpropanoid biosynthesis pathway is the most significant enrichment. The expression levels of transcription factors such as MYB, ERF, LBD, and bHLH in the fibers also changed, and most of MYB and ERF were up-regulated. Auxin-related genes IAA, GH3 and BIG GRAIN 1 were up-regulated, while ABPs were down-regulated, and the contents of 3 auxin metabolites were decreased. Our results provide important sphingolipid metabolites and regulatory pathways that influence fiber elongation.

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.

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.

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.

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.

Effect of the HCT Gene on Lignin Synthesis and Fiber Development in Gossypium barbadense

As one of the key enzymes in the metabolic pathway of phenylpropane, shikimate hydroxycinnamoyl transferase (HCT) is mainly involved in the biosynthesis of the plant secondary cell wall, which is closely related to cotton fiber quality. In this study, whole-genome identification and bioinformatics analysis of the HCT gene family were performed in G. barbadense. In the whole genome, we identified 136 GbHCT genes encoding 309-504 amino acids. Phylogenetic analysis divided the genome into 5 subfamilies, which were located on 25 chromosomes. Collinear analysis of polyploidization and tandem duplication events were the main driving forces for the rapid expansion and evolution of this family, and the genes underwent loose purifying selection constraints after duplication. Gene promoters identified a variety of cis-acting elements related to plant hormones and the stress response. Several members of the GbHCT family were highly expressed during the development of cotton fiber, and different members had different expression patterns in cotton fiber. After GbHCT114 gene silencing in cotton, the amount of stem surface trichomes and lignin content decreased, and the cell morphology and arrangement changed. After the GbHCT114 gene was overexpressed in Arabidopsis thaliana (L.) Heynh., the number of stem and leaf surface trichomes and the cross-sectional area of the secondary xylem duct cell wall increased. In addition, utilizing transcriptomic analysis, differentially expressed genes associated with lignin synthesis and fiber development were identified. Taken together, the results obtained in this study confirm that the GbHCT114 gene regulates plant trichome development, which lays a theoretical foundation for future studies on the function of GbHCT114 in cotton.

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.

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.

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.

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.