Skip to main content
Advertisement
  • Loading metrics

Whole Exome Sequencing Identifies TSC1/TSC2 Biallelic Loss as the Primary and Sufficient Driver Event for Renal Angiomyolipoma Development

  • Krinio Giannikou ,

    Contributed equally to this work with: Krinio Giannikou, Izabela A. Malinowska, Trevor J. Pugh

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Izabela A. Malinowska ,

    Contributed equally to this work with: Krinio Giannikou, Izabela A. Malinowska, Trevor J. Pugh

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Trevor J. Pugh ,

    Contributed equally to this work with: Krinio Giannikou, Izabela A. Malinowska, Trevor J. Pugh

    Affiliation Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • Rachel Yan,

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Yuen-Yi Tseng,

    Affiliation Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • Coyin Oh,

    Affiliation Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • Jaegil Kim,

    Affiliation Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • Magdalena E. Tyburczy,

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Yvonne Chekaluk,

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Yang Liu,

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Nicola Alesi,

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Geraldine A. Finlay,

    Affiliation Tufts New England Medical Center, Boston, Massachusetts, United States of America

  • Chin-Lee Wu,

    Affiliation Massachusetts General Hospital, Boston, Massachusetts, United States of America

  • Sabina Signoretti,

    Affiliation Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America

  • Matthew Meyerson,

    Affiliation Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • Gad Getz,

    Affiliation Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • Jesse S. Boehm,

    Affiliation Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • Elizabeth P. Henske ,

    EHENSKE@BWH.HARVARD.EDU (EPH); dk@rics.bwh.harvard.edu (DJK)

    Affiliations Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America, Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  •  [ ... ],
  • David J. Kwiatkowski

    EHENSKE@BWH.HARVARD.EDU (EPH); dk@rics.bwh.harvard.edu (DJK)

    Affiliations Division of Pulmonary and Critical Care Medicine and of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America, Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America

  • [ view all ]
  • [ view less ]

Abstract

Renal angiomyolipoma is a kidney tumor in the perivascular epithelioid (PEComa) family that is common in patients with Tuberous Sclerosis Complex (TSC) and Lymphangioleiomyomatosis (LAM) but occurs rarely sporadically. Though histologically benign, renal angiomyolipoma can cause life-threatening hemorrhage and kidney failure. Both angiomyolipoma and LAM have mutations in TSC2 or TSC1. However, the frequency and contribution of other somatic events in tumor development is unknown. We performed whole exome sequencing in 32 resected tumor samples (n = 30 angiomyolipoma, n = 2 LAM) from 15 subjects, including three with TSC. Two germline and 22 somatic inactivating mutations in TSC2 were identified, and one germline TSC1 mutation. Twenty of 32 (62%) samples showed copy neutral LOH (CN-LOH) in TSC2 or TSC1 with at least 8 different LOH regions, and 30 of 32 (94%) had biallelic loss of either TSC2 or TSC1. Whole exome sequencing identified a median of 4 somatic non-synonymous coding region mutations (other than in TSC2/TSC1), a mutation rate lower than nearly all other cancer types. Three genes with mutations were known cancer associated genes (BAP1, ARHGAP35 and SPEN), but they were mutated in a single sample each, and were missense variants with uncertain functional effects. Analysis of sixteen angiomyolipomas from a TSC subject showed both second hit point mutations and CN-LOH in TSC2, many of which were distinct, indicating that they were of independent clonal origin. However, three tumors had two shared mutations in addition to private somatic mutations, suggesting a branching evolutionary pattern of tumor development following initiating loss of TSC2. Our results indicate that TSC2 and less commonly TSC1 alterations are the primary essential driver event in angiomyolipoma/LAM, whereas other somatic mutations are rare and likely do not contribute to tumor development.

Author Summary

We performed comprehensive genome analysis of a kidney tumor called angiomyolipoma. These tumors are known to develop in most individuals who have Tuberous Sclerosis Complex (TSC) and those who have sporadic lymphangioleiomyomatosis (LAM), and are seen rarely in the general population. In these angiomyolipomas, we found consistent involvement of the TSC2 and TSC1 genes that are known to cause TSC, but very few (<5 on average) mutations elsewhere in the protein-coding regions. This is in stark contrast to other adult solid tumours that typically harbor hundreds to thousands of such mutations. Our results indicate that genetic alterations in TSC2/TSC1 are the primary and essential driver genetic events for development and progression of renal angiomyolipoma. Analysis of multiple angiomyolipomas from a single patient showed distinct genetic aberrations in the majority of samples, indicating that most of the tumors had developed independently. Branched clonal evolution was evident from the observation of three tumors that shared 2 mutations in addition to mutations private to each. Our results indicate that therapeutic approaches for treatment of patients with angiomyolipoma should focus on the consequences of TSC2/TSC1 loss, including but not limited to mTOR activation.

Introduction

