免费文献传递   相关文献

Application of DNA barcoding in cannabic identification

DNA条形码技术在毒品原植物大麻鉴定中的应用



全 文 :广 西 植 物 Guihaia Jul.2014,34(4):552-556           http://journal.gxzw.gxib.cn 
DOI:10.3969/j.issn.1000G3142.2014.04.024
宋炳轲,杨雪莹,倪萍娅,等.DNA条形码技术在毒品原植物大麻鉴定中的应用[J].广西植物,2014,34(4):552-556
SongBK,YangXY,NiPY,etal.ApplicationofDNAbarcodingincannabicidentification[J].Guihaia,2014,34(4):552-556
DNA条形码技术在毒品原植物大麻鉴定中的应用
宋炳轲1,杨雪莹2,倪萍娅3,裴 黎2∗,张 颖2,徐小玉2
(1.中国人民公安大学,北京100038;2.公安部物证鉴定中心,北京100038;3.杭州市公安局上城区分局,杭州311200)
摘 要:准确鉴定毒品原植物大麻的种属及品种具有重要的理论和实践意义.为了探讨DNA条形码技术用
于毒品原植物大麻种属鉴定及品种鉴定的可行性,该研究以60份大麻原植物(分别采自内蒙、黑龙江、陕西延
安、陕西榆林4个地区的栽培大麻雌雄各6株及新疆玛纳斯地区的野生大麻雌雄各6株)为材料,通过从其叶
片中提取的DNA为模版,利用核糖体DNA基因间隔区的通用引物ITS2和叶绿体DNA的通用引物psbAG
trnH进行PCR扩增,对扩增片段进行双向测序,将测序结果进行人工矫正和比对.结果显示:所有大麻样本
的ITS2扩增片段序列没有变异完全一致,但psbAGtrnH扩增片段变异较大共检测出8种cpDNA单倍型,用
MEGE5.1软件计算种间遗传距离,并构建NJ系统聚类树可以有效把这五个地区的大麻样本区别开来,因此
证明DNA条形码技术在毒品原植物大麻的种属鉴定方面具有可行性,但其用于大麻的种属鉴定的准确性、可
靠性及在其来源地鉴定及品种鉴定中的可能性还有待进一步深入地研究.
关键词:DNA条形码;大麻;种属鉴定
中图分类号:795  文献标识码:A  文章编号:1000G3142(2014)04G0552G05
ApplicationofDNAbarcodingincannabicidentification
SONGBingGKe1,YANGXueGYing2,NIPingGYa3,
PEILi2∗,ZHANGYing2,XUXiaoGYu2
(1.ChinesePeople’sPublicSecurityUniversity,Beijing100038,China;2.InstituteofForensicScience,MinistryofPublic
Security,Beijing100038,China;3.ShangChengDistrictBranchofHangzhouPublicSecurityBureau,Hangzhou311200,China)
Abstract:TostudythefeasibilityofDNAbarcodingoncannabicspeciesidentification,DNAwereextractedfromsixty
cannabisplants’leaves(cultivatedcannabisfromInnerMongolia,Heilongjiang,ShanxiYan’anandYulin,wildcannaG
bisfromXinjiangManas,sixmaleandfemalecannabisfromeachregion),usedribosomalDNAintergenicregion
ITS2universalprimersandcpDNApsbAGtrnHprimerforPCRamplification,thensequencedinbothdirections,the
resultscomparedtomanualcorrectionandfinalyforBlastcomparison.ITS2’amplifiedsequencesofalsampleswere
completelyconsistent,butpsbAGtrnH’amplifiedsequencesvariedgreatly,whichweredetectedeightkindsofcpDNA
haplotypes.MEGE5.1softwarewasusedtocalculatethegeneticdistancebetweenspecies,andbuildNJphylogenetic
treeswhichcouldefectivelyseparatethefiveregionsofcannabicsamples.SoDNAbarcodingoncannabicspecies
identificationwasprovedtohaveviability,buttheaccuracy,reliabilityandpossibilitiesofidentificationofitsoriginand
speciesneedfurtherresearch.
Keywords:DNAbarcoding;cannabis;speciesidentification
收稿日期:2014G01G19  修回日期:2014G03G08
基金项目:公安部物证鉴定中心基本科研业务费专项(2010JB002)
作者简介:宋炳轲(1987G),女,河南禹州市人,硕士研究生,主要从事法医遗传学的研究,(EGmail)sbk1987@126.com.
∗通讯作者:裴黎,主任法医师,硕士生导师,主要从事分子生物学与法医遗传学研究,(EGmail)peili128@sina.com.
  大麻(Cannabissativa)为被子植物门双子叶植
