陳宣, 云勇, 吳宇佳, 戚華沙, 楊立榮, 陳加利, 鄭道君*

海南島油茶種質(zhì)資源遺傳多樣性的SRAP分析

陳宣1, 云勇1, 吳宇佳2, 戚華沙1, 楊立榮1, 陳加利1, 鄭道君1*

(1. 海南省農(nóng)業(yè)科學(xué)院熱帶園藝研究所, 海南省熱帶特種經(jīng)濟(jì)植物種質(zhì)資源創(chuàng)新利用重點(diǎn)實(shí)驗(yàn)室，?？?571100; 2. 海南省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與土壤研究所, ?？?571100)

為了解海南島油茶()種質(zhì)資源的遺傳多樣性，采用SRAP分子標(biāo)記，對(duì)海南島油茶主要分布區(qū)的31個(gè)居群進(jìn)行了遺傳多樣性和親緣關(guān)系分析。結(jié)果表明，海南島油茶資源的遺傳多樣性低，物種水平的多態(tài)性百分率(PPB)為98.30%，Nei’s基因多樣性()為0.222 8，Shannon信息指數(shù)()為0.353 8；居群水平的PPB=40.96%, 觀(guān)測(cè)等位基因數(shù)(N)為1.409 6, 有效等位基因數(shù)(N)為1.237 1,=0.138 5,=0.208 3，這與海南島油茶豐富的表型多樣性水平不一致。海南島油茶資源遺傳分化較大，居群間基因交流有限，不同居群間的遺傳分化指數(shù)(G)為0.380，基因流(N)為0.813 9<1。遺傳變異主要發(fā)生在居群內(nèi)，有38.05%的變異存在居群間，61.95%存在于居群內(nèi)。遺傳距離為0.022 6~0.276 4，平均為0.107 7，居群間的親緣關(guān)系較近。UPGMA聚類(lèi)分析表明，在遺傳距離為0.11處，可將31個(gè)油茶居群聚為6類(lèi)，表現(xiàn)出明顯的行政區(qū)域性，而與地理距離關(guān)系不大。因此，海南島油茶資源遺傳多樣性低，親緣關(guān)系近可能導(dǎo)致自交或近交不親和，可能是海南油茶林分花量大而結(jié)實(shí)低的主要內(nèi)在原因。

油茶；SRAP標(biāo)記；遺傳多樣性；居群；海南島

油茶主要是指山茶科(Theaceae)山茶屬()植物中油脂含量較高且具有栽培經(jīng)濟(jì)價(jià)值的一類(lèi)植物的總稱(chēng)，是我國(guó)重要的木本食用油料樹(shù)種。油茶在海南俗稱(chēng)“山柚”。海南島為中國(guó)油茶資源分布的最南緣，各地均有野生種和栽培種，具有種植和利用油茶的悠久歷史[1]。海南島氣候獨(dú)特，海峽隔離，油茶又是異花授粉植物，在環(huán)境和遺傳控制的雙重作用下，海南島油茶資源表型變異極為豐富且獨(dú)特[2-4]。據(jù)海南省林業(yè)局統(tǒng)計(jì)，2015年底，全省油茶種植面積達(dá)到4 100 hm2，主要分布在瓊海、屯昌、澄邁、五指山、瓊中、定安、海口等市縣，2015年茶籽產(chǎn)量486.3 t，山柚油產(chǎn)量超過(guò)80 t，約占全國(guó)產(chǎn)量的0.1%。至2018年底，海南省已認(rèn)定油茶優(yōu)良品種20個(gè)，但尚未建立規(guī)范的良種采穗圃。海南島具有豐富而獨(dú)特的油茶資源、良好的水熱條件和產(chǎn)業(yè)發(fā)展勢(shì)頭，近年來(lái)海南島油茶資源的研究與開(kāi)發(fā)利用受到普遍關(guān)注，急需對(duì)海南島油茶資源開(kāi)展科學(xué)系統(tǒng)評(píng)價(jià)。

植物種質(zhì)資源鑒定與評(píng)價(jià)是資源有效保護(hù)與創(chuàng)新利用的前提。分子標(biāo)記技術(shù)不受季節(jié)和環(huán)境的限制，穩(wěn)定性高，并且分子標(biāo)記數(shù)量極多，遍及整個(gè)基因組，多態(tài)性高，是植物資源鑒定與評(píng)價(jià)的重要手段。其中，RAPD、ISSR和AFLP標(biāo)記技術(shù)已成功應(yīng)用于油茶資源的遺傳多樣性評(píng)價(jià)[5-7]、品種鑒定與分類(lèi)[8-10]、以及抗病性、果油率和產(chǎn)油量與分子標(biāo)記的相關(guān)性[11]，為油茶育種親本的選配和品種鑒定提供科學(xué)理論依據(jù)。相關(guān)序列擴(kuò)增多態(tài)性(sequence-related amplified polymorphism, SRAP)技術(shù)是利用特定引物對(duì)ORFs區(qū)域(open reading frames)進(jìn)行擴(kuò)增，具有易于操作、花費(fèi)少、多態(tài)性強(qiáng)及標(biāo)記分布均勻等特點(diǎn)[12]，已廣泛應(yīng)用于植物遺傳多樣性研究[13-14]。前人[15-17]應(yīng)用SRAP標(biāo)記技術(shù)分別對(duì)四川省名邛臺(tái)地油茶資源、全國(guó)19個(gè)居群的小果油茶資源和廣東省廣寧紅花油茶種質(zhì)資源進(jìn)行了遺傳多樣性評(píng)價(jià)。

