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.