The species translation challenge-a systems biology perspective on human and rat bronchial epithelial cells

Ex Vivo and In Vitro Proteomic Approach to Elucidate the Relevance of IL-4 and IL-10 in Intervertebral Disc Pathophysiology

Background

This study investigates the native presence and potential anabolic effects of interleukin (IL)-4 and IL-10 in the human intervertebral disc (IVD).

Methods

Human nucleus pulposus (NP) cells cultured in 3D from trauma and degenerate IVDs and NP explants were stimulated with 10 ng/mL IL-4, IL-10, or each in combination with 1 ng/mL IL-1β stimulation. The role of IL-4 and IL-10 in the IVD was evaluated using immunohistochemistry, gene expression, and Luminex multiplex immunoassay proteomics (73 secreted) and phosphoproteomics (21 phosphorylated proteins).

Results

IL-4, IL-4R, and IL-10R expression and localization in human cartilage endplate tissue were demonstrated for the first time. No significant gene expression changes were noted under IL-4 or IL-10 stimulation. However, IL-1β stimulation significantly increased MMP3, COX2, TIMP1, and TRPV4 expression in NP cells from trauma IVDs. Combined IL-4 and IL-1β treatment induced a significant increase in protein secretion of IL-1α, IL-7, IL-16, IL-17F, IL-18, IFNγ, TNF, ST2, PROK1, bFGF2, and stem cell factor exclusively in NP cells from degenerated IVDs. Conversely, the secretome profile of explants revealed an IL-4-mediated decrease in CXCL13 following treatment with IL-1β. Combined IL-10 and IL-1β treatment increased neurotrophic growth factor secretion compared with IL-10 baseline.

Conclusions

The NP cell phenotype affects the pleiotropic role of IL-4, which can induce a pro-inflammatory response in the presence of catabolic stimuli and enhance the effects of IL-1β in degenerated IVDs. Environmental factors, including 3D culture and hypoxia, may alter IL-4's role. Finally, IL-10's potential neurotrophic effects under catabolic stimuli warrant further investigation to clarify its role in IVD degeneration.

Systematic analysis of human colorectal cancer scRNA-seq revealed limited pro-tumoral IL-17 production potential in gamma delta T cells

Gamma delta T cells play a crucial role in anti-tumor immunity due to their cytotoxic properties. However, the role and extent of γδ T cells in production of pro-tumorigenic interleukin-17 (IL-17) within the tumor microenvironment of colorectal cancer (CRC) remains controversial. In this study, we re-analyzed nine published human CRC whole-tissue single-cell RNA sequencing datasets, identifying 18,483 γδ T cells out of 951,785 total cells, in the neoplastic or adjacent normal tissue of 165 human CRC patients. Our results confirm that tumor-infiltrating γδ T cells exhibit high cytotoxicity-related transcription in both tumor and adjacent normal tissues, but critically, none of the γδ T cell clusters showed IL-17 production potential. We also identified various γδ T cell subsets, including poised effector-like T cells, tissue-resident memory T cells, progenitor exhausted-like T cells, and exhausted T cells, and noted an increased expression of cytotoxic molecules in tumor-infiltrating γδ T cells compared to their normal area counterparts. We proposed anti-tumor γδ T effector cells may arise from tissue-resident progenitor cells based on the trajectory analysis. Our work demonstrates that γδ T cells in CRC primarily function as cytotoxic effector cells rather than IL-17 producers, mitigating the concerns about their potential pro-tumorigenic roles in CRC, highlighting the importance of accurately characterizing these cells for cancer immunotherapy research and the unneglectable cross-species discrepancy between the mouse and human immune system in the study of cancer immunology.

Systematic Analysis of Human Colorectal Cancer scRNA-seq Revealed Limited Pro-tumoral IL-17 Production Potential in Gamma Delta T Cells

Gamma delta (γδ) T cells play a crucial role in anti-tumor immunity due to their cytotoxic properties. However, the role and extent of γδ T cells in production of pro-tumorigenic interleukin- 17 (IL-17) within the tumor microenvironment (TME) of colorectal cancer (CRC) remains controversial. In this study, we re-analyzed nine published human CRC whole-tissue single-cell RNA sequencing (scRNA-seq) datasets, identifying 18,483 γδ T cells out of 951,785 total cells, in the neoplastic or adjacent normal tissue of 165 human CRC patients. Our results confirm that tumor-infiltrating γδ T cells exhibit high cytotoxicity-related transcription in both tumor and adjacent normal tissues, but critically, none of the γδ T cell clusters showed IL-17 production potential. We also identified various γδ T cell subsets, including Teff, TRM, Tpex, and Tex, and noted an increased expression of cytotoxic molecules in tumor-infiltrating γδ T cells compared to their normal area counterparts. Our work demonstrates that γδ T cells in CRC primarily function as cytotoxic effector cells rather than IL-17 producers, mitigating the concerns about their potential pro-tumorigenic roles in CRC, highlighting the importance of accurately characterizing these cells for cancer immunotherapy research and the unneglectable cross-species discrepancy between the mouse and human immune system in the study of cancer immunology.

Comparative toxicological assessment of cigarettes and new category products via an in vitro multiplex proteomics platform

Cigarette smoking is a risk factor for several diseases such as cancer, cardiovascular disease (CVD), and chronic obstructive pulmonary diseases (COPD), however, the underlying mechanisms are not fully understood. Alternative nicotine products with reduced risk potential (RRPs) including tobacco heating products (THPs), and e-cigarettes have recently emerged as viable alternatives to cigarettes that may contribute to the overall strategy of tobacco harm reduction due to the significantly lower levels of toxicants in these products' emissions as compared to cigarette smoke. Assessing the effects of RRPs on biological responses is important to demonstrate the potential value of RRPs towards tobacco harm reduction. Here, we evaluated the inflammatory and signaling responses of human lung epithelial cells to aqueous aerosol extracts (AqE) generated from the 1R6F reference cigarette, the glo™ THP, and the Vype ePen 3.0 e-cigarette using multiplex analysis of 37 inflammatory and phosphoprotein markers. Cellular exposure to the different RRPs and 1R6F AqEs resulted in distinct response profiles with 1R6F being the most biologically active followed by glo™ and ePen 3.0. 1R6F activated stress-related and pro-survival markers c-JUN, CREB1, p38 MAPK and MEK1 and led to the release of IL-1α. glo™ activated MEK1 and decreased IL-1β levels, whilst ePen 3.0 affected IL-1β levels but had no effect on the signaling activity compared to untreated cells. Our results demonstrated the reduced biological effect of RRPs and suggest that targeted analysis of inflammatory and cell signaling mediators is a valuable tool for the routine assessment of RRPs.

