車蘇容 張家源 盧偉 祁克明 魏藝聰

滴要：[目的]從基因表達的水平，初步分析草珊瑚葉和根之間次生代謝差異的分子機制，為二者之間臨床療效差異形成的分子機制分析提供信息。[方法]以福建省福州市的草珊瑚作為樣品，采用lluminaHiSeqTM高通量測序技術測定草珊瑚葉和根的轉錄組，然后經(jīng)過濾和Trinity組裝，得到的unigenes再通過blast與Nr、Nt、Pfam、KOG、Swiss-Prot、Kegg和GO進行比對注釋，并對葉和根的基因差異表達進行分析，尤其是對KEGG代謝通路富集的差異基因進行分析。[結果]轉錄組測序結果共獲得0.4億多個cleanreads，經(jīng)Trinity組裝后共得到508271個unigenes，其平均長度為740bp，最大長度為17.3kb?；赽last分析，共有148561個unigenes在七大功能注釋數(shù)據(jù)庫中得到成功注釋，占總基因數(shù)的58.80%。在分析基因表達水平差異時，發(fā)現(xiàn)草珊瑚葉和根的共同基因有93127個，葉和根的差異基因分別為36327個和52268個;同時還發(fā)現(xiàn)在29732種不同表達的unigenes中，有12511個上調(diào)基因和17221個下調(diào)基因;代謝相關KEGG具有顯著差異的通路有淀粉和蔗糖代謝、苯丙烷類生物合成、乙醛酸和二羧酸代謝、光合生物的固碳作用、吞噬體、谷胱甘肽代謝、光合作用、丙氨酸、天冬氨酸和谷氨酸代謝、倍半萜類和三萜類生物合成、卟啉和葉綠素代謝、氮素代謝、晝夜節(jié)律一植物、光合作用一天線蛋白、芪類、二芳基庚酸和姜酚生物合成、不飽和脂肪酸生物合成、檸檬烯和蒎烯降解、類胡蘿卜素生物合成、二萜類生物合成、類黃酮生物合成、脂肪酸延伸等。其中與藥效密切相關的次生代謝通路苯丙烷類、倍半萜類和三萜類、二萜類、類黃酮類生物合成等途徑分別有193個、82個、40個、35個差異表達基因，而上調(diào)倍半萜合酶、ent-kaur-16-烯合酶、黃酮醇合酶/黃烷酮3-羥化酶等基因和下調(diào)8-羥基香葉醇脫氫酶、vinorine合酶、角鯊烯合酶等關鍵酶基因差異顯著。[結論]草珊瑚葉和根中苯丙烷類、倍半萜類和三萜類、二萜類、類黃酮次生代謝途徑的相關基因差異最為顯著，其中差異顯著的關鍵酶基因可為分析其葉和根之間次生代謝差異的分子機制提供重要信息。

關鍵詞：草珊瑚;葉;根;Illumina測序;基因差異表達分析;次生代謝

中圖分類號：R932

文獻標志碼：A

文章編號：1008-0384（2020）060598-13

