DrDimont: explainable drug response prediction from differential analysis of multi-omics networks.

MOTIVATION: While it has been well established that drugs affect and help patients differently, personalized drug response predictions remain challenging. Solutions based on single omics measurements have been proposed, and networks provide means to incorporate molecular interactions into reasoning. However, how to integrate the wealth of information contained in multiple omics layers still poses a complex problem. RESULTS: We present DrDimont, Drug response prediction from Differential analysis of multi-omics networks. It allows for comparative conclusions between two conditions and translates them into differential drug response predictions. DrDimont focuses on molecular interactions. It establishes condition-specific networks from correlation within an omics layer that are then reduced and combined into heterogeneous, multi-omics molecular networks. A novel semi-local, path-based integration step ensures integrative conclusions. Differential predictions are derived from comparing the condition-specific integrated networks. DrDimont's predictions are explainable, i.e. molecular differences that are the source of high differential drug scores can be retrieved. We predict differential drug response in breast cancer using transcriptomics, proteomics, phosphosite and metabolomics measurements and contrast estrogen receptor positive and receptor negative patients. DrDimont performs better than drug prediction based on differential protein expression or PageRank when evaluating it on ground truth data from cancer cell lines. We find proteomic and phosphosite layers to carry most information for distinguishing drug response. AVAILABILITY AND IMPLEMENTATION: DrDimont is available on CRAN: https://cran.r-project.org/package=DrDimont. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Allergic airway inflammation delays glioblastoma progression and reinvigorates systemic and local immunity in mice.

BACKGROUND: Numerous patient-based studies have highlighted the protective role of immunoglobulin E-mediated allergic diseases on glioblastoma (GBM) susceptibility and prognosis. However, the mechanisms behind this observation remain elusive. Our objective was to establish a preclinical model able to recapitulate this phenomenon and investigate the role of immunity underlying such protection. METHODS: An immunocompetent mouse model of allergic airway inflammation (AAI) was initiated before intracranial implantation of mouse GBM cells (GL261). RAG1-KO mice served to assess tumor growth in a model deficient for adaptive immunity. Tumor development was monitored by MRI. Microglia were isolated for functional analyses and RNA-sequencing. Peripheral as well as tumor-associated immune cells were characterized by flow cytometry. The impact of allergy-related microglial genes on patient survival was analyzed by Cox regression using publicly available datasets. RESULTS: We found that allergy establishment in mice delayed tumor engraftment in the brain and reduced tumor growth resulting in increased mouse survival. AAI induced a transcriptional reprogramming of microglia towards a pro-inflammatory-like state, uncovering a microglia gene signature, which correlated with limited local immunosuppression in glioma patients. AAI increased effector memory T-cells in the circulation as well as tumor-infiltrating CD4+ T-cells. The survival benefit conferred by AAI was lost in mice devoid of adaptive immunity. CONCLUSION: Our results demonstrate that AAI limits both tumor take and progression in mice, providing a preclinical model to study the impact of allergy on GBM susceptibility and prognosis, respectively. We identify a potentiation of local and adaptive systemic immunity, suggesting a reciprocal crosstalk that orchestrates allergy-induced immune protection against GBM.

Machine learning for predicting neurodegenerative diseases in the general older population: a cohort study.

BACKGROUND: In the older general population, neurodegenerative diseases (NDs) are associated with increased disability, decreased physical and cognitive function. Detecting risk factors can help implement prevention measures. Using deep neural networks (DNNs), a machine-learning algorithm could be an alternative to Cox regression in tabular datasets with many predictive features. We aimed to compare the performance of different types of DNNs with regularized Cox proportional hazards models to predict NDs in the older general population. METHODS: We performed a longitudinal analysis with participants of the English Longitudinal Study of Ageing. We included men and women with no NDs at baseline, aged 60 years and older, assessed every 2 years from 2004 to 2005 (wave2) to 2016-2017 (wave 8). The features were a set of 91 epidemiological and clinical baseline variables. The outcome was new events of Parkinson's, Alzheimer or dementia. After applying multiple imputations, we trained three DNN algorithms: Feedforward, TabTransformer, and Dense Convolutional (Densenet). In addition, we trained two algorithms based on Cox models: Elastic Net regularization (CoxEn) and selected features (CoxSf). RESULTS: 5433 participants were included in wave 2. During follow-up, 12.7% participants developed NDs. Although the five models predicted NDs events, the discriminative ability was superior using TabTransformer (Uno's C-statistic (coefficient (95% confidence intervals)) 0.757 (0.702, 0.805). TabTransformer showed superior time-dependent balanced accuracy (0.834 (0.779, 0.889)) and specificity (0.855 (0.0.773, 0.909)) than the other models. With the CoxSf (hazard ratio (95% confidence intervals)), age (10.0 (6.9, 14.7)), poor hearing (1.3 (1.1, 1.5)) and weight loss 1.3 (1.1, 1.6)) were associated with a higher DNN risk. In contrast, executive function (0.3 (0.2, 0.6)), memory (0, 0, 0.1)), increased gait speed (0.2, (0.1, 0.4)), vigorous physical activity (0.7, 0.6, 0.9)) and higher BMI (0.4 (0.2, 0.8)) were associated with a lower DNN risk. CONCLUSION: TabTransformer is promising for prediction of NDs with heterogeneous tabular datasets with numerous features. Moreover, it can handle censored data. However, Cox models perform well and are easier to interpret than DNNs. Therefore, they are still a good choice for NDs.

