基于葉綠體DNA+PCRRFLP分析柿及其近緣種親緣關(guān)系

1、園藝學(xué)報(bào)，（）： 2014411224742480 http: / www. ahs. ac. cn acta horticulturae sinica e-mail: yuanyixuebao 收稿日期：20140724；修回日期：20141031 基金項(xiàng)目：國(guó)家十二五農(nóng)村領(lǐng)域科技計(jì)劃課題（2013bad14d0502） 基于葉綠體dna pcr-rflp分析柿及其近緣種親緣關(guān)系 梁晉軍，梁玉琴，傅建敏* （國(guó)家林業(yè)局泡桐研究開(kāi)發(fā)中心，中國(guó)林業(yè)科學(xué)研究院經(jīng)濟(jì)林研究開(kāi)發(fā)中心，鄭州 450003） 摘 要：應(yīng)用5對(duì)葉綠體dna（cpdna）特異引物對(duì)柿及其近緣的8個(gè)種共32個(gè)基因型的cpdna進(jìn)

2、行pcr擴(kuò)增。采用7種限制性?xún)?nèi)切酶對(duì)其擴(kuò)增產(chǎn)物酶切位點(diǎn)變異分析表明，在供試柿屬植物中存在豐富的種間多態(tài)性，但供試22個(gè)柿種內(nèi)基因型未檢測(cè)到變異位點(diǎn)。應(yīng)用ntsys中upgma法聚類(lèi)和主坐標(biāo)分析表明：柿與金棗柿最近緣，其次分別是油柿、云南野毛柿、浙江柿和君遷子，而與美洲柿和烏柿相距較遠(yuǎn)。野柿（浙江）為柿的變種，聚類(lèi)結(jié)果也表明野柿（浙江）與柿有一定差異。 關(guān)鍵詞：柿；cpdna；pcr-rflp；遺傳多樣性；親緣關(guān)系 中圖分類(lèi)號(hào)：s 665.2 文獻(xiàn)標(biāo)志碼：a 文章編號(hào)：0513-353x（2014）12-2474-07 genetic relationships of diospyros kak

3、i and related diospyros species using chloroplast dna pcr-rflp markers liang jin-jun，liang yu-qin，and fu jian-min* （china paulownia research center，non-timber forest research and development center，chinese academy of forestry，zhengzhou 450003，china） abstract：five special primer pairs of chloroplast

4、genome were used to amplified chloroplast dna （cpdna）regions in 8 diospyros species including 32 genotypes. the products were digested by seven restriction enzymes. the results showed high interspecific cpdna variations within the genus diospyros. twenty-two d. kaki genotypes were detected no visual

5、ly variation in this study. d. kaki had closely relationship withjinzaoshi，followed by d. oleifera，d. artrotricha，d. lotus and d. glaucifolia，but distant from d. virginiana and d. cathayensis in diagram based on the upgma cluster analysis by ntsyspc 2.10e program. wild persimmon（zhejiang）is a variet

6、y varietas of persimmon. the result showed that wild persimmon（zhejiang）is not closely similar to d. kaki. key words：diospyros kaki；cpdna；pcr-rflp；genetic diversity；genetic relationship 柿（diospyros kaki thunb.）品種多數(shù)為六倍體，2n = 6x = 90，但也有少量九倍體，2n = 9x = 135（zhuang et al.，1990）。前人已經(jīng)通過(guò)各種分子標(biāo)記對(duì)柿及其近緣種親緣關(guān)系進(jìn)行

7、研究，但結(jié)果差異較大。yuri和shozo（1994）通過(guò)對(duì)柿及其近緣種葉綠體dna（cpdna）和線(xiàn)粒體dna（mtdna） * 通信作者 author for correspondence（e-mail：fjm371） 12期 梁晉軍等：基于葉綠體dna pcr-rflp分析柿及其近緣種親緣關(guān)系 2475 的rflp分析，發(fā)現(xiàn)老鴉柿與柿最近緣；keizo等（1998）對(duì)柿屬植物cpdna pcr-rflp分析的結(jié)果表明柿與美洲柿最近緣；hu等（2008）、guo和luo（2011）分別用cpdna pcr-rflp和ssr分析表明柿與君遷子最近緣；keizo等（2008）用nrdna的it

