| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Priority Reports |
1 Cancer Genome Project, Wellcome Trust Sanger Institute, Hinxton, United Kingdom; 2 Molecular Pathology Unit, Brain Tumor Center, Neurosurgical Service and Center for Cancer Research, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts; 3 Cancer Research UK Genetic Epidemiology Unit, University of Cambridge, Cambridge, United Kingdom; 4 Johns Hopkins University School of Medicine, Baltimore, Maryland; and 5 Institute of Cancer Research, Sutton, Surrey, United Kingdom
Requests for reprints: P. Andrew Futreal, Cancer Genome Project, Wellcome Trust Sanger Institute, Hinxton CB10 1SA, United Kingdom. Phone: 44-1223-494857; Fax: 44-1223-494809; E-mail: paf{at}sanger.ac.uk.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
| Materials and Methods |
|---|
|
|
|---|
554 kb, of which 532 kb (98%) was successfully covered. Amplification of sequencing templates from the tumors, bidirectional sequencing, and confirmations were done as previously described (3). Greater than 90% of designed exons were successfully analyzed in >90% of samples. DNA from pretreatment tumor material within the H&E formalin-fixed paraffin-embedded slide from case PD1487a was prepared using the Gentra Puregene DNA Purification kit and sequenced using USB P-33 Thermosequenase kit. MGMT methylation analysis was done as detailed elsewhere (6, 7). In addition to the protein kinase domains, TP53 and PTEN were sequenced in the set (Supplementary Data). Genome-wide analysis of copy number was done on paired normal/tumor samples as previously described (8). To assess the presence of biological selection on the observed somatic mutations, analyses of the ratio of nonsynonymous to synonymous mutations were done as previously described (4). Additional statistical methods applied can be found in Supplementary Data. | Results |
|---|
|
|
|---|
The 48 nonsynonymous and splice site mutations were distributed among 44 genes. Three nonsynonymous mutations were found in BRD2 (one in PD1487a and two in PD1489a) and two were found in TRIM33 (in PD1487a and PD1489a). Two nonsynonymous mutations were also found in ERN1 (both in PD1487a). The remainder of the genes each had one or zero nonsynonymous mutations. Mutations in FGFR1, recently reported in a screen of 20 kinase domains (9), were not observed.
No frequently mutated protein kinase was identified in this screen. Nevertheless, there were individual mutations of potential interest. These include a mutation in the G-loop that mediates ATP binding (MLCK), mutations within the activation segment (MAPK8, PRKCB1, PTK2, and TLK1), a mutation of the canonical proline within the PSTAIRE motif of CDK2 important in interactions with cyclins/CDK inhibitor proteins (10), and a mutation in PDGFRA, amplification of which has previously been associated with glioma (11). To evaluate the significance of mutations detected, a follow-up screen was conducted of 21 glioblastomas through the kinase domain exons of 40 genes in which mutations had been found. Only two silent mutations were detected in the follow-up screen (Supplementary Data).
Six gliomas had no somatic mutations and one had a single mutation. In contrast to this generally low prevalence, two gliomas had numerous mutations, PD1487a with 34 and PD1489a with 30. Both PD1487a and PD1489a are early cultures derived from recurrences of glioblastomas. Analysis of a subset of these mutations in archival formalin-fixed, paraffin-embedded tissues showed that the mutations were present in the recurrences from which the cultures were derived (Supplementary Data).
The mutational patterns in PD1487a and PD1489a are shown in Tables 1 and 2 . Because they exhibit similar patterns, the data from these two tumors have been combined and their features are summarized together. All the somatic mutations detected were single base substitutions. In both PD1487a and PD1489a, most substitutions were C:G>T:A transitions (63 of 64, 98%). These C:G>T:A mutations occurred in a particular sequence context (P = 0.0006; Table 2). Considering the strand on which each C:G>T:A mutation is a C>T change, there was an excess of mutations at CpC dinucleotides where the 5' cytosine is mutated (50%) compared with that expected by chance (28%). There was also a reduction in guanine immediately 5' to the mutated cytosine (11% compared with 25%) and a weak strand bias for mutation of guanine on the nontranscribed strand (61% G>A to 39% C>T).
|
|
To identify the underlying causes of the mutator phenotype, several known DNA repair genes were evaluated. Inactivation of MGMT (O-6-methylguanine-DNA methyltransferase) by promoter methylation has been implicated in the response of gliomas to alkylating agent chemotherapy (7). MGMT status was evaluated by both sequencing and promoter methylation analysis. No mutations were observed in the tumors under study. Methylation of the MGMT promoter was observed in both the primary tumor recurrence and early culture of PD1489a and in one additional glioblastoma (PD1537), but not in PD1487a. Thus, MGMT status was not concordant with the mutator phenotype.
Somatic mutations were found in the MSH6 gene in both PD1487a and PD1489a. MSH6 heterodimerizes with MSH2 recognizing double-stranded DNA mismatches (12). Homozygous-inactivating germ line alleles of MSH6 have recently been reported in two children with glioma (13, 14) but somatic mutations of MSH6 have not previously been reported in gliomas. In PD1489a, a single base pair deletion, c.2425delG/p.809V>X, resulting in a translational frameshift was present. This protein truncating mutation was homozygous due to complete loss of the second copy of chromosome 2 (Supplementary Data). A somatic heterozygous nonsense mutation, c1453C>T/p.485Q>X, was detected in PD1487a. Loss of heterozygosity on chromosome 2 was not observed in this sample. However, a second somatic missense mutation, c.3907G>A/p.1303A>T, was detected. Alanine 1303 is conserved in nearly all MutS orthologues through prokaryotes and is within the ATPase domain (Supplementary Data), critical for the mismatch repair function of this family of proteins. Therefore, p.1303A>T is likely to abrogate MSH6 function. These somatic mutations were confirmed in DNA from the recurrent gliomas from which the early cultures were derived. Evaluation of the NIH-National Cancer Institute consensus markers and an additional seven dinucleotide repeats revealed no evidence of microsatellite instability in PD1487a, PD1489a, or in any other glioma (Supplementary Data).
A single H&E-stained section of homogeneous, densely cellular tumor was available from the pretreatment glioblastoma specimen of PD1487a. DNA was extracted from this section and PCR amplified to investigate the presence of the two MSH6 mutations found in the recurrent glioma. Neither of the two mutations was detected (Fig. 1 ). This result suggests that the MSH6 mutations were selected as a subclone of the primary glioma that subsequently expanded to constitute most of the posttreatment recurrence. No pretreatment material was available for PD1489a. MSH6 mutations were not found in the six primary malignant gliomas that had not previously been subject to therapy or in the glioma recurrent after radiotherapy alone (PD1536a). Screening of an additional series of 21 pretreatment malignant gliomas also failed to reveal any further MSH6 mutations.
|
|
| Discussion |
|---|
|
|
|---|
Exposure of a population of cells to a mutagen will result in a different set of mutations in each cell within the population. Such scattered mutations would not be detectable by sequencing bulk DNA made from the whole cell population. The fact that we can detect these somatic mutations therefore indicates that expansion of a clone of tumor cells underlies the emergence of these recurrences after alkylating agent chemotherapy. The observation of MSH6 mutations in the PD1487a recurrence but not in the pretreatment glioblastoma sample provides direct evidence in favor of this hypothesis. Therefore, the results show for the first time that, at least in some cases, recurrence after temozolomide therapy in patients with glioblastoma is through clonal evolution.
Alkylating agents induce adducts on DNA, notably O6-methylguanine. O6-methylguanine can mispair with thymine in double-stranded DNA (12). This mismatch is recognized by the MSH2/MSH6 dimer component of the mismatch repair system. Overwhelming numbers of adducts and G/T mismatches cause cells to enter apoptosis. This cell killing effect is dependent on the integrity of the mismatch repair system and specifically of MSH2/MSH6 dimers. Several studies in experimental models indicate that if the mismatch recognition mechanism is not functional, cells become resistant to the killing effects of alkylating agents (19, 20). Malignant gliomas that recur following alkylator chemotherapy typically do not respond to further alkylating agents (2). Given the observed MSH6 mutations in recurrent gliomas, we propose that somatic MSH6 mutations confer resistance to alkylating agents on clones of tumor cells that subsequently expand and manifest as clinical recurrence.
The etiology of the somatic MSH6 mutations is unclear. They could be due to endogenous mutagenic processes or could have themselves been the result of mutagenesis due to alkylating agents. In PD1489, the single base pair deletion in MSH6 together with loss of the other chromosome 2 copy are not typical of alkylating agent exposure. Conversely, both mutations in PD1487 were C:G>T:A changes, neither of which were at CpG dinucleotides, a pattern more compatible with the mutagenic effects of alkylating agents.
Mutational inactivation of MSH6 and the mismatch repair system seems to have had further consequences. In experimental models, cells with mismatch repair deficiency exposed to alkylating agents have an elevated mutation rate because G/T mismatches are no longer recognized and repaired (12). PD1487a and PD1489a exhibit the highest prevalence of somatic mutations per nucleotide screened of any primary cancers or cancer cell lines we have examined thus far by extensive sequencing (35). On the assumption that most of the somatic mutations observed are passenger mutations and hence that the prevalence of mutations is similar in the genome as a whole, these two cancers each have at least 200,000 somatic point mutations in their genomes,
1 per 15 kb. If distributed randomly, these would be expected to generate
1,500 amino acid changes in each cancer, affecting
8% protein coding genes. Thus, in addition to conferring resistance, inactivating mutations of MSH6 in the context of alkylating agent exposure foster a hypermutational process in gliomas, which is very likely to facilitate rapid evolution of clones with growth advantage and contribute to progression.
It has been shown in cell culture systems that alkylator chemoresistance emerges from mismatch repair deficiency (21). Our findings represent the first demonstration of this paradigm in vivo. They also reveal the development of rampant mutagenesis as a consequence of continuing alkylator exposure that is likely contributory to disease progression. Given the inherent difficulties in obtaining posttreatment brain tumor samples for molecular evaluation, a coordinated multicenter study is now indicated to explore in greater detail the role of MSH6 and potentially other mismatch repair defects in recurrent, chemoresistant malignant gliomas.
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank the patients who participated in this research, Wendy Haynes for administrative support, Nazneen Rahman for comments on the manuscript, and Darrell Bigner (Department of Pathology, Duke University Medical Center, Durham, NC) for the D-542MG cell line.
| Footnotes |
|---|
C. Hunter, R. Smith, and D.P. Cahill contributed equally to this work.
6 http://www.cbtrus.org/; http://seer.cancer.gov/. ![]()
7 http://www.sanger.ac.uk/genetics/CGP/Kinases/. ![]()
8 http://www.sanger.ac.uk/cosmic. ![]()
9 http://www-p53.iarc.fr/index.html. ![]()
Received 1/13/06. Revised 2/13/06. Accepted 2/24/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. W. Parsons, S. Jones, X. Zhang, J. C.-H. Lin, R. J. Leary, P. Angenendt, P. Mankoo, H. Carter, I-M. Siu, G. L. Gallia, et al. An Integrated Genomic Analysis of Human Glioblastoma Multiforme Science, September 26, 2008; 321(5897): 1807 - 1812. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Maxwell, S. P. Johnson, R. E. McLendon, D. W. Lister, K. S. Horne, A. Rasheed, J. A. Quinn, F. Ali-Osman, A. H. Friedman, P. L. Modrich, et al. Mismatch Repair Deficiency Does Not Mediate Clinical Resistance to Temozolomide in Malignant Glioma Clin. Cancer Res., August 1, 2008; 14(15): 4859 - 4868. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. M. Loriaux, R. L. Levine, J. W. Tyner, S. Frohling, C. Scholl, E. P. Stoffregen, G. Wernig, H. Erickson, C. A. Eide, R. Berger, et al. High-throughput sequence analysis of the tyrosine kinome in acute myeloid leukemia Blood, May 1, 2008; 111(9): 4788 - 4796. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Verbeek, T. D. Southgate, D. E. Gilham, and G. P. Margison O6-Methylguanine-DNA methyltransferase inactivation and chemotherapy Br. Med. Bull., March 1, 2008; 85(1): 17 - 33. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. H Scott, T. Homfray, N. L Huxter, S. G Mitton, R. Nash, M. N Potter, D. Lancaster, and N. Rahman Familial T-cell non-Hodgkin lymphoma caused by biallelic MSH2 mutations J. Med. Genet., July 1, 2007; 44(7): e83 - e83. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. P. Cahill, K. K. Levine, R. A. Betensky, P. J. Codd, C. A. Romany, L. B. Reavie, T. T. Batchelor, P. A. Futreal, M. R. Stratton, W. T. Curry, et al. Loss of the Mismatch Repair Protein MSH6 in Human Glioblastomas Is Associated with Tumor Progression during Temozolomide Treatment Clin. Cancer Res., April 1, 2007; 13(7): 2038 - 2045. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. K. Mirzoeva, T. Kawaguchi, and R. O. Pieper The Mre11/Rad50/Nbs1 complex interacts with the mismatch repair system and contributes to temozolomide-induced G2 arrest and cytotoxicity. Mol. Cancer Ther., November 1, 2006; 5(11): 2757 - 2766. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |