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

    甲狀腺癌中m6A甲基化調(diào)控因子的表達(dá)及其預(yù)后價(jià)值

    2021-09-09 06:35:18王文龍沈聰孫博韜白寧李新營陳嘉欣彭婀敏
    中國普通外科雜志 2021年8期
    關(guān)鍵詞:甲基化甲狀腺癌聚類

    王文龍,沈聰,孫博韜,白寧,李新營,陳嘉欣,彭婀敏

    (中南大學(xué)湘雅醫(yī)院1.甲狀腺外科3.國際醫(yī)療部,湖南長沙410008;2.中南大學(xué)湘雅醫(yī)學(xué)院,湖南長沙410013)

    過去幾十年,甲狀腺癌的發(fā)病率在全世界范圍內(nèi)迅速上升[1-2]。據(jù)國家癌癥中心數(shù)據(jù)顯示,甲狀腺癌已躋升至女性惡性腫瘤的第3 位,且發(fā)病出現(xiàn)年輕化趨勢(shì),已成為我國30 歲以下女性增長最快的惡性腫瘤[3]。甲狀腺乳頭狀癌是最常見的病理亞型,占85%以上[4-5]。手術(shù)是其主要治療方式,規(guī)范化手術(shù)治療后10年生存率達(dá)90%以上[6],但仍有1%~9%的患者在初次就診時(shí)就已出現(xiàn)遠(yuǎn)處轉(zhuǎn)移,嚴(yán)重影響了患者預(yù)后[7]。因此,闡述甲狀腺癌的發(fā)病機(jī)制對(duì)改善預(yù)后至關(guān)重要。

    m6A 甲基化修飾是指RNA 的腺苷酸的第6 位N處發(fā)生動(dòng)態(tài)可逆的甲基化過程,在調(diào)控mRNA 的轉(zhuǎn)錄、翻譯、穩(wěn)定、剪切、成熟等發(fā)揮了重要的生物學(xué)作用,其功能主要由甲基轉(zhuǎn)移酶(writer)、去甲基化酶(eraser)和識(shí)別蛋白(reader)決定[8-10]。越來越多的證據(jù)表明m6A 甲基化修飾在腫瘤的發(fā)生發(fā)展扮演了重要的角色[9,11-13]。METTL3 通過YTHDF2 依賴的途徑上調(diào)抑癌基因SOCS2 的m6A 甲基化水平,從而加速SOCS2 mRNA 的降解。METTL3 過表達(dá)促進(jìn)肝癌細(xì)胞的增殖和遷移,相反,敲低METTL3 顯著抑制腫瘤生長和轉(zhuǎn)移[14]。研究[15]揭示了IGF2BP3 在結(jié)直腸癌中高表達(dá),下調(diào)IGF2BP3 可以通過VEGF 和CCND1 的m6A 修飾來抑制血管生成和DNA 復(fù)制。以上研究證實(shí)m6A 甲基化調(diào)控因子高度參與癌癥發(fā)生發(fā)展,具有良好的潛在預(yù)后價(jià)值。然而,目前對(duì)于m6A 甲基化調(diào)控因子在甲狀腺癌致病調(diào)控機(jī)制尚不清楚。

    本研究利用癌癥基因組圖譜(The Cancer Genome Atlas,TCGA)數(shù)據(jù)庫中甲狀腺癌RNA-sep數(shù)據(jù),系統(tǒng)性分析了20 個(gè)m6A 甲基化調(diào)控因子在甲狀腺癌中的表達(dá)水平,再利用聚類分析和Lasso Cox 回歸分析篩選出預(yù)后相關(guān)基因及構(gòu)建風(fēng)險(xiǎn)預(yù)測(cè)模型并進(jìn)行危險(xiǎn)分層,為甲狀腺癌患者的臨床診療和預(yù)后改善提供了重要的理論參考依據(jù)。

    1 資料與方法

    1.1 數(shù)據(jù)的獲取與處理

    從TCGA 數(shù)據(jù)庫(https://portal.gdc.cancer.gov)下載了包含510 例甲狀腺癌組織和58 例正常組織樣本的RNA-seq 數(shù)據(jù)集及相應(yīng)臨床數(shù)據(jù)。臨床病理信息包括性別、年齡、腫瘤學(xué)多灶性、TNM 分期、淋巴結(jié)轉(zhuǎn)移、組織學(xué)類型、腫瘤位置和生存時(shí)間等,剔除RNA 序列以及臨床病理數(shù)據(jù)不全的樣本,共492 例癌組織和58 例正常組織樣本納入最終分析。所有患者的病理分期按照美國癌癥聯(lián)合委員會(huì)(AJCC)2018年1月頒布了第8 版甲狀腺癌TNM分期。根據(jù)之前m6A 甲基化在腫瘤中的研究[11,16],共篩選20 個(gè)m6A 甲基化調(diào)控因子,包括甲基轉(zhuǎn)移酶(WTAP、METTL3、METTL14、METTL16、KIAA1429、ZC3H13 和 RBM15)、識(shí) 別 蛋 白(HNRNPA2B1、IGF2BP1/2/3、YTHDC1/2、FMR1、LRPPRC 和YTHDF1/2/3)和去甲基化酶(ALKBH5和FTO)。應(yīng)用R 語言的“Limma”軟件包分析上述基因在癌與正常樣本的差異表達(dá)。

    1.2 聚類分析

    為了從功能上闡明甲狀腺癌中m6A 調(diào)控因子的生物學(xué)特性,使用“consensusclusterplus”軟件包對(duì)甲狀腺癌患者進(jìn)行一致性聚類分析并將患者分成cluster 1 和cluster 2 組。主成分分析(PCA)檢測(cè)基于m6A 調(diào)控因子表達(dá)是否能將兩組患者明顯區(qū)分,并探討不同的臨床病理因素和總體生存率在兩組聚類中是否存在差異。

    1.3 風(fēng)險(xiǎn)預(yù)測(cè)模型的建立與評(píng)估

    為了探討m6A 調(diào)控因子在甲狀腺癌預(yù)后的作用,通過“glmnet”軟件包進(jìn)行Lasso Cox 回歸分析篩選預(yù)后相關(guān)的基因,并計(jì)算基因的相關(guān)回歸系數(shù),根據(jù)回歸系數(shù)加權(quán)建立甲狀腺癌預(yù)后風(fēng)險(xiǎn)評(píng)估公式,risk score=Coei*Expi,Coei表示回歸系數(shù),Expi表示m6A 調(diào)控因子的表達(dá)量。據(jù)此公式計(jì)算每個(gè)患者的風(fēng)險(xiǎn)評(píng)分,以風(fēng)險(xiǎn)評(píng)分的中位數(shù)將患者分為高危組和低危組,比較兩組總體生存率是否存在差異,ROC 曲線下面積(AUC)評(píng)估模型的預(yù)測(cè)能力。

    1.4 統(tǒng)計(jì)學(xué)處理

    采用R 3.6.0 軟件進(jìn)行統(tǒng)計(jì)學(xué)分析,Wilcoxon 秩和檢驗(yàn)分析20 個(gè)m6A 甲基化調(diào)控因子在腫瘤和正常組織的差異表達(dá)。計(jì)量資料用t檢驗(yàn),計(jì)數(shù)資料用χ2檢驗(yàn)。采用Spearman 相關(guān)分析基因間的相關(guān)性,0.1≤|r|≤1.0 時(shí)定義為存在相關(guān);生存曲線采用Kaplan-Meier 法和Log-rank 檢驗(yàn),P<0.05 為差異有統(tǒng)計(jì)學(xué)意義。

    2 結(jié)果

    2.1 一般資料及m6A 甲基化調(diào)控因子的基因表達(dá)與相關(guān)性

    本研究共納入492 例甲狀腺癌患者,平均年齡47.29 歲,男女比列1∶2.76;臨床病理TNM 分期:I 期276 例,II 期52 例;III 期109 例,IV 期55 例;T 分期:T1期135例,T2期162例,T3期172例,T4期23例;44.9%的患者出現(xiàn)頸部淋巴結(jié)轉(zhuǎn)移;9 例患者出現(xiàn)遠(yuǎn)處轉(zhuǎn)移;228 例腫瘤呈現(xiàn)出多灶性,占46.3%;腫瘤位于峽部、左側(cè)、右側(cè)、雙側(cè)分別為22、172、214,84 例(表1)。

    表1 492例甲狀腺癌患者基本臨床病理資料Table 1 The clinicopathologic information of the 492 patients with thyroid cancer

    19 個(gè)m6A 甲基化調(diào)控因子在甲狀腺癌和正常組織表達(dá)具有統(tǒng)計(jì)學(xué)差異(P<0.05),其中HNRNPC(P<0.001)、IGF2BP2(P<0.001)、FMR1(P<0.05)在甲狀腺癌組織中明顯高表達(dá)。相反,ALKBH5(P<0.001)、FTO(P<0.001)、METTL3(P<0.001)、METTL14(P<0.001)、METTL16(P<0.05)、WTAP(P<0.001)、YTHDF1(P<0.001)、YTHDF3(P<0.001)、YTHDC2(P<0.001)、YTHDC1(P<0.001)、ZC3H13(P<0.001)、HNRNPA2B1(P<0.001)、RBM15(P<0.001)、IGF2BP1(P<0.001)、IGF2BP3(P<0.001)和LRPPRC(P<0.001)在甲狀腺癌標(biāo)本中表達(dá)明顯下調(diào)(圖1A)。此外,大部分m6A 甲基化調(diào)控因子都顯示正相關(guān),其中,METTL14 與YTHDC1、YTHDF3與ZC3H13 呈現(xiàn)最強(qiáng)正相關(guān)(r=0.73);而ALKBH5 與IGF2BP2呈明顯負(fù)相關(guān)(r=-0.33)(圖1B)。為了進(jìn)一步了解19 個(gè)m6A 甲基化調(diào)控因子相互關(guān)系,PPI 構(gòu)建了網(wǎng)絡(luò)相互作用網(wǎng)(圖1C)。

    圖1 m6A甲基化調(diào)控因子的基因表達(dá)及相關(guān)性A:20個(gè)m6A甲基化調(diào)控因子在甲狀腺癌和正常組織表達(dá)水平;B:Spearman相關(guān)分析m6A甲基化調(diào)控因子相關(guān)性;C:PPI網(wǎng)絡(luò)評(píng)估m(xù)6A甲基化調(diào)控因子相互作用關(guān)系Figure 1 Gene expression and correlations of the m6A methylation regulatorsA:The mRNA levels of the 20 m6A regulators in thyroid cancer; B:Spearman correlation analysis of co-expressions among the m6A regulators; C:PPI network of the m6A regulators

    2.2 一致性聚類分析m6A甲基化調(diào)控因子

    為了確定最佳的聚類數(shù),采用一致性聚類分析對(duì)492 例甲狀腺癌患者進(jìn)行聚類分組,共識(shí)矩陣k 值為2 時(shí),甲狀腺癌樣本之間不存在交叉,聚類穩(wěn)定性好(圖2A)。因此,492 例甲狀腺癌患者被分成cluster 1 和cluster 2 組,主成分分析可以將兩組患者明顯區(qū)分(圖2B)。聚類熱圖示,cluster 1組頸淋巴結(jié)轉(zhuǎn)移發(fā)生率明顯高于cluster 2(P<0.01),而TNM 分期、T 分期、M 分期、腫瘤多灶性、腫瘤位置、年齡、性別在兩組之間無差異(均P>0.05)(圖2C)。生存分析示,cluster 1 生存時(shí)間低于cluster 2(P=0.048)(圖2D)。

    圖2 一致性聚類分析m6A甲基化調(diào)控因子A:k值為2時(shí),甲狀腺癌樣本之間不存在交叉,聚類穩(wěn)定性好;B:主成分分析區(qū)分cluster 1和cluster 2;C:不同的聚類與臨床病理的關(guān)系;D:總體生存率在兩組聚類的情況Figure 2 Consensus clustering analysis of the m6A RNA methylation regulatorsA:The most appropriate selection with clustering stability when k=2 used; B:Principal component analysis to distinguish cluster 1 and cluster 2; C:Heatmap and clinicopathologic features of the two clusters;D:The overall survival of cluster 1 and cluster 2

    2.3 m6A甲基化調(diào)控因子構(gòu)建風(fēng)險(xiǎn)預(yù)后模型及檢驗(yàn)

    Lasso Cox 回歸分析篩選4 個(gè)基因(IGF2BP2、RBM15、YTHDF1、YTHDF3),風(fēng)險(xiǎn)評(píng)分=(-0.129×IGF2BP2)+(0.308×RBM15)+(0.331×YTHDF1)+(0.813×YTHDF3)(圖3A),根據(jù)回歸系數(shù)計(jì)算每個(gè)患者的風(fēng)險(xiǎn)評(píng)分,以風(fēng)險(xiǎn)評(píng)分的中位數(shù)將患者分為高危組和低危組,風(fēng)險(xiǎn)評(píng)分和生存狀態(tài)分布如圖3B 所示,高風(fēng)險(xiǎn)患者中死亡人數(shù)明顯增加。相比低風(fēng)險(xiǎn)組患者,高風(fēng)險(xiǎn)患者生存期明顯縮短(P=0.007)(圖3C),ROC 曲線下面積示該模型可預(yù)測(cè)甲狀腺癌患者的預(yù)后(AUC=0.731)(圖3D)。

    圖3 m6A甲基化調(diào)控因子構(gòu)建風(fēng)險(xiǎn)預(yù)后模型及檢驗(yàn)A:Lasso Cox回歸模型的m6A基因及相關(guān)系數(shù);B:風(fēng)險(xiǎn)評(píng)分和生存狀態(tài)分布;C:高、低危組的生存曲線;D:ROC曲線評(píng)估模型的預(yù)測(cè)效能Figure 3 Construction od prognostic risk model and validationA:Coefficients calculated by Lasso Cox regression analysis;B:Risk score and survival status distribution;C:Kaplan-Meier analysis of overall survival in high-risk and low-risk groups;D:Assessment of the prediction accuracy by ROC curve

    3 討論

    甲狀腺癌是一種常見的內(nèi)分泌惡性腫瘤,發(fā)病率在世界范圍內(nèi)呈快速上升趨勢(shì),絕大部分甲狀腺癌預(yù)后良好[17-18]。然而一些分化較差的甲狀腺癌和高危型甲狀腺乳頭狀癌具有高度侵襲性,易出現(xiàn)局部侵犯和遠(yuǎn)處轉(zhuǎn)移,且尚無特效治療方法[19-20]。因此闡述甲狀腺癌發(fā)生和演變進(jìn)展的分子機(jī)制對(duì)疾病的早期干預(yù)和預(yù)后預(yù)測(cè)具有重要的臨床應(yīng)用價(jià)值。

    m6A 修飾是近年來的研究熱點(diǎn),其作用主要參與轉(zhuǎn)錄后修飾,同時(shí)也參與了DNA 修復(fù)、細(xì)胞分化、細(xì)胞重編程、細(xì)胞應(yīng)激等生命活動(dòng)[21-22]。m6A 甲基化在不同的腫瘤中均有報(bào)道,在甲狀腺癌中m6A 相關(guān)基因的作用仍不清楚,Ye 等[23]發(fā)現(xiàn)IGF2BP2 在甲狀腺癌中高表達(dá),lncRNA MALAT1 通過競爭性結(jié)合miR-204,通過m6A 修飾識(shí)別,上調(diào)IGF2BP2,增強(qiáng)MYC 表達(dá),從而刺激甲狀腺癌細(xì)胞增殖、遷移和侵襲,同時(shí)減弱腫瘤生長和細(xì)胞凋亡。然而該研究樣本量有限無法驗(yàn)證lncRNA MALAT1/mir-204/IGF2BP2/m6A-MYC 軸在甲狀腺癌進(jìn)展機(jī)制。Wang 等[24]報(bào)道METTL3 作為致癌基因通過上調(diào)TCF1 的m6A 甲基化水平促進(jìn)甲狀腺癌的增值和侵襲。然而He 等[25]研究揭示METTL3 作為抑癌基因結(jié)合YTHDF2 通過c-Rel 和RelA 失活NF-κB 通路改變腫瘤相關(guān)的中性粒細(xì)胞浸潤來調(diào)節(jié)腫瘤生長。同樣,Hou 等[26]通過生物信息學(xué)證實(shí)METTL3 在甲狀腺癌組織中明顯低表達(dá)。上述研究證實(shí)m6A 相關(guān)基因表達(dá)水平和功能機(jī)制在甲狀腺癌中仍需進(jìn)一步闡述。本研究初步驗(yàn)證了METTL3 在正常組織中表達(dá)量明顯高于甲狀腺癌組織(P<0.001),提示METTL3 可能是一種抑癌基因,可改善甲狀腺癌患者的預(yù)后,這也進(jìn)一步驗(yàn)證了He 等[25-26]的研究結(jié)果。

    為了從功能上闡明甲狀腺癌中m6A 調(diào)控因子的生物學(xué)特性,本研究對(duì)492 例甲狀腺癌患者進(jìn)行一致性聚類分析后并分成cluster 1 和cluster 2,cluster 1 生存期低于cluster 2(P<0.05),相反,頸淋巴結(jié)轉(zhuǎn)移發(fā)生率明顯高于cluster 2(P<0.01)。暗示m6A 調(diào)控因子能預(yù)測(cè)頸淋巴結(jié)轉(zhuǎn)移并影響甲狀腺癌患者預(yù)后,具有非常重要的臨床實(shí)際意義。聚類分析在肝癌[27]、腎癌[16]、乳腺癌[11]、胃癌[28]等惡性腫瘤預(yù)后分析中均起到較好的作用。

    此外,為了探討m6A 調(diào)控因子在甲狀腺癌預(yù)后的作用,應(yīng)用Lasso Cox 回歸分析篩選的4 個(gè)基因(IGF2BP2、RBM15、YTHDF1、YTHDF3)構(gòu)建風(fēng)險(xiǎn)評(píng)估模型,相比IGF2BP2,其余3 個(gè)基因在甲狀腺癌的研究中甚少,RBM15 不具有甲基轉(zhuǎn)移酶的催化功能,但它可以結(jié)合METTL3 和WTAP,并將這兩個(gè)蛋白引導(dǎo)到特定的RNA 位點(diǎn)進(jìn)行m6A 修飾[8,29];而YTHDF1 和YTHDF3 均屬于YTH 域家族,YTHDF1 通過與起始因子相互作用增強(qiáng)mRNA 翻譯和蛋白質(zhì)合成,YTHDF3 分別通過與YTHDF1 相互作用增強(qiáng)RNA 翻譯,與YTHDF2 結(jié)合促進(jìn)RNA 降解[9,30]。相比低風(fēng)險(xiǎn)組患者,高風(fēng)險(xiǎn)患者生存期明顯縮短(P=0.007),ROC 曲線示該模型可以預(yù)測(cè)甲狀腺癌患者的預(yù)后(AUC=0.731)。相比TNM 分期和MACIS 預(yù)后評(píng)估系統(tǒng),該模型基于RNA 測(cè)序數(shù)據(jù),更符合當(dāng)今的精準(zhǔn)醫(yī)學(xué)理念。

    同時(shí),本研究也存在一些缺陷,首先,該研究基于生物信息學(xué)構(gòu)建的預(yù)后預(yù)測(cè)模型,無外部的實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證,臨床實(shí)際應(yīng)用價(jià)值有限。其次,該研究只是初步探討了m6A 調(diào)控因子在甲狀腺癌預(yù)后作用,不能完整的闡述甲狀腺癌發(fā)生和演變的分子機(jī)制。最后,本研究基于492 例甲狀腺癌患者,樣本量不大,后續(xù)需要進(jìn)行更大的獨(dú)立隊(duì)列研究來驗(yàn)證。

    綜上,本研究基于TCGA 數(shù)據(jù)庫分析甲狀腺癌m6A 調(diào)控因子表達(dá)水平,通過Lasso Cox 回歸分析篩選的4 個(gè)基因(IGF2BP2、RBM15、YTHDF1、YTHDF3)構(gòu)建風(fēng)險(xiǎn)評(píng)估模型,AUC 提示該模型具有良好的預(yù)測(cè)效能,將為甲狀腺癌臨床診治提供潛在的理論依據(jù)。

    猜你喜歡
    甲基化甲狀腺癌聚類
    分化型甲狀腺癌切除術(shù)后多發(fā)骨轉(zhuǎn)移一例
    分化型甲狀腺癌肺轉(zhuǎn)移的研究進(jìn)展
    基于DBSACN聚類算法的XML文檔聚類
    全甲狀腺切除術(shù)治療甲狀腺癌適應(yīng)證選擇及并發(fā)癥防治
    基于改進(jìn)的遺傳算法的模糊聚類算法
    一種層次初始的聚類個(gè)數(shù)自適應(yīng)的聚類方法研究
    鼻咽癌組織中SYK基因啟動(dòng)子區(qū)的甲基化分析
    胃癌DNA甲基化研究進(jìn)展
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    全甲基化沒食子兒茶素沒食子酸酯的制備
    叶爱在线成人免费视频播放| 在线观看国产h片| 久久韩国三级中文字幕| 久久精品夜色国产| 久久久久久久亚洲中文字幕| 色网站视频免费| av免费观看日本| 午夜久久久在线观看| 中文精品一卡2卡3卡4更新| 精品少妇内射三级| 中文字幕制服av| 欧美老熟妇乱子伦牲交| 边亲边吃奶的免费视频| 亚洲综合色网址| 久久精品国产综合久久久| 国产精品久久久久久精品电影小说| 美女视频免费永久观看网站| 精品一品国产午夜福利视频| 99精国产麻豆久久婷婷| 日韩中字成人| 丰满乱子伦码专区| 亚洲国产精品成人久久小说| 久久婷婷青草| 美女大奶头黄色视频| 18禁裸乳无遮挡动漫免费视频| 国产免费视频播放在线视频| 欧美成人精品欧美一级黄| 1024香蕉在线观看| 电影成人av| 十分钟在线观看高清视频www| 国产1区2区3区精品| 91精品伊人久久大香线蕉| 国产精品一二三区在线看| 一区二区日韩欧美中文字幕| 日韩人妻精品一区2区三区| 精品午夜福利在线看| 欧美日韩国产mv在线观看视频| 飞空精品影院首页| www.av在线官网国产| 国产精品无大码| 亚洲欧洲精品一区二区精品久久久 | 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久av不卡| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲欧美一区二区av| 国产极品粉嫩免费观看在线| 高清不卡的av网站| 九九爱精品视频在线观看| 捣出白浆h1v1| 99热国产这里只有精品6| 中国国产av一级| 99热全是精品| 亚洲,欧美,日韩| 国产极品天堂在线| 免费av中文字幕在线| a级毛片黄视频| 国产一区有黄有色的免费视频| 国产男女超爽视频在线观看| 国产成人免费观看mmmm| 精品卡一卡二卡四卡免费| 狠狠精品人妻久久久久久综合| 欧美+日韩+精品| 超色免费av| 欧美国产精品一级二级三级| 少妇熟女欧美另类| 亚洲,欧美,日韩| 黄片无遮挡物在线观看| 久久久久久久久免费视频了| 18禁动态无遮挡网站| 久久久久国产网址| 精品久久久久久电影网| 少妇猛男粗大的猛烈进出视频| 欧美日韩亚洲高清精品| 美女高潮到喷水免费观看| 99香蕉大伊视频| 色哟哟·www| 国产精品久久久久久精品电影小说| 青春草视频在线免费观看| 寂寞人妻少妇视频99o| 日韩三级伦理在线观看| 国产精品熟女久久久久浪| 男女午夜视频在线观看| 国产免费福利视频在线观看| 免费久久久久久久精品成人欧美视频| 亚洲欧美一区二区三区国产| 久久毛片免费看一区二区三区| 亚洲精品久久久久久婷婷小说| 最近手机中文字幕大全| 少妇精品久久久久久久| 中文天堂在线官网| 青春草亚洲视频在线观看| 国产日韩欧美视频二区| 久久鲁丝午夜福利片| 亚洲av国产av综合av卡| 日韩成人av中文字幕在线观看| 婷婷成人精品国产| 五月开心婷婷网| 国产一区二区在线观看av| 欧美日韩国产mv在线观看视频| 国产野战对白在线观看| 中文天堂在线官网| 观看美女的网站| 欧美成人精品欧美一级黄| 国产精品.久久久| 视频区图区小说| 一二三四中文在线观看免费高清| 免费黄频网站在线观看国产| 欧美精品一区二区大全| 久久久精品94久久精品| 国产成人欧美| 免费大片黄手机在线观看| 国产国语露脸激情在线看| 热re99久久国产66热| 极品少妇高潮喷水抽搐| 国产成人精品无人区| 日韩av免费高清视频| 国产 精品1| 午夜免费男女啪啪视频观看| 欧美bdsm另类| 精品国产国语对白av| 欧美av亚洲av综合av国产av | 日韩不卡一区二区三区视频在线| 免费在线观看视频国产中文字幕亚洲 | 91午夜精品亚洲一区二区三区| 日韩av免费高清视频| 久久人人爽人人片av| 亚洲男人天堂网一区| 波多野结衣一区麻豆| 国产成人a∨麻豆精品| 亚洲少妇的诱惑av| 久久久久久久久久久久大奶| 国产成人av激情在线播放| 欧美成人精品欧美一级黄| 亚洲av电影在线进入| 一本色道久久久久久精品综合| 国产精品嫩草影院av在线观看| 成年动漫av网址| 91午夜精品亚洲一区二区三区| 国产精品女同一区二区软件| 熟妇人妻不卡中文字幕| 日本-黄色视频高清免费观看| 精品久久蜜臀av无| 一区二区三区激情视频| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久 | 99久久精品国产国产毛片| 欧美在线黄色| 中文天堂在线官网| 久久 成人 亚洲| 涩涩av久久男人的天堂| 午夜日本视频在线| 久久精品国产综合久久久| 欧美97在线视频| 91午夜精品亚洲一区二区三区| 色视频在线一区二区三区| 亚洲精品国产av蜜桃| 大陆偷拍与自拍| 亚洲综合色网址| 青春草国产在线视频| 欧美激情极品国产一区二区三区| 爱豆传媒免费全集在线观看| 国产精品99久久99久久久不卡 | 电影成人av| 日韩av免费高清视频| 亚洲婷婷狠狠爱综合网| 国产精品国产av在线观看| 男人操女人黄网站| 免费日韩欧美在线观看| 亚洲第一青青草原| 极品人妻少妇av视频| www日本在线高清视频| 日韩欧美一区视频在线观看| 在线观看免费日韩欧美大片| 精品卡一卡二卡四卡免费| 男女无遮挡免费网站观看| 亚洲精品久久久久久婷婷小说| xxxhd国产人妻xxx| 国产色婷婷99| 波多野结衣一区麻豆| 婷婷色综合www| 国产不卡av网站在线观看| 婷婷色综合大香蕉| 99久国产av精品国产电影| 成人午夜精彩视频在线观看| 国产精品秋霞免费鲁丝片| 五月天丁香电影| 伦精品一区二区三区| 七月丁香在线播放| 久久久久久免费高清国产稀缺| 夫妻性生交免费视频一级片| 欧美国产精品一级二级三级| 国产日韩欧美亚洲二区| 欧美激情高清一区二区三区 | 高清不卡的av网站| 日韩一区二区视频免费看| 91aial.com中文字幕在线观看| 国语对白做爰xxxⅹ性视频网站| 久久这里有精品视频免费| 91aial.com中文字幕在线观看| 精品久久久精品久久久| 欧美+日韩+精品| 伦理电影免费视频| www.自偷自拍.com| 大香蕉久久成人网| 亚洲国产精品成人久久小说| 久久国产精品大桥未久av| 国产精品欧美亚洲77777| 18禁观看日本| 亚洲人成电影观看| 久久午夜福利片| 久久久精品免费免费高清| 国产极品粉嫩免费观看在线| 国产黄色视频一区二区在线观看| 看免费成人av毛片| 久久ye,这里只有精品| www.自偷自拍.com| 欧美日韩视频精品一区| 国产 精品1| 亚洲中文av在线| 国产亚洲午夜精品一区二区久久| 精品亚洲乱码少妇综合久久| 免费观看a级毛片全部| 老司机影院成人| 捣出白浆h1v1| 久久久久人妻精品一区果冻| 国产精品久久久av美女十八| av福利片在线| 一边摸一边做爽爽视频免费| 精品酒店卫生间| 亚洲熟女精品中文字幕| 精品少妇久久久久久888优播| 在线观看www视频免费| 亚洲精品一二三| 亚洲精品视频女| 日本猛色少妇xxxxx猛交久久| 五月天丁香电影| 国产成人精品一,二区| 国产高清不卡午夜福利| 中文字幕色久视频| 99久国产av精品国产电影| videos熟女内射| 久久女婷五月综合色啪小说| 日日撸夜夜添| 亚洲国产色片| h视频一区二区三区| 午夜福利在线观看免费完整高清在| 色婷婷av一区二区三区视频| 91精品伊人久久大香线蕉| 久久久久久人人人人人| 女性生殖器流出的白浆| 一本—道久久a久久精品蜜桃钙片| av卡一久久| 黑人巨大精品欧美一区二区蜜桃| 天堂中文最新版在线下载| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩av久久| 亚洲中文av在线| 国产精品 欧美亚洲| 欧美 日韩 精品 国产| 国产成人精品无人区| 卡戴珊不雅视频在线播放| 一级,二级,三级黄色视频| 久久国产亚洲av麻豆专区| 考比视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 飞空精品影院首页| 青草久久国产| 国产精品99久久99久久久不卡 | 国产黄色免费在线视频| 国产精品熟女久久久久浪| 国产一区二区激情短视频 | 欧美人与善性xxx| 久久久国产一区二区| 久久人人97超碰香蕉20202| 91午夜精品亚洲一区二区三区| 日韩欧美精品免费久久| 欧美激情极品国产一区二区三区| 久久久久久久久免费视频了| 国产又色又爽无遮挡免| 国产男女内射视频| 天天影视国产精品| 日韩大片免费观看网站| 日日啪夜夜爽| 精品国产一区二区久久| 黄色 视频免费看| 欧美精品一区二区大全| 老司机亚洲免费影院| www.av在线官网国产| 中文字幕最新亚洲高清| 欧美变态另类bdsm刘玥| av网站在线播放免费| 亚洲精品一区蜜桃| 日韩伦理黄色片| 少妇人妻精品综合一区二区| 老鸭窝网址在线观看| 亚洲精品国产av成人精品| 不卡视频在线观看欧美| 成人漫画全彩无遮挡| 国产色婷婷99| av又黄又爽大尺度在线免费看| 九九爱精品视频在线观看| 欧美精品一区二区免费开放| 精品少妇一区二区三区视频日本电影 | 亚洲视频免费观看视频| 午夜久久久在线观看| 日韩制服丝袜自拍偷拍| 99热国产这里只有精品6| 天堂俺去俺来也www色官网| 国产成人精品久久二区二区91 | www日本在线高清视频| 亚洲av电影在线进入| 在线观看美女被高潮喷水网站| 久久久精品免费免费高清| 色婷婷av一区二区三区视频| 欧美日韩亚洲高清精品| 婷婷色综合大香蕉| 晚上一个人看的免费电影| av在线观看视频网站免费| 免费在线观看视频国产中文字幕亚洲 | 只有这里有精品99| 人妻系列 视频| 边亲边吃奶的免费视频| 人妻人人澡人人爽人人| 黄色配什么色好看| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 有码 亚洲区| 久久精品国产鲁丝片午夜精品| 亚洲天堂av无毛| 一边摸一边做爽爽视频免费| 久久久亚洲精品成人影院| 2021少妇久久久久久久久久久| 久久精品人人爽人人爽视色| 啦啦啦中文免费视频观看日本| 熟女av电影| 久久久精品94久久精品| 亚洲人成网站在线观看播放| 多毛熟女@视频| 国产淫语在线视频| 欧美日韩亚洲国产一区二区在线观看 | 观看av在线不卡| 国产黄色视频一区二区在线观看| 天天操日日干夜夜撸| 亚洲欧美清纯卡通| 在线天堂最新版资源| 99久久综合免费| 国产成人一区二区在线| 欧美另类一区| 伊人久久大香线蕉亚洲五| 叶爱在线成人免费视频播放| 精品一区在线观看国产| 18禁观看日本| 久久免费观看电影| 国产免费一区二区三区四区乱码| 欧美变态另类bdsm刘玥| 亚洲美女视频黄频| 97精品久久久久久久久久精品| 人人妻人人爽人人添夜夜欢视频| 99热网站在线观看| 麻豆乱淫一区二区| 人体艺术视频欧美日本| 黑人欧美特级aaaaaa片| 午夜影院在线不卡| 久久精品国产自在天天线| 美国免费a级毛片| 成年女人在线观看亚洲视频| 如日韩欧美国产精品一区二区三区| 久久午夜福利片| 一区二区日韩欧美中文字幕| 亚洲av.av天堂| 一级黄片播放器| 国产白丝娇喘喷水9色精品| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久久久免| 久久久国产精品麻豆| 亚洲精品av麻豆狂野| 大香蕉久久网| 亚洲精品av麻豆狂野| 免费大片黄手机在线观看| 久久久久久久精品精品| 女性生殖器流出的白浆| 午夜老司机福利剧场| 亚洲成人av在线免费| 在线观看免费高清a一片| 欧美激情极品国产一区二区三区| 国产免费视频播放在线视频| 精品国产一区二区三区久久久樱花| 天天躁日日躁夜夜躁夜夜| 熟女电影av网| 男男h啪啪无遮挡| 亚洲国产日韩一区二区| 亚洲国产欧美网| 欧美最新免费一区二区三区| 人妻人人澡人人爽人人| 国产福利在线免费观看视频| 免费黄色在线免费观看| 新久久久久国产一级毛片| 久久狼人影院| 欧美中文综合在线视频| h视频一区二区三区| 亚洲欧洲日产国产| 亚洲精华国产精华液的使用体验| 日韩成人av中文字幕在线观看| 国产一区二区三区综合在线观看| 国产1区2区3区精品| 90打野战视频偷拍视频| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 最近的中文字幕免费完整| 汤姆久久久久久久影院中文字幕| 啦啦啦啦在线视频资源| 看非洲黑人一级黄片| 18禁观看日本| 99久久中文字幕三级久久日本| 老熟女久久久| xxx大片免费视频| 男人添女人高潮全过程视频| 激情五月婷婷亚洲| 欧美激情 高清一区二区三区| 9色porny在线观看| 女性生殖器流出的白浆| 高清欧美精品videossex| 国产免费福利视频在线观看| 亚洲欧美日韩另类电影网站| 亚洲美女黄色视频免费看| 国产成人欧美| 性少妇av在线| 国产淫语在线视频| 久久热在线av| 五月开心婷婷网| 亚洲综合精品二区| 天天躁夜夜躁狠狠久久av| 欧美成人午夜精品| 26uuu在线亚洲综合色| 97在线人人人人妻| 国产精品免费大片| 亚洲成av片中文字幕在线观看 | 亚洲人成网站在线观看播放| 国产日韩一区二区三区精品不卡| 欧美日本中文国产一区发布| 天天躁日日躁夜夜躁夜夜| 五月开心婷婷网| 国语对白做爰xxxⅹ性视频网站| 满18在线观看网站| 人妻少妇偷人精品九色| 在线观看免费日韩欧美大片| 国产在线免费精品| av国产久精品久网站免费入址| videosex国产| 高清欧美精品videossex| 国产成人精品久久二区二区91 | 纯流量卡能插随身wifi吗| 国产极品粉嫩免费观看在线| 日本91视频免费播放| 国产精品人妻久久久影院| 成年女人在线观看亚洲视频| 青青草视频在线视频观看| 黄色一级大片看看| 日韩免费高清中文字幕av| 成年美女黄网站色视频大全免费| 亚洲五月色婷婷综合| 91国产中文字幕| 自线自在国产av| 精品久久蜜臀av无| 18+在线观看网站| 91在线精品国自产拍蜜月| 午夜福利在线免费观看网站| 老司机影院成人| 亚洲一区二区三区欧美精品| 亚洲三区欧美一区| 一区在线观看完整版| 亚洲美女搞黄在线观看| 亚洲一级一片aⅴ在线观看| 男人添女人高潮全过程视频| 久久影院123| 免费大片黄手机在线观看| 亚洲国产最新在线播放| 日韩av不卡免费在线播放| 啦啦啦啦在线视频资源| 久久久久久伊人网av| 日本欧美视频一区| av天堂久久9| 女性生殖器流出的白浆| 欧美 亚洲 国产 日韩一| 精品人妻偷拍中文字幕| 国产一级毛片在线| 亚洲欧洲国产日韩| 欧美精品高潮呻吟av久久| 国产人伦9x9x在线观看 | 丰满少妇做爰视频| 男女啪啪激烈高潮av片| 欧美 日韩 精品 国产| 一本色道久久久久久精品综合| 秋霞伦理黄片| 啦啦啦在线免费观看视频4| 成人亚洲精品一区在线观看| 国产精品嫩草影院av在线观看| 国产黄色免费在线视频| 亚洲精品国产av蜜桃| 十分钟在线观看高清视频www| 日韩一区二区视频免费看| 91在线精品国自产拍蜜月| 最近中文字幕2019免费版| 人妻少妇偷人精品九色| 国产熟女午夜一区二区三区| 伊人久久国产一区二区| www.精华液| 人妻一区二区av| 国产精品久久久久久av不卡| 国产女主播在线喷水免费视频网站| 韩国高清视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 男人添女人高潮全过程视频| 国产成人av激情在线播放| 国产av精品麻豆| 午夜激情av网站| 一区在线观看完整版| 波多野结衣av一区二区av| 国产精品 欧美亚洲| 美女高潮到喷水免费观看| 午夜福利影视在线免费观看| 亚洲av国产av综合av卡| 亚洲美女搞黄在线观看| 国产精品免费大片| 女性生殖器流出的白浆| 9热在线视频观看99| 亚洲天堂av无毛| 午夜福利视频在线观看免费| 国产精品一区二区在线不卡| 波野结衣二区三区在线| 国产片特级美女逼逼视频| 欧美人与性动交α欧美精品济南到 | 亚洲国产日韩一区二区| 婷婷色综合www| 中文天堂在线官网| 亚洲国产成人一精品久久久| 自拍欧美九色日韩亚洲蝌蚪91| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久国产精品男人的天堂亚洲| 男女高潮啪啪啪动态图| 国产白丝娇喘喷水9色精品| 在线 av 中文字幕| 亚洲天堂av无毛| 久久国产精品男人的天堂亚洲| 午夜免费鲁丝| 欧美激情极品国产一区二区三区| a 毛片基地| 日本色播在线视频| 男人添女人高潮全过程视频| 黄色毛片三级朝国网站| 成人毛片60女人毛片免费| 性色av一级| 色婷婷av一区二区三区视频| 尾随美女入室| 99久国产av精品国产电影| 国产精品无大码| 高清在线视频一区二区三区| √禁漫天堂资源中文www| 90打野战视频偷拍视频| 一二三四在线观看免费中文在| 国产国语露脸激情在线看| 人妻 亚洲 视频| 亚洲精品av麻豆狂野| 久久久久久久久久久久大奶| 中文字幕亚洲精品专区| av网站在线播放免费| 成人毛片60女人毛片免费| 中文字幕制服av| 日本猛色少妇xxxxx猛交久久| 看十八女毛片水多多多| 国产黄色视频一区二区在线观看| 精品一区二区免费观看| 少妇被粗大的猛进出69影院| 性色avwww在线观看| 亚洲精品,欧美精品| 免费黄频网站在线观看国产| 亚洲精品久久午夜乱码| 国产女主播在线喷水免费视频网站| 女的被弄到高潮叫床怎么办| 午夜激情久久久久久久| 国产一区二区三区综合在线观看| 人妻少妇偷人精品九色| 精品一区二区免费观看| 一区福利在线观看| 男男h啪啪无遮挡| 1024视频免费在线观看| 免费黄频网站在线观看国产| 久久精品久久精品一区二区三区| 国产女主播在线喷水免费视频网站| 自线自在国产av| 亚洲精品国产色婷婷电影| 老司机影院毛片| 超碰97精品在线观看| 亚洲av福利一区| 91精品三级在线观看| 国产 精品1| 极品人妻少妇av视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美变态另类bdsm刘玥| 黄色怎么调成土黄色| 久久久久久久久久久久大奶| 亚洲av在线观看美女高潮| 青春草亚洲视频在线观看| 亚洲国产精品一区二区三区在线| 免费看av在线观看网站| av网站在线播放免费| 久久久国产精品麻豆| 久久久精品区二区三区| 亚洲av福利一区| 性高湖久久久久久久久免费观看| 免费黄色在线免费观看| 可以免费在线观看a视频的电影网站 | 美女高潮到喷水免费观看| 亚洲欧美精品综合一区二区三区 | 999精品在线视频|