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

    單細胞轉錄組測序分析結腸癌微環(huán)境中免疫細胞亞群與癌癥進程的關系

    2021-08-16 09:48:18周彥彤李春曉王勁松孫芳洲許東奎張雪燕袁芃錢海利
    解放軍醫(yī)學雜志 2021年7期
    關鍵詞:單細胞亞群結腸癌

    周彥彤,李春曉,王勁松,孫芳洲,許東奎,張雪燕,袁芃,錢海利*

    1國家癌癥中心/國家腫瘤臨床醫(yī)學研究中心/中國醫(yī)學科學院北京協(xié)和醫(yī)學院腫瘤醫(yī)院分子腫瘤學國家重點實驗室,北京 100021;2國家癌癥中心/國家腫瘤臨床醫(yī)學研究中心/中國醫(yī)學科學院北京協(xié)和醫(yī)學院腫瘤醫(yī)院腫瘤結直腸外科,北京 100021;3國家癌癥中心/國家腫瘤臨床醫(yī)學研究中心/中國醫(yī)學科學院北京協(xié)和醫(yī)學院腫瘤醫(yī)院腫瘤特需醫(yī)療部,北京 100021

    根據(jù)2020年世界衛(wèi)生組織(WHO)發(fā)布的《2020年世界癌癥進展報告》,結直腸癌發(fā)病率在世界范圍內居第3位,病死率居第2位,是嚴重威脅人類健康的疾病[1]。為了攻克癌癥,新的治療技術層出不窮,其中免疫治療展現(xiàn)出不俗的治療效果,于2013年被Science雜志評為該年度最重要的科學突破。為了使免疫治療這一手段讓更多的患者受益,針對腫瘤微環(huán)境中的免疫細胞與腫瘤細胞的相互調控及其影響腫瘤細胞行為的研究成為熱點[2]。有研究顯示,腫瘤微環(huán)境中的免疫細胞具有多種協(xié)助腫瘤發(fā)生發(fā)展的生物學行為,包括參與腫瘤局部血管生成、協(xié)助腫瘤細胞發(fā)生轉移及免疫逃逸等[3-6]。但目前在探究免疫細胞與腫瘤細胞相關調控行為的研究中,對免疫細胞表型的劃分較為籠統(tǒng),雖然基于單細胞測序及分析技術的提升,已經(jīng)能在結直腸癌免疫微環(huán)境中識別更細分的免疫亞群,但這些免疫亞群的分布及其在大樣本臨床隊列中的潛在價值尚不清楚。CIBERSORTx算法是一種機器學習方法,可用于估計細胞比例及推斷細胞類型特異性基因表達譜[7-8]。本研究以單細胞轉錄組測序(scRNAseq)結果為參考,利用CIBERSORTx算法構建特征表達譜,測算癌癥基因組圖譜(TCGA)結腸癌隊列中腫瘤的免疫細胞浸潤情況,并評價其對腫瘤患者預后的預測價值。

    1 資料與方法

    1.1資料數(shù)據(jù)來源 通過基因表達數(shù)據(jù)庫(GEO)獲得GSE146771[9]中43 817個單細胞轉錄組,分析CD45+細胞樣本的每千堿基記錄本(transcripts per kilobase million,TPM)數(shù)據(jù),調整樣本特征(nFeatures)數(shù)值,進一步篩選獲得更為勻質的13 193個細胞。通過10×測序平臺,獲得這些細胞的轉錄組數(shù)據(jù)。結腸癌單細胞轉錄組測序數(shù)據(jù)包含來自9例結腸癌患者的外周血、術中切除的正常腸黏膜組織、結腸腺瘤及腺癌組織的單細胞測序基因表達譜;從加州大學圣克魯茲分校癌癥功能基因組數(shù)據(jù)庫(UCSC XENA)獲取TCGA結直腸癌隊列,該隊列包含臨床結直腸癌樣本513例,收集其臨床資料如年齡、性別、TNM分期、隨訪時間及生存狀態(tài)等,并將每千個堿基的轉錄每百萬映射讀取片段(fragments per kilobase per million,F(xiàn)PKM)值的轉錄組測序(RNAseq)數(shù)據(jù)轉換為TPM數(shù)據(jù),以便與scRNA-seq組數(shù)據(jù)保持一致。

    1.2建立scRNA-seq特征矩陣 使用CIBERSORTx在線分析平臺(https://cibersortx.stanford.edu/),在不進行物理分離的情況下,分析細胞類型特異性基因表達譜。按照指令上傳單細胞表達矩陣,在默認參數(shù)保持不變的條件下,用CIBERSORTx計算獲得38個免疫細胞類型的簽名矩陣。

    1.3CIBERSORTx計算細胞組分 依據(jù)CIBERSORTx的提示,上傳TCGA及GEO,獲得結直腸癌表達譜數(shù)據(jù)集,選取之前以scRNA-seq為參考數(shù)據(jù)得到的簽名矩陣,選擇S模式進行批量校正。將排列數(shù)設置為500,其他參數(shù)保持默認值。運行CIBERSORTx后,獲得了每個樣本中38個腫瘤浸潤免疫細胞亞群的相對比例,將P<0.05的樣本納入進一步研究。

    1.4免疫細胞亞群與轉移的相關性分析 為了探究免疫細胞亞群與腫瘤轉移的關系,將513例TCGA結直腸癌隊列分為轉移組(216例)與未轉移組(297例),比較兩組免疫細胞亞群分數(shù)的差異,并分析免疫細胞亞群與淋巴結轉移數(shù)的相關性。

    1.5基于免疫細胞亞群的生存分析 為探究不同免疫細胞亞群對結腸癌預后的影響,使用Kaplan-Meier法繪制患者的生存曲線,用log-rankχ2檢驗分析患者總生存期與免疫細胞亞群細胞分數(shù)的關系。

    1.6免疫細胞亞群功能分析 采用功能富集分析數(shù)據(jù)庫(the database for annotation, visualization and integrated discovery,DAVID)以及單細胞轉錄組分析和注釋工具(single-cell transcriptome analysis of pathway activation signatures,scTPA)進行功能聚類。將CIBERSORTx分析獲得的CD4+T亞群的評分與其對應的表達譜的基因列表進行相關性分析,獲得的相關性列表剔除在先前評估細胞分數(shù)中使用的基因列表,在剔除后的列表中取正負相關性最高的前500個基因進行功能聚類。對非免疫反應通路進行比較,觀察不同CD4+T亞群對結腸癌細胞功能調控的傾向。在聚類分析結果中,除由于免疫細胞亞群差異導致的免疫功能相關通路的差異,還著重關注免疫細胞亞群差異對結腸癌腫瘤細胞生物功能的影響,關注其對腫瘤細胞周期、細胞生長、細胞增殖、能量代謝、氧化磷酸化、缺氧應激等方面的調控差異。

    1.7基于免疫細胞亞群的預后模型的構建 采用LASSO算法篩選免疫細胞亞群,并使用Cox回歸分析構建風險預后模型。

    1.8統(tǒng)計學處理 采用SPSS 23.0軟件進行統(tǒng)計學分析。正態(tài)分布的計量資料組間比較采用獨立樣本t檢驗。不符合正態(tài)分布的計量資料比較采用Wilcoxon秩和檢驗。P<0.05為差異有統(tǒng)計學意義。

    2 結 果

    2.138個免疫細胞亞群的簽名矩陣 取GSE146771檢測基因數(shù)目nFeatures_RNA值在800~1200的勻質細胞的測序數(shù)據(jù)作為后續(xù)分析的參考集,基于CIBERSORTx算法創(chuàng)建包括38個免疫細胞集群共3801個基因的特征矩陣。如圖1所示。

    2.2不同TNM分期的免疫細胞亞群的比例差異

    與TCGA結腸癌未轉移患者比較,轉移患者中hT04_CD4-TCF7分數(shù)上調(P<0.01),hT11_CD4-CTLA4、hT12_CD8-LEF1、hB01_PlasmaB-IgG、hM10_Macro-IL1B分數(shù)下調(P<0.05,圖2A),且上述hT04_CD4-TCF7、hT12_CD8-LEF1、hB01_PlasmaB-IgG、hM10_Macro-IL1B免疫細胞亞群與淋巴結轉移數(shù)存在明顯相關關系(P<0.05,圖2B)。

    2.3免疫細胞亞群與總生存期的相關性 如圖3所示,hT04_CD4-TCF7、hT03_CD4-GNLY、hT09_CD4-CXCL13、hT16_CD8-CD6、hT17_CD8-CD160、hT14_CD8-CX3CR1、hI02_NK-GZMK、hB03_FollicularB-IgD、hB01_PlasmaB-IgG、hB02_GALTB-IgA、hM01_Mast-TPSAB1、hM12_TAMC1QC、hM04_cDC1-BATF3、hM03_cDC2-CD1C、hM05_Mono-CD14與患者的總生存期明顯相關(P<0.05)。

    2.410個CD4+T細胞亞群對結腸癌細胞調控功能的異質性 CIBERSORTx分析共獲得10個CD4+T亞群的評分。其中CD4-GNLY、CD4-ANXA1、CD4-CXCR6對結腸癌腫瘤的調控作用優(yōu)先聚集在血管生成相關通路,CD4-IL23R優(yōu)先聚集在能量代謝相關通路,CD4-GZMK優(yōu)先聚集在缺氧應激相關通路,CD4-ANXA1、CD4-TCF7、CD4-CTLA4優(yōu)先聚集在細胞增殖相關通路,CD4-ANXA1、CD4-TCF7、CD4-CTLA4、CD4-CXCL13優(yōu)先聚集在細胞黏附及遷移相關通路(圖4A)。此外,CCR7、CXCR5兩個亞群富集結果顯示,這兩個亞群在促癌生物學過程中的作用較其他亞群不顯著。

    比較轉移組與非轉移組間存在明顯差異的免疫細胞亞群,CD4+T細胞亞群中T04_CD4-TCF7與轉移呈正相關,T11_CD4-CTLA4與轉移呈負相關。使用scTPA線上工具比較兩個亞群細胞內通路激活的差異,發(fā)現(xiàn)兩者在細胞黏附及免疫激活或抑制相關通路中的激活狀態(tài)存在明顯差異(P<0.05,圖4B)。生存相關性分析結果顯示,T04_CD4-TCF7與患者生存呈正相關關系,T09_CD4-CXCL13與患者生存呈負相關關系(P<0.05)。與CD4-TCF7相比,CD4-CXCL13在抗腫瘤及TCR相關通路中的激活上明顯上調(P<0.05,圖4C)。生存分析結果顯示,T17_CD8-CD160對TCR相關通路、Fc受體及免疫刺激細胞因子相關的通路激活明顯高于T16_CD8-CD6細胞亞群(P<0.05,圖4D)。

    2.5基于免疫細胞亞群反映患者免疫狀態(tài)的預后風險模型的構建 通過LASSO算法從38個免疫細胞亞群中選擇4類免疫細胞亞群CD160-CD8+T細胞、CX3CR1-CD8+T細胞、BATF3-cDC1細胞、IgG-PlasmaB細胞(圖5A)。其中CD160-CD8+T細胞(P=0.0064)、IgG-PlasmaB細胞與風險評分呈負相關,CX3CR1-CD8+T細胞(P=0.0098)、BATF3-cDC1細胞(P=0.011)與風險評分呈正相關(P<0.05)。采用Cox分析建立預后風險評分模型,根據(jù)整個TCGA結腸癌隊列獲得的風險評分,患者可分為高風險組與低風險組(圖5B)。Kaplan-Meier曲線表明,在訓練集、驗證集及整組中,高風險組患者的死亡風險較高(圖5C–E)。在免疫細胞亞群功能差異的分析中,CD160-CD8+T細胞比CD8+T細胞表現(xiàn)出更高的Fc受體(Fc receptors,F(xiàn)cR)激活水平及T細胞抗原受體(T cell receptor,TCR)激活水平(P<0.05)。高風險組患者的T細胞功能失調評分高于低風險組,差異有統(tǒng)計學意義(P<0.05,圖5F)。

    3 討 論

    免疫治療在結腸癌中的應用較少,尚處于較初級的探索階段。鑒別結腸癌的腫瘤微環(huán)境特征,特別是免疫細胞亞群及功能狀態(tài)的特點,以及可能為臨床精準實施免疫治療提供參考的可靠預測標志物,對于提高結腸癌患者的免疫治療效果,提出可行的輔助治療方案具有重要意義。近年來,單細胞測序技術的發(fā)展使研究者能在更高分辨水平研究腫瘤微環(huán)境中的各種細胞成分,其中包括免疫細胞亞群在基因表達及功能狀態(tài)上的異質性。一系列已經(jīng)發(fā)表的單細胞測序分析結果顯示,不同癌種的免疫細胞亞群譜存在差異[10-12],癌種間免疫治療效果的差異可能與其免疫細胞譜系的差異相關,目前相關研究主要通過不同的表面標記分子,結合免疫組化空間定位分析免疫細胞亞群的功能及預后價值[13-14]。免疫細胞調控腫瘤細胞功能及腫瘤發(fā)生發(fā)展的研究更多著眼于傳統(tǒng)的免疫細胞分群,而探討免疫微環(huán)境中免疫細胞異質性的單細胞測序分析多側重于免疫細胞本身的基因表達特性及功能特性。為進一步揭示結腸癌腫瘤微環(huán)境中各種免疫細胞亞群對結腸癌細胞功能調控的差異以及對腫瘤進展的影響,本研究利用CIBERSORTx算法,基于單細胞測序分析定義的高分辨結腸癌腫瘤免疫微環(huán)境免疫細胞譜對TCGA的結腸癌患者隊列進行研究,以期探究結腸癌微環(huán)境中特征免疫細胞譜對結腸癌腫瘤細胞的調控作用及其與腫瘤進程的關系,結果顯示,hT04_CD4-TCF7、hT11_CD4-CTLA4、hT12_CD8-LEF1、hB01_PlasmaB-IgG、hM10_Macro-IL1B細胞亞群與結腸癌淋巴結轉移相關,hT04_CD4-TCF7、hT03_CD4-GNLY、hT09_CD4-CXCL13、hT16_CD8-CD6、hT17_CD8-CD160、hT14_CD8-CX3CR1、hI02_NK-GZMK、hB03_FollicularBIgD、hB01_PlasmaB-IgG、hB02_GALTB-IgA、hM01_Mast-TPSAB1、hM12_TAM-C1QC、hM04_cDC1-BATF3、hM03_cDC2-CD1C、hM05_Mono-CD14細胞亞群與患者的總生存期相關。

    本研究還探索了不同CD4+T細胞亞群對結腸癌腫瘤細胞調控功能的差異,發(fā)現(xiàn)對于腫瘤血管生成、腫瘤細胞生長增殖、能量代謝、氧化磷酸化、缺氧應激等通路的調節(jié),不同CD4+T細胞亞群表現(xiàn)出不同的優(yōu)先調控方向,從而影響腫瘤的發(fā)展進程,導致患者在轉移和預后結局方面的不同。

    在轉移與預后分析中,本研究發(fā)現(xiàn),CD4+和CD8+T細胞的不同亞群對轉移及預后具有相反的相關關系,通過比較這些突出亞群的細胞內通路激活狀態(tài),發(fā)現(xiàn)這些細胞亞群間通路激活狀態(tài)的差異很好地解釋了其對腫瘤發(fā)生發(fā)展方向的影響。先前的分析只是單獨地看待這些免疫細胞亞群的功能差異,而經(jīng)過構建模型,可以看到這種顯著影響患者生存的腫瘤微環(huán)境差異主要來自于哪些免疫細胞亞群,從而可以從整體觀察免疫細胞亞群對患者疾病進程的影響。腫瘤部位基于免疫細胞亞群的免疫細胞構成的異質性及免疫狀態(tài),對于患者的生存具有較高的影響,該模型在訓練集、驗證集及整組中顯示出較好的基于免疫狀態(tài)的死亡風險預測能力,相較單純基于差異基因簇,基于免疫狀態(tài)預測患者的死亡風險具有更高的臨床參考價值,展現(xiàn)了患者檢測節(jié)點的免疫狀態(tài)。在這一基于免疫細胞亞群構建的預后模型中,筆者發(fā)現(xiàn),低風險組與高風險組中有4個免疫細胞亞群的構成存在差異,這可能表明某些免疫細胞亞群之間存在免疫應答過程中的相互作用,降低腫瘤微環(huán)境中BATF3-cDC1樹突細胞的浸潤可降低CX3CR1-CD8+T細胞的浸潤,或提升IgG-PlasmaB細胞的浸潤,進而提升CD160-CD8+T細胞的浸潤,改善CD8+T細胞群的功能活性,從而改善患者生存,或輔助免疫檢測點治療。

    綜上所述,結腸癌腫瘤微環(huán)境中的免疫細胞構成及免疫微環(huán)境狀態(tài)能夠通過對結腸癌細胞功能的差異調控影響結腸癌的進展,進而影響患者的預后。本研究基于單細胞測序數(shù)據(jù),使用CIBERSORTx算法在TCGA隊列中研究免疫細胞譜對結腸癌進展及預后的影響,為進一步明確各免疫細胞亞群對腫瘤微環(huán)境的影響提供了方向和線索。但本研究尚有一定的局限,如研究結果基于公共數(shù)據(jù)和算法分析,未能闡述各免疫細胞亞群之間的相互關系網(wǎng)絡,且需要進一步的細胞生物學實驗對結果進行驗證。

    猜你喜歡
    單細胞亞群結腸癌
    TB-IGRA、T淋巴細胞亞群與結核免疫的研究進展
    甲狀腺切除術后T淋巴細胞亞群的變化與術后感染的相關性
    人工智能助力微生物單細胞鑒定
    科學(2020年4期)2020-11-26 08:27:16
    MicroRNA-381的表達下降促進結腸癌的增殖與侵襲
    外周血T細胞亞群檢測在惡性腫瘤中的價值
    結腸癌切除術術后護理
    聯(lián)合SNaPshot和單倍型分析技術建立G6PD缺乏癥單細胞基因診斷體系
    單細胞測序技術研究進展
    疣狀胃炎與T淋巴細胞亞群的相關研究進展
    中西醫(yī)結合治療晚期結腸癌78例臨床觀察
    国内久久婷婷六月综合欲色啪| 在线观看舔阴道视频| 国产色视频综合| 国产精品日韩av在线免费观看 | 午夜福利在线观看吧| 国产精品一区二区三区四区久久 | 婷婷丁香在线五月| 露出奶头的视频| 国产av在哪里看| 九色亚洲精品在线播放| 久久午夜综合久久蜜桃| 这个男人来自地球电影免费观看| 精品一品国产午夜福利视频| 91九色精品人成在线观看| 精品无人区乱码1区二区| 激情视频va一区二区三区| 成人亚洲精品av一区二区| 在线观看午夜福利视频| 91成年电影在线观看| 97人妻精品一区二区三区麻豆 | bbb黄色大片| 99国产综合亚洲精品| 午夜福利欧美成人| 亚洲精品在线美女| 亚洲中文字幕日韩| 色综合婷婷激情| 国产成人啪精品午夜网站| 麻豆国产av国片精品| 欧美成人一区二区免费高清观看 | 亚洲欧美日韩无卡精品| 国产成人精品久久二区二区免费| 可以在线观看毛片的网站| 欧美久久黑人一区二区| 午夜久久久在线观看| 夜夜夜夜夜久久久久| 高清黄色对白视频在线免费看| 日韩欧美一区二区三区在线观看| 少妇粗大呻吟视频| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩亚洲国产一区二区在线观看| 99在线视频只有这里精品首页| 老汉色av国产亚洲站长工具| 成年版毛片免费区| 亚洲在线自拍视频| 精品无人区乱码1区二区| 最近最新中文字幕大全免费视频| 国产一卡二卡三卡精品| 一进一出抽搐gif免费好疼| 精品久久蜜臀av无| 国产成人精品久久二区二区91| 国产精品久久久久久精品电影 | 免费久久久久久久精品成人欧美视频| 香蕉丝袜av| 精品人妻在线不人妻| 妹子高潮喷水视频| 成人精品一区二区免费| 怎么达到女性高潮| 久久久久久亚洲精品国产蜜桃av| 亚洲性夜色夜夜综合| 亚洲视频免费观看视频| 一进一出抽搐gif免费好疼| 曰老女人黄片| 精品国产乱子伦一区二区三区| 久久久久久亚洲精品国产蜜桃av| 一级片免费观看大全| 一进一出好大好爽视频| 国产精品影院久久| 国产精品国产高清国产av| 国产99久久九九免费精品| 国产精品美女特级片免费视频播放器 | 人妻久久中文字幕网| 丝袜人妻中文字幕| 成人免费观看视频高清| 男人舔女人下体高潮全视频| 一边摸一边做爽爽视频免费| 热re99久久国产66热| 夜夜夜夜夜久久久久| 国产精品 国内视频| 久久精品人人爽人人爽视色| 一区福利在线观看| av天堂在线播放| 一级a爱视频在线免费观看| 伦理电影免费视频| 女人精品久久久久毛片| 久久精品成人免费网站| www.999成人在线观看| 老熟妇仑乱视频hdxx| 欧美日韩精品网址| 欧美最黄视频在线播放免费| 操美女的视频在线观看| 18禁裸乳无遮挡免费网站照片 | 女性被躁到高潮视频| av天堂久久9| 久久久久久亚洲精品国产蜜桃av| 又大又爽又粗| 欧美激情高清一区二区三区| 日本 av在线| 日韩欧美国产一区二区入口| 黑人巨大精品欧美一区二区mp4| 91大片在线观看| 在线观看日韩欧美| 一个人免费在线观看的高清视频| 国产精品久久久人人做人人爽| 国产精品免费视频内射| 中文字幕av电影在线播放| 亚洲精品国产一区二区精华液| 一级,二级,三级黄色视频| 天堂动漫精品| 日本 av在线| 男人舔女人的私密视频| 国产亚洲精品综合一区在线观看 | 人妻丰满熟妇av一区二区三区| 国产精华一区二区三区| 午夜免费观看网址| 国产真人三级小视频在线观看| 免费在线观看影片大全网站| 欧美性长视频在线观看| 脱女人内裤的视频| 美女扒开内裤让男人捅视频| 日本一区二区免费在线视频| 日日干狠狠操夜夜爽| 国产精品爽爽va在线观看网站 | 精品国内亚洲2022精品成人| 女人被狂操c到高潮| 在线观看日韩欧美| 波多野结衣一区麻豆| 视频在线观看一区二区三区| 成人亚洲精品一区在线观看| 欧美激情极品国产一区二区三区| 美女高潮喷水抽搐中文字幕| 少妇裸体淫交视频免费看高清 | 别揉我奶头~嗯~啊~动态视频| 大陆偷拍与自拍| 中文字幕人妻丝袜一区二区| 久久中文看片网| 在线永久观看黄色视频| 亚洲成人精品中文字幕电影| 少妇 在线观看| 天天一区二区日本电影三级 | 久热这里只有精品99| 巨乳人妻的诱惑在线观看| 亚洲少妇的诱惑av| 国产精品影院久久| 黄色 视频免费看| 老司机在亚洲福利影院| www日本在线高清视频| 香蕉久久夜色| 午夜日韩欧美国产| 热re99久久国产66热| 亚洲国产欧美网| 一级毛片女人18水好多| 最新美女视频免费是黄的| 亚洲免费av在线视频| 97超级碰碰碰精品色视频在线观看| 操出白浆在线播放| 免费搜索国产男女视频| 免费高清在线观看日韩| 久久久久久国产a免费观看| 亚洲激情在线av| 日本 av在线| 日韩 欧美 亚洲 中文字幕| 人人妻人人澡欧美一区二区 | 女性生殖器流出的白浆| 黄色视频,在线免费观看| 欧美激情高清一区二区三区| 99在线人妻在线中文字幕| 亚洲久久久国产精品| 国产一卡二卡三卡精品| 亚洲国产毛片av蜜桃av| 人人妻,人人澡人人爽秒播| 老汉色∧v一级毛片| 黄色 视频免费看| 免费在线观看影片大全网站| 亚洲欧美精品综合一区二区三区| www日本在线高清视频| 手机成人av网站| 亚洲av五月六月丁香网| 久久久久久久久中文| 咕卡用的链子| 国产亚洲精品一区二区www| 91在线观看av| 一区二区三区国产精品乱码| 亚洲中文日韩欧美视频| 怎么达到女性高潮| 午夜日韩欧美国产| 久久青草综合色| 中文字幕久久专区| 国产成人av激情在线播放| 日本三级黄在线观看| 国产成人啪精品午夜网站| 正在播放国产对白刺激| 最近最新中文字幕大全电影3 | 欧美激情极品国产一区二区三区| 免费av毛片视频| 男女下面进入的视频免费午夜 | 99久久99久久久精品蜜桃| 色综合婷婷激情| 香蕉久久夜色| 欧美乱码精品一区二区三区| 免费久久久久久久精品成人欧美视频| 很黄的视频免费| 久久久久久亚洲精品国产蜜桃av| 成熟少妇高潮喷水视频| 亚洲国产中文字幕在线视频| 黄色丝袜av网址大全| 制服人妻中文乱码| 久久影院123| 啦啦啦韩国在线观看视频| 日日夜夜操网爽| 亚洲一区二区三区不卡视频| 精品卡一卡二卡四卡免费| 久久 成人 亚洲| 欧美绝顶高潮抽搐喷水| 女性被躁到高潮视频| 欧美绝顶高潮抽搐喷水| 激情视频va一区二区三区| 亚洲电影在线观看av| 夜夜看夜夜爽夜夜摸| 国产一级毛片七仙女欲春2 | 亚洲av片天天在线观看| 一级毛片精品| 一本大道久久a久久精品| 欧美精品啪啪一区二区三区| 777久久人妻少妇嫩草av网站| 老汉色∧v一级毛片| 欧美亚洲日本最大视频资源| 精品国产超薄肉色丝袜足j| 12—13女人毛片做爰片一| 免费在线观看黄色视频的| 久久久精品国产亚洲av高清涩受| 免费在线观看亚洲国产| 99久久精品国产亚洲精品| 桃色一区二区三区在线观看| 国产精品久久电影中文字幕| 欧美久久黑人一区二区| 日韩大尺度精品在线看网址 | 国产不卡一卡二| 国产免费av片在线观看野外av| 露出奶头的视频| 少妇裸体淫交视频免费看高清 | 欧美中文日本在线观看视频| 亚洲av电影不卡..在线观看| 亚洲成国产人片在线观看| 欧美丝袜亚洲另类 | 成人特级黄色片久久久久久久| 精品一区二区三区视频在线观看免费| 亚洲电影在线观看av| 国产欧美日韩一区二区三区在线| 亚洲欧美日韩无卡精品| 国产xxxxx性猛交| 欧美午夜高清在线| 久久香蕉国产精品| 亚洲欧美一区二区三区黑人| 十八禁人妻一区二区| 很黄的视频免费| 国产激情久久老熟女| 成人18禁高潮啪啪吃奶动态图| 黑人欧美特级aaaaaa片| 国产亚洲精品av在线| 中文亚洲av片在线观看爽| 成人国产综合亚洲| 女人被躁到高潮嗷嗷叫费观| 久久国产精品男人的天堂亚洲| 免费少妇av软件| 免费在线观看完整版高清| 精品免费久久久久久久清纯| 亚洲七黄色美女视频| 色播亚洲综合网| 免费在线观看黄色视频的| 老司机午夜福利在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 久久 成人 亚洲| 成人国产综合亚洲| 亚洲精品国产区一区二| 最好的美女福利视频网| 黄色女人牲交| 国产精品久久久久久精品电影 | 人人妻,人人澡人人爽秒播| 久久亚洲精品不卡| 国产精品一区二区三区四区久久 | 黄网站色视频无遮挡免费观看| 天天躁夜夜躁狠狠躁躁| 男人舔女人下体高潮全视频| 亚洲男人天堂网一区| 国内毛片毛片毛片毛片毛片| 国产成人精品久久二区二区91| 一级作爱视频免费观看| 熟女少妇亚洲综合色aaa.| 91大片在线观看| 国产免费男女视频| 一区二区三区高清视频在线| 日本欧美视频一区| 99在线人妻在线中文字幕| 欧美色视频一区免费| 麻豆国产av国片精品| 国内精品久久久久久久电影| 国产成人av激情在线播放| 一边摸一边抽搐一进一出视频| 老熟妇仑乱视频hdxx| av欧美777| 叶爱在线成人免费视频播放| 麻豆av在线久日| av片东京热男人的天堂| 人妻丰满熟妇av一区二区三区| 亚洲欧美日韩高清在线视频| 在线十欧美十亚洲十日本专区| 丝袜在线中文字幕| 在线观看一区二区三区| 国产三级黄色录像| 91av网站免费观看| 香蕉国产在线看| 欧美中文日本在线观看视频| 一区二区三区国产精品乱码| 9热在线视频观看99| 国产精品影院久久| 一进一出好大好爽视频| 免费观看人在逋| 一边摸一边抽搐一进一小说| 很黄的视频免费| 久久久久久久久久久丰满 | 色综合色国产| av中文乱码字幕在线| 无遮挡黄片免费观看| 亚洲av免费在线观看| 亚洲精品粉嫩美女一区| www.www免费av| 成人国产麻豆网| 久久精品久久久久久噜噜老黄 | 日本与韩国留学比较| 亚洲综合色惰| 亚洲七黄色美女视频| 日本撒尿小便嘘嘘汇集6| 成人一区二区视频在线观看| 亚洲专区国产一区二区| www日本黄色视频网| 日韩欧美免费精品| 精品无人区乱码1区二区| 成人综合一区亚洲| 日本成人三级电影网站| 亚洲精品成人久久久久久| a级毛片a级免费在线| 美女cb高潮喷水在线观看| 成人亚洲精品av一区二区| 一个人看视频在线观看www免费| 亚洲 国产 在线| 极品教师在线免费播放| 欧美3d第一页| 天堂网av新在线| 国产精品乱码一区二三区的特点| 国产乱人视频| 床上黄色一级片| 高清毛片免费观看视频网站| 国产久久久一区二区三区| 国产亚洲精品av在线| 国产极品精品免费视频能看的| 男女下面进入的视频免费午夜| 免费大片18禁| 亚洲欧美激情综合另类| 中文在线观看免费www的网站| 午夜免费男女啪啪视频观看 | 白带黄色成豆腐渣| 亚洲成人精品中文字幕电影| a级毛片a级免费在线| 国产美女午夜福利| 又黄又爽又刺激的免费视频.| 小说图片视频综合网站| 日日撸夜夜添| 一边摸一边抽搐一进一小说| 男人狂女人下面高潮的视频| 搡老熟女国产l中国老女人| 亚洲真实伦在线观看| 一级毛片久久久久久久久女| 国产视频内射| 欧美人与善性xxx| 日本与韩国留学比较| 嫩草影院精品99| 非洲黑人性xxxx精品又粗又长| 国产精品人妻久久久久久| av在线观看视频网站免费| 亚洲性夜色夜夜综合| 美女高潮的动态| 伊人久久精品亚洲午夜| 国产一区二区在线观看日韩| 狠狠狠狠99中文字幕| 99久国产av精品| 老司机深夜福利视频在线观看| 五月伊人婷婷丁香| 男女之事视频高清在线观看| 久久精品国产99精品国产亚洲性色| 搡老熟女国产l中国老女人| 一级a爱片免费观看的视频| 少妇猛男粗大的猛烈进出视频 | 久久国内精品自在自线图片| 欧美日韩亚洲国产一区二区在线观看| 中文字幕高清在线视频| 国产爱豆传媒在线观看| 丰满乱子伦码专区| 色视频www国产| 国产亚洲精品综合一区在线观看| 亚洲真实伦在线观看| 国产一区二区三区在线臀色熟女| 热99re8久久精品国产| 内射极品少妇av片p| 村上凉子中文字幕在线| 两个人视频免费观看高清| 如何舔出高潮| 欧美激情国产日韩精品一区| 少妇人妻一区二区三区视频| 97热精品久久久久久| 亚洲,欧美,日韩| 国产伦精品一区二区三区视频9| 一进一出抽搐动态| 中文字幕av在线有码专区| 亚洲av中文av极速乱 | 欧美中文日本在线观看视频| 国产探花极品一区二区| 美女大奶头视频| 夜夜爽天天搞| 国内精品久久久久精免费| 色哟哟哟哟哟哟| 女生性感内裤真人,穿戴方法视频| 51国产日韩欧美| 色综合婷婷激情| av国产免费在线观看| 日韩av在线大香蕉| 日本在线视频免费播放| 国产高清不卡午夜福利| 深夜精品福利| 日本一二三区视频观看| av专区在线播放| 亚洲欧美日韩无卡精品| 国产精品久久久久久av不卡| 精品久久久噜噜| 一区二区三区高清视频在线| 国产精品女同一区二区软件 | 永久网站在线| 国产一区二区三区在线臀色熟女| 久久精品国产亚洲网站| 久久精品国产鲁丝片午夜精品 | av天堂中文字幕网| 国产精品人妻久久久影院| 日本撒尿小便嘘嘘汇集6| .国产精品久久| 国产精品日韩av在线免费观看| 久久草成人影院| 亚洲第一区二区三区不卡| 热99在线观看视频| 国产欧美日韩一区二区精品| 91狼人影院| 国产精品,欧美在线| 欧美xxxx性猛交bbbb| 久久久久国产精品人妻aⅴ院| 亚洲一级一片aⅴ在线观看| 午夜激情欧美在线| 99久久九九国产精品国产免费| 国产欧美日韩精品亚洲av| 免费看av在线观看网站| 在线观看66精品国产| 久久精品综合一区二区三区| 国产精品久久电影中文字幕| 国产女主播在线喷水免费视频网站 | 久久亚洲真实| 亚洲在线观看片| 精品久久久久久,| 午夜福利视频1000在线观看| 麻豆一二三区av精品| 色噜噜av男人的天堂激情| 亚洲欧美日韩高清专用| 久久精品国产亚洲av涩爱 | 最新在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 哪里可以看免费的av片| 国产69精品久久久久777片| 欧美区成人在线视频| 久久久精品欧美日韩精品| 日本三级黄在线观看| 人妻夜夜爽99麻豆av| 成人午夜高清在线视频| 成人欧美大片| 成人国产麻豆网| bbb黄色大片| 免费av不卡在线播放| 久久精品综合一区二区三区| 九九爱精品视频在线观看| 国产高清视频在线播放一区| 观看美女的网站| 亚洲在线观看片| 黄片wwwwww| av女优亚洲男人天堂| 亚洲va日本ⅴa欧美va伊人久久| 成人国产麻豆网| 欧美xxxx性猛交bbbb| 免费无遮挡裸体视频| 丰满乱子伦码专区| 国产又黄又爽又无遮挡在线| 亚洲成人精品中文字幕电影| 成年人黄色毛片网站| 久久久久国产精品人妻aⅴ院| 久久精品影院6| 精品久久久久久成人av| 又粗又爽又猛毛片免费看| 一边摸一边抽搐一进一小说| 欧美高清成人免费视频www| 我的女老师完整版在线观看| 又爽又黄无遮挡网站| 国产女主播在线喷水免费视频网站 | 久久久久九九精品影院| 久久久久免费精品人妻一区二区| 精品乱码久久久久久99久播| 桃色一区二区三区在线观看| 亚洲第一区二区三区不卡| 偷拍熟女少妇极品色| 日韩大尺度精品在线看网址| 午夜福利欧美成人| 在线看三级毛片| 日本黄色视频三级网站网址| а√天堂www在线а√下载| 天堂av国产一区二区熟女人妻| 国产精品av视频在线免费观看| 美女 人体艺术 gogo| 在线天堂最新版资源| 日本爱情动作片www.在线观看 | 免费黄网站久久成人精品| 日韩欧美国产一区二区入口| 嫁个100分男人电影在线观看| 不卡视频在线观看欧美| 麻豆国产97在线/欧美| 一本精品99久久精品77| 日本色播在线视频| 国内精品久久久久久久电影| 国产一级毛片七仙女欲春2| 亚洲七黄色美女视频| 韩国av一区二区三区四区| 变态另类丝袜制服| 级片在线观看| 免费看a级黄色片| 午夜爱爱视频在线播放| 熟女电影av网| 嫩草影院入口| 天天一区二区日本电影三级| 亚洲av日韩精品久久久久久密| 久久久久久久久久久丰满 | 国产蜜桃级精品一区二区三区| 美女 人体艺术 gogo| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区免费毛片| 网址你懂的国产日韩在线| 色精品久久人妻99蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久视频播放| 免费无遮挡裸体视频| av视频在线观看入口| 特级一级黄色大片| 欧美性猛交╳xxx乱大交人| 少妇的逼水好多| 国产一区二区亚洲精品在线观看| 啦啦啦韩国在线观看视频| 在线免费观看的www视频| 级片在线观看| 精品乱码久久久久久99久播| 乱码一卡2卡4卡精品| 国产私拍福利视频在线观看| 高清在线国产一区| 男人舔女人下体高潮全视频| 亚洲男人的天堂狠狠| 一a级毛片在线观看| 一级av片app| 2021天堂中文幕一二区在线观| 中文字幕精品亚洲无线码一区| 两个人的视频大全免费| 国产精华一区二区三区| 国产极品精品免费视频能看的| 最后的刺客免费高清国语| 免费av不卡在线播放| 欧美黑人巨大hd| 老司机深夜福利视频在线观看| 久久久精品欧美日韩精品| 午夜视频国产福利| 久久久久久国产a免费观看| 联通29元200g的流量卡| 欧美精品国产亚洲| 黄色欧美视频在线观看| 欧美黑人欧美精品刺激| 直男gayav资源| 久久精品国产亚洲av香蕉五月| 精品久久久久久成人av| 日本一本二区三区精品| 一级a爱片免费观看的视频| 啦啦啦啦在线视频资源| 无遮挡黄片免费观看| 日本黄色片子视频| 亚洲精品一区av在线观看| 国产中年淑女户外野战色| 最好的美女福利视频网| 免费看av在线观看网站| 黄片wwwwww| 亚洲男人的天堂狠狠| 久久精品国产自在天天线| 国产单亲对白刺激| 成人三级黄色视频| 女人被狂操c到高潮| 精品人妻视频免费看| 国产69精品久久久久777片| 无遮挡黄片免费观看| 又爽又黄无遮挡网站| 少妇人妻一区二区三区视频| 桃色一区二区三区在线观看| 色在线成人网| 亚洲精品456在线播放app | 少妇猛男粗大的猛烈进出视频 | 国产精品98久久久久久宅男小说| 夜夜爽天天搞| 一级黄色大片毛片| 久久九九热精品免费| 久久久久久九九精品二区国产|