Patient-derived organoids and orthotopic xenografts of primary and recurrent gliomas represent relevant patient avatars for precision oncology.

Patient-based cancer models are essential tools for studying tumor biology and for the assessment of drug responses in a translational context. We report the establishment a large cohort of unique organoids and patient-derived orthotopic xenografts (PDOX) of various glioma subtypes, including gliomas with mutations in IDH1, and paired longitudinal PDOX from primary and recurrent tumors of the same patient. We show that glioma PDOXs enable long-term propagation of patient tumors and represent clinically relevant patient avatars that retain histopathological, genetic, epigenetic, and transcriptomic features of parental tumors. We find no evidence of mouse-specific clonal evolution in glioma PDOXs. Our cohort captures individual molecular genotypes for precision medicine including mutations in IDH1, ATRX, TP53, MDM2/4, amplification of EGFR, PDGFRA, MET, CDK4/6, MDM2/4, and deletion of CDKN2A/B, PTCH, and PTEN. Matched longitudinal PDOX recapitulate the limited genetic evolution of gliomas observed in patients following treatment. At the histological level, we observe increased vascularization in the rat host as compared to mice. PDOX-derived standardized glioma organoids are amenable to high-throughput drug screens that can be validated in mice. We show clinically relevant responses to temozolomide (TMZ) and to targeted treatments, such as EGFR and CDK4/6 inhibitors in (epi)genetically defined subgroups, according to MGMT promoter and EGFR/CDK status, respectively. Dianhydrogalactitol (VAL-083), a promising bifunctional alkylating agent in the current clinical trial, displayed high therapeutic efficacy, and was able to overcome TMZ resistance in glioblastoma. Our work underscores the clinical relevance of glioma organoids and PDOX models for translational research and personalized treatment studies and represents a unique publicly available resource for precision oncology.

Single-cell transcriptomics reveals distinct inflammation-induced microglia signatures.

Microglia are specialized parenchymal-resident phagocytes of the central nervous system (CNS) that actively support, defend and modulate the neural environment. Dysfunctional microglial responses are thought to worsen CNS diseases; nevertheless, their impact during neuroinflammatory processes remains largely obscure. Here, using a combination of single-cell RNA sequencing and multicolour flow cytometry, we comprehensively profile microglia in the brain of lipopolysaccharide (LPS)-injected mice. By excluding the contribution of other immune CNS-resident and peripheral cells, we show that microglia isolated from LPS-injected mice display a global downregulation of their homeostatic signature together with an upregulation of inflammatory genes. Notably, we identify distinct microglial activated profiles under inflammatory conditions, which greatly differ from neurodegenerative disease-associated profiles. These results provide insights into microglial heterogeneity and establish a resource for the identification of specific phenotypes in CNS disorders, such as neuroinflammatory and neurodegenerative diseases.

Computational models for predicting drug responses in cancer research.

The computational prediction of drug responses based on the analysis of multiple types of genome-wide molecular data is vital for accomplishing the promise of precision medicine in oncology. This will benefit cancer patients by matching their tumor characteristics to the most effective therapy available. As larger and more diverse layers of patient-related data become available, further demands for new bioinformatics approaches and expertise will arise. This article reviews key strategies, resources and techniques for the prediction of drug sensitivity in cell lines and patient-derived samples. It discusses major advances and challenges associated with the different model development steps. This review highlights major trends in this area, and will assist researchers in the assessment of recent progress and in the selection of approaches to emerging applications in oncology.

Databases for lncRNAs: a comparative evaluation of emerging tools.

