• <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.

    国产黄片美女视频| 欧美xxxx黑人xx丫x性爽| 精品无人区乱码1区二区| 十八禁国产超污无遮挡网站| 变态另类成人亚洲欧美熟女| 国产精品福利在线免费观看| 免费看美女性在线毛片视频| 在线看三级毛片| 大又大粗又爽又黄少妇毛片口| 欧美日韩综合久久久久久| 欧美日韩在线观看h| 综合色av麻豆| 熟女人妻精品中文字幕| 亚洲国产欧洲综合997久久,| videossex国产| 人妻久久中文字幕网| 欧美日韩国产亚洲二区| 久久国产乱子免费精品| 日本爱情动作片www.在线观看 | 人人妻人人澡人人爽人人夜夜 | 国产麻豆成人av免费视频| 又爽又黄a免费视频| 熟妇人妻久久中文字幕3abv| 婷婷精品国产亚洲av在线| 精品久久国产蜜桃| 久久精品夜色国产| 成人欧美大片| 欧美又色又爽又黄视频| 我要搜黄色片| 色av中文字幕| 免费无遮挡裸体视频| 美女大奶头视频| 美女xxoo啪啪120秒动态图| 国产美女午夜福利| 中文资源天堂在线| 91av网一区二区| 在线观看美女被高潮喷水网站| 韩国av在线不卡| 亚洲人成网站在线观看播放| 亚洲精品一区av在线观看| 久久久久久久午夜电影| 久久精品国产亚洲网站| 亚洲,欧美,日韩| 蜜臀久久99精品久久宅男| 国国产精品蜜臀av免费| 午夜爱爱视频在线播放| 麻豆一二三区av精品| 国产精品综合久久久久久久免费| a级毛片a级免费在线| 亚洲人与动物交配视频| 不卡视频在线观看欧美| 国产伦在线观看视频一区| 真人做人爱边吃奶动态| 精品久久久噜噜| 中国美白少妇内射xxxbb| 三级男女做爰猛烈吃奶摸视频| 我的老师免费观看完整版| 国产精品无大码| 丰满人妻一区二区三区视频av| 欧美一区二区国产精品久久精品| 日日摸夜夜添夜夜爱| 日本欧美国产在线视频| 一级av片app| 国产色婷婷99| 欧美激情久久久久久爽电影| 男女做爰动态图高潮gif福利片| 亚洲av第一区精品v没综合| 免费看av在线观看网站| 免费av观看视频| ponron亚洲| 变态另类丝袜制服| 麻豆成人午夜福利视频| 国产欧美日韩精品一区二区| 亚洲激情五月婷婷啪啪| 久久久久久久午夜电影| 特大巨黑吊av在线直播| 日韩精品中文字幕看吧| 一进一出好大好爽视频| 成人性生交大片免费视频hd| 内射极品少妇av片p| 日韩欧美精品v在线| 欧美日韩乱码在线| 亚洲不卡免费看| 欧美日韩综合久久久久久| 国产精品免费一区二区三区在线| 国产人妻一区二区三区在| 美女被艹到高潮喷水动态| 国产高清视频在线播放一区| 国产黄a三级三级三级人| 国产人妻一区二区三区在| 赤兔流量卡办理| av天堂在线播放| 国产精品女同一区二区软件| 精品午夜福利视频在线观看一区| 99久久九九国产精品国产免费| 少妇丰满av| videossex国产| 丝袜喷水一区| ponron亚洲| 狂野欧美激情性xxxx在线观看| 免费在线观看成人毛片| 免费av不卡在线播放| 91麻豆精品激情在线观看国产| 最近手机中文字幕大全| 少妇的逼水好多| 午夜福利在线在线| 69人妻影院| 亚洲av不卡在线观看| 成年版毛片免费区| 日韩制服骚丝袜av| 欧美成人一区二区免费高清观看| 亚洲第一电影网av| 99久久久亚洲精品蜜臀av| 深夜精品福利| 成人午夜高清在线视频| 欧美性感艳星| 18禁在线播放成人免费| 午夜老司机福利剧场| 日本在线视频免费播放| 一区二区三区高清视频在线| 久久久久久大精品| 欧美成人免费av一区二区三区| 黄色日韩在线| 久久久久久大精品| 高清毛片免费看| 欧美日韩一区二区视频在线观看视频在线 | 欧美区成人在线视频| 免费不卡的大黄色大毛片视频在线观看 | 欧美日本视频| 久久久成人免费电影| 在现免费观看毛片| 久久久久久久久久成人| 99久久精品热视频| 99久久精品国产国产毛片| 成年免费大片在线观看| 久久精品国产亚洲av涩爱 | 国内揄拍国产精品人妻在线| 精品人妻熟女av久视频| 网址你懂的国产日韩在线| 国产精品,欧美在线| 国产精品一区www在线观看| 看黄色毛片网站| 欧美日韩在线观看h| 天天躁日日操中文字幕| 精品久久久久久久人妻蜜臀av| 神马国产精品三级电影在线观看| av在线播放精品| 欧美国产日韩亚洲一区| 老熟妇仑乱视频hdxx| 舔av片在线| 老司机午夜福利在线观看视频| 日韩国内少妇激情av| 九九久久精品国产亚洲av麻豆| 美女cb高潮喷水在线观看| 国产欧美日韩精品一区二区| 亚洲自拍偷在线| 亚洲国产精品成人综合色| 中文在线观看免费www的网站| 1000部很黄的大片| 亚洲欧美日韩卡通动漫| 99热精品在线国产| 久久这里只有精品中国| 午夜激情欧美在线| 插逼视频在线观看| 一级毛片久久久久久久久女| 午夜激情福利司机影院| 黄色配什么色好看| 久久精品综合一区二区三区| 亚洲图色成人| 亚洲精品在线观看二区| 国产精华一区二区三区| 久久久久国产网址| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 亚洲自拍偷在线| 日韩高清综合在线| 久久久久久久久久成人| 床上黄色一级片| 国产单亲对白刺激| 国产高潮美女av| 久久久久久久久久黄片| 国产高清有码在线观看视频| 国产伦在线观看视频一区| av天堂中文字幕网| 99热只有精品国产| 一级黄片播放器| 免费无遮挡裸体视频| 久久精品影院6| 久久久国产成人精品二区| 欧美区成人在线视频| 亚洲无线在线观看| 国产在线精品亚洲第一网站| 97人妻精品一区二区三区麻豆| 亚洲va在线va天堂va国产| 午夜免费男女啪啪视频观看 | 国内久久婷婷六月综合欲色啪| 成人av一区二区三区在线看| 亚洲婷婷狠狠爱综合网| 欧美激情国产日韩精品一区| 亚洲av免费在线观看| aaaaa片日本免费| 国产黄片美女视频| 97在线视频观看| 老司机福利观看| 身体一侧抽搐| 欧美激情在线99| 成人综合一区亚洲| 看非洲黑人一级黄片| 最近在线观看免费完整版| 成人永久免费在线观看视频| 天堂√8在线中文| 亚洲欧美精品综合久久99| 国产成人精品久久久久久| 99九九线精品视频在线观看视频| 联通29元200g的流量卡| 久久精品国产自在天天线| 成人性生交大片免费视频hd| 国产视频一区二区在线看| 乱人视频在线观看| 成人综合一区亚洲| 香蕉av资源在线| a级毛色黄片| 国产熟女欧美一区二区| 一区二区三区四区激情视频 | 久久精品人妻少妇| 亚洲精品国产av成人精品 | 欧美成人一区二区免费高清观看| 久久精品国产亚洲av涩爱 | 能在线免费观看的黄片| 欧美区成人在线视频| 乱人视频在线观看| 日韩一区二区视频免费看| avwww免费| 最近的中文字幕免费完整| 天天一区二区日本电影三级| 亚洲最大成人av| 久久精品国产鲁丝片午夜精品| 欧美又色又爽又黄视频| 亚洲国产精品成人久久小说 | 搡女人真爽免费视频火全软件 | 色哟哟哟哟哟哟| 亚洲av熟女| 精品久久久久久久人妻蜜臀av| 丰满的人妻完整版| 国产69精品久久久久777片| 亚洲国产高清在线一区二区三| 免费看美女性在线毛片视频| 波多野结衣巨乳人妻| 色综合亚洲欧美另类图片| 国产亚洲欧美98| 亚洲欧美精品自产自拍| 亚洲欧美日韩卡通动漫| 亚洲久久久久久中文字幕| 久久鲁丝午夜福利片| 精品一区二区三区人妻视频| 99久久中文字幕三级久久日本| 亚洲va在线va天堂va国产| 尾随美女入室| 色综合站精品国产| 观看免费一级毛片| 久久精品夜夜夜夜夜久久蜜豆| 国产美女午夜福利| 日本-黄色视频高清免费观看| 搡老妇女老女人老熟妇| 国产精品人妻久久久久久| 国产精品一及| 午夜福利成人在线免费观看| 久久热精品热| 嫩草影院新地址| av在线老鸭窝| 久久久成人免费电影| 久久久国产成人精品二区| 欧美一区二区精品小视频在线| 久久天躁狠狠躁夜夜2o2o| 91狼人影院| 麻豆国产97在线/欧美| av黄色大香蕉| 99久久精品热视频| 99热这里只有精品一区| 18禁黄网站禁片免费观看直播| 国产精品美女特级片免费视频播放器| 日本-黄色视频高清免费观看| 免费搜索国产男女视频| 欧美激情久久久久久爽电影| 日本爱情动作片www.在线观看 | 两个人视频免费观看高清| 小蜜桃在线观看免费完整版高清| 亚洲精品粉嫩美女一区| 人妻丰满熟妇av一区二区三区| 久久亚洲精品不卡| 成年av动漫网址| 最好的美女福利视频网| 非洲黑人性xxxx精品又粗又长| 国产成人精品久久久久久| 亚洲婷婷狠狠爱综合网| 一个人看的www免费观看视频| 免费大片18禁| 尾随美女入室| 俄罗斯特黄特色一大片| 中国美白少妇内射xxxbb| 国产黄a三级三级三级人| 最近最新中文字幕大全电影3| 黄片wwwwww| 国产69精品久久久久777片| 在线观看av片永久免费下载| 精品一区二区三区视频在线| 色哟哟哟哟哟哟| 最后的刺客免费高清国语| 国产毛片a区久久久久| 乱人视频在线观看| 欧美日本视频| 深爱激情五月婷婷| 国产伦在线观看视频一区| 亚洲无线在线观看| 久久久欧美国产精品| 欧美精品国产亚洲| 精品久久久久久久久久免费视频| 国产成人91sexporn| 精品99又大又爽又粗少妇毛片| 欧美日本亚洲视频在线播放| 成年免费大片在线观看| 亚洲自拍偷在线| 少妇熟女欧美另类| 欧美又色又爽又黄视频| 麻豆乱淫一区二区| 自拍偷自拍亚洲精品老妇| 国产精品一区www在线观看| 欧美三级亚洲精品| 春色校园在线视频观看| 精品一区二区三区视频在线观看免费| 18禁裸乳无遮挡免费网站照片| 国产白丝娇喘喷水9色精品| 搡老妇女老女人老熟妇| 少妇熟女aⅴ在线视频| 欧美性感艳星| 内射极品少妇av片p| 俄罗斯特黄特色一大片| 欧美+日韩+精品| 亚洲精品国产成人久久av| 99热精品在线国产| 成年女人永久免费观看视频| 午夜福利成人在线免费观看| 日本黄大片高清| 欧美性猛交╳xxx乱大交人| 亚洲精品456在线播放app| 日韩国内少妇激情av| 国产精品av视频在线免费观看| 美女免费视频网站| 国产精品久久久久久av不卡| 99热精品在线国产| www日本黄色视频网| 国产亚洲精品久久久com| 18+在线观看网站| 国产伦精品一区二区三区视频9| 久久久国产成人免费| 国产伦在线观看视频一区| 在线免费观看的www视频| 麻豆精品久久久久久蜜桃| 一级黄色大片毛片| 床上黄色一级片| 国产 一区 欧美 日韩| 菩萨蛮人人尽说江南好唐韦庄 | 毛片一级片免费看久久久久| 网址你懂的国产日韩在线| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| 我要搜黄色片| 国产精品国产高清国产av| 女人十人毛片免费观看3o分钟| 亚洲欧美日韩东京热| 欧美日韩乱码在线| 少妇高潮的动态图| 六月丁香七月| 在线天堂最新版资源| 午夜免费激情av| 全区人妻精品视频| 老司机影院成人| 国产免费一级a男人的天堂| 插阴视频在线观看视频| av.在线天堂| 午夜福利视频1000在线观看| 国产一级毛片七仙女欲春2| 给我免费播放毛片高清在线观看| 三级国产精品欧美在线观看| 欧美xxxx性猛交bbbb| 在线观看av片永久免费下载| 97超级碰碰碰精品色视频在线观看| 久久久午夜欧美精品| a级毛色黄片| 亚洲成人精品中文字幕电影| 成熟少妇高潮喷水视频| 丰满的人妻完整版| 亚洲欧美日韩高清在线视频| 女生性感内裤真人,穿戴方法视频| 美女被艹到高潮喷水动态| 色综合站精品国产| 国产精品亚洲美女久久久| 欧美日韩乱码在线| h日本视频在线播放| 国模一区二区三区四区视频| 亚洲国产高清在线一区二区三| 国内精品宾馆在线| 久久精品人妻少妇| 欧美中文日本在线观看视频| 国内揄拍国产精品人妻在线| 久久久a久久爽久久v久久| 亚洲国产欧洲综合997久久,| 赤兔流量卡办理| 99九九线精品视频在线观看视频| 久久久久国产网址| 性色avwww在线观看| 午夜福利在线观看免费完整高清在 | 国产精品美女特级片免费视频播放器| 久久欧美精品欧美久久欧美| 极品教师在线视频| 国产伦精品一区二区三区视频9| 中文字幕人妻熟人妻熟丝袜美| 99久久中文字幕三级久久日本| 国产成人a区在线观看| 久久人人爽人人爽人人片va| 91在线观看av| 69av精品久久久久久| 三级经典国产精品| 国产亚洲av嫩草精品影院| 国内精品美女久久久久久| 别揉我奶头~嗯~啊~动态视频| 在线观看av片永久免费下载| 又爽又黄无遮挡网站| 欧美bdsm另类| 国产高清三级在线| a级毛片免费高清观看在线播放| 久久亚洲国产成人精品v| 色哟哟·www| 亚洲人成网站高清观看| 国内揄拍国产精品人妻在线| 国产视频内射| 亚洲成人久久性| 免费人成视频x8x8入口观看| 麻豆一二三区av精品| 亚洲av免费在线观看| 又黄又爽又刺激的免费视频.| 91麻豆精品激情在线观看国产| 精品一区二区三区av网在线观看| 日韩强制内射视频| 又黄又爽又免费观看的视频| 日本一二三区视频观看| 悠悠久久av| 秋霞在线观看毛片| 亚洲精品日韩av片在线观看| 亚洲国产精品国产精品| 国产成人福利小说| 内射极品少妇av片p| 亚洲婷婷狠狠爱综合网| 两个人的视频大全免费| 麻豆精品久久久久久蜜桃| 久久精品国产自在天天线| 中文字幕人妻熟人妻熟丝袜美| 少妇的逼好多水| 搡老熟女国产l中国老女人| 91久久精品国产一区二区成人| 亚洲专区国产一区二区| 麻豆国产av国片精品| av中文乱码字幕在线| 亚洲人成网站在线播| 国产精品国产三级国产av玫瑰| 天堂av国产一区二区熟女人妻| 黄色视频,在线免费观看| 99久久九九国产精品国产免费| 最近手机中文字幕大全| 亚洲av中文av极速乱| 亚洲欧美精品综合久久99| 高清日韩中文字幕在线| 真人做人爱边吃奶动态| 亚洲欧美日韩卡通动漫| 日本黄大片高清| 国产日本99.免费观看| 午夜a级毛片| 精品久久久噜噜| 国产精品精品国产色婷婷| 亚洲av五月六月丁香网| 国产淫片久久久久久久久| av在线观看视频网站免费| 九九爱精品视频在线观看| 国产男人的电影天堂91| 欧美人与善性xxx| 嫩草影院入口| 久久热精品热| 亚洲av电影不卡..在线观看| 五月玫瑰六月丁香| av在线亚洲专区| videossex国产| 在线观看av片永久免费下载| 色av中文字幕| 天堂√8在线中文| 熟妇人妻久久中文字幕3abv| 春色校园在线视频观看| 国产男靠女视频免费网站| 国产精品久久电影中文字幕| 国产黄a三级三级三级人| 一个人看视频在线观看www免费| 最后的刺客免费高清国语| 嫩草影视91久久| 成人亚洲精品av一区二区| 亚洲精品色激情综合| 91狼人影院| 一边摸一边抽搐一进一小说| 日韩人妻高清精品专区| 内地一区二区视频在线| 精品久久国产蜜桃| 色综合色国产| 熟女人妻精品中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 99热精品在线国产| av视频在线观看入口| 日本-黄色视频高清免费观看| 99热这里只有是精品在线观看| 国产精品久久久久久精品电影| 午夜视频国产福利| 六月丁香七月| 亚洲av中文av极速乱| 亚洲人成网站在线播| 亚洲熟妇熟女久久| 日韩,欧美,国产一区二区三区 | 卡戴珊不雅视频在线播放| 成人av一区二区三区在线看| 人妻久久中文字幕网| АⅤ资源中文在线天堂| 黄色配什么色好看| 午夜爱爱视频在线播放| 成人国产麻豆网| a级一级毛片免费在线观看| 免费人成视频x8x8入口观看| 国产乱人视频| 久久久久久久久久久丰满| 成人特级av手机在线观看| 婷婷精品国产亚洲av| 三级毛片av免费| 久久人人精品亚洲av| 欧美3d第一页| 精品人妻视频免费看| 色视频www国产| 日产精品乱码卡一卡2卡三| 国产美女午夜福利| 午夜福利高清视频| 最近的中文字幕免费完整| 热99在线观看视频| 日韩欧美在线乱码| 午夜老司机福利剧场| av中文乱码字幕在线| 最近最新中文字幕大全电影3| 免费看美女性在线毛片视频| 永久网站在线| 特大巨黑吊av在线直播| 成人亚洲欧美一区二区av| 国产麻豆成人av免费视频| 99久久精品热视频| 91在线观看av| 久久国内精品自在自线图片| 美女免费视频网站| 成人欧美大片| 一区二区三区免费毛片| 最好的美女福利视频网| 精品久久久久久久久久免费视频| 狠狠狠狠99中文字幕| 一区二区三区高清视频在线| 国产精品久久电影中文字幕| 欧美性猛交黑人性爽| 国产精品亚洲一级av第二区| 亚洲一区高清亚洲精品| 国产一区亚洲一区在线观看| 看免费成人av毛片| 欧美xxxx性猛交bbbb| 亚洲av.av天堂| 国产精品99久久久久久久久| 麻豆国产av国片精品| 精品国产三级普通话版| 久久久久久久久中文| 一本久久中文字幕| 美女免费视频网站| a级一级毛片免费在线观看| 中国美女看黄片| 亚洲精品一区av在线观看| 国产成人a区在线观看| 欧美激情在线99| 亚洲精品亚洲一区二区| 国产 一区 欧美 日韩| а√天堂www在线а√下载| 久久热精品热| 一本精品99久久精品77| 嫩草影院新地址| 日韩人妻高清精品专区| 国产在视频线在精品| 一级毛片久久久久久久久女| 99久久久亚洲精品蜜臀av| 国产精品永久免费网站| 国产精品女同一区二区软件| 99久久久亚洲精品蜜臀av| 国产精品永久免费网站| 观看美女的网站| 老熟妇仑乱视频hdxx| 日韩欧美国产在线观看| 亚洲七黄色美女视频| 男女啪啪激烈高潮av片| 老女人水多毛片| 搡老妇女老女人老熟妇| 男女啪啪激烈高潮av片| 国产色爽女视频免费观看| 欧美激情在线99| 蜜臀久久99精品久久宅男| 成人综合一区亚洲| 久久天躁狠狠躁夜夜2o2o| 深夜a级毛片| 亚洲av中文字字幕乱码综合| 美女cb高潮喷水在线观看| 午夜a级毛片|