8、s區(qū)和cpdna matk序列分析的結(jié)果表明柿與油柿最近緣。此外，hu等（2008）和tang等（2014）分別用cpdna pcr-rflp和its，matk序列變異分析證明金棗柿為一個(gè)新種。 由于柿倍性復(fù)雜，其核基因組結(jié)構(gòu)和序列也是復(fù)雜多變，從核基因組開(kāi)發(fā)分子標(biāo)記有一定難度，相較于核基因，葉綠體基因組結(jié)構(gòu)和序列非常保守，所以葉綠體基因組能很好地用于種間甚至種內(nèi)系統(tǒng)進(jìn)化和親緣關(guān)系的研究（randall et al.，1998）。目前，基于cpdna的pcr-rflp的分子標(biāo)記已廣泛用于遺傳多樣性、系統(tǒng)進(jìn)化和植物種類(lèi)鑒定等研究（claudia et al.，2004；pharmawati et

9、 al.，2004；bao et al.，2005）。本研究中通過(guò)對(duì)4個(gè)柿近緣種（君遷子、浙江柿、金棗柿和油柿）cpdna全基因組測(cè)序結(jié)果對(duì)比分析，找出變異較大的4個(gè)區(qū)段（trnh-matk、rpob-petn、trns-trnfm和ndha）和1個(gè)功能基因（rbcl）片段，并設(shè)計(jì)特異引物對(duì)32份柿屬植物的5個(gè)cpdna區(qū)域進(jìn)行pcr-rflp分析，以期為柿及其近緣種親緣關(guān)系提供有效證據(jù)。 1 材料與方法 1.1 材料 試材為柿及其近緣種共8種，1個(gè)變種32個(gè)基因型（表1）。烏柿、美洲柿和油柿于2014年采自陜西楊陵國(guó)家柿種質(zhì)資源圃，其余于2013年采自國(guó)家林業(yè)局泡桐研究開(kāi)發(fā)中心柿種質(zhì)資源圃。

10、 1.2 cpdna區(qū)域選擇與引物設(shè)計(jì) 2013年6月采用kerstin等（2008）改良的蔗糖密度梯度離心分離法提取4個(gè)柿近緣種（金棗柿、君遷子、浙江柿和油柿）的葉綠體總dna；超聲波隨機(jī)打斷dna，電泳回收所需大小片段，加上接頭進(jìn)行文庫(kù)制備；最后用illumina-miseq測(cè)序儀完成規(guī)定數(shù)據(jù)的測(cè)序。采用calera assembler組裝軟件對(duì)測(cè)序結(jié)果進(jìn)行組裝。使用dogma軟件對(duì)基因進(jìn)行注釋?zhuān)詈蠼?jīng)過(guò)人工校正，用clustalx對(duì)4個(gè)全基因組序列進(jìn)行對(duì)比，找出差異較大的4個(gè)區(qū)段（trnh-matk、rpob-petn、trns-trnfm和ndha）和差異較小的1個(gè)功能基因（rbcl）

11、片段。 針對(duì)所選5個(gè)區(qū)段上、下游，利用primer 5.0軟件設(shè)計(jì)特異引物（表2），并用oligo7.0軟件對(duì)引物進(jìn)行評(píng)價(jià)，選出最優(yōu)特異引物，由生工生物工程（上海）股份有限公司合成。 1.3 dna提取、pcr擴(kuò)增和限制性酶切 1.3.1 dna提取和pcr擴(kuò)增 2013年6月采用ctab法（doyle & doyle，1987）從柿幼嫩葉片中提取植物總dna，2014年4月進(jìn)行補(bǔ)充試驗(yàn)。 2014年4月試材dna在bio-red公司pcr儀上進(jìn)行pcr擴(kuò)增。反應(yīng)總體積為25 µl，內(nèi)含1× pcr緩沖液（含1.5 mmol · l-1 mg2+），250

