Functional Enrichment Pathway Analysis Research Articles

Exploring hub genes and crucial pathways linked to oxidative stress in bipolar disorder depressive episodes through bioinformatics analysis.

Bipolar disorder (BD) is a complex and serious psychiatric condition primarily characterized by bipolar depression, with the underlying genetic determinants yet to be elucidated. There is a substantial body of literature linking psychiatric disorders, including BD, to oxidative stress (OS). Consequently, this study aims to assess the relationship between BD and OS by identifying key hub genes implicated in OS pathways. We acquired gene microarray data from GSE5392 through the Gene Expression Omnibus (GEO). Our approach encompassed differential expression analysis, weighted gene co-expression network analysis (WGCNA), and Protein-Protein Interaction (PPI) Network analysis to pinpoint hub genes associated with BD. Subsequently, we utilized Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) to identify hub genes relevant to OS. To evaluate the diagnostic accuracy of these hub genes, we performed receiver operating characteristic curve (ROC) analysis on both GSE5388 and GSE5389 datasets. Furthermore, we conducted a study involving ten BD patients and ten healthy controls (HCs) who met the special criteria, assessing the expression levels of these hub genes in their peripheral blood mononuclear cells (PBMCs). We identified 411 down-regulated genes and 69 up-regulated genes for further scrutiny. Through WGCNA, we obtained 22 co-expression modules, with the sienna3 module displaying the strongest association with BD. By integrating differential analysis with genes linked to OS, we identified 44 common genes. Subsequent PPI Network and WGCNA analyses confirmed three hub genes as potential biomarkers for BD. Functional enrichment pathway analysis revealed their involvement in neuronal signal transduction, oxidative phosphorylation, and metabolic obstacle pathways. Using the Cytoscape plugin "ClueGo assay," we determined that a majority of these targets regulate neuronal synaptic plasticity. ROC curve analysis underscored the excellent diagnostic value of these three hub genes. Quantitative reverse transcription-PCR (RT-qPCR) results indicated significant changes in the expression of these hub genes in the PBMCs of BD patients compared to HCs. We identified three hub genes (TAC1, MAP2K1, and MAP2K4) in BD associated with OS, potentially influencing the diagnosis and treatment of BD. Based on the GEO database, our study provides novel insights into the relationship between BD and OS, offering promising therapeutic targets.

Beyond the Rhythm: In Silico Identification of Key Genes and Therapeutic Targets in Atrial Fibrillation.

Atrial fibrillation (AF) is a prevalent cardiac arrhythmia worldwide and is characterized by a high risk of thromboembolism, ischemic stroke, and fatality. The precise molecular mechanisms of AF pathogenesis remain unclear. The purpose of this study was to use bioinformatics tools to identify novel key genes in AF, provide deeper insights into the molecular pathogenesis of AF, and uncover potential therapeutic targets. Four publicly available raw RNA-Seq datasets obtained through the ENA Browser, as well as proteomic analysis results, both derived from atrial tissues, were used in this analysis. Differential gene expression analysis was performed and cross-validated with proteomics results to identify common genes/proteins between them. A functional enrichment pathway analysis was performed. Cross-validation analysis revealed five differentially expressed genes, namely FGL2, IGFBP5, NNMT, PLA2G2A, and TNC, in patients with AF compared with those with sinus rhythm (SR). These genes play crucial roles in various cardiovascular functions and may be part of the molecular signature of AF. Furthermore, functional enrichment analysis revealed several pathways related to the extracellular matrix, inflammation, and structural remodeling. This study highlighted five key genes that constitute promising candidates for further experimental exploration as biomarkers as well as therapeutic targets for AF.

Identifying circRNA-miRNA-mRNA Regulatory Networks in Chemotherapy-Induced Peripheral Neuropathy.

