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

    熱應(yīng)激前后湖羊下丘腦差異表達新基因的篩選與注釋

    2019-09-10 07:22:44李隱俠馮小品張莉張俊錢勇孟春花王慧利仲躋峰曹少先
    江蘇農(nóng)業(yè)學(xué)報 2019年2期
    關(guān)鍵詞:熱應(yīng)激湖羊下丘腦

    李隱俠 馮小品 張莉 張俊 錢勇 孟春花 王慧利 仲躋峰 曹少先

    摘要:湖羊是中國南方唯一的綿羊品種,耐濕熱是其重要的特色性狀之一。但目前關(guān)于其耐濕熱機制尚不清楚。本研究以湖羊為研究對象,通過轉(zhuǎn)錄組測序技術(shù)鑒定熱應(yīng)激前后下丘腦組織中的新基因及其差異表達,并通過不同數(shù)據(jù)庫比對和分析,進行新基因功能注釋。結(jié)果表明,熱應(yīng)激前后在湖羊下丘腦組織中共檢測到1990個新轉(zhuǎn)錄本,通過功能注釋,發(fā)現(xiàn)1344個新基因。在熱應(yīng)激組和對照組湖羊中,差異表達新基因186個,其中130個新基因表達上調(diào),56個新基因表達下調(diào),通過比對發(fā)現(xiàn)這些基因分別與細胞免疫、繁殖、生長代謝及相關(guān)信號通路相關(guān)。Real-time PCR驗證試驗結(jié)果表明,10個熱應(yīng)激前后表達差異基因中有9個與轉(zhuǎn)錄組測序結(jié)果一致。這些結(jié)果為綿羊基因組的進一步完善及基因功能的挖掘提供了一定的數(shù)據(jù)基礎(chǔ),為全面解析綿羊熱應(yīng)激的分子機制提供更多的依據(jù)。

    關(guān)鍵詞:湖羊;熱應(yīng)激;下丘腦;轉(zhuǎn)錄組測序;real-time PCR

    中圖分類號:Q781

    文獻標(biāo)識碼:A

    文章編號:1000-4440(2019)02-0363-07

    熱應(yīng)激是指高溫環(huán)境下,動物機體對熱環(huán)境刺激產(chǎn)生的非特異性生理反應(yīng)的總稱。研究結(jié)果表明,影響家畜生產(chǎn)性能的主要環(huán)境因素有溫度、濕度、太陽輻射和風(fēng)速[2],在以高太陽輻射和高溫著稱的熱帶地區(qū),熱應(yīng)激是一個限制動物發(fā)育和生產(chǎn)的主要因素[3]。研究結(jié)果表明,急性熱應(yīng)激使雞生長性能下降,營養(yǎng)成分重分配,脂肪沉積增加[4-5];高溫環(huán)境破壞蛋雞內(nèi)分泌系統(tǒng)進而妨礙生殖活動[6];熱應(yīng)激影響綿羊的產(chǎn)奶量和繁殖性能[7-8]。

    熱應(yīng)激原因很簡單,但其影響和調(diào)控機制到目前為止尚不清楚。通過垂體將神經(jīng)系統(tǒng)和內(nèi)分泌系統(tǒng)連接起來的下丘腦組織是一個調(diào)節(jié)體溫的關(guān)鍵中樞[9]。眾所周知,下丘腦的前視區(qū)(POA)和前下丘腦區(qū)(AH)是體溫調(diào)節(jié)中心,下丘腦區(qū)域如室旁下丘腦核(PVN)、視神經(jīng)核(SO)和外側(cè)下丘腦在熱應(yīng)激狀態(tài)下被激活。熱應(yīng)激狀態(tài)下,下丘腦通過下丘腦-垂體-腎上腺軸、交感神經(jīng)系統(tǒng)和其他神經(jīng)內(nèi)分泌途徑保存或者消散熱量,從而調(diào)節(jié)動物體溫[10]。

    湖羊是中國南方主要的綿羊品種,具有非常重要的優(yōu)良特性,在中國低繁殖力綿羊品種遺傳改良和新品系培育中做出了一定的貢獻,其中耐濕熱是眾多優(yōu)良特性之一。本研究以湖羊下丘腦組織為研究對象,在熱應(yīng)激狀態(tài)和非熱應(yīng)激狀態(tài)下,通過轉(zhuǎn)錄組測序技術(shù)、參考基因組比對技術(shù)挖掘湖羊新轉(zhuǎn)錄本,解析調(diào)控湖羊熱應(yīng)激的差異表達新基因,并進行功能注釋,為進一步解析熱應(yīng)激調(diào)控的分子機制提供一定的理論依據(jù)。

    1 材料與方法

    1.1 試驗動物與試驗設(shè)計

    一周歲健康、體質(zhì)量[(33.00+1.58)kg]相近的湖羊母羊6只,急性熱應(yīng)激組和對照組各3只,自由采食和飲水。急性熱應(yīng)激組溫度42℃ 左右4 h,熱應(yīng)激模型和溫度參照文獻[11];對照組溫度22℃ 左右。基礎(chǔ)飼料組成及營養(yǎng)水平見文獻[12]。

    1.2 樣品采集及RNA提取

    熱應(yīng)激組和對照組湖羊處理后立即屠宰,取完整的下丘腦組織置于液氮中貯存。在液氮中將整個下丘腦組織研磨完全,用Trizol 法提取總RNA,并分別采用Nanodrop、Qubit 2.0和Agilent 2100方法檢測RNA樣品的純度、濃度和完整度等,符合測序要求的RNA送往北京百邁客生物科技有限公司進行轉(zhuǎn)錄組測序。

    1.3 cDNA文庫構(gòu)建和測序

    使用NEBNext UltraTM RNA文庫試劑盒(NEB,美國),按照其說明書進行文庫構(gòu)建,具體步驟參照文獻[13]。應(yīng)用llumina Hiseq 2500高通量測序平臺對構(gòu)建的eDNA文庫進行測序,測序讀取長度為125 bp(雙末端測序)。

    1.4 RNA-seq分析

    將原始數(shù)據(jù)進行數(shù)據(jù)過濾獲得高質(zhì)量的Cleandata,利用TopHat2將Clean reads 與參考基因組進行序列比對,并統(tǒng)計比對效率。使用Cufflinks軟件對新轉(zhuǎn)錄本和新基因的表達水平進行定量。同時對差異表達基因(DEG)進行統(tǒng)計,以差異倍數(shù)(FC)≥2且錯誤發(fā)現(xiàn)率(FDR)<0.01為標(biāo)準(zhǔn)。

    1.5 功能注釋

    使用BLAST軟件將差異表達新基因與Nr、Swiss-Prot、GO、COG、KOG、KEGG數(shù)據(jù)庫進行序列比對,使用KOBAS2.0得到新基因的KEGG orthology結(jié)果,預(yù)測完新基因編碼的蛋白質(zhì)氨基酸序列后使用HMMER軟件與Pfam數(shù)據(jù)庫比對,獲得新基因的注釋信息。

    1.6 Real-time PCR

    以轉(zhuǎn)錄組測序相對應(yīng)的RNA為模板,用逆轉(zhuǎn)錄試劑盒(TaKaRa,大連)逆轉(zhuǎn)錄為cDNA。Real-time PCR反應(yīng)體系為2.0 μleDNA模板,10 U的上下游引物各0.4 μl,2 ChamQ SYBR qPCR Master Mix10.0 μl。反應(yīng)程序為95℃ 預(yù)變性30s;95℃ 10s,60℃ 30s,40個循環(huán);95℃ 15s,60℃ 60s,95℃ 15s。每個樣本重復(fù)3次。

    2 結(jié)果與分析

    2.1 測序數(shù)據(jù)結(jié)果

    完成6只湖羊下丘腦組織的轉(zhuǎn)錄組分析,共獲得48.66 Gb Clean data,Q30堿基百分比在89.79%以上,得到大量高質(zhì)量原始數(shù)據(jù)。分別將各樣本的Cleanreads與綿羊參考基因組進行序列比對,比對效率從65.05%到71.20%不等。

    2.2 新基因預(yù)測及數(shù)據(jù)庫注釋

    通過與綿羊參考基因組(Oar_v4.0)進行序列比對,在2組湖羊下丘腦組織中共檢測到1990個新基因,其中綿羊以外數(shù)據(jù)庫被注釋的新基因為1344個。新基因在不同數(shù)據(jù)庫中進行注釋,結(jié)果發(fā)現(xiàn)Nr數(shù)據(jù)庫新基因數(shù)量最多,達1342個,而COG數(shù)據(jù)庫注釋的新基因數(shù)量最少,為78個(表1)。Nr數(shù)據(jù)庫中新基因同源物種分布結(jié)果顯示約99.98%新轉(zhuǎn)錄本注釋到最高的10個物種(圖1),即牛(Bostaurus)(31.22%)、野牦牛(Bos mutus)(21.16%)、綿羊(Ovisaries)(13.34%)、中國水牛(Bubalus bubalis)(5.51%)、山羊(Capra hircus)(5.44%)、藏羚羊(Pantholops hodgsonii)(4.84%)、野牛(Bison bison)(3.87%)、黑狐(Pteropus alecto)(1.12%)、海豚(Tursiops truncatus)(1.04%)和豬(Sus scrofa)(0.89%)。

    2.3 新基因GO分類

    GO分類發(fā)現(xiàn)690個新基因分別分類到50個基因本體(GO)條目中,細胞組分19個、分子功能13個和生物學(xué)過程18個(圖2),其中動物生長、繁殖、免疫系統(tǒng)、代謝及信號通路等生物學(xué)過程新基因所占比例最大,細胞組分次之,分子功能最少。

    2.4 熱應(yīng)激前后差異表達新基因的篩選和注釋

    對熱應(yīng)激和對照組湖羊的下丘腦組織差異表達新基因進行篩選。結(jié)果表明,兩組間共有186個新基因表達差異顯著,其中上調(diào)基因130個,下調(diào)基因56個。186個新基因只有146個在不同數(shù)據(jù)庫中被注釋,其中在Nr數(shù)據(jù)庫中被注釋的差異表達新基因數(shù)目最多(146個),在COG數(shù)據(jù)庫中最少(6個)(表2)。

    45個差異表達新基因參與到65個信號通路,包含氧化磷酸化、MAPK信號通路、鈣信號通路、代謝通路等。根據(jù)GO分類注釋信息發(fā)現(xiàn),其中參與湖羊生長調(diào)控的新基因有4個,參與繁殖調(diào)控的新基因有1個,參與免疫調(diào)控的新基因有2個,參與代謝調(diào)控的新基因有9個,參與氧化應(yīng)激調(diào)控的新基因有21個,對應(yīng)新基因編號及參與功能見表3。

    2.5 熱應(yīng)激前后差異表達新基因的Real-time PCR驗證

    為了進一步確認轉(zhuǎn)錄組測序得到的差異表達新基因在熱應(yīng)激前后表達量差異,設(shè)計12對引物(表4),以熱應(yīng)激前后湖羊下丘腦組織cDNA為模板,Real-time PCR驗證熱應(yīng)激前后基因表達差異。結(jié)果發(fā)現(xiàn)在10個表達有差異的新基因組中,9個基因的驗證結(jié)果都與轉(zhuǎn)錄組測序的結(jié)果一致,Ovis_aries_newGene_12971基因轉(zhuǎn)錄組測序結(jié)果顯示熱應(yīng)激后該基因表達量顯著下降,而Real-time PCR結(jié)果顯示該基因熱應(yīng)激后表達量下降,但是未達到顯著水平(圖3)。轉(zhuǎn)錄組測序熱應(yīng)激前后基因表達差異與Real-time PCR驗證熱應(yīng)激前后基因表達差異比較見表5。

    3 討論

    Marai等[13]根據(jù)溫濕指數(shù)(THI)將熱應(yīng)激劃分為不同的等級:22.2~23.3是普通熱應(yīng)激、23.4~25.6為嚴重?zé)釕?yīng)激>25.6屬于極嚴重?zé)釕?yīng)激。研究結(jié)果表明不同動物誘發(fā)熱應(yīng)激的溫度不同,對于羊而言,通常認為羊的臨界溫度為25~30℃ ,當(dāng)綿羊暴露于超過30℃時被認為發(fā)生熱應(yīng)激[14]。Srikandakumar等[15]在溫度為35.5~43.9℃ 時發(fā)現(xiàn)白毛澳大利亞美利奴產(chǎn)生熱應(yīng)激,呼吸頻率上升到每1min128次。DaCosta等在溫度為24~45℃ 的氣候室中觀察32只Polwarth母羊在不同級別熱應(yīng)激條件下飲食、反芻情況[11]。本研究根據(jù)Srikandaku-mar[15]和DaCosta[11]的研究模型選擇熱應(yīng)激時間和溫度對湖羊進行熱應(yīng)激研究。通過對對照組和熱應(yīng)激組湖羊下丘腦組織進行轉(zhuǎn)錄組測序,與之前的綿羊參考基因組序列相比,共檢測到1990個新轉(zhuǎn)錄本,其中有數(shù)據(jù)庫注釋的基因1344個。這些新轉(zhuǎn)錄本比對到GO、Nr、Swiss-Prot、KEGG數(shù)據(jù)庫的數(shù)量分別為689、1344、347和537,與這些數(shù)據(jù)庫的成功比對為進一步證實這些新基因的存在提供了理論依據(jù)。同源物種比對分析發(fā)現(xiàn),與牛的相似基因最多(31.22%),與綿羊的相似基因為13.34%,此結(jié)果一方面表明同為反芻動物的牛和羊在基因組序列中的高度保守性[1617],同時也表明現(xiàn)有的綿羊參考基因組的不完整性,本研究鑒定的綿羊新基因?qū)M一步完善綿羊參考基因組具有非常重要的意義。

    研究結(jié)果表明,當(dāng)動物經(jīng)歷熱應(yīng)激后,動物體溫上升,從而觸發(fā)進化保守的能調(diào)控多個細胞活性的熱應(yīng)激轉(zhuǎn)錄組反應(yīng)調(diào)節(jié)基因,包括蛋白折疊、蛋白講解轉(zhuǎn)運、代謝、DNA修復(fù)和復(fù)制等變化[1819],調(diào)控相關(guān)基因通路的變化,進而調(diào)控機體的氧化應(yīng)激[20]、生長[21]、繁殖[22]、代謝[23].免疫[24]等。本研究結(jié)果經(jīng)過與KEGG數(shù)據(jù)庫比對發(fā)現(xiàn),在熱應(yīng)激前后差異表達的新基因中有45個新基因參與到65個信號通路,包含氧化磷酸化、MAPK信號通路,鈣信號通路、代謝通路等。研究發(fā)現(xiàn)熱應(yīng)激能激活MAPK信號通路,進而阻斷細胞生長和引發(fā)凋亡[25],同時發(fā)現(xiàn)其參與熱應(yīng)激誘導(dǎo)的精子損傷[26],進而調(diào)控動物的繁殖。根據(jù)GO分類注釋信息,發(fā)現(xiàn)參與湖羊生長調(diào)控的新基因有4個,參與繁殖調(diào)控的新基因有1個,參與免疫調(diào)控的新基因有2個,參與代謝調(diào)控的新基因有9個,參與氧化應(yīng)激調(diào)控的新基因有21個。該結(jié)果表明,存在潛在的新基因,其在湖羊熱應(yīng)激后參與調(diào)控湖羊生長、繁殖、免疫、代謝和氧化應(yīng)激等功能。

    本試驗對湖羊熱應(yīng)激前后下丘腦組織差異表達新基因進行GO功能分類和信號通路分析,有助于進一步解析調(diào)控湖羊熱應(yīng)激的分子機制,為盡早挖掘出完備的湖羊耐濕熱特性的分子機理提供依據(jù)。

    參考文獻:

    [1]黃艷,周阿容,池文文,等.果蠅熱應(yīng)激生理機制研究進展[J].南方農(nóng)業(yè)學(xué)報,2018,49(5):912-921.

    [2]HULME P E.Adapting to climate change:is there scope for ecological management in the face of a global threat[J].Journal ofApplied Ecology,2005,42:784-294.

    [3]MCMANUS C,HERMUCHE P,PAIVA S R,et al.Geographical distribution of sheep breeds in brazil and their relationship with climatic and environmental factors as risk classification for conservation[J].Brazilian Journal Science and Technology,2014,1:1-15.

    [4]劉梅.急性熱應(yīng)激對肉仔雞生長性能及脂肪代謝的影響[J].動物營養(yǎng)學(xué)報,2011,23(5):862-868.

    [5]唐姣玉,謝宇潔,劉兆輝,等.大蒜素對熱應(yīng)激湘黃雞生產(chǎn)性能及營養(yǎng)物質(zhì)代謝的影響[J].江蘇農(nóng)業(yè)學(xué)報,2017,33(3):638-641.

    [6]TU W L,CHENGC Y,WANGS H,et al.Profiling of differential gene expression in the hypothalamus of broiler-type Taiwan country chickens in response to acute heat stress[J].Theriogenology,2016,85(3):483-494.

    [7]MARAI I F M,EI-DARAWANY A A,F(xiàn)ADIEL A,et al.Physiological traits as affected by heat stress in sheep-a review[J].Small Ruminant Research,2007,71:1-12.

    [8]FINOCCHIARO R,VAN KAAMJ BC H M,PORTOLANOB,et al.Effect of heat stress on production of Mediterranean dairy sheep[J].Journal of Dairy Science,2005,88:1855-1864.

    [9]HAMMEL H T,JACKSON D C,STOLWIJKJ A,et al.Temperature regulation by hypothalamic proportional control with an adjust-able set point[J].Journal of Applied Physiology,1963,18:1146-1154.

    [10]MILLER D B,O'CALLAGHAN J P.Neuroendocrine aspects of the response to stress[J].Metabolism,2002,51:5-10.

    [11]DA COSTA MJ,DA SOLVA RG,DE SOUZA R C. Effect of air temperature and humidity on ingestive behavior of sheep[J].International Journal of Biometeorology,1992,36(4):218-222.

    [12]李隱俠,王慧利,鄧唯,等.高溫環(huán)境下湖羊直腸溫度和血液生化指標(biāo)的監(jiān)測[J].中國畜牧獸醫(yī),2018,45(1):271-277.

    [13]MARAIIF M,AYYAT MS,ABD EI-MONEM U M.Growth performance and reproductive traits at first parity of New Zealand White female rabbits as affected by heat stress and its alleviation,under Egyptian conditions[J].Tropical Animal Health and Production,2001,33:451-462.

    [14]FUQUAY J W.Heat stress as it affects animal production[J].Journal of Animal Science,1981,52(1):164-174.

    [15]SRIKANDAKUMAR A,JOHNSON E H,MAHGOUB O.Efect of heat stress on respiratory rate,rectal temperature and blood chemistry in Omani and Australian Merino sheep[J].Small RuminantResearch,2003,49:193-198.

    [16]李隱俠,張俊,錢勇,等.烏骨綿羊NR5A2基因編碼區(qū)序列克隆及特征分析[J].江蘇農(nóng)業(yè)學(xué)報,2014,30(6):1375-1382.

    [17]張春香,張國林,郭麗娜,等.基于高通量轉(zhuǎn)錄組測序的山羊睪丸和附睪頭差異表達基因分析[J].畜牧獸醫(yī)學(xué)報,2014,45(3):391-401.

    [18]FEDER M E,HOFMANN G E.Heat-shock proteins,molecular chaperones,and the stress response:evolutionary and ecological physiology[J].Anual Review of Physiology,1999,61:243-282.

    [19]KREGEL K C.Heat shock proteins:modifying factors in physiological stress responses and acquired thermotolerance[J].Jourmal of Applied Physiology,2002,92(5):2177-2186.

    [20]ALEMUT W,PANDEY H O,SALILEW WONDIM D,et al.Oxidative and endoplasmic reticulum stress defense mechanisms of bovine granulosa cells exposed to heat stress[J].Theriogenology,2018,110:130-141.

    [21]PRAGNA P,SEJIAN V,BAGATH M,et al.Comparative assessment of growth performance of three different indigenous goatbreeds exposed to summer heat stress[J].Journal of Animal Physiology and Animal Nutrition (Berl),2018,102(4):825-836.

    [22]RAHMAN M B,SCHELLANDER K,LUCENO N L,et al.Heat stress responses in spermatozoa:Mechanisms and consequences forcattlefertility[J].Theriogenology,2018,113:102-112.

    [23]GANESAN S,SUMMERS C M,PEARCE S C,et al.Short-term heat stress altered metabolism and insulin signaling in skeletalmuscle[J].Journal of Animal Science,2018,96(1):154-167.

    [24]SAEED-ZIDANE M,LINDEN L,SALILEW-WONDIM D,et al.Cellular and exosome mediated molecular defense mechanism inbovine granulosa cells exposed to oxidative stress[J].PLoS ONE,2017,12(11):e0187569.

    [25]WADA T,PENNINGER J M.Mitogen-activated protein kinases in apoptosis regulation[J].Oncogene,2004,23(16):2838-2849.

    [26]RAHMAN M B,VANDAELE L,RIJSSELAERE T,et al.Bovine spermatozoa react to in vitro heat stress by activating the mitogen-activated protein kinase 14 signaling pathway[J].ReproductionFertility and Development,2014,26(2):245-257.

    猜你喜歡
    熱應(yīng)激湖羊下丘腦
    小而強大的下丘腦
    5G助力“智慧畜牧” 湖羊有了“健康碼”
    湖羊“地位”及“后時代”應(yīng)對策略
    湖羊的養(yǎng)殖探究——以金昌湖羊養(yǎng)殖為例
    科學(xué)家發(fā)現(xiàn)控制衰老開關(guān)
    中老年健康(2017年9期)2017-12-13 07:16:39
    中藥對下丘腦作用的研究進展
    中成藥(2017年6期)2017-06-13 07:30:34
    熱應(yīng)激對長毛兔生產(chǎn)性能及生理生化指標(biāo)的影響
    熱應(yīng)激對種豬的危害及防范措施
    湖羊巴氏桿菌病診治病例
    回藥阿夫忒蒙丸對失眠大鼠下丘腦5-HT和5-HIAA含量的影響
    久久香蕉精品热| 欧美一区二区亚洲| 99热6这里只有精品| 中文在线观看免费www的网站| 成人永久免费在线观看视频| 日本免费一区二区三区高清不卡| 高清毛片免费观看视频网站| 国产伦人伦偷精品视频| 一本久久中文字幕| 九色成人免费人妻av| a级毛片a级免费在线| 毛片一级片免费看久久久久 | 国产色爽女视频免费观看| 国产成人a区在线观看| 欧美高清成人免费视频www| 91狼人影院| 国产精品久久久久久久电影| 在线观看一区二区三区| 乱人视频在线观看| 午夜福利视频1000在线观看| 1024手机看黄色片| 99久国产av精品| 国产精品一区二区免费欧美| 亚洲va日本ⅴa欧美va伊人久久| 国产熟女欧美一区二区| 少妇人妻一区二区三区视频| 成年女人永久免费观看视频| 午夜福利在线观看吧| 成年免费大片在线观看| 99热精品在线国产| 欧美激情久久久久久爽电影| 日韩一本色道免费dvd| 男女做爰动态图高潮gif福利片| 国产男人的电影天堂91| 有码 亚洲区| 免费人成视频x8x8入口观看| 国产视频内射| 大又大粗又爽又黄少妇毛片口| 熟妇人妻久久中文字幕3abv| 国产高清视频在线播放一区| 亚洲av.av天堂| 露出奶头的视频| 夜夜爽天天搞| 国产视频内射| 赤兔流量卡办理| 网址你懂的国产日韩在线| 久久久久久九九精品二区国产| 18禁黄网站禁片免费观看直播| 天堂动漫精品| 狠狠狠狠99中文字幕| 亚洲最大成人手机在线| 亚洲人成网站在线播| 亚洲av中文字字幕乱码综合| 永久网站在线| 少妇被粗大猛烈的视频| 国产大屁股一区二区在线视频| 免费人成在线观看视频色| 一级黄色大片毛片| 成人国产一区最新在线观看| 美女免费视频网站| 99在线人妻在线中文字幕| 天美传媒精品一区二区| 在线看三级毛片| 淫妇啪啪啪对白视频| 国产毛片a区久久久久| 日本免费一区二区三区高清不卡| 他把我摸到了高潮在线观看| 国产伦人伦偷精品视频| 日本黄色片子视频| av在线老鸭窝| 国产成人a区在线观看| 亚洲最大成人中文| 国产女主播在线喷水免费视频网站 | 美女xxoo啪啪120秒动态图| 桃红色精品国产亚洲av| 久9热在线精品视频| 亚洲精品乱码久久久v下载方式| 国产亚洲精品综合一区在线观看| 免费看光身美女| 我的老师免费观看完整版| 在线观看一区二区三区| 国产91精品成人一区二区三区| 久久久精品欧美日韩精品| 日本一二三区视频观看| 特大巨黑吊av在线直播| 中文字幕免费在线视频6| 日本三级黄在线观看| 中出人妻视频一区二区| 在线播放国产精品三级| 黄色一级大片看看| 男女下面进入的视频免费午夜| 丰满乱子伦码专区| 桃红色精品国产亚洲av| 日韩 亚洲 欧美在线| xxxwww97欧美| 久久久久久久久久久丰满 | 黄色日韩在线| 精品久久久久久,| 亚洲熟妇中文字幕五十中出| 国产亚洲av嫩草精品影院| 欧美高清性xxxxhd video| av在线亚洲专区| 两个人的视频大全免费| 亚洲av免费在线观看| a级毛片免费高清观看在线播放| 亚洲五月天丁香| 久久久国产成人精品二区| 日日摸夜夜添夜夜添小说| 国语自产精品视频在线第100页| 国产一区二区三区av在线 | 在线a可以看的网站| 2021天堂中文幕一二区在线观| 在线观看一区二区三区| 中亚洲国语对白在线视频| 老司机午夜福利在线观看视频| 大又大粗又爽又黄少妇毛片口| 在线观看66精品国产| 亚洲欧美日韩无卡精品| 亚洲最大成人av| 给我免费播放毛片高清在线观看| 欧美另类亚洲清纯唯美| 91狼人影院| 男女啪啪激烈高潮av片| 精品久久国产蜜桃| 日韩 亚洲 欧美在线| 精品久久久久久久久久免费视频| 永久网站在线| 久久人人爽人人爽人人片va| 亚洲精品久久国产高清桃花| 色尼玛亚洲综合影院| 久久久久久久久久成人| 色在线成人网| 欧美一区二区亚洲| 久久精品国产99精品国产亚洲性色| 免费观看精品视频网站| 国产熟女欧美一区二区| 色吧在线观看| 日本欧美国产在线视频| 99热这里只有是精品50| 永久网站在线| 99在线人妻在线中文字幕| 日本一二三区视频观看| 亚洲天堂国产精品一区在线| 免费无遮挡裸体视频| 97人妻精品一区二区三区麻豆| 久久精品久久久久久噜噜老黄 | 亚洲欧美日韩东京热| 国产不卡一卡二| 一本久久中文字幕| 18+在线观看网站| 国产一级毛片七仙女欲春2| 国产精品,欧美在线| 中国美女看黄片| 亚洲国产精品久久男人天堂| 99视频精品全部免费 在线| 成人国产麻豆网| 一个人免费在线观看电影| 很黄的视频免费| av专区在线播放| 国产女主播在线喷水免费视频网站 | 深爱激情五月婷婷| 国产综合懂色| 精品人妻偷拍中文字幕| 欧美日韩乱码在线| 欧美极品一区二区三区四区| 国产精品98久久久久久宅男小说| 精品乱码久久久久久99久播| 精品欧美国产一区二区三| 最近最新中文字幕大全电影3| 12—13女人毛片做爰片一| 中文字幕熟女人妻在线| 成人亚洲精品av一区二区| 欧美三级亚洲精品| 欧美成人性av电影在线观看| 啦啦啦韩国在线观看视频| 人人妻人人澡欧美一区二区| 久久亚洲真实| 国产av麻豆久久久久久久| 亚洲成人精品中文字幕电影| 一个人免费在线观看电影| 国产精品久久久久久久久免| .国产精品久久| 黄色女人牲交| 久久午夜亚洲精品久久| 国产亚洲精品综合一区在线观看| 国产一区二区在线观看日韩| 听说在线观看完整版免费高清| 国产美女午夜福利| 欧美人与善性xxx| 欧美+亚洲+日韩+国产| 成人午夜高清在线视频| 成人欧美大片| 男女下面进入的视频免费午夜| 午夜福利在线观看吧| 午夜视频国产福利| ponron亚洲| 色哟哟·www| 国产中年淑女户外野战色| 久久精品国产鲁丝片午夜精品 | 91麻豆精品激情在线观看国产| 性色avwww在线观看| 九色国产91popny在线| 欧美三级亚洲精品| 国产一区二区在线av高清观看| 简卡轻食公司| 亚洲,欧美,日韩| 深夜a级毛片| 精品免费久久久久久久清纯| 久久精品国产鲁丝片午夜精品 | 露出奶头的视频| 最新在线观看一区二区三区| 国产探花极品一区二区| 亚洲国产欧美人成| 女人被狂操c到高潮| 少妇高潮的动态图| 色吧在线观看| 99精品久久久久人妻精品| 午夜精品一区二区三区免费看| 亚洲一区二区三区色噜噜| 久久人人爽人人爽人人片va| 午夜爱爱视频在线播放| 少妇的逼水好多| 中文字幕av成人在线电影| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线播| 国产高清不卡午夜福利| 搡女人真爽免费视频火全软件 | 五月玫瑰六月丁香| 国产精品美女特级片免费视频播放器| 18+在线观看网站| 色在线成人网| 两个人视频免费观看高清| 成熟少妇高潮喷水视频| 欧美精品啪啪一区二区三区| 国内揄拍国产精品人妻在线| 欧美成人性av电影在线观看| 日日夜夜操网爽| 最好的美女福利视频网| 久久6这里有精品| 校园春色视频在线观看| av.在线天堂| 九色成人免费人妻av| 热99re8久久精品国产| 欧美日本视频| 一a级毛片在线观看| 欧美成人性av电影在线观看| 真实男女啪啪啪动态图| 亚洲一区二区三区色噜噜| 久久久久久久精品吃奶| 高清在线国产一区| 少妇高潮的动态图| 久久久久久国产a免费观看| 两人在一起打扑克的视频| 伦精品一区二区三区| 亚洲精品成人久久久久久| 国模一区二区三区四区视频| 天堂av国产一区二区熟女人妻| av中文乱码字幕在线| 91久久精品国产一区二区成人| 国内久久婷婷六月综合欲色啪| 午夜福利在线观看免费完整高清在 | 亚洲欧美激情综合另类| 亚洲国产精品合色在线| 五月玫瑰六月丁香| 欧美3d第一页| 波野结衣二区三区在线| 亚洲av成人av| 中出人妻视频一区二区| 亚洲中文字幕一区二区三区有码在线看| 色综合站精品国产| 成熟少妇高潮喷水视频| 两个人的视频大全免费| 嫩草影院新地址| 如何舔出高潮| 国产一级毛片七仙女欲春2| 亚洲国产欧美人成| 精品午夜福利视频在线观看一区| 99国产极品粉嫩在线观看| 女人十人毛片免费观看3o分钟| 久久久久精品国产欧美久久久| 亚洲第一区二区三区不卡| 成年女人永久免费观看视频| 日本在线视频免费播放| 亚洲真实伦在线观看| 欧美在线一区亚洲| 成人综合一区亚洲| av国产免费在线观看| 春色校园在线视频观看| 国产午夜精品论理片| 黄色女人牲交| bbb黄色大片| 日本熟妇午夜| 美女黄网站色视频| 国产v大片淫在线免费观看| 日韩国内少妇激情av| 男女做爰动态图高潮gif福利片| 国产精品嫩草影院av在线观看 | 少妇丰满av| 日韩欧美精品免费久久| 小蜜桃在线观看免费完整版高清| 男女做爰动态图高潮gif福利片| 午夜精品在线福利| 在线天堂最新版资源| 久久久久久久精品吃奶| 黄色一级大片看看| 99国产精品一区二区蜜桃av| 亚洲va日本ⅴa欧美va伊人久久| 亚洲最大成人av| 成人欧美大片| 男女边吃奶边做爰视频| 久久久久性生活片| 亚洲精品亚洲一区二区| 亚洲成人久久爱视频| 在线天堂最新版资源| 久久久久久大精品| 啦啦啦观看免费观看视频高清| 国产精品一及| 老女人水多毛片| 男女下面进入的视频免费午夜| 女人被狂操c到高潮| 久久欧美精品欧美久久欧美| x7x7x7水蜜桃| 嫩草影视91久久| 亚洲国产色片| 国内少妇人妻偷人精品xxx网站| 嫁个100分男人电影在线观看| a级一级毛片免费在线观看| 成人毛片a级毛片在线播放| 日本爱情动作片www.在线观看 | 成人永久免费在线观看视频| 精品一区二区免费观看| 国产亚洲欧美98| av天堂中文字幕网| 成人特级av手机在线观看| 中文字幕免费在线视频6| 我的老师免费观看完整版| 观看美女的网站| 亚洲天堂国产精品一区在线| 人人妻,人人澡人人爽秒播| 99riav亚洲国产免费| 色综合婷婷激情| 久久亚洲精品不卡| 午夜激情欧美在线| 午夜福利在线观看免费完整高清在 | 国产色爽女视频免费观看| 国产三级在线视频| 国内精品宾馆在线| 国产乱人伦免费视频| 国产一区二区亚洲精品在线观看| 又爽又黄a免费视频| 自拍偷自拍亚洲精品老妇| 给我免费播放毛片高清在线观看| 又黄又爽又免费观看的视频| 精品久久国产蜜桃| 国产精品久久久久久久久免| 91久久精品电影网| 桃红色精品国产亚洲av| 欧美日韩瑟瑟在线播放| a级毛片a级免费在线| 干丝袜人妻中文字幕| 精品一区二区三区视频在线观看免费| 嫩草影院新地址| 日日摸夜夜添夜夜添小说| 精品一区二区免费观看| 免费看av在线观看网站| 我的老师免费观看完整版| 免费搜索国产男女视频| 国产亚洲精品av在线| 精品乱码久久久久久99久播| 内地一区二区视频在线| av在线观看视频网站免费| 日本一二三区视频观看| 国产在线男女| АⅤ资源中文在线天堂| 色尼玛亚洲综合影院| 午夜福利在线在线| 国产淫片久久久久久久久| 人妻制服诱惑在线中文字幕| 久久久久久久午夜电影| 国产大屁股一区二区在线视频| 18禁黄网站禁片免费观看直播| 日韩 亚洲 欧美在线| av天堂中文字幕网| 国内毛片毛片毛片毛片毛片| 日韩大尺度精品在线看网址| av国产免费在线观看| 熟女人妻精品中文字幕| 亚洲av二区三区四区| 淫秽高清视频在线观看| 久久久精品大字幕| www.色视频.com| 最新在线观看一区二区三区| 国产精品一区二区三区四区免费观看 | 2021天堂中文幕一二区在线观| 亚洲av免费高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看| a级毛片免费高清观看在线播放| 一本精品99久久精品77| av国产免费在线观看| 午夜爱爱视频在线播放| 日本黄色视频三级网站网址| 最近最新免费中文字幕在线| 国产精品99久久久久久久久| 噜噜噜噜噜久久久久久91| 精品日产1卡2卡| 欧美日韩瑟瑟在线播放| 亚洲av第一区精品v没综合| 午夜视频国产福利| 欧美精品国产亚洲| 直男gayav资源| 精品久久久久久,| 乱系列少妇在线播放| 麻豆成人av在线观看| 一个人免费在线观看电影| 日韩一区二区视频免费看| 国产伦精品一区二区三区视频9| 国产单亲对白刺激| 精品人妻1区二区| 桃红色精品国产亚洲av| 伦理电影大哥的女人| 色哟哟·www| 亚洲人成网站高清观看| 久久婷婷人人爽人人干人人爱| 日本黄色片子视频| 日本五十路高清| 亚洲成人久久爱视频| 日本黄大片高清| 性插视频无遮挡在线免费观看| 国产成人一区二区在线| 俄罗斯特黄特色一大片| 午夜免费男女啪啪视频观看 | 国产国拍精品亚洲av在线观看| 国产av一区在线观看免费| 亚洲欧美日韩无卡精品| 在线免费十八禁| 国产高清不卡午夜福利| 色精品久久人妻99蜜桃| 国产精品一区二区三区四区免费观看 | 88av欧美| aaaaa片日本免费| 亚洲一区二区三区色噜噜| 亚洲精品久久国产高清桃花| h日本视频在线播放| 嫩草影院新地址| 国产午夜精品论理片| or卡值多少钱| 黄色丝袜av网址大全| 一a级毛片在线观看| 老司机午夜福利在线观看视频| 99久久久亚洲精品蜜臀av| 一本久久中文字幕| 亚洲色图av天堂| 国产伦人伦偷精品视频| 欧美+日韩+精品| 淫秽高清视频在线观看| 真人一进一出gif抽搐免费| 老熟妇乱子伦视频在线观看| 亚洲av免费高清在线观看| 日韩精品青青久久久久久| 久久久久久久午夜电影| 亚洲va日本ⅴa欧美va伊人久久| 久久天躁狠狠躁夜夜2o2o| 国产午夜福利久久久久久| 午夜亚洲福利在线播放| 亚洲av中文字字幕乱码综合| 国产精品人妻久久久影院| 人妻久久中文字幕网| 亚洲美女黄片视频| 啦啦啦啦在线视频资源| 国产男靠女视频免费网站| 午夜激情福利司机影院| 精品福利观看| 少妇裸体淫交视频免费看高清| 色综合婷婷激情| 精品久久国产蜜桃| 免费观看精品视频网站| 日本精品一区二区三区蜜桃| 亚洲在线自拍视频| 天堂√8在线中文| 国产色婷婷99| 一a级毛片在线观看| 欧美丝袜亚洲另类 | 日韩欧美免费精品| 伦精品一区二区三区| or卡值多少钱| 天堂av国产一区二区熟女人妻| 国产一区二区在线av高清观看| 久久精品国产鲁丝片午夜精品 | 18禁黄网站禁片午夜丰满| 国产极品精品免费视频能看的| 一边摸一边抽搐一进一小说| 99精品在免费线老司机午夜| 午夜福利欧美成人| 精品久久久久久成人av| 国产一区二区在线观看日韩| 如何舔出高潮| 久久国产乱子免费精品| 美女高潮的动态| 欧美日本亚洲视频在线播放| 麻豆精品久久久久久蜜桃| 亚洲精品在线观看二区| aaaaa片日本免费| 2021天堂中文幕一二区在线观| 亚洲人成伊人成综合网2020| 波多野结衣高清无吗| 三级毛片av免费| 12—13女人毛片做爰片一| 国产精品嫩草影院av在线观看 | 成人美女网站在线观看视频| 亚洲精品粉嫩美女一区| 成人精品一区二区免费| 欧美黑人欧美精品刺激| 日日夜夜操网爽| 亚洲精品粉嫩美女一区| 欧美潮喷喷水| 亚洲av美国av| 人妻久久中文字幕网| 成年女人看的毛片在线观看| 中国美女看黄片| 91狼人影院| 男人舔女人下体高潮全视频| 久久精品国产99精品国产亚洲性色| 精华霜和精华液先用哪个| 如何舔出高潮| 最新中文字幕久久久久| 亚洲久久久久久中文字幕| 国产主播在线观看一区二区| 中文在线观看免费www的网站| 高清日韩中文字幕在线| 亚洲欧美日韩高清在线视频| АⅤ资源中文在线天堂| 天堂av国产一区二区熟女人妻| 国产伦精品一区二区三区视频9| 国产v大片淫在线免费观看| 久久精品国产自在天天线| 色尼玛亚洲综合影院| 免费电影在线观看免费观看| 亚洲三级黄色毛片| 久久精品综合一区二区三区| 韩国av一区二区三区四区| 亚洲av一区综合| 在线观看舔阴道视频| 黄色视频,在线免费观看| 看免费成人av毛片| 午夜福利在线观看吧| 99视频精品全部免费 在线| 国产亚洲av嫩草精品影院| 91久久精品电影网| 国产视频内射| 少妇丰满av| 亚洲国产精品久久男人天堂| 国产精品一及| 日日夜夜操网爽| 国产精品精品国产色婷婷| 国产精品99久久久久久久久| 99久久无色码亚洲精品果冻| 欧美激情国产日韩精品一区| 国产高清有码在线观看视频| 国产一区二区三区视频了| 99在线视频只有这里精品首页| 婷婷精品国产亚洲av| 中亚洲国语对白在线视频| 伦精品一区二区三区| 亚洲七黄色美女视频| 亚洲经典国产精华液单| 禁无遮挡网站| 国产乱人伦免费视频| 桃红色精品国产亚洲av| 国产高清视频在线播放一区| www.www免费av| 亚洲专区国产一区二区| 男人的好看免费观看在线视频| 久久精品国产亚洲av涩爱 | 精品日产1卡2卡| 精品国产三级普通话版| 亚洲图色成人| 88av欧美| 午夜日韩欧美国产| 91在线精品国自产拍蜜月| 精品不卡国产一区二区三区| 成人av在线播放网站| 日韩中文字幕欧美一区二区| 简卡轻食公司| 亚洲综合色惰| 国产精品野战在线观看| 国内毛片毛片毛片毛片毛片| 三级国产精品欧美在线观看| 国内精品久久久久久久电影| www日本黄色视频网| 赤兔流量卡办理| 狠狠狠狠99中文字幕| 三级男女做爰猛烈吃奶摸视频| 午夜免费激情av| 国产一区二区三区视频了| 国产成人aa在线观看| 亚洲乱码一区二区免费版| 亚洲国产精品久久男人天堂| 欧美精品国产亚洲| 国产亚洲精品综合一区在线观看| 日本a在线网址| 网址你懂的国产日韩在线| 亚洲最大成人手机在线| 亚洲精品亚洲一区二区| 中出人妻视频一区二区| 色噜噜av男人的天堂激情| 亚洲不卡免费看| 美女高潮喷水抽搐中文字幕| 亚洲最大成人中文| 最后的刺客免费高清国语| 人妻夜夜爽99麻豆av| 欧美丝袜亚洲另类 | 欧美性猛交╳xxx乱大交人| 免费看光身美女| 国产成人影院久久av| 欧美日韩精品成人综合77777|