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

Mendelian randomization of circulating proteome identifies actionable targets in heart failure

Background

Heart failure (HF) is a prevalent cause of mortality and morbidity. The molecular drivers of HF are still largely unknown.

Results

We aimed to identify circulating proteins causally associated with HF by leveraging genome-wide genetic association data for HF including 47,309 cases and 930,014 controls. We performed two-sample Mendelian randomization (MR) with multiple cis instruments as well as network and enrichment analysis using data from blood protein quantitative trait loci (pQTL) (2,965 blood proteins) measured in 3,301 individuals. Nineteen blood proteins were causally associated with HF, were not subject to reverse causality and were enriched in ligand-receptor and glycosylation molecules. Network pathway analysis of the blood proteins showed enrichment in NF-kappa B, TGF beta, lipid in atherosclerosis and fluid shear stress. Cross-phenotype analysis of HF identified genetic overlap with cardiovascular drugs, myocardial infarction, parental longevity and low-density cholesterol. Multi-trait MR identified causal associations between HF-associated blood proteins and cardiovascular outcomes. Multivariable MR showed that association of BAG3, MIF and APOA5 with HF were mediated by the blood pressure and coronary artery disease. According to the directional effect and biological action, 7 blood proteins are targets of existing drugs or are tractable for the development of novel therapeutics. Among the pathways, sialyl Lewis x and the activin type II receptor are potential druggable candidates.

Conclusions

Integrative MR analyses of the blood proteins identified causally-associated proteins with HF and revealed pleiotropy of the blood proteome with cardiovascular risk factors. Some of the proteins or pathway related mechanisms could be targeted as novel treatment approach in HF.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08811-2.

CCAS: One-stop and comprehensive annotation system for individual cancer genome at multi-omics level

