王丹丹,冷春玲
(遼東學院 農學院,遼寧丹東 118003)
石柱參植物分類學上隸屬于五加科(Araliaceae)人參屬(Panax)植物人參(PanaxginsengC.A.Mey.),以其干燥的根或根莖入藥;原產于遼寧省丹東市寬甸縣振江鄉石柱子村,以其蘆長、體靈、皮老紋深、須長須清、須有珍珠疙瘩、藥效高及形狀類似野山參等優點而在國際中藥市場享有較高名譽[1]。
人參皂苷是石柱參的主要藥效成分,屬于三萜皂苷[2-6],具有抗衰老、抗癌、抗炎癥、抗疲勞、治療心肌損傷、降血糖等作用。人參皂苷質量分數和組成是評價人參藥材質量的重要指標,但其在石柱參中質量分數較低,且人工合成尚未實現[7-11]。近年來,生物學專家對人參皂苷的體內生物合成途徑進行了大量研究,其中人參皂苷代謝途徑的分子機制研究進展卓著,但對道地產區石柱參人參皂苷合成的分子機制研究還較少[12]。為探究石柱參人參皂苷質量分數低的原因,本研究分別以5a生石柱參根、莖、葉組織為材料,測定石柱參6種人參皂苷(Rb1、Rb2、Rc、Re 、Rg1、Rg2)的質量分數,及整個生長時期中人參皂苷合成酶基因鯊烯合成酶(SS)、細胞色素P450酶—CYP82D47、CYP716A42的表達量,并對人參皂苷合成酶基因(CYP82D47)進行克隆,測序及生物信息學分析,為進一步分析人參皂苷合成酶基因表達量的動態變化與人參皂苷質量分數的關系,明確石柱參人參皂苷生物合成的分子生理機制,以期為石柱參栽培和人參藥材質量的提高奠定理論基礎。
分別以5 a生人參(對照組)、石柱參的根、莖、葉組織為試驗材料,2種試驗材料均于2017年5-9月采于丹東市寬甸滿族自治縣振江鎮大青村,每個生長時期均采集固定田壟的人參、石柱參各15株,采集后篩選出無病蟲害、個體差異小且生長勢良好的5株樣本進行后續試驗(表1),樣本液氮冷凍于-80 ℃,待用。
1.2.1 人參皂苷質量分數測定 稱取0.5 g 已干燥至恒質量的各樣品粉末,過24目篩,加入250 mL甲醇,于55 ℃水浴恒溫回流提取1 h,過濾提取液,濾渣重復提取3次,用適量體積溶液洗滌殘渣,合并濾液,60 ℃減壓濃縮,定容至5 mL,用0.45 μm濾膜過濾,取濾液為供試品;采用Agilent 1200 Series高效液相色譜儀對人參、石柱參不同部位的6種人參皂苷(Rb1、Rb2、Rc、Re、Rg1、Rg2)的質量分數進行測定[13],色譜柱為德國Agilent Eclipse XDB-C18(4.6 mm×250 mm, 5 μm);乙腈-水溶液為流動相,梯度洗脫,流速1.5 mL/min;A為水,B為乙腈;梯度洗脫程序: 0~16 min,56%B;16~20 min,56%→100% B; 20~40 min,100% B;40 min~47 min,56% B; 47~60 min,56% B;各組分均在70 min內出完,柱溫 35 ℃,檢測波長203 nm;各樣本人參皂苷質量分數的測定重復4次,取平均值。

表1 石柱參樣品采集信息Table 1 Collection information of Shizhu ginseng samples
注:為方便后續試驗記錄,生長時期順序依次以字母A~K記錄。
Note:Each growth period was recorded as the alphabet A-K,just in order to facilitate follow-up test records.
1.2.2 總RNA提取 分別取人參與石柱參根、莖、葉組織于液氮中單獨研磨,按照Takara公司總RNA提取試劑盒操作步驟提取總RNA,并檢測RNA完整性。并采用M-MuLV第1鏈cDNA合成試劑盒,以總RNA為模板,Oligo-dT為引物,逆轉錄合成cDNA[14-18],-20 ℃保存, 備用。
1.2.3 引物設計 根據GenBank上已報道的鯊烯合成酶(SS)、人參細胞色素P450的CYP82D47、CYP716A42基因序列為參考,采用DANMAN、ClustalW、Blast等軟件找出高度保守區設計特異性克隆引物和熒光定量PCR特異性引物(表2)。
1.2.4 人參皂苷合成酶基因表達量測定 以5 a生石柱參11個生長時期(表1,A~K時期)根、莖、葉組織的cDNA為模板,分別對石柱參SS、細胞色素酶P450的CYP82D47、CYP716A42基因進行qPCR擴增,每個反應體系重復3次。qPCR在ABI Biosystems的熒光定量PCR儀7500上完成的,RT-PCR體系用2xSG Fast QPCR Master Mix試劑盒,采用PrimerExpress軟件設計實時熒光定量PCR引物,同時以人參親環蛋白基因CCYP為qPCR內參基因(表2)。不同生長期的表達量以4月26日石柱參根組織表達量為對照[19-20]。qPCR擴增體系為20 μL:SYBR Premix ExTaq10 μL,Forward Primer 0.5 μL,Reverse Primer 0.5 μL,ROX Reference Dye 0.5 μL,cDNA模板 2 μL,ddH2O 6.5 μL;熒光定量qPCR反應程序:94 ℃預變性1 min; 94 ℃變性5 s, 55 ℃退火30 s,72 ℃延伸30 s,40個循環。

