朱弘, 蔡厚才, 李涌福, 陳萬(wàn)東, 陳林, 伊賢貴, 李蒙, 段一凡*, 王賢榮*

中國(guó)東部沿海水仙歸化群體的遺傳多樣性

朱弘1, 蔡厚才2, 李涌福1, 陳萬(wàn)東2, 陳林1, 伊賢貴1, 李蒙1, 段一凡1*, 王賢榮1*

(1. 南京林業(yè)大學(xué)南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心，生物與環(huán)境學(xué)院，亞熱帶森林生物多樣性保護(hù)國(guó)家林業(yè)局重點(diǎn)實(shí)驗(yàn)室，南京 210037；2. 南麂列島國(guó)家級(jí)海洋自然保護(hù)區(qū)管理局，浙江 溫州 325400)

為揭示我國(guó)東部歸化水仙(var.)的群體遺傳多樣性，利用2個(gè)葉綠體基因K和HA片段對(duì)采自滬、浙、閩的5個(gè)代表群體的49株水仙進(jìn)行了評(píng)估。結(jié)果表明，雙基因聯(lián)合序列的總長(zhǎng)為1 443 bp，共定義6個(gè)單倍型，各歸化群體的遺傳多樣性水平為DLSY>PTD>NJD=ZZ>CMD。AMOVA分析表明，群體內(nèi)變異為遺傳變異的主要來(lái)源(91.98%)，群體間的遺傳分化較低(F=0.08022)。群體物種水平上的譜系結(jié)構(gòu)不顯著(N=0.020<G=0.031;< 0.05)。Mantel檢驗(yàn)表明水仙群體間的遺傳距離與地理距離呈顯著的線(xiàn)性相關(guān)(=0.929，=0.02<0.05)。中性檢驗(yàn)和錯(cuò)配分布分析結(jié)果均暗示水仙群體背離了快速擴(kuò)張模型的假設(shè)。單倍型分布的中介網(wǎng)絡(luò)圖結(jié)合系統(tǒng)發(fā)育NJ樹(shù)均將所有群體劃分為2大分支。因此，我國(guó)東部沿海水仙歸化群體整體遺傳多樣性水平較低，各群體間遺傳分化較弱，遺傳變異主要來(lái)自群體內(nèi),物種近期未經(jīng)歷擴(kuò)張事件，可能是基因流受海島隔離、自身生物學(xué)特征、生境異質(zhì)性與及人為干擾的綜合作用影響。

水仙；中國(guó)東部沿海；歸化群體；葉綠體DNA；遺傳結(jié)構(gòu)

水仙(var.)為石蒜科(Amaryllidaceae)多年生鱗莖植物，是歐洲水仙()的變種，也是世界著名觀(guān)賞花卉和中國(guó)傳統(tǒng)十大名花之一，以其獨(dú)特的形態(tài)和香氣為世人所熟知。史籍考證與植物學(xué)證據(jù)表明該變種系地中海沿岸起源，主要?dú)w化分布于亞洲東部沿海地區(qū)或近海島嶼[1]，在我國(guó)栽培歷史悠久，具有重要的觀(guān)賞、經(jīng)濟(jì)與文化價(jià)值。隨著對(duì)水仙種質(zhì)資源數(shù)量與質(zhì)量需求的日益提高，前人已在抗性生理[2-3]、矮化機(jī)理[4]、組培擴(kuò)繁[5]、功能基因[6]等方面取得了一系列的研究。

掌握植物遺傳多樣性格局，是開(kāi)展資源保護(hù)與利用的前提基礎(chǔ)。探究植物群體遺傳結(jié)構(gòu)在時(shí)間和空間上的遺傳變異特征及分布格局，在反映物種進(jìn)化歷史、生境適應(yīng)等方面具有重要作用[7-9]。群體研究主要包括傳統(tǒng)的形態(tài)學(xué)和現(xiàn)代分子標(biāo)記2種手段,分子標(biāo)記具有不受環(huán)境或植株發(fā)育階段影響、規(guī)避測(cè)量誤差、可重復(fù)性高等優(yōu)勢(shì)。近10年來(lái)，隨著分子生物學(xué)發(fā)展，國(guó)內(nèi)外有關(guān)水仙屬遺傳多樣性的研究，已由早期的AFLP、ISSR、RAPD等標(biāo)記類(lèi)型[10-11]逐漸過(guò)渡到DNA測(cè)序。其中，基于葉綠體DNA (cpDNA)標(biāo)記和序列分析具有不受選擇影響、多態(tài)性豐富等優(yōu)點(diǎn)；同時(shí)，大多數(shù)被子植物cpDNA為母系遺傳，比雙親遺傳的nrDNA在亞種群間表現(xiàn)出更大的遺傳分化，能產(chǎn)生更清晰的種群歷史蹤跡[12]，因此，當(dāng)前cpDNA標(biāo)記已成為研究植物群體遺傳與譜系地理的主流手段之一，廣泛應(yīng)用于植物種間較低分類(lèi)階元和種內(nèi)譜系地理學(xué)研究。