[研究意義]草珊瑚，又名腫節(jié)風、九節(jié)風、接骨木、接骨金粟蘭等，是金粟蘭科草珊瑚屬植物草珊瑚Sarcandraglabra（Thunb.）Nakai的干燥全草叫。在2015年版《中國藥典》中的草珊瑚是以腫節(jié)風收載，其性平，味苦、辛，歸心、肝經(jīng)，具有清熱涼血，活血消斑，祛風通絡的功效，用于治療血熱發(fā)斑發(fā)疹，風濕痹痛，跌打損傷。民族用藥中，苗藥關于草珊瑚療效的記載和藥典大致相同，但特別加以記錄了臨床上不同部位的不同功效主治，如：莖、葉主治頭暈、夏季濕病;葉主治骨折，根主治痢疾、胃痛、肝風、跌打腫痛”。目前已從草珊瑚中分離出黃酮類、倍半萜類、香豆素類、有機酸類的活性成分，還得到了較為豐富的揮發(fā)油！、氨基酸和微量元素等成分，現(xiàn)代藥理研究顯示草珊瑚具有抗腫瘤、抗癌、抗菌、抗病毒、消炎鎮(zhèn)痛、對白細胞和血小板的影響、抗胃潰瘍、祛痰平喘、促進骨折愈合、調(diào)整機體免疫的作用因此，若能從分子水平研究草珊瑚不同部位次生代謝調(diào)控基因表達差異的相關信息，將更好地為草珊瑚藥用價值的分析提供線索，對有效利用草珊瑚資源具有重要意義。[前人研究進展]近年來，高通量轉錄組測序技術日益發(fā)展，在藥用植物轉錄組測序方面也已經(jīng)得到廣泛的應用?，F(xiàn)有的技術可以通過測定不同品種、不同采收期、不同產(chǎn)地、不同入藥部位等藥用植物轉錄組信息，整體地研究藥效成分相關酶基因的表達情況，進而發(fā)現(xiàn)其次生代謝產(chǎn)物的合成途徑及其調(diào)控機制凹。唐娟巧對三種不同品種的茯苓進行了轉錄組測序，研究了其基因的差異表達和可能與茯苓多糖含量有關的關鍵基因;鄧楠等叫研究了不同發(fā)芽時期麻黃種子在萌發(fā)時的相關基因及重要基因的表達情況;謝冬梅等5研究了道地藥材丹皮和非道地丹皮之間的基因表達差異，為道地性的形成闡明機制;陳延清等心建立了冬凌草葉和莖的轉錄組數(shù)據(jù)庫;朱孝軒對長春花的葉、花、根的轉錄組數(shù)據(jù)庫進行分析，得出了不同藥用部位中基因的表達差異。[本研究切人點]目前尚未見涉及轉錄組測序法研究草珊瑚不同部位基因差異表達的報道。[擬解決的關鍵問題]本試驗通過lluminaHiSeqTM高通量測序技術來獲得草珊瑚葉和根的轉錄組信息，從基因表達的角度，分析草珊瑚葉和根之間基因表達差異，為二者之間可能存在的次生代謝差異分析提供信息，為草珊瑚葉和根不同藥用價值分析提供參考，有利于草珊瑚葉和根資源的精準利用。

1材料與方法

1.1試驗材料

試驗樣品取自福建省福州市閩侯縣（東經(jīng)1191057.54"，北緯2604'36.10"），經(jīng)福建中醫(yī)藥大學盧偉教授鑒定為金粟蘭科植物草珊瑚Sarcandraglabra（Thunb.）Nakai。分別取林下種植、5年生的同一居群的3株草珊瑚葉和根樣本各3份，編號分別為L1、L2、L3和R1、R2、R3，經(jīng)過清洗、千燥后立即在液氮中冷凍并儲存于-80C備用。

1.2基因文庫的構建與RNA測序

使用TRIzol試劑盒分別提取草珊瑚樣品中葉和根的總RNA，經(jīng)檢測合格后，富集各自的mRNA，并以此模板進一步合成雙鏈cDNA，再經(jīng)純化處理、選擇合適的片段大小，最后通過PCR技術擴增，并純化PCR產(chǎn)物以獲得基因文庫。庫檢合格后，最終通過IlluminaHiSeqTM進行高通量測序。

1.3數(shù)據(jù)的篩選與轉錄本的組裝

使用Trinity軟件8將得到的rawreads進行過濾純化，去除帶接頭的reads、無法確定堿基信息比例大于10%的reads和低質(zhì)量的reads（質(zhì)量值Qphred《=20的堿基數(shù)占整個reads的50%以上的reads），得到cleanreads。再對cleanreads進行組裝，得到轉錄本序列。再使用Corset對轉錄本進行層次聚類，得到最終unigene。