The vast majority of the human transcriptome does not code for proteins. Advances in transcriptome arrays and deep sequencing are giving rise to a fast accumulation of large data sets, particularly of long noncoding RNAs (lncRNAs). Although it is clear that individual lncRNAs may play important and diverse biological roles, there is a large gap between the number of existing lncRNAs and their known relation to molecular/cellular function. This and related information have recently been gathered in several databases dedicated to lncRNA research. Here, we review the content of general and more specialized databases on lncRNAs. We evaluate these resources in terms of the quality of annotations, the reporting of validated or predicted molecular associations, and their integration with other resources and computational analysis tools. We illustrate our findings using known and novel cancer-related lncRNAs. Finally, we discuss limitations and highlight potential future directions for these databases to help delineating functions associated with lncRNAs.

Permanent culture of macrophages at physiological oxygen attenuates the antioxidant and immunomodulatory properties of dimethyl fumarate.

We hypothesized that O2 tension influences the redox state and the immunomodulatory responses of inflammatory cells to dimethyl fumarate (DMF), an activator of the nuclear factor Nrf2 that controls antioxidant genes expression. This concept was investigated in macrophages permanently cultured at either physiological (5% O2) or atmospheric (20% O2) oxygen levels and then treated with DMF or challenged with lipopolysaccharide (LPS) to induce inflammation. RAW 264.7 macrophages cultured at 20% O2 exhibited a pro-oxidant phenotype, reflected by a lower content of reduced glutathione, higher oxidized glutathione and increased production of reactive oxygen species when compared to macrophages continuously grown at 5% O2. At 20% O2, DMF induced a stronger antioxidant response compared to 5% O2 as evidenced by a higher expression of heme oxygenase-1, NAD(P)H:quinone oxydoreductase-1 and superoxide dismutase-2. After challenge of macrophages with LPS, several pro-inflammatory (iNOS, TNF-α, MMP-2, MMP-9), anti-inflammatory (arginase-1, IL-10) and pro-angiogenic (VEGF-A) mediators were evaluated in the presence or absence of DMF. All markers, with few interesting exceptions, were significantly reduced at 5% O2. This study brings new insights on the effects of O2 in the cellular adaptation to oxidative and inflammatory stimuli and highlights the importance of characterizing the effects of chemicals and drugs at physiologically relevant O2 tension. Our results demonstrate that the common practice of culturing cells at atmospheric O2 drives the endogenous cellular environment towards an oxidative stress phenotype, affecting inflammation and the expression of antioxidant pathways by exogenous modulators.

Therapeutic control and resistance of the EGFR-driven signaling network in glioblastoma.

The alteration of the epidermal growth factor receptor (EGFR)-driven signaling network is a characteristic feature of glioblastomas (GBM), and its inhibition represents a treatment strategy. However, EGFR-targeted interventions have been largely ineffective. Complex perturbations in this system are likely to be central to tumor cells with high adaptive capacity and resistance to therapies. We review key concepts and mechanisms relevant to EGFR-targeted treatment resistance at a systems level. Our understanding of treatment resistance as a systems-level phenomenon is necessary to develop effective therapeutic options for GBM patients. This is allowing us to go beyond the notion of therapeutic targets as single molecular components, into strategies that can weaken cancer signaling robustness and boost inherent network-level vulnerabilities.

Transcriptional response to cardiac injury in the zebrafish: systematic identification of genes with highly concordant activity across in vivo models.

BACKGROUND: Zebrafish is a clinically-relevant model of heart regeneration. Unlike mammals, it has a remarkable heart repair capacity after injury, and promises novel translational applications. Amputation and cryoinjury models are key research tools for understanding injury response and regeneration in vivo. An understanding of the transcriptional responses following injury is needed to identify key players of heart tissue repair, as well as potential targets for boosting this property in humans. RESULTS: We investigated amputation and cryoinjury in vivo models of heart damage in the zebrafish through unbiased, integrative analyses of independent molecular datasets. To detect genes with potential biological roles, we derived computational prediction models with microarray data from heart amputation experiments. We focused on a top-ranked set of genes highly activated in the early post-injury stage, whose activity was further verified in independent microarray datasets. Next, we performed independent validations of expression responses with qPCR in a cryoinjury model. Across in vivo models, the top candidates showed highly concordant responses at 1 and 3 days post-injury, which highlights the predictive power of our analysis strategies and the possible biological relevance of these genes. Top candidates are significantly involved in cell fate specification and differentiation, and include heart failure markers such as periostin, as well as potential new targets for heart regeneration. For example, ptgis and ca2 were overexpressed, while usp2a, a regulator of the p53 pathway, was down-regulated in our in vivo models. Interestingly, a high activity of ptgis and ca2 has been previously observed in failing hearts from rats and humans. CONCLUSIONS: We identified genes with potential critical roles in the response to cardiac damage in the zebrafish. Their transcriptional activities are reproducible in different in vivo models of cardiac injury.

