王明強， 湯詩杰， 方 彥， 薛佳宇

〔1. 江蘇省中國科學(xué)院植物研究所(南京中山植物園) 植物多樣性及系統(tǒng)演化研究中心， 江蘇 南京 210014；2. 南京森林警察學(xué)院， 江蘇 南京 210023〕

茅山地區(qū)南蒼術(shù)的特異性分子標記研究

王明強1，①， 湯詩杰1，①， 方 彥2， 薛佳宇1，②

〔1. 江蘇省中國科學(xué)院植物研究所(南京中山植物園) 植物多樣性及系統(tǒng)演化研究中心， 江蘇 南京 210014；2. 南京森林警察學(xué)院， 江蘇 南京 210023〕

以中國野生南蒼術(shù)〔Atractylodeslancea(Thunb.) DC.〕全部分布區(qū)9個居群(包括江蘇句容的茅山居群、射烏山居群和小九華山居群，安徽金寨的燕子河居群和天堂寨居群，河南嵩縣的五馬寺居群和紙坊居群，湖北的?？岛笃壕尤汉图t安七里坪居群)為研究對象，進行rbcL、matK、psbA-trnH、ITS和nadF片段序列比對及SRAP指紋圖譜比較來篩選江蘇句容茅山居群南蒼術(shù)的特異性分子標記。結(jié)果表明：通過5個DNA片段的序列比對，江蘇句容茅山居群南蒼術(shù)在nadF片段的第581位點處堿基為T，其余8個居群為A，可以將江蘇句容茅山居群與其余8個居群區(qū)分開。利用SRAP隨機引物篩選出1條江蘇句容茅山居群南蒼術(shù)特有的長度為323 bp的條帶，根據(jù)該條帶序列設(shè)計1對SCAR引物，正向引物序列為5′-ATCTTTGCTATTAACTATTAAGCTATTCTTCGCTATTCA-3′，反向引物序列為5′-CGATGTTGAAATTGCTGTCTCTGCTTTCTATTC-3′，該引物可在江蘇句容茅山居群所有南蒼術(shù)單株中擴增得到長度為268 bp的特異性條帶，而在其他居群中均不能擴增出該條帶。本研究篩選出的2種特異性標記可以為鑒定南蒼術(shù)的道地性提供可靠的技術(shù)支持。

南蒼術(shù)； 道地性； 茅山地區(qū)；nadF片段； SRAP-SCAR

南蒼術(shù)〔Atractylodeslancea(Thunb.) DC.〕，又名茅蒼術(shù)，為菊科(Asteraceae)蒼術(shù)屬(AtractylodesDC.)植物，主要分布于江蘇、安徽、湖北和河南等狹小區(qū)域，生于野生山坡草地、林下、灌叢及巖石縫隙中[1]。蒼術(shù)是中國傳統(tǒng)中藥材[2]，江蘇茅山地區(qū)的南蒼術(shù)質(zhì)量尤佳，美譽自古有之?！侗静萁?jīng)集注》[3]中記錄蒼術(shù)“今處處有之，以蔣山、白山、茅山者為勝”，即今江蘇的南京、句容等地，亦可知茅山南蒼術(shù)早已是道地藥材。

彭華勝等[4]觀察到，南蒼術(shù)葉及根狀莖的形態(tài)隨海拔升高呈現(xiàn)連續(xù)的變化，且在不同生長期和不同生境，同一居群不同個體間也有較大差異，不同居群的個體從形態(tài)上很難區(qū)分。在揮發(fā)油總含量上，茅山地區(qū)南蒼術(shù)較其他地區(qū)的南蒼術(shù)低[5]，且揮發(fā)油中的蒼術(shù)素和蒼術(shù)酮含量高，不含或含較少的茅術(shù)醇和β-桉葉醇[6-7]，但揮發(fā)油含量的測定周期較長、檢測成本較高，且個體間差異大[8]，因此，化學(xué)成分含量無法成為一種有效而穩(wěn)定的鑒別依據(jù)。

與形態(tài)和化學(xué)成分相比，植物的基因組DNA是穩(wěn)定的遺傳物質(zhì)，利用DNA分子標記鑒別道地藥材更為快速、準確。已有研究結(jié)果顯示，利用DNA序列差異可以鑒別道地藥材。余永邦等[9]對14個地區(qū)太子參〔Pseudostellariaheterophylla(Miq.) Pax〕的ITS序列進行了擴增和測序，結(jié)果表明ITS序列可以鑒別江蘇產(chǎn)區(qū)及其他非道地產(chǎn)區(qū)的太子參；張宏意等[10]也用ITS序列鑒別四川、廣東、湖北等地的何首烏〔Fallopiamultiflora(Thunb.) Harald.〕；孔德政等[11]用psbA-trnH序列分析牛膝(AchyranthesbidentataBlume)道地產(chǎn)區(qū)河南8個區(qū)域的野生居群，可以將新鄉(xiāng)九蓮山、南陽老界嶺、信陽靈山和南陽淮源的牛膝鑒別開。南蒼術(shù)也有類似的研究，如利用trnL-F和ITS[12-13]等序列可以將南蒼術(shù)與蒼術(shù)屬其他植物區(qū)分開，但尚無利用DNA片段序列差異鑒別道地產(chǎn)區(qū)與非道地產(chǎn)區(qū)南蒼術(shù)的相關(guān)研究報道。

此外，利用隨機引物擴增圖譜差異也可以鑒定道地藥材，如利用RAPD、SSR和ISSR圖譜等[14-17]可以區(qū)分不同來源的球花石斛(DendrobiumthyrsiflorumRchb. f.)、當歸〔Angelicasinensis(Oliv.) Diels〕、鐵皮石斛(DendrobiumofficinaleKimura et Migo)和貝母(Fritillariaspp.)等藥材，將隨機引物擴增出的特有條帶測序后轉(zhuǎn)化成特異性SCAR標記，能顯著提高隨機引物標記的穩(wěn)定性及可重復(fù)性，在川牛膝(CyathulaofficinalisK. C. Kuan)和寧夏枸杞(LyciumbarbarumLinn.)等[18-19]藥材的鑒別研究中均有報道。在針對南蒼術(shù)的研究中，任冰如等[20]利用擴增出的RAPD指紋圖譜差異，將南蒼術(shù)分為南蒼術(shù)類群和大別山類群；郭蘭萍等[21]利用RAPD探討不同產(chǎn)地南蒼術(shù)遺傳多樣性時，認為茅山南蒼術(shù)單獨為一類。但目前尚無利用隨機引物的指紋圖譜來開發(fā)用于南蒼術(shù)道地性鑒別的特異性SCAR分子標記。