Due to the explosion of cancer genome data and the urgent needs for cancer treatment, it is becoming increasingly important and necessary to easily and timely analyze and annotate cancer genomes. However, tumor heterogeneity is recognized as a serious barrier to annotate cancer genomes at the individual patient level. In addition, the interpretation and analysis of cancer multi-omics data rely heavily on existing database resources that are often located in different data centers or research institutions, which poses a huge challenge for data parsing. Here we present CCAS (Cancer genome Consensus Annotation System, https://ngdc.cncb.ac.cn/ccas/#/home), a one-stop and comprehensive annotation system for the individual patient at multi-omics level. CCAS integrates 20 widely recognized resources in the field to support data annotation of 10 categories of cancers covering 395 subtypes. Data from each resource are manually curated and standardized by using ontology frameworks. CCAS accepts data on single nucleotide variant/insertion or deletion, expression, copy number variation, and methylation level as input files to build a consensus annotation. Outputs are arranged in the forms of tables or figures and can be searched, sorted, and downloaded. Expanded panels with additional information are used for conciseness, and most figures are interactive to show additional information. Moreover, CCAS offers multidimensional annotation information, including mutation signature pattern, gene set enrichment analysis, pathways and clinical trial related information. These are helpful for intuitively understanding the molecular mechanisms of tumors and discovering key functional genes.

Identification of pathogenic genes associated with CKD: An integrated bioinformatics approach

Chronic kidney disease (CKD) is defined as a persistent abnormality in the structure and function of kidneys and leads to high morbidity and mortality in individuals across the world. Globally, approximately 8%–16% of the population is affected by CKD. Proper screening, staging, diagnosis, and the appropriate management of CKD by primary care clinicians are essential in preventing the adverse outcomes associated with CKD worldwide. In light of this, the identification of biomarkers for the appropriate management of CKD is urgently required. Growing evidence has suggested the role of mRNAs and microRNAs in CKD, however, the gene expression profile of CKD is presently uncertain. The present study aimed to identify diagnostic biomarkers and therapeutic targets for patients with CKD. The human microarray profile datasets, consisting of normal samples and treated samples were analyzed thoroughly to unveil the differentially expressed genes (DEGs). After selection, the interrelationship among DEGs was carried out to identify the overlapping DEGs, which were visualized using the Cytoscape program. Furthermore, the PPI network was constructed from the String database using the selected DEGs. Then, from the PPI network, significant modules and sub-networks were extracted by applying the different centralities methods (closeness, betweenness, stress, etc.) using MCODE, Cytohubba, and Centiserver. After sub-network analysis we identified six overlapped hub genes (RPS5, RPL37A, RPLP0, CXCL8, HLA-A, and ANXA1). Additionally, the enrichment analysis was undertaken on hub genes to determine their significant functions. Furthermore, these six genes were used to find their associated miRNAs and targeted drugs. Finally, two genes CXCL8 and HLA-A were common for Ribavirin drug (the gene-drug interaction), after docking studies HLA-A was selected for further investigation. To conclude our findings, we can say that the identified hub genes and their related miRNAs can serve as potential diagnostic biomarkers and therapeutic targets for CKD treatment strategies.

PlaqView 2.0: A comprehensive web portal for cardiovascular single-cell genomics

Single-cell RNA-seq (scRNA-seq) is a powerful genomics technology to interrogate the cellular composition and behaviors of complex systems. While the number of scRNA-seq datasets and available computational analysis tools have grown exponentially, there are limited systematic data sharing strategies to allow rapid exploration and re-analysis of single-cell datasets, particularly in the cardiovascular field. We previously introduced PlaqView, an open-source web portal for the exploration and analysis of published atherosclerosis single-cell datasets. Now, we introduce PlaqView 2.0 (www.plaqview.com), which provides expanded features and functionalities as well as additional cardiovascular single-cell datasets. We showcase improved PlaqView functionality, backend data processing, user-interface, and capacity. PlaqView brings new or improved tools to explore scRNA-seq data, including gene query, metadata browser, cell identity prediction, ad hoc RNA-trajectory analysis, and drug-gene interaction prediction. PlaqView serves as one of the largest central repositories for cardiovascular single-cell datasets, which now includes data from human aortic aneurysm, gene-specific mouse knockouts, and healthy references. PlaqView 2.0 brings advanced tools and high-performance computing directly to users without the need for any programming knowledge. Lastly, we outline steps to generalize and repurpose PlaqView's framework for single-cell datasets from other fields.

Transcriptomic data exploration of consensus genes and molecular mechanisms between chronic obstructive pulmonary disease and lung adenocarcinoma

Most current research has focused on chronic obstructive pulmonary disease (COPD) and lung adenocarcinoma (LUAD) alone; however, it is important to understand the complex mechanism of COPD progression to LUAD. This study is the first to explore the unique and jointly molecular mechanisms in the pathogenesis of COPD and LUAD across several datasets based on a variety of analysis methods. We used weighted correlation network analysis to search hub genes in two datasets from public databases: GSE10072 and GSE76925. We explored the unique and jointly molecular mechanistic signatures of the two diseases in pathogenesis through enrichment analysis, immune infiltration analysis, and therapeutic targets analysis. Finally, the results were confirmed using real-time quantitative reverse transcription PCR. Fifteen hub genes were identified: GPI, EZH2, EFNA4, CFB, ENO1, SH3PXD2B, SELL, CORIN, MAD2L1, CENPF, TOP2A, ASPM, IGFBP2, CDKN2A, and ELF3. For the first time, SELL, CORIN, GPI, and EFNA4 were found to play a role in the etiology of COPD and LUAD. The LUAD genes identified were primarily involved in the cell cycle and DNA replication processes; COPD genes we found were related to ubiquitin-mediated proteolysis, ribosome, and T/B-cell receptor signaling pathways. The tumor microenvironment of LUAD pathogenesis was influenced by CD4 + T cells, type 1 regulatory T cells, and T helper 1 cells. T follicular helper cells, natural killer T cells, and B cells all impact the immunological inflammation in COPD. The results of drug targets analysis suggest that cisplatin and tretinoin, as well as bortezomib and metformin may be potential targeted therapy for patients with COPD combined LUAD. These signatures may be provided a new direction for developing early interventions and treatments to improve the prognosis of COPD and LUAD.

In silico Analysis of ACE2 Receptor to Find Potential Herbal Drugs in COVID-19 Associated Neurological Dysfunctions

COVID-19 mainly causes the collapse of the pulmonary system thereby causing a dearth of oxygen in the human body. Patients infected with this viral disease have been reported to experience various signs and symptoms associated with brain dysfunction, from the feeling of vagueness to loss of smell and taste to severe strokes. These neurological problems have been reported by younger COVID-19 infected patients mainly in their thirties and forties. Various researchers from around the globe have discerned numerous other brain dysfunctions, such as headache, dizziness, numbness, major depressive disorder, anosmia, encephalitis, febrile seizures, and Guillain-Barre syndrome. The involvement of the CNS by this viral infection has been predicted to be for a longer period of time, even if the patient recovers from COVID-19. The neuronal cell damage caused by COVID-19 is a potent factor responsible for cognitive, behavioral, and psychological problems among its sufferers. The hypoxic conditions can also trigger the formation of beta-amyloid plaques and tau-tangles and thus the virus can even induce Alzheimer’s in patients in the near future. The virus affects the brain directly, thereby causing encephalitis. This pandemic has also been shown to have a negative psychological toll on people. This research aims to highlight the brain dysfunction associated with the ACE2 receptor that is known to be a crucial player in the COVID-19 pandemic using genetic networking approaches. Furthermore, we have identified herbal drug candidates that bind to the ACE2 receptor in order to identify potential treatments for the neurological manifestations of COVID-19.

Exploration and validation of metastasis-associated genes for skin cutaneous melanoma

Skin cutaneous melanoma is a malignant and highly metastatic skin tumor, and its morbidity and mortality are still rising worldwide. However, the molecular mechanisms that promote melanoma metastasis are unclear. Two datasets (GSE15605 and GSE46517) were retrieved to identify the differentially expressed genes (DEGs), including 23 normal skin tissues (N), 77 primary melanoma tissues (T) and 85 metastatic melanoma tissues (M). Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed to explore the functions of the DEGs. We constructed protein–protein interaction network using the STRING database and Cytoscape software. Using the cytoHubba plugin of Cytoscape, we identified the most significant hub genes by five analytical methods (Degree, Bottleneck, MCC, MNC, and EPC). Hub gene expression was validated using the UALCAN website. Clinical relevance was investigated using The Cancer Genome Atlas resources. Finally, we explored the association between metastasis-associated genes and immune infiltrates through the Tumor Immune Estimation Resource (TIMER) database and performed drug–gene interaction analysis using the Drug-Gene Interaction database. A total of 294 specific genes were related to melanoma metastasis and were mainly involved in the positive regulation of locomotion, mitotic cell cycle process, and epithelial cell differentiation. Four hub genes (CDK1, FOXM1, KIF11, and RFC4) were identified from the cytoHubba plugin of Cytoscape. CDK1 was significantly upregulated in metastatic melanoma compared with primary melanoma, and high CDK1 expression was positively correlated with worse overall survival. Immune infiltration analysis revealed that CDK1 expression negatively correlated with macrophage infiltration (Rho = − 0.164, P = 2.02e−03) and positively correlated with neutrophil cells (Rho = 0.269, P = 2.72e−07) in SKCM metastasis. In addition, we identified that CDK1 had a close interaction with 10 antitumor drugs. CDK1 was identified as a hub gene involved in the progression of melanoma metastasis and may be regarded as a therapeutic target for melanoma patients to improve prognosis and prevent metastasis in the future.

Identification of risk genes for Alzheimer's disease by gene embedding

Most disease-gene association methods do not account for gene-gene interactions, even though these play a crucial role in complex, polygenic diseases like Alzheimer's disease (AD). To discover new genes whose interactions may contribute to pathology, we introduce GeneEMBED. This approach compares the functional perturbations induced in gene interaction network neighborhoods by coding variants from disease versus healthy subjects. In two independent AD cohorts of 5,169 exomes and 969 genomes, GeneEMBED identified novel candidates. These genes were differentially expressed in post mortem AD brains and modulated neurological phenotypes in mice. Four that were differentially overexpressed and modified neurodegeneration in vivo are PLEC, UTRN, TP53, and POLD1. Notably, TP53 and POLD1 are involved in DNA break repair and inhibited by approved drugs. While these data show proof of concept in AD, GeneEMBED is a general approach that should be broadly applicable to identify genes relevant to risk mechanisms and therapy of other complex diseases.

Mutational Signature and Integrative Genomic Analysis of Human Papillomavirus-Associated Penile Squamous Cell Carcinomas from Latin American Patients

High-throughput DNA sequencing has allowed for the identification of genomic alterations and their impact on tumor development, progression, and therapeutic responses. In PSCC, for which the incidence has progressively increased worldwide, there are still limited data on the molecular mechanisms involved in the disease pathogenesis. In this study, we characterized the mutational signature of 30 human papillomavirus (HPV)-associated PSCC cases from Latin Americans, using whole-exome sequencing. Copy number variations (CNVs) were also identified and compared to previous array-generated data. Enrichment analyses were performed to reveal disrupted pathways and to identify alterations mapped to HPV integration sites (HPVis) and miRNA-mRNA hybridization regions. Among the most frequently mutated genes were NOTCH1, TERT, TTN, FAT1, TP53, CDKN2A, RYR2, CASP8, FBXW7, HMCN2, and ITGA8. Of note, 92% of these altered genes were localized at HPVis. We also found mutations in ten novel genes (KMT2C, SMARCA4, PTPRB, AJUBA, CR1, KMT2D, NBEA, FAM135B, GTF2I, and CIC), thus increasing our understanding of the potential HPV-disrupted pathways. Therefore, our study reveals innovative targets with potential therapeutic benefits for HPV-associated PSCCs. The CNV analysis by sequencing (CNV-seq) revealed five cancer-associated genes as the most frequent with gains (NOTCH1, MYC, NUMA1, PLAG1, and RAD21), while 30% of the tumors showed SMARCA4 with loss. Additionally, four cancer-associated genes (CARD11, CSMD3, KDR, and TLX3) carried untranslated regions (UTRs) variants, which may impact gene regulation by affecting the miRNAs hybridization regions. Altogether, these data contribute to the characterization of the mutational spectrum and its impact on cellular signaling pathways in PSCC, thus reinforcing the pivotal role of HPV infection in the molecular pathogenesis of these tumors.

Immune infiltration and clinical significance analyses of the coagulation-related genes in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common types of cancers and a global health challenge with a low early diagnosis rate and high mortality. The coagulation cascade plays an important role in the tumor immune microenvironment (TME) of HCC. In this study, based on the coagulation pathways collected from the KEGG database, two coagulation-related subtypes were distinguished in HCC patients. We demonstrated the distinct differences in immune characteristics and prognostic stratification between two coagulation-related subtypes. A coagulation-related risk score prognostic model was developed in the Cancer Genome Atlas (TCGA) cohort for risk stratification and prognosis prediction. The predictive values of the coagulation-related risk score in prognosis and immunotherapy were also verified in the TCGA and International Cancer Genome Consortium cohorts. A nomogram was also established to facilitate the clinical use of this risk score and verified its effectiveness using different approaches. Based on these results, we can conclude that there is an obvious correlation between the coagulation and the TME in HCC, and the risk score could serve as a robust prognostic biomarker, provide therapeutic benefits for chemotherapy and immunotherapy and may be helpful for clinical decision making in HCC patients.

Repurposing Drugs via Network Analysis: Opportunities for Psychiatric Disorders

Despite advances in pharmacology and neuroscience, the path to new medications for psychiatric disorders largely remains stagnated. Drug repurposing offers a more efficient pathway compared with de novo drug discovery with lower cost and less risk. Various computational approaches have been applied to mine the vast amount of biomedical data generated over recent decades. Among these methods, network-based drug repurposing stands out as a potent tool for the comprehension of multiple domains of knowledge considering the interactions or associations of various factors. Aligned well with the poly-pharmacology paradigm shift in drug discovery, network-based approaches offer great opportunities to discover repurposing candidates for complex psychiatric disorders. In this review, we present the potential of network-based drug repurposing in psychiatry focusing on the incentives for using network-centric repurposing, major network-based repurposing strategies and data resources, applications in psychiatry and challenges of network-based drug repurposing. This review aims to provide readers with an update on network-based drug repurposing in psychiatry. We expect the repurposing approach to become a pivotal tool in the coming years to battle debilitating psychiatric disorders.

SARS-CoV-2 potential drugs, drug targets, and biomarkers: a viral-host interaction network-based analysis

COVID-19 is a global pandemic impacting the daily living of millions. As variants of the virus evolve, a complete comprehension of the disease and drug targets becomes a decisive duty. The Omicron variant, for example, has a notably high transmission rate verified in 155 countries. We performed integrative transcriptomic and network analyses to identify drug targets and diagnostic biomarkers and repurpose FDA-approved drugs for SARS-CoV-2. Upon the enrichment of 464 differentially expressed genes, pathways regulating the host cell cycle were significant. Regulatory and interaction networks featured hsa-mir-93-5p and hsa-mir-17-5p as blood biomarkers while hsa-mir-15b-5p as an antiviral agent. MYB, RRM2, ERG, CENPF, CIT, and TOP2A are potential drug targets for treatment. HMOX1 is suggested as a prognostic biomarker. Enhancing HMOX1 expression by neem plant extract might be a therapeutic alternative. We constructed a drug-gene network for FDA-approved drugs to be repurposed against the infection. The key drugs retrieved were members of anthracyclines, mitotic inhibitors, anti-tumor antibiotics, and CDK1 inhibitors. Additionally, hydroxyquinone and digitoxin are potent TOP2A inhibitors. Hydroxyurea, cytarabine, gemcitabine, sotalol, and amiodarone can also be redirected against COVID-19. The analysis enforced the repositioning of fluorouracil and doxorubicin, especially that they have multiple drug targets, hence less probability of resistance.

Comprehensive Assessment of Indian Variations in the Druggable Kinome Landscape Highlights Distinct Insights at the Sequence, Structure and Pharmacogenomic Stratum

India confines more than 17% of the world’s population and has a diverse genetic makeup with several clinically relevant rare mutations belonging to many sub-group which are undervalued in global sequencing datasets like the 1000 Genome data (1KG) containing limited samples for Indian ethnicity. Such databases are critical for the pharmaceutical and drug development industry where diversity plays a crucial role in identifying genetic disposition towards adverse drug reactions. A qualitative and comparative sequence and structural study utilizing variant information present in the recently published, largest curated Indian genome database (IndiGen) and the 1000 Genome data was performed for variants belonging to the kinase coding genes, the second most targeted group of drug targets. The sequence-level analysis identified similarities and differences among different populations based on the nsSNVs and amino acid exchange frequencies whereas a comparative structural analysis of IndiGen variants was performed with pathogenic variants reported in UniProtKB Humsavar data. The influence of these variations on structural features of the protein, such as structural stability, solvent accessibility, hydrophobicity, and the hydrogen-bond network was investigated. In-silico screening of the known drugs to these Indian variation-containing proteins reveals critical differences imparted in the strength of binding due to the variations present in the Indian population. In conclusion, this study constitutes a comprehensive investigation into the understanding of common variations present in the second largest population in the world and investigating its implications in the sequence, structural and pharmacogenomic landscape. The preliminary investigation reported in this paper, supporting the screening and detection of ADRs specific to the Indian population could aid in the development of techniques for pre-clinical and post-market screening of drug-related adverse events in the Indian population.

Identify Multiple Gene-Drug Common Modules Via Constrained Graph Matching

Identifying gene-drug interactions is vital to understanding biological mechanisms and achieving precise drug repurposing. High-throughput technologies produce a large amount of pharmacological and genomic data, providing an opportunity to explore the associations between oncogenic genes and therapeutic drugs. However, most studies only focus on "one-to-one" or "one-to-many" interactions, ignoring the multivariate patterns between genes and drugs. In this article, a high-order graph matching model with hypergraph constraints is proposed to discover the gene-drug common regulatory modules. Moreover, the prior knowledge is formulated into hypergraph constraints to reveal their multiple correspondences, penalizing the tensor matching process. The experimental results on the synthetic data demonstrate the proposed model is robust to noise contamination and outlier corruption, achieving a better performance than four state-of-the-art methods. We then evaluate the statistical power of our proposed method on the pharmacogenomics data. Our identified gene-drug common modules not only show significantly enriched pathways associated with cancer but also manifest the highly close gene-drug interactions.

Identification and Validation of Autophagy-Related Genes in Primary Ovarian Insufficiency by Gene Expression Profile and Bioinformatic Analysis

Background

To investigate the relationship between primary ovarian insufficiency and autophagy, we detected and got the expression profile of human granulosa cell line SVOG, which was with or without LPS induced. The expression profile was analyzed with the focus on the autophagy genes, among which hub genes were identified.

Results

Totally, 6 genes were selected as candidate hub genes which might correlate with the process of primary ovarian insufficiency. The expression of hub genes was then validated by quantitative real-time PCR and two of them had significant expression change. Bioinformatics analysis was performed to observe the features of hub genes, including hub gene-RBP/TF/miRNA/drug network construction, functional analysis, and protein-protein interaction network. Pearson's correlation analysis was also performed to identify the correlation between hub genes and autophagy genes, among which there were four autophagy genes significantly correlated with hub genes, including ATG4B, ATG3, ATG13, and ULK1.

Conclusion

The results indicated that autophagy might play an essential role in the process and underlying molecular mechanism of primary ovarian insufficiency, which was revealed for the first time and may help to provide a molecular foundation for the development of diagnostic and therapeutic approaches for primary ovarian insufficiency.

Exploring the fuzzy border between senolytics and senomorphics with chemoinformatics and systems pharmacology

Senescent cells accumulate within tissues during aging and secrete an array of pro-inflammatory molecules known as senescent-associated secretory phenotype (SASP), which contribute to the appearance and progression of various chronic degenerative diseases. Novel pharmacological approaches aimed at modulating or eliminating senescent cells´ harmful effects have recently emerged: Senolytics are molecules that selectively eliminate senescent cells, while senomorphics modulate or decrease the inflammatory response to specific SASP. So far, the physicochemical, structural, and pharmacological properties that define these two kinds of pharmacological approaches remain unclear. Therefore, the identification and correct choice of molecules, based on their physicochemical, structural, and pharmacological properties, likely to exhibit the desired senotherapeutic activity is crucial for developing effective, selective, and safe senotherapies. Here we compared the physicochemical, structural, and pharmacological properties of 84 senolytics and 79 senomorphics using a chemoinformatic and systems pharmacology approach. We found great physicochemical, structural, and pharmacological similarities between them, also reflected in their cellular responses measured through transcriptome perturbations. The identified similarities between senolytics and senomorphics might explain the dual activity of some of those molecules. These findings will help design and discover new, more effective, and highly selective senotherapeutic agents.

Global analyses of mRNA expression in human sensory neurons reveal eIF5A as a conserved target for inflammatory pain

Nociceptors are a type of sensory neuron that are integral to most forms of pain. Targeted disruption of nociceptor sensitization affords unique opportunities to prevent pain. An emerging model for nociceptors are sensory neurons derived from human stem cells. Here, we subjected five groups to high-throughput sequencing: human induced pluripotent stem cells (hiPSCs) prior to differentiation, mature hiPSC-derived sensory neurons, mature co-cultures containing hiPSC-derived astrocytes and sensory neurons, mouse dorsal root ganglion (DRG) tissues, and mouse DRG cultures. Co-culture of nociceptors and astrocytes promotes expression of transcripts enriched in DRG tissues. Comparisons of the hiPSC models to tissue samples reveal that many key transcripts linked to pain are present. Markers indicative of a range of neuronal subtypes present in the DRG were detected in mature hiPSCs. Intriguingly, translation factors were maintained at consistently high expression levels across species and culture systems. As a proof of concept for the utility of this resource, we validated expression of eukaryotic initiation factor 5A (eIF5A) in DRG tissues and hiPSC samples. eIF5A is subject to a unique posttranslational hypusine modification required for its activity. Inhibition of hypusine biosynthesis prevented hyperalgesic priming by inflammatory mediators in vivo and diminished hiPSC activity in vitro. Collectively, our results illuminate the transcriptomes of hiPSC sensory neuron models. We provide a demonstration for this resource through our investigation of eIF5A. Our findings reveal hypusine as a potential target for inflammation associated pain in males.

Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease

More than 800 million people suffer from kidney disease, yet the mechanism of kidney dysfunction is poorly understood. In the present study, we define the genetic association with kidney function in 1.5 million individuals and identify 878 (126 new) loci. We map the genotype effect on the methylome in 443 kidneys, transcriptome in 686 samples and single-cell open chromatin in 57,229 kidney cells. Heritability analysis reveals that methylation variation explains a larger fraction of heritability than gene expression. We present a multi-stage prioritization strategy and prioritize target genes for 87% of kidney function loci. We highlight key roles of proximal tubules and metabolism in kidney function regulation. Furthermore, the causal role of SLC47A1 in kidney disease is defined in mice with genetic loss of Slc47a1 and in human individuals carrying loss-of-function variants. Our findings emphasize the key role of bulk and single-cell epigenomic information in translating genome-wide association studies into identifying causal genes, cellular origins and mechanisms of complex traits.

Inflammation subtypes in psychosis and their relationships with genetic risk for psychiatric and cardiometabolic disorders

Cardiometabolic disorders have known inflammatory implications, and peripheral measures of inflammation and cardiometabolic disorders are common in persons with psychotic disorders. Inflammatory signatures are also related to neurobiological and behavioral changes in psychosis. Relationships between systemic inflammation and cardiometabolic genetic risk in persons with psychosis have not been examined. Thirteen peripheral inflammatory markers and genome-wide genotyping were assessed in 122 participants (n ​= ​86 psychosis, n ​= ​36 healthy controls) of European ancestry. Cluster analyses of inflammatory markers classified higher and lower inflammation subgroups. Single-trait genetic risk scores (GRS) were constructed for each participant using previously reported GWAS summary statistics for the following traits: schizophrenia, bipolar disorder, major depressive disorder, coronary artery disease, type-2 diabetes, low-density lipoprotein, high-density lipoprotein, triglycerides, and waist-to-hip ratio. Genetic correlations across traits were quantified. Principal component (PC) analysis of the cardiometabolic GRSs generated six PC loadings used in regression models to examine associations with inflammation markers. Functional module discovery explored biological mechanisms of the inflammation association of cardiometabolic GRS genes. A subgroup of 38% persons with psychotic disorders was characterized with higher inflammation status. These higher inflammation individuals had lower BACS scores (p ​= ​0.038) compared to those with lower inflammation. The first PC of the cardiometabolic GRS matrix was related to higher inflammation status in persons with psychotic disorders (OR ​= ​2.037, p ​= ​0.001). Two of eight modules within the functional interaction network of cardiometabolic GRS genes were enriched for immune processes. Cardiometabolic genetic risk may predispose some individuals with psychosis to elevated inflammation which adversely impacts cognition associated with illness.

Metabolic modeling-based drug repurposing in Glioblastoma

The manifestation of intra- and inter-tumor heterogeneity hinders the development of ubiquitous cancer treatments, thus requiring a tailored therapy for each cancer type. Specifically, the reprogramming of cellular metabolism has been identified as a source of potential drug targets. Drug discovery is a long and resource-demanding process aiming at identifying and testing compounds early in the drug development pipeline. While drug repurposing efforts (i.e., inspecting readily available approved drugs) can be supported by a mechanistic rationale, strategies to further reduce and prioritize the list of potential candidates are still needed to facilitate feasible studies. Although a variety of ‘omics’ data are widely gathered, a standard integration method with modeling approaches is lacking. For instance, flux balance analysis is a metabolic modeling technique that mainly relies on the stoichiometry of the metabolic network. However, exploring the network’s topology typically neglects biologically relevant information. Here we introduce Transcriptomics-Informed Stoichiometric Modelling And Network analysis (TISMAN) in a recombinant innovation manner, allowing identification and validation of genes as targets for drug repurposing using glioblastoma as an exemplar.

Esophageal squamous cell carcinoma: Integrated bioinformatics analysis for differential gene expression with identification of hub genes and lncRNA

Background

Esophageal squamous cell carcinoma (ESCC) is a typical Gastro-Intestinal (GI) tract neoplasm. This study was conducted to know the Differential Expressed Genes (DEGs) profile of ESCC along with hub gene screening, lncRNA identification, and drug-genes interactions.

Methods

GSE161533, GSE20347, GSE45670 microarray datasets were retrieved from the NCBI Gene Expression Omnibus (GEO) database. GEO2R was used for the DEGs identification, whereas GO (Gene Ontology) and KEGG enrichment analysis were performed in DAVID. PPI network constructed using STRING and visualized with Cytoscape app with the help of MCODE. The top ten connectivity genes were selected as hub genes—further survival analysis was performed in the Kaplan-Meier plotter. Moreover, Boxplot, pathological stage plots were constructed using GEPIA (Gene Expression Profiling Interactive Analysis). The methylation heatmap assembled in the DiseaseMeth version 2.0. lncRNA (Long non-coding RNA) was identified comparing the list of genes in HUGO, and Gene-drug interactions were accumulated from the DgiDB platform.

Results

This experiment showed 16 upregulated, and 59 downregulated DEGs shared among the three datasets. Biological process analysis showed significant terms such as extracellular matrix disassembly and collagen catabolism. The extracellular region was detected as the most crucial cellular compartment. Notably, metalloen dopeptidease and serine-type endopeptidase activity showed significant molecular functions term. In contrast, transcriptional misregulation was a highly substantial KEGG pathway. Kaplan-Meier plotter showed higher expression of CXCL8, SPP1, MMP13, CXCL1, and TOP2A have a significant impact on the overall survival of the patients. Nine out of ten hub genes have significantly different expression levels than normal and cancer tissues. HYMAI was the only lncRNA commonly expressed upregulated among the three datasets. Drug-gene interaction showed multiple genes have no drug options exist till now.

•

GSE161533, GSE20347, and GSE45670 microarray datasets were analyzed.

•

16 upregulated and 59 downregulated DEGs shared among the three datasets.

•

CXCL8, SPP1, MMP13, CXCL1, and TOP2A have a significant impact on survival.

•

HYMAI was the only lncRNA commonly expressed.

•

Multiple genes have no drug options that exist.

Identification of hub genes and regulatory networks in histologically unstable carotid atherosclerotic plaque by bioinformatics analysis

Objective

This study identified underlying genetic molecules associated with histologically unstable carotid atherosclerotic plaques through bioinformatics analysis that may be potential biomarkers and therapeutic targets.

Methods

Three transcriptome datasets (GSE41571, GSE120521 and E-MTAB-2055) and one non-coding RNA dataset (GSE111794) that met histological grouping criteria of unstable plaque were downloaded. The common differentially expressed genes (co-DEGs) of unstable plaques identified from three mRNA datasets were annotated by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomics (KEGG). A protein–protein interaction (PPI) network was constructed to present the interaction between co-DEGs and screen out hub genes. MiRNet database and GSE111794 dataset were used to identify the miRNAs targeting hub genes. Associated transcription factors (TFs) and drugs were also predicted. These predicted results were used to construct miRNA/TFs-hub gene and drug-hub gene regulatory networks.

Results

A total of 105 co-DEGs were identified, including 42 up-regulated genes and 63 down-regulated genes, which were mainly enriched in collagen-containing extracellular matrix, focal adhesion, actin filament bundle, chemokine signaling pathway and regulates of actin cytoskeleton. Ten hub genes (up-regulated: HCK, C1QC, CD14, FCER1G, LCP1 and RAC2; down-regulated: TPM1, MYH10, PLS3 and FMOD) were screened. HCK and RAC2 were involved in chemokine signaling pathway, MYH10 and RAC2 were involved in regulation of actin cytoskeleton. We also predicted 12 miRNAs, top5 TFs and 25 drugs targeting hub genes. In the miRNA/TF-hub gene regulatory network, PLS3 was the most connected hub genes and was targeted by six miRNAs and all five screened TFs. In the drug-hub gene regulatory network, HCK was targeted by 20 drugs including 10 inhibitors.

Conclusions

We screened 10 hub genes and predicted miRNAs and TFs targeting them. These molecules may play a crucial role in the progression of histologically unstable carotid plaques and serve as potential biomarkers and therapeutic targets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01257-1.

urPTMdb/TeaProt: Upstream and Downstream Proteomics Analysis

We have developed the underrepresented post-translational modification (PTM) database (urPTMdb), a PTM gene set database to accelerate the discovery of enriched protein modifications in experimental data. urPTMdb provides curated lists of proteins reported to be substrates of underrepresented modifications. Their enrichment in proteomics datasets can reveal unexpected PTM regulations. urPTMdb can be implemented in existing workflows, or used in TeaProt, an online Shiny tool that integrates upstream transcription factor enrichment analysis with downstream pathway analysis through an easy-to-use interactive interface. TeaProt annotates user-uploaded data with drug-gene interactions, subcellular localizations, phenotypic functions, gene-disease associations, and enzyme-gene interactions. TeaProt enables gene set enrichment analysis (GSEA) to discover enrichments in gene sets from various resources, including MSigDB, CHEA, and urPTMdb. We demonstrate the utility of urPTMdb and TeaProt through the analysis of a previously published Western diet-induced remodeling of the tongue proteome, which revealed altered cellular processes associated with energy metabolism, interferon alpha/gamma response, adipogenesis, HMGylation substrate enrichment, and transcription regulation through PPARG and CEBPA. Additionally, we analyzed the interactome of ADP-ribose glycohydrolase TARG1, a key enzyme that removes mono-ADP-ribosylation. This analysis identified an enrichment of ADP-ribosylation, ribosomal proteins, and proteins localized in the nucleoli and endoplasmic reticulum. TeaProt and urPTMdb are accessible at https://tea.coffeeprot.com/.

A Knowledge Graph Embedding Based Approach to Predict the Adverse Drug Reactions Using a Deep Neural Network

Recently Artificial Intelligence(AI) has not only been used to diagnose the disease but also to cure the disease. Researchers started using AI for drug discovery. Predicting the Adverse Drug Reactions(ADRs) caused by the drug in the manufacturing stage or in the clinical trial stage is a very important problem in drug discovery. ADRs have become a major concern resulting in injuries and also becoming fatal sometimes. Drug safety has gained much importance over the years propelling to the forefront investigation of predicting the ADRs. Although prior studies have queried diverse approaches to predict ADRs, very few were found to be effective. Also, the problem of having fewer reports makes the prediction of ADRs more difficult. To tackle this problem effectively, a novel method has been proposed in this paper. The proposed method is based on Knowledge Graph(KG) embedding. Using the KG embedding, we designed and trained a custom-made Deep Neural Network(DNN) called KGDNN(Knowledge Graph DNN) for predicting the ADRs. A KG has been constructed with 6 types of entities: drugs, ADRs, target proteins, indications, pathways, and genes. Using the Node2Vec algorithm, each node has been embedded into a feature space. Using those embeddings, the ADRs are classified by the KGDNN model. The proposed method has obtained an AUROC score of 0.917 and significantly outperformed the existing methods. Two case studies on drugs causing liver injury and COVID-19 recommended drugs have been performed to illustrate the model efficacy.

Comprehensive bioinformatic analysis reveals a cancer-associated fibroblast gene signature as a poor prognostic factor and potential therapeutic target in gastric cancer

BACKGROUND

Gastric cancer is one of the deadliest cancers, currently available therapies have limited success. Cancer-associated fibroblasts (CAFs) are pivotal cells in the stroma of gastric tumors posing a great risk for progression and chemoresistance. The poor prognostic signature for CAFs is not clear in gastric cancer, and drugs that target CAFs are lacking in the clinic. In this study, we aim to identify a poor prognostic gene signature for CAFs, targeting which may increase the therapeutic success in gastric cancer.

METHODS

We analyzed four GEO datasets with a network-based approach and validated key CAF markers in The Cancer Genome Atlas (TCGA) and The Asian Cancer Research Group (ACRG) cohorts. We implemented stepwise multivariate Cox regression guided by a pan-cancer analysis in TCGA to identify a poor prognostic gene signature for CAF infiltration in gastric cancer. Lastly, we conducted a database search for drugs targeting the signature genes.

RESULTS

Our study revealed the COL1A1, COL1A2, COL3A1, COL5A1, FN1, and SPARC as the key CAF markers in gastric cancer. Analysis of the TCGA and ACRG cohorts validated their upregulation and poor prognostic significance. The stepwise multivariate Cox regression elucidated COL1A1 and COL5A1, together with ITGA4, Emilin1, and TSPAN9 as poor prognostic signature genes for CAF infiltration. The search on drug databases revealed collagenase clostridium histolyticum, ocriplasmin, halofuginone, natalizumab, firategrast, and BIO-1211 as the potential drugs for further investigation.

CONCLUSIONS

Our study demonstrated the central role of extracellular matrix components secreted and remodeled by CAFs in gastric cancer. The gene signature we identified in this study carries high potential as a predictive tool for poor prognosis in gastric cancer patients. Elucidating the mechanisms by which the signature genes contribute to poor patient outcomes can lead to the discovery of more potent molecular-targeted agents and increase the therapeutic success in gastric cancer.

Identification and Validation of Candidate Gene Module Along With Immune Cells Infiltration Patterns in Atherosclerosis Progression to Plaque Rupture via Transcriptome Analysis

Objective

To explore the differentially expressed genes (DEGs) along with infiltrating immune cells landscape and their potential mechanisms in the progression of atherosclerosis from onset to plaque rupture.

Methods

In this study, three atherosclerosis-related microarray datasets were downloaded from the NCBI-GEO database. The gene set enrichment analysis (GSEA) was performed for interpreting the biological insights of gene expression data. The CIBERSORTx algorithm was applied to infer the relative proportions of infiltrating immune cells of the atherosclerotic samples. DEGs of the datasets were screened using R. The protein interaction network was constructed via STRING. The cluster genes were analyzed by the Cytoscape software. Gene ontology (GO) enrichment was performed via geneontology.org. The least absolute shrinkage and selection operator (LASSO) logistic regression algorithm and receiver operating characteristics (ROC) analyses were performed to build machine learning models for differentiating atherosclerosis status. The Pearson correlation analysis was carried out to illustrate the relationship between cluster genes and immune cells. The expression levels of the cluster genes were validated in two external cohorts. Transcriptional factors and drug-gene interaction analysis were performed to investigate the promising targets for atherosclerosis intervention.

Results

Pathways related to immunoinflammatory responses were identified according to GSEA analysis, and the detailed fractions infiltrating immune cells were compared between the early and advanced atherosclerosis. Additionally, we identified 170 DEGs in atherosclerosis progression (|log2FC|≥1 and adjusted p < 0.05). They were mainly enriched in GO terms relating to inflammatory response and innate immune response. A cluster of nine genes, such as ITGB2, C1QC, LY86, CTSS, C1QA, CSF1R, LAPTM5, VSIG4, and CD163, were found to be significant, and their correlations with infiltrating immune cells were calculated. The cluster genes were also validated to be upregulated in two external cohorts. Moreover, C1QA and ITGB2 may exert pathogenic functions in the entire process of atherogenesis.

Conclusions

We reanalyzed the transcriptomic signature of atherosclerosis development from onset to plaque rupture along with the landscape of the immune cell, as well as revealed new insights and specific prospective DEGs for the investigation of disease-associated dynamic molecular processes and their regulations with immune cells.

Regulating DNA methylation could reduce neuronal ischemia response and apoptosis after ischemia-reperfusion injury

BACKGROUND

Ischemia-reperfusion injury (IRI) is an important pathophysiological condition that can cause cell injury and large-scale tissue injury in the nervous system. Previous studies have shown that epigenetic regulation may play a role in the pathogenesis of IRI.

METHODS

In this study, we isolated mouse cortical neurons and constructed an oxygen-glucose deprivation/reoxygenation (OGD) model to explore the change in DNA methylation and its effect on the expression of corresponding genes.

RESULTS

We found that DNA methylation in neurons increased with hypoxia duration and that hypermethylation of numerous promoters and 3'-untranslated regions increased. We performed Gene Ontology enrichment analysis to study gene function and Kyoto Encyclopedia of Genes and Genomes pathway analysis to identify the pathways associated with gene regulation. The results showed that hypermethylation-related genes expressed after OGD were related to physiological pathways such as neuronal projection, ion transport, growth and development, while hypomethylation-related genes were related to pathological pathways such as the external apoptosis signaling pathway, neuronal death regulation, and regulation of oxidative stress. However, the changes in DNA methylation were specific for certain genes and may have been related to OGD-induced neuronal damage. Importantly, we integrated transcription and DNA methylation data to identify several candidate target genes, including hypomethylated Apoe, Pax6, Bmp4, and Ptch1 and hypermethylated Adora2a, Crhr1, Stxbp1, and Tac1. This study further indicated the effect of DNA methylation on gene function in brain IRI from the perspective of epigenetics, and the identified genes may become new targets for achieving neuroprotection in the brain after IRI.

Regulatory T Cell-Related Gene Indicators in Pulmonary Hypertension

Objective: Regulatory T cells (Tregs) are critical immune modulators to maintain immune homeostasis and limit pulmonary hypertension (PH). This study was aimed to identify Treg-related genes (TRGs) in PH.

Methods: The gene expression profile from lungs of PH patients was retrieved from the Gene Expression Omnibus (GEO) database. The abundance of Tregs was estimated by the xCell algorithm, the correlation of which with differentially expressed genes (DEGs) was performed. DEGs with a |Pearson correlation coefficient| >0.4 were identified as TRGs. Functional annotation and the protein–protein interaction (PPI) network were analyzed. A gene signature for 25 hub TRGs (TRGscore) was generated by a single sample scoring method to determine its accuracy to distinguish PH from control subjects. TRGs were validated in datasets of transcriptional profiling of PH cohorts and in lung tissues of experimental PH mice.

Results: A total of 819 DEGs were identified in lungs of 58 PAH patients compared to that of 25 control subjects of dataset GSE117261. In total, 165 of all these DEGs were correlated with the abundance of Tregs and identified as TRGs, with 90 upregulated genes and 75 downregulated genes compared to that of control subjects. The upregulated TRGs were enriched in negative regulation of multiple pathways, such as cAMP-mediated signaling and I-kappaB kinase/NF-kappaB signaling, and regulated by multiple genes encoding transcriptional factors including HIF1A. Furthermore, 25 hub genes categorized into three clusters out of 165 TRGs were derived, and we identified 27 potential drugs targeting 10 hub TRGs. The TRGscore based on 25 hub TRGs was higher in PH patients and could distinguish PH from control subjects (all AUC >0.7). Among them, 10 genes including NCF2, MNDA/Ifi211, HCK, FGR, CSF3R, AQP9, S100A8, G6PD/G6pdx, PGD, and TXNRD1 were significantly reduced in lungs of severe PH patients of dataset GSE24988 as well as in lungs of hypoxic PH mice compared to corresponding controls.

Conclusion: Our finding will shed some light on the Treg-associated therapeutic targets in the progression of PH and emphasize on TRGscore as a novel indicator for PH.

Genetic association and single-cell transcriptome analyses reveal distinct features connecting autoimmunity with cancers

Autoimmune diseases (ADs) are at a significantly higher risk of cancers with unclear mechanism. By searching GWAS catalog database and Medline, susceptible genes for five common ADs, including systemic lupus erythematosus (SLE), rheumatoid arthritis, Sjögren syndrome, systemic sclerosis, and idiopathic inflammatory myopathies, were collected and then were overlapped with cancer driver genes. Single-cell transcriptome analysis was performed in the comparation between SLE and related cancer. We identified 45 carcinogenic autoimmune disease risk (CAD) genes, which were mainly enriched in T cell signaling pathway and B cell signaling pathway. Integrated single-cell analysis revealed immune cell signaling was significantly downregulated in renal cancer compared with SLE, while stemness signature was significantly enriched in both renal cancer or lymphoma and SLE in specific subpopulations. Drugs targeting CAD genes were shared between ADs and cancer. Our study highlights the common and specific features between ADs and related cancers, and sheds light on a new discovery of treatments.

Integrative multi-omic analysis identifies genetically influenced DNA methylation biomarkers for breast and prostate cancers

Aberrant DNA methylation has emerged as a hallmark in several cancers and contributes to risk, oncogenesis, progression, and prognosis. In this study, we performed imputation-based and conventional methylome-wide association analyses for breast cancer (BrCa) and prostate cancer (PrCa). The imputation-based approach identified DNA methylation at cytosine-phosphate-guanine sites (CpGs) associated with BrCa and PrCa risk utilising genome-wide association summary statistics (N_BrCa = 228,951, N_PrCa = 140,254) and prebuilt methylation prediction models, while the conventional approach identified CpG associations utilising TCGA and GEO experimental methylation data (N_BrCa = 621, N_PrCa = 241). Enrichment analysis of the association results implicated 77 and 81 genetically influenced CpGs for BrCa and PrCa, respectively. Furthermore, analysis of differential gene expression around these CpGs suggests a genome-epigenome-transcriptome mechanistic relationship. Conditional analyses identified multiple independent secondary SNP associations (P_cond < 0.05) around 28 BrCa and 22 PrCa CpGs. Cross-cancer analysis identified eight common CpGs, including a strong therapeutic target in SREBF1 (17p11.2)—a key player in lipid metabolism. These findings highlight the utility of integrative analysis of multi-omic cancer data to identify robust biomarkers and understand their regulatory effects on cancer risk.

Methylome-wide association studies identify genetically-influenced CpGs associated with breast and prostate cancer risk and (epi)genome-transcriptome mechanistic relationships, with lipid metabolism genes implicated as potential therapeutic targets.

An Integrated Mass Spectrometry-Based Glycomics-Driven Glycoproteomics Analytical Platform to Functionally Characterize Glycosylation Inhibitors

Cancer progression is linked to aberrant protein glycosylation due to the overexpression of several glycosylation enzymes. These enzymes are underexploited as potential anticancer drug targets and the development of rapid-screening methods and identification of glycosylation inhibitors are highly sought. An integrated bioinformatics and mass spectrometry-based glycomics-driven glycoproteomics analysis pipeline was performed to identify an N-glycan inhibitor against lung cancer cells. Combined network pharmacology and in silico screening approaches were used to identify a potential inhibitor, pictilisib, against several glycosylation-related proteins, such as Alpha1-6FucT, GlcNAcT-V, and Alpha2,6-ST-I. A glycomics assay of lung cancer cells treated with pictilisib showed a significant reduction in the fucosylation and sialylation of N-glycans, with an increase in high mannose-type glycans. Proteomics analysis and in vitro assays also showed significant upregulation of the proteins involved in apoptosis and cell adhesion, and the downregulation of proteins involved in cell cycle regulation, mRNA processing, and protein translation. Site-specific glycoproteomics analysis further showed that glycoproteins with reduced fucosylation and sialylation were involved in apoptosis, cell adhesion, DNA damage repair, and chemical response processes. To determine how the alterations in N-glycosylation impact glycoprotein dynamics, modeling of changes in glycan interactions of the ITGA5-ITGB1 (Integrin alpha 5-Integrin beta-1) complex revealed specific glycosites at the interface of these proteins that, when highly fucosylated and sialylated, such as in untreated A549 cells, form greater hydrogen bonding interactions compared to the high mannose-types in pictilisib-treated A549 cells. This study highlights the use of mass spectrometry to identify a potential glycosylation inhibitor and assessment of its impact on cell surface glycoprotein abundance and protein-protein interaction.

Drug repositioning in non-small cell lung cancer (NSCLC) using gene co-expression and drug–gene interaction networks analysis

Lung cancer is the most common cancer in men and women. This cancer is divided into two main types, namely non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC). Around 85 to 90 percent of lung cancers are NSCLC. Repositioning potent candidate drugs in NSCLC treatment is one of the important topics in cancer studies. Drug repositioning (DR) or drug repurposing is a method for identifying new therapeutic uses of existing drugs. The current study applies a computational drug repositioning method to identify candidate drugs to treat NSCLC patients. To this end, at first, the transcriptomics profile of NSCLC and healthy (control) samples was obtained from the GEO database with the accession number GSE21933. Then, the gene co-expression network was reconstructed for NSCLC samples using the WGCNA, and two significant purple and magenta gene modules were extracted. Next, a list of transcription factor genes that regulate purple and magenta modules' genes was extracted from the TRRUST V2.0 online database, and the TF–TG (transcription factors–target genes) network was drawn. Afterward, a list of drugs targeting TF–TG genes was obtained from the DGIdb V4.0 database, and two drug–gene interaction networks, including drug-TG and drug-TF, were drawn. After analyzing gene co-expression TF–TG, and drug–gene interaction networks, 16 drugs were selected as potent candidates for NSCLC treatment. Out of 16 selected drugs, nine drugs, namely Methotrexate, Olanzapine, Haloperidol, Fluorouracil, Nifedipine, Paclitaxel, Verapamil, Dexamethasone, and Docetaxel, were chosen from the drug-TG sub-network. In addition, nine drugs, including Cisplatin, Daunorubicin, Dexamethasone, Methotrexate, Hydrocortisone, Doxorubicin, Azacitidine, Vorinostat, and Doxorubicin Hydrochloride, were selected from the drug-TF sub-network. Methotrexate and Dexamethasone are common in drug-TG and drug-TF sub-networks. In conclusion, this study proposed 16 drugs as potent candidates for NSCLC treatment through analyzing gene co-expression, TF–TG, and drug–gene interaction networks.