在野外調(diào)查中，筆者發(fā)現(xiàn)當(dāng)前我國(guó)水仙歸化地普遍存在生境片斷化嚴(yán)重、種群規(guī)模和數(shù)量銳減等問(wèn)題，有關(guān)其在區(qū)域尺度上的遺傳多樣性現(xiàn)狀尚未報(bào)道。鑒于此，本研究采用2個(gè)水仙高變區(qū)的cpDNA標(biāo)記，首次對(duì)我國(guó)水仙集中分布的東部沿海地區(qū)的5個(gè)主要?dú)w化群體的序列變異進(jìn)行分析，旨在了解以水仙歸化群體為代表的海島植物地理隔絕狀態(tài)下的群體遺傳多樣性水平，以及分布格局的可能形成機(jī)制，以期為海島特色種質(zhì)資源的科學(xué)評(píng)價(jià)、合理保護(hù)與經(jīng)營(yíng)利用提供分子遺傳學(xué)依據(jù)。

1 材料和方法

1.1 材料

2017年冬季，個(gè)體樣本自北向南分別采自上海、浙江、福建的5個(gè)歸化水仙代表群體，覆蓋了水仙在中國(guó)東部沿海的主要分布區(qū)，共計(jì)49個(gè)樣本，記錄樣品編號(hào)、采集地點(diǎn)、經(jīng)緯度、花型等基本信息(表1), 其中平潭島(PTD)重瓣水仙為逸生的栽培品種‘玉玲瓏’。群體內(nèi)個(gè)體的采集距離至少相隔10 m以上，選擇當(dāng)年生新發(fā)的幼嫩葉片經(jīng)變色硅膠干燥處理后，置于-20℃冰箱備用。活體材料引種至南京中山植物園，憑證標(biāo)本存于南京林業(yè)大學(xué)樹(shù)木標(biāo)本館(NF)。

1.2 基因組DNA提取和PCR擴(kuò)增

植物基因組DNA采用柱式抽提純化法試劑盒(Tiangen, 北京)提取, 用1%瓊脂糖凝膠電泳和Tanon-2500全自動(dòng)數(shù)碼凝膠圖像分析儀(Tanon, 上海)檢測(cè)DNA質(zhì)量。參考已報(bào)道的水仙屬候選DNA條碼[13]進(jìn)行篩選，預(yù)試驗(yàn)結(jié)果表明，水仙cpDNA編碼區(qū)的K (F: 5?-CGATCTATTCATTCAATATTTC-3?, R: 5?-TCTAGCACACGAAAGTCGAAGT-3?)與非編碼區(qū)的H-A (F: 5?-GTTATGCATGAACGTAA- TGCTC-3?，R: 5?-CGCGCATGGTGGATTCACAAT- CC-3?)正反向引物片段在種內(nèi)群體水平表現(xiàn)出較高的變異，且測(cè)序成功率最高，故選擇上述序列進(jìn)行研究。PCR擴(kuò)增體系總體積為25L: 包含2×PCR MasterMix (TsingKe, 南京) 12.5L, 上下游引物各1.0L (GenScript, 南京)，DNA模板2.0L, 最后用ddH2O 8.5L補(bǔ)齊。K的PCR擴(kuò)增程序?yàn)?4℃預(yù)變性3 min，然后94℃變性30 s，41℃退火45 s，72℃延伸1 min，循環(huán)35次，最后72℃延伸10 min；H-A的PCR擴(kuò)增程序?yàn)?4℃預(yù)變性3 min，然后94℃變性30 s，47℃退火45 s，72℃延伸1 min，循環(huán)35次，最后72℃延伸10 min。產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測(cè)合格后送至南京金斯瑞生物科技公司進(jìn)行純化與測(cè)序，樣本基于ABI 3730XL測(cè)序平臺(tái)進(jìn)行正向測(cè)序。

表1 5個(gè)水仙歸化群體的采集信息

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

