Mapping the effects of drugs on the immune system

Pancreas Cancer Precision Treatment Using Avatar Mice from a Bioinformatics Perspective

Pancreatic ductal adenocarcinoma (PDAC) is a leading cause of cancer-related death among solid malignancies. Unfortunately, PDAC lethality has not substantially decreased over the past 20 years. This aggressiveness is related to the genomic complexity and heterogeneity of PDAC, but also to the absence of an effective screening for the detection of early-stage tumors and a lack of efficient therapeutic options. Therefore, there is an urgent need to improve the arsenal of anti-PDAC drugs for an effective treatment of these patients. Patient-derived xenograft (PDX) mouse models represent a promising strategy to personalize PDAC treatment, offering a bench testing of candidate treatments and helping to select empirical treatments in PDAC patients with no therapeutic targets. Moreover, bioinformatics-based approaches have the potential to offer systematic insights into PDAC etiology predicting putatively actionable tumor-specific genomic alterations, identifying novel biomarkers and generating disease-associated gene expression signatures. This review focuses on recent efforts to individualize PDAC treatments using PDX models. Additionally, we discuss the current understanding of the PDAC genomic landscape and the putative druggable targets derived from mutational studies. PDAC molecular subclassifications and gene expression profiling studies are reviewed as well. Finally, latest bioinformatics methodologies based on somatic variant detection and prioritization, in silico drug response prediction, and drug repositioning to improve the treatment of advanced PDAC tumors are also covered.

Comparing structural and transcriptional drug networks reveals signatures of drug activity and toxicity in transcriptional responses

We performed an integrated analysis of drug chemical structures and drug-induced transcriptional responses. We demonstrated that a network representing three-dimensional structural similarities among 5452 compounds can be used to automatically group together drugs with similar scaffolds, physicochemical parameters and mode-of-action. We compared the structural network to a network representing transcriptional similarities among a subset of 1309 drugs for which transcriptional response were available in the Connectivity Map data set. Analysis of structurally similar, but transcriptionally different drugs sharing the same MOA enabled us to detect and remove weak and noisy transcriptional responses, greatly enhancing the reliability of transcription-based approaches to drug discovery and drug repositioning. Cardiac glycosides exhibited the strongest transcriptional responses with a significant induction of pathways related to epigenetic regulation, which suggests an epigenetic mechanism of action for these drugs. Drug classes with the weakest transcriptional responses tended to induce expression of cytochrome P450 enzymes, hinting at drug-induced drug resistance. Analysis of transcriptionally similar, but structurally different drugs with unrelated MOA, led us to the identification of a 'toxic' transcriptional signature indicative of lysosomal stress (lysosomotropism) and lipid accumulation (phospholipidosis) partially masking the target-specific transcriptional effects of these drugs. We found that this transcriptional signature is shared by 258 compounds and it is associated to the activation of the transcription factor TFEB, a master regulator of lysosomal biogenesis and autophagy. Finally, we built a predictive Random Forest model of these 258 compounds based on 128 physicochemical parameters, which should help in the early identification of potentially toxic drug candidates.

Transcriptional responses to drug treatment can reveal mechanism of action and off-target effects thus enabling drug repositioning, but only if measured in the appropriate cells at clinically relevant concentrations. A team led by Diego di Bernardo and Diego Medina generated a network representing structural similarities among compounds to automatically group together drugs with similar scaffolds and MOA. By comparing the structural drug network with a transcriptional drug network based on similarities in transcriptional response, the team observed broad differences between the two. This observation led to the identification of a transcriptional signature related lysosomal stress and phospholipidosis, and a physicochemical model to identify such compounds. These results provide general guidelines to prevent erroneous conclusion when using transcriptional responses of small molecules for drug discovery and drug repositioning

Reversal of cancer gene expression correlates with drug efficacy and reveals therapeutic targets

