• <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.在线天堂| 国产高清不卡午夜福利| 免费大片黄手机在线观看| 22中文网久久字幕| 热99国产精品久久久久久7| 午夜福利网站1000一区二区三区| 国产精品偷伦视频观看了| 五月玫瑰六月丁香| 校园人妻丝袜中文字幕| 简卡轻食公司| 久久人人爽人人片av| 中文字幕免费在线视频6| 91精品一卡2卡3卡4卡| 99精国产麻豆久久婷婷| 欧美日韩精品成人综合77777| 大片电影免费在线观看免费| 美女福利国产在线| 精品国产露脸久久av麻豆| 男女啪啪激烈高潮av片| 免费av中文字幕在线| 久久久久久久精品精品| 亚洲精品中文字幕在线视频| 青春草视频在线免费观看| 精品人妻熟女毛片av久久网站| 波野结衣二区三区在线| 久久99蜜桃精品久久| 欧美日韩成人在线一区二区| 男女边摸边吃奶| 看十八女毛片水多多多| kizo精华| 欧美3d第一页| 一本久久精品| 中文精品一卡2卡3卡4更新| 三级国产精品片| 美女cb高潮喷水在线观看| 人妻人人澡人人爽人人| 亚洲欧美日韩卡通动漫| 国产国语露脸激情在线看| 一区二区三区精品91| 永久网站在线| 日本黄色日本黄色录像| 老女人水多毛片| 性色avwww在线观看| 九色成人免费人妻av| 亚洲一区二区三区欧美精品| 成人毛片a级毛片在线播放| 看十八女毛片水多多多| 久久久久久久久久久免费av| 黄色配什么色好看| 日韩成人伦理影院| 欧美激情极品国产一区二区三区 | 观看美女的网站| 久久久亚洲精品成人影院| 精品国产露脸久久av麻豆| 免费高清在线观看日韩| 黄色配什么色好看| 亚洲精品一二三| 熟女av电影| 欧美激情国产日韩精品一区| 免费看光身美女| 欧美日本中文国产一区发布| 亚洲精华国产精华液的使用体验| 少妇精品久久久久久久| 中国美白少妇内射xxxbb| 亚洲av日韩在线播放| 人人妻人人澡人人看| 精品一区在线观看国产| 久久国产亚洲av麻豆专区| 亚洲精品成人av观看孕妇| 亚洲国产精品国产精品| 亚洲丝袜综合中文字幕| 街头女战士在线观看网站| 中国三级夫妇交换| 亚洲精品国产av成人精品| 国产精品蜜桃在线观看| 国产又色又爽无遮挡免| 一级,二级,三级黄色视频| 97在线视频观看| 我要看黄色一级片免费的| 寂寞人妻少妇视频99o| 美女脱内裤让男人舔精品视频| 亚洲av成人精品一二三区| 校园人妻丝袜中文字幕| 如何舔出高潮| 精品久久久噜噜| 永久免费av网站大全| 日韩制服骚丝袜av| 美女内射精品一级片tv| 如何舔出高潮| 大香蕉久久网| 亚洲美女视频黄频| av在线播放精品| 一本一本综合久久| av在线app专区| 日本黄色日本黄色录像| 久久99精品国语久久久| 精品一区在线观看国产| 欧美成人精品欧美一级黄| 秋霞伦理黄片| 人成视频在线观看免费观看| 美女cb高潮喷水在线观看| 成年av动漫网址| 黑丝袜美女国产一区| 久久精品国产自在天天线| 人人澡人人妻人| 超碰97精品在线观看| 国产淫语在线视频| 两个人的视频大全免费| 国产精品三级大全| 伦理电影免费视频| 亚洲精品日韩av片在线观看| 美女福利国产在线| 秋霞伦理黄片| 亚洲精品亚洲一区二区| 综合色丁香网| 飞空精品影院首页| 精品酒店卫生间| 视频在线观看一区二区三区| 黄色毛片三级朝国网站| 亚洲图色成人| av网站免费在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区乱码不卡18| 18禁观看日本| 波野结衣二区三区在线| 亚洲av男天堂| 久久 成人 亚洲| 999精品在线视频| 另类亚洲欧美激情| 国产精品偷伦视频观看了| 中文字幕最新亚洲高清| 另类精品久久| 国产白丝娇喘喷水9色精品| 久久99一区二区三区| a级毛片黄视频| 91成人精品电影| 天堂中文最新版在线下载| 国产精品人妻久久久影院| 最后的刺客免费高清国语| 99久久人妻综合| 国产白丝娇喘喷水9色精品| 26uuu在线亚洲综合色| 丰满迷人的少妇在线观看| 青青草视频在线视频观看| 国产黄色免费在线视频| 日日啪夜夜爽| av线在线观看网站| 夜夜爽夜夜爽视频| 日韩电影二区| 精品亚洲成国产av| 人妻少妇偷人精品九色| 欧美性感艳星| 少妇猛男粗大的猛烈进出视频| 能在线免费看毛片的网站| 亚洲国产精品成人久久小说| 婷婷成人精品国产| 丁香六月天网| 欧美人与善性xxx| 亚洲人成77777在线视频| 满18在线观看网站| 免费观看a级毛片全部| 天堂8中文在线网| 国产黄色视频一区二区在线观看| 精品久久久噜噜| 亚洲国产精品一区二区三区在线| 日韩av免费高清视频| 免费日韩欧美在线观看| 久久av网站| 中文欧美无线码| 久久99热这里只频精品6学生| 婷婷成人精品国产| 人体艺术视频欧美日本| √禁漫天堂资源中文www| 亚洲精品456在线播放app| 99re6热这里在线精品视频| 国产精品一国产av| 最后的刺客免费高清国语| 免费黄网站久久成人精品| 亚洲精品国产av蜜桃| 性高湖久久久久久久久免费观看| 99久久中文字幕三级久久日本| 中文字幕av电影在线播放| 久久久久久久大尺度免费视频| 人体艺术视频欧美日本| √禁漫天堂资源中文www| 午夜激情福利司机影院| 国产在线视频一区二区| 各种免费的搞黄视频| 国产精品久久久久成人av| 国产精品麻豆人妻色哟哟久久| 成年人免费黄色播放视频| 26uuu在线亚洲综合色| 嫩草影院入口| 国产色爽女视频免费观看| 一区二区三区精品91| 18禁观看日本| 国产精品不卡视频一区二区| 亚洲怡红院男人天堂| 午夜老司机福利剧场| 亚洲国产毛片av蜜桃av| 纵有疾风起免费观看全集完整版| 大码成人一级视频| 婷婷成人精品国产| 精品国产露脸久久av麻豆| 国产男人的电影天堂91| 免费看av在线观看网站| 日韩成人伦理影院| 成人国产av品久久久| 久久精品久久久久久久性| 日本黄色日本黄色录像| 日本vs欧美在线观看视频| 精品久久蜜臀av无| 久久影院123| 国产老妇伦熟女老妇高清| 亚洲av在线观看美女高潮| 在线播放无遮挡| 午夜av观看不卡| 插逼视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产成人aa在线观看| 两个人免费观看高清视频| 午夜激情久久久久久久| 在线 av 中文字幕| 国产色婷婷99| 日本色播在线视频| 性色avwww在线观看| 日本欧美视频一区| 亚洲精品日韩在线中文字幕| 亚洲,一卡二卡三卡| 成人综合一区亚洲| 只有这里有精品99| 人人妻人人添人人爽欧美一区卜| 国语对白做爰xxxⅹ性视频网站| 国产成人a∨麻豆精品| 黑人巨大精品欧美一区二区蜜桃 | 22中文网久久字幕| 国产国语露脸激情在线看| 久久久久人妻精品一区果冻| 插阴视频在线观看视频| 日韩亚洲欧美综合| 99热这里只有精品一区| 亚洲丝袜综合中文字幕| 你懂的网址亚洲精品在线观看| av免费观看日本| 亚洲欧洲日产国产| 人妻系列 视频| 黄色配什么色好看| 国产亚洲欧美精品永久| 人妻系列 视频| 亚洲经典国产精华液单| 久久青草综合色| 男人爽女人下面视频在线观看| 午夜免费男女啪啪视频观看| 人人妻人人澡人人看| a级毛色黄片| 亚洲av二区三区四区| 亚洲精品色激情综合| 少妇熟女欧美另类| 人妻制服诱惑在线中文字幕| 七月丁香在线播放| 午夜激情av网站| 亚洲,欧美,日韩| 一级毛片 在线播放| 国产精品99久久99久久久不卡 | 国产免费福利视频在线观看| 黄色怎么调成土黄色| 国产精品国产三级国产专区5o| 成人亚洲精品一区在线观看| 国产精品免费大片| 久久精品国产鲁丝片午夜精品| 国产免费一级a男人的天堂| 午夜久久久在线观看| 天堂俺去俺来也www色官网| 欧美+日韩+精品| 国产精品人妻久久久影院| 啦啦啦啦在线视频资源| 满18在线观看网站| 国产免费视频播放在线视频| 多毛熟女@视频| 中文字幕最新亚洲高清| 秋霞伦理黄片| 99热国产这里只有精品6| 少妇熟女欧美另类| 亚洲欧美一区二区三区黑人 | av播播在线观看一区| 国产免费一区二区三区四区乱码| 午夜免费男女啪啪视频观看| 搡女人真爽免费视频火全软件| 老司机影院成人| 欧美日本中文国产一区发布| 精品久久久久久久久av| 丝瓜视频免费看黄片| 伊人亚洲综合成人网| 亚洲精品av麻豆狂野| 国产高清三级在线| 一区二区三区精品91| 亚洲高清免费不卡视频| av线在线观看网站| 美女主播在线视频| 黄色欧美视频在线观看| 视频在线观看一区二区三区| 亚洲国产精品成人久久小说| 乱码一卡2卡4卡精品| 老熟女久久久| 亚洲四区av| 亚洲精品成人av观看孕妇| 欧美精品高潮呻吟av久久| 国产免费现黄频在线看| 亚洲国产av新网站| 日日爽夜夜爽网站| 久久久久久人妻| 18+在线观看网站| 日本黄大片高清| 大片电影免费在线观看免费| 国产精品秋霞免费鲁丝片| 日日摸夜夜添夜夜爱| 亚洲国产欧美日韩在线播放| 夫妻性生交免费视频一级片| 国产国拍精品亚洲av在线观看| 校园人妻丝袜中文字幕| 亚洲人与动物交配视频| 91精品一卡2卡3卡4卡| 午夜av观看不卡| 青青草视频在线视频观看| 晚上一个人看的免费电影| 在线播放无遮挡| .国产精品久久| √禁漫天堂资源中文www| 一区二区三区免费毛片| 丝瓜视频免费看黄片| 黑丝袜美女国产一区| 国产一区二区在线观看av| 一本一本综合久久| 赤兔流量卡办理| 久久99蜜桃精品久久| a级毛片免费高清观看在线播放| a级毛色黄片| 欧美日本中文国产一区发布| 国产一区二区在线观看日韩| 欧美变态另类bdsm刘玥| 波野结衣二区三区在线| 少妇丰满av| 亚洲精品av麻豆狂野| 最近的中文字幕免费完整| 久久这里有精品视频免费| 日本免费在线观看一区| 蜜臀久久99精品久久宅男| 欧美另类一区| 成人二区视频| .国产精品久久| 成人18禁高潮啪啪吃奶动态图 | 91精品伊人久久大香线蕉| 中文字幕人妻熟人妻熟丝袜美| 男女边摸边吃奶| 天美传媒精品一区二区| 十八禁高潮呻吟视频| freevideosex欧美| 一级毛片我不卡| 成人综合一区亚洲| 男人操女人黄网站| 能在线免费看毛片的网站| 如何舔出高潮| 国产成人freesex在线| 欧美日韩在线观看h| 国产免费视频播放在线视频| 久久久久久久久久久丰满| 在线天堂最新版资源| 亚洲美女视频黄频| 亚洲国产色片| 欧美97在线视频| 插逼视频在线观看| 啦啦啦啦在线视频资源| 精品人妻一区二区三区麻豆| 边亲边吃奶的免费视频| 国产日韩一区二区三区精品不卡 | 国产精品久久久久久久久免| 国产成人一区二区在线| 欧美精品一区二区免费开放| 久久人人爽人人爽人人片va| 最近2019中文字幕mv第一页| 日韩一区二区视频免费看| 亚洲精品,欧美精品| 国产综合精华液| 美女视频免费永久观看网站| 欧美日韩一区二区视频在线观看视频在线| 又粗又硬又长又爽又黄的视频| 午夜激情av网站| av视频免费观看在线观看| 亚洲激情五月婷婷啪啪| 久久午夜综合久久蜜桃| 特大巨黑吊av在线直播| 亚洲成人av在线免费| 午夜激情久久久久久久| 国产日韩欧美视频二区| 国产探花极品一区二区| 草草在线视频免费看| 久久久久久久精品精品| 久久久a久久爽久久v久久| 欧美xxⅹ黑人| 国产一区亚洲一区在线观看| 国产 精品1| 日韩,欧美,国产一区二区三区| 美女主播在线视频| 免费不卡的大黄色大毛片视频在线观看| 欧美3d第一页| 赤兔流量卡办理| 国产极品天堂在线| 国产高清有码在线观看视频| 自线自在国产av| av有码第一页| 国产在线免费精品| 搡女人真爽免费视频火全软件| 好男人视频免费观看在线| 又黄又爽又刺激的免费视频.| 久久久久久久久久久久大奶| 免费观看a级毛片全部| 欧美少妇被猛烈插入视频| 人妻系列 视频| 亚洲精华国产精华液的使用体验| 亚洲色图 男人天堂 中文字幕 | av有码第一页| 最近最新中文字幕免费大全7| 久久久久国产精品人妻一区二区| 久久精品熟女亚洲av麻豆精品| 国产精品无大码| 久久午夜福利片| 久久女婷五月综合色啪小说| 亚洲欧美成人综合另类久久久| 一本一本综合久久| kizo精华| 国产精品久久久久久久电影| 国产欧美另类精品又又久久亚洲欧美| 国产成人aa在线观看| 在线观看人妻少妇| 纵有疾风起免费观看全集完整版| 一级二级三级毛片免费看| 免费黄频网站在线观看国产| a级毛片在线看网站| av在线播放精品| 精品久久久噜噜| 男人操女人黄网站| 最近手机中文字幕大全| 精品人妻一区二区三区麻豆| 日韩av在线免费看完整版不卡| 热re99久久国产66热| av专区在线播放| 亚洲精品乱码久久久久久按摩| 亚洲精品日韩在线中文字幕| 亚洲熟女精品中文字幕| 亚洲婷婷狠狠爱综合网| 两个人免费观看高清视频| www.色视频.com| av女优亚洲男人天堂| 久热这里只有精品99| 免费人妻精品一区二区三区视频| 欧美另类一区| 老熟女久久久| 久久热精品热| 亚洲精品av麻豆狂野| 精品久久久噜噜| 美女大奶头黄色视频| 综合色丁香网| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品无人区| 黄片播放在线免费| 天天影视国产精品| 婷婷色综合大香蕉| 欧美3d第一页| 免费观看a级毛片全部| 999精品在线视频| 制服诱惑二区| 国产视频内射| 精品人妻偷拍中文字幕| 久久国内精品自在自线图片| 夜夜看夜夜爽夜夜摸| 日韩欧美精品免费久久| 九九在线视频观看精品| 大香蕉97超碰在线| 国产女主播在线喷水免费视频网站| 亚洲精品日本国产第一区| 国产高清不卡午夜福利| 久久久国产一区二区| 伦理电影免费视频| 免费看光身美女| 亚洲欧美日韩卡通动漫| 亚洲国产日韩一区二区| 亚洲综合精品二区| 全区人妻精品视频| 一本一本综合久久| 国产午夜精品久久久久久一区二区三区| 美女福利国产在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av.av天堂| 成人综合一区亚洲| 丰满迷人的少妇在线观看| 人人妻人人澡人人爽人人夜夜| 午夜视频国产福利| 欧美性感艳星| 国产av码专区亚洲av| 精品人妻熟女av久视频| 美女主播在线视频| 精品久久久噜噜| 青青草视频在线视频观看| 美女脱内裤让男人舔精品视频| 少妇 在线观看| 日韩精品有码人妻一区| 免费高清在线观看日韩| 国内精品宾馆在线| 大香蕉久久成人网| 久久精品国产鲁丝片午夜精品| 一级a做视频免费观看| 国产精品99久久久久久久久| 一本大道久久a久久精品| 一区二区三区精品91| av在线app专区| 成人影院久久| 久久精品人人爽人人爽视色| 日韩 亚洲 欧美在线| 观看av在线不卡| 夜夜爽夜夜爽视频| 一区二区三区四区激情视频| 国产免费福利视频在线观看| 美女视频免费永久观看网站| 国产男人的电影天堂91| 亚洲国产日韩一区二区| 久久久久精品性色| 看非洲黑人一级黄片| 欧美精品一区二区免费开放| 中文字幕人妻熟人妻熟丝袜美| 中文字幕最新亚洲高清| 国产av码专区亚洲av| 建设人人有责人人尽责人人享有的| 亚洲欧美成人综合另类久久久| 一区二区日韩欧美中文字幕 | 一本大道久久a久久精品| 免费人妻精品一区二区三区视频| 在线观看人妻少妇| 国产成人av激情在线播放 | 亚洲高清免费不卡视频| 亚洲av国产av综合av卡| 伊人久久国产一区二区| 亚洲av电影在线观看一区二区三区| 波野结衣二区三区在线| 国产日韩欧美视频二区| 日韩大片免费观看网站| 99re6热这里在线精品视频| 99久国产av精品国产电影| 日本色播在线视频| 久久精品熟女亚洲av麻豆精品| 日日摸夜夜添夜夜爱| 国产精品 国内视频| 老熟女久久久| 成人18禁高潮啪啪吃奶动态图 | 天堂中文最新版在线下载| 秋霞在线观看毛片| 亚洲国产av影院在线观看| 日韩在线高清观看一区二区三区| 午夜福利影视在线免费观看| 国产精品一国产av| 美女中出高潮动态图| 国精品久久久久久国模美| 欧美精品人与动牲交sv欧美| 乱人伦中国视频| 国产日韩一区二区三区精品不卡 | 中文字幕亚洲精品专区| 中文精品一卡2卡3卡4更新| 九草在线视频观看| 国产成人aa在线观看| 国产成人91sexporn| 人妻一区二区av| 一本大道久久a久久精品| 国产精品一国产av| 日韩精品免费视频一区二区三区 | 中文乱码字字幕精品一区二区三区| av福利片在线| 国产av一区二区精品久久| 嫩草影院入口| 日本与韩国留学比较| 成人综合一区亚洲| 男男h啪啪无遮挡| 日产精品乱码卡一卡2卡三| 丁香六月天网| 91在线精品国自产拍蜜月| 日韩一本色道免费dvd| 中文字幕人妻熟人妻熟丝袜美| 色视频在线一区二区三区| 边亲边吃奶的免费视频| 亚洲av在线观看美女高潮| 人人澡人人妻人| 啦啦啦中文免费视频观看日本| av.在线天堂| 久久久久久久国产电影| 国产视频首页在线观看| 久久人妻熟女aⅴ| 久久久久久久国产电影| 久久午夜综合久久蜜桃| 日本黄大片高清| 黑人巨大精品欧美一区二区蜜桃 | 蜜桃国产av成人99| 亚洲欧洲日产国产| 精品国产露脸久久av麻豆| 如何舔出高潮| 亚洲精品国产色婷婷电影| 国产av码专区亚洲av| 欧美日韩亚洲高清精品| 亚洲欧洲精品一区二区精品久久久 | 国产成人精品久久久久久| 国产欧美日韩一区二区三区在线 | 亚洲av二区三区四区| 欧美丝袜亚洲另类| 九色成人免费人妻av| 99视频精品全部免费 在线| 欧美日韩亚洲高清精品| 国产精品一区二区在线不卡| 亚洲第一区二区三区不卡|