欒非時，盧丙洋,，周慧文,，王學(xué)征,

（1.農(nóng)業(yè)部東北地區(qū)園藝作物生物學(xué)與種質(zhì)創(chuàng)制重點實驗室，哈爾濱　150030；2.東北農(nóng)業(yè)大學(xué)園藝學(xué)院，哈爾濱　150030）

甜瓜遺傳連鎖圖譜構(gòu)建及果實相關(guān)性狀QTL定位

欒非時12，盧丙洋1,2，周慧文1,2，王學(xué)征1,2

（1.農(nóng)業(yè)部東北地區(qū)園藝作物生物學(xué)與種質(zhì)創(chuàng)制重點實驗室，哈爾濱150030；2.東北農(nóng)業(yè)大學(xué)園藝學(xué)院，哈爾濱150030）

試驗以厚皮甜瓜品系M4-5為母本材料，薄皮甜瓜品系M1-15為父本材料，配制雜交組合獲得F1、F2代，研究甜瓜果實相關(guān)性狀遺傳規(guī)律。對M4-5和M1-15兩親本間存在SNP突變位點序列作CAPS檢測，共設(shè)計CAPS標(biāo)記523對，CAPS標(biāo)記223對，多態(tài)率43.8%。利用多態(tài)性引物標(biāo)記F2代群體，構(gòu)建甜瓜遺傳連鎖圖譜，該圖譜包含195個標(biāo)記，12個連鎖群，1 702.55 cM，標(biāo)記間平均距離8.78 cM。28個，貢獻(xiàn)率介于3.4974%～17.9684%。其中與果實形狀相關(guān)位點14個，與產(chǎn)量相關(guān)位點9個，與可溶性固形物含量相關(guān)位點5個。

甜瓜；果實性狀；遺傳連鎖圖譜；QTL

網(wǎng)絡(luò)出版時間2016-8-24 15:04:00[URL]http://www.cnki.net/kcms/detail/23.1391.S.20160824.1504.004.html

欒非時,盧丙洋,周慧文,等.甜瓜遺傳連鎖圖譜構(gòu)建及果實相關(guān)性狀QTL定位[J].東北農(nóng)業(yè)大學(xué)學(xué)報,2016,47(8):9-20.

Luan Feishi,Lu Bingyang,Zhou Huiwen,et al.Construction of genetic linkage map and QTL location for fruit-related traits in melon[J].Journal of Northeast Agricultural University,2016,47(8):9-20.(in Chinese with English abstract)

甜瓜（Cucumis melon L.）是被子植物亞門、雙子葉植物綱、葫蘆目、葫蘆科、屬一年生蔓性草本植物。染色體2n=2x=24。我國現(xiàn)為世界甜瓜種植面積最大、產(chǎn)量最高的國家。果實形狀、體積、果重、可溶性固形物含量等性狀是影響果實品質(zhì)重要因素[1]。目前，果實形狀、果實重量、可溶性固形物含量等是甜瓜果實性狀基因定位研究主要對象。Périn等用同一個圓形果實甜瓜品系分別與兩個長形果實甜瓜品系雜交構(gòu)建兩個重組自交系群體，定位和分析控制甜瓜果形基因，兩個群體中果形性狀主要受果實長度影響，與果實形狀和子房形狀相關(guān)的QTL位點存在大量共區(qū)間現(xiàn)象，說明在子房發(fā)育過程中活躍著大量控制果實形狀基因，可能存在基因多效現(xiàn)象[2]。Monforte等利用F2和DH群體，共檢測到20個與甜瓜早熟性、果形及可溶性固形物含量相關(guān)QTL位點。在9號連鎖群檢測到1個與甜瓜外果皮顏色相關(guān)位點，在1號連鎖群檢測到和綠色果肉顏色相關(guān)位點?？刂瞥壬忸伾倪z傳因素相對復(fù)雜，在2號和12號連鎖群分別找到和橙色果肉顏色相關(guān)位點，擴(kuò)大群體規(guī)模有利于闡明控制橙色果肉遺傳因素[3]。果實重量是影響產(chǎn)量重要因素。Zalapa等在探究甜瓜產(chǎn)量相關(guān)性狀時共檢測到包括單株果重、單果平均重在內(nèi)5個性狀的37個QTL位點，試驗結(jié)果表明甜瓜種質(zhì)資源中高度分枝類型的基因更具增產(chǎn)潛力[1]。Obando等以RILs群體為基礎(chǔ)檢測到134個與甜瓜果實品質(zhì)性狀相關(guān)QTL位點，其中包括與形態(tài)學(xué)和外觀性狀相關(guān)位點52個，與顏色性狀相關(guān)位點69個，余下23個位點與果實風(fēng)味性狀相關(guān)[4]。Diaz等參考國際葫蘆科基因組計劃數(shù)據(jù)，將8個甜瓜圖譜整合為一個圖譜，分析果實相關(guān)性狀發(fā)現(xiàn)果實形狀遺傳可能性明顯高于果實重量和可溶性固形物含量[5]。Diaz等以中等大小、長果形的印度地方甜瓜品種PI1124112和果形較大、橢圓形西班牙栽培甜瓜品種PS為親本材料雜交獲得F2代群體，利用基于遺傳連鎖分析QTL定位，共檢測到10個和果實形狀相關(guān)QTL位點，分別與果實長度、果實直徑和果形相關(guān)，3個性狀中均檢測到貢獻(xiàn)率較大的QTL位點[6]。Ramamurthy等研究甜瓜果實形狀和子房形狀過程中，發(fā)現(xiàn)兩者之間存在顯著正相關(guān)，與Périn等結(jié)果相似，說明甜瓜果實形狀在開花期之前已確定。Ramamurthy等利用基于遺傳連鎖分析的QTL定位時，共檢測到31個與甜瓜果實品質(zhì)和形態(tài)性狀相關(guān)的QTL位點，其中23個位點得到證實，8個新定位的QTL位點涉及果實形狀、單株總果重、可溶性固形物含量、果肉顏色、子房形狀等[7]。

