Pathway Commons 2019 Update: integration, analysis and exploration of pathway data

A hybrid approach unveils drug repurposing candidates targeting an Alzheimer pathophysiology mechanism

The high number of failed pre-clinical and clinical studies for compounds targeting Alzheimer disease (AD) has demonstrated that there is a need to reassess existing strategies. Here, we pursue a holistic, mechanism-centric drug repurposing approach combining computational analytics and experimental screening data. Based on this integrative workflow, we identified 77 druggable modifiers of tau phosphorylation (pTau). One of the upstream modulators of pTau, HDAC6, was screened with 5,632 drugs in a tau-specific assay, resulting in the identification of 20 repurposing candidates. Four compounds and their known targets were found to have a link to AD-specific genes. Our approach can be applied to a variety of AD-associated pathophysiological mechanisms to identify more repurposing candidates.

Application of Transcriptomics for Predicting Protein Interaction Networks, Drug Targets and Drug Candidates

Protein interaction pathways and networks are critically-required for a vast range of biological processes. Improved discovery of candidate druggable proteins within specific cell, tissue and disease contexts will aid development of new treatments. Predicting protein interaction networks from gene expression data can provide valuable insights into normal and disease biology. For example, the resulting protein networks can be used to identify potentially druggable targets and drug candidates for testing in cell and animal disease models. The advent of whole-transcriptome expression profiling techniques—that catalogue protein-coding genes expressed within cells and tissues—has enabled development of individual algorithms for particular tasks. For example,: (i) gene ontology algorithms that predict gene/protein subsets involved in related cell processes; (ii) algorithms that predict intracellular protein interaction pathways; and (iii) algorithms that correlate druggable protein targets with known drugs and/or drug candidates. This review examines approaches, advantages and disadvantages of existing gene expression, gene ontology, and protein network prediction algorithms. Using this framework, we examine current efforts to combine these algorithms into pipelines to enable identification of druggable targets, and associated known drugs, using gene expression datasets. In doing so, new opportunities are identified for development of powerful algorithm pipelines, suitable for wide use by non-bioinformaticians, that can predict protein interaction networks, druggable proteins, and related drugs from user gene expression datasets.

Novel epigenetic network biomarkers for early detection of esophageal cancer

BACKGROUND

Early detection of esophageal cancer is critical to improve survival. Whilst studies have identified biomarkers, their interpretation and validity is often confounded by cell-type heterogeneity.

RESULTS

Here we applied systems-epigenomic and cell-type deconvolution algorithms to a discovery set encompassing RNA-Seq and DNA methylation data from esophageal adenocarcinoma (EAC) patients and matched normal-adjacent tissue, in order to identify robust biomarkers, free from the confounding effect posed by cell-type heterogeneity. We identify 12 gene-modules that are epigenetically deregulated in EAC, and are able to validate all 12 modules in 4 independent EAC cohorts. We demonstrate that the epigenetic deregulation is present in the epithelial compartment of EAC-tissue. Using single-cell RNA-Seq data we show that one of these modules, a proto-cadherin module centered around CTNND2, is inactivated in Barrett's Esophagus, a precursor lesion to EAC. By measuring DNA methylation in saliva from EAC cases and controls, we identify a chemokine module centered around CCL20, whose methylation patterns in saliva correlate with EAC status.

CONCLUSIONS

Given our observations that a CCL20 chemokine network is overactivated in EAC tissue and saliva from EAC patients, and that in independent studies CCL20 has been found to be overactivated in EAC tissue infected with the bacterium F. nucleatum, a bacterium that normally inhabits the oral cavity, our results highlight the possibility of using DNAm measurements in saliva as a proxy for changes occurring in the esophageal epithelium. Both the CTNND2/CCL20 modules represent novel promising network biomarkers for EAC that merit further investigation.

Metadichol®: A Novel Nanolipid Formulation That Inhibits SARS-CoV-2 and a Multitude of Pathological Viruses In Vitro

Increasing outbreaks of new pathogenic viruses have promoted the exploration of novel alternatives to time-consuming vaccines. Thus, it is necessary to develop a universal approach to halt the spread of new and unknown viruses as they are discovered. One such promising approach is to target lipid membranes, which are common to all viruses and bacteria. The ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has reaffirmed the importance of interactions between the virus envelope and the host cell plasma membrane as a critical mechanism of infection. Metadichol®, a nanolipid emulsion of long-chain alcohols, has been demonstrated as a strong candidate that inhibits the proliferation of SARS-CoV-2. Naturally derived substances, such as long-chain saturated lipid alcohols, reduce viral infectivity, including that of coronaviruses (such as SARS-CoV-2) by modifying their lipid-dependent attachment mechanism to human host cells. The receptor ACE2 mediates the entry of SARS-CoV-2 into the host cells, whereas the serine protease TMPRSS2 primes the viral S protein. In this study, Metadichol® was found to be 270 times more potent an inhibitor of TMPRSS2 (EC50 = 96 ng/mL) than camostat mesylate (EC50 = 26000 ng/mL). Additionally, it inhibits ACE with an EC50 of 71 ng/mL, but it is a very weak inhibitor of ACE2 at an EC50 of 31 μg/mL. Furthermore, the live viral assay performed in Caco-2 cells revealed that Metadichol® inhibits SARS-CoV-2 replication at an EC90 of 0.16 μg/mL. Moreover, Metadichol® had an EC90 of 0.00037 μM, making it 2081 and 3371 times more potent than remdesivir (EC50 = 0.77 μM) and chloroquine (EC50 = 1.14 μM), respectively.

Analyzing causal relationships in proteomic profiles using CausalPath

CausalPath (causalpath.org) evaluates proteomic measurements against prior knowledge of biological pathways and infers causality between changes in measured features, such as global protein and phospho-protein levels. It uses pathway resources to determine potential causality between observable omic features, which are called prior relations. The subset of the prior relations that are supported by the proteomic profiles are reported and evaluated for statistical significance. The end result is a network model of signaling that explains the patterns observed in the experimental dataset. For complete details on the use and execution of this protocol, please refer to Babur et al. (2021).

XDeathDB: a visualization platform for cell death molecular interactions

Lots of cell death initiator and effector molecules, signalling pathways and subcellular sites have been identified as key mediators in both cell death processes in cancer. The XDeathDB visualization platform provides a comprehensive cell death and their crosstalk resource for deciphering the signaling network organization of interactions among different cell death modes associated with 1461 cancer types and COVID-19, with an aim to understand the molecular mechanisms of physiological cell death in disease and facilitate systems-oriented novel drug discovery in inducing cell deaths properly. Apoptosis, autosis, efferocytosis, ferroptosis, immunogenic cell death, intrinsic apoptosis, lysosomal cell death, mitotic cell death, mitochondrial permeability transition, necroptosis, parthanatos, and pyroptosis related to 12 cell deaths and their crosstalk can be observed systematically by the platform. Big data for cell death gene-disease associations, gene-cell death pathway associations, pathway-cell death mode associations, and cell death-cell death associations is collected by literature review articles and public database from iRefIndex, STRING, BioGRID, Reactom, Pathway's commons, DisGeNET, DrugBank, and Therapeutic Target Database (TTD). An interactive webtool, XDeathDB, is built by web applications with R-Shiny, JavaScript (JS) and Shiny Server Iso. With this platform, users can search specific interactions from vast interdependent networks that occur in the realm of cell death. A multilayer spectral graph clustering method that performs convex layer aggregation to identify crosstalk function among cell death modes for a specific cancer. 147 hallmark genes of cell death could be observed in detail in these networks. These potential druggable targets are displayed systematically and tailoring networks to visualize specified relations is available to fulfil user-specific needs. Users can access XDeathDB for free at https://pcm2019.shinyapps.io/XDeathDB/ .

