北柴胡特異DNA條形碼篩選及種質資源鑒定

王 馨，尹光耀，張志飛，陳 穎，劉珊瑚，滿金輝，石 玥，黃鈺瑩，張曉芹，王曉暉，魏勝利

? 藥材與資源 ?

北柴胡特異DNA條形碼篩選及種質資源鑒定

王 馨1，尹光耀1，張志飛1，陳 穎1，劉珊瑚1，滿金輝1，石 玥1，黃鈺瑩1，張曉芹1，王曉暉2*，魏勝利1*

1. 北京中醫藥大學中藥學院，北京 102488 2. 北京中醫藥大學中藥研究院，北京 100029

基于北柴胡葉綠體基因組篩選特異性DNA條形碼，并利用特異DNA條形碼鑒定不同產地北柴胡種質資源。利用Illumina HiSeq X Ten平臺對北柴胡進行葉綠體基因組測序，利用mVISTA軟件和核酸多樣性分析篩選特異性片段，基于特異性片段進一步分析不同產區北柴胡樣品單倍型。不同產地3份北柴胡葉綠體基因組全長為155 557～155 959 bp，均呈現典型的環狀四分體結構，均編碼131個基因。、、、可以作為潛在的北柴胡種質資源鑒定的特異性DNA條形碼。基于擴增效率，選擇和對來自4省7個產地177份北柴胡樣品進行序列分析，結果表明和分別有91和78個變異位點，分別鑒定到28、29個單倍型；兩段序列聯合分析形成40個單倍型（Hap1～Hap40），各單倍型的遺傳距離為0～0.022。6個產地擁有特有單倍型，可以將不同產地的北柴胡種質資源進行區分。利用比較葉綠體基因組學篩選的特異DNA條形碼和可以用于北柴胡種內種質資源鑒定，為后續鑒定北柴胡的產地來源、種質資源保護利用和育種等工作奠定基礎。

北柴胡；葉綠體基因組；序列對比；種質資源；DNA條形碼

柴胡具有疏散退熱、疏肝解郁、升舉陽氣之功效，是最常用的大宗藥材之一。《中國藥典》2020年版[1]規定柴胡為傘形科柴胡屬藥用植物北柴胡DC.或狹葉柴胡（紅柴胡）Willd.的干燥根。北柴胡是中藥柴胡的主要來源，現代藥理學作用研究表明其具有多種藥理活性，如解熱、鎮靜催眠、抗炎保肝、抗病毒、抗腫瘤和正向調節免疫等功能[2-8]。因此對北柴胡進行種質資源鑒別，對于北柴胡的質量控制、安全合理利用及開展相關的研究具有重要意義。

葉綠體是一種擁有自身遺傳物質的多功能細胞器，葉綠體基因組與核基因組和線粒體基因組相比，依賴母系遺傳，相對分子質量較小，基因組成和結構相對保守，葉綠體基因組包含的大量遺傳信息被廣泛應用于植物分子進化及系統發育的研究[9-10]。隨著高通量測序技術和生物信息學的快速發展，多種植物的葉綠體基因組已經解析。目前柴胡屬的北柴胡、紫花闊葉柴胡H. Wolff、竹葉柴胡Wall. ex DC.、紅柴胡、長莖柴胡Wall. ex DC.等[11-15]物種的葉綠體基因組已有研究報道，但是種內葉綠體比較基因組學相關研究鮮有報道。

DNA條形碼是有關物種鑒定的新技術，通過用標準化的、較短的DNA序列作為條形碼，實現對物種進行快速、準確的鑒定[16]，是傳統性狀鑒定和理化鑒定方法的有效補充。如能夠鑒別大黃不同基原及混偽品[17]；而、、、等條形碼為柴胡屬植物的分類鑒定與系統學研究提供了重要依據[18]。隨著葉綠體基因組測序技術和葉綠體比較基因組學的發展，通過比較葉綠體基因組序列差異獲得不同物種或不同基原的特異DNA條形碼也變得切實可行，如利用葉綠體基因組及比較基因組學方法獲得鑒定唐古特大黃、藥用大黃和掌葉大黃的特異性條形碼為、、、以及L[19]。基于葉綠體基因組測序和比較基因組學獲得黃芩的突變位點最多的3個特異DNA條形碼、、，利用3段葉綠體基因組DNA的聯合分析對黃芩的野生居群展開研究，共發現29個變異位點，形成50個單倍型[20]。但是目前基于北柴胡葉綠體比較基因組學的特異DNA條形碼研究還未見報道，因此本研究選取不同產地的3個北柴胡樣品進行葉綠體測序，進行比較基因組學研究，篩選特異性DNA條形碼，并利用篩選的特異DNA條形碼對4省7產地177份北柴胡樣品進行種質資源鑒定，為后續北柴胡種質篩選、質量控制及相關的研究奠定基礎。

1. 材料與儀器

1.1 材料

用于葉綠體全基因組測序的北柴胡樣品分別采自遼寧省朝陽市、陜西省寶雞市和河南省洛陽市（表1）；用于種質資源鑒定的7個產地177份樣品來源及單倍型分布如表2所示。所用樣品經北京中醫藥大學魏勝利教授鑒定為北柴胡DC.。

