DNA條形碼在主要滇產石斛鑒別中的應用

張愛麗++高亞芳++黎太米++譚文紅 錢子剛

【摘要】目的：選擇云南境內應用面較廣、市場流通量較大的14個石斛品種，建立其 DNA條形碼鑒定方法。方法：結合Genbank數據庫序列，選擇ITS、psbA-trnH、matK、rbcL等4條常用序列對石斛DNA樣品進行擴增和測序，利用MEGA50計算種間種內 K2P遺傳距離并構建NJ樹，同時采用Taxon DNA軟件計算Barcoding gap（條形碼種內種間遺傳差異）。結果：ITS序列的K2P種內遺傳距離小于種間遺傳距離，NJ樹可以很好地將不同品種的石斛分開，且有較明顯的Barcoding gap。結論：ITS序列可以作為鑒別這14種滇產主要石斛的優選序列。

【關鍵詞】滇產石斛；DNA條形碼；分子鑒別

【中圖分類號】R197【文獻標志碼】 A【文章編號】1007-8517（2017）02-0011-06

Application of DNA Barcoding in species identification of Dendrobium widely used in Yunnan

ZHANG AiliGAO YafangLI TaimiTAN Wenhong*QING Zigang*

Yunnan University of Chinese Traditional Medicine，Kunming 650500，China

Abstract：Objective To identify 14 species of Dendrobium which widely used in Yunnan province by DNA barcoding. Methods Three chloroplast sequences including psbA-trnH， matK， rbcL and an ribosomal DNA ITS sequence of the 14 species were comparatively analyzed. K2P distances were calculated and NJ Tree was performed applying MEGA 50Meanwhile， Barcoding gap were calculated by Taxon DNA. Results The interspecific K2P distance of ITS were higher than the intraspecific one. And NJ tree could also divide species into different clades. The barcoding-gap were obvious between intra- and inter-species. Conclusions ITS sequence can be the optimal barcode to identify the 14 species of Dendrobium.

Keywords：Genus of Dendrobium in Yunnan； DNA Barcoding； Molecular Identification

石斛為我國著名中藥材，具有生津益胃、滋陰清熱、潤肺止咳等功效[1]。是歷代本草和歷版《中國藥典》收載的品種。《中國藥典》（2015年版）收載的石斛藥材為蘭科植物金釵石斛（Dendrobium nobile）、鼓槌石斛（Dendrobium chrysotoxum）、流蘇石斛（Dendrobium fimbriatum）的栽培品及同屬植物的近似種[2]。

石斛作為云南特有藥材之一，產量占全國的60%。在云南境內，除藥典品種鐵皮石斛以外，市場流通的石斛品種多達20幾種[3]。但這些品種沒有專屬的質量標準，使得這類石斛藥材質量得不到控制。因此，控制其質量、保證臨床用藥的安全和療效，必先對石斛藥材的來源進行嚴格的鑒定。但是，石斛屬 （Dendrobium SW）作為蘭科 （Orchidaceae） 最大的屬之一，其種類繁多，且許多種類形態非常相似，從外觀形態上對石斛屬植物進行鑒別有一定的困難[4]。

DNA條形碼鑒定是現階段一個精準、穩定并且易于操作的鑒定方法。本實驗收集了云南省境內應用面較廣、市場流通量較大的14個石斛品種，選擇ITS、psbA-trnH、matK、rbcL等4條常用條形碼片段對其進行了序列的分析比較，以期為鑒別這些主要滇產石斛品種提供分子依據，使得其基源鑒定更為準確可靠，為今后石斛的質量控制研究奠定基礎。

1儀器和材料

11樣品來源實驗通過走訪調查，文獻查閱，收集了云南省境內應用面廣、市場流通大的14個石斛品種，每個品種保證至少3個或以上重復樣品，實驗中不能獲得的序列信息，均從Genbank下載。具體信息見表1。