物纲荨麻目大麻科大麻属1~2年生雌雄异株的草
本植物.大麻属植物包括栽培大麻、印度大麻和野
生大麻三个主要品种及在长期生物进化过程中产生
的大麻变种和亚种(卢延旭等,2007).大麻植物所
含的四氢大麻酚(THC)具有致幻作用,所以大麻被
列为与海洛因、可卡因并列的三大毒品之一(谢幸媚
等,2000).不同品种大麻中THC的含量有很大的
差别且大麻是许多国家的重要经济作物,包括我国
在内的许多国家根据大麻体内THC含量的高低,
将大麻分为工业大麻、毒品大麻和中间型大麻,其中
后两种被列为严厉打击的对象.然而,在高额利润
的驱使下,全球范围的大麻种植和非法交易活动屡
禁不止.与此同时,在大麻犯罪案件中缴获成批未
加工的大麻植物(如整株植物、种子等)的情况也屡
见不鲜,在法庭科学中,虽然可用传统的形态学方法
对大麻进行鉴定,用理化方法如常规化学法、毛细管
电泳法、薄层色谱法等对大麻毒性成分进行测定,但
这些方法都有很大的局限性,对检材的要求较高且
很难或根本无法区分大麻的品种及来源.为了有效
打击大麻毒品犯罪,法庭科学急需寻找快速准确对
大麻进行种属鉴定和品种鉴定的新方法.
DNA条形码(DNAbarcoding)这一概念由加
拿大分类学家PaulHebert于2003年首次提出,随
后在世界范围内受到了广泛关注,成为生物学研究
新热点.它是通过对一个或多个标准目的基因片段
进行序列分析达到物种鉴定目的,是分类学中辅助
物种鉴定的新技术.目前该技术已成为生物分类学
研究中的一个新方向,在医药、生态、环境、食品等领
域有广泛应用(Nacirietal.,2012;Kressetal.,
2012).植物DNA条形码筛选区域主要集中在叶
绿体和核糖体内转录间隔区,主要包括 rpoB,
rpoCl,matK,rbcL,UP,trnHGpsbA,atpFGatpH,psG
bKGpsbI及核基因ITS等(Kressetal.,2012;Liet
al.2011;任保青等,2010).但目前尚未有获得广泛
认同的植物DNA条形码.2009年在墨西哥召开的
第三届DNA条形码国际学术大会上,与会者达成
共识一致建议选择进化速率较快的ITS和trnHG
psbA,探讨它们在陆地植物中的通用性、分辨率和
适合程度.ITS即核糖体DNA基因间隔区,组成
该片段的不同部分(ITSl、ITS2和5.8S)序列变异差
别较大,ITS1片段生物识别率高于ITS2,但研究发
现ITS2片段扩增成功率高于ITS1,而且植物叶绿
体DNA属单系遗传(Chenetal.,2010;Yaoetal.,
2010;Schochetal.,2012).通过全面综合考虑,我
们最终选择ITS2和psbAGtrnH作为初步研究我国
不同地区大麻变异情况的对象,对DNA条形码技
术在毒品原植物大麻鉴定中的应用可行性进行了初
步研究,以期为今后进一步的工作奠定基础.
1 材料与方法
1.1材料
采自内蒙古、黑龙江、新疆玛纳斯和陕西延安、
陕西榆林五个地区的大麻样本60份,每个地区雌雄
各六份,均由公安部物证鉴定中心提供.
1.2大麻全基因组DNA的提取与定量
取干叶10 mg,利用植物 DNA 提取试剂盒
(TiangenBiotechCo.,China)提取大麻全基因组
DNA,用 NanoDrop_NDG2000Spectrophotometer
(ThermoFisherScientific,Wilmington,DE)核酸
定量仪进行定量,定量结果见表1.
表1 大麻全基因组DNA定量结果
Table1  Quantitativeresults
样本采集地
Samplecolectingzone
OD260/280
平均浓度(ng/μL)
Meanconcentration
新疆玛纳斯 ManasiXinjiang 2.01 85.4
内蒙古InnerMongolia 1.98 110.2
黑龙江 Heilongjiang 1.75 93.6
陕西延安 Yan’an 1.81 100.1
陕西榆林ShanxiYulin 1.87 112.2
1.3PCR扩增与琼脂糖凝胶电泳
采用核糖体 DNA 基因间隔区的通用引物
ITS2(朱英杰等,2010)和叶绿体DNA的通用引物
psbAGtrnH(刘莹等,2013)进行扩增(引物序列见表
2,由上海英潍捷基公司合成),扩增反应在eppenG
dorf扩增仪上进行,PCR反应总体系为50μL,内含
2×PCRMastermix25μL,正反向引物(10mmol/
L)各1μL,去离子水23μL,模版DNA1μL.PCR
反应程序为94℃预变性5min;94℃变性30s,58
℃退火30s,72℃延伸45s,共30个循环;72℃终
延伸20min.PCR产物用2%琼脂糖进行凝胶电泳
检测.
1.4扩增产物测序
对扩增产物进行双向测序,测序反应由上海生
物工程公司完成.
3554期        宋炳轲等:DNA条形码技术在毒品原植物大麻鉴定中的应用
表2 ITS2及psbAGtrnH通用引物序列
Table2 Thesequencesofuniversal
primerITS2andpsbAGtrnH
引物编号PrimerNo. 引物序列Primersequences
ITS2 S2F
(正向引物) ATGCGATACTTGGTGTGAAT
S3R
(反向引物) GACGCTTCTCCAGACTACAAT
psbAGtrnH fwdPA
(正向引物) GTTATGCATGAACGTAATGCTA
RevTH
(反向引物)CGCGCAATGGTGGATTCACAATCC
1.5测序结果的拼接与分析
用DNAMANversion6软件及ChromasLITE
version2.01软件完成序列拼接,用 MEGA5.1软
件进行60份大麻样本的序列比对并进行变异位点
分析,然后登陆 http://blast.ncbi.nlm.nih.gov/
Blast.cgi网站,将所得序列进行blast比对,并利用
MEGA5.1进行多重比对,基于 KimuraG2GparameG
ter(K2P)双参数模型计算遗传距离,并用邻接法
(neighbourGjoining)构建系统聚类树(NJ树),同时
以bootstrap(自展支持率1000次)重复检验各分
支的支持率.
2 结果与分析
2.1扩增结果
实验结果表明ITS2和psbAGtrnH 对60份大
麻样本的扩增及测序成功率均为100%.
2.2测序及序列比对
通过测序得引物ITS2和psbAGtrnH的完整扩
增片段DNA序列,使用NCBI数据库中的Blast工
具进行搜索,比对后可以确定与该段序列相似度最
高的匹配均为大麻核糖体和叶绿体基因,匹配度分
别达到100%和99%,序列号分别为 KC292629.1
和KC578820.1.说明扩增序列即为所需要的大麻
核糖体和叶绿体基因的一部分.
2.3序列变异比较
应用 MEGE5.1进行序列对比,发现采自5个
不同地区的60份大麻样本的ITS2引物扩增区域
测序所得的序列长度均为446bp,且所有大麻样本
的序列完全一致,不含Poly结构,未发现变异位点;
引物psbAGtrnH 扩增片段的实际长度为404bp,共
发现11个单一多态性变异位点(分别为1、12、40、
47、202、284、330、377、381、391、404),占总分析位点
的2.72%,由于cpDNA常被看作是一个独立的遗
传单位,在细胞分裂时不经受重组(Schaaletal.,
1998).故对来自五个地区来源的60份大麻样本可
定义8种cpDNA单倍型(表3),其中新疆玛纳斯地
区的12份雌雄标本序列完全一致为 HA型,黑龙
江地区的六份雌性标本序列完全一致均为 HA型、
雄性均为HB型,内蒙地区的雌性序列完全一致均
为HC型、雄性均为HD型,陕西延安地区的雌性序
列完全一致为HE型、雄性为HF型,陕西榆林地区
的六份雌性标本序列完全一致为 HG型、雄性为
HH型.分析结果(图1).
表3 大麻的变异位点与单倍型
Table3 Variationpositionandhaplotypeofcannabissativa
引物序列长度
LengthofPrimer
sequences
H
A
H
B
H
C
H
D
H
E
H
F
H
G
H
H
1bp C C C A C A C C
12bp C T C C C C C T
40bp C T C C C C C C
47bp T A T T T T T T
202bp G G G G A G G G
284bp T A T T A T T T
330bp A A C C C A A A
377bp C C C G C C C C
381bp T T T G T T T T
391bp T G G G T T T T
404bp G G G G G T T G
图1 psbAGtrnH序列的变异位点
(“.”表示与黑龙江雌性碱基相同)
Fig.1 VariationpositionofpsbAGtrnHsequence
2.4遗传距离分析
将五个地区的大麻序列进行分析,根据碱基变
异位点的情况,用KimuraG2GParameter距离模式计
算的遗传距离(表4).五个地区的大麻样本遗传距
455 广 西 植 物                  34卷
表4 五个地区大麻遗传距离分析
Table4 GeneticdistanceanalysisofCannabicinfiveregions
黑龙江雌性
Heilongjiang
female
黑龙江雄性
Heilongjiang
male
内蒙雄性
InnerMongolia
male
内蒙雌性
InnerMongolia
female
新疆雌性
Xinjiang
female
新疆雄性
Xinjiang
male
延安雌性
Yan’an
female
延安雄性
Yan’an
male
榆林雌性
Yulin
female
黑龙江雌性
Heilongjiangfemale
0.0125
内蒙雌性
InnerMongoliafemale
0.005 0.0125
内蒙雄性
InnerMongoliamale
0.0125 0.0201 0.0075
新疆雌性
Xinjiangfemale
0 0.0125 0.005 0.0125
新疆雄性
Xinjiangmale
0 0.0125 0.005 0.0125 0
延安雌性
Yan’anfemale
0.005 0.0125 0.01 0.0175 0.005 0.005
延安雄性
Yan’anmale
0.005 0.0175 0.01 0.0125 0.005 0.005 0.01
榆林雌性
Yulinfemale
0.0025 0.015 0.0075 0.015 0.0025 0.0025 0.0075 0.0025
榆林雄性
Yulinmale
0.0025 0.01 0.0075 0.015 0.0025 0.0025 0.0075 0.0075 0.005
离在0.000~0.0201,平均遗传距离为0.008.其中
黑龙江的雄性大麻与内蒙的雄性大麻最大,为
0.0201,新疆的雄性雌性和黑龙江的雌性最小.
2.5系统发育邻接树分析
从图2可以看出,新疆的野生雌雄大麻与黑龙
江的雌性大麻聚为一支,陕西延安和榆林及黑龙江
地区的雄性大麻聚为一支,内蒙的雌性与雄性大麻
聚为一支,陕西延安和榆林的雌性大麻聚为一支.
说明sbAGtrnH条形码序列可用于不同地域大麻雌
雄间的鉴别分析.
图2 基于psbAGtrnH条形码序列构建的五个
地区雌雄大麻间的 NJ系统进化树
Fig.2 NJphylogenetictreeoffiveregionsbetweenmaleand
femalecannabisbasedonpsbAGtrnHbarcodesequences
3 讨论
大麻是一种集造纸、纺织、食用、药用等作用于
一身的经济价值很高的作物,具有抑制杂草生长、适
应环境能力强等特点,有利于戈壁滩涂和贫瘠土地
的绿化与开发利用.本实验结果表明,五个地区的
大麻样本的ITS2引物扩增区域测序所得的序列完
全一致,说明利用核糖体通用引物ITS2不能对不
同地域来源的大麻种属进行鉴定,虽然如此,但通过
blast比对却可以把大麻鉴定出来;而采用psbAG
trnH条形码序列对大麻研究的结果表明,其不但可
以把不同地域来源的大麻鉴定出来还可以把大麻的
雌性与雄性、野生大麻与栽培大麻分辨出来.实验
结果还表明,不同地区的雌雄大麻的cpDNA单配
体不同,而且由于大麻通常是雌雄异株,雄性植株的
纤维质量高于雌性植株,而雌性植株的有毒成分高
于雄性植株,雄株比雌株先成熟,对机械化的收割造
成不便;雌雄株在幼苗时期较难分辨,只有当其生殖
器官原基开始发育时才能正确识别(李仕金等,
2008).基于此,本研究一方面可以为大麻犯罪案件
中查找毒品来源、铲除毒品种植基地提供线索,有利
于打击个人及集团式的大麻贸易活动;另一方面也
为大麻育种提供了可能,而且有助于培育出低毒、高
纤维利用价值的雌雄同株的大麻新品种,以造福人
类和社会.
5554期        宋炳轲等:DNA条形码技术在毒品原植物大麻鉴定中的应用
4 结论
本实验是在有限的大麻标本中得出的结论,故
为了更好更精确地验证DNA条形码技术对毒品原
植物大麻的鉴定还需要一方面加大同一地区大麻雌
雄植株的采样,另一方面加大更多地区大麻品种的
研究;另外也需要进一步深入筛选更多更合适的条
形码序列.总的来说,DNA条形码技术有望成为将
来毒品案件刑事侦查和检验鉴定的新方法.
致谢 本论文由公安部物证鉴定中心基本科研
业务经费项目(2010JB002)资助,并得到了公安部
物证鉴定中心法医物证检验处的各位老师的指导和
帮助,在此深表感谢!
参考文献:
ChenSL,YaoH,HanJP,etal.2010.ValidationoftheITS2reG
gionasanovelDNAbarcodeforidentifyingmedicinalplantspeG
cies[J].PLoSone,5(1):8613
KressWJ,EricksDL.2012.DNABarcodes:MethodsandProtoG
cols[M].Berlin:Springer,858:3-8
LiDZ,LiuJQ,ChenZD,etal.2011.PlantDNAbarcodinginChiG
na[J].JSyst&Evol,49(3):165-168
LiSJ(李仕金),XinPY(辛培尧),ZhouJ(周军),etal.2008.ProG
gressesonrsearchofsexualdifferentiationofCannabissativaL.
(大麻性别研究进展)[J].PlantFiberSciChin(中国麻业
学),30(2):110-113
LuYX(卢延旭),DongP(董鹏),CuiXG(崔晓光),etal.2007.
DifferencebetweenindustrialhempandmarijuanahempandinG
dustrialhempcsusevalue(工业大麻与毒品大麻的区别及其可
利用价值)[J].ChinPharmacolBull(中国药理学通报),23
(8):1112-1114
LiuY(刘莹),ZhangAS(张安世),ZhaoLX(赵利新).2013.
Thephylogeneticrelationshipofsomebryophytesbasedon
trnHGPsbAgene(基于trnHGPsbA序列的部分苔藓植物亲
缘关系分析)[J].HubeiAgricSci(湖北农业科学),52
(20):5057-5059
NaciriY,CaetanoS,SalaminN.2012.PlantDNAbarcodesandthe
influenceofgeneflow[J].MolEcolResour,12(4):575-580
RenBQ(任保青),ChenZD(陈之端).2010.DNAbarcodingplant
life(植物DNA条形码技术)[J].ChinBullBot(植物学报),
45(1):1-12
SchochCL,SeifertKA,HuhndorfS,etal.2012.Nuclearribosomal
internaltranscribedspacer(ITS)regionasauniversalDNA
barcodemarkerforFungi[J].PNAS,109(16):6241-6246
SchaalBA,HayworthDA,OlsenKM,etal.1998.Phylogeographic
studiesinplants:problemsandprospects[J].MolEcol,7(4):
465-474
XieXM(谢幸媚),WuKL(吴坤龙),LiJY(李俊亿),etal.2000.
TheapplicationandevaluationofplantDNAtestingincriminal
investigations(植物DNA鉴定在犯罪侦查上之应用与评估)
[J].ForenSci(刑事科学),49:1-7
YaoH,SongJ,LiuC,etal.2010.UseofITS2Regionasthe
universal DNA Barcode for plants and animals [J].
PLoSone,5:13102
ZhuYJ(朱英杰),ChenSL(陈士林),YaoH(姚辉),etal.2010.
DNAbarcodingthemedicinalplantsofthegenusParis(重楼属
药用植物 DNA条形码鉴定研究)[J].ActaPharmSin(药学
学报),45(3):376-382
􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠􀤠
(上接第569页Continuefrompage569)
 Maxim.(狗枣猕猴桃果实发育的解剖学观察)[J].ChinBull