Multi-endpoint in vitro toxicological assessment of snus and tobacco-free nicotine pouch extracts

'Modern' oral tobacco-free nicotine pouches (NPs) are a nicotine containing product similar in appearance and concept to Swedish snus. A three-step approach was taken to analyse the biological effects of NPs and snus extracts in vitro. ToxTracker was used to screen for biomarkers for oxidative stress, cell stress, protein damage and DNA damage. Cytotoxicity, mutagenicity, and genotoxicity were assessed in the following respective assays: Neutral Red Uptake (NRU), Ames and Mouse Lymphoma Assay (MLA). Targeted analysis of phosphorylation signalling and inflammatory markers under non-toxic conditions was used to investigate any potential signalling pathways or inflammatory response. A reference snus (CRP1.1) and four NPs with various flavours and nicotine strengths were assessed. Test article extracts was generated by incubating one pouch in 20 mL of media (specific to each assay) with the inclusion of the pouch material. NP extracts did not induce any cytotoxicity or mutagenic response, genotoxic response was minimal and limited signalling or inflammatory markers were induced. In contrast, CRP1.1 induced a positive response in four toxicological endpoints in the absence of S9: Srxn1 (oxidative stress), Btg2 (cell stress), Ddit3 (protein damage) and Rtkn (DNA damage), and three endpoints in presence of S9: Srxn1, Ddit3 and Rtkn. CRP1.1 was genotoxic when assessed in MLA and activated signalling pathways involved in proliferation and cellular stress and specifically induced phosphorylation of c-JUN, CREB1, p53, p38 MAPK and to a lesser extent AKT1S1, GSK3α/β, ERK1/2 and RSK1 in a dose-dependent manner. CRP 1.1 extracts resulted in the release of several inflammatory mediators including cytokines IL-1α, IL5, IL6, IL8, IL-1RA, MIF and TNF-β, receptor IL-2RA, and growth factors FGF-basic, VEGF and M-CSF. In conclusion these assays contribute to the weight of evidence assessment of the potential comparative health risks of NPs and snus.

Protein kinase D drives the secretion of invasion mediators in triple-negative breast cancer cell lines

The protein kinase D (PKD) family members regulate the fission of cargo vesicles at the Golgi complex and play a pro-oncogenic role in triple-negative breast cancer (TNBC). Whether PKD facilitates the secretion of tumor-promoting factors in TNBC, however, is still unknown. Using the pharmacological inhibition of PKD activity and siRNA-mediated depletion of PKD2 and PKD3, we identified the PKD-dependent secretome of the TNBC cell lines MDA-MB-231 and MDA-MB-468. Mass spectrometry-based proteomics and antibody-based assays revealed a significant downregulation of extracellular matrix related proteins and pro-invasive factors such as LIF, MMP-1, MMP-13, IL-11, M-CSF and GM-CSF in PKD-perturbed cells. Notably, secretion of these proteins in MDA-MB-231 cells was predominantly controlled by PKD2 and enhanced spheroid invasion. Consistently, PKD-dependent secretion of pro-invasive factors was more pronounced in metastatic TNBC cell lines. Our study thus uncovers a novel role of PKD2 in releasing a pro-invasive secretome.

AutoTransOP: translating omics signatures without orthologue requirements using deep learning

The development of therapeutics and vaccines for human diseases requires a systematic understanding of human biology. Although animal and in vitro culture models can elucidate some disease mechanisms, they typically fail to adequately recapitulate human biology as evidenced by the predominant likelihood of clinical trial failure. To address this problem, we developed AutoTransOP, a neural network autoencoder framework, to map omics profiles from designated species or cellular contexts into a global latent space, from which germane information for different contexts can be identified without the typically imposed requirement of matched orthologues. This approach was found in general to perform at least as well as current alternative methods in identifying animal/culture-specific molecular features predictive of other contexts-most importantly without requiring homology matching. For an especially challenging test case, we successfully applied our framework to a set of inter-species vaccine serology studies, where 1-to-1 mapping between human and non-human primate features does not exist.

Α Humanized RANKL Transgenic Mouse Model of Progestin-Induced Mammary Carcinogenesis for Evaluation of Novel Therapeutics

Receptor activator of nuclear factor-κB ligand (RANKL) is critically involved in mammary gland pathophysiology, while its pharmaceutical inhibition is being currently investigated in breast cancer. Herein, we investigated whether the overexpression of human RANKL in transgenic mice affects hormone-induced mammary carcinogenesis, and evaluated the efficacy of anti-RANKL treatments, such as OPG-Fc targeting both human and mouse RANKL or Denosumab against human RANKL. We established novel MPA/DMBA-driven mammary carcinogenesis models in TgRANKL mice that express both human and mouse RANKL, as well as in humanized humTgRANKL mice expressing only human RANKL, and compared them to MPA/DMBA-treated wild-type (WT) mice. Our results show that TgRANKL and WT mice have similar levels of susceptibility to mammary carcinogenesis, while OPG-Fc treatment restored mammary ductal density, and prevented ductal branching and the formation of neoplastic foci in both genotypes. humTgRANKL mice also developed MPA/DMBA-induced tumors with similar incidence and burden to those of WT and TgRANKL mice. The prophylactic treatment of humTgRANKL mice with Denosumab significantly prevented the rate of appearance of mammary tumors from 86.7% to 15.4% and the early stages of carcinogenesis, whereas therapeutic treatment did not lead to any significant attenuation of tumor incidence or tumor burden compared to control mice, suggesting the importance of RANKL primarily in the initial stages of tumorigenesis. Overall, we provide unique genetic tools for investigating the involvement of RANKL in breast carcinogenesis, and allow the preclinical evaluation of novel therapeutics that target hormone-related breast cancers.

