Rakesh Nagarajan

Recurring Mutations Found by Sequencing an Acute Myeloid Leukemia Genome

BACKGROUND The full complement of DNA mutations that are responsible for the pathogenesis of acute myeloid leukemia (AML) is not yet known. METHODS We used massively parallel DNA sequencing to obtain a very high level of coverage (approximately 98%) of a primary, cytogenetically normal, de novo genome for AML with minimal maturation (AML-M1) and a matched normal skin genome. RESULTS We identified 12 acquired (somatic) mutations within the coding sequences of genes and 52 somatic point mutations in conserved or regulatory portions of the genome. All mutations appeared to be heterozygous and present in nearly all cells in the tumor sample. Four of the 64 mutations occurred in at least 1 additional AML sample in 188 samples that were tested. Mutations in NRAS and NPM1 had been identified previously in patients with AML, but two other mutations had not been identified. One of these mutations, in the IDH1 gene, was present in 15 of 187 additional AML genomes tested and was strongly associated with normal cytogenetic status; it was present in 13 of 80 cytogenetically normal samples (16%). The other was a nongenic mutation in a genomic region with regulatory potential and conservation in higher mammals; we detected it in one additional AML tumor. The AML genome that we sequenced contains approximately 750 point mutations, of which only a small fraction are likely to be relevant to pathogenesis. CONCLUSIONS By comparing the sequences of tumor and skin genomes of a patient with AML-M1, we have identified recurring mutations that may be relevant for pathogenesis.

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.

EGR2 Mutations in Inherited Neuropathies Dominant-Negatively Inhibit Myelin Gene Expression

The identification of EGR2 mutations in patients with neuropathies and the phenotype Egr2/Krox20(-/-) have demonstrated that the Egr2 transcription factor is critical for peripheral nerve myelination. However, the mechanism by which these mutations cause disease remains unclear, as most patients present with disease in the heterozygous state, whereas Egr2(+/-) mice are phenotypically normal. To understand the effect of aberrant Egr2 activity on Schwann cell gene expression, we performed microarray expression profiling to identify genes regulated by Egr2 in Schwann cells. These include genes encoding myelin proteins and enzymes required for synthesis of normal myelin lipids. Using these newly identified targets, we have shown that neuropathy-associated EGR2 mutants dominant-negatively inhibit wild-type Egr2-mediated expression of essential myelin genes to levels sufficiently low to result in the abnormal myelination observed in these patients.

Acquired copy number alterations in adult acute myeloid leukemia genomes

Cytogenetic analysis of acute myeloid leukemia (AML) cells has accelerated the identification of genes important for AML pathogenesis. To complement cytogenetic studies and to identify genes altered in AML genomes, we performed genome-wide copy number analysis with paired normal and tumor DNA obtained from 86 adult patients with de novo AML using 1.85 million feature SNP arrays. Acquired copy number alterations (CNAs) were confirmed using an ultra-dense array comparative genomic hybridization platform. A total of 201 somatic CNAs were found in the 86 AML genomes (mean, 2.34 CNAs per genome), with French-American-British system M6 and M7 genomes containing the most changes (10–29 CNAs per genome). Twenty-four percent of AML patients with normal cytogenetics had CNA, whereas 40% of patients with an abnormal karyotype had additional CNA detected by SNP array, and several CNA regions were recurrent. The mRNA expression levels of 57 genes were significantly altered in 27 of 50 recurrent CNA regions <5 megabases in size. A total of 8 uniparental disomy (UPD) segments were identified in the 86 genomes; 6 of 8 UPD calls occurred in samples with a normal karyotype. Collectively, 34 of 86 AML genomes (40%) contained alterations not found with cytogenetics, and 98% of these regions contained genes. Of 86 genomes, 43 (50%) had no CNA or UPD at this level of resolution. In this study of 86 adult AML genomes, the use of an unbiased high-resolution genomic screen identified many genes not previously implicated in AML that may be relevant for pathogenesis, along with many known oncogenes and tumor suppressor genes.