本試驗以厚皮甜瓜品系M4-5為母本材料，薄皮甜瓜品系M1-15為父本材料，配制雜交組合獲得F1、F2代，研究甜瓜果實相關(guān)性狀遺傳規(guī)律，利用CAPS標(biāo)記構(gòu)建遺傳連鎖圖譜，通過分子標(biāo)記找到控制果實相關(guān)性狀位點，為進(jìn)一步開展甜瓜果實相關(guān)基因定位和克隆研究提供理論依據(jù)。

1　材料與方法

1.1田間試驗材料

以遺傳和表型差異極大作為親本材料選擇依據(jù)，并配置雜交組合用以表型數(shù)據(jù)調(diào)查和圖譜構(gòu)建。母本材料為厚皮甜瓜伊麗莎白（M4-5），中早熟品系，果實圓形（果形指數(shù)0.9），果皮厚，果實大（單果重0.96 kg）。父本材料為薄皮甜瓜M1-15，早熟品系，果實長橢圓形（果形指數(shù)1.5），果皮薄，果實?。▎喂?.24 kg）。兩親本材料來源于東北農(nóng)業(yè)大學(xué)園藝學(xué)院西甜瓜分子育種研究室。利用母本M4-5和父本M1-15配置雜交組合，獲得F1代，F(xiàn)1代自交獲得F2代。

1.2田間試驗設(shè)計

2015年4月～2015年8月，在東北農(nóng)業(yè)大學(xué)設(shè)施園藝工程中心1號大棚播種P1、P2、F1和F2代群體，采用完全隨機(jī)區(qū)組設(shè)計，3次重復(fù)，親本P1、P2及F1各種植30株，每小區(qū)種植10株。隨機(jī)種植F2300株（見圖1），株行距35 cm×50 cm。待果實成熟后測定果實相關(guān)性狀。

圖1　田間布局Fig.1　Distribution in plastic greenhouse

1.3果實性狀調(diào)查方法

將收獲的父、母本，F(xiàn)1和F2群體中成熟甜瓜果實分株標(biāo)號，每株收獲1個果實，每個果實每個性狀測定3次，取測量3次平均值作結(jié)果，其中父、母本，F(xiàn)1代各測定30株，F(xiàn)2代測定284株，具體測定方法如下：

單果重：電子天平測量；

外果皮和食用果皮硬度：硬度計測量；

種腔、果實長度、寬度：直尺測量；

種腔與果實比例：種腔與果實長度比值；

可溶性固形物：手持糖度儀測定。

1.4分子標(biāo)記試驗內(nèi)容與方法

1.4.1CAPS分子標(biāo)記開發(fā)

以兩親本材料基因組重測序數(shù)據(jù)為依據(jù)，已發(fā)布甜瓜基因組數(shù)據(jù)為參考，利用東北農(nóng)業(yè)大學(xué)園藝學(xué)院西甜瓜實驗室自編Perl語言腳本提取SNP位點前后約500 bp堿基序列，利用SNP2CAPS[8]軟件查找存在差異的酶切位點候選序列，根據(jù)序列中突變位點所處位置，選擇適合CAPS標(biāo)記開發(fā)的序列設(shè)計引物。

1.4.2DNA提取

待M4-5（♀）、M1-15（♂）、F1代及其自交獲得的F2代群體定植緩苗后，采集幼嫩葉片，分袋保存，分別標(biāo)號，以免混淆，用以提取DNA。

