Genome sequencing reveals a deep intronic splicing ACVRL1 mutation hotspot in Hereditary Haemorrhagic Telangiectasia.

INTRODUCTION: Hereditary haemorrhagic telangiectasia (HHT) is a genetically heterogeneous disorder caused by mutations in the genes ENG, ACVRL1, and SMAD4. Yet the genetic cause remains unknown for some families even after exhaustive exome analysis. We hypothesised that non-coding regions of the known HHT genes may harbour variants that disrupt splicing in these cases. METHODS: DNA from 35 individuals with clinical findings of HHT and 2 healthy controls from 13 families underwent whole genome sequencing. Additionally, 87 unrelated cases suspected to have HHT were evaluated using a custom designed next-generation sequencing panel to capture the coding and non-coding regions of ENG, ACVRL1 and SMAD4. Individuals from both groups had tested negative previously for a mutation in the coding region of known HHT genes. Samples were sequenced on a HiSeq2500 instrument and data were analysed to identify novel and rare variants. RESULTS: Eight cases had a novel non-coding ACVRL1 variant that disrupted splicing. One family had an ACVRL1intron 9:chromosome 3 translocation, the first reported case of a translocation causing HHT. The other seven cases had a variant located within a ~300 bp CT-rich 'hotspot' region of ACVRL1intron 9 that disrupted splicing. CONCLUSIONS: Despite the difficulty of interpreting deep intronic variants, our study highlights the importance of non-coding regions in the disease mechanism of HHT, particularly the CT-rich hotspot region of ACVRL1intron 9. The addition of this region to HHT molecular diagnostic testing algorithms will improve clinical sensitivity.

Clinical utility of a next generation sequencing panel assay for Marfan and Marfan-like syndromes featuring aortopathy.

Aortopathy can be defined as aortic dilation, aneurysm, dissection, and tortuosity. Familial aortopathy may occur secondary to fibrillin-1 (FBN1) mutations in the setting of Marfan syndrome, or may occur as a result of other genetic defects with different, but occasionally overlapping, phenotypes. Because of the phenotypic overlap and genetic heterogeneity of disorders featuring aortopathy, we developed a next generation sequencing (NGS) assay and comparative genomic hybridization (CGH) array to detect mutations in 10 genes that cause thoracic aortic aneurysms (TAAs). Here, we report on the clinical and molecular findings in 175 individuals submitted for aortopathy panel testing at ARUP laboratories. Ten genes associated with heritable aortopathies were targeted using hybridization capture prior to sequencing. NGS results were analyzed, and variants were confirmed using Sanger sequencing. Array CGH was used to detect copy-number variation. Of 175 individuals, 18 had a pathogenic mutation and 32 had a variant of uncertain significance (VUS). Most pathogenic mutations (72%) were identified in FBN1. A novel large SMAD3 duplication and FBN1 deletion were identified. Over half who had TAAs or other aortic involvement tested negative for a mutation, suggesting that additional aortopathy genes exist. We anticipate that the clinical sensitivity of at least 10.3% will rise with VUS reclassification and as additional genes are identified and included in the panel. The aortopathy NGS panel aids in the timely molecular diagnosis of individuals with disorders featuring aortopathy and guides proper treatment.

Phospholipid isotope tracing suggests ß-catenin-driven suppression of phosphatidylcholine metabolism in hepatocellular carcinoma.

Activating mutations in the CTNNB1 gene encoding ß-catenin are among the most frequently observed oncogenic alterations in hepatocellular carcinoma (HCC). Profound alterations in lipid metabolism, including increases in fatty acid oxidation and transformation of the phospholipidome, occur in HCC with CTNNB1 mutations, but it is unclear what mechanisms give rise to these changes. We employed untargeted lipidomics and targeted isotope tracing to measure phospholipid synthesis activity in an inducible human liver cell line expressing mutant ß-catenin, as well as in transgenic zebrafish with activated ß-catenin-driven HCC. In both models, activated ß-catenin expression was associated with large changes in the lipidome including conserved increases in acylcarnitines and ceramides and decreases in triglycerides. Lipid isotope tracing analysis in human cells revealed a reduction in phosphatidylcholine (PC) production rates as assayed by choline incorporation. We developed lipid isotope tracing analysis for zebrafish tumors and observed reductions in phosphatidylcholine synthesis by both the CDP-choline and PEMT pathways. The observed changes in the ß-catenin-driven HCC phospholipidome suggest that zebrafish can recapitulate conserved features of HCC lipid metabolism and may serve as a model for identifying future HCC-specific lipid metabolic targets.