Novel Diagnostic and Therapeutic Options for KMT2A-Rearranged Acute Leukemias

The KMT2A (MLL) gene rearrangements (KMT2A-r) are associated with a diverse spectrum of acute leukemias. Although most KMT2A-r are restricted to nine partner genes, we have recently revealed that KMT2A-USP2 fusions are often missed during FISH screening of these genetic alterations. Therefore, complementary methods are important for appropriate detection of any KMT2A-r. Here we use a machine learning model to unravel the most appropriate markers for prediction of KMT2A-r in various types of acute leukemia. A Random Forest and LightGBM classifier was trained to predict KMT2A-r in patients with acute leukemia. Our results revealed a set of 20 genes capable of accurately estimating KMT2A-r. The SKIDA1 (AUC: 0.839; CI: 0.799-0.879) and LAMP5 (AUC: 0.746; CI: 0.685-0.806) overexpression were the better markers associated with KMT2A-r compared to CSPG4 (also named NG2; AUC: 0.722; CI: 0.659-0.784), regardless of the type of acute leukemia. Of importance, high expression levels of LAMP5 estimated the occurrence of all KMT2A-USP2 fusions. Also, we performed drug sensitivity analysis using IC50 data from 345 drugs available in the GDSC database to identify which ones could be used to treat KMT2A-r leukemia. We observed that KMT2A-r cell lines were more sensitive to 5-Fluorouracil (5FU), Gemcitabine (both antimetabolite chemotherapy drugs), WHI-P97 (JAK-3 inhibitor), Foretinib (MET/VEGFR inhibitor), SNX-2112 (Hsp90 inhibitor), AZD6482 (PI3Kβ inhibitor), KU-60019 (ATM kinase inhibitor), and Pevonedistat (NEDD8-activating enzyme (NAE) inhibitor). Moreover, IC50 data from analyses of ex-vivo drug sensitivity to small-molecule inhibitors reveals that Foretinib is a promising drug option for AML patients carrying FLT3 activating mutations. Thus, we provide novel and accurate options for the diagnostic screening and therapy of KMT2A-r leukemia, regardless of leukemia subtype.

