Comprehensive genetic and phenotype analysis of 95 individuals with mosaic tuberous sclerosis complex.

Tuberous sclerosis complex (TSC) is a neurogenetic disorder due to loss-of-function TSC1 or TSC2 variants, characterized by tumors affecting multiple organs, including skin, brain, heart, lung, and kidney. Mosaicism for TSC1 or TSC2 variants occurs in 10%-15% of individuals diagnosed with TSC. Here, we report comprehensive characterization of TSC mosaicism by using massively parallel sequencing (MPS) of 330 TSC samples from a variety of tissues and fluids from a cohort of 95 individuals with mosaic TSC. TSC1 variants in individuals with mosaic TSC are much less common (9%) than in germline TSC overall (26%) (p < 0.0001). The mosaic variant allele frequency (VAF) is significantly higher in TSC1 than in TSC2, in both blood and saliva (median VAF: TSC1, 4.91%; TSC2, 1.93%; p = 0.036) and facial angiofibromas (median VAF: TSC1, 7.7%; TSC2 3.7%; p = 0.004), while the number of TSC clinical features in individuals with TSC1 and TSC2 mosaicism was similar. The distribution of mosaic variants across TSC1 and TSC2 is similar to that for pathogenic germline variants in general TSC. The systemic mosaic variant was not present in blood in 14 of 76 (18%) individuals with TSC, highlighting the value of analysis of multiple samples from each individual. A detailed comparison revealed that nearly all TSC clinical features are less common in individuals with mosaic versus germline TSC. A large number of previously unreported TSC1 and TSC2 variants, including intronic and large rearrangements (n = 11), were also identified.

Genetic and transcriptional evolution alters cancer cell line drug response.

Human cancer cell lines are the workhorse of cancer research. Although cell lines are known to evolve in culture, the extent of the resultant genetic and transcriptional heterogeneity and its functional consequences remain understudied. Here we use genomic analyses of 106 human cell lines grown in two laboratories to show extensive clonal diversity. Further comprehensive genomic characterization of 27 strains of the common breast cancer cell line MCF7 uncovered rapid genetic diversification. Similar results were obtained with multiple strains of 13 additional cell lines. Notably, genetic changes were associated with differential activation of gene expression programs and marked differences in cell morphology and proliferation. Barcoding experiments showed that cell line evolution occurs as a result of positive clonal selection that is highly sensitive to culture conditions. Analyses of single-cell-derived clones demonstrated that continuous instability quickly translates into heterogeneity of the cell line. When the 27 MCF7 strains were tested against 321 anti-cancer compounds, we uncovered considerably different drug responses: at least 75% of compounds that strongly inhibited some strains were completely inactive in others. This study documents the extent, origins and consequences of genetic variation within cell lines, and provides a framework for researchers to measure such variation in efforts to support maximally reproducible cancer research.

Loss of heterozygosity does not occur in BRCA1/2 mutant pediatric solid and central nervous system tumors.

Utilization of tumor-only sequencing has expanded in pediatric cancer patients, which can lead to identification of pathogenic variants in genes that may be germline and/or have uncertain relevance to the tumor in question, such as the homologous recombination (HR) pathway genes BRCA1/2. We identified patients with pathogenic BRCA1/2 mutations from somatic tumor sequencing, and performed additional germline sequencing to assess for the presence of loss of heterozygosity (LOH). Of seven patients identified, four (57.1%) mutations were found in the germline and none had associated LOH. Our data suggest that BRCA1/2 mutations identified in this context are likely incidental findings.

Sporadic facial angiofibroma and sporadic angiomyolipoma mimicking tuberous sclerosis complex.

Tuberous sclerosis complex (TSC) is a genetic syndrome due to mutations in either TSC1 or TSC2, leading to the development of hamartomatous tumours at multiple body sites, including facial skin (facial angiofibroma (FAF)), brain (cortical tubers) and kidney (angiomyolipoma (AML)). In this report, we describe an individual with minimal TSC clinical features, who had 'no mutation identified' (NMI) by prior genetic testing in a clinical laboratory. Our massively parallel sequencing (MPS) analysis of multiple samples from different body sites and tumours (including blood, saliva, normal skin, AML and FAF) revealed an extraordinary situation in which FAF and AML had completely independent inactivating biallelic variants in TSC2, not present in other matched samples. This suggests that the two different lesions (AML and FAF) are not due to the same underlying germline or mosaic mutation, rather both are likely sporadic events. This case demonstrates the relevance of thorough clinical examination, high-coverage MPS of multiple tumours and matched normal tissues, and appropriate genetic counselling for individuals with marginal TSC features and possible TSC1 or TSC2 mosaicism.

