Genome-wide profiling of WRKY genes involved in flavonoid biosynthesis in Erigeron breviscapus.

The transcription factors of WRKY genes play essential roles in plant growth, stress responses, and metabolite biosynthesis. Erigeron breviscapus, a traditional Chinese herb, is abundant in flavonoids and has been used for centuries to treat cardiovascular and cerebrovascular diseases. However, the WRKY transcription factors that regulate flavonoid biosynthesis in E. breviscapus remain unknown. In this study, a total of 75 EbWRKY transcription factors were predicted through comprehensive genome-wide characterization of E. breviscapus and the chromosomal localization of each EbWRKY gene was investigated. RNA sequencing revealed transient responses of 74 predicted EbWRKY genes to exogenous abscisic acid (ABA), salicylic acid (SA), and gibberellin 3 (GA3) after 4 h of treatment. In contrast, the expression of key structural genes involved in flavonoid biosynthesis increased after 4 h in GA3 treatment. However, the content of flavonoid metabolites in leaves significantly increased at 12 h. The qRT-PCR results showed that the expression patterns of EbWRKY11, EbWRKY30, EbWRKY31, EbWRKY36, and EbWRKY44 transcription factors exhibited a high degree of similarity to the 11 structural genes involved in flavonoid biosynthesis. Protein-DNA interactions were performed between the key genes involved in scutellarin biosynthesis and candidate WRKYs. The result showed that F7GAT interacts with EbWRKY11, EbWRKY36, and EbWRKY44, while EbF6H has a self-activation function. This study provides comprehensive information on the regulatory control network of flavonoid accumulation mechanisms, offering valuable insights for breeding E. breviscapus varieties with enhanced scutellarin content.

Functional characterization of genes related to triterpene and flavonoid biosynthesis in Cyclocarya paliurus.

MAIN CONCLUSION: The oxidosqualene cyclases (OSCs) generating triterpenoid skeletons in Cyclocarya paliurus were identified for the first time, and two uridine diphosphate (UDP)-glycosyltransferases (UGTs) catalyzing the glycosylation of flavonoids were characterized. Cyclocarya paliurus, a native rare dicotyledonous plant in China, contains an abundance of triterpenoid saponins and flavonoid glycosides that exhibit valuable pharmaceutical effects in preventing hypertension, hyperlipidemia, and diabetes. However, the molecular mechanism explaining the biosynthesis of triterpenoid saponin and flavonoid glycoside in C. paliurus remains unclear. In this study, the triterpene content in different tissues and the expression pattern of genes encoding the key enzymes associated with triterpenoid saponin and flavonoid glycoside biosynthesis were studied using transcriptome and metabolome analysis. The eight upstream oxidosqualene cyclases (OSCs) involved in triterpenoid saponin biosynthesis were functionally characterized, among them CpalOSC6 catalyzed 2,3;22,23-dioxidosqualene to form 3-epicabraleadiol; CpalOSC8 cyclized 2,3-oxidosqualene to generate dammarenediol-II; CpalOSC2 and CpalOSC3 produced ß-amyrin and CpalOSC4 produced cycloartenol, while CpalOSC2-CpalOSC5, CpalOSC7, and CpalOSC8 all produced lanosterol. However, no catalytic product was detected for CpalOSC1. Moreover, two downstream flavonoid uridine diphosphate (UDP)-glycosyltransferases (UGTs) (CpalUGT015 and CpalUGT100) that catalyze the last step of flavonoid glycoside biosynthesis were functionally elucidated. These results uncovered the key genes involved in the biosynthesis of triterpenoid saponins and flavonoid glycosides in C. paliurus that could be applied to produce flavonoid glycosides and key triterpenoid saponins in the future via a synthetic strategy.

LipIDens: simulation assisted interpretation of lipid densities in cryo-EM structures of membrane proteins.