本文利用南蒼術(shù)不同居群基因組序列差異及特異性分子標記，建立南蒼術(shù)道地性鑒別的分子鑒別體系，以期為道地藥材的種質(zhì)保護提供簡便而穩(wěn)定的檢測方法。

1 材料和方法

1.1 材料

在中國野生南蒼術(shù)分布區(qū)江蘇、安徽、河南及湖北采集9個居群100個單株(表1)，每株取2～3枚幼嫩且無病蟲害葉片，清洗后保存于-80 ℃，用于總DNA的提取。

實驗材料經(jīng)江蘇省中國科學(xué)院植物研究所杭悅宇研究員鑒定，憑證標本存放于江蘇省中國科學(xué)院植物研究所標本館(NAS)。

1.2 方法

1.2.1 樣品總DNA提取 每個樣品取0.1 g葉片，用EasyPure Plant Genomic DNA Kit試劑盒(北京全式金生物技術(shù)有限公司)提取總DNA，溶于200 μL滅菌超純水中，-20 ℃保存、備用。

表1 南蒼術(shù)9個居群供試材料基本采集信息

Table 1 Basic collection information of materials tested from nine populations ofAtractylodeslancea(Thunb.) DC.

編號No．采集地Locality株數(shù)Numberofindividuals采集號No．ofcollectionAL1江蘇句容茅山MaoshaninJurongofJiangsu20201510021-201510040AL2江蘇句容射烏山ShewushaninJurongofJiangsu10201510011-201510020AL3江蘇句容小九華山XiaojiuhuashaninJurongofJiangsu10200606001-200606010AL4安徽金寨燕子河YanziheinJinzhaiofAnhui10200606021-200606030AL5安徽金寨天堂寨TiantangzhaiinJinzhaiofAnhui10201508011-201508020AL6河南嵩縣五馬寺WumasiinSongxianofHe’nan10201509001-201509010AL7河南嵩縣紙坊ZhifanginSongxianofHe’nan10201509011-201509020AL8湖北?？岛笃篐oupinginBaokangofHubei10200311001-200311010AL9湖北紅安七里坪QilipinginHonganofHubei10201510001-201510010

1.2.2 5個DNA片段序列測定

1.2.2.1 引物信息 對rbcL、matK、psbA-trnH、ITS和nadF5個DNA片段進行擴增。其中，前4個DNA片段為國際通用的條形碼片段，為區(qū)分物種或居群較常用的DNA片段；nadF片段為作者自行篩選的DNA片段，變異率較高，適合進行種下水平的區(qū)分。rbcL、matK、psbA-trnH和ITS片段的擴增引物選擇通用引物，根據(jù)GenBank數(shù)據(jù)庫中發(fā)布的菊科cpDNA序列，利用Oligo 6軟件自行設(shè)計nadF片段的引物，各引物信息見表2。

表2 南蒼術(shù)9個居群的PCR擴增片段引物信息

Table 2 Primer information of PCR amplified fragment from nine populations ofAtractylodeslancea(Thunb.) DC.

片段Fragment引物Primer引物序列(5′→3′)Primersequence(5′→3′)來源SourcerbcL1FATGTCACCACAAACAGAAACFa，etal．[22]724RTCGCATGTACCTGCAGTAGCmatKtrnK-3914FGCTCCAGAAGATGTTGATCGBayer，etal．[23]matK-1541RGGGGTTGCTAACTCAACGGpsbA-trnHpsbA-FGTTATGCATGAACGTAATGCTCSang，etal．[24]psbA-trnH-RCGCGCATGGTGGATTCACAATCCTate，etal．[25]ITSny47-FAACAAGGTTTCCGTAGGTGAMark，etal．[26]AST1TGAGGACGCTTCTCCAGACnadFrpl132-FCCGATTCACCAACTCTTATCTATTTCTArpl132-RCATTGGATGTGAAAGACATGTATTGTTC

1.2.2.2 序列擴增、測序及數(shù)據(jù)處理 每個居群選取5個單株用于DNA片段擴增和序列測定。擴增反應(yīng)在PerkinElmer GeneAmp PCR System 9600(美國PerkinElmer公司)上進行，試劑購自廣州東盛生物科技有限公司。PCR反應(yīng)體系總體積為50.0 μL，包含模板DNA 1.0 μL，2×Reaction Mix(含20 mmol·L-1Tris-HCl、100 mmol·L-1KCl、3 mmol·L-1MgCl2、400 μmol·L-1dNTPs和0.02 g·mL-1溴酚藍)25.0 μL，10 μmol·L-1引物各2.0 μL，2.5 U·μL-1TaqDNA聚合酶1.0 μL，最后用雙蒸水補至50.0 μL。matK、rbcL、psbA-trnH、ITS和nadF片段的擴增程序均為：95 ℃預(yù)變性5 min；95 ℃變性45 s，56 ℃退火45 s，70 ℃延伸90 s，35個循環(huán)；最后72 ℃延伸10 min。PCR產(chǎn)物經(jīng)質(zhì)量體積分數(shù)0.8%～1.0%瓊脂糖凝膠電泳后，溴化乙錠(EB)染色，用Gel Doc圖像分析儀(江蘇省捷達科技發(fā)展有限公司)觀察。

上述PCR擴增產(chǎn)物交由上海華大基因科技有限公司進行雙向測序。測序結(jié)果利用Sequencher 4.5軟件進行拼接，去除低質(zhì)量的序列；通過MEGA 5.1對序列進行分析，用NJ(Neighbor-Joining)法構(gòu)建系統(tǒng)進化樹評估南蒼術(shù)各居群間的親緣關(guān)系。

1.2.3 SRAP分析方法

表3 用于江蘇句容茅山居群南蒼術(shù)特異性條帶篩選的SRAP引物

Table 3 SRAP primers used for screening specific bands ofAtractylodeslancea(Thunb.) DC. from Maoshan population in Jurong of Jiangsu