采用改良CTAB法提取植物基因組DNA[9]，應(yīng)用氯仿、異戊醇（24 ? 1）抽提3次，以保證DNA純度。提取親本及F1的DNA用于CAPS標(biāo)記多態(tài)性篩選，F(xiàn)2單株DNA用于遺傳連鎖分析。

1.4.3CAPS標(biāo)記篩選與驗證

將引物母液按體積比1 ? 49稀釋后，采用降落式PCR（Touchdown PCR）反應(yīng)程序進(jìn)行引物PCR擴(kuò)增，PCR反應(yīng)程序參照束永俊等方法[10]，步驟略有改動，PCR產(chǎn)物在1%瓊脂糖凝膠電泳上檢驗。

參照限制性內(nèi)切酶操作說明書對酶切位點加以驗證，利用Eco RⅠ、HindⅢ、PstⅠ、Bam HⅠ 4種限制性核酸內(nèi)切酶分別對相應(yīng)PCR產(chǎn)物酶切驗證，酶切體系為15.3 μL，酶切反應(yīng)在水浴鍋中進(jìn)行，具體酶切溫度根據(jù)4種不同限制性內(nèi)切酶最適反應(yīng)溫度設(shè)定，均為37℃。1%瓊脂糖凝膠電泳酶切產(chǎn)物檢驗。

1.5遺傳連鎖圖譜構(gòu)建

利用IciMapping V3.3對具有多態(tài)性且應(yīng)用到F2代基因分型的分子標(biāo)記構(gòu)建遺傳圖譜，采用復(fù)合區(qū)間作圖法構(gòu)建遺傳連鎖圖譜。利用MapChart2.2軟件繪制遺傳連鎖圖譜，依次將連鎖群命名為LG1～LG12。

1.6QTL定位分析方法

QTL分析使用IciMapping V3.3軟件，采用復(fù)合區(qū)間作圖法，以1.0 cM步行速度在全基因組內(nèi)掃描。以性狀英文縮寫、連鎖群編號和QTL編號為依據(jù)命名其QTL。

2　結(jié)果與分析

2.1甜瓜果實相關(guān)性狀表型數(shù)據(jù)分析

2.1.1親本及F1果實表型數(shù)據(jù)分析

由表1可知，兩親本M4-5和M1-15果實各性狀存在差異，母本M4-5單果重、果實長、果實寬、種腔直徑均明顯大于父本M1-15，僅果形指數(shù)小于父本M1-15，可溶性固形物和種腔與果實比這兩個性狀兩親本差異不顯著。在F1代果實中，單果重、果實寬、果形指數(shù)、種腔直徑及種腔與果實比均介于兩親本之間，而果實長度存在超親現(xiàn)象，明顯大于兩親本值，可溶性固形物含量明顯小于兩親本。

表1　親本、F1及F2群體果實相關(guān)性狀表現(xiàn)Table 1　Fruit-related trait separation proportion of the parental materials,F1and F2generation

2.1.2果實單果重表型數(shù)據(jù)分析

由圖2可知，在單果重方面，母本M4-5單果重為（0.96±0.55）kg，父本M1-15單果重為（0.24± 0.13）kg，F(xiàn)1單果重為（0.60±0.68）kg，介于兩親本之間。在F2代果實中，單果重呈連續(xù)分離狀態(tài)，且存在超親現(xiàn)象，單果重分離范圍0.17～1.97 kg，整體分離趨勢基本符合正態(tài)分布，表明單果重受多個基因共同控制，同時存在主效基因，為數(shù)量性狀，可作QTL定位分析。

2.1.3果實長度表型數(shù)據(jù)分析

由圖3可知，在果實長度方面，母本M4-5果實長度為（11.6±0.21）cm，父本M1-15果實長度為（10.11±0.18）cm，F(xiàn)1果實長度為（13.36±0.57）cm，果實長度均高于兩親本。在F2代果實中，果實長度分離范圍7.00～20.80 cm，最高值和最低值均出現(xiàn)超親遺傳，分離趨勢基本符合正態(tài)分布，表明果實長度受多個基因共同控制，同時存在主效基因，為數(shù)量性狀，可作QTL定位分析。

圖2　果實單果重次數(shù)分布Fig.2　Fruit weight proportion distribution in F2population

圖3　果實長度次數(shù)分布Fig.3　Fruit length proportion distribution in F2population

2.1.4果實寬度表型數(shù)據(jù)分析

