基于DNA條形碼技術的維吾爾藥材桑葚基原鑒定研究

樊叢照 徐建國 李亞偉 李曉瑾

[摘要]該研究通過分析多途徑來源的新疆桑屬植物及藥材樣本的ITS2，psbAtrnH序列，為藥材分子鑒定提供依據。以51個新疆桑屬植物及藥材為樣本，對其ITS2，psbAtrnH序列進行PCR擴增和測序，用MEGA 60計算其種內、種間Kimura 2parameter（K2P）距離，分析變異位點，并構建NJ鑒別樹。ITS2序列分析結果顯示，桑Morus alba、韃靼桑M alba var tatarica、黑桑M nigra種內無變異；桑與韃靼桑種間無變異；桑與黑桑種間存在13個變異位點，種間平均KP2遺傳距離為004；桑與藥材樣本間無信息變異位點，NJ鑒別樹可將桑及韃靼桑與黑桑區分。psbAtrnH序列分析結果顯示，桑與黑桑種內各有1個變異位點，3種植物種間存在插入/缺失變異，可相互區分；種間變異與藥材樣本內變異一致。因此，ITS2序列可將來源于桑、韃靼桑的藥材樣本與黑桑區分，psbAtrnH序列可將來源于三者的藥材樣本區分，為維吾爾藥材真偽鑒別及市場監管提供依據。

[關鍵詞]維吾爾藥； 桑； DNA條形碼； ITS2； psbAtrnH

Identification of origin plant of Uygur medicine mulberry

based on DNA barcode

FAN Congzhao1， XU Jianguo1， LI Yawei2， LI Xiaojin1*

（1Xinjiang Institute of Chinese Materia Medica and Ethnical Materia State Administration of Traditional Chinese

Medicine， Key Laboratory of Traditional Chinese Medicine and Ethnic Medicine Resources， Urumqi 830002， China；

2 Testing Center of Xinjiang Entry Exit Inspection and Quarantine Portal， Urumqi 830063， China）

[Abstract]To provide molecular evidence for medical material identification， we analyzed the nucleotide sequence of ITS2， psbAtrnH gene in Morus genus plants and commercial products which were obtained from different places in Xinjiang The sequence of ITS2 and psbAtrnH in fiftyone samples were amplified and sequenced， MEGA 60 was used to analyze the intra and interspecific K2P distances， neighborjoining （NJ） tree was used to constructing clustering tree ITS2 sequence analyzed results showed that there is no intraspecific variation among Morus alba， M alba var tatarica and M nigra， but 13 variations sites were exist between M alba and M nigra and their interspecific K2P distances was 004， which indicated that there had significant variation in them We didn′t find informative variation sites between Morus genus plants and commercial products， and we also found that M nigra can be distinguished from other two species by NJ Tree PsbAtrnH analysis results showed there was only one variation site between M alba and M nigra， but insertion or deletion variation were remarkable evidence among M alba， M alba var tatarica and M Nigra Interspecific variation was accordance with intraspecific variation of commercial products So ITS2 and psbAtrnH gene were important marker for M alba， M alba var tatarica and M nigra identification This study provided important evidence for Uygur medicine identification and market supervision

[Key words]Uygur medicine； mulberry； DNA barcode； ITS2； psbAtrnHendprint

與漢民族“前不栽桑，后不栽柳”的習俗不同，維吾爾族民居具有庭院中栽種桑樹的傳統[1]，新疆桑樹種質資源豐富，不同產地性狀存在差異，種和品種之間的劃分界限模糊不清。新疆民間通常以桑椹的顏色或桑花的性別為依據，把桑樹分為黑桑、白桑、粉桑及公桑、藥桑5個類型[2]，《新疆植物志》將其分為白桑Morus alba、黑桑M nigra及變種韃靼桑M alba var tatarica[3]。桑屬植物具有重要的藥用價值，《中國藥典》及《新疆維吾爾自治區中藥維吾爾藥飲片炮制規范》均收載桑M alba為基原植物[45]，其藥性為一級濕熱，主治干寒性或黑膽質性疾病，黑桑雖未載入質量標準，亦作為維吾爾常用藥材使用，但其藥性為一級干二級寒，主治熱性或血液質性和膽液質性疾病，藥性與功能均與白桑存在一定的差異[67]，因此二者不能混用。

