Johnson, Sarah H.

Abstract Introduction: Identification of causal DNA copy number variations (CNVs) is critical to the understanding of cancer development and progression, and is important clinically for accurate prognosis and efficient treatment. T-cell lymphomas are replete with amplifications and deletions that result in over- and under-expression of genes. Array comparative genomic hybridization (aCGH) has been the technology of choice to detect CNVs but recently, it has been proposed that mate-pair Next Generation Sequencing (mpNGS) can produce similar results. Here, we report a comparison between the two techniques. Methods: We have developed an algorithm based on extreme value theory to identify CNVs across the genome at high resolution from mate-pair NGS data. The algorithm is based on plotting the distribution of read-counts in windows across the genome, and detecting deletions and gains by finding the windows that differ significantly from the distribution of the normal part of the genome. We compared this algorithm with aCGH data to detect CNVs. We ran five cases of T-cell lymphoma with both techniques and compared the results using concordance statistics. Results: With 50kb lower cutoff for a CNV, mpNGS detected 98% of the cumulative deleted genome and 94% of the amplified/gained genome that were called by aCGH. Most of the discordant CNVs that were detected by only aCGH were in regions of repetitive areas where mapping fails. Conversely, aCGH detected 75% of the cumulative deleted genome and 74% of the amplified/gained genome that were called by mpNGS. Conclusion: A detailed comparison of these techniques showed that both have pros and cons. In general, mpNGS can detect most CNVs detected by aCGH, as long as they are more than 50KB. Additionally, mpNGS offers confirmatory data, particularly for deletions, based on detection of aberrant mate pairs that map to the rejoined ends of deleted regional. Finally, the mate-pair technology can also detect balanced translocations and can also predict fusion genes. Citation Format: George Vasmatzis, Andrew L. Feldman, Sarah H. Johnson, Erik C. Thorland, Rafael Fonseca, Esteban Braggio, Troy J. Gliem. Copy number detection using genomics technologies: A comparison between aCGH and NGS. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3572. doi:10.1158/1538-7445.AM2014-3572

Abstract Background: Fusion genes, transcripts, and proteins are key drivers of hematopoietic oncogenesis and can serve as effective therapeutic targets. Among peripheral T-cell lymphomas (PTCLs), little is known about fusion genes except for those involving ALK. Discovery of fusion genes may provide therapeutic targets for PTCLs, which have poor outcomes with conventional therapy. We have developed methods for detecting candidate fusions using (1) mate-pair DNA sequencing (MPseq), an approach that spans the entire genome at a fraction of the cost of traditional whole-genome sequencing, and (2) RNA sequencing (RNAseq). Each technique has strengths, limitations, and false discovery rates. Here, we integrated results of both methods to discover and orthogonally validate expressed fusion genes in PTCLs. Methods: DNA and RNA were extracted from 49 frozen, OCT-embedded clinical PTCL specimens, including 25 anaplastic large cell lymphomas; 16 PTCLs, not otherwise specified; and 8 other PTCLs. MPseq and RNAseq were performed using Illumina library protocols and a HiSeq 2000. Detection of fusion genes and transcripts followed our published methods (Blood 2012;120:2280 and Nucleic Acids Res 2011;39:e100, respectively). Functional annotation of genes involved in fusions was performed using DAVID. Results: Of 2818 candidate chromosomal rearrangements identified by MPseq and 90 candidate fusion transcripts identified by RNAseq, 45 events were present in both datasets (mean/case, 0.9; range, 0-7). Fusions were inter-chromosomal in 16 and intra-chromosomal in 29. Annotation clusters with the highest enrichment scores related to SH2 domain-containing proteins and tyrosine kinases (including FER and MAP2K3). Seven genes were recurrently involved in fusions, including known fusion partners such as ALK and previously unreported fusion partners such as PPP1R8, a phosphatase inhibitor important for the activity of Polycomb group proteins such as EZH2. Detection of candidate events only by MPseq could be explained in part by breakpoints mapping to non-genic regions (79% of events). Detection only by RNAseq could be explained in part by unannotated transcripts spanning neighboring named genes (median intra-chromosomal distance, 0.09 MB vs. 12.7 MB for intra-chromosomal candidates identified by MPseq; p<0.0001, Wilcoxon test). Conclusions: Integrated DNA/RNA sequencing is an effective, clinically feasible tool for simultaneous discovery and orthogonal validation of expressed fusion genes. This approach yielded fusions in PTCLs involving genes with key functions in lymphocytes, including novel and recurrent fusions, some of which might be targetable with compounds currently in clinical trial. Candidate intra-chromosomal fusion transcripts not identified by MPseq may represent unrecognized isoforms or the results of novel read-through mechanisms, and merit further study. Citation Format: Andrew L. Feldman, George Vasmatzis, Yan W. Asmann, Sarah H. Johnson, Bruce W. Eckloff, Sumit Middha, Jaime I. Davila, Paul J. Kurtin, Brian K. Link, Stephen M. Ansell, James R. Cerhan. Integrated DNA/RNA sequencing for discovery and orthogonal validation of expressed fusion genes in peripheral T-cell lymphomas. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5585. doi:10.1158/1538-7445.AM2014-5585