Chemotherapy-induced peripheral neuropathy (CIPN) is a frequent and severe side effect of first-line chemotherapeutic agents. The association between circular RNAs (circRNAs) and CIPN remains unclear. In this study, CIPN models were constructed with Taxol, while 134 differentially expressed circRNAs, 353 differentially expressed long non-coding RNAs, and 86 differentially expressed messenger RNAs (mRNAs) were identified utilizing RNA sequencing. CircRNA-targeted microRNAs (miRNAs) were predicted using miRanda, and miRNA-targeted mRNAs were predicted using TargetScan and miRDB. The intersection of sequencing and mRNA prediction results was selected to establish the circRNA-miRNA-mRNA networks, which include 15 circRNAs, 18 miRNAs, and 11 mRNAs. Functional enrichment pathway analyses and immune infiltration analyses revealed that differentially expressed mRNAs were enriched in the immune system, especially in T cells, monocytes, and macrophages. Cdh1, Satb2, Fas, P2ry2, and Zfhx2 were further identified as hub genes and validated by RT-qPCR, correlating with macrophages, plasmacytoid dendritic cells, and central memory CD4 T cells in CIPN. Additionally, we predicted the associated diseases, 36 potential transcription factors (TFs), and 30 putative drugs for hub genes using the DisGeNET, TRRUST, and DGIdb databases, respectively. Our results indicated the crucial role of circRNAs, and the immune microenvironment played in CIPN, providing novel insights for further research.

Comprehensive analysis of lncRNA-mRNA co-expression networks in HPV-driven cervical cancer reveals the pivotal function of LINC00511-PGK1 in tumorigenesis

BackgroundMounting evidence suggests that noncoding RNAs (lncRNAs) were involved in various human cancers. However, the role of these lncRNAs in HPV-driven cervical cancer (CC) has not been extensively studied. Considering that HR-HPV infections contribute to cervical carcinogenesis by regulating the expression of lncRNAs, miRNAs and mRNAs, we aim to systematically analyze lncRNAs and mRNAs expression profile to identify novel lncRNAs-mRNAs co-expression networks and explore their potential impact on tumorigenesis in HPV-driven CC. MethodsLncRNA/mRNA microarray technology was utilized to identify the differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) in HPV-16 and HPV-18 cervical carcinogenesis compared to normal cervical tissues. Venn diagram and weighted gene co-expression network analysis (WGCNA) were used to identify the hub DElncRNAs/DEmRNAs that were both significantly correlated with HPV-16 and HPV-18 CC patients. LncRNA-mRNA correlation analysis and functional enrichment pathway analysis were performed on these key DElncRNAs/DEmRNAs in HPV-16 and HPV-18 CC patients to explore their mutual mechanism in HPV-driven CC. A lncRNA-mRNA co-expression score (CES) model was established and validated by using the Cox regression method. Afterward, the clinicopathological characteristics were analyzed between CES-high and CES-low groups. In vitro, functional experiments were performed to evaluate the role of LINC00511 and PGK1 in cell proliferation, migration and invasion in CC cells. To understand whether LINC00511 play as an oncogenic role partially via modulating the expression of PGK1, rescue assays were used. ResultsWe identified 81 lncRNAs and 211 mRNAs that were commonly differentially expressed in HPV-16 and HPV-18 CC tissues compared to normal tissues. The results of lncRNA-mRNA correlation analysis and functional enrichment pathway analysis showed that the LINC00511-PGK1 co-expression network may make an important contribution to HPV-mediated tumorigenesis and be closely associated with metabolism-related mechanisms. Combined with clinical survival data, the prognostic lncRNA-mRNA co-expression score (CES) model based on LINC00511 and PGK1 could precisely predict patients’ overall survival (OS). CES-high patients had a worse prognosis than CES-low patients and the enriched pathways and potential targets of applicable drugs were explored in CES-high patients. In vitro experiments confirmed the oncogenic functions of LINC00511 and PGK1 in the progression of CC, and revealed that LINC00511 functions in an oncogenic role in CC cells partially via modulating the expression of PGK1. ConclusionsTogether, these data identify co-expression modules that provide valuable information to understand the pathogenesis of HPV-mediated tumorigenesis, which highlights the pivotal function of the LINC00511-PGK1 co-expression network in cervical carcinogenesis. Furthermore, our CES model has a reliable predicting ability that could stratify CC patients into low- and high-risk groups of poor survival. This study provides a bioinformatics method to screen prognostic biomarkers which leads to lncRNA-mRNA co-expression network identification and construction for patients’ survival prediction and potential drug applications in other cancers.

Differential expression of MicroRNAs in Alzheimer's disease: a systematic review and meta-analysis.

Alzheimer's disease (AD) results in progressive cognitive decline owing to the accumulation of amyloid plaques and hyperphosphorylated tau. MicroRNAs (miRNAs) have attracted attention as a putative diagnostic and therapeutic target for neurodegenerative diseases. However, existing meta-analyses on AD and its association with miRNAs have produced inconsistent results. The primary objective of this study is to evaluate the magnitude and consistency of differences in miRNA levels between AD patients, mild cognitive impairment (MCI) patients and healthy controls (HC). Articles investigating miRNA levels in blood, brain tissue, or cerebrospinal fluid (CSF) of AD and MCI patients versus HC were systematically searched in PubMed/Medline from inception to February 16th, 2021. Fixed- and random-effects meta-analyses were complemented with the I2 statistic to measure the heterogeneity, assessment of publication bias, sensitivity subgroup analyses (AD severity, brain region, post-mortem versus ante-mortem specimen for CSF and type of analysis used to quantify miRNA) and functional enrichment pathway analysis. Of the 1512 miRNAs included in 61 articles, 425 meta-analyses were performed on 334 miRNAs. Fifty-six miRNAs were significantly upregulated (n = 40) or downregulated (n = 16) in AD versus HC and all five miRNAs were significantly upregulated in MCI versus HC. Functional enrichment analysis confirmed that pathways related to apoptosis, immune response and inflammation were statistically enriched with upregulated pathways in participants with AD relative to HC. This study confirms that miRNAs' expression is altered in AD and MCI compared to HC. These findings open new diagnostic and therapeutic perspectives for this disorder.

Identification of four key prognostic genes and three potential drugs in human papillomavirus negative head and neck squamous cell carcinoma

BackgroundHead and neck squamous cell carcinoma (HNSCC) is a common tumor worldwide with poor prognosis. The pathogenesis of human papillomavirus (HPV)-positive and HPV-negative HNSCCs differs. However, few studies have considered the HPV status when identifying biomarkers for HNSCC. Thus, the identification of biomarkers for HPV-positive and HPV-negative HNSCCs is urgently needed.MethodsThree microarray datasets from Gene Expression Omnibus (GEO) were analyzed, and the differentially expressed genes (DEGs) were obtained. Then, functional enrichment pathway analysis was performed and protein–protein interaction (PPI) networks were constructed. The expression of hub genes at both the mRNA and protein level was determined in Oncomine, The Cancer Genome Atlas (TCGA) and the Human Protein Atlas (HPA). In addition, survival analysis of the patient stratified by HPV status and the expression levels of key genes were performed based on TCGA data. The role of AREG, STAG3, CAV1 and C19orf57 in cancer were analyzed through Gene set enrichment analysis (GSEA). The top ten small molecule drugs were identified and the therapeutic value of zonisamide, NVP-AUY922, PP-2 and fostamatinib was further evaluated in six HPV-negative HNSCC cell lines. Finally, the therapeutic value of NVP-AUY922 was tested in vivo based on three HPV-negative HNSCC models, and statistical analysis was performed.ResultsIn total, 47 DEGs were obtained, 11 of which were identified as hub genes. Biological process analysis indicated that the hub genes were associated with the G1/S transition of the mitotic cell cycle. Survival analysis uncovered that the prognostic value of AREG, STAG3, C19orf57 and CAV1 differed between HPV-positive and HPV-negative patients. Gene set enrichment analysis (GSEA) showed the role of AREG, STAG3 and CAV1 in dysregulated pathways of tumor. Ten small molecules were identified as potential drugs specifically for HPV-positive or HPV-negative patients; three—NVP-AUY922, fostamatinib and PP-2—greatly inhibited the proliferation of six HPV-negative HNSCC cell lines in vitro, and NVP-AUY922 inhibited three HPV-negative HNSCC xenografts in vivo.ConclusionsIn conclusion, AREG, STAG3, C19orf57 and CAV1 are key prognostic factors and potential therapeutic targets in HPV-negative HNSCC. NVP-AUY922, fostamatinib and PP-2 may be effective drugs for HPV-negative HNSCC.

