BACKGROUND: The camphor tree (Cinnamomum camphora) is a valuable source of essential oils rich in monoterpenes and sesquiterpenes, which exhibit antibacterial, antioxidant, and insecticidal properties. The AP2/ERF transcription factor family is known for its roles in regulating secondary metabolite biosynthesis and stress responses in plants. However, the functions and regulatory mechanisms of AP2/ERF transcription factors in camphor tree remain largely unexplored. RESULTS: A total of 154 unique AP2/ERF genes were identified in the camphor tree genome and classified into four subfamilies (RAV, AP2, ERF, and DREB), among which the ERF subfamily comprises 75 members and represents the largest. The physicochemical properties of the proteins were diverse, with predicted nuclear localization. Phylogenetic analysis revealed distinct clustering patterns consistent with their subfamily classification. Gene structure analysis showed variation in exon-intron organization among subfamilies. A total of 5,910 cis-acting elements were identified in the promoters, with MeJA-responsive and stress responsiveness elements being most prevalent. Expression profile analysis indicated tissue-specific expression patterns, with most CcAP2/ERF genes highly expressed in the roots. Correlation analysis revealed that 44 CcAP2/ERF genes were strongly correlated with 6 terpene synthase (TPS) genes involved in monoterpene biosynthesis, suggesting regulatory functions. Strikingly, CcERF_104 was identified as a high-potential hub regulator, co-expressed two with borneol diphosphate synthase (BPPS) genes (CcTPS26 and CcTPS49). Moreover, a suite of regulators (including CcERF_27, CcERF_42, CcERF_83, CcERF_87, CcAP2_12, and CcAP2_13) showed strong positive correlation with CcTPS72, the most catalytically efficient BPPS synthase in C. camphora, highlighting sophisticated transcriptional control over monoterpene biosynthesis. CONCLUSION: This study provides a comprehensive analysis of the AP2/ERF transcription factor family in C. camphora and elucidates their potential roles in regulating terpene biosynthesis. The findings lay the foundation for future functional characterization of key CcAP2/ERF genes using genetic and biochemical approaches, with the goal of enhancing terpenoid content in C. camphora leaves through genetic improvement technologies.

Although Cinnamomum camphora's metabolic composition has been well studied, flavonoid distribution across its tissues remains poorly understood. This study combined transcriptome and metabolomic analyses on leaf, stem, and root tissues to uncover the synthesis pathway of flavonoids and to identify key regulatory genes. Metabolomic analysis revealed 2,893 metabolites, which can be divided into secondary metabolite 1,213(41.93

Background: Breast cancer treatment has been a global puzzle, and targeted strategies based on the hypoxic tumor microenvironment (TME) have attracted extensive attention. As a signature transcription factor overexpressed in hypoxia tumor, hypoxia-inducible factor-1 (HIF-1) contribute to cancer progression. Compound 7-(3-(2-chloro-1H-benzo[d]1midazole-1-yl) propoxy)-2-(3,4,5-trimethoxyphenyl)-4H-chromen-4-one, synthesized and named FB15 in our earlier research, a potential inhibitor of HIF-1 alpha signaling pathway, has been proved a promising drug candidate for many kinds of cancer chemotherapy. However, the poor solubility and undesirable pharmacokinetics of FB15 leads to limited treatment efficacy of tumor, which ultimately restricts its potential clinical applications. Carbonic anhydrase IX (CAIX), a tumor cell transmembrane protein, was overexpressed in hypoxia tumor site. Acetazolamide (AZA), a highly selective ligand targeting CAIX, can be utilized to delivery FB15 to hypoxia tumor site. Methods: In this study, we prepared and characterized FB15 loaded nano-mixed micelles with the AZA conjugated poloxamer 188 (AZA-P188) and D-a-Tocopherol Polyethylene 1000 Glycol Succinate (TPGS), denoted as, AZA-P188/TPGS@FB15. Its delivery efficiency in vitro and in vivo was assessed by in vitro drug release, cytotoxicity assay, cellular uptake, and in vivo pharmacokinetics and fluorescence imaging. Finally, therapeutic effect of AZA-P188/TPGS@FB15 was investigated using a preclinical breast cancer subcutaneous graft model in vivo. Results: In vitro studies revealed that AZA-P188/TPGS@FB15 could efficiently target breast cancer cells mediated by CAIX receptor, trigger FB15 release in response to acidic condition, and enhance cellular uptake and cytotoxicity against breast cancer cells. The pharmacokinetic studies showed that FB15-loaded AZA-functionalized micelles exhibited significantly increased AUC0-t over free FB15. In vivo imaging demonstrated that AZA-functionalized micelles significantly increased the drug distribution in the tumor site. In vivo experiments confirmed that AZA-P188/TPGS@FB15 exhibited superior inhibition of tumor growth in nude mice with good biosafety. Conclusion: AZA-P188/TPGS@FB15 hold promise as a potentially effective therapeutic way for breast cancer. Its targeted delivery system utilizing AZA as a carrier shows potential for improving the efficacy of FB15 in cancer therapy.

