SCoT分子標記對不同產(chǎn)地白芍親緣關系的研究△

席秀利，胡珊，鄔龍怡，鐘旺琴

(廣州市香雪制藥股份有限公司，廣東 廣州 510663)

白芍，別名芍藥，《中華人民共和國藥典》2015版載其原植物為毛茛科芍藥屬植物芍藥PaeonialactioraPall.的干燥根。其性味苦、酸，微寒。歸肝、脾經(jīng)。養(yǎng)血調(diào)經(jīng)，斂陰止汗，柔肝止痛，平抑肝陽。用于血虛萎黃，月經(jīng)不調(diào)，自汗，盜汗，脅痛，腹痛，四肢攣痛，頭痛眩暈等癥[1]。白芍主產(chǎn)于中國安徽、四川、浙江、山東等地[2]。浙江磐安、安徽亳州、四川中江和山東菏澤作為白芍的傳統(tǒng)產(chǎn)區(qū)，在長期的引種和交流中，種質(zhì)資源較為混雜，導致白芍種質(zhì)資源的遺傳結構有較大程度的差異，品種退化嚴重，栽培種質(zhì)混雜[3]。白芍是中國傳統(tǒng)常用大宗中藥材，需求量大，且品種繁多，而種質(zhì)資源的混雜嚴重影響了白芍規(guī)范化生產(chǎn)的發(fā)展。因此有必要從源頭即白芍種質(zhì)資源、種植規(guī)范方面嚴格要求，搞清白芍藥材的地道性，建立品種質(zhì)量標準，從而指導GAP中藥材生產(chǎn)質(zhì)量管理規(guī)范標準化生產(chǎn)。

目標起始密碼子多態(tài)性分子標記(Start codon targeted polymorphism，SCoT)是Collard和Mackill開發(fā)的基于SPAR(單引物擴增反應)的新的目的基因分子標記方法。其原理是根據(jù)植物基因中ATG翻譯起始位點側(cè)翼序列的保守性設計單引物，對基因組進行擴增，產(chǎn)生偏向候選功能基因區(qū)的顯性多態(tài)性標記[4]。SCoT標記結合了ISSR標記和RAPD標記的優(yōu)點，操作簡單，成本低，引物設計簡單且可以通用，多態(tài)性豐富，可用于分子輔助育種[5]。目前SCoT已被成功應用于柑橘[6]、鐵皮石斛[7]等作物，在種屬間、品種間的遺傳多樣性分析及雜交后代的遺傳分析等研究中取得較多進展。本研究利用SCoT分子標記對四川中江、浙江磐安、安徽亳州、山東菏澤四大傳統(tǒng)產(chǎn)地的22個居群的白芍樣品進行親緣關系研究和遺傳多樣性檢測。目的是從分子水平上研究22個白芍居群的親緣關系及遺傳多樣性，為白芍藥材的引種和栽培提供科學基礎。

1 材料與方法

1.1 材料

采集來自浙江，安徽，四川，山東共22個不同居群白芍樣品，經(jīng)廣州中醫(yī)藥大學黃海波副教授鑒定為白芍PaeonialactioraPall.，憑證標本保存于香雪制藥有限公司研發(fā)中心。詳細記錄樣品信息及編號，用75%乙醇清潔葉表面后保存于-80 ℃冰箱備用，具體樣品信息如下。

表1 樣品信息表

1.2 儀器和試劑

分析天平，核酸蛋白測定儀(Implen Nanophotometer)，PCR儀(BIO-RAD)，離心機(Eppendorf)，超低溫冰箱(Thermo 902-ULTS)，電泳儀、電泳槽(BIO-RAD)，制冰機，水浴鍋，凝膠成像儀(BIO-RAD)。

dNTPs、ExTaqDNA聚合酶、Mg2+、10x ExTaq Buffer、Loing Buffer(Takala)，Maker3(廣州東盛生物科技有限公司)，瓊脂糖(Biowest)，Golden View核酸染料(北京博凌科為生物科技有限公司)，TAE(50×)(生工)，植物基因組提取試劑盒(天根生化)，參考Collard和Mackill提供的36條SCoT引物序列。