1.2 儀器

GL-88B型渦旋器（海門市其林貝爾儀器制造有限公司）；HH-S4A型電熱恒溫水浴鍋（北京科偉永興儀器有限公司）；TGL-16型冷凍高速離心機（湘儀離心機儀器有限公司）；NanoDrop one型超微量分光光度計（Thermo Fisher）；JY300E型電泳儀（北京君意東方電泳設備有限公司）；WD-9413B型凝膠成像儀（北京六一生物科技有限公司）。

表1 北柴胡葉綠體全基因組測序樣品來源

Table 1 Source of chloroplast whole genome sequencing samples from B. chinense

名稱采集地點經度（E）緯度（N）海拔/m 北柴胡1遼寧省朝陽市朝陽縣臺子鎮孫丈子村120°16'46.44"41°40'91.13" 753 北柴胡2陜西省寶雞市106°18'00.00"33°35'00.00" 813 北柴胡3河南省洛陽市嵩縣車村鎮龍池峰村112°43'44.68"34°63'04.10"1539

表2 7個產地177份北柴胡樣品來源及單倍型分布

Table 2 Source and haplotype distribution of 177 samples of B. chinense from seven regions

編號樣品來源單倍型類型（數量） HNLY1～HNLY28河南省洛陽市Hap1(11)、Hap2(11)、Hap3(3)、Hap4(1)、Hap5(1)、Hap6(1) HBQHD1～HBQHD23河北省秦皇島市Hap7(9)、Hap8(9)、Hap9(2)、Hap10(2)、Hap11(1) SXLFHM1～SXLFHM53山西省臨汾市侯馬市Hap2(3)、Hap5(4)、Hap6(1)、Hap12(1)、Hap13(1)、Hap14(8)、Hap15(1)、Hap16(6)、Hap17(7)、Hap18(3)、Hap19(5)、Hap20(2)、Hap21(2)、Hap22(1)、Hap23(1)、Hap24(1)、Hap25(1)、Hap26(1)、Hap27(3)、Hap28(1) SXXZ1～SXXZ9山西省新絳縣西莊村Hap17(2)、Hap24(1)、Hap29(1)、Hap30(1)、Hap31(2)、Hap32(1)、Hap33(1) SXGJ1～SXGJ4山西省新絳縣古交村Hap17(1)、Hap27(1)、Hap30(1)、Hap34(1) SXBJ1～SXBJ18陜西省寶雞市Hap1(1)、Hap2(1)、Hap5(1)、Hap12(2)、Hap19(1)、Hap30(2)、Hap31(1)、Hap32(2)、Hap34(1)、Hap35(3)、Hap36(1)、Hap37(1)、Hap38(1) SXSL1～SXSL42陜西省商洛市商州區代街村Hap1(22)、Hap2(16)、Hap5(1)、Hap14(1)、Hap39(1)、Hap40(1)

2 方法

2.1 基因組DNA提取與測序

利用FastPure?Plant DNA Isolation Mini Kit（諾唯贊）提取植物總DNA，1.0%的瓊脂糖凝膠電泳檢測DNA質量，超微量分光光度計檢測DNA濃度。檢測合格的基因組總DNA構建插入片段長度約350 bp的文庫，利用Illumina HiSeq X Ten平臺進行序列讀長為150 bp的雙端測序，并對原始序列進行4步處理，分別是去除質量值連續≤20的堿基數達到40%的reads、去除含N的堿基數目總和達到10%的reads、去除adapter污染和去除duplication污染，得到高質量待分析序列（clean reads）。

2.2 葉綠體基因組組裝、拼接和注釋

使用NOVOPlasty將待分析序列組裝成完整的葉綠體基因組，采用PGA軟件進行組裝結果進行注釋（默認參數）。采用BWA將待分析序列比對北柴胡參考基因組（NC_046774）的葉綠體基因組序列，利用在線軟件tRNAscan-SE（http://lowelab. ucsc.edu/tRNAscan-SE/)確定所有tRNA基因的邊界，并通過CLC Sequence Viewer 8人工檢查確保組裝無誤。利用Organellar Genome DRAW（https:// chlorobox.mpimp-golm.mpg.de/OGDraw. html）在線繪制葉綠體全基因組圖譜。

2.3 重復序列分析

利用REPuter（https://bibiserv.cebitec.uni-bielefeld. de/reputer）檢測葉綠體全基因組序列中的分散重復序列，參數設置為最小重復序列長度為30 bp，重復序列間的相似度＞90%。串聯重復序列利用Tandem repeats finder軟件（https://tandem.bu.edu/ trf/trf.html）進行檢測，參數選擇默認值。利用MISA軟件（http://pgrc.ipk- gatersleben.de/misa/）檢測簡單重復序列位點類型與數目，對于單、二、三、四、五和六核苷酸重復序列，重復次數的最小值分別設置為10、5、4、3、3和3，2個SSR之間最小距離設置為100 bp。

2.4 葉綠體基因組比較分析

利用mVISTA在線軟件對測序獲得3份葉綠體基因組進行全局比對分析，利用Dna SP 6軟件檢測3條葉綠體基因組的核苷酸多樣性（nucleotide polymorphism，P）。

