HOXA9 Regulome and Pharmacological Interventions in Leukemia.

HOXA9, an important transcription factor (TF) in hematopoiesis, is aberrantly expressed in numerous cases of both acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL) and is a strong indicator of poor prognosis in patients. HOXA9 is a proto-oncogene which is both sufficient and necessary for leukemia transformation. HOXA9 expression in leukemia correlates with patient survival outcomes and response to therapy. Chromosomal transformations (such as NUP98-HOXA9), mutations, epigenetic dysregulation (e.g., MLL- MENIN -LEDGF complex or DOT1L/KMT4), transcription factors (such as USF1/USF2), and noncoding RNA (such as HOTTIP and HOTAIR) regulate HOXA9 mRNA and protein during leukemia. HOXA9 regulates survival, self-renewal, and progenitor cell cycle through several of its downstream target TFs including LMO2, antiapoptotic BCL2, SOX4, and receptor tyrosine kinase FLT3 and STAT5. This dynamic and multilayered HOXA9 regulome provides new therapeutic opportunities, including inhibitors targeting DOT1L/KMT4, MENIN, NPM1, and ENL proteins. Recent findings also suggest that HOXA9 maintains leukemia by actively repressing myeloid differentiation genes. This chapter summarizes the recent advances understanding biochemical mechanisms underlying HOXA9-mediated leukemogenesis, the clinical significance of its abnormal expression, and pharmacological approaches to treat HOXA9-driven leukemia.

Molecular insights into programmed cell death in esophageal squamous cell carcinoma.

Background: Esophageal squamous cell carcinoma (ESCC) is a deadly type of esophageal cancer. Programmed cell death (PCD) is an important pathway of cellular self-extermination and is closely involved in cancer progression. A detailed study of its mechanism may contribute to ESCC treatment. Methods: We obtained expression profiling data of ESCC patients from public databases and genes related to 12 types of PCD from previous studies. Hub genes in ESCC were screened from PCD-related genes applying differential expression analysis, machine learning analysis, linear support vector machine (SVM), random forest and Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis. In addition, based on the HTFtarget and TargetScan databases, transcription factors (TFs) and miRNAs interacting with the hub genes were selected. The relationship between hub genes and immune cells were analyzed using the CIBERSORT algorithm. Finally, to verify the potential impact of the screened hub genes on ESCC occurrence and development, a series of in vitro cell experiments were conducted. Results: We screened 149 PCD-related DEGs, of which five DEGs (INHBA, LRRK2, HSP90AA1, HSPB8, and EIF2AK2) were identified as the hub genes of ESCC. The area under the curve (AUC) of receiver operating characteristic (ROC) curve of the integrated model developed using the hub genes reached 0.997, showing a noticeably high diagnostic accuracy. The number of TFs and miRNAs regulating hub genes was 105 and 22, respectively. INHBA, HSP90AA1 and EIF2AK2 were overexpressed in cancer tissues and cells of ESCC. Notably, INHBA knockdown suppressed ECSS cell migration and invasion and altered the expression of important apoptotic and survival proteins. Conclusion: This study identified significant molecules with promising accuracy for the diagnosis of ESCC, which may provide a new perspective and experimental basis for ESCC research.

Transcriptional Control: A Directional Sign at the Crossroads of Adult Hepatic Progenitor Cells' Fates.

Hepatic progenitor cells (HPCs) have a bidirectional potential to differentiate into hepatocytes and bile duct epithelial cells and constitute a second barrier to liver regeneration in the adult liver. They are usually located in the Hering duct in the portal vein region where various cells, extracellular matrix, cytokines, and communication signals together constitute the niche of HPCs in homeostasis to maintain cellular plasticity. In various types of liver injury, different cellular signaling streams crosstalk with each other and point to the inducible transcription factor set, including FoxA1/2/3, YB-1, Foxl1, Sox9, HNF4α, HNF1α, and HNF1ß. These transcription factors exert different functions by binding to specific target genes, and their products often interact with each other, with diverse cascades of regulation in different molecular events that are essential for homeostatic regulation, self-renewal, proliferation, and selective differentiation of HPCs. Furthermore, the tumor predisposition of adult HPCs is found to be significantly increased under transcriptional factor dysregulation in transcriptional analysis, and the altered initial commitment of the differentiation pathway of HPCs may be one of the sources of intrahepatic tumors. Related transcription factors such as HNF4α and HNF1 are expected to be future targets for tumor treatment.

