Eckloff, Bruce W.

Abstract Next Generation Sequencing (NGS) technologies provide rapid genomic analyses of single nucleotide variants, RNA expression and DNA copy number. Application of these technologies to material isolated from formalin fixed paraffin embedded (FFPE) tissue and even degraded frozen material could provide powerful replication samples but remains challenging. We tested the nanoString platform to validate deep sequence analysis of gene expression and DNA copy number in degraded and FFPE material. Firstly, RNA from the Universal Human Reference RNA and a breast cancer cell line (MDA-MB-436) was artificially degraded to different degrees (RIN 1.2-6.8). We used the nanoString platform to simultaneously measure RNA expression across 226 genes in each degraded sample and the corresponding undegraded RNA. Secondly we isolated RNA and DNA from matched fresh frozen and FFPE tissues from nine breast cancer patients (3 HER2+/ER+/PR+, 2 HER2+/ER-/PR-, 2 HER2+/ER+/PR-, 2 HER2-/ER+/PR+) using the nanoString platform to compare expression and copy number across 226 and 86 genes respectively. Finally, we correlated expression and copy number data generated by nanoString with Illumina transcriptome and whole genome sequencing (WGS). NanoString log2 expression fold-change between all artificially degraded samples and their undegraded counterpart showed extremely high correlation (r2>0.91). NanoString DNA copy number between matched fresh-frozen and FFPE showed a high degree of correlation (r2=0.71). All gene amplifications with copy number ≥ 5 in DNA from fresh-frozen material (N=9) were successfully identified in DNA from FFPE material. We also observed good correlation of gene expression between whole transcriptome sequencing and the nanoString platform (r2 0.59 - 0.72) in FFPE and artificially degraded material and for DNA copy number between WGS and nanoString in DNA isolated from cancer cell lines (r2=0.96). The nanoString platform provides reliable data from highly degraded and FFPE material and correlates with sequence analysis of both expression and copy number from NGS platforms demonstrating potential for large-scale replication studies in FFPE material. Citation Format: Nadine Norton, Edith A. Perez, Yan W. Asmann, Jennifer M. Carr, Brian M. Necela, Jennifer M. Kachergus, Jin Jen, Bruce W. Eckloff, E Aubrey Thompson. Analysis of gene expression and copy number variation in breast tumors using both sequencing and hybridization-based platforms. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2005. doi:10.1158/1538-7445.AM2013-2005

Abstract Previous studies have analyzed patterns of transcription, transcription factor (TF) binding or mapped nucleosome occupancy across the genome. These suggest that the three aspects are genetically connected but the cause and effect relationships are still unknown. For example, physiologic TF binding studies involve many TFs, consequently, it is difficult to assign nucleosome reorganization to the binding site occupancy of any particular TF. Therefore, several aspects remain unclear: does TF binding influence nucleosome (re)organizations locally or impact the chromatin landscape at a more global level; are all or only a fraction of TF binding a result of reorganization in nucleosome occupancy and do all TF binding and associated changes in nucleosome occupancy result in altered gene expression? With these in mind, following characterization of two states (before and after induction of a single TF of choice) we determined: (i) genomic binding sites of the TF, (ii) promoter nucleosome occupancy and (iii) transcriptome profiles. Results demonstrated that promoter-proximal TF binding influenced expression of the target gene when it was coupled to nucleosome repositioning at or close to its binding site in most cases. In contrast, only in few cases change in target gene expression was found when TF binding occurred without local nucleosome reorganization.

