基于瓠瓜轉(zhuǎn)錄組測序的EST-SSR標(biāo)記的開發(fā)及其應(yīng)用

許端祥 杜文麗 陳中釤 趙瑞麗 徐同偉 高山

摘? 要：為豐富瓠瓜分子標(biāo)記庫、開發(fā)瓠瓜SSR標(biāo)記、尋找有效的瓠瓜種質(zhì)資源鑒定方法，利用MISA軟件對(duì)瓠瓜轉(zhuǎn)錄組測序獲得的87 518條unigene序列進(jìn)行篩選，共檢測出11 029個(gè)SSR位點(diǎn)，發(fā)生頻率為9.85%，分布距離為平均每8.3 kb有1個(gè)SSR 位點(diǎn)。其中只包含1個(gè)SSR位點(diǎn)的unigene 序列有6759條，其發(fā)生頻率為7.72%;有1858條unigene序列含有2個(gè)及2個(gè)以上的SSR位點(diǎn)，發(fā)生頻率為2.123%;有920條unigene序列屬于混合形式的SSR標(biāo)記，頻率為1.051%。瓠瓜轉(zhuǎn)錄組中優(yōu)勢重復(fù)基序?yàn)閱魏塑账?、三核苷酸和二核苷酸，分別占總SSR位點(diǎn)的55.51%、25.41%和17.07%。A/T、AG/CT和AAG/CTT分別是單核苷酸、二核苷酸和三核苷酸的優(yōu)勢重復(fù)基元，分別占總SSR重復(fù)類型的98.22%、55.39%和39.86%。用Primer 3.0共設(shè)計(jì)出8617對(duì)SSR引物，利用6對(duì)條帶清晰、多態(tài)性穩(wěn)定的SSR引物構(gòu)建了36個(gè)品種的指紋圖譜，部分位點(diǎn)具有高度的一致性，說明基因來源范圍較小。聚類分析顯示，在相似系數(shù)為0.68時(shí)將36個(gè)瓠瓜品種分為6類。其表明我國瓠瓜品種資源遺傳多樣性非常狹窄，利用現(xiàn)有瓠瓜資源開展品種選育潛力有限，應(yīng)加強(qiáng)對(duì)野生特異種質(zhì)的引進(jìn)與發(fā)掘利用。利用瓠瓜轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行SSR標(biāo)記開發(fā)，獲得的SSR位點(diǎn)能為其遺傳多樣性分析和遺傳圖譜構(gòu)建提供更豐富可靠的標(biāo)記選擇。

關(guān)鍵詞：瓠瓜;轉(zhuǎn)錄組;EST-SSR;指紋圖譜;遺傳多樣性

中圖分類號(hào)：S188+.1? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A

EST-SSR Marker Development and Utilization Based on Transcriptome of Lagenariasiceraria （Mol.） Stand.

XU Duanxiang， DU Wenli， CHEN Zhongshan， ZHAO Ruili， XU Tongwei， GAO Shan*

Fuzhou Institute of Vegetable Science， Fuzhou， Fujian 350111， China

Abstract： In order to enrich the molecular marker library， the SSR markers of Lagenariasiceraria （Mol.） Stand. 87 518 unigene sequences from bottle gourd transcriptome database were developed by the MISA software. A total of 11 029 SSR loci with an occurrence frequency of 9.85% were obtained from the transcriptome of Lagenariasiceraria （Mol.）， the distribution distance of the SSRs was about 1 SSR per 8.3 kb ESTs. Around 6759 unigene sequences contained a single SSR loci， and the occurrence frequency was 7.72%. Average of 1858 unigene contained 2 or more than 2 SSR loci， and the occurrence frequency was 2.123%. 920 unigene contained the mixed SSR loci， and the occurrence frequency was 1.051%. Among the SSR locis， mononucleotide， trinucleotide and dinucleotide were the major types， accounting for 55.51%， 25.41% and 17.07% of the total， respectively. A/T， AG/CT and AAG/CTT were the dominant elements among mononucleotide， dinucleotide and trinucleotide， accounting for 98.22%， 55.39% and 39.86% of the total SSRs. 8617 pairs of SSR primers were found by Primer 3.0， six primer pairs were selected and applied in PCR amplification of 36 varieties， based on the clear and repeatable polymorphism， the fingerprint of the cultivars was constructed， however， there were no differences in some loci probably due to narrow genes. The cluster analysis showed that 36 bottle gourd varieties were divided into six categories when the genetic similarity coefficient was 0.68. The results aslo showed that the genetic diversity basis of bottle gourd germplasm resources in China was very low， limiting the potential of variety development using the existing bottle gourd germplasm resources. The introduction and exploration of specific wild bottle gourd germplasm should be highly emphasized. The results indicated that SSR markers could be acquired in bottle gourd by transcriptome sequencing analysis and abundant SSR markers would provide more reliable markers for genetic diversity analysis and genetic mapping construction in bottle gourd.

