Search

Background: High-grade serous carcinoma (HGSC) gene expression subtypes are associated with differential survival. We characterized HGSC gene expression in Black individuals and considered whether gene expression differences by self-identified race may contribute to poorer HGSC survival among Black versus White individuals., Methods: We included newly generated RNA sequencing data from Black and White individuals and array-based genotyping data from four existing studies of White and Japanese individuals. We used K-means clustering, a method with no predefined number of clusters or dataset-specific features, to assign subtypes. Cluster- and dataset-specific gene expression patterns were summarized by moderated t-scores. We compared cluster-specific gene expression patterns across datasets by calculating the correlation between the summarized vectors of moderated t-scores. After mapping to The Cancer Genome Atlas-derived HGSC subtypes, we used Cox proportional hazards models to estimate subtype-specific survival by dataset., Results: Cluster-specific gene expression was similar across gene expression platforms and racial groups. Comparing the Black population with the White and Japanese populations, the immunoreactive subtype was more common (39% vs. 23%-28%) and the differentiated subtype was less common (7% vs. 22%-31%). Patterns of subtype-specific survival were similar between the Black and White populations with RNA sequencing data; compared with mesenchymal cases, the risk of death was similar for proliferative and differentiated cases and suggestively lower for immunoreactive cases [Black population HR = 0.79 (0.55, 1.13); White population HR = 0.86 (0.62, 1.19)]., Conclusions: Although the prevalence of HGSC subtypes varied by race, subtype-specific survival was similar., Impact: HGSC subtypes can be consistently assigned across platforms and self-identified racial groups., (©2024 American Association for Cancer Research.)

Advances in high-throughput microscopy have enabled the rapid acquisition of large numbers of high-content microscopy images. Whether by deep learning or classical algorithms, image analysis pipelines then produce single-cell features. To process these single-cells for downstream applications, we present Pycytominer, a user-friendly, open-source python package that implements the bioinformatics steps, known as "image-based profiling". We demonstrate Pycytominer's usefulness in a machine learning project to predict nuisance compounds that cause undesirable cell injuries., Competing Interests: Disclosure of competing interests The authors declare they have no competing interests relevant to this work.

The morphology of cells is dynamic and mediated by genetic and environmental factors. Characterizing how genetic variation impacts cell morphology can provide an important link between disease association and cellular function. Here, we combine genomic sequencing and high-content imaging approaches on iPSCs from 297 unique donors to investigate the relationship between genetic variants and cellular morphology to map what we term cell morphological quantitative trait loci (cmQTLs). We identify novel associations between rare protein altering variants in WASF2, TSPAN15, and PRLR with several morphological traits related to cell shape, nucleic granularity, and mitochondrial distribution. Knockdown of these genes by CRISPRi confirms their role in cell morphology. Analysis of common variants yields one significant association and nominate over 300 variants with suggestive evidence (P < 10 -6 ) of association with one or more morphology traits. We then use these data to make predictions about sample size requirements for increasing discovery in cellular genetic studies. We conclude that, similar to molecular phenotypes, morphological profiling can yield insight about the function of genes and variants., (© 2024. The Author(s).)

The field of high content imaging has steadily evolved and expanded substantially across many industry and academic research institutions since it was first described in the early 1990's. High content imaging refers to the automated acquisition and analysis of microscopic images from a variety of biological sample types. Integration of high content imaging microscopes with multiwell plate handling robotics enables high content imaging to be performed at scale and support medium- to high-throughput screening of pharmacological, genetic and diverse environmental perturbations upon complex biological systems ranging from 2D cell cultures to 3D tissue organoids to small model organisms. In this perspective article the authors provide a collective view on the following key discussion points relevant to the evolution of high content imaging: • Evolution and impact of high content imaging: An academic perspective • Evolution and impact of high content imaging: An industry perspective • Evolution of high content image analysis • Evolution of high content data analysis pipelines towards multiparametric and phenotypic profiling applications • The role of data integration and multiomics • The role and evolution of image data repositories and sharing standards • Future perspective of high content imaging hardware and software., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Competing Interests: Declaration of Conflicting Interests J.W-S is a consultant for InSphero AG. The authors declared no other potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Drug resistance is a challenge in anticancer therapy. In many cases, cancers can be resistant to the drug prior to exposure, that is, possess intrinsic drug resistance. However, we lack target-independent methods to anticipate resistance in cancer cell lines or characterize intrinsic drug resistance without a priori knowledge of its cause. We hypothesized that cell morphology could provide an unbiased readout of drug resistance. To test this hypothesis, we used HCT116 cells, a mismatch repair-deficient cancer cell line, to isolate clones that were resistant or sensitive to bortezomib, a well-characterized proteasome inhibitor and anticancer drug to which many cancer cells possess intrinsic resistance. We then expanded these clones and measured high-dimensional single-cell morphology profiles using Cell Painting, a high-content microscopy assay. Our imaging- and computation-based profiling pipeline identified morphological features that differed between resistant and sensitive cells. We used these features to generate a morphological signature of bortezomib resistance. We then employed this morphological signature to analyze a set of HCT116 clones (five resistant and five sensitive) that had not been included in the signature training dataset, and correctly predicted sensitivity to bortezomib in seven cases, in the absence of drug treatment. This signature predicted bortezomib resistance better than resistance to other drugs targeting the ubiquitin-proteasome system, indicating specificity for mechanisms of resistance to bortezomib. Our results establish a proof-of-concept framework for the unbiased analysis of drug resistance using high-content microscopy of cancer cells, in the absence of drug treatment., Competing Interests: MK, AB, DS, BC, YH, SS, AC, TK, GW No competing interests declared, (© 2023, Kelley et al.)