正確的藥材基原是保證藥材質量及療效的前提，維吾爾藥材質量控制水平整體較低，市售藥材基原不明極為常見[8]。維吾爾藥材桑葚多來源于民間貿易，不同植物來源的桑葚藥材形態特征相近，臨床應用時不易區分，具有潛在的隱患。利用傳統的基原鑒定和性狀鑒定難以準確快速地進行鑒別。DNA條形碼技術具有通用性強、鑒定結果可靠、重復性良好等優點，被廣泛用于中藥民族藥材鑒定中[9]。《中國藥典》2015 年版已將DNA條形碼ITS2及psbAtrnH序列作為候選序列收載[10]。本研究擬采用ITS2及psbAtrnH序列，通過分析多途徑來源的桑屬植物及藥材樣品，建立維吾爾醫常用藥材桑及其近緣種的鑒定方法，為維藥桑葚的質量控制及標準制定提供科學的依據，為正確合理利用維吾爾藥材資源提供理論基礎。

1材料

51份實驗樣品包括葉片30份、藥材樣品21份，基原植物由中國醫學科學院藥用植物研究所林余霖研究員及新疆維吾爾自治區中藥民族藥研究所王果平副研究員鑒定，憑證標本保存于新疆中藥民族藥研究所標本館（新疆XTNM），實驗所獲得基原植物的單倍型序列已提交至GenBank，見表1。

PCR擴增所用引物ITS2序列由上海生工生物工程合成，ITS2正向序列：5′ATGCGATACTTGGTGTGAAT3′，反向：5′GACGCTTCTCCAGACTACAAT3′；DNA提取試劑盒、DNA聚合酶及dNTP等均購自天根生物；psbAtrnH正向序列：5′GTTATGCATGAACGTAATGCTC3′，反向：5′CGCGCATGGTGGATTCACAATCC3′[11]。DNA提取試劑盒、DNA聚合酶及dNTP等均購自天根生物（Tiangen Biotech Co）。

DNA提取研磨儀GT100（GRINDER，中國），Anker TGL16C離心機（上海安亭科學儀器）； PCR擴增儀為070851（An Analytik Jena company，德國）。

2方法

21DNA提取及檢測稱取干燥葉片或果實30 mg，DNA提取研磨儀1 000 r·min-1研磨2 min，植物DNA提取試劑盒提取總DNA，1%瓊脂糖凝膠電泳檢測。

22PCR擴增及測序PCR反應體系、反應條件、序列拼接、序列注釋及分析對比參照Chen等研究結果[12]；使用MEGA 60分析序列，計算種內及種間遺傳距離，用NJ法構建聚類樹[13]。

3結果與分析

31序列信息ITS2序列分析結果顯示，桑為1個單倍型，長度236 bp，GC量631%；韃靼桑序列信息與桑一致，黑桑單倍型數為1，長度為236 bp，GC量為610%；藥材樣本有2個單倍型，序列長度為236 bp，GC量為610%，631%。PsbAtrnH序列分析結果表明，桑有2個單倍型，長度為452 bp，GC量為239%，241%；韃靼桑為1個單倍型，序列長度為470 bp，GC量234%；黑桑單倍型數為2，序列長度458 bp，GC量236%，234%；藥材樣本有3個單倍型，序列長度分別為452，458，470 bp，GC量分別為234%，236%，239%，見表2。

32種內及種間K2P遺傳距離在ITS2序列中，桑、韃靼桑及黑桑種內K2P遺傳距離均為0，藥材樣本內K2P遺傳距離為004；桑與韃靼桑種間無變異，兩者與黑桑種間K2P遺傳距離為004；桑、韃靼桑、黑桑與藥材樣本之間的K2P遺傳距離均為0～004，藥材與植物種內種間K2P遺傳距離差異不顯著。根據psbAtrnH序列分析結果，韃靼桑種內K2P遺傳距離為0，桑、黑桑種內K2P遺傳距離均為0002，藥材樣本內K2P遺傳距離為0～0002；桑、韃靼桑、黑桑種間K2P遺傳距離變化范圍為0001～0002，三者與藥材之間K2P遺傳距離變化范圍0～0004，種內種間K2P遺傳距離無顯著差異，見表3。

