Benchmarking bioinformatics approaches for tumour mutational burden evaluation from a large cancer panel against whole-exome sequencing.

Tumour mutational burden (TMB) is used to predict response to immunotherapies. Although several groups have proposed calculation methods for TMB, a clear consensus has not yet emerged. In this study, we explored TMB calculation approaches with a 586-gene cancer panel (1.75 Mb) benchmarked to TMB measured by whole-exome sequencing (WES), using 30 samples across a range of tumour types. We explored variant allelic fraction (VAF) cut-offs of 5% and 10%, population database filtering at 0.001, 0.0001 and 0.000025, as well as different combinations of synonymous, insertion/deletion and intronic (splice site) variants, as well as exclusion of hotspot mutations, and examined the effect on TMB correlation. Good correlation (Spearman, range 0.66-0.78) between WES and panel TMB was seen across all methods evaluated. Each method of TMB calculation evaluated showed good positive per cent agreement and negative per cent agreement using 10 mutations/Mb as a cut-off, suggesting that multiple TMB calculation approaches may yield comparable results.

Detection of mitochondrial DNA (mtDNA) mutations.

The maternally inherited mitochondrial DNA (mtDNA) is a circular 16,569bp double stranded DNA that encodes 37 genes, 24 of which (2 rRNAs and 22 tRNAs) are necessary for transcription and translation of 13 polypeptides that are all subunits of respiratory chain. Pathogenic mutations in mtDNA cause respiratory chain dysfunction, and are the underlying defect in an ever-increasing number of mtDNA-related encephalomyopathies with distinct phenotypes. In this chapter, we present an overview of mtDNA mutations and describe the molecular techniques currently employed in our laboratory to detect two types of mtDNA mutations: single-large-scale rearrangements and point mutations.

Microsatellite instability detection using a large next-generation sequencing cancer panel across diverse tumour types.

AIM: Microsatellite instability (MSI), a hallmark of DNA mismatch repair deficiency, is a key molecular biomarker with multiple clinical implications including the selection of patients for immunotherapy, identifying patients who may have Lynch syndrome and predicting prognosis in patients with colorectal tumours. Next-generation sequencing (NGS) provides the opportunity to interrogate large numbers of microsatellite loci concurrently with genomic variants. We sought to develop a method to detect MSI that would not require paired normal tissue and would leverage the sequence data obtained from a broad range of tumours tested using our 467-gene NGS Columbia Combined Cancer Panel (CCCP). METHODS: Altered mononucleotide and dinucleotide microsatellite loci across the CCCP region of interest were evaluated in clinical samples encompassing a diverse range of tumour types. The number of altered loci was used to develop a decision tree classifier model trained on the retrospectively collected cohort of 107 clinical cases sequenced by the CCCP assay. RESULTS: The classifier was able to correctly classify all cases and was then used to analyse a test set of clinical cases (n=112) and was able to correctly predict their MSI status with 100% sensitivity and specificity. Analysis of recurrently altered loci identified alterations in genes involved in DNA repair, signalling and transcriptional regulation pathways, many of which have been implicated in MSI tumours. CONCLUSION: This study highlights the utility of this approach, which should be applicable to laboratories performing similar testing.

Precision Medicine in Children and Young Adults with Hematologic Malignancies and Blood Disorders: The Columbia University Experience.

BACKGROUND: The advent of comprehensive genomic profiling has markedly advanced the understanding of the biology of pediatric hematological malignancies, however, its application to clinical care is still unclear. We present our experience integrating genomic data into the clinical management of children with high-risk hematologic malignancies and blood disorders and describe the broad impact that genomic profiling has in multiple aspects of patient care. METHODS: The Precision in Pediatric Sequencing Program at Columbia University Medical Center instituted prospective clinical next-generation sequencing (NGS) for high-risk malignancies and blood disorders. Testing included cancer whole exome sequencing (WES) of matched tumor-normal samples or targeted sequencing of 467 cancer-associated genes, when sample adequacy was a concern, and tumor transcriptome (RNA-seq). A multidisciplinary molecular tumor board conducted interpretation of results and final tiered reports were transmitted to the electronic medical record according to patient preferences. RESULTS: Sixty-nine samples from 56 patients with high-risk hematologic malignancies and blood disorders were sequenced. Patients carried diagnoses of myeloid malignancy (n = 25), lymphoid malignancy (n = 25), or histiocytic disorder (n = 6). Six patients had only constitutional WES, performed for a suspicion of an inherited predisposition for their disease. For the remaining 50 patients, tumor was sequenced with matched normal tissue when available. The mean number of somatic variants per sample was low across the different disease categories (2.85 variants/sample). Interestingly, a gene fusion was identified by RNA-seq in 58% of samples who had adequate RNA available for testing. Molecular profiling of tumor tissue led to clinically impactful findings in 90% of patients. Forty patients (80%) had at least one targetable gene variant or fusion identified in their tumor tissue; however, only seven received targeted therapy. Importantly, NGS findings contributed to the refinement of diagnosis and prognosis for 34% of patients. Known or likely pathogenic germline alterations were discovered in 24% of patients involving cancer predisposition genes in 12% of cases. CONCLUSION: Incorporating whole exome and transcriptome profiling of tumor and normal tissue into clinical practice is feasible, and the value that comprehensive testing provides extends beyond the ability to target-specific mutations.

