本發(fā)明屬于分子生物學(xué)領(lǐng)域，具體涉及大瀧六線魚單核苷酸多態(tài)性位點(diǎn)。
背景技術(shù)：
：大瀧六線魚(學(xué)名：hexagrammosotakii)為六線魚科六線魚屬的魚類。分布于朝鮮、日本以及中國東海、黃海、渤海等海域，系冷溫性近海底層魚類。大瀧六線魚的食性很廣，以底棲生物為主，也攝食小蝦、小魚。大瀧六線魚素以其味道鮮美肉質(zhì)細(xì)膩而聞名，素有“北方石斑”之稱，與其它魚類相比大瀧六線魚的含水量較低(73.85％)，可食部分比值高，蛋白質(zhì)含量極高(18.50％)，必需氨基酸含量高(7.25％)(kangbin)。大瀧六線魚是一種十分具有發(fā)展?jié)摿Φ慕?jīng)濟(jì)魚種，但大瀧六線魚生長發(fā)育十分緩慢，體長一般約為15-25厘。大瀧六線魚性成熟較早，性成熟時(shí)間一般在2～3齡，雄魚早于雌魚，性成熟后體重增長放緩，體重拐點(diǎn)在3.6齡。目前，國內(nèi)外有關(guān)大瀧六線魚的研究報(bào)道還比較少，主要集中在生物學(xué)地位、生態(tài)分布以及種群資源量，還有一些學(xué)者對其生理、生化等方面進(jìn)行了研究。對于大瀧六線魚遺傳學(xué)方面的研究相對較少，目前國外只有habib等人基于coi、coiii-nd3-nd4l和cytb對黃海和日本海附近群體遺傳結(jié)構(gòu)進(jìn)行研究，國內(nèi)的研究有劉奇和李瑩等人關(guān)于大瀧六線魚遺傳多樣性的研究，但取樣地都集中在黃渤海海域，且只使用了d-loop序列進(jìn)行遺傳多樣性與遺傳結(jié)構(gòu)的分析，并未使用cytb或coi基因進(jìn)行對比分析。未見大瀧六線魚單核苷酸多態(tài)性方面的研究。單核苷酸多態(tài)性標(biāo)記(snp：simplenucleotidepolymorphism)，是一種廣泛應(yīng)用的遺傳學(xué)領(lǐng)域的dna分子標(biāo)記。snp是廣泛存在真核生物基因組中的單核苷酸突變，一般有轉(zhuǎn)換、顛換、缺失、插入四種形式，其中動(dòng)物體中以轉(zhuǎn)換顛換最為常見。單核苷酸多態(tài)性具有多態(tài)性高，提供的遺傳信息多，pcr擴(kuò)增效果重復(fù)性好，以及在基因組中分散分布等優(yōu)點(diǎn)，常用于生物的遺傳多樣性分析，遺傳圖譜和雜交育種等。由于單核苷酸多態(tài)性研究所用dna的數(shù)量少，對樣品要求低，因此，單核苷酸多態(tài)性分析在遺傳學(xué)領(lǐng)域具有廣泛的應(yīng)用性。作為中國重要的漁業(yè)資源之一，開發(fā)大瀧六線魚與生長發(fā)育經(jīng)濟(jì)性狀相關(guān)的snp標(biāo)記研究是十分重要的。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的在于提供與大瀧六線魚生長發(fā)育相關(guān)的單核苷酸多態(tài)性位點(diǎn)，本發(fā)明提供的單核苷酸多態(tài)性位點(diǎn)選自：seqnoid：1所示序列第177位等位基因?yàn)閍/t，seqnoid：2所示序列第228位等位基因?yàn)間/a，seqnoid：3所示序列第228位等位基因?yàn)閏/t，seqnoid：4所示序列第269位等位基因?yàn)閍/g。本發(fā)明提供的單核苷酸多態(tài)性位點(diǎn)，其中與大瀧六線魚體長增長關(guān)聯(lián)顯著的單核苷酸多態(tài)性位點(diǎn)選自：seqnoid：1所示序列第177位等位基因?yàn)閍/t，seqnoid：3所示序列第228位等位基因?yàn)閏/tseqnoid：4所示序列第269位等位基因?yàn)閍/g。本發(fā)明提供的單核苷酸多態(tài)性位點(diǎn)，其中與大瀧六線魚體重增長關(guān)聯(lián)顯著的單核苷酸多態(tài)性位點(diǎn)選自：seqnoid：2所示序列第228位等位基因?yàn)間/a，seqnoid：3所示序列第228位等位基因?yàn)閏/t。本發(fā)明還提供所述單核苷酸多態(tài)性位點(diǎn)用于評估大瀧六線魚生長發(fā)育、大瀧六線魚選種或大瀧六線魚育種的用途。本發(fā)明提供一種大瀧六線魚單核苷酸多態(tài)性的檢測方法，包括以下步驟：1)基因組dna的提?。翰捎锰旄Ｑ髣?dòng)物組織dna提取試劑盒2)單核苷酸多態(tài)性pcr擴(kuò)增：利用assaydesign3.1軟件設(shè)計(jì)單核苷酸多態(tài)性引物擴(kuò)增大瀧六線魚基因組dna，獲得其擴(kuò)增產(chǎn)物。3)擴(kuò)增產(chǎn)物檢測：啟動(dòng)massarraynanodispenserrs1000點(diǎn)樣儀，將純化后的產(chǎn)物進(jìn)行點(diǎn)樣，并用maldi-tof質(zhì)譜儀分析，檢測結(jié)果使用typer4.0軟件進(jìn)行分析。4)遺傳多樣性分析：根據(jù)每個(gè)個(gè)體堿基類型取定基因型，利用popgene32計(jì)算遺傳多樣性參數(shù)。本發(fā)明體用提供能擴(kuò)增出所述單核苷酸多態(tài)性位點(diǎn)的的引物對，引物對具有選自以下序列對所示的序列：seqidno：5和seqidno：6；seqidno：7和seqidno：8；seqidno：9和seqidno：10；seqidno：11和seqidno：12。本發(fā)明提供所述引物對在大瀧六線魚單核苷酸多態(tài)性位點(diǎn)檢測中的用途。本發(fā)明提供了所述引物對在評估大瀧六線魚生長發(fā)育、大瀧六線魚選種或大瀧六線魚育種的用途。實(shí)施例一、單核苷酸多態(tài)性位點(diǎn)篩選首先采用2b-rad技術(shù)與高通量測序結(jié)合，對大瀧六線魚進(jìn)行簡化基因組測序，共發(fā)現(xiàn)158,733個(gè)snp，分布在27,404個(gè)序列標(biāo)簽中。通過blast比對篩選其中與生長發(fā)育相關(guān)且位于編碼區(qū)域的snp位點(diǎn)，共得到21個(gè)snp位點(diǎn)。通過sequenom實(shí)驗(yàn)平臺采用iplexgold技術(shù)及massarray系統(tǒng)檢測了21個(gè)snp位點(diǎn)在四個(gè)地理群體中的多樣性。四個(gè)地理群體所使用的樣本分別采自大連(n39°25′，e122°51′，野生群體，29尾)，舟山(n30°03′，e122°37′，野生群體，20尾)，瑯琊臺(n35°30′，e119°58′，野生群體，30尾)，即墨(即墨山東省海洋生物研究院養(yǎng)殖場，養(yǎng)殖群體，60尾)，取尾鰭和肌肉組織，存放于95％的無水乙醇中，放于-80℃冰箱中備用?；蚪Mdna的提?。菏褂帽本┨旄萍加邢薰镜暮Ｑ髣?dòng)物組織基因組dna提取試劑盒提取大瀧六線魚樣品基因組dna。具體流程如下：首先采用assaydesigner3.1軟件進(jìn)行引物設(shè)計(jì)，通過在四個(gè)地理群體120個(gè)樣本中進(jìn)行pcr擴(kuò)增，而后擴(kuò)增產(chǎn)物去鹽后4000rpm離心4分鐘，在芯片上點(diǎn)樣，再點(diǎn)calibrant，并將芯片放到scout板上并放回compact中，使用massarray進(jìn)行分析、質(zhì)量控制及報(bào)告輸出。報(bào)告結(jié)果采用popgene1.32計(jì)算等位基因數(shù)(observednumberofalleles，na)、有效等位基因數(shù)(effectivenumberofalleles，ne)、觀望雜合度(observedheterozygosity，ho)、期望雜合度(expectedheterozygosity，he)、nei’氏多樣性指數(shù)(nei′sgenediversityindex，nei)。多態(tài)信息含量(polymorphisminformationcontent，pic)計(jì)算公式為(式中pi和pj是第i和第j個(gè)等位基因頻率，n是等位基因數(shù)目)。分析結(jié)果顯示17個(gè)位點(diǎn)具有多態(tài)性，平均等位基因數(shù)為1.9524，平均有效等位基因數(shù)是1.4571，平均觀望雜合度0.2158，平均期望雜合度為0.2432，平均多態(tài)信息含量為0.2517。二、snp位點(diǎn)與生長發(fā)育性狀關(guān)聯(lián)性研究在關(guān)聯(lián)性研究中將多態(tài)性位點(diǎn)與大瀧六線魚體重體長增量的生長發(fā)育性狀進(jìn)行關(guān)聯(lián)。在即墨樣品中選取30尾魚體重體長相近，年齡約在1齡左右的大瀧六線魚進(jìn)行熒光魚體t型外掛標(biāo)簽標(biāo)記，分別記錄初體重與初體長數(shù)據(jù)，在山東省海洋生物研究院即墨鰲山衛(wèi)基地進(jìn)行為期125天的養(yǎng)殖，期間每天喂食兩次與養(yǎng)殖池內(nèi)的海水保持循環(huán)，記錄養(yǎng)殖結(jié)束后每個(gè)個(gè)體體重與體長數(shù)據(jù)。檢測養(yǎng)殖群體中具有多態(tài)性的snp位點(diǎn)，利用統(tǒng)計(jì)軟件spss對snp位點(diǎn)在30尾即墨養(yǎng)殖大瀧六線魚在不用基因型個(gè)體間的性狀進(jìn)行多重比較(lsd)，去除出現(xiàn)少于3次的基因型，保證每個(gè)位點(diǎn)分析價(jià)值。4個(gè)與生長相關(guān)的單核苷酸多態(tài)性位點(diǎn)pic值范圍從0.3457～0.6084，平均為0.4600，有限等位基因數(shù)范圍從1.8000～3.0958，平均為0.8682；觀測雜合度0.2933～0.5115，平均值0.4089；期望雜合度0.4444～0.6770，平均值0.5554。4個(gè)位點(diǎn)均未偏離h-w平衡(p＞0.05)，最小等位基因頻率均大于0.05，說明4個(gè)位點(diǎn)均未受到非隨機(jī)交配的影響，且突變率較高，適合進(jìn)行遺傳學(xué)方面的研究，具有進(jìn)一步研究的價(jià)值。與生長發(fā)育具有顯著性關(guān)系的snp位點(diǎn)及其兩翼序列如表1所示。位點(diǎn)單核苷酸多態(tài)性序列正向引物序列反向引物序列hos10seqidno：1seqidno：5seqidno：6hos13seqidno：2seqidno：7seqidno：8hos14seqidno：3seqidno：9seqidno：10hos17seqidno：4seqidno：11seqidno：124個(gè)單核苷酸多態(tài)性位點(diǎn)的多樣性結(jié)果如表2所示。具有顯著性的位點(diǎn)計(jì)算結(jié)果如表3所示，hos10中tt基因型的體長增長比較顯著。hos13中g(shù)g基因型與其他基因型相比，與體重增量有著較為顯著的關(guān)系。位點(diǎn)hos14中，體重增量上，位點(diǎn)缺失與其他基因型都沒有明顯差距：基因型tt比基因型cc、ct體重增長更快，差別極其顯著；體長增量上，與其余幾種基因型相比，tt基因型具有極其顯著的影響。hos17位點(diǎn)中g(shù)a基因型的體長增量明顯高于aa基因型，具有極顯著的影響?？傮w來看hos13對于體重增長有著顯著性影響，hos10與hos17對體長增長有著顯著性影響，hos14則對體重體長增長都用影響。一個(gè)標(biāo)記與多個(gè)性狀相關(guān)聯(lián)可能與體重與體長本身就具有一定聯(lián)系有關(guān)。開發(fā)出的多態(tài)性位點(diǎn)與生長發(fā)育相關(guān)的位點(diǎn)在今后大瀧六線魚的生長發(fā)育、人工繁殖等研究提供參考價(jià)值。表3不同位點(diǎn)的不同基因型個(gè)體間生長性狀的增量及多重性比較注：同列中不同字母數(shù)值間差異顯著(p＜0.05)，大寫字母表示差異及其顯著(p＜0.01)。三、單核苷酸多態(tài)性位點(diǎn)引物本發(fā)明根據(jù)基于簡化基因組測序(rad)利用二代測序技術(shù)所開發(fā)出的snp文庫，根據(jù)blast中生長相關(guān)堿基序列的比對結(jié)果，并根據(jù)相關(guān)位點(diǎn)，設(shè)計(jì)了50對的篩選引物，在大瀧六線魚的20個(gè)個(gè)體中進(jìn)行驗(yàn)證。單核苷酸多態(tài)性引物的驗(yàn)證擴(kuò)增位于兩端已知序列之間的dna區(qū)域。每個(gè)pcr循環(huán)，反應(yīng)混合體系都在模板變性、引物退火和引物延伸三種不同溫度下依序溫育。該過程可使用一種可編程序的溫度熱循環(huán)儀而達(dá)到自動(dòng)化。5μl體系：模板1μl，pcrmastermix4μl(water1.850μl，pcrbuffer0.625μl，mgcl20.325μl，dntp0.1000μl，primermix1.000μl，hotstartaq0.100μl)。將pcr產(chǎn)物用sap(shrimpalkalinephosphatase，堿性磷酸酶)處理，以去除體系中游離的dntps。反應(yīng)程序：94℃30s，[94℃5s，(52℃5s，80℃5s)5cycles，72℃3min]40cycles，4℃保存，產(chǎn)物用樹脂進(jìn)行純化。啟動(dòng)massarraynanodispenserrs1000點(diǎn)樣儀，將樹脂純化后的延伸產(chǎn)物移至384-wellspectrochipbioarray上。將點(diǎn)樣后的spectrochip芯片使用maldi-tof質(zhì)譜儀分析，檢測結(jié)果使用typer4.0軟件獲取原始數(shù)據(jù)及基因分型圖，檢查數(shù)據(jù)文件的完整性和正確性，將結(jié)果保存入相應(yīng)存儲媒介并遞交生物信息室分析。使用軟件popgene32進(jìn)行遺傳指標(biāo)(等位基因(na：observednumberofalleles)、有效等位基因(ne：effectivenumberofalleles)、觀測雜合度(observedheterozygosity)、期望雜合度(expectedheterozygosity)、香農(nóng)遺傳指數(shù)(shannon′sinformationindex(i))、nei氏多樣性指數(shù)(nei′sgenediversityindex(nei))、多態(tài)信息含量(polymorphisminformationcontent(pic)))計(jì)算。4個(gè)單核苷酸多態(tài)性位點(diǎn)的多樣性及引物見表4。表4其中tm是單核苷酸多態(tài)性的退火溫度，na為等位基因數(shù)，ne為有效等位基因數(shù)，obs_het觀望雜合度，exp-het為期望雜合度。pic為多態(tài)信息含量，i是香農(nóng)指數(shù)，nei為nei氏多樣性指數(shù)，pic為多態(tài)信息含量。本發(fā)明的4對引物對20個(gè)大瀧六線魚基因組dna進(jìn)行擴(kuò)增，遺傳多樣性分析結(jié)果表明每個(gè)單核苷酸多態(tài)性位點(diǎn)平均等位基因數(shù)為2.75，平均有效等位基因數(shù)是2.3381，平均觀望雜合度0.4089，平均期望雜合度為0.5554。遺傳雜合度(geneticheterozygosity)包括了觀望雜合度和期望雜合度，表示群體某一特定基因位點(diǎn)上雜合子的比例，反應(yīng)了群體遺傳變異程度的大小。觀測雜合度是指一個(gè)群體中觀察到的雜合子所占的比例；期望雜合度則是指在哈溫平衡假設(shè)下雜合度的期望值，也稱為基因多樣度(genediversity)，僅取決于等位基因頻率。雜合度的數(shù)值越低，遺傳一致性越高，遺傳多樣性越低，4個(gè)位點(diǎn)的雜合度處于相對較高的水平。4個(gè)位點(diǎn)均未偏離hw平衡，偏離hw平衡的原因可能是由于這些位點(diǎn)在繁殖時(shí)會受到選擇的壓力，進(jìn)而影響不同等位基因產(chǎn)生的頻率，或是由于非隨機(jī)交配所導(dǎo)致的平衡偏離。4個(gè)位點(diǎn)的最小等位基因頻率(minorallelefrequency，maf)均高于0.05突變率相對較高，較小的maf將會是統(tǒng)計(jì)效能降低，從而造成假陰性的結(jié)果，故而造成關(guān)聯(lián)性研究的偏差根據(jù)計(jì)算結(jié)果得出，4對單核苷酸多態(tài)性引物可用于大瀧六線魚遺傳多樣性檢測，數(shù)量性狀遺傳圖譜，遺傳連鎖做圖、以及家系和個(gè)體鑒定。本發(fā)明提供的4個(gè)單核苷酸多態(tài)性位點(diǎn)：seqnoid：1序列第177位等位基因a/t，seqnoid：2序列第228位等位基因g/a，seqnoid：3序列第228位等位基因c/t，seqnoid：4序列第269位等位基因a/g，均與大瀧六線魚生長發(fā)育性狀顯著性相關(guān)。同時(shí)驗(yàn)證了2b-rad技術(shù)開發(fā)大瀧六線魚單核苷酸多態(tài)性位點(diǎn)的可行性。本發(fā)明提供的snp位點(diǎn)能與大瀧六線魚生長發(fā)育相關(guān)，能有效用于大瀧六線魚評估大瀧六線魚生長發(fā)育、大瀧六線魚選種或大瀧六線魚育種。當(dāng)前第1頁12