12、 µmol · l-1 dntp，10 mmol · l-1 tris-hcl，50 mmol · l-1 kcl，0.2 µmol · l-1 primers，1.25 u taq dna polymerase，dna 40 ng。pcr程序?yàn)?4 預(yù)熱3 min，1個(gè)循環(huán)；94 變性45 s，52.7 56 退火1 min（表2），72 延伸1 min 15 s 2 min 30 s，35個(gè)循環(huán)；最后72 延伸5 min，4 保溫。pcr產(chǎn)物通過(guò)gelred染色的2.0%瓊脂糖凝膠電泳檢測(cè)。bio-red公司凝膠成像系統(tǒng)觀察并記錄

13、譜帶。 2476 園 藝 學(xué) 報(bào) 41卷 表1 試驗(yàn)材料 table 1 the materials in this study 學(xué)名 scientific name 倍性水平ploidy level 編號(hào) code 種/品種名 specie/cultivar name 脫澀類(lèi)型astringent type 來(lái)源 source diospyros lotus 2n = 2x = 30 1 君遷子（雌）date plum（female） 陜西 shaanxi 2n = 2x = 30 2 君遷子（雄）date plum（male） 陜西 shaanxi d. glaucifolia 2n =

14、2x = 30 3 浙江柿chekiang persimmon 浙江 zhejiang d. kaki thunb. var. silvestris 2n = 6x = 90 4 野柿（浙江）wild persimmon（zhejiang） 浙江 zhejiang d. oleifera 2n = 2x = 30 5 油柿oily persimmon 浙江 zhejiang d. virginiana l. 2n = 6x = 90 6 美洲柿（雌）common persimmon（female） 美國(guó) usa 2n = 6x = 90 7 美洲柿（雄）common persimmon（mal

15、e） 美國(guó) usa d. sp. 2n = 2x = 30 8 金棗柿jinzaoshi 浙江 zhejiang d. cathayensis 2n = 6x = 90 9 烏柿woolly-flowered persimmon 陜西 shaanxi d. artrotricha 2n = 2x = 30 10 云南野毛柿yunnan yemaoshi 云南 yunnan d. kaki thunb. 2n = 6x = 90 11 富有fuyuu pcna 日本 japan 12 次郎jirou pcna 日本 japan 13 晚御所oku-gosho pcna 日本 japan 14 陽(yáng)

16、豐youhou pcna 日本 japan 15 禪寺丸zenjimaru pvna 日本 japan 16 湖南鳥(niǎo)柿hunan niaoshi pca 湖南 hunan 17 黑柿blank persimmon pca 日本 japan 18 斤柿jin persimmon pca 日本 japan 19 正月 akagaki pvna 日本 japan 20 甜寶蓋baogai tianshi pcna 湖北 hubei 21 秋焰qiuyan pcna 湖北 hubei 22 羅田甜柿luotian tianshi pcna 湖北 hubei 23 磨盤(pán)柿 mopanshi pca 河北

17、 hebei 24 長(zhǎng)安懷抱月changan huaibaoyue pca 陜西 shaanxi 25 臺(tái)灣正柿taiwan zhengshi pca 臺(tái)灣 taiwan 26 牛頭柿niutoushi pca 福建 fujian 27 野柿1號(hào) wild persimmon 1 湖南 hunan 28 野柿2號(hào) wild persimmon 2 湖北 hubei 29 野柿3號(hào) wild persimmon 3 湖北 hubei 30 野柿4號(hào) wild persimmon 4 江蘇 jiangsu 31 野柿5號(hào) wild persimmon 5 湖北 hubei 2n = 9x = 13

18、5 32 平核無(wú)hiratanenashi pva 日本 japan 注：pcna：完全甜柿；pca：完全澀柿；pvna：不完全甜柿；pva：不完全澀柿；：未知。 note：pcna：pollination-constant and non-astringent；pca：pollination constant astringent；pvna：pollination variant non-astringent；pva：pollination variant astringent；：unknown. 表2 cpdna特異引物序列及擴(kuò)增條件 table 2 cpdna the sequence

19、and amplification conditions of cpdna special primer pairs in this study 區(qū)域 region 正向引物 squence forward（53） 反向引物 sequence reverse（53） 退火溫度/ annealing temperature 片段長(zhǎng)度/bp fragment lengthtrnh-matk attcttcgtcgccgtagtaaa ggcggatttggtatttgga 55 2 200 rpob-petn aaagttcttccgtcaagc tccccatactacgagtgaaa 52.7

