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

    Eight key long non-coding RNAs predict hepatitis virus positive hepatocellular carcinoma as prognostic targets

    2019-12-14 09:37:42ZiLinHuangWangLiQiFengChenPeiHongWuLuJunShen

    Zi-Lin Huang,Wang Li,Qi-Feng Chen,Pei-Hong Wu,Lu-Jun Shen

    Zi-Lin Huang,Wang Li,Qi-Feng Chen,Pei-Hong Wu,Lu-Jun Shen,Department of Medical Imaging and Interventional Radiology,Sun Yat-sen University Cancer Center,Guangzhou 510060,Guangdong Province,China

    Zi-Lin Huang,Wang Li,Qi-Feng Chen,Pei-Hong Wu,Lu-Jun Shen,State Key Laboratory of Oncology in South China,Collaborative Innovation Center for Cancer Medicine,Guangzhou 510060,Guangdong Province,China

    Abstract

    Key words: Long non-coding RNAs; Hepatitis virus; Hepatocellular carcinoma;Prognostic signature; Least absolute shrinkage and selection operator

    INTRODUCTION

    Hepatocellular carcinoma (HCC),one of the most frequent liver cancer types,is a severe worldwide health problem[1].Chronic hepatitis B virus (HBV) infection and/or hepatitis C virus (HCV) infection have been identified as the leading causes of HCC[2].Specifically,the prevalence of HCV infection is predominant in European countries,which accounts for about 60% of all HCC cases in Italy and Spain; by contrast,HBV infection is more commonly seen in eastern Asian countries such as China,with an incidence rate of over 2/3[3].Specifically,hepatitis viruses have been epidemiologically associated with the development of HCC[4,5].

    It has been extensively determined that hepatitis viruses are associated with HCC etiologically; nonetheless,it is still not completely clear about the molecular mechanism of hepatitis viruses in inducing liver cancer.An increasing number of mutation genes,like TP53,CTNNB1,and mTOR[6],are found to participate in HCC genesis as well as progression.Nonetheless,HCC is a highly heterogeneous disease,which can determine the altered biological progression of HCC and add to the difficulty in predicting prognosis.As a result,it is urgently needed to search for the efficient prognostic biomarkers for HCC.

    Long non-coding RNAs (LncRNAs) have been identified to be the superior prognostic biomarkers over other cancer hallmarks,which can be ascribed to their unique advantages,as shown below.First,lncRNA expression is greatly distinct among different tissues,in various diseases,as well as at different disease progression stages; as a result,they can better represent the disease characteristics[7].Second,lncRNAs have been found to regulate the expression of genes at the transcriptional,post-transcriptional,and epigenetic levels[8,9]; on this account,their levels as well as functions show a close association with cancer progression.As a result,increasing studies have been carried out to illustrate the clinical value of lncRNAs among cancers,such as HCC[10].

    In this study,hepatitis virus positive HCC (VHCC) patients were appropriately selected to explore the difference in the lncRNA expression profiles,so as to identify the candidate lncRNA biomarkers based on The Cancer Genome Atlas (TCGA)database.Moreover,the key lncRNAs were determined using the least absolute shrinkage and selection operator (LASSO) algorithm; thereafter,the risk score system for VHCC was established.This study aimed to establish a potent lncRNA signature based on the expression profiles through thorough genomic data analysis,thus improving the accuracy in predicting the prognosis of VHCC.

    MATERIALS AND METHODS

    Patient datasets as well as processing

    The data were downloaded from TCGA database.Gene expression,together with the clinical information of HCC patients (https://tcga-data.nci.nih.gov/),was downloaded using the Data Transfer Tool (provided by GDC Apps).Among them,151 patients were virus-related.Patients with unknown lncRNA expression were eliminated from this study (n= 2); as a result,a total of 149 HCC cases were finally enrolled for analysis.The flowchart of analysis is shown in Figure1.A total of 149 virus-induced HCC cases were enrolled for comprehensive integrated analysis.The data were publicly available,so approval was exempted from the Ethics Committee.Data were processed in accordance with the data access policies as well as the NIH TCGA human subject protection policies (http://cancergenome.nih.gov/publications/publicationguidelines).

    LncRNA expression in HCC patients was derived from Illumina HiSeq RNASeq platform,which was then standardized based on TCGA.Then,the R software“edgeR” package was employed to detect the differentially expressed RNAs,so as to identify the differentially expressed lncRNAs (DElncRNAs),and the thresholds were set at log2fold change > 2.0 as well as adjustedPvalue < 0.05.

    Construction of lncRNA signature

    The relationship of each lncRNA expression with the overall survival (OS) of patients was calculated using the univariate Cox model.Typically,lncRNAs would be deemed as statistically significant atP< 0.05 upon univariate Cox analysis.Subsequently,the chosen lncRNAs were screened and verified through LASSO regression with the R software “glmnet” package.Finally,the lncRNA-based prognostic risk score was constructed by linearly combining the products of expression level with the regression model (β) according to the formula below:Risk_score = βlncRNA1 ×lncRNA1 expression + βlncRNA2 × lncRNA2 expression + · ···· + βlncRNAn ×lncRNAn expression.

    Confirmation of lncRNA signature

    Patients together with their survival information were distributed according to the risk_score.Furthermore,patients were classified into a high-risk group and low-risk group according to the median risk_score.Then,the Kaplan-Meier survival curves were plotted,which could predict the high or low risk of patients.Thereafter,the univariate Cox proportional hazards regression analysis was performed.Then,the sensitivity and specificity regrading survival prediction were compared by the risk_score,and the receiver operating characteristic (ROC) curves with time dependence was then employed to evaluate the accuracy in predicting the 5-year prognosis.Moreover,multivariate Cox regression was also performed to examine whether lncRNA risk_score was predicted independently from other clinical factors.The R-package “party” “ctree” function was used to construct a conditional inference tree for better illustration.Besides,a nomogram was also established with the R project “rms” package.In addition,the predictive performance of the nomogram was verified for discrimination and calibration.Subsequently,subgroup analysis based on virus type was carried out.All the two-sidedP< 0.05 would be deemed to be statistically significant.The R software (version 3.4.1,R Foundation) was employed for all analyses.

    Figure1 Flowchart of study design.HCC:Hepatocellular carcinoma; TCGA:The Cancer Genome Atlas; lncRNA:Long non-coding RNAs.

    Functional analysis of lncRNAs for prognosis prediction

    Pearson correlation analysis was carried out for the screened lncRNAs as well as the protein coding genes in the TCGA database according to their different expression quantities.Typically,a coefficient of correlation > 0.4 andP< 0.001 would be deemed to show a significant correlation.Enrichment analysis was performed with the DAVID Bioinformatics Tool (version 6.8,https://david-d.ncifcrf.gov/).

    RESULTS

    DElncRNA identification

    The significant DElncRNAs were identified by comparing the tumor samples with non-tumor samples.Altogether,1420 DElncRNAs (including 1358 up-regulated and 62 down-regulated ones) were recognized using the R software “edgeR” package.Afterwards,the heat map (Figure2A) and volcano of top 20 DElncRNAs (Figure2B)were built.

    Construction of an lncRNA signature

    Univariate Cox regression analysis was performed between DElncRNAs and patients’survival; as suggested by the results,altogether 406 lncRNAs were significantly associated with OS (P< 0.05).Next,LASSO regression was employed for verifying further variables (Figure3).Eight lncRNAs were generated in this process,including AC005722.2,AC107959.3,AL353803.1,AL589182.1,AP000844.2,AP002478.1,FLJ36000,and NPSR1-AS1.Then,the prognostic risk_score was imputed as follows:(0.15103*AC005722.2 expression)+(0.07857*AC107959.3 expression)+(0.09283*AL353803.1 expression)+(0.17744*AL589182.1 expression) + (0.07857*AP000844.2 expression)+(0.09888*AP002478.1 expression)+(0.13460*FLJ36000 expression)+(0.06309*NPSR1-AS1 expression).

    Confirmation of the lncRNA signature

    Figure2 Heat map and volcano plot.A:Heat map demonstrating the differential expression of long non-coding RNAs (lncRNAs).The X-axis represents the samples,while the Y-axis denotes the differentially expressed lncRNAs.Green color is indicative of the down-regulated genes,while the red color stands for the upregulated genes; B:Volcano plot showing the differential expression of lncRNAs.The X-axis indicates the false discovery rate values after log transformation,whereas the Y-axis suggests the log2 (fold change) values for lncRNA expression.FDR:False discovery rate; FC:Fold change.

    The risk score of every case was calculated; then,patients were classified into a highrisk group or low-risk group based on the median threshold (Figure4A and B).The Kaplan-Meier curves for the high-risk as well as low-risk groups are presented in Figure4C.Patients in the high-risk group were associated with a shorter OS compared with those in the low-risk group (P= 2e-10).Meanwhile,the risk score hazard ratio (HR) upon univariate Cox proportional hazards regression was 5.360[95% confidence interval (CI):3.012-9.538] (Table1).Besides,multivariate Cox proportional hazards regression with adjusted clinical covariate also came to the same results (HR=1.94,95%CI:1.61-2.34) (Table2).Subsequently,the ROC curves with time dependency were adopted to assess the prognostic value of the eight-lncRNA signature.In addition,the areas under curves (AUCs) of the prognostic model constructed based on the eight-lncRNA signature were 0.73,0.758,and 0.788 for 1-,3-,and 5-year OS,respectively (Figure4D).

    The first node of the conditional inference tree divided the data according to the risk_score,and the lower risk_score (≤0.855) was related to a better OS (P< 0.001)(Figure5).For the lower risk_score (≤0.855),the branch was further divided depending on the stage (P= 0.006),and it was obvious that the risk_score ≤ 0.855 in stage I or II patients was associated with the best survival,while patients with the risk_score > 0.855 had the worst survival.

    A nomogram for prognosis was constructed based on the above-mentioned independent factors for prognosis for the sake of visualization as well as prediction,so as to construct an integrated predicting factor for HCC.According to the multivariate Cox regression analysis,the eight-lncRNA risk_score,TNM stage,and race were recognized as independent risk factors (P< 0.05) (Table2).Specifically,the nomogram was constituted by three factors (Figure6A),with a C-index of 0.763(95%CI:0.700-0.826).Moreover,the calibrated curves for the risk of 1-,3-,and 5-year recurrence were highly consistent between the nomogram predicted results and the observed results (Figure6B-D).

    Subgroup analysis by etiology

    Figure3 Plots of regression coefficient based on least absolute shrinkage and selection operator regression.A:Least absolute shrinkage and selection operator (LASSO) coefficient profiles of the fractions of the 406 significant long non-coding RNAs (lncRNAs) in univariate Cox regression analysis.The solid vertical line implies the optimal lambda value determined by cross-validation; B:Cross-validation for tuning parameter selection in the LASSO model.The vertical lines are drawn at the optimal values by minimum criteria and 1-standard error criteria.And the solid vertical line represents that eight lncRNAs were finally determined.

    Subsequently,the patients were further classified into three subgroups,namely,HBV+ HCV (n= 82),HBV alone (n= 58),and HCV alone (n= 9) groups.Kaplan-Meier curves for the high-risk and low-risk groups are displayed in Figure7.The OS of patients benefited from the low-risk score in the HBV + HCV group (P= 6e-5,Figure7A) and HBV alone group (P= 5e-6,Figure7B).At the same time,the AUCs regarding the 1-,3-,and 5-year OS in the HBV + HCV model were 0.681,0.759,and 0.724,respectively (Figure7D),while those in the HBV model were 0.787,0.720,and 0.826,respectively (Figure7E).Curves and correspondingPvalues for HCV patients are also shown in Figure7C and F.

    Functional analysis of prognostic lncRNAs

    LncRNA-related mRNAs were first filtered out by Pearson correlation analysis and then applied into online tool DAVID for pathway enrichment,so as to examine the possible biological functions of the above-mentioned identified lncRNAs.The top 10 pathways with the highest enrichment are presented in Figure8.Among those,the Wnt signaling pathway,angiogenesis,the p53 pathway,and the PI3 kinase pathway were recognized as the canonical pathways in the initiation as well as progression of HCC.

    DISCUSSION

    Few HCC prognostic models have been constructed on the basis of lncRNAs up to now.To our knowledge,little is known about lncRNA signature for hepatitis virus infection associated HCC patients.Meanwhile,the majority of studies analyzed virus combined with non-virus related HCC at the same time[11,12].Different from previous studies that carried out analysis regardless of etiology,we conducted the current research in virus positive patients only.The DElncRNAs with marked correlation with OS were comprehensively screened by applying biostatistical methods and univariate Cox analysis.Then,the LASSO regression was applied to analyze the lncRNA data from the carefully chosen patients from TCGA,and eight lncRNAs were identified,which were recognized to be the independent factors for the prognosis of VHCC.Afterwards,the identified significant lncRNAs were employed for constructing a prognostic model,and its prognostic value was verified through Kaplan-Meier analysis,time-dependent ROC curve,and Cox regression analysis.A nomogram was developed to present the constructed risk_score system,and it displayed a favorable accuracy in OS prediction for VHCC.Further functionalanalysis uncovered that several pathways play crucial roles in the effects of the eight lncRNAs on the progression of HCC.

    Table1 Stratified overall survival analysis

    Our research displayed obvious advantages when compared with other research on the prognostic value of lncRNAs in HCC.First,patient homogeneity was markedly enhanced through careful patient selection before analysis,since the ignorance of patient features might result in false relationships,thus invaliding the conclusions.Existing treatments can treat chronic viral hepatitis,and some research has been performed previously to domenstrate the influence of antiviral treatment on disease progression as well as HCC development[13].Moreover,many studies suggest that the survival of virus-related HCC patients can be prolonged through curative hepatectomy and/or local tumor ablation combined with adjuvant antiviral therapy[14-16].On the other hand,virus and non-virus HCCs have been recognized to have apparently different outcomes,and different therapeutic strategies were recommended for them.Therefore,the prognostic models should also be developed in a distinct way.All the 149 patients were virus positive in our study.Second,with regard to the methodology,the use of LASSO penalized regression could boost the accuracy of bioinformatic analysis.Different from the conventional stepwise regression applied in prior research,the LASSO algorithm contributed to analyzing all independent factors at the same time,which tended to screen the most significant variables[17].Then,the coefficients of the less significant variables would become zero following introduction into the penalty after regularization[18].Consequently,such formulation approach would display a higher accuracy than stepwise regression using the multivariate Cox model,especially in the presence of huge datasets,such as genomics[19].Third,an easy-to-use nomogram model was constructed on the basis of Cox analysis among the carefully screened patients with clinicopathological features,which would facilitate the prediction for individual patients.

    As hepatitis viruses vary in different countries[3],our cohort was also divided into three subgroups based on the virus type (namely,HBV alone,HCV alone,and combined).Subgroup survival analysis showed that the as-constructed eight-lncRNA signature performed well in patients with combined virus infection (5-year AUC =0.724),especially for those with HBV infection (5-year AUC = 0.826).No obvious OS heterogeneity was detected among the HCV patients (P= 0.51),which could be resulted from the small sample size (n= 9).

    No study had investigated the biological functions for the eight lncRNAs identified in HCC; however,they could affect the genesis as well as the progression of HCC by the Wnt signaling pathway[20-24],angiogenesis[25,26],p53 pathway[27,28],and PI3 kinase pathways[29-31],as suggested by the pathway enrichment results.More importantly,this research could provide new evidence that lncRNAs,whose biological functionshad not been reported before,might potentially serve as the potent markers to predict the prognosis of HCC.However,such findings should be further validated in future studies,and the molecular characteristics should also be investigated.

    Table2 Cox regression analysis for the overall survival-related clinicopathological features

    To our knowledge,this is the first study that constructs an lncRNA nomogram to predict the prognosis of VHCC.Nonetheless,some limitations should be noted in this study.First,no experimental studyin vitroorin vivowas performed to verify the effect of the eight-lncRNA signature,which was constructed based on online datasets by the bioinformatic approach,in predicting the prognosis of HCCs.Second,data regarding more valuable clinical features are not available,including viral hepatitis duration,state of virus,Milan criteria,HCC treatment,and anti-viral treatment.Moreover,our results should be further validated in future studies.

    To sum up,a novel lncRNA signature has been uncovered in this study based on comprehensive bioinformatic analysis combined with the expression patterns as well as the clinical information from the carefully selected cohort.Specifically,this risk model can be used as a biomarker to independently predict the prognosis of VHCC.Nonetheless,our results should be validated in studies that investigate the precise mechanisms regarding VHCC progression,as well as those that explore the functions of the eight identified lncRNAs.

    Figure4 Establishment and verification of a long non-coding RNA signature for predicting the prognosis of hepatitis virus positive hepatocellular carcinoma using the cut-off value 0.798.A:Distribution of the long non-coding RNA (lncRNA) risk scores; B:Survival status as well as overall survival (OS); C:Kaplan-Meier curve regarding the OS of the high-risk as well as low-risk group classified based on median risk score; D:Receiver operating characteristic curves regarding the 1-,3-,and 5-year survival discriminated by the lncRNA signature.AUC:Area under the curve.

    Figure5 Conditional interference tree for prediction.The nodes of the conditional inference tree were calculated automatically based on the influence power of variables.The first node of the tree was risk_score,while the second was stage.Node 3,4,or 5 was the corresponding group a certain patient would be assigned to.For example,a stage I patient with a risk_score ≤ 0.855 would be divided into group node 3,which was related to the best overall survival.

    Figure6 Overall survival nomogram predicting 1-,3-,and 5-year survival of hepatocellular carcinoma,and calibrated curves to predict the 1-year,3-years,and 5-year patient survival.A:Overall survival nomogram predicting 1-,3-,and 5-year survival of hepatocellular carcinoma; B:Calibrated curve to predict the 1-year patient survival; C:Calibrated curve to predict the 3-years patient survival; D:Calibrated curve to predict the 5-years patient survival.In calibrated curves,the X-axis represents the survival prediction by nomogram,and the Y-axis represents the actual observation which showed a good agreement with prediction.

    Figure7 Subgroup analysis regarding the etiology.A-C:Survival curves of risk_score for hepatitis B virus + hepatitis C virus (HBV + HCV) (A),HBV alone (B),and HCV alone (C) infected patients with HCC; D-F:Receiver operating characteristic curves for HBV + HCV (D),HBV alone (E),and HCV alone (F) infected patients with HCC.HBV:Hepatitis B virus; HCV:Hepatitis C virus.

    Figure8 Functional annotations for the eight long non-coding RNAs.The 10 most remarkable pathways enriched based on the co-expressed mRNAs of eight long non-coding RNAs (lncRNAs).Corresponding lncRNAs involved in the deregulated pathways are also presented.A:AC005722.2; B:AC107959.3; C:AL353803.1;D:AL589182.1; E:AP000844.2; F:AP002478.1; G:FLJ36000; H:NPSR1-AS1.

    ARTICLE HIGHLIGHTS

    Research background

    Hepatocellular carcinoma (HCC),one of the most frequent liver cancer subtype,has posed a serious health issue in the world.Hepatitis virus is recognized to be a major factor leading to HCC.Recently,a verity of genetic markers as well as prediction models are proposed to improve HCC treatment.At the same time,numerous statistical techniques are also utilized to mine data in numerous large-scale public databases.Thanks to the well-developed clinical approaches,prognosis models with a higher accuracy and robustness can be established for HCC.

    Research motivation

    This study aimed to establish a prognosis model on the basis of HCC molecular biomarkers.Typically,long non-coding RNAs (lncRNAs) have been recognized as new predictive factors.Great efforts have been made to establish lncRNA-based HCC models,however,the lncRNA features for hepatitis virus positive HCC (VHCC) are not available at present.

    Research objectives

    This study aimed to develop a prognostic lncRNA feature for VHCC based on candidate lncRNAs through analyzing data collected from The Cancer Genome Atlas (TCGA) database.

    Research methods

    Specifically,the least absolute shrinkage and selection operator (LASSO),the most advanced statistic algorithm,was used to establish the prediction model.This approach was carried out on the basis of typical lncRNAs selected according to the expression profiles of lncRNAs collected from the TCGA database.The as-established lncRNA feature was validated,and its clinical applicability was also examined.

    Research results

    Using LASSO,a risk score system was established to predict the overall survival (OS) for VHCC,which incorporated eight lncRNAs (including AC005722.2,AC107959.3,AL353803.1,AL589182.1,AP000844.2,AP002478.1,FLJ36000,and NPSR1-AS1).Notably,the as-constructed lncRNA feature was helpful in stratifying the risk for VHCC.To extend the model application range,a nomogram was also constructed,which involved both the lncRNA feature and other clinical features.Our results also suggested that the lncRNAs were markedly enriched in the Wnt signaling pathway,angiogenesis,the p53 pathway,and the PI3 kinase pathway.

    Research conclusions

    The as-constructed signature based on eight lncRNAs displayed a favorable capacity in predicting the prognosis for VHCC patients,which can also contribute to risk stratification and may provide more useful clinical advice for individual patients.

    Research perspectives

    Our results further support that lncRNAs can serve as possible functional regulating factors for the progression of VHCC.It is the future direction to search for the effective molecular biomarkers and predictive factors for the prognosis of HCC patients.

    国产成人影院久久av| 欧美xxⅹ黑人| 少妇粗大呻吟视频| 99精品久久久久人妻精品| 亚洲精品中文字幕一二三四区 | 亚洲精品中文字幕在线视频| 天堂8中文在线网| 久久久国产成人免费| 老熟妇仑乱视频hdxx| 成人国产av品久久久| 99久久人妻综合| 日本五十路高清| 日韩 欧美 亚洲 中文字幕| 91麻豆精品激情在线观看国产 | 午夜福利在线观看吧| 最近最新中文字幕大全免费视频| 久久女婷五月综合色啪小说| 丝袜喷水一区| 建设人人有责人人尽责人人享有的| 久久久精品国产亚洲av高清涩受| 精品亚洲乱码少妇综合久久| 精品国产超薄肉色丝袜足j| 男人爽女人下面视频在线观看| 一级黄色大片毛片| 国产深夜福利视频在线观看| 人妻 亚洲 视频| 国产精品国产三级国产专区5o| 99国产精品一区二区三区| 亚洲综合色网址| 在线天堂中文资源库| 欧美黄色片欧美黄色片| 老熟妇仑乱视频hdxx| 国产视频一区二区在线看| 中文字幕高清在线视频| 亚洲五月色婷婷综合| 亚洲熟女毛片儿| videos熟女内射| 91字幕亚洲| 无限看片的www在线观看| 亚洲国产欧美网| 亚洲国产av影院在线观看| 一级a爱视频在线免费观看| 人人妻人人添人人爽欧美一区卜| 亚洲第一欧美日韩一区二区三区 | 国产视频一区二区在线看| 精品久久久精品久久久| 99九九在线精品视频| 美国免费a级毛片| 国产成人av激情在线播放| 人妻久久中文字幕网| 亚洲人成电影免费在线| 操出白浆在线播放| 中国美女看黄片| 99精国产麻豆久久婷婷| 在线观看免费日韩欧美大片| 精品熟女少妇八av免费久了| 51午夜福利影视在线观看| 国产精品熟女久久久久浪| 91av网站免费观看| 99国产综合亚洲精品| 国产三级黄色录像| 久久中文字幕一级| 欧美精品一区二区免费开放| 操美女的视频在线观看| 国产成人欧美在线观看 | 午夜成年电影在线免费观看| 国产男女内射视频| 男女高潮啪啪啪动态图| 国产精品一区二区免费欧美 | av欧美777| 女性生殖器流出的白浆| 黄频高清免费视频| av网站免费在线观看视频| 亚洲精品国产精品久久久不卡| 又黄又粗又硬又大视频| 黄色视频不卡| 十八禁网站免费在线| 久久亚洲国产成人精品v| 啦啦啦 在线观看视频| 久久av网站| 国产精品香港三级国产av潘金莲| 电影成人av| 超碰97精品在线观看| 久久精品国产a三级三级三级| 啦啦啦啦在线视频资源| 日本91视频免费播放| 欧美黄色片欧美黄色片| 成人免费观看视频高清| 五月开心婷婷网| 亚洲av日韩精品久久久久久密| 欧美亚洲 丝袜 人妻 在线| 久久人妻熟女aⅴ| 亚洲精品美女久久久久99蜜臀| 久久狼人影院| 制服人妻中文乱码| 黑人猛操日本美女一级片| 真人做人爱边吃奶动态| 欧美日韩国产mv在线观看视频| 脱女人内裤的视频| 亚洲精品成人av观看孕妇| 欧美黑人精品巨大| 在线看a的网站| 亚洲伊人久久精品综合| 一级毛片精品| 日韩 亚洲 欧美在线| 三上悠亚av全集在线观看| 精品亚洲成国产av| 亚洲国产日韩一区二区| 日韩人妻精品一区2区三区| 热99re8久久精品国产| 国产极品粉嫩免费观看在线| 欧美亚洲 丝袜 人妻 在线| 久久精品久久久久久噜噜老黄| 亚洲伊人久久精品综合| 99热国产这里只有精品6| 老熟女久久久| 国产成人系列免费观看| 精品国产一区二区三区四区第35| 久久精品国产亚洲av香蕉五月 | 午夜福利在线免费观看网站| 日韩熟女老妇一区二区性免费视频| 国产欧美亚洲国产| 在线看a的网站| 国产成人精品久久二区二区91| 一进一出抽搐动态| 亚洲精品乱久久久久久| 91精品国产国语对白视频| 久久国产亚洲av麻豆专区| 欧美激情久久久久久爽电影 | 黄色片一级片一级黄色片| 欧美日韩亚洲国产一区二区在线观看 | 午夜两性在线视频| 一区二区三区激情视频| 日韩欧美国产一区二区入口| 亚洲 国产 在线| 91字幕亚洲| 国产精品麻豆人妻色哟哟久久| 一级黄色大片毛片| 丝袜美腿诱惑在线| 一本一本久久a久久精品综合妖精| 久久亚洲精品不卡| 久久狼人影院| 国产精品二区激情视频| 在线观看www视频免费| a在线观看视频网站| 免费人妻精品一区二区三区视频| 亚洲欧美成人综合另类久久久| 国产麻豆69| 亚洲人成电影免费在线| 老司机福利观看| av网站免费在线观看视频| 少妇 在线观看| 国产xxxxx性猛交| 老司机影院毛片| 1024视频免费在线观看| 免费观看av网站的网址| 亚洲一区二区三区欧美精品| 午夜精品久久久久久毛片777| 99热网站在线观看| 老司机午夜福利在线观看视频 | 少妇人妻久久综合中文| 国产精品 欧美亚洲| 久久午夜综合久久蜜桃| 久久精品亚洲熟妇少妇任你| 十八禁高潮呻吟视频| 人妻一区二区av| 精品国产一区二区三区四区第35| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品自产拍在线观看55亚洲 | 国产精品 国内视频| 十八禁人妻一区二区| 国产高清videossex| 一边摸一边做爽爽视频免费| 宅男免费午夜| tocl精华| 50天的宝宝边吃奶边哭怎么回事| 一级a爱视频在线免费观看| 精品一区二区三卡| 高清黄色对白视频在线免费看| 中文字幕色久视频| 99久久精品国产亚洲精品| 久久精品成人免费网站| 免费女性裸体啪啪无遮挡网站| 日韩,欧美,国产一区二区三区| 久久久久精品人妻al黑| 久久女婷五月综合色啪小说| 成人亚洲精品一区在线观看| 91av网站免费观看| 三级毛片av免费| 咕卡用的链子| 亚洲专区中文字幕在线| 国产一区二区三区综合在线观看| 成年美女黄网站色视频大全免费| 在线观看免费午夜福利视频| 丝袜美腿诱惑在线| 999久久久国产精品视频| 精品福利观看| 热re99久久精品国产66热6| www.999成人在线观看| 亚洲人成电影观看| 亚洲国产中文字幕在线视频| 午夜两性在线视频| 亚洲精品一二三| 精品国产一区二区三区久久久樱花| 午夜福利免费观看在线| 在线观看免费午夜福利视频| xxxhd国产人妻xxx| 18禁观看日本| 国产深夜福利视频在线观看| 男人添女人高潮全过程视频| 捣出白浆h1v1| 亚洲人成电影观看| 国产成人精品无人区| 夜夜夜夜夜久久久久| 亚洲精品国产av成人精品| 不卡一级毛片| 99久久人妻综合| 人妻 亚洲 视频| av视频免费观看在线观看| 五月开心婷婷网| 最黄视频免费看| 国产亚洲午夜精品一区二区久久| 国产精品香港三级国产av潘金莲| 精品久久久久久久毛片微露脸 | 欧美日韩成人在线一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 丰满少妇做爰视频| 精品亚洲乱码少妇综合久久| 99国产精品免费福利视频| 国产成人精品久久二区二区91| 午夜成年电影在线免费观看| 成年人午夜在线观看视频| 欧美精品高潮呻吟av久久| 亚洲国产精品一区二区三区在线| 国产欧美日韩一区二区精品| 欧美大码av| 亚洲成人免费电影在线观看| 欧美黑人精品巨大| 母亲3免费完整高清在线观看| 91麻豆av在线| 在线十欧美十亚洲十日本专区| 欧美+亚洲+日韩+国产| 人人妻人人爽人人添夜夜欢视频| 亚洲精华国产精华精| 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲 | 欧美黄色片欧美黄色片| 91国产中文字幕| 脱女人内裤的视频| 每晚都被弄得嗷嗷叫到高潮| 国产不卡av网站在线观看| 最新的欧美精品一区二区| 人人妻人人添人人爽欧美一区卜| 色婷婷久久久亚洲欧美| 丝袜喷水一区| 欧美在线一区亚洲| 日本猛色少妇xxxxx猛交久久| 老熟妇乱子伦视频在线观看 | 亚洲欧美一区二区三区黑人| 一级片'在线观看视频| 亚洲欧美日韩另类电影网站| 男女午夜视频在线观看| 性色av乱码一区二区三区2| 99久久国产精品久久久| 王馨瑶露胸无遮挡在线观看| 十八禁高潮呻吟视频| 老司机影院毛片| 爱豆传媒免费全集在线观看| 无限看片的www在线观看| 亚洲国产看品久久| 男女下面插进去视频免费观看| 国产精品免费大片| 少妇猛男粗大的猛烈进出视频| 91大片在线观看| 一边摸一边做爽爽视频免费| av视频免费观看在线观看| 视频在线观看一区二区三区| 首页视频小说图片口味搜索| 久久精品亚洲熟妇少妇任你| 日韩电影二区| 亚洲性夜色夜夜综合| 久久精品久久久久久噜噜老黄| 真人做人爱边吃奶动态| 精品亚洲乱码少妇综合久久| 国产免费福利视频在线观看| 精品欧美一区二区三区在线| 色老头精品视频在线观看| 亚洲伊人色综图| 捣出白浆h1v1| 日本vs欧美在线观看视频| 亚洲欧美清纯卡通| 精品一区二区三卡| 久久免费观看电影| 国产日韩一区二区三区精品不卡| 久久久久久免费高清国产稀缺| 亚洲五月色婷婷综合| 亚洲精品久久午夜乱码| 欧美成人午夜精品| 久久精品熟女亚洲av麻豆精品| 日韩大码丰满熟妇| 国产1区2区3区精品| 妹子高潮喷水视频| 国产男女内射视频| 日本wwww免费看| 欧美乱码精品一区二区三区| 国产又爽黄色视频| 男女下面插进去视频免费观看| 女警被强在线播放| 日日爽夜夜爽网站| 亚洲成人免费av在线播放| 精品国产乱码久久久久久男人| 91麻豆av在线| 色视频在线一区二区三区| 国产欧美日韩一区二区精品| 免费观看人在逋| 国产精品 国内视频| 国产精品秋霞免费鲁丝片| a 毛片基地| 成年人免费黄色播放视频| 人妻久久中文字幕网| 在线亚洲精品国产二区图片欧美| 男女之事视频高清在线观看| 久久99热这里只频精品6学生| 成年av动漫网址| 在线看a的网站| 老司机福利观看| 99久久99久久久精品蜜桃| 亚洲 国产 在线| 大陆偷拍与自拍| 狂野欧美激情性xxxx| 欧美激情久久久久久爽电影 | 美女脱内裤让男人舔精品视频| 中文字幕人妻熟女乱码| 久久人人爽av亚洲精品天堂| 狠狠精品人妻久久久久久综合| av视频免费观看在线观看| 最新的欧美精品一区二区| 国产激情久久老熟女| 亚洲精品中文字幕在线视频| 亚洲视频免费观看视频| 亚洲黑人精品在线| 成人影院久久| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| 国产亚洲欧美精品永久| 黄色 视频免费看| 亚洲精华国产精华精| 亚洲中文日韩欧美视频| 另类亚洲欧美激情| 69精品国产乱码久久久| 亚洲欧美清纯卡通| 亚洲欧美色中文字幕在线| 永久免费av网站大全| 久久久久视频综合| netflix在线观看网站| 天天躁夜夜躁狠狠躁躁| 美女脱内裤让男人舔精品视频| 后天国语完整版免费观看| 青春草亚洲视频在线观看| 久久人人97超碰香蕉20202| av超薄肉色丝袜交足视频| 欧美久久黑人一区二区| 久久久久久久大尺度免费视频| 日韩视频一区二区在线观看| 可以免费在线观看a视频的电影网站| 成人av一区二区三区在线看 | 亚洲av成人一区二区三| 午夜福利乱码中文字幕| 国产成人av教育| 亚洲国产精品成人久久小说| 国产免费av片在线观看野外av| 久久久久久免费高清国产稀缺| 一级a爱视频在线免费观看| 精品人妻1区二区| 国产伦人伦偷精品视频| 久久影院123| 国产又爽黄色视频| 69精品国产乱码久久久| 免费看十八禁软件| 日本猛色少妇xxxxx猛交久久| 另类精品久久| 亚洲国产精品一区二区三区在线| 中国国产av一级| 国产真人三级小视频在线观看| 中文字幕人妻熟女乱码| 高清欧美精品videossex| 日本精品一区二区三区蜜桃| 日本a在线网址| 18禁黄网站禁片午夜丰满| 国产精品久久久久成人av| 久久毛片免费看一区二区三区| 淫妇啪啪啪对白视频 | 国产精品久久久久久人妻精品电影 | 51午夜福利影视在线观看| 美女午夜性视频免费| 久久免费观看电影| 亚洲国产欧美一区二区综合| 大片电影免费在线观看免费| 不卡一级毛片| 欧美xxⅹ黑人| 青青草视频在线视频观看| 韩国精品一区二区三区| 大陆偷拍与自拍| 亚洲精品国产一区二区精华液| 成人国产一区最新在线观看| 国产有黄有色有爽视频| 欧美另类一区| 极品少妇高潮喷水抽搐| 高清黄色对白视频在线免费看| 女人精品久久久久毛片| 制服人妻中文乱码| 亚洲精品国产区一区二| 日韩欧美免费精品| 久久久精品国产亚洲av高清涩受| 久久精品国产亚洲av高清一级| 免费在线观看日本一区| 另类精品久久| 日本精品一区二区三区蜜桃| 久久香蕉激情| 亚洲自偷自拍图片 自拍| 美女中出高潮动态图| 国产精品秋霞免费鲁丝片| 亚洲精品国产区一区二| 狠狠狠狠99中文字幕| 国产男人的电影天堂91| 丝瓜视频免费看黄片| 亚洲欧美精品自产自拍| 18禁裸乳无遮挡动漫免费视频| 黄网站色视频无遮挡免费观看| 亚洲av日韩在线播放| 黄色a级毛片大全视频| 亚洲 国产 在线| 久久99一区二区三区| 久久久久久久精品精品| 女人高潮潮喷娇喘18禁视频| 久久精品成人免费网站| 亚洲自偷自拍图片 自拍| 久久精品aⅴ一区二区三区四区| 成年人黄色毛片网站| 久久精品亚洲熟妇少妇任你| 99re6热这里在线精品视频| 日韩有码中文字幕| 国产精品久久久久久精品古装| 亚洲国产毛片av蜜桃av| 国产在线观看jvid| 69av精品久久久久久 | 可以免费在线观看a视频的电影网站| 午夜精品国产一区二区电影| 亚洲av欧美aⅴ国产| 色播在线永久视频| 久久亚洲国产成人精品v| 在线av久久热| 99国产极品粉嫩在线观看| 欧美精品一区二区大全| 黑人操中国人逼视频| 亚洲久久久国产精品| 夫妻午夜视频| 天天影视国产精品| 亚洲国产av新网站| 国产亚洲一区二区精品| 满18在线观看网站| 9色porny在线观看| a级毛片黄视频| 国产av精品麻豆| 久久国产亚洲av麻豆专区| 国产91精品成人一区二区三区 | 国产亚洲午夜精品一区二区久久| 在线观看一区二区三区激情| 亚洲精品日韩在线中文字幕| 亚洲自偷自拍图片 自拍| 99精国产麻豆久久婷婷| 精品少妇内射三级| 日韩,欧美,国产一区二区三区| 精品国产乱子伦一区二区三区 | 久久久久精品国产欧美久久久 | 午夜视频精品福利| 久久久久国内视频| 手机成人av网站| 99精品欧美一区二区三区四区| √禁漫天堂资源中文www| 国产精品麻豆人妻色哟哟久久| 超色免费av| 精品国产乱子伦一区二区三区 | 国产人伦9x9x在线观看| 男人添女人高潮全过程视频| 爱豆传媒免费全集在线观看| 国产欧美日韩一区二区三区在线| 两个人免费观看高清视频| 亚洲伊人久久精品综合| 精品欧美一区二区三区在线| 99久久人妻综合| av有码第一页| 日韩人妻精品一区2区三区| 亚洲午夜精品一区,二区,三区| 亚洲欧美精品自产自拍| 中文字幕人妻熟女乱码| 婷婷色av中文字幕| 在线观看免费日韩欧美大片| 操美女的视频在线观看| 日韩欧美国产一区二区入口| 黄色片一级片一级黄色片| 亚洲国产欧美日韩在线播放| 成在线人永久免费视频| 免费一级毛片在线播放高清视频 | 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久成人av| 91成年电影在线观看| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 色老头精品视频在线观看| 国产一区二区激情短视频 | 亚洲综合色网址| 国产欧美日韩一区二区三区在线| 男女免费视频国产| 咕卡用的链子| 俄罗斯特黄特色一大片| 老司机福利观看| 久久狼人影院| 欧美精品人与动牲交sv欧美| 在线十欧美十亚洲十日本专区| 国产免费福利视频在线观看| 欧美精品啪啪一区二区三区 | 欧美黄色片欧美黄色片| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三区在线| 国产熟女午夜一区二区三区| 午夜两性在线视频| 在线观看免费日韩欧美大片| 久久久久久免费高清国产稀缺| 免费高清在线观看视频在线观看| 国产视频一区二区在线看| 久久免费观看电影| 黄色视频在线播放观看不卡| 久久亚洲精品不卡| 欧美黑人精品巨大| 男人操女人黄网站| 丝袜脚勾引网站| 国产精品久久久久久人妻精品电影 | 久久久国产精品麻豆| 国产高清国产精品国产三级| 精品亚洲成a人片在线观看| 精品国产乱码久久久久久小说| 欧美日韩亚洲国产一区二区在线观看 | 国产成人av激情在线播放| 91精品三级在线观看| 中文精品一卡2卡3卡4更新| 搡老岳熟女国产| 国产成人系列免费观看| 亚洲第一欧美日韩一区二区三区 | 老司机影院成人| 久久久久久久国产电影| 亚洲国产av新网站| 国产成人系列免费观看| 黄色视频不卡| 久久国产亚洲av麻豆专区| 91成年电影在线观看| 国产成+人综合+亚洲专区| 51午夜福利影视在线观看| 亚洲专区国产一区二区| 黄片大片在线免费观看| 另类精品久久| 女警被强在线播放| av免费在线观看网站| 欧美黑人欧美精品刺激| 精品国内亚洲2022精品成人 | 免费久久久久久久精品成人欧美视频| 巨乳人妻的诱惑在线观看| 精品国产一区二区久久| 女人高潮潮喷娇喘18禁视频| 成人国语在线视频| 激情视频va一区二区三区| 免费人妻精品一区二区三区视频| 秋霞在线观看毛片| 国产精品偷伦视频观看了| 亚洲国产成人一精品久久久| 久久久精品国产亚洲av高清涩受| 亚洲精品一区蜜桃| 国产成人免费观看mmmm| 国产又爽黄色视频| 国产亚洲欧美精品永久| 一区福利在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲精品久久成人aⅴ小说| 欧美日韩亚洲国产一区二区在线观看 | 国产精品九九99| 久久久久久久久久久久大奶| 少妇 在线观看| 麻豆国产av国片精品| 嫩草影视91久久| 欧美精品一区二区大全| 啦啦啦啦在线视频资源| 美女大奶头黄色视频| 国产精品一二三区在线看| 免费在线观看黄色视频的| 久久久精品94久久精品| 午夜福利,免费看| 在线观看免费高清a一片| 成人免费观看视频高清| 亚洲自偷自拍图片 自拍| 中国美女看黄片| 99热国产这里只有精品6| 韩国精品一区二区三区| 欧美黄色淫秽网站| 久久精品熟女亚洲av麻豆精品| 欧美国产精品一级二级三级| 国产又色又爽无遮挡免| 91老司机精品| 国产亚洲欧美在线一区二区| 中文字幕精品免费在线观看视频| 亚洲七黄色美女视频| 9热在线视频观看99| 交换朋友夫妻互换小说| 久久人人爽人人片av| 欧美日韩一级在线毛片| 午夜影院在线不卡|