The decreasing cost of genomic technologies has enabled the molecular characterization of large-scale clinical disease samples and of molecular changes upon drug treatment in various disease models. Exploring methods to relate diseases to potentially efficacious drugs through various molecular features is critically important in the discovery of new therapeutics. Here we show that the potency of a drug to reverse cancer-associated gene expression changes positively correlates with that drug's efficacy in preclinical models of breast, liver and colon cancers. Using a systems-based approach, we predict four compounds showing high potency to reverse gene expression in liver cancer and validate that all four compounds are effective in five liver cancer cell lines. The in vivo efficacy of pyrvinium pamoate is further confirmed in a subcutaneous xenograft model. In conclusion, this systems-based approach may be complementary to the traditional target-based approach in connecting diseases to potentially efficacious drugs.

Novel Therapeutics Identification for Fibrosis in Renal Allograft Using Integrative Informatics Approach

Chronic allograft damage, defined by interstitial fibrosis and tubular atrophy (IF/TA), is a leading cause of allograft failure. Few effective therapeutic options are available to prevent the progression of IF/TA. We applied a meta-analysis approach on IF/TA molecular datasets in Gene Expression Omnibus to identify a robust 85-gene signature, which was used for computational drug repurposing analysis. Among the top ranked compounds predicted to be therapeutic for IF/TA were azathioprine, a drug to prevent acute rejection in renal transplantation, and kaempferol and esculetin, two drugs not previously described to have efficacy for IF/TA. We experimentally validated the anti-fibrosis effects of kaempferol and esculetin using renal tubular cells in vitro and in vivo in a mouse Unilateral Ureteric Obstruction (UUO) model. Kaempferol significantly attenuated TGF-β1-mediated profibrotic pathways in vitro and in vivo, while esculetin significantly inhibited Wnt/β-catenin pathway in vitro and in vivo. Histology confirmed significantly abrogated fibrosis by kaempferol and esculetin in vivo. We developed an integrative computational framework to identify kaempferol and esculetin as putatively novel therapies for IF/TA and provided experimental evidence for their therapeutic activities in vitro and in vivo using preclinical models. The findings suggest that both drugs might serve as therapeutic options for IF/TA.

Transcriptional Profiling of Egg Allergy and Relationship to Disease Phenotype

Background

Egg allergy is one of the most common food allergies of childhood. There is a lack of information on the immunologic basis of egg allergy beyond the role of IgE.

Objective

To use transcriptional profiling as a novel approach to uncover immunologic processes associated with different phenotypes of egg allergy.

Methods

Peripheral blood mononuclear cells (PBMCs) were obtained from egg-allergic children who were defined as reactive (BER) or tolerant (BET) to baked egg, and from food allergic controls (AC) who were egg non-allergic. PBMCs were stimulated with egg white protein. Gene transcription was measured by microarray after 24 h, and cytokine secretion by multiplex assay after 5 days.

Results

The transcriptional response of PBMCs to egg protein differed between BER and BET versus AC subjects. Compared to the AC group, the BER group displayed increased expression of genes associated with allergic inflammation as well as corresponding increased secretion of IL-5, IL-9 and TNF-α. A similar pattern was observed for the BET group. Further similarities in gene expression patterns between BER and BET groups, as well as some important differences, were revealed using a novel Immune Annotation resource developed for this project. This approach identified several novel processes not previously associated with egg allergy, including positive associations with TLR4-stimulated myeloid cells and activated NK cells, and negative associations with an induced Treg signature. Further pathway analysis of differentially expressed genes comparing BER to BET subjects showed significant enrichment of IFN-α and IFN-γ response genes, as well as genes associated with virally-infected DCs.

Conclusions

Transcriptional profiling identified several novel pathways and processes that differed when comparing the response to egg allergen in BET, BER, and AC groups. We conclude that this approach is a useful hypothesis-generating mechanism to identify novel immune processes associated with allergy and tolerance to forms of egg.

In Search of 'Birth Month Genes': Using Existing Data Repositories to Locate Genes Underlying Birth Month-Disease Relationships