Chromosomal translocations leading to expression of abnormal fusion proteins play a major role in the pathogenesis of various hematologic malignancies. The recent development of high‐throughput, “deep” sequencing has allowed discovery of novel translocations leading to a rapid increase in understanding these diseases. Translocations involving the anaplastic lymphoma kinase (ALK) gene leading to ALK fusion proteins originally were discovered in anaplastic large cell lymphomas (ALCLs). Among ALCLs, NPM1‐ALK fusions are most common and lead to nuclear localization of the fusion protein. Here, we present a 50‐year‐old male with ALCL demonstrating cytoplasmic ALK immunoreactivity only, suggesting the presence of a non‐NPM1 fusion partner. We performed deep RNA sequencing of tumor tissue from this patient and identified a novel transcript fusing Exon 6 of TRAF1 to Exon 20 of ALK. The TRAF1‐ALK fusion transcript was confirmed at the mRNA level by Sanger sequencing and the fusion protein was visualized by Western blot. The discovery of this TRAF1‐ALK fusion expands the diversity of known ALK fusion partners and highlights the power of deep sequencing for fusion transcript discovery. © 2013 Wiley Periodicals, Inc.

Abstract Background: Colorectal cancer (CRC) in densely affected families without Lynch Syndrome may be due to mutations in undiscovered genetic loci. Familial linkage analyses have yielded disparate results; the use of exome sequencing in coding regions may identify novel segregating variants. Methods: We completed exome sequencing on 40 affected cases from 16 multicase pedigrees to identify novel loci. Variants shared among all sequenced cases within each family were identified and filtered to exclude common variants and single-nucleotide variants (SNV) predicted to be benign. Results: We identified 32 nonsense or splice-site SNVs, 375 missense SNVs, 1,394 synonymous or noncoding SNVs, and 50 indels in the 16 families. Of particular interest are two validated and replicated missense variants in CENPE and KIF23, which are both located within previously reported CRC linkage regions, on chromosomes 1 and 15, respectively. Conclusions: Whole-exome sequencing identified DNA variants in multiple genes. Additional sequencing of these genes in additional samples will further elucidate the role of variants in these regions in CRC susceptibility. Impact: Exome sequencing of familial CRC cases can identify novel rare variants that may influence disease risk. Cancer Epidemiol Biomarkers Prev; 22(7); 1239–51. ©2013 AACR.

Background— The primary role of natriuretic peptide receptor-3 (NPR3) or NPR-C is in the clearance of natriuretic peptides that play an important role in modulating intravascular volume and vascular tone. Genetic variation in NPR3 has been associated with variation in blood pressure and obesity. Despite the importance of NPR3, sequence variation in the gene has not been addressed using DNA from different ethnic populations. We set out to identify and functionally characterize genetic variation in NPR3 in 3 ethnic groups. Methods and Results— DNA samples from 96 European American, 96 African American, and 96 Han Chinese American healthy subjects were used to resequence NPR3 exons, splice junctions, and flanking regions. We identified 105 polymorphisms, 50 of which were novel, including 8 nonsynonymous single-nucleotide polymorphisms, 7 were novel. Expression constructs were created for the nonsynonymous single-nucleotide polymorphisms. HEK293 cells were transfected with constructs for wild type and variant allozymes; and recombinant proteins were measured by quantitative Western blot analysis. The most significant change in NPR3 protein was observed for the Arg146 variant allozyme, with 20% of wild-type protein, primarily because of autophagy-dependent degradation. NPR3 structural modeling confirmed that the Arg146 variant protein was not compatible with wild-type conformation and could result in protein misfolding or instability. Conclusions— Multiple novel NPR3 genetic polymorphisms were identified in 3 ethnic groups. The Arg146 allozyme displayed a significant decrease in protein quantity because of degradation mediated predominantly by autophagy. This genetic variation could have a significant effect on the metabolism of natriuretic peptides with potential clinical implications.

