黃連屬部分藥用植物遺傳多樣性與親緣關系的SCoT分析

陳大霞+王鈺+張雪+潘媛+李隆云

[摘要]采用SCoT分子標記分析4種黃連屬藥用植物的遺傳多樣性和親緣關系，并用POPGENE計算遺傳參數，TREECONW軟件進行聚類分析。結果表明：28條引物共檢測到434條帶，其中多態性條帶430條；各遺傳參數均表明黃連屬藥用植物種間的遺傳多樣性極為豐富（種間PPB=99.1%，N_a=1.990 6，N_e=1.329 3，H=0.212 2，I=0.337 8），但種內遺傳多樣性較低（種內PPB=16.8%，N_a=1.168 2，N_e=1.073 0，H=0.043 7，I=0.067 7）；遺傳分化系數G_st=0.794 0，說明79.4%的遺傳變異來源于種間；同種黃連屬藥用植物的不同植株均聚在一起，證實SCoT分子標記可作為黃連屬藥用植物親緣關系研究的有效方法之一。

[關鍵詞]黃連屬； SCoT； 藥用植物； 親緣關系

[Abstract]The genetic diversity and genetic relationship among four medicinal species of Coptis were detected by the approach of sequence-related amplified polymorphism （SCoT）. The associated genetic parameters were calculated by POPGENE1.31. The systematic diagram of genetic relationship were clustered by TREECONW. The results showed that a total of 434 bands were produced by using 28 primers， of which 430 were polymorphic loci. There was a high level of genetic diversity among species （PPB=99.1%，N_a=1.990 6，N_e=1.329 3，H=0.212 2，I=0.337 8）. However， genetic diversity was lower within species， the average of genetic parameters wasPPB=16.8%，N_a=1.168 2， N_e=1.073 0，H=0.043 7，I=0.067 7. The Nei′s genetic differentiation coefficient was 0.794 0， that indicated that most of the genetic variation existed among species. By clustering analysis， different individuals gathered in the same group. The results confirmed that SCoT marker can be used as one of the effective methods to reveal the genetic diversity and relationship among medicinal species of Coptis.

[Key words]SCoT； Coptis； medicinal plant； genetic relationship

黃連為名貴常用中藥之一，始載于《神農本草經》，在我國已有2 000多年的藥用歷史，具有清熱燥濕，瀉火解毒的功效。我國分布的黃連屬植物中黃連Coptis chinensis Franch.、三角葉黃連C. deltoidea C.Y.Cheng et Hsiao.及云南黃連C. teetoides C.Y.Cheng.是《中國藥典》各版本收錄藥材黃連的基原植物 [2]。此外，分布于四川峨眉、雅安與洪雅地區的峨嵋野連C. omeiemsis （Chen） C.Y.Cheng.、西藏南部昌都地區的西藏黃連C. teeta Wall.及四川馬邊一帶的線萼黃連C. linearisepala T.Z.Wang et C.K.hsieh.等多種栽培或野生黃連屬植物的根莖作為藥用?！度毡舅幘址健烦蛰d了《中國藥典》收錄藥材黃連的基原植物外，還收載了日本黃連C.japonica Makino.。因此，理清黃連屬法定用種及地方用種的親緣關系，準確進行基原鑒別對確保黃連的臨床用藥安全具有重要意義。

黃連屬藥用植物親緣關系的研究及鑒別早期主要側重于外部形態特征^[1]、數量分類^[2]及核型研究^[3-4]。近年來，將分子標記技術應用于黃連屬藥用植物遺傳多樣性與親緣關系的研究報道逐漸增多。本文以黃連屬的4種藥用植物為研究對象，采用SCoT（sequence-related amplified polymorphism，相關序列擴增多態性標記）分子標記技術^[5]研究黃連屬藥用植物種間的遺傳特征，探討物種間的親緣關系，為黃連屬藥用植物的鑒定建立新方法，為形態分類提供分子水平的依據。

1 材料和方法