本課題組前期研究表明，海南油茶資源具有豐富的表型多樣性，這是否是其豐富遺傳多樣性的體現(xiàn)？群體間的遺傳結(jié)構(gòu)和親緣關(guān)系是怎樣的模式？在悠久的引種與栽培過(guò)程中基因流是否通暢？本研究利用SRAP分子標(biāo)記技術(shù)，對(duì)海南島主要油茶分布區(qū)的代表性種質(zhì)資源進(jìn)行遺傳多樣性分析，探討其遺傳分化和親緣關(guān)系，為海南油茶分子標(biāo)記輔助育種，充分挖掘海南油茶優(yōu)良種質(zhì)資源和研究其引種栽培歷史提供科學(xué)依據(jù)。

1 材料和方法

1.1 材料

試驗(yàn)材料來(lái)自海南省海口、定安、文昌、瓊海、澄邁、屯昌、瓊中、五指山和臨高等9個(gè)市縣27個(gè)鄉(xiāng)鎮(zhèn)(國(guó)營(yíng)農(nóng)場(chǎng))的31個(gè)油茶居群，基本涵蓋了海南省的主要油茶分布區(qū)。居群信息見(jiàn)表1。

1.2 方法

總DNA的提取 采用彭邵鋒[8]的改良CTAB法。

SRAP-PCR擴(kuò)增 反應(yīng)在IcyclerTMThermal Cycler型PCR儀上進(jìn)行。SRAP-PCR反應(yīng)條件和反應(yīng)程序按鄭道君等[13]的方法，引物由上海英駿生物技術(shù)有限公司合成，每個(gè)引物對(duì)供試材料重復(fù)擴(kuò)增2次。取10L擴(kuò)增產(chǎn)物，采用2.0%瓊脂糖凝膠電泳檢測(cè)SRAP-PCR的產(chǎn)物，GoodView核酸染色, 在Gel Doc XR型凝膠成像分析系統(tǒng)下照相并記錄。以3S 100 bp DNA ladder作為對(duì)照分子量標(biāo)準(zhǔn)。

1.3 數(shù)據(jù)的統(tǒng)計(jì)分析

根據(jù)marker確定PCR條帶分子量，選擇清晰的擴(kuò)增譜帶用0-1二元矩陣表示，然后統(tǒng)計(jì)分析多態(tài)性條帶百分率(PPB)=多態(tài)位點(diǎn)數(shù)/位點(diǎn)總數(shù)。將0-1矩陣在NTSYSpc (2.1版)軟件中處理, 應(yīng)用SIMQUAL法計(jì)算Jaccard遺傳相似性系數(shù)(genetic similarity, GS)，GS=2N/(N+N)，式中,N指兩個(gè)個(gè)體共有的帶數(shù)，N+N指兩個(gè)個(gè)體所有帶數(shù)之和。采用SAHN程序和UPGMA (unweighted pair group mean average)方法進(jìn)行聚類(lèi)分析，通過(guò)Tree plot模塊構(gòu)建聚類(lèi)圖。應(yīng)用POPGENE 1.31軟件在假定種群處于Hardy-Weinberg平衡狀態(tài)下進(jìn)行遺傳參數(shù)分析, 分別計(jì)算觀(guān)測(cè)等位基因數(shù)(observed number of alleles,N)、有效等位基因數(shù)(effective number of alleles,N)、Nei’s基因多樣性指數(shù)(Nei’s gene diversity,)、Shannon’s信息指數(shù)(Shannon’s infor- mation index,)，并對(duì)不同居群間的遺傳分化系數(shù)(G)和基因流進(jìn)行分析。為了明確這些多樣性指數(shù)在各居群間的變異程度，對(duì)所有居群的各參數(shù)進(jìn)行變異系數(shù)(, %)分析。

表1 海南油茶居群

續(xù)表(Continued)

2 結(jié)果和分析

2.1 SRAP-PCR擴(kuò)增多態(tài)性

選取油茶來(lái)源地域及形態(tài)不同的5份種質(zhì)材料，對(duì)經(jīng)退火溫度初選的100對(duì)SRAP引物進(jìn)行篩選，共選出9對(duì)有效引物，然后分別對(duì)31個(gè)居群材料進(jìn)行SRAP-PCR擴(kuò)增，共擴(kuò)增出122條清晰的譜帶，每對(duì)引物平均擴(kuò)增13.56條，其中多態(tài)性條帶120條，擴(kuò)增的多態(tài)性條帶百分率(PPB)高達(dá)98.30% (表2)，多數(shù)條帶為200~2 500 bp。每對(duì)引物對(duì)不同材料擴(kuò)增的SRAP-PCR產(chǎn)物有一定差異,反映出不同地域與不同居群的種質(zhì)材料之間存在遺傳多樣性以及不同程度的親緣關(guān)系。