引物Primer引物序列(5′→3′)Primersequence(5′→3′)引物Primer引物序列(5′→3′)Primersequence(5′→3′)SRAP-F1TGAGTCCAAACCGGATASRAP-R5GACTGCGTACGAATTAACSRAP-F2TGAGTCCAAACCGGAGCSRAP-R6GACTGCGTACGAATTGCASRAP-F3TGAGTCCAAACCGGAATSRAP-R7GACTGCGTACGAATTCCASRAP-F4TGAGTCCAAACCGGACCSRAP-R8GACTGCGTACCAATTCGAGSRAP-F5TGAGTCCAAACCGGAAGSRAP-R9GACTGCGTACCAATTCGCCSRAP-F6TGAGTCCAAACCGGACGSRAP-R10GACTGCGTACCAATTCTCCSRAP-F7TGAGTCCAAACCGGACTSRAP-R11GACTGCGTACGAATTTCGSRAP-F8TGAGTCCAAACCGGAGGSRAP-R12GACTGCGTACGAATTGGTSRAP-R1GACTGCGTACGAATTAATSRAP-R13GGCTTGAACGAGTGACTGASRAP-R2GACTGCGTACGAATTTGCSRAP-R14GACTGCGTACGAATTCAASRAP-R3GACTGCGTACGAATTGACSRAP-R15GACTGCGTACGAATTCTGSRAP-R4GACTGCGTACGAATTTGASRAP-R16GACTGCGTACGAATTGTC

1.2.3.1 SRAP引物篩選 參照文獻[27-28]，用于茅山地區(qū)南蒼術(shù)特異性條帶篩選的SRAP引物(表3)由上海華大基因科技有限公司合成，其中正向引物8條，反向引物16條，兩兩組合共有128種組合方式。分別取每居群所有單株的20.0 μL總DNA混勻，作為SRAP引物篩選的模板。在Biometra TAdvanced型PCR儀(德國耶拿分析儀器股份公司)進行PCR反應(yīng)，反應(yīng)體系總體積為20.0 μL，包含DNA模板1.0 μL，2×Reaction Mix(含20 mmol·L-1Tris-HCl、100 mmol·L-1KCl、3 mmol·L-1MgCl2、400 μmol·L-1dNTPs和0.02 g·mL-1溴酚藍)10.0 μL，10 mmol·L-1正向引物和反向引物各0.7 μL，2.5 U·μL-1TaqDNA聚合酶0.4 μL，用雙蒸水補齊至20.0 μL。SRAP分子標記的擴增程序為：94 ℃預(yù)變性5 min；94 ℃變性1 min，35 ℃退火45 s，72 ℃延伸1 min，5個循環(huán)；94 ℃變性1 min，52 ℃退火45 s，72 ℃延伸1 min，35個循環(huán)；72 ℃延伸10 min。擴增產(chǎn)物經(jīng)質(zhì)量體積分數(shù)2.0%瓊脂糖凝膠電泳，EB染色，電壓60 V，電泳1.5 h后凝膠成像系統(tǒng)觀察、拍照，用2 000 bp DNA Ladder作為分子質(zhì)量標記。

1.2.3.2 SRAP-SCAR特異性條帶的克隆及測序選擇江蘇句容茅山居群南蒼術(shù)擴增出的特異性條帶，對特異性條帶進行切膠回收。將回收的目的DNA片段連接到pMD19-T載體〔寶生物工程(大連)有限公司〕，并將連接產(chǎn)物轉(zhuǎn)入大腸桿菌DH5α感受態(tài)細胞〔寶生物工程(大連)有限公司〕。在含有氨芐的細菌培養(yǎng)基上于37 ℃培養(yǎng)12 h。挑取單菌落進行PCR檢驗，將含有目的條帶的克隆交由上海華大基因科技有限公司進行測序。根據(jù)DNA序列分析結(jié)果，設(shè)計1對特異性引物用于后續(xù)實驗。

1.2.3.3 SRAP-SCAR標記的驗證 用合成的SCAR引物對南蒼術(shù)的9個居群100個單株進行驗證。在Biometra TAdvanced型PCR儀進行PCR反應(yīng)，反應(yīng)體系總體積為20.0 μL，包含DNA模板1.0 μL，2×Reaction Mix(含20 mmol·L-1Tris-HCl、100 mmol·L-1KCl、3 mmol·L-1MgCl2、400 μmol·L-1dNTPs和0.02 g·mL-1溴酚藍)10.0 μL，10 mmol·L-1正向引物和反向引物各0.7 μL，2.5 U·μL-1TaqDNA聚合酶0.4 μL，最后用雙蒸水補齊至20.0 μL。擴增程序為：95 ℃預(yù)變性5 min；95 ℃變性45 s，56 ℃退火45 s，72 ℃延伸30 s，30個循環(huán)；72 ℃延伸10 min。PCR擴增產(chǎn)物經(jīng)質(zhì)量體積分數(shù)1.0%瓊脂糖凝膠電泳，EB染色，電壓100 V，電泳0.5 h后凝膠成像系統(tǒng)觀察、拍照，采用2 000 bp DNA Ladder作為分子質(zhì)量標記。

2 結(jié)果和分析

2.1 5個DNA片段的分析結(jié)果

2.1.1 5個DNA片段的特征分析 南蒼術(shù)9個居群rbcL、matK、psbA-trnH、ITS和nadF片段的變異位點見表4。

由表4可以看出：南蒼術(shù)9個居群rbcL片段的變異位點較少，安徽2個居群南蒼術(shù)具有獨特的序列變異位點，在第617位點處安徽2個居群為T，而其余7個居群均為G。psbA-trnH片段在南蒼術(shù)居群間及居群內(nèi)雖然有變異位點，但只是個別單株的變異，無法代表所在地區(qū)的南蒼術(shù)。ITS片段的變異位點最多，但具有區(qū)分度的位點僅有1個，在第66位點處河南2個居群為A，而其余7個居群均為G，能將河南2個居群南蒼術(shù)鑒別出來；江蘇3個居群、安徽2個居群及湖北紅安七里坪居群在第9、第42、第55、第87、第115、第160、第197、第217、第403、第510和第597位點處分別為G、C、G、C、A、A、G、T、A、A、G，河南2個居群為A、T、A、T、G、C、C、C、T、T、A，可以將上述居群區(qū)分開。但是，matK、rbcL、psbA-trnH和ITS 4個片段均不能將江蘇句容茅山居群南蒼術(shù)與其余居群南蒼術(shù)區(qū)分開。

由表4還可以看出：江蘇句容茅山居群南蒼術(shù)在nadF片段的第581位點處為T，其余8個居群為A，可以將江蘇句容茅山居群南蒼術(shù)與其余8個居群的南蒼術(shù)區(qū)分開。另外，在nadF片段的第76和第343位點處，江蘇3個居群為A，其余8個居群為T。在nadF片段的第161位點處，江蘇3個居群和安徽金寨天堂寨居群為C，安徽金寨燕子河居群、河南2個居群和湖北2個居群為T。