將獲得的單基因序列分別在BioEditor 7.0.9軟件[14]中的ClustalW模塊進(jìn)行聯(lián)配與手工校正。使用SequenceMatrix 1.7.8軟件包[15]進(jìn)行雙基因序列的拼接。利用DNASP 5.10.1軟件[16]進(jìn)行聯(lián)合序列的特征分析，包括統(tǒng)計(jì)單倍型(haplotype)數(shù)量()、平均核苷酸差異數(shù)(average number of nucleotide,)、單倍型多樣性(haplotype diversity,d)、核苷酸多樣性(nucleotide diversity,)、遺傳分化系數(shù)(G)、基因流(N)等遺傳多樣性指數(shù)以及種群歷史動(dòng)態(tài)的中性檢驗(yàn)(neutrality test)和失配分析(mismatch distribution), 群體間基因流采用N=(1/G-1)/2[20]進(jìn)行估算。使用Arlequin 3.1軟件[17]進(jìn)行分子方差分析以判別變異的來(lái)源。使用PermutCpSSR 2.0軟件包(http://www. pierroton.inra.fr/genetics/labo/Software/Permut/)中的U統(tǒng)計(jì)法檢驗(yàn)(1 000次重復(fù)置換)比較遺傳分化系數(shù)G和N,以檢驗(yàn)水仙是否存在譜系地理結(jié)構(gòu)。使用GenAlEx 6.503軟件包[18]在Microsoft Excel 2007中進(jìn)行Mantal檢驗(yàn)，分析水仙群體遺傳距離與地理距離是否存在地理隔離(isolation by distance, IBD)作用。在Network 5.0.1.0軟件(http://www.fluxus- engineering.com/sharenet.htm/)中基于最大簡(jiǎn)約性原則的中介鄰接網(wǎng)絡(luò)法(median-joining network, MJ)獲得單倍型中介網(wǎng)絡(luò)圖，并在Adibe Illustrator CS6軟件中進(jìn)行優(yōu)化布局與重新繪制。并利用MEGA 5.0軟件[19]構(gòu)建以鄰接法(neighbor-joining, NJ)為算法的單倍型系統(tǒng)發(fā)育樹(shù), 并采用bootstrap法(1 000次重復(fù)置換)以檢驗(yàn)聚類(lèi)樹(shù)各分支置信度。

2 結(jié)果和分析

2.1 序列變異與單倍型分布

兩個(gè)cpDNA片段拼接后的序列總長(zhǎng)度為1 443 bp,G+C含量為33.10%，共檢測(cè)出8個(gè)變異位點(diǎn)(7個(gè)簡(jiǎn)約信息位點(diǎn)和1個(gè)自裔位點(diǎn))；共定義了6個(gè)單倍型(haplotypes, Hap1~Hap6)，其中K對(duì)應(yīng)的變異位點(diǎn)為673、767、801、807和809，H-A對(duì)應(yīng)的變異位點(diǎn)為906、1 333和1 383。根據(jù)溯祖理論, 貢獻(xiàn)單倍型(Hap2)位于單倍型中介網(wǎng)絡(luò)中心, 應(yīng)為5個(gè)地理群體的共有古老單倍型，而私自有單倍型則分散在各個(gè)群體，應(yīng)為近期輻射分化出來(lái)的年輕單倍型(圖1, 表2)。

2.2 遺傳多樣性評(píng)估

從表3可見(jiàn)，群體的總體單倍型多樣性(d)為0.363±0.086，核苷酸多樣性()為0.000 49±0.000 19；按遺傳多樣性水平，5個(gè)水仙歸化群體依次為大檑山嶼(DLSY)>平潭島(PTD)>南麂島(NJD)=漳州(ZZ)>崇明島(CMD)。

2.3 群體遺傳結(jié)構(gòu)與分化

從表4可見(jiàn)，兩兩群體間的基因流較強(qiáng), 且所有群體間的基因流都>4；分子方差分析(AMOVA) 結(jié)果表明，群體內(nèi)變異占了總變異的91.98%，為遺傳變異的主要來(lái)源；固定指數(shù)(F)為0.080 22，群體間的遺傳分化較弱(表5)；通過(guò)1 000次重復(fù)計(jì)算獲得的群體遺傳分化系數(shù)(N=0.020<G=0.031;<0.05), 表明水仙在物種水平上不具有明顯的譜系結(jié)構(gòu)。

圖1 中國(guó)東部沿海5個(gè)水仙歸化居群的葉綠體單倍型分布

表2 基于葉綠體matK+trnH-psbA序列的水仙6個(gè)單倍型中的變異位點(diǎn)

*: 該堿基變異位點(diǎn)與Hap1相同。

*: Variation site of the base is the same as Hap1.

表3 5個(gè)水仙歸化群體的遺傳多樣性

表4 5個(gè)水仙歸化群體間的遺傳分化系數(shù)(Gst, 對(duì)角線(xiàn)下)及基因流(Nm, 對(duì)角線(xiàn)上)

水仙歸化群體間的遺傳距離(G)與地理距離(K, 自然對(duì)數(shù)轉(zhuǎn)化)矩陣的Mantel檢驗(yàn)結(jié)果表明,兩者之間具有顯著的線(xiàn)性關(guān)系(=0.929,=0.02< 0.05) (圖2)。

2.4 群體擴(kuò)張事件檢測(cè)

選擇中性檢驗(yàn)和失配分析推測(cè)水仙群體的擴(kuò)張事件。通常用Tajima的指數(shù)[21]可以反映較大時(shí)間尺度上的種群事件，而Fu的*指數(shù)[22]對(duì)近期的事件比較敏感。本研究結(jié)果表明，中性檢驗(yàn)均無(wú)顯著性意義(D=-0.165 395,>0.05;*=0.590 83,> 0.10)。同時(shí)，錯(cuò)配分布分析(MDA)中，觀(guān)測(cè)值與期望值比較得到的雙峰曲線(xiàn)結(jié)果也表明水仙群體背離了快速擴(kuò)張模型的假設(shè)，即群體在進(jìn)化過(guò)程中未經(jīng)歷過(guò)快速擴(kuò)張(圖3)。

表5 5個(gè)水仙歸化群體的分子方差(AMOVA)分析

圖2 水仙歸化群體遺傳距離與地理距離的Mantel檢測(cè)

圖3 水仙49個(gè)個(gè)體錯(cuò)配分布分析

Fig. 3 Mismatch distribution analysis for 49 individuals ofvar.

2.5 單倍型系統(tǒng)發(fā)育關(guān)系

基于遺傳距離鄰接法構(gòu)建的單倍型系統(tǒng)發(fā)育樹(shù)(圖4)可知，與中介網(wǎng)絡(luò)圖結(jié)果一致(圖1)，水仙6個(gè)單倍型明顯被劃分為2個(gè)譜系分支，Hap4與Hap5聚為一大支，其余單倍型聚為另一大支，靴帶支持率(bootstrap)分別為80%和88%。其中Hap1、Hap3與Hap2又單獨(dú)聚為一小支, 從而表現(xiàn)出更近的親緣關(guān)系。

3 討論

開(kāi)展遺傳變異分析有助于評(píng)價(jià)物種的進(jìn)化潛力與環(huán)境適應(yīng)能力，也是開(kāi)展保護(hù)生物學(xué)研究的核心[23-25]?；赾pDNA序列變異的結(jié)果，從整體上看, 我國(guó)歸化水仙的遺傳多樣性水平較低，這與前人應(yīng)用RAPD等其他標(biāo)記手段的研究結(jié)果一致[10]; Mantel檢測(cè)的相關(guān)系數(shù)較高且達(dá)到顯著水平(= 0.929,<0.05)，符合地理隔離IBD模型，歸化群體遺傳距離更受地理距離影響顯著，表明海島間的地理隔離已經(jīng)成為該種基因流擴(kuò)散的主要屏障。細(xì)胞染色體研究表明，中國(guó)水仙主要為同源三倍體[26],在自然狀況下水仙只能采取以母鱗莖為中心的側(cè)球無(wú)性繁殖方式，大多呈聚集分布，缺乏長(zhǎng)距離的傳播機(jī)制[27]，使得居群內(nèi)交流遠(yuǎn)大于居群間。本研究中的分子方差分析結(jié)果亦顯示水仙群體變異主要來(lái)源于種群內(nèi)(91.98%)。此外，水仙棲息地嚴(yán)酷生境尤其是土壤的貧瘠及基鹽脅迫無(wú)疑降低了水仙種群的生態(tài)適應(yīng)性，而其物候花期寒冷的氣溫與強(qiáng)烈的海風(fēng)等因素或可導(dǎo)致長(zhǎng)距離的傳粉者缺失, 奠基者效應(yīng)的出現(xiàn)大大降低了群體遺傳多樣性。

圖4 水仙6個(gè)單倍型的系統(tǒng)發(fā)育樹(shù)。分支上數(shù)值為1 000次重復(fù)得到的靴帶支持率。

植物在長(zhǎng)期的遷移演化過(guò)程中為適應(yīng)不同生境形成了不同的地理生態(tài)類(lèi)群，生態(tài)隔離(isolation by ecology, IBE)作用往往可以加速物種的遺傳分化[28]。前人的研究均表明歸化水仙在我國(guó)存在一定的遺傳變異特征[29]，陳曉慧等[30]基于居群統(tǒng)計(jì)，分析比較了我國(guó)沿海4個(gè)島嶼8個(gè)群體的形態(tài)學(xué)特征差異, 認(rèn)為南麂產(chǎn)的水仙具有植株小、花朵大的特點(diǎn)，是一種珍稀遺傳資源。而結(jié)合本文的cpDNA單倍型分布來(lái)看(圖1)，各群體均含有1個(gè)以上的私有單倍型，暗示我國(guó)沿海島嶼的環(huán)境異質(zhì)性較高，水仙居群可能存在局域適應(yīng)(local adaptation)現(xiàn)象；水仙屬植物通常性喜溫暖、濕潤(rùn)及沙質(zhì)土壤，在我國(guó)東部歸化群體中，大檑山嶼(DLSY)含有2個(gè)私有單倍型導(dǎo)致其遺傳多樣性綜合排序第一，這可能是由于其地處南麂列島國(guó)家級(jí)自然保護(hù)區(qū)核心區(qū)的地理優(yōu)勢(shì)，在區(qū)系劃分上屬于華東植物區(qū)系南北過(guò)渡特征，亞熱帶海洋性季風(fēng)區(qū)及臺(tái)灣暖流帶來(lái)的水熱氣候條件相對(duì)充沛[31]；基于群落生態(tài)學(xué)研究亦表明,大檑山嶼具有適宜水仙生長(zhǎng)的排水較好、土壤鹽分適中與溫度適宜的微生境[27]，可能使得大檑山嶼群體產(chǎn)生了適應(yīng)性分化。有意思的是，該群體與緊鄰的南麂島(NJD)共享一個(gè)單倍型(Hap3)，且基因流(N=12.000)最高、遺傳分化(G=0.040)最小，相較其它群體而言顯示出最近的親緣關(guān)系；而雖然在物種水平上遺傳多樣性較低，但水仙各群體間的平均基因流均>4，說(shuō)明群體間的基因流持續(xù)發(fā)生(表5)，結(jié)合單倍型系統(tǒng)發(fā)育分析，聚類(lèi)結(jié)果也與地理距離不一致(圖4)。根據(jù)實(shí)際調(diào)查走訪(fǎng)得知，造成上述現(xiàn)象很可能與歷史上海島漁民人為多次移栽事件密切相關(guān)，即群體間可能經(jīng)常性交換遷徙個(gè)體; cpDNA單倍型中介網(wǎng)絡(luò)圖還顯示，平潭島(PTD)與其余群體的親緣關(guān)系相對(duì)最遠(yuǎn)，該群體是歸化水仙中為數(shù)不多的重瓣類(lèi)型‘玉玲瓏’，其副花冠和雄蕊花瓣?duì)?，黃白相間，與大部分單瓣、純白花色類(lèi)型截然不同[32],有關(guān)其花色及花形變異的生理與分子機(jī)制還有待進(jìn)一步研究，但無(wú)疑反映了人類(lèi)經(jīng)濟(jì)活動(dòng),對(duì)水仙遺傳分布格局帶來(lái)的持續(xù)影響。