Author-sourced capture of pathway knowledge in computable form using Biofactoid

Making the knowledge contained in scientific papers machine-readable and formally computable would allow researchers to take full advantage of this information by enabling integration with other knowledge sources to support data analysis and interpretation. Here we describe Biofactoid, a web-based platform that allows scientists to specify networks of interactions between genes, their products, and chemical compounds, and then translates this information into a representation suitable for computational analysis, search and discovery. We also report the results of a pilot study to encourage the wide adoption of Biofactoid by the scientific community.

Performance Assessment of the Network Reconstruction Approaches on Various Interactomes

Beyond the list of molecules, there is a necessity to collectively consider multiple sets of omic data and to reconstruct the connections between the molecules. Especially, pathway reconstruction is crucial to understanding disease biology because abnormal cellular signaling may be pathological. The main challenge is how to integrate the data together in an accurate way. In this study, we aim to comparatively analyze the performance of a set of network reconstruction algorithms on multiple reference interactomes. We first explored several human protein interactomes, including PathwayCommons, OmniPath, HIPPIE, iRefWeb, STRING, and ConsensusPathDB. The comparison is based on the coverage of each interactome in terms of cancer driver proteins, structural information of protein interactions, and the bias toward well-studied proteins. We next used these interactomes to evaluate the performance of network reconstruction algorithms including all-pair shortest path, heat diffusion with flux, personalized PageRank with flux, and prize-collecting Steiner forest (PCSF) approaches. Each approach has its own merits and weaknesses. Among them, PCSF had the most balanced performance in terms of precision and recall scores when 28 pathways from NetPath were reconstructed using the listed algorithms. Additionally, the reference interactome affects the performance of the network reconstruction approaches. The coverage and disease- or tissue-specificity of each interactome may vary, which may result in differences in the reconstructed networks.

COVID19 Disease Map, a computational knowledge repository of virus-host interaction mechanisms

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.

Low Levels of TRIM28-Interacting KRAB-ZNF Genes Associate with Cancer Stemness and Predict Poor Prognosis of Kidney Renal Clear Cell Carcinoma Patients

Simple Summary

This is the first report investigating the involvement of TRIM28-interacting KRAB-ZNFs in kidney cancer progression. We demonstrate a significant negative association between KRAB-ZNFs and cancer stemness followed by an attenuated immune-suppressive response and reveal the prognostic role for several KRAB-ZNFs. Our findings may help better understand the molecular basis of kidney cancer and ultimately pave the way to more appropriate prognostic tools and novel therapeutic strategies directly eradicating the dedifferentiated compartment of the tumor.

Abstract

Krüppel-associated box zinc finger (KRAB-ZNF) proteins are known to regulate diverse biological processes, such as embryonic development, tissue-specific gene expression, and cancer progression. However, their involvement in the regulation of cancer stemness-like phenotype acquisition and maintenance is scarcely explored across solid tumor types, and to date, there are no data for kidney renal clear cell cancer (KIRC). We have harnessed The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database transcriptomic data and used several bioinformatic tools (i.e., GEPIA2, GSCALite, TISIDB, GSEA, CIBERSORT) to verify the relation between the expression and genomic alterations in KRAB-ZNFs and kidney cancer, focusing primarily on tumor dedifferentiation status and antitumor immune response. Our results demonstrate a significant negative correlation between KRAB-ZNFs and kidney cancer dedifferentiation status followed by an attenuated immune-suppressive response. The transcriptomic profiles of high KRAB-ZNF-expressing kidney tumors are significantly enriched with stem cell markers and show a depletion of several inflammatory pathways known for favoring cancer stemness. Moreover, we show for the first time the prognostic role for several KRAB-ZNFs in kidney cancer. Our results provide new insight into the role of selected KRAB-ZNF proteins in kidney cancer development. We believe that our findings may help better understand the molecular basis of KIRC.

Drug Discovery of Spinal Muscular Atrophy (SMA) from the Computational Perspective: A Comprehensive Review

Spinal muscular atrophy (SMA), one of the leading inherited causes of child mortality, is a rare neuromuscular disease arising from loss-of-function mutations of the survival motor neuron 1 (SMN1) gene, which encodes the SMN protein. When lacking the SMN protein in neurons, patients suffer from muscle weakness and atrophy, and in the severe cases, respiratory failure and death. Several therapeutic approaches show promise with human testing and three medications have been approved by the U.S. Food and Drug Administration (FDA) to date. Despite the shown promise of these approved therapies, there are some crucial limitations, one of the most important being the cost. The FDA-approved drugs are high-priced and are shortlisted among the most expensive treatments in the world. The price is still far beyond affordable and may serve as a burden for patients. The blooming of the biomedical data and advancement of computational approaches have opened new possibilities for SMA therapeutic development. This article highlights the present status of computationally aided approaches, including in silico drug repurposing, network driven drug discovery as well as artificial intelligence (AI)-assisted drug discovery, and discusses the future prospects.

Biomolecule and Bioentity Interaction Databases in Systems Biology: A Comprehensive Review

Technological advances in high-throughput techniques have resulted in tremendous growth of complex biological datasets providing evidence regarding various biomolecular interactions. To cope with this data flood, computational approaches, web services, and databases have been implemented to deal with issues such as data integration, visualization, exploration, organization, scalability, and complexity. Nevertheless, as the number of such sets increases, it is becoming more and more difficult for an end user to know what the scope and focus of each repository is and how redundant the information between them is. Several repositories have a more general scope, while others focus on specialized aspects, such as specific organisms or biological systems. Unfortunately, many of these databases are self-contained or poorly documented and maintained. For a clearer view, in this article we provide a comprehensive categorization, comparison and evaluation of such repositories for different bioentity interaction types. We discuss most of the publicly available services based on their content, sources of information, data representation methods, user-friendliness, scope and interconnectivity, and we comment on their strengths and weaknesses. We aim for this review to reach a broad readership varying from biomedical beginners to experts and serve as a reference article in the field of Network Biology.

Comprior: facilitating the implementation and automated benchmarking of prior knowledge-based feature selection approaches on gene expression data sets

BACKGROUND

Reproducible benchmarking is important for assessing the effectiveness of novel feature selection approaches applied on gene expression data, especially for prior knowledge approaches that incorporate biological information from online knowledge bases. However, no full-fledged benchmarking system exists that is extensible, provides built-in feature selection approaches, and a comprehensive result assessment encompassing classification performance, robustness, and biological relevance. Moreover, the particular needs of prior knowledge feature selection approaches, i.e. uniform access to knowledge bases, are not addressed. As a consequence, prior knowledge approaches are not evaluated amongst each other, leaving open questions regarding their effectiveness.

RESULTS

