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

    Identification of genes related to tubulointerstitial injury in type 2 diabetic nephropathy based on bioinformatics and machine learning

    2022-12-12 11:52:28SUJiamingPENGJingCHENHaiminZHOUYingSHIYangDONGZhaoxiWENYaxuanLINZixuanLIUHongfang
    Journal of Hainan Medical College 2022年20期

    SU Jia-ming, PENG Jing, CHEN Hai-min, ZHOU Ying, SHI Yang, DONG Zhao-xi,WEN Ya-xuan, LIN Zi-xuan, LIU Hong-fang

    1. Beijing University of Traditional Chinese Medicine, Beijing 100029, China

    2. Dongzhimen Hospital Affiliated to Beijing University of Traditional Chinese Medicine, Beijing 100700, China

    Keywords:Diabetic nephropathy Renal tubulointerstitial injury Bioinformatics Machine Learning immune infiltration Core genes

    ABSTRACT Objective: To explore the genes related to renal tubulointerstitial injury in DN and to elucidate their underlying mechanism by using bioinformatics multi-chip joint analysis and machine learning technology, so as to provide new ideas for the diagnosis and treatment of DN. Methods: Four gene expression datasets of DN tubulointerstitial tissues were retrieved from the GEO database. GSE30122, GSE47185 and GSE99340 were used as the combined microarray datasets, and GSE104954 was used as the independent verification datasets. The differentially expressed genes (DEGs) were identified by R language, and Gene Ontology(GO) enrichment, KEGG pathway enrichment, Gene Set Enrichment Analysis (GSEA) and Immune Cell Infiltration Analysis were performed. Furthermore, LASSO regression, SVMRFE and RF machine learning algorithm were used to screen core genes, while external validation and Receiver Operating Curve (ROC) analysis as well as the model of prediction nomogram were performed. Finally, the influence of the clinical characteristics of DN patients was explored by Nephroseq. Results: A total of 107 DEGs were obtained, enrichment analysis revealed that the tubulointerstitial injury in DN was mainly involved in adaptive immune response, lymphocyte mediated immunity, regulation of immune effector process and immune-inflammatory pathways such as staphylococcus aureus infection, complement and coagulation cascades, phagosomes, and Th1 and Th2 cell differentiation. In addition,cell adhesion molecule, cytokine-cytokine receptor interaction and ECM-receptor interaction pathways were also significantly enriched. Memory resting CD4 T cells, γδδ T cells, resting mast cells and neutrophil cells were up-regulated, while CD8 T cells were down-regulated.Machine learning identified MARCKSL1, CX3CR1, FSTL1, AGR2, GADD45B as core genes with good diagnostic and predictive efficacy. Conclusion: The key pathological mechanism of tubulointerstitial injury in DN is immune disorder, inflammatory reaction, cytokine action and extracellular matrix deposition. Moreover, MARCKSL1, CX3CR1, FSTL1 may be the potential biomarkers for the diagnosis and prediction of DN.

    1. Introduction

    Currently, bioinformatics methods are widely used to develop biomarkers that are relevant to diagnosis or prognosis, but given the high false positive rate due to a limited sample size, it may be difficult to obtain reliable results in single-chip data analysis[4].In this study, a bioinformatics multi-chip approach was used to search the GEO database for gene sequencing chips in renal tubulointerstitial tissues of all existing DN patients, to search for differentially expressed genes (DEGs) which are specifically related to tubulointerstitial injury of DN, and to try to combine machine learning algorithm to screen the core genes, at the same time, the specific enrichment pathway and immune infiltration mechanism of related genes were studied, and the risk factors were evaluated by correlation analysis of clinical features, therefore, the molecular mechanism of renal tubulointerstitial injury associated with DN is revealed and the potential biomarkers are explored to provide theoretical reference and scientific basis for the early diagnosis and targeted therapy of DN.

    2. Materials and methods

    2.1 Filtering of Multi-chip Dataset

    The search term "Diabetic nephropathy" was used to retrieve gene chip data from the GEO database (http://www.ncbi.nlm.nih.gov/GEO). The following were the criteria for inclusion: (1) human mRNA expression data set; (2) all samples were renal tissue; (3)used a case-control design with sample sizes exceeding 2 for each group. After screening, the public datasets of the human DN sample Gene expression profiling, GSE30122, GSE47185 and GSE99340,were obtained as the combine chip datasets, and GSE104954 as the independent verification datasets, details were shown in Table 1 .

    2.2 Data quality and variance analysis

    We converts probe ids into“Entrez ID” based on platform annotation information. After combining the matrix files of GSE30122, GSE47185 and GSE99340 into the combine chip data set, the data set was corrected by Batch normalization in the Sva package of R software (version 4.2.0) , standardized Gene expression profiling were subjected to difference analysis by the“l(fā)imma” package, and DEGs were screened using Bayesian equation multiple-test correction with |log2fc|>1 and adjusted P<0.05 as criteria, the ggplot2 package is used to draw the volcano map of DEGs, and the pheatmap package is used to draw the heat map of DEGs.

    2.3 Functional annotation and enrichment analysis of DEGs

    In order to explore the main functions and approaches of DEGs, the Gene Ontology (GO) enrichment and KEGG (Kyotoencyclopedia of genes and genomes) pathway enrichment were applied by the“Bioconductor” package in R software, biological process (BP), cellular component (CC) , molecular function (MF) and signal pathway were screened by the criteria of P<0.05. In order to get more comprehensive results, gene set enrichment analysis (GSEA)method was used to select c5.go.v7.4.symbols.gmt and c2.cp.kegg.v7.4.symbols.gmt as the reference gene set, and GSEA software was used to perform 1 000 simulation analyses of the corrected gene expression matrix of DN patients in the combine chip data set,respectively, to obtain GO and KEGG enrichment analysis results.

    2.4 Immune infiltration analysis

    Using R software to perform immune infiltration analysis on the corrected gene expression matrix of the combined microarray data set obtained above. The gene expression data of 22 immune cells were downloaded from CIBERSOFT (https://cibersort.stanford.edu/), the ratio of 22 kinds of immune cells in each sample was calculated using R software“e1071” package, and the abundance map of immune infiltration was drawn. “corrplot” and “vioplot”were used to plot the correlation heatmap and violin plot respectively to analyze the correlation and difference of immune cell infiltration distribution. P<0.05 indicated that there was significant difference between the two groups.

    2.5 Machine learning screening of core genes

    Machine learning algorithms were used to obtain more refinedmodels, which has been widely used in the search for biomarkers. In this study, the minimum absolute convergence and selection operator(LASSO), support vector machine-recursive feature cancellation(SVM-RFE) and random forest (RF) are introduced[5]. Lasso is built from the “glmnet” R package, generalized linear model for variable selection and complexity adjustment while fitting, and SVMRFE is a supervised machine learning technique that sorts features recursively[6], RF uses a decision tree classifier model to iteratively score the classification variables, thus producing high-precision classification features and jointly screening the core genes in the joint microarray data set.

    Tab 1 Dataset details

    2.6 Core gene validation and diagnostic value analysis

    We used GSE104954 as an independent validation dataset and unpaired t-test with P<0.05 as the statistical difference, to verify the difference of expression of the selected core genes between the two groups. Subsequently, a receiver operating characteristic (ROC)curve was established and the area under the ROC curve (AUC) was calculated to evaluate the effectiveness of core gene diagnosis of DN.

    2.7 Construction of prediction model

    The core gene expression matrix in combine chip data set were integrated, and the logistic regression analysis was used to build prediction model using R software to visualize as nomograms, ROC curve was used to identify the model performance in predicting tubulointerstitial injury in DN patients.

    2.8 Analysis of clinical characteristics of core genes

    Core gene clinical profiling Nephroseq database (https://www.Nephroseq.org ) is a clinical database that stores gene expression data for kidney disease and its control groups, and is widely used in the study of kidney disease[7]. In this study, the core genes selected were verified by the database, pearson correlation analysis was used to explore the effect of the expression of core gene on glomerular filtration rate (GFR) , 24h(proteinuria/24h), serum creatinine (SCr)and blood urea nitrogen (BUN) in patients with DN, with P<0.05 as the statistical criterion.

    Mrs. Thompson took pains to open it in the middle of the other presents. Some of the children started to laugh when she found a rhinestone4() bracelet5 with some of the stones missing, and a bottle that was one-quaeter full of perfume. But she stifled6 the children s laughter when she exclaimed how pretty the brcelet was, putting it on, and dabbing7 some of the perfume on her wrist. Teddy stayed after school that day just long enough to say, Mrs. Thompson, today you smelled just like my mom used to.

    3. Results

    3.1 Identification of tubulointerstitial tissue-associated DEGs in DN

    A combined chip dataset with a total sample of 71 DN patients and 87 healthy controls was obtained, with a consistent distribution of data across datasets after batch correction. The results are shown in Figure 1. A total of 107 DN tubulointerstitial tissue-associated DEGs were obtained according to screening criteria compared with healthy controls, of which 26 genes were downregulated and 81 genes were upregulated. In order to show the change and cluster relation of DEGs, the volcano map and heat map are drawn, as shown in Figure 2.

    3.2 Results of enrichment analysis

    The 107 selected DEGs targets were imported into R software for GO enrichment analysis, the biological processes are mainly involved in the regulation of endopeptidase activity, the regulation of peptidase activity, the negative regulation of hydrolase activity,the positive regulation of Cell-mediated immunity and cytokine production, etc. The cellular components include collagen-containing extracellular matrix, secretory granules and cytoplasmic vacuoles.

    The molecular functions are mainly concentrated in the extracellular matrix structure, enzyme inhibitor activity and peptidase regulatory activity. There are 34 KEGG pathway enrichment results, the main related pathways are staphylococcus aureus infection, complement and coagulation cascade, cell adhesion molecule, phagosomes,ECM-receptor interaction, and Th1 and Th2 cell differentiation, as shown in Figure 3. Furthermore, GSEA enrichment analysis showed that in the tubulointerstitial gene expression matrix of DN patients,active GO functions are enriched in immune-related processes such as adaptive immune response, lymphoid Cell-mediated immunity,and regulation of immune effector processes; The active KEGG pathway was enriched in cell adhesion molecule, cytokine-Cytokine receptor interactions, and ECM receptor interactions, as shown in Figure 4.

    Fig 1 Standardization of chip data set

    Fig 2 Volcano plots(A) and heatmap (B) of differentially expressed genes

    3.3 Results of immune infiltration analysis

    A matrix of immune cell content in tubulointerstitial tissue was constructed from 87 healthy controls and 71 patients with DN,showing differences in the distribution of infiltration of 22 immune cells across different samples, as shown in Figure 5A. By immune cell correlation analysis, the interaction between eosinophil granulocyte and naive B cells was the most obvious and positively correlated (r=0.66) , as shown in Figure 5B. Compared with the healthy controls, the infiltration of 22 kinds of immune cells in renal tubulointerstitial tissue of DN patients was significant (P<0.05) ,memory resting CD4 T cells, γδT cells, resting mast cells and neutrophil cells were up-regulated, while CD8 T cells were downregulated, as shown in Figure 5.

    3.4 The core genes selected by machine learning

    LASSO Regression, SVM-RFE algorithm and RF algorithm were used to further screen the DEGs in tubulointerstitial tissue of DN,and LASSO regression model was constructed and to cross-verified.The minimum error value corresponded to 16 characteristic genes.The SVM-RFE algorithm selected 8 feature genes by 5 fold cross validation, and RF algorithm identified 10 feature genes. Five core genes, such as MARCKSL1, CX3CR1, FSTL1, AGR2 and GADD45B,were obtained by taking intersection.

    Fig 3 GO analysis(A) and KEGG analysis(B) of differentially expressed genes

    Fig 4 GSEA enrichment analysis of tubulointerstitial gene expression matrix in DN

    Fig 5 Infiltration of 22 kinds of immune cells

    3.5 Expression of core genes and verification of diagnostic ability

    External cross-validation using the GSE104954 validation data set showed that the expression of MARCKSL1, CX3CR1 and FSTL1 genes was significantly up-regulated in the tubulointerstitium of DN patients compared with healthy controls, while the expression of GADD45B gene was significantly down-regulated, there was no significant difference in AGR2 expression between the two groups(Figure 7A-E). At the same time, the ROC curves showed that the screened 5 core genes had higher diagnostic efficacy (AUC>0.7)in distinguishing DN patients from healthy controls within the validation data set (Figure 7F-J) .

    3.6 Predictive models of nomogram

    Based on the core gene expression matrix of the joint microarray data set, the prediction model was constructed by logistic regression and visualized as nomogram. The c-index of the prediction model was 0.994, which indicating a high correlation. In addition, ROC curves indicate that the combined nomogram model has the highest performance in predicting tubulointerstitial injury in patients with DN compared with other single-core gene models, as shown in Figure 8.

    Fig 6 The core genes of DN with tubulointerstitial injury were screened by machine learning algorithm

    Fig 7 The verification results of core genes

    3.7 Analysis of clinical characteristics of core genes

    The results of core gene clinical features analysis and the correlation between core gene expression and clinical features of DN were verified by Nephroseq. Four datasets including ERCB Nephrotic Syndrome Tublnt (10 DN patients) , Ju CKD TubInt (17 DN patients) , Schmid Diabetes Tublnt (9 DN patients) and Woroniecka Diabetes Tublnt (10 DN patients) were retrieved, pearson correlation analysis showed that the expressions of MARCKSL1, CX3CR1 and AGR2 were negatively correlated with GFR (mL/min) in DN patients (Figure 9 A-C) , and the expressions of MARCKSL1, FSTL1 and GADD45B were positively correlated with urine protein (g/d) in DN patients, see Figure 9 D- F. The remaining results are not statistically significant and were not shown in the Figure.

    Fig 8 Model of prediction nomogram

    Fig 9 The clinical characteristics results of core genes

    4. Discussion

    Tubulointerstitial injury in DN is not secondary to glomerular lesions, but rather exists early in DN and has a significant role in disease progression, making the research of relevant biomarkers one of the breakthroughs in the early diagnosis of DN[8, 9]. Based on a gene sequencing chip of tubulointerstitial tissue from DN patients included in the GEO database prior to January 2022, this work investigated the underlying biological pathways utilizing bioinformatics and machine learning. A total of 107 DEGs were found to be associated with tubulointerstitial lesions in DN patients and healthy controls. The results of DEGs enrichment analysis revealed that immune dysregulation, inflammatory reaction, cytokine action, and extracellular matrix deposition were the most prominent mechanisms in tubulointerstitial injury of DN. Renal interstitial fibrosis is a leading cause of progressive kidney failure in DN and the final consequence of DN progression to end-stage renal disease[10].ECM overdeposition determines the extent and progression of renal interstitial fibrosis[11] and immune inflammatory cells that infiltrate locally early in DN release multiple cytokines, growth factors, and vasoactive substances that cause kidney-intrinsic cell damage, which activate the transcription of many ECM proteins, such as fibrinogen and fibronectin, and infiltrate the damaged site, at the same time,fibroblasts and immune cells are cell adhesion molecule to the damaged site to facilitate repair or clearance of pathogens[12,13]. But with the long-term infiltration of an immune microinflammatory state and the persistence of high glucose stimulation, repair after injury extends to pathological changes[14], with excessive deposition of ECM in the tubulointerstitium and partial hydrolysis into biologically active fragments, stimulating peripheral cells to convert to fibrotic ECM such as collagen and fibronectin, which are difficult to degrade, results in the loss of normal kidney tissue structure and function and the formation of fibrotic damage in the renal interstitium[15,16]. Based on the gene enrichment pathway in this study, it is inferred that immune dysregulation and inflammatory reaction are the driving factors of tubulointerstitial injury in DN,and the interstitial fibrosis caused by ECM deposition is the final outcome.

    Many immune system components, including memory resting CD4 T cells, γδT cells, resting mast cells, neutrophil cells, and CD8 T cells, were found to be involved in tubulointerstitial injury in DN. Although DN is not a "immune-mediated" predominant kidney disease, numerous studies have demonstrated that dysregulation of the inflammatory response caused by innate and adaptive immune abnormalities leads to progressive renal impairment in patients with DN[17], including a variety of immune and inflammatory cells (T lymphocyte, monocyte macrophage, neutrophil, etc.),inflammatory factors [ Vascular endothelial growth factor (VEGF)],tumor necrosis factor (TNF-α), Transforming growth factor-β1(TGF-β1), interleukIL (IL), c-reactive proteCRP(CRP), (NF-κB),connective tissue growth fact(CTGF), monocyte chemoattractant protein (MCP-1), etc. ]. Through a variety of signaling pathways and interaction, this exacerbates the body's micro-inflammatory state and oxidative stress, which together promote ECM deposition, apoptosis,tubular sclerosis and changes in the Hemodynamics of the kidney,resulting in continuous deterioration of renal function. Therefore,DN is considered to be an immuno-inflammatory disease in recent years[18,19].

    To realize the early diagnosis and evaluate the prognosis of DN, a new molecular feature of tubulointerstitial injury of DN,including MARCKSL1, CX3CR1, FSTL1, AGR2, and GADD45B,was identified by machine learning, the ROC curve and nomogram prediction model demonstrated the excellent diagnostic efficacy.MARCKSL1 is involved in the control of numerous physiological processes, including cell migration, secretion, proliferation, and differentiation, and is expressed in numerous organs[20].MARCKSL1 plays an important regulatory role in the immune system, it can inhibit p38, JNK MAPKs and NF-B by inhibiting their phosphorylation, and reduce the levels of inflammatory cytokines such as TNF- and IL-6, it also affects the migration and adhesion of macrophages and neutrophil[21], but its function and role in immune infiltration and micro-inflammation of renal tubules in DN require further investigation. Multiple clinical studies have shown that CX3CR1 is elevated in DN patients' kidneys[22], and Song et al.[23] have suggested that CX3CR1 plays an important role in DN by upregulating ECM synthesis, the inhibition of CX3CR1 reduces ECM deposition in a mouse model of DN and improves renal macrophage infiltration and fibrosis, thus potentially becoming an effective target for the prevention and treatment of DN. FSTL1 is a fibroblastderived cytokine, it is closely associated with fibrosis in a variety of tissues and organs (kidney, liver, lung, etc.)[24,25], and has been shown to be up-regulated in patients with chronic kidney disease disease in both clinical and basic studies, it may be a new therapeutic target for chronic kidney disease to promote the process of renal fibrosis, inflammation and apoptosis. And AGR2 may be involved in cell proliferation and growth through a variety of pathways, there is increasing recognition of its role in the development, progression and targeted therapy of cancer[26] and no significant differences have been observed in the validation data set, however, Zhou et al.[27] also suggested that AGR2 is one of the key genes for tubulointerstitial lesions in DN. GADD45B is involved in cell cycle arrest, cell survival or apoptosis, and DNA damage and repair, a few studies have reported that high glucose can stimulate the high expression of GADD45B in kidney tissue of db/db mice and human proximal tubular epithelial cells, promoting tubular epithelial-mesenchymal transition and apoptosis, however, it has also been suggested that GADD45B may have anti-apoptotic effects in other cell lines and in different disease models[28], whereas its expression was found to be significantly down-regulated in the tubulointerstitial tissues of DN patients in this study.

    On the other hand, DN is characterized by progressive impairment of renal function such as worsening glomerular filtration rate,progressive proteinuria, elevated serum creatinine and urea nitrogen levels, and tubulointerstitial injury is more sensitive to deterioration of renal function[29], according to study. In current study, the expression of MARCKSL1, CX3CR1 and AGR2 were found to be associated with GFR and urinary protein levels in DN patients, it is suggested that this method may have some significance in judging the progress of DN. These results suggest that MARCKSL1, CX3CR1 and FSTL1 have remarkable diagnostic and predictive efficacy in tubulointerstitial injury of DN, and may even be potential therapeutic targets for DN, but AGR2 and GADD45B need to be further studied.

    In conclusion, based on the integration of gene maps of differential expression of DN tubulointerstitial tissues in GEO Database and the use of bioinformatics and machine learning methods, the biological significance of related biomarkers in renal tubulointerstitial injury of DN was elucidated, and a novel concept for the diagnosis and treatment of DN was proposed; however, the related detection methods and the range and accuracy of specific applications still necessitate further investigation.

    欧美激情国产日韩精品一区| 又粗又爽又猛毛片免费看| 日韩av不卡免费在线播放| 亚洲综合色惰| 美女cb高潮喷水在线观看| 内地一区二区视频在线| 网址你懂的国产日韩在线| 国产av在哪里看| 亚洲av.av天堂| 少妇猛男粗大的猛烈进出视频 | 伦理电影大哥的女人| 久久久午夜欧美精品| 国产成人福利小说| 国产麻豆成人av免费视频| 又黄又爽又刺激的免费视频.| 精品少妇黑人巨大在线播放 | 在线免费观看的www视频| 嫩草影院入口| 最近最新中文字幕大全电影3| 色综合色国产| 国产精品久久久久久精品电影| 免费看美女性在线毛片视频| 国产精品久久久久久久久免| 日韩,欧美,国产一区二区三区 | 免费看美女性在线毛片视频| 亚洲av不卡在线观看| 国产成人a区在线观看| 成人综合一区亚洲| 一区二区三区四区激情视频 | 伦理电影大哥的女人| 亚洲经典国产精华液单| 尾随美女入室| 在线观看免费视频日本深夜| 只有这里有精品99| www.色视频.com| 床上黄色一级片| 亚洲av一区综合| 国产精品久久久久久av不卡| 免费看美女性在线毛片视频| 日日摸夜夜添夜夜爱| 在线观看午夜福利视频| 亚洲精品粉嫩美女一区| 欧美性感艳星| 午夜福利在线观看吧| 狂野欧美白嫩少妇大欣赏| 夜夜夜夜夜久久久久| 国产av一区在线观看免费| 午夜精品一区二区三区免费看| 中文资源天堂在线| 国产极品精品免费视频能看的| 国产精品久久久久久亚洲av鲁大| 国产伦精品一区二区三区视频9| 精品免费久久久久久久清纯| 欧美另类亚洲清纯唯美| 在线观看免费视频日本深夜| 久久人人爽人人爽人人片va| 非洲黑人性xxxx精品又粗又长| 国产av麻豆久久久久久久| 国产高潮美女av| 深夜精品福利| 久久午夜亚洲精品久久| 99热这里只有精品一区| 波多野结衣高清作品| 波多野结衣高清作品| 久久6这里有精品| 99热6这里只有精品| 国产色爽女视频免费观看| 联通29元200g的流量卡| 赤兔流量卡办理| 午夜精品一区二区三区免费看| 欧美最黄视频在线播放免费| 亚洲精华国产精华液的使用体验 | 在线观看一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久大av| 国产精品日韩av在线免费观看| 亚洲av第一区精品v没综合| 99久久精品国产国产毛片| 欧美色欧美亚洲另类二区| 国产大屁股一区二区在线视频| 一本精品99久久精品77| 国产成人freesex在线| 国产高清不卡午夜福利| 久久国产乱子免费精品| 黄色日韩在线| 尤物成人国产欧美一区二区三区| 亚洲不卡免费看| 午夜久久久久精精品| 在线国产一区二区在线| av福利片在线观看| 久久久久久九九精品二区国产| 午夜a级毛片| 亚洲内射少妇av| 一个人观看的视频www高清免费观看| 99在线视频只有这里精品首页| 如何舔出高潮| 联通29元200g的流量卡| 亚洲欧美日韩东京热| 国产午夜福利久久久久久| 亚洲成人中文字幕在线播放| 国产成人福利小说| 日韩欧美国产在线观看| 麻豆久久精品国产亚洲av| 精品久久国产蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 伦理电影大哥的女人| 精品一区二区免费观看| 久久九九热精品免费| 国产成人精品一,二区 | 亚洲欧美日韩高清专用| 中文字幕av在线有码专区| 国内精品美女久久久久久| 久久久久久国产a免费观看| 日韩一区二区三区影片| 观看免费一级毛片| 亚洲精品久久久久久婷婷小说 | 插阴视频在线观看视频| 日韩欧美在线乱码| 日韩制服骚丝袜av| 热99re8久久精品国产| 色5月婷婷丁香| 精品久久久久久久久亚洲| 六月丁香七月| 精品久久久久久久久久免费视频| 久久欧美精品欧美久久欧美| 亚洲国产精品成人综合色| 别揉我奶头 嗯啊视频| 亚洲欧美精品自产自拍| 国产精品爽爽va在线观看网站| 成人鲁丝片一二三区免费| 国产精品.久久久| 久久中文看片网| 国产精品美女特级片免费视频播放器| 欧美激情国产日韩精品一区| 黑人高潮一二区| 一本一本综合久久| 亚洲最大成人中文| 丰满人妻一区二区三区视频av| 国产精品一区www在线观看| 欧美日本视频| 亚洲av.av天堂| 国产v大片淫在线免费观看| 波野结衣二区三区在线| 国产精品一区二区三区四区久久| 男人狂女人下面高潮的视频| 国产亚洲精品久久久久久毛片| 免费不卡的大黄色大毛片视频在线观看 | 免费观看精品视频网站| 国产免费一级a男人的天堂| 婷婷六月久久综合丁香| 国产精品无大码| 欧美日韩精品成人综合77777| 国产私拍福利视频在线观看| 日本爱情动作片www.在线观看| 国产老妇伦熟女老妇高清| 99久久成人亚洲精品观看| 少妇熟女aⅴ在线视频| 日韩国内少妇激情av| 麻豆久久精品国产亚洲av| www.色视频.com| 99九九线精品视频在线观看视频| 麻豆成人午夜福利视频| 国产精品永久免费网站| 国产精品野战在线观看| 亚洲一区二区三区色噜噜| 99久久中文字幕三级久久日本| 综合色丁香网| 国产中年淑女户外野战色| 国产久久久一区二区三区| 欧美+日韩+精品| 99国产极品粉嫩在线观看| 有码 亚洲区| 久久久色成人| 高清毛片免费看| 日本熟妇午夜| 欧美日韩在线观看h| 亚洲国产欧洲综合997久久,| 男人的好看免费观看在线视频| 久久久久久久久久久丰满| 精品久久久久久久久久免费视频| 国产成人精品久久久久久| 亚洲av男天堂| 禁无遮挡网站| 久久这里有精品视频免费| 亚洲精品日韩在线中文字幕 | 久久久久久久亚洲中文字幕| 日本-黄色视频高清免费观看| 国产极品天堂在线| 国产亚洲精品久久久久久毛片| 在线观看午夜福利视频| 成人鲁丝片一二三区免费| 日韩成人伦理影院| av在线观看视频网站免费| 亚洲欧洲国产日韩| 美女国产视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 嘟嘟电影网在线观看| 国产亚洲精品av在线| 国产精品1区2区在线观看.| 3wmmmm亚洲av在线观看| 少妇猛男粗大的猛烈进出视频 | 久久精品人妻少妇| 看非洲黑人一级黄片| 我要搜黄色片| 国产麻豆成人av免费视频| 蜜桃亚洲精品一区二区三区| 亚洲性久久影院| 哪个播放器可以免费观看大片| 国产一区二区亚洲精品在线观看| 亚洲国产色片| 大又大粗又爽又黄少妇毛片口| 熟女电影av网| 国产一区二区三区av在线 | 欧美日本视频| 亚洲精品久久久久久婷婷小说 | 国产私拍福利视频在线观看| 亚洲精品久久国产高清桃花| av黄色大香蕉| 欧美最新免费一区二区三区| 国产高清激情床上av| 亚洲精品影视一区二区三区av| 欧美另类亚洲清纯唯美| 久久久久国产网址| 日韩精品有码人妻一区| 日韩成人av中文字幕在线观看| 欧美丝袜亚洲另类| 人人妻人人看人人澡| 亚洲av中文av极速乱| 国国产精品蜜臀av免费| 国产麻豆成人av免费视频| 日日摸夜夜添夜夜添av毛片| 深夜a级毛片| 中国美女看黄片| 亚洲在线观看片| 久久久午夜欧美精品| 最近中文字幕高清免费大全6| 国产精品蜜桃在线观看 | 赤兔流量卡办理| 永久网站在线| 亚洲激情五月婷婷啪啪| 激情 狠狠 欧美| 久久久久久久久大av| 欧美日韩精品成人综合77777| 亚洲最大成人av| 99热全是精品| 亚洲国产精品国产精品| 在线观看66精品国产| 日韩欧美精品v在线| 日本免费a在线| 国产一级毛片在线| 欧美+日韩+精品| 日韩欧美精品v在线| 搡老妇女老女人老熟妇| 如何舔出高潮| 国产精品乱码一区二三区的特点| 五月伊人婷婷丁香| 国国产精品蜜臀av免费| 久久九九热精品免费| 99精品在免费线老司机午夜| 1000部很黄的大片| 久久久精品94久久精品| 秋霞在线观看毛片| 国产毛片a区久久久久| 国产成人91sexporn| 亚洲精品国产av成人精品| 国产欧美日韩精品一区二区| 精品人妻视频免费看| 六月丁香七月| 亚洲国产精品合色在线| 欧美一区二区亚洲| 精品国内亚洲2022精品成人| 亚洲av免费在线观看| 校园春色视频在线观看| 亚洲精品日韩av片在线观看| 人妻系列 视频| 亚洲人成网站在线播放欧美日韩| 男女视频在线观看网站免费| 中文字幕av在线有码专区| 亚洲精品色激情综合| 国产91av在线免费观看| 深夜a级毛片| 三级毛片av免费| 精品熟女少妇av免费看| 插阴视频在线观看视频| 在线a可以看的网站| 在线免费观看的www视频| 校园春色视频在线观看| 婷婷色综合大香蕉| 免费看日本二区| 国产三级中文精品| 午夜激情欧美在线| 天天躁夜夜躁狠狠久久av| 美女cb高潮喷水在线观看| 白带黄色成豆腐渣| 91久久精品国产一区二区成人| 中文字幕久久专区| 国产精品人妻久久久影院| 久久久久久久久久久丰满| 少妇人妻一区二区三区视频| 99热这里只有是精品在线观看| 欧美精品一区二区大全| 国产成人影院久久av| 国产午夜精品论理片| 久久久久久久久大av| 午夜激情福利司机影院| 18禁裸乳无遮挡免费网站照片| 人妻少妇偷人精品九色| 麻豆成人av视频| 日韩精品青青久久久久久| 18禁在线播放成人免费| 女人被狂操c到高潮| 久久久久久国产a免费观看| 一级av片app| 亚洲人成网站在线播放欧美日韩| 中文精品一卡2卡3卡4更新| 久久久久网色| 麻豆成人av视频| 国产精品免费一区二区三区在线| 美女高潮的动态| 日本成人三级电影网站| 26uuu在线亚洲综合色| 国产三级在线视频| 免费搜索国产男女视频| 国产在视频线在精品| 国产女主播在线喷水免费视频网站 | 一级毛片久久久久久久久女| 国产日本99.免费观看| 国产精品.久久久| 99热只有精品国产| 夜夜夜夜夜久久久久| 成年女人看的毛片在线观看| 黄色视频,在线免费观看| 在线免费观看不下载黄p国产| 久久久久久久久久久丰满| av天堂中文字幕网| 久久久成人免费电影| 日本av手机在线免费观看| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 亚洲色图av天堂| 看非洲黑人一级黄片| 在线观看午夜福利视频| 别揉我奶头 嗯啊视频| 国产亚洲精品久久久久久毛片| 五月玫瑰六月丁香| 亚洲成a人片在线一区二区| 美女脱内裤让男人舔精品视频 | 国产69精品久久久久777片| 天天一区二区日本电影三级| 国产精品免费一区二区三区在线| 给我免费播放毛片高清在线观看| 欧美日本视频| 欧美日本亚洲视频在线播放| 性插视频无遮挡在线免费观看| 嘟嘟电影网在线观看| 免费无遮挡裸体视频| 久久久久久久久久黄片| 可以在线观看的亚洲视频| 久久亚洲国产成人精品v| 日韩中字成人| av卡一久久| 99热这里只有是精品50| 国产成人精品久久久久久| 国国产精品蜜臀av免费| 亚洲一区高清亚洲精品| av.在线天堂| 一边亲一边摸免费视频| 久久人人爽人人片av| .国产精品久久| 嫩草影院精品99| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产亚洲av涩爱 | 青春草亚洲视频在线观看| 欧美xxxx性猛交bbbb| 欧美日韩国产亚洲二区| 91aial.com中文字幕在线观看| 91久久精品国产一区二区成人| 精品99又大又爽又粗少妇毛片| 午夜精品一区二区三区免费看| 久久热精品热| 99九九线精品视频在线观看视频| 51国产日韩欧美| 国产探花极品一区二区| 亚洲一级一片aⅴ在线观看| 欧美色视频一区免费| 看黄色毛片网站| 最好的美女福利视频网| 午夜免费激情av| 免费观看精品视频网站| 日韩一本色道免费dvd| 亚洲aⅴ乱码一区二区在线播放| 老师上课跳d突然被开到最大视频| 99国产极品粉嫩在线观看| 99riav亚洲国产免费| 亚洲三级黄色毛片| 亚洲精品自拍成人| 99九九线精品视频在线观看视频| 日本在线视频免费播放| 久久热精品热| 亚洲av免费在线观看| 给我免费播放毛片高清在线观看| 精品日产1卡2卡| 伦理电影大哥的女人| 中文欧美无线码| 两个人视频免费观看高清| 99热这里只有精品一区| 日韩国内少妇激情av| 亚洲真实伦在线观看| 国产精品女同一区二区软件| 国产成人aa在线观看| 一级毛片aaaaaa免费看小| 国产精品1区2区在线观看.| 男人舔奶头视频| 哪里可以看免费的av片| 亚洲高清免费不卡视频| 欧美高清性xxxxhd video| 欧美激情国产日韩精品一区| 免费看a级黄色片| 变态另类丝袜制服| 精品欧美国产一区二区三| 直男gayav资源| 色5月婷婷丁香| 97人妻精品一区二区三区麻豆| 色综合站精品国产| 插阴视频在线观看视频| 老司机福利观看| 99久久九九国产精品国产免费| 91麻豆精品激情在线观看国产| 男人和女人高潮做爰伦理| 天堂影院成人在线观看| 三级男女做爰猛烈吃奶摸视频| 久久精品影院6| 在现免费观看毛片| 欧美潮喷喷水| av国产免费在线观看| 亚洲va在线va天堂va国产| 国内久久婷婷六月综合欲色啪| 波野结衣二区三区在线| 亚洲av男天堂| 国产片特级美女逼逼视频| 中文字幕久久专区| 国产精品久久久久久亚洲av鲁大| 最后的刺客免费高清国语| 十八禁国产超污无遮挡网站| 99视频精品全部免费 在线| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品国产亚洲av香蕉五月| 国产精品久久久久久亚洲av鲁大| 99久久无色码亚洲精品果冻| 午夜免费男女啪啪视频观看| 男女下面进入的视频免费午夜| 可以在线观看的亚洲视频| 欧美激情久久久久久爽电影| 日韩制服骚丝袜av| 一级黄片播放器| 看免费成人av毛片| 亚州av有码| 麻豆精品久久久久久蜜桃| 国产精品国产高清国产av| 热99在线观看视频| 国产精品久久久久久精品电影小说 | 中文字幕久久专区| 波多野结衣高清无吗| 两个人视频免费观看高清| 国产精品永久免费网站| 精品一区二区三区人妻视频| 国产精品爽爽va在线观看网站| 免费看av在线观看网站| 婷婷精品国产亚洲av| 久久人人爽人人爽人人片va| 久久久久性生活片| 黄色一级大片看看| 亚洲精品粉嫩美女一区| 国内久久婷婷六月综合欲色啪| 久久久久久九九精品二区国产| 日日干狠狠操夜夜爽| 18+在线观看网站| 国产高清激情床上av| 热99re8久久精品国产| 成人亚洲欧美一区二区av| 免费看美女性在线毛片视频| 小蜜桃在线观看免费完整版高清| 日韩制服骚丝袜av| 人妻系列 视频| 成年女人看的毛片在线观看| 舔av片在线| 色尼玛亚洲综合影院| 边亲边吃奶的免费视频| 深爱激情五月婷婷| 国产成人a区在线观看| 少妇人妻一区二区三区视频| 亚洲高清免费不卡视频| 少妇人妻一区二区三区视频| 91久久精品国产一区二区三区| 欧美成人免费av一区二区三区| 神马国产精品三级电影在线观看| 在线观看午夜福利视频| 黄色一级大片看看| 亚洲欧洲日产国产| 人妻久久中文字幕网| 久久99精品国语久久久| 成人av在线播放网站| 美女脱内裤让男人舔精品视频 | 简卡轻食公司| 三级经典国产精品| a级毛片a级免费在线| 中文在线观看免费www的网站| 国产成人精品久久久久久| 在线观看午夜福利视频| av天堂中文字幕网| 一个人观看的视频www高清免费观看| 成人高潮视频无遮挡免费网站| 亚洲av第一区精品v没综合| 国产伦在线观看视频一区| 男女边吃奶边做爰视频| 久久久久久久久中文| 中文字幕熟女人妻在线| 欧美日韩国产亚洲二区| av视频在线观看入口| 亚洲经典国产精华液单| 国产精品久久久久久av不卡| 一级黄色大片毛片| 99riav亚洲国产免费| 国产精品1区2区在线观看.| 久久精品国产亚洲av香蕉五月| 久久久久久久久大av| 久久久国产成人精品二区| 在线观看av片永久免费下载| 久久午夜福利片| 国产伦在线观看视频一区| 97在线视频观看| 日韩一区二区视频免费看| 国产亚洲91精品色在线| 高清在线视频一区二区三区 | 亚洲欧美精品综合久久99| 亚洲激情五月婷婷啪啪| 亚洲精品国产av成人精品| 大香蕉久久网| 日本免费a在线| 老熟妇乱子伦视频在线观看| 亚洲五月天丁香| 久久精品国产亚洲av涩爱 | 国产精品永久免费网站| 国产精品无大码| 亚洲四区av| 三级经典国产精品| 狠狠狠狠99中文字幕| 国产色婷婷99| 国产伦精品一区二区三区视频9| 久久精品久久久久久久性| 国产精品一区二区性色av| 免费观看的影片在线观看| 亚洲国产高清在线一区二区三| 国产成人a∨麻豆精品| 中文字幕av在线有码专区| 亚洲精品色激情综合| 欧美另类亚洲清纯唯美| 神马国产精品三级电影在线观看| 午夜福利视频1000在线观看| 婷婷精品国产亚洲av| 校园人妻丝袜中文字幕| 看片在线看免费视频| 国产精品一区二区三区四区久久| 国产av一区在线观看免费| 美女国产视频在线观看| 男女那种视频在线观看| 国产伦精品一区二区三区视频9| 极品教师在线视频| 久久久午夜欧美精品| 在线免费十八禁| 国产成人精品久久久久久| 能在线免费看毛片的网站| 免费无遮挡裸体视频| 日本爱情动作片www.在线观看| 久久久久久久久中文| 欧美精品国产亚洲| 中文字幕人妻熟人妻熟丝袜美| av免费观看日本| 99热这里只有是精品在线观看| 三级男女做爰猛烈吃奶摸视频| 色哟哟·www| 人人妻人人看人人澡| 黑人高潮一二区| 色视频www国产| 国产女主播在线喷水免费视频网站 | 欧美xxxx黑人xx丫x性爽| 少妇熟女aⅴ在线视频| 在线观看美女被高潮喷水网站| 中文字幕av在线有码专区| a级毛片a级免费在线| 欧美人与善性xxx| av在线天堂中文字幕| 国产伦精品一区二区三区四那| 亚洲自偷自拍三级| 男女那种视频在线观看| 丰满的人妻完整版| 亚洲不卡免费看| 国内揄拍国产精品人妻在线| 女的被弄到高潮叫床怎么办| 国产精品无大码| 欧美3d第一页| 亚洲av中文av极速乱| 日韩一区二区三区影片| 欧美日韩精品成人综合77777| 女的被弄到高潮叫床怎么办| 色噜噜av男人的天堂激情| 99热这里只有精品一区| 男人舔女人下体高潮全视频| 欧美丝袜亚洲另类| 欧美+亚洲+日韩+国产| 麻豆乱淫一区二区| 亚洲精品成人久久久久久| 国产v大片淫在线免费观看| 亚洲自拍偷在线| 亚洲精品456在线播放app|