33種內及種間變異位點分析新疆桑屬植物及其藥材ITS2序列比較后，共有13個變異位點，見表4。包括9個轉換/顛換突變和4個插入/缺失變異，13個變異位點的堿基差異均可以將桑與黑桑區分開；藥材樣本中變異位點包括了桑與黑桑2種類型。

PsbAtrnH序列總變異位點數為27，見表5。包括3個轉換/顛換突變和24個插入/缺失變異，桑種內474位點存在A/C變異，黑桑種內73 bp位點存在G/T變異。66～71 bp位點的AATATT插入/缺失變異可將黑桑與其桑、韃靼桑區分，425～442 bp位點的插入缺失變異能將韃靼桑與桑、黑桑區分。

34桑屬植物NJ樹鑒別采用NJ樹法鑒定，基于ITS2序列，桑、韃靼桑可以與黑桑區分，但桑與韃靼桑無法區分，與編號YC125006為相同序列類型的藥材均與桑、韃靼桑聚在一起，與編號YC125001為相同序列類型的藥材均與黑桑聚在一起；基于psbAtrnH序列，黑桑單獨聚為一支，并與編號為YC125001的藥材聚在一起，見圖1，2。endprint

4討論

41DNA條形碼的選擇桑樹是一種集食用、藥用、生態保護等多種經濟價值為一體的經濟樹種[14]，在中國已經具有5 000多年的栽培歷史[15]，長期的栽培和廣泛的地理分布使其種間雜交嚴重，分類至今存在爭議[16]，因此傳統的形態學在區分不同種植資源存在一定的難度，桑M alba的干燥果實作為藥材，僅憑外觀性更是狀難以與桑屬其他種類的果實區分。諸多研究者通過遺傳多樣性、核糖體序列及葉綠體序列分子標記方法對不同來源的桑樹資源進行了鑒定[17]，雖然研究表明ITS序列區分不同種類的桑樹具有較高的效率，但ITS2序列在藥材樣品擴增中效率更高[18]，此外psbAtrnH序列是進化速率最快的葉綠體間隔區之一，作為植物DNA條形碼的候選片段，psbAtrnH序列具有引物通用性較好、擴增成功率較高，為單親遺傳，不存在序列雜合的影響[1922]等特點，因此本研究將ITS2及psbAtrnH序列作為候選DNA條形碼對不同來源的桑屬植物樣本和藥材進行鑒定。

42DNA條形碼可成功鑒定不同來源的藥材桑葚本研究選擇DNA條形碼ITS2，psbAtrnH等2條序列對不同來源的新疆桑屬植物樣本和藥材桑葚進行鑒定，結果表明，無論從葉片還是從藥材桑葚中均可成功提取基因組DNA，擴增效率均為100%。種間變異及NJ樹鑒定結果表明，ITS2序列可以鑒定藥材桑與黑桑，但不能將桑與韃靼桑區分開，作為變種其是否可以與桑同等入藥還未見此類報道；psbAtrnH序列中的插入/缺失變異能成功鑒定來源于桑、韃靼桑及黑桑的藥材，但插入缺失變異無法在NJ鑒別樹種表現出來。因此，ITS2與psbAtrnH片段結合使用，可成功鑒定來源于桑、韃靼桑和黑桑的藥材。

43DNA條形碼可為維吾爾藥材桑葚的正確應用提供依據相比中藥材，維吾爾藥材來源更為復雜，在長期的用藥過程中由于翻譯錯誤、鑒定錯誤、代用混用等原因，藥材質量標準整體偏低，藥材監管難度大[8]，諸多維吾爾藥材的來源無法追溯，考證難度較大。本研究所收集的藥材桑葚主要來源于維吾爾醫藥市場，鑒定結果表明，市場中的藥材包括桑、韃靼桑和黑桑，黑桑與桑屬于不同的藥材種類，藥效藥性也存在差異，長期誤用混用不僅影響了維吾爾藥材標準的權威性，也造成了潛在的用藥安全隱患。DNA條形碼技術可以追溯藥材的基原，快速準確的鑒定維吾爾藥材，在今后藥材來源、標準制定、藥材流通和市場管理等實際應用中具有重要的價值。