We present the Comprior benchmark tool, which facilitates the rapid development and effortless benchmarking of feature selection approaches, with a special focus on prior knowledge approaches. Comprior is extensible by custom approaches, offers built-in standard feature selection approaches, enables uniform access to multiple knowledge bases, and provides a customizable evaluation infrastructure to compare multiple feature selection approaches regarding their classification performance, robustness, runtime, and biological relevance.

CONCLUSION

Comprior allows reproducible benchmarking especially of prior knowledge approaches, which facilitates their applicability and for the first time enables a comprehensive assessment of their effectiveness.

The status of causality in biological databases: data resources and data retrieval possibilities to support logical modeling

Causal molecular interactions represent key building blocks used in computational modeling, where they facilitate the assembly of regulatory networks. Logical regulatory networks can be used to predict biological and cellular behaviors by system perturbations and in silico simulations. Today, broad sets of causal interactions are available in a variety of biological knowledge resources. However, different visions, based on distinct biological interests, have led to the development of multiple ways to describe and annotate causal molecular interactions. It can therefore be challenging to efficiently explore various resources of causal interaction and maintain an overview of recorded contextual information that ensures valid use of the data. This review lists the different types of public resources with causal interactions, the different views on biological processes that they represent, the various data formats they use for data representation and storage, and the data exchange and conversion procedures that are available to extract and download these interactions. This may further raise awareness among the targeted audience, i.e. logical modelers and other scientists interested in molecular causal interactions, but also database managers and curators, about the abundance and variety of causal molecular interaction data, and the variety of tools and approaches to convert them into one interoperable resource.

Proteo3Dnet: a web server for the integration of structural information with interactomics data

Proteo3Dnet is a web server dedicated to the analysis of mass spectrometry interactomics experiments. Given a flat list of proteins, its aim is to organize it in terms of structural interactions to provide a clearer overview of the data. This is achieved using three means: (i) the search for interologs with resolved structure available in the protein data bank, including cross-species remote homology search, (ii) the search for possibly weaker interactions mediated through Short Linear Motifs as predicted by ELM-a unique feature of Proteo3Dnet, (iii) the search for protein-protein interactions physically validated in the BioGRID database. The server then compiles this information and returns a graph of the identified interactions and details about the different searches. The graph can be interactively explored to understand the way the core complexes identified could interact. It can also suggest undetected partners to the experimentalists, or specific cases of conditionally exclusive binding. The interest of Proteo3Dnet, previously demonstrated for the difficult cases of the proteasome and pragmin complexes data is, here, illustrated in the context of yeast precursors to the small ribosomal subunits and the smaller interactome of 14-3-3zeta frequent interactors. The Proteo3Dnet web server is accessible at http://bioserv.rpbs.univ-paris-diderot.fr/services/Proteo3Dnet/.

Ultra-fast scalable estimation of single-cell differentiation potency from scRNA-Seq data

Motivation

An important task in the analysis of single-cell RNA-Seq data is the estimation of differentiation potency, as this can help identify stem-or-multipotent cells in non-temporal studies or in tissues where differentiation hierarchies are not well established. A key challenge in the estimation of single-cell potency is the need for a fast and accurate algorithm, scalable to large scRNA-Seq studies profiling millions of cells.

Results

Here, we present a single-cell potency measure, called Correlation of Connectome and Transcriptome (CCAT), which can return accurate single-cell potency estimates of a million cells in minutes, a 100-fold improvement over current state-of-the-art methods. We benchmark CCAT against 8 other single-cell potency models and across 28 scRNA-Seq studies, encompassing over 2 million cells, demonstrating comparable accuracy than the current state-of-the-art, at a significantly reduced computational cost, and with increased robustness to dropouts.

Availability and implementation

CCAT is part of the SCENT R-package, freely available from https://github.com/aet21/SCENT.

Supplementary information

Supplementary data are available at Bioinformatics online.

Statistical mechanics meets single-cell biology

Single-cell omics is transforming our understanding of cell biology and disease, yet the systems-level analysis and interpretation of single-cell data faces many challenges. In this Perspective, we describe the impact that fundamental concepts from statistical mechanics, notably entropy, stochastic processes and critical phenomena, are having on single-cell data analysis. We further advocate the need for more bottom-up modelling of single-cell data and to embrace a statistical mechanics analysis paradigm to help attain a deeper understanding of single-cell systems biology.

Reclassification of Facklamia ignava, Facklamia sourekii and Facklamia tabacinasalis as Falseniella ignava gen. nov., comb. nov., Hutsoniella sourekii gen. nov., comb. nov., and Ruoffia tabacinasalis gen. nov., comb. nov., and description of Ruoffia halotolerans sp. nov., isolated from hypersaline Inland Sea of Qatar

A Gram-stain-positive, non-pigmented, coccus-shaped, facultatively anaerobic and α-hemolytic bacterium designated as INB8^T was isolated from a hypersaline marine water sample collected at the Inland Sea of Qatar. The isolate was able to grow at 25-40 °C (optimum, 30 °C), at pH 5-11 and with 2-8% NaCl. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain INB8^T was placed within the family Aerococcaceae with the highest sequence similarity to Facklamia tabacinasalis CCUG 30090^T (99.5%), followed by Facklamia hominis CCUG 36813^T (93.9%), Facklamia sourekii Y17312^T (93.8%), Facklamia ignava CCUG 37419^T (93.6%), Facklamia miroungae CCUG 42728^T (93.5%), Suicoccus acidiformans ZY16052^T (93.5%), Facklamia languida CCUG 37842^T (93.2%), Ignavigranum ruoffiae (93.1%), and Dolosicoccus paucivorans DSM 15742^T (90.8%). Average nucleotide identity and digital DNA-DNA hybridization values between strain INB8^T and F. tabacinasalis CCUG 30090^T were determined to be 94.5% and 58.9% respectively, confirming strain INB8^T represents a novel species. The major fatty acids were C_14:0, C_16:0, C_18:0 and C_18:1 ω9c_. The G + C content of strain INB8^T determined from the genome was 36.3 mol%. Based on the phylogenetic, chemotaxonomic and phenotypic information, it is proposed that Facklamia tabacinasalis should be reclassified as Ruoffia tabacinasalis, Facklamia ignava be reclassified as Falseniella ignava, and Facklamia sourekii be reclassified Hutsoniella sourekii. It is further proposed that strain INB8^T should be classified as a species of the genus Ruoffia for which the name Ruoffia halotolerans sp. nov. is proposed. The type strain is INB8^T (= LMG 30291^T = CCUG 70701^T = QCC/B60/17^T).

A Machine Learning Approach to Differentiate Two Specific Breast Cancer Subtypes Using Androgen Receptor Pathway Genes

Triple-negative breast cancer is a heterogeneous disease with different molecular and histological subtypes. The Androgen receptor is expressed in a portion of triple-negative breast cancer cases and the activation of the androgen receptor pathway is thought to be a molecular subtyping signature as well as a therapeutic target for triple-negative breast cancer. Thus, identification of the androgen receptor pathway status is important for both molecular characterization andclinical management. In this study, we investigate the expression of the androgen receptor pathway in metaplastic breast cancer and luminal androgen receptor subtypes of triple-negative breast cancer and found that the androgen receptor pathway was downregulated in metaplastic breast cancer compared to luminal androgen receptor subtype. Using random forest, we found that the two subtypes of breast cancer can be molecularly classified with the gene expression of the androgen receptor pathway.

