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

    Identification of portal vein tumor thrombus with an independent clonal origin in hepatocellular carcinoma via multi-omics data analysis

    2019-03-23 05:16:52ShupengLiuZaixinZhouYinJiaJieXueZhiyongLiuKaiChengShuqunChengShanrongLiu
    Cancer Biology & Medicine 2019年1期

    Shupeng Liu, Zaixin Zhou, Yin Jia, Jie Xue, Zhiyong Liu, Kai Cheng, Shuqun Cheng, Shanrong Liu

    1Clinical Research Center, Changhai Hospital, Second Military Medical University, Shanghai 200433, China; 2Department of Laboratory Diagnostics, Changhai Hospital, Second Military Medical University, Shanghai 200433, China; 3Department of Hepatic Surgery VI, Eastern Hepatobiliary Surgery Hospital, Second Military Medical University, Shanghai 200433, China

    ABSTRACT Objective:Multiple mechanisms underlying the development of portal vein tumor thrombus (PVTT) in hepatocellular carcinoma(HCC) have been reported recently. However, the origins of PVTT remain unknown. Increasing multi-omics data on PVTTs in HCCs have made it possible to investigate whether PVTTs originate from the corresponding primary tumors (Ts).Methods:The clonal relationship between PVTTs and their corresponding primary Ts was investigated using datasets deposited in public databases. One DNA copy number variations dataset and three gene expression datasets were downloaded for the analyses.Clonality analysis was performed to investigate the clonal relationship between PVTTs and Ts from an individual patient.Differential gene expression analysis was applied to investigate the gene expression profiles of PVTTs and Ts.Results:One out of 19 PVTTs had no clonal relationship with its corresponding T, whereas the others did. The PVTTs with independent clonal origin showed different gene expression and enrichment in biological processes from the primary Ts. Based on the unique gene expression profiles, a gene signature including 24 genes was used to identify pairs of PVTTs and primary Ts without any clonal relationship. Validation in three datasets showed that these types of pairs of PVTTs and Ts can be identified by the 24-gene signature.Conclusions:Our findings show a direct evidence for PVTT origin and consolidate the heterogeneity of PVTTs observed in clinic.The results suggest that PVTT investigation at a molecular level is clinically necessary for diagnosis and treatment.

    KEYWORDS Hepatocellular carcinoma; portal vein tumor thrombus; clonal origin; copy number variation; bioinformatics

    Introduction

    Portal vein tumor thrombus (PVTT) is one of the main hepatocellular carcinoma (HCC) complications, found in approximately 50% of HCC patients1. The prognosis of HCC patients with PVTT is poor, and the median overall survival(OS) time of patients with PVTT was primarily less than 1 year2,3. Although clinicians in Europe and America recommend sorafenib alone, and multidisciplinary therapy has been adapted in Asia4,5, no worldwide consensus or guidelines have shown effectiveness for specific PVTT treatment, which might be due to a lack of knowledge about the nature of PVTT.

    Although currently, hypoxia, non-coding RNAs, and cancer stem cells have been found to contribute to PVTT development6-9, the clonal origins of PVTT still remain unknown. The identical gene expression profiles between PVTT and the corresponding primary HCCs demonstrated that PVTT could originate from HCC primary nodules by metastasis10,11. Contradictorily, in small HCC patients, some PVTTs are clinically defined, that are distant from the liver parenchyma tumor (PT) nodules, and some PVTTs are found without liver parenchyma tumor nodules12,13.Moreover, some PVTTs showed significantly different gene expression profiles compared with those of the paired primary tumors (Ts)11. These findings indicate that PVTTs display relatively high heterogeneity and may have different clonal origins from the paired primary Ts.

    Recently, high-throughput technology, such as microarray and sequencing, has been widely used to investigate the mechanisms underlying the development of Ts, including HCC. Among them, 3 published studies involved the genomic or transcriptomic datasets of HCC as well as PVTT9,11,14. In the present study, previously published datasets (GSE77509, GSE69164, GSE77275, and GSE74656),and freshly obtained samples were used for analysis, and we found a type of PVTT that has a different clonal origin from that of the corresponding primary Ts in HCC. Additionally,this type of PVTT was found to express a specific set of genes,which varied from that of the corresponding Ts.

    Materials and methods

    Data resources

    Four previously published datasets (Cytoscan HD array,GSE77275; RNA-seq, GSE77509 and GSE69164; and cDNA microarray, GSE74656) used in this study were downloaded from the NCBI Gene Expression Omnibus (GEO,http://www.ncbi.nlm.nih.gov/geo/). GSE77275 included the copy number variation (CNV) profiles of 60 matched samples (PVTT/T/PT from 20 patients) assessed by Cytoscan HD array. GSE77509 included the gene expression profiles of these 60 samples assessed by RNA-seq. DNA and total RNA from the 60 samples were used for CNV and gene expression analysis, respectively, as previously reported11.

    The samples of patient 14 (P14), in these two datasets,were not included in the present study because they were reported to be contaminated11. GSE69164 included the gene expression of 33 matched samples from 11 patients assessed by RNA-seq14. GSE74656 included the gene expression data of 15 matched samples from 5 patients assessed by cDNA analysis.

    Clonality testing

    To investigate the clonal relationship between PVTTs and Ts from an individual patient, the CNV profiles of 19 pairs of PVTTs and Ts were extracted from GSE77275 using R package Affy2sv and loaded to R package Clonality for further analysis. The CNV of the chromosome region in each patient was estimated using the circular binary segmentation(CBS) analysis based on the normalized weighted Log2ratio.A significant CNV in chromosome regions was determined according to the threshold value (median absolute deviations = 1.25). The clonal relationship between PVTTs and Ts from an individual patient was investigated according to the likelihood ratio 2 (LR2)15,16.

    Differential gene expression analysis

    The gene expression profiles of PVTTs and Ts from the datasets GSE77509, GSE69164, and GSE74656 were extracted using the appropriate methods11. EBSeq from NovelBrain BioCloud (https://cloud.novelbrain.com) was applied to analyze differentially expressed genes between PVTTs and Ts from individual patients in GSE77509. Genes with fold changes > 3 and FDR < 0.05 were identified as significantly differentially expressed.

    Gene ontology and KEGG pathway enrichment analyses

    The differentially expressed genes between PVTT13 and T13 were loaded to Metascape for Gene ontology (GO) and KEGG pathway enrichment analyses17. Protein-protein interaction (PPI) network analysis was also performed in Metascape, and the PPI network was exported using Cytoscape 3.4.

    Gene signature enrichment analysis

    Gene expression signature enrichment was evaluated using single-sample gene set enrichment analysis (ssGSEA) and nearest template prediction (NTP) analyses (Gene Pattern modules)18. Molecular Signature Database gene sets,(MSigDB, www.broadinstitute.org/msigdb ) including hematopoietic stem cells and liver cancer stem cells, and previously published gene-expression signatures representing different cell types including bile duct cells (BDC),hepatocytes, hepatic progenitor cells (HPC), and hepatic stellate cells (HSC)19(Supplementary Table S1), were tested using ssGSEA. Previously, published HCC molecular classifications were also analyzed using NTP analysis20,21.

    Tumor tissues

    Eight pairs of human PVTTs and Ts were obtained from HCC patients who underwent curative resection in the Eastern Hepatobiliary Surgery Hospital (Shanghai, China)after obtaining informed consent. All 8 pairs of PVTTs and Ts were used to perform real-time polymerase chain reaction(PCR) analysis. The study was approved by the Ethics Committee of Eastern Hepatobiliary Surgery Hospital.

    Quantitative real-time PCR

    Total RNA was isolated using TRIzol reagent (Invitrogen,CA, USA), and reverse-transcription (RT) was performed using the PrimeScript RT Reagent Kit (Takara, Kusatsu,Japan), according to manufacturer's instructions. Real-time PCR was performed with SYBR Premix Ex Taq (Takara,Kusatsu, Japan) using the StepOnePlus Real-Time PCR System (Applied Biosystems, CA, USA). The primers used are listed in Supplementary Table S2. The gene expression levels were calculated relative to the expression of β-actin.

    Cluster analysis

    All pairs of PVTT and Ts were classified using one minus Pearson correlation hierarchical cluster analysis by the 24-gene signatures (Morpheus). The 24 genes are listed in Supplementary Table S3.

    Statistical analysis

    All P values and FDR values were obtained by the appropriate methods used in the different gene expression and enrichment analyses. Statistical analysis of the gene signature enrichment score from ssGSEA was not performed due to only one PVTT-T pair identified to have no clonal relationship.

    Results

    Analysis of the clonal relationship between PVTTs and corresponding primary Ts

    To investigate the clonal relationship of PVTTs with corresponding Ts, DNA copy number variation (CNV)profiles of 19 paired PVTTs and Ts were analyzed using the R package Clonality22. The clonal-relatedness of the paired tissues was evaluated according to the LR2 value. Overall, 18 out of 19 PVTT-T pairs showed a significantly higher LR2 value (P < 0.001), whereas only one pair (patient 13, P13)showed a much lower LR2 value (P > 0.05) (Table 1). As a low LR2 value indicates no clonal relatedness, PVTT13 was identified to have an independent clonal origin from T13,whereas the other PVTTs had a clonal origin from the corresponding Ts. For further verification, the CNV of the chromosome regions and genes was investigated using the R package DNA copy23. The DNA CNV patterns of chromosome regions were estimated based on the threshold values (threshold = 0.3). Aberrant chromosome regions showed very different CNV patterns between PVTT13 and T13, whereas these differences were not observed in other pairs (Figure 1A). The region of 2q24.1-q31.1 was deleted in PVTT13 but amplified in T13. Regions of 5q13.2-q35.2 and 15q11.2-q21.1 were amplified in PVTT13 but deleted in T13.These three chromosome regions displayed similar CNV patterns between PVTTs and Ts from other patients (Figure 1A). CNV analysis of genes showed that 10 HCC driver genes have different copy numbers between the PVTTs and Ts from individual patients (Figure 1B). Less than 3 different mutant genes were observed in most PVTT-T pairs (n = 15),whereas 3 or more different mutant genes were observed in pairs from patient 6 (P6), patient 11 (P11), P13, and patient 22 (P22) (Figure 1B). Among these 4 PVTT-T pairs with more than 2 mutant genes with different copy numbers, only the pair from P13 had 2 mutant genes, with 2 and 3 different copies. Five copies of ADCY2 and 3 copies of CTNNB1 were observed in PVTT13, whereas 2 copies of ADCY2 and 5 copies of CTNNB1 were observed in T13 (Figure 1B).Additionally, Hoshida S1-S2-S3 classes of the PVTTs and Ts were evaluated by NTP analysis using their gene expression profiles20. All PVTTs and Ts were classified into one of the Hoshida S1-S2-S3 classes (Bonferroni P = 0.019, FDR < 0.01,Figure 1C). Among the 19 PVTT-T pairs, 6 pairs including PVTT13-T13 had their PVTTs and Ts in different subclasses(Figure 1C). PVTT13 was classified into the S3 class, which is associated with hepatocyte differentiation, and T13 was classified into the S1 class, which reflects an aberrant activation of the Wnt-signaling pathway (Figure 1C).Collectively, these data indicated that PVTT13 and T13 had different origins and that some PVTTs may not originate from the primary Ts.

    Table 1 Clonality analysis of PVTT and T (n = 19) according to copy number variations at the probe level

    Analysis of aberrantly expressed genes in PVTT13

    Figure 1 Analysis of the clonal relationship between PVTTs and corresponding primary Ts. (A) Copy number variations (CNVs) of chromosome regions in PVTTs and corresponding primary Ts. The specific CNV regions in PVTT13 and T13 were labeled by rectangles. The CNV threshold = 0.3; red, gain of CNV; blue, loss of CNV. (B) Copy numbers of HCC driver genes in PVTTs and Ts. (C) Molecular subtype of PVTTs and Ts according to Hoshida S1-S2-S3 classes.

    To investigate aberrantly expressed genes in PVTTs and Ts,EBSeq was used to analyze the differently expressed genes(DEGs) between PVTTs, Ts, and parenchyma tumor tissue(PT) from individual patients. Gene expression profiles of these tissues (RNA-seq, GSE77509) were downloaded from the GEO database. Genes with fold change > 3 and FDR< 0.05 were identified as having significantly different expression. More than 2,000 DEGs were observed in 5 PVTTT pairs including P13, whereas less than 1,000 DEGs were observed in the other pairs (Figure 2A). The Venn diagram in Supplementary Figure S1A shows that 512 and 3,257 DEGs were identified in PVTT13 and T13 compared with PT13, respectively. Among these 1,795 DEGs, 1,006 upregulated and 789 downregulated genes were identified as genes with specifically aberrant expression in PVTT13, and 1,540 DEGs with 595 upregulated and 945 downregulated genes were identified as genes with specifically aberrant expression in T13 (Supplementary Figure S1B). GO analysis showed that upregulated DEGs in PVTT13 were enriched in 15 biological processes (BPs) including the glutamate receptor signaling pathway, cell adhesion, and cell cycle(Figure 2B, Supplementary Table S4). Downregulated DEGs in PVTT13 were enriched in 16 BPs, including positive regulation of immune response, cytokine production, and peptidyl-tyrosine phosphorylation (Figure 2B,Supplementary Table S4). Upregulated DEGs in T13 were enriched in 15 BPs, including nuclear division, cell cycle checkpoints, and regulation of protein kinase activity (Figure 2C, Supplementary Table S4). Downregulated DEGs in T13 were enriched in 17 BPs, including inflammatory response,cell-cell adhesion, and negative regulation of cell proliferation (Figure 2C, Supplementary Table S4). KEGG enrichment analysis showed that upregulated DEGs in PVTT13 were enriched in pathways for nicotine addiction and type II diabetes mellitus, and The downregulated DEGs were enriched in signaling pathways, mainly PI3K-Akt, TGF-β, and p53 signaling (Figure 2D, Supplementary Table S4).The upregulated DEGs in T13 were enriched in signaling pathways, such as IL-17, p53 and cell cycle signaling, and the downregulated DEGs in T13 were enriched in signaling pathways including cGMP-PKG, Hippo and Rap1 (Figure 2E, Supplementary Table S4). The GO and KEGG analyses showed that DEGs in PVTT13 and T13 were enriched in different BPs and pathways. These data suggest that PVTT13 and T13 were different tumor types and that different mechanisms were involved in PVTT13 and T13 development.

    Enrichment of hepatocyte gene signature in PVTT13

    To further investigate the independent origin of PVTT13,enrichment of the gene signatures of the main cell types in liver tissues were analyzed using ssGSEA. The gene signatures of bile duct cells (BDCs), hepatocytes, hepatic progenitor cells (HPCs), and hepatic stellate cells (HSCs) were included due to their potential role in liver cancer24.

    The hepatocyte gene signature showed higher enrichment value in PVTT13 (7223.62) than that in T13 (5985.81), and the HPC gene signature showed lower enrichment value in PVTT13 (2608.81) than that in T13 (3212.07) (Figure 3A).Differences in the BDC and HSC gene signature enrichment values were not observed between PVTT13 and T13 (Figure 3A), whereas other PVTT-T pairs showed similar gene signature enrichment (Figure 3A). Further gene expression analysis showed a higher expression of hepatocyte-specific genes and lower expression of HSC-specific genes in PVTT13 than in T13 (Figure 3B). In addition, previous reports suggest that cancer stem cells (CSCs) and hematopoietic stem cells may be involved in PVTT development8,12. Our findings show that gene markers of these cells, such as NANOG and C-kit, were the nodes of the gene network or aberrantly expressed in PVTT13, further consolidating their potential roles in PVTT development (Supplementary Figure S2).Therefore, enrichment of gene signatures of hematopoietic stem cells and liver cancer stem cells were also investigated. A higher enrichment value of hematopoietic stem cell gene signatures was observed in PVTT13 (3114.07) than in T13(1031.52), whereas a lower enrichment value of the EpCAM+liver cancer cell gene signature was observed in PVTT13(-315.38) than in T13 (1401.61) (Figure 3C). Most PVTTs and Ts, except for PVTT6, PVTT11, PVTT19, T22, and pairs of P10, P16, and P25, had no enrichment of these signatures(Figure 3C). Gene expression analysis showed a higher expression of hematopoietic stem cell-specific genes, such as CD34 and KIT, and a lower expression of EpCAM+liver cancer cell-specific genes, such as EpCAM and CD44 in PVTT13 (Figure 3D). These data suggested that PVTT13 expressed high levels of hepatocyte- and hematopoietic stem cell-specific genes, and that these two cell types may be involved in PVTT13 development.

    Gene signatures identifying PVTTs with no clonal relatedness with Ts

    Figure 2 Analysis of aberrantly expressed genes in PVTT13 and T13. (A) Number of differently expressed genes between PVTTs and corresponding primary Ts. Top enriched biological processes in PVTT13 (B) and T13 (C) from GO enrichment analysis based on aberrantly expressed genes. Top enriched KEGG pathways in PVTT13 (D) and T13 (E) from KEGG pathway enrichment analysis based on aberrantly expressed genes.

    Having found that some PVTTs had no clonal relatedness with the corresponding Ts and that they have unique gene expression profiles, gene signatures identifying whether PVTTs have clonal relatedness with Ts may be very helpful for PVTT treatment. An integration analysis of CNV and gene expression profiles showed that 160 genes have differentially expressed copy number variations between PVTT13 and T13 (Supplementary Figure S3). The following PPI analysis found 24 genes to be the nodes in the network,suggesting that they may have crucial roles in PVTT13 development (Supplementary Figure S3). The ability of the 24-gene signature to identify PVTTs and Ts with no clonal relatedness was validated using three publicly available datasets with unsupervised clustering. In the dataset GSE77509, which includes 19 pairs of PVTT and Ts, the 24-gene signature classified PVTT13 and T13 into two different groups and classified other PVTTs and Ts from the same patient into the same group (Figure 4A). In the dataset GSE69164, including 11 pairs of PVTTs and Ts, PVTT9 and T9 were classified into two different groups, whereas other PVTTs and Ts from an individual patient were classified into the same group (Figure 4B). In the dataset GSE74656, which includes 5 pairs of PVTTs and Ts, all the PVTTs and Ts from an individual patient were classified into the same group(Figure 4C). In addition, the ability of the 24-gene signature to identify clonal relatedness was also validated using newly collected samples (n = 8). PVTT and T from patient 2 were,however, classified into two different groups (Figure 4D). All these data indicated that the 24-gene signature may be used to identify PVTTs and Ts with no clonal relatedness.

    Figure 3 Enrichment analysis of gene signatures in PVTTs and Ts. (A) Enrichment of gene signatures of bile duct cells (BDC), hepatocytes,hepatic progenitor cells (HPC), and hepatic stellate cells (HSC) in PVTTs and Ts by ssGSEA based on their gene expression profiles. (B)Expression of the genes included in gene signatures of different liver cells in PVTTs and Ts from the dataset GSE77509. (C) Enrichment of gene signatures of hematopoietic stem cells and cancer stem cells in PVTTs and Ts by ssGSEA. (D) Expression of the genes included in gene signatures of stem cells in PVTTs and Ts from the dataset GSE77509.

    Discussion

    Figure 4 Validation of identification ability of 24-gene signature cluster analysis of pairs of PVTTs and Ts from GSE77509 (A), GSE69164 (B),GSE74656 (C), and newly collected samples (D) based on expression of 24 genes assessed by RNA-seq, cDNA microarray, or real time PCR.Orange line indicates the two subgroups classified based on the 24-gene signature.

    In the present study, we found that PVTTs may have a different clonal origin than that of the corresponding HCC Ts. These types of PVTTs have different gene expression profiles from those of the corresponding Ts, whereas PVTTs originating from primary Ts have similar gene expression profiles to those of their origin Ts tissues. Previous studies reported the involvement of several factors, such as hypoxia and non-coding RNA in PVTT development, and similar gene expression profiles of PVTTs with Ts7-10, suggesting that PVTT may be a special type of HCC intrahepatic metastasis.Here, the clonality analysis, which was established to identify paired Ts with clonal origin, was applied to investigate the clonal relationship between PVTTs and Ts. Based on the clonality analysis and 24-gene signature identification, we found that most of the PVTTs (18 out of 19 in GSE77509, 10 out of 11 in GSE69164, and all 5 of GSE74656) possibly originate from the paired Ts (Table 1, Figure 4).Additionally, we also found another type of PVTTs, which received a very low LR2 value in the clonality analysis and were classified into a different group with the paired Ts in the cluster analysis (Table 1, Figure 4). These types of PVTTs were considered to represent tumor tissues having a different clonal origin than that of the paired Ts. These findings suggest the existence of two types of PVTTs according to their clonal origin, wherein one type has the same origin as the paired Ts and the other does not. Moreover, our findings can explain why PVTTs were observed in some patients without liver parenchyma tumor nodules13. The results can also explain why PVTTs from different patients showed very different gene expression profiles11. Of note is the observation that different CNV patterns were observed between PVTT13 and T13. These CNVs may be caused by gene translocation or genomic instability in the tumor cells.Genes transfer onto other chromosomes or extrachromosomal DNA (ecDNA) is amplified in tumor cells25. Genomic instability is a well-known characteristic of tumor cells, inducing alterations in the genomes26. Hepatitis virus infection is another important factor that causes genomic instability in HCC27. Hence, the mechanism underlying the different CNV patterns between PVTTs and Ts needs more investigation.

    Our study suggests that integrated analysis of the published datasets provides a better understanding of Ts, and that such analyses may lead to improved therapeutic strategies. With the development of high-throughput technology and multi-omics analyses, increasingly more omics data from Ts are uploaded into the public databases such as GEO and The Cancer Genome Atlas (TCGA)28. An integrated analysis of these datasets makes it possible to obtain a more comprehensive understanding of Ts at multidimensional levels. The comprehensive integrated analysis of HCC using six distinct data platforms determined subtypes of HCC associated with poorer prognosis, and the gene expression signatures correlated with poor survival and potential therapeutic targets for HCC29. Another multiplatform integrative analysis based on TCGA found that molecular signatures provide a higher accuracy for the clinical outcome prediction than that of the currently used tissue-of-origin-based classification in some cases30. By integrated analysis of CNV and gene expression profiles previously used to investigate aberrantly expressed lncRNAs11, we analyzed the origin of PVTTs and developed a gene signature to identify the type of PVTTs. It has been reported that the “cell-of-origin”-based features mainly dominate the molecular taxonomy of tumor types30.Therefore, the gene signatures developed in our study may provide a platform for a high-accuracy identification of PVTTs with independent origins. However, follow-up studies and additional samples are necessary to validate the findings reported here.

    Acknowledgements

    This work was supported by grants of China National Funds for Distinguished Young Scientists (Grant No. 81425019),National Natural Science Foundation of China (Grant No.81672899), the State Key Program of National Natural Science Foundation of China (Grant No. 81730076), Shanghai Science and Technology Committee Program (Grant No.18XD1405300) and Specially-Appointed Professor Fund of Shanghai (Grant No. GZ2015009).

    Conflict of interest statement

    No potential conflicts of interest are disclosed.

    亚洲av福利一区| 亚洲成色77777| 亚洲精华国产精华液的使用体验| 两个人免费观看高清视频| 久久人人97超碰香蕉20202| 欧美成人午夜精品| 欧美日韩国产mv在线观看视频| 亚洲成国产人片在线观看| 精品少妇黑人巨大在线播放| 在线观看一区二区三区激情| 成人黄色视频免费在线看| 国产成人精品在线电影| 亚洲成色77777| 国产成人精品久久二区二区91 | 欧美bdsm另类| 久久精品久久久久久噜噜老黄| 1024视频免费在线观看| 老汉色av国产亚洲站长工具| 欧美人与善性xxx| 国产成人一区二区在线| 看十八女毛片水多多多| 亚洲熟女精品中文字幕| 一二三四中文在线观看免费高清| 自线自在国产av| 亚洲精品美女久久久久99蜜臀 | av.在线天堂| 亚洲国产av影院在线观看| 亚洲精品av麻豆狂野| 建设人人有责人人尽责人人享有的| 久久狼人影院| 天堂中文最新版在线下载| 美女国产视频在线观看| 国产在视频线精品| 欧美在线黄色| 最近最新中文字幕大全免费视频 | 婷婷色综合大香蕉| 黄色毛片三级朝国网站| 亚洲精品第二区| 啦啦啦视频在线资源免费观看| 亚洲国产看品久久| 国产乱来视频区| 亚洲三区欧美一区| 少妇的逼水好多| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲经典国产精华液单| 色吧在线观看| av电影中文网址| av线在线观看网站| 欧美日韩一区二区视频在线观看视频在线| 亚洲成国产人片在线观看| 午夜福利影视在线免费观看| 日韩电影二区| 成人毛片60女人毛片免费| 日韩三级伦理在线观看| www日本在线高清视频| 精品亚洲成a人片在线观看| www.av在线官网国产| 午夜免费观看性视频| 欧美中文综合在线视频| 99久久中文字幕三级久久日本| av国产精品久久久久影院| 亚洲av成人精品一二三区| 成人国语在线视频| 视频区图区小说| 国产成人精品久久久久久| 日本免费在线观看一区| 国产精品一区二区在线观看99| 男女下面插进去视频免费观看| 亚洲美女黄色视频免费看| 国产精品 欧美亚洲| 人妻少妇偷人精品九色| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 色吧在线观看| 国产精品熟女久久久久浪| 午夜影院在线不卡| 2022亚洲国产成人精品| 美女xxoo啪啪120秒动态图| 亚洲精品中文字幕在线视频| 色94色欧美一区二区| 欧美亚洲 丝袜 人妻 在线| 国产男人的电影天堂91| av国产精品久久久久影院| 欧美精品亚洲一区二区| 男女无遮挡免费网站观看| 国产女主播在线喷水免费视频网站| 色哟哟·www| 国产成人aa在线观看| 国产不卡av网站在线观看| 在线观看免费日韩欧美大片| 国产成人免费无遮挡视频| 天天躁夜夜躁狠狠躁躁| 欧美日韩一区二区视频在线观看视频在线| 多毛熟女@视频| 亚洲五月色婷婷综合| 国产又色又爽无遮挡免| 欧美bdsm另类| 国语对白做爰xxxⅹ性视频网站| 国产亚洲精品第一综合不卡| 伊人亚洲综合成人网| 国产精品 欧美亚洲| 日韩一区二区三区影片| 女人久久www免费人成看片| 看非洲黑人一级黄片| 9191精品国产免费久久| 久久精品熟女亚洲av麻豆精品| 亚洲激情五月婷婷啪啪| 啦啦啦中文免费视频观看日本| 在线亚洲精品国产二区图片欧美| 欧美日韩一区二区视频在线观看视频在线| 中文字幕色久视频| 欧美日韩亚洲国产一区二区在线观看 | 精品国产乱码久久久久久小说| 国产精品女同一区二区软件| 人妻系列 视频| 国产午夜精品一二区理论片| 美女高潮到喷水免费观看| 亚洲精品,欧美精品| 男女国产视频网站| 丰满迷人的少妇在线观看| 午夜激情av网站| 久久精品国产a三级三级三级| 免费少妇av软件| 久久精品国产亚洲av涩爱| 男女下面插进去视频免费观看| 亚洲美女视频黄频| 日本av手机在线免费观看| av女优亚洲男人天堂| 丁香六月天网| 女人久久www免费人成看片| 黄片无遮挡物在线观看| 国产精品免费大片| 亚洲综合精品二区| 亚洲三区欧美一区| 18禁国产床啪视频网站| 永久网站在线| 精品一区二区免费观看| 久久人人爽av亚洲精品天堂| 亚洲成人一二三区av| 捣出白浆h1v1| 精品人妻偷拍中文字幕| 十分钟在线观看高清视频www| 日日爽夜夜爽网站| 啦啦啦啦在线视频资源| 午夜福利,免费看| 日产精品乱码卡一卡2卡三| 9色porny在线观看| 一边摸一边做爽爽视频免费| 国产精品久久久av美女十八| 精品国产露脸久久av麻豆| 精品久久久久久电影网| 在线观看国产h片| 99久久综合免费| 亚洲精品在线美女| 欧美变态另类bdsm刘玥| 成年美女黄网站色视频大全免费| 精品一区二区免费观看| 看免费av毛片| 伊人久久大香线蕉亚洲五| 日韩,欧美,国产一区二区三区| 免费观看a级毛片全部| 亚洲,一卡二卡三卡| 天天躁日日躁夜夜躁夜夜| 飞空精品影院首页| 18禁国产床啪视频网站| 肉色欧美久久久久久久蜜桃| 桃花免费在线播放| 国产视频首页在线观看| 韩国av在线不卡| 宅男免费午夜| 欧美日韩视频精品一区| 精品第一国产精品| 精品少妇一区二区三区视频日本电影 | 美女脱内裤让男人舔精品视频| 欧美bdsm另类| 色播在线永久视频| 亚洲美女搞黄在线观看| 亚洲中文av在线| 一区二区av电影网| 9191精品国产免费久久| 国产亚洲精品第一综合不卡| 宅男免费午夜| 一区在线观看完整版| 一区二区三区四区激情视频| 国产无遮挡羞羞视频在线观看| 国产乱人偷精品视频| 久久国产精品大桥未久av| 日韩精品有码人妻一区| 久久综合国产亚洲精品| 久久精品久久久久久久性| 国产日韩欧美在线精品| 人妻一区二区av| 视频在线观看一区二区三区| 男女高潮啪啪啪动态图| 一级黄片播放器| 秋霞伦理黄片| 久久99一区二区三区| 亚洲av成人精品一二三区| 亚洲 欧美一区二区三区| 亚洲精品中文字幕在线视频| 精品人妻熟女毛片av久久网站| 成人手机av| 亚洲精品国产色婷婷电影| 亚洲欧美日韩另类电影网站| 赤兔流量卡办理| freevideosex欧美| 国产人伦9x9x在线观看 | 亚洲第一区二区三区不卡| 亚洲成国产人片在线观看| 亚洲av免费高清在线观看| 一级爰片在线观看| 久久99一区二区三区| 午夜久久久在线观看| 侵犯人妻中文字幕一二三四区| 国产精品国产av在线观看| 国产黄频视频在线观看| 亚洲精品国产色婷婷电影| 欧美人与性动交α欧美精品济南到 | 久久国产精品男人的天堂亚洲| 国产精品国产av在线观看| 三级国产精品片| 丝袜喷水一区| 国产一区二区三区综合在线观看| 两性夫妻黄色片| 欧美中文综合在线视频| 国产精品免费视频内射| av在线老鸭窝| 黄频高清免费视频| 日本vs欧美在线观看视频| 国产精品久久久久久av不卡| 女人被躁到高潮嗷嗷叫费观| 人人妻人人爽人人添夜夜欢视频| 五月开心婷婷网| 哪个播放器可以免费观看大片| 中文字幕最新亚洲高清| 最近的中文字幕免费完整| 久久久a久久爽久久v久久| 看免费成人av毛片| 亚洲天堂av无毛| 综合色丁香网| 青春草国产在线视频| 亚洲精品美女久久久久99蜜臀 | 亚洲欧美色中文字幕在线| 日本色播在线视频| 国产乱来视频区| 制服丝袜香蕉在线| 精品久久久精品久久久| 欧美激情极品国产一区二区三区| www.熟女人妻精品国产| 免费少妇av软件| 最新中文字幕久久久久| 亚洲精品在线美女| 18在线观看网站| 性色avwww在线观看| 久久精品国产亚洲av涩爱| 免费久久久久久久精品成人欧美视频| 欧美成人午夜免费资源| 亚洲色图 男人天堂 中文字幕| 亚洲四区av| 大香蕉久久网| 国产成人av激情在线播放| 日韩视频在线欧美| 少妇的丰满在线观看| 久久久久久久国产电影| 欧美人与性动交α欧美精品济南到 | 亚洲第一av免费看| 亚洲欧美成人综合另类久久久| 亚洲欧美精品综合一区二区三区 | 亚洲 欧美一区二区三区| 在线观看一区二区三区激情| 日韩av不卡免费在线播放| 久久99精品国语久久久| 亚洲精品,欧美精品| 午夜日韩欧美国产| 国产精品国产av在线观看| 国产国语露脸激情在线看| 免费在线观看视频国产中文字幕亚洲 | 欧美 亚洲 国产 日韩一| 精品一区二区三卡| 午夜日本视频在线| 成人18禁高潮啪啪吃奶动态图| 亚洲精品中文字幕在线视频| 久久久久久久久久久免费av| 免费黄网站久久成人精品| 日韩电影二区| 国产高清不卡午夜福利| 国产精品人妻久久久影院| 水蜜桃什么品种好| 国产毛片在线视频| 少妇熟女欧美另类| 久久av网站| 免费观看a级毛片全部| 亚洲,一卡二卡三卡| 男人舔女人的私密视频| 最近中文字幕2019免费版| 免费不卡的大黄色大毛片视频在线观看| 亚洲av免费高清在线观看| 精品国产一区二区三区四区第35| 一级毛片 在线播放| 国产精品一二三区在线看| 国产成人精品无人区| 亚洲国产欧美网| 亚洲国产日韩一区二区| 我的亚洲天堂| 亚洲,欧美精品.| 国产男女内射视频| 成年动漫av网址| 桃花免费在线播放| 在线观看人妻少妇| 国产精品秋霞免费鲁丝片| 久热久热在线精品观看| 欧美人与性动交α欧美精品济南到 | 久久久久久久大尺度免费视频| 亚洲av中文av极速乱| 国产精品熟女久久久久浪| 老司机影院毛片| 精品第一国产精品| 中文字幕亚洲精品专区| 色94色欧美一区二区| xxx大片免费视频| 成人二区视频| 亚洲精品久久久久久婷婷小说| 色婷婷久久久亚洲欧美| 国产女主播在线喷水免费视频网站| 国产在视频线精品| 免费观看在线日韩| 午夜福利影视在线免费观看| 免费观看在线日韩| 老司机影院成人| 一级毛片黄色毛片免费观看视频| 又黄又粗又硬又大视频| 欧美激情高清一区二区三区 | 亚洲国产精品国产精品| 精品国产一区二区久久| 午夜精品国产一区二区电影| 看非洲黑人一级黄片| 美女视频免费永久观看网站| 少妇被粗大猛烈的视频| 国产精品蜜桃在线观看| 视频区图区小说| 亚洲国产成人一精品久久久| 看非洲黑人一级黄片| 久热久热在线精品观看| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 欧美日韩视频精品一区| 中文字幕亚洲精品专区| 久久久精品国产亚洲av高清涩受| 欧美成人午夜免费资源| 国产成人91sexporn| 国产在线一区二区三区精| 日韩不卡一区二区三区视频在线| 伊人久久国产一区二区| 国产精品亚洲av一区麻豆 | 狠狠婷婷综合久久久久久88av| 午夜福利视频精品| 欧美日韩综合久久久久久| 9191精品国产免费久久| 肉色欧美久久久久久久蜜桃| 亚洲精品日本国产第一区| 在线看a的网站| 久久国内精品自在自线图片| 日本免费在线观看一区| 亚洲一区中文字幕在线| 中文天堂在线官网| 99热全是精品| 街头女战士在线观看网站| 啦啦啦啦在线视频资源| 亚洲天堂av无毛| 久热久热在线精品观看| 亚洲av成人精品一二三区| 黑人巨大精品欧美一区二区蜜桃| 亚洲av.av天堂| 91成人精品电影| 久久这里只有精品19| 老熟女久久久| 精品国产乱码久久久久久小说| 亚洲欧美精品综合一区二区三区 | 男人添女人高潮全过程视频| 黑丝袜美女国产一区| 日韩中文字幕视频在线看片| 亚洲精品成人av观看孕妇| 国产一区二区三区综合在线观看| 亚洲精品乱久久久久久| 麻豆精品久久久久久蜜桃| 亚洲精品一二三| 国产97色在线日韩免费| 久久久精品区二区三区| 大话2 男鬼变身卡| 一区二区三区四区激情视频| 精品人妻在线不人妻| 亚洲欧洲国产日韩| 久久精品国产a三级三级三级| 成人黄色视频免费在线看| 飞空精品影院首页| 精品久久久久久电影网| 国产免费又黄又爽又色| 精品一区在线观看国产| 亚洲精品美女久久久久99蜜臀 | 少妇人妻 视频| 中文字幕亚洲精品专区| av免费在线看不卡| 国产色婷婷99| 男女国产视频网站| 久久av网站| 久久99热这里只频精品6学生| 欧美日本中文国产一区发布| 校园人妻丝袜中文字幕| 国产男女超爽视频在线观看| 成年人午夜在线观看视频| 天天操日日干夜夜撸| 97精品久久久久久久久久精品| 欧美日韩成人在线一区二区| 老汉色∧v一级毛片| av免费观看日本| 午夜精品国产一区二区电影| 中国三级夫妇交换| 80岁老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 午夜福利视频在线观看免费| 纵有疾风起免费观看全集完整版| 制服丝袜香蕉在线| 一级毛片电影观看| 久久这里有精品视频免费| 午夜久久久在线观看| 亚洲伊人色综图| 日韩av在线免费看完整版不卡| 男人爽女人下面视频在线观看| 狠狠婷婷综合久久久久久88av| 91精品伊人久久大香线蕉| 中文字幕人妻丝袜一区二区 | 欧美日韩亚洲国产一区二区在线观看 | 韩国av在线不卡| 亚洲一码二码三码区别大吗| 热re99久久国产66热| 捣出白浆h1v1| 欧美精品av麻豆av| 一级爰片在线观看| 哪个播放器可以免费观看大片| 人妻系列 视频| 一级片免费观看大全| 最近手机中文字幕大全| 欧美激情 高清一区二区三区| 欧美变态另类bdsm刘玥| 国产精品秋霞免费鲁丝片| 美女福利国产在线| 国产精品久久久久成人av| 最近中文字幕2019免费版| 王馨瑶露胸无遮挡在线观看| 久久久久久伊人网av| 国产在线一区二区三区精| 成人手机av| 亚洲国产色片| 制服丝袜香蕉在线| 99久国产av精品国产电影| 黄频高清免费视频| 午夜日韩欧美国产| 尾随美女入室| 久久久久国产一级毛片高清牌| 久久久久久久国产电影| 亚洲美女黄色视频免费看| 精品少妇一区二区三区视频日本电影 | 久久久久人妻精品一区果冻| 日韩精品免费视频一区二区三区| 中国国产av一级| 日本黄色日本黄色录像| 人妻人人澡人人爽人人| 天天操日日干夜夜撸| 国产亚洲一区二区精品| 亚洲国产av影院在线观看| 日韩视频在线欧美| 两个人看的免费小视频| 大香蕉久久网| 老司机亚洲免费影院| 汤姆久久久久久久影院中文字幕| 天天躁夜夜躁狠狠躁躁| 777米奇影视久久| a级片在线免费高清观看视频| 看免费av毛片| 免费女性裸体啪啪无遮挡网站| 777久久人妻少妇嫩草av网站| 如日韩欧美国产精品一区二区三区| 亚洲美女搞黄在线观看| 熟女少妇亚洲综合色aaa.| 亚洲第一青青草原| 美女福利国产在线| 国产女主播在线喷水免费视频网站| 天堂俺去俺来也www色官网| 男人操女人黄网站| 一个人免费看片子| 制服人妻中文乱码| av有码第一页| 国产精品熟女久久久久浪| 久久99精品国语久久久| 日日啪夜夜爽| 欧美中文综合在线视频| 久久鲁丝午夜福利片| 国语对白做爰xxxⅹ性视频网站| 一级,二级,三级黄色视频| 老汉色∧v一级毛片| 中国三级夫妇交换| a级片在线免费高清观看视频| 亚洲成人手机| 久久久亚洲精品成人影院| 久久久久久人人人人人| 日韩三级伦理在线观看| 亚洲,欧美,日韩| videossex国产| 美女大奶头黄色视频| 欧美亚洲日本最大视频资源| 国产一区亚洲一区在线观看| 国产又色又爽无遮挡免| 精品卡一卡二卡四卡免费| 人成视频在线观看免费观看| 亚洲国产成人一精品久久久| 久久久国产一区二区| 不卡视频在线观看欧美| 女性被躁到高潮视频| 波多野结衣一区麻豆| 中文字幕亚洲精品专区| 视频区图区小说| 久久精品国产鲁丝片午夜精品| 午夜福利乱码中文字幕| 久久久久精品久久久久真实原创| 美女视频免费永久观看网站| 国产精品成人在线| 丰满饥渴人妻一区二区三| 亚洲精品久久久久久婷婷小说| 麻豆乱淫一区二区| 视频在线观看一区二区三区| 午夜91福利影院| 宅男免费午夜| 人人澡人人妻人| 久久毛片免费看一区二区三区| 少妇人妻 视频| 亚洲精品一区蜜桃| 久久久精品94久久精品| 老司机影院毛片| 久久人人爽人人片av| 99热网站在线观看| 国产精品不卡视频一区二区| 美国免费a级毛片| 免费人妻精品一区二区三区视频| 97在线人人人人妻| 亚洲情色 制服丝袜| 99国产综合亚洲精品| 建设人人有责人人尽责人人享有的| 一区二区av电影网| 在线观看免费视频网站a站| av天堂久久9| 免费高清在线观看日韩| 97人妻天天添夜夜摸| 国产一区二区三区av在线| 国产精品麻豆人妻色哟哟久久| 精品第一国产精品| 丝袜美腿诱惑在线| 久久精品夜色国产| 久久久久精品久久久久真实原创| 成年女人在线观看亚洲视频| 亚洲精品av麻豆狂野| 丝袜美腿诱惑在线| 天堂俺去俺来也www色官网| 日本vs欧美在线观看视频| 精品人妻在线不人妻| 91精品国产国语对白视频| 亚洲一码二码三码区别大吗| 女人被躁到高潮嗷嗷叫费观| 自拍欧美九色日韩亚洲蝌蚪91| 在线天堂中文资源库| 不卡av一区二区三区| 午夜av观看不卡| 日本91视频免费播放| 中国国产av一级| 国产av码专区亚洲av| 国产亚洲一区二区精品| 成人国语在线视频| 伦精品一区二区三区| 超碰成人久久| 黄片无遮挡物在线观看| 午夜福利视频在线观看免费| 在线观看美女被高潮喷水网站| 人体艺术视频欧美日本| 亚洲精品成人av观看孕妇| 亚洲av男天堂| 人妻一区二区av| 亚洲精品国产av成人精品| 两性夫妻黄色片| 日本欧美视频一区| 久久97久久精品| 电影成人av| 亚洲伊人色综图| 男女国产视频网站| 亚洲经典国产精华液单| 狂野欧美激情性bbbbbb| 亚洲综合色惰| 国产成人一区二区在线| 国产精品嫩草影院av在线观看| 美女主播在线视频| 青春草视频在线免费观看| 久久国产亚洲av麻豆专区| 欧美最新免费一区二区三区| 欧美人与性动交α欧美软件| 欧美成人午夜精品| 1024视频免费在线观看| 飞空精品影院首页| 亚洲美女视频黄频| 国产日韩欧美亚洲二区| 日韩中文字幕视频在线看片| 久久婷婷青草| 国产一区二区 视频在线| 亚洲在久久综合| 美女福利国产在线| 在线天堂最新版资源| 国产老妇伦熟女老妇高清| 亚洲激情五月婷婷啪啪| 丰满饥渴人妻一区二区三|