Large Number Of Somatic Mutations Research Articles

Abstract 525: Higher mutation rate is associated with more frequent expression of cancer/testis antigens in human tumors

Abstract Tumors with high mutation rate such as those of melanoma and non-small cell lung cancers (NSCLC) have shown greatest response to the emerging immune therapies. These tumors harbor a large number of somatic mutations, some of which could be translated into neoantigens, potentially evoking host immune response. To what degree self antigens, many of which are cancer testis (CT) antigens, contribute to enhanced immune response in these tumors has never been systematically studied. Here we show that cancers with higher mutation rate frequently express a larger number of CT antigens than cancers with fewer mutations. In NSCLC, tumors with more than 300 non-synonymous mutations (75th percentile, n = 158) express more than 3 times as many CT antigens as those tumors with fewer than 100 mutations (25th percentile, n = 162), with a P-value of 1.6e-13. In head and neck cancers, tumors with more than 140 non-synonymous mutations (75th percentile, n = 79) express nearly twice as many CT antigens as those tumors with fewer than 60 mutations (25th percentile, n = 80), with a P-value of 0.004. Additionally, in NSCLC, tumors with high mutation rate are associated with higher expression of CD8B and granzyme B, markers representing cytotoxic lymphocytes. Tumors with more expressed CT antigens are also associated with higher expression of these markers. These data suggest that CT antigens, along with neoantigens, may contribute to a stronger host immune response in these hypermutative tumors. Citation Format: Junping Jing, Patrick Mayes, Yan Degenhardt, James R. Brown, Philippe Sanseau. Higher mutation rate is associated with more frequent expression of cancer/testis antigens in human tumors. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 525.

Merkel Cell Carcinoma Therapeutic Update.

Merkel cell carcinoma (MCC) is a rare and aggressive neuroendocrine tumor of the skin. Early-stage disease can be cured with surgical resection and radiotherapy (RT). Sentinel lymph node biopsy (SLNB) is an important staging tool, as a microscopic MCC is frequently identified. Adjuvant RT to the primary excision site and regional lymph node bed may improve locoregional control. However, newer studies confirm that patients with biopsy-negative sentinel lymph nodes may not benefit from regional RT. Advanced MCC currently lacks a highly effective treatment as responses to chemotherapy are not durable. Recent work suggests that immunotherapy targeting the programmed cell death receptor 1/programmed cell death ligand 1 (PD-1/PD-L1) checkpoint holds great promise in treating advanced MCC and may provide durable responses in a portion of patients. At the same time, high-throughput sequencing studies have demonstrated significant differences in the mutational profiles of tumors with and without the Merkel cell polyomavirus (MCV). An important secondary endpoint in the ongoing immunotherapy trials for MCC will be determining if there is a response difference between the virus-positive MCC tumors that typically lack a large mutational burden and the virus-negative tumors that have a large number of somatic mutations and predicted tumor neoantigens. Interestingly, sequencing studies have failed to identify a highly recurrent activated driver pathway in the majority of MCC tumors. This may explain why targeted therapies can demonstrate exceptional responses in case reports but fail when treating all comers with MCC. Ultimately, a precision medicine approach may be more appropriate for treating MCC, where identified driver mutations are used to direct targeted therapies. At a minimum, stratifying patients in future clinical trials based on tumor viral status should be considered as virus-negative tumors are more likely to harbor activating driver mutations.

Low frequency of broadly neutralizing HIV antibodies during chronic infection even in quaternary epitope targeting antibodies containing large numbers of somatic mutations.

