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

    產(chǎn)茸相關(guān)SNP位點在7個梅花鹿群體中的驗證

    2021-10-14 07:19:52田雪琪李洋劉匯濤王天驕邢秀梅
    特產(chǎn)研究 2021年5期
    關(guān)鍵詞:梅花鹿雜合多態(tài)性

    田雪琪,李洋,劉匯濤,王天驕,邢秀梅

    (中國農(nóng)業(yè)科學院特產(chǎn)研究所,特種經(jīng)濟動物分子生物重點實驗室,吉林 長春130112)

    2020年5月29日,以梅花鹿等16種特種畜禽列入《國家畜禽遺傳資源目錄》[1]。鹿茸作為傳統(tǒng)的中藥材在中國已有兩千多年的應(yīng)用歷史[2],產(chǎn)茸性能一直是評估梅花鹿優(yōu)劣的重要指標。梅花鹿的選種育種一直停留在常規(guī)表型選種階段,茸重、茸型、茸色和對稱性是梅花鹿茸用品種培育的主要選種依據(jù)。分子標記輔助選種仍處于起步階段,田萬年等[3]對雙陽梅花鹿胰島素樣生長因子-1(IGF-1)基因的研究,發(fā)現(xiàn)在824 bp處存在T→C突變,所產(chǎn)生的AB基因型個體的產(chǎn)茸性能顯著高于其他兩種基因型個體;杜智恒等[4]對東北梅花鹿生長激素基因(GH)進行研究,發(fā)現(xiàn)在第2內(nèi)含子中存在G→A突變,這個SNP位點不同基因型可能對梅花鹿產(chǎn)茸性狀有影響,尤其是對第5鋸的產(chǎn)茸量有顯著影響趨勢;熊家軍[5]以兩個地區(qū)梅花鹿種群作為研究對象,對褪黑素基因(MTNR1A)、生長激素基因(GH)以及雄激素受體基因(AR)的SNPs與產(chǎn)茸量相關(guān)性進行了分析,發(fā)現(xiàn)在MTNR1A第2外顯子中存在G→C的位點突變,對湖北荊門五三梅花鹿種群產(chǎn)茸量具有顯著影響;GH基因139 bp處存在C→T突變對吉林東豐地區(qū)梅花鹿種群產(chǎn)茸量具有顯著影響;AR第3外顯子存在的C→T突變對兩個種群的產(chǎn)茸量均有顯著影響;陳斌[6]對梅花鹿轉(zhuǎn)化生長因子-1基因(TGF-1)的研究,發(fā)現(xiàn)在371 bp處的SNP位點存在3種基因型,其中AA基因型與其他兩種基因型的產(chǎn)茸量存在顯著差異。由此可見,在IGF-1、GH、MTNR1A、AR、TGF-1基因上均存在著與產(chǎn)茸量相關(guān)的SNP位點,且不同群體中與產(chǎn)茸相關(guān)的SNP位點存在一定差異,見表1。前期實驗室在四平梅花鹿、東大梅花鹿兩個群體篩選得到8個與產(chǎn)茸性能相關(guān)的SNP位點。其中,SNP1427_255為極顯著相關(guān),其余SNP位點為顯著相關(guān)。本試驗應(yīng)用KASP分型技術(shù)在7個梅花鹿群體中進行基因分型,并對不同基因型與產(chǎn)茸量做關(guān)聯(lián)分析,以探究東北梅花鹿產(chǎn)茸量相關(guān)基因型或SNP組合,為早期選種提供依據(jù)。

    表1 各研究中與產(chǎn)茸相關(guān)SNP位點信息Table 1 The information of SNP locus in studies

    1 試驗材料和方法

    1.1 試驗材料

    試驗用雄性梅花鹿個體共計190頭,分別采集于吉林省遼源市東豐縣文福鹿場、吉林省通化市通化縣通化山寶鹿業(yè)、吉林省長春市雙陽區(qū)東鰲鹿業(yè)、吉林省遼源市東豐縣鑫利達鹿業(yè)、吉林省梅河口市大成加隆畜禽有限公司、吉林省四平市吉梅花鹿種鹿場和吉林省長春市雙陽區(qū)鹿業(yè)良種繁育有限公司。所有試驗個體使用抗凝真空采血管采集頸靜脈血約10 mL,帶回實驗室后﹣20℃保存,用于后續(xù)基因組DNA提取。同時記錄所有試驗個體的原始產(chǎn)茸記錄,用于后續(xù)產(chǎn)茸性能關(guān)聯(lián)分析。

    1.2 試驗方法

    1.2.1 基因組DNA的提取使用天根生化科技(北京)有限公司所生產(chǎn)的血液/細胞/組織基因組DNA提取試劑盒(離心柱型)提取試驗樣品的基因組DNA,相關(guān)操作參照試劑盒說明書。用1.0%的瓊脂糖凝膠電泳檢測DNA的完整性;用NanoDrop分光光度計測DNA的濃度,確定基因組DNA質(zhì)量合格后,置于﹣20℃保存?zhèn)溆谩?/p>

    表2 試驗樣品Table 2 The samples submitted to this study

    1.2.2 競爭性等位基因特異性PCR分型技術(shù)本實驗室提供8個SNP位點的序列信息,后續(xù)KASP分型實驗中引物設(shè)計、合成及位點檢測均由北京康普森公司完成。該技術(shù)是根據(jù)對終端熒光信號的讀取對位點進行判斷,每孔反應(yīng)都是采用雙色熒光檢測1個SNP位點的2種基因型,不同的SNP對應(yīng)著不同的熒光信號。通過掃描試驗用多重PCR板,采集不同的熒光信號得出檢測位點的具體信息。

    1.2.3 產(chǎn)茸性能估測在實際生產(chǎn)中,由于鹿場的管理模式、養(yǎng)殖人員更換等因素,可能會造成某些個體的產(chǎn)茸記錄部分缺失。李和平等[7]建立了根據(jù)整個群體和個體產(chǎn)茸記錄來準確估測公鹿產(chǎn)茸能力的數(shù)學模型,這個模型中假設(shè)存在一個包含若干頭公鹿的群體,其中第i頭公鹿有Ki次產(chǎn)茸記錄,則這頭公鹿i最有可能的產(chǎn)茸能力(Pi)為

    其中P指鹿群鋸平均產(chǎn)茸量;t指同一公鹿不同鋸別產(chǎn)茸記錄之間的相關(guān)(重復(fù)力);Xi指公鹿的鋸平均產(chǎn)茸量。本試驗根據(jù)這個模型對試驗個體進行產(chǎn)茸能力估測。

    1.2.4 SNP位點統(tǒng)計分析根據(jù)群體遺傳學的研究方法,運用Cervus軟件對群體的遺傳參數(shù)進行計算,其中包括基因頻率、基因型頻率、觀測雜合度(Ho)、期望雜合度(He)、有效等位基因數(shù)(Ne)、多態(tài)信息含量(PIC)和哈迪-溫伯格(H-W)平衡檢驗等[8],同時運用SPSS19.0軟件,采用一般線性模型對遺傳信息與產(chǎn)茸性能的關(guān)聯(lián)性進行最小二乘法分析。

    2 結(jié)果

    2.1 基因組DNA提取

    將所提取的基因組DNA經(jīng)1%瓊脂糖凝膠電泳和NanoDrop2000進行基因組DNA完整性、純度和濃度檢測,結(jié)果見圖1,電泳條帶單一且清晰明亮,大小位于15 000 kb Marker頂端條帶之上;Nano-Drop2000檢測結(jié)果顯示,所提取的DNA濃度在50~100ng/L,OD值介于1.8~2.0之間,紫外吸收峰峰型正常。

    圖1 基因組DNA瓊脂糖凝膠電泳圖Fig.1 Genomic DNA agarose gel electrophoresis

    2.2 引物設(shè)計結(jié)果及相關(guān)信息

    應(yīng)用Primer 5軟件,對8個位點所在序列進行分析,經(jīng)過篩選得出最適引物。本試驗中所用8個SNP位點的引物序列及相關(guān)信息見表3。

    表3 SNP位點引物信息Table 3 The information of SNP’s primer

    2.3 SNP位點分型結(jié)果

    運用KASP技術(shù),對190個樣本的8個SNP位點進行分型。研究表明,用于表型關(guān)聯(lián)分析的位點,基因型檢出率不低于94%[9]。本研究對8個SNP位點檢出率情況進行統(tǒng)計。結(jié)果顯示,位點的檢出率介于95.26%~98.42%,平均檢出率為97.43%,均大于94%,全部分型成功。同時還對這些位點進行最小等位基因頻率(MAF)計算分析,平均MAF值為0.188,結(jié)果見表4。

    表4 SNP位點MAF值計算結(jié)果Table 4 The calculated results of MAF value of SNP locus

    2.4 群體遺傳多樣性分析

    運用KASP對8個SNP位點在驗證群體進行基因分型,所得到的基因型信息使用Cervus軟件計算。通過得到的基因型信息對該群體的遺傳多樣性進行分析,結(jié)果見表5。8個SNP位點中,有7個位點具有多態(tài)性,其多態(tài)信息含量(PIC)介于0.078~0.369,平均值為0.231 6;最小基因等位基因頻率(MAF)介于0.044~0.396之間,平均值為0.188;各位點的觀測雜合度(Ho)介于0.090~0.426,平均值為0.265;期望雜合度(He)為0.082~0.490,平均值為0.282;有效等位基因數(shù)(Ne)介于1.2~2.8,平均值為1.9。觀測雜合度與期望雜合度不存在統(tǒng)計學上的顯著差異(P>0.05)。Hardy-Weinberg平衡分析結(jié)果顯示,SNP614_2585、SNP1881_130和SNP596_3774這3個SNP位點顯著偏離平衡(P<0.05)。

    表5 群體中9個SNP位點的遺傳參數(shù)信息Table 5 The genetic parameter information of 9 SNP locus in the population

    2.5 產(chǎn)茸性能估測結(jié)果

    運用茸重估測模型,計算驗證群體中每個個體的產(chǎn)茸估測值。計算結(jié)果顯示,試驗群體的產(chǎn)茸估測值介于1.12~8.55 kg。對所有個體的產(chǎn)茸估測值進行統(tǒng)計分析,根據(jù)相關(guān)研究,將產(chǎn)茸估測值分大于3.5 kg的高產(chǎn)組和小于2.5 kg的低產(chǎn)組[10]。對兩組數(shù)據(jù)進行單因素方差分析,結(jié)果顯示兩組產(chǎn)茸估測數(shù)據(jù)存在極顯著差異(P<0.01),證明本群體中存在高、低產(chǎn)個體。部分樣本的產(chǎn)茸記錄和估測值,見表6。

    表6 部分梅花鹿產(chǎn)茸性能估測Table 6 The estimation of the antler performance

    2.6 SNP位點與產(chǎn)茸性狀的關(guān)聯(lián)分析

    對190頭梅花鹿進行產(chǎn)茸性狀與SNP標記相關(guān)性分析。經(jīng)過計算,7個SNP位點的相關(guān)性P值介于0.264~0.849之間,均大于0.05,顯示所有位點SNP位點與產(chǎn)茸性狀均不存在顯著相關(guān)性(P>0.05),結(jié)果見表7。

    表7 SNP位點關(guān)聯(lián)性分析結(jié)果Table 7 The phenotypic interpretation rate of SNP Locus

    根據(jù)檢測結(jié)果顯示,位點SNP219_9123所有個體基因型均為CC,無多態(tài)性。SNP1427_255、SNP307_1564、SNP123_7847、SNP2027_423、SNP614_2585和SNP596_3774均存在3種基因型,但各基因型所占比例不同。SNP1881_130中僅存在兩種基因型,CC和CT,其中CC為優(yōu)勢基因型,基因型頻率為0.91。SNP1427_255中,純合CC為優(yōu)勢基因型,基因型頻率為0.56。SNP307_1564中,純合GG為優(yōu)勢基因型,基因型頻率為0.63。SNP123_7847中,純合TT為優(yōu)勢基因型,基因型頻率為0.56。SNP2027_423中,雜合AG為優(yōu)勢基因型,基因型頻率為0.43。SNP614_2585中,純合AA為優(yōu)勢基因型,基因型頻率為0.87。SNP596_3774中,純合TT為優(yōu)勢基因型,基因型頻率為0.82。對每個SNP位點不同基因型的產(chǎn)茸性估測值均值進行單因素方差分析。結(jié)果顯示,各位點的不同基因型間產(chǎn)茸估測值均值不存在顯著差異。進一步驗證了8個SNP位點在本驗證群中與產(chǎn)茸性能不存在顯著相關(guān)性。但在SNP307_156、SNP614_258和SNP596_377位點具有明顯的高產(chǎn)基因型。在SNP307_156中,存在與已知研究相同的高產(chǎn)CC基因型[11],分別高于CG基因型和GG基因型14.36%和17.85%。在該位點中,對具有CC基因型個體的產(chǎn)茸值與CG基因型和GG基因型個體進行比較,發(fā)現(xiàn)兩組產(chǎn)茸量具有顯著差異性(P<0.05)。還發(fā)現(xiàn)SNP614_258位點高產(chǎn)基因型為GG,SNP596_377位點高產(chǎn)基因型為CC。具體位點詳細信息見表8。

    表8 多態(tài)性位點的基因型信息Table 8 The genotype information of the polymorphic locus

    3 討論

    根據(jù)群體遺傳學的研究方法,運用Cervus軟件計算多態(tài)信息含量(PIC)、平均期望雜合度(He)、觀測雜合度(Ho)和Hardy-Weinberg群體平衡,以分析群體遺傳多樣性。多態(tài)信息含量這一指標是由Botsticn等[12]在1980年提出,該參數(shù)的高低用以描述家系連鎖分析標志基因(或標志序列)的多態(tài)性大小。Botsticn認為,當多態(tài)信息含量大于0.5時,該標志基因具有高度信息含量,小于0.5而大于0.25則為中度信息含量,小于0.25時具有低度信息含量。本試驗中所驗證的7個SNP位點的多態(tài)信息含量(PIC)介于0.078~0.369,其中SNP614_2585、SNP1881_130和SNP596_3774為低度多態(tài)性位點,SNP1427_255、SNP307_1564、SNP123_7847和SNP2027_423為中度多態(tài)性位點,沒有高度多態(tài)性位點。相關(guān)研究表明,與微衛(wèi)星標記相比,SNP的多態(tài)性較低,但因SNP位點為全基因組層面內(nèi)分布,可在一定程度上彌補其相較于微衛(wèi)星等其他分子標記方法中多態(tài)性較低的問題,使得SNP位點可用于相關(guān)性狀分析[13]。

    雜合度用來評價群體在檢測位點上出現(xiàn)雜合的頻率,是評價群體遺傳變異的重要參數(shù)。通過對比觀測雜合度(Ho)和期望雜合度(He)的相對大小來判定群體近交程度。當觀測雜合度小于期望雜合度時,群體遺傳多樣性低,為近親繁殖;當兩者相等時,群體內(nèi)個體為隨機交配;當觀測雜合度大于期望雜合度,群體遺傳多樣性較高,為遠親繁殖。本試驗群體的平均觀測雜合度(Ho)為0.265,平均期望雜合度(He)為0.282,呈觀測雜合度略小于期望雜合度,表明群體雜合度低,群體遺傳多樣性低。

    Hardy-Weinberg平衡是以完全隨機交配所組成的自然種群為基礎(chǔ)而做出的假設(shè)。因此,通過Hardy-Weinberg群體平衡卡方檢測來衡量群體遺傳平衡狀態(tài)。本研究群體經(jīng)過Hardy-Weinberg群體平衡檢測發(fā)現(xiàn),在具有多態(tài)性的7個位點中,有3個位點顯著偏離平衡狀態(tài)?,F(xiàn)有研究結(jié)果顯示,在小規(guī)模種群中,近親交配和無效等位基因等因素會導致雜合度不足,造成位點偏離Hardy-Weinberg平衡[14]。所以這3個顯著偏離Hardy-Weinberg平衡的位點可能與群體遺傳多樣性低、人工選育干預(yù)強度大的情況有關(guān)。東北地區(qū)的家養(yǎng)梅花鹿均為野生東北梅花鹿的后裔,結(jié)合本研究群體的多個遺傳信息結(jié)果可以看出,東北地區(qū)各梅花鹿種群間的遺傳背景十分相似。

    本試驗是對已知與產(chǎn)茸量相關(guān)的SNP位點進行驗證。在試驗群體中,這些SNP位點與本群體中的產(chǎn)茸量不存在顯著相關(guān)性。推測其原因,可能是由于在篩選群體中,主要是以四平梅花鹿和東大梅花鹿為主。這兩個品種的梅花鹿在已知品種中為高產(chǎn)品種。特別是東大梅花鹿群體,在產(chǎn)茸性能方面極顯著高于四平梅花鹿和其他群體。據(jù)統(tǒng)計,東大梅花鹿和四平梅花鹿的上鋸三杈鮮茸平均單產(chǎn)為4 kg和3.42 kg,明顯高于其他品種,如雙陽梅花鹿上鋸三杈鮮茸單產(chǎn)為2.9 kg[15],敖東梅花鹿上鋸三杈鮮茸單產(chǎn)為3.34kg[16]。究其原因可能與東大梅花鹿在育種期間引入馬鹿血源及高強度人工選擇有關(guān)[11]。本研究中的190頭梅花鹿群體經(jīng)全基因組SNP鑒別芯片檢測,判定為純種梅花鹿,純種梅花鹿在產(chǎn)茸量方面較雜交個體相差甚遠。因此,我們推測群體遺傳背景差異可能是造成與產(chǎn)茸相關(guān)性SNP位點在本試驗中不相關(guān)的主要原因。后續(xù)對梅花鹿產(chǎn)茸性能相關(guān)的SNP位點的篩選,建議擴大篩選群體,增加群體數(shù)量。在表型數(shù)據(jù)中,除依據(jù)產(chǎn)茸估測之外,還可以通過相同鋸齡的實際產(chǎn)茸量進行關(guān)聯(lián)分析,為吉林梅花鹿的產(chǎn)茸性能SNP位點提供有效參考。

    4 結(jié)論

    本研究運用競爭性等位基因特異性PCR分型技術(shù)對190個個體的8個SNP位點進行基因分型檢測,得出7個具有多態(tài)性的SNP位點。使用一般線性模型對這7個SNP位點遺傳信息與產(chǎn)茸量進行關(guān)聯(lián)分析,發(fā)現(xiàn)在本試驗中所建立的隨機驗證群體中,所有SNP位點與產(chǎn)茸量均不存在顯著關(guān)聯(lián),但存在SNP307_156、SNP614_258和SNP596_377位點具有明顯的高產(chǎn)基因型。本研究表明,前期基于部分群體篩選出的與產(chǎn)茸量相關(guān)SNP位點在本試驗群體的適用性不理想。

    猜你喜歡
    梅花鹿雜合多態(tài)性
    梅花鹿
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進展
    甘藍型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    貪玩的小梅花鹿
    貪玩的小梅花鹿
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    馬鈴薯cpDNA/mtDNA多態(tài)性的多重PCR檢測
    GlobalFiler~? PCR擴增試劑盒驗證及其STR遺傳多態(tài)性
    CYP3A4*1G基因多態(tài)性及功能的初步探討
    雄激素可調(diào)節(jié)的腎臟近端腎小管上皮細胞靶向雜合啟動子的優(yōu)化
    欧美人与善性xxx| 国产av国产精品国产| 少妇人妻精品综合一区二区| 成年女人在线观看亚洲视频| 国产熟女欧美一区二区| 精品国产乱码久久久久久小说| 国产亚洲最大av| 婷婷色综合www| 中文字幕制服av| 少妇的逼水好多| 亚洲av日韩在线播放| 日韩制服骚丝袜av| 一本久久精品| 丝瓜视频免费看黄片| 深夜精品福利| 日本色播在线视频| 国产日韩一区二区三区精品不卡| 久久人人爽人人片av| 天天影视国产精品| 成人毛片60女人毛片免费| 日韩大片免费观看网站| 一本久久精品| 少妇精品久久久久久久| 老司机影院毛片| 99香蕉大伊视频| 黄片播放在线免费| 欧美 日韩 精品 国产| 免费在线观看完整版高清| 亚洲精品视频女| 有码 亚洲区| 男人爽女人下面视频在线观看| 午夜福利视频在线观看免费| 男女免费视频国产| 欧美国产精品一级二级三级| 三级国产精品片| 免费看不卡的av| 久久免费观看电影| 成人漫画全彩无遮挡| 少妇熟女欧美另类| 18禁国产床啪视频网站| 欧美精品一区二区免费开放| 天堂8中文在线网| 99热这里只有是精品在线观看| 热99久久久久精品小说推荐| 成人手机av| 男女国产视频网站| 亚洲成人一二三区av| 人妻一区二区av| 国产精品偷伦视频观看了| 伦理电影免费视频| 久久久久久久久久久久大奶| 夫妻午夜视频| 水蜜桃什么品种好| 午夜视频国产福利| 这个男人来自地球电影免费观看 | 婷婷色综合大香蕉| 91在线精品国自产拍蜜月| 亚洲精品456在线播放app| 久久亚洲国产成人精品v| 亚洲av中文av极速乱| √禁漫天堂资源中文www| 久久久久网色| 亚洲丝袜综合中文字幕| 国产极品天堂在线| 国产av国产精品国产| 在线亚洲精品国产二区图片欧美| 伦理电影大哥的女人| 街头女战士在线观看网站| 日本午夜av视频| 国产精品三级大全| 亚洲,一卡二卡三卡| 黄色 视频免费看| 十八禁网站网址无遮挡| 亚洲国产av新网站| 久久 成人 亚洲| 久久久久久久国产电影| 熟女电影av网| 国产精品欧美亚洲77777| 欧美日韩av久久| 久久这里有精品视频免费| 欧美精品一区二区大全| 美女xxoo啪啪120秒动态图| 天堂俺去俺来也www色官网| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜日本视频在线| 多毛熟女@视频| 欧美xxxx性猛交bbbb| 午夜免费男女啪啪视频观看| 最黄视频免费看| 一二三四中文在线观看免费高清| 国产亚洲精品第一综合不卡 | 99热6这里只有精品| 亚洲精华国产精华液的使用体验| 美女主播在线视频| 精品卡一卡二卡四卡免费| 涩涩av久久男人的天堂| 欧美日韩一区二区视频在线观看视频在线| 日韩三级伦理在线观看| 人妻系列 视频| 久久久久精品人妻al黑| 国产极品粉嫩免费观看在线| 99久国产av精品国产电影| 亚洲国产看品久久| 制服诱惑二区| 美国免费a级毛片| 国产又爽黄色视频| 中文欧美无线码| 一区二区三区四区激情视频| 亚洲伊人色综图| 香蕉丝袜av| 中文精品一卡2卡3卡4更新| 国产片特级美女逼逼视频| 欧美变态另类bdsm刘玥| 久久久久精品久久久久真实原创| 我要看黄色一级片免费的| 日韩视频在线欧美| 如日韩欧美国产精品一区二区三区| 啦啦啦中文免费视频观看日本| 天天躁夜夜躁狠狠躁躁| 大话2 男鬼变身卡| 国产高清国产精品国产三级| 日韩制服丝袜自拍偷拍| 大话2 男鬼变身卡| 亚洲国产欧美在线一区| 18禁在线无遮挡免费观看视频| 另类精品久久| 巨乳人妻的诱惑在线观看| 国产高清不卡午夜福利| 国产亚洲午夜精品一区二区久久| 国产日韩欧美视频二区| 亚洲经典国产精华液单| 九色成人免费人妻av| a级毛色黄片| www.av在线官网国产| 久久久久久伊人网av| 中文乱码字字幕精品一区二区三区| 少妇被粗大的猛进出69影院 | 日韩 亚洲 欧美在线| 亚洲性久久影院| 国产精品秋霞免费鲁丝片| 夜夜骑夜夜射夜夜干| 国产淫语在线视频| 久久精品久久久久久噜噜老黄| 婷婷色综合www| av视频免费观看在线观看| 如日韩欧美国产精品一区二区三区| av卡一久久| 欧美精品高潮呻吟av久久| 日韩av不卡免费在线播放| 色视频在线一区二区三区| 十八禁高潮呻吟视频| 制服丝袜香蕉在线| 久久青草综合色| av黄色大香蕉| 欧美亚洲 丝袜 人妻 在线| 久久 成人 亚洲| 一级爰片在线观看| 日韩中字成人| 人妻少妇偷人精品九色| 成人亚洲精品一区在线观看| 91久久精品国产一区二区三区| 晚上一个人看的免费电影| 欧美日韩综合久久久久久| 午夜激情久久久久久久| 韩国高清视频一区二区三区| 午夜福利在线观看免费完整高清在| 啦啦啦在线观看免费高清www| 一级黄片播放器| 黄色配什么色好看| 国国产精品蜜臀av免费| 婷婷成人精品国产| 你懂的网址亚洲精品在线观看| 男女国产视频网站| 国产日韩欧美亚洲二区| 国产黄频视频在线观看| 99久国产av精品国产电影| 尾随美女入室| 欧美精品国产亚洲| 精品99又大又爽又粗少妇毛片| 肉色欧美久久久久久久蜜桃| 伦理电影免费视频| 亚洲欧美色中文字幕在线| 久久99蜜桃精品久久| 美女内射精品一级片tv| 午夜福利,免费看| www.色视频.com| 少妇 在线观看| 岛国毛片在线播放| 国产片内射在线| 亚洲图色成人| 亚洲精品视频女| 亚洲国产日韩一区二区| 大片电影免费在线观看免费| 久久精品国产亚洲av天美| 人成视频在线观看免费观看| av一本久久久久| 内地一区二区视频在线| 欧美日韩成人在线一区二区| 99久国产av精品国产电影| 高清视频免费观看一区二区| 久久人人爽人人爽人人片va| 伊人久久国产一区二区| 国产熟女午夜一区二区三区| 国产 精品1| av有码第一页| 亚洲情色 制服丝袜| 高清在线视频一区二区三区| 欧美日韩视频精品一区| 热re99久久国产66热| 免费人成在线观看视频色| 最黄视频免费看| 下体分泌物呈黄色| 久久久久久久久久人人人人人人| 狂野欧美激情性xxxx在线观看| 高清毛片免费看| 美女大奶头黄色视频| 久久 成人 亚洲| 成人午夜精彩视频在线观看| 久久久久精品久久久久真实原创| 性色av一级| 九九在线视频观看精品| 啦啦啦中文免费视频观看日本| 亚洲三级黄色毛片| 国产欧美日韩综合在线一区二区| a级毛色黄片| 免费观看在线日韩| 91精品国产国语对白视频| 久久婷婷青草| 黑人高潮一二区| 久久久久久久大尺度免费视频| 满18在线观看网站| 精品久久久精品久久久| 久久久久精品人妻al黑| 日本猛色少妇xxxxx猛交久久| 高清毛片免费看| 日产精品乱码卡一卡2卡三| 一级毛片我不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产亚洲欧美精品永久| 国产日韩欧美亚洲二区| 国产国语露脸激情在线看| 国产精品嫩草影院av在线观看| 捣出白浆h1v1| 青春草亚洲视频在线观看| 涩涩av久久男人的天堂| 精品国产一区二区久久| 亚洲精品456在线播放app| 伊人亚洲综合成人网| 欧美日韩国产mv在线观看视频| 97超碰精品成人国产| 一本色道久久久久久精品综合| 亚洲成人一二三区av| 欧美精品一区二区大全| 亚洲av免费高清在线观看| 岛国毛片在线播放| 久久ye,这里只有精品| 涩涩av久久男人的天堂| 亚洲国产av新网站| 免费黄网站久久成人精品| 欧美亚洲 丝袜 人妻 在线| 99久国产av精品国产电影| 精品酒店卫生间| 激情五月婷婷亚洲| 久久婷婷青草| 午夜激情久久久久久久| 日日啪夜夜爽| 99热国产这里只有精品6| 韩国av在线不卡| 午夜福利在线观看免费完整高清在| 一二三四中文在线观看免费高清| 久久亚洲国产成人精品v| 亚洲综合色网址| 熟妇人妻不卡中文字幕| 91aial.com中文字幕在线观看| 9热在线视频观看99| 欧美最新免费一区二区三区| 国精品久久久久久国模美| 考比视频在线观看| 夫妻午夜视频| 秋霞伦理黄片| 十分钟在线观看高清视频www| 亚洲av综合色区一区| 22中文网久久字幕| 在线 av 中文字幕| 亚洲第一区二区三区不卡| 亚洲,一卡二卡三卡| 蜜桃在线观看..| 国产男人的电影天堂91| 精品亚洲乱码少妇综合久久| 国产麻豆69| 亚洲国产欧美日韩在线播放| 国产黄频视频在线观看| 亚洲五月色婷婷综合| 制服丝袜香蕉在线| av不卡在线播放| 少妇高潮的动态图| 视频中文字幕在线观看| 欧美+日韩+精品| 夫妻午夜视频| 考比视频在线观看| 男女高潮啪啪啪动态图| 一区二区日韩欧美中文字幕 | 精品亚洲乱码少妇综合久久| 草草在线视频免费看| 啦啦啦啦在线视频资源| 免费高清在线观看视频在线观看| 韩国av在线不卡| 人成视频在线观看免费观看| 精品一区二区三区视频在线| 亚洲av日韩在线播放| 一本久久精品| av.在线天堂| 嫩草影院入口| av有码第一页| av卡一久久| 久久久久国产网址| 自拍欧美九色日韩亚洲蝌蚪91| 少妇熟女欧美另类| 婷婷色麻豆天堂久久| 久久 成人 亚洲| 一级爰片在线观看| 日日撸夜夜添| 国产亚洲av片在线观看秒播厂| 26uuu在线亚洲综合色| 啦啦啦啦在线视频资源| 亚洲色图 男人天堂 中文字幕 | 高清黄色对白视频在线免费看| 免费大片18禁| 极品少妇高潮喷水抽搐| 日本免费在线观看一区| 午夜激情久久久久久久| 99热这里只有是精品在线观看| 99国产综合亚洲精品| 欧美精品av麻豆av| 如何舔出高潮| 男女边摸边吃奶| 亚洲国产精品成人久久小说| 一级片免费观看大全| 欧美国产精品va在线观看不卡| 成人亚洲欧美一区二区av| 日日摸夜夜添夜夜爱| 亚洲成av片中文字幕在线观看 | 大片免费播放器 马上看| 国产在线免费精品| 精品少妇内射三级| 看十八女毛片水多多多| 王馨瑶露胸无遮挡在线观看| 国产白丝娇喘喷水9色精品| 一区二区三区四区激情视频| 亚洲中文av在线| 亚洲精品久久午夜乱码| 91精品三级在线观看| 中国国产av一级| 18禁国产床啪视频网站| 综合色丁香网| 国产在线免费精品| 亚洲高清免费不卡视频| 99国产综合亚洲精品| 丰满少妇做爰视频| 国产有黄有色有爽视频| 久久国内精品自在自线图片| 国产成人精品久久久久久| 国产精品久久久久成人av| 欧美日韩综合久久久久久| 你懂的网址亚洲精品在线观看| 日韩欧美精品免费久久| 欧美xxⅹ黑人| a级毛片黄视频| 男女免费视频国产| 国产黄色视频一区二区在线观看| 91午夜精品亚洲一区二区三区| 51国产日韩欧美| 在线看a的网站| 美国免费a级毛片| 国产熟女欧美一区二区| 欧美国产精品va在线观看不卡| 亚洲欧美一区二区三区黑人 | 亚洲欧美日韩卡通动漫| 国产日韩欧美亚洲二区| 我的女老师完整版在线观看| 欧美精品高潮呻吟av久久| 亚洲精品色激情综合| 国产精品久久久av美女十八| 九色成人免费人妻av| 午夜免费观看性视频| 有码 亚洲区| 午夜视频国产福利| 黄片播放在线免费| 国产精品.久久久| 卡戴珊不雅视频在线播放| 国产免费现黄频在线看| 夜夜骑夜夜射夜夜干| 午夜福利在线观看免费完整高清在| 欧美国产精品va在线观看不卡| 最近最新中文字幕大全免费视频 | 观看av在线不卡| 如何舔出高潮| 色婷婷av一区二区三区视频| 国产又爽黄色视频| 九色成人免费人妻av| 久久久精品免费免费高清| 夫妻午夜视频| 欧美人与善性xxx| 亚洲 欧美一区二区三区| 伦理电影大哥的女人| 99久久综合免费| 人人妻人人澡人人爽人人夜夜| 精品国产国语对白av| 成人亚洲欧美一区二区av| 春色校园在线视频观看| 成人毛片60女人毛片免费| 欧美精品亚洲一区二区| 久久99蜜桃精品久久| 如何舔出高潮| 女人精品久久久久毛片| 久久狼人影院| 久久午夜综合久久蜜桃| 黄色一级大片看看| 日韩欧美精品免费久久| 中文天堂在线官网| 欧美3d第一页| 久久这里只有精品19| 在线观看人妻少妇| 欧美国产精品一级二级三级| 99国产精品免费福利视频| 久久精品久久精品一区二区三区| www.av在线官网国产| 精品国产国语对白av| 欧美xxⅹ黑人| 国产日韩欧美视频二区| 99热国产这里只有精品6| 亚洲少妇的诱惑av| 一边亲一边摸免费视频| 国产精品不卡视频一区二区| 久久精品久久精品一区二区三区| 99国产综合亚洲精品| 色婷婷久久久亚洲欧美| 超色免费av| 如何舔出高潮| 18禁裸乳无遮挡动漫免费视频| 亚洲三级黄色毛片| 9色porny在线观看| 国产精品.久久久| 春色校园在线视频观看| 亚洲一区二区三区欧美精品| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 在线观看免费视频网站a站| 久久影院123| 国产在线视频一区二区| 国产在视频线精品| 亚洲av中文av极速乱| 国产精品久久久久久久久免| 我的女老师完整版在线观看| 欧美变态另类bdsm刘玥| 亚洲欧美日韩卡通动漫| 天堂中文最新版在线下载| 最近的中文字幕免费完整| 日本与韩国留学比较| 内地一区二区视频在线| 亚洲av福利一区| 宅男免费午夜| 啦啦啦中文免费视频观看日本| 久久久久精品人妻al黑| 久久人人爽av亚洲精品天堂| 赤兔流量卡办理| 久久久久精品久久久久真实原创| 大片免费播放器 马上看| 秋霞伦理黄片| 日日撸夜夜添| 免费黄频网站在线观看国产| 美女国产视频在线观看| 在线观看人妻少妇| 黑人高潮一二区| 国产成人精品一,二区| 亚洲国产精品专区欧美| 久久精品国产亚洲av天美| 欧美精品一区二区免费开放| 亚洲国产精品成人久久小说| 国产成人精品久久久久久| 在线观看一区二区三区激情| 2022亚洲国产成人精品| 亚洲人与动物交配视频| 这个男人来自地球电影免费观看 | 大码成人一级视频| 91成人精品电影| 国产精品秋霞免费鲁丝片| 精品99又大又爽又粗少妇毛片| 欧美日韩精品成人综合77777| 午夜精品国产一区二区电影| 免费高清在线观看视频在线观看| 日韩 亚洲 欧美在线| a级毛片在线看网站| 国产亚洲精品久久久com| 制服丝袜香蕉在线| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久久免费av| 久久韩国三级中文字幕| 国产一区二区激情短视频 | 黄色怎么调成土黄色| 看免费成人av毛片| 亚洲色图综合在线观看| 日韩,欧美,国产一区二区三区| 高清不卡的av网站| av有码第一页| 精品久久久精品久久久| 国产免费福利视频在线观看| av国产精品久久久久影院| 日韩,欧美,国产一区二区三区| 嫩草影院入口| 欧美日韩一区二区视频在线观看视频在线| 大香蕉97超碰在线| 欧美精品人与动牲交sv欧美| 亚洲欧美日韩另类电影网站| 香蕉国产在线看| 中文字幕制服av| 亚洲伊人久久精品综合| 在线观看一区二区三区激情| av一本久久久久| 不卡视频在线观看欧美| 欧美性感艳星| 国产成人一区二区在线| 亚洲精品久久午夜乱码| 久久精品国产亚洲av涩爱| 亚洲精品久久午夜乱码| 国产白丝娇喘喷水9色精品| 丝袜美足系列| 亚洲成色77777| 久久国内精品自在自线图片| 久久久久久久久久久久大奶| 日本欧美国产在线视频| 视频在线观看一区二区三区| 超碰97精品在线观看| 最新中文字幕久久久久| 丝袜喷水一区| 午夜久久久在线观看| 考比视频在线观看| 国产精品女同一区二区软件| 中国美白少妇内射xxxbb| 视频中文字幕在线观看| 国产一区二区三区综合在线观看 | 宅男免费午夜| 观看美女的网站| 国产免费现黄频在线看| 国内精品宾馆在线| a级毛片在线看网站| 热re99久久精品国产66热6| 99香蕉大伊视频| 亚洲av综合色区一区| 亚洲欧美成人综合另类久久久| 中文字幕亚洲精品专区| 国产国语露脸激情在线看| 少妇精品久久久久久久| 女性生殖器流出的白浆| 国产精品久久久久久av不卡| 国产av精品麻豆| 老司机亚洲免费影院| 久久久久久久久久久久大奶| 97精品久久久久久久久久精品| 18在线观看网站| 亚洲精品国产av蜜桃| 国产男人的电影天堂91| 午夜免费男女啪啪视频观看| 亚洲精品av麻豆狂野| 99视频精品全部免费 在线| 国产 精品1| 91aial.com中文字幕在线观看| 欧美 日韩 精品 国产| 草草在线视频免费看| 日本欧美视频一区| av天堂久久9| 国产白丝娇喘喷水9色精品| 最后的刺客免费高清国语| av网站免费在线观看视频| 建设人人有责人人尽责人人享有的| 国产欧美日韩综合在线一区二区| 欧美精品高潮呻吟av久久| 制服人妻中文乱码| 亚洲精华国产精华液的使用体验| 亚洲av电影在线观看一区二区三区| 欧美+日韩+精品| av卡一久久| 免费观看av网站的网址| 国产精品嫩草影院av在线观看| 国产免费一级a男人的天堂| 人妻系列 视频| 日本午夜av视频| 极品少妇高潮喷水抽搐| 亚洲精品国产色婷婷电影| 国产精品一区二区在线观看99| 欧美 亚洲 国产 日韩一| 国产极品天堂在线| 99久久人妻综合| 国产激情久久老熟女| 欧美激情国产日韩精品一区| 亚洲精品国产av成人精品| 国产女主播在线喷水免费视频网站| 少妇熟女欧美另类| 成人国产麻豆网| 久久精品熟女亚洲av麻豆精品| 久久久久久久精品精品| 99视频精品全部免费 在线| 少妇 在线观看| 久久精品国产自在天天线| 亚洲国产看品久久| 看免费av毛片| 国产av码专区亚洲av| 国产成人aa在线观看| 1024视频免费在线观看| 亚洲精品中文字幕在线视频| 美女国产视频在线观看| av视频免费观看在线观看| 国产成人免费无遮挡视频| 高清黄色对白视频在线免费看| 精品一区二区三区视频在线|