A primary obstacle in translating genetic associations with disease into therapeutic strategies is elucidating the cellular programs affected by genetic risk variants and effector genes. Here, we introduce LipocyteProfiler, a cardiometabolic-disease-oriented high-content image-based profiling tool that enables evaluation of thousands of morphological and cellular profiles that can be systematically linked to genes and genetic variants relevant to cardiometabolic disease. We show that LipocyteProfiler allows surveillance of diverse cellular programs by generating rich context- and process-specific cellular profiles across hepatocyte and adipocyte cell-state transitions. We use LipocyteProfiler to identify known and novel cellular mechanisms altered by polygenic risk of metabolic disease, including insulin resistance, fat distribution, and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse deep phenotypic profiling in lipocytes and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease., Competing Interests: J.C.F. has received consulting honoraria from Goldfinch Bio and Astra Zeneca and speaking honoraria from Novo Nordisk, Astra Zeneca, and Merck for research presentations over which he had full control of content. M.C. holds equity in Waypoint Bio, serves as a consultant for Pfizer, and is a member of the Nestle Scientific Advisory Board. The authors have filed a provisional patent application (63/218,656)., (© 2023.)

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community., Competing Interests: C.S.G.’s spouse was an employee of Alex’s Lemonade Stand Foundation, which was a sponsor of this research. J.A.S., C.L.S., C.J.B., S.J.S., and J.N.T. are or were employees of Alex’s Lemonade Stand Foundation, a sponsor of this research. A.J.W. is a member of the Scientific Advisory boards for Alexion and DayOne Biopharmaceuticals., (© 2023 The Author(s).)

Morphological and gene expression profiling can cost-effectively capture thousands of features in thousands of samples across perturbations by disease, mutation, or drug treatments, but it is unclear to what extent the two modalities capture overlapping versus complementary information. Here, using both the L1000 and Cell Painting assays to profile gene expression and cell morphology, respectively, we perturb human A549 lung cancer cells with 1,327 small molecules from the Drug Repurposing Hub across six doses, providing a data resource including dose-response data from both assays. The two assays capture both shared and complementary information for mapping cell state. Cell Painting profiles from compound perturbations are more reproducible and show more diversity but measure fewer distinct groups of features. Applying unsupervised and supervised methods to predict compound mechanisms of action (MOAs) and gene targets, we find that the two assays not only provide a partially shared but also a complementary view of drug mechanisms. Given the numerous applications of profiling in biology, our analyses provide guidance for planning experiments that profile cells for detecting distinct cell types, disease phenotypes, and response to chemical or genetic perturbations., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

A variational autoencoder (VAE) is a machine learning algorithm, useful for generating a compressed and interpretable latent space. These representations have been generated from various biomedical data types and can be used to produce realistic-looking simulated data. However, standard vanilla VAEs suffer from entangled and uninformative latent spaces, which can be mitigated using other types of VAEs such as β-VAE and MMD-VAE. In this project, we evaluated the ability of VAEs to learn cell morphology characteristics derived from cell images. We trained and evaluated these three VAE variants-Vanilla VAE, β-VAE, and MMD-VAE-on cell morphology readouts and explored the generative capacity of each model to predict compound polypharmacology (the interactions of a drug with more than one target) using an approach called latent space arithmetic (LSA). To test the generalizability of the strategy, we also trained these VAEs using gene expression data of the same compound perturbations and found that gene expression provides complementary information. We found that the β-VAE and MMD-VAE disentangle morphology signals and reveal a more interpretable latent space. We reliably simulated morphology and gene expression readouts from certain compounds thereby predicting cell states perturbed with compounds of known polypharmacology. Inferring cell state for specific drug mechanisms could aid researchers in developing and identifying targeted therapeutics and categorizing off-target effects in the future., Competing Interests: The authors have declared that no competing interests exist.

