Johnson, Sarah H.

AbstractCopy number variation (CNV) is a common form of structural variation detected in human genomes, occurring as both constitutional and somatic events. Cytogenetic techniques like chromosomal microarray (CMA) are widely used in analyzing CNVs. However, CMA techniques cannot resolve the full nature of these structural variations (i.e. the orientation and location of associated breakpoint junctions) and must be combined with other cytogenetic techniques, such as karyotyping or FISH, to do so. This makes the development of a next‐generation sequencing (NGS) approach capable of resolving both CNVs and breakpoint junctions desirable. Mate‐pair sequencing (MPseq) is a NGS technology designed to find large structural rearrangements across the entire genome. Here we present an algorithm capable of performing copy number analysis from mate‐pair sequencing data. The algorithm uses a step‐wise procedure involving normalization, segmentation, and classification of the sequencing data. The segmentation technique combines both read depth and discordant mate‐pair reads to increase the sensitivity and resolution of CNV calls. The method is particularly suited to MPseq, which is designed to detect breakpoint junctions at high resolution. This allows for the classification step to accurately calculate copy number levels at the relatively low read depth of MPseq. Here we compare results for a series of hematological cancer samples that were tested with CMA and MPseq. We demonstrate comparable sensitivity to the state‐of‐the‐art CMA technology, with the benefit of improved breakpoint resolution. The algorithm provides a powerful analytical tool for the analysis of MPseq results in cancer.

Abstract Motivation The ability to produce and analyze whole genome sequencing (WGS) data from samples with structural variations (SV) generated the need to visualize such abnormalities in simplified plots. Conventional two-dimensional representations of WGS data frequently use either circular or linear layouts. There are several diverse advantages regarding both these representations, but their major disadvantage is that they do not use the two-dimensional space very efficiently. We propose a layout, termed the Genome U-Plot, which spreads the chromosomes on a two-dimensional surface and essentially quadruples the spatial resolution. We present the Genome U-Plot for producing clear and intuitive graphs that allows researchers to generate novel insights and hypotheses by visualizing SVs such as deletions, amplifications, and chromoanagenesis events. The main features of the Genome U-Plot are its layered layout, its high spatial resolution and its improved aesthetic qualities. We compare conventional visualization schemas with the Genome U-Plot using visualization metrics such as number of line crossings and crossing angle resolution measures. Based on our metrics, we improve the readability of the resulting graph by at least 2-fold, making apparent important features and making it easy to identify important genomic changes. Results A whole genome visualization tool with high spatial resolution and improved aesthetic qualities. Availability and implementation An implementation and documentation of the Genome U-Plot is publicly available at https://github.com/gaitat/GenomeUPlot. Supplementary information Supplementary data are available at Bioinformatics online.

Abstract TMPRSS2-ERG gene fusions occur in over 50% of prostate cancers, but their impact on clinical outcomes is not well understood. Retention of interstitial genes between TMPRSS2 and ERG has been reported to influence tumor progression in an animal model. In this study, we analyzed the status of TMPRSS2-ERG fusion genes and interstitial genes in tumors from a large cohort of men treated surgically for prostate cancer, associating alterations with biochemical progression. Through whole-genome mate pair sequencing, we mapped and classified rearrangements driving ETS family gene fusions in 133 cases of very low-, low-, intermediate-, and high-risk prostate cancer from radical prostatectomy specimens. TMPRSS2-ERG gene fusions were observed in 44% of cases, and over 90% of these fusions occurred in ERG exons 3 or 4. ERG fusions retaining interstitial sequences occurred more frequently in very low-risk tumors. These tumors also frequently displayed ERG gene fusions involving alternative 5′-partners to TMPRSS2, specifically SLC45A3 and NDRG1 and other ETS family genes, which retained interstitial TMPRSS2/ERG sequences. Lastly, tumors displaying TMPRSS2-ERG fusions that retained interstitial genes were less likely to be associated with biochemical recurrence (P = 0.028). Our results point to more favorable clinical outcomes in patients with ETS family fusion-positive prostate cancers, which retain potential tumor-suppressor genes in the interstitial regions between TMPRSS2 and ERG. Identifying these patients at biopsy might improve patient management, particularly with regard to active surveillance. Cancer Res; 77(22); 6157–67. ©2017 AACR.

