Corrigendum: Rare cell variability and drug-induced reprogramming as a mode of cancer drug resistance.

This corrects the article DOI: 10.1038/nature22794.

Rare cell variability and drug-induced reprogramming as a mode of cancer drug resistance.

Therapies that target signalling molecules that are mutated in cancers can often have substantial short-term effects, but the emergence of resistant cancer cells is a major barrier to full cures. Resistance can result from secondary mutations, but in other cases there is no clear genetic cause, raising the possibility of non-genetic rare cell variability. Here we show that human melanoma cells can display profound transcriptional variability at the single-cell level that predicts which cells will ultimately resist drug treatment. This variability involves infrequent, semi-coordinated transcription of a number of resistance markers at high levels in a very small percentage of cells. The addition of drug then induces epigenetic reprogramming in these cells, converting the transient transcriptional state to a stably resistant state. This reprogramming begins with a loss of SOX10-mediated differentiation followed by activation of new signalling pathways, partially mediated by the activity of the transcription factors JUN and/or AP-1 and TEAD. Our work reveals the multistage nature of the acquisition of drug resistance and provides a framework for understanding resistance dynamics in single cells. We find that other cell types also exhibit sporadic expression of many of these same marker genes, suggesting the existence of a general program in which expression is displayed in rare subpopulations of cells.

A comprehensive characterisation of phaeochromocytoma and paraganglioma tumours through histone protein profiling, DNA methylation and transcriptomic analysis genome wide.

BACKGROUND: Phaeochromocytomas and paragangliomas (PPGLs) are rare neuroendocrine tumours. Pathogenic variants have been identified in more than 15 susceptibility genes; associated tumours are grouped into three Clusters, reinforced by their transcriptional profiles. Cluster 1A PPGLs have pathogenic variants affecting enzymes of the tricarboxylic acid cycle, including succinate dehydrogenase. Within inherited PPGLs, these are the most common. PPGL tumours are known to undergo epigenetic reprograming, and here, we report on global histone post-translational modifications and DNA methylation levels, alongside clinical phenotypes. RESULTS: Out of the 25 histone post-translational modifications examined, Cluster 1A PPGLs were distinguished from other tumours by a decrease in hyper-acetylated peptides and an increase in H3K4me2. DNA methylation was compared between tumours from individuals who developed metastatic disease versus those that did not. The majority of differentially methylated sites identified tended to be completely methylated or unmethylated in non-metastatic tumours, with low inter-sample variance. Metastatic tumours by contrast consistently had an intermediate DNA methylation state, including the ephrin receptor EPHA4 and its ligand EFNA3. Gene expression analyses performed to identify genes involved in metastatic tumour behaviour pin-pointed a number of genes previously described as mis-regulated in Cluster 1A tumours, as well as highlighting the tumour suppressor RGS22 and the pituitary tumour-transforming gene PTTG1. CONCLUSIONS: Combined transcriptomic and DNA methylation analyses revealed aberrant pathways, including ones that could be implicated in metastatic phenotypes and, for the first time, we report a decrease in hyper-acetylated histone marks in Cluster 1 PPGLs.

Correction: Anastopoulos et al. Multi-Drug Featurization and Deep Learning Improve Patient-Specific Predictions of Adverse Events. Int. J. Environ. Res. Public Health 2021, 18, 2600.

In the original publication [...].

Multi-Drug Featurization and Deep Learning Improve Patient-Specific Predictions of Adverse Events.

While the clinical approval process is able to filter out medications whose utility does not offset their adverse drug reaction profile in humans, it is not well suited to characterizing lower frequency issues and idiosyncratic multi-drug interactions that can happen in real world diverse patient populations. With a growing abundance of real-world evidence databases containing hundreds of thousands of patient records, it is now feasible to build machine learning models that incorporate individual patient information to provide personalized adverse event predictions. In this study, we build models that integrate patient specific demographic, clinical, and genetic features (when available) with drug structure to predict adverse drug reactions. We develop an extensible graph convolutional approach to be able to integrate molecular effects from the variable number of medications a typical patient may be taking. Our model outperforms standard machine learning methods at the tasks of predicting hospitalization and death in the UK Biobank dataset yielding an R2 of 0.37 and an AUC of 0.90, respectively. We believe our model has potential for evaluating new therapeutic compounds for individualized toxicities in real world diverse populations. It can also be used to prioritize medications when there are multiple options being considered for treatment.

Using Transcriptional Signatures to Find Cancer Drivers with LURE.

