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

    菲律賓蛤仔EST_SSR標(biāo)記與生長(zhǎng)性狀的相關(guān)分析

    2015-03-15 07:36:08牛泓博聶鴻濤朱德鵬閆喜武
    生態(tài)學(xué)報(bào) 2015年6期
    關(guān)鍵詞:蛤仔微衛(wèi)星雜合

    牛泓博,聶鴻濤,朱德鵬,楊 鳳,閆喜武

    大連海洋大學(xué),遼寧省貝類良種繁育工程技術(shù)研究中心, 大連 116023

    菲律賓蛤仔EST_SSR標(biāo)記與生長(zhǎng)性狀的相關(guān)分析

    牛泓博,聶鴻濤,朱德鵬,楊 鳳,閆喜武*

    大連海洋大學(xué),遼寧省貝類良種繁育工程技術(shù)研究中心, 大連 116023

    研究利用20 個(gè)微衛(wèi)星標(biāo)記對(duì)菲律賓蛤仔斑馬蛤F2代家系107 個(gè)個(gè)體進(jìn)行遺傳多樣性分析,并對(duì)標(biāo)記位點(diǎn)與生長(zhǎng)相關(guān)性狀進(jìn)行分析。在20 個(gè)微衛(wèi)星位點(diǎn)共檢測(cè)到41 個(gè)等位基因,各位點(diǎn)等位基因數(shù)為2—3 個(gè),等位基因片段大小為109—430 bp,平均等位基因數(shù)為2.05 個(gè)。平均有效等位基因數(shù)為1.71個(gè),觀測(cè)雜合度平均值為0.504,期望雜合度的平均值為0.431,平均多態(tài)信息含量為0.324。經(jīng)卡方檢驗(yàn),3 個(gè)位點(diǎn)SSR11,SSR164和SSR213的基因型分布顯著偏離了孟德爾定律(P<0.01)。運(yùn)用SPSS 20.0對(duì)20 個(gè)微衛(wèi)星位點(diǎn)與菲律賓蛤仔斑馬蛤家系生長(zhǎng)性狀的相關(guān)性(殼長(zhǎng)、殼寬、殼高和體重)進(jìn)行連鎖顯著性檢驗(yàn)。結(jié)果表明,SSR9位點(diǎn)與殼高存在顯著的相關(guān)關(guān)系(P<0.05),SSR135和SSR164位點(diǎn)與殼寬呈顯著相關(guān)(P<0.05),SSR142位點(diǎn)與體重呈顯著性相關(guān)(P<0.05)。研究結(jié)果可為菲律賓蛤仔的分子標(biāo)記輔助選育提供參考。

    菲律賓蛤仔; 微衛(wèi)星標(biāo)記; 生長(zhǎng)性狀; 相關(guān)分析

    菲律賓蛤仔(Ruditapesphilippinarum)隸屬于軟體動(dòng)物門(Mollusca),雙殼綱(Bivalvia),異齒亞綱(Heterodonta),簾蛤目(Veneroida),簾蛤科(Veneridae),蛤仔屬(Ruditapes),是我國傳統(tǒng)四大養(yǎng)殖貝類之一,廣泛分布在我國南北沿海。據(jù)統(tǒng)計(jì),2010年菲律賓蛤仔世界產(chǎn)量達(dá)到360多萬t[1],我國年產(chǎn)量在300 萬t左右,占世界菲律賓蛤仔產(chǎn)量的90%以上[2]。由于其營養(yǎng)豐富、味道鮮美,近年來市場(chǎng)需求量也日益增大。然而由于人工苗種的累代養(yǎng)殖,菲律賓蛤仔養(yǎng)殖群體的種質(zhì)逐漸下降,養(yǎng)殖產(chǎn)業(yè)迫切需要菲律賓蛤仔的優(yōu)良品種。近年來,菲律賓蛤仔的遺傳改良技術(shù)取得了一定的成績(jī),如選擇育種和雜交育種等[2-6]。傳統(tǒng)育種方法對(duì)目標(biāo)性狀缺乏準(zhǔn)確可靠的早期預(yù)選手段,造成育種盲目性大、效率低。而緊密連鎖的DNA標(biāo)記為實(shí)現(xiàn)水產(chǎn)動(dòng)物目標(biāo)性狀的早期選擇提供了有效途徑。因此,研究與水產(chǎn)動(dòng)物重要經(jīng)濟(jì)性狀緊密連鎖的分子標(biāo)記,特別是對(duì)缺乏遺傳連鎖圖譜的水產(chǎn)動(dòng)物DNA標(biāo)記,就顯得十分重要。

    微衛(wèi)星(microsatellites)又稱簡(jiǎn)單序列重復(fù)(Simple Sequence Repeat, SSR),一般由1—6 bp的簡(jiǎn)單序列重復(fù)排列組成。由于微衛(wèi)星具有多態(tài)性豐富、共顯性遺傳、遵循孟德爾分離定律、易于(Polymerase Chain Reaction) PCR擴(kuò)增等優(yōu)點(diǎn),是分子遺傳學(xué)研究的理想分子標(biāo)記,已廣泛應(yīng)用于水產(chǎn)動(dòng)物的重要經(jīng)濟(jì)性狀相關(guān)分子標(biāo)記篩選、(Quantitative Trait Locus) QTL定位等方面[7-15]。在海洋貝類,蝦夷扇貝[16]和馬氏珠母貝[17]等已報(bào)道了生長(zhǎng)性狀與微衛(wèi)星標(biāo)記的相關(guān)性分析。目前,對(duì)菲律賓蛤仔的研究,主要集中在生物學(xué)、生態(tài)學(xué)和養(yǎng)殖技術(shù)等方面[2]。然而,有關(guān)菲律賓蛤仔重要經(jīng)濟(jì)性狀的分子基礎(chǔ)研究還剛剛起步,標(biāo)記與生長(zhǎng)性狀的相關(guān)分析迄今還未見報(bào)道。

    本實(shí)驗(yàn)利用20 對(duì)微衛(wèi)星標(biāo)記對(duì)菲律賓蛤仔斑馬蛤家系進(jìn)行遺傳多樣性分析,并對(duì)標(biāo)記位點(diǎn)與生長(zhǎng)性狀相關(guān)性進(jìn)行分析,篩選出對(duì)斑馬蛤生長(zhǎng)有利的基因型,實(shí)現(xiàn)分子標(biāo)記和生長(zhǎng)性狀的連鎖分析。為菲律賓蛤仔優(yōu)良品系的進(jìn)一步選育,以及遺傳圖譜的構(gòu)建和QTL定位等提供參考。

    1 材料與方法

    1.1 實(shí)驗(yàn)材料

    實(shí)驗(yàn)所用斑馬蛤材料為2011年大連石河菲律賓蛤仔野生群體中挑選的斑馬蛤?yàn)橛H本,交配產(chǎn)生F1代群體,上選F1代群體中大的個(gè)體,采用巢式設(shè)計(jì)建立了30 個(gè)F1家系,本實(shí)驗(yàn)所用家系為第8 組家系,隨機(jī)選取其中107 個(gè)個(gè)體,對(duì)所有個(gè)體的殼長(zhǎng)、殼寬、殼高和體重進(jìn)行測(cè)量。取部分足尖組織,固定于95%乙醇中用于DNA的提取。

    1.2 DNA提取及微衛(wèi)星分析

    DNA的提取采用苯酚-氯仿-異戊醇法抽提DNA,并用1%瓊脂糖凝膠電泳進(jìn)行檢測(cè)。DNA保存于-20 ℃?zhèn)溆谩CR反應(yīng)體系:反應(yīng)總體積為10 μL,其中DNA模板1 μL、上下游引物(10 μmol/L) 各0.4 μL、easytaq MIX 5 μL、ddH2O 3.2 μL。實(shí)驗(yàn)中所用微衛(wèi)星引物的核心序列及退火溫度見表1。PCR反應(yīng)程序?yàn)椋?4 ℃ 5 min,94 ℃ 40 s,退火溫度40 s,72 ℃ 40 s,30 個(gè)循環(huán);72 ℃ 延伸10 min,4 ℃保存。PCR擴(kuò)增產(chǎn)物經(jīng)12%聚丙烯酰胺凝膠電泳檢測(cè),銀染顯色。凝膠在掃描儀上成像,拍照保存后對(duì)電泳譜帶進(jìn)行分析。等位基因大小用100 bp marker為參照標(biāo)準(zhǔn)對(duì)電泳條帶進(jìn)行判讀。

    表1 菲律賓蛤仔20個(gè)微衛(wèi)星引物序列和退火溫度Table 1 Characteristics of the 20 polymorphic microsatellite loci from Ruditapes philippinarum

    1.3 統(tǒng)計(jì)分析

    用PopGene3.2軟件進(jìn)行數(shù)據(jù)處理、遺傳多樣性分析, 計(jì)算等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、觀測(cè)雜合度(Ho)及期望雜合度(He)。根據(jù)Botstein等[18]的方法可計(jì)算多態(tài)信息含量(PIC)。利用PIC-CALC計(jì)算軟件,根據(jù)等位基因頻率來計(jì)算每個(gè)位點(diǎn)的多態(tài)信息含量,其計(jì)算公式如下:

    式中,n為某一位點(diǎn)上的等位基因數(shù);pi、pj分別為群體中第i和第j個(gè)位點(diǎn)的等位基因頻率,j=i+1。

    用SPSS20.0軟件中的一般線性模型(GLM) 對(duì)菲律賓蛤仔斑馬蛤家系的生長(zhǎng)相關(guān)性狀與20個(gè)微衛(wèi)星標(biāo)記相關(guān)性進(jìn)行分析。其線性模型如下:

    yij=μ+gi+εij

    式中,yij為斑馬蛤某性狀第i個(gè)標(biāo)記第j個(gè)個(gè)體表型值,μ為群體平均值,gi為第i個(gè)標(biāo)記的效應(yīng)值,εij為隨機(jī)誤差。

    2 結(jié)果

    2.1 性狀的表型值分布

    經(jīng)測(cè)量得到的殼長(zhǎng)、殼寬、殼高和體重4 個(gè)生長(zhǎng)性狀,均符合正態(tài)分布,可直接進(jìn)行相關(guān)性分析。正態(tài)分布信息見表2。

    表2 殼長(zhǎng)、殼寬、殼高和體重的表型數(shù)據(jù)及正態(tài)分布統(tǒng)計(jì)檢驗(yàn)Table 2 The phenotypic data of shell lenth, width, height, body weight and normal distribution test

    2.2 微衛(wèi)星位點(diǎn)基因型

    圖1 引物SSR372和SSR373在斑馬蛤F2 家系部分個(gè)體的聚丙烯電泳圖譜Fig.1 Electrophoresis pattern of PCR product of partial individuals in the F2 family of Ruditapes philippinarum amplified with microsatellite marker SSR372 and SSR373

    利用父母本和8個(gè)子代進(jìn)行引物篩選,經(jīng)PCR 擴(kuò)增、聚丙烯酰胺凝膠電泳,從32 對(duì)微衛(wèi)星引物中篩選得到20 對(duì)擴(kuò)增條帶清晰、穩(wěn)定,具有多態(tài)性的位點(diǎn)。檢測(cè)了父母本和107 個(gè)F2子代個(gè)體在20 個(gè)微衛(wèi)星位點(diǎn)的基因型(表3)。SSR372和SSR373位點(diǎn)在菲律賓蛤仔斑馬蛤家系部分個(gè)體中的擴(kuò)增結(jié)果見圖1。

    20 個(gè)微衛(wèi)星位點(diǎn)共檢測(cè)到41 個(gè)等位基因,片段大小為109—430 bp,各位點(diǎn)的等位基因數(shù)2—3 個(gè),平均等位基因數(shù)2.05 個(gè),平均有效等位基因數(shù)為1.71 個(gè)。各位點(diǎn)的期望雜合度、觀測(cè)雜合度、基因型分布和多態(tài)信息含量等值見表3。20 個(gè)位點(diǎn)的期望雜合度平均值為0.431,觀測(cè)雜合度平均值為0.504,多態(tài)信息含量平均值為0.324。在20 個(gè)位點(diǎn)中,有2 個(gè)位點(diǎn)SSR213和SSR261的觀測(cè)雜合度低于期望雜合度,表現(xiàn)出純合子過剩、雜合子不足現(xiàn)象。另外18 個(gè)位點(diǎn)的觀測(cè)雜合度則高于期望雜合度,其中SSR11和SSR164的觀測(cè)雜合度與期望雜合度的差異較大。經(jīng)卡方檢驗(yàn),有17 個(gè)位點(diǎn)的基因型比例符合孟德爾分離定律,另外3 個(gè)位點(diǎn)SSR11,SSR164和SSR213的基因型分布顯著偏離了孟德爾定律(P<0.01)(表3)。

    表3 菲律賓蛤仔20個(gè)微衛(wèi)星位點(diǎn)在斑馬蛤F2家系中的統(tǒng)計(jì)信息Table 3 Statistic information for 20 microsatellite loci in an F2 family of Ruditapes philippinarum

    Na: The number of alleles;Ne: Number of effective alleles;Ho: Observed heterozygosity;He: Expected heterozygosity; PIC: Polymorphism information content;P為基因型分離比的卡方檢驗(yàn);加粗P值代表偏分離的基因型分離比 (P<0.01)

    2.3 微衛(wèi)星標(biāo)記與生長(zhǎng)性狀的相關(guān)性分析

    用SPSS 20.0軟件中的一般線性模型(General Linear Model, GLM),對(duì)微衛(wèi)星位點(diǎn)基因型與斑馬蛤主要生長(zhǎng)性狀(殼長(zhǎng)、殼寬、殼高和體重)之間的相關(guān)性進(jìn)行分析,經(jīng)方差分析檢驗(yàn)呈顯著性差異的位點(diǎn)。在20 對(duì)微衛(wèi)星位點(diǎn)中,共找出4 個(gè)與斑馬蛤生長(zhǎng)相關(guān)的位點(diǎn)(表4)。SSR9位點(diǎn)與殼高呈顯著關(guān)聯(lián),SSR135和SSR164與殼寬顯著相關(guān),SSR142位點(diǎn)與體重呈顯著性相關(guān)。由于相關(guān)分析得到的微衛(wèi)星位點(diǎn)除SSR9號(hào)引物外其他3 個(gè)位點(diǎn)基因型個(gè)數(shù)都小于3,故使用SPSS中compare means過程進(jìn)行不同位點(diǎn)表型值的比較。在SSR9位點(diǎn)上,AA基因型的殼長(zhǎng),殼寬,殼高和體重4 項(xiàng)生長(zhǎng)指標(biāo)平均值均高于其他基因型個(gè)體,說明AA基因型對(duì)斑馬蛤生長(zhǎng)起正效應(yīng)。另外,SSR135位點(diǎn)上的BB基因型、SSR142位點(diǎn)上的AA基因型和SSR164位點(diǎn)上的AB基因型,在4 項(xiàng)生長(zhǎng)指標(biāo)中的平均值均高于各自位點(diǎn)中的其他基因型,說明4個(gè)基因型在各自位點(diǎn)上對(duì)斑馬蛤的生長(zhǎng)具有正效應(yīng)。

    表4 20個(gè)微衛(wèi)星位點(diǎn)不同基因型主要生長(zhǎng)性狀的平均值及多重比較Table 4 Means and multiple comparisons of major growth traits with different genotypes at 20 microsatellites

    上標(biāo)字母表示在同一位點(diǎn)中不同基因型之間差異顯著(P<0.05),無字母表示基因型間差異不顯著,*表示與生長(zhǎng)相關(guān)的位點(diǎn)

    3 討論

    群體的遺傳多樣性主要表現(xiàn)在雜合度、等位基因數(shù)和多態(tài)信息含量三個(gè)方面。雜合度大小可以反映出遺傳結(jié)構(gòu)變異程度的高低,雜合度期望越高,群體的遺傳結(jié)構(gòu)越復(fù)雜[19]。本實(shí)驗(yàn)中,觀測(cè)雜合度和期望雜合度的平均值分別為0.504和0.431,說明實(shí)驗(yàn)家系的遺傳變異和遺傳多樣性較豐富。從等位基因來看,有效等位基因越接近所檢測(cè)到的等位基因的絕對(duì)值,說明等位基因在群體中的分布越均勻[20]。本研究在20 個(gè)微衛(wèi)星位點(diǎn)共檢測(cè)到41 個(gè)等位基因,其中有效等位基因數(shù)為34 個(gè),接近于實(shí)際等位基因數(shù),表明所檢測(cè)的微衛(wèi)星位點(diǎn)的等位基因在群體中分布較均勻。多態(tài)信息含量可作為一個(gè)衡量遺傳標(biāo)記所包含的或能提供的遺傳信息容量的指標(biāo)。一般認(rèn)為,在某一群體中,當(dāng)PIC> 0.5時(shí)該位點(diǎn)表現(xiàn)為高度多態(tài),0.25

    在許多海洋貝類中,經(jīng)常會(huì)出現(xiàn)偏分離的等位基因,這似乎是海洋貝類中普遍存在的現(xiàn)象。在本實(shí)驗(yàn)中,發(fā)現(xiàn)3 個(gè)微衛(wèi)星位點(diǎn)的基因型分離比顯著偏離孟德爾遺傳定律(P<0.01),部分位點(diǎn)偏分離的原因可能與生活環(huán)境和人工選擇有關(guān),再有可能是由隱性致死基因、取樣或基因型判讀錯(cuò)誤等因素引起[22]。大部分微衛(wèi)星位點(diǎn)的期望雜合度與觀測(cè)雜合度基本接近,說明群體內(nèi)的基因型頻率未發(fā)生大的改變。但有2 個(gè)位點(diǎn)SSR11和SSR164的觀測(cè)雜合度明顯大于期望雜合度,在SSR164發(fā)現(xiàn)了一種純合子完全缺失的現(xiàn)象(基因型AA為0),說明該位點(diǎn)附近可能存在著與之緊密連鎖的隱性致死基因[23]。

    菲律賓蛤仔的殼長(zhǎng)、殼寬、殼高、體重等性狀屬于數(shù)量性狀,數(shù)量性狀受多個(gè)基因的控制,其遺傳基礎(chǔ)比較復(fù)雜且容易受環(huán)境的影響。標(biāo)記與性狀之間的連鎖分析,是根據(jù)標(biāo)記位點(diǎn)的基因型以及數(shù)量性狀的表型對(duì)個(gè)體進(jìn)行顯著性檢驗(yàn),差異顯著則說明標(biāo)記與數(shù)量性狀存在關(guān)聯(lián)[24]。因此,如果一個(gè)群體的性狀差異顯著,或兩個(gè)群體的差異很大,就可以通過標(biāo)記與性狀的相關(guān)分析,找出性狀與一個(gè)或多個(gè)標(biāo)記的遺傳相關(guān),一旦發(fā)現(xiàn)顯著相關(guān),即可認(rèn)為存在一個(gè)數(shù)量性狀位點(diǎn),從而實(shí)現(xiàn)從表型到基因型選擇育種的轉(zhuǎn)變[25]。借助微衛(wèi)星標(biāo)記進(jìn)行性狀連鎖分析或 QTL 定位,在水產(chǎn)動(dòng)物育種中已也有較多報(bào)道[10-17]。

    目前,菲律賓蛤仔遺傳標(biāo)記的開發(fā)較少, 不足以滿足遺傳作圖和QTL 定位的需要。所以,有必要應(yīng)用有限的微衛(wèi)星標(biāo)記對(duì)菲律賓蛤仔生長(zhǎng)相關(guān)性狀進(jìn)行連鎖分析。在本研究中,共找到了4 個(gè)與生長(zhǎng)性狀相關(guān)的位點(diǎn),SSR9位點(diǎn)與殼高存在顯著的相關(guān)關(guān)系(P<0.05),SSR135和SSR164與殼寬顯著相關(guān)(P<0.05),SSR142與體重顯著相關(guān)(P<0.05)。在這些關(guān)聯(lián)中,出現(xiàn)了幾個(gè)標(biāo)記同一個(gè)性狀相關(guān),說明這些位點(diǎn)存在多因一效的現(xiàn)象,并符合數(shù)量性狀的相關(guān)理論。對(duì)上述位點(diǎn)不同基因型的多重比較表明,SSR9位點(diǎn)的AA基因型、SSR135位點(diǎn)的BB基因型、SSR142位點(diǎn)的AA基因型和SSR164位點(diǎn)的AB基因型都與殼長(zhǎng)、殼寬、殼高和體重性狀呈正相關(guān)。目前,分子標(biāo)記在菲律賓蛤仔中的應(yīng)用主要集中于遺傳多樣性研究[26-27],有關(guān)菲律賓蛤仔性狀相關(guān)分子標(biāo)記的篩選迄今未見報(bào)道。本文篩選出菲律賓蛤仔生長(zhǎng)性狀相關(guān)的微衛(wèi)星標(biāo)記,為菲律賓蛤仔分子標(biāo)記輔助選育提供參考。

    [1] FAO. Yearbook of fishery statistics summary tables. 2011, Available from ftp://ftp.fao.org/fi/stat/summary/default.htm#aqua.

    [2] 張國范, 閆喜武. 蛤仔養(yǎng)殖學(xué). 北京: 科學(xué)出版社, 2010.

    [3] 閆喜武, 張躍環(huán), 霍忠明, 孫煥強(qiáng), 楊鳳, 張國范. 不同殼色菲律賓蛤仔品系F2的表型性狀. 水產(chǎn)學(xué)報(bào), 2010, 34(6): 701-709.

    [4] 閆喜武, 張躍環(huán), 孫煥強(qiáng), 霍忠明, 孫欣, 楊鳳, 張國范. 菲律賓蛤仔兩道紅與白斑馬品系的三元雜交. 水產(chǎn)學(xué)報(bào), 2010, 34(8): 1190-1197.

    [5] Yan X W, Huo Z M, Yang F, Zhang G F. Heritability of larval and juvenile growth for two stocks of Manila clamRuditapesphilippinarum. Aquaculture Research, 2013, doi: 10.1111/j.1365-2109.2012.03250.x.

    [6] Zhao L Q, Yan X W, Huo Z M, Yang F, Zhang G F. Divergent selection for shell length in the Manila clam,Ruditapesphilippinarum. Journal of the World Aquaculture Society, 2012, 43(6): 878-884.

    [7] Cnaani A, Hallerman E M, Ron M, Weller J I, Indelman M, Kashi Y, Graham A E G, Hulata G. Detection of a chromosomal region with two quantitative trait loci, affecting cold tolerance and fish size, in an F2tilapia hybrid. Aquaculture, 2003, 223(1/4): 117-128.

    [8] Rodriguez M F, LaPatra S, Williams S, Famula T, May B. Genetic markers associated with resistance to infectious hematopoietic necrosis in rainbow and steelhead trout (Oncorhynchusmykiss) backcrosses. Aquaculture, 2004, 241(1/4): 93-115.

    [9] 高國強(qiáng), 常玉梅, 韓啟霞, 池炳杰, 李明云, 薛良義, 梁利群. 大黃魚耐低溫性狀相關(guān)微衛(wèi)星標(biāo)記的篩選. 遺傳, 2010, 32(3): 248-253.

    [10] 張?zhí)鞎r(shí), 劉萍, 李健, 孔杰, 王清印. 中國對(duì)蝦與生長(zhǎng)性狀相關(guān)微衛(wèi)星DNA分子標(biāo)記的初步研究. 海洋水產(chǎn)研究, 2006, 27(5): 34-38.

    [11] 樊佳佳, 白俊杰, 李小慧, 何小平, 何小燕, 李勝杰, 葉星, 吳立新. 大口黑鱸生長(zhǎng)性狀的微衛(wèi)星DNA標(biāo)記篩選. 遺傳, 2009, 31(5): 515-522.

    [12] 張研, 梁利群, 常玉梅, 侯寧, 魯翠云, 孫效文. 鯉魚體長(zhǎng)性狀的QTL定位及其遺傳效應(yīng)分析. 遺傳, 2007, 29(10): 1243-1248.

    [13] 王桂興, 劉永新, 孫效文, 劉海金. 牙鲆微衛(wèi)星分子標(biāo)記與生長(zhǎng)性狀的相關(guān)性分析. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2009, 40(7): 77-84.

    [14] 劉賢德, 韋信鍵, 蔡明夷, 劉洋, 王志勇. 大黃魚22個(gè)微衛(wèi)星標(biāo)記在F1家系中的分離方式及與生長(zhǎng)性狀的相關(guān)分析. 水產(chǎn)學(xué)報(bào), 2012, 36(9): 1322-1330.

    [15] 劉磊, 李健, 劉萍, 趙法箴, 高保全, 杜盈, 馬春艷. 三疣梭子蟹微衛(wèi)星標(biāo)記與生長(zhǎng)相關(guān)性狀的相關(guān)性分析. 水產(chǎn)學(xué)報(bào),2012, 36(7): 1034-1041.

    [16] 張曼, 陳蒙, 姬南京. 蝦夷扇貝家系群體遺傳結(jié)構(gòu)及其微衛(wèi)星標(biāo)記與經(jīng)濟(jì)性狀相關(guān)性的分析研究. 中國農(nóng)學(xué)通報(bào), 2012, 28(20): 125-130.

    [17] 鄧岳文, 高遠(yuǎn)鎮(zhèn), 王學(xué)穎, 杜曉東. 馬氏珠母貝生長(zhǎng)性狀與EST-SSR標(biāo)記的關(guān)聯(lián)分析. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2013, 21(1): 77-88.

    [18] Botstein D, White R L, Skolnick M, Davis R W. Construction of genetic linkage map in man using restriction fragment length polymorphisms. American Journal of Human Genetics, 1980, 32(3): 314-331.

    [19] Yang M S, Hertel H, Schneck V. Genetic diversity and population structure ofRobiniapseudoacaciaprovenances from middle Europe. Acta Ecologica Sinica, 2004, 24(12): 2700-2706.

    [20] Hines H C, ZiKakis J P, Haenlein G F, Kiddy C G, Trowbridge C L. Linkage relationships among loci of polymorphism in blood and milk of cattle. Journal of Dairy Science, l 981, 64(1): 71-76.

    [21] Wright S. Evolution and the Genetics of Populations: Variability Within and Among Natural Population. Chicago: University of Chicago Press, 1978.

    [22] Launey S, Hedgecock D. High genetic load in the Pacific oysterCrassostreagigas. Genetics, 2001, 159(1): 255-265.

    [23] Nie H T, Li Q, Kong L. Centromere mapping in the Pacific abalone (Haliotisdiscushannai) through half-tetrad analysis in gynogenetic diploid families. Animal Genetics, 2012, 43(3): 290-297.

    [25] 鐘金城, 陳智華. 分子遺傳學(xué)與動(dòng)物育種. 成都: 四川大學(xué)出版社, 2001.

    [26] 虞志飛, 閆喜武, 楊霏, 王金海, 張躍環(huán), 楊鳳, 張國范. 菲律賓蛤仔大連群體不同世代的遺傳多樣性. 生態(tài)學(xué)報(bào), 2011, 31(15): 4199-4206.

    [27] 閆喜武, 虞志飛, 秦艷杰, 楊霏, 王金海, 張躍環(huán), 楊鳳, 張國范. 菲律賓蛤仔EST-SSR標(biāo)記開發(fā)及不同地理群體遺傳多樣性. 生態(tài)學(xué)報(bào), 2011, 31(15): 4190-4198.

    Identification of EST_SSR markers associated with growth-related traits in the Manila clamRuditapesphilippinarum

    NIU Hongbo, NIE Hongtao, ZHU Depeng, YANG Feng, YAN Xiwu*

    EngineeringandTechnologyResearchCenterofShellfishBreedingofLiaoningProvince,DalianOceanUniversity,Dalian116023,China

    The Manila clam,Ruditapesphilippinarum, which is widely distributed along the coasts of China, is an economically important marine bivalve species in China′s aquaculture industry. The world production of this species was 3.6 million metric tons in 2010. China is the first largest country in the world in terms of production of the Manila clam, producing about 3.0 million metric tons annually, which accounts for about 90% of global production. This species has several pedigrees including White, Zebra, Liangdao Red and Marine Red distributing in the coastal areas in North China. Microsatellites or simple sequence repeats (SSRs) are tandemly repeated motifs of 1—6 genetic base pairs. Microsatellite markers are powerful molecular markers due to their high polymorphism, stability, and co-dominance, and they are used widely in studies of genetic diversity, parentage assignment, genetic linkage map construction, and trait-related marker screening. In this study, 20 microsatellite DNA markers were used to analyze the genetic diversity of 107 individuals of a Zebra F2 pedigree ofR.philippinarum. Forty-one alleles were detected, and the number of alleles (Na) was 2—3 at each locus (average, 2.05). The effective number of alleles (Ne) was 1.71, and the DNA fragment length was 109—430 base pairs. The mean values of observed heterozygosity (Ho), expected heterozygosity (He), and polymorphism information content (PIC) were 0.504, 0.431, and 0.324, respectively. The probability value of the chi-square test showed that three loci significantly deviated from Mendelian segregation (P<0.01), which suggested that these loci might link with the adaptive gene, and two loci (SSR11 and SSR164) may link with recessive homozygous lethal genes. The general linear model procedure in SPSS20.0 was used to analyze the correlation between the 20 microsatellites and growth-related traits ofR.philippinarum(i.e., shell length, shell width, shell height, and body weight). Four loci were significantly related to the growth traits (P<0.05): SSR135 and SSR164 were significantly related to shell width (P<0.05), SSR9 was significantly related to shell height (P<0.05), and SSR142 was significantly related to total weight (P<0.05). Favorable genotypes for each growth trait were identified by a multiple comparison among the loci. Allele A for SSR9 had a significant impact on shell height (P<0.05) and had the highest phenotype value; thus, it can be used as a molecular marker for selective breeding. SSR135 and SSR164 had a significant impact on shell width and total weight (P<0.05), respectively, and allele B for SSR135 and allele A for SSR142 had a positive effect on growth-related traits. The following four genotypes of these loci had a favorable effect on growth-related traits: AA for SSR9, BB for SSR135, AA for SSR142, and AB for SSR164. The trait-related microsatellite loci identified in this study will be valuable for marker-assisted breeding ofR.philippinarum.

    Ruditapesphilippinarum; microsatellites; growth-related traits; correlation analysis

    國家自然科學(xué)基金(31302183); 國家高技術(shù)研究發(fā)展計(jì)劃“863”項(xiàng)目(2012AA10A400); 現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-48)

    2013-05-16;

    日期:2014-04-25

    10.5846/stxb201305161083

    *通訊作者Corresponding author.E-mail: yanxiwu@dlou.edu

    牛泓博,聶鴻濤,朱德鵬,楊鳳,閆喜武.菲律賓蛤仔EST_SSR標(biāo)記與生長(zhǎng)性狀的相關(guān)分析.生態(tài)學(xué)報(bào),2015,35(6):1910-1916.

    Niu H B, Nie H T, Zhu D P, Yang F, Yan X W.Identification of EST_SSR markers associated with growth-related traits in the Manila clamRuditapesphilippinarum.Acta Ecologica Sinica,2015,35(6):1910-1916.

    猜你喜歡
    蛤仔微衛(wèi)星雜合
    蛤仔山東萊州群體與福建漳浦群體生長(zhǎng)存活比較研究
    梅花鹿基因組微衛(wèi)星分布特征研究
    遼寧蛤仔產(chǎn)業(yè)發(fā)展初探
    甘藍(lán)型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國“一箭雙星”成功將“遙感衛(wèi)星二十一號(hào)”與“天拓二號(hào)視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    不同海藻餌料對(duì)菲律賓蛤仔生長(zhǎng)的影響
    雄激素可調(diào)節(jié)的腎臟近端腎小管上皮細(xì)胞靶向雜合啟動(dòng)子的優(yōu)化
    成人鲁丝片一二三区免费| 国产久久久一区二区三区| 亚洲国产欧美人成| 国产主播在线观看一区二区| 亚洲 国产 在线| 亚洲成人免费电影在线观看| 一进一出抽搐gif免费好疼| 在线a可以看的网站| 首页视频小说图片口味搜索| 精品久久久久久成人av| 色尼玛亚洲综合影院| 国产精品久久久久久久久免 | 韩国av一区二区三区四区| 久久久精品大字幕| 99精品在免费线老司机午夜| 亚洲在线观看片| 国产精品野战在线观看| 香蕉久久夜色| 国产伦一二天堂av在线观看| 黄色成人免费大全| 老司机在亚洲福利影院| 一级黄色大片毛片| 女同久久另类99精品国产91| 午夜激情福利司机影院| 午夜福利成人在线免费观看| 国产精品久久电影中文字幕| 真实男女啪啪啪动态图| 熟妇人妻久久中文字幕3abv| 欧美+日韩+精品| 欧美最黄视频在线播放免费| 两个人看的免费小视频| 人妻丰满熟妇av一区二区三区| 一级毛片高清免费大全| 欧美日韩亚洲国产一区二区在线观看| 国产单亲对白刺激| 欧美国产日韩亚洲一区| 尤物成人国产欧美一区二区三区| 少妇的逼水好多| 国产欧美日韩精品一区二区| 欧美一级a爱片免费观看看| 丁香六月欧美| 一本综合久久免费| 精品一区二区三区视频在线观看免费| 国产精品香港三级国产av潘金莲| 最近视频中文字幕2019在线8| 美女cb高潮喷水在线观看| 日本黄色片子视频| 真实男女啪啪啪动态图| 国产探花极品一区二区| 国产乱人伦免费视频| 国产精品久久久久久久电影 | 国内揄拍国产精品人妻在线| 亚洲精品色激情综合| 国产av不卡久久| 国产精品女同一区二区软件 | 日韩欧美国产一区二区入口| 午夜免费观看网址| 搡老妇女老女人老熟妇| 最后的刺客免费高清国语| 午夜激情欧美在线| 国内精品久久久久精免费| 香蕉久久夜色| 中文字幕av成人在线电影| 黄色女人牲交| 久久久久久大精品| 女人十人毛片免费观看3o分钟| 成人特级av手机在线观看| 亚洲最大成人手机在线| 成人无遮挡网站| 午夜免费激情av| 麻豆国产av国片精品| 免费在线观看成人毛片| 欧美极品一区二区三区四区| 欧美大码av| 亚洲成av人片在线播放无| 国产一区二区三区视频了| 1000部很黄的大片| 在线观看av片永久免费下载| 12—13女人毛片做爰片一| 亚洲av成人不卡在线观看播放网| 少妇人妻一区二区三区视频| 超碰av人人做人人爽久久 | 三级毛片av免费| 成年免费大片在线观看| 99久久综合精品五月天人人| 一级毛片高清免费大全| 香蕉久久夜色| 1000部很黄的大片| 亚洲精品成人久久久久久| 亚洲熟妇中文字幕五十中出| 女人十人毛片免费观看3o分钟| 成年女人毛片免费观看观看9| 国产高清videossex| 18+在线观看网站| 亚洲国产精品合色在线| 麻豆久久精品国产亚洲av| 亚洲成人精品中文字幕电影| 天天添夜夜摸| 久久久久九九精品影院| 女警被强在线播放| 蜜桃亚洲精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 99热精品在线国产| 国产精品国产高清国产av| 成人特级黄色片久久久久久久| 精品一区二区三区视频在线观看免费| 精品久久久久久,| 亚洲性夜色夜夜综合| 欧美xxxx黑人xx丫x性爽| 午夜激情欧美在线| 99久久精品国产亚洲精品| 中文字幕av在线有码专区| 老熟妇乱子伦视频在线观看| 国产精品亚洲一级av第二区| 精品不卡国产一区二区三区| 两个人的视频大全免费| 少妇丰满av| 亚洲欧美日韩卡通动漫| 给我免费播放毛片高清在线观看| 最新在线观看一区二区三区| 天堂影院成人在线观看| 国产精品亚洲一级av第二区| 亚洲精品一区av在线观看| 在线免费观看的www视频| 观看美女的网站| 18禁裸乳无遮挡免费网站照片| 精品午夜福利视频在线观看一区| 日韩大尺度精品在线看网址| 日韩国内少妇激情av| 国产精品一区二区三区四区免费观看 | 中出人妻视频一区二区| 蜜桃亚洲精品一区二区三区| 国产精品乱码一区二三区的特点| 麻豆成人av在线观看| 村上凉子中文字幕在线| 国产私拍福利视频在线观看| 久久久久久国产a免费观看| 真人一进一出gif抽搐免费| 国产探花在线观看一区二区| 亚洲国产高清在线一区二区三| 不卡一级毛片| 男人舔女人下体高潮全视频| 欧美性猛交╳xxx乱大交人| 一级作爱视频免费观看| 99热这里只有是精品50| 国产激情欧美一区二区| 国产老妇女一区| 国产野战对白在线观看| 老汉色∧v一级毛片| 村上凉子中文字幕在线| 九色成人免费人妻av| 99国产综合亚洲精品| 99精品久久久久人妻精品| 高清在线国产一区| 欧美高清成人免费视频www| 国产视频一区二区在线看| 人人妻,人人澡人人爽秒播| 99久久精品热视频| 国产成年人精品一区二区| 叶爱在线成人免费视频播放| 无人区码免费观看不卡| 日本免费a在线| 禁无遮挡网站| 两人在一起打扑克的视频| 尤物成人国产欧美一区二区三区| 看黄色毛片网站| h日本视频在线播放| 久久欧美精品欧美久久欧美| 男女下面进入的视频免费午夜| 十八禁网站免费在线| 免费搜索国产男女视频| 国内精品久久久久久久电影| www国产在线视频色| 蜜桃久久精品国产亚洲av| 在线国产一区二区在线| 此物有八面人人有两片| 欧美日韩乱码在线| 免费在线观看成人毛片| 怎么达到女性高潮| 欧美区成人在线视频| 人人妻人人澡欧美一区二区| 成人特级av手机在线观看| 亚洲黑人精品在线| 99视频精品全部免费 在线| 日韩 欧美 亚洲 中文字幕| 午夜福利免费观看在线| 好男人在线观看高清免费视频| 国产麻豆成人av免费视频| 中文字幕久久专区| 变态另类成人亚洲欧美熟女| 国产免费av片在线观看野外av| 天天躁日日操中文字幕| 老司机午夜十八禁免费视频| 最近最新中文字幕大全电影3| 日韩欧美在线乱码| 日韩免费av在线播放| 久久久久久久久大av| 51午夜福利影视在线观看| 在线观看免费午夜福利视频| 亚洲av成人精品一区久久| 我的老师免费观看完整版| 欧美在线黄色| 午夜影院日韩av| 亚洲av电影在线进入| 久久精品91无色码中文字幕| tocl精华| 国产激情偷乱视频一区二区| 久久精品人妻少妇| 免费电影在线观看免费观看| 啪啪无遮挡十八禁网站| 国产单亲对白刺激| 国产精品电影一区二区三区| 国产精品久久久久久久电影 | 亚洲精品影视一区二区三区av| 久久久国产成人精品二区| 人人妻人人看人人澡| 在线观看av片永久免费下载| 欧美日韩黄片免| 午夜影院日韩av| 女生性感内裤真人,穿戴方法视频| 亚洲黑人精品在线| 99热6这里只有精品| 最近最新中文字幕大全电影3| 国产精品1区2区在线观看.| 亚洲 欧美 日韩 在线 免费| 欧美在线黄色| 免费搜索国产男女视频| 午夜久久久久精精品| 日日夜夜操网爽| 久久精品91蜜桃| 91麻豆精品激情在线观看国产| 欧美xxxx黑人xx丫x性爽| 操出白浆在线播放| 国内精品久久久久精免费| 人妻夜夜爽99麻豆av| 午夜免费成人在线视频| 久久婷婷人人爽人人干人人爱| 久久久久久久亚洲中文字幕 | 午夜福利成人在线免费观看| 又黄又爽又免费观看的视频| 欧美av亚洲av综合av国产av| 国产黄a三级三级三级人| 夜夜爽天天搞| 内地一区二区视频在线| 国产毛片a区久久久久| 午夜福利欧美成人| 欧美性猛交╳xxx乱大交人| 午夜激情欧美在线| 美女被艹到高潮喷水动态| 欧美乱码精品一区二区三区| 中文亚洲av片在线观看爽| 日韩欧美 国产精品| 亚洲18禁久久av| 一进一出抽搐动态| 99热精品在线国产| 久久欧美精品欧美久久欧美| 中文字幕久久专区| 淫妇啪啪啪对白视频| 麻豆久久精品国产亚洲av| 亚洲熟妇中文字幕五十中出| 一边摸一边抽搐一进一小说| 国产 一区 欧美 日韩| 色老头精品视频在线观看| 国产主播在线观看一区二区| av女优亚洲男人天堂| 日本熟妇午夜| 男女做爰动态图高潮gif福利片| 老鸭窝网址在线观看| 色综合欧美亚洲国产小说| 国产亚洲精品久久久久久毛片| 亚洲av免费在线观看| 免费在线观看成人毛片| 黄色日韩在线| 99热只有精品国产| 国内精品久久久久久久电影| 又黄又爽又免费观看的视频| 欧美日本亚洲视频在线播放| 九九热线精品视视频播放| 国产精品自产拍在线观看55亚洲| 18禁裸乳无遮挡免费网站照片| 一级毛片女人18水好多| 亚洲真实伦在线观看| 国产欧美日韩精品亚洲av| 国产精品免费一区二区三区在线| 叶爱在线成人免费视频播放| 亚洲欧美激情综合另类| 国产伦精品一区二区三区四那| 国产精品久久久久久亚洲av鲁大| 亚洲精品粉嫩美女一区| 免费在线观看成人毛片| 国产伦一二天堂av在线观看| 色老头精品视频在线观看| 国产蜜桃级精品一区二区三区| 午夜免费激情av| av专区在线播放| 国产视频一区二区在线看| av国产免费在线观看| 欧美性猛交╳xxx乱大交人| 色视频www国产| h日本视频在线播放| 波多野结衣高清作品| 一本一本综合久久| av黄色大香蕉| 免费看日本二区| 美女被艹到高潮喷水动态| 三级国产精品欧美在线观看| 久久这里只有精品中国| 99久久无色码亚洲精品果冻| aaaaa片日本免费| 色老头精品视频在线观看| 18禁国产床啪视频网站| 久久久国产成人免费| 欧美乱码精品一区二区三区| 男人舔女人下体高潮全视频| 观看美女的网站| 色视频www国产| 一级毛片女人18水好多| 免费无遮挡裸体视频| 欧美黑人欧美精品刺激| 老司机福利观看| 又紧又爽又黄一区二区| a级毛片a级免费在线| 免费观看精品视频网站| 国产精品久久久久久精品电影| 国产亚洲av嫩草精品影院| 小蜜桃在线观看免费完整版高清| 国产高清激情床上av| 国产亚洲av嫩草精品影院| 少妇的逼水好多| 啦啦啦免费观看视频1| 午夜免费激情av| 小蜜桃在线观看免费完整版高清| 久久亚洲精品不卡| 日韩有码中文字幕| 亚洲av五月六月丁香网| 19禁男女啪啪无遮挡网站| 国产亚洲欧美98| 国语自产精品视频在线第100页| 国产蜜桃级精品一区二区三区| 国产淫片久久久久久久久 | 久久久国产精品麻豆| 午夜福利在线观看免费完整高清在 | 一个人观看的视频www高清免费观看| 精品久久久久久成人av| 国产亚洲精品久久久久久毛片| 综合色av麻豆| 精品久久久久久久久久久久久| 免费人成在线观看视频色| 婷婷丁香在线五月| 国产aⅴ精品一区二区三区波| 久久香蕉精品热| 久久久久久国产a免费观看| 欧美zozozo另类| 久久久久国内视频| 一级黄片播放器| 国产老妇女一区| 亚洲av熟女| 国产激情偷乱视频一区二区| 99热只有精品国产| 亚洲第一电影网av| www日本黄色视频网| xxxwww97欧美| 好男人电影高清在线观看| 日本免费a在线| 91在线观看av| 欧美一区二区国产精品久久精品| 99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频| 国产精品 欧美亚洲| 女人高潮潮喷娇喘18禁视频| 亚洲av第一区精品v没综合| 色综合婷婷激情| 亚洲熟妇中文字幕五十中出| 99久久成人亚洲精品观看| 人人妻,人人澡人人爽秒播| 嫁个100分男人电影在线观看| 久久精品亚洲精品国产色婷小说| 亚洲av美国av| 99久久综合精品五月天人人| 在线观看av片永久免费下载| 日韩欧美精品免费久久 | 9191精品国产免费久久| 午夜免费观看网址| 成人鲁丝片一二三区免费| 麻豆久久精品国产亚洲av| 丝袜美腿在线中文| 国产av麻豆久久久久久久| 一区福利在线观看| 国产欧美日韩一区二区精品| 国产伦一二天堂av在线观看| 亚洲片人在线观看| 欧美成人一区二区免费高清观看| 国产精品野战在线观看| 国产高清视频在线播放一区| 精品免费久久久久久久清纯| 熟女少妇亚洲综合色aaa.| 亚洲人成网站高清观看| 看免费av毛片| 婷婷精品国产亚洲av在线| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 美女高潮喷水抽搐中文字幕| 欧美色欧美亚洲另类二区| 熟女电影av网| 成人午夜高清在线视频| 久久久久久久久中文| 久久婷婷人人爽人人干人人爱| 亚洲av熟女| 国产欧美日韩一区二区精品| 国产免费男女视频| 波多野结衣高清作品| 久99久视频精品免费| 日韩欧美在线二视频| 亚洲精品在线观看二区| 午夜精品在线福利| aaaaa片日本免费| 亚洲欧美日韩高清在线视频| 国产成人系列免费观看| 老司机午夜福利在线观看视频| av天堂在线播放| 51国产日韩欧美| 身体一侧抽搐| 国产精品嫩草影院av在线观看 | 色综合站精品国产| 国产精品久久电影中文字幕| 国产精品乱码一区二三区的特点| 亚洲色图av天堂| 长腿黑丝高跟| 亚洲男人的天堂狠狠| 97碰自拍视频| a在线观看视频网站| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品999在线| 18禁裸乳无遮挡免费网站照片| 高清毛片免费观看视频网站| 午夜老司机福利剧场| 最近最新中文字幕大全电影3| 99久久久亚洲精品蜜臀av| 国产亚洲av嫩草精品影院| 中文亚洲av片在线观看爽| 午夜福利在线观看吧| 少妇高潮的动态图| 狠狠狠狠99中文字幕| 免费在线观看日本一区| 岛国视频午夜一区免费看| 亚洲天堂国产精品一区在线| 18禁国产床啪视频网站| 尤物成人国产欧美一区二区三区| 欧美bdsm另类| 老司机福利观看| 村上凉子中文字幕在线| 99在线视频只有这里精品首页| 日韩中文字幕欧美一区二区| 18美女黄网站色大片免费观看| 免费无遮挡裸体视频| 啪啪无遮挡十八禁网站| 在线十欧美十亚洲十日本专区| 中文字幕熟女人妻在线| 窝窝影院91人妻| 波多野结衣高清作品| 中文字幕人妻熟人妻熟丝袜美 | 亚洲精品影视一区二区三区av| 国产精品99久久久久久久久| 99国产综合亚洲精品| 久久精品国产综合久久久| 国产精品三级大全| 怎么达到女性高潮| 变态另类丝袜制服| 一夜夜www| 最近最新中文字幕大全电影3| 国产高清激情床上av| 亚洲精品乱码久久久v下载方式 | 亚洲欧美一区二区三区黑人| 在线播放无遮挡| 久久久国产精品麻豆| 美女被艹到高潮喷水动态| 十八禁网站免费在线| 精品久久久久久久久久久久久| 男女下面进入的视频免费午夜| 亚洲欧美日韩卡通动漫| 亚洲国产日韩欧美精品在线观看 | 国内精品久久久久精免费| 久久精品国产亚洲av香蕉五月| 少妇的丰满在线观看| 最好的美女福利视频网| 18禁黄网站禁片午夜丰满| 精品一区二区三区人妻视频| 亚洲五月婷婷丁香| 日韩中文字幕欧美一区二区| 国产伦人伦偷精品视频| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 午夜免费男女啪啪视频观看 | 成人av在线播放网站| 99久久久亚洲精品蜜臀av| 一本一本综合久久| 大型黄色视频在线免费观看| 午夜福利成人在线免费观看| 久久精品影院6| 内射极品少妇av片p| 成人无遮挡网站| 在线观看av片永久免费下载| 熟女少妇亚洲综合色aaa.| 亚洲国产欧洲综合997久久,| 老司机深夜福利视频在线观看| 18美女黄网站色大片免费观看| 日韩人妻高清精品专区| 欧美一级毛片孕妇| 一二三四社区在线视频社区8| 一本综合久久免费| 男插女下体视频免费在线播放| av天堂在线播放| 高清在线国产一区| a级毛片a级免费在线| 欧美日韩精品网址| 在线观看免费午夜福利视频| avwww免费| 脱女人内裤的视频| 国产高清有码在线观看视频| 亚洲成人久久性| 午夜激情欧美在线| 亚洲自拍偷在线| 99久久综合精品五月天人人| 搡女人真爽免费视频火全软件 | 美女被艹到高潮喷水动态| 乱人视频在线观看| 久久精品国产亚洲av涩爱 | 日韩欧美一区二区三区在线观看| 女人被狂操c到高潮| 亚洲在线观看片| 18禁国产床啪视频网站| 淫妇啪啪啪对白视频| 一个人免费在线观看的高清视频| 少妇的逼水好多| 欧美成人性av电影在线观看| 欧美区成人在线视频| 亚洲国产色片| 人人妻,人人澡人人爽秒播| 久久草成人影院| 国产精品香港三级国产av潘金莲| 精品人妻一区二区三区麻豆 | 亚洲,欧美精品.| 波多野结衣巨乳人妻| 国产伦人伦偷精品视频| 欧美绝顶高潮抽搐喷水| 亚洲精品影视一区二区三区av| 亚洲精品在线美女| 一进一出抽搐动态| 日本五十路高清| 成人亚洲精品av一区二区| 国产伦精品一区二区三区视频9 | 欧美日韩综合久久久久久 | 久久欧美精品欧美久久欧美| 久久草成人影院| 国产免费av片在线观看野外av| 天堂网av新在线| 少妇的逼好多水| 亚洲国产欧美网| 免费大片18禁| 老司机在亚洲福利影院| 久久性视频一级片| 亚洲精品在线观看二区| 蜜桃亚洲精品一区二区三区| 在线观看一区二区三区| 男女那种视频在线观看| 国产黄a三级三级三级人| 亚洲午夜理论影院| 色播亚洲综合网| 亚洲中文字幕一区二区三区有码在线看| 热99re8久久精品国产| 搡老岳熟女国产| 在线播放国产精品三级| 欧美黄色片欧美黄色片| 免费在线观看日本一区| 亚洲av电影在线进入| 欧美精品啪啪一区二区三区| 中亚洲国语对白在线视频| 亚洲精品一卡2卡三卡4卡5卡| 1000部很黄的大片| 真人一进一出gif抽搐免费| 一进一出抽搐gif免费好疼| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 亚洲精品乱码久久久v下载方式 | xxxwww97欧美| 色噜噜av男人的天堂激情| 最近最新中文字幕大全电影3| 国产一区二区在线av高清观看| netflix在线观看网站| 老熟妇仑乱视频hdxx| 日本撒尿小便嘘嘘汇集6| 18禁国产床啪视频网站| 少妇熟女aⅴ在线视频| 三级毛片av免费| 亚洲,欧美精品.| 又黄又爽又免费观看的视频| 一级黄片播放器| 国产成年人精品一区二区| 成年版毛片免费区| 非洲黑人性xxxx精品又粗又长| 美女高潮喷水抽搐中文字幕| 少妇高潮的动态图| 国产视频内射| 天堂√8在线中文| 变态另类成人亚洲欧美熟女| 国产单亲对白刺激| 久久精品综合一区二区三区| 久久精品91蜜桃| 亚洲国产中文字幕在线视频| 久久久久久久午夜电影| 动漫黄色视频在线观看| 精品午夜福利视频在线观看一区| a级毛片a级免费在线| 日本黄大片高清|