Common fragile sites (CFS) are chromosome regions that are prone to form gaps or breaks in response to DNA replication stress. They are often found as hotspots for sister chromatid exchanges, deletions, and amplifications in different cancers. Many of the CFS regions are found to span genes whose genomic sequence is greater than 1 Mb, some of which have been demonstrated to function as important tumor suppressors. CFS regions are also hotspots for human papillomavirus (HPV) integrations in cervical cancer. We used mate‐pair sequencing to examine HPV integration events and chromosomal structural variations in 34 oropharyngeal squamous cell carcinoma (OPSCC). We used endpoint PCR and Sanger sequencing to validate each HPV integration event and found HPV integrations preferentially occurred within CFS regions similar to what is observed in cervical cancer. We also found that many of the chromosomal alterations detected also occurred at or near the cytogenetic location of CFSs. Several large genes were also found to be recurrent targets of rearrangements, independent of HPV integrations, including CSMD1 (2.1Mb), LRP1B (1.9Mb), and LARGE1 (0.7Mb). Sanger sequencing revealed that the nucleotide sequences near to identified junction sites contained repetitive and AT‐rich sequences that were shown to have the potential to form stem‐loop DNA secondary structures that might stall DNA replication fork progression during replication stress. This could then cause increased instability in these regions which could lead to cancer development in human cells. Our findings suggest that CFSs and some specific large genes appear to play important roles in OPSCC. © 2016 Wiley Periodicals, Inc.

Abstract PURPOSE: Previous mutational analyses have failed to identify a predictive or prognostic marker for pulmonary carcinoid tumors (PCT). The heterogeneity of PCTs coupled with a lack of detailed information about the tumor biology, impedes patient stratification into groups based on tumor phenotypes or treatment response. We sought to use high-throughput next-generation sequencing to improve prognostication. METHODS: We stratified 37 patients with resected PCTs as low risk (n = 17) and high risk (n = 20) based on a 5 year long follow up for recurrence. Patients with no recurrence within 5 years from initial treatment were deemed low risk. Macrodissection was followed by genomic DNA isolation 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: Carcinoids commonly harbor rearrangements (median 3.5, range 1-167). Large chromosomal gains or losses were found in 24 of the cases (64.8%). We developed a prognostic score that included the total amount of genetic alterations (deletions and translocations). ROC analysis revealed an AUC of 0.76. High risk patients had a mean score of 20 and low risk patients had a mean score of 5.45 (p value = 0.03, Welch Two Sample t-test). Further analysis of the specific genetic alterations harbored by high and low risk PCTs and correlation with the pathologic and immunohistochemical information is currently underway. CONCLUSION: Large chromosomal rearrangements and aneuploidy portend a poor prognosis for patients with PCTs. Mate Pair sequencing of PCTs provides valuable prognostic information for this highly heterogeneous group of cancers. Citation Format: Konstantinos Leventakos, Aaron S. Mansfield, Stephen J. Murphy, Sarah H. Johnson, Sarah E. Kerr, Marie Christine Aubry, George Vasmatzis, Julian R. Molina. A prognostic model for pulmonary carcinoid tumors based on large chromosomal alterations. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 4513.