表2 9對(duì)SRAP引物對(duì)31個(gè)海南油茶居群種質(zhì)的擴(kuò)增

2.2 遺傳多樣性分析

根據(jù)SRAP-PCR譜帶統(tǒng)計(jì)結(jié)果，利用PopGene軟件分析31個(gè)居群的遺傳多樣性(表3)。在物種水平上，PPB為98.30%，N為1.975 4，N為1.352 8，為0.222 8，為0.353 8，可見(jiàn)海南島油茶資源具有較高的遺傳多樣性水平。但是，在居群水平上的遺傳變異卻比較低，PPB為13.93%~59.84%，平均40.96%；N為1.139 3~1.598 4，平均1.409；N為1.083 4~1.357 9，平均1.237；為0.14 8~0.307 7, 平均0.208 3；為0.107 7~0.193 5，平均0.138 5。31個(gè)居群間的遺傳多樣性差異較明顯，、和PPB的變異系數(shù)()分別為23.82%、23.67%和24.26%, 其中sjgx居群的遺傳多樣性水平最低，zxlw居群的遺傳多樣性水平最高。

表3 海南島油茶居群的遺傳多樣性指數(shù)

2.3 居群的遺傳結(jié)構(gòu)和遺傳變異分析

對(duì)31個(gè)居群進(jìn)行基因分化分析，計(jì)算不同居群間的遺傳分化系數(shù)(G)，結(jié)果表明，居群總的遺傳多樣性(H)為0.223 6±0.028 2，居群內(nèi)的基因多樣性(H)為0.138 5±0.010 3，基因流(N)僅為0.813 9<1，表明海南島油茶遺傳多樣性水平較低，G為0.380 5,可見(jiàn)居群間具有較大的遺傳分化，同時(shí)也表明有38.05%的遺傳變異存在居群間，而61.95%的遺傳變異存在于居群內(nèi)，其遺傳變異以居群內(nèi)變異為主。

2.4 聚類(lèi)分析

根據(jù)SRAP擴(kuò)增條帶分析結(jié)果，采用UPGMA法對(duì)31個(gè)不同居群進(jìn)行聚類(lèi)分析(圖1)。不同居群間的遺傳距離為0.022 6~0.276 4，平均遺傳距離為0.107 7。其中l(wèi)kfh與lkly遺傳關(guān)系最近，遺傳距離僅0.022 6，均來(lái)自定安嶺口鎮(zhèn)；sjgx與tnlf雖然地理距離不遠(yuǎn)，但遺傳距離卻最遠(yuǎn)，為0.276 4。在遺傳距離為0.11處，可將31個(gè)居群明顯聚為6類(lèi),其中類(lèi)群A由來(lái)自瓊海的所有居群sbwl、yjmp、hszj、lj、zyxc和海口的部分居群jzcc、ssjx、yxrz、sjqm組成；B類(lèi)群包括澄邁的所有居群sfbl、wral、zxlw、jllm，臨高的所有居群htww、hsct、nbgj和定安的所有居群lkfh和lkly；C類(lèi)群由瓊中的所有居群wlby、czxz、nknl，屯昌的hlsy，以及五指山的cthw和csfs組成；D類(lèi)群僅由文昌的2個(gè)居群jslj和tnlf組成；E類(lèi)包括五指山的smxc和?？诘膉zbs; F類(lèi)是?？诘膕jgx和hq。除了來(lái)自?？诘木尤汉臀逯干骄尤翰荒芨髯跃蹫橐活?lèi)外，其他市縣來(lái)源的居群均按行政區(qū)域聚為一類(lèi)，表現(xiàn)出明顯的區(qū)域性。

圖1 31個(gè)海南油茶居群間的聚類(lèi)分析圖

3 討論

3.1 海南島油茶資源的遺傳多樣性

