Ionizing radiation is a potent carcinogen, inducing cancer through DNA damage. The signatures of mutations arising in human tissues following in vivo exposure to ionizing radiation have not been documented. Here, we searched for signatures of ionizing radiation in 12 radiation-associated second malignancies of different tumour types. Two signatures of somatic mutation characterize ionizing radiation exposure irrespective of tumour type. Compared with 319 radiation-naive tumours, radiation-associated tumours carry a median extra 201 deletions genome-wide, sized 1-100 base pairs often with microhomology at the junction. Unlike deletions of radiation-naive tumours, these show no variation in density across the genome or correlation with sequence context, replication timing or chromatin structure. Furthermore, we observe a significant increase in balanced inversions in radiation-associated tumours. Both small deletions and inversions generate driver mutations. Thus, ionizing radiation generates distinctive mutational signatures that explain its carcinogenic potential.

BRAF (v-raf murine sarcoma viral oncogene homolog B) inhibitors elicit a transient anti-tumor response in ∼ 80% of BRAF(V600)-mutant melanoma patients that almost uniformly precedes the emergence of resistance. Here we used a mouse model of melanoma in which melanocyte-specific expression of Braf(V618E) (analogous to the human BRAF(V600E) mutation) led to the development of skin hyperpigmentation and nevi, as well as melanoma formation with incomplete penetrance. Sleeping Beauty insertional mutagenesis in this model led to accelerated and fully penetrant melanomagenesis and synchronous tumor formation. Treatment of Braf(V618E) transposon mice with the BRAF inhibitor PLX4720 resulted in tumor regression followed by relapse. Analysis of transposon insertions identified eight genes including Braf, Mitf, and ERas (ES-cell expressed Ras) as candidate resistance genes. Expression of ERAS in human melanoma cell lines conferred resistance to PLX4720 and induced hyperphosphorylation of AKT (v-akt murine thymoma viral oncogene homolog 1), a phenotype reverted by combinatorial treatment with PLX4720 and the AKT inhibitor MK2206. We show that ERAS expression elicits a prosurvival signal associated with phosphorylation/inactivation of BAD, and that the resistance of hepatocyte growth factor-treated human melanoma cells to PLX4720 can be reverted by treatment with the BAD-like BH3 mimetic ABT-737. Thus, we define a role for the AKT/BAD pathway in resistance to BRAF inhibition and illustrate an in vivo approach for finding drug resistance genes.

Adamantinomatous craniopharyngiomas (ACPs) are clinically challenging tumours, the majority of which have activating mutations in CTNNB1. They are histologically complex, showing cystic and solid components, the latter comprised of different morphological cell types (e.g. β-catenin-accumulating cluster cells and palisading epithelium), surrounded by a florid glial reaction with immune cells. Here, we have carried out RNA sequencing on 18 ACP samples and integrated these data with an existing ACP transcriptomic dataset. No studies so far have examined the patterns of gene expression within the different cellular compartments of the tumour. To achieve this goal, we have combined laser capture microdissection with computational analyses to reveal groups of genes that are associated with either epithelial tumour cells (clusters and palisading epithelium), glial tissue or immune infiltrate. We use these human ACP molecular signatures and RNA-Seq data from two ACP mouse models to reveal that cell clusters are molecularly analogous to the enamel knot, a critical signalling centre controlling normal tooth morphogenesis. Supporting this finding, we show that human cluster cells express high levels of several members of the FGF, TGFB and BMP families of secreted factors, which signal to neighbouring cells as evidenced by immunostaining against the phosphorylated proteins pERK1/2, pSMAD3 and pSMAD1/5/9 in both human and mouse ACP. We reveal that inhibiting the MAPK/ERK pathway with trametinib, a clinically approved MEK inhibitor, results in reduced proliferation and increased apoptosis in explant cultures of human and mouse ACP. Finally, we analyse a prominent molecular signature in the glial reactive tissue to characterise the inflammatory microenvironment and uncover the activation of inflammasomes in human ACP. We validate these results by immunostaining against immune cell markers, cytokine ELISA and proteome analysis in both solid tumour and cystic fluid from ACP patients. Our data support a new molecular paradigm for understanding ACP tumorigenesis as an aberrant mimic of natural tooth development and opens new therapeutic opportunities by revealing the activation of the MAPK/ERK and inflammasome pathways in human ACP.