Implementation of a prostate cancer-specific targeted sequencing panel for credentialing of patient-derived cell lines and genomic characterization of patient samples.

BACKGROUND: Primary and metastatic prostate cancers have low mutation rates and recurrent alterations in a small set of genes, enabling targeted sequencing of prostate cancer-associated genes as an efficient approach to characterizing patient samples (compared to whole-exome and whole-genome sequencing). For example, targeted sequencing provides a flexible, rapid, and cost-effective method for genomic assessment of patient-derived cell lines to evaluate fidelity to initial patient tumor samples. METHODS: We developed a prostate cancer-specific targeted next-generation sequencing (NGS) panel to detect alterations in 62 prostate cancer-associated genes as well as recurring gene fusions with ETS family members, representing the majority of common alterations in prostate cancer. We tested this panel on primary prostate cancer tissues and blood biopsies from patients with metastatic prostate cancer. We generated patient-derived cell lines from primary prostate cancers using conditional reprogramming methods and applied targeted sequencing to evaluate the fidelity of these cell lines to the original patient tumors. RESULTS: The prostate cancer-specific panel identified biologically and clinically relevant alterations, including point mutations in driver oncogenes and ETS family fusion genes, in tumor tissues from 29 radical prostatectomy samples. The targeted panel also identified genomic alterations in cell-free DNA and circulating tumor cells (CTCs) from patients with metastatic prostate cancer, and in standard prostate cancer cell lines. We used the targeted panel to sequence our set of patient-derived cell lines; however, no prostate cancer-specific mutations were identified in the tumor-derived cell lines, suggesting preferential outgrowth of normal prostate epithelial cells. CONCLUSIONS: We evaluated a prostate cancer-specific targeted NGS panel to detect common and clinically relevant alterations (including ETS family gene fusions) in prostate cancer. The panel detected driver mutations in a diverse set of clinical samples of prostate cancer, including fresh-frozen tumors, cell-free DNA, CTCs, and cell lines. Targeted sequencing of patient-derived cell lines highlights the challenge of deriving cell lines from primary prostate cancers and the importance of genomic characterization to credential candidate cell lines. Our study supports that a prostate cancer-specific targeted sequencing panel provides an efficient, clinically feasible approach to identify genetic alterations across a spectrum of prostate cancer samples and cell lines.

Adverse prognostic impact of the loss of STAG2 protein expression in patients with newly diagnosed localised Ewing sarcoma: A report from the Children's Oncology Group.

BACKGROUND: Ewing sarcoma (EWS) is an aggressive sarcoma with no validated molecular biomarkers. We aimed to determine the frequency of STAG2 protein loss by immunohistochemistry (IHC) and whether loss of expression is associated with outcome. METHODS: We performed a retrospective cohort study of patients with EWS enrolled to Children's Oncology Group studies. We obtained unstained slides from 235 patients and DNA for sequencing from 75 patients. STAG2 expression was tested for association with clinical features and survival was estimated using Kaplan-Meier methods with log-rank tests. RESULTS: In total, 155 cases passed quality control for STAG2 IHC. STAG2 expression in 20/155 cases could not be categorised with the limited available tissue, leaving 135 patients with definitive STAG2 IHC. In localised and metastatic disease, STAG2 was lost in 29/108 and 6/27 cases, respectively. Among patients with IHC and sequencing, 0/17 STAG2 expressing cases had STAG2 mutations, and 2/7 cases with STAG2 loss had STAG2 mutations. Among patients with localised disease, 5-year event-free survival was 54% (95% CI 34-70%) and 75% (95% CI 63-84%) for patients with STAG2 loss vs. expression (P = 0.0034). CONCLUSION: STAG2 loss of expression is identified in a population of patients without identifiable STAG2 mutations and carries a poor prognosis.

Genomic analyses of PMBL reveal new drivers and mechanisms of sensitivity to PD-1 blockade.

