• <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片天天在线观看| 变态另类成人亚洲欧美熟女 | 在线观看免费视频日本深夜| 91成人精品电影| 天堂√8在线中文| 这个男人来自地球电影免费观看| www.精华液| 欧美日韩国产mv在线观看视频| 久久婷婷成人综合色麻豆| 精品一区二区三卡| 亚洲一码二码三码区别大吗| 首页视频小说图片口味搜索| 日韩免费高清中文字幕av| 大型黄色视频在线免费观看| 国产精品.久久久| 欧美乱妇无乱码| 午夜成年电影在线免费观看| 老司机在亚洲福利影院| 国产精品香港三级国产av潘金莲| 久久99一区二区三区| 国产成人免费观看mmmm| 无限看片的www在线观看| 韩国精品一区二区三区| 午夜福利在线免费观看网站| 日本撒尿小便嘘嘘汇集6| 亚洲精品av麻豆狂野| 老司机靠b影院| 亚洲,欧美精品.| av超薄肉色丝袜交足视频| 欧美一级毛片孕妇| 中文字幕精品免费在线观看视频| 精品亚洲成a人片在线观看| 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 9191精品国产免费久久| 午夜影院日韩av| 一个人免费在线观看的高清视频| 99久久精品国产亚洲精品| 91精品国产国语对白视频| 国产一区二区三区视频了| 91麻豆精品激情在线观看国产 | 欧美另类亚洲清纯唯美| 亚洲专区国产一区二区| 久久亚洲精品不卡| 久久精品成人免费网站| 丰满迷人的少妇在线观看| 国产精品1区2区在线观看. | 欧美日韩亚洲国产一区二区在线观看 | 国产成人免费无遮挡视频| 亚洲精品在线美女| 精品久久久久久,| 成年人免费黄色播放视频| 制服诱惑二区| 免费一级毛片在线播放高清视频 | 777久久人妻少妇嫩草av网站| 91成人精品电影| 日韩视频一区二区在线观看| 夜夜夜夜夜久久久久| 9191精品国产免费久久| 亚洲精品在线观看二区| 午夜福利免费观看在线| 美女 人体艺术 gogo| 久久国产亚洲av麻豆专区| 美女视频免费永久观看网站| 99re6热这里在线精品视频| 亚洲五月婷婷丁香| 久久 成人 亚洲| 999久久久国产精品视频| 美女高潮到喷水免费观看| 99久久国产精品久久久| 黄色毛片三级朝国网站| 在线视频色国产色| 男女下面插进去视频免费观看| 老司机午夜福利在线观看视频| 欧美人与性动交α欧美软件| 两个人免费观看高清视频| 亚洲黑人精品在线| 久久影院123| 中文亚洲av片在线观看爽 | 男男h啪啪无遮挡| 亚洲全国av大片| 国产亚洲欧美在线一区二区| 两人在一起打扑克的视频| 黑人操中国人逼视频| 婷婷精品国产亚洲av在线 | 国产精品1区2区在线观看. | 中文字幕最新亚洲高清| 亚洲精品一卡2卡三卡4卡5卡| av欧美777| 国产精品久久视频播放| 久久影院123| 建设人人有责人人尽责人人享有的| 日韩 欧美 亚洲 中文字幕| 亚洲av日韩精品久久久久久密| 一进一出抽搐动态| 好看av亚洲va欧美ⅴa在| 午夜影院日韩av| 久久精品国产亚洲av高清一级| 成熟少妇高潮喷水视频| 亚洲中文字幕日韩| av有码第一页| 精品欧美一区二区三区在线| 亚洲午夜精品一区,二区,三区| 80岁老熟妇乱子伦牲交| 欧美日韩视频精品一区| 在线观看舔阴道视频| 动漫黄色视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产一区二区三区视频了| 日韩三级视频一区二区三区| 国产亚洲欧美98| 欧美人与性动交α欧美软件| a级片在线免费高清观看视频| 桃红色精品国产亚洲av| 国精品久久久久久国模美| 男人舔女人的私密视频| 性色av乱码一区二区三区2| 精品欧美一区二区三区在线| 老司机在亚洲福利影院| 亚洲成人国产一区在线观看| bbb黄色大片| 在线国产一区二区在线| 国产又色又爽无遮挡免费看| 久久九九热精品免费| 免费一级毛片在线播放高清视频 | 中文亚洲av片在线观看爽 | 国产成+人综合+亚洲专区| 国产亚洲精品久久久久久毛片 | 丁香六月欧美| 久久狼人影院| tocl精华| 免费日韩欧美在线观看| 美女高潮喷水抽搐中文字幕| 亚洲熟女精品中文字幕| 热99国产精品久久久久久7| 日韩精品免费视频一区二区三区| 最近最新免费中文字幕在线| 欧美最黄视频在线播放免费 | av在线播放免费不卡| 高清毛片免费观看视频网站 | 久久精品亚洲熟妇少妇任你| 亚洲成人国产一区在线观看| 精品一区二区三区四区五区乱码| 少妇 在线观看| 捣出白浆h1v1| 精品久久蜜臀av无| 亚洲精品在线美女| 日韩视频一区二区在线观看| 国产高清国产精品国产三级| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 免费av中文字幕在线| 久久久精品区二区三区| 欧美日韩亚洲综合一区二区三区_| 熟女少妇亚洲综合色aaa.| 一区二区三区激情视频| 欧美日韩亚洲综合一区二区三区_| 国产精品免费大片| 韩国精品一区二区三区| 丝袜人妻中文字幕| 咕卡用的链子| 久久精品熟女亚洲av麻豆精品| 青草久久国产| 中文字幕精品免费在线观看视频| 精品人妻熟女毛片av久久网站| 国产麻豆69| 国产精品乱码一区二三区的特点 | av在线播放免费不卡| 一区二区三区精品91| 啦啦啦视频在线资源免费观看| 精品亚洲成国产av| 久久精品91无色码中文字幕| 日日夜夜操网爽| 欧美精品高潮呻吟av久久| a在线观看视频网站| 啦啦啦在线免费观看视频4| 午夜影院日韩av| 免费看十八禁软件| 欧美人与性动交α欧美软件| 亚洲精品粉嫩美女一区| 黄色丝袜av网址大全| 亚洲欧美精品综合一区二区三区| 90打野战视频偷拍视频| 色尼玛亚洲综合影院| 在线观看66精品国产| 欧美日本中文国产一区发布| 少妇的丰满在线观看| 一区二区三区激情视频| 亚洲一码二码三码区别大吗| 制服人妻中文乱码| 丁香欧美五月| 久久精品亚洲av国产电影网| 免费在线观看视频国产中文字幕亚洲| 精品人妻熟女毛片av久久网站| 亚洲精品乱久久久久久| ponron亚洲| 人成视频在线观看免费观看| 看黄色毛片网站| 久久久精品免费免费高清| 国产熟女午夜一区二区三区| 水蜜桃什么品种好| 欧美激情极品国产一区二区三区| 黄色丝袜av网址大全| 韩国精品一区二区三区| 久久精品亚洲av国产电影网| 下体分泌物呈黄色| 久久中文字幕一级| 久久国产精品人妻蜜桃| 巨乳人妻的诱惑在线观看| 大香蕉久久成人网| 9191精品国产免费久久| a级片在线免费高清观看视频| 欧美日韩福利视频一区二区| 黄色丝袜av网址大全| netflix在线观看网站| 国产精品秋霞免费鲁丝片| 国产精品一区二区免费欧美| 亚洲成人免费电影在线观看| 一进一出好大好爽视频| 久久久久国内视频| 国产精品一区二区在线不卡| av线在线观看网站| 18禁美女被吸乳视频| 午夜成年电影在线免费观看| 超碰成人久久| 少妇被粗大的猛进出69影院| 一级作爱视频免费观看| 啦啦啦在线免费观看视频4| 国产深夜福利视频在线观看| 少妇裸体淫交视频免费看高清 | 99国产精品免费福利视频| 热99re8久久精品国产| 久久精品亚洲熟妇少妇任你| 他把我摸到了高潮在线观看| 大香蕉久久成人网| 国产精品电影一区二区三区 | 日本撒尿小便嘘嘘汇集6| 岛国在线观看网站| 亚洲国产毛片av蜜桃av| 俄罗斯特黄特色一大片| 欧美日韩一级在线毛片| 在线观看免费高清a一片| 成年动漫av网址| 国产成人一区二区三区免费视频网站| 老汉色∧v一级毛片| 天天操日日干夜夜撸| 99国产精品99久久久久| 人人妻人人爽人人添夜夜欢视频| 在线播放国产精品三级| 久久精品熟女亚洲av麻豆精品| 视频区图区小说| 麻豆乱淫一区二区| 色精品久久人妻99蜜桃| 手机成人av网站| 岛国毛片在线播放| 亚洲成国产人片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 黄色视频不卡| 久久草成人影院| 亚洲熟女毛片儿| 亚洲,欧美精品.| 成人国语在线视频| 悠悠久久av| 亚洲黑人精品在线| 欧美日韩瑟瑟在线播放| 欧美日韩乱码在线| 中文字幕av电影在线播放| 午夜免费观看网址| 黑人猛操日本美女一级片| 色在线成人网| 夜夜夜夜夜久久久久| 欧美 日韩 精品 国产| 1024视频免费在线观看| 午夜福利影视在线免费观看| 亚洲精品国产精品久久久不卡| 国产人伦9x9x在线观看| 精品人妻在线不人妻| 99热网站在线观看| 国产精品综合久久久久久久免费 | 久久精品国产亚洲av香蕉五月 | 亚洲,欧美精品.| 成熟少妇高潮喷水视频| 国产在线观看jvid| 久久久国产成人精品二区 | 午夜福利免费观看在线| 国产1区2区3区精品| 亚洲,欧美精品.| 国产激情久久老熟女| 亚洲美女黄片视频| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人巨大精品欧美一区二区mp4| 18禁美女被吸乳视频| 亚洲国产欧美网| 天堂√8在线中文| 国产高清videossex| 午夜福利欧美成人| 国产成人欧美在线观看 | 色综合欧美亚洲国产小说| 18禁观看日本| 久久久国产成人免费| 少妇 在线观看| 满18在线观看网站| 精品久久久精品久久久| 亚洲av日韩精品久久久久久密| 黄片播放在线免费| 深夜精品福利| 精品一区二区三区av网在线观看| 一进一出好大好爽视频| 欧美精品av麻豆av| 91麻豆精品激情在线观看国产 | 免费日韩欧美在线观看| 亚洲视频免费观看视频| 精品无人区乱码1区二区| 黄色成人免费大全| 国产精品av久久久久免费| 最新美女视频免费是黄的| 色在线成人网| 国产欧美日韩一区二区精品| 热99国产精品久久久久久7| 99在线人妻在线中文字幕 | av线在线观看网站| 在线观看日韩欧美| 丝袜在线中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩黄片免| 十八禁人妻一区二区| 国产精品1区2区在线观看. | www日本在线高清视频| 国产男女超爽视频在线观看| 91九色精品人成在线观看| 岛国毛片在线播放| 在线观看免费高清a一片| a级片在线免费高清观看视频| 免费av中文字幕在线| 午夜福利影视在线免费观看| 国产精品.久久久| xxxhd国产人妻xxx| 精品人妻熟女毛片av久久网站| 成人影院久久| 国产不卡av网站在线观看| 欧美国产精品一级二级三级| 久久久精品免费免费高清| 十分钟在线观看高清视频www| 91字幕亚洲| 91国产中文字幕| 9热在线视频观看99| 男女下面插进去视频免费观看| av不卡在线播放| 99国产综合亚洲精品| 9191精品国产免费久久| 欧美激情高清一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 久久人妻福利社区极品人妻图片| 两性午夜刺激爽爽歪歪视频在线观看 | 国产男女超爽视频在线观看| 法律面前人人平等表现在哪些方面| 亚洲色图 男人天堂 中文字幕| 欧美中文综合在线视频| 国产精品乱码一区二三区的特点 | 两人在一起打扑克的视频| 又黄又粗又硬又大视频| 纯流量卡能插随身wifi吗| 午夜福利免费观看在线| 欧美丝袜亚洲另类 | 老司机靠b影院| 满18在线观看网站| 在线观看舔阴道视频| 国产三级黄色录像| 精品第一国产精品| 午夜福利欧美成人| 极品人妻少妇av视频| 久久婷婷成人综合色麻豆| 天堂中文最新版在线下载| 久久午夜综合久久蜜桃| 久久国产乱子伦精品免费另类| 麻豆av在线久日| 美女国产高潮福利片在线看| 女人久久www免费人成看片| 正在播放国产对白刺激| 精品一区二区三卡| 美女 人体艺术 gogo| 欧美 亚洲 国产 日韩一| 久久久久精品人妻al黑| 国产片内射在线| 亚洲专区中文字幕在线| 国产精品 国内视频| 两性午夜刺激爽爽歪歪视频在线观看 | 王馨瑶露胸无遮挡在线观看| 亚洲精品一二三| 国产97色在线日韩免费| 一区二区三区激情视频| 老鸭窝网址在线观看| 精品高清国产在线一区| 少妇被粗大的猛进出69影院| 亚洲熟女精品中文字幕| 美女国产高潮福利片在线看| www.999成人在线观看| 777米奇影视久久| 老司机靠b影院| a在线观看视频网站| 无限看片的www在线观看| 99国产精品免费福利视频| 99久久综合精品五月天人人| 国产高清视频在线播放一区| 99国产精品一区二区蜜桃av | 亚洲成人国产一区在线观看| av网站在线播放免费| 亚洲少妇的诱惑av| 操出白浆在线播放| 精品福利观看| 国内毛片毛片毛片毛片毛片| 久久 成人 亚洲| 日韩欧美一区二区三区在线观看 | 女同久久另类99精品国产91| 母亲3免费完整高清在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品粉嫩美女一区| 精品国产乱子伦一区二区三区| 精品一品国产午夜福利视频| 日本wwww免费看| 女人被躁到高潮嗷嗷叫费观| 一级,二级,三级黄色视频| 欧美性长视频在线观看| 黄色视频,在线免费观看| 黄色毛片三级朝国网站| 亚洲av成人av| 久久人妻av系列| 国产免费av片在线观看野外av| 中文字幕色久视频| 日本五十路高清| 日本一区二区免费在线视频| 激情视频va一区二区三区| 日韩欧美国产一区二区入口| 日韩欧美一区二区三区在线观看 | 真人做人爱边吃奶动态| 亚洲男人天堂网一区| 亚洲精品国产色婷婷电影| 国产真人三级小视频在线观看| 国产亚洲欧美98| 一进一出抽搐gif免费好疼 | 精品国产一区二区久久| 亚洲精品国产区一区二| 在线观看午夜福利视频| 精品久久久精品久久久| 午夜福利,免费看| www.精华液| av中文乱码字幕在线| 99香蕉大伊视频| 亚洲午夜精品一区,二区,三区| 国产男女超爽视频在线观看| 人妻丰满熟妇av一区二区三区 | 18禁裸乳无遮挡免费网站照片 | 亚洲国产精品合色在线| videosex国产| 精品少妇久久久久久888优播| 亚洲片人在线观看| bbb黄色大片| 亚洲欧美日韩另类电影网站| 老司机影院毛片| 国精品久久久久久国模美| 午夜福利视频在线观看免费| 一级,二级,三级黄色视频| 亚洲自偷自拍图片 自拍| 午夜福利一区二区在线看| 欧美丝袜亚洲另类 | 最新美女视频免费是黄的| 女同久久另类99精品国产91| 男女之事视频高清在线观看| 中文字幕最新亚洲高清| 亚洲av成人不卡在线观看播放网| 色精品久久人妻99蜜桃| 精品高清国产在线一区| 99国产极品粉嫩在线观看| 男女下面插进去视频免费观看| 国产成人精品久久二区二区免费| 久久精品熟女亚洲av麻豆精品| 男人操女人黄网站| 好看av亚洲va欧美ⅴa在| 在线视频色国产色| 久久亚洲真实| 一本综合久久免费| 国产一卡二卡三卡精品| 日韩欧美在线二视频 | 久9热在线精品视频| 又黄又粗又硬又大视频| 亚洲精品一二三| 久久这里只有精品19| 日韩中文字幕欧美一区二区| 色94色欧美一区二区| 国产精品免费大片| 新久久久久国产一级毛片| 国产精品成人在线| 在线观看免费视频日本深夜| 黄片大片在线免费观看| 动漫黄色视频在线观看| 91国产中文字幕| 女同久久另类99精品国产91| а√天堂www在线а√下载 | 黑丝袜美女国产一区| 亚洲色图综合在线观看| 日韩欧美三级三区| 亚洲精品久久午夜乱码| 亚洲第一青青草原| 夜夜躁狠狠躁天天躁| 嫩草影视91久久| 性少妇av在线| 9191精品国产免费久久| 亚洲国产欧美日韩在线播放| 夜夜爽天天搞| 一本综合久久免费| 欧美日韩乱码在线| 又黄又粗又硬又大视频| 久久香蕉国产精品| 91大片在线观看| 热99久久久久精品小说推荐| 国产单亲对白刺激| 精品亚洲成a人片在线观看| 国产成人一区二区三区免费视频网站| 欧美日韩黄片免| 久久精品成人免费网站| 一级黄色大片毛片| 国产亚洲精品久久久久5区| 亚洲精品久久午夜乱码| av视频免费观看在线观看| 精品欧美一区二区三区在线| 91字幕亚洲| 国产精品成人在线| 老熟妇乱子伦视频在线观看| 少妇裸体淫交视频免费看高清 | 捣出白浆h1v1| 亚洲欧美日韩高清在线视频| 欧美在线黄色| x7x7x7水蜜桃| 欧美+亚洲+日韩+国产| 国产成人免费观看mmmm| 久久久久久久国产电影| 精品国产一区二区久久| 色在线成人网| 国产精品免费大片| 18禁国产床啪视频网站| 国产成人av激情在线播放| 大香蕉久久网| 天天躁日日躁夜夜躁夜夜| 91av网站免费观看| 看黄色毛片网站| 51午夜福利影视在线观看| 成年版毛片免费区| 夜夜爽天天搞| 一级黄色大片毛片| 中亚洲国语对白在线视频| 国产精品98久久久久久宅男小说| 男人的好看免费观看在线视频 | 亚洲欧美激情综合另类| videos熟女内射| 91麻豆av在线| 少妇被粗大的猛进出69影院| 丰满迷人的少妇在线观看| 老汉色av国产亚洲站长工具| 真人做人爱边吃奶动态| 大陆偷拍与自拍| 久久精品亚洲熟妇少妇任你| av网站免费在线观看视频| 国产亚洲av高清不卡| 国产av又大| 免费在线观看影片大全网站| 中文字幕人妻熟女乱码| 丝瓜视频免费看黄片| 老司机午夜十八禁免费视频| 午夜福利视频在线观看免费| 高清av免费在线| 91麻豆精品激情在线观看国产 | 欧美最黄视频在线播放免费 | 国产精品1区2区在线观看. | 亚洲成人手机| 免费观看人在逋| 丝袜美足系列| 狠狠狠狠99中文字幕| 亚洲五月色婷婷综合| 国产精华一区二区三区| 久久久精品区二区三区| 大香蕉久久成人网| 精品一区二区三区视频在线观看免费 | 欧美 亚洲 国产 日韩一| 亚洲国产欧美一区二区综合| 在线观看66精品国产| 欧美日韩视频精品一区| 亚洲av成人不卡在线观看播放网| 国产成人欧美在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 人妻一区二区av| 看黄色毛片网站| 亚洲精品中文字幕一二三四区| 一级黄色大片毛片| 午夜精品久久久久久毛片777| 极品少妇高潮喷水抽搐| 美女扒开内裤让男人捅视频| 黄片小视频在线播放| 啦啦啦在线免费观看视频4| 国产成人精品久久二区二区91| 老熟女久久久| 淫妇啪啪啪对白视频| 亚洲av熟女| 午夜免费成人在线视频| 高潮久久久久久久久久久不卡| 一级片免费观看大全| 色在线成人网| 啦啦啦视频在线资源免费观看| 国产成人影院久久av| 国产真人三级小视频在线观看| 丝袜美足系列| 制服诱惑二区| 国产精品久久久人人做人人爽| 精品久久久精品久久久| 中文字幕人妻熟女乱码|