Drug-Gene Interaction Database Research Articles

Exploring potential therapeutic strategy for hepatocellular carcinoma and COVID-19 using bioinformatics analysis

BackgroundHepatocellular carcinoma (HCC) constitutes an important contributor to fatalities. Coronavirus disease 2019 (COVID-19) frequently presents with complications such as respiratory distress, systemic inflammatory responses, and damage to various organs. Several studies have investigated the relationship between COVID-19 and mortality in patients with liver cancer, but there are few research on the relationship between them. This study is to explore the correlation between the two diseases and drugs treating them.MethodsThe Gene Expression Omnibus (GEO) database provides gene datasets of COVID-19 patients and HCC patients. Through differential gene analysis and weighted gene co-expression network analysis, we determined 223 genes represented in HCC and COVID-19. We then used functional annotation, protein–protein interaction network construction, predictive model development and verification, prognostic value analysis, and miRNA–gene network construction. Besides, we created a drug–hub–gene network by predicting possible medications that interact with hub genes using the Drug–Gene Interaction Database (DGIdb). Ultimately, we applied immunohistochemistry to ascertain the hub genes expression.ResultsThis study revealed that eight core genes (RRM2, TPX2, DTL, CDT1, TYMS, CDCA5, CDC25C, and HJURP) co-existed in both HCC and COVID-19 and were differentially expressed in both HCC and normal tissues.CDC25C, RRM2, CDCA5, and HJURP had diagnostic value (AUC > 0.8) and prognostic value (adjusted P-value < 0.05). Genome enrichment analysis indicated that eight genes may function in liver cancer through engagement in the cell cycle, DNA replication, etc. In liver cancer samples, these genes were significantly and adversely associated with plasma cells while RRM2 was positively associated with neutrophil and NK cell activation and with dendritic cell resting. Using the miRNAnet database and DGIdb, 9 transcription factors, 7 miRNAs, and 51 drugs or molecular compounds were predicted to interact with the hub genes. Finally, RRM2 expression showed significant variation in clinical specimens, and analysis of the association of RRM2 with immunomodulators indicated that RRM2 was closely connected to MICB and CD276.ConclusionsOur study revealed several metabolic genes related to HCC and COVID-19. Moreover, potential drugs related to central genes were predicted. These findings may provide new ideas for treating COVID-19 and HCC.

Identification of essential genes and drug discovery in bladder cancer and inflammatory bowel disease via text mining and bioinformatics analysis.

Bladder cancer (BCa) is the most common malignancy of the urinary system. Inflammation is critical in the occurrence and development of BCa. The purpose of this study was to identify key genes and pathways of inflammatory bowel disease in BCa through text mining technology and bioinformatics technology, and to explore potential therapeutic drugs for BCa. Genes associated with BCa and Crohn's disease (CD) were detected using the text mining tool GenClip3, and analyzed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). A protein-protein interaction (PPI) network was constructed by STRING and visualized in Cytoscape, and modular analysis was performed using the Molecular Complex Detection plugin (MCODE). Finally, the genes clustered in the first two modules were selected as core genes, and the drug-gene interaction database was used to discover potential therapeutic drugs. We identified 796 genes shared by "Bladder cancer" and " Crohn's disease " by text mining. Gene function enrichment analysis yielded 18 enriched GO terms and the 6 most relevant KEGG pathways. A PPI network with 758 nodes and 4014 edges was constructed, and 20 gene modules were obtained using MCODE. We selected the top two gene clusters as core candidate genes. We found that 3 out of 55 selected core genes could be targeted by 26 existing drugs. The results indicated that CXCL12, FGF2 and FSCN1 are potential key genes involved in CD with BCa. Additionally, 26 drugs were identified as potential therapeutics for BCa treatment and management.

Harnessing transcriptomic signals for amyotrophic lateral sclerosis to identify novel drugs and enhance risk prediction

IntroductionAmyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease. This study integrates common genetic association results from the latest ALS genome-wide association study (GWAS) summary statistics with functional genomic annotations with the aim of providing mechanistic insights into ALS risk loci, inferring drug repurposing opportunities, and enhancing prediction of ALS risk and clinical characteristics. MethodsGenes associated with ALS were identified using GWAS summary statistic methodology including SuSiE SNP-based fine-mapping, and transcriptome- and proteome-wide association study (TWAS/PWAS) analyses. Using several approaches, gene associations were integrated with the DrugTargetor drug-gene interaction database to identify drugs that could be repurposed for the treatment of ALS. Furthermore, ALS gene associations from TWAS were combined with observed blood expression in two external ALS case-control datasets to calculate polytranscriptomic scores and evaluate their utility for prediction of ALS risk and clinical characteristics, including site of onset, age at onset, and survival. ResultsSNP-based fine-mapping, TWAS and PWAS identified 118 genes associated with ALS, with TWAS and PWAS providing novel mechanistic insights. Drug repurposing analyses identified six drugs significantly enriched for interactions with ALS associated genes, though directionality could not be determined. Additionally, drug class enrichment analysis showed gene signatures linked to calcium channel blockers may reduce ALS risk, whereas antiepileptic drugs may increase ALS risk. Across the two observed expression target samples, ALS polytranscriptomic scores significantly predicted ALS risk (R2 = 5.1 %; p-value = 3.2 × 10−27) and clinical characteristics. ConclusionsFunctionally-informed analyses of ALS GWAS summary statistics identified novel mechanistic insights into ALS aetiology, highlighted several therapeutic research avenues, and enabled statistically significant prediction of ALS risk.