遺傳多樣性是物種長(zhǎng)期生存和保持進(jìn)化的基礎(chǔ)，同時(shí)也植物良種培育的遺傳物質(zhì)基礎(chǔ)。PPB、、作為評(píng)價(jià)種群內(nèi)和種群間遺傳多樣性的指標(biāo), 其數(shù)值越大，表明種群的遺傳多樣性越高[18]。SRAP、RAPD和ISSR一樣，均為顯性標(biāo)記，研究結(jié)果具有可比性。Nybom[19]統(tǒng)計(jì)認(rèn)為植物的平均為0.22 (基于RAPD和ISSR的平均值)，其中37種多年生植物的為0.25。比較而言，海南油茶資源物種水平上的值(0.222 8)僅與上述物種的平均值相近，而居群水平的平均H僅為0.138 5 (0.048 6~0.207)。此外，黃勇等[16]利用SRAP分子標(biāo)記對(duì)貴州、湖南、江西等8個(gè)小果油茶分布區(qū)19個(gè)居群進(jìn)行遺傳多樣性分析，認(rèn)為小果油茶具有豐富的遺傳多樣性, 19個(gè)居群的、、PPB分別為0.321 3~0.388 7、0.473 2~ 0.567 6和84.96%~95.58%；物種水平的、、PPB分別為0.422 3、0.612 6和100%。謝薈等[17]利用SRAP分子標(biāo)記建立了廣東省廣寧紅花油茶()種質(zhì)資源4個(gè)種群的遺傳多樣性評(píng)價(jià)體系,結(jié)果表明，種群的、、PPB分別為0.217 9~0.356 3、0.321 7~0.524 8和58.25%~92.72%，平均值分別為0.307 8、0.456 3和83.62%。張婷等[20]采用SRAP標(biāo)記對(duì)湖北省油茶5個(gè)居群的遺傳多樣性進(jìn)行了研究，結(jié)果表明，居群的平均N、、分別為1.516 1、0.274 8和0.389 9。因此，本研究中海南島油茶資源表現(xiàn)出較低的遺傳多樣性水平，物種水平的和僅分別為0.222 8和0.353 8；居群的遺傳變異較物種水平更低，和僅為0.048 6~0.207和0.073~ 0.288 3，這與海南油茶在表型方面豐富的多樣性水平不一致。楊立榮等[21]對(duì)海南島21個(gè)老油茶林果實(shí)和種子的10個(gè)數(shù)量性狀指標(biāo)進(jìn)行了分析，結(jié)果表明海南油茶種實(shí)表型多樣性十分豐富, 10個(gè)種實(shí)數(shù)量性狀的變異系數(shù)()為7.87%~52.93%，平均變異系數(shù)為29.36%，變異程度均高于小果油茶()、浙江紅花油茶()、滇西騰沖紅花油茶()、廣東省的廣寧紅花油茶和越南油茶()；平均Shannon- Wiener多樣性指數(shù)為2.749 9，接近或大于小果油茶(主要分布于7個(gè)省份)。可見(jiàn)，海南油茶低水平的遺傳變異不能解析豐富的表型變異。

油茶組是山茶屬中多倍體出現(xiàn)頻率較高的類(lèi)群，有2、3、4、5、6、7和8，且核型分析存在多樣性，其中奇數(shù)多倍體為栽培品種，包括油茶()、小果油茶、越南油茶和攸縣油茶()均是多倍體[22-24]，這些物種除小果油茶外其余物種均存在種內(nèi)多倍性現(xiàn)象，如越南油茶核型2=120=86m+26sm+8st+2sat[25], 或2=105=61m+37sm+7st+2sat，或2=120=86m+26sm+8st+ 2sat[26]；普通油茶核型為2=30=24m+4sm+2st+ 2sat, 或2=90=60m+29sm+1st+3sat，或2=90= 60m+22sm+8s[27]。倍性變異會(huì)影響個(gè)體、花朵和果實(shí)大小等農(nóng)藝性狀[28]。海南油茶的核型多樣性是否豐富？是否在海南島油茶資源豐富的表型多樣性過(guò)程中扮演重要角色？從本文中海南島油茶資源較低水平的遺傳多樣性結(jié)果推測(cè)，核型的變異或許不能使海南島油茶資源間的遺傳背景或遺傳多樣性變得豐富。目前，本課題組正在進(jìn)行海南油茶資源的核型多樣性分析，并研究其與形態(tài)多樣性和遺傳多樣性的相關(guān)性，對(duì)上述推測(cè)進(jìn)行驗(yàn)證。

油茶在海南島有悠久的栽培歷史，現(xiàn)存的油茶居群均是人工栽種的半野生居群，海南島油茶資源表現(xiàn)出物種水平高于居群水平的遺傳多樣性，表明半野生海南油茶居群種源來(lái)源比較復(fù)雜，而同一居群種源來(lái)源相對(duì)穩(wěn)定。同時(shí)居群水平存在一定的遺傳變異，究其原因是這些居群均是以實(shí)生苗建成, 油茶是較為嚴(yán)格的蟲(chóng)媒異花授粉植物，F(xiàn)1代實(shí)生苗分離的結(jié)果。

3.2 海南島油茶群體遺傳結(jié)構(gòu)與親緣關(guān)系