Keywords： bottle gourd; transcriptome; EST-SSR; fingerprint; genetic diversity

DOI： 10.3969/j.issn.1000-2561.2020.03.015

瓠瓜[Lagenaria siceraria （Mol.） Stand.]隸屬于葫蘆科葫蘆屬1年生草本植物，是世界上重要的葫蘆科作物之一。瓠瓜原產(chǎn)于非洲南部低地，現(xiàn)廣泛分布于非洲、印度、東南亞等地區(qū)。我國是世界上栽培瓠瓜歷史最悠久的國家之一，目前廣泛栽培于長江流域及以南地區(qū)，現(xiàn)已成為我國南方地區(qū)春、夏季重要的瓜類蔬菜作物之一，在北方地區(qū)栽培面積也日趨增大[1-2]。瓠瓜種質(zhì)資源眾多，但各品種間親緣關(guān)系的鑒定和分類工作尚未形成體系，僅在形態(tài)上的劃分不足以揭示瓠瓜種質(zhì)資源的遺傳多樣性，盡管瓠瓜在生產(chǎn)上占有重要地位，但其遺傳學(xué)基礎(chǔ)研究卻遠(yuǎn)遠(yuǎn)落后于黃瓜、西瓜、甜瓜等其他瓜類作物，目前在分子水平上檢測瓠瓜不同種質(zhì)間遺傳多樣性的研究還比較欠缺，嚴(yán)重阻礙了育種工作的進(jìn)行。

近幾年隨著分子生物學(xué)技術(shù)飛速發(fā)展，作為重復(fù)性好、特異性強(qiáng)和多態(tài)性高的共顯性遺傳標(biāo)記的微衛(wèi)星序列即簡單重復(fù)序列（simple sequence repeat， SSR）分子標(biāo)記早已成為植物品種鑒定、親緣關(guān)系鑒定和遺傳圖譜構(gòu)建等育種工作的有力手段。通過構(gòu)建基因組文庫、核酸數(shù)據(jù)庫篩選等傳統(tǒng)方法獲得SSR引物存在的工作量大、成功率低、多態(tài)性差等諸多問題，隨著第2代高通量測序技術(shù)的迅速發(fā)展，利用轉(zhuǎn)錄組測序的表達(dá)序列獲得EST-SSR標(biāo)記比基于基因組序列開發(fā)進(jìn)而獲得的SSR標(biāo)記的方法具有更高的種間通用性，使SSR位點(diǎn)克隆篩選的工作效率大大提高[3]。迄今，在番茄[4]、黃瓜[5]、黃秋葵[6]等一些蔬菜作物中均已有轉(zhuǎn)錄組數(shù)據(jù)庫成功開發(fā)出EST-SSR標(biāo)記并應(yīng)用于遺傳多樣性分析及指紋圖譜構(gòu)建的報(bào)道。瓠瓜的分類并沒有其他主要糧食作物那樣有著較明顯的劃分，由于人為選擇和長期的引種等原因，育種工作者只是根據(jù)不同的瓠瓜種質(zhì)資源的表型特性進(jìn)行粗略歸類，凸顯了在分子水平上檢測瓠瓜遺傳多樣性的重要性。瓠瓜的基因組序列信息和位點(diǎn)特異性分子標(biāo)記（SSR、SNP等）長期以來少有報(bào)道，僅見王莎[1]利用SSR分子標(biāo)記對(duì)瓠瓜種質(zhì)資源的遺傳多樣性和親緣關(guān)系進(jìn)行研究，因此，有必要在深度和廣度上進(jìn)行拓展。