1.1 供試黃連 供試的黃連屬4種藥用植物見表1，其中黃連為長期人工栽培種，峨眉黃連與三角葉黃連為就地保護撫育與栽培，日本黃連為引進栽培。采集植株的幼嫩葉片洗凈晾干后，放入裝有足量硅膠的自封袋快速干燥?；氐綄嶒炇液?，取出葉片提取基因組DNA。每種藥用植物均采集24株。

1.2 儀器與試劑 PCR儀（S1000 Thermal Cycler，美國BIO-RAD）、凝膠成像系統（Gel Doc XR，美國BIO-RAD）、超微量分光光度計（ND-2000，美國Thermo Scientific）、電泳儀（DYY-6C型，北京六一）、電泳槽（DYCP-31F，北京六一）。

SCoT引物（上海生工）、植物基因組DNA提取試劑盒（北京天根）、Trans2K Plus DNA Marker （北京全式金）、TaqDNA聚合酶（日本TOYOBO）、dNTPs（日本TOYOBO）、GelRed核酸凝膠染料（美國Biotium）、瓊脂糖（進口分裝），其余試劑為國產分析純。

1.3 基因組DNA的提取 采用植物基因組DNA提取試劑盒提取4種黃連屬藥用植物，共計96個植株的基因組DNA，濃度和純度采用超微量分光光度計進行測定。擴增時，將DNA濃度稀釋至40 mg·L^-1。

1.4 PCR擴增與電泳檢測 從Collard和Mackill^[5]開發的SCoT引物中篩選出28條擴增條帶清晰的引物用于SCoT-PCR擴增，引物序列見表2。SCoT-PCR擴增反應總體積為15 μL，內含1×PCRbuffer，200 μmol·L^-1 dNTP，1.5 mmol·L^-1 Mg²⁺，0.4 μmol·L^-1引物，1 U TaqDNA聚合酶和40 ng模板，不足部分用滅菌的超純水補足。PCR擴增程序為：94 ℃預變性4 min；然后94 ℃變性1 min，50 ℃退火1 min，72 ℃延伸2 min，36 個循環；最后72 ℃延伸10 min，4 ℃保存。

PCR擴增產物采用瓊脂糖凝膠電泳進行檢測。取5 μL PCR擴增產物與2.5 μL 6×DNA 上樣緩沖液（內含GelRed核酸凝膠染料）混勻后，以1×TAE為電泳緩沖液在1.5%的瓊脂糖凝膠上電泳分離，電壓150 V。當溴酚藍指示劑距凝膠前沿大約2～3 cm時停止電泳，取出凝膠在成像系統上觀察、照相與保存。

1.5 數據的統計與分析 根據擴增條帶的電泳遷移率及其有無統計數據，有帶記“1”，無帶記“0”。聚類分析采用TREECONW軟件進行，遺傳參數采用POPGENE version軟件計算，包括多態位點百分率（percentage of polymorphic bands，PPB）、表觀等位基因數量（observed number of alleles， N_a）、等位基因有效值（Effective number of alleles，N_e）、Nei′s基因多樣度指數（Nei′s gene diversity，H）、Shannon′s多樣性信息指數（I）、總的基因多樣度（Whole genetic diversity，H_t）、種內基因多樣度（genetiv diversity within population，H_s）、種間基因多樣性（genetic diversity among populations，D_st）、基因分化系數（genetic differentiation，G_st）、基因流（gene flow，N_m）。

2 結果與分析

2.1 SCoT標記的擴增多態性 篩選出28條擴增譜帶清晰并呈現多態性的引物用于黃連屬4種藥用植物的親緣關系分析， SP21的擴增圖譜見圖1。28條SCoT引物共擴增出434條帶，其中有430條為多態性條帶，占總擴增條帶的99.1%。每個引物擴增條帶數在7～22條，平均15.5條，其中多態性條帶數7～22條，平均15.4條。從以上SCoT標記的擴增多態性數據可看出供試的4種黃連屬藥用植物遺傳多樣性極為豐富。

