Luo, Shujun

Abstract Few methods exist via which to assess treatment response with high fidelity in gliomas. While CSF has previously been used as a source of cfDNA, it remains unknown how cfDNA abundance and sequencing changes throughout treatment and inevitable recurrence in gliomas. We utilized our biobank of longitudinal intracranial CSF specimens to determine the potential utility of CSF cfDNA as a monitoring tool for treatment response in patients with gliomas. Longitudinal CSF specimens were acquired via Ommaya reservoirs (NCT04692337) or ventriculoperitoneal shunts (NCT04692324) in patients with gliomas. cfDNA was extracted from 1-5 mL of CSF, and next-Generation sequencing or low-pass whole genome sequencing (LPWGS) was performed by Predicine, Inc. CSF acquired prior to versus after resection demonstrated an increase in quantified cfDNA (2.97x, range: 1.58-5.26x), consistent with the impact of increased parenchymal disruption and closer contact with CSF. Thereafter, CSF cfDNA continued to decrease throughout chemoradiation and adjuvant treatment and increased at recurrence, including in patients co-enrolled in immunotherapy trials for IL-7 agonists and pembrolizumab. Copy number burden (CNB) increased by the time of recurrence in three patients with known disease progression, and continued to decrease throughout immunotherapy treatment in one patient despite concerns of pseudoprogression. In one patient with a hypermutated glioblastoma, CSF cfDNA revealed over 59 genomic alterations, including MAP2K1, KIT, and PDGFRA, excluding even more variants of unknown significance, each of which had been detected in sequencing of the tissue acquired at resection and subsequently decreased with chemoradiation. In another patient with a known progressing EGFR-amplified GBM, over 200 new variants were identified by the patient’s final sample and EGFR copy number increased from 2 to over 30. In conclusion, cfDNA analysis over a patient’s disease course, including resection and treatment with standard-of-care and experimental therapies, may be of use for disease monitoring and treatment response.

Abstract Next-generation sequencing is now commonly applied to measure genomic modifications, providing the opportunity for a more accurate view of cancer progression. That being said, several challenges including tumor purity, heterogeneity and the broad set of cancer associated cellular functions, increase the difficulty of interpreting cancer genomes and make analyzing them distinct from analysis of germline diseases. The ability to detect low allelic mutations in the broad set of cancer-associated genes is critically important. To address these issues, we developed an extended cancer panel, a targeted enrichment sequencing platform that includes over 1,300 cancer genes and 200 miRNAs that improves coverage in traditionally difficult to sequence NGS regions. In this study, we apply an integrated RNA/DNA analysis approach as well as our extended cancer panel to explore mutations in metastatic renal cell carcinoma (RCC) tissue samples before and after treatment. Our goal with this study was to test the importance of integrating DNA and RNA analysis, which allows for detection of mutational events that cannot be observed by either method alone. With RNA analysis it is possible to detect gene fusion events, expression of small variants (SNVs/indels), and gene expression levels while DNA analysis allows for detection of loss of heterozygosity (LOH), copy number variants (CNVs), and small variants. By integrating both the RNA and DNA analysis, it is possible to identify unique variant classes, such as unexpressed variants, allele specific expression (ASE), effects of CNVs on gene expression levels, and how fusion gene integration sites affect CNVs. In our study, each variant class was interpreted pre- and post-treatment, allowing for optimal analysis of treatment effectiveness. By applying both our RNA/DNA integration approach and our extended cancer panel to investigate metastatic RCC samples, we were able to classify many genomic alterations that drive cancer progression. For example, we found that anywhere from 5% to 30% of small variant mutations in driver genes called in the DNA reside in genes that are unexpressed in the RNA. Likewise, we noted multiple cases of ASE, where high profile variants observed in DNA were entirely absent in the RNA despite high expression levels. Strong correlations between CNVs observed in DNA and gene expression changes found in RNA were also detected. We also observed that about half of the RCC samples we tested were biallelic for VHL mutations, suggesting inherited predisposition combined with second hit somatic mutations, facilitating cancer progression. This analysis demonstrates how combining an integrated RNA/DNA approach with a cancer focused augmented enrichment panel allows for detection of both low allelic representation variants and unique variant classes, both of which are critical for accurate interpretation of cancer samples. Citation Format: Sean M. Boyle, Michael J. Clark, Elena Helman, Alexander S. Parker, Thai Ho, Shujun Luo, Scott Kirk, Parin Sripakdeevong, Mirian Karbelashvili, Deanna M. Church, Michael Snyder, John West, Rich Chen. Integrating RNA/DNA analysis with a comprehensive cancer panel to improve interpretations of stage four metastatic renal cell carcinoma. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 3899. doi:10.1158/1538-7445.AM2015-3899

Abstract Background: Estrogen receptor (ER) expression in breast caner is an important biomarker for targeted therapy and outcome prediction. Nevertheless, the molecular mechanism that controls the phenotypic difference is little known. We used deep sequence technology to profile the transcriptome, gene copy number, and CpG island methylation simultaneously in eight commonly used breast cell lines to develop a model for how these genomic features are integrated in ER+ and ER- breast cancer. Materials and methods: Total mRNA sequence (mRNA-seq), gene copy number (DNA-seq), and genomic CpG island methylation (Methyl-seq) were carried out using the Illumina Genome Analyzer. Sequences were mapped to the human genome (hg18) to obtain digitized gene expression data, DNA copy number in reference to the non-tumor cell line (MCF10A), and methylation status of 21,570 CpG islands; differentially expressed genes between ER+ and ER- cell lines were selected and then correlated with methylation status of these genes’ CpG islands within 5kb of transcript start or copy number changes. The genes that appeared controlled by methylation or copy number changes in cell lines were further evaluated for their expression in the dataset of 129 primary breast tumors. Results: ER+ breast cancer cells had very different gene expression pattern from ER- cancer cells. The two cell line types formed two distinct clusters in unsupervised cluster and 1,873 genes were differentially expressed by moderated t statistics. We identified 149 differentially expressed genes that exhibited differential methylation of one or more CpG islands within 5kb of the 5’ end of the gene and for which mRNA abundance was inversely correlated with CpG island methylation status. Eighty nine of these 149 genes were also differentially expressed in the primary tumor samples and 84 of them were consistent with cell line gene expression and methylation data. The set of 149 genes were significantly enriched in the ER+ tumors by Gene Set Enrichment Analysis. Part of chromosome 8 was deleted in all ER- cells and part of chromosome 17 amplified in all ER+ cells. These loci encoded 30 genes that were overexpressed in ER+ cells; however, only 9 of these genes were overexpressed in ER+ tumors. Many of the methylation or CNA affected genes were functionally significant and some are known to be associated with breast cancer outcome such as GATA3 and LYN. Conclusions: A global pattern of differential CpG island methylation influences expression of a cohort of genes that contribute to the transcriptome landscape of ER+ and ER- breast cancer cells and tumors. The role of gene amplification/deletion in defining the transcriptome landscape of ER+ and ER- cells appears to more modest, although several potentially significant genes appear to be globally regulated by copy number aberrations. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4908. doi:10.1158/1538-7445.AM2011-4908