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

    杜洛克豬全基因組連鎖不平衡分析

    2016-12-19 06:19:15刁淑琪羅元宇陳桂華陳贊謀李加琪
    廣東農(nóng)業(yè)科學(xué) 2016年11期
    關(guān)鍵詞:杜洛克連鎖等位基因

    刁淑琪,羅元宇,蔡 迪,陳桂華,陳贊謀,張 豪,李加琪,張 哲

    (華南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院,廣東 廣州 510642)

    杜洛克豬全基因組連鎖不平衡分析

    刁淑琪,羅元宇,蔡 迪,陳桂華,陳贊謀,張 豪,李加琪,張 哲

    (華南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院,廣東 廣州 510642)

    利用豬Illumina Porcine SNP60K芯片對(duì)福建某核心種豬場(chǎng)杜洛克豬216個(gè)個(gè)體進(jìn)行基因型檢測(cè),基于該高密度SNP芯片數(shù)據(jù),運(yùn)用Haploview軟件計(jì)算全基因組連鎖不平衡并構(gòu)建杜洛克豬連鎖不平衡圖譜。結(jié)果表明,該杜洛克豬群體不同染色體上相鄰標(biāo)記間r2存在波動(dòng),波動(dòng)范圍為0.46~0.59,相鄰標(biāo)記間的平均連鎖不平衡程度r2為0.52,SSC10的r2最低(平均為0.46),SSC6的r2最高(平均為0.59),連鎖不平衡水平隨著標(biāo)記間距的增加而衰減、變異程度隨之減小。該研究結(jié)果可為杜洛克豬遺傳分析及全基因組選擇研究提供參考。

    豬;連鎖不平衡;杜洛克;SNP

    連鎖不平衡(Linkage Disequilibrium,LD)又稱等位基因關(guān)聯(lián),是指同一條染色體上,等位基因間廣泛存在的非隨機(jī)組合現(xiàn)象[1]。連鎖不平衡是群體遺傳學(xué)研究的重要內(nèi)容之一,在分子生物學(xué)研究與動(dòng)物遺傳育種領(lǐng)域,群體基因組連鎖不平衡信息是數(shù)量性狀基因座定位、群體遺傳結(jié)構(gòu)分析和全基因組選擇的基礎(chǔ)。連鎖不平衡受選擇、突變和遺傳漂變等因素的影響。隨著科學(xué)技術(shù)的發(fā)展,覆蓋全基因組的高密度遺傳標(biāo)記已廣泛運(yùn)用于人和動(dòng)植物的連鎖不平衡研究,并取得了很大的研究進(jìn)展。在人類基因組學(xué)研究中,將單核苷酸多態(tài)性(Single Nucleotide Polymorphisms,SNP)作為遺傳標(biāo)記來(lái)深入研究連鎖不平衡已成為一種趨勢(shì)。結(jié)合關(guān)聯(lián)分析和功能研究,將包括復(fù)雜疾病及其中間表型在內(nèi)的人類各種表型特征,與基因組上以SNP為主的多態(tài)性位點(diǎn)聯(lián)系起來(lái),這對(duì)疾病預(yù)測(cè)、預(yù)防和治療等的突破具有重要作用[2]。

    目前,在畜禽方面有很多基于連鎖不平衡的研究。在牛的連鎖不平衡研究中,Mckay等[3]利用2 670個(gè)SNP標(biāo)記構(gòu)建了8個(gè)不同品種的連鎖不平衡圖譜,并發(fā)現(xiàn)其LD區(qū)域長(zhǎng)度不超過(guò)500 kb。Lu等[4]利用SNP標(biāo)記對(duì)安格斯牛、夏洛萊牛和雜交牛3個(gè)群體進(jìn)行LD分析,發(fā)現(xiàn)當(dāng)標(biāo)記間距小于30 kb時(shí),3個(gè)品種的LD分別為0.29、0.22和0.21,認(rèn)為不同品種和不同染色體間LD強(qiáng)度具有顯著差異。Ross等[5]以r2作為L(zhǎng)D的度量指標(biāo),利用SNP標(biāo)記估計(jì)了澳大利亞荷斯坦牛、澳洲安格斯牛、新西蘭黑白花奶牛和娟姍牛4個(gè)群體的LD,其大小分別為0.35、0.25、0.22和0.14。

    在豬的連鎖不平衡研究方面,Amaral等[6]通過(guò)分析中國(guó)家豬、歐洲家豬和歐洲野豬的371個(gè)SNP,發(fā)現(xiàn)歐洲豬種的LD程度高于中國(guó)豬種。Badke等[7]以r2作為度量指標(biāo)對(duì)長(zhǎng)白豬、約克夏、杜洛克和漢普夏的LD進(jìn)行比較,發(fā)現(xiàn)在標(biāo)記間距較小時(shí),杜洛克豬的LD較高;而標(biāo)記間距較大時(shí),漢普夏豬的LD較低。此外,Ai等[8]比較了中國(guó)地方豬和西方家豬的LD,發(fā)現(xiàn)西方家豬的LD高于中國(guó)地方豬。

    本研究以杜洛克豬為研究對(duì)象,基于SNP芯片基因組數(shù)據(jù)分析杜洛克豬群體的連鎖不平衡水平及其衰減規(guī)律,旨在為杜洛克豬的連鎖不平衡研究提供一定的參考。

    1 材料與方法

    1.1 試驗(yàn)材料

    試驗(yàn)對(duì)象為福建某核心種豬場(chǎng)的216頭杜洛克豬(公豬20頭,母豬196頭),試驗(yàn)豬均在常規(guī)條件下飼養(yǎng),體況正常。使用OMEGA Tissue DNA Kit對(duì)所有個(gè)體的耳組織抽提基因組DNA。采用Illumina Porcine SNP60K芯片[9]進(jìn)行SNP基因分型[10],共得到61 565個(gè)SNP。

    1.2 基因型質(zhì)量控制

    采用美國(guó)農(nóng)業(yè)部和愛(ài)荷華州立大學(xué)于2010年9月聯(lián)合組裝的豬基因組(Sus scrofa Build 10)作為參考序列(http://www.animalgenome. org/blast/),利用Plink軟件[11]對(duì)供試的216個(gè)個(gè)體和61 565個(gè)SNPs進(jìn)行基因型質(zhì)量控制。質(zhì)控標(biāo)準(zhǔn)為:個(gè)體檢出率(call rate)>0.95,SNP檢出率>0.95,最小等位基因頻率(Minor allele frequency,MAF)>0.05和哈迪-溫伯格平衡P值(Hardy-Weinberg equilibrium,HWE)≥10-6。經(jīng)過(guò)質(zhì)控后,剩余215個(gè)個(gè)體和32 186個(gè)有效SNP用于后續(xù)研究。

    1.3 LD度量

    連鎖不平衡程度通常用D'[12]和r2[13]進(jìn)行度量,本研究選用r2作為L(zhǎng)D的度量標(biāo)準(zhǔn)。r2表示兩個(gè)位點(diǎn)統(tǒng)計(jì)上和遺傳上相關(guān)程度的關(guān)系(0< r2<1),其表現(xiàn)較穩(wěn)定r2,對(duì)基因頻率的變化不敏感。r2的計(jì)算公式為:

    式中,PA1和PB1是兩個(gè)標(biāo)記位點(diǎn)上第1個(gè)等位基因的頻率,PA1B1是等位基因之間形成的單倍型頻率。

    采用Haploview軟件[15]計(jì)算各SNP之間的r2值,采用R軟件進(jìn)行后續(xù)的統(tǒng)計(jì)分析。

    2 結(jié)果與分析

    圖1 杜洛克豬質(zhì)量控制后SNP最小等位基因頻率分布

    2.1 基因型質(zhì)量控制

    供試的216個(gè)個(gè)體和61 565個(gè)SNPs進(jìn)行基因型質(zhì)量控制后,剩余215個(gè)個(gè)體和32 186個(gè)有效SNP進(jìn)行后續(xù)分析。質(zhì)控后的SNP最小等位基因頻率(MAF)分布情況見(jiàn)圖1。從圖1可以看出,MAF呈均勻分布,本研究中稀有突變位點(diǎn)極少。其中,70%以上SNP的MAF>0.20,說(shuō)明本研究中杜洛克豬群體中的SNP多態(tài)性較高,分型數(shù)據(jù)適用于LD的計(jì)算和圖譜的構(gòu)建。

    2.2 杜洛克豬各染色體LD

    杜洛克豬18條常染色體和X染色體r2的衰減趨勢(shì)如圖2所示。從圖2可以看出,相鄰標(biāo)記間的LD水平隨標(biāo)記間物理距離的增大而衰減。此外,除了SSC6、SSC14、SSC16、SSC17和SSC18等染色體外,其他染色體上個(gè)別間距較大的標(biāo)記之間也存在高r2現(xiàn)象。

    2.3 杜洛克豬全基因組LD及LD圖譜

    圖2 杜洛克豬各染色體的r2分布

    杜洛克豬全基因組r2的衰減趨勢(shì)如圖3所示,在全基因組水平上,LD變化趨勢(shì)與各染色體LD一致,即LD水平隨標(biāo)記間距的增加而衰減,LD的變異程度也隨之減小。但在某些距離較遠(yuǎn)的SNP標(biāo)記之間,LD程度也很高。

    本研究還構(gòu)建了杜洛克豬全基因組連鎖不平衡圖譜(圖4),圖4中顏色越深代表LD程度越強(qiáng)。結(jié)果表明,除了在相同染色體上的標(biāo)記間存在LD,不同染色體上的標(biāo)記間也存在一定程度的LD。

    2.4 杜洛克豬不同標(biāo)記間距的LD分布

    本研究統(tǒng)計(jì)了杜洛克豬各染色體上相鄰標(biāo)記間的r2分布情況,從表1可見(jiàn),杜洛克豬不同染色體上相鄰標(biāo)記間的r2存在波動(dòng),波動(dòng)范圍為0.46~0.59。

    對(duì)杜洛克豬的r2統(tǒng)計(jì)發(fā)現(xiàn),不同染色體上標(biāo)記間的LD程度存在差異。本研究分別統(tǒng)計(jì)了杜洛克豬在標(biāo)記間距為50 kb、500 kb、1 Mb和5 Mb時(shí)各染色體的LD分布,結(jié)果(表1)表明,SSC10的連鎖不平衡程度較低,SSC6的連鎖不平衡程度較高。

    圖3 杜洛克豬全基因組的r2分布

    圖4 杜洛克豬全基因組連鎖不平衡圖譜

    3 結(jié)論與討論

    本研究利用Illumina Porcine SNP60K芯片,對(duì)216頭杜洛克豬進(jìn)行了全基因組LD分析,揭示了該杜洛克豬群體中LD的變化規(guī)律,成功構(gòu)建了杜洛克豬的高密度連鎖不平衡圖譜,為杜洛克豬QTL定位、群體遺傳結(jié)構(gòu)分析及基因組選擇等研究提供了參考。

    在相鄰標(biāo)記(標(biāo)記間距平均為82 kb)之間,該杜洛克豬群體的平均LD程度為0.52,結(jié)果與Badke等[16]對(duì)4個(gè)外種豬的LD研究結(jié)果(r2=0.46,相鄰標(biāo)記間距約為70 kb)相似。在每條染色體之間,杜洛克的LD程度波動(dòng)范圍為0.46~0.59。這說(shuō)明杜洛克擁有較高的連鎖不平衡水平,推測(cè)其在進(jìn)化過(guò)程中受到高強(qiáng)度的人工選擇,長(zhǎng)期的選育也導(dǎo)致了杜洛克豬種的高度純化,從而使其LD程度較高。此外,本研究比較了標(biāo)記間距為50 kb、500 kb、1 Mb和5 Mb時(shí)各染色體的r2分布,發(fā)現(xiàn)SSC10的r2較小,SSC6的r2較大。這與Uimari等[17]對(duì)芬蘭長(zhǎng)白豬和芬蘭約克夏豬LD進(jìn)行分析的研究中SSC10 r2偏小的結(jié)果相同。同時(shí),也與Lei等[18]對(duì)3個(gè)丹麥豬種進(jìn)行LD分析研究得到SSC10的r2偏小的結(jié)果相同。

    表1 杜洛克豬各染色體上標(biāo)記間的r2分布情況

    LD隨著標(biāo)記間距的增加呈衰減趨勢(shì),標(biāo)記間r2逐漸降低,本研究結(jié)果與上述規(guī)律一致。但在本研究中發(fā)現(xiàn)某些距離較遠(yuǎn)的標(biāo)記間也存在強(qiáng)LD,這可能和參考基因組序列有關(guān)[19-20]。通過(guò)比較不同物種間相鄰標(biāo)記間的連鎖不平衡發(fā)現(xiàn),杜洛克豬的LD程度高于人類的LD水平[21-22],也高于綿羊的LD水平[23]。參考GWAS和GS的研究,將有用的LD界定為相鄰標(biāo)記間平均r2>0.3[24],這在本研究中SNP密度是充足的,但質(zhì)控后的SNP在基因組上呈非均勻分布,這可能會(huì)影響到SNP芯片在相關(guān)研究中的應(yīng)用效果。

    [1]Hill W G,Robertson A. Linkage disequilibrium in finite populations[J]. TAG. Theoretical and applied genetics. Theoretische und angewandteGenetik,1968,38(6):226-231.

    [2]Wall J D,Pritchard J K. Haplotype blocks and linkage disequilibrium in the human genome[J]. Nat Rev Genet,2003,4(8):587-597.

    [3]McKay S D,Schnabel R D,Murdoch B M,et al. Whole genome linkage disequilibrium maps in cattle[J]. BMC GENET,2007,8(1):74.

    [4]Lu D,Sargolzaei M,Kelly M,et al. Linkage disequilibrium in Angus,Charolais,and Crossbred beef cattle[J]. Frontiers in Genetics,2012,3.

    [5]de Roos Al APWD,Hayes B J,Spelman R J ,et al. Linkage Disequilibrium and Persistence of Phase in Holsten-Friesian,Jersey and Angus Cattle[J]. Genetics,2008,179(3):1503-1512.

    [6]Amaral A J,Megens H J,Crooijmans R P M A,et al. Linkage Disequilibrium Decay and Haplotype Block Structure in the Pig[J]. Genetics,2008,179(1):569-579.

    [7]Badke Y M,Bates R O,Ernst C W,et al Estimation of linkage disequilibrium in four US pig breeds[J]. Bmc Genomics,2012,1(13):1-10.

    [8]Ai H S,Huang L S,Ren J. Genetic Diversity,Linkage Disequilibrium and Selection Signatures in Chinese and Western Pigs Revealed by Genome-Wide SNP Markers[J]. Plos One,2013,2(8):e56001.

    [9]Ramos A M,Crooijmans R P M A,Affara N A,et al. Design of a High Density SNP Genotyping Assay in the Pig Using SNPs Identified and Characterized by Next Generation Sequencing Technology[J]. Plos One,2009,8(4):e6524.

    [10]Heaton M P,Harhay G P,Bennett G L,et al. Selection and use of SNP markers for animal identification and paternity analysis in U.S. beef cattle[J]. Mamm Genome,2002,13(5):272-281.

    [11]Purcell S,Neale B,Todd-Brown K,et al. PLINK:A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses[J]. The American Journal of Human Genetics,2007,81(3):559-575.

    [12]Lewontin R C. The Interaction of Selection and Linkage. I. General Considerations;Heterotic Models[J]. Gnentics,1964,1(49):49-67.

    [13]Hill W G. Estimation of linkage disequilibrium in randomly mating populations.[J]. Heredity,1974,33(2):229-239.

    [14]Zhao H,Nettleton D,Soller M,et al. Evaluation of linkage disequilibrium measures between multi-allelic markers as predictors of linkage disequilibrium between markers and QTL[J]. Genet Research,2005,86(1):77-87.

    [15]Barrett J C,F(xiàn)ry B,Maller J,et al. Haploview: analysis and visualization of LD and haplotype maps[J]. Bioinformatics,2005,21(2):263-265.

    [16]Badke Y M,Bates R O,Ernst C W,et al. Estimation of linkage disequilibrium in four US pig breeds[J]. Bmc Genomics,2012,1(13):1-10.

    [17]Uimari P,Tapio M. Extent of linkage disequilibrium and effective population size in Finnish Landrace and Finnish Yorkshire pig breeds[J]. J Anim Sci,2011,89(3):609-614.

    [18]Wang L,Rensen P S,Janss L,et al. Genomewide and local pattern of linkage disequilibrium and persistence of phase for 3 Danish pig breeds[J]. Bmc Genet,2013,1(14):1-11.

    [19]Farnir F,Coppieters W,Arranz J J,et al. Extensive Genome-wide Linkage Disequilibrium in Cattle[J]. Genome Res,2000,10(2):220-227.

    [20]Smith E M,Wang X,Littrell J,et al. Comparison of linkage disequilibrium patterns between the HapMap CEPH samples and a family-based cohort of Northern European descent[J]. Genomics,2006,88(4):407-414.

    [21]Dunning A M,Durocher F,Healey C S,et al. The extent of linkage disequilibrium in four populations with distinct demographic histories[J]. Am J Hum Genet,2000,67(6):1544-1554.

    [22]Kruglyak L. Prospects for whole-genome linkage disequilibrium mapping of common disease genes.[J]. Nat Genet,1999,2(22):139-144.

    [23]Zhao F P,Wang G K,Zeng T,et al. Estimations of genomic linkage disequilibrium and effective population sizes in three sheep populations[J]. Livest Sci,2014,170:22-29.

    [24]Ai H S,Huang L S,Ren J. Genetic Diversity,Linkage Disequilibrium and Selection Signatures in Chinese and Western Pigs Revealed by Genome-Wide SNP Markers[J]. Plos One,2013,8:e560012.

    (責(zé)任編輯 崔建勛)

    Genome-wide linkage disequilibrium analysis in Duroc pigherd

    DIAO Shu-qi,LUO Yuan-yu,CAI Di,CHEN Gui-hua,CHEN Zan-mou,ZHANG Hao,LI Jia-qi,ZHANG Zhe
    (College of Animal Science,South China Agricultural University,Guangzhou 510642,China)

    In this study,we collected 216 Duroc pigs from a breeding herd in Fujian province and genotyped them with Illumina Porcine SNP60K Bead Chip. Based on the high-density SNP data,the genome-wide LD was calculated with Haploview. Results showed that LD reduced as the marker intervals increasing. The r2between adjacent markers ranged from 0.46 to 0.59 in the Duroc population. The mean r2was 0.52 among adjacent markers across the genome,with the minimum in SSC10 (r2=0.46),and the maximum in SSC6 (r2=0.59). This study provides useful information for further genome-wide analysis of Duroc pig.

    pig;linkage disequilibrium;Duroc;SNP

    S828.2;Q343.1

    A

    1004-874X(2016)11-0116-06

    2016-08-10

    廣東省自然科學(xué)基金(2014A03031345);國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系項(xiàng)目(CARS-36);華南農(nóng)業(yè)大學(xué)大學(xué)生創(chuàng)新訓(xùn)練項(xiàng)目(201410564155)

    刁淑琪(1993-),女,在讀碩士生,E-mail:saradiao@126.com

    張哲(1984-),男,博士,副教授,E-mail:zhezhang@scau.edu.cn

    刁淑琪,羅元宇,蔡迪,等.杜洛克豬全基因組連鎖不平衡分析[J].廣東農(nóng)業(yè)科學(xué),2016,43(11):116-121.

    猜你喜歡
    杜洛克連鎖等位基因
    親子鑒定中男性個(gè)體Amelogenin基因座異常1例
    智慧健康(2021年17期)2021-07-30 14:38:32
    《全國(guó)優(yōu)秀公豬站》手冊(cè)(二)
    專注零售連鎖空間打造
    《全國(guó)優(yōu)秀公豬站》手冊(cè)(一)
    杜洛克與納雍糯谷豬雜交母系選育效果觀察
    庫(kù)里受傷的連鎖效應(yīng)
    NBA特刊(2018年7期)2018-06-08 05:48:32
    布拉格Burrito Loco連鎖快餐店
    WHOHLA命名委員會(huì)命名的新等位基因HLA-A*24∶327序列分析及確認(rèn)
    DXS101基因座稀有等位基因的確認(rèn)1例
    有壹手——重新定義快修連鎖
    午夜免费鲁丝| 成人亚洲欧美一区二区av| 97精品久久久久久久久久精品| 欧美日韩视频高清一区二区三区二| 国产欧美另类精品又又久久亚洲欧美| 搡女人真爽免费视频火全软件| 一区在线观看完整版| 亚洲天堂av无毛| 免费高清在线观看视频在线观看| 久久鲁丝午夜福利片| 国产伦精品一区二区三区四那| 精品亚洲成a人片在线观看| 黑人巨大精品欧美一区二区蜜桃 | 精品人妻偷拍中文字幕| 亚洲欧美一区二区三区黑人 | 麻豆成人av视频| 久久久欧美国产精品| 亚洲精华国产精华液的使用体验| 久久这里有精品视频免费| 亚洲不卡免费看| 久久久久久久国产电影| 久久精品国产亚洲av天美| 超碰97精品在线观看| 国产精品99久久99久久久不卡 | 免费人妻精品一区二区三区视频| 高清黄色对白视频在线免费看 | 精品一区在线观看国产| 99久久精品热视频| 免费人妻精品一区二区三区视频| 一级,二级,三级黄色视频| 2022亚洲国产成人精品| 日韩av在线免费看完整版不卡| 亚洲第一区二区三区不卡| 日本黄大片高清| 大陆偷拍与自拍| 亚洲成色77777| 大陆偷拍与自拍| 黑丝袜美女国产一区| 国产成人精品福利久久| 欧美老熟妇乱子伦牲交| 国产精品人妻久久久久久| 精品99又大又爽又粗少妇毛片| 中文字幕免费在线视频6| 久久久久久久大尺度免费视频| 国产爽快片一区二区三区| 日本91视频免费播放| 精品亚洲成国产av| 成人二区视频| 纯流量卡能插随身wifi吗| 男人狂女人下面高潮的视频| 在线观看www视频免费| 韩国av在线不卡| 丰满饥渴人妻一区二区三| 国产美女午夜福利| av不卡在线播放| 欧美日韩视频高清一区二区三区二| 久久青草综合色| 视频中文字幕在线观看| 国产高清有码在线观看视频| 一二三四中文在线观看免费高清| 人人澡人人妻人| 欧美日韩视频高清一区二区三区二| 免费观看性生交大片5| 国产高清有码在线观看视频| 午夜福利视频精品| 久久免费观看电影| 亚洲成人av在线免费| 黄色一级大片看看| 亚洲性久久影院| 国产在线男女| 男人添女人高潮全过程视频| 欧美一级a爱片免费观看看| 成人午夜精彩视频在线观看| 人妻 亚洲 视频| 国产精品一区二区性色av| av网站免费在线观看视频| 亚洲av男天堂| 久久97久久精品| 欧美丝袜亚洲另类| 中文在线观看免费www的网站| 欧美+日韩+精品| 91久久精品电影网| 亚洲av.av天堂| 91久久精品国产一区二区成人| 三级国产精品片| 成年女人在线观看亚洲视频| 国产熟女欧美一区二区| 欧美日韩精品成人综合77777| 亚洲国产欧美在线一区| 自拍欧美九色日韩亚洲蝌蚪91 | 丰满饥渴人妻一区二区三| 国产精品嫩草影院av在线观看| 国产淫片久久久久久久久| 日韩 亚洲 欧美在线| 少妇人妻精品综合一区二区| 建设人人有责人人尽责人人享有的| 晚上一个人看的免费电影| 熟妇人妻不卡中文字幕| 日本91视频免费播放| 日韩一区二区视频免费看| 国产精品人妻久久久久久| av有码第一页| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区欧美精品| 欧美日韩精品成人综合77777| av免费在线看不卡| 精品午夜福利在线看| a级毛片免费高清观看在线播放| 国产精品女同一区二区软件| 偷拍熟女少妇极品色| 校园人妻丝袜中文字幕| 一本久久精品| 97精品久久久久久久久久精品| 国产精品久久久久久精品古装| 妹子高潮喷水视频| 国产乱人偷精品视频| av视频免费观看在线观看| 91精品一卡2卡3卡4卡| 女性生殖器流出的白浆| 九九在线视频观看精品| 天天躁夜夜躁狠狠久久av| 草草在线视频免费看| 多毛熟女@视频| 日本爱情动作片www.在线观看| 人人妻人人添人人爽欧美一区卜| 久久青草综合色| 一本久久精品| 99热全是精品| 成年av动漫网址| 丝袜喷水一区| 尾随美女入室| 内射极品少妇av片p| 亚洲经典国产精华液单| av.在线天堂| 亚洲精品乱码久久久久久按摩| 国产精品国产三级专区第一集| 成人黄色视频免费在线看| 亚洲精品久久久久久婷婷小说| 国产欧美日韩综合在线一区二区 | 欧美精品国产亚洲| 国产视频首页在线观看| 亚洲,一卡二卡三卡| 在线观看免费高清a一片| 国产精品久久久久久av不卡| 久久国内精品自在自线图片| 久久国内精品自在自线图片| 秋霞在线观看毛片| 国产精品久久久久久久久免| 久久久国产欧美日韩av| 亚洲不卡免费看| 亚洲三级黄色毛片| a级片在线免费高清观看视频| 国产黄片美女视频| 亚洲欧美一区二区三区黑人 | 亚洲中文av在线| 国产成人精品久久久久久| 国产黄色免费在线视频| 99久久精品热视频| 久久久久人妻精品一区果冻| 精品卡一卡二卡四卡免费| 建设人人有责人人尽责人人享有的| 亚洲人成网站在线播| a级一级毛片免费在线观看| 中国美白少妇内射xxxbb| 欧美日韩综合久久久久久| 成人毛片a级毛片在线播放| 男人舔奶头视频| 亚洲伊人久久精品综合| 亚洲内射少妇av| videossex国产| 一级av片app| 国产亚洲av片在线观看秒播厂| 国产伦在线观看视频一区| 午夜福利在线观看免费完整高清在| 国产精品无大码| 久久久久久久久大av| 国产成人精品婷婷| 国产黄片视频在线免费观看| 婷婷色综合www| av播播在线观看一区| 99热国产这里只有精品6| 高清av免费在线| 亚洲成人手机| 麻豆精品久久久久久蜜桃| 国产一区有黄有色的免费视频| 伦精品一区二区三区| 黑丝袜美女国产一区| 久久99一区二区三区| 成人黄色视频免费在线看| 亚洲欧洲精品一区二区精品久久久 | 国产精品福利在线免费观看| 美女国产视频在线观看| 国产高清有码在线观看视频| av不卡在线播放| 国产精品伦人一区二区| 成人综合一区亚洲| 亚洲成色77777| 久久99精品国语久久久| √禁漫天堂资源中文www| 美女cb高潮喷水在线观看| 欧美 日韩 精品 国产| 免费看不卡的av| 天天操日日干夜夜撸| 国产成人免费观看mmmm| 一级毛片久久久久久久久女| 黄色视频在线播放观看不卡| 新久久久久国产一级毛片| 久久鲁丝午夜福利片| 肉色欧美久久久久久久蜜桃| 日韩不卡一区二区三区视频在线| 99九九线精品视频在线观看视频| 亚洲四区av| 国产91av在线免费观看| 日韩在线高清观看一区二区三区| 国产白丝娇喘喷水9色精品| 深夜a级毛片| 精品久久国产蜜桃| 国模一区二区三区四区视频| 男人狂女人下面高潮的视频| 亚洲无线观看免费| 亚洲av国产av综合av卡| 免费av中文字幕在线| 熟女人妻精品中文字幕| 欧美精品高潮呻吟av久久| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品日韩av片在线观看| 久久人人爽av亚洲精品天堂| 亚洲欧美日韩另类电影网站| 十八禁高潮呻吟视频 | 欧美另类一区| 国产成人a∨麻豆精品| 亚洲精品456在线播放app| 久久久精品免费免费高清| 高清黄色对白视频在线免费看 | 国产日韩欧美视频二区| 成人综合一区亚洲| 国产av精品麻豆| 日韩av不卡免费在线播放| 又大又黄又爽视频免费| 成人特级av手机在线观看| 男人爽女人下面视频在线观看| 国产在线一区二区三区精| 少妇高潮的动态图| 如日韩欧美国产精品一区二区三区 | 欧美日韩av久久| 91成人精品电影| 国产精品99久久久久久久久| 男人舔奶头视频| 久久午夜综合久久蜜桃| 曰老女人黄片| 大又大粗又爽又黄少妇毛片口| 黑丝袜美女国产一区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产男女超爽视频在线观看| 最黄视频免费看| 久久ye,这里只有精品| 国产在视频线精品| 狠狠精品人妻久久久久久综合| a级片在线免费高清观看视频| 交换朋友夫妻互换小说| 哪个播放器可以免费观看大片| 国产精品三级大全| 亚洲国产色片| 美女中出高潮动态图| 成人毛片a级毛片在线播放| 少妇的逼好多水| 人妻少妇偷人精品九色| 国产欧美日韩精品一区二区| 国产女主播在线喷水免费视频网站| 国产男人的电影天堂91| 一本一本综合久久| 又爽又黄a免费视频| 熟女电影av网| 亚洲av男天堂| 一二三四中文在线观看免费高清| 免费看av在线观看网站| 国产成人精品一,二区| 国产深夜福利视频在线观看| 男女国产视频网站| 午夜福利,免费看| 久久午夜福利片| 一级毛片黄色毛片免费观看视频| 亚洲人成网站在线观看播放| 欧美xxⅹ黑人| 这个男人来自地球电影免费观看 | 久久精品熟女亚洲av麻豆精品| 亚洲av中文av极速乱| 99热国产这里只有精品6| 一级毛片久久久久久久久女| 天堂俺去俺来也www色官网| 亚洲精华国产精华液的使用体验| 亚洲av国产av综合av卡| 色哟哟·www| 51国产日韩欧美| 蜜桃在线观看..| 91aial.com中文字幕在线观看| 又粗又硬又长又爽又黄的视频| 亚洲欧美一区二区三区国产| av专区在线播放| a级毛片在线看网站| 91久久精品国产一区二区三区| 少妇人妻精品综合一区二区| 日韩成人伦理影院| 成人黄色视频免费在线看| 国产av国产精品国产| 91精品国产九色| 精品酒店卫生间| 99久久精品热视频| 春色校园在线视频观看| 美女中出高潮动态图| 伊人亚洲综合成人网| 国国产精品蜜臀av免费| 欧美日韩av久久| 国产精品人妻久久久久久| 只有这里有精品99| 人人澡人人妻人| 国语对白做爰xxxⅹ性视频网站| 香蕉精品网在线| 久久婷婷青草| 欧美xxⅹ黑人| 伦理电影免费视频| 91精品国产九色| 五月开心婷婷网| 桃花免费在线播放| 亚洲国产精品成人久久小说| 蜜桃在线观看..| 69精品国产乱码久久久| av卡一久久| 久久精品国产a三级三级三级| 亚洲电影在线观看av| 亚洲国产欧美在线一区| 97精品久久久久久久久久精品| 国产成人一区二区在线| 久久精品国产自在天天线| av国产久精品久网站免费入址| 丝袜喷水一区| 欧美精品一区二区大全| 丰满乱子伦码专区| 中文字幕制服av| 肉色欧美久久久久久久蜜桃| 国产亚洲欧美精品永久| 国产中年淑女户外野战色| 高清视频免费观看一区二区| 少妇人妻 视频| 在现免费观看毛片| 赤兔流量卡办理| 婷婷色综合大香蕉| 一区二区三区乱码不卡18| 另类亚洲欧美激情| 亚洲精品日韩av片在线观看| 精品久久久精品久久久| 亚洲三级黄色毛片| 国产黄频视频在线观看| 少妇猛男粗大的猛烈进出视频| 免费人成在线观看视频色| 亚洲欧美精品专区久久| 亚洲精品久久午夜乱码| 一个人免费看片子| 午夜日本视频在线| 亚洲,一卡二卡三卡| av黄色大香蕉| 汤姆久久久久久久影院中文字幕| 精品久久久精品久久久| 成人午夜精彩视频在线观看| 亚洲av成人精品一二三区| 下体分泌物呈黄色| 麻豆成人av视频| 亚洲av.av天堂| 精品人妻熟女av久视频| av有码第一页| 一级a做视频免费观看| 午夜精品国产一区二区电影| 成人亚洲欧美一区二区av| 欧美另类一区| 久久久久久久久久久丰满| 久久免费观看电影| 精品少妇内射三级| 美女中出高潮动态图| 亚洲熟女精品中文字幕| 精品一区二区三区视频在线| 国产成人一区二区在线| 极品人妻少妇av视频| 秋霞在线观看毛片| 久久久国产精品麻豆| 伦理电影免费视频| 国产在线一区二区三区精| 男女啪啪激烈高潮av片| 18禁在线无遮挡免费观看视频| av专区在线播放| 亚洲怡红院男人天堂| 国产成人aa在线观看| 岛国毛片在线播放| 97在线视频观看| 国产91av在线免费观看| 纵有疾风起免费观看全集完整版| 啦啦啦中文免费视频观看日本| 少妇精品久久久久久久| 免费高清在线观看视频在线观看| 国语对白做爰xxxⅹ性视频网站| 日韩在线高清观看一区二区三区| 亚洲欧美成人综合另类久久久| 最后的刺客免费高清国语| 久久久欧美国产精品| 日韩制服骚丝袜av| 国产淫片久久久久久久久| 五月开心婷婷网| xxx大片免费视频| 女的被弄到高潮叫床怎么办| 男女免费视频国产| 亚洲精品久久午夜乱码| 久久久精品免费免费高清| 亚州av有码| 超碰97精品在线观看| 午夜福利在线观看免费完整高清在| 嫩草影院新地址| 成年女人在线观看亚洲视频| 少妇被粗大猛烈的视频| 大香蕉97超碰在线| 国产精品99久久久久久久久| 777米奇影视久久| 精品99又大又爽又粗少妇毛片| 边亲边吃奶的免费视频| 久久精品久久精品一区二区三区| 成年女人在线观看亚洲视频| 久久精品国产自在天天线| 人妻少妇偷人精品九色| 国产成人精品福利久久| 男女国产视频网站| 极品少妇高潮喷水抽搐| 自拍偷自拍亚洲精品老妇| 国产亚洲精品久久久com| 国产精品99久久久久久久久| 国产精品三级大全| 欧美成人午夜免费资源| 日日啪夜夜撸| 亚洲四区av| h日本视频在线播放| 黄色欧美视频在线观看| 日本与韩国留学比较| 人妻人人澡人人爽人人| 亚洲在久久综合| 欧美少妇被猛烈插入视频| 成年av动漫网址| 亚洲精品久久久久久婷婷小说| 免费观看的影片在线观看| 亚洲久久久国产精品| 欧美国产精品一级二级三级 | av网站免费在线观看视频| 日日爽夜夜爽网站| av视频免费观看在线观看| 综合色丁香网| 春色校园在线视频观看| 日韩成人伦理影院| 五月玫瑰六月丁香| 欧美 日韩 精品 国产| 大陆偷拍与自拍| av在线播放精品| 秋霞伦理黄片| 一本—道久久a久久精品蜜桃钙片| 嘟嘟电影网在线观看| 日韩av在线免费看完整版不卡| 插阴视频在线观看视频| 国产探花极品一区二区| 久久久精品94久久精品| 成人美女网站在线观看视频| 精品一区二区免费观看| 天堂俺去俺来也www色官网| 欧美精品一区二区大全| 99久国产av精品国产电影| 久久久久久久久久人人人人人人| 一级二级三级毛片免费看| 午夜影院在线不卡| 国产一区二区在线观看av| 超碰97精品在线观看| 男女啪啪激烈高潮av片| 国产老妇伦熟女老妇高清| 成人亚洲欧美一区二区av| 亚洲,一卡二卡三卡| kizo精华| 一级a做视频免费观看| 亚洲国产毛片av蜜桃av| 国产伦精品一区二区三区视频9| 国产亚洲一区二区精品| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人久久小说| 99视频精品全部免费 在线| 日本爱情动作片www.在线观看| 欧美高清成人免费视频www| 日本猛色少妇xxxxx猛交久久| 又黄又爽又刺激的免费视频.| 搡女人真爽免费视频火全软件| 三级国产精品片| 午夜日本视频在线| 欧美日韩av久久| 亚洲av日韩在线播放| 日韩免费高清中文字幕av| 伦理电影大哥的女人| 国产精品女同一区二区软件| 午夜免费男女啪啪视频观看| 免费久久久久久久精品成人欧美视频 | 女的被弄到高潮叫床怎么办| 国产午夜精品久久久久久一区二区三区| 少妇人妻精品综合一区二区| 制服丝袜香蕉在线| 免费观看的影片在线观看| 国产欧美亚洲国产| 久久97久久精品| 精品酒店卫生间| 男男h啪啪无遮挡| 伊人久久精品亚洲午夜| 欧美性感艳星| 91精品一卡2卡3卡4卡| 国产亚洲最大av| 久久狼人影院| 国产亚洲欧美精品永久| 国产色爽女视频免费观看| 日韩制服骚丝袜av| 18禁动态无遮挡网站| 精品国产乱码久久久久久小说| 丝袜在线中文字幕| 国产精品无大码| 国产精品伦人一区二区| 黄片无遮挡物在线观看| 热re99久久精品国产66热6| 亚洲精品自拍成人| 一级黄片播放器| 久久久午夜欧美精品| 性色av一级| 深夜a级毛片| 18禁动态无遮挡网站| 免费播放大片免费观看视频在线观看| 超碰97精品在线观看| 五月伊人婷婷丁香| 中文字幕制服av| 久久久久精品久久久久真实原创| 熟妇人妻不卡中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品456在线播放app| 婷婷色综合大香蕉| 纯流量卡能插随身wifi吗| 如日韩欧美国产精品一区二区三区 | 欧美日韩精品成人综合77777| 亚洲国产欧美日韩在线播放 | a级毛色黄片| 午夜av观看不卡| 如何舔出高潮| 国内揄拍国产精品人妻在线| 国产欧美日韩一区二区三区在线 | 久久久久久久久大av| 十八禁高潮呻吟视频 | av在线app专区| 国产一区二区在线观看av| 国产乱人偷精品视频| 99久久中文字幕三级久久日本| 高清在线视频一区二区三区| 女人精品久久久久毛片| av在线老鸭窝| 国产男人的电影天堂91| 日韩一区二区视频免费看| 多毛熟女@视频| 欧美老熟妇乱子伦牲交| 亚洲av电影在线观看一区二区三区| 高清视频免费观看一区二区| 久久精品久久精品一区二区三区| 人妻制服诱惑在线中文字幕| 少妇猛男粗大的猛烈进出视频| 老司机影院毛片| 亚洲四区av| 在线观看av片永久免费下载| 久久精品国产亚洲网站| 国产真实伦视频高清在线观看| 国产精品伦人一区二区| 欧美日韩精品成人综合77777| 国产欧美日韩一区二区三区在线 | 亚洲自偷自拍三级| 免费人成在线观看视频色| 搡老乐熟女国产| 亚洲无线观看免费| 久久97久久精品| 欧美日韩综合久久久久久| 在线 av 中文字幕| 午夜日本视频在线| 免费看av在线观看网站| 最新中文字幕久久久久| 涩涩av久久男人的天堂| 爱豆传媒免费全集在线观看| 精品视频人人做人人爽| 久久久久人妻精品一区果冻| 免费看光身美女| 亚洲高清免费不卡视频| 国产精品久久久久久久电影| 中文乱码字字幕精品一区二区三区| 我的女老师完整版在线观看| 97在线人人人人妻| 亚州av有码| 80岁老熟妇乱子伦牲交| 久久人人爽av亚洲精品天堂| 精品酒店卫生间| 日韩免费高清中文字幕av| 天天操日日干夜夜撸| 一级黄片播放器| av福利片在线| 新久久久久国产一级毛片| 啦啦啦在线观看免费高清www| 亚洲丝袜综合中文字幕| 欧美性感艳星| 国产白丝娇喘喷水9色精品| 大香蕉97超碰在线| 亚洲真实伦在线观看| 亚洲不卡免费看| 免费不卡的大黄色大毛片视频在线观看| 新久久久久国产一级毛片| 少妇的逼好多水| 久久久久久久久久成人| 久久午夜福利片|