Abstract Background 6-Methoxydihydrosanguinarine (6-MDS) has shown promising potential in fighting against a variety of malignancies. Yet, its anti‑lung adenocarcinoma (LUAD) effect and the underlying mechanism remain largely unexplored. This study sought to explore the targets and the probable mechanism of 6-MDS in LUAD through network pharmacology and experimental validation. Methods The proliferative activity of human LUAD cell line A549 was evaluated by Cell Counting Kit-8 (CCK8) assay. LUAD related targets, potential targets of 6-MDS were obtained from databases. Venn plot analysis were performed on 6-MDS target genes and LUAD related genes to obtain potential target genes for 6-MDS treatment of LUAD. The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database was utilized to perform a protein-protein interaction (PPI) analysis, which was then visualized by Cytoscape. The hub genes in the network were singled out by CytoHubba. Metascape was employed for GO and KEGG enrichment analyses. molecular docking was carried out using AutoDock Vina 4.2 software. Gene expression levels, overall survival of hub genes were validated by the GEPIA database. Protein expression levels, promotor methylation levels of hub genes were confirmed by the UALCAN database. Timer database was used for evaluating the association between the expression of hub genes and the abundance of infiltrating immune cells. Furthermore, correlation analysis of hub genes expression with immune subtypes of LUAD were performed by using the TISIDB database. Finally, the results of network pharmacology analysis were validated by qPCR. Results Experiments in vitro revealed that 6-MDS significantly reduced tumor growth. A total of 33 potential targets of 6-MDS in LUAD were obtained by crossing the LUAD related targets with 6-MDS targets. Utilizing CytoHubba, a network analysis tool, the top 10 genes with the highest centrality measures were pinpointed, including MMP9, CDK1, TYMS, CCNA2, ERBB2, CHEK1, KIF11, AURKB, PLK1 and TTK. Analysis of KEGG enrichment hinted that these 10 hub genes were located in the cell cycle signaling pathway, suggesting that 6-MDS may mainly inhibit the occurrence of LUAD by affecting the cell cycle. Molecular docking analysis revealed that the binding energies between 6-MDS and the hub proteins were all higher than − 6 kcal/Mol with the exception of AURKB, indicating that the 9 targets had strong binding ability with 6-MDS.These results were corroborated through assessments of mRNA expression levels, protein expression levels, overall survival analysis, promotor methylation level, immune subtypes andimmune infiltration. Furthermore, qPCR results indicated that 6-MDS can significantly decreased the mRNA levels of CDK1, CHEK1, KIF11, PLK1 and TTK. Conclusions According to our findings, it appears that 6-MDS could possibly serve as a promising option for the treatment of LUAD. Further investigations in live animal models are necessary to confirm its potential in fighting cancer and to delve into the mechanisms at play.