表4 南蒼術(shù)9個居群5個DNA片段變異位點的比較1)

Table 4 Comparison on variable sites in five DNA fragments from nine populations ofAtractylodeslancea(Thunb.) DC.1)

編號2)No．2)rbcLmatKpsbA-trnH617783445416056466737168341052113383858688199257258259260261262284285286287288312313350371372AL1-1GGGCAAGGTAGAAAAT———————————ATTT—AL1-2GGGCAAGGTAGAAAAT———————————AGT——AL1-3GGGCAAGGTAGAAAAT——————CAATTATGT—AL1-4GGGCAAGGTAGTTTTT———————————AGT——AL1-5GGGCAAGGTAGAAAAT———————————MTTT—AL2-1GGGCAAGGTAGAAAAT———————————ATT——AL2-2GGGCAAGGTAGAAAAT———————————ATT——AL2-3GGGCAAGGTAGAAAAT———————————ATT——AL2-4GGGCAAGGTAGAAAAT———————————ATT——AL2-5GGGCAAGGTAGAAAAT———————————ATT——AL3-1GGGCAAGGTAGAAAAT———————————ATTT—AL3-2GGGCAAGGTAGAAAAT———————————ATTT—AL3-3GGGCAAGGTAGAAAAT———————————ATTT—AL3-4GGGCAAGGTAGAAAAT———————————ATT——AL3-5GGGCAAGGTAGAAAAT———————————ATTT—AL4-1TGGCAAGGTAAAAAAT———————————ATTT—AL4-2TGGCGACGTAAAAAAT———————————ATTT—AL4-3TGGCAAGATAAAAAAT———————————ATTT—AL4-4TGGCAAGGTAAAAAAT———————————ATTT—AL4-5TGGCAAGGTAGAAAAT———————————ATTT—AL5-1TGGCAAGGTAAAAAAT———————————ATTT—AL5-2TGGCAAGGTAAAAAAT———————————ATTT—AL5-3TGGCAAGGTAAAAAAT———————————ATTT—AL5-4TGGCAAGGTAAAAAAT———————————ATTT—AL5-5TGGCAAGGTAGAAAAT———————————ATTT—AL6-1GGGCAAGGTAGAAAAT———————————ATTT—AL6-2GGGCAAGGTAGAAAAT———————————ATT——AL6-3GGGCAAGGTAGAAAAT———————————ATT——AL6-4GGGCAAGGTAGAAAAT———————————ATTT—AL6-5GGGCAAGGTAGAAAAT———————————ATTT—AL7-1GGGCAAGGTAGAAAAT———————————ATTT—AL7-2GGGCAAGGTAGAAAATATATGA—————ATTTTAL7-3GGGCAAGGTAGAAAAC———————————ATTT—AL7-4GGGCAAGGTAGAAAATATATGA—————ATTTTAL7-5GGGCAAGGTAGAAAAC———————————ATTT—AL8-1GGGCAAGGTAGAAAAT———————————ATTT—AL8-2GGGCAGGGTAGAAAAT———————————AATT—AL8-3GGGGAAGGTAGAAAAT———————————ATTT—AL8-4GGGCAAGGTAGAAAAT———————————ATT——AL8-5GGGCAAGGTAGAAAAT———————————ATTT—AL9-1GTACAAGGTAGAAAAT———————————ATTT—AL9-2GGGCAAGGTAGAAAAT———————————ATTT—AL9-3GGGCAAGGTGGAAAAT———————————ATTT—AL9-4GGGCAAGGTAGAAAAT———————————ATTT—AL9-5GGGCAAGGGAGAAAAT———————————ATTT—

續(xù)表4 Table 4 (Continued)

編號2)No．2)ITSnadF942556683871061151601941972174034644675105977576154161343395411468473581628704705AL1-1GCGGCCCAAGGTACTAGTACCATTACTTTTAL1-2GCGGCCCAAGGTACTAGTACCATTACTTTTAL1-3GCGGYCCAAGGTAYTAGTACCAAT—CTTT—AL1-4GCGGCCCAAGGTAYTAGTACCAAT—CTAT—AL1-5GCGGYCCAAGGTACTAGTACCAAT—CTAT—AL2-1GCGGCCCAAGGTACTAGTACCATT—CAAT—AL2-2GCGGCCCAAGGTACTAGTACCAATACAATTAL2-3GCGGCCCAAGGTACTAGTACCAATACAATTAL2-4GCGGCCCAAGGTACTAGTACCAATACAATTAL2-5GCGGCCCAAGGTACTAGTACCAAT—CAAT—AL3-1GCGGCCCAAGGTACTAGTACCAAT—TAAT—AL3-2GCGGCCCAAGGTACTAGTACCAATATAAT—AL3-3GCGGCCCAAGGTACTAGTACCAATACATT—AL3-4GCGGCCCAAGGTACTAGTACCAAT—CAATTAL3-5GCGGCCCAAGGTACTAGTACCAAT—CATTTAL4-1GCGGCCCAAGGTACTAGTTCTTAT—CATT—AL4-2GCGGCCCAAGGTACTAGTTCTTAT—CATT—AL4-3GCGGCCCAAGGTACTAGTTCTTAT—CATT—AL4-4GCGGCCCAAGGTACTAGTTCTTAT—CAT——AL4-5GCGGCCCAAGGTACTAGTTCTTAT—CATT—AL5-1GCGGCCCAAGGTACTAGTTCCTAT—CATT—AL5-2GCGGCCCAAGGTACTAGTTCCTAT—CAT——AL5-3GCGGCCCAAGGTACTAGTTCCTAT—CATT—AL5-4GCGGCCCAAGGTACTAGTTCCTAT—CATT—AL5-5GCGGCCCAAGGTACTAGTTCCTAT—CATT—AL6-1ATAACTCGCGCCTCCTATTGTTAT—CAT——AL6-2ATAACTCGCTCCTCTTATTCTTAT—CATTTAL6-3ATAACTCGCGCCTCCTATTCTTAT—CAT——AL6-4ATAACTCGCGCCTCCTAGTCTTAT—CATT—AL6-5ATAACTCGCGCCTCCTAGTCTTAT—CATT—AL7-1ATAACTCGCGCCTCCTAGTCTTAT—CAT——AL7-2ATAACTCGCGCCTCCTATTGTTAT—CATT—AL7-3ATAACTCGCGCCTCCTATTCTTAT—CATTTAL7-4ATAACTCGCGCCTCCTATTCTTAT—CAT——AL7-5ATAACTCGCGCCTCCTATTCTTAG—CAT——AL8-1GCGGCCCAAGGTACTAGTTCTTAT—CATT—AL8-2RYRGCYCRMGSYWCYWRTTCTTAT—CATT—AL8-3RYRGCYCRMGSYWCYWRTTCTTAT—CAT——AL8-4RYRGCYCRMGSYWCYWRTTCTTAT—CAA——AL8-5ATAGCTCGCGCCTCCTATTCTTAT—CAA——AL9-1GCGGCCCAAGGTACTAGTTCTTAT—CAAT—AL9-2GCGGCCCAAGGTACTAGTTCTTAT—CAAT—AL9-3GCGGCCCAAGGTACTAGTTCTTAT—CAT——AL9-4GCGGCCCAAGGTACTAGTTCTTAT—CAT——AL9-5GCGGCCCAAGGTACTAGTTCTTAT—CATT—