Neutralizing antibodies (Abs) are thought to be a critical component of an appropriate HIV vaccine response. It has been proposed that Abs recognizing conformationally dependent quaternary epitopes on the HIV envelope (Env) trimer may be necessary to neutralize diverse HIV strains. A number of recently described broadly neutralizing monoclonal Abs (mAbs) recognize complex and quaternary epitopes. Generally, many such Abs exhibit extensive numbers of somatic mutations and unique structural characteristics. We sought to characterize the native antibody (Ab) response against circulating HIV focusing on such conformational responses, without a prior selection based on neutralization. Using a capture system based on VLPs incorporating cleaved envelope protein, we identified a selection of B cells that produce quaternary epitope targeting Abs (QtAbs). Similar to a number of broadly neutralizing Abs, the Ab genes encoding these QtAbs showed extensive numbers of somatic mutations. However, when expressed as recombinant molecules, these Abs failed to neutralize virus or mediate ADCVI activity. Molecular analysis showed unusually high numbers of mutations in the Ab heavy chain framework 3 region of the variable genes. The analysis suggests that large numbers of somatic mutations occur in Ab genes encoding HIV Abs in chronically infected individuals in a non-directed, stochastic, manner.

Mining the coding and non-coding genome for cancer drivers

Progress in next-generation sequencing provides unprecedented opportunities to fully characterize the spectrum of somatic mutations of cancer genomes. Given the large number of somatic mutations identified by such technologies, the prioritization of cancer-driving events is a consistent bottleneck. Most bioinformatics tools concentrate on driver mutations in the coding fraction of the genome, those causing changes in protein products. As more non-coding pathogenic variants are identified and characterized, the development of computational approaches to effectively prioritize cancer-driving variants within the non-coding fraction of human genome is becoming critical. After a short summary of methods for coding variant prioritization, we here review the highly diverse non-coding elements that may act as cancer drivers and describe recent methods that attempt to evaluate the deleteriousness of sequence variation in these elements. With such tools, the prioritization and identification of cancer-implicated regulatory elements and non-coding RNAs is becoming a reality.

Whole-exome sequencing of fibroblast and its iPS cell lines derived from a patient diagnosed with xeroderma pigmentosum.

Cells from a patient with a DNA repair-deficiency disorder are anticipated to bear a large number of somatic mutations. Because such mutations occur independently in each cell, there is a high degree of mosaicism in patients' tissues. While major mutations that have been expanded in many cognate cells are readily detected by sequencing, minor ones are overlaid with a large depth of non-mutated alleles and are not detected. However, cell cloning enables us to observe such cryptic mutations as well as major mutations. In the present study, we focused on a fibroblastic cell line that is derived from a patient diagnosed with xeroderma pigmentosum (XP), which is an autosomal recessive disorder caused by a deficiency in nucleotide excision repair. By making a list of somatic mutations, we can expect to see a characteristic pattern of mutations caused by the hereditary disorder. We cloned a cell by generating an iPS cell line and performed a whole-exome sequencing analysis of the progenitor and its iPS cell lines. Unexpectedly, we failed to find causal mutations in the XP-related genes, but we identified many other mutations including homozygous deletion of GSTM1 and GSTT1. In addition, we found that the long arm of chromosome 9 formed uniparental disomy in the iPS cell line, which was also confirmed by a structural mutation analysis using a SNP array. Type and number of somatic mutations were different from those observed in XP patients. Taken together, we conclude that the patient might be affected by a different type of the disorder and that some of the mutations that we identified here may be responsible for exhibiting the phenotype. Sequencing and SNP-array data have been submitted to SRA and GEO under accession numbers SRP059858 and GSE55520, respectively.

PD-1 Blockade in Tumors with Mismatch-Repair Deficiency

