Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads.

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Nanopore sequencing of 1000 Genomes Project samples to build a comprehensive catalog of human genetic variation.

Less than half of individuals with a suspected Mendelian condition receive a precise molecular diagnosis after comprehensive clinical genetic testing. Improvements in data quality and costs have heightened interest in using long-read sequencing (LRS) to streamline clinical genomic testing, but the absence of control datasets for variant filtering and prioritization has made tertiary analysis of LRS data challenging. To address this, the 1000 Genomes Project ONT Sequencing Consortium aims to generate LRS data from at least 800 of the 1000 Genomes Project samples. Our goal is to use LRS to identify a broader spectrum of variation so we may improve our understanding of normal patterns of human variation. Here, we present data from analysis of the first 100 samples, representing all 5 superpopulations and 19 subpopulations. These samples, sequenced to an average depth of coverage of 37x and sequence read N50 of 54 kbp, have high concordance with previous studies for identifying single nucleotide and indel variants outside of homopolymer regions. Using multiple structural variant (SV) callers, we identify an average of 24,543 high-confidence SVs per genome, including shared and private SVs likely to disrupt gene function as well as pathogenic expansions within disease-associated repeats that were not detected using short reads. Evaluation of methylation signatures revealed expected patterns at known imprinted loci, samples with skewed X-inactivation patterns, and novel differentially methylated regions. All raw sequencing data, processed data, and summary statistics are publicly available, providing a valuable resource for the clinical genetics community to discover pathogenic SVs.

Identification of therapeutic targets in chordoma through comprehensive genomic and transcriptomic analyses.

Chordoma is a rare, orphan cancer arising from embryonal precursors of bone. Surgery and radiotherapy (RT) provide excellent local control, often at the price of significant morbidity because of the structures involved and the need for relatively high doses of RT; however, recurrence remains high. Although our understanding of the genetic changes that occur in chordoma is evolving rapidly, this knowledge has yet to translate into treatments. We performed comprehensive DNA (paired tumor/normal whole-exome and shallow whole-genome) and RNA (tumor whole-transcriptome) next-generation sequencing analyses of archival sacral and clivus chordoma specimens. Incorporation of transcriptomic data enabled the identification of gene overexpression and expressed DNA alterations, thus providing additional support for potential therapeutic targets. In three patients, we identified alterations that may be amenable to off-label FDA-approved treatments for other tumor types. These alterations include FGFR1 overexpression (ponatinib, pazopanib) and copy-number duplication of CDK4 (palbociclib) and ERBB3 (gefitinib). In a third patient, germline DNA demonstrated predicted pathogenic changes in CHEK2 and ATM, which may have predisposed the patient to developing chordoma at a young age and may also be associated with potential sensitivity to PARP inhibitors because of homologous recombination repair deficiency. Last, in the fourth patient, a missense mutation in IGF1R was identified, suggesting potential activity for investigational anti-IGF1R strategies. Our findings demonstrate that chordoma patients present with aberrations in overlapping pathways. These results provide support for targeting the IGF1R/FGFR/EGFR and CDK4/6 pathways as treatment strategies for chordoma patients. This study underscores the value of comprehensive genomic and transcriptomic analysis in the development of rational, individualized treatment plans for chordoma.

Whole Exome Library Construction for Next Generation Sequencing.

Whole exome sequencing (WES) is a DNA sequencing strategy that provides a survey of base substitutions across coding genomic locations and other regions of interest. As the coding portion of the genome encompasses only 1-2% of the entire genome, this approach represents a more cost-effective strategy to detect DNA alterations that may alter protein function, compared to whole genome sequencing. Although the research community has and is currently delineating the functional implications of sequence changes in noncoding regions of the genome, WES is a currently available assay that provides valuable information for both discovery research and precision medicine applications. In this chapter, we present a WES library preparation protocol using the KAPA Hyper Prep Kit with Agilent SureSelect Human All Exon V5+UTR probes that demonstrates high DNA-to-library conversion efficiency for sequencing on the Illumina HiSeq platform.

Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma.

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver BRAF/NRAS/NF1 mutations. In contrast with CM, we observed PAK1 copy gains in 15% of patients, and somatic TERT translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of TERT in ALM tumorigenesis and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.