目前在我國(guó)水仙市場(chǎng)上，栽培品種中除了‘玉玲瓏’，還有‘金三角’、‘金邊’和‘云香’等，但都分布不均，栽培面積極少且產(chǎn)地不明，品種類(lèi)群十分單一；而且由于水仙很難通過(guò)有性繁殖同其他物種進(jìn)行雜交獲得新的基因，長(zhǎng)期的無(wú)性繁殖，勢(shì)必導(dǎo)致種群衰退、品質(zhì)降低等問(wèn)題，這將嚴(yán)重制約水仙花卉產(chǎn)業(yè)的可持續(xù)發(fā)展[33]。因此，開(kāi)展我國(guó)歸化水仙資源種群動(dòng)態(tài)和遺傳變異的研究進(jìn)而獲得基礎(chǔ)的遺傳學(xué)背景信息是非常必要的。本研究綜合利用了植物分子譜系地理學(xué)的方法，首次對(duì)我國(guó)東部沿海水仙5個(gè)代表歸化群體49個(gè)樣本進(jìn)行遺傳多樣性檢測(cè)及分析，初步闡明了該物種在區(qū)域尺度上的遺傳多樣性水平與空間分布，這無(wú)疑有助于下一步水仙衰退種群的恢復(fù)保育工作的開(kāi)展。同樣以南麂島產(chǎn)水仙尤其是大檑山嶼為代表，基于群體遺傳cpDNA標(biāo)記的分子證據(jù)，該歸化群體應(yīng)為中國(guó)東部沿海一個(gè)特殊且重要的基因型或生態(tài)型分化類(lèi)群，目前雖然已經(jīng)列入核心保護(hù)區(qū)，但由于靠近旅游區(qū)仍然面臨人類(lèi)經(jīng)濟(jì)活動(dòng)的日益影響，應(yīng)結(jié)合沿海島嶼特殊的生態(tài)系統(tǒng)進(jìn)行綜合管理，做好持續(xù)的宣傳、保護(hù)與管理工作。短期而言，就地保護(hù)可能是最為有效的措施，可以在保護(hù)區(qū)內(nèi)增加護(hù)欄、控制伴生入侵雜草及定期巡護(hù)工作；長(zhǎng)期而言，一方面通過(guò)原地建立固定樣地，開(kāi)展對(duì)種群及物候的動(dòng)態(tài)監(jiān)測(cè)，另一方面通過(guò)離體培養(yǎng)、人工種苗擴(kuò)繁、自然回歸試驗(yàn)等措施，從而促進(jìn)其種群的正向更新。同時(shí)也應(yīng)發(fā)揮其重要的科學(xué)價(jià)值，將該區(qū)域優(yōu)先選為核心種質(zhì)苗圃或遺傳基因庫(kù)建設(shè)的遺傳材料來(lái)源地，或可從中選擇適宜的親本以提高雜交育種效率?？傊?本研究可為海島特色種質(zhì)資源的綜合評(píng)價(jià)、科學(xué)保護(hù)與開(kāi)發(fā)利用提供理論依據(jù)。