Purpose Distinguishing independent primary tumors from intrapulmonary metastases in non–small-cell carcinoma remains a clinical dilemma with significant clinical implications. Using next-generation DNA sequencing, we developed a chromosomal rearrangement–based approach to differentiate multiple primary tumors from metastasis. Methods Tumor specimens from patients with known independent primary tumors and metastatic lesions were used for lineage test development, which was then applied to multifocal tumors. Laser capture microdissection was performed separately for each tumor. Genomic DNA was isolated using direct in situ whole-genome amplification methodology, and next-generation sequencing was performed using an Illumina mate-pair library protocol. Sequence reads were mapped to the human genome, and primers spanning the fusion junctions were used for validation polymerase chain reaction. Results A total of 41 tumor samples were sequenced (33 adenocarcinomas [ADs] and eight squamous cell carcinomas [SQCCs]), with a range of three to 276 breakpoints per tumor identified. Lung tumors predicted to be independent primary tumors based on different histologic subtype did not share any genomic rearrangements. In patients with lung primary tumors and paired distant metastases, shared rearrangements were identified in all tumor pairs, emphasizing the patient specificity of identified breakpoints. Multifocal AD and SQCC samples were reviewed independently by two pulmonary pathologists. Concordance between histology and genomic data occurred in the majority of samples. Discrepant tumor samples were resolved by genome sequencing. Conclusion A diagnostic lineage test based on genomic rearrangements from mate-pair sequencing demonstrates promise for distinguishing independent primary from metastatic disease in lung cancer.

AbstractThe development of adenocarcinoma of the lung is believed to proceed from in situ disease (adenocarcinoma in situ, AIS) to minimally invasive disease with prominent lepidic growth (minimally invasive adenocarcinoma, MIA), then to fully invasive adenocarcinoma (AD), but direct evidence for this model has been lacking. Because some lung adenocarcinomas show prominent lepidic growth (AD-L), we designed a study to address the lineage relationship between the lepidic (noninvasive) component (L) and the adjacent nonlepidic growth component representing invasive disease within individual tumors. Lineage relationships were evaluated by next-generation DNA sequencing to define large genomic rearrangements in microdissected tissue specimens collected by laser capture. We found a strong lineage relationship between the majority of adjacent lepidic and invasive components, supporting a putative AIS–AD transition. Notably, many rearrangements were detected in the less aggressive lepidic component, although the invasive component exhibited an overall higher rate of genomic rearrangement. Furthermore, a significant number of genomic rearrangements were present in histologically normal lung adjacent to tumor, but not in host germline DNA, suggesting field defects restricted to zonal regions near a tumor. Our results offer a perspective on the genetic pathogenesis underlying adenocarcinoma development and its clinical management. Cancer Res; 74(11); 3157–67. ©2014 AACR.

