DNA sequencing of a cytogenetically normal acute myeloid leukaemia genome

Acute myeloid leukaemia is a highly malignant haematopoietic tumour that affects about 13,000 adults in the United States each year. The treatment of this disease has changed little in the past two decades, because most of the genetic events that initiate the disease remain undiscovered. Whole-genome sequencing is now possible at a reasonable cost and timeframe to use this approach for the unbiased discovery of tumour-specific somatic mutations that alter the protein-coding genes. Here we present the results obtained from sequencing a typical acute myeloid leukaemia genome, and its matched normal counterpart obtained from the same patient's skin. We discovered ten genes with acquired mutations; two were previously described mutations that are thought to contribute to tumour progression, and eight were new mutations present in virtually all tumour cells at presentation and relapse, the function of which is not yet known. Our study establishes whole-genome sequencing as an unbiased method for discovering cancer-initiating mutations in previously unidentified genes that may respond to targeted therapies.     

Mutations in kelch-like 3 and cullin 3 cause hypertension and electrolyte abnormalities

Hypertension affects one billion people and is a principal reversible risk factor for cardiovascular disease. Pseudohypoaldosteronism type II (PHAII), a rare Mendelian syndrome featuring hypertension, hyperkalaemia and metabolic acidosis, has revealed previously unrecognized physiology orchestrating the balance between renal salt reabsorption and K(+) and H(+) excretion. Here we used exome sequencing to identify mutations in kelch-like 3 (KLHL3) or cullin 3 (CUL3) in PHAII patients from 41 unrelated families. KLHL3 mutations are either recessive or dominant, whereas CUL3 mutations are dominant and predominantly de novo. CUL3 and BTB-domain-containing kelch proteins such as KLHL3 are components of cullin-RING E3 ligase complexes that ubiquitinate substrates bound to kelch propeller domains. Dominant KLHL3 mutations are clustered in short segments within the kelch propeller and BTB domains implicated in substrate and cullin binding, respectively. Diverse CUL3 mutations all result in skipping of exon 9, producing an in-frame deletion. Because dominant KLHL3 and CUL3 mutations both phenocopy recessive loss-of-function KLHL3 mutations, they may abrogate ubiquitination of KLHL3 substrates. Disease features are reversed by thiazide diuretics, which inhibit the Na-Cl cotransporter in the distal nephron of the kidney; KLHL3 and CUL3 are expressed in this location, suggesting a mechanistic link between KLHL3 and CUL3 mutations, increased Na-Cl reabsorption, and disease pathogenesis. These findings demonstrate the utility of exome sequencing in disease gene identification despite the combined complexities of locus heterogeneity, mixed models of transmission and frequent de novo mutation, and establish a fundamental role for KLHL3 and CUL3 in blood pressure, K(+) and pH homeostasis.     

Somatic mutations affect key pathways in lung adenocarcinoma

Determining the genetic basis of cancer requires comprehensive analyses of large collections of histopathologically well-classified primary tumours. Here we report the results of a collaborative study to discover somatic mutations in 188 human lung adenocarcinomas. DNA sequencing of 623 genes with known or potential relationships to cancer revealed more than 1,000 somatic mutations across the samples. Our analysis identified 26 genes that are mutated at significantly high frequencies and thus are probably involved in carcinogenesis. The frequently mutated genes include tyrosine kinases, among them the EGFR homologue ERBB4; multiple ephrin receptor genes, notably EPHA3; vascular endothelial growth factor receptor KDR; and NTRK genes. These data provide evidence of somatic mutations in primary lung adenocarcinoma for several tumour suppressor genes involved in other cancers--including NF1, APC, RB1 and ATM--and for sequence changes in PTPRD as well as the frequently deleted gene LRP1B. The observed mutational profiles correlate with clinical features, smoking status and DNA repair defects. These results are reinforced by data integration including single nucleotide polymorphism array and gene expression array. Our findings shed further light on several important signalling pathways involved in lung adenocarcinoma, and suggest new molecular targets for treatment.     

Whole-genome analysis informs breast cancer response to aromatase inhibition