本研究應(yīng)用Illumina高通量測序技術(shù)獲得豐富的瓠瓜轉(zhuǎn)錄組數(shù)據(jù)，通過搜索unigene中的SSR位點(diǎn)信息，分析其分布及組成特征并進(jìn)行可用性評(píng)價(jià)，針對(duì)SSR位點(diǎn)的2側(cè)序列開發(fā)設(shè)計(jì)了大量SSR引物，并以福州市蔬菜科學(xué)研究所種質(zhì)資源庫保存的具有不同優(yōu)良性狀的共36個(gè)瓠瓜品種進(jìn)行引物的有效性驗(yàn)證，進(jìn)而對(duì)36個(gè)瓠瓜品種進(jìn)行了DNA指紋圖譜的構(gòu)建和親緣關(guān)系分析，以期為瓠瓜種質(zhì)鑒定、親緣關(guān)系分析、基因的定位及分子標(biāo)記輔助選擇育種提供理論支撐，也為我國瓠瓜遺傳改良、育種親本選配、種質(zhì)保護(hù)、核心種質(zhì)庫構(gòu)建等提供科學(xué)依據(jù)。

1? 材料與方法

1.1? 轉(zhuǎn)錄組數(shù)據(jù)來源

從福州市蔬菜科學(xué)研究所瓠瓜種質(zhì)資源庫中選取地方特有品種福州芋瓠作為轉(zhuǎn)錄測序材料。于2018年4月中旬播種，5月上旬定植于福州市蔬菜科學(xué)研究所基地，于花后6 d剪去嫩葉轉(zhuǎn)移至液氮速凍，?80 ℃保存?zhèn)溆?。測序時(shí)每個(gè)單株構(gòu)建1個(gè)文庫，委托上海凌恩生物科技有限公司采用Illumina TruseqTM RNA sample prep Kit方法構(gòu)建文庫進(jìn)行無參轉(zhuǎn)錄組分析，經(jīng)過濾后得到高質(zhì)量的測序數(shù)據(jù)525 959 048個(gè)reads片段，共計(jì)77 337 803 108 bp（77.34 Gb）個(gè)核苷酸，采用Trinity進(jìn)行序列組裝產(chǎn)生了87 518條unigene，平均長度為1027.79 bp。76.03%的unigene長度主要集中在1～1000 bp，堿基Q20大于98%，Q30大于94%，GC含量為56.54%。該測序結(jié)果可用于后續(xù)分析。

1.2? 材料來源及其DNA提取

用于SSR引物篩選和可用性評(píng)價(jià)的36個(gè)瓠瓜品種均來源于福州蔬菜科學(xué)研究所中心資源庫（表1），包括了31個(gè)國內(nèi)品種及引自泰國5個(gè)品種，其均代表目前市場上栽培的主栽品種。采用CTAB法對(duì)瓠瓜基因組DNA進(jìn)行提取。

1.3? 轉(zhuǎn)錄組SSR位點(diǎn)鑒別及SSR引物設(shè)計(jì)

使用MISA 程序（http：//pgrc.ikp-gaters le ben.de/misa）搜索瓠瓜SSR 位點(diǎn)搜索，包括2、3、4、5、6核苷酸最少重復(fù)次數(shù)分別為6、5、5、4和4次為判斷標(biāo)準(zhǔn)進(jìn)行搜索，同時(shí)篩選復(fù)合型SSR，篩選標(biāo)準(zhǔn)為打斷SSR的堿基間隔小于或等于100 bp。使用Primer 3.0軟件對(duì)SSR重復(fù)單元前后的序列進(jìn)行評(píng)價(jià)并分別設(shè)計(jì)3條引物。引物序列長度18～27 bp，GC含量為40%～60%，退火溫度57～63 ℃，上、下游引物的Tm值相差≤2 ℃，預(yù)期擴(kuò)增產(chǎn)物長度100～500 bp，且無2級(jí)結(jié)構(gòu)和2聚體。

1.4? EST-SSR 引物篩選和指紋圖譜代碼構(gòu)建

