Reactome knowledgebase of human biological pathways and processes

Inflammatory gene regulation by Cdc42 in airway epithelial cells

Cytokine release from airway epithelial cells is a key immunological process that coordinates an immune response in the lungs. We propose that the Rho GTPase, Cdc42, regulates both transcription and trafficking of cytokines, ultimately affecting the essential process of cytokine release and subsequent inflammation in the lungs. Here, we examined the pro-inflammatory transcriptional profile that occurs in bronchial epithelial cells (BEAS-2B) in response to TNF-α using RNA-Seq and differential gene expression analysis. To interrogate the role of Cdc42 in inflammatory gene expression, we used a pharmacological inhibitor of Cdc42, ML141, and determined changes in the transcriptomic profile induced by Cdc42 inhibition. Our results indicated that Cdc42 inhibition with ML141 resulted in a unique inflammatory phenotype concomitant with increased gene expression of ER stress genes, Golgi membrane and vesicle transport genes. To further interrogate the inflammatory pathways regulated by Cdc42, we made BEAS-2B knockdown strains for the signaling targets TRIB3, DUSP5, SESN2 and BMP4, which showed high differential expression in response to Cdc42 inhibition. Depletion of DUSP5 and TRIB3 reduced the pro-inflammatory response triggered by Cdc42 inhibition as shown by a reduction in cytokine transcript levels. Depletion of SESN2 and BMP4 did not affect cytokine transcript level, however, Golgi fragmentation was reduced. These results provide further evidence that in airway epithelial cells, Cdc42 is part of a signaling network that controls inflammatory gene expression and secretion by regulating Golgi integrity. Summary sentence:We define the Cdc42-regulated gene networks for inflammatory signaling in airway epithelial cells which includes regulation of ER stress response and vesicle trafficking pathways.

CCPA: cloud-based, self-learning modules for consensus pathway analysis using GO, KEGG and Reactome