2.5 高變區域擴增、遺傳距離分析與進化樹構建

為利用高變區基因對不同產區的種質資源進行分析，以高變區基因為模板設計引物，對收集到的7個產區177份樣品進行PCR擴增分析。PCR混合體系共50 μL，包括ddH2O 37 μL，10×Buffer 5 μL,dNTP 4 μL，引物-F 1 μL，引物-R 1 μL，DNA模板0.5 μL，0.1%BSA 1 μL，TaKaRa Taq 0.5 μL。擴增程序為93 ℃、3 min，50 ℃、2 min，30個循環（93 ℃、30 s，44.4或36.8 ℃、45 s，70 ℃、45 s），70 ℃、5 min，4 ℃保存，同時根據不同的引物設定不同的退火溫度（表3）。通過1.0%的瓊脂糖凝膠電泳和FastPure?Gel DNA Extraction Mini Kit（諾唯贊）對PCR產物進行純化，純化后的產物送北京六合華大公司進行雙向測序。

利用DNAMAN和Chromas軟件對測序結果進行核對和單倍型匯總。將匯總后的單倍型序列通過mafft軟件進行比對。比對后的序列通過MEGA X軟件計算單倍型間的遺傳距離，構建鄰接法系統進化樹，設置bootstrap重復值為1000。

表3 北柴胡2個基因引物信息

Table 3 Primer information of two genes of B. chinense

基因名引物名稱堿基序列退火溫度/℃ petN_psbMpetN_psbM-F5’-AGGTACTACGAATTGATTGAGGAATC-3’44.4 petN_psbM-R5’-TTAATTTTAATTAATGTTTACTTCTG-3’ ndhF_rpl32ndhF_rpl32-F5’-AAAAAAAGTTTTTTTTCTTAATTAAT-3’36.8 ndhF_rpl32-R5’-TTTAAAATGAAGAGGTTACTCGTTG-3’

3 結果與分析

3.1 葉綠體全基因組測序、序列組裝與拼接注釋

將北柴胡1、北柴胡2和北柴胡3的測序結果過濾去除低質量序列和接頭序列等雜質,分別得到30 038 384條（4.50 Gb）、29 684 332條（4.45 Gb）和27 841 256條（4.18 Gb）clean reads。將3份北柴胡葉綠體全基因組組裝拼接后均得到完整的環狀四分體結構，序列總長度分別為155 557、155 959和155 959 bp。環狀四分體結構由1個大單拷貝區（large single copy area，LSC）、1個小單拷貝區（small single copy region，SSC）和2個反向重復區（inverted repeat regions，IRA和IRB）組成，其中LSC長度分別為85 431、85 831、85 831 bp，SSC長度均為17 546 bp，IRs長度分別為26 290、26 291和26 291 bp（圖1）。北柴胡2和北柴胡3的葉綠體全基因組序列及LSC、SSC和IRs的GC含量均相同，分別為37.68%、35.80%、31.48%和42.81%。北柴胡1的葉綠體全基因組序列及LSC、SSC和IRs的GC含量依次為37.69%、35.83%、31.43%和42.81%（表4）。

外圈基因逆時針轉錄，內圈基因順時針轉錄；彩色條表示不同的功能組；內圈較深的灰色區域表示GC含量，而較淺的灰色對應于基因組的AT含量。LSC-大單拷貝區 SSC-小單拷貝區域 IRA和IRB-反向重復區

表4 3份北柴胡的葉綠體基因組特征概述

Table 4 Summary of chloroplast genome characteristics of three B. chinense

樣本全基因組LSCSSCIRsCDS 長度/bp(G＋C)/%長度/bp(G＋C)/%長度/bp(G＋C)/%長度/bp(G＋C)/%長度/bp(G＋C)/% 北柴胡1155 55737.6985 43135.8317 54631.4352 58042.8177 79338.20 北柴胡2155 95937.6885 83135.8017 54631.4852 58242.8177 79338.21 北柴胡3155 95937.6885 83135.8017 54631.4852 58242.8177 79338.21

3個不同產地的北柴胡葉綠體全基因組序列都注釋了131個基因（表5），包括86個蛋白質編碼基因、37個tRNA基因和8個rRNA基因。其中，SSC區含有11個蛋白質編碼基因（、、、、、、、、、和）和1個tRNA（）；有17個基因分別在2個IR區出現1次，包含6個蛋白質編碼基因（、、、、和），7個tRNA（、、、、、和）和4個rRNA（、、和）；LSC區含有61個蛋白質編碼基因和22個tRNA。此外，蛋白質編碼基因和分別橫跨LSC/IRb邊界和SSC/IRA邊界。含1個內含子的基因有16個（、、、、、、、、、、、、、、和），含2個內含子的基因有2個（和）。

表5 3份北柴胡葉綠體基因組的基因組成

Table 5 Gene composition of chloroplast genomes of three B. chinense