A pan-cancer survey of cell line tumor similarity by feature-weighted molecular profiles

Patient-derived cell lines are often used in pre-clinical cancer research, but some cell lines are too different from tumors to be good models. Comparison of genomic and expression profiles can guide the choice of pre-clinical models, but typically not all features are equally relevant. We present TumorComparer, a computational method for comparing cellular profiles with higher weights on functional features of interest. In this pan-cancer application, we compare ∼600 cell lines and ∼8,000 tumor samples of 24 cancer types, using weights to emphasize known oncogenic alterations. We characterize the similarity of cell lines and tumors within and across cancers by using multiple datum types and rank cell lines by their inferred quality as representative models. Beyond the assessment of cell lines, the weighted similarity approach is adaptable to patient stratification in clinical trials and personalized medicine.

TSCCA: A tensor sparse CCA method for detecting microRNA-gene patterns from multiple cancers

Existing studies have demonstrated that dysregulation of microRNAs (miRNAs or miRs) is involved in the initiation and progression of cancer. Many efforts have been devoted to identify microRNAs as potential biomarkers for cancer diagnosis, prognosis and therapeutic targets. With the rapid development of miRNA sequencing technology, a vast amount of miRNA expression data for multiple cancers has been collected. These invaluable data repositories provide new paradigms to explore the relationship between miRNAs and cancer. Thus, there is an urgent need to explore the complex cancer-related miRNA-gene patterns by integrating multi-omics data in a pan-cancer paradigm. In this study, we present a tensor sparse canonical correlation analysis (TSCCA) method for identifying cancer-related miRNA-gene modules across multiple cancers. TSCCA is able to overcome the drawbacks of existing solutions and capture both the cancer-shared and specific miRNA-gene co-expressed modules with better biological interpretations. We comprehensively evaluate the performance of TSCCA using a set of simulated data and matched miRNA/gene expression data across 33 cancer types from the TCGA database. We uncover several dysfunctional miRNA-gene modules with important biological functions and statistical significance. These modules can advance our understanding of miRNA regulatory mechanisms of cancer and provide insights into miRNA-based treatments for cancer.

MicroRNAs (miRNAs) are a class of small non-coding RNAs. Previous studies have revealed that miRNA-gene regulatory modules play key roles in the occurrence and development of cancer. However, little has been done to discover miRNA-gene regulatory modules from a pan-cancer view. Thus, it is urgently needed to develop new methods to explore the complex cancer-related miRNA-gene patterns by integrating multi-omics data of multi-cancers. To build the connections between miRNA-gene regulatory modules across different cancer types, we propose a tensor sparse canonical correlation analysis (TSCCA) method. Our specific contributions are two-fold: (1) We propose a sparse statistical learning model TSCCA and an efficient block-coordinate descent algorithm to solve it. (2) We apply TSCCA to a multi-omics data set of 33 cancer types from TCGA and identify some cancer-related miRNA-gene modules with important biological functions and statistical significance.

Bioinformatics, Molecular Docking and Experiments In Vitro Analyze the Prognostic Value of CXC Chemokines in Breast Cancer

The increasing incidence and mortality rate of Breast cancer (BC) make it a major public health problem around the world. CXC chemokines can mediate the migration of immune cells and regulate apoptosis in tumor. However, the expression and prognostic value of them in BC and their targeted drugs have not been clarified. Therefore, in this study, ONCOMINE, GEPIA2.0, UALCAN, Venny2.1.0, cBioPortal, STRING, Gene MANIA, Pathway Commons, DAVID6.8, Omicshare, Cytoscape3.6.1, TIMER2.0, Drug Bank, TCMSP, RSCBPDB, PubChem, pkCSM, Chem Draw, AutoDockTools-1.5.6 and PyMOL were utilized for analysis. The expression of CXCL1-3, CXCL9-13 between BC and normal tissues was significantly different in all the three databases. And the expression of CXCL1-2, CXCL12-13 was correlated with the stages of BC. But only CXCL1-3 were prone to mutation, and negatively correlated with survival and prognosis of BC patients. Taken together, CXCL1-2 might be therapeutic targets and biomarkers for BC patients. In addition, both of them were associated with immune infiltration. The results of molecular docking showed that Quercetin was most likely to be developed as drugs that interacted directly with CXCL1-2. And GLU29 of CXCL1, ASP-1, PRO-96, TRP-47 and LEU-45 of CXCL2 were the most potential sites, which provided valuable reference for further study of pharmacodynamics and mechanism. In addition, the inhibitory effect of Quercetin on proliferation and promoting apoptosis of BC related cell lines were confirmed in vitro. Western blot and Real-Time PCR confirmed that it increased the expression of CXCL1-2 in MDA-MB-231 and MCF-7 cells.

Drug Repurposing for COVID-19 Treatment by Integrating Network Pharmacology and Transcriptomics

Since coronavirus disease 2019 (COVID-19) is a serious new worldwide public health crisis with significant morbidity and mortality, effective therapeutic treatments are urgently needed. Drug repurposing is an efficient and cost-effective strategy with minimum risk for identifying novel potential treatment options by repositioning therapies that were previously approved for other clinical outcomes. Here, we used an integrated network-based pharmacologic and transcriptomic approach to screen drug candidates novel for COVID-19 treatment. Network-based proximity scores were calculated to identify the drug-disease pharmacological effect between drug-target relationship modules and COVID-19 related genes. Gene set enrichment analysis (GSEA) was then performed to determine whether drug candidates influence the expression of COVID-19 related genes and examine the sensitivity of the repurposing drug treatment to peripheral immune cell types. Moreover, we used the complementary exposure model to recommend potential synergistic drug combinations. We identified 18 individual drug candidates including nicardipine, orantinib, tipifarnib and promethazine which have not previously been proposed as possible treatments for COVID-19. Additionally, 30 synergistic drug pairs were ultimately recommended including fostamatinib plus tretinoin and orantinib plus valproic acid. Differential expression genes of most repurposing drugs were enriched significantly in B cells. The findings may potentially accelerate the discovery and establishment of an effective therapeutic treatment plan for COVID-19 patients.

Histamine, Metabolic Remodelling and Angiogenesis: A Systems Level Approach

Histamine is a highly pleiotropic biogenic amine involved in key physiological processes including neurotransmission, immune response, nutrition, and cell growth and differentiation. Its effects, sometimes contradictory, are mediated by at least four different G-protein coupled receptors, which expression and signalling pathways are tissue-specific. Histamine metabolism conforms a very complex network that connect many metabolic processes important for homeostasis, including nitrogen and energy metabolism. This review brings together and analyses the current information on the relationships of the "histamine system" with other important metabolic modules in human physiology, aiming to bridge current information gaps. In this regard, the molecular characterization of the role of histamine in the modulation of angiogenesis-mediated processes, such as cancer, makes a promising research field for future biomedical advances.

A Detailed Catalogue of Multi-Omics Methodologies for Identification of Putative Biomarkers and Causal Molecular Networks in Translational Cancer Research

