Person:

Beroukhim, Rameen

Loss of heterozygosity (LOH) of chromosomal regions bearing tumor suppressors is a key event in the evolution of epithelial and mesenchymal tumors. Identification of these regions usually relies on genotyping tumor and counterpart normal DNA and noting regions where heterozygous alleles in the normal DNA become homozygous in the tumor. However, paired normal samples for tumors and cell lines are often not available. With the advent of oligonucleotide arrays that simultaneously assay thousands of single-nucleotide polymorphism (SNP) markers, genotyping can now be done at high enough resolution to allow identification of LOH events by the absence of heterozygous loci, without comparison to normal controls. Here we describe a hidden Markov model-based method to identify LOH from unpaired tumor samples, taking into account SNP intermarker distances, SNP-specific heterozygosity rates, and the haplotype structure of the human genome. When we applied the method to data genotyped on 100 K arrays, we correctly identified 99% of SNP markers as either retention or loss. We also correctly identified 81% of the regions of LOH, including 98% of regions greater than 3 megabases. By integrating copy number analysis into the method, we were able to distinguish LOH from allelic imbalance. Application of this method to data from a set of prostate samples without paired normals identified known regions of prevalent LOH. We have developed a method for analyzing high-density oligonucleotide SNP array data to accurately identify of regions of LOH and retention in tumors without the need for paired normal samples.

Amplification, deletion, and loss of heterozygosity of genomic DNA are hallmarks of cancer. In recent years a variety of studies have emerged measuring total chromosomal copy number at increasingly high resolution. Similarly, loss-of-heterozygosity events have been finely mapped using high-throughput genotyping technologies. We have developed a probe-level allele-specific quantitation procedure that extracts both copy number and allelotype information from single nucleotide polymorphism (SNP) array data to arrive at allele-specific copy number across the genome. Our approach applies an expectation-maximization algorithm to a model derived from a novel classification of SNP array probes. This method is the first to our knowledge that is able to (a) determine the generalized genotype of aberrant samples at each SNP site (e.g., CCCCT at an amplified site), and (b) infer the copy number of each parental chromosome across the genome. With this method, we are able to determine not just where amplifications and deletions occur, but also the haplotype of the region being amplified or deleted. The merit of our model and general approach is demonstrated by very precise genotyping of normal samples, and our allele-specific copy number inferences are validated using PCR experiments. Applying our method to a collection of lung cancer samples, we are able to conclude that amplification is essentially monoallelic, as would be expected under the mechanisms currently believed responsible for gene amplification. This suggests that a specific parental chromosome may be targeted for amplification, whether because of germ line or somatic variation. An R software package containing the methods described in this paper is freely available at http://genome.dfci.harvard.edu/~tlaframb/PLASQ.

Background: Single nucleotide polymorphisms (SNPs) are the most common genetic variations in the human genome and are useful as genomic markers. Oligonucleotide SNP microarrays have been developed for high-throughput genotyping of up to 900,000 human SNPs and have been used widely in linkage and cancer genomics studies. We have previously used Hidden Markov Models (HMM) to analyze SNP array data for inferring copy numbers and loss-of-heterozygosity (LOH) from paired normal and tumor samples and unpaired tumor samples. Results: We proposed and implemented major copy proportion (MCP) analysis of oligonucleotide SNP array data. A HMM was constructed to infer unobserved MCP states from observed allele-specific signals through emission and transition distributions. We used 10 K, 100 K and 250 K SNP array datasets to compare MCP analysis with LOH and copy number analysis, and showed that MCP performs better than LOH analysis for allelic-imbalanced chromosome regions and normal contaminated samples. The major and minor copy alleles can also be inferred from allelic-imbalanced regions by MCP analysis. Conclusion: MCP extends tumor LOH analysis to allelic imbalance analysis and supplies complementary information to total copy numbers. MCP analysis of mixing normal and tumor samples suggests the utility of MCP analysis of normal-contaminated tumor samples. The described analysis and visualization methods are readily available in the user-friendly dChip software.

Background: Protein tyrosine kinases are important regulators of cellular homeostasis with tightly controlled catalytic activity. Mutations in kinase-encoding genes can relieve the autoinhibitory constraints on kinase activity, can promote malignant transformation, and appear to be a major determinant of response to kinase inhibitor therapy. Missense mutations in the EGFR kinase domain, for example, have recently been identified in patients who showed clinical responses to EGFR kinase inhibitor therapy. Methods and Findings: Encouraged by the promising clinical activity of epidermal growth factor receptor (EGFR) kinase inhibitors in treating glioblastoma in humans, we have sequenced the complete EGFR coding sequence in glioma tumor samples and cell lines. We identified novel missense mutations in the extracellular domain of EGFR in 13.6% (18/132) of glioblastomas and 12.5% (1/ 8) of glioblastoma cell lines. These EGFR mutations were associated with increased EGFR gene dosage and conferred anchorage-independent growth and tumorigenicity to NIH-3T3 cells. Cells transformed by expression of these EGFR mutants were sensitive to small-molecule EGFR kinase inhibitors. Conclusions: Our results suggest extracellular missense mutations as a novel mechanism for oncogenic EGFR activation and may help identify patients who can benefit from EGFR kinase inhibitors for treatment of glioblastoma.