[1] CHEN X Q, WU Y X. Additional notes on Chinese sacred lily, with a discussion with Mr. Lu [J]. J Wuhan Bot Res, 1991, 9(1): 70-74. 陳心啟, 吳應(yīng)祥. 中國(guó)水仙續(xù)考——與盧履范同志商榷 [J]. 武漢植物學(xué)研究, 1991, 9(1): 70-74.

[2] LI X F, SHAO X H, DENG X J, et al. Necessity of high temperature for the dormancy release ofvar.[J]. J Plant Physiol, 2012, 169(14): 1340-1347. doi: 10.1016/j.jplph.2012.05.017.

[3] BIAN A A, PAN D M. Effects of salt stress on growth and inorganic ion distribution inL. var.Roem. seedlings [J]. Hort Sci, 2018, 53(8): 1152-1156. doi: 10.21273/hortsci13145-18.

[4] LI Y, CAI P, JIANG L L, et al. Propagation of wild daffodils in Suzhou City, Jiangsu Province [J]. Jiangsu J Agric Sci, 2010(6): 276-278. 李穎, 蔡平, 姜麗麗, 等. 蘇州野生水仙繁殖技術(shù)研究 [J]. 江蘇農(nóng)業(yè)科學(xué), 2010(6): 276-278.

[5] ZHANG X Q, GAO J, PENG Z H. Callus induction, differentiation and histological observation ofvar.floral organs [J]. For Res, 2013, 26(3): 320-325. 張曉晴, 高健, 彭鎮(zhèn)華. 中國(guó)水仙花器官愈傷組織誘導(dǎo)、分化及組織學(xué)觀(guān)察 [J]. 林業(yè)科學(xué)研究, 2013, 26(3): 320-325.