Necroptosis in recurrent implantation failure: A bioinformatics analysis of key genes and therapeutic targets.

Recurrent implantation failure (RIF), a major issue in assisted reproductive technology (ART), may be influenced by necroptosis, a form of cell death linked to several diseases. This study was aimed at investigating the involvement of necroptosis in RIF. Using RNA-sequencing data from the Gene Expression Omnibus database, we identified differentially expressed necroptosis-related genes (DENRGs) in RIF patients compared with those in controls. Functional enrichment, protein-protein interaction (PPI) networks, and transcription factor (TF) regulatory networks were analyzed to identify key genes. Immune cell infiltration was analyzed using the single-sample gene set enrichment analysis (ssGSEA) algorithm. Finally, potential therapeutic drugs targeting key genes were explored using a Drug Gene Interaction Database. In total, 20 DENRGs associated with RIF were identified, with a focus on 6 key genes (MLKL, FASLG, XIAP, CASP1, BIRC3, and TLR3) implicated in necroptosis and immune processes. These genes were used to develop a predictive model for RIF, which was validated in 2 datasets. The model and TF network analysis underscored the importance of TLR3. Immune infiltration analysis showed reduced levels of 16 immune cells in RIF patients, highlighting immune system alterations. Several drugs targeting CASP1, such as nivocasan and emricasan, were identified as potential treatments. The study sheds light on the role of necroptosis in RIF, identifying key genes and immune alterations that could serve as biomarkers and therapeutic targets. These findings pave the way for future experimental research and clinical applications targeting necroptosis in RIF treatment.

Integrative Bioinformatics Analysis for Targeting Hub Genes in Hepatocellular Carcinoma Treatment

Background: The damage in the liver and hepatocytes is where the primary liver cancer begins, and this is referred to as Hepatocellular Carcinoma (HCC). One of the best methods for detecting changes in gene expression of hepatocellular carcinoma is through bioinformatics approaches. Objective: This study aimed to identify potential drug target(s) hubs mediating HCC progression using computational approaches through gene expression and protein-protein interaction datasets. Methodology: Four datasets related to HCC were acquired from the GEO database, and Differentially Expressed Genes (DEGs) were identified. Using Evenn, the common genes were chosen. Using the Fun Rich tool, functional associations among the genes were identified. Further, protein-protein interaction networks were predicted using STRING, and hub genes were identified using Cytoscape. The selected hub genes were subjected to GEPIA and Shiny GO analysis for survival analysis and functional enrichment studies for the identified hub genes. The up-regulating genes were further studied for immunohistopathological studies using HPA to identify gene/protein expression in normal vs HCC conditions. Drug Bank and Drug Gene Interaction Database were employed to find the reported drug status and targets. Finally, STITCH was performed to identify the functional association between the drugs and the identified hub genes. Results: The GEO2R analysis for the considered datasets identified 735 upregulating and 284 downregulating DEGs. Functional gene associations were identified through the Fun Rich tool. Further, PPIN network analysis was performed using STRING. A comparative study was carried out between the experimental evidence and the other seven data evidence in STRING, revealing that most proteins in the network were involved in protein-protein interactions. Further, through Cytoscape plugins, the ranking of the genes was analyzed, and densely connected regions were identified, resulting in the selection of the top 20 hub genes involved in HCC pathogenesis. The identified hub genes were: KIF2C, CDK1, TPX2, CEP55, MELK, TTK, BUB1, NCAPG, ASPM, KIF11, CCNA2, HMMR, BUB1B, TOP2A, CENPF, KIF20A, NUSAP1, DLGAP5, PBK, and CCNB2. Further, GEPIA and Shiny GO analyses provided insights into survival ratios and functional enrichment studied for the hub genes. The HPA database studies further found that upregulating genes were involved in changes in protein expression in Normal vs HCC tissues. These findings indicated that hub genes were certainly involved in the progression of HCC. STITCH database studies uncovered that existing drug molecules, including sorafenib, regorafenib, cabozantinib, and lenvatinib, could be used as leads to identify novel drugs, and identified hub genes could also be considered as potential and promising drug targets as they are involved in the gene-chemical interaction networks. Conclusion: The present study involved various integrated bioinformatics approaches, analyzing gene expression and protein-protein interaction datasets, resulting in the identification of 20 topranked hubs involved in the progression of HCC. They are KIF2C, CDK1, TPX2, CEP55, MELK, TTK, BUB1, NCAPG, ASPM, KIF11, CCNA2, HMMR, BUB1B, TOP2A, CENPF, KIF20A, NUSAP1, DLGAP5, PBK, and CCNB2. Gene-chemical interaction network studies uncovered that existing drug molecules, including sorafenib, regorafenib, cabozantinib, and lenvatinib, can be used as leads to identify novel drugs, and the identified hub genes can be promising drug targets. The current study underscores the significance of targeting these hub genes and utilizing existing molecules to generate new molecules to combat liver cancer effectively and can be further explored in terms of drug discovery research to develop treatments for HCC.