Cinnamomum burmannii, a medicinal plant with a history of traditional use in Chinese folk medicine, utilizes various plant parts, including bark, root bark, leaves, and stems, for therapeutic purposes. Recent investigations into the essential oil extracted from the twigs and leaves of C. burmannii have revealed a notable abundance of volatile monoterpenes, especially D-borneol, α-pinene, and camphene. In this study, an extensive chemical profiling on the essential oil of the roots, stems, and leaves of C. burmannii was conducted. The analysis results indicated that the root and leaf components exhibit the most diverse and the most abundant of volatile monoterpenes, respectively. To elucidate the biosynthesis of monoterpenes in C. burmannii, candidate genes with monoterpene synthase activities were identified through transcriptome sequencing. Subsequently, function characterization on three mono-terpene synthases (TPSs), designated as CbTPS1, CbTPS2, and CbTPS3, were conducted using phylogenetic analysis and heterogeneous expression in Escherichia coli. The primary enzymatic products were identified as 3-carene, α-phellandrene, and bornyl diphosphate (BPP), respectively. Additionally, a confocal laser microscopy assay suggested the chloroplast localization of these mono-TPSs through transient expression in tobacco. Further validation of their functionality was confirmed through eukaryotic expression in tobacco. In conclusion, this study has unveiled critical enzymes responsible for the biosynthesis of major monoterpenes in C. burmannii. These findings provide essential elements for future studies in synthetic biology, facilitating a deep understanding of the biosynthetic pathways and potential applications in medicinal plant engineering.

Benign prostatic hyperplasia (BPH) as a common geriatric disease in urology, the incidence and prevalence are rapidly increasing with the aging society, prompting an urgent need for effective prevention and treatment of BPH. However, limited therapeutic efficacy and higher risk of complications result in the treatment of BPH remaining challenging. The unclear pathogenic mechanism also hampers further exploration of therapeutic approaches for BPH. In this study, we used multi-omics methods to integrate genomics, transcriptomics, immunomics, and metabolomics data and identify biomolecules associated with BPH. We performed transcriptomic imputation, summary data-based Mendelian randomization (SMR), joint/conditional analysis, colocalization analysis, and FOCUS to explore high-confidence genes associated with BPH in blood and prostate tissue. Subsequently, three-step SMR was used to identify the DNA methylation sites regulating high-confidence genes to improve the pathogenic pathways of BPH. We also used cis-instruments of druggable genes to conduct SMR analysis to find potential drug targets for BPH. Finally, we used MR analysis to explore the immune pathways and metabolomics related to BPH. Multiple analytical methods identified BTN3A2 (Blood: TWAS Z score = 5.02912, TWAS P = 4.93 × 10-7; Prostate: TWAS Z score = 4.89, TWAS P = 1.01 × 10-6) and C4A (Blood: TWAS Z score = 4.90754, TWAS P = 9.22 × 10-7; Prostate: TWAS Z score = 5.084, TWAS P = 3.70 × 10-7) as high-confidence genes for BPH and identified the cg14345882-BTN3A2-BPH pathogenic pathway. We also used druggable gene data to identify 30 promising therapeutic target genes, including BTN3A2 and C4A. For MR analysis of immune pathways, we identified immune cell surface molecules as well as the inflammatory factor IL-17 (OR = 1.25, 95% CI = 1.09-1.43, PFDR = 0.12, Maximum likelihood) as risk factors for BPH. In addition, we found that disulfide levels of cysteinylglycine (OR = 1.11, 95% CI = 1.05-1.18, P = 5.18 × 10-4, Weighted median), oxidation levels of cysteinylglycine (OR = 1.09, 95% CI = 1.04-1.14, P = 3.87 × 10-4, Weighted median), and sebacate levels (OR = 1.05, 95% CI = 1.02-1.08, P = 3.0 × 10-4, Maximum likelihood) increase the risk of BPH. This multi-omics study explored biomolecules associated with BPH, improved the pathogenic pathways of BPH, and identified promising therapeutic targets. Our results provide evidence for future studies aimed at developing appropriate therapeutic interventions.