2 方法

2.1 基因組DNA提取與檢測

采用天根植物基因組試劑盒提取所有材料的基因組DNA。用微量紫外核酸儀檢測DNA濃度和純度。將合格樣品DNA濃度稀釋至20 ng·μL-1后置于-20 ℃保存?zhèn)溆谩?/p>

2.2 SCoT-PCR 反應體系正交優(yōu)化

選用L16(54)正交試驗設計(表2)，對Mg2+、dNTPs濃度、ExTaq DNA聚合酶、引物和DNA用量進行5因素4水平的正交試驗。以山東編號為J的DNA為模板，隨機挑選S6號引物作為擴增引物，反應總體積20 μL，含有10×PCR Buffer 2 μL，每個處理重復2次。

表2 PCR反應正交體系濃度表

2.3 SCoT-PCR引物篩選

根據(jù)正交試驗結果選擇最佳反應體系對36條引物進行篩選，設平行組。

2.4 退火溫度篩選

對篩選后的引物進行溫度篩選，在PCR儀上設定50 ℃～60 ℃溫度區(qū)間，即50 ℃，50.9 ℃，52.3 ℃，54 ℃，56.3 ℃，58.1 ℃，59.3 ℃，60 ℃，8個溫度梯度進行退火溫度篩選。

2.5 SCoT-PCR反應體系的驗證與多態(tài)性引物篩選

選擇浙江、四川、安徽三個不同產(chǎn)地的白芍DNA為模板，按照上述優(yōu)化的SCoT-PCR體系進行SCoT分子標記，篩選出多態(tài)性好、條帶清晰的引物。

2.6 SCoT-PCR 擴增與檢測

PCR擴增反應在PCR儀上進行，擴增程序為：94 ℃預變性5 min；94 ℃變性45 s，55 ℃退火45 s，72 ℃延伸90 s，35個循環(huán)；72 ℃延伸10 min，最后4 ℃保存。擴增反應結束后，取10 μL擴增產(chǎn)物加入2 μL 6×loading buffer用2.0%瓊脂糖凝膠，在1XTAE電極緩沖液，120 V，30 min電泳下檢測條帶，凝膠成像儀下拍照。

2.7 數(shù)據(jù)分析

根據(jù)PCR擴增產(chǎn)物的電泳結果，對清晰且易于辨認的條帶采用“0/1”系統(tǒng)記錄其位置，建立SCoT標記的0/1矩陣。用NTSYS-pc2.10 e軟件計算樣品間的相似系數(shù)，并用不加權成對算術平均法(UPGMA)進行聚類分析。

3 實驗結果

3.1 DNA提取結果

將提取的所有樣品在超微量紫外分光光度計下測量其濃度和純度，得到其A260/A280均在1.7～2.0之間，說明DNA純度較高，符合PCR技術要求。

3.2 SCoT-PCR反應體系正交優(yōu)化結果分析

使用引物S6對反應體系優(yōu)化的電泳結果見圖1，參考桂騰琴[14]直觀評價法依次對16個實驗組結果多態(tài)性進行評分，條帶數(shù)量豐富、清晰的計16分，涂布或無條帶的計1分，反應1～16組體系的平均整數(shù)計分依次為：4、8、5、9、7、12、7、6、8、9、3、10、5、15、16、14分。依照計分原則，以第15組計分最高，條帶最為清晰明亮。因此考慮體系最終為1×PCR buffer，3 mmol·L-1Mg2+，0.2 mol·L-1dNTPs，0.6 μmol·L-1primer，0.75 U ExTaq酶，30ngDNA模板。

注：1～16分別為正交1～16組不同PCR反應體系。圖1 SCoT正交體系優(yōu)化電泳圖

3.3 引物篩選及退火溫度優(yōu)化