由圖4可知，在果實寬度方面，母本M4-5果實寬度為（12.95±0.29）cm，父本M1-15果實寬度為（6.64±0.12）cm，F(xiàn)1果實寬度為（9.59±0.44）cm，F(xiàn)1代表型數(shù)據(jù)介于父本和母本之間。在F2代果實中，果實寬度為6.50～16.00 cm，與兩親本相比存在超親現(xiàn)象，分離趨勢呈連續(xù)正態(tài)分布，表明果實寬度受多個基因共同控制，同時存在主效基因，為數(shù)量性狀，可作QTL定位分析。

2.1.5果形指數(shù)表型數(shù)據(jù)分析

由圖5可知，在果形指數(shù)方面，母本M4-5果形指數(shù)為0.90±0.13，果實為近圓形，父本M1-15果形指數(shù)為1.53±0.18，果實為長橢圓形，F(xiàn)1代果形指數(shù)為1.40±0.25，果形指數(shù)介于兩親本之間，且偏向父本M1-15，為長橢圓形。

在F2代果實中，果形指數(shù)分離范圍0.89～ 1.98，既有親本類型的近圓形和長橢圓形，也存在兩者之間過渡類型且有超親現(xiàn)象。分離情況呈連續(xù)正態(tài)分布，表明果形指數(shù)為數(shù)量性狀，可作QTL定位分析。

圖4　果實寬度次數(shù)分布Fig.4　Fruit width proportion distribution in F2population

圖5　果形指數(shù)次數(shù)分布Fig.5　Fruit shape index proportion distribution in F2population

2.1.6種腔表型數(shù)據(jù)分析

由圖6可知，母本M4-5種腔直徑為（6.55± 0.17）cm，父本M1-15種腔直徑為（3.78±0.69）cm，父母本差異顯著。F1代種腔直徑為（5.16±0.24）cm，介于兩親本之間，且偏向于父本。F2代種腔直徑3.40～9.50 cm，存在超親現(xiàn)象，分離情況呈正態(tài)分布，為數(shù)量性狀，可作QTL定位分析。

2.1.7種腔與果實比表型數(shù)據(jù)分析

由圖7可知，種腔與果實比可顯示果實可食用率，通過種腔與果實比數(shù)據(jù)顯示，M4-5和M1-15可食用率基本相同，F(xiàn)1代數(shù)據(jù)介于兩親本之間，F(xiàn)2代群體中可食用率出現(xiàn)較大變異，種腔與果實最大比值為0.76，最小比值為0.40，呈正態(tài)分布，為數(shù)量性狀。

2.1.8可溶性固形物表型數(shù)據(jù)分析

由圖8可知，母本M4-5可溶性固形物含量為（9.24±0.43）%，父本M-15可溶性固形物含量為（10.69±0.26）%，F(xiàn)1代果實可溶性固形物含量低于兩親本，F(xiàn)2代群體可溶性固形物含量4.00%～14.00%，存在較大幅度變異，均出現(xiàn)超親遺傳現(xiàn)象。F2代群體呈正態(tài)單峰分布，表明可溶性固形物含量受多個基因共同控制，同時存在主效基因，為數(shù)量性狀。

圖6　種腔次數(shù)分布Fig.6　Seed cavity proportion distribution in F2population

圖7　種腔與果實比次數(shù)分布Fig.7　Ratio of seed cavity and fruit proportion distribution in F2population

圖8　可溶性固形物次數(shù)分布Fig.8　Brix content proportion distribution in F2population

2.2甜瓜遺傳連鎖圖譜構(gòu)建

2.2.1甜瓜基因組DNA提取

利用改良CTAB法提取基因組DNA，經(jīng)1%瓊脂糖凝膠電泳檢測，電泳結(jié)果見圖9，父母本、F1及F2代分離群體DNA濃度和純度均可用于分子標(biāo)記檢測，無RNA和蛋白質(zhì)，對分子標(biāo)記檢測不存在影響。

圖9　甜瓜基因組DNAFig.9　Genomic DNA of melon

2.2.2CAPS引物篩選與驗證

通過測序結(jié)果，根據(jù)HindⅢ、Eco RⅠ、Bam HⅠ和PstⅠ酶切位點信息對M4-5和M1-15兩親本間有存在SNP突變位點序列作CAPS檢測，利用Primer 6設(shè)計并選用523對引物。利用提取的親本及F1代DNA，對設(shè)計的523對CAPS引物作PCR和酶切反應(yīng)，驗證引物多態(tài)性和共顯性。圖10為部分有多態(tài)性引物經(jīng)PCR和酶切反應(yīng)后產(chǎn)物電泳條帶。驗證結(jié)果顯示，523對CAPS引物中有229對存在多態(tài)性和共顯性，多態(tài)率達(dá)43.8%。

圖10　CAPS引物在親本及F1代多態(tài)性表現(xiàn)Fig.10　Polymorphism of CAPS markers among parental materials and F1generation