Computational discrete models of tissue growth and regeneration.

Tissue growth and regeneration are fundamental processes underpinning crucial physiological and pathological conditions: ranging from normal blood vessel network development, response to stem cells therapy and cancers. Modelling of such biological phenomena has been addressed through mathematical and algorithmic approaches. The former implements continuous representations based on differential equations. The latter exploit operational descriptions in the form of computing programs to represent and execute the models. Within this area, models that define the cell as the fundamental unit of model development, as well as discrete representations of different model entities, are important to plan in vitro experiments and to generate new testable hypotheses. This article reviews the application of algorithmic discrete models, with a focus on tissue growth and regeneration phenomena in the context of health and disease. The review begins with an overview of basic concepts, problems and approaches of computational discrete models. This will include a discussion of basic assumptions and design principles. An overview of key cell-driven approaches and examples of applications in tissue growth and regeneration is provided. The specification, implementation and analysis of a model are illustrated with a hypothetical example, which mimics the branching and sprouting patterns observed in blood vessel network development. The article concludes with a discussion of current challenges and recommendations.

Gene set analysis in the cloud.

UNLABELLED: Cloud computing offers low cost and highly flexible opportunities in bioinformatics. Its potential has already been demonstrated in high-throughput sequence data analysis. Pathway-based or gene set analysis of expression data has received relatively less attention. We developed a gene set analysis algorithm for biomarker identification in the cloud. The resulting tool, YunBe, is ready to use on Amazon Web Services. Moreover, here we compare its performance to those obtained with desktop and computing cluster solutions. AVAILABILITY AND IMPLEMENTATION: YunBe is open-source and freely accessible within the Amazon Elastic MapReduce service at s3n://lrcv-crp-sante/app/yunbe.jar. Source code and user's guidelines can be downloaded from http://tinyurl.com/yunbedownload.

Prognostic transcriptional association networks: a new supervised approach based on regression trees.

MOTIVATION: The application of information encoded in molecular networks for prognostic purposes is a crucial objective of systems biomedicine. This approach has not been widely investigated in the cardiovascular research area. Within this area, the prediction of clinical outcomes after suffering a heart attack would represent a significant step forward. We developed a new quantitative prediction-based method for this prognostic problem based on the discovery of clinically relevant transcriptional association networks. This method integrates regression trees and clinical class-specific networks, and can be applied to other clinical domains. RESULTS: Before analyzing our cardiovascular disease dataset, we tested the usefulness of our approach on a benchmark dataset with control and disease patients. We also compared it to several algorithms to infer transcriptional association networks and classification models. Comparative results provided evidence of the prediction power of our approach. Next, we discovered new models for predicting good and bad outcomes after myocardial infarction. Using blood-derived gene expression data, our models reported areas under the receiver operating characteristic curve above 0.70. Our model could also outperform different techniques based on co-expressed gene modules. We also predicted processes that may represent novel therapeutic targets for heart disease, such as the synthesis of leucine and isoleucine. AVAILABILITY: The SATuRNo software is freely available at http://www.lsi.us.es/isanepo/toolsSaturno/.

Use of circulating microRNAs to diagnose acute myocardial infarction.

BACKGROUND: Rapid and correct diagnosis of acute myocardial infarction (MI) has an important impact on patient treatment and prognosis. We compared the diagnostic performance of high-sensitivity cardiac troponin T (hs-cTnT) and cardiac enriched microRNAs (miRNAs) in patients with MI. METHODS: Circulating concentrations of cardiac-enriched miR-208b and miR-499 were measured by quantitative PCR in a case-control study of 510 MI patients referred for primary mechanical reperfusion and 87 healthy controls. RESULTS: miRNA-208b and miR-499 were highly increased in MI patients (>10(5)-fold, P < 0.001) and nearly undetectable in healthy controls. Patients with ST-elevation MI (n= 397) had higher miRNA concentrations than patients with non-ST-elevation MI (n = 113) (P < 0.001). Both miRNAs correlated with peak concentrations of creatine kinase and cTnT (P < 10(-9)). miRNAs and hs-cTnT were already detectable in the plasma 1 h after onset of chest pain. In patients who presented <3 h after onset of pain, miR-499 was positive in 93% of patients and hs-cTnT in 88% of patients (P= 0.78). Overall, miR-499 and hs-cTnT provided comparable diagnostic value with areas under the ROC curves of 0.97. The reclassification index of miR-499 to a clinical model including several risk factors and hs-cTnT was not significant (P = 0.15). CONCLUSION: Circulating miRNAs are powerful markers of acute MI. Their usefulness in the establishment of a rapid and accurate diagnosis of acute MI remains to be determined in unselected populations of patients with acute chest pain.