以最佳反應體系對36條SCoT通用引物進行篩選，重復2次初篩選出14條引物見下圖2，進行8個退火溫度梯度篩選，根據(jù)信號強度判斷條帶，采用區(qū)間溫度對篩選的引物進行匹配最終獲得在53 ℃，56 ℃，57 ℃，59 ℃四個溫度擴增下的13條引物。其中S29和S30號引物的退火溫度梯度篩選如下圖3。

注：1～36分別為S1～S36號引物。圖2 36條引物初篩電泳圖

注：1～8分別為50 ℃～60 ℃梯度溫度。圖3 S29和S30在8個溫度梯度下擴增電泳圖

3.4 SCoT擴增的多態(tài)性分析

以浙江，四川，安徽三個地域差異大的樣品對退火優(yōu)化后的13條引物進行多態(tài)性篩選，最終獲得12條擴增帶型完整、清晰、多態(tài)性豐畗的引物(表3)。利用這12條引物對24份材料進行SCoT擴增，共獲得總條帶121條，其中多態(tài)性條帶共77條，多態(tài)位點百分率共占63.6%，多態(tài)性比率較高，在分子水平證實了白芍具有豐富的種內(nèi)遺傳多樣性。

注：a.安徽白芍；b.四川白芍；c.浙江白芍。圖4 部分SCoT引物多態(tài)性篩選

表3 SCOT引物擴增結果

注：1～22為樣品表中依次排序的A1～J號樣品。圖5 S6引物對22份材料SCoT-PCR擴增電泳圖

注：1～22為樣品表中依次排序的A1～J號樣品。圖6 S25引物對22份材料SCoT-PCR擴增電泳圖

3.5 SCoT標記的遺傳相似性及聚類分析

利用NTSYS軟件構建全部材料基于SCoT數(shù)據(jù)的UPGMA樹狀圖，結果如圖7。22份材料兩兩之間的相似系數(shù)在0.78～0.94之間，以相似系數(shù)0.82為閾值可將22份材料劃分為5組，浙江的4個樣品Z1-Z4聚為一支；四川的4個樣品C3、C4、C5、C6聚為一支；安徽A1、A2與四川C1、C2聚為一支，說明安徽與四川白芍親緣關系較近，從地理位置分析可能由于白芍地方引種和雜交導致，因此安徽白芍與四川白芍差異不大；山東B、C、J三個樣品聚為一支，且與四川和安徽白芍親緣關系較近，可能因山東白芍多來源于栽培，而栽培品種來源混雜，從聚類樹上可以初步認為部分山東白芍由四川或者安徽引種得來；剩下7個的山東樣品則聚為一組。

圖7 22個不同產(chǎn)地白芍樣品的聚類分析圖

4 討論

近年來，采用分子生物技術對藥材鑒別已成為趨勢，目前白芍及品種的分子鑒定僅出現(xiàn)AFLP[9]、RAPD[10-11]、ISSR[12-13]、DNA條形碼[14]的研究，然而尚未有關SCoT標記在白芍親緣關系鑒定中應用的研究報道。其中周佐斌[11]使用ISSR分析對白芍種質(zhì)資源的遺傳多態(tài)性性分析，主要對浙江磐安白芍為主要研究對象，并以花色特征的種內(nèi)變異為研究內(nèi)容，發(fā)現(xiàn)浙江磐安白芍與外省居群分支明顯。徐攀[13]等人采用DNA條形碼技術和隨機多態(tài)性RAPD分子標記技術對12個不同品種的白芍種質(zhì)資源進行遺傳多樣性研究，發(fā)現(xiàn)DNA條形碼分子標記在種內(nèi)水平中并不能很好的區(qū)分相關的品種和品系，而RAPD分子標記的結果也證明了現(xiàn)有的白芍種質(zhì)資源中存在較為豐富的遺傳多樣性。