Phospholipid isotope tracing reveals ß-catenin-driven suppression of phosphatidylcholine metabolism in hepatocellular carcinoma.

Background and Aims: Activating mutations in the CTNNB1 gene encoding ß-catenin are among the most frequently observed oncogenic alterations in hepatocellular carcinoma (HCC). HCC with CTNNB1 mutations show profound alterations in lipid metabolism including increases in fatty acid oxidation and transformation of the phospholipidome, but it is unclear how these changes arise and whether they contribute to the oncogenic program in HCC. Methods: We employed untargeted lipidomics and targeted isotope tracing to quantify phospholipid production fluxes in an inducible human liver cell line expressing mutant ß-catenin, as well as in transgenic zebrafish with activated ß-catenin-driven HCC. Results: In both models, activated ß-catenin expression was associated with large changes in the lipidome including conserved increases in acylcarnitines and ceramides and decreases in triglycerides. Lipid flux analysis in human cells revealed a large reduction in phosphatidylcholine (PC) production rates as assayed by choline tracer incorporation. We developed isotope tracing lipid flux analysis for zebrafish and observed similar reductions in phosphatidylcholine synthesis flux accomplished by sex-specific mechanisms. Conclusions: The integration of isotope tracing with lipid abundances highlights specific lipid class transformations downstream of ß-catenin signaling in HCC and suggests future HCC-specific lipid metabolic targets.

Intratumor Heterogeneity in Hepatocellular Carcinoma: Challenges and Opportunities.

Hepatocellular carcinoma (HCC) represents a leading cause of cancer-related death, but it remains difficult to treat. Intratumor genetic and phenotypic heterogeneity are inherent properties of breast, skin, lung, prostate, and brain tumors, and intratumor heterogeneity (ITH) helps define prognosis and therapeutic response in these cancers. Several recent studies estimate that ITH is inherent to HCC and attribute the clinical intractability of HCC to this heterogeneity. In this review, we examine the evidence for genomic, phenotypic, and tumor microenvironment ITH in HCC, with a focus on two of the top molecular drivers of HCC: ß-catenin (CTNNB1) and Telomerase reverse transcriptase (TERT). We discuss the influence of ITH on HCC diagnosis, prognosis, and therapy, while highlighting the gaps in knowledge and possible future directions.

Large-Scale Identification of Clonal Hematopoiesis and Mutations Recurrent in Blood Cancers.

Clonal hematopoiesis of indeterminate potential (CHIP) is characterized by detectable hematopoietic-associated gene mutations in a person without evidence of hematologic malignancy. We sought to identify additional cancer-presenting mutations useable for CHIP detection by performing a data mining analysis of 48 somatic mutation studies reporting mutations at diagnoses of 7,430 adult and pediatric patients with hematologic malignancies. Following extraction of 20,141 protein-altering mutations, we identified 434 significantly recurrent mutation hotspots, 364 of which occurred at loci confidently assessable for CHIP. We then performed an additional large-scale analysis of whole exome sequencing data from 4,538 persons belonging to three non-cancer cohorts for clonal mutations. We found the combined cohort prevalence of CHIP with mutations identical to those reported at blood cancer mutation hotspots to be 1.8%, and that some of these CHIP mutations occurred in children. Our findings may help to improve CHIP detection and pre-cancer surveillance for both children and adults.

NF1 Somatic Mutation in Dystrophic Scoliosis.

Scoliosis is a common manifestation of neurofibromatosis type 1, causing significant morbidity. The etiology of dystrophic scoliosis in neurofibromatosis type 1 is not fully understood and therapies are lacking. Somatic mutations in NF1 have been shown in tibial pseudarthrosis providing rationale for similar processes in neurofibromatosis type 1-associated dystrophic scoliosis. Spinal samples from surgical procedures with matched peripheral blood of two individuals with neurofibromatosis type 1 and dystrophic scoliosis were obtained and DNA extracted. Next generation sequencing of various spinal sections as well as the germline/blood sample were performed using a RASopathy gene panel (includes the NF1 gene). Variants were compared between the spinal tissue samples and the germline data. In addition, the next generation sequencing allele frequency data were used to detect somatic loss of heterozygosity. All samples had a detected potentially inactivating NF1 germline mutation. Both individuals demonstrated an allelic imbalance inclusive of NF1 in the next generation sequencing data. In addition, for the same two individuals, there was an increase in the % variant reads for the germline mutation in some of the surgical spinal samples corresponding to the allelic imbalance. Contra analysis did not show any deletion in Chromosome 17 next generation sequencing data. Microarray analysis verified somatic copy neutral loss of heterozygosity for these two individuals for the majority of the chromosome 17 q-arm, inclusive of the NF1 gene. These results suggest that the cause of dystrophic scoliosis is multifactorial and that a somatic NF1 mutation contributes to the etiology.