2.2.3甜瓜遺傳圖譜構(gòu)建

利用篩選出的多態(tài)性引物在F2代群體中作分析，IciMapping V3.3軟件對多態(tài)性CAPS標(biāo)記連鎖作圖，排除條帶不清與偏分離標(biāo)記，用于作圖引物共計195對，如圖11所示，構(gòu)建一張包含12個連鎖群，覆蓋基因組長度1 702.55 cM，標(biāo)記間平均距離8.78 cM的甜瓜F2代臨時群體遺傳圖譜，圖譜中來源于同一條染色體CAPS標(biāo)記均位于相同連鎖群上，將連鎖群依次命名為LG1～LG12，連鎖群長度介于105.16～187.34 cM。該圖譜最長連鎖群為LG2，長度為187.34 cM，包含18個CAPS標(biāo)記，最短連鎖群為LG10，長度105.16 cM，包含12個CAPS標(biāo)記（見表2和圖11），包含最多標(biāo)記連鎖群為LG4，共23個標(biāo)記，包含標(biāo)記最少連鎖群為LG10，共12個標(biāo)記。連鎖群LG10為包含標(biāo)記最少，長度最短連鎖群。

2.3果實相關(guān)性狀QTL定位分析

2.3.1單果重QTL分析

對于單果重，如表3所示，共定位到3個QTL，均為微效QTL位點，分別位于連鎖群LG5、LG6和LG11上。其中FW5.1和FW6.1對單果重起增加作用，而FW11.1起減小作用。這3個QTL位點中貢獻(xiàn)率最大的是位于連鎖群LG5上的FW5.1，貢獻(xiàn)率為6.8534%，定位在標(biāo)記M5-37和M5-22之間，距離M5-37距離1.02 cM，距離M5-22距離0.33 cM，距離兩標(biāo)記距離均相對較近。貢獻(xiàn)率最小的為定位在連鎖群LG11的FW11.1，貢獻(xiàn)率僅3.7775%，位于標(biāo)記M11-2和M11-22之間，距離兩標(biāo)記距離分別為9.65和0.28 cM，該位點只與一端標(biāo)記距離近，而與另一端標(biāo)記距離則較遠(yuǎn)。

2.3.2果實長度QTL分析

通過田間數(shù)據(jù)與分子數(shù)據(jù)結(jié)合分析，共檢測到與果實長度相關(guān)的QTL位點5個（見表4），其中包含1個貢獻(xiàn)率>10%的QTL位點，其余QTL位點貢獻(xiàn)率3.4974%～4.6475%。這5個QTL分別位于連鎖群LG3、LG5、LG6、LG9和LG11上。除FL5.1對增加甜瓜果實長度表現(xiàn)為正的加性效應(yīng)外，其余4個QTL位點對增加甜瓜果實長度表現(xiàn)為負(fù)的加性效應(yīng)，即降低果實長度。定位在連鎖群LG3 的FL3.1，定位在連鎖群LG6的FL6.1和定位在連鎖群LG11的FL11.1這3個QTL位點均與一側(cè)標(biāo)記距離較近，分別為0.04、0.39和0.28 cM，而距另一側(cè)標(biāo)記相對較遠(yuǎn)，分別為13、10.93和9.65 cM。貢獻(xiàn)率最大的FL5.1定位距離與其他QTL位點相比，與相鄰標(biāo)記間距離相對較遠(yuǎn)，分別為14和10.01 cM。

圖11　甜瓜遺傳連鎖圖譜及果皮與果實相關(guān)性狀定位位點Fig.11　Genetic linkage map and locus for pericarp-related and fruit-related traits of melon

表2　甜瓜遺傳連鎖圖譜基本參數(shù)Table 2　Parameter of melon genetic linkage map

表3　甜瓜單果重QTL遺傳效應(yīng)分析Table 3　Genetic effects analysis of QTLs for fruit weight in melon

2.3.3果實寬度QTL分析

如表5所示，關(guān)于果實寬度僅定位到兩個QTL位點，分別位于連鎖群LG5和LG6上，兩標(biāo)記貢獻(xiàn)率分別達(dá)7.1447%和9.4046%，同時對于果實寬度均起增加作用。在連鎖群LG5上定位到的FWID5.1距離標(biāo)記M5-48距離為3.01 cM，距離標(biāo)記M5-49距離為0.64 cM。在連鎖群LG6上定位到的FWID6.1距離兩標(biāo)記距離分別為2.99和3.07 cM。與定位到的關(guān)于果實長度QTL位點相比，關(guān)于果實寬度的QTL位點定位距離相對較近。

2.3.4果形指數(shù)QTL分析