Software has provided cell biologists the power to quantify specific cellular features in cell images at scale. Before long, these biologists also recognized the potential to extract much more biological information from the same images. From here, the field of image-based profiling, the process of extracting unbiased representations that capture morphological cell state, was born. We are still in the early days of image-based profiling, and it is clear that the many opportunities to interrogate biological systems come with significant challenges. These challenges include building expressive and biologically-relevant representations, adjusting for technical noise, writing generalizable software infrastructure, continuing to foster a culture of open science, and promoting FAIR (findable, accessible, interoperable, and reusable) data. We present a workshop at the Pacific Symposium on Biocomputing 2022 to introduce the field of image-based profiling to the broader computational biology community. In the following document, we introduce image-based profiling, discuss current state-of-the-art methods and limitations, and provide rationale for why now is the perfect time for the field to expand. We also introduce our invited speakers and agenda, which together provide an introduction to the field complemented by in-depth application areas in industry and academia. We also include five lightning talks to complement the invited speakers on various methodological and discovery advances.

Evolving in sync with the computation revolution over the past 30 years, computational biology has emerged as a mature scientific field. While the field has made major contributions toward improving scientific knowledge and human health, individual computational biology practitioners at various institutions often languish in career development. As optimistic biologists passionate about the future of our field, we propose solutions for both eager and reluctant individual scientists, institutions, publishers, funding agencies, and educators to fully embrace computational biology. We believe that in order to pave the way for the next generation of discoveries, we need to improve recognition for computational biologists and better align pathways of career success with pathways of scientific progress. With 10 outlined steps, we call on all adjacent fields to move away from the traditional individual, single-discipline investigator research model and embrace multidisciplinary, data-driven, team science., Competing Interests: The authors have declared that no competing interests exist.

Genetic and chemical perturbations impact diverse cellular phenotypes, including multiple indicators of cell health. These readouts reveal toxicity and antitumorigenic effects relevant to drug discovery and personalized medicine. We developed two customized microscopy assays, one using four targeted reagents and the other three targeted reagents, to collectively measure 70 specific cell health phenotypes including proliferation, apoptosis, reactive oxygen species, DNA damage, and cell cycle stage. We then tested an approach to predict multiple cell health phenotypes using Cell Painting, an inexpensive and scalable image-based morphology assay. In matched CRISPR perturbations of three cancer cell lines, we collected both Cell Painting and cell health data. We found that simple machine learning algorithms can predict many cell health readouts directly from Cell Painting images, at less than half the cost. We hypothesized that these models can be applied to accurately predict cell health assay outcomes for any future or existing Cell Painting dataset. For Cell Painting images from a set of 1500+ compound perturbations across multiple doses, we validated predictions by orthogonal assay readouts. We provide a web app to browse predictions: http://broad.io/cell-health-app. Our approach can be used to add cell health annotations to Cell Painting datasets.

Background: Bone accrual impacts lifelong skeletal health, but genetic discovery has been primarily limited to cross-sectional study designs and hampered by uncertainty about target effector genes. Here, we capture this dynamic phenotype by modeling longitudinal bone accrual across 11,000 bone scans in a cohort of healthy children and adolescents, followed by genome-wide association studies (GWAS) and variant-to-gene mapping with functional follow-up., Results: We identify 40 loci, 35 not previously reported, with various degrees of supportive evidence, half residing in topological associated domains harboring known bone genes. Of several loci potentially associated with later-life fracture risk, a candidate SNP lookup provides the most compelling evidence for rs11195210 (SMC3). Variant-to-gene mapping combining ATAC-seq to assay open chromatin with high-resolution promoter-focused Capture C identifies contacts between GWAS loci and nearby gene promoters. siRNA knockdown of gene expression supports the putative effector gene at three specific loci in two osteoblast cell models. Finally, using CRISPR-Cas9 genome editing, we confirm that the immediate genomic region harboring the putative causal SNP influences PRPF38A expression, a location which is predicted to coincide with a set of binding sites for relevant transcription factors., Conclusions: Using a new longitudinal approach, we expand the number of genetic loci putatively associated with pediatric bone gain. Functional follow-up in appropriate cell models finds novel candidate genes impacting bone accrual. Our data also raise the possibility that the cell fate decision between osteogenic and adipogenic lineages is important in normal bone accrual.

