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

    基于條件關(guān)聯(lián)互補(bǔ)基因的乳腺癌預(yù)后分析

    2020-09-10 12:29:38張屹閆雙雙張成王冠方黃海云韓珊珊
    關(guān)鍵詞:乳腺癌

    張屹 閆雙雙 張成 王冠方 黃海云 韓珊珊

    摘 要:為了提高乳腺癌患者的生存率,改善病人的臨床治療效果,從分子機(jī)制上研究了乳腺癌的致病基因。首先對(duì)113個(gè)正常組織和1 109個(gè)癌癥組織的表達(dá)量進(jìn)行差異分析,然后對(duì)差異表達(dá)的基因采用條件聯(lián)合分析方式對(duì)互補(bǔ)基因進(jìn)行分組,并用逐步Cox回歸挑選出一組基因擬合預(yù)后模型。研究結(jié)果顯示:VWCE,SPDYC,CRYBG3,DEFB1,SEL1L2,NMNAT2 6個(gè)基因?qū)颊呱媛适怯泻Φ?,AMZ1,GJB2,CXCL2,ALDOC 4個(gè)基因?qū)颊呱媛适怯欣?,最終確定10個(gè)基因的預(yù)后模型能夠顯著地將樣本分為高風(fēng)險(xiǎn)組和低風(fēng)險(xiǎn)組,并且對(duì)乳腺癌患者5年和10年的生存率進(jìn)行了預(yù)測(cè),依賴時(shí)間的AUC值均可達(dá)0.7以上。所提方法能夠利用基因與基因之間的關(guān)聯(lián)性,很好地對(duì)高維數(shù)據(jù)進(jìn)行降維,消除基因與基因之間的共線性問(wèn)題,10個(gè)基因的預(yù)后模型可以對(duì)患者的臨床預(yù)測(cè)提供幫助。

    關(guān)鍵詞:生物數(shù)學(xué);乳腺癌;條件關(guān)聯(lián)基因;預(yù)后模型;臨床預(yù)測(cè)

    中圖分類(lèi)號(hào):Q786?文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1008-1542(2020)04-0349-07

    doi:10.7535/hbkd.2020yx04008

    乳腺癌是女性常見(jiàn)的一種癌癥,對(duì)女性身心健康有非常嚴(yán)重的影響。世界癌癥研究基金會(huì)數(shù)據(jù)顯示,2018年,全球有超過(guò)200萬(wàn)新增乳腺癌患者[1]。其中有研究表明被診斷為Ⅰ期和Ⅱ期(腫瘤比較小,或只見(jiàn)局部擴(kuò)散)的乳腺癌患者,5年的生存率大約為80%以上;但對(duì)于Ⅲ和Ⅳ期的晚期乳腺癌患者,5年的存活率僅有20%左右[2-3]??梢?jiàn),對(duì)乳腺癌患者的治療效果還有待提高。手術(shù)和藥物治療對(duì)乳腺癌患者的治療有所改善,但是近一半的乳腺癌患者仍會(huì)復(fù)發(fā),并因復(fù)發(fā)而導(dǎo)致死亡[4-5]。傳統(tǒng)治療是通過(guò)腫瘤大小、分期等風(fēng)險(xiǎn)因素作為乳腺癌患者的風(fēng)險(xiǎn)評(píng)估和確定治療計(jì)劃的依據(jù)。但是這些臨床病理風(fēng)險(xiǎn)因素?zé)o法將高風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)的乳腺癌患者區(qū)分開(kāi),亦無(wú)法預(yù)測(cè)哪類(lèi)患者更能受益于化療的治療方案。

    隨著高通量技術(shù)的不斷提高,開(kāi)始采用測(cè)序和基因芯片的數(shù)據(jù)來(lái)查找致癌基因或確立疾病的生物標(biāo)志物。這種數(shù)據(jù)的格式為在n個(gè)樣本中測(cè)得m個(gè)基因的基因表達(dá)量,結(jié)果應(yīng)為n×m維的矩陣,xij表示第i個(gè)樣本中測(cè)的第j個(gè)基因的值。針對(duì)同一基因在n個(gè)樣本中的測(cè)量數(shù)據(jù)所組成的向量稱為該基因的基因表達(dá)向量,可以表示為xi=[x1i,x2i,…,xni]?;虮磉_(dá)數(shù)據(jù)最明顯的特性就是樣本少而基因(特征)多,這種高維度的基因表達(dá)數(shù)據(jù)增加了鑒定出與癌癥相關(guān)基因的挑戰(zhàn)難度,嚴(yán)重影響了對(duì)疾病的正確預(yù)測(cè)。而且基因與基因之間也有著復(fù)雜的聯(lián)系,不同基因之間的關(guān)聯(lián)有可能會(huì)影響病人的生存期。國(guó)內(nèi)外對(duì)如何從眾多基因中挑選出與癌癥預(yù)后相關(guān)的基因一直都有研究。LIU等[6]對(duì)差異表達(dá)的基因,先使用單變量Cox比例回歸模型選出與預(yù)后顯著相關(guān)的,再使用逐步多元Cox回歸確立了最終的7個(gè)基因的預(yù)后模型。BIERMANN等[7]在單變量Cox回歸顯著基因的基礎(chǔ)上,將貝葉斯迭代模型平均應(yīng)用于多變量Cox回歸模型,確定了最終18個(gè)基因的預(yù)后模型。SHI等[8]通過(guò)構(gòu)建基因共表達(dá)網(wǎng)絡(luò)和挖掘與生存相關(guān)的模塊來(lái)鑒定可作為乳腺癌預(yù)后標(biāo)志物的基因。關(guān)于最后的預(yù)后模型的建立,通常采用基于基本的統(tǒng)計(jì)方法Cox回歸來(lái)選擇最后的基因,采用一致性指數(shù)(C-Index)來(lái)評(píng)估模型擬合,用風(fēng)險(xiǎn)評(píng)分(線性預(yù)測(cè)因子)將不同風(fēng)險(xiǎn)組的患者分開(kāi),最后通過(guò)ROC曲線函數(shù)下的時(shí)間依賴區(qū)域AUC對(duì)臨床預(yù)測(cè)能力進(jìn)行驗(yàn)證[9]。

    目前已發(fā)現(xiàn)與乳腺癌相關(guān)的某些分子標(biāo)記,多個(gè)基因生物標(biāo)記物對(duì)于癌癥的預(yù)后比單個(gè)生物標(biāo)記物更為準(zhǔn)確[10-13],但當(dāng)前臨床預(yù)測(cè)和治療的結(jié)果仍不能令人滿意。主要的難點(diǎn)在于發(fā)現(xiàn)基因與基因之間復(fù)雜的聯(lián)系并有效地減少基因個(gè)數(shù),為建立良好的預(yù)后模型奠定基礎(chǔ)。因此,本文借鑒了YANG 等[14]研究SNP對(duì)復(fù)雜性狀影響中采用的條件聯(lián)合多SNP分析的部分方法,先尋找互補(bǔ)基因,再對(duì)互補(bǔ)基因進(jìn)行Cox回歸,選擇與乳腺癌預(yù)后相關(guān)的基因,確定模型能夠很好地將樣本分為高風(fēng)險(xiǎn)組和低風(fēng)險(xiǎn)組,通過(guò)了對(duì)秩數(shù)的檢驗(yàn),并對(duì)乳腺

    癌患者5年和10年的生存率進(jìn)行了預(yù)測(cè),檢驗(yàn)了此預(yù)后模型的準(zhǔn)確性。這種從分子機(jī)制上確定疾病的分子標(biāo)志物對(duì)病人治療的臨床結(jié)果有了很大改善,研發(fā)分子診斷標(biāo)記物來(lái)監(jiān)測(cè)和預(yù)測(cè)各種癌癥的療效是有實(shí)際意義的[15]。

    通過(guò)確立出的生物標(biāo)志物能夠很好地將不同風(fēng)險(xiǎn)的個(gè)體區(qū)分開(kāi),采取不同的治療方案,改善治療效果,減少病痛,延長(zhǎng)患者壽命。

    1?材料及方法

    針對(duì)基因表達(dá)量的數(shù)據(jù)特點(diǎn),提出了基于條件篩選互補(bǔ)基因并進(jìn)行Cox回歸擬合的預(yù)后模型,對(duì)確定的預(yù)后模型進(jìn)行測(cè)試驗(yàn)證,算法過(guò)程如圖1所示。

    1.1?數(shù)據(jù)的下載與預(yù)處理

    從TCGA(the cancer genome atlas)數(shù)據(jù)庫(kù)中下載乳腺癌的基因表達(dá)量數(shù)據(jù)和對(duì)應(yīng)的臨床數(shù)據(jù),由于乳腺癌的基因表達(dá)量數(shù)據(jù)包含1 222個(gè)樣本(其中正常組織樣本113個(gè),癌癥組織樣本1 109個(gè)),樣本量較大,不能從網(wǎng)頁(yè)直接下載基因表達(dá)量的數(shù)據(jù),需下載官方下載工具gdc-client并在終端命令行中下載。將下載后的每個(gè)樣本中的各個(gè)基因的表達(dá)量數(shù)據(jù)進(jìn)行合并,合并后的基因表達(dá)量數(shù)據(jù)為矩陣形式,行代表的是1 222個(gè)樣本,列代表的是5萬(wàn)多個(gè)基因,基因表達(dá)量數(shù)據(jù)有3種格式,本文中用的是counts值格式的數(shù)據(jù)。對(duì)ENSEMBL格式的基因名字進(jìn)行轉(zhuǎn)換,轉(zhuǎn)換成只對(duì)人類(lèi)基因命名的Official Gene Symbol格式,并對(duì)樣本中80%以上表達(dá)量為0的低質(zhì)量基因進(jìn)行剔除,減少基因個(gè)數(shù)。整理好基因表達(dá)量數(shù)據(jù)后,對(duì)113個(gè)正常組織和1 109個(gè)癌癥組織的表達(dá)量進(jìn)行差異分析。首先對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化,然后創(chuàng)建分組矩陣和設(shè)計(jì)矩陣,查看上調(diào)基因和下調(diào)基因,篩選出與癌癥顯著相關(guān)的基因,降低維數(shù),以利于下一步模型的計(jì)算與構(gòu)建。其中差異表達(dá)基因的篩選標(biāo)準(zhǔn)為log FC>2和FDR<0.05,將滿足條件的基因與臨床數(shù)據(jù)合并,其中的臨床數(shù)據(jù)包括生存時(shí)間、生存狀態(tài)、年齡、組織學(xué)亞型等。

    1.2?尋找多組關(guān)聯(lián)互補(bǔ)基因

    對(duì)已有的高維基因表達(dá)量數(shù)據(jù)進(jìn)行差異分析,雖然很大程度降低了維數(shù),但是并不能直接對(duì)此時(shí)的差異基因進(jìn)行預(yù)后模型估計(jì),此時(shí)的基因不能滿足Cox回歸模型的前提假設(shè),存在多重共線性問(wèn)題。因此提出了一種新的方式尋找多組互補(bǔ)關(guān)聯(lián)基因。

    假設(shè)有多個(gè)基因?qū)ι鏁r(shí)間有線性影響,采用式(1)來(lái)表示此模型:

    1.3?Cox回歸

    對(duì)篩選出來(lái)的每一組基因分別進(jìn)行多元Cox回歸時(shí),并不是一組里的每個(gè)基因都是有用的。因此,將數(shù)據(jù)分為訓(xùn)練集和測(cè)試集,在訓(xùn)練集中對(duì)每一組數(shù)據(jù)進(jìn)行逐步Cox回歸,逐步Cox回歸采取AIC信息準(zhǔn)則作為衡量模型擬合的優(yōu)劣,AIC越小,說(shuō)明模型擬合得越好,多個(gè)基因的Cox回歸模型如式(18)所示:

    其中h(t,X)表示的是在時(shí)刻t時(shí),m個(gè)基因影響下的危險(xiǎn)率,h0(t)稱為基準(zhǔn)危險(xiǎn)率。w1是自變量的偏回歸系數(shù),一般取RRj=exp(wj)為相對(duì)危險(xiǎn)度,通過(guò)相對(duì)危險(xiǎn)度來(lái)直觀地解釋Xj的取值是否對(duì)h(t,X)的取值產(chǎn)生影響。如果RRj>1,說(shuō)明第j個(gè)基因?qū)Π┌Y樣本有著危害影響;如果RRj<1,說(shuō)明第j個(gè)基因?qū)Π┌Y樣本有著保護(hù)影響;如果RRj=1,說(shuō)明第j個(gè)基因基本不會(huì)對(duì)癌癥樣本產(chǎn)生影響。因此,可以通過(guò)擬合好的模型提取對(duì)乳腺癌病人生存時(shí)間有著主要影響的基因。Cox比例風(fēng)險(xiǎn)模型一般采用一致性指數(shù)即C-Index來(lái)評(píng)價(jià)模型的預(yù)測(cè)能力,一致性指數(shù)越接近于1,說(shuō)明模型的預(yù)測(cè)能力越強(qiáng),預(yù)測(cè)的準(zhǔn)確率越高。

    1.4?風(fēng)險(xiǎn)評(píng)分及預(yù)后分析

    對(duì)訓(xùn)練集和測(cè)試集的基因表達(dá)數(shù)據(jù)利用確定好的多個(gè)基因的Cox回歸模型對(duì)風(fēng)險(xiǎn)評(píng)分進(jìn)行預(yù)測(cè),風(fēng)險(xiǎn)評(píng)分表達(dá)式如式(19)所示:

    預(yù)測(cè)出的風(fēng)險(xiǎn)評(píng)分的中位數(shù)為閾值,將每個(gè)樣本的風(fēng)險(xiǎn)評(píng)分與閾值進(jìn)行比較,大于閾值的樣本為高風(fēng)險(xiǎn)組,小于樣本閾值的為低風(fēng)險(xiǎn)組。對(duì)訓(xùn)練集、測(cè)試集和整個(gè)數(shù)據(jù)集分別繪制K-M生存曲線,并用對(duì)數(shù)秩檢驗(yàn),檢驗(yàn)兩組患者的風(fēng)險(xiǎn)是否有顯著區(qū)別,如果檢驗(yàn)的p值<0.05,說(shuō)明通過(guò)了顯著性檢驗(yàn),選出來(lái)的基因確定的Cox回歸預(yù)測(cè)模型能夠很好地預(yù)測(cè)患者的風(fēng)險(xiǎn)。然后再對(duì)患者5年和10年的臨床預(yù)測(cè)能力繪制ROC,并計(jì)算出AUC值,越接近于1,說(shuō)明臨床預(yù)測(cè)能力越強(qiáng),準(zhǔn)確率越高。

    2?實(shí)驗(yàn)結(jié)果及分析

    2.1?差異分析結(jié)果

    下載后的基因表達(dá)量完成數(shù)據(jù)的預(yù)處理后為1 222個(gè)樣本和21 487個(gè)基因,對(duì)113個(gè)正常組織和1 109個(gè)癌癥組織進(jìn)行差異分析,使用R語(yǔ)言的limma包,將閾值設(shè)定為log FC>2和FDR<0.05,對(duì)2個(gè)條件同時(shí)滿足的差異表達(dá)基因繪制火山圖,如圖2所示,差異基因中有278個(gè)上調(diào)基因和868個(gè)下調(diào)基因。

    2.2?預(yù)后基因的確定

    1)將差異分析后得到的顯著基因與臨床數(shù)據(jù)合并,分別進(jìn)行單基因Cox回歸和單基因的一元線性回歸分析,選出與乳腺癌患者總生存期之間顯著的基因,檢驗(yàn)每個(gè)基因回歸系數(shù)值,2個(gè)回歸中以p<0.05作為閾值,篩選顯著基因并進(jìn)行合并,其中線性回歸中顯著基因有211個(gè),Cox回歸中顯著基因有438個(gè),對(duì)2種單基因回歸分析中顯著基因取并集后的612個(gè)基因按關(guān)聯(lián)互補(bǔ)方式選擇開(kāi)始進(jìn)行分組。

    2)提取單基因的一元線性回歸系數(shù)及顯著性檢驗(yàn)的p值, p值最小的一個(gè)作為第一組的第一個(gè)入選基因,剩下的基因叫做備選組,計(jì)算在入選基因組的基礎(chǔ)上再?gòu)膫溥x組中選擇一個(gè)基因加入的條件閾值,如果有p<0.05,選擇最小的p值加入,以此類(lèi)推,加入下一個(gè)基因,直至沒(méi)有p值小于0.05,這樣就選出了第一組基因,共4個(gè)。

    3)重復(fù)步驟2),至所有的基因都分了組,共分47個(gè)組。通過(guò)驗(yàn)證證明了組中基因個(gè)數(shù)小于15,進(jìn)行擬合Cox回歸,C-Index均小于0.7,所以將組中基因個(gè)數(shù)小于15的組刪掉,剩下15個(gè)組。

    4)完成數(shù)據(jù)分組后,對(duì)935個(gè)樣本數(shù)據(jù)集采用計(jì)算機(jī)生成隨機(jī)數(shù)的方式劃分訓(xùn)練集735個(gè)樣本和測(cè)試集200個(gè)樣本,對(duì)每一組基因分別進(jìn)行逐步Cox回歸,最終選出了第8組的10個(gè)基因作為乳腺癌的預(yù)后基因。 通過(guò)這10個(gè)基因確定了乳腺癌的預(yù)后模型,Cox回歸結(jié)果如表1所示。

    最終多變量Cox回歸模型的10個(gè)基因中有6個(gè)基因(VWCE,SPDYC,CRYBG3,DEFB1,SEL1L2,NMNAT2)的風(fēng)險(xiǎn)比值均大于1,說(shuō)明對(duì)乳腺癌患者的生存率有著危害的影響,剩下的4個(gè)基因(AMZ1,GJB2,CXCL2,ALDOC)的風(fēng)險(xiǎn)比值均小于1,說(shuō)明對(duì)乳腺癌患者的生存率有著有利的影響。對(duì)10個(gè)基因的回歸系數(shù)進(jìn)行顯著性檢驗(yàn),8個(gè)p值均小于0.05,剩下2個(gè)p值小于0.1,由于對(duì)多變量模型的預(yù)測(cè)能力強(qiáng)而被保留在模型中。

    2.3?預(yù)后模型的測(cè)試

    通過(guò)10個(gè)基因在訓(xùn)練集中估算乳腺癌患者的生存風(fēng)險(xiǎn)評(píng)分:

    最終的計(jì)算結(jié)果表明,風(fēng)險(xiǎn)評(píng)分的中位數(shù)為1.013,平均數(shù)為1.473,最小值為0.052,最大值為26.788。采用訓(xùn)練集中風(fēng)險(xiǎn)評(píng)分的中位數(shù)將訓(xùn)練集、測(cè)試集和整個(gè)數(shù)據(jù)集中的乳腺癌患者分為高風(fēng)險(xiǎn)組和低風(fēng)險(xiǎn)組,繪制生存風(fēng)險(xiǎn)曲線。訓(xùn)練集中高風(fēng)險(xiǎn)組和低風(fēng)險(xiǎn)組中樣本個(gè)數(shù)分別為367個(gè)和368個(gè),測(cè)試集中高風(fēng)險(xiǎn)組和低風(fēng)險(xiǎn)組中樣本個(gè)數(shù)分別為88個(gè)和112個(gè)。K-M分析顯示這2組的生存曲線顯著不同,與低風(fēng)險(xiǎn)組相比,高風(fēng)險(xiǎn)組的患者表現(xiàn)出更短的生存時(shí)間,且通過(guò)了對(duì)秩數(shù)檢驗(yàn),p值都分別小于0.05,如圖3所示。

    為了評(píng)估這10個(gè)基因預(yù)后模型的準(zhǔn)確率,分別對(duì)訓(xùn)練集、測(cè)試集和整個(gè)數(shù)據(jù)集繪制了時(shí)間依賴性的ROC曲線,

    如圖4所示,在訓(xùn)練集中,5年和10年的AUC值分別為 0.75和0.81。

    如圖5所示,在測(cè)試集中,5年和10年的AUC值分別為0.77和0.78。以上這些結(jié)果表明,風(fēng)險(xiǎn)評(píng)分是乳腺癌患者臨床結(jié)果的有力預(yù)測(cè)指標(biāo)。

    3?結(jié)?論

    研究并設(shè)計(jì)出能夠區(qū)分基因與基因之間復(fù)雜關(guān)系、更好地降低維數(shù)、更為簡(jiǎn)單高效的癌癥預(yù)后模型,具有十分重要的理論意義和臨床應(yīng)用價(jià)值。本研究使用了一種新的方式來(lái)尋找關(guān)聯(lián)互補(bǔ)基因,對(duì)每組基因分別進(jìn)行Cox回歸,得到了一組穩(wěn)健的預(yù)后模型,并分別通過(guò)了訓(xùn)練集和測(cè)試集的驗(yàn)證,通過(guò)10個(gè)預(yù)后基因,可以將高風(fēng)險(xiǎn)組和低風(fēng)險(xiǎn)組顯著分開(kāi)。

    通過(guò)在GeneCard上查找已經(jīng)選擇出來(lái)的基因的分子功能,發(fā)現(xiàn)VWCE是β-catenin信號(hào)傳導(dǎo)途徑中的調(diào)控元件,也是化學(xué)預(yù)防細(xì)胞癌的靶標(biāo)。SPDYC通過(guò)結(jié)合和激活CDK1與CDK2,能夠促進(jìn)整個(gè)細(xì)胞周期的進(jìn)展。GJB2表達(dá)于耳蝸,在表皮的上基底層和乳腺及子宮內(nèi)膜的上皮細(xì)胞中弱表達(dá)。防御素β1(DEFB1)基因與COPD相關(guān),在COPD患者中,人DEFB1的表達(dá)上調(diào)會(huì)影響肺功能的下降[16]。CXCL2是一種趨化因子,炎癥介質(zhì),趨化因子在協(xié)調(diào)免疫細(xì)胞向炎癥或損傷部位的協(xié)調(diào)募集中起作用[17]。比如其他趨化因子CCL2已被報(bào)道能夠負(fù)面調(diào)節(jié)管腔B乳腺癌細(xì)胞的過(guò)程,包括自噬和壞死[18]。有文獻(xiàn)證明CXCL2通過(guò)抑制ERK1/2信號(hào)傳導(dǎo)通路,減弱成骨細(xì)胞分化[19]。NMNAT2是維持健康軸突的重要生存因素[20]。

    總之,確定的10個(gè)乳腺癌的預(yù)后生物標(biāo)記物基因被證明具有高預(yù)測(cè)能力,且會(huì)隨著時(shí)間的推移保持穩(wěn)定研究結(jié)果有助于開(kāi)發(fā)更有效的預(yù)后工具,最終改善患者預(yù)后。得到的這些標(biāo)志性基因在其他文獻(xiàn)中雖然沒(méi)有顯示出與乳腺癌直接關(guān)聯(lián),但是這些基因與疾病相關(guān),說(shuō)明有致病影響,還需要作進(jìn)一步的實(shí)驗(yàn)研究,了解這些基因之間的相互作用,對(duì)臨床結(jié)果進(jìn)行預(yù)測(cè),制定治療方案,從而為高危乳腺癌患者提供更好的治療選擇,對(duì)低?;颊邷p少過(guò)度治療。此方法也可作為篩選其他癌癥預(yù)后標(biāo)志物的方法,結(jié)果有待進(jìn)一步測(cè)試驗(yàn)證。

    參考文獻(xiàn)/References:

    [1]?BRAY F, FERLAY J, SOERJOMATARAM I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA: A Cancer Journal for Clinicians, 2018, 68(6):394-424.

    [2]?CHEN L, LINDEN H M, ANDERSON B O, et al. Trends in 5-year survival rates among breast cancer patients by hormone receptor status and stage[J]. Breast Cancer Research and Treatment, 2014, 147(3):609-616.

    [3]?SAADATMAND S, BRETVELD R, SIESLING S, et al. Influence of tumour stage at breast cancer detection on survival in modern times: Population based study in 173,797 patients[J]. Nederlands Tijdschrift Voor Geneeskunde, 2016,160: A9800.

    [4]?BREWSTER A M, HORTOBAGYI G N, BROGLIO K R, et al. Residual risk of breast cancer recurrence 5 years after adjuvant therapy[J]. Journal of the National Cancer Institute, 2008, 100(16):1179-1183.

    [5]?ADES F, TRYFONIDIS K, ZARDAVAS D. The past and future of breast cancer treatment from the papyrus to individualised treatment approaches[J]. Ecancer Medical Science, 2017.doi:10.3332/ecancer.2017.746.

    [6]?LIU L, CHEN Z, SHI W, et al. Breast cancer survival prediction using seven prognostic biomarker genes[J]. Oncology Letters, 2019, 18(3):2907-2916.

    [7]?BIERMANN J, NEMES S, PARRIS T Z, et al. A novel 18-marker panel predicting clinical outcome in breast cancer[J].Cancer Epidemiology and Prevention Biomarkers, 2017, 26(11): 1619-1628.

    [8]?MALLETT S, ROYSTON P, WATERS R, et al. Reporting performance of prognostic models in cancer: A review[J]. BMC Medicine, 2010, 8(1):21.

    [9]?HEAGERTY P J, LUMLEY T, PEPE M S. Time-dependent ROC curves for censored survival data and a diagnostic marker[J]. Biometrics, 2000, 56(2): 337-344.

    [10]閆麗娜, 覃婷, 王彤. LASSO 方法在 Cox 回歸模型中的應(yīng)用[J]. 中國(guó)衛(wèi)生統(tǒng)計(jì), 2012, 29(1):58-60.

    YAN Lina, QIN Ting, WANG Tong. The application of LASSO in the Cox model[J]. Chinese Journal of Health Stats, 2012, 29(1):58-60.

    [11]線云開(kāi), 孫明立, 于兆進(jìn), 等. 基于TCGA數(shù)據(jù)庫(kù)篩選乳腺癌不良預(yù)后相關(guān)mi RNAs及風(fēng)險(xiǎn)評(píng)估[J]. 解剖科學(xué)進(jìn)展, 2019, 25(1):38-40.

    [12]孫景波, 陳嘉煒, 王植治, 等. NUF2 基因在乳腺癌中的表達(dá)及臨床意義[J]. 南方醫(yī)科大學(xué)學(xué)報(bào), 2019, 39(5): 591-597.

    SUN Jingbo, CHEN Jiawei, WANG Zhizhi, et al. Expression of NUF2 in breast cancer and its clinical significance[J]. Journal of Southern Medical University, 2019,39(5): 591-597.

    [13]SHI H J, ZHANG L, QU Y J, et al. Prognostic genes of breast cancer revealed by gene coexpression networkanalysis[J]. Oncology Letters, 2017,14(4): 4535-4542.

    [14]YANG J, FERREIRA T, MORRIS A P, et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits[J]. Nature Genetics, 2012, 44(4):369-375.

    [15]CH O, WILLIAM C S. Molecular diagnostics for monitoring and predicting therapeutic effect in cancer[J]. Expert Review of Molecular Diagnostics, 2011, 11(1):9-12.

    [16]ELLEN A, GNTHER G, BULLWINKEL J, et al. Increased expression of beta-defensin 1(DEFB1) in chronic obstructive pulmonary disease[J]. Plos One, 2011. doi:10.1371/journal.pone.0021898.

    [17]BORO M, BALAJI K N. CXCL1 and CXCL2 regulate NLRP3 inflammasome activation via G-protein-coupled receptor CXCR2[J]. The Journal of Immunology, 2017, 199(5):1660-1671.

    [18]FANG W B, YAO M, JOKAR I, et al. The CCL2 chemokine is a negative regulator of autophagy and necrosis in luminal B breast cancer cells[J]. Breast Cancer Research and Treatment, 2015, 150(2):309-320.

    [19]YANG Y, ZHOU X Y, LI Y J, et al. CXCL2 attenuates osteoblast differentiation by inhibiting the ERK1/2 signaling pathway[J]. Journal of Cell Science, 2019, 132(16): jcs230490.

    [20]GILLEY J, COLEMAN M P, BARRES B A. Endogenous Nmnat2 is an essential survival factor for maintenance of healthy axons[J]. Plos Biology, 2010, 8(1): e1000300.

    猜你喜歡
    乳腺癌
    絕經(jīng)了,是否就離乳腺癌越來(lái)越遠(yuǎn)呢?
    中老年保健(2022年6期)2022-08-19 01:41:48
    中醫(yī)治療乳腺癌的研究進(jìn)展
    乳腺癌的認(rèn)知及保健
    甘肅科技(2020年20期)2020-04-13 00:30:42
    乳腺癌是吃出來(lái)的嗎
    胸大更容易得乳腺癌嗎
    男人也得乳腺癌
    防治乳腺癌吃什么:禽比獸好
    幸福家庭(2019年14期)2019-01-06 09:15:38
    別逗了,乳腺癌可不分男女老少!
    祝您健康(2018年5期)2018-05-16 17:10:16
    PI3K在復(fù)發(fā)乳腺癌中的表達(dá)及意義
    CD47與乳腺癌相關(guān)性的研究進(jìn)展
    av黄色大香蕉| 国产黄片视频在线免费观看| 97精品久久久久久久久久精品| 女人久久www免费人成看片| 亚洲不卡免费看| a级一级毛片免费在线观看| 久久精品国产亚洲av涩爱| 看免费成人av毛片| 黄色一级大片看看| 成人毛片a级毛片在线播放| 亚洲欧美日韩无卡精品| 最黄视频免费看| 国产视频首页在线观看| 高清欧美精品videossex| 韩国av在线不卡| 一个人免费看片子| 欧美高清成人免费视频www| 夫妻性生交免费视频一级片| 三级国产精品欧美在线观看| 国产精品一区二区性色av| 免费观看a级毛片全部| 亚洲伊人久久精品综合| 欧美精品国产亚洲| 国产亚洲最大av| 综合色丁香网| 免费观看性生交大片5| 欧美成人a在线观看| 大香蕉97超碰在线| 亚洲四区av| 热re99久久精品国产66热6| 九九爱精品视频在线观看| 又爽又黄a免费视频| 晚上一个人看的免费电影| 天美传媒精品一区二区| 九色成人免费人妻av| 久久久精品94久久精品| 亚洲国产欧美人成| 午夜免费鲁丝| 欧美xxxx黑人xx丫x性爽| 丝袜脚勾引网站| 国产v大片淫在线免费观看| 亚洲图色成人| 免费人成在线观看视频色| 人妻一区二区av| 亚洲最大成人中文| 成年美女黄网站色视频大全免费 | 精品久久久久久久末码| 亚洲综合精品二区| 亚洲av成人精品一区久久| 国产精品久久久久久久电影| 亚洲欧洲日产国产| 人妻少妇偷人精品九色| 人妻一区二区av| 91aial.com中文字幕在线观看| 国产精品一区二区在线不卡| 久久精品久久久久久久性| 熟女av电影| 免费高清在线观看视频在线观看| 91久久精品国产一区二区成人| 国产精品蜜桃在线观看| av福利片在线观看| 我的女老师完整版在线观看| 日韩亚洲欧美综合| 国产精品国产av在线观看| 亚洲精品视频女| 狂野欧美激情性xxxx在线观看| 少妇人妻精品综合一区二区| 成人高潮视频无遮挡免费网站| 26uuu在线亚洲综合色| 91精品国产国语对白视频| 国产91av在线免费观看| 欧美三级亚洲精品| 高清av免费在线| 欧美精品亚洲一区二区| 久久av网站| 少妇人妻一区二区三区视频| 精品久久久久久久末码| 精品亚洲成国产av| 欧美精品一区二区大全| 黄色怎么调成土黄色| 蜜桃亚洲精品一区二区三区| 久久女婷五月综合色啪小说| 日本猛色少妇xxxxx猛交久久| 精品人妻熟女av久视频| 熟女av电影| 欧美日韩视频精品一区| 欧美成人一区二区免费高清观看| 乱码一卡2卡4卡精品| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看| 高清欧美精品videossex| 直男gayav资源| 国产毛片在线视频| 亚洲欧美日韩无卡精品| 春色校园在线视频观看| 久久毛片免费看一区二区三区| 免费高清在线观看视频在线观看| 午夜福利视频精品| 黄色欧美视频在线观看| 国产真实伦视频高清在线观看| 亚洲精品久久午夜乱码| 观看免费一级毛片| 久久99精品国语久久久| 美女内射精品一级片tv| 国产欧美日韩精品一区二区| 少妇人妻一区二区三区视频| 国产亚洲精品久久久com| 久久这里有精品视频免费| 91午夜精品亚洲一区二区三区| 成人影院久久| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 国产成人91sexporn| 久久国内精品自在自线图片| 日日摸夜夜添夜夜添av毛片| 日本wwww免费看| av不卡在线播放| 国产日韩欧美亚洲二区| 蜜臀久久99精品久久宅男| 久久综合国产亚洲精品| 精品亚洲乱码少妇综合久久| 久久精品人妻少妇| 大片免费播放器 马上看| 免费黄频网站在线观看国产| 男女免费视频国产| 成人18禁高潮啪啪吃奶动态图 | 欧美一区二区亚洲| av在线app专区| 免费观看在线日韩| 黄色怎么调成土黄色| 亚洲av在线观看美女高潮| 熟女av电影| av国产免费在线观看| 亚洲婷婷狠狠爱综合网| 日韩一区二区视频免费看| 这个男人来自地球电影免费观看 | 久久久久人妻精品一区果冻| 国产日韩欧美亚洲二区| 国产精品福利在线免费观看| 毛片女人毛片| 欧美日本视频| 高清黄色对白视频在线免费看 | 一本久久精品| 日本欧美国产在线视频| 久久99热6这里只有精品| 国产成人a区在线观看| 少妇高潮的动态图| 亚洲美女黄色视频免费看| 麻豆国产97在线/欧美| 欧美xxxx黑人xx丫x性爽| 国产精品久久久久久久久免| 亚洲丝袜综合中文字幕| 插逼视频在线观看| 精品国产乱码久久久久久小说| 搡老乐熟女国产| 99久国产av精品国产电影| 九九爱精品视频在线观看| 狠狠精品人妻久久久久久综合| 自拍欧美九色日韩亚洲蝌蚪91 | 国产高清国产精品国产三级 | 永久网站在线| 最新中文字幕久久久久| 我要看日韩黄色一级片| 在线 av 中文字幕| 老师上课跳d突然被开到最大视频| 女人久久www免费人成看片| 成人黄色视频免费在线看| 精品久久国产蜜桃| 亚洲精品国产av成人精品| 国产乱人偷精品视频| 亚洲成人av在线免费| 日韩不卡一区二区三区视频在线| 黄色视频在线播放观看不卡| 国产黄片美女视频| 青春草亚洲视频在线观看| 亚洲精品一区蜜桃| 久久99蜜桃精品久久| 97热精品久久久久久| 男女无遮挡免费网站观看| 黄片无遮挡物在线观看| 成人无遮挡网站| 人妻 亚洲 视频| 哪个播放器可以免费观看大片| 精品久久久精品久久久| 中国美白少妇内射xxxbb| 九九久久精品国产亚洲av麻豆| 中文乱码字字幕精品一区二区三区| 亚洲性久久影院| 国产精品一及| 国产精品伦人一区二区| 午夜福利网站1000一区二区三区| 国产成人免费无遮挡视频| 国产成人精品一,二区| 国产精品无大码| 高清欧美精品videossex| 国产伦在线观看视频一区| 国产成人91sexporn| 少妇人妻精品综合一区二区| 精品熟女少妇av免费看| 午夜福利视频精品| 伦精品一区二区三区| 七月丁香在线播放| 99热这里只有是精品在线观看| .国产精品久久| av国产免费在线观看| 久久韩国三级中文字幕| 一级毛片aaaaaa免费看小| 欧美日韩一区二区视频在线观看视频在线| 成人影院久久| 男女无遮挡免费网站观看| 婷婷色综合大香蕉| 久久人人爽人人片av| 国产白丝娇喘喷水9色精品| 亚洲国产欧美人成| 精品亚洲成a人片在线观看 | 亚洲人成网站高清观看| 最新中文字幕久久久久| 欧美3d第一页| 看十八女毛片水多多多| 青青草视频在线视频观看| tube8黄色片| 精品人妻视频免费看| 99久久精品国产国产毛片| 国产探花极品一区二区| 免费av中文字幕在线| 免费黄网站久久成人精品| 我的老师免费观看完整版| 中文精品一卡2卡3卡4更新| 啦啦啦在线观看免费高清www| 在线观看一区二区三区| 精品少妇久久久久久888优播| 亚洲精品456在线播放app| 中文资源天堂在线| 日日啪夜夜撸| 美女xxoo啪啪120秒动态图| 干丝袜人妻中文字幕| 欧美人与善性xxx| 国产精品国产三级专区第一集| 夫妻午夜视频| 免费在线观看成人毛片| 伦精品一区二区三区| 男女下面进入的视频免费午夜| 欧美zozozo另类| 纯流量卡能插随身wifi吗| 大香蕉久久网| 26uuu在线亚洲综合色| 久久6这里有精品| 午夜精品国产一区二区电影| 成人亚洲精品一区在线观看 | 色哟哟·www| 韩国av在线不卡| 在线精品无人区一区二区三 | 国产精品蜜桃在线观看| 亚洲av二区三区四区| 777米奇影视久久| 女的被弄到高潮叫床怎么办| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美清纯卡通| 网址你懂的国产日韩在线| 久久久久人妻精品一区果冻| 亚洲av中文字字幕乱码综合| 九色成人免费人妻av| 伊人久久国产一区二区| 26uuu在线亚洲综合色| 99热全是精品| 久久久久久九九精品二区国产| 最后的刺客免费高清国语| 亚洲在久久综合| 亚洲内射少妇av| 又爽又黄a免费视频| 哪个播放器可以免费观看大片| 国产黄色视频一区二区在线观看| 国产免费福利视频在线观看| 夜夜骑夜夜射夜夜干| 日韩成人伦理影院| av国产免费在线观看| 美女福利国产在线 | 色综合色国产| av天堂中文字幕网| 草草在线视频免费看| 日韩av在线免费看完整版不卡| 日本午夜av视频| 毛片一级片免费看久久久久| 亚洲成人一二三区av| 少妇高潮的动态图| 丰满人妻一区二区三区视频av| 亚洲国产精品成人久久小说| 日韩在线高清观看一区二区三区| 成人黄色视频免费在线看| 一二三四中文在线观看免费高清| 久久久久久久久久久丰满| 成年美女黄网站色视频大全免费 | 国产伦在线观看视频一区| 亚洲欧洲国产日韩| 我要看黄色一级片免费的| 大又大粗又爽又黄少妇毛片口| 蜜臀久久99精品久久宅男| 男男h啪啪无遮挡| 成人一区二区视频在线观看| 各种免费的搞黄视频| 久久久久久人妻| 中文欧美无线码| 少妇人妻精品综合一区二区| 亚洲精品乱码久久久久久按摩| 我要看黄色一级片免费的| 波野结衣二区三区在线| 亚洲高清免费不卡视频| 国产成人精品福利久久| 少妇猛男粗大的猛烈进出视频| 欧美+日韩+精品| 97在线视频观看| 欧美97在线视频| 七月丁香在线播放| 精品亚洲成国产av| 日本vs欧美在线观看视频 | 黄色怎么调成土黄色| 成人美女网站在线观看视频| 亚洲一级一片aⅴ在线观看| 国产深夜福利视频在线观看| 草草在线视频免费看| 欧美极品一区二区三区四区| 久久久久精品久久久久真实原创| 国产精品久久久久久精品古装| 男女边摸边吃奶| 国产精品嫩草影院av在线观看| 插逼视频在线观看| 国产精品久久久久成人av| 国产白丝娇喘喷水9色精品| 久久99热这里只频精品6学生| 国产精品不卡视频一区二区| 搡老乐熟女国产| 亚洲精品国产av成人精品| 久久女婷五月综合色啪小说| 国模一区二区三区四区视频| 国产熟女欧美一区二区| 2021少妇久久久久久久久久久| 插逼视频在线观看| 亚洲va在线va天堂va国产| 国产男女超爽视频在线观看| 久久韩国三级中文字幕| av福利片在线观看| 日韩免费高清中文字幕av| 三级经典国产精品| 国精品久久久久久国模美| 亚洲无线观看免费| 人人妻人人添人人爽欧美一区卜 | 有码 亚洲区| 青春草视频在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 久久久a久久爽久久v久久| 高清日韩中文字幕在线| 最近最新中文字幕大全电影3| 久久久久久久久大av| 亚洲精品一区蜜桃| 97超视频在线观看视频| 丝袜喷水一区| 99久久中文字幕三级久久日本| av黄色大香蕉| 国产精品人妻久久久影院| 在线观看人妻少妇| 啦啦啦在线观看免费高清www| 免费少妇av软件| 日韩av不卡免费在线播放| 在线免费十八禁| 人妻夜夜爽99麻豆av| 在线 av 中文字幕| 国产91av在线免费观看| 精品国产乱码久久久久久小说| 亚洲国产成人一精品久久久| 最新中文字幕久久久久| 免费看日本二区| 一区在线观看完整版| 国产无遮挡羞羞视频在线观看| 久久99蜜桃精品久久| 妹子高潮喷水视频| 精品国产乱码久久久久久小说| 国产91av在线免费观看| 久久久久久伊人网av| 97在线人人人人妻| 一本色道久久久久久精品综合| tube8黄色片| 在线观看三级黄色| 精品久久久噜噜| 国产亚洲91精品色在线| 亚洲丝袜综合中文字幕| 成人18禁高潮啪啪吃奶动态图 | 蜜桃亚洲精品一区二区三区| 久久久久精品久久久久真实原创| 国产乱人偷精品视频| av免费在线看不卡| 精品视频人人做人人爽| 中国国产av一级| 国产精品99久久久久久久久| 亚洲色图综合在线观看| 国产av精品麻豆| 亚洲国产最新在线播放| 国产免费一级a男人的天堂| a 毛片基地| 男女边摸边吃奶| a级毛色黄片| 黄片无遮挡物在线观看| 亚洲人成网站在线播| 国产精品秋霞免费鲁丝片| 极品少妇高潮喷水抽搐| 国产亚洲欧美精品永久| 赤兔流量卡办理| 狂野欧美激情性xxxx在线观看| 久久av网站| 91在线精品国自产拍蜜月| 国产精品一区二区三区四区免费观看| 一区二区三区乱码不卡18| 国产免费福利视频在线观看| 秋霞在线观看毛片| 精品人妻视频免费看| 国产欧美日韩一区二区三区在线 | 日韩中字成人| 中文在线观看免费www的网站| av又黄又爽大尺度在线免费看| 精品久久久噜噜| 国国产精品蜜臀av免费| 高清视频免费观看一区二区| 毛片一级片免费看久久久久| 永久网站在线| 免费在线观看成人毛片| 精品国产乱码久久久久久小说| av线在线观看网站| 又粗又硬又长又爽又黄的视频| 性色av一级| 免费观看a级毛片全部| 亚洲美女搞黄在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 插逼视频在线观看| 高清不卡的av网站| 免费看不卡的av| 欧美一级a爱片免费观看看| 身体一侧抽搐| 免费观看无遮挡的男女| 欧美变态另类bdsm刘玥| 高清视频免费观看一区二区| 国产欧美亚洲国产| 各种免费的搞黄视频| 亚洲精品456在线播放app| 亚洲精品自拍成人| 亚洲精品日本国产第一区| 老司机影院毛片| 黄色怎么调成土黄色| 97超碰精品成人国产| 国产男女内射视频| 丝袜喷水一区| 亚洲国产精品专区欧美| 欧美极品一区二区三区四区| 亚洲欧美日韩无卡精品| 麻豆精品久久久久久蜜桃| 日韩av不卡免费在线播放| 国产亚洲精品久久久com| 国产又色又爽无遮挡免| 久久ye,这里只有精品| 欧美+日韩+精品| 美女内射精品一级片tv| 亚洲精品aⅴ在线观看| 久久99热6这里只有精品| 日韩精品有码人妻一区| tube8黄色片| 精品一区在线观看国产| 色综合色国产| 精品一区二区三区视频在线| 少妇人妻一区二区三区视频| 97在线人人人人妻| a级毛片免费高清观看在线播放| 亚洲久久久国产精品| av国产久精品久网站免费入址| 亚洲精品aⅴ在线观看| 热re99久久精品国产66热6| 国产精品秋霞免费鲁丝片| 日韩视频在线欧美| 午夜免费鲁丝| av.在线天堂| 最近2019中文字幕mv第一页| 久久精品国产自在天天线| 国产 一区 欧美 日韩| 丰满乱子伦码专区| 国产白丝娇喘喷水9色精品| 亚洲av免费高清在线观看| 男女啪啪激烈高潮av片| 特大巨黑吊av在线直播| xxx大片免费视频| 亚洲无线观看免费| 18禁裸乳无遮挡免费网站照片| 日韩av在线免费看完整版不卡| 欧美日韩亚洲高清精品| 亚洲精品乱久久久久久| 国产毛片在线视频| 亚洲综合精品二区| 久久99精品国语久久久| av在线app专区| 精品一区在线观看国产| 久久女婷五月综合色啪小说| 欧美精品亚洲一区二区| 又大又黄又爽视频免费| 黄色欧美视频在线观看| 亚洲精品久久午夜乱码| 永久网站在线| 亚洲精品456在线播放app| 国产精品国产三级国产av玫瑰| 欧美变态另类bdsm刘玥| 欧美成人a在线观看| 国产爱豆传媒在线观看| 99热这里只有精品一区| 少妇的逼好多水| 久久久久国产网址| 国产精品久久久久久精品电影小说 | 亚洲精品日韩在线中文字幕| 高清毛片免费看| 色哟哟·www| 国产成人精品福利久久| 国产人妻一区二区三区在| av免费观看日本| 夜夜爽夜夜爽视频| 国产av国产精品国产| 日本免费在线观看一区| 妹子高潮喷水视频| av又黄又爽大尺度在线免费看| 欧美人与善性xxx| 在线观看免费高清a一片| 五月开心婷婷网| 王馨瑶露胸无遮挡在线观看| 中文字幕av成人在线电影| 晚上一个人看的免费电影| 免费播放大片免费观看视频在线观看| 精品亚洲乱码少妇综合久久| 多毛熟女@视频| 女性生殖器流出的白浆| 狂野欧美白嫩少妇大欣赏| 欧美bdsm另类| 亚洲av中文字字幕乱码综合| 天天躁日日操中文字幕| 日韩av不卡免费在线播放| 国产久久久一区二区三区| 国国产精品蜜臀av免费| av天堂中文字幕网| 青春草国产在线视频| 我的女老师完整版在线观看| 欧美+日韩+精品| 国产精品嫩草影院av在线观看| www.av在线官网国产| 毛片女人毛片| 婷婷色综合大香蕉| 日韩一区二区视频免费看| 天天躁日日操中文字幕| 午夜福利在线在线| 免费黄频网站在线观看国产| 日韩成人av中文字幕在线观看| 永久网站在线| 纯流量卡能插随身wifi吗| 夫妻午夜视频| 大陆偷拍与自拍| 久久精品国产亚洲网站| 免费不卡的大黄色大毛片视频在线观看| 欧美+日韩+精品| 国产老妇伦熟女老妇高清| av国产免费在线观看| 秋霞伦理黄片| 亚洲国产最新在线播放| 欧美精品一区二区免费开放| 大又大粗又爽又黄少妇毛片口| 国精品久久久久久国模美| 中文字幕av成人在线电影| 尤物成人国产欧美一区二区三区| 国产精品一区二区在线不卡| 免费观看av网站的网址| 嘟嘟电影网在线观看| 亚洲国产欧美在线一区| 亚洲最大成人中文| 久久久久精品性色| 日本wwww免费看| 久久av网站| 精品一区二区三区视频在线| 日本wwww免费看| 亚洲性久久影院| 男人爽女人下面视频在线观看| 精品国产三级普通话版| 日韩成人伦理影院| 五月天丁香电影| 国产精品秋霞免费鲁丝片| 亚洲无线观看免费| 国产淫片久久久久久久久| 在线观看一区二区三区| 亚洲精品一二三| 毛片女人毛片| 亚洲熟女精品中文字幕| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久av不卡| freevideosex欧美| 国产日韩欧美在线精品| 91精品伊人久久大香线蕉| 一区二区三区精品91| 男人添女人高潮全过程视频| 亚洲婷婷狠狠爱综合网| 国产真实伦视频高清在线观看| 久久这里有精品视频免费| 一级毛片我不卡| 亚洲成人中文字幕在线播放| 亚洲人成网站在线观看播放| 久久国内精品自在自线图片| 国内揄拍国产精品人妻在线| 中文资源天堂在线| 下体分泌物呈黄色| 国产精品一及| 国产在线免费精品| 国产视频内射| 国产精品免费大片| 99视频精品全部免费 在线| 女的被弄到高潮叫床怎么办| 成人影院久久| 成人毛片60女人毛片免费|