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

    Integrated bioinformatics analysis of potential biomarkers and candidate drugs of esophageal squamous cell carcinoma

    2022-08-23 01:03:00MingQianHaoTingZhaoDaYongChenFengXianZhangYiWenZhangQiTengDingShuWenSunJinPingZhangLingDongChuanBoDingWenCongLiu
    Medical Data Mining 2022年3期

    Ming-Qian Hao ,Ting Zhao ,Da-Yong Chen ,Feng-Xian Zhang ,Yi-Wen Zhang ,Qi-Teng Ding ,Shu-Wen Sun,Jin-Ping Zhang,Ling Dong,Chuan-Bo Ding,Wen-Cong Liu*

    1College of Traditional Chinese Medicine,Jilin Agricultural Science and Technology College,Jilin 132109,China.2School of Chinese Medicinal Materials,Jilin Agricultural University,Changchun 130118,China.3Jilin Province Center for Drug Inspection,Changchun 130062,China.

    Abstract Esophageal squamous cell carcinoma (ESCC),the major subtype of esophageal carcinoma(ESCA),is one of the most lethal malignancies worldwide.This study aimed to identify potential biomarkers and/or therapeutic targets for ESCC.The datasets GSE44021,GSE77861,GSE20347,and GSE29001 retrieved from the Gene Expression Omnibus (GEO)database contained 117 ESCC tissues and 109 normal tissues.Differentially expressed genes(DEGs) associated with ESCC were identified using the GEO2R tool.Dysregulated pathways associated with ESCC mainly included mitotic regulation,cell cycle,ECM-receptor interaction,DNA replication,etc.The protein-protein interaction (PPI) network of overlapping DEGs was constructed and nine key genes (KGs) were identified from the complex interaction network using Degree,maximum neighborhood component (MNC),and maximal clique centrality(MCC)algorithms.Expression patterns of KGs at the transcriptional and translational levels were validated using ESCC-related data from the Cancer Genome Atlas (TCGA),Oncomine,and Human Protein Atlas (HPA) databases.Genetic alterations calculation,immune cell infiltrates evaluation,methylation analysis,prognostic analysis,transcription factors (TFs) and miRNAs regulatory networks construction,and targeted drug prediction were further performed.It was also found that the knockout of these KGs affected the survival of more than two types of ESCC cells by genome-wide CRISPR-Cas9 dropout screens.In conclusion,we identified KGs,TFs,and miRNAs with biomarker potential (e.g.,NDC80,BUB1,TOP2A,AURKA,AURKB,TTK,UBE2C,TPX2,BUB1B,E2F1,and hsa-miR-483-5p) and 23 candidate targeted drugs for ESCC by utilizing an integrated multi-omics approach.These findings provide additional insights into uncovering the molecular mechanism and improving the efficiency of clinical diagnosis and treatment for ESCC.

    Keywords: esophageal squamous cell carcinoma;biomarker;transcription factor;miRNA;CRISPR-Cas9;candidate drug

    Background

    ESCC,the major pathological type of ESCA,is a global malignant tumor with poor diagnosis,prognosis,and survival rate [1].For the lack of effective preventive and diagnostic methods,most patients with ESCC are only diagnosed and detected at an advanced stage with a 5-year overall survival rate of less than 20% [2].In addition,although clinical therapies such as immunotherapy,neoadjuvant chemoradiotherapy,and minimally invasive esophagectomy have achieved some promising results,the risk of mortality from ESCC is still high due to population increase and aging,as well as incomplete treatment and tumor metastasis [3].Therefore,it is necessary to explore promising biomarkers and valuable therapeutic drugs for ESCC.

    High-throughput sequencing technology and bioinformatics method have become powerful tools to monitor genome-wide gene expression changes and provide an indispensable platform for integrating public databases and exploring the molecular mechanisms of various human diseases [4,5].Wang et al.re-utilized the published transcriptome data to identify DEGs and signaling pathways associated with ESCC,and further found that aberrant expression of MAPK1,ACOX1,and SCP2 affected the overall survival time of ESCC patients[6].Han et al.demonstrated that AURKA was an abnormally methylated DEGs in ESCC,and its overexpression was associated with poor overall survival through an integrated bioinformatics analysis [7].Moreover,a recent multi-omics report identified 10 hub nodes of ESCC by expression regulatory network analysis of abnormal genes and investigated their potential as diagnostic and therapeutic targets for ESCC [8].Another bioinformatics research based on large-scale sequencing data reported that miR-139-3p was significantly down-regulated in the serum of patients with ESCC/ESCA and might be a candidate biomarker for predicting ESCC [9].The above bioinformatics research mainly focuses on whether single or multiple genes affect tumorigenesis and prognosis.However,the related work on gene function and stability,tumor immune microenvironment,and drug candidates is insufficient and needs to be further implemented.

    Based on the above characteristics,we performed a comprehensive multi-omics task,including DEGs identification,PPI network construction,functional enrichment analysis,KGs identification,genetic alteration and methylation analysis,fitness gene prediction,tumor microenvironment and immune cell infiltration assessment,miRNAs-and TFs-KGs regulatory networks mapping,prognostic analysis,and targeted drug discovery,to reprocess large-scale sequencing data of ESCC from multiple human tumor-related databases.We believed this study could provide insights into the pathogenesis,clinical diagnosis,prognosis,and treatment of ESCC based on the identified biological mechanisms,biomarker targets,and drug candidates.

    Methods

    DEGs identification

    Four gene expression profiles of ESCC (GSE44021,GSE77861,GSE29001,and GSE20347) were randomly searched from the GEO(https://www.ncbi.nlm.nih.gov/geo/) database and DEGs between ESCC tissues and normal tissues were identified by using GEO2R tool(https://www.ncbi.nlm.nih.gov/geo/geo2r/)with a screening criteria of adjustedP-value (adj.P.Val) <0.05 and |log2 fold change (logFC)|>1 [10].Overlapping DEGs from different datasets were obtained using a VennDiagram package in R language (version 3.6.2) [11].Table S1 listed the details of the above datasets.

    PPI network construction and KGs identification

    PPI network was integrated using the STRING (version 11.0,https://string-db.org/) database with an interaction score of 0.7,and then visualized by the cytoscape software (version 3.7.2) [12-13].KGs were identified from the co-expression network using Degree,MNC and MCC algorithms of the CytoHubba plugin [14].

    Expression verification,methylation assessment and prognostic analysis of KGs

    The UALCAN (http://ualcan.path.uab.edu/index.html) and Oncomine(https://www.oncomine.org) databases were used to verify the transcriptional expression of KGs [15,16].Differences of these genes between normal and tumor tissues at the translational level were evaluated by the HPA (https://www.proteinatlas.org) database [17].We also assessed the methylation degree and prognostic value of KGs using the UALCAN tool.P< 0.05 was considered statistically significant.

    Identification of fitness genes in ESCC cell lines

    Fitness genes were defined as genes required for cell survival and growth.The Project Score database(https://score.depmap.sanger.ac.uk/),a public resource that contains genome-wide CRISPR-Cas9 dropout screening data in 323 cancer cell lines,was applied to evaluate fitness genes that promote the growth of three or more ESCC cell lines based on fitness scores [18].Corrected log fold change <-1 means that gene depletion prevents cell growth or viability.

    Functional enrichment analysis

    To elucidate the biological characteristics of overlapping DEGs,we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis using the Cluster Profiler package of R [19].Statistical significance was set atP<0.05.

    Evaluation of immune infiltration

    Tumor Immune Estimation Resource (TIMER,https://cistrome.shinyapps.io/timer/) is an integrated platform for analyzing 32 cancer types’immune infiltration data[20].We explored the link between abnormal expression of KGs and the number of tumor-infiltrating immune cells (B Cell,CD4+T Cell,CD8+T Cell,Macrophage,Neutrophil,and Dendritic Cell) in ESCA samples.Correlation analysis between KGs was performed by the “Correlation”module.

    Analysis of genetic alterations of KGs

    cBioPortal (http://cbioportal.org) facilitates the exploration and visualization of complex cancer genomics data[21].In this study,227 samples from 2 ESCC reports in the cBioPortal database were used to explore genetic alterations in KGs.

    Prediction of TFs and miRNAs-KGs interactions

    To understand the transcriptional regulation process of KGs,we constructed a TFs-KGs interaction network using TRRUST platform(version 2,https://www.grnpedia.org/trrust/) [22].In addition,the miRNAs-mRNAs network was constructed to further explore the function of KGs and the regulation of miRNAs in ESCC.Briefly,we identified overlapping differentially expressed miRNAs (DEmiRNAs)from datasets GSE114110 and GSE43732 utilizing the GEO2R and VennDiagram tools with criteria of|logFC|>1 and adj.P.Val <0.05.The interactions between KGs and overlapping DEmiRNAs were obtained from the miRDB (http://mirdb.org/),miRWalk(http://mirwalk.umm.uni-heidelberg.de/),and miRTarBase (http://mirtarbase.cuhk.edu.cn/php/index.php) databases [23-25].

    Targeted drugs prediction

    The DGIdb tool (http://www.dgidb.org/) can redesign and mine clinically relevant drug-gene interactions for personalized medicine[26].Hence,we further re-analyzed the compounds in the DGIdb database to identify anti-neoplastic drugs targeting KGs in ESCC.

    Results

    Identification of DEGs in ESCC

    In the section on DEGs identification,we randomly selected four datasets involving a total of 117 ESCC samples and 109 paired normal samples.The data boxplot showed that the distribution of samples in each individual dataset tended to be consistent and could be used for further analysis (Figure S1).As presented in Figure 1A-D,1048 DEGs in GSE44021 (537 up-and 511 down-regulated),1260 DEGs in GSE77861 (684 up-and 576 down-regulated),4858 DEGs in GSE29001(2952 up-and 1906 down-regulated),and 1428 DEGs(660 up-and 768 down-regulated) in GSE20347 were identified.To avoid the bias of individual studies,we performed a VennDiagram analysis and found 462 overlapping DEGs (212 up-and 250 down-regulated)from the above datasets (Figure 1E and F,Table S2).

    PPI network integration and KGs identification

    A PPI network composed of overlapping DEGs was shown in Figure 2A,containing 255 nodes and 1151 interactions.Subsequently,nine KGs,namely NDC80,BUB1,TOP2A,AURKA,AURKB,TTK,UBE2C,TPX2,and BUB1B,were identified based on the Degree,MNC,and MCC algorithms (Figure 2B-E).Notably,these KGs were all significantly up-regulated in malignant lesion tissues and might be inseparable from the pathogenesis and progression of ESCC.

    Expression verification,methylation assessment and prognostic analysis of KGs

    Based on clinical data and pan-cancer reports retrieved from the TCGA and oncomine databases,we found that these KGs in ESCC were abnormally up-regulated in most human cancers (e.g.,Brain and CNS Cancer,ESCA,Breast Cancer,Lung Cancer,Bladder Cancer,and Sarcoma,Figure 2F and Figure S2A).

    Figure 1 Identification of overlapping differentially expressed genes (DEGs) from four esophageal squamous cell carcinoma (ESCC)datasets.(A-D) Volcano plots of DEGs in GSE44021,GSE77861,GSE29001,and GSE20347,respectively.Red means up-regulated DEGs,green means down-regulated DEGs,and black means no significant difference.(E) Venn diagram of up-regulated DEGs.(F) Venn diagram of down-regulated DEGs.DEGs,differentially expressed genes;ESCC,esophageal squamous cell carcinoma.

    Figure 2 Protein-protein interaction (PPI) network construction and key genes (KGs) identification.(A) PPI network of overlapping DEGs.Red diamonds: up-regulated genes;Blue ellipses: down-regulated genes.(B-D) Top 10 genes identified from the PPI network via MCC,MNC,and Degree algorithms,respectively.(E) Venn diagram for identifying KGs.(F) Heatmap of KGs expression in esophageal carcinoma (ESCA) data from the cancer genome atlas(TCGA)database.PPI,protein-protein interaction;KGs,key genes;ESCA,esophageal carcinoma;TCGA,the cancer genome atlas.

    Further,we verified the expression trends of nine KGs in ESCC tissues by analyzing relevant samples from the TCGA database (Figure S2B-J).Interestingly,these genes were also over-expressed in esophageal adenocarcinoma (P<0.001).We then compared the expression differences of KGs at the translational level in ESCC and normal tissues using immunohistochemical staining data from the HPA database.TOP2A is over-expressed in ESCC tissues and moderately expressed in normal tissues (Figure 3A and B).AURKA showed medium expression in tumor samples but was negative in adjacent normal tissues (Figure 3C and D).AURKB displayed medium expression in both ESCC and adjacent tissues (Figure 3E and F).Medium expression of TPX2 was found in the nucleus of ESCC cells and the cytoplasm/membrane of normal tissue cells (Figure 3G and H).UBE2C was found to have high expression in normal tissues but moderately expressed in ESCC tissues,which might be relevant to the relatively small sample size of ESCC included in the HPA database(Figure 3I and J).Unfortunately,we did not find the expression data of NDC80,BUB1,TTK,and BUB1B proteins in ESCC from the HPA database.In summary,we predict that the abnormal expression of these proteins is highly likely to participate in the pathogenesis and progression of ESCC,which needs further experimental confirmation.

    Furthermore,we analyzed the methylation data of ESCC included in the TCGA database and found that the methylation levels of BUB1,UBE2C,and BUB1B were significantly higher in ESCC tissues than in normal tissues (P<0.01 orP<0.001);while the methylation levels of NDC80,BUB1,TOP2A,and TTK in esophageal adenocarcinoma were higher than those in normal tissues (P<0.05 orP<0.001,Figure S3).These findings indicated that KGs might be associated with the prognosis of ESCC,but current data showed that there was no significant relationship between the abnormal expression of these KGs and the overall survival of patients with ESCC (P>0.05,Figure S4).Therefore,we recognized that further evidence was needed for prognostic analysis of ESCC.

    Figure 3 Expression validation of key genes (KG) at the translational level by the Human Protein Atlas database.(A) TOP2A protein in normal tissues (staining: medium;intensity: strong;quantity: <25%;location: nuclear).(B) TOP2A protein in tumor tissues (staining: high;intensity:strong;quantity: 75%-25%;location:nuclear).(C)AURKA protein in normal tissues(staining:not detected;intensity:negative;quantity:none;location: none).(D) AURKA protein in tumor tissues (staining: medium;intensity: moderate;quantity: 75%-25%;location:cytoplasmic/membranous,nuclear).(E)AURKB protein in normal tissues(staining:medium;intensity:strong;quantity:<25%;location:nuclear).(F) AURKB protein in tumor tissues (staining: medium;intensity: strong;quantity: <25%;location: nuclear).(G) TPX2 protein in normal tissues(staining: medium;intensity: moderate;quantity: 75%-25%;location: cytoplasmic/membranous).(H) TPX2 protein in tumor tissues (staining:medium;intensity: strong;quantity: <25%;location: nuclear).(I) UBE2C protein in normal tissues (staining: high;intensity: strong;quantity: >75%;location: nuclear).(J) UBE2C protein in tumor tissues (staining: medium;intensity: moderate;quantity: >75%;location: nuclear).KG,key genes.

    Evaluation of fitness genes

    The impacts of knockout KGs on the growth and cell viability of ESCC cell lines were further investigated using genome-scale CRISPR-Cas9 dropout screening data.In the results,all nine KGs could be defined as fitness genes of ESCC,that is,the deletion of these genes hindered the survival and growth of more than two ESCC cell lines.Knockout of TOP2A or AURKB,but not other KGs,most significantly inhibited the growth of 19 ESCC cell lines.In addition,KYSE-450,TE-10,and T-T cells,but not other ESCC cell lines,were more sensitive to KGs deletion,although there was no statistical significance (Figure 4).

    GO and KEGG enrichment analysis

    GO and KEGG enrichment analysis was further performed to elaborate the biological annotations of 462 overlapping DEGs in ESCC.GO enrichment mainly included three categories: biological process (BP),cell component (CC),and molecular function (MF) (Figure 5A).For the BP category,overlapping DEGs were mainly enriched in cell cycle-related terms,such as extracellular matrix organization,cell cycle checkpoint,mitotic nuclear division,and positive regulation of cell cycle,which was consistent with the biological features of the rapid proliferation of ESCC cells.CC enrichment analysis indicated that overlapping DEGs were related to the collagen-containing extracellular matrix,chromosomal region,endoplasmic reticulum lumen,and extracellular matrix component.In the MF category,overlapping DEGs were markedly involved in growth factor binding,extracellular matrix structural constituent,cell adhesion molecule binding,DNA replication origin binding,and so on.KEGG pathway analysis revealed that 462 overlapping DEGs were primarily linked to DNA replication,cell cycle,transcriptional mis-regulation in cancer,ECM-receptor interaction,protein digestion and absorption,and proteoglycans in cancer (Figure 5B).Furthermore,KGs such as BUB1B,TTK,and BUB1 mainly participated in the cell cycle.

    These results were consistent with the understanding that abnormalities in cell growth regulatory factors and cell cycle were the main reasons for tumor formation,which further indicated that KGs and overlapping DEGs played crucial roles in the pathogenesis and tumor microenvironment of ESCC [27].

    Characterization of immune infiltration and genetic changes of key genes

    Tumor tissues contain not only tumor cells but also tumor-infiltrating immune cells closely related to tumor progression.In this study,nine KGs were correlated with purity and six types of immune-infiltrating cells (B Cell,CD4+T Cell,CD8+T Cell,Macrophage,Neutrophil,and Dendritic Cell) to varying degrees(Figure S5).Furthermore,there was a positive correlation in expression patterns among nine KGs in ESCC tissues,showing that these genes may be essential functional partners (Figure S6).

    Figure 4 Effects of knockout of key genes (KGs) on cells of esophageal squamous cell carcinoma (ESCC).Corrected log fold change <-1 means that gene depletion prevented cell growth or viability.KGs,key genes;ESCC,esophageal squamous cell carcinoma.

    Figure 5 Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis for overlapping differentially expressed genes (DEGs).(A) GO enrichment analysis of overlapping DEGs (Top 10 of each category).BP,biological process;CC,cellular component;MF,molecular function.The color represents P-value and the size of the spots represents the gene number.(B) KEGG pathway enrichment analysis.GO,gene ontology;KEGG,kyoto encyclopedia of genes and genomes;DEGs,differentially expressed genes.

    Genetic alterations of 9 KGs in 227 ESCC patients were analyzed using the cBioPortal database.As a result,we observed that AURKA and BUB1B produced 0.9% alterations due to missense mutations in ESCC samples,respectively,and the genetic alteration rate of the remaining genes was 0% (Figure 6A).Therefore,these KGs were considered low-frequency mutant genes.

    TFs-and miRNAs-KGs regulatory networks construction

    TFs-KGs network containing 7 KGs (TOP2A,NDC80,AURKA,AURKB,TTK,UBE2C,and BUB1B) and 12 TFs(E2F1,E2F3,E2F4,MED1,etc.)was constructed (Figure 6B).We found that other TFs excluding ATF4 and OTX2 were significantly up-regulated in ESCC samples compared to normal tissues (P<0.001,Figure S7).Interestingly,E2F1 not only activated AURKA and TOP2A genes but also affected the transcriptional process of AURKB.E2F3,MED1,and OTX2 as TFs all had potential activating effects on AURKA.In addition,ATF4 and CREB1 might activate the transcriptional process of NDC80,and ZNF143 and MYC had an activating effect on the expression of BUB1B and UBE2C,respectively.Conversely,E2F4 might repress TTK expression at the transcriptional level.

    Fifty-seven overlapping DEmiRNAs (36 up-and 21 down-regulated)were identified and obtained from datasets GSE114110 and GSE43732(Table S3).A miRNAs-KGs regulatory network comprising 7 KGs,16 overlapping DEmiRNAs (13 up-and 3 down-regulated),and 48 interactions was subsequently constructed (Figure 6C).Survival analysis indicated that the high expression of hsa-miR-483 was negatively correlated with the overall survival of patients with ESCC(P<0.0001,Figure 6D).Other miRNAs (hsa-miR-183,hsa-miR-145,hsa-miR-139,hsa-miR-708,hsa-miR-7,hsa-miR-96,etc.) presented similar prognostic outcomes in patients with ESCC,although there was no statistical difference (Figure S8).

    Targeted drugs analysis

    A total of 23 potential anti-ESCC compounds targeting KGs were identified,and most small molecular drugs (20/23) might target TOP2A,AURKA,AURKB,and TTK for the treatment of ESCC (Table 1).Co-inhibitors of AURKA and AURKB included AMG-900,AT-9283,barasertib,danusertib,GSK-1070916,MK-5108,PF-03814735,and SNS-314.

    Table 1 Candidate drugs targeting key genes (KGs) in esophageal squamous cell carcinoma (ESCC)

    Discussion

    Although numerous studies on the molecular mechanism,early diagnosis,and clinical treatment of ESCC have been performed in the past few years,the morbidity and mortality of ESCC are still increasing globally [3].The main reason for the above may be the false positive or negative rates of a single independent analysis.In this study,we used four independent microarray datasets from different reports in the GEO database to obtain more accurate and reliable information.Four hundred and sixty-two overlapping DEGs (212 upand 250 down-regulated) were identified for in-depth exploration of the biological mechanisms and KGs of ESCC (Figure 1E and F).Cell cycle,transcriptional mis-regulation in cancer,ECM-receptor interaction,DNA replication,protein digestion and absorption,proteoglycans in cancer,and biosynthesis of unsaturated fatty acids have been previously shown to be involved in ESCA [27],and we found that they may also participate in the carcinogenesis and progression of ESCC(Figure 5B).Numerous publications have pointed out that abnormal DNA replication and cell cycle are important signals for the initiation and progression of various cancers,including ESCC[28].ECM is an essential ingredient of the cancer cell niche and the first barrier against tumor invasion and metastasis [29].Consistent with previous studies,overlapping DEGs associated with the collagen family in this work (e.g.,COL1A1,COL1A2,COL4A2,and COL4A1)were enriched in the ECM-receptor interaction pathway (Figure 5B),which undoubtedly exacerbated malignancy growth,angiogenesis,and invasion [30].Remodeling of energy metabolism is a symptom of malignancy.Therefore,it is inevitable that ESCC cells create a more suitable environment for growth and invasion by activating pathways related to substance synthesis and energy metabolism,such as protein digestion and absorption,proteoglycans in cancer,and biosynthesis of unsaturated fatty acids [31].

    Figure 6 Genetic alterations analysis,and exploration of Transcription factors (TFs) and miRNAs regulatory networks of nine key genes(KGs).(A) Genetic alterations evaluation of KGs.(B) TFs-KGs regulatory network.Red rhomboids: KGs.Red triangles: up-regulated TFs in ESCC tissues compared with normal tissues.Gray triangles: TFs with no significant expression changes in ESCC tissues compared with normal tissues.Delta-shaped arrows:activation of KGs by TFs.Half circular arrows:repression of KGs by TFs.(C)miRNAs-KGs regulatory network.Red rhomboids:KGs.Red rectangles: up-regulated DEmiRNAs.Blue rectangles: down-regulated DEmiRNAs.(D) Prognostic analysis of hsa-miR-483.TFs,transcription factors;KGs,key genes.

    NDC80,BUB1,TOP2A,AURKA,AURKB,TTK,UBE2C,TPX2,and BUB1B were identified as KGs in ESCC using multiple machine algorithms,and their expression changes in diseased tissues were validated using a large number of ESCC-related data recorded in the Oncomine and TCGA databases (Figure 2 and Figure S2).The expression of TOP2A and AURKA in ESCC tissues was significantly higher compared to normal tissues at the translational level (Figure 3A-D).The methylation levels of BUB1,UBE2C,and BUB1B were higher in ESCC than those in normal samples (Figure S3).We also evaluated the inhibitory effect of KGs knockout on the survival and growth of 19 ESCC cell lines using genome-wide CRISPR-Cas9 dropout screens (Figure 4).In addition,nine KGs were significantly associated with immune cell infiltration,and they were also strongly positively correlated with each other (Figure S5 and S6),demonstrating that these KGs might be functional complexes that play a vital coordinating role in manipulating the immune microenvironment of ESCC.

    We performed a literature review of nine promising biomarkers for ESCC.NDC80 is a mitotic protein that interacts with other proteins to regulate the cell cycle.Studies have shown that NDC80 is over-expressed in colorectal cancer and clear cell renal cell carcinoma[32,33].Based on current findings,we speculate that NDC80 plays a potentially critical role in the carcinogenesis of ESCC.BUB1 and BUB1B have similar functions in the cell cycle.Currently,the role of BUB1 as an oncogene in human cancers including gastric cancer and liver cancer has been observed [34].The clinical value of BUB1B as a potential target for the diagnosis,treatment,and prognosis of hepatocellular carcinoma was previously demonstrated [35].TOP2A is considered a potential biomarker for the diagnosis and prognosis of patients with malignant tumors [36].Both AURKA and AURKB are mitotic regulators,which participate in the regulation of chromosome arrangement and separation during mitosis and meiosis.The over-expression of AURKA is related to the differentiation degree,invasive ability,distant lymph node metastasis,and poor prognosis of ESCC [7].AURKB has been reported to affect the occurrence of ESCC[37].The above studies have pointed out the application direction of AURKA and AURKB,although there is no clear indication that they can be used as clinical diagnosis or treatment targets for ESCC.TTK is a core unit of the spindle assembly checkpoint.Several reports indicate that the abnormally expressed TTK gene is a potential prognostic and/or therapeutic target for human malignancies such as thyroid cancer,breast cancer,liver cancer,and glioblastoma [38].Data from the HPA database showed that UBE2C protein was over-expressed in normal esophageal squamous epithelial cells and moderately expressed in ESCC tissues.Conversely,Palumbo et al.found that UBE2C protein was highly expressed in ESCC samples,and its silencing reduced the malignant phenotype of ESCC cells [39].In recent years,increasing reports have found that high expression of TPX2 further leads to more active proliferation of tumor cells [40].Therefore,it is believed that KGs may be potential oncogenes or driver genes of ESCC,and suppressing their expression can directly or indirectly inhibit the growth and invasion of ESCC cells.

    During the exploration of the TFs-KGs network,we identified 12 TFs targeting KGs,of which E2F1,E2F3,E2F4,and MED1 had the most target KGs (Figure 6B).Numerous reports indicate that E2F family members are essential for coordinating cell cycle progression [41].E2F1 is a transcriptional activator of AURKA and AURKB,and its over-expression is a vital indicator of tumor progression and prognosis in patients with ESCC [42].Furthermore,AURKA over-expression increases the stability and transcriptional activity of E2F1 protein in cancer samples by inhibiting proteasome-dependent protein degradation[43].The current study showed that E2F1,as an activator of AURKA,AURKB,and TOP2A genes,aggravated the deterioration of ESCC.E2F4 had a potential reversal effect on ESCC by inhibiting the expression of TTK.Elevated MED1 is a key molecular event related to the occurrence of prostate cancer[44].Similarly,we found that MED1 was abnormally up-regulated in ESCC tissues and might have an activating function for AURKA.

    Candidate miRNAs,as oncogenes or tumor suppressors,play important roles in the detection,treatment,and prognosis of malignant tumors [45].Jiang et al.validated that miR-29a-5p could down-regulate the abnormal TPX2 gene,thereby inhibiting the invasion and proliferation of endometrial cancer-derived cells and promoting apoptosis [46].We explored the potential regulatory relationships between KGs and miRNAs to identify miRNA biomarkers related to ESCC (Figure 6C).Consistent with previous studies,we found that hsa-miR-183-5p,hsa-miR-21-5p,hsa-miR-483-5p and hsa-miR-196a-5p had potential diagnostic and/or prognostic biomarker roles for ESCC [47-49].As far as we know,other dysregulated DEmiRNAs might also directly or indirectly affect the etiology and progression of ESCC,and more experiments were needed to support these hypotheses.

    Studies have proved that alisertib can inhibit gastrointestinal cancer and glioblastoma by inhibiting the expression of AURKA [50,51].Danusertib is also an effective pan-Aurora kinase inhibitor,which results in cell cycle arrest and programmed cell death in tumor cells[52].Dexrazoxane,an inhibitor of TOP2A,inhibits cancer progression by inducing DNA fragmentation and DNA damage in cancer cells[53].As far as we know,the above drugs are rarely used in the treatment of ESCC,so the data may provide new clues for targeted therapy of ESCC.We speculate that the remaining five KGs may be novel therapeutic targets for ESCC in the future,although no drug information related to them was obtained in this study.

    Conclusions

    In conclusion,we identified nine KGs(NDC80,BUB1,TOP2A,AURKA,AURKB,TTK,UBE2C,TPX2,and BUB1B)and some significant TFs and miRNAs from large-scale ESCC-related expression data,which may be promising biomarkers for the diagnosis,treatment,and prognosis of ESCC.We also screened 23 potential anti-neoplastic drugs.These data provide more useful diagnostic and therapeutic strategies for clinical research in ESCC.

    亚洲激情在线av| 两个人视频免费观看高清| 麻豆成人av在线观看| 成年人黄色毛片网站| 亚洲精品日韩av片在线观看| 久久人妻av系列| 欧美成人一区二区免费高清观看| 少妇人妻一区二区三区视频| 啦啦啦观看免费观看视频高清| 久久久久久久午夜电影| 男女下面进入的视频免费午夜| 国产黄片美女视频| 精品国内亚洲2022精品成人| 丁香欧美五月| 亚洲av不卡在线观看| 欧美zozozo另类| 在线观看舔阴道视频| 午夜精品在线福利| 国产成年人精品一区二区| 少妇高潮的动态图| 如何舔出高潮| 久久热精品热| 丰满人妻熟妇乱又伦精品不卡| 亚洲av日韩精品久久久久久密| 亚洲成人免费电影在线观看| 两人在一起打扑克的视频| 国产一区二区激情短视频| 国内精品美女久久久久久| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 天堂影院成人在线观看| 欧美性猛交黑人性爽| 国内精品美女久久久久久| 久久久色成人| 国产精品久久久久久精品电影| 日本熟妇午夜| 一进一出抽搐动态| 九九久久精品国产亚洲av麻豆| a级一级毛片免费在线观看| 中国美女看黄片| 久久人人精品亚洲av| 脱女人内裤的视频| 精品人妻熟女av久视频| 哪里可以看免费的av片| 永久网站在线| 久久伊人香网站| 天堂网av新在线| 亚洲国产精品成人综合色| 午夜福利欧美成人| 人妻制服诱惑在线中文字幕| 97超级碰碰碰精品色视频在线观看| 亚洲熟妇中文字幕五十中出| 一个人看视频在线观看www免费| 简卡轻食公司| 高清日韩中文字幕在线| 国产黄片美女视频| 亚洲欧美精品综合久久99| 亚洲av一区综合| 日韩中字成人| 欧美最新免费一区二区三区 | 欧美精品国产亚洲| 黄色日韩在线| 色av中文字幕| 亚洲欧美日韩卡通动漫| 村上凉子中文字幕在线| 内射极品少妇av片p| 91麻豆av在线| 亚洲狠狠婷婷综合久久图片| 亚洲成人精品中文字幕电影| 国产成+人综合+亚洲专区| 亚洲精品乱码久久久v下载方式| 亚洲欧美清纯卡通| 国产精品,欧美在线| 日本黄大片高清| av黄色大香蕉| 成年免费大片在线观看| 一个人免费在线观看电影| 国产不卡一卡二| 欧美日本亚洲视频在线播放| 精品人妻1区二区| 亚洲av成人不卡在线观看播放网| 人妻制服诱惑在线中文字幕| 欧美国产日韩亚洲一区| 亚洲人与动物交配视频| 国产男靠女视频免费网站| 国内精品一区二区在线观看| 久久精品国产99精品国产亚洲性色| 日本 欧美在线| 亚洲国产欧美人成| 99久久九九国产精品国产免费| 在线观看一区二区三区| 中文字幕av成人在线电影| 人妻久久中文字幕网| a级毛片a级免费在线| 国产免费一级a男人的天堂| 色尼玛亚洲综合影院| 亚洲国产精品成人综合色| 一级a爱片免费观看的视频| 熟妇人妻久久中文字幕3abv| 久久6这里有精品| 国产黄片美女视频| 淫妇啪啪啪对白视频| 免费电影在线观看免费观看| 日本五十路高清| 成人高潮视频无遮挡免费网站| 性欧美人与动物交配| 国产伦人伦偷精品视频| 91麻豆精品激情在线观看国产| 亚洲国产高清在线一区二区三| 亚洲无线观看免费| 亚洲美女搞黄在线观看 | 色播亚洲综合网| 午夜福利欧美成人| 欧美三级亚洲精品| 国产亚洲精品av在线| 久久精品国产亚洲av天美| 免费大片18禁| 欧美午夜高清在线| 97碰自拍视频| 国产精品1区2区在线观看.| 小蜜桃在线观看免费完整版高清| 亚洲人成伊人成综合网2020| 俄罗斯特黄特色一大片| 男人舔女人下体高潮全视频| 国产高清视频在线播放一区| 少妇高潮的动态图| 国内精品久久久久精免费| 桃色一区二区三区在线观看| 丝袜美腿在线中文| 88av欧美| 国产乱人视频| 亚洲经典国产精华液单 | 精品国内亚洲2022精品成人| 1024手机看黄色片| 久久久久久大精品| 久久国产乱子免费精品| 成人毛片a级毛片在线播放| 成人高潮视频无遮挡免费网站| 欧美国产日韩亚洲一区| 亚洲欧美清纯卡通| 一个人免费在线观看电影| 久久人人爽人人爽人人片va | 小蜜桃在线观看免费完整版高清| 国产麻豆成人av免费视频| 热99在线观看视频| 午夜a级毛片| 午夜a级毛片| 国产主播在线观看一区二区| 国产三级在线视频| 亚洲经典国产精华液单 | 成人国产综合亚洲| 成人欧美大片| 欧美日韩福利视频一区二区| 在线观看美女被高潮喷水网站 | 国产精品电影一区二区三区| 观看免费一级毛片| 91麻豆av在线| 99久国产av精品| 精品熟女少妇八av免费久了| 欧美日韩乱码在线| 噜噜噜噜噜久久久久久91| 变态另类丝袜制服| 国产一区二区亚洲精品在线观看| 波多野结衣高清作品| 久久精品影院6| 一二三四社区在线视频社区8| 成人av一区二区三区在线看| 我的老师免费观看完整版| 999久久久精品免费观看国产| 国内精品美女久久久久久| 婷婷精品国产亚洲av| 看十八女毛片水多多多| 成人国产综合亚洲| 757午夜福利合集在线观看| 久久久久免费精品人妻一区二区| 亚洲第一电影网av| 老司机午夜十八禁免费视频| 亚洲av熟女| 久久亚洲真实| 日本黄色视频三级网站网址| 久久精品国产亚洲av涩爱 | 国产 一区 欧美 日韩| 可以在线观看的亚洲视频| 禁无遮挡网站| 欧美成人性av电影在线观看| 天堂√8在线中文| 亚洲美女视频黄频| 欧美一区二区精品小视频在线| 伦理电影大哥的女人| 欧美最黄视频在线播放免费| 免费在线观看亚洲国产| 99在线视频只有这里精品首页| 久久久国产成人免费| 午夜老司机福利剧场| 狠狠狠狠99中文字幕| av女优亚洲男人天堂| 老司机午夜十八禁免费视频| 久久6这里有精品| 一区二区三区四区激情视频 | 亚洲男人的天堂狠狠| 国产欧美日韩一区二区精品| 少妇人妻精品综合一区二区 | 国产av在哪里看| 国产三级中文精品| 男女下面进入的视频免费午夜| 亚洲av.av天堂| 国产精品av视频在线免费观看| 色噜噜av男人的天堂激情| 日本与韩国留学比较| 久久国产乱子伦精品免费另类| av中文乱码字幕在线| 一级a爱片免费观看的视频| 在线观看一区二区三区| 美女黄网站色视频| 天堂av国产一区二区熟女人妻| 在线播放国产精品三级| 亚洲久久久久久中文字幕| 亚洲天堂国产精品一区在线| 国内精品久久久久精免费| 国产激情偷乱视频一区二区| 美女高潮的动态| 亚洲精品乱码久久久v下载方式| 国产av一区在线观看免费| 日本一本二区三区精品| 欧洲精品卡2卡3卡4卡5卡区| 99精品久久久久人妻精品| 桃红色精品国产亚洲av| 久99久视频精品免费| 成年免费大片在线观看| 如何舔出高潮| 欧美3d第一页| 国产亚洲欧美98| 国产av一区在线观看免费| 国内精品久久久久精免费| 亚洲不卡免费看| 国产av麻豆久久久久久久| 日韩欧美在线乱码| 亚洲人成网站在线播放欧美日韩| 动漫黄色视频在线观看| 1000部很黄的大片| 国产精品久久久久久亚洲av鲁大| 最好的美女福利视频网| 99国产极品粉嫩在线观看| 波多野结衣高清无吗| 午夜免费成人在线视频| 久久精品国产亚洲av香蕉五月| 精品一区二区免费观看| 国产成人福利小说| 国产精品精品国产色婷婷| 亚洲国产精品成人综合色| 在线观看一区二区三区| 久久这里只有精品中国| 欧美午夜高清在线| 国产精品av视频在线免费观看| 五月玫瑰六月丁香| 亚洲av电影在线进入| 99久久无色码亚洲精品果冻| 男人和女人高潮做爰伦理| 少妇人妻精品综合一区二区 | 级片在线观看| 国产精品亚洲av一区麻豆| 亚洲av不卡在线观看| 女人被狂操c到高潮| 国内久久婷婷六月综合欲色啪| 免费在线观看日本一区| 国产免费男女视频| av国产免费在线观看| 日韩中文字幕欧美一区二区| 悠悠久久av| 美女大奶头视频| 色综合亚洲欧美另类图片| 九色成人免费人妻av| 免费观看人在逋| 午夜福利视频1000在线观看| 丰满的人妻完整版| 精品午夜福利在线看| 亚洲熟妇中文字幕五十中出| 90打野战视频偷拍视频| 一区二区三区免费毛片| 国产在视频线在精品| 中文字幕人妻熟人妻熟丝袜美| 国产大屁股一区二区在线视频| 麻豆一二三区av精品| 人人妻人人看人人澡| 变态另类成人亚洲欧美熟女| av天堂在线播放| 88av欧美| 国产黄色小视频在线观看| 禁无遮挡网站| 夜夜夜夜夜久久久久| 尤物成人国产欧美一区二区三区| 免费在线观看日本一区| 国产乱人伦免费视频| 亚洲精品亚洲一区二区| 亚洲中文字幕一区二区三区有码在线看| 村上凉子中文字幕在线| 99国产综合亚洲精品| 婷婷丁香在线五月| 熟妇人妻久久中文字幕3abv| 无遮挡黄片免费观看| 最近最新免费中文字幕在线| 级片在线观看| 亚洲中文字幕日韩| 黄片小视频在线播放| 桃红色精品国产亚洲av| 精品久久久久久久久久久久久| 黄色配什么色好看| 又紧又爽又黄一区二区| 欧美一区二区精品小视频在线| 成人三级黄色视频| 亚洲熟妇中文字幕五十中出| av视频在线观看入口| 高清毛片免费观看视频网站| 欧美黄色淫秽网站| 亚洲久久久久久中文字幕| 亚洲第一电影网av| 成人欧美大片| h日本视频在线播放| 日韩欧美国产一区二区入口| 国产乱人视频| 亚洲av免费在线观看| 狠狠狠狠99中文字幕| 一进一出抽搐动态| 极品教师在线视频| 色综合婷婷激情| 女同久久另类99精品国产91| 精品久久久久久,| 欧美xxxx黑人xx丫x性爽| 久久国产精品人妻蜜桃| 国产精品国产高清国产av| 成年版毛片免费区| 国产综合懂色| 中文字幕人妻熟人妻熟丝袜美| 成人美女网站在线观看视频| 亚洲美女视频黄频| 日本 av在线| 少妇丰满av| 日韩成人在线观看一区二区三区| 天堂动漫精品| 亚洲aⅴ乱码一区二区在线播放| 毛片女人毛片| 91麻豆av在线| 少妇的逼水好多| 草草在线视频免费看| 草草在线视频免费看| 深爱激情五月婷婷| 久久性视频一级片| 亚洲人成网站在线播| 看黄色毛片网站| 成人毛片a级毛片在线播放| av欧美777| 他把我摸到了高潮在线观看| www.999成人在线观看| a级毛片a级免费在线| 听说在线观看完整版免费高清| 久久精品人妻少妇| 亚洲成人免费电影在线观看| 亚洲一区二区三区色噜噜| 国产免费一级a男人的天堂| 网址你懂的国产日韩在线| 成人一区二区视频在线观看| 黄色丝袜av网址大全| 老司机午夜福利在线观看视频| 免费看美女性在线毛片视频| 亚洲av第一区精品v没综合| 国产三级中文精品| 99热精品在线国产| 欧美绝顶高潮抽搐喷水| 在线观看免费视频日本深夜| 色5月婷婷丁香| 99久久久亚洲精品蜜臀av| 最好的美女福利视频网| 欧美区成人在线视频| 中文字幕久久专区| 国产探花在线观看一区二区| 中文字幕免费在线视频6| 午夜视频国产福利| 国语自产精品视频在线第100页| 久久人妻av系列| 亚洲av.av天堂| 欧美日韩国产亚洲二区| 三级国产精品欧美在线观看| 精品久久久久久久人妻蜜臀av| 级片在线观看| 两人在一起打扑克的视频| 色尼玛亚洲综合影院| 成人美女网站在线观看视频| 尤物成人国产欧美一区二区三区| 免费黄网站久久成人精品 | 人人妻人人澡欧美一区二区| 国产成人欧美在线观看| 日韩亚洲欧美综合| 久久久色成人| 1000部很黄的大片| 中文在线观看免费www的网站| 欧美激情久久久久久爽电影| 99精品在免费线老司机午夜| 免费观看精品视频网站| 久久国产精品影院| 成人美女网站在线观看视频| 欧美激情国产日韩精品一区| 亚洲专区国产一区二区| 俄罗斯特黄特色一大片| 亚洲成人免费电影在线观看| 欧美一级a爱片免费观看看| 亚洲av五月六月丁香网| 在线免费观看的www视频| 成年女人永久免费观看视频| 露出奶头的视频| 深夜精品福利| 亚洲人成电影免费在线| 国产精品久久电影中文字幕| 日韩精品中文字幕看吧| 九九久久精品国产亚洲av麻豆| 国产精品爽爽va在线观看网站| 欧美黑人巨大hd| 热99re8久久精品国产| 小蜜桃在线观看免费完整版高清| av在线蜜桃| 久久久久亚洲av毛片大全| 国产成人福利小说| or卡值多少钱| 嫁个100分男人电影在线观看| av女优亚洲男人天堂| 最近最新免费中文字幕在线| 在线免费观看不下载黄p国产 | 久久久久久久亚洲中文字幕 | h日本视频在线播放| 国产三级黄色录像| 国产精品久久久久久精品电影| bbb黄色大片| 国产午夜福利久久久久久| 国产av不卡久久| 很黄的视频免费| 12—13女人毛片做爰片一| 久久久久九九精品影院| 伊人久久精品亚洲午夜| 女生性感内裤真人,穿戴方法视频| 国产精品久久电影中文字幕| 无遮挡黄片免费观看| 国产中年淑女户外野战色| 国产极品精品免费视频能看的| av视频在线观看入口| 欧美日本视频| 亚洲18禁久久av| 国产久久久一区二区三区| 久久欧美精品欧美久久欧美| 一本久久中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 免费高清视频大片| 亚洲精品乱码久久久v下载方式| 无遮挡黄片免费观看| 99riav亚洲国产免费| 人妻制服诱惑在线中文字幕| 亚洲av日韩精品久久久久久密| 精品一区二区三区视频在线| 一卡2卡三卡四卡精品乱码亚洲| 久久精品人妻少妇| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 免费av毛片视频| 色视频www国产| 日韩精品青青久久久久久| 久久午夜亚洲精品久久| 精品福利观看| 男女床上黄色一级片免费看| 别揉我奶头 嗯啊视频| 国产免费一级a男人的天堂| 欧美成人一区二区免费高清观看| 少妇高潮的动态图| 在线观看美女被高潮喷水网站 | 中文字幕精品亚洲无线码一区| 中文亚洲av片在线观看爽| 欧美激情在线99| 99国产精品一区二区三区| 噜噜噜噜噜久久久久久91| 日日摸夜夜添夜夜添小说| 免费高清视频大片| www.熟女人妻精品国产| 18禁黄网站禁片免费观看直播| 亚洲人成伊人成综合网2020| 日韩欧美 国产精品| 最近在线观看免费完整版| 精品不卡国产一区二区三区| 亚洲黑人精品在线| 少妇的逼水好多| 九色成人免费人妻av| 精品一区二区三区av网在线观看| 一区福利在线观看| 亚洲欧美精品综合久久99| 精品久久久久久,| 亚洲一区高清亚洲精品| 成人av一区二区三区在线看| a在线观看视频网站| 亚洲色图av天堂| 狠狠狠狠99中文字幕| 91九色精品人成在线观看| 少妇的逼好多水| 精品久久久久久久人妻蜜臀av| 有码 亚洲区| 美女被艹到高潮喷水动态| 一进一出好大好爽视频| 日本三级黄在线观看| 国产不卡一卡二| 亚洲国产精品久久男人天堂| 亚洲美女视频黄频| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品日韩av片在线观看| 久久天躁狠狠躁夜夜2o2o| 国产精品人妻久久久久久| 国内毛片毛片毛片毛片毛片| 在线播放无遮挡| av在线天堂中文字幕| 特级一级黄色大片| 国产精品女同一区二区软件 | 久久精品国产亚洲av涩爱 | 51午夜福利影视在线观看| 国产激情偷乱视频一区二区| 亚洲精品粉嫩美女一区| 啦啦啦观看免费观看视频高清| 久久99热这里只有精品18| 免费看日本二区| 少妇人妻一区二区三区视频| 国产单亲对白刺激| 久久国产乱子免费精品| 国产成人欧美在线观看| 欧美在线黄色| 日本三级黄在线观看| www.色视频.com| 99久久成人亚洲精品观看| 91麻豆精品激情在线观看国产| 午夜久久久久精精品| 亚洲激情在线av| 亚洲国产精品久久男人天堂| 国产精品av视频在线免费观看| 尤物成人国产欧美一区二区三区| 久久精品国产清高在天天线| 国产大屁股一区二区在线视频| 丰满人妻一区二区三区视频av| 乱码一卡2卡4卡精品| 国产免费av片在线观看野外av| 黄色丝袜av网址大全| www.色视频.com| 精品人妻一区二区三区麻豆 | 免费观看的影片在线观看| 中文亚洲av片在线观看爽| 国产亚洲精品久久久com| 久久久久久九九精品二区国产| 欧美xxxx黑人xx丫x性爽| 国产视频一区二区在线看| 欧美zozozo另类| .国产精品久久| 免费大片18禁| netflix在线观看网站| 亚洲美女视频黄频| 麻豆成人午夜福利视频| 日韩免费av在线播放| 亚洲欧美激情综合另类| 97人妻精品一区二区三区麻豆| 日本黄色片子视频| 色在线成人网| 国产亚洲欧美98| 亚洲成人中文字幕在线播放| 内地一区二区视频在线| 香蕉av资源在线| 淫妇啪啪啪对白视频| 欧美一区二区国产精品久久精品| 中文资源天堂在线| av视频在线观看入口| 波野结衣二区三区在线| 真人一进一出gif抽搐免费| 久久热精品热| 亚洲成a人片在线一区二区| 美女大奶头视频| 亚洲av不卡在线观看| 久久国产乱子免费精品| 人妻久久中文字幕网| 国产蜜桃级精品一区二区三区| 久久人妻av系列| 精品久久国产蜜桃| 一本综合久久免费| 日韩av在线大香蕉| 熟女电影av网| 欧美乱妇无乱码| 国产午夜福利久久久久久| 亚洲人成电影免费在线| 日本五十路高清| 色av中文字幕| 高清日韩中文字幕在线| 99热精品在线国产| 乱人视频在线观看| 欧美激情国产日韩精品一区| 日日干狠狠操夜夜爽| 亚洲专区国产一区二区| 国产淫片久久久久久久久 | 国产精品精品国产色婷婷| 亚洲欧美日韩东京热| 久久久久国内视频| 久久香蕉精品热| www.www免费av| 我的老师免费观看完整版| 亚洲成a人片在线一区二区| 欧美bdsm另类| 观看免费一级毛片| a级毛片免费高清观看在线播放| 久久久久久久久久成人| 国模一区二区三区四区视频| 亚洲成人免费电影在线观看| 亚洲国产色片| 91午夜精品亚洲一区二区三区 | xxxwww97欧美| 国产一区二区在线观看日韩| 天堂影院成人在线观看| 免费电影在线观看免费观看| 中文字幕久久专区| 舔av片在线|