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

    97人妻天天添夜夜摸| 午夜精品国产一区二区电影| 亚洲 国产 在线| 成人午夜精彩视频在线观看| 国产有黄有色有爽视频| 亚洲成国产人片在线观看| 91精品国产国语对白视频| 久久久精品区二区三区| 国产真人三级小视频在线观看| 亚洲av日韩精品久久久久久密 | 高清欧美精品videossex| 超碰97精品在线观看| 国产精品久久久人人做人人爽| 精品免费久久久久久久清纯 | 无遮挡黄片免费观看| 国产欧美日韩精品亚洲av| 色94色欧美一区二区| 亚洲精品日韩在线中文字幕| 国产人伦9x9x在线观看| 九色亚洲精品在线播放| 久久毛片免费看一区二区三区| 国产免费福利视频在线观看| 色综合欧美亚洲国产小说| 91九色精品人成在线观看| 在线亚洲精品国产二区图片欧美| 人人妻人人爽人人添夜夜欢视频| 后天国语完整版免费观看| 久久午夜综合久久蜜桃| 侵犯人妻中文字幕一二三四区| 日韩电影二区| 性色av乱码一区二区三区2| 日韩一卡2卡3卡4卡2021年| 国产1区2区3区精品| 超碰97精品在线观看| 午夜福利,免费看| 丝袜人妻中文字幕| svipshipincom国产片| 色婷婷久久久亚洲欧美| 国产精品香港三级国产av潘金莲 | 丰满饥渴人妻一区二区三| 久久精品熟女亚洲av麻豆精品| 黄频高清免费视频| 一边亲一边摸免费视频| 97精品久久久久久久久久精品| 久久精品熟女亚洲av麻豆精品| 亚洲一码二码三码区别大吗| 男女国产视频网站| 亚洲中文字幕日韩| 久久久久精品人妻al黑| 一区二区av电影网| 国产精品久久久久成人av| 十八禁高潮呻吟视频| 成人亚洲精品一区在线观看| 七月丁香在线播放| 国语对白做爰xxxⅹ性视频网站| 新久久久久国产一级毛片| 亚洲中文av在线| a级毛片在线看网站| 亚洲欧美成人综合另类久久久| 国产片内射在线| 99国产精品一区二区三区| 亚洲图色成人| 国产欧美日韩精品亚洲av| 久久久亚洲精品成人影院| 国产一区二区在线观看av| 精品少妇久久久久久888优播| 亚洲欧美一区二区三区国产| 人人妻,人人澡人人爽秒播 | 亚洲欧美精品自产自拍| 国产福利在线免费观看视频| 日韩av在线免费看完整版不卡| 国产国语露脸激情在线看| 99久久精品国产亚洲精品| 国产一区二区 视频在线| 免费在线观看视频国产中文字幕亚洲 | 男女边吃奶边做爰视频| 97人妻天天添夜夜摸| 日韩视频在线欧美| 欧美亚洲 丝袜 人妻 在线| 黑人欧美特级aaaaaa片| 久久久久久免费高清国产稀缺| 免费看不卡的av| 国产黄色免费在线视频| 日韩一卡2卡3卡4卡2021年| av电影中文网址| 免费黄频网站在线观看国产| 老鸭窝网址在线观看| 欧美国产精品一级二级三级| 女人被躁到高潮嗷嗷叫费观| 老汉色av国产亚洲站长工具| 精品久久久久久电影网| 免费日韩欧美在线观看| 91麻豆精品激情在线观看国产 | 女警被强在线播放| 我要看黄色一级片免费的| 午夜影院在线不卡| 国产黄色免费在线视频| 一级毛片女人18水好多 | 国产成人免费观看mmmm| 曰老女人黄片| 亚洲久久久国产精品| 日本av免费视频播放| av又黄又爽大尺度在线免费看| 色婷婷久久久亚洲欧美| 欧美黑人欧美精品刺激| 国产极品粉嫩免费观看在线| 9191精品国产免费久久| 80岁老熟妇乱子伦牲交| 婷婷色综合大香蕉| 高清欧美精品videossex| avwww免费| 亚洲人成电影观看| 久久久久久久精品精品| videos熟女内射| 男女之事视频高清在线观看 | 极品少妇高潮喷水抽搐| 性色av乱码一区二区三区2| 亚洲国产中文字幕在线视频| 在线观看www视频免费| 精品人妻熟女毛片av久久网站| 亚洲色图 男人天堂 中文字幕| 亚洲,一卡二卡三卡| 亚洲av成人不卡在线观看播放网 | 可以免费在线观看a视频的电影网站| 国产老妇伦熟女老妇高清| 两个人免费观看高清视频| 18在线观看网站| 我的亚洲天堂| 亚洲少妇的诱惑av| 日本a在线网址| 国产亚洲午夜精品一区二区久久| 一区二区三区激情视频| 亚洲av日韩精品久久久久久密 | 亚洲av成人不卡在线观看播放网 | 亚洲国产中文字幕在线视频| 亚洲欧美激情在线| 老司机影院毛片| av网站在线播放免费| 麻豆乱淫一区二区| 成在线人永久免费视频| 国产视频首页在线观看| 亚洲图色成人| 天天操日日干夜夜撸| 国产视频首页在线观看| 涩涩av久久男人的天堂| 欧美日韩成人在线一区二区| 亚洲自偷自拍图片 自拍| 午夜福利乱码中文字幕| 久久精品久久久久久噜噜老黄| 又黄又粗又硬又大视频| 亚洲成av片中文字幕在线观看| 国产欧美日韩综合在线一区二区| www日本在线高清视频| 欧美变态另类bdsm刘玥| 菩萨蛮人人尽说江南好唐韦庄| 超碰成人久久| 91麻豆精品激情在线观看国产 | 欧美乱码精品一区二区三区| 亚洲视频免费观看视频| 18禁黄网站禁片午夜丰满| 午夜福利影视在线免费观看| 久久毛片免费看一区二区三区| 婷婷色综合大香蕉| 亚洲欧美精品综合一区二区三区| 亚洲成色77777| 熟女少妇亚洲综合色aaa.| 国产精品国产三级国产专区5o| 性高湖久久久久久久久免费观看| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品电影小说| 亚洲第一青青草原| 少妇人妻久久综合中文| 天堂俺去俺来也www色官网| 久久女婷五月综合色啪小说| www.精华液| 另类精品久久| 亚洲一码二码三码区别大吗| 色视频在线一区二区三区| 国产深夜福利视频在线观看| 男女免费视频国产| 日日爽夜夜爽网站| 成人午夜精彩视频在线观看| 欧美日韩一级在线毛片| 国精品久久久久久国模美| 午夜视频精品福利| 亚洲精品美女久久av网站| 国产熟女午夜一区二区三区| www.精华液| 黄色视频在线播放观看不卡| 国产野战对白在线观看| 亚洲精品美女久久av网站| av福利片在线| 看免费av毛片| 伊人久久大香线蕉亚洲五| 看免费成人av毛片| 久久久久久亚洲精品国产蜜桃av| 亚洲三区欧美一区| 午夜福利免费观看在线| 久久综合国产亚洲精品| 婷婷色av中文字幕| 国产免费又黄又爽又色| 国产精品亚洲av一区麻豆| 青春草亚洲视频在线观看| av片东京热男人的天堂| 黑人欧美特级aaaaaa片| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美一区二区三区国产| 国产一区有黄有色的免费视频| 精品少妇黑人巨大在线播放| 亚洲久久久国产精品| 亚洲欧美中文字幕日韩二区| 精品熟女少妇八av免费久了| 亚洲欧洲日产国产| 亚洲欧美日韩高清在线视频 | 大香蕉久久成人网| 青青草视频在线视频观看| 欧美人与性动交α欧美软件| 国产成人系列免费观看| 香蕉丝袜av| 在线av久久热| 热99国产精品久久久久久7| 天天影视国产精品| 你懂的网址亚洲精品在线观看| 久久久久久亚洲精品国产蜜桃av| 一本综合久久免费| 女人高潮潮喷娇喘18禁视频| 狠狠婷婷综合久久久久久88av| 欧美日韩视频高清一区二区三区二| 人人妻,人人澡人人爽秒播 | 又黄又粗又硬又大视频| 日韩熟女老妇一区二区性免费视频| 黑人猛操日本美女一级片| 国产在线一区二区三区精| 18在线观看网站| 夜夜骑夜夜射夜夜干| 国产成人一区二区在线| 免费av中文字幕在线| 十八禁高潮呻吟视频| 91国产中文字幕| 欧美国产精品va在线观看不卡| 亚洲国产日韩一区二区| 9色porny在线观看| av视频免费观看在线观看| 色婷婷av一区二区三区视频| 丝袜在线中文字幕| 老鸭窝网址在线观看| 9191精品国产免费久久| 人人妻人人爽人人添夜夜欢视频| 男女高潮啪啪啪动态图| 亚洲国产精品国产精品| 麻豆国产av国片精品| 美女国产高潮福利片在线看| 国产成人精品在线电影| 国产成人欧美在线观看 | 一区二区三区激情视频| 欧美老熟妇乱子伦牲交| 波多野结衣一区麻豆| 黑人欧美特级aaaaaa片| 久久精品国产a三级三级三级| av网站免费在线观看视频| 真人做人爱边吃奶动态| 成年人免费黄色播放视频| 一级毛片电影观看| 国产免费又黄又爽又色| 一本色道久久久久久精品综合| 久久精品久久久久久久性| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 1024香蕉在线观看| 黑人猛操日本美女一级片| 一个人免费看片子| 亚洲国产av新网站| 好男人电影高清在线观看| 只有这里有精品99| 久久人人爽人人片av| 你懂的网址亚洲精品在线观看| 国精品久久久久久国模美| 免费看不卡的av| 91麻豆av在线| 亚洲国产精品国产精品| 国产精品av久久久久免费| 日韩中文字幕视频在线看片| 久久久久国产精品人妻一区二区| 大香蕉久久网| 久久久久视频综合| 国产精品一二三区在线看| 国产欧美日韩精品亚洲av| 国产精品一国产av| 秋霞在线观看毛片| 免费看十八禁软件| 一级a爱视频在线免费观看| 男女之事视频高清在线观看 | 两个人免费观看高清视频| 亚洲精品乱久久久久久| 精品亚洲乱码少妇综合久久| 欧美日韩一级在线毛片| 老司机靠b影院| 最新的欧美精品一区二区| 2021少妇久久久久久久久久久| 国产三级黄色录像| 高清不卡的av网站| 国产在线免费精品| 日日摸夜夜添夜夜爱| 午夜福利乱码中文字幕| 大码成人一级视频| 女人高潮潮喷娇喘18禁视频| 日韩一卡2卡3卡4卡2021年| 黄色 视频免费看| 日韩av在线免费看完整版不卡| 欧美国产精品va在线观看不卡| 国产男女超爽视频在线观看| 中文字幕av电影在线播放| 国产精品久久久人人做人人爽| 亚洲av男天堂| 一级黄色大片毛片| 亚洲视频免费观看视频| 午夜免费成人在线视频| 成人亚洲精品一区在线观看| 亚洲av成人精品一二三区| 久久久久精品国产欧美久久久 | 成人手机av| 国产成人免费无遮挡视频| 欧美97在线视频| 欧美日韩成人在线一区二区| 丝袜美足系列| 欧美在线黄色| 后天国语完整版免费观看| 不卡av一区二区三区| 欧美老熟妇乱子伦牲交| 亚洲国产av影院在线观看| 国产黄色视频一区二区在线观看| 大香蕉久久网| 男女边吃奶边做爰视频| 三上悠亚av全集在线观看| 精品免费久久久久久久清纯 | 久久久久久人人人人人| 美女大奶头黄色视频| 2021少妇久久久久久久久久久| 99精国产麻豆久久婷婷| 又黄又粗又硬又大视频| 大话2 男鬼变身卡| 在现免费观看毛片| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| 婷婷色麻豆天堂久久| 国产黄色视频一区二区在线观看| 久久精品成人免费网站| 天堂俺去俺来也www色官网| 这个男人来自地球电影免费观看| 一级a爱视频在线免费观看| 久久久久久久久久久久大奶| 又紧又爽又黄一区二区| 国产黄频视频在线观看| 黄色 视频免费看| 午夜av观看不卡| 97人妻天天添夜夜摸| av有码第一页| 欧美激情 高清一区二区三区| 考比视频在线观看| 国产亚洲欧美在线一区二区| 色综合欧美亚洲国产小说| 女人被躁到高潮嗷嗷叫费观| 欧美精品一区二区大全| 色视频在线一区二区三区| 欧美av亚洲av综合av国产av| 午夜av观看不卡| 老司机影院成人| www.999成人在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕亚洲精品专区| 交换朋友夫妻互换小说| 亚洲色图综合在线观看| 亚洲国产看品久久| 99香蕉大伊视频| 日韩 亚洲 欧美在线| 天天添夜夜摸| 国产男女内射视频| 自线自在国产av| 国产欧美日韩综合在线一区二区| 亚洲国产精品国产精品| 少妇的丰满在线观看| 日本91视频免费播放| 各种免费的搞黄视频| 久久久国产欧美日韩av| 婷婷丁香在线五月| 另类亚洲欧美激情| 高潮久久久久久久久久久不卡| 久久精品国产亚洲av涩爱| 看免费av毛片| 免费高清在线观看日韩| 无限看片的www在线观看| 精品免费久久久久久久清纯 | 亚洲国产精品999| 国产成人av激情在线播放| 十八禁网站网址无遮挡| 亚洲成人手机| 国产熟女欧美一区二区| 美女国产高潮福利片在线看| 久久精品成人免费网站| 丝袜在线中文字幕| 国产一区二区在线观看av| 丁香六月天网| 日本猛色少妇xxxxx猛交久久| 赤兔流量卡办理| 成人18禁高潮啪啪吃奶动态图| 高清视频免费观看一区二区| 真人做人爱边吃奶动态| 免费女性裸体啪啪无遮挡网站| 人妻人人澡人人爽人人| 在线看a的网站| 国产片特级美女逼逼视频| 一级片免费观看大全| 亚洲欧美激情在线| 在线观看免费午夜福利视频| 人人妻,人人澡人人爽秒播 | 极品少妇高潮喷水抽搐| 又粗又硬又长又爽又黄的视频| 午夜免费观看性视频| 久久久精品94久久精品| 久久精品国产亚洲av涩爱| 国产免费视频播放在线视频| 人体艺术视频欧美日本| 丝袜脚勾引网站| 国产在线免费精品| 91精品国产国语对白视频| 19禁男女啪啪无遮挡网站| 日日爽夜夜爽网站| 中文字幕人妻丝袜制服| 久久国产精品影院| 如日韩欧美国产精品一区二区三区| 老司机靠b影院| 精品国产一区二区久久| 久久国产精品男人的天堂亚洲| 婷婷丁香在线五月| 另类亚洲欧美激情| 一区二区三区四区激情视频| 久久国产精品大桥未久av| 最近手机中文字幕大全| 日韩av在线免费看完整版不卡| 国产亚洲av高清不卡| 两性夫妻黄色片| 亚洲 国产 在线| 亚洲精品乱久久久久久| 国产精品国产av在线观看| 高潮久久久久久久久久久不卡| 国产免费现黄频在线看| 欧美性长视频在线观看| 亚洲人成电影观看| 91字幕亚洲| 精品一区二区三区四区五区乱码 | 欧美成狂野欧美在线观看| 飞空精品影院首页| 欧美黄色淫秽网站| 午夜福利乱码中文字幕| xxxhd国产人妻xxx| 精品国产乱码久久久久久男人| 在现免费观看毛片| 久久精品亚洲av国产电影网| 天天躁夜夜躁狠狠躁躁| 777米奇影视久久| 人妻 亚洲 视频| 丰满迷人的少妇在线观看| 女性生殖器流出的白浆| 在线 av 中文字幕| 亚洲美女黄色视频免费看| 亚洲av成人不卡在线观看播放网 | 9热在线视频观看99| 久久性视频一级片| av在线app专区| 午夜激情久久久久久久| 国产精品二区激情视频| 美国免费a级毛片| 天天操日日干夜夜撸| 久久精品久久久久久久性| 久久这里只有精品19| 美女中出高潮动态图| 久久午夜综合久久蜜桃| 国产精品一区二区在线观看99| 超碰成人久久| 久久精品亚洲av国产电影网| 午夜免费成人在线视频| 激情五月婷婷亚洲| 视频区欧美日本亚洲| 成年人午夜在线观看视频| 亚洲精品自拍成人| 亚洲国产成人一精品久久久| 亚洲成国产人片在线观看| 热99久久久久精品小说推荐| 精品高清国产在线一区| 伊人亚洲综合成人网| 每晚都被弄得嗷嗷叫到高潮| 日韩大片免费观看网站| 欧美激情 高清一区二区三区| 亚洲成人手机| 亚洲国产欧美一区二区综合| 久久久久久久国产电影| 美女扒开内裤让男人捅视频| 黄色视频不卡| 不卡av一区二区三区| 国产黄频视频在线观看| 丁香六月天网| 亚洲av电影在线观看一区二区三区| www.av在线官网国产| 无限看片的www在线观看| 亚洲精品久久久久久婷婷小说| 日韩制服丝袜自拍偷拍| 欧美另类一区| 国产成人精品在线电影| 久久久久国产精品人妻一区二区| 亚洲欧洲国产日韩| 久热这里只有精品99| 啦啦啦视频在线资源免费观看| 男女免费视频国产| 婷婷成人精品国产| 久久久久国产一级毛片高清牌| 深夜精品福利| 国产精品久久久久成人av| 大香蕉久久成人网| 亚洲欧美精品自产自拍| 在线观看www视频免费| 亚洲免费av在线视频| 赤兔流量卡办理| 91精品伊人久久大香线蕉| 欧美精品人与动牲交sv欧美| 国产爽快片一区二区三区| 一本—道久久a久久精品蜜桃钙片| 久久毛片免费看一区二区三区| 中文字幕高清在线视频| 纵有疾风起免费观看全集完整版| 国产欧美日韩一区二区三区在线| 午夜久久久在线观看| 国产在线视频一区二区| 不卡av一区二区三区| 后天国语完整版免费观看| 午夜影院在线不卡| 看免费av毛片| 国产女主播在线喷水免费视频网站| 欧美 日韩 精品 国产| 一级片'在线观看视频| 黄色视频不卡| 水蜜桃什么品种好| 久久这里只有精品19| 久久久欧美国产精品| 国产成人精品久久二区二区免费| 母亲3免费完整高清在线观看| 午夜视频精品福利| 夫妻午夜视频| 美女高潮到喷水免费观看| 国产1区2区3区精品| 成人亚洲欧美一区二区av| av不卡在线播放| 狠狠婷婷综合久久久久久88av| 国产在视频线精品| 五月开心婷婷网| 亚洲国产欧美网| 国产伦人伦偷精品视频| 国产精品久久久久久精品电影小说| 天天躁夜夜躁狠狠久久av| 一级黄片播放器| 一区二区三区精品91| 国产高清视频在线播放一区 | 午夜老司机福利片| videos熟女内射| 男女高潮啪啪啪动态图| 黄色视频不卡| 在现免费观看毛片| 十八禁人妻一区二区| 黄色 视频免费看| 精品一区二区三区四区五区乱码 | 男女下面插进去视频免费观看| 亚洲专区国产一区二区| 丝袜人妻中文字幕| 啦啦啦啦在线视频资源| 国产成人一区二区三区免费视频网站 | 伊人久久大香线蕉亚洲五| 欧美黄色淫秽网站| 人妻 亚洲 视频| xxxhd国产人妻xxx| 亚洲av电影在线进入| 9191精品国产免费久久| 国产女主播在线喷水免费视频网站| 国产视频一区二区在线看| 久久精品久久久久久噜噜老黄| 热99久久久久精品小说推荐| 最近手机中文字幕大全| 69精品国产乱码久久久| 80岁老熟妇乱子伦牲交| 午夜福利影视在线免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人免费观看高清视频| 国语对白做爰xxxⅹ性视频网站| 男人操女人黄网站| 人妻 亚洲 视频| 美国免费a级毛片| 观看av在线不卡| 午夜免费男女啪啪视频观看| 99久久99久久久精品蜜桃| 在线看a的网站| 久久精品久久久久久久性| 国产高清不卡午夜福利| 极品人妻少妇av视频| 在线观看人妻少妇| 国产成人免费观看mmmm| 欧美在线黄色| 久久热在线av| 夫妻性生交免费视频一级片| 一本大道久久a久久精品| 久久性视频一级片| 晚上一个人看的免费电影| 首页视频小说图片口味搜索 | 中文字幕人妻熟女乱码| 国产深夜福利视频在线观看| 一级片'在线观看视频| 亚洲男人天堂网一区| 精品国产国语对白av|