• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    “鲆優(yōu)1號(hào)”牙鲆和親本群體的表型及遺傳性狀變異分析

    2014-06-01 12:30:10田永勝齊文山姜靜王磊張英平劉萬(wàn)軍陳紅林陳松林
    海洋學(xué)報(bào) 2014年6期
    關(guān)鍵詞:牙鲆微衛(wèi)星雜合

    田永勝,齊文山,2,姜靜,2,王磊,3,張英平,2,劉萬(wàn)軍,3,陳紅林,2,陳松林*

    (1.中國(guó)水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所,海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室,山東青島 266071;2.上海海洋大學(xué)水產(chǎn)與生命學(xué)院,上海,201306;3.中國(guó)海洋大學(xué)海洋生命學(xué)院,山東青島 266003)

    “鲆優(yōu)1號(hào)”牙鲆和親本群體的表型及遺傳性狀變異分析

    田永勝1,齊文山1,2,姜靜1,2,王磊1,3,張英平1,2,劉萬(wàn)軍1,3,陳紅林1,2,陳松林1*

    (1.中國(guó)水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所,海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室,山東青島 266071;2.上海海洋大學(xué)水產(chǎn)與生命學(xué)院,上海,201306;3.中國(guó)海洋大學(xué)海洋生命學(xué)院,山東青島 266003)

    牙鲆是我國(guó)及大西洋沿岸天然捕撈和養(yǎng)殖的主導(dǎo)品種,對(duì)其遺傳資源保護(hù)及養(yǎng)殖新品種開(kāi)發(fā)具有重要的意義。本文對(duì)人工培育的牙鲆養(yǎng)殖新品種“鲆優(yōu)1號(hào)”(ZJ)、其父本韓國(guó)牙鲆群體(KS)、母本F0750家系(JX)的18個(gè)表型性狀進(jìn)行測(cè)量,利用SSPS程序進(jìn)行方差分析和多重比較,發(fā)現(xiàn)體長(zhǎng)/體寬,頭長(zhǎng)/眼徑,頭長(zhǎng)/眼間距、眼徑/眼間距,側(cè)線鱗和尾鰭條6個(gè)性狀與父本或母本具有顯著性差異(p<0.05),可做為鑒別“鲆優(yōu)1號(hào)”的標(biāo)志性表型性狀?!蚌覂?yōu)1號(hào)”與母本無(wú)顯著差異(p>0.05)性狀占50%。與父本無(wú)顯著性差異(p>0.05)性狀占33.33%。篩選了24個(gè)多態(tài)性微衛(wèi)星引物對(duì)KS、JX、ZJ和FJ的遺傳多樣性進(jìn)行分析,總共擴(kuò)增出143個(gè)等位基因。4個(gè)群體的平均等位基因數(shù)(Na)大于平均有效等位基因數(shù)(Nae),平均觀測(cè)雜合度大于平均期望雜合度,多態(tài)信息含量(PICa)分析顯示:KS>JX>FJ>ZJ,說(shuō)明人工選育導(dǎo)致了遺傳雜合度的降低。24個(gè)微衛(wèi)星位點(diǎn)在4個(gè)牙鲆群體中的遺傳偏離指數(shù)(D)在-0.759 7到0.188 4之間,在KS、JX、ZJ和FJ群體中分別有10、16、11和17個(gè)位點(diǎn)顯著(p<0.05)或極顯著(p<0.01)偏離Hardy-Weinberg平衡。24個(gè)基因位點(diǎn)的遺傳分化系數(shù)(Fst)在0.076 7~0.258 6之間,平均為0.146 8,表明14.68%的變異來(lái)自種群間,85.32%的變異來(lái)自種群內(nèi)。Nei's遺傳距離顯著,KS和FJ遺傳距離最近(0.282 5),遺傳相似系數(shù)最大(0.753 9);ZJ 與JX遺傳距離最大(0.420 5),遺傳相似系數(shù)最?。?.656 7)。本文研究表明:“鲆優(yōu)1號(hào)”具有明顯的表型和遺傳特征,與父母本及反交群體比較發(fā)生了一定程度的種內(nèi)和種群間遺傳分化,為“鲆優(yōu)1號(hào)”的種質(zhì)鑒別提供了理論依據(jù)。

    牙鲆,鲆優(yōu)1號(hào),表型性狀,衛(wèi)微星標(biāo)記,遺傳多樣性

    1 引言

    牙鲆Paralichthys olivaceus是我國(guó)主要養(yǎng)殖和捕撈魚(yú)類,年產(chǎn)量達(dá)到26 477.5 t(2009年),排在鲆鰈魚(yú)類產(chǎn)量第二位[1]。但是由于牙鲆苗種養(yǎng)殖成活率低、病害多,極大的限制了養(yǎng)殖產(chǎn)業(yè)的發(fā)展。為了改善牙鲆的養(yǎng)殖現(xiàn)狀,培育生長(zhǎng)快、抗病力強(qiáng)的牙鲆新品種,從21世紀(jì)初開(kāi)始不斷從黃渤海海域收集了當(dāng)?shù)匾吧丽胰后w,利用人工選擇方法篩選出抗鰻弧菌牙鲆群體[2],并從日本、韓國(guó)引進(jìn)了牙鲆群體,構(gòu)建了牙鲆育種基礎(chǔ)群體。通過(guò)家系建立和篩選、抗病力測(cè)定、后裔鑒定、家系間雜交等方法[3],于2010年培育出一個(gè)生長(zhǎng)快、成活率高的牙鲆新品種“鲆優(yōu)1號(hào)”[4]。但是對(duì)“鲆優(yōu)1號(hào)”牙鲆的生物學(xué)性狀及分子遺傳性狀的研究至今還未有研究,為了查明“鲆優(yōu)1號(hào)”標(biāo)志性表型特征,以及與其父母本的遺傳區(qū)別,本文對(duì)“鲆優(yōu)1號(hào)”群體、韓國(guó)牙鲆群體、優(yōu)良家系F0750群體及其反交后代群體的表型性狀和微衛(wèi)星標(biāo)記遺傳性狀進(jìn)行了研究。

    表型性狀被大量應(yīng)用于魚(yú)類傳統(tǒng)分類,但近年來(lái)相關(guān)研究報(bào)告極少,在牙鲆形態(tài)性狀研究方面,劉永新和劉金海對(duì)牙鲆家系早期形態(tài)性狀進(jìn)行了分析[5],在魚(yú)類分類和種質(zhì)鑒別方面,近年來(lái)人們大都傾向于利用分子標(biāo)記對(duì)魚(yú)類種群進(jìn)行鑒定和分析。其中微衛(wèi)星標(biāo)記被廣泛應(yīng)用于遺傳圖譜構(gòu)建、QTL定位、親子鑒定、群體遺傳結(jié)構(gòu)分析等方面,趙瑩瑩等利用微衛(wèi)星標(biāo)記對(duì)6個(gè)虹鱒群體的遺傳多樣性進(jìn)行了分析[6],全迎春等利用微衛(wèi)星標(biāo)記對(duì)4個(gè)鯉魚(yú)群體的遺傳多樣性進(jìn)行了分析[7],Kang等利用180個(gè)微衛(wèi)星標(biāo)記和31個(gè)EST分子標(biāo)記構(gòu)建了牙鲆遺傳連鎖圖譜[8],Song等利用1 624個(gè)微衛(wèi)星標(biāo)記構(gòu)建了高密度的牙鲆遺傳連鎖圖譜[9]。Fuji等利用微衛(wèi)星標(biāo)記選育出一個(gè)抗淋巴囊腫病的牙鲆新品種[10]??梢?jiàn)微衛(wèi)星標(biāo)記在魚(yú)類種質(zhì)鑒別、遺傳分析和新品種輔助選育方面具有重要的應(yīng)用價(jià)值,本文試圖將牙鲆育種群體的表型性狀與微衛(wèi)星標(biāo)記性狀結(jié)合起來(lái),對(duì)牙鲆育種群體進(jìn)行分析,為牙鲆新品種的培育和鑒定探索新的技術(shù)途徑。

    2 材料和方法

    2.1 樣本來(lái)源

    實(shí)驗(yàn)用牙鲆4個(gè)群體樣本均取自山東海陽(yáng)黃海水產(chǎn)有限公司,“鲆優(yōu)1號(hào)”樣本(ZJ)為2010年建立的優(yōu)良家系,其父本為韓國(guó)牙鲆群體(KS),母本為F0750家系群體(JX);“鲆優(yōu)1號(hào)”具有生長(zhǎng)快、成活率高的優(yōu)點(diǎn),近年來(lái)已經(jīng)在山東威海、海陽(yáng),福建東山等地大面積推廣養(yǎng)殖。韓國(guó)群體(KS)為2008年從韓國(guó)濟(jì)州島收集的牙鲆群體。F0750家系群體(JX)為2007年建立和篩選出的生長(zhǎng)快優(yōu)良家系。鲆優(yōu)1號(hào)反交家系(FJ)是利用KS為母本,F(xiàn)0750為父本建立的全同胞家系。

    2.2 牙鲆群體表型性狀測(cè)量和分析

    2010年10月對(duì)生長(zhǎng)到2齡的F0750、4齡的韓國(guó)牙鲆及當(dāng)年的“鲆優(yōu)1號(hào)”群體的18個(gè)生物學(xué)性狀進(jìn)行了測(cè)定,生物學(xué)性狀如表2。利用SSPS統(tǒng)計(jì)軟件對(duì)獲得的表型性狀數(shù)據(jù)進(jìn)行單因素方差分析(One-Way ANOVA),采用Student-Newman-Keuls進(jìn)行多重比較和差異顯著性分析,p<0.05表示具有顯著性差異,p>0.05表示差異不顯著,表型性狀值為平均值±標(biāo)準(zhǔn)差。

    2.3 4個(gè)牙鲆群體遺傳材料的收集和DNA的提取

    以上4個(gè)群體各隨機(jī)選取30個(gè)個(gè)體,剪取鰭條后,分別保存于無(wú)水乙醇中,-20℃下保存?zhèn)溆谩;蚪MDNA提取按常規(guī)酚/氯仿提取方法,將提取的DNA調(diào)整濃度至100 ng/μL備用。

    2.4 牙鲆基因組DNA的微衛(wèi)星分析

    從牙鲆微衛(wèi)星引物中挑選24條多態(tài)性較高的引物(見(jiàn)表1),由華大基因合成。采用15μL反應(yīng)體系,含1μL模板DNA,上下游引物各0.4μL,其中1.5 μL 10×PCR緩沖液(TIANGEN,含Mg2+),0.5μL dNTP(CWBIO,2.5 mmol/L),0.1μLTaq DNA聚合酶(TIANGEN,5 U/μL),用滅菌的雙蒸水補(bǔ)充至15 μL。在PCR儀(AB,Alpha-SE)上進(jìn)行PCR反應(yīng),反應(yīng)程序?yàn)椋?5℃預(yù)變性5 min,35個(gè)循環(huán)(95℃變性30 s,按每對(duì)引物的實(shí)際退火溫度反應(yīng)30~40 s,72℃延伸40 s),72℃延伸10 min,4℃保溫。PCR產(chǎn)物經(jīng)6%的聚丙烯酰胺凝膠電泳后,銀染顯色。

    2.5 遺傳數(shù)據(jù)分析與處理

    對(duì)照PBR322分子量標(biāo)準(zhǔn)確定每個(gè)個(gè)體的基因型,利用Popgene軟件計(jì)算4個(gè)群體的等位基因數(shù)(No)、有效等位基因數(shù)(Ne)、平均等位基因數(shù)(Na)、平均有效等位基因數(shù)(Nae)、觀測(cè)雜合度(Ho)和期望雜合度(He)。同時(shí)計(jì)算了各群體的遺傳相似指數(shù)、遺傳距離、近交系數(shù)(Fis)、遺傳分化指數(shù)(Fst)和基因流(Nm)、香農(nóng)指數(shù)(I)和多態(tài)信息含量(PIC),并進(jìn)行Hardy-Weinberg遺傳偏離指數(shù)(D)的計(jì)算。采取UPGMA方法構(gòu)建進(jìn)化樹(shù)。

    3 結(jié)果

    3.1 牙鲆不同群體生物學(xué)性狀比較

    在測(cè)量的18個(gè)生物學(xué)性狀中“鲆優(yōu)1號(hào)”體長(zhǎng)/體寬顯著大于父母本,頭長(zhǎng)/眼徑顯著小于父母本(p<0.05),頭長(zhǎng)/眼間距和眼徑/眼間距顯著大于父母本(p<0.05),側(cè)線鱗和尾鰭條顯著多于父母本(p<0.05)??梢詫⑦@幾個(gè)性狀做為鲆優(yōu)1號(hào)的特征性狀(表2)。

    “鲆優(yōu)1號(hào)”的體長(zhǎng)/尾柄長(zhǎng)、尾柄長(zhǎng)/尾柄寬、背鰭條3個(gè)性狀與母本F0750無(wú)顯著性差異(p> 0.05),但與父本韓國(guó)牙鲆有顯著性差異(p<0.05)。其他性狀與父母本都無(wú)顯著性差異,“鲆優(yōu)1號(hào)”有9個(gè)性狀與母本無(wú)顯著差異(p>0.05),占全部性狀的50%?!蚌覂?yōu)1號(hào)”有6個(gè)性狀與父本無(wú)顯著性差異(p>0.05),占全部性狀的33.33%。說(shuō)明“鲆優(yōu)1號(hào)”與母本的相似性大于父本。

    表1 本文中24對(duì)SSR位點(diǎn)信息

    表2 鲆優(yōu)1號(hào)(ZJ)與F0750家系(JX)、韓國(guó)牙鲆(KS)群體生物學(xué)性狀比較

    續(xù)表2

    續(xù)表2

    3.2 微衛(wèi)星位點(diǎn)的多態(tài)性

    試驗(yàn)所用24對(duì)微衛(wèi)星引物均能在的DNA樣品能穩(wěn)定地?cái)U(kuò)增出相應(yīng)條帶,引物scaffold773_63332在4個(gè)牙鲆群體中的擴(kuò)增結(jié)果見(jiàn)圖1。各微衛(wèi)星引物在24個(gè)微衛(wèi)星位點(diǎn)中有19個(gè)在4個(gè)牙鲆群體中均顯示出了多態(tài)性(見(jiàn)表3),而位點(diǎn)5和21在ZJ群、位點(diǎn)9 和19在FJ群、位點(diǎn)15在JX群中僅有1個(gè)等位基因,主要原因是由于采集樣本中大部分為全同胞家系個(gè)體造成的。24個(gè)微衛(wèi)星位點(diǎn)在KS群中均表現(xiàn)出多態(tài)性。24對(duì)微衛(wèi)星標(biāo)記總共擴(kuò)增出143個(gè)等位基因,等位基因最多的是位點(diǎn)10和21,有9個(gè)等位基因,而位點(diǎn)19等位基因最少,有2個(gè)等位基因,KS群、JX群、ZJ群和FJ群平均等位基因數(shù)(Na)分別為5.833 3、3.583 3、3.291 7和2.791 7;平均有效等位基因數(shù)(Nae)分別為3.693 0、2.738 3、2.276 4和2.428 9(見(jiàn)表4),所有位點(diǎn)的有效等位基因數(shù)都不大于觀察到的等位基因數(shù)。

    圖1 位點(diǎn)3(scaffold773_63332)在4個(gè)群體中的擴(kuò)增結(jié)果

    3.3 牙鲆4個(gè)群體遺傳多樣性

    在24個(gè)微衛(wèi)星位點(diǎn)中,平均觀測(cè)雜合度在0.189 2~1.000 0之間,平均期望雜合度在0.241 6~0.841 5之間(見(jiàn)表5)。KS、JX、ZJ和FJ群體平均觀測(cè)雜合度分別為3.693 0、2.738 3、2.276 4和2.428 9,平均期望雜合度分別為0.690 4、0.577 3、0.490 0 和0.527 2。4個(gè)群體的Shannon's多樣性指數(shù)Ia分別為3.693 0、1.000 0、0.827 0和0.859 9(見(jiàn)表4);24個(gè)位點(diǎn)在4個(gè)牙鲆群體中大部分表現(xiàn)為中度或高度多態(tài),只有位點(diǎn)19在KS群體,2、7、18位點(diǎn)在ZJ群中表現(xiàn)為低度多態(tài)。4個(gè)群體的平均多態(tài)信息含量PICa分別為0.637 0、0.525 9、0.460 1和0.487 9(見(jiàn)表4),4個(gè)群體都屬于高度多態(tài),但“鲆優(yōu)1號(hào)”多態(tài)性相對(duì)較低。以上數(shù)據(jù)表明,經(jīng)過(guò)人工選擇后的3個(gè)群體(JX,ZJ和FJ)基因多樣性要遠(yuǎn)低于韓國(guó)群體(KS),且正交和反交群體的基因多樣性要低于父母本群體。

    表4 4個(gè)牙鲆群體在24個(gè)微衛(wèi)星基因座上的平均等位基因數(shù)(Na)、平均有效等位基因數(shù)(Nae)、觀測(cè)雜合度(Hao)、期望雜合度(Hae)、Nei氏期望雜合度(Nei)、Shannon's多樣性指數(shù)(Ia)以及多態(tài)信息含量(PICa)

    表5 24個(gè)微衛(wèi)星位點(diǎn)在4個(gè)牙鲆群體的觀測(cè)雜合度(Ho)、期望雜合度(He)、Hardy-Weinberg遺傳偏離指數(shù)(D)、遺傳分化信息(Fis,F(xiàn)st,Nm)

    3.4 基因位點(diǎn)遺傳變異

    24個(gè)基因位點(diǎn)的遺傳分化系數(shù)(Fst)和基因流(Nm)分析結(jié)果顯示(表5),F(xiàn)st在0.076 7~0.258 6之間,表明群體間遺傳分化在中等水平以上,其中在位點(diǎn)1遺傳分化最低,在位點(diǎn)12和21遺傳分化最高,各位點(diǎn)的平均遺傳分化指數(shù)為0.146 8,表明僅有14.68%的變異由種群間分化導(dǎo)致,而85.32%的變異來(lái)自種群內(nèi)分化。各位點(diǎn)的Nm在0.716 8~3.008 0之間,其數(shù)值越大,表明遺傳分化程度越低,因此,這同樣說(shuō)明位點(diǎn)1在各群體間的遺傳分化程度最低,在位點(diǎn)12和21遺傳分化最高。

    3.5 Hardy-Weinberg平衡分析

    對(duì)24個(gè)微衛(wèi)星位點(diǎn)分別進(jìn)行等位基因頻率Hardy-Weinberg平衡卡方檢驗(yàn)和遺傳偏離指數(shù)(D)計(jì)算。在KS、JX、ZJ和FJ群體中分別有10、16、11和 17個(gè)位點(diǎn)顯著(p<0.05)或極顯著(p<0.01)偏離Hardy-Weinberg平衡,F(xiàn)0750(JX)和反交群體(FJ)偏離位點(diǎn)最多。只有位點(diǎn)2(scaffold912_69534)在各群體中符合Hardy-Weinberg平衡(表6)。24個(gè)微衛(wèi)星位點(diǎn)在4個(gè)牙鲆群體中的遺傳偏離指數(shù)(D)在-0.759 7到0.188 4之間,其中12個(gè)位點(diǎn)表現(xiàn)為雜合子過(guò)剩(d>0),其余位點(diǎn)均表現(xiàn)為雜合子缺失(d<0)(見(jiàn)表5)。

    表6 4個(gè)牙鲆群體在24個(gè)微衛(wèi)星位點(diǎn)Hardy-Weinberg平衡的卡方檢驗(yàn)

    3.6 牙鲆群體間遺傳距離

    各群體間Nei's遺傳距離在0.282 5~0.420 5,平均0.336 7;相似性指數(shù)在0.656 7~0.753 9,平均0.714 8(見(jiàn)表7)。KS群體和FJ群體的遺傳距離最近為0.282 5,遺傳相似系數(shù)最大為0.753 9;ZJ群體與JX群體的遺傳距離最遠(yuǎn)達(dá)到0.420 5,其相似性遺傳指數(shù)最小0.656 7。根據(jù)各群體之間的遺傳距離,采用UPGMA法對(duì)4個(gè)牙鲆群體進(jìn)行聚類(見(jiàn)圖2)。結(jié)果表明KS群體首先和FJ群體聚為一支,它們?cè)倥cJX聚為一支,最后與ZJ聚合。

    表7 4個(gè)牙鲆群體間的Nei's遺傳距離(下)及相似性指數(shù)(上)

    圖2 牙鲆4個(gè)群體UPGMA聚類圖

    4 討論

    牙鲆主要分布于渤海、黃海、東海、南海及朝鮮、日本、俄國(guó)遠(yuǎn)東沿岸[11],是這些國(guó)家主要捕撈和養(yǎng)殖魚(yú)類。分布在我國(guó)海域的牙鲆被認(rèn)為有兩個(gè)生態(tài)種群,北種群(黃渤海群)1-2月份分布于33°30′~37° 30′N和122°30′~124°0′E,南種群(東海群)1-2月份分布于27°~27°30′N和121°30′~122°30′E[12]。在牙鲆雜交品種“鲆優(yōu)1號(hào)”的培育過(guò)程中,最初利用黃海牙鲆種群,日本牙鲆種群,韓國(guó)濟(jì)州島牙鲆種群,以及利用人工感染鰻弧菌后篩選出的抗病牙鲆群體[2—3]。分布于太平洋西岸的牙鲆在分類上屬于同一個(gè)種,但由于地理上的分割,不同的地理群體在外形特征上表現(xiàn)出一定的差異。濟(jì)州島牙鲆群體的體色近似于沙灘色,日本牙鲆群體體背具有較大的白色花瓣?duì)罨y,而黃海牙鲆體背具有白色點(diǎn)狀花紋。形態(tài)上的差異是長(zhǎng)期地理隔離和環(huán)境選擇的結(jié)果。利用不同地理群體進(jìn)行選擇和雜交育種,使不同地理群體的優(yōu)良性狀得到重組,從而培育出符合人們需求的養(yǎng)殖新品種,是魚(yú)類育種主要途徑之一?!蚌覂?yōu)1號(hào)”的培育充分利用了分布在太平洋西岸不同環(huán)境下牙鲆種群的優(yōu)良性狀,經(jīng)過(guò)大量家系建立和篩選,從而培育出一個(gè)具有生長(zhǎng)快、成活率高的牙鲆優(yōu)良養(yǎng)殖品種。經(jīng)過(guò)幾年來(lái)在山東海陽(yáng)、福建東山島等養(yǎng)殖公司的大量推廣養(yǎng)殖,“鲆優(yōu)1號(hào)”體現(xiàn)出生長(zhǎng)快、成活高的優(yōu)良特點(diǎn),與當(dāng)?shù)仄胀ㄑ丽蚁啾容^,可將苗種養(yǎng)殖生長(zhǎng)速度提高20%以上,成活率提高30%以上[4],深受養(yǎng)殖用戶的好評(píng)。

    表型性狀是判斷一個(gè)物種或品種區(qū)別于其他品種最直觀的特征,在魚(yú)類表型性狀的描述中為了避免不同年齡造成的生長(zhǎng)性狀差異,采用了體長(zhǎng)/體寬、頭長(zhǎng)/眼徑等比例性狀,可有效的減小年齡因素造成的誤差,本文中分析牙鲆4個(gè)群體的表型性狀時(shí)也采用了這一傳統(tǒng)的方法。表型性狀對(duì)于魚(yú)類育種群體的篩選和應(yīng)用可提供直觀的信息,例如對(duì)牙鲆11個(gè)家系后裔不同時(shí)期生長(zhǎng)性狀的分析發(fā)現(xiàn),全長(zhǎng)貢獻(xiàn)率是形態(tài)性狀中最高的[5]。對(duì)大菱鲆4個(gè)地理群體的生長(zhǎng)性能進(jìn)行分析,發(fā)現(xiàn)在生長(zhǎng)速度上依次為法國(guó)、英國(guó)、丹麥和挪威群體[13]。對(duì)90日齡大菱鲆全長(zhǎng)、體長(zhǎng)、頭長(zhǎng)、吻長(zhǎng)、體高、尾柄高、尾柄長(zhǎng)和體質(zhì)量等8項(xiàng)性狀進(jìn)行分析,發(fā)現(xiàn)其中7個(gè)測(cè)量性狀與體重的相關(guān)系數(shù)達(dá)到極顯著水平[14],為育種方案的制定提供了一定的依據(jù)。本文對(duì)“鲆優(yōu)1號(hào)”及其父本韓國(guó)牙鲆、母本F0750(日本牙鲆♂×抗病牙鲆♀的雜交后代)的18個(gè)生物學(xué)性狀進(jìn)行測(cè)量和比較分析,由于韓國(guó)牙鲆和F0750都為篩選出的優(yōu)良親本,在育種中相當(dāng)珍貴,在測(cè)量時(shí)只采集了10尾樣本,但經(jīng)過(guò)分析發(fā)現(xiàn)其可量性狀和可數(shù)性狀的變動(dòng)已很小,說(shuō)明此樣本量已能較準(zhǔn)確的反映群體的表型性狀。研究發(fā)現(xiàn)“鲆優(yōu)1號(hào)”體長(zhǎng)/體寬,頭長(zhǎng)/眼徑,頭長(zhǎng)/眼間距、眼徑/眼間距,側(cè)線鱗和尾鰭條6個(gè)性狀與父本或母本具有顯著性差異,從此為“鲆優(yōu)1號(hào)”的鑒別找到了標(biāo)志性的表型性狀。

    分子標(biāo)記在魚(yú)類種質(zhì)鑒別中的應(yīng)用愈來(lái)愈廣;利用RAPD技術(shù)分析了奧利亞羅非魚(yú)♀×尼羅羅非魚(yú)♂正反雜交子代與親本的遺傳相似系數(shù),顯示其反交后代與母本極其相似,而其正交群體的性狀介于親本之間[15]。4個(gè)尼羅羅非魚(yú)群體Oreochromis niloticus的種群內(nèi)遺傳相似度為0.745 8~0.815 8[16]。本文利用牙鲆18個(gè)表型性狀通過(guò)多重比較和顯著性分析,發(fā)現(xiàn)“鲆優(yōu)1號(hào)”與母本F0750有50%的性狀無(wú)顯著差異(p>0.05),與父本有33.33%的性狀無(wú)顯著性差異(p>0.05),說(shuō)明鲆優(yōu)1號(hào)與母本的相似性大于父本。遺傳距離分析和聚類圖顯示反交群體與父本韓國(guó)牙鲆的遺傳距離最近,其次是母本F0750,而正交群體“鲆優(yōu)1號(hào)”與父本遺傳距離相對(duì)較遠(yuǎn),與母本遺傳距離相對(duì)較近,這一點(diǎn)與表型性狀研究結(jié)果具有相似性。利用24個(gè)微衛(wèi)星標(biāo)記分析顯示“鲆優(yōu)1號(hào)”與母本為遺傳相似性指數(shù)反而小于父本。相反“鲆優(yōu)1號(hào)”親本的反交后代與父本韓國(guó)牙鲆群體的遺傳距離最近,與母本F0750遺傳距離較遠(yuǎn)。是由于分析方法的不同,還是所采用的微衛(wèi)星標(biāo)記與表型性狀相關(guān)性較低造成,有待于進(jìn)一步研究。

    利用16個(gè)微衛(wèi)星位點(diǎn)對(duì)丹東、北戴河、威海、青島、榮成5個(gè)牙鲆群體進(jìn)行遺傳分析,顯示5個(gè)群體遺傳多樣性差異不顯著,群體間基因分化系數(shù)(GST)為0.099 1,各群體之間存在中度遺傳分化[17]。利用10對(duì)微衛(wèi)星引物對(duì)牙鲆普通群體、感病群體和抗病群體進(jìn)行遺傳分析,顯示2個(gè)選擇性養(yǎng)殖群體遺傳多樣性降低[18]。利用9對(duì)多態(tài)性微衛(wèi)星引物對(duì)韓國(guó)沿海5個(gè)野生牙鲆群體和3個(gè)養(yǎng)殖群體的遺傳多樣性進(jìn)行了分析,顯示野生群體與養(yǎng)殖群體具有明顯的遺傳分離,東、西和南部的野生群體聚為一支,而養(yǎng)殖群體聚為一支[19]。利用11個(gè)微衛(wèi)星位點(diǎn)對(duì)日本沿海周圍的7個(gè)牙鲆群體進(jìn)行分析,7個(gè)群體的遺傳分化系數(shù)相當(dāng)?shù)停‵st=0.002 5)[20],可以認(rèn)為是同一個(gè)地理群體。從以上的遺傳分析可以認(rèn)為本文應(yīng)用韓國(guó)牙鲆群體、日本牙鲆群體應(yīng)當(dāng)為代表性群體。而本文中涉及4個(gè)牙鲆群體的遺傳分化系數(shù)為0.076 7~0.258 6,分化系數(shù)相對(duì)較高,主要原因是由于“鲆優(yōu)1號(hào)”父本為韓國(guó)群體,母本F0750的父母本分別為日本群體和從養(yǎng)殖群體中人工篩選而來(lái)的抗病群體,其遺傳背景較為豐富,4個(gè)群體之間產(chǎn)生了中度的遺傳分化,分化水平分別高于日本牙鲆和中國(guó)沿海養(yǎng)殖群體。

    觀測(cè)雜合度(Ho)和期望雜合度(He)是分析種群內(nèi)遺傳變異的2個(gè)重要指標(biāo),日本海區(qū)7個(gè)牙鲆種群的平均等位基因No為15.2~18.2,期望雜合度He為0.74~0.78[20],韓國(guó)沿海牙鲆野生種群平均有效等位基因AR為10.9~16.1,觀測(cè)雜合度Ho為0.820 ~0.888[19]。中國(guó)北部沿海5個(gè)牙鲆養(yǎng)殖群體的觀測(cè)雜合度Ho為0.220 0~0.800 0,期望雜合度He為0.206 1~0.818 7[17],本文4個(gè)群體的觀測(cè)雜合度Ho為0.538 3~0.619 3,期望雜合度He為0.490 0~0.690 4,本文中牙鲆群體的遺傳多樣性低于日本和韓國(guó)自然群體,但遺傳多樣性高于國(guó)內(nèi)沿海養(yǎng)殖群體,主要原因是由于本文中群體是從以上群體選育而來(lái),因此遺傳多樣性低于原產(chǎn)地養(yǎng)殖群體,但經(jīng)過(guò)幾代選育,和本地養(yǎng)殖群體相比遺傳多樣性有所提高。

    等位基因數(shù)(Na)或有效等位基因數(shù)(Ne)易受到外部環(huán)境(選擇、漂變、遺傳瓶頸等)因素的影響[22],適合于評(píng)價(jià)群體遺傳變異[23]。本文4個(gè)牙鲆群體的平均有效等位基因數(shù)(Nae)小于平均等位基因數(shù)(Na),其原因主要是由于群體選育及針對(duì)性的家系選育造成的。Locus5、Locus21在ZJ群體中,Locus9 和Locus19在FJ群體中出現(xiàn)了只有1個(gè)等位基因的現(xiàn)象,主要是由于ZJ和FJ群體為一對(duì)父母本交配產(chǎn)生的全同胞家系群體。而韓國(guó)群體為引進(jìn)的養(yǎng)殖群體,因此利用24個(gè)衛(wèi)微星標(biāo)記進(jìn)行檢測(cè)時(shí),韓國(guó)群體在每個(gè)位點(diǎn)表現(xiàn)為多態(tài)性,而在其ZJ和FJ家系群體中出現(xiàn)了單態(tài)現(xiàn)象。

    多態(tài)信息含量是衡量位點(diǎn)多樣性的較好指標(biāo),一般認(rèn)為在某一群體中,當(dāng)PIC>0.5時(shí)該位點(diǎn)表現(xiàn)為高度多態(tài),當(dāng)0.25<PIC<0.5時(shí)該位點(diǎn)表現(xiàn)為中度多態(tài),當(dāng)PIC<0.25時(shí)該位點(diǎn)表現(xiàn)為低度多態(tài)[21]。中國(guó)北部沿海牙鲆野生群體的30個(gè)位點(diǎn)多態(tài)信息含量為PIC為0.709~0.892,平均為0.832[23]。本文中4個(gè)牙鲆群體的PICa為0.460 1~0.637 0,平均為0.527 7,相對(duì)于野生群體較低,但都表現(xiàn)為高度多態(tài),表明這4個(gè)牙鲆群體種群內(nèi)遺傳變異較大,信息含量較高,多態(tài)信息含量從大到小依次為:韓國(guó)牙鲆、F0750、鲆優(yōu)1號(hào)反交群體和鲆優(yōu)1號(hào)。

    Hardy-Weinberg遺傳偏離指數(shù)(d)反映了Ho和He兩者之間的平衡關(guān)系[24]。北部沿海野生牙鲆群體的遺傳偏離指數(shù)d變化范圍從-0.247到0.512[23],本文中24個(gè)微衛(wèi)星位點(diǎn)在4個(gè)牙鲆群體中的遺傳偏離指數(shù)(D)在-0.759 7到0.188 4之間,其中50%的位點(diǎn)雜合子過(guò)剩(d>0),其余為雜合子缺失(d<0)。4個(gè)群體中F0750(JX)和反交群體(FJ)顯著的偏離位點(diǎn)最多,分別達(dá)到16和17個(gè)。主要原因是由于奠基者效應(yīng)(founder effect)和瓶頸效應(yīng)(bottle effect)造成[25]。韓國(guó)牙鲆(KS)是從國(guó)外引進(jìn),種群數(shù)量有限,選育群體F0750(JX)是一個(gè)生長(zhǎng)快成活率高的家系,父母本為一對(duì)個(gè)體。鲆優(yōu)1號(hào)(ZJ)和其反交群體(FJ)樣本都采自建立的家系群,因而發(fā)生了大部分基因位點(diǎn)的偏離。

    綜上所述,人工培育的牙鲆雜交新品種“鲆優(yōu)1號(hào)”與其父母本相比較,在6個(gè)表型性狀上具有顯著性差異;遺傳性狀分析顯示4個(gè)群體之間產(chǎn)生了中度的遺傳分化,而且種群內(nèi)也產(chǎn)生了較大的遺傳變異。說(shuō)明人工培育的“鲆優(yōu)1號(hào)”在表型和遺傳上具有一定的特質(zhì),這些特點(diǎn)為鑒別“鲆優(yōu)1號(hào)”提供表型和遺傳標(biāo)志,同時(shí)為其優(yōu)良的生產(chǎn)性狀(生長(zhǎng)快、成活率高)提供了一定的遺傳依據(jù)。

    [1] 楊正勇,王春曉.全球視野下中國(guó)鲆鰈類養(yǎng)殖業(yè)的發(fā)展[J].中國(guó)漁業(yè)經(jīng)濟(jì),2009,6(27):115-127.

    [2] 陳松林,田永勝,徐田軍,等.牙鲆抗病群體和家系的建立及其生長(zhǎng)和抗病性能初步測(cè)定[J].水產(chǎn)學(xué)報(bào),2008,32(5):665-674.

    [3] 田永勝,徐田軍,陳松林,等.三個(gè)牙鲆育種群體親本效應(yīng)及遺傳參數(shù)估計(jì)[J].海洋學(xué)報(bào),2009,31(6):119-129.

    [4] 陳松林,田永勝.全國(guó)水產(chǎn)技術(shù)推廣總站-2010水產(chǎn)新品種推廣指南[M].北京:中國(guó)農(nóng)業(yè)出版社,2011:242-268.

    [5] 劉永新,劉海金.牙鲆不同家系早期形態(tài)性狀差異比較[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2008,39(8):82-88.

    [6] 趙瑩瑩,朱曉琛,孫效文.虹鱒6個(gè)養(yǎng)殖群體遺傳多樣性的微衛(wèi)星分析[J].遺傳,2006,28(8):956-962.

    [7] 全迎春,孫效文,梁利群.應(yīng)用微衛(wèi)星多態(tài)分析四個(gè)鯉魚(yú)群體的遺傳多樣性[J].動(dòng)物學(xué)研究,2005,26(6):595—602.

    [8] Kang J H,Kim W J,Lee W J.Genetic linkage map of olive flounder,Paralichthys olivaceus[J].Int J Biol Sci,2008,4(3):143-149.

    [9] Song W T,Pang R Y,Niu Y Z,et al.Construction of high-density genetic linkage maps and mapping of growth-related quantitative trail loci in the Japanese flounder(Paralichthys olivaceus)[J].Plos One,2012,7(11):e50404.

    [10] FujiK,Hasegawa O,HondaK,et al.Marker-assisted breeding of a lymphocystis disease-resistant Japanese flounder(Paralichthysolivaceus)[J]. Aquaculture,2007,272:291-295.

    [11] 田永勝,陳松林,嚴(yán)安生.牙鲆胚胎玻璃化冷凍技術(shù)研究[J].高技術(shù)通訊,2005,15(3):105-110.

    [12] 李思忠,王惠民.中國(guó)動(dòng)物志:硬骨魚(yú)綱-鰈形目[M].北京:科學(xué)出版社,1995:116-123.

    [13] 王新安,馬愛(ài)軍,侯仕營(yíng),等.大菱鲆4個(gè)不同地理群體生長(zhǎng)性能的比較[J].漁業(yè)科學(xué)進(jìn)展,2010,31(1):34-39.

    [14] 于飛,張慶文,孔杰,等.大菱鲆測(cè)量性狀對(duì)體重的影響效果分析[J].海洋水產(chǎn)研究,2008,29(60):33—40.

    [15] 夏德全,曹螢,楊弘,等.羅非魚(yú)雜交F代與親本的遺傳關(guān)系及其雜種優(yōu)勢(shì)的利用[J].水產(chǎn)科學(xué),1999,6(4):29-32.

    [16] 杜誠(chéng),盧邁新,葉星,等.5個(gè)羅非魚(yú)群體的遺傳分析與分子標(biāo)記[J].大連海洋學(xué)院學(xué)報(bào),2005,20(1):25-28.

    [17] 劉海金,朱曉琛,孫效文,等.牙鲆5個(gè)養(yǎng)殖群體的遺傳多樣性分析[J].中國(guó)水產(chǎn)科學(xué),2008,15(1):30-38.

    [18] 劉云國(guó),陳松林,李八方,等.牙鲆選擇性養(yǎng)殖群體遺傳結(jié)構(gòu)的微衛(wèi)星分析[J].高技術(shù)通訊,2006,16(1):94-99.

    [19] Kim W J,KimKK,Han H S,et al.Population structure of the olive flounder(Paralichthys olivaceus)inKorea inferred from microsatellite marker analysis[J].J Fish Biol,2010,76:1958-1971.

    [20] Sekino M,Hara M.Application of microsatellite markers to population genetics studies of Japanese flounder Paralichthysolivaceus[J].Mar Biotechnol(NY),2001,3:572-589.

    [21] Li D,Kang D,Yin Q,et al.Microsatellite DNA marker analysis of genetic diversity in wild common carp(Cyprinuscarpio L.)populations[J].J Genet Genomics,2007,34:984-993.

    [22] Nei M,Maruyama T,Chakraborty R.The bottleneck effect and genetic variability in populations[J].Evolution,1975,29(1):1-10.

    [23] Wang L,Zhang L D,Wan Y M,et al.Isolation of microsatellite markers from Paralichthysolivaceus and its application in genetic structure analysis[J].Yi Chuan,2010,32:1057-1064.

    [24] 李鵬飛,劉萍,柳學(xué)周,等.漠斑牙鲆引進(jìn)種群同工酶的遺傳多態(tài)性分析[J].中國(guó)水產(chǎn)科學(xué),2006,13(1):13-19.

    [25] Suci A,Uthairat N N,WorawutK.Study of genetic diversity of orange spotted grouper,Epinepheluscoioides from Thailand and Indonesia using microsatellite markers[J].Mar Biotechnol,2005,8:17-26.

    Analysis of phenotype and genetic variation of on“Flounder No.1”Paralichthys Olivaceus and their parental populations

    Tian Yongsheng1,Qi Wenshan1,2,Jiang Jing1,2,Wang Lei1,3,Zhang Yingping1,2,Liu Wanjun1,3,Chen Honglin1,2,Chen Songlin1
    (1.Key Laboratory of Sustainable Development of Marine Fisheries,Ministry of Agriculture,Yellow Sea Fisheries ResearchInstitute,Chinese Academy of Fishery Sciences,Qingdao 266071,China;2.College of Fisheries andLife Science,Shanghai Ocean University,Shanghai 201306,China;3.School of Marine life Sciences,0cean University of China,Qingdao 266003,China)

    Japanese flounder(Paralichthys olivaceus)is one of the most popularspecies for fishing and breeding in China and the Atlantic coast,so it is of great importance for the protection of genetic resources and the development of new species.In present study,18 phenotypic traits of“Flounder No.1”(ZJ),a newly breeding species of Japanese flounder,as well as its male(theKorean Japanese flounder stock,KS)and female parents(F0750 family,JX),were measured.The data were compared by ANOVA and SPSS software.The results showed that there were significant differences(p<0.05)between“flounder No.1”and its parents regarding 6 traits,e.g.body length/body width,head length/eye diameter,head length/eye spacing,eye diameter/eye spacing,lateral line scales and caudal fin ray which can be used as typical characters for identification of“Flounder No.1”.In“Flounder No.1”,50%of the traits showed no significant differences(p>0.05)bcompared to its female parent,while the value decreased to 33.3%between“flounder No.1”and its male parent.A total of 24 microsatellite primers were designed to screen the genetic diversity ofKS,JX,ZJand FJ(reciprocal cross offspring),and 143 alleles were identified.For four populations,KS,JX,ZJ and FJ,their mean allelic number(Na)was higher than the effective mean number of alleles(Nae),the mean heterozygosity was higher than the predicted value.The polymorphism information content (PICa)showed thatKS>JX>FJ>ZJ,indicating that the artificial breeding caused the reduction of the genetic heterozygosity.Among the four populations,the genetic deviation index(D)of 24 above mentioned microsatellite loci ranged from-0.759 7 to 0.188 4.And inKS,JX,ZJ and FJ,10,16,11 and 17 loci were significantly(p<0.05)or extremely(p<0.01)deviated from Hardy-Weinberg equilibrium,respectively.The genetic-differentiation-index(Fst)of 24 loci ranged from 0.076 7 to 0.258 6,with the mean value 0.1468,which suggests that 14.68%of the variations result from the interspecies,while 85.32%are intraspecific variations.Furthermore,there was significant difference for Nei's genetic distance among four populations.The nearest genetic distance occured betweenKS and FJ(0.282 5),togehter with the highest genetic similarity coefficient(0.753 9).While the furthest(0.420 5)one was found between ZJ and JX with the lowest(0.656 7)genetic similarity coefficient.In this study,the observation indicated:“flounder No.1”has the obvious phenotypic and genetic characteristics,shows a certain degree of genetic differentiation in intraspecies and interspecies in comparison to its parents and reciprocal cross offspring,which provides the theoretical basis for the germplasm identification of“Flounder No.1”.

    Paralichthys Olivaceus;“Flounder No.1”;phenotypic traits;microsatellite markers;genetic diversity

    Q959.486

    A

    0253-4193(2014)06-0075-12

    2013-04-18;

    2014-02-17。

    “863”高技術(shù)研究發(fā)展計(jì)劃(2012AA10A408);山東省泰山學(xué)者建設(shè)工程專項(xiàng)資助。

    田永勝(1964-),男,甘肅省會(huì)寧縣人,研究員,博士,主要從事魚(yú)類低溫生物學(xué)及遺傳育種研究。E-mail:tianys@ysfri.ac.cn

    *通信作者:陳松林(1960-),男,研究員,博士,主要從事水產(chǎn)生物技術(shù)研究。E-mail:chensl@ysfri.ac.cn

    田永勝,齊文山,姜靜,等.“鲆優(yōu)1號(hào)”牙鲆和親本群體的表型及遺傳性狀變異分析[J].海洋學(xué)報(bào),2014,36(6):75—86,

    10.3969/j.issn.0253-4193.2014.06.010

    Tian Yongsheng,Qi Wenshan,Jiang Jing,et al.Analysis of phenotype and genetic variation of on“Flounder No.1”Paralichthys Olivaceus and their parental populations[J].Acta Oceanologica Sinica(in Chinese),2014,36(6):75—86,doi:10.3969/j.issn.0253-4193.2014.06.010

    猜你喜歡
    牙鲆微衛(wèi)星雜合
    梅花鹿基因組微衛(wèi)星分布特征研究
    2種規(guī)格全雌牙鲆池塘養(yǎng)殖效果分析與評(píng)價(jià)
    甘藍(lán)型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見(jiàn)的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    雌核發(fā)育、野生及養(yǎng)殖牙鲆形態(tài)差異分析
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國(guó)“一箭雙星”成功將“遙感衛(wèi)星二十一號(hào)”與“天拓二號(hào)視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    不同生長(zhǎng)環(huán)境褐牙鲆的外觀性狀與肌肉營(yíng)養(yǎng)成分的比較分析
    雄激素可調(diào)節(jié)的腎臟近端腎小管上皮細(xì)胞靶向雜合啟動(dòng)子的優(yōu)化
    亚洲欧美日韩高清在线视频 | 少妇 在线观看| 久久久久久人人人人人| 亚洲av美国av| 人成视频在线观看免费观看| 亚洲精品乱久久久久久| 精品少妇内射三级| 亚洲精品中文字幕在线视频| 欧美黄色片欧美黄色片| 久久性视频一级片| 成人影院久久| 777久久人妻少妇嫩草av网站| 国产精品麻豆人妻色哟哟久久| 久久久国产精品麻豆| 成人黄色视频免费在线看| 多毛熟女@视频| 国产老妇伦熟女老妇高清| 亚洲专区中文字幕在线| 999久久久国产精品视频| 999久久久国产精品视频| 国产野战对白在线观看| 久久精品熟女亚洲av麻豆精品| 乱人伦中国视频| 免费看av在线观看网站| 精品国产超薄肉色丝袜足j| 91麻豆av在线| 亚洲七黄色美女视频| 午夜福利一区二区在线看| 一级毛片女人18水好多 | 国产1区2区3区精品| av国产精品久久久久影院| 国产成人免费观看mmmm| 波野结衣二区三区在线| 男女午夜视频在线观看| 国产精品免费视频内射| 亚洲欧美日韩高清在线视频 | 国产国语露脸激情在线看| 无遮挡黄片免费观看| 国产亚洲午夜精品一区二区久久| 久久亚洲国产成人精品v| 国产av国产精品国产| 亚洲九九香蕉| 水蜜桃什么品种好| 后天国语完整版免费观看| 日韩制服丝袜自拍偷拍| 一本色道久久久久久精品综合| 一区福利在线观看| 欧美久久黑人一区二区| 欧美日本中文国产一区发布| 亚洲国产欧美网| cao死你这个sao货| 亚洲精品在线美女| 午夜老司机福利片| 亚洲国产看品久久| 天堂8中文在线网| 大陆偷拍与自拍| 老司机亚洲免费影院| 亚洲精品久久久久久婷婷小说| 免费在线观看影片大全网站 | 欧美中文综合在线视频| 久久亚洲国产成人精品v| 99久久综合免费| 人成视频在线观看免费观看| 超碰97精品在线观看| 99久久精品国产亚洲精品| 亚洲精品久久成人aⅴ小说| 我的亚洲天堂| 欧美日韩亚洲综合一区二区三区_| 一二三四在线观看免费中文在| 免费观看a级毛片全部| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区 视频在线| 大陆偷拍与自拍| 麻豆乱淫一区二区| 国产成人一区二区三区免费视频网站 | 免费人妻精品一区二区三区视频| 免费在线观看影片大全网站 | 人人澡人人妻人| 你懂的网址亚洲精品在线观看| 久久久精品区二区三区| 成人国产av品久久久| 成年美女黄网站色视频大全免费| 亚洲中文字幕日韩| 桃花免费在线播放| 国产不卡av网站在线观看| 男男h啪啪无遮挡| 母亲3免费完整高清在线观看| 久久久久久久大尺度免费视频| 可以免费在线观看a视频的电影网站| 熟女av电影| 我要看黄色一级片免费的| 成人黄色视频免费在线看| 美女午夜性视频免费| 高清不卡的av网站| 久久久久久久久久久久大奶| 亚洲五月婷婷丁香| 亚洲国产成人一精品久久久| 国产精品免费大片| 久久鲁丝午夜福利片| 亚洲国产精品成人久久小说| av福利片在线| 人妻 亚洲 视频| 9191精品国产免费久久| 在线观看免费日韩欧美大片| 在线观看一区二区三区激情| 两个人看的免费小视频| 久久女婷五月综合色啪小说| 99热全是精品| 91字幕亚洲| 狂野欧美激情性bbbbbb| 黄片播放在线免费| 99九九在线精品视频| 汤姆久久久久久久影院中文字幕| 在线观看人妻少妇| 人人澡人人妻人| 1024视频免费在线观看| 天天影视国产精品| 脱女人内裤的视频| 国产成人精品久久二区二区91| 欧美 日韩 精品 国产| 纯流量卡能插随身wifi吗| 老司机影院成人| 91精品伊人久久大香线蕉| 国产成人av激情在线播放| 日本猛色少妇xxxxx猛交久久| 性少妇av在线| 日韩,欧美,国产一区二区三区| 亚洲熟女毛片儿| 欧美少妇被猛烈插入视频| 国产伦人伦偷精品视频| 男人舔女人的私密视频| 天堂8中文在线网| 国产精品偷伦视频观看了| 欧美在线黄色| 国语对白做爰xxxⅹ性视频网站| 久久精品亚洲av国产电影网| 日韩 欧美 亚洲 中文字幕| 一级毛片我不卡| 欧美日韩福利视频一区二区| 91精品三级在线观看| 国产精品免费视频内射| 国产成人啪精品午夜网站| 香蕉丝袜av| 91成人精品电影| 久久99热这里只频精品6学生| 欧美亚洲 丝袜 人妻 在线| 成人午夜精彩视频在线观看| 天天躁日日躁夜夜躁夜夜| 视频区图区小说| 在线天堂中文资源库| 黄色 视频免费看| 啦啦啦中文免费视频观看日本| 最新的欧美精品一区二区| 美女中出高潮动态图| 欧美精品啪啪一区二区三区 | 久久国产精品人妻蜜桃| 国产精品久久久久久人妻精品电影 | 欧美少妇被猛烈插入视频| 一级片免费观看大全| 久久国产精品男人的天堂亚洲| 国产成人精品久久二区二区免费| 精品熟女少妇八av免费久了| 亚洲国产日韩一区二区| 午夜精品国产一区二区电影| 一区二区三区四区激情视频| 免费观看a级毛片全部| 成年人免费黄色播放视频| 免费在线观看影片大全网站 | 一级毛片黄色毛片免费观看视频| 国产日韩欧美亚洲二区| 夜夜骑夜夜射夜夜干| 18禁裸乳无遮挡动漫免费视频| 成年美女黄网站色视频大全免费| 国产又爽黄色视频| a 毛片基地| 免费少妇av软件| 国产有黄有色有爽视频| 亚洲成色77777| 波多野结衣av一区二区av| 天天躁夜夜躁狠狠躁躁| 精品国产国语对白av| 十分钟在线观看高清视频www| 免费高清在线观看日韩| 97在线人人人人妻| 欧美人与性动交α欧美软件| 久久久久久亚洲精品国产蜜桃av| 天堂中文最新版在线下载| 人妻 亚洲 视频| 免费看av在线观看网站| 亚洲色图综合在线观看| 少妇人妻久久综合中文| 国产熟女欧美一区二区| 精品视频人人做人人爽| 99热全是精品| 国产av国产精品国产| 性少妇av在线| 99国产精品一区二区蜜桃av | 曰老女人黄片| 男女高潮啪啪啪动态图| 日本av免费视频播放| 人人妻人人澡人人爽人人夜夜| 亚洲成人免费av在线播放| 9热在线视频观看99| 国产日韩一区二区三区精品不卡| 男人操女人黄网站| 欧美日韩综合久久久久久| 日韩电影二区| 天堂俺去俺来也www色官网| 国产av国产精品国产| 婷婷色综合www| 欧美 亚洲 国产 日韩一| 亚洲欧美一区二区三区久久| 国产精品久久久久成人av| 色婷婷久久久亚洲欧美| 建设人人有责人人尽责人人享有的| 欧美国产精品一级二级三级| 婷婷色综合大香蕉| 最近最新中文字幕大全免费视频 | 91麻豆av在线| 9191精品国产免费久久| 精品亚洲乱码少妇综合久久| 90打野战视频偷拍视频| 久久九九热精品免费| 黑丝袜美女国产一区| 在线观看免费视频网站a站| 每晚都被弄得嗷嗷叫到高潮| 欧美人与性动交α欧美软件| 久久综合国产亚洲精品| 最新在线观看一区二区三区 | 爱豆传媒免费全集在线观看| 黄色怎么调成土黄色| av国产久精品久网站免费入址| a级毛片黄视频| 国产无遮挡羞羞视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产免费又黄又爽又色| 看十八女毛片水多多多| 午夜免费观看性视频| av天堂在线播放| 久久久久久免费高清国产稀缺| av不卡在线播放| 丁香六月欧美| 亚洲欧洲精品一区二区精品久久久| 麻豆国产av国片精品| 大话2 男鬼变身卡| 亚洲中文av在线| 国产不卡av网站在线观看| 在线观看人妻少妇| 黄色 视频免费看| 老司机在亚洲福利影院| 精品卡一卡二卡四卡免费| 男女高潮啪啪啪动态图| 日韩大片免费观看网站| 日本色播在线视频| 亚洲三区欧美一区| 亚洲自偷自拍图片 自拍| av国产精品久久久久影院| 国产激情久久老熟女| 精品久久久久久电影网| 久久性视频一级片| 精品免费久久久久久久清纯 | 晚上一个人看的免费电影| 精品福利观看| 一区二区三区乱码不卡18| 午夜免费鲁丝| 亚洲欧洲日产国产| 菩萨蛮人人尽说江南好唐韦庄| 免费人妻精品一区二区三区视频| 亚洲七黄色美女视频| 久久久国产精品麻豆| 欧美黄色片欧美黄色片| 日韩精品免费视频一区二区三区| 亚洲图色成人| 亚洲精品国产区一区二| 9热在线视频观看99| 亚洲欧美清纯卡通| 中文字幕av电影在线播放| 久久毛片免费看一区二区三区| 亚洲国产毛片av蜜桃av| 伊人久久大香线蕉亚洲五| 国产亚洲av高清不卡| 国产片特级美女逼逼视频| 亚洲欧美一区二区三区久久| 一区二区三区四区激情视频| 中国国产av一级| 亚洲精品乱久久久久久| 一区二区av电影网| 久久性视频一级片| 日韩av在线免费看完整版不卡| 少妇裸体淫交视频免费看高清 | 成年女人毛片免费观看观看9 | 欧美+亚洲+日韩+国产| 婷婷成人精品国产| 亚洲三区欧美一区| 侵犯人妻中文字幕一二三四区| 国产片特级美女逼逼视频| 国产在视频线精品| 精品卡一卡二卡四卡免费| 亚洲精品av麻豆狂野| 50天的宝宝边吃奶边哭怎么回事| 男女高潮啪啪啪动态图| 亚洲欧洲日产国产| 久久久国产一区二区| 少妇的丰满在线观看| 亚洲伊人久久精品综合| 又大又黄又爽视频免费| 婷婷色综合大香蕉| 亚洲天堂av无毛| 一区在线观看完整版| 国产精品国产三级专区第一集| 麻豆av在线久日| 亚洲av国产av综合av卡| a级毛片在线看网站| 一级,二级,三级黄色视频| 久久免费观看电影| 久久免费观看电影| 麻豆av在线久日| 18在线观看网站| 大香蕉久久成人网| 久久99精品国语久久久| 18在线观看网站| 免费人妻精品一区二区三区视频| 汤姆久久久久久久影院中文字幕| 日韩一本色道免费dvd| 中国美女看黄片| 欧美亚洲 丝袜 人妻 在线| 一区二区三区精品91| 亚洲色图 男人天堂 中文字幕| 亚洲熟女精品中文字幕| 脱女人内裤的视频| 亚洲伊人色综图| 欧美黄色片欧美黄色片| 天天躁夜夜躁狠狠久久av| 国产成人欧美在线观看 | 狂野欧美激情性xxxx| av网站在线播放免费| www.熟女人妻精品国产| 精品一区二区三区av网在线观看 | 少妇人妻久久综合中文| 在线av久久热| 欧美中文综合在线视频| 国产精品 欧美亚洲| av天堂久久9| 免费在线观看日本一区| 精品少妇内射三级| 高清欧美精品videossex| 亚洲 国产 在线| 国产成人精品久久二区二区91| 欧美国产精品一级二级三级| 女性被躁到高潮视频| 国产日韩欧美在线精品| 黑丝袜美女国产一区| 亚洲av成人精品一二三区| 高清不卡的av网站| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲欧美在线一区二区| 久久精品国产亚洲av高清一级| 只有这里有精品99| 99热国产这里只有精品6| 少妇的丰满在线观看| 亚洲精品国产色婷婷电影| 国产av国产精品国产| 一本大道久久a久久精品| 国产精品一区二区精品视频观看| 免费看av在线观看网站| 亚洲av在线观看美女高潮| 日韩 欧美 亚洲 中文字幕| 亚洲第一青青草原| 久久国产精品大桥未久av| 18禁裸乳无遮挡动漫免费视频| 一级毛片电影观看| 日日爽夜夜爽网站| 欧美国产精品va在线观看不卡| 97精品久久久久久久久久精品| 国产高清视频在线播放一区 | 一边摸一边抽搐一进一出视频| 激情视频va一区二区三区| 日本五十路高清| 欧美变态另类bdsm刘玥| 亚洲中文av在线| 交换朋友夫妻互换小说| 国产精品成人在线| 久久这里只有精品19| 国产精品 欧美亚洲| 国产成人免费无遮挡视频| 男人爽女人下面视频在线观看| 亚洲精品在线美女| 久久久国产欧美日韩av| 黄片小视频在线播放| 丝袜美腿诱惑在线| 一区二区三区精品91| av片东京热男人的天堂| 久久久久精品国产欧美久久久 | 精品免费久久久久久久清纯 | 国产欧美日韩精品亚洲av| 少妇裸体淫交视频免费看高清 | 亚洲成人免费电影在线观看 | 国产成人影院久久av| 国产成人91sexporn| 亚洲精品国产色婷婷电影| netflix在线观看网站| 一边摸一边做爽爽视频免费| 十八禁高潮呻吟视频| 黄片小视频在线播放| 免费女性裸体啪啪无遮挡网站| 久久九九热精品免费| 色综合欧美亚洲国产小说| 国产av国产精品国产| 日韩 亚洲 欧美在线| 精品人妻1区二区| 欧美黑人欧美精品刺激| 亚洲欧美一区二区三区黑人| 波多野结衣一区麻豆| 欧美中文综合在线视频| 中文字幕高清在线视频| 亚洲五月色婷婷综合| 欧美精品高潮呻吟av久久| 看免费成人av毛片| 麻豆av在线久日| 亚洲免费av在线视频| 国产99久久九九免费精品| 黄色视频在线播放观看不卡| 日本五十路高清| 十八禁网站网址无遮挡| 大型av网站在线播放| 性少妇av在线| 亚洲欧美日韩高清在线视频 | 国产91精品成人一区二区三区 | 国产亚洲午夜精品一区二区久久| 99久久99久久久精品蜜桃| 婷婷色av中文字幕| videos熟女内射| 大香蕉久久网| 国产色视频综合| 欧美激情高清一区二区三区| 久9热在线精品视频| 成人影院久久| 欧美日韩av久久| 国产野战对白在线观看| 九草在线视频观看| 青草久久国产| av网站在线播放免费| 国产在视频线精品| 自线自在国产av| 久久久久精品人妻al黑| 这个男人来自地球电影免费观看| 首页视频小说图片口味搜索 | 赤兔流量卡办理| 国产一级毛片在线| 午夜福利,免费看| 人人妻人人澡人人看| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品一区二区免费欧美 | 国产免费视频播放在线视频| 精品国产超薄肉色丝袜足j| 国产午夜精品一二区理论片| 在线看a的网站| 国产精品人妻久久久影院| 一本久久精品| 在线观看人妻少妇| 99国产综合亚洲精品| 亚洲伊人久久精品综合| 亚洲成人手机| 男人舔女人的私密视频| www.av在线官网国产| 制服诱惑二区| 国产精品熟女久久久久浪| 美女高潮到喷水免费观看| 亚洲综合色网址| 国产高清视频在线播放一区 | 最近中文字幕2019免费版| 丰满人妻熟妇乱又伦精品不卡| 18禁观看日本| 欧美在线一区亚洲| 咕卡用的链子| 国产精品麻豆人妻色哟哟久久| 肉色欧美久久久久久久蜜桃| 国产精品九九99| 精品国产乱码久久久久久男人| 中文字幕人妻熟女乱码| 久久99热这里只频精品6学生| 婷婷色av中文字幕| 久久久久精品人妻al黑| 又紧又爽又黄一区二区| av欧美777| 伊人亚洲综合成人网| 久久久国产一区二区| 丝袜脚勾引网站| 另类亚洲欧美激情| 这个男人来自地球电影免费观看| 国产精品一区二区在线不卡| 国产精品一二三区在线看| 高清视频免费观看一区二区| 国产精品成人在线| 亚洲色图综合在线观看| 亚洲 国产 在线| 美女脱内裤让男人舔精品视频| 桃花免费在线播放| www日本在线高清视频| 80岁老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 欧美久久黑人一区二区| 国产高清视频在线播放一区 | 国产精品久久久久久精品古装| 中文字幕人妻丝袜一区二区| 欧美日韩一级在线毛片| 色播在线永久视频| 搡老岳熟女国产| 亚洲精品自拍成人| 精品一区二区三区av网在线观看 | 丝瓜视频免费看黄片| 中文字幕另类日韩欧美亚洲嫩草| 又粗又硬又长又爽又黄的视频| 狂野欧美激情性bbbbbb| 中文字幕人妻丝袜一区二区| 人人妻人人澡人人看| 超碰成人久久| 热99久久久久精品小说推荐| 男人操女人黄网站| 精品国产国语对白av| 男女无遮挡免费网站观看| 嫁个100分男人电影在线观看 | 日韩 欧美 亚洲 中文字幕| 亚洲精品久久久久久婷婷小说| 免费在线观看日本一区| 国产野战对白在线观看| 美女大奶头黄色视频| 精品一区二区三区av网在线观看 | 国产日韩欧美亚洲二区| 欧美精品人与动牲交sv欧美| 色视频在线一区二区三区| 亚洲少妇的诱惑av| 男女边摸边吃奶| 免费在线观看日本一区| 波多野结衣av一区二区av| 婷婷色麻豆天堂久久| 少妇猛男粗大的猛烈进出视频| 日韩视频在线欧美| 国产一区二区三区av在线| 91九色精品人成在线观看| 午夜福利一区二区在线看| 亚洲人成电影免费在线| 老司机影院成人| 亚洲精品美女久久久久99蜜臀 | 建设人人有责人人尽责人人享有的| 国产亚洲欧美在线一区二区| 2018国产大陆天天弄谢| 蜜桃国产av成人99| 欧美xxⅹ黑人| 日本av手机在线免费观看| 国产在线观看jvid| 欧美成狂野欧美在线观看| 夜夜骑夜夜射夜夜干| 亚洲,欧美精品.| 亚洲免费av在线视频| 丰满饥渴人妻一区二区三| 一级毛片黄色毛片免费观看视频| 一本久久精品| 欧美xxⅹ黑人| 亚洲欧美日韩另类电影网站| 日本a在线网址| 女人爽到高潮嗷嗷叫在线视频| 美女高潮到喷水免费观看| 91老司机精品| 欧美老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 国产精品香港三级国产av潘金莲 | 丝袜喷水一区| 别揉我奶头~嗯~啊~动态视频 | 三上悠亚av全集在线观看| 性少妇av在线| 无遮挡黄片免费观看| 国产成人欧美在线观看 | 大片免费播放器 马上看| a级片在线免费高清观看视频| 成人国语在线视频| 午夜精品国产一区二区电影| 久久天躁狠狠躁夜夜2o2o | 熟女av电影| 久久久国产一区二区| 亚洲伊人色综图| 制服诱惑二区| 亚洲专区国产一区二区| 免费观看av网站的网址| 亚洲色图综合在线观看| 热99久久久久精品小说推荐| 国产一区亚洲一区在线观看| 91字幕亚洲| 亚洲美女黄色视频免费看| 国产亚洲一区二区精品| 国产黄频视频在线观看| 女人被躁到高潮嗷嗷叫费观| 在线观看免费午夜福利视频| 欧美成人精品欧美一级黄| 观看av在线不卡| 久久久久久久大尺度免费视频| 久久九九热精品免费| 啦啦啦啦在线视频资源| 制服诱惑二区| 中文字幕亚洲精品专区| 国产国语露脸激情在线看| 极品少妇高潮喷水抽搐| 亚洲精品乱久久久久久| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕亚洲精品专区| 如日韩欧美国产精品一区二区三区| 亚洲专区国产一区二区| 亚洲av欧美aⅴ国产| 国产精品熟女久久久久浪| 欧美日韩亚洲综合一区二区三区_| 国产成人免费无遮挡视频| 久久久久精品人妻al黑| 精品一区二区三区av网在线观看 | 人人妻人人爽人人添夜夜欢视频|