A Combination of Conformation-Specific RAF Inhibitors Overcome Drug Resistance Brought about by RAF Overexpression

Cancer cells often adapt to targeted therapies, yet the molecular mechanisms underlying adaptive resistance remain only partially understood. Here, we explore a mechanism of RAS/RAF/MEK/ERK (MAPK) pathway reactivation through the upregulation of RAF isoform (RAFs) abundance. Using computational modeling and in vitro experiments, we show that the upregulation of RAFs changes the concentration range of paradoxical pathway activation upon treatment with conformation-specific RAF inhibitors. Additionally, our data indicate that the signaling output upon loss or downregulation of one RAF isoform can be compensated by overexpression of other RAF isoforms. We furthermore demonstrate that, while single RAF inhibitors cannot efficiently inhibit ERK reactivation caused by RAF overexpression, a combination of two structurally distinct RAF inhibitors synergizes to robustly suppress pathway reactivation.

Crowdsourced benchmarking of taxonomic metagenome profilers: lessons learned from the sbv IMPROVER Microbiomics challenge

Background

Selection of optimal computational strategies for analyzing metagenomics data is a decisive step in determining the microbial composition of a sample, and this procedure is complex because of the numerous tools currently available. The aim of this research was to summarize the results of crowdsourced sbv IMPROVER Microbiomics Challenge designed to evaluate the performance of off-the-shelf metagenomics software as well as to investigate the robustness of these results by the extended post-challenge analysis. In total 21 off-the-shelf taxonomic metagenome profiling pipelines were benchmarked for their capacity to identify the microbiome composition at various taxon levels across 104 shotgun metagenomics datasets of bacterial genomes (representative of various microbiome samples) from public databases. Performance was determined by comparing predicted taxonomy profiles with the gold standard.

Results

Most taxonomic profilers performed homogeneously well at the phylum level but generated intermediate and heterogeneous scores at the genus and species levels, respectively. kmer-based pipelines using Kraken with and without Bracken or using CLARK-S performed best overall, but they exhibited lower precision than the two marker-gene-based methods MetaPhlAn and mOTU. Filtering out the 1% least abundance species—which were not reliably predicted—helped increase the performance of most profilers by increasing precision but at the cost of recall. However, the use of adaptive filtering thresholds determined from the sample’s Shannon index increased the performance of most kmer-based profilers while mitigating the tradeoff between precision and recall.

Conclusions

kmer-based metagenomic pipelines using Kraken/Bracken or CLARK-S performed most robustly across a large variety of microbiome datasets. Removing non-reliably predicted low-abundance species by using diversity-dependent adaptive filtering thresholds further enhanced the performance of these tools. This work demonstrates the applicability of computational pipelines for accurately determining taxonomic profiles in clinical and environmental contexts and exemplifies the power of crowdsourcing for unbiased evaluation.

Supplementary Information

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

Stress-induced inflammation evoked by immunogenic cell death is blunted by the IRE1α kinase inhibitor KIRA6 through HSP60 targeting

Mounting evidence indicates that immunogenic therapies engaging the unfolded protein response (UPR) following endoplasmic reticulum (ER) stress favor proficient cancer cell-immune interactions, by stimulating the release of immunomodulatory/proinflammatory factors by stressed or dying cancer cells. UPR-driven transcription of proinflammatory cytokines/chemokines exert beneficial or detrimental effects on tumor growth and antitumor immunity, but the cell-autonomous machinery governing the cancer cell inflammatory output in response to immunogenic therapies remains poorly defined. Here, we profiled the transcriptome of cancer cells responding to immunogenic or weakly immunogenic treatments. Bioinformatics-driven pathway analysis indicated that immunogenic treatments instigated a NF-κB/AP-1-inflammatory stress response, which dissociated from both cell death and UPR. This stress-induced inflammation was specifically abolished by the IRE1α-kinase inhibitor KIRA6. Supernatants from immunogenic chemotherapy and KIRA6 co-treated cancer cells were deprived of proinflammatory/chemoattractant factors and failed to mobilize neutrophils and induce dendritic cell maturation. Furthermore, KIRA6 significantly reduced the in vivo vaccination potential of dying cancer cells responding to immunogenic chemotherapy. Mechanistically, we found that the anti-inflammatory effect of KIRA6 was still effective in IRE1α-deficient cells, indicating a hitherto unknown off-target effector of this IRE1α-kinase inhibitor. Generation of a KIRA6-clickable photoaffinity probe, mass spectrometry, and co-immunoprecipitation analysis identified cytosolic HSP60 as a KIRA6 off-target in the IKK-driven NF-κB pathway. In sum, our study unravels that HSP60 is a KIRA6-inhibitable upstream regulator of the NF-κB/AP-1-inflammatory stress responses evoked by immunogenic treatments. It also urges caution when interpreting the anti-inflammatory action of IRE1α chemical inhibitors.

Machine learning applications for therapeutic tasks with genomics data

Thanks to the increasing availability of genomics and other biomedical data, many machine learning algorithms have been proposed for a wide range of therapeutic discovery and development tasks. In this survey, we review the literature on machine learning applications for genomics through the lens of therapeutic development. We investigate the interplay among genomics, compounds, proteins, electronic health records, cellular images, and clinical texts. We identify 22 machine learning in genomics applications that span the whole therapeutics pipeline, from discovering novel targets, personalizing medicine, developing gene-editing tools, all the way to facilitating clinical trials and post-market studies. We also pinpoint seven key challenges in this field with potentials for expansion and impact. This survey examines recent research at the intersection of machine learning, genomics, and therapeutic development.