Dynamic BTB-domain filaments promote clustering of ZBTB proteins.

The formation of dynamic protein filaments contributes to various biological functions by clustering individual molecules together and enhancing their binding to ligands. We report such a propensity for the BTB domains of certain proteins from the ZBTB family, a large eukaryotic transcription factor family implicated in differentiation and cancer. Working with Xenopus laevis and human proteins, we solved the crystal structures of filaments formed by dimers of the BTB domains of ZBTB8A and ZBTB18 and demonstrated concentration-dependent higher-order assemblies of these dimers in solution. In cells, the BTB-domain filamentation supports clustering of full-length human ZBTB8A and ZBTB18 into dynamic nuclear foci and contributes to the ZBTB18-mediated repression of a reporter gene. The BTB domains of up to 21 human ZBTB family members and two related proteins, NACC1 and NACC2, are predicted to behave in a similar manner. Our results suggest that filamentation is a more common feature of transcription factors than is currently appreciated.

Polymerization of ZBTB transcription factors regulates chromatin occupancy.

BCL6, an oncogenic transcription factor (TF), forms polymers in the presence of a small-molecule molecular glue that stabilizes a complementary interface between homodimers of BCL6's broad-complex, tramtrack, and bric-à-brac (BTB) domain. The BTB domains of other proteins, including a large class of TFs, have similar architectures and symmetries, raising the possibility that additional BTB proteins self-assemble into higher-order structures. Here, we surveyed 189 human BTB proteins with a cellular fluorescent reporter assay and identified 18 ZBTB TFs that show evidence of polymerization. Through biochemical and cryoelectron microscopy (cryo-EM) studies, we demonstrate that these ZBTB TFs polymerize into filaments. We found that BTB-domain-mediated polymerization of ZBTB TFs enhances chromatin occupancy within regions containing homotypic clusters of TF binding sites, leading to repression of target genes. Our results reveal a role of higher-order structures in regulating ZBTB TFs and suggest an underappreciated role for TF polymerization in modulating gene expression.

Genome-wide identification of R2R3-MYB family in Eriobotrya japonica and functional analysis of EjMYB5 involved in proanthocyanidin biosynthesis.

Loquat (Eriobotrya japonica Lindl.) is a popular fruit and medicinal plant. Proanthocyanidins (PAs), as one of the main types of flavonoids, are the key components of loquat fruit quality and medicinal properties. However, the identification of transcription factors (TFs) involved in PA accumulation in loquat remains limited. R2R3-MYB TFs play key regulatory role in PA accumulation in plants. In this study, 190 R2R3-MYB TFs were identified in loquat genome. Combined with transcriptome data, R2R3-MYB TF EjMYB5 involved in PA accumulation in loquat was isolated. EjMYB5 was transcriptional activator localized to nucleus. Expression of EjMYB5 was closely related to PA accumulation in loquat fruits. Heterogenous overexpression of EjMYB5 in tomato (Solanum lycopersicum) inhibited anthocyanin accumulation and promoted PA accumulation. Additionally, transient overexpression of EjMYB5 in tobacco (Nicotiana benthamiana) leaves promoted PA accumulation by upregulating flavonoid biosynthesis genes (NtDFR, NtANS, and NtLAR). Transcriptome analysis of EjMYB5-overexpressing tomato fruits suggested that EjMYB5 was involved in several biological pathways, including lipid metabolism, MAPK signaling, phenylpropanoid biosynthesis, and flavonoid biosynthesis. Collectively, our findings provided basic data for further analysis the function of R2R3-MYB TFs in loquat, and revealed that EjMYB5 functioned as PA accumulation in loquat.

Eggplant NAC domain transcription factor SmNST1 as an activator promotes secondary cell wall thickening.

