Transcriptome analysis of Stephania yunnanensis and functional validation of CYP80s involved in benzylisoquinoline alkaloids biosynthesis

Transcriptome analysis of Stephania yunnanensis and functional validation of CYP80s involved in benzylisoquinoline alkaloids biosynthesis | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Transcriptome analysis of Stephania yunnanensis and functional validation of CYP80s involved in benzylisoquinoline alkaloids biosynthesis Wenlong Shi, Qishuang Li, Xinyi Li, Jingyi Gan, Ying Ma, Jian Wang, and 6 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5384973/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract The medicinal plant Stephania yunnanensis is rich in aporphine alkaloids, a type of benzylisoquinoline alkaloids (BIAs), with aporphine being the representative and most abundant compound, but our understanding on the biosynthesis of BIA alkaloids in this plant have been relatively limited. Previous research has reported the genome of S. yunnanensis and preliminarily identified the upstream gene norcoclaurine synthase (NCS) in the BIA biosynthetic pathways. However, the key genes promoting the formation of the aporphine skeleton have not yet been reported. In this study, based on the differences in the content of crebanine and several other BIAs in different tissues, we conducted transcriptome sequencing of roots, stems, and leaves. We then identified candidate genes through functional annotation and sequence alignment, followed by transcriptomic and genomic analyses. Based on this analysis, we identified three CYP80 enzymes (SyCYP80Q5-1, SyCYP80Q5-3, and SyCYP80G6), which exhibited different activities towards ( S )- and ( R )-configured substrates in S. yunnanensis and demonstrated strict stereoselectivity enroute to aporphine. This study provides metabolomic and transcriptomic information on the biosynthesis of BIAs in S. yunnanensis and offers valuable insights into the elucidation of BIA biosynthesis, and lays the foundation for the complete analysis of pathways for more aporphine alkaloids. Molecular Biology Plant Molecular Biology and Genetics aporphine alkaloids biosynthetic gene clusters biosynthesis CYP80 Stephania yunnanensis transcriptome Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Key Message This study found that BIA biosynthetic genes in Stephania yunnanensis are not all highly expressed in roots and are clustered on chromosomes. Two CYP80s catalyze key steps in aporphine biosynthesis with strict conformational selectivity. 1. Introduction Stephania yunnanensis is one of the source plants of the traditional Chinese medicine ‘Shanwugui’, which is a plant from the genus Stephania in the family Menispermaceae of the order Ranales (Heming & Xianrui, 1980 ). Traditional Chinese medicine theory believes that S. yunnanensis has the effects of vomiting phlegm and foods, treating malaria, and detoxifying sores (Zhao et al., 2021 ). Previous chemical research has shown that the main active components of S. yunnanensis were four types of benzylisoquinoline alkaloids (BIAs), including 1-benzylisoquinolines (1-BIAs), aporphines, morphinans, and protoberberines (Dai et al., 2012 ; Hu et al., 2010 ; Wang et al., 2022 ), such as salutaridine, tetrahydropalmatine, roemerine, stephanine, and crebanine, etc. Crebanine, a kind of aporphine, was major active compound found in S. yunnanensis (Peng, 2014 ), recent studies have indicated that crebanine has neuroprotective (Yang et al., 2023 ), cardioprotective(Wang et al., 2016 ), anticancer (Tan et al., 2023 ; Yeh et al., 2024 ), anti-inflammatory and analgesic activities (Cui et al., 2022 ). Since the extraction of BIAs from plants was an expensive process and the content was also affected by the growth conditions of the plants, this low yield and variability limited the commercial development of BIAs. The complexity of chemical synthesis further hindered the large-scale production of BIAs. The demand for BIAs ultimately promoted the development of BIAs metabolic engineering and synthetic biology. The biosynthetic pathways of BIAs were initiated with the enzymatic reaction of norcoclaurine synthase (NCS), which catalyzes the synthesis of norcoclaurine from dopamine and 4-hydroxyphenylacetic acid (4-HPAA) (Ghirga et al., 2017 ). Subsequent biochemical transformations, included hydroxylation, methylation, and isomerization processes involving an array of cytochrome P450s, methyltransferases, and reductases, such as norcoclaurine 6- O -methyltransferase (6OMT) (Menendez-Perdomo & Facchini, 2020 ), coclaurine N -methyltransferase (CNMT), ( S )- N -methylcoclaurine 3’-hydroxylase (NMCH) (Liu et al., 2021 ; Pauli & Kutchan, 1998 ), and 3'-hydroxy- N -methylcoclaurine 4'- O -methyltransferase (4'OMT) (Gurkok et al., 2016 ). These enzymes catalyzed the production of a series of crucial 1-BIA intermediates, including coclaurine (COC), N -methylcoclaurine (NMC), 3'-hydroxyl- N -methylcoclaurine (HNMC), and reticuline, which subsequently served as substrates for the biosynthesis of a diverse array of BIAs (Fig. 4 ). Aporphine alkaloids like crebanine were the major compounds in S. yunnanensis , which had an unclear biosynthetic pathway. Recent studies have discovered that corytuberine synthase (CTS) like CjCYP80G2 from Coptis japonica (Ikezawa et al., 2008 ) and StCYP80G6 (Li et al., 2024 ) from Stephania tetrandra could catalytic C-C coupling of ( S )-reticuline to the type Ⅰ aporphine ( S )-corytuberine. Additionally, NnCYP80Q1 or StCYP80Q5 could catalytic C-C coupling of ( R )-coclaurine and ( R )- N -methylcoclaurine to the protoaporphines ( R )-crotoflorine and ( R )-glaziovine, two important precursors of type Ⅱ aporphines (Li et al., 2024 ; Pyne et al., 2023 ). Subsequently, the proaporphines crotoflorine and glaziovine underwent unknown oxidative rearrangements, followed by a series of hydroxylations, methylations, and the formation of methylenedioxy bridges, leading to the production of type II aporphines, including crebanine and nelumboferine (Li et al., 2024 ). These CYP80s were considered to be the starting points of the aporphines biosynthetic pathways, however, no CYP80s with similar catalytic C-C coupling functions had been found in S. yunnanensis . The combination of transcriptomics and metabolomics is currently the most common method for studying the biosynthesis and regulation of plants secondary metabolites (Lai et al., 2024 ; Li et al., 2022 ; Li et al., 2023 ; Tong et al., 2022 ). To date, de novo transcriptomes of multiple BIA-producing plants have been reported, focusing on the Papaveraceae, Ranunculaceae, and Menispermaceae families within the Ranunculales (Hagel et al., 2015 ). Papaveraceae plants have been a primary focus due to their unique accumulation of morphinan alkaloids (Park et al., 2024 ; Pei et al., 2021 ; Xu et al., 2022 ; Zhao et al., 2019 ). As the pharmacological effects of other types BIAs have been discovered, the transcriptomes of other BIA-producing plants have also been reported to enable more in-depth research into BIA biosynthesis. The genome of S. yunnanensis has been previously reported (Leng et al., 2024 ), providing valuable genomic data for investigating its BIA biosynthetic pathways. However, for a deeper understanding of the BIA biosynthesis and its regulation, it is also necessary to have information on gene expression levels and metabolite profiles. Unfortunately, there are currently no reported transcriptome data available for S. yunnanensis . In this study, we determined the relative concentrations of four BIAs in different tissues of S. yunnanensis to investigate the differential accumulation patterns of these BIAs. We also generated a de novo transcriptome of S. yunnanensis and identified hundreds of genes potentially related to BIA biosynthesis. Subsequently, using differential expression analysis and phylogenetic analysis, we identified some candidate genes involved in BIA biosynthesis in S. yunnanensis and conducted genomic-level analysis on them. We then validated the functions of CYP80s genes of S. yunnanensis in vitro . Our research will contribute to further in-depth studies on BIA biosynthetic pathways and their evolution in Stephania species, aiming to elucidate the complete biosynthetic pathways of BIAs represented by crebanine. 2. Materials and methods 2.1 Plant materials, chemicals, reagents and strains The S. yunnanensis plant samples were collected from Yunnan University of Traditional Chinese Medicine in Kunming City, Yunnan Province, China. Authentic standards of ( S )- N -methylcoclaurine and ( R )- N -methylcoclaurine were obtained through commercial chemical synthesis by WuXi LabNetwork. ( S )-norcoclaurine, ( R )-norcoclaurine, ( S )-coclaurine and ( R )-coclaurine were derived from the chiral separation of commercial racemic standard. The racemic norcoclaurine and coclaurine were purchased from Baoji Herbest Bio-Tech Co., Ltd. Other standards were purchased from Shanghai yuanye Bio-Technology Co., Ltd. ( R )-crotoflorine was prepared in previous work (Li et al., 2024 ). YPD medium was composed of 20 g L − 1 peptone (OXOID), 10 g L − 1 yeast extract (OXOID), and 20 g L − 1 glucose (Solarbio), serving as a standard medium for cultivating yeast and preparing competent cells. YPL medium, used to induce gene expression in yeast, uses galactose (Solarbio) instead of glucose compared to YPD medium. Synthetic dropout minus uracil medium (SD-Ura, FunGenome) with 20 g L − 1 glucose was employed to select positive colonies transformed with the pESC-Ura vector. For plates preparation, agar (DING GUO) at a concentration of 20 g L − 1 was added as necessary. WAT11 yeast strain with endogenous cytochrome P450 reductase replaced by one from Arabidopsis thaliana was maintained by this laboratory (Urban et al., 1997 ). pESC-Ura was purchased from Agilent Technologies. Escherichia coli Trans1 T1 strain from TransGen Biotech was used for routine plasmid assembly. 2.2 Alkaloid extraction and composition analysis 50 mg root, stem, and leaf freeze-dried powders of S. yunnanensis were mixed to 2 mL methanol and then sonicated extracts for 0.5 h. After centrifugation, the supernatants were filtered with a nylon syringe filter (0.22 µm). Quantitative analysis was conducted on a UPLC-QTOF-MS system (Waters Technologies). The Acquity UPLC was carried out using a T3 column (Waters Technologies, 2.1 mm × 100 mm, 2.7 µm particle size) at 38 ℃. The mobile phases consisted of eluent acetonitrile (A) and 0.1% aqueous formic acid (B) with a flow rate of 0.1 mL min − 1 with the linear gradient elution program: 5% − 12% A from 0.0 min to 1.0 min, 12% A from 1.0 min to 15.0 min, 12% − 17% A from 15.0 min to 16.0 min, 17% A from 16.0 min to 66.0 min, 17% ~ 90% A from 66.0 min to 67.0 min, 90% A from 67.0 min to 72.0 min, 90% − 5% A from 72.0 min to 72.5 min, and 5% A from 72.5 min to 76.0 min. The Acquity UPLC system was coupled to a Waters Xevo G2-S QTOF mass spectrometer equipped with electrospray ionization (ESI). The instrument was operated in positive ion mode to perform full scan monitoring in the range of m/z 50–800. The other operating parameters were set as follows: capillary voltage of 0.5 kV; sample cone voltage of 40 V; extraction cone voltage of 4 V; source temperature of 100 ℃; desolvation temperature 300 ℃; and desolvation gas flow of 800 L h − 1 . The trap collision energy of low energy function was set at 6 eV, while ramp trap collision energy of high energy function was set at 30–50 eV (Xiao et al., 2018 ). Data acquisition and processing were performed using MassLynx software. 2.3 Transcriptome analysis S. yunnanensis RNA was extracted using the R6827 Plant RNA Kit (Omega Bio-tek) according to the manufacturer's protocol. The RNA was quality-checked using Nanodrop, Qubit, and gel electrophoresis. Once the RNA samples were verified to be of acceptable quality, they were randomly fragmented using a Covaris sonicator. The entire library preparation process then involved end repair, A-tailing, adapter ligation, purification, and PCR amplification. Specifically, fragmented mRNA was used as a template with random oligonucleotides as primers to synthesize the first strand of cDNA in an M-MuLV reverse transcriptase system. The RNA strand was degraded with RNaseH, and the second strand of cDNA was synthesized using DNA polymerase I with dNTPs as substrates. The double-stranded cDNA underwent end repair, A-tailing, and sequencing adapter ligation. AMPure XP beads were used to select cDNA fragments of approximately 250–300 bp. After PCR amplification, the library was purified again using AMPure XP beads. Quality control of the library involved quantification with a Qubit 2.0 Fluorometer, followed by size verification with an Agilent 2100 bioanalyzer. Once the insert size met expectations, the effective library concentration was accurately quantified using qRT-PCR (with the effective library concentration exceeding 2 nM). The library preparation kit used was the NEBNext® Ultra™ RNA Library Prep Kit for Illumina®, and sequencing was performed on an Illumina Novaseq 6000 platform. To obtain comprehensive functional information of the transcripts, functional annotation was performed using seven databases, including Nr (Yangyang et al., 2006 ), Pfam (Finn et al., 2014 ), Uniprot (Apweiler et al., 2004 ), KEGG (Kanehisa et al., 2004 ), GO (Ashburner et al., 2000 ), KOG/COG (Tatusov et al., 2003 ), and PATHWAY (Kanehisa et al., 2004 ). Functional annotation primarily employed two methods: sequence similarity search and motif similarity search. For sequence similarity search, the protein sequences encoded by the transcripts were compared with existing protein databases, such as Uniprot, Nr (Yangyang et al., 2006 ), and the metabolic pathway database KEGG (Kanehisa et al., 2004 ), using diamond BLASTp (version: 2.0.6.144; parameters: –evalue 1e − 5 ) (Buchfink et al., 2015 ) to obtain functional information and potential metabolic pathway information. KEGG annotation was performed using KOBAS (version: 3.0) (Xie et al., 2011 ) to associate with KEGG ORTHOLOGY and PATHWAY. The Uniprot (Apweiler et al., 2004 ) database records the correspondence between each protein family and functional nodes in Gene Ontology (Ashburner et al., 2000 ), allowing prediction of the biological functions executed by the protein sequences encoded by the transcripts. Based on the associations between databases, KOG/COG (Tatusov et al., 2003 ) annotation results were obtained, followed by classification statistics and plotting for KOG/COG (Tatusov et al., 2003 ). For motif similarity search, proteins typically consist of one or more functional regions, commonly referred to as domains. The different combinations of these domains result in a variety of proteins; thus, identifying protein domains is crucial for analyzing protein functions. Domain prediction was performed using hmmscan (version: 3.3.2; parameters: e-value 0.01) (Finn et al., 2011 ) to obtain conserved sequences, motifs, and domains of the proteins. The Pfam (Finn et al., 2014 ) database is a large collection of protein families based on multiple sequence alignments and hidden Markov models. The expression level of transcripts was assessed by calculating their FPKM value. If the samples all contain biological replicates, used DESeq2 for differential expression analysis; otherwise, used edgeR. Herein, the differentially expressed genes (DEGs) were screened with the thresholds were a significance level of corrected value of padj 1. Perform GO enrichment and KEGG pathway enrichment analysis on DEGs, with the padj < 0.05, which was considered significantly enriched among DEGs, respectively. 2.4 Analysis of candidate genes in the BIA biosynthetic pathways To identify the candidate genes associated with BIA biosynthetic pathways of S. yunnanensis , local BLASTp by BioEdit, version 7.1.3.0 (Hall, 1999 ) was performed against S. yunnanensis sequences using query sequences of BIA-producing plants obtained from GenBank databases. The resulting transcripts with high degrees of identity were selected as candidate genes. A heatmap was generated using Tbtools with row scaling. Based on the reported genome of S. yunnanensis (Leng et al., 2024 ), we studied the chromosomal localization of these candidate genes and conducted synteny analysis of three species using TbTools. Based on the reference genes of BIA-producing plants and the obtained candidate genes, phylogenetic analysis was conducted using MEGA11 (Tamura et al., 2021 ) through the maximum likelihood estimation. 2.5 Cloning of candidate genes and eukaryotic expression of recombinant plasmids All candidate genes were amplified from S. yunnanensis root or leaf cDNA using the primers listed in the supplementary information (Table S4) and PrimeSTAR® Max DNA Polymerase (TaKaRa), according to the manufacturer's instructions. All CYP450 genes were purified and constructed into the pESC-Ura vector (digested with BamH Ⅰ) using the Gibson Assembly Kit (TransGen Biotech). The recombinant vectors were transformed into competent E. coli Trans1 T1 and identified by colony PCR (TransGen Biotech) and Sanger sequencing (Ruiotechnology). The identified positive recombinant plasmids were extracted from Trans1 T1 using the Plasmid Purification Kit (Magen Biotech). The recombinant plasmids pESC-Ura carrying candidate CYP450 genes were each transformed into WAT11 using the Frozen-EZ Yeast Transformation II Kit (Zymo Research). WAT11 transformed with empty pESC-Ura was employed as a control. The positive transformants were screened on SD-Ura medium containing 20 g L − 1 glucose. The positive clones were incubated with shaking at 30°C for 20–24 hours until the OD600 reached 2–3. Cells were centrifuged to remove the SD-Ura medium. The cells were then resuspended in YPL induction medium and grown overnight at 30°C to induce recombinant protein expression. 2.6 Microsomes extraction, enzymatic activity assay and LC-MS analysis Microsomes of recombinant yeast were prepared as previously described (Pompon et al., 1996 ), which has proven to be successful for previous such characterization. In vitro activity assays were performed in a 250 µL reaction system that included 100 mM Tris-HCl (pH 7.5, sangong) and 500 µM NADPH (Solarbio), 0.1 g microsomal protein, and 20 µM of the substrate. The reactions were incubated at 30 ℃ for 2 hours with 180 rpm shaking, then 5 µL ammonium hydroxide of was added to the reaction system to make the pH to about 10, and then extracted with 250 µL of ethyl acetate. After the ethyl acetate solution of the reaction product was concentrated to dryness, 150 µL of methanol was added to redissolve and centrifuged at 20,000 g for 15 min before LC-MS analysis. All enzymatic reaction products were detected using UPLC-QTOF-MS system (Waters Technologies). The Acquity UPLC was carried out as previously described. The linear gradient elution program: 5% − 30% A from 0.0 min to 6.0 min, 30% − 60% A from 6.0 min to 8.0 min, 60% − 90% A from 8.0 min to 8.5 min, 90% A from 8.5 min to 9.5 min, 90% − 5% A from 9.5 min to 10.0 min, and 5% A from 10 min to 11 min. Data acquisition and processing were performed using MassLynx software. 3. Results 3.1 Determination of BIAs in S. yunnanensis S. yunnanensis is a medicinal plant with a long history of use, with BIAs as its main active components (Dai et al., 2012 ; Hu et al., 2010 ). Due to its strict environmental requirements and long growth cycle, it had not yet been cultivated on a large scale. Recent studies have reported the presence of various aporphine alkaloids in S. yunnanensis (Wang et al., 2022 ) (Fig. S1). We analyzed three tissues of S. yunnanensis using liquid chromatography-mass spectrometry (LC-MS) and researched the relative contents of salutaridine, roemerine, stephanine, and crebanine to investigate the accumulation patterns of these four representative BIAs in S. yunnanensis (Fig. 1 , Fig. S2). Surprisingly, we detected four representative BIAs in the roots, each corresponding to the purchased reference standards, but in the stems and leaves, there were hardly any or only minimal amounts detected (Fig. S2). Furthermore, crebanine was the most abundant among the four compounds (Fig. 1 ), which is consistent with previous studies (Peng, 2014 ). Among these compounds, salutaridine was identified as a morphinan, another major compound reported in S. yunnanensis (Peng, 2014 ); Additionally, roemerine and stephanine were considered important intermediates, along with the representative product crebanine, all of which belong to type Ⅱ aporphines (Wang et al., 2022 ). We also attempted to detect additional aporphines that have been reported in other Stephania species (Peng, 2014 ; Wang et al., 2022 ). However, we did not detect magnoflorine, a representative type Ⅰ aporphine, suggesting that S. yunnanensis may contain more type Ⅱ aporphines, with type Ⅰ aporphines being less prevalent. This result might explain why S. yunnanensis was primarily used for its tuberous roots in medicine rather than its aerial parts. However, under the current extraction and detection conditions, we did not find any representative bisbenzylisoquinoline cepharanthine, which contradicts another study that reported the presence of cepharanthine and other bisbenzylisoquinolines in the metabolome of S. yunnanensis (Leng et al., 2024 ). In summary, the analysis of these compounds indicated that numerous structurally diverse BIAs are present in S. yunnanensis , specifically accumulating in the roots, with aporphines, represented by crebanine, being the most abundant. 3.2 Transcriptome sequencing, assembly, and analysis of S. yunnanensis . To explore the biosynthesis of aporphine alkaloids in S. yunnanensis , we prepared nine samples from the roots, stems, and leaves for transcriptome sequencing. RNA was extracted and cDNA libraries were constructed using Illumina and PacBio Sequel II for sequencing, resulting in 41,640,198 raw reads. De novo assembly yielded 50,119 transcripts, with an N50 length of 2,041 bp and an average length of 1,681 bp (Table 1 ). Table 1 Assembly results of de novo transcriptome of S. yunnanensis Item Number Seq Num 50,119 Seq Base (bp) 91,278,165 N50 (bp) 2,041 Max Length (bp) 7,046 Min Length (bp) 108 Average Length (bp) 1,821.23 Mean Length (bp) 1,681.00 Benchmarking Universal Single-Copy Orthologs (BUSCO) assessment (Manni et al., 2021 ) indicated that the S. yunnanensis transcriptome contained 87.3% complete BUSCOs (Fig. S3a, Table S1). Annotation was performed for all 50,119 (100%) transcripts across various databases, with an annotated transcript N50 length of 1,518 bp (Table S2). Gene Ontology (GO) classification (Ashburner et al., 2000 ) was conducted to describe gene functions in three categories: cellular component, molecular function, and biological process (Fig. S3b). To investigate the distribution of assembled transcripts in each tissue, their relative expression levels were determined by calculating the fragment per kilobase of transcript per million fragments mapped (FPKM) values of assembled transcripts. To compare the differentially expressed genes (DEGs) between roots and other tissues, we conducted a comparative analysis of the expression levels of assembled transcripts (Fig. 2 ). The results indicated that the number of DEGs between roots and leaves was the highest compared to other groups, with 1,780 up-regulated and 932 down-regulated (Fig. 2 a & 2 b). To further analyze the DEGs between roots and leaves, GO and KEGG clustering analyses were performed (Figs. 2 c & 2 d, Fig. S4). GO clustering revealed that most DEGs between roots and leaves were classified into “monooxygenase activity-related” categories, except for photosynthesis-related categories (Fig. 2 c). On the other hand, KEGG clustering showed that a total of 685 DEGs were considered to be involved in secondary metabolism, with only a small number (11 DEGs, with 4 down-regulated and 7 up-regulated) clustered into the “isoquinoline alkaloid biosynthesis” category (Fig. 2 d). 3.3 Analysis of BIA biosynthesis genes Apart from the two reported SyNCS4, SyNCS5 (Leng et al., 2024 ), the biosynthesis of BIAs in S. yunnanensis has not yet been elucidated. The upstream pathways, from dopamine and 4-HPAA to reticuline, were largely conserved among BIA-producing plants (Lee & Facchini, 2010 ). Analysis of the annotation results revealed a total of 141 OMT, 117 NMT, 45 NCSs, and 505 P450s (Tables S5-S8). In recent studies, CYP80s involved in BIA biosynthesis have gained increasing attention for their roles in the synthesis of 1-BIAs, as well as bisbenzylisoquinoline alkaloids and aporphine alkaloids (Hao et al., 2024 ; Li et al., 2024 ). Notably, among the P450s, 25 were further annotated into the CYP80 family, which are likely involved in BIA biosynthesis. To investigate the BIA biosynthesis in S. yunnanensis , we identified candidate genes with high similarity to previously reported BIA biosynthesis genes from S. tetrandra (Li et al., 2022 ; Li et al., 2020 ; Li et al., 2023 ). We selected sequences with homology greater than 55% and coverage greater than 90% as candidate genes, and a total of 26 candidate transcripts for the BIA biosynthesis were identified from the S. yunnanensis transcriptome. Phylogenetic analysis and domain analysis indicated that these candidate transcripts have structural similarities to the reference genes (Fig. 3 , Table S3). The sequence alignment revealed that these candidate genes share similar sequences with the reference genes (Fig. S5-S8). To analyze the expression patterns of these candidate transcripts in different tissues of S. yunnanensis , a heatmap was generated using the FPKM of the candidate genes (Fig. S9). The heatmap showed that, except for NCS, 4'OMT, and CYP80G6, most candidate transcripts did not exhibit distinct differential expression patterns and maintained high expression levels across multiple tissue parts. This may explain why only a few DEGs were clustered into the "isoquinoline alkaloid biosynthesis" category. The disparity between gene expression and compound accumulation suggests that BIA compounds transport might occur in S. yunnanensis . On the other hand, this disparity also suggests that the downstream genes of the BIA biosynthesis might not be highly expressed only in the roots, an important consideration for gene screening. Furthermore, based on the reported S. yunnanensis genome (Leng et al., 2024 ), we examined the distribution and location of these candidate transcripts on the chromosomes (Fig. 5 a). Sequence alignment revealed that 26 candidate transcripts corresponded to 18 genes in genome, they were arranged in descending order of homology and named according to their respective functional genes before being mapped onto the chromosomes. Notably, among these candidate genes, one SyNMCH, two SyCNMT and two Sy6OMT were located on pseudochromosome 3 in close proximity, suggesting that this region might form a BIA biosynthetic gene cluster (BGC). These were the genes involved in the three consecutive steps following NCS in the upstream BIA biosynthetic pathways. Similarly, another candidate gene SyCYP80G6 co-localized with two SyCNMT on pseudochromosome 11. Additionally, the candidate genes for CYP80Q5 were found on three different pseudochromosomes. Specifically, four were located on pseudochromosome 1, one was located on pseudochromosome 9, and another was co-localized with Sy4'OMT and SyNMCH2 on pseudochromosome 13. Among them, SyCYP80Q5-1 and SyCYP80Q5-2 shared over 98% identity, indicating they are tandem repeats, whereas they had 90% identity with SyCYP80Q5-3 (Fig. S10). The other three genes, SyCYP80Q5-4, SyCYP80Q5-5, and SyCYP80Q5-6, had low identity with other SyCYP80Q5 genes, approximately 59–65%. These chromosomal locations preliminarily suggested that the BIA biosynthetic genes in S. yunnanensis may exhibit a certain degree of clustered distribution. To investigate the distribution and association of these candidate genes in Stephania species, we conducted a synteny analysis using the reported genomes of two other plants, S. japonica and S. cepharantha (Leng et al., 2024 ) (Fig. 5 b). Interestingly, the two tandemly duplicated genes of SyCYP80Q5 on chromosome 1 of S. yunnanensis , SyCYP80Q5-1 and SyCYP80Q5-2 had only one ortholog in S. cepharantha and none in S. japonica , but SyCYP80Q5-3 was conserved across the three species. Notably, the clusters observed in S. yunnanensis , were present on different chromosomes in the three Stephania species. This indicated that the clustering patterns of BIA biosynthetic genes, were conserved across the three species, suggesting that such clustering is likely conserved within Stephania as well. In summary, chromosome localization and synteny analysis have provided insights into the genomic-level associations of BIA biosynthetic genes in S. yunnanensis . 3.4 Functional verification and phylogenetic analysis of CYP80s To elucidate the biosynthesis of the BIAs of S. yunnanensis , we conducted in vitro functional validation of the two key CYP80s involved in aporphine skeleton formation: CYP80G6 (CTS) and CYP80Q5. These enzymes specifically catalyzed the ( S ) or ( R )- types of 1-BIA substrates to produce the corresponding type Ⅰ aporphines or type Ⅱ aporphines (protoaporphines). The further rearrangement and modification of these intermediates ultimately led to the formation of aporphines such as crebanine. Since the high identity between SyCYP80Q5-1 and SyCYP80Q5-2, we cloned only SyCYP80Q5-1 and another four, SyCYP80Q5-3, SyCYP80Q5-4, SyCYP80Q5-5, SyCYP80Q5-6, as well as SyCYP80G6, and validated their functions of them in vitro . Recombinant plasmids containing these candidate genes were constructed and expressed in Saccharomyces cerevisiae (WAT11), in which the endogenous cytochrome P450 reductase was replaced by one from Arabidopsis thaliana . Microsomes extracted from the recombinant yeast were used for in vitro enzyme assays and were detected by LC-MS. We found that only SyCYP80Q5-1, SyCYP80Q5-3 and SyCYP80G6 exhibited catalytic functions (Fig. 6 , Figs. S11 & S12). Specifically, SyCYP80G6 catalyzes the conversion of ( S )-reticuline and ( S )-NMC into ( S )-corytuberine and ( S )-glaziovine (Figs. 6 a- 6 d, Figs. S11a & S11c). Both SyCYP80Q5-1 and SyCYP80Q5-3 were capable of catalyzing the conversion of ( R )-COC and ( R )-NMC into the aporphine alkaloids ( R )-crotoflorine and ( R )-glaziovine (Figs. 6 e- 6 h, Figs. S11b & S12). Furthermore, the catalytic efficiency of SyCYP80Q5-1 was slightly higher than that of SyCYP80Q5-3 when using either ( R )-COC or ( R )-NMC as substrates (Figs. 6 f & 6 h). Additionally, we discussed the activity of SyCYP80Q5-1 and SyCYP80G6 on a broader range of BIA substrates and other configurations (Fig. S13 & S14). As expected, SyCYP80Q5 and SyCYP80G6 exhibited strict substrate specificity and stereoselectivity, which was consistent with the results of other CYP80s (Li et al., 2024 ). We performed a phylogenetic analysis of the CYP80 family, and the results showed that the three SyCYP80s clustered within the CYP80 family (Fig. 7), and they were distributed across various clades of the CYP80 family. The functions of these enzymes matched the previously reported enzymes (Li et al., 2024 ), indicating that these three enzymes play similar roles in the BIA biosynthesis in S. yunnanensis . 4. Discussion The biosynthetic pathways of BIAs have not been fully elucidated, mainly due to the complex transformations from 1-BIAs to various types of BIAs scaffolds and the subsequent downstream modifications, including C-C coupling, C-O coupling, hydroxylation, and methylation, etc. (Li et al., 2024 ; Meng et al., 2024 ). This study reported the de novo transcriptome of S. yunnanensis , and through the combined analysis of the transcriptome, metabolome, and genome, identified 18 candidate genes involved in BIA biosynthesis in S. yunnanensis . The functions of three CYP80s were validated in vitro , they specifically catalyzed the C-C coupling of ( S )- or ( R )-configured 1-BIA substrates, resulting in the formation of aporphines or protoaporphines. This is an important step in the biosynthetic pathways of aporphines (Li et al., 2024 ). In medicinal plants where roots or rhizomes are used, the active components are often highly concentrated in the underground parts, with corresponding functional genes and regulatory factors showing significantly higher expression levels compared to the aboveground parts (He et al., 2018 ; Liu et al., 2017 ; Liu et al., 2023 ; Tong et al., 2022 ; Wang et al., 2015 ; Zhan et al., 2019 ). In this study, we analyzed the representative BIAs of S. yunnanensis and found that the major active component, crebanine and three other BIAs were highly accumulated in the roots, similar to other root-based medicinal plants. Notably, no upstream 1-BIAs were detected in the metabolomes of any tissues, possibly due to their complete consumption or concentrations below the detection limit. Furthermore, under the extraction and detection conditions of this study, the representative bisbenzylisoquinoline alkaloid, cepharanthine, was not detected. Then we sequenced the transcriptomes of S. yunnanensis using next-generation sequencing, followed by de novo assembly and functional annotation. Previous studies have shown that P450s are involved in almost the entire biosynthetic pathways of BIAs, specifically CYP80, which participates in hydroxylation and downstream C-O and C-C coupling reactions (An et al., 2024 ; Hao et al., 2024 ; Li et al., 2024 ; Meng et al., 2024 ). CYP719 catalyzed the formation of the methylenedioxy bridge (Hori et al., 2018 ; Ikezawa et al., 2003 ; Li et al., 2024 ; Menendez-Perdomo & Facchini, 2023 ), a characteristic group of many active BIA components, including crebanine. We identified 25 transcripts annotated as “CYP80” and 11 were annotated as “CYP719” in the de novo transcriptome of S. yunnanensis . Through annotation and sequence alignment, we identified 26 candidate transcripts involved in 10 steps of the BIA biosynthetic pathways. The heatmap showed that, except for a few genes, the majority were not specifically expressed in the roots but were highly expressed across multiple tissues. This pattern of metabolite accumulation and differential gene expression levels, while uncommon, is not unprecedented in BIA-producing plants. Similar mechanisms involving transporter proteins facilitating the movement of compounds within the plant have been reported previously in C. japonica (Sakai et al., 2002 ; Shitan et al., 2003 )d somniferum (Dastmalchi et al., 2019 ). We speculate that similar mechanisms may also exist in S. yunnanensis , facilitating the internal transfer of compounds through processes such as endocytosis and exocytosis. In summary, the transcriptome data obtained in this study are valuable for elucidating the biosynthetic pathways of BIAs, including crebanine. BGCs, composed of tightly arranged genes that collectively participate in the biosynthesis of specific metabolites, allow organisms to coordinate and efficiently regulate gene expression and synergistic metabolic reactions (Li et al., 2017 ). With the reported high-quality genome of S. yunnanensis (Leng et al., 2024 ), this study also examined the chromosomal localization of the 18 candidate genes and found a tendency for several candidate genes to cluster. Further collinearity analysis showed that this clustering phenomenon is conserved among different Stephania species, indicating important functions in their evolutionary processes. Notably, the tandem duplication of SyCYP80Q5-1 and SyCYP80Q5-2 on chromosome 1 of S. yunnanensis was not conserved among different species, whereas another SyCYP80Q5-3 gene on chromosome 13 was widely conserved. This might explain why all Stephania species can produce aporphines and protoaporphines, but differ in the types and amounts they contain. Furthermore, the number of CYP80Q5 copies in S. yunnanensis was significantly higher than CYP80G6, which was consistent with the finding that there are more type II aporphines present in S. yunnanensis . In conclusion, the differences and conservation of the BIA biosynthetic pathways in Stephania suggest that there are still variations within the genus. This may explain why there are significant differences in the types and contents of BIAs among different Stephania species. As an important gene family in the BIA biosynthetic pathways, CYP80, this study further validated the functions of three CYP80s from S. yunnanensis in vitro , showing that they all had the ability to catalyze substrate C-C coupling with configurational selectivity. Specifically, SyCYP80G6 specifically catalyzed the production of corresponding aporphines or protoaporphines from ( S )-type substrates; SyCYP80Q5-1 and SyCYP80Q5-3 specifically catalyzed the production of corresponding protoaporphines, the precursor of type Ⅱ aporphines, from ( R )-type substrates. The functions of CYP80 are complex and diverse, and there might be some evolutionary correlation. It is noteworthy that, although SyCYP80G6 can also catalyze ( S )-type substrates to produce aporphines and protoaporphines, similar to other CYP80Gs, we did not detect any representative reported ( S )-type aporphines, such as magnoflorine or mecambroline, in S. yunnanensis . However, the expression level of SyCYP80G6 is not low. This discrepancy between metabolites and gene expression suggests that there may be mechanisms in S. yunnanensis 's aporphine biosynthesis that we have yet to understand. It is possible that the metabolic flux is diverted towards (R)-reticuline, leading to the biosynthesis of ( R )-type aporphines, morphinans, and protoberberines. In conclusion, the combined approach of metabolite analysis and transcriptome sequencing identified 26 candidate transcripts responsible for the biosynthesis of BIAs, including crebanine, in S. yunnanensis , these genes might be responsible for BIAs biosynthesis in S. yunnanensis . Furthermore, genome analysis of the 18 candidate genes showed clustering on chromosomes, suggesting the presence of related BGCs, and collinearity analysis indicated these BGCs are conserved among different Stephania species. Finally, we identified three CYP80s from the S. yunnanensis transcriptome data as involved in BIAs biosynthesis using in vitro enzymatic reaction. Overall, our work provides valuable genetic information on S. yunnanensis and reveals the biosynthesis of BIAs in this medicinal plant. Declarations Competing Interests The authors have no relevant financial or non-financial interests to disclose. Funding This work was supported by the National Key R&D Program of China (2020YFA0908000), the National Natural Science Foundation of China (82011530137, 31961133007), Scientific and technological innovation project of CACMS (CI2023D002, CI2023E002) and Key project at central government level: The ability to establish sustainable use of valuable Chinese medicine resources (2060302). Author contributions LH and JGuo conceived and designed the entire research plans; WS and QL performed most of the experiments; LL, XL, JGan, YM, TC, and PS participated in some of the experiments; WS analyzed results of all the experiments; WS, JGuo, JW, YZ, XM, and LH wrote the manuscript; QL and YM revised the manuscript. All authors read and approved the manuscript. Data availability The raw data from the transcriptome sequencing have been uploaded to public database under the accession number XXXXXXXX. All other subsequent data are provided in the supplementary materials. References An Z, Gao R, Chen S, Tian Y, Li Q, Tian L, Zhang W, Kong L, Zheng B, Hao L, Xin T, Yao H, Wang Y, Song W, Hua X, Liu C, Song J, Fan H, Sun W, Chen S, Xu Z (2024) Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots. Adv Sci (Weinh) 11(19):e2309990. https://doi.org/10.1002/advs.202309990 Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LS (2004) UniProt: the Universal Protein knowledgebase. Nucleic Acids Res 32(Database issue) D115-119. https://doi.org/10.1093/nar/gkh131 Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25(1):25–29. https://doi.org/10.1038/75556 Buchfink B, Xie C, Huson DH (2015) Fast and sensitive protein alignment using DIAMOND. Nat Methods 12(1):59–60. https://doi.org/10.1038/nmeth.3176 Cui L, Peng C, Li J, Cheng X, Fan X, Li J, Yang Z, Zhao Y, Ma Y (2022) The anti-inflammatory and analgesic activities of 2Br-Crebanine and Stephanine from Stephania yunnanenses H. S.Lo. Front Pharmacol 13:1092583. https://doi.org/10.3389/fphar.2022.1092583 Dai X, Hu R, Sun C, Pan Y (2012) Comprehensive separation and analysis of alkaloids from Stephania yunnanensis by counter-current chromatography coupled with liquid chromatography tandem mass spectrometry analysis. J Chromatogr A 1226:18–23. https://doi.org/10.1016/j.chroma.2011.10.022 Dastmalchi M, Chang L, Chen R, Yu L, Chen X, Hagel JM, Facchini PJ (2019) Purine Permease-Type Benzylisoquinoline Alkaloid Transporters in Opium Poppy. Plant Physiol 181(3):916–933. https://doi.org/10.1104/pp.19.00565 Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer EL, Tate J, Punta M (2014) Pfam: the protein families database. Nucleic Acids Res 42(Database issue):D222–230. https://doi.org/10.1093/nar/gkt1223 Finn RD, Clements J, Eddy SR (2011) HMMER web server: interactive sequence similarity searching. Nucleic Acids Res W29–37, 39(Web Server issue). https://doi.org/10.1093/nar/gkr367 Ghirga F, Bonamore A, Calisti L, D'Acquarica I, Mori M, Botta B, Boffi A, Macone A (2017) Green Routes for the Production of Enantiopure Benzylisoquinoline Alkaloids. Int J Mol Sci 18(11). https://doi.org/10.3390/ijms18112464 Gurkok T, Ozhuner E, Parmaksiz I, Ozcan S, Turktas M, Ipek A, Demirtas I, Okay S, Unver T (2016) Functional Characterization of 4'OMT and 7OMT Genes in BIA Biosynthesis. Front Plant Sci 7:98. https://doi.org/10.3389/fpls.2016.00098 Hagel JM, Morris JS, Lee EJ, Desgagne-Penix I, Bross CD, Chang L, Chen X, Farrow SC, Zhang Y, Soh J, Sensen CW, Facchini PJ (2015) Transcriptome analysis of 20 taxonomically related benzylisoquinoline alkaloid-producing plants. BMC Plant Biol 15:227. https://doi.org/10.1186/s12870-015-0596-0 Hall TA (1999) BIOEDIT: A USER-FRIENDLY BIOLOGICAL SEQUENCE ALIGNMENT EDITOR AND ANALYSIS PROGRAM FOR WINDOWS 95/98/ NT Hao C, Yu Y, Liu Y, Liu A, Chen S (2024) The CYP80A and CYP80G Are Involved in the Biosynthesis of Benzylisoquinoline Alkaloids in the Sacred Lotus (Nelumbo nucifera). Int J Mol Sci 25(2). https://doi.org/10.3390/ijms25020702 He SM, Liang YL, Cong K, Chen G, Zhao X, Zhao QM, Zhang JJ, Wang X, Dong Y, Yang JL, Zhang GH, Qian ZL, Fan W, Yang SC (2018) Identification and Characterization of Genes Involved in Benzylisoquinoline Alkaloid Biosynthesis in Coptis Species. Front Plant Sci 9:731. https://doi.org/10.3389/fpls.2018.00731 Heming Y, Xianrui L (1980) BOTANICAL AND PHARMACOGNOSTICAL STUDIES OF THE CHINESE DRUG SHAN-WU-GUI. Acta Pharm Sinica (11), 674–683 Hori K, Yamada Y, Purwanto R, Minakuchi Y, Toyoda A, Hirakawa H, Sato F (2018) Mining of the Uncharacterized Cytochrome P450 Genes Involved in Alkaloid Biosynthesis in California Poppy Using a Draft Genome Sequence. Plant Cell Physiol 59(2):222–233. https://doi.org/10.1093/pcp/pcx210 Hu R, Dai X, Lu Y, Pan Y (2010) Preparative separation of isoquinoline alkaloids from Stephania yunnanensis by pH-zone-refining counter-current chromatography. J Chromatogr B Analyt Technol Biomed Life Sci 878(21):1881–1884. https://doi.org/10.1016/j.jchromb.2010.05.005 Ikezawa N, Iwasa K, Sato F (2008) Molecular cloning and characterization of CYP80G2, a cytochrome P450 that catalyzes an intramolecular C-C phenol coupling of (S)-reticuline in magnoflorine biosynthesis, from cultured Coptis japonica cells. J Biol Chem 283(14):8810–8821. https://doi.org/10.1074/jbc.M705082200 Ikezawa N, Tanaka M, Nagayoshi M, Shinkyo R, Sakaki T, Inouye K, Sato F (2003) Molecular cloning and characterization of CYP719, a methylenedioxy bridge-forming enzyme that belongs to a novel P450 family, from cultured Coptis japonica cells. J Biol Chem 278(40):38557–38565. https://doi.org/10.1074/jbc.M302470200 Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M (2004) The KEGG resource for deciphering the genome. Nucleic Acids Res 32(Database issue). https://doi.org/10.1093/nar/gkh063 . D277-280 Lai Y, Ma J, Zhang X, Xuan X, Zhu F, Ding S, Shang F, Chen Y, Zhao B, Lan C, Unver T, Huo G, Li X, Wang Y, Liu Y, Lu M, Pan X, Yang D, Li M, Zhang B, Zhang D (2024) High-quality chromosome-level genome assembly and multi-omics analysis of rosemary (Salvia rosmarinus) reveals new insights into the environmental and genome adaptation. Plant Biotechnol J 22(7):1833–1847. https://doi.org/10.1111/pbi.14305 Lee EJ, Facchini P (2010) Norcoclaurine synthase is a member of the pathogenesis-related 10/Bet v1 protein family. Plant Cell 22(10):3489–3503. https://doi.org/10.1105/tpc.110.077958 Leng L, Xu Z, Hong B, Zhao B, Tian Y, Wang C, Yang L, Zou Z, Li L, Liu K, Peng W, Liu J, An Z, Wang Y, Duan B, Hu Z, Zheng C, Zhang S, Li X, Li M, Liu Z, Bi Z, He T, Liu B, Fan H, Song C, Tong Y, Chen S (2024) Cepharanthine analogs mining and genomes of Stephania accelerate anti-coronavirus drug discovery. Nat Commun 15(1):1537. https://doi.org/10.1038/s41467-024-45690-5 Li K, Chen X, Zhang J, Wang C, Xu Q, Hu J, Kai G, Feng Y (2022) Transcriptome Analysis of Stephania tetrandra and Characterization of Norcoclaurine-6-O-Methyltransferase Involved in Benzylisoquinoline Alkaloid Biosynthesis. Front Plant Sci 13:874583. https://doi.org/10.3389/fpls.2022.874583 Li L, Jiang W, Lu Y (2017) New strategies and approaches for engineering biosynthetic gene clusters of microbial natural products. Biotechnol Adv 35(8):936–949. https://doi.org/10.1016/j.biotechadv.2017.03.007 Li Q, Bu J, Ma Y, Yang J, Hu Z, Lai C, Xu Y, Tang J, Cui G, Wang Y, Zhao Y, Jin B, Shen Y, Guo J, Huang L (2020) Characterization of O-methyltransferases involved in the biosynthesis of tetrandrine in Stephania tetrandra. J Plant Physiol 250:153181. https://doi.org/10.1016/j.jplph.2020.153181 Li Q, Jiao X, Li X, Shi W, Ma Y, Tan X, Gan J, Liu J, Yang J, Wang J, Jin B, Chen T, Su P, Zhao Y, Zhang Y, Tang J, Cui G, Chen Y, Guo J, Huang L (2024) Identification of the cytochrome P450s responsible for the biosynthesis of two types of aporphine alkaloids and their de novo biosynthesis in yeast. J Integr Plant Biol. https://doi.org/10.1111/jipb.13724 Li X, Li Q, Jiao X, Tang H, Cheng Y, Ma Y, Cui G, Tang J, Chen Y, Guo J, Huang L (2023) Phylogenetic analysis and functional characterization of norcoclaurine synthase involved in benzylisoquinoline alkaloids biosynthesis in Stephania tetrandra. J Cell Physiol. https://doi.org/10.1002/jcp.31065 Liu S, Liu M, Wang S, Lin Y, Zhang H, Wang Q, Zhao Y (2017) Analysis of the Panax ginseng stem/leaf transcriptome and gene expression during the leaf expansion period. Mol Med Rep 16(5):6396–6404. https://doi.org/10.3892/mmr.2017.7377 Liu X, Bu J, Ma Y, Chen Y, Li Q, Jiao X, Hu Z, Cui G, Tang J, Guo J, Huang L (2021) Functional characterization of (S)-N-methylcoclaurine 3'-hydroxylase (NMCH) involved in the biosynthesis of benzylisoquinoline alkaloids in Corydalis yanhusuo. Plant Physiol Biochem 168:507–515. https://doi.org/10.1016/j.plaphy.2021.09.042 Liu X, Shao P, Wang Y, Chen Y, Cui S (2023) Anti-inflammatory mechanism of the optimized active ingredients of Sargentodoxa cuneata and Patrinia villosa. Int Immunopharmacol 120:110337. https://doi.org/10.1016/j.intimp.2023.110337 Manni M, Berkeley MR, Seppey M, Simao FA, Zdobnov EM (2021) BUSCO Update: Novel and Streamlined Workflows along with Broader and Deeper Phylogenetic Coverage for Scoring of Eukaryotic, Prokaryotic, and Viral Genomes. Mol Biol Evol 38(10):4647–4654. https://doi.org/10.1093/molbev/msab199 Menendez-Perdomo IM, Facchini PJ (2020) Isolation and characterization of two O-methyltransferases involved in benzylisoquinoline alkaloid biosynthesis in sacred lotus (Nelumbo nucifera). J Biol Chem 295(6):1598–1612. https://doi.org/10.1074/jbc.RA119.011547 Menendez-Perdomo IM, Facchini PJ (2023) Elucidation of the (R)-enantiospecific benzylisoquinoline alkaloid biosynthetic pathways in sacred lotus (Nelumbo nucifera). Sci Rep 13(1):2955. https://doi.org/10.1038/s41598-023-29415-0 Meng F, Zhang S, Su J, Zhu B, Pan X, Qiu X, Cui X, Wang C, Niu L, Li C, Lu S (2024) Characterization of two CYP80 enzymes provides insights into aporphine alkaloid skeleton formation in Aristolochia contorta. Plant J 118(5):1439–1454. https://doi.org/10.1111/tpj.16686 Park S, An B, Park S (2024) Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae). BMC Plant Biol 24(1):303. https://doi.org/10.1186/s12870-024-05025-4 Pauli HH, Kutchan TM (1998) Molecular cloning and functional heterologous expression of two alleles encoding (S)-N-methylcoclaurine 3'-hydroxylase (CYP80B1), a new methyl jasmonate-inducible cytochrome P-450-dependent mono-oxygenase of benzylisoquinoline alkaloid biosynthesis. Plant J 13(6):793–801. https://doi.org/10.1046/j.1365-313x.1998.00085.x Pei L, Wang B, Ye J, Hu X, Fu L, Li K, Ni Z, Wang Z, Wei Y, Shi L, Zhang Y, Bai X, Jiang M, Wang S, Ma C, Li S, Liu K, Li W, Cong B (2021) Genome and transcriptome of Papaver somniferum Chinese landrace CHM indicates that massive genome expansion contributes to high benzylisoquinoline alkaloid biosynthesis. Hortic Res 8(1):5. https://doi.org/10.1038/s41438-020-00435-5 Peng F (2014) Chemical constituents from Stephania yunnannensis [master, Yunnan University of Chinese Medicine] Pompon D, Louerat B, Bronine A, Urban P (1996) Yeast expression of animal and plant P450s in optimized redox environments. Methods Enzymol 272:51–64. https://doi.org/10.1016/s0076-6879(96)72008-6 Pyne ME, Gold ND, Martin VJJ (2023) Pathway elucidation and microbial synthesis of proaporphine and bis-benzylisoquinoline alkaloids from sacred lotus (Nelumbo nucifera). Metab Eng 77:162–173. https://doi.org/10.1016/j.ymben.2023.03.010 Sakai K, Shitan N, Sato F, Ueda K, Yazaki K (2002) Characterization of berberine transport into Coptis japonica cells and the involvement of ABC protein. J Exp Bot 53(376):1879–1886. https://doi.org/10.1093/jxb/erf052 Shitan N, Bazin I, Dan K, Obata K, Kigawa K, Ueda K, Sato F, Forestier C, Yazaki K (2003) Involvement of CjMDR1, a plant multidrug-resistance-type ATP-binding cassette protein, in alkaloid transport in Coptis japonica. Proc Natl Acad Sci U S A 100(2):751–756. https://doi.org/10.1073/pnas.0134257100 Tamura K, Stecher G, Kumar S (2021) MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol Biol Evol 38(7):3022–3027. https://doi.org/10.1093/molbev/msab120 Tan J, Xiang Y, Xiong Y, Zhang Y, Qiao B, Zhang H (2023) Crebanine induces ROS-dependent apoptosis in human hepatocellular carcinoma cells via the AKT/FoxO3a signaling pathway. Front Pharmacol 14:1069093. https://doi.org/10.3389/fphar.2023.1069093 Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA (2003) The COG database: an updated version includes eukaryotes. BMC Bioinformatics 4:41. https://doi.org/10.1186/1471-2105-4-41 Tong Q, Zhang C, Tu Y, Chen J, Li Q, Zeng Z, Wang F, Sun L, Huang D, Li M, Qiu S, Chen W (2022) Biosynthesis-based spatial metabolome of Salvia miltiorrhiza Bunge by combining metabolomics approaches with mass spectrometry-imaging. Talanta 238(Pt 2):123045. https://doi.org/10.1016/j.talanta.2021.123045 Urban P, Mignotte C, Kazmaier M, Delorme F, Pompon D (1997) Cloning, yeast expression, and characterization of the coupling of two distantly related Arabidopsis thaliana NADPH-cytochrome P450 reductases with P450 CYP73A5. J Biol Chem 272(31):19176–19186. https://doi.org/10.1074/jbc.272.31.19176 Wang H, Cheng X, Kong S, Yang Z, Wang H, Huang Q, Li J, Chen C, Ma Y (2016) Synthesis and Structure-Activity Relationships of a Series of Aporphine Derivatives with Antiarrhythmic Activities and Acute Toxicity. Molecules 21(12). https://doi.org/10.3390/molecules21121555 Wang K, Jiang S, Sun C, Lin Y, Yin R, Wang Y, Zhang M (2015) The Spatial and Temporal Transcriptomic Landscapes of Ginseng, Panax ginseng C. A. Meyer. Sci Rep 5:18283. https://doi.org/10.1038/srep18283 Wang M, Zhang XM, Fu X, Zhang P, Hu WJ, Yang BY, Kuang HX (2022) Alkaloids in genus stephania (Menispermaceae): A comprehensive review of its ethnopharmacology, phytochemistry, pharmacology and toxicology. J Ethnopharmacol 293:115248. https://doi.org/10.1016/j.jep.2022.115248 Xiao J, Song N, Lu T, Pan Y, Song J, Chen G, Sun L, Li N (2018) Rapid characterization of TCM Qianjinteng by UPLC-QTOF-MS and its application in the evaluation of three species of Stephania. J Pharm Biomed Anal 156:284–296. https://doi.org/10.1016/j.jpba.2018.04.044 Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L (2011) KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res W316–322, 39(Web Server issue). https://doi.org/10.1093/nar/gkr483 Xu T, Yang X, Jia Y, Li Z, Tang G, Li X, Wang B, Wang T, Lin J, Guo L, Ye K (2022) A global survey of the transcriptome of the opium poppy (Papaver somniferum) based on single-molecule long-read isoform sequencing. Plant J 110(2):607–620. https://doi.org/10.1111/tpj.15689 Yang Y, Hao T, Yao X, Che Y, Liu Y, Fang M, Wang Y, Zhou D, Chai H, Li N, Hou Y (2023) Crebanine ameliorates ischemia-reperfusion brain damage by inhibiting oxidative stress and neuroinflammation mediated by NADPH oxidase 2 in microglia. Phytomedicine 120:155044. https://doi.org/10.1016/j.phymed.2023.155044 Yangyang D, Jianqi LI, Songfeng WU, Yunping Z, Yaowen C, Fuchu HE (2006) Integrated nr Database in Protein Annotation System and Its Localization. Comput Eng 32(5):71–72 Yeh PS, Liu CT, Yu CY, Chang YC, Lin SY, Li YC, Luan YZ, Sung WW (2024) Crebanine, an aporphine alkaloid, induces cancer cell apoptosis through PI3K-Akt pathway in glioblastoma multiforme. Front Pharmacol 15:1419044. https://doi.org/10.3389/fphar.2024.1419044 Zhan Z, Fang W, Ma X, Chen T, Cui G, Ma Y, Kang L, Nan T, Lin H, Tang J, Zhang Y, Lai C, Ren Z, Wang Y, Zhao Y, Shen Y, Wang L, Zeng W, Guo J, Huang L (2019) Metabolome and transcriptome analyses reveal quality change in the orange-rooted Salvia miltiorrhiza (Danshen) from cultivated field. Chin Med 14:42. https://doi.org/10.1186/s13020-019-0265-6 Zhao Y, Cui L, Yang XX, Sun X, Liu Y, Yang Z, Zhu L, Peng C, Li D, Cai J, Ma Y (2021) Sinoacutine inhibits inflammatory responses to attenuates acute lung injury by regulating NF-kappaB and JNK signaling pathways. BMC Complement Med Ther 21(1):284. https://doi.org/10.1186/s12906-021-03458-0 Zhao Y, Zhang Z, Li M, Luo J, Chen F, Gong Y, Li Y, Wei Y, Su Y, Kong L (2019) Transcriptomic profiles of 33 opium poppy samples in different tissues, growth phases, and cultivars. Sci Data 6(1):66. https://doi.org/10.1038/s41597-019-0082-x Additional Declarations The authors declare no competing interests. Supplementary Files Supplementaryfile1.docx Supplementary file 1 Supplementaryfile2.xlsx Supplementary file 2 Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-5384973","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":373641676,"identity":"966614f1-eeef-4ad6-8f82-db2556078275","order_by":0,"name":"Wenlong Shi","email":"","orcid":"https://orcid.org/0009-0003-6333-7137","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Wenlong","middleName":"","lastName":"Shi","suffix":""},{"id":373643015,"identity":"adc1a8fa-33de-4671-87d3-ab7a0f7dd818","order_by":1,"name":"Qishuang Li","email":"","orcid":"https://orcid.org/0000-0002-2148-7917","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Qishuang","middleName":"","lastName":"Li","suffix":""},{"id":373643016,"identity":"9b85c366-e90d-4b80-abe9-26a83872174c","order_by":2,"name":"Xinyi Li","email":"","orcid":"https://orcid.org/0009-0003-1144-5363","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Xinyi","middleName":"","lastName":"Li","suffix":""},{"id":373643017,"identity":"ef59d7e6-3ef6-4ff5-b5ec-13f066194dc6","order_by":3,"name":"Jingyi Gan","email":"","orcid":"https://orcid.org/0009-0004-9261-6159","institution":"College of Chinese Materia Medica, Yunnan Key Laboratory of Southern Medicinal Utilization, Yunnan University of Chinese Medicine","correspondingAuthor":false,"prefix":"","firstName":"Jingyi","middleName":"","lastName":"Gan","suffix":""},{"id":373643018,"identity":"334f3e81-a5af-4d4a-a5d9-99df69d09da1","order_by":4,"name":"Ying Ma","email":"","orcid":"https://orcid.org/0009-0005-7484-1720","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Ying","middleName":"","lastName":"Ma","suffix":""},{"id":373643019,"identity":"fd0f98d1-3301-4822-9435-ed82d0115ebd","order_by":5,"name":"Jian Wang","email":"","orcid":"https://orcid.org/0000-0002-9386-7986","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Jian","middleName":"","lastName":"Wang","suffix":""},{"id":373643020,"identity":"3bcf5c08-f7bb-4e0c-8c1e-bd3e86fe2680","order_by":6,"name":"Tong Chen","email":"","orcid":"https://orcid.org/0000-0003-3134-3113","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Tong","middleName":"","lastName":"Chen","suffix":""},{"id":373643021,"identity":"b45afb4d-4878-4225-80e4-c4283d03fe81","order_by":7,"name":"Yifeng Zhang","email":"","orcid":"https://orcid.org/0000-0002-6380-4431","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Yifeng","middleName":"","lastName":"Zhang","suffix":""},{"id":373643022,"identity":"6f2ae9f0-1e2a-45de-b614-fd5a8543abf0","order_by":8,"name":"Ping Su","email":"","orcid":"https://orcid.org/0000-0002-5187-4889","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Ping","middleName":"","lastName":"Su","suffix":""},{"id":373643023,"identity":"d1362569-e3b0-4d4b-b8a9-3b76ea283d1d","order_by":9,"name":"Xiaohui Ma","email":"","orcid":"https://orcid.org/0000-0001-5961-3168","institution":"College of Chinese Materia Medica, Yunnan Key Laboratory of Southern Medicinal Utilization, Yunnan University of Chinese Medicine","correspondingAuthor":false,"prefix":"","firstName":"Xiaohui","middleName":"","lastName":"Ma","suffix":""},{"id":373643024,"identity":"06b01c55-0975-441d-8522-a45f9a3e0cbc","order_by":10,"name":"Juan Guo","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAyklEQVRIiWNgGAWjYJACAzDJ3tj48ANpWngONxtLkGaXRHqbAA8xCuUjkg8U89TcsTe4+bCNQYLBTk63gYAWwxtpCcY8x54xG9xObHtQwJBsbHaAkJYZOQbGPGyH2YBa2g0kGA4kbiNOy7/DPAY3D7ZJ8BCjRV4CqIW37bCEwQ1GIrUY8DxLMJzbd9hA8kwiMJANiPCLfHvyMYM33w7b8x0//vDhhwo7OYJaDA4wsIGjUgGs0oCAcrAtDQzMD6CMUTAKRsEoGAXYAQDTWUM3SkdbkQAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0000-0002-1327-4428","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":true,"prefix":"","firstName":"Juan","middleName":"","lastName":"Guo","suffix":""},{"id":373643025,"identity":"a8fb0423-e232-491a-ab44-2b0716221973","order_by":11,"name":"Luqi Huang","email":"","orcid":"https://orcid.org/0000-0002-2070-4318","institution":"State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences","correspondingAuthor":false,"prefix":"","firstName":"Luqi","middleName":"","lastName":"Huang","suffix":""}],"badges":[],"createdAt":"2024-11-04 05:34:09","currentVersionCode":1,"declarations":{"humanSubjects":false,"vertebrateSubjects":false,"conflictsOfInterestStatement":false,"humanSubjectEthicalGuidelines":false,"humanSubjectConsent":false,"humanSubjectClinicalTrial":false,"humanSubjectCaseReport":false,"vertebrateSubjectEthicalGuidelines":false},"doi":"10.21203/rs.3.rs-5384973/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5384973/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":68251499,"identity":"69de907e-a642-4163-bd88-9291f35c1e38","added_by":"auto","created_at":"2024-11-05 10:02:04","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":199551,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003ea \u003c/strong\u003eDetermination of BIAs in \u003cem\u003eS. yunnanensis\u003c/em\u003e.\u003cstrong\u003e \u003c/strong\u003eLC-MS chromatogram of \u003cem\u003eS. yunnanensis\u003c/em\u003e roots, with peaks of four compounds labeled with corresponding numbers. The concentration of all four chemical standards was 0.02 mM. \u003cstrong\u003eb\u003c/strong\u003e A bar chart of relative content based on the average peak area, with four compounds represented by different colors. \u003cstrong\u003ec\u003c/strong\u003e The structures of the four main BIAs\u003c/p\u003e","description":"","filename":"image1.png","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/c03513cd8b0d745cb7fec3d3.png"},{"id":68252214,"identity":"b5bd5a58-97fc-404b-add0-255cd8d0304f","added_by":"auto","created_at":"2024-11-05 10:10:05","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":360702,"visible":true,"origin":"","legend":"\u003cp\u003eFunctional annotation of DEGs. \u003cstrong\u003ea\u003c/strong\u003e DEGs of root vs. leaf, root vs. stem, and stem vs. leaf in \u003cem\u003eS. yunnanensis\u003c/em\u003e, and the Venn diagram (\u003cstrong\u003eb\u003c/strong\u003e). \u003cstrong\u003ec\u003c/strong\u003e GO enrichment analysis of DEGs in root vs. leaf. \u003cstrong\u003ed\u003c/strong\u003e Scatterplot of KEGG pathway enrichment of DEGs in root vs. leaf\u003c/p\u003e","description":"","filename":"image2.png","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/34e245ef4a72ed54f3d3278a.png"},{"id":68251502,"identity":"3384fa68-e37a-4607-bba4-34c4ef5826a5","added_by":"auto","created_at":"2024-11-05 10:02:05","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":500576,"visible":true,"origin":"","legend":"\u003cp\u003ePhylogenetic and domain analysis of genes involved in the BIA biosynthesis. Phylogenetic tree of functional genes and candidate transcripts from \u003cem\u003eS. yunnanensis\u003c/em\u003e, constructed using the maximum likelihood method, and protein domain analysis\u003c/p\u003e","description":"","filename":"image3.png","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/b1462c292e3f7eb041537bdc.png"},{"id":68251501,"identity":"f2df7b0a-3182-455d-8b14-98a20a67ced5","added_by":"auto","created_at":"2024-11-05 10:02:05","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":381725,"visible":true,"origin":"","legend":"\u003cp\u003eChromosomal localization of \u003cem\u003eS. yunnanensis\u003c/em\u003e gene and synteny analysis of the genomes of three Stephania species.\u003cstrong\u003e a\u003c/strong\u003eChromosomal localization of candidate genes for BIA biosynthesis in \u003cem\u003eS. yunnanensis\u003c/em\u003e. Pseudomolecules are represented as long bars, with colors indicating gene density: blue for low density and red for high density. \u003cstrong\u003eb\u003c/strong\u003eSynteny analysis of the genomes of three Stephania species, with candidate genes in \u003cem\u003eS. yunnanensis\u003c/em\u003e highlighted\u003c/p\u003e","description":"","filename":"image4.png","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/ed471d144592ab3cf1c33434.png"},{"id":68251506,"identity":"4fc05d60-43a6-4c37-876a-e889bb74f503","added_by":"auto","created_at":"2024-11-05 10:02:05","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":411635,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eIn vitro\u003c/em\u003e functional validation of the CYP80s from \u003cem\u003eS. yunnanensis\u003c/em\u003e. SyCYP80G6 was responsible for catalyzing (\u003cem\u003eS\u003c/em\u003e)-type substrates, SyCYP80Q5 was responsible for catalyzing (\u003cem\u003eR\u003c/em\u003e)-type substrates. The color of products' names corresponds to the color of chromatographic peaks. \u003cstrong\u003ea, b\u003c/strong\u003e SYCYP80G6 catalyzes the conversion of (\u003cem\u003eS\u003c/em\u003e)-reticuline to (\u003cem\u003eS\u003c/em\u003e)-corytuberine. \u003cstrong\u003ec, d\u003c/strong\u003e SyCYP80G6 catalyzes the conversion of (\u003cem\u003eS\u003c/em\u003e)-NMC to (\u003cem\u003eS\u003c/em\u003e)-glaziovine. \u003cstrong\u003ee, f\u003c/strong\u003e SyCYP80Q5-1 and SyCYP80Q5-3 catalyzes the conversion of (\u003cem\u003eR\u003c/em\u003e)-NMC to (\u003cem\u003eR\u003c/em\u003e)-glaziovine. \u003cstrong\u003eg, h\u003c/strong\u003e SyCYP80Q5-1 and SyCYP80Q5-3 catalyzes the conversion of (\u003cem\u003eR\u003c/em\u003e)-COC to (\u003cem\u003eR\u003c/em\u003e)-crotoflorine\u003c/p\u003e","description":"","filename":"image5.png","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/c7924f405171b854a8e9e09a.png"},{"id":68251504,"identity":"8fab35a6-5ada-4f55-aca7-91e9006fc44f","added_by":"auto","created_at":"2024-11-05 10:02:05","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":920683,"visible":true,"origin":"","legend":"\u003cp\u003ePhylogenetic tree of CYP80s.\u003cstrong\u003e \u003c/strong\u003eThe red font represent P450s identified in this study\u003c/p\u003e","description":"","filename":"image6.png","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/3a86d38a1ea0ea7f572dbec8.png"},{"id":68252310,"identity":"ba5d7517-c897-4db3-b1e8-35e0aa92f6cc","added_by":"auto","created_at":"2024-11-05 10:18:08","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2752298,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/95dd06d8-3264-4275-87e0-c188c120c77d.pdf"},{"id":68252215,"identity":"78dca138-ca11-45fb-9998-2eb8175356b0","added_by":"auto","created_at":"2024-11-05 10:10:05","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":7827832,"visible":true,"origin":"","legend":"\u003cp\u003eSupplementary file 1\u003c/p\u003e","description":"","filename":"Supplementaryfile1.docx","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/3ca15cb3fa1828fc0162ce41.docx"},{"id":68252213,"identity":"32b8bb70-6632-456d-9e05-14d6b7911bef","added_by":"auto","created_at":"2024-11-05 10:10:05","extension":"xlsx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":79396,"visible":true,"origin":"","legend":"\u003cp\u003eSupplementary file 2\u003c/p\u003e","description":"","filename":"Supplementaryfile2.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-5384973/v1/89b4afa7f104872e433c5d9a.xlsx"}],"financialInterests":"The authors declare no competing interests.","formattedTitle":"\u003cp\u003eTranscriptome analysis of \u003cem\u003eStephania yunnanensis\u003c/em\u003e and functional validation of CYP80s involved in benzylisoquinoline alkaloids biosynthesis\u003c/p\u003e","fulltext":[{"header":"Key Message","content":"\u003cp\u003eThis study found that BIA biosynthetic genes in \u003cem\u003eStephania yunnanensis\u003c/em\u003e are not all highly expressed in roots and are clustered on chromosomes. Two CYP80s catalyze key steps in aporphine biosynthesis with strict conformational selectivity.\u003c/p\u003e"},{"header":"1. Introduction","content":"\u003cp\u003e \u003cem\u003eStephania yunnanensis\u003c/em\u003e is one of the source plants of the traditional Chinese medicine \u0026lsquo;Shanwugui\u0026rsquo;, which is a plant from the genus \u003cem\u003eStephania\u003c/em\u003e in the family \u003cem\u003eMenispermaceae\u003c/em\u003e of the order \u003cem\u003eRanales\u003c/em\u003e (Heming \u0026amp; Xianrui, \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e1980\u003c/span\u003e). Traditional Chinese medicine theory believes that \u003cem\u003eS. yunnanensis\u003c/em\u003e has the effects of vomiting phlegm and foods, treating malaria, and detoxifying sores (Zhao et al., \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Previous chemical research has shown that the main active components of \u003cem\u003eS. yunnanensis\u003c/em\u003e were four types of benzylisoquinoline alkaloids (BIAs), including 1-benzylisoquinolines (1-BIAs), aporphines, morphinans, and protoberberines (Dai et al., \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Hu et al., \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Wang et al., \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), such as salutaridine, tetrahydropalmatine, roemerine, stephanine, and crebanine, etc. Crebanine, a kind of aporphine, was major active compound found in \u003cem\u003eS. yunnanensis\u003c/em\u003e (Peng, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), recent studies have indicated that crebanine has neuroprotective (Yang et al., \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), cardioprotective(Wang et al., \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), anticancer (Tan et al., \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Yeh et al., \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2024\u003c/span\u003e), anti-inflammatory and analgesic activities (Cui et al., \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eSince the extraction of BIAs from plants was an expensive process and the content was also affected by the growth conditions of the plants, this low yield and variability limited the commercial development of BIAs. The complexity of chemical synthesis further hindered the large-scale production of BIAs. The demand for BIAs ultimately promoted the development of BIAs metabolic engineering and synthetic biology. The biosynthetic pathways of BIAs were initiated with the enzymatic reaction of norcoclaurine synthase (NCS), which catalyzes the synthesis of norcoclaurine from dopamine and 4-hydroxyphenylacetic acid (4-HPAA) (Ghirga et al., \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Subsequent biochemical transformations, included hydroxylation, methylation, and isomerization processes involving an array of cytochrome P450s, methyltransferases, and reductases, such as norcoclaurine 6-\u003cem\u003eO\u003c/em\u003e-methyltransferase (6OMT) (Menendez-Perdomo \u0026amp; Facchini, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2020\u003c/span\u003e), coclaurine \u003cem\u003eN\u003c/em\u003e-methyltransferase (CNMT), (\u003cem\u003eS\u003c/em\u003e)-\u003cem\u003eN\u003c/em\u003e-methylcoclaurine 3\u0026rsquo;-hydroxylase (NMCH) (Liu et al., \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Pauli \u0026amp; Kutchan, \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e1998\u003c/span\u003e), and 3'-hydroxy-\u003cem\u003eN\u003c/em\u003e-methylcoclaurine 4'-\u003cem\u003eO\u003c/em\u003e-methyltransferase (4'OMT) (Gurkok et al., \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). These enzymes catalyzed the production of a series of crucial 1-BIA intermediates, including coclaurine (COC), \u003cem\u003eN\u003c/em\u003e-methylcoclaurine (NMC), 3'-hydroxyl-\u003cem\u003eN\u003c/em\u003e-methylcoclaurine (HNMC), and reticuline, which subsequently served as substrates for the biosynthesis of a diverse array of BIAs (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eAporphine alkaloids like crebanine were the major compounds in \u003cem\u003eS. yunnanensis\u003c/em\u003e, which had an unclear biosynthetic pathway. Recent studies have discovered that corytuberine synthase (CTS) like CjCYP80G2 from \u003cem\u003eCoptis japonica\u003c/em\u003e (Ikezawa et al., \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2008\u003c/span\u003e) and StCYP80G6 (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e) from \u003cem\u003eStephania tetrandra\u003c/em\u003e could catalytic C-C coupling of (\u003cem\u003eS\u003c/em\u003e)-reticuline to the type Ⅰ aporphine (\u003cem\u003eS\u003c/em\u003e)-corytuberine. Additionally, NnCYP80Q1 or StCYP80Q5 could catalytic C-C coupling of (\u003cem\u003eR\u003c/em\u003e)-coclaurine and (\u003cem\u003eR\u003c/em\u003e)-\u003cem\u003eN\u003c/em\u003e-methylcoclaurine to the protoaporphines (\u003cem\u003eR\u003c/em\u003e)-crotoflorine and (\u003cem\u003eR\u003c/em\u003e)-glaziovine, two important precursors of type Ⅱ aporphines (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Pyne et al., \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Subsequently, the proaporphines crotoflorine and glaziovine underwent unknown oxidative rearrangements, followed by a series of hydroxylations, methylations, and the formation of methylenedioxy bridges, leading to the production of type II aporphines, including crebanine and nelumboferine (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). These CYP80s were considered to be the starting points of the aporphines biosynthetic pathways, however, no CYP80s with similar catalytic C-C coupling functions had been found in \u003cem\u003eS. yunnanensis\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eThe combination of transcriptomics and metabolomics is currently the most common method for studying the biosynthesis and regulation of plants secondary metabolites (Lai et al., \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Li et al., \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Li et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Tong et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). To date, \u003cem\u003ede novo\u003c/em\u003e transcriptomes of multiple BIA-producing plants have been reported, focusing on the Papaveraceae, Ranunculaceae, and Menispermaceae families within the Ranunculales (Hagel et al., \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Papaveraceae plants have been a primary focus due to their unique accumulation of morphinan alkaloids (Park et al., \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Pei et al., \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Xu et al., \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Zhao et al., \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). As the pharmacological effects of other types BIAs have been discovered, the transcriptomes of other BIA-producing plants have also been reported to enable more in-depth research into BIA biosynthesis. The genome of \u003cem\u003eS. yunnanensis\u003c/em\u003e has been previously reported (Leng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2024\u003c/span\u003e), providing valuable genomic data for investigating its BIA biosynthetic pathways. However, for a deeper understanding of the BIA biosynthesis and its regulation, it is also necessary to have information on gene expression levels and metabolite profiles. Unfortunately, there are currently no reported transcriptome data available for \u003cem\u003eS. yunnanensis\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eIn this study, we determined the relative concentrations of four BIAs in different tissues of \u003cem\u003eS. yunnanensis\u003c/em\u003e to investigate the differential accumulation patterns of these BIAs. We also generated a \u003cem\u003ede novo\u003c/em\u003e transcriptome of \u003cem\u003eS. yunnanensis\u003c/em\u003e and identified hundreds of genes potentially related to BIA biosynthesis. Subsequently, using differential expression analysis and phylogenetic analysis, we identified some candidate genes involved in BIA biosynthesis in \u003cem\u003eS. yunnanensis\u003c/em\u003e and conducted genomic-level analysis on them. We then validated the functions of CYP80s genes of \u003cem\u003eS. yunnanensis in vitro\u003c/em\u003e. Our research will contribute to further in-depth studies on BIA biosynthetic pathways and their evolution in Stephania species, aiming to elucidate the complete biosynthetic pathways of BIAs represented by crebanine.\u003c/p\u003e"},{"header":"2. Materials and methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 Plant materials, chemicals, reagents and strains\u003c/h2\u003e \u003cp\u003eThe \u003cem\u003eS. yunnanensis\u003c/em\u003e plant samples were collected from Yunnan University of Traditional Chinese Medicine in Kunming City, Yunnan Province, China.\u003c/p\u003e \u003cp\u003eAuthentic standards of (\u003cem\u003eS\u003c/em\u003e)-\u003cem\u003eN\u003c/em\u003e-methylcoclaurine and (\u003cem\u003eR\u003c/em\u003e)-\u003cem\u003eN\u003c/em\u003e-methylcoclaurine were obtained through commercial chemical synthesis by WuXi LabNetwork. (\u003cem\u003eS\u003c/em\u003e)-norcoclaurine, (\u003cem\u003eR\u003c/em\u003e)-norcoclaurine, (\u003cem\u003eS\u003c/em\u003e)-coclaurine and (\u003cem\u003eR\u003c/em\u003e)-coclaurine were derived from the chiral separation of commercial racemic standard. The racemic norcoclaurine and coclaurine were purchased from Baoji Herbest Bio-Tech Co., Ltd. Other standards were purchased from Shanghai yuanye Bio-Technology Co., Ltd. (\u003cem\u003eR\u003c/em\u003e)-crotoflorine was prepared in previous work (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eYPD medium was composed of 20 g L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e peptone (OXOID), 10 g L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e yeast extract (OXOID), and 20 g L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e glucose (Solarbio), serving as a standard medium for cultivating yeast and preparing competent cells. YPL medium, used to induce gene expression in yeast, uses galactose (Solarbio) instead of glucose compared to YPD medium. Synthetic dropout minus uracil medium (SD-Ura, FunGenome) with 20 g L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e glucose was employed to select positive colonies transformed with the pESC-Ura vector. For plates preparation, agar (DING GUO) at a concentration of 20 g L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e was added as necessary.\u003c/p\u003e \u003cp\u003eWAT11 yeast strain with endogenous cytochrome P450 reductase replaced by one from Arabidopsis thaliana was maintained by this laboratory (Urban et al., \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e1997\u003c/span\u003e). pESC-Ura was purchased from Agilent Technologies. \u003cem\u003eEscherichia coli\u003c/em\u003e Trans1 T1 strain from TransGen Biotech was used for routine plasmid assembly.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2 Alkaloid extraction and composition analysis\u003c/h2\u003e \u003cp\u003e50 mg root, stem, and leaf freeze-dried powders of \u003cem\u003eS. yunnanensis\u003c/em\u003e were mixed to 2 mL methanol and then sonicated extracts for 0.5 h. After centrifugation, the supernatants were filtered with a nylon syringe filter (0.22 \u0026micro;m). Quantitative analysis was conducted on a UPLC-QTOF-MS system (Waters Technologies). The Acquity UPLC was carried out using a T3 column (Waters Technologies, 2.1 mm \u0026times; 100 mm, 2.7 \u0026micro;m particle size) at 38 ℃. The mobile phases consisted of eluent acetonitrile (A) and 0.1% aqueous formic acid (B) with a flow rate of 0.1 mL min\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e with the linear gradient elution program: 5% \u0026minus;\u0026thinsp;12% A from 0.0 min to 1.0 min, 12% A from 1.0 min to 15.0 min, 12% \u0026minus;\u0026thinsp;17% A from 15.0 min to 16.0 min, 17% A from 16.0 min to 66.0 min, 17% ~ 90% A from 66.0 min to 67.0 min, 90% A from 67.0 min to 72.0 min, 90% \u0026minus;\u0026thinsp;5% A from 72.0 min to 72.5 min, and 5% A from 72.5 min to 76.0 min. The Acquity UPLC system was coupled to a Waters Xevo G2-S QTOF mass spectrometer equipped with electrospray ionization (ESI). The instrument was operated in positive ion mode to perform full scan monitoring in the range of m/z 50\u0026ndash;800. The other operating parameters were set as follows: capillary voltage of 0.5 kV; sample cone voltage of 40 V; extraction cone voltage of 4 V; source temperature of 100 ℃; desolvation temperature 300 ℃; and desolvation gas flow of 800 L h\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. The trap collision energy of low energy function was set at 6 eV, while ramp trap collision energy of high energy function was set at 30\u0026ndash;50 eV (Xiao et al., \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). Data acquisition and processing were performed using MassLynx software.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003e2.3 Transcriptome analysis\u003c/h2\u003e \u003cp\u003e\u003cem\u003eS. yunnanensis\u003c/em\u003e RNA was extracted using the R6827 Plant RNA Kit (Omega Bio-tek) according to the manufacturer's protocol. The RNA was quality-checked using Nanodrop, Qubit, and gel electrophoresis. Once the RNA samples were verified to be of acceptable quality, they were randomly fragmented using a Covaris sonicator. The entire library preparation process then involved end repair, A-tailing, adapter ligation, purification, and PCR amplification. Specifically, fragmented mRNA was used as a template with random oligonucleotides as primers to synthesize the first strand of cDNA in an M-MuLV reverse transcriptase system. The RNA strand was degraded with RNaseH, and the second strand of cDNA was synthesized using DNA polymerase I with dNTPs as substrates. The double-stranded cDNA underwent end repair, A-tailing, and sequencing adapter ligation. AMPure XP beads were used to select cDNA fragments of approximately 250\u0026ndash;300 bp. After PCR amplification, the library was purified again using AMPure XP beads. Quality control of the library involved quantification with a Qubit 2.0 Fluorometer, followed by size verification with an Agilent 2100 bioanalyzer. Once the insert size met expectations, the effective library concentration was accurately quantified using qRT-PCR (with the effective library concentration exceeding 2 nM). The library preparation kit used was the NEBNext\u0026reg; Ultra\u0026trade; RNA Library Prep Kit for Illumina\u0026reg;, and sequencing was performed on an Illumina Novaseq 6000 platform.\u003c/p\u003e \u003cp\u003eTo obtain comprehensive functional information of the transcripts, functional annotation was performed using seven databases, including Nr (Yangyang et al., \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2006\u003c/span\u003e), Pfam (Finn et al., \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), Uniprot (Apweiler et al., \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2004\u003c/span\u003e), KEGG (Kanehisa et al., \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2004\u003c/span\u003e), GO (Ashburner et al., \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2000\u003c/span\u003e), KOG/COG (Tatusov et al., \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2003\u003c/span\u003e), and PATHWAY (Kanehisa et al., \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2004\u003c/span\u003e). Functional annotation primarily employed two methods: sequence similarity search and motif similarity search. For sequence similarity search, the protein sequences encoded by the transcripts were compared with existing protein databases, such as Uniprot, Nr (Yangyang et al., \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2006\u003c/span\u003e), and the metabolic pathway database KEGG (Kanehisa et al., \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2004\u003c/span\u003e), using diamond BLASTp (version: 2.0.6.144; parameters: \u0026ndash;evalue 1e\u003csup\u003e\u0026minus;\u0026thinsp;5\u003c/sup\u003e) (Buchfink et al., \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2015\u003c/span\u003e) to obtain functional information and potential metabolic pathway information. KEGG annotation was performed using KOBAS (version: 3.0) (Xie et al., \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2011\u003c/span\u003e) to associate with KEGG ORTHOLOGY and PATHWAY. The Uniprot (Apweiler et al., \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2004\u003c/span\u003e) database records the correspondence between each protein family and functional nodes in Gene Ontology (Ashburner et al., \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2000\u003c/span\u003e), allowing prediction of the biological functions executed by the protein sequences encoded by the transcripts. Based on the associations between databases, KOG/COG (Tatusov et al., \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2003\u003c/span\u003e) annotation results were obtained, followed by classification statistics and plotting for KOG/COG (Tatusov et al., \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2003\u003c/span\u003e). For motif similarity search, proteins typically consist of one or more functional regions, commonly referred to as domains. The different combinations of these domains result in a variety of proteins; thus, identifying protein domains is crucial for analyzing protein functions. Domain prediction was performed using hmmscan (version: 3.3.2; parameters: e-value 0.01) (Finn et al., \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2011\u003c/span\u003e) to obtain conserved sequences, motifs, and domains of the proteins. The Pfam (Finn et al., \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) database is a large collection of protein families based on multiple sequence alignments and hidden Markov models.\u003c/p\u003e \u003cp\u003eThe expression level of transcripts was assessed by calculating their FPKM value. If the samples all contain biological replicates, used DESeq2 for differential expression analysis; otherwise, used edgeR. Herein, the differentially expressed genes (DEGs) were screened with the thresholds were a significance level of corrected value of padj\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and |log2FoldChange| \u0026gt; 1. Perform GO enrichment and KEGG pathway enrichment analysis on DEGs, with the padj\u0026thinsp;\u0026lt;\u0026thinsp;0.05, which was considered significantly enriched among DEGs, respectively.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003e2.4 Analysis of candidate genes in the BIA biosynthetic pathways\u003c/h2\u003e \u003cp\u003eTo identify the candidate genes associated with BIA biosynthetic pathways of \u003cem\u003eS. yunnanensis\u003c/em\u003e, local BLASTp by BioEdit, version 7.1.3.0 (Hall, \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e1999\u003c/span\u003e) was performed against \u003cem\u003eS. yunnanensis\u003c/em\u003e sequences using query sequences of BIA-producing plants obtained from GenBank databases. The resulting transcripts with high degrees of identity were selected as candidate genes. A heatmap was generated using Tbtools with row scaling. Based on the reported genome of \u003cem\u003eS. yunnanensis\u003c/em\u003e (Leng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2024\u003c/span\u003e), we studied the chromosomal localization of these candidate genes and conducted synteny analysis of three species using TbTools.\u003c/p\u003e \u003cp\u003eBased on the reference genes of BIA-producing plants and the obtained candidate genes, phylogenetic analysis was conducted using MEGA11 (Tamura et al., \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) through the maximum likelihood estimation.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003e2.5 Cloning of candidate genes and eukaryotic expression of recombinant plasmids\u003c/h2\u003e \u003cp\u003eAll candidate genes were amplified from \u003cem\u003eS. yunnanensis\u003c/em\u003e root or leaf cDNA using the primers listed in the supplementary information (Table S4) and PrimeSTAR\u0026reg; Max DNA Polymerase (TaKaRa), according to the manufacturer's instructions. All CYP450 genes were purified and constructed into the pESC-Ura vector (digested with BamH Ⅰ) using the Gibson Assembly Kit (TransGen Biotech). The recombinant vectors were transformed into competent \u003cem\u003eE. coli\u003c/em\u003e Trans1 T1 and identified by colony PCR (TransGen Biotech) and Sanger sequencing (Ruiotechnology). The identified positive recombinant plasmids were extracted from Trans1 T1 using the Plasmid Purification Kit (Magen Biotech).\u003c/p\u003e \u003cp\u003eThe recombinant plasmids pESC-Ura carrying candidate CYP450 genes were each transformed into WAT11 using the Frozen-EZ Yeast Transformation II Kit (Zymo Research). WAT11 transformed with empty pESC-Ura was employed as a control. The positive transformants were screened on SD-Ura medium containing 20 g L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e glucose. The positive clones were incubated with shaking at 30\u0026deg;C for 20\u0026ndash;24 hours until the OD600 reached 2\u0026ndash;3. Cells were centrifuged to remove the SD-Ura medium. The cells were then resuspended in YPL induction medium and grown overnight at 30\u0026deg;C to induce recombinant protein expression.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003e2.6 Microsomes extraction, enzymatic activity assay and LC-MS analysis\u003c/h2\u003e \u003cp\u003eMicrosomes of recombinant yeast were prepared as previously described (Pompon et al., \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e1996\u003c/span\u003e), which has proven to be successful for previous such characterization. \u003cem\u003eIn vitro\u003c/em\u003e activity assays were performed in a 250 \u0026micro;L reaction system that included 100 mM Tris-HCl (pH 7.5, sangong) and 500 \u0026micro;M NADPH (Solarbio), 0.1 g microsomal protein, and 20 \u0026micro;M of the substrate. The reactions were incubated at 30 ℃ for 2 hours with 180 rpm shaking, then 5 \u0026micro;L ammonium hydroxide of was added to the reaction system to make the pH to about 10, and then extracted with 250 \u0026micro;L of ethyl acetate. After the ethyl acetate solution of the reaction product was concentrated to dryness, 150 \u0026micro;L of methanol was added to redissolve and centrifuged at 20,000 g for 15 min before LC-MS analysis. All enzymatic reaction products were detected using UPLC-QTOF-MS system (Waters Technologies). The Acquity UPLC was carried out as previously described. The linear gradient elution program: 5% \u0026minus;\u0026thinsp;30% A from 0.0 min to 6.0 min, 30% \u0026minus;\u0026thinsp;60% A from 6.0 min to 8.0 min, 60% \u0026minus;\u0026thinsp;90% A from 8.0 min to 8.5 min, 90% A from 8.5 min to 9.5 min, 90% \u0026minus;\u0026thinsp;5% A from 9.5 min to 10.0 min, and 5% A from 10 min to 11 min. Data acquisition and processing were performed using MassLynx software.\u003c/p\u003e \u003c/div\u003e"},{"header":"3. Results","content":"\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Determination of BIAs in \u003cem\u003eS. yunnanensis\u003c/em\u003e\u003c/h2\u003e \u003cp\u003e \u003cem\u003eS. yunnanensis\u003c/em\u003e is a medicinal plant with a long history of use, with BIAs as its main active components (Dai et al., \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Hu et al., \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Due to its strict environmental requirements and long growth cycle, it had not yet been cultivated on a large scale. Recent studies have reported the presence of various aporphine alkaloids in \u003cem\u003eS. yunnanensis\u003c/em\u003e (Wang et al., \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) (Fig. S1). We analyzed three tissues of \u003cem\u003eS. yunnanensis\u003c/em\u003e using liquid chromatography-mass spectrometry (LC-MS) and researched the relative contents of salutaridine, roemerine, stephanine, and crebanine to investigate the accumulation patterns of these four representative BIAs in \u003cem\u003eS. yunnanensis\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, Fig. S2). Surprisingly, we detected four representative BIAs in the roots, each corresponding to the purchased reference standards, but in the stems and leaves, there were hardly any or only minimal amounts detected (Fig. S2). Furthermore, crebanine was the most abundant among the four compounds (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e), which is consistent with previous studies (Peng, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Among these compounds, salutaridine was identified as a morphinan, another major compound reported in \u003cem\u003eS. yunnanensis\u003c/em\u003e (Peng, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2014\u003c/span\u003e); Additionally, roemerine and stephanine were considered important intermediates, along with the representative product crebanine, all of which belong to type Ⅱ aporphines (Wang et al., \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). We also attempted to detect additional aporphines that have been reported in other Stephania species (Peng, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Wang et al., \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). However, we did not detect magnoflorine, a representative type Ⅰ aporphine, suggesting that \u003cem\u003eS. yunnanensis\u003c/em\u003e may contain more type Ⅱ aporphines, with type Ⅰ aporphines being less prevalent. This result might explain why \u003cem\u003eS. yunnanensis\u003c/em\u003e was primarily used for its tuberous roots in medicine rather than its aerial parts. However, under the current extraction and detection conditions, we did not find any representative bisbenzylisoquinoline cepharanthine, which contradicts another study that reported the presence of cepharanthine and other bisbenzylisoquinolines in the metabolome of \u003cem\u003eS. yunnanensis\u003c/em\u003e (Leng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). In summary, the analysis of these compounds indicated that numerous structurally diverse BIAs are present in \u003cem\u003eS. yunnanensis\u003c/em\u003e, specifically accumulating in the roots, with aporphines, represented by crebanine, being the most abundant.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003e3.2 Transcriptome sequencing, assembly, and analysis of \u003cem\u003eS. yunnanensis\u003c/em\u003e.\u003c/h2\u003e \u003cp\u003eTo explore the biosynthesis of aporphine alkaloids in \u003cem\u003eS. yunnanensis\u003c/em\u003e, we prepared nine samples from the roots, stems, and leaves for transcriptome sequencing. RNA was extracted and cDNA libraries were constructed using Illumina and PacBio Sequel II for sequencing, resulting in 41,640,198 raw reads. \u003cem\u003eDe novo\u003c/em\u003e assembly yielded 50,119 transcripts, with an N50 length of 2,041 bp and an average length of 1,681 bp (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eAssembly results of \u003cem\u003ede novo\u003c/em\u003e transcriptome of \u003cem\u003eS. yunnanensis\u003c/em\u003e\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"2\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eItem\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNumber\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSeq Num\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e50,119\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSeq Base (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e91,278,165\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eN50 (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2,041\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eMax Length (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e7,046\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eMin Length (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e108\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAverage Length (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1,821.23\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eMean Length (bp)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e1,681.00\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eBenchmarking Universal Single-Copy Orthologs (BUSCO) assessment (Manni et al., \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) indicated that the \u003cem\u003eS. yunnanensis\u003c/em\u003e transcriptome contained 87.3% complete BUSCOs (Fig. S3a, Table S1). Annotation was performed for all 50,119 (100%) transcripts across various databases, with an annotated transcript N50 length of 1,518 bp (Table S2). Gene Ontology (GO) classification (Ashburner et al., \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2000\u003c/span\u003e) was conducted to describe gene functions in three categories: cellular component, molecular function, and biological process (Fig. S3b).\u003c/p\u003e \u003cp\u003eTo investigate the distribution of assembled transcripts in each tissue, their relative expression levels were determined by calculating the fragment per kilobase of transcript per million fragments mapped (FPKM) values of assembled transcripts. To compare the differentially expressed genes (DEGs) between roots and other tissues, we conducted a comparative analysis of the expression levels of assembled transcripts (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). The results indicated that the number of DEGs between roots and leaves was the highest compared to other groups, with 1,780 up-regulated and 932 down-regulated (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea \u0026amp; \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb). To further analyze the DEGs between roots and leaves, GO and KEGG clustering analyses were performed (Figs.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec \u0026amp; \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed, Fig. S4). GO clustering revealed that most DEGs between roots and leaves were classified into \u0026ldquo;monooxygenase activity-related\u0026rdquo; categories, except for photosynthesis-related categories (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec). On the other hand, KEGG clustering showed that a total of 685 DEGs were considered to be involved in secondary metabolism, with only a small number (11 DEGs, with 4 down-regulated and 7 up-regulated) clustered into the \u0026ldquo;isoquinoline alkaloid biosynthesis\u0026rdquo; category (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003e3.3 Analysis of BIA biosynthesis genes\u003c/h2\u003e \u003cp\u003eApart from the two reported SyNCS4, SyNCS5 (Leng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2024\u003c/span\u003e), the biosynthesis of BIAs in \u003cem\u003eS. yunnanensis\u003c/em\u003e has not yet been elucidated. The upstream pathways, from dopamine and 4-HPAA to reticuline, were largely conserved among BIA-producing plants (Lee \u0026amp; Facchini, \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Analysis of the annotation results revealed a total of 141 OMT, 117 NMT, 45 NCSs, and 505 P450s (Tables S5-S8). In recent studies, CYP80s involved in BIA biosynthesis have gained increasing attention for their roles in the synthesis of 1-BIAs, as well as bisbenzylisoquinoline alkaloids and aporphine alkaloids (Hao et al., \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). Notably, among the P450s, 25 were further annotated into the CYP80 family, which are likely involved in BIA biosynthesis. To investigate the BIA biosynthesis in \u003cem\u003eS. yunnanensis\u003c/em\u003e, we identified candidate genes with high similarity to previously reported BIA biosynthesis genes from \u003cem\u003eS. tetrandra\u003c/em\u003e (Li et al., \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Li et al., \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Li et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). We selected sequences with homology greater than 55% and coverage greater than 90% as candidate genes, and a total of 26 candidate transcripts for the BIA biosynthesis were identified from the \u003cem\u003eS. yunnanensis\u003c/em\u003e transcriptome. Phylogenetic analysis and domain analysis indicated that these candidate transcripts have structural similarities to the reference genes (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, Table S3). The sequence alignment revealed that these candidate genes share similar sequences with the reference genes (Fig. S5-S8). To analyze the expression patterns of these candidate transcripts in different tissues of \u003cem\u003eS. yunnanensis\u003c/em\u003e, a heatmap was generated using the FPKM of the candidate genes (Fig. S9). The heatmap showed that, except for NCS, 4'OMT, and CYP80G6, most candidate transcripts did not exhibit distinct differential expression patterns and maintained high expression levels across multiple tissue parts. This may explain why only a few DEGs were clustered into the \"isoquinoline alkaloid biosynthesis\" category. The disparity between gene expression and compound accumulation suggests that BIA compounds transport might occur in \u003cem\u003eS. yunnanensis\u003c/em\u003e. On the other hand, this disparity also suggests that the downstream genes of the BIA biosynthesis might not be highly expressed only in the roots, an important consideration for gene screening.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eFurthermore, based on the reported \u003cem\u003eS. yunnanensis\u003c/em\u003e genome (Leng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2024\u003c/span\u003e), we examined the distribution and location of these candidate transcripts on the chromosomes (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea). Sequence alignment revealed that 26 candidate transcripts corresponded to 18 genes in genome, they were arranged in descending order of homology and named according to their respective functional genes before being mapped onto the chromosomes. Notably, among these candidate genes, one SyNMCH, two SyCNMT and two Sy6OMT were located on pseudochromosome 3 in close proximity, suggesting that this region might form a BIA biosynthetic gene cluster (BGC). These were the genes involved in the three consecutive steps following NCS in the upstream BIA biosynthetic pathways. Similarly, another candidate gene SyCYP80G6 co-localized with two SyCNMT on pseudochromosome 11. Additionally, the candidate genes for CYP80Q5 were found on three different pseudochromosomes. Specifically, four were located on pseudochromosome 1, one was located on pseudochromosome 9, and another was co-localized with Sy4'OMT and SyNMCH2 on pseudochromosome 13. Among them, SyCYP80Q5-1 and SyCYP80Q5-2 shared over 98% identity, indicating they are tandem repeats, whereas they had 90% identity with SyCYP80Q5-3 (Fig. S10). The other three genes, SyCYP80Q5-4, SyCYP80Q5-5, and SyCYP80Q5-6, had low identity with other SyCYP80Q5 genes, approximately 59\u0026ndash;65%. These chromosomal locations preliminarily suggested that the BIA biosynthetic genes in \u003cem\u003eS. yunnanensis\u003c/em\u003e may exhibit a certain degree of clustered distribution. To investigate the distribution and association of these candidate genes in Stephania species, we conducted a synteny analysis using the reported genomes of two other plants, \u003cem\u003eS. japonica\u003c/em\u003e and \u003cem\u003eS. cepharantha\u003c/em\u003e (Leng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2024\u003c/span\u003e) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eb). Interestingly, the two tandemly duplicated genes of SyCYP80Q5 on chromosome 1 of \u003cem\u003eS. yunnanensis\u003c/em\u003e, SyCYP80Q5-1 and SyCYP80Q5-2 had only one ortholog in \u003cem\u003eS. cepharantha\u003c/em\u003e and none in \u003cem\u003eS. japonica\u003c/em\u003e, but SyCYP80Q5-3 was conserved across the three species. Notably, the clusters observed in \u003cem\u003eS. yunnanensis\u003c/em\u003e, were present on different chromosomes in the three Stephania species. This indicated that the clustering patterns of BIA biosynthetic genes, were conserved across the three species, suggesting that such clustering is likely conserved within Stephania as well. In summary, chromosome localization and synteny analysis have provided insights into the genomic-level associations of BIA biosynthetic genes in \u003cem\u003eS. yunnanensis\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003e3.4 Functional verification and phylogenetic analysis of CYP80s\u003c/h2\u003e \u003cp\u003eTo elucidate the biosynthesis of the BIAs of \u003cem\u003eS. yunnanensis\u003c/em\u003e, we conducted \u003cem\u003ein vitro\u003c/em\u003e functional validation of the two key CYP80s involved in aporphine skeleton formation: CYP80G6 (CTS) and CYP80Q5. These enzymes specifically catalyzed the (\u003cem\u003eS\u003c/em\u003e) or (\u003cem\u003eR\u003c/em\u003e)- types of 1-BIA substrates to produce the corresponding type Ⅰ aporphines or type Ⅱ aporphines (protoaporphines). The further rearrangement and modification of these intermediates ultimately led to the formation of aporphines such as crebanine. Since the high identity between SyCYP80Q5-1 and SyCYP80Q5-2, we cloned only SyCYP80Q5-1 and another four, SyCYP80Q5-3, SyCYP80Q5-4, SyCYP80Q5-5, SyCYP80Q5-6, as well as SyCYP80G6, and validated their functions of them \u003cem\u003ein vitro\u003c/em\u003e. Recombinant plasmids containing these candidate genes were constructed and expressed in \u003cem\u003eSaccharomyces cerevisiae\u003c/em\u003e (WAT11), in which the endogenous cytochrome P450 reductase was replaced by one from \u003cem\u003eArabidopsis thaliana\u003c/em\u003e. Microsomes extracted from the recombinant yeast were used for \u003cem\u003ein vitro\u003c/em\u003e enzyme assays and were detected by LC-MS. We found that only SyCYP80Q5-1, SyCYP80Q5-3 and SyCYP80G6 exhibited catalytic functions (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e, Figs. S11 \u0026amp; S12). Specifically, SyCYP80G6 catalyzes the conversion of (\u003cem\u003eS\u003c/em\u003e)-reticuline and (\u003cem\u003eS\u003c/em\u003e)-NMC into (\u003cem\u003eS\u003c/em\u003e)-corytuberine and (\u003cem\u003eS\u003c/em\u003e)-glaziovine (Figs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ea-\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ed, Figs. S11a \u0026amp; S11c). Both SyCYP80Q5-1 and SyCYP80Q5-3 were capable of catalyzing the conversion of (\u003cem\u003eR\u003c/em\u003e)-COC and (\u003cem\u003eR\u003c/em\u003e)-NMC into the aporphine alkaloids (\u003cem\u003eR\u003c/em\u003e)-crotoflorine and (\u003cem\u003eR\u003c/em\u003e)-glaziovine (Figs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ee-\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eh, Figs. S11b \u0026amp; S12). Furthermore, the catalytic efficiency of SyCYP80Q5-1 was slightly higher than that of SyCYP80Q5-3 when using either (\u003cem\u003eR\u003c/em\u003e)-COC or (\u003cem\u003eR\u003c/em\u003e)-NMC as substrates (Figs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ef \u0026amp; \u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eh).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eAdditionally, we discussed the activity of SyCYP80Q5-1 and SyCYP80G6 on a broader range of BIA substrates and other configurations (Fig. S13 \u0026amp; S14). As expected, SyCYP80Q5 and SyCYP80G6 exhibited strict substrate specificity and stereoselectivity, which was consistent with the results of other CYP80s (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). We performed a phylogenetic analysis of the CYP80 family, and the results showed that the three SyCYP80s clustered within the CYP80 family (Fig.\u0026nbsp;7), and they were distributed across various clades of the CYP80 family. The functions of these enzymes matched the previously reported enzymes (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e), indicating that these three enzymes play similar roles in the BIA biosynthesis in \u003cem\u003eS. yunnanensis\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eThe biosynthetic pathways of BIAs have not been fully elucidated, mainly due to the complex transformations from 1-BIAs to various types of BIAs scaffolds and the subsequent downstream modifications, including C-C coupling, C-O coupling, hydroxylation, and methylation, etc. (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Meng et al., \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). This study reported the \u003cem\u003ede novo\u003c/em\u003e transcriptome of \u003cem\u003eS. yunnanensis\u003c/em\u003e, and through the combined analysis of the transcriptome, metabolome, and genome, identified 18 candidate genes involved in BIA biosynthesis in \u003cem\u003eS. yunnanensis\u003c/em\u003e. The functions of three CYP80s were validated \u003cem\u003ein vitro\u003c/em\u003e, they specifically catalyzed the C-C coupling of (\u003cem\u003eS\u003c/em\u003e)- or (\u003cem\u003eR\u003c/em\u003e)-configured 1-BIA substrates, resulting in the formation of aporphines or protoaporphines. This is an important step in the biosynthetic pathways of aporphines (Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn medicinal plants where roots or rhizomes are used, the active components are often highly concentrated in the underground parts, with corresponding functional genes and regulatory factors showing significantly higher expression levels compared to the aboveground parts (He et al., \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Liu et al., \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Liu et al., \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Tong et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Wang et al., \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Zhan et al., \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). In this study, we analyzed the representative BIAs of \u003cem\u003eS. yunnanensis\u003c/em\u003e and found that the major active component, crebanine and three other BIAs were highly accumulated in the roots, similar to other root-based medicinal plants. Notably, no upstream 1-BIAs were detected in the metabolomes of any tissues, possibly due to their complete consumption or concentrations below the detection limit. Furthermore, under the extraction and detection conditions of this study, the representative bisbenzylisoquinoline alkaloid, cepharanthine, was not detected.\u003c/p\u003e \u003cp\u003eThen we sequenced the transcriptomes of \u003cem\u003eS. yunnanensis\u003c/em\u003e using next-generation sequencing, followed by \u003cem\u003ede novo\u003c/em\u003e assembly and functional annotation. Previous studies have shown that P450s are involved in almost the entire biosynthetic pathways of BIAs, specifically CYP80, which participates in hydroxylation and downstream C-O and C-C coupling reactions (An et al., \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Hao et al., \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Meng et al., \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). CYP719 catalyzed the formation of the methylenedioxy bridge (Hori et al., \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Ikezawa et al., \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2003\u003c/span\u003e; Li et al., \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; Menendez-Perdomo \u0026amp; Facchini, \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), a characteristic group of many active BIA components, including crebanine. We identified 25 transcripts annotated as \u0026ldquo;CYP80\u0026rdquo; and 11 were annotated as \u0026ldquo;CYP719\u0026rdquo; in the \u003cem\u003ede novo\u003c/em\u003e transcriptome of \u003cem\u003eS. yunnanensis\u003c/em\u003e. Through annotation and sequence alignment, we identified 26 candidate transcripts involved in 10 steps of the BIA biosynthetic pathways. The heatmap showed that, except for a few genes, the majority were not specifically expressed in the roots but were highly expressed across multiple tissues. This pattern of metabolite accumulation and differential gene expression levels, while uncommon, is not unprecedented in BIA-producing plants. Similar mechanisms involving transporter proteins facilitating the movement of compounds within the plant have been reported previously in \u003cem\u003eC. japonica\u003c/em\u003e (Sakai et al., \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2002\u003c/span\u003e; Shitan et al., \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2003\u003c/span\u003e)d \u003cem\u003esomniferum\u003c/em\u003e (Dastmalchi et al., \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). We speculate that similar mechanisms may also exist in \u003cem\u003eS. yunnanensis\u003c/em\u003e, facilitating the internal transfer of compounds through processes such as endocytosis and exocytosis. In summary, the transcriptome data obtained in this study are valuable for elucidating the biosynthetic pathways of BIAs, including crebanine.\u003c/p\u003e \u003cp\u003eBGCs, composed of tightly arranged genes that collectively participate in the biosynthesis of specific metabolites, allow organisms to coordinate and efficiently regulate gene expression and synergistic metabolic reactions (Li et al., \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). With the reported high-quality genome of \u003cem\u003eS. yunnanensis\u003c/em\u003e (Leng et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2024\u003c/span\u003e), this study also examined the chromosomal localization of the 18 candidate genes and found a tendency for several candidate genes to cluster. Further collinearity analysis showed that this clustering phenomenon is conserved among different Stephania species, indicating important functions in their evolutionary processes. Notably, the tandem duplication of SyCYP80Q5-1 and SyCYP80Q5-2 on chromosome 1 of \u003cem\u003eS. yunnanensis\u003c/em\u003e was not conserved among different species, whereas another SyCYP80Q5-3 gene on chromosome 13 was widely conserved. This might explain why all Stephania species can produce aporphines and protoaporphines, but differ in the types and amounts they contain. Furthermore, the number of CYP80Q5 copies in \u003cem\u003eS. yunnanensis\u003c/em\u003e was significantly higher than CYP80G6, which was consistent with the finding that there are more type II aporphines present in \u003cem\u003eS. yunnanensis\u003c/em\u003e. In conclusion, the differences and conservation of the BIA biosynthetic pathways in Stephania suggest that there are still variations within the genus. This may explain why there are significant differences in the types and contents of BIAs among different Stephania species.\u003c/p\u003e \u003cp\u003eAs an important gene family in the BIA biosynthetic pathways, CYP80, this study further validated the functions of three CYP80s from \u003cem\u003eS. yunnanensis in vitro\u003c/em\u003e, showing that they all had the ability to catalyze substrate C-C coupling with configurational selectivity. Specifically, SyCYP80G6 specifically catalyzed the production of corresponding aporphines or protoaporphines from (\u003cem\u003eS\u003c/em\u003e)-type substrates; SyCYP80Q5-1 and SyCYP80Q5-3 specifically catalyzed the production of corresponding protoaporphines, the precursor of type Ⅱ aporphines, from (\u003cem\u003eR\u003c/em\u003e)-type substrates. The functions of CYP80 are complex and diverse, and there might be some evolutionary correlation. It is noteworthy that, although SyCYP80G6 can also catalyze (\u003cem\u003eS\u003c/em\u003e)-type substrates to produce aporphines and protoaporphines, similar to other CYP80Gs, we did not detect any representative reported (\u003cem\u003eS\u003c/em\u003e)-type aporphines, such as magnoflorine or mecambroline, in \u003cem\u003eS. yunnanensis\u003c/em\u003e. However, the expression level of SyCYP80G6 is not low. This discrepancy between metabolites and gene expression suggests that there may be mechanisms in \u003cem\u003eS. yunnanensis\u003c/em\u003e's aporphine biosynthesis that we have yet to understand. It is possible that the metabolic flux is diverted towards (R)-reticuline, leading to the biosynthesis of (\u003cem\u003eR\u003c/em\u003e)-type aporphines, morphinans, and protoberberines.\u003c/p\u003e \u003cp\u003eIn conclusion, the combined approach of metabolite analysis and transcriptome sequencing identified 26 candidate transcripts responsible for the biosynthesis of BIAs, including crebanine, in \u003cem\u003eS. yunnanensis\u003c/em\u003e, these genes might be responsible for BIAs biosynthesis in \u003cem\u003eS. yunnanensis\u003c/em\u003e. Furthermore, genome analysis of the 18 candidate genes showed clustering on chromosomes, suggesting the presence of related BGCs, and collinearity analysis indicated these BGCs are conserved among different Stephania species. Finally, we identified three CYP80s from the \u003cem\u003eS. yunnanensis\u003c/em\u003e transcriptome data as involved in BIAs biosynthesis using \u003cem\u003ein vitro\u003c/em\u003e enzymatic reaction. Overall, our work provides valuable genetic information on \u003cem\u003eS. yunnanensis\u003c/em\u003e and reveals the biosynthesis of BIAs in this medicinal plant.\u003c/p\u003e"},{"header":"Declarations","content":"\u003ch2\u003eCompeting Interests\u003c/h2\u003e \u003cp\u003eThe authors have no relevant financial or non-financial interests to disclose.\u003c/p\u003e\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eThis work was supported by the National Key R\u0026amp;D Program of China (2020YFA0908000), the National Natural Science Foundation of China (82011530137, 31961133007), Scientific and technological innovation project of CACMS (CI2023D002, CI2023E002) and Key project at central government level: The ability to establish sustainable use of valuable Chinese medicine resources (2060302).\u003c/p\u003e\u003ch2\u003eAuthor contributions\u003c/h2\u003e \u003cp\u003eLH and JGuo conceived and designed the entire research plans; WS and QL performed most of the experiments; LL, XL, JGan, YM, TC, and PS participated in some of the experiments; WS analyzed results of all the experiments; WS, JGuo, JW, YZ, XM, and LH wrote the manuscript; QL and YM revised the manuscript. All authors read and approved the manuscript.\u003c/p\u003e\u003ch2\u003eData availability\u003c/h2\u003e \u003cp\u003eThe raw data from the transcriptome sequencing have been uploaded to public database under the accession number XXXXXXXX. All other subsequent data are provided in the supplementary materials.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eAn Z, Gao R, Chen S, Tian Y, Li Q, Tian L, Zhang W, Kong L, Zheng B, Hao L, Xin T, Yao H, Wang Y, Song W, Hua X, Liu C, Song J, Fan H, Sun W, Chen S, Xu Z (2024) Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots. Adv Sci (Weinh) 11(19):e2309990. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/advs.202309990\u003c/span\u003e\u003cspan address=\"10.1002/advs.202309990\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eApweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LS (2004) UniProt: the Universal Protein knowledgebase. Nucleic Acids Res 32(Database issue) D115-119. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gkh131\u003c/span\u003e\u003cspan address=\"10.1093/nar/gkh131\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAshburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25(1):25\u0026ndash;29. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/75556\u003c/span\u003e\u003cspan address=\"10.1038/75556\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBuchfink B, Xie C, Huson DH (2015) Fast and sensitive protein alignment using DIAMOND. Nat Methods 12(1):59\u0026ndash;60. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nmeth.3176\u003c/span\u003e\u003cspan address=\"10.1038/nmeth.3176\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCui L, Peng C, Li J, Cheng X, Fan X, Li J, Yang Z, Zhao Y, Ma Y (2022) The anti-inflammatory and analgesic activities of 2Br-Crebanine and Stephanine from Stephania yunnanenses H. S.Lo. Front Pharmacol 13:1092583. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fphar.2022.1092583\u003c/span\u003e\u003cspan address=\"10.3389/fphar.2022.1092583\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDai X, Hu R, Sun C, Pan Y (2012) Comprehensive separation and analysis of alkaloids from Stephania yunnanensis by counter-current chromatography coupled with liquid chromatography tandem mass spectrometry analysis. J Chromatogr A 1226:18\u0026ndash;23. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.chroma.2011.10.022\u003c/span\u003e\u003cspan address=\"10.1016/j.chroma.2011.10.022\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDastmalchi M, Chang L, Chen R, Yu L, Chen X, Hagel JM, Facchini PJ (2019) Purine Permease-Type Benzylisoquinoline Alkaloid Transporters in Opium Poppy. Plant Physiol 181(3):916\u0026ndash;933. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1104/pp.19.00565\u003c/span\u003e\u003cspan address=\"10.1104/pp.19.00565\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFinn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer EL, Tate J, Punta M (2014) Pfam: the protein families database. Nucleic Acids Res 42(Database issue):D222\u0026ndash;230. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gkt1223\u003c/span\u003e\u003cspan address=\"10.1093/nar/gkt1223\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFinn RD, Clements J, Eddy SR (2011) HMMER web server: interactive sequence similarity searching. Nucleic Acids Res W29\u0026ndash;37, 39(Web Server issue). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gkr367\u003c/span\u003e\u003cspan address=\"10.1093/nar/gkr367\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGhirga F, Bonamore A, Calisti L, D'Acquarica I, Mori M, Botta B, Boffi A, Macone A (2017) Green Routes for the Production of Enantiopure Benzylisoquinoline Alkaloids. Int J Mol Sci 18(11). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/ijms18112464\u003c/span\u003e\u003cspan address=\"10.3390/ijms18112464\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGurkok T, Ozhuner E, Parmaksiz I, Ozcan S, Turktas M, Ipek A, Demirtas I, Okay S, Unver T (2016) Functional Characterization of 4'OMT and 7OMT Genes in BIA Biosynthesis. Front Plant Sci 7:98. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fpls.2016.00098\u003c/span\u003e\u003cspan address=\"10.3389/fpls.2016.00098\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHagel JM, Morris JS, Lee EJ, Desgagne-Penix I, Bross CD, Chang L, Chen X, Farrow SC, Zhang Y, Soh J, Sensen CW, Facchini PJ (2015) Transcriptome analysis of 20 taxonomically related benzylisoquinoline alkaloid-producing plants. BMC Plant Biol 15:227. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s12870-015-0596-0\u003c/span\u003e\u003cspan address=\"10.1186/s12870-015-0596-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHall TA (1999) BIOEDIT: A USER-FRIENDLY BIOLOGICAL SEQUENCE ALIGNMENT EDITOR AND ANALYSIS PROGRAM FOR WINDOWS 95/98/ NT\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHao C, Yu Y, Liu Y, Liu A, Chen S (2024) The CYP80A and CYP80G Are Involved in the Biosynthesis of Benzylisoquinoline Alkaloids in the Sacred Lotus (Nelumbo nucifera). Int J Mol Sci 25(2). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/ijms25020702\u003c/span\u003e\u003cspan address=\"10.3390/ijms25020702\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHe SM, Liang YL, Cong K, Chen G, Zhao X, Zhao QM, Zhang JJ, Wang X, Dong Y, Yang JL, Zhang GH, Qian ZL, Fan W, Yang SC (2018) Identification and Characterization of Genes Involved in Benzylisoquinoline Alkaloid Biosynthesis in Coptis Species. Front Plant Sci 9:731. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fpls.2018.00731\u003c/span\u003e\u003cspan address=\"10.3389/fpls.2018.00731\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHeming Y, Xianrui L (1980) BOTANICAL AND PHARMACOGNOSTICAL STUDIES OF THE CHINESE DRUG SHAN-WU-GUI. Acta Pharm Sinica (11), 674\u0026ndash;683\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHori K, Yamada Y, Purwanto R, Minakuchi Y, Toyoda A, Hirakawa H, Sato F (2018) Mining of the Uncharacterized Cytochrome P450 Genes Involved in Alkaloid Biosynthesis in California Poppy Using a Draft Genome Sequence. Plant Cell Physiol 59(2):222\u0026ndash;233. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/pcp/pcx210\u003c/span\u003e\u003cspan address=\"10.1093/pcp/pcx210\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHu R, Dai X, Lu Y, Pan Y (2010) Preparative separation of isoquinoline alkaloids from Stephania yunnanensis by pH-zone-refining counter-current chromatography. J Chromatogr B Analyt Technol Biomed Life Sci 878(21):1881\u0026ndash;1884. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jchromb.2010.05.005\u003c/span\u003e\u003cspan address=\"10.1016/j.jchromb.2010.05.005\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIkezawa N, Iwasa K, Sato F (2008) Molecular cloning and characterization of CYP80G2, a cytochrome P450 that catalyzes an intramolecular C-C phenol coupling of (S)-reticuline in magnoflorine biosynthesis, from cultured Coptis japonica cells. J Biol Chem 283(14):8810\u0026ndash;8821. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M705082200\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M705082200\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIkezawa N, Tanaka M, Nagayoshi M, Shinkyo R, Sakaki T, Inouye K, Sato F (2003) Molecular cloning and characterization of CYP719, a methylenedioxy bridge-forming enzyme that belongs to a novel P450 family, from cultured Coptis japonica cells. J Biol Chem 278(40):38557\u0026ndash;38565. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M302470200\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M302470200\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M (2004) The KEGG resource for deciphering the genome. Nucleic Acids Res 32(Database issue). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gkh063\u003c/span\u003e\u003cspan address=\"10.1093/nar/gkh063\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e. D277-280\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLai Y, Ma J, Zhang X, Xuan X, Zhu F, Ding S, Shang F, Chen Y, Zhao B, Lan C, Unver T, Huo G, Li X, Wang Y, Liu Y, Lu M, Pan X, Yang D, Li M, Zhang B, Zhang D (2024) High-quality chromosome-level genome assembly and multi-omics analysis of rosemary (Salvia rosmarinus) reveals new insights into the environmental and genome adaptation. Plant Biotechnol J 22(7):1833\u0026ndash;1847. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/pbi.14305\u003c/span\u003e\u003cspan address=\"10.1111/pbi.14305\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLee EJ, Facchini P (2010) Norcoclaurine synthase is a member of the pathogenesis-related 10/Bet v1 protein family. Plant Cell 22(10):3489\u0026ndash;3503. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1105/tpc.110.077958\u003c/span\u003e\u003cspan address=\"10.1105/tpc.110.077958\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLeng L, Xu Z, Hong B, Zhao B, Tian Y, Wang C, Yang L, Zou Z, Li L, Liu K, Peng W, Liu J, An Z, Wang Y, Duan B, Hu Z, Zheng C, Zhang S, Li X, Li M, Liu Z, Bi Z, He T, Liu B, Fan H, Song C, Tong Y, Chen S (2024) Cepharanthine analogs mining and genomes of Stephania accelerate anti-coronavirus drug discovery. Nat Commun 15(1):1537. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-024-45690-5\u003c/span\u003e\u003cspan address=\"10.1038/s41467-024-45690-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi K, Chen X, Zhang J, Wang C, Xu Q, Hu J, Kai G, Feng Y (2022) Transcriptome Analysis of Stephania tetrandra and Characterization of Norcoclaurine-6-O-Methyltransferase Involved in Benzylisoquinoline Alkaloid Biosynthesis. Front Plant Sci 13:874583. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fpls.2022.874583\u003c/span\u003e\u003cspan address=\"10.3389/fpls.2022.874583\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi L, Jiang W, Lu Y (2017) New strategies and approaches for engineering biosynthetic gene clusters of microbial natural products. Biotechnol Adv 35(8):936\u0026ndash;949. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.biotechadv.2017.03.007\u003c/span\u003e\u003cspan address=\"10.1016/j.biotechadv.2017.03.007\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi Q, Bu J, Ma Y, Yang J, Hu Z, Lai C, Xu Y, Tang J, Cui G, Wang Y, Zhao Y, Jin B, Shen Y, Guo J, Huang L (2020) Characterization of O-methyltransferases involved in the biosynthesis of tetrandrine in Stephania tetrandra. J Plant Physiol 250:153181. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jplph.2020.153181\u003c/span\u003e\u003cspan address=\"10.1016/j.jplph.2020.153181\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi Q, Jiao X, Li X, Shi W, Ma Y, Tan X, Gan J, Liu J, Yang J, Wang J, Jin B, Chen T, Su P, Zhao Y, Zhang Y, Tang J, Cui G, Chen Y, Guo J, Huang L (2024) Identification of the cytochrome P450s responsible for the biosynthesis of two types of aporphine alkaloids and their de novo biosynthesis in yeast. J Integr Plant Biol. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/jipb.13724\u003c/span\u003e\u003cspan address=\"10.1111/jipb.13724\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi X, Li Q, Jiao X, Tang H, Cheng Y, Ma Y, Cui G, Tang J, Chen Y, Guo J, Huang L (2023) Phylogenetic analysis and functional characterization of norcoclaurine synthase involved in benzylisoquinoline alkaloids biosynthesis in Stephania tetrandra. J Cell Physiol. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/jcp.31065\u003c/span\u003e\u003cspan address=\"10.1002/jcp.31065\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu S, Liu M, Wang S, Lin Y, Zhang H, Wang Q, Zhao Y (2017) Analysis of the Panax ginseng stem/leaf transcriptome and gene expression during the leaf expansion period. Mol Med Rep 16(5):6396\u0026ndash;6404. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3892/mmr.2017.7377\u003c/span\u003e\u003cspan address=\"10.3892/mmr.2017.7377\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu X, Bu J, Ma Y, Chen Y, Li Q, Jiao X, Hu Z, Cui G, Tang J, Guo J, Huang L (2021) Functional characterization of (S)-N-methylcoclaurine 3'-hydroxylase (NMCH) involved in the biosynthesis of benzylisoquinoline alkaloids in Corydalis yanhusuo. Plant Physiol Biochem 168:507\u0026ndash;515. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.plaphy.2021.09.042\u003c/span\u003e\u003cspan address=\"10.1016/j.plaphy.2021.09.042\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu X, Shao P, Wang Y, Chen Y, Cui S (2023) Anti-inflammatory mechanism of the optimized active ingredients of Sargentodoxa cuneata and Patrinia villosa. Int Immunopharmacol 120:110337. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.intimp.2023.110337\u003c/span\u003e\u003cspan address=\"10.1016/j.intimp.2023.110337\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eManni M, Berkeley MR, Seppey M, Simao FA, Zdobnov EM (2021) BUSCO Update: Novel and Streamlined Workflows along with Broader and Deeper Phylogenetic Coverage for Scoring of Eukaryotic, Prokaryotic, and Viral Genomes. Mol Biol Evol 38(10):4647\u0026ndash;4654. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/molbev/msab199\u003c/span\u003e\u003cspan address=\"10.1093/molbev/msab199\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMenendez-Perdomo IM, Facchini PJ (2020) Isolation and characterization of two O-methyltransferases involved in benzylisoquinoline alkaloid biosynthesis in sacred lotus (Nelumbo nucifera). J Biol Chem 295(6):1598\u0026ndash;1612. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.RA119.011547\u003c/span\u003e\u003cspan address=\"10.1074/jbc.RA119.011547\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMenendez-Perdomo IM, Facchini PJ (2023) Elucidation of the (R)-enantiospecific benzylisoquinoline alkaloid biosynthetic pathways in sacred lotus (Nelumbo nucifera). Sci Rep 13(1):2955. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41598-023-29415-0\u003c/span\u003e\u003cspan address=\"10.1038/s41598-023-29415-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMeng F, Zhang S, Su J, Zhu B, Pan X, Qiu X, Cui X, Wang C, Niu L, Li C, Lu S (2024) Characterization of two CYP80 enzymes provides insights into aporphine alkaloid skeleton formation in Aristolochia contorta. Plant J 118(5):1439\u0026ndash;1454. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/tpj.16686\u003c/span\u003e\u003cspan address=\"10.1111/tpj.16686\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePark S, An B, Park S (2024) Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae). BMC Plant Biol 24(1):303. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s12870-024-05025-4\u003c/span\u003e\u003cspan address=\"10.1186/s12870-024-05025-4\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePauli HH, Kutchan TM (1998) Molecular cloning and functional heterologous expression of two alleles encoding (S)-N-methylcoclaurine 3'-hydroxylase (CYP80B1), a new methyl jasmonate-inducible cytochrome P-450-dependent mono-oxygenase of benzylisoquinoline alkaloid biosynthesis. Plant J 13(6):793\u0026ndash;801. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1046/j.1365-313x.1998.00085.x\u003c/span\u003e\u003cspan address=\"10.1046/j.1365-313x.1998.00085.x\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePei L, Wang B, Ye J, Hu X, Fu L, Li K, Ni Z, Wang Z, Wei Y, Shi L, Zhang Y, Bai X, Jiang M, Wang S, Ma C, Li S, Liu K, Li W, Cong B (2021) Genome and transcriptome of Papaver somniferum Chinese landrace CHM indicates that massive genome expansion contributes to high benzylisoquinoline alkaloid biosynthesis. Hortic Res 8(1):5. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41438-020-00435-5\u003c/span\u003e\u003cspan address=\"10.1038/s41438-020-00435-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePeng F (2014) \u003cem\u003eChemical constituents from Stephania yunnannensis\u003c/em\u003e [master, Yunnan University of Chinese Medicine]\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePompon D, Louerat B, Bronine A, Urban P (1996) Yeast expression of animal and plant P450s in optimized redox environments. Methods Enzymol 272:51\u0026ndash;64. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/s0076-6879(96)72008-6\u003c/span\u003e\u003cspan address=\"10.1016/s0076-6879(96)72008-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePyne ME, Gold ND, Martin VJJ (2023) Pathway elucidation and microbial synthesis of proaporphine and bis-benzylisoquinoline alkaloids from sacred lotus (Nelumbo nucifera). Metab Eng 77:162\u0026ndash;173. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ymben.2023.03.010\u003c/span\u003e\u003cspan address=\"10.1016/j.ymben.2023.03.010\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSakai K, Shitan N, Sato F, Ueda K, Yazaki K (2002) Characterization of berberine transport into Coptis japonica cells and the involvement of ABC protein. J Exp Bot 53(376):1879\u0026ndash;1886. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/jxb/erf052\u003c/span\u003e\u003cspan address=\"10.1093/jxb/erf052\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShitan N, Bazin I, Dan K, Obata K, Kigawa K, Ueda K, Sato F, Forestier C, Yazaki K (2003) Involvement of CjMDR1, a plant multidrug-resistance-type ATP-binding cassette protein, in alkaloid transport in Coptis japonica. Proc Natl Acad Sci U S A 100(2):751\u0026ndash;756. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1073/pnas.0134257100\u003c/span\u003e\u003cspan address=\"10.1073/pnas.0134257100\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTamura K, Stecher G, Kumar S (2021) MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol Biol Evol 38(7):3022\u0026ndash;3027. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/molbev/msab120\u003c/span\u003e\u003cspan address=\"10.1093/molbev/msab120\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTan J, Xiang Y, Xiong Y, Zhang Y, Qiao B, Zhang H (2023) Crebanine induces ROS-dependent apoptosis in human hepatocellular carcinoma cells via the AKT/FoxO3a signaling pathway. Front Pharmacol 14:1069093. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fphar.2023.1069093\u003c/span\u003e\u003cspan address=\"10.3389/fphar.2023.1069093\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA (2003) The COG database: an updated version includes eukaryotes. BMC Bioinformatics 4:41. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/1471-2105-4-41\u003c/span\u003e\u003cspan address=\"10.1186/1471-2105-4-41\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTong Q, Zhang C, Tu Y, Chen J, Li Q, Zeng Z, Wang F, Sun L, Huang D, Li M, Qiu S, Chen W (2022) Biosynthesis-based spatial metabolome of Salvia miltiorrhiza Bunge by combining metabolomics approaches with mass spectrometry-imaging. Talanta 238(Pt 2):123045. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.talanta.2021.123045\u003c/span\u003e\u003cspan address=\"10.1016/j.talanta.2021.123045\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eUrban P, Mignotte C, Kazmaier M, Delorme F, Pompon D (1997) Cloning, yeast expression, and characterization of the coupling of two distantly related Arabidopsis thaliana NADPH-cytochrome P450 reductases with P450 CYP73A5. J Biol Chem 272(31):19176\u0026ndash;19186. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.272.31.19176\u003c/span\u003e\u003cspan address=\"10.1074/jbc.272.31.19176\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang H, Cheng X, Kong S, Yang Z, Wang H, Huang Q, Li J, Chen C, Ma Y (2016) Synthesis and Structure-Activity Relationships of a Series of Aporphine Derivatives with Antiarrhythmic Activities and Acute Toxicity. Molecules 21(12). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/molecules21121555\u003c/span\u003e\u003cspan address=\"10.3390/molecules21121555\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang K, Jiang S, Sun C, Lin Y, Yin R, Wang Y, Zhang M (2015) The Spatial and Temporal Transcriptomic Landscapes of Ginseng, Panax ginseng C. A. Meyer. Sci Rep 5:18283. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/srep18283\u003c/span\u003e\u003cspan address=\"10.1038/srep18283\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang M, Zhang XM, Fu X, Zhang P, Hu WJ, Yang BY, Kuang HX (2022) Alkaloids in genus stephania (Menispermaceae): A comprehensive review of its ethnopharmacology, phytochemistry, pharmacology and toxicology. J Ethnopharmacol 293:115248. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jep.2022.115248\u003c/span\u003e\u003cspan address=\"10.1016/j.jep.2022.115248\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXiao J, Song N, Lu T, Pan Y, Song J, Chen G, Sun L, Li N (2018) Rapid characterization of TCM Qianjinteng by UPLC-QTOF-MS and its application in the evaluation of three species of Stephania. J Pharm Biomed Anal 156:284\u0026ndash;296. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jpba.2018.04.044\u003c/span\u003e\u003cspan address=\"10.1016/j.jpba.2018.04.044\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L (2011) KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res W316\u0026ndash;322, 39(Web Server issue). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gkr483\u003c/span\u003e\u003cspan address=\"10.1093/nar/gkr483\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXu T, Yang X, Jia Y, Li Z, Tang G, Li X, Wang B, Wang T, Lin J, Guo L, Ye K (2022) A global survey of the transcriptome of the opium poppy (Papaver somniferum) based on single-molecule long-read isoform sequencing. Plant J 110(2):607\u0026ndash;620. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/tpj.15689\u003c/span\u003e\u003cspan address=\"10.1111/tpj.15689\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang Y, Hao T, Yao X, Che Y, Liu Y, Fang M, Wang Y, Zhou D, Chai H, Li N, Hou Y (2023) Crebanine ameliorates ischemia-reperfusion brain damage by inhibiting oxidative stress and neuroinflammation mediated by NADPH oxidase 2 in microglia. Phytomedicine 120:155044. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.phymed.2023.155044\u003c/span\u003e\u003cspan address=\"10.1016/j.phymed.2023.155044\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYangyang D, Jianqi LI, Songfeng WU, Yunping Z, Yaowen C, Fuchu HE (2006) Integrated nr Database in Protein Annotation System and Its Localization. Comput Eng 32(5):71\u0026ndash;72\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYeh PS, Liu CT, Yu CY, Chang YC, Lin SY, Li YC, Luan YZ, Sung WW (2024) Crebanine, an aporphine alkaloid, induces cancer cell apoptosis through PI3K-Akt pathway in glioblastoma multiforme. Front Pharmacol 15:1419044. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fphar.2024.1419044\u003c/span\u003e\u003cspan address=\"10.3389/fphar.2024.1419044\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhan Z, Fang W, Ma X, Chen T, Cui G, Ma Y, Kang L, Nan T, Lin H, Tang J, Zhang Y, Lai C, Ren Z, Wang Y, Zhao Y, Shen Y, Wang L, Zeng W, Guo J, Huang L (2019) Metabolome and transcriptome analyses reveal quality change in the orange-rooted Salvia miltiorrhiza (Danshen) from cultivated field. Chin Med 14:42. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s13020-019-0265-6\u003c/span\u003e\u003cspan address=\"10.1186/s13020-019-0265-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao Y, Cui L, Yang XX, Sun X, Liu Y, Yang Z, Zhu L, Peng C, Li D, Cai J, Ma Y (2021) Sinoacutine inhibits inflammatory responses to attenuates acute lung injury by regulating NF-kappaB and JNK signaling pathways. BMC Complement Med Ther 21(1):284. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s12906-021-03458-0\u003c/span\u003e\u003cspan address=\"10.1186/s12906-021-03458-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao Y, Zhang Z, Li M, Luo J, Chen F, Gong Y, Li Y, Wei Y, Su Y, Kong L (2019) Transcriptomic profiles of 33 opium poppy samples in different tissues, growth phases, and cultivars. Sci Data 6(1):66. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41597-019-0082-x\u003c/span\u003e\u003cspan address=\"10.1038/s41597-019-0082-x\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[{"identity":"742dea68-d79f-44bd-b065-a3b1fe25e394","identifier":"10.13039/501100012165","name":"Key Technologies Research and Development Program","awardNumber":"2020YFA0908000","order_by":0},{"identity":"48fb446e-d56c-4cb0-9855-a3f9b3cf9a3c","identifier":"10.13039/501100001809","name":"National Natural Science Foundation of China","awardNumber":"82011530137","order_by":1},{"identity":"909a580c-15f1-4239-8e08-9ef7551f7613","identifier":"10.13039/501100001809","name":"National Natural Science Foundation of China","awardNumber":"31961133007","order_by":2},{"identity":"55cc0940-7252-46eb-9eab-89e9ea9fe75c","identifier":"10.13039/501100005892","name":"China Academy of Chinese Medical Sciences","awardNumber":"CI2023D002","order_by":3},{"identity":"cab6b96a-2793-4dd6-a657-61c90e96fe95","identifier":"10.13039/501100005892","name":"China Academy of Chinese Medical Sciences","awardNumber":"CI2023E002","order_by":4}],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"China Academy of Chinese Medical Sciences","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"aporphine alkaloids, biosynthetic gene clusters, biosynthesis, CYP80, Stephania yunnanensis, transcriptome","lastPublishedDoi":"10.21203/rs.3.rs-5384973/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5384973/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThe medicinal plant \u003cem\u003eStephania yunnanensis\u003c/em\u003e is rich in aporphine alkaloids, a type of benzylisoquinoline alkaloids (BIAs), with aporphine being the representative and most abundant compound, but our understanding on the biosynthesis of BIA alkaloids in this plant have been relatively limited. Previous research has reported the genome of \u003cem\u003eS. yunnanensis\u003c/em\u003e and preliminarily identified the upstream gene norcoclaurine synthase (NCS) in the BIA biosynthetic pathways. However, the key genes promoting the formation of the aporphine skeleton have not yet been reported. In this study, based on the differences in the content of crebanine and several other BIAs in different tissues, we conducted transcriptome sequencing of roots, stems, and leaves. We then identified candidate genes through functional annotation and sequence alignment, followed by transcriptomic and genomic analyses. Based on this analysis, we identified three CYP80 enzymes (SyCYP80Q5-1, SyCYP80Q5-3, and SyCYP80G6), which exhibited different activities towards (\u003cem\u003eS\u003c/em\u003e)- and (\u003cem\u003eR\u003c/em\u003e)-configured substrates in \u003cem\u003eS. yunnanensis\u003c/em\u003eand demonstrated strict stereoselectivity enroute to aporphine. This study provides metabolomic and transcriptomic information on the biosynthesis of BIAs in \u003cem\u003eS. yunnanensis\u003c/em\u003e and offers valuable insights into the elucidation of BIA biosynthesis, and lays the foundation for the complete analysis of pathways for more aporphine alkaloids.\u003c/p\u003e","manuscriptTitle":"Transcriptome analysis of Stephania yunnanensis and functional validation of CYP80s involved in benzylisoquinoline alkaloids biosynthesis","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-11-05 10:02:00","doi":"10.21203/rs.3.rs-5384973/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"eee6e76e-85c5-4a96-b127-b050997f2504","owner":[],"postedDate":"November 5th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":39763098,"name":"Molecular Biology"},{"id":39763099,"name":"Plant Molecular Biology and Genetics"}],"tags":[],"updatedAt":"2024-11-05T10:02:00+00:00","versionOfRecord":[],"versionCreatedAt":"2024-11-05 10:02:00","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-5384973","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5384973","identity":"rs-5384973","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}