Background: Unsupervised compression algorithms applied to gene expression data extract latent or hidden signals representing technical and biological sources of variation. However, these algorithms require a user to select a biologically appropriate latent space dimensionality. In practice, most researchers fit a single algorithm and latent dimensionality. We sought to determine the extent by which selecting only one fit limits the biological features captured in the latent representations and, consequently, limits what can be discovered with subsequent analyses., Results: We compress gene expression data from three large datasets consisting of adult normal tissue, adult cancer tissue, and pediatric cancer tissue. We train many different models across a large range of latent space dimensionalities and observe various performance differences. We identify more curated pathway gene sets significantly associated with individual dimensions in denoising autoencoder and variational autoencoder models trained using an intermediate number of latent dimensionalities. Combining compressed features across algorithms and dimensionalities captures the most pathway-associated representations. When trained with different latent dimensionalities, models learn strongly associated and generalizable biological representations including sex, neuroblastoma MYCN amplification, and cell types. Stronger signals, such as tumor type, are best captured in models trained at lower dimensionalities, while more subtle signals such as pathway activity are best identified in models trained with more latent dimensionalities., Conclusions: There is no single best latent dimensionality or compression algorithm for analyzing gene expression data. Instead, using features derived from different compression models across multiple latent space dimensionalities enhances biological representations.

Understanding the aberrant transcriptional landscape of neuroblastoma is necessary to provide insight to the underlying influences of the initiation, progression and persistence of this developmental cancer. Here, we present chromatin immunoprecipitation sequencing (ChIP-Seq) data for the oncogenic transcription factors, MYCN and MYC, as well as regulatory histone marks H3K4me1, H3K4me3, H3K27Ac, and H3K27me3 in ten commonly used human neuroblastoma-derived cell line models. In addition, for all of the profiled cell lines we provide ATAC-Seq as a measure of open chromatin. We validate specificity of global MYCN occupancy in MYCN amplified cell lines and functional redundancy of MYC occupancy in MYCN non-amplified cell lines. Finally, we show with H3K27Ac ChIP-Seq that these cell lines retain expression of key neuroblastoma super-enhancers (SE). We anticipate this dataset, coupled with available transcriptomic profiling on the same cell lines, will enable the discovery of novel gene regulatory mechanisms in neuroblastoma.

Recent work has highlighted the tumor microenvironment as a central player in cancer. In particular, interactions between tumor and immune cells may help drive the development of brain tumors such as glioblastoma multiforme (GBM). Despite significant research into the molecular classification of glioblastoma, few studies have characterized in a comprehensive manner the immune infiltrate in situ and within different GBM subtypes.In this study, we use an unbiased, automated immunohistochemistry-based approach to determine the immune phenotype of the four GBM subtypes (classical, mesenchymal, neural and proneural) in a cohort of 98 patients. Tissue Micro Arrays (TMA) were stained for CD20 (B lymphocytes), CD5, CD3, CD4, CD8 (T lymphocytes), CD68 (microglia), and CD163 (bone marrow derived macrophages) antibodies. Using automated image analysis, the percentage of each immune population was calculated with respect to the total tumor cells. Mesenchymal GBMs displayed the highest percentage of microglia, macrophage, and lymphocyte infiltration. CD68 + and CD163 + cells were the most abundant cell populations in all four GBM subtypes, and a higher percentage of CD163 + cells was associated with a worse prognosis. We also compared our results to the relative composition of immune cell type infiltration (using RNA-seq data) across TCGA GBM tumors and validated our results obtained with immunohistochemistry with an external cohort and a different method. The results of this study offer a comprehensive analysis of the distribution and the infiltration of the immune components across the four commonly described GBM subgroups, setting the basis for a more detailed patient classification and new insights that may be used to better apply or design immunotherapies for GBM.

A major driver of multiple myeloma (MM) is thought to be aberrant signaling, yet no kinase inhibitors have proven successful in the clinic. Here, we employed an integrated, systems approach combining phosphoproteomic and transcriptome analysis to dissect cellular signaling in MM to inform precision medicine strategies. Unbiased phosphoproteomics initially revealed differential activation of kinases across MM cell lines and that sensitivity to mammalian target of rapamycin (mTOR) inhibition may be particularly dependent on mTOR kinase baseline activity. We further noted differential activity of immediate downstream effectors of Ras as a function of cell line genotype. We extended these observations to patient transcriptome data in the Multiple Myeloma Research Foundation CoMMpass study. A machine-learning-based classifier identified surprisingly divergent transcriptional outputs between NRAS- and KRAS-mutated tumors. Genetic dependency and gene expression analysis revealed mutated Ras as a selective vulnerability, but not other MAPK pathway genes. Transcriptional analysis further suggested that aberrant MAPK pathway activation is only present in a fraction of RAS-mutated vs wild-type RAS patients. These high-MAPK patients, enriched for NRAS Q61 mutations, have inferior outcomes, whereas RAS mutations overall carry no survival impact. We further developed an interactive software tool to relate pharmacologic and genetic kinase dependencies in myeloma. Collectively, these predictive models identify vulnerable signaling signatures and highlight surprising differences in functional signaling patterns between NRAS and KRAS mutants invisible to the genomic landscape. These results will lead to improved stratification of MM patients in precision medicine trials while also revealing unexplored modes of Ras biology in MM., (© 2019 by The American Society of Hematology.)

