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

    胃癌相關(guān)核心基因的生物信息學(xué)分析

    2019-08-06 08:29:02陳秀瓊孟凡橋王雅麗周洋媚唐雯華鄒燕梅
    癌變·畸變·突變 2019年4期
    關(guān)鍵詞:差異基因胃癌通路

    陳秀瓊,孟凡橋,熊 華,王雅麗,周洋媚,唐雯華,鄒燕梅*

    (華中科技大學(xué)同濟(jì)醫(yī)學(xué)院附屬同濟(jì)醫(yī)院腫瘤中心,湖北 武漢 430030)

    胃癌是最常見的癌癥之一,在癌癥相關(guān)死亡原因中排位第2[1-2],大多數(shù)胃癌患者僅在疾病晚期才確診,中位生存期小于1年[3],胃癌確診的金標(biāo)準(zhǔn)是通過內(nèi)鏡或手術(shù)病理活檢,多種血清腫瘤標(biāo)志物如CA199、CA724、CEA、CA125等早已用于胃癌的早期篩查及預(yù)后評(píng)價(jià)[4],但靈敏度及特異度均欠佳。多個(gè)基因靶點(diǎn)也已被挖掘用于胃癌靶向治療藥物的研發(fā),如EGFR、VEGFR、HER2等。但由于早診率低、癌細(xì)胞增殖快、分化程度低、侵襲性強(qiáng),多種檢查及治療手段的聯(lián)合仍未達(dá)到預(yù)想效果。目前,高通量技術(shù)用于數(shù)據(jù)挖掘的運(yùn)用,能有效、快速的篩選出在癌組織與正常組織中顯著差異表達(dá)的基因,并可利用數(shù)據(jù)庫中的臨床資料評(píng)估靶基因的預(yù)后及診斷價(jià)值。因此,在臨床可應(yīng)用于挖掘新的特異性的早期診斷標(biāo)志物,尋找新的治療靶點(diǎn)及預(yù)后評(píng)判生物標(biāo)志物。

    本研究運(yùn)用生物信息分析的方法,從GEO數(shù)據(jù)庫下載3個(gè)不同的胃癌組織芯片表達(dá)譜GSE2685、GSE33335、GSE81948,整理芯片數(shù)據(jù)并從數(shù)據(jù)中獲取在胃癌組織及正常組織中顯著差異表達(dá)的基因,對(duì)差異基因進(jìn)行GO功能注釋及KEGG通路富集分析,篩選出參與胃癌發(fā)生發(fā)展機(jī)制的差異基因。通過構(gòu)建蛋白互作網(wǎng)絡(luò)(protein-proteininteractionnetwork,PPI)選出核心基因,最后對(duì)核心基因進(jìn)行生存分析,找出可能提示胃癌患者預(yù)后、具有診斷與治療意義的靶基因。

    1 材料與方法

    1.1 文件下載及處理

    從GEO數(shù)據(jù)庫(GeneExpressionOmnibus,GEO,http://www.ncbi.nlm.nih.gov/geo/)下載胃癌相關(guān)基因表達(dá)譜芯片(編號(hào)為GSE2685、GSE33335、GSE81948)及其相對(duì)應(yīng)的平臺(tái)信息文件(GPL80、GPL5175、GPL6244),三者包含的樣本量均不低于20例,樣本分類明確,僅包含癌組織與癌旁組織。GSE2685芯片包含8例正常胃組織及22例胃癌組織,GSE33335芯片包含25例正常胃組織及25例胃癌組織,GSE81948芯片包含5例正常胃組織及15例胃癌組織。利用平臺(tái)信息文件將表達(dá)譜芯片中的探針矩陣轉(zhuǎn)換為基因矩陣。

    1.2 差異表達(dá)基因的篩選

    利用R語言軟件(https://www.r-project.org/)的Limma包運(yùn)行已處理的基因矩陣文件,設(shè)置差異表達(dá)基因(differentiallyexpressedgenes,DEGs)篩選標(biāo)準(zhǔn)為P<0.05,差異倍數(shù)(foldchange,F(xiàn)C)對(duì)數(shù)值的絕對(duì)值|Log2(FC)|≥2,分別得到3個(gè)芯片譜中正常組織與胃癌組織間差異表達(dá)的所有基因包括上調(diào)基因與下調(diào)基因,并用火山圖體現(xiàn),用Funrich軟件取三者差異基因的交集,并制作韋恩圖。

    1.3 GO功能注釋和KEGG通路富集分析

    基因本體分析(geneontologyanalysis,GO)是一種常用于解釋基因和基因產(chǎn)物以及識(shí)別高通量基因組或轉(zhuǎn)錄組數(shù)據(jù)特征生物學(xué)屬性的方法[5];京都基因與基因組百科全書(Kyotoencyclopediaofgenesandgenomes,KEGG)是一個(gè)處理基因組、生物途徑、疾病、藥物和化學(xué)物質(zhì)數(shù)據(jù)庫的集合[6]。在基因功能注釋在線網(wǎng)絡(luò)工具(DatabaseforAnnotationVisualizationand IntegratedDiscovery,DAVID)(https://david.ncifcrf.gov/)中,對(duì)上述所得的差異表達(dá)基因進(jìn)行GO功能注釋,包括生物學(xué)過程(biologicalprocess,BP)、細(xì)胞組成(cellular component, CC)、 分 子 功 能 (molecular function,MF)[7],選取校正后P值小于0.05的GO,同時(shí)也進(jìn)行KEGG信號(hào)通路富集分析,用R軟件對(duì)富集分析的結(jié)果進(jìn)行可視化展示。

    1.4 PPI構(gòu)建

    交互基因檢索工具(SearchToolfortheRetrieval ofInteractingGenes,STRING)是一種用于評(píng)價(jià)蛋白質(zhì)相互作用信息的在線工具[8]。將所得的顯著表達(dá)差異基因?qū)隨TRING網(wǎng)絡(luò)在線工具中(http://string-db.org),設(shè)置置信度閾值大于0.4,運(yùn)行STRING,得到差異基因表達(dá)的PPI,隱藏其中未連接點(diǎn)。將STRING中得到的PPI文件導(dǎo)入Cytoscape軟件中,使用CytoHubba網(wǎng)絡(luò)分析插件篩選核心基因,核心基因表達(dá)的蛋白為具有調(diào)節(jié)生理功能的蛋白。

    1.5 核心基因的生存分析

    KM數(shù)據(jù)庫使用10461個(gè)癌癥樣本評(píng)估54675個(gè)基因?qū)ι娴挠绊?,包?143個(gè)乳腺癌,1816個(gè)卵巢癌,2437個(gè)肺癌和1065個(gè)胃癌樣本[9]。利用KM數(shù)據(jù)庫(Kaplan-Meierdatabase,KM)對(duì)核心基因進(jìn)行生存分析,將排名前10的核心基因置于KM數(shù)據(jù)庫中,以檢查其基因表達(dá)水平的改變與患者的5年存活率之間的關(guān)聯(lián)。計(jì)算95%置信區(qū)間(confidence interval,CI)和對(duì)數(shù)秩P值的風(fēng)險(xiǎn)比(hazardratio,HR)并將其可視化。

    1.6 核心基因的診斷價(jià)值及相關(guān)性分析

    整合所有的胃癌樣本及正常組織,篩選出核心基因在所有樣本中的表達(dá)量,將其分為實(shí)驗(yàn)組與對(duì)照組。將兩組數(shù)據(jù)導(dǎo)入GraphPad軟件中,通過繪制ROC曲線(ReceiverOperatingCharacteristicCurve,ROC)將結(jié)果進(jìn)行可視化展示,并計(jì)算出核心基因的ROC曲線下面積(AreaUnderCurve,AUC)。AUC>0.5表示該基因具有良好的診斷價(jià)值。利用GEPIA數(shù)據(jù)庫(Gene ExpressionProfilingInteractiveAnalysis,GEPIA)對(duì)核心基因進(jìn)行皮爾遜相關(guān)性分析,相關(guān)系數(shù)用r表示,r>0.4代表兩者具有較高的相關(guān)性,r的絕對(duì)值越大表示相關(guān)性越強(qiáng),同時(shí)P<0.05表示結(jié)果具有統(tǒng)計(jì)學(xué)意義。

    2 結(jié)果

    2.1 差異表達(dá)基因

    在R軟件的Limma包中設(shè)置差異基因篩選條件為校正后P值<0.05,基因表達(dá)差異倍數(shù)為2,得到DEGs有1839個(gè),上調(diào)基因851個(gè),下調(diào)基因988個(gè),來自GSE2685的上調(diào)基因287個(gè),下調(diào)基因238個(gè),來自GSE33335的上調(diào)基因104個(gè),下調(diào)基因165個(gè),來自GSE81948的上調(diào)基因460個(gè),下調(diào)基因585個(gè),見圖1;取交集后,在3個(gè)表達(dá)譜芯片中均有顯著差異表達(dá)的基因有66個(gè),上調(diào)基因24個(gè),下調(diào)基因42個(gè),見圖2。

    2.2 差異基因的GO功能注釋

    選取校正后P值即FDR<0.05,富集程度最高的前5條GO功能,結(jié)果顯示:差異基因所富集的功能主要表現(xiàn)在細(xì)胞外空間、細(xì)胞外外泌體、消化、細(xì)胞外基質(zhì)組織、膠原纖維組織。如圖3、表1中可見以細(xì)胞外空間與細(xì)胞外外泌體功能區(qū)富集的基因數(shù)目最多。

    2.3 差異基因的KEGG信號(hào)通路富集分析

    DAVID在線工具進(jìn)行差異基因KEGG通路富集分析結(jié)果顯示,按差異基因富集數(shù)目進(jìn)行排序,選取富集程度最高的前9條通路,包括蛋白質(zhì)消化和吸收、胃酸分泌、氮代謝、ECM-受體相互作用、礦物質(zhì)吸收、AGE-RAGE信號(hào)通路等,如圖4。表1中可見差異基因富集程度最高的通路為蛋白質(zhì)消化和吸收與胃酸分泌。

    圖1 差異表達(dá)基因火山圖

    圖2 3個(gè)表達(dá)譜芯片取交集的韋恩圖

    圖3 差異表達(dá)基因的GO功能注釋圖

    圖4 為差異基因的KEGG通路富集分析圖

    2.4 PPI構(gòu)建及核心基因

    將STRING在線數(shù)據(jù)庫構(gòu)建的蛋白互作網(wǎng)絡(luò)PPI數(shù)據(jù)文件導(dǎo)入Cytoscape軟件中,得到的PPI網(wǎng)絡(luò)由48個(gè)基因和186條邊組成,如圖5,圖中節(jié)點(diǎn)為差異基因,節(jié)點(diǎn)之間的邊表示基因之間的聯(lián)系,與節(jié)點(diǎn)相連的邊越多,基因所起的作用越大。CytoHubba網(wǎng)絡(luò)分析插件通過計(jì)算與節(jié)點(diǎn)連接的邊的數(shù)量篩選出核心基因并做出顏色標(biāo)記,顏色越紅,該基因在PPI網(wǎng)絡(luò)中所起的作用越大,由此得出相關(guān)度最高的前10個(gè)核心基因即 TOP2A、COL1A1、COL1A2、ATP4A、COL3A1、COL6A3、FN1、LUM、SPARC、SST,如表2所示,SST和ATP4A在胃癌組織中表達(dá)上調(diào)(較正常組織),TOP2A、COL1A1、COL1A2、COL3A1、COL6A3、FN1、LUM、SPARC在胃癌組織中表達(dá)下調(diào)(較正常組織)。10個(gè)核心基因在胃癌組織中表達(dá)的差異均具有統(tǒng)計(jì)學(xué)意義即FDR均<0.05。

    表1 胃癌組織中差異表達(dá)基因的GO功能注釋和KEGG富集分析

    2.5 核心基因與生存率的關(guān)系

    圖5 PPI網(wǎng)絡(luò)圖

    篩選出的核心基因在胃癌組織與正常組織的差異表達(dá)顯著,為進(jìn)一步挖掘核心基因與胃癌預(yù)后的關(guān)系,我們將核心基因的表達(dá)分為低表達(dá)組和高表達(dá)組,利用KM繪圖儀在線工具繪出每個(gè)核心基因與胃癌患者生存率關(guān)系的KM曲線,并計(jì)算出P值、HR,求得95%CI。在10個(gè)核心基因中,僅COL1A1(1556499_s_at等)、FN1(210495_x_at等)的表達(dá)水平與胃癌患者總生存率的關(guān)系具有顯著的統(tǒng)計(jì)學(xué)意義[COL1A1,HR=1.49,95%CI(1.22, 1.81),P<0.05;FN1,HR=1.46,95%CI(1.23,1.74),P<0.05)],如圖6所示。COL1A1及FN1低表達(dá)組的預(yù)后明顯優(yōu)于高表達(dá)組。

    表2 胃癌基因芯片中10個(gè)核心基因表達(dá)特點(diǎn)

    圖6 生存曲線圖

    2.6 核心基因的診斷價(jià)值

    當(dāng)AUC>0.5時(shí),提示該指標(biāo)具有較好的診斷價(jià)值。利用GraphPad分別計(jì)算核心基因COL1A1與FN1兩者的AUC,得出AUCCOL1A1=0.90,AUCFN1=0.93,如圖7所示,兩者的P值均<0.01。結(jié)果提示,COL1A1和FN1兩者在胃癌患者與正常人群中的差異表達(dá)具有臨床診斷意義,可用于胃癌患者的診斷及預(yù)后的提示。皮爾遜法相關(guān)系數(shù)分析得出COL1A1與FN1兩者之間的r=0.59,P=0.00,結(jié)果表明兩者具有較高的相關(guān)性,在胃癌的發(fā)生發(fā)展中可能具有協(xié)同作用。

    3 討論

    近年來,隨著基因芯片技術(shù)的快速發(fā)展,生物信息學(xué)分析越來越多地用于尋找各種癌癥新的治療靶點(diǎn)和診斷標(biāo)志物,實(shí)現(xiàn)惡性腫瘤的預(yù)后評(píng)價(jià)和精準(zhǔn)診療[10]。本研究通過分析來自GEO數(shù)據(jù)庫的胃癌基因芯片數(shù)據(jù)表達(dá)譜,獲取胃癌組織與正常組織之間的顯著差異表達(dá)基因,包含的851個(gè)上調(diào)基因與988個(gè)下調(diào)基因說明胃癌發(fā)生發(fā)展是由多個(gè)基因表達(dá)異?;蛞种苹蚴Щ畹囊粋€(gè)復(fù)雜的生物學(xué)過程。

    圖7 核心基因的ROC曲線圖

    差異基因的GO功能注釋表明差異基因功能主要表現(xiàn)在細(xì)胞外空間、細(xì)胞外外泌體、消化、細(xì)胞外基質(zhì)組織、膠原纖維組織,KEGG富集分析顯示差異基因主要富集在蛋白質(zhì)消化和吸收、胃酸分泌、氮代謝、ECM-受體相互作用、礦物質(zhì)吸收等。由上可見,差異基因富集的通路及功能主要作用于腫瘤微環(huán)境。近年來,隨著對(duì)腫瘤微環(huán)境研究的深入,大量研究結(jié)果表明,腫瘤微環(huán)境是一個(gè)以免疫抑制機(jī)制為主,免疫激活機(jī)制為輔的位點(diǎn),在腫瘤的免疫治療,治療耐藥,治療敏感性,腫瘤復(fù)發(fā)轉(zhuǎn)移等過程中起著關(guān)鍵作用[11],操縱腫瘤微環(huán)境中不同的信號(hào)通路或基因靶點(diǎn)可以創(chuàng)造比以往任何時(shí)候都更有效的治療。因此,尋找有效的靶點(diǎn)及通路至關(guān)重要。外泌體是由細(xì)胞分泌的小囊泡,被發(fā)現(xiàn)可介導(dǎo)相鄰或遠(yuǎn)端細(xì)胞之間的信號(hào)轉(zhuǎn)導(dǎo)[12-13],而介導(dǎo)腫瘤細(xì)胞與微環(huán)境中基質(zhì)細(xì)胞之間信號(hào)傳導(dǎo)的外泌體被稱為腫瘤細(xì)胞介導(dǎo)的外泌體(tumor-derivedexosomes,TDEs),已有研究表明TDEs通過攜帶不同的miRNA促進(jìn)多種惡性腫瘤的復(fù)發(fā)轉(zhuǎn)移,如外泌體介導(dǎo)的miR-105通過提高血管內(nèi)皮細(xì)胞的通透性促進(jìn)乳腺癌細(xì)胞的遠(yuǎn)處轉(zhuǎn)移,外泌介導(dǎo)的miR-25-3p促進(jìn)結(jié)腸癌細(xì)胞的肝肺轉(zhuǎn)移[14-15],胃癌來源的外泌體通過引起間皮屏障破壞和腹膜纖維化介導(dǎo)腫瘤定向轉(zhuǎn)移[16],如通過傳遞表皮生長(zhǎng)因子(EGFR)調(diào)控肝臟微環(huán)境,促進(jìn)胃癌肝轉(zhuǎn)移[17]。有研究表明胃癌患者血清和組織中細(xì)胞外基質(zhì)蛋白(ECM)通過ECM-受體相互作用誘導(dǎo)ITGB4/FAK/SOX2/HIF-1α信號(hào)通路調(diào)節(jié)胃癌細(xì)胞轉(zhuǎn)移和葡萄糖代謝,對(duì)于預(yù)防腫瘤轉(zhuǎn)移和復(fù)發(fā)的治療靶標(biāo)的發(fā)展具有重要意義[18]。這些差異基因富集的GO功能與KEGG通路為胃癌的起始和進(jìn)展的分子機(jī)制提供了深入的見解,尋找操控胃癌發(fā)生發(fā)展的信號(hào)通路及靶點(diǎn)、開發(fā)新的治療策略能為胃癌患者的生存預(yù)后帶來福音。

    為進(jìn)一步明確差異基因之間的相互作用關(guān)系,我們通過構(gòu)建差異基因的蛋白互作網(wǎng)絡(luò)PPI,運(yùn)用CytoHubba插件篩選出10個(gè)核心基因,在10個(gè)核心基因中,我們發(fā)現(xiàn)COL1A1、FN1與胃癌患者的預(yù)后明顯相關(guān)。COL1A1又稱I型膠原α1鏈,與I型膠原α2鏈(COL1A2)組成典型的I型膠原,I型膠原存在于大多數(shù)結(jié)締組織和胚胎組織中,是膠原家族的重要成員,是細(xì)胞外基質(zhì)的關(guān)鍵結(jié)構(gòu)成分[19-20]。已有研究指出膠原蛋白在胃癌組織中降解水平升高是腫瘤細(xì)胞侵襲周圍組織的關(guān)鍵步驟[21],發(fā)現(xiàn)胃癌組織中COL1A1的上調(diào)與胃癌侵襲性顯著相關(guān),COL1A1的敲除可抑制癌細(xì)胞的增殖,遷移和侵襲能力,同時(shí)也發(fā)現(xiàn)miRNAlet-7i、miR-129-5p表達(dá)上調(diào)通過靶向下調(diào)COL1A1促進(jìn)胃癌侵襲轉(zhuǎn)移[22-23]。生存分析提示,COL1A1的表達(dá)在胃癌組織中顯著降低,低表達(dá)組患者預(yù)后較差(P<0.05),同時(shí),ROC分析顯示AUC=0.90,表明COL1A1在胃癌患者的診斷中具有一定的價(jià)值。因此,我們推測(cè)COL1A1可作為胃癌預(yù)后評(píng)價(jià)、診斷、靶向治療研究對(duì)象的靶標(biāo)之一。FN1,又稱fibronectin1,編碼纖連蛋白,一種以血漿中可溶性二聚體形式存在的糖蛋白,以及細(xì)胞表面和細(xì)胞外基質(zhì)中的二聚體或多聚體形式。纖連蛋白參與細(xì)胞粘附和遷移過程,通過激活MMP2/MMP9通路促進(jìn)胃癌細(xì)胞增殖、遷移和侵襲[24],有研究結(jié)果指出胃癌組織中FN1蛋白的陽性表達(dá)明顯高于正常組織,尤其在淋巴結(jié)轉(zhuǎn)移陽性、分期較晚的胃癌患者中,同時(shí)也發(fā)現(xiàn)miRNA-200c表達(dá)上調(diào)通過靶向下調(diào)FN1促進(jìn)胃癌細(xì)胞的增值、侵襲、轉(zhuǎn)移[25-26]。AUC=0.93提示FN1對(duì)胃癌患者具有良好的診斷價(jià)值,通過生存分析顯示,與正常胃組織相比,胃癌組織中的FN1明顯降低,F(xiàn)N1低表達(dá)組預(yù)后差(P<0.05),Cai等[27]的研究表明,在結(jié)腸癌(CRC)中,低表達(dá)的FN1通過與ITGA5相互作用促進(jìn)細(xì)胞凋亡,抑制CRC細(xì)胞的存活、侵襲和遷移。而FN1在胃癌中的抗腫瘤作用仍需進(jìn)一步研究加以驗(yàn)證。皮爾遜相關(guān)系數(shù)R=0.59,表明COL1A1與FN1在胃癌的發(fā)生發(fā)展過程中具有協(xié)同作用。目前尚無相關(guān)基礎(chǔ)研究對(duì)兩者的協(xié)同作用深入探討,因此,設(shè)計(jì)進(jìn)一步的實(shí)驗(yàn)進(jìn)行研究可能帶來新的突破。

    綜上,生物信息分析方法可為未來胃癌基因組個(gè)體化診斷和治療提供有力證據(jù),利用基因芯片表達(dá)譜數(shù)據(jù)分析得到的核心基因所富集的功能與通路說明胃癌的發(fā)生發(fā)展是一個(gè)多基因突變、表達(dá)異常導(dǎo)致細(xì)胞無限增殖的復(fù)雜過程。核心基因COL1A1、FN1在胃癌組織中的表達(dá)與患者預(yù)后明顯相關(guān),未來我們將進(jìn)行進(jìn)一步的基礎(chǔ)實(shí)驗(yàn)及臨床研究證實(shí)其作為評(píng)判預(yù)后、分子靶向治療靶標(biāo)的價(jià)值。

    猜你喜歡
    差異基因胃癌通路
    ICR鼠肝和腎毒性損傷生物標(biāo)志物的篩選
    基于RNA 測(cè)序研究人參二醇對(duì)大鼠心血管內(nèi)皮細(xì)胞基因表達(dá)的影響 (正文見第26 頁)
    P53及Ki67在胃癌中的表達(dá)及其臨床意義
    Kisspeptin/GPR54信號(hào)通路促使性早熟形成的作用觀察
    胃癌組織中LKB1和VEGF-C的表達(dá)及其意義
    proBDNF-p75NTR通路抑制C6細(xì)胞增殖
    胃癌組織中VEGF和ILK的表達(dá)及意義
    通路快建林翰:對(duì)重模式應(yīng)有再認(rèn)識(shí)
    SSH技術(shù)在絲狀真菌功能基因篩選中的應(yīng)用
    中醫(yī)辨證結(jié)合化療治療中晚期胃癌50例
    免费不卡的大黄色大毛片视频在线观看 | 韩国av在线不卡| 熟女人妻精品中文字幕| 精品久久久久久成人av| 看片在线看免费视频| 日本熟妇午夜| 中文字幕av成人在线电影| 亚洲欧美成人综合另类久久久 | 久久欧美精品欧美久久欧美| 亚洲欧洲国产日韩| 成人特级av手机在线观看| 天天躁日日操中文字幕| 99热这里只有是精品在线观看| 国产亚洲av片在线观看秒播厂 | 国产 一区精品| 亚洲综合色惰| 免费搜索国产男女视频| 久久人人爽人人片av| 免费黄网站久久成人精品| av专区在线播放| 岛国毛片在线播放| 日本与韩国留学比较| 午夜福利在线在线| 深夜精品福利| videossex国产| 精品人妻熟女av久视频| 亚洲在线自拍视频| 国语自产精品视频在线第100页| 午夜福利在线观看免费完整高清在 | 免费av毛片视频| 偷拍熟女少妇极品色| 日本一二三区视频观看| 欧洲精品卡2卡3卡4卡5卡区| 在线播放无遮挡| 国语自产精品视频在线第100页| 99国产极品粉嫩在线观看| 国产精品免费一区二区三区在线| 天天躁日日操中文字幕| 级片在线观看| 天堂中文最新版在线下载 | 少妇熟女欧美另类| 成人午夜精彩视频在线观看| 天天躁夜夜躁狠狠久久av| 变态另类丝袜制服| 日本欧美国产在线视频| 亚洲人成网站在线播| 在线国产一区二区在线| 成人毛片60女人毛片免费| 欧美bdsm另类| 午夜亚洲福利在线播放| 麻豆精品久久久久久蜜桃| 国产高清三级在线| 在线观看免费视频日本深夜| 99久久精品热视频| 两个人视频免费观看高清| 精品人妻熟女av久视频| 99热精品在线国产| 国产一级毛片七仙女欲春2| 最新中文字幕久久久久| 99国产精品一区二区蜜桃av| 99热这里只有是精品在线观看| 啦啦啦观看免费观看视频高清| 亚洲色图av天堂| 欧美日本亚洲视频在线播放| 97在线视频观看| videossex国产| 午夜激情福利司机影院| 欧美最黄视频在线播放免费| 麻豆一二三区av精品| 在线播放国产精品三级| 一本精品99久久精品77| 搞女人的毛片| ponron亚洲| 国产精品,欧美在线| 91久久精品国产一区二区成人| 美女 人体艺术 gogo| 亚洲在线自拍视频| av.在线天堂| 欧美bdsm另类| 69av精品久久久久久| 国产69精品久久久久777片| 午夜久久久久精精品| 欧美变态另类bdsm刘玥| 三级国产精品欧美在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av在线天堂中文字幕| 成年免费大片在线观看| 成人欧美大片| 国产爱豆传媒在线观看| 天堂√8在线中文| 97热精品久久久久久| 精品日产1卡2卡| 亚洲成av人片在线播放无| 国产成人精品久久久久久| 国产伦精品一区二区三区视频9| 97在线视频观看| av免费观看日本| 日本撒尿小便嘘嘘汇集6| 欧美精品国产亚洲| 男人舔女人下体高潮全视频| 亚洲成人中文字幕在线播放| 日本色播在线视频| 两个人的视频大全免费| 亚洲欧美精品自产自拍| 搞女人的毛片| 久久鲁丝午夜福利片| h日本视频在线播放| 欧美区成人在线视频| 国产精品1区2区在线观看.| av福利片在线观看| 啦啦啦观看免费观看视频高清| 国产色婷婷99| 97在线视频观看| 亚洲精品成人久久久久久| 日韩av不卡免费在线播放| 亚洲欧美精品专区久久| 最新中文字幕久久久久| 日韩制服骚丝袜av| 国产成人福利小说| 成人无遮挡网站| 热99re8久久精品国产| 国产欧美日韩精品一区二区| 国产成人精品久久久久久| 国产高清激情床上av| 乱人视频在线观看| 麻豆成人午夜福利视频| 麻豆av噜噜一区二区三区| 99热这里只有精品一区| 亚洲,欧美,日韩| 12—13女人毛片做爰片一| 菩萨蛮人人尽说江南好唐韦庄 | 精品一区二区三区人妻视频| 男人舔女人下体高潮全视频| 美女xxoo啪啪120秒动态图| 亚洲av免费高清在线观看| 亚洲美女搞黄在线观看| 欧美高清成人免费视频www| 婷婷亚洲欧美| 91麻豆精品激情在线观看国产| 99久国产av精品国产电影| 美女被艹到高潮喷水动态| 免费看a级黄色片| 成人午夜高清在线视频| 大又大粗又爽又黄少妇毛片口| 男女啪啪激烈高潮av片| 日韩强制内射视频| av黄色大香蕉| 亚洲真实伦在线观看| 欧美zozozo另类| 国产综合懂色| 自拍偷自拍亚洲精品老妇| 少妇猛男粗大的猛烈进出视频 | 偷拍熟女少妇极品色| 日本三级黄在线观看| 日本黄色视频三级网站网址| 色哟哟·www| a级毛色黄片| 看免费成人av毛片| 美女黄网站色视频| 成人鲁丝片一二三区免费| 大香蕉久久网| av又黄又爽大尺度在线免费看 | 亚洲av电影不卡..在线观看| 两个人的视频大全免费| 一个人观看的视频www高清免费观看| 国产又黄又爽又无遮挡在线| 少妇人妻一区二区三区视频| 美女cb高潮喷水在线观看| 国产激情偷乱视频一区二区| 草草在线视频免费看| 国产精品久久久久久精品电影小说 | 日本熟妇午夜| 亚洲欧洲日产国产| 中文字幕制服av| 全区人妻精品视频| or卡值多少钱| 日本一本二区三区精品| 精品人妻视频免费看| 欧美xxxx黑人xx丫x性爽| 国产 一区 欧美 日韩| 一个人看视频在线观看www免费| 99热网站在线观看| 久久精品久久久久久久性| av天堂中文字幕网| 国产精品精品国产色婷婷| 又爽又黄无遮挡网站| 国产黄片美女视频| 少妇的逼好多水| 91精品国产九色| 免费观看在线日韩| 免费人成视频x8x8入口观看| 日韩成人伦理影院| 亚洲一级一片aⅴ在线观看| 国产麻豆成人av免费视频| 欧美激情久久久久久爽电影| 日韩强制内射视频| 国产日韩欧美在线精品| 午夜福利高清视频| 18禁在线无遮挡免费观看视频| 精品国内亚洲2022精品成人| 亚洲综合色惰| 亚洲欧美精品综合久久99| av.在线天堂| 一级毛片久久久久久久久女| 一级毛片aaaaaa免费看小| 亚洲精品成人久久久久久| 久久午夜亚洲精品久久| 黄色配什么色好看| 国产一区二区三区在线臀色熟女| 国产色婷婷99| 久久久久久久久大av| 国产成人aa在线观看| 国产精品国产高清国产av| 国产亚洲91精品色在线| 亚洲不卡免费看| 欧美+亚洲+日韩+国产| 欧美在线一区亚洲| 日日摸夜夜添夜夜添av毛片| 欧美精品国产亚洲| 免费黄网站久久成人精品| 日韩av不卡免费在线播放| 性插视频无遮挡在线免费观看| 99热精品在线国产| 久久中文看片网| 干丝袜人妻中文字幕| 你懂的网址亚洲精品在线观看 | 又爽又黄无遮挡网站| 不卡一级毛片| 日韩欧美精品免费久久| 日韩,欧美,国产一区二区三区 | 亚洲国产精品成人久久小说 | 国产私拍福利视频在线观看| 人体艺术视频欧美日本| 中文欧美无线码| 欧美性猛交黑人性爽| 欧美丝袜亚洲另类| 亚洲av不卡在线观看| 中文字幕久久专区| 亚洲欧美成人精品一区二区| 在线观看免费视频日本深夜| 国产一区二区在线观看日韩| 欧美zozozo另类| 高清毛片免费看| 亚洲成人av在线免费| www.色视频.com| 亚洲五月天丁香| 精品久久久久久久久久久久久| 人人妻人人澡人人爽人人夜夜 | 99riav亚洲国产免费| 亚洲美女搞黄在线观看| 免费av不卡在线播放| eeuss影院久久| 嫩草影院入口| 日本一本二区三区精品| 夜夜看夜夜爽夜夜摸| 黄色视频,在线免费观看| 又爽又黄a免费视频| 中文欧美无线码| 国产淫片久久久久久久久| 看十八女毛片水多多多| 欧美日韩乱码在线| 亚洲国产日韩欧美精品在线观看| 哪里可以看免费的av片| 亚洲精品乱码久久久v下载方式| ponron亚洲| 午夜老司机福利剧场| 亚洲,欧美,日韩| 亚洲欧洲日产国产| 午夜久久久久精精品| 91av网一区二区| 国产真实乱freesex| 夫妻性生交免费视频一级片| av视频在线观看入口| 午夜视频国产福利| 免费电影在线观看免费观看| 亚洲国产精品国产精品| 精品国产三级普通话版| 欧美成人精品欧美一级黄| 日韩人妻高清精品专区| 日韩高清综合在线| 亚洲经典国产精华液单| 成人午夜高清在线视频| 久久久久免费精品人妻一区二区| 97在线视频观看| 国产在线精品亚洲第一网站| 丰满的人妻完整版| 蜜桃亚洲精品一区二区三区| 免费看a级黄色片| 亚洲欧美中文字幕日韩二区| 国产午夜福利久久久久久| 嫩草影院入口| 国语自产精品视频在线第100页| 在线观看美女被高潮喷水网站| 国产伦精品一区二区三区四那| 色综合亚洲欧美另类图片| 日本撒尿小便嘘嘘汇集6| 欧美精品一区二区大全| 国产精品嫩草影院av在线观看| av福利片在线观看| 18禁在线无遮挡免费观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 乱人视频在线观看| 欧美日本视频| av在线天堂中文字幕| 国产激情偷乱视频一区二区| 亚洲精品粉嫩美女一区| 亚洲丝袜综合中文字幕| 有码 亚洲区| 国产成人精品一,二区 | 91狼人影院| 色哟哟·www| 久久久精品大字幕| 中文精品一卡2卡3卡4更新| 天堂影院成人在线观看| 免费av毛片视频| 国产91av在线免费观看| 午夜a级毛片| 国产精品伦人一区二区| 哪个播放器可以免费观看大片| 在现免费观看毛片| 欧美日韩一区二区视频在线观看视频在线 | 99久国产av精品| 国产乱人视频| 在线免费观看不下载黄p国产| 1000部很黄的大片| 在现免费观看毛片| 悠悠久久av| 毛片女人毛片| 欧美三级亚洲精品| 99riav亚洲国产免费| 国产精品1区2区在线观看.| 成人午夜高清在线视频| 伦理电影大哥的女人| 日本av手机在线免费观看| 日韩欧美在线乱码| 亚洲精品粉嫩美女一区| 国产在线男女| 黄片wwwwww| 国产淫片久久久久久久久| 国产午夜福利久久久久久| 高清毛片免费观看视频网站| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久av不卡| 国产白丝娇喘喷水9色精品| 日本爱情动作片www.在线观看| 亚洲成av人片在线播放无| 丰满乱子伦码专区| 亚洲成av人片在线播放无| 久久久久久大精品| 免费av观看视频| 久久久久久大精品| 中文资源天堂在线| 深爱激情五月婷婷| 成人永久免费在线观看视频| 观看美女的网站| 国产成人一区二区在线| 最近中文字幕高清免费大全6| 久久久久久久久久成人| 综合色丁香网| 少妇高潮的动态图| 久久韩国三级中文字幕| www.av在线官网国产| av视频在线观看入口| 久久鲁丝午夜福利片| 岛国在线免费视频观看| 久久久久久久久久成人| 日日干狠狠操夜夜爽| 久久精品国产亚洲av涩爱 | 波野结衣二区三区在线| 亚洲精品久久国产高清桃花| 久久精品91蜜桃| 男插女下体视频免费在线播放| 成人亚洲欧美一区二区av| 国模一区二区三区四区视频| 你懂的网址亚洲精品在线观看 | 中文字幕久久专区| 一级毛片电影观看 | 伦理电影大哥的女人| 卡戴珊不雅视频在线播放| 亚洲国产精品合色在线| 毛片一级片免费看久久久久| av福利片在线观看| 日本欧美国产在线视频| 日韩一本色道免费dvd| 国产黄色视频一区二区在线观看 | a级毛片a级免费在线| 免费一级毛片在线播放高清视频| 激情 狠狠 欧美| 色5月婷婷丁香| 天美传媒精品一区二区| 国产精品一区二区三区四区免费观看| 毛片一级片免费看久久久久| 亚洲国产高清在线一区二区三| 亚洲av第一区精品v没综合| 级片在线观看| 久久人人精品亚洲av| av卡一久久| av在线老鸭窝| 国产精品福利在线免费观看| 久久久国产成人免费| 人体艺术视频欧美日本| 亚洲最大成人av| 尤物成人国产欧美一区二区三区| 国产三级在线视频| 亚洲av一区综合| 国产大屁股一区二区在线视频| 国产爱豆传媒在线观看| 亚洲aⅴ乱码一区二区在线播放| 中文字幕人妻熟人妻熟丝袜美| 两个人的视频大全免费| 亚洲精品影视一区二区三区av| 欧美+日韩+精品| 天堂√8在线中文| 欧美高清性xxxxhd video| ponron亚洲| 春色校园在线视频观看| 国产伦在线观看视频一区| 国产精品免费一区二区三区在线| 成熟少妇高潮喷水视频| av在线蜜桃| 亚洲精品影视一区二区三区av| 亚洲精华国产精华液的使用体验 | 国产老妇伦熟女老妇高清| 久久鲁丝午夜福利片| 91久久精品国产一区二区三区| 内地一区二区视频在线| 亚洲激情五月婷婷啪啪| 在线观看午夜福利视频| 久久久久久久久久成人| 国产高清视频在线观看网站| a级毛色黄片| 亚洲中文字幕日韩| 免费观看人在逋| 一级毛片aaaaaa免费看小| 国产伦精品一区二区三区视频9| 麻豆精品久久久久久蜜桃| 国产极品精品免费视频能看的| 亚洲av一区综合| 可以在线观看毛片的网站| 精品国产三级普通话版| 久久99热6这里只有精品| 亚洲精品亚洲一区二区| 简卡轻食公司| 成人毛片a级毛片在线播放| 国产亚洲精品久久久com| 在线观看66精品国产| 亚洲五月天丁香| 亚洲,欧美,日韩| 国产精品日韩av在线免费观看| 别揉我奶头 嗯啊视频| 亚洲成人精品中文字幕电影| 国产精品嫩草影院av在线观看| 免费无遮挡裸体视频| 成人漫画全彩无遮挡| 联通29元200g的流量卡| 99热网站在线观看| 欧美成人精品欧美一级黄| 欧美又色又爽又黄视频| 亚洲精品国产av成人精品| 白带黄色成豆腐渣| 国产精品女同一区二区软件| 日韩大尺度精品在线看网址| 国内揄拍国产精品人妻在线| 成人亚洲精品av一区二区| 免费观看a级毛片全部| 如何舔出高潮| 免费不卡的大黄色大毛片视频在线观看 | 男插女下体视频免费在线播放| 一区二区三区高清视频在线| 夜夜夜夜夜久久久久| 亚洲欧美成人综合另类久久久 | 一区二区三区高清视频在线| 成人二区视频| 亚洲一级一片aⅴ在线观看| 国产在线男女| 欧美在线一区亚洲| 两个人的视频大全免费| 一级av片app| 丰满的人妻完整版| 赤兔流量卡办理| 18禁在线播放成人免费| 色综合色国产| 老司机影院成人| 日本五十路高清| 精品久久久久久成人av| 少妇的逼好多水| 国产午夜精品论理片| 69av精品久久久久久| 亚洲中文字幕一区二区三区有码在线看| 国产美女午夜福利| 国产不卡一卡二| 亚洲在线观看片| 欧美一区二区亚洲| 一级黄色大片毛片| 在线免费观看不下载黄p国产| 久久精品国产自在天天线| 亚洲在线观看片| 一级毛片久久久久久久久女| 在线天堂最新版资源| 中文字幕精品亚洲无线码一区| 少妇丰满av| 色尼玛亚洲综合影院| 五月玫瑰六月丁香| av在线亚洲专区| 一级毛片我不卡| 男女啪啪激烈高潮av片| 亚洲美女搞黄在线观看| 国产精品日韩av在线免费观看| 一个人看视频在线观看www免费| 高清毛片免费观看视频网站| 成人午夜高清在线视频| 精品一区二区三区人妻视频| 中出人妻视频一区二区| 亚洲熟妇中文字幕五十中出| 搞女人的毛片| 久久久色成人| 五月伊人婷婷丁香| 亚洲欧美日韩高清在线视频| av在线天堂中文字幕| 女同久久另类99精品国产91| 日韩一区二区三区影片| 91麻豆精品激情在线观看国产| 欧美高清成人免费视频www| 美女内射精品一级片tv| 亚洲av中文av极速乱| 色尼玛亚洲综合影院| a级毛片免费高清观看在线播放| 网址你懂的国产日韩在线| 国产精品一区www在线观看| 欧美bdsm另类| 免费观看精品视频网站| 欧美zozozo另类| 欧美一区二区精品小视频在线| 中文字幕制服av| 床上黄色一级片| 少妇高潮的动态图| 日韩 亚洲 欧美在线| 久久这里只有精品中国| 欧美高清成人免费视频www| 国模一区二区三区四区视频| av在线观看视频网站免费| 中文在线观看免费www的网站| 床上黄色一级片| 99久久九九国产精品国产免费| 精品一区二区三区视频在线| 男人舔女人下体高潮全视频| 午夜视频国产福利| 97热精品久久久久久| 免费看a级黄色片| 十八禁国产超污无遮挡网站| 久久综合国产亚洲精品| 亚洲中文字幕一区二区三区有码在线看| 国产精品伦人一区二区| 哪里可以看免费的av片| 看十八女毛片水多多多| 99热6这里只有精品| av免费在线看不卡| 搞女人的毛片| 亚洲最大成人手机在线| 男的添女的下面高潮视频| 欧美3d第一页| 亚洲欧洲国产日韩| 乱系列少妇在线播放| 日本三级黄在线观看| www.色视频.com| 国产大屁股一区二区在线视频| 精品午夜福利在线看| 午夜视频国产福利| 国产精品日韩av在线免费观看| 国产激情偷乱视频一区二区| 亚洲av熟女| 女同久久另类99精品国产91| 国产精品久久久久久精品电影小说 | 99在线视频只有这里精品首页| 国产亚洲5aaaaa淫片| 国产v大片淫在线免费观看| 国产精品久久久久久精品电影| 联通29元200g的流量卡| 舔av片在线| 国产精品乱码一区二三区的特点| 男女视频在线观看网站免费| 三级男女做爰猛烈吃奶摸视频| 美女内射精品一级片tv| 欧美日韩一区二区视频在线观看视频在线 | 国产成人精品婷婷| 在线播放无遮挡| 男女视频在线观看网站免费| 边亲边吃奶的免费视频| 啦啦啦啦在线视频资源| 久久亚洲精品不卡| 国产激情偷乱视频一区二区| 日日干狠狠操夜夜爽| av.在线天堂| 久久精品国产亚洲网站| av国产免费在线观看| 女的被弄到高潮叫床怎么办| 久久久久久久久久黄片| 97超碰精品成人国产| 国产成人a∨麻豆精品| 12—13女人毛片做爰片一| 干丝袜人妻中文字幕| 丰满人妻一区二区三区视频av| 国产v大片淫在线免费观看| 五月伊人婷婷丁香| 男女那种视频在线观看| 美女xxoo啪啪120秒动态图| 国产精品一区二区在线观看99 | 成人毛片a级毛片在线播放| 欧美3d第一页| 精品一区二区免费观看| 亚洲人成网站在线观看播放| 一区二区三区四区激情视频 | 国产真实乱freesex| 亚洲成人中文字幕在线播放| 久久久久久久亚洲中文字幕|