Purpose: Gestational diabetes mellitus (GDM) is a common metabolic disorder during pregnancy that is associated with placental inflammation and adverse pregnancy outcomes. However, the mechanisms of inflammation in GDM are still unclear. Methods: Bulk transcriptome, single-cell transcriptome, clinical information, and samples were collected from GSE154414, GSE70493, GSE173193 and a retrospective cohort. Bioinformatics prediction was used to explore the mechanisms of placental inflammation, and multiplex immunofluorescence was used to validate the results. Results: First, we found that GDM is characterized by low-grade inflammation and is linked to several adverse pregnancy outcomes, as supported by our collected clinical data. Additionally, we identified ten hub genes (FCGR3B, CXCR1, MMP9, ITGAX, CCL5, GZMB, S100A8, LCN2, TGFB1, and LTF) as potential therapy targets and confirmed the binding of corresponding predictive therapeutic agents by molecular docking. Transcriptome sequencing analysis has shown that macrophages are primarily responsible for the emergence of placental inflammation, and that M1 macrophage polarization increased while M2 macrophage polarization decreased in GDM when compared to the control sample. Multiplex immunofluorescence staining of CD68, CD80, and ACSL4 was performed and suggested that ferroptosis of macrophages may contribute to placental inflammation in GDM. Conclusion: In conclusion, our findings provide a better understanding of the mechanisms of inflammation in GDM and suggest for this condition.

Background: Presenilin enhancer-2 (PSENEN, PEN-2), one of the four components of the γ-secretase complex, has been increasingly revealed to be significant in cancer types such as pancreatic cancer, gastric cancer, and breast cancer. However, the pan-cancer clinical relevance of PSENEN remains unclear. Methods: Raw data on PSENEN expression in normal and cancer tissues were obtained from The Cancer Genome Atlas (TGCA) database and Genotype-Tissue Expression Project (GTEx). The difference of PSENEN expression was analyzed using data from the TCGA repository and TIMER2 database. Meanwhile, Cox regression analysis and KM plotter were used to analyze the pan-cancer prognostic significance of PSENEN. We also analyzed the correlation between PSENEN expression and tumor immune infiltration using the TIMER and XCELL algorithms. Moreover, we used the TISIDB database to determine PSENEN expression in different immune and molecular subtypes of human cancers. Pan-cancer analysis of genetic alteration of PSENEN was performed using online tools such as cBioPortal and UALCAN. Finally, six Gene Expression Omnibus datasets from the Gene Expression Profiling Interactive Analysis database were used to validate the expression level of PSENEN in lung adenocarcinoma (LUAD). Results: Contrary to nonmalignant tissues, the expression level of PSENEN was significantly upregulated in the cancerous tissues in 22 cancer types. Elevated PSENEN expression was correlated with worse overall survival in 7 cancer types and with worse disease-specific survival in 8 of them. By using xCell algorithm, TIMER algorithm, and Spearman's correlation analysis, we found that PSENEN expression was closely correlated with the infiltration of immune cells across all cancer types. Moreover, aberrant PSENEN expression was associated with 60 immune checkpoint pathway-related genes, microsatellite instability, and tumor mutation burden in various cancer types. Besides, PSENEN expression was significantly associated with different immune subtypes and molecular subtypes of various human cancers. Notably, ovarian epithelial tumor samples demonstrated the highest frequency of PSENEN mutation among all cancer types. Conclusions: Our pan-cancer analysis demonstrated that PSENEN might serve as a prognostic biomarker and revealed the importance of an in-depth understanding of the oncogenic role of PSENEN in various malignancies.

Blumea balsamifera, a wooden plant belonging to the family Asteraceae, is a medicinal herb with anticancer, antiviral, and multiple pharmacological effects, which are believed to be caused by its essential oil. The essential oil from B. balsamifera is comprised of mono- and sesqui-terpenes as the majority. Unfortunately, this plant has been facing the challenge of resource shortage, which could be effectively alleviated by biological engineering. Therefore, the identification of key elements involved in the biosynthesis of active ingredients becomes an indispensable prerequisite. In this study, candidate genes encoding monoterpene synthase were screened by transcriptome sequencing combined with metabolomics profiling in the roots, stems, and leaves of B. balsamifera. Then, these candidates were successfully cloned and verified by heterologous expression and in vitro enzyme activity assays. As a result, six candidate BbTPS genes were isolated from B. balsamifera, of which three encoded single-product monoterpene synthases and one encoded a multi-product monoterpene synthase. Among them, BbTPS1, BbTPS3, and BbTPS4 could catalyze the formation of D-limonene, α-phellandrene, and L-borneol, respectively. Meanwhile, BbTPS5 functioned in catalyzing GPP into terpinol, β-phellandrene, β-myrcene, D-limonene, and 2-carene in vitro. In general, our results provided important elements for the synthetic biology of volatile terpenes in B. balsamifera, which laid a foundation for subsequent heterologous production of these terpenoids through metabolic engineering and increasing their yield, as well as promoting sustainable development and utilization of B. balsamifera.