1)rbcL：rbcL片段變異位點的堿基Bases at variable sites inrbcLfragment；matK：matK片段變異位點的堿基Bases at variable sites inmatKfragment；psbA-trnH：psbA-trnH片段變異位點的堿基Bases at variable sites inpsbA-trnHfragment； ITS： ITS片段變異位點的堿基Bases at variable sites in ITS fragment；nadF：nadF片段變異位點的堿基Bases at variable sites innadFfragment. —： 堿基缺失Base missing.

2)AL1-1，AL1-2，AL1-3，AL1-4，AL1-5： 江蘇句容茅山居群5個單株Five individuals from Maoshan population in Jurong of Jiangsu； AL2-1，AL2-2，AL2-3，AL2-4，AL2-5： 江蘇句容射烏山居群5個單株Five individuals from Shewushan population in Jurong of Jiangsu； AL3-1，AL3-2，AL3-3，AL3-4，AL3-5： 江蘇句容小九華山居群5個單株Five individuals from Xiaojiuhuashan population in Jurong of Jiangsu； AL4-1，AL4-2，AL4-3，AL4-4，AL4-5： 安徽金寨燕子河居群5個單株Five individuals from Yanzihe population in Jinzhai of Anhui； AL5-1，AL5-2，AL5-3，AL5-4，AL5-5： 安徽金寨天堂寨居群5個單株Five individuals from Tiantangzhai population in Jinzhai of Anhui； AL6-1，AL6-2，AL6-3，AL6-4，AL6-5： 河南嵩縣五馬寺居群5個單株Five individuals from Wumasi population in Songxian of He’nan； AL7-1，AL7-2，AL7-3，AL7-4，AL7-5： 河南嵩縣紙坊居群5個單株Five individuals from Zhifang population in Songxian of He’nan； AL8-1，AL8-2，AL8-3，AL8-4，AL8-5： 湖北?？岛笃壕尤?個單株Five individuals from Houping population in Baokang of Hubei； AL9-1，AL9-2，AL9-3，AL9-4，AL9-5： 湖北紅安七里坪居群5個單株Five individuals from Qiliping population in Hongan of Hubei.

南蒼術(shù)9個居群rbcL、matK、psbA-trnH、ITS和nadF5個DNA片段的分析結(jié)果見表5。由表5可以看出：這5個片段的可匹配序列長度分別為643、1 312、455、686和1 169 bp。當空缺(gap)始終作缺失(missing)處理時，rbcL片段有642個保守位點，保守性最高；psbA-trnH片段有434個保守位點，保守性最低；psbA-trnH片段的變異位點最多，有21個，變異性最大；ITS片段的信息位點最多，有13個。

表5 南蒼術(shù)9個居群5個DNA片段信息1)

Table 5 Information of five DNA fragments from nine populations ofAtractylodeslancea(Thunb.) DC.1)

片段FragmentASL/bpNCSNVSNPISrbcL64364211matK13121302101psbA-trnH455434210ITS6866691713nadF1169116094

1)ASL： 可匹配序列長度Alignment sequence length； NCS： 保守位點數(shù)Number of conserved sites； NVS： 變異位點數(shù)Number of variable sites； NPIS： 信息位點數(shù)Number of parsim-informative sites.

AL1： 江蘇句容茅山居群Maoshan population in Jurong of Jiangsu； AL2： 江蘇句容射烏山居群Shewushan population in Jurong of Jiangsu； AL3： 江蘇句容小九華山居群Xiaojiuhuashan population in Jurong of Jiangsu； AL4： 安徽金寨燕子河居群Yanzihe population in Jinzhai of Anhui； AL5： 安徽金寨天堂寨居群Tiantangzhai population in Jinzhai of Anhui； AL6： 河南嵩縣五馬寺居群Wumasi population in Songxian of He’nan； AL7： 河南嵩縣紙坊居群Zhifang population in Songxian of He’nan； AL8： 湖北保康后坪居群Houping population in Baokang of Hubei； AL9： 湖北紅安七里坪居群Qiliping population in Hongan of Hubei. 分支上數(shù)據(jù)為支持率Datums above the branches are support rate.圖1 基于5個DNA片段聯(lián)合構(gòu)建的南蒼術(shù)9個居群的NJ系統(tǒng)進化樹Fig. 1 NJ phylogenetic tree of nine populations of Atractylodes lancea (Thunb.) DC. based on combining five DNA fragments

2.1.2 基于5個DNA片段聯(lián)合的南蒼術(shù)居群間聚類分析 基于5個DNA片段聯(lián)合構(gòu)建的南蒼術(shù)9個居群間的NJ系統(tǒng)進化樹(圖1)顯示：與江蘇句容茅山居群南蒼術(shù)親緣關(guān)系較近的是江蘇句容的小九華山居群和射烏山居群，其次分別是安徽2個居群、湖北2個居群、河南2個居群，推測江蘇句容茅山居群南蒼術(shù)與其余居群的親緣關(guān)系與南蒼術(shù)分布的地理距離成正相關(guān)。

2.2 SRAP-SCAR分析結(jié)果

2.2.1 與茅山地區(qū)南蒼術(shù)特異性相關(guān)的SRAP分子標記篩選 應(yīng)用128對SRAP引物分別對供試南蒼術(shù)9個居群的總DNA混合樣品進行PCR擴增，根據(jù)指紋圖譜進行篩選，結(jié)果顯示：在引物對SRAP-F8/R14的PCR擴增圖譜中，江蘇句容茅山居群在約300 bp處有1條特有條帶，而其余8個居群均無此條帶(圖2)。