[6] WU J C, WANG G Q, ANWAR M, et al. Cloning and functional analysis of R2R3-MYB geneinvar.[J]. Acta Hort Sin, 2018, 45(7): 1327-1337. doi: 10.16420/j.issn.0513- 353x.2017-0889.吳嘉誠(chéng), 王桂青, ANWAR M, 等. 中國(guó)水仙R2R3-MYB基因的克隆和功能研究 [J]. 園藝學(xué)報(bào), 2018, 45(7): 1327-1337. doi: 10.16420/j.issn.0513-353x.2017-0889.

[7] QU R Z, HOU L, Lü H L, et al. The gene flow of population genetic structure [J]. Hereditas, 2004, 26(3): 377-382. 曲若竹, 侯林, 呂紅麗, 等. 群體遺傳結(jié)構(gòu)中的基因流 [J]. 遺傳, 2004, 26(3): 377-382.

[8] WANG Z F, LIAN J Y, HUANG G M, et al. Genetic groups in the common plant speciesand their associations with topographic habitats [J]. Oikos, 2012, 121(12): 2044-2051. doi: 10.1111/j.1600-0706.2012.20483.x.

[9] LIU P P, LIAO H B, JIANG Q C, et al. Spatial genetic structure ofin Wuguishan, Zhongshan, Guangdong Province [J]. J Trop Subtrop Bot, 2019, 27(1): 65-73. doi: 10.11926/jtsb.3917.劉盼盼, 廖浩斌, 蔣謙才, 等. 廣東省中山市五桂山土沉香空間遺傳結(jié)構(gòu) [J]. 熱帶亞熱帶植物學(xué)報(bào), 2019, 27(1): 65-73. doi: 10.11926/ jtsb.3917.

[10] LIU X M, ZHANG X Z, SHI Y M, et al. Genetic diversity analysis of ninebased on morphological characteristics and random amplified polymorphic DNA markers [J]. Hort Sci, 2017, 52(2): 212- 220. doi: 10.21273/HORTSCI11171-16.

[11] CHEN L J, TIAN H Q, WU J. The study on RAPD fingerprints ofin China and Europe [J]. J Trop Subtrop Bot, 2003, 11(2): 177-180. doi: 10.3969/j.issn.1005-3395.2003.2.019.陳林姣, 田惠橋, 武劍. 中國(guó)水仙與歐洲水仙品種RAPD指紋的研究 [J]. 熱帶亞熱帶植物學(xué)報(bào), 2003, 11(2): 177-180. doi: 10.3969/j. issn.1005-3395.2003.2.019.