定位到與果形指數(shù)相關(guān)的QTL位點共計7個（見表6），分別位于連鎖群LG1、LG2、LG4、LG7、LG9和LG10。其中貢獻(xiàn)率為17.9684%的主效QTL，位于連鎖群LG7上，該位點被定位在標(biāo)記M7-39和M7-24之間，距離兩標(biāo)記分別為1.06 和1.80 cM。在連鎖群LG4上，共檢測到兩個與果形指數(shù)有關(guān)QTL位點，兩位點貢獻(xiàn)率為7個QTL中最小。這7個關(guān)于果形指數(shù)的QTL貢獻(xiàn)率累計達(dá)41.55%。在LG2上定位到的FSI2.1與標(biāo)記M2-12緊密連鎖，距離僅為0.16 cM。在定位到與果形指數(shù)相關(guān)的QTL位點中，定位結(jié)果與果實長度QTL位點一樣，也存在QTL位點只與一側(cè)標(biāo)記緊密連鎖，與另一標(biāo)記距離相對較遠(yuǎn)的情況。

表5　 甜瓜果實寬度QTL遺傳效應(yīng)分析Table 5　Genetic effects analysis of QTLs for fruit width in melon

表6　甜瓜果形指數(shù)QTL遺傳效應(yīng)分析Table 6　Genetic effects analysis of QTLs for fruit shape index in melon

2.3.5種腔QTL分析

種腔體積既影響果實食用率同時也影響產(chǎn)量，種腔越小，食用率越高，產(chǎn)量也越高。在種腔定位中，共定位到3個QTL位點，分別位于連鎖群LG4、LG5和LG6（見表7），其中定位在連鎖群LG5上的QTL位點SC5.1與相鄰兩標(biāo)記距離最近，分別為0.67和0.95 cM。定位在LG6上的QTL位點SC6.1貢獻(xiàn)率最大，達(dá)7.8495%。但并未定位到關(guān)于種腔的主效QTL位點。同時發(fā)現(xiàn)種腔和果實寬度共同定位在標(biāo)記M6-39和M6-22之間。

2.3.6果實可溶性固形物含量QTL分析

如表8所示，共定位到關(guān)于可溶性固形物含量QTL位點5個，連鎖群LG2上存在1個，LG6上存在1個，LG7上存在2個，LG9上存在1個。5個QTL貢獻(xiàn)率達(dá)3.8933%-5.7898%，BRC6.1和BRC9.1加性效應(yīng)值為負(fù)值，表明對果實可溶性固形物含量起減小作用，其余QTL加性效應(yīng)值為正值，對果實可溶性固形物含量起增加作用。

2.3.7種腔與果實比QTL分析

種腔與果實比可從側(cè)面反映果皮厚度及可食用率，比值越小可食用率越大。結(jié)合表型數(shù)據(jù)與分子數(shù)據(jù)，共計定位到3個QTL位點，分別位于連鎖群LG4、LG8和LG9上。其中CFR4.1加性效應(yīng)為正，表明增加種腔與果實比，即對可食用率起降低作用，而其余位點加性效應(yīng)為負(fù)，表明對可食用率起增加作用。

在定位距離上，CFR4.1和CFR8.1出現(xiàn)與一側(cè)標(biāo)記緊密連鎖的情況，與另一側(cè)標(biāo)記距離相對較遠(yuǎn)。CFR9.1與兩標(biāo)記間距離分別為4.87和2.08 cM（見表9）。

表7　 甜瓜種腔QTL遺傳效應(yīng)分析Table 7　Genetic effects analysis of QTLs for seed cavity in melon

表8　 甜瓜可溶性固形物QTL遺傳效應(yīng)分析Table 8　Genetic effects analysis of QTLs for brix content in melon

表9　甜瓜果實與種腔比QTL遺傳效應(yīng)分析Table 9　Genetic effects analysis of QTLs for seed cavity and fruit ratio in melon

3　討　論

3.1基于重測序的分子標(biāo)記開發(fā)

本研究用于構(gòu)建遺傳連鎖圖譜的195個CAPS分子標(biāo)記均源于甜瓜基因組重測序數(shù)據(jù)開發(fā)，是該圖譜最顯著的特點。

與其他標(biāo)記相比，CAPS標(biāo)記具有以下優(yōu)勢：①兩親本材料序列比對后設(shè)計的CAPS引物具有更強(qiáng)針對性和更高多態(tài)性，本試驗中CAPS引物多態(tài)性達(dá)43.8%。②用瓊脂糖凝膠電泳即能對CAPS標(biāo)記檢測，方便且易于觀測。③基因組中數(shù)量龐大的單堿基突變?yōu)闃?gòu)建高飽和遺傳連鎖圖譜提供更大可能性。④多數(shù)單堿基突變位點位于基因編碼區(qū)及功能區(qū)，利于對相關(guān)性狀精細(xì)定位。⑤所構(gòu)建遺傳連鎖圖譜與實際染色體相對應(yīng)，基因定位更有實際意義。