Cadherin-11, Sparc-related modular calcium binding protein-2, and Pigment epithelium-derived factor are promising non-invasive biomarkers of kidney fibrosis

Kidney fibrosis constitutes the shared final pathway of nearly all chronic nephropathies, but biomarkers for the non-invasive assessment of kidney fibrosis are currently not available. To address this, we characterize five candidate biomarkers of kidney fibrosis: Cadherin-11 (CDH11), Sparc-related modular calcium binding protein-2 (SMOC2), Pigment epithelium-derived factor (PEDF), Matrix-Gla protein, and Thrombospondin-2. Gene expression profiles in single-cell and single-nucleus RNA-sequencing (sc/snRNA-seq) datasets from rodent models of fibrosis and human chronic kidney disease (CKD) were explored, and Luminex-based assays for each biomarker were developed. Plasma and urine biomarker levels were measured using independent prospective cohorts of CKD: the Boston Kidney Biopsy Cohort, a cohort of individuals with biopsy-confirmed semiquantitative assessment of kidney fibrosis, and the Seattle Kidney Study, a cohort of patients with common forms of CKD. Ordinal logistic regression and Cox proportional hazards regression models were used to test associations of biomarkers with interstitial fibrosis and tubular atrophy and progression to end-stage kidney disease and death, respectively. Sc/snRNA-seq data confirmed cell-specific expression of biomarker genes in fibroblasts. After multivariable adjustment, higher levels of plasma CDH11, SMOC2, and PEDF and urinary CDH11 and PEDF were significantly associated with increasing severity of interstitial fibrosis and tubular atrophy in the Boston Kidney Biopsy Cohort. In both cohorts, higher levels of plasma and urinary SMOC2 and urinary CDH11 were independently associated with progression to end-stage kidney disease. Higher levels of urinary PEDF associated with end-stage kidney disease in the Seattle Kidney Study, with a similar signal in the Boston Kidney Biopsy Cohort, although the latter narrowly missed statistical significance. Thus, we identified CDH11, SMOC2, and PEDF as promising non-invasive biomarkers of kidney fibrosis.

Prediction of combination therapies based on topological modeling of the immune signaling network in multiple sclerosis

Background

Multiple sclerosis (MS) is a major health problem, leading to a significant disability and patient suffering. Although chronic activation of the immune system is a hallmark of the disease, its pathogenesis is poorly understood, while current treatments only ameliorate the disease and may produce severe side effects.

Methods

Here, we applied a network-based modeling approach based on phosphoproteomic data to uncover the differential activation in signaling wiring between healthy donors, untreated patients, and those under different treatments. Based in the patient-specific networks, we aimed to create a new approach to identify drug combinations that revert signaling to a healthy-like state. We performed ex vivo multiplexed phosphoproteomic assays upon perturbations with multiple drugs and ligands in primary immune cells from 169 subjects (MS patients, n=129 and matched healthy controls, n=40). Patients were either untreated or treated with fingolimod, natalizumab, interferon-β, glatiramer acetate, or the experimental therapy epigallocatechin gallate (EGCG). We generated for each donor a dynamic logic model by fitting a bespoke literature-derived network of MS-related pathways to the perturbation data. Last, we developed an approach based on network topology to identify deregulated interactions whose activity could be reverted to a “healthy-like” status by combination therapy. The experimental autoimmune encephalomyelitis (EAE) mouse model of MS was used to validate the prediction of combination therapies.

Results

Analysis of the models uncovered features of healthy-, disease-, and drug-specific signaling networks. We predicted several combinations with approved MS drugs that could revert signaling to a healthy-like state. Specifically, TGF-β activated kinase 1 (TAK1) kinase, involved in Transforming growth factor β-1 proprotein (TGF-β), Toll-like receptor, B cell receptor, and response to inflammation pathways, was found to be highly deregulated and co-druggable with all MS drugs studied. One of these predicted combinations, fingolimod with a TAK1 inhibitor, was validated in an animal model of MS.

Conclusions

Our approach based on donor-specific signaling networks enables prediction of targets for combination therapy for MS and other complex diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00925-8.

An interspecies translation model implicates integrin signaling in infliximab-resistant inflammatory bowel disease

Anti-tumor necrosis factor (anti-TNF) therapy resistance is a major clinical challenge in inflammatory bowel disease (IBD), due, in part, to insufficient understanding of disease-site, protein-level mechanisms. Although proteomics data from IBD mouse models exist, data and phenotype discrepancies contribute to confounding translation from preclinical animal models of disease to clinical cohorts. We developed an approach called translatable components regression (TransComp-R) to overcome interspecies and trans-omic discrepancies between mouse models and human subjects. TransComp-R combines mouse proteomic data with patient pretreatment transcriptomic data to identify molecular features discernable in the mouse data that are predictive of patient response to therapy. Interrogating the TransComp-R models revealed activated integrin pathway signaling in patients with anti-TNF-resistant colonic Crohn's disease (cCD) and ulcerative colitis (UC). As a step toward validation, we performed single-cell RNA sequencing (scRNA-seq) on biopsies from a patient with cCD and analyzed publicly available immune cell proteomics data to characterize the immune and intestinal cell types contributing to anti-TNF resistance. We found that ITGA1 was expressed in T cells and that interactions between these cells and intestinal cell types were associated with resistance to anti-TNF therapy. We experimentally showed that the α₁ integrin subunit mediated the effectiveness of anti-TNF therapy in human immune cells. Thus, TransComp-R identified an integrin signaling mechanism with potential therapeutic implications for overcoming anti-TNF therapy resistance. We suggest that TransComp-R is a generalizable framework for addressing species, molecular, and phenotypic discrepancies between model systems and patients to translationally deliver relevant biological insights.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

Using human in vitro transcriptome analysis to build trustworthy machine learning models for prediction of animal drug toxicity