Cryo-electron microscopy (cryo-EM) enables the determination of membrane protein structures in native-like environments. Characterising how membrane proteins interact with the surrounding membrane lipid environment is assisted by resolution of lipid-like densities visible in cryo-EM maps. Nevertheless, establishing the molecular identity of putative lipid and/or detergent densities remains challenging. Here we present LipIDens, a pipeline for molecular dynamics (MD) simulation-assisted interpretation of lipid and lipid-like densities in cryo-EM structures. The pipeline integrates the implementation and analysis of multi-scale MD simulations for identification, ranking and refinement of lipid binding poses which superpose onto cryo-EM map densities. Thus, LipIDens enables direct integration of experimental and computational structural approaches to facilitate the interpretation of lipid-like cryo-EM densities and to reveal the molecular identities of protein-lipid interactions within a bilayer environment. We demonstrate this by application of our open-source LipIDens code to ten diverse membrane protein structures which exhibit lipid-like densities.

Comparative genomics reveals the diversification of triterpenoid biosynthesis and origin of ocotillol-type triterpenes in Panax.

Gene duplication is assumed to be the major force driving the evolution of metabolite biosynthesis in plants. Freed from functional burdens, duplicated genes can mutate toward novelties until fixed due to selective fitness. However, the extent to which this mechanism has driven the diversification of metabolite biosynthesis remains to be tested. Here we performed comparative genomics analysis and functional characterization to evaluate the impact of gene duplication on the evolution of triterpenoid biosynthesis using Panax species as models. We found that whole-genome duplications (WGDs) occurred independently in Araliaceae and Apiaceae lineages. Comparative genomics revealed the evolutionary trajectories of triterpenoid biosynthesis in plants, which was mainly promoted by WGDs and tandem duplication. Lanosterol synthase (LAS) was likely derived from a tandem duplicate of cycloartenol synthase that predated the emergence of Nymphaeales. Under episodic diversifying selection, the LAS gene duplicates produced by γ whole-genome triplication have given rise to triterpene biosynthesis in core eudicots through neofunctionalization. Moreover, functional characterization revealed that oxidosqualene cyclases (OSCs) responsible for synthesizing dammarane-type triterpenes in Panax species were also capable of producing ocotillol-type triterpenes. Genomic and biochemical evidence suggested that Panax genes encoding the above OSCs originated from the specialization of one OSC gene duplicate produced from a recent WGD shared by Araliaceae (Pg-ß). Our results reveal the crucial role of gene duplication in diversification of triterpenoid biosynthesis in plants and provide insight into the origin of ocotillol-type triterpenes in Panax species.

BMF-AS1/BMF Promotes Diabetic Vascular Calcification and Aging both In Vitro and In Vivo.

Vascular calcification and aging often increase morbidity and mortality in patients with diabetes mellitus (DM); however, the underlying mechanisms are still unknown. In the present study, we found that Bcl-2 modifying factor (BMF) and BMF antisense RNA 1 (BMF-AS1) were significantly increased in high glucose-induced calcified and senescent vascular smooth muscle cells (VSMCs) as well as artery tissues from diabetic mice. Inhibition of BMF-AS1 and BMF reduced the calcification and senescence of VSMCs, whereas overexpression of BMF-AS1 and BMF generates the opposite results. Mechanistic analysis showed that BMF-AS1 interacted with BMF directly and up-regulated BMF at both mRNA and protein levels, but BMF did not affect the expression of BMF-AS1. Moreover, knocking down BMF-AS1 and BMF suppressed the calcification and senescence of VSMCs, and BMF knockout (BMF-/-) diabetic mice presented less vascular calcification and aging compared with wild type diabetic mice. In addition, higher coronary artery calcification scores (CACs) and increased plasma BMF concentration were found in patients with DM, and there was a positive correlation between CACs and plasma BMF concentration. Thus, BMF-AS1/BMF plays a key role in promoting high glucose-induced vascular calcification and aging both in vitro and in vivo. BMF-AS1 and BMF represent potential therapeutic targets in diabetic vascular calcification and aging.

Mass spectrometry captures biased signalling and allosteric modulation of a G-protein-coupled receptor.