Abstract Summary: Mate pair library sequencing is an effective and economical method for detecting genomic structural variants and chromosomal abnormalities. Unfortunately, the mapping and alignment of mate-pair read pairs to a reference genome is a challenging and time-consuming process for most next-generation sequencing alignment programs. Large insert sizes, introduction of library preparation protocol artifacts (biotin junction reads, paired-end read contamination, chimeras, etc.) and presence of structural variant breakpoints within reads increase mapping and alignment complexity. We describe an algorithm that is up to 20 times faster and 25% more accurate than popular next-generation sequencing alignment programs when processing mate pair sequencing. Availability: http://bioinformaticstools.mayo.edu/research/bima/ Contact: vasm@mayo.edu or vasmatzis.george@mayo.edu Supplementary Information: Supplementary data are available at Bioinformatics online.

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.

We used high resolution mate‐pair sequencing (HRMPS) in 15 patients with primary myelofibrosis (PMF): eight with normal karyotype and seven with PMF‐characteristic cytogenetic abnormalities, including der(6)t(1;6)(q21–23;p21.3) (n = 4), der(7)t(1;7)(q10;p10) (n = 2), del(20)(q11.2q13.3) (n = 3), and complex karyotype (n = 1). We describe seven novel deletions/translocations in five patients (including two with normal karyotype) whose breakpoints were PCR‐validated and involved MACROD2, CACNA2D4, TET2, SGMS2, LRBA, SH3D19, INTS3, FOP (CHTOP), SCLT1, and PHF17. Deletions with breakpoints involving MACROD2 (lysine deacetylase; 20p12.1) were recurrent and found in two of the 15 study patients. A novel fusion transcript was found in one of the study patients (INTS3‐CHTOP), and also in an additional non‐study patient with PMF. In two patients with der(6)t(1;6)(q21–23;p21.3), we were able to map the precise translocation breakpoints, which involved KCNN3 and GUSBP2 in one case and HYDIN2 in another. This study demonstrates the utility of HRMPS in uncovering submicroscopic deletions/translocations/fusions, and precise mapping of breakpoints in those with overt cytogenetic abnormalities. The overall results confirm the genetic heterogeneity of PMF, given the low frequency of recurrent specific abnormalities, identified by this screening strategy. Currently, we are pursuing the pathogenetic relevance of some of the aforementioned findings. Am. J. Hematol. 88:741–746, 2013. © 2013 Wiley Periodicals, Inc.

Abstract Gleason score 7 (GS7) prostate cancer [tumors with both Gleason patterns 3 (GP3) and 4 (GP4)] portends a significantly more aggressive tumor than Gleason score 6 (GS6). It is, therefore, critical to understand the molecular relationship of adjacent GP3 and GP4 tumor cell populations and relate molecular abnormalities to disease progression. To decipher molecular relatedness, we used laser capture microdissection (LCM) and whole-genome amplification (WGA) to separately collect and amplify DNA from adjacent GP3 and GP4 cell populations from 14 cases of GS7 prostate cancer. We then carried out massively parallel mate-pair next generation sequencing (NGS) to examine the landscape of large chromosomal alterations. We identified four to 115 DNA breakpoints in GP3 and 17 to 480 in GP4. Our findings indicate that while GP3 and GP4 from the same tumor each possess unique breakpoints, they also share identical ones, indicating a common origin. Approximately 300 chromosomal breakpoints were localized to the regions affected in at least two tumors, whereas more than 3,000 were unique within the set of 14 tumors. TMPRSS2–ERG was the most recurrent rearrangement present in eight cases, in both GP3 and GP4. PTEN rearrangements were found in five of eight TMPRSS2–ERG fusion–positive cases in both GP3 and GP4. Hierarchical clustering analysis revealed that GP3 has greater breakpoint similarity to its partner GP4 compared with GP3 from different patients. We show evidence that LCM, WGA, and NGS of adjacent tumor regions provide an important tool in deciphering lineage relationships and discovering chromosomal alterations associated with tumor progression. Cancer Res; 73(11); 3275–84. ©2013 AACR.

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 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.