類別基因群基因名稱 光合系統基因光系統IpsaA、psaB、psaC、psaI、psaJ 光系統IIpsbA、psbB、psbC、psbD、psbE、psbF、psbH、psbI、psbJ、psbK、psbL、psbM、psbN、psbT、psbZ NADH脫氫酶ndhA*、ndhB*(2)、ndhC、ndhD、ndhE、ndhF、ndhG、ndhH、ndhI、ndhJ、ndhK 細胞色素b/f復合體petA、petB*、petD*、petG、petL、petN ATP合成酶atpA、atpB、atpE、atpF*、atpH、atpI 二磷酸核酮糖羧化酶大亞基rbcL 遺傳系統基因核糖體蛋白大亞基rpl14、rpl16*、rpl2*(2)、rpl20、rpl22、rpl23(2)、rpl32、rpl33、rpl36 核糖體蛋白小亞基rps11、rps12、rps12*(2)、rps14、rps15、rps16*、rps18、rps19、rps2、rps3、rps4、rps7(2)、rps8 RNA聚合酶rpoA、rpoB、rpoC1*、rpoC2 核糖體RNArrn16(2)、rrn23(2)、rrn4.5(2)、rrn5(2) 轉運RNAtrnA-UGC*(2)、trnC-GCA、trnD-GUC、trnE-UUC、trnF-GAA、trnG-GCC、trnG-UCC*、trnH-GUG、trnI-CAU(2)、trnI-GAU*(2)、trnK-UUU*、trnL-CAA(2)、trnL-UAA*、trnL-UAG、trnM-CAU、trnN-GUU(2)、trnP-UGG、trnQ-UUG、trnR-ACG(2)、trnR-UCU、trnS-GCU、trnS-GGA、trnS-UGA、trnT-GGU、trnT-UGU、trnV-GAC(2)、trnV-UAC*、trnW-CCA、trnY-GUA、trnfM-CAU 其他基因成熟酶matK 蛋白酶clpP** 胞膜蛋白cemA 乙酰輔酶A羧化酶accD 細胞色素C合成酶ccsA 翻譯起始因子infA 未知功能基因假定葉綠體開放性閱讀框ycf1、ycf2(2)、ycf3**、ycf4

*-含有1個內含子的基因**-含有2個內含子的基因 (2)-有2個拷貝的基因

*-genes with one intron**-genes with two introns (2)-genes with two copies

3.2 重復序列分析

不同產地的北柴胡樣品葉綠體全基因組序列中北柴胡1發現66條分散重復序列，北柴胡2和北柴胡3均發現75條分散重復序列。其中北柴胡2和北柴胡3包括29條正向重復（forward repeats）、12條反向重復（reverse repeats）、29條回文重復（palindromic repeats）和5條互補重復（complement repeats）；北柴胡1包括30條正向重復、4條反向重復、28條回文重復和4條互補重復序列（圖2）。同時檢測到3份北柴胡樣品中長度為9～49 bp、重復次數為2～4次的串聯重復序列分別為37、38和38條。對不同產地北柴胡樣品的葉綠體全基因組簡單重復序列（SSR）的數量、類型及其在葉綠體基因組中的分布情況分析顯示，SSRs主要分布于葉綠體基因組的LSC區（71.83%～72.31%）和非編碼的基因間區序列中。編碼基因序列總長度占葉綠體全基因組的49.88%～50.00%，而編碼基因序列中分布的SSRs數量僅占其總數的15.38%～15.49%（表6、7）。

圖2 3份北柴胡葉綠體基因組中重復序列的數量

3.3 比較基因組學分析

為檢測3份北柴胡葉綠體基因組的種內變異情況，以注釋過的北柴胡樣本葉綠體基因組（NC_046774）為參考，利用mVISTA在線軟件將獲得的其余3個葉綠體基因組進行全局對比（圖3）。結果顯示，葉綠體基因組序列中蛋白編碼區域變異低于非編碼區變異，且IRs區變異明顯小于其他LSC區和SSC區。此外，絕大多數基因的相似度在90%以上，rRNA基因（、、、）高度保守沒有變異。基因區、、、和基因間區、、、、、等基因有不同程度的變異。同時，對3份北柴胡樣品的P進行分析顯示種內變異度較小，P值范圍為0～0.017 78，LSC、SSC和IRs區的平均P值分別為0.001 32、0.001 63和0.000 24，與mVISTA分析結果一致（圖4）。基因區（0.008 89）和基因間區（0.006 67）、0.017 78）、（0.006 67）的P值較高。綜上可知，北柴胡種內基因的保守程度很高，基因間區變異程度高于基因區，結合mVISTA和P結果，顯示基因區和基因間區、、可作為潛在的分子高變區域。

表6 SSRs在3份北柴胡葉綠體基因組中的分布

Table 6 Distribution of SSRs in chloroplast genomes of three B. chinense

樣本SSRs總量SSRs數量（占比/%） LSCIRsSSCCDS 北柴胡16547（72.31）17（26.15）1（1.54）10（15.38） 北柴胡27151（71.83）16（22.54）4（5.63）11（15.49） 北柴胡37151（71.83）16（22.54）4（5.63）11（15.49）

表7 3份北柴胡葉綠體基因組中SSRs種類及其數量

Table 7 Types and numbers of SSRs in chloroplast genomes of three B. chinense