[12] XU G B. Intraspecific phylogeography and its application in conser- vation strategies of plant genetic diversity [J]. J CS Univ For Technol, 2011, 31(12): 1-6. 徐剛標(biāo). 種內(nèi)譜系地理學(xué)及在植物遺傳多樣性保護(hù)中的應(yīng)用 [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2011, 31(12): 1-6.

[13] ZHANG C F, YU W T, LI M, et al. Identification of DNA barcodes in[J]. J Fujian Agric For Univ (Nat Sci), 2018, 47(4): 452-460. 張程飛, 于文濤, 李敏, 等. 水仙屬DNA條形碼鑒定技術(shù) [J].福建農(nóng)林大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 47(4): 452-460.

[14] HALL T A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT [J]. Nucl Acids Symp Ser, 1999, 41: 95-98.

[15] VAIDYA G, LOHMAN D J, MEIER R. SequenceMatrix: Concate- nation software for the fast assembly of multi-gene datasets with character set and codon information [J]. Cladistics, 2011, 27(2): 171- 180. doi: 10. 1111/j.1096-0031.2010.00329.x.

[16] LIBRADO P, ROZAS J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data [J]. Bioinformatics, 2009, 25(11): 1451-1452. doi: 10.1093/bioinformatics/btp187.

[17] EXCOFFIER L, LAVAL G, SCHNEIDER S. Arlequin (Version 3.0): An integrated software package for population genetics data analysis [J]. Evol Bioinform, 2005, 1: 47-50. doi: 10.1177/117693430500100003.

[18] PEAKALL R, SMOUSE P E. GenAlEx 6. 5: Genetic analysis in Excel. Population genetic software for teaching and research: An update [J]. Bioinformatics, 2012, 28(19): 2537-2539. doi: 10.1093/bioinformatics/ bts460.

[19]TAMURA K, PETERSON D, PETERSON N, et al. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods [J]. Mol Biol Evol, 2011, 28(10): 2731-2739. doi: 10.1093/molbev/msr121.

[20] BAO Y X, CHENG H Y, ZHOU X W, et al. Genetic diversity of three populations of the black muntjac () [J]. Acta Ecol Sin, 2008, 28(8): 4030-4036. 鮑毅新, 程宏毅, 周襄武, 等. 黑麂() 3個(gè)種群的遺傳多樣性 [J]. 生態(tài)學(xué)報(bào), 2008, 28(8): 4030-4036.

[21] TAJIMA F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism [J]. Genetics, 1989, 123(3): 585-595.

[22] FU Y X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection [J]. Genetics, 1997, 147 (2): 915-925.

[23] MENEZES T, ROMEIRAS M M, DE SEQUEIRA M M, et al. Phylo- genetic relationships and phylogeography of relevant lineages within the complex Campanulaceae family in Macaronesia [J]. Ecol Evol, 2018, 8(1): 88-108. doi: 10.1002/ece3.3640.

[24] JIAN H Y, LI S F, GUO J L, et al. High genetic diversity and differen- tiation of an extremely narrowly distributed and critically endangered decaploid rose (): Implications for its conservation [J]. Conserv Genet, 2018, 19(4): 761-776. doi: 10.1007/s10592-018-1052-0.

[25] ZHANG L W, HAN G X. A review on the relationships between plant genetic diversity and ecosystem functioning [J]. Chin J Plant Ecol, 2018, 42(10): 977-989. doi: 10.17521/cjpe.2018.0013.張儷文, 韓廣軒. 植物遺傳多樣性與生態(tài)系統(tǒng)功能關(guān)系研究進(jìn)展 [J]. 植物生態(tài)學(xué)報(bào), 2018, 42(10): 977-989. doi: 10.17521/cjpe.2018.0013.

[26] HONG L, LIANG G L. Research advance in genetic and biological characteristics of[J]. Subtrop Plant Sci, 2011, 40(2): 89-96. 洪林, 梁國(guó)魯. 水仙屬植物遺傳學(xué)及生物學(xué)特性研究進(jìn)展 [J]. 亞熱帶植物科學(xué), 2011, 40(2): 89-96.

[27] ZHU H, CAI H C, YOU L X, et al. Analyses on species diversity, environmental interpretation and spatial distribution pattern of natural population ofvar.on Daleishan Island of Nanji Islands in Zhejiang Province [J]. J Plant Resour Environ, 2017, 26(3): 100-108. doi: 10.3969/j.issn.1674-7895.2017.03.13. 朱弘, 蔡厚才, 尤祿祥, 等. 浙江南麂列島大檑山嶼水仙自然居群的物種多樣性、環(huán)境解釋及空間分布格局分析 [J]. 植物資源與環(huán)境學(xué)報(bào), 2017, 26(3): 100-108. doi: 10.3969/j.issn.1674-7895.2017.03.13.

[28] LI M, WANG S X, GAO B J. Analysis of genetic diversity of Chinese pine () natural secondary forest populations and correlation with theirs habitat ecological factors [J]. Acta Ecol Sin, 2013, 33(12): 3602-3610. 李明, 王樹(shù)香, 高寶嘉. 油松天然次生林居群遺傳多樣性及與產(chǎn)地地理氣候因子的關(guān)聯(lián)分析 [J]. 生態(tài)學(xué)報(bào), 2013, 33(12): 3602-3610.

[29] CHEN D F, GAO J, PENG Z H. Advances and perspectives in studies onplants [J]. Sci Silv Sin, 2008, 44(3): 140-146. 陳段芬, 高健, 彭鎮(zhèn)華. 水仙屬植物研究進(jìn)展 [J]. 林業(yè)科學(xué), 2008, 44(3): 140-146.

[30] CHEN X H, ZHENG K, CHEN Z H, et al. Morphological variation in wild populations ofvar.and taxonomy position ofin Nanji Islands [J]. J Anhui Agric Sci, 2012, 40(33): 16041-16044. 陳曉慧, 鄭康, 陳珍慧, 等. 中國(guó)水仙野生居群形態(tài)變異及其南麂水仙分類(lèi)地位的研究 [J]. 安徽農(nóng)業(yè)科學(xué), 2012, 40(33): 16041-16044.

[31] ZHU H, KU W P, RONG J T, et al. Species diversity and floristic characteristics of vascular plants in Nanji Island, Zhejiang Provience [J]. Plant Div Resour, 2015, 37(6): 713-720. doi: 10.7677/ynzwyj2015 15050.朱弘, 庫(kù)偉鵬, 戎建濤, 等. 浙江南麂島陸生維管束植物多樣性及區(qū)系特征 [J]. 植物分類(lèi)與資源學(xué)報(bào), 2015, 37(6): 713-720. doi: 10. 7677/ynzwyj201515050.

[32] LIN X H. Morphological and anatomical study on threevar.variant strains (strains) [J]. J CS Univ For Technol, 2012, 32(10): 135-140. 林曉紅. 三個(gè)中國(guó)水仙變異品系(株系)的形態(tài)解剖研究 [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2012, 32(10): 135-140.

[33] SHI S J, YE R M. Plant classification ofand current status of utilization and study on Chinesegermplasm resource [J]. Fujian Agric Sci Technol, 2016, 47(4): 52-58. 施世強(qiáng), 葉如夢(mèng). 水仙屬植物分類(lèi)及中國(guó)水仙資源利用研究現(xiàn)狀 [J]. 福建農(nóng)業(yè)科技, 2016, 47(4): 52-58.

Genetic Diversity ofvar.Populations in Eastern Coast of China

ZHU Hong1, CAI Hou-cai2, LI Yong-fu1, CHEN Wan-dong2, CHEN Lin1, YI Xian-gui1, LI Meng1, DUAN Yi-fan1*, WANG Xian-rong1*

(1. Co-Innovation Center for Sustainable Forestry in Southern China, College of Biology and the Environment, Key Laboratory of State Forestry Administration on Subtropical Forest Biodiversity Conservation,Nanjing Forestry University,Nanjing 210037, China; 2. Nanji Islands National Marine Natural Reserve Administration Bureau,Wenzhou 325401, Zhejiang, China)

To reveal the genetic diversity of naturalizedvar.populations in eastern coast of China, 49 individuals from five representative populations distributed in Shanghai, Zhejiang and Fujian Proviences were evaluated by using two chloroplast DNAK andHA fragments. The results showed that total of six haplotypes were identi?ed from the combined two sequences with an alignment length of 1 443 bp. The genetic diversity of naturalized population was in order of DLSY>PTD>NJD=ZZ>CMD. AMOVA analysis showed that the variation within populations was the main source of the genetic variation (91.98%) and low genetic differentiation between populations (F=0.080 22). Insigni?cant phylogeographic structure was detected at species level (Nt=0.020<G=0.031;<0.05). Mantel test showed a significant linear correlation between genetic distance and geographical distance (=0.929,=0.02<0.05). The neutrality test combined bimodalcurve in mismatch distribution analysis both rejected the hypothesis of the demographic expansion of these populations. Analysis of median network and phylogenetic NJ tree based on the distribution of haplotype both divided all the populations into two major branches. These indicated a low level of genetic diversity with a weak genetic differentiation amongvar.populations in eastern coast of China. Genetic variation mainly occurs in populations, while species had experienced no demographic expansion. These might likely owing to the gene flows which were caused by the combination of islands segregation, biological characteristics, habitat heterogeneity as well as human disturbance.

var.;Eastern coast ofChina; Naturalized population; cpDNA; Genetic structure

10.11926/jtsb.4045

2019-01-21

2019-03-15

南麂列島國(guó)家級(jí)海洋自然保護(hù)區(qū)研究項(xiàng)目(NJKJ-2016-006)資助

This work was supported by the research program of Nanji Islands national marine natural reserve (Grant No. NJKJ-2016-006).

朱弘(1990~)，男，博士研究生，主要從事亞熱帶森林生態(tài)與譜系地理學(xué)研究。E-mail: 1107401987@qq.com

Corresponding author. E-mail: wangxianrong66@nifu.edu.cn; 67181278@qq.com