Recent advances in sequencing and biotechnological methodologies have led to the generation of large volumes of molecular data of different omics layers, such as genomics, transcriptomics, proteomics and metabolomics. Integration of these data with clinical information provides new opportunities to discover how perturbations in biological processes lead to disease. Using data-driven approaches for the integration and interpretation of multi-omics data could stably identify links between structural and functional information and propose causal molecular networks with potential impact on cancer pathophysiology. This knowledge can then be used to improve disease diagnosis, prognosis, prevention, and therapy. This review will summarize and categorize the most current computational methodologies and tools for integration of distinct molecular layers in the context of translational cancer research and personalized therapy. Additionally, the bioinformatics tools Multi-Omics Factor Analysis (MOFA) and netDX will be tested using omics data from public cancer resources, to assess their overall robustness, provide reproducible workflows for gaining biological knowledge from multi-omics data, and to comprehensively understand the significantly perturbed biological entities in distinct cancer types. We show that the performed supervised and unsupervised analyses result in meaningful and novel findings.

An Update on MRMAssayDB: A Comprehensive Resource for Targeted Proteomics Assays in the Community

Precise multiplexed quantification of proteins in biological samples can be achieved by targeted proteomics using multiple or parallel reaction monitoring (MRM/PRM). Combined with internal standards, the method achieves very good repeatability and reproducibility enabling excellent protein quantification and allowing longitudinal and cohort studies. A laborious part of performing such experiments lies in the preparation steps dedicated to the development and validation of individual protein assays. Several public repositories host information on targeted proteomics assays, including NCI’s Clinical Proteomic Tumor Analysis Consortium assay portals, PeptideAtlas SRM Experiment Library, SRMAtlas, PanoramaWeb, and PeptideTracker, with all offering varying levels of details. We introduced MRMAssayDB in 2018 as an integrated resource for targeted proteomics assays. The Web-based application maps and links the assays from the repositories, includes comprehensive up-to-date protein and sequence annotations, and provides multiple visualization options on the peptide and protein level. We have extended MRMAssayDB with more assays and extensive annotations. Currently it contains >828 000 assays covering >51 000 proteins from 94 organisms, of which >17 000 proteins are present in >2400 biological pathways, and >48 000 mapping to >21 000 Gene Ontology terms. This is an increase of about four times the number of assays since introduction. We have expanded annotations of interaction, biological pathways, and disease associations. A newly added visualization module for coupled molecular structural annotation browsing allows the user to interactively examine peptide sequence and any known PTMs and disease mutations, and map all to available protein 3D structures. Because of its integrative approach, MRMAssayDB enables a holistic view of suitable proteotypic peptides and commonly used transitions in empirical data. Availability: http://mrmassaydb.proteincentre.com.

Artificial Intelligence and Early Detection of Pancreatic Cancer: 2020 Summative Review

ABSTRACT

Despite considerable research efforts, pancreatic cancer is associated with a dire prognosis and a 5-year survival rate of only 10%. Early symptoms of the disease are mostly nonspecific. The premise of improved survival through early detection is that more individuals will benefit from potentially curative treatment. Artificial intelligence (AI) methodology has emerged as a successful tool for risk stratification and identification in general health care. In response to the maturity of AI, Kenner Family Research Fund conducted the 2020 AI and Early Detection of Pancreatic Cancer Virtual Summit (www.pdac-virtualsummit.org) in conjunction with the American Pancreatic Association, with a focus on the potential of AI to advance early detection efforts in this disease. This comprehensive presummit article was prepared based on information provided by each of the interdisciplinary participants on one of the 5 following topics: Progress, Problems, and Prospects for Early Detection; AI and Machine Learning; AI and Pancreatic Cancer-Current Efforts; Collaborative Opportunities; and Moving Forward-Reflections from Government, Industry, and Advocacy. The outcome from the robust Summit conversations, to be presented in a future white paper, indicate that significant progress must be the result of strategic collaboration among investigators and institutions from multidisciplinary backgrounds, supported by committed funders.

Automating parameter selection to avoid implausible biological pathway models

A common way to integrate and analyze large amounts of biological "omic" data is through pathway reconstruction: using condition-specific omic data to create a subnetwork of a generic background network that represents some process or cellular state. A challenge in pathway reconstruction is that adjusting pathway reconstruction algorithms' parameters produces pathways with drastically different topological properties and biological interpretations. Due to the exploratory nature of pathway reconstruction, there is no ground truth for direct evaluation, so parameter tuning methods typically used in statistics and machine learning are inapplicable. We developed the pathway parameter advising algorithm to tune pathway reconstruction algorithms to minimize biologically implausible predictions. We leverage background knowledge in pathway databases to select pathways whose high-level structure resembles that of manually curated biological pathways. At the core of this method is a graphlet decomposition metric, which measures topological similarity to curated biological pathways. In order to evaluate pathway parameter advising, we compare its performance in avoiding implausible networks and reconstructing pathways from the NetPath database with other parameter selection methods across four pathway reconstruction algorithms. We also demonstrate how pathway parameter advising can guide reconstruction of an influenza host factor network. Pathway parameter advising is method agnostic; it is applicable to any pathway reconstruction algorithm with tunable parameters.

GeneWalk identifies relevant gene functions for a biological context using network representation learning

A bottleneck in high-throughput functional genomics experiments is identifying the most important genes and their relevant functions from a list of gene hits. Gene Ontology (GO) enrichment methods provide insight at the gene set level. Here, we introduce GeneWalk ( github.com/churchmanlab/genewalk ) that identifies individual genes and their relevant functions critical for the experimental setting under examination. After the automatic assembly of an experiment-specific gene regulatory network, GeneWalk uses representation learning to quantify the similarity between vector representations of each gene and its GO annotations, yielding annotation significance scores that reflect the experimental context. By performing gene- and condition-specific functional analysis, GeneWalk converts a list of genes into data-driven hypotheses.

CellMiner Cross-Database (CellMinerCDB) version 1.2: Exploration of patient-derived cancer cell line pharmacogenomics

CellMiner Cross-Database (CellMinerCDB, discover.nci.nih.gov/cellminercdb) allows integration and analysis of molecular and pharmacological data within and across cancer cell line datasets from the National Cancer Institute (NCI), Broad Institute, Sanger/MGH and MD Anderson Cancer Center (MDACC). We present CellMinerCDB 1.2 with updates to datasets from NCI-60, Broad Cancer Cell Line Encyclopedia and Sanger/MGH, and the addition of new datasets, including NCI-ALMANAC drug combination, MDACC Cell Line Project proteomic, NCI-SCLC DNA copy number and methylation data, and Broad methylation, genetic dependency and metabolomic datasets. CellMinerCDB (v1.2) includes several improvements over the previously published version: (i) new and updated datasets; (ii) support for pattern comparisons and multivariate analyses across data sources; (iii) updated annotations with drug mechanism of action information and biologically relevant multigene signatures; (iv) analysis speedups via caching; (v) a new dataset download feature; (vi) improved visualization of subsets of multiple tissue types; (vii) breakdown of univariate associations by tissue type; and (viii) enhanced help information. The curation and common annotations (e.g. tissues of origin and identifiers) provided here across pharmacogenomic datasets increase the utility of the individual datasets to address multiple researcher question types, including data reproducibility, biomarker discovery and multivariate analysis of drug activity.

Multi-omics integration in biomedical research - A metabolomics-centric review