樣本SSRs總量核苷酸類型占比（數量） 單核苷酸/%二核苷酸/%三核苷酸/%四核苷酸/%五核苷酸/%六核苷酸/% 北柴胡16540（61.54）12（18.46）6（9.23）4（6.15）2（3.08）1（1.54） 北柴胡27149（69.01）10（14.08）7（9.86）5（7.04）// 北柴胡37149（69.01）10（14.08）7（9.86）5（7.04）//

3.4 不同產地北柴胡樣品的petN_psbM和ndhF_rpl32的擴增

根據（643 bp）、（656 bp）、（1102 bp）、（1220 bp）序列設計引物，進行PCR擴增，發現（643 bp）、（656 bp）擴增效率不高，、的擴增效率為100%，因此本研究選擇、為北柴胡種內鑒定的特異DNA條形碼，進行后續研究。以7個產地177份北柴胡樣品總DNA為模板對篩選的2個葉綠體基因組高變區進行PCR擴增，并對測序結果進行分析。利用設計的引物進行PCR擴增，全部樣品擴增到1000～2000 bp（圖5），對PCR產物純化后測序，與北柴胡1的序列對比可判斷其為序列；利用設計的引物進行PCR擴增，全部樣品擴增到1000～2000 bp（圖6），對PCR產物純化后測序，與北柴胡1的序列對比可判斷其為序列。

圖4 3份北柴胡葉綠體基因組核苷酸多樣性

3.5 不同產地樣品DNA條形碼分析

基因共檢測到91個變異位點，共形成28個單倍型（圖7）。其中位于249、396、576、612、837和926 bp的6個位點突變類型為插入；位于932～976 bp的45個位點突變類型為缺失，同時942～944 bp和950 bp伴有單點突變；其他40個位點為單點突變，分別位于219、228、303、327、331、362、401、403、428、445、485、536、578、581、607、620、621、626、628～631 bp、633、636、639、641～644、646、652、682、708、716、749、869、873、883、928、1013 bp。

M-Marker 1-HNLY10 (Hap3) 2-HBQHD1 (Hap7) 3-SXLFHM19 (Hap20) 4-SXXZ2 (Hap29) 5-SXGJ1 (Hap34) 6-SXBJ1 (Hap35) 7-SXSL7 (Hap39) 8-陰性對照，圖6同

圖6 北柴胡樣品ndhF_rpl32序列PCR擴增

基因共檢測到78個變異位點，共形成29個單倍型（圖8）。其中位于255、256、280、331、749、794、803、807、950、1001bp的10個位點突變類型為插入，此外，749 bp位點也存在缺失突變和單點突變；位于258、750、775～780、819～827、853～856 bp的21個位點突變類型為缺失，同時853 bp伴有單點突變；其他47個位點為單點突變，分別位于132、134、167、178、179、184、207、223、224、262、279、286、317、318、343、353、356、423、449、462、472、478、482、518、570、574、577、594、647、656、713、741、742、790、796、809、839、874、876、877、886、960、1 011、1014、1015、1022和1044 bp。

對2基因聯合分析共形成了40個單倍型（表2），占比最多分布最廣的單倍型為Hap1，占全部樣品的19.21%，其次是Hap2，占比為17.51%（圖9），而特異的單倍型可作為該產地特有的種質資源加以擴繁，也可以作為鑒定北柴胡產地來源的DNA分子標簽。河南省洛陽市的DNA分子標簽為Hap3和Hap4，河北秦皇島市的DNA分子標簽為Hap7、Hap8、Hap9、Hap10和Hap11，山西省臨汾市侯馬市的DNA分子標簽為Hap13、Hap15、Hap16、Hap18、Hap20、Hap21、Hap22、Hap23、Hap25、Hap26和Hap28，山西省新絳縣西莊村的DNA分子標簽為Hap29和Hap33，陜西省寶雞市的DNA分子標簽為Hap35、Hap36、Hap37和Hap38，陜西省商洛市商州區代街村的DNA分子標簽為Hap39和Hap40。

*-核酸與Hap1相同 “-”-核酸缺失，圖8同

*-nucleic acid is the same as Hap1 “-”-nucleic acid deletion, same as fig.8

圖7基因的單倍型變異位點

Fig. 7 Haplotypes variation outliers ofgene

3.6 遺傳距離分析

使用Mega X對2個基因聯合分析的單倍型進行遺傳距離分析，結果顯示40個單倍型的遺傳距離為0～0.022，平均遺傳距離為0.010，表明北柴胡樣品間親緣關系差別不大。最大遺傳距離0.022存在于Hap3和Hap20之間，親緣關系最遠；最小遺傳距離為0.000存在于多個單倍型之間，分別在Hap6、Hap16、Hap17和Hap33之間，Hap14、Hap23和Hap24之間，Hap26和Hap27之間，Hap35和Hap38之間，親緣關系最近。占比較多的Hap1和Hap2之間遺傳距離為0.017，親緣關系較遠。結合表2中每個產地樣品的單倍型發現，Hap35和Hap38為陜西省寶雞市的特異單倍型，Hap16為山西省臨汾市侯馬市特異單倍型，Hap33為山西省新絳縣西莊村特異單倍型。此外，陜西省寶雞市的另2個特異單倍型Hap36和Hap37之間遺傳距離為0.001，親緣關系較近。山西省臨汾市侯馬市的另一特異單倍型Hap15和山西省新絳縣西莊村的另一特異單倍型Hap29之間遺傳距離為0.003，親緣關系較近。河北省秦皇島市的2個特異單倍型Hap7和Hap11之間遺傳距離為0.001，親緣關系較近。陜西省商洛市商州區代街村的特異單倍型Hap40和河南省洛陽市的特異單倍型Hap4之間遺傳距離為0.002，親緣關系較近。