Liver cirrhosis is associated with decreased hepatic cytochrome P4503A (CYP3A) activity but the pathogenesis of this phenomenon is not well elucidated. In this study, we examined if certain microRNAs (miRNA) are associated with decreased hepatic CYP3A activity in cirrhosis.

Obeticholic acid (OCA), a potent farnesoid X receptor agonist, was studied as monotherapy in an international, randomized, double-blind, placebo-controlled phase 2 study in patients with primary biliary cholangitis who were then followed for up to 6 years. The goals of the study were to assess the benefit of OCA in the absence of ursodeoxycholic acid, which is relevant for patients who are intolerant of ursodeoxycholic acid and at higher risk of disease progression. Patients were randomized and dosed with placebo (n = 23), OCA 10 mg (n = 20), or OCA 50 mg (n = 16) given as monotherapy once daily for 3 months (1 randomized patient withdrew prior to dosing). The primary endpoint was the percent change in alkaline phosphatase from baseline to the end of the double-blind phase of the study. Secondary and exploratory endpoints included change from baseline to month 3/early termination in markers of cholestasis, hepatocellular injury, and farnesoid X receptor activation. Efficacy and safety continue to be monitored through an ongoing 6-year open-label extension (N = 28). Alkaline phosphatase was reduced in both OCA groups (median% [Q1, Q3], OCA 10 mg -53.9% [-62.5, -29.3], OCA 50 mg -37.2% [-54.8, -24.6]) compared to placebo (-0.8% [-6.4, 8.7]; P < 0.0001) at the end of the study, with similar reductions observed through 6 years of open-label extension treatment. OCA improved many secondary and exploratory endpoints (including γ-glutamyl transpeptidase, alanine aminotransferase, conjugated bilirubin, and immunoglobulin M). Pruritus was the most common adverse event; 15% (OCA 10 mg) and 38% (OCA 50 mg) discontinued due to pruritus.

The genetics of renal cancer is dominated by inactivation of the VHL tumour suppressor gene in clear cell carcinoma (ccRCC), the commonest histological subtype. A recent large-scale screen of ∼3,500 genes by PCR-based exon re-sequencing identified several new cancer genes in ccRCC including UTX (also known as KDM6A), JARID1C (also known as KDM5C) and SETD2 (ref. 2). These genes encode enzymes that demethylate (UTX, JARID1C) or methylate (SETD2) key lysine residues of histone H3. Modification of the methylation state of these lysine residues of histone H3 regulates chromatin structure and is implicated in transcriptional control. However, together these mutations are present in fewer than 15% of ccRCC, suggesting the existence of additional, currently unidentified cancer genes. Here, we have sequenced the protein coding exome in a series of primary ccRCC and report the identification of the SWI/SNF chromatin remodelling complex gene PBRM1 (ref. 4) as a second major ccRCC cancer gene, with truncating mutations in 41% (92/227) of cases. These data further elucidate the somatic genetic architecture of ccRCC and emphasize the marked contribution of aberrant chromatin biology.