Recent advances in high-throughput technologies have enabled the profiling of multiple layers of a biological system, including DNA sequence data (genomics), RNA expression levels (transcriptomics), and metabolite levels (metabolomics). This has led to the generation of vast amounts of biological data that can be integrated in so-called multi-omics studies to examine the complex molecular underpinnings of health and disease. Integrative analysis of such datasets is not straightforward and is particularly complicated by the high dimensionality and heterogeneity of the data and by the lack of universal analysis protocols. Previous reviews have discussed various strategies to address the challenges of data integration, elaborating on specific aspects, such as network inference or feature selection techniques. Thereby, the main focus has been on the integration of two omics layers in their relation to a phenotype of interest. In this review we provide an overview over a typical multi-omics workflow, focusing on integration methods that have the potential to combine metabolomics data with two or more omics. We discuss multiple integration concepts including data-driven, knowledge-based, simultaneous and step-wise approaches. We highlight the application of these methods in recent multi-omics studies, including large-scale integration efforts aiming at a global depiction of the complex relationships within and between different biological layers without focusing on a particular phenotype.

Drug2ways: Reasoning over causal paths in biological networks for drug discovery

Elucidating the causal mechanisms responsible for disease can reveal potential therapeutic targets for pharmacological intervention and, accordingly, guide drug repositioning and discovery. In essence, the topology of a network can reveal the impact a drug candidate may have on a given biological state, leading the way for enhanced disease characterization and the design of advanced therapies. Network-based approaches, in particular, are highly suited for these purposes as they hold the capacity to identify the molecular mechanisms underlying disease. Here, we present drug2ways, a novel methodology that leverages multimodal causal networks for predicting drug candidates. Drug2ways implements an efficient algorithm which reasons over causal paths in large-scale biological networks to propose drug candidates for a given disease. We validate our approach using clinical trial information and demonstrate how drug2ways can be used for multiple applications to identify: i) single-target drug candidates, ii) candidates with polypharmacological properties that can optimize multiple targets, and iii) candidates for combination therapy. Finally, we make drug2ways available to the scientific community as a Python package that enables conducting these applications on multiple standard network formats.

Measurement Recorder: developing a useful tool for making species descriptions that produces computable phenotypes

To use published phenotype information in computational analyses, there have been efforts to convert descriptions of phenotype characters from human languages to ontologized statements. This postpublication curation process is not only slow and costly, it is also burdened with significant intercurator variation (including curator-author variation), due to different interpretations of a character by various individuals. This problem is inherent in any human-based intellectual activity. To address this problem, making scientific publications semantically clear (i.e. computable) by the authors at the time of publication is a critical step if we are to avoid postpublication curation. To help authors efficiently produce species phenotypes while producing computable data, we are experimenting with an author-driven ontology development approach and developing and evaluating a series of ontology-aware software modules that would create publishable species descriptions that are readily useable in scientific computations. The first software module prototype called Measurement Recorder has been developed to assist authors in defining continuous measurements and reported in this paper. Two usability studies of the software were conducted with 22 undergraduate students majoring in information science and 32 in biology. Results suggest that participants can use Measurement Recorder without training and they find it easy to use after limited practice. Participants also appreciate the semantic enhancement features. Measurement Recorder's character reuse features facilitate character convergence among participants by 48% and have the potential to further reduce user errors in defining characters. A set of software design issues have also been identified and then corrected. Measurement Recorder enables authors to record measurements in a semantically clear manner and enriches phenotype ontology along the way. Future work includes representing the semantic data as Resource Description Framework (RDF) knowledge graphs and characterizing the division of work between authors as domain knowledge providers and ontology engineers as knowledge formalizers in this new author-driven ontology development approach.

Phosphoproteomic quantitation and causal analysis reveal pathways in GPVI/ITAM-mediated platelet activation programs

Platelets engage cues of pending vascular injury through coordinated adhesion, secretion, and aggregation responses. These rapid, progressive changes in platelet form and function are orchestrated downstream of specific receptors on the platelet surface and through intracellular signaling mechanisms that remain systematically undefined. This study brings together cell physiological and phosphoproteomics methods to profile signaling mechanisms downstream of the immunotyrosine activation motif (ITAM) platelet collagen receptor GPVI. Peptide tandem mass tag (TMT) labeling, sample multiplexing, synchronous precursor selection (SPS), and triple stage tandem mass spectrometry (MS3) detected >3000 significant (false discovery rate < 0.05) phosphorylation events on >1300 proteins over conditions initiating and progressing GPVI-mediated platelet activation. With literature-guided causal inference tools, >300 site-specific signaling relations were mapped from phosphoproteomics data among key and emerging GPVI effectors (ie, FcRγ, Syk, PLCγ2, PKCδ, DAPP1). Through signaling validation studies and functional screening, other less-characterized targets were also considered within the context of GPVI/ITAM pathways, including Ras/MAPK axis proteins (ie, KSR1, SOS1, STAT1, Hsp27). Highly regulated GPVI/ITAM targets out of context of curated knowledge were also illuminated, including a system of >40 Rab GTPases and associated regulatory proteins, where GPVI-mediated Rab7 S72 phosphorylation and endolysosomal maturation were blocked by TAK1 inhibition. In addition to serving as a model for generating and testing hypotheses from omics datasets, this study puts forth a means to identify hemostatic effectors, biomarkers, and therapeutic targets relevant to thrombosis, vascular inflammation, and other platelet-associated disease states.

Hypomethylation of GDNF family receptor alpha 1 promotes epithelial-mesenchymal transition and predicts metastasis of colorectal cancer

Tumor metastasis is the major cause of poor prognosis and mortality in colorectal cancer (CRC). However, early diagnosis of highly metastatic CRC is currently difficult. In the present study, we screened for a novel biomarker, GDNF family receptor alpha 1 (GFRA1) based on the expression and methylation data in CRC patients from The Cancer Genome Altlas (TCGA), followed by further analysis of the correlation between the GFRA1 expression, methylation, and prognosis of patients. Our results show DNA hypomethylation-mediated upregulation of GFRA1 in invasive CRC, and it was found to be correlated with poor prognosis of CRC patients. Furthermore, GFRA1 methylation-modified sequences were found to have potential as methylation diagnostic markers of highly metastatic CRC. The targeted demethylation of GFRA1 by dCas9-TET1CD and gRNA promoted CRC metastasis in vivo and in vitro. Mechanistically, demethylation of GFRA1 induces epithelial-mesenchymal transition (EMT) by promoting AKT phosphorylation and increasing c-Jun expression in CRC cells. Collectively, our findings indicate that GFRA1 hypomethylation can promote CRC invasion via inducing EMT, and thus, GFRA1 methylation can be used as a biomarker for the early diagnosis of highly metastasis CRC.

Abnormal DNA methylation, one of important characteristics in tumor cells, is exploited as biomarkers for cancer diagnosis and prognosis prediction. Early diagnosis of highly metastatic CRC will be helpful for the clinical management, thus prolongs patient survival. However, it is currently difficult to make early diagnosis of highly metastatic CRC in clinical practice. Currently, we screened a novel biomarker gene, GFRA1, which associated with the invasion and poor prognosis of CRC. The targeted demethylation of GFRA1 exerted a significant promoting effect on CRC metastasis, and GFRA1 methylation-modified sequences are valuable diagnostic biomarker for CRC metastasis risk assessment. Mechanically, demethylation of GFRA1 induced epithelial-mesenchymal transition (EMT) by upregulating AKT phosphorylation and c-Jun expression in CRC cells. Our results demonstrate the promoting effect of GFRA1 demethylation on CRC invasion and GFRA1 methylation may be a potential prognostic marker for predicting metastasis of CRC.