Low levels of vascular endothelial growth factor B predict left ventricular remodeling after acute myocardial infarction.

BACKGROUND: Left ventricular (LV) remodeling is a prognostically important development after acute myocardial infarction (AMI). We recently reported that vascular endothelial growth factor B (VEGFB) may be a potential new biomarker of LV remodeling. This potential biomarker was evaluated in the present study. METHODS AND RESULTS: Patients with AMI (n = 290) and healthy volunteers (n = 42) were included. Plasma VEGFB levels were assessed before discharge. LV remodeling was determined by echocardiography at 6 months' follow-up. Levels of VEGFB were elevated in AMI patients compared with healthy volunteers (1.5-fold; P = .001). Mean plasma levels of VEGFB were 64% higher (P < .001) in patients in whom LV end-diastolic volume (EDV) decreased during follow-up (ΔEDV ≤ 0; n = 144; reverse remodeling) compared with patients in whom ΔEDV increased (ΔEDV > 0; n = 146; remodeling). Using logistic regression models, independent relationships were found between VEGFB (odds ratio [OR] 0.8, 95% confidence interval [CI] 0.7-0.9; P = .0007) and infarct territory (OR 1.7, 95% CI 1.1-2.8; P = .02). Patients with anterior MI and low levels of VEGFB had the highest risk of remodeling. VEFGB outperformed N-terminal pro-B-type natriuretic peptide to predict LV remodeling, and low levels of VEGFB (<100 pg/mL) provided a specificity of 90%. Adding VEGFB to a clinical model involving age, sex, smoking habit, and infarct territory resulted in a net reclassification index of 11.7%. CONCLUSIONS: Plasma levels of VEGFB increase after AMI and correlate with preservation of cardiac function. Low levels of VEGFB accurately predict LV remodeling. Therefore, circulating VEGFB may have clinical utility in the identification of patients at high risk of remodeling after AMI.

Computational biology for cardiovascular biomarker discovery.

Computational biology is essential in the process of translating biological knowledge into clinical practice, as well as in the understanding of biological phenomena based on the resources and technologies originating from the clinical environment. One such key contribution of computational biology is the discovery of biomarkers for predicting clinical outcomes using 'omic' information. This process involves the predictive modelling and integration of different types of data and knowledge for screening, diagnostic or prognostic purposes. Moreover, this requires the design and combination of different methodologies based on statistical analysis and machine learning. This article introduces key computational approaches and applications to biomarker discovery based on different types of 'omic' data. Although we emphasize applications in cardiovascular research, the computational requirements and advances discussed here are also relevant to other domains. We will start by introducing some of the contributions of computational biology to translational research, followed by an overview of methods and technologies used for the identification of biomarkers with predictive or classification value. The main types of 'omic' approaches to biomarker discovery will be presented with specific examples from cardiovascular research. This will include a review of computational methodologies for single-source and integrative data applications. Major computational methods for model evaluation will be described together with recommendations for reporting models and results. We will present recent advances in cardiovascular biomarker discovery based on the combination of gene expression and functional network analyses. The review will conclude with a discussion of key challenges for computational biology, including perspectives from the biosciences and clinical areas.

Ontology- and graph-based similarity assessment in biological networks.

SUMMARY: A standard systems-based approach to biomarker and drug target discovery consists of placing putative biomarkers in the context of a network of biological interactions, followed by different 'guilt-by-association' analyses. The latter is typically done based on network structural features. Here, an alternative analysis approach in which the networks are analyzed on a 'semantic similarity' space is reported. Such information is extracted from ontology-based functional annotations. We present SimTrek, a Cytoscape plugin for ontology-based similarity assessment in biological networks. AVAILABILITY: http://rosalind.infj.ulst.ac.uk/SimTrek.html CONTACT: francisco.azuaje@crp-sante.lu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.