The near-complete genome assembly of Reynoutria multiflora reveals the genetic basis of stilbene and anthraquinone biosynthesis

Abstract

Reynoutria multiflora is a widely used medicinal plant in China. Its medicinal compounds are mainly stilbenes and anthraquinones which possess important pharmacological activities in anti-aging, anti-inflammatory and anti-oxidation, but their biosynthetic pathways are still largely unresolved. Here, we reported a near-complete genome assembly of R. multiflora consisting of 1.39 Gb with a contig N50 of 122.91 Mb and only one gap left. Genome evolution analysis revealed that two recent bursts of long terminal repeats (LTRs) contributed significantly to the increased genome size of R. multiflora, and numerous large chromosome rearrangements were observed between R. multiflora and Fagopyrum tataricum genomes. Comparative genomics analysis revealed that a recent whole-genome duplication specific to Polygonaceae led to a significant expansion of gene families associated with disease tolerance and the biosynthesis of stilbenes and anthraquinones in R. multiflora. Combining transcriptomic and metabolomic analyses, we elucidated the molecular mechanisms underlying the dynamic changes in content of medicinal ingredients in R. multiflora roots across different growth years. Additionally, we identified several putative key genes responsible for anthraquinone and stilbene biosynthesis. We identified a stilbene synthase gene PM0G05131 highly expressed in roost, which may exhibit an important role in the accumulation of stilbenes in R. multiflora. These genomic data will expedite the discovery of anthraquinone and stilbenes biosynthesis pathways in medicinal plants.

Keywords

comparative genomics / genome assembly / Reynoutria multiflora / stilbene and anthraquinone biosynthesis

References

[1] AramakiT,Blanc-Mathieu R,EndoH,OhkuboK,Kanehisa M,GotoS,OgataH. 2020. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics 36:2251–2252.

[2] BaiC,YangJ, CaoB,XueY, GaoP,LiangH, LiG. 2020. Growth years and post-harvest processing methods have critical roles on the contents of medicinal active ingredients of Scutellaria baicalensis. Industrial Crops and Products 158:112985.

[3] BinghangL,ShiY, YuanJ,Galaxy Y,ZhangH,LiN,LiZ, ChenY,Fan W. 2013. Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects. Quantitative Biology 35:62–67.

[4] CarboneroP,Iglesias-Fernández R,Vicente-CarbajosaJ. 2017. The AFL subfamily of B3 transcription factors: Evolution and function in angiosperm seeds. Journal of Experimental Botany 68:871–880.

[5] ChanY-C,ChengF-C, WangM-F. 2002. Beneficial effects of different Polygonum multiflorum Thunb. extracts on memory and hippocampus morphology. Journal of Nutritional Science and Vitaminology 48:491–497.

[6] ChangY-M,LinH-H, LiuW-Y,Yu C-P,ChenH-J,WartiniPP,KaoY-Y, WuY-H,Lin J-J,LuM-YJ,TuS-L,WuS-H, ShiuS-H,Ku MSB,LiW-H 2019. Comparative transcriptomics method to infer gene coexpression networks and its applications to maize and rice leaf transcriptomes. Proceedings of the National Academy of Sciences of the United States of America 116:3091–3099.

[7] ChenW,GongL, GuoZ,WangW, ZhangH,Liu X,YuS,XiongL,LuoJ. 2013. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: Application in the study of rice metabolomics. Molecular Plant 6:1769–1780.

[8] ChengH,Concepcion GT,FengX,ZhangH,LiH. 2021. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nature Methods 18:170–175.

[9] ChengY,LiuH, TongX,Liu Z,ZhangX,LiD,JiangX, YuX. 2020. Identification and analysis of CYP450 and UGT supergene family members from the transcriptome of Aralia elata (Miq.) seem reveal candidate genes for triterpenoid saponin biosynthesis. BMC Plant Biology 20:214.

[10] DoleželJ,Greilhuber J,SudaJ. 2007. Estimation of nuclear DNA content in plants using flow cytometry. Nature Protocols 2:2233–2244.

[11] DudchenkoO,BatraSS, OmerAD,Nyquist SK,HoegerM,DurandNC,ShamimMS, MacholI,Lander ES,AidenAP,AidenEL. 2017. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science 356:92–95.

[12] DurandNC,ShamimMS, MacholI,Rao SSP,HuntleyMH,LanderES,AidenEL. 2016. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Systems 3:95–98.

[13] EddyS. 2009. A new generation of homology search tools based on probabilistic inference. Genome Informatics. International Conference on Genome Informatics 23:205–211.