Pathway information extracted from 25 years of pathway figures

Thousands of pathway diagrams are published each year as static figures inaccessible to computational queries and analyses. Using a combination of machine learning, optical character recognition, and manual curation, we identified 64,643 pathway figures published between 1995 and 2019 and extracted 1,112,551 instances of human genes, comprising 13,464 unique NCBI genes, participating in a wide variety of biological processes. This collection represents an order of magnitude more genes than found in the text of the same papers, and thousands of genes missing from other pathway databases, thus presenting new opportunities for discovery and research.

Uncovering cancer gene regulation by accurate regulatory network inference from uninformative data

The interactions among the components of a living cell that constitute the gene regulatory network (GRN) can be inferred from perturbation-based gene expression data. Such networks are useful for providing mechanistic insights of a biological system. In order to explore the feasibility and quality of GRN inference at a large scale, we used the L1000 data where ~1000 genes have been perturbed and their expression levels have been quantified in 9 cancer cell lines. We found that these datasets have a very low signal-to-noise ratio (SNR) level causing them to be too uninformative to infer accurate GRNs. We developed a gene reduction pipeline in which we eliminate uninformative genes from the system using a selection criterion based on SNR, until reaching an informative subset. The results show that our pipeline can identify an informative subset in an overall uninformative dataset, allowing inference of accurate subset GRNs. The accurate GRNs were functionally characterized and potential novel cancer-related regulatory interactions were identified.

OTUB1 Recruits Tumor Infiltrating Lymphocytes and Is a Prognostic Marker in Digestive Cancers

Background

The deubiquitinating enzyme (DUB) OTUB1 can regulate the process of ubiquitination, but the influence of OTUB1 on immunity, apoptosis, autophagy, and the prognosis of digestive cancers requires further exploration.

Methods

OTUB1 expression was analyzed with the Oncomine and TIMER database. Kaplan-Meier plotter was used to calculate the association between OTUB1 and clinical prognosis. The regulation of OTUB1 on cancer immunocyte infiltration was determined by the TIMER database. The interaction between OTUB1 and immune genes, gene expression profiling (GEP), key genes of apoptosis and autophagy were analyzed via GEPIA. Protein-protein interaction (PPI), gene expression profiling (GEP), and functional pathway enrichment were also performed with the STRING and Pathway Common databases, respectively.

Results

High OTUB1 expression was found in CHOL, LIHC, READ, ESCA, and COAD, which was significantly associated with the poorer OS of LIHC (HR = 2.07, 95% CI = 1.30-3.30, P = 0.002), with modifications by sex, stage, grade, and mutant burden. OTUB1 can promote the recruitment of B cells, CD8 + T cells, macrophages in ESCA, B cells, and neutrophils in LIHC. We determined a significant interaction between OTUB1 and USP8, RNF128, LRIG1, UBB, UBC, STAM2, RNF41, EGFR, RPS27A, and HGS by PPI. This functional pathway indicates the regulatory role of OTUB1on immune, apoptosis, and autophagy through its interaction with TP53 and ATG.

Conclusions

OTUB1 performed as a molecular indicator of poor prognosis in digestive cancers, regulated the infiltration of tumor immunocytes, and exerted a significant influence on apoptosis and autophagy. OTUB1 is a potential antitumor target for digestive tumors.

netboxr: Automated discovery of biological process modules by network analysis in R

SUMMARY

Large-scale sequencing projects, such as The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC), have generated high throughput sequencing and molecular profiling data sets, but it is still challenging to identify potentially causal changes in cellular processes in cancer as well as in other diseases in an automated fashion. We developed the netboxr package written in the R programming language, which makes use of the NetBox algorithm to identify candidate cancer-related functional modules. The algorithm makes use of a data-driven, network-based approach that combines prior knowledge with a network clustering algorithm, obviating the need for and the limitation of independently curated functionally labeled gene sets. The method can combine multiple data types, such as mutations and copy number alterations, leading to more reliable identification of functional modules. We make the tool available in the Bioconductor R ecosystem for applications in cancer research and cell biology.

AVAILABILITY AND IMPLEMENTATION

The netboxr package is free and open-sourced under the GNU GPL-3 license R package available at https://www.bioconductor.org/packages/release/bioc/html/netboxr.html.

Functional in vivo and in vitro effects of 20q11.21 genetic aberrations on hPSC differentiation

Human pluripotent stem cells (hPSCs) have promising therapeutic applications due to their infinite capacity for self-renewal and pluripotency. Genomic stability is imperative for the clinical use of hPSCs; however, copy number variation (CNV), especially recurrent CNV at 20q11.21, may contribute genomic instability of hPSCs. Furthermore, the effects of CNVs in hPSCs at the whole-transcriptome scale are poorly understood. This study aimed to examine the functional in vivo and in vitro effects of frequently detected CNVs at 20q11.21 during early-stage differentiation of hPSCs. Comprehensive transcriptome profiling of abnormal hPSCs revealed that the differential gene expression patterns had a negative effect on differentiation potential. Transcriptional heterogeneity identified by single-cell RNA sequencing (scRNA-seq) of embryoid bodies from two different isogenic lines of hPSCs revealed alterations in differentiated cell distributions compared with that of normal cells. RNA-seq analysis of 22 teratomas identified several differentially expressed lineage-specific markers in hPSCs with CNVs, consistent with the histological results of the altered ecto/meso/endodermal ratio due to CNVs. Our results suggest that CNV amplification contributes to cell proliferation, apoptosis, and cell fate specification. This work shows the functional consequences of recurrent genetic abnormalities and thereby provides evidence to support the development of cell-based applications.

Comprehensive analysis of the expression and prognostic value of CXC chemokines in colorectal cancer

The C-X-C motif (CXC) chemokines play an important role in inflammatory processes and angiogenesis and are also associated with tumor development, progression and metastasis. They can be either promoting or inhibiting factors in colorectal cancers (CRC). The expression patterns and prognostic values of the CXC family still need further investigation. In this study, we investigated data related to transcription, translation, survival and tumor immune infiltration for CXC chemokines in patients with CRC from the ONCOMINE, GEPIA, cBioPortal, HPA and TIMER databases. We found that the expression levels of CXCL1-3, CXCL5, and CXCL8 were higher in CRC tissues than in colorectal tissues. Expression among stages significantly varied for CXCL1-3 and CXCL9-11. The survival analysis revealed that high transcriptional levels of CXCL4 and CXCL9-11 could serve as positive prognostic factors for patients with CRC. CXCL9-11 were highly associated with CD8+ T cells and natural killer (NK) cells in the tumor immune infiltration analysis, indicating their role in the antitumor immune response. This study implies that CXCL1-3, CXCL5, and CXCL8 are important factors during CRC oncogenesis and that CXCL9-11 could be new biomarkers for the prognosis of CRC.

The tissue proteome in the multi-omic landscape of kidney disease