Accelerating cures for children with cancer remains an immediate challenge as a result of extensive oncogenic heterogeneity between and within histologies, distinct molecular mechanisms evolving between diagnosis and relapsed disease, and limited therapeutic options. To systematically prioritize and rationally test novel agents in preclinical murine models, researchers within the Pediatric Preclinical Testing Consortium are continuously developing patient-derived xenografts (PDXs)-many of which are refractory to current standard-of-care treatments-from high-risk childhood cancers. Here, we genomically characterize 261 PDX models from 37 unique pediatric cancers; demonstrate faithful recapitulation of histologies and subtypes; and refine our understanding of relapsed disease. In addition, we use expression signatures to classify tumors for TP53 and NF1 pathway inactivation. We anticipate that these data will serve as a resource for pediatric oncology drug development and will guide rational clinical trial design for children with cancer., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Investigators often interpret genome-wide data by analyzing the expression levels of genes within pathways. While this within-pathway analysis is routine, the products of any one pathway can affect the activity of other pathways. Past efforts to identify relationships between biological processes have evaluated overlap in knowledge bases or evaluated changes that occur after specific treatments. Individual experiments can highlight condition-specific pathway-pathway relationships; however, constructing a complete network of such relationships across many conditions requires analyzing results from many studies., Results: We developed PathCORE-T framework by implementing existing methods to identify pathway-pathway transcriptional relationships evident across a broad data compendium. PathCORE-T is applied to the output of feature construction algorithms; it identifies pairs of pathways observed in features more than expected by chance as functionally co-occurring . We demonstrate PathCORE-T by analyzing an existing eADAGE model of a microbial compendium and building and analyzing NMF features from the TCGA dataset of 33 cancer types. The PathCORE-T framework includes a demonstration web interface, with source code, that users can launch to (1) visualize the network and (2) review the expression levels of associated genes in the original data. PathCORE-T creates and displays the network of globally co-occurring pathways based on features observed in a machine learning analysis of gene expression data., Conclusions: The PathCORE-T framework identifies transcriptionally co-occurring pathways from the results of unsupervised analysis of gene expression data and visualizes the relationships between pathways as a network. PathCORE-T recapitulated previously described pathway-pathway relationships and suggested experimentally testable additional hypotheses that remain to be explored., Competing Interests: Not applicableNot applicableThe authors declare that they have no competing interests.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFβ signaling, p53 and β-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy., (Copyright © 2018. Published by Elsevier Inc.)

Precision oncology uses genomic evidence to match patients with treatment but often fails to identify all patients who may respond. The transcriptome of these "hidden responders" may reveal responsive molecular states. We describe and evaluate a machine-learning approach to classify aberrant pathway activity in tumors, which may aid in hidden responder identification. The algorithm integrates RNA-seq, copy number, and mutations from 33 different cancer types across The Cancer Genome Atlas (TCGA) PanCanAtlas project to predict aberrant molecular states in tumors. Applied to the Ras pathway, the method detects Ras activation across cancer types and identifies phenocopying variants. The model, trained on human tumors, can predict response to MEK inhibitors in wild-type Ras cell lines. We also present data that suggest that multiple hits in the Ras pathway confer increased Ras activity. The transcriptome is underused in precision oncology and, combined with machine learning, can aid in the identification of hidden responders., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