表2 引物信息Table 2 The information of primers
1.2.5 重復性試驗 從“1.2.4”中5 a生石柱參不同生長時期根、莖、葉組織的相同部位分別取3份樣品,分別在同等條件下進行擴增,每份樣品平均做3個平行試驗,作為生物學重復,分析其組內及組間差異,以檢測該熒光定量PCR方法的重復性及穩定性。
1.2.6 石柱參CYP82D47基因擴增 分別以石柱參綠果期根、莖、葉為材料,以逆轉錄的cDNA為模板進行PCR反應,克隆CYP82D47基因。反應體系[12]為20 μL:ddH2O 7.8 μL,Mix 10 μL,Forward Primer 0.6 BμL,Reverse Primer 0.6 μL,DNA模板 1 μL; PCR反應程序:94 ℃預變性8 min;94 ℃變性30 s, 55 ℃退火30 s, 72 ℃延伸3 min,35個循環;72 ℃延伸10 min,產物4 ℃保存。
1.2.7 PCR產物回收、連接、轉化與測序 采用SanPrep柱式DNA膠回收試劑盒對PCR產物進行回收、純化,采用Takara公司Pmd20-T試劑盒,將得到的較高純度DNA連接到pMD20-T載體上,后將連接產物轉化到Top10 DH5α大腸桿菌感受態細胞中,以菌液為模板進行PCR擴增[22],最終將陽性克隆的菌液送至上海生工生物工程有限公司進行測序。
1.2.8 系統進化樹構建 將核苷酸序列在Blast在線數據庫進行同源性比對,以此檢測克隆結果是否正確;并在NCBI核苷酸數據庫和氨基酸數據庫中搜索近緣種的細胞色素P450基因序列信息,采用DNAMAN軟件分析序列的開放閱讀框、起始密碼子和終止密碼子,并翻譯出氨基酸序列,最后通過MEGA 5.0構建Neighbor-joining系統進化樹。
采用MxPro與SPSS 19.0對試驗數據進行方差分析和線性回歸分析,以P<0.05表示差異顯著,P<0.01表示差異極顯著。
石柱參3個組織的總RNA電泳圖譜如圖1,28S rRNA和18S rRNA的條帶完整且清晰,二者亮度與寬度比值均接近2∶1,表明提取的總RNA完整性較好;經檢測OD260/OD280為1.9~2.1,表明總RNA質量較高,無污染,可用于后續試驗。

M.2 kb plus DNA ladder;1~3.石柱參根、莖、葉 Root,stem and leaves of Shizhu ginseng;4~6.人參根、莖、葉 Root,stem and leaves of ginseng
圖1 石柱參總RNA的提取
Fig.1 Total RNA extraction from shizhu ginseng
采用高效液相色譜法測定人參、石柱參綠果期不同部位人參皂苷,試驗重復測定4次,取平均值(表3)。可見,人參不同部位根、莖、葉人參皂苷總和的質量分數分別為1.174 4%、7.623 7%、8.933 5%,均明顯高于石柱參根莖葉中人參皂苷的質量分數0.584 4%、1.895 2%、3.307 0%,可能是由于人參皂苷合成途徑中某種人參皂苷合成酶基因的表達量過低造成的;此外,人參、石柱參莖、葉中人參皂苷質量分數相對較高,尤其是人參中Re達到4.013 8%(P<0.05)。分別取配好的人參皂苷標準品儲備液稀釋得到一系列不同質量濃度混合標準操作液,以質量濃度(mg/mL)為橫坐標x,峰面積為縱坐標y,獲得6種人參皂苷線性回歸方程,由表4可見,人參皂苷Rb1、Rb2、Rc、Re 、Rg1、Rg2在標準液質量濃度為0.016~ 3.500 mg/mL呈良好的線性關系。
采用實時熒光定量PCR技術,以人參Cyclophilin作為內參基因,檢測5 a生石柱參11個生長時期(A~K)根、莖、葉組織中人參皂苷合成酶基因(SS、CYP82D47、CYP716A42)的表達量。
SS基因在5 a生石柱參根、莖、葉中不同生長時期的表達量見圖2。可見,石柱參根組織SS基因表達量從展葉期到紅果期無明顯變化,隨后至果熟期降為最低值,而在根生長期表達量有回升,但不高;莖組織SS基因表達量從展葉期到果熟期逐漸降低,尤其在根生長期幾乎不表達,而在枯萎期有少量回升;SS基因在葉組織展葉期表達量較高,從始花期至綠果初期表達量穩步上升,達到峰值,約為根部表達量的3倍,在綠果末期至枯萎期其表達量逐步下降,而此時莖、葉組織的SS基因幾乎不表達,僅在根部有少量表達,說明SS基因的表達具有明顯的組織特異性,在石柱參根、葉中表達活躍,在莖中表達量有限。可見,石柱參SS基因隨著生長周期的結束,其活躍程度降低,趨向于不表達狀態,經統計學分析后,各組間差異極顯著(P<0.01)。