對表2中每個產地的特異單倍型之間遺傳距離分析發現，河南省洛陽市的特異單倍型Hap3與河北省秦皇島市的特異單倍型Hap7之間遺傳距離為0.021；Hap3與山西省臨汾市侯馬市特異單倍型Hap20之間遺傳距離為0.022，親緣關系最遠。河南省洛陽市的另一特異單倍型Hap4與山西省新絳縣西莊村的特異單倍型Hap29之間遺傳距離為0.018；Hap7與Hap29之間的遺傳距離為0.019；Hap20與Hap29之間的遺傳距離為0.020，親緣關系較遠。Hap4與陜西省寶雞市的特異單倍型Hap35之間遺傳距離為0.019；Hap7與Hap35之間的遺傳距離為0.020；Hap20與Hap35之間的遺傳距離為0.021，親緣關系較遠。Hap4與陜西省商洛市商州區代街村的特異單倍型Hap39之間遺傳距離為0.019；Hap7與Hap39之間遺傳距離為0.020；Hap20與Hap39之間遺傳距離為0.021；Hap29與陜西省商洛市商州區代街村的另一特異單倍型Hap40之間遺傳距離為0.017；Hap35與Hap40之間遺傳距離為0.017，親緣關系較遠。這些結果表明每個產地的特異單倍型之間的遺傳距離較遠。

圖9 2個基因聯合分析的單倍型占比圖

3.7 系統進化分析

將2個基因聯合分析后的單倍型序列對比后構建NJ Tree，從而進一步分析各單倍型的親緣關系。如圖10所示，Hap1、Hap3、Hap14、Hap18、Hap19、Hap23、Hap24、Hap28、Hap35和Hap38聚在同一支上，其他單倍型聚為一支，同一支的單倍型之間親緣關系較近。占比較多的Hap1和Hap2分別屬于2個不同的單系分支，親緣關系較遠。

結合表2分析各產地的特異單倍型發現，部分單倍型之間親緣關系較近，與遺傳距離分析結果一致。陜西省寶雞市的特異單倍型Hap35和Hap38聚為同一支，支持率為62%；陜西省寶雞市的另2個特異單倍型Hap36和Hap37聚為同一支，支持率為98%；河北省秦皇島市的2個特異單倍型Hap7和Hap11聚為同一支，支持率為96%；陜西省商洛市商州區代街村的特異單倍型Hap40和河南省洛陽市的特異單倍型Hap4聚為同一支，支持率為92%；山西省臨汾市侯馬市的特異單倍型Hap16與山西省新絳縣西莊村的特異單倍型Hap33聚為同一支，支持率為23%；山西省臨汾市侯馬市的另一特異單倍型Hap15和山西省新絳縣西莊村的另一特異單倍型Hap29聚為同一支，支持率為62%。

結合表2中每個產地的特異單倍型發現，與遺傳距離分析結果一致，遺傳距離較大的特異單倍型均不聚在同一支，親緣關系較遠。

圖10 2個基因聯合分析的單倍型NJ Tree

4 討論

據報道北柴胡葉綠體基因組為典型的環狀四分體結構，總長度為155 458～155 869 bp，GC含量為37.68%～37.70%，IR長度為52 596～52 620 bp，LSC長度為85 343～85 772 bp，SSC長度為17 495～17 505 bp；共包含129～133個基因，其中蛋白質編碼基因為19～84個，tRNA基因為30～37個，rRNA基因為4-8個[11, 21-22]。北柴胡及其同屬植物的分散重復序列為33～49條，串聯重復序列為22～30條，簡單重復序列為57～72條[21-22]。本研究所測的3個北柴胡葉綠體全基因組大小和結構與上述北柴胡研究結果基本相符。

