基于ITS序列分析技術鑒定川牛膝與常見偽品麻牛膝

陳詩晴，蒲沁琳，陳小軍，宋瑱

·品種品質·

基于ITS序列分析技術鑒定川牛膝與常見偽品麻牛膝

陳詩晴1，蒲沁琳1，陳小軍2，宋瑱1

目的：比較川牛膝與其常見偽品麻牛膝之間的ITS序列差異及規律。為川牛膝與麻牛膝的DNA條形碼鑒別提供適合的分子標記。方法：收集川牛膝及麻牛膝成品藥材并提取純化其基因組DNA，經PCR擴增得到ITS序列（包括ITS1、5.8S nrDNA、ITS2）并進行T-A克隆后測序，分析兩者序列差異。結果：PCR擴增獲得兩者ITS序列，經多序列對比分析得出川牛膝與偽品麻牛膝的ITS序列存在明顯差異。結論：ITS序列分析可以用作鑒定川牛膝與麻牛膝藥材。

川牛膝；麻牛膝；PCR；ITS序列；DNA測序

中藥材川牛膝是莧科植物川牛膝（Cyathμla officinalis Kuan）的干燥根，為《中華人民共和國藥典》1990年版收載，是四川省著名的道地藥材，具有逐瘀通經、通利關節、利尿通淋等功能[1,2]。近年來，現代藥理研究證明川牛膝多糖成分具有抗腫瘤、抗生育等新的作用[3]，使其再次受到關注。但是與川牛膝同科植物的麻牛膝入藥部位外觀與川牛膝相似，用傳統鑒別方法（如外觀、形態、質地等）很難準確地把兩種藥材區別開來；中醫臨床經驗認為麻牛膝的性味功效與川牛膝不同，不宜作川牛膝藥用。因此，如何簡便、快速和準確地鑒定川牛膝的真偽是我們廣泛應用道地藥材川牛膝面臨的最大難題。

ITS序列是真核生物細胞基因組中的核糖體rRNA基因的轉錄間隔序列。18S nrDNA序列、ITS1間隔序列、5.8S nrDNA序列、ITS2間隔序列和26S nrDNA依次緊密連接在一起，形成一個完整的轉錄單位，該轉錄單位是典型的高度重復序列，在植物基因組中拷貝數可達500～40000。在技術上易于擴增，其中的ITS1序列和ITS2序列的可變性與物種遺傳關系密切，成為了眾多研究者對于藥用植物的鑒別，揭示藥用植物的遺傳差異和鑒定中藥材品種的道地性的新方法和手段[4～7]。

圖1 ITS完整序列示意圖

1 材料和方法

1.1 材料

1.1.1 藥材樣品 分別從四川雅安天全縣、西昌、寶興等地采集川牛膝藥材10份，西昌會理采集及成都藥材市場購買麻牛膝藥材5份。道地藥材為雅安川牛膝GAP基地產品。采集后置于干燥環境保存。

1.1.2 主要試劑 CTAB、Tris、蛋白酶K購自sigma公司；HOTSTART taq DNA聚合酶，pEASY-T5 zero質粒載體購自北京全式金公司；凝膠回收試劑盒購自北京博邁德公司，磁珠法細菌基因組DNA提取試劑盒及擴增引物購自上海生工公司。

表1 擴增引物及其序列

1.2 方法

1.2.1 基因組DNA分離純化 采用本實驗室經過預試驗改良后的磁珠法提取，以75%酒精消毒川牛膝及麻牛膝藥材，剔去表面后液氮研磨成粉末，10%CTAB/NaCl裂解液裂解樣品粉末后加入吸附性磁珠。以磁力架吸附包含DNA的磁珠，加入70%的乙醇洗滌磁珠，揮干乙醇后加入50μL洗脫液洗脫獲得基因組DNA。將獲得的DNA用0.7%瓊脂糖凝膠電泳鑒定完整性，紫外分光光度儀檢測OD260值計算濃度。