表3 人參、石柱參不同藥用部位中人參皂苷質量分數Table 3 Mass fraction of ginsenosides in ginseng and shizhu ginseng %

表4 6種人參皂苷的線性回歸方程Table 4 Linear regression equation of six ginsenosides
CYP82D47基因在5 a生石柱參根、莖、葉組織中不同生長時期的表達量(圖3)。可見,CYP82D47基因在根組織中的表達量除展葉期較高外,其余各生長期的表達均處于較低水平;在莖組織中,從展葉期至根生長期Ⅰ的CYP82D47基因表達量逐漸升高,并達到峰值,隨后下降,但仍高于同時期根組織的表達量;在葉組織中,從展葉期至綠果初期,其表達量穩步上升,之后均具有較高的表達量,高于同時期根、莖的表達量,可見,CYP82D47基因在石柱參中的表達部位主要是莖、葉組織,經統計學分析后,各組間差異極顯著(P<0.01)。

“**”表示差異極顯著(P<0.01),“*”表示差異顯著(P<0.05);下圖同 “**”the difference was statistically very significant(P<0.01) ,“*”the difference was statistically significant (P<0.05);the same below
圖2 石柱參人參皂苷合成酶SS基因qPCR表達分析
Fig.2 The qPCR relative expression analysis ofSSgene in shizhu ginseng

圖3 石柱參人參皂苷合成酶 CYP82D47基因qPCR表達分析Fig.3 The qPCR relative expression analysis of CYP82D47 gene in shizhu ginseng
CYP716A42基因在5 a生石柱參根、莖、葉組織中不同生長時期的表達量(圖4)。可見,CYP716A42基因在根部的表達與CYP82D47基因完全不同,其表達量先增高后降低,且在綠果末期達到峰值,約為展葉期的3倍(P<0.05),隨后逐漸下降至枯萎期幾乎不表達;在莖組織中的表達不活躍,在結果期和綠果末期其表達量較高,隨后各生長期幾乎不表達;在葉組織中表達量變化較大,展葉期較高,至盛花期降低至不表達,隨后至綠果初期表達量逐步增加并出現峰值,隨后降低,但同期相比,葉組織中的CYP82D47基因在枯萎期表達水平遠高于根、莖組織,經統計學分析后,各組間差異極顯著(P<0.01)。
綜上所述,人參皂苷合成酶基因在石柱參綠果期的葉組織中相對表達量較高;所以,本試驗檢測同時期人參葉片中人參皂苷合成酶基因的表達量為對照組,結果如圖5,人參、石柱參的人參皂苷合成酶基因SS、CYP716A42在葉片組織綠果期的表達量差異極顯著(P<0.01),基因CYP82D47在葉片組織綠果期的表達量差異極顯著(P<0.01)。
熒光定量PCR對3個平行試驗檢測的結果計算標準差及變異系數表明,同一份樣品的組內重復變異系數為0.12%~1.03%,同一植株相同部位3份樣品間的組間重復變異系數為 0.19%~ 1.92%(表5),表明該方法具有較好的組內和組間重復性。因此,說明本試驗建立的實時熒光定量PCR方法不僅靈敏度和準確性高,而且穩定性較好。

圖4 石柱參人參皂苷合成酶 CYP716A42基因qPCR表達分析Fig.4 The qPCR relative expression analysis of CYP716A42 gene in shizhu ginseng
由于綠果期CYP82D47基因的表達量較高,所以分別以石柱參綠果期根、莖、葉為材料,以逆轉錄cDNA為模板,采用細胞色素CYP82D47基因特異性引物擴增該基因,片段如圖所示(圖6)。