Mutation Burden Analysis of Six Common Mental Disorders in African Americans by Whole Genome Sequencing

Mental disorders present a global health concern, with limited treatment options. In today's medical practice, medications such as antidepressants are prescribed not only for depression, but also for conditions such as anxiety and attention deficit hyperactivity disorder (ADHD). Therefore, identifying gene targets for specific disorders is important and offers improved precision. In this study, we performed a genetic analysis of six common mental disorders, ADHD, anxiety, depression, delays in mental developments, intellectual disabilities (ID), and speech/language disorder in the ethnic minority of African Americans (AA) using whole genome sequencing (WGS). WGS data was generated from blood-derived DNA from 4178 AA individuals, including 1384 patients with the diagnosis of at least one mental disorder. Mutation burden analysis was applied based on rare and deleterious mutations in the AA population between cases and controls, and further analyzed in the context of patients with single mental disorder diagnosis. Certain genes uncovered demonstrated significant p values in mutation burden analysis. In addition, exclusive recurrences in specific type of disorder were scanned through gene-drug interaction databases to assess for availability of potential medications. We uncovered 15 genes harboring deleterious mutations, including HMGCR and UST for ADHD; FNTB for anxiety, XIRP2, NPPC, , STK33, PANX1 and NTS for depression; RUNX3, TACR1, and NDUFS7 for delays in mental developments; HPN for ID; COL6A3, DDB1, and NDUFA11 for speech/language disorder. Taken together, we have established critical insight into the development of new precision medicine approaches for mental disorders in African Americans.