Renal angiomyolipoma is a pathologically benign mesenchymal kidney tumor, characterized by vascular, smooth muscle, and adipocyte elements [1, 2]. Angiomyolipomas are rare in the general population (autopsy prevalence about 1 in 1,000 [3]) but are seen in > 70% adults with Tuberous Sclerosis Complex (TSC (MIM: 191100, 613254 [4]) where they are usually multifocal and bilateral [5]. TSC is a genetic disorder characterized by seizures, tumor development in the brain, heart, kidney, and skin, and a distinctive set of neurodevelopmental syndromes known as TSC-associated neurologic disorder (TAND) [6, 7]. TSC is due to inactivating heterozygous or mosaic mutations in either TSC1 (~21%) or TSC2 (~79%) [6]. Lymphangioleiomyomatosis (LAM) is a pathologic condition related to angiomyolipoma in which a similar cellular proliferation occurs in the lung in the form of small (< 5 mm) nodules, causing progressive cystic lung disease that can be fatal [8]. LAM occurs almost exclusively in women, and is seen at much higher frequency in TSC than in the general population, with up to 80% of women with TSC having evidence of cystic lung disease [9]. Both angiomyolipoma and LAM belong to the family of Perivascular Epithelioid Cell tumors (PEComa) [10, 11].

Previous genetic studies have shown that angiomyolipomas arising in TSC patients occur due to biallelic inactivation of either TSC2 or TSC1 [1]. The first mutational event (‘hit’) in TSC2 or TSC1 is the germline mutation that is the cause of TSC in the individual. The genetic “second hit” event leads to hyperactivation of mTORC1 (mammalian Target Of Rapamycin Complex 1) and contributes to tumor development [12, 13]. Sporadic, non-TSC associated, renal angiomyolipoma and LAM are due nearly exclusively to mutations in TSC2 [14, 15]. Concurrent analysis of angiomyolipomas and LAM lung samples from women with sporadic LAM but not TSC has demonstrated that these two lesions have identical point mutations in TSC2, providing strong evidence that they are clonal and derived from a common cell [16]. However, TSC2 mutations have not been identified in all cases of LAM, possibly due to technical limitations related to the small amount of tumor tissue available [14]. Angiomyolipoma samples from patients with either sporadic or TSC-associated tumors often display loss of the wild type allele (loss of heterozygosity, LOH) for TSC2 [17, 18] [19], as well as evidence of activation of mTORC1 [12, 13, 19].

It is unknown whether additional genetic events, beyond loss of TSC1 or TSC2, contribute to angiomyolipoma development, particularly in large tumors requiring surgical resection. We report here the first exome-wide genetic analysis of these tumors.

Results

Samples

Thirty-two tumor samples (n = 30 renal angiomyolipomas; n = 2 LAM) from 15 patients were examined for somatic mutations by performing whole exome sequencing (Table 1, Fig 1). For one subject with LAM but not TSC (P1), we analyzed both a resected abdominal LAM tumor and a LAM cell cluster isolated from chylous pleural fluid. Nineteen angiomyolipoma samples were available from 3 TSC subjects, including one subject undergoing bilateral kidney removal (P13), from whom 16 macroscopically-distinct renal angiomyolipomas (S14-S29) were available, with 12 tumors from the left kidney and 4 from the right (S1 Fig). Single angiomyolipoma samples were collected from 11 subjects without TSC (P2, P3, P5-P12, P15), of whom two had LAM. All samples were operative samples reflecting the clinical situation that these tumors were relatively large necessitating surgical intervention.

thumbnail
Table 1. LIst of patient samples available, and mutation findings in TSC2 and TSC1.

https://doi.org/10.1371/journal.pgen.1006242.t001

thumbnail
Fig 1. The diagram shows the different kinds of tumor samples analyzed from the different kinds of patients.

S#: Sample number, P#: Patient number.

https://doi.org/10.1371/journal.pgen.1006242.g001

Mutation and LOH findings in TSC2 and TSC1

Inactivating mutations in TSC2 were found in 30 of 32 (94%) tumors from 14 of 15 patients (Table 1). The remaining 2 tumors (6%, S30 and S31) were from one patient (P14), who had a germ line TSC1 nonsense mutation (c.2074C>T; p.R692*). Among the 22 somatic mutations in TSC2, there were 6 small deletions (4 single nucleotide, one 18 nucleotide), 5 splice site mutations, 7 nonsense mutations, 2 small insertions (1 nucleotide each), and 2 large-scale deletions that removed the entire gene in one sample (S7) (homozygous deletion region 50kb, single copy loss region 150kb) (Fig 2). Apart from the absence of missense mutations, this mutation distribution is similar to that for TSC2 germ-line mutations [20]. See S1 Table for a glossary of genetic terms used in this publication.

thumbnail
Fig 2. Map of 20 somatic TSC2 mutations detected in angiomyolipoma and LAM specimens.

The GAP domain of TSC2 is indicated. Novel variants (n = 13) are in black font; variants previously reported (n = 5) are in blue font. Circle colors reflect different types of mutation, as indicated. Note that in two instances, the mutations differ by a single amino acid position, and hence their circles overlap (p.R749delfs/p.L750delfs, p.Q1588*/p.P1589delfs). Two mutations were seen twice each in two samples.

https://doi.org/10.1371/journal.pgen.1006242.g002

Of the 11 angiomyolipoma samples from subjects without TSC, six (S6, S8, S9, S10, S11, S12) had two distinct TSC2 somatic mutations, indicative of biallelic inactivation. One sample (S7) showed a homozygous somatic deletion (50kb) of TSC2, as noted above, confirmed by multiplex ligation-dependent probe amplification (MLPA) and CapSeg analysis (S2 Fig) [21]. The remaining 4 of 11 non-TSC angiomyolipoma samples (S3, S4, S13, S32) showed copy neutral LOH. The chylous fluid and abdominal LAM samples (S1, S2) from patient P1 who had sporadic LAM, had the same TSC2 single nucleotide mutation (c.2319delA, p.L773delfs) and copy neutral LOH. Thus all 13 sporadic angiomyolipoma and LAM specimens showed evidence of biallelic TSC2 inactivation.

Analysis of the 16 renal angiomyolipomas from one patient (P13), who has TSC, revealed a non-mosaic germline TSC2 deletion mutation (c.5135delC) that was also seen in normal tissue (Table 1). Two tumor samples (S26, S27) from this subject had TSC2 somatic “second hit” mutations. Twelve of the remaining 14 angiomyolipomas showed TSC2 copy neutral LOH, while neither second hit mutation nor LOH was detected in two of 16 samples (S16, S28) (Table 1).

LOH mapping using informative SNPs for 18 different samples from 6 patients permitted the delineation of at least 8 different regions of chromosome 16p LOH (Fig 3; S3 Fig). In all of these tumors the LOH region extended from the first SNP identified near 16pter to a region more centromeric on chromosome 16p, and encompassed TSC2 at g.16:2,095,990–2,140,713 (hg19). The smallest region of LOH was ~ 3Mb, and the largest was ~ 30Mb. The largest LOH region included the centromeric region of chromosome 16 but none extended to include all of chromosome 16, suggesting that mitotic recombination had occurred somewhere on the p arm or centromere in each case to replace the wild type TSC2 locus with the pathogenic mutation. Note that this mechanism preserves diploid copy number for all other genes on 16p. Copy neutral LOH for all of 9q, including TSC1 at g.9:135,764,735–135,822,020 (hg19) on 9q34, was seen in the two samples (S30, S31) from patient P14 with the germline TSC1 p.R692* mutation. Copy number analysis using ABSOLUTE [22] also identified a subclonal, 2.4-fold cancer cell fraction copy number for chromosome 9q in sample S15, although the target of this event is unknown. Apart from these alterations centered on TSC1 and TSC2, no other chromosomal gains or losses were identified in our cohort.

thumbnail
Fig 3. Copy neutral-LOH on chromosome 16p.

At least 8 different regions of copy neutral LOH were seen in 18 tumor samples with 16p LOH. The blue bars reflect the region of copy neutral LOH for each sample extending from the first to the last SNP with skewed allele frequency (AF <0.4 or >0.6). The gray bars represent the interval between the last SNP with normal AF (0.4 < AF < 0.6) and the first SNP with skewed AF on each side of the region with LOH, and reflect regions with no informative SNP markers to assess LOH.

https://doi.org/10.1371/journal.pgen.1006242.g003

Whole exome sequencing findings

Whole exome sequence analysis of 23 tumor normal pairs uncovered a median of 4 coding mutations per tumor (not including those identified in TSC1/TSC2, range 0–12), with 3 of 78 (4%) frame-shift deletions, 3 of 78 (4%) in-frame deletions, 2 of 78 (3%) nonsense, 2 of 78 (3%) splice site, and 68 of 78 (87%) missense mutations (Table 2). In addition, among the 68 missense mutations, 10 (15%) were neutral, 21 (31%) were of low potential, 31 (46%) were of medium potential, and 6 (9%) were of high potential for having functional effects according to Mutation Assessor [23, 24]. This corresponds to a mutation rate of 0.12 (range 0–0.36) non-synonymous mutations/Mb (excluding TSC1/TSC2). Three angiomyolipomas from P13 (S20, S23, and S24) had identical mutations in MAATS1 and NCF1 (Table 2), suggesting that they were derived from a single clone, as discussed in more detail below. Two samples from P13 had different mutations in TRIP12 (S24 c.5470-2A>G, S27 p.R1595Q) (Table 2), suggesting that alteration of this gene may contribute to development of these tumors. TRIP12 is an ubiquitin ligase for ARF, which is reported to suppress cell growth by activating p53 [25], and frequently harbors inactivating somatic point mutations (100 of 380 (26%) somatic mutations are inactivating [26, 27]) in diverse cancer types. However, although the p.R1595Q somatic mutation has been seen in a bladder cancer and a stomach cancer, the functional significance of that missense mutation is unknown.

thumbnail
Table 2. Somatic genetic alterations other than TSC1/TSC2 identified in 23 angiomyolipoma/LAM samples subject to tumor-normal paired exome sequencing.

https://doi.org/10.1371/journal.pgen.1006242.t002

Apart from MAATS1, NCF1, and TRIP12, no other genes were mutated in more than one tumor. Gene set enrichment analysis (GSEA) showed no enrichment of genes with mutations within hallmark gene sets, computational gene sets (including cancer gene neighborhoods, and cancer modules), GO gene sets (including GO cellular component, molecular function and GO molecular function gene dataset), or oncogenic signatures genesets [28] [29]. Furthermore, KEGG and Protein Interaction Network Module analyses found no enrichment of genes with mutations within any biological pathway or functional gene network [30] [31, 32]. Next, we examined each mutated gene for its potential genetic or physical interaction with TSC2 or TSC1, using BioGRID, Coremine, and Esyn databases. Only GSK3B was identified as an interactor with TSC2, based on their common involvement in PI3K-AKT-mTOR signaling [33] [34].

Singleton mutations were found in 3 of the 263 genes recurrently mutated across cancer as annotated by the ‘PanCan’ dataset [35, 36]: BAP1 (p.T517M in S23 from P13), ARHGAP35 (p.E1273A in S17 from P13), and SPEN (p.A1813T in S7 from P6) (Table 1). The BAP1 mutation has been previously seen in two cancer cell lines, derived from a Ewing sarcoma, and a lung adenocarcinoma [26, 27]. The other two mutations have not been reported previously [26, 27]. All are missense changes of uncertain functional effect.

In a search for potential structural rearrangements affecting these tumors, we performed whole genome sequencing of two tumor-normal pairs (S1 and S3). However, there was no evidence for genomic rearrangement in these samples using dRanger [37], and no new exonic mutations were identified by this analysis beyond the previous whole exome analysis.

Phylogenetic analysis of the 16 angiomyolipomas from TSC patient P13

LOH mapping and somatic mutations were used to develop a phylogenetic tree for the 16 renal angiomyolipomas from P13 (Fig 4). There were 5 distinct regions of LOH on 16p in these tumors (Fig 3, S3 Fig). Samples S20, S23, and S24 all had the same region of chromosome 16 LOH, and identical mutations in NCF1 and MAATS1, but also had unique mutations not seen the other two samples (Table 2). These observations suggest that these three tumors arose from a common precursor cell, sustaining these two early mutations prior to clonal expansion, followed by acquisition of additional mutations after dispersion (Fig 4). Five other tumors, S14, S15, S17, S18, S19, also had identical regions of LOH to the resolution limits of exome sequencing. However, these tumors did not share any sequence mutations, suggesting that mutations were acquired later in tumor development, or that they are not of common clonal ancestry. Five tumors, S21, S25, S26, S27, S29, had either unique regions of LOH (n = 3) or unique second hit TSC2 point mutations (n = 2), indicating that each was clonally distinct. Two angiomyolipomas (S16 and S28) from this patient had no evidence of TSC2 LOH, nor any confirmed somatic mutations, suggesting alternate mechanisms of TSC2 loss not detected by our analysis, or potential low tumor purity. Sample S22 also had low tumor purity (~10%) but showed evidence of TSC2 LOH without somatic mutations.

thumbnail
Fig 4. Phylogenetic tree of 16 angiomyolipomas from one TSC patient.

A model of angiomyolipoma development from a progenitor cell in the kidney is shown, including 3 tumors (S20, S23, S24) that had a definite common truncal precursor, and 5 tumors (S14, S15, S17, S18, S19) that may have had a common truncal precursor. Note that the extent of copy neutral LOH in S22 is uncertain due to low tumor purity. CN-LOH denotes copy neutral LOH.

https://doi.org/10.1371/journal.pgen.1006242.g004

Discussion

Consistent with extensive previous data and Knudson’s two-hit hypothesis [38], we identified biallelic inactivation of TSC2 or TSC1 in 30 of 32 angiomyolipomas and LAM tumors including point mutations, small indels, large genomic deletions, and copy neutral LOH. Interestingly, biallelic small mutations were more common in subjects without TSC (7 of 13, 54%) than in subjects with TSC (3 of 19, 12%; p = 0.05 Fisher exact test). Copy neutral LOH events were much more common in those with TSC. Using SNP allele frequency from our sequencing data, we were able to localize the site of aberrant cross-over events on 16p, and determine that at least 8 different regions were sites of mitotic recombination that led to copy number neutral LOH of TSC2.

We found that angiomyolipomas have very few non-synonymous exonic mutations (4 on average), and almost no recurrently mutated genes beyond TSC1 and TSC2. Although a number of genes with singleton mutations in this data set are mutated in other cancers, we suspect that most are background non-functional mutations that happened to be present in the initiating cells of these tumors, and are not tumor driver events. In support of this notion, 31 of 82 (38%) of the mutations appeared to be subclonal with variant AF < 0.15, and 20 of 82 (24%) had variant AF < 0.1 (Table 2). In addition, 87% were missense mutations, and of those 46% were neutral or had low potential for functional effects. Furthermore, we found no aggregation of the genes with mutations in any particular pathway using multiple statistical approaches and candidate gene sets. While it is possible that the mutations in ARHGAP35, BAP1, and SPEN are significant, all were missense mutations and the functional effects are uncertain. Two mutations in TRIP12 were also seen in different tumors, and while one is located in a splice site, the second (p.R1595Q) is of unknown functional effect.

Given its role in the mTOR signaling pathway [39], the singleton loss-of-function mutation in GSK3B (c.1195+1G>A, p.D399fs seen at clonal AF, 0.293, in S7) is arguably the strongest candidate to be a secondary driver mutation in this cohort. However, further study is required to examine the functional effect of this mutation in angiomyolipoma development.

Our data provide important information on clonality and tumor development for this class of tumors. Comparison of an abdominal LAM mass and a LAM cell cluster isolated from chylous fluid of a single patient (P1) revealed identical mutations in TSC2 and somatic alterations in other genes indicating a common clonal origin, and no candidate genetic event to account for metastasis and survival of the LAM cell cluster. Analysis of 16 angiomyolipomas (S14 –S29) derived from one patient (P13) identified at least 7 tumors with independent clonal origin with unique regions of 16p copy neutral LOH or unique TSC2 somatic mutations. These data strongly support the concept that the multiple angiomyolipomas that develop in most TSC adults are due to independent second hit mutations affecting TSC2 or TSC1 occurring in distinct progenitor cells, similar to what we have previously reported in TSC-associated RCC [40].

In P13, we also found evidence that single clones can seed multiple tumour masses, as 3 angiomyolipomas shared a common region of 16p LOH as well as two identical somatic mutations. These tumours each contained a small set of unique somatic mutations, indicative of further, independent clonal evolution. In addition, five tumor samples had the same region of 16p copy neutral LOH as well as other unique somatic mutations, and may or may not have been derived from a common precursor cell that became dispersed during tumor development. We also note that 3 of 16 angiomyolipomas from P13 had no somatic mutation findings, including two without second hit mutations in TSC2. Although these angiomyolipomas may be completely silent at the genomic level, it is also possible that the samples analyzed had a low tumor representation masking identification of both LOH and somatic mutations, or that epigenetic or other types of cancer genome variation underlie these tumors.

Our findings of the low somatic mutation rate of angiomyolipoma contrast sharply with cumulative cancer genome sequencing studies that have shown that adult malignancies have a wide range of non-synonymous mutations per cancer, ranging from 8 (chronic lymphocytic leukemia) to over 1,000 (colorectal cancer with microsatellite instability; other cancers with POLE mutations) [35, 41]. Angiomyolipoma, with a median of 4 non-synonymous mutations per tumor (range of 0 to 12) is well below the range of these common cancer types. However, infant and pediatric malignancies have a much lower mutation rate in general with a median of 4 to 13 among 5 pediatric cancers (acute lymphoblastic leukemia (ALL), medulloblastoma, neuroblastoma, glioblastoma, and rhabdoid cancers) [41]. In addition, several infant or pediatric malignancies appear to be due to a single genetic event: MLL-AF4 fusion in infant acute lymphoblastic leukemia [42], homozygous inactivation of SMARCB1 in pediatric rhabdoid tumors [43], and C11orf95-RELA fusion in brain ependymoma [44]. As noted by others [45], each of these three mutations affect chromatin remodeling genes that may be predicted to have wide-ranging effects on gene transcription, driving tumor development and growth in a pleiotypic manner. In addition, these singleton genetic events are seen in unique and rare malignancies likely derived from specific cell types at specific developmental stages of early childhood, suggesting that each (different) cell of origin is sensitive to the global chromatin and transcriptional regulatory derangements induced by the singleton genetic event.

Our findings suggest that angiomyolipomas fit this general model of pediatric tumor development, with few somatic mutations and a single critical target of inactivating mutation (in the case of angiomyolipoma, this is usually TSC2 and less commonly TSC1) that initiates and drives tumor development. Consistent with this model is the observation that small angiomyolipomas can be seen in the kidneys of young children with TSC and are present in up to 80% of TSC children by the age of 10 [46]. In addition, these tumors are typically indolent, growing slowly over a number of years, and usually do not require surgical intervention or other treatment, with a slowing in tumor growth in older adults. Furthermore these observations suggest that there is likely a specific cell type (currently unknown) that is resident in the early kidney that is sensitive to the tumor and growth promoting effects of TSC2 or TSC1 loss. Loss of TSC2 or TSC1 results in unregulated activation of mTORC1 with wide-ranging signaling, metabolomic, and transcriptional effects [39, 47, 48], which likely enable this singleton transforming event. In addition, there may be other clinically important consequences of TSC protein complex loss beyond mTORC1 activation that are not well understood.

Very limited similar studies have been performed for tumors occurring in other tumor suppressor gene syndromes with low growth potential. Whole exome sequencing of 7 neurofibromas from a patient with Neurofibromatosis type I led to identification of second hit mutations in NF1 in 5 of 7 tumors, and median of 0 (maximum 1) non-synonymous mutations in the rest of the exome [49]. Whole exome sequencing of 4 renal cell carcinomas (RCC) occurring in a 32 yo man with Von Hippel Lindau syndrome led to identification of unique second hit mutations in VHL in each RCC, as well as 13–23 other somatic non-synonymous mutations throughout the genome [50]. This increased level of somatic mutation in VHL RCCs is consistent with the more malignant behavior of these tumors compared to angiomyolipoma, but is less than that seen in RCC in those without VHL syndrome [50].

The remarkably low mutation rate and genomic stability observed in angiomyolipoma is consistent with continuing long-term responses (> 5 years in some cases) of these tumors to treatment with rapamycin and related drugs that inhibit mTORC1[51]. However, it is notable that that complete responses are virtually unknown and even partial responses, according to RECIST criteria, are rare.

Materials and Methods

Ethics statement

Informed written consent was obtained from all patients whose samples were studied. The study was approved by the Partners Human Research Committee of Brigham and Women’s Hospital (2007P000699, 2010P001818).

Patient samples

Thirty-two samples, comprising 30 fresh resected tumor tissues, one chylous pleural fluid sample, and one short-term culture from a resected angiomyolipoma, were collected from 15 patients with angiomyolipoma (n = 11) or concurrent angiomyolipoma and LAM (n = 4) (Table 1). All samples were collected from patients undergoing surgery as part of routine clinical care, and had been seen either at one of our hospitals (Brigham and Women’s Hospital or Massachusetts General Hospital), or in one case was referred from another institution. Normal adjacent tissue, saliva, and/or blood were also collected as control samples for 24 of the samples.

Sample processing

Genomic DNA was extracted from fresh frozen tumor specimens using Puregene DNA Purification Kit (QIAGEN, formerly Gentra/Puregene) according to the manufacturer’s protocol. The chylous fluid sample was subjected to filtering on a 40 μm nylon mesh (Cell Strainer from BD Biosciences). LAM cell clusters were collected and cultured in Lonza Clonetics SmGM-2 Smooth Muscle Growth Medium (containing 5% fetal bovine serum) in a humidified incubator at 5% CO2 [52]. Clusters of LAM cells were visually identified and a single large cluster was aspirated, and subjected to DNA extraction as described [52], and Phi29-based whole genome amplification using the GE Healthcare illustra GenomiPhi V2 DNA Amplification Kit.

Deep coverage sequence analysis of TSC1 and TSC2

To identify mutations in TSC1 and TSC2 with high sensitivity, we performed targeted sequence analysis of long-range PCR products at high read depth (>500x) as described previously [53]. This assay was used on the first six samples (subject P1-P6). Sequencing data output was analyzed using the Genome Analysis Toolkit [54], as well as custom software as described [53] to enable the detection of all sequence variants at allele frequency ≥1%.

Whole exome and genome sequencing

Exome sequencing was performed by the Broad Institute Genomics Platform and analyzed using a standard analytic pipeline deployed in the Firehose environment. Exome capture targeted 33 Mb in 193,094 exons in 18,863 genes [55]. Briefly, reads were aligned using bwa, followed by indel realignment and quality score recalibration using the Genome Analysis Toolkit [54]. Somatic mutations were identified from tumor-normal pairs using MuTect [56] and Indelocator [57]. This sequence data was also examined for mutations in TSC1 and TSC2 using custom software as described [53] to enhance detection of sequence variants. For the first 13 samples (Table 1), the median read depth was 86x for the targeted region, with median of 89% of target bases having a read depth of ≥ 20x, and median of 70% of target nt having a read depth of ≥ 50x. For the latter 19 angiomyolipoma samples (Table 1), median read depth was 191x for the targeted region, with median of 95.5% of target nt having a read depth of ≥ 20x, and median of 90% of target nt having a read depth of ≥ 50x.

Whole genome sequencing (WGS) was also performed at the Broad Institute following similar methods, with median read depth of 33x. Sequence variants were detected using the tools noted above and structural rearrangements were detected using dRanger [37].

Copy number changes and LOH mapping

Copy number profiles were derived from fractional coverage values of each exon compared with a panel of normals using Cap-Seg [58] and Allelic-CapSeg in the Clonal Evolution Exome Suite in the Firehose environment. Somatic LOH was quantified using ABSOLUTE [22] applied to variant and copy number calls derived from exome sequence data.

Detailed LOH mapping was performed for the paired tumor samples by: 1) identification of heterozygous SNPs in the normal sample using the following filters: read depth > 19, and variant AF between 0.40 and 0.60; and 2) determination of the variant AF for these SNPs in the tumor sample. In unpaired tumor samples, likely SNPs were identified by filtering for variants with variant AF between 0.05 and 0.95. Allele frequencies for these SNPs for each sample were graphed according to nucleotide position using Excel.