The emergence of high-throughput, next-generation sequencing technologies has dramatically altered the way we assess genomes in population genetics and in cancer genomics. Currently, there are four commonly used whole-genome sequencing platforms on the market: Illumina's HiSeq2000, Life Technologies' SOLiD 4 and its completely redesigned 5500xl SOLiD, and Complete Genomics' technology. A number of earlier studies have compared a subset of those sequencing platforms or compared those platforms with Sanger sequencing, which is prohibitively expensive for whole genome studies. Here we present a detailed comparison of the performance of all currently available whole genome sequencing platforms, especially regarding their ability to call SNVs and to evenly cover the genome and specific genomic regions. Unlike earlier studies, we base our comparison on four different samples, allowing us to assess the between-sample variation of the platforms. We find a pronounced GC bias in GC-rich regions for Life Technologies' platforms, with Complete Genomics performing best here, while we see the least bias in GC-poor regions for HiSeq2000 and 5500xl. HiSeq2000 gives the most uniform coverage and displays the least sample-to-sample variation. In contrast, Complete Genomics exhibits by far the smallest fraction of bases not covered, while the SOLiD platforms reveal remarkable shortcomings, especially in covering CpG islands. When comparing the performance of the four platforms for calling SNPs, HiSeq2000 and Complete Genomics achieve the highest sensitivity, while the SOLiD platforms show the lowest false positive rate. Finally, we find that integrating sequencing data from different platforms offers the potential to combine the strengths of different technologies. In summary, our results detail the strengths and weaknesses of all four whole-genome sequencing platforms. It indicates application areas that call for a specific sequencing platform and disallow other platforms. This helps to identify the proper sequencing platform for whole genome studies with different application scopes.

Exposure to extracellular 5'-adenosine triphosphate (ATP) is known to induce membrane blebbing. In this study, we investigated the subcellular distribution of the cytoskeletal adaptor protein paxillin in primary bovine osteoblasts upon stimulation with ATP. Cells expressing a fusion protein of green fluorescent protein (GFP) and paxillin were followed by time-lapse video-microscopy after stimulation with 100 μM ATP. Within 100 s, GFP-paxillin became incorporated in numerous de novo formed focal aggregates localized at the cell periphery. The assembly of individual paxillin-containing aggregates occurred with a mean half-life time of <60 s, whereas their disassembly lasted twice as long. Despite the ongoing presence of ATP, the formation of paxillin aggregates was self-limiting within 25 min. Paxillin clustering was preceded by a transient rise in cytoplasmic calcium transients, which peaked already 20 s after adding ATP. The high mobility of paxillin was confirmed by measuring the dissociation rate of GFP-paxillin at mature focal adhesions, demonstrating the presence of a highly mobile fraction with a mean recovery half-life of 8.2 ± 1.2 s, followed by a slower phase (53 ± 20 s). Thus, both the exchange of paxillin at mature focal adhesions and the increase in intracellular calcium concentrations upon ATP stimulation are very rapid processes, which override the time course of ATP-induced paxillin membrane clustering by one to two orders of magnitude. Our data demonstrate that the transient recruitment of paxillin in membrane protuberances is based on the high intracytoplasmic mobility of unbound paxillin molecules and their rapid focal accumulation.

In a systematic sequencing screen of the coding exons of the X chromosome in 250 families with X-linked mental retardation (XLMR), we identified two nonsense mutations and one consensus splice-site mutation in the AP1S2 gene on Xp22 in three families. Affected individuals in these families showed mild-to-profound mental retardation. Other features included hypotonia early in life and delay in walking. AP1S2 encodes an adaptin protein that constitutes part of the adaptor protein complex found at the cytoplasmic face of coated vesicles located at the Golgi complex. The complex mediates the recruitment of clathrin to the vesicle membrane. Aberrant endocytic processing through disruption of adaptor protein complexes is likely to result from the AP1S2 mutations identified in the three XLMR-affected families, and such defects may plausibly cause abnormal synaptic development and function. AP1S2 is the first reported XLMR gene that encodes a protein directly involved in the assembly of endocytic vesicles.

Transmissible cancers are clonal lineages that spread through populations via contagious cancer cells. Although rare in nature, two facial tumor clones affect Tasmanian devils. Here we perform comparative genetic and functional characterization of these lineages. The two cancers have similar patterns of mutation and show no evidence of exposure to exogenous mutagens or viruses. Genes encoding PDGF receptors have copy number gains and are present on extrachromosomal double minutes. Drug screening indicates causative roles for receptor tyrosine kinases and sensitivity to inhibitors of DNA repair. Y chromosome loss from a male clone infecting a female host suggests immunoediting. These results imply that Tasmanian devils may have inherent susceptibility to transmissible cancers and present a suite of therapeutic compounds for use in conservation.