During the development of new drugs or compounds there is a requirement for preclinical trials, commonly involving animal tests, to ascertain the safety of the compound prior to human trials. Machine learning techniques could provide an in-silico alternative to animal models for assessing drug toxicity, thus reducing expensive and invasive animal testing during clinical trials, for drugs that are most likely to fail safety tests. Here we present a machine learning model to predict kidney dysfunction, as a proxy for drug induced renal toxicity, in rats. To achieve this, we use inexpensive transcriptomic profiles derived from human cell lines after chemical compound treatment to train our models combined with compound chemical structure information. Genomics data due to its sparse, high-dimensional and noisy nature presents significant challenges in building trustworthy and transparent machine learning models. Here we address these issues by judiciously building feature sets from heterogenous sources and coupling them with measures of model uncertainty achieved through Gaussian Process based Bayesian models. We combine the use of insight into the feature-wise contributions to our predictions with the use of predictive uncertainties recovered from the Gaussian Process to improve the transparency and trustworthiness of the model.

Signatures of cell death and proliferation in perturbation transcriptomics data-from confounding factor to effective prediction

Transcriptional perturbation signatures are valuable data sources for functional genomics. Linking perturbation signatures to screenings opens the possibility to model cellular phenotypes from expression data and to identify efficacious drugs. We linked perturbation transcriptomics data from the LINCS-L1000 project with cell viability information upon genetic (Achilles project) and chemical (CTRP screen) perturbations yielding more than 90 000 signature–viability pairs. An integrated analysis showed that the cell viability signature is a major factor underlying perturbation signatures. The signature is linked to transcription factors regulating cell death, proliferation and division time. We used the cell viability–signature relationship to predict viability from transcriptomics signatures, and identified and validated compounds that induce cell death in tumor cell lines. We showed that cellular toxicity can lead to unexpected similarity of signatures, confounding mechanism of action discovery. Consensus compound signatures predicted cell-specific drug sensitivity, even if the signature is not measured in the same cell line, and outperformed conventional drug-specific features. Our results can help in understanding mechanisms behind cell death and removing confounding factors of transcriptomic perturbation screens. To interactively browse our results and predict cell viability in new gene expression samples, we developed CEVIChE (CEll VIability Calculator from gene Expression; https://saezlab.shinyapps.io/ceviche/).

From expression footprints to causal pathways: contextualizing large signaling networks with CARNIVAL

While gene expression profiling is commonly used to gain an overview of cellular processes, the identification of upstream processes that drive expression changes remains a challenge. To address this issue, we introduce CARNIVAL, a causal network contextualization tool which derives network architectures from gene expression footprints. CARNIVAL (CAusal Reasoning pipeline for Network identification using Integer VALue programming) integrates different sources of prior knowledge including signed and directed protein-protein interactions, transcription factor targets, and pathway signatures. The use of prior knowledge in CARNIVAL enables capturing a broad set of upstream cellular processes and regulators, leading to a higher accuracy when benchmarked against related tools. Implementation as an integer linear programming (ILP) problem guarantees efficient computation. As a case study, we applied CARNIVAL to contextualize signaling networks from gene expression data in IgA nephropathy (IgAN), a condition that can lead to chronic kidney disease. CARNIVAL identified specific signaling pathways and associated mediators dysregulated in IgAN including Wnt and TGF-β, which we subsequently validated experimentally. These results demonstrated how CARNIVAL generates hypotheses on potential upstream alterations that propagate through signaling networks, providing insights into diseases.

Computational translation of genomic responses from experimental model systems to humans

The high failure rate of therapeutics showing promise in mouse models to translate to patients is a pressing challenge in biomedical science. Though retrospective studies have examined the fidelity of mouse models to their respective human conditions, approaches for prospective translation of insights from mouse models to patients remain relatively unexplored. Here, we develop a semi-supervised learning approach for inference of disease-associated human differentially expressed genes and pathways from mouse model experiments. We examined 36 transcriptomic case studies where comparable phenotypes were available for mouse and human inflammatory diseases and assessed multiple computational approaches for inferring human biology from mouse datasets. We found that semi-supervised training of a neural network identified significantly more true human biological associations than interpreting mouse experiments directly. Evaluating the experimental design of mouse experiments where our model was most successful revealed principles of experimental design that may improve translational performance. Our study shows that when prospectively evaluating biological associations in mouse studies, semi-supervised learning approaches, combining mouse and human data for biological inference, provide the most accurate assessment of human in vivo disease processes. Finally, we proffer a delineation of four categories of model system-to-human “Translation Problems” defined by the resolution and coverage of the datasets available for molecular insight translation and suggest that the task of translating insights from model systems to human disease contexts may be better accomplished by a combination of translation-minded experimental design and computational approaches.

Empirical comparison of genomic responses in mouse models and human disease contexts is not sufficient for addressing the challenge of prospective translation from mouse models to human disease contexts. We address this challenge by developing a semi-supervised machine learning approach that combines supervised modeling of mouse datasets with unsupervised modeling of human disease-context datasets to predict human in vivo differentially expressed genes and enriched pathways. Semi-supervised training of a feed forward neural network was the most efficacious model for translating experimentally derived mouse biological associations to the human in vivo disease context. We find that computational generalization of signaling insights substantially improves upon direct generalization of mouse experimental insights and argue that such approaches can facilitate more clinically impactful translation of insights from preclinical studies in model systems to patients.

Phosphoprotein patterns predict trametinib responsiveness and optimal trametinib sensitisation strategies in melanoma

Malignant melanoma is a highly aggressive form of skin cancer responsible for the majority of skin cancer-related deaths. Recent insight into the heterogeneous nature of melanoma suggests more personalised treatments may be necessary to overcome drug resistance and improve patient care. To this end, reliable molecular signatures that can accurately predict treatment responsiveness need to be identified. In this study, we applied multiplex phosphoproteomic profiling across a panel of 24 melanoma cell lines with different disease-relevant mutations, to predict responsiveness to MEK inhibitor trametinib. Supported by multivariate statistical analysis and multidimensional pattern recognition algorithms, the responsiveness of individual cell lines to trametinib could be predicted with high accuracy (83% correct predictions), independent of mutation status. We also successfully employed this approach to case specifically predict whether individual melanoma cell lines could be sensitised to trametinib. Our predictions identified that combining MEK inhibition with selective targeting of c-JUN and/or FAK, using siRNA-based depletion or pharmacological inhibitors, sensitised resistant cell lines and significantly enhanced treatment efficacy. Our study indicates that multiplex proteomic analyses coupled with pattern recognition approaches could assist in personalising trametinib-based treatment decisions in the future.

Dissecting RAF Inhibitor Resistance by Structure-based Modeling Reveals Ways to Overcome Oncogenic RAS Signaling

Clinically used RAF inhibitors are ineffective in RAS mutant tumors because they enhance homo- and heterodimerization of RAF kinases, leading to paradoxical activation of ERK signaling. Overcoming enhanced RAF dimerization and the resulting resistance is a challenge for drug design. Combining multiple inhibitors could be more effective, but it is unclear how the best combinations can be chosen. We built a next-generation mechanistic dynamic model to analyze combinations of structurally different RAF inhibitors, which can efficiently suppress MEK/ERK signaling. This rule-based model of the RAS/ERK pathway integrates thermodynamics and kinetics of drug-protein interactions, structural elements, posttranslational modifications, and cell mutational status as model rules to predict RAF inhibitor combinations for inhibiting ERK activity in oncogenic RAS and/or BRAFV600E backgrounds. Predicted synergistic inhibition of ERK signaling was corroborated by experiments in mutant NRAS, HRAS, and BRAFV600E cells, and inhibition of oncogenic RAS signaling was associated with reduced cell proliferation and colony formation.

The sbv IMPROVER Systems Toxicology Computational Challenge: Identification of Human and Species-Independent Blood Response Markers as Predictors of Smoking Exposure and Cessation Status

Cigarette smoking entails chronic exposure to a mixture of harmful chemicals that trigger molecular changes over time, and is known to increase the risk of developing diseases. Risk assessment in the context of 21^st century toxicology relies on the elucidation of mechanisms of toxicity and the identification of exposure response markers, usually from high-throughput data, using advanced computational methodologies. The sbv IMPROVER Systems Toxicology computational challenge (Fall 2015-Spring 2016) aimed to evaluate whether robust and sparse (≤40 genes) human (sub-challenge 1, SC1) and species-independent (sub-challenge 2, SC2) exposure response markers (so called gene signatures) could be extracted from human and mouse blood transcriptomics data of current (S), former (FS) and never (NS) smoke-exposed subjects as predictors of smoking and cessation status. Best-performing computational methods were identified by scoring anonymized participants' predictions. Worldwide participation resulted in 12 (SC1) and six (SC2) final submissions qualified for scoring. The results showed that blood gene expression data were informative to predict smoking exposure (i.e. discriminating smoker versus never or former smokers) status in human and across species with a high level of accuracy. By contrast, the prediction of cessation status (i.e. distinguishing FS from NS) remained challenging, as reflected by lower classification performances. Participants successfully developed inductive predictive models and extracted human and species-independent gene signatures, including genes with high consensus across teams. Post-challenge analyses highlighted "feature selection" as a key step in the process of building a classifier and confirmed the importance of testing a gene signature in independent cohorts to ensure the generalized applicability of a predictive model at a population-based level. In conclusion, the Systems Toxicology challenge demonstrated the feasibility of extracting a consistent blood-based smoke exposure response gene signature and further stressed the importance of independent and unbiased data and method evaluations to provide confidence in systems toxicology-based scientific conclusions.

Crowd-Sourced Verification of Computational Methods and Data in Systems Toxicology: A Case Study with a Heat-Not-Burn Candidate Modified Risk Tobacco Product

Systems toxicology intends to quantify the effect of toxic molecules in biological systems and unravel their mechanisms of toxicity. The development of advanced computational methods is required for analyzing and integrating high throughput data generated for this purpose as well as for extrapolating predictive toxicological outcomes and risk estimates. To ensure the performance and reliability of the methods and verify conclusions from systems toxicology data analysis, it is important to conduct unbiased evaluations by independent third parties. As a case study, we report here the results of an independent verification of methods and data in systems toxicology by crowdsourcing. The sbv IMPROVER systems toxicology computational challenge aimed to evaluate computational methods for the development of blood-based gene expression signature classification models with the ability to predict smoking exposure status. Participants created/trained models on blood gene expression data sets including smokers/mice exposed to 3R4F (a reference cigarette) or noncurrent smokers/Sham (mice exposed to air). Participants applied their models on unseen data to predict whether subjects classify closer to smoke-exposed or nonsmoke exposed groups. The data sets also included data from subjects that had been exposed to potential modified risk tobacco products (MRTPs) or that had switched to a MRTP after exposure to conventional cigarette smoke. The scoring of anonymized participants' predictions was done using predefined metrics. The top 3 performers' methods predicted class labels with area under the precision recall scores above 0.9. Furthermore, although various computational approaches were used, the crowd's results confirmed our own data analysis outcomes with regards to the classification of MRTP-related samples. Mice exposed directly to a MRTP were classified closer to the Sham group. After switching to a MRTP, the confidence that subjects belonged to the smoke-exposed group decreased significantly. Smoking exposure gene signatures that contributed to the group separation included a core set of genes highly consistent across teams such as AHRR, LRRN3, SASH1, and P2RY6. In conclusion, crowdsourcing constitutes a pertinent approach, in complement to the classical peer review process, to independently and unbiasedly verify computational methods and data for risk assessment using systems toxicology.

Hypothiocyanite produced by human and rat respiratory epithelial cells inactivates extracellular H1N2 influenza A virus

OBJECTIVE AND DESIGN

Our aim was to study whether an extracellular, oxidative antimicrobial mechanism inherent to tracheal epithelial cells is capable of inactivating influenza H1N2 virus.

MATERIAL OR SUBJECTS

Epithelial cells were isolated from tracheas of male Sprague-Dawley rats. Both primary human and rat tracheobronchial epithelial cells were differentiated in air-liquid interface cultures.

TREATMENT

A/swine/Illinois/02860/09 (swH1N2) influenza A virions were added to the apical side of airway cells for 1 h in the presence or absence of lactoperoxidase or thiocyanate.

METHODS

Characterization of rat epithelial cells (morphology, Duox expression) occurred via western blotting, PCR, hydrogen peroxide production measurement and histology. The number of viable virions was determined by plaque assays. Statistical difference of the results was analyzed by ANOVA and Tukey's test.

RESULTS

Our data show that rat tracheobronchial epithelial cells develop a differentiated, polarized monolayer with high transepithelial electrical resistance, mucin production and expression of dual oxidases. Influenza A virions are inactivated by human and rat epithelial cells via a dual oxidase-, lactoperoxidase- and thiocyanate-dependent mechanism.

CONCLUSIONS

Differentiated air-liquid interface cultures of rat tracheal epithelial cells provide a novel model to study airway epithelium-influenza interactions. The dual oxidase/lactoperoxidase/thiocyanate extracellular oxidative system producing hypothiocyanite is a fast and potent anti-influenza mechanism inactivating H1N2 viruses prior to infection of the epithelium.

Animal models of Duchenne muscular dystrophy: from basic mechanisms to gene therapy

Duchenne muscular dystrophy (DMD) is a progressive muscle-wasting disorder. It is caused by loss-of-function mutations in the dystrophin gene. Currently, there is no cure. A highly promising therapeutic strategy is to replace or repair the defective dystrophin gene by gene therapy. Numerous animal models of DMD have been developed over the last 30 years, ranging from invertebrate to large mammalian models. mdx mice are the most commonly employed models in DMD research and have been used to lay the groundwork for DMD gene therapy. After ~30 years of development, the field has reached the stage at which the results in mdx mice can be validated and scaled-up in symptomatic large animals. The canine DMD (cDMD) model will be excellent for these studies. In this article, we review the animal models for DMD, the pros and cons of each model system, and the history and progress of preclinical DMD gene therapy research in the animal models. We also discuss the current and emerging challenges in this field and ways to address these challenges using animal models, in particular cDMD dogs.

Trans-species learning of cellular signaling systems with bimodal deep belief networks

MOTIVATION

Model organisms play critical roles in biomedical research of human diseases and drug development. An imperative task is to translate information/knowledge acquired from model organisms to humans. In this study, we address a trans-species learning problem: predicting human cell responses to diverse stimuli, based on the responses of rat cells treated with the same stimuli.

RESULTS

We hypothesized that rat and human cells share a common signal-encoding mechanism but employ different proteins to transmit signals, and we developed a bimodal deep belief network and a semi-restricted bimodal deep belief network to represent the common encoding mechanism and perform trans-species learning. These 'deep learning' models include hierarchically organized latent variables capable of capturing the statistical structures in the observed proteomic data in a distributed fashion. The results show that the models significantly outperform two current state-of-the-art classification algorithms. Our study demonstrated the potential of using deep hierarchical models to simulate cellular signaling systems.

AVAILABILITY AND IMPLEMENTATION

The software is available at the following URL: http://pubreview.dbmi.pitt.edu/TransSpeciesDeepLearning/. The data are available through SBV IMPROVER website, https://www.sbvimprover.com/challenge-2/overview, upon publication of the report by the organizers.

CONTACT

xinghua@pitt.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Network-Based Analysis of Nutraceuticals in Human Hepatocellular Carcinomas Reveals Mechanisms of Chemopreventive Action

Chronic inflammation is associated with the development of human hepatocellular carcinoma (HCC), an essentially incurable cancer. Anti-inflammatory nutraceuticals have emerged as promising candidates against HCC, yet the mechanisms through which they influence the cell signaling machinery to impose phenotypic changes remain unresolved. Herein we implemented a systems biology approach in HCC cells, based on the integration of cytokine release and phospoproteomic data from high-throughput xMAP Luminex assays to elucidate the action mode of prominent nutraceuticals in terms of topology alterations of HCC-specific signaling networks. An optimization algorithm based on SigNetTrainer, an Integer Linear Programming formulation, was applied to construct networks linking signal transduction to cytokine secretion by combining prior knowledge of protein connectivity with proteomic data. Our analysis identified the most probable target phosphoproteins of interrogated compounds and predicted translational control as a new mechanism underlying their anticytokine action. Induced alterations corroborated with inhibition of HCC-driven angiogenesis and metastasis.

Network reconstruction based on proteomic data and prior knowledge of protein connectivity using graph theory

Modeling of signal transduction pathways is instrumental for understanding cells’ function. People have been tackling modeling of signaling pathways in order to accurately represent the signaling events inside cells’ biochemical microenvironment in a way meaningful for scientists in a biological field. In this article, we propose a method to interrogate such pathways in order to produce cell-specific signaling models. We integrate available prior knowledge of protein connectivity, in a form of a Prior Knowledge Network (PKN) with phosphoproteomic data to construct predictive models of the protein connectivity of the interrogated cell type. Several computational methodologies focusing on pathways’ logic modeling using optimization formulations or machine learning algorithms have been published on this front over the past few years. Here, we introduce a light and fast approach that uses a breadth-first traversal of the graph to identify the shortest pathways and score proteins in the PKN, fitting the dependencies extracted from the experimental design. The pathways are then combined through a heuristic formulation to produce a final topology handling inconsistencies between the PKN and the experimental scenarios. Our results show that the algorithm we developed is efficient and accurate for the construction of medium and large scale signaling networks. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGF/TNFA stimulation against made up experimental data. To avoid the possibility of erroneous predictions, we performed a cross-validation analysis. Finally, we validate that the introduced approach generates predictive topologies, comparable to the ILP formulation. Overall, an efficient approach based on graph theory is presented herein to interrogate protein–protein interaction networks and to provide meaningful biological insights.

Modeling signal transduction from protein phosphorylation to gene expression

BACKGROUND

Signaling networks are of great importance for us to understand the cell’s regulatory mechanism. The rise of large-scale genomic and proteomic data, and prior biological knowledge has paved the way for the reconstruction and discovery of novel signaling pathways in a data-driven manner. In this study, we investigate computational methods that integrate proteomics and transcriptomic data to identify signaling pathways transmitting signals in response to specific stimuli. Such methods can be applied to cancer genomic data to infer perturbed signaling pathways.

METHOD

We proposed a novel Bayesian Network (BN) framework to integrate transcriptomic data with proteomic data reflecting protein phosphorylation states for the purpose of identifying the pathways transmitting the signal of diverse stimuli in rat and human cells. We represented the proteins and genes as nodes in a BN in which edges reflect the regulatory relationship between signaling proteins. We designed an efficient inference algorithm that incorporated the prior knowledge of pathways and searched for a network structure in a data-driven manner.

RESULTS

We applied our method to infer rat and human specific networks given gene expression and proteomic datasets. We were able to effectively identify sparse signaling networks that modeled the observed transcriptomic and proteomic data. Our methods were able to identify distinct signaling pathways for rat and human cells in a data-driven manner, based on the facts that rat and human cells exhibited distinct transcriptomic and proteomics responses to a common set of stimuli. Our model performed well in the SBV IMPROVER challenge in comparison to other models addressing the same task. The capability of inferring signaling pathways in a data-driven fashion may contribute to cancer research by identifying distinct aberrations in signaling pathways underlying heterogeneous cancers subtypes.

Understanding the limits of animal models as predictors of human biology: lessons learned from the sbv IMPROVER Species Translation Challenge

Motivation: Inferring how humans respond to external cues such as drugs, chemicals, viruses or hormones is an essential question in biomedicine. Very often, however, this question cannot be addressed because it is not possible to perform experiments in humans. A reasonable alternative consists of generating responses in animal models and ‘translating’ those results to humans. The limitations of such translation, however, are far from clear, and systematic assessments of its actual potential are urgently needed. sbv IMPROVER (systems biology verification for Industrial Methodology for PROcess VErification in Research) was designed as a series of challenges to address translatability between humans and rodents. This collaborative crowd-sourcing initiative invited scientists from around the world to apply their own computational methodologies on a multilayer systems biology dataset composed of phosphoproteomics, transcriptomics and cytokine data derived from normal human and rat bronchial epithelial cells exposed in parallel to 52 different stimuli under identical conditions. Our aim was to understand the limits of species-to-species translatability at different levels of biological organization: signaling, transcriptional and release of secreted factors (such as cytokines). Participating teams submitted 49 different solutions across the sub-challenges, two-thirds of which were statistically significantly better than random. Additionally, similar computational methods were found to range widely in their performance within the same challenge, and no single method emerged as a clear winner across all sub-challenges. Finally, computational methods were able to effectively translate some specific stimuli and biological processes in the lung epithelial system, such as DNA synthesis, cytoskeleton and extracellular matrix, translation, immune/inflammation and growth factor/proliferation pathways, better than the expected response similarity between species.

Contact:pmeyerr@us.ibm.com or Julia.Hoeng@pmi.com

Supplementary information:Supplementary data are available at Bioinformatics online.

Inter-species inference of gene set enrichment in lung epithelial cells from proteomic and large transcriptomic datasets

Motivation: Translating findings in rodent models to human models has been a cornerstone of modern biology and drug development. However, in many cases, a naive ‘extrapolation’ between the two species has not succeeded. As a result, clinical trials of new drugs sometimes fail even after considerable success in the mouse or rat stage of development. In addition to in vitro studies, inter-species translation requires analytical tools that can predict the enriched gene sets in human cells under various stimuli from corresponding measurements in animals. Such tools can improve our understanding of the underlying biology and optimize the allocation of resources for drug development.

Results: We developed an algorithm to predict differential gene set enrichment as part of the sbv IMPROVER (systems biology verification in Industrial Methodology for Process Verification in Research) Species Translation Challenge, which focused on phosphoproteomic and transcriptomic measurements of normal human bronchial epithelial (NHBE) primary cells under various stimuli and corresponding measurements in rat (NRBE) primary cells. We find that gene sets exhibit a higher inter-species correlation compared with individual genes, and are potentially more suited for direct prediction. Furthermore, in contrast to a similar cross-species response in protein phosphorylation states 5 and 25 min after exposure to stimuli, gene set enrichment 6 h after exposure is significantly different in NHBE cells compared with NRBE cells. In spite of this difference, we were able to develop a robust algorithm to predict gene set activation in NHBE with high accuracy using simple analytical methods.

Availability and implementation: Implementation of all algorithms is available as source code (in Matlab) at http://bhanot.biomaps.rutgers.edu/wiki/codes_SC3_Predicting_GeneSets.zip, along with the relevant data used in the analysis. Gene sets, gene expression and protein phosphorylation data are available on request.

Contact:hormoz@kitp.ucsb.edu

Proteomics for systems toxicology

Current toxicology studies frequently lack measurements at molecular resolution to enable a more mechanism-based and predictive toxicological assessment. Recently, a systems toxicology assessment framework has been proposed, which combines conventional toxicological assessment strategies with system-wide measurement methods and computational analysis approaches from the field of systems biology. Proteomic measurements are an integral component of this integrative strategy because protein alterations closely mirror biological effects, such as biological stress responses or global tissue alterations. Here, we provide an overview of the technical foundations and highlight select applications of proteomics for systems toxicology studies. With a focus on mass spectrometry-based proteomics, we summarize the experimental methods for quantitative proteomics and describe the computational approaches used to derive biological/mechanistic insights from these datasets. To illustrate how proteomics has been successfully employed to address mechanistic questions in toxicology, we summarized several case studies. Overall, we provide the technical and conceptual foundation for the integration of proteomic measurements in a more comprehensive systems toxicology assessment framework. We conclude that, owing to the critical importance of protein-level measurements and recent technological advances, proteomics will be an integral part of integrative systems toxicology approaches in the future.