2.2 種間與種內的SCoT遺傳多樣性 POPGENE計算的各種遺傳參數表明黃連屬4種藥用植物種間的遺傳多樣性非常豐富，見表3。種間的PPB=99.1%，N_a=1.990 6，N_e=1.329 3，H=0.212 2，I=0.337 8。但種內的遺傳多樣性極低，且變化幅度不大：PPB在10.4%～24.9%，平均為16.8%；N_a變化幅度在1.103 5～1.249 4，平均為1.168 2；N_e變化幅度在1.041 3～1.093 4，平均為1.073 0；H變化幅度在0.026 0～0.057 4，平均為0.043 7；I變化幅度在0.042 1～0.090 5，平均為0.067 7）。

2.3 種群遺傳結構分析 在假設遺傳平衡條件下計算出供試材料總基因多樣性H_t=0.212 2，其中種內基因多樣性H_s=0.043 7，種間的基因多樣度性（D_st=H_t-H_s）為0.168 5，大于H_s；Nei′s的基因分化系數G_st=0.794 0，表明有79.4%的遺傳變異存在于種間，20.6%的遺傳變異存在于種內；基因流[N_m=0.5（1-G_st）/G_st]為0.129 7，表明種間基因交流處于極低等水平，交流極少。以上POPGENE計算出的遺傳參數表明：供試黃連屬藥用植物種間的遺傳多樣性遠大于種內的遺傳多樣性，種間存在著高度的遺傳分化。

2.4 親緣關系與聚類分析 POPGENE計算出供試黃連屬藥用植物種間遺傳一致度在0.730 8～0.792 9，平均為0.765 0，遺傳距離分布在0.232 1～0.313 6，平均為0.268 4，見表4。黃連、三角葉黃連與峨眉黃連3種黃連屬藥用植物之間的遺傳一致度、遺傳距離差異不大，表明三者之間的親緣關緣較近。日本黃連與其他3種黃連屬藥用植物的遺傳一致度較低，遺傳距離較遠，表明親緣關系較遠。

用TREECONW軟件計算供試材料個體間的遺傳距離，按UPGMA法進行聚類分析，建立聚類分支樹狀圖，見圖2。從聚類圖中可看出：峨眉黃連與三角葉黃連首先聚為一支，再與黃連聚在一起；日本黃連處于單獨一個進化枝上，與其他分子界限明顯。個體聚類的置信值均為100，可信度極高，均聚在同種藥用植物內。

3 討論

3.1 基于分子標記的黃連屬藥用植物遺傳多樣性 從已有文獻報道看，揭示黃連屬藥用植物遺傳多樣性主要采用RAPD標記，DNA條形碼多用于探索黃連屬藥用植物的鑒別。常艷波^[6]采用RAPD標記揭示黃連屬藥用植物（黃連、三角葉黃連、云南黃連、峨眉黃連、線萼黃連、西藏黃連、日本黃連）在物種水平上多態位點比率（PPB）為97.46%，同一種C. chinesis居群間的多態位點比率為10.24%，顯著低于物種間的多態位點比率，表明黃連屬植物的遺傳多樣性是豐富的。武敬亮^[7]用RAPD標記揭示黃連屬藥用價值較大的4個物種（黃連、三角葉黃連、云連、峨眉黃連）之間多態性比例65.3%，說明各物種具有較為豐富的遺傳多樣性。韓凡^[8]用RAPD標記揭示4種黃連屬藥用植物（黃連、三角葉黃連、峨眉黃連和線萼黃連）的多態百分率為89.37%，說明黃連屬樣本之間遺傳多樣性較為豐富。本文的研究結果表明供試的4種黃連屬藥用植物物種水平上的遺傳多樣性極為豐富（多態性百分率達99.1%），但同一種內個體之間的遺傳多樣性均較低（多態性百分率10.4%～24.9%，平均僅為16.8%），其中，峨眉黃連、三角葉黃連的遺傳多樣性水平較低（多態性百分率分別為10.4%，12.5%），黃連、日本黃連的遺傳多樣性水平相對較高（19.5%，24.9%）。由于供試材料、分子標記的不同，所得到的結果有一定差異，但均從不同角度揭示了黃連屬藥用植物較為豐富的遺傳多樣性，不同種之間存在特有的遺傳特性。