M： DNA marker； AL1： 江蘇句容茅山居群Maoshan population in Jurong of Jiangsu； AL2： 江蘇句容射烏山居群Shewushan population in Jurong of Jiangsu； AL3： 江蘇句容小九華山居群Xiaojiuhuashan population in Jurong of Jiangsu； AL4： 安徽金寨燕子河居群Yanzihe population in Jinzhai of Anhui； AL5： 安徽金寨天堂寨居群Tiantangzhai population in Jinzhai of Anhui； AL6： 河南嵩縣五馬寺居群Wumasi population in Songxian of He’nan； AL7： 河南嵩縣紙坊居群Zhifang population in Songxian of He’nan； AL8： 湖北?？岛笃壕尤篐ouping population in Baokang of Hubei； AL9： 湖北紅安七里坪居群Qiliping population in Hongan of Hubei. 箭頭示江蘇句容茅山居群特有條帶The arrow indicates specific bands of population at Maoshan in Jurong of Jiangsu.圖2 引物對SRAP-F8R14的南蒼術(shù)9個居群總DNA的PCR擴增圖譜Fig. 2 PCR amplification pattern of total DNA from nine populations of Atractylodes lancea (Thunb.) DC. by primer pair SRAP-F8R14

2.2.2 茅山地區(qū)南蒼術(shù)特有條帶的序列測定及與茅山地區(qū)南蒼術(shù)相關(guān)的SCAR引物設(shè)計 江蘇句容茅山居群南蒼術(shù)特有條帶所包含的序列是該居群的特有序列，根據(jù)該序列設(shè)計的引物僅在江蘇句容茅山居群南蒼術(shù)單株中擴增得到條帶，而在其余8個居群南蒼術(shù)單株中則無法得到。

對引物對SRAP-F8R14擴增片段的測序結(jié)果顯示該片段長度為323 bp，其核苷酸序列見圖3。

下劃線部分為設(shè)計的SCAR引物The underline parts are the designed SCAR primers.圖3 引物對SRAP-F8R14在江蘇句容茅山居群南蒼術(shù)中擴增出特有片段的DNA序列Fig. 3 DNA sequence of specific fragment amplified by primer pair SRAP-F8R14 from Atractylodes lancea (Thunb.) DC. from Maoshan population in Jurong of Jiangsu

根據(jù)引物設(shè)計的一般原則，設(shè)計1對SCAR引物，正向引物序列為5′-ATCTTTGCTATTAACTATTAAGCTATTCTTCGCTATTCA-3′，反向引物序列為5′-CGATGTTGAAATTGCTGTCTCTGCTTTCTATTC-3′，此對引物擴增出的序列長度為268 bp。

2.2.3 SRAP-SCAR標記在南蒼術(shù)各居群單株中的驗證 用設(shè)計的SCAR引物對南蒼術(shù)9個居群100個單株進行驗證，結(jié)果見圖4。除江蘇句容茅山居群南蒼術(shù)的20個單株在268 bp處擴增出條帶外，其余8個居群的所有單株均未在該處擴增出條帶。

AL1： 江蘇句容茅山居群Maoshan population in Jurong of Jiangsu； AL2： 江蘇句容射烏山居群Shewushan population in Jurong of Jiangsu； AL3： 江蘇句容小九華山居群Xiaojiuhuashan population in Jurong of Jiangsu； AL4： 安徽金寨燕子河居群Yanzihe population in Jinzhai of Anhui； AL5： 安徽金寨天堂寨居群Tiantangzhai population in Jinzhai of Anhui； AL6： 河南嵩縣五馬寺居群Wumasi population in Songxian of He’nan； AL7： 河南嵩縣紙坊居群Zhifang population in Songxian of He’nan； AL8： 湖北?？岛笃壕尤篐ouping population in Baokang of Hubei； AL9： 湖北紅安七里坪居群Qiliping population in Hongan of Hubei.圖4 SRAP-SCAR標記在南蒼術(shù)9個居群單株中的驗證Fig. 4 Identification of individuals from nine populations of Atractylodes lancea (Thunb.) DC. by SRAP-SCAR marker

3 討 論

3.1 DNA片段鑒別道地產(chǎn)區(qū)藥材的可行性及片段篩選

藥材道地性形成的根本原因是長期地理隔離導(dǎo)致的與其他區(qū)域居群產(chǎn)生的遺傳差異，實質(zhì)上是藥材基因組DNA序列出現(xiàn)了變異，因此，不同區(qū)域同一植物反映在DNA片段上的穩(wěn)定差異實際上是各區(qū)域該植物的特質(zhì)，就此產(chǎn)生了鑒別的依據(jù)。但是，不同DNA片段演化速率的快慢存在差異，作為種下水平的鑒別，必須篩選到演化速率快的片段。盡管不同片段在不同類群植物中演化速率也有差異，但總體上看，rbcL和matK片段的演化速率較慢，往往適用于種或者種以上水平的鑒別，適用于種下水平的往往序列較小，鑒定效果不佳；相對而言，ITS和psbA-trnH等片段更適合遺傳差異小的種下水平的鑒定，但是也并非一定有效[29-31]。

鄭麗[32]利用rbcL、matK、psbA-trnH和ITS共4個DNA片段作為標記，對蒼術(shù)屬5種植物進行了分子鑒別，認為僅ITS片段能將該屬供試的5個種區(qū)分開，而其他3個片段則不行，說明該屬遺傳差異較小、序列較保守，常規(guī)片段很有可能無法進行種下水平的區(qū)分，還需要尋找更具鑒別度的快速演化片段。本研究選用的5個片段中，matK和rbcL片段的變異位點較少，rbcL片段僅第617位點能鑒別安徽2個居群南蒼術(shù)，但安徽2個居群間則無法區(qū)分。雖然psbA-trnH片段的變異位點較多，如河南嵩縣五馬寺居群和江蘇句容茅山居群均有少數(shù)單株存在小片段插入，但只是居群內(nèi)個別單株間的變異，無法代表所在地區(qū)的南蒼術(shù)。ITS片段的變異位點最多，但具有區(qū)分度的僅有1個，能將河南2個居群南蒼術(shù)鑒別出來，但河南2個居群間也無法區(qū)分。nadF片段的變異位點較多，鑒別效果最佳，不僅能將道地產(chǎn)區(qū)江蘇句容茅山居群的南蒼術(shù)區(qū)分開，而且還可以將江蘇句容的射烏山居群和小九華山居群與其他居群分開，但不能區(qū)分射烏山居群和小九華山居群。nadF片段由本研究自行篩選得到，位于葉綠體基因組的基因間區(qū)，變異率較高，適合進行種下水平的區(qū)分。此前的分子鑒別、遺傳演化等研究中極少使用nadF片段，該片段為單拷貝片段，長度穩(wěn)定(1 000 bp左右)，兩端序列保守(分別為蛋白編碼基因nadF和rpl32)[33]，引物設(shè)計容易，有望被開發(fā)為進行種下研究的熱點片段。nadF片段在江蘇句容茅山居群南蒼術(shù)的鑒定效率達到100%。