20、 2 000 trns-trnfm gctatcaaccactcagccatct gaggtcacgggttcaaatc 56 1 000 ndha attcagtttgataacctgctac tttggctatcgtcttgttcta 53 1 500 rbcl ggagggatttatgtcaccac ttgtattcggctcaatcctt 55 1 500 12期 梁晉軍等：基于葉綠體dna pcr-rflp分析柿及其近緣種親緣關(guān)系 2477 1.3.2 酶切反應(yīng) 應(yīng)用apo、hinf、taq、mse、mbo、ava和dra 7種限制性?xún)?nèi)切酶（購(gòu)自neb公司）對(duì)擴(kuò)增出的cpdna片段

21、進(jìn)行酶切。酶切反應(yīng)體積為15 µl，內(nèi)含6 µl pcr產(chǎn)物，2 4 u限制性?xún)?nèi)切酶和1× 限制性?xún)?nèi)切酶反應(yīng)緩沖液，酶切時(shí)間為10 h，除taq反應(yīng)溫度為65 外，其余均為37 。酶切產(chǎn)物通過(guò)gelred染色的2.5%瓊脂糖凝膠電泳檢測(cè)。bio-red凝膠成像系統(tǒng)觀察并記錄譜帶。 根據(jù)擴(kuò)增片段酶切產(chǎn)物的瓊脂糖電泳圖譜上條帶的有或無(wú)標(biāo)記為1或0，建成原始矩陣。應(yīng)用ntsys 2.10，依據(jù)upgma（unweighted pair group method arithmetic averages）法進(jìn)行聚類(lèi)分析。 2 結(jié)果與分析 2.1 柿及其近緣種cpdna pc

22、r-rflp的多態(tài)性 設(shè)計(jì)的5對(duì)特異引物均能在供試32份柿及其近緣種中擴(kuò)增出一條譜帶，通過(guò)2%瓊脂糖凝膠電泳，檢測(cè)出擴(kuò)增產(chǎn)物無(wú)明顯長(zhǎng)度差異。 5個(gè)pcr擴(kuò)增后的產(chǎn)物經(jīng)過(guò)apo、hinf、taq、mse、mbo、ava和dra 7種限制性?xún)?nèi)切酶酶切，35對(duì)特異引物/限制性?xún)?nèi)切酶組合中有22對(duì)（62.9%）檢測(cè)到明顯的種間差異，共獲得181個(gè)片段，其中105個(gè)片段（58%）呈現(xiàn)多態(tài)性。在rpob-petn、trns-trnfm、ndha、rbcl區(qū)域均未檢測(cè)到ava的酶切位點(diǎn)，trnh-matk、rpob-petn、ndha區(qū)域與其余6種限制性?xún)?nèi)切酶組合均存在豐富的種間多態(tài)性，說(shuō)明本試驗(yàn)所選區(qū)域種

23、間存在較大變異；而rbcl區(qū)域僅有一個(gè)酶切組合存在多態(tài)性，說(shuō)明在關(guān)鍵功能基因片段（rbcl）處柿近緣種存在較少變異，以此反推所選近緣種均與柿親緣關(guān)系較近。在所有35對(duì)酶切組合中未檢測(cè)到22個(gè)柿種內(nèi)試材有酶切位點(diǎn)差異，表明柿種內(nèi)葉綠體基因組十分保守。 由圖1可以看出，依據(jù)trnh-matk/dra、rpob-petn/mbo兩種區(qū)域/酶切組合可明確將8個(gè)種，1個(gè)變種，分為9組：1和2：君遷子；3：浙江柿；4：野柿（浙江）；5：油柿；6和7：美洲柿；8：金棗柿；16：烏柿；10：云南野毛柿；11 15：柿品種。在22對(duì)多態(tài)性酶切組合中，僅trns-trnfm/apo組合能檢測(cè)出君遷子（1和2）的雌