隨機(jī)挑選并合成50對(duì)引物（華大生物科技有限公司）。PCR反應(yīng)體系為20 μL，其中10 mmol/L dNTPs 0.4 μL，5 U Taq酶0.3 μL，100 ng DNA 1.5 μL，10 μmol /L的上、下游引物各1 μL，10× Buffer（Mg2+）2.5 μL，加dd H2O補(bǔ)至20 μL（購自大連寶生生物有限公司）。PCR擴(kuò)增程序?yàn)椋?5 ℃預(yù)變性5 min;37個(gè)循環(huán)，每個(gè)循環(huán)包括94 ℃變性30 s，60 ℃退火30 s（退火溫度因不同引物而異），72 ℃延伸1 min;最后72 ℃延伸7 min。PCR擴(kuò)增產(chǎn)物在1%瓊脂糖凝膠電泳初篩，再用6%變性聚丙烯酰胺凝膠檢測，160 V電壓，3 h，銀染顯色后觀察拍照。指紋圖譜代碼構(gòu)建方法借鑒馬紅勃等[7]及段艷鳳等[8]的方法，對(duì)清晰、可重復(fù)的條帶進(jìn)行記錄。按同一位點(diǎn)條帶有無分別賦值，有帶記為“1”，無帶記為“0”，模板缺失記為“9”，指紋圖譜的構(gòu)建按照多態(tài)性條帶分子量從大到小依次進(jìn)行記錄。

1.5? 數(shù)據(jù)處理

SSR 發(fā)生頻率為SSR的unigene數(shù)量與總unigene數(shù)量之比，SSR的出現(xiàn)頻率為SSR的個(gè)數(shù)與總unigene的數(shù)量比，SSR平均分布距離為1 kb以上的unigene的堿基數(shù)與SSR數(shù)量之比。指紋圖譜代碼統(tǒng)計(jì)結(jié)果錄入Excel軟件，利用Power Marker軟件計(jì)算PIC值，通過 NTSYSpc 2.10e 進(jìn)行種質(zhì)聚類分析。

2? 結(jié)果與分析

2.1? 瓠瓜轉(zhuǎn)錄組中SSR的分布及結(jié)構(gòu)特點(diǎn)

瓠瓜轉(zhuǎn)錄組經(jīng)測序組裝共產(chǎn)生了87 518條unigene，平均長度為1027.79 bp。用MISA軟件對(duì)1 kb以上的unigene進(jìn)行檢索發(fā)現(xiàn)8617條unigene序列中存在11 029個(gè)SSR位點(diǎn)，SSR位點(diǎn)在瓠瓜轉(zhuǎn)錄組中總的發(fā)生頻率為9.846%，平均每8.3 Kb出現(xiàn)1個(gè)SSR;其中只包含1個(gè)SSR位點(diǎn)的unigene序列有6759條，其發(fā)生頻率為7.72%，有1858條unigene序列含有2個(gè)及2個(gè)以上的SSR位點(diǎn)，發(fā)生頻率為2.123%，有920條unigene序列屬于混合形式的SSR標(biāo)記，頻率為1.051%。

SSR重復(fù)類型包含2核苷酸至6核苷酸重復(fù)5大類。其中3核苷酸比例最占優(yōu)勢，占總SSR的65.02%，其次是2核苷酸占總SSR的31.29%;4核苷酸、5核苷酸和6核苷酸重復(fù)類型數(shù)量較少，分別只占總數(shù)的2.24%、0.75%和0.70%（表2）。

瓠瓜轉(zhuǎn)錄組SSR重復(fù)單元的重復(fù)次數(shù)分布在5～28次之間，其中5～10次的SSR共有5582個(gè)，占總數(shù)的92.77%;其次為11～28次的SSR共有435個(gè)，占總數(shù)的7.23%。瓠瓜轉(zhuǎn)錄組SSR的長度集中在11～311 bp之間，其中長度在11～20 bp的SSR達(dá)6621個(gè);占總數(shù)的65.55%，長度在20 bp以上的SSR有3488個(gè)，占34.50%。

2.2? 瓠瓜轉(zhuǎn)錄組中SSR基序重復(fù)類型和頻率特征

