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

    Screening and comprehensive analysis of key genes in liver hepatocellular carcinoma based on bioinformatics

    2022-02-04 10:32:00WANGRuiBAIHaotianLIYalanYANGJing
    Journal of Hainan Medical College 2022年21期

    WANG Rui,BAI Hao-tian,LI Ya-lan,YANG Jing

    1.Heilongjiang University of Traditional Chinese Medicine, Harbin 150040, China

    Keywords:

    ABSTRACT

    1. Introduction

    Hepatocellular carcinoma belongs to liver cancer. It is one of the malignant tumors of digestive system that pose a great threat to human life and health all over the world. In recent years, the global cancer statistics show that liver cancer ranks the top of the list of malignant tumors, from lung cancer, breast cancer and stomach cancer, from the perspective of incidence rate or mortality rate[1]. The incidence rate of incidence of cancer increased year by year[2], the incidence rate of liver cancer in 2007 was about 700 thousand. By 2017, the incidence rate of liver cancer had reached 950 thousand. The incidence of female was generally higher than that of men, and the mortality rate of men was about 2 times[3].Worldwidely, more than half of the incidence rate is from China,which shows that liver cancer is a serious threat to people's health.Liver cancer is generally divided into primary liver cancer and metastatic liver cancer[4]. There is no essential difference between primary liver cancer and hepatocellular carcinoma, which is difficult to treat, difficult to find, easy to metastasis and high mortality. Due to the hidden incidence of hepatocellular carcinoma, most of the patients found clinically are late stage patients[5]. At present, there are many clinical treatment methods for hepatocellular carcinoma,but the treatment effect is not very good and the survival rate is low.Even in developed countries, the survival rate is only about 3%[6].The widely accepted pathogenesis of hepatocellular carcinoma is a complex process involving multiple genes and affected by both heredity and environment. Among these pathogenic factors, the most important are hepatitis B virus and hepatitis C virus[7]. Therefore, the study of the pathogenesis of hepatocellular carcinoma and the related genes of hepatocellular carcinoma can promote the early diagnosis and treatment of hepatocellular carcinoma, and effectively alleviate the current situation of hepatocellular carcinoma. Existing literature reports show that the molecular mechanisms of hepatocellular carcinoma mainly include proto oncogene activation and tumor suppressor gene inactivation, abnormal activation of multiple molecular signaling pathways, abnormal expression of related proteins and growth factors in vivo[8]. The molecular mechanism of hepatocellular carcinoma is complex. At present, the most dominant signal pathway has not been found. Therefore, in recent years, more and more biological analysis methods have been applied to cancer research. Mining from gene chip through bioinformation technology can explore the pathogenesis of hepatocellular carcinoma from the molecular level.

    In this paper, we use GEO to search hepatocellular carcinoma genes, and then screen hepatocellular carcinoma related expression genes. Then, we analyze go function and KEGG function, construct differential gene characteristics and protein interaction network,screen key genes, and clarify the differential expression of key genes, prognosis analysis and immune infiltration. It provides a new idea for the selection of prognostic biomarkers of hepatocellular carcinoma.

    2. Materials and methods

    2.1 Microarray data analysis

    At GEO(https://www.ncbi.nlm.nih.gov/geo/)Input the search term"liver hepatecellular carcinoma" into the database to download the gene expression profile of gse109903. The expression profile is from Affymetrix GPL19057 platform and provided by Fong el.GSE109903 data included 28 samples, including 14 hepatocellular carcinoma samples and 14 normal control samples.

    2.2 Screening differential genes

    The downloaded sample data were grouped and named as hepatocellular carcinoma group and control group respectively.Biojupiter was introduced to process and screen the differential genes. Principal component analysis is a statistical method used to identify the global pattern of concentrated high-dimensional data. It is usually used to explore the high similarity of biological samples in RNA SEQ data sets[9]. The processed data were further screened to identify differentially expressed genes (DEGs). Set the condition log2 fold change 2, and the corrected data set condition P < 0.05 is used as the screening condition.

    2.3 Functional analysis and pathway analysis of differential genes

    Go analysis process exists in large-scale function enrichment and is a common research method in function enrichment research. Gene function can be divided into the following three steps: biological process (BP), molecular function (MF) and cellular component (CC).KEGG database is also widely used in the process of biomedical analysis, which contains a large number of data on genome,biological processes and diseases. Then, the David database is used for annotation analysis and pathway enrichment analysis of DEGs in this study. Among them, the screening conditions set in go analysis and KEGG analysis were P<0.05 and gene count>5.

    2.4 Visualization of differential gene expression

    L1000FWD database(https://maayanlab.cloud/l1000fwd/),it provides interactive visualization of gene expression characteristics induced by more than 16000 drugs and small molecules. It also supports coloring the characteristics through different attributes(such as cell type, time point, concentration) and drug attributes, so as to analyze the function of small molecules and the mode of action of drugs, so as to infer the affinity of drugs to the above different genes[10]. The obtained drug data were ranked from high to low according to the similarity score.

    2.5 Construction of PPI interaction network

    String database(https://string-db.org/)It is a website that can visualize protein interactions[11] and can be used to predict the analysis database of gene interaction network in EMBL[12]. The relationship between different genes can be visually analyzed. The protein network interaction map can be obtained by submitting the studied genes to the string database for protein interaction data analysis, setting a high confidence score (0.500) and hiding the nodes without interaction relationship. Screen the proteins related to hepatocellular carcinoma for further research.

    2.6 Differential expression of key genes and prognostic analysis

    GEPIA (http://GEPIA.cancer-pku.cn/index.html)It is an analysis tool using standard processing pathway, which is composed of thousands of tumor and normal tissue sample data[13]. In our study,GEPIA was used to compare the differential expression of genes between tumor and normal tissues, analyze the pathological stage,and finally analyze the related prognosis. It is generally considered that there is a significant difference between P<0.05.

    Kaplan-Meier-plotter(http://kmplot.com/analysis/index.php?p=service&cancer= liver _ rnaseq)The mapping tool can carry out survival analysis for more than 54000 genes (mRNA, miRNA,protein) in 21 kinds of tumors, including breast cancer, ovarian cancer, hepatocellular carcinoma and gastric cancer. The mRNA expression of core genes screened in patients with hepatocellular carcinoma was analyzed by Kaplan Meier plotter. It was statistically significant when P<0.05[14].

    2.7 Function prediction of key genes

    GeneMANIA (http://www.genemania.org)It is a website with gene information and prioritization of genes for functional analysis, with high accuracy prediction algorithm[15]. Input the screened key genes into the platform, analyze the functions of the screened key genes,and infer the predictive value of the screened genes.

    2.8 Relationship between key gene expression and immune infiltration level of hepatocellular carcinoma

    TIMER2.0(https://cistrome.shinyapps.io/TIMER/),the infiltration of different immune cells and their clinical effects can be systematically evaluated through six main analysis modules[16]. The selected key genes were input through the "gene module" to generate a scatter diagram to observe the relationship between their expression and the level of immune invasion of hepatocellular carcinoma.

    3. Results

    3.1 Data processing and differential gene screening results

    In principal component analysis, each point represents a selected sample, and different colors represent the control group and model group, thus forming an interactive three-dimensional scatter diagram.If the sample has a more similar gene expression profile, the distance will be closer in the three-dimensional space (Fig. 1). As shown in Fig 1, the expressions of the control group and the model group are independent of each other, so the next analysis can be carried out. After processing the obtained gene text, 1509 differentially expressed genes were screened, including 637 up-regulated genes and 872 down-regulated genes. Finally, the screened results were visually analyzed and volcanic map and cluster map were drawn(Fig. 2).

    Fig1 Principal component analysis results

    3.2 Functional enrichment of differential gene go and analysis results of KEGG pathway

    After functional enrichment analysis and pathway analysis, the results are visualized (Fig. 3). After analyzing 1509 samples, the following results were obtained. Go enrichment analysis mainly includes three components. Firstly, the biological process (BP)enriched by differential genes includes positive regulation of transcription (DNA template) with DNA as template and positive regulation of transcription from RNA polymerase II promoter,Regulation of transcription (DNA templated), nucleosome assembly,chromatin assembly, nucleosome organization, mitochondrial ATP synthesis electron transport, Protein DNA complex assembly,respiratory electron transport chain, etc; The component analysis(CC) of differential genes mainly focused on nuclear body, nuclear chromatin, chromatin, focal adhesionx, nuclear chromatome part,nuclear speck, nuclear plast part, mitochondrial inner membrane,Mitochondrial respiratory chain complex I, nuclear periphery, etc;The molecular function (MF) of differential genes was analyzed as RNA binding, transcription coactivator activity, cadherin binding,RNA polymerase II core promoter proximal region sequence specific DNA binding, Protein kinase binding, RNA polymerase II transcription cofactor activity, activating transcription factor binding;Secondly, the pathways of differentially expressed genes were analyzed. The results are as follows: the main enrichment pathways of differentially expressed genes are related to systemic lupus erythematosus, alcoholism and other diseases, RNA polymerase I promoter opening and DNA methylation, Activated pkn1 stimulates the transcription of AR regulated genes klk2 and Klk3. The enriched genes and joint scores of each pathway are shown in the table(Tab.1).

    Fig2 Volcano map and heat map of differential gene clustering

    Fig3 GO enrichment analysis results of differential genes

    3.3 Interactive visualization results of differential gene characteristic expression

    L1000FWD platform can visualize the expression characteristics of more than 16000 drugs and small molecule induced genes.After introducing the differential gene into l1000FWD platform,the similarity results of small molecule and drug characteristics are obtained, as shown in Tab.2. Drugs with similar characteristicsof differential gene mainly include CDC inhibitor, prostaglandin,serotonin receptor antagonist, BAF transcription repressor, tyrosine phosphatase inhibitor, etc. After the data is processed, the shape represents the time point, the color represents the type of cells, and the fireworks diagram is drawn (Fig. 4).

    Tab1 KEGG pathway analysis

    3.4 Visualization results of differential gene PPI network

    out by using string database to explore the potential interaction between them. As expected, 160 nodes and 200 connections are obtained in the PPI network diagram. After removing the nodes without interaction relationship, the string connection diagram is drawn (Fig. 5). As shown in the Fig.5, the genes closely related to the occurrence and development of hepatocellular carcinoma are mainly EEF1A1, HK2, FAM38A, LAMB3, etc.

    PPI network analysis of differentially expressed genes was carried

    Fig4 Visual pyrotechnic map of the expression of differentially expressed genes

    Tab2 Main drugs with similar characteristics with different genes

    Fig5 Network diagram of different gene protein interactions

    3.5 Different expression of EEF1A1and HK2 in patients with hepatocellular carcinoma

    It can be seen from Fig. 5 that EEF1A1 is connected with 14 related proteins, with a total of 26 lines; HK2 was linked to 13 related proteins, with a total of 24 connections; Fam83a was connected with 8 related proteins, with a total of 22 lines; LAMB3 is connected with 5 related proteins, with a total of 18 lines. By comparing the degree of its internal network, we further analyzed the selected genes EEF1A1 and HK2, which are closely related to the occurrence and development of hepatocellular carcinoma.

    Using the GEPIA data set, we compared the mRNA expression of EEF1A1 and HK2 in LIHC and normal hepatocytes. The results showed that the expression levels of EEF1A1and HK2 in hepatocellular carcinoma were higher than those in normal tissues(Fig.6). Next, we continued to evaluate the correlation between the expression difference of EEF1A1and HK2 and the pathological stage of patients with hepatocellular carcinoma. According to the violin diagram (Fig.7), the P values of EEF1A1and HK2 groups are less than 0.05, and the results are meaningful. The central points of EEF1A1and HK2 are not on the same straight line, and they all change significantly in the fourth stage, which proves that genes play an important role in this stage. The diagram shows that EEF1A1and HK2 may play an important role in the occurrence and development of hepatocellular carcinoma.

    Fig6 Expression of EEF1A1 and HK2 in patients with liver hepatocellular carcinoma

    3.6 Prognostic value of EEF1A1and HK2 expression in patients with hepatocellular carcinoma

    Fig7 Correlation between EEF1A1 and HK2 expression and tumor stage in patients with liver hepatocellular carcinoma

    In order to evaluate the value of the expression of EEF1A1 and HK2 in the process of colorectal cancer, GEPIA was used to analyze the correlation between EEF1A1, HK2 and clinical results.(Fig. 8)shows the curve of disease-free survival (DFS) and overall survival (OS). Patients with high transcription level of EEF1A1(P= 0.036) were significantly correlated with short-term disease-free survival. Similarly, patients with high transcription level of HK2 (P= 0.02) were also significantly correlated with short-term diseasefree survival. It is proved that both EEF1A1 and HK2 are related to survival.

    In addition, we also performed subgroup analysis of the model using the Kaplan Meier plotter mapping tool to obtain the prognostic value of EEF1A1 and HK2 in patients with hepatocellular carcinoma (Fig. 9). The expression of EEF1A1(HR = 0.8, P =0.03) and HK2 (HR = 1.82, P = 0.00086) mRNA in patients with hepatocellular carcinoma was significantly correlated with shortterm overall survival, which was the same as that obtained by GEPIA.

    3.7 Prediction results of differential gene function analysis

    Fig8 Prognostic value of EEF1A1 and HK2 mRNA expression in liver hepatocellular carcinoma (GEPIA)

    Fig9 Prognostic value of EEF1A1 and HK2 in liver hepatocellular carcinoma (Kaplan-Meier plotter)

    EEF1A1and HK2 genes were input into the GeneMIANA platform respectively, and the interaction network diagram was obtained after processing (Fig.10). It can be seen from the Fig.10 that the results reflect the functions of differentially expressed genes EEF1A1and HK2 and the functions of genes similar to these two genes. For example, CHRM4, EEF1A2, ICA1 and GATA1 are also reflected in the string network interaction diagram. The Fig.10 shows the functionally similar parts of the selected genes and their similar genes. It is mainly related to translation factor activity, chaperone mediated autophagy, translation regulator activity, protein targeted endoplasmic reticulum and carbohydrate catabolism, purine nucleoside diphosphate metabolism, ADP metabolism, etc. Physical interaction accounted for 77.64%, and intergenic interaction accounted for 2.87%.

    Fig10 Cluster analysis of key genes EEF1A1, HK2 and their similar genes

    3.8 Immune cell infiltration of EEF1A1and HK2 in patients with liver cancer

    The level of immune cells is related to the proliferation and development of cancer cells. In this study, we used TIMER2.0 database to explore the correlation between EEF1A1and HK2 immune cell infiltration. Input EEF1A1and HK2 into TIMER2.0 platform, and the results are shown in Fig.11. The expression of EEF1A1was positively correlated with the infiltration of B cells,CD8 + T cells, CD4 + T cells, macrophages and neutrophils in hepatocellular carcinoma. The expression of HK2 was positively correlated with the infiltration of macrophages and neutrophils in hepatocellular carcinoma, but not with the infiltration of B cells,CD8 + T cells, CD4 cells and dendritic cells .

    Fig11 Correlation between differentially expressed genes EEF1A1 and HK2 and immune cell infiltration

    4. Discussion

    As one of the malignant tumors that cause great harm to global health, hepatocellular carcinoma has shown an upward trend year by year. Due to its high degree of malignancy and easy recurrence,traditional treatment methods are difficult to improve the survival rate of patients with hepatocellular carcinoma, and due to the influence of genetic factors, the treatment methods of hepatocellular carcinoma have attracted extensive attention[17-18], Therefore, it is urgent to study the pathogenesis and related genes of hepatocellular carcinoma. With the development of bioinformatics, it provides great possibility and great convenience for studying the pathogenesis of hepatocellular carcinoma, exploring new targets, inventing new treatment methods and improving the survival rate of patients. In this study, using the rich data information of geo database, a total of 1509 differentially expressed genes were screened by bioinformatics technology, and the number of up-regulated and down-regulated genes were 637 and 872 respectively. PPI network visual analysis screened four genes closely related to hepatocellular carcinoma,including EEF1A1, HK2, fAM38A and LAMB3. Through node relationship comparison, it is clear that the core genes significantly related to hepatocellular carcinoma are EEF1A1and HK2. According to GEPIA and Kaplan-Meier plotter databases, the screened genes are closely related to the pathological stage and prognosis of hepatocellular carcinoma. GeneMANIA database shows that key genes are mainly related to translation factor activity, chaperone mediated autophagy, translation regulator activity, protein targeted endoplasmic reticulum and carbohydrate catabolism, purine nucleoside diphosphate metabolism, ADP metabolism, etc. Finally,through the analysis of TIMER2.0 database, it was found that the two genes were positively correlated with different immune cell infiltration in cancerous tissues. The differential expression,prognostic value and immune cell infiltration of EEF1A1 and HK2 in patients with hepatocellular carcinoma were further clarified. The results showed that EEF1A1 and HK2 mRNA were significantly expressed in HCC cells than in normal tissues, which may be potential prognostic biomarkers for the survival of HCC patients.

    EEF1A1 is a trigger for protein translation extension[19-20],.which can participate in cell biological behavior, express significantly in a variety of cancer cells, and significantly inhibit the proliferation,invasion and migration of tumor cells. Zhou Shuting et al.[21]studied the effect of EEF1A1 on HCC HCA-P cells through RNA interference technology related experiments, and further analyzed the mechanism involved in regulation. The results showed that shRNA targeting EEF1A1could reduce the expression of EEF1A1, and low expression of EEF1A1could induce apoptosis and inhibit proliferation and migration of HCC cells. The regulation mechanism showed that EEF1A1 was closely related to ANXA7, It can participate in the regulation of tumorigenesis. The results also showed that the expression of EEF1A2 decreased with the decrease of EEF1A1,and the two were highly correlated. Qiu FN et al.[22] detected the expression of EEF1A2 in 62 HCC tissue samples and hepatocellular carcinoma by real-time PCR and immunohistochemistry. The results showed that EEF1A2 could effectively silence the proliferation,migration and invasion of cancer cells, promote the increase of apoptosis of liver cancer cells and induce cell cycle arrest.EEF1A2 reduces the occurrence of hepatocellular carcinoma by inhibiting P13K/Akt/ NFKB signaling pathway. HK2 is the first key rate limiting enzyme in the glycolysis pathway of the body. Its combination with mitochondrial outer membrane voltage dependent cation channel protein can promote the occurrence of glycolysis and the formation of tumor cells[23]. Inhibiting the expression of HK2 can inhibit the increase of tumor in tumor cells and animal models.Zhang Yu et al.[24] detected the expression of HK2 in 48 liver cancer tissues by immunohistochemistry, A series of experiments such as statistical analysis explored the relationship between HK2 and hepatocellular carcinoma. The results showed that HK2 was highly expressed in hepatocellular carcinoma, which could significantly reduce the value-added activity of hepatoma cells, make the cell cycle stagnate in S phase and promote the apoptosis of hepatoma cells. Finally, it was proved that HK2 showed strong antitumor effect and had certain clinical significance, It also lays a foundation for our further study on the survival relationship between HKs family and hepatocellular carcinoma.

    In conclusion, this experiment uses bioinformatics and gene chip to explore the differentially expressed genes of hepatocellular carcinoma. After preliminary analysis of the genes, we get the signal pathways and functions involved in the genes, and further screen the key genes EEF1A1 and HK2, so as to provide potential prognostic biomarkers for the survival of patients with hepatocellular carcinoma, To provide theoretical support, new ideas and schemes for the clinical treatment and mechanism research of hepatocellular carcinoma, but more rigorous and reliable experiments are still needed for verification and analysis.

    Description of the author's contribution: the design idea of this paper is provided by professors Wang Rui and Yang Jing; Bai Haotian is responsible for target screening, pathway analysis and differential gene expression analysis; Li Yalan is responsible for data visualization and data processing.

    国产男靠女视频免费网站| 热99re8久久精品国产| 欧美一区二区国产精品久久精品| 免费看a级黄色片| 三级国产精品欧美在线观看| 好看av亚洲va欧美ⅴa在| 亚洲男人的天堂狠狠| tocl精华| 18禁美女被吸乳视频| 久久99热这里只有精品18| 国产探花极品一区二区| 免费观看的影片在线观看| 久久久久久久精品吃奶| 九九热线精品视视频播放| 精品乱码久久久久久99久播| 搡老岳熟女国产| 三级国产精品欧美在线观看| 精品一区二区三区人妻视频| 国语自产精品视频在线第100页| 久久精品国产综合久久久| 欧美大码av| 高潮久久久久久久久久久不卡| 国产精品永久免费网站| 免费一级毛片在线播放高清视频| 亚洲欧美日韩东京热| 九九在线视频观看精品| 国产乱人视频| avwww免费| 国产精品一区二区三区四区久久| 我要搜黄色片| 中国美女看黄片| 午夜精品久久久久久毛片777| 亚洲一区高清亚洲精品| 法律面前人人平等表现在哪些方面| 黄色丝袜av网址大全| 两个人看的免费小视频| 9191精品国产免费久久| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩精品网址| 亚洲精品在线美女| 午夜福利高清视频| 老熟妇仑乱视频hdxx| 久久久久亚洲av毛片大全| 精品久久久久久,| 亚洲欧美日韩无卡精品| 国产伦精品一区二区三区视频9 | 欧美高清成人免费视频www| bbb黄色大片| 亚洲国产精品久久男人天堂| 看黄色毛片网站| 久久精品国产亚洲av涩爱 | 国产av一区在线观看免费| 亚洲欧美日韩高清在线视频| 亚洲人成网站高清观看| 亚洲午夜理论影院| 国产色婷婷99| 免费av观看视频| 国产成+人综合+亚洲专区| 99久久九九国产精品国产免费| av福利片在线观看| 12—13女人毛片做爰片一| 久久国产乱子伦精品免费另类| av欧美777| 国产精品,欧美在线| 国产真人三级小视频在线观看| 久久欧美精品欧美久久欧美| 亚洲美女视频黄频| 变态另类丝袜制服| 嫩草影院精品99| 免费看十八禁软件| 免费在线观看成人毛片| av黄色大香蕉| 欧美高清成人免费视频www| 激情在线观看视频在线高清| 每晚都被弄得嗷嗷叫到高潮| 97人妻精品一区二区三区麻豆| 国产精品久久久久久精品电影| 天堂动漫精品| 性色av乱码一区二区三区2| 给我免费播放毛片高清在线观看| 在线免费观看不下载黄p国产 | 在线视频色国产色| 日本 av在线| 宅男免费午夜| 欧美一区二区亚洲| 给我免费播放毛片高清在线观看| 99热6这里只有精品| 在线播放国产精品三级| 国产精品久久电影中文字幕| 看免费av毛片| 熟妇人妻久久中文字幕3abv| 欧美三级亚洲精品| 波多野结衣高清无吗| 精华霜和精华液先用哪个| 国产av不卡久久| 亚洲av电影在线进入| 性色avwww在线观看| 黄色片一级片一级黄色片| 人妻久久中文字幕网| 成人18禁在线播放| 欧美高清成人免费视频www| 在线观看66精品国产| av中文乱码字幕在线| 午夜免费观看网址| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全免费视频| 99在线视频只有这里精品首页| 免费看日本二区| 男女那种视频在线观看| 国产精品亚洲美女久久久| 亚洲国产高清在线一区二区三| 国产精品久久久久久亚洲av鲁大| 色吧在线观看| 一夜夜www| 精品一区二区三区视频在线 | 亚洲成人久久性| 成人三级黄色视频| 天堂动漫精品| 精品久久久久久久末码| 少妇人妻精品综合一区二区 | 一二三四社区在线视频社区8| 中文字幕熟女人妻在线| 俺也久久电影网| 亚洲欧美激情综合另类| 最近最新中文字幕大全电影3| 欧美高清成人免费视频www| 黄色成人免费大全| 精品人妻1区二区| 宅男免费午夜| 欧美日韩乱码在线| 精品欧美国产一区二区三| 色视频www国产| 日本一本二区三区精品| 国产精品国产高清国产av| ponron亚洲| 亚洲美女视频黄频| 欧美性猛交╳xxx乱大交人| 中国美女看黄片| 亚洲欧美激情综合另类| 国产成人a区在线观看| 99国产精品一区二区三区| 999久久久精品免费观看国产| 美女免费视频网站| 中文字幕人妻丝袜一区二区| 免费大片18禁| 中国美女看黄片| 国内精品久久久久久久电影| 国产在视频线在精品| 日本熟妇午夜| 亚洲一区二区三区色噜噜| 在线观看一区二区三区| 免费在线观看影片大全网站| 热99在线观看视频| 国产成人福利小说| 中文字幕熟女人妻在线| 黄片大片在线免费观看| 亚洲男人的天堂狠狠| 亚洲国产欧美网| 九九在线视频观看精品| 五月伊人婷婷丁香| 亚洲av电影在线进入| 国产欧美日韩一区二区精品| 内地一区二区视频在线| 国产亚洲精品一区二区www| 国产一区二区三区视频了| 啦啦啦免费观看视频1| 国产成人福利小说| 久久久久国内视频| 国产v大片淫在线免费观看| 国产精品野战在线观看| 国产一区二区三区视频了| 久久这里只有精品中国| 国内毛片毛片毛片毛片毛片| 欧美最新免费一区二区三区 | 国产成人av激情在线播放| 又黄又粗又硬又大视频| 日韩精品青青久久久久久| 亚洲av美国av| 色精品久久人妻99蜜桃| 性欧美人与动物交配| av专区在线播放| 国产极品精品免费视频能看的| 一区福利在线观看| 日韩欧美免费精品| 亚洲人成电影免费在线| 欧洲精品卡2卡3卡4卡5卡区| www日本在线高清视频| 午夜福利18| 一二三四社区在线视频社区8| 久久精品91无色码中文字幕| 精品日产1卡2卡| 91在线观看av| 国产免费一级a男人的天堂| www.色视频.com| 亚洲av一区综合| 欧美3d第一页| 国产激情偷乱视频一区二区| 99久久九九国产精品国产免费| 日韩欧美国产一区二区入口| 国产精品一区二区三区四区久久| 在线观看午夜福利视频| 亚洲国产中文字幕在线视频| 无限看片的www在线观看| 亚洲精品亚洲一区二区| 少妇的逼好多水| 亚洲欧美日韩高清专用| 中文字幕人妻熟人妻熟丝袜美 | 久久这里只有精品中国| 久久香蕉国产精品| 亚洲五月婷婷丁香| 制服人妻中文乱码| 在线十欧美十亚洲十日本专区| 一区二区三区激情视频| www.色视频.com| 亚洲欧美一区二区三区黑人| 热99re8久久精品国产| 国产激情偷乱视频一区二区| 我的老师免费观看完整版| 国产精品一区二区三区四区久久| 性色avwww在线观看| 人妻久久中文字幕网| 国产熟女xx| 亚洲乱码一区二区免费版| 亚洲欧美一区二区三区黑人| 日本黄色片子视频| 日日摸夜夜添夜夜添小说| 欧美激情久久久久久爽电影| 国产男靠女视频免费网站| 久久香蕉国产精品| 在线观看午夜福利视频| 亚洲国产中文字幕在线视频| 激情在线观看视频在线高清| 村上凉子中文字幕在线| 精品国产美女av久久久久小说| 老司机午夜福利在线观看视频| 757午夜福利合集在线观看| 欧美+日韩+精品| 午夜免费成人在线视频| 欧美一区二区国产精品久久精品| 久久久国产成人免费| 国产在线精品亚洲第一网站| 女同久久另类99精品国产91| 久久精品国产亚洲av香蕉五月| 精品久久久久久久毛片微露脸| 亚洲人成伊人成综合网2020| 中文字幕高清在线视频| 一级黄色大片毛片| 久久精品夜夜夜夜夜久久蜜豆| 男女之事视频高清在线观看| 欧美一区二区国产精品久久精品| 好男人在线观看高清免费视频| av黄色大香蕉| 女同久久另类99精品国产91| 亚洲美女视频黄频| 午夜老司机福利剧场| 欧美大码av| 精品一区二区三区视频在线观看免费| 在线免费观看的www视频| 亚洲成av人片免费观看| 99热这里只有是精品50| 久久久久久久久中文| av在线蜜桃| 麻豆国产97在线/欧美| 99在线视频只有这里精品首页| 亚洲在线观看片| 国产午夜福利久久久久久| 久久精品国产清高在天天线| 亚洲精品影视一区二区三区av| 日本在线视频免费播放| 欧美日韩综合久久久久久 | 日韩欧美三级三区| 岛国在线免费视频观看| 一级黄片播放器| a级一级毛片免费在线观看| 一进一出抽搐gif免费好疼| 国产成人福利小说| 国产精品美女特级片免费视频播放器| 成人欧美大片| 我的老师免费观看完整版| 亚洲精品色激情综合| 久久午夜亚洲精品久久| 啦啦啦观看免费观看视频高清| 两个人的视频大全免费| 19禁男女啪啪无遮挡网站| 国产精品嫩草影院av在线观看 | 色在线成人网| 成年女人永久免费观看视频| 国产美女午夜福利| 真实男女啪啪啪动态图| 国产一区二区三区视频了| 久久精品综合一区二区三区| 老司机福利观看| 亚洲人成网站在线播放欧美日韩| 18禁国产床啪视频网站| 久久九九热精品免费| 99热精品在线国产| 国产精品久久久久久亚洲av鲁大| 少妇裸体淫交视频免费看高清| 成年免费大片在线观看| 搞女人的毛片| 男女下面进入的视频免费午夜| 欧美激情在线99| 亚洲天堂国产精品一区在线| 精品国产超薄肉色丝袜足j| 亚洲精品美女久久久久99蜜臀| 无人区码免费观看不卡| 国产亚洲精品av在线| 男人舔女人下体高潮全视频| 午夜福利在线观看吧| 国产久久久一区二区三区| 全区人妻精品视频| 十八禁网站免费在线| 可以在线观看的亚洲视频| 中文字幕高清在线视频| 2021天堂中文幕一二区在线观| 欧美激情在线99| xxx96com| 亚洲av成人不卡在线观看播放网| 久久久久久久午夜电影| 很黄的视频免费| 免费看光身美女| 一区福利在线观看| 美女cb高潮喷水在线观看| 99久久99久久久精品蜜桃| 日韩成人在线观看一区二区三区| 亚洲精品影视一区二区三区av| 制服人妻中文乱码| 成人av一区二区三区在线看| 日韩欧美免费精品| 中文字幕人妻熟人妻熟丝袜美 | 精品无人区乱码1区二区| 一个人看的www免费观看视频| 有码 亚洲区| 蜜桃久久精品国产亚洲av| 亚洲成人久久性| 日日干狠狠操夜夜爽| 99热这里只有是精品50| 又粗又爽又猛毛片免费看| 一个人免费在线观看电影| 国产精品永久免费网站| 国产三级中文精品| 午夜亚洲福利在线播放| 日本三级黄在线观看| 精品人妻偷拍中文字幕| 成人三级黄色视频| 亚洲精华国产精华精| 99精品在免费线老司机午夜| 亚洲国产高清在线一区二区三| 精品乱码久久久久久99久播| av女优亚洲男人天堂| 亚洲午夜理论影院| 3wmmmm亚洲av在线观看| 日韩欧美精品免费久久 | 搞女人的毛片| 精品一区二区三区人妻视频| 久久久久久久午夜电影| 欧美色欧美亚洲另类二区| 亚洲色图av天堂| 国产高清视频在线观看网站| 欧美另类亚洲清纯唯美| 免费在线观看影片大全网站| 制服丝袜大香蕉在线| 日本在线视频免费播放| 国内精品美女久久久久久| 精品久久久久久久毛片微露脸| 一卡2卡三卡四卡精品乱码亚洲| 变态另类成人亚洲欧美熟女| АⅤ资源中文在线天堂| 美女 人体艺术 gogo| 757午夜福利合集在线观看| АⅤ资源中文在线天堂| 99热6这里只有精品| 又黄又粗又硬又大视频| 国产91精品成人一区二区三区| 国产乱人视频| 88av欧美| 久久6这里有精品| 久久国产精品人妻蜜桃| 夜夜看夜夜爽夜夜摸| 国产精品精品国产色婷婷| 欧美成人免费av一区二区三区| 午夜两性在线视频| 九九热线精品视视频播放| 欧美另类亚洲清纯唯美| 成人国产一区最新在线观看| 男女之事视频高清在线观看| www.熟女人妻精品国产| 成年女人看的毛片在线观看| 深爱激情五月婷婷| 国产高清videossex| 欧美一区二区精品小视频在线| 俄罗斯特黄特色一大片| 天堂影院成人在线观看| 欧美又色又爽又黄视频| 色精品久久人妻99蜜桃| 亚洲中文日韩欧美视频| 色吧在线观看| 特大巨黑吊av在线直播| 一级黄色大片毛片| 久久精品综合一区二区三区| 精品久久久久久成人av| 国产老妇女一区| 久久久精品欧美日韩精品| 亚洲国产中文字幕在线视频| 精品一区二区三区视频在线观看免费| 少妇的逼好多水| 欧美极品一区二区三区四区| 男人舔奶头视频| 1024手机看黄色片| 女同久久另类99精品国产91| 成人特级av手机在线观看| 亚洲天堂国产精品一区在线| 国内揄拍国产精品人妻在线| 美女cb高潮喷水在线观看| 9191精品国产免费久久| 一区福利在线观看| 别揉我奶头~嗯~啊~动态视频| 特级一级黄色大片| 高清日韩中文字幕在线| 免费在线观看成人毛片| 黄色日韩在线| 午夜a级毛片| 久久久国产成人免费| 久久精品国产亚洲av涩爱 | 免费av观看视频| 久久精品国产清高在天天线| 欧美丝袜亚洲另类 | 亚洲精华国产精华精| 人人妻人人看人人澡| 黄色女人牲交| 免费在线观看亚洲国产| 国产老妇女一区| 男女之事视频高清在线观看| 国产色婷婷99| 国产精品av视频在线免费观看| 黄色片一级片一级黄色片| 婷婷精品国产亚洲av| 欧美日韩中文字幕国产精品一区二区三区| 99久久成人亚洲精品观看| 久久人人精品亚洲av| 亚洲精品久久国产高清桃花| 99国产精品一区二区三区| 听说在线观看完整版免费高清| 国产高清三级在线| 久久亚洲精品不卡| av片东京热男人的天堂| 成年女人看的毛片在线观看| av天堂中文字幕网| 久久99热这里只有精品18| 99久久精品国产亚洲精品| 成人av一区二区三区在线看| 一区二区三区免费毛片| 啦啦啦韩国在线观看视频| 国产蜜桃级精品一区二区三区| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 青草久久国产| 日韩精品中文字幕看吧| 欧美日韩中文字幕国产精品一区二区三区| www.熟女人妻精品国产| 亚洲乱码一区二区免费版| 亚洲欧美日韩高清专用| 欧美最新免费一区二区三区 | 91久久精品电影网| 中亚洲国语对白在线视频| 国产伦精品一区二区三区四那| www.色视频.com| 欧美一级毛片孕妇| 99国产极品粉嫩在线观看| 久久久久九九精品影院| 白带黄色成豆腐渣| 免费观看人在逋| 在线播放无遮挡| 九色成人免费人妻av| 在线播放国产精品三级| 3wmmmm亚洲av在线观看| 国产69精品久久久久777片| 人妻夜夜爽99麻豆av| 在线观看午夜福利视频| 亚洲最大成人中文| 超碰av人人做人人爽久久 | 成人国产一区最新在线观看| 久久久色成人| 久久国产精品影院| 久久亚洲真实| 日韩精品中文字幕看吧| 我要搜黄色片| 欧美乱色亚洲激情| 亚洲狠狠婷婷综合久久图片| 欧美一级a爱片免费观看看| 国产私拍福利视频在线观看| 蜜桃久久精品国产亚洲av| 亚洲成人精品中文字幕电影| 欧美极品一区二区三区四区| 亚洲精品亚洲一区二区| 午夜精品在线福利| 久久香蕉精品热| 日韩欧美精品免费久久 | 一个人看的www免费观看视频| 制服丝袜大香蕉在线| 亚洲国产日韩欧美精品在线观看 | 啦啦啦观看免费观看视频高清| 精品国产美女av久久久久小说| svipshipincom国产片| 久久99热这里只有精品18| 亚洲电影在线观看av| 欧美成人性av电影在线观看| 国产精品电影一区二区三区| 国产私拍福利视频在线观看| 亚洲久久久久久中文字幕| 久久国产乱子伦精品免费另类| 午夜激情福利司机影院| 99久久久亚洲精品蜜臀av| 亚洲最大成人手机在线| 精品不卡国产一区二区三区| 久久久久久久精品吃奶| 又黄又粗又硬又大视频| 国产一区二区亚洲精品在线观看| 真实男女啪啪啪动态图| 久久99热这里只有精品18| 婷婷丁香在线五月| 99国产极品粉嫩在线观看| 日韩有码中文字幕| 成人国产一区最新在线观看| 国产真实乱freesex| 99久久久亚洲精品蜜臀av| 亚洲av免费高清在线观看| 成人av一区二区三区在线看| 亚洲av免费高清在线观看| 动漫黄色视频在线观看| 嫩草影院入口| 一个人看的www免费观看视频| 久久6这里有精品| 欧美乱色亚洲激情| 免费看美女性在线毛片视频| 欧美高清成人免费视频www| 欧美性猛交╳xxx乱大交人| 国产亚洲精品av在线| 欧美一区二区国产精品久久精品| 精品免费久久久久久久清纯| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 久久精品国产99精品国产亚洲性色| 欧美日韩综合久久久久久 | 他把我摸到了高潮在线观看| 18禁黄网站禁片午夜丰满| 国产免费一级a男人的天堂| 欧美性感艳星| 无限看片的www在线观看| 免费观看人在逋| 亚洲av中文字字幕乱码综合| 欧美一区二区国产精品久久精品| 色播亚洲综合网| 网址你懂的国产日韩在线| 舔av片在线| 一个人看视频在线观看www免费 | 一卡2卡三卡四卡精品乱码亚洲| 一级a爱片免费观看的视频| 好男人在线观看高清免费视频| 国产高清视频在线播放一区| 叶爱在线成人免费视频播放| 久久精品国产亚洲av涩爱 | 欧美3d第一页| 美女大奶头视频| 蜜桃久久精品国产亚洲av| 久久这里只有精品中国| 99久久成人亚洲精品观看| 精品乱码久久久久久99久播| 午夜精品一区二区三区免费看| 免费高清视频大片| 亚洲中文日韩欧美视频| 国产精品爽爽va在线观看网站| 免费av观看视频| 久久久精品欧美日韩精品| 欧美日本视频| 亚洲av第一区精品v没综合| 波野结衣二区三区在线 | 天堂网av新在线| 亚洲av免费在线观看| av专区在线播放| 一个人免费在线观看的高清视频| 亚洲一区二区三区不卡视频| 免费观看人在逋| 少妇高潮的动态图| 国产高清视频在线观看网站| 麻豆一二三区av精品| 深爱激情五月婷婷| 国产伦在线观看视频一区| АⅤ资源中文在线天堂| 99久久精品一区二区三区| 免费电影在线观看免费观看| 欧美最黄视频在线播放免费| 亚洲成人精品中文字幕电影| 亚洲久久久久久中文字幕| 国产免费av片在线观看野外av| 老司机福利观看| 男女做爰动态图高潮gif福利片| 美女高潮的动态| 精品国产三级普通话版| 国产麻豆成人av免费视频| 国产老妇女一区| 亚洲av免费高清在线观看| 人妻夜夜爽99麻豆av| 午夜福利免费观看在线| 亚洲不卡免费看| 国产成年人精品一区二区| 窝窝影院91人妻| 少妇的逼好多水| 好看av亚洲va欧美ⅴa在| 国产一区二区激情短视频| 9191精品国产免费久久| 国产成年人精品一区二区| 9191精品国产免费久久| 欧美激情在线99| 搞女人的毛片|