基于多序列比對(duì)的甜瓜SSR標(biāo)記開(kāi)發(fā)及應(yīng)用

王盼喬，周亞峰，許彥賓，胡建斌，楊路明，孫守如

(河南農(nóng)業(yè)大學(xué)園藝學(xué)院，河南 鄭州 450002)

基于多序列比對(duì)的甜瓜SSR標(biāo)記開(kāi)發(fā)及應(yīng)用

王盼喬，周亞峰，許彥賓，胡建斌，楊路明，孫守如

(河南農(nóng)業(yè)大學(xué)園藝學(xué)院，河南 鄭州 450002)

采用多序列比對(duì)的方法，從甜瓜EST序列中發(fā)現(xiàn)約1 100個(gè)SSR差異位點(diǎn)，并隨機(jī)挑選50個(gè)位點(diǎn)設(shè)計(jì)成SSR引物，對(duì)47份遺傳背景不同的甜瓜種質(zhì)進(jìn)行PCR擴(kuò)增。結(jié)果表明，41對(duì)引物可產(chǎn)生清晰的帶譜，有效擴(kuò)增率達(dá)82.00%，且均具有多態(tài)性。共檢測(cè)到130個(gè)等位基因，平均3.17個(gè)。Na、Ne、Ho、He、PIC等參數(shù)表明，這些標(biāo)記具有一定程度的多態(tài)性。本研究開(kāi)發(fā)的SSR標(biāo)記可用于甜瓜的遺傳多樣性和群體結(jié)構(gòu)分析，基于多序列比對(duì)的標(biāo)記開(kāi)發(fā)方法亦可在基因組信息豐富的物種中應(yīng)用。

甜瓜；序列比對(duì)；EST-SSR；遺傳多樣性；群體結(jié)構(gòu)

甜瓜(CucumismeloL.)是世界重要的園藝作物，其生產(chǎn)規(guī)模在世界10大水果中居第9位[1]。甜瓜是中國(guó)重要的設(shè)施蔬菜種類(lèi)之一，其生產(chǎn)規(guī)模僅次于西瓜，播種面積和產(chǎn)量歷年穩(wěn)居世界第一[2]。近年來(lái)，甜瓜基因組領(lǐng)域的研究進(jìn)展十分迅速，目前已經(jīng)構(gòu)建了至少20張遺傳圖譜，基因組信息不斷增加[3]。2007年，西班牙、美國(guó)、中國(guó)等14個(gè)實(shí)驗(yàn)室聯(lián)合成立了國(guó)際葫蘆科基因組計(jì)劃(International Cucurbit Genomics Initiative, ICuGI)，啟動(dòng)了包括甜瓜在內(nèi)的葫蘆科植物全基因組測(cè)序工作。2012年，甜瓜全基因組測(cè)序結(jié)果公開(kāi)發(fā)表[4]，是繼黃瓜、西瓜之后完成全基因組測(cè)序的葫蘆科植物。甜瓜基因組學(xué)的快速發(fā)展產(chǎn)生了大量的數(shù)據(jù)信息，目前GenBank中甜瓜的EST序列已經(jīng)超過(guò)13萬(wàn)條，為甜瓜的遺傳研究(如標(biāo)記開(kāi)發(fā))提供了豐富的數(shù)據(jù)信息。植物EST序列中一般含有SSR位點(diǎn)，其出現(xiàn)頻率一般為1/10～1/5[5]，這些SSR位點(diǎn)可用來(lái)開(kāi)發(fā)SSR標(biāo)記。通過(guò)SSR搜尋工具(如SSRHunter軟件)可將含有SSR的EST序列找出，在SSR位點(diǎn)側(cè)翼設(shè)計(jì)引物，便可開(kāi)發(fā)成SSR標(biāo)記。這種方法充分利用了現(xiàn)有的基因組信息，無(wú)需構(gòu)建文庫(kù)和測(cè)序，被認(rèn)為是一種“免費(fèi)”開(kāi)發(fā)SSR標(biāo)記的途徑，目前已經(jīng)應(yīng)用于多種植物[6]。然而，許多研究表明，通過(guò)這種方法開(kāi)發(fā)出的SSR標(biāo)記，其多態(tài)性標(biāo)記的比例不高。如白菜中多態(tài)性EST-SSR標(biāo)記的比例為46.7%[7]，蘿卜中多態(tài)性EST-SSR標(biāo)記的比例為35.0%[8]，黃瓜中多態(tài)性EST-SSR標(biāo)記比例僅為24.3%[9]。本研究采用多序列比對(duì)的方法開(kāi)發(fā)甜瓜SSR標(biāo)記，期望進(jìn)一步提高多態(tài)性標(biāo)記的比例，并對(duì)這些標(biāo)記的可用性進(jìn)行分析。