油茶是較為嚴(yán)格的蟲(chóng)媒異花授粉植物。Nybom[19]對(duì)31種雜交植物進(jìn)行統(tǒng)計(jì)，平均G為0.22；張德全等[29]的統(tǒng)計(jì)表明，蟲(chóng)媒植物的平均G為0.264。本研究中G為0.380 5，高于上述平均值，也比我國(guó)小果油茶遺傳分化值要高(=0.134 3)[16]，可見(jiàn)海南島油茶居群間具有較大的遺傳分化。此外，在整體水平上居群間的基因流為0.813 9，說(shuō)明海南島油茶資源不同居群間存在較小的基因交流。遺傳分化系數(shù)顯示，海南油茶總的遺傳變異中有38.05%的變異存在于居群間，有61.95%的遺傳變異存在于居群內(nèi)；在小果油茶總的遺傳變異中，居群內(nèi)占86.58%，而居群間僅占13.42%[16]；Kaundun[30]利用RAPD-PCR標(biāo)記對(duì)韓國(guó)、日本和臺(tái)灣27個(gè)優(yōu)良茶樹(shù)品種的多樣性進(jìn)行了研究，顯示71%的變異位于組內(nèi)，29%位于組間。本研究結(jié)論與其一致，其遺傳變異均以居群內(nèi)變異為主。海南油茶資源31個(gè)居群間的遺傳距離為0.022 6~0.276 4，平均為0.107 7, 顯示較近的親緣關(guān)系，僅略高于廣布于我國(guó)南方的小果油茶(19個(gè)居群的遺傳距離為0.021 6~0.164 5，平均為0.098 8)[16]。聚類(lèi)分析表明，在遺傳距離為0.11處，可將31個(gè)海南油茶居群聚為6類(lèi)，聚類(lèi)結(jié)果與行政區(qū)域有一定的相關(guān)性，但與各居群的地理距離不相關(guān)，除了來(lái)自海口的居群和五指山居群不能各自聚為一類(lèi)外，其他市縣來(lái)源的居群均按行政區(qū)域聚為一類(lèi)，表現(xiàn)出明顯的行政區(qū)域性。這與劉揚(yáng)[31]應(yīng)用ISSR標(biāo)記對(duì)來(lái)自海南、廣西、江西的55份油茶資源進(jìn)行多態(tài)性和聚類(lèi)分析的結(jié)果基本一致，也同小果油茶的聚類(lèi)結(jié)果相一致[16]?？梢?jiàn),海南島油茶親緣關(guān)系較近，且表現(xiàn)出較強(qiáng)的行政區(qū)域性。我們推測(cè)，海南油茶資源居群的聚類(lèi)關(guān)系可能與海南油茶悠久的引種栽培歷史, 以及計(jì)劃經(jīng)濟(jì)體制下各行政區(qū)域的油茶推廣方式密切相關(guān)。

油茶為異花授粉植物，自交或近交不親和，結(jié)實(shí)率極低[32]。我們的調(diào)查表明，現(xiàn)有海南島油茶林分個(gè)體的開(kāi)花率高，可達(dá)100%，花量也極大，但結(jié)果個(gè)體僅為20%左右，甚至更低，以致產(chǎn)量極低，這種現(xiàn)象在近5~10年栽植的油茶林分中也表現(xiàn)得更為突出[33]，極大地影響了農(nóng)民種植油茶的信心和積極性。本研究表明海南島油茶資源親緣關(guān)系近,依據(jù)油茶自交不結(jié)實(shí)或結(jié)實(shí)率低的特性，因此，我們認(rèn)為海南島油茶資源親緣關(guān)系近可能造成自交或近交不親和，是致使海南油茶林分花量大而結(jié)實(shí)低的主要內(nèi)在原因之一。本課題組下一步將從小氣候?qū)ǚ刍钚院蛡鞣壅咝袨榈挠绊憽⑼寥鲤B(yǎng)分供給對(duì)結(jié)實(shí)的影響，以及林分個(gè)體間的遺傳親緣關(guān)系等方面進(jìn)行分析和論證。

[1] Zhou Y D. On history and culture of[J]. Trop For, 2017, 45(2): 42-48. doi: 10.3969/j.issn.1672-0938.2017.02.013.周亞?wèn)|. 淺談?dòng)筒栉幕?[J]. 熱帶林業(yè), 2017, 45(2): 42-48. doi: 10.3969/j.issn.1672-0938.2017.02.013.

[2] Yuan J, Han Z Q, He S Y, et al. Investigation and cluster analysis of main morphological and economical characters for oiltea resource in Hainan Province [J]. J Plant Genet Resour, 2014, 15(6): 1380-1384. doi: 10.13430/j.cnki.jpgr.2014.06.031.袁軍, 韓志強(qiáng), 賀舍予, 等. 海南省油茶資源主要形態(tài)和經(jīng)濟(jì)性狀及聚類(lèi)分析 [J]. 植物遺傳資源學(xué)報(bào), 2014, 15(6): 1380-1384. doi: 10.13430/j.cnki.jpgr.2014.06.031.

[3] ZHENG D J, PAN X Z, ZHANG D M, et al. Survey and analysis on tea-oil camellia resource in Hainan[J]. J NW For Univ, 2016, 31(1): 130-135. doi: 10.3969/j.issn.1001-7461.2016.01.24.鄭道君, 潘孝忠, 張冬明, 等. 海南油茶資源調(diào)查與分析 [J]. 西北林學(xué)院學(xué)報(bào), 2016, 31(1): 130-135. doi: 10.3969/j.issn.1001-7461. 2016.01.24.