3.2果實相關(guān)性狀分析及QTL定位

通過相關(guān)性分析結(jié)果顯示，甜瓜果實性狀大部分均顯著相關(guān)，在定位分析中，這些性狀也均在相同連鎖群上，有的甚至定位在相同區(qū)域。這種現(xiàn)象存在于甜瓜中[7]，也存在于其他作物中，如西瓜、茄子、辣椒、番茄、玉米、紅麻等[11-19]。Ramamurthy等研究表明，與甜瓜果實長度、果實寬度、果形指數(shù)、單果重、果肉厚度等相關(guān)QTL位點富集于連鎖群LG8上[7]，本研究結(jié)果也出現(xiàn)QTL簇，QTL聚集在連鎖群LG6和LG9上，在LG6上同時定位到關(guān)于果實形狀、產(chǎn)量及食用果皮厚度和裂果等相關(guān)性狀位點，在LG9上定位到果實形狀、產(chǎn)量以及外果皮類型等相關(guān)性狀位點。對于這種現(xiàn)象，有學(xué)者認(rèn)為與基因多因一效和一因多效遺傳機(jī)制有關(guān)[20-21]，為聚合多種優(yōu)良性狀分子育種手段奠定理論基礎(chǔ)。

4　結(jié)論

4.1甜瓜遺傳圖譜構(gòu)建

利用生物信息學(xué)分析軟件挖掘SNP位點，對M4-5和M1-15兩親本間存在SNP突變位點序列作CAPS檢測，共設(shè)計CAPS標(biāo)記523對，篩選出在親本間有多態(tài)性CAPS標(biāo)記223對，多態(tài)率43.8%。利用多態(tài)性引物對F2代群體作標(biāo)記分析，構(gòu)建甜瓜遺傳連鎖圖譜，該圖譜包含195個標(biāo)記，12個連鎖群，覆蓋基因組長度1 702.55 cM，標(biāo)記間平均距離8.78 cM。該圖譜最長連鎖群為LG2，最短連鎖群為LG10。

4.2甜瓜果實相關(guān)性狀遺傳分析

果實相關(guān)性狀包括果實形狀、產(chǎn)量、可溶性固形物，研究表明這些與果實相關(guān)性狀均為數(shù)量性狀。通過遺傳分析和相關(guān)性研究發(fā)現(xiàn)果皮及果實相關(guān)性狀密切相關(guān)。結(jié)合表型數(shù)據(jù)與分子數(shù)據(jù)，共計定位到果皮及果實相關(guān)性狀位點28個，并明確其在連鎖群上分布情況及位點遺傳效應(yīng)。在果實相關(guān)性狀定位中，與形狀相關(guān)位點共計14個，與產(chǎn)量相關(guān)位點共計9個，與可溶性固形物含量相關(guān)位點5個。

[1]Zalapa J E,Staub J E,Mccreight J D,et al.Detection of QTL for yield-related traits using recombinant inbred lines derived from exotic and elite US western shipping melon germplasm[J].Theo?retical&Applied Genetics,2007,114(7):1185-1201.

[2]Périn C,Hagen L S,Giovinazzo N,et al.Genetic control of fruit shape acts prior to anthesis in melon(Cucumis melo L.)[J].Molec?ular Genetics&Genomics Mgg,2002,266(6):933-941.

[3]Monforte A J,Oliver M,Gonzalo M J,et al.Identification of quan?titative trait loci involved in fruit quality traits in melon(Cucumis melo L.)[J].Tag Theoretical&Applied Genetics Theoretische Und Angewandte Genetik,2004,108(4):750-758.

[4]Obando J,Fernández-Trujillo J P,Martinez J A,et al.Identifica?tion of melon fruit quality quantitative trait loci using near-isogen?ic lines[J].Journal of the American Society for Horticultural Sci?ence,2008,133(1):139-151.

[5]Díaz A,Fergany M,Formisano G,et al.A consensus linkage map for molecular markers and quantitative trait loci associated with economically important traits in melon(Cucumis melo L.)[J]. BMC Plant Biology,2011,11(2):121-136.

[6]Díaz A,Zarouri B,Fergany M,et al.Mapping and introgression of QTL involved in fruit shape transgressive segregation into‘piel de sapo,melon(Cucucumis melo L.)[J].PloS One,2014,9(8): e104188.

[7]Ramamurthy R K,Waters B M.Identification of fruit quality and morphology QTLs in melon(Cucumis melo)using a population de?rived from flexuosus and cantalupensis botanical groups[J].Eu?phytica,2015,204(1):1-15.