NAC-domain transcription factors (TFs) are plant-specific transcriptional regulators playing crucial roles in plant secondary cell wall (SCW) biosynthesis. SCW is important for plant growth and development, maintaining plant morphology, providing rigid support, ensuring material transportation and participating in plant stress responses as a protective barrier. However, the molecular mechanisms underlying SCW in eggplant have not been thoroughly explored. In this study, the NAC domain TFs SmNST1 and SmNST2 were cloned from the eggplant line 'Sanyue qie'. SmNST1 and SmNST2 expression levels were the highest in the roots and stems. Subcellular localization analysis showed that they were localized in the cell membrane and nucleus. Their overexpression in transgenic tobacco showed that SmNST1 promotes SCW thickening. The expression of a set of SCW biosynthetic genes for cellulose, xylan and lignin, which regulate SCW formation, was increased in transgenic tobacco. Bimolecular fluorescence and luciferase complementation assays showed that SmNST1 interacted with SmNST2 in vivo. Yeast one-hybrid, electrophoretic mobility shift assay (EMSA) and Dual-luciferase reporter assays showed that SmMYB26 directly bound to the SmNST1 promoter and acted as an activator. SmNST1 and SmNST2 interact with the SmMYB108 promoter and repress SmMYB108 expression. Altogether, we showed that SmNST1 positively regulates SCW formation, improving our understanding of SCW biosynthesis transcriptional regulation.

Early B-Cell Factor 1: An Archetype for a Lineage-Restricted Transcription Factor Linking Development to Disease.

The development of highly specialized blood cells from hematopoietic stem cells (HSCs) in the bone marrow (BM) is dependent upon a stringently orchestrated network of stage- and lineage-restricted transcription factors (TFs). Thus, the same stem cell can give rise to various types of differentiated blood cells. One of the key regulators of B-lymphocyte development is early B-cell factor 1 (EBF1). This TF belongs to a small, but evolutionary conserved, family of proteins that harbor a Zn-coordinating motif and an IPT/TIG (immunoglobulin-like, plexins, transcription factors/transcription factor immunoglobulin) domain, creating a unique DNA-binding domain (DBD). EBF proteins play critical roles in diverse developmental processes, including body segmentation in the Drosophila melanogaster embryo, and retina formation in mice. While several EBF family members are expressed in neuronal cells, adipocytes, and BM stroma cells, only B-lymphoid cells express EBF1. In the absence of EBF1, hematopoietic progenitor cells (HPCs) fail to activate the B-lineage program. This has been attributed to the ability of EBF1 to act as a pioneering factor with the ability to remodel chromatin, thereby creating a B-lymphoid-specific epigenetic landscape. Conditional inactivation of the Ebf1 gene in B-lineage cells has revealed additional functions of this protein in relation to the control of proliferation and apoptosis. This may explain why EBF1 is frequently targeted by mutations in human leukemia cases. This chapter provides an overview of the biochemical and functional properties of the EBF family proteins, with a focus on the roles of EBF1 in normal and malignant B-lymphocyte development.

Forkhead box O1 in metabolic dysfunction-associated fatty liver disease: molecular mechanisms and drug research.

Metabolic dysfunction-associated fatty liver disease (MAFLD) is a chronic liver disease that progresses from hepatic steatosis to non-alcoholic steatohepatitis, cirrhosis, and liver cancer, posing a huge burden on human health. Existing research has confirmed that forkhead box O1 (FOXO1), as a member of the FOXO transcription factor family, is upregulated in MAFLD. Its activity is closely related to nuclear-cytoplasmic shuttling and various post-translational modifications including phosphorylation, acetylation, and methylation. FOXO1 mediates the progression of MAFLD by regulating glucose metabolism, lipid metabolism, insulin resistance, oxidative stress, hepatic fibrosis, hepatocyte autophagy, apoptosis, and immune inflammation. This article elaborates on the regulatory role of FOXO1 in MAFLD, providing a summary and new insights for the current status of drug research and targeted therapies for MAFLD.

Chromatin accessibility and transcriptional landscape in PK-15 cells during early exposure to Aflatoxin B1.

Aflatoxin B1 (AFB1) not only causes significant losses in livestock production but also poses a serious threat to human health. It is the most carcinogenic among known chemicals. Pigs are more susceptible to AFB1 and experience a higher incidence. However, the molecular mechanism of the toxic effect of AFB1 remains unclear. In this study, we used assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA-seq to uncover chromatin accessibility and gene expression dynamics in PK-15 cells during early exposure to AFB1. We observed that the toxic effects of AFB1 involve signaling pathways such as p53, PI3K-AKT, Hippo, MAPK, TLRs, apoptosis, autophagy, and cancer pathways. Basic leucine zipper (bZIP) transcription factors (TFs), including AP-1, Fos, JunB, and Fra2, play a crucial role in regulating the biological processes involved in AFB1 challenge. Several new TFs, such as BORIS, HNF1b, Atf1, and KNRNPH2, represent potential targets for the toxic mechanism of AFB1. In addition, it is crucial to focus on the concentration of intracellular zinc ions. These findings will contribute to a better understanding of the mechanisms underlying AFB1-induced nephrotoxicity and offer new molecular targets.