瓠瓜轉(zhuǎn)錄組SSR種類豐富，6種SSR均有出現(xiàn)，但出現(xiàn)頻率不同，其中3核苷酸最多，2核苷酸次之，4核苷酸、5核苷酸、6核苷酸重復(fù)的數(shù)量較少。瓠瓜轉(zhuǎn)錄組中4907個(gè)SSR位點(diǎn)共含170種重復(fù)基序，2、3、4、5、6核苷酸重復(fù)分別有4、10、25、43和99種。從分布頻率來看（表3），以3核苷酸重復(fù)類型AAG/CTT出現(xiàn)最多，占總SSR的22.76%，占3核苷酸重復(fù)基序的總數(shù)的39.86%。其次是2核苷酸重復(fù)類型AG/CT，占總SSR的21.26%，占2核苷酸總數(shù)的55.39%。4核苷酸、5核苷酸尤其是6核苷酸重復(fù)基序分布較為分散，出現(xiàn)的次數(shù)在40以下，所占比例極低。

2.3? 瓠瓜轉(zhuǎn)錄組SSR引物設(shè)計(jì)、篩選及多態(tài)性分析

利用Primer 3.0軟件對(duì)含SSR位點(diǎn)的8617條unigene序列進(jìn)行引物設(shè)計(jì)，共設(shè)計(jì)出8088對(duì)引物。隨機(jī)挑選50對(duì)不同重復(fù)單元（2、3、4、5、6核苷酸）的引物合成后對(duì)本所收集的36份瓠瓜DNA混樣進(jìn)行SSR-PCR擴(kuò)增，以驗(yàn)證其有效性。根據(jù)引物設(shè)計(jì)要求，無帶和條帶高于750 bp的引物都不能進(jìn)行后續(xù)實(shí)驗(yàn)，50對(duì)引物中有4對(duì)引物無帶或者帶不清晰，其余46對(duì)PCR擴(kuò)增產(chǎn)物與預(yù)期大小符合，有效擴(kuò)增率達(dá)到92%，說明引物設(shè)計(jì)效果基本合格。在46對(duì)有效擴(kuò)增引物中經(jīng)過多次試驗(yàn)確定最適退火溫度，篩選出6對(duì)多態(tài)性引物能將36份瓠瓜品種區(qū)分開，占有效擴(kuò)增率的13%;共獲得86條條帶，其中多態(tài)性位點(diǎn)為59個(gè)，多態(tài)性比率為68.6%（圖1）;其中引物L(fēng)sSSR25、LsSSR424和LsSSR3擴(kuò)增位點(diǎn)數(shù)最多，分別為13、11和10個(gè)（表4）;每對(duì)引物有效多態(tài)性位點(diǎn)數(shù)在8～13個(gè)，每對(duì)引物平均擴(kuò)增出9.8多態(tài)性位點(diǎn);多態(tài)信息含量PIC值為0.299～0.839，其中YHSSR42的PIC值達(dá)到0.87，平均多態(tài)信息含量為0.513，說明選擇的引物適用于試驗(yàn)瓠瓜樣品的指紋圖譜構(gòu)建。

2.4? DNA數(shù)字指紋圖譜構(gòu)建

根據(jù)篩選的6對(duì)SSR引物在36個(gè)瓠瓜品種的擴(kuò)增結(jié)果，構(gòu)建出36個(gè)品種資源的DNA數(shù)字

指紋圖譜（表5），可見部分位點(diǎn)在多數(shù)品種中呈現(xiàn)較高的一致性，說明36份瓠瓜的基因來源范圍較窄。但36個(gè)品種的DNA指紋圖譜均獨(dú)一無二，因此可用于瓠瓜品種資源的分類與鑒定。

2.5? 利用SSR標(biāo)記對(duì)瓠瓜品種資源聚類分析