[8]Thiel T,Kota R,Grosse I,et al.SNP2CAPS:A SNP and INDEL analysis tool for CAPS marker development[J].Nucleic Acids Re?search,2004,32(1):5.

[9]Luan F,Delannay I,Staub J E.Chinese melon(Cucumis melo L.) diversity analyses provide strategies for germplasm curation,ge?netic improvement,and evidentiary support of domestication pat? terns[J].Euphytica,2008,164(2):445-461.

[10]束永俊,李勇,吳娜拉胡,等.大豆EST-SNP的挖掘、鑒定及其CAPS標(biāo)記的開發(fā)[J].作物學(xué)報,2010,36(4):574-579.

[11]Chaim A B,Paran I,Grube R C,et al.QTL mapping of fruit-relat?ed traits in pepper(Capsicum annuum)[J].Theoretical&Applied Genetics,2001,102(6-7):1016-1028.

[12]程瑤,劉識,朱子成,等.西瓜糖含量相關(guān)性狀的QTL分析[J].中國蔬菜,2016(2):24-31.

[13]喬軍,陳鈺輝,王利英,等.茄子果形的QTL定位[J].園藝學(xué)報, 2012,39(6):1115-1122.

[14]Saliba-Colombani V,Causse M,Langlois D,et al.Genetic analy?sis of organoleptic quality in fresh market tomato.1.Mapping QTLs for physical and chemical traits[J].Theoretical&Applied Genetics,2001,102(2-3):259-272.

[15]Prothro J,Sandlin K,Abdel-Haleem H,et al.Main and epistatic quantitative trait loci associated with seed size in watermelon[J]. Journal of the American Society for Horticultural Science Ameri?can Society for Horticultural Science,2012,137(6):452-457.

[16]Katherine S,Jason P,Adam H,et al.Comparative mapping in wa?termelon[Citrullus lanatus(Thunb.)Matsum.et Nakai][J].Tag Theoretical&Applied Genetics Theoretische Und Angewandte Genetik,2012,125(8):1603-1618.

[17]Ren Y,Mcgregor C,Zhang Y,et al.An integrated genetic map based on four mapping populations and quantitative trait loci asso?ciated with economically important traits in watermelon(Citrullus lanatus)[J].BMC Plant Biology,2014,14(1):1.

[18] 趙璞,劉瑞響,李成璞,等.基于掖478導(dǎo)入系的玉米產(chǎn)量性狀QTL鑒定[J].中國農(nóng)業(yè)科學(xué),2011,44(17):3508-3519.

[19] 陳美霞,祁建民,方平平,等.紅麻6個重要產(chǎn)量性狀的QTL定位[J].中國農(nóng)業(yè)科學(xué),2011,44(5):874-883.

[20]Katherine S,Jason P,Adam H,et al.Comparative mapping in wa?termelon[Citrullus lanatus(Thunb.)Matsum.et Nakai][J].Tag Theoretical&Applied Genetics Theoretische Und Angewandte Genetik,2012,125(8):1603-1618.

[21]Zhang X,Rhodes B,Wang M,et al.Genes controlling watermelon seed size[J].Cucurbitaceae,1994,94:144-147.

Construction of genetic linkage map and QTL location for fruit-related

traits in melon/LUAN Feishi1,2,LU Bingyang1,2,ZHOU Huiwen1,2,WANG Xuezheng1,2(1.Ministry of Agriculture KeyLaboratory of Biology and Germplasm Enhancement of Horticultural Crops in Northeast China,Harbin 150030,China;2.School of Horticulture,Northeast Agricultural University, Harbin 150030,China)

The present study was to develop SNPs based on the resequencing data and bioinformatics methods of M4-5 and M1-15 and convert them into cleaved amplified polymorphic sequence (CAPS)markers.In this study,we developed 523 CAPS markers totally and 223 markers showed polymorphism between parents,the polymorphism rate was 43.8%.We markered F2generation using CAPS markers and constructed a genetic map.The map contained 195 markers and grouped on 12 linkage groups,spanning a total length of 1 702.55 cM,with an average of 8.78 cM between markers.A total of 28 locus for friut-related traits which contained 14 locus for fruit,nine locus for yield and five locus brix content with contribution rate ranged from 3.4974%to 17.9684%.

melon;fruit trait;genetic linkage;QTL

S652

A

1005-9369（2016）08-0009-12

2016-04-11

國家西甜瓜產(chǎn)業(yè)技術(shù)體系（CARS-026-02)

欒非時（1964-），女，教授，博士，研究方向為西瓜甜瓜遺傳育種。E-mail:luanfeishi@sina.com