Review of variant calls

All sequence variants identified by Mutect and Indelocator were reviewed. Those identified in intronic and intergenic regions, and in ncRNA; or with a total number of variant reads < 3 or reads present only in one direction; or those seen at an allele frequency of < 5% were not considered further. Remaining variants were reviewed manually using Integrative Genomics Viewer (IGV) [59, 60], and were examined in multiple samples including the tumor and normal control. Artifacts, misaligned reads, and synonymous variants were discarded. The remaining variants were validated by one of three methods. First, both whole exome and whole genome sequencing was performed on two samples (S1 and S3, Table 1), and variants seen concordantly by the two analyses were considered confirmed. Second, RNA-Seq was performed on 5 samples (Table 1), and variants seen in both whole exome sequencing and RNA-Seq analyses were also considered confirmed (full details on RNA-Seq will be reported elsewhere). Third, variants were confirmed by Sanger sequencing for those seen at > 10% allele frequency, and by SNaPshot single nucleotide sequencing for those seen at allele frequency > 5% and < 10%. Sanger and SNaPshot extension products were analyzed on an ABI 3100 sequencer (Applied Biosystems, Carlsbad, CA, USA). SNaPshot experiments were run in duplicate.

Pathway analysis

Integrated pathway analysis was performed for genes in which somatic mutations were detected in two ways: one analysis for mutations found in the 16 distinct tumors from P13, and another for all somatic mutations found in 15 subjects. We used Gene Set Enrichment analysis (GSEA) [28, 29] and WebGestalt [31, 32] tools to search for pathway enrichment among the genes with somatic mutations in these tumors. We also explored each mutant gene using BioGRID, Coremine, and Esyn, to examine potential interactions between each gene and TSC2/TSC1 [33, 34]. We also searched COSMIC [61, 62], cBioPortal for Cancer Genomics [26, 27], that currently contains data from 105 cancer genomics studies and TumorPortal [36] to assess whether each of the identified genes and confirmed mutant variants had been seen previously in any type of cancer.

