Identification of new key genes and their association with breast cancer occurrence and poor survival using in silico and in vitro methods

Breast cancer is one of the most prevalent types of cancer diagnosed globally and continues to have a significant impact on the global number of cancer deaths. Despite all efforts of epidemiological and experimental research, therapeutic concepts in cancer are still unsatisfactory. Gene expression datasets are widely used to discover the new biomarkers and molecular therapeutic targets in diseases. In the present study, we analyzed four datasets using R packages with accession number GSE29044, GSE42568, GSE89116, and GSE109169 retrieved from NCBI-GEO and differential expressed genes (DEGs) were identified. Protein–protein interaction (PPI) network was constructed to screen the key genes. Subsequently, the GO function and KEGG pathways were analyzed to determine the biological function of key genes. Expression profile of key genes was validated in MCF-7 and MDA-MB-231 human breast cancer cell lines using qRT-PCR. Overall expression level and stage wise expression pattern of key genes was determined by GEPIA. The bc-GenExMiner was used to compare expression level of genes among groups of patients with respect to age factor. OncoLnc was used to analyze the effect of expression levels of LAMA2, TIMP4, and TMTC1 on the survival of breast cancer patients. We identified nine key genes, of which COL11A1, MMP11, and COL10A1 were found up-regulated and PCOLCE2, LAMA2, TMTC1, ADAMTS5, TIMP4, and RSPO3 were found down-regulated. Similar expression pattern of seven among nine genes (except ADAMTS5 and RSPO3) was observed in MCF-7 and MDA-MB-231 cells. Further, we found that LAMA2, TMTC1, and TIMP4 were significantly expressed among different age groups of patients. LAMA2 and TIMP4 were found significantly associated and TMTC1 was found less correlated with breast cancer occurrence. We found that the expression level of LAMA2, TIMP4, and TMTC1 was abnormal in all TCGA tumors and significantly associated with poor survival.Indian Council of Medical Research | Ref. BMI/11/(35)/2020MICINN | Ref. RYC-2017-2289

In silico analyses for potential key genes associated with gastric cancer

Background Understanding hub genes involved in gastric cancer (GC) metastasis could lead to effective approaches to diagnose and treat cancer. In this study, we aim to identify the hub genes and investigate the underlying molecular mechanisms of GC. Methods To explore potential therapeutic targets for GC,three expression profiles (GSE54129, GSE33651, GSE81948) of the genes were extracted from the Gene Expression Omnibus (GEO) database. The GEO2R online tool was applied to screen out differentially expressed genes (DEGs) between GC and normal gastric samples. Database for Annotation, Visualization and Integrated Discovery was applied to perform Gene Ontology (GO) and KEGG pathway enrichment analysis. The protein-protein interaction (PPI) network of these DEGs was constructed using a STRING online software. The hub genes were identified by the CytoHubba plugin of Cytoscape software. Then, the prognostic value of these identified genes was verified by gastric cancer database derived from Kaplan-Meier plotter platform. Results A total of 85 overlapped upregulated genes and 44 downregulated genes were identified. The majority of the DEGs were enriched in extracellular matrix organization, endodermal cell differentiation, and endoderm formation. Moreover, five KEGG pathways were significantly enriched, including ECM-receptor interaction, amoebiasis, AGE-RAGE signaling pathway in diabetic complications, focal adhesion, protein digestion and absorption. By combining the results of PPI network and CytoHubba, a total of nine hub genes including COL1A1, THBS1, MMP2, CXCL8, FN1, TIMP1, SPARC, COL4A1, and ITGA5 were selected. The Kaplan-Meier plotter database confirmed that overexpression levels of these genes were associated with reduced overall survival, except for THBS1 and CXCL8. Conclusions Our study suggests that COL1A1, MMP2, FN1, TIMP1, SPARC, COL4A1, and ITGA5 may be potential biomarkers and therapeutic targets for GC. Further study is needed to assess the effect of THBS1 and CXCL8 on GC

Common Core Genes Play Vital Roles in Gastric Cancer With Different Stages