Comprehensive analysis of genomic mutation signature and tumor mutation burden for prognosis of intrahepatic cholangiocarcinoma

BackgroundIntrahepatic cholangiocarcinoma (iCCA) is a highly lethal malignancy of the biliary tract. Analysis of somatic mutational profiling can reveal new prognostic markers and actionable treatment targets. In this study, we explored the utility of genomic mutation signature and tumor mutation burden (TMB) in predicting prognosis in iCCA patients.MethodsWhole-exome sequencing and corresponding clinical data were collected from the ICGC portal and cBioPortal database to detect the prognostic mutated genes and determine TMB values. To identify the hub prognostic mutant signature, we used Cox regression and Lasso feature selection. Mutation-related signature (MRS) was constructed using multivariate Cox regression. The predictive performances of MRS and TMB were assessed using Kaplan–Meier (KM) analysis and receiver operating characteristic (ROC). We performed a functional enrichment pathway analysis using gene set enrichment analysis (GSEA) for mutated genes. Based on the MRS, TMB, and the TNM stage, a nomogram was constructed to visualize prognosis in iCCA patients.ResultsThe mutation landscape illustrated distributions of mutation frequencies and types in iCCA, and generated a list of most frequently mutated genes (such as Tp53, KRAS, ARID1A, and IDH1). Thirty-two mutated genes associated with overall survival (OS) were identified in iCCA patients. We obtained a six-gene signature using the Lasso and Cox method. AUCs for the MRS in the prediction of 1-, 3-, and 5-year OS were 0.759, 0.732, and 0.728, respectively. Kaplan–Meier analysis showed a significant difference in prognosis for patients with iCCA having a high and low MRS score (P < 0.001). GSEA was used to show that several signaling pathways, including MAPK, PI3K-AKT, and proteoglycan, were involved in cancer. Conversely, survival analysis indicated that TMB was significantly associated with prognosis. GSEA indicated that samples with high MRS or TMB also showed an upregulated expression of pathways involved in tumor signaling and the immune response. Finally, the predictive nomogram (that included MRS, TMB, and the TNM stage) demonstrated satisfactory performance in predicting survival in patients with iCCA.ConclusionsMutation-related signature and TMB were associated with prognosis in patients with iCCA. Our study provides a valuable prognostic predictor for determining outcomes in patients with iCCA.

975. Why Sex Make a Difference in HIV Clinical Course? Bioinformatics Analysis of Differential Expressed Gene in Females and Males with HIV Disease

BackgroundHuman immunodeficiency virus (HIV) disease progression are different among genders, in which women usually progress to acquired immunodeficiency syndrome (AIDS) faster than men. The mechanisms resulting in the gender biases of HIV progression are unclear. We conducted a bioinformatics analysis of differentially expressed genes (DEGs) in women and men with HIV disease to understand the sex-based differences in HIV pathogenesis.MethodsWe obtained microarray data from the Gene Expression Omnibus (GEO) database using our pre-defined search strategy and analyzed data using the GEO2R platform. The t-test was done to compare DEGs between females and males with HIV diseases. The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was implemented to systematically extract biological features and processes of retrieving DEGs via gene ontology (GO) analysis. A Systemic search was performed to evaluate each DEG function and its possible association with HIV.ResultsOne gene expression profiling data were retrieved: GSE 140713, composed of 40 males and 10 females with HIV1 infected samples. A GEO2R analysis yielded 19 DEGs (Table 1). The GO analysis result was demonstrated in Tables 2 and 3. Following a systemic search, we found two DEGs, which have previous studies reported an association with HIV: DDX3X (20 studies) and PDS5 (1 study). We proposed DDX3X (t 5.3, p 0.0037) is responsible for gender inequalities of HIV progression because of: 1. DDX3X is needed in the HIV1 life cycle. 2. Several studies confirmed a positive correlation between DDX3X expression and HIV1 replication. 3. Our study found an up-regulated DDX3X expression in women corresponded to the fact that women progress to AIDS faster than men. 4. Our GO analysis showed female up-regulated genes were enriched in positive regulation of the gene expression pathway, which can be explained by DDX3X and its underlying mechanism.Table 1: DEGs in women and men with HIV1 disease Table 2: GO functional enrichment pathway analyses of overall retrieving DEGs Table 3: GO functional enrichment pathway analyses of down- and up-regulated clusters of DEGs ConclusionAberrant DDX3X expression may contribute to sex-based differences in HIV disease. Drugs modifying DDX3X gene expression will be beneficial in the treatment of HIV especially resolving the HIV drug resistance problem because current anti-HIV drugs target viral components posed the risk of viral mutation.Disclosures All Authors: No reported disclosures