Analysis of multiple databases identifies crucial genes correlated with prognosis of hepatocellular carcinoma

Despite advancements made in the therapeutic strategies on hepatocellular carcinoma (HCC), the survival rate of HCC patient is not satisfactory enough. Therefore, there is an urgent need for the valuable prognostic biomarkers in HCC therapy. In this study, we aimed to screen hub genes correlated with prognosis of HCC via multiple databases. 117 HCC-related genes were obtained from the intersection of the four databases. We subsequently identify 10 hub genes (JUN, IL10, CD34, MTOR, PTGS2, PTPRC, SELE, CSF1, APOB, MUC1) from PPI network by Cytoscape software analysis. Significant differential expression of hub genes between HCC tissues and adjacent tissues were observed in UALCAN, HCCDB and HPA databases. These hub genes were significantly associated with immune cell infiltrations and immune checkpoints. The hub genes were correlated with clinical parameters and survival probability of HCC patients. 147 potential targeted therapeutic drugs for HCC were identified through the DGIdb database. These hub genes could be used as novel prognostic biomarkers for HCC therapy.

Correlation of Immunological and Molecular Profiles with Response to Crizotinib in Alveolar Soft Part Sarcoma: An Exploratory Study Related to the EORTC 90101 "CREATE" Trial

Alveolar soft part sarcoma (ASPS) is a rare subtype of soft tissue sarcoma characterized by an unbalanced translocation, resulting in ASPSCR1-TFE3 fusion that transcriptionally upregulates MET expression. The European Organization for Research and Treatment of Cancer (EORTC) 90101 “CREATE” phase II trial evaluated the MET inhibitor crizotinib in ASPS patients, achieving only limited antitumor activity. We performed a comprehensive molecular analysis of ASPS tissue samples collected in this trial to identify potential biomarkers correlating with treatment outcome. A tissue microarray containing 47 ASPS cases was used for the characterization of the tumor microenvironment using multiplex immunofluorescence. DNA isolated from 34 available tumor samples was analyzed to detect recurrent gene copy number alterations (CNAs) and mutations by low-coverage whole-genome sequencing and whole-exome sequencing. Pathway enrichment analysis was used to identify diseased-associated pathways in ASPS sarcomagenesis. Kaplan–Meier estimates, Cox regression, and the Fisher’s exact test were used to correlate histopathological and molecular findings with clinical data related to crizotinib treatment, aiming to identify potential factors associated with patient outcome. Tumor microenvironment characterization showed the presence of PD-L1 and CTLA-4 in 10 and 2 tumors, respectively, and the absence of PD-1 in all specimens. Apart from CD68, other immunological markers were rarely expressed, suggesting a low level of tumor-infiltrating lymphocytes in ASPS. By CNA analysis, we detected a number of broad and focal alterations. The most common alteration was the loss of chromosomal region 1p36.32 in 44% of cases. The loss of chromosomal regions 1p36.32, 1p33, 1p22.2, and 8p was associated with shorter progression-free survival. Using whole-exome sequencing, 13 cancer-associated genes were found to be mutated in at least three cases. Pathway enrichment analysis identified genetic alterations in NOTCH signaling, chromatin organization, and SUMOylation pathways. NOTCH4 intracellular domain dysregulation was associated with poor outcome, while inactivation of the beta-catenin/TCF complex correlated with improved outcome in patients receiving crizotinib. ASPS is characterized by molecular heterogeneity. We identify genetic aberrations potentially predictive of treatment outcome during crizotinib therapy and provide additional insights into the biology of ASPS, paving the way to improve treatment approaches for this extremely rare malignancy.

