Integration of the Drug-Gene Interaction Database (DGIdb 4.0) with open crowdsource efforts

mRNA-miRNA bipartite networks reconstruction in different tissues of bladder cancer based on gene co-expression network analysis

Bladder cancer (BC) is one of the most important cancers worldwide, and if it is diagnosed early, its progression in humans can be prevented and long-term survival will be achieved accordingly. This study aimed to identify novel micro-RNA (miRNA) and gene-based biomarkers for diagnosing BC. The microarray dataset of BC tissues (GSE13507) listed in the GEO database was analyzed for this purpose. The gene expression data from three BC tissues including 165 primary bladder cancer (PBC), 58 normal looking-bladder mucosae surrounding cancer (NBMSC), and 23 recurrent non-muscle invasive tumor tissues (RNIT) were used to reconstruct gene co-expression networks. After preprocessing and normalization, deferentially expressed genes (DEGs) were obtained and used to construct the weighted gene co-expression network (WGCNA). Gene co-expression modules and low-preserved modules were extracted among BC tissues using network clustering. Next, the experimentally validated mRNA-miRNA interaction information were used to reconstruct three mRNA-miRNA bipartite networks. Reactome pathway database and Gene ontology (GO) was subsequently performed for the extracted genes of three bipartite networks and miRNAs, respectively. To further analyze the data, ten hub miRNAs (miRNAs with the highest degree) were selected in each bipartite network to reconstruct three bipartite subnetworks. Finally, the obtained biomarkers were comprehensively investigated and discussed in authentic studies. The obtained results from our study indicated a group of genes including PPARD, CST4, CSNK1E, PTPN14, ETV6, and ADRM1 as well as novel miRNAs (e.g., miR-16-5p, miR-335-5p, miR-124-3p, and let-7b-5p) which might be potentially associated with BC and could be a potential biomarker. Afterward, three drug-gene interaction networks were reconstructed to explore candidate drugs for the treatment of BC. The hub miRNAs in the mRNA-miRNA bipartite network played a fundamental role in BC progression; however, these findings need further investigation.

The Diagnostic and Prognostic Values of HOXA Gene Family in Kidney Clear Cell Renal Cell Carcinoma

Kidney renal clear cell carcinoma (KIRC) is one of the most common cancers with high mortality worldwide. As members of the homeobox (HOX) family, homeobox-A (HOXA) genes have been reported to play an increasingly important role in tumorigenesis and the progression of multiple cancers. However, limited studies have investigated the potential diagnostic and prognostic roles of HOXA genes in KIRC. In this research, we explored the expression pattern of the HOXA gene family in KIRC progression by differential analysis of expression profiles from The Cancer Genome Atlas (TCGA). By using univariate Cox analysis and lasso regression analysis, we comprehensively evaluated the prognostic value of HOXA genes and eventually identified a prognostic risk model consisting of five HOXA genes (HOXA2, HOXA3, HOXA7, HOXA11, and HOXA13). The risk model was further validated as a novel independent prognostic factor for KIRC patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses, and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, to explore the potential mechanism of tumorigenesis and clinical application of KIRC, we also developed the HOXA-based competing endogenous RNA (ceRNA) regulatory network and machine learning classification model. Valproic acid and tretinoin were predicted to be the most promising small molecules to adjuvant treatment of KIRC by mining the CMAP and DGIdb drug database. Subsequently, pathway and functional enrichment analyses provided us with new ways to search for a possible mechanism of action of drugs. Taken together, our study demonstrated the nonnegligible role of HOXA genes in KIRC and constructed an effective prognostic and diagnostic model, which offers novel insights into KIRC prognosis.

Cleavage Stimulation Factor Subunit 2: Function Across Cancers and Potential Target for Chemotherapeutic Drugs

The cleavage stimulation factor subunit complex is involved in the cleavage and polyadenylation of 3′-end pre-mRNAs that regulate mRNA formation and processing. However, cleavage stimulation factor subunit 2 (CSTF2) was found to play a more critical regulatory role across cancers. General cancer data sets from The Cancer Genome Atlas and Genotype-Tissue Expression project were thus downloaded for differential analysis, and the possible functions and mechanisms of CSTF2 in general cancer were analyzed using the Compartments database, cBioPortal database, Tumor Immune Single-cell Hub database, and Comparative Toxigenomics database using gene set enrichment analysis and R software. The results showed that CSTF2 could affect DNA repair and methylation in tumor cells. In addition, CSTF2 was associated with multiple tumor immune infiltrates in a wide range of cancers, and its high expression was associated with multiple immune checkpoints; therefore, it could serve as a potential target for many drug molecules. We also proved that CSTF2 promotes oral cell proliferation and migration. The high diagnostic efficacy of CSTF2 suggested that this gene may act as a new biomarker and personalized therapeutic target for a variety of tumors.

The Prognostic Value and Immunological Role of STEAP1 in Pan-Cancer: A Result of Data-Based Analysis

Purpose. This study is aimed at systematically analyzing the expression, function, and prognostic value of six transmembrane epithelial antigen of the prostate 1 (STEAP1) in various cancers. Methods. The expressions of STEAP1 between normal and tumor tissues were analyzed using TCGA and GTEx. Clinicopathologic data was collected from GEPIA and TCGA. Prognostic analysis was conducted by Cox proportional hazard regression and Kaplan-Meier survival. DNA methylation, mutation features, and molecular subtypes of cancers were also investigated. The top-100 coexpressed genes with STEAP1 were involved in functional enrichment analysis. ESTIMATE algorithm was used to analyze the correlation between STEAP1 and immunity value. The relationships of STEAP1 and biomarkers including tumor mutational burden (TMB), microsatellite instability (MSI), and stemness score as well as chemosensitivity were also illustrated. Results. Among 33 cancers, STEAP1 was overexpressed in 19 cancers such as cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), colon adenocarcinoma, and lymphoid neoplasm diffuse large B cell lymphoma while was downregulated in 5 cancers such as adrenocortical carcinoma, breast invasive carcinoma (BRCA), and kidney chromophobe renal cell carcinoma. STEAP1 has significant prognostic relationships in multiple cancers. 15 cancers exhibited differences of DNA methylation including bladder urothelial carcinoma, BRCA, and CESC. STEAP1 expression was positively correlated to immune molecules especially in thyroid carcinoma and negatively especially in uveal melanoma. STEAP1 was associated with TMB and MSI in certain cancers. In addition, STEAP1 was connected with increased chemosensitivity of drugs such as trametinib and pimasertib. Conclusions. STEAP1 was an underlying target for prognostic prediction in different cancer types and a potential biomarker of TMB, MSI, tumor microenvironment, and chemosensitivity.

A Computational Text Mining-Guided Meta-Analysis Approach to Identify Potential Xerostomia Drug Targets

Xerostomia (subjective complaint of dry mouth) is commonly associated with salivary gland hypofunction. Molecular mechanisms associated with xerostomia pathobiology are poorly understood, thus hampering drug development. Our objectives were to (i) use text-mining tools to investigate xerostomia and dry mouth concepts, (ii) identify associated molecular interactions involving genes as candidate drug targets, and (iii) determine how drugs currently used in clinical trials may impact these genes and associated pathways. PubMed and PubMed Central were used to identify search terms associated with xerostomia and/or dry mouth. Search terms were queried in pubmed2ensembl. Protein-protein interaction (PPI) networks were determined using the gene/protein network visualization program search tool for recurring instances of neighboring genes (STRING). A similar program, Cytoscape, was used to determine PPIs of overlapping gene sets. The drug-gene interaction database (DGIdb) and the clinicaltrials.gov database were used to identify potential drug targets from the xerostomia/dry mouth PPI gene set. We identified 64 search terms in common between xerostomia and dry mouth. STRING confirmed PPIs between identified genes (CL = 0.90). Cytoscape analysis determined 58 shared genes, with cytokine-cytokine receptor interaction representing the most significant pathway (p = 1.29 × 10-23) found in the Kyoto encyclopedia of genes and genomes (KEGG). Fifty-four genes in common had drug interactions, per DGIdb analysis. Eighteen drugs, targeting the xerostomia/dry mouth PPI network, have been evaluated for xerostomia, head and neck cancer oral complications, and Sjögren's Syndrome. The PPI network genes IL6R, EGFR, NFKB1, MPO, and TNFSF13B constitute a possible biomarker signature of xerostomia. Validation of the candidate biomarkers is necessary to better stratify patients at the genetic and molecular levels to facilitate drug development or to monitor response to treatment.

A Hybrid Machine Learning and Network Analysis Approach Reveals Two Parkinson's Disease Subtypes from 115 RNA-Seq Post-Mortem Brain Samples

