Exome profiling of primary, metastatic and recurrent ovarian carcinomas in a BRCA1-positive patient
© Zhang et al.; licensee BioMed Central Ltd. 2013
Received: 17 October 2012
Accepted: 13 March 2013
Published: 22 March 2013
Ovarian carcinoma is a common, and often deadly, gynecological cancer. Mutations in BRCA1 and BRCA2 genes are present in at least a fifth of patients. Uncovering other genes that become mutated subsequent to BRCA1/BRCA2 inactivation during cancer development will be helpful for more effective treatments.
We performed exome sequencing on the blood, primary tumor, omental metastasis and recurrence following therapy with carboplatin and paclitaxel, from a patient carrying a BRCA1 S1841R mutation.
We observed loss of heterozygosity in the BRCA1 mutation in the primary and subsequent tumors, and somatic mutations in the TP53 and NF1 genes were identified, suggesting their role along with BRCA1 driving the tumor development. Notably, we show that exome sequencing is effective in detecting large chromosomal rearrangements such as deletions and amplifications in cancer. We found that a large deletion was present in the three tumors in the regions containing BRCA1, TP53, and NF1 mutations, and an amplification in the regions containing MYC. We did not observe the emergence of any new mutations among tumors from diagnosis to relapse after chemotherapy, suggesting that mutations already present in the primary tumor contributed to metastases and chemotherapy resistance.
Our findings suggest that exome sequencing of matched samples from one patient is a powerful method of detecting somatic mutations and prioritizing their potential role in the development of the disease.
KeywordsDriver mutations Gynecological cancer Hereditary cancer Next generation sequencing Tumor suppressor genes Chromosomal rearrangements
Ovarian carcinoma (OC) is the leading cause of death from gynecological cancer in western countries. The most important predisposing factors are germline mutations in inherited cancer susceptibility genes, most notably BRCA1, BRCA2, RAD51C, RAD51D and the mismatch repair genes [1, 2]. Recently, next generation (exome) sequencing of 316 OC revealed that over 20 percent of these cancers carried either somatic or germline inactivating mutations in either BRCA1 or BRCA2, thus emphasizing the importance of these two genes in the pathogenesis of OC . Notably, about a quarter of women diagnosed with OC in their fifth decade will carry a BRCA1 or BRCA2 mutation . Several studies have observed that BRCA1 and BRCA2 mutation carriers tend to have a better outcome than stage-matched non-carriers, and that this better outcome is largely attributable to the combination of BRCA mutation status and DNA-damaging chemotherapeutic drugs such as cisplatinum . There have also been case reports of rare cures achieved in BRCA1/2 carriers with ovarian and other cancers following other, older treatments such as melphalan . Together, these findings suggest that optimal alignment of chemotherapeutic agents with both host and tumor genetic events is possible and is in fact required to achieve improved outcomes. To further understand the interaction between treatment, host genetics and tumor-specific mutations, we extracted DNA from four sources obtained from a single patient carrying a deleterious mutation in BRCA1 (blood, primary tumor, omental metastasis and relapse (recurrence) following standard post-operative therapy with carboplatin and paclitaxel). These four DNA samples were then subjected to whole exome sequencing, thus allowing us to identify tumor-specific variants and to determine potential changes in allele frequencies and emergence of new variants in the different tumor samples.
Tumor samples used for exome sequencing
Tumor samples were kept at −80 degrees Celsius. All examined tumor blocks contained poorly differentiated serous adenocarcinoma (Figure 2). The histiotype was ascertained in routine histological slides obtained from the same tumor which was fixed in formalin and sections were obtained from paraffin-embedded tissue. This was done because cell morphology was not preserved well enough to provide information on the histiotype of the malignant cells. The serous histiotype was further demonstrated by immunohistochemistry: the neoplastic cells of all tumor samples stained strongly and diffusely for CA-125, p16, TP53, Ki-67 and WTI. They failed to stain for caldesmon, fascin and only very weakly and focally for B-cadherin. This immunohistochemical profile is consistent with serous differentiation.
Exome sequencing and SNP/small indel detection
Exome sequencing was applied on the primary tumor, the omental metastasis, the tumor present at relapse, and the blood from the patient to identify somatic mutations. Exomes were captured from a total of 3 μg of genomic DNA, using the Illumina TruSeq exome enrichment kit, according to manufacturer’s protocols. Samples were sequenced using one lane of paired-end, 100 bp reads on Illumina Hiseq for each sample. We ensured that only read pairs with both mates present were subsequently used. Adaptor sequences and quality trimmed reads were removed using Fastx toolkit (http://hannonlab.cshl.edu/fastx_toolkit/). Reads that passed quality control were aligned to the UCSC hg19 reference genome with BWA . Duplicate reads were marked using Picard (http://picard.sourceforge.net/) and were excluded from downstream analyses. SAMtools was used to call SNV and indel variants . Next, we applied additional quality control measures to all identified raw variants based on the following criteria: 1) The Phred-like score is no less than 20 for SNPs and 50 for indels; 2) the read coverage of no less than three reads per base; 3) at least three and 10% of covering reads had to support the alternate base for the primary tumor sample. Finally, we used Annovar to identify SNVs and indels that located in protein coding regions as well as variants affecting canonical splice sites .
Copy number variant detection
Copy number variant (CNV) detection was done by comparing normalized read coverage or read-depth between the blood and each of the primary, metastatic, and recurrent tumors, using an algorithm based on ExomeCNV . Read-depth was normalized to Reads Per Kilobase of exon model per Million mapped reads (RPKM)  for each exon, and the log ratio of RPKM values were calculated. Log ratios serve as input for DNAcopy, which segments chromosomal regions based on similar log ratios . In this study, because the use of exome sequencing data is still not well proven in CNV detection, we refrained from attempting to identify small structural variants and concentrated on larger segments, which we can detect with high confidence. In order to identify large scale rearrangements, the DNAcopy outputs were smoothed by removing small CNV calls and merging adjacent segments. Some large CNVs may be represented by more than one segment because they span regions where exonic data are unavailable.
In our present analysis, the identified regions contain at least 100 exons which collectively show deviation from the expected. The probability that all of these show the same deviation by random chance is negligible (i.e. if p-values for each exon within the segments are combined using Fisher’s Method, the resulting p-value approaches zero).
Results and discussion
Numbers of variant calls from exome sequencing results
Variants after quality check
Rare variants after filtering
Validated somatic variants
We proceeded to examine the change in frequency of the BRCA1 missense mutation (chr17, 41197766, S1841R) and observed increasing allele frequencies of this mutation: 0.48 in the blood, 0.57 in the primary tumor, 0.76 in the metastasis, and 0.72 in the recurrence. Upon validation using Sanger sequencing, this mutation showed consistent increase in frequency: 0.39 in the blood, 0.50 in the primary tumor, 0.68 in the metastasis, and 0.78 in recurrence. We note that the measurements from exome appear more accurate than from Sanger sequencing, because the allele frequency from exome sequencing of the inherited BRCA1 mutation in the blood sample was closer to the expected 0.5, representing heterozygosity. Although we observed increase in frequency of this mutation from blood to tumor samples, we did not observe complete loss of the wild-type allele in the tumors. Based on previous investigations of series of BRCA1 mutation-positive patients  the primary, metastatic and recurrent tumors will frequently exhibit complete loss of heterozygosity (LOH), and therefore the mutant allele frequency in the tumors should be close to 1, instead of 0.57 - 0.76, suggesting that the tumor samples may contain considerable proportion of non-malignant tissue. Allowing for sampling issues, it does appear that the frozen primary tissue (equivalent paraffin section images shown in Figure 2A) contains a considerable amount of non-malignant tissue, whereas, as shown in Figure 2B, the percentage of malignant tissue in the omental biopsy is higher (fat cells take up some of the sample, top of the figure). This is even more evident in Figure 2C, where there appears to be very little non-malignant tissue present. Further corroborating these data, CNV detection results showed that the allelic frequency of all the identified large deletions/duplications is increased from primary tumor to metastatic and recurrent tumors. Concurrently, we find no evidence for de novo alleles in the primary tumor that are absent in the subsequent tumors – which would have indicated that the primary tumor contained a mixture of different malignant clones. Thus, we hypothesize that the primary tumor sample we obtained for sequencing contained a relatively larger proportion of normal tissue than the metastases. The increased mutant allele frequencies among tumor samples are likely to reflect a more pure tumor sample, rather than a selection process. Moreover, CNV detection suggested that the region (17q11-17q21) containing BRCA1 gene was deleted in all tumors, including the primary. This result is consistent with LOH, and that in this patient, the inherited mutation and the somatic deletion in BRCA1 together initiated the tumor growth.
The above observation implies that most of the detected mutations were present in the primary tumor and that very little, if any selection has occurred thereafter. The most compelling hypothesis regarding the origin of BRCA1-related high-grade serous ovarian carcinoma is that in fact the majority of them arise in the fallopian tube . Our findings suggest that most of the critical tumor-driving clonal evolution occurs very early in the life of BRCA1-related highgrade serous carcinomas. One can reasonably speculate that the three tumors we studied here were all in fact “secondary” to the primary origin of the tumor and metastases from the now obscured primary tumor, likely in the fallopian tube. Surgery and chemotherapy failed to eradicate the original clone. Furthermore, when taking into account the relatively lower purity of the primary tumor, it is highly likely that most of the somatic mutations detected in this study were already present at high allelic frequency and high level of clonality in the tumor arising in the ovary. In agreement with our data, Castellarin et al. have recently suggested that in high-grad serous carcinoma patients, most somatic mutations found in recurrent tumors during platinum-based chemotherapy were present in primary tumors . Our data thus suggests that little genetic evolution of the tumor has taken place from time of diagnosis to relapse following three courses of highly-active chemotherapy. It is possible that the 2.5 fold increase in allele frequency of the NF1 mutation from the primary tumor to the metastasis (Table 1) indicates that this mutation appeared in the primary tumor later than for example, TP53 mutation but was required for the full metastatic phenotype. It is likely that the primary tumor that is detected in patients is descended from cells that already contain a significant and potentially lethal mutational load.
Sanger sequencing confirmed somatic mutations with increased frequencies in tumor samples
Mutant allele frequency from exome sequencing
Mutant allele frequency from sanger sequencing
splice site SNV
splice site SNV
Copy number variants (CNVs) that were detected in primary, metastatic and recurrent tumors
CNV segements indicating deletion/amplification
Mean log ratio
Mean log ratio
Mean log ratio
17p + 17q11-17q21
This work used whole exome capture and massively parallel DNA sequencing to study targeted candidate mutations in selected genes, as well as performing a “hypothesis-free” analysis where we aimed to identify potential driver mutations by identifying variants with increased proportion of mutant alleles. Genetic evolution of tumors from diagnosis to relapse after chemotherapy was not observed. Instead, we suggest that most of the critical tumor-driving and chemotherapy resistant mutations were already present in the primary tumor. We show that high-throughput sequencing is effective in detecting large chromosomal rearrangements such as deletions and amplifications that occur in cancer. It is notable that the patient responded very poorly to platinum-based therapy; relapse after only 3 course of therapy usually betokens a very poor survival. This early platinum failure is somewhat less common in BRCA1-related cancer than in non-hereditary ovarian cancer , and it seems unlikely that this failure is related to type of mutation (i.e. missense mutation) that was present in this patient. The large number of deleterious somatic mutations present in the primary tumor likely contributed to the rapid progression of the disease. It will be important to conduct studies such as ours in large numbers of patients to establish whether specific exomic profiles at initial diagnosis are associated with subsequent resistance to standard chemotherapy. In these situations, alternative forms of first-line therapy may be chosen. As many similar studies are going to be carried out in the near future, correlation of such candidate lists across patients will provide unprecedented information regarding recurrent mutations in specific genes responsible for metastasis and resistance to therapy. In addition, pathway analysis of the mutated genes will allow definition of the functional pathways involved in the above processes.
We thank Ms. Sonya Zaor for her clinical care and Rachel Silva-Smith for help with preparing the Figures. This work was supported by the Marsha Rivkin Centre for Ovarian Cancer Research (WDF, JM); the Weekend to End Women’s Cancer (WDF, WHG); and the Réseau de Médecine Génétique Appliquée (JZ); and the Canadian Institute for Health Research (EL). JM is a recipient of the Canada Research Chair. We would also like to acknowledge the Genome Quebec High Throughput Sequencing Platform for performing the exome sequencing.
- Loveday C, Turnbull C, Ramsay E, Hughes D, Ruark E, Frankum JR, Bowden G, Kalmyrzaev B, Warren-Perry M, Snape K, Adlard JW, Barwell J, Berg J, Brady AF, Brewer C, Brice G, Chapman C, Cook J, Davidson R, Donaldson A, Douglas F, Greenhalgh L, Henderson A, Izatt L, Kumar A, Lalloo F, Miedzybrodzka Z, Morrison PJ, Paterson J, Porteous M, Rogers MT, Shanley S, Walker L, Eccles D, Evans DG, Renwick A, Seal S, Lord CJ, Ashworth A, Reis-Filho J, Antoniou AC, Rahman N: Germline mutations in RAD51D confer susceptibility to ovarian cancer. Nat Genet. 2011, 43 (9): 879-882. 10.1038/ng.893.View ArticlePubMedPubMed CentralGoogle Scholar
- Walsh T, Casadei S, Lee MK, Pennil CC, Nord AS, Thornton AM, Roeb W, Agnew KJ, Stray SM, Wickramanayake A, Norquist B, Pennington KP, Garcia RL, King M, Swisher EM: Mutations in 12 genes for inherited ovarian, fallopian tube, and peritoneal carcinoma identified by massively parallel sequencing. PNAS. 2011, 108 (44): 18032-18037. .View ArticlePubMedPubMed CentralGoogle Scholar
- The Cancer Genome Atlas Research Network: Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474 (7353): 609-615. 10.1038/nature10166.View ArticlePubMed CentralGoogle Scholar
- Zhang S, Royer R, Li S, McLaughlin JR, Rosen B, Risch HA, Fan I, Bradley L, Shaw PA, Narod SA: Frequencies of BRCA1 and BRCA2 mutations among 1,342 unselected patients with invasive ovarian cancer. Gynecol Oncol. 2011, 121 (2): 353-357. 10.1016/j.ygyno.2011.01.020.View ArticlePubMedGoogle Scholar
- Bolton KL, Chenevix- Trench G, Goh C, Sadetzki S, Ramus SJ, Karlan BY, Lambrechts D, Despierre E, Barrowdale D, McGuffog L, Healey S, Easton DF, Sinilnikova O, BenÃtez J, GarcÃa MJ, Neuhausen S, Gail MH, Hartge P, Peock S, Frost D, Evans DG, Eeles R, Godwin AK, Daly MB, Kwong A, Ma ESK, LÃ¡zaro C, Blanco I, Montagna M: Association between BRCA1 and BRCA2 mutations and survival in women with invasive epithelial ovarian cancer. JAMA. 2012, 307 (4): 382-389. 10.1001/jama.2012.20.View ArticlePubMedPubMed CentralGoogle Scholar
- Osher DJ, Kushner YB, Arseneau J, Foulkes WD: Melphalan as a treatment for BRCA-related ovarian carcinoma: can you teach an old drug new tricks?. J Clin Pathol. 2011, 64 (10): 924-926. .View ArticlePubMedGoogle Scholar
- Lee MS, Green R, Marsillac SM, Coquelle N, Williams RS, Yeung T, Foo D, Hau DD, Hui B, Monteiro ANA, Glover JNM: Comprehensive analysis of missense variations in the BRCT domain of BRCA1 by structural and functional assays. Cancer Res. 2010, 70 (12): 4880-4890. 10.1158/0008-5472.CAN-09-4563.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: 1000 Genome project data, processing subgroup: the sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38 (16): e164-e164. 10.1093/nar/gkq603.View ArticlePubMedPubMed CentralGoogle Scholar
- Ge B, Gurd S, Gaudin T, Dore C, Lepage P, Harmsen E, Hudson TJ, Pastinen T: Survey of allelic expression using EST mining. Genome Res. 2005, 15 (11): 1584-1591. 10.1101/gr.4023805.View ArticlePubMedPubMed CentralGoogle Scholar
- Sathirapongsasuti JF, Lee H, Horst BAJ, Brunner G, Cochran AJ, Binder S, Quackenbush J, Nelson SF: Exome sequencing-based copy-number variation and loss of heterozygosity detection: ExomeCNV. Bioinformatics. 2011, 27 (19): 2648-2654. .View ArticlePubMedPubMed CentralGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.View ArticleGoogle Scholar
- Venkatraman ES, Olshen AB: A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics. 2007, 23 (6): 657-663. 10.1093/bioinformatics/btl646.View ArticlePubMedGoogle Scholar
- Piek JMJ, van Diest PJ, Zweemer RP, Jansen JW, Poort-Keesom R, Menko FH, Gille JJP, Jongsma APM, Pals G, Kenemans P, Verheijen RHM: Dysplastic changes in prophylactically removed Fallopian tubes of women predisposed to developing ovarian cancer. J Pathol. 2001, 195 (4): 451-456. 10.1002/path.1000.View ArticlePubMedGoogle Scholar
- Castellarin M, Milne K, Zeng T, Tse K, Mayo M, Zhao Y, Webb JR, Watson PH, Nelson BH, Holt RA: Clonal evolution of high-grade serous ovarian carcinoma from primary to recurrent disease. J Pathol. 2012, 229 (4): 515-524.View ArticlePubMedGoogle Scholar
- Sangha N, Wu R, Kuick R, Powers S, Mu D, Fiander D, Yuen K, Katabuchi H, Tashiro H, Fearon ER, Cho KR: Neurofibromin 1 (NF1) defects are common in human ovarian serous carcinomas and co-occur with TP53 mutations. Neoplasia. 2008, 10 (12): 1362-1372.View ArticlePubMedPubMed CentralGoogle Scholar
- Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR: A method and server for predicting damaging missense mutations. Nat Meth. 2010, 7 (4): 248-249. 10.1038/nmeth0410-248.View ArticleGoogle Scholar
- Jones PA, Baylin SB: The epigenomics of cancer. Cell. 2007, 128 (4): 683-692. 10.1016/j.cell.2007.01.029.View ArticlePubMedPubMed CentralGoogle Scholar
- Downward J: Targeting RAS signalling pathways in cancer therapy. Nat Rev Cancer. 2003, 3 (1): 11-22.View ArticlePubMedGoogle Scholar
- Lord CJ, Ashworth A: The DNA damage response and cancer therapy. Nature. 2012, 481 (7381): 287-294. 10.1038/nature10760.View ArticlePubMedGoogle Scholar
- Brown LA, Kalloger SE, Miller MA, Shih I, McKinney SE, Santos JL, Swenerton K, Spellman PT, Gray J, Gilks CB, Huntsman DG: Amplification of 11q13 in ovarian carcinoma. Genes Chromosom Cancer. 2008, 47 (6): 481-489. 10.1002/gcc.20549.View ArticlePubMedGoogle Scholar
- Johannsdottir H, Johannesdottir G, Agnarsson B, Eerola H, Arason A, Johannsson O, Heikkilä P, Egilsson V, Olsson H, Borg A, Nevanlinna H, Barkardottir R: Deletions on chromosome 4 in sporadic and BRCA mutated tumors and association with pathological variables. Anticancer Res. 2004, 24 (5): 2681-2688.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/146/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.