1 材料與方法

1.1材料

供試材料選自河南農(nóng)業(yè)大學(xué)園藝學(xué)院瓜類(lèi)育種實(shí)驗(yàn)室已構(gòu)建的甜瓜核心種質(zhì)庫(kù)[10]，共47份，其中厚皮種質(zhì)(subsp.melo) 28份，薄皮種質(zhì)(subsp.agrestis) 19份。所有材料均已培育成高代自交系，材料名稱(chēng)及來(lái)源見(jiàn)表1。

1.2方法

1.2.1 DNA提取 取各種質(zhì)幼葉約300 mg，采用CTAB法[11]提取基因組DNA，用無(wú)菌ddH2O稀釋至250 mg·L-1，備用。

1.2.2 SSR位點(diǎn)發(fā)掘及引物設(shè)計(jì) 從國(guó)際葫蘆科基因組(ICuGI)網(wǎng)站(http://www.icugi.org)下載含有SSR位點(diǎn)的Unigene(MU43177～MU67602)，摒棄現(xiàn)有報(bào)道中已用于開(kāi)發(fā)SSR標(biāo)記的序列，將剩余的Unigene與GenBank中甜瓜EST序列進(jìn)行同源性比對(duì)，選取適宜的序列進(jìn)行引物設(shè)計(jì)。選擇含有SSR的序列的條件為：1) 序列比對(duì)的相似性為85%以上;2) 至少有2條EST序列與該序列存在SSR位點(diǎn)差異。

采用Primer 5在符合條件的目標(biāo)序列(含有SSR差異位點(diǎn)的序列)側(cè)翼設(shè)計(jì)引物，引物設(shè)計(jì)的標(biāo)準(zhǔn)為：Tm 值為55～60 ℃，引物長(zhǎng)度為18～24 bp，PCR產(chǎn)物大小為100～350 bp，GC含量為40%～60%。引物由上海生工生物工程有限公司合成。

表1 試驗(yàn)材料和來(lái)源Table 1 Experimental materials and their origins

表2 試驗(yàn)所設(shè)計(jì)的50對(duì)引物的信息Table 2 Information of 50 primer pairs designed in the present study

續(xù)表Continuing Table

注：‘-’ 不能確定標(biāo)記的位置。

Note: ‘-’indicates that the marker is not located.

1.2.3 PCR擴(kuò)增及電泳檢測(cè) 以47份甜瓜種質(zhì)基因組DNA為模板，采用新合成的50對(duì)引物進(jìn)行PCR擴(kuò)增。擴(kuò)增體系為15 μL，其中2 μL 50 mg·L-1模板DNA，1.5 μL 10×PCR buffer (含Mg2+)，0.4 μL 2.5 mmol·L-1dNTPs，0.2 μL 2 500 U·L-1TaqDNA聚合酶，0.2 μL 10 μmol·L-1引物，8.6 μL ddH2O。PCR擴(kuò)增在Mastercycler PCR 儀上進(jìn)行，反應(yīng)程序?yàn)椋?4 ℃預(yù)變性4 min；94 ℃變性40 s，55～60 ℃退火40 s，72 ℃延伸1 min，共循環(huán)35次；最后72 ℃延伸8 min。采用本實(shí)驗(yàn)室建立的PAGE 凝膠制備方法和銀染的技術(shù)[12]，對(duì)PCR擴(kuò)增產(chǎn)物進(jìn)行檢測(cè)、記錄。

1.3SSR的多態(tài)性檢測(cè)與數(shù)據(jù)統(tǒng)計(jì)

根據(jù)電泳帶譜，在相同遷移位置有帶記作1，無(wú)帶記作0，建立基因分型的0/1矩陣。采用PopGene 1.32軟件計(jì)算等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、期望雜合度(He)和觀(guān)測(cè)雜合度(Ho)。根據(jù)BOTSTEIN等[13]的方法計(jì)算多態(tài)性信息量(Polymorphism information content, PIC)，

(1)

式中:Xi為第i種基因型出現(xiàn)的頻率。