[4] CHEN Y Z, CHEN L S, LI R F, et al. Survey and the industry develop- ment proposals ofsp. resources in Hainan [J]. Trop For, 2017, 45(1): 49-52. doi: 10.3969/j.issn.1672-0938.2017.01.015.陳永忠, 陳隆升, 李儒法, 等. 海南油茶資源調(diào)研及產(chǎn)業(yè)發(fā)展建議 [J]. 熱帶林業(yè), 2017, 45(1): 49-52. doi: 10.3969/j.issn.1672-0938. 2017.01.015.

[5] HUANG Y F, CHEN X M, ZHUANG X Y, et al. Analysis of genetic diversity ingermplasms [J]. Sci Silv Sin, 2006, 42(4): 38-43. doi: 10.3321/j.issn:1001-7488.2006.04.007.黃永芳, 陳錫沐, 莊雪影, 等. 油茶種質(zhì)資源遺傳多樣性分析 [J]. 林業(yè)科學(xué), 2006, 42(4): 38-43. doi: 10.3321/j.issn:1001-7488.2006.04. 007.

[6] YU X Y, YU F Y, LIU J B, et al. Identification and genetic diversity analysis ofvarieties using ISSR marker [J]. J Nan- jing For Univ (Nat Sci), 2013, 37(1): 61-66. 于小玉, 喻方圓, 劉建兵, 等. ISSR在油茶品種鑒別和遺傳多樣性分析中的應(yīng)用 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 37(1): 61-66.

[7] ZHANG T, LIU S Q, MEI H, et al. AFLP analysis on the genetic diversity offrom different regions of Hubei Province [J]. J Anhui Agric Sci, 2011, 39(23): 14070-14071,14075. doi: 10. 3969/j.issn.0517-6611.2011.23.067.張婷, 劉雙青, 梅輝, 等. 湖北省不同地區(qū)油茶遺傳多樣性的AFLP分析 [J]. 安徽農(nóng)業(yè)科學(xué), 2011, 39(23): 14070-14071,14075. doi: 10. 3969/j.issn.0517-6611.2011.23.067.