一個物種遺傳多樣性的高低與其分布范圍、繁育方式等有關。從物種的分布范圍看：黃連主要分布在長江中下游各省市，是黃連屬藥用植物中用藥最廣泛、栽培面積最廣、產量最大的物種；峨眉黃連常生長在陸峭的巖縫間，主要分布在四川峨眉山脈，居群小而分散，屬于狹域種；三角葉黃連僅分布于僅分布于四川峨眉山、洪雅一帶，多為栽培，生長區域狹窄；日本黃連為從日本引進物種。從以上分析推斷，黃連屬藥用植物遺傳多樣性水平與物種的分布區域及小生境有著一定的聯系。從繁育方式上看，三角葉黃連染色體為三倍體（2n=3x=27），偶爾有種子結實，但其大多不育，生產上采用“扯秧苗栽秧”的無性繁殖方式進行繁殖^[3]，因此其遺傳多樣性水平低于種子繁殖物種。

3.2 黃連屬常用藥用植物的親緣關系研究 中藥材黃連的來源植物除《中國藥典》收載的黃連、云連、三角葉黃連之外，我國各地尚有多種黃連屬植物如峨眉黃連、西藏黃連、短萼黃連等的根莖作為黃連入藥。對于某些種的分類地位無論是形態標記還是分子標記得到的結果是一致的，如云連，但某些種的分類地位頗有爭議，爭議主要集中在黃連、三角葉黃連與峨眉野連三者的親緣關系。研究者采用不同方法，多角度對黃連屬藥用植物的親緣關系進行研究，早期多采用形態標記，后期引入分子標記，但不同的研究方法得到的結果并非完全一致。形態分類主要以花瓣形狀、萼片的形狀及寬窄、葉形差異為黃連屬植物分類的依據，認為三角葉黃連與黃連親緣關系更近^[1]。武敬亮采用RAPD分子標記，也認為三角葉黃連與黃連親緣關系更近。但韓凡等^[8]采用RAPD分子標記與本文采用的SCoT標記，認為三角葉黃連與峨眉黃連親緣關系更近。

此外，對三角葉黃連的栽培起源爭議也頗大，一種認為三角葉黃連為黃連與峨眉野連的雜交種，另一種認為三角葉黃連是峨眉黃連種內多倍化的結果。前者的主要證據為劉曉光^[9]采用DNA條形碼研究發現三角葉黃連分別與黃連和峨眉黃連具有較多共有基因型，由此可初步推斷三角葉黃連是由黃連和峨眉黃連雜交起源而來。后者主要證據源于黃驥等^[3]對國產黃連屬植物的染色體核型分析，表明三角葉黃連為三倍體，對稱程度較高，核型為3A型，其染色體大小與峨眉黃連最接近，結合這2個種分布區相重疊的實際，可以推測三角葉黃連是峨眉黃連種內多倍化并以營養繁殖方式固定下來，即三角葉黃連應栽培起源于峨眉黃連。本文的聚類結果支持后一種推測。

[參考文獻]

[1]王天志，方忻平，張浩，等. 藥用黃連的原植物研究[J].中藥通報，1988，13（1）：7.

[2]孫紅祥，張昌禧，吳遠文. 黃連屬藥用植物的數量分類學研究[J].中國藥學雜志，1995，30（1）：7.

[3]黃驥，陳浙，夏志華，等.國產黃連屬植物的染色體核型分析[J].西北植物學報，2013，33（5）：931.

[4]向剛，傅體華，范巧佳.中藥黃連3種基原植物的核型研究[J].中國中藥雜志，2010， 35 （11）：1367.

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

[6]常艷波. 黃連屬藥用植物RAPD鑒定與親緣關系研究[D].成都：四川大學，2005.

[7]武敬亮. 4種黃連屬（Coptis）植物的RAPD鑒定以及黃連的居群多樣性分析[D].南充：西華師范大學，2005.

[8]韓凡，鄒嘉欣，宋良科，等. 峨眉黃連及其近緣種的遺傳多樣性分析[J].安徽農業科學，2012，40（30）：14648.

[9]劉曉光. 基于群體遺傳學的黃連屬DNA條形碼研究[D].濟南：山東中醫藥大學，2012.

[責任編輯 呂冬梅]