12儀器SW-CJ-1D潔凈工作臺（江蘇蘇潔凈化設備廠）；高速冷凍離心機（HITACHI）；電泳儀（Bio-Rad）；微型電泳槽（Bio-Rad）；凝膠成像系統（Bio- Rad）；PCR反應擴增儀（Bio-Rad）；移液器（范圍100～1000μL，20～200μL，05～10μL（Eppordorf）。

[JP3]13試劑4S Red Plus核酸染色劑[BF]（上海生工，生產批號：A606695）；瓊脂糖（Biowest，生產批號：111860）；Mix Tag酶（Takara公司，生產批號：A3201A）。[BFQ][JP]

2方法

21DNA提取采集的新鮮石斛葉片用變色硅膠干燥， 采用植物基因組DNA快速提取試劑盒 （Bioteke， 北京），提取石斛樣品的總DNA。

22PCR擴增PCR 體系總體積為 25 μL， 其中包含： 2×Taq PCR Mix 125 μL、DNA 模板 2 μL、ddH2O 85 μL、引物各 1 μL。PCR 反應條件及引物序列見表2。擴增結束后，用1% 瓊脂糖凝膠電泳檢測，凝膠成像分析系統分析。

23純化測序PCR產物送上海生工進行純化后，測序。測序工作由上海生物的ABI-PRISM3730DNA自動測序儀完成。每個產物樣品采用雙向測序。

24數據處理測序峰圖采用Contigexpress序列拼接軟件校對拼接，去除低質量序列及引物區，利用MEGA 50對拼接后序列進行手工排序矯正后，比對分析；同時從 GenBank 數據庫中搜集并下載未獲得的ITS、 psbA-trnH、rbcL、matK序基因序列；共計221條（實驗獲得和GenBank下載）；運用 MEGA 50 軟件，以密花石豆蘭 （Bulbophyllum odoratissimum） 為外類群[5]，并基于K2P 距離模型構建NJ系統發育樹[6]，計算種間種內的K2P遺傳距離；采用Taxon DNA軟件做出各候選序列的barcoding gap圖。

3結果與分析

31序列信息、PCR 擴增效率和測序成功率PCR 擴增效率和測序成功率是評價DNA 條形碼的一個重要指標，本實驗對4個候選序列的 PCR擴增效率（出現明顯 PCR 條帶即判定成功）和測序成功率（測序重復2次后獲得高質量序列即為成功）進行了統計。其中rbcL序列由于長度較長，其擴增和測序成功率都非常底，不到70%和60%。psbA-trnH序列的擴增成功率為100%，但其測序成功率較低，僅為67%。matK序列的擴增成功率為100%，其測序成功率為89%。ITS序列無論擴增還是測序其成功率都最高，都接近100%。另外，各序列 GC 含量差異較大，ITS序列最高，平均含量達到 537%，matK序列GC 含量較低，平均僅為326%。詳細信息見表3。

32不同DNA條形碼序列種類種間差距種內種間差異也是評價條形碼序列的重要指標。合適的條形碼序列其種間變異要大于種內變異，以便準確鑒別同屬不同種的物種[7]。研究基于K2P模型，計算了14個石斛品種的種間種內遺傳距離 （表4）。結果顯示，4條序列中，ITS序列的種間遺傳距離為001～020，種內遺傳距離為000～008，除疊鞘石斛（Daurantiacum）以外，其余品種的種間遺傳距離都大于種內遺傳距離。而psbA-trnH、matK及rbcL序列的種間遺傳距離基本與種內遺傳距離相一致。見表4。

33不同序列 barcoding gap檢驗 Barcoding gap是指序列種內與種間遺傳變異由于顯著差異而形成的間隔區[8-9]，如圖1為4條序列的barcoding gap 圖。matK、psbA-trnH及rbcL序列均無明顯的間隔區。相對其他序列，ITS序列有較為明顯的間隔區，種內與種間變異明顯分為兩大部分，且其分布較高的遺傳變異也可明顯區分，因此在barcoding gap 分析方面，ITS序列為鑒別不同石斛較理想條形碼序列。

