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.

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.

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.

We have identified truncating mutations in the human DLG3 (neuroendocrine dlg) gene in 4 of 329 families with moderate to severe X-linked mental retardation. DLG3 encodes synapse-associated protein 102 (SAP102), a member of the membrane-associated guanylate kinase protein family. Neuronal SAP102 is expressed during early brain development and is localized to the postsynaptic density of excitatory synapses. It is composed of three amino-terminal PDZ domains, an src homology domain, and a carboxyl-terminal guanylate kinase domain. The PDZ domains interact directly with the NR2 subunits of the NMDA glutamate receptor and with other proteins responsible for NMDA receptor localization, immobilization, and signaling. The mutations identified in this study all introduce premature stop codons within or before the third PDZ domain, and it is likely that this impairs the ability of SAP102 to interact with the NMDA receptor and/or other proteins involved in downstream NMDA receptor signaling pathways. NMDA receptors have been implicated in the induction of certain forms of synaptic plasticity, such as long-term potentiation and long-term depression, and these changes in synaptic efficacy have been proposed as neural mechanisms underlying memory and learning. The disruption of NMDA receptor targeting or signaling, as a result of the loss of SAP102, may lead to altered synaptic plasticity and may explain the intellectual impairment observed in individuals with DLG3 mutations.

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.

Mental retardation (MR) is the most frequent handicap among children and young adults. Although a large proportion of X-linked MR genes have been identified, only four genes responsible for autosomal-recessive nonsyndromic MR (AR-NSMR) have been described so far. Here, we report on two genes involved in autosomal-recessive and X-linked NSMR. First, autozygosity mapping in two sibs born to first-cousin French parents led to the identification of a region on 8p22-p23.1. This interval encompasses the gene N33/TUSC3 encoding one subunit of the oligosaccharyltransferase (OTase) complex, which catalyzes the transfer of an oligosaccharide chain on nascent proteins, the key step of N-glycosylation. Sequencing N33/TUSC3 identified a 1 bp insertion, c.787_788insC, resulting in a premature stop codon, p.N263fsX300, and leading to mRNA decay. Surprisingly, glycosylation analyses of patient fibroblasts showed normal N-glycan synthesis and transfer, suggesting that normal N-glycosylation observed in patient fibroblasts may be due to functional compensation. Subsequently, screening of the X-linked N33/TUSC3 paralog, the IAP gene, identified a missense mutation (c.932T-->G, p.V311G) in a family with X-linked NSMR. Recent studies of fucosylation and polysialic-acid modification of neuronal cell-adhesion glycoproteins have shown the critical role of glycosylation in synaptic plasticity. However, our data provide the first demonstration that a defect in N-glycosylation can result in NSMR. Together, our results demonstrate that fine regulation of OTase activity is essential for normal cognitive-function development, providing therefore further insights to understand the pathophysiological bases of MR.

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.

We have previously described a novel pathway controlling drug resistance, epithelial-to-mesenchymal transition (EMT) and stemness in breast cancer cells. Upstream in the pathway, three miRs (miR-106b, miR-93 and miR-25) target EP300, a transcriptional activator of E-cadherin. Upregulation of these miRs leads to the downregulation of EP300 and E-cadherin with initiation of an EMT. However, miRs regulate the expression of many genes, and the contribution to EMT by miR targets other than EP300 cannot be ruled out.

Linkage analysis and DNA sequencing in a family exhibiting an X-linked mental retardation (XLMR) syndrome, characterized by microcephaly, epilepsy, ataxia, and absent speech and resembling Angelman syndrome, identified a deletion in the SLC9A6 gene encoding the Na(+)/H(+) exchanger NHE6. Subsequently, other mutations were found in a male with mental retardation (MR) who had been investigated for Angelman syndrome and in two XLMR families with epilepsy and ataxia, including the family designated as having Christianson syndrome. Therefore, mutations in SLC9A6 cause X-linked mental retardation. Additionally, males with findings suggestive of unexplained Angelman syndrome should be considered as potential candidates for SLC9A6 mutations.

Cancer therapies directed at specific molecular targets in signaling pathways of cancer cells, such as tamoxifen, aromatase inhibitors and trastuzumab, have proven useful for treatment of advanced breast cancers. However, increased risk of endometrial cancer with long-term tamoxifen administration and of bone fracture due to osteoporosis in postmenopausal women undergoing aromatase inhibitor therapy are recognized side effects. These side effects as well as drug resistance make it necessary to search for novel molecular targets for drugs on the basis of well-characterized mechanisms of action.