DNA damage repair (DDR) pathways modulate cancer risk, progression, and therapeutic response. We systematically analyzed somatic alterations to provide a comprehensive view of DDR deficiency across 33 cancer types. Mutations with accompanying loss of heterozygosity were observed in over 1/3 of DDR genes, including TP53 and BRCA1/2. Other prevalent alterations included epigenetic silencing of the direct repair genes EXO5, MGMT, and ALKBH3 in ∼20% of samples. Homologous recombination deficiency (HRD) was present at varying frequency in many cancer types, most notably ovarian cancer. However, in contrast to ovarian cancer, HRD was associated with worse outcomes in several other cancers. Protein structure-based analyses allowed us to predict functional consequences of rare, recurrent DDR mutations. A new machine-learning-based classifier developed from gene expression data allowed us to identify alterations that phenocopy deleterious TP53 mutations. These frequent DDR gene alterations in many human cancers have functional consequences that may determine cancer progression and guide therapy., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Deep learning describes a class of machine learning algorithms that are capable of combining raw inputs into layers of intermediate features. These algorithms have recently shown impressive results across a variety of domains. Biology and medicine are data-rich disciplines, but the data are complex and often ill-understood. Hence, deep learning techniques may be particularly well suited to solve problems of these fields. We examine applications of deep learning to a variety of biomedical problems-patient classification, fundamental biological processes and treatment of patients-and discuss whether deep learning will be able to transform these tasks or if the biomedical sphere poses unique challenges. Following from an extensive literature review, we find that deep learning has yet to revolutionize biomedicine or definitively resolve any of the most pressing challenges in the field, but promising advances have been made on the prior state of the art. Even though improvements over previous baselines have been modest in general, the recent progress indicates that deep learning methods will provide valuable means for speeding up or aiding human investigation. Though progress has been made linking a specific neural network's prediction to input features, understanding how users should interpret these models to make testable hypotheses about the system under study remains an open challenge. Furthermore, the limited amount of labelled data for training presents problems in some domains, as do legal and privacy constraints on work with sensitive health records. Nonetheless, we foresee deep learning enabling changes at both bench and bedside with the potential to transform several areas of biology and medicine., (© 2018 The Authors.)

Differential expression experiments or other analyses often end in a list of genes. Pathway enrichment analysis is one method to discern important biological signals and patterns from noisy expression data. However, pathway enrichment analysis may perform suboptimally in situations where there are multiple implicated pathways - such as in the case of genes that define subtypes of complex diseases. Our simulation study shows that in this setting, standard overrepresentation analysis identifies many false positive pathways along with the true positives. These false positives hamper investigators' attempts to glean biological insights from enrichment analysis. We develop and evaluate an approach that combines community detection over functional networks with pathway enrichment to reduce false positives. Our simulation study demonstrates that a large reduction in false positives can be obtained with a small decrease in power. Though we hypothesized that multiple communities might underlie previously described subtypes of high-grade serous ovarian cancer and applied this approach, our results do not support this hypothesis. In summary, applying community detection before enrichment analysis may ease interpretation for complex gene sets that represent multiple distinct pathways.

The Cancer Genome Atlas (TCGA) has profiled over 10,000 tumors across 33 different cancer-types for many genomic features, including gene expression levels. Gene expression measurements capture substantial information about the state of each tumor. Certain classes of deep neural network models are capable of learning a meaningful latent space. Such a latent space could be used to explore and generate hypothetical gene expression profiles under various types of molecular and genetic perturbation. For example, one might wish to use such a model to predict a tumor's response to specific therapies or to characterize complex gene expression activations existing in differential proportions in different tumors. Variational autoencoders (VAEs) are a deep neural network approach capable of generating meaningful latent spaces for image and text data. In this work, we sought to determine the extent to which a VAE can be trained to model cancer gene expression, and whether or not such a VAE would capture biologically-relevant features. In the following report, we introduce a VAE trained on TCGA pan-cancer RNA-seq data, identify specific patterns in the VAE encoded features, and discuss potential merits of the approach. We name our method "Tybalt" after an instigative, cat-like character who sets a cascading chain of events in motion in Shakespeare's "Romeo and Juliet". From a systems biology perspective, Tybalt could one day aid in cancer stratification or predict specific activated expression patterns that would result from genetic changes or treatment effects.

Genome-wide association studies (GWAS) have contributed significantly to the understanding of complex disease genetics. However, GWAS only report association signals and do not necessarily identify culprit genes. As most signals occur in non-coding regions of the genome, it is often challenging to assign genomic variants to the underlying causal mechanism(s). Topologically associating domains (TADs) are primarily cell-type-independent genomic regions that define interactome boundaries and can aid in the designation of limits within which an association most likely impacts gene function. We describe and validate a computational method that uses the genic content of TADs to prioritize candidate genes. Our method, called 'TAD&#95;Pathways', performs a Gene Ontology (GO) analysis over genes that reside within TAD boundaries corresponding to GWAS signals for a given trait or disease. Applying our pipeline to the bone mineral density (BMD) GWAS catalog, we identify 'Skeletal System Development' (Benjamini-Hochberg adjusted P=1.02x10 -5 ) as the top-ranked pathway. In many cases, our method implicated a gene other than the nearest gene. Our molecular experiments describe a novel example: ACP2, implicated near the canonical 'ARHGAP1' locus. We found ACP2 to be an important regulator of osteoblast metabolism, whereas ARHGAP1 was not supported. Our results via BMD, for example, demonstrate how basic principles of three-dimensional genome organization can define biologically informed association windows.