[14] EllinghausD,KurtzS, WillhoeftU. 2008. LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons. BMC Bioinformatics 9:18.

[15] EmmsDM,KellyS. 2019. OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biology 20:238.

[16] GoschC,Halbwirth H,SchneiderB,HölscherD,StichK. 2010. Cloning and heterologous expression of glycosyltransferases from Malus x domestica and Pyrus communis, which convert phloretin to phloretin 2’-O-glucoside (phloridzin). Plant Science 178:299–306.

[17] GremmeG,BrendelV, SparksME,Kurtz S. 2005. Engineering a software tool for gene structure prediction in higher organisms. Information and Software Technology 47:965–978.

[18] HaasBJ,Papanicolaou A,YassourM,GrabherrM,BloodPD, BowdenJ,Couger MB,EcclesD,LiB,LieberM, MacmanesMD,Ott M,OrvisJ,PochetN,StrozziF, WeeksN,Westerman R,WilliamT,DeweyCN,Henschel R,LeDucRD,FriedmanN,RegevA. 2013. De novo transcript sequence reconstruction from RNA-Seq using the Trinity platform for reference generation and analysis. Nature protocols 8:1494–1512.

[19] HaasBJ,Salzberg SL,ZhuW,PerteaM,AllenJE, OrvisJ,White O,BuellCR,WortmanJR. 2008. Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments. Genome Biology 9:R7.

[20] HanMV,ThomasGWC, Lugo-MartinezJ,HahnMW. 2013. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Molecular Biology and Evolution 30:1987–1997.

[21] HanY,ChenP, ChenY. 2021. Prospects for planting Polygonum multiflorum under forests in Hainan. Forest Science and Technology 5:66–68.

[22] HeM,HeY, ZhangK,Lu X,ZhangX,GaoB,FanY, ZhaoH,Jha R,HudaMN,TangY,WangJ, YangW,Yan M,ChengJ,RuanJ,DullooE, ZhangZ,Georgiev MI,ChapmanMA,ZhouM. 2022. Comparison of buckwheat genomes reveals the genetic basis of metabolomic divergence and ecotype differentiation. New Phytologist 235:1927–1943.

[23] HeQ,MaD, LiW,XingL, ZhangH,Wang Y,DuC,LiX,JiaZ, LiX,LiuJ, LiuZ,MiaoY, FengR,Lv Y,WangM,LuH,LiX, XiaoY,Wang R,LiangH,ZhouQ,ZhangL, LiangC,Du H. 2023. High-quality Fagopyrum esculentum genome provides insights into the flavonoid accumulation among different tissues and self-incompatibility. Journal of Integrative Plant Biology 65:1423–1441.

[24] IsahT. 2019. Stress and defense responses in plant secondary metabolites production. Biological Research 52:39.

[25] IyerLM,KooninEV, AravindL. 2001. Adaptations of the helix-grip fold for ligand binding and catalysis in the START domain superfamily. Proteins: Structure, Function, and Genetics 43:134–144.

[26] KangM,FuR, ZhangP,Lou S,YangX,ChenY,MaT, ZhangY,Xi Z,LiuJ. 2021. A chromosome-level Camptotheca acuminata genome assembly provides insights into the evolutionary origin of camptothecin biosynthesis. Nature Communications 12:3531.

[27] KangS-H,PandeyRP, LeeC-M,Sim J-S,JeongJ-T,ChoiB-S,JungM, GinzburgD,Zhao K,WonSY,OhT-J,YuY, KimN-H,Lee OR,LeeT-H,BashyalP,KimT-S, LeeW-H,Hawkins C,KimC-K,KimJS,AhnBO, RheeSY,Sohng JK. 2021. Genome-enabled discovery of anthraquinone biosynthesis in Senna tora. Nature Communication 11:5875.

[28] KarppinenK,Hokkanen J,MattilaS,NeubauerP,HohtolaA. 2008. Octaketide-producing type III polyketide synthase from Hypericum perforatum is expressed in dark glands accumulating hypericins. The FEBS Journal 275:4329–4342.

[29] KimCK,KimYK. 2018. The multipartite mitochondrial genome of Fallopia multiflora (Caryophyllales: Polygonaceae). Mitochondrial DNA Part B 3:155–156.

[30] KimD,PaggiJM, ParkC,Bennett C,SalzbergSL. 2019. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nature Biotechnology 37:907–915.

[31] LangfelderP,HorvathS. 2009. WGCNA: An R package for weighted correlation network analysis. BMC Bioinformatics 9:559.

