• <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)分析鑒定法系大白豬繁殖性狀基因組區(qū)域和候選基因

    2022-08-13 02:57:16李秀領(lǐng)陳健梅吳姿儀韓雪蕾喬瑞敏王克君王獻(xiàn)偉李新建
    中國畜牧雜志 2022年8期
    關(guān)鍵詞:遺傳變異遺傳力表型

    李秀領(lǐng),陳健梅,吳姿儀,韓雪蕾,喬瑞敏,王克君,楊 峰,王獻(xiàn)偉,李新建*

    (1.河南農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院,河南鄭州 450046;2.河南省畜牧總站,河南鄭州 450008)

    母豬的繁殖能力在養(yǎng)豬業(yè)的經(jīng)濟(jì)效益構(gòu)成中起著重要作用。為了獲取更高的收益,提高母豬的繁殖性能尤為重要??偖a(chǎn)仔數(shù)(Total Number Born,TNB)、產(chǎn)活仔數(shù)(Number Born Alive,NBA)和出生窩重(Litter Weight Born Alive,LWB)是評(píng)估母豬繁殖性能的關(guān)鍵指標(biāo)。然而,繁殖性狀的遺傳結(jié)構(gòu)較為復(fù)雜,由于其低遺傳力特性(約0.1),傳統(tǒng)的育種計(jì)劃獲得的遺傳改良非常緩慢。因此,標(biāo)記輔助選擇(Marker-Assisted Selection,MAS)和基因組選擇(Genomic Selection,GS)被廣泛應(yīng)用于繁殖性狀的遺傳改良。

    近幾十年來,隨著分子育種技術(shù)的發(fā)展,全基因組關(guān)聯(lián)分析(Genome-Wide Association Study,GWAS)已被廣泛用于識(shí)別復(fù)雜性狀中的數(shù)量性狀基因座(Quan titative Trait Loci,QTL)和數(shù)量性狀核苷酸(Quantita tive Trait Nucleotides,QTN)。在豬相關(guān)的研究中,GWAS 已經(jīng)成為定位與復(fù)雜性狀相關(guān)的SNPs 和候選基因的有力方法之一。在GWAS 的統(tǒng)計(jì)方法中,單步全基因組關(guān)聯(lián)分析(Single-Step GWAS,ssGWAS)已被用于研究豬的經(jīng)濟(jì)性狀。與傳統(tǒng)GWAS 相比,ssGWAS 可以將個(gè)體的基因型、表型和系譜信息都整合進(jìn)混合線性模型中。因此,ssGWAS 可以提高低遺傳力性狀候選基因定位的精度。

    此外,繁殖性狀受許多微效基因控制,但到目前為止,只有少數(shù)候選基因被確定,而且只是解釋了一小部分遺傳變異。因此,繁殖性狀的具體遺傳機(jī)制仍然需要進(jìn)一步研究。本研究的目的是通過在法系大白豬群體中進(jìn)行ssGWAS 檢測影響繁殖性狀(TNB、NBA 和LWB)的遺傳變異并確定候選基因。

    1 材料與方法

    1.1 動(dòng)物和表型數(shù)據(jù) 在本研究中,收集了河南某核心養(yǎng)豬場392 頭法系大白豬的系譜信息以及2020 年1—10 月的表型信息。其中,表型包括TNB、NBA 和LWB,所有表型記錄均為母豬第1 胎次產(chǎn)仔記錄。所有個(gè)體均被分組在標(biāo)準(zhǔn)的商業(yè)豬舍中并且每個(gè)個(gè)體都是在相同的飼料和飼養(yǎng)管理?xiàng)l件下飼養(yǎng)的。TNB、NBA和LWB 的表型記錄的平均值、標(biāo)準(zhǔn)差、最小值和最大值見表1。

    表1 繁殖性狀的描述性統(tǒng)計(jì)

    1.2 基因型數(shù)據(jù)和質(zhì)控 使用組織基因組DNA 提取試劑盒從豬耳組織樣中提取DNA,并使用紐勤公司GeneSeek Genomic Profiler(GGP)豬50K SNP 芯 片進(jìn)行基因分型。數(shù)據(jù)質(zhì)量控制由PLINK v1.90進(jìn)行。剔除基因型檢出率<0.90、次等位基因頻率<0.05 和Hardy-Weinberg 平衡檢驗(yàn)<10的SNP 標(biāo)記,刪 除SNP 標(biāo)記檢出率<0.90 的個(gè)體。所有位于未映射區(qū)域的SNPs 也被排除。最終保留了43 098 個(gè)SNPs 和392 個(gè)個(gè)體用于后續(xù)分析。

    1.3 遺傳參數(shù)估計(jì) 使用GCTA v1.93.2進(jìn)行限制性最大似然算法來估計(jì)所有性狀的遺傳參數(shù)。所有的常染色體SNPs 被用來計(jì)算基因組關(guān)系矩陣(GRM),然后計(jì)算主成分。用單變量動(dòng)物模型來估計(jì)遺傳力,用雙變量GREML 分析來估計(jì)性狀間的遺傳相關(guān)。統(tǒng)計(jì)模型為:

    其中,是表型觀測值的向量;代表固定效應(yīng),包括前3 個(gè)主成分、產(chǎn)仔季節(jié)和妊娠天數(shù),是的相關(guān)矩陣;是隨機(jī)多基因效應(yīng)的向量;是殘差效應(yīng)的向量。

    1.4 ssGWAS GWAS 是使用BLUPF90 軟件家族進(jìn)行的。RENUMF90 模塊用于準(zhǔn)備表型和系譜數(shù)據(jù),BLUPF90 模塊用于計(jì)算基因組估計(jì)育種值(Genomic Estimated Breeding Values,GEBV),POSTGSF90 模塊用于將GEBVs 轉(zhuǎn)換為SNPs 效應(yīng),并計(jì)算相鄰窗口所解釋的方差。

    在本研究中,單性狀動(dòng)物模型被用于GWAS 分析。統(tǒng)計(jì)學(xué)模型可以寫成:

    其中,是基于所有動(dòng)物系譜的分子關(guān)系矩陣,是基因分型個(gè)體的分子關(guān)系矩陣的子集,是基因組關(guān)系矩陣。上標(biāo)-1 表示它是相關(guān)矩陣的逆數(shù)。VanRaden(2008)構(gòu)建G 矩陣的方法如下:

    其中,Z 是根據(jù)等位基因頻率編碼的標(biāo)記矩陣,D 是SNP 效應(yīng)方差的權(quán)重對(duì)角矩陣(最初D=I)。ssGWAS的計(jì)算方法如下:

    ①構(gòu)建G 矩陣;

    ②計(jì)算所有動(dòng)物的GEBVs;

    ⑤歸一化SNP 權(quán)重以保持總方差不變;

    ⑥根據(jù)SNP 權(quán)重重新計(jì)算矩陣;=';

    ⑦循環(huán)步驟②。

    從第②步到第⑦步重復(fù)一次迭代過程,按照Wang等的描述,計(jì)算第i 個(gè)連續(xù)SNP(SNP 窗口)解釋的遺傳變異百分比:

    Beissinger 等報(bào)告說,與其他較大的窗口尺寸相比,5 個(gè)或10 個(gè)SNP 組成的滑動(dòng)窗口的檢測率與假陽性率都是最有利的。因此,本研究計(jì)算了由10 個(gè)相鄰SNP組成的滑動(dòng)窗口所解釋的加性遺傳變異的比例。

    1.5 候選基因組區(qū)域注釋 解釋的遺傳變異大于1%的基因組窗口被認(rèn)為是與性狀相關(guān)的。在11.1基因組(http://www.ensembl.org/)上搜索了候選窗口內(nèi)的基因或部分重疊的基因。使用大鼠基因組數(shù)據(jù)庫(RGD,https://rgd.mcw.edu/rgdweb/homepage/)進(jìn)行人類表型本體論(Human Phenotype Ontology,HPO)分析,選擇Bonferroni 校正后<0.05 的條目進(jìn)行后續(xù)分析。隨后,使用KOBAS 在線網(wǎng)站對(duì)與生殖系統(tǒng)相關(guān)條目中涉及的候選基因進(jìn)行分析,以確定基因參與的KEGG 信號(hào)通路。<0.05 的KEGG 通路被選中。使用PigQTLdb 數(shù)據(jù)庫(http://www.animalgenome.org/cgi-bin/QTLdb/SS/index)對(duì)豬基因組中先前確定的QTL 進(jìn)行評(píng)估。

    2 結(jié)果

    2.1 表型統(tǒng)計(jì) 表1 列出了3 個(gè)繁殖性狀的描述性統(tǒng)計(jì)。其中,TNB 和NBA 的范圍是0~24,平均值分別為14.91 和13.11。LWB 平均值為15.14。此外,所有性狀的變異系數(shù)值都較高,這表明數(shù)據(jù)在某種程度上是離散的。另外,3 個(gè)性狀之間高度正相關(guān),表型相關(guān)系數(shù)分別為0.86、0.67 和0.81(表2)。

    2.2 遺傳力和遺傳相關(guān) 本研究對(duì)3 個(gè)繁殖性狀的遺傳力和遺傳相關(guān)都進(jìn)行了估計(jì)(表2)。3 個(gè)性狀的估計(jì)遺傳力在0.06~0.12 之間,其中,NBA 遺傳力最高,為0.12,TNB 遺傳力最低,為0.06。3 個(gè)性狀兩兩間的遺傳相關(guān)從0.62 到0.83,都呈現(xiàn)高度正相關(guān)。

    表2 繁殖性狀的遺傳力(對(duì)角線,粗體)、遺傳相關(guān)(上三角)和表型相關(guān)系數(shù)(下三角)

    2.3 GWAS 本研究將連續(xù)10 個(gè)SNPs 劃分為一個(gè)窗口后,共有43 098 個(gè)窗口用于后續(xù)的GWAS 分析,其中最后9 個(gè)窗口包含的SNPs 少于10 個(gè)。

    本研究中3 個(gè)性狀分別鑒定到95、72、109 個(gè)相關(guān)的窗口,它們分別解釋了TNB、NBA 和LWB 總遺傳變異的1% 或以上(圖1)。共有40 個(gè)基因組區(qū)域與繁殖性狀相關(guān),其中TNB、NBA 和LWB 分別包括15、14、11 個(gè)區(qū)域。對(duì)于TNB,共有15 個(gè)基因組區(qū)域,分布在染色體(SSC)1、3、5、7、13、15、16、17 和X 上,其中最高的區(qū)域解釋了4.794%的遺傳變異。有14 個(gè)基因組區(qū)域與NBA 顯著相關(guān),它們分布在SSC 1、2、3、6、7、16、17 和X 上,最高的區(qū)域解釋了2.214%的遺傳變異。對(duì)于LWB,109 個(gè)窗口串聯(lián)成11 個(gè)基因組區(qū)域,這些基因組區(qū)域分布在SSC 2、3、6、10、13、16 和X 上,它們所解釋的最大遺傳變異為2.364%(表3)。

    表3 繁殖性狀相關(guān)基因組區(qū)域

    圖1 繁殖性狀ssGWAS 曼哈頓圖

    2.4 候選基因鑒定 經(jīng)過注釋,共得到532 個(gè)基因,包括381 個(gè)編碼蛋白的基因。隨后,381 個(gè)基因被用來進(jìn)行HPO 分析,結(jié)果富集在以下3 個(gè)條目:泌尿生殖系統(tǒng)表型、心血管系統(tǒng)表型和神經(jīng)系統(tǒng)表型(圖2A)。在泌尿生殖系統(tǒng)表型中,共得到81 個(gè)候選基因(圖2B,2C)。有趣的是,大多數(shù)候選基因都位于2 號(hào)染色體上。隨后,對(duì)候選基因進(jìn)行了進(jìn)一步的KEGG 分析,得到了82 條顯著的通路。在82 條顯著的通路中共發(fā)現(xiàn)7 條與生殖相關(guān)的KEGG 通路,包括5 條與生殖激素相關(guān)的信號(hào)通路和2 條卵細(xì)胞減數(shù)分裂及成熟相關(guān)的通路(表4)。對(duì)于所有的候選基因,本研究查詢了以前的研究是否報(bào)道過與生殖性狀有關(guān)的文獻(xiàn),結(jié)果列于表5。

    表4 候選基因參與的顯著的繁殖相關(guān)KEGG 通路

    表5 候選基因中先前報(bào)道過的與繁殖相關(guān)的基因

    圖2 候選基因鑒定

    3 討 論

    產(chǎn)仔性狀包含仔豬的TNB、NBA 和LWB,與母豬的繁殖性能密切相關(guān),影響豬的生產(chǎn)效率和經(jīng)濟(jì)利益。因此,探索豬產(chǎn)仔性狀的遺傳結(jié)構(gòu)是極為重要的。此外,探索產(chǎn)仔性狀的遺傳機(jī)制非常復(fù)雜,整個(gè)繁殖過程受多種激素和生理途徑的調(diào)節(jié)。GWAS 是檢測與目標(biāo)性狀相關(guān)的候選基因最有效方法之一。本研究利用ssGWAS方法探索了392 頭大白豬的TNB、NBA 和LWB 的遺傳變異,通過HPO 和KEGG 分析進(jìn)一步挖掘了與生殖相關(guān)的基因以及這些基因涉及的通路。

    繁殖性狀的遺傳力通常在0.1 左右。在本研究中,3個(gè)繁殖性狀的估計(jì)遺傳力為0.06~0.12,與前人研究結(jié)果一致。此外,3 個(gè)性狀之間的遺傳和表型相關(guān)都呈現(xiàn)出高度正相關(guān),這些結(jié)果與Wolf 等的研究一致。

    本研究發(fā)現(xiàn)了40 個(gè)解釋1%以上遺傳變異的基因組區(qū)域。在這些區(qū)域中,還發(fā)現(xiàn)了一些以前報(bào)道的與繁殖有關(guān)的QTL。在3 個(gè)性狀中SSC13 上的48.65~52.85 Mb 區(qū)域包含了最多解釋超過1% 的遺傳變異的窗口,并且在這個(gè)基因組區(qū)域中有9 個(gè)先前報(bào)道的相關(guān)QTLs。這些QTLs 包括3 個(gè)黃體數(shù)量QTLs(493、8 826、24 285),2 個(gè)卵巢重量QTLs(18 354、18 356),1 個(gè)青春期年齡QTLs(8 824),1 個(gè)子宮角重量QTLs(18 409),1 個(gè)總產(chǎn)仔數(shù)QTLs(15 1333)和1 個(gè)死胎數(shù)QTLs(57 525)。此外,SSC14 的104.44~106.90 Mb 區(qū)域解釋了2.87%的TNB 性狀的遺傳變異,并且在此區(qū)域注釋到了基因。據(jù)報(bào)道,基因在胚胎發(fā)育中起著不可或缺的作用。

    基于關(guān)聯(lián)分析鑒定到的顯著基因組區(qū)域內(nèi)的蛋白編碼基因,本研究進(jìn)行了HPO 分析,發(fā)現(xiàn)這些基因聚成了3 類:泌尿生殖系統(tǒng)、心血管系統(tǒng)和神經(jīng)系統(tǒng)。共有7 種表型涉及泌尿生殖系統(tǒng),包括原發(fā)性閉經(jīng)、閉經(jīng)、月經(jīng)周期異常、女性生殖系統(tǒng)生理功能異常、泌尿生殖系統(tǒng)異常、生殖系統(tǒng)異常以及生殖系統(tǒng)生理學(xué)異常。這7 個(gè)表型共包含81 個(gè)候選基因,且大多聚集在SSC 2 上。經(jīng)過搜索文獻(xiàn),在SSC2 上獲得了10 個(gè)先前報(bào)道的與生殖有關(guān)的基因。其中,畸形表皮自動(dòng)調(diào)節(jié)因子1(DEAF1)編碼一個(gè)多域轉(zhuǎn)錄因子。相關(guān)結(jié)合位點(diǎn)的變化可能與山羊的產(chǎn)仔數(shù)有關(guān)?;騾⑴c磷酸戊糖途徑(PPP)的第2 階段,并在糖酵解和PPP 之間提供一個(gè)鏈接。缺陷的患者在新生兒期出現(xiàn)低出生體重,甚至在嬰兒期早期死亡。胰島素樣生長因子2()基因是一個(gè)父系印記基因,已知它能增加豬的肌肉生長并減少脂肪沉積。然而,也有報(bào)道稱對(duì)母豬的繁殖力有影響。此外,基因已被確認(rèn)為人類、小鼠和牛的一個(gè)印記基因。印記基因與早期胚胎和胎盤的發(fā)育有關(guān)系。這些基因的表達(dá)水平的變化可以導(dǎo)致異常的生長和胚胎死亡。據(jù)報(bào)道,在豬的心臟、卵巢等部位有母體表達(dá),在肺和腎有雙側(cè)表達(dá)。因此,基因在豬身上的印記還有待驗(yàn)證?;蚴莝irtuin家族的一個(gè)重要成員,在調(diào)節(jié)細(xì)胞生理和衰老方面發(fā)揮了廣泛的作用。Gorczyca 等發(fā)現(xiàn)基因在豬卵巢卵泡細(xì)胞中表達(dá)。

    為進(jìn)一步了解候選基因的功能和這些基因之間的聯(lián)系,本研究還進(jìn)行了KEGG 分析,在所有顯著的KEGG 信號(hào)通路中得到了7 條與繁殖性狀相關(guān)的通路。這7 條通路分別是催產(chǎn)素信號(hào)通路、催乳素信號(hào)通路、GnRH 信號(hào)通路、松弛素信號(hào)通路、雌激素信號(hào)通路、母細(xì)胞減數(shù)分裂和孕激素介導(dǎo)卵細(xì)胞成熟通路。其中催產(chǎn)素、催乳素、促性腺激素釋放激素(GnRH)、松弛素和雌激素都屬于生殖激素,與生殖活動(dòng)如排卵、懷孕、胚胎發(fā)育、分娩和泌乳直接相關(guān)。參與生殖激素相關(guān)途徑的共有7 個(gè)基因,即和CREB3L3。其中,MAP2K2編碼的蛋白質(zhì)屬于MAPK 激酶家族,基因?qū)μケP形成期間合體滋養(yǎng)層的正常發(fā)育有貢獻(xiàn)?;?qū)儆赗as 腫瘤基因家族,其編碼的產(chǎn)物在信號(hào)轉(zhuǎn)導(dǎo)途徑中發(fā)揮作用。有報(bào)道稱miR-210 通過和調(diào)節(jié)水牛排卵前的顆粒細(xì)胞功能?;騾⑴c調(diào)控血管生長和基礎(chǔ)血壓。Zou 等研究表明基因在胎盤血管的形成中有重要作用?;蚓幋aAMP 依賴性轉(zhuǎn)錄因子家族,它可以調(diào)節(jié)脂質(zhì)代謝并介導(dǎo)對(duì)感染的反應(yīng)。相關(guān)研究表明通過控制GH 信號(hào)傳導(dǎo),對(duì)身體的生長有負(fù)向調(diào)節(jié)作用。此外,有2 個(gè)基因參與卵母細(xì)胞減數(shù)分裂和成熟相關(guān)通路,分別是和。Aurora-A(AURKA)是一種絲氨酸/蘇氨酸蛋白激酶。Yao 等研究結(jié)果顯示基因在豬卵母細(xì)胞減數(shù)分裂成熟、受精和早期胚胎有絲分裂過程中對(duì)微管的組裝起著關(guān)鍵的調(diào)節(jié)作用?;蚓幋a胰島素,這是一種多肽激素,可調(diào)節(jié)生理學(xué)的各個(gè)方面,是決定后生動(dòng)物雌性生殖的關(guān)鍵因素。據(jù)報(bào)道,胰島素參與了奶牛的卵泡成熟。

    4 結(jié) 論

    本研究使用50K SNP 芯片對(duì)法系大白豬群體中的TNB、NBA 和LWB 表型進(jìn)行了ssGWAS 分析,鑒定到了40 個(gè)基因組區(qū)域。區(qū)域內(nèi)注釋到的所有編碼蛋白的基因被用于進(jìn)行HPO 分析,發(fā)現(xiàn)了81 個(gè)與繁殖有關(guān)的候選基因,這些基因主要涉及GnRH 信號(hào)通路、雌性激素信號(hào)通路和卵細(xì)胞成熟通路等。其中,本研究還發(fā)現(xiàn)一些與產(chǎn)仔數(shù)有關(guān)的基因,包括、、、、等。盡管篩選出了一些與繁殖相關(guān)的候選基因,但部分基因在畜禽中還沒有被研究過,因此這些候選基因的功能還需要進(jìn)一步探討。下一步,將利用高密度測序數(shù)據(jù)對(duì)這些區(qū)域進(jìn)行精細(xì)定位,以進(jìn)一步確定與產(chǎn)仔數(shù)相關(guān)的基因并篩選出致因突變,為豬的遺傳機(jī)制分析和育種提供參考。

    猜你喜歡
    遺傳變異遺傳力表型
    某國家生豬核心育種場加系大白母豬繁殖性狀遺傳參數(shù)估計(jì)
    美系大白種豬生長性狀和繁殖性狀遺傳力估計(jì)
    養(yǎng)豬(2022年4期)2022-08-17 07:07:02
    先導(dǎo)編輯技術(shù)可編輯近90%的人類遺傳變異
    基于改進(jìn)遺傳變異算子的海島算法
    電子制作(2019年24期)2019-02-23 13:22:18
    建蘭、寒蘭花表型分析
    火力楠子代遺傳變異分析及優(yōu)良家系選擇
    通過胎盤效率改良母豬繁殖性能的研究進(jìn)展
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    慢性乙型肝炎患者HBV基因表型與血清學(xué)測定的臨床意義
    72例老年急性白血病免疫表型分析
    精品国产美女av久久久久小说| 校园春色视频在线观看| 男女高潮啪啪啪动态图| 精品久久久精品久久久| 嫩草影视91久久| 最好的美女福利视频网| 国产免费现黄频在线看| 亚洲精品中文字幕一二三四区| 国产欧美日韩一区二区精品| 免费av中文字幕在线| 成熟少妇高潮喷水视频| 丝袜美足系列| 日本免费一区二区三区高清不卡 | 日韩视频一区二区在线观看| 欧美日韩乱码在线| 嫁个100分男人电影在线观看| 国产精品久久久久久人妻精品电影| 一级片免费观看大全| 亚洲国产看品久久| 色哟哟哟哟哟哟| 国产精品亚洲一级av第二区| 午夜福利影视在线免费观看| 高清欧美精品videossex| 三级毛片av免费| 这个男人来自地球电影免费观看| 性色av乱码一区二区三区2| 久久天堂一区二区三区四区| 精品少妇一区二区三区视频日本电影| 欧美乱码精品一区二区三区| 三级毛片av免费| 免费日韩欧美在线观看| 老司机亚洲免费影院| 一边摸一边做爽爽视频免费| 久久久国产欧美日韩av| 欧美老熟妇乱子伦牲交| 成人免费观看视频高清| 国产欧美日韩综合在线一区二区| 正在播放国产对白刺激| 黄色毛片三级朝国网站| 欧美精品亚洲一区二区| 亚洲美女黄片视频| 亚洲精品粉嫩美女一区| 啪啪无遮挡十八禁网站| 久久影院123| 精品久久久精品久久久| 欧美激情久久久久久爽电影 | 色老头精品视频在线观看| 不卡一级毛片| 国产不卡一卡二| www.999成人在线观看| 久久中文字幕人妻熟女| 欧美一级毛片孕妇| 中文字幕最新亚洲高清| 色老头精品视频在线观看| 欧美性长视频在线观看| 悠悠久久av| 在线观看一区二区三区激情| 国产免费现黄频在线看| 男女床上黄色一级片免费看| av超薄肉色丝袜交足视频| 欧美日韩亚洲国产一区二区在线观看| 一二三四在线观看免费中文在| 日韩欧美一区视频在线观看| 夜夜爽天天搞| 91老司机精品| 身体一侧抽搐| 亚洲久久久国产精品| 大香蕉久久成人网| 亚洲色图综合在线观看| 激情视频va一区二区三区| 国产成人系列免费观看| 日韩精品免费视频一区二区三区| 又大又爽又粗| 欧美最黄视频在线播放免费 | 日韩成人在线观看一区二区三区| 亚洲九九香蕉| tocl精华| 亚洲欧洲精品一区二区精品久久久| 黄色怎么调成土黄色| 精品国产乱子伦一区二区三区| 国产亚洲av高清不卡| 国产成人精品久久二区二区免费| 欧美日本亚洲视频在线播放| 国产三级黄色录像| 久久久精品国产亚洲av高清涩受| 在线观看舔阴道视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美日韩无卡精品| 12—13女人毛片做爰片一| 久热这里只有精品99| cao死你这个sao货| 成人18禁高潮啪啪吃奶动态图| 又黄又粗又硬又大视频| 亚洲专区字幕在线| 男人舔女人的私密视频| 中文字幕高清在线视频| 亚洲欧美一区二区三区久久| av中文乱码字幕在线| 精品国产美女av久久久久小说| 欧美日韩乱码在线| 久久这里只有精品19| 国产在线精品亚洲第一网站| 亚洲午夜精品一区,二区,三区| 成人国产一区最新在线观看| 免费高清在线观看日韩| 国产亚洲欧美精品永久| 久久香蕉精品热| 国产成年人精品一区二区 | 美女福利国产在线| 亚洲色图av天堂| 国产一区二区三区视频了| av欧美777| 男女下面插进去视频免费观看| e午夜精品久久久久久久| 国产日韩一区二区三区精品不卡| 国产亚洲av高清不卡| 美女午夜性视频免费| 母亲3免费完整高清在线观看| 国产极品粉嫩免费观看在线| 欧美日韩亚洲国产一区二区在线观看| 欧美性长视频在线观看| 国产又色又爽无遮挡免费看| 国产精品香港三级国产av潘金莲| 18美女黄网站色大片免费观看| 性欧美人与动物交配| 亚洲一区中文字幕在线| 亚洲久久久国产精品| 中国美女看黄片| 精品一区二区三区av网在线观看| 国产成人精品久久二区二区免费| 国产亚洲精品一区二区www| 老司机在亚洲福利影院| 精品少妇一区二区三区视频日本电影| 不卡一级毛片| 伊人久久大香线蕉亚洲五| 久久国产精品男人的天堂亚洲| 免费少妇av软件| 国产91精品成人一区二区三区| 丰满的人妻完整版| 久久国产精品影院| 国产精品98久久久久久宅男小说| 丝袜在线中文字幕| 午夜免费成人在线视频| 亚洲成人久久性| 午夜福利欧美成人| 亚洲自拍偷在线| 99re在线观看精品视频| 啦啦啦 在线观看视频| 成人三级黄色视频| 一本综合久久免费| 怎么达到女性高潮| 男人舔女人的私密视频| 黑丝袜美女国产一区| 男人操女人黄网站| 香蕉国产在线看| 精品久久久久久成人av| 亚洲专区中文字幕在线| 亚洲中文av在线| 久久精品91无色码中文字幕| 女人精品久久久久毛片| 无限看片的www在线观看| 一级a爱片免费观看的视频| 黄色视频,在线免费观看| 国产黄色免费在线视频| 97超级碰碰碰精品色视频在线观看| 国产精品电影一区二区三区| 国产欧美日韩一区二区精品| 最近最新中文字幕大全电影3 | 大型黄色视频在线免费观看| 亚洲av日韩精品久久久久久密| 69精品国产乱码久久久| 精品久久久久久久久久免费视频 | 精品国产亚洲在线| 在线观看一区二区三区| 亚洲一区二区三区欧美精品| 一级片'在线观看视频| 老司机深夜福利视频在线观看| 欧美一区二区精品小视频在线| 黄色视频不卡| 黄色片一级片一级黄色片| 亚洲欧美日韩高清在线视频| 亚洲国产精品合色在线| 88av欧美| 亚洲国产毛片av蜜桃av| 99国产精品一区二区蜜桃av| 国产人伦9x9x在线观看| 老司机靠b影院| 国产精品98久久久久久宅男小说| 精品乱码久久久久久99久播| 电影成人av| 正在播放国产对白刺激| 国产精品成人在线| 国产野战对白在线观看| 五月开心婷婷网| 中国美女看黄片| 国产精品久久电影中文字幕| 在线十欧美十亚洲十日本专区| 麻豆成人av在线观看| 精品国产一区二区三区四区第35| 在线天堂中文资源库| 久久精品91无色码中文字幕| 亚洲avbb在线观看| 欧美日韩瑟瑟在线播放| 一级a爱片免费观看的视频| 亚洲九九香蕉| 日本撒尿小便嘘嘘汇集6| 人人妻人人添人人爽欧美一区卜| 亚洲欧美精品综合一区二区三区| 在线观看一区二区三区激情| 啦啦啦 在线观看视频| 日韩欧美免费精品| 精品熟女少妇八av免费久了| 日日爽夜夜爽网站| 每晚都被弄得嗷嗷叫到高潮| a级片在线免费高清观看视频| av超薄肉色丝袜交足视频| 老司机福利观看| 91在线观看av| 日本五十路高清| av视频免费观看在线观看| 少妇裸体淫交视频免费看高清 | 少妇裸体淫交视频免费看高清 | 桃红色精品国产亚洲av| 国产精品久久电影中文字幕| www日本在线高清视频| 91麻豆精品激情在线观看国产 | 国产真人三级小视频在线观看| 午夜精品在线福利| 手机成人av网站| 亚洲精品美女久久久久99蜜臀| 国产单亲对白刺激| 18禁观看日本| 青草久久国产| 热99国产精品久久久久久7| 亚洲成av片中文字幕在线观看| 亚洲免费av在线视频| 在线观看www视频免费| 18禁美女被吸乳视频| av国产精品久久久久影院| 国产一区在线观看成人免费| 亚洲人成电影免费在线| 侵犯人妻中文字幕一二三四区| 国产欧美日韩精品亚洲av| 老司机福利观看| 三上悠亚av全集在线观看| 免费搜索国产男女视频| 99国产精品99久久久久| 1024香蕉在线观看| 亚洲男人的天堂狠狠| 自线自在国产av| 在线观看午夜福利视频| 99国产精品一区二区三区| 国产成年人精品一区二区 | 美女高潮喷水抽搐中文字幕| 91麻豆av在线| 99久久久亚洲精品蜜臀av| 亚洲欧美一区二区三区黑人| 国产激情欧美一区二区| 99热只有精品国产| 熟女少妇亚洲综合色aaa.| 好男人电影高清在线观看| 国产黄色免费在线视频| 国产精品国产高清国产av| 成人免费观看视频高清| 精品一区二区三区av网在线观看| 久久精品aⅴ一区二区三区四区| 国产麻豆69| 丰满饥渴人妻一区二区三| 亚洲国产精品一区二区三区在线| 免费一级毛片在线播放高清视频 | 99riav亚洲国产免费| 国产成+人综合+亚洲专区| 在线免费观看的www视频| 婷婷丁香在线五月| 亚洲少妇的诱惑av| 99在线视频只有这里精品首页| 国产乱人伦免费视频| 90打野战视频偷拍视频| 日日爽夜夜爽网站| 天堂动漫精品| 久久热在线av| 国产一区在线观看成人免费| 国产精品亚洲av一区麻豆| 女人被躁到高潮嗷嗷叫费观| 夜夜爽天天搞| 久久久久国产精品人妻aⅴ院| 国产成人精品在线电影| 真人做人爱边吃奶动态| 深夜精品福利| 精品久久久久久成人av| 国产真人三级小视频在线观看| 久久久精品国产亚洲av高清涩受| 纯流量卡能插随身wifi吗| www.www免费av| 欧美精品啪啪一区二区三区| 亚洲中文av在线| 久久影院123| 国产又色又爽无遮挡免费看| 久久草成人影院| 伊人久久大香线蕉亚洲五| 电影成人av| 精品国产乱子伦一区二区三区| 亚洲色图av天堂| 国产精品一区二区免费欧美| 一区二区三区激情视频| 两个人免费观看高清视频| 一边摸一边抽搐一进一小说| 天堂影院成人在线观看| 每晚都被弄得嗷嗷叫到高潮| 天堂动漫精品| 如日韩欧美国产精品一区二区三区| 免费在线观看日本一区| 激情视频va一区二区三区| 免费在线观看完整版高清| 欧美在线一区亚洲| 久9热在线精品视频| 高清av免费在线| 高清欧美精品videossex| 久久精品国产亚洲av香蕉五月| 久久性视频一级片| 真人做人爱边吃奶动态| 午夜免费成人在线视频| 久久久久久大精品| 两人在一起打扑克的视频| 午夜福利影视在线免费观看| 精品一品国产午夜福利视频| 午夜免费激情av| 亚洲av成人av| 91麻豆av在线| 成人三级黄色视频| 老司机午夜福利在线观看视频| 午夜免费成人在线视频| 亚洲欧美激情综合另类| 国产精品影院久久| 操美女的视频在线观看| 91成人精品电影| 性少妇av在线| www.自偷自拍.com| 99re在线观看精品视频| 国产成人欧美在线观看| 亚洲国产精品sss在线观看 | 亚洲欧美日韩高清在线视频| 搡老岳熟女国产| 国产区一区二久久| 好男人电影高清在线观看| 别揉我奶头~嗯~啊~动态视频| 久久久久久免费高清国产稀缺| 中国美女看黄片| 97人妻天天添夜夜摸| 满18在线观看网站| 又黄又爽又免费观看的视频| 久久亚洲真实| 国产精品免费视频内射| 一边摸一边抽搐一进一小说| 亚洲五月色婷婷综合| 亚洲国产精品一区二区三区在线| 久久久久国产精品人妻aⅴ院| 久久人妻熟女aⅴ| 国产伦一二天堂av在线观看| 久久人妻熟女aⅴ| 国产成+人综合+亚洲专区| 黄色视频,在线免费观看| 亚洲av成人一区二区三| 免费观看人在逋| 成人永久免费在线观看视频| 一夜夜www| 国产免费男女视频| 亚洲一码二码三码区别大吗| 色老头精品视频在线观看| 一夜夜www| 两性午夜刺激爽爽歪歪视频在线观看 | 神马国产精品三级电影在线观看 | av免费在线观看网站| 亚洲精品美女久久久久99蜜臀| 成年女人毛片免费观看观看9| 国产一区二区三区视频了| 亚洲黑人精品在线| 亚洲一区二区三区不卡视频| 99国产精品免费福利视频| av网站在线播放免费| 日韩免费av在线播放| 高清黄色对白视频在线免费看| 欧美久久黑人一区二区| 巨乳人妻的诱惑在线观看| 亚洲黑人精品在线| 天堂动漫精品| 热99re8久久精品国产| 亚洲中文字幕日韩| 久久精品人人爽人人爽视色| 99久久精品国产亚洲精品| 黑人巨大精品欧美一区二区mp4| av在线天堂中文字幕 | 国产精华一区二区三区| 久热爱精品视频在线9| 免费搜索国产男女视频| 亚洲精品成人av观看孕妇| 又大又爽又粗| 午夜两性在线视频| 欧美精品一区二区免费开放| 国产精品综合久久久久久久免费 | 色婷婷av一区二区三区视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品美女特级片免费视频播放器 | 少妇粗大呻吟视频| 国产成人一区二区三区免费视频网站| 村上凉子中文字幕在线| 成熟少妇高潮喷水视频| av视频免费观看在线观看| 国产高清国产精品国产三级| 国产欧美日韩精品亚洲av| 国产1区2区3区精品| 国产精品爽爽va在线观看网站 | 99热国产这里只有精品6| 波多野结衣高清无吗| 国产亚洲欧美98| 一级a爱片免费观看的视频| 高清av免费在线| 久热这里只有精品99| 国产精品一区二区在线不卡| 一二三四社区在线视频社区8| 人妻丰满熟妇av一区二区三区| 亚洲精品国产精品久久久不卡| 中文字幕av电影在线播放| 成人亚洲精品一区在线观看| 亚洲av成人一区二区三| 最近最新中文字幕大全免费视频| 亚洲欧美日韩另类电影网站| 99在线视频只有这里精品首页| 大型黄色视频在线免费观看| 桃色一区二区三区在线观看| 视频区图区小说| 嫩草影视91久久| 国产区一区二久久| 久久国产精品男人的天堂亚洲| 校园春色视频在线观看| 国产精品免费一区二区三区在线| 丁香六月欧美| 欧美激情高清一区二区三区| 九色亚洲精品在线播放| 男人操女人黄网站| 女人高潮潮喷娇喘18禁视频| 久久久国产一区二区| 久久天堂一区二区三区四区| 可以在线观看毛片的网站| 一个人免费在线观看的高清视频| 欧美激情久久久久久爽电影 | 99热国产这里只有精品6| 日本黄色视频三级网站网址| 在线观看免费视频日本深夜| 啦啦啦 在线观看视频| 亚洲精品国产色婷婷电影| 久久久国产欧美日韩av| 中文字幕人妻熟女乱码| 国产精品1区2区在线观看.| 天天添夜夜摸| 国产成人av激情在线播放| 亚洲美女黄片视频| 18禁美女被吸乳视频| 久久久国产一区二区| 午夜精品国产一区二区电影| 日日干狠狠操夜夜爽| 亚洲 欧美一区二区三区| 精品久久久久久,| 亚洲aⅴ乱码一区二区在线播放 | 精品熟女少妇八av免费久了| 女同久久另类99精品国产91| 精品免费久久久久久久清纯| 日日摸夜夜添夜夜添小说| 亚洲成人国产一区在线观看| 亚洲成国产人片在线观看| 国产精品综合久久久久久久免费 | 亚洲男人的天堂狠狠| 久久人妻av系列| 19禁男女啪啪无遮挡网站| 搡老熟女国产l中国老女人| 人成视频在线观看免费观看| 国产一区二区三区视频了| 一进一出抽搐动态| 久久国产精品男人的天堂亚洲| 正在播放国产对白刺激| 色综合站精品国产| 欧美精品啪啪一区二区三区| 亚洲人成电影观看| 亚洲 国产 在线| 久久久久国产一级毛片高清牌| 80岁老熟妇乱子伦牲交| 黄色怎么调成土黄色| 日韩视频一区二区在线观看| 超碰成人久久| 这个男人来自地球电影免费观看| 亚洲 国产 在线| 色综合欧美亚洲国产小说| 久久狼人影院| www.精华液| 日本五十路高清| 久热这里只有精品99| 热99re8久久精品国产| 日韩中文字幕欧美一区二区| 亚洲五月天丁香| 久久精品91蜜桃| 国产色视频综合| 亚洲国产欧美日韩在线播放| 欧美日韩视频精品一区| 免费在线观看影片大全网站| 无人区码免费观看不卡| 亚洲欧美日韩无卡精品| 欧美日韩亚洲国产一区二区在线观看| 国产av一区二区精品久久| 久久久久国产一级毛片高清牌| 啦啦啦免费观看视频1| 人人妻人人澡人人看| 嫩草影视91久久| 久久久久久亚洲精品国产蜜桃av| 搡老岳熟女国产| 欧美另类亚洲清纯唯美| 久久天堂一区二区三区四区| 国产高清国产精品国产三级| 欧美午夜高清在线| 精品久久久久久久久久免费视频 | 中文字幕人妻丝袜制服| av国产精品久久久久影院| 88av欧美| 人人澡人人妻人| 深夜精品福利| 免费一级毛片在线播放高清视频 | videosex国产| 精品日产1卡2卡| 黑人巨大精品欧美一区二区蜜桃| 日韩大码丰满熟妇| 侵犯人妻中文字幕一二三四区| 亚洲av第一区精品v没综合| 国产伦一二天堂av在线观看| 天堂√8在线中文| 怎么达到女性高潮| 日本黄色日本黄色录像| 免费在线观看日本一区| 十八禁网站免费在线| 国产精品电影一区二区三区| 国产精品日韩av在线免费观看 | 淫妇啪啪啪对白视频| 老汉色av国产亚洲站长工具| 国产99久久九九免费精品| 亚洲精品国产精品久久久不卡| 日日摸夜夜添夜夜添小说| 99国产精品免费福利视频| 伦理电影免费视频| 欧美老熟妇乱子伦牲交| 久久午夜综合久久蜜桃| 欧美老熟妇乱子伦牲交| 一级毛片高清免费大全| 久久精品aⅴ一区二区三区四区| 丁香六月欧美| 久久亚洲真实| 乱人伦中国视频| 欧美大码av| 精品国产乱码久久久久久男人| 精品一区二区三区视频在线观看免费 | 国产1区2区3区精品| 一级片免费观看大全| 乱人伦中国视频| 国产高清国产精品国产三级| 少妇裸体淫交视频免费看高清 | 这个男人来自地球电影免费观看| 国产男靠女视频免费网站| 久久人妻av系列| 性色av乱码一区二区三区2| 免费在线观看影片大全网站| 欧美中文日本在线观看视频| 婷婷六月久久综合丁香| 亚洲在线自拍视频| 黄色视频不卡| 99精品在免费线老司机午夜| 嫩草影院精品99| 亚洲一区中文字幕在线| 国产成人影院久久av| 丝袜在线中文字幕| 99热国产这里只有精品6| 最好的美女福利视频网| 国产成人啪精品午夜网站| 十八禁网站免费在线| 日日夜夜操网爽| 国产一区二区激情短视频| 亚洲国产欧美网| 亚洲av熟女| 动漫黄色视频在线观看| 久久久国产成人免费| 精品福利观看| 久久精品亚洲av国产电影网| 乱人伦中国视频| www.熟女人妻精品国产| 久久久国产精品麻豆| 欧美精品亚洲一区二区| 久9热在线精品视频| 欧美av亚洲av综合av国产av| 黑人猛操日本美女一级片| 人人妻人人爽人人添夜夜欢视频| 在线观看舔阴道视频| 免费看十八禁软件| 日本免费一区二区三区高清不卡 | 伦理电影免费视频| 淫秽高清视频在线观看| 国产精品1区2区在线观看.| 在线观看日韩欧美| 亚洲成人精品中文字幕电影 | 女性被躁到高潮视频| 在线观看午夜福利视频| 欧美日韩亚洲高清精品| 国产区一区二久久| 80岁老熟妇乱子伦牲交| xxxhd国产人妻xxx| 久久人妻福利社区极品人妻图片| 国产精品1区2区在线观看.| tocl精华| 国产野战对白在线观看| 国内久久婷婷六月综合欲色啪|