We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

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.

Liver receptor homologue 1 (LRH-1) is an orphan nuclear receptor that has been implicated in the progression of breast, pancreatic and colorectal cancer (CRC). To determine mechanisms underlying growth promotion by LRH-1 in CRC, we undertook global expression profiling following siRNA-mediated LRH-1 knockdown in HCT116 cells, which require LRH-1 for growth and in HT29 cells, in which LRH-1 does not regulate growth. Interestingly, expression of the cell cycle inhibitor p21 (CDKN1A) was regulated by LRH-1 in HCT116 cells. p21 regulation was not observed in HT29 cells, where p53 is mutated. p53 dependence for the regulation of p21 by LRH-1 was confirmed by p53 knockdown with siRNA, while LRH-1-regulation of p21 was not evident in HCT116 cells where p53 had been deleted. We demonstrate that LRH-1-mediated p21 regulation in HCT116 cells does not involve altered p53 protein or phosphorylation, and we show that LRH-1 inhibits p53 recruitment to the p21 promoter, likely through a mechanism involving chromatin remodelling. Our study suggests an important role for LRH-1 in the growth of CRC cells that retain wild-type p53.

The BioMart Community Portal (www.biomart.org) is a community-driven effort to provide a unified interface to biomedical databases that are distributed worldwide. The portal provides access to numerous database projects supported by 30 scientific organizations. It includes over 800 different biological datasets spanning genomics, proteomics, model organisms, cancer data, ontology information and more. All resources available through the portal are independently administered and funded by their host organizations. The BioMart data federation technology provides a unified interface to all the available data. The latest version of the portal comes with many new databases that have been created by our ever-growing community. It also comes with better support and extensibility for data analysis and visualization tools. A new addition to our toolbox, the enrichment analysis tool is now accessible through graphical and web service interface. The BioMart community portal averages over one million requests per day. Building on this level of service and the wealth of information that has become available, the BioMart Community Portal has introduced a new, more scalable and cheaper alternative to the large data stores maintained by specialized organizations.

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.

Cancer is driven by somatically acquired point mutations and chromosomal rearrangements, conventionally thought to accumulate gradually over time. Using next-generation sequencing, we characterize a phenomenon, which we term chromothripsis, whereby tens to hundreds of genomic rearrangements occur in a one-off cellular crisis. Rearrangements involving one or a few chromosomes crisscross back and forth across involved regions, generating frequent oscillations between two copy number states. These genomic hallmarks are highly improbable if rearrangements accumulate over time and instead imply that nearly all occur during a single cellular catastrophe. The stamp of chromothripsis can be seen in at least 2%-3% of all cancers, across many subtypes, and is present in ∼25% of bone cancers. We find that one, or indeed more than one, cancer-causing lesion can emerge out of the genomic crisis. This phenomenon has important implications for the origins of genomic remodeling and temporal emergence of cancer.

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.

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.

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.

Trained immunity describes the capacity of innate immune cells to develop heterologous memory in response to certain exogenous exposures. This phenomenon mediates, at least in part, the beneficial off-target effects of the BCG vaccine. Using an in vitro model of trained immunity, we show that BCG exposure induces a persistent change in active histone modifications, DNA methylation, transcription, and adenosine-to-inosine RNA modification in human monocytes. By profiling DNA methylation of circulating monocytes from infants in the MIS BAIR clinical trial, we identify a BCG-associated DNA methylation signature that persisted more than 12 months after neonatal BCG vaccination. Genes associated with this epigenetic signature are involved in viral response pathways, consistent with the reported off-target protection against viral infections in neonates, adults, and the elderly. Our findings indicate that the off-target effects of BCG in infants are accompanied by epigenetic remodeling of circulating monocytes that lasts more than 1 year.

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.

COSMIC (http://www.sanger.ac.uk/cosmic) curates comprehensive information on somatic mutations in human cancer. Release v48 (July 2010) describes over 136,000 coding mutations in almost 542,000 tumour samples; of the 18,490 genes documented, 4803 (26%) have one or more mutations. Full scientific literature curations are available on 83 major cancer genes and 49 fusion gene pairs (19 new cancer genes and 30 new fusion pairs this year) and this number is continually increasing. Key amongst these is TP53, now available through a collaboration with the IARC p53 database. In addition to data from the Cancer Genome Project (CGP) at the Sanger Institute, UK, and The Cancer Genome Atlas project (TCGA), large systematic screens are also now curated. Major website upgrades now make these data much more mineable, with many new selection filters and graphics. A Biomart is now available allowing more automated data mining and integration with other biological databases. Annotation of genomic features has become a significant focus; COSMIC has begun curating full-genome resequencing experiments, developing new web pages, export formats and graphics styles. With all genomic information recently updated to GRCh37, COSMIC integrates many diverse types of mutation information and is making much closer links with Ensembl and other data resources.

Genome editing by CRISPR-Cas9 technology allows large-scale screening of gene essentiality in cancer. A confounding factor when interpreting CRISPR-Cas9 screens is the high false-positive rate in detecting essential genes within copy number amplified regions of the genome. We have developed the computational tool CRISPRcleanR which is capable of identifying and correcting gene-independent responses to CRISPR-Cas9 targeting. CRISPRcleanR uses an unsupervised approach based on the segmentation of single-guide RNA fold change values across the genome, without making any assumption about the copy number status of the targeted genes.

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.