采用NTSYSpc 2.10軟件進(jìn)行數(shù)據(jù)分析，計(jì)算Jaccard相似系數(shù)，采用UPGMA法對(duì)47份材料的基因型數(shù)據(jù)進(jìn)行聚類(lèi)分析，并通過(guò)Treeplot程序生成聚類(lèi)圖。為了估測(cè)供試種質(zhì)的群體遺傳結(jié)構(gòu)，應(yīng)用STRUCTURE 2.3.1軟件對(duì)多態(tài)性位點(diǎn)的分子數(shù)據(jù)進(jìn)行分析，計(jì)算各種質(zhì)的Q值(第i個(gè)材料基因組變異源于第K亞群的概率)。先設(shè)定群體數(shù)(K)為1～10，將MCMC(markov chain monte carlo)開(kāi)始時(shí)的不作數(shù)迭代(length of burn-in period)設(shè)為10 000次，再將不作數(shù)迭代后的MCMC設(shè)為100 000次，依據(jù)似然值最大的原則選取1個(gè)合適的K值，繪制群體遺傳結(jié)構(gòu)圖。

2 結(jié)果與分析

2.1SSR位點(diǎn)識(shí)別及擴(kuò)增效率

在ICuGI網(wǎng)站中發(fā)現(xiàn)有3 263個(gè)unigene至少含有1個(gè)SSR位點(diǎn)，約 70% 的unigene能在GenBank中找到同源序列。采用Clustalx軟件對(duì)unigene及其同源序列進(jìn)行多序列比對(duì)，結(jié)果發(fā)現(xiàn)約1 100個(gè)SSR位點(diǎn)在多序列比對(duì)中表現(xiàn)出位點(diǎn)差異。其中，unigene MU45323的序列比對(duì)結(jié)果見(jiàn)圖1，該序列的重復(fù)基元CTC和TCT與其他EST序列存在差異。本試驗(yàn)隨機(jī)選取了50個(gè)符合條件的差異SSR位點(diǎn)，設(shè)計(jì)50對(duì)EST-SSR引物(表2)，根據(jù)ICuGI網(wǎng)站公布的最新的甜瓜遺傳圖譜信息，將除HNM09和HNM27外的48個(gè)標(biāo)記定位在除LG3和LG11之外的所有連鎖群上。采用這些標(biāo)記對(duì)47份遺傳背景不同的甜瓜種質(zhì)進(jìn)行基因分型，結(jié)果發(fā)現(xiàn)，41對(duì)引物產(chǎn)生清晰可見(jiàn)的帶譜，引物的有效擴(kuò)增率為82.00%。在這41對(duì)引物中，4對(duì)引物(HNM01、HNM10、HNM35和HNM50)擴(kuò)增出非預(yù)期大小的片段，實(shí)際片段比期望長(zhǎng)度多出50 ～ 200 bp。編號(hào)為HNM27的引物同時(shí)檢測(cè)到2個(gè)SSR位點(diǎn)的變異，因而可分為HNM27a (～300 bp)和HNM27b (～ 200 bp) 2個(gè)標(biāo)記。剩余36對(duì)引物擴(kuò)增出預(yù)期產(chǎn)物大小的片段，其中引物HNM08的擴(kuò)增帶譜見(jiàn)圖2。

圖1 MU45323的多序列比對(duì)結(jié)果Fig.1 The result of multiple sequence alignment of MU45323

2.2SSR標(biāo)記的多態(tài)性驗(yàn)證

41對(duì)引物共檢測(cè)到130個(gè)等位基因，每對(duì)引物檢測(cè)到的等位基因數(shù)為2～9，平均3.167個(gè)(表3)。由于每對(duì)引物在47份種質(zhì)中均能檢測(cè)到等位變異，因而多態(tài)性引物的比例為100%。所有引物的有效等位基因數(shù)為1.403～5.111，HNM7最大，HNM49最小，平均2.410。期望雜合度為0.088～0.250，HNM01最大，HNM15最小，平均0.191。觀(guān)測(cè)雜合度為0.000～0.095，平均0.023。多態(tài)性信息量為0.246～0.775，HNM07最大，HNM49最小，平均0.470。這些參數(shù)說(shuō)明，供試種質(zhì)的純合度較高，新開(kāi)發(fā)的引物均具有一定程度的多態(tài)性。所有引物的擴(kuò)增產(chǎn)物呈“single-locus”變異，表型出典型的共顯性特征(圖2)，為電泳數(shù)據(jù)記錄提供了便利。

