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

    基于RNA-Seq對(duì)心肌肌鈣蛋白IR193H突變限制性心肌病小鼠轉(zhuǎn)錄組的分析

    2019-01-17 02:40:56干意封葉趙唯安劉振國(guó)張蕾田杰
    解放軍醫(yī)學(xué)雜志 2018年11期
    關(guān)鍵詞:差異基因測(cè)序顯著性

    干意,封葉,趙唯安,劉振國(guó),張蕾,田杰

    限制性心肌病(restrictive cardiomyopathy,RCM)是兒童心肌病的一種類(lèi)型,占小兒心肌病的2%~5%[1],是一種由于心肌僵硬度升高導(dǎo)致以舒張功能?chē)?yán)重受損為主要特征的心肌病,表現(xiàn)為心室舒張末容積正?;蚩s小,心室壁厚度正常或輕度增加,而收縮功能大多正常或僅有輕度受損[2-3]。兒童RCM多為特發(fā)性,病因不明,約30%病例有家族發(fā)病傾向,多為常染色體顯性遺傳。已有研究發(fā)現(xiàn)心肌肌小節(jié)蛋白、結(jié)蛋白、心肌肌鈣蛋白I(cardiac troponin I,cTnI)、β-肌球蛋白重鏈等基因突變與RCM有關(guān)[4],其中,cTnI作為重要的心肌結(jié)構(gòu)蛋白,其突變與原發(fā)性心肌病關(guān)系尤為密切[5-6]。研究發(fā)現(xiàn),6種類(lèi)型的cTnI突變可導(dǎo)致RCM(L144Q、R145W、A171T、K178E、D190G、R192H),其中兩種類(lèi)型的突變(K178E、R192H)臨床癥狀最為嚴(yán)重[7]。本課題組擁有黃旭培教授成功構(gòu)建的來(lái)源于人cTnI基因突變的RCM小鼠動(dòng)物模型(在小鼠系193位點(diǎn),cTnI R193H),該模型小鼠與攜帶該突變的RCM患者(在人192位點(diǎn),cTnI R192H)有相似的病理改變、癥狀和表型[8-9],其在6~8周時(shí),超聲即可檢測(cè)到舒張功能的改變[10],但具體致病機(jī)制尚不明確。本研究利用高通量測(cè)序的方法,分析RCM小鼠與正常小鼠差異表達(dá)的基因,篩選出可能參與致病過(guò)程的分子及信號(hào)通路,為RCM發(fā)病機(jī)制研究提供新的思路和方向。

    1 材料與方法

    1.1 主要試劑和設(shè)備 R193H(美國(guó)佛羅里達(dá)大西洋大學(xué)黃旭培教授贈(zèng)送)C57和野生型(WT)C57 3月齡小鼠(購(gòu)自重慶醫(yī)科大學(xué)實(shí)驗(yàn)動(dòng)物中心,飼養(yǎng)于重慶醫(yī)科大學(xué)附屬兒童醫(yī)院SPF動(dòng)物中心),RNA提取試劑盒(BioTeke,北京),反轉(zhuǎn)錄試劑盒(TaKaRa,日本),SYBR Green試劑盒(KAPA,美國(guó))。實(shí)驗(yàn)室常用設(shè)備由重慶醫(yī)科大學(xué)兒研所提供。動(dòng)物實(shí)驗(yàn)符合倫理要求,并經(jīng)重慶醫(yī)科大學(xué)醫(yī)學(xué)研究倫理委員會(huì)批準(zhǔn)。

    1.2 小鼠心臟組織總RNA提取 準(zhǔn)備3月齡R193H和WT小鼠,利用CO2麻醉處死后,剪開(kāi)小鼠胸腔,灌洗心臟,取心臟組織,利用高純度總RNA快速提取試劑盒提取心臟組織總RNA??俁NA樣本質(zhì)量由上海裂冰生物醫(yī)藥有限公司采用Agilent2200儀器進(jìn)行檢測(cè)。RNA分子完整數(shù)(RNA integrity number,RIN)>8.0的RNA用于cDNA文庫(kù)構(gòu)建。樣本和RNA均放入–80℃冰箱保存。

    1.3 基于RNA-Seq技術(shù)的差異表達(dá)基因篩選基因的表達(dá)量計(jì)算主要采用RPKM方法(Reads per kb per million reads),計(jì)算公式為:RPKM=106C/(NL/103)。其中RPKM為基因表達(dá)量,C為唯一比對(duì)到該基因上讀段(reads)數(shù),N為唯一比對(duì)到參考基因上的總讀段數(shù),L為該基因編碼區(qū)的堿基數(shù)。RPKM法消除了基因長(zhǎng)度和測(cè)序量差異對(duì)計(jì)算基因表達(dá)的影響,計(jì)算得到的基因表達(dá)量可直接用于比較不同樣品間的基因表達(dá)差異[11]。采用國(guó)際公認(rèn)算法EBSeq軟件包對(duì)表達(dá)量(Counts)進(jìn)行差異篩選。以差異倍數(shù)(fold change,F(xiàn)C)≥1.2或≤0.833、P<0.05為標(biāo)準(zhǔn)篩選差異表達(dá)基因,其中錯(cuò)誤檢出率(false discovery rate,F(xiàn)DR)<0.05。利用Cluster軟件,以歐氏距離為距離距陣計(jì)算公式,對(duì)差異表達(dá)基因和樣本同時(shí)進(jìn)行分層聚類(lèi)分析,聚類(lèi)結(jié)果用Java Treeview 1.1.6 R2軟件顯示。

    1.4 基因功能分析 將分析得到的差異基因基于數(shù)據(jù)庫(kù)NCBI(http://www.ncbi.nlm.nih.gov/)、UniProt(http://www.uniprot.org/)、GOC(http://www.geneontology.org/),分別從生物學(xué)過(guò)程(biological process,BP)、細(xì)胞組分(celluar components,CC)、分子功能(molecular function,MF)這三個(gè)層面進(jìn)行基因本體(gene ontology,GO)注釋?zhuān)玫交騾⑴c的所有GO,采用Fisher檢驗(yàn)計(jì)算每個(gè)GO的顯著性水平(P值),從而篩選出差異基因富集的顯著性GO。其中,P<0.05的值用灰白相間的表格標(biāo)出。

    1.5 功能調(diào)控網(wǎng)絡(luò)構(gòu)建 基于GO的層次結(jié)構(gòu),利用NovelBio云計(jì)算平臺(tái)將所有GO之間的相互調(diào)控及從屬關(guān)系整理成數(shù)據(jù)庫(kù),構(gòu)建功能關(guān)系網(wǎng)絡(luò),總結(jié)實(shí)驗(yàn)影響的功能群體,以及顯著性功能的內(nèi)在從屬關(guān)系。以差異基因所做GO分析中的顯著性GO條目(P<0.05)為研究對(duì)象進(jìn)行功能調(diào)控分析,構(gòu)建功能調(diào)控網(wǎng)絡(luò)。

    1.6 信號(hào)通路分析(pathway-analysis) 將篩選出的差異基因于KEGG(Kyoto encyclopedia of genes and genomes)數(shù)據(jù)庫(kù)進(jìn)行通路注釋?zhuān)玫讲町惢騾⑴c的所有通路條目,采用Fisher檢驗(yàn)計(jì)算通路的顯著性水平(P值),從而篩選出差異基因富集的顯著性通路條目。P<0.05的值用灰白相間的表格標(biāo)出。

    1.7 信號(hào)通路相互作用關(guān)系分析 選用顯著性通路,根據(jù)KEGG數(shù)據(jù)庫(kù)中通路之間的關(guān)聯(lián)關(guān)系整合成顯著性通路之間的信號(hào)轉(zhuǎn)導(dǎo)網(wǎng)絡(luò)。在此以差異基因所做的通路分析中的顯著性通路條目(P<0.05)進(jìn)行信號(hào)通路相互關(guān)系網(wǎng)絡(luò)構(gòu)建。

    1.8 引物的設(shè)計(jì)和合成 在NCBI的GenBank數(shù)據(jù)庫(kù)(https://www.ncbi.nlm.nih.gov/genbank/)中檢索Tifab、Ifi202b、Ctss、Itgb2、Pik3ap1、Arhgap21、Myocd、MT-ND5基因的mRNA序列。采用Primer Premier 5.0軟件針對(duì)這些基因CDS核心編碼區(qū)設(shè)計(jì)特異性引物,引物的堿基序列由上海生工生物工程股份有限公司合成。引物序列為:Tifab,正向5'-GCTGCCACAGTTGTCTCGCTAC-3',反向5'-ATG CCAGAGAAGGAGATTCTGTTGATG-3';Ifi202b,正向5'-ACTGACCAGTGTCAACTGTGAGATTG-3',反向5'-TGACCTCCATGTAACTGTACCTCGTAG-3';Ctss,正向5'-GCCAGCCAヰCCTCCヰCヰCヰC-3',反向5'-TCTTGCCATCCGAATGTATCCTTGATC-3';Itgb2,正向5'-GCAGCAGAAGGACGGAAGGAA C-3',反向5'-ATGACCAGGAGGAGGACACCAATC-3';Pik3ap1,正向5'-GCAGCCAACCCAGTACAGTT-3',反向5'-GGGACAAATCCAGCCATAGA-3';Arhgap21,正向5'-GAGAGCAGACAACGCCATCAGAAG-3',反向5'-TCCGCCTCTTGGTCTCACTGTAC-3';Myocd,正向5'-CTGCCAGACACCヰCACCGATG-3',反向5'-CCACCAGCATCTTGTCCTTCTCAG-3';MT-ND5,正向5'-CCTGGCAGACGAACAAGACATC C-3',反向5'-CGAGGCヰCCGAヰACTAGGCATG-3'。

    1.9 Q-PCR檢測(cè)基因表達(dá)量 為驗(yàn)證RNA-Seq結(jié)果的可靠性,隨機(jī)選擇8個(gè)差異表達(dá)基因進(jìn)行Q-PCR檢測(cè),分別從R193H和WT小鼠中提取總RNA,并利用反轉(zhuǎn)錄試劑盒(TaKaRa,日本)將其反轉(zhuǎn)錄成cDNA,利用SYBR Green試劑盒(Qiagen,德國(guó))進(jìn)行Real-time PCR,以β-actin為內(nèi)參基因,利用實(shí)時(shí)熒光定量PCR儀(Bio-Rad,美國(guó))測(cè)量吸光度,mRNA相對(duì)表達(dá)量分析采用2–ΔΔCt法。實(shí)驗(yàn)數(shù)據(jù)采用GraphPad Prism 6軟件進(jìn)行獨(dú)立樣本t檢驗(yàn)并制圖。P<0.05為差異有統(tǒng)計(jì)學(xué)意義。

    2 結(jié) 果

    2.1 差異表達(dá)基因篩選 對(duì)R193H和WT兩組共6個(gè)樣本(生物學(xué)重復(fù)3個(gè))進(jìn)行分析,得出差異表達(dá)的基因共52個(gè)(表1),其中28個(gè)基因上調(diào),24個(gè)基因下調(diào)。對(duì)差異表達(dá)基因和樣本同時(shí)進(jìn)行分層聚類(lèi)分析(圖1),可見(jiàn)表達(dá)模式相似的基因具有功能相關(guān)性。根據(jù)差異表達(dá)的基因構(gòu)建火山圖,對(duì)樣本之間和基因之間的差異進(jìn)行直觀(guān)展示(圖2)。

    表1 RNA-Seq分析野生型和突變型R193H小鼠差異表達(dá)基因Tab.1 Analysis of the differentially expressed genes of wild type and mutant R193H mice by RNA-Seq

    圖1 R193H與WT基因差異表達(dá)模式聚類(lèi)圖Fig.1 Clustering diagram of R193H vs. WT gene expression pattern

    圖2 R193H與WT基因差異表達(dá)火山圖Fig.2 R193H vs. WT genetic differentially expressed volcanic maps

    2.2 利用生物信息學(xué)方法進(jìn)行的基因功能分析從BP、CC、MF這3個(gè)層面進(jìn)行GO注釋(圖3),基于通路的分析進(jìn)一步了解基因的生物學(xué)功能。將差異基因進(jìn)行GO顯著性富集(圖4),差異表達(dá)基因在生物過(guò)程分組中主要富集在翻譯延伸(translational elongation)、固有免疫應(yīng)答(innate immune response)、ATP合成耦合電子傳輸(ATP synthesis coupled electron transport)這三個(gè)生物學(xué)過(guò)程中,涉及轉(zhuǎn)錄翻譯、免疫以及能量代謝相關(guān)基因。層次樹(shù)狀圖顯示出顯著性功能的內(nèi)在從屬關(guān)系(圖5)。

    圖3 GO功能分類(lèi)Fig.3 GO function classification

    圖4 KEGG富集散點(diǎn)圖Fig.4 KEGG enriched scatter plot

    圖5 顯著性功能間的層次樹(shù)形關(guān)系圖Fig.5 The hierarchical tree diagram for significant functions

    2.3 顯著性差異基因的KEGG分析 根據(jù)差異表達(dá)的KEGG功能注釋結(jié)果進(jìn)行通路顯著性富集分析(圖6),共得到16條顯著富集的通路,涉及了氨?;?tRNA生物合成通路(aminoacyl-tRNA biosynthesis)、NF-κB信號(hào)通路(NF-kappa B signaling pathway)、Toll樣受體信號(hào)通路(Toll-like receptor signaling pathway)等。對(duì)顯著性差異的通路條目進(jìn)行調(diào)控網(wǎng)絡(luò)構(gòu)建,可以直觀(guān)顯示顯著性功能之間的關(guān)系(圖7)。

    圖6 顯著通路富集度柱狀圖Fig.6 Significant pathway enrichment histogram

    圖7 顯著性功能間的網(wǎng)絡(luò)構(gòu)建關(guān)系圖Fig.7 Network construction diagram of significant functions

    2.4 差異表達(dá)基因驗(yàn)證 隨機(jī)選取8個(gè)差異表達(dá)的基因,每組取4個(gè)樣本進(jìn)行實(shí)時(shí)熒光定量PCR驗(yàn)證,分別為T(mén)ifab、Ifi202b、Ctss、Itgb2、Pik3ap1、Arhgap21、Myocd、MT-ND5。結(jié)果顯示,突變組Tifab、Ifi202b、Ctss、Itgb2、Pik3ap1基因mRNA表達(dá)量與WT組比較明顯升高,Arhgap21、Myocd、MT-ND5基因mRNA表達(dá)量與WT組比較明顯降低(P<0.05),該結(jié)果與測(cè)序結(jié)果一致(圖8)。

    圖8 野生型和突變型相關(guān)基因的mRNA水平Fig.8 mRNA levels of wild-type and mutant

    3 討 論

    目前研究認(rèn)為,RCM可能是特發(fā)性的、家族性的或繼發(fā)于其他條件發(fā)生的疾病[12]。最新分子遺傳技術(shù)的進(jìn)展明顯加深了對(duì)該疾病的理解[5],在編碼肌節(jié)蛋白的相關(guān)基因中已經(jīng)鑒定了幾個(gè)RCM特異性突變,其中cTnI作為重要的心肌結(jié)構(gòu)蛋白,其突變可以導(dǎo)致RCM和家族性肥厚性心肌病[13]。本課題組擁有來(lái)源于人cTnI基因192位點(diǎn)突變的RCM小鼠動(dòng)物模型(在小鼠系193位點(diǎn),cTnI R193H),模型小鼠與攜帶該突變的RCM患者有相似的病理改變、癥狀和表型[10,14],從而驗(yàn)證了cTnI突變與RCM的關(guān)系。目前,RCM相關(guān)機(jī)制的研究頗多。有研究表明,cTnI突變可能引起心肌細(xì)胞鈣敏性改變而導(dǎo)致RCM[15],也有學(xué)者認(rèn)為cTnI突變可能導(dǎo)致與其相互作用的基因或蛋白發(fā)生改變,從而導(dǎo)致心臟功能的改變[16],但cTnI突變引起RCM的機(jī)制尚不明確。

    轉(zhuǎn)錄組是特定組織或細(xì)胞在某一發(fā)育階段或功能狀態(tài)下轉(zhuǎn)錄出來(lái)的所有RNA的集合。RNA-Seq利用高通量測(cè)序技術(shù)對(duì)組織或細(xì)胞中所有RNA反轉(zhuǎn)錄而成的cDNA文庫(kù)進(jìn)行測(cè)序,通過(guò)統(tǒng)計(jì)相關(guān)reads數(shù)計(jì)算出不同RNA的表達(dá)量[17],尋找出兩組樣本中差異表達(dá)的基因[18]。這種方法改變了轉(zhuǎn)錄序列的研究方式,主要應(yīng)用于轉(zhuǎn)錄組學(xué)的研究。本課題組利用RNA-Seq技術(shù)探究cTnI在RCM發(fā)病機(jī)制中的生物學(xué)作用[19]。通過(guò)檢測(cè)cTnI突變導(dǎo)致的RCM模型鼠和同月齡WT小鼠基因表達(dá)譜的差異,找出52個(gè)差異表達(dá)基因,隨機(jī)挑選出8個(gè)差異表達(dá)基因利用Q-PCR進(jìn)行驗(yàn)證,其結(jié)果初步驗(yàn)證了測(cè)序結(jié)果的準(zhǔn)確性。對(duì)差異表達(dá)的基因進(jìn)行GO分析發(fā)現(xiàn),上調(diào)和下調(diào)的基因在細(xì)胞組分、分子功能、生物過(guò)程三個(gè)方面差異有統(tǒng)計(jì)學(xué)意義。其中上調(diào)的基因及產(chǎn)物主要分布于細(xì)胞表面和質(zhì)膜外側(cè)面,分子水平主要是具有結(jié)合活性,多參與到固有免疫應(yīng)答的過(guò)程,涉及的基因有Cxcr3、Fcgr1、Icam1、Clec12a、Ctss、Lbp等。而下調(diào)的基因及產(chǎn)物主要分布于細(xì)胞骨架和線(xiàn)粒體中,分子功能主要是催化活性,參與的生物學(xué)過(guò)程主要是翻譯延長(zhǎng)和ATP合成耦聯(lián)電子運(yùn)輸,涉及的基因有D2hgdh、ND5、Myom1、Ubr4、Spta1、Lrpprc、Arhgap21。這一結(jié)果表明,cTnI突變可能伴隨免疫相關(guān)分子的上調(diào)和能量代謝的抑制。通過(guò)通路顯著性富集分析,共得到16條顯著富集的通路,涉及氨?;?tRNA生物合成、NF-κB信號(hào)通路、Toll樣受體信號(hào)通路。但是由于差異基因的數(shù)量并不大,通路分析的結(jié)果無(wú)統(tǒng)計(jì)學(xué)意義。因此,本研究后續(xù)重點(diǎn)分析了差異基因參與的生物學(xué)進(jìn)程及參與的信號(hào)通路。通過(guò)分析該測(cè)序結(jié)果發(fā)現(xiàn),差異表達(dá)的基因主要參與免疫應(yīng)答、線(xiàn)粒體能量代謝、心肌細(xì)胞分化等過(guò)程,如Ifi202b、ND5、Itgb2、Myocd等基因。既往研究表明,cTnI某些序列是免疫調(diào)節(jié)序列,能通過(guò)免疫炎癥反應(yīng)導(dǎo)致擴(kuò)張性心肌疾病[20],而本研究結(jié)果顯示,在cTnI突變所致的RCM動(dòng)物模型中參與免疫應(yīng)答的一系列基因表達(dá)量發(fā)生改變,提示免疫功能的改變可能在RCM的發(fā)病中具有重要作用。另有研究表明,通過(guò)轉(zhuǎn)運(yùn)蛋白和結(jié)合蛋白之間復(fù)雜的相互作用,調(diào)節(jié)Ca2+信號(hào)通路來(lái)調(diào)節(jié)線(xiàn)粒體的代謝活性,從而導(dǎo)致能量代謝受損可能是RCM發(fā)生的潛在機(jī)制[21]。本研究發(fā)現(xiàn),線(xiàn)粒體能量代謝相關(guān)基因差異表達(dá)顯著,提示在RCM發(fā)病中cTnI突變可能伴隨著能量代謝改變。NADH泛醌氧化還原酶亞基5(ND5或MT-ND5,EC 1.6.5.3)由位于線(xiàn)粒體DNA的鳥(niǎo)嘌呤富集重鏈中的最大基因編碼(mtDNA坐標(biāo):12337-14148),參與的相關(guān)途徑包括呼吸電子傳遞,通過(guò)化學(xué)滲透耦合進(jìn)行ATP合成,以及通過(guò)解耦聯(lián)蛋白產(chǎn)生熱量,該基因的表達(dá)降低可能與線(xiàn)粒體能量代謝受損相關(guān)[22],但其具體內(nèi)在機(jī)制尚需深入研究。差異表達(dá)的基因中也發(fā)現(xiàn)有PI3KAP1基因,結(jié)合課題組前期研究證實(shí)cTnI與多種PI3K-Akt信號(hào)通路中的蛋白相結(jié)合[23],表明PI3K-Akt信號(hào)通路可能通過(guò)cTnI參與RCM的發(fā)生。

    綜上所述,利用測(cè)序技術(shù)和特有的RCM模型,本研究初步篩選出了可能參與RCM發(fā)病機(jī)制的基因及其相關(guān)通路,可以通過(guò)尋找對(duì)差異表達(dá)基因有調(diào)控作用的藥物,提高對(duì)該疾病的治療效果。cTnI作為一種重要的心肌蛋白,作用機(jī)制十分復(fù)雜,其功能仍需進(jìn)一步深入研究。

    猜你喜歡
    差異基因測(cè)序顯著性
    ICR鼠肝和腎毒性損傷生物標(biāo)志物的篩選
    杰 Sir 帶你認(rèn)識(shí)宏基因二代測(cè)序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    二代測(cè)序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    基于RNA 測(cè)序研究人參二醇對(duì)大鼠心血管內(nèi)皮細(xì)胞基因表達(dá)的影響 (正文見(jiàn)第26 頁(yè))
    基于顯著性權(quán)重融合的圖像拼接算法
    電子制作(2019年24期)2019-02-23 13:22:26
    基于視覺(jué)顯著性的視頻差錯(cuò)掩蓋算法
    一種基于顯著性邊緣的運(yùn)動(dòng)模糊圖像復(fù)原方法
    論商標(biāo)固有顯著性的認(rèn)定
    基因捕獲測(cè)序診斷血癌
    單細(xì)胞測(cè)序技術(shù)研究進(jìn)展
    嫩草影院精品99| 亚洲18禁久久av| 日韩亚洲欧美综合| 欧美一区二区国产精品久久精品| 非洲黑人性xxxx精品又粗又长| 日产精品乱码卡一卡2卡三| 桃色一区二区三区在线观看| 国产大屁股一区二区在线视频| 99久久精品一区二区三区| 久久综合国产亚洲精品| 精品熟女少妇av免费看| 男人的好看免费观看在线视频| 乱码一卡2卡4卡精品| 99久久中文字幕三级久久日本| 欧美激情国产日韩精品一区| 精品久久久久久久久久久久久| 国产片特级美女逼逼视频| 免费人成在线观看视频色| 一级毛片aaaaaa免费看小| 精品一区二区三区人妻视频| 97超碰精品成人国产| 春色校园在线视频观看| 亚洲av电影不卡..在线观看| 国产午夜精品论理片| 婷婷色av中文字幕| 亚洲高清免费不卡视频| 色综合站精品国产| 国产真实伦视频高清在线观看| 国产探花极品一区二区| 亚洲av电影不卡..在线观看| 国产av码专区亚洲av| 啦啦啦韩国在线观看视频| 偷拍熟女少妇极品色| 黄色日韩在线| 国产乱来视频区| 亚洲内射少妇av| 国产精品一区二区在线观看99 | 嘟嘟电影网在线观看| 欧美bdsm另类| 国产又色又爽无遮挡免| 亚洲自拍偷在线| 亚洲成av人片在线播放无| 综合色av麻豆| 熟妇人妻久久中文字幕3abv| 99九九线精品视频在线观看视频| 国产高潮美女av| 午夜免费男女啪啪视频观看| 国产成人aa在线观看| av在线天堂中文字幕| АⅤ资源中文在线天堂| 麻豆久久精品国产亚洲av| 久久久国产成人免费| 色哟哟·www| 国内精品美女久久久久久| 午夜精品在线福利| 日韩中字成人| 国产亚洲精品av在线| 欧美成人一区二区免费高清观看| 一级av片app| 亚洲18禁久久av| 可以在线观看毛片的网站| 日韩三级伦理在线观看| 免费看av在线观看网站| 欧美成人a在线观看| 黄片wwwwww| 搞女人的毛片| 亚洲av熟女| 午夜视频国产福利| 亚洲高清免费不卡视频| 国产黄色小视频在线观看| av免费观看日本| 婷婷色麻豆天堂久久 | 哪个播放器可以免费观看大片| 免费av不卡在线播放| 中文字幕久久专区| 啦啦啦啦在线视频资源| 99久久精品热视频| 老司机影院毛片| 99国产精品一区二区蜜桃av| 久久热精品热| 亚洲色图av天堂| 午夜久久久久精精品| 观看美女的网站| 亚洲aⅴ乱码一区二区在线播放| 欧美三级亚洲精品| 日韩成人av中文字幕在线观看| 搡女人真爽免费视频火全软件| 欧美日本亚洲视频在线播放| 国产精品一区二区三区四区久久| 欧美日韩精品成人综合77777| 人体艺术视频欧美日本| 亚洲欧美成人精品一区二区| 午夜久久久久精精品| av线在线观看网站| 欧美一区二区国产精品久久精品| 亚洲欧洲国产日韩| 国产极品精品免费视频能看的| 一级毛片aaaaaa免费看小| 老司机福利观看| 亚洲成人久久爱视频| 国产免费视频播放在线视频 | 亚洲国产精品国产精品| 欧美激情在线99| 热99在线观看视频| 纵有疾风起免费观看全集完整版 | 国产精品国产高清国产av| 婷婷色av中文字幕| 欧美bdsm另类| 男女那种视频在线观看| 婷婷六月久久综合丁香| 你懂的网址亚洲精品在线观看 | 久久这里只有精品中国| 中文字幕av成人在线电影| 亚洲精华国产精华液的使用体验| 久久精品久久久久久久性| 69人妻影院| 国产亚洲91精品色在线| 久久精品夜夜夜夜夜久久蜜豆| 我的老师免费观看完整版| 日韩欧美在线乱码| 1024手机看黄色片| 69人妻影院| 综合色丁香网| 黑人高潮一二区| 久久99蜜桃精品久久| 精品国产露脸久久av麻豆 | a级毛片免费高清观看在线播放| 国产伦精品一区二区三区视频9| 亚洲欧洲国产日韩| eeuss影院久久| 日韩一本色道免费dvd| 夫妻性生交免费视频一级片| 18+在线观看网站| 精品国内亚洲2022精品成人| 国产精品久久电影中文字幕| 青春草国产在线视频| 久久午夜福利片| 免费观看人在逋| 免费大片18禁| 久久精品久久精品一区二区三区| 男女边吃奶边做爰视频| 欧美日韩精品成人综合77777| 日韩av在线免费看完整版不卡| 国产毛片a区久久久久| 熟妇人妻久久中文字幕3abv| 国产精品人妻久久久久久| 久久精品熟女亚洲av麻豆精品 | 亚洲欧洲日产国产| 亚洲精品乱久久久久久| 日本与韩国留学比较| 精品国产露脸久久av麻豆 | 久久久精品欧美日韩精品| 成人毛片60女人毛片免费| 18+在线观看网站| 一级av片app| 免费黄色在线免费观看| 成人午夜高清在线视频| 国产黄色小视频在线观看| 久久精品影院6| 最新中文字幕久久久久| 少妇被粗大猛烈的视频| 日韩av在线免费看完整版不卡| 久久久精品大字幕| 非洲黑人性xxxx精品又粗又长| 国产精华一区二区三区| 国内精品宾馆在线| or卡值多少钱| 国产精品久久久久久av不卡| 最近视频中文字幕2019在线8| 99热网站在线观看| 久久久久久伊人网av| 午夜福利网站1000一区二区三区| 中文亚洲av片在线观看爽| 国产爱豆传媒在线观看| 国产伦在线观看视频一区| 大又大粗又爽又黄少妇毛片口| 久久人妻av系列| 精品人妻一区二区三区麻豆| 国产精品电影一区二区三区| 尾随美女入室| 精品久久久久久久人妻蜜臀av| a级毛色黄片| 成年av动漫网址| 少妇熟女欧美另类| 三级毛片av免费| 非洲黑人性xxxx精品又粗又长| 能在线免费看毛片的网站| 18禁裸乳无遮挡免费网站照片| 日本av手机在线免费观看| 在线观看一区二区三区| 色噜噜av男人的天堂激情| 国产成年人精品一区二区| 国产精品一区二区三区四区久久| 亚洲av成人av| 深夜a级毛片| 国产真实伦视频高清在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久午夜欧美精品| www.av在线官网国产| 九九爱精品视频在线观看| 日日干狠狠操夜夜爽| 久久精品国产自在天天线| 精品午夜福利在线看| 日韩成人av中文字幕在线观看| 综合色丁香网| 久久久久久久久大av| 哪个播放器可以免费观看大片| 汤姆久久久久久久影院中文字幕 | 丰满乱子伦码专区| 国产亚洲午夜精品一区二区久久 | 亚洲av男天堂| 爱豆传媒免费全集在线观看| 观看美女的网站| 国模一区二区三区四区视频| 麻豆乱淫一区二区| 国产精品久久电影中文字幕| 一二三四中文在线观看免费高清| 亚洲国产欧美人成| 欧美成人a在线观看| 久久国内精品自在自线图片| 亚洲va在线va天堂va国产| 久久精品国产亚洲网站| 国产真实乱freesex| 一区二区三区乱码不卡18| 少妇猛男粗大的猛烈进出视频 | 我的女老师完整版在线观看| 精品不卡国产一区二区三区| 国产精品麻豆人妻色哟哟久久 | 精华霜和精华液先用哪个| 国产午夜精品久久久久久一区二区三区| 男的添女的下面高潮视频| 国产伦在线观看视频一区| 国产亚洲精品av在线| 久久久久国产网址| 国产成人一区二区在线| 亚洲,欧美,日韩| 大香蕉久久网| 免费在线观看成人毛片| 国产精品国产三级国产av玫瑰| 七月丁香在线播放| 18禁在线播放成人免费| 久久久久性生活片| 18禁裸乳无遮挡免费网站照片| 免费观看的影片在线观看| 亚洲高清免费不卡视频| 国产成人a∨麻豆精品| 久久久久久久午夜电影| 偷拍熟女少妇极品色| 亚洲综合精品二区| 成年免费大片在线观看| 中文在线观看免费www的网站| 麻豆成人av视频| 欧美日韩一区二区视频在线观看视频在线 | 欧美性猛交黑人性爽| 97超视频在线观看视频| 最近视频中文字幕2019在线8| 亚洲四区av| 午夜精品国产一区二区电影 | 国产极品精品免费视频能看的| 精华霜和精华液先用哪个| 国产成人福利小说| 性插视频无遮挡在线免费观看| 国产精品熟女久久久久浪| 嘟嘟电影网在线观看| 欧美潮喷喷水| 国产熟女欧美一区二区| 国产一区二区在线观看日韩| 九九热线精品视视频播放| 亚洲欧美日韩高清专用| 日韩一区二区三区影片| 最近最新中文字幕大全电影3| 人妻夜夜爽99麻豆av| 精品国产露脸久久av麻豆 | 亚洲天堂国产精品一区在线| 亚洲人成网站在线观看播放| 成人亚洲欧美一区二区av| 久久久亚洲精品成人影院| 在线免费观看的www视频| 两个人视频免费观看高清| 久99久视频精品免费| 级片在线观看| 高清视频免费观看一区二区 | 麻豆成人av视频| 国产精品综合久久久久久久免费| 久久久久精品久久久久真实原创| 国产亚洲午夜精品一区二区久久 | 国产美女午夜福利| 婷婷色av中文字幕| 久久人妻av系列| 欧美区成人在线视频| 国产av在哪里看| 少妇丰满av| 国产精品日韩av在线免费观看| 最后的刺客免费高清国语| 精品久久久久久成人av| 久久久久免费精品人妻一区二区| 男人舔奶头视频| 天堂√8在线中文| 久久国内精品自在自线图片| 中文字幕人妻熟人妻熟丝袜美| 三级男女做爰猛烈吃奶摸视频| 男女下面进入的视频免费午夜| 天堂√8在线中文| 国国产精品蜜臀av免费| 偷拍熟女少妇极品色| 午夜激情欧美在线| 99久久精品一区二区三区| 国产精品一区二区性色av| av免费在线看不卡| 色综合站精品国产| 亚洲国产色片| 精品久久久久久久久久久久久| 亚洲,欧美,日韩| 亚洲精品日韩在线中文字幕| 免费av毛片视频| 欧美成人一区二区免费高清观看| 久久久国产成人免费| 亚洲人成网站在线播| 51国产日韩欧美| 欧美丝袜亚洲另类| 亚洲无线观看免费| eeuss影院久久| 黑人高潮一二区| 国产在视频线精品| 欧美xxxx性猛交bbbb| 欧美成人精品欧美一级黄| 国产成人一区二区在线| 欧美一区二区精品小视频在线| 免费播放大片免费观看视频在线观看 | 免费av不卡在线播放| 久久99蜜桃精品久久| 日本猛色少妇xxxxx猛交久久| 亚洲欧美日韩高清专用| 午夜久久久久精精品| 免费av不卡在线播放| 中文字幕av成人在线电影| 久久精品国产自在天天线| 国产视频首页在线观看| 欧美xxxx性猛交bbbb| 日韩视频在线欧美| 一级黄色大片毛片| 色吧在线观看| 能在线免费看毛片的网站| 观看免费一级毛片| 国产毛片a区久久久久| 久久久久久伊人网av| 成人午夜高清在线视频| 午夜精品在线福利| 直男gayav资源| 欧美精品一区二区大全| 69av精品久久久久久| 六月丁香七月| 男女视频在线观看网站免费| 亚洲国产精品成人综合色| 国产一区二区亚洲精品在线观看| 能在线免费观看的黄片| 夜夜爽夜夜爽视频| 色网站视频免费| 日韩中字成人| 99久久精品国产国产毛片| av专区在线播放| 国产精品综合久久久久久久免费| 欧美精品国产亚洲| 欧美zozozo另类| 好男人视频免费观看在线| 直男gayav资源| 国产精品国产高清国产av| 少妇人妻一区二区三区视频| 日日干狠狠操夜夜爽| 亚洲欧美精品综合久久99| 国产精品国产高清国产av| 久久久久久伊人网av| 亚洲精品国产av成人精品| 黄色欧美视频在线观看| 又粗又爽又猛毛片免费看| 男人舔女人下体高潮全视频| 综合色av麻豆| 国产乱来视频区| 久久精品国产亚洲av涩爱| 亚洲国产精品sss在线观看| 国内少妇人妻偷人精品xxx网站| 国产精品国产三级国产av玫瑰| 久久久色成人| 日韩欧美在线乱码| 黄片wwwwww| 99久久精品国产国产毛片| 国产欧美日韩精品一区二区| 国产乱人视频| 欧美性猛交黑人性爽| 久久精品夜色国产| 国产成人午夜福利电影在线观看| 91aial.com中文字幕在线观看| 最近中文字幕高清免费大全6| 一级av片app| 伦理电影大哥的女人| 日本午夜av视频| 国语对白做爰xxxⅹ性视频网站| 亚洲一区高清亚洲精品| 国产白丝娇喘喷水9色精品| 色综合亚洲欧美另类图片| 乱系列少妇在线播放| 欧美日本视频| 波多野结衣高清无吗| 久久精品久久精品一区二区三区| 久久精品国产亚洲av天美| 汤姆久久久久久久影院中文字幕 | 一边摸一边抽搐一进一小说| 国产乱人偷精品视频| 日本免费一区二区三区高清不卡| 深爱激情五月婷婷| av黄色大香蕉| 国产精品一区二区性色av| 久久久久久大精品| 超碰97精品在线观看| 看十八女毛片水多多多| 亚洲av.av天堂| 欧美zozozo另类| 成人美女网站在线观看视频| 91久久精品电影网| 成人特级av手机在线观看| 又爽又黄无遮挡网站| 成人亚洲欧美一区二区av| 亚洲国产欧美人成| 秋霞在线观看毛片| 搞女人的毛片| 七月丁香在线播放| 午夜激情欧美在线| 成人亚洲欧美一区二区av| 久久99热这里只有精品18| 日韩欧美精品免费久久| 久久精品熟女亚洲av麻豆精品 | 国产美女午夜福利| 波多野结衣巨乳人妻| 国产一区二区在线观看日韩| 老司机影院成人| 久久精品国产亚洲av天美| 国产色爽女视频免费观看| 美女cb高潮喷水在线观看| 国产色婷婷99| 国产真实伦视频高清在线观看| 18禁在线无遮挡免费观看视频| videos熟女内射| 亚洲真实伦在线观看| 欧美成人一区二区免费高清观看| 亚洲综合色惰| 简卡轻食公司| 国产精品乱码一区二三区的特点| 国产精品国产三级国产av玫瑰| 国产美女午夜福利| 好男人在线观看高清免费视频| 午夜爱爱视频在线播放| 亚洲最大成人中文| 久久人人爽人人爽人人片va| 乱系列少妇在线播放| 国语自产精品视频在线第100页| 久久精品影院6| 国产精品久久久久久久久免| 国内精品美女久久久久久| 免费观看人在逋| 黄色欧美视频在线观看| 99久久精品热视频| 少妇人妻精品综合一区二区| 一个人看的www免费观看视频| 久久人人爽人人爽人人片va| 日本三级黄在线观看| 日本五十路高清| 亚洲综合精品二区| 欧美不卡视频在线免费观看| 国产精品女同一区二区软件| 国产又黄又爽又无遮挡在线| 久久热精品热| 性色avwww在线观看| 麻豆一二三区av精品| 哪个播放器可以免费观看大片| 1024手机看黄色片| .国产精品久久| 狂野欧美激情性xxxx在线观看| 高清毛片免费看| 在线观看av片永久免费下载| 亚洲高清免费不卡视频| 欧美高清性xxxxhd video| 日本-黄色视频高清免费观看| 亚洲av一区综合| 丝袜喷水一区| 在线观看av片永久免费下载| 99国产精品一区二区蜜桃av| 成人综合一区亚洲| 大香蕉久久网| 亚洲18禁久久av| 国产91av在线免费观看| 久久欧美精品欧美久久欧美| 精品久久久久久久末码| av专区在线播放| 久久久成人免费电影| 51国产日韩欧美| 九九在线视频观看精品| 国产三级中文精品| 亚洲精品影视一区二区三区av| 日韩欧美在线乱码| АⅤ资源中文在线天堂| 亚洲国产精品久久男人天堂| 亚洲激情五月婷婷啪啪| 性色avwww在线观看| 国产午夜精品论理片| 国产成年人精品一区二区| 国产v大片淫在线免费观看| 在线观看av片永久免费下载| 欧美成人一区二区免费高清观看| 日韩成人伦理影院| 亚洲精品国产成人久久av| 亚洲天堂国产精品一区在线| 中文欧美无线码| 一本一本综合久久| 老司机影院毛片| 99国产精品一区二区蜜桃av| 色5月婷婷丁香| 一级毛片电影观看 | 欧美极品一区二区三区四区| 成年女人看的毛片在线观看| 听说在线观看完整版免费高清| 欧美成人午夜免费资源| 成人毛片a级毛片在线播放| 欧美激情在线99| 99久久精品一区二区三区| 热99在线观看视频| 麻豆av噜噜一区二区三区| 一区二区三区高清视频在线| 亚洲国产最新在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 一边亲一边摸免费视频| 内地一区二区视频在线| 免费看a级黄色片| 一本一本综合久久| 小说图片视频综合网站| 亚洲精品影视一区二区三区av| 日韩国内少妇激情av| 舔av片在线| 成人一区二区视频在线观看| 99热网站在线观看| 桃色一区二区三区在线观看| 国产成人精品久久久久久| 插阴视频在线观看视频| 日本熟妇午夜| 免费大片18禁| 麻豆成人av视频| 舔av片在线| 国产亚洲精品av在线| 亚洲av福利一区| 亚洲不卡免费看| 深夜a级毛片| 一级黄片播放器| av天堂中文字幕网| 久久综合国产亚洲精品| 久久这里有精品视频免费| 中文精品一卡2卡3卡4更新| 婷婷色av中文字幕| 日韩一本色道免费dvd| 亚洲av免费高清在线观看| 夫妻性生交免费视频一级片| 久久精品久久精品一区二区三区| 尾随美女入室| 中国国产av一级| 成人欧美大片| 日本猛色少妇xxxxx猛交久久| 国产日韩欧美在线精品| 国产精品,欧美在线| 又爽又黄a免费视频| 国产精品av视频在线免费观看| 特大巨黑吊av在线直播| 亚洲国产精品sss在线观看| 久久99热6这里只有精品| 神马国产精品三级电影在线观看| av国产久精品久网站免费入址| 黄色一级大片看看| 狠狠狠狠99中文字幕| 国产又色又爽无遮挡免| 亚洲18禁久久av| 联通29元200g的流量卡| 国产av不卡久久| 日韩一区二区视频免费看| 免费电影在线观看免费观看| 亚洲精品456在线播放app| 3wmmmm亚洲av在线观看| 国产v大片淫在线免费观看| 久久国产乱子免费精品| 美女高潮的动态| 国产av在哪里看| 女人被狂操c到高潮| av在线观看视频网站免费| 青春草国产在线视频| 国产精品国产三级专区第一集| 国产美女午夜福利| 国产乱人视频| 人妻少妇偷人精品九色| 永久免费av网站大全| 在线免费十八禁| 两个人的视频大全免费| 午夜免费激情av| 午夜福利在线观看免费完整高清在| 午夜福利视频1000在线观看| 人体艺术视频欧美日本| 日韩欧美精品免费久久| 你懂的网址亚洲精品在线观看 | 午夜精品一区二区三区免费看| 日韩三级伦理在线观看| 欧美97在线视频| 嫩草影院精品99| 国产精品一区二区三区四区免费观看| 亚洲av熟女| 一级黄色大片毛片| 久久久久九九精品影院| 国产午夜福利久久久久久| 国产老妇伦熟女老妇高清| 成人毛片a级毛片在线播放| 午夜激情欧美在线|