BackgroundSomatic mutations have the potential to encode “non-self” immunogenic antigens. We hypothesized that tumors with a large number of somatic mutations due to mismatch-repair defects may be susceptible to immune checkpoint blockade.MethodsWe conducted a phase 2 study to evaluate the clinical activity of pembrolizumab, an anti–programmed death 1 immune checkpoint inhibitor, in 41 patients with progressive metastatic carcinoma with or without mismatch-repair deficiency. Pembrolizumab was administered intravenously at a dose of 10 mg per kilogram of body weight every 14 days in patients with mismatch repair–deficient colorectal cancers, patients with mismatch repair–proficient colorectal cancers, and patients with mismatch repair–deficient cancers that were not colorectal. The coprimary end points were the immune-related objective response rate and the 20-week immune-related progression-free survival rate.ResultsThe immune-related objective response rate and immune-related progression-free survival rate were 40% (4 of 10 patients) and 78% (7 of 9 patients), respectively, for mismatch repair–deficient colorectal cancers and 0% (0 of 18 patients) and 11% (2 of 18 patients) for mismatch repair–proficient colorectal cancers. The median progression-free survival and overall survival were not reached in the cohort with mismatch repair–deficient colorectal cancer but were 2.2 and 5.0 months, respectively, in the cohort with mismatch repair–proficient colorectal cancer (hazard ratio for disease progression or death, 0.10 [P<0.001], and hazard ratio for death, 0.22 [P=0.05]). Patients with mismatch repair–deficient noncolorectal cancer had responses similar to those of patients with mismatch repair–deficient colorectal cancer (immune-related objective response rate, 71% [5 of 7 patients]; immune-related progression-free survival rate, 67% [4 of 6 patients]). Whole-exome sequencing revealed a mean of 1782 somatic mutations per tumor in mismatch repair–deficient tumors, as compared with 73 in mismatch repair–proficient tumors (P=0.007), and high somatic mutation loads were associated with prolonged progression-free survival (P=0.02).ConclusionsThis study showed that mismatch-repair status predicted clinical benefit of immune checkpoint blockade with pembrolizumab. (Funded by Johns Hopkins University and others; ClinicalTrials.gov number, NCT01876511.)

PD-1 blockade in tumors with mismatch-repair deficiency

Background Somatic mutations have the potential to encode “non-self” immunogenic antigens. We hypothesized that tumors with a large number of somatic mutations due to mismatch-repair defects may be susceptible to immune checkpoint blockade.

Abstract 5350: Aberrant activation-induced cytidine deaminase (AID) expression is induced in the process of hepatitis C virus-associated hepatocarcinogenesis

Abstract Activation-induced cytidine deaminase (AID) is a DNA mutator enzyme, which is encoded by Aicda gene. Under physiological condition, AID is essential for somatic hypermutation of immunogloburin genes in activated B lymphocytes and involved in immune diversity. Previous in vitro studies demonstrated that in cultured adult hepatocytes, both the response to chronic inflammation and hepatitis C virus (HCV) infection itself triggers aberrant AID upregulation, and due to the mutagenic activity of AID, nucleotide alterations gradually accumulate in tumor-related genes. In this study, we aimed to obtain a direct evidence of AID expression induced by HCV infection in vivo using human and mice subjects. Immunostaining of human clinical specimen revealed endogenous AID expression was detectable in both the hepatocellular carcinoma (HCC) and the non-tumorous inflamed liver tissues in patients with chronic HCV infection while no AID expression was detected in normal liver tissues without HCV infection. In mice, as endogenous AID expression level of liver is not sufficient to detect by immustaining, we utilized the transgenic mouse (Aicda-cre transgenic mouse), which was generated by inserting the coding sequence of Cre recombinase into the locus of the Aicda gene, crossing with a genetically-marked reporter mouse (Rosa-tdRFP mouse). In this Aicda-cre/Rosa-tdRFP mouse model, AID-expressing cells can be detected as RFP-positive. We injected HCV core antigen-inserted plasmid to the liver of Aicda-cre/Rosa-tdRFP mice by “hydrodynamic tail vein injection” procedure. 14days after injection, RFP-positive hepatocytes were detected in the HCV antibody-labelled liver cells whereas no RFP-positive cells were found in control mice with physiological saline injection, indicating HCV core antigen induced aberrant AID expression in hepatocytes. Next, to examine the effect of constitutive AID expression in hepatic-lineage cells, we established a transplantation model of hepatic progenitor cells derived from AID Tg mice. We enriched hepatic progenitor cells from the E13.5 fetal liver of AID Tg mice and transplanted those cells into recipient mice. HCC developed in 7 of 11 (63.6%) recipient mice receiving enriched hepatic progenitor cells from AID Tg mice. On the other hand, no tumor was found in recipient mice transplanted with hepatic progenitor cells derived from wild-type mice (control). In addition, whole exome sequencing of tumor cells clarified the landscape of the accumulated genetic alterations during tumorigenesis and detected a large number of somatic mutations in each tumor tissue originating from the transplanted hepatic progenitor cells of AID Tg mice. These results suggested aberrant AID production induced by HCV infection could contribute to hepatocarcinogenesis via the accumulation of genetic alterations by its mutagenic activity. Citation Format: Soo Ki Kim, Hiroyuki Marusawa, Tsutomu Chiba. Aberrant activation-induced cytidine deaminase (AID) expression is induced in the process of hepatitis C virus-associated hepatocarcinogenesis. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5350. doi:10.1158/1538-7445.AM2014-5350

