Expanding the clinical phenotype of FGFR1 internal tandem duplication.

Closed spinal dysraphism (SD) is a type of neural tube defect originating during early embryonic development whereby the neural tissue of the spinal defect remains covered by skin, often coinciding with markers of cutaneous stigmata. It is hypothesized that these events are caused by multifactorial processes, including genetic and environmental causes. We present an infant with a unique congenital midline lesion associated with a closed SD. Through comprehensive molecular profiling of the intraspinal lesion and contiguous skin lesion, an internal tandem duplication (ITD) of the kinase domain of the fibroblast growth factor receptor 1 (FGFR1) gene was found. This ITD variant is somatic mosaic in nature as supported by a diminished variant allele frequency in the lesional tissue and by its absence in peripheral blood. FGFR1 ITD results in constitutive activation of the receptor tyrosine kinase to promote cell growth, differentiation, and survival through RAS/MAPK signaling. Identification of FGFR1 ITD outside of central nervous system tumors is exceedingly rare, and this report broadens the phenotypic spectrum of somatic mosaic FGFR1-related disease.

Long-read whole genome sequencing reveals HOXD13 alterations in synpolydactyly.

Synpolydactyly 1, also called syndactyly type II (SDTY2), is a genetic limb malformation characterized by polydactyly with syndactyly involving the webbing of the third and fourth fingers, and the fourth and fifth toes. It is caused by heterozygous alterations in HOXD13 with incomplete penetrance and phenotypic variability. In our study, a five-generation family with an SPD phenotype was enrolled in our Rare Disease Genomics Protocol. A comprehensive examination of three generations using Illumina short-read whole-genome sequencing (WGS) did not identify any causative variants. Subsequent WGS using Pacific Biosciences (PacBio) long-read HiFi Circular Consensus Sequencing (CCS) revealed a heterozygous 27-bp duplication in the polyalanine tract of HOXD13. Sanger sequencing of all available family members confirmed that the variant segregates with affected individuals. Reanalysis of an unrelated family with a similar SPD phenotype uncovered a 21-bp (7-alanine) duplication in the same region of HOXD13. Although ExpansionHunter identified these events in most individuals in a retrospective analysis, low sequence coverage due to high GC content in the HOXD13 polyalanine tract makes detection of these events challenging. Our findings highlight the value of long-read WGS in elucidating the molecular etiology of congenital limb malformation disorders.

Molecular classification of a complex structural rearrangement of the RB1 locus in an infant with sporadic, isolated, intracranial, sellar region retinoblastoma.

Retinoblastoma is a childhood cancer of the retina involving germline or somatic alterations of the RB Transcriptional Corepressor 1 gene, RB1. Rare cases of sellar-suprasellar region retinoblastoma without evidence of ocular or pineal tumors have been described. A nine-month-old male presented with a sellar-suprasellar region mass. Histopathology showed an embryonal tumor with focal Flexner-Wintersteiner-like rosettes and loss of retinoblastoma protein (RB1) expression by immunohistochemistry. DNA array-based methylation profiling confidently classified the tumor as pineoblastoma group A/intracranial retinoblastoma. The patient was subsequently enrolled on an institutional translational cancer research protocol and underwent comprehensive molecular profiling, including paired tumor/normal exome and genome sequencing and RNA-sequencing of the tumor. Additionally, Pacific Biosciences (PacBio) Single Molecule Real Time (SMRT) sequencing was performed from comparator normal and disease-involved tissue to resolve complex structural variations. RNA-sequencing revealed multiple fusions clustered within 13q14.1-q21.3, including a novel in-frame fusion of RB1-SIAH3 predicted to prematurely truncate the RB1 protein. SMRT sequencing revealed a complex structural rearrangement spanning 13q14.11-q31.3, including two somatic structural variants within intron 17 of RB1. These events corresponded to the RB1-SIAH3 fusion and a novel RB1 rearrangement expected to correlate with the complete absence of RB1 protein expression. Comprehensive molecular analysis, including DNA array-based methylation profiling and sequencing-based methodologies, were critical for classification and understanding the complex mechanism of RB1 inactivation in this diagnostically challenging tumor.

Characterization and mitigation of fragmentation enzyme-induced dual stranded artifacts.

High-throughput short-read sequencing relies on fragmented DNA for optimal sampling of input nucleic acid. Several vendors now offer proprietary enzyme cocktails as a cheaper and more streamlined method of fragmentation when compared to acoustic shearing. We have discovered that these enzymes induce the formation of library molecules containing regions of nearby DNA from opposite strands. Sequencing reads derived from these molecules can lead to artifact-derived variant calls appearing at variant allele frequencies <5%. We present Fragmentation Artifact Detection and Elimination (FADE), software to remove these artifacts from mapped reads and mitigate artifact-related effects on downstream analysis. We find that the artifacts principally affect downstream analyses that are sensitive to a 1-3% artifact bias in the sequencing reads, such as targeted resequencing and rare variant discovery.

Transcriptionally Active Androgen Receptor Splice Variants Promote Hepatocellular Carcinoma Progression.