目前已經有多篇文獻報道[19, 23-25]利用葉綠體測序和比較基因組進行篩選高變區作為特異DNA條形碼用于物種的鑒定，但是用篩選的特異DNA條形碼驗證的比較少。如通過比較葉綠體基因組學研究有學者推薦、、、、等葉綠體基因組高變區作為鑒定石豆蘭屬藥用植物的特異DNA條形碼[23]；、、、、是通過葉綠體基因組學分析篩選的龍膽屬植物鑒定的高變片段[24]。目前也有少量研究對葉綠體基因組高變區用于鑒定物種進行了驗證，如基于葉綠體基因組篩選的秦艽的特異性DNA片段和2個高變區，能夠用于長梗秦艽、全萼秦艽的鑒別[25]；高變區、、、與能作為特異DNA條形碼鑒定大黃藥材3個基原[19]。但是也有相關報道證明利用葉綠體基因組測序獲得特異DNA條形碼后擴增效果不理想，如豆蔻屬的特異性DNA片段擴增效果不理想，不適合豆蔻屬的鑒定[26]。綜合文獻研究表明[25-26]，在篩選高變區域后，需要經過實際擴增效率驗證。近期已有學者通過比較基因組學分析柴胡屬植物葉綠體基因組高變區為、、、、、為柴胡屬種間高變區，但是沒有相應的實驗驗證[22]。本研究基于mVISTA和Pi結果分析北柴胡葉綠體基因組的高變區，結果表明、、、是北柴胡種內的潛在DNA條形碼，結合以前的文獻[22]報道結果表明、既可以作為北柴胡種內特異性DNA片段也可以作為種間鑒定的DNA條形碼。基于擴增效率，選擇、用于后續北柴胡種質資源的鑒定。

DNA條形碼技術在藥用植物物種鑒定及種內資源鑒定方面起到重要的作用，具有簡便、高效、客觀等特點。目前DNA條形碼技術也在柴胡屬藥用植物鑒定方面取得顯著進展，如趙晴等[27]利用ITS序列鑒定了北柴胡、紅柴胡、三島柴胡的種子；曾偉萍[28]考察了DNA條形碼候選序列、、、在柴胡屬藥用植物的鑒定能力，發現序列和序列可作為柴胡屬藥用植物鑒定的候選序列。但是目前關于北柴胡種內資源鑒別的研究比較少。本研究基于葉綠體高變區篩選的2條DNA條形碼和對北柴胡的遺傳多樣性進行分析，發現了大量的變異位點和較全面的遺傳信息。其中6個產地的特異單倍型作為該產地特有的種質資源，為后續鑒別不同產地北柴胡樣品奠定基礎；各單倍型之間遺傳距離較小，北柴胡種內親緣關系較近，但主流單倍型Hap1和Hap2之間親緣關系相對較遠，部分產地特異單倍型之間親緣關系較遠。因此本研究基于比較葉綠體基因組學篩選的北柴胡鑒定的特異性片段和可以有效鑒定北柴胡種質資源，為北柴胡種質資源保護和后續育種工作奠定基礎。

利益沖突 所有作者均聲明不存在利益沖突

[1] 中國藥典[S]. 一部. 2020: 293.

[2] 孫世君. 柴胡的藥理學分析以及臨床應用 [J]. 中國醫藥指南, 2010, 8(29): 210-211.

[3] 金順姬. 柴胡的藥理作用及臨床應用 [J]. 現代醫藥衛生, 2009, 25(7): 1074-1075.

[4] 王占一, 南極星. 北柴胡對對乙酰氨基酚所致小鼠急性肝損傷的保護作用 [J]. 中國藥師, 2008, 11(7): 747-749.

[5] 胡繼鷹, 許湘, 潘克英, 等. 保康北柴胡解熱抗炎作用的藥效學研究 [J]. 中醫藥學刊, 2005, 23(4): 631-632.

[6] 劉曉斌, 高燕, 劉永仙, 等. 北柴胡提取組分對小鼠淋巴細胞活性的影響 [J]. 細胞與分子免疫學雜志, 2002, 18(6): 600-601.

[7] 馮煦, 王鳴, 趙友誼, 等. 北柴胡莖葉總黃酮抗流感病毒的作用 [J]. 植物資源與環境學報, 2002, 11(4): 15-18.

[8] 安玉明, 接傳胤, 祝世偉. 柴胡的研究進展 [J]. 人參研究, 2001, 13(1): 11-13.

[9] 楊小英, 劉暢, 曾憲法, 等. 八角楓及其亞種葉綠體基因組序列結構及系統發育分析 [J/OL]. 藥學學報, 2022: [2022-07-18]. https://kns.cnki.net/kcms/ detail/11. 2163.R.20220716.1931.002.html.

[10] 李卓蔚, 邱遷, 郎佳琪, 等. 尖刀唇石斛和翅梗石斛葉綠體全基因組分析 [J]. 中草藥, 2022, 53(16): 5159-5169.

[11] Zhang F, Zhao Z Y, Yuan Q J,. The complete chloroplast genome sequence ofDC. (Apiaceae) [J]., 2019, 4(2): 3665-3666.

[12] Wu Y, Zhang T Z, Qiu D Y,. Complete plastid genome of, an endemic herb plant in Western China [J]., 2018, 10(4): 635-637.

[13] Deng X D, Liu H Z, Yang Y H,. The first complete chloroplast genome sequence of the medicinal plant(Apiaceae) [J]., 2020, 5(2): 1836-1838.

[14] Zhang F, Yang Z, Wang Z,. The complete chloroplast genome sequence ofWilld. (Apiaceae) [J]., 2020, 5(2): 1998-1999.

[15] Yang L C, Xiong F, Xiao Y M,. The complete chloroplast genome ofvar., an annual herb endemic to China [J]., 2020, 5(1): 899-901.