Recent updating of the World Health Organization (WHO) classification of central nervous system (CNS) tumors in 2016 demonstrates the first organized effort to restructure brain tumor classification by incorporating histomorphologic features with recurrent molecular alterations. Revised CNS tumor diagnostic criteria also attempt to reduce interobserver variability of histological interpretation and provide more accurate stratification related to clinical outcome. As an example, diffuse gliomas (WHO grades II-IV) are now molecularly stratified based upon isocitrate dehydrogenase 1 or 2 (IDH) mutational status, with gliomas of WHO grades II and III being substratified according to 1p/19q codeletion status. For now, grading of diffuse gliomas is still dependent upon histological parameters. Independent of WHO classification criteria, multidimensional scaling analysis of molecular signatures for diffuse gliomas from The Cancer Genome Atlas (TCGA) has identified distinct molecular subgroups, and allows for their visualization in 2-dimensional (2D) space. Using the web-based platform Oncoscape as a tool, we applied multidimensional scaling-derived molecular groups to the 2D visualization of the 2016 WHO classification of diffuse gliomas. Here we show that molecular multidimensional scaling of TCGA data provides 2D clustering that represents the 2016 WHO classification of diffuse gliomas. Additionally, we used this platform to successfully identify and define novel copy-number alteration-based molecular subtypes, which are independent of WHO grading, as well as predictive of clinical outcome. The prognostic utility of these molecular subtypes was further validated using an independent data set of the German Glioma Network prospective glioblastoma patient cohort.

Combinations of anti-cancer drugs can overcome resistance and provide new treatments1,2. The number of possible drug combinations vastly exceeds what could be tested clinically. Efforts to systematically identify active combinations and the tissues and molecular contexts in which they are most effective could accelerate the development of combination treatments. Here we evaluate the potency and efficacy of 2,025 clinically relevant two-drug combinations, generating a dataset encompassing 125 molecularly characterized breast, colorectal and pancreatic cancer cell lines. We show that synergy between drugs is rare and highly context-dependent, and that combinations of targeted agents are most likely to be synergistic. We incorporate multi-omic molecular features to identify combination biomarkers and specify synergistic drug combinations and their active contexts, including in basal-like breast cancer, and microsatellite-stable or KRAS-mutant colon cancer. Our results show that irinotecan and CHEK1 inhibition have synergistic effects in microsatellite-stable or KRAS-TP53 double-mutant colon cancer cells, leading to apoptosis and suppression of tumour xenograft growth. This study identifies clinically relevant effective drug combinations in distinct molecular subpopulations and is a resource to guide rational efforts to develop combinatorial drug treatments.

Loss-of-function mutations in the progranulin gene (GRN) are one of the major causes of familial frontotemporal lobar degeneration. Our objective was to determine the rates and trajectories of lobar cortical atrophy from longitudinal structural magnetic resonance imaging in both asymptomatic and symptomatic GRN mutation carriers. Individuals in this study were from the ADRC and LEFFTDS studies at the Mayo Clinic. We identified 13 GRN mutation carriers (8 asymptomatic, 5 symptomatic) and noncarriers (n = 10) who had at least 2 serial T1-weighted structural magnetic resonance images and were followed annually with a median of 3 years (range 1.0-9.8 years). Longitudinal changes in lobar cortical volume were analyzed using the tensor-based morphometry with symmetric normalization (TBM-SyN) algorithm. Linear mixed-effect models were used to model cortical volume change over time among 3 groups. The annual rates of frontal (p < 0.05) and parietal (p < 0.01) lobe cortical atrophy were higher in asymptomatic GRN mutation carriers than noncarriers. The symptomatic GRN mutation carriers also had increased rates of atrophy in the frontal (p < 0.01) and parietal lobe (p < 0.01) cortices than noncarriers. In addition, greater rates of cortical atrophy were observed in the temporal lobe cortices of symptomatic GRN mutation carriers than noncarriers (p < 0.001). We found that a decline in frontal and parietal lobar cortical volume occurs in asymptomatic GRN mutation carriers and continues in the symptomatic GRN mutation carriers, whereas an increased rate of temporal lobe cortical atrophy is observed only in symptomatic GRN mutation carriers. This sequential pattern of cortical involvement in GRN mutation carriers has important implications for using imaging biomarkers of neurodegeneration as an outcome measure in potential treatment trials involving GRN mutation carriers.