[32] LeeJ,WaminalNE, ChoiH-I,Perumal S,LeeS-C,NguyenVB,JangW, KimN-H,Gao L,YangT-J. 2017. Rapid amplification of four retrotransposon families promoted speciation and genome size expansion in the genus Panax. Scientific Reports 7:9045.

[33] LeeW,KuS-K, BaeJ-S. 2013. Emodin-6-O-β-d-glucoside down-regulates endothelial protein C receptor shedding. Archives of Pharmacal Research 36:1160–1165.

[34] LiG,WangL, YangJ,He H,JinH,LiX,RenT, RenZ,LiF, HanX,ZhaoX, DongL,Li Y,SongZ,YanZ,ZhengN, ShiC,WangZ, YangS,Xiong Z,ZhangM,SunG,ZhengX, GouM,JiC, DuJ,ZhengH, DoleželJ,DengXW,SteinN, YangQ,Zhang K,WangD. 2021. A high-quality genome assembly highlights rye genomic characteristics and agronomically important genes. Nature Genetics 53:574–584.

[35] LiH,DurbinR. 2009. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics 25:1754–1760.

[36] LiHL,WuL, DongZ,Jiang Y,JiangS,XingH,LiQ, LiuG,TianS, WuZ,WuB, LiZ,ZhaoP, ZhangY,Tang J,XuJ,HuangK,LiuX, ZhangW,Liao Q,RenY,HuangX,LiQ, LiC,WangY, Xavier-RaviB,LiH,LiuY, WanT,LiuQ, ZouY,JianJ, XiaQ,LiuY. 2021. Haplotype-resolved genome of diploid ginger (Zingiber officinale) and its unique gingerol biosynthetic pathway. Horticulture Research 8:189.

[37] LiP,QuanX, JiaG,XiaoJ, CloutierS,You FM. 2016. RGAugury: A pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants. BMC Genomics 17:852.

[38] LiY,NiuZ, ZhuM,WangZ, XuR,LiM, ZhengZ,Lu Z,DongC,HuH,YangY, WuY,WangD, YangJ,Zhang J,WanD,AbbottR,LiuJ, YangY. 2023. Multi-omics data provide insight into the adaptation of the glasshouse plant Rheum nobile to the alpine subnival zone. Communications Biology 6(1):906.

[39] LinL,NiB, LinH,ZhangM, LiX,YinX, QuC,NiJ. 2015. Traditional usages, botany, phytochemistry, pharmacology and toxicology of Polygonum multiflorum Thunb.: A review. Journal of Ethnopharmacology 159:158–183.

[40] LiuCJ,HuhmanD, SumnerLW,Dixon RA. 2003. Regiospecific hydroxylation of isoflavones by cytochrome p450 81E enzymes from Medicago truncatula. The Plant Journal 36:471–484.

[41] LiuJ,LiuY, WangY,Abozeid A,ZuYG,ZhangXN,TangZH. 2017. GC-MS metabolomic analysis to reveal the metabolites and biological pathways involved in the developmental stages and tissue response of Panax ginseng. Molecules 22:496.

[42] LiuY,WangB, ShuS,LiZ, SongC,Liu D,NiuY,LiuJ,ZhangJ, LiuH,HuZ, HuangB,Liu X,LiuW,JiangL,AlamiMM, ZhouY,Ma Y,HeX,YangY,ZhangT, HuH,BarkerMS, ChenS,Wang X,NieJ. 2021. Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids. Nature Communications 12:3276.

[43] LvL,GuX, TangJ,Ho C. 2007. Antioxidant activity of stilbene glycoside from Polygonum multiflorum Thunb in vivo. Food Chemistry 104:1678–1681.

[44] MizuuchiY,ShiS-P, WanibuchiK,Kojima A,MoritaH,NoguchiH,AbeI. 2009. Novel type III polyketide synthases from Aloe arborescens. The FEBS Journal 276:2391–2401.

[45] NawrockiEP,KolbeDL, EddySR. 2009. Infernal 1.0: Inference of RNA alignments. Bioinformatics 25:1335–1337.

[46] OuS,ChenJ, JiangN. 2018. Assessing genome assembly quality using the LTR Assembly Index (LAI). Nucleic Acids Research 46:e126.

[47] OuS,JiangN. 2018. LTR_retriever: A highly accurate and sensitive program for identification of long terminal repeat retrotransposons. Plant Physiology 176:1410–1422.

[48] OudelaarAM,HiggsDR. 2021. The relationship between genome structure and function. Nature Reviews Genetics 22:154–168.