Precision medicine emphasizes fine-grained diagnostics, taking individual variability into account to enhance treatment effectiveness. Parkinson’s disease (PD) heterogeneity among individuals proves the existence of disease subtypes, so subgrouping patients is vital for better understanding disease mechanisms and designing precise treatment. The purpose of this study was to identify PD subtypes using RNA-Seq data in a combined pipeline including unsupervised machine learning, bioinformatics, and network analysis. Two hundred and ten post mortem brain RNA-Seq samples from PD (n = 115) and normal controls (NCs, n = 95) were obtained with systematic data retrieval following PRISMA statements and a fully data-driven clustering pipeline was performed to identify PD subtypes. Bioinformatics and network analyses were performed to characterize the disease mechanisms of the identified PD subtypes and to identify target genes for drug repurposing. Two PD clusters were identified and 42 DEGs were found (p adjusted ≤ 0.01). PD clusters had significantly different gene network structures (p < 0.0001) and phenotype-specific disease mechanisms, highlighting the differential involvement of the Wnt/β-catenin pathway regulating adult neurogenesis. NEUROD1 was identified as a key regulator of gene networks and ISX9 and PD98059 were identified as NEUROD1-interacting compounds with disease-modifying potential, reducing the effects of dopaminergic neurodegeneration. This hybrid data analysis approach could enable precision medicine applications by providing insights for the identification and characterization of pathological subtypes. This workflow has proven useful on PD brain RNA-Seq, but its application to other neurodegenerative diseases is encouraged.

Computational Drug Repurposing Based on a Recommendation System and Drug-Drug Functional Pathway Similarity

Drug repurposing identifies new clinical indications for existing drugs. It can be used to overcome common problems associated with cancers, such as heterogeneity and resistance to established therapies, by rapidly adapting known drugs for new treatment. In this study, we utilized a recommendation system learning model to prioritize candidate cancer drugs. We designed a drug-drug pathway functional similarity by integrating multiple genetic and epigenetic alterations such as gene expression, copy number variation (CNV), and DNA methylation. When compared with other similarities, such as SMILES chemical structures and drug targets based on the protein-protein interaction network, our approach provided better interpretable models capturing drug response mechanisms. Furthermore, our approach can achieve comparable accuracy when evaluated with other learning models based on large public datasets (CCLE and GDSC). A case study about the Erlotinib and OSI-906 (Linsitinib) indicated that they have a synergistic effect to reduce the growth rate of tumors, which is an alternative targeted therapy option for patients. Taken together, our computational method characterized drug response from the viewpoint of a multi-omics pathway and systematically predicted candidate cancer drugs with similar therapeutic effects.

Prognostic tools and candidate drugs based on plasma proteomics of patients with severe COVID-19 complications

COVID-19 complications still present a huge burden on healthcare systems and warrant predictive risk models to triage patients and inform early intervention. Here, we profile 893 plasma proteins from 50 severe and 50 mild-moderate COVID-19 patients, and 50 healthy controls, and show that 375 proteins are differentially expressed in the plasma of severe COVID-19 patients. These differentially expressed plasma proteins are implicated in the pathogenesis of COVID-19 and present targets for candidate drugs to prevent or treat severe complications. Based on the plasma proteomics and clinical lab tests, we also report a 12-plasma protein signature and a model of seven routine clinical tests that validate in an independent cohort as early risk predictors of COVID-19 severity and patient survival. The risk predictors and candidate drugs described in our study can be used and developed for personalized management of SARS-CoV-2 infected patients.

Prognostic markers for patients with COVID-19 are of critical importance in determining the course of SARS-CoV-2 infection and patient handling. Here the authors determine and apply a prognostic proteomic panel for risk and drug prediction in the management of SARS-CoV-2 infected patients.

Deciphering the Interactions of SARS-CoV-2 Proteins with Human Ion Channels Using Machine-Learning-Based Methods

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is accountable for the protracted COVID-19 pandemic. Its high transmission rate and pathogenicity led to health emergencies and economic crisis. Recent studies pertaining to the understanding of the molecular pathogenesis of SARS-CoV-2 infection exhibited the indispensable role of ion channels in viral infection inside the host. Moreover, machine learning (ML)-based algorithms are providing a higher accuracy for host-SARS-CoV-2 protein–protein interactions (PPIs). In this study, PPIs of SARS-CoV-2 proteins with human ion channels (HICs) were trained on the PPI-MetaGO algorithm. PPI networks (PPINs) and a signaling pathway map of HICs with SARS-CoV-2 proteins were generated. Additionally, various U.S. food and drug administration (FDA)-approved drugs interacting with the potential HICs were identified. The PPIs were predicted with 82.71% accuracy, 84.09% precision, 84.09% sensitivity, 0.89 AUC-ROC, 65.17% Matthews correlation coefficient score (MCC) and 84.09% F1 score. Several host pathways were found to be altered, including calcium signaling and taste transduction pathway. Potential HICs could serve as an initial set to the experimentalists for further validation. The study also reinforces the drug repurposing approach for the development of host directed antiviral drugs that may provide a better therapeutic management strategy for infection caused by SARS-CoV-2.

Development of Gene Expression-Based Random Forest Model for Predicting Neoadjuvant Chemotherapy Response in Triple-Negative Breast Cancer

Simple Summary

Only 20–50% of patients with triple negative breast cancer achieve a pathological complete response from neoadjuvant chemotherapy, a strong indicator of patient survival. Therefore, there is an urgent need for a reliable predictive model of the patient’s pathological complete response prior to actual treatment. The purpose of this study was to develop such a model based on random forest recursive feature elimination and to benchmark the performance of the proposed model against existing predictive models. Our study suggests that an 86-gene-based random forest model associated to DNA repair and cell cycle mechanisms can provide reliable predictions of neoadjuvant chemotherapy response in patients with triple negative breast cancer.

Abstract

Neoadjuvant chemotherapy (NAC) response is an important indicator of patient survival in triple negative breast cancer (TNBC), but predicting chemosensitivity remains a challenge in clinical practice. We developed an 86-gene-based random forest (RF) classifier capable of predicting neoadjuvant chemotherapy response (pathological Complete Response (pCR) or Residual Disease (RD)) in TNBC patients. The performance of pCR classification of the proposed model was evaluated by Receiver Operating Characteristic (ROC) curve and Precision Recall (PR) curve. The AUROC and AUPRC of the proposed model on the test set were 0.891 and 0.829, respectively. At a predefined specificity (>90%), the proposed model shows a superior sensitivity compared to the best performing reported NAC response prediction model (69.2% vs. 36.9%). Moreover, the predicted pCR status by the model well explains the distance recurrence free survival (DRFS) of TNBC patients. In addition, the pCR probabilities of the proposed model using the expression profiles of the CCLE TNBC cell lines show a high Spearman rank correlation with cyclophosphamide sensitivity in the TNBC cell lines (SRCC =0.697, p-value =0.031). Associations between the 86 genes and DNA repair/cell cycle mechanisms were provided through function enrichment analysis. Our study suggests that the random forest-based prediction model provides a reliable prediction of the clinical response to neoadjuvant chemotherapy and may explain chemosensitivity in TNBC.

Correlation of age of onset and clinical severity in Niemann-Pick disease type C1 with lysosomal abnormalities and gene expression

Niemann-Pick disease type C1 (NPC1) is a rare, prematurely fatal lysosomal storage disorder which exhibits highly variable severity and disease progression as well as a wide-ranging age of onset, from perinatal stages to adulthood. This heterogeneity has made it difficult to obtain prompt diagnosis and to predict disease course. In addition, small NPC1 patient sample sizes have been a limiting factor in acquiring genome-wide transcriptome data. In this study, primary fibroblasts from an extensive cohort of 41 NPC1 patients were used to validate our previous findings that the lysosomal quantitative probe LysoTracker can be used as a predictor for age of onset and disease severity. We also examined the correlation between these clinical parameters and RNA expression data from primary fibroblasts and identified a set of genes that were significantly associated with lysosomal defects or age of onset, in particular neurological symptom onset. Hierarchical clustering showed that these genes exhibited distinct expression patterns among patient subgroups. This study is the first to collect transcriptomic data on such a large scale in correlation with clinical and cellular phenotypes, providing a rich genomic resource to address NPC1 clinical heterogeneity and discover potential biomarkers, disease modifiers, or therapeutic targets.

A Practical Strategy for Exploring the Pharmacological Mechanism of Luteolin Against COVID-19/Asthma Comorbidity: Findings of System Pharmacology and Bioinformatics Analysis