Innovative in Silico Approaches for Characterization of Genes and Proteins

Bioinformatics is an amalgamation of biology, mathematics and computer science. It is a science which gathers the information from biology in terms of molecules and applies the informatic techniques to the gathered information for understanding and organizing the data in a useful manner. With the help of bioinformatics, the experimental data generated is stored in several databases available online like nucleotide database, protein databases, GENBANK and others. The data stored in these databases is used as reference for experimental evaluation and validation. Till now several online tools have been developed to analyze the genomic, transcriptomic, proteomics, epigenomics and metabolomics data. Some of them include Human Splicing Finder (HSF), Exonic Splicing Enhancer Mutation taster, and others. A number of SNPs are observed in the non-coding, intronic regions and play a role in the regulation of genes, which may or may not directly impose an effect on the protein expression. Many mutations are thought to influence the splicing mechanism by affecting the existing splice sites or creating a new sites. To predict the effect of mutation (SNP) on splicing mechanism/signal, HSF was developed. Thus, the tool is helpful in predicting the effect of mutations on splicing signals and can provide data even for better understanding of the intronic mutations that can be further validated experimentally. Additionally, rapid advancement in proteomics have steered researchers to organize the study of protein structure, function, relationships, and dynamics in space and time. Thus the effective integration of all of these technological interventions will eventually lead to steering up of next-generation systems biology, which will provide valuable biological insights in the field of research, diagnostic, therapeutic and development of personalized medicine.