[參考文獻]

[1]盧紅， 左少純， 吳麗莉， 等 新疆維吾爾庭院桑文化內涵與文化經濟模式思考[J]. 蠶學通訊， 2014， 34（3）： 8.

[2]張思春， 朱美武 新疆的桑樹資源[J]. 中國蠶業， 1995， 12（1）： 44.

[3]新疆植物志編輯委員會 新疆植物志第1卷[M] 烏魯木齊： 新疆科技衛生出版社， 1992： 218.

[4]中國藥典一部[S] 2015： 300.

[5]新疆維吾爾自治區食品藥品監督管理局 新疆維吾爾自治區中藥維吾爾藥飲片炮制規范[S] 2010： 214.

[6]國家中醫藥管理局中華本草編委會 中華本草·維吾爾藥卷[M] 上海： 上海科學技術出版社， 2005： 263.

[7]劉勇民 維吾爾藥志[M] 烏魯木齊： 新疆科技衛生出版社， 1999： 701.

[8]蘇來曼·哈力克， 孫磊， 沈曉麗，等 維吾爾藥質量標準現狀分析及發展建議[J]. 中國藥事， 2015（12）： 1256.

[9]陳士林， 姚輝， 韓建萍， 等 中藥材DNA條形碼分子鑒定指導原則[J]. 中國中藥雜志， 2013， 38（2）： 141.

[10]中國藥典四部[S] 2015：383.

[11]陳士林 中藥DNA條形碼分子鑒定[M] 北京： 人民衛生出版社， 2012： 14.

[12]陳士林， 龐曉慧， 姚輝， 等 中藥DNA條形碼鑒定體系及研究方向[J]. 世界科學技術——中醫藥現代化， 2011（5）： 747.

[13]Koetschan C， Hackl T， Müller T， et al. ITS2 Database Ⅳ： interactive taxon sampling for internal transcribed spacer 2 based phylogenies[J]. Mol Phylogenet Evol， 2012， 63（3）： 585.

[14]Priya S Medicinal values of mulberryan overview[J]. J Pharmacy Res， 2012， 5： 3588.

[15]He N， Zhang C， Qi X， et al. Draft genome sequence of the mulberry tree Morus notabilis[J]. Nat Commun， 2013， doi： 101038/ncomms3445.

[16]Burgess K S， Morgan M， Deverno L， et al. Asymmetrical introgression between two Morus species （M alba， M rubra） that differ in abundance[J]. Mol Ecol， 2005， 14： 3471.

[17]Zeng Q， Chen H， Zhang C， et al. Definition of eight mulberry species in the genus Morus by internal transcribed spacer based phylogeny[J]. PLoS ONE， 2015， 10（8）： e0135411.

[18]Fan C Z， Li X J， Zhu J， et al. Endangered Uyghur medicinal plant ferula identification through the second internal transcribed spacer[J]. EvidBased Compl Alt， 2015， 2015（1）： 1.

[19]Li X W， Hu Z G， Lin X H， et al. Highthroughput pyrosequencing of the complete chloroplast genome of Magnolia officinalis and its application in species identification[J]. Acta Pharm Sin， 2012， 47： 124.

[20]Shaw J， Lickey E B， Beck J T， et al The tortoise and the hare Ⅱ：Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis[J]. Am J Bot， 2005， 92（1）： 142.

[21]Kress W J， Efickson D L A twolocus global DNA barcode for land plants： the coding rbcL gene complements the noncoding trnHpsbA spacer region[J]. PLoS ONE， 2007， 2（6）： e508.

[22]Fazekas A J， Burgess K S Kesanakurti P R， et al Multiple muhilocus DNA barcodes from the plastid genome discriminate plant species equally well[J]. PLoS ONE， 2008， 3（7）： e2802

[責任編輯呂冬梅]endprint