Background: Owing to complex molecular mechanisms in gastric cancer (GC) oncogenesis and progression, existing biomarkers and therapeutic targets could not significantly improve diagnosis and prognosis. This study aims to identify the key genes and signaling pathways related to GC oncogenesis and progression using bioinformatics and meta-analysis methods.Methods: Eligible microarray datasets were downloaded and integrated using the meta-analysis method. According to the tumor stage, GC gene chips were classified into three groups. Thereafter, the three groups’ differentially expressed genes (DEGs) were identified by comparing the gene data of the tumor groups with those of matched normal specimens. Enrichment analyses were conducted based on common DEGs among the three groups. Then protein–protein interaction (PPI) networks were constructed to identify relevant hub genes and subnetworks. The effects of significant DEGs and hub genes were verified and explored in other datasets. In addition, the analysis of mutated genes was also conducted using gene data from The Cancer Genome Atlas database.Results: After integration of six microarray datasets, 1,229 common DEGs consisting of 1,065 upregulated and 164 downregulated genes were identified. Alpha-2 collagen type I (COL1A2), tissue inhibitor matrix metalloproteinase 1 (TIMP1), thymus cell antigen 1 (THY1), and biglycan (BGN) were selected as significant DEGs throughout GC development. The low expression of ghrelin (GHRL) is associated with a high lymph node ratio (LNR) and poor survival outcomes. Thereafter, we constructed a PPI network of all identified DEGs and gained 39 subnetworks and the top 20 hub genes. Enrichment analyses were performed for common DEGs, the most related subnetwork, and the top 20 hub genes. We also selected 61 metabolic DEGs to construct PPI networks and acquired the relevant hub genes. Centrosomal protein 55 (CEP55) and POLR1A were identified as hub genes associated with survival outcomes.Conclusion: The DEGs, hub genes, and enrichment analysis for GC with different stages were comprehensively investigated, which contribute to exploring the new biomarkers and therapeutic targets

INTEGRATIVE ANALYSIS OF OMICS DATA IN ADULT GLIOMA AND OTHER TCGA CANCERS TO GUIDE PRECISION MEDICINE

Transcriptomic profiling and gene expression signatures have been widely applied as effective approaches for enhancing the molecular classification, diagnosis, prognosis or prediction of therapeutic response towards personalized therapy for cancer patients. Thanks to modern genome-wide profiling technology, scientists are able to build engines leveraging massive genomic variations and integrating with clinical data to identify “at risk” individuals for the sake of prevention, diagnosis and therapeutic interventions. In my graduate work for my Ph.D. thesis, I have investigated genomic sequencing data mining to comprehensively characterise molecular classifications and aberrant genomic events associated with clinical prognosis and treatment response, through applying high-dimensional omics genomic data to promote the understanding of gene signatures and somatic molecular alterations contributing to cancer progression and clinical outcomes. Following this motivation, my dissertation has been focused on the following three topics in translational genomics. 1) Characterization of transcriptomic plasticity and its association with the tumor microenvironment in glioblastoma (GBM). I have integrated transcriptomic, genomic, protein and clinical data to increase the accuracy of GBM classification, and identify the association between the GBM mesenchymal subtype and reduced tumorpurity, accompanied with increased presence of tumor-associated microglia. Then I have tackled the sole source of microglial as intrinsic tumor bulk but not their corresponding neurosphere cells through both transcriptional and protein level analysis using a panel of sphere-forming glioma cultures and their parent GBM samples.FurthermoreI have demonstrated my hypothesis through longitudinal analysis of paired primary and recurrent GBM samples that the phenotypic alterations of GBM subtypes are not due to intrinsic proneural-to-mesenchymal transition in tumor cells, rather it is intertwined with increased level of microglia upon disease recurrence. Collectively I have elucidated the critical role of tumor microenvironment (Microglia and macrophages from central nervous system) contributing to the intra-tumor heterogeneity and accurate classification of GBM patients based on transcriptomic profiling, which will not only significantly impact on clinical perspective but also pave the way for preclinical cancer research. 2) Identification of prognostic gene signatures that stratify adult diffuse glioma patientsharboring1p/19q co-deletions. I have compared multiple statistical methods and derived a gene signature significantly associated with survival by applying a machine learning algorithm. Then I have identified inflammatory response and acetylation activity that associated with malignant progression of 1p/19q co-deleted glioma. In addition, I showed this signature translates to other types of adult diffuse glioma, suggesting its universality in the pathobiology of other subset gliomas. My efforts on integrative data analysis of this highly curated data set usingoptimizedstatistical models will reflect the pending update to WHO classification system oftumorsin the central nervous system (CNS). 3) Comprehensive characterization of somatic fusion transcripts in Pan-Cancers. I have identified a panel of novel fusion transcripts across all of TCGA cancer types through transcriptomic profiling. Then I have predicted fusion proteins with kinase activity and hub function of pathway network based on the annotation of genetically mobile domains and functional domain architectures. I have evaluated a panel of in -frame gene fusions as potential driver mutations based on network fusion centrality hypothesis. I have also characterised the emerging complexity of genetic architecture in fusion transcripts through integrating genomic structure and somatic variants and delineating the distinct genomic patterns of fusion events across different cancer types. Overall my exploration of the pathogenetic impact and clinical relevance of candidate gene fusions have provided fundamental insights into the management of a subset of cancer patients by predicting the oncogenic signalling and specific drug targets encoded by these fusion genes. Taken together, the translational genomic research I have conducted during my Ph.D. study will shed new light on precision medicine and contribute to the cancer research community. The novel classification concept, gene signature and fusion transcripts I have identified will address several hotly debated issues in translational genomics, such as complex interactions between tumor bulks and their adjacent microenvironments, prognostic markers for clinical diagnostics and personalized therapy, distinct patterns of genomic structure alterations and oncogenic events in different cancer types, therefore facilitating our understanding of genomic alterations and moving us towards the development of precision medicine

