陳月娟+劉文簡+陳丹娜+詹賢福+蔣林

[摘要] 應(yīng)用Cytb序列對燕窩進行DNA條形碼（DNA barcodes）研究，為燕窩的溯源研究及全面評價燕窩質(zhì)量提供理論依據(jù)。以馬來西亞、印度尼西亞、越南、泰國等地的燕窩，共計39個樣本為研究對象，提取基因組DNA，擴增Cytb序列并雙向測序。采用MEGA 7.0軟件，對39條樣本序列及36條從Gen Bank中下載5種金絲燕的Cytb序列進行序列特征分析，基于Kimura雙參數(shù)法計算序列種內(nèi)、種間遺傳距離，構(gòu)建NJ和UPGMA系統(tǒng)聚類樹進行分析。研究結(jié)果表明，39個樣本共來源于3個不同的金絲燕物種。Cytb序列種間變異顯著大于其種內(nèi)變異，可以準確地區(qū)分不同種類的燕窩。所構(gòu)建的NJ和UPGMA系統(tǒng)發(fā)育樹，也能將全部物種區(qū)分開。表明基于Cytb序列的DNA條形碼技術(shù)可用于準確快速地鑒別燕窩的生物基原。

[關(guān)鍵詞] 燕窩； DNA條形碼； Cytb序列； 物種鑒定

[Abstract] To provide theoretical basis for the traceability and quality evaluation of edible bird′s nests （EBNs）， the Cytb sequence was applied to identify the origin of EBNs. A total of 39 experiment samples were collected from Malaysia， Indonesia， Vietnam and Thailand. Genomic DNA was extracted for the PCR reaction. The amplified products were sequenced. 36 sequences were downloaded from Gen Bank including edible nest swiftlet， black nest swiftlet， mascarene swiftlet， pacific swiftlet and germain′s swiftlet. MEGA 7.0 was used to analyze the distinction of sequences by the method of calculating the distances in intraspecific and interspecific divergences and constructing NJ and UPMGA phylogenetic tree based on Kimera-2-parameter model. The results showed that 39 samples were from three kinds of EBNs. Interspecific divergences were significantly greater than the intraspecific one. Samples could be successfully distinguished by NJ and UPMGA phylogenetic tree. In conclusion， Cytb sequence could be used to distinguish the origin of EBNs and it is efficient for tracing the origin species of EBNs.

[Key words] edible bird′s nests； DNA barcodes； Cytb； species identification

燕窩為雨燕科Apodidae金絲燕屬Aerodramus和侏金絲燕屬Collocalia的幾種燕類分泌的唾液與其絨羽混合筑成的巢窩，按筑巢的地方可分為“屋燕”和“洞燕”，主要產(chǎn)于印度尼西亞、馬來西亞等東南亞國家[1-2]。燕窩，性平、味甘，歸心、肺、腎三經(jīng)，是一種宜藥宜膳的高級保健品[3-4]，具有豐富的營養(yǎng)價值和藥用價值。中醫(yī)認為燕窩具有養(yǎng)陰潤燥、益氣補中、化痰止咳之功效[5]。古代醫(yī)書《本草求真》中記載[6]：“燕窩入肺生氣，入腎滋水，入胃補中，俾其補不致燥，潤不致滯，而為藥中至平至美之味者也”。燕窩中含有大量蛋白質(zhì)及生物活性物質(zhì)[7]，近年來研究表明，燕窩具有抗衰老、提高免疫力、抗病毒、強骨作用等功效[8-10]。

燕窩價格昂貴，市場需求量巨大，利益豐厚。它的價格和品質(zhì)不僅受燕窩顏色、成分、形狀等影響，而且不同種屬金絲燕所產(chǎn)的燕窩種類不同，在價格和質(zhì)量上也差異巨大[11]。傳統(tǒng)意義上認為燕窩的質(zhì)量與產(chǎn)地有密切的關(guān)系，如越南會安的燕窩優(yōu)于印尼、泰國等地的燕窩，這可能與不同產(chǎn)地金絲燕物種不同有關(guān)[4]。鳥類分類學(xué)家將產(chǎn)燕窩的金絲燕分為兩大類：一是雨燕科侏金絲燕屬Collocalia，包括白腹金絲燕C.esculenta、穴金絲燕C.linchi；二是雨燕科金絲燕屬Aerodramus，包括爪哇金絲燕A. fuciphagus、戈式金絲燕A.germani、大金絲燕A.maximus、印度金絲燕A.sunicolor、小灰腰金絲燕A.francicus[12]。目前國內(nèi)外對燕窩生物基原方面的研究報道較少，市場上以劣充優(yōu)的現(xiàn)象時有發(fā)生。本研究利用DNA條形碼鑒定技術(shù)，對39個未知來源的毛燕窩進行鑒定，為燕窩的溯源研究及全面評價燕窩質(zhì)量提供參考。

