Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

Integrated Proteogenomic Characterization across Major Histological Types of Pediatric Brain Cancer.

We report a comprehensive proteogenomics analysis, including whole-genome sequencing, RNA sequencing, and proteomics and phosphoproteomics profiling, of 218 tumors across 7 histological types of childhood brain cancer: low-grade glioma (n = 93), ependymoma (32), high-grade glioma (25), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16), and atypical teratoid rhabdoid tumor (12). Proteomics data identify common biological themes that span histological boundaries, suggesting that treatments used for one histological type may be applied effectively to other tumors sharing similar proteomics features. Immune landscape characterization reveals diverse tumor microenvironments across and within diagnoses. Proteomics data further reveal functional effects of somatic mutations and copy number variations (CNVs) not evident in transcriptomics data. Kinase-substrate association and co-expression network analysis identify important biological mechanisms of tumorigenesis. This is the first large-scale proteogenomics analysis across traditional histological boundaries to uncover foundational pediatric brain tumor biology and inform rational treatment selection.

Kinome-Wide Virtual Screening by Multi-Task Deep Learning.

Deep learning is a machine learning technique to model high-level abstractions in data by utilizing a graph composed of multiple processing layers that experience various linear and non-linear transformations. This technique has been shown to perform well for applications in drug discovery, utilizing structural features of small molecules to predict activity. Here, we report a large-scale study to predict the activity of small molecules across the human kinome-a major family of drug targets, particularly in anti-cancer agents. While small-molecule kinase inhibitors exhibit impressive clinical efficacy in several different diseases, resistance often arises through adaptive kinome reprogramming or subpopulation diversity. Polypharmacology and combination therapies offer potential therapeutic strategies for patients with resistant diseases. Their development would benefit from a more comprehensive and dense knowledge of small-molecule inhibition across the human kinome. Leveraging over 650,000 bioactivity annotations for more than 300,000 small molecules, we evaluated multiple machine learning methods to predict the small-molecule inhibition of 342 kinases across the human kinome. Our results demonstrated that multi-task deep neural networks outperformed classical single-task methods, offering the potential for conducting large-scale virtual screening, predicting activity profiles, and bridging the gaps in the available data.

Machine and deep learning approaches for cancer drug repurposing.

Knowledge of the underpinnings of cancer initiation, progression and metastasis has increased exponentially in recent years. Advanced "omics" coupled with machine learning and artificial intelligence (deep learning) methods have helped elucidate targets and pathways critical to those processes that may be amenable to pharmacologic modulation. However, the current anti-cancer therapeutic armamentarium continues to lag behind. As the cost of developing a new drug remains prohibitively expensive, repurposing of existing approved and investigational drugs is sought after given known safety profiles and reduction in the cost barrier. Notably, successes in oncologic drug repurposing have been infrequent. Computational in-silico strategies have been developed to aid in modeling biological processes to find new disease-relevant targets and discovering novel drug-target and drug-phenotype associations. Machine and deep learning methods have especially enabled leaps in those successes. This review will discuss these methods as they pertain to cancer biology as well as immunomodulation for drug repurposing opportunities in oncologic diseases.

LINCS Data Portal 2.0: next generation access point for perturbation-response signatures.