RUNX1 regulates MCM2/CDC20 to promote COAD progression modified by deubiquitination of USP31.

Colon adenocarcinoma (COAD) is the second leading cause of cancer death, and there is still a lack of diagnostic biomarkers and therapeutic targets. In this study, bioinformatics analysis of the TCGA database was used to obtain RUNX1, a gene with prognostic value in COAD. RUNX1 plays an important role in many malignancies, and its molecular regulatory mechanisms in COAD remain to be fully understood. To explore the physiological role of RUNX1, we performed functional analyses, such as CCK-8, colony formation and migration assays. In addition, we investigated the underlying mechanisms using transcriptome sequencing and chromatin immunoprecipitation assays. RUNX1 is highly expressed in COAD patients and significantly correlates with survival. Silencing of RUNX1 significantly slowed down the proliferation and migratory capacity of COAD cells. Furthermore, we demonstrate that CDC20 and MCM2 may be target genes of RUNX1, and that RUNX1 may be physically linked to the deubiquitinating enzyme USP31, which mediates the upregulation of RUNX1 protein to promote transcriptional function. Our results may provide new insights into the mechanism of action of RUNX1 in COAD and reveal potential therapeutic targets for this disease.

Durable Objective Response to Lurbinectedin in Small Cell Bladder Cancer with TP53 Mutation: A Molecular-Directed Strategy.

Small cell bladder cancer (SCBC) is a rare and aggressive disease, often treated with platinum/etoposide-based chemotherapy. Key molecular drivers include the inactivation of onco-suppressor genes (TP53, RB1) and amplifications in proto-oncogenes (MYC). We report a patient with SCBC who achieved an objective and prolonged response to lurbinectedin, which has been approved for metastatic small cell lung cancer, after developing disease progression on cisplatin/etoposide and nivolumab/ipilimumab. A genomic analysis of a metastatic biopsy prior to lurbinectedin initiation revealed a TP53 mutation and amplification of the cell cycle regulators E2F3 and MYCL. A repeat biopsy following the development of lurbinectedin resistance showed a new actionable ERBB2 alteration without significant change in the tumor mutation burden (six mutations/Mb). The present report suggests that lurbinectedin may be active and should be further explored in SCBC harboring TP53 mutations and amplifications in E2F3 and MYC family complexes.

Emerging role of MYB transcription factors in cancer drug resistance.

Decades ago, the viral myeloblastosis oncogene v-myb was identified as a gene responsible for the development of avian leukemia. However, the relevance of MYB proteins for human cancer diseases, in particular for solid tumors, remained basically unrecognized for a very long time. The human family of MYB transcription factors comprises MYB (c-MYB), MYBL2 (b-MYB), and MYBL1 (a-MYB), which are overexpressed in several cancers and are associated with cancer progression and resistance to anticancer drugs. In addition to overexpression, the presence of activated MYB-fusion proteins as tumor drivers was described in certain cancers. The identification of anticancer drug resistance mediated by MYB proteins and their underlying mechanisms are of great importance in understanding failures of current therapies and establishing new and more efficient therapy regimens. In addition, new drug candidates targeting MYB transcription factor activity and signaling have emerged as a promising class of potential anticancer therapeutics that could tackle MYB-dependent drug-resistant cancers in a more selective way. This review describes the correlation of MYB transcription factors with the formation and persistence of cancer resistance to various approved and investigational anticancer drugs.

Genome-Wide Analysis of BBX Gene Family in Three Medicago Species Provides Insights into Expression Patterns under Hormonal and Salt Stresses.

BBX protein is a class of zinc finger transcription factors that have B-box domains at the N-terminus, and some of these proteins contain a CCT domain at the C-terminus. It plays an important role in plant growth, development, and metabolism. However, the expression pattern of BBX genes in alfalfa under hormonal and salt stresses is still unclear. In this study, we identified a total of 125 BBX gene family members by the available Medicago reference genome in diploid alfalfa (Medicago sativa spp. Caerulea), a model plant (M. truncatula), and tetraploid alfalfa (M. sativa), and divided these members into five subfamilies. We found that the conserved motifs of BBXs of the same subfamily reveal similarities. We analyzed the collinearity relationship and duplication mode of these BBX genes and found that the expression pattern of BBX genes is specific in different tissues. Analysis of the available transcriptome data suggests that some members of the BBX gene family are involved in multiple abiotic stress responses, and the highly expressed genes are often clustered together. Furthermore, we identified different expression patterns of some BBX genes under salt, ethylene, salt and ethylene, salicylic acid, and salt and salicylic acid treatments, verified by qRT-PCR, and analyzed the subcellular localization of MsBBX2, MsBBX17, and MsBBX32 using transient expression in tobacco. The results showed that BBX genes were localized in the nucleus. This study systematically analyzed the BBX gene family in Medicago plants, which provides a basis for the study of BBX gene family tolerance to abiotic stresses.