A systematic profile of clinical inhibitors responsive to EGFR somatic amino acid mutations in lung cancer: implication for the molecular mechanism of drug resistance and sensitivity

Human epidermal growth factor receptor (EGFR) has become a well-established target for the treatment of patients with non-small cell lung cancer (NSCLC). However, a large number of somatic mutations in such protein have been observed to cause drug resistance or sensitivity during pathological progression, limiting the application of reversible EGFR tyrosine kinase inhibitor therapy in NSCLC. In the current work, we describe an integration of in silico analysis and in vitro assay to profile six representative EGFR inhibitors against a panel of 71 observed somatic mutations in EGFR tyrosine kinase domain. In the procedure, the changes in interaction free energy of inhibitors with EGFR upon various mutations were calculated one by one using a rigorous computational scheme, which was preoptimized based on a set of structure-solved, affinity-known samples to improve its performance in characterizing the EGFR-inhibitor system. This method was later demonstrated to be effective in inferring drug response to the classical L858R and G719S mutations that confer constitutive activation for the EGFR kinase. It is found that the Staurosporine, a natural product isolated from the bacterium Streptomyces staurosporeus, exhibits selective inhibitory activity on the T790M and T790M/L858R mutants. This finding was subsequently solidified by in vitro kinase assay experiment; the inhibitory IC50 values of Staurosporine against wild-type, T790M and T790M/L858R mutant EGFR were measured to be 937, 12 and 3 nM, respectively.

Mutational signature of aristolochic acid exposure as revealed by whole-exome sequencing.

In humans, exposure to aristolochic acid (AA) is associated with urothelial carcinoma of the upper urinary tract (UTUC). Exome sequencing of UTUCs from 19 individuals with documented exposure to AA revealed a remarkably large number of somatic mutations and an unusual mutational signature attributable to AA. Most of the mutations (72%) in these tumors were A:T-to-T:A transversions, located predominantly on the nontranscribed strand, with a strong preference for deoxyadenosine in a consensus sequence (T/CAG). This trinucleotide motif overlaps the canonical splice acceptor site, possibly accounting for the excess of splice site mutations observed in these tumors. The AA mutational fingerprint was found frequently in oncogenes and tumor suppressor genes in AA-associated UTUC. The AA mutational signature was observed in one patient's tumor from a UTUC cohort without previous indication of AA exposure. Together, these results directly link an established environmental mutagen to cancer through genome-wide sequencing and highlight its power to reveal individual exposure to carcinogens.

Abstract 36: The OncoNetwork Consortium: A European Collaborative Research study on the development of an Ion AmpliSeq gene panel targeting hot spots in colon and lung cancers.