SCoT標記作為一種新型的目的基因分子標記，結合了ISSR和RAPD的優(yōu)點，操作簡單、重復性好、引物通用性強，實驗成本較低。本實驗從36條SCoT通用引物中最終篩選出12條清晰、條帶豐富的作為不同產(chǎn)地白芍的SCoT分子標記的引物，驗證了其穩(wěn)定性及可靠性。引物擴增多態(tài)位點百分率占63.6%，在分子水平證實了不同產(chǎn)地白芍具有豐富的種內(nèi)遺傳多樣性，本研究通過NTSYS軟件進行相似系數(shù)和UPGMA聚類分析，實現(xiàn)了對不同產(chǎn)地白芍22份樣品的初步區(qū)分。不同地區(qū)白芍樣品在聚類樹上基本聚為一類，浙江、安徽白芍在聚類圖上的聚類相對明顯，且同產(chǎn)地內(nèi)的白芍相似系數(shù)高。而安徽白芍與四川部分白芍則出現(xiàn)混種現(xiàn)象，同時部分山東白芍樣品與安徽，四川兩地部分樣品親緣關系近，這可能由于白芍經(jīng)過短期引種和馴化導致品種間混雜現(xiàn)象。這對比研究結果表明SCoT技術能從分子水平上檢測不同產(chǎn)地和居群之間白芍的基因差異，可以從基因水平把不同地區(qū)植株外型相似的白芍品種區(qū)分開，對新品種的選育有很大的意義。

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

[2] 查良平，楊俊，彭華勝，等.四大產(chǎn)地白芍的種質(zhì)調(diào)查.中藥材.2011，34(7)：1037-1040.

[3] 徐攀，梁衛(wèi)青，周潔，等.白芍種質(zhì)資源的DNA分子標記研究 [J].中華中醫(yī)藥學刊，2017，35(4)：49-53.

[4] Collard B C Y，Mackill D J.Start codon targeted(SCoT)polymorphism：a simple，novel DNA marker technique for generating genetargeted markers in plants[J].Plant Mol Biol Rep，2009，27(1)：86-93.

[5] 熊發(fā)前，唐榮華，陳忠良，等.目標起始密碼子多態(tài)性(SCoT)：一種基于翻譯起始位點目的基因標記技術[J].分子植物育種，2009，7(3)：635-638.

[6] 蔣巧巧，龍桂友，李武文，等.SCoT結合克隆測序鑒別湖南甜橙變異類型[J].中國農(nóng)學通報，2011，27(6)：148-154.

[7] 趙瑞強，高燕會，章曉玲，等.鐵皮石斛SCoT-PCR 反應體系構建及優(yōu)化[J].核農(nóng)學報，2012，26(4)：648-655.

[8] 朱紅霞.牡丹、芍藥品種DNA指紋圖譜繪制的初步研究[D].北京：北京林業(yè)大學，2004.

[9] 陳丙鑾.杭悅宇，周義峰，等.白芍傳統(tǒng)品種原植物遺傳多樣性及親緣關系分析[J].植物資源與環(huán)境學報，2003，12(2)：17-21.

[10] 來明達，陳瑩，朱振洪，等.不同居群白芍的DNA隨機擴增多態(tài)性研究 [J].中國現(xiàn)代醫(yī)藥雜志，2012，14(12)：4-6.

[11] 周佐斌.白芍原植物種質(zhì)資源遺傳多樣性的ISSR分析[D].南昌：南昌大學，2007.

[12] 蔣雨晗，張麗萍，周學剛，等.ISSR分子標記在白芍親緣關系的應用研究[J].世界科學技術，2011，18(3)：429-433.

[13] 徐攀，梁衛(wèi)青，周潔，等.白芍種質(zhì)資源的DNA分子標記研究[J].中華中醫(yī)藥學刊，2017，35(4)：429-433.

[14] 桂騰琴，孫敏，喬愛民，等.正交設計優(yōu)化果梅 ISSR 反應體系[J].果樹學報，2009，26(1)：108-112.