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

    7個(gè)地方山羊品種遺傳多樣性及遺傳結(jié)構(gòu)分析

    2020-07-10 06:42:14張樂超劉月琴段春輝張英杰王泳郭云霞
    生物技術(shù)通報(bào) 2020年6期
    關(guān)鍵詞:遺傳變異絨山羊多態(tài)

    張樂超 劉月琴 段春輝 張英杰 王泳 郭云霞

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

    山羊是人類馴化較早的動(dòng)物,為人們提供肉、絨、毛、皮、奶等多種畜產(chǎn)品,在畜牧生產(chǎn)中占據(jù)重要地位。我國山羊品種資源豐富,有著悠久的馴化和養(yǎng)殖歷史,就河北省及周邊地區(qū)就存在多個(gè)優(yōu)良地方山羊品種,如武安山羊、承德無角山羊、太行山羊、燕山絨山羊、遼寧絨山羊、內(nèi)蒙古絨山羊和濟(jì)寧青山羊等。這些地方品種,無論是在促進(jìn)當(dāng)?shù)亟?jīng)濟(jì)發(fā)展還是在促進(jìn)我國出口創(chuàng)匯中,均占有重要地位。隨著國外羊品種的大量引進(jìn)及肉用綿羊的大量養(yǎng)殖,這些優(yōu)秀的地方山羊品種的養(yǎng)殖規(guī)模和存欄數(shù)量不斷減少,作為生物多樣性的重要組成部分,亟需加強(qiáng)保護(hù)和利用。

    有許多研究表明,分子標(biāo)記可用于基因多樣性分析。到目前為止,有各種分子標(biāo)記研究遺傳多樣性,例如限制性片段長(zhǎng)度多態(tài)性(Restriction fragment length polymorphi,RFLP)、 隨 機(jī) 擴(kuò) 增 多態(tài) 性 DNA(Random amplified polymorphic DNA,RAPD)、擴(kuò)增片段長(zhǎng)度多態(tài)性(Rmplified fragment length polymorphism,AFLP)、簡(jiǎn)單重復(fù)序列(Simple sequence repeat,SSR)和單核苷酸多態(tài)性(Single nucleotide polymorphisms,SNP) 等[1-5]。SSR 又 稱為微衛(wèi)星標(biāo)記,由于其共顯性遺傳、高程度的等位基因多樣化、相對(duì)多態(tài)性豐富和廣泛的基因組分布,可以有效地促進(jìn)遺傳關(guān)系的建立,是分析生物遺傳多樣性和群體遺傳結(jié)構(gòu)的常用方法[6]。

    為充分保護(hù)和利用地方山羊遺傳資源的多樣性,本研究選取15對(duì)微衛(wèi)星位點(diǎn)分析河北省及周邊地區(qū)7個(gè)地方山羊群體的遺傳多樣性和群體結(jié)構(gòu),旨在為地方山羊品種的種質(zhì)資源保護(hù)和利用提供理論依據(jù)。

    1 材料與方法

    1.1 材料

    采集30只燕山絨山羊(青龍利紅絨山羊技術(shù)服務(wù)中心)、23只濟(jì)寧青山羊(灤平虎什哈養(yǎng)殖場(chǎng))、30只承德無角山羊(寬城立東養(yǎng)殖有限公司)、30只太行山羊(阜平縣銀洞山肉羊養(yǎng)殖專業(yè)合作社)、36只武安山羊(武安市馬家莊太行飼養(yǎng)場(chǎng))、31只遼寧絨山羊(遼寧絨山羊科技示范場(chǎng))的血液樣品,以及33只內(nèi)蒙古絨山羊(阿拉善型)(內(nèi)蒙古億維白絨山羊有限責(zé)任公司)的耳組織,213份樣品均為隨機(jī)取樣,且均為沒有親緣關(guān)系的母羊。利用血液基因組DNA提取試劑盒(DP348)天根生化科技有限公司(北京)提取基因組DNA(圖1)。

    圖1 DNA瓊脂糖凝膠電泳圖

    1.2 方法

    1.2.1 微衛(wèi)星位點(diǎn)的引物設(shè)計(jì) 參照聯(lián)合國糧農(nóng)組織(FAO)及國際遺傳學(xué)會(huì)(ISAG)的推薦并查閱相關(guān)文獻(xiàn),選用15對(duì)熒光標(biāo)記微衛(wèi)星引物(表1),引物由通用生物系統(tǒng)有限公司合成(安徽)。

    1.2.2 PCR擴(kuò)增 PCR擴(kuò)增的反應(yīng)體系為:ddH2O 14.8 μL,dNTP 0.4 μL,Buffer 2 μL,F(xiàn) 0.3 μL(20μmol/L),R 0.3 μL(20 μmol/L),DNA 模 板 2 μL,Taq 0.2 μL,共 20 μL。PCR 反應(yīng)條件為:94℃預(yù)變性5 min;94℃變性30 s,54℃退火35 s,72℃延伸40 s,共35個(gè)循環(huán);最終72℃延伸3 min(圖2)。

    1.2.3 毛細(xì)管電泳 將甲酰胺與分子量?jī)?nèi)標(biāo)按100∶1的體積比混勻后,取15 μL加入上樣板中,再加入1 μL稀釋10倍的PCR產(chǎn)物,然后使用3730XL測(cè)序儀進(jìn)行毛細(xì)管電泳(圖3)。

    1.2.4 遺傳變異分析 利用Genemarker中的Fragment(Plant)片段分析軟件對(duì)測(cè)序儀得到的原始數(shù)據(jù)進(jìn)行分析,將各泳道內(nèi)分子量?jī)?nèi)標(biāo)的位置與各樣品峰值的位置做比較分析,得到片段大小。把15個(gè)微衛(wèi)星位點(diǎn)的等位基因片段大小進(jìn)行編號(hào),由大到小依次為1-21。利用Popgene32計(jì)算等位基因頻率、有效等位基因數(shù)(Effective numbers of alleles,Ne)、期望雜 合 度(Expected heterozygosity,He)、觀測(cè)雜合度(Observed heterozygosity,Ho)、F統(tǒng)計(jì)量(Fit和Fst)、基因流(Gene flow,Nm)、Shannon指數(shù)(I)等指標(biāo)。利用PIC-Calc軟件計(jì)算多態(tài)信息含量(Polymorphism information content,PIC)。利用Dispan軟件統(tǒng)計(jì)各群體Nei’s遺傳距離,并聚類分析。利用MEGA4.0軟件構(gòu)建群體UPGMA系統(tǒng)發(fā)育樹。

    表1 15個(gè)微衛(wèi)星位點(diǎn)的引物序列

    圖2 PCR產(chǎn)物瓊脂糖凝膠電泳圖

    圖3 毛細(xì)管電泳檢測(cè)部分結(jié)果

    2 結(jié)果

    2.1 7個(gè)地方山羊品種群體內(nèi)遺傳變異分析

    7個(gè)地方山羊品種15個(gè)微衛(wèi)星位點(diǎn)的等位基因頻率分析結(jié)果見表2。由表2可知,15個(gè)位點(diǎn)共檢測(cè)到了172個(gè)等位基因。其中,BM6526位點(diǎn)的等位基因最多,為21個(gè),片段大小為152-191 bp;INRA063位點(diǎn)的等位基因最少,為6個(gè),片段大小為169-179 bp。在所檢測(cè)到的等位基因中ILSTS005位點(diǎn)181 bp片段基因頻率最高為0.689 6,其次為BM1225位點(diǎn)221 bp片段,基因頻率為0.669;而OarFCB48標(biāo)記的163 bp片段、OarFCB20標(biāo)記的112 bp片段、CSRD247標(biāo)記的248 bp片段、ILSTS087標(biāo)記的164 bp片段、TGLA53標(biāo)記的158 bp片段、BM1225標(biāo)記的252 bp片段和BMS574標(biāo)記的122 bp片段的基因頻率最低,均為0.002 3。

    7個(gè)地方山羊品種15個(gè)微衛(wèi)星位點(diǎn)的多態(tài)性分析結(jié)果見表3。由表3可知,所有微衛(wèi)星位點(diǎn)的有效等位基因數(shù)(Ne)范圍為1.928 1(ILSTS005)-10.257(BM6526), 觀 察 雜 合 度(Ho) 范 圍 為 0.436 6(BM1225)- 0.887 3(BMS574),期望雜合度(He)范圍為 0.482 5(ILST005)-0.904 7(BM6526),平均多態(tài)信息含量(PIC)范圍為0.439 9(ILSTS005)-0.894 9(BM6526),均屬于高度多態(tài)。

    7個(gè)地方山羊品種15個(gè)微衛(wèi)星位點(diǎn)的遺傳變異分析結(jié)果見表4。由表4可知,7個(gè)地方山羊群體15個(gè)微衛(wèi)星位點(diǎn)Shannon指數(shù)(I)均大于1,有效等位基因數(shù)在3-5之間,等位基因數(shù)較多,有較高的觀察雜合度和較高的期望雜合度,平均多態(tài)信息含量范圍為0.637 1-0.748 5,為高度多態(tài)。

    2.2 微衛(wèi)星位點(diǎn)的遺傳分化

    各微衛(wèi)星位點(diǎn)F-統(tǒng)計(jì)量分析見表5。由表5可知,7個(gè)山羊群體的平均群體間基因分化系數(shù)(Fst)為0.544 2,平均基因流(Nm)為0.209 4,總?cè)后w近交系數(shù)(Fit)為0.088 3。

    表2 7個(gè)地方品種山羊15個(gè)微衛(wèi)星位點(diǎn)的等位基因頻率

    2.3 7個(gè)地方山羊品種遺傳距離及聚類分析

    由表6可知,7個(gè)地方山羊品種的相似系數(shù)范圍為0.764 8-0.927 4,遺傳距離范圍為0.075 4-0268 2。武安山羊和承德無角山羊相似系數(shù)最大(0.927 4),遺傳距離最近(0.075 4);太行山羊和燕山絨山羊相似系數(shù)最小(0.764 8),遺傳距離最遠(yuǎn)(0.268 2)。同時(shí),依據(jù)遺傳相似系數(shù)對(duì)7個(gè)地方山羊群體進(jìn)行聚類分析,結(jié)果見圖4,7個(gè)山羊群體被聚為2類,濟(jì)寧青山羊、承德無角山羊、武安山羊、內(nèi)蒙古絨山羊和太行山羊聚為一類,燕山絨山羊和遼寧絨山羊聚為一類。

    表3 7個(gè)地方山羊品種15個(gè)微衛(wèi)星位點(diǎn)多態(tài)信息指標(biāo)

    表4 7個(gè)山羊群體在15個(gè)微衛(wèi)星位點(diǎn)遺傳變異結(jié)果

    表5 15個(gè)微衛(wèi)星位點(diǎn)F-統(tǒng)計(jì)量分析表

    表6 7個(gè)地方山羊群體間的Nei’s遺傳距離和相似系數(shù)

    3 討論

    3.1 7個(gè)地方山羊品種群體內(nèi)遺傳變異分析

    等位基因頻率是衡量動(dòng)物群體遺傳結(jié)構(gòu)、估計(jì)遺傳變異的重要指標(biāo)[7],群體中頻率最高、片段最大的等位基因是該物種最保守的基因,而其余等位基因則是在進(jìn)化過程中由該基因突變產(chǎn)生的[8]。王可等[9]利用12個(gè)微衛(wèi)星標(biāo)記研究濟(jì)寧青山羊群體遺傳多樣性,發(fā)現(xiàn)SRCRSP8位點(diǎn)的230 bp和232 bp片段的基因頻率為0.281 3和0.322 9,遠(yuǎn)大于236 bp和240 bp的基因頻率0.010 4,認(rèn)為前兩個(gè)片段可能為原始基因,后兩個(gè)片段可能是突變基因。本試驗(yàn)也發(fā)現(xiàn)了基因頻率相差較大的片段,如ILSTS005位點(diǎn)181 bp片段基因頻率為0.689 6,173 bp片段基因頻率為0.002 4,推測(cè)181 bp片段可能為原始基因,173 bp片段可能為突變基因。有效等位基因數(shù)可反映群體遺傳變異,其數(shù)值越接近等位基因數(shù),等位基因在群體中分布越均勻[10]。本試驗(yàn)發(fā)現(xiàn)BM6526位點(diǎn)等位基因數(shù)最多(21個(gè))片段大小為152-191 bp,說明BM6526位點(diǎn)變異程度大,對(duì)研究山羊選育更有價(jià)值。武志娟等[11]利用14個(gè)微衛(wèi)星標(biāo)記分析5個(gè)黑山羊群體,發(fā)現(xiàn)其中有11個(gè)位點(diǎn)的等位基因數(shù)和有效等位基因數(shù)差異大,認(rèn)為這種不平衡可能是由于地理隔離和人工選擇導(dǎo)致的。Lin等[12]研究亞洲10個(gè)山羊群體的遺傳多樣性和遺傳結(jié)構(gòu)發(fā)現(xiàn),群體離馴養(yǎng)中心越遠(yuǎn),其遺傳多樣性漸漸下降。王可等[13]研究濟(jì)寧青山羊4個(gè)群體同樣發(fā)現(xiàn)等位基因數(shù)要高于Ne。Hines等[14]認(rèn)為等位基因分布不均可能與近交有關(guān)。本試驗(yàn)中15個(gè)微衛(wèi)星位點(diǎn)等位基因數(shù)與Ne相差較多,片段大小不同,頻率分布不均勻,說明7個(gè)地方山羊群體有豐富的遺傳多樣性。

    基因雜合度與品種的變異程度具有正相關(guān)關(guān)系,即基因雜合度越大,品種的變異程度越大[15]。若一個(gè)群體的平均雜合度大于0.5,表明該群體的遺傳多樣性豐富[16]。郭丹等[17]利用7個(gè)微衛(wèi)星位點(diǎn)分析遼寧絨山羊遺傳多樣性,發(fā)現(xiàn)群體的平均雜合度為0.638,說明遼寧絨山羊品種內(nèi)的遺傳變異處于較高水平,選擇潛力大。趙艷紅[18]分析阿拉善型內(nèi)蒙古絨山羊的遺傳多樣性,發(fā)現(xiàn)群體的平均雜合度為0.613 6,同樣說明內(nèi)蒙古絨山羊遺傳多樣性豐富。本研究中,所檢測(cè)的群體內(nèi)平均雜合度比較高,說明群體的遺傳變異較大,適應(yīng)環(huán)境的能力較強(qiáng)。并且,不同微衛(wèi)星在同一山羊群體,同一微衛(wèi)星在不同山羊群體中變異程度也有很大差別。BM6526位點(diǎn)的He最高為0.904 7,ILSTS005位點(diǎn)相對(duì)其它微衛(wèi)星位點(diǎn)較低為0.482 5;燕山絨山羊最低為0.647 9,內(nèi)蒙古絨山羊的He最高為0.760 2。7個(gè)群體中僅太行山羊Ho > He,其他6個(gè)群體均是Ho < He,說明除太行山羊外,其他6個(gè)群體均存在不同程度的近交,承德無角山羊近交程度高于其他群體。王可等[9]研究濟(jì)寧青山羊遺傳多樣性發(fā)現(xiàn),各標(biāo)記的觀測(cè)雜合度在0.428 6-0.904 8之間,期望雜合度在

    0.598 1-0.839 7之間,平均期望雜合度為0.757 4,屬于高度雜合,遺傳變異大。I是研究物種多樣性的指標(biāo)之一,與群體的遺傳變異程度有正相關(guān)關(guān)系,即I值越大,群體的遺傳變異程度越高[19]。本試驗(yàn)中,7個(gè)地方山羊群體I值均大于1,表明試驗(yàn)群體遺傳變異程度較高,受近交的影響較小。多態(tài)信息含量(PIC)反映了微衛(wèi)星標(biāo)記的變異程度和估計(jì)群體內(nèi)的遺傳變異程度,PIC < 0.25時(shí)為低度多態(tài),0.25 0.5時(shí)為高度多態(tài)[20],PIC > 0.7時(shí)認(rèn)為是最優(yōu)的遺傳標(biāo)記[21]。本試驗(yàn)的15個(gè)微衛(wèi)星位點(diǎn)僅有ILSTS005為中度多態(tài),其余均為高度多態(tài),有11個(gè)微衛(wèi)星位點(diǎn)大于0.7,說明試驗(yàn)地方山羊群體具有豐富的遺傳基礎(chǔ),有很大的保種潛力。

    圖4 7個(gè)地方山羊群體聚類系統(tǒng)發(fā)育樹

    3.2 7個(gè)地方山羊品種的遺傳距離及聚類分析

    Fst反映了群體間的遺傳分化程度,F(xiàn)st在0-0.05時(shí),各群體間遺傳分化很小,在0.05-0.15時(shí),為中等程度的遺傳分化,在0.15-0.25時(shí),為高度遺傳分化;Nm數(shù)值小于1時(shí),群體可能向遺傳分化的方向發(fā)展,Nm數(shù)值大于1時(shí),群體能有效地抑制由遺傳漂變引起的遺傳分化[22]。Radhika等[23]利用 10 個(gè)微衛(wèi)星位點(diǎn)分析印度本地山羊和雜交山羊品種的遺傳多樣性和群體結(jié)構(gòu)發(fā)現(xiàn),二者均具有高的遺傳多樣性,并且Fst值為0.02,說明群體間遺傳分化很小。本試驗(yàn)發(fā)現(xiàn)14個(gè)微衛(wèi)星位點(diǎn)在7個(gè)山羊群體的Fst值均在0.484以上,說明群體間的遺傳分化程度很大。Nm值均在0.162 7以上,說明7個(gè)群體間的基因交流較少。Fit反應(yīng)了總?cè)后w的近交程度,表示整個(gè)群體中攜帶的一對(duì)等位基因是同源的概率[24]。本試驗(yàn)Fit值為0.088 3,表明群體間存在8.83%的遺傳變異來自遠(yuǎn)交,占的比例相對(duì)較小。

    遺傳距離表示不同種群或品種間的基因差異程度,通常利用等位基因頻率計(jì)算[25],是研究物種遺傳多樣性的基礎(chǔ),可反映群體的系統(tǒng)進(jìn)化,群體分化時(shí)間越短,遺傳距離就越小[26]。趙艷紅等[27]分析中國6個(gè)山羊群體遺傳多樣性發(fā)現(xiàn)遼寧絨山羊和內(nèi)蒙古絨山羊(阿拉善型)遺傳距離為0.134,與本試驗(yàn)的結(jié)果(0.1589)相似。楊清芳[28]利用15個(gè)微衛(wèi)星研究5個(gè)山羊品種遺傳多樣性發(fā)現(xiàn),遼寧絨山羊與承德無角山羊的遺傳距離為0.2270,與本試驗(yàn)研究結(jié)果(0.2341)相似。劉恩民[29]利用SNP芯片分析中國地方山羊品種的遺傳結(jié)構(gòu)發(fā)現(xiàn)內(nèi)蒙古絨山羊和太行山羊聚在一大類,與本試驗(yàn)結(jié)果相一致。Alatiyat等[30]利用17個(gè)微衛(wèi)星位點(diǎn)UPGMA聚類分析3種沙特山羊發(fā)現(xiàn),Jabli、Bishi和Tohami聚為一大類。本研究結(jié)果顯示,燕山絨山羊與遼寧絨山羊聚為一類,說明這兩個(gè)群體遺傳距離較近,燕山絨山羊與太行山羊遺傳距離為0.268 2,相對(duì)于其他群體最遠(yuǎn)。

    4 結(jié)論

    燕山絨山羊、遼寧絨山羊、承德無角山羊、濟(jì)寧青山羊、太行山羊、武安山羊和內(nèi)蒙古絨山羊7個(gè)地方山羊群體遺傳多樣性豐富,群體間遺傳分化程度大,基因交流少,受近交程度影響小,具有較高的利用價(jià)值和潛力。

    猜你喜歡
    遺傳變異絨山羊多態(tài)
    新華指數(shù) 新增岢嵐絨山羊
    分層多態(tài)加權(quán)k/n系統(tǒng)的可用性建模與設(shè)計(jì)優(yōu)化
    參差多態(tài)而功不唐捐
    先導(dǎo)編輯技術(shù)可編輯近90%的人類遺傳變異
    基于改進(jìn)遺傳變異算子的海島算法
    電子制作(2019年24期)2019-02-23 13:22:18
    遼寧絨山羊產(chǎn)業(yè)現(xiàn)狀及市場(chǎng)分析
    秋冬季節(jié)絨山羊的管理注意事項(xiàng)
    火力楠子代遺傳變異分析及優(yōu)良家系選擇
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    人多巴胺D2基因啟動(dòng)子區(qū)—350A/G多態(tài)位點(diǎn)熒光素酶表達(dá)載體的構(gòu)建與鑒定及活性檢測(cè)
    女人久久www免费人成看片 | 91午夜精品亚洲一区二区三区| 精品久久久久久久久亚洲| 高清日韩中文字幕在线| 亚洲va在线va天堂va国产| 秋霞在线观看毛片| 亚洲在线观看片| 欧美区成人在线视频| 欧美bdsm另类| or卡值多少钱| 51国产日韩欧美| 内地一区二区视频在线| 老司机影院毛片| 看非洲黑人一级黄片| 免费观看的影片在线观看| 亚洲不卡免费看| 国产成人freesex在线| 亚洲真实伦在线观看| 一级毛片电影观看 | 波野结衣二区三区在线| 高清日韩中文字幕在线| 一本一本综合久久| 久久精品国产鲁丝片午夜精品| 岛国毛片在线播放| 99久久精品热视频| 春色校园在线视频观看| 国产色婷婷99| 亚洲av中文av极速乱| 床上黄色一级片| 男插女下体视频免费在线播放| 日本三级黄在线观看| 国产精品国产高清国产av| 国产 一区精品| 舔av片在线| 亚洲欧美一区二区三区国产| 欧美日本视频| 日本午夜av视频| 免费观看人在逋| 国产精品国产三级国产专区5o | 69av精品久久久久久| 国产单亲对白刺激| 女人久久www免费人成看片 | 欧美日本视频| 黄片wwwwww| 一夜夜www| 干丝袜人妻中文字幕| 久久久午夜欧美精品| 男女下面进入的视频免费午夜| 嫩草影院新地址| 美女高潮的动态| 97超碰精品成人国产| 男人舔奶头视频| 一级毛片电影观看 | 日本黄色视频三级网站网址| 天堂√8在线中文| 在现免费观看毛片| 青春草视频在线免费观看| 亚洲欧美精品自产自拍| 一边亲一边摸免费视频| 成人av在线播放网站| 欧美97在线视频| 国产成人精品久久久久久| 国产精品国产三级国产av玫瑰| 一级毛片我不卡| 午夜久久久久精精品| 欧美日韩综合久久久久久| 只有这里有精品99| 校园人妻丝袜中文字幕| 中文字幕人妻熟人妻熟丝袜美| av福利片在线观看| 国产av不卡久久| 免费观看性生交大片5| 亚洲欧美精品综合久久99| 天天一区二区日本电影三级| av国产免费在线观看| 日本爱情动作片www.在线观看| 久久久国产成人精品二区| 国产精品无大码| 国产精品爽爽va在线观看网站| 欧美xxxx性猛交bbbb| a级一级毛片免费在线观看| 国产精品久久久久久久电影| 七月丁香在线播放| 我要看日韩黄色一级片| 久久99热这里只频精品6学生 | 免费av观看视频| 汤姆久久久久久久影院中文字幕 | 国产成人a∨麻豆精品| 亚洲精品乱码久久久v下载方式| 白带黄色成豆腐渣| 国产精品人妻久久久影院| 久久久久网色| 国产激情偷乱视频一区二区| 22中文网久久字幕| 欧美性感艳星| 久久草成人影院| 国产毛片a区久久久久| 成人国产麻豆网| 国语对白做爰xxxⅹ性视频网站| 观看免费一级毛片| 麻豆久久精品国产亚洲av| 高清av免费在线| 国产淫片久久久久久久久| 国产激情偷乱视频一区二区| 国产不卡一卡二| 日日干狠狠操夜夜爽| 99热这里只有精品一区| 寂寞人妻少妇视频99o| 国产黄a三级三级三级人| 亚洲欧美中文字幕日韩二区| 国产伦理片在线播放av一区| 国产成人精品久久久久久| 久久久久精品久久久久真实原创| 深爱激情五月婷婷| 99在线视频只有这里精品首页| av专区在线播放| 神马国产精品三级电影在线观看| 欧美一区二区亚洲| 色尼玛亚洲综合影院| av.在线天堂| 一边亲一边摸免费视频| 国产精品久久久久久久久免| 你懂的网址亚洲精品在线观看 | 特级一级黄色大片| 99久久中文字幕三级久久日本| 亚洲精品国产av成人精品| 夜夜爽夜夜爽视频| 日韩精品青青久久久久久| 18禁在线播放成人免费| 亚洲精品一区蜜桃| 亚洲国产色片| 97超碰精品成人国产| 97人妻精品一区二区三区麻豆| 91久久精品国产一区二区成人| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 久久久久久伊人网av| 亚洲欧美日韩高清专用| 久久精品国产亚洲av涩爱| 少妇丰满av| 插阴视频在线观看视频| 亚洲国产日韩欧美精品在线观看| 国产伦精品一区二区三区视频9| 亚洲国产最新在线播放| 中文字幕av在线有码专区| 自拍偷自拍亚洲精品老妇| 一区二区三区四区激情视频| 国产色婷婷99| 久久国产乱子免费精品| 老司机影院成人| 亚洲精华国产精华液的使用体验| 哪个播放器可以免费观看大片| 六月丁香七月| 午夜免费男女啪啪视频观看| 看非洲黑人一级黄片| 男人的好看免费观看在线视频| 可以在线观看毛片的网站| 免费看a级黄色片| 99久久无色码亚洲精品果冻| 精品久久久久久久人妻蜜臀av| 亚洲自拍偷在线| 一区二区三区免费毛片| 97超碰精品成人国产| 色噜噜av男人的天堂激情| 久久久久九九精品影院| 欧美+日韩+精品| 美女国产视频在线观看| 桃色一区二区三区在线观看| 欧美日韩一区二区视频在线观看视频在线 | 一区二区三区免费毛片| 欧美高清成人免费视频www| 日韩中字成人| 熟女人妻精品中文字幕| 国产精品1区2区在线观看.| 天堂√8在线中文| 久久久国产成人精品二区| 黄色一级大片看看| 2021天堂中文幕一二区在线观| 成年女人永久免费观看视频| 午夜爱爱视频在线播放| 久久久久免费精品人妻一区二区| 国产一区亚洲一区在线观看| 久久精品久久久久久久性| 亚洲成色77777| 国产黄a三级三级三级人| 久久精品国产99精品国产亚洲性色| 国产精品一区二区在线观看99 | 观看美女的网站| 尤物成人国产欧美一区二区三区| 自拍偷自拍亚洲精品老妇| 欧美成人一区二区免费高清观看| 亚洲,欧美,日韩| 国产成人免费观看mmmm| 免费看av在线观看网站| 成人毛片a级毛片在线播放| 我要搜黄色片| 国产精品野战在线观看| 99热网站在线观看| 国产一级毛片在线| 久久久亚洲精品成人影院| 国产高清有码在线观看视频| 看片在线看免费视频| 黄片无遮挡物在线观看| 狠狠狠狠99中文字幕| 亚洲人成网站高清观看| 日韩欧美精品v在线| 一区二区三区免费毛片| 日韩国内少妇激情av| 丰满人妻一区二区三区视频av| 中文字幕精品亚洲无线码一区| 亚洲,欧美,日韩| 日韩三级伦理在线观看| 久久久久久久国产电影| 午夜福利在线在线| 精品欧美国产一区二区三| 久久久久久久久久成人| a级毛片免费高清观看在线播放| 成年av动漫网址| ponron亚洲| 丝袜美腿在线中文| 日韩成人伦理影院| 又粗又硬又长又爽又黄的视频| 成人二区视频| 在线观看av片永久免费下载| 18禁裸乳无遮挡免费网站照片| 日韩在线高清观看一区二区三区| АⅤ资源中文在线天堂| 中文字幕免费在线视频6| 纵有疾风起免费观看全集完整版 | 九九久久精品国产亚洲av麻豆| 欧美激情在线99| 99久久成人亚洲精品观看| 你懂的网址亚洲精品在线观看 | 欧美97在线视频| 国产成人精品久久久久久| 99在线人妻在线中文字幕| 人妻少妇偷人精品九色| 天天一区二区日本电影三级| 婷婷色麻豆天堂久久 | 久久精品久久精品一区二区三区| 如何舔出高潮| 熟女电影av网| 国产精品,欧美在线| 国产美女午夜福利| 欧美高清成人免费视频www| 国产午夜精品久久久久久一区二区三区| av又黄又爽大尺度在线免费看 | 亚洲精品一区蜜桃| 亚洲中文字幕日韩| 色哟哟·www| 51国产日韩欧美| 亚洲性久久影院| 亚洲真实伦在线观看| 欧美一区二区亚洲| 中文字幕免费在线视频6| 夫妻性生交免费视频一级片| 高清日韩中文字幕在线| 久久久色成人| 大又大粗又爽又黄少妇毛片口| 久久热精品热| 嫩草影院精品99| 亚洲一区高清亚洲精品| 中国美白少妇内射xxxbb| 国产精品1区2区在线观看.| 看非洲黑人一级黄片| 毛片一级片免费看久久久久| 欧美3d第一页| 国产色爽女视频免费观看| 国产成人精品久久久久久| 99热网站在线观看| 一区二区三区四区激情视频| 综合色av麻豆| 中文字幕亚洲精品专区| 国产亚洲精品av在线| 日本一本二区三区精品| 国产成人a∨麻豆精品| 欧美一级a爱片免费观看看| 亚洲激情五月婷婷啪啪| 色哟哟·www| 国产一级毛片在线| 国产免费一级a男人的天堂| 亚洲av一区综合| 一级爰片在线观看| 国产一区二区三区av在线| 亚洲精品影视一区二区三区av| 亚洲av成人精品一区久久| 免费看a级黄色片| 国产一区二区亚洲精品在线观看| 欧美色视频一区免费| 中文字幕亚洲精品专区| 国产精品一区二区性色av| 久久99热这里只频精品6学生 | 国产精品乱码一区二三区的特点| 国产精品人妻久久久影院| 免费看av在线观看网站| 国产午夜福利久久久久久| 超碰av人人做人人爽久久| 最近最新中文字幕大全电影3| 91午夜精品亚洲一区二区三区| 久久6这里有精品| 精品国内亚洲2022精品成人| 久久精品国产自在天天线| 精品一区二区三区视频在线| 少妇的逼水好多| 成人亚洲欧美一区二区av| 99久久中文字幕三级久久日本| 一区二区三区免费毛片| 国产淫语在线视频| 男的添女的下面高潮视频| or卡值多少钱| 日韩欧美国产在线观看| 国产一区二区三区av在线| 午夜福利在线观看吧| 久久久a久久爽久久v久久| 国产午夜福利久久久久久| av在线天堂中文字幕| 国产精品永久免费网站| 亚洲真实伦在线观看| 男人和女人高潮做爰伦理| 一个人看视频在线观看www免费| 少妇丰满av| 少妇人妻精品综合一区二区| 在线观看66精品国产| 国产高清视频在线观看网站| 超碰97精品在线观看| 最近中文字幕2019免费版| 中国国产av一级| 18禁动态无遮挡网站| 边亲边吃奶的免费视频| 欧美成人a在线观看| 可以在线观看毛片的网站| 久久精品久久久久久噜噜老黄 | 亚洲av中文av极速乱| 免费观看的影片在线观看| 日韩国内少妇激情av| 亚洲伊人久久精品综合 | 真实男女啪啪啪动态图| 身体一侧抽搐| 国产高清有码在线观看视频| 欧美变态另类bdsm刘玥| 亚洲美女视频黄频| 国产精品久久久久久精品电影小说 | 性插视频无遮挡在线免费观看| 国产高清三级在线| 精品国产露脸久久av麻豆 | 日韩制服骚丝袜av| 亚洲婷婷狠狠爱综合网| av黄色大香蕉| 边亲边吃奶的免费视频| 一边摸一边抽搐一进一小说| 伦精品一区二区三区| 亚洲国产欧美人成| 中文在线观看免费www的网站| 国产亚洲5aaaaa淫片| 免费观看的影片在线观看| 色综合色国产| 精品一区二区免费观看| 久热久热在线精品观看| 汤姆久久久久久久影院中文字幕 | 色综合站精品国产| 别揉我奶头 嗯啊视频| 人妻夜夜爽99麻豆av| 男女视频在线观看网站免费| 日韩一区二区视频免费看| 日韩精品有码人妻一区| 99热全是精品| 纵有疾风起免费观看全集完整版 | 精品熟女少妇av免费看| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 日本色播在线视频| 最近的中文字幕免费完整| 美女脱内裤让男人舔精品视频| 我要搜黄色片| 国产色婷婷99| 精品一区二区免费观看| 18禁动态无遮挡网站| 日韩av不卡免费在线播放| 2021少妇久久久久久久久久久| av在线亚洲专区| 久久久久精品久久久久真实原创| 国产黄色视频一区二区在线观看 | 精品久久久久久久末码| 中国美白少妇内射xxxbb| 国产一区亚洲一区在线观看| 欧美成人午夜免费资源| 黄片无遮挡物在线观看| 亚洲欧美中文字幕日韩二区| 亚洲色图av天堂| 免费观看性生交大片5| 97超碰精品成人国产| 国产男人的电影天堂91| 中文字幕免费在线视频6| 日本与韩国留学比较| 午夜福利在线观看吧| 午夜精品一区二区三区免费看| 国内揄拍国产精品人妻在线| 久久人人爽人人片av| 身体一侧抽搐| 亚洲va在线va天堂va国产| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av.av天堂| 亚洲国产精品成人久久小说| 老女人水多毛片| 亚洲成人中文字幕在线播放| 三级毛片av免费| 亚洲中文字幕一区二区三区有码在线看| 久久精品久久精品一区二区三区| 亚洲人成网站在线观看播放| 国产成人免费观看mmmm| 欧美另类亚洲清纯唯美| 日本五十路高清| 午夜福利视频1000在线观看| 国产精品无大码| 国产精品久久久久久av不卡| 国产成人一区二区在线| 亚洲av男天堂| 中文字幕亚洲精品专区| 日韩精品青青久久久久久| 午夜精品在线福利| 在线a可以看的网站| 国产亚洲最大av| 成人欧美大片| 国产欧美另类精品又又久久亚洲欧美| 国产乱人偷精品视频| 国内揄拍国产精品人妻在线| 又爽又黄无遮挡网站| 日韩一本色道免费dvd| 18禁裸乳无遮挡免费网站照片| 欧美+日韩+精品| АⅤ资源中文在线天堂| 亚洲国产成人一精品久久久| 美女大奶头视频| 99久国产av精品国产电影| .国产精品久久| 在线观看美女被高潮喷水网站| 国产免费又黄又爽又色| 亚洲欧洲国产日韩| 久久久久性生活片| 国产老妇女一区| 一个人免费在线观看电影| 免费人成在线观看视频色| 日本wwww免费看| 亚洲精品影视一区二区三区av| 国产亚洲精品av在线| 欧美精品国产亚洲| 99热全是精品| 国产免费一级a男人的天堂| 亚洲婷婷狠狠爱综合网| 色5月婷婷丁香| 国产成人福利小说| 免费少妇av软件| 又大又黄又爽视频免费| 亚洲av电影在线进入| 97人妻天天添夜夜摸| 久久人妻熟女aⅴ| 美国免费a级毛片| 成人免费观看视频高清| 在线观看国产h片| 午夜福利乱码中文字幕| 爱豆传媒免费全集在线观看| 国产日韩欧美亚洲二区| av免费在线看不卡| 丝袜喷水一区| 亚洲国产欧美在线一区| 天天影视国产精品| 国产欧美另类精品又又久久亚洲欧美| 国产精品国产三级国产专区5o| 精品福利永久在线观看| 在现免费观看毛片| 午夜福利在线观看免费完整高清在| 亚洲高清免费不卡视频| 亚洲成色77777| 美女福利国产在线| 国产精品.久久久| 国产成人a∨麻豆精品| 亚洲av免费高清在线观看| 亚洲av电影在线进入| 狂野欧美激情性xxxx在线观看| 日韩成人伦理影院| 美女xxoo啪啪120秒动态图| 久久精品国产亚洲av涩爱| 亚洲 欧美一区二区三区| 免费看不卡的av| 久久午夜综合久久蜜桃| 尾随美女入室| 免费久久久久久久精品成人欧美视频 | 最近2019中文字幕mv第一页| 满18在线观看网站| a级毛片黄视频| 天天躁夜夜躁狠狠久久av| 51国产日韩欧美| 2022亚洲国产成人精品| 欧美成人午夜精品| 只有这里有精品99| 婷婷色av中文字幕| 中文字幕人妻熟女乱码| 欧美成人午夜免费资源| 人人妻人人澡人人看| 精品99又大又爽又粗少妇毛片| 亚洲熟女精品中文字幕| 99热这里只有是精品在线观看| 2022亚洲国产成人精品| 国产精品99久久99久久久不卡 | 久久99热6这里只有精品| 一级,二级,三级黄色视频| 男人爽女人下面视频在线观看| 亚洲人成77777在线视频| 春色校园在线视频观看| 伊人亚洲综合成人网| 一二三四在线观看免费中文在 | 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 欧美精品av麻豆av| 久久av网站| 一级片免费观看大全| 在线 av 中文字幕| 天天操日日干夜夜撸| 国语对白做爰xxxⅹ性视频网站| 99国产精品免费福利视频| 丰满迷人的少妇在线观看| 成年人免费黄色播放视频| 亚洲欧美日韩卡通动漫| 亚洲av电影在线进入| 国语对白做爰xxxⅹ性视频网站| 免费观看av网站的网址| 亚洲国产精品成人久久小说| 午夜福利网站1000一区二区三区| 国产精品99久久99久久久不卡 | 亚洲av在线观看美女高潮| 国产成人91sexporn| 国产xxxxx性猛交| 91久久精品国产一区二区三区| 熟女av电影| 啦啦啦视频在线资源免费观看| 99精国产麻豆久久婷婷| 曰老女人黄片| 国产黄色免费在线视频| 国产日韩欧美视频二区| 免费观看a级毛片全部| 纯流量卡能插随身wifi吗| 99热网站在线观看| 欧美97在线视频| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩亚洲高清精品| 成人国产麻豆网| 飞空精品影院首页| 亚洲av中文av极速乱| 下体分泌物呈黄色| 久久99精品国语久久久| 如何舔出高潮| 色婷婷av一区二区三区视频| 国产高清三级在线| 国国产精品蜜臀av免费| 久久久欧美国产精品| 中国国产av一级| 欧美日本中文国产一区发布| 亚洲高清免费不卡视频| 日本-黄色视频高清免费观看| 老司机亚洲免费影院| 男女啪啪激烈高潮av片| 热re99久久精品国产66热6| 久久久久精品久久久久真实原创| 亚洲国产精品国产精品| 亚洲精品乱久久久久久| 日韩视频在线欧美| av天堂久久9| 欧美激情国产日韩精品一区| 91在线精品国自产拍蜜月| 99国产综合亚洲精品| 成人无遮挡网站| 亚洲国产精品一区二区三区在线| 久久精品夜色国产| 制服诱惑二区| 黑人巨大精品欧美一区二区蜜桃 | videos熟女内射| 一级片'在线观看视频| 久久毛片免费看一区二区三区| 在线观看国产h片| 国产成人精品福利久久| 美国免费a级毛片| 久久久亚洲精品成人影院| 亚洲精品中文字幕在线视频| 99热国产这里只有精品6| 在线观看人妻少妇| 男女边吃奶边做爰视频| 在线精品无人区一区二区三| 亚洲高清免费不卡视频| 日本-黄色视频高清免费观看| 亚洲五月色婷婷综合| 丝袜人妻中文字幕| 飞空精品影院首页| 久久这里只有精品19| 免费人成在线观看视频色| 久久人人爽av亚洲精品天堂| 国产白丝娇喘喷水9色精品| 春色校园在线视频观看| 午夜日本视频在线| 中文字幕亚洲精品专区| 最近最新中文字幕免费大全7| 亚洲国产看品久久| 在线观看国产h片| 人人妻人人澡人人爽人人夜夜| 精品国产国语对白av| 男女下面插进去视频免费观看 | 黄片无遮挡物在线观看| 日韩制服丝袜自拍偷拍| 黑人高潮一二区| 久久久欧美国产精品| 黄色毛片三级朝国网站| 在线看a的网站| 国产亚洲一区二区精品| 亚洲国产av新网站| 成年动漫av网址| 日韩成人av中文字幕在线观看|