Protein S is a cofactor in the tissue factor pathway inhibitor (TFPI) anticoagulant pathway. It enhances TFPIα-mediated inhibition of factor (F)Xa activity and generation. The enhancement is dependent on a TFPIα-protein S interaction involving TFPIα Kunitz 3 and protein S laminin G-type (LG)-1. C4b binding protein (C4BP), which binds to protein S LG1, almost completely abolishes its TFPI cofactor function. However, neither the amino acids involved in TFPIα enhancement nor the mechanisms underlying the reduced TFPI cofactor function of C4BP-bound protein S are known. To screen for functionally important regions within protein S LG1, we generated 7 variants with inserted N-linked glycosylation attachment sites. Protein S D253T and Q427N/K429T displayed severely reduced TFPI cofactor function while showing normal activated protein C (APC) cofactor function and C4BP binding. Based on these results, we designed 4 protein S variants in which 4 to 6 surface-exposed charged residues were substituted for alanine. One variant, protein S K255A/E257A/D287A/R410A/K423A/E424A, exhibited either abolished or severely reduced TFPI cofactor function in plasma and FXa inhibition assays, both in the presence or absence of FV-short, but retained normal APC cofactor function and high-affinity C4BP binding. The C4BP β-chain was expressed to determine the mechanisms behind the reduced TFPI cofactor function of C4BP-bound protein S. Like C4BP-bound protein S, C4BP β-chain-bound protein S had severely reduced TFPI cofactor function. These results show that protein S Lys255, Glu257, Asp287, Arg410, Lys423, and Glu424 are critical for protein S-mediated enhancement of TFPIα and that binding of the C4BP β-chain blocks this function.

Alzheimer's disease prevalence will reach epidemic proportions in coming decades. There is a need for impactful educational materials to help patients, families, medical practitioners, and policy makers understand the nature and impact of the disease. Defining an effective workflow to create such models from existing segmentation tools will be a valuable contribution in creating these patient-specific models.

The incidence of acute myeloid leukaemia (AML) increases with age and mortality exceeds 90% when diagnosed after age 65. Most cases arise without any detectable early symptoms and patients usually present with the acute complications of bone marrow failure1. The onset of such de novo AML cases is typically preceded by the accumulation of somatic mutations in preleukaemic haematopoietic stem and progenitor cells (HSPCs) that undergo clonal expansion2,3. However, recurrent AML mutations also accumulate in HSPCs during ageing of healthy individuals who do not develop AML, a phenomenon referred to as age-related clonal haematopoiesis (ARCH)4-8. Here we use deep sequencing to analyse genes that are recurrently mutated in AML to distinguish between individuals who have a high risk of developing AML and those with benign ARCH. We analysed peripheral blood cells from 95 individuals that were obtained on average 6.3 years before AML diagnosis (pre-AML group), together with 414 unselected age- and gender-matched individuals (control group). Pre-AML cases were distinct from controls and had more mutations per sample, higher variant allele frequencies, indicating greater clonal expansion, and showed enrichment of mutations in specific genes. Genetic parameters were used to derive a model that accurately predicted AML-free survival; this model was validated in an independent cohort of 29 pre-AML cases and 262 controls. Because AML is rare, we also developed an AML predictive model using a large electronic health record database that identified individuals at greater risk. Collectively our findings provide proof-of-concept that it is possible to discriminate ARCH from pre-AML many years before malignant transformation. This could in future enable earlier detection and monitoring, and may help to inform intervention.

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed "kataegis," was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

