• <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例
    有壹手——重新定義快修連鎖
    精品少妇内射三级| 一级片免费观看大全| 午夜福利视频在线观看免费| 麻豆乱淫一区二区| 黄色视频在线播放观看不卡| www日本在线高清视频| 日本av手机在线免费观看| 啦啦啦在线观看免费高清www| 在线观看免费午夜福利视频| 1024视频免费在线观看| 国产精品av久久久久免费| 免费高清在线观看日韩| 捣出白浆h1v1| 日韩,欧美,国产一区二区三区| 欧美变态另类bdsm刘玥| 成人午夜精彩视频在线观看| 国产xxxxx性猛交| 午夜福利视频在线观看免费| 亚洲少妇的诱惑av| 欧美亚洲日本最大视频资源| 一级黄片播放器| 水蜜桃什么品种好| 九色亚洲精品在线播放| 亚洲中文日韩欧美视频| 啦啦啦在线观看免费高清www| 菩萨蛮人人尽说江南好唐韦庄| 欧美人与性动交α欧美精品济南到| 亚洲一区二区三区欧美精品| 精品久久久久久久毛片微露脸 | 色播在线永久视频| 久久ye,这里只有精品| 国产视频首页在线观看| 9色porny在线观看| 亚洲情色 制服丝袜| 国产日韩一区二区三区精品不卡| 久久国产精品人妻蜜桃| 欧美国产精品一级二级三级| 大码成人一级视频| 亚洲av男天堂| 国产国语露脸激情在线看| 久久精品亚洲av国产电影网| 久久九九热精品免费| 大香蕉久久成人网| 久久性视频一级片| 中文字幕另类日韩欧美亚洲嫩草| 美女脱内裤让男人舔精品视频| 美女脱内裤让男人舔精品视频| 少妇裸体淫交视频免费看高清 | 国产xxxxx性猛交| 美女午夜性视频免费| 日本欧美视频一区| 久久国产精品男人的天堂亚洲| av欧美777| 亚洲精品久久成人aⅴ小说| 美女午夜性视频免费| 精品一品国产午夜福利视频| 亚洲av电影在线观看一区二区三区| 在线av久久热| 国产麻豆69| 少妇人妻 视频| 精品国产一区二区久久| 十分钟在线观看高清视频www| 国产一区亚洲一区在线观看| 亚洲国产中文字幕在线视频| 中国国产av一级| 91老司机精品| 大片免费播放器 马上看| 国产成人av激情在线播放| 婷婷丁香在线五月| 亚洲免费av在线视频| 桃花免费在线播放| 日韩av免费高清视频| 欧美另类一区| 在线观看国产h片| 中国国产av一级| 亚洲 欧美一区二区三区| 国产免费一区二区三区四区乱码| 欧美人与性动交α欧美软件| 国产黄色视频一区二区在线观看| 黄网站色视频无遮挡免费观看| 亚洲三区欧美一区| 欧美日韩黄片免| 男女午夜视频在线观看| 久久精品熟女亚洲av麻豆精品| 免费在线观看完整版高清| 国产成人精品久久久久久| 亚洲视频免费观看视频| 亚洲精品美女久久久久99蜜臀 | 欧美国产精品一级二级三级| 久久影院123| 久久青草综合色| 在现免费观看毛片| 久久久精品区二区三区| 伊人久久大香线蕉亚洲五| 涩涩av久久男人的天堂| 巨乳人妻的诱惑在线观看| 国产精品国产av在线观看| 欧美人与善性xxx| 国产有黄有色有爽视频| 你懂的网址亚洲精品在线观看| 亚洲精品日韩在线中文字幕| 欧美日韩福利视频一区二区| 成人国语在线视频| av网站在线播放免费| 成人影院久久| 真人做人爱边吃奶动态| 亚洲精品日本国产第一区| 欧美精品亚洲一区二区| 天堂俺去俺来也www色官网| 韩国精品一区二区三区| 在线观看一区二区三区激情| 亚洲,欧美精品.| 国产精品一区二区在线观看99| 久久精品国产亚洲av高清一级| 久久人妻福利社区极品人妻图片 | 久久ye,这里只有精品| 啦啦啦在线免费观看视频4| 免费久久久久久久精品成人欧美视频| 日韩大片免费观看网站| 后天国语完整版免费观看| 桃花免费在线播放| 久久影院123| 亚洲欧美清纯卡通| 国产视频一区二区在线看| 自线自在国产av| 国产午夜精品一二区理论片| 免费高清在线观看日韩| 日韩一本色道免费dvd| 免费在线观看黄色视频的| 日韩av免费高清视频| 日日摸夜夜添夜夜爱| 国产精品久久久久久精品电影小说| 秋霞在线观看毛片| 老熟女久久久| 免费看不卡的av| 亚洲成色77777| 亚洲av欧美aⅴ国产| 99香蕉大伊视频| 亚洲精品国产一区二区精华液| 久久青草综合色| 欧美性长视频在线观看| 国产精品国产三级国产专区5o| 伊人久久大香线蕉亚洲五| 国产真人三级小视频在线观看| 在线精品无人区一区二区三| 老司机影院毛片| 欧美日韩成人在线一区二区| 国产成人91sexporn| 啦啦啦在线观看免费高清www| av在线播放精品| 老司机影院成人| 日本五十路高清| kizo精华| 最近中文字幕2019免费版| 亚洲熟女精品中文字幕| 美女扒开内裤让男人捅视频| 看免费av毛片| 欧美性长视频在线观看| 日韩中文字幕欧美一区二区 | 欧美人与性动交α欧美精品济南到| 两人在一起打扑克的视频| 久久国产精品人妻蜜桃| 亚洲精品日本国产第一区| 亚洲国产成人一精品久久久| 美女午夜性视频免费| 中文字幕人妻丝袜一区二区| 最近最新中文字幕大全免费视频 | 亚洲精品第二区| 人人妻,人人澡人人爽秒播 | 麻豆av在线久日| 久久这里只有精品19| 免费一级毛片在线播放高清视频 | 高清欧美精品videossex| 国产精品99久久99久久久不卡| av线在线观看网站| 尾随美女入室| 在线看a的网站| 欧美日韩福利视频一区二区| 色综合欧美亚洲国产小说| 日日摸夜夜添夜夜爱| 久热爱精品视频在线9| 在线观看一区二区三区激情| 黄色一级大片看看| 国产成人精品久久二区二区91| 99精品久久久久人妻精品| 丰满人妻熟妇乱又伦精品不卡| 王馨瑶露胸无遮挡在线观看| 欧美日韩黄片免| 80岁老熟妇乱子伦牲交| 老司机影院毛片| 国产精品久久久av美女十八| 国产成人精品在线电影| 男女下面插进去视频免费观看| 日韩一卡2卡3卡4卡2021年| 人人澡人人妻人| netflix在线观看网站| 免费看av在线观看网站| 亚洲国产毛片av蜜桃av| 国产免费现黄频在线看| 1024视频免费在线观看| 国产成人91sexporn| 青青草视频在线视频观看| 9热在线视频观看99| 男女午夜视频在线观看| 男女床上黄色一级片免费看| 麻豆av在线久日| 日韩av免费高清视频| 精品卡一卡二卡四卡免费| 亚洲专区中文字幕在线| 伊人久久大香线蕉亚洲五| 精品视频人人做人人爽| 国产精品国产av在线观看| 欧美黄色片欧美黄色片| 18在线观看网站| av一本久久久久| 女性被躁到高潮视频| tube8黄色片| 老司机影院毛片| 国产精品偷伦视频观看了| 欧美精品av麻豆av| 侵犯人妻中文字幕一二三四区| 波多野结衣av一区二区av| 亚洲成色77777| 亚洲人成电影免费在线| 丰满迷人的少妇在线观看| 国产精品一二三区在线看| 国产欧美日韩一区二区三 | 搡老岳熟女国产| 成在线人永久免费视频| 一区在线观看完整版| 色视频在线一区二区三区| 久久人人爽人人片av| 日韩欧美一区视频在线观看| 男女之事视频高清在线观看 | 在线观看免费高清a一片| 麻豆乱淫一区二区| 赤兔流量卡办理| 美女大奶头黄色视频| 精品欧美一区二区三区在线| 一本综合久久免费| 新久久久久国产一级毛片| 成人亚洲欧美一区二区av| 国产精品秋霞免费鲁丝片| 精品少妇久久久久久888优播| 精品国产超薄肉色丝袜足j| 18禁国产床啪视频网站| 这个男人来自地球电影免费观看| 久久久精品94久久精品| 99国产精品一区二区三区| 国产免费又黄又爽又色| 少妇裸体淫交视频免费看高清 | 欧美日韩亚洲国产一区二区在线观看 | 97在线人人人人妻| 日韩熟女老妇一区二区性免费视频| 亚洲天堂av无毛| 精品少妇久久久久久888优播| 午夜视频精品福利| 麻豆乱淫一区二区| 菩萨蛮人人尽说江南好唐韦庄| 一级毛片电影观看| 成人国语在线视频| 精品欧美一区二区三区在线| 十分钟在线观看高清视频www| av一本久久久久| 午夜激情av网站| 考比视频在线观看| 高清黄色对白视频在线免费看| 两人在一起打扑克的视频| 丰满少妇做爰视频| 精品人妻一区二区三区麻豆| 精品国产乱码久久久久久男人| 菩萨蛮人人尽说江南好唐韦庄| 国产成人精品久久二区二区免费| 欧美 日韩 精品 国产| 亚洲,欧美,日韩| 免费在线观看黄色视频的| 亚洲免费av在线视频| netflix在线观看网站| 晚上一个人看的免费电影| 一二三四社区在线视频社区8| 久久精品亚洲av国产电影网| 精品一区二区三区av网在线观看 | 亚洲第一av免费看| 免费观看人在逋| 成年女人毛片免费观看观看9 | 成人国产av品久久久| 婷婷成人精品国产| 亚洲成人免费av在线播放| 亚洲欧美精品自产自拍| 久久久久国产一级毛片高清牌| 日本欧美国产在线视频| 久久国产精品男人的天堂亚洲| 国产av一区二区精品久久| 久久免费观看电影| 青春草亚洲视频在线观看| 老司机靠b影院| 人人妻人人澡人人爽人人夜夜| 国产成人91sexporn| 亚洲一区中文字幕在线| 成人三级做爰电影| 国产一区亚洲一区在线观看| 无限看片的www在线观看| 国产xxxxx性猛交| 美女大奶头黄色视频| 男女国产视频网站| 婷婷成人精品国产| av有码第一页| 国产精品亚洲av一区麻豆| 在现免费观看毛片| 蜜桃在线观看..| 欧美日韩视频精品一区| 我的亚洲天堂| 99久久人妻综合| 国产男女超爽视频在线观看| e午夜精品久久久久久久| 国产精品熟女久久久久浪| 午夜影院在线不卡| 久久免费观看电影| 两个人看的免费小视频| 丝袜在线中文字幕| 亚洲人成77777在线视频| 在线av久久热| 下体分泌物呈黄色| 亚洲第一青青草原| 国产一区二区三区av在线| 亚洲成人手机| 十八禁高潮呻吟视频| 久久99一区二区三区| 国产精品一区二区免费欧美 | 一边摸一边做爽爽视频免费| 久久鲁丝午夜福利片| 国产色视频综合| 青青草视频在线视频观看| www.自偷自拍.com| 欧美精品一区二区大全| xxxhd国产人妻xxx| 好男人电影高清在线观看| 亚洲欧美精品自产自拍| 男女国产视频网站| 18禁观看日本| 亚洲欧美精品自产自拍| 国产精品 欧美亚洲| 丝袜在线中文字幕| 国产av精品麻豆| 多毛熟女@视频| 亚洲中文字幕日韩| 午夜福利在线免费观看网站| 又粗又硬又长又爽又黄的视频| 97人妻天天添夜夜摸| 人人妻人人爽人人添夜夜欢视频| 99热国产这里只有精品6| 亚洲国产最新在线播放| av欧美777| 一区在线观看完整版| 亚洲欧美清纯卡通| 人妻一区二区av| 国产成人91sexporn| 久久人人爽人人片av| 国产精品九九99| 日韩av免费高清视频| 久久久久久人人人人人| 精品国产乱码久久久久久男人| 成人三级做爰电影| 国产精品一二三区在线看| 老司机影院毛片| 男女国产视频网站| 老熟女久久久| 日韩人妻精品一区2区三区| 高潮久久久久久久久久久不卡| 欧美av亚洲av综合av国产av| 丝袜人妻中文字幕| 午夜福利视频精品| 精品国产一区二区三区四区第35| 亚洲熟女精品中文字幕| 国产精品久久久久成人av| 亚洲国产精品999| 一级毛片 在线播放| 日韩一区二区三区影片| 黑丝袜美女国产一区| 青春草亚洲视频在线观看| 亚洲 国产 在线| 欧美精品一区二区大全| 啦啦啦在线免费观看视频4| 欧美在线黄色| 观看av在线不卡| 曰老女人黄片| 五月开心婷婷网| 日本猛色少妇xxxxx猛交久久| 国产高清视频在线播放一区 | 久久久久久久久久久久大奶| 脱女人内裤的视频| 搡老岳熟女国产| 日本av免费视频播放| 2021少妇久久久久久久久久久| 尾随美女入室| 日本午夜av视频| 中文乱码字字幕精品一区二区三区| www.av在线官网国产| 欧美日韩亚洲国产一区二区在线观看 | 极品人妻少妇av视频| 婷婷色综合www| 国产一区二区三区综合在线观看| 视频区图区小说| 欧美日韩av久久| 国产三级黄色录像| 多毛熟女@视频| 亚洲av成人不卡在线观看播放网 | 国产精品久久久久久人妻精品电影 | 亚洲熟女毛片儿| 国产亚洲精品久久久久5区| 欧美大码av| 亚洲国产毛片av蜜桃av| 制服诱惑二区| 精品福利观看| 精品久久久久久久毛片微露脸 | 欧美日韩一级在线毛片| 精品福利观看| 久久久国产精品麻豆| 伦理电影免费视频| 久久人人97超碰香蕉20202| 亚洲黑人精品在线| 国产成人免费无遮挡视频| 考比视频在线观看| 女人精品久久久久毛片| 午夜两性在线视频| 亚洲天堂av无毛| 男人爽女人下面视频在线观看| 大香蕉久久网| 欧美大码av| 午夜影院在线不卡| 日韩伦理黄色片| 日韩一区二区三区影片| 久久人妻福利社区极品人妻图片 | 午夜福利一区二区在线看| 七月丁香在线播放| 亚洲精品国产区一区二| 国产女主播在线喷水免费视频网站| 国精品久久久久久国模美| av视频免费观看在线观看| 国产成人免费观看mmmm| 欧美日韩黄片免| 久久久久久久久久久久大奶| 天天躁日日躁夜夜躁夜夜| 十八禁网站网址无遮挡| 午夜久久久在线观看| videos熟女内射| 一区在线观看完整版| 亚洲国产精品成人久久小说| 岛国毛片在线播放| www.熟女人妻精品国产| 欧美av亚洲av综合av国产av| 中文欧美无线码| 日韩 亚洲 欧美在线| 黄频高清免费视频| 日日爽夜夜爽网站| 宅男免费午夜| 国产精品人妻久久久影院| www.自偷自拍.com| 青草久久国产| 精品福利观看| 国产熟女欧美一区二区| 午夜福利视频精品| 18禁国产床啪视频网站| 亚洲图色成人| 18禁裸乳无遮挡动漫免费视频| 久久久久久免费高清国产稀缺| 国产高清不卡午夜福利| 久久综合国产亚洲精品| 热99久久久久精品小说推荐| 性少妇av在线| 飞空精品影院首页| 欧美精品人与动牲交sv欧美| 咕卡用的链子| 欧美日本中文国产一区发布| av天堂久久9| 欧美变态另类bdsm刘玥| 精品欧美一区二区三区在线| 国产爽快片一区二区三区| 纵有疾风起免费观看全集完整版| 久久久久久久精品精品| a级片在线免费高清观看视频| 啦啦啦在线观看免费高清www| 极品少妇高潮喷水抽搐| 观看av在线不卡| 曰老女人黄片| 久久精品国产亚洲av高清一级| 好男人视频免费观看在线| 国产在线免费精品| 青青草视频在线视频观看| 久久毛片免费看一区二区三区| 大片电影免费在线观看免费| 丁香六月天网| 午夜老司机福利片| 亚洲 国产 在线| bbb黄色大片| 久久久精品区二区三区| 亚洲精品久久久久久婷婷小说| 97精品久久久久久久久久精品| 久久中文字幕一级| 色播在线永久视频| 看十八女毛片水多多多| 天堂8中文在线网| 欧美黄色片欧美黄色片| 免费av中文字幕在线| 男人添女人高潮全过程视频| 爱豆传媒免费全集在线观看| 无遮挡黄片免费观看| 狂野欧美激情性xxxx| av在线app专区| 一级毛片 在线播放| 精品人妻一区二区三区麻豆| 菩萨蛮人人尽说江南好唐韦庄| 激情视频va一区二区三区| 亚洲国产欧美一区二区综合| 曰老女人黄片| 最近最新中文字幕大全免费视频 | 五月开心婷婷网| 国产免费视频播放在线视频| 国产亚洲欧美在线一区二区| 啦啦啦中文免费视频观看日本| 一二三四社区在线视频社区8| 人人澡人人妻人| 蜜桃在线观看..| 青春草亚洲视频在线观看| 亚洲成人国产一区在线观看 | 国产不卡av网站在线观看| 男女国产视频网站| 亚洲人成网站在线观看播放| svipshipincom国产片| 精品视频人人做人人爽| 婷婷丁香在线五月| 欧美97在线视频| 成人国语在线视频| 国产深夜福利视频在线观看| 一区在线观看完整版| 亚洲av综合色区一区| 又大又爽又粗| 国产视频一区二区在线看| 汤姆久久久久久久影院中文字幕| 少妇 在线观看| 汤姆久久久久久久影院中文字幕| 国产精品国产三级国产专区5o| 一区二区三区精品91| 久久国产精品影院| 丝袜人妻中文字幕| 久久久精品区二区三区| 国产一区二区三区av在线| 国产视频一区二区在线看| 免费高清在线观看视频在线观看| 亚洲成人免费电影在线观看 | 欧美精品啪啪一区二区三区 | 国产精品国产av在线观看| 国产精品 欧美亚洲| 午夜免费男女啪啪视频观看| 久久影院123| www.999成人在线观看| 人人妻人人添人人爽欧美一区卜| 午夜福利,免费看| 亚洲 欧美一区二区三区| 精品人妻熟女毛片av久久网站| 国产99久久九九免费精品| 欧美激情极品国产一区二区三区| 啦啦啦在线免费观看视频4| 国产亚洲欧美在线一区二区| 王馨瑶露胸无遮挡在线观看| 国产有黄有色有爽视频| 无限看片的www在线观看| 午夜免费成人在线视频| 性少妇av在线| 日本欧美视频一区| 午夜激情久久久久久久| 美女中出高潮动态图| 亚洲欧洲日产国产| 久久亚洲精品不卡| 中文字幕最新亚洲高清| 深夜精品福利| 国精品久久久久久国模美| 黄色一级大片看看| av国产精品久久久久影院| 老司机影院毛片| 久久久久网色| 久久天躁狠狠躁夜夜2o2o | 久久久久精品人妻al黑| 啦啦啦在线免费观看视频4| 国产精品一区二区在线观看99| 亚洲国产最新在线播放| 亚洲欧美精品自产自拍| 男人爽女人下面视频在线观看| 日韩人妻精品一区2区三区| 秋霞在线观看毛片| 香蕉国产在线看| 韩国精品一区二区三区| 十分钟在线观看高清视频www| 久久狼人影院| 欧美精品av麻豆av| 亚洲国产日韩一区二区| 精品国产国语对白av| 国产免费福利视频在线观看| 久热爱精品视频在线9| 成年人黄色毛片网站| 亚洲精品国产av蜜桃| 色婷婷av一区二区三区视频| 欧美在线黄色| 国产一区二区 视频在线| 国产成人一区二区在线| 日日夜夜操网爽| 一本一本久久a久久精品综合妖精| 亚洲国产欧美在线一区| 国产又爽黄色视频| 亚洲精品国产色婷婷电影| 国产成人一区二区在线| 精品高清国产在线一区| 黑丝袜美女国产一区| 国产老妇伦熟女老妇高清| 免费人妻精品一区二区三区视频|