Owing to the marked sexual dimorphism of hepatocellular carcinoma (HCC), sex hormone receptor signaling has been implicated in numerous aspects of liver cancer pathogenesis. We sought to reconcile the clear contribution of androgen receptor (AR) activity that has been established in preclinical models of HCC with the clinical failure of AR antagonists in patients with advanced HCC by evaluating potential resistance mechanisms to AR-targeted therapy. The AR locus was interrogated for resistance-causing genomic modifications using publicly available primary HCC datasets (1,019 samples). Analysis of HCC tumor and cell line RNA-seq data revealed enriched expression of constitutively active, treatment-refractory AR splice variants (AR-SV). HCC cell lines expressed C-terminal-truncated AR-SV; 28 primary HCC samples abundantly expressed AR-SV. Low molecular weight AR species were nuclear localized and constitutively active. Furthermore, AR/AR-SV signaling promoted AR-mediated HCC cell progression and conferred resistance to AR antagonists. Ligand-dependent and -independent AR signaling mediated HCC epithelial-to-mesenchymal transition by regulating the transcription factor SLUG. These data suggest that AR-SV expression in HCC drives HCC progression and resistance to traditional AR antagonists. Novel therapeutic approaches that successfully target AR-SVs may be therapeutically beneficial for HCC. SIGNIFICANCE: Treatment-refractory, constitutively active androgen receptor splice variants promote hepatocellular carcinoma progression by regulating the epithelial-to-mesenchymal transition pathway.

Challenges in proteogenomics: a comparison of analysis methods with the case study of the DREAM proteogenomics sub-challenge.

BACKGROUND: Proteomic measurements, which closely reflect phenotypes, provide insights into gene expression regulations and mechanisms underlying altered phenotypes. Further, integration of data on proteome and transcriptome levels can validate gene signatures associated with a phenotype. However, proteomic data is not as abundant as genomic data, and it is thus beneficial to use genomic features to predict protein abundances when matching proteomic samples or measurements within samples are lacking. RESULTS: We evaluate and compare four data-driven models for prediction of proteomic data from mRNA measured in breast and ovarian cancers using the 2017 DREAM Proteogenomics Challenge data. Our results show that Bayesian network, random forests, LASSO, and fuzzy logic approaches can predict protein abundance levels with median ground truth-predicted correlation values between 0.2 and 0.5. However, the most accurately predicted proteins differ considerably between approaches. CONCLUSIONS: In addition to benchmarking aforementioned machine learning approaches for predicting protein levels from transcript levels, we discuss challenges and potential solutions in state-of-the-art proteogenomic analyses.

BRCA1/2 Functional Loss Defines a Targetable Subset in Leiomyosarcoma.

BACKGROUND: Soft-tissue sarcomas (STS) describe a heterogeneous group of mesenchymal tumors with limited treatment options. Targeted therapies exist for BRCA1/2 gene alterations, but their prevalence and role have not been fully described in STS. Here, we present the largest effort to characterize the frequency of homologous recombination (HR) DNA repair pathway alterations in STS subtypes and highlight the unique nature of leiomyosarcoma (LMS). MATERIALS AND METHODS: DNA sequencing data were analyzed for HR pathway alterations for 1,236 patients with STS. DNA sequencing data from an additional 1,312 patients were used to confirm the prevalence of HR pathway alterations in LMS. Four uterine LMS (uLMS) patients with functional BRCA2 loss were evaluated for response to poly (ADP-ribose) polymerase (PARP) inhibition. RESULTS: In an unselected STS study population, BRCA2 alterations were identified in 15 (1%) patients, and homozygous BRCA2 loss was detected in 9 (<1%). However, subset analysis revealed that these BRCA2 alterations were concentrated in uLMS as compared with any other STS subtype. Notably, 10% of uLMS tumors had a BRCA2 alteration. We further report that PARP inhibitors had demonstrated durable clinical benefit in four uLMS patients with BRCA2 loss. CONCLUSION: HR pathway alterations are rare in most STS. However, we identify uLMS to be enriched for BRCA2 loss and report the positive outcomes of a series of patients treated with PARP inhibitors. Our data suggest that patients with uLMS should be considered for somatic BRCA2 profiling. Prospective trials are necessary to confirm the efficacy of PARP inhibition in uLMS. IMPLICATIONS FOR PRACTICE: Soft-tissue sarcomas are a highly morbid, diverse set of tumors with limited treatment options. This study identifies an increased prevalence of functional BRCA1/2 loss in patients with uterine leiomyosarcoma (uLMS). It also presents four patients with uLMS and BRCA2 loss who achieved durable clinical benefit from poly (ADP-ribose) polymerase inhibition. These data suggest that patients with uLMS in particular should be screened for BRCA1/2 alterations and may benefit from treatment targeted to these alterations.

Validation of a Targeted RNA Sequencing Assay for Kinase Fusion Detection in Solid Tumors.