根據(jù)6對(duì)引物擴(kuò)增的59個(gè)多態(tài)性位點(diǎn)，對(duì)36份瓠瓜供試品種資源進(jìn)行了聚類分析，其遺傳相似系數(shù)變幅為0.60～0.97。在遺傳距離0.60處，36份瓠瓜劃分成2大類，在遺傳距離0.68處劃分為6類（圖2）：第1類包含有19份瓠瓜品種資源，按瓜皮色種質(zhì)區(qū)分，除2份白綠色外其余均為綠色;其中長棒形的有6份，棒形的有8份，短棒形的有5份;6份來自中國福建，6份來自中國湖北，3份來自中國廣東，2份來自中國浙江，還有2份來自泰國。第2類有來源于中國湖北、中國廣東和泰國的3份瓠瓜品種資源;其中長棒形1份，短棒形2份;2份綠瓜皮色，1份綠白皮色。第3類有4份瓠瓜種質(zhì)，瓜型、瓜皮色一致，且品種資源均來自中國福建。第4類包含來自福建和泰國2份瓠瓜品種資源，瓜型均為長棒形，瓜皮色分別為亮綠和綠色。第5類包含7份瓠瓜品種資源，其中長棒形1份，棒形4份，短棒形2份;4份瓜皮綠色，2份瓜皮綠白色，1份瓜皮亮綠色;來自福建和廣東的各有3份，1份來自泰國。來自廣東的瓜形短棒形、瓜皮色為綠白的瓠瓜被單獨(dú)列為一類。

3? 討論

本研究基于瓠瓜轉(zhuǎn)錄組獲得87 518條uni gene序列進(jìn)行SSR位點(diǎn)搜尋，共檢測出11 029個(gè)SSR位點(diǎn)，總發(fā)生頻率為9.85%，分布距離為平均每8.3 kb有1個(gè)SSR位點(diǎn)，此結(jié)果略高于大王杜鵑（9.32%）[9]、花楸（6.86%）[10]、芒果（7.57%）[11]、枇杷（8.55%）[12]等的頻率范圍，但顯著低于蘿卜（20.37%）[13]、桃（20.76%）[14]和黃秋葵（20.59%）[6]的位點(diǎn)頻率，表明不同物種中轉(zhuǎn)錄組的SSR發(fā)生頻率并沒有特定的規(guī)律，在進(jìn)化過程中不同的經(jīng)歷都有可能影響SSR的位點(diǎn)發(fā)生頻率，同時(shí)受轉(zhuǎn)錄組SSR挖掘工具及篩選的條件影響[15]。但從總體上看瓠瓜轉(zhuǎn)錄組中SSR數(shù)量較豐富。

SSR重復(fù)類型也是轉(zhuǎn)錄組的重要特征之一。瓠瓜SSR重復(fù)類型，包含2核苷酸至6核苷酸重復(fù)5大類，其中3核苷酸比例最占優(yōu)勢，占總SSR的65.02%;其次是2核苷酸占總SSR的31.29%;4核苷酸、5核苷酸和6核苷酸重復(fù)類型數(shù)量較少，分別占總數(shù)的2.24%、0.75%和0.70%。其與杜鵑[16]、菜薹[17]、黃秋葵[6]等作物的研究結(jié)果不盡相同，猜測原因是由于2、4和5核苷酸重復(fù)主要與非翻譯區(qū)（UTR）相關(guān)。而3、6核苷酸重復(fù)主要發(fā)生在編碼序列（CDS）內(nèi)[18]，瓠瓜中3核苷酸重復(fù)比例最高，說明其轉(zhuǎn)錄組unigene中可能包含更多的CDS信息。進(jìn)一步對(duì)重復(fù)基序進(jìn)行分析，發(fā)現(xiàn)以3核苷酸重復(fù)類型AAG/CTT出現(xiàn)最多，占總SSR的22.76%;其次是2核苷酸重復(fù)類型AG/CT，占總SSR的21.26%。這與黃瓜[19]、蓖麻[20]、芝麻[21]、黃秋葵[6]等雙子葉植物現(xiàn)象相類似，高比例AAG/CTT和AG/CT基序的出現(xiàn)可能是雙子葉植物SSR的特征之一。