Terpene synthase (TPS) plays an important role in terpenoids biosynthesis. Cinnamomum camphora (camphor tree) contains dozens of terpenoids with medicinal value, especially borneol, which has been widely used since ancient times. However, limited information is available regarding the genome-wide identification and characterization of the TPS family in the C. camphora. In this study, 82 CcTPS genes were identified from the camphor tree genome (CTG). Gene cluster and sequence syntenic analysis suggested that tandem duplication occurred within the TPS family of the CTG, especially for the TPS-b subfamily. The chemotype-specific gene expression analysis showed significantly differential expression patterns among six chemotypes. It is worth noting that three genes (CcTPS26, CcTPS49 and CcTPS72) exhibited relatively high expression in the borneol-type camphor tree, compared to the other five chemotypes. Further functional characterization of them indicated that they were all bornyl diphosphate synthases (BPPSs), which function in catalyzing GPP into BPP and then undergoes dephosphorylation to yield borneol. This is the first report that multiple BPPSs exist within a single species. Intriguingly, CcTPS49 and CcTPS72 lead to the generation of dextral-borneol, while CcTPS26 contributes to the biosynthesis of levo-borneol. In addition, the functional characterization of another six CcTPSs suggested that they are responsible for the biosynthesis of linalool, eucalyptol and several other monoterpenes in camphor tree. In conclusion, these novel results provide a foundation for further exploration of the role of the CcTPS gene family and shed light on a better understanding of the biosynthesis and accumulation of monoterpenes in camphor tree.

Background Buddleja lindleyana Fort., which belongs to the Loganiaceae with a distribution throughout the tropics, is widely used as an ornamental plant in China. There are several morphologically similar species in the genus Buddleja, but the lack of comprehensive molecular and phylogenetic studies makes it difficult to distinguish related species, which hinders further studies of this genus. Objective Using molecular biology techniques to sequence and analyze the complete chloroplast (cp) genome of B. lindleyana. Methods After sequencing of the genomic DNA using next-generation sequencing, a series of bioinformatics software were used to assemble and analyze the molecular structure of the cp genome of B. lindleyana. Results The complete cp genome of B. lindleyana is a circular 154 487-bp-long molecule with a GC (Guanine and Cytosine) content of 38.1%. It has a quadripartite structure, including a LSC region (85 489 base pair (bp)), a small single-copy region (17 898 bp), and a pair of inverted repeat regions (25 550 bp). A total of 133 genes were identified in this genome, including 86 protein-coding genes, 37 tRNA (transfer Ribonucleic Acid) genes, eight rRNA (ribosomal Ribonucleic Acid) genes, and two pseudogenes. Conclusion These results suggest that the B. lindelyana cp genome could be used as a potential genomic resource to resolve the phylogenetic positions and relationships of Loganiaceae, and will offer valuable information for future research in the identification of Buddleja species and will conduce to genomic investigations into these species.

Helicteres angustifolia has long been used in Chinese traditional medicine. It has multiple pharmacological benefits, including anti-inflammatory, anti-viral and anti-tumor effects. Its main active chemicals include betulinic acid, oleanolic acid, helicteric acid, helicterilic acid, and other triterpenoid sapions. It is worth noting that the acylation of triterpenoids, such as helicteric acid and helicterilic acid, are characteristic components of Helicteres and are relatively rare among other plants. However, reliance on natural plants as the only sources of these restricts the potential use of H. angustifolia. Therefore, engineering of its metabolic pathway is of high research value for enhancing the production of secondary metabolites. Unfortunately, little is known of the biosynthesis of acylated triterpenoids, hindering its further investigation. Here, the RNAs of different groups treated by metabolic stimulation were sequenced with an Illumina high-throughput sequencing platform, resulting in 121 gigabases of data. A total of 424,824 unigenes were obtained after the trimming and assembly of the raw data, and 22,430 unigenes were determined to be differentially expressed. In addition, three oxidized squalene cyclases (OSCs) and four Cytochrome P450 (CYP450s) were screened, of which one OSC (HaOSC1) and one CYP450 (HaCYPi3) achieved functional verification, suggesting that they could catalyze the production of lupeol and oleanolic acid, respectively. At the same time, we also screened two triterpenoid acetyl transferases (TATs) and one triterpenoid benzoyl transferase (TBT), as subsequent structural modificatory genes, their preliminary study laid a foundation for the acylation of triterpenoids. In a nutshell, these results shed light on the regulation of acylated triterpenoids biosynthesis.