[16] Hebert P D N, Cywinska A, Ball S L,. Biological identifications through DNA barcodes [J]., 2003, 270(1512): 313-321.

[17] 張曉芹, 劉春生, 閆興麗, 等. 多基原藥材大黃葉綠體基因序列分析及鑒定研究 [J]. 藥學學報, 2013, 48(11): 1722-1728.

[18] 袁伯川. 柴胡屬藥用植物的分子鑒定及柴胡藥材的質量考察 [D]. 北京: 北京中醫藥大學, 2017.

[19] 李冉郡, 武立偉, 辛天怡, 等. 大黃藥材基原物種葉綠體基因組分析與特異DNA條形碼開發 [J]. 藥學學報, 2022, 57(5): 1495-1505.

[20] 姜丹. 黃岑道地性的遺傳和化學物質基礎研究 [D]. 北京: 北京中醫藥大學, 2018.

[21] 張明英, 張雨曲, 李依民, 等. 北柴胡、紫花闊葉柴胡葉綠體全基因組解析及柴胡屬藥用植物葉綠體基因組比較與系統發育分析 [J]. 藥學學報, 2021, 56(2): 618-629.

[22] Li J, Xie D F, Guo X L,. Comparative analysis of the complete plastid genome of fivespecies and new insights into DNA barcoding and phylogenetic relationship [J]., 2020, 9(4): 543.

[23] 楊嘉鵬, 朱紫樂, 范雅娟, 等. 三種石豆蘭屬藥用植物的葉綠體基因組比較分析及其在物種鑒定中的意義 [J]. 藥學學報, 2020, 55(11): 2736-2745.

[24] 胡海粟, 張德全. 幾種滇產龍膽屬藥用植物的DNA超級條形碼研究 [J]. 中國中藥雜志, 2021, 46(20): 5260-5269.

[25] 董博然, 趙志禮, 倪梁紅, 等. 基于葉綠體基因組的長梗秦艽、全萼秦艽物種DNA分子標記與鑒定 [J]. 藥學學報, 2021, 56(9): 2584-2591.

[26] 崔英賢. 藥食兩用藥材砂仁、枸杞、山楂和姜基原植物葉綠體基因組結構解析 [D]. 北京: 北京協和醫學院, 2020.

[27] 趙晴, 謝紅波, 央拉, 等. 基于DNA條形碼技術的北柴胡種子分子鑒定 [J]. 中國實驗方劑學雜志, 2020, 26(14): 182-189.

[28] 曾偉萍. 柴胡屬藥用植物的條形碼研究[D]. 廣州: 南方醫科大學, 2016.

Screening of specific DNA barcoding and identification of germplasm resources of

WANG Xin1, YIN Guang-yao1, ZHANG Zhi-fei1, CHEN Ying1, LIU Shan-hu1, MAN Jin-hui1, SHI Yue1, HUANG Yu-ying1, ZHANG Xiao-qin1, WANG Xiao-hui2, WEI Sheng-li1

1. School of ChineseMateria Medica, Beijing University of Chinese Medicine, Beijing 102488, China 2. Institute of Chinese Medicine, Beijing University of Chinese Medicine, Beijing 100029, China

To screen the specific DNA barcodes based on Beichaihu () chloroplast genome, so as to identify the plasm resource offrom different producing areas.The Illumina HiSeq X Ten platform was used to sequence the chloroplast genome of three samples of. Specific DNA barcodes was analyzed with mVISTA soft and nucleotide polymorphism analysis. The haplotype offrom different producing areas was identified with specific DNA barcodes.The total length of chloroplast genomes of three samples from different habitats was 155 557—155 959 bp, all of which had typical circular tetrad structure and coded 131 genes.,,andcan be used as specific DNA barcodes for identification of potential germplasm resources of. Based on amplification efficiency,andwere selected to analyse the sequence of 177 samples offrom seven places in four provinces. The results showed thatandhad 91 and 78 mutation loci, and 28 and 29 haplotypes were identified; The two sequences were analyzed together to form 40 haplotypes (Hap1—Hap40), and the genetic distance of each haplotype was 0—0.022. Six producing areas had unique haplotypes, which can distinguish the germplasm resources offrom different producing areas.The specific DNA barcodesand, identified bychloroplast genome analysis,could be used to demonstrate the genotypes of, which provides values for the identification ofproducing area, the protection and utilization of germplasm resources and breeding work.

DC.; chloroplast genome; sequence comparison; germplasm resources; DNA barcoding

R286.12

A

0253 - 2670(2023)17 - 5703 - 13

10.7501/j.issn.0253-2670.2023.17.023

2023-02-03

北京市科學技術委員會基金項目（Z201100005420005）；北柴胡精準藥材批次分子防偽技術研究項目（2020071720419）

王 馨（1999—），碩士，主要從事中藥資源與分子生藥學研究。Tel: 18813062869 E-mail: nz18813062869@163.com

王曉暉，副研究員，主要從事中藥活性成分生物合成調控研究。Tel: 18810722975 E-mail: wangxhui2014@163.com

魏勝利，教授，博士生導師，主要從事中藥資源研究。Tel: 13683336930 E-mail: wsl7491@126.com

[責任編輯 時圣明]