Pooled CRISPR Interference Screening Identifies Crucial Transcription Factors in Gas-Fermenting Clostridium ljungdahlii.

Gas-fermenting Clostridium species hold tremendous promise for one-carbon biomanufacturing. To unlock their full potential, it is crucial to unravel and optimize the intricate regulatory networks that govern these organisms; however, this aspect is currently underexplored. In this study, we employed pooled CRISPR interference (CRISPRi) screening to uncover a wide range of functional transcription factors (TFs) in Clostridium ljungdahlii, a representative species of gas-fermenting Clostridium, with a special focus on TFs associated with the utilization of carbon resources. Among the 425 TF candidates, we identified 75 and 68 TF genes affecting the heterotrophic and autotrophic growth of C. ljungdahlii, respectively. We focused our attention on two of the screened TFs, NrdR and DeoR, and revealed their pivotal roles in the regulation of deoxyribonucleoside triphosphates (dNTPs) supply, carbon fixation, and product synthesis in C. ljungdahlii, thereby influencing the strain performance in gas fermentation. Based on this, we proceeded to optimize the expression of deoR in C. ljungdahlii by adjusting its promoter strength, leading to an improved growth rate and ethanol synthesis of C. ljungdahlii when utilizing syngas. This study highlights the effectiveness of pooled CRISPRi screening in gas-fermenting Clostridium species, expanding the horizons for functional genomic research in these industrially important bacteria.

Unveiling the potential of FOXO3 in lung cancer: From molecular insights to therapeutic prospects.

Lung cancer poses a significant challenge regarding molecular heterogeneity, as it encompasses a wide range of molecular alterations and cancer-related pathways. Recent discoveries made it feasible to thoroughly investigate the molecular mechanisms underlying lung cancer, giving rise to the possibility of novel therapeutic strategies relying on molecularly targeted drugs. In this context, forkhead box O3 (FOXO3), a member of forkhead transcription factors, has emerged as a crucial protein commonly dysregulated in cancer cells. The regulation of the FOXO3 in reacting to external stimuli plays a key role in maintaining cellular homeostasis as a component of the molecular machinery that determines whether cells will survive or dies. Indeed, various extrinsic cues regulate FOXO3, affecting its subcellular location and transcriptional activity. These regulations are mediated by diverse signaling pathways, non-coding RNAs (ncRNAs), and protein interactions that eventually drive post-transcriptional modification of FOXO3. Nevertheless, while it is no doubt that FOXO3 is implicated in numerous aspects of lung cancer, it is unclear whether they act as tumor suppressors, promotors, or both based on the situation. However, FOXO3 serves as an intriguing possible target in lung cancer therapeutics while widely used anti-cancer chemo drugs can regulate it. In this review, we describe a summary of recent findings on molecular mechanisms of FOXO3 to clarify that targeting its activity might hold promise in lung cancer treatment.

Multiomic analysis of monocyte-derived alveolar macrophages in idiopathic pulmonary fibrosis.