Breast cancer is a complex disease consisting of four distinct molecular subtypes. DNA methylation-based (DNAm) studies in tumors are complicated further by disease heterogeneity. In the present study, we compared DNAm in breast tumors with normal-adjacent breast samples from The Cancer Genome Atlas (TCGA). We constructed models stratified by tumor stage and PAM50 molecular subtype and performed cell-type reference-free deconvolution to control for cellular heterogeneity. We identified nineteen differentially methylated gene regions (DMGRs) in early stage tumors across eleven genes (AGRN, C1orf170, FAM41C, FLJ39609, HES4, ISG15, KLHL17, NOC2L, PLEKHN1, SAMD11, WASH5P). These regions were consistently differentially methylated in every subtype and all implicated genes are localized to the chromosomal cytoband 1p36.3. Seventeen of these DMGRs were independently validated in a similar analysis of an external data set. The identification and validation of shared DNAm alterations across tumor subtypes in early stage tumors advances our understanding of common biology underlying breast carcinogenesis and may contribute to biomarker development. We also discuss evidence of the specific importance and potential function of 1p36 in cancer.

Purpose of Review: Only recently has it become clear that epithelial ovarian cancer (EOC) is comprised of such distinct histotypes--with different cells of origin, morphology, molecular features, epidemiologic factors, clinical features, and survival patterns-that they can be thought of as different diseases sharing an anatomical location. Herein, we review opportunities and challenges in studying EOC heterogeneity., Recent Findings: The 2014 World Health Organization diagnostic guidelines incorporate accumulated evidence that high- and low-grade serous tumors have different underlying pathogenesis, and that, on the basis of shared molecular features, most high grade tumors, including some previously classified as endometrioid, are now considered to be high-grade serous. At the same time, several studies have reported that high-grade serous EOC, which is the most common histotype, is itself made up of reproducible subtypes discernable by gene expression patterns., Summary: These major advances in understanding set the stage for a new era of research on EOC risk and clinical outcomes with the potential to reduce morbidity and mortality. We highlight the need for multidisciplinary studies with pathology review using the current guidelines, further molecular characterization of the histotypes and subtypes, inclusion of women of diverse racial/ethnic and socioeconomic backgrounds, and updated epidemiologic and clinical data relevant to current generations of women at risk of EOC., Competing Interests: Conflict of Interest Jennifer Anne Doherty, Lauren Cole Peres, Chen Wang, and Gregory P. Way each declare no potential conflicts of interest. Casey S. Greene reports grants from Gordon and Betty Moore Foundation during the conduct of the study and personal fees from SomaLogic outside the submitted work. Joellen M. Schildkraut reports grants from National Cancer Institute during the conduct of the study.

Background: We have identified molecules that exhibit synthetic lethality in cells with loss of the neurofibromin 1 (NF1) tumor suppressor gene. However, recognizing tumors that have inactivation of the NF1 tumor suppressor function is challenging because the loss may occur via mechanisms that do not involve mutation of the genomic locus. Degradation of the NF1 protein, independent of NF1 mutation status, phenocopies inactivating mutations to drive tumors in human glioma cell lines. NF1 inactivation may alter the transcriptional landscape of a tumor and allow a machine learning classifier to detect which tumors will benefit from synthetic lethal molecules., Results: We developed a strategy to predict tumors with low NF1 activity and hence tumors that may respond to treatments that target cells lacking NF1. Using RNAseq data from The Cancer Genome Atlas (TCGA), we trained an ensemble of 500 logistic regression classifiers that integrates mutation status with whole transcriptomes to predict NF1 inactivation in glioblastoma (GBM). On TCGA data, the classifier detected NF1 mutated tumors (test set area under the receiver operating characteristic curve (AUROC) mean = 0.77, 95% quantile = 0.53 - 0.95) over 50 random initializations. On RNA-Seq data transformed into the space of gene expression microarrays, this method produced a classifier with similar performance (test set AUROC mean = 0.77, 95% quantile = 0.53 - 0.96). We applied our ensemble classifier trained on the transformed TCGA data to a microarray validation set of 12 samples with matched RNA and NF1 protein-level measurements. The classifier's NF1 score was associated with NF1 protein concentration in these samples., Conclusions: We demonstrate that TCGA can be used to train accurate predictors of NF1 inactivation in GBM. The ensemble classifier performed well for samples with very high or very low NF1 protein concentrations but had mixed performance in samples with intermediate NF1 concentrations. Nevertheless, high-performing and validated predictors have the potential to be paired with targeted therapies and personalized medicine.

