Cancer is a complex collection of diseases that are to some degree unique to each patient. Precision oncology aims to identify the best drug treatment regime using molecular data on tumor samples. While omics-level data is becoming more widely available for tumor specimens, the datasets upon which computational learning methods can be trained vary in coverage from sample to sample and from data type to data type. Methods that can 'connect the dots' to leverage more of the information provided by these studies could offer major advantages for maximizing predictive potential. We introduce a multi-view machinelearning strategy called PLATYPUS that builds 'views' from multiple data sources that are all used as features for predicting patient outcomes. We show that a learning strategy that finds agreement across the views on unlabeled data increases the performance of the learning methods over any single view. We illustrate the power of the approach by deriving signatures for drug sensitivity in a large cancer cell line database. Code and additional information are available from the PLATYPUS website https://sysbiowiki.soe.ucsc.edu/platypus.

Abstract Combining heterogeneous sources of data is essential for accurate prediction of protein function. The task is complicated by the fact that while sequence-based features can be readily compared across species, most other data are species-specific. In this paper, we present a multi-view extension to GOstruct, a structured-output framework for function annotation of proteins. The extended framework can learn from disparate data sources, with each data source provided to the framework in the form of a kernel. Our empirical results demonstrate that the multi-view framework is able to utilize all available information, yielding better performance than sequence-based models trained across species and models trained from collections of data within a given species. This version of GOstruct participated in the recent Critical Assessment of Functional Annotations (CAFA) challenge; since then we have significantly improved the natural language processing component of the method, which now provides performance that is on par with that provided by sequence information. The GOstruct framework is available for download at http://strut.sourceforge.net.

The purpose of this study was to investigate changes in the lipidome of patients with sepsis to identify signaling lipids associated with poor outcomes that could be linked to future therapies. Adult patients with sepsis were enrolled within 24h of sepsis recognition. Patients meeting Sepsis-3 criteria were enrolled from the emergency department or intensive care unit and blood samples were obtained. Clinical data were collected and outcomes of rapid recovery, chronic critical illness (CCI), or early death were adjudicated by clinicians. Lipidomic analysis was performed on two platforms, the Sciex™ 5500 device to perform a lipidomic screen of 1450 lipid species and a targeted signaling lipid panel using liquid-chromatography tandem mass spectrometry. For the lipidomic screen, there were 274 patients with sepsis: 192 with rapid recovery, 47 with CCI, and 35 with early deaths. CCI and early death patients were grouped together for analysis. Fatty acid (FA) 12:0 was decreased in CCI/early death, whereas FA 17:0 and 20:1 were elevated in CCI/early death, compared to rapid recovery patients. For the signaling lipid panel analysis, there were 262 patients with sepsis: 189 with rapid recovery, 45 with CCI, and 28 with early death. Pro-inflammatory signaling lipids from ω-6 poly-unsaturated fatty acids (PUFAs), including 15-hydroxyeicosatetraenoic (HETE), 12-HETE, and 11-HETE (oxidation products of arachidonic acid [AA]) were elevated in CCI/early death patients compared to rapid recovery. The pro-resolving lipid mediator from ω-3 PUFAs, 14(S)-hydroxy docosahexaenoic acid (14S-HDHA), was also elevated in CCI/early death compared to rapid recovery. Signaling lipids of the AA pathway were elevated in poor-outcome patients with sepsis and may serve as targets for future therapies.

Background

Patient stratification to identify subtypes with different disease manifestations, severity, and expected survival time is a critical task in cancer diagnosis and treatment. While stratification approaches using various biomarkers (including high-throughput gene expression measurements) for patient-to-patient comparisons have been successful in elucidating previously unseen subtypes, there remains an untapped potential of incorporating various genotypic and phenotypic data to discover novel or improved groupings.

Methods

Here, we present HOCUS, a unified analytical framework for patient stratification that uses a community detection technique to extract subtypes out of sparse patient measurements. HOCUS constructs a patient-to-patient network from similarities in the data and iteratively groups and reconstructs the network into higher order clusters. We investigate the merits of using higher-order correlations to cluster samples of cancer patients in terms of their associations with survival outcomes.

Results

In an initial test of the method, the approach identifies cancer subtypes in mutation data of glioblastoma, ovarian, breast, prostate, and bladder cancers. In several cases, HOCUS provides an improvement over using the molecular features directly to compare samples. Application of HOCUS to glioblastoma images reveals a size and location classification of tumors that improves over human expert-based stratification.

Conclusions

Subtypes based on higher order features can reveal comparable or distinct groupings. The distinct solutions can provide biologically- and treatment-relevant solutions that are just as significant as solutions based on the original data.