The Nuclear Receptor (NR) superfamily of transcription factors comprises 48 members, several of which have been implicated in breast cancer. Most important is estrogen receptor-α (ERα), which is a key therapeutic target. ERα action is facilitated by co-operativity with other NR and there is evidence that ERα function may be recapitulated by other NRs in ERα-negative breast cancer. In order to examine the inter-relationships between nuclear receptors, and to obtain evidence for previously unsuspected roles for any NRs, we undertook quantitative RT-PCR and bioinformatics analysis to examine their expression in breast cancer. While most NRs were expressed, bioinformatic analyses differentiated tumours into distinct prognostic groups that were validated by analyzing public microarray data sets. Although ERα and progesterone receptor were dominant in distinguishing prognostic groups, other NR strengthened these groups. Clustering analysis identified several family members with potential importance in breast cancer. Specifically, RORγ is identified as being co-expressed with ERα, whilst several NRs are preferentially expressed in ERα-negative disease, with TLX expression being prognostic in this subtype. Functional studies demonstrated the importance of TLX in regulating growth and invasion in ERα-negative breast cancer cells.

The nature and pace of genome mutation is largely unknown. Because standard methods sequence DNA from populations of cells, the genetic composition of individual cells is lost, de novo mutations in cells are concealed within the bulk signal and per cell cycle mutation rates and mechanisms remain elusive. Although single-cell genome analyses could resolve these problems, such analyses are error-prone because of whole-genome amplification (WGA) artefacts and are limited in the types of DNA mutation that can be discerned. We developed methods for paired-end sequence analysis of single-cell WGA products that enable (i) detecting multiple classes of DNA mutation, (ii) distinguishing DNA copy number changes from allelic WGA-amplification artefacts by the discovery of matching aberrantly mapping read pairs among the surfeit of paired-end WGA and mapping artefacts and (iii) delineating the break points and architecture of structural variants. By applying the methods, we capture DNA copy number changes acquired over one cell cycle in breast cancer cells and in blastomeres derived from a human zygote after in vitro fertilization. Furthermore, we were able to discover and fine-map a heritable inter-chromosomal rearrangement t(1;16)(p36;p12) by sequencing a single blastomere. The methods will expedite applications in basic genome research and provide a stepping stone to novel approaches for clinical genetic diagnosis.

Melanoma is the most deadly form of skin cancer. Expression of oncogenic BRAF or NRAS, which are frequently mutated in human melanomas, promote the formation of nevi but are not sufficient for tumorigenesis. Even with germline mutated p53, these engineered melanomas present with variable onset and pathology, implicating additional somatic mutations in a multi-hit tumorigenic process.

The stromal microenvironment controls response to injury and inflammation, and is also an important determinant of cancer cell behavior. However, our understanding of its modulation by miRNA (miR) and their respective targets is still sparse. Here, we identified the miR-25-93-106b cluster and two new target genes as critical drivers for metastasis and immune evasion of cancer cells. Using miR-25-93-106b knockout mice or antagomiRs, we demonstrated regulation of the production of the chemoattractant CXCL12 controlling bone marrow metastasis. Moreover, we identified the immune checkpoint PD-L1 (CD274) as a novel miR-93/106b target playing a central role in diminishing tumor immunity. Eventually, upregulation of miR-93 and miR-106b via miR-mimics or treatment with an epigenetic reader domain (BET) inhibitor resulted in diminished expression of CXCL12 and PD-L1. These data suggest a potential new therapeutic rationale for use of BET inhibitors for dual targeting of cancers with strong immunosuppressive and metastatic phenotypes.

Clear cell renal cell carcinoma (ccRCC) is the most common form of adult kidney cancer, characterized by the presence of inactivating mutations in the VHL gene in most cases, and by infrequent somatic mutations in known cancer genes. To determine further the genetics of ccRCC, we have sequenced 101 cases through 3,544 protein-coding genes. Here we report the identification of inactivating mutations in two genes encoding enzymes involved in histone modification-SETD2, a histone H3 lysine 36 methyltransferase, and JARID1C (also known as KDM5C), a histone H3 lysine 4 demethylase-as well as mutations in the histone H3 lysine 27 demethylase, UTX (KMD6A), that we recently reported. The results highlight the role of mutations in components of the chromatin modification machinery in human cancer. Furthermore, NF2 mutations were found in non-VHL mutated ccRCC, and several other probable cancer genes were identified. These results indicate that substantial genetic heterogeneity exists in a cancer type dominated by mutations in a single gene, and that systematic screens will be key to fully determining the somatic genetic architecture of cancer.