G-protein-coupled receptors signal through cognate G proteins. Despite the widespread importance of these receptors, their regulatory mechanisms for G-protein selectivity are not fully understood. Here we present a native mass spectrometry-based approach to interrogate both biased signalling and allosteric modulation of the ß1-adrenergic receptor in response to various ligands. By simultaneously capturing the effects of ligand binding and receptor coupling to different G proteins, we probed the relative importance of specific interactions with the receptor through systematic changes in 14 ligands, including isoprenaline derivatives, full and partial agonists, and antagonists. We observed enhanced dynamics of the intracellular loop 3 in the presence of isoprenaline, which is capable of acting as a biased agonist. We also show here that endogenous zinc ions augment the binding in receptor-Gs complexes and propose a zinc ion-binding hotspot at the TM5/TM6 intracellular interface of the receptor-Gs complex. Further interrogation led us to propose a mechanism in which zinc ions facilitate a structural transition of the intermediate complex towards the stable state.

Genome-wide identification of bHLH transcription factors: Discovery of a candidate regulator related to flavonoid biosynthesis in Erigeron breviscapus.

Erigeron breviscapus is a Compositae plant, and its rich flavonoids have shown strong preventative and curative effects in the treatment of cardio- and cerebrovascular diseases. bHLH genes play a crucial role in plant growth and development. There are 116 EbbHLH genes in E. breviscapus, and each gene has been named based on its chromosome location. Our phylogenetic analysis divided these genes into 18 subfamilies. To further investigate its function, EbbHLH80 was isolated from E. breviscapus leaves. Next, transcriptomic and metabolomic analyses of tobacco leaves were performed. Among 421 differentially accumulated compounds, 98 flavonoids were identified. In addition, differentially expressed genes were identified using RNA-seq, and further analysis suggested that EbbHLH80-OE could not only regulate the expression of some structural genes in the flavonoid biosynthesis pathway to achieve flavonoid accumulation but also be involved in the regulation of a series of downstream pathways, such as stress response, ABA and ethylene signal transduction, to affect plant growth and development. The results of our analysis provide new insights into the function of EbbHLH80 and lay the foundation for future functional studies on E. breviscapus.

EbMYBP1, a R2R3-MYB transcription factor, promotes flavonoid biosynthesis in Erigeron breviscapus.

Erigeron breviscapus, a traditional Chinese medicinal plant, is enriched in flavonoids that are beneficial to human health. While we know that R2R3-MYB transcription factors (TFs) are crucial to flavonoid pathway, the transcriptional regulation of flavonoid biosynthesis in E. breviscapus has not been fully elucidated. Here, EbMYBP1, a R2R3-MYB transcription factor, was uncovered as a regulator involved in the regulation of flavonoid accumulation. Transcriptome and metabolome analysis revealed that a large group of genes related to flavonoid biosynthesis were significantly changed, accompanied by significantly increased concentrations of the flavonoid in EbMYBP1-OE transgenic tobacco compared with the wild-type (WT). In vitro and in vivo investigations showed that EbMYBP1 participated in flavonoid biosynthesis, acting as a nucleus-localized transcriptional activator and activating the transcription of flavonoid-associated genes like FLS, F3H, CHS, and CHI by directly binding to their promoters. Collectively, these new findings are advancing our understanding of the transcriptional regulation that modulates the flavonoid biosynthesis.

Hepatocyte cholesterol content modulates glucagon receptor signalling.