Network-based genetic profiling reveals cellular pathway differences between follicular thyroid carcinoma and follicular thyroid adenoma

Molecular mechanisms underlying the pathogenesis and progression of malignant thyroid cancers, such as follicular thyroid carcinomas (FTCs), and how these differ from benign thyroid lesions, are poorly understood. In this study, we employed network-based integrative analyses of FTC and benign follicular thyroid adenoma (FTA) lesion transcriptomes to identify key genes and pathways that differ between them. We first analysed a microarray gene expression dataset (Gene Expression Omnibus GSE82208, n = 52) obtained from FTC and FTA tissues to identify differentially expressed genes (DEGs). Pathway analyses of these DEGs were then performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) resources to identify potentially important pathways, and protein-protein interactions (PPIs) were examined to identify pathway hub genes. Our data analysis identified 598 DEGs, 133 genes with higher and 465 genes with lower expression in FTCs. We identified four significant pathways (one carbon pool by folate, p53 signalling, progesterone-mediated oocyte maturation signalling, and cell cycle pathways) connected to DEGs with high FTC expression; eight pathways were connected to DEGs with lower relative FTC expression. Ten GO groups were significantly connected with FTC-high expression DEGs and 80 with low-FTC expression DEGs. PPI analysis then identified 12 potential hub genes based on degree and betweenness centrality; namely, TOP2A, JUN, EGFR, CDK1, FOS, CDKN3, EZH2, TYMS, PBK, CDH1, UBE2C, and CCNB2. Moreover, transcription factors (TFs) were identified that may underlie gene expression differences observed between FTC and FTA, including FOXC1, GATA2, YY1, FOXL1, E2F1, NFIC, SRF, TFAP2A, HINFP, and CREB1. We also identified microRNA (miRNAs) that may also affect transcript levels of DEGs; these included hsa-mir-335-5p, -26b-5p, -124-3p, -16-5p, -192-5p, -1-3p, -17-5p, -92a-3p, -215-5p, and -20a-5p. Thus, our study identified DEGs, molecular pathways, TFs, and miRNAs that reflect molecular mechanisms that differ between FTC and benign FTA. Given the general similarities of these lesions and common tissue origin, some of these differences may reflect malignant progression potential, and include useful candidate biomarkers for FTC and identifying factors important for FTC pathogenesis

Screening Driving Transcription Factors in the Processing of Gastric Cancer