24、雄位點(diǎn)變異（圖1，d）；僅trnh-matk/dra（圖1，a）組合能檢測(cè)出柿變種野柿（浙江）與柿位點(diǎn)變異（圖1，a中4與11 15）。ndha/mbo（圖1，c）組合可能有存在柿種內(nèi)些微長(zhǎng)度差異（未記錄新條帶），尚待進(jìn)一步驗(yàn)證。 材料16 32全部為無(wú)差異位點(diǎn)柿種內(nèi)品種，未單獨(dú)列出。 2.2 柿及其近緣種間的親緣關(guān)系 根據(jù)所得到酶切位點(diǎn)變異數(shù)據(jù)建成的原始矩陣，利用ntsys-pc version 2.10軟件計(jì)算sm相似指數(shù)，upgma法進(jìn)行聚類(lèi)（圖2），柿種內(nèi)22個(gè)品種無(wú)差異歸為柿種下；同時(shí)進(jìn)行8個(gè)種種間主坐標(biāo)分析（圖3），雌雄（君遷子和美洲柿）歸為一種，柿種內(nèi)22個(gè)品種和1個(gè)變種野柿（浙

25、江）歸為柿種下。在圖2和圖3中，六倍體柿與同為多倍體的美洲柿和烏柿相似性最差，分布位置最遠(yuǎn)，表明柿與烏柿和美洲柿親緣關(guān)系相對(duì)較遠(yuǎn)，盡管柿與烏柿、美洲柿倍性相同，但cpdna序列變異較大，表明兩者有不同的起源祖先，這與前人研究結(jié)果（nakasuka et al.，2002；choi et al.，2003；guo & luo，2006；guo et al.，2006；hu et al.，2008）一致；柿與金棗柿、油柿和云南野毛柿聚為一類(lèi)，且柿與二倍體金棗柿相似度最高，親緣關(guān)系最近，該結(jié)果與hu等（2008）和guo等（2011）關(guān)于柿與君遷子最近緣的結(jié)果不同。 2478 園 藝 學(xué) 報(bào)

26、 41卷 圖1 部分柿及其近緣種cpdna pcr產(chǎn)物經(jīng)限制性?xún)?nèi)切酶酶切后2.5%瓊脂糖凝膠電泳圖譜 1 16見(jiàn)表1，ma為dl2000，mb為50 bp ladder。 fig. 1 restriction patterns of amplified and digested products of chloroplast genomic dna from diospyros spp. 116 refer to the numbers list in table 1. ma：dl2000；mb：50 bp ladder. 圖2 pcr-rflp分析的upgma聚類(lèi)圖 fig. 2 dendr

27、ogram of all genotypes from the upgma cluster analysis based on pcr-rflp data 12期 梁晉軍等：基于葉綠體dna pcr-rflp分析柿及其近緣種親緣關(guān)系 2479 圖3 柿屬植物種間pcr-rflp主坐標(biāo)分析散點(diǎn)圖 fig. 3 diagram showing the relationships among species based on principal coordinates analysis by pcr-rflp 3 討論 前人通過(guò)ty1-copia反轉(zhuǎn)錄轉(zhuǎn)座子southern雜交（nakasuka e

28、t al.，2002）；檢測(cè)dna重復(fù)位點(diǎn)序列同源性（choi et al.，2003）；srap、irap、remap（guo & luo，2006；guo et al.，2006）；cpdna pcr-rflp多態(tài)性分析（hu et al.，2008）均表明美洲柿與柿親緣關(guān)系較遠(yuǎn)，本試驗(yàn)結(jié)果與之一致。但keizo等（1998）通過(guò)對(duì)東南亞部分柿屬植物包括24個(gè)種的限制性?xún)?nèi)切酶酶切研究，用鄰接法分析表明柿與美洲柿、君遷子聚成一組，導(dǎo)致結(jié)果不完全相同的原因是所選材料和擴(kuò)增區(qū)段不同。 有研究用cpdna pcr-rflp和ssr分析表明柿與君遷子親緣關(guān)系最近（hu et al.，2008