To correlate the variable clinical features of oestrogen-receptor-positive breast cancer with somatic alterations, we studied pretreatment tumour biopsies accrued from patients in two studies of neoadjuvant aromatase inhibitor therapy by massively parallel sequencing and analysis. Eighteen significantly mutated genes were identified, including five genes (RUNX1, CBFB, MYH9, MLL3 and SF3B1) previously linked to haematopoietic disorders. Mutant MAP3K1 was associated with luminal A status, low-grade histology and low proliferation rates, whereas mutant TP53 was associated with the opposite pattern. Moreover, mutant GATA3 correlated with suppression of proliferation upon aromatase inhibitor treatment. Pathway analysis demonstrated that mutations in MAP2K4, a MAP3K1 substrate, produced similar perturbations as MAP3K1 loss. Distinct phenotypes in oestrogen-receptor-positive breast cancer are associated with specific patterns of somatic mutations that map into cellular pathways linked to tumour biology, but most recurrent mutations are relatively infrequent. Prospective clinical trials based on these findings will require comprehensive genome sequencing.     

Oncogenically active MYD88 mutations in human lymphoma

The activated B-cell-like (ABC) subtype of diffuse large B-cell lymphoma (DLBCL) remains the least curable form of this malignancy despite recent advances in therapy. Constitutive nuclear factor (NF)-κB and JAK kinase signalling promotes malignant cell survival in these lymphomas, but the genetic basis for this signalling is incompletely understood. Here we describe the dependence of ABC DLBCLs on MYD88, an adaptor protein that mediates toll and interleukin (IL)-1 receptor signalling, and the discovery of highly recurrent oncogenic mutations affecting MYD88 in ABC DLBCL tumours. RNA interference screening revealed that MYD88 and the associated kinases IRAK1 and IRAK4 are essential for ABC DLBCL survival. High-throughput RNA resequencing uncovered MYD88 mutations in ABC DLBCL lines. Notably, 29% of ABC DLBCL tumours harboured the same amino acid substitution, L265P, in the MYD88 Toll/IL-1 receptor (TIR) domain at an evolutionarily invariant residue in its hydrophobic core. This mutation was rare or absent in other DLBCL subtypes and Burkitt's lymphoma, but was observed in 9% of mucosa-associated lymphoid tissue lymphomas. At a lower frequency, additional mutations were observed in the MYD88 TIR domain, occurring in both the ABC and germinal centre B-cell-like (GCB) DLBCL subtypes. Survival of ABC DLBCL cells bearing the L265P mutation was sustained by the mutant but not the wild-type MYD88 isoform, demonstrating that L265P is a gain-of-function driver mutation. The L265P mutant promoted cell survival by spontaneously assembling a protein complex containing IRAK1 and IRAK4, leading to IRAK4 kinase activity, IRAK1 phosphorylation, NF-κB signalling, JAK kinase activation of STAT3, and secretion of IL-6, IL-10 and interferon-β. Hence, the MYD88 signalling pathway is integral to the pathogenesis of ABC DLBCL, supporting the development of inhibitors of IRAK4 kinase and other components of this pathway for the treatment of tumours bearing oncogenic MYD88 mutations.     

Clonal evolution in relapsed acute myeloid leukaemia revealed by whole-genome sequencing

Most patients with acute myeloid leukaemia (AML) die from progressive disease after relapse, which is associated with clonal evolution at the cytogenetic level. To determine the mutational spectrum associated with relapse, we sequenced the primary tumour and relapse genomes from eight AML patients, and validated hundreds of somatic mutations using deep sequencing; this allowed us to define clonality and clonal evolution patterns precisely at relapse. In addition to discovering novel, recurrently mutated genes (for example, WAC, SMC3, DIS3, DDX41 and DAXX) in AML, we also found two major clonal evolution patterns during AML relapse: (1) the founding clone in the primary tumour gained mutations and evolved into the relapse clone, or (2) a subclone of the founding clone survived initial therapy, gained additional mutations and expanded at relapse. In all cases, chemotherapy failed to eradicate the founding clone. The comparison of relapse-specific versus primary tumour mutations in all eight cases revealed an increase in transversions, probably due to DNA damage caused by cytotoxic chemotherapy. These data demonstrate that AML relapse is associated with the addition of new mutations and clonal evolution, which is shaped, in part, by the chemotherapy that the patients receive to establish and maintain remissions.     

The landscape of cancer genes and mutational processes in breast cancer