Exploration of the Shared Gene Signatures and Molecular Mechanisms between Chronic Bronchitis and Antineutrophil Cytoplasmic Antibody-associated Glomerulonephritis: Evidence from Transcriptome Data

Background: Chronic Bronchitis (CB) is a recurrent and persistent pulmonary inflammation disease. Growing evidence suggests an association between CB and Anti-neutrophil Cytoplasmic Antibody-associated Glomerulonephritis (ANCA-GN). However, the precise mechanisms underlying their association remain unclear. Aims: The purpose of this study was to further explore the molecular mechanism of the occurrence of chronic bronchitis (CB) associated with anti-neutrophil cytoplasmic antibody-associated glomerulonephritis (ANCA- GN). Objective: Our study aimed to investigate the potential shared pathogenesis of CB-associated ANCA-GN. Methods: Datasets of ANCA (GSE108113 and GSE104948) and CB (GSE151052 and GSE162635) were obtained from the Gene Expression Omnibus (GEO) datasets. Firstly, GSE108113 and GSE151052 were analyzed to identify common differentially expressed genes (DEGs) by Limma package. Based on common DEGs, protein-protein interaction (PPI) network and functional enrichment analyses, including GO, KEGG, and GSEA, were performed. Then, hub genes were identified by degree algorithm and validated in GSE104948 and GSE162635. Further PPI network and functional enrichment analyses were performed on hub genes. Additionally, a competitive ceRNA network was constructed through miRanda and spongeScan. Transcription factors (TFs) were predicted and verified using the TRRUST database. Furthermore, the CIBERSORT algorithm was employed to explore immune cell infiltration. The Drug Gene Interaction Database (DGIDB) was utilized to predict small-molecular compounds of CB and ANCA-GN. Result: A total of 963 DEGs were identified in the integrated CB dataset, and 610 DEGs were identified in the integrated ANCA-GN dataset. Totally, we identified 22 common DEGs, of which 10 hub genes (LYZ, IRF1, PIK3CG, IL2RG, NT5E, ARG2, HBEGF, NFATC2, ALPL, and FKBP5) were primarily involved in inflammation and immune responses. Focusing on hub genes, we constructed a ceRNA network composed of 323 miRNAs and 348 lncRNAs. Additionally, five TFs (SP1, RELA, NFKB1, HIF1A, and SP3) were identified to regulate the hub genes. Furthermore, immune cell infiltration results revealed immunoregulation in CB and ANCA-GN. Finally, some small-molecular compounds (Daclizumab, Aldesleukin, and NT5E) were predicted to predominantly regulate inflammation and immunity, especially IL-2. Conclusion: Our study explores the inflammatory-immune pathways underlying CB-associated ANCA-GN and emphasizes the importance of NETs and lymphocyte differentiation, providing novel insights into the shared pathogenesis and therapeutic targets.

Potential Mechanism and Therapeutic Targets of Polymyositis.

To investigate the variability in the expression profile of genes associated with polymyositis (PM), explore the potential molecular mechanisms underlying PM, and predict novel targets for intervention. Descriptive study. Place and Duration of the Study: Department of Rheumatology, Taizhou Municipal Hospital, Taizhou, China, from August to November 2023. Three microarray datasets (GSE3112, GSE39454, and GSE128470) were extracted from the gene expression omnibus (GEO). The analysis of this research involved identifying the differentially expressed genes (DEGs) in PM compared to normal samples. Enrichment analysis, gene-microRNA, gene-transcription factor (TF), and protein-protein interaction (PPI) network studies were conducted to identify hub genes and relevant pathways. Additionally, the drug-gene interaction database (DGIdb) was used to predict therapeutic medications. Eighty-eight DEGs were identified. The enrichment analysis results highlighted the significant involvement of downregulated DEGs in antigen processing and presentation. Based on the PPI networks, seven hub genes with high connectivity degrees were selected including a cluster of differentiation 74 (CD74), human leukocyte antigen (HLA)-DPA1, HLA-B, guanylate-binding protein 1 (GBP1), recombinant 2', 5'-oligoadenylate synthetase 1 (OAS1), HLA-C, and HLA-E. This research screened-out core genes, projected prospective therapeutic medications, discovered DEGs between PM and normal samples, and offered fresh perspectives for additional research into the possible mechanism and therapeutic targets of PM. Polymyositis, DEGs, Hub genes, Bioinformatics, Potential therapeutic agents.

Comprehensive analysis of cuproptosis-related genes involved in immune infiltration and their use in the diagnosis of hepatic ischemia‒reperfusion injury: an experimental Study.