This manuscript describes the development of a resource module that is part of a learning platform named 'NIGMS Sandbox for Cloud-based Learning' (https://github.com/NIGMS/NIGMS-Sandbox). The module delivers learning materials on Cloud-based Consensus Pathway Analysis in an interactive format that uses appropriate cloud resources for data access and analyses. Pathway analysis is important because it allows us to gain insights into biological mechanisms underlying conditions. But the availability of many pathway analysis methods, the requirement of coding skills, and the focus of current tools on only a few species all make it very difficult for biomedical researchers to self-learn and perform pathway analysis efficiently. Furthermore, there is a lack of tools that allow researchers to compare analysis results obtained from different experiments and different analysis methods to find consensus results. To address these challenges, we have designed a cloud-based, self-learning module that provides consensus results among established, state-of-the-art pathway analysis techniques to provide students and researchers with necessary training and example materials. The training module consists of five Jupyter Notebooks that provide complete tutorials for the following tasks: (i) process expression data, (ii) perform differential analysis, visualize and compare the results obtained from four differential analysis methods (limma, t-test, edgeR, DESeq2), (iii) process three pathway databases (GO, KEGG and Reactome), (iv) perform pathway analysis using eight methods (ORA, CAMERA, KS test, Wilcoxon test, FGSEA, GSA, SAFE and PADOG) and (v) combine results of multiple analyses. We also provide examples, source code, explanations and instructional videos for trainees to complete each Jupyter Notebook. The module supports the analysis for many model (e.g. human, mouse, fruit fly, zebra fish) and non-model species. The module is publicly available at https://github.com/NIGMS/Consensus-Pathway-Analysis-in-the-Cloud. This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [1] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

Transcriptomic Profile of Breast Tissue of Premenopausal Women Following Treatment with Progesterone Receptor Modulator: Secondary Outcomes of a Randomized Controlled Trial

Progesterone receptor antagonism is gaining attention due to progesterone's recognized role as a major mitogen in breast tissue. Limited but promising data suggest the potential efficacy of antiprogestins in breast cancer prevention. The present study presents secondary outcomes from a randomized controlled trial and examines changes in breast mRNA expression following mifepristone treatment in healthy premenopausal women. We analyzed 32 paired breast biopsies from 16 women at baseline and after two months of mifepristone treatment. In total, 27 differentially expressed genes were identified, with enriched biological functions related to extracellular matrix remodeling. Notably, the altered gene signature induced by mifepristone in vivo was rather similar to the in vitro signature. Furthermore, this gene expression signature was linked to breast carcinogenesis and notably linked with progesterone receptor expression status in breast cancer, as validated in The Cancer Genome Atlas dataset using the R2 platform. The present study is the first to explore the breast transcriptome following mifepristone treatment in normal breast tissue in vivo, enhancing the understanding of progesterone receptor antagonism and its potential protective effect against breast cancer.

Network pharmacology: an efficient but underutilized approach in oral, head and neck cancer therapy-a review

The application of network pharmacology (NP) has advanced our understanding of the complex molecular mechanisms underlying diseases, including neck, head, and oral cancers, as well as thyroid carcinoma. This review aimed to explore the therapeutic potential of natural network pharmacology using compounds and traditional Chinese medicines for combating these malignancies. NP serves as a pivotal tool that provides a comprehensive view of the interactions among compounds, genes, and diseases, thereby contributing to the advancement of disease treatment and management. In parallel, this review discusses the significance of publicly accessible databases in the identification of oral, head, and neck cancer-specific genes. These databases, including those for head and neck oral cancer, head and neck cancer, oral cancer, and genomic variants of oral cancer, offer valuable insights into the genes, miRNAs, drugs, and genetic variations associated with these cancers. They serve as indispensable resources for researchers, clinicians, and drug developers, contributing to the pursuit of precision medicine and improved treatment of these challenging malignancies. In summary, advancements in NP could improve the globalization and modernization of traditional medicines and prognostic targets as well as aid in the development of innovative drugs. Furthermore, this review will be an eye-opener for researchers working on drug development from traditional medicines by applying NP approaches.

A critical assessment of Traditional Chinese Medicine databases as a source for drug discovery

Traditional Chinese Medicine (TCM) has been used for thousands of years to treat human diseases. Recently, many databases have been devoted to studying TCM pharmacology. Most of these databases include information about the active ingredients of TCM herbs and their disease indications. These databases enable researchers to interrogate the mechanisms of action of TCM systematically. However, there is a need for comparative studies of these databases, as they are derived from various resources with different data processing methods. In this review, we provide a comprehensive analysis of the existing TCM databases. We found that the information complements each other by comparing herbs, ingredients, and herb-ingredient pairs in these databases. Therefore, data harmonization is vital to use all the available information fully. Moreover, different TCM databases may contain various annotation types for herbs or ingredients, notably for the chemical structure of ingredients, making it challenging to integrate data from them. We also highlight the latest TCM databases on symptoms or gene expressions, suggesting that using multi-omics data and advanced bioinformatics approaches may provide new insights for drug discovery in TCM. In summary, such a comparative study would help improve the understanding of data complexity that may ultimately motivate more efficient and more standardized strategies towards the digitalization of TCM.

Predicting pathways for old and new metabolites through clustering

The diverse metabolic pathways are fundamental to all living organisms, as they harvest energy, synthesize biomass components, produce molecules to interact with the microenvironment, and neutralize toxins. While the discovery of new metabolites and pathways continues, the prediction of pathways for new metabolites can be challenging. It can take vast amounts of time to elucidate pathways for new metabolites; thus, according to HMDB (Human Metabolome Database), only 60% of metabolites get assigned to pathways. Here, we present an approach to identify pathways based on metabolite structure. We extracted 201 features from SMILES annotations and identified new metabolites from PubMed abstracts and HMDB. After applying clustering algorithms to both groups of features, we quantified correlations between metabolites, and found the clusters accurately linked 92% of known metabolites to their respective pathways. Thus, this approach could be valuable for predicting metabolic pathways for new metabolites.

Benchmark Dataset for Training Machine Learning Models to Predict the Pathway Involvement of Metabolites

Metabolic pathways are a human-defined grouping of life sustaining biochemical reactions, metabolites being both the reactants and products of these reactions. But many public datasets include identified metabolites whose pathway involvement is unknown, hindering metabolic interpretation. To address these shortcomings, various machine learning models, including those trained on data from the Kyoto Encyclopedia of Genes and Genomes (KEGG), have been developed to predict the pathway involvement of metabolites based on their chemical descriptions; however, these prior models are based on old metabolite KEGG-based datasets, including one benchmark dataset that is invalid due to the presence of over 1500 duplicate entries. Therefore, we have developed a new benchmark dataset derived from the KEGG following optimal standards of scientific computational reproducibility and including all source code needed to update the benchmark dataset as KEGG changes. We have used this new benchmark dataset with our atom coloring methodology to develop and compare the performance of Random Forest, XGBoost, and multilayer perceptron with autoencoder models generated from our new benchmark dataset. Best overall weighted average performance across 1000 unique folds was an F1 score of 0.8180 and a Matthews correlation coefficient of 0.7933, which was provided by XGBoost binary classification models for 11 KEGG-defined pathway categories.

From network analysis to experimental validation: identification of regulators of non-muscle myosin II contractility using the folded-gastrulation signaling pathway

The morphogenetic process of apical constriction, which relies on non-muscle myosin II (NMII) generated constriction of apical domains of epithelial cells, is key to the development of complex cellular patterns. Apical constriction occurs in almost all multicellular organisms, but one of the most well-characterized systems is the Folded-gastrulation (Fog)-induced apical constriction that occurs in Drosophila. The binding of Fog to its cognizant receptors Mist/Smog results in a signaling cascade that leads to the activation of NMII-generated contractility. Despite our knowledge of key molecular players involved in Fog signaling, we sought to explore whether other proteins have an undiscovered role in its regulation. We developed a computational method to predict unidentified candidate NMII regulators using a network of pairwise protein-protein interactions called an interactome. We first constructed a Drosophila interactome of over 500,000 protein-protein interactions from several databases that curate high-throughput experiments. Next, we implemented several graph-based algorithms that predicted 14 proteins potentially involved in Fog signaling. To test these candidates, we used RNAi depletion in combination with a cellular contractility assay in Drosophila S2R + cells, which respond to Fog by contracting in a stereotypical manner. Of the candidates we screened using this assay, two proteins, the serine/threonine phosphatase Flapwing and the putative guanylate kinase CG11811 were demonstrated to inhibit cellular contractility when depleted, suggestive of their roles as novel regulators of the Fog pathway.

Benchmark dataset for training machine learning models to predict the pathway involvement of metabolites

Metabolic pathways are a human-defined grouping of life sustaining biochemical reactions, metabolites being both the reactants and products of these reactions. But many public datasets include identified metabolites whose pathway involvement is unknown, hindering metabolic interpretation. To address these shortcomings, various machine learning models, including those trained on data from the Kyoto Encyclopedia of Genes and Genomes (KEGG), have been developed to predict the pathway involvement of metabolites based on their chemical descriptions; however, these prior models are based on old metabolite KEGG-based datasets, including one benchmark dataset that is invalid due to the presence of over 1500 duplicate entries. Therefore, we have developed a new benchmark dataset derived from the KEGG following optimal standards of scientific computational reproducibility and including all source code needed to update the benchmark dataset as KEGG changes. We have used this new benchmark dataset with our atom coloring methodology to develop and compare the performance of Random Forest, XGBoost, and multilayer perceptron with autoencoder models generated from our new benchmark dataset. Best overall weighted average performance across 1000 unique folds was an F1-score of 0.8180 and Matthews correlation coefficient of 0.7933, which was provided by XGBoost binary classification models for 11 KEGG-defined pathway categories.

Epigenome-wide association study of diabetic chronic kidney disease progression in the Korean population: the KNOW-CKD study

Since the etiology of diabetic chronic kidney disease (CKD) is multifactorial, studies on DNA methylation for kidney function deterioration have rarely been performed despite the need for an epigenetic approach. Therefore, this study aimed to identify epigenetic markers associated with CKD progression based on the decline in the estimated glomerular filtration rate in diabetic CKD in Korea. An epigenome-wide association study was performed using whole blood samples from 180 CKD recruited from the KNOW-CKD cohort. Pyrosequencing was also performed on 133 CKD participants as an external replication analysis. Functional analyses, including the analysis of disease-gene networks, reactome pathways, and protein-protein interaction networks, were conducted to identify the biological mechanisms of CpG sites. A phenome-wide association study was performed to determine the associations between CpG sites and other phenotypes. Two epigenetic markers, cg10297223 on AGTR1 and cg02990553 on KRT28 indicated a potential association with diabetic CKD progression. Based on the functional analyses, other phenotypes (blood pressure and cardiac arrhythmia for AGTR1) and biological pathways (keratinization and cornified envelope for KRT28) related to CKD were also identified. This study suggests a potential association between the cg10297223 and cg02990553 and the progression of diabetic CKD in Koreans. Nevertheless, further validation is needed through additional studies.

A systematic review of biologically-informed deep learning models for cancer: fundamental trends for encoding and interpreting oncology data

BACKGROUND

There is an increasing interest in the use of Deep Learning (DL) based methods as a supporting analytical framework in oncology. However, most direct applications of DL will deliver models with limited transparency and explainability, which constrain their deployment in biomedical settings.

METHODS

This systematic review discusses DL models used to support inference in cancer biology with a particular emphasis on multi-omics analysis. It focuses on how existing models address the need for better dialogue with prior knowledge, biological plausibility and interpretability, fundamental properties in the biomedical domain. For this, we retrieved and analyzed 42 studies focusing on emerging architectural and methodological advances, the encoding of biological domain knowledge and the integration of explainability methods.

RESULTS

We discuss the recent evolutionary arch of DL models in the direction of integrating prior biological relational and network knowledge to support better generalisation (e.g. pathways or Protein-Protein-Interaction networks) and interpretability. This represents a fundamental functional shift towards models which can integrate mechanistic and statistical inference aspects. We introduce a concept of bio-centric interpretability and according to its taxonomy, we discuss representational methodologies for the integration of domain prior knowledge in such models.

CONCLUSIONS

The paper provides a critical outlook into contemporary methods for explainability and interpretability used in DL for cancer. The analysis points in the direction of a convergence between encoding prior knowledge and improved interpretability. We introduce bio-centric interpretability which is an important step towards formalisation of biological interpretability of DL models and developing methods that are less problem- or application-specific.

E3 ubiquitin ligase ASB8 promotes selinexor-induced proteasomal degradation of XPO1

Selinexor (KPT-330), a small-molecule inhibitor of exportin-1 (XPO1, CRM1) with potent anticancer activity, has recently been granted FDA approval for treatment of relapsed/refractory multiple myeloma and diffuse large B-cell lymphoma (DLBCL), with a number of additional indications currently under clinical investigation. Since selinexor has often demonstrated synergy when used in combination with other drugs, notably bortezomib and dexamethasone, a more comprehensive approach to uncover new beneficial interactions would be of great value. Moreover, stratifying patients, personalizing therapeutics and improving clinical outcomes requires a better understanding of the genetic vulnerabilities and resistance mechanisms underlying drug response. Here, we used CRISPR-Cas9 loss-of-function chemogenetic screening to identify drug-gene interactions with selinexor in chronic myeloid leukemia, multiple myeloma and DLBCL cell lines. We identified the TGFβ-SMAD4 pathway as an important mediator of resistance to selinexor in multiple myeloma cells. Moreover, higher activity of this pathway correlated with prolonged progression-free survival in multiple myeloma patients treated with selinexor, indicating that the TGFβ-SMAD4 pathway is a potential biomarker predictive of therapeutic outcome. In addition, we identified ASB8 (ankyrin repeat and SOCS box containing 8) as a shared modulator of selinexor sensitivity across all tested cancer types, with both ASB8 knockout and overexpression resulting in selinexor hypersensitivity. Mechanistically, we showed that ASB8 promotes selinexor-induced proteasomal degradation of XPO1. This study provides insight into the genetic factors that influence response to selinexor treatment and could support both the development of predictive biomarkers as well as new drug combinations.

ABL1/2 and DDR1 Drive MEKi Resistance in NRAS-Mutant Melanomas by Stabilizing RAF/MYC/ETS1 and Promoting RAF Homodimerization

Melanomas harboring NRAS mutations are a particularly aggressive and deadly subtype. If patients cannot tolerate or the melanomas are insensitive to immune checkpoint blockade, there are no effective 2nd-line treatment options. Drugs targeting the RAF/MEK/ERK pathway, which are used for BRAF-mutant melanomas, do little to increase progression-free survival (PFS). Here, using both loss-of-function and gain-of-function approaches, we show that ABL1/2 and DDR1 are critical nodes during NRAS-mutant melanoma intrinsic and acquired MEK inhibitor (MEKi) resistance. In some acquired resistance cells, ABL1/2 and DDR1 cooperate to stabilize RAF proteins, activate ERK cytoplasmic and nuclear signaling, repress p27/KIP1 expression, and drive RAF homodimerization. In contrast, other acquired resistance cells depend solely on ABL1/2 for their survival, and are sensitive to highly specific allosteric ABL1/2 inhibitors, which prevent β-catenin nuclear localization and destabilize MYC and ETS1 in an ERK-independent manner. Significantly, targeting ABL1/2 and DDR1 with an FDA-approved anti-leukemic drug, reverses intrinsic MEKi resistance, delays acquisition of acquired resistance, and doubles the survival time in a NRAS-mutant mouse model. These data indicate that repurposing FDA-approved drugs targeting ABL1/2 and DDR1 may be a novel and effective strategy for treating patients with treatment-refractory NRAS-driven melanomas.

A comprehensive survey of the approaches for pathway analysis using multi-omics data integration

Pathway analysis has been widely used to detect pathways and functions associated with complex disease phenotypes. The proliferation of this approach is due to better interpretability of its results and its higher statistical power compared with the gene-level statistics. A plethora of pathway analysis methods that utilize multi-omics setup, rather than just transcriptomics or proteomics, have recently been developed to discover novel pathways and biomarkers. Since multi-omics gives multiple views into the same problem, different approaches are employed in aggregating these views into a comprehensive biological context. As a result, a variety of novel hypotheses regarding disease ideation and treatment targets can be formulated. In this article, we review 32 such pathway analysis methods developed for multi-omics and multi-cohort data. We discuss their availability and implementation, assumptions, supported omics types and databases, pathway analysis techniques and integration strategies. A comprehensive assessment of each method's practicality, and a thorough discussion of the strengths and drawbacks of each technique will be provided. The main objective of this survey is to provide a thorough examination of existing methods to assist potential users and researchers in selecting suitable tools for their data and analysis purposes, while highlighting outstanding challenges in the field that remain to be addressed for future development.

Applying differential network analysis to longitudinal gene expression in response to perturbations

Differential Network (DN) analysis is a method that has long been used to interpret changes in gene expression data and provide biological insights. The method identifies the rewiring of gene networks in response to external perturbations. Our study applies the DN method to the analysis of RNA-sequencing (RNA-seq) time series datasets. We focus on expression changes: (i) in saliva of a human subject after pneumococcal vaccination (PPSV23) and (ii) in primary B cells treated ex vivo with a monoclonal antibody drug (Rituximab). The DN method enabled us to identify the activation of biological pathways consistent with the mechanisms of action of the PPSV23 vaccine and target pathways of Rituximab. The community detection algorithm on the DN revealed clusters of genes characterized by collective temporal behavior. All saliva and some B cell DN communities showed characteristic time signatures, outlining a chronological order in pathway activation in response to the perturbation. Moreover, we identified early and delayed responses within network modules in the saliva dataset and three temporal patterns in the B cell data.

A transcriptomic analysis of the effects of macrophage polarization and endotoxin tolerance on the response to Salmonella

Salmonella is an intracellular pathogen causing significant morbidity and mortality. Its ability to grow inside macrophages is important to virulence, and is dependent on the activation state of the macrophages. Classically activated M1 macrophages are non-permissive for Salmonella growth, while alternatively activated M2 macrophages are permissive for Salmonella growth. Here we showed that endotoxin-primed macrophages (MEP), such as those associated with sepsis, showed similar levels of Salmonella resistance to M1 macrophages after 2 hr of intracellular infection, but at the 4 hr and 24 hr time points were susceptible like M2 macrophages. To understand this mechanistically, transcriptomic sequencing, RNA-Seq, was performed. This showed that M1 and MEP macrophages that had not been exposed to Salmonella, demonstrated a process termed here as primed activation, in expressing relatively higher levels of particular anti-infective genes and pathways, including the JAK-STAT (Janus kinase-signal transducer and activator of transcription) pathway. In contrast, in M2 macrophages these genes and pathways were largely expressed only in response to infection. Conversely, in response to infection, M1 macrophages, but not MEP macrophages, modulated additional genes known to be associated with susceptibility to Salmonella infection, possibly contributing to the differences in resistance at later time points. Application of the JAK inhibitor Ruxolitinib before infection reduced resistance in M1 macrophages, supporting the importance of early JAK-STAT signalling in M1 resistance to Salmonella.

Exploring the multidimensional heterogeneities of glioblastoma multiforme based on sample-specific edge perturbation in gene interaction network

Glioblastoma multiforme (GBM) is the most malignant brain cancer with great heterogeneities in many aspects, such as prognosis, clinicopathological features, immune landscapes, and immunotherapeutic responses. Considering that gene interaction network is relatively stable in a healthy state but widely perturbed in cancers, we sought to explore the multidimensional heterogeneities of GBM through evaluating the degree of network perturbations. The gene interaction network perturbations of GBM samples (TCGA cohort) and normal samples (GTEx database) were characterized by edge perturbations, which were quantized through evaluating the change in relative gene expression value. An unsupervised consensus clustering analysis was performed to identify edge perturbation-based clusters of GBM samples. Results revealed that the edge perturbation of GBM samples was stronger than that of normal samples. Four edge perturbation-based clusters of GBM samples were identified and showed prominent heterogeneities in prognosis, clinicopathological features, somatic genomic alterations, immune landscapes, and immunotherapeutic responses. In addition, a sample-specific perturbation of gene interaction score (SPGIScore) was constructed based on the differently expressed genes (DEGs) among four clusters, and exhibited a robust ability to predict prognosis. In conclusion, the bioinformatics approach based on sample-specific edge perturbation in gene interaction network provided a new perspective to understanding the multidimensional heterogeneities of GBM.

Biolinguistics: A Scientometric Analysis of Research on (Children's) Molecular Genetics of Speech and Language (Disorders)

There are numerous children and adolescents throughout the world who are either diagnosed with speech and language disorders or manifest any of them as a result of another disorder. Meanwhile, since the emergence of language as an innate capability, the question of whether it constitutes a behaviour or an innate ability has been debated for decades. There have been several theories developed that support and demonstrate the biological foundations of human language. Molecular evidence of the biological basis of language came from the FOXP2 gene, also known as the language gene. Taking a closer look at both human language and biology, biolinguistics is at the core of these inquiries-attempting to understand the aetiologies of the genetics of speech and language disorders in children and adolescents. This paper presents empirical evidence based on both scientometrics and bibliometrics. We collected data between 1935 and 2022 from Scopus, WOS, and Lens. A total of 1570 documents were analysed from Scopus, 1440 from the WOS, and 5275 from Lens. Bibliometric analysis was performed using Excel based on generated reports from these three databases. CiteSpace 5.8.R3 and VOSviewer 1.6.18 were used to conduct the scientometric analysis. Eight bibliometric and eight scientometric indicators were used to measure the development of the field of biolinguistics, including but not limited to the production size of knowledge, the most examined topics, and the most frequent concepts and variables. A major finding of our study is identifying the most examined topics in the genetics of speech and language disorders. These included: gestural communication, structural design, cultural evolution, neural network, language tools, human language faculty, evolutionary biology, molecular biology, and theoretical perspective on language evolution.

COVID-19-specific transcriptomic signature detectable in blood across multiple cohorts

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading across the world despite vast global vaccination efforts. Consequently, many studies have looked for potential human host factors and immune mechanisms associated with the disease. However, most studies have focused on comparing COVID-19 patients to healthy controls, while fewer have elucidated the specific host factors distinguishing COVID-19 from other infections. To discover genes specifically related to COVID-19, we reanalyzed transcriptome data from nine independent cohort studies, covering multiple infections, including COVID-19, influenza, seasonal coronaviruses, and bacterial pneumonia. The identified COVID-19-specific signature consisted of 149 genes, involving many signals previously associated with the disease, such as induction of a strong immunoglobulin response and hemostasis, as well as dysregulation of cell cycle-related processes. Additionally, potential new gene candidates related to COVID-19 were discovered. To facilitate exploration of the signature with respect to disease severity, disease progression, and different cell types, we also offer an online tool for easy visualization of the selected genes across multiple datasets at both bulk and single-cell levels.

Application of modular response analysis to medium- to large-size biological systems

The development of high-throughput genomic technologies associated with recent genetic perturbation techniques such as short hairpin RNA (shRNA), gene trapping, or gene editing (CRISPR/Cas9) has made it possible to obtain large perturbation data sets. These data sets are invaluable sources of information regarding the function of genes, and they offer unique opportunities to reverse engineer gene regulatory networks in specific cell types. Modular response analysis (MRA) is a well-accepted mathematical modeling method that is precisely aimed at such network inference tasks, but its use has been limited to rather small biological systems so far. In this study, we show that MRA can be employed on large systems with almost 1,000 network components. In particular, we show that MRA performance surpasses general-purpose mutual information-based algorithms. Part of these competitive results was obtained by the application of a novel heuristic that pruned MRA-inferred interactions a posteriori. We also exploited a block structure in MRA linear algebra to parallelize large system resolutions.

Evaluating the predictive accuracy of curated biological pathways in a public knowledgebase

Abstract

Reactome is a database of human biological pathways manually curated from the primary literature and peer-reviewed by experts. To evaluate the utility of Reactome pathways for predicting functional consequences of genetic perturbations, we compared predictions of perturbation effects based on Reactome pathways against published empirical observations. Ten cancer-relevant Reactome pathways, representing diverse biological processes such as signal transduction, cell division, DNA repair and transcriptional regulation, were selected for testing. For each pathway, root input nodes and key pathway outputs were defined. We then used pathway-diagram-derived logic graphs to predict, either by inspection by biocurators or using a novel algorithm MP-BioPath, the effects of bidirectional perturbations (upregulation/activation or downregulation/inhibition) of single root inputs on the status of key outputs. These predictions were then compared to published empirical tests. In total, 4968 test cases were analyzed across 10 pathways, of which 847 were supported by published empirical findings. Out of the 847 test cases, curators’ predictions agreed with the experimental evidence in 670 and disagreed in 177 cases, resulting in ∼81% overall accuracy. MP-BioPath predictions agreed with experimental evidence for 625 and disagreed for 222 test cases, resulting in ∼75% overall accuracy. The expected accuracy of random guessing was 33%. Per-pathway accuracy did not correlate with the number of pathway edges nor the number of pathway nodes but varied across pathways, ranging from 56% (curator)/44% (MP-BioPath) for ‘Mitotic G1 phase and G1/S transition’ to 100% (curator)/94% (MP-BioPath) for ‘RAF/MAP kinase cascade’. This study highlights the potential of pathway databases such as Reactome in modeling genetic perturbations, promoting standardization of experimental pathway activity readout and supporting hypothesis-driven research by revealing relationships between pathway inputs and outputs that have not yet been directly experimentally tested.

Database URL

www.reactome.org

HyMM: hybrid method for disease-gene prediction by integrating multiscale module structure

MOTIVATION

Identifying disease-related genes is an important issue in computational biology. Module structure widely exists in biomolecule networks, and complex diseases are usually thought to be caused by perturbations of local neighborhoods in the networks, which can provide useful insights for the study of disease-related genes. However, the mining and effective utilization of the module structure is still challenging in such issues as a disease gene prediction.

RESULTS

We propose a hybrid disease-gene prediction method integrating multiscale module structure (HyMM), which can utilize multiscale information from local to global structure to more effectively predict disease-related genes. HyMM extracts module partitions from local to global scales by multiscale modularity optimization with exponential sampling, and estimates the disease relatedness of genes in partitions by the abundance of disease-related genes within modules. Then, a probabilistic model for integration of gene rankings is designed in order to integrate multiple predictions derived from multiscale module partitions and network propagation, and a parameter estimation strategy based on functional information is proposed to further enhance HyMM's predictive power. By a series of experiments, we reveal the importance of module partitions at different scales, and verify the stable and good performance of HyMM compared with eight other state-of-the-arts and its further performance improvement derived from the parameter estimation.

CONCLUSIONS

The results confirm that HyMM is an effective framework for integrating multiscale module structure to enhance the ability to predict disease-related genes, which may provide useful insights for the study of the multiscale module structure and its application in such issues as a disease-gene prediction.

The importance of enhancer methylation for epigenetic regulation of tumorigenesis in squamous lung cancer

Lung squamous cell carcinoma (LUSC) is a subtype of non-small cell lung cancer (NSCLC). LUSC occurs at the bronchi, shows a squamous appearance, and often occurs in smokers. To determine the epigenetic regulatory mechanisms of tumorigenesis, we performed a genome-wide analysis of DNA methylation in tumor and adjacent normal tissues from LUSC patients. With the Infinium Methylation EPIC Array, > 850,000 CpG sites, including ~350,000 CpG sites for enhancer regions, were profiled, and the differentially methylated regions (DMRs) overlapping promoters (pDMRs) and enhancers (eDMRs) between tumor and normal tissues were identified. Dimension reduction based on DMR profiles revealed that eDMRs alone and not pDMRs alone can differentiate tumors from normal tissues with the equivalent performance of total DMRs. We observed a stronger negative correlation of LUSC-specific gene expression with methylation for enhancers than promoters. Target genes of eDMRs rather than pDMRs were found to be enriched for tumor-associated genes and pathways. Furthermore, DMR methylation associated with immune infiltration was more frequently observed among enhancers than promoters. Our results suggest that methylation of enhancer regions rather than promoters play more important roles in epigenetic regulation of tumorigenesis and immune infiltration in LUSC.

Identifying individualized risk subpathways reveals pan-cancer molecular classification based on multi-omics data

Cancer is a highly heterogeneous disease with different functional disorders among individuals. The initiation and progression of cancer is usually related to dysregulation of local regions within pathways. Identification of individualized risk pathways is crucial for revealing the mechanisms of tumorigenesis and heterogeneity. However, approach that focused on mining patient-specific risk subpathway regions is still lacking. Here, we developed an individualized cancer risk subpathway identification method that was referred as InCRiS by integrating multi-omics data. Then, the method was applied to nearly 3000 samples across 9 TCGA cancer types and its robustness and reliability were evaluated. Dissecting dysregulated subpathways in these tumor samples revealed several key regions which may play oncogenic roles in cancer. The construction of risk subpathway dysregulation profile of pan-cancers revealed 11 pan-cancer molecular classification (InCRiS subtypes) with significantly different clinical outcomes. Moreover, subpathway regions that tend to be disordered in individuals of specific subtypes were examined for understanding the pathogenesis of tumor and some key genes such as CTNNB1, EP300 and PRKDC were nominated in different subtypes. In summary, the proposed method and resulting data presented useful resources to study the mechanism of tumor heterogeneity and to discovery novel therapeutic targets for precise treatment of cancer.

Bioinformatic analyses of hydroxylated polybrominated diphenyl ethers toxicities on impairment of adrenocortical secretory function

BACKGROUND

Polybrominated diphenyl ethers (PBDEs) and their metabolites have severe impact on human health, but few studies focus on their nephrotoxicity. This study was conceived to explore hub genes that may be involved in two hydroxylated polybrominated diphenyl ethers toxicities on impairment of adrenocortical secretory function.

METHODS

Gene dataset was obtained from Gene Expression Omnibus (GEO). Principal component analysis and correlation analysis were used to classify the samples. Differentially expressed genes (DEGs) were screened using the limma package in RStudio (version 4.1.0). Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome enrichment analyses of DEGs were conducted. Protein-protein interaction (PPI) network was established using STRING network, and genes were filtered by Cytoscape (version 3.8.2). Finally, the hub genes were integrated by plug-in CytoHubba and RobustRankAggreg, and were preliminarily verified by the Comparative Toxicogenomics Database (CTD).

RESULTS

GSE8588 dataset was selected in this study. About 190 upregulated and 224 downregulated DEGs in 2-OH-BDE47 group, and 244 upregulated and 276 downregulated DEGs in 2-OH-BDE85 group. Functional enrichment analyses in the GO, KEGG and Reactome indicated the potential involvement of DEGs in endocrine metabolism, oxidative stress mechanisms, regulation of abnormal cell proliferation, apoptosis, DNA damage and repair. 2-OH-BDE85 is more cytotoxic in a dose-dependent manner than 2-OH-BDE47. A total of 98 hub genes were filtered, and 91 nodes and 359 edges composed the PPI network. Besides, 9 direct-acting genes were filtered for the intersection of hub genes by CTD.

CONCLUSIONS

OH-PBDEs may induce H295R adrenocortical cancer cells in the disorder of endocrine metabolism, regulation of abnormal cell proliferation, DNA damage and repair. The screened hub genes may play an important role in this dysfunction.

RDFizing the biosynthetic pathway of E.coli O-antigen to enable semantic sharing of microbiology data

Background

The abundance of glycomics data that have accumulated has led to the development of many useful databases to aid in the understanding of the function of the glycans and their impact on cellular activity. At the same time, the endeavor for data sharing between glycomics databases with other biological databases have contributed to the creation of new knowledgebases. However, different data types in data description have impeded the data sharing for knowledge integration. To solve this matter, Semantic Web techniques including Resource Description Framework (RDF) and ontology development have been adopted by various groups to standardize the format for data exchange. These semantic data have contributed to the expansion of knowledgebases and hold promises of providing data that can be intelligently processed. On the other hand, bench biologists who are experts in experimental finding are end users and data producers. Therefore, it is indispensable to reduce the technical barrier required for bench biologists to manipulate their experimental data to be compatible with standard formats for data sharing.

Results

There are many essential concepts and practical techniques for data integration but there is no method to enable researchers to easily apply Semantic Web techniques to their experimental data. We implemented our procedure on unformatted information of E.coli O-antigen structures collected from the web and show how this information can be expressed as formatted data applicable to Semantic Web standards. In particular, we described the E-coli O-antigen biosynthesis pathway using the BioPAX ontology developed to support data exchange between pathway databases.

Conclusions

The method we implemented to semantically describe O-antigen biosynthesis should be helpful for biologists to understand how glycan information, including relevant pathway reaction data, can be easily shared. We hope this method can contribute to lower the technical barrier that is required when experimental findings are formulated into formal representations and can lead bench scientists to readily participate in the construction of new knowledgebases that are integrated with existing ones. Such integration over the Semantic Web will enable future work in artificial intelligence and machine learning to enable computers to infer new relationships and hypotheses in the life sciences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02384-y.

A systematic genome-wide mapping of oncogenic mutation selection during CRISPR-Cas9 genome editing

Recent studies have reported that genome editing by CRISPR–Cas9 induces a DNA damage response mediated by p53 in primary cells hampering their growth. This could lead to a selection of cells with pre-existing p53 mutations. In this study, employing an integrated computational and experimental framework, we systematically investigated the possibility of selection of additional cancer driver mutations during CRISPR-Cas9 gene editing. We first confirm the previous findings of the selection for pre-existing p53 mutations by CRISPR-Cas9. We next demonstrate that similar to p53, wildtype KRAS may also hamper the growth of Cas9-edited cells, potentially conferring a selective advantage to pre-existing KRAS-mutant cells. These selective effects are widespread, extending across cell-types and methods of CRISPR-Cas9 delivery and the strength of selection depends on the sgRNA sequence and the gene being edited. The selection for pre-existing p53 or KRAS mutations may confound CRISPR-Cas9 screens in cancer cells and more importantly, calls for monitoring patients undergoing CRISPR-Cas9-based editing for clinical therapeutics for pre-existing p53 and KRAS mutations.

CRISPR-Cas9 gene editing can induce a p53 mediated damage response. Here the authors investigate the possibility of selection of pre-existing cancer driver mutations during CRISPR-Cas9 knockout based gene editing and identify KRAS mutants that may confer a selected advantage to edited cells.

DysPIA: A Novel Dysregulated Pathway Identification Analysis Method

Differential co-expression-based pathway analysis is still limited and not widely used. In most current methods, the pathways were considered as gene sets, but the gene regulation relationships were not considered, and the computational speed was slow. In this article, we proposed a novel Dysregulated Pathway Identification Analysis (DysPIA) method to overcome these shortcomings. We adopted the idea of Correlation by Individual Level Product into analysis and performed a fast enrichment analysis. We constructed a combined gene-pair background which was much more sufficient than the background used in Edge Set Enrichment Analysis. In simulation study, DysPIA was able to identify the causal pathways with high AUC (0.9584 to 0.9896). In p53 mutation data, DysPIA obtained better performance than other methods. It obtained more potential dysregulated pathways that could be literature verified, and it ran much faster (∼1,700-8,000 times faster than other methods when 10,000 permutations). DysPIA was also applied to breast cancer relapse dataset and breast cancer subtype dataset. The results show that DysPIA is effective and has a great biological significance. R packages "DysPIA" and "DysPIAData" are constructed and freely available on R CRAN (https://cran.r-project.org/web/packages/DysPIA/index.html and https://cran.r-project.org/web/packages/DysPIAData/index.html), and on GitHub (https://github.com/lemonwang2020).

Transcriptome analysis of MBD5-associated neurodevelopmental disorder (MAND) neural progenitor cells reveals dysregulation of autism-associated genes

MBD5-associated neurodevelopmental disorder (MAND) is an autism spectrum disorder (ASD) characterized by intellectual disability, motor delay, speech impairment and behavioral problems; however, the biological role of methyl-CpG-binding domain 5, MBD5, in neurodevelopment and ASD remains largely undefined. Hence, we created neural progenitor cells (NPC) derived from individuals with chromosome 2q23.1 deletion and conducted RNA-seq to identify differentially expressed genes (DEGs) and the biological processes and pathways altered in MAND. Primary skin fibroblasts from three unrelated individuals with MAND and four unrelated controls were converted into induced pluripotent stem cell (iPSC) lines, followed by directed differentiation of iPSC to NPC. Transcriptome analysis of MAND NPC revealed 468 DEGs (q < 0.05), including 20 ASD-associated genes. Comparison of DEGs in MAND with SFARI syndromic autism genes revealed a striking significant overlap in biological processes commonly altered in neurodevelopmental phenotypes, with TGFβ, Hippo signaling, DNA replication, and cell cycle among the top enriched pathways. Overall, these transcriptome deviations provide potential connections to the overlapping neurocognitive and neuropsychiatric phenotypes associated with key high-risk ASD genes, including chromatin modifiers and epigenetic modulators, that play significant roles in these disease states.

Spatial Distribution of Private Gene Mutations in Clear Cell Renal Cell Carcinoma

Simple Summary

Tumours consist of multiple groups of similar cells resulting from differing evolutionary trajectories, i.e., subclones. These subclones are prevalent in clear cell renal cell carcinoma (ccRCC). The aim of this study is to determine how similar or dissimilar the subclones in 89 ccRCC tumours are from one another regarding their gene mutations and expression profiles, i.e., the extent of intra-tumour heterogeneity. The implications of these alterations with respect to signalling pathways is also assessed. Deep sequencing allows for the identification of mutations with low-allele frequencies, providing a more comprehensive view of the heterogeneity present in the tumours. With an average of 62% of mutations having been identified in only one of the two biopsies, some of which in turn are found to impact gene expression, the complex makeup of ccRCC tumours is evident, and this can drastically influence treatment outcome.

Abstract

Intra-tumour heterogeneity is the molecular hallmark of renal cancer, and the molecular tumour composition determines the treatment outcome of renal cancer patients. In renal cancer tumourigenesis, in general, different tumour clones evolve over time. We analysed intra-tumour heterogeneity and subclonal mutation patterns in 178 tumour samples obtained from 89 clear cell renal cell carcinoma patients. In an initial discovery phase, whole-exome and transcriptome sequencing data from paired tumour biopsies from 16 ccRCC patients were used to design a gene panel for follow-up analysis. In this second phase, 826 selected genes were targeted at deep coverage in an extended cohort of 89 patients for a detailed analysis of tumour heterogeneity. On average, we found 22 mutations per patient. Pairwise comparison of the two biopsies from the same tumour revealed that on average, 62% of the mutations in a patient were detected in one of the two samples. In addition to commonly mutated genes (VHL, PBRM1, SETD2 and BAP1), frequent subclonal mutations with low variant allele frequency (<10%) were observed in TP53 and in mucin coding genes MUC6, MUC16, and MUC3A. Of the 89 ccRCC tumours, 87 (~98%) harboured private mutations, occurring in only one of the paired tumour samples. Clonally exclusive pathway pairs were identified using the WES data set from 16 ccRCC patients. Our findings imply that shared and private mutations significantly contribute to the complexity of differential gene expression and pathway interaction and might explain the clonal evolution of different molecular renal cancer subgroups. Multi-regional sequencing is central for the identification of subclones within ccRCC.

High-content phenotypic and pathway profiling to advance drug discovery in diseases of unmet need

Conventional thinking in modern drug discovery postulates that the design of highly selective molecules which act on a single disease-associated target will yield safer and more effective drugs. However, high clinical attrition rates and the lack of progress in developing new effective treatments for many important diseases of unmet therapeutic need challenge this hypothesis. This assumption also impinges upon the efficiency of target agnostic phenotypic drug discovery strategies, where early target deconvolution is seen as a critical step to progress phenotypic hits. In this review we provide an overview of how emerging phenotypic and pathway-profiling technologies integrate to deconvolute the mechanism-of-action of phenotypic hits. We propose that such in-depth mechanistic profiling may support more efficient phenotypic drug discovery strategies that are designed to more appropriately address complex heterogeneous diseases of unmet need.

Lead DEAD/H box helicase biomarkers with the therapeutic potential identified by integrated bioinformatic approaches in lung cancer

DEAD/H box helicases are implicated in lung cancer but have not been systematically investigated for their clinical significance and function. In this study, we aimed to evaluate the potential of DEAD/H box helicases as prognostic biomarkers and therapeutic targets in lung cancer by integrated bioinformatic analysis of multivariate large-scale databases. Survival and differential expression analysis of these helicases enabled us to identify four biomarkers with the most significant alterations. These were found to be the negative prognostic factors DDX11, DDX55 and DDX56, and positive prognostic factor DDX5. Pathway enrichment analysis indicates that MYC signalling is negatively associated with expression levels of the DDX5 gene while positively associated with that of DDX11, DDX55 and DDX56. High expression levels of the DDX5 gene is associated with low mutation levels of TP53 and MUC16, the two most frequently mutated genes in lung cancer. In contrast, high expression levels of DDX11, DDX55 and DDX56 genes are associated with high levels of TP53 and MUC16 mutation. The tumour-infiltrated CD8 + T and B cells positively correlate with levels of DDX5 gene expression, while negatively correlate with that of the other three DEAD box helicases, respectively. Moreover, the DDX5-associated miRNA profile is distinguished from the miRNA profiles of DDX11, DDX55 and DDX56, although each DDX has a different miRNA signature. The identification of these four DDX helicases as biomarkers will be valuable for prognostic prediction and targeted therapeutic development in lung cancer.

Frailty subtypes and recovery in older survivors of acute respiratory failure: a pilot study

BACKGROUND

Identifying subtypes of acute respiratory failure survivors may facilitate patient selection for post-intensive care unit (ICU) follow-up clinics and trials.

METHODS

We conducted a single-centre prospective cohort study of 185 acute respiratory failure survivors, aged ≥ 65 years. We applied latent class modelling to identify frailty subtypes using frailty phenotype and cognitive impairment measurements made during the week before hospital discharge. We used Fine-Gray competing risks survival regression to test associations between frailty subtypes and recovery, defined as returning to a basic Activities of Daily Living disability count less than or equal to the pre-hospitalisation count within 6 months. We characterised subtypes by pre-ICU frailty (Clinical Frailty Scale score ≥ 5), the post-ICU frailty phenotype, and serum inflammatory cytokines, hormones and exosome proteomics during the week before hospital discharge.

RESULTS

We identified five frailty subtypes. The recovery rate decreased 49% across each subtype independent of age, sex, pre-existing disability, comorbidity and Acute Physiology and Chronic Health Evaluation II score (recovery rate ratio: 0.51, 95% CI 0.41 to 0.63). Post-ICU frailty phenotype prevalence increased across subtypes, but pre-ICU frailty prevalence did not. In the subtype with the slowest recovery, all had cognitive impairment. The three subtypes with the slowest recovery had higher interleukin-6 levels (p=0.03) and a higher prevalence of ≥ 2 deficiencies in insulin growth factor-1, dehydroepiandrostersone-sulfate, or free-testosterone (p=0.02). Exosome proteomics revealed impaired innate immunity in subtypes with slower recovery.

CONCLUSIONS

Frailty subtypes varied by prehospitalisation frailty and cognitive impairment at hospital discharge. Subtypes with the slowest recovery were similarly characterised by greater systemic inflammation and more anabolic hormone deficiencies at hospital discharge.

multiGSEA: a GSEA-based pathway enrichment analysis for multi-omics data

Background

Gaining biological insights into molecular responses to treatments or diseases from omics data can be accomplished by gene set or pathway enrichment methods. A plethora of different tools and algorithms have been developed so far. Among those, the gene set enrichment analysis (GSEA) proved to control both type I and II errors well. In recent years the call for a combined analysis of multiple omics layers became prominent, giving rise to a few multi-omics enrichment tools. Each of these has its own drawbacks and restrictions regarding its universal application.

Results

Here, we present the multiGSEA package aiding to calculate a combined GSEA-based pathway enrichment on multiple omics layers. The package queries 8 different pathway databases and relies on the robust GSEA algorithm for a single-omics enrichment analysis. In a final step, those scores will be combined to create a robust composite multi-omics pathway enrichment measure. multiGSEA supports 11 different organisms and includes a comprehensive mapping of transcripts, proteins, and metabolite IDs.

Conclusions

With multiGSEA we introduce a highly versatile tool for multi-omics pathway integration that minimizes previous restrictions in terms of omics layer selection, pathway database availability, organism selection and the mapping of omics feature identifiers. multiGSEA is publicly available under the GPL-3 license at https://github.com/yigbt/multiGSEA and at bioconductor: https://bioconductor.org/packages/multiGSEA.

Transcriptomic profile of Pea3 family members reveal regulatory codes for axon outgrowth and neuronal connection specificity

PEA3 transcription factor subfamily is present in a variety of tissues with branching morphogenesis, and play a particularly significant role in neural circuit formation and specificity. Many target genes in axon guidance and cell-cell adhesion pathways have been identified for Pea3 transcription factor (but not for Erm or Er81); however it was not so far clear whether all Pea3 subfamily members regulate same target genes, or whether there are unique targets for each subfamily member that help explain the exclusivity and specificity of these proteins in neuronal circuit formation. In this study, using transcriptomics and qPCR analyses in SH-SY5Y neuroblastoma cells, hypothalamic and hippocampal cell line, we have identified cell type-specific and subfamily member-specific targets for PEA3 transcription factor subfamily. While Pea3 upregulates transcription of Sema3D and represses Sema5B, for example, Erm and Er81 upregulate Sema5A and Er81 regulates Unc5C and Sema4G while repressing EFNB3 in SH-SY5Y neuroblastoma cells. We furthermore present a molecular model of how unique sites within the ETS domain of each family member can help recognize specific target motifs. Such cell-context and member-specific combinatorial expression profiles help identify cell-cell and cell-extracellular matrix communication networks and how they establish specific connections.

Genome-wide analysis of therapeutic response uncovers molecular pathways governing tamoxifen resistance in ER+ breast cancer

Background

Prioritization of breast cancer patients based on the risk of resistance to tamoxifen plays a significant role in personalized therapeutic planning and improving disease course and outcomes.

Methods

In this work, we demonstrate that a genome-wide pathway-centric computational framework elucidates molecular pathways as markers of tamoxifen resistance in ER+ breast cancer patients. In particular, we associated activity levels of molecular pathways with a wide spectrum of response to tamoxifen, which defined markers of tamoxifen resistance in patients with ER+ breast cancer.

Findings

We identified five biological pathways as markers of tamoxifen failure and demonstrated their ability to predict the risk of tamoxifen resistance in two independent patient cohorts (Test cohort1: log-rank p-value = 0.02, adjusted HR = 3.11; Test cohort2: log-rank p-value = 0.01, adjusted HR = 4.24). We have shown that these pathways are not markers of aggressiveness and outperform known markers of tamoxifen response. Furthermore, for adoption into clinic, we derived a list of pathway read-out genes and their associated scoring system, which assigns a risk of tamoxifen resistance for new incoming patients.

Interpretation

We propose that the identified pathways and their read-out genes can be utilized to prioritize patients who would benefit from tamoxifen treatment and patients at risk of tamoxifen resistance that should be offered alternative regimens.

Funding

This work was supported by the Rutgers SHP Dean's research grant, Rutgers start-up funds, Libyan Ministry of Higher Education and Scientific Research, and Katrina Kehlet Graduate Award from The NJ Chapter of the Healthcare Information Management Systems Society.

A systematic study of critical miRNAs on cells proliferation and apoptosis by the shortest path

BACKGROUND

MicroRNAs are a class of important small noncoding RNAs, which have been reported to be involved in the processes of tumorigenesis and development by targeting a few genes. Existing studies show that the imbalance between cell proliferation and apoptosis is closely related to the initiation and development of cancers. However, the impact of miRNAs on this imbalance has not been studied systematically.

RESULTS

In this study, we first construct a cell fate miRNA-gene regulatory network. Then, we propose a systematical method for calculating the global impact of miRNAs on cell fate genes based on the shortest path. Results on breast cancer and liver cancer datasets show that most of the cell fate genes are perturbed by the differentially expressed miRNAs. Most of the top-identified miRNAs are verified in the Human MicroRNA Disease Database (HMDD) and are related to breast and liver cancers. Function analysis shows that the top 20 miRNAs regulate multiple cell fate related function modules and interact tightly based on their functional similarity. Furthermore, more than half of them can promote sensitivity or induce resistance to some anti-cancer drugs. Besides, survival analysis demonstrates that the top-ranked miRNAs are significantly related to the overall survival time in the breast and liver cancers group.

CONCLUSION

In sum, this study can help to systematically study the important role of miRNAs on proliferation and apoptosis and thereby uncover the key miRNAs during the process of tumorigenesis. Furthermore, the results of this study will contribute to the development of clinical therapy based miRNAs for cancers.

Transcriptomic study of lipopolysaccharide-induced sepsis damage in a mouse heart model

Sepsis is an emergency systemic illness caused by pathogen infection and the combined result of the underactivity and overactivity of a patient's own immune system. However, the molecular mechanism of this illness remains largely unknown. Lipopolysaccharide (LPS) was injected to establish a sepsis model, and heart tissue was used to analyze transcriptome changes in mice. LPS injection was used to develop a sepsis model, which resulted in cardiac tissue rearrangement and inflammatory response activation. An RNA-sequencing-based transcriptome assay using mouse heart tissue with or without LPS injection showed that 3,326 and 1,769 genes were upregulated and downregulated, respectively (>2-fold changes; P<0.05). Furthermore, these differentially expressed genes were classified into 20 pathways, including ‘Wnt signaling pathway’, ‘VEGF signaling pathway’ and ‘TGF-β signaling pathway’, and these altered genes were enriched in 41 Gene Ontology terms. The application of Wnt3a inhibited the activation of the LPS-induced inflammatory response and activated Wnt signaling, as well as protecting against LPS-mediated cardiac tissue damage in mice. In contrast, inhibition of Wnt signaling by injection of its inhibitor IWR induced plasminogen activator inhibitor-1 expression and resulted in cardiac structure derangement, which was similar to the symptoms caused by injection of LPS, suggesting that LPS-induced damage to heart tissue may be via inhibition of Wnt signaling. The present analyses showed that Wnt signaling serves a pivotal role in sepsis development and may improve our understanding of the molecular basis underlying sepsis.

Transcriptome network of the papillary thyroid carcinoma radiation marker CLIP2

Background

We present a functional gene association network of the CLIP2 gene, generated by de-novo reconstruction from transcriptomic microarray data. CLIP2 was previously identified as a potential marker for radiation induced papillary thyroid carcinoma (PTC) of young patients in the aftermath of the Chernobyl reactor accident. Considering the rising thyroid cancer incidence rates in western societies, potentially related to medical radiation exposure, the functional characterization of CLIP2 is of relevance and contributes to the knowledge about radiation-induced thyroid malignancies.

Methods

We generated a transcriptomic mRNA expression data set from a CLIP2-perturbed thyroid cancer cell line (TPC-1) with induced CLIP2 mRNA overexpression and siRNA knockdown, respectively, followed by gene-association network reconstruction using the partial correlation-based approach GeneNet. Furthermore, we investigated different approaches for prioritizing differentially expressed genes for network reconstruction and compared the resulting networks with existing functional interaction networks from the Reactome, Biogrid and STRING databases. The derived CLIP2 interaction partners were validated on transcript and protein level.

Results

The best reconstructed network with regard to selection parameters contained a set of 20 genes in the 1st neighborhood of CLIP2 and suggests involvement of CLIP2 in the biological processes DNA repair/maintenance, chromosomal instability, promotion of proliferation and metastasis. Peptidylprolyl Isomerase Like 3 (PPIL3), previously identified as a potential direct interaction partner of CLIP2, was confirmed in this study by co-expression at the transcript and protein level.

Conclusion

In our study we present an optimized preselection approach for genes subjected to gene-association network reconstruction, which was applied to CLIP2 perturbation transcriptome data of a thyroid cancer cell culture model. Our data support the potential carcinogenic role of CLIP2 overexpression in radiation-induced PTC and further suggest potential interaction partners of the gene.

The Role of lncRNAs TAPIR-1 and -2 as Diagnostic Markers and Potential Therapeutic Targets in Prostate Cancer

In search of new biomarkers suitable for the diagnosis and treatment of prostate cancer, genome-wide transcriptome sequencing was carried out with tissue specimens from 40 prostate cancer (PCa) and 8 benign prostate hyperplasia patients. We identified two intergenic long non-coding transcripts, located in close genomic proximity, which are highly expressed in PCa. Microarray studies on a larger cohort comprising 155 patients showed a profound diagnostic potential of these transcripts (AUC~0.94), which we designated as tumor associated prostate cancer increased lncRNA (TAPIR-1 and -2). To test their therapeutic potential, knockdown experiments with siRNA were carried out. The knockdown caused an increase in the p53/TP53 tumor suppressor protein level followed by downregulation of a large number of cell cycle- and DNA-damage repair key regulators. Furthermore, in radiation therapy resistant tumor cells, the knockdown leads to a renewed sensitization of these cells to radiation treatment. Accordingly, in a preclinical PCa xenograft model in mice, the systemic application of nanoparticles loaded with siRNA targeting TAPIR-1 significantly reduced tumor growth. These findings point to a crucial role of TAPIR-1 and -2 in PCa.

Omics Integration for Mitochondria Systems Biology

Significance: Elucidation of the central importance of mitophagy in homeostasis of cells and organisms emphasizes that mitochondrial functions extend far beyond short-term needs for energy production. In mitochondria systems biology, the mitochondrial genome, proteome, and metabolome operate as a functional network in coordination of cell activities. Organization occurs through subnetworks that are interconnected by membrane potential, transport activities, allosteric and cooperative interactions, redox signaling mechanisms, rheostatic control by post-translational modifications, and metal ion homeostasis. These subnetworks enable use of varied energy precursors, defense against environmental stressors, and macromolecular rewiring to titrate energy production, biosynthesis, and detoxification according to cell-specific needs. Rewiring mechanisms, termed mitochondrial reprogramming, enhance fitness to respond to metabolic resources and challenges from the environment. Maladaptive responses can cause cell death. Maladaptive rewiring can cause disease. In cancer, adaptive rewiring can interfere with effective treatment. Recent Advances: Many recent advances have been facilitated by the development of new omics tools, which create opportunities to use data-driven analysis of omics data to address these complex adaptive and maladaptive mechanisms of mitochondrial reprogramming in human disease. Critical Issues: Application of omics integration to model systems reveals a critical role for metal ion homeostasis broadly impacting mitochondrial reprogramming. Importantly, data show that trans-omics associations are more robust and biologically relevant than single omics associations. Future Directions: Application of omics integration to mitophagy research creates new opportunities to link the complex, interactive functions of mitochondrial form and function in mitochondria systems biology.

Simultaneous kinase inhibition with ibrutinib and BCL2 inhibition with venetoclax offers a therapeutic strategy for acute myeloid leukemia

Acute myeloid leukemia (AML) results from the enhanced proliferation and impaired differentiation of hematopoietic stem and progenitor cells. Using an ex vivo functional screening assay, we identified that the combination of the BTK inhibitor ibrutinib and BCL2 inhibitor venetoclax (IBR + VEN), currently in clinical trials for chronic lymphocytic leukemia (CLL), demonstrated enhanced efficacy on primary AML patient specimens, AML cell lines, and in a mouse xenograft model of AML. Expanded analyses among a large cohort of hematologic malignancies (n = 651 patients) revealed that IBR + VEN sensitivity associated with selected genetic and phenotypic features in both CLL and AML specimens. Among AML samples, 11q23 MLL rearrangements were highly sensitive to IBR + VEN. Analysis of differentially expressed genes with respect to IBR + VEN sensitivity indicated pathways preferentially enriched in patient samples with reduced ex vivo sensitivity, including IL-10 signaling. These findings suggest that IBR + VEN may represent an effective therapeutic option for patients with AML.

Similarities and Differences of Blood N-Glycoproteins in Five Solid Carcinomas at Localized Clinical Stage Analyzed by SWATH-MS

Cancer is mostly incurable when diagnosed at a metastatic stage, making its early detection via blood proteins of immense clinical interest. Proteomic changes in tumor tissue may lead to changes detectable in the protein composition of circulating blood plasma. Using a proteomic workflow combining N-glycosite enrichment and SWATH mass spectrometry, we generate a data resource of 284 blood samples derived from patients with different types of localized-stage carcinomas and from matched controls. We observe whether the changes in the patient's plasma are specific to a particular carcinoma or represent a generic signature of proteins modified uniformly in a common, systemic response to many cancers. A quantitative comparison of the resulting N-glycosite profiles discovers that proteins related to blood platelets are common to several cancers (e.g., THBS1), whereas others are highly cancer-type specific. Available proteomics data, including a SWATH library to study N-glycoproteins, will facilitate follow-up biomarker research into early cancer detection.

2018 George Lyman Duff Memorial Lecture: Genetics and Genomics of Coronary Artery Disease: A Decade of Progress

Recent studies have led to a broader understanding of the genetic architecture of coronary artery disease and demonstrate that it largely derives from the cumulative effect of multiple common risk alleles individually of small effect size rather than rare variants with large effects on coronary artery disease risk. The tools applied include genome-wide association studies encompassing over 200 000 individuals complemented by bioinformatic approaches including imputation from whole-genome data sets, expression quantitative trait loci analyses, and interrogation of ENCODE (Encyclopedia of DNA Elements), Roadmap Epigenetic Project, and other data sets. Over 160 genome-wide significant loci associated with coronary artery disease risk have been identified using the genome-wide association studies approach, 90% of which are situated in intergenic regions. Here, I will describe, in part, our research over the last decade performed in collaboration with a series of bright trainees and an extensive number of groups and individuals around the world as it applies to our understanding of the genetic basis of this complex disease. These studies include computational approaches to better understand missing heritability and identify causal pathways, experimental approaches, and progress in understanding at the molecular level the function of the multiple risk loci identified and potential applications of these genomic data in clinical medicine and drug discovery.

Application of metabolomics technologies toward cancer prognosis and therapy

Altered metabolism is one of the defining features of cancer. Since the discovery of the Warburg effect in 1924, research into the metabolic aspects of cancer has been reinvigorated over the past decade. Metabolomics is an invaluable tool for gaining insights into numerous biochemical processes including those related to cancer metabolism and metabolic aspects of other diseases. The combination of untargeted and targeted metabolomics approaches has greatly facilitated the discovery of many cancer biomarkers with prognostic potential. Using mass spectrometry-based stable isotope-resolved metabolomics (SIRM) with isotopic labeling, a powerful tool used in pathway analysis, researchers have discovered novel cancer metabolic pathways and metabolic targets for therapeutic application. Metabolomics technologies provide invaluable metabolic insights reflecting cancer progression in coordination with genomics and proteomics aspects. The systematic study of metabolite levels in the metabolome and their dynamics within a biological organism has been, in recent years, applied across a wide range of fields. Metabolomics technologies have been applied to both early clinical trials and pre-clinical research in several essential aspects of human health. This chapter will give an overview of metabolomics technologies and their application in the discovery of novel pathways using isotopic labeled and non-labeled metabolomics.

Astrocyte-specific transcriptome analysis using the ALDH1L1 bacTRAP mouse reveals novel biomarkers of astrogliosis in response to neurotoxicity

Neurotoxicology is hampered by the inability to predict regional and cellular targets of toxicant-induced damage. Evaluating astrogliosis overcomes this problem because reactive astrocytes highlight the location of toxicant-induced damage. While enhanced expression of glial fibrillary acidic protein is a hallmark of astrogliosis, few other biomarkers have been identified. However, bacterial artificial chromosome - translating ribosome affinity purification (bacTRAP) technology allows for characterization of the actively translating transcriptome of a particular cell type; use of this technology in aldehyde dehydrogenase 1 family member L1 (ALDH1L1) bacTRAP mice can identify genes selectively expressed in astrocytes. The aim of this study was to characterize additional biomarkers of neurotoxicity-induced astrogliosis using ALDH1L1 bacTRAP mice. The known dopaminergic neurotoxicant 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP; 12.5 mg/kg s.c.) was used to induce astrogliosis. Striatal tissue was obtained 12, 24, and 48 h following exposure for the isolation of actively translating RNA. Subsequently, MPTP-induced changes in this RNA pool were analyzed by microarray and 184 statistically significant, differentially expressed genes were identified. The dataset was interrogated by gene ontology, pathway, and co-expression network analyses, which identified novel genes, as well as those with known immune and inflammatory functions. Using these analyses, we were directed to several genes associated with reactive astrocytes. Of these, TIMP1 and miR-147 were identified as candidate biomarkers because of their robust increased expression following both MPTP and trimethyl tin exposures. Thus, we have demonstrated that bacTRAP can be used to identify new biomarkers of astrogliosis and aid in the characterization of astrocyte phenotypes induced by toxicant exposures. OPEN SCIENCE BADGES: This article has received a badge for *Open Materials* because it provided all relevant information to reproduce the study in the manuscript. The complete Open Science Disclosure form for this article can be found at the end of the article. More information about the Open Practices badges can be found at https://cos.io/our-services/open-science-badges/. Cover Image for this issue: doi: 10.1111/jnc.14518.

A Multi-Cohort and Multi-Omics Meta-Analysis Framework to Identify Network-Based Gene Signatures

Although massive amounts of condition-specific molecular profiles are being accumulated in public repositories every day, meaningful interpretation of these data remains a major challenge. In an effort to identify the biomarkers that describe the key biological phenomena for a given condition, several approaches have been developed over the past few years. However, the majority of these approaches either (i) do not consider the known intermolecular interactions, or (ii) do not integrate molecular data of multiple types (e.g., genomics, transcriptomics, proteomics, epigenomics, etc.), and thus potentially fail to capture the true biological changes responsible for complex diseases (e.g., cancer). In addition, these approaches often ignore the heterogeneity and study bias present in independent molecular cohorts. In this manuscript, we propose a novel multi-cohort and multi-omics meta-analysis framework that overcomes all three limitations mentioned above in order to identify robust molecular subnetworks that capture the key dynamic nature of a given biological condition. Our framework integrates multiple independent gene expression studies, unmatched DNA methylation studies, and protein-protein interactions to identify methylation-driven subnetworks. We demonstrate the proposed framework by constructing subnetworks related to two complex diseases: glioblastoma and low-grade gliomas. We validate the identified subnetworks by showing their ability to predict patients' clinical outcome on multiple independent validation cohorts.

Loss of RNA-Binding Protein Sfpq Causes Long-Gene Transcriptopathy in Skeletal Muscle and Severe Muscle Mass Reduction with Metabolic Myopathy

Growing evidences are suggesting that extra-long genes in mammals are vulnerable for full-gene length transcription and dysregulation of long genes is a mechanism underlying human genetic disorders. How long-distance transcription is achieved is a fundamental question to be elucidated. In previous study, we had discovered that RNA-binding protein SFPQ preferentially binds to long pre-mRNAs and specifically regulates the cluster of neuronal genes >100 kbp. Here we investigated the roles of SFPQ for long gene expression, target specificities, and also physiological functions in skeletal muscle. Loss of Sfpq selectively downregulated genes >100 kbp including Dystrophin, which is 2.26 Mbp in length. Sfpq knockout (KO) mice showed progressive muscle mass reduction and metabolic myopathy characterized by glycogen accumulation and decreased abundance of mitochondrial oxidative phosphorylation complexes. Functional clustering analysis identified energy metabolism pathway genes as SFPQ's targets. These findings indicate target gene specificities and tissue-specific physiological functions of SFPQ in skeletal muscle.

•

SFPQ is essential for long gene expression, including Dystrophin, in skeletal muscle

•

Disruption of Sfpq caused severe muscle mass reduction and premature death

•

SFPQ is required for metabolic pathway gene expression in skeletal muscle

•

Loss of Sfpq decreased OXPHOS complexes and caused glycogen accumulation