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

    孟德爾隨機(jī)化在胰腺癌研究中的應(yīng)用現(xiàn)狀與展望

    2024-12-31 00:00:00杜凱豪侯立朝羅蘭明慧東小鴿蔣威王展
    臨床肝膽病雜志 2024年10期

    摘要:胰腺癌的發(fā)病隱匿性強(qiáng),治療困難,早期診斷和治療存在諸多局限性。本文總結(jié)了孟德爾隨機(jī)化(MR)在胰腺癌風(fēng)險(xiǎn)因素探索中的應(yīng)用進(jìn)展,特別分析了腸道微生態(tài)、生活方式、代謝性疾病等因素的因果關(guān)系。通過大規(guī)?;蚪M關(guān)聯(lián)研究(GWAS)數(shù)據(jù),MR分析揭示了若干與胰腺癌風(fēng)險(xiǎn)相關(guān)的生物標(biāo)志物。雙樣本MR是當(dāng)前研究中的常用方法,包括逆方差加權(quán)法、加權(quán)中位數(shù)、MR-Egger法等,有助于從遺傳角度解釋疾病的因果網(wǎng)絡(luò)。盡管MR策略可以為理解胰腺癌的病因?qū)W提供新視角,但其在數(shù)據(jù)合成、工具變量選擇及多效性評(píng)估方面仍需謹(jǐn)慎。新興分析模型如貝葉斯加權(quán)MR(BWMR)、CAUSE和多變量MR(MVMR)等,為綜合評(píng)估多重風(fēng)險(xiǎn)因素及其相互關(guān)系帶來新的可能。未來,結(jié)合上述方法及累積增多的遺傳流行病學(xué)資料,MR分析預(yù)期能夠?yàn)榘l(fā)現(xiàn)胰腺癌的潛在治療靶點(diǎn)和制訂預(yù)防策略提供更為堅(jiān)實(shí)的證據(jù)基礎(chǔ)。

    關(guān)鍵詞:胰腺腫瘤;孟德爾隨機(jī)化分析;因果律

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(82160131);青海省科技廳項(xiàng)目(2021-ZJ-963Q)

    Current status and prospects of the application of Mendelian randomization in pancreatic cancer research

    DU Kaihao1,HOU Lizhao2a,LUO Lanminghui1,DONG Xiaoge1,JIANG Wei1,WANG Zhan2b.(1.Clinical Medical College,Qinghai University,Xining 810000,China;2.a.Department of General Surgery,b.Department of Medical Engineering Integration and Translational Application,The Affiliated Hospital of Qinghai University,Xining 810000,China)

    Corresponding author:WANG Zhan,ufofu01@163.com(ORCID:0000-0002-4225-8136)

    Abstract:Pancreatic cancer often has an insidious onset and difficulties in treatment,with various limitations in early diagnosis and treatment.This article reviews the application of Mendelian randomization(MR)in exploring the risk factors for pancreatic cancer,with a special focus on the causal relationships of factors such as gut microbiota,lifestyle,and metabolic diseases.Leveraging data from large-scale genome-wide association studies(GWAS),MR analysis has revealed several biomarkers associated with the risk of pancreatic cancer.The two-sample MR approach is commonly used in current research,including the methods such as Inverse Variance Weighted,Weighted Median,and MR-Egger,which helps to explain the causal network of the disease from a genetic perspective.While MR strategy provides a new perspective for understanding the etiology of pancreatic cancer,caution is still needed in data synthesis,selection of instrumental variables,and pleiotropy assessment.The use of emerging analytical models such as BWMR,CAUSE,and MVMR offers new possibilities for the comprehensive evaluation of multiple risk factors and their interaction.In the future,with the combination of these methods and the ever-increasing genetic epidemiological data,MR analysis is expected to provide more solid evidence for identifying potential therapeutic targets for pancreatic cancer and formulating prevention strategies.

    Key words:Pancreatic Neoplasms;Mendelian Randomization Analysis;Causality

    Research funding:National Natural Science Foundation of China(82160131);Qinghai Provincial Department of Science and Technology Project(2021-ZJ-963Q)

    1背景介紹

    1.1胰腺癌胰腺癌是消化系統(tǒng)中常見的惡性腫瘤之一,全球致死率極高,在常見癌癥發(fā)病率和死亡率中分別排名第12位和第7位。2018年數(shù)據(jù)[1-2]顯示,胰腺癌死亡人數(shù)與新發(fā)病例數(shù)近乎相當(dāng),且男性發(fā)病率稍高。研究[3]預(yù)測(cè),到2025年,胰腺癌將超過乳腺癌,成為癌癥死亡的第三大原因。胰腺癌患者5年生存率不足5%。在我國(guó),胰腺癌發(fā)病率、死亡率均位居前列[4-6]。

    由于胰腺深藏于腹腔并接鄰多個(gè)器官,胰腺癌的早期癥狀往往不明顯,易與其他消化道疾病混淆,導(dǎo)致誤診或漏診。目前,胰腺癌的臨床診斷主要依賴影像學(xué)、實(shí)驗(yàn)室檢查和病理學(xué)檢查等手段。盡管醫(yī)療技術(shù)已取得長(zhǎng)足進(jìn)步,但對(duì)于胰腺癌的早期篩查能力依舊有限。手術(shù)切除仍為最有效的治療方法[7-8]。胰腺癌的發(fā)生風(fēng)險(xiǎn)因素包括吸煙、肥胖、遺傳和糖尿病等,但其成因機(jī)制尚未完全明確[9-12]。因此,研究其與其他因素間的因果關(guān)系,對(duì)揭示胰腺癌的深層病理機(jī)制、制訂臨床干預(yù)手段和公共衛(wèi)生政策具有重要意義。

    1.2孟德爾隨機(jī)化孟德爾隨機(jī)化(mendelian randomization,MR)是利用與風(fēng)險(xiǎn)因素相關(guān)的遺傳變異作為工具變量(instrumental variable,IV),以推斷風(fēng)險(xiǎn)因素與疾病之間因果關(guān)系的方法。該方法基于孟德爾遺傳定律,認(rèn)為影響特定性狀的等位基因隨減數(shù)分裂遺傳給子代時(shí)隨機(jī)分配,從而避免了反向因果和混雜因素的影響,增強(qiáng)了研究結(jié)果的可靠性[13-14]。此方法的建立基于3個(gè)核心假設(shè)(圖1),同時(shí)也是IV選擇的原則:(1)相關(guān)性假設(shè)(假設(shè)一)——IV[如基因單核苷酸多態(tài)性(single nucleotide polymorphism,SNP)]必須與暴露因素(因果變量)強(qiáng)相關(guān)(利用SNP與暴露因素相關(guān)性的P值及弱IV的F值來保證IV與暴露因素的強(qiáng)相關(guān)性,以Plt;5×10-8、Fgt;10為默認(rèn)值);(2)獨(dú)立性假設(shè)(假設(shè)二)——IV與混淆因素?zé)o關(guān),這意味著IV不應(yīng)該與可能影響結(jié)果的未觀察到的因素相關(guān)聯(lián);(3)排他性假設(shè)(假設(shè)三)——IV與結(jié)果(即結(jié)局)無關(guān),這確保了IV僅通過影響暴露因素來影響結(jié)果,沒有直接影響結(jié)局[15]。同時(shí),在進(jìn)行IV的篩選時(shí)也需要完成連鎖不平衡以確保IV的可靠性(連鎖不平衡是實(shí)際遺傳過程中的一個(gè)現(xiàn)象,其中2個(gè)基因位點(diǎn)不是獨(dú)立遺傳的)[16]。在IV的選擇過程中,應(yīng)盡量避免連鎖不平衡,以確保IV的有效性??梢允褂胷2值和窗口(kb)來評(píng)估連鎖不平衡。隨著r2的減小和kb的增大,可以去除更多存在連鎖不平衡的SNP。

    采用MR-Egger、MR多效性殘差和離群值(MR-PRESSO)檢驗(yàn)水平多效性與離群值[17],應(yīng)用MR-Egger鑒定水平多效性的存在,Pgt;0.05提示水平多效性不顯著。與MR-Egger相比,MR-PRESSO在識(shí)別水平多效性和離群值方面具有較高的準(zhǔn)確率。多效性的評(píng)估加強(qiáng)了IV的可靠性。

    目前多采用兩樣本MR,常用方法有逆方差加權(quán)法(inverse variance weighted,IVW)、加權(quán)中位數(shù)法(weighted median method,WM)和MR-Egger法等,通過IV與疾病結(jié)局之間的關(guān)系來推測(cè)暴露因素的因果效應(yīng)[18]。MR有效規(guī)避了觀察性研究和隨機(jī)對(duì)照試驗(yàn)中的限制,但仍存在一些不足,如IV的選擇和生物多樣性假設(shè)驗(yàn)證困難,以及結(jié)果轉(zhuǎn)化應(yīng)用上的難點(diǎn)。因此,盡管處理了一些偏倚,確定因果關(guān)系的完整性仍有待多方法驗(yàn)證和隨機(jī)對(duì)照試驗(yàn)結(jié)果支持[15,19-20]。

    IVW是MR匯總數(shù)據(jù)的標(biāo)準(zhǔn)方法[21],可將每一個(gè)IV的Wald比率估計(jì)合并到元分析中,相當(dāng)于實(shí)現(xiàn)了IV與結(jié)果之間關(guān)聯(lián)的加權(quán)線性回歸,其IV的截距被限制為0。IVW的優(yōu)勢(shì)在于可以無偏估計(jì)狀態(tài)而沒有水平多效性。然而,IVW方法也容易受到弱工具偏倚的影響,有模擬研究[22]結(jié)果表明,IVW的弱工具偏倚大小可通過F統(tǒng)計(jì)量進(jìn)行量化。此外,WM可以減少第一類錯(cuò)誤的發(fā)生,并允許一些遺傳變異無效。當(dāng)大多數(shù)具有相似因果估計(jì)的IV有效時(shí),即使有部分IV不滿足MR方法的因果推斷要求,WM仍然是可信的。WM方法對(duì)遺傳變異與暴露因素之間的相關(guān)性進(jìn)行加權(quán),使得與暴露因素相關(guān)性較強(qiáng)的遺傳變異在效應(yīng)估計(jì)中占更大比重,從而降低弱IV(即與暴露因素相關(guān)性較弱的遺傳變異)對(duì)結(jié)果的影響[23],這可能導(dǎo)致弱IV對(duì)總體估計(jì)的貢獻(xiàn)較小。雖然這有助于減少弱IV可能引入的偏差,但也可能導(dǎo)致總體估計(jì)的統(tǒng)計(jì)效能降低,尤其是在沒有足夠強(qiáng)IV的情況下。MR-Egger方法依據(jù)inSIDE的假設(shè),主要反映IV與結(jié)果之間的劑量關(guān)系,并考慮了一定程度的多效性[18]。MR-Egger放寬多效性假設(shè)的代價(jià)是其違反NOME假設(shè),帶來的偏倚比IVW更大,且特別容易受到弱工具偏倚的影響[24]。另外,只有在基因多效性具有方向性的時(shí)候(即多效性具有非零平均值),MR-Egger回歸才能檢測(cè)出多效性[18]。當(dāng)幾種常用MR方法結(jié)果不一致時(shí),仍以經(jīng)典的IVW結(jié)果為主。

    1.3工具變量的獲取與過濾IV的數(shù)據(jù)可以通過多個(gè)公共數(shù)據(jù)庫獲取,這些數(shù)據(jù)庫包括但不限于IEU Open GWAS(https://gwas.mrcieu.ac.uk/)、FINNGEN(https://www.finngen.fi/)、UK Biobank(https://www.ukbiobank.ac.uk/)及GWAS Catalog(https://www.ebi.ac.uk/gwas/)。這些資源提供了豐富的遺傳變異信息,涵蓋SNP的基因型頻率、物理位置和已知的生物學(xué)功能。選擇與研究暴露因素相關(guān)的生物路徑中的IV至關(guān)重要,因此,研究人員需要對(duì)相關(guān)文獻(xiàn)和先前的遺傳關(guān)聯(lián)研究有深入的理解。

    在選擇IV時(shí),研究人員可以遵循但不限于以下步驟[25-28](TwoSampleMR、MRPRESSO包):(1)應(yīng)優(yōu)先選擇那些暴露與結(jié)局人群一致、發(fā)表時(shí)間較新的研究數(shù)據(jù)集,因?yàn)檫@些數(shù)據(jù)集更可能反映當(dāng)前的科研進(jìn)展和最新的科學(xué)發(fā)現(xiàn)。(2)初始設(shè)定與暴露相關(guān)的SNP的全基因組顯著性閾值5×10?8(默認(rèn))。如果符合此條件的IV數(shù)量有限,可以考慮采用更寬松的P值閾值,但這可能會(huì)減弱IV與暴露之間的相關(guān)性,即降低第一假設(shè)的強(qiáng)度。(3)利用1 000 Genomes Project的歐洲樣本數(shù)據(jù)作為參考,計(jì)算SNP之間的連鎖不平衡(LD),默認(rèn)LD參數(shù)設(shè)為r2lt;0.001,遺傳間距設(shè)為10 000 kb。r2用于衡量2個(gè)SNP之間的獨(dú)立性,其值越接近1,表示LD越強(qiáng),即2個(gè)SNP之間的遺傳相關(guān)性越強(qiáng)。遺傳間距是指在計(jì)算LD時(shí)考慮的最大物理距離,在此范圍內(nèi)的所有SNP都會(huì)被納入考慮,以確定哪些SNP違反了獨(dú)立性標(biāo)準(zhǔn)。如果符合此條件的IV數(shù)量有限,也可以采用更寬松的閾值,但需謹(jǐn)慎行事,以免將遺傳上相關(guān)的SNP錯(cuò)誤地視為獨(dú)立IV,影響MR分析的有效性。(4)為確保所選IV與暴露顯著相關(guān),使用F統(tǒng)計(jì)量評(píng)估IV的強(qiáng)度。F值的計(jì)算公式為F統(tǒng)計(jì)量,F(xiàn)值計(jì)算公式為F=(R2/(1-R2))×((N-K-1)/K)[29](R2的計(jì)算公式為:R2=2×EAF×(1-EAF)×β2 2×β2×EAF(1-EAF)+2×SE2×N×EAF×(1-EAF),其中EAF為效應(yīng)等位基因頻率,β為等位基因效應(yīng)值,SE為標(biāo)準(zhǔn)誤差,N為樣本量,K為自由度等于1),如果計(jì)算得到Flt;10,表明IV與暴露之間的聯(lián)系較弱,應(yīng)從分析中排除這類IV。(5)為確保所選IV與結(jié)局的風(fēng)險(xiǎn)因素?zé)o關(guān),需要查閱相關(guān)文獻(xiàn),并使用PhenoScanner數(shù)據(jù)庫或LDlink數(shù)據(jù)庫對(duì)這些IV進(jìn)行詳細(xì)分析,排除可能與這些風(fēng)險(xiǎn)因素存在潛在關(guān)聯(lián)的SNP。(6)在完成上述篩選步驟后,進(jìn)入數(shù)據(jù)提取與合并階段。從結(jié)局?jǐn)?shù)據(jù)中提取與篩選出的IV相對(duì)應(yīng)的SNP信息,確保分析聚焦于那些與暴露有顯著關(guān)聯(lián)且獨(dú)立的IV。(7)在處理遺傳數(shù)據(jù)時(shí),需特別注意回文SNP的影響?;匚腟NP在DNA序列中正反兩個(gè)方向的讀取序列相同,可能導(dǎo)致分析混淆。為確保分析的準(zhǔn)確性,設(shè)置去除回文SNP的參數(shù)action=2,代表中度強(qiáng)度的處理。具體而言,action=1處理低強(qiáng)度回文SNP,action=3處理高強(qiáng)度回文SNP,這些SNP應(yīng)立即排除。通過這些參數(shù)設(shè)置和處理步驟,可以有效地識(shí)別和去除回文SNP,提高研究的準(zhǔn)確性和可靠性。(8)在初步數(shù)據(jù)清洗和篩選后,使用MR-PRESSO工具進(jìn)一步識(shí)別和排除潛在的異常SNP。MR-PRESSO能夠檢測(cè)并處理MR分析中的異常值或偏離預(yù)期的SNP,這些異常SNP可能由于基因型測(cè)量誤差、樣本分層、基因多效性等原因表現(xiàn)出異常。通過MR-PRESSO的細(xì)致分析,可以精確地識(shí)別并排除這些SNP,提高數(shù)據(jù)集的純凈度和統(tǒng)計(jì)分析的可靠性。

    通過以上步驟,研究人員可以確保所選IV的質(zhì)量和相關(guān)性,為后續(xù)的MR分析提供堅(jiān)實(shí)的數(shù)據(jù)基礎(chǔ)。

    2 MR在胰腺癌研究中的應(yīng)用現(xiàn)狀

    目前,MR已在一定范圍內(nèi)應(yīng)用于胰腺癌研究,研究范圍包括腸道菌群、飲食來源的循環(huán)抗氧化劑、生物鐘類型、2型糖尿病、BMI、C反應(yīng)蛋白(CRP)、血液代謝物、炎癥性腸病及IL-1β與胰腺癌的因果關(guān)系。

    2.1腸道菌群與胰腺癌人體腸道內(nèi)存在著豐富的微生物群,廣泛參與人體健康與疾病的多種生理與病理過程[30],如肥胖、2型糖尿病、肝脂肪變性等[31]。當(dāng)前研究[32-33]已揭示胰腺癌的發(fā)生和發(fā)展與遺傳突變、飲食習(xí)慣及環(huán)境因素有關(guān)。隨著宏基因組學(xué)和16S核糖體RNA基因測(cè)序技術(shù)的進(jìn)步,腸道微生物組(GM)與包括胰腺癌在內(nèi)的胰腺疾病之間的相關(guān)性已被發(fā)現(xiàn)[34]。GM失調(diào)不僅直接影響腸道疾?。?5-36],還影響腸外器官,如胰臟和肝臟[37]。研究[38]表明,與健康個(gè)體相比,胰腺癌患者的口腔、胃腸道和胰腺中的微生物群發(fā)生多重變化,凸顯了GM在胰腺癌中的作用。GM的這種改變通過多種機(jī)制對(duì)胰腺癌產(chǎn)生影響。具體而言,GM可以與胰腺癌的各種風(fēng)險(xiǎn)因素相互作用,如肥胖和糖尿病、炎癥反應(yīng)等。這種炎癥反應(yīng)的激活可以通過不同途徑促進(jìn)GM及其分泌物進(jìn)入胰臟,包括血液傳播、淋巴轉(zhuǎn)移和胰腺逆流等[39]。此外,GM通過調(diào)節(jié)炎癥反應(yīng)、免疫細(xì)胞浸潤(rùn)等機(jī)制,影響胰腺癌的發(fā)展和進(jìn)展[40-42]。但具體的因果效應(yīng)尚未被完全明了。有研究[25]利用GM和胰腺癌的GWAS匯總數(shù)據(jù)行雙樣本MR,以評(píng)估GM和胰腺癌之間的因果效應(yīng),共分析119個(gè)屬水平GM,包含1 198個(gè)SNP,其中IVW被選為MR分析的主要方法,鑒定出Senegalimassilia是胰腺癌的保護(hù)因素,Odoribacter、Ruminiclostridium 9、Ruminococcaceae(UCG011)、Streptococcus是胰腺癌的易感因素(表1)。另一項(xiàng)MR研究[26]的IVW分析表明,Lentisphaerae、Lentisphaeria、Victivallales和Lachnospiraceae UCG004與胰腺癌風(fēng)險(xiǎn)呈負(fù)相關(guān)。此外,Erysipelotrichia、Erysipelotrichales、Erysipelotrichaceae、Flavonifractor和Terrisporobacter與胰腺癌的風(fēng)險(xiǎn)呈正相關(guān)(表1)。上述2項(xiàng)研究均提到其研究局限在于IV的閾值非常低(Plt;1×10?5),其生物學(xué)作用尚不明確,這可能導(dǎo)致多效性現(xiàn)象,即單個(gè)IV影響多個(gè)途徑或?qū)傩裕瑥亩で搜芯拷Y(jié)果。這可能是造成兩項(xiàng)研究結(jié)果不同的原因。上述研究填補(bǔ)了GM與胰腺癌之間的因果關(guān)系研究領(lǐng)域的空白,并為未來胰腺癌的靶向治療相關(guān)研究提供了重要方向。

    2.2飲食、營(yíng)養(yǎng)因素與胰腺癌氧化應(yīng)激在許多腫瘤的發(fā)生和發(fā)展中起著關(guān)鍵作用,而膳食抗氧化劑可能在其中扮演重要角色[43-46]。既往諸多研究探討了維生素(如維生素A、E、C)、類胡蘿卜素與消化系統(tǒng)腫瘤之間的關(guān)系,但由于研究方法限制和混雜因素,研究結(jié)果并不一致。MR研究[47]證實(shí),飲食來源的循環(huán)抗氧化劑不能降低胰腺癌的風(fēng)險(xiǎn)。因此,對(duì)于不缺乏相關(guān)抗氧化劑的人群,不建議額外服用維生素補(bǔ)充劑來預(yù)防胰腺癌的發(fā)生。另一項(xiàng)MR研究[27]證實(shí),飲食來源的抗氧化劑與胰腺癌無顯著相關(guān)性。研究還發(fā)現(xiàn),基因預(yù)測(cè)的飲食來源的更高抗壞血酸和視黃醇絕對(duì)濃度與胰腺癌的風(fēng)險(xiǎn)略有降低有關(guān)。上述研究的IV較少,結(jié)果可靠性不高,因此仍需其他研究方法加以驗(yàn)證。

    2.3生物節(jié)律與胰腺癌生物鐘類型(或稱晨昏節(jié)律類型),反映的是生物體在一天中自然傾向于活動(dòng)與休息的時(shí)間,其2種極端表現(xiàn)為偏好早晨(早起型)與偏好晚上(晚起型)[48-49]。研究[50]表明,晚起型個(gè)體可能因與不良生活方式相關(guān)聯(lián),從而增加癌癥發(fā)生風(fēng)險(xiǎn)。一項(xiàng)使用宏基因組測(cè)序分析的研究[51]發(fā)現(xiàn),早晨型生物鐘與艾氏菌屬的減少和棘桿菌屬的增加有關(guān)。這兩個(gè)腸道微生物屬可能在胃腸道癌癥的發(fā)展中發(fā)揮作用[52-53]。盡管有少數(shù)研究將生物鐘類型與消化道癌癥聯(lián)系起來,但尚無研究進(jìn)一步探索這一關(guān)聯(lián)。MR研究[54]結(jié)果證實(shí)了早起型與胰腺癌無關(guān)。

    2.4代謝與胰腺癌隨著全球糖尿病患者數(shù)量的飆升,2型糖尿病與多種惡性腫瘤之間的相關(guān)性逐漸受到關(guān)注[55]。流行病學(xué)研究[56]揭示了2型糖尿病與胰腺癌風(fēng)險(xiǎn)增加之間的相關(guān)性,考慮到胰腺癌不容樂觀的5年生存率,深入探究這一關(guān)聯(lián)至關(guān)重要。相關(guān)MR結(jié)果顯示,2型糖尿病與胰腺癌風(fēng)險(xiǎn)增加之間存在密切關(guān)聯(lián)[57]。在東亞人群中,2型糖尿病與胰腺癌風(fēng)險(xiǎn)增加24%有關(guān)[58](表2)。因此,針對(duì)2型糖尿病患者,特別是年輕發(fā)病人群,實(shí)施早期和密集的隨訪,及時(shí)控制血糖,可能是減少胰腺癌發(fā)生的有效預(yù)防策略。

    當(dāng)前研究[59]強(qiáng)調(diào)了肥胖與胰腺癌風(fēng)險(xiǎn)增加之間的聯(lián)系,肥胖通常與慢性炎癥狀態(tài)有關(guān),可能通過促炎因子的分泌促進(jìn)腫瘤生長(zhǎng)。研究[60-64]表明,全身性的輕度炎癥,特別是CRP水平提高,是肥胖與胰腺癌之間關(guān)系的潛在中介。MR研究[65]顯示,BMI、CRP與胰腺癌之間存在正向的因果關(guān)系(β3、β2),BMI和CRP之間也存在正向的因果關(guān)系(β1);控制CRP之后,BMI和胰腺癌之間的因果關(guān)系不再顯著(β5);控制BMI之后,CRP和胰腺癌之間仍然存在正向的因果關(guān)系(β4)(表2),利用因果效應(yīng)β1和β4計(jì)算CRP的中介效應(yīng)β5=30%(圖2),提示BMI對(duì)胰腺癌的總效應(yīng)中有30%是通過CRP實(shí)現(xiàn)的。這項(xiàng)研究利用中介MR的方法,在一定程度上證實(shí)了BMI、CRP和胰腺癌之間的因果關(guān)系,以及CRP介導(dǎo)的BMI對(duì)胰腺癌的重要性,對(duì)于指導(dǎo)臨床工作以及人群的飲食和生活方式的決策具有重要意義。

    2.5血液代謝物與胰腺癌代謝重編程被認(rèn)為是癌癥發(fā)展中不可或缺的一環(huán),特定代謝物甚至可能成為新型疾病標(biāo)志物,有助于癌癥風(fēng)險(xiǎn)評(píng)估和早期診斷[66-67]。既往研究[68-69]揭示,多種血清代謝物的變化與胰腺病變緊密相關(guān),且這些變化同時(shí)受腸道菌群的影響。在控制混雜因素的情況下,一項(xiàng)雙樣本MR研究[70]分析血漿中代謝物對(duì)胰腺導(dǎo)管腺癌(pancreatic ductaladenocarcinoma,PDAC)風(fēng)險(xiǎn)的潛在因果效應(yīng),通過表型全基因組MR(Phe-MR)分析進(jìn)一步評(píng)估以確定的代謝物為靶向的PDAC干預(yù)的潛在副作用,識(shí)別出44種獨(dú)特的代謝物與PDAC風(fēng)險(xiǎn)顯著相關(guān),這些代謝物可能有助于解釋PDAC的機(jī)制。

    2.6炎癥、免疫因素與胰腺癌炎癥性腸?。╥nflammatory bowel disease,IBD)作為一種慢性胃腸道炎癥性疾病,在患者長(zhǎng)期隨訪研究[71-72]中顯示可能增加多種癌癥的發(fā)生風(fēng)險(xiǎn)。值得關(guān)注的是,一些研究[73-75]評(píng)估了IBD與胰腺癌之間的潛在聯(lián)系,但結(jié)果不盡一致。一項(xiàng)MR研究[28]利用IWM、WM、MR-Egger等方法分析IBD與胰腺癌之間的因果關(guān)系,結(jié)果顯示,潰瘍性結(jié)腸炎與胰腺癌不存在因果關(guān)系,克羅恩病與胰腺癌的發(fā)生存在正相關(guān)(表3)。

    實(shí)驗(yàn)研究[76]發(fā)現(xiàn),調(diào)節(jié)IL-1β活性對(duì)緩解急性胰腺炎的嚴(yán)重程度及降低病死率發(fā)揮著關(guān)鍵作用。IL-1β過量表達(dá)可能導(dǎo)致慢性胰腺炎的發(fā)展[77],并在胰腺炎向胰腺癌進(jìn)展過程中發(fā)揮促進(jìn)作用,其中一個(gè)途徑是通過增加免疫抑制性B淋巴細(xì)胞數(shù)量,促進(jìn)腫瘤的生成[78]。IL-1β在肥胖與胰腺癌的關(guān)聯(lián)中也發(fā)揮作用[59],并與腫瘤增長(zhǎng)[79]和癌細(xì)胞遷移[80]有關(guān)。盡管有研究[81]觀察到IL-1β與某些胰腺腫瘤類型存在相關(guān)性,但是人群中IL-1與胰腺疾病風(fēng)險(xiǎn)之間的因果關(guān)系尚有待進(jìn)一步明確。這些關(guān)系的準(zhǔn)確評(píng)估不僅有助于加強(qiáng)對(duì)胰腺疾病的病理基礎(chǔ)的理解,還能夠?yàn)殚_發(fā)治療藥物(如抑制核糖體IL-1藥物阿那曲拉等)提供科學(xué)支撐。MR研究[82]證實(shí),IL-1RA與急、慢性胰腺炎及胰腺癌風(fēng)險(xiǎn)較低有關(guān),IL-1RA水平越高,疾病風(fēng)險(xiǎn)越低(表3)。上述結(jié)果揭示了IL-1RA在胰腺疾病發(fā)展中的作用,為IL-1抑制劑(如阿那白滯素)用于治療急、慢性胰腺炎及胰腺癌提供了一定依據(jù),未來還需更多臨床數(shù)據(jù)進(jìn)一步證實(shí)。

    3 MR在胰腺癌研究中的應(yīng)用潛力

    在胰腺癌的研究中,MR分析作為探討潛在因果關(guān)系的強(qiáng)有力工具,已經(jīng)取得了一定的前沿成果。目前,主流的策略依賴于經(jīng)典的兩樣本MR分析,該分析框架通過多種方法(例如IVW、WM、MR-Egger、Simple mode、Weighted mode)探索暴露因素與臨床結(jié)局之間的相互作用。特別是IVW,通過將每個(gè)IV的Wald比率整合為一次加權(quán)線性回歸,提供了一種無偏的狀態(tài)估計(jì),且在沒有水平多效性的情況下尤為有效[18]。WM則即使在高達(dá)50%的IV失效的情況下,依舊能夠提供可靠的因果關(guān)系估算。MR-Egger進(jìn)一步考慮了IV可能存在的多效性。與簡(jiǎn)單模式方法不同,Weighted mode方法在減少偏差及無效遺傳變異對(duì)結(jié)果的影響方面,提供了每個(gè)IV不同的權(quán)重。MR分析的敏感性檢驗(yàn)主要包括:MR-Egger方法評(píng)估方向性、多效性,留一法分析來評(píng)價(jià)單個(gè)IV的影響強(qiáng)度,MR-PRESSO評(píng)估水平多效性[17]。同時(shí),Cochran’sQ統(tǒng)計(jì)量被用于檢驗(yàn)SNP間的異質(zhì)性[83]。

    然而,盡管經(jīng)過細(xì)致的方法分析和敏感性檢驗(yàn),MR結(jié)果的可靠性仍受制于一些局限性因素,包括:多基因性帶來的多個(gè)弱SNP-暴露效應(yīng)的不確定性,觀察到的SNP同時(shí)影響多個(gè)特征(即多效性或水平多效性),以及在運(yùn)用數(shù)據(jù)進(jìn)行MR分析時(shí)可能遇到的連鎖不平衡等風(fēng)險(xiǎn)[19]。為了應(yīng)對(duì)這些挑戰(zhàn),Zhao等[19]提出了一種貝葉斯加權(quán)孟德爾隨機(jī)化方法(BWMR)。這種方法明確考慮了由于多基因性而產(chǎn)生的弱效應(yīng)相關(guān)的不確定性,并通過使用貝葉斯加權(quán)的異常值檢測(cè)來解決由于多效性導(dǎo)致的IV假設(shè)違反問題[19]。相較于傳統(tǒng)MR方法,BWMR通過考慮多基因性和多效性的影響,能夠提供更準(zhǔn)確和穩(wěn)健的因果推斷,尤其是在處理復(fù)雜疾病遺傳學(xué)研究中常見的弱效應(yīng)和異常值時(shí)。有研究[84]在分析免疫細(xì)胞與乳腺癌關(guān)系時(shí)使用BWMR檢驗(yàn)MR研究的結(jié)果,同樣的,在胰腺癌研究利用上述手段亦可以加強(qiáng)結(jié)果的可靠性。

    隨著全基因組關(guān)聯(lián)研究數(shù)量的激增,MR分析的機(jī)遇與挑戰(zhàn)并存?,F(xiàn)有的MR方法并非沒有弱點(diǎn),其中嚴(yán)格的假設(shè)常被現(xiàn)實(shí)打破,致使假陽性頻發(fā)。特別是水平多效性問題,當(dāng)一個(gè)遺傳變異同時(shí)影響多個(gè)特征時(shí),該問題即尤為復(fù)雜。為了解決這一問題,Morrison等[85]學(xué)者提出了一種新的MR策略——CAUSE分析,其能夠有效避免由水平多效性帶來的假陽性結(jié)論。Baranova等[86]在研究阿爾茨海默病與新型冠狀病毒感染因果影響時(shí),利用CAUSE評(píng)估MR估計(jì)中潛在的相關(guān)水平多效性,結(jié)果顯示,重癥新型冠狀病毒感染與阿爾茨海默病之間的因果模型優(yōu)于非因果模型。Micha?lsson等[87]同樣在研究中加入了CAUSE-MR以加強(qiáng)結(jié)果的可靠性。相較于傳統(tǒng)MR方法,CAUSE-MR允許在評(píng)估因果關(guān)系時(shí)存在一定程度的多效性,即允許遺傳IV同時(shí)影響多個(gè)暴露因素。這與傳統(tǒng)的MR方法(如IVW)不同,后者假設(shè)IV僅通過影響單一暴露因素來影響結(jié)局,在存在多效性時(shí)可能會(huì)導(dǎo)致偏差[85]。因此,CAUSE-MR如果應(yīng)用于胰腺癌研究,可輔助評(píng)估傳統(tǒng)MR結(jié)果,加強(qiáng)結(jié)果的可靠性。

    考慮到傳統(tǒng)MR方法的局限性,研究不能僅停留在單一遺傳變異與單一結(jié)果的關(guān)聯(lián)分析。多變量MR(MVMR)通過包含模型中的其他暴露條件分析,能夠同時(shí)估算多個(gè)相關(guān)暴露的獨(dú)立效應(yīng),更準(zhǔn)確地揭示復(fù)雜的因果關(guān)系網(wǎng)絡(luò)[88]。在某些情況下,遺傳變異會(huì)多效性地與多個(gè)相關(guān)表型相關(guān)聯(lián)。例如,與脂蛋白代謝相關(guān)的遺傳變異很少只與某一特定脂質(zhì)分?jǐn)?shù)相關(guān)聯(lián)[89-90]。傳統(tǒng)MR方法很可能會(huì)因?yàn)榇嬖谶@種水平多效性而導(dǎo)致關(guān)于因果關(guān)系的誤導(dǎo)性結(jié)論。多變量MR能夠通過使用與多個(gè)暴露因素相關(guān)的IV來共同估計(jì)每個(gè)風(fēng)險(xiǎn)因素對(duì)結(jié)果的獨(dú)立因果效應(yīng),從而克服這個(gè)問題[91-93]。在胰腺癌研究中,可能需要分析胰腺癌與血糖水平、胰島素抵抗、肥胖和脂質(zhì)代謝等多個(gè)因素的關(guān)系,而多變量MR可適用于同時(shí)評(píng)估這些代謝因素對(duì)胰腺癌風(fēng)險(xiǎn)的獨(dú)立影響。

    4 MR的解讀

    在進(jìn)行MR敏感性分析時(shí),多效性檢測(cè)和評(píng)估是至關(guān)重要的一環(huán)。多效性是指一個(gè)SNP同時(shí)影響多個(gè)表型(如暴露因素和疾病結(jié)果)。這可能導(dǎo)致對(duì)因果關(guān)系的誤判,因?yàn)槎嘈钥赡苁沟肐V與結(jié)果之間存在非因果的關(guān)聯(lián)。因此,在MR分析中,剔除具有多效性的離群值對(duì)于確保研究結(jié)果的準(zhǔn)確性和可靠性至關(guān)重要。MR-Egger、MR多效性殘差和離群值(MR-PRESSO)檢驗(yàn)是常用方法,用于檢測(cè)和評(píng)估多效性及其對(duì)結(jié)果的影響。MR-Egger回歸可通過引入一個(gè)截距項(xiàng)來檢測(cè)水平多效性的存在,如果截距顯著不為0,則表明存在多效性[18]。MR-PRESSO方法則可以通過迭代地識(shí)別并剔除具有較大效應(yīng)估計(jì)的離群值,評(píng)估多效性對(duì)結(jié)果的影響[17]。剔除離群值的目的是減少多效性對(duì)因果估計(jì)的干擾,從而提高研究結(jié)果的準(zhǔn)確性。在剔除離群值時(shí),需要權(quán)衡剔除帶來的潛在偏差和提高結(jié)果準(zhǔn)確性之間的利弊。過度剔除可能導(dǎo)致結(jié)果的過度校正,從而忽視一些真實(shí)的因果關(guān)系。然而,不剔除離群值可能會(huì)使研究結(jié)果受到多效性的影響,從而導(dǎo)致錯(cuò)誤的因果推斷。因此,在進(jìn)行MR分析時(shí),應(yīng)綜合考慮多效性檢測(cè)結(jié)果、生物學(xué)機(jī)制及現(xiàn)有研究證據(jù),以決定是否剔除離群值以及如何解釋剔除后的結(jié)果。

    在MR的應(yīng)用中也會(huì)出現(xiàn)結(jié)果不一致的情況,可能是由于多種因素造成的。首先,IV的選擇和強(qiáng)度可能影響結(jié)果的可靠性,弱IV或與多個(gè)途徑相關(guān)的多效性IV可能導(dǎo)致結(jié)果的偏倚。其次,遺傳數(shù)據(jù)的來源、樣本量大小及潛在的混雜因素也可能影響MR分析的準(zhǔn)確性。此外,不同的分析方法(如IVW、WM)可能對(duì)多效性的敏感性不同,從而產(chǎn)生不同的結(jié)論。作為臨床工作者,在解讀MR結(jié)果時(shí),應(yīng)綜合考慮這些潛在的偏差來源。首先,評(píng)估所選IV的相關(guān)性和有效性,即是否滿足獨(dú)立性、排他性和相關(guān)性假設(shè)。其次,考慮使用多種MR方法并進(jìn)行敏感性分析,以檢驗(yàn)結(jié)果的穩(wěn)健性。最后,結(jié)合生物學(xué)機(jī)制和現(xiàn)有的臨床知識(shí)來解釋MR結(jié)果,不應(yīng)僅僅依賴單一的統(tǒng)計(jì)分析。在臨床實(shí)踐中,MR分析提供的證據(jù)應(yīng)被視為補(bǔ)充傳統(tǒng)臨床研究的額外信息,而非替代性證據(jù)。通過這種方式,臨床工作者可以更全面地理解疾病機(jī)制,并為患者提供基于證據(jù)的醫(yī)療服務(wù)。

    5小結(jié)與展望

    總而言之,MR作為一種解析潛在病因關(guān)系的工具,在胰腺癌研究中的應(yīng)用展現(xiàn)了其獨(dú)特的價(jià)值,特別是在解決傳統(tǒng)觀察性研究中的因果推斷難題方面,取得了令人矚目的成果。MR利用遺傳變異作為IV,這些變異在生命早期隨機(jī)分配,因此可以作為一種自然實(shí)驗(yàn),減少傳統(tǒng)研究中常見的混雜因素和反向因果關(guān)系的影響。在胰腺癌的研究中,MR方法能夠幫助研究者探索多種潛在的風(fēng)險(xiǎn)因素,如遺傳變異、生活方式因素(如飲食習(xí)慣)、代謝性疾?。ㄈ?型糖尿?。┘把装Y標(biāo)志物等與胰腺癌之間的因果關(guān)系。MR分析同時(shí)揭示了腸道菌群失調(diào)、循環(huán)抗氧化劑水平、炎癥性腸病及IL-1β等與胰腺癌風(fēng)險(xiǎn)之間的關(guān)系,這些發(fā)現(xiàn)對(duì)于理解胰腺癌的復(fù)雜病理機(jī)制至關(guān)重要。此外,MR分析還提供了一種評(píng)估這些風(fēng)險(xiǎn)因素在預(yù)防和治療胰腺癌中潛在作用的方法,有助于指導(dǎo)臨床決策和公共衛(wèi)生政策的制訂。面對(duì)越來越多的GWAS數(shù)據(jù),MR分析策略的優(yōu)化和創(chuàng)新正在不斷發(fā)展。新的算法和策略,如BWMR、CAUSE和MVMR,為克服傳統(tǒng)MR方法的限制開辟了道路,并進(jìn)一步挖掘出潛在風(fēng)險(xiǎn)因素之間的復(fù)雜因果鏈條。目前,BWMR、CAUSE、MVMR等方法尚未有運(yùn)用到胰腺癌研究的報(bào)道出現(xiàn),上述方法的合理運(yùn)用,將在傳統(tǒng)MR分析的基礎(chǔ)上,為胰腺癌與風(fēng)險(xiǎn)因子之間的因果關(guān)系提供更強(qiáng)有力的證據(jù),從而促進(jìn)潛在治療靶點(diǎn)的研究,豐富對(duì)腫瘤生物學(xué)的認(rèn)識(shí)。此外,需要認(rèn)識(shí)到的是,盡管先進(jìn)的統(tǒng)計(jì)模型為研究提供了強(qiáng)有力的工具,真實(shí)的生物學(xué)機(jī)制和多樣化的病因構(gòu)成了更復(fù)雜的挑戰(zhàn)。期待未來通過多學(xué)科協(xié)作研究,整合臨床、基因組學(xué)和流行病學(xué)數(shù)據(jù),為胰腺癌的早期預(yù)防、精準(zhǔn)診療和個(gè)性化醫(yī)療提供新的策略。

    利益沖突聲明:本文不存在任何利益沖突。

    作者貢獻(xiàn)聲明:杜凱豪負(fù)責(zé)課題設(shè)計(jì),論文的撰寫與修改;羅蘭明慧、東小鴿、蔣威參與文獻(xiàn)搜集,修改論文;王展、侯立朝負(fù)責(zé)擬定寫作思路,指導(dǎo)撰寫文章,并最后定稿。

    參考文獻(xiàn):

    [1]ARNOLD M,ABNET CC,NEALE RE,et al.Global burden of 5 majortypes of gastrointestinal cancer[J].Gastroenterology,2020,159(1):335-349.e15.DOI:10.1053/j.gastro.2020.02.068.

    [2]ILIC M,ILIC I.Epidemiology of pancreatic cancer[J].World J Gas?troenterol,2016,22(44):9694-9705.DOI:10.3748/wjg.v22.i44.9694.

    [3]SUNG H,F(xiàn)ERLAY J,SIEGEL RL,et al.Global cancer statistics 2020:GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J].CA Cancer J Clin,2021,71(3):209-249.DOI:10.3322/caac.21660.

    [4]FERLAY J,SOERJOMATARAM I,DIKSHIT R,et al.Cancer inci?dence and mortality worldwide:Sources,methods and major patterns in GLOBOCAN 2012[J].Int J Cancer,2015,136(5):E359-E386.DOI:10.1002/ijc.29210.

    [5]HIDALGO M,CASCINU S,KLEEFF J,et al.Addressing the chal?lenges of pancreatic cancer:Future directions for improving out?comes[J].Pancreatology,2015,15(1):8-18.DOI:10.1016/j.pan.2014.10.001.

    [6]General Office of National Health Commission.Standard for diagno?sis and treatment of pancreatic cancer(2022 edition)[J].J Clin Hepatol,2022,38(5):1006-1015.DOI:10.3969/j.issn.1001-5256.2022.05.007.

    國(guó)家衛(wèi)生健康委辦公廳.胰腺癌診療指南(2022年版)[J].臨床肝膽病雜志,2022,38(5):1006-1015.DOI:10.3969/j.issn.1001-5256.2022.05.007.

    [7]MILLIKAN KW,DEZIEL DJ,SILVERSTEIN JC,et al.Prognostic fac?tors associated with resectable adenocarcinoma of the head of the pancreas[J].Am Surg,1999,65(7):618-623;discussion 623-624.

    [8]KOLBEINSSON HM,CHANDANA S,WRIGHT GP,et al.Pancre?atic cancer:A review of current treatment and novel therapies[J].J Invest Surg,2023,36(1):2129884.DOI:10.1080/08941939.2022.2129884.

    [9]PARKIN DM,BOYD L,WALKER LC.16.The fraction of cancer attrib?utable to lifestyle and environmental factors in the UK in 2010[J].Br J Cancer,2011,105(Suppl 2):S77-S81.DOI:10.1038/bjc.2011.489.

    [10]WILLETT WC.Diet and cancer[J].Oncol,2000,5(5):393-404.DOI:10.1634/theoncologist.5-5-393.

    [11]BOSETTI C,BERTUCCIO P,NEGRI E,et al.Pancreatic cancer:Overview of descriptive epidemiology[J].Mol Carcinog,2012,51(1):3-13.DOI:10.1002/mc.20785.

    [12]ANAND P,KUNNUMAKKARA AB,SUNDARAM C,et al.Cancer is a preventable disease that requires major lifestyle changes[J].Pharm Res,2008,25(9):2097-2116.DOI:10.1007/s11095-008-9661-9.

    [13]WANG LN,ZHANG ZF.Mendelian randomization approach,used for causal inferences[J].Chin J Epidemiol,2017,38(4):547-552.DOI:10.3760/cma.j.issn.0254-6450.2017.04.027.

    王莉娜,ZHANG ZuoFeng.孟德爾隨機(jī)化法在因果推斷中的應(yīng)用[J].中華流行病學(xué)雜志,2017,38(4):547-552.DOI:10.3760/cma.j.issn.0254-6450.2017.04.027.

    [14]DAVEY SMITH G,HEMANI G.Mendelian randomization:Genetic an?chors for causal inference in epidemiological studies[J].Hum Mol Genet,2014,23(R1):R89-R98.DOI:10.1093/hmg/ddu328.

    [15]ZHENG J,BAIRD D,BORGES MC,et al.Recent developments in Mendelian randomization studies[J].Curr Epidemiol Rep,2017,4(4):330-345.DOI:10.1007/s40471-017-0128-6.

    [16]LAWLOR DA,HARBORD RM,STERNE JA,et al.Mendelian random?ization:Using genes as instruments for making causal inferences in epidemiology[J].Stat Med,2008,27(8):1133-1163.DOI:10.1002/sim.3034.

    [17]VERBANCK M,CHEN CY,NEALE B,et al.Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases[J].Nat Genet,2018,50(5):693-698.DOI:10.1038/s41588-018-0099-7.

    [18]BOWDEN J,DAVEY SMITH G,BURGESS S.Mendelian randomiza?tion with invalid instruments:Effect estimation and bias detectionthrough Egger regression[J].Int J Epidemiol,2015,44(2):512-525.DOI:10.1093/ije/dyv080.

    [19]ZHAO J,MING JS,HU XH,et al.Bayesian weighted Mendelian ran?domization for causal inference based on summary statistics[J].Bioinformatics,2020,36(5):1501-1508.DOI:10.1093/bioinformatics/btz749.

    [20]LAWLOR DA.Commentary:Two-sample Mendelian randomization:Opportunities and challenges[J].Int J Epidemiol,2016,45(3):908-915.DOI:10.1093/ije/dyw127.

    [21]BURGESS S,BUTTERWORTH A,THOMPSON SG.Mendelian ran?domization analysis with multiple genetic variants using summarized data[J].Genet Epidemiol,2013,37(7):658-665.DOI:10.1002/gepi.21758.

    [22]YU TQ,XU WT,SU YN,et al.Basic principles,methods and limita?tions of Mendelian randomization studies[J].Chin J Evid-Based Med,2021,21(10):1227-1234.DOI:10.7507/1672-2531.202107008.

    于天琦,徐文濤,蘇雅娜,等.孟德爾隨機(jī)化研究基本原理、方法和局限性[J].中國(guó)循證醫(yī)學(xué)雜志,2021,21(10):1227-1234.DOI:10.7507/1672-2531.202107008.

    [23]BOWDEN J,DAVEY SMITH G,HAYCOCK PC,et al.Consistent esti?mation in Mendelian randomization with some invalid instruments us?ing a weighted Median estimator[J].Genet Epidemiol,2016,40(4):304-314.DOI:10.1002/gepi.21965.

    [24]BOWDEN J,DEL GRECO MF,MINELLI C,et al.Assessing the suit?ability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression:The role of the I2 statistic[J].Int J Epidemiol,2016,45(6):1961-1974.DOI:10.1093/ije/dyw220.

    [25]JIANG ZC,MOU YP,WANG HJ,et al.Causal effect between gut mi?crobiota and pancreatic cancer:A two-sample Mendelian random?ization study[J].BMC Cancer,2023,23(1):1091.DOI:10.1186/s12885-023-11493-y.

    [26]ZHU SJ,DING Z.Association between gut microbiota and sevengastrointestinal diseases:A Mendelian randomized study[J].J Gene Med,2024,26(1):e3623.DOI:10.1002/jgm.3623.

    [27]ZHANG XN,ZHAO H,MAN JY,et al.Investigating causal associa?tions of diet-derived circulating antioxidants with the risk of digestive system cancers:A Mendelian randomization study[J].Nutrients,2022,14(15):3237.DOI:10.3390/nu14153237.

    [28]MIN Y,LIU ZR,LI RD,et al.Association between inflammatory bowel disease and pancreatic cancer:Results from the two-sample Mendelian randomization study[J].Front Oncol,2023,13:1155123.DOI:10.3389/fonc.2023.1155123.

    [29]PIERCE BL,AHSAN H,VANDERWEELE TJ.Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants[J].Int J Epidemiol,2011,40(3):740-752.DOI:10.1093/ije/dyq151.

    [30]ADAK A,KHAN MR.An insight into gut microbiota and its function?alities[J].Cell Mol Life Sci,2019,76(3):473-493.DOI:10.1007/s00018-018-2943-4.

    [31]de VOS WM,TILG H,van HUL M,et al.Gut microbiome and health:Mechanistic insights[J].Gut,2022,71(5):1020-1032.DOI:10.1136/gutjnl-2021-326789.

    [32]LIU MY.Arid1a:A gatekeeper in the development of pancreatic can?cer from a rare precursor lesion[J].Gastroenterology,2022,163(2):371-373.DOI:10.1053/j.gastro.2022.05.046.

    [33]CAI J,CHEN HD,LU M,et al.Advances in the epidemiology of pan?creatic cancer:Trends,risk factors,screening,and prognosis[J].Cancer Lett,2021,520:1-11.DOI:10.1016/j.canlet.2021.06.027.

    [34]TONG Y,GAO HR,QI QC,et al.High fat diet,gut microbiome and gastrointestinal cancer[J].Theranostics,2021,11(12):5889-5910.DOI:10.7150/thno.56157.

    [35]COHEN LJ,CHO JH,GEVERS D,et al.Genetic factors and the in?testinal microbiome guide development of microbe-based therapiesfor inflammatory bowel diseases[J].Gastroenterology,2019,156(8):2174-2189.DOI:10.1053/j.gastro.2019.03.017.

    [36]WONG SH,YU J.Gut microbiota in colorectal cancer:Mechanisms of action and clinical applications[J].Nat Rev Gastroenterol Hepa?tol,2019,16(11):690-704.DOI:10.1038/s41575-019-0209-8.

    [37]MA C,HAN MJ,HEINRICH B,et al.Gut microbiome-mediated bile acid metabolism regulates liver cancer via NKT cells[J].Science,2018,360(6391):eaan5931.DOI:10.1126/science.aan5931.

    [38]ERTZ-ARCHAMBAULT N,KEIM P,von HOFF D.Microbiome and pancreatic cancer:A comprehensive topic review of literature[J].World J Gastroenterol,2017,23(10):1899-1908.DOI:10.3748/wjg.v23.i10.1899.

    [39]YANG QY,ZHANG JH,ZHU Y.Potential roles of the gut microbiota in pancreatic carcinogenesis and therapeutics[J].Front Cell Infect Microbiol,2022,12:872019.DOI:10.3389/fcimb.2022.872019.

    [40]PUSHALKAR S,HUNDEYIN M,DALEY D,et al.The pancreatic can?cer microbiome promotes oncogenesis by induction of innate and adaptive immune suppression[J].Cancer Discov,2018,8(4):403-416.DOI:10.1158/2159-8290.CD-17-1134.

    [41]SETHI V,KURTOM S,TARIQUE M,et al.Gut microbiota promotes tumor growth in mice by modulating immune response[J].Gastro?enterology,2018,155(1):33-37.e6.DOI:10.1053/j.gastro.2018.04.001.

    [42]YU Q,NEWSOME RC,BEVERIDGE M,et al.Intestinal microbiota modulates pancreatic carcinogenesis through intratumoral natural killer cells[J].Gut Microbes,2022,14(1):2112881.DOI:10.1080/19490976.2022.2112881.

    [43]CARINI F,MAZZOLA M,RAPPA F,et al.Colorectal carcinogenesis:Role of oxidative stress and antioxidants[J].Anticancer Res,2017,37(9):4759-4766.DOI:10.21873/anticanres.11882.

    [44]SHIMIZU Y,TAMURA T,KEMMOCHI A,et al.Oxidative stress andLiver X Receptor agonist induce hepatocellular carcinoma in Non-al?coholic steatohepatitis model[J].J Gastroenterol Hepatol,2021,36(3):800-810.DOI:10.1111/jgh.15239.

    [45]AKBARI A,MAJD HM,RAHNAMA R,et al.Cross-talk between oxi?dative stress signaling and microRNA regulatory systems in carci?nogenesis:Focused on gastrointestinal cancers[J].Biomedecine Pharmacother,2020,131:110729.DOI:10.1016/j.biopha.2020.110729.

    [46]KLAUNIG JE.Oxidative stress and cancer[J].Curr Pharm Des,2018,24(40):4771-4778.DOI:10.2174/1381612825666190215121712.

    [47]YIN LL,YAN HH,CHEN KD,et al.Diet-derived circulating antioxi?dants and risk of digestive system tumors:A Mendelian randomizationstudy[J].Nutrients,2022,14(16):3274.DOI:10.3390/nu14163274.

    [48]SEMPERE-RUBIO N,AGUAS M,F(xiàn)AUBEL R.Association betweenchronotype,physical activity and sedentary behaviour:A systematic review[J].Int J Environ Res Public Health,2022,19(15):9646.DOI:10.3390/ijerph19159646.

    [49]MAZRI FH,MANAF ZA,SHAHAR S,et al.The association between chronotype and dietary pattern among adults:A scoping review[J].Int J Environ Res Public Health,2019,17(1):68.DOI:10.3390/ijerph17010068.

    [50]SROUR B,PLANCOULAINE S,ANDREEVA VA,et al.Circadian nutri?tional behaviours and cancer risk:New insights from the NutriNet-santéprospective cohort study:Disclaimers[J].Int J Cancer,2018,143(10):2369-2379.DOI:10.1002/ijc.31584.

    [51]CARASSO S,F(xiàn)ISHMAN B,LASK LS,et al.Metagenomic analysis re?veals the signature of gut microbiota associated with human chronotypes[J].FASEB J,2021,35(11):e22011.DOI:10.1096/fj.202100857RR.

    [52]MENG CT,BAI CM,BROWN TD,et al.Human gut microbiota and gastrointestinal cancer[J].Genomics Proteomics Bioinformatics,2018,16(1):33-49.DOI:10.1016/j.gpb.2017.06.002.

    [53]FENG Q,LIANG SS,JIA HJ,et al.Gut microbiome development along the colorectal adenoma-carcinoma sequence[J].Nat Com?mun,2015,6:6528.DOI:10.1038/ncomms7528.

    [54]YUAN S,MASON AM,TITOVA OE,et al.Morning chronotype and di?gestive tract cancers:Mendelian randomization study[J].Int J Can?cer,2023,152(4):697-704.DOI:10.1002/ijc.34284.

    [55]RISK FACTOR COLLABORATION(NCD-RISC)NCD.Worldwide trends in diabetes since 1980:A pooled analysis of 751 population-based studies with 4.4 million participants[J].Lancet,2016,387(10027):1513-1530.DOI:10.1016/S0140-6736(16)00618-8.

    [56]TSILIDIS KK,KASIMIS JC,LOPEZ DS,et al.Type 2 diabetes and cancer:Umbrella review of meta-analyses of observational studies[J].BMJ,2015,350:g7607.DOI:10.1136/bmj.g7607.

    [57]PEARSON-STUTTARD J,PAPADIMITRIOU N,MARKOZANNES G,et al.Type 2 diabetes and cancer:An umbrella review of observa?tional and Mendelian randomization studies[J].Cancer Epidemiol Biomarkers Prev,2021,30(6):1218-1228.DOI:10.1158/1055-9965.EPI-20-1245.

    [58]SHEN BY,LI YY,SHENG CS,et al.Association between age at dia?betes onset or diabetes duration and subsequent risk of pancreatic cancer:Results from a longitudinal cohort and Mendelian random?ization study[J].Lancet Reg Health West Pac,2023,30:100596.DOI:10.1016/j.lanwpc.2022.100596.

    [59]INCIO J,LIU H,SUBOJ P,et al.Obesity-induced inflammation and desmoplasia promote pancreatic cancer progression and resistance to chemotherapy[J].Cancer Discov,2016,6(8):852-869.DOI:10.1158/2159-8290.CD-15-1177.

    [60]PARK JH,AMERI AH,DEMPSEY KE,et al.Nuclear IL-33/SMAD sig?naling axis promotes cancer development in chronic inflammation[J].EMBO J,2021,40(7):e106151.DOI:10.15252/embj.2020106151.

    [61]YUDKIN JS,STEHOUWER CD,EMEIS JJ,et al.C-reactive protein in healthy subjects:Associations with obesity,insulin resistance,and endothelial dysfunction:A potential role for cytokines originating from adipose tissue?[J].Arterioscler Thromb Vasc Biol,1999,19(4):972-978.DOI:10.1161/01.atv.19.4.972.

    [62]XU HY,BARNES GT,YANG Q,et al.Chronic inflammation in fat plays a crucial role in the development of obesity-related insulin re?sistance[J].J Clin Invest,2003,112(12):1821-1830.DOI:10.1172/JCI19451.

    [63]SINGH A,MAYENGBAM SS,YADUVANSHI H,et al.Obesity pro?grams macrophages to support cancer progression[J].Cancer Res,2022,82(23):4303-4312.DOI:10.1158/0008-5472.CAN-22-1257.

    [64]XIA B,HE QS,PAN YH,et al.Metabolic syndrome and risk of pan?creatic cancer:A population-based prospective cohort study[J].Int J Cancer,2020,147(12):3384-3393.DOI:10.1002/ijc.33172.

    [65]LI ZQ,JIN LQ,XIA L,et al.Body mass index,C-reactive protein,and pancreatic cancer:A Mendelian randomization analysis to in?vestigate causal pathways[J].Front Oncol,2023,13:1042567.DOI:10.3389/fonc.2023.1042567.

    [66]SON J,LYSSIOTIS CA,YING HQ,et al.Glutamine supports pancre?atic cancer growth through a KRAS-regulated metabolic pathway[J].Nature,2013,496(7443):101-105.DOI:10.1038/nature12040.

    [67]CHEUNG PK,MA MH,TSE HF,et al.The applications of metabolo?mics in the molecular diagnostics of cancer[J].Expert Rev Mol Di?agn,2019,19(9):785-793.DOI:10.1080/14737159.2019.1656530.

    [68]XIE GX,LU LG,QIU YP,et al.Plasma metabolite biomarkers for the detection of pancreatic cancer[J].J Proteome Res,2015,14(2):1195-1202.DOI:10.1021/pr501135f.

    [69]LI QX,JIN M,LIU YH,et al.Gut microbiota:Its potential roles in pan?creatic cancer[J].Front Cell Infect Microbiol,2020,10:572492.DOI:10.3389/fcimb.2020.572492.

    [70]ZHONG H,LIU S,ZHU JJ,et al.Elucidating the role of blood me?tabolites on pancreatic cancer risk using two-sample Mendelian ran?domization analysis[J].Int J Cancer,2024,154(5):852-862.DOI:10.1002/ijc.34771.

    [71]NADEEM MS,KUMAR V,AL-ABBASI FA,et al.Risk of colorectal cancer in inflammatory bowel diseases[J].Semin Cancer Biol,2020,64:51-60.DOI:10.1016/j.semcancer.2019.05.001.

    [72]ZHANG HM,ZHANG MM,CHEN XF,et al.Risk of malignancy in pa?tients with inflammatory bowel disease:A population-based cohort study from China[J].Int J Cancer,2022,150(11):1770-1778.DOI:10.1002/ijc.33932.

    [73]JUNG YS,HAN M,PARK S,et al.Cancer risk in the early stages of inflammatory bowel disease in Korean patients:A nationwide population-based study[J].J Crohns Colitis,2017,11(8):954-962.DOI:10.1093/ecco-jcc/jjx040.

    [74]EVERHOV?H,ERICHSEN R,SACHS MC,et al.Inflammatory bowel disease and pancreatic cancer:A Scandinavian register-based co?hort study 1969-2017[J].Aliment Pharmacol Ther,2020,52(1):143-154.DOI:10.1111/apt.15785.

    [75]LO B,ZHAO M,VIND I,et al.The risk of extraintestinal cancer in in?flammatory bowel disease:A systematic review and meta-analysis of population-based cohort studies[J].Clin Gastroenterol Hepatol,2021,19(6):1117-1138.e19.DOI:10.1016/j.cgh.2020.08.015.

    [76]DENHAM W,YANG J,F(xiàn)INK G,et al.Gene targeting demonstrates additive detrimental effects of interleukin 1 and tumor necrosis factorduring pancreatitis[J].Gastroenterology,1997,113(5):1741-1746.DOI:10.1053/gast.1997.v113.pm9352880.

    [77]MARRACHE F,TU SP,BHAGAT G,et al.Overexpression of interleu?kin-1beta in the murine pancreas results in chronic pancreatitis[J].Gastroenterology,2008,135(4):1277-1287.DOI:10.1053/j.gastro.2008.06.078.

    [78]TAKAHASHI R,MACCHINI M,SUNAGAWA M,et al.Interleukin-1β-induced pancreatitis promotes pancreatic ductal adenocarcinoma via B lymphocyte-mediated immune suppression[J].Gut,2021,70(2):330-341.DOI:10.1136/gutjnl-2019-319912.

    [79]ZHANG Y,CHEN XY,WANG HM,et al.Innate immune mediator,In?terleukin-1 receptor accessory protein(IL1RAP),is expressed and pro-tumorigenic in pancreatic cancer[J].J Hematol Oncol,2022,15(1):70.DOI:10.1186/s13045-022-01286-4.

    [80]TJOMSLAND V,SANDNES D,POMIANOWSKA E,et al.The TGFβ-SMAD3 pathway inhibits IL-1“induced interactions between human pancreatic stellate cells and pancreatic carcinoma cells and re?stricts cancer cell migration[J].J Exp Clin Cancer Res,2016,35(1):122.DOI:10.1186/s13046-016-0400-5.

    [81]MAKER AV,KATABI N,QIN LX,et al.Cyst fluid interleukin-1beta(IL1beta)levels predict the risk of carcinoma in intraductal papillary mucinous neoplasms of the pancreas[J].Clin Cancer Res,2011,17(6):1502-1508.DOI:10.1158/1078-0432.CCR-10-1561.

    [82]YUAN S,MIAO YY,RUAN XX,et al.Therapeutic role of interleukin-1 receptor antagonist in pancreatic diseases:Mendelian randomization study[J].Front Immunol,2023,14:1240754.DOI:10.3389/fimmu.2023.1240754.

    [83]GRECO M FD,MINELLI C,SHEEHAN NA,et al.Detecting pleiot?ropy in Mendelian randomisation studies with summary data and a continuous outcome[J].Stat Med,2015,34(21):2926-2940.DOI:10.1002/sim.6522.

    [84]WANG X,GAO HY,ZENG YY,et al.A Mendelian analysis of the rela?tionships between immune cells and breast cancer[J].Front Oncol,2024,14:1341292.DOI:10.3389/fonc.2024.1341292.

    [85]MORRISON J,KNOBLAUCH N,MARCUS JH,et al.Mendelian ran?domization accounting for correlated and uncorrelated pleiotropic effects using genome-wide summary statistics[J].Nat Genet,2020,52(7):740-747.DOI:10.1038/s41588-020-0631-4.

    [86]BARANOVA A,CAO HB,ZHANG FQ.Causal effect of COVID-19 on Alzheimer’s disease:A Mendelian randomization study[J].J Med Virol,2023,95(1):e28107.DOI:10.1002/jmv.28107.

    [87]MICHA?LSSON M,YUAN S,MELHUS H,et al.The impact and causal directions for the associations between diagnosis of ADHD,socioeconomic status,and intelligence by use of a bi-directional two-sample Mendelian randomization design[J].BMC Med,2022,20(1):106.DOI:10.1186/s12916-022-02314-3.

    [88]SANDERSON E.Multivariable Mendelian randomization and media?tion[J].Cold Spring Harb Perspect Med,2021,11(2):a038984.DOI:10.1101/cshperspect.a038984.

    [89]WILLER CJ,SCHMIDT EM,SENGUPTA S,et al.Discovery and re?finement of loci associated with lipid levels[J].Nat Genet,2013,45(11):1274-1283.DOI:10.1038/ng.2797.

    [90]WüRTZ P,KANGAS AJ,SOININEN P,et al.Lipoprotein subclass profiling reveals pleiotropy in the genetic variants of lipid risk factors for coronary heart disease:A note on Mendelian randomization stud?ies[J].J Am Coll Cardiol,2013,62(20):1906-1908.DOI:10.1016/j.jacc.2013.07.085.

    [91]BURGESS S,THOMPSON SG.Multivariable Mendelian randomiza?tion:The use of pleiotropic genetic variants to estimate causal ef?fects[J].Am J Epidemiol,2015,181(4):251-260.DOI:10.1093/aje/kwu283.

    [92]BURGESS S,DUDBRIDGE F,THOMPSON SG.Re:“Multivariable Mendelian randomization:The use of pleiotropic genetic variants to estimate causal effects”[J].Am J Epidemiol,2015,181(4):290-291.DOI:10.1093/aje/kwv017.

    [93]KEMP JP,SAYERS A,SMITH GD,et al.Using Mendelian randomiza?tion to investigate a possible causal relationship between adiposity and increased bone mineral density at different skeletal sites in chil?dren[J].Int J Epidemiol,2016,45(5):1560-1572.DOI:10.1093/ije/dyw079.

    收稿日期:2024-01-21;錄用日期:2024-04-19

    本文編輯:邢翔宇

    引證本文:DU KH, HOU LZ, LUO LMH, et al. Current status and prospects of the application of Mendelian randomization in pancreatic cancer research[J]. J Clin Hepatol, 2024, 40 (10): 2127-2136.

    杜凱豪, 侯立朝, 羅蘭明慧, 等 . 孟德爾隨機(jī)化在胰腺癌研究中的 應(yīng)用現(xiàn)狀與展望[J]. 臨床肝膽病雜志, 2024, 40(10): 2127- 2136

    午夜a级毛片| 久久久精品欧美日韩精品| 国产精品爽爽va在线观看网站 | 亚洲精品国产一区二区精华液| 可以免费在线观看a视频的电影网站| 免费久久久久久久精品成人欧美视频| 黄色毛片三级朝国网站| 亚洲自偷自拍图片 自拍| av免费在线观看网站| 999久久久精品免费观看国产| 日本免费一区二区三区高清不卡 | 一本久久中文字幕| 很黄的视频免费| 久久亚洲真实| 大陆偷拍与自拍| 国产成人精品在线电影| 男人的好看免费观看在线视频 | 国产亚洲精品第一综合不卡| 1024香蕉在线观看| av在线天堂中文字幕| 狠狠狠狠99中文字幕| 一级片免费观看大全| 真人一进一出gif抽搐免费| 美女 人体艺术 gogo| 大码成人一级视频| 窝窝影院91人妻| 欧美性长视频在线观看| 伊人久久大香线蕉亚洲五| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品,欧美在线| 在线天堂中文资源库| 亚洲七黄色美女视频| 欧美成狂野欧美在线观看| 99久久精品国产亚洲精品| 欧美另类亚洲清纯唯美| 国产精品一区二区免费欧美| 欧美 亚洲 国产 日韩一| 99国产精品一区二区三区| 精品一区二区三区视频在线观看免费| 侵犯人妻中文字幕一二三四区| 中文字幕高清在线视频| 日韩中文字幕欧美一区二区| 国产在线观看jvid| 深夜精品福利| 香蕉国产在线看| 久久久精品国产亚洲av高清涩受| 看免费av毛片| 少妇 在线观看| 亚洲全国av大片| 精品一品国产午夜福利视频| 一个人免费在线观看的高清视频| 欧美成人一区二区免费高清观看 | 性欧美人与动物交配| 国产片内射在线| 久久国产亚洲av麻豆专区| 亚洲国产精品合色在线| 日韩免费av在线播放| 美女免费视频网站| 天堂动漫精品| 男女下面插进去视频免费观看| 后天国语完整版免费观看| 香蕉国产在线看| 国产日韩一区二区三区精品不卡| 欧美成狂野欧美在线观看| 黄色成人免费大全| 国产激情久久老熟女| 在线av久久热| 国产欧美日韩综合在线一区二区| 国产精品亚洲av一区麻豆| 热99re8久久精品国产| 桃红色精品国产亚洲av| 变态另类丝袜制服| 国产乱人伦免费视频| 咕卡用的链子| 黄色a级毛片大全视频| 国产国语露脸激情在线看| 成人手机av| 狠狠狠狠99中文字幕| 亚洲av电影在线进入| 日韩欧美一区二区三区在线观看| 老汉色∧v一级毛片| 在线观看66精品国产| 国产免费男女视频| 免费观看精品视频网站| 人妻久久中文字幕网| 亚洲人成伊人成综合网2020| tocl精华| 中文字幕最新亚洲高清| www日本在线高清视频| 神马国产精品三级电影在线观看 | 啦啦啦观看免费观看视频高清 | 18美女黄网站色大片免费观看| 成人手机av| 精品日产1卡2卡| 丝袜美足系列| 色哟哟哟哟哟哟| 国产主播在线观看一区二区| 亚洲人成77777在线视频| 国产熟女午夜一区二区三区| 亚洲人成电影免费在线| 成人18禁在线播放| 精品福利观看| 中文字幕人成人乱码亚洲影| 欧美日韩亚洲综合一区二区三区_| 69av精品久久久久久| 欧美成人一区二区免费高清观看 | 国产成人精品无人区| 久久九九热精品免费| 女人被狂操c到高潮| 精品第一国产精品| 757午夜福利合集在线观看| 18禁观看日本| 视频区欧美日本亚洲| 人人妻人人澡人人看| 国产精品一区二区免费欧美| av免费在线观看网站| 在线视频色国产色| 欧美激情极品国产一区二区三区| 久久精品成人免费网站| 成熟少妇高潮喷水视频| 91九色精品人成在线观看| 不卡一级毛片| 长腿黑丝高跟| 欧美中文日本在线观看视频| 制服人妻中文乱码| 激情在线观看视频在线高清| 村上凉子中文字幕在线| 99国产综合亚洲精品| 操出白浆在线播放| 国产区一区二久久| 久久久久久久久中文| 精品福利观看| 午夜福利18| 久久精品影院6| 黄片大片在线免费观看| 高潮久久久久久久久久久不卡| 激情在线观看视频在线高清| 久久精品国产亚洲av香蕉五月| 人成视频在线观看免费观看| а√天堂www在线а√下载| 欧美日韩福利视频一区二区| 免费一级毛片在线播放高清视频 | 午夜免费观看网址| 一边摸一边抽搐一进一小说| 热99re8久久精品国产| 超碰成人久久| 女人高潮潮喷娇喘18禁视频| 在线观看午夜福利视频| 久久国产精品人妻蜜桃| 国产野战对白在线观看| 久久久国产精品麻豆| 一区在线观看完整版| 亚洲欧洲精品一区二区精品久久久| 窝窝影院91人妻| 精品不卡国产一区二区三区| 国产欧美日韩一区二区三| 国产av精品麻豆| 国产欧美日韩综合在线一区二区| √禁漫天堂资源中文www| 国产精品乱码一区二三区的特点 | 久久久久久久精品吃奶| 国产亚洲欧美98| 久久 成人 亚洲| 成人特级黄色片久久久久久久| 国产又色又爽无遮挡免费看| 性色av乱码一区二区三区2| 成年人黄色毛片网站| 色av中文字幕| 黄网站色视频无遮挡免费观看| 熟妇人妻久久中文字幕3abv| 大型黄色视频在线免费观看| 国产精品爽爽va在线观看网站 | 18禁美女被吸乳视频| 国产亚洲av嫩草精品影院| 色av中文字幕| 免费在线观看视频国产中文字幕亚洲| 大型黄色视频在线免费观看| 亚洲成人国产一区在线观看| 久久久久久久精品吃奶| 一级,二级,三级黄色视频| www.999成人在线观看| 色播亚洲综合网| 一进一出抽搐动态| 国产成人av激情在线播放| 国产亚洲av高清不卡| 在线永久观看黄色视频| 久久久国产精品麻豆| 黄色片一级片一级黄色片| 啦啦啦观看免费观看视频高清 | 国产伦一二天堂av在线观看| 精品久久久久久,| 日本免费a在线| 久久这里只有精品19| 91麻豆av在线| 国产精品 欧美亚洲| 日韩精品中文字幕看吧| 99re在线观看精品视频| 啪啪无遮挡十八禁网站| 精品日产1卡2卡| 免费在线观看视频国产中文字幕亚洲| 97人妻精品一区二区三区麻豆 | 琪琪午夜伦伦电影理论片6080| 亚洲男人天堂网一区| 午夜免费激情av| 级片在线观看| 国产精品98久久久久久宅男小说| 一区福利在线观看| 久久精品91无色码中文字幕| 两性夫妻黄色片| 91大片在线观看| 岛国在线观看网站| 丁香欧美五月| 电影成人av| www国产在线视频色| 在线观看舔阴道视频| 日韩欧美一区二区三区在线观看| 制服诱惑二区| 99久久99久久久精品蜜桃| 亚洲av成人不卡在线观看播放网| 黄色视频,在线免费观看| 最新在线观看一区二区三区| 国产乱人伦免费视频| 美女大奶头视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美绝顶高潮抽搐喷水| 亚洲国产看品久久| 亚洲成av人片免费观看| 国产精品久久电影中文字幕| 一本大道久久a久久精品| 久久人人精品亚洲av| 亚洲人成77777在线视频| 精品熟女少妇八av免费久了| 国产区一区二久久| 中文亚洲av片在线观看爽| 国产精品,欧美在线| 日韩精品免费视频一区二区三区| 黑丝袜美女国产一区| 97超级碰碰碰精品色视频在线观看| 久久精品亚洲精品国产色婷小说| 男女下面进入的视频免费午夜 | 亚洲专区国产一区二区| 国产精品久久久久久精品电影 | 亚洲av美国av| tocl精华| 波多野结衣巨乳人妻| 999久久久精品免费观看国产| 啦啦啦 在线观看视频| 精品熟女少妇八av免费久了| 国产亚洲精品av在线| 两个人看的免费小视频| 精品人妻1区二区| 禁无遮挡网站| 亚洲av片天天在线观看| 神马国产精品三级电影在线观看 | 成人三级做爰电影| 精品欧美一区二区三区在线| 国产成人欧美在线观看| netflix在线观看网站| 午夜免费成人在线视频| 亚洲欧美一区二区三区黑人| 午夜久久久在线观看| 手机成人av网站| 天天躁狠狠躁夜夜躁狠狠躁| 88av欧美| av福利片在线| 好看av亚洲va欧美ⅴa在| 亚洲成人精品中文字幕电影| 日韩 欧美 亚洲 中文字幕| 99热只有精品国产| 午夜老司机福利片| 一边摸一边抽搐一进一出视频| 国产成+人综合+亚洲专区| 国产成人啪精品午夜网站| 成人精品一区二区免费| 久久人人精品亚洲av| 亚洲av熟女| 在线播放国产精品三级| 欧美另类亚洲清纯唯美| 亚洲精品中文字幕在线视频| 丝袜美腿诱惑在线| 欧洲精品卡2卡3卡4卡5卡区| 国产精品电影一区二区三区| 亚洲全国av大片| 国产精华一区二区三区| 99国产精品99久久久久| 国产精品影院久久| 91成年电影在线观看| e午夜精品久久久久久久| 欧美 亚洲 国产 日韩一| 国产主播在线观看一区二区| 香蕉国产在线看| 亚洲色图综合在线观看| 男女午夜视频在线观看| 天天添夜夜摸| 两人在一起打扑克的视频| 老司机深夜福利视频在线观看| 在线观看一区二区三区| 一区在线观看完整版| 日日摸夜夜添夜夜添小说| av天堂在线播放| 51午夜福利影视在线观看| 日韩欧美在线二视频| 看片在线看免费视频| 无遮挡黄片免费观看| 黑人巨大精品欧美一区二区蜜桃| 欧美午夜高清在线| 非洲黑人性xxxx精品又粗又长| 精品国产国语对白av| 91在线观看av| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成狂野欧美在线观看| 淫妇啪啪啪对白视频| 12—13女人毛片做爰片一| 大陆偷拍与自拍| 精品午夜福利视频在线观看一区| 精品不卡国产一区二区三区| 深夜精品福利| 老司机福利观看| 亚洲欧美日韩另类电影网站| 国产精品二区激情视频| 91国产中文字幕| 亚洲中文日韩欧美视频| 久久婷婷成人综合色麻豆| 国产精品一区二区三区四区久久 | 久久久久国产一级毛片高清牌| 一本大道久久a久久精品| 午夜福利一区二区在线看| 男人的好看免费观看在线视频 | 纯流量卡能插随身wifi吗| 视频区欧美日本亚洲| 最新在线观看一区二区三区| 亚洲电影在线观看av| 一区二区日韩欧美中文字幕| 国产熟女午夜一区二区三区| 9热在线视频观看99| 他把我摸到了高潮在线观看| 欧美色欧美亚洲另类二区 | 不卡av一区二区三区| av欧美777| 亚洲精品中文字幕在线视频| 国内毛片毛片毛片毛片毛片| 精品日产1卡2卡| 老汉色∧v一级毛片| 国产伦人伦偷精品视频| 国产精华一区二区三区| 国产免费男女视频| 在线国产一区二区在线| av欧美777| 色播亚洲综合网| 日日爽夜夜爽网站| 亚洲自拍偷在线| 日韩av在线大香蕉| 精品国内亚洲2022精品成人| 亚洲中文字幕日韩| 国产一区二区三区综合在线观看| 久久国产乱子伦精品免费另类| 午夜免费激情av| av电影中文网址| 国产91精品成人一区二区三区| 波多野结衣一区麻豆| 久久久久精品国产欧美久久久| 麻豆久久精品国产亚洲av| 国产一区二区在线av高清观看| 在线av久久热| 久久伊人香网站| 欧美一区二区精品小视频在线| 国产精品爽爽va在线观看网站 | videosex国产| 91麻豆av在线| av免费在线观看网站| 99在线人妻在线中文字幕| 国产成人欧美在线观看| 不卡一级毛片| 69av精品久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 国产高清videossex| 欧美绝顶高潮抽搐喷水| 国产精品九九99| www国产在线视频色| 久久精品91蜜桃| 亚洲狠狠婷婷综合久久图片| 老司机福利观看| 久久人人97超碰香蕉20202| 国产亚洲精品久久久久5区| 成人三级做爰电影| 在线观看66精品国产| 国产精品自产拍在线观看55亚洲| 99国产综合亚洲精品| 久久国产精品人妻蜜桃| 男女床上黄色一级片免费看| 日本撒尿小便嘘嘘汇集6| 夜夜躁狠狠躁天天躁| 亚洲美女黄片视频| 香蕉丝袜av| 91大片在线观看| www.自偷自拍.com| 亚洲一卡2卡3卡4卡5卡精品中文| 久久人人爽av亚洲精品天堂| 亚洲精品久久成人aⅴ小说| 制服丝袜大香蕉在线| 少妇的丰满在线观看| 精品一区二区三区视频在线观看免费| 国产成人精品久久二区二区91| 99精品在免费线老司机午夜| 亚洲成人免费电影在线观看| 色播亚洲综合网| 91精品国产国语对白视频| 男女做爰动态图高潮gif福利片 | 999久久久国产精品视频| 国产精品爽爽va在线观看网站 | 午夜影院日韩av| 日韩成人在线观看一区二区三区| 久久久久久久午夜电影| 久久午夜亚洲精品久久| 国产国语露脸激情在线看| 亚洲中文日韩欧美视频| 又紧又爽又黄一区二区| 97人妻精品一区二区三区麻豆 | 黑人操中国人逼视频| www.自偷自拍.com| 美女高潮到喷水免费观看| 久久亚洲真实| 国产激情欧美一区二区| 久久欧美精品欧美久久欧美| 1024视频免费在线观看| 欧美国产精品va在线观看不卡| 欧美黑人欧美精品刺激| 精品高清国产在线一区| 亚洲男人天堂网一区| 国产精品亚洲av一区麻豆| 好看av亚洲va欧美ⅴa在| 久久香蕉国产精品| 9191精品国产免费久久| 午夜福利免费观看在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产精品999在线| 人人妻人人澡人人看| 97人妻天天添夜夜摸| 成人免费观看视频高清| 欧美一级a爱片免费观看看 | 女警被强在线播放| 亚洲午夜理论影院| 亚洲熟妇熟女久久| 国产aⅴ精品一区二区三区波| 97人妻精品一区二区三区麻豆 | 人人妻人人澡欧美一区二区 | 欧美日韩亚洲国产一区二区在线观看| 国产午夜福利久久久久久| 十八禁网站免费在线| 成人欧美大片| 欧美在线黄色| 99精品久久久久人妻精品| 国产精品国产高清国产av| 亚洲中文av在线| 久久热在线av| 黄色女人牲交| 中出人妻视频一区二区| 亚洲 欧美 日韩 在线 免费| 国产精品 国内视频| 757午夜福利合集在线观看| 色在线成人网| 麻豆国产av国片精品| 制服诱惑二区| 在线观看免费视频日本深夜| 91精品三级在线观看| 男人的好看免费观看在线视频 | 女性被躁到高潮视频| 丝袜美腿诱惑在线| 成人永久免费在线观看视频| 亚洲男人天堂网一区| 成熟少妇高潮喷水视频| 亚洲精品一区av在线观看| 如日韩欧美国产精品一区二区三区| 精品不卡国产一区二区三区| 亚洲七黄色美女视频| 国产99久久九九免费精品| 18禁国产床啪视频网站| 亚洲精品国产一区二区精华液| 国产乱人伦免费视频| 黄频高清免费视频| 欧美日韩福利视频一区二区| 欧美乱码精品一区二区三区| 露出奶头的视频| 欧美另类亚洲清纯唯美| 999久久久精品免费观看国产| 一级,二级,三级黄色视频| 免费在线观看视频国产中文字幕亚洲| 90打野战视频偷拍视频| 丝袜人妻中文字幕| 免费搜索国产男女视频| 欧美日本视频| 日韩视频一区二区在线观看| 夜夜夜夜夜久久久久| 精品高清国产在线一区| 国产一区二区在线av高清观看| 久久人妻av系列| 日本免费一区二区三区高清不卡 | 国产真人三级小视频在线观看| 少妇的丰满在线观看| 久久久久国产精品人妻aⅴ院| 亚洲性夜色夜夜综合| 麻豆一二三区av精品| АⅤ资源中文在线天堂| 性欧美人与动物交配| 亚洲va日本ⅴa欧美va伊人久久| 国产精品精品国产色婷婷| 日本三级黄在线观看| 超碰成人久久| av网站免费在线观看视频| 9热在线视频观看99| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 国产又爽黄色视频| 免费人成视频x8x8入口观看| 老司机在亚洲福利影院| 一二三四社区在线视频社区8| 高潮久久久久久久久久久不卡| 在线免费观看的www视频| 国产极品粉嫩免费观看在线| 别揉我奶头~嗯~啊~动态视频| 麻豆国产av国片精品| 麻豆国产av国片精品| 成年版毛片免费区| 在线天堂中文资源库| 亚洲免费av在线视频| 亚洲av五月六月丁香网| 久久精品国产亚洲av香蕉五月| 午夜精品在线福利| 一本大道久久a久久精品| 黄色视频,在线免费观看| 国产亚洲精品综合一区在线观看 | 成在线人永久免费视频| 成人精品一区二区免费| 丰满的人妻完整版| 欧美中文综合在线视频| 精品卡一卡二卡四卡免费| 波多野结衣高清无吗| 操出白浆在线播放| 一区二区三区高清视频在线| 男人操女人黄网站| 动漫黄色视频在线观看| 国产熟女午夜一区二区三区| 国产精品综合久久久久久久免费 | 国产区一区二久久| 久久人人精品亚洲av| 色综合站精品国产| 欧美国产日韩亚洲一区| 欧美黑人精品巨大| 亚洲欧洲精品一区二区精品久久久| 在线观看一区二区三区| 日韩欧美在线二视频| 一区二区日韩欧美中文字幕| 亚洲av电影不卡..在线观看| avwww免费| 极品人妻少妇av视频| 女人高潮潮喷娇喘18禁视频| 亚洲成a人片在线一区二区| 欧美av亚洲av综合av国产av| 午夜福利在线观看吧| 精品不卡国产一区二区三区| 亚洲精华国产精华精| 午夜激情av网站| 亚洲精华国产精华精| 少妇 在线观看| 免费人成视频x8x8入口观看| 亚洲电影在线观看av| 欧美人与性动交α欧美精品济南到| 女人被躁到高潮嗷嗷叫费观| 日本一区二区免费在线视频| 精品人妻1区二区| x7x7x7水蜜桃| 丰满的人妻完整版| 淫妇啪啪啪对白视频| 91成年电影在线观看| 亚洲专区字幕在线| 亚洲国产精品成人综合色| 久久久久国产精品人妻aⅴ院| 亚洲第一av免费看| www国产在线视频色| 国产免费av片在线观看野外av| www.www免费av| 麻豆成人av在线观看| 久久香蕉激情| 久久久久国产精品人妻aⅴ院| 麻豆av在线久日| 手机成人av网站| 高清毛片免费观看视频网站| 午夜两性在线视频| 在线观看日韩欧美| 午夜激情av网站| 国产亚洲av嫩草精品影院| 欧美精品亚洲一区二区| 亚洲成av人片免费观看| 无限看片的www在线观看| 最近最新中文字幕大全免费视频| 侵犯人妻中文字幕一二三四区| 成年人黄色毛片网站| 成人国产一区最新在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 一二三四社区在线视频社区8| 精品欧美国产一区二区三| 一进一出抽搐gif免费好疼| 久久这里只有精品19| 97人妻精品一区二区三区麻豆 | 黑人操中国人逼视频| 国产蜜桃级精品一区二区三区| 亚洲国产日韩欧美精品在线观看 | 欧美最黄视频在线播放免费| 视频区欧美日本亚洲| 高潮久久久久久久久久久不卡| 久久久久久久久久久久大奶| 欧美精品亚洲一区二区| 在线av久久热| 亚洲成人精品中文字幕电影| 亚洲第一电影网av|