BACKGROUND: Monocyte-derived alveolar macrophages (Mo_AMs) are increasingly recognised as potential pathogenic factors for idiopathic pulmonary fibrosis (IPF). While scRNAseq analysis has proven valuable in the transcriptome profiling of Mo_AMs, the integration analysis of multi-omics may provide additional dimensions of understanding of these cellular populations. METHODS: We performed multi-omics analysis on 116 scRNAseq, 119 bulkseq and five scATACseq lung tissue samples from IPF. We built a large-scale IPF scRNAseq atlas and conducted the Monocle 2/3 as well as the Cellchat to explore the developmental path and intercellular communication on Mo_AMs. We also reported the difference in metabolisms, tissue repair and phagocytosis between Mo_AMs and tissue-resident alveolar macrophages (TRMs). To determine whether Mo_AMs affected pulmonary function, we projected clinical phenotypes (FVC%pred) from the bulkseq dataset onto the scRNAseq atlas. Finally, we used scATATCseq to uncover the upstream regulatory mechanisms and determine key drivers in Mo_AMs. RESULTS: We identified three Mo_AMs clusters and the trajectory analysis further validated the origin of these clusters. Moreover, via the Cellchat analysis, the CXCL12/CXCR4 axis was found to be involved in the molecular basis of reciprocal interactions between Mo_AMs and fibroblasts through the activation of the ERK pathway in Mo_AMs. SPP1_RecMacs (RecMacs, recruited macrophages) were higher in the low-FVC group than in the high-FVC group. Specifically, compared with TRMs, the functions of lipid and energetic metabolism as well as tissue repair were higher in Mo_AMs than TRMs. But, TRMs may have higher level of phagocytosis than TRMs. SPIB (PU.1), JUNB, JUND, BACH2, FOSL2, and SMARCC1 showed stronger association with open chromatin of Mo_AMs than TRMs. Significant upregulated expression and deep chromatin accessibility of APOE were observed in both SPP1_RecMacs and TRMs. CONCLUSION: Through trajectory analysis, it was confirmed that SPP1_RecMacs derived from Monocytes. Besides, Mo_AMs may influence FVC% pred and aggravate pulmonary fibrosis through the communication with fibroblasts. Furthermore, distinctive transcriptional regulators between Mo_AMs and TRMs implied that they may depend on different upstream regulatory mechanisms. Overall, this work provides a global overview of how Mo_AMs govern IPF and also helps determine better approaches and intervention therapies.

Epigenomic signatures of sarcomatoid differentiation to guide the treatment of renal cell carcinoma.

Renal cell carcinoma with sarcomatoid differentiation (sRCC) is associated with poor survival and a heightened response to immune checkpoint inhibitors (ICIs). Two major barriers to improving outcomes for sRCC are the limited understanding of its gene regulatory programs and the low diagnostic yield of tumor biopsies due to spatial heterogeneity. Herein, we characterized the epigenomic landscape of sRCC by profiling 107 epigenomic libraries from tissue and plasma samples from 50 patients with RCC and healthy volunteers. By profiling histone modifications and DNA methylation, we identified highly recurrent epigenomic reprogramming enriched in sRCC. Furthermore, CRISPRa experiments implicated the transcription factor FOSL1 in activating sRCC-associated gene regulatory programs, and FOSL1 expression was associated with the response to ICIs in RCC in two randomized clinical trials. Finally, we established a blood-based diagnostic approach using detectable sRCC epigenomic signatures in patient plasma, providing a framework for discovering epigenomic correlates of tumor histology via liquid biopsy.

Comparative transcriptome analysis reveals genes involved in trichome development and metabolism in tobacco.

BACKGROUND: The glandular trichomes of tobacco (Nicotiana tabacum) can efficiently produce secondary metabolites. They act as natural bioreactors, and their natural products function to protect plants against insect-pests and pathogens and are also components of industrial chemicals. To clarify the molecular mechanisms of tobacco glandular trichome development and secondary metabolic regulation, glandular trichomes and glandless trichomes, as well as other different developmental tissues, were used for RNA sequencing and analysis. RESULTS: By comparing glandless and glandular trichomes with other tissues, we obtained differentially expressed genes. They were obviously enriched in KEGG pathways, such as cutin, suberine, and wax biosynthesis, flavonoid and isoflavonoid biosynthesis, terpenoid biosynthesis, and plant-pathogen interaction. In particular, the expression levels of genes related to the terpenoid, flavonoid, and wax biosynthesis pathway mainly showed down-regulation in glandless trichomes, implying that they lack the capability to synthesize certain exudate compounds. Among the differentially expressed genes, 234 transcription factors were found, including AP2-ERFs, MYBs, bHLHs, WRKYs, Homeoboxes (HD-ZIP), and C2H2-ZFs. These transcription factor and genes that highly expressed in trichomes or specially expressed in GT or GLT. Following the overexpression of R2R3-MYB transcription factor Nitab4.5_0011760g0030.1 in tobacco, an increase in the number of branched glandular trichomes was observed. CONCLUSIONS: Our data provide comprehensive gene expression information at the transcriptional level and an understanding of the regulatory pathways involved in glandular trichome development and secondary metabolism.