Background. Construction of the transcriptional regulatory network can provide additional clues on the regulatory mechanisms and therapeutic applications in gastric cancer. Methods. Gene expression profiles of gastric cancer were downloaded from GEO database for integrated analysis. All of DEGs were analyzed by GO enrichment and KEGG pathway enrichment. Transcription factors were further identified and then a global transcriptional regulatory network was constructed. Results. By integrated analysis of the six eligible datasets (340 cases and 43 controls), a bunch of 2327 DEGs were identified, including 2100 upregulated and 227 downregulated DEGs. Functional enrichment analysis of DEGs showed that digestion was a significantly enriched GO term for biological process. Moreover, there were two important enriched KEGG pathways: cell cycle and homologous recombination. Furthermore, a total of 70 differentially expressed TFs were identified and the transcriptional regulatory network was constructed, which consisted of 566 TF-target interactions. The top ten TFs regulating most downstream target genes were BRCA1, ARID3A, EHF, SOX10, ZNF263, FOXL1, FEV, GATA3, FOXC1, and FOXD1. Most of them were involved in the carcinogenesis of gastric cancer. Conclusion. The transcriptional regulatory network can help researchers to further clarify the underlying regulatory mechanisms of gastric cancer tumorigenesis

Exploring novel prognostic biomarkers and biologic processes involved in NASH, cirrhosis and HCC based on survival analysis using systems biology approach

There is an unmet need to develop medications or drug combinations which can stop advancement of NASH to liver cirrhosis and HCC. Therefore, identifying key biomarkers based on overall survival and the exploring biological processes involved in the pathogenesis and progression of NASH toward cirrhosis and HCC to improve therapeutic interventions is necessary. The microarray dataset from the GPL13667 platform was downloaded from the Gene Expression Omnibus (GEO). The inclusion criteria for the DEGs included an adjusted p-value <0.05 and a log(2) fold change >1. In the first step, protein-protein interaction (PPI) network of differentially expressed genes (DEGs) in every three groups (NASH, Cirrhosis and HCC) was constructed using STRING online database. In the second step, the DEGs of each group were imported to Cytoscape software separately, and the genes with Degree >3 were selected. In the third step, the genes with Degree >3 were imported to Gephi software (0.9.2 version) and the genes with betweenness centrality >0 were selected. Venn diagram of these genes was depicted for the NASH, cirrhotic and HCC groups. According to venn diagram, 96 (NASH), 30 (cirrhotic) and 213 (HCC) genes were specifically upregulated (based on inclusion criteria). Among specific genes, 22 (NASH), 5 (cirrhotic) and 82 (HCC) genes were with poor overall survival. The overlap among the 3 groups (NASH, Cirrhosis and HCC) contained 4 upregulated genes HLA-F, HLA-DPA1, TPM1 and YWHAZ. From 4 upregulated genes, only YWHAZ gene was with poor overall survival. The present study detected new candidate genes and key biological processes in NASH, cirrhosis and HCC based on overall survival and using in silico analysis. Therefore, performing in vitro and in vivo researches to verify the results is necessary.Comment: 20 pages, 5 figures, 1 table

Identification of Key Biomarkers in Bladder Cancer: Evidence from a Bioinformatics Analysis

Bladder cancer (BCa) is one of the most common malignancies and has a relatively poor outcome worldwide. However, the molecular mechanisms and processes of BCa development and progression remain poorly understood. Therefore, the present study aimed to identify candidate genes in the carcinogenesis and progression of BCa. Five GEO datasets and TCGA-BLCA datasets were analyzed by statistical software R, FUNRICH, Cytoscape, and online instruments to identify differentially expressed genes (DEGs), to construct protein‒protein interaction networks (PPIs) and perform functional enrichment analysis and survival analyses. In total, we found 418 DEGs. We found 14 hub genes, and gene ontology (GO) analysis revealed DEG enrichment in networks and pathways related to cell cycle and proliferation, but also in cell movement, receptor signaling, and viral carcinogenesis. Compared with noncancerous tissues, TPM1, CRYAB, and CASQ2 were significantly downregulated in BCa, and the other hub genes were significant upregulated. Furthermore, MAD2L1 and CASQ2 potentially play a pivotal role in lymph nodal metastasis. CRYAB and CASQ2 were both significantly correlated with overall survival (OS) and disease-free survival (DFS). The present study highlights an up to now unrecognized possible role of CASQ2 in cancer (BCa). Furthermore, CRYAB has never been described in BCa, but our study suggests that it may also be a candidate biomarker in BCa