進一步分析這四個代謝途徑中差異倍數(shù)前10名的基因，可以發(fā)現(xiàn)在苯丙烷類代謝途徑上，差異倍數(shù)較大的上調(diào)基因有莽草酸氧羥基肉桂?；D移酶、過氧化物酶、B一葡萄糖苷酶，下調(diào)基因有8-羥基香葉醇脫氫酶、過氧化物酶、甘露醇脫氫酶;在倍半萜類和三萜類代謝途徑上，差異倍數(shù)較大的上調(diào)基因有倍半萜合酶、鍺烷烯-D合酶，下調(diào)基因有（-）-鍺烷烯D合酶、valencene合酶、8-鎘烯合酶同工酶、角鯊烯合酶;在二萜類代謝途徑上，差異倍數(shù)較大的上調(diào)基因有赤霉素2-B-雙加氧酶、ent-kaur-16-烯合酶，下調(diào)基因有en-kaurene氧化酶、戊二烯丙基二磷酸合酶、ent-kaurenoicacidOXidase、赤霉素2-B-雙加氧酶;在類黃酮代謝途徑上，差異倍數(shù)較大的上調(diào)基因有黃酮醇合酶/黃烷酮3-羥化酶、莽草酸氧-羥基肉桂?；D移酶、查爾酮合成酶、黃烷酮3-羥化酶，下調(diào)基因有vinorine合酶、咖啡酰輔酶AO-甲基轉移酶、4-香豆酸-3-羥化酶、莽草酸氧-羥基肉桂酰基轉移酶、咖啡酰輔酶A3-0-甲基轉移酶。結果如表6～9。

3討論與結論

在本試驗中，我們采用IluminaHiSeq”^高通量測序技術對草珊瑚葉和根的轉錄組進行測序，去除冗余后各獲得了0.4億多cleanreads，其中，各樣品的Q20、Q30值均大于97.60%、93.10%，GC量均大于45.20%，說明建立的整個文庫質(zhì)量較高，再將這些讀數(shù)通過使用Trinity組裝成508271個轉錄物，其中Transcripts平均長度為740bp，最大長度為17.3kb，N50=l178bp，說明組裝的轉錄本質(zhì)量較高，測序質(zhì)量良好?；赽last分析，共有148561個unigenes可以在功能注釋數(shù)據(jù)庫（Nr、Nt、Pfam、KOG、Swiss-Prot、Kegg和GO），上匹配，占總基因數(shù)的58.80%，剩余的unigenes由于與未知功能的蛋白質(zhì)匹配，可能被認為是新的轉錄物或替代體，這些unigenes有待進一步的研究。

在分析基因表達水平差異時，采用DESeq2篩選出合格的readcount數(shù)據(jù)，發(fā)現(xiàn)草珊瑚葉和根的共同基因有93127個，葉和根的差異基因分別為36327個和52268個，同時還發(fā)現(xiàn)在29732種不同表達的unigenes中，有12511個上調(diào)基因和17221個下調(diào)基因，在葉中被大幅度上調(diào)，與根表達差異倍數(shù)最大的基因有抗壞血酸過氧化酶、植酰苯醌甲基轉移酶、腺苷硫酸鹽還原酶等;在根中大幅度下調(diào)，與葉中表達差異倍數(shù)最大的基因有葡萄糖脫氫酶、類枯草素蛋白酶、谷胱甘肽轉移酶等，這與Mizra-chi等人2的研究結果相似。

在差異基因KEGG顯著性富集通路中，與藥效密切相關的次生代謝通路包括苯丙烷類生物合成基因193個，其中差異倍數(shù)較大的上調(diào)基因有莽草酸氧一羥基肉桂?；D移酶、過氧化物酶、B-葡萄糖苷酶，下調(diào)基因有8-羥基香葉醇脫氫酶、過氧化物酶、甘露醇脫氫酶;倍半萜類和三萜類生物合成基因82個，其中差異倍數(shù)較大的上調(diào)基因有倍半萜合酶、鍺烷烯-D合酶，下調(diào)基因有（一）一鍺烷烯D合酶、valencene合酶、8-鎘烯合酶同工酶、角鯊烯合酶：萜類生物合成基因40個，其中差異倍數(shù)較大的上調(diào)基因有赤霉素2-B-雙加氧酶、ent-kaur-16-烯合酶，下調(diào)基因有en-kaurene氧化酶、戊二烯丙基二磷酸合酶、ent-kaurenoicacidoxidase、赤霉素2-B-雙加氧酶;類黃酮生物合成基因35個，其中差異倍數(shù)較大的上調(diào)基因有黃酮醇合酶/黃烷酮3-羥化酶、莽草酸氧羥基肉桂?；D移酶、查爾酮合成酶、黃烷酮3-羥化酶，下調(diào)基因有vinorine合酶、咖啡酰輔酶AO-甲基轉移酶、4-香豆酸-3-羥化酶、莽草酸氧-羥基肉桂?；D移酶、咖啡酰輔酶A3-0-甲基轉移酶。