Prenatal and perinatal exposures vary seasonally (e.g., sunlight, allergens) and many diseases are linked with variance in exposure. Epidemiologists often measure these changes using birth month as a proxy for seasonal variance. Likewise, Genome-Wide Association Studies have associated or implicated these same diseases with many genes. Both disparate data types (epidemiological and genetic) can provide key insights into the underlying disease biology. We developed an algorithm that links 1) epidemiological data from birth month studies with 2) genetic data from published gene-disease association studies. Our framework uses existing data repositories - PubMed, DisGeNET and Gene Ontology - to produce a bipartite network that connects enriched seasonally varying biofactorss with birth month dependent diseases (BMDDs) through their overlapping developmental gene sets. As a proof-of-concept, we investigate 7 known BMDDs and highlight three important biological networks revealed by our algorithm and explore some interesting genetic mechanisms potentially responsible for the seasonal contribution to BMDDs.

The roles of inflammation and immune mechanisms in Alzheimer's disease

The Alzheimer's Association's Research roundtable met in April 2015 to explore the role of neuroinflammatory mechanisms in the progression of Alzheimer's disease (AD). The ability of innate immune cells, particularly microglia and astrocytes, to mediate neuroinflammation in AD has been implicated as a significant contributor to disease pathogenesis. Adaptive immunity, which plays an important role in responding to injury and some diseases of the central nervous system, may contribute to neuroinflammation in AD as well. Communication between the central and peripheral immune systems may also be important in AD. An increased understanding of the physiology of the innate immune system may aid the identification of new therapeutic targets or mechanisms. The development of predictive animal models and translatable neuroinflammation biomarkers for AD would also facilitate the advancement of novel treatments for innate immunity. Important challenges impeding the advancement of new therapeutic agents and strategies to overcome them were discussed.

Quantitative Systems Pharmacology: A Framework for Context

Quantitative Systems Pharmacology (QSP) is receiving increased attention. As the momentum builds and the expectations grow it is important to (re)assess and formalize the basic concepts and approaches. In this short review, I argue that QSP, in addition to enabling the rational integration of data and development of complex models, maybe more importantly, provides the foundations for developing an integrated framework for the assessment of drugs and their impact on disease within a broader context expanding the envelope to account in great detail for physiology, environment and prior history. I articulate some of the critical enablers, major obstacles and exciting opportunities manifesting themselves along the way. Charting such overarching themes will enable practitioners to identify major and defining factors as the field progressively moves towards personalized and precision health care delivery.

Exploiting Gene-Expression Deconvolution to Probe the Genetics of the Immune System

Sequence variation can affect the physiological state of the immune system. Major experimental efforts targeted at understanding the genetic control of the abundance of immune cell subpopulations. However, these studies are typically focused on a limited number of immune cell types, mainly due to the use of relatively low throughput cell-sorting technologies. Here we present an algorithm that can reveal the genetic basis of inter-individual variation in the abundance of immune cell types using only gene expression and genotyping measurements as input. Our algorithm predicts the abundance of immune cell subpopulations based on the RNA levels of informative marker genes within a complex tissue, and then provides the genetic control on these predicted immune traits as output. A key feature of the approach is the integration of predictions from various sets of marker genes and refinement of these sets to avoid spurious signals. Our evaluation of both synthetic and real biological data shows the significant benefits of the new approach. Our method, VoCAL, is implemented in the freely available R package ComICS.

Quantitative trait locus (QTL) studies have identified a plethora of genetic variants that lead to inter-individual variation in the abundance of immune cell subpopulations, both in normal and disease states. Cell sorting is an effective method of monitoring immune cell type quantities; however, owing to the large number of possible immune cell subsets, it can be difficult to apply this method to each cell type over multiple individuals. Recent QTL studies dealt with this difficulty by focusing on an a priori selection of one or a few cell subsets. Here we introduce VoCAL, a deconvolution-based method that utilizes transcriptome data to infer the quantities of immune cell types, and then uses these quantitative traits to uncover the underlying DNA loci. Our results in synthetic data and lung cohorts show that the VoCAL method outperforms other alternatives in revealing the genetic basis of immune physiology.