Hepatic ischemia reperfusion injury (HIRI) is a common injury not only during liver transplantation but also during major hepatic surgery. HIRI causes severe complications and affects the prognosis and survival of patients. Cuproptosis, a newly identified form of cell death, plays an important role in a variety of illnesses. However, its role in HIRI remains unknown. The GSE151648 dataset was mined from the Gene Expression Omnibus (GEO) database, and differences were analyzed for intersections. Based on the differentially expressed genes (DEGs), functional annotation, differentially expressed cuproptosis-related genes (DE-CRGs) identification and lasso logistic regression were conducted. Correlation analysis of DE-CRGs and immune infiltration was further conducted, and DE-CRGs were applied to construct an HIRI diagnostic model. The hierarchical clustering method was used to classify the specimens of HIRI, and functional annotation was conducted to verify the accuracy of these DE-CRGs in predicting HIRI progression. The GSE14951 microarray dataset and GSE171539 single-cell sequencing dataset were chosen as validation datasets. At the same time, the significance of DE-CRGs was verified using a mouse model of HIRI with cuproptosis inhibitors and inducers. Finally, a network of transcription-factor-DE-CRGs and miRNA-DE-CRGs was constructed to reveal the regulation mechanisms. And potential drugs for DE-CRGs were predicted using Drug Gene Interaction Database (DGIdb). Overall, 2390 DEGs and 19 DE-CRGs were identified. Through machine learning algorithms, 8 featured DE-CRGs (GNL3, ALAS1, TSC22D2, KLF5, GTF2B, DNTTIP2, SLFN11 and HNRNPU) were screened, and 2 cuproptosis-related subclusters were defined. Based on the 8 DE-CRGs obtained from the HIRI model (AUC=0.97), the nomogram model demonstrated accuracy in predicting HIRI. Eight DE-CRGs were highly expressed in HIRI samples and were negatively related to immune cell infiltration. A higher level of immune infiltration and expression of CRG group B was found in the HIRI population. Differences in cell death and immune regulation were found between the 2 groups. The diagnostic value of the 8 DE-CRGs was confirmed in the validation of two datasets. The identification of 7 DE-CRGs (SLFN11 excluded) by HIRI animal model experiments was also confirmed. Using hTFtarget, miRWalk and DGIDB database, we predicted that 17 transcription factors, 192 miRNAs and 10 drugs might interact with the DE-CRGs. This study shows that cuproptosis may occur in HIRI and is correlated with immune infiltration. Additionally, a cuproptosis-related predictive model was constructed for studying the causes of HIRI and developing targeted treatment options for HIRI.

Genetic associations in ankylosing spondylitis: circulating proteins as drug targets and biomarkers.

Ankylosing spondylitis (AS) is a complex condition with a significant genetic component. This study explored circulating proteins as potential genetic drug targets or biomarkers to prevent AS, addressing the need for innovative and safe treatments. We analyzed extensive data from protein quantitative trait loci (pQTLs) with up to 1,949 instrumental variables (IVs) and selected the top single-nucleotide polymorphism (SNP) associated with AS risk. Utilizing a two-sample Mendelian randomization (MR) approach, we assessed the causal relationships between identified proteins and AS risk. Colocalization analysis, functional enrichment, and construction of protein-protein interaction networks further supported these findings. We utilized phenome-wide MR (phenMR) analysis for broader validation and repurposing of drugs targeting these proteins. The Drug-Gene Interaction database (DGIdb) was employed to corroborate drug associations with potential therapeutic targets. Additionally, molecular docking (MD) techniques were applied to evaluate the interaction between target protein and four potential AS drugs identified from the DGIdb. Our analysis identified 1,654 plasma proteins linked to AS, with 868 up-regulated and 786 down-regulated. 18 proteins (AGER, AIF1, ATF6B, C4A, CFB, CLIC1, COL11A2, ERAP1, HLA-DQA2, HSPA1L, IL23R, LILRB3, MAPK14, MICA, MICB, MPIG6B, TNXB, and VARS1) that show promise as therapeutic targets for AS or biomarkers, especially MAPK14, supported by evidence of colocalization. PhenMR analysis linked these proteins to AS and other diseases, while DGIdb analysis identified potential drugs related to MAPK14. MD analysis indicated strong binding affinities between MAPK14 and four potential AS drugs, suggesting effective target-drug interactions. This study underscores the utility of MR analysis in AS research for identifying biomarkers and therapeutic drug targets. The involvement of Th17 cell differentiation-related proteins in AS pathogenesis is particularly notable. Clinical validation and further investigation are essential for future applications.

Identifying potential drug targets for varicose veins through integration of GWAS and eQTL summary data.

Varicose veins (VV) are a common chronic venous disease that is influenced by multiple factors. It affects the quality of life of patients and imposes a huge economic burden on the healthcare system. This study aimed to use integrated analysis methods, including Mendelian randomization analysis, to identify potential pathogenic genes and drug targets for VV treatment. This study conducted Summary-data-based Mendelian Randomization (SMR) analysis and colocalization analysis on data collected from genome-wide association studies and cis-expression quantitative trait loci databases. Only genes with PP.H4 > 0.7 in colocalization were chosen from the significant SMR results. After the above analysis, we screened 12 genes and performed Mendelian Randomization (MR) analysis on them. After sensitivity analysis, we identified four genes with potential causal relationships with VV. Finally, we used transcriptome-wide association studies and The Drug-Gene Interaction Database data to identify and screen the remaining genes and identified four drug targets for the treatment of VV. We identified four genes significantly associated with VV, namely, KRTAP5-AS1 [Odds ratio (OR) = 1.08, 95% Confidence interval (CI): 1.05-1.11, p = 1.42e-10] and PLEKHA5 (OR = 1.13, 95% CI: 1.06-1.20, p = 6.90e-5), CBWD1 (OR = 1.05, 95% CI: 1.01-1.11, p = 1.42e-2) and CRIM1 (OR = 0.87, 95% CI: 0.81-0.95, p = 3.67e-3). Increased expression of three genes, namely, KRTAP5-AS1, PLEKHA5, and CBWD1, was associated with increased risk of the disease, and increased expression of CRIM1 was associated with decreased risk of the disease. These four genes could be targeted for VV therapy. We identified four potential causal proteins for varicose veins with MR. A comprehensive analysis indicated that KRTAP5-AS1, PLEKHA5, CBWD1, and CRIM1 might be potential drug targets for varicose veins.