對于苯丙烷類生物合成代謝，穆紅梅等巧總結了植物中苯丙烷類物質(zhì)生物合成的主要途徑，其下游分支有形成木質(zhì)素類的成分、形成花青素、單寧等物質(zhì)和形成黃酮類成分（包括查爾酮、黃烷酮、總黃酮）。其中，黃酮類物質(zhì)作為草珊瑚的主要化學成分，具有抗癌、抗腫瘤、消炎抗菌、抗氧化、降血脂的功效。王敦清等26研究發(fā)現(xiàn)，在草珊瑚不同部位中，葉的總黃酮含量相對較高，本試驗的結果不僅從基因的角度證明了王敦清等的研究，還進一步發(fā)現(xiàn)在葉中差異最大的基因是黃酮醇合酶/黃烷酮3-羥化酶基因，該類基因顯著上調(diào)，其下游產(chǎn)物：高良姜素、楊梅素、山奈素、槲皮素的含量顯著增多。因此，在臨床上治療癌定、腫瘤、心血管、肝臟、神經(jīng)系統(tǒng)等疾病時往往可以優(yōu)先選擇草珊瑚的葉作為主要原料部位。在苯丙烷類物質(zhì)中也包括木質(zhì)素類（香豆素類）成分，如現(xiàn)行的藥典規(guī)定草珊瑚含量測定的指標性成分之一的異嗪皮啶！。包俠萍56發(fā)現(xiàn)，草珊瑚中異嗪皮啶的含量在根中最高，其次是莖、葉中最低;本試驗發(fā)現(xiàn)，在根中差異倍數(shù)最大的基因是8-羥基香葉醇脫氫酶，其下游產(chǎn)物芥子醇參與木質(zhì)素類生物合成途徑”，故進一步佐證了合成異嗪皮啶成分的基因在葉和根中的差異較為顯著，認為差異與8～羥基香葉醇脫氫酶基因顯著上調(diào)有關。而對于植物萜類物質(zhì)生物合成代謝，近年來已有多數(shù)研究學者對其進行詳盡的歸納總結8-39，藥典中規(guī)定草珊瑚以全草人藥，其化學成分中含量最豐富的是倍半萜類的成分，在差異基因KEGG富集通路排名靠前的結果中，我們可以發(fā)現(xiàn)草珊瑚葉和根在二萜類生物合成代謝上差異較大，以根中的表達較高，差異最大的基因是en-kaurene氧化酶，其下游產(chǎn)物是ent-Kaur-16-en-19-oate，該基因的下調(diào)可能是二者療效差異的指標之一;而倍半萜類和三萜類生物合成代謝上基因差異性不是很顯著，差異基因倍半萜合酶和（-）-鍺烷烯D合酶的下游產(chǎn)物均為（-）-鍺烷烯，說明葉和根倍半萜類化合物含量相差不大，以葉中含量較高，并且僅在根中發(fā)現(xiàn)角鯊烯合酶基因表現(xiàn)為上周，說明根還具有一定的抗氧化、保肝、促進血液循環(huán)等作用，對治療肝和心血管方面的疾病具有優(yōu)勢。因此，本試驗可以表明草珊瑚的葉和根之間基因表達的顯著差異，主要表現(xiàn)在苯丙烷類、倍半萜類和三萜類、二萜類、類黃酮生物合成次生代謝途徑的差異。其中上調(diào)倍半萜合酶、ent-kaur-16-烯合酶、黃酮醇合酶/黃烷酮3～羥化酶等基因和下調(diào)8-羥基香葉醇脫氫酶、vinorine合酶、角鯊烯合酶等基因，是葉和根療效差異的主要分子機制。