Kinase gene fusions are important drivers of oncogenic transformation and can be inhibited with targeted therapies. Clinical grade diagnostics using RNA sequencing to detect gene rearrangements in solid tumors are limited, and the few that are available require prior knowledge of fusion break points. To address this, we have analytically validated a targeted RNA sequencing assay (OSU-SpARKFuse) for fusion detection that interrogates complete transcripts from 93 kinase and transcription factor genes. From a total of 74 positive and 36 negative control samples, OSU-SpARKFuse had 93.3% sensitivity and 100% specificity for fusion detection. Assessment of repeatability and reproducibility revealed 96.3% and 94.4% concordance between intrarun and interrun technical replicates, respectively. Application of this assay on prospective patient samples uncovered OLFM4 as a novel RET fusion partner in a small-bowel cancer and led to the discovery of a KLK2-FGFR2 fusion in a patient with prostate cancer who subsequently underwent treatment with a pan-fibroblast growth factor receptor inhibitor. Beyond fusion detection, OSU-SpARKFuse has built-in capabilities for discovery research, including gene expression analysis, detection of single-nucleotide variants, and identification of alternative splicing events.

Landscape of Microsatellite Instability Across 39 Cancer Types.

PURPOSE: Microsatellite instability (MSI) is a pattern of hypermutation that occurs at genomic microsatellites and is caused by defects in the mismatch repair system. Mismatch repair deficiency that leads to MSI has been well described in several types of human cancer, most frequently in colorectal, endometrial, and gastric adenocarcinomas. MSI is known to be both predictive and prognostic, especially in colorectal cancer; however, current clinical guidelines only recommend MSI testing for colorectal and endometrial cancers. Therefore, less is known about the prevalence and extent of MSI among other types of cancer. METHODS: Using our recently published MSI-calling software, MANTIS, we analyzed whole-exome data from 11,139 tumor-normal pairs from The Cancer Genome Atlas and Therapeutically Applicable Research to Generate Effective Treatments projects and external data sources across 39 cancer types. Within a subset of these cancer types, we assessed mutation burden, mutational signatures, and somatic variants associated with MSI. RESULTS: We identified MSI in 3.8% of all cancers assessed-present in 27 of tumor types-most notably adrenocortical carcinoma (ACC), cervical cancer (CESC), and mesothelioma, in which MSI has not yet been well described. In addition, MSI-high ACC and CESC tumors were observed to have a higher average mutational burden than microsatellite-stable ACC and CESC tumors. CONCLUSION: We provide evidence of as-yet-unappreciated MSI in several types of cancer. These findings support an expanded role for clinical MSI testing across multiple cancer types as patients with MSI-positive tumors are predicted to benefit from novel immunotherapies in clinical trials.

Performance evaluation for rapid detection of pan-cancer microsatellite instability with MANTIS.

In current clinical practice, microsatellite instability (MSI) and mismatch repair deficiency detection is performed with MSI-PCR and immunohistochemistry. Recent research has produced several computational tools for MSI detection with next-generation sequencing (NGS) data; however a comprehensive analysis of computational methods has not yet been performed. In this study, we introduce a new MSI detection tool, MANTIS, and demonstrate its favorable performance compared to the previously published tools mSINGS and MSISensor. We evaluated 458 normal-tumor sample pairs across six cancer subtypes, testing classification performance on variable numbers of target loci ranging from 10 to 2539. All three computational methods were found to be accurate, with MANTIS exhibiting the highest accuracy with 98.91% of samples from all six diseases classified correctly. MANTIS displayed superior performance among the three tools, having the highest overall sensitivity (MANTIS 97.18%, MSISensor 96.48%, mSINGS 76.06%) and specificity (MANTIS 99.68%, mSINGS 99.68%, MSISensor 98.73%) across six cancer types, even with loci panels of varying size. Additionally, MANTIS also had the lowest resource consumption (<1% of the space and <7% of the memory required by mSINGS) and fastest running times (49.6% and 8.7% of the running times of MSISensor and mSINGS, respectively). This study highlights the potential utility of MANTIS in classifying samples by MSI-status, allowing its incorporation into existing NGS pipelines.

Cancer Driver Log (CanDL): Catalog of Potentially Actionable Cancer Mutations.

Massively parallel sequencing technologies have enabled characterization of genomic alterations across multiple tumor types. Efforts have focused on identifying driver mutations because they represent potential targets for therapy. However, because of the presence of driver and passenger mutations, it is often challenging to assign the clinical relevance of specific mutations observed in patients. Currently, there are multiple databases and tools that provide in silico assessment for potential drivers; however, there is no comprehensive resource for mutations with functional characterization. Therefore, we created an expert-curated database of potentially actionable driver mutations for molecular pathologists to facilitate annotation of cancer genomic testing. We reviewed scientific literature to identify variants that have been functionally characterized in vitro or in vivo as driver mutations. We obtained the chromosome location and all possible nucleotide positions for each amino acid change and uploaded them to the Cancer Driver Log (CanDL) database with associated literature reference indicating functional driver evidence. In addition to a simple interface, the database allows users to download all or selected genes as a comma-separated values file for incorporation into their own analysis pipeline. Furthermore, the database includes a mechanism for third-party contributions to support updates for novel driver mutations. Overall, this freely available database will facilitate rapid annotation of cancer genomic testing in molecular pathology laboratories for mutations.