Functional Investigation and Two-sample Mendelian Randomization Study of Inguinal Hernia Hub Genes Obtained by Bioinformatics Analysis.

Inguinal hernia in adults is a common and frequent disease in surgery, prone to occur in the elderly or in those with a weak abdominal wall. Despite its prevalence, Molecular mechanisms underlying inguinal hernia formation are unclear. This study aims to identify potential gene markers for inguinal hernia and available drugs. Pubmed2Ensembl text mining was used to identify genes related to "inguinal hernia" keywords. The GeneCodis system was used to specify GO biological process terms and KEGG pathways defined in the Kyoto Encyclopedia of Genes and Genomes (KEGG). The STRING tool was used to construct protein-protein interaction networks, which were then visualized using Cytoscape.CytoHubba and Molecular Complex Detection were utilized to analyze the module (MCODE). A GO and KEGG analysis of gene modules was conducted using the DAVID platform database. Hub genes are those that are concentrated in prominent modules. The druggene interaction database was also used to identify potential drugs for inguinal hernia patients based on their interactions between the hub genes. Finally, a Mendelian randomization study was conducted based on genome-wide association studies to determine whether hub genes cause inguinal hernias. The identification of 96 genes associated with inguinal hernia was carried out using text mining techniques. It was constructed using PPI networks with 80 nodes and 476 edges, and the sequence of the genes was performed using CytoHubba. MCODE analysis identified three gene modules. Three modules contain 37 genes clustered as hub candidate genes associated with inguinal hernia patients. The PI3K-Akt, MAPK, AGE-RAGE, and HIF-1 pathways were found to be enriched in signaling pathways. Sixteen of the 37 genes were found to be targetable by 30 existing drugs. The relationship between hub genes and inguinal hernia was examined using Mendelian randomization. The research revealed nine genes that may be connected with inguinal hernia, such as POMC, CD40LG, TFRC, VWF, LOX, IGF2, BRCA1, TNF, and HGF in the plasma. By inverse variance weighting, ALB was associated with an increased risk of inguinal hernia with an OR of 1.203 (OR [95%] = 1,04 [1.012 to 1.089], p = 0.008). We identified potential hub genes for inguinal hernia, predicted potential drugs for inguinal hernia, and reverse-validated potential genes by Mendelian randomization. This may provide further insights into asymptomatic pre-diagnostic methods and contribute to studies to understand the molecular mechanisms of risk genes associated with inguinal hernia.

Network-based identification of key proteins and repositioning of drugs for non-small cell lung cancer.

NSCLC is a lethal cancer that is highly prevalent and accounts for 85% of cases of lung cancer. Conventional cancer treatments, such as chemotherapy and radiation, frequently exhibit limited efficacy and notable adverse reactions. Therefore, a drug repurposing method is proposed for effective NSCLC treatment. This study aims to evaluate candidate drugs that are effective for NSCLC at the clinical level using a systems biology and network analysis approach. Differentially expressed genes in transcriptomics data were identified using the systems biology and network analysis approaches. A network of gene co-expression was developed with the aim of detecting two modules of gene co-expression. Following that, the Drug-Gene Interaction Database was used to find possible drugs that target important genes within two gene co-expression modules linked to non-small cell lung cancer (NSCLC). The use of Cytoscape facilitated the creation of a drug-gene interaction network. Finally, gene set enrichment analysis was done to validate candidate drugs. Unlike previous research on repositioning drugs for NSCLC, which uses a gene co-expression network, this project is the first to research both gene co-expression and co-occurrence networks. And the co-occurrence network also accounts for differentially expressed genes in cancer cells and their adjacent normal cells. For effective management of non-small cell lung cancer (NSCLC), drugs that show higher gene regulation and gene affinity within the drug-gene interaction network are thought to be important. According to the discourse, NSCLC genes have a lot of control over medicines like vincristine, fluorouracil, methotrexate, clotrimazole, etoposide, tamoxifen, sorafenib, doxorubicin, and pazopanib. Hence, there is a possibility of repurposing these drugs for the treatment of non-small-cell lung cancer.

Identification of potential cell death-related biomarkers for diagnosis and treatment of osteoporosis