The Library of Integrated Network-Based Cellular Signatures (LINCS) is an NIH Common Fund program with the goal of generating a large-scale and comprehensive catalogue of perturbation-response signatures by utilizing a diverse collection of perturbations across many model systems and assay types. The LINCS Data Portal (LDP) has been the primary access point for the compendium of LINCS data and has been widely utilized. Here, we report the first major update of LDP (http://lincsportal.ccs.miami.edu/signatures) with substantial changes in the data architecture and APIs, a completely redesigned user interface, and enhanced curated metadata annotations to support more advanced, intuitive and deeper querying, exploration and analysis capabilities. The cornerstone of this update has been the decision to reprocess all high-level LINCS datasets and make them accessible at the data point level enabling users to directly access and download any subset of signatures across the entire library independent from the originating source, project or assay. Access to the individual signatures also enables the newly implemented signature search functionality, which utilizes the iLINCS platform to identify conditions that mimic or reverse gene set queries. A newly designed query interface enables global metadata search with autosuggest across all annotations associated with perturbations, model systems, and signatures.

DrugCentral 2018: an update.

DrugCentral is a drug information resource (http://drugcentral.org) open to the public since 2016 and previously described in the 2017 Nucleic Acids Research Database issue. Since the 2016 release, 103 new approved drugs were updated. The following new data sources have been included: Food and Drug Administration (FDA) Adverse Event Reporting System (FAERS), FDA Orange Book information, L1000 gene perturbation profile distance/similarity matrices and estimated protonation constants. New and existing entries have been updated with the latest information from scientific literature, drug labels and external databases. The web interface has been updated to display and query new data. The full database dump and data files are available for download from the DrugCentral website.

Data Portal for the Library of Integrated Network-based Cellular Signatures (LINCS) program: integrated access to diverse large-scale cellular perturbation response data.

The Library of Integrated Network-based Cellular Signatures (LINCS) program is a national consortium funded by the NIH to generate a diverse and extensive reference library of cell-based perturbation-response signatures, along with novel data analytics tools to improve our understanding of human diseases at the systems level. In contrast to other large-scale data generation efforts, LINCS Data and Signature Generation Centers (DSGCs) employ a wide range of assay technologies cataloging diverse cellular responses. Integration of, and unified access to LINCS data has therefore been particularly challenging. The Big Data to Knowledge (BD2K) LINCS Data Coordination and Integration Center (DCIC) has developed data standards specifications, data processing pipelines, and a suite of end-user software tools to integrate and annotate LINCS-generated data, to make LINCS signatures searchable and usable for different types of users. Here, we describe the LINCS Data Portal (LDP) (http://lincsportal.ccs.miami.edu/), a unified web interface to access datasets generated by the LINCS DSGCs, and its underlying database, LINCS Data Registry (LDR). LINCS data served on the LDP contains extensive metadata and curated annotations. We highlight the features of the LDP user interface that is designed to enable search, browsing, exploration, download and analysis of LINCS data and related curated content.

Ontological representation, integration, and analysis of LINCS cell line cells and their cellular responses.

BACKGROUND: Aiming to understand cellular responses to different perturbations, the NIH Common Fund Library of Integrated Network-based Cellular Signatures (LINCS) program involves many institutes and laboratories working on over a thousand cell lines. The community-based Cell Line Ontology (CLO) is selected as the default ontology for LINCS cell line representation and integration. RESULTS: CLO has consistently represented all 1097 LINCS cell lines and included information extracted from the LINCS Data Portal and ChEMBL. Using MCF 10A cell line cells as an example, we demonstrated how to ontologically model LINCS cellular signatures such as their non-tumorigenic epithelial cell type, three-dimensional growth, latrunculin-A-induced actin depolymerization and apoptosis, and cell line transfection. A CLO subset view of LINCS cell lines, named LINCS-CLOview, was generated to support systematic LINCS cell line analysis and queries. In summary, LINCS cell lines are currently associated with 43 cell types, 131 tissues and organs, and 121 cancer types. The LINCS-CLO view information can be queried using SPARQL scripts. CONCLUSIONS: CLO was used to support ontological representation, integration, and analysis of over a thousand LINCS cell line cells and their cellular responses.

Epigenetic pathways and glioblastoma treatment: insights from signaling cascades.

There is an urgent need to identify novel therapies for glioblastoma (GBM) as most therapies are ineffective. A first step in this process is to identify and validate targets for therapeutic intervention. Epigenetic modulators have emerged as attractive drug targets in several cancers including GBM. These epigenetic regulators affect gene expression without changing the DNA sequence. Recent studies suggest that epigenetic regulators interact with drivers of GBM cell and stem-like cell proliferation. These drivers include components of the Notch, Hedgehog, and Wingless (WNT) pathways. We highlight recent studies connecting epigenetic and signaling pathways in GBM. We also review systems and big data approaches for identifying patient specific therapies in GBM. Collectively, these studies will identify drug combinations that may be effective in GBM and other cancers.

Support of the histaminergic hypothesis in Tourette syndrome: association of the histamine decarboxylase gene in a large sample of families.

BACKGROUND: Gilles de la Tourette Syndrome is a neurodevelopmental disorder that is caused by the interaction of environment with a complex genetic background. The genetic etiology of the disorder remains, so far, elusive, although multiple promising leads have been recently reported. The recent implication of the histamine decarboxylase (HDC) gene, the key enzyme in histamine production, raises the intriguing hypothesis of a possible role of histaminergic dysfunction leading to TS onset. METHODS: Following up on the finding of a nonsense mutation in a single family with TS, we investigated variation across the HDC gene for association with TS. As a result of a collaborative international effort, we studied a large sample of 520 nuclear families originating from seven European populations (Greek, Hungarian, Italian, Polish, German, Albanian, Spanish) as well as a sample collected in Canada. RESULTS AND CONCLUSIONS: Interrogating 12 tagging SNPs (tSNP) across the HDC region, we find strong over-transmission of alleles at two SNPs (rs854150 and rs1894236) in the complete sample, as well as a statistically significant associated haplotypes. Analysis of individual populations also reveals signals of association in the Canadian, German and Italian samples. Our results provide strong support for the histaminergic hypothesis in TS etiology and point to a possible role of histamine pathways in neuronal development.

Mass spectrometry database of lacrimal gland adenoid cystic carcinoma and normal lacrimal gland tissue identifies extracellular matrix remodeling in these tumors.

Adenoid cystic carcinoma of the lacrimal gland (LGACC) is a slow-growing but aggressive orbital malignancy. Due to the rarity of LGACC, it is poorly understood, which makes diagnosing, treating, and monitoring disease progression difficult. The aim is to understand the molecular drivers of LGACC further to identify potential targets for treating this cancer. Mass spectrometry was performed on LGACC and normal lacrimal gland samples to examine the differentially expressed proteins to understand this cancer's proteomic characteristics. Downstream gene ontology and pathway analysis revealed the extracellular matrix is the most upregulated process in LGACC. This data serves as a resource for further understanding LGACC and identifying potential treatment targets. This dataset is publicly available.

Spatial Transcriptomics Identifies Expression Signatures Specific to Lacrimal Gland Adenoid Cystic Carcinoma Cells.

Although primary tumors of the lacrimal gland are rare, adenoid cystic carcinoma (ACC) is the most common and lethal epithelial lacrimal gland malignancy. Traditional management of lacrimal gland adenoid cystic carcinoma (LGACC) involves the removal of the eye and surrounding socket contents, followed by chemoradiation. Even with this radical treatment, the 10-year survival rate for LGACC is 20% given the propensity for recurrence and metastasis. Due to the rarity of LGACC, its pathobiology is not well-understood, leading to difficulties in diagnosis, treatment, and effective management. Here, we integrate bulk RNA sequencing (RNA-seq) and spatial transcriptomics to identify a specific LGACC gene signature that can inform novel targeted therapies. Of the 3499 differentially expressed genes identified by bulk RNA-seq, the results of our spatial transcriptomic analysis reveal 15 upregulated and 12 downregulated genes that specifically arise from LGACC cells, whereas fibroblasts, reactive fibrotic tissue, and nervous and skeletal muscle account for the remaining bulk RNA-seq signature. In light of the analysis, we identified a transitional state cell or stem cell cluster. The results of the pathway analysis identified the upregulation of PI3K-Akt signaling, IL-17 signaling, and multiple other cancer pathways. This study provides insights into the molecular and cellular landscape of LGACC, which can inform new, targeted therapies to improve patient outcomes.

Correction: The Genetic Landscape of Ocular Adnexa MALT Lymphoma Reveals Frequent Aberrations in NFAT and MEF2B Signaling Pathways.

[This corrects the article DOI: 10.1158/2767-9764.CRC-21-0022.][This corrects the article DOI: 10.1158/2767-9764.CRC-21-0022.].

Exploring genomic structure differences and similarities between the Greek and European HapMap populations: implications for association studies.

Studies of the genomic structure of the Greek population and Southeastern Europe are limited, despite the central position of the area as a gateway for human migrations into Europe. HapMap has provided a unique tool for the analysis of human genetic variation. Europe is represented by the CEU (Northwestern Europe) and the TSI populations (Tuscan Italians from Southern Europe), which serve as reference for the design of genetic association studies. Furthermore, genetic association findings are often transferred to unstudied populations. Although initial studies support the fact that the CEU can, in general, be used as reference for the selection of tagging SNPs in European populations, this has not been extensively studied across Europe. We set out to explore the genomic structure of the Greek population (56 individuals) and compare it to the HapMap TSI and CEU populations. We studied 1112 SNPs (27 regions, 13 chromosomes). Although the HapMap European populations are, in general, a good reference for the Greek population, regions of population differentiation do exist and results should not be light-heartedly generalized. We conclude that, perhaps due to the individual evolutionary history of each genomic region, geographic proximity is not always a perfect guide for selecting a reference population for an unstudied population.

Connecting omics signatures and revealing biological mechanisms with iLINCS.

There are only a few platforms that integrate multiple omics data types, bioinformatics tools, and interfaces for integrative analyses and visualization that do not require programming skills. Here we present iLINCS ( http://ilincs.org ), an integrative web-based platform for analysis of omics data and signatures of cellular perturbations. The platform facilitates mining and re-analysis of the large collection of omics datasets (>34,000), pre-computed signatures (>200,000), and their connections, as well as the analysis of user-submitted omics signatures of diseases and cellular perturbations. iLINCS analysis workflows integrate vast omics data resources and a range of analytics and interactive visualization tools into a comprehensive platform for analysis of omics signatures. iLINCS user-friendly interfaces enable execution of sophisticated analyses of omics signatures, mechanism of action analysis, and signature-driven drug repositioning. We illustrate the utility of iLINCS with three use cases involving analysis of cancer proteogenomic signatures, COVID 19 transcriptomic signatures and mTOR signaling.

A multiparametric pharmacogenomic strategy for drug repositioning predicts therapeutic efficacy for glioblastoma cell lines.

BACKGROUND: Poor prognosis of glioblastoma patients and the extensive heterogeneity of glioblastoma at both the molecular and cellular level necessitates developing novel individualized treatment modalities via genomics-driven approaches. METHODS: This study leverages numerous pharmacogenomic and tissue databases to examine drug repositioning for glioblastoma. RNA-seq of glioblastoma tumor samples from The Cancer Genome Atlas (TCGA, n = 117) were compared to "normal" frontal lobe samples from Genotype-Tissue Expression Portal (GTEX, n = 120) to find differentially expressed genes (DEGs). Using compound gene expression data and drug activity data from the Library of Integrated Network-Based Cellular Signatures (LINCS, n = 66,512 compounds) CCLE (71 glioma cell lines), and Chemical European Molecular Biology Laboratory (ChEMBL) platforms, we employed a summarized reversal gene expression metric (sRGES) to "reverse" the resultant disease signature for GBM and its subtypes. A multiparametric strategy was employed to stratify compounds capable of blood-brain barrier penetrance with a favorable pharmacokinetic profile (CNS-MPO). RESULTS: Significant correlations were identified between sRGES and drug efficacy in GBM cell lines in both ChEMBL(r = 0.37, P < .001) and Cancer Therapeutic Response Portal (CTRP) databases (r = 0.35, P < 0.001). Our multiparametric algorithm identified two classes of drugs with highest sRGES and CNS-MPO: HDAC inhibitors (vorinostat and entinostat) and topoisomerase inhibitors suitable for drug repurposing. CONCLUSIONS: Our studies suggest that reversal of glioblastoma disease signature correlates with drug potency for various GBM subtypes. This multiparametric approach may set the foundation for an early-phase personalized -omics clinical trial for glioblastoma by effectively identifying drugs that are capable of reversing the disease signature and have favorable pharmacokinetic and safety profiles.

The Genetic Landscape of Ocular Adnexa MALT Lymphoma Reveals Frequent Aberrations in NFAT and MEF2B Signaling Pathways.

A comprehensive constellation of somatic non-silent mutations and copy number (CN) variations in ocular adnexa marginal zone lymphoma (OAMZL) is unknown. By utilizing whole-exome sequencing in 69 tumors we define the genetic landscape of OAMZL. Mutations and CN changes in CABIN1 (30%), RHOA (26%), TBL1XR1 (22%), and CREBBP (17%) and inactivation of TNFAIP3 (26%) were among the most common aberrations. Candidate cancer driver genes cluster in the B-cell receptor (BCR), NFkB, NOTCH and NFAT signaling pathways. One of the most commonly altered genes is CABIN1, a calcineurin inhibitor acting as a negative regulator of the NFAT and MEF2B transcriptional activity. CABIN1 deletions enhance BCR-stimulated NFAT and MEF2B transcriptional activity, while CABIN1 mutations enhance only MEF2B transcriptional activity by impairing binding of mSin3a to CABIN1. Our data provide an unbiased identification of genetically altered genes that may play a role in the molecular pathogenesis of OAMZL and serve as therapeutic targets.