Implementation of next generation sequencing into pediatric hematology-oncology practice: moving beyond actionable alterations.

BACKGROUND: Molecular characterization has the potential to advance the management of pediatric cancer and high-risk hematologic disease. The clinical integration of genome sequencing into standard clinical practice has been limited and the potential utility of genome sequencing to identify clinically impactful information beyond targetable alterations has been underestimated. METHODS: The Precision in Pediatric Sequencing (PIPseq) Program at Columbia University Medical Center instituted prospective clinical next generation sequencing (NGS) for pediatric cancer and hematologic disorders at risk for treatment failure. We performed cancer whole exome sequencing (WES) of patient-matched tumor-normal samples and RNA sequencing (RNA-seq) of tumor to identify sequence variants, fusion transcripts, relative gene expression, and copy number variation (CNV). A directed cancer gene panel assay was used when sample adequacy was a concern. Constitutional WES of patients and parents was performed when a constitutionally encoded disease was suspected. Results were initially reviewed by a molecular pathologist and subsequently by a multi-disciplinary molecular tumor board. Clinical reports were issued to the ordering physician and posted to the patient's electronic medical record. RESULTS: NGS was performed on tumor and/or normal tissue from 101 high-risk pediatric patients. Potentially actionable alterations were identified in 38% of patients, of which only 16% subsequently received matched therapy. In an additional 38% of patients, the genomic data provided clinically relevant information of diagnostic, prognostic, or pharmacogenomic significance. RNA-seq was clinically impactful in 37/65 patients (57%) providing diagnostic and/or prognostic information for 17 patients (26%) and identified therapeutic targets in 15 patients (23%). Known or likely pathogenic germline alterations were discovered in 18/90 patients (20%) with 14% having germline alternations in cancer predisposition genes. American College of Medical Genetics (ACMG) secondary findings were identified in six patients. CONCLUSIONS: Our results demonstrate the feasibility of incorporating clinical NGS into pediatric hematology-oncology practice. Beyond the identification of actionable alterations, the ability to avoid ineffective/inappropriate therapies, make a definitive diagnosis, and identify pharmacogenomic modifiers is clinically impactful. Taking a more inclusive view of potential clinical utility, 66% of cases tested through our program had clinically impactful findings and samples interrogated with both WES and RNA-seq resulted in data that impacted clinical decisions in 75% of cases.

Isoform-level brain expression profiling of the spermidine/spermine N1-Acetyltransferase1 (SAT1) gene in major depression and suicide.

Low brain expression of the spermidine/spermine N-1 acetyltransferase (SAT1) gene, the rate-limiting enzyme involved in catabolism of polyamines that mediate the polyamine stress response (PSR), has been reported in depressed suicides. However, it is unknown whether this effect is associated with depression or with suicide and whether all or only specific isoforms expressed by SAT1, such as the primary 171 amino acid protein-encoding transcript (SSAT), or an alternative splice variant (SSATX) that is involved in SAT1 regulated unproductive splicing and transcription (RUST), are involved. We applied next generation sequencing (RNA-seq) to assess gene-level, isoform-level, and exon-level SAT1 expression differences between healthy controls (HC, N = 29), DSM-IV major depressive disorder suicides (MDD-S, N = 21) and MDD non-suicides (MDD, N = 9) in the dorsal lateral prefrontal cortex (Brodmann Area 9, BA9) of medication-free individuals postmortem. Using small RNA-seq, we also examined miRNA species putatively involved in SAT1 post-transcriptional regulation. A DSM-IV diagnosis was made by structured interview. Toxicology and history ruled out recent psychotropic medication. At the gene-level, we found low SAT1 expression in both MDD-S (vs. HC, p = 0.002) and MDD (vs. HC, p = 0.002). At the isoform-level, reductions in MDD-S (vs. HC) were most pronounced in four transcripts including SSAT and SSATX, while reductions in MDD (vs. HC) were pronounced in three transcripts, one of which was reduced in MDD relative to MDD-S (all p < 0.1 FDR corrected). We did not observe evidence for differential exon-usage (i.e. splicing) nor differences in miRNA expression. Results replicate the finding of low SAT1 brain expression in depressed suicides in an independent sample and implicate low SAT1 brain expression in MDD independent of suicide. Low expressions of both SSAT and SATX isoforms suggest that shared transcriptional mechanisms involved in RUST may account for low SAT1 brain expression in depressed suicides. Future studies are required to understand the functions and regulation of SAT1 isoforms, and how they relate to the pathogenesis of MDD and suicide.