Abstract Circulating tumor DNA (ctDNA) provides great potential for the non-invasive clinical assessment of tumor burden. We developed a monitoring protocol using next generation mate-pair sequencing to identify chromosomal rearrangement signatures in primary tumors to follow ctDNA in blood. Primary tumors of 10 patients with ovarian cancer were sequenced, and individualized tumor-specific panels were designed to detect the junctions in ctDNA. A novel algorithm was developed and applied to data generated by quantitative PCR to calculate the percent of ctDNA. Selected junctions, some with potential as therapeutic targets, were detected in pre-surgically drawn blood in eight out of ten patients. Of these eight, three demonstrated the continued presence of circulating tumor DNA post-surgery, and five had undetectable levels, consistent with their documented presence or absence of disease. The detection of somatic junctions derived from ctDNA could be an effective way to monitor ovarian cancer patients for relapse and therapy response. Citation Format: Faye R. Harris, Irina Kovtun, James Smadbeck, Francesco Multinu, Aminah Jatoi, Kimberly R. Kalli, Stephen J. Murphy, Geoffrey C. Halling, Farhad Kosari, Sarah H. Johnson, Andrea Mariani, George Vasmatzis. Quantification of somatic chromosomal rearrangements in circulating cell-free DNA from ovarian cancers. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 440.

Abstract Peripheral T-cell lymphomas (PTCLs) represent a heterogeneous group of T-cell malignancies that generally demonstrate aggressive clinical behavior, often are refractory to standard therapy, and remain significantly understudied. The most common World Health Organization subtype is PTCL, not otherwise specified (NOS), essentially a “wastebasket” category because of inadequate understanding to assign cases to a more specific diagnostic entity. Identification of novel fusion genes has contributed significantly to improving the classification, biologic understanding, and therapeutic targeting of PTCLs. Here, we integrated mate-pair DNA and RNA next-generation sequencing to identify chromosomal rearrangements encoding expressed fusion transcripts in PTCL, NOS. Two of 11 cases had novel fusions involving VAV1, encoding a truncated form of the VAV1 guanine nucleotide exchange factor important in T-cell receptor signaling. Fluorescence in situ hybridization studies identified VAV1 rearrangements in 10 of 148 PTCLs (7%). These were observed exclusively in PTCL, NOS (11%) and anaplastic large cell lymphoma (11%). In vitro, ectopic expression of a VAV1 fusion promoted cell growth and migration in a RAC1-dependent manner. This growth was inhibited by azathioprine, a clinically available RAC1 inhibitor. We also identified novel kinase gene fusions, ITK-FER and IKZF2-ERBB4, as candidate therapeutic targets that show similarities to known recurrent oncogenic ITK-SYK fusions and ERBB4 transcript variants in PTCLs, respectively. Additional novel and potentially clinically relevant fusions also were discovered. Together, these findings identify VAV1 fusions as recurrent and targetable events in PTCLs and highlight the potential for clinical sequencing to guide individualized therapy approaches for this group of aggressive malignancies.

Abstract Many somatic mutations have been detected in pancreatic ductal adenocarcinoma (PDAC), leading to the identification of some key drivers of disease progression, but the involvement of large genomic rearrangements has often been overlooked. In this study, we performed mate pair sequencing (MPseq) on genomic DNA from 24 PDAC tumors, including 15 laser-captured microdissected PDAC and 9 patient-derived xenografts, to identify genome-wide rearrangements. Large genomic rearrangements with intragenic breakpoints altering key regulatory genes involved in PDAC progression were detected in all tumors. SMAD4, ZNF521, and FHIT were among the most frequently hit genes. Conversely, commonly reported genes with copy number gains, including MYC and GATA6, were frequently observed in the absence of direct intragenic breakpoints, suggesting a requirement for sustaining oncogenic function during PDAC progression. Integration of data from MPseq, exome sequencing, and transcriptome analysis of primary PDAC cases identified limited overlap in genes affected by both rearrangements and point mutations. However, significant overlap was observed in major PDAC-associated signaling pathways, with all PDAC exhibiting reduced SMAD4 expression, reduced SMAD-dependent TGFβ signaling, and increased WNT and Hedgehog signaling. The frequent loss of SMAD4 and FHIT due to genomic rearrangements strongly implicates these genes as key drivers of PDAC, thus highlighting the strengths of an integrated genomic and transcriptomic approach for identifying mechanisms underlying disease initiation and progression. Cancer Res; 76(3); 749–61. ©2015 AACR.