圖5 人參、石柱參人參皂苷合成酶基因表達量的比較Fig.5 Comparison of ginsenoside synthase gene expression in ginseng and shizhu ginseng

M.Marker;1~3.依次為以石柱參根、莖、葉組織cDNA為模板克隆的CYP82D47基因片段 Amplified fragment ofCYP82D47gene from Shizhu ginseng’s root,stem and leaves cDNA as the template
圖6CYP82D47基因克隆圖譜
Fig.6 Cloning map ofCYP82D47gene

表5 重復性試驗結果Table 5 Results of repeated experiments of real-time quantitative PCR assay
將純化得到的目的片段連接到pMD20-T載體上并轉化到大腸桿菌Top10 DH5α的感受態細胞中,采用菌落PCR方法檢測陽性克隆菌液,電泳得到約1 524 bp目的基因條帶,隨后對陽性克隆的菌液進行測序,采用DNAMAN軟件對序列進行拼接,獲得細胞色素P450CYP82D47基因cDNA全長序列,片段大小為1 524 bp,起始密碼子位于114 bp,終止密碼子位于1 205 bp,開放閱讀框總長1 091 bp,共編碼氨基酸324個,序列如圖7。
根據石柱參細胞色素P450CYP82D47基因的氨基酸序列,在NCBI數據庫中檢索與其氨基酸序列相似度較高的序列建立系統進化樹(圖8)。可見,14個遺傳關系較近的品種聚為4個類群,第一類群2個品種,即楤木和刺楸;第二類群5個品種,即通草、羽葉三七、姜狀三七、三七和扣子七;第三類群4個品種,即西洋參、竹節參、人參和石柱參;第四類群3個品種,即無梗五加、細柱五加和刺五加。其中石柱參CYP82D47基因的氨基酸序列與人參、竹節參和西洋參親緣關系較近,說明該結果與物種進化程度相符合。

圖7 CYP82D47基因cDNA序列Fig.7 cDNA sequence of CYP82D47 gene

圖8 石柱參 CYP82D47基因系統發育樹Fig.8 Phylogenetic tree of CYP82D47 gene in shizhu ginseng
人參皂苷質量分數是評價石柱參品質的關鍵,其質量分數受環境和人參皂苷合成酶的協同影響,人參皂苷合成酶基因的表達影響人參皂苷合成途徑中各物質流向。曲正義等[22]對不同生長時期人參根、莖、葉中的人參皂苷質量分數進行測定分析,發現人參中總皂苷質量分數為葉>根>莖,人參葉片中質量分數較高的單體皂苷是Re、Rd、Rb1和Rg1,人參根中質量分數較高的單體皂苷是Rg1、Re和Rb1,人參莖中質量分數較高的單體皂苷是Rg1和Re,牛云云等[23]發現三七人參皂苷合成酶基因在莖、葉中的表達量明顯高于根中,均與本文研究結果一致;此外,人參皂苷合成酶基因的表達在根、莖、葉組織中具有較高的組織差異性,比如SS基因在石柱參根、葉中表達活躍,而莖中表達量則有限,且石柱參SS基因隨著生長周期的結束,其活躍程度降低,趨向于不表達狀態;CYP82D47基因主要表達部位是莖、葉組織;而CYP716A42基因在根、莖、葉中表達較活躍,尤其是根、葉中;人參皂苷合成酶基因的表達量在綠果期葉片中質量分數較高,但由于基因表達量不穩定且帶有波動性,所以以人參為對照組,檢測二者綠果期葉片中各人參皂苷合成酶基因的表達量,結果發現相同條件下,石柱參人參皂苷質量分數較人參差了近3倍,Han等[24]研究發現,人參CYP82D47在茉莉酸甲酯調控下表達量上調,而CYP82D47是一個尚未明確功能的P450羥化酶,所以造成人參皂苷質量分數較低的原因之一可能是人參皂苷合成酶基因CYP82D47的表達量相對較低,影響了人參皂苷合成過程中量的積累,該基因主要催化達瑪烷型原二醇到原人參二醇的C-12位羥基化,推進了人參皂苷的生物合成。此外,克隆CYP82D47基因cDNA全長序列,進行生物信息學分析,由CYP82D47基因氨基酸序列及檢索相似度較高的序列建立石柱參系統進化樹,可見,石柱參與人參、竹節參、西洋參等的親緣關系較近,結果與物種進化程度相符合。CYP716A42和CYP82D47同為細胞色素P450超基因家族成員,深入研究CYP450s,并進一步驗證其功能和作用位點,闡述人參皂苷合成酶基因的表達如何調控人參皂苷生物合成的生理生態機制,對今后在微生物中重建人參皂苷合成途徑、通過基因工程手段培育高質量分數人參皂苷的石柱參具有重要的實踐意義。