The progression of chronic liver disease to hepatocellular carcinoma is caused by the acquisition of somatic mutations that affect 20-30 cancer genes1-8. Burdens of somatic mutations are higher and clonal expansions larger in chronic liver disease9-13 than in normal liver13-16, which enables positive selection to shape the genomic landscape9-13. Here we analysed somatic mutations from 1,590 genomes across 34 liver samples, including healthy controls, alcohol-related liver disease and non-alcoholic fatty liver disease. Seven of the 29 patients with liver disease had mutations in FOXO1, the major transcription factor in insulin signalling. These mutations affected a single hotspot within the gene, impairing the insulin-mediated nuclear export of FOXO1. Notably, six of the seven patients with FOXO1S22W hotspot mutations showed convergent evolution, with variants acquired independently by up to nine distinct hepatocyte clones per patient. CIDEB, which regulates lipid droplet metabolism in hepatocytes17-19, and GPAM, which produces storage triacylglycerol from free fatty acids20,21, also had a significant excess of mutations. We again observed frequent convergent evolution: up to fourteen independent clones per patient with CIDEB mutations and up to seven clones per patient with GPAM mutations. Mutations in metabolism genes were distributed across multiple anatomical segments of the liver, increased clone size and were seen in both alcohol-related liver disease and non-alcoholic fatty liver disease, but rarely in hepatocellular carcinoma. Master regulators of metabolic pathways are a frequent target of convergent somatic mutation in alcohol-related and non-alcoholic fatty liver disease.

Genome-wide DNA sequencing was used to decrypt the phylogeny of multiple samples from distinct areas of cancer and morphologically normal tissue taken from the prostates of three men. Mutations were present at high levels in morphologically normal tissue distant from the cancer, reflecting clonal expansions, and the underlying mutational processes at work in morphologically normal tissue were also at work in cancer. Our observations demonstrate the existence of ongoing abnormal mutational processes, consistent with field effects, underlying carcinogenesis. This mechanism gives rise to extensive branching evolution and cancer clone mixing, as exemplified by the coexistence of multiple cancer lineages harboring distinct ERG fusions within a single cancer nodule. Subsets of mutations were shared either by morphologically normal and malignant tissues or between different ERG lineages, indicating earlier or separate clonal cell expansions. Our observations inform on the origin of multifocal disease and have implications for prostate cancer therapy in individual cases.

Multiple somatic rearrangements are often found in cancer genomes; however, the underlying processes of rearrangement and their contribution to cancer development are poorly characterized. Here we use a paired-end sequencing strategy to identify somatic rearrangements in breast cancer genomes. There are more rearrangements in some breast cancers than previously appreciated. Rearrangements are more frequent over gene footprints and most are intrachromosomal. Multiple rearrangement architectures are present, but tandem duplications are particularly common in some cancers, perhaps reflecting a specific defect in DNA maintenance. Short overlapping sequences at most rearrangement junctions indicate that these have been mediated by non-homologous end-joining DNA repair, although varying sequence patterns indicate that multiple processes of this type are operative. Several expressed in-frame fusion genes were identified but none was recurrent. The study provides a new perspective on cancer genomes, highlighting the diversity of somatic rearrangements and their potential contribution to cancer development.

The application of paired-end next generation sequencing approaches has made it possible to systematically characterize rearrangements of the cancer genome to base-pair level. Utilizing this approach, we report the first detailed analysis of ovarian cancer rearrangements, comparing high-grade serous and clear cell cancers, and these histotypes with other solid cancers. Somatic rearrangements were systematically characterized in eight high-grade serous and five clear cell ovarian cancer genomes and we report here the identification of > 600 somatic rearrangements. Recurrent rearrangements of the transcriptional regulator gene, TSHZ3, were found in three of eight serous cases. Comparison to breast, pancreatic and prostate cancer genomes revealed that a subset of ovarian cancers share a marked tandem duplication phenotype with triple-negative breast cancers. The tandem duplication phenotype was not linked to BRCA1/2 mutation, suggesting that other common mechanisms or carcinogenic exposures are operative. High-grade serous cancers arising in women with germline BRCA1 or BRCA2 mutation showed a high frequency of small chromosomal deletions. These findings indicate that BRCA1/2 germline mutation may contribute to widespread structural change and that other undefined mechanism(s), which are potentially shared with triple-negative breast cancer, promote tandem chromosomal duplications that sculpt the ovarian cancer genome.