Background The large genus Ficus comprises approximately 800 species, most of which possess high ornamental and ecological values. However, its evolutionary history remains largely unknown. Plastome (chloroplast genome) analysis had become an essential tool for species identification and for unveiling evolutionary relationships between species, genus and other rank groups. In this work we present the plastomes of ten Ficus species. Results The complete chloroplast (CP) genomes of eleven Ficus specimens belonging to ten species were determined and analysed. The full length of the Ficus plastome was nearly 160 kbp with a similar overall GC content, ranging from 35.88 to 36.02%. A total of 114 unique genes, distributed in 80 protein-coding genes, 30 tRNAs, and 4 rRNAs, were annotated in each of the Ficus CP genome. In addition, these CP genomes showed variation in their inverted repeat regions (IR). Tandem repeats and mononucleotide simple sequence repeat (SSR) are widely distributed across the Ficus CP genome. Comparative genome analysis showed low sequence variability. In addition, eight variable regions to be used as potential molecular markers were proposed for future Ficus species identification. According to the phylogenetic analysis, these ten Ficus species were clustered together and further divided into three clades based on different subgenera. Simultaneously, it also showed the relatedness between Ficus and Morus . Conclusion The chloroplast genome structure of 10 Ficus species was similar to that of other angiosperms, with a typical four-part structure. Chloroplast genome sizes vary slightly due to expansion and contraction of the IR region. And the variation of noncoding regions of the chloroplast genome is larger than that of coding regions. Phylogenetic analysis showed that these eleven sampled CP genomes were divided into three clades, clustered with species from subgenus Urostigma , Sycomorus , and Ficus , respectively. These results support the Berg classification system, in which the subgenus Ficus was further decomposed into the subgenus Sycomorus . In general, the sequencing and analysis of Ficus plastomes, especially the ones of species with no or limited sequences available yet, contribute to the study of genetic diversity and species evolution of Ficus , while providing useful information for taxonomic and phylogenetic studies of Ficus .

Necroptosis is a promising novel target for cervical cancer therapy. Nevertheless, differentially expressed necroptosis-related genes (NRGs) in cervical cancer and their associations with prognosis are far from fully clarified. In this study, differentially expressed NRGs (DE-NRGs) were screened out and their bio-function was elucidated. Subsequently, a prognostic scoring model based on the regression coefficients of the screened out NRGs and their corresponding mRNA expressions were constructed and validated. Finally, the survival probability of cervical cancer patients based on the constructed prognostic scoring model in 3 and 5 years was predicted and assessed. We found 17 DE-NRGs in cervical cancer tissues which were closely related to cancer progression, and most of them were significantly highly expressed. Furthermore, 3 NRG were confirmed as the prognostic signature genes from 17 DE-NRGs by regression analysis. Overall survival predicted through our prognostic scoring model was lower in the high-risk group than in the low-risk group (p < 0.05) in both the TCGA cohort and the external GEO44001 validation cohort. What's more, the prediction performance of our prognostic scoring models well verified by the ROC curve, and the risk score calculated could act as an independent prognostic factor for cervical cancer patients. The calibration curve and C-index (0.776) of the nomogram analysis suggested that the predictive performance of the nomogram was satisfactory. Our study identified and validated a necroptosis-related prognostic signature in cervical cancer, which could well predict the prognosis for cervical cancer patients.