Abstract Chromosome rearrangements that result in gene fusions have important roles in the initial steps of tumorigenesis especially in leukemias and lymphomas, but the biological and clinical impact of gene fusions in the common solid tumor has been less appreciated. We previously generated RNAseq data from 12 oropharyngeal tumors along with matched normal tissue from the same patients. We utilized and compared two bioinformatic pipelines, SnowShoes-FTD and FusionHunter to identify two novel fusion events involving of MGST3-ZMAT5 and MS4A7-DKFZp586P0123. These fusions were validated utilizing. Reverse transcriptase polymerase chain reaction (RT-PCR) and Sanger sequencing validated. In the first translocation the 5′ UTR of MGST3 was fused with exon 1 of ZMAT5, and in the second translocation the exon 2 of MS4A7 was fused with 3′ UTR of DKFZp586P0123. These precise fusions were not present in the matched normal tissue, however, they were only found in the single tumor in which they were identified. Given a growing and powerful set of high-throughput massively parallel sequencing technologies, the future for genome-wide analysis of chromosome rearrangement appears open in many directions. Ultimately, we hope that this integrative methodology will elucidate key aspects of tumor biology as well as generate novel targets for cancer diagnostics and therapy. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 3100. doi:1538-7445.AM2012-3100

Abstract Peripheral T-cell lymphomas (PTCLs) are aggressive malignancies of mature T lymphocytes with 5-year overall survival rates of only ∼ 35%. Improvement in outcomes has been stymied by poor understanding of the genetics and molecular pathogenesis of PTCL, with a resulting paucity of molecular targets for therapy. We developed bioinformatic tools to identify chromosomal rearrangements using genome-wide, next-generation sequencing analysis of mate-pair DNA libraries and applied these tools to 16 PTCL patient tissue samples and 6 PTCL cell lines. Thirteen recurrent abnormalities were identified, of which 5 involved p53-related genes (TP53, TP63, CDKN2A, WWOX, and ANKRD11). Among these abnormalities were novel TP63 rearrangements encoding fusion proteins homologous to ΔNp63, a dominant-negative p63 isoform that inhibits the p53 pathway. TP63 rearrangements were seen in 11 (5.8%) of 190 PTCLs and were associated with inferior overall survival; they also were detected in 2 (1.2%) of 164 diffuse large B-cell lymphomas. As TP53 mutations are rare in PTCL compared with other malignancies, our findings suggest that a constellation of alternate genetic abnormalities may contribute to disruption of p53-associated tumor suppressor function in PTCL.

Abstract Several platforms for genome sequencing and methods for exon enrichment are now available and continue to evolve quickly. In this study, we compared the performance of whole exome sequencing using two instruments, Illumina's GAIIx and HiSeq 2000, comparing single versus duplicate runs. Eight runs of eight samples each were performed with genomic DNA derived from whole blood and sample libraries prepared with the Agilent SureSelect capture array. One run was conducted in duplicate on the GAIIx, and, for some analyses, the combined data were used. We used TREAT (Targeted RE-sequencing and Annotation Tool), developed in-house for data analysis, including sequence alignment (MAQ, BWA), local re-alignment (GATK), variant calling (utilizing MAQ), annotation (SIFT and Seattle Seq), and visualization. Overall, duplicate runs on the GAIIx platform (compared to a single run) increased the number of total reads by approximately two-fold: from ∼66 M to ∼130 M. For both single and duplicate runs, ∼90% of reads mapped to the reference sequence, and ∼55% of reads mapped on-target. Of importance, the percent coverage of the target region increased substantially with duplicate runs; ∼72% of the target region was covered at 40-fold or more when run in duplicate compared to only ∼52% when run once. At 30x, 20x, and 10x, there was a 1.3- (62% v. 78%), 1.2- (72% v. 84%), and 1.1-fold (84% v. 91%) increase respectively, in the percent coverage of the target region for samples run twice. A similar number of filtered, on-target SNPs per sample (∼24 K) was found for both single and duplicate run analyses. However, a 1.3-fold increase in the number of on-target indels was seen in the duplicate run (∼1,200) compared to the single run (∼900). When the combined data for samples run twice on a GAIIx were compared to the samples run using a HiSeq, the total number of reads were similar (∼130 M), although several of the samples had substantially more reads on the HiSeq platform. Samples run on a HiSeq had an increase in the percent coverage of the target region at 40x (72% v. 81%), 30x (78% v. 86%), 20x (84% v. 91%) and 10x (91% v. 95%) compared to the two-run approach on a GAIIx. A similar number of filtered on-target SNPs (∼24 K/sample) and indels (∼1200) was found with the HiSeq compared to two runs on a GAIIx. In summary, our results demonstrate increases in the total number of reads and in overall coverage with the HiSeq 2000 instrument. The total number of SNPs and indel that mapped on target for the version 1 Agilent SureSelect capture array is ∼25 K and 1,200 per sample, respectively. As capture is of variable efficiency for individual runs, higher average coverage of the target region is necessary for sufficient coverage of poorly captured regions. Additional analyses comparing Agilent SureSelect versions 1 and 2 arrays are underway. 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 2742. doi:10.1158/1538-7445.AM2011-2742