Bot(植物学通报),19(4):469-476
ChenYP(陈亚萍),ZhangYC(章英才),JingHX(景红霞),etal.
2013.Thecharacteristicofphotosynthesisandthegrowthoffruit
ofZiziphusjujubaMilcv.lingwuchangzao(灵武长枣光合特性
及果实生长发育特征)[J].NorthernHortic(北方园艺),(19):
1-4
LiuGC(刘国成),MaHY(马怀宇),LvDG(吕德国),etal.2012.
Researchonthechangeof‘Hanfu’applefruitanatomicstructure
andqualityduringstorage(“寒富”苹果贮藏期果实解剖结构及
品质变化研究)[J].NorthernHortic(北方园艺),(15):1-4
WangXL(王勋陵).1974.ApreliminarystudyofthedevelopmenG
talanatomyofChinesejujubefruit(枣果实发育解剖学的初步
观察)[J].ActaBotSin(植物学报),16(2):161-169
XieMH(颉敏华),ZhangHM(张永茂),LiSQ(李守强),etal.
2008.PostharestfreshGkeepingstoragecharacteristicsof‘lingG
wuchangzao︐(ZiziphusjujubaMil.)(灵武长枣采后保鲜贮藏
特性研究)[J].ActaBotBorGOccSin(西北植物学报),28(5):
1031-1035
YangJ(杨军),ZhangYC(章英才),SuWD(苏伟东),etal.2011.
IsolationandassayingofpolysaccharideinZiziphusjujuba Mil
cv.lingwuchangzao(灵武长枣多糖含量的测定研究)[J].
NorthernHortic(北方园艺),(14):35-37
YangSJ(杨淑娟),ZhangYC(章英才),ZhengGQ(郑国琦).2010.
Comparativestudyondissectedstructuresofnormalfruitand
crackingonesofLingwuLongGjujube(灵武长枣正常果与裂果
解剖结构的比较研究)[J].NorthernHortic(北方园艺),(22):
15-18
YangSJ(杨淑娟),ZhengGQ(郑国琦),ZhangYC(章英才).2011.
UltracytochemicallocalizationofCa2+ inthenormalfruitsand
crackingonesofLingwuLongGjujube(灵武长枣正常果及裂果
中Ca2+的细胞化学定位研究)[J].ActaBotBorGOccSin(西北
植物学报),31(1):0084-0088
YuJF(喻菊芳),WeiTJ(魏天军),ChenWJ(陈卫军).2008.The
researchofgermplasmresourcesandthestudyonvarietyselecG
tionofZiziphusjujubaMilcv.lingwuchangzao(灵武长枣种质
资源调查和品种选优研究)[J].ChinFruits(中国果树),(1):
56-57
YuJF(喻菊芳),ZhuLC(朱连成),WeiTJ(魏天军).2004.
Ningxia(Lingwu)LongGjujubetextualresearch(宁夏(灵武)长
枣考证)[J].NingxiaJAgricForSciTechnol(宁夏农林科
技),(5):31-32
ZhangYC(章英才),SuWD(苏伟东),YangJ(杨军).2012.The
storelocalizationcharacteristicsofpolysaccharideinZiziphusjuG
jubaMilcv.lingwuchangzao(灵武长枣多糖积累分布特征研
究)[J].NorthernHortic(北方园艺),(21):7-11
655 广 西 植 物                  34卷