1.2.2 PCR擴增ITS序列 參照NCBI上牛膝屬ITS完整序列（包括ITS1、5.8S nrDNA、ITS2，GenBank: DQ497186.1）設計引物ITS1及ITS4。以上一步所得的DNA為模板構建川牛膝ITS序列PCR反應體系（50μL）：10×Buffer，200mM dNTP，0.25uM引物，1U HOTSTART taq DNA聚合酶，100ng模板，補超純水至50μL。擴增反應參數為：94℃預變性5分鐘；94℃變性30秒；52℃退火30秒；72℃延伸1分鐘；循環30次后72℃延伸7分鐘。1.5%瓊脂糖凝膠電泳檢測PCR產物。

1.2.3 PCR產物的T-A克隆及序列分析 將PCR產物經紫外分光定量后按摩爾比7: 1加入pEASY-T5 zero質粒載體，25℃連接10分鐘，42℃熱激30秒轉化JM109感受態細胞，在含氨芐青霉素（100mg/L）的SOC平板上培養24小時后進行菌落PCR篩選，獲得陽性克隆。將篩選后的陽性克隆送上海生工公司測序鑒定。

2 結果

2.1 川牛膝及麻牛膝藥材（根莖）的基因組DNA的鑒定

取“1.2.1”中獲得牛膝藥材DNA溶液5μL上樣于0.7%瓊脂糖凝膠電泳結果見圖2。圖中可以分析出DNA分子完整，無明顯降解，且濃度較高，可以滿足下游PCR實驗。紫外分光檢測A260/A280均為1.7-2.0范圍內，無明顯蛋白污染。

圖2 川牛膝及麻牛膝基因組DNA的電泳鑒定結果

2.2 PCR擴增ITS序列的電泳結果

以1.2.2方法擴增16份樣品ITS序列經電泳檢測均獲得明顯的特異性產物，大小670bp。見圖3。

圖3 川牛膝及麻牛膝ITS序列PCR產物的電泳結果

2.3 川牛膝及麻牛膝ITS序列測序結果分析

經測序后，參照NCBI中莧科植物ITS序列經BlAST比對，以及DNAMAN軟件多序列分析，確定11株川牛膝ITS序列長度約為591bp，G+C%約為53.3-53.7%。其中ITS1序列約232bp、5.8S nrDNA序列約162bp、ITS2序列約197bp。5株麻牛膝ITS序列長度約為587bp，G+C%約為51.6-52.1%。其中ITS1序列約230bp、5.8S nrDNA序列約162bp、ITS2序列約195bp。所有樣品5.8S nrDNA序列都高度保守，ITS1及ITS2序列在同種間只有少量差異，在川牛膝與麻牛膝間有明顯的差異30余處，見圖4。

圖4 川牛膝及麻牛膝ITS序列差異對比分析結果

3 結論與討論

本次研究對于樣品的采樣部位是成品藥材即根莖部位，而由于植物根莖的DNA分離純化難度較高，所以現有眾多研究者采樣部位主要是采用了新鮮葉片提取DNA，方法以CTAB法為主[8,9]。但此方法不適合對藥材成品的DNA提取。本研究采用了優化后的磁珠法。磁珠法提取生物樣本DNA在方法上有快速、樣品純度及濃度高的優點，但對于植物根莖的細胞裂解不及CTAB法效率高。本次研究綜合利用磁珠法和CTAB法的優點，先利用CTAB裂解液處理樣品裂解細胞，再使用磁珠法收集純化DNA，獲得了良好的效果。

對11株川牛膝及5株麻牛膝ITS序列比對分析發現，在同種的川牛膝ITS序列中5.8S nrDNA序列高度保守，未發現有突變位點，ITS1及ITS2序列在隨產地不同只有少數幾個堿基有差異，道地藥材和其余10株不同產地川牛膝沒有明顯穩定的差異。所以，ITS序列分析不能用作鑒別川牛膝道地藥材，需要另尋其他分子標記。