Background The borneol chemotype of Cinnamomum camphora (BCC), a monoterpene-rich woody plant species, is the sole source prescribed by the Chinese Pharmacopoeia for the production of natural D-borneol, a major monoterpene in BCC used for millennia as a topical analgesic in China. Nevertheless, the possible gene-regulatory roles of transcription factors (TFs) in BCC’s monoterpenoid biosynthesis remained unknown. Here, a joint analysis of the transcriptome and terpenoid metabolome of BCC induced by mechanical damage (MD) was used to comprehensively explore the interaction between TFs and terpene synthase (TPS) unigenes that might participate in monoterpene biosynthesis in BCC. Results Gas chromatography–mass spectrometry analysis detected 14 monoterpenes and seven sesquiterpenes. All but two monoterpenes underwent a significantly increased accumulation after the MD treatment. RNA sequencing data revealed that 10 TPS, 82 MYB, 70 AP2/ERF, 38 BHLH, 31 WRKY, and 29 bZIP unigenes responded to the MD treatment. A correlation analysis revealed that three monoterpene synthase genes (CcTPS1, CcTPS3, CcTPS4) highly correlated with multiple monoterpenes, namely D-borneol, camphor, and bornyl acetate, which could be responsible for monoterpenoid biosynthesis in BCC. Furthermore, five WRKY, 15 MYB, 10 ERF/AP2, five bZIP, and two BHLH genes had strong, positive correlations with CcTPS1 or CcTPS4, judging by their high coefficient values (R2 > 0.8). The bioinformatics results were verified by quantitative real-time PCR. Conclusion This study provides insight into the genes involved in the biosynthesis and regulation of monoterpene in BCC and thus provides a pool of candidate genes for future mechanistic analyses of how monoterpenes accumulate in BCC.

Ardisia Sw. (Primulaceae) is naturally distributed in tropical and subtropical areas. Most of them possess edible and medicinal values and are popular in clinical and daily use in China. However, ambiguous species delineation and genetic information limit the development and utilization of this genus. In this study, the chloroplast genomes of four Ardisia species, namely A. gigantifolia Stapf, A. crenata Sims, A. villosa Roxb. and A. mamillata Hance, were sequenced, annotated, and analyzed comparatively. All the four chloroplast genomes possess a typical quadripartite structure, and each of the genomes is about 156 Kb in size. The structure and gene content of the Ardisia plastomes were conservative and showed low sequence divergence. Furthermore, we identified five mutation hotspots as candidate DNA barcodes for Ardisia, namely, trnT-psbD, ndhF-rpl32, rpl32-ccsA, ccsA-ndhD and ycf1. Phylogenetic analysis based on the whole-chloroplast genomes data showed that Ardisia was sister to Tapeinosperma Hook. f. In addition, the results revealed a great topological profile of Ardisia's with strong support values, which matches their geographical distribution patterns. Summarily, our results provide useful information for investigations on taxonomic differences, molecular identification, and phylogenetic relationships of Ardisia plants.

We investigated the antitumor effects of four fractions of Dendrobium officinale Kimura & Migo (D. officinale) polysaccharides with different molecular weights (Mw), Astragalus membranaceus polysaccharides (APS) and Lentinus edodes polysaccharides (LNT) on colorectal cancer (CRC) using a zebrafish xenograft model. Transcriptome sequencing was performed to further explore the possible antitumor mechanisms of D. officinale polysaccharides. Fractions of D. officinale polysaccharides, LNT, and APS could significantly inhibit the growth of HT-29 cells in a zebrafish xenograft model. One fraction of D. officinale polysaccharides called DOPW-1 (Mw of 389.98 kDa) exhibited the strongest tumor inhibition. Compared with the control group, RNA-seq revealed that the DOPW-1–treated experimental group had 119 differentially expressed genes (DEGs), of which 45 had upregulated expression and 74 had downregulated expression. Analyses using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes suggested that the pathway “apoptosis-multiple species” was the most significantly enriched. Our data indicated that 1) fractions of D. officinale polysaccharides of Mw 389.98 kDa were most suitable against CRC; 2) DOPW-1 could be developed into a clinical agent against CRC; and 3) an apoptosis pathway is important for DOPW-1 to inhibit the proliferation of HT-29 cells.