BackgroundThis study aimed to identify potential biomarkers for the diagnosis and treatment of osteoporosis (OP).MethodsData sets were downloaded from the Gene Expression Omnibus database, and differentially programmed cell death-related genes were screened. Functional analyses were performed to predict the biological processes associated with these genes. Least absolute shrinkage and selection operator (LASSO), support vector machine (SVM), and random forest (RF) machine learning algorithms were used to screen for characteristic genes, and receiver operating characteristics were used to evaluate the diagnosis of disease characteristic gene values. Gene set enrichment analysis (GSEA) and single-sample GSEA were conducted to analyze the correlation between characteristic genes and immune infiltrates. Cytoscape and the Drug Gene Interaction Database (DGIdb) were used to construct the mitochondrial RNA-mRNA-transcription factor network and explore small-molecule drugs. Reverse transcription real-time quantitative PCR (RT-qPCR) analysis was performed to evaluate the expression of biomarker genes in clinical samples.ResultsIn total, 25 differential cell death genes were identified. Among these, two genes were screened using the LASSO, SVM, and RF algorithms as characteristic genes, including BRSK2 and VPS35. In GSE56815, the area under the receiver operating characteristic curve of BRSK2 was 0.761 and that of VPS35 was 0.789. In addition, immune cell infiltration analysis showed that BRSK2 positively correlated with CD56dim natural killer cells and negatively correlated with central memory CD4 + T cells. Based on the data from DGIdb, hesperadin was associated with BRSK2, and melagatran was associated with VPS35. BRSK2 and VPS35 were expectably upregulated in OP group compared with controls (all p < 0.05).ConclusionsBRSK2 and VPS35 may be important diagnostic biomarkers of OP.

Abstract 3482: Single-cell transcriptome analysis reveals that PLK4 signature genes are highly associated with cancer proliferation and hotspot mutations in non-small cell lung cancer (NSCLC)