Identification of a Novel TP53 Cancer Susceptibility Mutation Through Whole-Genome Sequencing of a Patient With Therapy-Related AML

CONTEXT The identification of patients with inherited cancer susceptibility syndromes facilitates early diagnosis, prevention, and treatment. However, in many cases of suspected cancer susceptibility, the family history is unclear and genetic testing of common cancer susceptibility genes is unrevealing. OBJECTIVE To apply whole-genome sequencing to a patient without any significant family history of cancer but with suspected increased cancer susceptibility because of multiple primary tumors to identify rare or novel germline variants in cancer susceptibility genes. DESIGN, SETTING, AND PARTICIPANT: Skin (normal) and bone marrow (leukemia) DNA were obtained from a patient with early-onset breast and ovarian cancer (negative for BRCA1 and BRCA2 mutations) and therapy-related acute myeloid leukemia (t-AML) and analyzed with the following: whole-genome sequencing using paired-end reads, single-nucleotide polymorphism (SNP) genotyping, RNA expression profiling, and spectral karyotyping. MAIN OUTCOME MEASURES Structural variants, copy number alterations, single-nucleotide variants, and small insertions and deletions (indels) were detected and validated using the described platforms. RESULTS; Whole-genome sequencing revealed a novel, heterozygous 3-kilobase deletion removing exons 7-9 of TP53 in the patients normal skin DNA, which was homozygous in the leukemia DNA as a result of uniparental disomy. In addition, a total of 28 validated somatic single-nucleotide variations or indels in coding genes, 8 somatic structural variants, and 12 somatic copy number alterations were detected in the patients leukemia genome. CONCLUSION Whole-genome sequencing can identify novel, cryptic variants in cancer susceptibility genes in addition to providing unbiased information on the spectrum of mutations in a cancer genome.

Distinct Genetic Signatures among Pilocytic Astrocytomas Relate to Their Brain Region Origin

Pilocytic astrocytomas (PAs) are the most common glioma in children. Whereas many PAs are slow-growing or clinically indolent, others exhibit more aggressive features with tumor recurrence and death. To identify genetic signatures that might predict PA clinical behavior, we did gene expression profiling on 41 primary PAs arising sporadically and in patients with neurofibromatosis type 1 (NF1). Whereas no expression signature was found that could discriminate clinically aggressive or recurrent tumors from more indolent cases, PAs arising in patients with NF1 did exhibit a unique gene expression pattern. In addition, we identified a gene expression signature that stratified PAs by location (supratentorial versus infratentorial). Lastly, we also identified a gene expression pattern common to PAs and normal mouse astrocytes and neural stem cells from these distinct brain regions as well as a gene expression pattern shared between PAs and another human glial tumor (ependymoma) arising supratentorially compared with those originating in the posterior fossa. These results suggest that glial tumors share an intrinsic, lineage-specific molecular signature that reflects the brain region in which their nonmalignant predecessors originated.

Validation of a Next-Generation Sequencing Assay for Clinical Molecular Oncology

Currently, oncology testing includes molecular studies and cytogenetic analysis to detect genetic aberrations of clinical significance. Next-generation sequencing (NGS) allows rapid analysis of multiple genes for clinically actionable somatic variants. The WUCaMP assay uses targeted capture for NGS analysis of 25 cancer-associated genes to detect mutations at actionable loci. We present clinical validation of the assay and a detailed framework for design and validation of similar clinical assays. Deep sequencing of 78 tumor specimens (≥ 1000× average unique coverage across the capture region) achieved high sensitivity for detecting somatic variants at low allele fraction (AF). Validation revealed sensitivities and specificities of 100% for detection of single-nucleotide variants (SNVs) within coding regions, compared with SNP array sequence data (95% CI = 83.4-100.0 for sensitivity and 94.2-100.0 for specificity) or whole-genome sequencing (95% CI = 89.1-100.0 for sensitivity and 99.9-100.0 for specificity) of HapMap samples. Sensitivity for detecting variants at an observed 10% AF was 100% (95% CI = 93.2-100.0) in HapMap mixes. Analysis of 15 masked specimens harboring clinically reported variants yielded concordant calls for 13/13 variants at AF of ≥ 15%. The WUCaMP assay is a robust and sensitive method to detect somatic variants of clinical significance in molecular oncology laboratories, with reduced time and cost of genetic analysis allowing for strategic patient management.