Abstract The knowledge of molecular alterations involved in colorectal carcinoma (CRC) and non-small-cell lung carcinoma (NSCLC) has significantly increased in the past few years. Molecular subgroups of tumors carrying mutations in druggable genes have been identified in CRC and NSCLC. In this respect, the availability of approaches that are compatible with small amounts of input DNA from FFPE tissue and that allow a rapid screening of a large number of somatic mutations is definitely needed. Ion AmpliSeq™ technology introduces a groundbreaking workflow enabling the rapid sequencing of hundreds of mutations with low allele frequency using as little as 10ng of DNA per reaction. The OncoNetwork Consortium is a European collaborative effort of 8 Cancer Translational Research Institutions to evaluate a custom panel targeting hotspots mutations in 22 genes implicated in CRC and NSCLC. In particular, 7 labs will test the gene panel in 3 phases using 180 retrospective FFPE samples previously tested with orthogonal technologies, and 1 lab will perform confirmatory testing on novel SNPs. According to the established workflow, 10 ng of genomic DNA are used with a single primer pool that includes 86 pairs of primers spanning 15.7kb of targets. Individual libraries are barcoded during the ligation step, quantified and pooled 5 by 5 for amplification on spheres using the Ion OneTouch™ 200 Template Kit v2. Spheres are loaded on Ion 316 chip for sequencing with the Ion PGM™ Sequencer using the Ion PGM™ 200 Sequencing Kit. Results are analyzed using the Ion Reporter™ Software v1.2. In the first phase of the project, 7 labs received the same 5 blind, control FFPE samples. The average coverage obtained for the genes included in the panel ranged between 1.000x and 6.000x. The percentage of reads on target was 88-93% for all samples except for a sample showing a mean of 77%, with an average per base accuracy &amp;gt;98% for all the samples. Data analysis confirmed that all laboratories identified the expected mutations, when using a mapping quality value (MQV) of 4 as first filter. In the second phase of the project, each lab submitted to the consortium 10 previously tested CRC or NSCLC FFPE samples and received back 10 blind samples for sequencing. These samples were expected to contain 5 to 10 variants including point mutations and indels. The six labs that completed this phase were able to detect all the expected variants. Therefore, a 100% true positive rate was found in controls and blind study. The third phase of the project is ongoing and involves the assessment in each participating lab of 10 NSCLC and 5 CRC samples with a low percentage of tumor cells. In conclusion, the preliminary results suggest that the Onconetwork Ion AmpliSeq Colon and Lung Cancer panel is able to detect with high sensitivity and reproducibility the most frequent mutations of CRC and NSCLC by using a just 10 ng of DNA from FFPE samples. Citation Format: Nicola Normanno, Rosella Petraroli, Alain Rico, Henriette Kurth, Bastiaan Tops, Eliana Amato, Andrea Mafficini, Delphine Le Corre, Anna Maria Rachiglio, Anne Reiman, Orla Sheils, Christoph Noppen, Nora Rieber, Jonathan Mangion, Chrysanthi Ainali, Ludovic Lacroix, Ian Cree, Marjolijn Ligtenberg, Andy Felton, Aldo Scarpa, Pierre Laurent Puig. The OncoNetwork Consortium: A European Collaborative Research study on the development of an Ion AmpliSeq gene panel targeting hot spots in colon and lung cancers. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 36. doi:10.1158/1538-7445.AM2013-36

An empirical Bayesian framework for somatic mutation detection from cancer genome sequencing data