Chromosome 8q24 is the most commonly amplified region across multiple cancer types, and the typical length of the amplification suggests that it may target additional genes to MYC. To explore the roles of the genes most frequently included in 8q24 amplifications, we analyzed the relation between copy number alterations and gene expression in three sets of endometrial cancers (N = 252); and in glioblastoma, ovarian, and breast cancers profiled by TCGA. Among the genes neighbouring MYC, expression of the bromodomain-containing gene ATAD2 was the most associated with amplification. Bromodomain-containing genes have been implicated as mediators of MYC transcriptional function, and indeed ATAD2 expression was more closely associated with expression of genes known to be upregulated by MYC than was MYC itself. Amplifications of 8q24, expression of genes downstream from MYC, and overexpression of ATAD2 predicted poor outcome and increased from primary to metastatic lesions. Knockdown of ATAD2 and MYC in seven endometrial and 21 breast cancer cell lines demonstrated that cell lines that were dependent on MYC also depended upon ATAD2. These same cell lines were also the most sensitive to the histone deacetylase (HDAC) inhibitor Trichostatin-A, consistent with prior studies identifying bromodomain-containing proteins as targets of inhibition by HDAC inhibitors. Our data indicate high ATAD2 expression is a marker of aggressive endometrial cancers, and suggest specific inhibitors of ATAD2 may have therapeutic utility in these and other MYC-dependent cancers.

Human cancer is caused by the accumulation of genetic alterations in cells. Of special importance are changes that occur early during malignant transformation because they may result in oncogene addiction and thus represent promising targets for therapeutic intervention. We have previously described a computational approach, called Retracing the Evolutionary Steps in Cancer (RESIC), to determine the temporal sequence of genetic alterations during tumorigenesis from cross-sectional genomic data of tumors at their fully transformed stage. Since alterations within a set of genes belonging to a particular signaling pathway may have similar or equivalent effects, we applied a pathway-based systems biology approach to the RESIC methodology. This method was used to determine whether alterations of specific pathways develop early or late during malignant transformation. When applied to primary glioblastoma (GBM) copy number data from The Cancer Genome Atlas (TCGA) project, RESIC identified a temporal order of pathway alterations consistent with the order of events in secondary GBMs. We then further subdivided the samples into the four main GBM subtypes and determined the relative contributions of each subtype to the overall results: we found that the overall ordering applied for the proneural subtype but differed for mesenchymal samples. The temporal sequence of events could not be identified for neural and classical subtypes, possibly due to a limited number of samples. Moreover, for samples of the proneural subtype, we detected two distinct temporal sequences of events: (i) RAS pathway activation was followed by TP53 inactivation and finally PI3K2 activation, and (ii) RAS activation preceded only AKT activation. This extension of the RESIC methodology provides an evolutionary mathematical approach to identify the temporal sequence of pathway changes driving tumorigenesis and may be useful in guiding the understanding of signaling rearrangements in cancer development.

Aberrant NF-κB activation is frequently observed in human cancers. Genome characterization efforts have identified genetic alterations in multiple components of the NF-κB pathway, some of which have been shown to be essential for cancer initiation and tumor maintenance. Here using patient tumors and cancer cell lines, we identify the NF-κB regulator, TRAF2 as an oncogene that is recurrently amplified and rearranged in 15% of human epithelial cancers. Suppression of TRAF2 in cancer cells harboring TRAF2 copy number gain inhibits proliferation, NF-κB activation, anchorage-independent growth and tumorigenesis. Cancer cells that are dependent on TRAF2 also require NF-κB for survival. The phosphorylation of TRAF2 at serine 11 is essential for the survival of cancer cells harboring TRAF2 amplification. Together these observations identify TRAF2 as a frequently amplified oncogene.

The incidence of esophageal adenocarcinoma (EAC) has risen 600% over the last 30 years. With a five-year survival rate of 15%, identification of new therapeutic targets for EAC is greatly important. We analyze the mutation spectra from whole exome sequencing of 149 EAC tumors/normal pairs, 15 of which have also been subjected to whole genome sequencing. We identify a mutational signature defined by a high prevalence of A to C transversions at AA dinucleotides. Statistical analysis of exome data identified significantly mutated 26 genes. Of these genes, four (TP53, CDKN2A, SMAD4, and PIK3CA) have been previously implicated in EAC. The novel significantly mutated genes include chromatin modifying factors and candidate contributors: SPG20, TLR4, ELMO1, and DOCK2. Functional analyses of EAC-derived mutations in ELMO1 reveal increased cellular invasion. Therefore, we suggest a new hypothesis about the potential activation of the RAC1 pathway to be a contributor to EAC tumorigenesis.

Meningiomas are the most common primary nervous system tumor. The tumor suppressor NF2 is disrupted in approximately half of meningiomas1 but the complete spectrum of genetic changes remains undefined. We performed whole-genome or whole-exome sequencing on 17 meningiomas and focused sequencing on an additional 48 tumors to identify and validate somatic genetic alterations. Most meningiomas exhibited simple genomes, with fewer mutations, rearrangements, and copy-number alterations than reported in other adult tumors. However, several meningiomas harbored more complex patterns of copy-number changes and rearrangements including one tumor with chromothripsis. We confirmed focal NF2 inactivation in 43% of tumors and found alterations in epigenetic modifiers among an additional 8% of tumors. A subset of meningiomas lacking NF2 alterations harbored recurrent oncogenic mutations in AKT1 (E17K) and SMO (W535L) and exhibited immunohistochemical evidence of activation of their pathways. These mutations were present in therapeutically challenging tumors of the skull base and higher grade. These results begin to define the spectrum of genetic alterations in meningiomas and identify potential therapeutic targets.