Deciphering peripheral nerve myelination by using Schwann cell expression profiling.

Although mutations in multiple genes are associated with inherited demyelinating neuropathies, the molecular components and pathways crucial for myelination remain largely unknown. To approach this question, we performed genome-wide expression analysis in several paradigms where the status of peripheral nerve myelination is dynamically changing. Anchor gene correlation analysis, a form of microarray analysis that integrates functional information, using correlation-based clustering, with a statistically rigorous test, the Westfall and Young step-down algorithm, was applied to this data set. Biological pathways active in myelination, genes encoding proteins involved in myelin synthesis, and genes whose mutation results in myelination defects were identified. Many known genes and previously uncharacterized ESTs not heretofore associated with myelination were also identified. One of these ESTs, MASR (myelin-associated SUR4 protein), encodes a member of the SUR4 family of fatty acid desaturases, enzymes involved in elongation of very long chain fatty acids. Its specific localization in myelinating Schwann cells indicates a crucial role for MASR in normal myelin lipid synthesis.

Nab proteins are essential for peripheral nervous system myelination

Mutations that disrupt Egr2 transcriptional activity cause severe demyelinating peripheral neuropathies. Here we provide evidence that Nab1 and Nab2 proteins are critical transcriptional modulators of Egr2 in myelinating Schwann cells. Like Egr2, these proteins are essential for Schwann cell differentiation into the myelinating state. Mice lacking both Nab1 and Nab2 show severe congenital hypomyelination of peripheral nerves, with Schwann cell development arresting at the promyelinating stage, despite elevated Egr2 expression. As observed for Egr2, Nab proteins are necessary for Schwann cells to exit the cell cycle, downregulate suppressed cAMP-inducible protein (SCIP) expression and upregulate expression of critical myelination genes. The mRNA expression signature of Schwann cells deficient in both Nab1 and Nab2 is highly similar to that of Egr2-deficient Schwann cells, further indicating that the Egr2/Nab protein complex is a key regulator of the Schwann cell myelination program and that disruption of this transcriptional complex is likely to result in Schwann cell dysfunction in patients with Egr2 mutations.

High throughput digital quantification of mRNA abundance in primary human acute myeloid leukemia samples

Acute promyelocytic leukemia (APL) is characterized by the t(15;17) chromosomal translocation, which results in fusion of the retinoic acid receptor alpha (RARA) gene to another gene, most commonly promyelocytic leukemia (PML). The resulting fusion protein, PML-RARA, initiates APL, which is a subtype (M3) of acute myeloid leukemia (AML). In this report, we identify a gene expression signature that is specific to M3 samples; it was not found in other AML subtypes and did not simply represent the normal gene expression pattern of primary promyelocytes. To validate this signature for a large number of genes, we tested a recently developed high throughput digital technology (NanoString nCounter). Nearly all of the genes tested demonstrated highly significant concordance with our microarray data (P < 0.05). The validated gene signature reliably identified M3 samples in 2 other AML datasets, and the validated genes were substantially enriched in our mouse model of APL, but not in a cell line that inducibly expressed PML-RARA. These results demonstrate that nCounter is a highly reproducible, customizable system for mRNA quantification using limited amounts of clinical material, which provides a valuable tool for biomarker measurement in low-abundance patient samples.