All cancers carry somatic mutations in their genomes. A subset, known as driver mutations, confer clonal selective advantage on cancer cells and are causally implicated in oncogenesis, and the remainder are passenger mutations. The driver mutations and mutational processes operative in breast cancer have not yet been comprehensively explored. Here we examine the genomes of 100 tumours for somatic copy number changes and mutations in the coding exons of protein-coding genes. The number of somatic mutations varied markedly between individual tumours. We found strong correlations between mutation number, age at which cancer was diagnosed and cancer histological grade, and observed multiple mutational signatures, including one present in about ten per cent of tumours characterized by numerous mutations of cytosine at TpC dinucleotides. Driver mutations were identified in several new cancer genes including AKT2, ARID1B, CASP8, CDKN1B, MAP3K1, MAP3K13, NCOR1, SMARCD1 and TBX3. Among the 100 tumours, we found driver mutations in at least 40 cancer genes and 73 different combinations of mutated cancer genes. The results highlight the substantial genetic diversity underlying this common disease.     

Novel mutations target distinct subgroups of medulloblastoma

Medulloblastoma is a malignant childhood brain tumour comprising four discrete subgroups. Here, to identify mutations that drive medulloblastoma, we sequenced the entire genomes of 37 tumours and matched normal blood. One-hundred and thirty-six genes harbouring somatic mutations in this discovery set were sequenced in an additional 56 medulloblastomas. Recurrent mutations were detected in 41 genes not yet implicated in medulloblastoma; several target distinct components of the epigenetic machinery in different disease subgroups, such as regulators of H3K27 and H3K4 trimethylation in subgroups 3 and 4 (for example, KDM6A and ZMYM3), and CTNNB1-associated chromatin re-modellers in WNT-subgroup tumours (for example, SMARCA4 and CREBBP). Modelling of mutations in mouse lower rhombic lip progenitors that generate WNT-subgroup tumours identified genes that maintain this cell lineage (DDX3X), as well as mutated genes that initiate (CDH1) or cooperate (PIK3CA) in tumorigenesis. These data provide important new insights into the pathogenesis of medulloblastoma subgroups and highlight targets for therapeutic development.     

Genome-wide analysis of genetic alterations in acute lymphoblastic leukaemia

Chromosomal aberrations are a hallmark of acute lymphoblastic leukaemia (ALL) but alone fail to induce leukaemia. To identify cooperating oncogenic lesions, we performed a genome-wide analysis of leukaemic cells from 242 paediatric ALL patients using high-resolution, single-nucleotide polymorphism arrays and genomic DNA sequencing. Our analyses revealed deletion, amplification, point mutation and structural rearrangement in genes encoding principal regulators of B lymphocyte development and differentiation in 40% of B-progenitor ALL cases. The PAX5 gene was the most frequent target of somatic mutation, being altered in 31.7% of cases. The identified PAX5 mutations resulted in reduced levels of PAX5 protein or the generation of hypomorphic alleles. Deletions were also detected in TCF3 (also known as E2A), EBF1, LEF1, IKZF1 (IKAROS) and IKZF3 (AIOLOS). These findings suggest that direct disruption of pathways controlling B-cell development and differentiation contributes to B-progenitor ALL pathogenesis. Moreover, these data demonstrate the power of high-resolution, genome-wide approaches to identify new molecular lesions in cancer.     

Two levels of protection for the B cell genome during somatic hypermutation

Somatic hypermutation introduces point mutations into immunoglobulin genes in germinal centre B cells during an immune response. The reaction is initiated by cytosine deamination by the activation-induced deaminase (AID) and completed by error-prone processing of the resulting uracils by mismatch and base excision repair factors. Somatic hypermutation represents a threat to genome integrity and it is not known how the B cell genome is protected from the mutagenic effects of somatic hypermutation nor how often these protective mechanisms fail. Here we show, by extensive sequencing of murine B cell genes, that the genome is protected by two distinct mechanisms: selective targeting of AID and gene-specific, high-fidelity repair of AID-generated uracils. Numerous genes linked to B cell tumorigenesis, including Myc, Pim1, Pax5, Ocab (also called Pou2af1), H2afx, Rhoh and Ebf1, are deaminated by AID but escape acquisition of most mutations through the combined action of mismatch and base excision repair. However, approximately 25% of expressed genes analysed were not fully protected by either mechanism and accumulated mutations in germinal centre B cells. Our results demonstrate that AID acts broadly on the genome, with the ultimate distribution of mutations determined by a balance between high-fidelity and error-prone DNA repair.     