Supporting Information

S1 Fig. Kidneys resected from TSC patient P13 due to massive angiomyolipoma involvement.

Left kidney had weight 1.465 kg and size 28.4 x 11.9 x 11.5 cm, and right kidney had weight 0.515 kg and size 15.6 x 10.8 x 5.4 cm.

https://doi.org/10.1371/journal.pgen.1006242.s001

(DOCX)

S2 Fig. Homozygous genomic loss of TSC2 in sample S7. A. MLPA graph.

Top, analysis of S7 with a probe set covering TSC2. Control probes from other genomic sites have values from 0.81–1.25, while probes from TSC2 have values from 0.29 to 0.45, indicative of homozygous loss with some normal cell contamination. A control sample is shown at bottom. Note that 05ex refers to exon 5 of TSC2, etc. Probes are sorted by size. B. Capseg analysis visualized using IGV. A 199 kb region surrounding TSC2 is shown in IGV. Note that samples S4, S5, S6 show no copy number loss, while sample S7 shows homozygous copy number loss—dark blue, 0.42 copies—in a 50kb region; and single copy loss—light blue, 1.47 copies—in flanking regions. Gray regions are those with indeterminate copy number due to low coverage; white regions have copy number of 2.

https://doi.org/10.1371/journal.pgen.1006242.s002

