An interactive resource to identify cancer genetic and lineage dependencies targeted by small molecules.

The high rate of clinical response to protein-kinase-targeting drugs matched to cancer patients with specific genomic alterations has prompted efforts to use cancer cell line (CCL) profiling to identify additional biomarkers of small-molecule sensitivities. We have quantitatively measured the sensitivity of 242 genomically characterized CCLs to an Informer Set of 354 small molecules that target many nodes in cell circuitry, uncovering protein dependencies that: (1) associate with specific cancer-genomic alterations and (2) can be targeted by small molecules. We have created the Cancer Therapeutics Response Portal (http://www.broadinstitute.org/ctrp) to enable users to correlate genetic features to sensitivity in individual lineages and control for confounding factors of CCL profiling. We report a candidate dependency, associating activating mutations in the oncogene ß-catenin with sensitivity to the Bcl-2 family antagonist, navitoclax. The resource can be used to develop novel therapeutic hypotheses and to accelerate discovery of drugs matched to patients by their cancer genotype and lineage.

Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

The effect of external stimulation on functional networks in the aging healthy human brain.

Understanding the brain changes occurring during aging can provide new insights for developing treatments that alleviate or reverse cognitive decline. Neurostimulation techniques have emerged as potential treatments for brain disorders and to improve cognitive functions. Nevertheless, given the ethical restrictions of neurostimulation approaches, in silico perturbation protocols based on causal whole-brain models are fundamental to gaining a mechanistic understanding of brain dynamics. Furthermore, this strategy could serve to identify neurophysiological biomarkers differentiating between age groups through an exhaustive exploration of the global effect of all possible local perturbations. Here, we used a resting-state fMRI dataset divided into middle-aged (N =310, <65 years) and older adults (N =310, $\geq $65) to characterize brain states in each group as a probabilistic metastable substate (PMS) space. We showed that the older group exhibited a reduced capability to access a metastable substate that overlaps with the rich club. Then, we fitted the PMS to a whole-brain model and applied in silico stimulations in each node to force transitions from the brain states of the older- to the middle-aged group. We found that the precuneus was the best stimulation target. Overall, these findings could have important implications for designing neurostimulation interventions for reversing the effects of aging on whole-brain dynamics.

Whole-Brain Dynamics in Aging: Disruptions in Functional Connectivity and the Role of the Rich Club.

Normal aging causes disruptions in the brain that can lead to cognitive decline. Resting-state functional magnetic resonance imaging studies have found significant age-related alterations in functional connectivity across various networks. Nevertheless, most of the studies have focused mainly on static functional connectivity. Studying the dynamics of resting-state brain activity across the whole-brain functional network can provide a better characterization of age-related changes. Here, we employed two data-driven whole-brain approaches based on the phase synchronization of blood-oxygen-level-dependent signals to analyze resting-state fMRI data from 620 subjects divided into two groups (middle-age group (n = 310); age range, 50-64 years versus older group (n = 310); age range, 65-91 years). Applying the intrinsic-ignition framework to assess the effect of spontaneous local activation events on local-global integration, we found that the older group showed higher intrinsic ignition across the whole-brain functional network, but lower metastability. Using Leading Eigenvector Dynamics Analysis, we found that the older group showed reduced ability to access a metastable substate that closely overlaps with the so-called rich club. These findings suggest that functional whole-brain dynamics are altered in aging, probably due to a deficiency in a metastable substate that is key for efficient global communication in the brain.

Obesity-associated deficits in inhibitory control are phenocopied to mice through gut microbiota changes in one-carbon and aromatic amino acids metabolic pathways.

BACKGROUND: Inhibitory control (IC) is critical to keep long-term goals in everyday life. Bidirectional relationships between IC deficits and obesity are behind unhealthy eating and physical exercise habits. METHODS: We studied gut microbiome composition and functionality, and plasma and faecal metabolomics in association with cognitive tests evaluating inhibitory control (Stroop test) and brain structure in a discovery (n=156), both cross-sectionally and longitudinally, and in an independent replication cohort (n=970). Faecal microbiota transplantation (FMT) in mice evaluated the impact on reversal learning and medial prefrontal cortex (mPFC) transcriptomics. RESULTS: An interplay among IC, brain structure (in humans) and mPFC transcriptomics (in mice), plasma/faecal metabolomics and the gut metagenome was found. Obesity-dependent alterations in one-carbon metabolism, tryptophan and histidine pathways were associated with IC in the two independent cohorts. Bacterial functions linked to one-carbon metabolism (thyX,dut, exodeoxyribonuclease V), and the anterior cingulate cortex volume were associated with IC, cross-sectionally and longitudinally. FMT from individuals with obesity led to alterations in mice reversal learning. In an independent FMT experiment, human donor's bacterial functions related to IC deficits were associated with mPFC expression of one-carbon metabolism-related genes of recipient's mice. CONCLUSION: These results highlight the importance of targeting obesity-related impulsive behaviour through the induction of gut microbiota shifts.

Germline variants and somatic mutation signatures of breast cancer across populations of African and European ancestry in the US and Nigeria.

Somatic mutation signatures may represent footprints of genetic and environmental exposures that cause different cancer. Few studies have comprehensively examined their association with germline variants, and none in an indigenous African population. SomaticSignatures was employed to extract mutation signatures based on whole-genome or whole-exome sequencing data from female patients with breast cancer (TCGA, training set, n = 1,011; Nigerian samples, validation set, n = 170), and to estimate contributions of signatures in each sample. Association between somatic signatures and common single nucleotide polymorphisms (SNPs) or rare deleterious variants were examined using linear regression. Nine stable signatures were inferred, and four signatures (APOBEC C>T, APOBEC C>G, aging and homologous recombination deficiency) were highly similar to known COSMIC signatures and explained the majority (60-85%) of signature contributions. There were significant heritable components associated with APOBEC C>T signature (h2 = 0.575, p = 0.010) and the combined APOBEC signatures (h2 = 0.432, p = 0.042). In TCGA dataset, seven common SNPs within or near GNB5 were significantly associated with an increased proportion (beta = 0.33, 95% CI = 0.21-0.45) of APOBEC signature contribution at genome-wide significance, while rare germline mutations in MTCL1 was also significantly associated with a higher contribution of this signature (p = 6.1 × 10-6 ). This is the first study to identify associations between germline variants and mutational patterns in breast cancer across diverse populations and geography. The findings provide evidence to substantiate causal links between germline genetic risk variants and carcinogenesis.

Nuclear factor I/B is an oncogene in small cell lung cancer.

Small cell lung cancer (SCLC) is an aggressive cancer often diagnosed after it has metastasized. Despite the need to better understand this disease, SCLC remains poorly characterized at the molecular and genomic levels. Using a genetically engineered mouse model of SCLC driven by conditional deletion of Trp53 and Rb1 in the lung, we identified several frequent, high-magnitude focal DNA copy number alterations in SCLC. We uncovered amplification of a novel, oncogenic transcription factor, Nuclear factor I/B (Nfib), in the mouse SCLC model and in human SCLC. Functional studies indicate that NFIB regulates cell viability and proliferation during transformation.

Melanoma genome sequencing reveals frequent PREX2 mutations.

Melanoma is notable for its metastatic propensity, lethality in the advanced setting and association with ultraviolet exposure early in life. To obtain a comprehensive genomic view of melanoma in humans, we sequenced the genomes of 25 metastatic melanomas and matched germline DNA. A wide range of point mutation rates was observed: lowest in melanomas whose primaries arose on non-ultraviolet-exposed hairless skin of the extremities (3 and 14 per megabase (Mb) of genome), intermediate in those originating from hair-bearing skin of the trunk (5-55 per Mb), and highest in a patient with a documented history of chronic sun exposure (111 per Mb). Analysis of whole-genome sequence data identified PREX2 (phosphatidylinositol-3,4,5-trisphosphate-dependent Rac exchange factor 2)--a PTEN-interacting protein and negative regulator of PTEN in breast cancer--as a significantly mutated gene with a mutation frequency of approximately 14% in an independent extension cohort of 107 human melanomas. PREX2 mutations are biologically relevant, as ectopic expression of mutant PREX2 accelerated tumour formation of immortalized human melanocytes in vivo. Thus, whole-genome sequencing of human melanoma tumours revealed genomic evidence of ultraviolet pathogenesis and discovered a new recurrently mutated gene in melanoma.

The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity.

The systematic translation of cancer genomic data into knowledge of tumour biology and therapeutic possibilities remains challenging. Such efforts should be greatly aided by robust preclinical model systems that reflect the genomic diversity of human cancers and for which detailed genetic and pharmacological annotation is available. Here we describe the Cancer Cell Line Encyclopedia (CCLE): a compilation of gene expression, chromosomal copy number and massively parallel sequencing data from 947 human cancer cell lines. When coupled with pharmacological profiles for 24 anticancer drugs across 479 of the cell lines, this collection allowed identification of genetic, lineage, and gene-expression-based predictors of drug sensitivity. In addition to known predictors, we found that plasma cell lineage correlated with sensitivity to IGF1 receptor inhibitors; AHR expression was associated with MEK inhibitor efficacy in NRAS-mutant lines; and SLFN11 expression predicted sensitivity to topoisomerase inhibitors. Together, our results indicate that large, annotated cell-line collections may help to enable preclinical stratification schemata for anticancer agents. The generation of genetic predictions of drug response in the preclinical setting and their incorporation into cancer clinical trial design could speed the emergence of 'personalized' therapeutic regimens.

Functional genomics reveal that the serine synthesis pathway is essential in breast cancer.

Cancer cells adapt their metabolic processes to drive macromolecular biosynthesis for rapid cell growth and proliferation. RNA interference (RNAi)-based loss-of-function screening has proven powerful for the identification of new and interesting cancer targets, and recent studies have used this technology in vivo to identify novel tumour suppressor genes. Here we developed a method for identifying novel cancer targets via negative-selection RNAi screening using a human breast cancer xenograft model at an orthotopic site in the mouse. Using this method, we screened a set of metabolic genes associated with aggressive breast cancer and stemness to identify those required for in vivo tumorigenesis. Among the genes identified, phosphoglycerate dehydrogenase (PHGDH) is in a genomic region of recurrent copy number gain in breast cancer and PHGDH protein levels are elevated in 70% of oestrogen receptor (ER)-negative breast cancers. PHGDH catalyses the first step in the serine biosynthesis pathway, and breast cancer cells with high PHGDH expression have increased serine synthesis flux. Suppression of PHGDH in cell lines with elevated PHGDH expression, but not in those without, causes a strong decrease in cell proliferation and a reduction in serine synthesis. We find that PHGDH suppression does not affect intracellular serine levels, but causes a drop in the levels of α-ketoglutarate, another output of the pathway and a tricarboxylic acid (TCA) cycle intermediate. In cells with high PHGDH expression, the serine synthesis pathway contributes approximately 50% of the total anaplerotic flux of glutamine into the TCA cycle. These results reveal that certain breast cancers are dependent upon increased serine pathway flux caused by PHGDH overexpression and demonstrate the utility of in vivo negative-selection RNAi screens for finding potential anticancer targets.

COT drives resistance to RAF inhibition through MAP kinase pathway reactivation.

Oncogenic mutations in the serine/threonine kinase B-RAF (also known as BRAF) are found in 50-70% of malignant melanomas. Pre-clinical studies have demonstrated that the B-RAF(V600E) mutation predicts a dependency on the mitogen-activated protein kinase (MAPK) signalling cascade in melanoma-an observation that has been validated by the success of RAF and MEK inhibitors in clinical trials. However, clinical responses to targeted anticancer therapeutics are frequently confounded by de novo or acquired resistance. Identification of resistance mechanisms in a manner that elucidates alternative 'druggable' targets may inform effective long-term treatment strategies. Here we expressed ∼600 kinase and kinase-related open reading frames (ORFs) in parallel to interrogate resistance to a selective RAF kinase inhibitor. We identified MAP3K8 (the gene encoding COT/Tpl2) as a MAPK pathway agonist that drives resistance to RAF inhibition in B-RAF(V600E) cell lines. COT activates ERK primarily through MEK-dependent mechanisms that do not require RAF signalling. Moreover, COT expression is associated with de novo resistance in B-RAF(V600E) cultured cell lines and acquired resistance in melanoma cells and tissue obtained from relapsing patients following treatment with MEK or RAF inhibitors. We further identify combinatorial MAPK pathway inhibition or targeting of COT kinase activity as possible therapeutic strategies for reducing MAPK pathway activation in this setting. Together, these results provide new insights into resistance mechanisms involving the MAPK pathway and articulate an integrative approach through which high-throughput functional screens may inform the development of novel therapeutic strategies.

The landscape of somatic copy-number alteration across human cancers.

A powerful way to discover key genes with causal roles in oncogenesis is to identify genomic regions that undergo frequent alteration in human cancers. Here we present high-resolution analyses of somatic copy-number alterations (SCNAs) from 3,131 cancer specimens, belonging largely to 26 histological types. We identify 158 regions of focal SCNA that are altered at significant frequency across several cancer types, of which 122 cannot be explained by the presence of a known cancer target gene located within these regions. Several gene families are enriched among these regions of focal SCNA, including the BCL2 family of apoptosis regulators and the NF-kappaBeta pathway. We show that cancer cells containing amplifications surrounding the MCL1 and BCL2L1 anti-apoptotic genes depend on the expression of these genes for survival. Finally, we demonstrate that a large majority of SCNAs identified in individual cancer types are present in several cancer types.

CDYL bridges REST and histone methyltransferases for gene repression and suppression of cellular transformation.

The neuronal gene repressor REST/NRSF recruits corepressors, including CoREST, to modify histones and repress transcription. REST also functions as a tumor suppressor, but the mechanism remains unclear. We identified chromodomain on Y-like (CDYL) as a REST corepressor that physically bridges REST and the histone methylase G9a to repress transcription. Importantly, RNAi knockdown of REST, CDYL, and G9a, but not CoREST, induced oncogenic transformation of immortalized primary human cells and derepression of the proto-oncogene TrkC. Significantly, transgenic expression of TrkC also induced transformation. This implicates CDYL-G9a, but not CoREST, in REST suppression of transformation, possibly by oncogene repression. CDYL knockdown also augments transformation in a cell culture model of cervical cancer, where loss of heterozygosity of the CDYL locus occurs. These findings demonstrate molecular strategies by which REST carries out distinct biological functions via different corepressors and provide critical insights into the role of histone-modifying complexes in regulating cellular transformation.

High-throughput oncogene mutation profiling in human cancer.

Systematic efforts are underway to decipher the genetic changes associated with tumor initiation and progression. However, widespread clinical application of this information is hampered by an inability to identify critical genetic events across the spectrum of human tumors with adequate sensitivity and scalability. Here, we have adapted high-throughput genotyping to query 238 known oncogene mutations across 1,000 human tumor samples. This approach established robust mutation distributions spanning 17 cancer types. Of 17 oncogenes analyzed, we found 14 to be mutated at least once, and 298 (30%) samples carried at least one mutation. Moreover, we identified previously unrecognized oncogene mutations in several tumor types and observed an unexpectedly high number of co-occurring mutations. These results offer a new dimension in tumor genetics, where mutations involving multiple cancer genes may be interrogated simultaneously and in 'real time' to guide cancer classification and rational therapeutic intervention.

CDK8 is a colorectal cancer oncogene that regulates beta-catenin activity.

Aberrant activation of the canonical WNT/beta-catenin pathway occurs in almost all colorectal cancers and contributes to their growth, invasion and survival. Although dysregulated beta-catenin activity drives colon tumorigenesis, further genetic perturbations are required to elaborate full malignant transformation. To identify genes that both modulate beta-catenin activity and are essential for colon cancer cell proliferation, we conducted two loss-of-function screens in human colon cancer cells and compared genes identified in these screens with an analysis of copy number alterations in colon cancer specimens. One of these genes, CDK8, which encodes a member of the mediator complex, is located at 13q12.13, a region of recurrent copy number gain in a substantial fraction of colon cancers. Here we show that the suppression of CDK8 expression inhibits proliferation in colon cancer cells characterized by high levels of CDK8 and beta-catenin hyperactivity. CDK8 kinase activity was necessary for beta-catenin-driven transformation and for expression of several beta-catenin transcriptional targets. Together these observations suggest that therapeutic interventions targeting CDK8 may confer a clinical benefit in beta-catenin-driven malignancies.

Systematic investigation of genetic vulnerabilities across cancer cell lines reveals lineage-specific dependencies in ovarian cancer.

A comprehensive understanding of the molecular vulnerabilities of every type of cancer will provide a powerful roadmap to guide therapeutic approaches. Efforts such as The Cancer Genome Atlas Project will identify genes with aberrant copy number, sequence, or expression in various cancer types, providing a survey of the genes that may have a causal role in cancer. A complementary approach is to perform systematic loss-of-function studies to identify essential genes in particular cancer cell types. We have begun a systematic effort, termed Project Achilles, aimed at identifying genetic vulnerabilities across large numbers of cancer cell lines. Here, we report the assessment of the essentiality of 11,194 genes in 102 human cancer cell lines. We show that the integration of these functional data with information derived from surveying cancer genomes pinpoints known and previously undescribed lineage-specific dependencies across a wide spectrum of cancers. In particular, we found 54 genes that are specifically essential for the proliferation and viability of ovarian cancer cells and also amplified in primary tumors or differentially overexpressed in ovarian cancer cell lines. One such gene, PAX8, is focally amplified in 16% of high-grade serous ovarian cancers and expressed at higher levels in ovarian tumors. Suppression of PAX8 selectively induces apoptotic cell death of ovarian cancer cells. These results identify PAX8 as an ovarian lineage-specific dependency. More generally, these observations demonstrate that the integration of genome-scale functional and structural studies provides an efficient path to identify dependencies of specific cancer types on particular genes and pathways.

Stratified analyses refine association between TLR7 rare variants and severe COVID-19.

Despite extensive global research into genetic predisposition for severe COVID-19, knowledge on the role of rare host genetic variants and their relation to other risk factors remains limited. Here, 52 genes with prior etiological evidence were sequenced in 1,772 severe COVID-19 cases and 5,347 population-based controls from Spain/Italy. Rare deleterious TLR7 variants were present in 2.4% of young (<60 years) cases with no reported clinical risk factors (n=378), compared to 0.24% of controls (odds ratio (OR)=12.3, p=1.27x10-10). Incorporation of the results of either functional assays or protein modeling led to a pronounced increase in effect size (ORmax=46.5, p=1.74x10-15). Association signals for X-chromosomal TLR7 were also detected in the female-only subgroup, suggesting the existence of additional mechanisms beyond X-linked recessive inheritance in males. Additionally, supporting evidence was generated for a contribution to severe COVID-19 of the previously implicated genes IFNAR2, IFIH1 and TBK1. Our results refine the genetic contribution of rare TLR7 variants to severe COVID-19, and strengthen evidence for the etiological relevance of genes in the interferon signaling pathway.

Integrative analysis of the melanoma transcriptome.

Global studies of transcript structure and abundance in cancer cells enable the systematic discovery of aberrations that contribute to carcinogenesis, including gene fusions, alternative splice isoforms, and somatic mutations. We developed a systematic approach to characterize the spectrum of cancer-associated mRNA alterations through integration of transcriptomic and structural genomic data, and we applied this approach to generate new insights into melanoma biology. Using paired-end massively parallel sequencing of cDNA (RNA-seq) together with analyses of high-resolution chromosomal copy number data, we identified 11 novel melanoma gene fusions produced by underlying genomic rearrangements, as well as 12 novel readthrough transcripts. We mapped these chimeric transcripts to base-pair resolution and traced them to their genomic origins using matched chromosomal copy number information. We also used these data to discover and validate base-pair mutations that accumulated in these melanomas, revealing a surprisingly high rate of somatic mutation and lending support to the notion that point mutations constitute the major driver of melanoma progression. Taken together, these results may indicate new avenues for target discovery in melanoma, while also providing a template for large-scale transcriptome studies across many tumor types.