M：分子標(biāo)量marker pUC18 DNA/Msp1；泳道1～ 47為甜瓜種質(zhì)，序號(hào)同表1。 M: marker pUC18 DNA/Msp1; Lanes 1～ 47 mean the melon germplasms and their codes are the same to what is in table 1.

標(biāo)記名稱(chēng)Name等位基因數(shù)/個(gè)Allelenumber有效等位基因數(shù)/個(gè)Effectiveallelenumber期望雜合度Expectedheterozygosity觀(guān)測(cè)雜合度Observedheterozygosity多態(tài)性信息量PolymorphisminformationcontentHNM122.0000.2500.0210.375HNM242.9520.1650.0000.593HNM421.5410.1760.0000.290HNM521.9380.2420.0000.367HNM775.1110.1150.0650.775HNM821.6670.2000.0000.320HNM921.9960.2490.0210.375HNM1043.4050.1770.0230.653HNM1132.5960.2050.0430.533HNM1274.4490.1110.0950.742HNM1342.2890.1410.0000.503HNM1432.3180.1900.0000.477HNM1594.7300.0880.0000.761HNM1642.6010.1540.0650.549HNM1721.9130.2390.0210.363HNM1821.9770.2470.0440.372HNM1921.9570.2440.0000.369HNM2021.9570.2440.0000.369HNM2231.9530.1630.0430.388HNM2321.6580.1980.0000.318HNM2432.7360.2120.0430.558HNM2521.9960.2500.0210.375HNM2632.6650.2080.0000.547HNM27a42.4680.1490.0000.543HNM27b32.3490.1910.0210.502HNM2821.5840.1840.0000.301HNM3153.1050.1360.0640.629HNM3221.9410.2420.0000.367HNM3353.3730.1410.0460.646HNM3421.9760.2470.0000.372HNM3521.9860.2480.0210.373HNM3632.0060.1670.0000.412HNM3732.5880.2050.0640.533HNM3832.0620.1720.0210.459HNM3921.9570.2440.0000.369HNM4042.7740.1600.0000.589HNM4132.1460.1780.0480.469HNM4232.2020.1820.0250.484HNM4532.1200.1760.0510.441HNM4632.5580.2030.0430.526HNM4921.4030.1440.0440.246HNM5032.1980.1820.0000.486平均Mean3.1672.4100.1910.0230.470

2.3遺傳多樣性分析

為了進(jìn)一步評(píng)價(jià)所開(kāi)發(fā)的SSR標(biāo)記的可用性，對(duì)41對(duì)SSR引物在47份甜瓜種質(zhì)中的擴(kuò)增數(shù)據(jù)進(jìn)行聚類(lèi)分析。所有種質(zhì)的Jaccard相似系數(shù)為0.41～0.94，相似系數(shù)變化范圍較廣，說(shuō)明參試種質(zhì)遺傳背景差異較大。內(nèi)蒙古的‘金巴齒’和河南的‘牛角蜜’的相似系數(shù)最大，兩者在果實(shí)性狀上存在明顯的相似性；河北的薄皮種質(zhì)‘小翠瓜’和土耳其厚皮種質(zhì)‘PI-169358’相似系數(shù)最小，兩者在果實(shí)性狀和地理分布上差異明顯(表1)。

47份種質(zhì)在相似系數(shù)為0.52的水平上可明顯分為2類(lèi)(I和II) (圖3)。第I類(lèi)包括26份厚皮種質(zhì)，第II類(lèi)包括剩余的19份薄皮種質(zhì)和2份厚皮種質(zhì)。在第I類(lèi)中，不同地理來(lái)源地厚皮種質(zhì)聚在一起，說(shuō)明它們的遺傳關(guān)系緊密，可能具有相同的起源中心；但是河南的‘Horneseoneen’與其他種質(zhì)遺傳關(guān)系較遠(yuǎn)。 所有的薄皮種質(zhì)(中國(guó)17份，韓國(guó)和日本各1份)均聚在第II類(lèi)，平均相似系數(shù)高達(dá)0.78，其遺傳關(guān)系十分緊密。然而2份厚皮種質(zhì)(中國(guó)臺(tái)灣的‘PI-321004-S2’和伊朗的‘PI-137852’)也被劃分在第II類(lèi)，可能由于它們含有薄皮甜瓜的血緣，或因所用標(biāo)記數(shù)量的限制而不能將它們分開(kāi)；但在0.77的相似系數(shù)水平上，它們可與薄皮種質(zhì)完全區(qū)別開(kāi)。該研究結(jié)果表明，通過(guò)序列比對(duì)開(kāi)發(fā)的SSR標(biāo)記能較準(zhǔn)確地將甜瓜厚皮和薄皮種質(zhì)區(qū)分開(kāi)，可用于甜瓜種質(zhì)資源的遺傳多樣性研究。