OBJECTIVE: To determine whether glucagon receptor (GCGR) actions are modulated by cellular cholesterol levels. METHODS: We determined the effects of experimental cholesterol depletion and loading on glucagon-mediated cAMP production, ligand internalisation and glucose production in human hepatoma cells, mouse and human hepatocytes. GCGR interactions with lipid bilayers were explored using coarse-grained molecular dynamic simulations. Glucagon responsiveness was measured in mice fed a high cholesterol diet with or without simvastatin to modulate hepatocyte cholesterol content. RESULTS: GCGR cAMP signalling was reduced by higher cholesterol levels across different cellular models. Ex vivo glucagon-induced glucose output from mouse hepatocytes was enhanced by simvastatin treatment. Mice fed a high cholesterol diet had increased hepatic cholesterol and a blunted hyperglycaemic response to glucagon, both of which were partially reversed by simvastatin. Simulations identified likely membrane-exposed cholesterol binding sites on the GCGR, including a site where cholesterol is a putative negative allosteric modulator. CONCLUSIONS: Our results indicate that cellular cholesterol content influences glucagon sensitivity and indicate a potential molecular basis for this phenomenon. This could be relevant to the pathogenesis of non-alcoholic fatty liver disease, which is associated with both hepatic cholesterol accumulation and glucagon resistance.

Integrated Metabolomic and Transcriptomic Analysis and Identification of Dammarenediol-II Synthase Involved in Saponin Biosynthesis in Gynostemma longipes.

Gynostemma longipes contains an abundance of dammarane-type ginsenosides and gypenosides that exhibit extensive pharmacological activities. Increasing attention has been paid to the elucidation of cytochrome P450 monooxygenases (CYPs) and UDP-dependent glycosyltransferases (UGTs) that participate downstream of ginsenoside biosynthesis in the Panax genus. However, information on oxidosqualene cyclases (OSCs), the upstream genes responsible for the biosynthesis of different skeletons of ginsenoside and gypenosides, is rarely reported. Here, an integrative study of the metabolome and the transcriptome in the leaf, stolon, and rattan was conducted and the function of GlOSC1 was demonstrated. In total, 46 triterpenes were detected and found to be highly abundant in the stolon, whereas gene expression analysis indicated that the upstream OSC genes responsible for saponin skeleton biosynthesis were highly expressed in the leaf. These findings indicated that the saponin skeletons were mainly biosynthesized in the leaf by OSCs, and subsequently transferred to the stolon via CYPs and UGTs biosynthesis to form various ginsenoside and gypenosides. Additionally, a new dammarane-II synthase (DDS), GlOSC1, was identified by bioinformatics analysis, yeast expression assay, and enzyme assays. The results of the liquid chromatography-mass spectrometry (LC-MS) analysis proved that GlOSC1 could catalyze 2,3-oxidosqualene to form dammarenediol-II via cyclization. This work uncovered the biosynthetic mechanism of dammarenediol-II, an important starting substrate for ginsenoside and gypenosides biosynthesis, and may achieve the increased yield of valuable ginsenosides and gypenosides produced under excess substrate in a yeast cell factory through synthetic biology strategy.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations.

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

[Identification and expression profiling of R2R3-MYB transcription factors in Erigeron breviscapus].

R2 R3-MYB transcription factors are ubiquitous in plants, playing a role in the regulation of plant growth, development, and secondary metabolism. In this paper, the R2 R3-MYB transcription factors were identified by bioinformatics analysis of the genomic data of Erigeron breviscapus, and their gene sequences, structures, physical and chemical properties were analyzed. The functions of R2 R3-MYB transcription factors were predicted by cluster analysis. Meanwhile, the expression patterns of R2 R3-MYB transcription factors in response to hormone treatments were analyzed. A total of 108 R2 R3-MYB transcription factors, named EbMYB1-EbMYB108, were identified from the genome of E. breviscapus. Most of the R2 R3-MYB genes carried 2-4 exons. The phylogenetic tree of MYBs in E. breviscapus and Arabidopsis thaliala was constructed, which classified 234 MYBs into 30 subfamilies. The MYBs in the five MYB subfamilies of A.thaliala were clustered into independent clades, and those in E. breviscapus were clustered into four clades. The transcriptome data showed that MYB genes were differentially expressed in different tissues of E. breviscapus and in response to the treatments with exogenous hormones such as ABA, SA, and GA for different time. The transcription of 13 R2 R3-MYB genes did not change significantly, and the expression patterns of some genes were up-regulated or down-regulated with the extension of hormone treatment time. This study provides a theoretical basis for revealing the mechanisms of R2 R3-MYB transcription factors in regulating the growth and development, stress(hormone) response, and active ingredient accumulation in E. breviscapus.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations.

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins.

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