Recent advances in high-throughput sequencing technologies have enabled a comprehensive dissection of the cancer genome clarifying a large number of somatic mutations in a wide variety of cancer types. A number of methods have been proposed for mutation calling based on a large amount of sequencing data, which is accomplished in most cases by statistically evaluating the difference in the observed allele frequencies of possible single nucleotide variants between tumours and paired normal samples. However, an accurate detection of mutations remains a challenge under low sequencing depths or tumour contents. To overcome this problem, we propose a novel method, Empirical Bayesian mutation Calling (https://github.com/friend1ws/EBCall), for detecting somatic mutations. Unlike previous methods, the proposed method discriminates somatic mutations from sequencing errors based on an empirical Bayesian framework, where the model parameters are estimated using sequencing data from multiple non-paired normal samples. Using 13 whole-exome sequencing data with 87.5–206.3 mean sequencing depths, we demonstrate that our method not only outperforms several existing methods in the calling of mutations with moderate allele frequencies but also enables accurate calling of mutations with low allele frequencies (≤10%) harboured within a minor tumour subpopulation, thus allowing for the deciphering of fine substructures within a tumour specimen.

MS1-3: Tumor Heterogeneity and Metastasis.

Abstract Massively parallel sequencing using Next Gen sequencing approaches have opened a new window into our understanding of the molecular architecture of breast cancer. A significant advantage of these technologies are that they are digital in nature because the measurement technology begins with single DNA molecules attached to a substrate. As a result, not only do these technologies provide an efficient means to identify somatic mutations, information is also provided on the frequency with which any given mutation is observed within the population of DNA molecules in the sample under analysis. The first and most striking feature of Next Gen analysis of breast cancer is the large numbers of somatic mutations that have been identified. Genome wide, thousands of mutations can be present in a single tumor. While this finding generates immense complexity from the perspective of defining biological drivers versus passengers, it does open up a new opportunity to study variations in mutation frequency. In a perfectly monoclonal tumor with a diploid genome and no contaminating normal cell DNA the mutation frequencies generated by digital sequencing will be 50% for a heterozygous mutation or 100% for homozygous mutations. Very few breast cancers subjected to whole genome sequencing exhibit this pattern, rather mutations occur with a wide variation of frequencies from the detection limit, around 5%, through to 100%. Furthermore mutation frequency clustering often occurs, a phenomenon best explained by the presence of a repertoire of mutations present in a founder clone that are present in all tumors, combined with additional less common mutations that represent a subdominant or minority population of cells that have arisen through additional clonal outgrowth. In a multiclonal model one of the most critical questions is “which clone determines the prognosis?”. One way to begin to answer this question is to sequence tumors before and after therapy and to compare primary tumors versus metastases. In our initial experience of a patient with a basal-like breast cancer in which a whole cancer genome was generated from her primary, her subsequent brain metastasis and a mouse xenograft generated from the breast primary, there was clear evidence of a shift in mutation frequency primary to metastasis, suggesting the brain metastasis arose from a minor sub-clone. The mutational profile of the xenograft more closely resembled the metastasis rather than the primary, suggesting the grafting process captured the metastatic clone (1). In this invited lecture we will present an extension of our sequencing experiences with further human progenitor-mouse xenograft comparisons, comparisons between human primaries and paired metastases and also examples of ER+ breast cancer genomes obtained before and after neoadjuvant endocrine therapy. In all these experiments multiclonality appears to be the rule rather than the exception. The clinical and biological implications of these findings will be discussed. 1. Ding L, Ellis MJ, Li S, Larson DE, Chen K, Wallis JW, et al. Genome remodelling in a basal-like breast cancer metastasis and xenograft. Nature. 2010;464(7291):999–1005. PMCID: 2872544. Citation Information: Cancer Res 2011;71(24 Suppl):Abstract nr MS1-3.

ProKinO: An Ontology for Integrative Analysis of Protein Kinases in Cancer

BackgroundProtein kinases are a large and diverse family of enzymes that are genomically altered in many human cancers. Targeted cancer genome sequencing efforts have unveiled the mutational profiles of protein kinase genes from many different cancer types. While mutational data on protein kinases is currently catalogued in various databases, integration of mutation data with other forms of data on protein kinases such as sequence, structure, function and pathway is necessary to identify and characterize key cancer causing mutations. Integrative analysis of protein kinase data, however, is a challenge because of the disparate nature of protein kinase data sources and data formats.ResultsHere, we describe ProKinO, a protein kinase-specific ontology, which provides a controlled vocabulary of terms, their hierarchy, and relationships unifying sequence, structure, function, mutation and pathway information on protein kinases. The conceptual representation of such diverse forms of information in one place not only allows rapid discovery of significant information related to a specific protein kinase, but also enables large-scale integrative analysis of protein kinase data in ways not possible through other kinase-specific resources. We have performed several integrative analyses of ProKinO data and, as an example, found that a large number of somatic mutations (∼288 distinct mutations) associated with the haematopoietic neoplasm cancer type map to only 8 kinases in the human kinome. This is in contrast to glioma, where the mutations are spread over 82 distinct kinases. We also provide examples of how ontology-based data analysis can be used to generate testable hypotheses regarding cancer mutations.ConclusionWe present an integrated framework for large-scale integrative analysis of protein kinase data. Navigation and analysis of ontology data can be performed using the ontology browser available at: http://vulcan.cs.uga.edu/prokino.

The structural impact of cancer-associated missense mutations in oncogenes and tumor suppressors

BackgroundCurrent large-scale cancer sequencing projects have identified large numbers of somatic mutations covering an increasing number of different cancer tissues and patients. However, the characterization of these mutations at the structural and functional level remains a challenge.ResultsWe present results from an analysis of the structural impact of frequent missense cancer mutations using an automated method. We find that inactivation of tumor suppressors in cancer correlates frequently with destabilizing mutations preferably in the core of the protein, while enhanced activity of oncogenes is often linked to specific mutations at functional sites. Furthermore, our results show that this alteration of oncogenic activity is often associated with mutations at ATP or GTP binding sites.ConclusionsWith our findings we can confirm and statistically validate the hypotheses for the gain-of-function and loss-of-function mechanisms of oncogenes and tumor suppressors, respectively. We show that the distinct mutational patterns can potentially be used to pre-classify newly identified cancer-associated genes with yet unknown function.

Important role of indels in somatic mutations of human cancer genes.

BackgroundCancer is clonal proliferation that arises owing to mutations in a subset of genes that confer growth advantage. More and more cancer related genes are found to have accumulated somatic mutations. However, little has been reported about mutational patterns of insertions/deletions (indels) in these genes.ResultsWe analyzed indels' abundance and distribution, the relative ratio between indels and somatic base substitutions and the association between those two forms of mutations in a large number of somatic mutations in the Catalogue of Somatic Mutations in Cancer database. We found a strong correlation between indels and base substitutions in cancer-related genes and showed that they tend to concentrate at the same locus in the coding sequences within the same samples. More importantly, a much higher proportion of indels were observed in somatic mutations, as compared to meiotic ones. Furthermore, our analysis demonstrated a great diversity of indels at some loci of cancer-related genes. Particularly in the genes with abundant mutations, the proportion of 3n indels in oncogenes is 7.9 times higher than that in tumor suppressor genes.ConclusionsThere are three distinct patterns of indel distribution in somatic mutations: high proportion, great abundance and non-random distribution. Because of the great influence of indels on gene function (e.g., the effect of frameshift mutation), these patterns indicate that indels are frequently under positive selection and can often be the 'driver mutations' in oncogenesis. Such driver forces can better explain why much less frameshift mutations are in oncogenes while much more in tumor suppressor genes, because of their different function in oncogenesis. These findings contribute to our understanding of mutational patterns and the relationship between indels and cancer.

Prediction of Cancer Driver Mutations in Protein Kinases

A large number of somatic mutations accumulate during the process of tumorigenesis. A subset of these mutations contribute to tumor progression (known as "driver" mutations) whereas the majority of these mutations are effectively neutral (known as "passenger" mutations). The ability to differentiate between drivers and passengers will be critical to the success of upcoming large-scale cancer DNA resequencing projects. Here we show a method capable of discriminating between drivers and passengers in the most frequently cancer-associated protein family, protein kinases. We apply this method to multiple cancer data sets, validating its accuracy by showing that it is capable of identifying known drivers, has excellent agreement with previous statistical estimates of the frequency of drivers, and provides strong evidence that predicted drivers are under positive selection by various sequence and structural analyses. Furthermore, we identify particular positions in protein kinases that seem to play a role in oncogenesis. Finally, we provide a ranked list of candidate driver mutations.

A Hypermutation Phenotype and Somatic MSH6 Mutations in Recurrent Human Malignant Gliomas after Alkylator Chemotherapy

Malignant gliomas have a very poor prognosis. The current standard of care for these cancers consists of extended adjuvant treatment with the alkylating agent temozolomide after surgical resection and radiotherapy. Although a statistically significant increase in survival has been reported with this regimen, nearly all gliomas recur and become insensitive to further treatment with this class of agents. We sequenced 500 kb of genomic DNA corresponding to the kinase domains of 518 protein kinases in each of nine gliomas. Large numbers of somatic mutations were observed in two gliomas recurrent after alkylating agent treatment. The pattern of mutations in these cases showed strong similarity to that induced by alkylating agents in experimental systems. Further investigation revealed inactivating somatic mutations of the mismatch repair gene MSH6 in each case. We propose that inactivating somatic mutations of MSH6 confer resistance to alkylating agents in gliomas in vivo and concurrently unleash accelerated mutagenesis in resistant clones as a consequence of continued exposure to alkylating agents in the presence of defective mismatch repair. The evidence therefore suggests that when MSH6 is inactivated in gliomas, alkylating agents convert from induction of tumor cell death to promotion of neoplastic progression. These observations highlight the potential of large scale sequencing for revealing and elucidating mutagenic processes operative in individual human cancers.

Tracking of the V4-34 (VH4-21) gene in human tonsil reveals clonal isotype switch events and a highly variable degree of somatic hypermutation.

The V4-34 (VH4-21) gene has been found to encode certain IgM autoantibodies, and is mandatory for pathological IgM anti-erythrocyte antibodies of I/i specificity. The gene is also commonly used by normal IgM-positive B lymphocytes, but its involvement in B cells which have undergone class switching to IgG or IgA is less clear. In order to track V4-34 gene usage and class switching events during a normal immune response, we have probed RNA in a limited area of human tonsil. Results indicate that the V4-34 gene undergoes class switching to IgG or IgA, with the progeny either remaining unmutated or containing large numbers of somatic mutations. Mutational patterns indicate possible 'hot spots', and some mutations appear deleterious. At the level of individual B cells, we have tracked a clonal isotype switch event from IgM to IgA, with each retaining close to germ-line configuration. In addition, we have followed a clonal switch from a mutated IgM to IgG, with no further accumulation of somatic mutations. These data indicate that the V4-34 gene is involved in a maturing immune response, and that the routes to production of IgG or IgA antibodies are various.

Cytosine methylation determines hot spots of DNA damage in the human P53 gene.

In the P53 tumor suppressor gene, a remarkably large number of somatic mutations are found at methylated CpG dinucleotides. We have previously mapped the distribution of (+/-) anti-7beta,8alpha-dihydroxy-9alpha,10alpha-epoxy -7,8,9,10-tetrahydrobenzo[a]pyrene (BPDE) adducts along the human P53 gene [Denissenko, M. F., Pao, A., Tang, M.-s. & Pfeifer, G. P. (1996) Science 274, 430-432]. Strong and selective formation of adducts occurred at guanines in CpG sequences of codons 157, 248, and 273, which are the major mutational hot spots in lung cancer. Chromatin structure was not involved in preferential modification of these sites by BPDE. To investigate other possible mechanisms underlying the selectivity of BPDE binding, we have mapped the adducts in plasmid DNA containing genomic P53 sequences. The adduct profile obtained was different from that in genomic DNA. However, when cytosines at CpG sequences were converted to 5-methylcytosines by the CpG-specific methylase SssI and the DNA was subsequently treated with BPDE, adduct hot spots were created which were similar to those seen in genomic DNA where all CpGs are methylated. A strong positive effect of 5-methylcytosine on BPDE adduct formation at CpG sites was also documented with sequences of the PGK1 gene derived from an active or inactive human X chromosome and having differential methylation patterns. These results show that methylated CpG dinucleotides, in addition to being an endogenous promutagenic factor, may represent a preferential target for exogenous chemical carcinogens. The data open new avenues concerning the reasons that the majority of mutational hot spots in human genes are at CpGs.