34NJ系統樹的構建與分析NJ 樹是基于 K2P 遺傳距離建立，通過分析不同物種間的親緣關系遠近以鑒別物種的系統發育樹，該方法可直觀展示不同物種的系統關系。如圖（2-5）所示，以密花石豆蘭 （Bulbophyllum odoratissimum） 為外類群[5]，基于K2P 距離模型構建NJ系統發育樹[6]，結果顯示ITS序列可以將不同物種分到不同的分支上，將不同品種的石斛區分開來，而matK、 psbA-trnH及rbcL序列均不能將不同品種的石斛完全具為1個分支。

4討論

DNA 條形碼是指利用有足夠變異、易擴增且相對較短的標準DNA片段，在種內特異性與種間多樣性中建立的一種新的生物身份識別系統，從而實現對物種進行快速、準確地識別和鑒定[10]。理想的條形碼序列應具有足夠的種間差異（種間差異應小于種內差異），易于擴增，且引物通用性好，能夠通過雙向測序得到高質量的序列[11]。本實驗選取了當前在植物學界推薦的4個熱點候選序列（ITS、matK、rbcL、psbA-trnH序列）對主要的14個滇產石斛品種進行鑒別[12]，四個片段的鑒別能力各有差別。

41matK序列對不同品種石斛的鑒定能力評價matK及片段有較高的的擴增和測序成功率，分別為100%和89%。計算K2P遺傳距離表明，其種內和種間差異保守。barcoding gap分析顯示，并沒有明顯的barcoding gap。NJ樹也不能將不同品種的石斛很好的聚為一支（圖3）。同時，其作為條形碼的最主要爭議是引物通用性差[13-14]，在不同物種間的通用性較小。鑒于此，本實驗中不建議選用該片段作為鑒定片段。

42psbA-trnH序列對不同品種石斛的鑒定能力評價psbA-trnH片段的PCR擴增和測序成功率，分別為100%和67%。根據barcoding gap分析顯示，其分布較高的遺傳變異與分布較低的遺傳變異可以明顯區分。但K2P遺傳距離表明，其種間遺傳距離沒有明顯大于種內遺傳距離，不能將不同種間的差異表現出來。同時，NJ樹也不能將不同品種的石斛很好的聚為一支。鑒于此，本實施不推薦psbA-trnH序列作為區分這14個石斛品種的理想條形碼片段。

43rbcL序列對不同品種石斛的鑒定能力評價與其它三個片段相比較，rbcL的序列長度最長，在本實驗中，其最長片段長度接近1400bp。因此，其擴增和測序的成功率都不高，僅為70%和59%。計算其K2P遺傳距離表明，其種間種內遺傳距離無明顯區別（000～001）。基于K2P遺傳距離構建的NJ樹也不能將不同品種的石斛很好的聚為一支。文獻研究報道[15]，rbcL的變異主要存在于種以上的水平，其鑒別主要應用在科屬及以上水平。鑒于此，rbcL也并不是鑒別該14個品種石斛的最佳條形碼。

44ITS序列對不同品種石斛的鑒定能力評價對于ITS片段，由于具有序列長度較短，易于擴增，對樣本要求較低。在本實驗中，其擴增和測序成功率均接近100%。根據實驗數據分析，除了疊鞘石斛（D aurantiacum）的K2P種內遺傳距離為006，大于種間遺傳距離，其余各品種的種間遺傳距離均大于種內遺傳距離。使用Mega 50軟件進行1000次自展計算，對14個石斛品種進行NJ樹的構建，由圖2可以看出，以密花石豆蘭 （Bulbophyllum odoratissimum）為外類群[6]，基于ITS序列構建的NJ系統聚類樹，不同的石斛品種基本上都聚在一起。同時barcoding gap分析顯示，其種內與種間變異區分為兩大部分，有明顯的barcoding gap，其分布較高的遺傳變異也可以明顯區分。綜合以上分析，ITS序列可以作為鑒別不同石斛物種的優選序列。

綜上，本實驗篩選了4個熱點DNA條形碼片段，用以區分云南市場上流通量較大的14個石斛品種。實驗結果表明，ITS序列是鑒別這14個石斛品種較為理想的DNA條形碼片段。ITS序列作為國際上承認的通用核心條形碼，在植物類的鑒別中有著廣泛的應用前景。

彭祿[16]等通過ITS序列鑒定了獨活的17個種，得到四帶芹類可作為獨活的替補品的結論。劉靜[17]等研究也發現，ITS序列可將17種藥用石斛完全鑒別開來。

ITS條形碼在中藥鑒定領域的廣泛應用將實現對中藥原植物及其藥材和飲片的準確快速鑒。

參考文獻

[1]郁美娟，孟慶華，黃德音，等. 石斛屬植物有效成分及藥理作用研究 [J]. 中成藥，2003，25（11）：918-921.

[2]國家藥典委員會.中華人民共和國藥典（第一部）[S]. 北京：中國醫藥科技出版社，2015：93-94.

[3]云南省藥材公司.云南中藥資源名錄 [M]. 北京：科學出版社，1993：669-672.

[4]Adams P.B. Systematics of Dendrobiinae （Orchidaceae）， with special reference to Australian taxa [J]. Bot J Linn Soc， 2011， 166： 105-126.

[5]Xiang X. G.， Schuiteman A.， Li D.Z.， et al. Molecular systematic of Dendrobium （Orchidaceae， Dendrobieae） from mainland Asia based on plastid and nuclear sequences [J]. Mol Phylogenet Evol， 2013， 69： 950-960.

[6]Chen S. L.， Yao H.， Han J. P.， et al. Validation of the ITS2 regionas a novel DNA barcode for identifying medicinal plant species [J]. PLoS One， 2010， 5： e8613.

[7]鄭司浩，李亞康，任偉光，等.基于DNA條形碼的百合屬分子鑒定[J].藥學學報，2014，49（12）：1730-1738.

[8]Meyer， C. P.， Paulay， G. DNA barcoding： error rates based on comprehensive sampling [J]. PLoS Biol， 2005， 3： e422.

[9]Lahaye， R.， Bank， M.， Bogarin， D.， et al. DNA barcoding the floras of biodiversity hotspots [J]. Proc Natl Acad Sci USA， 2008， 105： 2923-2928.[ZK）]

[10]Yuan， Y. W. Olmstead， R. G. Evolution and phylogenetic utility of the PHOT gene duplicates in Verbena complex （Verbenaceae）： dramtic intron size variation and footprint of abcestral recombination [J]. American jouranal of bantany， 2008， 95（9）， 1166-1176.

[11]陳士林，宋經元，姚輝. 藥用植物DNA 條形碼鑒定策略及關鍵技術分析 [J]. 中國天然藥物，2009，7 （5）：322-327.

[12]陳士林，姚輝，宋經元，等. 基于DNA barcoding （條形碼） 技術的中藥材鑒定 [J]. 世界科學技術-中醫藥現代化，2007，9（3）：7-12.

[13]Chase， M. W.， Cowan， R. S.， Hollingsworth， P.M.， et al. Hollingsworth ML. A proposal for a standardized protocol to barcode all land plants [J]. Taxon. 2007， 56： 295-299.

[14]Hollingsworth， P. M. DNA barcoding plants in biodiversity hot spots： Progress and outstanding questions [J].Heredity， 2008， 101： 1-2.

[15]Newmaster， S. G.， Fazekas， A. J.， Ragupathy， S. DNA barcoding in land plants： evaluation of rbcL in a multigenetiered approach [J]. Can J Bot，2008， 84： 335-341.

[16] 彭祿，余言，王志新，等. 基于ITS序列對獨活17個種的分子鑒定 [J]. 中草藥，2013，43（12）：1648-1653.

[17] 劉靜，何濤，淳澤. 基于ITS序列的中國藥用石斛及其混偽品的分子鑒定 [J]. 中國中藥雜志，2009，34（22）：2853-2856.

（編輯：梁志慶）