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

    通過加權基因共表達網(wǎng)絡和LASSO分析識別兒童呼吸道合胞病毒的診斷基因

    2024-12-31 00:00:00董海麗王國勝谷雙龍宋南萍
    上海醫(yī)藥 2024年19期
    關鍵詞:中性粒細胞

    摘 要 目的:運用生物信息學方法挖掘兒童呼吸道合胞病毒(RSV)的診斷標志物,為RSV的靶向治療提供依據(jù)。方法:從GEO中獲得RSV的基因表達數(shù)據(jù)集。運用差異表達分析、免疫浸潤分析、WGCNA和LASSO確定RSV的關鍵生物標志物。同時,繪制ROC曲線對關鍵基因的診斷性能進行評估。結(jié)果:共識別出592個差異基因和10種差異浸潤的免疫細胞,例如中性粒細胞。WGCNA和LASSO分析篩選出7個關鍵的生物標志物,包括ADM、LTF、UBE2J1、PDK4、S100A4、BST1、CAMP。此外,ROC曲線分析發(fā)現(xiàn)這些基因具有較高的診斷價值。結(jié)論:篩選的7個基因可作為RSV的診斷生物標志物。

    關鍵詞 呼吸道合胞病毒 診斷基因 中性粒細胞 加權基因共表達網(wǎng)絡分析

    中圖分類號:R446.9; R725.1 文獻標志碼:A 文章編號:1006-1533(2024)19-0063-09

    引用本文 董海麗, 王國勝, 谷雙龍, 等. 通過加權基因共表達網(wǎng)絡和LASSO分析識別兒童呼吸道合胞病毒的診斷基因[J]. 上海醫(yī)藥, 2024, 45(19): 63-71.

    Identification of diagnostic genes related to respiratory syncytial virus in children based on WGCNA and LASSO algorithms

    DONG Haili, WANG Guosheng, GU Shuanglong, SONG Nanping

    (Department of Pediatrics, Zhengzhou 460 Hospital, Zhengzhou 450000, China)

    ABSTRACT Objective: The diagnostic biomarkers of respiratory syncytial virus (RSV) in children were explored by bioinformatics analysis to provide basis for targeted therapy of RSV. Methods: Gene expression profiles of RSV were obtained from the GEO database. Differential expression analysis, immune infiltration analysis, WGCNA and LASSO were applied to identify the key biomarkers of RSV. Meanwhile, ROC curve was utilized to evaluate the diagnostic performance of key genes. Results: A total of 592 differentially expressed genes and 10 differential immune cells (such as neutrophils) were identified between RSV and healthy samples. Further, WGCNA and LASSO analysis screened out seven key biomarkers, including ADM, LTF, UBE2J1, PDK4, S100A4, BST1, and CAMP. In addition, ROC curve showed that these genes had high diagnostic value. Conclusion: The identified seven biomarkers may serve as diagnostic biomarkers for RSV.

    KEY WORDS respiratory syncytial virus; diagnostic genes; neutrophils; weighted gene co-expression network analysis

    呼吸道合胞病毒(respiratory syncytial virus,RSV)是全球兒童呼吸道感染的主要原因[1]。據(jù)統(tǒng)計,每年有210萬5歲以下的兒童因RSV感染需要門診治療或住院,并且在全世界28 d至6個月的嬰兒中,RSV疾病的總死亡率負擔為1/28[2-3]。因此,迫切需要有效的RSV預防策略來解決這一重大公共衛(wèi)生問題和負擔。帕利珠單抗廣泛的用于預防嬰兒RSV疾病[4]。盡管其藥效已經(jīng)得到證實,但是尚未引進國內(nèi)臨床應用[5]。迄今為止,RSV疫苗仍然處于開發(fā)階段。此外,診斷RSV感染的有效方法也十分有限。所以,了解RSV感染的分子發(fā)病機制對于指導未來的診斷和治療至關重要。目前,已有研究旨在尋找RSV的診斷標志物,揭示其潛在的分子機制[6-7]。但是,識別RSV診斷生物標志物的相關研究仍然缺乏。本研究利用一系列生物信息學分析旨在探討與RSV中免疫浸潤細胞相關的關鍵生物標志物。本工作將有望為提高RSV患者的診斷以及治療靶點的開發(fā)提供幫助。

    1 材料和方法

    1.1 表達譜數(shù)據(jù)的收集和預處理

    從NCBI GEO數(shù)據(jù)庫中下載6個兒童RSV芯片數(shù)據(jù)集,包括GSE188427、GSE105450、GSE103842、GSE103119、GSE80179、GSE77087。其中,GSE188427包括353個樣本,來自于51個健康志愿者和302個RSV患者的全血樣本,在本研究中作為訓練集。其余數(shù)據(jù)集作為模型驗證數(shù)據(jù)集。另外,使用R語言中的affy包對下載的數(shù)據(jù)進行預處理,并使用rma函數(shù)進行歸一化處理。

    1.2 基因差異表達分析

    通過limma包對健康和RSV患者中的基因進行差異表達分析,計算基因相應的P值和logFC值。同時,利用FDR方法對P值進行多重檢驗校正,得到P-adjust(Benjamini-Hochberg校正)。差異表達基因(differential expressed genes,DEGs)的閾值為:P-adjust<0.05和|logFC|>1.2。

    1.3 GSEA功能分析

    為了探索疾病組和正常組之間的差異調(diào)節(jié)通路,使用GSEA(V4.3.2)對DEGs進行功能富集分析,包括GO-BP、KEGG和Reactome。顯著富集的閾值為P-adjust<0.05。

    1.4 免疫微環(huán)境分析

    基于患者和對照樣本的基因表達數(shù)據(jù),利用CIBERSORT對樣本中22種免疫細胞的比例進行計算。隨后,運用Wilcox檢驗比較了兩組間免疫細胞浸潤水平的差異。

    1.5 WGCNA分析

    使用WGCNA篩選與RSV有關的關鍵基因。通過篩選中位數(shù)絕對偏差(median absolute deviation,MAD)排名前75%且MAD值>0.01的基因進行共表達分析(13 963個基因)。隨后,計算每個模塊與差異的免疫細胞浸潤程度之間的相關性。獲得與免疫細胞最相關的模塊后(又稱為hub模塊),將模塊基因與DEGs取交集,從而得到與免疫細胞浸潤相關的關鍵基因用于后續(xù)分析。

    1.6 關鍵基因的功能富集分析

    利用DAVID 2021對這些基因進行GO功能和KEGG通路富集分析。選取P-adjust<0.05作為顯著富集的閾值。

    1.7 診斷基因篩選

    運用R4.3.1語言中的glmnet包(V4.1-7)從關鍵基因中篩選診斷基因。首先,進行單因素logistic回歸分析,選取P<0.05的基因進行后續(xù)分析。接下來,利用LASSO分析來篩選RSV的診斷基因。LASSO回歸分析在R軟件中運行10次K交叉驗證,用于對包含的變量進行居中和歸一化處理,然后選擇最佳的l值。

    1.8 列線圖的構(gòu)建

    通過rms包繪制包括診斷基因的列線圖來判斷RSV患病的風險率,并且繪制校正曲線來驗證模型的準確性。

    1.9 外部數(shù)據(jù)集驗證

    在5個驗證隊列中,對診斷基因在健康對照和RSV樣本中的表達水平進行比較。隨后,利用pROC包繪制診斷基因的ROC曲線,通過計算AUC值評估基因的預測性能。

    2 結(jié)果

    2.1 DEGs篩選和GSEA富集分析

    在RSV和對照中共篩選到592個DEGs,包括453個上調(diào)和139個下調(diào)基因。GSEA結(jié)果揭示這些基因主要富集在380個GO-BP類目,57個KEGG通路和149個Reactome通路。圖1展示了NES排名前5的結(jié)果。

    2.2 RSV和健康對照患者呈現(xiàn)不同的免疫景觀

    10種免疫細胞的浸潤水平在兩組中有顯著差異(圖2A和2B)。在RSV樣本中,活化樹突狀細胞(activated dendritic cells)、M0型巨噬細胞(M0 macrophages)、活化肥大細胞(activated mast cells)、中性粒細胞(neutrophils)、活化CD4記憶T細胞(activated T cells CD4 memory)的豐度遠高于對照樣本,而native B細胞(native B cells)、嗜酸性粒細胞(eosinophils)、靜息肥大細胞(resting mast cells)、靜息自然殺傷細胞(resting NK cells)、天然CD4 T細胞(native T CD4 cells)的豐度低于對照樣本。

    2.3 利用WGCNA識別關鍵模塊

    如圖3A所示,當軟閾值為6時(對應的R2為0.75),滿足網(wǎng)絡的無標度拓撲,且平均連通性高。通過基因相關性構(gòu)建層次聚類樹狀圖,共識別出20個相似的基因模塊(圖3B)。隨后,以免疫浸潤程度作為表型,分析模塊特征值與10種差異免疫細胞之間的相關性。包含991個基因的“MEYellow”模塊與中性粒細胞具有較強的相關性,皮爾森相關系數(shù)為0.84(圖3C)。因此認為“MEYellow”模塊為hub模塊,其中包含的基因用于后續(xù)分析。

    2.4 RSV中關鍵差異基因的篩選和功能富集分析

    將DEGs和“MEYellow”模塊基因整合,得到232個重疊的基因(圖4A)。對這些關鍵基因進行功能富集分析。GO富集分析結(jié)果顯示,基因主要影響細胞因子產(chǎn)生的正向調(diào)節(jié)(positive regulation of cytokine production)(GO-BP)、分泌顆粒腔(secretory granule lumen)(GO-CC)、模式識別受體活性(pattern recognition receptor activity)(GO-MF)(圖4B~D)。KEGG通路分析發(fā)現(xiàn)基因在中性粒細胞胞外陷阱的形成(neutrophil extracellular trap formation)等通路中顯著富集(圖4E)。

    2.5 RSV中診斷基因的鑒定

    為了進一步識別核心的診斷基因,基于lambda值進行了LASSO回歸分析。共識別出7個RSV中最具代表性的診斷基因,包括ADM、LTF、UBE2J1、PDK4、S100A4、 BST1、CAMP(圖5A~C)。隨后,基于這7個診斷生物標志物,建立了預測RSV的列線圖模型(圖5D)。校準曲線顯示疾病的預測概率和實際發(fā)生概率之間的誤差很小,表明列線圖對RSV患者的預測準確性較高(圖5E)。

    2.6 生物標志物的驗證

    利用5個外部驗證數(shù)據(jù)集對7個特征基因的表達和診斷效能進行驗證(圖6)。與對照組相比,RSV組中7個基因的mRNA水平呈現(xiàn)上調(diào)趨勢,尤其是S100A4、CAMP、ADM、BST1(均P<0.05)。在5個數(shù)據(jù)集中,診斷基因聯(lián)合模型的AUC值均大于0.9,表明模型具有較高的分類性能。此外,還對每個基因進行了ROC分析。在5個驗證集中,ADM、UBE2J1、S100A4、BST1、CAMP的AUC值均大于0.7。這些結(jié)果表明篩選的特征標志物在RSV中顯示出良好的診斷效能。

    3 討論

    有證據(jù)表明天然免疫應答的不同組成部分參與RSV感染的控制,且適當?shù)墓逃忻庖邞鹪赗SV感染的消退中具有重要的作用[8-9]。所以,RSV的免疫病理學機制的研究有助于疫苗的研發(fā)。

    本研究發(fā)現(xiàn)中性粒細胞在RSV中的浸潤水平顯著增加。已有研究報道了RSV感染與中性粒細胞之間的相關性。臨床研究發(fā)現(xiàn),RSV感染后,約80%浸潤的細胞是中性粒細胞,并且其強烈的內(nèi)流是RSV誘導的急性下呼吸道感染的氣道炎癥的常見特征[10-11]。另外,RSV感染中的中性粒細胞氣道炎癥會導致黏液生成過多、上皮損傷和氣道重塑,進而引起肺功能損傷和肺部疾病[12]。因此,在RSV感染期間靶向中性粒細胞可能是改善患兒的不良肺部預后的潛在治療方法。后續(xù)的分析中,本研究重點關注了中性粒細胞相關的關鍵基因。

    結(jié)合WGCNA和LASSO回歸分析,最終確定了7個具有診斷性能的生物標志物,包括ADM、LTF、UBE2J1、PDK4、S100A4、BST1、CAMP。ROC分析發(fā)現(xiàn)這些基因均具有良好的診斷性能。氣道ADM的異常表達可能導致上皮和黏膜功能失調(diào),與哮喘的發(fā)生和發(fā)展有關,并可能是預測社區(qū)獲得性肺炎死亡率的一個關鍵的生物標志物[13-14]。UBE2J1在內(nèi)質(zhì)網(wǎng)上可降解錯誤折疊的蛋白質(zhì),其過表達可增強登革熱病毒感染,并且該基因可通過影響泛素化過程在重癥COVID-19患者中起關鍵作用[15-16]。S100A4在炎癥過程中發(fā)揮重要作用。已有證據(jù)表明其與淋巴細胞計數(shù)呈正相關,因此該基因參與哮喘氣道炎癥的發(fā)病機制[17]。BST1可抑制中性粒細胞通過內(nèi)皮層的浸潤,發(fā)揮抗炎特性[18]。CAMP(又稱為LL-37)具有抗菌、細胞趨化和炎癥反應調(diào)節(jié)功能[19]。RSV感染后可誘導氣道上皮細胞表達LL-37,并且該基因可能是先天性宿主抵抗RSV感染的重要靶點[20]。LTF編碼的蛋白在中性粒細胞的次級顆粒中發(fā)現(xiàn),具有已知的抗病毒活性,包括RSV[21]。兒童膿毒癥心肌病患者的血清中PDK4水平顯著升高,并且與疾病的嚴重程度和器官損傷密切相關[22]。總而言之,這些研究證實篩選出來的基因均與炎癥相關的疾病有關,可間接或直接地參與RVS的發(fā)生。

    本研究還發(fā)現(xiàn)基因主要富集在“中性粒細胞胞外陷阱的形成”通路中。在RSV患兒的氣道中觀察到大量的中性粒細胞募集,而激活的中性粒細胞可釋放中性粒細胞胞外誘捕網(wǎng)(neutrophil extracellular trap,NETs)[23]。研究發(fā)現(xiàn)NETs能夠捕獲RSV,從而阻止病毒顆粒與靶細胞結(jié)合,防止感染。然而,NETs的大量產(chǎn)生會阻塞氣道并損害肺功能,從而加重嬰幼兒感染的炎癥癥狀[24]。因此,必須嚴格調(diào)控NETs的形成[25]。然而,NETs在RVS感染下的有益和有害影響之間的界限仍未揭示,所以尋找一種靶向NETs及其相關蛋白的理想療法是亟待解決的一個挑戰(zhàn)。

    我們的工作首次利用WGCNA結(jié)合LASSO方法識別了RSV的特異性診斷生物標志物,加深了對其發(fā)病機制的認識。然而,這些發(fā)現(xiàn)均是通過生物信息學分析得到的,還需要通過體內(nèi)和體外實驗來探索這些基因在RSV中的具體作用。

    參考文獻

    [1] Weinberg GA. Respiratory syncytial virus mortality among young children[J]. Lancet Glob Health, 2017, 5(10): e951-e952.

    [2] Li Y, Wang X, Blau DM, et al. Global, regional, and national disease burden estimates of acute lower respiratory infections due to respiratory syncytial virus in children younger than 5 years in 2019: a systematic analysis[J]. Lancet, 2022, 399(10340): 2047-2064.

    [3] Hall CB, Weinberg GA, Blumkin AK, et al. Respiratory syncytial virus-associated hospitalizations among children less than 24 months of age[J]. Pediatrics, 2013, 132(2): e341-e348.

    [4] Garegnani L, Styrmisdóttir L, Roson Rodriguez P, et al. Palivizumab for preventing severe respiratory syncytial virus(RSV) infection in children[J]. Cochrane Database Syst Rev, 2021, 11(11): CD013757.

    [5] 李秀東. 小兒呼吸道合胞病毒感染性肺炎防治研究進展[J]. 齊齊哈爾醫(yī)學院學報, 2015, 36(26): 3990-3991.

    [6] Dapat C, Kumaki S, Sakurai H, et al. Gene signature of children with severe respiratory syncytial virus infection[J]. Pediatr Res, 2021, 89(7): 1664-1672.

    [7] Li L, Ni YA, Song Z, et al. Identification of pathogenic genes and transcription factors in respiratory syncytial virus[J]. BMC Pediatr, 2021, 21(1): 27.

    [8] Kim HY, Yun KW, Cheong HJ, et al. Respiratory syncytial virus infection and the need for immunization in Korea[J]. Expert Rev Vaccines, 2023, 22(1): 327-340.

    [9] Luo J, Qin H, Lei L, et al. Virus-like particles containing a prefusion-stabilized F protein induce a balanced immune response and confer protection against respiratory syncytial virus infection in mice[J]. Front Immunol, 2022, 13: 1054005.

    [10] Jones AC, Anderson D, Galbraith S, et al. Personalized transcriptomics reveals heterogeneous immunophenotypes in children with viral bronchiolitis[J]. Am J Respir Crit Care Med, 2019, 199(12): 1537-1549.

    [11] 劉金玲, 陳志敏. 呼吸道合胞病毒毛細支氣管炎患兒中性粒細胞活化及蛋白酶平衡的研究[J]. 中華醫(yī)學雜志, 2009, 89(6): 397-399.

    [12] Tahamtan A, Besteman S, Samadizadeh S, et al. Neutrophils in respiratory syncytial virus infection: from harmful effects to therapeutic opportunities[J]. Br J Pharmacol, 2021, 178(3): 515-530.

    [13] Hagner S, Welz H, Kicic A, et al. Suppression of adrenomedullin contributes to vascular leakage and altered epithelial repair during asthma[J]. Allergy, 2012, 67(8): 998-1006.

    [14] Torres A, Ramirez P, Montull B, et al. Biomarkers and community-acquired pneumonia: tailoring management with biological data[J]. Semin Respir Crit Care Med, 2012, 33(3): 266-271.

    [15] Feng T, Deng L, Lu X, et al. Ubiquitin-conjugating enzyme UBE2J1 negatively modulates interferon pathway and promotes RNA virus infection[J]. Virol J, 2018, 15(1): 132.

    [16] Huang J, Li R, Cheng J, et al. Transcriptome analysis of peripheral blood mononuclear cells response in patients with severe COVID-19 reveals crucial genes regulating protein ubiquitination[J]. Med Sci Monit, 2022, 28: e937532.

    [17] Huang X, Qu D, Liang Y, et al. Elevated S100A4 in asthmatics and an allergen-induced mouse asthma model[J]. J Cell Biochem, 2019, 120(6): 9667-9676.

    [18] Wang LJ, Huang YC, Lin PY, et al. BST-1 as a serum protein biomarker involved in neutrophil infiltration in schizophrenia[J]. World J Biol Psychiatry, 2022, 23(7): 537-547.

    [19] Zhang L, Wu WKK, Gallo RL, et al. Critical role of antimicrobial peptide cathelicidin for controlling Helicobacter pylori survival and infection[J]. J Immunol, 2016, 196(4): 1799-1809.

    [20] Currie SM, Findlay EG, Mchugh BJ, et al. The human cathelicidin LL-37 has antiviral activity against respiratory syncytial virus[J]. PLoS One, 2013, 8(8): e73659.

    [21] Van Der Strate BW, Beljaars L, Molema G, et al. Antiviral activities of lactoferrin[J]. Antiviral Res, 2001, 52(3): 225-239.

    [22] Chen T, Ye L, Zhu J, et al. Inhibition of pyruvate dehydrogenase kinase 4 attenuates myocardial and mitochondrial injury in sepsis-induced cardiomyopathy[J]. J Infect Dis, 2023, 25: jiad365.

    [23] Funchal GA, Jaeger N, Czepielewski RS, et al. Respiratory syncytial virus fusion protein promotes TLR-4-dependent neutrophil extracellular trap formation by human neutrophils[J]. PLoS One, 2015, 10(4): e0124082.

    [24] Cortjens B, De Boer OJ, De Jong R, et al. Neutrophil extracellular traps cause airway obstruction during respiratory syncytial virus disease[J]. J Pathol, 2015, 238(3): 401-411.

    [25] Porto BN, Stein RT. Neutrophil extracellular traps in pulmonary diseases: too much of a good thing?[J]. Front Immunol, 2016, 7: 311.

    猜你喜歡
    中性粒細胞
    中性粒細胞/淋巴細胞比值、淋巴細胞/中性粒細胞比值與系統(tǒng)性紅斑狼瘡活動及狼瘡腎炎的關系
    不同類型突發(fā)性聾患者外周血相關細胞比值和 治療療效的關系探討
    布魯桿菌感染患者外周血中性粒細胞/淋巴細胞比值和血小板/淋巴細胞比值的臨床意義
    紅細胞分布寬度和中性粒細胞/淋巴細胞比值在系統(tǒng)性紅斑狼瘡中的應用價值
    小細胞肺癌患者化療后中性粒細胞減少與含鉑類方案化療療效的相關性研究
    術前NIR比值對非肌層浸潤性膀胱癌術后的預后評估價值
    奧曲肽對急性胰腺炎患者外周血中性粒細胞凋亡和炎癥因子的影響
    中性粒細胞和淋巴細胞與鼻咽癌療效的相關性研究
    活性氧自由基在試驗性乳腺炎大鼠發(fā)病機制中的作用
    中性粒細胞/淋巴細胞比率對小兒熱性驚厥的價值
    黄色片一级片一级黄色片| 一级黄色大片毛片| 丁香欧美五月| 久久久久久大精品| 久久性视频一级片| 国产高清videossex| 国产伦在线观看视频一区| 老熟妇仑乱视频hdxx| 色综合亚洲欧美另类图片| 国产亚洲精品久久久久5区| 99精品久久久久人妻精品| 中文在线观看免费www的网站 | 天天一区二区日本电影三级| 人人妻人人澡欧美一区二区| 久久久精品欧美日韩精品| 久久香蕉激情| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| 国产不卡一卡二| 久久久精品欧美日韩精品| 午夜精品在线福利| 性欧美人与动物交配| 熟女少妇亚洲综合色aaa.| 国产高清videossex| 亚洲av片天天在线观看| 国产精品99久久99久久久不卡| 亚洲精品国产区一区二| av中文乱码字幕在线| 国产精品免费一区二区三区在线| 黑丝袜美女国产一区| 禁无遮挡网站| 亚洲精品国产一区二区精华液| 人人妻人人看人人澡| 日韩中文字幕欧美一区二区| 国产亚洲精品第一综合不卡| 宅男免费午夜| av在线天堂中文字幕| 亚洲av熟女| 亚洲欧美一区二区三区黑人| 美女大奶头视频| 国产欧美日韩精品亚洲av| 88av欧美| 一区二区日韩欧美中文字幕| 国产av在哪里看| 日韩国内少妇激情av| 精品少妇一区二区三区视频日本电影| 伊人久久大香线蕉亚洲五| 婷婷精品国产亚洲av在线| 午夜a级毛片| 亚洲中文av在线| 夜夜爽天天搞| 麻豆成人午夜福利视频| 欧美中文综合在线视频| 国产1区2区3区精品| 女人爽到高潮嗷嗷叫在线视频| 色老头精品视频在线观看| 午夜激情av网站| 久久久久免费精品人妻一区二区 | 国产午夜福利久久久久久| 亚洲 国产 在线| 女人高潮潮喷娇喘18禁视频| 亚洲国产看品久久| 国产又色又爽无遮挡免费看| 国产私拍福利视频在线观看| 日韩三级视频一区二区三区| 免费看美女性在线毛片视频| 一边摸一边抽搐一进一小说| 亚洲人成伊人成综合网2020| 国产1区2区3区精品| 久久精品人妻少妇| 90打野战视频偷拍视频| 成人国语在线视频| 欧美成人免费av一区二区三区| 99久久99久久久精品蜜桃| 午夜福利在线在线| 99国产综合亚洲精品| 精品一区二区三区四区五区乱码| 男男h啪啪无遮挡| 亚洲精品在线观看二区| 国产一区二区三区在线臀色熟女| 国产区一区二久久| 欧美国产精品va在线观看不卡| 欧美在线黄色| 色精品久久人妻99蜜桃| 国产av一区在线观看免费| 黑人巨大精品欧美一区二区mp4| 午夜免费观看网址| 日韩大码丰满熟妇| 成人欧美大片| 国产男靠女视频免费网站| 黄色视频,在线免费观看| 国产单亲对白刺激| 老司机在亚洲福利影院| 亚洲成人国产一区在线观看| 亚洲五月婷婷丁香| 高潮久久久久久久久久久不卡| 999精品在线视频| 久久久水蜜桃国产精品网| 别揉我奶头~嗯~啊~动态视频| 一卡2卡三卡四卡精品乱码亚洲| www.熟女人妻精品国产| 国产精品久久久av美女十八| 叶爱在线成人免费视频播放| bbb黄色大片| 最新美女视频免费是黄的| 亚洲国产精品合色在线| 成人av一区二区三区在线看| 香蕉av资源在线| 国产黄色小视频在线观看| 国产真实乱freesex| 亚洲第一欧美日韩一区二区三区| 91大片在线观看| 日日夜夜操网爽| netflix在线观看网站| 欧美激情久久久久久爽电影| 久热爱精品视频在线9| 12—13女人毛片做爰片一| 亚洲真实伦在线观看| 久久亚洲精品不卡| 久久香蕉国产精品| 免费搜索国产男女视频| 三级毛片av免费| 午夜久久久在线观看| 日日夜夜操网爽| 999久久久精品免费观看国产| 香蕉av资源在线| av视频在线观看入口| 久久久精品国产亚洲av高清涩受| 高潮久久久久久久久久久不卡| 别揉我奶头~嗯~啊~动态视频| 精品国产一区二区三区四区第35| 99在线人妻在线中文字幕| 国产精品 国内视频| 亚洲人成伊人成综合网2020| 麻豆一二三区av精品| 999精品在线视频| 天堂动漫精品| 中亚洲国语对白在线视频| 国产av又大| 中出人妻视频一区二区| 国产又黄又爽又无遮挡在线| 9191精品国产免费久久| 妹子高潮喷水视频| 久久精品夜夜夜夜夜久久蜜豆 | 免费在线观看影片大全网站| 日韩 欧美 亚洲 中文字幕| 亚洲狠狠婷婷综合久久图片| 黄片大片在线免费观看| 久久精品国产99精品国产亚洲性色| 国产精品电影一区二区三区| 国产精品亚洲一级av第二区| 亚洲国产看品久久| 亚洲九九香蕉| 热99re8久久精品国产| 日韩欧美国产在线观看| 97碰自拍视频| cao死你这个sao货| 高潮久久久久久久久久久不卡| 亚洲真实伦在线观看| 免费搜索国产男女视频| 午夜免费鲁丝| 最好的美女福利视频网| 亚洲国产欧美一区二区综合| 免费在线观看视频国产中文字幕亚洲| 国产精品,欧美在线| av福利片在线| 欧美激情久久久久久爽电影| 波多野结衣av一区二区av| 91av网站免费观看| 精品久久久久久久久久久久久 | 18禁国产床啪视频网站| 岛国视频午夜一区免费看| 日韩精品免费视频一区二区三区| 美女免费视频网站| 哪里可以看免费的av片| 18美女黄网站色大片免费观看| 夜夜夜夜夜久久久久| 免费在线观看日本一区| 日本成人三级电影网站| 18禁黄网站禁片免费观看直播| 在线天堂中文资源库| 一本久久中文字幕| 亚洲avbb在线观看| 熟妇人妻久久中文字幕3abv| 大型黄色视频在线免费观看| 男女床上黄色一级片免费看| 国产精品久久视频播放| 婷婷丁香在线五月| 欧美乱码精品一区二区三区| 日本免费一区二区三区高清不卡| 免费人成视频x8x8入口观看| 黄色毛片三级朝国网站| 日本 av在线| 久久久久久久久久黄片| 女人被狂操c到高潮| 又紧又爽又黄一区二区| 久久中文字幕人妻熟女| xxx96com| 午夜影院日韩av| 久久久久九九精品影院| 中文字幕精品免费在线观看视频| 亚洲熟妇中文字幕五十中出| 亚洲第一欧美日韩一区二区三区| 不卡一级毛片| 香蕉丝袜av| 91大片在线观看| av天堂在线播放| 午夜福利在线观看吧| 看片在线看免费视频| 国内揄拍国产精品人妻在线 | 白带黄色成豆腐渣| a级毛片a级免费在线| 日本三级黄在线观看| 久久久久久九九精品二区国产 | 母亲3免费完整高清在线观看| 亚洲男人天堂网一区| 久久国产精品男人的天堂亚洲| 亚洲av美国av| 久9热在线精品视频| 精品日产1卡2卡| 亚洲一区二区三区不卡视频| 18禁观看日本| 久久久久久国产a免费观看| 两人在一起打扑克的视频| 国产精品野战在线观看| 在线观看66精品国产| 国产成人系列免费观看| 久久精品国产亚洲av高清一级| 亚洲一区高清亚洲精品| 韩国av一区二区三区四区| 久久精品人妻少妇| 欧美性猛交黑人性爽| 亚洲成人久久爱视频| 成人一区二区视频在线观看| 老汉色∧v一级毛片| 99热只有精品国产| 99riav亚洲国产免费| 国产人伦9x9x在线观看| 麻豆成人av在线观看| 久久99热这里只有精品18| 18美女黄网站色大片免费观看| 精品久久久久久成人av| 婷婷精品国产亚洲av在线| 桃红色精品国产亚洲av| 欧美不卡视频在线免费观看 | 亚洲黑人精品在线| 亚洲国产日韩欧美精品在线观看 | 一级毛片高清免费大全| 宅男免费午夜| 亚洲精品国产一区二区精华液| 91字幕亚洲| 亚洲精品久久国产高清桃花| 在线十欧美十亚洲十日本专区| 国产av不卡久久| 校园春色视频在线观看| 村上凉子中文字幕在线| 国产精品美女特级片免费视频播放器 | 两个人视频免费观看高清| 丁香六月欧美| 亚洲,欧美精品.| 91麻豆av在线| 国产三级在线视频| 麻豆成人午夜福利视频| 日本黄色视频三级网站网址| 最近最新中文字幕大全免费视频| netflix在线观看网站| 亚洲精品色激情综合| 婷婷丁香在线五月| 村上凉子中文字幕在线| 久久这里只有精品19| 成人特级黄色片久久久久久久| 日本a在线网址| 啦啦啦韩国在线观看视频| 草草在线视频免费看| 国产亚洲欧美98| 欧美亚洲日本最大视频资源| 99久久无色码亚洲精品果冻| 国产伦一二天堂av在线观看| 一个人免费在线观看的高清视频| 中文字幕久久专区| 美女高潮喷水抽搐中文字幕| 亚洲第一电影网av| 国产不卡一卡二| 国产日本99.免费观看| 999久久久精品免费观看国产| 免费人成视频x8x8入口观看| 不卡av一区二区三区| 美女高潮到喷水免费观看| 午夜久久久在线观看| 亚洲美女黄片视频| 日本撒尿小便嘘嘘汇集6| 亚洲一区二区三区色噜噜| 一级片免费观看大全| 一本综合久久免费| 天天添夜夜摸| 亚洲av片天天在线观看| 成人免费观看视频高清| 怎么达到女性高潮| 男女床上黄色一级片免费看| xxx96com| 日韩欧美三级三区| 欧美国产日韩亚洲一区| 国产主播在线观看一区二区| 精品欧美一区二区三区在线| 亚洲精品一区av在线观看| 真人一进一出gif抽搐免费| 国产乱人伦免费视频| 国产精品自产拍在线观看55亚洲| 最近最新中文字幕大全电影3 | 18美女黄网站色大片免费观看| 国产精品乱码一区二三区的特点| 久久人妻福利社区极品人妻图片| 精品国产国语对白av| 欧美激情极品国产一区二区三区| 99精品在免费线老司机午夜| 国产精品一区二区三区四区久久 | 精品国产一区二区三区四区第35| 欧美一级毛片孕妇| www日本黄色视频网| 日韩一卡2卡3卡4卡2021年| 少妇粗大呻吟视频| 午夜精品久久久久久毛片777| 女同久久另类99精品国产91| 午夜福利在线观看吧| 日韩精品中文字幕看吧| 满18在线观看网站| 色老头精品视频在线观看| 午夜亚洲福利在线播放| 不卡av一区二区三区| 俄罗斯特黄特色一大片| 成人永久免费在线观看视频| 久久中文字幕人妻熟女| 欧美乱色亚洲激情| 精品国产乱码久久久久久男人| 热re99久久国产66热| 亚洲精品粉嫩美女一区| 国内精品久久久久精免费| 在线免费观看的www视频| 久久久久精品国产欧美久久久| 欧美在线黄色| 欧美中文综合在线视频| 国产精品久久久人人做人人爽| 国产成人啪精品午夜网站| 中文字幕人成人乱码亚洲影| 欧美另类亚洲清纯唯美| 国产精品久久久久久亚洲av鲁大| 夜夜爽天天搞| 丁香欧美五月| 国产黄a三级三级三级人| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 日韩欧美三级三区| 久久九九热精品免费| 亚洲av片天天在线观看| 一区福利在线观看| 久久天躁狠狠躁夜夜2o2o| 久久 成人 亚洲| 2021天堂中文幕一二区在线观 | 成人欧美大片| 久久精品夜夜夜夜夜久久蜜豆 | 欧美日韩乱码在线| 两个人看的免费小视频| 999精品在线视频| 免费人成视频x8x8入口观看| 啦啦啦韩国在线观看视频| 国产一区在线观看成人免费| 国产欧美日韩一区二区精品| 一本大道久久a久久精品| 国产欧美日韩一区二区三| 99热6这里只有精品| 亚洲五月色婷婷综合| 久久中文字幕一级| 国产在线精品亚洲第一网站| 国产av在哪里看| 久久午夜亚洲精品久久| 一级毛片精品| 久久狼人影院| 国产成人av教育| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜精品久久久久久| 欧美另类亚洲清纯唯美| 国产视频一区二区在线看| 日韩三级视频一区二区三区| 中出人妻视频一区二区| 亚洲精品国产区一区二| 男女之事视频高清在线观看| 精品第一国产精品| 亚洲精品一区av在线观看| 一本久久中文字幕| 国产成+人综合+亚洲专区| 国产私拍福利视频在线观看| 国产精品亚洲av一区麻豆| 妹子高潮喷水视频| 搡老妇女老女人老熟妇| 国产成人av激情在线播放| 99re在线观看精品视频| 国产熟女午夜一区二区三区| 女性被躁到高潮视频| 亚洲性夜色夜夜综合| 欧美久久黑人一区二区| 日韩精品免费视频一区二区三区| 日韩一卡2卡3卡4卡2021年| av天堂在线播放| 在线av久久热| 日韩 欧美 亚洲 中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 宅男免费午夜| 亚洲精品av麻豆狂野| 国产精品九九99| 精品高清国产在线一区| 国产黄色小视频在线观看| 日韩精品中文字幕看吧| 99riav亚洲国产免费| 妹子高潮喷水视频| 久久久精品国产亚洲av高清涩受| 丁香欧美五月| 丰满人妻熟妇乱又伦精品不卡| 国产精品免费一区二区三区在线| 可以在线观看毛片的网站| 90打野战视频偷拍视频| 日本熟妇午夜| 精品久久久久久久末码| 妹子高潮喷水视频| 非洲黑人性xxxx精品又粗又长| 久久人人精品亚洲av| 国产v大片淫在线免费观看| 精品国产一区二区三区四区第35| 久久午夜综合久久蜜桃| АⅤ资源中文在线天堂| 国产成人精品久久二区二区免费| 黄频高清免费视频| 亚洲五月色婷婷综合| 久久午夜亚洲精品久久| 婷婷精品国产亚洲av在线| 久久草成人影院| 国产不卡一卡二| 精品久久久久久久久久久久久 | 日韩欧美在线二视频| 免费电影在线观看免费观看| 国产精品电影一区二区三区| 久久中文字幕一级| av免费在线观看网站| 国产色视频综合| 亚洲 国产 在线| www国产在线视频色| 婷婷六月久久综合丁香| 男女那种视频在线观看| 国产精品爽爽va在线观看网站 | 90打野战视频偷拍视频| 欧美日韩瑟瑟在线播放| 亚洲五月婷婷丁香| 日韩欧美国产一区二区入口| 国产亚洲欧美精品永久| 亚洲av美国av| 国产亚洲精品av在线| 黄色片一级片一级黄色片| 日韩欧美国产一区二区入口| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三| 免费观看精品视频网站| 听说在线观看完整版免费高清| 女性被躁到高潮视频| 俺也久久电影网| 午夜视频精品福利| 老司机午夜福利在线观看视频| 一区福利在线观看| 精品国产乱子伦一区二区三区| 国产精品日韩av在线免费观看| 国产伦一二天堂av在线观看| 中国美女看黄片| 身体一侧抽搐| 国产97色在线日韩免费| 国产主播在线观看一区二区| 免费在线观看日本一区| 欧美日韩乱码在线| 亚洲成人免费电影在线观看| 非洲黑人性xxxx精品又粗又长| 欧美在线黄色| 黄色a级毛片大全视频| a级毛片a级免费在线| 亚洲国产精品合色在线| 欧美成狂野欧美在线观看| 99国产极品粉嫩在线观看| 91麻豆av在线| 天堂动漫精品| 一本大道久久a久久精品| 一边摸一边抽搐一进一小说| 成人国产一区最新在线观看| 精品久久久久久成人av| 久久久水蜜桃国产精品网| 久久精品亚洲精品国产色婷小说| 日韩高清综合在线| 欧美又色又爽又黄视频| 日本在线视频免费播放| av有码第一页| 又黄又爽又免费观看的视频| 真人一进一出gif抽搐免费| 91成人精品电影| 欧美国产日韩亚洲一区| 久久久久久久午夜电影| 亚洲色图av天堂| 欧美人与性动交α欧美精品济南到| 国产野战对白在线观看| 国产av一区在线观看免费| 亚洲中文av在线| 国产精品国产高清国产av| 香蕉国产在线看| 村上凉子中文字幕在线| 最近最新中文字幕大全电影3 | 一级毛片女人18水好多| 一夜夜www| 亚洲 欧美 日韩 在线 免费| 国产97色在线日韩免费| 精品乱码久久久久久99久播| 99久久99久久久精品蜜桃| 日本 欧美在线| 欧美一区二区精品小视频在线| 精品欧美一区二区三区在线| 久久狼人影院| 成人18禁在线播放| 免费无遮挡裸体视频| 99re在线观看精品视频| 老司机午夜福利在线观看视频| 国产三级在线视频| 美女 人体艺术 gogo| 一边摸一边做爽爽视频免费| 丝袜人妻中文字幕| 午夜两性在线视频| 麻豆国产av国片精品| 每晚都被弄得嗷嗷叫到高潮| 97超级碰碰碰精品色视频在线观看| 国产真人三级小视频在线观看| 成人18禁在线播放| 欧美黄色片欧美黄色片| 亚洲 欧美一区二区三区| 亚洲狠狠婷婷综合久久图片| 国产亚洲精品av在线| 中文字幕另类日韩欧美亚洲嫩草| 亚洲午夜精品一区,二区,三区| 18禁裸乳无遮挡免费网站照片 | av在线播放免费不卡| 国产97色在线日韩免费| 国产不卡一卡二| 淫妇啪啪啪对白视频| 国产一区二区在线av高清观看| 搡老熟女国产l中国老女人| 国产成人精品无人区| 欧美日韩精品网址| 嫁个100分男人电影在线观看| 18禁观看日本| 免费在线观看亚洲国产| 国产又爽黄色视频| 最新美女视频免费是黄的| 女生性感内裤真人,穿戴方法视频| 欧美日本视频| 久久人妻福利社区极品人妻图片| 免费无遮挡裸体视频| 欧美色视频一区免费| 成人国产一区最新在线观看| aaaaa片日本免费| 国产一区在线观看成人免费| 级片在线观看| 亚洲熟妇熟女久久| 国产亚洲精品综合一区在线观看 | 在线观看免费日韩欧美大片| 色综合婷婷激情| 欧美成人性av电影在线观看| av免费在线观看网站| 男人的好看免费观看在线视频 | 日本黄色视频三级网站网址| 午夜精品在线福利| 亚洲av成人不卡在线观看播放网| 国产国语露脸激情在线看| 亚洲成人精品中文字幕电影| 性欧美人与动物交配| 亚洲国产精品sss在线观看| 欧美中文综合在线视频| 欧美午夜高清在线| 日韩大尺度精品在线看网址| 91成年电影在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲片人在线观看| 欧美久久黑人一区二区| 美女高潮到喷水免费观看| 成人永久免费在线观看视频| 母亲3免费完整高清在线观看| 亚洲国产高清在线一区二区三 | 中文资源天堂在线| 亚洲一码二码三码区别大吗| 国产亚洲精品久久久久久毛片| 欧美日韩福利视频一区二区| 91成人精品电影| svipshipincom国产片| 一二三四社区在线视频社区8| 波多野结衣巨乳人妻| 18禁美女被吸乳视频| 亚洲人成网站在线播放欧美日韩| 日韩有码中文字幕| 久久国产精品男人的天堂亚洲| 亚洲天堂国产精品一区在线| 视频区欧美日本亚洲| 亚洲最大成人中文| 丝袜人妻中文字幕| 波多野结衣av一区二区av| 成人欧美大片| 日韩大码丰满熟妇| 啦啦啦免费观看视频1| 人人妻人人看人人澡| 亚洲成av人片免费观看| 99re在线观看精品视频| 精品卡一卡二卡四卡免费| 午夜精品在线福利| 日韩欧美一区视频在线观看| 亚洲av成人一区二区三| 亚洲色图 男人天堂 中文字幕|