圖3 基于SSR數(shù)據(jù)構(gòu)建的47份甜瓜種質(zhì)的聚類(lèi)圖Fig.3 Dendrogram of 47 melon germplasms based on the SSR data

2.4群體結(jié)構(gòu)分析

為了進(jìn)一步分析參試種質(zhì)的群體結(jié)構(gòu)，試驗(yàn)采用STRUCTURE軟件對(duì)47份種質(zhì)的基因型數(shù)據(jù)進(jìn)行基于Bayesian算法的群體結(jié)構(gòu)劃分，以推定的K所對(duì)應(yīng)的最大似然值LnP(D)來(lái)選擇合適的群體數(shù)目。由于LnP(D)隨K的增大持續(xù)增加，難以確定真實(shí)的K值，因此參照EVANNO等[14]的方法通過(guò)ΔK來(lái)確定K值。當(dāng)K=2時(shí)，ΔK出現(xiàn)明顯的峰值(圖4)，即參試種質(zhì)明顯存在2個(gè)亞群(圖5)。19份薄皮種質(zhì)和一份厚皮種質(zhì)(‘PI-321004-S2’)被劃分為亞群Pop1，27份厚皮種質(zhì)被劃分為亞群Pop2，基于Bayesian算法基本能將參試種質(zhì)按厚皮和薄皮亞種區(qū)分開(kāi)。供試種質(zhì)中，39份(95.12%)的種質(zhì)的Q值大于0.90，表明其血緣組成相對(duì)單一，2份厚皮種質(zhì)(‘PI-321004-S2’和‘PI-137852’)的Q值為0.30～0.70，具有混合血緣。本研究結(jié)果與UPGMA聚類(lèi)方法一致。

K,群體數(shù)目。K,The number of population.

圖5 47份甜瓜種質(zhì)的群體遺傳結(jié)構(gòu)Fig.5 Population structure of 47 melon germplasms

3 結(jié)論與討論

隨著新一代的測(cè)序技術(shù)快速發(fā)展，海量的序列信息(主要是EST序列)被釋放到公共數(shù)據(jù)庫(kù)中，成為開(kāi)發(fā)SSR、InDel、SNP等標(biāo)記的寶貴資源。由于SSR標(biāo)記具有共顯性、穩(wěn)定性、多等位變異等特點(diǎn)，是基于EST序列開(kāi)發(fā)的主要對(duì)象。目前，基于EST序列的SSR標(biāo)記開(kāi)發(fā)在植物[6]、動(dòng)物[15]、微生物[16]中均有報(bào)道，所開(kāi)發(fā)的EST-SSR標(biāo)記已廣泛用于遺傳多樣性分析[17]、遺傳圖譜構(gòu)建[18]、比較作圖[19]、基因定位[20]等領(lǐng)域。然而，眾多的研究表明，從EST序列中開(kāi)發(fā)的SSR標(biāo)記主要存在兩個(gè)問(wèn)題：一是標(biāo)記的有效擴(kuò)增率不高；二是標(biāo)記的多態(tài)性較低[21]。

從基因組文庫(kù)開(kāi)發(fā)SSR標(biāo)記是傳統(tǒng)方法，由于該方法是在基因組DNA上搜尋SSR位點(diǎn)，因而所開(kāi)發(fā)的標(biāo)記對(duì)基因組的擴(kuò)增是有效的，有效擴(kuò)增率可達(dá)100%[22]。而EST-SSR標(biāo)記是來(lái)自于基因組的編碼區(qū)，如果EST-SSR引物恰好落在2個(gè)外顯子上，或者上下游引物之間存在一個(gè)較長(zhǎng)的內(nèi)含子(超過(guò)了DNA聚合酶的擴(kuò)增長(zhǎng)度)，標(biāo)記則不會(huì)有擴(kuò)增產(chǎn)物。因此，EST-SSR標(biāo)記的有效擴(kuò)增率一般不會(huì)達(dá)到100%。本試驗(yàn)中，甜瓜SSR標(biāo)記的有效擴(kuò)增率為82.00%，略低于黃瓜的89.19%[9]和蘿卜的86.88%[8]，高于辣椒的60.44%[23]。隨著一些植物的全基因組序列的公開(kāi)釋放，標(biāo)記開(kāi)發(fā)者可先將目的EST序列與該物種的全基因組序列進(jìn)行比對(duì)，排除內(nèi)含子存在引物結(jié)合處或引物之間的情況，可提高標(biāo)記的擴(kuò)增效率。