Lack of heterogeneity in antihapten antibodies of a phylogenetically primitive shark

Individual mammals have the capacity to express at least one million distinct antigen binding specificities, implying a high degree of structural heterogeneity in the variable heavy and light chain (VH and VL) portions of the antibody molecules. Studies of higher vertebrate species suggest that this heterogeneity is created both through a sizeable repertoire of germ-line VH and VL genes and through random rearrangements of V and joining genes. Additional somatic mechanisms probably also contribute to the ultimate heterogeneity; one-third of murine plasmacytomas producing lambda 1 immunoglobulin carry a somatically mutated Ig1-V gene. The relative contributions of these various mechanisms to the overall imunoglobulin variability are difficult to evaluate. The production of different antibodies to a defined determinant in different individuals of an inbred mouse strain [for example, (3-iodo-4-hydroxy-5-nitrophenyl) acetyl (NIP) in CBA mice] suggests the involvement of somatic mutations or rearrangement but does not rule out the possibility that each individual CBA mouse expresses only a small random fraction from a large germ-line repertoire of V genes determining different anti-NIP binding sites. The opposite finding, that different individuals produce nearly identical antibodies to a defined determinant, would suggest the presence and expression of a limited number of germ-line genes without somatic alterations. Data presented here suggest that primitive sharks (Heterodontus fransisci) produce such antibodies to the hapten furyloxazolone.     

Mutations of the BRAF gene in human cancer

Cancers arise owing to the accumulation of mutations in critical genes that alter normal programmes of cell proliferation, differentiation and death. As the first stage of a systematic genome-wide screen for these genes, we have prioritized for analysis signalling pathways in which at least one gene is mutated in human cancer. The RAS RAF MEK ERK MAP kinase pathway mediates cellular responses to growth signals. RAS is mutated to an oncogenic form in about 15% of human cancer. The three RAF genes code for cytoplasmic serine/threonine kinases that are regulated by binding RAS. Here we report BRAF somatic missense mutations in 66% of malignant melanomas and at lower frequency in a wide range of human cancers. All mutations are within the kinase domain, with a single substitution (V599E) accounting for 80%. Mutated BRAF proteins have elevated kinase activity and are transforming in NIH3T3 cells. Furthermore, RAS function is not required for the growth of cancer cell lines with the V599E mutation. As BRAF is a serine/threonine kinase that is commonly activated by somatic point mutation in human cancer, it may provide new therapeutic opportunities in malignant melanoma.     

Patterns of somatic mutation in human cancer genomes

Cancers arise owing to mutations in a subset of genes that confer growth advantage. The availability of the human genome sequence led us to propose that systematic resequencing of cancer genomes for mutations would lead to the discovery of many additional cancer genes. Here we report more than 1,000 somatic mutations found in 274 megabases (Mb) of DNA corresponding to the coding exons of 518 protein kinase genes in 210 diverse human cancers. There was substantial variation in the number and pattern of mutations in individual cancers reflecting different exposures, DNA repair defects and cellular origins. Most somatic mutations are likely to be 'passengers' that do not contribute to oncogenesis. However, there was evidence for 'driver' mutations contributing to the development of the cancers studied in approximately 120 genes. Systematic sequencing of cancer genomes therefore reveals the evolutionary diversity of cancers and implicates a larger repertoire of cancer genes than previously anticipated.     

The role of clonal selection and somatic mutation in autoimmunity

Polyclonal activation has been proposed as the reason that autoantibodies are produced during autoimmune disease. This model denies a role for specific antigen selection of B cells and predicts instead a multiclonal population of unmutated or randomly mutated autoantibodies. We have found that the genetic features and clonal composition of spontaneously derived immunoglobulin G (IgG) antiself-IgG (rheumatoid factor (RF] autoantibodies derived from the autoimmune MRL/lpr mouse strain are inconsistent with both the predictions of this model and the actual outcome of experimental polyclonal activation. Instead we have found that MRL/lpr RFs are oligoclonal or even monoclonal in origin. They harbour numerous somatic mutations which are distributed in a way that suggests immunoglobulin-receptor-dependent selection of these mutations. In this sense, the MRL/lpr RFs resemble antibodies elicited by exogenous antigens after secondary immunization. The parallels suggest that, like secondary immune responses, antigen stimulation is important in the generation of MRL/lpr RFs.