Search

We performed an integrated genomic, transcriptomic and proteomic characterization of 373 endometrial carcinomas using array- and sequencing-based technologies. Uterine serous tumours and ~25% of high-grade endometrioid tumours had extensive copy number alterations, few DNA methylation changes, low oestrogen receptor/progesterone receptor levels, and frequent TP53 mutations. Most endometrioid tumours had few copy number alterations or TP53 mutations, but frequent mutations in PTEN, CTNNB1, PIK3CA, ARID1A and KRAS and novel mutations in the SWI/SNF chromatin remodelling complex gene ARID5B. A subset of endometrioid tumours that we identified had a markedly increased transversion mutation frequency and newly identified hotspot mutations in POLE. Our results classified endometrial cancers into four categories: POLE ultramutated, microsatellite instability hypermutated, copy-number low, and copy-number high. Uterine serous carcinomas share genomic features with ovarian serous and basal-like breast carcinomas. We demonstrated that the genomic features of endometrial carcinomas permit a reclassification that may affect post-surgical adjuvant treatment for women with aggressive tumours., National Institutes of Health (U.S.) (Grant 5U24CA143799-04), National Institutes of Health (U.S.) (Grant 5U24CA143835-04), National Institutes of Health (U.S.) (Grant 5U24CA143840-04), National Institutes of Health (U.S.) (Grant 5U24CA143843-04), National Institutes of Health (U.S.) (Grant 5U24CA143845-04), National Institutes of Health (U.S.) (Grant 5U24CA143848-04), National Institutes of Health (U.S.) (Grant 5U24CA143858-04), National Institutes of Health (U.S.) (Grant 5U24CA143866-04), National Institutes of Health (U.S.) (Grant 5U24CA143867-04), National Institutes of Health (U.S.) (Grant 5U24CA143882-04), National Institutes of Health (U.S.) (Grant 5U24CA143883-04), National Institutes of Health (U.S.) (Grant 5U24CA144025-04), National Institutes of Health (U.S.) (Grant U54HG003067-11), National Institutes of Health (U.S.) (Grant U54HG003079-10), National Institutes of Health (U.S.) (Grant U54HG003273-10)

The most common pediatric brain tumors are low-grade gliomas (LGGs). We used whole-genome sequencing to identify multiple new genetic alterations involving BRAF, RAF1, FGFR1, MYB, MYBL1 and genes with histone-related functions, including H3F3A and ATRX, in 39 LGGs and low-grade glioneuronal tumors (LGGNTs). Only a single non-silent somatic alteration was detected in 24 of 39 (62%) tumors. Intragenic duplications of the portion of FGFR1 encoding the tyrosine kinase domain (TKD) and rearrangements of MYB were recurrent and mutually exclusive in 53% of grade II diffuse LGGs. Transplantation of Trp53-null neonatal astrocytes expressing FGFR1 with the duplication involving the TKD into the brains of nude mice generated high-grade astrocytomas with short latency and 100% penetrance. FGFR1 with the duplication induced FGFR1 autophosphorylation and upregulation of the MAPK/ERK and PI3K pathways, which could be blocked by specific inhibitors. Focusing on the therapeutically challenging diffuse LGGs, our study of 151 tumors has discovered genetic alterations and potential therapeutic targets across the entire range of pediatric LGGs and LGGNTs.

The genetic basis of hypodiploid acute lymphoblastic leukemia (ALL), a subtype of ALL characterized by aneuploidy and poor outcome, is unknown. Genomic profiling of 124 hypodiploid ALL cases, including whole-genome and exome sequencing of 40 cases, identified two subtypes that differ in the severity of aneuploidy, transcriptional profiles and submicroscopic genetic alterations. Near-haploid ALL with 24-31 chromosomes harbor alterations targeting receptor tyrosine kinase signaling and Ras signaling (71%) and the lymphoid transcription factor gene IKZF3 (encoding AIOLOS; 13%). In contrast, low-hypodiploid ALL with 32-39 chromosomes are characterized by alterations in TP53 (91.2%) that are commonly present in nontumor cells, IKZF2 (encoding HELIOS; 53%) and RB1 (41%). Both near-haploid and low-hypodiploid leukemic cells show activation of Ras-signaling and phosphoinositide 3-kinase (PI3K)-signaling pathways and are sensitive to PI3K inhibitors, indicating that these drugs should be explored as a new therapeutic strategy for this aggressive form of leukemia.

SummaryWe report the results of whole-genome and transcriptome sequencing of tumor and adjacent normal tissue samples from 17 patients with non-small cell lung carcinoma (NSCLC). We identified 3,726 point mutations and more than 90 indels in the coding sequence, with an average mutation frequency more than 10-fold higher in smokers than in never-smokers. Novel alterations in genes involved in chromatin modification and DNA repair pathways were identified, along with DACH1, CFTR, RELN, ABCB5, and HGF. Deep digital sequencing revealed diverse clonality patterns in both never-smokers and smokers. All validated EFGR and KRAS mutations were present in the founder clones, suggesting possible roles in cancer initiation. Analysis revealed 14 fusions, including ROS1 and ALK, as well as novel metabolic enzymes. Cell-cycle and JAK-STAT pathways are significantly altered in lung cancer, along with perturbations in 54 genes that are potentially targetable with currently available drugs.

Massively parallel sequencing technology and the associated rapidly decreasing sequencing costs have enabled systemic analyses of somatic mutations in large cohorts of cancer cases. Here we introduce a comprehensive mutational analysis pipeline that uses standardized sequence-based inputs along with multiple types of clinical data to establish correlations among mutation sites, affected genes and pathways, and to ultimately separate the commonly abundant passenger mutations from the truly significant events. In other words, we aim to determine the Mutational Significance in Cancer (MuSiC) for these large data sets. The integration of analytical operations in the MuSiC framework is widely applicable to a broad set of tumor types and offers the benefits of automation as well as standardization. Herein, we describe the computational structure and statistical underpinnings of the MuSiC pipeline and demonstrate its performance using 316 ovarian cancer samples from the TCGA ovarian cancer project. MuSiC correctly confirms many expected results, and identifies several potentially novel avenues for discovery.

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.

A variety of microbial communities and their genes (the microbiome) exist throughout the human body, with fundamental roles in human health and disease. The National Institutes of Health (NIH)-funded Human Microbiome Project Consortium has established a population-scale framework to develop metagenomic protocols, resulting in a broad range of quality-controlled resources and data including standardized methods for creating, processing and interpreting distinct types of high-throughput metagenomic data available to the scientific community. Here we present resources from a population of 242 healthy adults sampled at 15 or 18 body sites up to three times, which have generated 5,177 microbial taxonomic profiles from 16S ribosomal RNA genes and over 3.5 terabases of metagenomic sequence so far. In parallel, approximately 800 reference strains isolated from the human body have been sequenced. Collectively, these data represent the largest resource describing the abundance and variety of the human microbiome, while providing a framework for current and future studies.

Summary To correlate the variable clinical features of estrogen receptor positive (ER+) breast cancer with somatic alterations, we studied pre-treatment tumour biopsies accrued from patients in a study of neoadjuvant aromatase inhibitor (AI) 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 hematopoietic disorders. Mutant MAP3K1 was associated with Luminal A status, low grade histology and low proliferation rates whereas mutant TP53 associated with the opposite pattern. Moreover, mutant GATA3 correlated with suppression of proliferation upon AI treatment. Pathway analysis demonstrated mutations in MAP2K4, a MAP3K1 substrate, produced similar perturbations as MAP3K1 loss. Distinct phenotypes in ER+ breast cancer are associated with specific patterns of somatic mutations that map into cellular pathways linked to tumor biology but most recurrent mutations are relatively infrequent. Prospective clinical trials based on these findings will require comprehensive genome sequencing.

Retinoblastoma is an aggressive childhood cancer of the developing retina that is initiated by the biallelic loss of RB1. Tumours progress very quickly following RB1 inactivation but the underlying mechanism is not known. Here we show that the retinoblastoma genome is stable, but that multiple cancer pathways can be epigenetically deregulated. To identify the mutations that cooperate with RB1 loss, we performed whole-genome sequencing of retinoblastomas. The overall mutational rate was very low; RB1 was the only known cancer gene mutated. We then evaluated the role of RB1 in genome stability and considered non-genetic mechanisms of cancer pathway deregulation. For example, the proto-oncogene SYK is upregulated in retinoblastoma and is required for tumour cell survival. Targeting SYK with a small-molecule inhibitor induced retinoblastoma tumour cell death in vitro and in vivo. Thus, retinoblastomas may develop quickly as a result of the epigenetic deregulation of key cancer pathways as a direct or indirect result of RB1 loss.

Early T-cell precursor acute lymphoblastic leukaemia (ETP ALL) is an aggressive malignancy of unknown genetic basis. We performed whole-genome sequencing of 12 ETP ALL cases and assessed the frequency of the identified somatic mutations in 94 T-cell acute lymphoblastic leukaemia cases. ETP ALL was characterized by activating mutations in genes regulating cytokine receptor and RAS signalling (67% of cases; NRAS, KRAS, FLT3, IL7R, JAK3, JAK1, SH2B3 and BRAF), inactivating lesions disrupting haematopoietic development (58%; GATA3, ETV6, RUNX1, IKZF1 and EP300) and histone-modifying genes (48%; EZH2, EED, SUZ12, SETD2 and EP300). We also identified new targets of recurrent mutation including DNM2, ECT2L and RELN. The mutational spectrum is similar to myeloid tumours, and moreover, the global transcriptional profile of ETP ALL was similar to that of normal and myeloid leukaemia haematopoietic stem cells. These findings suggest that addition of myeloid-directed therapies might improve the poor outcome of ETP ALL.

BACKGROUND The genetic alterations responsible for an adverse outcome in most patients with acute myeloid leukemia (AML) are unknown. METHODS Using massively parallel DNA sequencing, we identified a somatic mutation in DNMT3A, encoding a DNA methyltransferase, in the genome of cells from a patient with AML with a normal karyotype. We sequenced the exons of DNMT3A in 280 additional patients with de novo AML to define recurring mutations. RESULTS A total of 62 of 281 patients (22.1%) had mutations in DNMT3A that were predicted to affect translation. We identified 18 different missense mutations, the most common of which was predicted to affect amino acid R882 (in 37 patients). We also identified six frameshift, six nonsense, and three splice-site mutations and a 1.5-Mbp deletion encompassing DNMT3A. These mutations were highly enriched in the group of patients with an intermediate-risk cytogenetic profile (56 of 166 patients, or 33.7%) but were absent in all 79 patients with a favorable-risk cytogenetic profile (P

Massively parallel DNA sequencing technologies provide an unprecedented ability to screen entire genomes for genetic changes associated with tumor progression. Here we describe the genomic analyses of four DNA samples from an African-American patient with basal-like breast cancer: peripheral blood, the primary tumor, a brain metastasis, and a xenograft derived from the primary tumor. The metastasis contained two de novo mutations and a large deletion not present in the primary tumor, and was significantly enriched for 20 shared mutations. The xenograft retained all primary tumor mutations, and displayed a mutation enrichment pattern that paralleled the metastasis (16 of 20 genes). Two overlapping large deletions, encompassing CTNNA1, were present in all three tumor samples. The differential mutation frequencies and structural variation patterns in metastasis and xenograft compared to the primary tumor suggest that secondary tumors may arise from a minority of cells within the primary.

Lay Summary Acute myeloid leukemia is a highly malignant hematopoietic tumor that affects about 13,000 adults yearly in the United States. The treatment of this disease has changed little in the past two decades, since most of the genetic events that initiate the disease remain undiscovered. Whole genome sequencing is now possible at a reasonable cost and timeframe to utilize this approach for unbiased discovery of tumor-specific somatic mutations that alter the protein-coding genes. Here we show the results obtained by sequencing a typical acute myeloid leukemia genome and its matched normal counterpart, obtained from the patient’s skin. We discovered 10 genes with acquired mutations; two were previously described mutations thought to contribute to tumor progression, and 8 were novel mutations present in virtually all tumor cells at presentation and relapse, whose function is not yet known. Our study establishes whole genome sequencing as an unbiased method for discovering initiating mutations in cancer genomes, and for identifying novel genes that may respond to targeted therapies. We used massively parallel sequencing technology to sequence the genomic DNA of tumor and normal skin cells obtained from a patient with a typical presentation of FAB M1 Acute Myeloid Leukemia (AML) with normal cytogenetics. 32.7-fold ‘haploid’ coverage (98 billion bases) was obtained for the tumor genome, and 13.9-fold coverage (41.8 billion bases) was obtained for the normal sample. Of 2,647,695 well-supported Single Nucleotide Variants (SNVs) found in the tumor genome, 2,588,486 (97.7%) also were detected in the patient’s skin genome, limiting the number of variants that required further study. For the purposes of this initial study, we restricted our downstream analysis to the coding sequences of annotated genes: we found only eight heterozygous, non-synonymous somatic SNVs in the entire genome. All were novel, including mutations in protocadherin/cadherin family members (CDH24 and PCLKC), G-protein coupled receptors (GPR123 and EBI2), a protein phosphatase (PTPRT), a potential guanine nucleotide exchange factor (KNDC1), a peptide/drug transporter (SLC15A1), and a glutamate receptor gene (GRINL1B). We also detected previously described, recurrent somatic insertions in the FLT3 and NPM1 genes. Based on deep readcount data, we determined that all of these mutations (except FLT3) were present in nearly all tumor cells at presentation, and again at relapse 11 months later, suggesting that the patient had a single dominant clone containing all of the mutations. These results demonstrate the power of whole genome sequencing to discover novel cancer-associated mutations.

Human cancer cells typically harbour multiple chromosomal aberrations, nucleotide substitutions and epigenetic modifications that drive malignant transformation. The Cancer Genome Atlas ( TCGA) pilot project aims to assess the value of large- scale multi- dimensional analysis of these molecular characteristics in human cancer and to provide the data rapidly to the research community. Here we report the interim integrative analysis of DNA copy number, gene expression and DNA methylation aberrations in 206 glioblastomas - the most common type of primary adult brain cancer - and nucleotide sequence aberrations in 91 of the 206 glioblastomas. This analysis provides new insights into the roles of ERBB2, NF1 and TP53, uncovers frequent mutations of the phosphatidylinositol- 3- OH kinase regulatory subunit gene PIK3R1, and provides a network view of the pathways altered in the development of glioblastoma. Furthermore, integration of mutation, DNA methylation and clinical treatment data reveals a link between MGMT promoter methylation and a hypermutator phenotype consequent to mismatch repair deficiency in treated glioblastomas, an observation with potential clinical implications. Together, these findings establish the feasibility and power of TCGA, demonstrating that it can rapidly expand knowledge of the molecular basis of cancer.

In this work, we present the Genome Modeling System (GMS), an analysis information management system capable of executing automated genome analysis pipelines at a massive scale. The GMS framework provides detailed tracking of samples and data coupled with reliable and repeatable analysis pipelines. The GMS also serves as a platform for bioinformatics development, allowing a large team to collaborate on data analysis, or an individual researcher to leverage the work of others effectively within its data management system. Rather than separating ad-hoc analysis from rigorous, reproducible pipelines, the GMS promotes systematic integration between the two. As a demonstration of the GMS, we performed an integrated analysis of whole genome, exome and transcriptome sequencing data from a breast cancer cell line (HCC1395) and matched lymphoblastoid line (HCC1395BL). These data are available for users to test the software, complete tutorials and develop novel GMS pipeline configurations. The GMS is available at https://github.com/genome/gms.

The combined insight gained from laboratory and computational experiments has historically provided detailed information about the complex interactions present in many heterogeneous catalytic systems. Although this information has been effectively used to interpret and explain experimental observations, it unfortunately has not, in general, been incorporated into kinetic models of the reacting systems. In particular, catalytic nonuniformity is one such phenomenon which is often ignored when kinetic models are constructed despite there being much experimental and computational evidence indicating its presence. Catalytic nonuniformity has been therefore incorporated into both microkinetic models and dynamic Monte Carlo simulations of catalytic reacting systems, and the performance and applicability of each solution method has been investigated for several model systems exhibiting different causes of nonuniform behavior.

In recent years, the use of mechanistic modeling to identify the underlying kinetics of complex systems has increased greatly. One of the challenges to kinetic modeling is the construction of a model which can capture the essential chemistry of a system while a manageable size is retained. The rate-based generation of mechanistic models is an attractive approach because kinetically significant species can be determined and selectively included in the final mechanism. An algorithm for the rate-based generation of reaction mechanisms developed previously1 was improved and used to construct a compact mechanistic model for low-pressure tetradecane pyrolysis. Though thousands of species and reactions were generated, only a small portion of these (2% of species and 20% of reactions) was deemed necessary and incorporated into the final model. Experimental data were used to determine frequency factors for a subset of the reaction families, while all other kinetic parameters were set on the basis of literature val...

Recently, the synthesis of ammonia over ruthenium-based catalysts has become an industrially viable process. Unfortunately, investigations of ammonia synthesis over ruthenium are scarce, particularly in comparison to the number of studies carried out over iron. To begin to fill this void, we have performed a series of electronic density-functional theory (DFT) calculations to investigate the effect of particle size and surface structure on ammonia synthesis over ruthenium. Our study has focused on the dissociative adsorption of dinitrogen, which is thought to be the rate-determining step in the synthesis, on both single-crystal surfaces and spherical clusters of ruthenium. The equilibrium adsorbate geometries were remarkably similar on both the single-crystal surfaces and the spherical clusters studied. The binding energy of dinitrogen in the end-on state exhibited a strong dependence on ruthenium surface atom coordination, being much stronger on atoms with low coordination. The main difference between the two single-crystal surfaces studied was the ability of the open Ru(1120) face to stabilize a low-energy side-on dinitrogen state, while the close-packed Ru(0001) face could not. It is likely that this stable side-on state provides a low-energy dissociation pathway.

Nonuniformity, the variation in heats of adsorption and kinetic parameters with changes in coverage or adsorption site, is present in many heterogeneous catalytic reaction systems. Unfortunately, there have been few investigations addressing the incorporation of nonuniformity into kinetic models of catalytic systems. The main reason for this disparity is the difficulty encountered in establishing a quantitative assessment of the effect of nonuniformity on reaction kinetics. Indeed, many systems in which nonuniformity is known to be present are well described by uniform kinetic models. To begin to get a better understanding of the effect nonuniformity has on reaction kinetics, both model and real systems in which nonuniformity is known to be present have been studied using microkinetic modeling. The model system was a simple condensation mechanism involving only one surface reaction. Two types of nonuniformity, adsorbate-adsorbate interactions and biographic heterogeneity, were incorporated into the kinetic parameters of the model, and two sets of nonuniform data were generated. Uniform models were fit to transient nonuniform data, and their accuracy was assessed. The uniform models were then used to predict steady-state reaction rates. While the uniform models performed fairly well in fitting the transient nonuniform data, they overpredicted steady-state reaction rates by roughly a factor of 20 regardless of the type of nonuniformity used to generate the data. The insights gained from the model system were used to model the reaction of carbon monoxide and hydrogen to methane and water (methanation). The nonuniform methanation model included adsorbate-adsorbate repulsive interactions and sites created by the adsorption of carbon monoxide upon which only hydrogen could adsorb. The nonuniform model did an excellent job of fitting and predicting transient methanation data. While it could not predict batch reactor methanation rates within an order of magnitude, it did capture the trends in rate with changing temperature and pressure which had eluded other microkinetic models of methanation reported in the literature.

Methylidyne is a key intermediate in hydrocarbon growth reactions over transition metal surfaces. However, experimental data characterizing its interaction with various surfaces are scarce. Therefore, to deepen our understanding of the chemisorption of methylidyne and to quantify its role in the catalytic formation of hydrocarbons, we have calculated the binding energy of methylidyne on Ni(111) and Co(0001) surfaces using density-functional theory within the generalized gradient approximation and the full-potential linear augmented planewave (FP-LAPW) method. The dependence of the binding energy on both the adsorption site and the surface coverage was significant on both surfaces. The binding energy of CH adsorbed in the fcc threefold hollow site of Ni(111) decreased by 1.1 eV as the surface coverage increased from 0.25 to 1.0 ML, whereas the increase of 0.74 eV observed over Co(0001) for the hcp threefold hollow site was slightly less pronounced. The density-of-state plots revealed that rehybridization of the carbon atom interacting with the surface occurred, and the electronic changes induced upon chemisorption were similar for both surfaces.

The Cancer Genome Atlas (TCGA) Research Network has profiled and analyzed large numbers of human tumors to discover molecular aberrations at the DNA, RNA, protein and epigenetic levels. The resulting rich data provide a major opportunity to develop an integrated picture of commonalities, differences and emergent themes across tumor lineages. The Pan-Cancer initiative compares the first 12 tumor types profiled by TCGA. Analysis of the molecular aberrations and their functional roles across tumor types will teach us how to extend therapies effective in one cancer type to others with a similar genomic profile., National Cancer Institute (U.S.), National Human Genome Research Institute (U.S.)

The myelodysplastic syndromes are a group of hematologic disorders that often evolve into secondary acute myeloid leukemia (AML). The genetic changes that underlie progression from the myelodysplastic syndromes to secondary AML are not well understood.We performed whole-genome sequencing of seven paired samples of skin and bone marrow in seven subjects with secondary AML to identify somatic mutations specific to secondary AML. We then genotyped a bone marrow sample obtained during the antecedent myelodysplastic-syndrome stage from each subject to determine the presence or absence of the specific somatic mutations. We identified recurrent mutations in coding genes and defined the clonal architecture of each pair of samples from the myelodysplastic-syndrome stage and the secondary-AML stage, using the allele burden of hundreds of mutations.Approximately 85% of bone marrow cells were clonal in the myelodysplastic-syndrome and secondary-AML samples, regardless of the myeloblast count. The secondary-AML samples contained mutations in 11 recurrently mutated genes, including 4 genes that have not been previously implicated in the myelodysplastic syndromes or AML. In every case, progression to acute leukemia was defined by the persistence of an antecedent founding clone containing 182 to 660 somatic mutations and the outgrowth or emergence of at least one subclone, harboring dozens to hundreds of new mutations. All founding clones and subclones contained at least one mutation in a coding gene.Nearly all the bone marrow cells in patients with myelodysplastic syndromes and secondary AML are clonally derived. Genetic evolution of secondary AML is a dynamic process shaped by multiple cycles of mutation acquisition and clonal selection. Recurrent gene mutations are found in both founding clones and daughter subclones. (Funded by the National Institutes of Health and others.).

The Human Microbiome Project will establish a reference data set for analysis of the microbiome of healthy adults by surveying multiple body sites from 300 people and generating data from over 12,000 samples. To characterize these samples, the participating sequencing centers evaluated and adopted 16S rDNA community profiling protocols for ABI 3730 and 454 FLX Titanium sequencing. In the course of establishing protocols, we examined the performance and error characteristics of each technology, and the relationship of sequence error to the utility of 16S rDNA regions for classification- and OTU-based analysis of community structure. The data production protocols used for this work are those used by the participating centers to produce 16S rDNA sequence for the Human Microbiome Project. Thus, these results can be informative for interpreting the large body of clinical 16S rDNA data produced for this project.

Author Manuscript date: 2013 February 05., Studies of the human microbiome have revealed that even healthy individuals differ remarkably in the microbes that occupy habitats such as the gut, skin and vagina. Much of this diversity remains unexplained, although diet, environment, host genetics and early microbial exposure have all been implicated. Accordingly, to characterize the ecology of human-associated microbial communities, the Human Microbiome Project has analysed the largest cohort and set of distinct, clinically relevant body habitats so far. We found the diversity and abundance of each habitat’s signature microbes to vary widely even among healthy subjects, with strong niche specialization both within and among individuals. The project encountered an estimated 81–99% of the genera, enzyme families and community configurations occupied by the healthy Western microbiome. Metagenomic carriage of metabolic pathways was stable among individuals despite variation in community structure, and ethnic/racial background proved to be one of the strongest associations of both pathways and microbes with clinical metadata. These results thus delineate the range of structural and functional configurations normal in the microbial communities of a healthy population, enabling future characterization of the epidemiology, ecology and translational applications of the human microbiome.

SummaryMost mutations in cancer genomes are thought to be acquired after the initiating event, which may cause genomic instability and drive clonal evolution. However, for acute myeloid leukemia (AML), normal karyotypes are common, and genomic instability is unusual. To better understand clonal evolution in AML, we sequenced the genomes of M3-AML samples with a known initiating event (PML-RARA) versus the genomes of normal karyotype M1-AML samples and the exomes of hematopoietic stem/progenitor cells (HSPCs) from healthy people. Collectively, the data suggest that most of the mutations found in AML genomes are actually random events that occurred in HSPCs before they acquired the initiating mutation; the mutational history of that cell is “captured” as the clone expands. In many cases, only one or two additional, cooperating mutations are needed to generate the malignant founding clone. Cells from the founding clone can acquire additional cooperating mutations, yielding subclones that can contribute to disease progression and/or relapse.

Motivation: The sequencing of tumors and their matched normals is frequently used to study the genetic composition of cancer. Despite this fact, there remains a dearth of available software tools designed to compare sequences in pairs of samples and identify sites that are likely to be unique to one sample. Results: In this article, we describe the mathematical basis of our SomaticSniper software for comparing tumor and normal pairs. We estimate its sensitivity and precision, and present several common sources of error resulting in miscalls. Availability and implementation: Binaries are freely available for download at http://gmt.genome.wustl.edu/somatic-sniper/current/, implemented in C and supported on Linux and Mac OS X. Contact: delarson@wustl.edu; lding@wustl.edu Supplementary information: Supplementary data are available at Bioinformatics online.

To characterize somatic alterations in colorectal carcinoma, we conducted a genome-scale analysis of 276 samples, analysing exome sequence, DNA copy number, promoter methylation and messenger RNA and microRNA expression. A subset of these samples (97) underwent low-depth-of-coverage whole-genome sequencing. In total, 16% of colorectal carcinomas were found to be hypermutated: three-quarters of these had the expected high microsatellite instability, usually with hypermethylation and MLH1 silencing, and one-quarter had somatic mismatch-repair gene and polymerase ε (POLE) mutations. Excluding the hypermutated cancers, colon and rectum cancers were found to have considerably similar patterns of genomic alteration. Twenty-four genes were significantly mutated, and in addition to the expected APC, TP53, SMAD4, PIK3CA and KRAS mutations, we found frequent mutations in ARID1A, SOX9 and FAM123B. Recurrent copy-number alterations include potentially drug-targetable amplifications of ERBB2 and newly discovered amplification of IGF2. Recurrent chromosomal translocations include the fusion of NAV2 and WNT pathway member TCF7L1. Integrative analyses suggest new markers for aggressive colorectal carcinoma and an important role for MYC-directed transcriptional activation and repression., National Institutes of Health (U.S.) (Grant U24CA143799), National Institutes of Health (U.S.) (Grant U24CA143835), National Institutes of Health (U.S.) (Grant U24CA143840), National Institutes of Health (U.S.) (Grant U24CA143843), National Institutes of Health (U.S.) (Grant U24CA143845), National Institutes of Health (U.S.) (Grant U24CA143848), National Institutes of Health (U.S.) (Grant U24CA143858), National Institutes of Health (U.S.) (Grant U24CA143866), National Institutes of Health (U.S.) (Grant U24CA143867), National Institutes of Health (U.S.) (Grant U24CA143882), National Institutes of Health (U.S.) (Grant U24CA143883), National Institutes of Health (U.S.) (Grant U24CA144025), National Institutes of Health (U.S.) (Grant U54HG003067), National Institutes of Health (U.S.) (Grant U54HG003079), National Institutes of Health (U.S.) (Grant U54HG003273)

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders that often progress to chemotherapy-resistant secondary acute myeloid leukemia (sAML). We used whole-genome sequencing to perform an unbiased comprehensive screen to discover the somatic mutations in a sample from an individual with sAML and genotyped the loci containing these mutations in the matched MDS sample. Here we show that a missense mutation affecting the serine at codon 34 (Ser34) in U2AF1 was recurrently present in 13 out of 150 (8.7%) subjects with de novo MDS, and we found suggestive evidence of an increased risk of progression to sAML associated with this mutation. U2AF1 is a U2 auxiliary factor protein that recognizes the AG splice acceptor dinucleotide at the 3' end of introns, and the alterations in U2AF1 are located in highly conserved zinc fingers of this protein. Mutant U2AF1 promotes enhanced splicing and exon skipping in reporter assays in vitro. This previously unidentified, recurrent mutation in U2AF1 implicates altered pre-mRNA splicing as a potential mechanism for MDS pathogenesis.

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.

Summary Genomic structural variants (SVs) are abundant in humans, differing from other variation classes in extent, origin, and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (i.e., copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analyzing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.

The full complement of DNA mutations that are responsible for the pathogenesis of acute myeloid leukemia (AML) is not yet known.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.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.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.

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.

Massively parallel sequencing instruments enable rapid and inexpensive DNA sequence data production. Because these instruments are new, their data require characterization with respect to accuracy and utility. To address this, we sequenced a Caernohabditis elegans N2 Bristol strain isolate using the Solexa Sequence Analyzer, and compared the reads to the reference genome to characterize the data and to evaluate coverage and representation. Massively parallel sequencing facilitates strain-to-reference comparison for genome-wide sequence variant discovery. Owing to the short-read-length sequences produced, we developed a revised approach to determine the regions of the genome to which short reads could be uniquely mapped. We then aligned Solexa reads from C. elegans strain CB4858 to the reference, and screened for single-nucleotide polymorphisms (SNPs) and small indels. This study demonstrates the utility of massively parallel short read sequencing for whole genome resequencing and for accurate discovery of genome-wide polymorphisms.

The Drug-Gene Interaction database (DGIdb) mines existing resources that generate hypotheses about how mutated genes might be targeted therapeutically or prioritized for drug development. It provides an interface for searching lists of genes against a compendium of drug-gene interactions and potentially 'druggable' genes. DGIdb can be accessed at http://dgidb.org/.

Ruthenium has long been known to be an effective catalyst for ammonia synthesis. However, compared to the traditional iron-based catalysts, studies on ruthenium-based catalysts are limited. The rate determining step of ammonia synthesis, the dissociative adsorption of dinitrogen, has been shown to be extremely structure sensitive on both iron and ruthenium catalysts. To study this structure sensitivity on ruthenium, density functional theory calculations were performed on Ru(001) and Ru(110) clusters. End-on, side-on, and dissociated adsorption states were investigated on both surfaces. While the Ru(110) cluster could stabilize all three adsorption modes, a minimum energy structure for the side-on adsorption on Ru(001) could not be found. It is likely that this side-on mode can provide a low energy pathway to the dissociated state, thereby resulting in faster dissociative adsorption on Ru(110).

Background: Estrogen receptors are over-expressed in around 70% of breast cancer cases. The genetic changes that occur during aromatase inhibitor (AI) treatment are not well understood and may differ depending upon the patient's response phenotype. Methods: We performed whole genome sequencing (WGS) of matched blood, pre-treatment, and post-treatment biopsy samples from 22 estrogen receptor positive breast cancer patients treated with neoadjuvant aromatase inhibitors. For 5 cases, we performed the whole genome sequencing (WGS) on patients’ matched normal, two pre AI-treatment, and two post AI-treatment DNA isolates from biopsy samples. We validated all putative coding and non-coding somatic mutations using deep sequencing. By comparing the validated somatic mutations from pre- and post- AI treatment biopsy samples, we were able to determine the alterations in the tumor genomes. In every case we defined the clonal architecture of each pair of pre-treatment and post-treatment biopsy samples by comparing the variant allele frequencies from thousands of validated somatic mutations. Results: Comparisons of the two pre AI-treatment biopsy samples from the same patient indicates that the variant allele frequencies of mutations showed high concordances in all 5 cases, 0.74 to 0.95 range of correlation coefficient. Only a small percentage of somatic mutations were detected in one pre-treatment sample and not the other (4.65% overall). In comparing the somatic variations between pre-treatment and matched post-treatment biopsy samples in 22 cases, we found that patients with good clinical response to AI treatment retained known driver mutations only in their pre-treatment tumors. Conversely, those patients with poor clinical response presented new driver mutations in their post-treatment samples. Furthermore, the variant allele frequency for most mutated genes decreased in post AI treatment samples for patients with good AI treatment response; on the contrary, the variant allele frequency increased for patients with poor clinical response. Conclusions: From WGS of matched normal, pre-treatment, and post-treatment biopsy samples, we identified new driver genes mutated in patients with poor clinical response, while patients with good clinical response had lost mutated driver genes in their post-treatment biopsy samples. The genetic landscape revealed by WGS of pre-treatment and post-treatment biopsy samples reveals mutational repertoires are remodeled by AI therapy. This finding suggests deep sequencing of AI treated samples will be necessary to reveal the complete complement of mutations present in a patient's tumor. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr LB-423. doi:1538-7445.AM2012-LB-423

Abstract 404 To characterize the genomic events associated with distinct subtypes of AML, we used whole genome sequencing to compare 24 tumor/normal sample pairs from patients with normal karyotype (NK) M1-AML (12 cases) and t(15;17)-positive M3-AML (12 cases). All single nucleotide variants (SNVs), small insertions and deletions (indels), and cryptic structural variants (SVs) identified by whole genome sequencing (average coverage 28x) were validated using sample-specific custom Nimblegen capture arrays, followed by Illumina sequencing; an average coverage of 972 reads per somatic variant yielded 10,597 validated somatic variants (average 421/genome). Of these somatic mutations, 308 occurred in 286 unique genes; on average, 9.4 somatic mutations per genome had translational consequences. Several important themes emerged: 1) AML genomes contain a diverse range of recurrent mutations. We assessed the 286 mutated genes for recurrency in an additional 34 NK M1-AML cases and 9 M3-AML cases. We identified 51 recurrently mutated genes, including 37 that had not previously been described in AML; on average, each genome had 3 recurrently mutated genes (M1 = 3.2; M3 = 2.8, p = 0.32). 2) Many recurring mutations cluster in mutually exclusive pathways, suggesting pathophysiologic importance. The most commonly mutated genes were: FLT3 (36%), NPM1 (25%), DNMT3A (21%), IDH1 (18%), IDH2 (10%), TET2 (10%), ASXL1 (6%), NRAS (6%), TTN (6%), and WT1 (6%). In total, 3 genes (excluding PML-RARA) were mutated exclusively in M3 cases. 22 genes were found only in M1 cases (suggestive of alternative initiating mutations which occurred in methylation, signal transduction, and cohesin complex genes). 25 genes were mutated in both M1 and M3 genomes (suggestive of common progression mutations relevant for both subtypes). A single mutation in a cell growth/signaling gene occurred in 38 of 67 cases (FLT3, NRAS, RUNX1, KIT, CACNA1E, CADM2, CSMD1); these mutations were mutually exclusive of one another, and many of them occurred in genomes with PML-RARA, suggesting that they are progression mutations. We also identified a new leukemic pathway: mutations were observed in all four genes that encode members of the cohesin complex (STAG2, SMC1A, SMC3, RAD21), which is involved in mitotic checkpoints and chromatid separation. The cohesin mutations were mutually exclusive of each other, and collectively occur in 10% of non-M3 AML patients. 3) AML genomes also contain hundreds of benign “passenger” mutations. On average 412 somatic mutations per genome were translationally silent or occurred outside of annotated genes. Both M1 and M3 cases had similar total numbers of mutations per genome, similar mutation types (which favored C>T/G>A transitions), and a similar random distribution of variants throughout the genome (which was affected neither by coding regions nor expression levels). This is consistent with our recent observations of random “passenger” mutations in hematopoietic stem cell (HSC) clones derived from normal patients (Ley et al manuscript in preparation), and suggests that most AML-associated mutations are not pathologic, but pre-existed in the HSC at the time of initial transformation. In both studies, the total number of SNVs per genome correlated positively with the age of the patient (R2 = 0.48, p = 0.001), providing a possible explanation for the increasing incidence of AML in elderly patients. 4) NK M1 and M3 AML samples are mono- or oligo-clonal. By comparing the frequency of all somatic mutations within each sample, we could identify clusters of mutations with similar frequencies (leukemic clones) and determined that the average number of clones per genome was 1.8 (M1 = 1.5; M3 = 2.2; p = 0.04). 5) t(15;17) is resolved by a non-homologous end-joining repair pathway, since nucleotide resolution of all 12 t(15;17) breakpoints revealed inconsistent micro-homologies (0 – 7 bp). Summary: These data provide a genome-wide overview of NK and t(15;17) AML and provide important new insights into AML pathogenesis. AML genomes typically contain hundreds of random, non-genic mutations, but only a handful of recurring mutated genes that are likely to be pathogenic because they cluster in mutually exclusive pathways; specific combinations of recurring mutations, as well as rare and private mutations, shape the leukemia phenotype in an individual patient, and help to explain the clinical heterogeneity of this disease. Disclosures: Westervelt: Novartis: Speakers Bureau.

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 patient's 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 patient's 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.

Abstract 99 Whole genome sequencing with next generation technologies represents a new, unbiased approach for discovering somatic variations in cancer genomes. Our group recently reported the DNA sequence and analysis of the genomes of two patients with normal karyotype acute myeloid leukemia (AML). Improvements in next generation sequencing technologies (principally, paired-end sequencing) led us to reevaluate the first case (Ley et al, Nature 456:66–72, 2008) with deeper sequence coverage. We discovered a novel frameshift mutation in DNMT3A, one of the three genes in humans (DNMT1, DNMT3A, and DNMT3B) that encodes a DNA methyltransferase that catalyzes the addition of methyl groups to cytosine within CpG dinucleotides. We then sequenced all the coding exons of this gene in 280 additional de novo cases of AML to define recurring mutations. 62/281 de novo AML cases (22%) had mutations with translational effects in the DNMT3A gene. 18 different missense mutations were identified, the most common of which was at amino acid R882 (37 cases). Frameshifts (n=6), nonsense mutations (n=6), splice site mutations (n=3), and a 1.5 Mbp deletion that included the DNMT3A gene were also identified. DNMT3A mutations were highly enriched in cases with intermediate risk cytogenetics (56/166=33.7%; p Disclosures: Westervelt: Novartis: Honoraria; Celgene: Honoraria, Speakers Bureau. DiPersio:Genzyme: Honoraria.

Abstract 299 Given the poor prognosis of secondary AML evolving from prior MDS and the limited knowledge of mutations that occur during transformation, we set out to comprehensively discover gene mutations that occur during MDS to AML transformation. To achieve this goal, we performed whole genome sequencing (WGS) of paired DNA samples from normal (skin) and tumor (bone marrow) specimens from a patient with MDS-derived AML. The patient was a 65 year old man who presented with pancytopenia and MDS (FAB RA, 4% myeloblasts, del(5q), -17, del(20q)). Three years later he became RBC transfusion dependent and had skin and bone marrow samples banked at Washington University after providing informed consent. His bone marrow biopsy was consistent with RAEB MDS (6% myeloblasts, del(5q), -17, del(20q)), and he subsequently received 4 cycles of decitibine and 2 cycles of lenalidomide before developing AML 2 years later (69% myeloblasts, del(5q), -17, del(20q)). He underwent an allogeneic BMT and died 1 year later. DNA libraries were prepared from the normal (skin) and flow sorted (CD45 dim, low side scatter) AML samples. Using 75–100 bp paired-end reads on the Illumina platform, we generated 126.2 Gb and 101.7 Gb of sequence from the normal and tumor libraries, respectively. Aligned, deduplicated sequence provided 28.7× (normal) and 26.5× (tumor) haploid coverage of the genome, and >96% diploid coverage, using informative SNPs as a metric. Analysis of paired tumor and normal genome sequences allowed us to discriminate between inherited and acquired sequence variants. Indels and structural variants, including copy number alterations, inversions, and translocations were identified using a combination of SAMTools, Pindel, GATK, and Breakdancer algorithms. The cytogenetically visible monosomy 17 and del(20q) lesions were confirmed, and the del(5q) was resolved into 11 distinct interstitial deletions using WGS data. Single nucleotide resolution assembly was possible for 6 of the deleted segments on chromosome 5 and the chromosome 20 deletion which revealed microhomology at the breakpoints, implicating the error prone nonhomologous end joining repair pathway. Interestingly, the tumor suppressor genes APC and DCC were in separate micro-retained regions within larger deletions on chromosome 5, suggesting a selective pressure to retain two copies of these genes during AML progression. Of the 4.8M single nucleotide variants (SNVs) detected in the tumor genome, 33,006 were not found in databases or this patient's normal genome and are, therefore, potential somatic mutations. We identified 46 tier 1 high confidence (HC) SNVs (coding and splice site variants), 388 tier 2 HC SNVs (conserved non-coding variants), 2,185 tier 3 HC SNVs (variants in non-repetitive regions), and 5,862 tier 4 HC SNVs (all other variants). Tier 1 SNVs were prioritized for validation given their potential to be functionally significant. 22/46 Tier 1 HC SNVs were validated as acquired mutations in the AML sample using the 454 platform (5 synonymous, 1 nonsense, and 16 missense mutations). Only 1/22 SNVs has previously been found in myeloid cancers (WT1). The mutant allele frequency of the 22 mutations ranged from 9–100% in the tumor sample based on the 454 deep read counts, suggesting that the underlying heterogeneity that is common in de novo MDS persists during transformation to AML. By sequencing the 22 somatic SNVs in the MDS sample, we identified 10 mutations that were present in both the MDS and AML samples and 12 AML-specific mutations. The latter are likely to be important for MDS to AML evolution. By analyzing the mutant allele frequencies of all 22 SNVs, we predict that 1 dominant clone and 2 sub-clones coexist in the AML sample. All 22 genes with mutations are expressed in the AML sample but none of the 22 SNVs have recurrent mutations at the same nucleotide position in an additional 150 de novo MDS samples. Finally, the residual non-deleted allele of ACTR5, located in the 2.5Mb minimally deleted region (MDR) on chromosome 20, is mutated during AML transformation, implicating it as a potential 20q tumor suppressor gene. Collectively, analysis of WGS data identified mutations in 21 genes not previously implicated in MDS or AML, identified 12 genes as potential drivers of evolution from MDS to AML, allowed for fine mapping of deletion breakpoints on chromosomes 5 and 20, and identified a potential tumor suppressor on the del(20q) MDR. Disclosures: No relevant conflicts of interest to declare.

Abstract 580 Therapy-related acute myeloid leukemia/myelodysplasia (t-AML/t-MDS) accounts for 10–20% of all cases of AML, and its incidence is rising. Treatment options are limited and the prognosis very poor, highlighting the need for new therapies in t-AML/t-MDS. However, the genetic mutations contributing to transformation in t-AML/t-MDS are largely unknown, limiting the development of novel targeted therapeutics. Our group previously reported the sequence of the first two cancer genomes, both in patients with de novo AML (Nature 456:66, 2008; NEJM 361:1058, 2009). Herein, we report the sequence of the cancer genome of a patient with t-AML. The patient presented with early-onset breast, then ovarian cancer (age 96% of heterozygous SNPs were detected. Similar data were obtained for the skin sample. A total of 27 validated somatic single nucleotide variants or indels were detected in coding sequences. None of these mutations have been previously reported in de novo AML. Eight novel chromosomal translocations were identified and the breakpoints defined. One translocation, t(3;4)(q27.3;p15.32), resulted in the production of an in frame fusion transcript of DGKG (diacylglycerol kinase gamma) with BST1 (bone marrow cell stromal antigen 1). Studies are underway to characterize the effect of this fusion gene on hematopoietic cell growth and differentiation. In addition to -7 and del(5q), somatic copy number alterations on chromosome 3 and 12 were identified. There is controversy whether haploinsufficiency of genes on chromosomes 7 and 5q is sufficient to contribute to transformation, or whether further mutations lead to loss of heterozygosity of one or more genes in these regions. In the present case, careful review of the sequence and array data revealed no ‘homozygous' somatic single nucleotide variants, indels, or copy number alterations of coding genes on the remaining copy of chromosome 5 or 7. The patient's clinical presentation strongly suggested genetic cancer susceptibility. Analysis of the skin genome of this patient identified a heterozygous deletion of exons 7–9 of TP53, likely contributing to the early onset of her breast cancer; a uniparental disomy event resulted in the deletion being homozygous in the leukemia sample. Interestingly, the mutant TP53 allele is expressed, and it is predicted to produce a truncated p53 protein lacking most of its DNA binding domain. Functional studies of the mutant p53 protein are underway. Of note, based on a detailed family history and genotyping of the patient's mother, we suspect that the TP53 deletion occurred spontaneously. Ongoing whole genome sequencing studies in a large number of t-AML samples should identify novel somatic mutations and germline variants that contribute to t-AML. Disclosures: No relevant conflicts of interest to declare.

The genetic events responsible for the initiation and progression of AML are unknown for most patients. Although many important mutations for pathogenesis have been discovered, unbiased genomic approaches will be required to discover all relevant mutations. The development of “Next Generation” massively parallel sequencing approaches (454 and Solexa) have substantially reduced the cost of whole genome sequencing, and have allowed us to initiate a study designed to identify all potentially relevant mutations in a case of M1 AML. The first patient selected for study was a previously healthy 57-year-old Caucasian female who presented with a white blood count of 102,500 (91% blasts); the bone marrow biopsy revealed 100% blasts. Cytogenetics were normal. Tumor and skin (normal) samples were banked with informed consent, using our Washington University IRB-approved protocol that specifically permits whole genome sequencing. Resequencing of 14 genes frequently mutated in AML revealed somatic mutations in FLT3 (ITD) and NPM. Oligomer array-based comparative genomic hybridization using the 2.1 million-oligomer NimbleGen array revealed only one small (

// Justina McEvoy 1,* , Panduka Nagahawatte 2,* , David Finkelstein 2 , Jennifer Richards-Yutz 6 , Marcus Valentine 13 , Jing Ma 14 , Charles Mullighan 14 , Guangchun Song 14 , Xiang Chen 2 , Matthew Wilson 4 , Rachel Brennan 12 , Stanley Pounds 3 , Jared Becksfort 2 , Robert Huether 2 , Charles Lu 7 , Robert S. Fulton 7,8 , Lucinda L. Fulton 7,8 , Xin Hong 7,8 , David J. Dooling 7,8 , Kerri Ochoa 7,8 , Elaine R. Mardis 7,8,9 , Richard K.Wilson 7,8,10 , John Easton 2 , Jinghui Zhang 2 , James R. Downing 14 , Arupa Ganguly 5,6,* and Michael A. Dyer 1,4,11 for the St. Jude Children’s Research Hospital – Washington University Pediatric Cancer Genome Project 1 Departments of Developmental Neurobiology, St. Jude Children’s Research Hospital, Memphis, TN, USA. 2 Computational Biology and Bioinformatics, St. Jude Children’s Research Hospital, Memphis, TN, USA. 3 Biostatistics, St. Jude Children’s Research Hospital, Memphis, TN, USA. 4 Department of Ophthalmology, University of Tennessee Health Science Center, Memphis, TN 5 Department of Genetics, School of Medicine, Philadelphia, PA, USA 6 Genetic Diagnostic Laboratory at University of Pennsylvania, School of Medicine, Philadelphia, PA, USA. 7 The Genome Institute, Washington University School of Medicine in St Louis, St Louis, Missouri, USA. 8 Department of Genetics, Washington University School of Medicine in St Louis, St Louis, Missouri, USA. 9 Siteman Cancer Center, Washington University School of Medicine in St Louis, St Louis, Missouri, USA. 10 Department of Medicine, Washington University School of Medicine in St Louis, St Louis, Missouri, USA 11 Howard Hughes Medical Institute, Chevy Chase, MD 12 Oncology, St. Jude Children’s Research Hospital, Memphis, TN, USA. 13 Cytogenetics, St. Jude Children’s Research Hospital, Memphis, TN, USA. 14 Pathology, St. Jude Children’s Research Hospital, Memphis, TN, USA. * These authors contributed equally to this work. Correspondence: Michael A. Dyer , email: // Keywords : chromothripsis, retinoblastoma, RB1, MYCN Received : December 13, 2013 Accepted : January 7, 2014 Published : January 11, 2014 Abstract Retinoblastoma is a rare childhood cancer of the developing retina. Most retinoblastomas initiate with biallelic inactivation of the RB1 gene through diverse mechanisms including point mutations, nucleotide insertions, deletions, loss of heterozygosity and promoter hypermethylation. Recently, a novel mechanism of retinoblastoma initiation was proposed. Gallie and colleagues discovered that a small proportion of retinoblastomas lack RB1 mutations and had MYCN amplification [ 1 ]. In this study, we identified recurrent chromosomal, regional and focal genomic lesions in 94 primary retinoblastomas with their matched normal DNA using SNP 6.0 chips. We also analyzed the RB1 gene mutations and compared the mechanism of RB1 inactivation to the recurrent copy number variations in the retinoblastoma genome. In addition to the previously described focal amplification of MYCN and deletions in RB1 and BCOR , we also identified recurrent focal amplification of OTX2 , a transcription factor required for retinal photoreceptor development. We identified 10 retinoblastomas in our cohort that lacked RB1 point mutations or indels. We performed whole genome sequencing on those 10 tumors and their corresponding germline DNA. In one of the tumors, the RB1 gene was unaltered, the MYCN gene was amplified and RB1 protein was expressed in the nuclei of the tumor cells. In addition, several tumors had complex patterns of structural variations and we identified 3 tumors with chromothripsis at the RB1 locus. This is the first report of chromothripsis as a mechanism for RB1 gene inactivation in cancer.

Abstract 68 Early T-cell precursor acute lymphoblastic leukemia (ETP ALL) is characterized by an immature T-lineage immunophenotype (cCD3+, CD1a-, CD8- and CD5dim) aberrant expression of myeloid and stem cell markers, a distinct gene expression profile and very poor outcome. The underlying genetic basis of this form of leukemia is unknown. Here we report results of whole genome sequencing (WGS) of tumor and normal DNA from 12 children with ETP ALL. Genomes were sequenced to 30-fold haploid coverage using the Illumina GAIIx platform, and all putative somatic sequence and structural variants were validated. The frequency of mutations in 43 genes was assessed in a recurrence cohort of 52 ETP and 42 non-ETP T-ALL samples from patients enrolled in St Jude, Children's Oncology Group and AEIOP trials. Transcriptomic resequencing was performed for two WGS cases, and whole exome sequencing for three ETP ALL cases in the recurrence cohort. We identified 44 interchromosomal translocations (mean 4 per patient, range 0–12), 32 intrachromosomal translocations (mean 3, 0–7), 53 deletions (mean 4, 0–10) and 16 insertions (mean 1, 0–5). Three cases exhibited a pattern of complex rearrangements suggestive of a single cellular catastrophe (“chromothripsis”), two of which had mutations targeting mismatch and DNA repair (MLH3 and DCLRE1C). While no single chromosomal alteration was present in all cases, 10 of 12 ETP ALLs harbored chromosomal rearrangements, several of which involved complex multichromosomal translocations and resulted in the expression of chimeric in-frame novel fusion genes disrupting hematopoietic regulators, including ETV6-INO80D, NAP1L1-MLLT10, RUNX1-EVX1 and NUP214-SQSTM1, each occurring in a single case. An additional ETP case with the ETV6-INO80D fusion was identified in the recurrence cohort. Additionally, 51% of structural variants had breakpoints in genes, including those with roles in hematopoiesis and leukemogenesis, and genes also targeted by mutation in other cases (MLH3, SUZ12, RUNX1). We identified a high frequency of activating mutations in genes regulating cytokine receptor and Ras signalling in ETP ALL (67.2% of ETP compared to 19% of non-ETP T-ALL) including NRAS (17%), FLT3 (14%), JAK3 (9%), SH2B3 (or LNK; 9%), IL7R (8%), JAK1 (8%), KRAS (3%), and BRAF (2%). Seven cases (5 ETP, 2 non-ETP) harbored in frame insertion mutations in the transmembrane domain of IL7R, which were transforming when expressed in the murine cell lines, and resulted in enhanced colony formation when expressed in primary murine hematopoietic cells. The IL7R mutations resulted in constitutive Jak-Stat activation in these cell lines and primary leukemic cells expressing these mutations. Fifty-eight percent of ETP cases (compared to 17% of non-ETP cases) harbored mutations known or predicted to disrupt hematopoietic and lymphoid development, including ETV6 (33%), RUNX1 (16%), IKZF1 (14%), GATA3 (10%), EP300 (5%) and GATA2 (2%). GATA3 regulates early T cell development, and mutations in this gene were observed exclusively in ETP ALL. The mutations were commonly biallelic, and were clustered at R276, a residue critical for binding of GATA3 to DNA. Strikingly, mutations disrupting chromatin modifying genes were also highly enriched in ETP ALL. Genes encoding the the polycomb repressor complex 2 (EZH2, SUZ12 and EED), that mediates histone 3 lysine 27 (H3K27) trimethylation were deleted or mutated in 42% of ETP ALL compared to 12% of non-ETP T-ALL. In addition, alterations of the H3K36 trimethylase SETD2 were observed in 5 ETP cases, but not in non-ETP ALL. We also identified recurrent mutations in genes that have not previously been implicated in hematopoietic malignancies including RELN, DNM2, ECT2L, HNRNPA1 and HNRNPR. Using gene set enrichment analysis we demonstrate that the gene expression profile of ETP ALL shares features not only with normal human hematopoietic stem cells, but also with leukemic initiating cells (LIC) purified from patients with acute myeloid leukemia (AML). These results indicate that mutations that drive proliferation, impair differentiation and disrupt histone modification cooperate to induce an aggressive leukemia with an aberrant immature phenotype. The similarity of the gene expression pattern with that observed in the LIC of AML raises the possibility that myeloid-directed therapies might improve the outcome of ETP ALL. Disclosures: Evans: St. Jude Children's research Hospital: Employment, Patents & Royalties; NIH & NCI: Research Funding; Aldagen: Membership on an entity's Board of Directors or advisory committees.

We analysed primary breast cancers by genomic DNA copy number arrays, DNA methylation, exome sequencing, messenger RNA arrays, microRNA sequencing and reverse-phase protein arrays. Our ability to integrate information across platforms provided key insights into previously defined gene expression subtypes and demonstrated the existence of four main breast cancer classes when combining data from five platforms, each of which shows significant molecular heterogeneity. Somatic mutations in only three genes (TP53, PIK3CA and GATA3) occurred at >10% incidence across all breast cancers; however, there were numerous subtype-associated and novel gene mutations including the enrichment of specific mutations in GATA3, PIK3CA and MAP3K1 with the luminal A subtype. We identified two novel protein-expression-defined subgroups, possibly produced by stromal/microenvironmental elements, and integrated analyses identified specific signalling pathways dominant in each molecular subtype including a HER2/phosphorylated HER2/EGFR/phosphorylated EGFR signature within the HER2-enriched expression subtype. Comparison of basal-like breast tumours with high-grade serous ovarian tumours showed many molecular commonalities, indicating a related aetiology and similar therapeutic opportunities. The biological finding of the four main breast cancer subtypes caused by different subsets of genetic and epigenetic abnormalities raises the hypothesis that much of the clinically observable plasticity and heterogeneity occurs within, and not across, these major biological subtypes of breast cancer.

The unparalleled efficiency of next-generation sequencing (NGS) has prompted widespread adoption, but significant problems remain in the use of NGS data for whole genome assembly. We explore the advantages and disadvantages of chicken genome assemblies generated using a variety of sequencing and assembly methodologies. NGS assemblies are equivalent in some ways to a Sanger-based assembly yet deficient in others. Nonetheless, these assemblies are sufficient for the identification of the majority of genes and can reveal novel sequences when compared to existing assembly references.