Asthma patients may increase their susceptibility to SARS-CoV-2 infection and the poor prognosis of coronavirus disease 2019 (COVID-19). However, anti-COVID-19/asthma comorbidity approaches are restricted on condition. Existing evidence indicates that luteolin has antiviral, anti-inflammatory, and immune regulation capabilities. We aimed to evaluate the possibility of luteolin evolving into an ideal drug and explore the underlying molecular mechanisms of luteolin against COVID-19/asthma comorbidity. We used system pharmacology and bioinformatics analysis to assess the physicochemical properties and biological activities of luteolin and further analyze the binding activities, targets, biological functions, and mechanisms of luteolin against COVID-19/asthma comorbidity. We found that luteolin may exert ideal physicochemical properties and bioactivity, and molecular docking analysis confirmed that luteolin performed effective binding activities in COVID-19/asthma comorbidity. Furthermore, a protein–protein interaction network of 538 common targets between drug and disease was constructed and 264 hub targets were obtained. Then, the top 6 hub targets of luteolin against COVID-19/asthma comorbidity were identified, namely, TP53, AKT1, ALB, IL-6, TNF, and VEGFA. Furthermore, the enrichment analysis suggested that luteolin may exert effects on virus defense, regulation of inflammation, cell growth and cell replication, and immune responses, reducing oxidative stress and regulating blood circulation through the Toll-like receptor; MAPK, TNF, AGE/RAGE, EGFR, ErbB, HIF-1, and PI3K–AKT signaling pathways; PD-L1 expression; and PD-1 checkpoint pathway in cancer. The possible “dangerous liaison” between COVID-19 and asthma is still a potential threat to world health. This research is the first to explore whether luteolin could evolve into a drug candidate for COVID-19/asthma comorbidity. This study indicated that luteolin with superior drug likeness and bioactivity has great potential to be used for treating COVID-19/asthma comorbidity, but the predicted results still need to be rigorously verified by experiments.

Regulatory Variant rs2535629 in ITIH3 Intron Confers Schizophrenia Risk By Regulating CTCF Binding and SFMBT1 Expression

Genome-wide association studies have identified 3p21.1 as a robust risk locus for schizophrenia. However, the underlying molecular mechanisms remain elusive. Here a functional regulatory variant (rs2535629) is identified that disrupts CTCF binding at 3p21.1. It is confirmed that rs2535629 is also significantly associated with schizophrenia in Chinese population and the regulatory effect of rs2535629 is validated. Expression quantitative trait loci analysis indicates that rs2535629 is associated with the expression of three distal genes (GLT8D1, SFMBT1, and NEK4) in the human brain, and CRISPR-Cas9-mediated genome editing confirmed the regulatory effect of rs2535629 on GLT8D1, SFMBT1, and NEK4. Interestingly, differential expression analysis of GLT8D1, SFMBT1, and NEK4 suggested that rs2535629 may confer schizophrenia risk by regulating SFMBT1 expression. It is further demonstrated that Sfmbt1 regulates neurodevelopment and dendritic spine density, two key pathological characteristics of schizophrenia. Transcriptome analysis also support the potential role of Sfmbt1 in schizophrenia pathogenesis. The study identifies rs2535629 as a plausibly causal regulatory variant at the 3p21.1 risk locus and demonstrates the regulatory mechanism and biological effect of this functional variant, indicating that this functional variant confers schizophrenia risk by altering CTCF binding and regulating expression of SFMBT1, a distal gene which plays important roles in neurodevelopment and synaptic morphogenesis.

An Integrative Network Approach to Identify Common Genes for the Therapeutics in Tuberculosis and Its Overlapping Non-Communicable Diseases

Tuberculosis (TB) is the leading cause of death from a single infectious agent. The estimated total global TB deaths in 2019 were 1.4 million. The decline in TB incidence rate is very slow, while the burden of noncommunicable diseases (NCDs) is exponentially increasing in low- and middle-income countries, where the prevention and treatment of TB disease remains a great burden, and there is enough empirical evidence (scientific evidence) to justify a greater research emphasis on the syndemic interaction between TB and NCDs. The current study was proposed to build a disease-gene network based on overlapping TB with NCDs (overlapping means genes involved in TB and other/s NCDs), such as Parkinson's disease, cardiovascular disease, diabetes mellitus, rheumatoid arthritis, and lung cancer. We compared the TB-associated genes with genes of its overlapping NCDs to determine the gene-disease relationship. Next, we constructed the gene interaction network of disease-genes by integrating curated and experimentally validated interactions in humans and find the 13 highly clustered modules in the network, which contains a total of 86 hub genes that are commonly associated with TB and its overlapping NCDs, which are largely involved in the Inflammatory response, cellular response to cytokine stimulus, response to cytokine, cytokine-mediated signaling pathway, defense response, response to stress and immune system process. Moreover, the identified hub genes and their respective drugs were exploited to build a bipartite network that assists in deciphering the drug-target interaction, highlighting the influential roles of these drugs on apparently unrelated targets and pathways. Targeting these hub proteins by using drugs combination or drug repurposing approaches will improve the clinical conditions in comorbidity, enhance the potency of a few drugs, and give a synergistic effect with better outcomes. Thus, understanding the Mycobacterium tuberculosis (Mtb) infection and associated NCDs is a high priority to contain its short and long-term effects on human health. Our network-based analysis opens a new horizon for more personalized treatment, drug-repurposing opportunities, investigates new targets, multidrug treatment, and can uncover several side effects of unrelated drugs for TB and its overlapping NCDs.

Identification of Molecular Subtypes and Potential Small-Molecule Drugs for Esophagus Cancer Treatment Based on m6A Regulators

BACKGROUND

Esophagus cancer (ESCA) is the sixth most frequent cancer in males, with 5-year overall survival of 15%-25%. RNA modifications function critically in cancer progression, and m⁶A regulators are associated with ESCA prognosis. This study further revealed correlations between m⁶A and ESCA development.

METHODS

Univariate Cox regression analysis and consensus clustering were applied to determine molecular subtypes. Functional pathways and gene ontology terms were enriched by gene set enrichment analysis. Protein-protein interaction (PPI) analysis on differentially expressed genes (DEGs) was conducted for hub gene screening. Public drug databases were employed to study the interactions between hub genes and small molecules.

RESULTS

Three molecular subtypes related to ESCA prognosis were determined. Based on multiple analyses among molecular subtypes, 146 DEGs were screened, and a PPT network of 15 hub genes was visualized. Finally, 8 potential small-molecule drugs (BMS-754807, gefitinib, neratinib, zuclopenthixol, puromycin, sulfasalazine, and imatinib) were identified for treating ESCA.

CONCLUSIONS

This study applied a new approach to analyzing the relation between m⁶A and ESCA prognosis, providing a reference for exploring potential targets and drugs for ESCA treatment.

Integrating gene expression and clinical data to identify drug repurposing candidates for hyperlipidemia and hypertension

Discovering novel uses for existing drugs, through drug repurposing, can reduce the time, costs, and risk of failure associated with new drug development. However, prioritizing drug repurposing candidates for downstream studies remains challenging. Here, we present a high-throughput approach to identify and validate drug repurposing candidates. This approach integrates human gene expression, drug perturbation, and clinical data from publicly available resources. We apply this approach to find drug repurposing candidates for two diseases, hyperlipidemia and hypertension. We screen >21,000 compounds and replicate ten approved drugs. We also identify 25 (seven for hyperlipidemia, eighteen for hypertension) drugs approved for other indications with therapeutic effects on clinically relevant biomarkers. For five of these drugs, the therapeutic effects are replicated in the All of Us Research Program database. We anticipate our approach will enable researchers to integrate multiple publicly available datasets to identify high priority drug repurposing opportunities for human diseases.

Priority index: database of genetic targets in immune-mediated disease

We describe a comprehensive and unique database 'Priority index' (Pi; http://pi.well.ox.ac.uk) of prioritized genes encoding potential therapeutic targets that encompasses all major immune-mediated diseases. We provide targets at the gene level, each receiving a 5-star rating supported by: genomic evidence arising from disease genome-wide associations and functional immunogenomics, annotation evidence using ontologies restricted to genes with genomic evidence, and network evidence from protein interactions. Target genes often act together in related molecular pathways. The underlying Pi approach is unique in identifying a network of highly rated genes that mediate pathway crosstalk. In the Pi website, disease-centric pages are specially designed to enable the users to browse a complete list of prioritized genes and also a manageable list of nodal genes at the pathway crosstalk level; both switchable by clicks. Moreover, target genes are cross-referenced and supported using additional information, particularly regarding tractability, including druggable pockets viewed in 3D within protein structures. Target genes highly rated across diseases suggest drug repurposing opportunity, while genes in a particular disease reveal disease-specific targeting potential. To facilitate the ease of such utility, cross-disease comparisons involving multiple diseases are also supported. This facility, together with the faceted search, enhances integrative mining of the Pi resource to accelerate early-stage therapeutic target identification and validation leveraging human genetics.

Inflammation Subtypes and Translating Inflammation-Related Genetic Findings in Schizophrenia and Related Psychoses: A Perspective on Pathways for Treatment Stratification and Novel Therapies

Dysregulation of immunological and inflammatory processes is frequently observed in psychotic disorders. Numerous studies have examined the complex components of innate and adaptive immune processes in schizophrenia and related psychoses. Elevated inflammation in these conditions is related to neurobiological phenotypes and associated with both genetics and environmental exposures. Recent studies have utilized multivariate cytokine approaches to identify what appears to be a subset of individuals with elevated inflammation. The degree to which these findings represent a general process of dysregulated inflammation or whether there are more refined subtypes remains unclear. Brain-imaging studies have attempted to establish the link between peripheral inflammation and gray matter disruption, white matter abnormalities, and neuropsychological phenotypes. However, the interplay between peripheral inflammation and neuroinflammation, as well as the consequences of this interplay, in the context of psychosis remains unclear and requires further investigation. This Perspectives article reviews the following elements of immune dysregulation and its clinical and therapeutic implications: (1) evidence supporting inflammation and immune dysregulation in schizophrenia and related psychoses; (2) recent advances in approaches to characterizing subgroups of patients with elevated inflammation; (3) relationships between peripheral inflammation and brain-imaging indicators of neuroinflammation; (4) convergence of large-scale genetic findings and peripheral inflammation findings; and (5) therapeutic implications: anti-inflammation interventions leveraging genetic findings for drug discovery and repurposing. We offer perspectives and examples of how multiomics technologies may be useful for constructing and studying immunogenetic signatures. Advancing research in this area will facilitate biomarker discovery, disease subtyping, and the development of etiological treatments for immune dysregulation in psychosis.

Application of deep learning algorithm on whole genome sequencing data uncovers structural variants associated with multiple mental disorders in African American patients

Mental disorders present a global health concern, while the diagnosis of mental disorders can be challenging. The diagnosis is even harder for patients who have more than one type of mental disorder, especially for young toddlers who are not able to complete questionnaires or standardized rating scales for diagnosis. In the past decade, multiple genomic association signals have been reported for mental disorders, some of which present attractive drug targets. Concurrently, machine learning algorithms, especially deep learning algorithms, have been successful in the diagnosis and/or labeling of complex diseases, such as attention deficit hyperactivity disorder (ADHD) or cancer. In this study, we focused on eight common mental disorders, including ADHD, depression, anxiety, autism, intellectual disabilities, speech/language disorder, delays in developments, and oppositional defiant disorder in the ethnic minority of African Americans. Blood-derived whole genome sequencing data from 4179 individuals were generated, including 1384 patients with the diagnosis of at least one mental disorder. The burden of genomic variants in coding/non-coding regions was applied as feature vectors in the deep learning algorithm. Our model showed ~65% accuracy in differentiating patients from controls. Ability to label patients with multiple disorders was similarly successful, with a hamming loss score less than 0.3, while exact diagnostic matches are around 10%. Genes in genomic regions with the highest weights showed enrichment of biological pathways involved in immune responses, antigen/nucleic acid binding, chemokine signaling pathway, and G-protein receptor activities. A noticeable fact is that variants in non-coding regions (e.g., ncRNA, intronic, and intergenic) performed equally well as variants in coding regions; however, unlike coding region variants, variants in non-coding regions do not express genomic hotspots whereas they carry much more narrow standard deviations, indicating they probably serve as alternative markers.

Bioinformatics and Network Pharmacology Identify the Therapeutic Role and Potential Mechanism of Melatonin in AD and Rosacea

Rosacea is significantly associated with dementia, particularly Alzheimer’s disease (AD). However, the common underlying molecular mechanism connecting these two diseases remains limited. This study aimed to reveal the common molecular regulatory networks and identify the potential therapeutic drugs for rosacea and AD. There were 747 overlapped DEGs (ol-DEGs) that were detected in AD and rosacea, enriched in inflammation-, metabolism-, and apoptosis-related pathways. Using the TF regulatory network analysis, 37 common TFs and target genes were identified as hub genes. They were used to predict the therapeutic drugs for rosacea and AD using the DGIdb/CMap database. Among the 113 predicted drugs, melatonin (MLT) was co-associated with both RORA and IFN-γ in AD and rosacea. Subsequently, network pharmacology analysis identified 19 pharmacological targets of MLT and demonstrated that MLT could help in treating AD/rosacea partly by modulating inflammatory and vascular signaling pathways. Finally, we verified the therapeutic role and mechanism of MLT on rosacea in vivo and in vitro. We found that MLT treatment significantly improved rosacea-like skin lesion by reducing keratinocyte-mediated inflammatory cytokine secretion and repressing the migration of HUVEC cells. In conclusion, this study contributes to common pathologies shared by rosacea and AD and identified MLT as an effective treatment strategy for rosacea and AD via regulating inflammation and angiogenesis.

A comprehensive gene-centric pleiotropic association analysis for 14 psychiatric disorders with GWAS summary statistics

BACKGROUND

Recent genome-wide association studies (GWASs) have revealed the polygenic nature of psychiatric disorders and discovered a few of single-nucleotide polymorphisms (SNPs) associated with multiple psychiatric disorders. However, the extent and pattern of pleiotropy among distinct psychiatric disorders remain not completely clear.

METHODS

We analyzed 14 psychiatric disorders using summary statistics available from the largest GWASs by far. We first applied the cross-trait linkage disequilibrium score regression (LDSC) to estimate genetic correlation between disorders. Then, we performed a gene-based pleiotropy analysis by first aggregating a set of SNP-level associations into a single gene-level association signal using MAGMA. From a methodological perspective, we viewed the identification of pleiotropic associations across the entire genome as a high-dimensional problem of composite null hypothesis testing and utilized a novel method called PLACO for pleiotropy mapping. We ultimately implemented functional analysis for identified pleiotropic genes and used Mendelian randomization for detecting causal association between these disorders.

RESULTS

We confirmed extensive genetic correlation among psychiatric disorders, based on which these disorders can be grouped into three diverse categories. We detected a large number of pleiotropic genes including 5884 associations and 2424 unique genes and found that differentially expressed pleiotropic genes were significantly enriched in pancreas, liver, heart, and brain, and that the biological process of these genes was remarkably enriched in regulating neurodevelopment, neurogenesis, and neuron differentiation, offering substantial evidence supporting the validity of identified pleiotropic loci. We further demonstrated that among all the identified pleiotropic genes there were 342 unique ones linked with 6353 drugs with drug-gene interaction which can be classified into distinct types including inhibitor, agonist, blocker, antagonist, and modulator. We also revealed causal associations among psychiatric disorders, indicating that genetic overlap and causality commonly drove the observed co-existence of these disorders.

CONCLUSIONS

Our study is among the first large-scale effort to characterize gene-level pleiotropy among a greatly expanded set of psychiatric disorders and provides important insight into shared genetic etiology underlying these disorders. The findings would inform psychiatric nosology, identify potential neurobiological mechanisms predisposing to specific clinical presentations, and pave the way to effective drug targets for clinical treatment.

Combining Human Genetics of Multiple Sclerosis with Oxidative Stress Phenotype for Drug Repositioning

In multiple sclerosis (MS), oxidative stress (OS) is implicated in the neurodegenerative processes that occur from the beginning of the disease. Unchecked OS initiates a vicious circle caused by its crosstalk with inflammation, leading to demyelination, axonal damage and neuronal loss. The failure of MS antioxidant therapies relying on the use of endogenous and natural compounds drives the application of novel approaches to assess target relevance to the disease prior to preclinical testing of new drug candidates. To identify drugs that can act as regulators of intracellular oxidative homeostasis, we applied an in silico approach that links genome-wide MS associations and molecular quantitative trait loci (QTLs) to proteins of the OS pathway. We found 10 drugs with both central nervous system and oral bioavailability, targeting five out of the 21 top-scoring hits, including arginine methyltransferase (CARM1), which was first linked to MS. In particular, the direction of brain expression QTLs for CARM1 and protein kinase MAPK1 enabled us to select BIIB021 and PEITC drugs with the required target modulation. Our study highlights OS-related molecules regulated by functional MS variants that could be targeted by existing drugs as a supplement to the approved disease-modifying treatments.

Exploring metformin as a candidate drug for rosacea through network pharmacology and experimental validation

Rosacea is a common chronic inflammatory disease that affects the middle of the face. Due to the unclear pathogenesis, the effective treatment options for rosacea remain limited. In this study, weighted gene co-expression network analyses (WGCNA) identified three rosacea-related hub modules, which were involved in immune-, metabolic- and development- related signaling pathways. Next, the key genes from green and brown modules were submitted to CMap database for drug prediction and metformin was identified as a candidate drug for rosacea. Moreover, network pharmacology analysis identified pharmacological targets of metformin and demonstrated that metformin could help in treating rosacea partly by modulating inflammatory and angiogenesis signaling pathways. Finally, we verified the therapeutic role and mechanism of metformin on rosacea in vivo and vitro. We found that metformin treatment significantly improved rosacea-like skin lesions including immune cells infiltration, cytokines/chemokines expression and angiogenesis. Moreover, metformin suppressed LL37- and TNF-α-induced the ROS production and MAPK-NF-κB signal activation in keratinocytes cells. In conclusion, our findings identified and verified metformin as a novel therapeutic candidate for rosacea, and it alleviates the pathological symptoms, possibly by suppressing inflammatory responses, angiogenesis in rosacea.

Histopathological and Molecular Profiling of Clear Cell Sarcoma and Correlation with Response to Crizotinib: An Exploratory Study Related to EORTC 90101 "CREATE" Trial

Simple Summary

Clear cell sarcoma (CCSA) is a rare subtype of soft tissue sarcoma characterized by EWSR1 rearrangement and subsequent MET upregulation. The European Organisation for Research and Treatment of Cancer 90101 phase II trial evaluated the MET inhibitor crizotinib in CCSA but resulted in only sporadic responses. The aim of this exploratory study was to identify the molecular alterations potentially relevant for the treatment outcome by using archival CCSA samples and trial-related clinical data. We characterized MET signaling and revealed an infrequent activation of MET, which may explain the lack of response to crizotinib in the disease cohort. Based on sequencing analyses, we discovered copy number alterations, mutations and dysregulated pathways with potentially predictive or prognostic values for patients’ outcomes. This work describes the molecular heterogeneity in CCSA and provides deep insight into the biology of this ultra-rare malignancy, which may potentially lead to better therapeutic approaches.

Abstract

Clear cell sarcoma (CCSA) is characterized by a chromosomal translocation leading to EWSR1 rearrangement, resulting in aberrant transcription of multiple genes, including MET. The EORTC 90101 phase II trial evaluated the MET inhibitor crizotinib in CCSA but resulted in only sporadic responses. We performed an in-depth histopathological and molecular analysis of archival CCSA samples to identify alterations potentially relevant for the treatment outcome. Immunohistochemical characterization of MET signaling was performed using a tissue microarray constructed from 32 CCSA cases. The DNA from 24 available tumor specimens was analyzed by low-coverage whole-genome sequencing and whole-exome sequencing for the detection of recurrent copy number alterations (CNAs) and mutations. A pathway enrichment analysis was performed to identify the pathways relevant for CCSA tumorigenesis. Kaplan–Meier estimates and Fisher’s exact test were used to correlate the molecular findings with the clinical features related to crizotinib treatment, aiming to assess a potential association with the outcomes. The histopathological analysis showed the absence of a MET ligand and MET activation, with the presence of MET itself in most of cases. However, the expression/activation of MET downstream molecules was frequently observed, suggesting the role of other receptors in CCSA signal transduction. Using sequencing, we detected a number of CNAs at the chromosomal arm and region levels. The most common alteration was a gain of 8q24.21, observed in 83% of the cases. The loss of chromosomes 9q and 12q24 was associated with shorter survival. Based on exome sequencing, 40 cancer-associated genes were found to be mutated in more than one sample, with SRGAP3 and KMT2D as the most common alterations (each in four cases). The mutated genes encoded proteins were mainly involved in receptor tyrosine kinase signaling, polymerase-II transcription, DNA damage repair, SUMOylation and chromatin organization. Disruption in chromatin organization was correlated with longer progression-free survival in patients receiving crizotinib. Conclusions: The infrequent activation of MET may explain the lack of response to crizotinib observed in the majority of cases in the clinical trial. Our work describes the molecular heterogeneity in CCSA and provides further insight into the biology of this ultra-rare malignancy, which may potentially lead to better therapeutic approaches for CCSA.

Proteome-wide Association Study Provides Insights Into the Genetic Component of Protein Abundance in Psychiatric Disorders

BACKGROUND

Genome-wide association studies have identified multiple risk variants for psychiatric disorders. Nevertheless, how the risk variants confer risk of psychiatric disorders remains largely unknown.

METHODS

We performed proteome-wide association studies to identify genes whose cis-regulated protein abundance change in the human brain were associated with psychiatric disorders.

RESULTS

By integrating genome-wide associations of four common psychiatric disorders and two independent brain proteomes (n = 376 and n = 152, respectively) from the dorsolateral prefrontal cortex, we identified 61 genes (including 48 genes for schizophrenia, 12 genes for bipolar disorder, 5 genes for depression, and 2 genes for attention-deficit/hyperactivity disorder) whose genetically regulated protein abundance levels were associated with risk of psychiatric disorders. Comparison with transcriptome-wide association studies identified 18 overlapping genes that showed significant associations with psychiatric disorders at both proteome-wide and transcriptome-wide levels, suggesting that genetic risk variants likely confer risk of psychiatric disorders by regulating messenger RNA expression and protein abundance of these genes.

CONCLUSIONS

Our study not only provides new insights into the genetic component of protein abundance in psychiatric disorders but also highlights several high-confidence risk proteins (including CNNM2 and CTNND1) for schizophrenia and depression. These high-confidence risk proteins represent promising therapeutic targets for future drug development.

MIKB: A manually curated and comprehensive knowledge base for myocardial infarction

Myocardial infarction knowledge base (MIKB; http://www.sysbio.org.cn/mikb/; latest update: December 31, 2020) is an open-access and manually curated database dedicated to integrating knowledge about MI to improve the efficiency of translational MI research. MIKB is an updated and expanded version of our previous MI Risk Knowledge Base (MIRKB), which integrated MI-related risk factors and risk models for providing help in risk assessment or diagnostic prediction of MI. The updated MIRKB includes 9701 records with 2054 single factors, 209 combined factors, 243 risk models, 37 MI subtypes and 3406 interactions between single factors and MIs collected from 4817 research articles. The expanded functional module, i.e. MIGD, is a database including not only MI associated genetic variants, but also the other multi-omics factors and the annotations for their functional alterations. The goal of MIGD is to provide a multi-omics level understanding of the molecular pathogenesis of MI. MIGD includes 1782 omics factors, 28 MI subtypes and 2347 omics factor-MI interactions as well as 1253 genes and 6 chromosomal alterations collected from 2647 research articles. The functions of MI associated genes and their interaction with drugs were analyzed. MIKB will be continuously updated and optimized to provide precision and comprehensive knowledge for the study of heterogeneous and personalized MI.

Drug repurposing for coronavirus (SARS-CoV-2) based on gene co-expression network analysis

Severe acute respiratory syndrome (SARS) is a highly contagious viral respiratory illness. This illness is spurred on by a coronavirus known as SARS-associated coronavirus (SARS-CoV). SARS was first detected in Asia in late February 2003. The genome of this virus is very similar to the SARS-CoV-2. Therefore, the study of SARS-CoV disease and the identification of effective drugs to treat this disease can be new clues for the treatment of SARS-Cov-2. This study aimed to discover novel potential drugs for SARS-CoV disease in order to treating SARS-Cov-2 disease based on a novel systems biology approach. To this end, gene co-expression network analysis was applied. First, the gene co-expression network was reconstructed for 1441 genes, and then two gene modules were discovered as significant modules. Next, a list of miRNAs and transcription factors that target gene co-expression modules' genes were gathered from the valid databases, and two sub-networks formed of transcription factors and miRNAs were established. Afterward, the list of the drugs targeting obtained sub-networks' genes was retrieved from the DGIDb database, and two drug-gene and drug-TF interaction networks were reconstructed. Finally, after conducting different network analyses, we proposed five drugs, including FLUOROURACIL, CISPLATIN, SIROLIMUS, CYCLOPHOSPHAMIDE, and METHYLDOPA, as candidate drugs for SARS-CoV-2 coronavirus treatment. Moreover, ten miRNAs including miR-193b, miR-192, miR-215, miR-34a, miR-16, miR-16, miR-92a, miR-30a, miR-7, and miR-26b were found to be significant miRNAs in treating SARS-CoV-2 coronavirus.

Network-guided identification of cancer-selective combinatorial therapies in ovarian cancer

Each patient's cancer consists of multiple cell subpopulations that are inherently heterogeneous and may develop differing phenotypes such as drug sensitivity or resistance. A personalized treatment regimen should therefore target multiple oncoproteins in the cancer cell populations that are driving the treatment resistance or disease progression in a given patient to provide maximal therapeutic effect, while avoiding severe co-inhibition of non-malignant cells that would lead to toxic side effects. To address the intra- and inter-tumoral heterogeneity when designing combinatorial treatment regimens for cancer patients, we have implemented a machine learning-based platform to guide identification of safe and effective combinatorial treatments that selectively inhibit cancer-related dysfunctions or resistance mechanisms in individual patients. In this case study, we show how the platform enables prediction of cancer-selective drug combinations for patients with high-grade serous ovarian cancer using single-cell imaging cytometry drug response assay, combined with genome-wide transcriptomic and genetic profiles. The platform makes use of drug-target interaction networks to prioritize those combinations that warrant further preclinical testing in scarce patient-derived primary cells. During the case study in ovarian cancer patients, we investigated (i) the relative performance of various ensemble learning algorithms for drug response prediction, (ii) the use of matched single-cell RNA-sequencing data to deconvolute cell population-specific transcriptome profiles from bulk RNA-seq data, (iii) and whether multi-patient or patient-specific predictive models lead to better predictive accuracy. The general platform and the comparison results are expected to become useful for future studies that use similar predictive approaches also in other cancer types.

Target Discovery for Host-Directed Antiviral Therapies: Application of Proteomics Approaches

Current epidemics, such as AIDS or flu, and the emergence of new threatening pathogens, such as the one causing the current coronavirus disease 2019 (COVID-19) pandemic, represent major global health challenges. While vaccination is an important part of the arsenal to counter the spread of viral diseases, it presents limitations and needs to be complemented by efficient therapeutic solutions. Intricate knowledge of host-pathogen interactions is a powerful tool to identify host-dependent vulnerabilities that can be exploited to dampen viral replication. Such host-directed antiviral therapies are promising and are less prone to the development of drug-resistant viral strains. Here, we first describe proteomics-based strategies that allow the rapid characterization of host-pathogen interactions. We then discuss how such data can be exploited to help prioritize compounds with potential host-directed antiviral activity that can be tested in preclinical models.

Genomic Sub-Classification of Ovarian Clear Cell Carcinoma Revealed by Distinct Mutational Signatures

Ovarian clear cell carcinoma (OCCC) is characterized by dismal prognosis, partially due to its low sensitivity to standard chemotherapy regimen. It is also well-known for presenting unique molecular features in comparison to other epithelial ovarian cancer subtypes. Here, we aim to identify potential subgroups of patients in order to (1) determine their molecular features and (2) characterize their mutational signature. Furthermore, we sought to perform the investigation based on a potentially clinically relevant setting. To that end, we assessed the mutational profile and genomic instability of 55 patients extracted from the Gynecologic Cancer Database (DGCD) by using a panel comprised of 409 cancer-associated genes and a microsatellite assay, respectively; both are currently used in our routine environment. In accordance with previous findings, ARID1A and PIK3CA were the most prevalent mutations, present in 49.1% and 41.8%, respectively. From those, the co-occurrence of ARID1A and PIK3CA mutations was observed in 36.1% of subjects, indicating that this association might be a common feature of OCCC. The microsatellite instability frequency was low across samples. An unbiased assessment of signatures identified the presence of three subgroups, where "PIK3CA" and "Double hit" (with ARID1A and PIK3CA double mutation) subgroups exhibited unique signatures, whilst "ARID1A" and "Undetermined" (no mutations on ARID1A nor PIK3CA) subgroups showed similar profiles. Those differences were further indicated by COSMIC signatures. Taken together, the current findings suggest that OCCC presents distinct mutational landscapes within its group, which may indicate different therapeutic approaches according to its subgroup. Although encouraging, it is noteworthy that the current results are limited by sample size, and further investigation on a larger group would be crucial to better elucidate them.

Integrative Network Pharmacology of Moringa oleifera Combined with Gemcitabine against Pancreatic Cancer

Gemcitabine (GEM) is the first-line chemotherapy drug for patients with advanced pancreatic cancer. Moringa oleifera (MO) exhibited various biological activities, including anticancer effects. Nevertheless, the effectiveness of their combination against pancreatic cancer has not yet been explored. This study evaluates the effect of MO and GEM against pancreatic cancer through network pharmacology. TCMSP, TCMID, and PubMed were used to identify and screen MO bioactive compounds. MO and GEM genes were predicted through DGIdb, CTD, and DrugBank. Pancreatic cancer genes were retrieved from OMIM and MalaCards. Protein–protein interaction (PPI) and compound-target-pathway network were established via STRING and Cytoscape. Gene ontology (GO) and pathway enrichment analysis were conducted using DAVID Bioinformatic Tools. Catechin, kaempferol, quercetin, and epicatechin that met the drug screening requirements, and three additional compounds, glucomoringin, glucoraphanin, and moringinine, were identified as bioactive compounds in MO. Catechin was found to be the main hub compound in MO. TP53, AKT1, VEGFA, and CCND1 from PPI network were discovered as hub genes to have biological importance in pancreatic cancer. GO and pathway analysis revealed that MO and GEM combination was mainly associated with cancer, including pancreatic cancer, through regulation of apoptosis. Combination therapy between MO and GEM might provide insight in pancreatic cancer treatment.

CTR-DB, an omnibus for patient-derived gene expression signatures correlated with cancer drug response

To date, only some cancer patients can benefit from chemotherapy and targeted therapy. Drug resistance continues to be a major and challenging problem facing current cancer research. Rapidly accumulated patient-derived clinical transcriptomic data with cancer drug response bring opportunities for exploring molecular determinants of drug response, but meanwhile pose challenges for data management, integration, and reuse. Here we present the Cancer Treatment Response gene signature DataBase (CTR-DB, http://ctrdb.ncpsb.org.cn/), a unique database for basic and clinical researchers to access, integrate, and reuse clinical transcriptomes with cancer drug response. CTR-DB has collected and uniformly reprocessed 83 patient-derived pre-treatment transcriptomic source datasets with manually curated cancer drug response information, involving 28 histological cancer types, 123 drugs, and 5139 patient samples. These data are browsable, searchable, and downloadable. Moreover, CTR-DB supports single-dataset exploration (including differential gene expression, receiver operating characteristic curve, functional enrichment, sensitizing drug search, and tumor microenvironment analyses), and multiple-dataset combination and comparison, as well as biomarker validation function, which provide insights into the drug resistance mechanism, predictive biomarker discovery and validation, drug combination, and resistance mechanism heterogeneity.

GWAS Meta-Analysis Reveals Shared Genes and Biological Pathways between Major Depressive Disorder and Insomnia

Major depressive disorder (MDD) is one of the most prevalent and disabling mental disorders worldwide. Among the symptoms of MDD, sleep disturbance such as insomnia is prominent, and the first reason patients may seek professional help. However, the underlying pathophysiology of this comorbidity is still elusive. Recently, genome-wide association studies (GWAS) have begun to unveil the genetic background of several psychiatric disorders, including MDD and insomnia. Identifying the shared genomic risk loci between comorbid psychiatric disorders could be a valuable strategy to understanding their comorbidity. This study seeks to identify the shared genes and biological pathways between MDD and insomnia based on their shared genetic variants. First, we performed a meta-analysis based on the GWAS summary statistics of MDD and insomnia obtained from Psychiatric Genomics Consortium and UK Biobank, respectively. Next, we associated shared genetic variants to genes using two gene mapping strategies: (a) positional mapping based on genomic proximity and (b) expression quantitative trait loci (eQTL) mapping based on gene expression linkage across multiple tissues. As a result, a total of 719 shared genes were identified. Over half (51%) of them are protein-coding genes. Functional enrichment analysis shows that the most enriched biological pathways are related to epigenetic modification, sensory perception, and immunologic signatures. We also identified druggable targets using a network approach. Together, these results may provide insights into understanding the genetic predisposition and underlying biological pathways of comorbid MDD and insomnia symptoms.

Genome-Wide Meta-Analysis Identifies Two Novel Risk Loci for Epilepsy

Epilepsy (affects about 70 million people worldwide) is one of the most prevalent brain disorders and imposes a huge economic burden on society. Epilepsy has a strong genetic component. In this study, we perform the largest genome-wide meta-analysis of epilepsy (N = 8,00,869 subjects) by integrating four large-scale genome-wide association studies (GWASs) of epilepsy. We identified three genome-wide significant (GWS) (p < 5 × 10^–8) risk loci for epilepsy. The risk loci on 7q21.11 [lead single nucleotide polymorphism (SNP) rs11978015, p = 9.26 × 10^–9] and 8p23.1 (lead SNP rs28634186, p = 4.39 × 10^–8) are newly identified in the present study. Of note, rs11978015 resides in upstream of GRM3, which encodes glutamate metabotropic receptor 3. GRM3 has pivotal roles in neurotransmission and is involved in most aspects of normal brain function. In addition, we also identified three genes (TTC21B, RP11-375N15.2, and TNKS) whose cis-regulated expression level are associated with epilepsy, indicating that risk variants may confer epilepsy risk through regulating the expression of these genes. Our study not only provides new insights into genetic architecture of epilepsy but also prioritizes potential molecular targets (including GRM3 and TTC21B) for development of new drugs and therapeutics for epilepsy.

A high-resolution temporal atlas of the SARS-CoV-2 translatome and transcriptome

COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which infected >200 million people resulting in >4 million deaths. However, temporal landscape of the SARS-CoV-2 translatome and its impact on the human genome remain unexplored. Here, we report a high-resolution atlas of the translatome and transcriptome of SARS-CoV-2 for various time points after infecting human cells. Intriguingly, substantial amount of SARS-CoV-2 translation initiates at a novel translation initiation site (TIS) located in the leader sequence, termed TIS-L. Since TIS-L is included in all the genomic and subgenomic RNAs, the SARS-CoV-2 translatome may be regulated by a sophisticated interplay between TIS-L and downstream TISs. TIS-L functions as a strong translation enhancer for ORF S, and as translation suppressors for most of the other ORFs. Our global temporal atlas provides compelling insight into unique regulation of the SARS-CoV-2 translatome and helps comprehensively evaluate its impact on the human genome.

Here, Kim et al. apply various sequencing techniques (RPF-seq, QTI-seq, mRNA-seq, sRNA-seq) to unravel the high-resolution, longitudinal translatome and transcriptome of SARS-CoV-2. They identify a translation initiation site in the leader sequence of all genomic and subgenomic RNAs and show its relevance for the SARS-CoV-2 translatome.

Exploring Chemical Information in PubChem

PubChem (https://pubchem.ncbi.nlm.nih.gov) is a public chemical database that serves scientific communities as well as the general public. This database collects chemical information from hundreds of data sources and organizes them into multiple data collections, including Substance, Compound, BioAssay, Protein, Gene, Pathway, and Patent. These collections are interlinked with each other, allowing users to discover related records in the various collections (e.g., drugs targeting a protein or genes modulated by a chemical). PubChem can be searched by keyword (e.g., a chemical, protein, or gene name) as well as by chemical structure. The input structure can be provided using popular line notations or drawn with the PubChem Sketcher. PubChem supports various types of structure searches, including identity search, 2‐D and 3‐D similarity searches, and substructure and superstructure searches. Results from multiple searches can be combined using Boolean operators (i.e., AND, OR, and NOT) to formulate complex queries. PubChem allows the user to quickly retrieve a list of records annotated with a particular classification or ontological term. This paper provides step‐by‐step instructions on how to explore PubChem data with examples of commonly requested tasks. © 2021. This article is a U.S. Government work and is in the public domain in the USA. Current Protocols published by Wiley Periodicals LLC.

Basic Protocol 1: Finding genes and proteins that interact with a given compound

Basic Protocol 2: Finding drug‐like compounds similar to a query compound through a two‐dimensional (2‐D) similarity search

Basic Protocol 3: Finding compounds similar to a query compound through a three‐dimensional (3‐D) similarity search

Support Protocol: Computing similarity scores between compounds

Basic Protocol 4: Getting the bioactivity data for the hit compounds from substructure search

Basic Protocol 5: Finding drugs that target a particular gene

Basic Protocol 6: Getting bioactivity data of all chemicals tested against a protein.

Basic Protocol 7: Finding compounds annotated with classifications or ontological terms

Basic Protocol 8: Finding stereoisomers and isotopomers of a compound through identity search

RNA-Protein Interaction Analysis of SARS-CoV-2 5' and 3' Untranslated Regions Reveals a Role of Lysosome-Associated Membrane Protein-2a during Viral Infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive-strand RNA virus. The viral genome is capped at the 5′ end, followed by an untranslated region (UTR). There is a poly(A) tail at the 3′ end, preceded by a UTR. The self-interaction between the RNA regulatory elements present within the 5′ and 3′ UTRs and their interaction with host/virus-encoded proteins mediate the function of the 5′ and 3′ UTRs. Using an RNA-protein interaction detection (RaPID) assay coupled to liquid chromatography with tandem mass spectrometry, we identified host interaction partners of SARS-CoV-2 5′ and 3′ UTRs and generated an RNA-protein interaction network. By combining these data with the previously known protein-protein interaction data proposed to be involved in virus replication, we generated the RNA-protein-protein interaction (RPPI) network, likely to be essential for controlling SARS-CoV-2 replication. Notably, bioinformatics analysis of the RPPI network revealed the enrichment of factors involved in translation initiation and RNA metabolism. Lysosome-associated membrane protein-2a (Lamp2a), the receptor for chaperone-mediated autophagy, is one of the host proteins that interact with the 5′ UTR. Further studies showed that the Lamp2 level is upregulated in SARS-CoV-2-infected cells and that the absence of the Lamp2a isoform enhanced the viral RNA level whereas its overexpression significantly reduced the viral RNA level. Lamp2a and viral RNA colocalize in the infected cells, and there is an increased autophagic flux in infected cells, although there is no change in the formation of autophagolysosomes. In summary, our study provides a useful resource of SARS-CoV-2 5′ and 3′ UTR binding proteins and reveals the role of Lamp2a protein during SARS-CoV-2 infection.

IMPORTANCE Replication of a positive-strand RNA virus involves an RNA-protein complex consisting of viral genomic RNA, host RNA(s), virus-encoded proteins, and host proteins. Dissecting out individual components of the replication complex will help decode the mechanism of viral replication. 5′ and 3′ UTRs in positive-strand RNA viruses play essential regulatory roles in virus replication. Here, we identified the host proteins that associate with the UTRs of SARS-CoV-2, combined those data with the previously known protein-protein interaction data (expected to be involved in virus replication), and generated the RNA-protein-protein interaction (RPPI) network. Analysis of the RPPI network revealed the enrichment of factors involved in translation initiation and RNA metabolism, which are important for virus replication. Analysis of one of the interaction partners of the 5′-UTR (Lamp2a) demonstrated its role in reducing the viral RNA level in SARS-CoV-2-infected cells. Collectively, our study provides a resource of SARS-CoV-2 UTR-binding proteins and identifies an important role for host Lamp2a protein during viral infection.

Identification of hub genes associated with prognosis, diagnosis, immune infiltration and therapeutic drug in liver cancer by integrated analysis

BACKGROUND

Liver cancer is one of the most common cancers and causes of cancer death worldwide. The objective was to elucidate novel hub genes which were benefit for diagnosis, prognosis, and targeted therapy in liver cancer via integrated analysis.

METHODS

GSE84402, GSE101685, and GSE112791 were filtered from the Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) were identified by using the GEO2R. The GO and KEGG pathway of DEGs were analyzed in the DAVID. PPI and TF network of the DEGs were constructed by using the STRING, TRANSFAC, and Harmonizome. The relationship between hub genes and prognoses in liver cancer was analyzed in UALCAN based on The Cancer Genome Atlas (TCGA). The diagnostic value of hub genes was evaluated by ROC. The relationship between hub genes and tumor-infiltrate lymphocytes was analyzed in TIMER. The protein levels of hub genes were verified in HPA. The interaction between the hub genes and the drug were identified in DGIdb.

RESULTS

In total, 108 upregulated and 60 downregulated DEGs were enriched in 148 GO terms and 20 KEGG pathways. The mRNA levels and protein levels of CDK1, HMMR, PTTG1, and TTK were higher in liver cancer tissues compared to normal tissues, which showed excellent diagnostic and prognostic value. CDK1, HMMR, PTTG1, and TTK were positively correlated with tumor-infiltrate lymphocytes, which might involve tumor immune response. The CDK1, HMMR, and TTK had close interaction with anticancer agents.

CONCLUSIONS

The CDK1, HMMR, PTTG1, and TTK were hub genes in liver cancer; hence, they might be potential biomarkers for diagnosis, prognosis, and targeted therapy of liver cancer.

A Methodological Framework to Discover Pharmacogenomic Interactions Based on Random Forests

The identification of genomic alterations in tumor tissues, including somatic mutations, deletions, and gene amplifications, produces large amounts of data, which can be correlated with a diversity of therapeutic responses. We aimed to provide a methodological framework to discover pharmacogenomic interactions based on Random Forests. We matched two databases from the Cancer Cell Line Encyclopaedia (CCLE) project, and the Genomics of Drug Sensitivity in Cancer (GDSC) project. For a total of 648 shared cell lines, we considered 48,270 gene alterations from CCLE as input features and the area under the dose-response curve (AUC) for 265 drugs from GDSC as the outcomes. A three-step reduction to 501 alterations was performed, selecting known driver genes and excluding very frequent/infrequent alterations and redundant ones. For each model, we used the concordance correlation coefficient (CCC) for assessing the predictive performance, and permutation importance for assessing the contribution of each alteration. In a reasonable computational time (56 min), we identified 12 compounds whose response was at least fairly sensitive (CCC > 20) to the alteration profiles. Some diversities were found in the sets of influential alterations, providing clues to discover significant drug-gene interactions. The proposed methodological framework can be helpful for mining pharmacogenomic interactions.

Identification and validation of inferior prognostic genes associated with immune signatures and chemotherapy outcome in acute myeloid leukemia

Acute myeloid leukemia (AML) is a group of heterogeneous hematological malignancies. We identified key genes as ITGAM and lncRNA ITGB2-AS1 through different bioinformatics tools. Furthermore, qPCR was performed to verify the expression level of essential genes in clinical samples. Retrospective research on 179 AML cases was used to investigate the relationship between the expression of ITGAM and the characteristics of AML. The critical gene relationship with immune infiltration in AML was estimated. The clinical validation and prognostic investigation showed that ITGAM, PPBP, and ITGB2-AS1 are highly expressed in AML (P < 0.001) and significantly associated with the overall survival in AML. Moreover, the retrospective research on 179 clinical cases showed that positive expression of ITGAM is substantially related to AML classification (P < 0.001), higher count of white blood cells (P < 0.01), and poor chemotherapy outcome (P < 0.05). Furthermore, based on grouping ITGAM as the high and low expression in TCGA-LAML profile, we found that genes in the highly expressed ITGAM group are mainly involved in immune infiltration and inflammation-related signaling pathways. Finally, we discovered that the expression level of ITGAM and lncRNA ITGB2-AS1 are not just closely related to the immune score and stromal score (P < 0.001) but also significantly positively correlated with various Immune signatures in AML (P < 0.001), indicating the association of these genes with immunosuppression in AML. The prediction of candidate drugs indicated that certain immunosuppressive drugs have potential therapeutic effects for AML. The critical genes could be used as potential biomarkers to evaluate the survival and prognosis of AML.

The association between hypertension and nonalcoholic fatty liver disease (NAFLD): literature evidence and systems biology analysis

Nonalcoholic fatty liver disease (NAFLD) has become a major public health issue as its progression increases risks of multisystem morbidity and mortality. Recent evidence indicates a more complex relationship between hypertension and NAFLD than previously thought. In this study, a comprehensive literature search was used to gather information supporting the comorbidity phenomenon of hypertension and NAFLD. Then, systems biology approach was applied to identify the potential genes and mechanisms simultaneously associated with hypertension and NAFLD. With the help of protein-protein interaction network-based algorithm, we found that the distance between hypertension and NAFLD was much less than random ones. Sixty-four shared genes of hypertension and NAFLD modules were identified as core genes. Kyoto Encyclopedia of Genes and Genomes(KEGG) enrichment analysis indicated that some inflammatory, metabolic and endocrine signals were related to the potential biological functions of core genes. More importantly, drugs used to treat cardiovascular diseases, hypertension, hyperlipidemia, inflammatory diseases and depression could be potential therapeutics against hypertension-NAFLD co-occurrence. After analyzing public OMICs data, ALDH1A1 was identified as a potential therapeutic target, without being affected by reverse causality. These findings give a clue for the potential mechanisms of comorbidity of hypertension and NAFLD and highlight the multiple target-therapeutic strategy of NAFLD for future clinical research.

NIDM: network impulsive dynamics on multiplex biological network for disease-gene prediction

The prediction of genes related to diseases is important to the study of the diseases due to high cost and time consumption of biological experiments. Network propagation is a popular strategy for disease-gene prediction. However, existing methods focus on the stable solution of dynamics while ignoring the useful information hidden in the dynamical process, and it is still a challenge to make use of multiple types of physical/functional relationships between proteins/genes to effectively predict disease-related genes. Therefore, we proposed a framework of network impulsive dynamics on multiplex biological network (NIDM) to predict disease-related genes, along with four variants of NIDM models and four kinds of impulsive dynamical signatures (IDSs). NIDM is to identify disease-related genes by mining the dynamical responses of nodes to impulsive signals being exerted at specific nodes. By a series of experimental evaluations in various types of biological networks, we confirmed the advantage of multiplex network and the important roles of functional associations in disease-gene prediction, demonstrated superior performance of NIDM compared with four types of network-based algorithms and then gave the effective recommendations of NIDM models and IDS signatures. To facilitate the prioritization and analysis of (candidate) genes associated to specific diseases, we developed a user-friendly web server, which provides three kinds of filtering patterns for genes, network visualization, enrichment analysis and a wealth of external links (http://bioinformatics.csu.edu.cn/DGP/NID.jsp). NIDM is a protocol for disease-gene prediction integrating different types of biological networks, which may become a very useful computational tool for the study of disease-related genes.

RIFS2D: A two-dimensional version of a randomly restarted incremental feature selection algorithm with an application for detecting low-ranked biomarkers

The era of big data introduces both opportunities and challenges for biomedical researchers. One of the inherent difficulties in the biomedical research field is to recruit large cohorts of samples, while high-throughput biotechnologies may produce thousands or even millions of features for each sample. Researchers tend to evaluate the individual correlation of each feature with the class label and use the incremental feature selection (IFS) strategy to select the top-ranked features with the best prediction performance. Recent experimental data showed that a subset of continuously ranked features randomly restarted from a low-ranked feature (an RIFS block) may outperform the subset of top-ranked features. This study proposed a feature selection Algorithm RIFS2D by integrating multiple RIFS blocks. A comprehensive comparative experiment was conducted with the IFS, RIFS and existing feature selection algorithms and demonstrated that a subset of low-ranked features may also achieve promising prediction performance. This study suggested that a prediction model with promising performance may be trained by low-ranked features, even when top-ranked features did not achieve satisfying prediction performance. Further comparative experiments were conducted between RIFS2D and t-tests for the detection of early-stage breast cancer. The data showed that the RIFS2D-recommended features achieved better prediction accuracy and were targeted by more drugs than the t-test top-ranked features.

A COVID-19 Drug Repurposing Strategy through Quantitative Homological Similarities Using a Topological Data Analysis-Based Framework

Since its emergence in March 2020, the SARS-CoV-2 global pandemic has produced more than 116 million cases and 2.5 million deaths worldwide. Despite the enormous efforts carried out by the scientific community, no effective treatments have been developed to date. We applied a novel computational pipeline aimed to accelerate the process of identifying drug repurposing candidates which allows us to compare three-dimensional protein structures. Its use in conjunction with two in silico validation strategies (molecular docking and transcriptomic analyses) allowed us to identify a set of potential drug repurposing candidates targeting three viral proteins (3CL viral protease, NSP15 endoribonuclease, and NSP12 RNA-dependent RNA polymerase), which included rutin, dexamethasone, and vemurafenib. This is the first time that a topological data analysis (TDA)-based strategy has been used to compare a massive number of protein structures with the final objective of performing drug repurposing to treat SARS-CoV-2 infection.

Quantitative Proteomic Approach Reveals Altered Metabolic Pathways in Response to the Inhibition of Lysine Deacetylases in A549 Cells under Normoxia and Hypoxia

Growing evidence is showing that acetylation plays an essential role in cancer, but studies on the impact of KDAC inhibition (KDACi) on the metabolic profile are still in their infancy. Here, we analyzed, by using an iTRAQ-based quantitative proteomics approach, the changes in the proteome of KRAS-mutated non-small cell lung cancer (NSCLC) A549 cells in response to trichostatin-A (TSA) and nicotinamide (NAM) under normoxia and hypoxia. Part of this response was further validated by molecular and biochemical analyses and correlated with the proliferation rates, apoptotic cell death, and activation of ROS scavenging mechanisms in opposition to the ROS production. Despite the differences among the KDAC inhibitors, up-regulation of glycolysis, TCA cycle, oxidative phosphorylation and fatty acid synthesis emerged as a common metabolic response underlying KDACi. We also observed that some of the KDACi effects at metabolic levels are enhanced under hypoxia. Furthermore, we used a drug repositioning machine learning approach to list candidate metabolic therapeutic agents for KRAS mutated NSCLC. Together, these results allow us to better understand the metabolic regulations underlying KDACi in NSCLC, taking into account the microenvironment of tumors related to hypoxia, and bring new insights for the future rational design of new therapies.

Patient-Specific Network for Personalized Breast Cancer Therapy with Multi-Omics Data

The development of new computational approaches that are able to design the correct personalized drugs is the crucial therapeutic issue in cancer research. However, tumor heterogeneity is the main obstacle to developing patient-specific single drugs or combinations of drugs that already exist in clinics. In this study, we developed a computational approach that integrates copy number alteration, gene expression, and a protein interaction network of 73 basal breast cancer samples. 2509 prognostic genes harboring a copy number alteration were identified using survival analysis, and a protein-protein interaction network considering the direct interactions was created. Each patient was described by a specific combination of seven altered hub proteins that fully characterize the 73 basal breast cancer patients. We suggested the optimal combination therapy for each patient considering drug-protein interactions. Our approach is able to confirm well-known cancer related genes and suggest novel potential drug target genes. In conclusion, we presented a new computational approach in breast cancer to deal with the intra-tumor heterogeneity towards personalized cancer therapy.

Screening and identification of potential biomarkers and therapeutic drugs in melanoma via integrated bioinformatics analysis

Melanoma is a highly aggressive malignant skin tumor with a high rate of metastasis and mortality. In this study, a comprehensive bioinformatics analysis was used to clarify the hub genes and potential drugs. Download the GSE3189, GSE22301, and GSE35388 microarray datasets from the Gene Expression Omnibus (GEO), which contains a total of 33 normal samples and 67 melanoma samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) approach analyze DEGs based on the DAVID. Use STRING to construct protein-protein interaction network, and use MCODE and cytoHubba plug-ins in Cytoscape to perform module analysis and identified hub genes. Use Gene Expression Profile Interactive Analysis (GEPIA) to assess the prognosis of genes in tumors. Finally, use the Drug-Gene Interaction Database (DGIdb) to screen targeted drugs related to hub genes. A total of 140 overlapping DEGs were identified from the three microarray datasets, including 59 up-regulated DEGs and 81 down-regulated DEGs. GO enrichment analysis showed that these DEGs are mainly involved in the biological process such as positive regulation of gene expression, positive regulation of cell proliferation, positive regulation of MAP kinase activity, cell migration, and negative regulation of the apoptotic process. The cellular components are concentrated in the membrane, dendritic spine, the perinuclear region of cytoplasm, extracellular exosome, and membrane raft. Molecular functions include protein homodimerization activity, calmodulin-binding, transcription factor binding, protein binding, and cytoskeletal protein binding. KEGG pathway analysis shows that these DEGs are mainly related to protein digestion and absorption, PPAR signaling pathway, signaling pathways regulating stem cells' pluripotency, and Retinol metabolism. The 23 most closely related DEGs were identified from the PPI network and combined with the GEPIA prognostic analysis, CDH3, ESRP1, FGF2, GBP2, KCNN4, KIT, SEMA4D, and ZEB1 were selected as hub genes, which are considered to be associated with poor prognosis of melanoma closely related. Besides, ten related drugs that may have therapeutic effects on melanoma were also screened. These newly discovered genes and drugs provide new ideas for further research on melanoma.