Kidney research is entering an era of 'big data' and molecular omics data can provide comprehensive insights into the molecular footprints of cells. In contrast to transcriptomics, proteomics and metabolomics generate data that relate more directly to the pathological symptoms and clinical parameters observed in patients. Owing to its complexity, the proteome still holds many secrets, but has great potential for the identification of drug targets. Proteomics can provide information about protein synthesis, modification and degradation, as well as insight into the physical interactions between proteins, and between proteins and other biomolecules. Thus far, proteomics in nephrology has largely focused on the discovery and validation of biomarkers, but the systematic analysis of the nephroproteome can offer substantial additional insights, including the discovery of mechanisms that trigger and propagate kidney disease. Moreover, proteome acquisition might provide a diagnostic tool that complements the assessment of a kidney biopsy sample by a pathologist. Such applications are becoming increasingly feasible with the development of high-throughput and high-coverage technologies, such as versatile mass spectrometry-based techniques and protein arrays, and encourage further proteomics research in nephrology.

Transcriptomic Profiling for the Autophagy Pathway in Colorectal Cancer

The role of autophagy in colorectal cancer (CRC) pathogenesis appears to be crucial. Autophagy acts both as a tumor suppressor, by removing redundant cellular material, and a tumor-promoting factor, by providing access to components necessary for growth, metabolism, and proliferation. To date, little is known about the expression of genes that play a basal role in the autophagy in CRC. In this study, we aimed to compare the expression levels of 46 genes involved in the autophagy pathway between tumor-adjacent and tumor tissue, employing large RNA sequencing (RNA-seq) and microarray datasets. Additionally, we verified our results using data on 38 CRC cell lines. Gene set enrichment analysis revealed a significant deregulation of autophagy-related gene sets in CRC. The unsupervised clustering of tumors using the mRNA levels of autophagy-related genes revealed the existence of two major clusters: microsatellite instability (MSI)-enriched and -depleted. In cluster 1 (MSI-depleted), ATG9B and LAMP1 genes were the most prominently expressed, whereas cluster 2 (MSI-enriched) was characterized by DRAM1 upregulation. CRC cell lines were also clustered according to MSI-enriched/-depleted subgroups. The moderate deregulation of autophagy-related genes in cancer tissue, as compared to adjacent tissue, suggests a prominent field cancerization or early disruption of autophagy. Genes differentiating these clusters are promising candidates for CRC targeting therapy worthy of further investigation.

Designing custom CRISPR libraries for hypothesis-driven drug target discovery

Over the last decade Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) has been developed into a potent molecular biology tool used to rapidly modify genes or their expression in a multitude of ways. In parallel, CRISPR-based screening approaches have been developed as powerful discovery platforms for dissecting the genetic basis of cellular behavior, as well as for drug target discovery. CRISPR screens can be designed in numerous ways. Here, we give a brief background to CRISPR screens and discuss the pros and cons of different design approaches, including unbiased genome-wide screens that target all known genes, as well as hypothesis-driven custom screens in which selected subsets of genes are targeted (Fig. 1). We provide several suggestions for how a custom screen can be designed, which could broadly serve as inspiration for any experiment that includes candidate gene selection. Finally, we discuss how results from CRISPR screens could be translated into drug development, as well as future trends we foresee in the rapidly evolving CRISPR screen field.

CCTs as new biomarkers for the prognosis of head and neck squamous cancer

The chaperonin-containing T-complex protein 1 (CCT) subunits participate in diverse diseases. However, little is known about their expression and prognostic values in human head and neck squamous cancer (HNSC). This article aims to evaluate the effects of CCT subunits regarding their prognostic values for HNSC. We mined the transcriptional and survival data of CCTs in HNSC patients from online databases. A protein-protein interaction network was constructed and a functional enrichment analysis of target genes was performed. We observed that the mRNA expression levels of CCT1/2/3/4/5/6/7/8 were higher in HNSC tissues than in normal tissues. Survival analysis revealed that the high mRNA transcriptional levels of CCT3/4/5/6/7/8 were associated with a low overall survival. The expression levels of CCT4/7 were correlated with advanced tumor stage. And the overexpression of CCT4 was associated with higher N stage of patients. Validation of CCTs' differential expression and prognostic values was achieved by the Human Protein Atlas and GEO datasets. Mechanistic exploration of CCT subunits by the functional enrichment analysis suggests that these genes may influence the HNSC prognosis by regulating PI3K-Akt and other pathways. This study implies that CCT3/4/6/7/8 are promising biomarkers for the prognosis of HNSC.

The first 10 years of the international coordination network for standards in systems and synthetic biology (COMBINE)

This paper presents a report on outcomes of the 10th Computational Modeling in Biology Network (COMBINE) meeting that was held in Heidelberg, Germany, in July of 2019. The annual event brings together researchers, biocurators and software engineers to present recent results and discuss future work in the area of standards for systems and synthetic biology. The COMBINE initiative coordinates the development of various community standards and formats for computational models in the life sciences. Over the past 10 years, COMBINE has brought together standard communities that have further developed and harmonized their standards for better interoperability of models and data. COMBINE 2019 was co-located with a stakeholder workshop of the European EU-STANDS4PM initiative that aims at harmonized data and model standardization for in silico models in the field of personalized medicine, as well as with the FAIRDOM PALs meeting to discuss findable, accessible, interoperable and reusable (FAIR) data sharing. This report briefly describes the work discussed in invited and contributed talks as well as during breakout sessions. It also highlights recent advancements in data, model, and annotation standardization efforts. Finally, this report concludes with some challenges and opportunities that this community will face during the next 10 years.

Comparative Pathway Integrator: A Framework of Meta-Analytic Integration of Multiple Transcriptomic Studies for Consensual and Differential Pathway Analysis

Pathway enrichment analysis provides a knowledge-driven approach to interpret differentially expressed genes associated with disease status. Many tools have been developed to analyze a single study. However, when multiple studies of different conditions are jointly analyzed, novel integrative tools are needed. In addition, pathway redundancy introduced by combining multiple public pathway databases hinders interpretation and knowledge discovery. We present a meta-analytic integration tool, Comparative Pathway Integrator (CPI), to address these issues using adaptively weighted Fisher's method to discover consensual and differential enrichment patterns, a tight clustering algorithm to reduce pathway redundancy, and a text mining algorithm to assist interpretation of the pathway clusters. We applied CPI to jointly analyze six psychiatric disorder transcriptomic studies to demonstrate its effectiveness, and found functions confirmed by previous biological studies as well as novel enrichment patterns. CPI's R package is accessible online on Github metaOmics/MetaPath.

Probing Protein Interaction Networks by Combining MS-Based Proteomics and Structural Data Integration

Protein-protein interactions play a major role in the molecular machinery of life, and various techniques such as AP-MS are dedicated to their identification. However, those techniques return lists of proteins devoid of organizational structure, not detailing which proteins interact with which others. Proposing a hierarchical view of the interactions between the members of the flat list becomes highly tedious for large data sets when done by hand. To help hierarchize this data, we introduce a new bioinformatics protocol that integrates information of the multimeric protein 3D structures available in the Protein Data Bank using remote homology detection, as well as information related to Short Linear Motifs and interaction data from the BioGRID. We illustrate on two unrelated use-cases of different complexity how our approach can be useful to decipher the network of interactions hidden in the list of input proteins, and how it provides added value compared to state-of-the-art resources such as Interactome3D or STRING. Particularly, we show the added value of using homology detection to distinguish between orthologs and paralogs, and to distinguish between core obligate and more facultative interactions. We also demonstrate the potential of considering interactions occurring through Short Linear Motifs.