Modulation of adenosine A2a receptor oligomerization by receptor activation and PIP2 interactions.

GPCRs have been shown to form oligomers, which generate distinctive signaling outcomes. However, the structural nature of the oligomerization process remains uncertain. We have characterized oligomeric configurations of the adenosine A2a receptor (A2aR) by combining large-scale molecular dynamics simulations with Markov state models. These oligomeric structures may also serve as templates for studying oligomerization of other class A GPCRs. Our simulation data revealed that receptor activation results in enhanced oligomerization, more diverse oligomer populations, and a more connected oligomerization network. The active state conformation of the A2aR shifts protein-protein association interfaces to those involving intracellular loop ICL3 and transmembrane helix TM6. Binding of PIP2 to A2aR stabilizes protein-protein interactions via PIP2-mediated association interfaces. These results indicate that A2aR oligomerization is responsive to the local membrane lipid environment. This, in turn, suggests a modulatory effect on A2aR whereby a given oligomerization profile favors the dynamic formation of specific supramolecular signaling complexes.

Identification of two UDP-glycosyltransferases involved in the main oleanane-type ginsenosides in Panax japonicus var. major.

MAIN CONCLUSION: Two UDP-glycosyltransferases from Panax japonicus var. major were identified, and the biosynthetic pathways of three oleanane-type ginsenosides (chikusetsusaponin IVa, ginsenoside Ro, zingibroside R1) were elucidated. Chikusetsusaponin IVa and ginsenoside Ro are primary active components formed by stepwise glycosylation of oleanolic acid in five medicinal plants of the genus Panax. However, the key UDP-glycosyltransferases (UGTs) in the biosynthetic pathway of chikusetsusaponin IVa and ginsenoside Ro are still unclear. In this study, two UGTs (PjmUGT1 and PjmUGT2) from Panax japonicus var. major involved in the biosynthesis of chikusetsusaponin IVa and ginsenoside Ro were identified based on bioinformatics analysis, heterologous expression and enzyme assays. The results show that PjmUGT1 can transfer a glucose moiety to the C-28 carboxyl groups of oleanolic acid 3-O-ß-D-glucuronide and zingibroside R1 to form chikusetsusaponin IVa and ginsenoside Ro, respectively. Meanwhile, PjmUGT2 can transfer a glucose moiety to oleanolic acid 3-O-ß-D-glucuronide and chikusetsusaponin IVa to form zingibroside R1 and ginsenoside Ro. This work uncovered the biosynthetic mechanism of chikusetsusaponin IVa and ginsenoside Ro, providing the rational production of valuable saponins through synthetic biology strategy.

Improved de novo Assembly of the Achlorophyllous Orchid Gastrodia elata.

Achlorophyllous plants are full mycoheterotrophic plants with no chlorophyll and they obtain their nutrients from soil fungi. Gastrodia elata is a perennial, achlorophyllous orchid that displays distinctive evolutionary strategy of adaptation to the non-photosynthetic lifestyle. Here in this study, the genome of G. elata was assembled to 1.12 Gb with a contig N50 size of 110 kb and a scaffold N50 size of 1.64 Mb so that it helped unveil the genetic basics of those adaptive changes. Based on the genomic data, key genes related to photosynthesis, leaf development, and plastid division pathways were found to be lost or under relaxed selection during the course of evolution. Thus, the genome sequence of G. elata provides a good resource for future investigations of the evolution of orchids and other achlorophyllous plants.

The Glycosphingolipid GM3 Modulates Conformational Dynamics of the Glucagon Receptor.

