Caroline Conway

Using next-generation sequencing for high resolution multiplex analysis of copy number variation from nanogram quantities of DNA from formalin-fixed paraffin-embedded specimens

The use of next-generation sequencing technologies to produce genomic copy number data has recently been described. Most approaches, however, reply on optimal starting DNA, and are therefore unsuitable for the analysis of formalin-fixed paraffin-embedded (FFPE) samples, which largely precludes the analysis of many tumour series. We have sought to challenge the limits of this technique with regards to quality and quantity of starting material and the depth of sequencing required. We confirm that the technique can be used to interrogate DNA from cell lines, fresh frozen material and FFPE samples to assess copy number variation. We show that as little as 5 ng of DNA is needed to generate a copy number karyogram, and follow this up with data from a series of FFPE biopsies and surgical samples. We have used various levels of sample multiplexing to demonstrate the adjustable resolution of the methodology, depending on the number of samples and available resources. We also demonstrate reproducibility by use of replicate samples and comparison with microarray-based comparative genomic hybridization (aCGH) and digital PCR. This technique can be valuable in both the analysis of routine diagnostic samples and in examining large repositories of fixed archival material.

Gene Expression Profiling of Paraffin-Embedded Primary Melanoma Using the DASL Assay Identifies Increased Osteopontin Expression as Predictive of Reduced Relapse-Free Survival

Purpose: Gene expression studies in melanoma have been few because tumors are small and cryopreservation is rarely possible. The purpose of this study was to evaluate the Illumina DASL Array Human Cancer Panel for gene expression studies in formalin-fixed melanoma primary tumors and to identify prognostic biomarkers. Experimental Design: Primary tumors from two studies were sampled using a tissue microarray needle. Study 1: 254 tumors from a melanoma cohort recruited from 2000 to 2006. Study 2: 218 tumors from a case-control study of patients undergoing sentinel node biopsy. Results: RNA was obtained from 76 of blocks; 1.4 of samples failed analysis (transcripts from <250 of the 502 genes on the DASL chip detected). Increasing age of the block and increased melanin in the tumor were associated with reduced number of genes detected. The gene whose expression was most differentially expressed in association with relapse-free survival in study 1 was osteopontin (SPP1; P = 2.11 106) and supportive evidence for this was obtained in study 2 used as a validation set (P = 0.006; unadjusted data). Osteopontin level in study 1 remained a significant predictor of relapse-free survival when data were adjusted for age, sex, tumor site, and histologic predictors of relapse. Genes whose expression correlated most strongly with osteopontin were PBX1, BIRC5 (survivin), and HLF. Conclusion: Expression data were obtained from 74 of primary melanomas and provided confirmatory evidence that osteopontin expression is a prognostic biomarker. These results suggest that predictive biomarker studies may be possible using stored blocks from mature clinical trials. (Clin Cancer Res 2009;15(22):693946)

Patterns of Expression of DNA Repair Genes and Relapse From Melanoma

Purpose: To use gene expression profiling of formalin-fixed primary melanoma samples to detect expression patterns that are predictive of relapse and response to chemotherapy. Experimental Design: Gene expression profiles were identified in samples from two studies (472 tumors). Gene expression data for 502 cancer-related genes from these studies were combined for analysis. Results: Increased expression of DNA repair genes most strongly predicted relapse and was associated with thicker tumors. Increased expression of RAD51 was the most predictive of relapse-free survival in unadjusted analysis (hazard ratio, 2.98; P = 8.80 × 10−6). RAD52 (hazard ratio, 4.73; P = 0.0004) and TOP2A (hazard ratio, 3.06; P = 0.009) were independent predictors of relapse-free survival in multivariable analysis. These associations persisted when the analysis was further adjusted for demographic and histologic features of prognostic importance (RAD52 P = 0.01; TOP2A P = 0.02). Using principal component analysis, expression of DNA repair genes was summarized into one variable. Genes whose expression correlated with this variable were predominantly associated with the cell cycle and DNA repair. In 42 patients treated with chemotherapy, DNA repair gene expression was greater in tumors from patients who progressed on treatment. Further data supportive of a role for increased expression of DNA repair genes as predictive biomarkers are reported, which were generated using multiplex PCR. Conclusions: Overexpression of DNA repair genes (predominantly those involved in double-strand break repair) was associated with relapse. These data support the hypothesis that melanoma progression requires maintenance of genetic stability and give insight into mechanisms of melanoma drug resistance and potential therapies. Clin Cancer Res; 16(21); 5211–21. ©2010 AACR.

Next-generation sequencing for simultaneous determination of human papillomavirus load, subtype, and associated genomic copy number changes in tumors.

Human papillomavirus (HPV) infection in cases of squamous cell carcinoma of the oropharynx is a powerful predictive and prognostic biomarker. We describe how the use of next-generation sequencing can provide a novel method for the detection of HPV in DNA isolated from formalin-fixed paraffin-embedded tissues. Using this methodology in a cohort of 44 head and neck tumors, we identified the samples that contained HPV sequences, the viral subtype involved, and a direct readout of viral load. Specificity of HPV detection by sequencing compared to traditional detection methods using either PCR or p16 immunohistochemistry was 100%. Sensitivity was 50% when either compared to PCR [confidence interval (CI) = 29% to 71%] or 75% when compared to p16 (CI = 47% to 91%). In addition, we demonstrate the ability of next-generation sequencing to detect other HPV subtypes that would not have been detected by traditional methods, and we demonstrated the ability to apply this method to any tumor and any virus in a panel of eight human cancer cell lines. This methodology also provides a tumor genomic copy number karyogram, and in the samples analyzed here, a lower level of chromosome instability was detected in HPV-positive tumors compared to HPV-negative tumors, as observed in previous studies. Thus, the use of next-generation sequencing for the detection of HPV provides a multiplicity of data with clinical significance in a single test.

Deletion at chromosome arm 9p in relation to BRAF/NRAS mutations and prognostic significance for primary melanoma.

We report an investigation of gene dosage at 9p21.3 and mutations in BRAF and NRAS, as predictors of relapse and histological markers of poor melanoma prognosis. Formalin‐fixed primary melanomas from 74 relapsed and 42 nonrelapsed patients were sequenced for common BRAF and NRAS mutations (N = 71 results) and gene dosage at 9p21.3 including the genes CDKN2A (which encodes CDKN2A and P14ARF), CDKN2B (CDKN2B), and MTAP was measured using multiplexed ligation‐dependant probe amplification (MLPA), (N = 75 results). BRAF/NRAS mutations were detected in 77% of relapsers and 58% of nonrelapsers (Fishers exact P = 0.17), and did not predict ulceration or mitotic rate. There was no relationship between BRAF/NRAS mutations and gene dosage at 9p21.3. Reduced gene dosage at MTAP showed a borderline association with BRAF mutation (P = 0.04) and reduced gene dosage at the interferon gene cluster was borderline associated with wild type NRAS (P = 0.05). Reduced gene dosage in the CDKN2A regions coding for CDKN2A was associated with an increased risk of relapse (P = 0.03). Reduced gene dosage across 9p21.3 was associated with increased tumor thickness, mitotic rate, and ulceration (P = 0.02, 0.02, and 0.002, respectively), specifically in coding regions impacting on CDKN2B and P14ARF and CDKN2A. Loss at MTAP (P = 0.05) and the interferon gene cluster (P = 0.03) on 9p21 was also associated with tumor ulceration. There was no association between reduced gene dosage at 9p21.3 and subtype or site of tumor. This study presents supportive evidence that CDKN2B, P14ARF, and CDKN2A may all play a tumor suppressor role in melanoma progression.

Melanoma sentinel node biopsy and prediction models for relapse and overall survival

Background:To optimise predictive models for sentinal node biopsy (SNB) positivity, relapse and survival, using clinico-pathological characteristics and osteopontin gene expression in primary melanomas.Methods:A comparison of the clinico-pathological characteristics of SNB positive and negative cases was carried out in 561 melanoma patients. In 199 patients, gene expression in formalin-fixed primary tumours was studied using Illuminas DASL assay. A cross validation approach was used to test prognostic predictive models and receiver operating characteristic curves were produced.Results:Independent predictors of SNB positivity were Breslow thickness, mitotic count and tumour site. Osteopontin expression best predicted SNB positivity (P=2.4 × 10−7), remaining significant in multivariable analysis. Osteopontin expression, combined with thickness, mitotic count and site, gave the best area under the curve (AUC) to predict SNB positivity (72.6%). Independent predictors of relapse-free survival were SNB status, thickness, site, ulceration and vessel invasion, whereas only SNB status and thickness predicted overall survival. Using clinico-pathological features (thickness, mitotic count, ulceration, vessel invasion, site, age and sex) gave a better AUC to predict relapse (71.0%) and survival (70.0%) than SNB status alone (57.0, 55.0%). In patients with gene expression data, the SNB status combined with the clinico-pathological features produced the best prediction of relapse (72.7%) and survival (69.0%), which was not increased further with osteopontin expression (72.7, 68.0%).Conclusion:Use of these models should be tested in other data sets in order to improve predictive and prognostic data for patients.