1 材料

ESP-600型電泳儀，C1000TM PCR儀，GBOX-F3凝膠成像儀，Centrifuge 5417R臺式高速冷凍離心機，Nanodrop2000微量定量儀口腔拭子基因組DNA提取試劑盒[天根生化科技（北京）有限公司]，引物（上海捷瑞生物工程有限公司廣州分公司），2×Taq PCR StarMix with Loading Dye、瓊脂糖、50×TAE（北京康潤誠業(yè)生物科技有限公司），Gold View I型核酸染色劑，DNA Marker（上海捷瑞生物工程有限公司），其余試劑均為分析純。endprint

本研究所有樣品均經(jīng)中山大學(xué)藥學(xué)院蔣林研究員鑒定，憑證樣本保存于中山大學(xué)藥學(xué)院。實驗樣本信息見表1，從Genbank下載36條Cytb序列，Gen Bank登錄號等信息見表2。

2 方法

2.1 總DNA提取 將燕窩樣品于50 ℃烘48 h，研磨成細粉，過120目篩，密封保存于-20 ℃冰箱備用。采用口腔拭子基因組DNA提取試劑盒提取總DNA。稱取燕窩樣品約40 mg，按照試劑盒操作說明進行總DNA的提取，利用微量定量儀對提取DNA的濃度和純度進行檢測。

2.2 PCR擴增及擴增產(chǎn)物的檢測 擴增序列的引物為Genbank中上傳最多的Cytb基因序列[13-14]，見表3。PCR反應(yīng)體系為20 μL：其中DNA模板1 μL，正向引物（10 μmol·L-1） 1 μL，反向引物（10 μmol·L-1） 1 μL，2×Taq PCR StarMix 10 μL， ddH2O 7 μL。PCR反應(yīng)條件的設(shè)置：94 ℃，2 min；94 ℃，30 s；55～65 ℃，30 s；72 ℃，1 min，35個循環(huán)；72 ℃，5 min。取5 μL的PCR 產(chǎn)物，用1.5%的瓊脂糖凝膠進行電泳。Marker采用上海捷瑞DL 2501，Marker相對分子質(zhì)量標準從上至下為2 000，1 000，750，500，250，100 bp。將檢測合格的PCR擴增產(chǎn)物送至生工生物工程有限公司進行雙向測序。

2.3 數(shù)據(jù)處理 測序結(jié)果用DNASTAR軟件包中的Seqman和Editseq進行校對拼接，去除引物區(qū)和低質(zhì)量的序列，然后通過Clustal X V2.1[15]校對，MEGA 7.0[16]對齊刪減，最終選取所有樣品長度均為530 bp的一段Cytb序列。從Genbank下載的36條Cytb序列，與拼接好的樣品序列比對后，對齊刪減獲得相對應(yīng)的區(qū)段。所有序列應(yīng)用MEGA 7.0軟件比對并計算K2P（Kimura-2-parameter）遺傳距離[17]，建立NJ系統(tǒng)聚類樹[18]和UPGMA系統(tǒng)聚類樹[19]。在Gen Bank上通過BLAST相似性搜索法[20]，進行序列的比對分析，驗證各樣品的生物基原。

3 結(jié)果與分析

3.1 樣品序列特征分析 利用MEGA 7.0軟件，將以Cytb為引物所獲得的39條序列進行比對，最終選取對齊后的530 bp的序列進行分析。所有序列不含插入或缺失堿基以及終止密碼子。結(jié)果顯示，序列中保守位點516個，變異位點14個，簡約信息位點7個，單一多態(tài)位點7個。A，T，C，G堿基平均量分別為27.3%，23.2%，35.1%，14.4%。

3.2 種內(nèi)種間遺傳距離分析 在軟件MEGA 7.0中應(yīng)用Kimura雙參數(shù)法計算39條樣品序列和36條從Genbank下載的5種金絲燕的Cytb序列的種內(nèi)遺傳距離和種間遺傳距離。結(jié)果顯示，大金絲燕和小灰腰金絲燕與其他金絲燕種間遺傳距離在1.85～1.92，其他3種金絲燕種間遺傳變異在0.02～0.08，種內(nèi)遺傳距離在0.000 3～0.023 6，見表4。種間遺傳距離明顯高于種內(nèi)遺傳距離，表明各序列的種間變異均顯著大于種內(nèi)變異，這與條形碼種間變異大而種內(nèi)變異小的要求是相一致的[18]，說明Cytb條形碼序列適用于燕窩生物基原的鑒定，并且可以對其進行有效的區(qū)分。

3.3 NJ系統(tǒng)聚類樹和UPGMA系統(tǒng)聚類樹 基于Cytb序列構(gòu)建的NJ系統(tǒng)聚類樹和UPGMA系統(tǒng)聚類樹見圖1，2，結(jié)果顯示2種方法獲得的物種聚集信息分布基本一致。所有序列被聚為兩大支，各個序列的建樹效果良好，基本上能將不同種類燕窩區(qū)分開。

UPGMA系統(tǒng)聚類樹節(jié)點的支持率明顯高于NJ系統(tǒng)聚類樹，以UPGMA系統(tǒng)聚類樹進行分析。爪哇金絲燕、白腰雨燕和戈氏金絲燕聚為一大支，節(jié)點支持率為99%。白腰雨燕單獨聚為一小支，爪哇金絲燕和戈氏金絲燕聚為一小支，同時又各自聚為不同的小支。大金絲燕和小灰腰金絲燕聚為一大支，節(jié)點支持率均為99%，大金絲燕和小灰腰金絲燕又各自聚為一小支。06號樣品與小灰腰金絲燕聚為一大支，節(jié)點支持率為89%，表明樣品生物基原為小灰腰金絲燕。39號樣品與戈式金絲燕聚為一支，節(jié)點支持率為56%，其生物基原為戈式金絲燕。其余樣品均與爪哇金絲燕聚為一支，表明其生物基原為爪哇金絲燕。以上結(jié)果與Blast搜索結(jié)果一致，表明Cytb序列可對不同種類的燕窩生物基原進行準確鑒別。

4 討論

DNA條形碼是近年來生物分類和鑒定的研究熱點和方向，是利用基因組中一段公認的標準短序列來進行物種鑒定的分子診斷新技術(shù)[21]。DNA條形碼研究的主要目的是用于物種的鑒定和新種的發(fā)現(xiàn)[22-23]，特別是在鑒定未知來源樣品方面有著極大 的優(yōu)勢[24]。近年來在經(jīng)過加工的肉類、海鮮食品的物種溯源方面得到成功的應(yīng)用[25-26]，為食品質(zhì)量監(jiān)督提供保障。燕窩從采摘到加工成品，需要經(jīng)過復(fù)雜的挑毛、塑型等過程，成品的外形十分相似，無法通過形態(tài)特征等進行物種溯源鑒定。本文采用DNA條形碼技術(shù)，對39份燕窩樣品進行溯源研究。實驗設(shè)計Cytb序列作為擴增引物，擴增和測序成功的效率高達100%。對樣本Cytb序列進行BLAST比對，結(jié)果顯示其相似度為99%～100%，E值為0，為鑒定的準確性提供保障。

《中國藥典》2010年版[27]和2015年版[28]已正式收載了DNA分子鑒別方法及中藥材DNA條形碼分子鑒定法指導(dǎo)原則。DNA條形碼技術(shù)可以為動植物來源的中藥材提供快速、準確的鑒定[29]。它通過構(gòu)建NJ，UPGMA系統(tǒng)聚類樹、計算遺傳距離、分析變異位點等方法相結(jié)合準確鑒定出中藥材基原。本研究基于Cytb序列將DNA條形碼技術(shù)應(yīng)用于燕窩樣品生物基原的鑒定。結(jié)果表明：Cytb序列種間變異顯著大于其種內(nèi)變異，符合理想的DNA條形碼 基因片段的特點[30]。所構(gòu)建的NJ和UPGMA系統(tǒng)發(fā)育樹結(jié)果顯示所有序列被聚為兩大支，各個序列的建樹效果良好，且2種方法獲得的物種聚集信息分布基本一致，可以將全部物種區(qū)分開。DNA條形碼技術(shù)對燕窩樣品的鑒定成功為燕窩的溯源研究及全面評價燕窩質(zhì)量提供理論參考，同時也為合理開發(fā)燕窩資源，規(guī)范燕窩市場提供理論依據(jù)。endprint

[致謝] 感謝鐘元、周王曉云、郭襄、林國超、柯海燕、許旭、寶迪、周利和、丁香、陶敏等提供燕窩樣品。

[參考文獻]

[1] 于?；ǎ於孛?，周昱，等. 燕窩的研究現(xiàn)狀[J]. 食品安全質(zhì)量檢測學(xué)報，2015，6（1）：197.

[2] 蔣林. 解開燕窩密碼[M]. 廣州：廣東省地圖出版社，2013：6.

[3] 南京中醫(yī)藥大學(xué).中藥大辭典[M]. 上海：上?？茖W(xué)技術(shù)出版社，2006.

[4] 林潔茹，周華，賴小平，等. 燕窩DNA提取方法研究[J]. 世界科學(xué)技術(shù)——中醫(yī)藥現(xiàn)代化，2010，12（2）：202.