從差異位點(diǎn)的檢測(cè)效率來(lái)講，EST-SSR的多態(tài)性一般低于基因組SSR。HU等[24]比較分析了黃瓜基因組SSR和EST-SSR的多態(tài)性，發(fā)現(xiàn)基因組SSR檢測(cè)到的等位基因數(shù)和PIC值均明顯高于EST-SSR。本試驗(yàn)開(kāi)發(fā)的41個(gè)SSR標(biāo)記平均檢測(cè)到的等位基因3.17個(gè)，PIC值0.470，而蘿卜的EST-SSR標(biāo)記僅能檢測(cè)2.5個(gè)等位基因[7]。與基因組的非編碼區(qū)相比，編碼區(qū)在植物進(jìn)化過(guò)程中承受的選擇壓力較大，其序列變異相對(duì)較小，因而EST-SSR標(biāo)記檢測(cè)到的等位變異的幾率不會(huì)很高[21]。TEMNYKH等[25]認(rèn)為SSR的多態(tài)性與重復(fù)基元的長(zhǎng)度有關(guān)，選擇長(zhǎng)度大于20 bp的SSR可提高標(biāo)記的多態(tài)性。但在百合[26]、甜瓜[27]、刺梨[28]等植物中，一些長(zhǎng)度超過(guò)20 bp的EST-SSR位點(diǎn)也不具有多態(tài)性，本研究中一些低重復(fù)的SSR (長(zhǎng)度小于20 bp)卻在不同種質(zhì)間表現(xiàn)等位變異。本研究是針對(duì)具有SSR差異的序列(簇)開(kāi)發(fā)標(biāo)記，從源頭上保證了標(biāo)記的多態(tài)性潛能，在以遺傳背景差異明顯的材料作為試材的前提下，表現(xiàn)多態(tài)性的標(biāo)記比例可達(dá)到100%。而針對(duì)單個(gè)EST序列開(kāi)發(fā)的SSR標(biāo)記，如在上述蘿卜、黃瓜、辣椒等植物中，并非全部具有多態(tài)性。因此，在基因組信息豐富的物種中，通過(guò)序列比對(duì)法開(kāi)發(fā)SSR標(biāo)記，其優(yōu)勢(shì)十分明顯。目前，基因組測(cè)序的成本越來(lái)越低，一些植物還完成了不同遺傳背景材料的重測(cè)序。在這種情況下，電子PCR (Electronic PCR, e-PCR)技術(shù)的優(yōu)勢(shì)十分突出。它通過(guò)借助相應(yīng)的序列分析軟件，將設(shè)計(jì)的PCR引物與所查詢(xún)樣品的序列進(jìn)行比對(duì)，以評(píng)價(jià)引物的擴(kuò)增效率和多態(tài)性，可用于PCR 引物輔助篩選[29]。

此外，本研究還將新開(kāi)發(fā)的SSR標(biāo)記用于47份甜瓜種質(zhì)的遺傳多樣性和群體結(jié)構(gòu)分析，結(jié)果發(fā)現(xiàn)，基于新開(kāi)發(fā)標(biāo)記的聚類(lèi)結(jié)果和群體劃分與甜瓜分類(lèi)基本一致。由于薄皮甜瓜主要分布在東亞地區(qū)(中國(guó)中、東部，朝鮮半島和日本)，厚皮種質(zhì)主要分布在中亞、西亞、歐洲和美洲[30]，因此聚類(lèi)分析和群體劃分的結(jié)果也與甜瓜的地理分布基本相符，說(shuō)明所開(kāi)發(fā)的標(biāo)記可用于甜瓜遺傳多樣性和群體結(jié)構(gòu)研究。