[49] ParajuliP,PandeyRP, NguyenTHT,Dhakal D,SohngJK. 2018. Substrate scope of o-methyltransferase from Streptomyces peucetius for biosynthesis of diverse natural products methoxides. Applied Biochemistry and Biotechnology 184:1404–1420.

[50] PeninAA,Kasianov AS,KlepikovaAV,KirovIV,Gerasimov ES,FesenkoAN,LogachevaMD. 2021. High-resolution transcriptome atlas and improved genome assembly of common buckwheat,Fagopyrum esculentum. Frontiers in Plant Science 12:612382.

[51] PerteaM,PerteaGM, AntonescuCM,Chang T-C,MendellJT,SalzbergSL. 2015. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nature Biotechnology 33:290–295.

[52] ServantN,Varoquaux N,LajoieBR,ViaraE,ChenC-J, VertJ-P,Heard E,DekkerJ,BarillotE. 2015. HiC-Pro: An optimized and flexible pipeline for Hi-C data processing. Genome Biology 16:259.

[53] SettS,PrasadA, PrasadM. 2022. Resistance genes on the verge of plant-virus interaction. Trends in Plant Science 27:1242–1252.

[54] ShannonP,MarkielA, OzierO,Baliga NS,WangJT,RamageD,AminN, SchwikowskiB,IdekerT. 2003. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Research 13:2498–2504.

[55] ShengS-J,LiuZ-Y, ZhaoW,Shao L,ZhaoS-J. 2010. Molecular analysis of a type III polyketide synthase gene in Fallopia multiflora. Biologia 65:939–946.

[56] SinghA,PanwarR, MittalP,Hassan MI,SinghIK. 2021. Plant cytochrome P450s: Role in stress tolerance and potential applications for human welfare. International Journal of Biological Macromolecules 184:874–886.

[57] SongT,YaoM, YangY,Liu Z,ZhangL,LiW. 2023. Integrative identification by Hi-C revealed distinct advanced structural variations in lung adenocarcinoma tissue. Phenomics 3:390–407.

[58] StamatakisA. 2014. RAxML Version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30:1312–1313.

[59] StankeM,WaackS. 2003. Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics 19:ii215–ii225.

[60] SteeleML,TruongJ, GovindaraghavanS,OoiL,SucherNJ, MünchG. 2013. Cytoprotective properties of traditional Chinese medicinal herbal extracts in hydrogen peroxide challenged human U373 astroglia cells. Neurochemistry International 62:522–529.

[61] SunXQ,BaiMM, YaoH,GuoJL, LiMM,HangYY. 2013. DNA barcoding of populations of Fallopia multiflora, an indigenous herb in China. Genetics and Molecular Research 12:4078–4089.

[62] TangY-C,LiuY-J, HeG-R,Cao Y-W,BiM-M,SongM,YangP-P, XuL-F,Ming J. 2021. Comprehensive analysis of secondary metabolites in the extracts from different lily bulbs and their antioxidant ability. Antioxidants 10:1634.

[63] Tarailo-GraovacM,Chen N. 2009. Using RepeatMasker to identify repetitive elements in genomic sequences. Current Protocols in Bioinformatics unit:4.10.1-4.10.14.

[64] TengHM,FangMF, CaiX,HuZH. 2009. Localization and dynamic change of saponin in vegetative organs of Polygala tenuifolia. Journal of Integrative Plant Biology 51:529–536.

[65] Van de PeerY,MaereS, MeyerA. 2009. The evolutionary significance of ancient genome duplications. Nature Reviews Genetics 10:725–732.

[66] WangJ,XuS, MeiY,CaiS, GuY,SunM, LiangZ,Xiao Y,ZhangM,YangS. 2021. A high-quality genome assembly of Morinda officinalis, a famous native southern herb in the Lingnan region of southern China. Horticulture Research 8:135.

[67] WangL,WangY. 2019. Transcription factor VqERF114 regulates stilbene synthesis in Chinese wild Vitis quinquangularis by interacting with VqMYB35. Plant Cell Reports 38:1347–1360.

[68] WangM,ZhaoR, WangW,Mao X,YuJ. 2012. Lipid regulation effects of Polygoni Multiflori Radix, its processed products and its major substances on steatosis human liver cell line L02. Journal of Ethnopharmacology 139:287–293.

[69] WangY,TangH, DebarryJD,Tan X,LiJ,WangX,LeeT, JinH,MarlerB, GuoH,Kissinger JC,PatersonAH. 2012. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Research 40:e49.

[70] WatanabeK,Praseuth AP,WangCC. 2007. A comprehensive and engaging overview of the type III family of polyketide synthases. Current Opinion in Chemical Biology 11:279–286.