由于前人未對草珊瑚開展過全面的基因組測序，可供參考的遺傳信息較少，在一定程度上限制了草珊瑚功能基因的發(fā)掘和研究。本研究獲得了草珊瑚葉和根的轉錄組信息，通過比較不同部位的差異基因表達，尤其是葉和根在次生代謝中相關的基因差異性，可為進一步分析草珊瑚葉和根之間可能的次生代謝差異分析提供基因表達方面的信息，為草珊瑚葉和根不同藥用價值分析提供了重要信息與線索，有利于提升草珊瑚葉和根資源的精準利用。

參考文獻：

[1]福建省科學技術委員會《福建植物志》編寫組.福建植物志：第四卷[M].福州：福建科學技術出版社，1990.

[2]中國科學院《中國植物志》編輯委員會.中國植物志[M].北京：科學出版社，2005：287-290.

[3]國家藥典委員會.中華人民共和國藥典：2015年版一部[M].北京：中國醫(yī)藥科技出版社，2015：234-235.

[4]賈敏如，李星煒.中國民族藥志要[M].北京：中國醫(yī)藥科技出版社，2005：365366.

[5]楊榮平，王賓豪，勵娜，等.GC-MS法分析腫節(jié)風葉中揮發(fā)油化學成分[J].中成藥，2008，30（11）：1703-1704.

[6]徐麗麗.全緣金粟蘭及草珊瑚的化學成分研究[D].昆明：云南中醫(yī)學院，2016.

[7]徐艷琴，劉小麗，黃小方，等.草珊瑚的研究現(xiàn)狀與展望[J].中草藥，

2011，42（12）：2552-2559.

[8]周斌，劉可越，常軍，等.中藥腫節(jié)風的化學成分和藥理作用研究進展[J].中國現(xiàn)代應用藥學，2009，26（12）：982-986.

[9]姜伶，李景輝.中藥腫節(jié)風的抗腫瘤作用研究進展[].中國執(zhí)業(yè)藥師，2014，11（04）：29-31，35.

[10]梅全喜，胡瑩.腫節(jié)風的藥理作用及臨床應用研究進展[].時珍國醫(yī)國藥，2011，22（1）：230-232.

[11]孫洪計，魏慧君.RNA一Seq技術在轉錄組研究中的應用].中外醫(yī)學研究，2018，16（20）：184-187.

[12]王媛媛，楊美青.藥用植物轉錄組的研究進展[J].安徽農(nóng)學通報，2019，25（8）：13-15，52.

[13]唐娟.三個茯苓品種的品質(zhì)特性及轉錄組分析[D].長沙：湖南農(nóng)業(yè)大學，2016.

[14]鄧楠，史勝青，常二梅，等.膜果麻黃種子不同發(fā)育時期的轉錄組測序分析[J].東北林業(yè)大學學報，2015，43（2）：28-32.

[15]謝冬梅，俞年軍，黃璐琦，等.基于高通量測序的藥用植物“鳳丹根皮的轉錄組分析[J].中國中藥雜志，2017，42（15）：2954-2961.

[16]陳延清，胡志剛，劉大會，等.藥用植物冬凌草高通量轉錄組測序與分析[J].中國現(xiàn)代中藥，2018，20（12）：1476-1482.

[17]朱孝軒.長春花轉錄組與萜類吲哚生物堿代謝途徑研究[D].北京：北京協(xié)和醫(yī)學院，2015.

[18] GRABHERR M G， HAAS B J， YASSOUR M， et al. Full-lengthtranscriptome assembly from RNA-Seq data without a reference genome [J]. Nature Biotechnology， 2011，29 （7） ： 644-652.

[19] GOTZ S， GARCIA-GOMEZ J M， TEROL J， et al. High-throughputfunctional annotation and data mining with the Blast2GO suite [J]. Nucleic Acids Research， 2008， 36 （10） ： 3420-3435.

[20] LOVE M I， HUBER W， ANDERS S. Moderated estimation of foldchange and dispersion for RNA-seq data with DESeq2 [J]. Genome Biology， 2014， 15 （12） ： 550.

[21] YOUNG M D， WAKEFIELD MJ， SMYTH G K， et al. Gene ontologyanalysis for RNA-seq： accounting for selection bias [J]. Genome Biology， 2010， 11 （2） ：R14.