Utilization of Whole-Exome Next-Generation Sequencing Variant Read Frequency for Detection of Lesion-Specific, Somatic Loss of Heterozygosity in a Neurofibromatosis Type 1 Cohort with Tibial Pseudarthrosis.

A subset of neurofibromatosis type 1 patients develop tibial dysplasia, which can lead to pseudarthrosis. The tissue from the tibial pseudarthrosis region commonly has a somatic second hit in NF1: single-nucleotide variants, small deletions, or loss of heterozygosity (LOH). We used exome next-generation sequencing (NGS) variant frequency data (allelic imbalance analysis) to detect somatic LOH in pseudarthrosis tissue from three individuals with clinically and diagnostically confirmed neurofibromatosis type 1, and verified the results with microarray. The variant files were parsed and plotted using python scripts, and the NGS variant frequencies between the affected tissue and blood sample were compared. Individuals without somatic single-nucleotide variants or small insertions/deletions were tested for somatic LOH using the NGS variant allele frequencies. One individual's NGS data indicated no LOH in chromosome 17. The other two individuals demonstrated somatic LOH inclusive of NF1: one had an LOH region of approximately one million bases and Contra (NGS copy number program) indicated a somatic deletion and the other individual had LOH for most of chromosome 17q and Contra indicated no copy number change (microarray data verified this sample as copy neutral somatic LOH). Both LOH and copy number variation detected by NGS data correlated with microarray data, demonstrating the somatic LOH second hit can be detected directly from the NGS data.

Hereditary hemorrhagic telangiectasia: genetics and molecular diagnostics in a new era.

Hereditary hemorrhagic telangiectasia (HHT) is a vascular dysplasia characterized by telangiectases and arteriovenous malformations (AVMs) in particular locations described in consensus clinical diagnostic criteria published in 2000. Two genes in the transforming growth factor-beta (TGF-ß) signaling pathway, ENG and ACVRL1, were discovered almost two decades ago, and mutations in these genes have been reported to cause up to 85% of HHT. In our experience, approximately 96% of individuals with HHT have a mutation in these two genes, when published (Curaçao) diagnostic criteria for HHT are strictly applied. More recently, two additional genes in the same pathway, SMAD4 and GDF2, have been identified in a much smaller number of patients with a similar or overlapping phenotype to HHT. Yet families still exist with compelling evidence of a hereditary telangiectasia disorder, but no identifiable mutation in a known gene. Recent availability of whole exome and genome testing has created new opportunities to facilitate gene discovery, identify genetic modifiers to explain clinical variability, and potentially define an increased spectrum of hereditary telangiectasia disorders. An expanded approach to molecular diagnostics for inherited telangiectasia disorders that incorporates a multi-gene next generation sequencing (NGS) HHT panel is proposed.

Next-generation sequencing of custom amplicons to improve coverage of HaloPlex multigene panels.

Next-generation sequencing (NGS) of multigene panels performed for genetic clinical diagnostics requires 100% coverage of all targeted genes. In the genetic diagnostics laboratory, coverage gaps are typically filled with Sanger sequencing after NGS data are collected and analyzed. Libraries prepared using the hybridization-based custom capture HaloPlex method are covered at ~98% and include gaps in coverage because of the location of the restriction enzyme sites used for fragmentation and differences in the designed and actual library insert size. We describe a method for improving the coverage of HaloPlex libraries by generating a set of amplicons spanning known low-coverage regions that are pooled, indexed by sample, and sequenced together with the HaloPlex libraries. This approach reduces the number of post-NGS Sanger sequencing reactions required and complements any NGS library preparation method when complete gene coverage is necessary.