S1: Patient Samples and Associated Genomic Information Tumor Samples Were Ascertained from Three Independent Cohorts
Total Page:16
File Type:pdf, Size:1020Kb
S1: Patient samples and associated genomic information Tumor samples were ascertained from three independent cohorts. 1) The TCGA cohort of HGSC was collated by investigators from international hospital and universities, mostly from North America (1). Germline, somatic mutation and methylation status of HRR pathway members was available to us for 316 tumors and patients. From these we selected cases that had either somatic or germline pathogenic BRCA1/2 mutation or no apparent disruption of the BRCA pathway (wild-type) for gene expression and copy number analysis, including: • 280 cases with gene expression profiles, including 210 cases for which the molecular subtypes were identified in a previous publication (2) (27 BRCA1 mutated, 28 BRCA2 mutated, 10 BRCA1 methylated and 145 wild type). • 204 cases with copy number data (34 BRCA1 mutated, 30 BRCA2 mutated, 140 wild-type). • Gene expression profiles of additional 196 samples for clinical validation of the BRCA1/2 classifier. Mutation status of these samples was not available at the time of this study. 2) The Australian Ovarian Cancer Study (AOCS) is a population based case-control cohort of ovarian cancer patients ascertained at diagnosis between 2002 and 2006 (3). Research associated with the use of AOCS samples and clinical data was approved by the Human Research Ethics Committees at the Peter MacCallum Cancer Centre. A cohort of high grade serous samples was selected from tumors where germline BRCA1 and BRCA2 mutations (4) and Affymetrix U133 2.0 data were available (3). As described below, all 132 tumors were screened for somatic BRCA1/2 mutations using a high-resolution melt analysis (5), and for methylation of the BRCA1, PALB2 and FANCF gene promoters using methylation-sensitive high-resolution melting technology (6, 7) after bisulfite conversion (EpiTect Bisulfite Kit, Qiagen). One hundred and eleven of the 132 tumors profiled could be classified into one of the four molecular subtypes of HGSC previously described (3) and were used in the subsequent gene expression analysis. 3) Gene expression data for a cohort of 61 ovarian cancer tumors with mixed histologies, ascertained from the Memorial Sloan-Kettering Cancer Centre, and generated using in-house cDNA microarray chips manufactured at the National Cancer Institute Microarry Facility, was publically available as supplementary data provided by the journal 1 (8). 1 http://jnci.oxfordjournals.org/content/94/13/990/suppl/DC1 Study ID Gene Exon Nomenclature Location 2476 BRCA1 11 c.3302G>A Germline 3958 BRCA1 2 c.80+3A>C Germline 2410 BRCA2 11 c.[4094G>A] Germline 4256 BRCA2 3 c.[68-7T>A] Germline 1096 BRCA2 22 c.[8830A>T] Germline 1033 BRCA2 15 c.7504C>T Germline 6496 BRCA2 18 c.[7985C>A] Germline 66211 BRCA1 15 c.4669G>C Germline 20392 BRCA1 7 c.[305C>G] Germline 32023 BRCA1/BRCA2 17/11 c.[4987-20A>G/c.[5446A>C] Germline Supplementary Table 1 : Germline unclassified sequence variants identified in 132 AOCS tumors. Pathogenicity of sequence variants was assessed by reference to the Breast Information Core database (BIC), protein function prediction algorithims (SIFT, PolyPhen 2.0, AGVGD), splice site prediction algorithms (NNSPLICE, MaxEntScan, Splice Site Finder-Like, Human Splicing Finder) and literature searches as necessary. This process was standardised by use of the Alamut program (Interactive Biosoftware) as described in (4). 1 also found to be methylated at the BRCA1 gene promoter. 2 co-existing BRCA1 pathogenic mutation 3 co-existing BRCA2 pathogenic mutation, and BRCA1 gene promoter methylation S2: Somatic BRCA1/2 mutation detection using High-Resolution Melting (HRM) Analysis DNA extracted from fresh-frozen primary tumor tissue was screened for mutations in all coding exons, and intron-exon boundaries, of BRCA1 and BRCA2 using high resolution melting (HRM) analysis, as previously described (5) with the LightCycler®480 thermocycler (Roche). The success of the melting curve analysis relies on a heteroduplex being formed by the wild-type and mutant DNA sequences during the PCR step. To account for potential loss of heterozygosity at either BRCA1 or BRCA2 locus, which would result in a duplex being formed between two similar DNA strands (i.e. two mutant strands in the absence of a wild-type strand), needle macrodissected or high tumor content samples (≥80%tumor) were mixed with a 4ng/ul wild-type control DNA in a 1:1 ratio before being added to the assay. Data output was analysed with the LightCycler®480 Software (Version 1.5.0.39; Roche) and using the Gene Scanning workflow. The melting curves for each amplicon were normalised, and the temperature-shift differences plotted. Samples in which a melt profile differed from that of the wild- type control were selected for Sanger DNA sequence analysis. The HRM product of selected amplicons was PCR amplified, and sequenced using the BigDye Terminator v3.1 (Applied Biosystems) assay. After purification and re-suspension in HiDi Formamide solution (Applied Biosystems) the samples underwent capillary electrophoresis on an ABI3730 sequencer (Applied Biosystems). Mutant sequences were imported into the SequencherTM 4.10.1 software (Gene Codes) for comparison to a wild-type sequence. Mutations were confirmed in a second independent analysis by sequencing of both the forward and reverse strands. S3: Methylation-sensitive HRM analysis Tumor DNA that was used for HRM analysis also underwent a methylation-sensitive high-resolution melting analysis as described previously (6, 7). Briefly, 200ng of DNA was bisulfite converted using the EpiTect Bisulfite Kit (Qiagen) as per the manufacturer’s instructions. Bisulfite converted DNA was added to a mastermix utilising HotStarTaq DNA polymerase containing forward and reverse primers for each assay (Supplementary Table 2) and 5 μmol/L SYTO 9 (Life Technologies, Carlsbad, CA). Assays were performed on the Rotor-Gene 6000 (Corbett, Sydney, Australia), with slight variations in cycling conditions depending on the target (Supplementary Table 3). Temperature was increased by 0.2°C per second during the HRM step. Samples were assayed in duplicate. Appropriate DNA methylation standard controls were included in each assay (in duplicate) for a comparison of extent of methylation during the analysis; 100%, 50%, 25%,10%, 5%, 1% or 0% of fully methylated DNA in unmethylated DNA, as well as a no template control. HRM analysis was performed with the software provided for the Rotor-Gene 6000 (Corbett). For each sample the negative first derivative of the fluorescence over temperature were plotted; the resulting peaks representing the ratio of both methylated and unmethylated DNA in the amplified product. Melting curves were then compared to those of the DNA methylation standard controls, and the extent of methylation estimated based on the shared features between the samples and the controls (6, 7). Gene Forward Primer Reverse Primer Product Size BRCA1 TTGTTGTTTAGCGGTAGTTTTTTGGTT* CAATCGCAATTTTAATTTATCTATAATTCCC* 81bp PALB2 TTTTCGGTTTAGGGTTAATTGGGTT CACCTTTTCCTTCTCCTCACAACTAAA 135bp FANCF ATTGATATGTATTTCGATTAATAGTATTGT ATCCAAATACTACAAAAAAAATTCCATAAA# 149bp Supplementary Table 2: Primers for the methylation-sensitive HRM analysis. * These primers were published in Wong et al. (2011) (9). # This primer was published in Taniguchi et. al (2003) (10). BRCA1 Assay PALB2 Assay FANCF Assay Temperature Time Cycles Temperature Time Cycles Temperature Time Cycles 95°C 15 minutes 1 95°C 15 minutes 1 95°C 15 minutes 1 95°C 10 seconds 50 95°C 20 seconds 50 95°C 20 seconds 50 61°C 10 seconds 62°C 20 seconds 56°C 25 seconds 72°C 20 seconds 72°C 30 seconds 72°C 20 seconds 95°C 1 minute 1 95°C1 minute197°C1 minute1 65-95°C 1 70-95°C 1 65-95°C 1 Supplementary Table 3: Cycling conditions for the methylation-sensitive HRM analysis Twenty-three cases showed evidence of methylation in the analysis, although in two cases this was deemed to be not biologically significant due to the fact that both cases carried a germline BRCA2 pathogenic mutation (Supplementary Table 4). Study ID Estimated % of methylated alleles in Estimated % of tumor material in DNA Other BRCA DNA event 958 <5% MD BRCA2 GL PATH 3202 ~ 50% 50% BRCA2 GL PATH 7492 ~ 50% 60% - 3102 ~15-20% 60% - 3133 ~10% 50% - 5154 ~20% 80% - 1024 ~ 50% 90% - 7329 ~ 50% 40% - 451 ~90% 80% - 1806 ~50% 60-70% - 1978 ~60% MD - 3612 ~100% 90% - 3793 ~50% 80% - 4465 ~50% 70% - 9152 ~50% 85-90% - 6483 ~80% 50-60% - 5349 ~25-30% 90-100% - 4051 ~30% 70% - 2738 ~30% MD - 6621 ~30% 50% BRCA1 UV 3961 ~20% MD - 3843 ~40% 60-80% - 6977 ~30% 70% - Supplementary Table 4: BRCA1 gene promoter methylation detected in the 132 AOCS cases by methylation-sensitive HRM analysis. GL = germline; UV = unclassified variant; MD = macro-dissected GL BRCA1/2 Somatic BRCA1/2 BRCA1 methylation Wild-Type positive positive positive Total Cases 21 8 21 82 Age at Diagnosis (years) Mean 55.01 62.48 55.57 62.79 Standard Deviation 8.71 7.17 9.31 9.76 Primary Site Ovary 81.0% 62.5% 85.7% 74.4% Fallopian Tube 4.8% - 4.8% 4.9% Peritoneum 14.3% 37.5% 9.5% 20.7% Histology Serous 100% 100% 100% 100% FIGO Stage I 4.8% - 14.3% 2.4% II - - 14.3% 2.4% III 81.0% 100% 57.1% 82.9% IV 9.5% - 14.3% 6.1% Not Known 4.8% - - 6.1% Tumour Grade 1 - - - 3.7% 2 4.8% - 23.8% 19.5% 3 81.0% 100% 76.2% 69.5% Not Known 14.3% - - 7.3% Residual Disease Nil macroscopic 28.6% - 33.3% 22.0% ≤ 1cm 52.4% 62.5% 38.1% 32.9% > 1cm 19.0% 37.5% 28.6% 39.0% Not Known/Size Not Known - - - 6.1% Time from diagnosis to disease progression (months) Median 18.10 14.60 23.40 12.30 Range 8.51 – 41.98 10.85 – 20.71 5.00 – 37.31 0.76 – 55.59 Time from diagnosis to death (months) Median 47.1 59.60 62.10 37.10 Range 19.20 – 84.00 20.25 – 71.34 14.56 – 66.94 0.76 – 83.47 Supplementary Table 5: Clinico-pathological features of the AOCS cohort used in this analysis S4: Bioinformatic analyses S4.1 BRCA mutation status and molecular subtypes A sample was considered to have abnormal BRCA status if the sample had a pathogenic germline or somatic mutation in either BRCA gene, or the BRCA1 promoter was methylated.