Four gene expression subtypes of high-grade serous ovarian cancer (HGSC) have been previously described. In these early studies, a fraction of samples that did not fit well into the four subtype classifications were excluded. Therefore, we sought to systematically determine the concordance of transcriptomic HGSC subtypes across populations without removing any samples. We created a bioinformatics pipeline to independently cluster the five largest mRNA expression datasets using k-means and nonnegative matrix factorization (NMF). We summarized differential expression patterns to compare clusters across studies. While previous studies reported four subtypes, our cross-population comparison does not support four. Because these results contrast with previous reports, we attempted to reproduce analyses performed in those studies. Our results suggest that early results favoring four subtypes may have been driven by the inclusion of serous borderline tumors. In summary, our analysis suggests that either two or three, but not four, gene expression subtypes are most consistent across datasets., (Copyright © 2016 Way et al.)

A behavioral syndrome exists when specific behaviors interact under different contexts. Zebrafish have been test subjects in recent studies and it is important to standardize protocols to ensure proper analyses and interpretations. In our previous studies, we have measured boldness by monitoring a series of behaviors (time near surface, latency in transitions, number of transitions, and darts) in a 1.5 L trapezoidal tank. Likewise, we quantified aggression by observing bites, lateral displays, darts, and time near an inclined mirror in a rectangular 19 L tank. By dividing a 76 L tank into thirds, we also examined shoaling preferences. The shoaling assay is a highly customizable assay and can be tailored for specific hypotheses. However, protocols for this assay also must be standardized, yet flexible enough for customization. In previous studies, end chambers were either empty, contained 5 or 10 zebrafish, or 5 pearl danios (D. albolineatus). In the following manuscript, we present a detailed protocol and representative data that accompany successful applications of the protocol, which will allow for replication of behavioral syndrome experiments.

A behavioral syndrome is considered present when individuals consistently express correlated behaviors across two or more axes of behavior. These axes of behavior are shy-bold, exploration-avoidance, activity, aggression, and sociability. In this study we examined aggression, boldness and sociability (shoaling) within a juvenile convict cichlid brood (Amatitlania nigrofasciatus). Because young convict cichlids are social, we used methodologies commonly used by ethologists studying social fishes. We did not detect an aggression-boldness behavioral syndrome, but we did find that the aggression, boldness, and possibly the exploration behavioral axes play significant roles in shaping the observed variation in individual convict cichlid behavior. While juvenile convict cichlids did express a shoaling preference, this social preference was likely convoluted by aggressive interactions, despite the small size and young age of the fish. There is a need for the development of behavioral assays that allow for more reliable measurement of behavioral axes in juvenile neo-tropical cichlids., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Aggression is a psychological construct that is commonly used to classify zebrafish behavior. Aggression is a complex trait that can be difficult to accurately measure. The literature on fish behavior describes many different methodologies to examine aggression, which, we believe, have not been compared in a formal manner. In this study we observed 19 individual zebrafish (Danio rerio) and quantified bites, lateral displays, charges, darts, and time near the stimulus in six common assays used to measure aggression. The methodologies included an inclined mirror assay, two flat mirror assays with different acclimation periods, a live conspecific assay, a clay model stimulus assay, and a video recording assay. Our results indicate high repeatability in most aggressive behaviors over time, which confirms the value of each assay to measure personality. However, our results also indicate significant differences between the assays. Specifically, assays using a flat mirror or live conspecific as a stimulus for aggression elicited more attempted bites than an inclined mirror, a clay model stimulus, or a video recording stimulus. Furthermore, the inclined mirror stimulus provoked more darts than any other assay. The results suggest the need for researchers to consider specific research goals when selecting the appropriate stimulus to provoke aggression in zebrafish.

A behavioral syndrome is observed in a population when specific behaviors overlap at the individual level in different contexts. Here, we explore boldness and aggression personality spectra, the repeatability of shoaling, and possible associated correlations between the behaviors in a population of lab-reared zebrafish (Danio rerio). Our findings describe a sex-specific boldness-shoaling behavioral syndrome, as a link between boldness and shoaling behaviors is detected. The results indicate that bold males are likely to have a stronger shoaling propensity than shy males for unfamiliar conspecifics. Conversely, bold females are more likely to shoal than shy females, but only when presented with heterospecific individuals. Additionally, aggression does not correlate with boldness or shoaling propensity for either sex. A positive relationship between boldness and shoaling that differs by sex is contrary to most of the present literature, but could help to explain population dynamics and may also have evolutionary implications., (Copyright © 2015 Elsevier B.V. All rights reserved.)