Primary mediastinal large B-cell lymphomas (PMBLs) are aggressive tumors that typically present as large mediastinal masses in young women. PMBLs share clinical, transcriptional, and molecular features with classical Hodgkin lymphoma (cHL), including constitutive activation of nuclear factor κB (NF-κB), JAK/STAT signaling, and programmed cell death protein 1 (PD-1)-mediated immune evasion. The demonstrated efficacy of PD-1 blockade in relapsed/refractory PMBLs led to recent approval by the US Food and Drug Administration and underscored the importance of characterizing targetable genetic vulnerabilities in this disease. Here, we report a comprehensive analysis of recurrent genetic alterations -somatic mutations, somatic copy number alterations, and structural variants-in a cohort of 37 newly diagnosed PMBLs. We identified a median of 9 genetic drivers per PMBL, including known and newly identified components of the JAK/STAT and NF-κB signaling pathways and frequent B2M alterations that limit major histocompatibility complex class I expression, as in cHL. PMBL also exhibited frequent, newly identified driver mutations in ZNF217 and an additional epigenetic modifier, EZH2. The majority of these alterations were clonal, which supports their role as early drivers. In PMBL, we identified several previously uncharacterized molecular features that may increase sensitivity to PD-1 blockade, including high tumor mutational burden, microsatellite instability, and an apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) mutational signature. The shared genetic features between PMBL and cHL provide a framework for analyzing the mechanism of action of PD-1 blockade in these related lymphoid malignancies.

Plasma cell-free DNA variant analysis compared with methylated DNA analysis in renal cell carcinoma.

PURPOSE: Plasma cell-free DNA (cfDNA) variant analysis is commonly used in many cancer subtypes. Cell-free methylated DNA immunoprecipitation sequencing (cfMeDIP-seq) has shown high sensitivity for cancer detection. To date, studies have not compared the sensitivity of both methods in a single cancer subtype. METHODS: cfDNA from 40 metastatic RCC (mRCC) patients was subjected to targeted panel variant analysis. For 34 of 40, cfMeDIP-seq was also performed. A separate cohort of 38 mRCC patients were used in cfMeDIP-seq analysis to train an RCC classifier. RESULTS: cfDNA variant analysis detected 21 candidate variants in 11 of 40 mRCC patients (28%), after exclusion of 2 germline variants and 6 variants reflecting clonal hematopoiesis. Among 23 patients with parallel tumor sequencing, cfDNA analysis alone identified variants in 9 patients (39%), while cfDNA analysis focused on tumor sequencing variant findings improved the sensitivity to 52%. In 34 mRCC patients undergoing cfMeDIP-seq, cfDNA variant analysis identified variants in 7 (21%), while cfMeDIP-seq detected all mRCC cases (100% sensitivity) with 88% specificity in 34 control subjects. In 5 patients with cfDNA variants and serial samples, variant frequency correlated with response to therapy. CONCLUSION: cfMeDIP-seq is significantly more sensitive for mRCC detection than cfDNA variant analysis. However, cfDNA variant analysis may be useful for monitoring response to therapy.

Correction: Detection of circulating tumour DNA is associated with inferior outcomes in Ewing sarcoma and osteosarcoma: a report from the Children's Oncology Group.

The authors have noticed that the final paragraph of the Results section contains errors in the number of patients involved. The correct number of patients is included in the text below. These errors do not affect the Figure referenced.In osteosarcoma, we focused on 8q gain as a specific biological feature of interest. Among the 41 patients with detectable ctDNA in the osteosarcoma cohort, 8q gain was detected in 73.2% (30/41). The 3-year EFS for patients with 8q gain (n = 30) in ctDNA was 60.0% (95% CI 40.5-75.0) compared to 80.8 (95% CI 42.4-94.9) in patients without 8q gain (n = 11) in ctDNA (p = 0.18; Fig. 3).

Detection of Mutations in Barrett's Esophagus Before Progression to High-Grade Dysplasia or Adenocarcinoma.