轉(zhuǎn)錄組測序的出現(xiàn)，使得EST數(shù)量迅速增加，許多植物中EST-SSR標(biāo)記已被開發(fā)應(yīng)用[22]，但瓠瓜EST-SSR標(biāo)記鮮少報(bào)道。本研究利用基于瓠瓜轉(zhuǎn)錄組數(shù)據(jù)開發(fā)的SSR引物，篩選出6對(duì)多態(tài)性引物將36份瓠瓜品種資源區(qū)分開來，多態(tài)性比率為68.6%。多態(tài)性信息含量（PIC）是衡量引物多態(tài)性信息含量水平的重要指標(biāo)（高：PIC>0.5、適中：0.25

根據(jù)6對(duì)引物擴(kuò)增的59個(gè)多態(tài)性位點(diǎn)，為36份供試瓠瓜品種資源進(jìn)行了聚類分析，其遺傳相似系數(shù)變幅為0.60～0.97。在遺傳相似系數(shù)為0.60處36份瓠瓜品種資源可分為2個(gè)類群。在遺傳相似系數(shù)為0.68處可分為6類：第1類包含有19份品種資源，來源地主要集中在福建和湖北，還有2份來自泰國;瓜皮色綠色的占89%，說明此類可以通過地域和果皮顏色不同來區(qū)分。第2類包含來源于湖北、廣東和泰國3份品種資源。第3類包含4份瓠瓜品種資源，瓜型均為長棒形，瓜皮色綠色，且種質(zhì)資源均來自福建。第4類包含2份瓠瓜品種資源，來自福建和泰國，瓜型均為長棒形，瓜皮色分別為亮綠和綠色。第5類包含7份瓠瓜品種資源，其中來自福建和廣東的各有3份。第6類為來自廣東的短棒形、綠白瓜皮色的瓠瓜品種資源。通過聚類分析的結(jié)果表明，利用篩選得到的6對(duì)SSR引物可將36份瓠瓜品種資源中的大部分區(qū)分開，36份瓠瓜品種資源的聚類結(jié)果與種源地、瓜型和瓜皮色有很大的相關(guān)性。而王莎[1]利用14個(gè)SSR標(biāo)記對(duì)44份瓠瓜品種的多樣性進(jìn)行分析，表明相同瓜型的品種有聚集在一起的趨勢，不同品種資源的來源地分類卻不明顯。雖然本研究中仍有部分瓠瓜品種資源未能區(qū)分，其可能原因與市場上主栽瓠瓜品種間親緣關(guān)系近、不同地域之間的瓠瓜相互引種及同種異名等有關(guān)，且在瓠瓜的瓜皮色判斷上因栽培方式、人為因素等原因也容易出現(xiàn)混淆。

本研究結(jié)果表明，利用瓠瓜轉(zhuǎn)錄組數(shù)據(jù)開發(fā)出的SSR標(biāo)記具有可行性，為瓠瓜品種資源遺傳圖譜構(gòu)建、基因定位、克隆等育種工作奠定了基礎(chǔ)。現(xiàn)有瓠瓜核心品種資源間親緣關(guān)系近，在今后的育種工作中應(yīng)加大國外品種資源或野生種質(zhì)的引進(jìn)，發(fā)掘更多的瓠瓜優(yōu)異資源提高遺傳多樣性。

參考文獻(xiàn)

王? 莎. 瓠瓜SSR和InDel分子標(biāo)記的開發(fā)及應(yīng)用[D]. 金華： 浙江師范大學(xué)， 2013.

Decker-Walters D S， Wilkins-Ellert M， Staub C J E. Discovery and Genetic Assessment of Wild Bottle Gourd [Lagenaria Siceraria （Mol.） Standley; Cucurbitaceae] from Zimbabwe[J]. Economic Botany， 2004， 58（4）： 501-508.