3.2 茅山地區(qū)南蒼術(shù)鑒別的SCAR分子標記的開發(fā)意義

利用隨機引物的指紋圖譜進一步開發(fā)的特異性分子標記是專一性、特異性的分子標記，不僅穩(wěn)定，而且操作十分簡單，已經(jīng)成為檢驗藥材道地性的重要手段。與DNA序列的差異相比，由隨機引物進一步開發(fā)的SCAR分子標記更適用于遺傳背景相近的對象的鑒定。由于隨機引物擴增產(chǎn)生的多態(tài)性條帶能給予更大的篩選范圍，更容易發(fā)現(xiàn)基因組的差異，因此，種下水平的分子標記篩選常用隨機引物的多態(tài)性條帶來開發(fā)，甚至在同一居群不同個體間也能用隨機引物來發(fā)現(xiàn)差異，進而開發(fā)成SCAR分子標記。如韋陽連等[34]采用RAPD引物找到了茅山南蒼術(shù)不同個體間與蒼術(shù)素含量相關(guān)的分子標記。

本研究開發(fā)的江蘇句容茅山居群南蒼術(shù)SRAP-SCAR標記源自不同居群南蒼術(shù)通過SRAP隨機引物擴增后的指紋圖譜的比較篩選，但這種標記不用進行序列比對分析，僅用PCR產(chǎn)物進行凝膠電泳檢測便可以達到鑒別的目的，操作更加簡便，判斷也更加直觀，對操作者的專業(yè)知識要求更低，因此，在應(yīng)用推廣上更具優(yōu)勢。

總體而言，本研究基于nadF片段的序列差異和SRAP的指紋圖譜開發(fā)的茅山地區(qū)南蒼術(shù)特有的分子標記均可以有效地將茅山地區(qū)南蒼術(shù)與其余居群南蒼術(shù)區(qū)分開。茅山地區(qū)南蒼術(shù)作為南蒼術(shù)的道地藥材，其專一性、特異性鑒別的分子標記的開發(fā)對其他道地藥材的保護性鑒別也提供了理論和科學(xué)依據(jù)。

[1] 中國科學(xué)院中國植物志編輯委員會. 中國植物志： 第七十八卷第一分冊[M]. 北京： 科學(xué)出版社， 1987： 23.

[2] 國家藥典委員會. 中國人民共和國藥典： 2015年版(一部)[M]. 北京： 中國醫(yī)藥科技出版社， 2015： 161.

[3] 陶弘景. 本草經(jīng)集注[M]. 輯校本. 北京： 人民衛(wèi)生出版社， 1994： 196-197.

[4] 彭華勝， 王德群. 安徽省術(shù)屬藥用植物過渡類群的居群生物學(xué)研究[J]. 中國中藥雜志， 2007， 32(9)： 793-797.

[5] 朱曉琴， 賀善安. 不同產(chǎn)地蒼術(shù)藥材化學(xué)成分的比較[J]. 植物資源與環(huán)境， 1994， 3(4)： 18-22.

[6] 劉曉寧， 侯芳潔， 谷 巍， 等. 不同產(chǎn)地蒼術(shù)揮發(fā)油特征性成分分析[J]. 南京中醫(yī)藥大學(xué)學(xué)報， 2009， 25(1)： 51-53.

[7] 沈曉婧， 楊文國， 石 瑩， 等. 茅蒼術(shù)與北蒼術(shù)揮發(fā)油的統(tǒng)計分析[J]. 科技創(chuàng)新導(dǎo)報， 2013(34)： 67-68.

[8] 吳寶成， 韋陽連， 高 興， 等. 茅山蒼術(shù)揮發(fā)油組分特征及與根莖鮮質(zhì)量相關(guān)性研究[J]. 林產(chǎn)化學(xué)與工業(yè)， 2008， 28(5)： 84-88.

[9] 余永邦， 秦民堅， 梁之桃， 等. 不同產(chǎn)區(qū)太子參的rDNA ITS區(qū)序列的比較[J]. 植物資源與環(huán)境學(xué)報， 2003， 12(4)： 1-5.

[10] 張宏意， 石祥剛. 不同產(chǎn)地何首烏的ITS序列研究[J]. 中草藥， 2007， 38(6)： 911-914.

[11] 孔德政， 張 偉， 李 永. 道地產(chǎn)區(qū)不同牛膝種群的遺傳多樣性分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版)， 2014， 42(11)： 151-155.

[12] 葛燕芬， 杭悅宇， 夏 冰， 等. 5種蒼術(shù)屬藥用植物的trnL-F序列測定及種間遺傳關(guān)系分析[J]. 植物資源與環(huán)境學(xué)報， 2007， 16(2)： 12-16.

[13] PENG H S， YUAN Q J， LI Q Q， et al. Molecular systematics of genusAtractylodes(Compositae， Cardueae)： evidence from internal transcribed spacer (ITS) andtrnL-Fsequences[J]. International Journal of Molecular Sciences， 2012， 13： 14623-14633.

[14] 袁英惠. 利用微衛(wèi)星標記分析球花石斛的遺傳多樣性并鑒定其藥材的地理來源[D]. 南京： 南京師范大學(xué)生命科學(xué)學(xué)院， 2011： 37-42.

[15] 高文遠， 秦恩強， 肖小河， 等. 當歸藥材道地性的RAPD分析[J]. 中草藥， 2001， 32(10)： 926-929.

[16] 沈 潔， 丁小余， 丁 鴿， 等. 鐵皮石斛居群差異的研究Ⅱ： ISSR指紋標記方法的建立與優(yōu)化[J]. 中國中藥雜志， 2006， 31(4)： 291-294.

[17] 施云燕， 黃海飛， 朱振洪. 貝母類藥材隨機引物擴增多態(tài)性分析研究[J]. 現(xiàn)代中藥研究與實踐， 2014， 28(4)： 10-13.