BACKGROUND & AIMS: Barrett's esophagus (BE) is the greatest risk factor for esophageal adenocarcinoma (EAC), but only a small proportion of patients with BE develop cancer. Biomarkers might be able to identify patients at highest risk of progression. We investigated genomic differences in surveillance biopsies collected from patients whose BE subsequently progressed compared to patients whose disease did not progress. METHODS: We performed a retrospective case-control study of 24 patients with BE that progressed to high-grade dysplasia (HGD, n = 14) or EAC (n = 10). The control group (n = 73, called non-progressors) comprised patients with BE and at least 5 years of total endoscopic biopsy surveillance without progression to HGD or EAC. From each patient, we selected a single tissue sample obtained more than 1 year before progression (cases) or more than 2 years before the end of follow-up (controls). Pathogenic mutations, gene copy numbers, and ploidy were compared between samples from progressors and non-progressors. RESULTS: TP53 mutations were detected in 46% of samples from progressors and 5% of non-progressors. In this case-control sample set, TP53 mutations in BE tissues increased the adjusted risk of progression 13.8-fold (95% confidence interval, 3.2-61.0) (P < .001). We did not observe significant differences in ploidy or copy-number profile between groups. We identified 147 pathogenic mutations in 57 distinct genes-the average number of pathogenic mutations was higher in samples from progressors (n = 2.5) than non-progressors (n = 1.2) (P < .001). TP53 and other somatic mutations were recurrently detected in samples with limited copy-number changes (aneuploidy). CONCLUSIONS: In genomic analyses of BE tissues from patients with or without later progression to HGD or EAC, we found significantly higher numbers of TP53 mutations in BE from patients with subsequent progression. These mutations were frequently detected before the onset of dysplasia or substantial changes in copy number.

Low-level mosaicism in tuberous sclerosis complex: prevalence, clinical features, and risk of disease transmission.

PURPOSE: To examine the prevalence and spectrum of mosaic variant allele frequency (MVAF) in tuberous sclerosis complex (TSC) patients with low-level mosaicism and correlate genetic findings with clinical features and transmission risk. METHODS: Massively parallel sequencing was performed on 39 mosaic TSC patients with 170 different tissue samples. RESULTS: TSC mosaic patients (MVAF: 0-10%, median 1.7% in blood DNA) had a milder and distinct clinical phenotype in comparison with other TSC series, with similar facial angiofibromas (92%) and kidney angiomyolipomas (83%), and fewer seizures, cortical tubers, and multiple other manifestations (p < 0.0001 for six features). MVAF of TSC1/TSC2 pathogenic variants was highly variable in different tissue samples. Remarkably, skin lesions were the most reliable tissue for variant identification, and 6 of 39 (15%) patients showed no evidence of the variant in blood. Semen analysis showed absence of the variant in 3 of 5 mosaic men. The expected distribution of MVAF in comparison with that observed here suggests that there is a considerable number of individuals with low-level mosaicism for a TSC2 pathogenic variant who are not recognized clinically. CONCLUSION: Our findings provide information on variability in MVAF and risk of transmission that has broad implications for other mosaic genetic disorders.

Germline cancer susceptibility gene variants, somatic second hits, and survival outcomes in patients with resected pancreatic cancer.

PURPOSE: Germline variants in double-strand DNA damage repair (dsDDR) genes (e.g., BRCA1/2) predispose to pancreatic adenocarcinoma (PDAC) and may predict sensitivity to platinum-based chemotherapy and poly(ADP) ribose polymerase (PARP) inhibitors. We sought to determine the prevalence and significance of germline cancer susceptibility gene variants in PDAC with paired somatic and survival analyses. METHODS: Using a customized next-generation sequencing panel, germline/somatic DNA was analyzed from 289 patients with resected PDAC ascertained without preselection for high-risk features (e.g., young age, personal/family history). All identified variants were assessed for pathogenicity. Outcomes were analyzed using multivariable-adjusted Cox proportional hazards regression. RESULTS: We found that 28/289 (9.7%; 95% confidence interval [CI] 6.5-13.7%) patients carried pathogenic/likely pathogenic germline variants, including 21 (7.3%) dsDDR gene variants (3 BRCA1, 4 BRCA2, 14 other dsDDR genes [ATM, BRIP1, CHEK2, NBN, PALB2, RAD50, RAD51C]), 3 Lynch syndrome, and 4 other genes (APC p.I1307K, CDKN2A, TP53). Somatic sequencing and immunohistochemistry identified second hits in the tumor in 12/27 (44.4%) patients with germline variants (1 failed sequencing). Compared with noncarriers, patients with germline dsDDR gene variants had superior overall survival (hazard ratio [HR] 0.54; 95% CI 0.30-0.99; P = 0.05). CONCLUSION: Nearly 10% of PDAC patients harbor germline variants, although the majority lack somatic second hits, the therapeutic significance of which warrants further study.

Unique, dual-indexed sequencing adapters with UMIs effectively eliminate index cross-talk and significantly improve sensitivity of massively parallel sequencing.

