• <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含量的影響
    国产97色在线日韩免费| 黄色怎么调成土黄色| 热re99久久国产66热| 国产精品免费视频内射| 少妇精品久久久久久久| 国产熟女欧美一区二区| bbb黄色大片| 精品第一国产精品| 精品少妇一区二区三区视频日本电影 | 久久 成人 亚洲| 少妇 在线观看| 搡老岳熟女国产| 女性生殖器流出的白浆| 男人操女人黄网站| 欧美久久黑人一区二区| 丰满饥渴人妻一区二区三| 久久久精品免费免费高清| 免费女性裸体啪啪无遮挡网站| 久久这里只有精品19| 精品少妇久久久久久888优播| 秋霞在线观看毛片| 国产一卡二卡三卡精品 | 天天躁日日躁夜夜躁夜夜| 天堂俺去俺来也www色官网| 两个人免费观看高清视频| 亚洲欧美精品自产自拍| 亚洲成国产人片在线观看| 国产免费一区二区三区四区乱码| 97在线人人人人妻| 啦啦啦 在线观看视频| 国产精品无大码| 欧美精品av麻豆av| 亚洲欧美精品综合一区二区三区| 男女无遮挡免费网站观看| 男女之事视频高清在线观看 | 汤姆久久久久久久影院中文字幕| 国产精品久久久久久久久免| 国产精品99久久99久久久不卡 | 18禁国产床啪视频网站| 一本一本久久a久久精品综合妖精| 精品人妻一区二区三区麻豆| 波多野结衣av一区二区av| 久久精品熟女亚洲av麻豆精品| 亚洲人成网站在线观看播放| 在线天堂中文资源库| 制服丝袜香蕉在线| 亚洲精品久久久久久婷婷小说| 欧美日韩一级在线毛片| 色吧在线观看| 国产日韩一区二区三区精品不卡| 毛片一级片免费看久久久久| 曰老女人黄片| 日韩制服丝袜自拍偷拍| 国产精品二区激情视频| 哪个播放器可以免费观看大片| 91老司机精品| 国产爽快片一区二区三区| 精品少妇黑人巨大在线播放| 最黄视频免费看| 精品一区在线观看国产| 人人妻人人添人人爽欧美一区卜| 老司机深夜福利视频在线观看 | 免费黄网站久久成人精品| 亚洲少妇的诱惑av| 在线观看人妻少妇| 成年动漫av网址| 女性生殖器流出的白浆| 久久久亚洲精品成人影院| 最黄视频免费看| 欧美日韩亚洲国产一区二区在线观看 | 欧美成人精品欧美一级黄| 国产日韩欧美在线精品| 在线天堂中文资源库| 大话2 男鬼变身卡| 永久免费av网站大全| 日韩精品有码人妻一区| 亚洲国产最新在线播放| 天天影视国产精品| 亚洲精品自拍成人| 久久青草综合色| 久久久久久人人人人人| 欧美黑人欧美精品刺激| 国产精品麻豆人妻色哟哟久久| xxx大片免费视频| 日韩不卡一区二区三区视频在线| 精品久久蜜臀av无| 天堂8中文在线网| 老司机深夜福利视频在线观看 | 国产麻豆69| 亚洲av电影在线观看一区二区三区| 精品卡一卡二卡四卡免费| 亚洲精品第二区| tube8黄色片| 波野结衣二区三区在线| 欧美日韩亚洲国产一区二区在线观看 | av卡一久久| 亚洲伊人色综图| 精品国产露脸久久av麻豆| 一级毛片 在线播放| 我要看黄色一级片免费的| 国产一区有黄有色的免费视频| 亚洲七黄色美女视频| 免费久久久久久久精品成人欧美视频| 一区二区三区激情视频| 国产一区二区激情短视频 | 成人黄色视频免费在线看| 亚洲色图综合在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 男女边摸边吃奶| 国产日韩一区二区三区精品不卡| 亚洲美女搞黄在线观看| 免费观看av网站的网址| 亚洲av电影在线观看一区二区三区| 妹子高潮喷水视频| 国产男女超爽视频在线观看| 久久久国产一区二区| 日韩 亚洲 欧美在线| 黄片小视频在线播放| 亚洲色图 男人天堂 中文字幕| 久久97久久精品| 精品人妻一区二区三区麻豆| 亚洲精品自拍成人| 男男h啪啪无遮挡| 国产亚洲精品第一综合不卡| 国产精品 国内视频| 亚洲欧洲日产国产| 在线亚洲精品国产二区图片欧美| 色94色欧美一区二区| 久久久久视频综合| 久久精品久久久久久久性| 国产在线免费精品| 中文字幕色久视频| 性色av一级| 伦理电影大哥的女人| 国产97色在线日韩免费| 亚洲三区欧美一区| av有码第一页| 亚洲欧洲日产国产| 国产爽快片一区二区三区| 你懂的网址亚洲精品在线观看| 国产精品三级大全| 九九爱精品视频在线观看| 大陆偷拍与自拍| 在线观看三级黄色| 欧美亚洲日本最大视频资源| 久久人妻熟女aⅴ| 伦理电影免费视频| 亚洲av成人不卡在线观看播放网 | 老司机影院毛片| 精品久久蜜臀av无| 激情视频va一区二区三区| 女性生殖器流出的白浆| 国产成人啪精品午夜网站| 热re99久久精品国产66热6| 美女高潮到喷水免费观看| 中国三级夫妇交换| 建设人人有责人人尽责人人享有的| 国产成人精品福利久久| 国产精品人妻久久久影院| 亚洲欧美一区二区三区国产| 色94色欧美一区二区| 亚洲综合色网址| 欧美精品亚洲一区二区| 高清视频免费观看一区二区| 免费高清在线观看视频在线观看| 亚洲国产av影院在线观看| 狂野欧美激情性bbbbbb| 亚洲欧洲精品一区二区精品久久久 | 精品久久久久久电影网| 久久久精品94久久精品| 在现免费观看毛片| 国产在线一区二区三区精| 岛国毛片在线播放| 只有这里有精品99| √禁漫天堂资源中文www| 国产在视频线精品| av线在线观看网站| 在线天堂中文资源库| 男人操女人黄网站| 观看美女的网站| 日韩中文字幕欧美一区二区 | 黄色视频不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久久免费视频了| 免费日韩欧美在线观看| 欧美日本中文国产一区发布| 熟女少妇亚洲综合色aaa.| 国产无遮挡羞羞视频在线观看| 大陆偷拍与自拍| 丰满乱子伦码专区| 黑人巨大精品欧美一区二区蜜桃| 色精品久久人妻99蜜桃| 在线观看免费高清a一片| 天天影视国产精品| 午夜老司机福利片| 国产乱来视频区| 国产成人精品在线电影| 久久久精品94久久精品| 亚洲av成人不卡在线观看播放网 | 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区国产| 精品国产露脸久久av麻豆| 中国国产av一级| 熟妇人妻不卡中文字幕| 久久性视频一级片| 亚洲精品日本国产第一区| 另类亚洲欧美激情| 国产女主播在线喷水免费视频网站| 99精国产麻豆久久婷婷| 少妇精品久久久久久久| 免费在线观看完整版高清| 亚洲av国产av综合av卡| 国产精品欧美亚洲77777| 下体分泌物呈黄色| 日韩精品免费视频一区二区三区| 纵有疾风起免费观看全集完整版| 亚洲国产欧美在线一区| 欧美 亚洲 国产 日韩一| 丝袜喷水一区| 国产成人精品无人区| 色综合欧美亚洲国产小说| 欧美人与性动交α欧美软件| 国产精品久久久av美女十八| 国产免费福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 伊人久久国产一区二区| 欧美日韩视频精品一区| 欧美日韩av久久| 午夜福利免费观看在线| 观看美女的网站| 国产免费又黄又爽又色| 亚洲欧美成人综合另类久久久| 国产 精品1| 国产精品一区二区在线观看99| 成人影院久久| 欧美97在线视频| 国产欧美日韩综合在线一区二区| 久久久精品区二区三区| 久久天躁狠狠躁夜夜2o2o | videos熟女内射| 欧美亚洲 丝袜 人妻 在线| 午夜福利影视在线免费观看| 久久久久久久久久久免费av| 久久影院123| 天天影视国产精品| 69精品国产乱码久久久| 精品卡一卡二卡四卡免费| 国产精品一国产av| 色婷婷av一区二区三区视频| xxxhd国产人妻xxx| 亚洲av在线观看美女高潮| 日韩,欧美,国产一区二区三区| 男人操女人黄网站| 亚洲成国产人片在线观看| 一级黄片播放器| 欧美激情 高清一区二区三区| 国产精品一国产av| 日日摸夜夜添夜夜爱| 如何舔出高潮| av.在线天堂| 国产激情久久老熟女| 欧美成人午夜精品| 啦啦啦 在线观看视频| 精品一区二区免费观看| 免费在线观看视频国产中文字幕亚洲 | 欧美精品一区二区大全| 夫妻午夜视频| 亚洲精品在线美女| 色94色欧美一区二区| 一边亲一边摸免费视频| 青春草国产在线视频| 99国产综合亚洲精品| 亚洲国产欧美在线一区| 亚洲国产成人一精品久久久| 午夜av观看不卡| 亚洲欧美一区二区三区国产| 交换朋友夫妻互换小说| 老司机深夜福利视频在线观看 | 亚洲成av片中文字幕在线观看| 中文精品一卡2卡3卡4更新| 老熟女久久久| 亚洲欧美成人综合另类久久久| 久久久精品区二区三区| 免费高清在线观看视频在线观看| 一二三四在线观看免费中文在| 亚洲成人av在线免费| 亚洲色图综合在线观看| 捣出白浆h1v1| 久久人人爽av亚洲精品天堂| 宅男免费午夜| 久久久久人妻精品一区果冻| 久热这里只有精品99| av片东京热男人的天堂| 97在线人人人人妻| 精品一区二区免费观看| 999久久久国产精品视频| 亚洲精品,欧美精品| 一边摸一边抽搐一进一出视频| 欧美日韩综合久久久久久| 啦啦啦在线免费观看视频4| 久久久欧美国产精品| xxx大片免费视频| 十八禁高潮呻吟视频| 操出白浆在线播放| 国产无遮挡羞羞视频在线观看| 久久免费观看电影| 777久久人妻少妇嫩草av网站| 精品国产乱码久久久久久男人| 美女福利国产在线| 大香蕉久久网| 一二三四中文在线观看免费高清| 国产一区二区三区综合在线观看| 久久人人爽av亚洲精品天堂| 国产成人一区二区在线| 青青草视频在线视频观看| 男女无遮挡免费网站观看| 成人国语在线视频| 精品人妻一区二区三区麻豆| 国产精品国产三级专区第一集| 一二三四中文在线观看免费高清| 亚洲精品美女久久久久99蜜臀 | 91老司机精品| 91精品三级在线观看| 亚洲国产精品一区三区| 尾随美女入室| 国产成人av激情在线播放| 国产又色又爽无遮挡免| 制服人妻中文乱码| 免费少妇av软件| 女人精品久久久久毛片| 一个人免费看片子| 久久久久国产一级毛片高清牌| 国产日韩欧美亚洲二区| 国产成人免费观看mmmm| 天天躁狠狠躁夜夜躁狠狠躁| 可以免费在线观看a视频的电影网站 | 国产 一区精品| 亚洲三区欧美一区| 中文欧美无线码| 亚洲综合精品二区| av有码第一页| 精品第一国产精品| 国产色婷婷99| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩视频高清一区二区三区二| 天堂中文最新版在线下载| 美女主播在线视频| 黄色视频在线播放观看不卡| 一级毛片黄色毛片免费观看视频| 午夜福利视频精品| 亚洲国产精品999| 亚洲国产中文字幕在线视频| av线在线观看网站| 十八禁人妻一区二区| 99久久精品国产亚洲精品| 国产精品一二三区在线看| 男女免费视频国产| 久久久久人妻精品一区果冻| 在线观看免费视频网站a站| 国产欧美日韩综合在线一区二区| 少妇 在线观看| 久久精品国产亚洲av高清一级| 久久精品国产综合久久久| 国产精品国产av在线观看| 中文欧美无线码| 不卡av一区二区三区| 最近2019中文字幕mv第一页| 乱人伦中国视频| 国产成人精品久久久久久| 久久综合国产亚洲精品| 在线观看人妻少妇| 亚洲美女黄色视频免费看| 午夜激情av网站| 国产成人精品在线电影| 国产高清国产精品国产三级| 高清在线视频一区二区三区| 汤姆久久久久久久影院中文字幕| av不卡在线播放| 亚洲第一青青草原| 99久久精品国产亚洲精品| 久久精品久久久久久久性| 亚洲av日韩在线播放| 黄色毛片三级朝国网站| 免费看不卡的av| 中文字幕色久视频| 熟妇人妻不卡中文字幕| 亚洲精品久久久久久婷婷小说| 亚洲欧美清纯卡通| 久久鲁丝午夜福利片| 2018国产大陆天天弄谢| 99久久人妻综合| 亚洲色图 男人天堂 中文字幕| 亚洲欧美中文字幕日韩二区| 色94色欧美一区二区| 丰满乱子伦码专区| 久久久国产欧美日韩av| 中文字幕人妻熟女乱码| 亚洲精品成人av观看孕妇| 精品免费久久久久久久清纯 | 嫩草影院入口| 一级毛片黄色毛片免费观看视频| 美女国产高潮福利片在线看| 色精品久久人妻99蜜桃| 超碰97精品在线观看| 日本色播在线视频| 老熟女久久久| 国产精品人妻久久久影院| 日韩精品免费视频一区二区三区| 久久影院123| 久久精品国产亚洲av涩爱| 久久 成人 亚洲| 老司机影院毛片| 婷婷色综合www| 制服丝袜香蕉在线| 亚洲精品第二区| 黄片无遮挡物在线观看| 日韩制服丝袜自拍偷拍| 国产精品人妻久久久影院| 秋霞伦理黄片| 色视频在线一区二区三区| 九草在线视频观看| 精品久久久精品久久久| 欧美精品人与动牲交sv欧美| 久久精品亚洲熟妇少妇任你| 免费高清在线观看日韩| 观看美女的网站| h视频一区二区三区| 国产97色在线日韩免费| 欧美中文综合在线视频| 一区二区三区四区激情视频| 巨乳人妻的诱惑在线观看| 国产熟女午夜一区二区三区| 在线免费观看不下载黄p国产| 十八禁网站网址无遮挡| 国产视频首页在线观看| 涩涩av久久男人的天堂| 午夜福利网站1000一区二区三区| 日本av手机在线免费观看| 日韩视频在线欧美| 欧美精品一区二区免费开放| 中文乱码字字幕精品一区二区三区| 性高湖久久久久久久久免费观看| 久久热在线av| av视频免费观看在线观看| 大片电影免费在线观看免费| 高清在线视频一区二区三区| 精品久久蜜臀av无| 少妇被粗大的猛进出69影院| 天美传媒精品一区二区| 成人亚洲欧美一区二区av| www日本在线高清视频| 久久天躁狠狠躁夜夜2o2o | 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品国产欧美久久久 | 成年av动漫网址| 成人国语在线视频| 国产又色又爽无遮挡免| 国产免费视频播放在线视频| 国产免费福利视频在线观看| 麻豆精品久久久久久蜜桃| 少妇被粗大猛烈的视频| 午夜影院在线不卡| 久久热在线av| 成人亚洲欧美一区二区av| a级毛片黄视频| 大香蕉久久成人网| 九九爱精品视频在线观看| 国产一区亚洲一区在线观看| 亚洲,欧美精品.| 中文字幕av电影在线播放| 大陆偷拍与自拍| 97人妻天天添夜夜摸| 欧美精品一区二区免费开放| 涩涩av久久男人的天堂| 免费av中文字幕在线| 亚洲国产日韩一区二区| 亚洲四区av| 免费高清在线观看视频在线观看| 黄色 视频免费看| 免费高清在线观看视频在线观看| 亚洲三区欧美一区| 菩萨蛮人人尽说江南好唐韦庄| 日日啪夜夜爽| 一级毛片 在线播放| 十八禁高潮呻吟视频| 亚洲国产中文字幕在线视频| 看免费成人av毛片| 99久国产av精品国产电影| 视频在线观看一区二区三区| 自线自在国产av| 高清av免费在线| 午夜福利,免费看| 9191精品国产免费久久| 老鸭窝网址在线观看| 婷婷色av中文字幕| 一本—道久久a久久精品蜜桃钙片| 午夜免费男女啪啪视频观看| 久久韩国三级中文字幕| 电影成人av| 亚洲国产成人一精品久久久| 啦啦啦在线观看免费高清www| 国产精品二区激情视频| 欧美人与性动交α欧美软件| 街头女战士在线观看网站| 伊人久久国产一区二区| 亚洲精品国产色婷婷电影| 亚洲av日韩精品久久久久久密 | 久久久久久人人人人人| 又粗又硬又长又爽又黄的视频| 国产精品二区激情视频| 欧美 亚洲 国产 日韩一| 久久精品亚洲av国产电影网| 国产亚洲av片在线观看秒播厂| 中文字幕最新亚洲高清| 午夜福利免费观看在线| 久久久久国产精品人妻一区二区| 少妇人妻 视频| 久久精品aⅴ一区二区三区四区| 中文字幕色久视频| 国产成人午夜福利电影在线观看| 大片免费播放器 马上看| 色精品久久人妻99蜜桃| 亚洲图色成人| svipshipincom国产片| 国产欧美亚洲国产| a级毛片黄视频| 欧美日韩一级在线毛片| 成年动漫av网址| 黄色毛片三级朝国网站| 电影成人av| 亚洲美女搞黄在线观看| 男女高潮啪啪啪动态图| 看十八女毛片水多多多| 我要看黄色一级片免费的| 欧美精品一区二区大全| 久热爱精品视频在线9| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区欧美精品| 曰老女人黄片| 日韩av在线免费看完整版不卡| 欧美成人午夜精品| 国产精品成人在线| 婷婷色av中文字幕| 哪个播放器可以免费观看大片| 亚洲精品在线美女| 69精品国产乱码久久久| 免费黄网站久久成人精品| 日韩精品有码人妻一区| 伊人久久大香线蕉亚洲五| 精品国产一区二区久久| 建设人人有责人人尽责人人享有的| 伊人亚洲综合成人网| 亚洲五月色婷婷综合| 国产精品偷伦视频观看了| 国产精品一二三区在线看| 亚洲av中文av极速乱| 久久久精品区二区三区| 国产精品一国产av| 18禁动态无遮挡网站| 少妇人妻 视频| av在线app专区| 又粗又硬又长又爽又黄的视频| 午夜福利乱码中文字幕| 亚洲欧洲日产国产| 日韩一区二区视频免费看| 国产片内射在线| 亚洲,欧美,日韩| 啦啦啦视频在线资源免费观看| 夫妻性生交免费视频一级片| 嫩草影院入口| av.在线天堂| 国产精品久久久久久精品电影小说| 亚洲av电影在线进入| 又黄又粗又硬又大视频| 成年美女黄网站色视频大全免费| 五月天丁香电影| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美清纯卡通| 亚洲一区二区三区欧美精品| 超碰97精品在线观看| 少妇的丰满在线观看| 19禁男女啪啪无遮挡网站| 国产精品免费视频内射| 最新在线观看一区二区三区 | 国产伦理片在线播放av一区| 99久国产av精品国产电影| 看非洲黑人一级黄片| 欧美变态另类bdsm刘玥| 婷婷色综合www| 最黄视频免费看| 国产精品国产av在线观看| 我要看黄色一级片免费的| 看非洲黑人一级黄片| 亚洲视频免费观看视频| 2021少妇久久久久久久久久久| 久久精品久久久久久噜噜老黄| 考比视频在线观看| 中文字幕精品免费在线观看视频| 欧美日韩一级在线毛片| 蜜桃在线观看..| 老司机影院成人| 精品国产一区二区久久| 亚洲成av片中文字幕在线观看| 丝袜脚勾引网站| 少妇人妻久久综合中文| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美中文字幕日韩二区| 日日撸夜夜添| 19禁男女啪啪无遮挡网站| 韩国av在线不卡| av又黄又爽大尺度在线免费看| 色94色欧美一区二区| 婷婷色av中文字幕| 久久久精品94久久精品| 亚洲欧美一区二区三区黑人|