(DOCX)

S3 Fig. LOH Mapping in renal angiomyolipoma and LAM tumor samples.

Graphs of allele fractions are shown for SNPs on chromosome 16 (Samples S1 –S29, S32) and chromosome 9 (Samples S30, S31).

https://doi.org/10.1371/journal.pgen.1006242.s003

(PPTX)

S1 Table. Glossary of genetic terms used in this publication.

https://doi.org/10.1371/journal.pgen.1006242.s004

(DOCX)

Acknowledgments

The authors thank the patients involved in this study for their participation.

Author Contributions

Conceived and designed the experiments: EPH DJK. Performed the experiments: KG IAM RY YYT CO MET YC YL. Analyzed the data: KG IAM TJP RY CO MET YC. Contributed reagents/materials/analysis tools: TJP JK NA GAF CLW SS MM GG JSB DJK. Wrote the paper: KG IAM TJP CO JK MET YC YL SS MM EPH DJK.

References

  1. 1. Henske EP, Neumann HP, Scheithauer BW, Herbst EW, Short MP, Kwiatkowski DJ. Loss of heterozygosity in the tuberous sclerosis (TSC2) region of chromosome band 16p13 occurs in sporadic as well as TSC-associated renal angiomyolipomas. Genes, Chromosomes & Cancer. 1995;13(4):295–8.
  2. 2. Karbowniczek M, Yu J, Henske EP. Renal angiomyolipomas from patients with sporadic lymphangiomyomatosis contain both neoplastic and non-neoplastic vascular structures. Am J Pathol. 2003;162(2):491–500. pmid:12547707.
  3. 3. Kozlowska J, Okon K. Renal tumors in postmortem material. Pol J Pathol. 2008;59(1):21–5. Epub 2008/07/29. pmid:18655367.
  4. 4. http://www.omim.org
  5. 5. Henske EP, Rasooly R, Siroky B, Bissler J. Tuberous sclerosis complex, mTOR, and the kidney: report of an NIDDK-sponsored workshop. Am J Physiol Renal Physiol. 2014;306(3):F279–83. Epub 2013/11/15. ajprenal.00525.2013 [pii] 10.1152/ajprenal.00525.2013. pmid:24226526; PubMed Central PMCID: PMC3920052.
  6. 6. Kwiatkowski DJ, Thiele EA, Whittemore VH. Tuberous sclerosis complex. Weinheim, Germany: Wiley-VCH; 2010. 409 p.
  7. 7. de Vries PJ, Whittemore VH, Leclezio L, Byars AW, Dunn D, Ess KC, et al. Tuberous sclerosis associated neuropsychiatric disorders (TAND) and the TAND Checklist. Pediatr Neurol. 2015;52(1):25–35. pmid:25532776; PubMed Central PMCID: PMC4427347.
  8. 8. Henske EP, McCormack FX. Lymphangioleiomyomatosis—a wolf in sheep's clothing. J Clin Invest. 2012;122(11):3807–16. Epub 2012/11/02. 58709 [pii] 10.1172/JCI58709. pmid:23114603.
  9. 9. Cudzilo CJ, Szczesniak RD, Brody AS, Rattan M, Krueger DA, Bissler JJ, et al. Lymphangioleiomyomatosis Screening in Women with Tuberous Sclerosis. Chest. 2013. Epub 2013/03/30. 1672803 [pii] pmid:23539171.
  10. 10. Thway K, Fisher C. PEComa: morphology and genetics of a complex tumor family. Ann Diagn Pathol. 2015;19(5):359–68. Epub 2015/07/07. S1092-9134(15)00090-8 [pii] 10.1016/j.anndiagpath.2015.06.003. pmid:26144278.
  11. 11. Bonetti F, Pea M, Martignoni G, Doglioni C, Zamboni G, Capelli P, et al. Clear cell ("sugar") tumor of the lung is a lesion strictly related to angiomyolipoma—the concept of a family of lesions characterized by the presence of the perivascular epithelioid cells (PEC). Pathology. 1994;26(3):230–6. Epub 1994/07/01. pmid:7991275.
  12. 12. El-Hashemite N, Zhang H, Henske EP, Kwiatkowski DJ. Mutation in TSC2 and activation of mammalian target of rapamycin signalling pathway in renal angiomyolipoma. Lancet. 2003;361(9366):1348–9. pmid:12711473.
  13. 13. Kenerson H, Folpe AL, Takayama TK, Yeung RS. Activation of the mTOR pathway in sporadic angiomyolipomas and other perivascular epithelioid cell neoplasms. Hum Pathol. 2007;38(9):1361–71. pmid:17521703.
  14. 14. Badri KR, Gao L, Hyjek E, Schuger N, Schuger L, Qin W, et al. Exonic mutations of TSC2/TSC1 are common but not seen in all sporadic pulmonary lymphangioleiomyomatosis. Am J Respir Crit Care Med. 2013;187(6):663–5. Epub 2013/03/19. 187/6/663 [pii] 10.1164/ajrccm.187.6.663. pmid:23504366.
  15. 15. Qin W BV, Malinowska I, Lu X, MacConaill L, Wu C-L, Kwiatkowski DJ. Angiomyolipoma have common mutations in TSC2 but no other common genetic events. PLOS One. 2011;in press.
  16. 16. Carsillo T, Astrinidis A, Henske EP. Mutations in the tuberous sclerosis complex gene TSC2 are a cause of sporadic pulmonary lymphangioleiomyomatosis. Proc Natl Acad Sci U S A. 2000;97(11):6085–90. pmid:10823953
  17. 17. Henske EP, Scheithauer BW, Short MP, Wollmann R, Nahmias J, Hornigold N, et al. Allelic loss is frequent in tuberous sclerosis kidney lesions but rare in brain lesions. Am J Hum Genet. 1996;59:400–6. pmid:8755927
  18. 18. Cai X, Pacheco-Rodriguez G, Fan QY, Haughey M, Samsel L, El-Chemaly S, et al. Phenotypic characterization of disseminated cells with TSC2 loss of heterozygosity in patients with lymphangioleiomyomatosis. Am J Respir Crit Care Med. 2010;182(11):1410–8. pmid:20639436; PubMed Central PMCID: PMC3029931.
  19. 19. Pan CC, Chung MY, Ng KF, Liu CY, Wang JS, Chai CY, et al. Constant allelic alteration on chromosome 16p (TSC2 gene) in perivascular epithelioid cell tumour (PEComa): genetic evidence for the relationship of PEComa with angiomyolipoma. J Pathol. 2008;214(3):387–93. pmid:18085521.
  20. 20. Kwiatkowski DJ. Genetics of Tuberous Sclerosis Complex. In: Kwiatkowski DJ, Whittemore VH, Thiele EA, editors. Tuberous Sclerosis Complex. Weinheim, Germany: Wiley-VCH; 2010. p. 29–60.
  21. 21. Kozlowski P, Roberts P, Dabora S, Franz D, Bissler J, Northrup H, et al. Identification of 54 large deletions/duplications in TSC1 and TSC2 using MLPA, and genotype-phenotype correlations. Hum Genet. 2007;121(3–4):389–400. pmid:17287951.
  22. 22. Carter SL, Cibulskis K, Helman E, McKenna A, Shen H, Zack T, et al. Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol. 2012;30(5):413–21. pmid:22544022; PubMed Central PMCID: PMC4383288.
  23. 23. Reva B, Antipin Y, Sander C. Predicting the functional impact of protein mutations: application to cancer genomics. Nucleic Acids Res. 2011;39(17):e118. pmid:21727090; PubMed Central PMCID: PMC3177186.
  24. 24. http://mutationassessor.org/r2/
  25. 25. Park Y, Yoon SK, Yoon JB. The HECT domain of TRIP12 ubiquitinates substrates of the ubiquitin fusion degradation pathway. J Biol Chem. 2009;284(3):1540–9. Epub 2008/11/26. M807554200 [pii] pmid:19028681.
  26. 26. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer discovery. 2012;2(5):401–4. pmid:22588877; PubMed Central PMCID: PMC3956037.
  27. 27. http://www.cbioportal.org
  28. 28. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. pmid:16199517; PubMed Central PMCID: PMC1239896.
  29. 29. http://software.broadinstitute.org/gsea/index.jsp
  30. 30. Kanehisa M. The KEGG database. Novartis Found Symp. 2002;247:91–101; discussion -3, 19–28, 244–52. Epub 2003/01/24. pmid:12539951.
  31. 31. Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 2013;41(Web Server issue):W77–83. pmid:23703215; PubMed Central PMCID: PMC3692109.
  32. 32. http://bioinfo.vanderbilt.edu/webgestalt/
  33. 33. Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M. BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006;34(Database issue):D535–9. Epub 2005/12/31. 34/suppl_1/D535 [pii] 10.1093/nar/gkj109. pmid:16381927; PubMed Central PMCID: PMC1347471.
  34. 34. http://thebiogrid.org/
  35. 35. Lawrence MS, Stojanov P, Mermel CH, Robinson JT, Garraway LA, Golub TR, et al. Discovery and saturation analysis of cancer genes across 21 tumour types. Nature. 2014;505(7484):495–501. pmid:24390350; PubMed Central PMCID: PMC4048962.
  36. 36. http://www.tumorportal.org
  37. 37. http://www.broadinstitute.org/cancer/cga/dranger
  38. 38. Knudson AG Jr. Mutation and cancer: statistical study of retinoblastoma. Proc Natl Acad Sci U S A. 1971;68:820–3. pmid:5279523
  39. 39. Tee AR, Sampson JR, Pal DK, Bateman JM. The role of mTOR signalling in neurogenesis, insights from tuberous sclerosis complex. Semin Cell Dev Biol. 2016;52:12–20. pmid:26849906.
  40. 40. Tyburczy ME, Jozwiak S, Malinowska IA, Chekaluk Y, Pugh TJ, Wu CL, et al. A shower of second hit events as the cause of multifocal renal cell carcinoma in tuberous sclerosis complex. Hum Mol Genet. 2015;24(7):1836–42. pmid:25432535; PubMed Central PMCID: PMC4355019.
  41. 41. Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA Jr., Kinzler KW. Cancer genome landscapes. Science. 2013;339(6127):1546–58. pmid:23539594; PubMed Central PMCID: PMC3749880.
  42. 42. Andersson AK, Ma J, Wang J, Chen X, Gedman AL, Dang J, et al. The landscape of somatic mutations in infant MLL-rearranged acute lymphoblastic leukemias. Nat Genet. 2015;47(4):330–7. pmid:25730765; PubMed Central PMCID: PMC4553269.
  43. 43. Lee RS, Stewart C, Carter SL, Ambrogio L, Cibulskis K, Sougnez C, et al. A remarkably simple genome underlies highly malignant pediatric rhabdoid cancers. J Clin Invest. 2012;122(8):2983–8. pmid:22797305; PubMed Central PMCID: PMC3408754.
  44. 44. Parker M, Mohankumar KM, Punchihewa C, Weinlich R, Dalton JD, Li Y, et al. C11orf95-RELA fusions drive oncogenic NF-kappaB signalling in ependymoma. Nature. 2014;506(7489):451–5. pmid:24553141; PubMed Central PMCID: PMC4050669.
  45. 45. Greaves M. When one mutation is all it takes. Cancer Cell. 2015;27(4):433–4. pmid:25873168.
  46. 46. Ewalt DH, Sheffield E, Sparagana SP, Delgado MR, Roach ES. Renal lesion growth in children with tuberous sclerosis complex. J Urol. 1998;160(1):141–5. pmid:9628635
  47. 47. Dibble CC, Manning BD. Signal integration by mTORC1 coordinates nutrient input with biosynthetic output. Nat Cell Biol. 2013;15(6):555–64. pmid:23728461; PubMed Central PMCID: PMC3743096.
  48. 48. Lamming DW, Sabatini DM. A Central role for mTOR in lipid homeostasis. Cell Metab. 2013;18(4):465–9. pmid:23973332; PubMed Central PMCID: PMC3818790.
  49. 49. Emmerich D, Zemojtel T, Hecht J, Krawitz P, Spielmann M, Kuhnisch J, et al. Somatic neurofibromatosis type 1 (NF1) inactivation events in cutaneous neurofibromas of a single NF1 patient. Eur J Hum Genet. 2015;23(6):870–3. pmid:25293717; PubMed Central PMCID: PMC4795057.
  50. 50. Fisher R, Horswell S, Rowan A, Salm MP, de Bruin EC, Gulati S, et al. Development of synchronous VHL syndrome tumors reveals contingencies and constraints to tumor evolution. Genome Biol. 2014;15(8):433. pmid:25159823; PubMed Central PMCID: PMC4166471.
  51. 51. Bissler JJ, Kingswood JC, Radzikowska E, Zonnenberg BA, Frost M, Belousova E, et al. Everolimus for renal angiomyolipoma in patients with tuberous sclerosis complex or sporadic lymphangioleiomyomatosis: extension of a randomized controlled trial. Nephrol Dial Transplant. 2016;31(1):111–9. pmid:26156073.
  52. 52. Malinowska-Kolodziej I, Finlay G, Campbell G, Garcia-Moliner ML, Weinstein B, Doughty CA, et al. Lymphangioleiomyomatosis: cause of a malignant chylous pleural effusion. J Clin Oncol. 2010;28(1):e4–6. Epub 2009/11/11. JCO.2009.23.3890 [pii] 10.1200/JCO.2009.23.3890. pmid:19901112.
  53. 53. Tyburczy ME, Dies KA, Glass J, Camposano S, Chekaluk Y, Thorner AR, et al. Mosaic and Intronic Mutations in TSC1/TSC2 Explain the Majority of TSC Patients with No Mutation Identified by Conventional Testing. PLoS Genet. 2015;11(11):e1005637. pmid:26540169; PubMed Central PMCID: PMC4634999.
  54. 54. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. Epub 2010/07/21. gr.107524.110 [pii] 10.1101/gr.107524.110. pmid:20644199; PubMed Central PMCID: PMC2928508.
  55. 55. Chapman MA, Lawrence MS, Keats JJ, Cibulskis K, Sougnez C, Schinzel AC, et al. Initial genome sequencing and analysis of multiple myeloma. Nature. 2011;471(7339):467–72. pmid:21430775; PubMed Central PMCID: PMC3560292.
  56. 56. Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31(3):213–9. Epub 2013/02/12. nbt.2514 [pii] 10.1038/nbt.2514. pmid:23396013.
  57. 57. http://www.broadinstitute.org/cancer/cga/indelocator
  58. 58. https://github.com/aaronmck/CapSeg
  59. 59. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6. Epub 2011/01/12. nbt.1754 [pii] 10.1038/nbt.1754. pmid:21221095; PubMed Central PMCID: PMC3346182.
  60. 60. https://www.broadinstitute.org/igv
  61. 61. Forbes SA, Tang G, Bindal N, Bamford S, Dawson E, Cole C, et al. COSMIC (the Catalogue of Somatic Mutations in Cancer): a resource to investigate acquired mutations in human cancer. Nucleic Acids Res. 2010;38(Database issue):D652–7. Epub 2009/11/13. gkp995 [pii] 10.1093/nar/gkp995. pmid:19906727; PubMed Central PMCID: PMC2808858.
  62. 62. http://cancer.sanger.ac.uk/cosmic/