本研究采用多序列比對(duì)的方法從甜瓜EST序列中發(fā)現(xiàn)約1 100個(gè)SSR差異位點(diǎn)，并開(kāi)發(fā)了50對(duì)SSR引物，48個(gè)標(biāo)記分布在除LG3和LG11之外的其他連鎖群。對(duì)47份遺傳背景不同的甜瓜種質(zhì)進(jìn)行分子標(biāo)記分析，41對(duì)引物產(chǎn)生清晰的帶譜，且均具有多態(tài)性，共檢測(cè)到130個(gè)等位基因。Na、Ne、Ho、He、PIC等遺傳參數(shù)顯示，這些標(biāo)記具有一定程度的多態(tài)性。UPGMA聚類(lèi)方法和Bayesian群體劃分方法均能按照甜瓜植物學(xué)分類(lèi)和地理分布將所有種質(zhì)分為2類(lèi)或亞群。本研究開(kāi)發(fā)的甜瓜SSR標(biāo)記不僅進(jìn)一步豐富了甜瓜DNA標(biāo)記資源，而且對(duì)于基因組豐富物種的標(biāo)記開(kāi)發(fā)，也具有十分重要的啟示意義。

[1] 劉君璞, 許勇, 孫小武, 等. 我國(guó)西瓜甜瓜產(chǎn)業(yè)“十一五”的展望及建議[J]. 中國(guó)瓜菜, 2006, 19(1): 1-3.

[2] 王志丹. 中國(guó)甜瓜產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展研究[D].北京：中國(guó)農(nóng)業(yè)科學(xué)院, 2014.

[3] ARGYRIS J M, PUJOL M. MARTIN-HERNANDEZ A M, et al. Combined use of genetic and genomics resources to understand virus resistance and fruit quality traits in melon[J]. Physiologia plantarum, 2015, 155(1): 4-11.

[4] GARCIA-MAS J, BENJAK A, SANSEVERINO W, et al. The genome of melon (CucumismeloL.) [J]. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(29): 11872-11877.

[5] GAO L, TANG J, LI H, et al. Analysis of microsatellites in major crops assessed by computational and experimental approaches[J]. Molecular Breeding, 2003, 12(3): 245-261.

[6] KALIA R K, RAI M K, KALIA S, et al. Microsatellite markers: an overview of the recent progress in plants[J]. Euphytica, 2011, 177(3): 309-334.

[7] 忻雅, 崔海瑞, 盧美貞, 等. 白菜EST-SSR信息分析與標(biāo)記的建立[J]. 園藝學(xué)報(bào),2006, 33(3): 549-554.

[8] 崔娜, 邱楊, 李錫香, 等. 蘿卜EST資源的SSR信息分析及EST-SSRs標(biāo)記開(kāi)發(fā)[J]. 園藝學(xué)報(bào), 2012, 39(7): 1303-1312.

[9] 關(guān)媛, 李效尊, 潘俊松, 等. 黃瓜果實(shí)EST-SSR標(biāo)記的開(kāi)發(fā)[J]. 分子植物育種, 2007, 5(5): 725-728.

[10] 胡建斌, 馬雙武, 王吉明, 等. 基于表型性狀的甜瓜核心種質(zhì)構(gòu)建[J]. 果樹(shù)學(xué)報(bào), 2013, 30(3): 404-411.

[11] PATERSON A H, BRUBAKER C L, WENDEL J F. A rapid method for extraction of cotton (Gossypium spp.) genomic DNA suitable for RFLP or PCR analysis[J]. Plant Molecular Biology Reporter, 1993, 11(2): 122-127.

[12] 胡建斌, 李靜, 袁鳴, 等. PAGE凝膠制備方法和銀染方法的改進(jìn)[J]. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào), 2009, 31(4): 742-745.

[13] BOTSTEIN D, WHITE R L, SKOLNICK M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. American journal of human genetics, 1980, 32(3): 314-331.

[14] EVANNO G, REGNAUT S, GOUDET J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study[J]. Molecular Ecology, 2005, 14(8): 2611-2620.

[15] 王解香, 于凌云, 白俊杰, 等. 草魚(yú)EST-SSR標(biāo)記及5個(gè)不同地域群體的遺傳結(jié)構(gòu)分析[J]. 動(dòng)物學(xué)雜志, 2011, 46(5): 24-32.

[16] 朱振東, 霍云龍, 王曉鳴, 等. 從大豆疫霉菌ESTs中發(fā)掘SSR標(biāo)記[J]. 科學(xué)通報(bào), 2004, 49(17): 1749-1754.

