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

    華西牛胴體性狀GWAS 分析

    2022-11-17 03:52:46王澤昭李宏偉安炳星李海鵬徐凌洋蔡文濤張路培高會江李俊雅
    中國畜牧雜志 2022年11期
    關鍵詞:華西胴體表型

    高 翰,葛 菲,王澤昭,李宏偉,安炳星,李海鵬,徐凌洋,朱 波,蔡文濤,張路培,高 雪,高會江,陳 燕,李俊雅

    (中國農(nóng)業(yè)科學院北京畜牧獸醫(yī)研究所,北京 100193)

    我國是肉牛生產(chǎn)大國,也是牛肉消費大國[1]。隨著居民生活水平的提高,消費者對牛肉的需求持續(xù)增長,供需缺口呈擴大趨勢[2]。胴體性狀作為肉牛生產(chǎn)實踐中最重要的經(jīng)濟性狀之一,是典型的由多基因調(diào)控的數(shù)量性狀,遺傳調(diào)控機制復雜[3-4]。Risch 于1996 年首次提出全基因組關聯(lián)分析(Genome-Wide Association Study,GWAS)的應用[5],隨著測序技術的不斷發(fā)展和測序成本的降低,GWAS 被廣泛應用于研究家畜數(shù)量性狀的研究中[6-10]。利用GWAS,在肉牛中鑒定出了PLAG1[11]、IGF2[12]、MTPN[13]等影響胴體性狀的重要基因。大部分GWAS 都是基于單一性狀與SNPs 的關聯(lián)分析,即使同時分析多個相關表型,也只是逐個性狀分別進行單一性狀GWAS[14-15]。Bolormaa 等[16]研究指出,對于具有高估計育種值(Estimated Breeding Value,EBV)準確性的性狀,多性狀GWAS 的效力至少會與單性狀GWAS 的效力一樣高。如果性狀間高度相關,則多性狀分析相較于單性狀分析更具有優(yōu)勢。因為多性狀GWAS 只進行一次統(tǒng)計檢驗,同時考慮多個性狀的性狀內(nèi)和性狀間方差組分[17],降低了多重檢驗造成的誤差[18],提高了檢驗功效[19-20]和參數(shù)估計的精度[21]。因此,與逐個性狀與遺傳位點間的關聯(lián)分析相比,多個性狀聯(lián)合信息的多性狀GWAS 通常具有更好的關聯(lián)分析效果[22-24]。

    “華西?!笔侵袊r(nóng)業(yè)科學院北京畜牧獸醫(yī)研究所歷經(jīng)40 余年潛心培育的肉牛新品種,2021 年通過國家畜禽遺傳資源委員會審定,并獲得國家畜禽新品種證書,其具有生長速度快、飼料轉化率高、凈肉率高、產(chǎn)肉和繁殖性能好、抗逆性強等特性。與國際同類型肉牛品種相比,“華西?!钡娜赵鲋?、屠宰率、凈肉率處于國際先進水平。本研究以內(nèi)蒙古烏拉蓋管理區(qū)建立的華西牛資源群體為研究對象,基于Illumina BovineHD 芯片基因型數(shù)據(jù),對胴體重(Carcass Weight,CW)、胴體長(Carcass Length,CL)、胴體胸深(Chest Depth of Carcass,CD)3 個胴體性狀分別進行單性狀與多性狀GWAS 分析,旨在挖掘影響華西牛胴體性狀的顯著SNPs 位點及候選基因,為全基因組選擇育種提供有效的分子標記,為進一步研究華西牛胴體相關性狀的潛在遺傳機制提供有價值的參考。

    1 材料與方法

    1.1 實驗數(shù)據(jù)收集 實驗動物群體選自中國農(nóng)業(yè)科學院北京畜牧獸醫(yī)研究所牛遺傳育種創(chuàng)新團隊在內(nèi)蒙古錫林郭勒盟烏拉蓋管理區(qū)組建的華西牛資源群體,該群體2008 年開始建立,每年進行擴群,并且使用Illumina BovineHD 芯片對每個個體進行基因分型。胴體重、胴體長和胴體胸深3 個胴體性狀表型數(shù)據(jù)嚴格按照GB/T 27643-2011《牛胴體及鮮肉分割》[25]的要求進行測定。

    1.2 表型處理和基因型質(zhì)控 對表型數(shù)據(jù)進行處理,包括剔除表型缺失值、三倍標準差以外的異常值。采用一般線性模型評估育肥場、屠宰場、屠宰年份、屠宰季節(jié)以及年齡對性狀影響的顯著性,將顯著的效應作為協(xié)變量加入到模型當中。利用PLINK 軟件[26]對群體進行主成分分析(Principal Components Analysis,PCA),使用R 包ggplot2 繪制群體PCA 圖。利用EIGENSTRAT計算各主成分是否具有顯著統(tǒng)計學意義,將達到顯著水平P<0.05 的主成分作為協(xié)變量加入到模型當中。使用PLINK 軟件對相應個體的基因型數(shù)據(jù)進行質(zhì)量控制,標準如下:剔除SNP 缺失率大于10%的個體以及檢出率低于90%、最小哈德溫伯格平衡小于1.0×10-6、最小等位基因頻率小于0.05 的SNPs 位點。最后,使用R包ASReml[27]中的多性狀動物模型計算各性狀之間的表型相關和遺傳相關。

    1.3 關聯(lián)統(tǒng)計分析

    1.3.1 單性狀GWAS 使用GEMMA 軟件[21]進行全基因組關聯(lián)分析,采用考慮親緣關系的混合線性模型(Linear Mixed Model,LMM),計算公式如下:

    式(1)中,y為所有個體表型性狀的n×1 向量;W為協(xié)變量矩陣,α為包括截距在內(nèi)對應系數(shù)的一列向量;x為標記基因型向量,β表示標記位點效應的大??;u為個體的隨機效應向量,u~N(0,KVg),K代表由SNP 標記計算出的已知的n×n遺傳關系矩陣,Vg為加性遺傳方差;ε代表n×1 隨機誤差向量,符合ε~N(0,IVe)分布,I代表n×n身份矩陣,Ve代表殘差組分。本次研究中主要利用Wald 測驗作為GWAS 分析顯著統(tǒng)計指標。

    1.3.2 多性狀GWAS 使用GEMMA 軟件[21]中的多性狀混合線性模型(Multivariate Linear Mixed Model,mvLMM)進行多個性狀表型值和SNP 標記的GWAS,計算公式如下:

    式(2)中,Y為n×t維表型矩陣,n是樣本數(shù)目,t是分析性狀個數(shù);W為n×c維協(xié)變量矩陣,A為c維相應系數(shù)行向量,c是包含截距項在內(nèi)的行變量個數(shù);x為n維當前檢驗SNP 的基因型指示變量列向量,β為n×t維不包括當前檢驗SNP 的剩余多基因效應矩陣;G~MN(0,Vg?K),K是由全基因組SNP 標記構建的基因組親緣關系矩陣,Vg是t×t維剩余多基因方差-協(xié)方差矩陣;E為n×t維誤差項矩陣,E~MN(0,Ve?I),Ve是t×t維誤差方差-協(xié)方差矩陣,I是單位矩陣。最終,全基因組上每個SNP 都能得到1 個一因多效Wald 統(tǒng)計量。應用Bonferroni 方法進行多重校正,確定染色體水平顯著性閾值為1/N,N 為用于分析的SNPs 數(shù)目。最后,用R 包ggplot2 繪制曼哈頓圖和QQ 圖。

    1.4 位點版本轉化及基因注釋 Illumina BovineHD 芯片中的SNPs 位置信息為Bos taurus UMD 3.1 版本,使用在線網(wǎng)站UCSC 中的liftover 工具(http://genome.ucsc.edu/cgi-bin/hgLiftOver/)將其轉化成最新的?;蚪M版本ARS-UCD 1.2。根據(jù)得到的顯著SNPs 位點的物理位置,利用Ensemble 在線數(shù)據(jù)庫的biomart 模塊(http://www.ensembl.org/biomart/)搜索每個顯著SNP位點上下游50 kb 區(qū)域尋找距離最近的候選基因,作為與目標性狀關聯(lián)的重要功能基因。

    2 結果

    2.1 主成分分析 由圖1 可知,實驗群體存在分層現(xiàn)象。利用EIGENSTRAT 計算發(fā)現(xiàn),有2 個達到P<0.05 顯著水平的主成分,因此在后續(xù)分析中將前2 個主成分作為協(xié)變量加入到模型中。

    圖1 群體主成分分析圖

    2.2 表型及基因型描述性統(tǒng)計

    2.2.1 表型值的描述性統(tǒng)計 分別對840 頭華西牛公牛的CW、CL 和CD 3 個性狀開展單性狀GWAS 與多性狀GWAS 分析。表1 為各性狀的原始表型描述性統(tǒng)計分析結果。3 個性狀的表型頻率分布圖見圖2,可以看出,各性狀的表型值基本都符合正態(tài)分布。此外,對各表型值進行了表型相關分析與遺傳相關分析(表2)。結果顯示,胴體重、胴體長、胴體胸深3 個性狀間均具有較高的表型相關(r>0.6);胴體長與胴體胸深具有中等遺傳相關性(r=0.339 2),胴體重與胴體長、胴體胸深具有強遺傳相關性(r>0.4)。

    表1 胴體重、胴體長、胴體胸深的描述性統(tǒng)計

    表2 性狀間表型及遺傳相關程度

    圖2 CW、CL 和CD 的頻率分布

    2.2.2 基因型的描述性統(tǒng)計 按篩選原則對基因型數(shù)據(jù)進行質(zhì)控后,共剩余607 198 個SNPs 標記。由圖3 可知,各染色體上的SNP 位點分布較均勻,可以用于華西牛群體的GWAS 分析。

    圖3 SNPs 標記在29 條染色體的分布情況

    2.3 GWAS 分析 使用GEMMA 軟件中考慮親緣關系的混合線性模型(LMM)對3 個性狀分別做單性狀GWAS 分析。曼哈頓和QQ 圖如圖4 所示,在CW 和CL 性狀中都沒有檢測到顯著的遺傳變異位點,僅在CD 性狀中發(fā)現(xiàn)了6 個顯著的SNPs 位點,分別位于4、7、10、23、24、25 號染色體上,共注釋到了5 個基因(表3)。

    表3 單性狀GWAS 和多性狀GWAS 基因注釋結果

    圖4 單性狀GWAS 曼哈頓和QQ 圖

    使用mvLMM 分別進行了CW-CL、CW-CD、CLCD、CW-CL-CD 4 個組合的多性狀GWAS 分析(圖5)。在CW-CL 組合的多性狀GWAS 中,發(fā)現(xiàn)1 個超過染色體顯著水平的SNP 位點,位于1 號染色體上,注釋到了1 個基因;在CW-CD 組合的多性狀GWAS 中,發(fā)現(xiàn)5 個超過染色體顯著水平的SNPs 位點,分別位于4、7、23、25 號染色體上,共注釋到3 個基因;在CL-CD組合的多性狀GWAS 中,發(fā)現(xiàn)2 個超過染色體顯著水平的SNPs 位點,分別位于4、25 號染色體上,注釋到了1 個基因;在CW-CL-CD 組合的多性狀GWAS 中,發(fā)現(xiàn)3 個超過染色體顯著水平的SNPs 位點,分別位于4、7、23 號染色體上,共注釋到3 個基因。合并多性狀GWAS 的4 個組合的結果發(fā)現(xiàn),有4 個SNPs 位點在多性狀GWAS 中被重復檢測到,其中Hapmap36353-SCAFFOLD29708_3468 在CW-CD、CL-CD、CW-CLCD 3 個組合中被重復檢測到,BovineHD0700018227和BovineHD2300004986均在CW-CD、CW-CL-CD 2個組合中被重復檢測到,BovineHD2500006960 在CLCD、CW-CD 2 個組合中被重復檢測到。去除重復位點,多性狀GWAS 共找到6 個一因多效的SNPs 位點。

    圖5 多性狀GWAS 曼哈頓和QQ 圖

    2.4 重要候選基因 對CW、CL、CD 3 個性狀分別進行單性狀GWAS 和不同性狀之間組合的多性狀GWAS結果進行匯總統(tǒng)計。結果得到,單性狀與多性狀GWAS共檢測到9 個與華西牛胴體性狀顯著相關的SNPs 位點,根據(jù)顯著SNPs 位點的物理位置,在這9 個顯著SNPs 位點上下游50 kb 的區(qū)間內(nèi)尋找距離最近的候選基因,共在8 個位點附近找到6 個候選基因,分別是突觸結合蛋白5L(Syntaxin Binding Protein 5 Like,STXBP5L)、磷酸二酯酶1C(Phosphodiesterase 1C,PDE1C)、過氧化物酶體增殖物激活受體γ輔激活因子-1β(Peroxisome Proliferator-Activated Receptor Gamma Coactivator 1-Beta,PPARGC1B)、(Metaxin 3,MTX3)、鈣調(diào)磷酸酶調(diào)節(jié)因子2(Regulator Of Calcineurin 2,RCAN2)、??康鞍?(Docking Protein 6,DOK6)。其中,PDE1C、PPARGC1B、RCAN23 個基因在單性狀GWAS 與多性狀GWAS 中被重復檢測到。

    由圖6 可知,單性狀與多性狀GWAS 重復檢測到的3 個顯著SNPs 位點分別是Hapmap36353-SCAFFOLD29708_3468、BovineHD2300004986 和Bovine HD2500006960。其中,單性狀GWAS 與多性狀GWAS中CW-CD、CL-CD 組合共同檢測到的BovineHD25000 06960 未注釋到相關基因。此外,存在不同策略檢測到的多個SNPs 位點注釋到同一個基因的情況,具體為:在CD 單性狀GWAS 檢測到的BovineHD0700018210 與CW-CD、CW-CL-CD 2 個多性狀GWAS 共同檢測到的BovineHD0700018227,同時注釋到基因PPARGC1B;CW-CD 多性狀GWAS 檢測到的BovineHD2300004985與單性狀GWAS 和CW-CD、CW-CL-CD 2 個多性狀GWAS 共同檢測到的BovineHD2300004986,同時注釋到基因RCAN2。PPARGC1B和RCAN2將作為后續(xù)重點關注的與胴體性狀關聯(lián)的重要候選基因。

    圖6 單性狀GWAS 與多性狀GWAS 結果韋恩圖

    3 討 論

    在肉牛產(chǎn)業(yè)中,胴體性狀被認為是影響肉牛生產(chǎn)的重要經(jīng)濟性狀,對于提高生產(chǎn)效率和效益起著至關重要的作用[28-29]。然而,肉牛胴體性狀的測定只能在個體屠宰后進行,這就導致胴體性狀的相關表型數(shù)據(jù)收集困難且成本高昂[30]。因此,闡明胴體性狀的遺傳調(diào)控機制,提供可靠的分子遺傳標記,可以對肉牛胴體相關性狀進行早期遺傳評估,從而降低育種成本,提高生產(chǎn)效益。

    GWAS 方法已經(jīng)廣泛應用于動植物復雜性狀的遺傳解析[31-33]。本研究基于Illumina BovineHD 芯片基因分型,對華西牛CW、CL、CD 3 個胴體性狀進行了關聯(lián)分析,綜合考慮單性狀和多性狀GWAS 分析得到的結果,共得到9 個顯著的SNPs 位點和6 個相關基因,分別是STXBP5L、PDE1C、PPARGC1B、MTX3、RCAN2和DOK6。其中,多性狀與單性狀策略共同鑒定到3 個基因,分別為PDE1C、PPARGC1B和RCAN2。PDE1C基因在小鼠胚胎早期發(fā)育過程中中樞神經(jīng)系統(tǒng)內(nèi)遷移的神經(jīng)元細胞中表達[34]。在人類疾病的研究中,PDE1C基因的異常會導致兒童發(fā)育遲緩[35]。PPARGC1B基因是一種調(diào)控基因轉錄的共激活劑,對轉錄因子和核受體活動、脂肪氧化、糖酵解以及能量代謝等多種生物學過程具有調(diào)控作用[36-39]。Jiang 等[40]研究了秦川牛犢牛和成年發(fā)育階段牛脂肪組織中circRNA 的表達,發(fā)現(xiàn)circFUT10 在脂肪細胞增殖和分化中有重要作用,circFUT10 間接促進PPARGC1B的轉錄,然后調(diào)節(jié)細胞的增殖和分化,且PPARGC1B在成年牛脂肪組織中表達水平較高,說明該基因對胴體性狀的表型差異可能存在重要影響。在小鼠模型中的相關研究表明,RCAN2基因對體重調(diào)節(jié)有重要作用[41-43]。RCAN2在成纖維細胞、心臟、腦、肝臟和骨骼肌中表達,能夠抑制鈣調(diào)磷酸酶活性,在骨骼發(fā)育和維護中有重要作用[44]。由此可見,該基因可能是影響華西牛胴體性狀的重要候選基因。MTX3與DOK6是單性狀GWAS 分析中鑒定到的與CD 性狀相關的2 個特有基因。MTX3在蛋白質(zhì)轉運到線粒體這一生物過程中發(fā)揮重要作用[45-46]。DOK6能促進RET 介導的神經(jīng)突生長,可能在大腦發(fā)育或維護中發(fā)揮作用[47]。Jiao 等[48]對杜洛克豬群體進行GWAS 分析,提出DOK6基因為平均日采食量、育肥期日增重、活體背膘厚3 個重要經(jīng)濟性狀最有可能的候選基因之一。STXBP5L基因為CW-CL 組合的多性狀GWAS 中單獨鑒定到的基因,它在維持葡萄糖穩(wěn)態(tài)[49]、囊泡運輸和胞吐作用抑制[50]等過程中發(fā)揮著重要作用。

    對單性狀GWAS 與多性狀GWAS 結果進行比較,結果發(fā)現(xiàn),在相同的顯著性閾值下,兩者檢測到的基因數(shù)目無明顯差別。但是單一性狀GWAS 分析中,CW、CL 性狀都沒有檢測到顯著SNPs 位點及相關基因。通過將CW 與CL、CD 兩兩組合進行多性狀GWAS,在CW-CL、CW-CD、CL-CD、CW-CL-CD 4 個不同組合中均檢測到了不同數(shù)目的顯著SNPs 位點和相關基因。同時,在單性狀GWAS 以及4 組多性狀GWAS 中,Hapmap36353-SCAFFOLD29708_3468、BovineHD070 0018227、BovineHD2300004986、BovineHD2500006960共4 個SNPs 位點,以及PDE1C、PPARGC1B、RCAN2共3 個基因被重復檢測到,推測認為這些位點及基因極大可能是潛在的一因多效位點和一因多效基因??傮w而言,多性狀GWAS 可以提高對位點的檢測效率和準確性,這與高進等[51]、Bolormaa 等[16]的研究結果相一致。在實際育種當中,CW、CL 和CD 的表型通常是同時被選擇的,由于多性狀GWAS 分析檢測到的SNPs 位點具有一因多效的優(yōu)越性,因此綜合考慮多性狀GWAS 分析在生物學和實際應用層面更有意義。

    4 結 論

    本研究基于Illumina BovineHD 高密度基因分型結果,利用840 頭華西牛胴體重、胴體長和胴體胸深表型數(shù)據(jù),采用單變量和多變量線性混合模型方法進行單性狀和多性狀的GWAS 分析,共檢測到9 個顯著關聯(lián)的SNPs 位點,注釋到STXBP5L、PDE1C、PPARGC1B、MTX3、RCAN2、DOK6共6個候選基因,其中PPARGC1B、RCAN2被不同方法檢測到的多個SNPs 位點同時注釋到,且參與體重調(diào)節(jié)、脂肪細胞增殖和分化等過程,推測其可能是影響胴體性狀的重要候選基因。

    猜你喜歡
    華西胴體表型
    影響豬胴體瘦肉率的因素及提高措施
    建蘭、寒蘭花表型分析
    百年精誠 譽從信來——走進四川大學華西眼視光之一
    在華西人與晚清軍事技術近代化
    華西追蹤:“戰(zhàn)時狀態(tài)”不見了
    內(nèi)洛爾公牛的飼料轉化率及其與胴體品質(zhì)、非胴體品質(zhì)和肉品質(zhì)的關系
    飼料博覽(2016年5期)2016-07-12 11:58:39
    GABABR2基因遺傳變異與肥胖及代謝相關表型的關系
    藍塘豬與長白豬正反交F1代胴體性狀和肉品質(zhì)的比較
    慢性乙型肝炎患者HBV基因表型與血清學測定的臨床意義
    72例老年急性白血病免疫表型分析
    日本三级黄在线观看| 国产一区二区三区在线臀色熟女| 亚洲在线自拍视频| 级片在线观看| 午夜福利在线观看吧| 女生性感内裤真人,穿戴方法视频| 国产精品久久久久久精品电影 | 丝袜美足系列| 天堂影院成人在线观看| 性欧美人与动物交配| 久久性视频一级片| 18禁美女被吸乳视频| 女性生殖器流出的白浆| 男女床上黄色一级片免费看| 黄色女人牲交| 好男人电影高清在线观看| 黄网站色视频无遮挡免费观看| 国产亚洲精品av在线| 久久午夜综合久久蜜桃| 法律面前人人平等表现在哪些方面| 夜夜夜夜夜久久久久| 欧美色视频一区免费| av视频在线观看入口| 国产亚洲欧美在线一区二区| 88av欧美| 午夜免费观看网址| 亚洲精品久久成人aⅴ小说| 久久久国产欧美日韩av| 91字幕亚洲| 黑人巨大精品欧美一区二区mp4| 99国产精品一区二区三区| 一区二区三区国产精品乱码| xxx96com| 亚洲国产中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 久久午夜亚洲精品久久| 亚洲精华国产精华精| 又大又爽又粗| 国内毛片毛片毛片毛片毛片| av在线天堂中文字幕| 91成年电影在线观看| 亚洲精品美女久久久久99蜜臀| 国产91精品成人一区二区三区| 久久狼人影院| 老熟妇仑乱视频hdxx| 欧美日韩精品网址| 欧美成人一区二区免费高清观看 | 在线天堂中文资源库| 亚洲成人国产一区在线观看| 亚洲激情在线av| 一级a爱片免费观看的视频| 最好的美女福利视频网| 乱人伦中国视频| 每晚都被弄得嗷嗷叫到高潮| 成人av一区二区三区在线看| 国产99白浆流出| 国产精品免费视频内射| av免费在线观看网站| 一边摸一边抽搐一进一出视频| 欧美乱码精品一区二区三区| 99久久99久久久精品蜜桃| 精品一区二区三区av网在线观看| 国产人伦9x9x在线观看| 色尼玛亚洲综合影院| 校园春色视频在线观看| 欧美另类亚洲清纯唯美| 男女下面插进去视频免费观看| 亚洲精品久久成人aⅴ小说| 一级a爱视频在线免费观看| 亚洲欧洲精品一区二区精品久久久| 麻豆av在线久日| av视频免费观看在线观看| 9热在线视频观看99| 成人18禁高潮啪啪吃奶动态图| 午夜福利视频1000在线观看 | 欧美日韩一级在线毛片| 美女高潮到喷水免费观看| 亚洲自拍偷在线| 手机成人av网站| 多毛熟女@视频| 国产精品亚洲av一区麻豆| 亚洲av电影不卡..在线观看| 中文字幕最新亚洲高清| 50天的宝宝边吃奶边哭怎么回事| 精品久久久久久久人妻蜜臀av | 麻豆成人av在线观看| av免费在线观看网站| 天天躁夜夜躁狠狠躁躁| 一a级毛片在线观看| 最好的美女福利视频网| 国产三级在线视频| 亚洲少妇的诱惑av| 国产野战对白在线观看| 成人国产一区最新在线观看| 国产精品影院久久| 亚洲自拍偷在线| 男女下面进入的视频免费午夜 | 中文亚洲av片在线观看爽| 欧美日韩黄片免| 久久欧美精品欧美久久欧美| 波多野结衣av一区二区av| 男女做爰动态图高潮gif福利片 | 免费搜索国产男女视频| 男人操女人黄网站| 免费在线观看黄色视频的| e午夜精品久久久久久久| 国产成人av教育| www国产在线视频色| 正在播放国产对白刺激| 亚洲精品av麻豆狂野| 欧美日韩瑟瑟在线播放| 黑人巨大精品欧美一区二区mp4| 精品不卡国产一区二区三区| 午夜视频精品福利| 国产三级在线视频| 少妇被粗大的猛进出69影院| 国产黄a三级三级三级人| 一二三四社区在线视频社区8| 亚洲伊人色综图| 久热爱精品视频在线9| 国产亚洲精品一区二区www| 国产高清videossex| 操美女的视频在线观看| 日本 欧美在线| 精品久久久久久久毛片微露脸| 国产主播在线观看一区二区| ponron亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲片人在线观看| 妹子高潮喷水视频| 可以在线观看毛片的网站| 狂野欧美激情性xxxx| 午夜福利视频1000在线观看 | 自拍欧美九色日韩亚洲蝌蚪91| 自拍欧美九色日韩亚洲蝌蚪91| 久久人妻熟女aⅴ| 桃红色精品国产亚洲av| 国产麻豆成人av免费视频| 久久草成人影院| 精品第一国产精品| 色综合婷婷激情| 亚洲欧美精品综合一区二区三区| 午夜两性在线视频| 波多野结衣av一区二区av| 国产一卡二卡三卡精品| 啦啦啦 在线观看视频| 两个人视频免费观看高清| 国产成人影院久久av| 欧美乱妇无乱码| 99久久精品国产亚洲精品| 国产精品 欧美亚洲| netflix在线观看网站| 国产精品一区二区精品视频观看| 黄色 视频免费看| av在线播放免费不卡| 十八禁网站免费在线| 久久精品国产亚洲av高清一级| 午夜免费激情av| 欧美黄色淫秽网站| 最新美女视频免费是黄的| 精品卡一卡二卡四卡免费| 久久国产精品影院| 精品久久久久久,| 啦啦啦 在线观看视频| 国产成人精品无人区| 麻豆国产av国片精品| 精品久久蜜臀av无| 啦啦啦观看免费观看视频高清 | 亚洲无线在线观看| 欧美成人一区二区免费高清观看 | 日韩欧美国产在线观看| 无限看片的www在线观看| 午夜两性在线视频| 美女免费视频网站| 99精品欧美一区二区三区四区| 欧美激情高清一区二区三区| 国产精华一区二区三区| 欧美中文日本在线观看视频| 怎么达到女性高潮| 色婷婷久久久亚洲欧美| 天天躁夜夜躁狠狠躁躁| 他把我摸到了高潮在线观看| 日韩欧美免费精品| 欧美成狂野欧美在线观看| 亚洲欧美激情综合另类| 欧美国产日韩亚洲一区| 香蕉久久夜色| 97人妻天天添夜夜摸| 国产欧美日韩一区二区三区在线| 欧美黑人精品巨大| 成人免费观看视频高清| 十八禁人妻一区二区| 亚洲欧美激情综合另类| 9色porny在线观看| 99国产精品一区二区三区| 老熟妇仑乱视频hdxx| 国产私拍福利视频在线观看| 黄频高清免费视频| 99精品久久久久人妻精品| 国产av一区在线观看免费| 少妇的丰满在线观看| 久久久精品国产亚洲av高清涩受| 视频在线观看一区二区三区| 性色av乱码一区二区三区2| av有码第一页| 久久久久久国产a免费观看| 在线永久观看黄色视频| 成人三级黄色视频| 乱人伦中国视频| ponron亚洲| 一区二区三区高清视频在线| 国产午夜福利久久久久久| 国产亚洲精品久久久久久毛片| 757午夜福利合集在线观看| 不卡av一区二区三区| 在线观看舔阴道视频| 欧美日本视频| 亚洲成人免费电影在线观看| 午夜福利一区二区在线看| 91av网站免费观看| av福利片在线| 国产aⅴ精品一区二区三区波| 一本大道久久a久久精品| 黑人巨大精品欧美一区二区蜜桃| 岛国视频午夜一区免费看| 亚洲aⅴ乱码一区二区在线播放 | 人人妻,人人澡人人爽秒播| www.999成人在线观看| 丝袜美腿诱惑在线| 非洲黑人性xxxx精品又粗又长| 我的亚洲天堂| 在线观看66精品国产| 日韩成人在线观看一区二区三区| 欧美日本中文国产一区发布| 麻豆av在线久日| 亚洲国产精品久久男人天堂| 亚洲午夜精品一区,二区,三区| 久久欧美精品欧美久久欧美| 久久香蕉国产精品| 久9热在线精品视频| 欧洲精品卡2卡3卡4卡5卡区| 美国免费a级毛片| 亚洲一区二区三区色噜噜| 一本大道久久a久久精品| 婷婷精品国产亚洲av在线| 国产成人精品在线电影| 国产欧美日韩一区二区精品| 午夜久久久在线观看| 精品久久久精品久久久| 身体一侧抽搐| 男女下面进入的视频免费午夜 | 手机成人av网站| 日本五十路高清| 91精品三级在线观看| av超薄肉色丝袜交足视频| 99久久综合精品五月天人人| 久久精品国产综合久久久| 亚洲国产欧美日韩在线播放| 人人妻人人爽人人添夜夜欢视频| 一级片免费观看大全| 国产精品亚洲av一区麻豆| 9191精品国产免费久久| 精品电影一区二区在线| 欧美激情 高清一区二区三区| av视频在线观看入口| 极品人妻少妇av视频| 99精品在免费线老司机午夜| 精品国产亚洲在线| 色综合亚洲欧美另类图片| 色播亚洲综合网| 亚洲欧美精品综合久久99| 国产亚洲精品一区二区www| 日韩欧美国产一区二区入口| АⅤ资源中文在线天堂| 男人的好看免费观看在线视频 | 不卡一级毛片| 日韩精品免费视频一区二区三区| 久久人妻福利社区极品人妻图片| 美女 人体艺术 gogo| 精品卡一卡二卡四卡免费| 伊人久久大香线蕉亚洲五| 国产午夜精品久久久久久| 久9热在线精品视频| 黄色毛片三级朝国网站| 伊人久久大香线蕉亚洲五| 成人亚洲精品一区在线观看| 欧美日韩福利视频一区二区| 国产成人免费无遮挡视频| 午夜精品国产一区二区电影| 免费高清视频大片| 国产av又大| xxx96com| 欧美日韩黄片免| 长腿黑丝高跟| 一区在线观看完整版| 一级毛片精品| 此物有八面人人有两片| 国产成人精品在线电影| x7x7x7水蜜桃| 亚洲av成人一区二区三| 在线观看一区二区三区| a级毛片在线看网站| 国产单亲对白刺激| 久久精品亚洲熟妇少妇任你| 午夜福利18| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 欧美不卡视频在线免费观看 | 91大片在线观看| 在线观看66精品国产| 91麻豆精品激情在线观看国产| 这个男人来自地球电影免费观看| 国产成人av激情在线播放| 少妇粗大呻吟视频| 久久久久久大精品| 十分钟在线观看高清视频www| 欧美日韩瑟瑟在线播放| 国产主播在线观看一区二区| 麻豆av在线久日| 91老司机精品| 欧美+亚洲+日韩+国产| 久久人妻福利社区极品人妻图片| 国产精品亚洲av一区麻豆| 成人国产综合亚洲| 日本 av在线| 色精品久久人妻99蜜桃| 亚洲在线自拍视频| 精品电影一区二区在线| 99国产综合亚洲精品| 狂野欧美激情性xxxx| 国产精品一区二区精品视频观看| 一二三四社区在线视频社区8| 女警被强在线播放| 搞女人的毛片| 精品人妻在线不人妻| 日韩高清综合在线| 丁香六月欧美| 久久久国产欧美日韩av| 俄罗斯特黄特色一大片| 午夜福利成人在线免费观看| 亚洲欧美日韩另类电影网站| 88av欧美| 国产欧美日韩综合在线一区二区| 老司机深夜福利视频在线观看| 午夜福利欧美成人| 视频在线观看一区二区三区| 黑人欧美特级aaaaaa片| 亚洲国产欧美日韩在线播放| 琪琪午夜伦伦电影理论片6080| 热99re8久久精品国产| 手机成人av网站| 亚洲精品美女久久久久99蜜臀| 国产亚洲欧美98| 亚洲午夜精品一区,二区,三区| 亚洲久久久国产精品| 欧美色欧美亚洲另类二区 | 久久中文字幕人妻熟女| 欧美日韩亚洲国产一区二区在线观看| 精品国产亚洲在线| 色综合站精品国产| www国产在线视频色| 黑人巨大精品欧美一区二区mp4| 国产亚洲欧美在线一区二区| 亚洲第一欧美日韩一区二区三区| 精品不卡国产一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 一区在线观看完整版| 搡老熟女国产l中国老女人| av欧美777| 成年人黄色毛片网站| 窝窝影院91人妻| 久久精品91无色码中文字幕| 亚洲在线自拍视频| 国产精品影院久久| 中亚洲国语对白在线视频| 国产成人免费无遮挡视频| 久久狼人影院| 一区在线观看完整版| 国产亚洲精品一区二区www| 91麻豆精品激情在线观看国产| 在线观看免费午夜福利视频| 色婷婷久久久亚洲欧美| 亚洲av美国av| 老司机午夜福利在线观看视频| 一区二区三区激情视频| 91成年电影在线观看| 亚洲欧洲精品一区二区精品久久久| 91字幕亚洲| 久久久久久久久久久久大奶| 黄色毛片三级朝国网站| 男女做爰动态图高潮gif福利片 | 午夜日韩欧美国产| 一区二区三区高清视频在线| 搡老妇女老女人老熟妇| 高清在线国产一区| 不卡一级毛片| 日本a在线网址| 国产99久久九九免费精品| 亚洲成人精品中文字幕电影| 亚洲黑人精品在线| 黄色a级毛片大全视频| 国产野战对白在线观看| 久久性视频一级片| 日韩国内少妇激情av| 99国产精品免费福利视频| 色av中文字幕| 欧美日韩黄片免| 国产成人精品久久二区二区91| 国产精品电影一区二区三区| 在线观看一区二区三区| 在线观看日韩欧美| 中文字幕高清在线视频| 变态另类成人亚洲欧美熟女 | 悠悠久久av| 欧美精品亚洲一区二区| 好看av亚洲va欧美ⅴa在| 天堂影院成人在线观看| 亚洲国产中文字幕在线视频| 久久精品亚洲熟妇少妇任你| 国产精品一区二区三区四区久久 | 日韩欧美免费精品| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| 1024视频免费在线观看| 亚洲午夜理论影院| 亚洲第一青青草原| 男人操女人黄网站| 曰老女人黄片| 琪琪午夜伦伦电影理论片6080| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲精品第一综合不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品亚洲一区二区| 国产精品爽爽va在线观看网站 | 这个男人来自地球电影免费观看| 亚洲国产欧美网| 久久精品国产清高在天天线| 国产成人欧美在线观看| 色哟哟哟哟哟哟| 大型av网站在线播放| 亚洲中文日韩欧美视频| 欧美午夜高清在线| 少妇的丰满在线观看| 国产欧美日韩一区二区三区在线| 黄网站色视频无遮挡免费观看| 国产精品久久久人人做人人爽| 夜夜夜夜夜久久久久| 身体一侧抽搐| 免费人成视频x8x8入口观看| 亚洲精品av麻豆狂野| 免费av毛片视频| 成人国语在线视频| 999久久久精品免费观看国产| 色播亚洲综合网| 久久精品国产综合久久久| 亚洲av成人av| 在线播放国产精品三级| 久久人人精品亚洲av| 国产欧美日韩一区二区精品| 久热爱精品视频在线9| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦观看免费观看视频高清 | 女人被狂操c到高潮| 男人舔女人的私密视频| 午夜福利在线观看吧| 韩国av一区二区三区四区| 99国产精品免费福利视频| 成人18禁在线播放| 国产在线观看jvid| 午夜老司机福利片| 亚洲精品在线美女| 电影成人av| 亚洲自偷自拍图片 自拍| 琪琪午夜伦伦电影理论片6080| 人人澡人人妻人| 一区福利在线观看| 亚洲欧美激情综合另类| 国产99白浆流出| 午夜福利在线观看吧| 亚洲国产精品合色在线| 黑丝袜美女国产一区| 亚洲精品av麻豆狂野| 国产精品久久久人人做人人爽| 中出人妻视频一区二区| 人人妻,人人澡人人爽秒播| 视频在线观看一区二区三区| 不卡av一区二区三区| 大陆偷拍与自拍| 一区二区三区国产精品乱码| 亚洲成av片中文字幕在线观看| 色尼玛亚洲综合影院| 涩涩av久久男人的天堂| 国产精品美女特级片免费视频播放器 | 国产精品美女特级片免费视频播放器 | 久久 成人 亚洲| 久久中文看片网| 国产精品 欧美亚洲| 97人妻精品一区二区三区麻豆 | 黄频高清免费视频| 久久久久久久久中文| 夜夜爽天天搞| 午夜福利影视在线免费观看| 人人妻人人爽人人添夜夜欢视频| 日韩欧美国产在线观看| 亚洲片人在线观看| 亚洲成国产人片在线观看| 好男人电影高清在线观看| 一级,二级,三级黄色视频| 成人永久免费在线观看视频| av有码第一页| 高潮久久久久久久久久久不卡| 久久亚洲精品不卡| 99久久综合精品五月天人人| 精品久久久久久,| 亚洲欧美激情综合另类| 91字幕亚洲| 久99久视频精品免费| 亚洲电影在线观看av| 一个人观看的视频www高清免费观看 | 嫁个100分男人电影在线观看| 精品福利观看| 亚洲av第一区精品v没综合| 看黄色毛片网站| 国产亚洲欧美在线一区二区| 91老司机精品| 日韩欧美三级三区| 伦理电影免费视频| 一区二区日韩欧美中文字幕| 国产成人啪精品午夜网站| 国产真人三级小视频在线观看| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 性少妇av在线| 亚洲精品在线观看二区| 国产麻豆成人av免费视频| 午夜精品国产一区二区电影| 后天国语完整版免费观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲激情在线av| 校园春色视频在线观看| 美女午夜性视频免费| 大码成人一级视频| 日韩成人在线观看一区二区三区| 国内精品久久久久久久电影| 嫩草影视91久久| 18禁黄网站禁片午夜丰满| 伊人久久大香线蕉亚洲五| 欧美人与性动交α欧美精品济南到| 欧美亚洲日本最大视频资源| 高清毛片免费观看视频网站| 国语自产精品视频在线第100页| 久久天躁狠狠躁夜夜2o2o| 在线十欧美十亚洲十日本专区| 久久人妻熟女aⅴ| www日本在线高清视频| 69av精品久久久久久| 色老头精品视频在线观看| av视频在线观看入口| 美女国产高潮福利片在线看| 最近最新中文字幕大全免费视频| 色婷婷久久久亚洲欧美| 精品久久久久久久毛片微露脸| 午夜日韩欧美国产| av在线天堂中文字幕| 国产精品亚洲美女久久久| 免费女性裸体啪啪无遮挡网站| 91在线观看av| 亚洲九九香蕉| 欧美乱色亚洲激情| 日本免费a在线| 咕卡用的链子| 午夜福利免费观看在线| 丝袜美足系列| 国内毛片毛片毛片毛片毛片| 亚洲精品av麻豆狂野| 欧美日本中文国产一区发布| 久久青草综合色| 国产精品日韩av在线免费观看 | 色综合亚洲欧美另类图片| 亚洲精品国产色婷婷电影| 他把我摸到了高潮在线观看| xxx96com| 国产一区二区激情短视频| 黑人操中国人逼视频| 日日干狠狠操夜夜爽| 777久久人妻少妇嫩草av网站| 啪啪无遮挡十八禁网站| 国产精品一区二区精品视频观看| 久久天躁狠狠躁夜夜2o2o| 亚洲精品在线美女| 99精品久久久久人妻精品| 国产精品久久久久久亚洲av鲁大| 精品无人区乱码1区二区| 女人高潮潮喷娇喘18禁视频| 成年女人毛片免费观看观看9| 亚洲片人在线观看| 制服人妻中文乱码| 99在线人妻在线中文字幕| 在线观看舔阴道视频| 啪啪无遮挡十八禁网站| 国产99白浆流出| 97超级碰碰碰精品色视频在线观看| 男女下面进入的视频免费午夜 | 三级毛片av免费| 久久国产乱子伦精品免费另类| 精品欧美一区二区三区在线| 午夜a级毛片| 久久人妻av系列| 久久香蕉国产精品| 国产伦一二天堂av在线观看| 亚洲国产精品成人综合色| 亚洲午夜精品一区,二区,三区| 黄片大片在线免费观看| 激情视频va一区二区三区| 午夜精品久久久久久毛片777|