Cancer genome projects have produced multidimensional datasets on thousands of samples. Yet, depending on the tumor type, 5-50% of samples have no known driving event. We introduce a semi-supervised method called Learning UnRealized Events (LURE) that uses a progressive label learning framework and minimum spanning analysis to predict cancer drivers based on their altered samples sharing a gene expression signature with the samples of a known event. We demonstrate the utility of the method on the TCGA Pan-Cancer Atlas dataset for which it produced a high-confidence result relating 59 new connections to 18 known mutation events including alterations in the same gene, family, and pathway. We give examples of predicted drivers involved in TP53, telomere maintenance, and MAPK/RTK signaling pathways. LURE identifies connections between genes with no known prior relationship, some of which may offer clues for targeting specific forms of cancer. Code and Supplemental Material are available on the LURE website: https://sysbiowiki.soe.ucsc.edu/lure.

CODEX2: full-spectrum copy number variation detection by high-throughput DNA sequencing.

High-throughput DNA sequencing enables detection of copy number variations (CNVs) on the genome-wide scale with finer resolution compared to array-based methods but suffers from biases and artifacts that lead to false discoveries and low sensitivity. We describe CODEX2, as a statistical framework for full-spectrum CNV profiling that is sensitive for variants with both common and rare population frequencies and that is applicable to study designs with and without negative control samples. We demonstrate and evaluate CODEX2 on whole-exome and targeted sequencing data, where biases are the most prominent. CODEX2 outperforms existing methods and, in particular, significantly improves sensitivity for common CNVs.

BRCA locus-specific loss of heterozygosity in germline BRCA1 and BRCA2 carriers.

Complete loss of BRCA1 or BRCA2 function is associated with sensitivity to DNA damaging agents. However, not all BRCA1 and BRCA2 germline mutation-associated tumors respond. Herein we report analyses of 160 BRCA1 and BRCA2 germline mutation-associated breast and ovarian tumors. Retention of the normal BRCA1 or BRCA2 allele (absence of locus-specific loss of heterozygosity (LOH)) is observed in 7% of BRCA1 ovarian, 16% of BRCA2 ovarian, 10% of BRCA1 breast, and 46% of BRCA2 breast tumors. These tumors have equivalent homologous recombination deficiency scores to sporadic tumors, significantly lower than scores in tumors with locus-specific LOH (ovarian, P = 0.0004; breast P < 0.0001, two-tailed Student's t-test). Absence of locus-specific LOH is associated with decreased overall survival in ovarian cancer patients treated with platinum chemotherapy (P = 0.01, log-rank test). Locus-specific LOH may be a clinically useful biomarker to predict primary resistance to DNA damaging agents in patients with germline BRCA1 and BRCA2 mutations.Most tumours associated with germline BRCA1/BRCA2 loss of function mutations respond to DNA damaging agents, however, some do not. Herein, the authors identify that a subset of breast/ovarian tumors retain a normal allele, which is associated with decreased overall survival after DNA damage-inducing platinum chemotherapy.

Genetic and Genomic Characterization of 462 Melanoma Patient-Derived Xenografts, Tumor Biopsies, and Cell Lines.

Tumor-sequencing studies have revealed the widespread genetic diversity of melanoma. Sequencing of 108 genes previously implicated in melanomagenesis was performed on 462 patient-derived xenografts (PDXs), cell lines, and tumors to identify mutational and copy number aberrations. Samples came from 371 unique individuals: 263 were naive to treatment, and 108 were previously treated with targeted therapy (34), immunotherapy (54), or both (20). Models of all previously reported major melanoma subtypes (BRAF, NRAS, NF1, KIT, and WT/WT/WT) were identified. Multiple minor melanoma subtypes were also recapitulated, including melanomas with multiple activating mutations in the MAPK-signaling pathway and chromatin-remodeling gene mutations. These well-characterized melanoma PDXs and cell lines can be used not only as reagents for a large array of biological studies but also as pre-clinical models to facilitate drug development.

A Comprehensive Patient-Derived Xenograft Collection Representing the Heterogeneity of Melanoma.

Therapy of advanced melanoma is changing dramatically. Following mutational and biological subclassification of this heterogeneous cancer, several targeted and immune therapies were approved and increased survival significantly. To facilitate further advancements through pre-clinical in vivo modeling, we have established 459 patient-derived xenografts (PDX) and live tissue samples from 384 patients representing the full spectrum of clinical, therapeutic, mutational, and biological heterogeneity of melanoma. PDX have been characterized using targeted sequencing and protein arrays and are clinically annotated. This exhaustive live tissue resource includes PDX from 57 samples resistant to targeted therapy, 61 samples from responders and non-responders to immune checkpoint blockade, and 31 samples from brain metastasis. Uveal, mucosal, and acral subtypes are represented as well. We show examples of pre-clinical trials that highlight how the PDX collection can be used to develop and optimize precision therapies, biomarkers of response, and the targeting of rare genetic subgroups.