Suppression of colorectal cancer cell growth by combined treatment of 6-gingerol and γ-tocotrienol via alteration of multiple signalling pathways.

Our previous study reported that combined treatment of γ-tocotrienol with 6-gingerol showed promising anticancer effects by synergistically inhibiting proliferation of human colorectal cancer cell lines. This study aimed to identify and elucidate molecular mechanisms involved in the suppression of SW837 colorectal cancer cells modulated by combined treatment of γ-tocotrienol and 6-gingerol. Total RNA from both untreated and treated cells was prepared for transcriptome analysis using RNA sequencing techniques. We performed high-throughput sequencing at approximately 30-60million coverage on both untreated and 6G + γT3-treated cells. The results showed that cancer-specific differential gene expression occurred and functional enrichment pathway analysis suggested that more than one pathway was modulated in 6G + γT3-treated cells. Combined treatment with 6G + γT3 augmented its chemotherapeutic effect by interfering with the cell cycle process, downregulating the Wnt signalling pathway and inducing apoptosis mainly through caspase-independent programmed cell death through mitochondrial dysfunction, activation of ER-UPR, disruption of DNA repair mechanisms and inactivation of the cell cycle process through the downregulation of main genes in proliferation such as FOXM1 and its downstream genes. The combined treatment exerted its cytotoxic effect through upregulation of genes in stress response activation and cytostatic effects demonstrated by downregulation of main regulator genes in the cell cycle. Selected genes involved in particular pathways including ATF6, DDIT3, GADD34, FOXM1, CDK1 and p21 displayed concordant patterns of gene expression between RNA sequencing and RT-qPCR. This study provides new insights into combined treatment with bioactive compounds not only in terms of its pleiotropic effects that enhance multiple pathways but also specific target genes that could be exploited for therapeutic purposes, especially in suppressing cancer cell growth.

CDH13 and HCRTR2 May Be Associated with Hypersomnia Symptom of Bipolar Depression: A Genome-Wide Functional Enrichment Pathway Analysis.

Although bipolar disorder is highly heritable, the identification of specific genetic variations is limited because of the complex traits underlying the disorder. We performed a genome-wide association study of bipolar disorder using a subphenotype that shows hypersomnia symptom during a major depressive episode. We investigated a total of 2,191 cases, 1,434 controls, and 703,012 single nucleotide polymorphisms (SNPs) in the merged samples obtained from the Translational Genomics Institute and the Genetic Association Information Network. The gene emerging as the most significant by statistical analysis was rs1553441 (odds ratio=0.4093; p=1.20×10-5; Permuted p=6.0×10-6). However, the 5×0-8 threshold for statistical significance required in a genome-wide association study was not achieved. The functional enrichment pathway analysis showed significant enrichments in the adhesion, development-related, synaptic transmission-related, and cell recognition-related pathways. For further evaluation, each gene of the enriched pathways was reviewed and matched with genes that were suggested to be associated with psychiatric disorders by previous genetic studies. We found that the cadherin 13 and hypocretin (orexin) receptor 2 genes may be involved in the hypersomnia symptom during a major depressive episode of bipolar disorder.