Transcriptional Profiling of Hippocampus Identifies Network Alterations in Alzheimer’s Disease

Alzheimer’s disease (AD) is a neurodegenerative disease characterized by rapid brain cell degeneration affecting different areas of the brain. Hippocampus is one of the earliest involved brain regions in the disease. Modern technologies based on high-throughput data have identified transcriptional profiling of several neurological diseases, including AD, for a better comprehension of genetic mechanisms of the disease. In this study, we investigated differentially expressed genes (DEGs) from six Gene Expression Omnibus (GEO) datasets of hippocampus of AD patients. The identified DEGs were submitted to Weighted correlation network analysis (WGCNA) and ClueGo to explore genes with a higher degree centrality and to comprehend their biological role. Subsequently, MCODE was used to identify subnetworks of interconnected DEGs. Our study found 40 down-regulated genes and 36 up-regulated genes as consensus DEGs. Analysis of the co-expression network revealed ACOT7, ATP8A2, CDC42, GAD1, GOT1, INA, NCALD, and WWTR1 to be genes with a higher degree centrality. ClueGO revealed the pathways that were mainly enriched, such as clathrin coat assembly, synaptic vesicle endocytosis, and DNA damage response signal transduction by p53 class mediator. In addition, we found a subnetwork of 12 interconnected genes (AMPH, CA10, CALY, NEFL, SNAP25, SNAP91, SNCB, STMN2, SV2B, SYN2, SYT1, and SYT13). Only CA10 and CALY are targets of known drugs while the others could be potential novel drug targets.

Central and Peripheral Immune Dysregulation in Posttraumatic Stress Disorder: Convergent Multi-Omics Evidence

Posttraumatic stress disorder (PTSD) is a chronic and multifactorial disorder with a prevalence ranging between 6-10% in the general population and ~35% in individuals with high lifetime trauma exposure. Growing evidence indicates that the immune system may contribute to the etiology of PTSD, suggesting the inflammatory dysregulation as a hallmark feature of PTSD. However, the potential interplay between the central and peripheral immune system, as well as the biological mechanisms underlying this dysregulation remain poorly understood. The activation of the HPA axis after trauma exposure and the subsequent activation of the inflammatory system mediated by glucocorticoids is the most common mechanism that orchestrates an exacerbated immunological response in PTSD. Recent high-throughput analyses in peripheral and brain tissue from both humans with and animal models of PTSD have found that changes in gene regulation via epigenetic alterations may participate in the impaired inflammatory signaling in PTSD. The goal of this review is to assess the role of the inflammatory system in PTSD across tissue and species, with a particular focus on the genomics, transcriptomics, epigenomics, and proteomics domains. We conducted an integrative multi-omics approach identifying TNF (Tumor Necrosis Factor) signaling, interleukins, chemokines, Toll-like receptors and glucocorticoids among the common dysregulated pathways in both central and peripheral immune systems in PTSD and propose potential novel drug targets for PTSD treatment.

Multivariate genome-wide association study on tissue-sensitive diffusion metrics highlights pathways that shape the human brain

The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discover 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% are validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development are found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders.

Integrated computational approaches to screen gene expression data to determine key genes and therapeutic targets for type-2 diabetes mellitus