29、；guo & luo，2011）；而本試驗(yàn)結(jié)果表明柿與君遷子、浙江柿親緣關(guān)系較遠(yuǎn)，柿與金棗柿、油柿、云南野毛柿聚為一類(lèi)，親緣關(guān)系較與君遷子和浙江柿更近（圖2，圖3），導(dǎo)致兩種不同結(jié)果最可能的原因是選用擴(kuò)增區(qū)域和限制性?xún)?nèi)切酶的不同。keizo等（2008）通過(guò)用nrdna its區(qū)和cpdna matk序列分析表明柿與油柿最近緣，這與本試驗(yàn)中柿與金棗柿、油柿、云南野毛柿聚為一類(lèi)的結(jié)果（圖2）相一致，就柿與油柿的葉片表觀形態(tài)而言，柿也與油柿相近。本試驗(yàn)結(jié)果表明柿種內(nèi)無(wú)酶切位點(diǎn)差異與前人（yuri & shozo，1994；keizo et al.，1998；hu et al.，20

30、08）用同一方法得到的結(jié)果一致。 金棗柿曾被學(xué)術(shù)界認(rèn)定為柿的一個(gè)二倍體品種（楊勇 等，1999），后經(jīng)hu等（2008）和tang等（2014）分別用cpdna pcr-rflp和its，matk序列變異分析證明金棗柿為一個(gè)新種。本試驗(yàn)中柿與金棗柿能完全分開(kāi)，支持金棗柿為一新種的結(jié)論。 references bao y，lu b r，ge s. 2005. identification of genomic constitutions of oryza species with the b and c genomes by the pcr-rflp method. genetic resour

31、ces and crop evolution，52：6976. choi y a，tao r，keizo y，sugiura a. 2003. genomic distribution of three repetitive dnas in cultivated hexaploid diospyros spp.（d. kaki and d. virginiana）and their wild relatives. genes genet syst，78：301308. claudia c，muhamment t，phillip d g. 2004. discrimination of dipl

32、oid crucifer species using pcr-rflp of chloroplast dna. hortscience，39 (7)：15751577. doyle j j，doyle j l. 1987. a rapid dna isolation procedure for small quantities of fresh leaf tissue. phytochem bul，19：1115. 2480 園 藝 學(xué) 報(bào) 41卷 guo d l，luo z r. 2006. genetic relationships of some pcna persimmons（dios

33、pyros kaki thunb.）from china and japan revealed by srap analysis. genet res crop evol，53：15971603. guo d l，luo z r. 2011. genetic relationships of the japanese persimmon diospyros kaki（ebenaceae）and related species revealed by ssr analysis. genet mol res，10 (2)：10601068. guo d l，zhang h q，luo z r. 2

34、006. genetic relationships of diospyros kaki thunb and related species revealed by irap and remap analysis. plant sci，170：528533. hu d c，zhang q l，luo z r. 2008. phylogenetic analysis in some diospyros spp.（ebenaceae）and japanese persimmon using chloroplast dna pcr-rflp markers. scientia horticultur

35、ae，117：3238. keizo y，chitose h，shinya k，hitofumi i，ayako l，akira k，akira s，dan e p. 2008. sequence analyses of the its regions and the matk gene for determining phylogenetic relationships of diospyros kaki（persimmon）with other wild diospyros（ebenaceae）species. tree genetics & genomes，(4)：149158.

36、 keizo y，shinya k，dan e p，naoki u，suranant s，akira s. 1998. phylogenetic relationship of diospyros kaki（persimmon）to diospyros spp. （ebenaceae）of thailand and four temperate zone diospyros spp. from an analysis of rflp variation in amplified cpdna. genome，41： 173182. kerstin d，trevor r h，ecelyn f，su

37、sanne b. 2008. an optimized chloroplast dna extraction protocol for grasses（poaceae）proves suitable for whole plastid genome sequencing and snp detection. plos one，3 (7)：e2813. nakasuka a，iwami n，matusmoto s，itamura h，yamagishi m. 2002. ty1-copia group retrotransposons in persimmon（diospyros kaki th

38、unb）. genes genet syst，77：131136. pharmawati m，yan g，sedgley r，finnegan p m. 2004. chloroplast dna inheritance and variation in leucadendron species（proteaceae）as revealed by pcr-rflp. theor appl genet，109：16941701. randall l s，julie a r，richard c c，tosak s，jonathan f w. 1998. the tortoise and the hare：choosing between non coding plastome and nuclear adh sequences for phylogeny reconstruction in a recently diverged plant group. ame j bot，85：13011315. tang d l，hu y，zhang q l，yang y，luo z r. 2014. discriminant analysis of“jinzaoshi”from persimmon