• <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)記在櫛孔扇貝中的通用性研究*

    2018-08-31 08:19:08張廣明孫秀俊楊愛國劉志鴻周麗青劉寒苗
    漁業(yè)科學(xué)進(jìn)展 2018年4期
    關(guān)鍵詞:通用性扇貝微衛(wèi)星

    張廣明 孫秀俊 吳 彪 楊愛國 劉志鴻 周麗青 劉寒苗 趙 慶

    ?

    蝦夷扇貝EST-SSR標(biāo)記在櫛孔扇貝中的通用性研究*

    張廣明1,2,3孫秀俊1,2吳 彪1,2楊愛國1,2①劉志鴻1,2周麗青1,2劉寒苗1,2,3趙 慶1,2,3

    (1. 農(nóng)業(yè)農(nóng)村部海洋漁業(yè)可持續(xù)發(fā)展重點實驗室 中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所 青島 266071; 2. 青島海洋科學(xué)與技術(shù)試點國家實驗室海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實驗室 青島 266071; 3. 上海海洋大學(xué)水產(chǎn)與生命學(xué)院 上海 201306)

    為研究蝦夷扇貝()和櫛孔扇貝()的群體遺傳多樣性并分析其親緣關(guān)系,對在蝦夷扇貝中開發(fā)的60個EST-SSR位點在櫛孔扇貝中的通用性進(jìn)行了研究。結(jié)果顯示,21個位點可在櫛孔扇貝中擴增出特異性條帶,通用性比例達(dá)35.00%,其中,17個位點具有多態(tài)性,多態(tài)性比例達(dá)28.33%。在蝦夷扇貝和櫛孔扇貝群體中,平均等位基因數(shù)(N)分別為2.7647和2.3529;平均有效等位基因數(shù)(N)分別為1.9487和1.6350;平均觀測雜合度(H)為0.6314和0.3333;平均期望雜合度(H)為0.4569和0.3139;平均多態(tài)信息含量(PIC)為0.3726和0.2597;Nei's (1973)基因多樣性指數(shù)()分別為0.4493和0.3087;Shannon信息指數(shù)分別為0.7176和0.5041;固定指數(shù)(F)檢測發(fā)現(xiàn),7個位點偏離了Hardy-Weinberg平衡;遺傳分化系數(shù)(F)為0.2398。本研究開發(fā)的21對通用性EST-SSR標(biāo)記,為進(jìn)一步開展蝦夷扇貝和櫛孔扇貝遺傳多樣性分析、分子輔助育種、基因發(fā)掘和種質(zhì)資源評價奠定了基礎(chǔ)。

    蝦夷扇貝;櫛孔扇貝;EST-SSRs;微衛(wèi)星;多態(tài)性;通用性

    櫛孔扇貝()主要分布于我國北方近海、日本、韓國和俄羅斯等沿海。隨著養(yǎng)殖規(guī)模的擴大,出現(xiàn)了養(yǎng)殖產(chǎn)量下降和病害頻發(fā)等問題(王如才等, 2004)。而蝦夷扇貝()屬冷水性雙殼貝類,原產(chǎn)于日本北海道及本洲北部、俄羅斯千島群島的南部海域及朝鮮附近(常亞青等, 2008; 楊愛國等, 2004)。近年來,隨著養(yǎng)殖規(guī)模的擴大,蝦夷扇貝養(yǎng)殖群體出現(xiàn)了苗種成活率低、成體小型化、高死亡率和病害頻發(fā)等問題(張明明等, 2008; 楊培民等, 2007; 李云峰等, 2011)。蝦夷扇貝和櫛孔扇貝在育苗和養(yǎng)殖上均遇到了瓶頸和無法單獨解決的難題,而應(yīng)用遺傳學(xué)手段進(jìn)行良種培育,選育抗病、抗逆和生長快的新品種(品系)是實現(xiàn)扇貝養(yǎng)殖業(yè)可持續(xù)發(fā)展的有效途徑(Li, 2005)。

    EST-SSR側(cè)翼序列在物種之間高度保守,因此,EST-SSR引物可在物種之間通用。種間通用性分子標(biāo)記作為重要的遺傳學(xué)研究工具,對種間進(jìn)化關(guān)系研究、種間比較作圖、通用性遺傳圖譜構(gòu)建等具有重要的意義(王曉澗, 2012)。蝦夷扇貝和櫛孔扇貝標(biāo)記開發(fā)及應(yīng)用的工作已有大量研究,現(xiàn)已有ISSR技術(shù) (呂振明, 2006; 何斌等, 2007)、RAPD技術(shù)(滕麗莉等, 2005)、SRAP技術(shù)(程寧寧等, 2009)和AFLP技術(shù) (于濤等, 2011)研究了蝦夷扇貝、櫛孔扇貝的遺傳特性,但由于操作復(fù)雜或準(zhǔn)確率低等缺點限制了分子標(biāo)記的進(jìn)一步發(fā)展,而基因組EST(表達(dá)序列標(biāo)簽,Expressed Sequence Tags,ESTs)兼具基因組微衛(wèi)星標(biāo)記的孟德爾遺傳共顯性的遺傳特征和直接標(biāo)記功能基因的優(yōu)點,同時,相比其他分子標(biāo)記,可靠性更高,且擴增的目標(biāo)序列在基因的編碼區(qū)開發(fā)成本較低,遺傳信息含量豐富等優(yōu)點(劉芳,2006)。因此,本研究對60對蝦夷扇貝EST引物在櫛孔扇貝基因組中的通用性進(jìn)行了分析,對櫛孔扇貝和蝦夷扇貝的群體遺傳多樣性和群體遺傳結(jié)構(gòu)進(jìn)行了分析,以期為進(jìn)一步的分子輔助育種提供理論依據(jù)。

    1 材料和方法

    1.1 材料

    蝦夷扇貝(殼長為 10~12 cm)和櫛孔扇貝(殼長為7~8 cm)均于2015年6月取自山東煙臺長島養(yǎng)殖區(qū),從2個群體中分別隨機選取30個個體作為研究材料,活體運回實驗室,–80℃保存?zhèn)溆谩?/p>

    1.2 引物的設(shè)計與合成

    從已發(fā)表的文獻(xiàn)(Zhang, 2016; Sun, 2016)中選擇多態(tài)性高且擴增效果好的60對蝦夷扇貝微衛(wèi)星引物,引物由生工生物工程(上海)有限公司合成。

    1.3 模板DNA的制備

    模板DNA提取自扇貝的閉殼肌。取其閉殼肌組織100 mg剪碎,加入500 μl勻漿緩沖液(10 mmol/L Tris-HCl,pH 8.0;100 mmol/L EDTA,pH 8.0),混勻后加入終濃度為0.5%的SDS和50 μg/ml的蛋白酶K,55℃消化 3 h,然后分別用等體積的酚∶氯仿(1:1),氯仿∶異戊醇(24:1)抽提,2倍體積乙醇沉淀,ddH2O溶解?;蚪MDNA濃度由GENEQUANTpro (Eppendorf) RNA/DNA分析定量儀和瓊脂糖凝膠電泳檢測其質(zhì)量。模板DNA用滅菌的超純水稀釋成50 ng/μl,–20℃保存?zhèn)溆谩?/p>

    1.4 PCR擴增與電泳檢驗

    PCR反應(yīng)體系總體積為10 μl,包括10×Buffer 1 μl、Mg2+(25 mmol/L)0.6 μl、dNTPs (各2 mmol/L) 1 μl、上下游引物(10 μmol/L)各1μl、模板DNA 1μl、DNA聚合酶(Promega) 0.5 U,加超純水補足10 μl。反應(yīng)在Eppendorf擴增儀上進(jìn)行。PCR反應(yīng)程序: 95℃預(yù)變性2 min,然后95℃ 45 s,退火45 s,72℃延伸45 s,共35個循環(huán);最后72℃延伸10 min。擴增產(chǎn)物在 1%的瓊脂糖凝膠上電泳進(jìn)行預(yù)檢測,然后在8%的變性聚丙烯酰氨凝膠上進(jìn)行電泳分離,電泳結(jié)果銀染顯色。待凝膠干燥后使用掃描儀記錄電泳圖譜。

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

    根據(jù)每個個體產(chǎn)生的條帶位置確定其基因型,利用Cervus 3.0軟件計算多態(tài)信息含量(Polymorphism information content, PIC)、觀測雜合度(Observed heterozygosity,H)、期望雜合度(Expected heterozygosity, H)、等位基因數(shù)(Number of allele,N)、有效等位基因數(shù)(Effective number of allele,N)、Nei群體間的相似性系數(shù)和群體間遺傳距離以及統(tǒng)計量進(jìn)行遺傳多樣性的分析,統(tǒng)計顯著性水平用Bonferroni校正。

    有效等位基因數(shù):

    式中,p為第個等位基因的頻率。

    多態(tài)信息含量的計算方法:

    式中,為等位基因數(shù)目,pp分別為第和第個等位基因的頻率。

    平均雜合度觀測值:H= 觀察到的雜合個體總數(shù)/觀察到的個體總數(shù);

    固定指數(shù)(F):F=1 – (o/e)

    -統(tǒng)計量(F):F= σ2/(1 –)

    式中,為某等位基因在整個群體中的平均頻率; σ2為該等位基因在分群體之間的方差。

    Nei群體間的相似性系數(shù):

    式中,X、Y分別為和群體第個位點的等位基因頻率。

    Nei群體間遺傳距離:

    D= –ln

    加上偏差矯正后為:

    式中,X、Y分別為和群體第個位點的等

    位基因頻率。

    2 結(jié)果

    2.1 扇貝屬間通用性引物的篩選

    將60對蝦夷扇貝微衛(wèi)星引物在蝦夷扇貝和櫛孔扇貝基因組DNA中進(jìn)行PCR擴增。根據(jù)電泳條帶優(yōu)化反應(yīng)條件,調(diào)整退火溫度及時間,最終有21對引物可穩(wěn)定擴增,條帶之間較為清晰,引物通用性比例達(dá)35.00%,其中,17對引物的擴增產(chǎn)物具有多態(tài)性,占總數(shù)的28.33%(表1)。

    表1 在蝦夷扇貝和櫛孔扇貝基因組中擴增出多態(tài)性條帶的17對EST-SSRs

    Tab.1 Seventeen polymorphic EST-SSRs were amplified from the genomes of the scallops P. yessoensis and C. farreri.

    2.2 群體遺傳結(jié)構(gòu)

    17對EST引物在蝦夷扇貝、櫛孔扇貝群體中均能檢測到較高的多態(tài)性,但多態(tài)性水平有所不同(表2)。蝦夷扇貝、櫛孔扇貝群體等位基因數(shù)變化范圍都為2.00~4.00,平均等位基因數(shù)(N)分別為2.7647和2.3529;平均有效等位基因數(shù)(N)分別為1.9487和1.6350;觀測雜合度(H)范圍為0.0667~0.9667和0~0.9667,平均為0.6314和0.3333;期望雜合度(H)范圍為0.1588~ 0.6576和0.0333~0.6672,平均為0.4569和0.3139;多態(tài)信息含量(PIC)范圍為 0.1500~0.5724 和 0.0323~ 0.5814,平均為0.3726和0.2597。

    表2 17個EST-SSR位點在蝦夷扇貝、櫛孔扇貝各30個個體中的遺傳特征

    Tab.2 Characterization of 17 EST-SSRs in 30 scallops P. yessoensis and C. farreri

    F1:蝦夷扇貝群體;F2:櫛孔扇貝群體。統(tǒng)計顯著性水平用Bonferroni校正(Rice, 1989;=17),*:<0.05;**:<0.01

    F1:; F2:. Table-wide significance levels were applied, using the sequential Bonferroni technique (Rice, 1989;=17), *:<0.05; **:<0.01

    2.3 扇貝群體遺傳變異性分析

    蝦夷扇貝和櫛孔扇貝群體遺傳結(jié)構(gòu)變化和21對EST引物在2個群體內(nèi)偏離Hardy-Weinberg平衡程度見表2。從表2可以看出,從F參數(shù)看,2個群體中7對引物偏離Hardy-Weinberg平衡,其中,位點comp88965_c0在2個群體中均偏離了Hardy-Weinberg平衡。從-檢驗的數(shù)據(jù)來看(表 3),配對比較F值(F<0.05)。通過17對微衛(wèi)星位點計算總的遺傳分化系數(shù)(F)為0.2398,表明有23.98%的遺傳變異來自于群體之間,76.02%來自個體之間。此外,對F值的計算表明,除了comp100440_c1、contig_963940和contig_988858的F>0,其他14個位點F<0,呈現(xiàn)雜合子過剩狀態(tài)。2個扇貝群體在comp91359_c0、comp100440_c1、comp91909_c0、contig_1217920、contig_1443641、contig_293601、contig_963940、contig_ 915269、contig_763705、contig_988858 10個引物呈現(xiàn)一定程度的雜合子缺失。

    表3 扇貝2個群體17個微衛(wèi)星位點的-檢驗分析

    Tab.3 F-statistics for two scallop populations at seventeen microsatellite loci

    3 討論

    3.1 蝦夷扇貝EST-SSR位點的通用性研究

    EST-SSR標(biāo)記來源于基因的編碼區(qū),具有開發(fā)成本低、信息含量豐富、擴增效率高等優(yōu)點(Wang, 2009; 李云峰等, 2011)。本研究從蝦夷扇貝的EST文庫中分離和篩選了60對多態(tài)性位點,多態(tài)性比例達(dá)40.0%,這與太平洋牡蠣()EST-SSR的多態(tài)性比例相似(王艷紅等, 2007)。而劉博等(2011)研究表明,縊蟶()EST-SSR的多態(tài)性比例達(dá)35.0%,與本研究結(jié)果基本一致。由于微衛(wèi)星側(cè)翼序列在種間或?qū)匍g有高度保守性,因此,開發(fā)通用性位點應(yīng)用于近源物種成為可能。有學(xué)者研究表明,文蛤()微衛(wèi)星引物在彩虹明櫻蛤()和青蛤()的引物通用性分別為10.3%和6.9%(李曉英等, 2010),美洲紫海膽()EST-SSR在光棘球海膽()中的引物通用性為22.50% (李晶晶,2008)。而本研究中,蝦夷扇貝EST-SSR在櫛孔扇貝中的引物通用性比例達(dá)35.00%,這與楊璞等(2008)和李曉英等(2010)的研究結(jié)果相一致,但低于櫛孔扇貝 cSNP 引物在蝦夷扇貝中52.27%的通用性比例(李晶晶, 2008),原因可能是SNP標(biāo)記相比SSR重復(fù)序列具有更高的遺傳穩(wěn)定性,且更容易對數(shù)據(jù)進(jìn)行自動化和規(guī)模化分析,王曉澗等(2012)的研究結(jié)果與本研究結(jié)果相一致。

    3.2 蝦夷扇貝EST-SSR位點的開發(fā)及群體遺傳多樣性分析

    本研究中2個扇貝群體中的F大部分為負(fù)值,表明這些基因座存在遠(yuǎn)交或雜交現(xiàn)象,2個種群均處于雜合度過剩狀態(tài),這可能是純合致死或是群體內(nèi)自然選擇引起的(張海濱, 2005; 王宇等, 2012)。Hardy- Weinberg平衡和多態(tài)性信息含量是經(jīng)常用來衡量群體遺傳結(jié)構(gòu)變異程度高低的重要參數(shù)。本研究發(fā)現(xiàn),有7個位點偏離了平衡,可能是人為因素的干擾,導(dǎo)致種群內(nèi)基因缺失造成雜合子缺失。其中,6個位點為高度多態(tài)性基因座(PIC>0.5),占總位點數(shù)的35.29%;11個位點為中度多態(tài)性基因座(0.25H和H分別為0.6314和0.4569,高于在櫛孔扇貝群體中的H(0.3333)和o(0.3139),由此可知,蝦夷扇貝群體的遺傳多樣性略高于櫛孔扇貝群體。同時,這組數(shù)據(jù)高于利用AFLP研究不同地理群體蝦夷扇貝遺傳多樣性的平均雜合度,可能是SSR比AFLP技術(shù)更加靈敏導(dǎo)致,劉芳(2006)、薛明等(2006)的研究也說明這一點。

    Chang YQ, Chen XX, Ding J,. Genetic diversity in five scallop populations of the Japanese scallop(). Acta Ecologica Sinica, 2007, 27(3): 1145–1152 [常亞青,陳曉霞, 丁君, 等. 蝦夷扇貝()5個群體的遺傳多樣性. 生態(tài)學(xué)報, 2007, 27(3): 1145–1152]

    Chang YQ, Zhang CS, Cao XB,. Effect of morphometrical traits on weight traits in one-year old Yesso scallop. Journal of Dalian Fisheries University, 2008, 23(5): 330–334 [常亞青, 張存善, 曹學(xué)彬, 等. 1齡蝦夷扇貝形態(tài)性狀對重量性狀的影響效果分析. 大連海洋大學(xué)學(xué)報, 2008, 23(5): 330–334]

    Cheng NN, Yang AG, Liu ZH,. Analysis on heterosis of×by SRAP marker. Marine Sciences, 2009, 33(10): 107–111 [程寧寧, 楊愛國, 劉志鴻, 等. 櫛孔扇貝(♀)×蝦夷扇貝()子一代雜種優(yōu)勢的SRAP分析. 海洋科學(xué), 2009, 33(10): 107–111]

    He B, Yang AG, Wang QY,. ISSR analysis of the F1hybrids of scallop♀×. Journal of Dalian Fisheries University, 2007, 22(4): 273– 277 [何斌, 楊愛國, 王清印, 等. 櫛孔扇貝♀×蝦夷扇貝♂單對雜交子一代的ISSR分析. 大連海洋大學(xué)學(xué)報, 2007, 22(4): 273–277]

    Li JJ. Isolation and characterization of microsatellite markers in two species of bivalve and sea urchin. Doctoral Dissertation of Ocean University of China, 2008 [李晶晶. 兩種雙殼貝類和光棘球海膽微衛(wèi)星標(biāo)記的分離及特性研究.中國海洋大學(xué)博士研究生學(xué)位論文, 2008]

    Li L, Xiang JH, Liu X,. Construction of AFLP-based genetic linkage map for Zhikong scallop,Jones et Preston and mapping of sex-linked markers. Aquaculture, 2005, 245(1–4): 63–73

    Li XY, Dong ZG, Yan BL. A Study on the universality of partial shellfish microsatellites in GenBank inand. First International Conference on Cellular, Molecular Biology, Biophysics and Bioengineering, 2010 [李曉英, 董志國, 閻斌倫, 等. GenBank中部分貝類微衛(wèi)星標(biāo)記在青蛤和彩虹明櫻蛤的通用性研究. 第I屆細(xì)胞生物學(xué), 分子生物學(xué), 生物物理與生物工程國際會議, 2010]

    Li YF, Liu WD, Gao XG,. Construction of cDNA libraries from mantle and kidney of Japanese scallop() and ESTs analysis. Academic Annual Meeting of China Fisheries Society, 2011, 578–585 [李云峰, 劉衛(wèi)東, 高祥剛, 等. 蝦夷扇貝外套膜和腎臟組織cDNA文庫構(gòu)建以及EST的初步分析. 中國水產(chǎn)學(xué)會學(xué)術(shù)年會, 2011, 578–585]

    Liu B, Shao YQ, Teng SS,. Characterization, development and utilization of EST-derived microsatellities in. China Zoological Society, Shellfish China Institute of Oceanology and Limnology Society of the Ninth Congress and the Fifteenth National Symposium on Shellfish Abstract Set, 2011 [劉博, 邵艷卿, 滕爽爽, 等. 縊蟶EST-SSR分布特征及引物開發(fā)利用. 中國動物學(xué)會、中國海洋湖沼學(xué)會貝類學(xué)分會第九次會員代表大會暨第十五次學(xué)術(shù)討論會會議摘要集, 2011]

    Liu F. Genetic diversity of different populations ofusing AFLP and SSR molecular markers. Master¢s Thesis of Liaoning Normal University, 2006 [劉芳. 利用AFLP和SSR分子標(biāo)記研究不同地理蝦夷扇貝群體的遺傳多樣性. 遼寧師范大學(xué)碩士研究生學(xué)位論文, 2006]

    Liu XL, Chang YQ, Xiang JH,. Analysis of effects of shell size characters on live weight in Chinese scallop. Oceanologia et Limnologia Sinica, 2002, 33(6): 673–678 [劉小林, 常亞青, 相建海, 等. 櫛孔扇貝殼尺寸性狀對活體重的影響效果分析. 海洋與湖沼, 2002, 33(6): 673–678]

    Lv ZM. Cytogenetic and molecular analysis of hybridization betweenand. Doctoral Dissertation of Ocean University of China, 2006[呂振明. 櫛孔扇貝和蝦夷扇貝雜交過程中的細(xì)胞與分子遺傳學(xué)分析. 中國海洋大學(xué)博士研究生學(xué)位論文, 2006]

    Qi XY, Dong YH, Yao HH,. Identification of 30 microsatellite markers inand their transferability inand. Journal of Fisheries of China, 2013, 37(8): 1147–1154 [齊曉艷, 董迎輝, 姚韓韓, 等. 文蛤30個微衛(wèi)星標(biāo)記的開發(fā)及在斧文蛤和簾文蛤中的通用性檢測. 水產(chǎn)學(xué)報, 2013, 37(8): 1147–1154]

    Sun XJ, Zhang GM, Zhou LQ,. Development of genome-wide microsatellite genetic resources by using RAD tag sequencing in Yesso scallop. Conservation Genetic Resource, 2016, 8(2): 192–194

    Teng LL, Yang AG, Zhao F,. Analysis on heterosis of×by RAPD marker. Chinese High Technology Letters, 2005, 15(6): 97–101 [滕麗莉, 楊愛國, 趙峰, 等. 櫛孔扇貝×蝦夷扇貝子一代雜種優(yōu)勢的RAPD分析. 高技術(shù)通訊, 2005, 15(6): 97–101]

    Wang RC, Zheng XD.Progress of marine shellfishes culture in China and its prospect. Periodical of Ocean University of China(Natural Science), 2004, 34(5): 775–780 [王如才, 鄭小東. 我國海產(chǎn)貝類養(yǎng)殖進(jìn)展及發(fā)展前景. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 2004, 34(5): 775–780]

    Wang XL, Song B, Qiu XM,. Development of EST-SSRs in scallop () from sequence database. Conservation Genetics, 2009, 10(4): 1129–1131

    Wang XJ. Development of cSNP markers in Zhikong scallop () and their cross-species applications. Master¢s Thesis of Ocean University of China, 2012[王曉澗. 櫛孔扇貝()cSNP標(biāo)記的開發(fā)及其種間通用性研究. 中國海洋大學(xué)碩士研究生學(xué)位論文, 2012]

    Wang Y, Li L, Zhang SD,. Microsatellite segregation distortion analysis of the out bred, inbred and selfed families of the Bay scallop(). Marine Sciences, 2012, 36(8): 109–115 [王宇, 李莉, 張守都, 等. 海灣扇貝雜交、近交和自交家系的微衛(wèi)星標(biāo)記偏分離分析. 海洋科學(xué), 2012, 36(8): 109–115]

    Wang YH, Ren R, Yu ZN. Exploitation of EST-SSR markers for. China Zoological Society, Shellfish China Institute of Oceanology and Limnology society of the Eighth Congress and the Thirteenth National Symposium on Shellfish Abstract Set, 2007 [王艷紅, 任睿, 喻子牛. 太平洋牡蠣()的EST-SSR標(biāo)記開發(fā). 中國動物學(xué)會、中國海洋湖沼學(xué)會貝類學(xué)分會第八次會員代表大會暨第十三次全國貝類學(xué)術(shù)討論會論文摘要集, 2007]

    Xue M, Du XD, Huang RL,. Biochemicalgenetic variation within three wild populations of. Marine Science Bulletin, 2006, 25(1): 38–43 [薛明, 杜曉東, 黃榮蓮, 等. 文蛤三個野生種群的生化遺傳變異. 海洋通報, 2006, 25(1): 38–43]

    Yang AG, Wang QY, Liu ZH,. The hybrid between the scallopsandand the inheritance characteristics of its first filial generation. Marine Fisheries Research, 2004, 25(5): 1–5 [楊愛國, 王清印, 劉志鴻, 等. 櫛孔扇貝與蝦夷扇貝雜交及子一代的遺傳性狀. 漁業(yè)科學(xué)進(jìn)展, 2004, 25(5): 1–5]

    Yang P, Yang AG, Liu ZH,. The selection of universal microsatellite primers betweenandand its application in hybrid identification. Journal of Anhui Agricultural Sciences, 2008, 36(19): 8287–8289 [楊璞, 楊愛國, 劉志鴻, 等. 櫛孔扇貝和蝦夷扇貝通用微衛(wèi)星引物的篩選及其在雜種鑒定中的應(yīng)用. 安徽農(nóng)業(yè)科學(xué), 2008, 36(19): 8287–8289]

    Yang PM, Yang AG, Liu ZH,. Cryopreservation ofsperm and its application to hybridization. Journal of Shanghai Ocean University, 2007, 16(4): 351–356 [楊培民, 楊愛國, 劉志鴻, 等. 蝦夷扇貝精子的冷凍保存及其雜交試驗應(yīng)用研究. 上海海洋大學(xué)學(xué)報, 2007, 16(4): 351–356]

    Yu T, Yang AG, Zhou LQ,. Genetic variation of,and their hybrids. Journal of Fishery Sciences of China, 2011, 18(3): 574–580 [于濤, 楊愛國, 周麗青, 等. 櫛孔扇貝、蝦夷扇貝及其雜交子代的群體遺傳多樣性分析. 中國水產(chǎn)科學(xué), 2011, 18(3): 574–580]

    Zhang GM, Sun XJ, Zhou LQ,. Isolation and characterization of 44 microsatellite markers in Yesso scallopby transcriptome mining. Conservation Genetic Resource, 2016, 8(2): 186–188

    Zhang HB. Inbreeding effect and genetic improvement of Bay scallop,(Lamarck) .Master¢s Thesis of University of Chinese Academy of Sciences (Institute of Oceanology), 2005 [張海濱. 海灣扇貝近交生物學(xué)效應(yīng)和遺傳改良研究. 中國科學(xué)院大學(xué)研究生院(海洋研究所)碩士研究生學(xué)位論文, 2005]

    Zhang MM, Zhao W. Causes of death ofin China and control measures. China Fisheries, 2008(2): 65–66 [張明明, 趙文. 我國蝦夷扇貝死亡原因的探討及控制對策. 中國水產(chǎn), 2008(2): 65–66]

    (編輯 陳嚴(yán))

    Transferability of EST-SSR frominto

    ZHANG Guangming1,2,3, SUN Xiujun1,2, WU Biao1,2, YANG Aiguo1,2①, LIU Zhihong1,2, ZHOU Liqing1,2, LIU Hanmiao1,2,3, ZHAO Qing1,2,3

    (1. Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071; 2. Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266071; 3. College of Fisheries and Life Science, Shanghai Ocean University, Shanghai 201306)

    In this study, expressed sequence tag (EST)-SSR markers were developed to investigate the genetic relationship between two commercially important scallop species, Yesso scallopand Zhikong scallopA total of 60 EST-SSRs previously developed fromwere selected, and their cross-species amplification inwas analyzed. As a result, 21 pairs of EST-SSR primers showed unique PCR products. The interspecies transferability was calculated to be 35.00%. Among the 21 EST-SSR primers, 17 were polymorphic in the studied populations, which resulted in 28.33% transferability. In the two populations ofand, the number of alleles ranged from 2.00 to 4.00, with mean allele numbers of 2.7647 and 2.3529, respectively. The mean effective number of alleles was estimated to be 1.9487 and 1.6350, while the mean observed heterozygosity was 0.6314 and 0.3333, respectively. Similar to the observed heterozygosity, the mean expected heterozygosity was higher in(0.4569) than in(0.3139). For the two populations, the diversity index was consistently higher inthan inThe mean polymorphism information content was estimated to be 0.3726 and 0.2597, and Nei’s (1973) genediversity index was calculated to be 0.4493 and 0.3087, respectively. Furthermore, Shannon’s information index in the two populations was 0.7176 and 0.5041, respectively. The genetic identity between the two species was calculated to be 0.619, with a high genetic divergence between them (0.480). Among the polymorphic markers, seven loci significantly deviated from the Hardy-Weinberg equilibrium according to the average fixation index (F). The genetic differentiation index (F) between the two populations was estimated to be 0.2398. The EST-SSR markers developed for cross-species amplification ofandare important resources for the study of genetic diversity, marker-assisted breeding, gene discovery, and genetic evaluation of germplasm resources.

    ;; EST-SSRs; Microsatellites; Polymorphism; Transferability

    YANG Aiguo, E-mail: yangag@ysfri.ac.cn

    10.19663/j.issn2095-9869.20170503002

    S968.31

    A

    2095-9869(2018)04-0139-08

    * 國家自然科學(xué)基金(31602153)、青島海洋科學(xué)與技術(shù)國家實驗室海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實驗室開放基金 (2016LMFS-B02)、海水養(yǎng)殖教育部重點實驗室開放基金(KLM2018004)和山東省自然科學(xué)基金(ZR2016CQ32)共同資助[This work was supported by National Natural Science Foundation of China (31602153), Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, P. R. China (2016LMFS-B02), the Key Laboratory of Mariculture, Ministry of Education, Ocean University of China (KLM2018004), and Natural Science Foundation of Shandong Province (ZR2016CQ32)]. 張廣明,E-mail: zgm985623@163.com

    楊愛國,研究員,E-mail: yangag@ysfri.ac.cn

    2017-05-03,

    2017-05-16

    張廣明, 孫秀俊, 吳彪, 楊愛國, 劉志鴻, 周麗青, 劉寒苗, 趙慶. 蝦夷扇貝EST-SSR標(biāo)記在櫛孔扇貝中的通用性研究. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(4): 139–146

    Zhang GM, Sun XJ, Wu B, Yang AG, Liu ZH, Zhou LQ, Liu HM, Zhao Q. Transferability of EST-SSR frominto. Progress in Fishery Sciences, 2018, 39(4): 139–146

    猜你喜歡
    通用性扇貝微衛(wèi)星
    Improving polyp detection at colonoscopy: Non-technological techniques
    梅花鹿基因組微衛(wèi)星分布特征研究
    扇貝的眼睛在哪里
    基于元模型的通用性列控仿真平臺基礎(chǔ)環(huán)境研究
    烤扇貝:貝殼反復(fù)利用,殼比肉還貴
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    四川動物(2017年4期)2017-07-31 23:54:19
    拋丸機吊具的通用性設(shè)計以及拋丸器的布置
    基于轉(zhuǎn)錄組測序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國“一箭雙星”成功將“遙感衛(wèi)星二十一號”與“天拓二號視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    扇貝吃了一頭大象
    日韩人妻精品一区2区三区| 夜夜骑夜夜射夜夜干| 欧美另类亚洲清纯唯美| 正在播放国产对白刺激| 国产成人欧美在线观看 | 亚洲精品自拍成人| 久久香蕉激情| 国产人伦9x9x在线观看| 757午夜福利合集在线观看| 日韩中文字幕欧美一区二区| 日韩 欧美 亚洲 中文字幕| 午夜福利乱码中文字幕| 老司机影院毛片| 最近最新中文字幕大全免费视频| 欧美日韩福利视频一区二区| 午夜成年电影在线免费观看| 一本大道久久a久久精品| 国产精品亚洲一级av第二区| 久久久久久久大尺度免费视频| 国产在线一区二区三区精| 日韩视频在线欧美| 久久中文字幕人妻熟女| 日本av免费视频播放| 国产精品亚洲av一区麻豆| 国产aⅴ精品一区二区三区波| 久久精品国产综合久久久| 亚洲欧美日韩另类电影网站| 久久国产精品男人的天堂亚洲| 成年人免费黄色播放视频| 国产亚洲精品久久久久5区| 国产一区二区三区在线臀色熟女 | 在线永久观看黄色视频| 黄频高清免费视频| 亚洲熟女精品中文字幕| 亚洲情色 制服丝袜| 老司机深夜福利视频在线观看| 动漫黄色视频在线观看| 亚洲伊人色综图| 他把我摸到了高潮在线观看 | 在线十欧美十亚洲十日本专区| 99国产极品粉嫩在线观看| 日韩三级视频一区二区三区| 午夜老司机福利片| 丝袜美足系列| 少妇精品久久久久久久| 男女下面插进去视频免费观看| 国产国语露脸激情在线看| 一级片'在线观看视频| 日韩成人在线观看一区二区三区| 亚洲专区字幕在线| 一边摸一边抽搐一进一出视频| 国产1区2区3区精品| 国产高清videossex| 啪啪无遮挡十八禁网站| 大香蕉久久网| 国产精品九九99| 宅男免费午夜| 1024香蕉在线观看| 人妻一区二区av| 国产一区二区三区综合在线观看| 日本黄色日本黄色录像| 精品一区二区三区av网在线观看 | 国产欧美日韩一区二区精品| 成人三级做爰电影| 亚洲国产欧美网| 亚洲精品中文字幕在线视频| 99re6热这里在线精品视频| 人成视频在线观看免费观看| 可以免费在线观看a视频的电影网站| 久久 成人 亚洲| 国产老妇伦熟女老妇高清| 韩国精品一区二区三区| 欧美黑人欧美精品刺激| 18在线观看网站| 777米奇影视久久| 久久精品国产综合久久久| 精品一区二区三区视频在线观看免费 | 免费观看人在逋| 99热国产这里只有精品6| 法律面前人人平等表现在哪些方面| 国产亚洲精品一区二区www | 国产高清视频在线播放一区| 亚洲成人国产一区在线观看| 亚洲精品中文字幕一二三四区 | 午夜福利影视在线免费观看| 成人国语在线视频| 少妇精品久久久久久久| 午夜福利视频精品| 日日夜夜操网爽| 99国产精品99久久久久| 99九九在线精品视频| 一级a爱视频在线免费观看| 超碰成人久久| 精品国产亚洲在线| 欧美亚洲日本最大视频资源| 美女视频免费永久观看网站| 国产男女内射视频| 亚洲熟女精品中文字幕| 中文字幕人妻丝袜制服| 亚洲av国产av综合av卡| 天堂中文最新版在线下载| 亚洲视频免费观看视频| 在线永久观看黄色视频| 久久狼人影院| kizo精华| 岛国毛片在线播放| 久久中文字幕人妻熟女| 精品国产一区二区三区久久久樱花| 51午夜福利影视在线观看| 免费看十八禁软件| 久久午夜亚洲精品久久| 午夜老司机福利片| 2018国产大陆天天弄谢| 我的亚洲天堂| 手机成人av网站| av电影中文网址| 欧美日韩精品网址| 欧美激情 高清一区二区三区| 母亲3免费完整高清在线观看| 汤姆久久久久久久影院中文字幕| 欧美国产精品一级二级三级| 午夜免费成人在线视频| 一区在线观看完整版| 国产激情久久老熟女| 亚洲 欧美一区二区三区| 汤姆久久久久久久影院中文字幕| 亚洲全国av大片| 国产成人免费无遮挡视频| 一边摸一边抽搐一进一小说 | 在线观看一区二区三区激情| 一进一出抽搐动态| 大型黄色视频在线免费观看| 脱女人内裤的视频| 一边摸一边抽搐一进一小说 | 久热爱精品视频在线9| 午夜成年电影在线免费观看| 最近最新免费中文字幕在线| 捣出白浆h1v1| 免费看a级黄色片| 亚洲欧美一区二区三区久久| 黄片大片在线免费观看| 国产精品.久久久| 好男人电影高清在线观看| 人人澡人人妻人| 人人妻人人澡人人看| 免费观看av网站的网址| 1024香蕉在线观看| 国产高清视频在线播放一区| 在线观看舔阴道视频| 多毛熟女@视频| 一区在线观看完整版| 久久久久久久精品吃奶| 精品乱码久久久久久99久播| 韩国精品一区二区三区| 女性被躁到高潮视频| 丰满迷人的少妇在线观看| bbb黄色大片| 日本a在线网址| 99re6热这里在线精品视频| 日韩免费av在线播放| 怎么达到女性高潮| 日韩中文字幕视频在线看片| 69精品国产乱码久久久| av网站免费在线观看视频| 亚洲av美国av| 美女午夜性视频免费| 亚洲国产欧美日韩在线播放| 亚洲专区字幕在线| 无限看片的www在线观看| 久久99热这里只频精品6学生| 一进一出抽搐动态| 日韩大码丰满熟妇| 一夜夜www| 国产一区二区三区综合在线观看| 精品亚洲乱码少妇综合久久| 国产精品一区二区精品视频观看| 91九色精品人成在线观看| 国产人伦9x9x在线观看| 欧美国产精品一级二级三级| 在线看a的网站| 老司机影院毛片| 久久精品aⅴ一区二区三区四区| 无限看片的www在线观看| 高清毛片免费观看视频网站 | 欧美黑人欧美精品刺激| 美女高潮喷水抽搐中文字幕| kizo精华| 亚洲,欧美精品.| 亚洲中文av在线| 亚洲九九香蕉| 日韩大码丰满熟妇| 一进一出抽搐动态| 性少妇av在线| 男女高潮啪啪啪动态图| 精品高清国产在线一区| 欧美日韩av久久| 亚洲美女黄片视频| 国产在线视频一区二区| 精品少妇黑人巨大在线播放| 久久中文字幕一级| 亚洲一区二区三区欧美精品| 国产欧美日韩一区二区精品| 一本色道久久久久久精品综合| 久久人人爽av亚洲精品天堂| 黑人操中国人逼视频| 中文亚洲av片在线观看爽 | 老司机靠b影院| 十八禁高潮呻吟视频| kizo精华| 极品教师在线免费播放| 国产黄色免费在线视频| 操出白浆在线播放| aaaaa片日本免费| 亚洲综合色网址| 久久精品成人免费网站| 男女之事视频高清在线观看| 在线观看免费日韩欧美大片| 99热网站在线观看| 久久午夜综合久久蜜桃| 免费一级毛片在线播放高清视频 | 美女福利国产在线| av免费在线观看网站| 三上悠亚av全集在线观看| 啦啦啦免费观看视频1| 岛国毛片在线播放| 久久99一区二区三区| 久9热在线精品视频| 2018国产大陆天天弄谢| 99精品在免费线老司机午夜| 九色亚洲精品在线播放| av线在线观看网站| 蜜桃在线观看..| 肉色欧美久久久久久久蜜桃| 亚洲色图 男人天堂 中文字幕| 91麻豆av在线| 99国产精品99久久久久| 变态另类成人亚洲欧美熟女 | 自线自在国产av| 久热这里只有精品99| 国内毛片毛片毛片毛片毛片| 亚洲五月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看 | 亚洲色图综合在线观看| 国产亚洲精品一区二区www | 精品国产亚洲在线| 国产又爽黄色视频| 国产精品一区二区在线观看99| 亚洲人成77777在线视频| 欧美激情极品国产一区二区三区| 宅男免费午夜| 国产伦理片在线播放av一区| 日本欧美视频一区| 天天躁夜夜躁狠狠躁躁| 国产亚洲av高清不卡| 9热在线视频观看99| 一本大道久久a久久精品| 欧美日韩国产mv在线观看视频| 久久热在线av| 免费久久久久久久精品成人欧美视频| 在线观看免费日韩欧美大片| av网站免费在线观看视频| 夜夜爽天天搞| av网站免费在线观看视频| 天天操日日干夜夜撸| 69av精品久久久久久 | 50天的宝宝边吃奶边哭怎么回事| 变态另类成人亚洲欧美熟女 | 久久久国产精品麻豆| 欧美精品一区二区免费开放| 黑人猛操日本美女一级片| 欧美 亚洲 国产 日韩一| av天堂在线播放| 亚洲av美国av| 男女边摸边吃奶| 777米奇影视久久| 夜夜骑夜夜射夜夜干| 99国产极品粉嫩在线观看| 亚洲中文av在线| 国产成人一区二区三区免费视频网站| 精品第一国产精品| 老熟妇乱子伦视频在线观看| 少妇精品久久久久久久| 亚洲九九香蕉| 国产成人av教育| 在线观看舔阴道视频| 午夜福利视频在线观看免费| 男女免费视频国产| 最新在线观看一区二区三区| 成人精品一区二区免费| 飞空精品影院首页| 亚洲中文av在线| 人人妻人人澡人人看| 大片免费播放器 马上看| www.自偷自拍.com| 国产精品免费一区二区三区在线 | 一区福利在线观看| 久久亚洲真实| www.精华液| 国产精品98久久久久久宅男小说| 国产在线观看jvid| 肉色欧美久久久久久久蜜桃| 日韩中文字幕视频在线看片| 欧美激情久久久久久爽电影 | 日韩欧美一区视频在线观看| 久久亚洲真实| 久久久久精品国产欧美久久久| 下体分泌物呈黄色| 热99久久久久精品小说推荐| 欧美精品一区二区大全| 欧美成人免费av一区二区三区 | 国产区一区二久久| 91九色精品人成在线观看| 丁香欧美五月| 高清黄色对白视频在线免费看| 免费久久久久久久精品成人欧美视频| 色老头精品视频在线观看| 看免费av毛片| 女同久久另类99精品国产91| 99精品欧美一区二区三区四区| 成年人黄色毛片网站| 亚洲 欧美一区二区三区| 大型黄色视频在线免费观看| 成在线人永久免费视频| 狂野欧美激情性xxxx| 色婷婷久久久亚洲欧美| 色在线成人网| 亚洲九九香蕉| 老司机福利观看| 97人妻天天添夜夜摸| 1024香蕉在线观看| 久久午夜亚洲精品久久| 大片免费播放器 马上看| 精品人妻1区二区| 电影成人av| 欧美大码av| 999久久久精品免费观看国产| 精品人妻熟女毛片av久久网站| 后天国语完整版免费观看| 国产在线精品亚洲第一网站| 12—13女人毛片做爰片一| 俄罗斯特黄特色一大片| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 99国产精品免费福利视频| 操出白浆在线播放| 一本综合久久免费| 久久久精品免费免费高清| 母亲3免费完整高清在线观看| 国产日韩一区二区三区精品不卡| 午夜精品国产一区二区电影| 午夜福利免费观看在线| 亚洲成a人片在线一区二区| 激情在线观看视频在线高清 | 亚洲欧美色中文字幕在线| 欧美精品人与动牲交sv欧美| 免费av中文字幕在线| 婷婷成人精品国产| 久久中文字幕一级| 2018国产大陆天天弄谢| 999久久久精品免费观看国产| 岛国毛片在线播放| 成人国语在线视频| 色老头精品视频在线观看| 热99久久久久精品小说推荐| 黄片大片在线免费观看| 在线观看免费高清a一片| 国产野战对白在线观看| 男女之事视频高清在线观看| 丰满人妻熟妇乱又伦精品不卡| 日本黄色视频三级网站网址 | 午夜福利欧美成人| 国产av精品麻豆| 日本撒尿小便嘘嘘汇集6| 成人影院久久| 天天躁夜夜躁狠狠躁躁| 免费女性裸体啪啪无遮挡网站| 午夜久久久在线观看| 婷婷成人精品国产| 18禁美女被吸乳视频| 国产高清videossex| 女人被躁到高潮嗷嗷叫费观| 香蕉丝袜av| 久久人人爽av亚洲精品天堂| 国产亚洲av高清不卡| 精品视频人人做人人爽| 国产成人欧美| 国产成人精品久久二区二区免费| 黄片大片在线免费观看| 捣出白浆h1v1| 在线观看免费视频日本深夜| 999精品在线视频| 久久久久精品人妻al黑| 操美女的视频在线观看| 欧美激情久久久久久爽电影 | 欧美老熟妇乱子伦牲交| 成人黄色视频免费在线看| 中文字幕精品免费在线观看视频| 午夜福利一区二区在线看| 99热国产这里只有精品6| 亚洲午夜精品一区,二区,三区| 日韩三级视频一区二区三区| 亚洲第一欧美日韩一区二区三区 | 欧美黄色片欧美黄色片| 日韩有码中文字幕| 黄色视频在线播放观看不卡| 国产亚洲欧美在线一区二区| 午夜久久久在线观看| 十八禁高潮呻吟视频| 欧美日韩成人在线一区二区| 亚洲欧美一区二区三区久久| 亚洲精品一二三| 色婷婷av一区二区三区视频| 一级a爱视频在线免费观看| 国产在线免费精品| 少妇的丰满在线观看| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| svipshipincom国产片| 男人舔女人的私密视频| 亚洲国产欧美一区二区综合| 啦啦啦视频在线资源免费观看| 欧美黑人精品巨大| 日本wwww免费看| 国产亚洲av高清不卡| 国产免费现黄频在线看| 在线播放国产精品三级| 亚洲精品在线观看二区| 麻豆av在线久日| 色播在线永久视频| 国产伦理片在线播放av一区| 久久香蕉激情| av免费在线观看网站| 大陆偷拍与自拍| 男女高潮啪啪啪动态图| 无遮挡黄片免费观看| 成年人黄色毛片网站| 亚洲自偷自拍图片 自拍| 男女床上黄色一级片免费看| 亚洲第一欧美日韩一区二区三区 | 免费在线观看完整版高清| 亚洲精品乱久久久久久| 在线观看www视频免费| 日韩欧美免费精品| 久久人人爽av亚洲精品天堂| 国产在线观看jvid| 国产熟女午夜一区二区三区| 久久精品国产a三级三级三级| 久久国产精品男人的天堂亚洲| 高清毛片免费观看视频网站 | 精品视频人人做人人爽| 如日韩欧美国产精品一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品九九99| 多毛熟女@视频| 97在线人人人人妻| 亚洲色图 男人天堂 中文字幕| 下体分泌物呈黄色| 国产成人精品在线电影| 欧美久久黑人一区二区| 欧美精品高潮呻吟av久久| 亚洲色图 男人天堂 中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区三区综合在线观看| 国产亚洲午夜精品一区二区久久| 欧美+亚洲+日韩+国产| 国产黄频视频在线观看| 岛国在线观看网站| 一个人免费在线观看的高清视频| 精品第一国产精品| 亚洲第一av免费看| 一进一出抽搐动态| 婷婷成人精品国产| 免费看十八禁软件| 欧美性长视频在线观看| 午夜福利,免费看| 精品久久蜜臀av无| 飞空精品影院首页| 久久久久久久久久久久大奶| 精品国产国语对白av| 一二三四社区在线视频社区8| 色综合欧美亚洲国产小说| www.自偷自拍.com| 国产精品成人在线| 水蜜桃什么品种好| 香蕉丝袜av| 真人做人爱边吃奶动态| 80岁老熟妇乱子伦牲交| 国产精品亚洲av一区麻豆| 狠狠精品人妻久久久久久综合| 午夜福利在线观看吧| 国产精品熟女久久久久浪| 亚洲av电影在线进入| 国产精品电影一区二区三区 | 精品视频人人做人人爽| 热99re8久久精品国产| 欧美亚洲日本最大视频资源| 午夜福利在线观看吧| 成年版毛片免费区| 肉色欧美久久久久久久蜜桃| 男男h啪啪无遮挡| 久久av网站| 国产亚洲午夜精品一区二区久久| 免费看a级黄色片| 亚洲人成电影免费在线| 久久久久久免费高清国产稀缺| 高清毛片免费观看视频网站 | 深夜精品福利| 免费久久久久久久精品成人欧美视频| 少妇精品久久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 欧美 日韩 精品 国产| 亚洲精品自拍成人| 老司机亚洲免费影院| 国产三级黄色录像| 黄色怎么调成土黄色| 国产精品久久久av美女十八| 日本撒尿小便嘘嘘汇集6| 久久久久久久国产电影| 在线观看免费日韩欧美大片| 午夜激情av网站| 亚洲七黄色美女视频| 午夜91福利影院| videosex国产| 国产免费视频播放在线视频| 久久久水蜜桃国产精品网| 欧美久久黑人一区二区| 亚洲视频免费观看视频| 久久久精品94久久精品| av在线播放免费不卡| 亚洲人成电影免费在线| www日本在线高清视频| 国产男女超爽视频在线观看| 国产av国产精品国产| www.999成人在线观看| av天堂久久9| 国产日韩欧美在线精品| 狂野欧美激情性xxxx| 一个人免费看片子| aaaaa片日本免费| 久久精品国产a三级三级三级| 亚洲国产av新网站| 日本五十路高清| 久久久久久亚洲精品国产蜜桃av| 婷婷成人精品国产| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产区一区二| 久久精品91无色码中文字幕| 一本综合久久免费| 叶爱在线成人免费视频播放| 午夜日韩欧美国产| 久久精品人人爽人人爽视色| 男男h啪啪无遮挡| 老司机影院毛片| 午夜日韩欧美国产| 一个人免费在线观看的高清视频| 国产1区2区3区精品| 国产欧美日韩精品亚洲av| 50天的宝宝边吃奶边哭怎么回事| 国产精品免费一区二区三区在线 | 1024视频免费在线观看| 又紧又爽又黄一区二区| 色精品久久人妻99蜜桃| 五月天丁香电影| 天天操日日干夜夜撸| a级毛片黄视频| 国产精品久久久久久人妻精品电影 | 欧美+亚洲+日韩+国产| 国产成人av教育| 99国产极品粉嫩在线观看| 99久久精品国产亚洲精品| 波多野结衣av一区二区av| 大码成人一级视频| 成年动漫av网址| √禁漫天堂资源中文www| 久久毛片免费看一区二区三区| 成人影院久久| 精品人妻熟女毛片av久久网站| 啦啦啦中文免费视频观看日本| 国产午夜精品久久久久久| 一级毛片精品| 久久久国产精品麻豆| 在线观看舔阴道视频| 黄频高清免费视频| 女性被躁到高潮视频| 亚洲一码二码三码区别大吗| 水蜜桃什么品种好| 高清毛片免费观看视频网站 | 国产免费现黄频在线看| netflix在线观看网站| 国产伦人伦偷精品视频| 亚洲一区中文字幕在线| 最近最新免费中文字幕在线| 国产精品亚洲一级av第二区| 香蕉国产在线看| 丝瓜视频免费看黄片| 久久久精品国产亚洲av高清涩受| av天堂在线播放| 精品久久蜜臀av无| 欧美午夜高清在线| 色视频在线一区二区三区| 欧美黑人欧美精品刺激| 黄色怎么调成土黄色| 性高湖久久久久久久久免费观看| 女人爽到高潮嗷嗷叫在线视频| 黄频高清免费视频| 91麻豆av在线| 久久久国产欧美日韩av| 2018国产大陆天天弄谢| av视频免费观看在线观看| 人人妻人人添人人爽欧美一区卜| 视频区欧美日本亚洲| 国产亚洲精品一区二区www | 精品国产乱码久久久久久男人| 老司机亚洲免费影院| 久久精品成人免费网站|