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

    兩樣本孟德爾隨機(jī)化分析4種常量營養(yǎng)素相對攝入量與腎結(jié)石的因果關(guān)系

    2024-12-13 00:00:00王聰席羽佳王振興賈宜君郝川
    現(xiàn)代泌尿外科雜志 2024年12期

    Two sample Mendelian randomization analysis of the causal relationship between relative intake of four constant nutrients and kidney stones

    WANG Cong1,2,XI Yujia1,2,WANG Zhenxing2,JIA Yijun1,2,HAO Chuan2

    (1.The Second Clinical Medical College,2.Department of Urology,The Second Hospital,Shanxi Medical University,Taiyuan 030001,China)

    ABSTRACT:Objective To investigate the causal relationship between the relative intake of dietary components (including carbohydrate,added sugar,protein,and fat) and kidney stones using two sample Mendelian randomization (MR).Methods The genetic instrumental variables (IVs) from genome-wide association study (GWAS) statistics,exposure-related summary data from the SSGAC genome study,and kidney stone data from the Finn Gen Consortium were collected,which were then analyzed with inverse variance weighting (IVW).The MR results and sensitivity were assessed with MR-Egger,weighted median,MR-PRESSO test,Cochran’s Q test,leave-one analysis,and Steiger filter.Results IVW results showed that the relative intake of carbohydrate (OR=2.01,95%CI:1.02-3.96,P=0.042) and added sugar (OR=2.07,95%CI:1.00-4.29,P=0.049) had a causal relationship and were risk factors for the development of kidney stones.Relative intake of protein (OR=0.70,95%CI:0.38-1.29,P=0.249) and fat (OR=0.79,95%CI:0.27-2.35,P=0.671) were not associated with kidney stones.Sensitivity analysis showed no heterogeneity or pleiotropy (Pgt;0.05).Conclusion Relative intake of carbohydrate and added sugar are risk factors for kidney stones,suggesting that limiting carbohydrate and added sugar intake may prevent the development of kidney stones.

    KEY WORDS:Mendelian randomization; kidney stones; diet; carbohydrate; added sugar

    摘要:目的 通過兩樣本孟德爾隨機(jī)化(MR)探討飲食成分相對攝入量(即碳水化合物、添加糖、蛋白質(zhì)和脂肪)與腎結(jié)石(KSD)之間的因果關(guān)系,為KSD的預(yù)防提供參考。方法 檢索整理全基因組關(guān)聯(lián)研究(GWAS)統(tǒng)計數(shù)據(jù),暴露組數(shù)據(jù)來自社會科學(xué)遺傳協(xié)會聯(lián)合會的基因組研究,KSD數(shù)據(jù)選自Finn Gen生物樣本庫。使用逆方差加權(quán)法(IVW)作為兩樣本MR分析的主要方法,采用MR-Egger、加權(quán)中位數(shù)法、MR-多效性殘差和離群值(MR-PRESSO)檢驗(yàn)、Cochran’s Q檢驗(yàn)、“留一法”分析和Steiger濾波作為MR結(jié)果和敏感性分析的補(bǔ)充。結(jié)果 IVW結(jié)果顯示碳水化合物相對攝入量(OR=2.01,95%CI:1.02~3.96,P=0.042)、添加糖相對攝入量(OR=2.07,95%CI:1.00~4.29,P=0.049)與KSD之間存在因果關(guān)系,是其發(fā)生發(fā)展的危險因素。未發(fā)現(xiàn)蛋白質(zhì)相對攝入量(OR=0.70,95%CI:0.38~1.29,P=0.249)和脂肪相對攝入量(OR=0.79,95%CI:0.27~2.35,P=0.671)與KSD存在關(guān)聯(lián)。敏感性分析顯示上述結(jié)果無異質(zhì)性和多效性(P>0.05)。結(jié)論 碳水化合物相對攝入量與添加糖相對攝入量可能是KSD的危險因素,提示控制碳水化合物與添加糖的攝入可預(yù)防KSD的發(fā)生。

    關(guān)鍵詞:孟德爾隨機(jī)化;腎結(jié)石;飲食;碳水化合物;添加糖

    中圖分類號:R692.4 文獻(xiàn)標(biāo)志碼:A

    DOI:10.3969/j.issn.1009-8291.2024.12.011

    收稿日期:2024-03-07 修回日期:2024-06-02

    基金項目:山西省醫(yī)師協(xié)會 (No.YSXH-QL2023MN002)

    通信作者:郝川,碩士研究生導(dǎo)師,主任醫(yī)師。

    E-mail:haochuan2003@sina.com

    作者簡介:王聰,碩士研究生在讀。研究方向:泌尿系結(jié)石。

    E-mail:wc18834563080@163.com

    腎結(jié)石(kidney stone,KSD)是泌尿系統(tǒng)最常見的疾病之一,其臨床表現(xiàn)包括腰部或上腹部陣發(fā)性劇烈絞痛、血尿、腎積水等。這些臨床表現(xiàn)不但給患者帶來了巨大的身體和精神負(fù)擔(dān)[1],而且頗高的復(fù)發(fā)率也增加了治療的經(jīng)濟(jì)衛(wèi)生成本[2]。因此,針對KSD的誘發(fā)因素制定詳細(xì)的預(yù)防及治療策略是一個全球性的醫(yī)療保健問題。有研究表明通過調(diào)控飲食和生活方式,KSD發(fā)病率可降低50%以上[3]。先后有隨機(jī)對照研究及觀察性研究報道,4種常量營養(yǎng)素(即碳水化合物、蛋白質(zhì)、脂肪、添加糖)相對攝入量的改變會增加或降低KSD發(fā)生發(fā)展的風(fēng)險。然而上述研究采用傳統(tǒng)研究方法受混雜因素和反向因果關(guān)系的影響,結(jié)論在很大程度上存在偏倚性。

    為了得到更明確的結(jié)果,本研究利用兩樣本孟德爾隨機(jī)化(mendelian randomization,MR)方法來推斷上述4種飲食的攝入量與KSD發(fā)生的潛在因果關(guān)系。與傳統(tǒng)對照試驗(yàn)相比,此方法利用了樣本量龐大的數(shù)據(jù)庫,避免了長期隨訪這一耗時且高成本的環(huán)節(jié)。此外,MR選用了單核苷酸多態(tài)位點(diǎn)(single nucleotide polymorphism,SNP)作為工具變量(instrumental variables,IVs),而基因在受孕時具有隨機(jī)分配性,基因變異先于疾病的發(fā)展,且不受外界環(huán)境、社會行為等因素的影響,因此在推斷因果關(guān)聯(lián)時具有合理的時序性,很大程度上減少了傳統(tǒng)觀察性研究中存在的混雜因素和反向因果關(guān)系影響[4],并且本研究使用的暴露與結(jié)局樣本無重疊,保證了研究結(jié)果的準(zhǔn)確性,為KSD的飲食預(yù)防提供更精確的證據(jù)。

    1 資料與方法

    1.1 研究設(shè)計 采用兩樣本MR模型評價飲食成分相對攝入量(即碳水化合物、添加糖、脂肪、蛋白質(zhì))與腎結(jié)石之間的潛在因果關(guān)系。IVs滿足以下3個假設(shè)[5]:①相關(guān)性假設(shè):IVs與暴露(即4種營養(yǎng)素相對攝入量)間密切相關(guān),P<5×10-8為強(qiáng)相關(guān);②排他性假設(shè):確保IVs與結(jié)局(腎結(jié)石)無直接聯(lián)系,在結(jié)局變量的全基因組關(guān)聯(lián)研究(genome-wide association study,GWAS)數(shù)據(jù)中查找并篩選出IVs,同時剔除與結(jié)局因素高度相關(guān)的IVs,以確保IVs不直接與結(jié)局產(chǎn)生關(guān)聯(lián);③獨(dú)立性假設(shè):IVs與暴露和結(jié)局之間潛在的混雜因素?zé)o關(guān)。本研究的流程圖見圖1。

    1.2 數(shù)據(jù)來源 關(guān)于KSD的IVs,本研究選擇了Finn Gen(https://www.finngen.fi/en)生物樣本庫中的GWAS數(shù)據(jù)。在調(diào)整年齡、性別、遺傳相關(guān)性、基因分型批次和前10個主要成分后,收集到關(guān)于KSD的相關(guān)數(shù)據(jù)包括376 406名芬蘭成人受試者,其中病例有9713例,對照組有366 693例[6]。

    暴露組的數(shù)據(jù)來源于社會科學(xué)遺傳協(xié)會聯(lián)合會(Social Science Genetic Associations Consortium,SSGAC)由MEDDENS等[7]科研團(tuán)隊整理的GWAS研究數(shù)據(jù),包括碳水化合物相對攝入量、蛋白質(zhì)相對攝入量、脂肪相對攝入量、添加糖相對攝入量,該研究對超過 235 000 名具有歐洲血統(tǒng)的人進(jìn)行問卷案調(diào)查,調(diào)查所得的24 h飲食回憶問卷的信息被用于英國生物庫(UK Biobank,UKB)隊列中進(jìn)行分析,食物頻率問卷的結(jié)果用于復(fù)制隊列?!俺A繝I養(yǎng)素密度”是通過將常量營養(yǎng)素攝入量除以總攝入量而獲得的。由于常量營養(yǎng)素的攝入量與總攝入量呈非線性關(guān)系,作者使用修正因子β進(jìn)行修正,即

    相對攝入量=常量營養(yǎng)素的攝入量/總攝入量×β

    最終采用上式來衡量脂肪、蛋白質(zhì)、碳水化合物和添加糖的相對攝入量。

    MR研究要求SNP暴露和SNP結(jié)局關(guān)聯(lián)的數(shù)據(jù)具有獨(dú)立性[8],本研究暴露與結(jié)局IVs來源于不同的生物樣本庫,樣本之間基本不存在重疊,并且受試者的倫理批準(zhǔn)和書面知情同意情況在原始出版物中均有說明。

    1.3 篩選IVs IVs篩選標(biāo)準(zhǔn):IVs與暴露間強(qiáng)相關(guān)P<5×10-8。為減少偏倚結(jié)果的發(fā)生,以r2<0.001、窗口大小為10 000 kb的等位基因?yàn)榧s束條件,消除具有強(qiáng)連鎖不平衡(linkage disequilibrium,LD)的SNPs[9]。使用F統(tǒng)計量來確定是否存在弱 IVs偏倚[10]。Steiger濾波用于解釋暴露和結(jié)局的方差比,明確單個SNP的方向性[11]。通過Phenosanner數(shù)據(jù)庫(http://www.phenoscanner.medschl.cam.ac.uk/)篩查相關(guān)混雜因素[12],排除有混雜關(guān)系的SNPs。

    1.4 因果效應(yīng)分析 逆方差加權(quán)(inverse variance weighting,IVW)法作為主要分析方法來評估4種飲食成分相對攝入量與KSD之間的關(guān)系,其結(jié)果在IVs不存在水平多效性的情況下無偏倚[13]。其中水平多效性通過MR-Egger法進(jìn)行檢驗(yàn),此方法在計算分析中考慮了截距項的存在,當(dāng)截距項恰好為零時,MR-Egger法的結(jié)果等同于IVW法,本研究用MR-Egger法作為補(bǔ)充方法以評估IVW法結(jié)果的穩(wěn)健性。由于加權(quán)中位數(shù)法(weighted median method,WMM)在高達(dá)50%對結(jié)局不產(chǎn)生影響的IVs中,仍可以得出更接近真實(shí)值的估計結(jié)果,故本研究也用其作為IVW法的一個有效補(bǔ)充方法[11]。

    1.5 質(zhì)量控制 使用R軟件(版本4.2.0)進(jìn)行分析,R包為“Two-sample MR”和“MR-PRESSO”。MR多效性殘差和離群值(MR-pleiotropy residual sum and outlier,MR-PRESSO)方法可以用來檢測異常值,并校正水平多效性[14]。采用Cochran’s Q檢驗(yàn)對異質(zhì)性進(jìn)行量化,P<0.05表明分析結(jié)果存在顯著的異質(zhì)性[15],通過散點(diǎn)圖將暴露與結(jié)局SNP的聯(lián)系可視化。MR-Egger回歸用于分析評價無效IVs和水平多效性效應(yīng)引起的偏倚,如果截距項的Pgt;0.05,則表明所選IVs不存在水平多效性[16]。通過Steiger tests判斷暴露與結(jié)局因果關(guān)系的方向[11]。此外,應(yīng)用“留一法”檢驗(yàn)結(jié)果是否穩(wěn)健,通過逐個剔除SNPs后計算剩余SNPs的合并效應(yīng),以評估單個SNP對暴露和結(jié)局變量之間關(guān)聯(lián)的影響[17]

    2 結(jié) 果

    2.1 IVs的提取 通過嚴(yán)格把控篩選標(biāo)準(zhǔn),并且排除混雜因素相關(guān)的SNPs,最終篩選出18個SNP。具體信息見表1,P值均遠(yuǎn)小于5×10-8,說明這些SNPs與暴露因素高度相關(guān)。r2均小于0.001,說明SNPs之間不具有LD。F值均大于10,說明納入的SNP均不是弱IVs。

    2.2 4種常量營養(yǎng)素相對攝入量與KDS因果關(guān)系的MR分析結(jié)果 使用IVW、WMM、MR-Egger、MR-PRESSO法逐一分析4種常量營養(yǎng)素相對攝入量與KSD的因果關(guān)系,結(jié)果見圖2。IVW法分析顯示碳水化合物相對攝入量(OR=2.01,95%CI=1.02~3.96,P=0.042),添加糖相對攝入量(OR=2.07,95%CI=1.00~4.29,P=0.049)是KSD的危險因素。蛋白質(zhì)相對攝入量(OR=0.70,95%CI=0.38~1.29,P=0.249),脂肪相對攝入量(OR=0.79,95%CI=0.27~2.35,P=0.671)與KSD無顯著因果關(guān)系。

    2.3 結(jié)果可視化 運(yùn)用TwoSample包可視化得出相關(guān)SNPs對KSD因果影響的散點(diǎn)圖(圖3),圖上的每一個點(diǎn)代表一個SNP,每個點(diǎn)上的線實(shí)際反映的是95%CI,彩色的線表示MR的擬合結(jié)果。

    從圖3C、D可以得出,IVW彩色線隨著暴露的升高,結(jié)局的發(fā)病風(fēng)險也在升高。MR可視化森林圖(圖2)中的每一條水平實(shí)線反映的是對應(yīng)的暴露對象用不同方法估計出來的結(jié)果:從圖中可以得出添加糖、碳水化合物實(shí)線完全在紅色虛線的右邊,說明由這兩種暴露對象增加能升高結(jié)局的發(fā)病風(fēng)險。綜上兩種可視化圖可以得出,添加糖、碳水化合物相對攝入量的增加會提高KSD的發(fā)病風(fēng)險。

    2.4 異質(zhì)性分析 上述有意義的兩個分析結(jié)果經(jīng)MR-PRESSO Global 檢驗(yàn)P均 gt; 0.05(碳水化合物P=0.746、添加糖P=0.645),表明IVs不存在水平多效性,未發(fā)現(xiàn)離群值。Cochran’s Q檢驗(yàn)(碳水化合物IVW:P=0.519,MR-Egger:P=0.501;添加糖IVW:P=0.549,MR-Egger:P=0.356)表明結(jié)果無異質(zhì)性。MR-Egger截距試驗(yàn)結(jié)果無統(tǒng)計學(xué)意義(碳水化合物:截距=0.038,P=0.404;添加糖:截距=-0.010,P=0.850),表明沒有觀察到多效性。Steiger檢驗(yàn)方向分析示暴露工具與結(jié)局保持正確的因果關(guān)系方向(表2)。漏斗圖由于分析所篩選的SNP數(shù)量較少,無法有效觀察是否基本對稱。

    2.5 敏感性分析 對分析結(jié)果有意義的IVs的敏感性進(jìn)行“留一法”分析如圖4所示,每個黑點(diǎn)代表1個SNP,均位于0的右側(cè),表示無論去除哪個SNP都不會對結(jié)果產(chǎn)生根本影響,表明結(jié)果是穩(wěn)健的。

    3 討 論

    KSD患病率逐年上升,飲食習(xí)慣在疾病發(fā)生發(fā)展中發(fā)揮著一定的作用,通過控制飲食可以有效減少結(jié)石的產(chǎn)生率及復(fù)發(fā)率[18]。在本項MR分析中,IVs從大規(guī)模、非重疊的GWAS的匯總統(tǒng)計中提取。通過使用遺傳基因替代碳水化合物、添加糖、脂肪和蛋白質(zhì)相對攝入量作為暴露工具,結(jié)局為KSD,應(yīng)用兩樣本MR分析分別評估上述4種飲食相對攝入量與KSD風(fēng)險的潛在因果關(guān)系。本研究結(jié)果表明,碳水化合物及添加糖的相對攝入量與KSD風(fēng)險呈正相關(guān),且通過對有意義的結(jié)果進(jìn)行異質(zhì)性、多效性的評估,結(jié)果不存在明顯偏差。

    先前有隨機(jī)對照研究報道,KSD患者和對照組的碳水化合物攝入量相似[19-20]。另一項健康人群與特發(fā)性鈣結(jié)石的對照試驗(yàn)表明結(jié)石組碳水化合物攝入量更高[21]。BARGHOUTHY等[22]在一篇關(guān)于營養(yǎng)與KSD的綜述中總結(jié)了碳水化合物攝入量的改變可通過多種途徑影響KSD的發(fā)生風(fēng)險。TRINCHIERI等[23]的一項回顧性研究發(fā)現(xiàn)單純尿酸KSD患者以及混合尿酸/草酸鈣KSD患者的碳水化合物攝入量顯著高于健康對照組。一項前瞻性研究報道年輕女性蔗糖攝入量與結(jié)石形成的風(fēng)險呈正相關(guān)[24]。另一項關(guān)于年輕女性的大型流行病學(xué)研究也提到,高蔗糖飲食攝入量會提高KSD的發(fā)病率及復(fù)發(fā)率[25],然而在男性的觀察性研究中未發(fā)現(xiàn)兩者存在關(guān)聯(lián)[13,24]。NGUYEN等[26]研究發(fā)現(xiàn)簡單碳水化合物食用負(fù)荷增加會導(dǎo)致尿液中鈣鹽、草酸鹽的增加,精制碳水化合物攝入量的增加同樣會提高尿液中鈣排泄量水平[23],這是草酸鈣結(jié)石復(fù)發(fā)的主要尿液決定因素[27]。LEMANN等[28]的隨機(jī)對照研究發(fā)現(xiàn),健康對照組與草酸鈣結(jié)石患者在口服100 g葡萄糖或蔗糖后,尿鈣排泄率有不同程度的升高。蔗糖攝入增加會提高人體內(nèi)胰島素水平,轉(zhuǎn)而抑制腎臟對鈣的重吸收,可能是內(nèi)在的機(jī)制之一[29]

    飲食中的添加糖主要成分為果糖,果糖攝入量增加顯著提高了結(jié)石形成的風(fēng)險[30],這與本研究結(jié)果一致。果糖可以通過影響尿酸代謝從而增加結(jié)石形成的風(fēng)險[31],這可能是因?yàn)楣谴x顯著抑制了三磷酸腺苷(adenosine triphosphate,ATP)的再生,使尿酸產(chǎn)生的前體成分—磷酸腺苷(adenosine monophosphate,AMP)含量升高,從而進(jìn)一步促進(jìn)了嘌呤合成和尿酸水平的提高[32]。NGUYEN等[26]在研究中也提到,果糖的攝入增加同樣會提高尿液中草酸鹽及尿鈣水平。此外,提高果糖攝入量會增加血管升壓素的釋放[33-34],血清和尿液滲透壓升高,排尿量隨之減少,同樣增加了KSD的發(fā)生風(fēng)險。綜上所述,高碳水化合物攝入以及膳食中過多的添加糖會提高KSD的發(fā)生風(fēng)險。

    蛋白質(zhì)、脂肪同樣在飲食結(jié)構(gòu)中占據(jù)不小的比例,但在本研究中未發(fā)現(xiàn)蛋白質(zhì)相對攝入量、脂肪相對攝入量與KSD存在關(guān)聯(lián)。一項大型流行病學(xué)研究報道,高蛋白攝入會增加男性KSD發(fā)病風(fēng)險或復(fù)發(fā)風(fēng)險[35]。然而,PEDERSEN等[36]的綜述認(rèn)為蛋白質(zhì)攝入與腎結(jié)石風(fēng)險之間的關(guān)系尚無定論。ASOUDEH等[37]通過一項薈萃分析得出總蛋白質(zhì)攝入量包括動物蛋白與KSD的發(fā)生風(fēng)險之間不存在顯著關(guān)聯(lián)。TAYLOR等[38]卻認(rèn)為動物蛋白攝入量僅在BMI<25的男性人群中與KSD形成存在因果關(guān)系,其相關(guān)原因可能是尿液鈣平衡失調(diào)、尿液酸堿度降低以及尿液排出的其他成分變化[39]。很少有報道提到乳制品蛋白及植物蛋白攝入增加會增加KSD的風(fēng)險[40]??紤]到本研究未將蛋白質(zhì)攝入的來源詳細(xì)分類,可能是導(dǎo)致結(jié)果陰性的一個原因。關(guān)于脂肪攝入量和尿路結(jié)石形成風(fēng)險之間因果聯(lián)系的研究較少,而且結(jié)論不一致。兩項隨機(jī)對照研究顯示結(jié)石患者和對照組脂肪攝入量無明顯區(qū)別[19-20]。NAYA等[41]通過分析58例特發(fā)性草酸鈣結(jié)石患者的24 h尿樣和膳食記錄,結(jié)果表明膳食中脂肪酸的含量與尿草酸排泄量顯著相關(guān),可能會增加草酸鹽結(jié)石的發(fā)生率。

    本研究的優(yōu)點(diǎn)在于基于大樣本人群,首次系統(tǒng)地從遺傳學(xué)角度探討不同飲食成分相對攝入量與KSD發(fā)生的潛在因果關(guān)聯(lián)。通過運(yùn)用遺傳變異作為IVs,有效降低了傳統(tǒng)觀察性研究中混雜因素和反向因果關(guān)系的影響,并且研究中使用的暴露和結(jié)局?jǐn)?shù)據(jù)樣本沒有重疊,確保了研究結(jié)果的準(zhǔn)確性。此外還采用了IVW、WMM、MR-Egger回歸分析及MR-PRESSO法等多種統(tǒng)計學(xué)方法驗(yàn)證了結(jié)果的穩(wěn)定性。但是該研究也存在一定的局限性,例如①為了減少人口分層的風(fēng)險,納入的GWAS數(shù)據(jù)來源于歐洲人群,所得結(jié)果和結(jié)論可能不適用于其他種族人群。因此,有必要在其他種族人群中進(jìn)一步研究。②由于飲食數(shù)據(jù)本身的原因,4項MR分析納入的SNP數(shù)量較少,并且研究中P未進(jìn)行校正,可能會導(dǎo)致統(tǒng)計學(xué)效能較低,得出的結(jié)論存在偏倚。但是本研究通過嚴(yán)格的標(biāo)準(zhǔn)篩選了IVs變量,保證其都是強(qiáng)IVs,且研究結(jié)果不存在水平多效性和異質(zhì)性,為研究結(jié)論的可靠性提供了理論支撐。③本次研究未將暴露與結(jié)局IVs進(jìn)一步詳細(xì)分類,例如,膳食攝入的蛋白質(zhì)包括動物蛋白、植物蛋白,脂肪攝入分為飽和脂肪酸與不飽和脂肪酸,而KSD包括鈣鹽結(jié)石、尿酸結(jié)石等不同亞組,未來仍需要更進(jìn)一步的研究。

    總而言之,本研究發(fā)現(xiàn)了強(qiáng)有力的遺傳證據(jù),證明碳水化合物相對攝入量、添加糖相對攝入量是增加KSD發(fā)生風(fēng)險的危險性因素,蛋白相對攝入量、脂肪相對攝入量與KSD發(fā)生風(fēng)險無明顯相關(guān)性,該結(jié)論對KSD的飲食預(yù)防策略有重要意義,值得在更全面的隨機(jī)臨床試驗(yàn)中進(jìn)行驗(yàn)證。

    參考文獻(xiàn):

    [1] RAJA A,HEKMATI Z,JOSHI HB.How do urinary calculi influence health-related quality of life and patient treatment preference:a systematic review[J].J Endourol,2016,30(7):727-743.

    [2] GERAGHTY RM,COOK P,WALKER V,et al.Evaluation of the economic burden of kidney stone disease in the UK:a retrospective cohort study with a mean follow-up of 19 years[J].BJU Int,2020,125(4):586-594.

    [3] FERRARO PM,TAYLOR EN,GAMBARO G,et al.Dietary and lifestyle risk factors associated with incident kidney stones in men and women[J].J Urol,2017,198(4):858-863.

    [4] SKRIVANKOVA VW,RICHMOND RC,WOOLF BAR,et al.Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR):explanation and elaboration[J].BMJ,2021,375:n2233.

    [5] BURGESS S,BUTTERWORTH A,THOMPSON SG.Mendelian randomization analysis with multiple genetic variants using summarized data[J].Genet Epidemiol,2013,37(7):658-665.

    [6] KURKI MI,KARJALAINEN J,PALTA P,et al.FinnGen provides genetic insights from a well-phenotyped isolated population[J].Nature,2023,613(7944):508-518.

    [7] MEDDENS SFW,DE VLAMING R,BOWERS P,et al.Genomic analysis of diet composition finds novel loci and associations with health and lifestyle[J].Mol Psychiatry,2021,26(6):2056-2069.

    [8] DAVEY SMITH G,HEMANI G.Mendelian randomization:genetic anchors for causal inference in epidemiological studies[J].Hum Mol Genet,2014,23(R1):R89-R98.

    [9] YAVORSKA OO,BURGESS S.Mendelian randomization:an R package for performing Mendelian randomization analyses using summarized data[J].Int J Epidemiol,2017,46(6):1734-1739.

    [10] PRITCHARD JK,PRZEWORSKI M.Linkage disequilibrium in humans:models and data[J].Am J Hum Genet,2001,69(1):1-14.

    [11] BOWDEN J,DAVEY SMITH G,BURGESS S.Mendelian randomization with invalid instruments:effect estimation and bias detection through Egger regression[J].Int J Epidemiol,2015,44(2):512-525.

    [12] KAMAT MA,BLACKSHAW JA,YOUNG R,et al.PhenoScanner V2:an expanded tool for searching human genotype-phenotype associations[J].Bioinformatics,2019,35(22):4851-4853.

    [13] BURGESS S,DAVEY SMITH G,DAVIES NM,et al.Guidelines for performing Mendelian randomization investigations:update for summer 2023[J].Wellcome Open Res,2019,4:186.

    [14] ONG JS,MACGREGOR S.Implementing MR-PRESSO and GCTA-GSMR for pleiotropy assessment in Mendelian randomization studies from a practitioner’s perspective[J].Genet Epidemiol,2019,43(6):609-616.

    [15] DAVIES NM,HOLMES MV,DAVEY SMITH G.Reading Mendelian randomisation studies:a guide,glossary,and checklist for clinicians[J].BMJ,2018,362:k601.

    [16] BURGESS S,THOMPSON SG.Interpreting findings from Mendelian randomization using the MR-Egger method[J].Eur J Epidemiol,2017,32(5):377-389.

    [17] RUSHING C,BULUSU A,HURWITZ HI,et al.A leave-one-out cross-validation SAS macro for the identification of markers associated with survival[J].Comput Biol Med,2015,57:123-129.

    [18] BORGHI L,MESCHI T,MAGGIORE U,et al.Dietary therapy in idiopathic nephrolithiasis[J].Nutr Rev,2006,64(7 Pt 1):301-312.

    [19] POWER C,BARKER DJ,NELSON M,et al.Diet and renal stones:a case-control study[J].Br J Urol,1984,56(5):456-459.

    [20] FELLSTR?M B,DANIELSON BG,KARLSTR?M B,et al.Dietary habits in renal stone patients compared with healthy subjects[J].Br J Urol,1989,63(6):575-580.

    [21] AL ZAHRANI H,NORMAN RW,THOMPSON C,et al.The dietary habits of idiopathic calcium stone-formers and normal control subjects[J].BJU Int,2000,85(6):616-20.

    [22] BARGHOUTHY Y,CORRALES M,SOMANI B.The relationship between modern fad diets and kidney stone disease:a systematic review of literature[J].Nutrients,2021,13(12):e186-e195.

    [23] TRINCHIERI A,MONTANARI E.Biochemical and dietary factors of uric acid stone formation[J].Urolithiasis,2018,46(2):167-172.

    [24] CURHAN GC,WILLETT WC,KNIGHT EL,et al.Dietary factors and the risk of incident kidney stones in younger women:Nurses’ Health Study Ⅱ[J].Arch Intern Med,2004,164(8):885-891.

    [25] CURHAN GC,WILLETT WC,SPEIZER FE,et al.Comparison of dietary calcium with supplemental calcium and other nutrients as factors affecting the risk for kidney stones in women[J].Ann Int Med,1997,126(7):497-504.

    [26] NGUYEN NU,DUMOULIN G,HENRIET MT,et al.Carbohydrate metabolism and urinary excretion of calcium and oxalate after ingestion of polyol sweeteners[J].J Clin Endocrinol Metab,1993,77(2):388-392.

    [27] SIENER R,GLATZ S,NICOLAY C,et al.Prospective study on the efficacy of a selective treatment and risk factors for relapse in recurrent calcium oxalate stone patients[J].Eur Urol,2003,44(4):467-474.

    [28] LEMANN JJR,PIERING WF,LENNON EJ.Possible role of carbohydrate-induced calciuria in calcium oxalate kidney-stone formation[J].N Engl J Med,1969,280(5):232-237.

    [29] HOLL MG,ALLEN LH.Sucrose ingestion,insulin response and mineral metabolism in humans[J].J Nutr,1987,117(7):1229-1233.

    [30] LIN BB,LIN ME,HUANG RH,et al.Dietary and lifestyle factors for primary prevention of nephrolithiasis:a systematic review and meta-analysis[J].BMC Nephrol,2020,21(1):267.

    [31] JOHNSON RJ,PEREZ-POZO SE,LILLO JL,et al.Fructose increases risk for kidney stones:potential role in metabolic syndrome and heat stress[J].BMC Nephrol,2018,19(1):315.

    [32] RHO YH,ZHU Y,CHOI HK.The epidemiology of uric acid and fructose[J].Semin Nephrol,2011,31(5):410-419.

    [33] SONG Z,RONCAL-JIMENEZ CA,LANASPA-GARCIA MA,et al.Role of fructose and fructokinase in acute dehydration-induced vasopressin gene expression and secretion in mice[J].J Neurophysiol,2017,117(2):646-654.

    [34] GARCíA-ARROYO FE,CRISTóBAL M,ARELLANO-BUENDíA AS,et al.Rehydration with soft drink-like beverages exacerbates dehydration and worsens dehydration-associated renal injury[J].Am J Physiol Regul Integr Comp Physiol,2016,311(1):R57-R65.

    [35] CURHAN GC,WILLETT WC,RIMM EB,et al.A prospective study of dietary calcium and other nutrients and the risk of symptomatic kidney stones[J].N Engl J Med,1993,328(12):833-838.

    [36] PEDERSEN AN,KONDRUP J,B?RSHEIM E.Health effects of protein intake in healthy adults:a systematic literature review[J].Food Nutr Res,2013,30:57.

    [37] ASOUDEH F,TALEBI S,JAYEDI A,et al.Associations of total protein or animal protein intake and animal protein sources with risk of kidney stones:a systematic review and dose-response meta-analysis[J].Adv Nutr,2022,13(3):821-832.

    [38] TAYLOR EN,STAMPFER MJ,CURHAN GC.Dietary factors and the risk of incident kidney stones in men:new insights after 14 years of follow-up[J].J Am Soc Nephrol,2004,15(12):3225-3232.

    [39] REDDY ST,WANG CY,SAKHAEE K,et al.Effect of low-carbohydrate high-protein diets on acid-base balance,stone-forming propensity,and calcium metabolism[J].Am J Kidney Dis,2002,40(2):265-274.

    [40] FERRARO PM,BARGAGLI M,TRINCHIERI A,et al.Risk of kidney stones:influence of dietary factors,dietary patterns,and vegetarian-vegan diets[J].Nutrients,2020,12(3):3561-3571.

    [41] NAYA Y,ITO H,MASAI M,et al.Association of dietary fatty acids with urinary oxalate excretion in calcium oxalate stone-formers in their fourth decade[J].BJU Int,2002,89(9):842-846.

    (編輯 郭楚君)

    黑人猛操日本美女一级片| 悠悠久久av| 久久这里只有精品19| 巨乳人妻的诱惑在线观看| 亚洲一区二区三区欧美精品| 又紧又爽又黄一区二区| 久久精品91无色码中文字幕| 亚洲国产欧美日韩在线播放| 法律面前人人平等表现在哪些方面| 国产av一区二区精品久久| 国产精品久久久久成人av| 久久久精品国产亚洲av高清涩受| 在线观看一区二区三区激情| 国产老妇伦熟女老妇高清| 久久影院123| 欧美日韩av久久| 99国产精品一区二区三区| 天天操日日干夜夜撸| 欧美黄色片欧美黄色片| 在线 av 中文字幕| 成年动漫av网址| 亚洲精品av麻豆狂野| 亚洲成人免费电影在线观看| 精品国产乱码久久久久久男人| 一级黄色大片毛片| 咕卡用的链子| 色在线成人网| 国产极品粉嫩免费观看在线| 涩涩av久久男人的天堂| 日日夜夜操网爽| 久久亚洲真实| 亚洲av美国av| 亚洲avbb在线观看| 欧美日韩精品网址| 十八禁人妻一区二区| 天天躁日日躁夜夜躁夜夜| 成人免费观看视频高清| 在线观看免费高清a一片| 十八禁人妻一区二区| 久热这里只有精品99| 成人国语在线视频| 99久久人妻综合| 亚洲欧美一区二区三区久久| 麻豆av在线久日| 国产真人三级小视频在线观看| 免费久久久久久久精品成人欧美视频| 久久狼人影院| 亚洲 欧美一区二区三区| 一级片免费观看大全| 精品一区二区三区视频在线观看免费 | 成人av一区二区三区在线看| 国产在线精品亚洲第一网站| 亚洲精品美女久久av网站| 激情视频va一区二区三区| 人妻久久中文字幕网| 久久天躁狠狠躁夜夜2o2o| 国产91精品成人一区二区三区 | 亚洲伊人色综图| 这个男人来自地球电影免费观看| 国产精品国产高清国产av | 一级毛片女人18水好多| 99精品久久久久人妻精品| 日韩熟女老妇一区二区性免费视频| 欧美亚洲日本最大视频资源| 久久青草综合色| 欧美日韩一级在线毛片| 黄频高清免费视频| 好男人电影高清在线观看| 一区二区三区国产精品乱码| 日韩欧美三级三区| 大片免费播放器 马上看| 99久久精品国产亚洲精品| 亚洲国产毛片av蜜桃av| 无限看片的www在线观看| 日本a在线网址| 精品国产一区二区三区四区第35| 一区二区av电影网| 少妇被粗大的猛进出69影院| 不卡一级毛片| 免费观看a级毛片全部| 波多野结衣av一区二区av| 无限看片的www在线观看| 日韩中文字幕视频在线看片| 成人亚洲精品一区在线观看| 午夜福利视频精品| 午夜激情av网站| 中文字幕人妻熟女乱码| 又紧又爽又黄一区二区| 黄片播放在线免费| 国产精品欧美亚洲77777| 国产在线视频一区二区| 99国产精品免费福利视频| 大型黄色视频在线免费观看| 热99国产精品久久久久久7| 99精品在免费线老司机午夜| 欧美中文综合在线视频| 中文字幕人妻丝袜制服| 1024视频免费在线观看| 黑人欧美特级aaaaaa片| 色婷婷av一区二区三区视频| 日本a在线网址| 日韩欧美三级三区| 国产高清激情床上av| 成人av一区二区三区在线看| 麻豆成人av在线观看| 一本久久精品| 国产1区2区3区精品| 黄片播放在线免费| 丰满人妻熟妇乱又伦精品不卡| 亚洲综合色网址| 桃红色精品国产亚洲av| 人成视频在线观看免费观看| 老汉色∧v一级毛片| 精品高清国产在线一区| 国产成人啪精品午夜网站| 欧美黑人欧美精品刺激| 在线观看人妻少妇| 自线自在国产av| 激情视频va一区二区三区| 久久久精品94久久精品| 汤姆久久久久久久影院中文字幕| 欧美精品人与动牲交sv欧美| 最近最新免费中文字幕在线| 中国美女看黄片| 满18在线观看网站| 欧美黑人欧美精品刺激| 国产一区二区三区视频了| 中文字幕色久视频| 久久精品91无色码中文字幕| 十八禁网站网址无遮挡| 在线观看一区二区三区激情| 国产亚洲一区二区精品| 757午夜福利合集在线观看| 亚洲精品国产区一区二| 黑人操中国人逼视频| 制服人妻中文乱码| 亚洲avbb在线观看| 99久久国产精品久久久| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美一区二区三区黑人| 男女之事视频高清在线观看| 蜜桃国产av成人99| 亚洲一区中文字幕在线| 国产精品久久久久成人av| 男男h啪啪无遮挡| 日日摸夜夜添夜夜添小说| 欧美另类亚洲清纯唯美| 中文字幕色久视频| 99热国产这里只有精品6| 18禁美女被吸乳视频| 亚洲熟妇熟女久久| 欧美成人午夜精品| 成人18禁高潮啪啪吃奶动态图| 天天操日日干夜夜撸| 人人妻,人人澡人人爽秒播| 欧美成人午夜精品| 老司机在亚洲福利影院| 少妇 在线观看| 日日夜夜操网爽| 在线播放国产精品三级| 免费久久久久久久精品成人欧美视频| 国产97色在线日韩免费| 久久精品国产亚洲av高清一级| 国产成人精品无人区| 蜜桃在线观看..| 变态另类成人亚洲欧美熟女 | 黄片大片在线免费观看| 狠狠精品人妻久久久久久综合| 黄色成人免费大全| 变态另类成人亚洲欧美熟女 | 欧美黄色淫秽网站| 久久精品91无色码中文字幕| 成在线人永久免费视频| 99国产精品一区二区蜜桃av | 动漫黄色视频在线观看| 国产成人欧美| 18禁观看日本| 色综合婷婷激情| 精品少妇内射三级| 757午夜福利合集在线观看| 搡老熟女国产l中国老女人| 王馨瑶露胸无遮挡在线观看| 最新美女视频免费是黄的| 欧美中文综合在线视频| 女警被强在线播放| 中文欧美无线码| 国产色视频综合| 黄片大片在线免费观看| 午夜精品久久久久久毛片777| 99香蕉大伊视频| 成年人黄色毛片网站| 欧美在线一区亚洲| 制服诱惑二区| 午夜老司机福利片| 又大又爽又粗| 午夜日韩欧美国产| 日本wwww免费看| 18在线观看网站| 91精品三级在线观看| 99精国产麻豆久久婷婷| 80岁老熟妇乱子伦牲交| 日韩一卡2卡3卡4卡2021年| 成人18禁在线播放| 日韩制服丝袜自拍偷拍| 亚洲色图av天堂| 成年版毛片免费区| 中文字幕av电影在线播放| 亚洲欧美日韩高清在线视频 | 久久精品国产99精品国产亚洲性色 | 别揉我奶头~嗯~啊~动态视频| 久久性视频一级片| 黄片小视频在线播放| 亚洲色图 男人天堂 中文字幕| cao死你这个sao货| 搡老岳熟女国产| 丁香六月天网| 免费黄频网站在线观看国产| 91麻豆av在线| 美女扒开内裤让男人捅视频| 中文字幕色久视频| 久久久久精品国产欧美久久久| 国产成人av激情在线播放| 99国产精品99久久久久| 亚洲av片天天在线观看| 色综合婷婷激情| 亚洲欧美一区二区三区久久| 热99国产精品久久久久久7| 亚洲午夜精品一区,二区,三区| 亚洲自偷自拍图片 自拍| 又黄又粗又硬又大视频| 国产精品国产高清国产av | 女人久久www免费人成看片| 午夜福利欧美成人| 性高湖久久久久久久久免费观看| 国产成人精品无人区| 757午夜福利合集在线观看| 91av网站免费观看| 高清av免费在线| 中文亚洲av片在线观看爽 | 精品久久久精品久久久| 波多野结衣av一区二区av| 天天添夜夜摸| 亚洲av成人不卡在线观看播放网| 99re6热这里在线精品视频| 大型av网站在线播放| 欧美精品人与动牲交sv欧美| 亚洲七黄色美女视频| 一级,二级,三级黄色视频| 在线观看免费日韩欧美大片| av一本久久久久| 怎么达到女性高潮| 青青草视频在线视频观看| 五月天丁香电影| 99re在线观看精品视频| 一级片免费观看大全| av又黄又爽大尺度在线免费看| 成人永久免费在线观看视频 | 最黄视频免费看| 伦理电影免费视频| 亚洲av成人一区二区三| 纯流量卡能插随身wifi吗| videos熟女内射| av欧美777| 一个人免费在线观看的高清视频| 熟女少妇亚洲综合色aaa.| 黄色成人免费大全| 亚洲成av片中文字幕在线观看| 精品国产一区二区三区久久久樱花| 热99re8久久精品国产| 91精品国产国语对白视频| 亚洲成人免费av在线播放| 久久中文字幕人妻熟女| 国产av又大| 99久久精品国产亚洲精品| 欧美 亚洲 国产 日韩一| 天天躁狠狠躁夜夜躁狠狠躁| 这个男人来自地球电影免费观看| 黄片播放在线免费| 欧美 日韩 精品 国产| 亚洲中文av在线| 极品人妻少妇av视频| 精品卡一卡二卡四卡免费| 美女国产高潮福利片在线看| 久久天堂一区二区三区四区| 国产av又大| 国产1区2区3区精品| 亚洲国产欧美日韩在线播放| 国产成人影院久久av| 又黄又粗又硬又大视频| 亚洲精品久久成人aⅴ小说| av超薄肉色丝袜交足视频| 黑人猛操日本美女一级片| 亚洲欧美日韩另类电影网站| 人妻一区二区av| 国产成人影院久久av| 久久国产精品人妻蜜桃| 久久久久久久国产电影| 亚洲av国产av综合av卡| 国产成人免费无遮挡视频| 亚洲精品在线美女| 一个人免费在线观看的高清视频| 亚洲精品乱久久久久久| 一二三四在线观看免费中文在| 热99久久久久精品小说推荐| 狂野欧美激情性xxxx| 无限看片的www在线观看| 老熟女久久久| 大香蕉久久网| 亚洲精品国产色婷婷电影| 国产精品98久久久久久宅男小说| 国产精品1区2区在线观看. | 国产精品二区激情视频| 男人操女人黄网站| 老司机亚洲免费影院| 亚洲熟妇熟女久久| 欧美乱妇无乱码| 精品熟女少妇八av免费久了| 中文字幕最新亚洲高清| av天堂在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品电影一区二区三区 | 91字幕亚洲| 男女边摸边吃奶| av国产精品久久久久影院| 午夜免费成人在线视频| 久久九九热精品免费| 在线播放国产精品三级| 精品一品国产午夜福利视频| 男人操女人黄网站| 成人永久免费在线观看视频 | 日韩大码丰满熟妇| 热re99久久国产66热| 少妇精品久久久久久久| 久久精品91无色码中文字幕| 日本一区二区免费在线视频| 欧美中文综合在线视频| 亚洲中文字幕日韩| 午夜福利,免费看| 精品国产乱码久久久久久小说| 国产成人精品久久二区二区免费| 老熟妇仑乱视频hdxx| 国产亚洲精品久久久久5区| 脱女人内裤的视频| 成人国产一区最新在线观看| 久9热在线精品视频| 国产精品欧美亚洲77777| 国产在视频线精品| 日韩成人在线观看一区二区三区| 亚洲七黄色美女视频| 欧美中文综合在线视频| 欧美亚洲日本最大视频资源| 欧美激情高清一区二区三区| 中文字幕色久视频| 怎么达到女性高潮| 国产又色又爽无遮挡免费看| 黄色视频不卡| 亚洲精品国产一区二区精华液| 国产视频一区二区在线看| 啦啦啦 在线观看视频| 亚洲精品乱久久久久久| 精品少妇一区二区三区视频日本电影| 久久中文看片网| 如日韩欧美国产精品一区二区三区| 淫妇啪啪啪对白视频| 丁香欧美五月| 日本av免费视频播放| 亚洲av成人一区二区三| 1024香蕉在线观看| 黄片大片在线免费观看| 国产精品电影一区二区三区 | 久久精品国产99精品国产亚洲性色 | 老司机午夜福利在线观看视频 | 99久久人妻综合| 亚洲男人天堂网一区| 国产亚洲午夜精品一区二区久久| 午夜两性在线视频| 久久精品成人免费网站| 99久久99久久久精品蜜桃| 又紧又爽又黄一区二区| 久久婷婷成人综合色麻豆| 中文字幕精品免费在线观看视频| 一二三四在线观看免费中文在| 天天躁狠狠躁夜夜躁狠狠躁| 精品亚洲成国产av| 久久av网站| 少妇裸体淫交视频免费看高清 | 不卡一级毛片| 伊人久久大香线蕉亚洲五| 国产av国产精品国产| 女同久久另类99精品国产91| 午夜福利在线免费观看网站| 久久久国产一区二区| 999久久久精品免费观看国产| 亚洲人成77777在线视频| 黑人欧美特级aaaaaa片| 国产一区有黄有色的免费视频| 亚洲天堂av无毛| 欧美 亚洲 国产 日韩一| 老汉色av国产亚洲站长工具| 成人黄色视频免费在线看| 99久久99久久久精品蜜桃| 久久久国产欧美日韩av| av网站免费在线观看视频| 伦理电影免费视频| 飞空精品影院首页| 90打野战视频偷拍视频| 久久毛片免费看一区二区三区| 国产熟女午夜一区二区三区| 一区二区三区国产精品乱码| 午夜福利欧美成人| 欧美精品一区二区免费开放| 一级毛片精品| 国产精品香港三级国产av潘金莲| 成人黄色视频免费在线看| 中文字幕精品免费在线观看视频| 蜜桃在线观看..| 免费不卡黄色视频| 日韩一区二区三区影片| 啦啦啦 在线观看视频| 国产精品亚洲av一区麻豆| 日本五十路高清| 国产精品亚洲av一区麻豆| videos熟女内射| 亚洲成a人片在线一区二区| 日本五十路高清| 老司机午夜十八禁免费视频| 自线自在国产av| 国产欧美日韩综合在线一区二区| 91字幕亚洲| 精品一区二区三区四区五区乱码| 波多野结衣av一区二区av| 欧美黄色淫秽网站| 男女无遮挡免费网站观看| 国产精品欧美亚洲77777| 国产精品国产av在线观看| 久久国产精品人妻蜜桃| 久久久水蜜桃国产精品网| 欧美久久黑人一区二区| 亚洲精品国产精品久久久不卡| 99国产精品99久久久久| 少妇裸体淫交视频免费看高清 | 天天添夜夜摸| 中文字幕精品免费在线观看视频| 欧美日韩黄片免| 欧美日韩亚洲高清精品| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 久久久久精品国产欧美久久久| 国产不卡一卡二| 制服诱惑二区| 麻豆乱淫一区二区| 69av精品久久久久久 | 欧美黑人精品巨大| 久久这里只有精品19| 国产亚洲一区二区精品| 一夜夜www| 免费一级毛片在线播放高清视频 | 天天躁夜夜躁狠狠躁躁| 啪啪无遮挡十八禁网站| 国产精品久久久久久精品电影小说| 99精国产麻豆久久婷婷| 亚洲av片天天在线观看| videosex国产| 日日摸夜夜添夜夜添小说| 国产亚洲欧美在线一区二区| 午夜福利视频精品| 久久久久精品人妻al黑| 黄色片一级片一级黄色片| 欧美日韩av久久| a级片在线免费高清观看视频| 亚洲精品国产一区二区精华液| 日韩 欧美 亚洲 中文字幕| 亚洲综合色网址| 亚洲国产看品久久| 少妇的丰满在线观看| 精品亚洲乱码少妇综合久久| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩中文字幕欧美一区二区| 超色免费av| 麻豆乱淫一区二区| 久久亚洲真实| 欧美激情久久久久久爽电影 | 亚洲av国产av综合av卡| 757午夜福利合集在线观看| 叶爱在线成人免费视频播放| 午夜日韩欧美国产| videosex国产| 久久 成人 亚洲| 国产午夜精品久久久久久| 亚洲精品久久午夜乱码| 一本一本久久a久久精品综合妖精| 亚洲精品国产色婷婷电影| 一区二区日韩欧美中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 一个人免费看片子| 三级毛片av免费| 99久久人妻综合| 日本wwww免费看| www.自偷自拍.com| bbb黄色大片| 999久久久国产精品视频| 别揉我奶头~嗯~啊~动态视频| 欧美日韩黄片免| 国产熟女午夜一区二区三区| 欧美精品高潮呻吟av久久| av天堂久久9| 亚洲av日韩精品久久久久久密| 日日摸夜夜添夜夜添小说| 人人妻人人澡人人爽人人夜夜| 欧美激情久久久久久爽电影 | 午夜免费成人在线视频| 成人18禁在线播放| 黄色视频不卡| 一边摸一边抽搐一进一出视频| 久久午夜亚洲精品久久| tube8黄色片| 国产日韩一区二区三区精品不卡| 黑丝袜美女国产一区| av免费在线观看网站| 亚洲五月婷婷丁香| 热re99久久精品国产66热6| 国产成人一区二区三区免费视频网站| 国产精品香港三级国产av潘金莲| 在线播放国产精品三级| 五月开心婷婷网| 手机成人av网站| 99精品久久久久人妻精品| 亚洲成人免费av在线播放| 制服人妻中文乱码| 叶爱在线成人免费视频播放| 亚洲三区欧美一区| 日韩大片免费观看网站| 亚洲男人天堂网一区| 一区二区三区乱码不卡18| 这个男人来自地球电影免费观看| 人人妻,人人澡人人爽秒播| 男女高潮啪啪啪动态图| 露出奶头的视频| 国产精品免费视频内射| 国产成人一区二区三区免费视频网站| 精品人妻在线不人妻| 美女福利国产在线| 久久天躁狠狠躁夜夜2o2o| 欧美+亚洲+日韩+国产| 夫妻午夜视频| 一个人免费在线观看的高清视频| 亚洲欧美日韩另类电影网站| 99热网站在线观看| 黄片小视频在线播放| 一本综合久久免费| 欧美成狂野欧美在线观看| 欧美精品高潮呻吟av久久| 黄色a级毛片大全视频| 18禁美女被吸乳视频| 精品国产亚洲在线| 精品久久蜜臀av无| 国产一区二区三区综合在线观看| 色老头精品视频在线观看| 精品国产乱码久久久久久小说| 一区福利在线观看| 精品少妇内射三级| 日韩大码丰满熟妇| 一个人免费看片子| 精品少妇久久久久久888优播| 亚洲av日韩在线播放| 久久久久久久大尺度免费视频| 美女扒开内裤让男人捅视频| 日韩熟女老妇一区二区性免费视频| av免费在线观看网站| 精品久久久久久电影网| 黄片小视频在线播放| 欧美 亚洲 国产 日韩一| 一区在线观看完整版| av在线播放免费不卡| 久久国产亚洲av麻豆专区| 久久精品91无色码中文字幕| 极品少妇高潮喷水抽搐| 青青草视频在线视频观看| 狠狠狠狠99中文字幕| 免费在线观看影片大全网站| 亚洲国产中文字幕在线视频| 丝袜喷水一区| 日本av手机在线免费观看| 亚洲av第一区精品v没综合| 国产高清视频在线播放一区| 午夜成年电影在线免费观看| 黑丝袜美女国产一区| 精品国产亚洲在线| 久久久国产一区二区| 午夜福利免费观看在线| 亚洲精品成人av观看孕妇| 18禁裸乳无遮挡动漫免费视频| 精品人妻熟女毛片av久久网站| 9热在线视频观看99| 极品人妻少妇av视频| 一本一本久久a久久精品综合妖精| 国产精品熟女久久久久浪| 国产精品美女特级片免费视频播放器 | 十八禁人妻一区二区| 亚洲国产毛片av蜜桃av| 亚洲精品av麻豆狂野| 一边摸一边抽搐一进一小说 | 天堂中文最新版在线下载| 在线亚洲精品国产二区图片欧美| 高清在线国产一区| 91av网站免费观看| 在线亚洲精品国产二区图片欧美| 久久久久精品国产欧美久久久| 肉色欧美久久久久久久蜜桃| 久久久久久人人人人人| 正在播放国产对白刺激| 一夜夜www| 丝瓜视频免费看黄片| 国产aⅴ精品一区二区三区波| 在线看a的网站|