[8] PENG S F, ZHANG D Q, CHEN Y Z, et al. Genetic diversity analysis of elite cultivars ofby SRAP [J]. J CS Univ For Technol, 2011, 31(1): 80-85. doi: 10.3969/j.issn.1673-923X.2011.01.013.彭邵鋒, 張黨權(quán), 陳永忠, 等. 14個(gè)油茶良種遺傳多樣性的SRAP分析 [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2011, 31(1): 80-85. doi: 10.3969/ j.issn.1673-923X.2011.01.013.

[9] ZHANG G W, ZHONG W B, WUYunTaNa, et al. Identification of oil tea () superior clones by ISSR molecular marker [J]. For Res, 2007, 20(2): 278-282. doi: 10.3321/j.issn:1001-1498. 2007.02.024.張國(guó)武, 鐘文斌, 烏云塔娜, 等. 油茶優(yōu)良無(wú)性系ISSR分子鑒別 [J]. 林業(yè)科學(xué)研究, 2007, 20(2): 278-282. doi: 10.3321/j.issn:1001- 1498.2007.02.024.

[10] WU Y Y. Analysis of genetic diversity and identification ofgermplasms [D]. Nanjing: Nanjing Forestry University, 2011: 36-40. 吳鶯鶯. 油茶種質(zhì)資源遺傳多樣性分析及品種鑒定 [D]. 南京: 南京林業(yè)大學(xué), 2011: 36-40.

[11] LI L. Analysis on genetic diversity and quality ofresources in Hubei Province [D]. Wuhan: Huazhong Agricultural University, 2010: 18-34. 李良. 湖北省油茶優(yōu)良單株遺傳多樣性及品質(zhì)初步分析 [D]. 武漢: 華中農(nóng)業(yè)大學(xué), 2010: 18-34.

[12] 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 in[J]. Theor Appl Genet, 2001, 103: 455-461. doi: 10.1007/s001220100570.

[13] ZHENG D J, YUN T H, ZHANG Z L, et al. Study on genetic diversity and relationship for the Hainan island landraces of[J]. J Nucl Agric Sci, 2016, 30(5): 869-877. doi: 10.11869/j.issn.100- 8551.2016.05.0869.鄭道君, 云天海, 張治禮, 等. 中國(guó)南瓜海南農(nóng)家品種資源遺傳多樣性和親緣關(guān)系研究 [J]. 核農(nóng)學(xué)報(bào), 2016, 30(5): 869-877. doi: 10.11869/j.issn.100-8551.2016.05.0869.

[14] QI L, WANG X N, ZHANG J Z, et al. Study on genetic diversity and differentiation of three wild rice species in Hainan using SRAP markers [J]. J Plant Genet Resour, 2013, 14(3): 402-406. 齊蘭, 王效寧, 張吉貞, 等. 利用SRAP標(biāo)記研究海南野生稻的遺傳多樣性與遺傳分化 [J]. 植物遺傳資源學(xué)報(bào), 2013, 14(3): 402- 406.

[15] ZHOU L J, KAO A D, LIU J, et al. Genetic diversity analysis of wild’s germplasm resources based on SRAP markers [J]. Mol Plant Breed, 2017, 15(10): 4242-4247. 周莉君, 考安都, 劉靜, 等. 基于SRAP標(biāo)記的野生油茶種質(zhì)遺傳多樣性分析 [J]. 分子植物育種, 2017, 15(10): 4242-4247.

[16] HUANG Y. Analysis ofgenetic diversity based on SRAP markers [J]. Sci Silv Sin, 2013, 49(3): 43-50. doi: 10.11707/j. 1001-7488.20130306.黃勇. 基于SRAP分子標(biāo)記的小果油茶遺傳多樣性分析 [J]. 林業(yè)科學(xué), 2013, 49(3): 43-50. doi: 10.11707/j.1001-7488.20130306.

[17] XIE H, ZHU W J, ZHANG Y Z, et al. Genetic diversity ofchi germplasm based on SRAP markers [J]. J CS Univ For Technol, 2017, 37(8): 93-97,113. 謝薈, 祝文娟, 張應(yīng)中, 等. 基于SRAP分子標(biāo)記的廣寧紅花油茶種質(zhì)資源遺傳多樣性研究 [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2017, 37(8): 93-97,113.

[18] WU Z X, WANG H Z, SHI N N, et al. The genetic diversity ofby ISSR [J]. Hereditas, 2008, 30(5): 627-632. doi: 10.3321/ j.issn:0253-9772.2008.05.016.吳振興, 王慧中, 施農(nóng)農(nóng), 等. 蘭屬植物ISSR遺傳多樣性分析 [J]. 遺傳, 2008, 30(5): 627-632. doi: 10.3321/j.issn:0253- 9772.2008.05.016.

[19] NYBOM H. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants [J]. Mol Ecol, 2004, 13(5): 1143-1155. doi: 10.1111/j.1365-294X.2004.02141.x.

[20] ZHANG T, LIU S Q, MEI H, et al. Genetic diversity offrom Hubei province revealed by SRAP analysis [J]. J Hua- zhong Norm Univ (Nat Sci), 2011, 45(3): 477-479,484. 張婷, 劉雙青, 梅輝, 等. 湖北省油茶資源遺傳多樣性的SRAP分析 [J]. 華中師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 45(3): 477-479,484.

[21] YANG L R, ZHANG Z L, CHEN J L, et al. The quantitative characters and diversity of oiltea fruit in Hainan Province [J]. Nonwood For Res, 2018, 36(3): 69-76. doi: 10.14067/j.cnki.1003-8981.2018.03.011.楊立榮, 張治禮, 陳加利, 等. 海南油茶種實(shí)數(shù)量性狀多樣性評(píng)價(jià) [J]. 經(jīng)濟(jì)林研究, 2018, 36(3): 69-76. doi: 10.14067/j.cnki.1003-8981. 2018.03.011.

[22] LI G T, LIANG T. Chromosome numbers and karyotypes in the genus[J]Guihaia, 1990, 10(2): 127-138. 李光濤, 梁濤. 山茶屬植物的染色體數(shù)目和核型 [J]. 廣西植物, 1990, 10(2): 127-138.

[23] GU Z J, SUN X F. A karyomorphological study of seventeen species of Chinese[J]. Acta Bot Yunnan, 1997, 19(2): 159-170. 顧志建, 孫先鳳. 山茶屬17個(gè)種的核形態(tài)學(xué)研究 [J]. 云南植物研究, 1997, 19(2): 159-170.

[24] ZHANG W J, MIN T L. A cytogeographical study of, Sect.[J]. Acta Bot Yunnan, 1998, 20(3): 321-328. 張文駒, 閔天祿. 山茶屬山茶組的細(xì)胞地理學(xué)研究 [J]. 云南植物研究, 1998, 20(3): 321-328.

[25] MO Z Q. Karyotype analysis in[J]. Guihaia, 1990, 10(1): 31-32. 莫澤乾. 越南油茶的核型分析 [J]. 廣西植物, 1990, 10(1): 31-32.

[26] XIAO T, GU Z, XIA L. A karyomophological study of ten species of Chinese[J]. La Kromosomo, 1991, 2(61): 2051-2085.

[27] LI M Q. Karyotype and banding patterns of Giemsa C-bands inand[J]. For Sci Technol, 1981 (4): 9-12. 黎麥秋. 普通油茶、板栗染色體組型和Giemsa C-帶的帶型研究 [J]. 林業(yè)科技通訊, 1981(4): 9-12.

[28] QIN S Y, RONG J, ZHANG W J, et al. Cultivation history ofand genetic resources in the Yangtze River Basin [J]. Biodiv Sci, 2018, 26(4): 384-395. doi: 10.17520/biods.2017254.秦聲遠(yuǎn), 戎俊, 張文駒, 等. 油茶栽培歷史與長(zhǎng)江流域油茶遺傳資源 [J]. 生物多樣性, 2018, 26(4): 384-395. doi: 10.17520/biods.2017254.

[29] ZHANG D Q, YANG Y P. A statistical and comparative analysis of genetic diversity detected by different molecular markers [J]. Acta Bot Yunnan, 2008, 30(2): 159-167. 張德全, 楊永平. 幾種常用分子標(biāo)記遺傳多樣性參數(shù)的統(tǒng)計(jì)分析 [J]. 云南植物研究, 2008, 30(2): 159-167.

[30] KAUNDUN S S, ZHYVOLOUP A, PARK Y G. Evaluation of the genetic diversity among elite tea (var.) accessions using RAPD markers [J]. Euphytica, 2000, 115(1): 7-16. doi: 10.1023/A:100393912.

[31] LIU Y, PENG Y, HE R, et al. Genetic diversity ofin Hainan by ISSR [J]. Mol Plant Breeding, 2016, 14(2): 517-523. 劉揚(yáng), 彭赟, 賀瑞, 等. 基于ISSR標(biāo)記的海南油茶種質(zhì)資源遺傳多樣性分析 [J]. 分子植物育種, 2016, 14(2): 517-523.

[32] DENG Y Y, YU X L, LUO Y B. The role of native bees on the reproductive success ofin Hunan Province, central south China [J]. Acta Ecol Sin, 2010, 30(16): 4427-4436.鄧園藝, 喻勛林, 羅毅波. 傳粉昆蟲(chóng)對(duì)我國(guó)中南地區(qū)油茶結(jié)實(shí)和結(jié)籽的作用 [J]. 生態(tài)學(xué)報(bào), 2010, 30(16): 4427-4436.

[33] ZHENG D J, PAN X Z, XIE L S, et al. Investigation and analysis ofindustry development status in Hainan Province [J]. Nonwood For Res, 2015, 33(1): 131-135. doi: 10.14067/j.cnki.1003- 8981.2015.01.025.鄭道君, 潘孝忠, 謝良商, 等. 海南省油茶產(chǎn)業(yè)發(fā)展現(xiàn)狀調(diào)查與分析 [J]. 經(jīng)濟(jì)林研究, 2015, 33(1): 131-135. doi: 10.14067/j.cnki.1003- 8981.2015.01.025.

Genetic Diversity Analysis ofResources Based on SRAP Markers in Hainan Island

CHEN Xuan1, YUN Yong1, WU Yu-jia2, QI Hua-sha1, YANG Li-rong1, CHEN Jia-li1, ZHENG Dao-jun1*

(1.Institute of Tropical Horticulture, Hainan Academy of Agricultural Sciences,Hainan Key Laboratory for Innovative Development and Utiliation of Tropical Special Economic Plants,Haikou 571100, China; 2. Institute of Agricultural Environment and Soil, Hainan Academy of Agricultural Sciences, Haikou 571100, China)

In order to evaluate genetic diversity ofresources in Hainan Island, the genetic diversity and relationship among 31 populations wereanalyzedby using SRAP molecular marker. The results showed that the genetic diversity ofresources in Hainan Island was low, the percentage of polymorphic bands (PPB), Nei’s genetic diversity index () and Shannon’s information index () were 98.30%, 0.222 8 and 0.353 8 at the species level, and 40.96%, 0.138 5 and 0.208 3 at the population level, respectively, which was inconsistent with their rich phenotypic diversity in Hainan. The coefficient of genetic differentiation among populations (G) was 0.380 5, and the gene flow (N) was 0.813 9 between different populations, indicating that genetic differentiation was high with limited gene flow among populations. The genetic variation mainly existed within populations (61.95%) and few among them (38.05%). The genetic distance ranged from 0.022 6 to 0.276 4, with the mean of 0.107 7 between 31 populations, which revealed the relatively close genetic relationship among the populations. UPGMA cluster analysis showed that 31 populations could be divided into 6 groups at the genetic distance of 0.11, showing obvious administrative area, but had little to do with geographical distance. So, it was suggested that the relatively close genetic relationship ofresources in Hainan Island might cause self-incompatibility, which could be the main internal cause for blooming but low yield offorests in Hainan Island.

; SRAP marker; Genetic diversity; Population; Hainan Island

10.11926/jtsb.4027

2018-12-04

2019-03-25

國(guó)家自然科學(xué)基金項(xiàng)目(31860082); 海南省省屬科研院所技術(shù)開(kāi)發(fā)專(zhuān)項(xiàng)(KYYS-2018-14); 海南省農(nóng)業(yè)科學(xué)院2018年度農(nóng)業(yè)科技創(chuàng)新專(zhuān)項(xiàng)資助

This work was supported by the National Natural Science Foundation of China (Grant No. 31860082), the Project for Technology Development of Hainan Provincial Scientific Research Institutions (Grant No. KYYS-2018-14), and the Project for Agricultural Science and Technology Innovation of Hainan Academy of Agricultural Sciences in 2018.

陳宣(1983~ )，男，碩士，副研究員。研究方向?yàn)闊釒Ы?jīng)濟(jì)作物遺傳育種。E-mail: chenxuan06@sina.com

Corresponding author. E-mail: daojunzh@163.com