G-B-T 43271-2023 木材鉴定 DNA条形码方法（正式版）

ICS79.010木材鉴定DNA条形码方法国家标准化管理委员会国家市场监督管理总局发布国家标准化管理委员会GB/T43271—2023 I Ⅱ 2规范性引用文件 2 2 6 附录A(资料性)木材DNA条形码鉴定方法的常用条形码类型及通用引物 附录B(资料性)木材DNA条形码的序列一致性分析 8附录C(资料性)木材DNA条形码序列的邻接树构建 9附录D(资料性)木材DNA条形码序列的最大似然树构建 IGB/T43271—2023本文件按照GB/T1.1—2020《标准化工作导则第1部分：标准化文件的结构和起草规则》的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由国家林业和草原局提出。本文件由全国木材标准化技术委员会(SAC/TC41)归口。科委司法鉴定中心、天津工业大学、中国检验认证集团天津有限公司、东阳市千百年红木家具有限公司、北京天怡利华沉香技术研究院、山东巧夺天工家具有限公司、浙江精典木材检测服务有限公司、浙江省木雕红木家具产品质量检验中心、山西紫檀臻品红木家具有限公司。ⅡGB/T43271—2023易执法和国内市场监督管理的迫切需求。DNA条形码方法为木材“种”水平鉴定提供了技术依据。业链监管水平，进一步推动我国林业产业健康可持续发展，特制定本文件。1GB/T43271—2023木材鉴定DNA条形码方法定报告。2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文本文件。GB/T6682分析实验室用水规格和试验方法3.13.2DNA条形码DNAbarcode生物体基因组中用于物种鉴定的一段标准的、相对较短的、易扩增的DNA(deoxyribonucleicacid,脱氧核糖核酸)序列。3.4亲缘关系相近的种的集合。从树木主干或枝条部位采集的以木质部为主的包含完整树木分类学及采集信息的植物样本。3.6刚伐倒后未经干燥的木材。2GB/T43271—2023种是心材和边材区别明显的树种。木材干燥wooddrying;woodseasoning在热能作用下以蒸发或沸腾方式排除木材水分的处理过程。3.10聚合酶链式反应polymerasechainreaction;PCR用于扩增位于两段已知序列之间DNA的分子生物学技术。3.11与特定DNA片段两端序列互补的人工合成的寡核苷酸短片段，用于聚合酶链式反应。3.123.13两条序列在同一位点核苷酸完全相同的部分所占的比例。3.14间系统发生关系的树状图。4原理DNA条形码方法是通过从未知样本(待检样本)中获取标准化DNA区域的核苷酸序列，然后使用生物信息学方法将该序列与参考数据库中的已知序列进行对比，最终根据对比结果将样本鉴定到物种水平。采用该方法的前提是条形码序列的种内差异5试剂5.1所使用的水均应符合GB/T6682的灭菌双蒸水。除另有规定外，所用试剂溶液宜大体积配制、小3GB/T43271—20235.2次氯酸钠(sodiumhypochlorite,NaClO)。5.3十六烷基三甲基溴化铵[hexadecyltrimethylammoniumbromide,CTAB,C₁₆H₃₃(CH₃)₃NBr]。5.4乙二胺四乙酸(ethylenediaminetetraaceticacid,EDTA,C₁₀H₁₆N₂O₈)。5.5三羟甲基氨基甲烷[tris(hydroxymethyl)aminomethane,C₄H₁₁NO₃]。5.6盐酸(hydrochloricacid,HCl)。5.7聚乙烯吡咯烷酮[polyvinylpyrrolidone,PVP-40,(C₆H₉NO)40]。5.8二硫苏糖醇(dithiothreitol,DTT,C₄H₁₀O₂S₂)。5.9乙酸钠(sodiumacetate,CH₃COONa)。5.10三氯甲烷(chloroform,CHCl₃)。5.11异戊醇(isoamylalcohol,C₅H₁₂O)。5.12三氯甲烷：异戊醇(体积比24:1)。5.13异丙醇(isopropanol,C₃H₈O)。5.1470%乙醇(ethanol,C₂H₆O)。5.15蛋白酶K(proteinaseK)。5.16PCR缓冲液。5.17TaqDNA聚合酶(TaqDNApolymerase)。5.18氯化镁(magnesiumchloride,MgCl₂)。5.19脱氧核糖核苷三磷酸(deoxyribonucleosidetriphosphate,dNTP)。5.20琼脂糖(agarose)。5.21电泳缓冲液TAE[由三羟甲基氨基甲烷、乙酸(aceticacid,C₂H₄O₂)和乙二胺四乙酸组成]。5.22核酸染料。6.1超净工作台(平均风速0.25m/s～0.6m/s)。6.4冰箱(—28℃~4℃)。6.6恒温水浴锅(工作范围30℃~70℃)。6.8低温高速离心机(可实现4℃离心，离心力9500g以上)。6.9分光光度计。6.10PCR仪。6.12凝胶成像系统。6.13DNA测序仪。7试验步骤7.1.1应采用规范的实验室操作措施有效避免或消除外源污染及交4GB/T43271—20237.1.3根据样品加工处理条件、保存状态、存储时间和取样部位的差异，通过优化DNA提取参数以提高DNA提取质量和效率。经干燥或长期存储的样品宜采用十六烷基三甲基溴化铵(CTAB)方法，生材样品宜采用DNA提取试剂盒方法。7.1.4木材DNA提取和PCR扩增试验步骤中应设置阴性对照。水率为12%时)。样品首先采用3%次氯酸钠浸泡30min,然后切除样品外层表面(1mm～3mm)以去一般来说，样品采集后尽快进行DNA条形码分析，否则宜储存在-80℃或液氮(-196℃)条件下。7.3DNA提取根据离心柱法或磁珠法DNA提取试剂盒说明书进行。操作步骤如下。a)细胞裂解：取约300mg木粉装入离心管中，分别加入3mL木材DNA裂解液[0.03g/mLCTAB,100mmol/LTris-HCl(pH8.0),1.4mol/LNaCl,20mmol/LEDTA,0.02g/mL~b)抽提：水浴结束后，向离心管中加入等体积三氯甲烷：异戊醇(24:1),颠倒混匀5min;9000g、4℃离心10min;吸取上层水相于另一无菌离心管中。重复此操作1次。c)沉淀：向装有抽提后水相的离心管中加入等体积预冷的异丙醇(1/10体积3mol/L的乙酸钠和终浓度为0.02mg/mL～1mg/mL的糖原),上下颠倒混匀，此时可能出现白色絮状沉淀，在—20℃冰箱中放置30min。d)洗涤：9000g、4℃离心10min,倒掉上清，用70%乙醇洗涤沉淀2次，置于室温下晾干。f)纯化：利用DNA模板纯化试剂盒进行木材DNA纯化。7.4DNA定量采用紫外分光光度法或荧光定量法等对木材DNA提取液进行浓度测定。7.5引物设计与PCR扩增7.5.1DNA条形码选择一般可选择rbcL、matK、psbA-trnH、trnL-trnF、ycfl等叶绿体DNA片段和细胞核核糖体DNA内转录间隔区ITS作为通用DNA条形码。对于特定树种类群，也可使用针对该类群开发的专属DNA条形码。5GB/T43271—2023引物序列应位于高度保守区，与非扩增区无同源序列；引物长度一般为15bp～30bp;碱基尽可能随机分布，GC含量宜为40%～60%;引物内部应避免形成稳定的二聚体或发夹结构；引物序列的3’端应避免出现碱基A及3个以上的连续碱基。常用DNA条形码类型及通用引物实例见附录A。PCR操作流程如下。a)反应体系：反应体系的体积可选择20μL～50pL。当PCR反应体系为20μL时，各组分的终浓度如表1所示。反应组分原浓度终浓度反应体积(20μL)ddH₂O10×PCR缓冲液MgCl₂25mmol/L2.5mmol/LDNTPs2.5mmol/L0.15mmol/LTaq酶正向引物反向引物b)反应程序：反应程序通常可采用：94℃预变性为90s,94℃变性5s,Tm-5℃退火30s(Tm为引物熔解温PCR完成后，可将PCR产物临时存放于4℃冰箱，对于需长期保存的(一周以上)PCR产物应置于—20℃冰箱保存。琼脂糖凝胶电泳检测流程如下：a)称取适量的琼脂糖并加入电泳缓冲液中，配制成质量浓度1%～2%的溶液，加热溶解；7.6DNA测序与序列拼接PCR产物纯化后，以PCR扩增引物作为测序引物对DNA片段进行测序。6对测序结果进行质量评估，去除序列两端Phred质量分数小于20的碱基和测序引物序列，然后将正向和反向互补序列组装成一个完整序列。7.7DNA条形码序列分析将获得的DNA序列与参考数据库进行比对，通过两两序列局部比对或者搜索短的核苷酸字符串来查询数据库中最匹配的树种序列。参考数据库应已包含相关树种的序列。序列一致性分析实例见附将获得的DNA序列与参考数据库序列进行多重序列比对，根据优选的建树方法构建分子系统发育树，并用重复抽样检验方法评估分子系统发育树的可靠性。分子系统发育树构建实例见附录C和附录D。当采用其他方法进行DNA序列分析时，应说明该方法的具体步骤。8判定当采用其他的序列分析方法时，应说明具体判定依据和判定结果。当DNA条形码序列分析无法判定相应树种，则应给出科、属等更高分类阶元的鉴定结果。9鉴定报告7GB/T43271—2023(资料性)木材DNA条形码鉴定方法的常用条形码类型及通用引物常用DNA条形码类型及通用引物见表A.1。表A.1木材DNA条形码鉴定方法的常用条形码类型及通用引物引物名称引物序列(5'-3')适用范围rbcLrbcLbFAGACCTWTTTGAAGAAGGTTCWGTrbcLbRTCGGTYAGAGCRGGCATRTGCCAmatKmatK472FCCCRTYCATCTGGAAATCTTGGTTC被子植物(含阔叶树)matK1248RGCTRTRATAATGAGAAAGATTTCTGCGym_F1AATYGYRCTTTTATGTTTACARGC裸子植物(含针叶树)Gym_R1ATCAYCCGGARATTTTGGTTCGpsbA-trnHtrnHf_05CGCGCATGGTGGATTCACAATCCpsbA3fGTTATGCATGAACGTAATGCTCycf1ycflbFTCTCGACGAAAATCAGATTGTTGTGAAT被子植物(含阔叶树)ycflbRATACATGTCAAAGTGATGGAAAAycflgFTGAAAGCTCTAAGCAATGGATCYCC裸子植物(含针叶树)ycflgFATACGACCAATATTTTTRGCTATTATATGCGATACTTGGTGTGAATGACGCTTCTCCAGACTACAATCCTTATCAYTTAGAGGAAGGAGYGACTCTCGGCAACGGATAGCGTTCAAAGAYTCGATGRTTCRGTTTCTTTTCCTCCGCTTA8GB/T43271—2023(资料性)木材DNA条形码的序列一致性分析以常用DNA数据库NCBI(NationalCenterforBiotechnologyInformation)为例。采用基于局部对比算法搜索工具(BLAST)的相似性搜索方法分析时，进入BLAST页面，将所获得的木材样品DNA条形码序列作为查询序列置于“输入查询序列(EnterQuerySequence)”窗口中；选择“选择搜索设置(ChooseSearchSet)”中“标准数据库(Standarddatabases)”,默认搜索范围为“核苷酸数据库(Nucleotidecollection)”;程序选择(ProgramSelection)中推荐优先选择“高度相似序列(Highlysimilarsequences)”;点击“BLAST”,获得最大一致性，确定匹配唯一性，完成序列一致性分析(见图B.1)。检素结果序列一致性DescriptionScientificCormmonNameNameMaxTotalQueryEPer.AAAcoessionvoucherJY10smallsubunitribosomalRNAgene.partialsequencs;intema…Phoebe.….NA9431011097210972100%0.0100.00%5941MT628602.1PhoabeboumeivoucherSCH08smallsubunitribosomalRNAgene,partialsaquenca;intern...Phoebe...NA460784109641096499%PhaebehuivoucherZF49smallsubunitribosomalRNAgene_partialsequence;intemaltrans..Phaebe...N4A60785108591085999%PhaabespZF55smallubunitribasoialRNAgen9,partialsaquence;internaltranscribeds..Phaabe...455541082910829100%PhoebeshearerivoucherWH01smallsubunitribosomalRNAgeng,.partialsequence;intema…BPhoebe...0789108131081399%PhoabeshearenivoucherZF45smallsubunitribosomalRNAgeng,partialsequance;intemal..Phogbe...N4A60789108111081199%CinnamomumcaudiferumvoucherFD24smallsubunitribosomalRNAgene,.patialsequenc...Cinnam.NA2745578105381053899%CinnarmomumparthenoxylonvoucherYB04smallsubunitibosomalRNAgene,partialseque.Cinnarn.NA714455105211052199%CinnamomumcamphoravoucherKZ05smallsubunitibosomalRNAgene.patialsecuence:..inam.campho.….13429104531045399%图B.1基于木材DNA条形码的序列一致性分析示例0.099.98%5941MT628616.10.099.66%5942MT628645.10.099.58%5943MT628648.10.099.56%5935MT628621.10.099.56%5934MT628643.10.098.71%5935MT628595.10.098.65%5835MT628624.10.098.47%5928MT628606.19GB/T43271—2023(资料性)木材DNA条形码序列的邻接树构建以分子进化遗传分析软件(MolecularEvolutionaryGeneticsAnalysis,MEGA¹)中邻接法(Neigh- bor-joiningmethod)为例。采用基于邻接法的分子系统发育树方法分析时，进入MEGA界面，从“文件 (File)”菜单中选择“打开文件(OpenaFile/Session)”,点击MEGA主窗口“系统发育(Phylogeny)”菜单选择“构建/测试邻接树(Construct/TestNeighbor-JoiningTree)”,选择“数据类型(DataType)”“选择遗传密码(SelectGeneticCode)”相应参数后进入“分析参数选择(AnalysisPreferences)”窗口，在“系统发育测试(PhylogenyTest)”菜单中选择“自展法(Bootstrapmethod)”,数值一般设置为1000;选择“模型/方法(Model/Method)”“替换类型(SubstitutionsType)”“位点进化速率(RatesamongSites)”和“空位处理(Gaps/MissingDataTreatment)”等下拉菜单相应参数后，点击“计算(Compute)”,确定聚类关系，完成分子系统发育树构建(见图C.1)。P.zhōn190_01P.zhōn191_02P.zhōnn8n_23.zhie/inan_25F.zliennan_27F.zhiennan_imH033832P.chekiangans3_01Pchekiangensis_0247P:hekiargansis_0:hekiargersis_D4F.hekiargmasis0E90Pshcarcr_02Pshcaror_03btPshGarr_0bt—P.bcu!mei_KY348512P.hGinsstha_MH384354P.neurantho_AH394255e.newanths_A-3943530.neUgnthe_A4394252pauIa_AH178403:1pawia_44Moauia_a29pauha_02janangnsis_KT54551850Miepicysp:05的Micpicp:07Micoisphpl02原icoiophpl_01原icoiophpl_0128氧unbcgi_03硫bunb6mi_64硫tounb₆gi_65试tounbsgiiMn178404CinnaromUn?csniphcrs_MH050970 图C.1基于木材DNA条形码的邻接树构建示例GB/T43271—2023(资料性)木材DNA条形码序列的最大似然树构建以MEGA软件中最大似然法(MaximumLikelihoodmethod)为例。采用基于最大似然法的分子系统发育树方法分析时，进入MEGA界面，从“文件(File)”菜单中选择“打开文件(OpenaFileSession)”,点击MEGA主窗口“系统发育(Phylogeny)”菜单选择“构建/测试最大似然树(Construct/TestMaximumLikelihoodTree)”,选择“数据类型(DataType)”“选择遗传密码(SelectGeneticCode)”相应参数后进入“分析参数选择(AnalysisPreferences)”窗口，在“系统发育测试(PhylogenyTest)”菜单中选择“自展法(Bootstrapmethod)”,数值一般设置为1000;选择“模型/方法(Model/Method)”“替换类型(SubstitutionsType)”“位点进化速率(RatesamongSites)”和“空位处理分子系统发育树构建(见图D.1)。Pzhennan_01F.zhannan_021P.2herwur_05P.2fcratr_07lpzhennaniH033832Pshg:srenr_(F.shgaar_02P.shearai_05P.bowxi_KYJ46512P.neuranira_43943