在川牛膝和麻牛膝不同種牛膝間發現，其ITS1及ITS2序列在42位、50位、88位、90位等30多處位點有明顯的堿基替換、缺失、插入等差異。并且這些差異具有較好的穩定性。經過重復雙向測序驗證，這些結果準確可靠。

通過本研究發現ITS序列測序后經過序列比對能較好的反映出川牛膝和麻牛膝的分子差異。因此可利用這種技術鑒定川牛膝和麻牛膝藥材。

[1] 張敬杰,蘆衛紅,史旭霞.牛膝的草本學研究[J].現代中西結合雜志,1994,7(11):696.

[2] 胡世林.中國道地藥材論叢[M].北京:北京中醫古籍出版社, 1997:62.

[3] 陳紅,劉友平.川牛膝多糖抗腫瘤作用初探[J].成都中醫藥大學學報,2001,24(1):16.

[4] Gottschling M,plotner J.plotner.Secondary structure models ofthe nuclear intemal -transcribed spacer regions and 5.8S rRNA in CaleiodinelIoideae(Peridiniaeeae) and other dinofagellates[J].Nucleic Acids Res,2004,32(1):307．

[5] Baldwin BG,Sandersonm J,Porter JM．TheITSIegion of nuclear ribosomal DNA:A valuable soμl'ce of evidence on angiosperm phylogeny[J].Ann Missouri Bot Gard,1995,82:247.

[6] Hyosig Won,Susanne S,Renner.The internal transcribed spacer of nuclear ribosomal DNA in the gymnosperm Gnetum[J]. Molecμlar Phylogenetics and Evolution,2005,36:581.

[7] 趙歡,吳衛.核糖體DNA ITS序列分析在藥用植物研究中的應用[J].時珍國醫國藥,2009,20(4):959.

[8] 丁平,方琴.巴戟天與常見混偽品的rDNA-ITS序列分析及其分子鑒定[J].中草藥,2005,36(6):908.

[9] 趙永亮,傅體華. 川牛膝(Cyathula oficinalis Kuan)RAPD和AFLP標記的多態性聚類分析[J].安徽農業科學, 2008,36(16): 6682.

（責任編輯：蔣淼）

The identifcation of Chuanniuxi and adulterants Maniuxi based on ITS sequence analysis

CHEN Shi-qing1,PU Qinglin1,CHEN Xiao-jun2,SONG Zhen1//(1. Chengdu University of Traditional Chinese Medicine, Chengdu 611137;2. Sichuan Lu Tian Hua Co., Ltd. Luzhou 646300)

Objective:To fnd the difference of ITS sequences and regularity between Chuanniuxi and its common adulterants Maniuxi, and provide suitable molecular markers for the DNA barcoding identification of Chuanniuxi and Maniuxi. Method:The samples of Chuanniuxi and Maniuxi were collected, and genomic DNA were extracted to get the ITS sequences by PCR amplification (including ITS1, 5.8S nrDNA and ITS2). Then ITS sequence and sequence divergence were analyzed after T-A cloning.Result:The ITS sequence were got after PCR amplifcation and signifcant differences between Chuanniuxi and Maniuxi was found by multiple sequence alignment analysis. Conclusion:ITS sequence analysis can be used to identify the raw material Chuanniuxiu and Maniuxi.

Chuanniuxi; Maniuxi; PCR; ITS sequences; DNA sequence analysis

R 282.5

A

1674-926X(2014)04-002-03

成都中醫藥大學科技發展基金(編號ZRYB2010)；成都中醫藥大學大學生科研實踐創新課題（編號2012-90）

1.成都中醫藥大學醫學技術學院，四川 成都611137；2.四川省瀘天化股份有限公司中化室，四川 瀘州 646300

陳詩晴（1993-），女，本科，成都中醫藥大學醫學技術學院生物技術專業Tel: 18215525190 Email:502590353@qq.com

宋瑱（1979-），男，醫學碩士，主管檢驗技師，主要從事醫學分子生物學研究Email:red_cell@qq.com

2013-10-24