The extracellular domain (ECD) of class B1 G-protein-coupled receptors (GPCRs) plays a central role in signal transduction and is uniquely positioned to sense both the extracellular and membrane environments. Although recent studies suggest a role for membrane lipids in the modulation of class A and class F GPCR signaling properties, little is known about the effect of lipids on class B1 receptors. In this study, we employed multiscale molecular dynamics simulations to access the dynamics of the glucagon receptor (GCGR) ECD in the presence of native-like membrane bilayers. Simulations showed that the ECD could move about a hinge region formed by residues Q122-E126 to adopt both closed and open conformations relative to the transmembrane domain. ECD movements were modulated by binding of the glycosphingolipid GM3. These large-scale fluctuations in ECD conformation may affect the ligand binding and receptor activation properties. We also identify a unique phosphatidylinositol (4,5)-bisphosphate (PIP2) interaction profile near intracellular loop (ICL) 2/TM3 at the G-protein-coupling interface, suggesting a mechanism of engaging G-proteins that may have a distinct dependence on PIP2 compared with class A GPCRs. Given the structural conservation of class B1 GPCRs, the modulatory effects of GM3 and PIP2 on GCGR may be conserved across these receptors, offering new insights into potential therapeutic targeting.

Lipid-Dependent Regulation of Ion Channels and G Protein-Coupled Receptors: Insights from Structures and Simulations.

Ion channels and G protein-coupled receptors (GPCRs) are regulated by lipids in their membrane environment. Structural studies combined with biophysical and molecular simulation investigations reveal interaction sites for specific lipids on membrane protein structures. For K channels, PIP2 plays a key role in regulating Kv and Kir channels. Likewise, several recent cryo-EM structures of TRP channels have revealed bound lipids, including PIP2 and cholesterol. Among the pentameric ligand-gated ion channel family, structural and biophysical studies suggest the M4 TM helix may act as a lipid sensor, e.g., forming part of the binding sites for neurosteroids on the GABAA receptor. Structures of GPCRs have revealed multiple cholesterol sites, which may modulate both receptor dynamics and receptor oligomerization. PIP2 also interacts with GPCRs and may modulate their interactions with G proteins. Overall, it is evident that multiple lipid binding sites exist on channels and receptors that modulate their function allosterically and are potential druggable sites.

Transcriptome analysis of Panax zingiberensis identifies genes encoding oleanolic acid glucuronosyltransferase involved in the biosynthesis of oleanane-type ginsenosides.

MAIN CONCLUSION: Oleanolic acid glucuronosyltransferase (OAGT) genes synthesizing the direct precursor of oleanane-type ginsenosides were discovered. The four recombinant proteins of OAGT were able to transfer glucuronic acid at C-3 of oleanolic acid that yields oleanolic acid 3-O-ß-glucuronide. Ginsenosides are the primary active components in the genus Panax, and great efforts have been made to elucidate the mechanisms underlying dammarane-type ginsenoside biosynthesis. However, there is limited information on oleanane-type ginsenosides. Here, high-performance liquid chromatography analysis demonstrated that oleanane-type ginsenosides (particularly ginsenoside Ro and chikusetsusaponin IV and IVa) are the abundant ginsenosides in Panax zingiberensis, an extremely endangered Panax species in southwest China. These ginsenosides are derived from oleanolic acid 3-O-ß-glucuronide, which may be formed from oleanolic acid catalyzed by an unknown oleanolic acid glucuronosyltransferase (OAGT). Transcriptomic analysis of leaves, stems, main roots, and fibrous roots of P. zingiberensis was performed, and a total of 46,098 unigenes were obtained, including all the identified homologous genes involved in ginsenoside biosynthesis. The most upstream genes were highly expressed in the leaves, and the UDP-glucosyltransferase genes were highly expressed in the roots. This finding indicated that the precursors of ginsenosides are mainly synthesized in the leaves and transported to different parts for the formation of particular ginsenosides. For the first time, enzyme activity assay characterized four genes (three from P. zingiberensis and one from P. japonicus var. major, another Panax species with oleanane-type ginsenosides) encoding OAGT, which particularly transfer glucuronic acid at C-3 of oleanolic acid to form oleanolic acid 3-O-ß-glucuronide. Taken together, our study provides valuable genetic information for P. zingiberensis and the genes responsible for synthesizing the direct precursor of oleanane-type ginsenosides.