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

    全基因組關(guān)聯(lián)分析及其在豬育種中的研究進(jìn)展

    2016-02-01 20:02:25宋志芳曹洪戰(zhàn)蘆春蓮
    豬業(yè)科學(xué) 2016年12期
    關(guān)鍵詞:表型基因組關(guān)聯(lián)

    宋志芳,曹洪戰(zhàn),2*,蘆春蓮,2

    (1.河北農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院,河北 保定 071000;2.河北農(nóng)業(yè)大學(xué)豬業(yè)科學(xué)研究所,河北 保定 071000)

    全基因組關(guān)聯(lián)分析及其在豬育種中的研究進(jìn)展

    宋志芳1,曹洪戰(zhàn)1,2*,蘆春蓮1,2

    (1.河北農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院,河北 保定 071000;2.河北農(nóng)業(yè)大學(xué)豬業(yè)科學(xué)研究所,河北 保定 071000)

    近年來,隨著高通量單核苷酸芯片和基因分型技術(shù)的不斷發(fā)展,利用全基因組關(guān)聯(lián)分析豬的性狀成為可能。全基因組關(guān)聯(lián)分析是一種新興的遺傳分析方法,能有效進(jìn)行復(fù)雜疾病和性狀的研究。國內(nèi)外相關(guān)研究人員針對豬性狀進(jìn)行全基因組關(guān)聯(lián)分析,積累了大量的單核苷酸多態(tài)性(SNP)標(biāo)記、候選基因以及數(shù)量性狀位點(diǎn),為豬分子育種提供基礎(chǔ)。該文主要對全基因組關(guān)聯(lián)分析的基本原理、分析方法以及對豬性狀的研究進(jìn)展進(jìn)行綜述。

    全基因組關(guān)聯(lián)分析;數(shù)量性狀位點(diǎn);SNP標(biāo)記;多重檢驗(yàn)

    我國是一個(gè)豬生產(chǎn)和消費(fèi)大國,國內(nèi)對豬肉的需求量巨大。最新數(shù)據(jù)顯示,2015年我國生豬出欄數(shù)為7.082 5億頭,豬肉產(chǎn)量5 487萬t。預(yù)計(jì)2016年我國生豬出欄規(guī)模是6.15億頭,農(nóng)業(yè)部4月份公布了我國能繁母豬存欄數(shù)是3 771萬頭。如此大規(guī)模的飼養(yǎng)和出欄規(guī)模,如果能夠很好地針對豬的重要經(jīng)濟(jì)性狀進(jìn)行遺傳改良,必然會(huì)帶來良好的經(jīng)濟(jì)效益。然而,有些數(shù)量性狀如生長性狀遺傳力低或難以準(zhǔn)確度量,利用傳統(tǒng)的育種方法耗時(shí)且進(jìn)展緩慢,準(zhǔn)確度不高。豬的60k高密SNP芯片已經(jīng)實(shí)現(xiàn)了商業(yè)化應(yīng)用,能使基因組基因定位的精確度提高,并且對豬的體型外貌[1]、生長[2-4]、肉質(zhì)[5-8]、繁殖[9-11]及行為[12]等多種性狀進(jìn)行全基因組關(guān)聯(lián)分析(GWAS),能夠進(jìn)一步深入了解性狀的遺傳基礎(chǔ)。目前全基因組關(guān)聯(lián)分析方法在豬性狀的研究方面已經(jīng)取得了一定的研究成果,積累了豐富的經(jīng)驗(yàn)和大量數(shù)據(jù),對豬育種工作奠定了基礎(chǔ)。

    1 全基因組關(guān)聯(lián)分析的提出

    關(guān) 聯(lián) 分 析 最 早 由N.Risch和Merikangas等人于1996年提出[13],認(rèn)為關(guān)聯(lián)分析的檢測效率要明顯高于連鎖分析,且檢測力更強(qiáng),并預(yù)測極有可能會(huì)應(yīng)用關(guān)聯(lián)分析大規(guī)模檢測復(fù)雜疾病的遺傳機(jī)制。在GWAS發(fā)展的初級階段,關(guān)聯(lián)分析多數(shù)被應(yīng)用于醫(yī)學(xué)上復(fù)雜疾病的研究。而GWAS是研究復(fù)雜疾病遺傳變異的最行之有效的方法,對疾病的遺傳機(jī)制在基因水平上有更深入地了解。Klein R J等人于2005年在《Science》雜志上報(bào)道了第1例與年齡相關(guān)的黃斑變性全基因組關(guān)聯(lián)分析研究[14],這是第一次利用GWAS發(fā)現(xiàn)復(fù)雜疾病的分子標(biāo)記,之后更多的人應(yīng)用GWAS方法研究人類復(fù)雜疾病和動(dòng)植物復(fù)雜性狀的遺傳機(jī)制,試圖探究遺傳機(jī)理為疾病的治療或性狀的改善提供科學(xué)的理論依據(jù)。利用GWAS方法在人類復(fù)雜疾病的遺傳研究方面日漸成熟,研究結(jié)果有助于人們逐漸了解人類復(fù)雜疾病和性狀的遺傳基礎(chǔ),再進(jìn)一步深入研究所檢測到的關(guān)聯(lián)位點(diǎn)和基因,以揭開人類復(fù)雜疾病和性狀的遺傳奧秘。

    2 全基因組關(guān)聯(lián)分析的基本原理

    GWAS是指在人類全基因組范圍內(nèi)找出存在的序列變異,即單核苷酸多態(tài)性(SNP),從中篩選出與人類疾病或動(dòng)植物復(fù)雜性狀相關(guān)的SNPs,可以利用其在全基因組水平上對疾病和復(fù)雜性狀的遺傳變異進(jìn)行關(guān)聯(lián)分析[15]?;驹硎牵阂赃B鎖不平衡(LD)為基礎(chǔ),通過識(shí)別數(shù)百個(gè)或數(shù)千個(gè)個(gè)體的定位群體中高密度的分子標(biāo)記,一般是上萬個(gè)甚至上百萬個(gè)SNPs標(biāo)記,剔除無關(guān)SNPs,篩選出與復(fù)雜性狀表型變異相關(guān)的分子標(biāo)記,然后以SNP為分子遺傳標(biāo)記,進(jìn)行全基因組水平上的相關(guān)分析,從而成為發(fā)現(xiàn)影響復(fù)雜性狀基因變異的一種新方案。動(dòng)物復(fù)雜性狀GWAS分析原理是以SNP為輔助,進(jìn)行總體關(guān)聯(lián)分析,首先利用SNP芯片在全基因組范圍內(nèi)進(jìn)行基因分型,比較個(gè)體之間每個(gè)遺傳變異及其頻率的差異,然后統(tǒng)計(jì)分析全部變異與研究性狀之間的關(guān)聯(lián)性,最后對最相關(guān)的遺傳變異加以驗(yàn)證,并根據(jù)驗(yàn)證結(jié)果確定其與研究性狀之間的關(guān)聯(lián)程度。

    GWAS的一般流程大致為:確定研究群體和樣本數(shù)量(樣本數(shù)量越多,分析結(jié)果越可靠),提取樣本DNA,利用芯片進(jìn)行SNP分型,對基因型數(shù)據(jù)進(jìn)行質(zhì)控,采用合適的數(shù)學(xué)統(tǒng)計(jì)模型對相關(guān)性狀表型值和SNPs進(jìn)行關(guān)聯(lián)分析,對關(guān)聯(lián)分析結(jié)果進(jìn)行數(shù)據(jù)分析及驗(yàn)證等。

    《Epidemiology》《Genetic》《Biometrics》等著名雜志從遺傳統(tǒng)計(jì)學(xué)角度對GWAS進(jìn)行了數(shù)據(jù)統(tǒng)計(jì)學(xué)方向的探討和研究,期望低成本、高效地找到遺傳標(biāo)記與疾病之間的聯(lián)系,GWAS分析過程中出現(xiàn)的假陽性問題得到解決。在表型選擇時(shí),要選擇遺傳度較高的疾病或者表型進(jìn)行檢測可提升遺傳學(xué)關(guān)聯(lián)研究的把握度[16],同時(shí)選擇表型的測定要盡量簡單、準(zhǔn)確。

    3 全基因組關(guān)聯(lián)分析的研究方法

    按照研究對象分類,可以分為基于無關(guān)個(gè)體的關(guān)聯(lián)分析和基于家系的關(guān)聯(lián)分析?;跓o關(guān)個(gè)體的關(guān)聯(lián)分析主要采用病例-對照研究設(shè)計(jì),大多用來研究質(zhì)量性狀;基于隨機(jī)人群的關(guān)聯(lián)分析,主要用來研究數(shù)量性狀;在研究基于家系的樣本時(shí),采用傳遞不平衡檢驗(yàn)(TDT)分析遺傳標(biāo)記與疾病數(shù)量表型和質(zhì)量表型的關(guān)聯(lián)可以排除人群混雜對于關(guān)聯(lián)分析的影響,但其在發(fā)現(xiàn)陽性關(guān)聯(lián)的檢驗(yàn)方面不如相同樣本量的病例對照研究有效。FBAT是運(yùn)用十分廣泛的基于家系的統(tǒng)計(jì)分析工具,能夠分析質(zhì)量性狀及數(shù)量性狀、調(diào)整混雜因素、分析基因-環(huán)境相互作用、分析單倍型、調(diào)整多重比較等。

    按照研究階段分類,分為單個(gè)階段研究和兩個(gè)或多個(gè)階段研究。單個(gè)階段研究即在有足夠大量的病例和對照樣本數(shù)量后,統(tǒng)一對全部選中的SNP進(jìn)行基因分型,然后分析每個(gè)SNP與疾病的基本關(guān)聯(lián),計(jì)算其關(guān)聯(lián)強(qiáng)度和OR值(比值比)。然而,因需要大量的樣本數(shù)量,單個(gè)階段研究基因分型耗費(fèi)的資金巨大。兩個(gè)或多個(gè)階段研究即采用小樣本數(shù)量進(jìn)行第一階段的全基因組范圍內(nèi)的SNP基因分型,統(tǒng)計(jì)分析后篩選少量陽性SNPs,在第二階段用更大數(shù)量的樣本對這些陽性SNPs進(jìn)行基因分型,最后整合兩個(gè)階段的結(jié)果進(jìn)行分析。研究證明DNA pool和微陣列試劑盒都能減少基因分型的工作量,能夠進(jìn)行低成本、高效益的SNP篩選。以GWAS作為研究工具能夠鑒定出與性狀相關(guān)聯(lián)的基因區(qū)域和DNA標(biāo)記。利用其先初步定位候選QTL,再借助其他研究手段進(jìn)一步精確定位數(shù)量性狀核苷酸(QTN)位點(diǎn)。因?yàn)檫M(jìn)行全基因組關(guān)聯(lián)分析研究的根本目的是要尋找影響性狀表型的QTN,QTN是QTL內(nèi)對數(shù)量性狀變異真正發(fā)揮作用的核苷酸多態(tài),能夠具體到DNA核苷酸水平,有助于深入地了解豬等動(dòng)物性狀的遺傳機(jī)制,揭示性狀差異的根本原因。

    4 全基因組關(guān)聯(lián)分析在豬育種中的應(yīng)用情況

    4.1 國外有關(guān)豬的GWAS研究現(xiàn)狀

    2005年公布了家豬的基因組序列[17],為以后的基因研究提供參考。2009年推出了豬的Illumina Porcine SNP60芯片,一經(jīng)推出國內(nèi)外科研人員利用豬的高密度芯片對豬重要性狀進(jìn)行全基因組關(guān)聯(lián)分析,發(fā)現(xiàn)了許多重要的SNPs、QTLs和候選基因,這些研究成果有助于未來豬分子育種的順利推進(jìn)。2012年《Nature》雜志報(bào)道了杜洛克豬的全基因組測序結(jié)果[18],為其性狀的改善提供參考。Fa等[1]利用豬的IlluminaPorcineSNP60芯片對820頭商品母豬的體型構(gòu)造、背膘厚、肢體穩(wěn)固性狀進(jìn)行定位,發(fā)現(xiàn)黑素皮質(zhì)素受體-4(MC4R)是影響背膘厚的候選基因,而IGF2是影響眼肌面積的重要候選基因,同時(shí)定位到了新的候選基因。Becker等[19]對192頭瑞典大白公豬的重要經(jīng)濟(jì)性狀進(jìn)行全基因組關(guān)聯(lián)分析分析,在14號(hào)和2號(hào)染色體上分別鑒定出1個(gè)pH1和胴體長的QTL。Wilson等[20]對豬的咬尾行為進(jìn)行全基因組關(guān)聯(lián)分析分析,在有咬尾行為的豬上定位到2個(gè)顯著的SNP。W.Luo等[7]以大白×民豬雜交豬為試驗(yàn)群體對其肉質(zhì)性狀進(jìn)行全基因組關(guān)聯(lián)分析,發(fā)現(xiàn)45個(gè)顯著SNPs與肉質(zhì)性狀有關(guān)。L.Fontanesi等[21]對意大利大白豬的背膘厚進(jìn)行全基因組關(guān)聯(lián)分析,鑒定出4個(gè)顯著SNPs與其背膘厚度有關(guān)。D.Becker等[22]對瑞士大白公豬進(jìn)行全基因組關(guān)聯(lián)分析,發(fā)現(xiàn)了4個(gè)QTLs可能與豬的肉質(zhì)、外觀、繁殖力和生長性狀緊密相關(guān)。

    4.2 中國國內(nèi)有關(guān)豬的GWAS研究現(xiàn)狀

    劉晨龍等人[23]對610頭杜長大豬、336頭二花臉豬和333頭萊蕪豬3個(gè)群體25種血液性狀進(jìn)行全基因組關(guān)聯(lián)分析,通過全基因組關(guān)聯(lián)分析和Meta分析共得到610個(gè)顯著影響杜長大豬、二花臉豬和萊蕪豬3個(gè)群體25種血液性狀的SNP位點(diǎn),初步確定F13A1、SPTA1、DBNL、SLC25A28和CTSC基因分別是嗜堿性粒細(xì)胞百分比(BASR)、紅細(xì)胞壓積(HCT)、淋巴細(xì)胞數(shù)(LYM)、平均紅細(xì)胞血紅蛋白濃度(MCHC)和中性粒細(xì)胞百分比(NEUR)的重要位置候選基因,為分析豬種的血液性狀或免疫性疾病提供重要參考。張哲等人[24]對191頭杜洛克豬的全基因組高密度SNP基因型數(shù)據(jù)以及生長性狀表型數(shù)據(jù)進(jìn)行全基因組關(guān)聯(lián)分析,結(jié)果表明在體重達(dá)100 kg日齡(D100)性狀上檢測到1個(gè)顯著關(guān)聯(lián)的SNP,在活體眼肌面積(LMA)性狀上檢測到6個(gè)染色體水平顯著關(guān)聯(lián)的SNP,都位于5號(hào)染色體上,分析表明BTG1和EFCAB6是影響生長性狀的兩個(gè)候選基因。楊慧等人[25]對288頭白色杜洛克×二花臉F2群體母豬殺嬰行為性狀進(jìn)行全基因組關(guān)聯(lián)分析,根據(jù)整個(gè)家系和60K SNP芯片信息推測出目標(biāo)個(gè)體的基因型,結(jié)果表明質(zhì)控后有效SNPs有34 591個(gè),同時(shí)在8號(hào)染色體檢測到了49個(gè)與母豬殺嬰行為顯著相關(guān)的SNPs。劉欣[26]選擇大白豬×民豬F2資源群體為研究對象,對其胴體性狀進(jìn)行全基因組關(guān)聯(lián)分析,共檢測到頭重、蹄重、板油重、屠宰重、胴體長以及背膘厚的顯著相關(guān)SNP位點(diǎn)數(shù)分別是63、84、56、8、78和35個(gè)。梁晶[27]通過GWAS在5號(hào)染色體上定位了與豬耳面積性狀顯著相關(guān)的區(qū)間,發(fā)現(xiàn)WIF1和LEMD3基因是影響豬耳性狀的有利候選基因。已經(jīng)利用全基因組關(guān)聯(lián)分析方法研究了豬的多種性狀,如生長性狀、胴體性狀、血液性狀、豬耳型性狀等,研究方法也日漸成熟,發(fā)現(xiàn)并定位了影響這些性狀的QTLs和SNPs,找到了與之相關(guān)的候選基因和主基因。

    5 全基因組關(guān)聯(lián)分析存在的問題

    雖然GWAS在人類疾病和動(dòng)植物相關(guān)性狀的研究成果有很多,也得到了一些重要的理論成果,有益于疾病的治療和性狀的改良。具有精確性高、統(tǒng)計(jì)效率高、比較范圍廣等優(yōu)點(diǎn),但是在應(yīng)用過程中也會(huì)不可避免地存在一些問題:1)在大樣本研究中,人群混雜是導(dǎo)致出現(xiàn)假陽性、假陰性結(jié)果的重要原因之一,數(shù)據(jù)結(jié)果存在誤差。利用統(tǒng)計(jì)分析手段控制人群混雜的影響、采用基于家系的研究等都能降低人群混雜對關(guān)聯(lián)分析結(jié)果的影響。為了降低假陽性位點(diǎn)出現(xiàn)的概率,要對群體進(jìn)行分層,主成分分析法、基因組控制法和結(jié)構(gòu)關(guān)聯(lián)法是處理群體分層的主要方法。2)使用GWAS時(shí)必須進(jìn)行多重檢驗(yàn),提高結(jié)果的可信度。常用的方法有Bonferroni校正法、控制錯(cuò)誤發(fā)生率法和Permutation法。3)對基因-變異-環(huán)境之間的相互關(guān)系進(jìn)行解釋,需要更多與變異相關(guān)的微效基因。4)不能獲取功能基因的相關(guān)信息,只能得到基因與表型的統(tǒng)計(jì)信息,具體相關(guān)基因的功能還需要后續(xù)分析。

    6 GWAS的展望和應(yīng)用前景

    利用GWAS能夠幫助人們更多地了解人類和動(dòng)植物復(fù)雜疾病和性狀的遺傳基礎(chǔ),發(fā)現(xiàn)大量的SNP標(biāo)記,篩選出很多性狀的相關(guān)候選基因,在染色體上定位QTL。利用GWAS能找出全基因組范圍內(nèi)所有變異的等位基因頻率。另外,也能發(fā)現(xiàn)許多從未發(fā)現(xiàn)的基因及染色體區(qū)域,為研究復(fù)雜疾病的發(fā)病機(jī)制和相關(guān)性狀的遺傳改良提供更多的線索。大量有關(guān)人類和動(dòng)植物上應(yīng)用GWAS進(jìn)行相關(guān)研究的文章被發(fā)表,涉及成百上千種性狀和疾病。但是利用GWAS并不能找到真正的功能基因,清楚這些基因作用于性狀產(chǎn)生多大的影響,故弄清引起變異的根本原因存在一定的困難,還需要對發(fā)現(xiàn)的QTLs和分子標(biāo)記進(jìn)行更多地研究和分析。此外,利用其進(jìn)行數(shù)據(jù)分析時(shí),結(jié)果存在一定的誤差,故其應(yīng)用還存在局限性。隨著基因芯片技術(shù)的發(fā)展和基因組研究的深入,相信在很長一段時(shí)間內(nèi)GWAS會(huì)得到越來越多的應(yīng)用,GWAS研究和計(jì)算方法也會(huì)得到不斷更新。

    [1] FAN B,ONTERU S K,DU Z Q,et al.Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs[J].PLoS One,2011,6(2):e14726.

    [2] FERN NDEZ A I,P REZ -MONTARELO D,BARRAG N C,et al.Genome-wide linkage analysis of QTL for growth and body composition employing the PorcineSNP60 BeadChip [J].BMC Genetics,2012,13(1):41.

    [3] FONTANESI L,SCHIAVO G,GALIMBERTI G,et al.A genome wide association study for backfat thickness in Italian Large White pigs highlights new regions affecting fat deposition includingneuronal genes[J]. BMC Genomics,2012,13(1):583.

    [4] SAHANA G,KADLECOV V,HORNSH J H,et al.A genomewide association scan in pig identifies novel regions associated with feed efficiency trait[J].Journal of animal science,2013,91(3):1041-1050.

    [5] RAMAYO-CALDAS Y, MERCAD A,CASTELL A,et al.Genome-wide association study for intramuscular fatty acid composition in an Iberian× Landrace cross[J].Journal of animal science,2012,90(9):2883-2893.

    [6] DUIJVESTEIJN N,KNOL E F,MERKS J W,et al.A genome-wide association study on androstenone levels in pigs reveals a cluster of candidate genes on chromosome 6[J].BMC Genetics,2010,11(1):1-11.

    [7] LUO W,CHENG D,CHEN S,et al.Genome-wide association analysis of meat quality traits in a porcine Large White× Minzhu intercross population[J].Int J Biol Sci,2012,8(4):580.

    [8] BECKER D,WIMMERS K,LUTHER H,et al.A genome wide association study to detect QTL for commercially important traits in Swiss Large White boars[J].PLoS One,2013,8(2):e55951.

    [9] ONTERU S K,F(xiàn)AN B,DU Z Q,et al.A whole genome association study for pig reproductive traits[J].Animal Genetics,2012,43(1):18-26.

    [10] ONTERU S K,GORBACH D M,YOUNG J M,et al.Whole genome association studies of residual feed intake and related traits in the pig[J].PLoS One,2013,8(6): e61756.

    [11] UIMARI P,SIRONEN A,SEV NAIMONEN M L.Whole-genome SNP association analysis of reproduction traits in the Finnish Landrace pig breed[J].Genetics Selection Evolution,2011,43(1):544-553.

    [12] WILSON K,ZANELLA R,VENTURA C,et al.Identification of chromosomal locations associated with tail biting and being a victim of tail-biting behaviour in the domestic pig(Sus scrofa domesticus)[J].Journal of applied genetics,2012,53(4): 449-456.

    [13] RISCH N,MERIKANGAS K.The future of genetic studies of comple×human diseases[J].Science,1996,273(5281):1516-1517.以下參考文獻(xiàn)若有需要可與作者聯(lián)系

    2016-10-11)

    宋志芳,E-mail:18730285576@163.com

    猜你喜歡
    表型基因組關(guān)聯(lián)
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    “一帶一路”遞進(jìn),關(guān)聯(lián)民生更緊
    奇趣搭配
    建蘭、寒蘭花表型分析
    智趣
    讀者(2017年5期)2017-02-15 18:04:18
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    慢性乙型肝炎患者HBV基因表型與血清學(xué)測定的臨床意義
    72例老年急性白血病免疫表型分析
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    天天操日日干夜夜撸| 夜夜爽天天搞| 亚洲国产成人一精品久久久| 色播在线永久视频| 多毛熟女@视频| 亚洲人成电影免费在线| 在线看a的网站| 一本一本久久a久久精品综合妖精| 亚洲第一av免费看| 好男人电影高清在线观看| 欧美国产精品va在线观看不卡| 精品国产乱子伦一区二区三区| 老司机午夜福利在线观看视频 | 热99久久久久精品小说推荐| 丝瓜视频免费看黄片| 国产精品久久久久久精品电影小说| 啦啦啦视频在线资源免费观看| 18禁裸乳无遮挡动漫免费视频| 久久九九热精品免费| 大陆偷拍与自拍| 黄网站色视频无遮挡免费观看| 精品一区二区三区av网在线观看 | 亚洲av片天天在线观看| 精品人妻熟女毛片av久久网站| 精品亚洲成国产av| 一级黄色大片毛片| 精品卡一卡二卡四卡免费| 在线播放国产精品三级| e午夜精品久久久久久久| 免费人妻精品一区二区三区视频| 女性生殖器流出的白浆| 午夜福利欧美成人| 亚洲专区中文字幕在线| 久久国产精品影院| 一本综合久久免费| 一级毛片女人18水好多| 成人三级做爰电影| 桃花免费在线播放| 在线观看免费日韩欧美大片| 国产成人精品久久二区二区91| 亚洲精品美女久久久久99蜜臀| 日韩免费高清中文字幕av| 制服人妻中文乱码| 老司机福利观看| 国产精品美女特级片免费视频播放器 | 制服诱惑二区| 最近最新中文字幕大全免费视频| 视频区欧美日本亚洲| 99国产精品免费福利视频| 一区二区日韩欧美中文字幕| 丝袜喷水一区| av线在线观看网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲专区国产一区二区| av不卡在线播放| 精品卡一卡二卡四卡免费| 真人做人爱边吃奶动态| 熟女少妇亚洲综合色aaa.| 99热网站在线观看| 午夜老司机福利片| 午夜老司机福利片| 国产成人av教育| 国产国语露脸激情在线看| 欧美日韩视频精品一区| cao死你这个sao货| 精品国产一区二区三区久久久樱花| 深夜精品福利| 久久亚洲精品不卡| 午夜福利视频在线观看免费| 亚洲精品在线美女| 成人18禁在线播放| 色精品久久人妻99蜜桃| 色94色欧美一区二区| 宅男免费午夜| 另类亚洲欧美激情| 日韩欧美一区视频在线观看| 欧美在线黄色| 免费看a级黄色片| 国产精品影院久久| 国产男靠女视频免费网站| 精品亚洲成a人片在线观看| 大型黄色视频在线免费观看| 久久久久久久国产电影| 午夜福利视频精品| 久久国产亚洲av麻豆专区| 国产av国产精品国产| 国产人伦9x9x在线观看| 午夜成年电影在线免费观看| 中文字幕色久视频| 国产一区二区在线观看av| 亚洲人成电影免费在线| 一区在线观看完整版| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品一二三| 精品一品国产午夜福利视频| 国产不卡一卡二| 精品视频人人做人人爽| 亚洲国产看品久久| 天天影视国产精品| 午夜福利欧美成人| av欧美777| 欧美在线黄色| 好男人电影高清在线观看| 日韩有码中文字幕| 亚洲人成电影免费在线| 久久九九热精品免费| 欧美成狂野欧美在线观看| 欧美成狂野欧美在线观看| 久久久国产成人免费| 免费在线观看黄色视频的| 如日韩欧美国产精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 在线看a的网站| 90打野战视频偷拍视频| 男女床上黄色一级片免费看| 久久毛片免费看一区二区三区| avwww免费| av天堂在线播放| 亚洲美女黄片视频| 亚洲成人免费av在线播放| 欧美国产精品va在线观看不卡| 日韩免费高清中文字幕av| 无限看片的www在线观看| 久久影院123| 国产精品 欧美亚洲| 国产片内射在线| 最近最新中文字幕大全电影3 | 天天操日日干夜夜撸| 天堂中文最新版在线下载| 免费观看人在逋| 叶爱在线成人免费视频播放| 一夜夜www| 男女免费视频国产| 成人免费观看视频高清| 免费在线观看完整版高清| 91精品三级在线观看| 大香蕉久久网| 一二三四社区在线视频社区8| 美女高潮到喷水免费观看| 国产在线免费精品| 他把我摸到了高潮在线观看 | 夜夜爽天天搞| 中亚洲国语对白在线视频| 久久精品成人免费网站| 如日韩欧美国产精品一区二区三区| 99精品久久久久人妻精品| 两性夫妻黄色片| 纵有疾风起免费观看全集完整版| 99国产精品一区二区蜜桃av | 亚洲精品美女久久av网站| 国产精品免费一区二区三区在线 | 免费高清在线观看日韩| 日韩视频一区二区在线观看| 俄罗斯特黄特色一大片| 精品福利永久在线观看| 久久免费观看电影| 91国产中文字幕| 18禁黄网站禁片午夜丰满| 日韩人妻精品一区2区三区| 一个人免费在线观看的高清视频| 叶爱在线成人免费视频播放| 国产精品.久久久| 人妻 亚洲 视频| 国产在视频线精品| 啦啦啦视频在线资源免费观看| 成年女人毛片免费观看观看9 | 桃花免费在线播放| 露出奶头的视频| videos熟女内射| 亚洲国产毛片av蜜桃av| 97在线人人人人妻| 亚洲国产毛片av蜜桃av| 国产xxxxx性猛交| 日本五十路高清| 自线自在国产av| 亚洲欧美色中文字幕在线| 久久精品熟女亚洲av麻豆精品| 国产精品久久久人人做人人爽| 人妻一区二区av| 久久久国产精品麻豆| 国产男女超爽视频在线观看| av有码第一页| 欧美精品一区二区大全| 90打野战视频偷拍视频| 国产av精品麻豆| 久久久久久久精品吃奶| 青草久久国产| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人av教育| 最黄视频免费看| 国产精品成人在线| 日本一区二区免费在线视频| 欧美日韩亚洲国产一区二区在线观看 | 最新美女视频免费是黄的| 免费不卡黄色视频| 欧美成狂野欧美在线观看| 亚洲人成伊人成综合网2020| 欧美中文综合在线视频| 高潮久久久久久久久久久不卡| 人人妻人人添人人爽欧美一区卜| 国产激情久久老熟女| 正在播放国产对白刺激| 最近最新中文字幕大全电影3 | 久久免费观看电影| 久久99一区二区三区| 女警被强在线播放| 在线观看免费午夜福利视频| 大片电影免费在线观看免费| 人人妻人人添人人爽欧美一区卜| 久久精品亚洲精品国产色婷小说| 老司机在亚洲福利影院| 亚洲欧美精品综合一区二区三区| 欧美精品亚洲一区二区| 在线亚洲精品国产二区图片欧美| 日韩视频在线欧美| 国产欧美亚洲国产| 久久中文字幕人妻熟女| cao死你这个sao货| 国产熟女午夜一区二区三区| 首页视频小说图片口味搜索| 9色porny在线观看| 欧美乱码精品一区二区三区| 亚洲午夜理论影院| 女性被躁到高潮视频| 精品一区二区三区视频在线观看免费 | 97在线人人人人妻| 久久精品aⅴ一区二区三区四区| 亚洲中文av在线| 日韩人妻精品一区2区三区| 免费在线观看黄色视频的| 香蕉国产在线看| 黄色片一级片一级黄色片| 免费久久久久久久精品成人欧美视频| 怎么达到女性高潮| 亚洲一区中文字幕在线| 黄色 视频免费看| 亚洲精品久久成人aⅴ小说| 午夜福利免费观看在线| 国产精品美女特级片免费视频播放器 | 欧美日韩黄片免| 亚洲视频免费观看视频| 男女下面插进去视频免费观看| 成年女人毛片免费观看观看9 | 一级片免费观看大全| 精品久久蜜臀av无| 亚洲 欧美一区二区三区| 在线亚洲精品国产二区图片欧美| 欧美老熟妇乱子伦牲交| 18禁美女被吸乳视频| 亚洲男人天堂网一区| 精品一区二区三区av网在线观看 | 自拍欧美九色日韩亚洲蝌蚪91| 精品国产乱子伦一区二区三区| 久久久精品国产亚洲av高清涩受| 十分钟在线观看高清视频www| 正在播放国产对白刺激| 中文字幕精品免费在线观看视频| 国产精品1区2区在线观看. | 亚洲成国产人片在线观看| 中文字幕最新亚洲高清| 亚洲自偷自拍图片 自拍| 黑人巨大精品欧美一区二区mp4| 国产精品影院久久| av网站在线播放免费| 亚洲一区二区三区欧美精品| e午夜精品久久久久久久| 欧美久久黑人一区二区| www.熟女人妻精品国产| 丝袜美足系列| 免费日韩欧美在线观看| 久久国产精品影院| 国产一区有黄有色的免费视频| 久久国产精品大桥未久av| 日韩欧美免费精品| 亚洲人成电影免费在线| 成人18禁在线播放| 一区二区三区乱码不卡18| 999精品在线视频| 精品国产一区二区三区四区第35| www.999成人在线观看| 精品福利永久在线观看| 日韩有码中文字幕| 亚洲精品国产色婷婷电影| 亚洲精品av麻豆狂野| 一边摸一边做爽爽视频免费| 在线观看一区二区三区激情| 日韩欧美三级三区| 亚洲av日韩精品久久久久久密| av网站免费在线观看视频| 视频区图区小说| 久久国产精品人妻蜜桃| a级毛片黄视频| 久久精品国产a三级三级三级| 精品乱码久久久久久99久播| 69av精品久久久久久 | 热99re8久久精品国产| svipshipincom国产片| 搡老熟女国产l中国老女人| 欧美 日韩 精品 国产| 日韩欧美一区视频在线观看| 成人免费观看视频高清| 最近最新免费中文字幕在线| 国产黄色免费在线视频| 欧美黄色淫秽网站| 我的亚洲天堂| 亚洲色图综合在线观看| 亚洲欧美激情在线| 国产高清激情床上av| 亚洲成人国产一区在线观看| 中文字幕另类日韩欧美亚洲嫩草| 99精品久久久久人妻精品| 久久午夜综合久久蜜桃| 欧美国产精品va在线观看不卡| 国产精品久久久久久精品古装| 精品乱码久久久久久99久播| 中亚洲国语对白在线视频| 日韩中文字幕视频在线看片| 精品国产乱子伦一区二区三区| 黄色成人免费大全| 12—13女人毛片做爰片一| 1024香蕉在线观看| av网站免费在线观看视频| 久久久久久亚洲精品国产蜜桃av| 窝窝影院91人妻| 国产精品1区2区在线观看. | 国产国语露脸激情在线看| 国产精品电影一区二区三区 | 老司机靠b影院| 久久精品国产亚洲av香蕉五月 | 手机成人av网站| 新久久久久国产一级毛片| 女人久久www免费人成看片| 制服诱惑二区| 国产aⅴ精品一区二区三区波| 黄色成人免费大全| 国产日韩欧美在线精品| 国产一卡二卡三卡精品| 久久久久久久久久久久大奶| aaaaa片日本免费| 色视频在线一区二区三区| 一二三四在线观看免费中文在| 午夜精品国产一区二区电影| 欧美激情 高清一区二区三区| 成在线人永久免费视频| 男人操女人黄网站| 男女高潮啪啪啪动态图| 一区二区三区国产精品乱码| 日本黄色视频三级网站网址 | 老司机午夜十八禁免费视频| 大片免费播放器 马上看| 亚洲国产成人一精品久久久| 18禁国产床啪视频网站| 一级片免费观看大全| av天堂久久9| 水蜜桃什么品种好| 天天影视国产精品| 国产精品自产拍在线观看55亚洲 | 嫩草影视91久久| 亚洲精品在线观看二区| 日韩 欧美 亚洲 中文字幕| 肉色欧美久久久久久久蜜桃| 精品熟女少妇八av免费久了| 国产片内射在线| 中文字幕精品免费在线观看视频| 人人妻,人人澡人人爽秒播| 精品一品国产午夜福利视频| 大型av网站在线播放| 最新在线观看一区二区三区| netflix在线观看网站| 欧美日韩黄片免| 国产欧美日韩一区二区三| 亚洲黑人精品在线| 中文字幕色久视频| 黄色怎么调成土黄色| 免费在线观看影片大全网站| 国产一区二区三区综合在线观看| 十八禁人妻一区二区| 最近最新中文字幕大全电影3 | 香蕉丝袜av| 亚洲国产av新网站| 国产成+人综合+亚洲专区| 久久久精品区二区三区| 亚洲精品美女久久久久99蜜臀| 狂野欧美激情性xxxx| 麻豆成人av在线观看| 亚洲精品粉嫩美女一区| 亚洲成国产人片在线观看| 九色亚洲精品在线播放| www.999成人在线观看| 国产成人免费观看mmmm| 黄色成人免费大全| 国产在线免费精品| 成人三级做爰电影| 动漫黄色视频在线观看| 91成年电影在线观看| 法律面前人人平等表现在哪些方面| 999久久久精品免费观看国产| 久久国产精品大桥未久av| 亚洲精品国产一区二区精华液| 久久婷婷成人综合色麻豆| 日韩视频一区二区在线观看| 高清av免费在线| 欧美亚洲日本最大视频资源| 免费女性裸体啪啪无遮挡网站| 国产aⅴ精品一区二区三区波| 日本五十路高清| 男女边摸边吃奶| 亚洲av成人一区二区三| 另类亚洲欧美激情| 亚洲国产av影院在线观看| 亚洲专区字幕在线| 91精品三级在线观看| 婷婷丁香在线五月| 国产一区二区三区视频了| 国产亚洲精品久久久久5区| 亚洲 欧美一区二区三区| 成年动漫av网址| 波多野结衣一区麻豆| 曰老女人黄片| 精品人妻熟女毛片av久久网站| 久久久国产一区二区| 老熟妇仑乱视频hdxx| 久久免费观看电影| 精品国产国语对白av| 宅男免费午夜| 久久久久国产一级毛片高清牌| 飞空精品影院首页| 亚洲国产欧美一区二区综合| 久久99一区二区三区| 中文字幕色久视频| 每晚都被弄得嗷嗷叫到高潮| 欧美精品一区二区大全| 亚洲国产欧美网| a级片在线免费高清观看视频| av片东京热男人的天堂| 老熟妇乱子伦视频在线观看| 色94色欧美一区二区| 黑人欧美特级aaaaaa片| 欧美av亚洲av综合av国产av| 久9热在线精品视频| 国产精品免费视频内射| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久久免费视频了| 无遮挡黄片免费观看| 中文字幕人妻熟女乱码| 一二三四在线观看免费中文在| 国内毛片毛片毛片毛片毛片| 人人妻人人添人人爽欧美一区卜| 国产不卡一卡二| 日本撒尿小便嘘嘘汇集6| 欧美成人午夜精品| 香蕉丝袜av| 丁香六月欧美| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 18禁观看日本| av超薄肉色丝袜交足视频| 亚洲欧美一区二区三区久久| 久久久久精品人妻al黑| 午夜久久久在线观看| 国产高清视频在线播放一区| 中亚洲国语对白在线视频| 高清欧美精品videossex| 日本黄色日本黄色录像| 高清在线国产一区| 国产一区有黄有色的免费视频| 国产精品美女特级片免费视频播放器 | 最近最新免费中文字幕在线| 热re99久久国产66热| 99精品欧美一区二区三区四区| 欧美国产精品一级二级三级| 国产深夜福利视频在线观看| tocl精华| 成人特级黄色片久久久久久久 | 在线播放国产精品三级| 91成人精品电影| 国产麻豆69| 在线观看www视频免费| 极品教师在线免费播放| 欧美久久黑人一区二区| 精品视频人人做人人爽| 免费不卡黄色视频| 不卡一级毛片| 欧美亚洲 丝袜 人妻 在线| 欧美性长视频在线观看| av视频免费观看在线观看| 99热网站在线观看| 少妇 在线观看| 人成视频在线观看免费观看| 99九九在线精品视频| 91精品三级在线观看| 丰满迷人的少妇在线观看| 久久亚洲精品不卡| 亚洲免费av在线视频| 黑人巨大精品欧美一区二区mp4| 亚洲av日韩在线播放| 美国免费a级毛片| 女人久久www免费人成看片| 中文字幕最新亚洲高清| 80岁老熟妇乱子伦牲交| 一级片'在线观看视频| 一个人免费看片子| 国产成人av教育| 久久婷婷成人综合色麻豆| 汤姆久久久久久久影院中文字幕| 色婷婷久久久亚洲欧美| 露出奶头的视频| 国产免费视频播放在线视频| 国产一区二区三区综合在线观看| 一进一出好大好爽视频| 大片电影免费在线观看免费| 日本欧美视频一区| 久久久久国产一级毛片高清牌| 亚洲美女黄片视频| 色婷婷av一区二区三区视频| 久久久久久久国产电影| 香蕉国产在线看| 亚洲精品国产区一区二| 久久午夜亚洲精品久久| 亚洲欧美精品综合一区二区三区| 欧美精品人与动牲交sv欧美| 亚洲午夜精品一区,二区,三区| 高清黄色对白视频在线免费看| 久久精品亚洲熟妇少妇任你| 黄频高清免费视频| 色播在线永久视频| 俄罗斯特黄特色一大片| 91九色精品人成在线观看| 99精国产麻豆久久婷婷| 黑人操中国人逼视频| 久久国产精品大桥未久av| 精品少妇内射三级| 大陆偷拍与自拍| 丝袜美足系列| 天堂中文最新版在线下载| 黄网站色视频无遮挡免费观看| av片东京热男人的天堂| 午夜福利,免费看| 亚洲专区字幕在线| 欧美成狂野欧美在线观看| 日本五十路高清| 国产aⅴ精品一区二区三区波| 99国产综合亚洲精品| 黄色视频在线播放观看不卡| 一区二区av电影网| 考比视频在线观看| 欧美日韩国产mv在线观看视频| 色精品久久人妻99蜜桃| 欧美 日韩 精品 国产| 亚洲国产欧美在线一区| 岛国毛片在线播放| 国产精品1区2区在线观看. | 十分钟在线观看高清视频www| 国产伦理片在线播放av一区| 欧美人与性动交α欧美软件| 汤姆久久久久久久影院中文字幕| 无遮挡黄片免费观看| 午夜免费鲁丝| 大型av网站在线播放| 999久久久国产精品视频| 黑人猛操日本美女一级片| 国产亚洲欧美在线一区二区| 久久国产精品影院| 久久久久久亚洲精品国产蜜桃av| a级片在线免费高清观看视频| 亚洲专区字幕在线| 欧美人与性动交α欧美精品济南到| 男女无遮挡免费网站观看| 啦啦啦视频在线资源免费观看| 中亚洲国语对白在线视频| 精品人妻1区二区| 国产精品一区二区免费欧美| 男女免费视频国产| 大型黄色视频在线免费观看| 成人18禁在线播放| 欧美日韩亚洲综合一区二区三区_| 动漫黄色视频在线观看| 人妻久久中文字幕网| 性高湖久久久久久久久免费观看| 在线观看免费视频日本深夜| 亚洲精品乱久久久久久| 19禁男女啪啪无遮挡网站| 亚洲精品一卡2卡三卡4卡5卡| 操美女的视频在线观看| 99re6热这里在线精品视频| 久久久久网色| 国产麻豆69| 人人妻人人爽人人添夜夜欢视频| 在线 av 中文字幕| 久久久精品国产亚洲av高清涩受| 欧美激情 高清一区二区三区| 热99国产精品久久久久久7| 国产欧美日韩综合在线一区二区| 久久天堂一区二区三区四区| 老司机深夜福利视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 80岁老熟妇乱子伦牲交| 国产在视频线精品| 欧美日韩一级在线毛片| 国产不卡av网站在线观看| 99国产精品一区二区三区| bbb黄色大片| 黄片小视频在线播放| 国产老妇伦熟女老妇高清| 国产精品美女特级片免费视频播放器 | 久久久精品国产亚洲av高清涩受| 国产精品免费一区二区三区在线 | 欧美日韩精品网址| 日韩有码中文字幕| 久久中文看片网| 两性午夜刺激爽爽歪歪视频在线观看 | 美女午夜性视频免费| 精品卡一卡二卡四卡免费| 亚洲五月婷婷丁香| 一区福利在线观看|