Abstract There are a number of human genes which span very large chromosomal regions. Forty genes span greater than 1 megabase and 200 are 500 Kb or larger. We and others have demonstrated that many of the largest genes are derived from within regions of profound genomic instability known as common fragile sites (CFSs). We have been studying oropharyngeal cancers which are derived from the base of tongue, tonsil and larynx. While the traditional etiologic factors underlying the development of this aggressive cancer is a history of smoking and drinking there are an increasing number of younger patients developing this cancer without these risk factors. The most likely cause of these cancers is the presence of human papillomavirus which we detect in over 90% of oropharyngeal cancers derived from the Mayo Clinic. In an attempt to better understand the molecular alterations that underlie oropharyngeal cancers arising from different risk factors we utilized the power of Next Generation sequencing on the Illumina GA IIx to compare the transcriptomes of 12 oropharyngeal cancers to patient matched normal tissue. We first observe a division of patients into two groups of oropharyngeal cancers: those caused by mutations and other alterations induced by a history of smoking and drinking and those associated with HPV infection caused by viral-induced genomic instability. We then examined the expression of the 200 largest human genes in the oropharyngeal tumors and find that each tumor has inactivated the expression of at least 20 different very large genes. There is a distinct pattern of inactivation of these genes in different oropharyngeal cancers depending upon the underlying etiology. Some of the genes most frequently inactivated include Dlg2 (the human homolog of the large discs tumor suppressor gene from drosophila and a mutation target in a number of other cancers); ERBB4 (Her4 and a member of the EGFR family), BAI3 (a brain specific angiogenesis inhibitor that is also mutated in multiple cancers), and ANK2 (mutated in long QT patients). Each of these genes may be functioning as tumor suppressors as has already been demonstrated for the two large CFS genes, FHIT and WWOX. The inactivation of so many of these genes in each oropharyngeal tumor may have a profound effect upon the phenotype of the resulting tumor. When we examine the expression of the large genes most frequently inactivated in each of the two groups of oropharyngeal cancers in panels of breast or lung cancer samples we find that we do not detect significant loss of expression of these genes. We believe that genomic instability promotes alterations in the unstable CFS regions but that there is selective inactivation of different large CFS genes in different cancers. The signature of which specific genes are inactivated may have phenotypic and clinical significance which we are just beginning to investigate. 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 3885. doi:10.1158/1538-7445.AM2011-3885

Abstract The genetics of peripheral T-cell lymphomas are poorly understood. The most well-characterized abnormalities are translocations involving ALK, occurring in approximately half of anaplastic large cell lymphomas (ALCLs). To gain insight into the genetics of ALCLs lacking ALK translocations, we combined mate-pair DNA library construction, massively parallel (“Next Generation”) sequencing, and a novel bioinformatic algorithm. We identified a balanced translocation disrupting the DUSP22 phosphatase gene on 6p25.3 and adjoining the FRA7H fragile site on 7q32.3 in a systemic ALK-negative ALCL. Using fluorescence in situ hybridization, we demonstrated that the t(6;7)(p25.3;q32.3) was recurrent in ALK-negative ALCLs. Furthermore, t(6;7)(p25.3;q32.3) was associated with down-regulation of DUSP22 and up-regulation of MIR29 microRNAs on 7q32.3. These findings represent the first recurrent translocation reported in ALK-negative ALCL and highlight the utility of massively parallel genomic sequencing to discover novel translocations in lymphoma and other cancers.