BACKGROUND: Sample index cross-talk can result in false positive calls when massively parallel sequencing (MPS) is used for sensitive applications such as low-frequency somatic variant discovery, ancient DNA investigations, microbial detection in human samples, or circulating cell-free tumor DNA (ctDNA) variant detection. Therefore, the limit-of-detection of an MPS assay is directly related to the degree of index cross-talk. RESULTS: Cross-talk rates up to 0.29% were observed when using standard, combinatorial adapters, resulting in 110,180 (0.1% cross-talk rate) or 1,121,074 (0.29% cross-talk rate) misassigned reads per lane in non-patterned and patterned Illumina flow cells, respectively. Here, we demonstrate that using unique, dual-matched indexed adapters dramatically reduces index cross-talk to ≤1 misassigned reads per flow cell lane. While the current study was performed using dual-matched indices, using unique, dual-unrelated indices would also be an effective alternative. CONCLUSIONS: For sensitive downstream analyses, the use of combinatorial indices for multiplexed hybrid capture and sequencing is inappropriate, as it results in an unacceptable number of misassigned reads. Cross-talk can be virtually eliminated using dual-matched indexed adapters. These results suggest that use of such adapters is critical to reduce false positive rates in assays that aim to identify low allele frequency events, and strongly indicate that dual-matched adapters be implemented for all sensitive MPS applications.

Clinical and mutational spectrum of highly differentiated, paired box 3:forkhead box protein o1 fusion-negative rhabdomyosarcoma: A report from the Children's Oncology Group.

BACKGROUND: Pediatric paired box 3:forkhead box protein O1 fusion-negative (PF-) rhabdomyosarcoma (RMS) represents a diverse spectrum of tumors with marked differences in histology, myogenic differentiation, and clinical behavior. METHODS: This study sought to evaluate the clinical and mutational spectrum of 24 pediatric PF- human RMS tumors with high levels of myogenic differentiation. Tumors were sequenced with OncoPanel v.2, a panel consisting of the coding regions of 504 genes previously linked to human cancer. RESULTS: Most of the tumors (19 of 24) arose at head/neck or genitourinary sites, and the overall survival rate was 100% with a median follow-up time of 4.6 years (range, 1.4-8.6 years). RAS pathway gene mutations were the most common mutations in PF-, highly differentiated RMS tumors. In addition, Hedgehog (Hh) and mechanistic target of rapamycin (mTOR) gene mutations with evidence for functional relevance (high-impact) were identified in subsets of tumors. The presence of Hh and mTOR pathway gene mutations was mutually exclusive and was associated with high-impact RAS pathway gene mutations in 3 of 4 Hh-mutated tumors and in 1 of 6 mTOR-mutated tumors. CONCLUSIONS: Interestingly, Hh and mTOR gene mutations were previously associated with rhabdomyomas, which are also known to preferentially arise at head/neck and genitourinary sites. Findings from this study further support the idea that PF-, highly differentiated RMS tumors and rhabdomyomas may represent a continuous spectrum of tumors. Cancer 2018;124:1973-81. © 2018 American Cancer Society.

Detection of circulating tumour DNA is associated with inferior outcomes in Ewing sarcoma and osteosarcoma: a report from the Children's Oncology Group.

BACKGROUND: New prognostic markers are needed to identify patients with Ewing sarcoma (EWS) and osteosarcoma unlikely to benefit from standard therapy. We describe the incidence and association with outcome of circulating tumour DNA (ctDNA) using next-generation sequencing (NGS) assays. METHODS: A NGS hybrid capture assay and an ultra-low-pass whole-genome sequencing assay were used to detect ctDNA in banked plasma from patients with EWS and osteosarcoma, respectively. Patients were coded as positive or negative for ctDNA and tested for association with clinical features and outcome. RESULTS: The analytic cohort included 94 patients with EWS (82% from initial diagnosis) and 72 patients with primary localised osteosarcoma (100% from initial diagnosis). ctDNA was detectable in 53% and 57% of newly diagnosed patients with EWS and osteosarcoma, respectively. Among patients with newly diagnosed localised EWS, detectable ctDNA was associated with inferior 3-year event-free survival (48.6% vs. 82.1%; p = 0.006) and overall survival (79.8% vs. 92.6%; p = 0.01). In both EWS and osteosarcoma, risk of event and death increased with ctDNA levels. CONCLUSIONS: NGS assays agnostic of primary tumour sequencing results detect ctDNA in half of the plasma samples from patients with newly diagnosed EWS and osteosarcoma. Detectable ctDNA is associated with inferior outcomes.

Targetable genetic features of primary testicular and primary central nervous system lymphomas.

