• <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)化
    波野结衣二区三区在线| 久久久色成人| 免费人成在线观看视频色| 男女国产视频网站| 纵有疾风起免费观看全集完整版 | 麻豆久久精品国产亚洲av| 可以在线观看毛片的网站| 韩国av在线不卡| 久久韩国三级中文字幕| 亚洲国产欧美在线一区| 在线观看一区二区三区| 国产精品国产三级专区第一集| 国产av码专区亚洲av| 久久久久久久亚洲中文字幕| 亚洲成人久久爱视频| 一边亲一边摸免费视频| 亚洲精品自拍成人| 日本三级黄在线观看| av国产久精品久网站免费入址| 一级片'在线观看视频| 永久免费av网站大全| 国产免费又黄又爽又色| 99热全是精品| 亚洲aⅴ乱码一区二区在线播放| 久久久精品94久久精品| 婷婷色麻豆天堂久久| 国产淫语在线视频| 久99久视频精品免费| 亚洲18禁久久av| 午夜日本视频在线| 中文字幕av在线有码专区| 亚洲欧美中文字幕日韩二区| 日韩伦理黄色片| 熟女电影av网| av在线亚洲专区| av专区在线播放| 国产色婷婷99| 国产亚洲一区二区精品| 国产精品av视频在线免费观看| 国产成人精品久久久久久| 久久久精品免费免费高清| 午夜老司机福利剧场| 一级片'在线观看视频| 九草在线视频观看| 国产精品国产三级国产专区5o| 国产精品一区二区三区四区久久| 身体一侧抽搐| 91精品国产九色| 丰满少妇做爰视频| 日本猛色少妇xxxxx猛交久久| 日韩欧美 国产精品| 韩国av在线不卡| 哪个播放器可以免费观看大片| 婷婷色av中文字幕| 国产精品日韩av在线免费观看| 乱码一卡2卡4卡精品| 免费av毛片视频| 国产精品国产三级专区第一集| 欧美不卡视频在线免费观看| 国产成人a区在线观看| 午夜激情欧美在线| 午夜福利视频1000在线观看| 97在线视频观看| 亚洲国产精品sss在线观看| 九草在线视频观看| 成年女人在线观看亚洲视频 | 色哟哟·www| 亚洲欧美精品专区久久| 国产成人freesex在线| 日韩国内少妇激情av| 蜜桃久久精品国产亚洲av| 亚洲国产最新在线播放| 国国产精品蜜臀av免费| 美女大奶头视频| 亚洲av日韩在线播放| 99re6热这里在线精品视频| av在线老鸭窝| 啦啦啦啦在线视频资源| 男人爽女人下面视频在线观看| 国产一区二区亚洲精品在线观看| 中文欧美无线码| 在线a可以看的网站| 中国美白少妇内射xxxbb| 极品少妇高潮喷水抽搐| 少妇人妻精品综合一区二区| 久久97久久精品| 最近视频中文字幕2019在线8| 精华霜和精华液先用哪个| 精品酒店卫生间| 一级毛片aaaaaa免费看小| 欧美日韩亚洲高清精品| 老司机影院成人| 日本av手机在线免费观看| 一边亲一边摸免费视频| 91av网一区二区| 国产永久视频网站| 日韩在线高清观看一区二区三区| 亚洲av电影不卡..在线观看| 亚洲国产精品成人久久小说| 久久97久久精品| 青青草视频在线视频观看| 男人舔女人下体高潮全视频| 一区二区三区高清视频在线| 黄色配什么色好看| 大又大粗又爽又黄少妇毛片口| 久久精品国产亚洲av天美| 最近的中文字幕免费完整| 色播亚洲综合网| 久久久精品94久久精品| 精品国产三级普通话版| av播播在线观看一区| 国内精品美女久久久久久| 女人久久www免费人成看片| 国产在线男女| 国产在视频线在精品| 亚洲一区高清亚洲精品| 久久久久性生活片| 偷拍熟女少妇极品色| 国产一区亚洲一区在线观看| 精品久久久久久电影网| 深夜a级毛片| 国产女主播在线喷水免费视频网站 | 2018国产大陆天天弄谢| 免费av毛片视频| 深夜a级毛片| 日本与韩国留学比较| 亚洲欧美精品自产自拍| 你懂的网址亚洲精品在线观看| 久久久久国产网址| 亚洲人成网站在线观看播放| 又黄又爽又刺激的免费视频.| 蜜桃久久精品国产亚洲av| 免费看日本二区| 欧美成人午夜免费资源| 韩国av在线不卡| 简卡轻食公司| 日韩亚洲欧美综合| 欧美成人a在线观看| 听说在线观看完整版免费高清| 久久久国产一区二区| 国内揄拍国产精品人妻在线| 18禁在线无遮挡免费观看视频| 国产精品一区二区三区四区久久| 26uuu在线亚洲综合色| 深夜a级毛片| 国产永久视频网站| 一个人看的www免费观看视频| 色5月婷婷丁香| 免费少妇av软件| 极品教师在线视频| 亚洲欧美一区二区三区国产| 午夜福利视频1000在线观看| 女人十人毛片免费观看3o分钟| 国产精品久久久久久精品电影小说 | 日本与韩国留学比较| 久久99精品国语久久久| a级毛片免费高清观看在线播放| 国产老妇伦熟女老妇高清| 国产免费又黄又爽又色| 成人毛片60女人毛片免费| 色吧在线观看| 一本久久精品| 久久久久久久国产电影| 久热久热在线精品观看| 蜜桃久久精品国产亚洲av| 亚洲va在线va天堂va国产| 亚洲人成网站在线播| 欧美不卡视频在线免费观看| 久久精品人妻少妇| 2018国产大陆天天弄谢| 丰满少妇做爰视频| 在线观看人妻少妇| 亚洲电影在线观看av| av免费观看日本| 男女国产视频网站| 在现免费观看毛片| 久久精品国产亚洲av天美| 黄色配什么色好看| 日韩精品有码人妻一区| 97热精品久久久久久| 久久久久久久久久久丰满| 国产成人freesex在线| 男女下面进入的视频免费午夜| 国产精品久久久久久av不卡| 女的被弄到高潮叫床怎么办| 日日摸夜夜添夜夜添av毛片| 国产av码专区亚洲av| 国产精品国产三级专区第一集| 亚洲精品视频女| 能在线免费观看的黄片| 国产精品久久视频播放| 日韩亚洲欧美综合| 简卡轻食公司| 免费黄色在线免费观看| 99热这里只有是精品在线观看| 国产精品伦人一区二区| 女的被弄到高潮叫床怎么办| 日本三级黄在线观看| 中文字幕人妻熟人妻熟丝袜美| 日日啪夜夜爽| 最近最新中文字幕大全电影3| 成人一区二区视频在线观看| 美女被艹到高潮喷水动态| 在线播放无遮挡| 一个人看视频在线观看www免费| 少妇熟女aⅴ在线视频| 夫妻性生交免费视频一级片| 欧美区成人在线视频| 亚洲精品自拍成人| av播播在线观看一区| 久久国产乱子免费精品| 国产伦精品一区二区三区视频9| 日韩制服骚丝袜av| 男女边吃奶边做爰视频| 欧美日韩视频高清一区二区三区二| 亚洲高清免费不卡视频| 人妻少妇偷人精品九色| 天堂影院成人在线观看| 男人狂女人下面高潮的视频| 久久精品久久久久久久性| 亚洲av不卡在线观看| 99久久人妻综合| 国产日韩欧美在线精品| 尾随美女入室| 99久久九九国产精品国产免费| 精品不卡国产一区二区三区| 亚洲经典国产精华液单| 寂寞人妻少妇视频99o| 2022亚洲国产成人精品| 日本色播在线视频| 99久久中文字幕三级久久日本| 日本免费a在线| 能在线免费看毛片的网站| 国产黄色视频一区二区在线观看| 国产69精品久久久久777片| av线在线观看网站| 99热这里只有是精品50| 联通29元200g的流量卡| 免费大片黄手机在线观看| 亚洲av二区三区四区| 成人午夜精彩视频在线观看| 国产单亲对白刺激| 亚洲va在线va天堂va国产| a级毛片免费高清观看在线播放| 最新中文字幕久久久久| 国产成人精品久久久久久| 日韩亚洲欧美综合| 亚洲av.av天堂| 国产一级毛片七仙女欲春2| 国产激情偷乱视频一区二区| 久久99热6这里只有精品| 嫩草影院入口| 国产黄色小视频在线观看| 99热网站在线观看| 亚洲18禁久久av| av在线亚洲专区| 久久久国产一区二区| 午夜激情久久久久久久| 久久精品久久精品一区二区三区| 久久久精品欧美日韩精品| 欧美zozozo另类| 午夜福利视频1000在线观看| 国产麻豆成人av免费视频| 91精品国产九色| 国产亚洲av片在线观看秒播厂 | 能在线免费观看的黄片| 日本黄色片子视频| 国产亚洲av片在线观看秒播厂 | 一个人看视频在线观看www免费| 最近中文字幕2019免费版| 一级毛片黄色毛片免费观看视频| ponron亚洲| 国产精品福利在线免费观看| 九九久久精品国产亚洲av麻豆| 免费不卡的大黄色大毛片视频在线观看 | 99久久精品热视频| 一个人观看的视频www高清免费观看| 69av精品久久久久久| 亚洲最大成人手机在线| av国产久精品久网站免费入址| 丰满少妇做爰视频| 欧美xxxx性猛交bbbb| 国产一级毛片在线| 最近最新中文字幕大全电影3| 97超碰精品成人国产| 国产精品一区二区三区四区免费观看| 在线a可以看的网站| 亚洲欧美一区二区三区黑人 | 视频中文字幕在线观看| 一个人观看的视频www高清免费观看| 男人爽女人下面视频在线观看| 国产精品美女特级片免费视频播放器| 日本欧美国产在线视频| 亚洲精品国产av成人精品| a级毛色黄片| 国产在视频线在精品| 黄色日韩在线| 亚洲激情五月婷婷啪啪| 日韩中字成人| 国产一区二区三区av在线| 在线免费观看的www视频| 一夜夜www| 国产精品美女特级片免费视频播放器| 成人亚洲精品一区在线观看 | 国产高清不卡午夜福利| 国产真实伦视频高清在线观看| 69av精品久久久久久| 天堂影院成人在线观看| 22中文网久久字幕| 国国产精品蜜臀av免费| 亚洲一区高清亚洲精品| 97超碰精品成人国产| 搡女人真爽免费视频火全软件| 亚洲最大成人手机在线| 人妻一区二区av| 国产一区二区在线观看日韩| 亚洲欧美日韩卡通动漫| 国产老妇伦熟女老妇高清| 国产成人91sexporn| 99九九线精品视频在线观看视频| 一本久久精品| 欧美成人精品欧美一级黄| 欧美日本视频| 亚洲欧美一区二区三区黑人 | 欧美性感艳星| 亚洲国产色片| 熟妇人妻不卡中文字幕| 18+在线观看网站| 高清午夜精品一区二区三区| 日本与韩国留学比较| 亚洲精品色激情综合| 亚洲va在线va天堂va国产| 亚洲成人一二三区av| 一级爰片在线观看| 人人妻人人澡人人爽人人夜夜 | 免费av毛片视频| 免费观看性生交大片5| 99热全是精品| 精品国内亚洲2022精品成人| 一个人看的www免费观看视频| 国产三级在线视频| 国产欧美另类精品又又久久亚洲欧美| 嫩草影院新地址| 在线 av 中文字幕| 又爽又黄a免费视频| 人妻制服诱惑在线中文字幕| 男女边摸边吃奶| 男人舔女人下体高潮全视频| 午夜亚洲福利在线播放| 国产中年淑女户外野战色| 国产色爽女视频免费观看| 亚洲成色77777| 国产一区二区三区av在线| 久久久国产一区二区| 久久人人爽人人片av| 一级爰片在线观看| 国产白丝娇喘喷水9色精品| 亚洲av中文av极速乱| xxx大片免费视频| 成人无遮挡网站| 国产91av在线免费观看| 熟妇人妻久久中文字幕3abv| 久久热精品热| 国产精品人妻久久久影院| 国模一区二区三区四区视频| 亚洲av不卡在线观看| 亚洲无线观看免费| 日日啪夜夜撸| 一级a做视频免费观看| 能在线免费看毛片的网站| 亚洲人成网站在线播| 国产一区二区三区综合在线观看 | 国产男女超爽视频在线观看| 深夜a级毛片| 国产精品精品国产色婷婷| 亚洲欧美一区二区三区国产| 午夜精品一区二区三区免费看| 肉色欧美久久久久久久蜜桃 | 国产爱豆传媒在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产中年淑女户外野战色| 人人妻人人澡人人爽人人夜夜 | 国产亚洲av嫩草精品影院| 久久久久久久大尺度免费视频| 在线观看美女被高潮喷水网站| 久久精品国产亚洲网站| 中文字幕av在线有码专区| 国产av在哪里看| 午夜亚洲福利在线播放| 狠狠精品人妻久久久久久综合| 色综合色国产| 国产成人91sexporn| 2021少妇久久久久久久久久久| 亚洲天堂国产精品一区在线| 色网站视频免费| 国产av不卡久久| 亚洲国产色片| 免费av观看视频| 亚洲av国产av综合av卡| 国产精品女同一区二区软件| 在线免费观看不下载黄p国产| 免费av观看视频| 人人妻人人澡欧美一区二区| 天美传媒精品一区二区| 国产久久久一区二区三区| 国产精品久久久久久久电影| 久久人人爽人人爽人人片va| 狂野欧美白嫩少妇大欣赏| 91狼人影院| 国产 一区精品| 日本av手机在线免费观看| 熟妇人妻不卡中文字幕| 精品99又大又爽又粗少妇毛片| 亚洲色图av天堂| 美女黄网站色视频| 精品欧美国产一区二区三| 精品久久久久久久人妻蜜臀av| 亚洲精品,欧美精品| 最近中文字幕高清免费大全6| 欧美成人精品欧美一级黄| 国产男女超爽视频在线观看| 最近最新中文字幕免费大全7| 男人和女人高潮做爰伦理| 成人美女网站在线观看视频| 日韩精品青青久久久久久| 美女xxoo啪啪120秒动态图| 特大巨黑吊av在线直播| 精品国产露脸久久av麻豆 | 国产精品久久久久久精品电影| 亚洲最大成人手机在线| 80岁老熟妇乱子伦牲交| av播播在线观看一区| 午夜福利视频精品| 中文在线观看免费www的网站| 中文字幕av成人在线电影| 亚洲精品,欧美精品| eeuss影院久久| 午夜福利视频1000在线观看| 久久这里有精品视频免费| av播播在线观看一区| 天天躁日日操中文字幕| 免费少妇av软件| 蜜臀久久99精品久久宅男| 国产成人a∨麻豆精品| 国产成人a区在线观看| 亚洲精品久久久久久婷婷小说| 亚洲国产精品成人综合色| 乱码一卡2卡4卡精品| 婷婷色麻豆天堂久久| 国产免费视频播放在线视频 | 日本与韩国留学比较| 啦啦啦中文免费视频观看日本| 观看免费一级毛片| 亚洲精品乱码久久久久久按摩| 免费看美女性在线毛片视频| 亚洲美女搞黄在线观看| 日本免费在线观看一区| a级毛片免费高清观看在线播放| 干丝袜人妻中文字幕| 亚洲成人av在线免费| 日韩av免费高清视频| 中文字幕亚洲精品专区| 亚洲av日韩在线播放| 日本免费a在线| 中文字幕人妻熟人妻熟丝袜美| av网站免费在线观看视频 | 午夜爱爱视频在线播放| 熟妇人妻久久中文字幕3abv| 国产成人免费观看mmmm| 老司机影院成人| 国产成人一区二区在线| 久久精品人妻少妇| 激情五月婷婷亚洲| 免费看av在线观看网站| av天堂中文字幕网| 亚洲18禁久久av| 国产人妻一区二区三区在| 三级经典国产精品| 精品久久国产蜜桃| 免费电影在线观看免费观看| 国产精品久久视频播放| 女人久久www免费人成看片| 国产高清有码在线观看视频| 日韩国内少妇激情av| 听说在线观看完整版免费高清| 老师上课跳d突然被开到最大视频| 黄片无遮挡物在线观看| 婷婷色综合www| 三级经典国产精品| 免费不卡的大黄色大毛片视频在线观看 | 精品一区二区三区人妻视频| 国产精品女同一区二区软件| 日本三级黄在线观看| 自拍偷自拍亚洲精品老妇| 一本一本综合久久| 夫妻性生交免费视频一级片| 欧美性猛交╳xxx乱大交人| 啦啦啦啦在线视频资源| 全区人妻精品视频| 国产成人福利小说| 人妻系列 视频| 亚洲经典国产精华液单| 国产亚洲91精品色在线| 18禁动态无遮挡网站| 男女啪啪激烈高潮av片| 亚洲欧美清纯卡通| 亚洲欧美日韩东京热| 国产精品国产三级国产av玫瑰| 免费在线观看成人毛片| 精品酒店卫生间| 亚洲欧洲日产国产| 中文字幕制服av| 22中文网久久字幕| 美女内射精品一级片tv| 可以在线观看毛片的网站| 五月伊人婷婷丁香| 国产片特级美女逼逼视频| 高清午夜精品一区二区三区| 一级二级三级毛片免费看| 午夜精品一区二区三区免费看| 少妇人妻一区二区三区视频| 99久国产av精品国产电影| 国产毛片a区久久久久| 舔av片在线| 国产亚洲精品av在线| 中国国产av一级| 亚洲无线观看免费| 日韩强制内射视频| 国产成人午夜福利电影在线观看| xxx大片免费视频| 老师上课跳d突然被开到最大视频| 亚洲欧洲国产日韩| 国产成年人精品一区二区| 久久99热6这里只有精品| 97人妻精品一区二区三区麻豆| 久久久久久久大尺度免费视频| 亚洲欧美成人综合另类久久久| 晚上一个人看的免费电影| 最新中文字幕久久久久| 极品少妇高潮喷水抽搐| 亚洲av福利一区| 精品久久久噜噜| 亚洲内射少妇av| 毛片女人毛片| 一个人看的www免费观看视频| 亚洲欧美日韩卡通动漫| 高清视频免费观看一区二区 | 不卡视频在线观看欧美| 少妇人妻一区二区三区视频| 好男人在线观看高清免费视频| 联通29元200g的流量卡| 99久久精品热视频| 色视频www国产| 欧美最新免费一区二区三区| 国产精品精品国产色婷婷| 国产日韩欧美在线精品| ponron亚洲| videossex国产| 亚洲欧美成人精品一区二区| 日韩精品青青久久久久久| 99久久九九国产精品国产免费| 免费电影在线观看免费观看| 成人二区视频| 国产乱人视频| 最后的刺客免费高清国语| 赤兔流量卡办理| 国产不卡一卡二| 国产精品嫩草影院av在线观看| 欧美日本视频| 又大又黄又爽视频免费| 日韩欧美精品免费久久| 国产又色又爽无遮挡免| 深夜a级毛片| 亚洲综合精品二区| 成人一区二区视频在线观看| 欧美日韩在线观看h| 成人无遮挡网站| 久久精品国产亚洲av涩爱| 高清日韩中文字幕在线| 在线免费观看不下载黄p国产| 免费av观看视频| 少妇被粗大猛烈的视频| 美女内射精品一级片tv| 久久久久久久久久久丰满| av在线蜜桃| 国语对白做爰xxxⅹ性视频网站| 日韩,欧美,国产一区二区三区| 日韩强制内射视频| 亚洲丝袜综合中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 麻豆乱淫一区二区| 国内少妇人妻偷人精品xxx网站| 男插女下体视频免费在线播放| 国产黄片美女视频| 欧美性感艳星| 久久精品国产亚洲av涩爱| 丝袜美腿在线中文| 欧美日韩一区二区视频在线观看视频在线 | 成人亚洲精品av一区二区| 人人妻人人澡欧美一区二区| 麻豆av噜噜一区二区三区| av免费观看日本| 亚洲av电影不卡..在线观看| 亚洲成人中文字幕在线播放| 老司机影院成人| 我的老师免费观看完整版| 好男人在线观看高清免费视频| 国产成人精品久久久久久| 国内精品宾馆在线| 天堂中文最新版在线下载 | av免费观看日本| 欧美成人精品欧美一级黄| 男人舔女人下体高潮全视频| 欧美成人午夜免费资源| 国产 亚洲一区二区三区 |