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

    馬氏珠母貝生長性狀相關(guān)QTL在兩個群體中的驗證

    2017-11-01 06:58:23韋國建管云雁劉文廣林堅士何毛賢
    海洋通報 2017年5期
    關(guān)鍵詞:母貝基因型雜交

    韋國建,管云雁,劉文廣,林堅士,何毛賢

    (1.中國科學(xué)院南海海洋研究所熱帶海洋生物資源與生態(tài)重點實驗室廣東省應(yīng)用海洋生物學(xué)重點實驗室,廣東廣州510301;2.中國科學(xué)院大學(xué),北京100049)

    馬氏珠母貝生長性狀相關(guān)QTL在兩個群體中的驗證

    韋國建1,2,管云雁1,劉文廣1,林堅士1,何毛賢1

    (1.中國科學(xué)院南海海洋研究所熱帶海洋生物資源與生態(tài)重點實驗室廣東省應(yīng)用海洋生物學(xué)重點實驗室,廣東廣州510301;2.中國科學(xué)院大學(xué),北京100049)

    對兩個通過遺傳圖譜構(gòu)建定位的馬氏珠母貝(Pinctada ucata)生長性狀相關(guān)的QTL在兩個不同群體中進(jìn)行驗證分析。結(jié)果顯示:QTL471489在深圳群體中與殼高、殼寬顯著相關(guān)(P<0.05),其CT基因型個體的殼高、殼寬和體質(zhì)量顯著大于CC基因型,CT為優(yōu)勢基因型;而在雜交群體中該QTL與生長性狀關(guān)聯(lián)不顯著(P>0.05),但CT基因型的生長性狀值仍最大。QTL410206在深圳群體中與生長性狀關(guān)聯(lián)不顯著(P>0.05),而在雜交群體中該QTL與殼高、殼長、殼寬和體質(zhì)量顯著相關(guān)(P<0.05),AA為優(yōu)勢基因型?;谶@兩個QTL在這2個群體中分別構(gòu)建二倍型,發(fā)現(xiàn)在深圳群體中C1(CTAT)型和雜交群體中D1(CTAA)、D2(TTAA)、D3(CCAT)、D4(TTAT)、D5(CTAA)型對馬氏珠母貝的生長最為有利,可作為優(yōu)選的基因型組合。

    馬氏珠母貝;數(shù)量性狀位點;生長性狀;驗證

    水產(chǎn)養(yǎng)殖動物的大多數(shù)經(jīng)濟(jì)性狀不僅由多個基因控制,而且容易受環(huán)境等因素的影響,如體質(zhì)量、體長和生長率等,它們都是數(shù)量性狀(Gje原drem et al,2010)。挖掘數(shù)量性狀位點(Quantita原tive trait locus,QTL)的最終目的就是利用與目標(biāo)經(jīng)濟(jì)性狀緊密連鎖的分子標(biāo)記,通過對分子標(biāo)記基因型的檢測,就能獲知目標(biāo)基因的選擇,發(fā)展分子標(biāo)記輔助選擇育種(Molecular marker-assistedse原lection,MAS),即實現(xiàn)基因型選擇育種,從而促進(jìn)對性狀的遺傳改良(Naish et al,2008)。MAS與傳統(tǒng)的育種方法相比,它可以大大提高選擇的效率和可靠性(Ozaki et al,2012),因此已成為遺傳育種研究的主要發(fā)展方向。盡管國內(nèi)外在水產(chǎn)養(yǎng)殖動物的遺傳作圖、QTL定位以及性狀關(guān)聯(lián)分析方面取得了較大的進(jìn)展,但是基于QTL結(jié)果進(jìn)行指導(dǎo)育種的研究還很少,其主要原因就是定位到的QTL很少被鑒定,這些QTL的真實性和通用性還不清楚,因而限制了其應(yīng)用。因此對QTL定位結(jié)果進(jìn)行驗證和鑒定是實施MAS的重要前提(Yue,2014)。目前僅在牙鲆(Paralichthys olivaceus)(Fuji et al,2007,Wang et al,2014)、大西洋鮭魚(Salmo salar)(Moen et al,2009;Tsai et al,2015)、虹鱒(On原corhynchus mykiss)(Wringe et al,2010;Baerwald et al,2011;Vallejo et al,2014)、尖吻鱸(Lates calcar原ifer)(Wang et al,2008)、大黃魚(Larimichthyscrocea)(劉賢德等,2013)、鏡鯉(Cyprinus carpio)(魯翠云等,2012)等少數(shù)魚類的抗病和生長標(biāo)記進(jìn)行了驗證,并且MAS已經(jīng)成功應(yīng)用于牙鲆抗淋巴囊腫?。↙D)(Fuji et al,2007)、大西洋鮭魚抗傳染性胰臟壞死癥(IPN)(Moen et al,2009)、半滑舌鰨(Cynoglossus semilaevis)的早期性別鑒定(Chen et al,2008)、鏡鯉優(yōu)良性狀的培育(孫效文等,2009)等。在海洋貝類中,還極少見有對QTL定位結(jié)果進(jìn)行鑒定和應(yīng)用的報道。

    馬氏珠母貝(Pinctada fucata)是世界上重要的海水珍珠生產(chǎn)貝類,具有極高的經(jīng)濟(jì)價值。與大多數(shù)水產(chǎn)養(yǎng)殖動物一樣,目前馬氏珠母貝的育種方法主要是基于表型的傳統(tǒng)選擇,為了提高其育種效率,有必要開展基因型選擇育種研究。最近,Shi等(2014)利用2b-RAD測序技術(shù)構(gòu)建了馬氏珠母貝高密度遺傳連鎖圖譜,定位到6個與殼長、鉸合線長、體質(zhì)量、珍珠層厚度相關(guān)的QTL。Li等(2014)利用RAD測序技術(shù)獲得了馬氏珠母貝1373個SNP標(biāo)記,構(gòu)建其高密度遺傳連鎖圖譜,定位到與殼高、殼長、殼寬、絞合線等7個生長性狀相關(guān)的QTL,其中定位到QTL471489與軟體部重顯著相關(guān),QTL410206與殼長、總重和軟體部重顯著相關(guān)。但是這些QTL定位結(jié)果都是基于單個家系構(gòu)建遺傳連鎖圖譜獲得的,而且樣本量小,導(dǎo)致獲得的QTL結(jié)果可能并不適用于其他家系或群體。鑒于育種的繁殖群體多為來源于大樣本隨機(jī)群體,因此將獲得的QTL結(jié)果在不同遺傳背景的群體中進(jìn)行驗證十分必要。目前,僅有少數(shù)QTL標(biāo)記被驗證(韋國建等,2015)。本研究用馬氏珠母貝兩個不同遺傳背景的群體對其中的兩個QTL標(biāo)記(471489和410206)進(jìn)行驗證(Li et al,2014),評估QTL標(biāo)記與生長性狀的關(guān)聯(lián)性,為開展基于馬氏珠母貝QTL結(jié)果的分子標(biāo)記輔助選擇育種奠定基礎(chǔ)。

    1 材料與方法

    1.1 實驗群體

    實驗所用的驗證分析群體為深圳群體和1個雜交群體(其親本為深圳群體和湛江養(yǎng)殖群體)。其中深圳養(yǎng)殖群體為8月齡貝,雜交群體為2月齡貝,均養(yǎng)殖于深圳大鵬澳海區(qū)。從深圳養(yǎng)殖群體和雜交群體中分別隨機(jī)抽樣170只和118只。每只貝的殼高、殼長和殼寬用游標(biāo)卡尺測量,精確到0.01mm;體質(zhì)量用電子天平稱量,精確到0.01 g。測量完成后活體取樣,即從活體的貝中取下一小塊的鰓組織置于90%乙醇,于-20益保存。

    1.2 基因組DNA的提取和基因分型

    使用DNA提取試劑盒(Hipure Universal DNA Kit)提取鰓組織的基因組DNA,相關(guān)操作參照試劑盒說明書。提取DNA后,于1.0%的瓊脂糖凝膠電泳檢測其完整性;并用紫外分光光度計檢測其純度及濃度。提取的DNA于–20益保存。

    將獲得的QTL標(biāo)記(471489和410206)所在的短序列作為源序列(Li et al,2014),在已公布的馬氏珠母貝基因組草圖(http://marinegenomics.oist.jp/)和本實驗已有的馬氏珠母貝轉(zhuǎn)錄組數(shù)據(jù)庫中進(jìn)行Blast比對,找出與之匹配的序列。利用基質(zhì)輔助激光解吸電離飛行時間質(zhì)譜法(MALDI-TOF-MS)對QTL標(biāo)記SNP進(jìn)行基因分型,由北京華諾時代科技有限公司完成。

    1.3 數(shù)據(jù)統(tǒng)計與分析

    利用PIC_CALC軟件計算QTL標(biāo)記的多態(tài)信息含量(polymorphism content information,PIC)。利用PopGene1.3軟件計算QTL在群體中期望雜合度(expected heterozygosity,He)和觀測雜合度(observed heterozygosity,Ho),并進(jìn)行哈迪-溫伯格平衡(Hardy-Weinberg equilibrium,HWE)的檢驗。利用SPSS19.0軟件中的一般線性模型(GLM)進(jìn)行標(biāo)記與生長性狀的相關(guān)性分析,模型如下:Yij=怎+Gior Di+eij,其中Yij表示基因型為i的第j個個體的指標(biāo),怎表示群體均值,Gi表示第i種基因型的效應(yīng),Di表示第i種二倍型的效應(yīng),eij表示隨機(jī)誤差變量;基因型對應(yīng)性狀均值之間的差異用Duancan進(jìn)行多重比較分析,顯著性差異設(shè)為P<0.05。

    2 結(jié)果

    2.1 QTL標(biāo)記的遺傳參數(shù)

    利用MALDI-TOF-MS法對深圳群體和雜交群體個體的QTL標(biāo)記SNP進(jìn)行了基因分型,將獲得QTL471489和410206標(biāo)記的基因分型數(shù)據(jù),利用Popgene1.3軟件和PIC_CALC軟件分析這兩個群體的2個QTL標(biāo)記SNP的遺傳多樣性參數(shù)(見表1)。結(jié)果表明,QTL471489和410206標(biāo)記在深圳群體中Ho和He分布區(qū)間分別為0.195 2~0.512 8和0.249 8~0.440 5,在雜交群體中Ho和He分布區(qū)間分別為0.283 1~0.613 5和0.252 4~0.500 1。根據(jù)中度多態(tài)性原則0.25

    表1 QTL標(biāo)記的遺傳參數(shù)

    2.2 QTL標(biāo)記與生長性狀的關(guān)聯(lián)分析

    生長性狀分析表明,深圳群體和雜交群體的殼高、殼長、殼寬和總重均符合正態(tài)分布。經(jīng)SPSS19.0軟件的一般線性模型分析(GLM),QTL471489標(biāo)記在深圳群體中與殼高、殼寬顯著相關(guān)(P<0.05)(表2),基因型多重比較分析發(fā)現(xiàn),CT基因型的個體在殼高、殼寬和體質(zhì)量顯著大于TT基因型,分別高出4.90%、5.13%和10.22%,CT為優(yōu)勢基因型(表3)。而在雜交群體中該QTL與生長性狀關(guān)聯(lián)不顯著(P>0.05)(表2),但是經(jīng)過不同基因型多重比較分析發(fā)現(xiàn),CT基因型個體的生長性狀均值仍最大,CT基因型的個體在殼寬顯著大于CC基因型(表3)。

    經(jīng)過GLM模型分析,結(jié)果表明QTL410206標(biāo)記在雜交群體中與殼高、殼長、殼寬和體質(zhì)量顯著相關(guān)(P<0.05)(表2);AA、AT基因型個體的殼高、殼長、殼寬和體質(zhì)量顯著高于TT基因,AA和AT基因型的生長性狀值沒有顯著性差異,且AA基因型的生長性狀值最大,AA為優(yōu)勢基因型,AA基因型個體在殼高、殼長、殼寬和體質(zhì)量比TT基因型高出8.88%、9.70%、8.29%和27.30%(表3)。而在深圳群體中該QTL與生長性狀關(guān)聯(lián)不顯著(P>0.05),AT基因型個體的生長性狀值大于TT基因型。

    表2 馬氏珠母貝QTL標(biāo)記與生長性狀的GLM分析

    表3 馬氏珠母貝2個QTL標(biāo)記不同基因型多重比較

    2.3 QTL標(biāo)記SNP二倍型構(gòu)建及與生長性狀的關(guān)聯(lián)分析

    為了進(jìn)一步研究這兩個QTL標(biāo)記SNP與生長性狀的作用,基于這兩個QTL在這兩個群體分別構(gòu)建二倍型。經(jīng)過SPSS19.0軟件的GLM分析發(fā)現(xiàn),這兩個QTL標(biāo)記在這兩個群體中不存在交互作用(P>0.05)。但不同基因型在生長性狀上差異顯著,通過基因型多重比較分析發(fā)現(xiàn),在深圳群體中二倍型C1(CTAT)個體在殼高、殼長、殼寬和體質(zhì)量顯著高于C6(TTAT)(P<0.05);在雜交群體中二倍型D1(CTAA)、D2(TTAA)、D3(CCAT)、D4(TTAT)、D5(CTAA)個體在殼高、殼長、殼寬和體質(zhì)量顯著大于二倍型D7(TTTT)、D8(CCTT)(P<0.05),二倍型D1、D2、D3、D4和D5的生長性狀值沒有顯著差異(P>0.05)(表4)。

    表4 QTL標(biāo)記SNP二倍型與馬氏珠母貝生長性狀的關(guān)聯(lián)分析

    3 討論

    Andrew(2010)認(rèn)為在一個理想的群體中各個等位基因應(yīng)當(dāng)是處于Hardy-Weinberg平衡狀態(tài),基因型的分布頻率是穩(wěn)定的,而且其雜合度之間的差異不顯著。本研究中字2檢驗發(fā)現(xiàn)QTL471489標(biāo)記在深圳群體和雜交群體中顯著偏離了Hardy-Weinberg平衡狀態(tài)。分析原因主要是由于人工選擇和非隨機(jī)交配等因素導(dǎo)致的(Waples et al,2015;Baird et al,2015)。QTL471489標(biāo)記SNP的C等位基因頻率在深圳群體為51.83%,而在雜交群體中的頻率只有25.85%,這可能是由于人工選擇育種造成的(Sun et al,2010)。

    研究表明,QTL在不同的遺傳背景、不同的養(yǎng)殖環(huán)境以及不同的發(fā)育時期中驗證應(yīng)該有不完全相同的結(jié)果。Baerwald等(2011)在虹鱒魚(On原corhynchus mykiss)中定位到與眩暈癥緊密連鎖的QTL,在大樣本的實驗群體和驗證家系中,獲得了一個共同QTL區(qū)間位于染色體Omy9,家系間可解釋的表型差異在50~86%間。Wang等(2008)用兩個不同遺傳背景的家系對尖吻鱸(Lates calcarifer)生長性狀QTL進(jìn)行驗證,發(fā)現(xiàn)連鎖群LG2上的QTL標(biāo)記Lca371在兩個家系中與體重具有顯著的相關(guān)性,而其他QTL僅在一個家系表現(xiàn)出與性狀的相關(guān)性。Laghari等(2013)定位到與鯉魚(Cyprinus carpio)3個不同階段體重相關(guān)的QTL,在第6連鎖群上的1個QTL與3個階段的體重都相關(guān)的QTL,而不同階段也有特異的QTL存在。在本研究中,用兩個不同遺傳背景的馬氏珠母貝群體對QTL471489和410206標(biāo)記進(jìn)行驗證,發(fā)現(xiàn)QTL471489標(biāo)記在深圳群體中與殼高、殼寬顯著相關(guān),CT為優(yōu)勢基因型;在雜交群體中與生長性狀關(guān)聯(lián)不顯著,但CT基因型仍是生長性狀值最大。QTL410206在雜交群體中與殼高、殼長、殼寬和體質(zhì)量顯著相關(guān),而在深圳群體中與生長性狀關(guān)聯(lián)不顯著,這可能是由于AA基因型頻率過低,導(dǎo)致關(guān)聯(lián)分析結(jié)果被低估或高估(Ferraz et al,2009)。這兩個QTL標(biāo)記在兩個群體中驗證的結(jié)果不一樣的原因可能是由于不同遺傳背景下QTL的等位基因頻率不一樣以及QTL間上位性效應(yīng)的影響(Frankel et al,1996)。Mackay(2014)在用果蠅(Drosophila melanogaster)表型性狀與基因間的相互作用時發(fā)現(xiàn)上位性效應(yīng)對性狀表型值具有重要的作用。本研究用到兩個群體處于不同的發(fā)育階段,也會影響QTL在群體中驗證的結(jié)果。孫效文等(2009)認(rèn)為同一物種但具有不同遺傳背景的群體或家系應(yīng)有不完全相同的QTL結(jié)果。因此,從已有的研究和本實驗結(jié)果來看,影響QTL驗證結(jié)果的因素是很多的、也很復(fù)雜,如遺傳背景、發(fā)育階段、群體數(shù)量、標(biāo)記間的互作等等,要想篩選到一個能廣泛應(yīng)用的關(guān)聯(lián)QTL標(biāo)記,需進(jìn)行大量的驗證工作。

    從本實驗結(jié)果可以看到,這兩個QTL標(biāo)記均與多個生長性狀相關(guān),同時這兩個標(biāo)記位點又與同一生長性狀相關(guān),表現(xiàn)出一因多效或多因一效的現(xiàn)象,說明這些生長性狀可能是由一個以上的QTL所控制的,該現(xiàn)象也符合QTL的定義(Liu et al,2004)。Li等(2014)定位到QTL471489標(biāo)記與軟體部質(zhì)量顯著相關(guān),QTL410206標(biāo)記與殼長、體質(zhì)量和軟體部質(zhì)量顯著相關(guān)。雖然有一些性狀指標(biāo)在該研究中沒有檢測,但在本研究中顯著關(guān)聯(lián)的性狀仍與這篇文獻(xiàn)的報道存在一定的差異。這是因為這些定位到的QTL一般是通過特定的小樣本量家系作圖獲得的,獲得的QTL結(jié)果可能不適合用于另外一個家系或群體。可見,對QTL定位結(jié)果進(jìn)行驗證十分必要。

    單獨的SNP位點進(jìn)行性狀的關(guān)聯(lián)分析提供的信息往往是有限的。不同SNPs之間存在相互作用,因而通過2個或2個以上的SNP位點基因型組合可以提供更豐富的遺傳信息來彌補(bǔ)單獨SNP的缺陷(He et al,2008;He et al,2012)。Cong等(2014)的研究發(fā)現(xiàn),太平洋牡蠣(Crassostrea gigas)IRR基因上的2個SNP標(biāo)記構(gòu)建的二倍型D3(GGTT)的生長性狀表型值比其他二倍型大,對太平洋牡蠣的生長最為有利。在本研究中,基于這兩個QTL標(biāo)記在深圳群體和雜交群體中構(gòu)建二倍型,發(fā)現(xiàn)在深圳群體中二倍型C1(CTAT)對馬氏珠母貝生長最為有利,可作為優(yōu)選的基因型組合;在雜交群體中二倍型D7(TTTT)、D8(CCTT)顯著小于D1(CTAA)、D2(TTAA)、D3(CCAT)、D4(TTAT)和D5(CTAA)型,在選育過程中可以避免選擇D7、D8型。這兩個QTL的相互作用需要實驗進(jìn)一步驗證。

    致謝:本實驗使用的馬氏珠母貝材料在中國科學(xué)院大亞灣海洋生物綜合試驗站養(yǎng)殖培育。

    Andrews C,2010.The Hardy-Weinberg Principle.Nature Education Knowledge,1(8):65.

    Baerwald M R,Petersen J L,Hedrick R P,et al,2011.A major effect quantitative trait locus for whirling disease resistance identified in rainbow trout(Oncorhynchusmykiss).Heredity,106(6):920-926.

    Baird S J E,2015.Exploring linkage disequilibrium.Molecular Ecology Resources,15(5):1017-1019.

    Chen S L,Deng S P,Ma h Y,et al,2008.Molecular marker-assisted sex control in half-smooth tongue sole(Cynoglossussemilaevis).Aquaculture,283(1):7-12.

    Cong R,Kong L,Yu H,et al,2014.Association between polymorphism in the insulin receptor-related receptor gene and growth traits in the Pacific oyster Crassostreagigas.Biochemical Systematics and Ecolo原gy,54(1):144-149.

    Ferraz J B S,Pinto L F B,Meirelles F V,et al,2009.Association of single nucleotide polymorphisms with carcass traits in Nellorecattle.Genet Mol Res,8(4),1360-1366.

    Frankel W N,Schork N J,1996.Who's afraid of epistasis?.Nature Genet-ics,14(4):371-373.

    Fuji K,Hasegawa O,Honda K,et al,2007.Marker-assisted breeding of a lymphocystisdisease-resistantJapaneseflounder(Paralichthysolivaceus).Aquaculture,272(1):291-295.

    Gjedrem T,Baranski M,2010.Selective breeding in aquaculture:an in-troduction.Springer Science&Business Media,155-159.

    He F,Wen H S,Dong S L,et al,2008.Identification of single nucleotide polymorphism cytochrome P450-c19a and its relation to reproduc-tive traits in Japanese flounder(Paralichthysolivaceus).Aquacul-ture,279(1):177-181.

    He X P,Xia J H,Wang C M,et al,2012.Significant associations of poly-morphisms in the prolactin gene with growth traits in Asian seabass(Latescalcarifer).Animal Genetics,43(2):233-236.

    Laghari M Y,Zhang Y,Lashari P,et al,2013.Quantitativetrait loci(QTL)associatedwithgrowthratetraitincommoncarp(Cyprinuscar-pio).AquacultureInternational,21(6):1 373-1 379.

    Li Y,He M,2014.Genetic mapping and QTL analysis of growth-related traits in Pinctadafucata using restriction-site associated DNA se-quencing.PLoS One,9(11):e111707.

    Liu Z J,Cordes J F,2004.DNA marker technologies and their applica-tions in aquaculture genetics.Aquaculture,238(1):1-37.

    Ma L,Qu Y J,Huai Y T,et al,2011.Polymorphisms identification and associations of KLF7 gene with cattle growth traits.Livestock Sci原ence,135(1):1-7.

    Mackay T F C,2014.Epistasis and quantitative traits:using model organ-isms to study gene-gene interactions.Nature Reviews Genetics,15(1):22-33.

    Moen T,Baranski M,Sonesson A K,et al,2009.Confirmation and finemapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon(Salmosalar):population-level asso-ciations between markers and trait.BMC Genomics,10(1):368.

    NAISH K A,HARD J J,2008.Bridging the gap between the genotype and the phenotype:linking genetic variation,selection and adaptation in fishes.Fish and Fisheries,9(4):396-422.

    Ozaki A,Araki K,Okamoto H,et al,2012.Progress of DNA markerassisted breeding in mariculturedfinfish.Bulletin of Fisheries Research Agency,35(1):31-37.

    Shi Y H,Wang S,Gu Z F,et al,2014.High-density single nucleotide polymorphisms linkage and quantitative trait locus mapping of the pearl oyster,Pinctadafucatamartensii Dunker.Aquaculture,434(3):376-384.

    Sun J,Zhang C,Chen H,et al,2010.Two novel SNPs in coding region of the caprine Fat-inducing transcript gene and their association with growth traits.Molecular Biology Reports,37(1):485-490.

    Tsai H Y,Hamilton A,Guy D R,et al,2015.Verification of SNPs Associated with Growth Traits in Two Populations of Farmed Atlantic Salmon.International Journal of Molecular Sciences,17(1):5.

    Vallejo R L,Palti Y,Liu S,et al,2014.Validation of linked QTL for bacterial cold water disease resistance and spleen size on rainbow trout chromosome Omy19.Aquaculture,432(3):139-143.

    Wang C M,Lo L C,Feng F,et al,2008.Identification andverification of QTL associated with growth traits in twogenetic backgrounds of Barramundi(Latescalcarifer).Animal Genetics,39(1):34-39.

    Wang L,Fan C,Liu Y,et al,2014.A genome scan for quantitative trait loci associated with Vibrio anguillarum infection resistance in Japanese flounder(Paralichthysolivaceus)by bulked segregantan alysis.Marine Biotechnology,16(5):513-521.

    Waples R S,2015.Testing for Hardy-Weinberg proportions:have we lost the plot?.Journal of Heredity,106(1):esu062.

    Wringe B F,Devlin R H,Ferguson M M,et al,2010.Growth-related quantitativetraitlociindomesticandwildrainbowtrout(Oncorhynchusmykiss).BMC Genetics,11(1):63.

    Yue G H,2014.Recent advances of genome mapping and markerassisted selection in aquaculture.Fish and Fisheries,15(3):376-396.

    劉賢德,隋班良,王志勇,等.2013.大黃魚快長相關(guān)微衛(wèi)星標(biāo)記的篩選與驗證.水生生物學(xué)報,37(6):1 036-1 043.

    魯翠云,曹頂臣,李超,等.2012.生長相關(guān)的微衛(wèi)星標(biāo)記在鏡鯉繁殖群體中的驗證.天津農(nóng)學(xué)院學(xué)報,19(3):13-18.

    孫效文,魯翠云,曹頂臣,等.2009.鏡鯉體重相關(guān)分子標(biāo)記與優(yōu)良子代的篩選和培育.水產(chǎn)學(xué)報,33(2):177-181.

    孫效文,魯翠云,賈智英,等.2009.水產(chǎn)動物分子育種研究進(jìn)展.中國水產(chǎn)科學(xué),16(6):981-990.

    韋國建,劉文廣,林堅士,等.2015.馬氏珠母貝兩個與生長性狀相關(guān)QTL的驗證.海洋科學(xué),39(11):13-19.

    Verification of QTL associated with growth traits in two populations of Pearl oyster Pinctada fucata

    WEI Guo-jian1,2,GUAN Yun-yan1,LIU Wen-guang1,LIN Jian-shi1,HE Mao-xian1

    (1.CAS Key Laboratory of Tropical Marine Bio-resources and Ecology,Guangdong Provincial Key Laboratory of Applied Marine Biology,South China SeaInstitute of Oceanology,Chinese Academy of Sciences,Guangzhou 510301,China;2.University of Chinese Academy of Sciences,Beijing 100049,China)

    Quantitative trait loci(QTL)affecting growth traits have been mapped in genetic map of Pinctadafucata.In this study,two of these QTL were verified in two populations with the different genetic backgrounds.The results showed that the QTL471489 was significantly associated with shell height and shell width in Shenzhen population(P<0.05),and oysters with genotype CT were bigger than those with the genotype TT in shell height,shell width and body mass,CT was the superior genotype.However,it was not significantly associated with growth traits in cross-population(P>0.05),and CT was still be the genotype with the maximal average of growth traits.QTL410206 was not significantly associated with growth traits in Shenzhen population(P>0.05),but it was significantly associated with shell height,shell length,shell width and body mass in cross-population(P<0.05),and AA was the superior genotype.Diplotypes were constructed based on the two QTLs,and association analyses revealed that C1(CTAT)in Shenzhen population and D1(CTAA),D2(TTAA),D3(CCAT),D4(TTAT)and D5(CTAA)in cross-population might be the most advantageous diplotypes for growth traits.

    Pinctadafucata;Quantitative trait locus(QTL);growth trait;verification

    Q341

    A

    1001原6932(圓園17)05原園547原07

    10.11840/j.issn.1001-6392.2017.05.010

    2015-09-12;

    2016-11-07

    廣東省漁業(yè)科技與發(fā)展專項(Z2014014;Z2015014);廣東省科技計劃項目(2014B030301064)。

    韋國建(1988-),碩士研究生,從事馬氏珠母貝分子育種。電子郵箱:weiguojian2014@163.com。

    何毛賢,博士,研究員,從事貝類遺傳育種與分子生物學(xué)。電子郵箱:hmx@scsio.ac.cn。

    (本文編輯:袁澤軼)

    猜你喜歡
    母貝基因型雜交
    風(fēng)韻之夢
    鳳凰生活(2018年1期)2018-01-09 17:49:48
    高等植物雜交染色體及其雜交基因表達(dá)的性狀——三論高等植物染色體雜交
    6年生雜交桉無性系對比試驗
    馬氏珠母貝紅色閉殼肌F1代的家系選育及家系評定
    溫度和鹽度對4種殼色馬氏珠母貝谷胱甘肽過氧化物酶基因表達(dá)的影響
    再論高等植物染色體雜交
    西安地區(qū)育齡婦女MTHFRC677T基因型分布研究
    雜交牛
    小說月刊(2014年11期)2014-04-18 14:12:27
    電子背散射衍射法研究馬氏珠母貝珍珠層中文石擇優(yōu)取向
    BAMBI基因敲除小鼠的繁育、基因型鑒定
    久久久久久人妻| 精品一品国产午夜福利视频| 天堂中文最新版在线下载| 99精品久久久久人妻精品| 日韩av免费高清视频| 欧美另类一区| 中文字幕av电影在线播放| 我要看黄色一级片免费的| 久久人人97超碰香蕉20202| 国产又爽黄色视频| 一区二区日韩欧美中文字幕| 日韩人妻精品一区2区三区| 一本色道久久久久久精品综合| 亚洲 欧美一区二区三区| 欧美 亚洲 国产 日韩一| 亚洲欧美一区二区三区国产| 天天躁夜夜躁狠狠躁躁| 国产伦理片在线播放av一区| 女人久久www免费人成看片| 精品国产一区二区三区久久久樱花| 黄频高清免费视频| av又黄又爽大尺度在线免费看| 国产99久久九九免费精品| av国产久精品久网站免费入址| 老司机亚洲免费影院| 建设人人有责人人尽责人人享有的| 日韩精品免费视频一区二区三区| 校园人妻丝袜中文字幕| 精品国产超薄肉色丝袜足j| av视频免费观看在线观看| 免费女性裸体啪啪无遮挡网站| 一区二区av电影网| 色播在线永久视频| 亚洲美女搞黄在线观看| 国产精品一区二区在线观看99| 99精品久久久久人妻精品| 天堂8中文在线网| 亚洲精品第二区| tube8黄色片| 1024香蕉在线观看| 天天添夜夜摸| 丝袜人妻中文字幕| 伦理电影免费视频| 最近2019中文字幕mv第一页| 国产一卡二卡三卡精品 | 亚洲免费av在线视频| 精品第一国产精品| av不卡在线播放| 精品人妻一区二区三区麻豆| 午夜福利在线免费观看网站| 国产精品亚洲av一区麻豆 | 亚洲av欧美aⅴ国产| 日韩中文字幕视频在线看片| a 毛片基地| 日韩 欧美 亚洲 中文字幕| 久久人妻熟女aⅴ| 一级毛片 在线播放| 好男人视频免费观看在线| 久久久久久人妻| 国产精品.久久久| 久久精品久久久久久噜噜老黄| 日本爱情动作片www.在线观看| 亚洲精品在线美女| 天天躁狠狠躁夜夜躁狠狠躁| 制服丝袜香蕉在线| 成年女人毛片免费观看观看9 | 一边摸一边抽搐一进一出视频| 在线亚洲精品国产二区图片欧美| 久久精品熟女亚洲av麻豆精品| 最近中文字幕2019免费版| 国产精品成人在线| 我的亚洲天堂| 国产精品免费视频内射| 中文字幕人妻丝袜一区二区 | av在线app专区| 亚洲人成77777在线视频| 日本wwww免费看| 国产成人啪精品午夜网站| 女性被躁到高潮视频| 国产成人系列免费观看| 男女边摸边吃奶| www.熟女人妻精品国产| 免费黄色在线免费观看| 亚洲国产日韩一区二区| 一本—道久久a久久精品蜜桃钙片| 日韩欧美精品免费久久| 欧美黑人精品巨大| 国产一区二区三区av在线| 午夜免费鲁丝| 女人高潮潮喷娇喘18禁视频| 欧美在线一区亚洲| 日韩制服骚丝袜av| 毛片一级片免费看久久久久| 在线观看一区二区三区激情| 亚洲av中文av极速乱| 高清av免费在线| 大片免费播放器 马上看| 免费看av在线观看网站| 在线观看一区二区三区激情| 国产伦人伦偷精品视频| 一区在线观看完整版| 欧美少妇被猛烈插入视频| 9191精品国产免费久久| 精品国产一区二区三区四区第35| 视频在线观看一区二区三区| 青青草视频在线视频观看| 久久人人爽人人片av| 一区福利在线观看| 亚洲熟女精品中文字幕| 亚洲综合精品二区| 国产xxxxx性猛交| 欧美精品一区二区大全| 各种免费的搞黄视频| 国产成人a∨麻豆精品| 国产深夜福利视频在线观看| 午夜福利免费观看在线| 欧美人与性动交α欧美精品济南到| 国产av国产精品国产| 三上悠亚av全集在线观看| 久久精品人人爽人人爽视色| 又粗又硬又长又爽又黄的视频| 一二三四在线观看免费中文在| 久久亚洲国产成人精品v| 亚洲精品视频女| 欧美国产精品一级二级三级| 亚洲三区欧美一区| 成人午夜精彩视频在线观看| 精品免费久久久久久久清纯 | 日日撸夜夜添| 精品久久久精品久久久| 18禁观看日本| 看十八女毛片水多多多| 午夜av观看不卡| 亚洲精品久久午夜乱码| 老汉色∧v一级毛片| 免费看不卡的av| 久久天躁狠狠躁夜夜2o2o | 久久精品国产综合久久久| 亚洲av成人不卡在线观看播放网 | 久久婷婷青草| 免费人妻精品一区二区三区视频| 自拍欧美九色日韩亚洲蝌蚪91| 丰满少妇做爰视频| 久热爱精品视频在线9| 国产精品国产三级国产专区5o| 欧美亚洲日本最大视频资源| 免费黄网站久久成人精品| 亚洲国产最新在线播放| 免费看av在线观看网站| 免费黄频网站在线观看国产| 国产探花极品一区二区| 亚洲成人av在线免费| 69精品国产乱码久久久| 中文精品一卡2卡3卡4更新| 欧美日韩视频高清一区二区三区二| 婷婷色综合www| 久久精品国产a三级三级三级| h视频一区二区三区| 亚洲欧美精品自产自拍| 色视频在线一区二区三区| 乱人伦中国视频| 国产成人精品久久二区二区91 | 巨乳人妻的诱惑在线观看| 日韩av在线免费看完整版不卡| 黄色怎么调成土黄色| 午夜精品国产一区二区电影| 电影成人av| 综合色丁香网| av在线老鸭窝| 精品免费久久久久久久清纯 | 考比视频在线观看| 国产成人精品福利久久| 亚洲激情五月婷婷啪啪| 国产精品人妻久久久影院| 热re99久久国产66热| 亚洲国产欧美日韩在线播放| 久久狼人影院| 免费观看人在逋| 久久精品久久久久久噜噜老黄| 啦啦啦在线观看免费高清www| 精品一区在线观看国产| 大片免费播放器 马上看| 99热全是精品| 亚洲成人免费av在线播放| 午夜福利,免费看| 狠狠婷婷综合久久久久久88av| 少妇被粗大猛烈的视频| 美女扒开内裤让男人捅视频| 一二三四中文在线观看免费高清| 老司机影院成人| 国产精品麻豆人妻色哟哟久久| 90打野战视频偷拍视频| 久久久久视频综合| 一区二区三区激情视频| 国产有黄有色有爽视频| av国产精品久久久久影院| 天天躁狠狠躁夜夜躁狠狠躁| 精品国产乱码久久久久久小说| 欧美 亚洲 国产 日韩一| 欧美激情 高清一区二区三区| 久久这里只有精品19| 国产精品一区二区精品视频观看| 女性生殖器流出的白浆| 美女扒开内裤让男人捅视频| 精品人妻一区二区三区麻豆| av免费观看日本| 精品卡一卡二卡四卡免费| 亚洲av成人不卡在线观看播放网 | 最近中文字幕高清免费大全6| 黄色视频在线播放观看不卡| 啦啦啦 在线观看视频| 各种免费的搞黄视频| 在线观看www视频免费| 在线天堂最新版资源| 成年av动漫网址| 亚洲精品一区蜜桃| 国产熟女午夜一区二区三区| 女性被躁到高潮视频| 在线精品无人区一区二区三| 国产成人啪精品午夜网站| 在线亚洲精品国产二区图片欧美| 欧美亚洲 丝袜 人妻 在线| 精品国产国语对白av| 日韩视频在线欧美| av.在线天堂| av又黄又爽大尺度在线免费看| 无限看片的www在线观看| av在线app专区| 777米奇影视久久| 精品少妇久久久久久888优播| 在线免费观看不下载黄p国产| 久久久久久人人人人人| 欧美黑人欧美精品刺激| 日本av免费视频播放| 亚洲精品国产色婷婷电影| 亚洲三区欧美一区| 欧美97在线视频| 51午夜福利影视在线观看| 99精品久久久久人妻精品| 国产亚洲最大av| 日韩一区二区三区影片| 男女免费视频国产| 天天躁日日躁夜夜躁夜夜| 欧美av亚洲av综合av国产av | 丝袜美足系列| 午夜福利一区二区在线看| 日韩一卡2卡3卡4卡2021年| 下体分泌物呈黄色| 精品人妻在线不人妻| av福利片在线| 无限看片的www在线观看| 精品国产一区二区三区四区第35| av在线观看视频网站免费| 国精品久久久久久国模美| 久久久久久人妻| 免费观看av网站的网址| 亚洲国产欧美在线一区| 亚洲精品第二区| 国产精品二区激情视频| 日本色播在线视频| 黑人欧美特级aaaaaa片| a 毛片基地| 国产野战对白在线观看| 国产男女超爽视频在线观看| 99香蕉大伊视频| 日日爽夜夜爽网站| 麻豆乱淫一区二区| 国产一区有黄有色的免费视频| 久久女婷五月综合色啪小说| 午夜福利,免费看| 日韩大码丰满熟妇| 男女下面插进去视频免费观看| 国产视频首页在线观看| 婷婷成人精品国产| 亚洲欧美成人精品一区二区| 国产精品 国内视频| 国产一区二区 视频在线| 国产乱人偷精品视频| 在线看a的网站| 欧美另类一区| 999精品在线视频| 看免费成人av毛片| 国产黄频视频在线观看| 人人澡人人妻人| 韩国高清视频一区二区三区| 成年美女黄网站色视频大全免费| 精品国产乱码久久久久久小说| 两个人看的免费小视频| 婷婷成人精品国产| 精品人妻熟女毛片av久久网站| 国产黄频视频在线观看| 丁香六月欧美| 久久久国产欧美日韩av| 一级毛片黄色毛片免费观看视频| 夫妻午夜视频| 一区二区三区激情视频| 国产精品一区二区在线观看99| 一级毛片 在线播放| 亚洲精品国产av成人精品| 国产成人系列免费观看| 最近中文字幕2019免费版| 日韩一卡2卡3卡4卡2021年| 一本大道久久a久久精品| 日日啪夜夜爽| 满18在线观看网站| 各种免费的搞黄视频| 国产精品嫩草影院av在线观看| 一级爰片在线观看| 伊人亚洲综合成人网| 国产亚洲午夜精品一区二区久久| 亚洲一区中文字幕在线| 精品一品国产午夜福利视频| 麻豆av在线久日| 国产1区2区3区精品| 日本91视频免费播放| 国产精品偷伦视频观看了| 欧美乱码精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产成人av激情在线播放| 久久精品aⅴ一区二区三区四区| 操美女的视频在线观看| 亚洲av在线观看美女高潮| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 久久99精品国语久久久| 亚洲精品第二区| 国产精品久久久久久人妻精品电影 | 欧美 日韩 精品 国产| 亚洲国产欧美一区二区综合| 国产成人av激情在线播放| 乱人伦中国视频| 汤姆久久久久久久影院中文字幕| 日日摸夜夜添夜夜爱| 精品一区二区三卡| 少妇的丰满在线观看| 人人妻,人人澡人人爽秒播 | 成年女人毛片免费观看观看9 | 999久久久国产精品视频| 日韩av不卡免费在线播放| 欧美 亚洲 国产 日韩一| 精品国产一区二区三区久久久樱花| 少妇猛男粗大的猛烈进出视频| 高清视频免费观看一区二区| 日韩成人av中文字幕在线观看| 亚洲第一青青草原| 天堂8中文在线网| 黄片小视频在线播放| 精品一区在线观看国产| 宅男免费午夜| 成人漫画全彩无遮挡| 巨乳人妻的诱惑在线观看| 日本午夜av视频| 欧美成人午夜精品| 在线观看免费日韩欧美大片| 精品国产超薄肉色丝袜足j| 香蕉丝袜av| 人人妻人人添人人爽欧美一区卜| 91精品三级在线观看| 精品第一国产精品| 亚洲久久久国产精品| 久久热在线av| 在线观看免费高清a一片| 麻豆精品久久久久久蜜桃| 一级,二级,三级黄色视频| 下体分泌物呈黄色| 亚洲精品视频女| 97精品久久久久久久久久精品| 伊人亚洲综合成人网| 欧美日韩av久久| 午夜福利免费观看在线| 青春草亚洲视频在线观看| 黄色 视频免费看| 在线天堂中文资源库| 亚洲av男天堂| 亚洲精品中文字幕在线视频| 欧美日韩综合久久久久久| 国产男人的电影天堂91| 好男人视频免费观看在线| tube8黄色片| 久久99精品国语久久久| 亚洲精品国产av成人精品| 国产精品久久久人人做人人爽| 亚洲av成人精品一二三区| 97精品久久久久久久久久精品| 欧美最新免费一区二区三区| 亚洲精品国产色婷婷电影| 极品少妇高潮喷水抽搐| 超色免费av| 欧美国产精品va在线观看不卡| 日本黄色日本黄色录像| 国产精品免费大片| 嫩草影视91久久| 亚洲精品久久久久久婷婷小说| 美国免费a级毛片| 日本黄色日本黄色录像| 国产激情久久老熟女| 丝瓜视频免费看黄片| 十八禁人妻一区二区| 在线 av 中文字幕| 亚洲在久久综合| √禁漫天堂资源中文www| 十八禁高潮呻吟视频| 成人黄色视频免费在线看| 中文乱码字字幕精品一区二区三区| 啦啦啦 在线观看视频| 国产男女内射视频| 日日啪夜夜爽| 亚洲精品日韩在线中文字幕| 久久精品国产综合久久久| av在线老鸭窝| av福利片在线| 免费黄色在线免费观看| 在线天堂最新版资源| 看免费av毛片| 国产精品国产av在线观看| 亚洲免费av在线视频| 如何舔出高潮| a级毛片在线看网站| 人妻人人澡人人爽人人| 成年人免费黄色播放视频| 亚洲欧洲精品一区二区精品久久久 | 欧美变态另类bdsm刘玥| 久久久精品国产亚洲av高清涩受| 国产日韩欧美视频二区| 中文字幕人妻丝袜制服| 亚洲精品久久午夜乱码| av视频免费观看在线观看| 成年女人毛片免费观看观看9 | 又粗又硬又长又爽又黄的视频| 日日啪夜夜爽| 99香蕉大伊视频| 99热全是精品| 亚洲av欧美aⅴ国产| 色婷婷久久久亚洲欧美| 亚洲国产欧美一区二区综合| 国产无遮挡羞羞视频在线观看| 精品国产乱码久久久久久小说| 国产乱人偷精品视频| 精品久久蜜臀av无| 丝袜在线中文字幕| 亚洲av中文av极速乱| 午夜福利一区二区在线看| 一级片'在线观看视频| 美女中出高潮动态图| 欧美精品av麻豆av| 精品少妇黑人巨大在线播放| 亚洲成人手机| 男女之事视频高清在线观看 | 中国国产av一级| 久久久久精品人妻al黑| 高清不卡的av网站| 最近中文字幕2019免费版| 久久久久网色| 欧美亚洲 丝袜 人妻 在线| 国产精品蜜桃在线观看| av片东京热男人的天堂| 日韩电影二区| 国产高清国产精品国产三级| 欧美老熟妇乱子伦牲交| 在线观看免费高清a一片| 免费观看人在逋| 国产片特级美女逼逼视频| 日韩熟女老妇一区二区性免费视频| 最近中文字幕高清免费大全6| 99国产精品免费福利视频| 亚洲成人手机| 精品福利永久在线观看| 久久久精品区二区三区| 精品一区二区免费观看| 久热这里只有精品99| 国产精品三级大全| 日韩精品免费视频一区二区三区| 亚洲一级一片aⅴ在线观看| 亚洲精品国产av蜜桃| 99精品久久久久人妻精品| 999精品在线视频| 日本色播在线视频| 久久久久久久大尺度免费视频| 男人操女人黄网站| 欧美亚洲日本最大视频资源| 国产一区二区三区av在线| 欧美亚洲日本最大视频资源| 最新在线观看一区二区三区 | 国产精品久久久久久精品电影小说| 免费高清在线观看视频在线观看| 国产深夜福利视频在线观看| 国产男女内射视频| avwww免费| 青青草视频在线视频观看| 亚洲男人天堂网一区| 午夜免费鲁丝| 国产探花极品一区二区| 最近中文字幕高清免费大全6| 国产一卡二卡三卡精品 | 一区二区三区激情视频| 精品国产乱码久久久久久小说| 欧美精品高潮呻吟av久久| 国产精品一区二区在线不卡| 亚洲成色77777| 9热在线视频观看99| a级毛片在线看网站| 免费黄色在线免费观看| 欧美日韩一级在线毛片| 精品一区二区三卡| 看免费av毛片| 亚洲国产精品一区二区三区在线| 亚洲精品国产区一区二| 最近最新中文字幕免费大全7| 亚洲欧美色中文字幕在线| 美女扒开内裤让男人捅视频| 亚洲欧美成人精品一区二区| 99九九在线精品视频| 男女边摸边吃奶| 免费黄网站久久成人精品| 亚洲精品国产一区二区精华液| 国产一卡二卡三卡精品 | 国产精品蜜桃在线观看| 日韩精品有码人妻一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久综合国产亚洲精品| 肉色欧美久久久久久久蜜桃| 国产精品国产三级国产专区5o| 亚洲国产av影院在线观看| 青春草视频在线免费观看| 最近中文字幕高清免费大全6| 男女之事视频高清在线观看 | 亚洲av在线观看美女高潮| 18禁裸乳无遮挡动漫免费视频| 丝袜喷水一区| 免费不卡黄色视频| 亚洲伊人色综图| 国产成人欧美在线观看 | 高清视频免费观看一区二区| av国产久精品久网站免费入址| 中文天堂在线官网| 久久久久久人人人人人| 一级片'在线观看视频| 99国产精品免费福利视频| 咕卡用的链子| 久久天躁狠狠躁夜夜2o2o | 国产亚洲精品第一综合不卡| 亚洲一级一片aⅴ在线观看| 午夜影院在线不卡| 美女扒开内裤让男人捅视频| 人人妻人人爽人人添夜夜欢视频| 一二三四中文在线观看免费高清| 日韩熟女老妇一区二区性免费视频| 亚洲国产最新在线播放| 国产黄色免费在线视频| 日韩制服丝袜自拍偷拍| av网站在线播放免费| 欧美人与善性xxx| 国语对白做爰xxxⅹ性视频网站| 亚洲av中文av极速乱| 老熟女久久久| svipshipincom国产片| 啦啦啦 在线观看视频| 久久精品亚洲熟妇少妇任你| 亚洲激情五月婷婷啪啪| 侵犯人妻中文字幕一二三四区| 一级,二级,三级黄色视频| 亚洲欧美日韩另类电影网站| 1024视频免费在线观看| 一二三四在线观看免费中文在| 老鸭窝网址在线观看| 久久精品亚洲熟妇少妇任你| 亚洲欧美一区二区三区黑人| 最近手机中文字幕大全| 欧美xxⅹ黑人| 中文乱码字字幕精品一区二区三区| 国产福利在线免费观看视频| 国产精品久久久久久精品电影小说| 老鸭窝网址在线观看| 韩国av在线不卡| 久久精品亚洲av国产电影网| 国精品久久久久久国模美| 黄色视频不卡| 午夜91福利影院| 亚洲av日韩在线播放| 亚洲五月色婷婷综合| av免费观看日本| 狠狠精品人妻久久久久久综合| 欧美日韩福利视频一区二区| 欧美成人精品欧美一级黄| 在线观看免费午夜福利视频| 丝瓜视频免费看黄片| 18禁裸乳无遮挡动漫免费视频| 男人添女人高潮全过程视频| 欧美少妇被猛烈插入视频| 高清视频免费观看一区二区| 啦啦啦在线免费观看视频4| 亚洲成色77777| 国产精品熟女久久久久浪| 色精品久久人妻99蜜桃| 麻豆乱淫一区二区| 99香蕉大伊视频| 高清在线视频一区二区三区| 各种免费的搞黄视频| 午夜福利一区二区在线看| 国产成人av激情在线播放| 色网站视频免费| 曰老女人黄片| 亚洲熟女精品中文字幕| videos熟女内射| 伦理电影免费视频| 国产在线免费精品| 亚洲,一卡二卡三卡| 欧美日韩一区二区视频在线观看视频在线| 国产熟女午夜一区二区三区| 飞空精品影院首页| 爱豆传媒免费全集在线观看| 亚洲精华国产精华液的使用体验| 啦啦啦中文免费视频观看日本| 国产有黄有色有爽视频| 精品久久久久久电影网|