Primary central nervous system lymphomas (PCNSLs) and primary testicular lymphomas (PTLs) are extranodal large B-cell lymphomas (LBCLs) with inferior responses to current empiric treatment regimens. To identify targetable genetic features of PCNSL and PTL, we characterized their recurrent somatic mutations, chromosomal rearrangements, copy number alterations (CNAs), and associated driver genes, and compared these comprehensive genetic signatures to those of diffuse LBCL and primary mediastinal large B-cell lymphoma (PMBL). These studies identify unique combinations of genetic alterations in discrete LBCL subtypes and subtype-selective bases for targeted therapy. PCNSLs and PTLs frequently exhibit genomic instability, and near-uniform, often biallelic, CDKN2A loss with rare TP53 mutations. PCNSLs and PTLs also use multiple genetic mechanisms to target key genes and pathways and exhibit near-uniform oncogenic Toll-like receptor signaling as a result of MYD88 mutation and/or NFKBIZ amplification, frequent concurrent B-cell receptor pathway activation, and deregulation of BCL6. Of great interest, PCNSLs and PTLs also have frequent 9p24.1/PD-L1/PD-L2 CNAs and additional translocations of these loci, structural bases of immune evasion that are shared with PMBL.

Diffuse large B-cell lymphoma patient-derived xenograft models capture the molecular and biological heterogeneity of the disease.

Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous disease defined by transcriptional classifications, specific signaling and survival pathways, and multiple low-frequency genetic alterations. Preclinical model systems that capture the genetic and functional heterogeneity of DLBCL are urgently needed. Here, we generated and characterized a panel of large B-cell lymphoma (LBCL) patient-derived xenograft (PDX) models, including 8 that reflect the immunophenotypic, transcriptional, genetic, and functional heterogeneity of primary DLBCL and 1 that is a plasmablastic lymphoma. All LBCL PDX models were subjected to whole-transcriptome sequencing to classify cell of origin and consensus clustering classification (CCC) subtypes. Mutations and chromosomal rearrangements were evaluated by whole-exome sequencing with an extended bait set. Six of the 8 DLBCL models were activated B-cell (ABC)-type tumors that exhibited ABC-associated mutations such as MYD88, CD79B, CARD11, and PIM1. The remaining 2 DLBCL models were germinal B-cell type, with characteristic alterations of GNA13, CREBBP, and EZH2, and chromosomal translocations involving IgH and either BCL2 or MYC Only 25% of the DLBCL PDX models harbored inactivating TP53 mutations, whereas 75% exhibited copy number alterations of TP53 or its upstream modifier, CDKN2A, consistent with the reported incidence and type of p53 pathway alterations in primary DLBCL. By CCC criteria, 6 of 8 DLBCL PDX models were B-cell receptor (BCR)-type tumors that exhibited selective surface immunoglobulin expression and sensitivity to entospletinib, a recently developed spleen tyrosine kinase inhibitor. In summary, we have established and characterized faithful PDX models of DLBCL and demonstrated their usefulness in functional analyses of proximal BCR pathway inhibition.

Mosaic and Intronic Mutations in TSC1/TSC2 Explain the Majority of TSC Patients with No Mutation Identified by Conventional Testing.

Tuberous sclerosis complex (TSC) is an autosomal dominant tumor suppressor gene syndrome due to germline mutations in either TSC1 or TSC2. 10-15% of TSC individuals have no mutation identified (NMI) after thorough conventional molecular diagnostic assessment. 53 TSC subjects who were NMI were studied using next generation sequencing to search for mutations in these genes. Blood/saliva DNA including parental samples were available from all subjects, and skin tumor biopsy DNA was available from six subjects. We identified mutations in 45 of 53 subjects (85%). Mosaicism was observed in the majority (26 of 45, 58%), and intronic mutations were also unusually common, seen in 18 of 45 subjects (40%). Seventeen (38%) mutations were seen at an allele frequency < 5%, five at an allele frequency < 1%, and two were identified in skin tumor biopsies only, and were not seen at appreciable frequency in blood or saliva DNA. These findings illuminate the extent of mosaicism in TSC, indicate the importance of full gene coverage and next generation sequencing for mutation detection, show that analysis of TSC-related tumors can increase the mutation detection rate, indicate that it is not likely that a third TSC gene exists, and enable provision of genetic counseling to the substantial population of TSC individuals who are currently NMI.

BreaKmer: detection of structural variation in targeted massively parallel sequencing data using kmers.

Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for 'targeted' resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a 'kmer' strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.