杜改改， 孫? 鵬， 索玉靜， 等. 基于柿雌雄花芽轉(zhuǎn)錄組測序的SSR和SNP多態(tài)性分析[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2017， 22（10）： 45-55.

韓明利， 崔? 娜， 于志海， 等. 番茄果實(shí)EST資源SSR信息分析[J]. 華北農(nóng)學(xué)報(bào)， 2011， 26（4）： 213-217.

胡建斌， 李建吾. 黃瓜EST-SSR位點(diǎn)信息[J]. 江西農(nóng)業(yè)學(xué)報(bào)， 2008， 20（11）： 58-60， 68.

李永平， 劉建汀， 陳敏氡， 等. 利用黃秋葵轉(zhuǎn)錄組信息挖掘SSR標(biāo)記及用于種質(zhì)分析[J]. 園藝學(xué)報(bào)， 2018， 45（3）： 579-590.

馬紅勃， 許旭明， 韋新宇， 等. 基于SSR標(biāo)記的福建省若干水稻品種DNA指紋圖譜構(gòu)建及遺傳多樣性分析[J]. 福建農(nóng)業(yè)學(xué)報(bào)， 2010， 25（1）： 33-38.

段艷鳳， 劉? 杰， 卞春松， 等. 中國88個(gè)馬鈴薯審定品種SSR指紋圖譜構(gòu)建與遺傳多樣性分析[J]. 作物學(xué)報(bào)， 2009， 35（8）： 1451-1457.

蘇應(yīng)雄. 云南轎子山大王杜鵑的種子生理及遺傳多樣性研究[D]. 昆明： 云南大學(xué)， 2017.

張? 敏， 周彩紅， 陳? 燾， 等. 榧樹轉(zhuǎn)錄組SSR信息分析及其分子標(biāo)記開發(fā)[J]. 果樹學(xué)報(bào)， 2017， 34（10）： 1258-1265.

羅? 純， 武紅霞， 姚全勝， 等. 芒果轉(zhuǎn)錄組中SSR位點(diǎn)信息分析與引物篩選[J]. 熱帶作物學(xué)報(bào)， 2015， 36（7）： 1261-1266.

鄭婷婷， 魏偉淋， 楊向暉， 等. 基于枇杷轉(zhuǎn)錄組序列的SSR分子標(biāo)記引物開發(fā)[J]. 亞熱帶植物科學(xué)， 2015， 44（4）： 274-278.

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

何? 平， 李林光， 王海波，等. 桃（Prunus persica， [L.] Batsch）轉(zhuǎn)錄組SSR信息分析及其分子標(biāo)記開發(fā)[J]. 分子植物育種， 2016， 14（11）： 3130-3135.

楊? 彬， 許薔薇， 牛明月， 等. 云錦杜鵑轉(zhuǎn)錄組SSR分析及其分子標(biāo)記開發(fā)[J]. 核農(nóng)學(xué)報(bào)， 2018， 32（12）： 2335-2345.

李太強(qiáng)， 劉雄芳， 萬友名， 等. 滇東南瀕危植物長梗杜鵑轉(zhuǎn)錄組微衛(wèi)星特征分析[J]. 林業(yè)科學(xué)研究， 2017， 30（4）： 533-541.

李榮華， 王直亮， 陳靜芳， 等. 菜薹轉(zhuǎn)錄組中SSR信息與可用性分析[J]. 園藝學(xué)報(bào)， 2016， 43（9）： 1816-1824.

Cavagnaro P F， Senalik D A， Yang L， et al. Genome-wide characterization of simple sequence repeats in cucumber （Cucumis sativus L.）[J]. BMC Genomics， 2010， 11（1）： 569-587.

Wu J， Cai C F， Cheng F Y， et al. Characterisation and development of EST-SSR markers in tree peony using transcriptome sequences[J]. Molecular Breeding， 2014， 34（4）： 1853-1866.

Qiu L J， Yang C， Tian B， et al. Exploiting EST databases for the development and characterization of EST-SSR markers in castor bean （Ricinus communis L.）[J]. BMC Plant Biology， 2010， 10（1）： 278-288.

Wei W L， Qi X Q， Wang L H， et al. Characterization of the sesame （Sesamum indicum L.） global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers[J]. BMC Genomics， 2011， 12（1）： 451-464.

You Y M， Liu D C， Liu H B， et al. Development and characterisation of EST-SSR markers by transcriptome sequencing in taro （Colocasia esculenta （L.） Schoot）[J]. Molecular Breeding， 2015， 35（6）： 134.

徐玉仙， 張微微， 莫海波， 等. 基于EST-SSR標(biāo)記的蓮屬種質(zhì)資源遺傳多樣性分析[J]. 植物分類與資源學(xué)報(bào)， 2015， 37（5）： 595-604.