Trans-isopentenyl diphosphate synthases (TIDSs) genes are known to be important determinants for terpene diversity and the accumulation of terpenoids. The essential oil of Cinnamomum camphora , which is rich in monoterpenes, sesquiterpenes, and other aromatic compounds, has a wide range of pharmacological activities and has therefore attracted considerable interest. However, the TIDS gene family, and its relationship to the camphor tree ( C. camphora L. Presl.), has not yet been characterized. In this study, we identified 10 TIDS genes in the genome of the C. camphora borneol chemotype that were unevenly distributed on chromosomes. Synteny analysis revealed that the TIDS gene family in this species likely expanded through segmental duplication events. Furthermore, cis-element analyses demonstrated that C. camphora TIDS (CcTIDS) genes can respond to multiple abiotic stresses. Finally, functional characterization of eight putative short-chain TIDS proteins revealed that CcTIDS3 and CcTIDS9 exhibit farnesyl diphosphate synthase (FPPS) activity, while CcTIDS1 and CcTIDS2 encode geranylgeranyl diphosphate synthases (GGPPS). Although, CcTIDS8 and CcTIDS10 were found to be catalytically inactive alone, they were able to bind to each other to form a heterodimeric functional geranyl diphosphate synthase (GPPS) in vitro , and this interaction was confirmed using a yeast two-hybrid assay. Furthermore, transcriptome analysis revealed that the CcTIDS3 , CcTIDS8 , CcTIDS9 , and CcTIDS10 genes were found to be more active in C. camphora roots as compared to stems and leaves, which were verified by quantitative real-time PCR (qRT-PCR). These novel results provide a foundation for further exploration of the role of the TIDS gene family in camphor trees, and also provide a potential mechanism by which the production of camphor tree essential oil could be increased for pharmacological purposes through metabolic engineering.

Background Dextrorotatory borneol (D-borneol), a cyclic monoterpene, is widely used in traditional Chinese medicine as an efficient topical analgesic drug. Fresh leaves of Cinnamomum trees, e.g., C. burmannii and C. camphor, are the main sources from which D-borneol is extracted by steam distillation, yet with low yields. Insufficient supply of D-borneol has hampered its clinical use and production of patent remedies for a long time. Biological synthesis of D-borneol offers an additional approach; however, mechanisms of D-borneol biosynthesis remain mostly unresolved. Hence, it is important and necessary to elucidate the biosynthetic pathway of D-borneol. Results Comparative analysis on the gene expression patterns of different D-borneol production C. burmannii samples facilitates elucidation on the underlying biosynthetic pathway of D-borneol. Herein, we collected three different chemotypes of C. burmannii, which harbor different contents of D-borneol.A total of 100,218 unigenes with an N50 of 1,128 bp were assembled de novo using Trinity from a total of 21.21 Gb clean bases. We used BLASTx analysis against several public databases to annotate 45,485 unigenes (45.38%) to at least one database, among which 82 unigenes were assigned to terpenoid biosynthesis pathways by KEGG annotation. In addition, we defined 8,860 unigenes as differentially expressed genes (DEGs), among which 13 DEGs were associated with terpenoid biosynthesis pathways. One 1-deoxy-D-xylulose-5-phosphate synthase (DXS) and two monoterpene synthase, designated as CbDXS9, CbTPS2 and CbTPS3, were up-regulated in the high-borneol group compared to the low-borneol and borneol-free groups, and might be vital to biosynthesis of D-borneol in C. burmannii. In addition, we identified one WRKY, two BHLH, one AP2/ERF and three MYB candidate genes, which exhibited the same expression patterns as CbTPS2 and CbTPS3, suggesting that these transcription factors might potentially regulate D-borneol biosynthesis. Finally, quantitative real-time PCR was conducted to detect the actual expression level of those candidate genes related to the D-borneol biosynthesis pathway, and the result showed that the expression patterns of the candidate genes related to D-borneol biosynthesis were basically consistent with those revealed by transcriptome analysis. Conclusions We used transcriptome sequencing to analyze three different chemotypes of C. burmannii, identifying three candidate structural genes (one DXS, two monoterpene synthases) and seven potential transcription factor candidates (one WRKY, two BHLH, one AP2/ERF and three MYB) involved in D-borneol biosynthesis. These results provide new insight into our understanding of the production and accumulation of D-borneol in C. burmannii.