[71] XiaW,RuiW, ZhaoW,Sheng S,LeiL,FengY,ZhaoS. 2018. Stable isotope labeling and 2, 3, 5, 4’-tetrahydroxystilbene-2-O-β-D-glucopyranoside biosynthetic pathway characterization in Fallopia multiflora. Planta 247:613–623.

[72] XuZ,WangH. 2007. LTR_FINDER: An efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Research 35:W265–W268.

[73] YangX. 2005. The effects of refined Polygonum multiflorum Thunb. polysaccharide on learning memory and the activities of enzymes in the brain for the experimental mice with dementia. Progress in Pharmaceutical Sciences 29:557–559.

[74] YinQ,HanX, HanZ,ChenQ, ShiY,GaoH, ZhangT,Dong G,XiongC,SongC,SunW, ChenS. 2020. Genome-wide analyses reveals a glucosyltransferase involved in rutin and emodin glucoside biosynthesis in tartary buckwheat. Food Chemistry 318:126478.

[75] YuG,WangL-G, HanY,HeQ-Y. 2012. clusterProfiler: An R package for comparing biological themes among gene clusters. Omics: A Journal of Integrative Biology 16:284–287.

[76] YuanZ,DongF, PangZ,Fallah N,ZhouY,LiZ,HuC. 2022. Integrated metabolomics and transcriptome analyses unveil pathways involved in sugar content and rind color of two sugarcane varieties. Frontiers in Plant Science 13:921536.

[77] ZhangD,LiW, XiaE,ZhangQ, LiuY,ZhangY, TongY,Zhao Y,NiuY,XuJ,GaoL. 2017. The medicinal herb Panax notoginseng genome provides insights into ginsenoside biosynthesis and genome evolution. Molecular Plant 10:903–907.

[78] ZhangL,LiX, MaB,GaoQ, DuH,HanY, LiY,CaoY, QiM,ZhuY, LuH,MaM, LiuL,ZhouJ, NanC,QinY, WangJ,Cui L,LiuH,LiangC,QiaoZ. 2017. The tartary buckwheat genome provides insights into rutin biosynthesis and abiotic stress tolerance. Molecular Plant 10:1224–1237.

[79] ZhangLS,ChenCL, MaoXX,Liu JQ. 2022. Phylogenomics and evolutionary diversification of the subfamily Polygonoideae. Journal of Systematics and Evolution 61:587–598.

[80] ZhangY,ZhengL, ZhengY,Zhou C,HuangP,XiaoX,ZhaoY, HaoX,HuZ, ChenQ,Li H,WangX,FukushimaK,WangG, LiC. 2019. Assembly and annotation of a draft genome of the medicinal plant Polygonum cuspidatum. Frontiers in Plant Science 10:1274.

[81] ZhaoW,XiaW, LiJ,ShengS, LeiL,ZhaoS. 2014. Transcriptome profiling and digital gene expression analysis of Fallopia multiflora to discover putative genes involved in the biosynthesis of 2, 3, 5, 4’-tetrahydroxy stilbene-2-O-β-D-glucoside. Gene 547:126–135.

[82] ZhaoY,ChuS, GuiS,QinY, XuR,ShanT, PengH. 2021. Tissue-specific metabolite profiling of Fallopia multiflora (Heshouwu) and Fallopia multiflora var. angulata by mass spectrometry imaging and laser microdissection combined with UPLC-Q/TOF-MS. Journal of Pharmaceutical and Biomedical Analysis 200:114070.

[83] ZhengY,JiaoC, SunH,RosliHG, PomboMA,Zhang P,BanfM,DaiX,MartinGB, GiovannoniJJ,ZhaoPX,RheeSY, FeiZ. 2016. iTAK: A program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Molecular Plant 9:1667–1670.

[84] ZhongJ,LiuS-J, MaSS,YangW, HuYL,WuQ, LinZP. 2004. Effect of matrix attachment regions on resveratrol production in tobacco with transgene of stilbene synthase from Parthenocissus henryana. Journal of Integrative Plant Biology 46:948–954.

[85] ZhuZ-J,SchultzAW, WangJ,Johnson CH,YannoneSM,PattiGJ,SiuzdakG. 2013. Liquid chromatography quadrupole time-of-flight mass spectrometry characterization of metabolites guided by the METLIN database. Nature Protocols 8:451–460.

RIGHTS & PERMISSIONS

2024 The Authors. Journal of Systematics and Evolution published by John Wiley & Sons Australia, Ltd on behalf of Institute of Botany, Chinese Academy of Sciences.