Evidence from numerous cancers suggests that increased aggressiveness is accompanied by up-regulation of signaling pathways and acquisition of properties common to stem cells. It is unclear if different subtypes of late-stage cancer vary in stemness properties and whether or not these subtypes are transcriptionally similar to normal tissue stem cells. We report a gene signature specific for human prostate basal cells that is differentially enriched in various phenotypes of late-stage metastatic prostate cancer. We FACS-purified and transcriptionally profiled basal and luminal epithelial populations from the benign and cancerous regions of primary human prostates. High-throughput RNA sequencing showed the basal population to be defined by genes associated with stem cell signaling programs and invasiveness. Application of a 91-gene basal signature to gene expression datasets from patients with organ-confined or hormone-refractory metastatic prostate cancer revealed that metastatic small cell neuroendocrine carcinoma was molecularly more stem-like than either metastatic adenocarcinoma or organ-confined adenocarcinoma. Bioinformatic analysis of the basal cell and two human small cell gene signatures identified a set of E2F target genes common between prostate small cell neuroendocrine carcinoma and primary prostate basal cells. Taken together, our data suggest that aggressive prostate cancer shares a conserved transcriptional program with normal adult prostate basal stem cells.

We introduce a novel method called Prophetic Granger Causality (PGC) for inferring gene regulatory networks (GRNs) from protein-level time series data. The method uses an L1-penalized regression adaptation of Granger Causality to model protein levels as a function of time, stimuli, and other perturbations. When combined with a data-independent network prior, the framework outperformed all other methods submitted to the HPN-DREAM 8 breast cancer network inference challenge. Our investigations reveal that PGC provides complementary information to other approaches, raising the performance of ensemble learners, while on its own achieves moderate performance. Thus, PGC serves as a valuable new tool in the bioinformatics toolkit for analyzing temporal datasets. We investigate the general and cell-specific interactions predicted by our method and find several novel interactions, demonstrating the utility of the approach in charting new tumor wiring.

Recent genomic analyses have provided substantial evidence for past periods of gene flow from polar bears (Ursus maritimus) into Alaskan brown bears (Ursus arctos), with some analyses suggesting a link between climate change and genomic introgression. However, because it has mainly been possible to sample bears from the present day, the timing, frequency, and evolutionary significance of this admixture remains unknown. Here, we analyze genomic DNA from three additional and geographically distinct brown bear populations, including two that lived temporally close to the peak of the last ice age. We find evidence of admixture in all three populations, suggesting that admixture between these species has been common in their recent evolutionary history. In addition, analyses of ten fossil bears from the now-extinct Irish population indicate that admixture peaked during the last ice age, whereas brown bear and polar bear ranges overlapped. Following this peak, the proportion of polar bear ancestry in Irish brown bears declined rapidly until their extinction. Our results support a model in which ice age climate change created geographically widespread conditions conducive to admixture between polar bears and brown bears, as is again occurring today. We postulate that this model will be informative for many admixing species pairs impacted by climate change. Our results highlight the power of paleogenomics to reveal patterns of evolutionary change that are otherwise masked in contemporary data.

Vast amounts of molecular data are being collected on tumor samples, which provide unique opportunities for discovering trends within and between cancer subtypes. Such cross-cancer analyses require computational methods that enable intuitive and interactive browsing of thousands of samples based on their molecular similarity. We created a portal called TumorMap to assist in exploration and statistical interrogation of high-dimensional complex "omics" data in an interactive and easily interpretable way. In the TumorMap, samples are arranged on a hexagonal grid based on their similarity to one another in the original genomic space and are rendered with Google's Map technology. While the important feature of this public portal is the ability for the users to build maps from their own data, we pre-built genomic maps from several previously published projects. We demonstrate the utility of this portal by presenting results obtained from The Cancer Genome Atlas project data. Cancer Res; 77(21); e111-4. ©2017 AACR.

This is a study of lipid metabolic gene expression patterns to discover precision medicine for sepsis.

Objectives

Sepsis patients experience poor outcomes including chronic critical illness (CCI) or early death (within 14 d). We investigated lipid metabolic gene expression differences by outcome to discover therapeutic targets.

Design setting and particitpants

Secondary analysis of samples from prospectively enrolled sepsis patients (first 24 hr) and a zebrafish endotoxemia model for drug discovery. Patients were enrolled from the emergency department or ICU at an urban teaching hospital. Enrollment samples from sepsis patients were analyzed. Clinical data and cholesterol levels were recorded. Leukocytes were processed for RNA sequencing and reverse transcriptase polymerase chain reaction. A lipopolysaccharide zebrafish endotoxemia model was used for confirmation of human transcriptomic findings and drug discovery.

Main outcomes and measures

The derivation cohort included 96 patients and controls (12 early death, 13 CCI, 51 rapid recovery, and 20 controls) and the validation cohort had 52 patients (6 early death, 8 CCI, and 38 rapid recovery).

Results

The cholesterol metabolism gene 7-dehydrocholesterol reductase (DHCR7) was significantly up-regulated in both derivation and validation cohorts in poor outcome sepsis compared with rapid recovery patients and in 90-day nonsurvivors (validation only) and validated using RT-qPCR analysis. Our zebrafish sepsis model showed up-regulation of dhcr7 and several of the same lipid genes up-regulated in poor outcome human sepsis (dhcr24, sqlea, cyp51, msmo1, and ldlra) compared with controls. We then tested six lipid-based drugs in the zebrafish endotoxemia model. Of these, only the Dhcr7 inhibitor AY9944 completely rescued zebrafish from lipopolysaccharide death in a model with 100% lethality.

Conclusions

DHCR7, an important cholesterol metabolism gene, was up-regulated in poor outcome sepsis patients warranting external validation. This pathway may serve as a potential therapeutic target to improve sepsis outcomes.