We have identified one frameshift mutation, one splice-site mutation, and two missense mutations in highly conserved residues in ZDHHC9 at Xq26.1 in 4 of 250 families with X-linked mental retardation (XLMR). In three of the families, the mental retardation phenotype is associated with a Marfanoid habitus, although none of the affected individuals meets the Ghent criteria for Marfan syndrome. ZDHHC9 is a palmitoyltransferase that catalyzes the posttranslational modification of NRAS and HRAS. The degree of palmitoylation determines the temporal and spatial location of these proteins in the plasma membrane and Golgi complex. The finding of mutations in ZDHHC9 suggests that alterations in the concentrations and cellular distribution of target proteins are sufficient to cause disease. This is the first XLMR gene to be reported that encodes a posttranslational modification enzyme, palmitoyltransferase. Furthermore, now that the first palmitoyltransferase that causes mental retardation has been identified, defects in other palmitoylation transferases become good candidates for causing other mental retardation syndromes.

Multivalent second-generation TRAIL-R2 agonists are currently in late preclinical development and early clinical trials. Herein, we use a representative second-generation agent, MEDI3039, to address two major clinical challenges facing these agents: lack of predictive biomarkers to enable patient selection and emergence of resistance. Genome-wide CRISPR knockout screens were notable for the lack of resistance mechanisms beyond the canonical TRAIL-R2 pathway (caspase-8, FADD, BID) as well as p53 and BAX in TP53 wild-type models, whereas a CRISPR activatory screen identified cell death inhibitors MCL-1 and BCL-XL as mechanisms to suppress MEDI3039-induced cell death. High-throughput drug screening failed to identify genomic alterations associated with response to MEDI3039; however, transcriptomics analysis revealed striking association between MEDI3039 sensitivity and expression of core components of the extrinsic apoptotic pathway, most notably its main apoptotic effector caspase-8 in solid tumor cell lines. Further analyses of colorectal cell lines and patient-derived xenografts identified caspase-8 expression ratio to its endogenous regulator FLIP(L) as predictive of sensitivity to MEDI3039 in several major solid tumor types and a further subset indicated by caspase-8:MCL-1 ratio. Subsequent MEDI3039 combination screening of TRAIL-R2, caspase-8, FADD, and BID knockout models with 60 compounds with varying mechanisms of action identified two inhibitor of apoptosis proteins (IAP) that exhibited strong synergy with MEDI3039 that could reverse resistance only in BID-deleted models. In summary, we identify the ratios of caspase-8:FLIP(L) and caspase-8:MCL-1 as potential predictive biomarkers for second-generation TRAIL-R2 agonists and loss of key effectors such as FADD and caspase-8 as likely drivers of clinical resistance in solid tumors.

Transmission-blocking vaccines (TBVs) that target sexual stage parasite development could be an integral part of measures for malaria elimination. Pfs25 is a leading TBV candidate, and previous studies conducted in animals demonstrated an improvement of its functional immunogenicity after conjugation to EPA, a recombinant, detoxified ExoProtein A from Pseudomonas aeruginosa. In this report, we describe results of an open-label, dose-escalating Phase 1 trial to assess the safety and immunogenicity of Pfs25-EPA conjugates formulated with Alhydrogel®. Thirty malaria-naïve healthy adults received up to four doses of the conjugate vaccine, with 8, 16, or 47 μg of conjugated Pfs25 mass, at 0, 2, 4, and 10 months. Vaccinations were generally well tolerated. The majority of solicited adverse events were mild in severity with pain at the injection site the most common complaint. Anemia was the most common laboratory abnormality, but was considered possibly related to the study in only a minority of cases. No vaccine-related serious adverse events occurred. The peak geometric mean anti-Pfs25 antibody level in the highest dose group was 88 (95% CI 53, 147) μg/mL two weeks after the 4th vaccination, and declined to near baseline one year later. Antibody avidity increased over successive vaccinations. Transmission blocking activity demonstrated in a standard membrane feeding assay (SMFA) also increased from the second to the third dose, and correlated with antibody titer and, after the final dose, with antibody avidity. These results support the further evaluation of Pfs25-EPA/Alhydrogel® in a malaria-endemic population.