There is a rapid rise in cases of Type-2-diabetes mellitus (T2DM) globally, irrespective of the geography, ethnicity or any other variable factors. The molecular mechanisms that could cause the condition of T2DM need to be more thoroughly analysed to understand the clinical manifestations and to derive better therapeutic regimes. Tools in bioinformatics are used to trace out key gene elements and to identify the key causative gene elements and their possible therapeutic agents. Microarray datasets were retrieved from the Gene expression omnibus database and studied using R to derive different expressed gene (DEG) elements. With the comparison of the expressed genes with disease specific genes in DisGeNET, the final annotated genes were taken for analysis. Gene Ontology studies, Protein-protein interaction (PPI), Co-expression analysis, Gene-drug interactions were performed to scale down the hub genes and to identify the novelty across the genes analysed so far. In vivo and invitro analysis of key genes and the trace of interaction pathway is crucial to better understand the unique outcomes from the novel genes, forming the basis to understand the pathway that ends up causing T2DM. Afterwards, docking was executed enabling recognition of interacting residues involved in inhibition. The complex CCL5-265 and CD8A-40585 thus docked showed best results as is evident from its PCA analysis and MMGBSA calculation. There is now scope for deriving candidate drugs that could possibly detect personalized therapies for T2DM.

In silico identification of novel biomarkers for key players in transition from normal colon tissue to adenomatous polyps

Adenomatous polyps of the colon are the most common neoplastic polyps. Although most of adenomatous polyps do not show malign transformation, majority of colorectal carcinomas originate from neoplastic polyps. Therefore, understanding of this transformation process would help in both preventive therapies and evaluation of malignancy risks. This study uncovers alterations in gene expressions as potential biomarkers that are revealed by integration of several network-based approaches. In silico analysis performed on a unified microarray cohort, which is covering 150 normal colon and adenomatous polyp samples. Significant gene modules were obtained by a weighted gene co-expression network analysis. Gene modules with similar profiles were mapped to a colon tissue specific functional interaction network. Several clustering algorithms run on the colon-specific network and the most significant sub-modules between the clusters were identified. The biomarkers were selected by filtering differentially expressed genes which also involve in significant biological processes and pathways. Biomarkers were also validated on two independent datasets based on their differential gene expressions. To the best of our knowledge, such a cascaded network analysis pipeline was implemented for the first time on a large collection of normal colon and polyp samples. We identified significant increases in TLR4 and MSX1 expressions as well as decrease in chemokine profiles with mostly pro-tumoral activities. These biomarkers might appear as both preventive targets and biomarkers for risk evaluation. As a result, this research proposes novel molecular markers that might be alternative to endoscopic approaches for diagnosis of adenomatous polyps.

Epithelial-mesenchymal transition-related genes in coronary artery disease

Epithelial-mesenchymal transition (EMT) is critical in the development of coronary artery disease (CAD). However, landscapes of EMT-related genes have not been fully established in CAD. We identified the differentially expressed mRNAs and lncRNAs (DElncRNAs) from the Gene Expression Omnibus database. Pearson’s correlation analysis, the least absolute shrinkage and selection operator regression, and support vector machine reverse feature elimination algorithms were used to screen EMT-related lncRNAs. The cis–trans regulatory networks were constructed based on EMT-related lncRNAs. Quantitative real-time polymerase chain reaction was performed to validate the expression of EMT-related genes in a cohort of six patients with CAD and six healthy controls. We further estimated the infiltration of the immune cells in CAD patients with five algorithms, and the correlation between EMT-related genes and infiltrating immune cells was analyzed. We identified eight EMT-related lncRNAs in CAD. The area under curve value was greater than 0.95. The immune analysis revealed significant CD8 T cells, monocytes, and NK cells in CAD and found that EMT-related lncRNAs were correlated with these immune cell subsets. Moreover, SNAI2, an EMT-TF gene, was found in the trans-regulatory network of EMT-related lncRNAs. Further, we found SNAI2 as a biomarker for the diagnosis of CAD but it also had a close correlation with immune cell subsets in CAD. Eight EMT-related lncRNAs and SNAI2 have important significance in the diagnosis of CAD patients.

From pharmacogenetics to pharmaco-omics: Milestones and future directions

The origins of pharmacogenetics date back to the 1950s, when it was established that inter-individual differences in drug response are partially determined by genetic factors. Since then, pharmacogenetics has grown into its own field, motivated by the translation of identified gene-drug interactions into therapeutic applications. Despite numerous challenges ahead, our understanding of the human pharmacogenetic landscape has greatly improved thanks to the integration of tools originating from disciplines as diverse as biochemistry, molecular biology, statistics, and computer sciences. In this review, we discuss past, present, and future developments of pharmacogenetics methodology, focusing on three milestones: how early research established the genetic basis of drug responses, how technological progress made it possible to assess the full extent of pharmacological variants, and how multi-dimensional omics datasets can improve the identification, functional validation, and mechanistic understanding of the interplay between genes and drugs. We outline novel strategies to repurpose and integrate molecular and clinical data originating from biobanks to gain insights analogous to those obtained from randomized controlled trials. Emphasizing the importance of increased diversity, we envision future directions for the field that should pave the way to the clinical implementation of pharmacogenetics.

A conditional gene-based association framework integrating isoform-level eQTL data reveals new susceptibility genes for schizophrenia

Linkage disequilibrium and disease-associated variants in the non-coding regions make it difficult to distinguish the truly associated genes from the redundantly associated genes for complex diseases. In this study, we proposed a new conditional gene-based framework called eDESE that leveraged an improved effective chi-squared statistic to control the type I error rates and remove the redundant associations. eDESE initially performed the association analysis by mapping variants to genes according to their physical distance. We further demonstrated that the isoform-level eQTLs could be more powerful than the gene-level eQTLs in the association analysis using a simulation study. Then the eQTL-guided strategies, that is, mapping variants to genes according to their gene/isoform-level variant-gene cis-eQTLs associations, were also integrated with eDESE. We then applied eDESE to predict the potential susceptibility genes of schizophrenia and found that the potential susceptibility genes were enriched with many neuronal or synaptic signaling-related terms in the Gene Ontology knowledgebase and antipsychotics-gene interaction terms in the drug-gene interaction database (DGIdb). More importantly, seven potential susceptibility genes identified by eDESE were the target genes of multiple antipsychotics in DrugBank. Comparing the potential susceptibility genes identified by eDESE and other benchmark approaches (i.e., MAGMA and S-PrediXcan) implied that strategy based on the isoform-level eQTLs could be an important supplement for the other two strategies (physical distance and gene-level eQTLs). We have implemented eDESE in our integrative platform KGGSEE (http://pmglab.top/kggsee/#/) and hope that eDESE can facilitate the prediction of candidate susceptibility genes and isoforms for complex diseases in a multi-tissue context.

The Shared Mechanism and Candidate Drugs of Multiple Sclerosis and Sjögren's Syndrome Analyzed by Bioinformatics Based on GWAS and Transcriptome Data

Objective

This study aimed to explore the shared mechanism and candidate drugs of multiple sclerosis (MS) and Sjögren's syndrome (SS).

Methods

MS- and SS-related susceptibility genes and differentially expressed genes (DEGs) were identified by bioinformatics analysis based on genome-wide association studies (GWAS) and transcriptome data from GWAS catalog and Gene Expression Omnibus (GEO) database. Pathway enrichment, Gene Ontology (GO) analysis, and protein-protein interaction analysis for susceptibility genes and DEGs were performed. The drugs targeting common pathways/genes were obtained through Comparative Toxicogenomics Database (CTD), DrugBank database, and Drug-Gene Interaction (DGI) Database. The target genes of approved/investigational drugs for MS and SS were obtained through DrugBank and compared with the common susceptibility genes.

Results

Based on GWAS data, we found 14 hub common susceptibility genes (HLA-DRB1, HLA-DRA, STAT3, JAK1, HLA-B, HLA-DQA1, HLA-DQA2, HLA-DQB1, HLA-DRB5, HLA-DPA1, HLA-DPB1, TYK2, IL2RA, and MAPK1), with 8 drugs targeting two or more than two genes, and 28 common susceptibility pathways, with 15 drugs targeting three or more than three pathways. Based on transcriptome data, we found 3 hub common DEGs (STAT1, GATA3, PIK3CA) with 3 drugs and 10 common risk pathways with 435 drugs. "JAK-STAT signaling pathway" was included in common susceptibility pathways and common risk pathways at the same time. There were 133 overlaps including JAK-STAT inhibitors between agents from GWAS and transcriptome data. Besides, we found that IL2RA and HLA-DRB1, identified as hub common susceptibility genes, were the targets of daclizumab and glatiramer that were used for MS, indicating that daclizumab and glatiramer may be therapeutic for SS.

Conclusion

We observed the shared mechanism of MS and SS, in which JAK-STAT signaling pathway played a vital role, which may be the genetic and molecular bases of comorbidity of MS with SS. Moreover, JAK-STAT inhibitors were potential therapies for MS and SS, especially for their comorbidity.