[22] MIZRACHI E， HEFER C A， RANIK M， et al. De novo assembledexpressed gene catalog of a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq [J]. BMC Genomics， 2010， 11 （1） ： 681.

[23]穆紅梅，杜秀菊，張秀省，等.植物MYB轉錄因子調(diào)控苯丙烷類生物合成研究[J].北方園藝，2015（24）：171-174.

[24]邵佳.草珊瑚總黃酮提取純化及藥理作用研究[D].貴陽：貴州大學，2008.

[25]鄭永標，許小萍，鄒先文，等.草珊瑚藥材抗氧化活性化學成分研究[J].福建師范大學學報（自然科學版），2016，32（3）：98-102.

[26]王敦清，李先春.草珊瑚根莖葉中總黃酮成分的研究[J].中草藥，1996，27（6）：337-338.

[27]陳鄭，哈文波.高良姜素抗腫瘤作用機制的研究進展[J].醫(yī)學綜述，2017，23（9）：1752-1756.

[28]張旭光，尹航，陳峰，等.高良姜素藥理活性的研究進展[J].中國現(xiàn)代中藥，2016，18（11）：1532-1536.

[29]吳少花，陳君，劉亞萌，等.楊梅素抗腫瘤活性研究進展[].吉林醫(yī)藥學院學報，2015，36（5）：381-383.

[30]王潞，周云英.楊梅素抗感染、抗炎及抗氧化活性研究進展[].中草藥，2019，50（3）：778-784.

[31]林國鋇，謝燕，李國文.楊梅素的研究進展[J].國際藥學研究雜志，2012，39（6）：483-487.

[32]翟廣玉，馬海英，郜蕾.槲皮素及其衍生物的抗腫瘤活性研究進展[J].化學試劑，2015，37（2）：97-103.

[33]陳振華，胡曉艷，趙滕，等.槲皮素對心血管系統(tǒng)疾病的影響及其新劑型研究進展].時珍國醫(yī)國藥，2019，30（2）：440-443.

[34]張超，李昌平.槲皮素治療非酒精性脂肪性肝病的作用機制研究進展[J].中藥新藥與臨床藥理，2015，26（5）：718-721.

[35]石玥，梁曉春.槲皮素防治神經(jīng)退行性疾病的機制研究進展[J].中國中西醫(yī)結合雜志，2012，32（10）：1432-1435.

[36]包俠萍.草珊瑚不同部位及不同采收期有效成分含量的考察[J].海峽藥學，2014，26（12）：42-44.

[37]馮鶴翔，涂軼.木質(zhì)素生物合成的研究[J].青島大學學報（自然科學版）2018，31（1）：46-54.

[38]時敏，王瑤，周偉，等.藥用植物萜類化合物的生物合成與代謝調(diào)控研究進展[J].中國科學：生命科學，2018，48（4）：352-364.

[39]蘇文炳，蔣園園，白昀鷺，等.轉錄因子調(diào)控植物萜類化合物生物合成研究進展[J].農(nóng)業(yè)生物技術學報，2019，27（5）：919-926.

[40]金祖漢，金捷，毛培江，等.角鯊烯對四氯化碳急性肝損傷模型小鼠的抗氧化和護肝作用研究[J].浙江中醫(yī)藥大學學報，2015，39（9）：666-670.

[41]繆云萍，陳愛瑛，夏志國，等.角鯊烯對小鼠急性酒精性肝損傷的保護作用[J].食品工業(yè)科技，2015，36（16）：364-365，377.

[42] KATSELOU M G， MATRALIS A N， KOUROUNAKIS A P.Developing potential agents against atherosclerosis： Design， synthesis and pharmacological evaluation of novel dual inhibitors of oxidative stress and Squalene Synthase activity [J]. European Journal of Medicinal Chemistry， 2017， 138： 748-760.

（責任編輯：吳宇琳）

車蘇容，張家源，盧偉，等.轉錄組測序法研究草珊瑚葉和根的基因差異表達[J].福建農(nóng)業(yè)學報，2020，35（6）：598-610.