[17] 劉俊青, 李佩芳, 胡建斌, 等. 用EST-SSR標(biāo)記評(píng)價(jià)21份黃瓜種質(zhì)的遺傳多樣性[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 45(6): 657-661.

[18] 陳海梅, 李林志, 衛(wèi)憲云, 等. 小麥EST-SSR標(biāo)記的開(kāi)發(fā)、染色體定位和遺傳作圖[J]. 科學(xué)通報(bào), 2005, 50(20): 2208-2216.

[19] YU J K, DAKE T M, SINGH S, et al. Development and mapping of EST-derived simple sequence repeat markers for hexaploid wheat[J].Genome, 2004,47(5): 805-818.

[20] 張寧, 張顯, 張勇, 等. 甜瓜果實(shí)糖含量相關(guān)性狀QTL分析[J]. 西北植物學(xué)報(bào), 2015, 35(2): 252-257.

[21] VARSHNEY R K, GRANER A, SORRELLS M E. Genic microsatellite markers in plants: features and applications[J]. Trends in Biotechnology, 2005, 23(1): 48-55.

[22] 張?jiān)龃? 侯喜林. SSR分子標(biāo)記開(kāi)發(fā)策略及評(píng)價(jià)[J]. 遺傳, 2004, 26(5): 763-768.

[23] 劉峰, 王運(yùn)生, 田雪亮, 等. 辣椒轉(zhuǎn)錄組SSR挖掘及其多態(tài)性分析[J]. 園藝學(xué)報(bào), 2012, 39(1): 168-174.

[24] HU J, WANG L, LI J. Comparison of genomic SSR and EST-SSR markers for estimating genetic diversity in cucumber[J]. Biologia Plantarum, 2011, 55(3): 577-580.

[25] TEMNYKH S. Computational and Experimental Analysis of Microsatellites in Rice (Oryzasativa L.): Frequency, Length Variation, Transposon Associations, and Genetic Marker Potential [J]. Genome Research, 2001, 11(8): 1441-1452.

[26] 楊素麗, 明軍, 劉春, 等. 基于EST信息的百合SSR標(biāo)記的建立[J]. 園藝學(xué)報(bào), 2008, 35(7): 1069-1074.

[27] 胡建斌, 李建吾. 甜瓜EST-SSR位點(diǎn)信息及標(biāo)記開(kāi)發(fā)[J]. 園藝學(xué)報(bào), 2009, 36(4): 513-520.

[28] 鄢秀芹, 魯敏, 安華明. 刺梨轉(zhuǎn)錄組SSR信息分析及其分子標(biāo)記開(kāi)發(fā)[J]. 園藝學(xué)報(bào), 2015, 42(2): 341-349.

[29] 張潔, 陸海峰, 李有志. 電子PCR [J]. 分子植物育種, 2004, 2(1): 139-145.

[30] 林德佩. 中國(guó)栽培甜瓜植物的起源、分類(lèi)及進(jìn)化[J]. 中國(guó)瓜菜, 2010, 23(4): 34-36.

(責(zé)任編輯：李瑩)

DevelopmentofSSRmarkersbasedonmultiplesequencealignmentandtheirapplicationsinmelon

WANG Panqiao, ZHOU Yafeng, XU Yanbin, HU Jianbin, YANG Luming, SUN Shouru

(College of Horticulture, Henan Agricultural University, Zhengzhou 450002, China)

Multiple sequence alignment was used to develop SSR markers in melon, which identified about 1 100 different SSR loci from the melon ESTs. Of these loci, 50 were randomly selected for primer design and the next PCR amplification. The genomic DNA of 47 melon germplasms with different genetic backgrounds was analyzed with these primers. The results showed that 41 primer pairs generated clear band patterns, with the amplification efficiency of 82.00%, and all the products were polymorphic. A total of 130 alleles were detected with a mean of 3.17. The genetic parameters, Na, Ne, Ho, He, and PIC, verified a certain level of polymorphism of these newly-developed markers. Therefore, the SSR markers developed herein can be used to analyze genetic diversity and population structure in melon. Also, the method of SSR development based on multiple sequence alignment can be well applied in the species with abundant genomic information.

melon; sequence alignment; EST-SSR; genetic diversity; population structure

S652

：A

2015-09-07

國(guó)家自然科學(xué)基金項(xiàng)目(31101544)；河南省高等學(xué)校重點(diǎn)科研項(xiàng)目(15A210011)

王盼喬(1990-) ，男，河南郟縣人，碩士研究生，從事甜瓜育種研究。

胡建斌(1976-) ，男，湖北武漢人，副教授，博士。

1000-2340(2016)02-0189-09