[18] 田孟良， 官 宇， 劉 帆， 等. 基于RAPD標記的SCAR分子標記技術(shù)鑒定川牛膝及其混淆品[J]. 中國中藥雜志， 2010， 35(8)： 953-956.

[19] SZE S C W， SONG J X， WONG R N S， et al. Application of SCAR (sequence characterized amplified region) analysis to authenticateLyciumbarbarum(wolfberry) and its adulterants[J]. Biotechnology and Applied Biochemistry， 2008， 51： 15-21.

[20] 任冰如， 賀善安， 於 虹， 等. 用RAPD技術(shù)評估蒼術(shù)居群間的親緣關(guān)系[J]. 中草藥， 2000， 31(6)： 458-461.

[21] 郭蘭萍， 黃璐琦， 王 敏， 等. 南北蒼術(shù)的RAPD分析及其劃分的初步探討[J]. 中國中藥雜志， 2001， 26(3)： 156-158.

[22] FAY M F， SWENSEN S M， CHASE M W. Taxonomic affinities ofMedusagyneoppositifolia(Medusagynaceae)[J]. Kew Bulletin， 1997， 52： 111-120.

[23] BAYER R J， GREBER D G， BAGNALL N H. Phylogeny of Australian Gnaphalieae (Asteraceae) based on chloroplast and nuclear sequences， thetrnLintro，trnL/trnFintergenic spacer，matK， and ETS[J]. Systematic Botany， 2002， 27： 801-814.

[24] SANG T. Utility of low-copy nuclear gene sequences in plant phylogenetics[J]. Critical Reviews in Biochemistry and Molecular Biology， 2002， 37： 121-147.

[25] TATE J A， SIMPSON B B. Paraphyly ofTarasa(Malvaceae) and diverse origins of the polyploid species[J]. Systematic Botany， 2003， 28： 723-737.

[26] CHASE M W， COWAN R S， HOLLINGSWORTH P M， et al. A proposal for a standardised protocol to barcode all land plants[J]. Taxon， 2007， 56： 295-299.

[27] GENG J F， ZHU C S， ZHANG X W， et al. A genetic linkage map of non-heading Chinese cabbage[J]. Journal of the American Society for Horticultural Science， 2007， 132： 816-823.

[28] LI G， QUIROS C F. Sequence-related amplified polymorphism (SRAP)， a new marker system based on a simple PCR reaction： its application to mapping and gene tagging inBrassica[J]. Theoretical and Applied Genetics， 2001， 103： 455-461.

[29] CBOL Plant Working Group， JANZEN D H. A DNA barcode for land plants[J]. Proceedings of the National Academy of Sciences of the United States of America， 2009， 106： 12794-12797.

[30] KRESS W J， ERICKSON D L. DNA barcodes： genes， genomics， and bioinformatics[J]. Proceedings of the National Academy of Sciences of the United States of America， 2008， 105： 2761-2762.

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

[32] 鄭 麗. 蒼術(shù)屬系統(tǒng)發(fā)育及藥用植物白術(shù)的群體遺傳和栽培起源研究[D]. 杭州： 浙江大學(xué)生命科學(xué)學(xué)院， 2013： 124-125.

[33] DANIELL H， LIN C S， YU M， et al. Chloroplast genomes： diversity， evolution， and applications in genetic engineering[J]. Genome Biology， 2016， 17： 134.

[34] 韋陽連， 吳寶成， 周義鋒. 茅山蒼術(shù)與蒼術(shù)素相關(guān)的RAPD標記研究[J]. 安徽農(nóng)業(yè)科學(xué)， 2008， 36(24)： 10368-10370.

(責(zé)任編輯： 張明霞)

Study on specific molecular marker ofAtractylodeslanceain Maoshan area

WANG Mingqiang1，①， TANG Shijie1，①， FANG Yan2， XUE Jiayu1，②

(1. Center for Plant Diversity and System Evolution， Institute of Botany， Jiangsu Province and Chinese Academy of Sciences， Nanjing 210014， China； 2. Nanjing Forest Police College， Nanjing 210023， China)，J.PlantResour. &Environ.， 2017， 26(2)： 17-26

The study focuses on nine populations from all distribution areas of wildAtractylodeslancea(Thunb.) DC. in China， including Maoshan population， Shewushan population and Xiaojiuhuashan population in Jurong of Jiangsu， Yanzihe population and Tiantangzhai population in Jinzhai of Anhui， Wumasi population and Zhifang population in Songxian of He’nan， and Houping population in Baokang and Qiliping population in Hongan of Hubei. Taking these nine populations as objects， specific molecular marker ofA.lanceafrom Maoshan population in Jurong of Jiangsu is screened by using sequence alignment of fragments ofrbcL，matK，psbA-trnH， ITS andnadFand comparison on SRAP fingerprints. The results show that through sequence alignment of five DNA fragments， the base at the 581th locus innadFfragment ofA.lanceafrom Maoshan population in Jurong of Jiangsu is A， and that from other eight populations is T， which can distinguish Maoshan population in Jurong of Jiangsu from other eight populations. A specific band with length of 323 bp inA.lanceafrom Maoshan population in Jurong of Jiangsu is screened by SRAP random primers， and a pair of SCAR primers is designed according to the sequence of this band.The forward primer sequence is 5′-ATCTTTGCTATTAACTATTAAGCTATTCTTCGCTATTCA-3′，and the reverse primer sequence is 5′-CGATGTTGAAATTGCTGTCTCTGCTTTCTATTC-3′， which can amplify specific band with length of 268 bp from all individuals ofA.lanceafrom Maoshan population in Jurong of Jiangsu， and cannot amplify any band from other populations. The two specific markers screened in this study can provide reliable technical support for identifying the genuineness ofA.lancea.

Atractylodeslancea(Thunb.) DC.； genuineness； Maoshan area；nadFfragment； SRAP-SCAR

2016-12-08

江蘇省林業(yè)三新工程項目(LYSX[2015]22)

王明強(1992—)，男，安徽阜陽人，碩士研究生，主要從事植物系統(tǒng)與演化方面的研究。 湯詩杰(1969—)，男，安徽東至人，博士，副研究員，主要從事園林植物與植物多樣性保護方面的研究。

Q946-33； S567.21+1

A

1674-7895(2017)02-0017-10

10.3969/j.issn.1674-7895.2017.02.03

①共同第一作者

②通信作者E-mail： xuejy@cnbg.net