Abstract Background: Polo-like kinase 4 (PLK4) is a member of the polo-like kinase family known for its precise regulation during the cell cycle and its key role in centriole duplication, which is essential for proper cell division. Dysregulation of PLK4 is associated with abnormal centriole number, genomic instability and cancer development. In our previous study, we proposed PLK4 as a potential therapeutic target and biomarker for LUAD patients with TP53 mutations by bulk RNA-seq analysis and confirmed its association with drug resistance. In this study, we analyzed the expression of PLK4 and its association with cell cycle at the cellular level by scRNA-seq in NSCLC data. Methods: In a previous study, we analyzed NSCLC data from The Cancer Genome Atlas (TCGA) database using the UCSC Xena data portal. In this study, we analyzed E-MTAB-6653 dataset from the single cell expression atlas(https://www.ebi.ac.uk/gxa/sc/home). Results: In the previous study, we demonstrated that PLK4 expression is higher in primary tumors compared to normal tissues in NSCLC. Higher PLK4 expression in patients is associated with an unfavorable prognosis. Using genes highly correlated with PLK4, we confirmed a significant association between PLK4 and TP53 mutation by statistical analysis. Using the Drug Gene Interaction database, we observed an association between the drug resistance-related genes PBK and BIRC5 with PLK4, and these genes were also correlated with TP53 mutation. In this study, we investigated the association between PLK4 and the cell cycle using single cell RNA-seq. First, we performed data pre-processing to selectively isolate NSCLC cells. Using a set of 87 PLK4 signature genes obtained from previous research, we stratified the samples into PLK4 high and low groups based on the expression of these genes. We then analyzed cell cycle related genes in both groups and observed a higher expression of G2/M phase related genes in the PLK4 high group compared to the PLK4 low group. We also examined the expression of PLK4 signature genes based on pathological differences and found higher expression in tumor tissue compared to normal tissue. There was also an increasing expression gradient from the tumor margin to the tumor core. Finally, using Gene Ontology analysis of both groups, we confirmed the upregulation of genes associated with transcription and translation, as indicated by the GO terms obtained. Conclusion: In a previous study, the association between PLK4 expression, TP53 mutation and drug resistance was identified by bulk RNA-seq. In the current study, we used single cell RNA-seq to verify the correlation between PLK4 and cell cycle in NSCLC. In addition, we observed that the expression of PLK4 signature genes increased towards the tumor core based on pathological differences. Citation Format: Youngtaek Kim, Joon Yeon Hwang, YoungJoon Park, Dong Kwon Kim, Kwangmin Na, Seul Lee, Sujeong Baek, Seong-san Kang, Seung Min Yang, Mi Hyun Kim, Heekyung Han, Chai Young Lee, Yu Jin Han, Min Hee Hong, Jii Bum Lee, Sun Min Lim, Kyoung-Ho Pyo, Byoung Chul Cho. Single-cell transcriptome analysis reveals that PLK4 signature genes are highly associated with cancer proliferation and hotspot mutations in non-small cell lung cancer (NSCLC) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 3482.

Identification of a novel glioblastoma multiforme molecular subtype with poor prognosis and high immune infiltration based on oxidative stress-related genes.

Glioblastoma multiforme (GBM) is a highly malignant primary brain tumor with a poor prognosis. Reactive oxygen species that accumulate during tumorigenesis can cause oxidative stress (OS), which plays a crucial role in cancer cell survival. Clinical and transcriptome data of TCGA-GBM dataset from UCSC Xena database were analyzed. Consensus clustering analysis was conducted to identify OS-related molecular subtypes for GBM. The immune infiltrate level between subtypes were characterized by ESTIMATE algorithm. Differentially expressed genes (DEGs) between subtypes were screened by DESeq2 package. Two OS-related molecular subtypes of GBM were identified, and cluster 2 had poorer overall survival and higher immune infiltration levels than cluster 1. Enrichment analysis showed that 54 DEGs in cluster 2 were significantly enriched in cytokine/chemokine-related functions or pathways. Ten hub genes (CSF2, CSF3, CCL7, LCN2, CXCL6, MMP8, CCR8, TNFSF11, IL22RA2, and ORM1) were identified in GBM subtype 2 through protein-protein interaction network, most of which were positively correlated with immune factors and immune checkpoints. A total of 55 small molecule drugs obtained from drug gene interaction database (DGIdb) may have potential therapeutic effects in GBM subtype 2 patients. Our study identified 10 hub genes as potential therapeutic targets in GBM subtype 2 patients, who have poorer overall survival and higher immune infiltration levels. These findings could pave the way for new treatments for this aggressive form of brain cancer.

Searching for Essential Genes and Targeted Drugs Common to Breast Cancer and Osteoarthritis.

It is documented that osteoarthritis can promote the progression of breast cancer (BC). This study aims to search for the essential genes associated with breast cancer (BC) and osteoarthritis (OA), explore the relationship between epithelial-mesenchymal transition (EMT)- related genes and the two diseases, and identify the candidate drugs. The genes related to both BC and OA were determined by text mining. Protein-protein Interaction (PPI) analysis was carried out, and as a result, the exported genes were found to be related to EMT. PPI and the correlation of mRNA of these genes were also analyzed. Different kinds of enrichment analyses were performed on these genes. A prognostic analysis was performed on these genes for examining their expression levels at different pathological stages, in different tissues, and in different immune cells. Drug-gene interaction database was employed for potential drug discovery. A total number of 1422 genes were identified as common to BC and OA and 58 genes were found to be related to EMT. We found that HDAC2 and TGFBR1 were significantly poor in overall survival. High expression of HDAC2 plays a vital role in the increase of pathological stages. Four immune cells might play a role in this process. Fifty-seven drugs were identified that could potentially have therapeutic effects. EMT may be one of the mechanisms by which OA affects BC. Using the drugs can have potential therapeutic effects, which may benefit patients with both diseases and broaden the indications for drug use.

Abstract A001: Transcriptomic and microRNA profiling revealed novel therapeutic targets for breast cancer subtypes

Abstract Numerous studies suggested breast cancer (BC) from the Middle East and North Africa (MENA) to be presented at younger age with advanced tumor grade, suggesting underlying biological differences. Given the limited transcriptomic profiling data on BC from the MENA region, we performed mRNA and microRNA (miRNA) transcriptomic profiling on a local cohort of BC (n=96) from Qatar. Our data revealed the differentially expressed genes and miRNAs in relation to BC molecular subtypes (HR+, HER2+, HER2+HR+, and TNBC), tumor grade (high vs low), patients’ age (young (≤ 40) vs old (&amp;gt; 40)), and ethnicity (MENA vs non-MENA). Our data revealed close similarity between TNBC and HER2+, while the transcriptome of HER2+HR+ tumor was similar to that from HR+ tumors. Integrated network analysis identified complex miRNA-mRNA regulatory networks in each BC molecular subtype, in high vs low grade tumors, in tumors from young vs old, and in tumors from MENA vs non-MENA, thus implicating miRNA-mediated gene regulation as an key mechanism in shaping the transcriptome of BC. Integration of our transcriptomic data with CRISPR-Cas9 functional screen data and the Drug-Gene Interaction Database identified numerous dependencies and therapeutic vulnerabilities for each BC molecular subtype. Among those, CDC123 was functionally validated as potential therapeutic target for TNBC. Survival analyses identified mRNA and miRNA-based signatures predicative of relapse free survival (RFS), which were validated in larger BC cohorts. Our data provides comprehensive transcriptomic profiling and unraveled the miRNA-mRNA regulatory networks in BC patients from the region and identified numerous actionable gene targets, employing integrated approach. Findings from the current study have potential implications to improve the current standard-of-care for BC from the MENA as well as patients from other ethnicities. Citation Format: Ramesh Elango, Radhakrishnan Vishnubalaji, Sameera Rashid, Nehad Alajez. Transcriptomic and microRNA profiling revealed novel therapeutic targets for breast cancer subtypes [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Breast Cancer Research; 2023 Oct 19-22; San Diego, California. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_1):Abstract nr A001.

Transcriptomic Analysis of Metastatic Colorectal Tumor with Low Mutational Burden

Objective: To identify potential drug targets for metastasis colorectal cancer (CRC) patients with low mutational burden by examining differences in immune-related gene expression. Methods: For this study, 623 samples were collected from The Cancer Genome Atlas (TCGA) database, comprising tumor mutational burden (TMB), RNA sequencing (RNA-Seq), and clinical data. Differential gene expression analysis, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the identified genes were conducted using the R package. Additionally, a comparative analysis of immune infiltrating cell composition in metastatic and non-metastatic groups was performed. Hub genes, exhibiting high levels of interaction, were selected using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database. The Drug Gene Interaction Database (DGIdb) was then utilized to estimate drugs targeting the identified hub genes. Results: The transcriptome data of 326 colorectal cancer patients with low TMB were analyzed, comprising 58 patients with metastasis and 268 patients without metastasis. Among the differential expression in 1,111 genes for patients with metastasis compared to those without metastasis, 733 genes were upregulated, and 378 genes were downregulated. KEGG and GO enrichment analysis indicated significant differences in gene expression in CRC metastatic patients with low TMB compared to non-metastasis patients with low TMB. Enriched pathways included humoral immuneresponse, immunoglobulin production, and regulation of AMPA receptor activity. Two genes related to interleukin-12 were identified through secondary enrichment for immune-related genes. Analysis of tumor-infiltrating immune cell data revealed significant differences in memory-activated T cell CD4 and T cell CD8. Conclusions: This analysis of RNA sequencing data and immune-filtrating cell data revealed significant differences between metastatic colorectal cancer patients with low TMB and their non-metastatic counterparts. These distinctions suggest the possibility of identifying more effective drugs or therapies for metastatic colorectal cancer patients with low TMB.

Identification of CCL20 as a Prognostic Predictor for Severe Fever With Thrombocytopenia Syndrome Based on Plasma Proteomics.

Severe fever with thrombocytopenia syndrome (SFTS), a lethal tick-borne hemorrhagic fever, prompted our investigation into prognostic predictors and potential drug targets using plasma Olink Proteomics. Employing the Olink assay, we analyzed 184 plasma proteins in 30 survivors and 8 nonsurvivors of SFTS. Validation was performed in a cohort of 154 patients with SFTS via enzyme-linked immunosorbent assay. We utilized the Drug-Gene Interaction Database to identify protein-drug interactions. Nonsurvivors exhibited 110 differentially expressed proteins as compared with survivors, with functional enrichment in the cell chemotaxis-related pathway. Thirteen differentially expressed proteins-including C-C motif chemokine 20 (CCL20), calcitonin gene-related peptide alpha, and pleiotrophin-were associated with multiple-organ dysfunction syndrome. CCL20 emerged as the top predictor of death, demonstrating an area under the curve of 1 (P = .0004) and 0.9033 (P < .0001) in the discovery and validation cohorts, respectively. Patients with CCL20 levels exceeding 45.74 pg/mL exhibited a fatality rate of 45.65%, while no deaths occurred in those with lower CCL20 levels. Furthermore, we identified 202 Food and Drug Administration-approved drugs targeting 37 death-related plasma proteins. Distinct plasma proteomic profiles characterize SFTS cases with different outcomes, with CCL20 emerging as a novel, sensitive, accurate, and specific biomarker for predicting SFTS prognosis.

Whole transcriptome sequencing identifies key lncRNAs,circRNAs, and mRNAs for exploring the pathogenesis and therapeutic target of mouse pneumoconiosis

BackgroundPneumoconiosis is a kind of lung dysfunction caused by the inhalation of mineral dust. However, the potential molecular mechanism of pneumoconiosis have not been fully elucidated. MethodsIn this study, the silica-treated pneumoconiosis mice model was constructed and the transcriptome sequencing data including lncRNA, circRNA, and mRNA were obtained. Firstly, differentially expressed lncRNA, circRNA, and mRNA (DElncRNA, DEcircRNA, DEGs) between control and pneumoconiosis/silicosis samples were screened, the target miRNAs (co-pre-miRNAs) were obtained by intersecting the miRNAs predicted by DElncRNA and DEcircRNA, respectively, and the target mRNAs (co-mRNA) were obtained by intersecting the mRNAs predicted by target miRNA and DEGs. Then, the lncRNA/circRNA-miRNA-mRNA networks were constructed by Cytoscape. Next, the key mRNAs were obtained by protein–protein interaction (PPI) analysis, and the key lncRNAs/circRNAs were selected by correlation analysis. Moreover, the expression of the key lncRNAs, circRNAs and mRNAs on chromosome were studied by the “circlize” package. Furthermore, the TFs-miRNA-mRNA network was constructed and the function of DEGs were explored by Ingenuity Pathway Analysis (IPA). To demonstrate the feasibility and value of the constructed ceRNA networks, we validated key genes and mmu-miR-682 pathway. Finally, We used the Drug-Gene Interaction database to predict potential drugs that could interfere with key genes,which may help to find promising treatment. ResultsThere were 427 DElncRNAs, 107 DEcircRNAs and 1,597 DEGs between silicosis and control groups. Totals of 77 co-pre-miRNAs and 96 co-mRNA were screened, and the lncRNA/circRNA-miRNA-mRNA networks were constructed with 27 lncRNA/25 circRNAs, 74 miRNAs and 96 mRNAs. Then, 6 key mRNAs including Igf1, Klf4, Ptgs2, Epas1, Gnao1, and Il1a were obtained by PPI, and all of these key mRNAs and 10 key lncRNAs and 8 circRNAs were significantly different between the pneumoconiosis and normal groups, in which 10 lncRNAs and 9 circRNA that have not been previously studied in pneumoconiosis/silicosis can be used as new potential therapeutic targets. Moreover, the TFs-miRNA-mRNA network were constructed with 11 TFs, 1 key miRNA (mmu-miR-682) and 3 key mRNAs (Igf1, Epas1, Ptgs2). And the validation of key genes revealing by RNA-seq through experimental approaches shows the the predictive power of this study. Finally, IPA results indicated that 41 pathways were activated and 2 pathways were suppressed in pneumoconiosis/silicosis groups, and Pathogen Induced Cytokine Storm Signaling Pathway was the most significant pathway affected by pneumoconiosis/silicosis. In addition, 93 drugs were screened out by Drug-Gene Interaction database. Among them, Hydroxychloroquine was a kind of drug which associated with Il1a and Ptgs2, may be a promising treatment. ConclusionThis study constructed the lncRNA/circRNA-miRNA-mRNA and TFs-miRNA-mRNA networks, which could deepen the potential molecular regulatory mechanism of pneumoconiosis/silicosis.