[5] 盧端萍，程佳華，陳碩，等. 不同產(chǎn)地燕窩中唾液酸含量的高效液相色譜法測定[J]. 時珍國醫(yī)國藥，2016，27（2）：371.

[6] 黃宮繡. 本草求真[M]. 上海：上?？茖W(xué)技術(shù)出版社，1959：25.

[7] Yang M， Cheung S H， Li S C， et al. Establishment of a holistic and scientific protocol for the authentication and quality assurance of edible bird′s nest[J]. Food Chem， 2014，151：271.

[8] Hou Z P， Imam M U， Ismail M， et al. Effects of edible bird′s on hippocampal and cortical neurodegeneration in ovariectomized rats[J]. Food Funct， 2015， 6： 1701.

[9] Guo C T， Takahashi T， Bukawa W， et al. Edible bird′s nest extract inhibits influenza virus infection[J]. Antivir Res， 2006， 70： 140.

[10] Chua K H， Lee T H， Nagandran K， et al. Edible bird′s nest extract as a chondro-protective agent for human chondrocytes isolated from osteoarthritic knee： in vitro study[J]. Bmc Complem Altern M， 2013， 13： 19.

[11] Koon L C， Cranbrook E. Swiftlets of Borneo： builders of edible nests[M]. Borneo： Natural History Publications， 2003.

[12] 蔣林.走進燕窩世界[M]. 廣州：廣東省地圖出版社，2016：9.

[13] Lin J R， Zhou H， Lai X P， et al. Genetic identification of edible bird′s nest based on mitochondrial DNA sequences[J]. Food Res Int， 2009， 42（8）： 1054 .

[14] 林潔茹.燕窩DNA基源鑒定及抗病毒作用研究[D]. 廣州：廣州中醫(yī)藥大學(xué)，2010.

[15] Larkin M A， Blackshields G， Brown N P， et al. Clustal W and Clustal X version 2.0[J]. Bioinform， 2007， 23： 2947.

[16] 王孟虎，許亮，康廷國，等. 基于COI序列的哈蟆油基原動物DNA條形碼鑒定研究[J]. 中國中藥雜志，2017，42（8）：1572.

[17] Laopichienpong N， Muangmai N， Supikamolseni A， et al. Assessment of snake DNA barcodes based on mitochondrial and Cytb genes revealed multiple putative cryptic species in Thailand[J]. Gene， 2016， 594（2）： 238.

[18] Khaksar R， Carlson T， Schaffner D W， et al. Unmasking seafood mislabeling in U.S. markets： DNA barcoding as a unique technology for food authentication and quality control[J]. Food Control， 2015，56： 71.

[19] 張爽，劉宇靖，吳沿勝，等. 罌粟屬植物核心DNA條形碼的篩選[J]. 中國中藥雜志，2015，8，40（15）：2964.

[20] Xiong X， Guardone L， Cornax M J， et al. DNA barcoding reveals substitution of sablefish （Anoplopoma fimbria） with Patagonian and Antarctic toothfish （Dissostichus eleginoides and Dissostichus mawsoni） in online market in China. How mislabeling opens door to IUU fishing[J]. Food Control， 2016， 70： 380.endprint

[21] 陳士林，郭寶林，張貴君，等. 中藥鑒定學(xué)新技術(shù)新方法研究進展[J]. 中國中藥雜志，2012，37（8）：1043.

[22] Hebert P D， Stoeckle M， Zemlak T， et al. Identification of birds through DNA barcodes [J]. PLoS Biol， 2004， 2： 1657.

[23] Hebert P D N， Penton E H， Burns J M， et al. Ten species in one： DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator [J]. Natl Acad Sci USA， 2004， 101（41）： 14812.

[24] 陳士林，宋經(jīng)元，姚輝，等. 藥用植物 DNA 條形碼鑒定策略及關(guān)鍵技術(shù)分析[J]. 中國天然藥物，2009，7（5）：322.

[25] Willette D A， Simmonds S E， Cheng S H， et al. Using DNA barcoding to track seafood mislabeling in Los Angeles restaurants [J]. Conserv Biol， 2017， 4： 1.

[26] Quinto C A， Tinoco R， Hellberg R S. DNA barcoding reveals mislabeling of game meat species on the U.S. commercial market [J]. Food Control， 2016， 59： 386.

[27] 中國藥典. 一部[S].2010.

[28] 中國藥典. 一部[S].2015.

[29] 熊波，趙志禮，倪梁紅，等. DNA條形碼技術(shù)在中藥中的應(yīng)用、局限性與展望[J]. 中藥材，2015，38（10）：2202.

[30] 楊鵬，沈文華，石建明，等. 唇形科常見藥用植物DNA條形碼的鑒定研究[J]. 中草藥，2017，48（7）：1397.

[責(zé)任編輯 呂冬梅]endprint