The clinicopathological and gene expression patterns associated with ulceration of primary melanoma

Ulceration of primary melanomas is associated with poor prognosis yet is reported to predict benefit from adjuvant interferon. To better understand the biological processes involved, clinicopathological factors associated with ulceration were determined in 1804 patients. From this cohort, 348 primary tumor blocks were sampled to generate gene expression data using a 502‐gene cancer panel and 195 blocks were used for immunohistochemistry to detect macrophage infiltration and vessel density. Gene expression results were validated using a whole genome array in two independent sample sets. Ulceration of primary melanomas was associated with more proliferative tumors, tumor vessel invasion, and increased microvessel density. Infiltration of tumors with greater number of macrophages and gene expression pathways associated with wound healing and up‐regulation of pro‐inflammatory cytokines suggests that ulceration is associated with tumor‐related inflammation. The relative benefit from interferon reported in patients with ulcerated tumors may reflect modification of signaling pathways involved in inflammation.

Clinicopathologic features of V600E and V600K melanoma--letter.

The presence of a BRAF mutation in a melanoma tumor predicts response to BRAF inhibitors; however, the biological characteristics of tumors with different BRAF mutations have not been investigated until recently. Menzies and colleagues reported that the rarer V600K BRAF mutations were found more

A computational index derived from whole-genome copy number analysis is a novel tool for prognosis in early stage lung squamous cell carcinoma.

Squamous cell carcinoma of the lung is remarkable for the extent to which the same chromosomal abnormalities are detected in individual tumours. We have used next generation sequencing at low coverage to produce high resolution copy number karyograms of a series of 89 non-small cell lung tumours specifically of the squamous cell subtype. Because this methodology is able to create karyograms from formalin-fixed paraffin-embedded material, we were able to use archival stored samples for which survival data were available and correlate frequently occurring copy number changes with disease outcome. No single region of genomic change showed significant correlation with survival. However, adopting a whole-genome approach, we devised an algorithm that relates to total genomic damage, specifically the relative ratios of copy number states across the genome. This algorithm generated a novel index, which is an independent prognostic indicator in early stage squamous cell carcinoma of the lung.

The clonal relationships between pre‐cancer and cancer revealed by ultra‐deep sequencing

The study of the relationships between pre‐cancer and cancer and identification of early driver mutations is becoming increasingly important as the value of molecular markers of early disease and personalised drug targets is recognized, especially now the extent of clonal heterogeneity in fully invasive disease is being realized. It has been assumed that pre‐cancerous lesions exhibit a fairly passive progression to invasive disease; the degree to which they, too, are heterogeneous is unknown. We performed ultra‐deep sequencing of thousands of selected mutations, together with copy number analysis, from multiple, matched pre‐invasive lesions, primary tumours and metastases from five patients with oral cancer, some with multiple primary tumours presenting either synchronously or metachronously, totalling 75 samples. This allowed the clonal relationships between the samples to be observed for each patient. We expose for the first time the unexpected variety and complexity of the relationships between this group of oral dysplasias and their associated carcinomas and, ultimately, the diversity of processes by which tumours are initiated, spread and metastasize. Instead of a series of genomic precursors of their adjacent invasive disease, we have shown dysplasia to be a distinct dynamic entity, refuting the belief that pre‐cancer and invasive tumours with a close spatial relationship always have linearly related genomes. We show that oral pre‐cancer exhibits considerable subclonal heterogeneity in its own right, that mutational changes in pre‐cancer do not predict the onset of invasion, and that the genomic pathway to invasion is neither unified nor predictable. Sequence data from this study have been deposited in the European Nucleotide Archive, Accession No. PRJEB6588. Copyright