Evaluation of ultra-deep targeted sequencing for personalized breast cancer care

Prevalence and spectrum of AKT1, PIK3CA, PTEN and TP53 somatic mutations in Chinese breast cancer patients

Breast cancer, one of the most frequently occurring cancers worldwide, is the leading cause of cancer-related death among women. AKT1, PIK3CA, PTEN and TP53 mutations were common observed in breast cancer representing potential clinical biomarkers for cancer classification and treatment. A comprehensive knowledge of AKT1, PIK3CA, PTEN and TP53 mutations in breast cancer was still insufficient in Chinese population. In this study, the complete coding regions and exon-intron boundaries of AKT1, PIK3CA, PTEN and TP53 genes were sequenced in paired breast tumor and normal tissues from 313 Chinese breast cancer patients using microfluidic PCR-based target enrichment and next-generation sequencing technology. Total 120 somatic mutations were identified in 190 of the 313 patients (60.7%), with the mutation frequency of AKT1 as 3.2%, PIK3CA as 36.4%, PTEN as 4.8%, and TP53 as 33.9%. Among these mutations, 1 in PIK3CA (p.I69N), 3 in PTEN (p.K62X, c.635-12_636delTTAACCATGCAGAT and p.N340IfsTer4) and 5 in TP53 (p.Q136AfsTer5, p.K139_P142del, p.Y234dup, p.V274LfsTer31 and p.N310TfsTer35) were novel. Notably, PIK3CA somatic mutations were significantly associated with ER-positive or PR-positive tumors. TP53 somatic mutations were significantly associated with ER-negative, PR-negative, HER2-positive, BRCA1 mutation, Ki67 high expression and basal-like tumors. Our findings provided a comprehensive mutation profiling of AKT1, PIK3CA, PTEN and TP53 genes in Chinese breast cancer patients, which have potential implications in clinical management.

Genetic variants of prospectively demonstrated phenocopies in BRCA1/2 kindreds

Background

In kindreds carrying path_BRCA1/2 variants, some women in these families will develop cancer despite testing negative for the family’s pathogenic variant. These families may have additional genetic variants, which not only may increase the susceptibility of the families’ path_BRCA1/2, but also be capable of causing cancer in the absence of the path_BRCA1/2 variants. We aimed to identify novel genetic variants in prospectively detected breast cancer (BC) or gynecological cancer cases tested negative for their families’ pathogenic BRCA1/2 variant (path_BRCA1 or path_BRCA2).

Methods

Women with BC or gynecological cancer who had tested negative for path_BRCA1 or path_BRCA2 variants were included. Forty-four cancer susceptibility genes were screened for genetic variation through a targeted amplicon-based sequencing assay. Protein- and RNA splicing-dedicated in silico analyses were performed for all variants of unknown significance (VUS). Variants predicted as the ones most likely affecting pre-mRNA splicing were experimentally analyzed in a minigene assay.

Results

We identified 48 women who were tested negative for their family’s path_BRCA1 (n = 13) or path_BRCA2 (n = 35) variants. Pathogenic variants in the ATM, BRCA2, MSH6 and MUTYH genes were found in 10% (5/48) of the cases, of whom 15% (2/13) were from path_BRCA1 and 9% (3/35) from path_BRCA2 families. Out of the 26 unique VUS, 3 (12%) were predicted to affect RNA splicing (APC c.721G > A, MAP3K1 c.764A > G and MSH2 c.815C > T). However, by using a minigene, assay we here show that APC c.721G > A does not cause a splicing defect, similarly to what has been recently reported for the MAP3K1 c.764A > G. The MSH2 c.815C > T was previously described as causing partial exon skipping and it was identified in this work together with the path_BRCA2 c.9382C > T (p.R3128X).

Conclusion

All women in breast or breast/ovarian cancer kindreds would benefit from being offered genetic testing irrespective of which causative genetic variants have been demonstrated in their relatives.

Electronic supplementary material

The online version of this article (10.1186/s13053-018-0086-0) contains supplementary material, which is available to authorized users.

Exome and deep sequencing of clinically aggressive neuroblastoma reveal somatic mutations that affect key pathways involved in cancer progression

The spectrum of somatic mutation of the most aggressive forms of neuroblastoma is not completely determined. We sought to identify potential cancer drivers in clinically aggressive neuroblastoma.

Whole exome sequencing was conducted on 17 germline and tumor DNA samples from high-risk patients with adverse events within 36 months from diagnosis (HR-Event3) to identify somatic mutations and deep targeted sequencing of 134 genes selected from the initial screening in additional 48 germline and tumor pairs (62.5% HR-Event3 and high-risk patients), 17 HR-Event3 tumors and 17 human-derived neuroblastoma cell lines.

We revealed 22 significantly mutated genes, many of which implicated in cancer progression. Fifteen genes (68.2%) were highly expressed in neuroblastoma supporting their involvement in the disease. CHD9, a cancer driver gene, was the most significantly altered (4.0% of cases) after ALK.

Other genes (PTK2, NAV3, NAV1, FZD1 and ATRX), expressed in neuroblastoma and involved in cell invasion and migration were mutated at frequency ranged from 4% to 2%.

Focal adhesion and regulation of actin cytoskeleton pathways, were frequently disrupted (14.1% of cases) thus suggesting potential novel therapeutic strategies to prevent disease progression.

Notably BARD1, CHEK2 and AXIN2 were enriched in rare, potentially pathogenic, germline variants.

In summary, whole exome and deep targeted sequencing identified novel cancer genes of clinically aggressive neuroblastoma. Our analyses show pathway-level implications of infrequently mutated genes in leading neuroblastoma progression.

Evolving Therapeutic Strategies to Exploit Chromosome Instability in Cancer

Cancer is a devastating disease that claims over 8 million lives each year. Understanding the molecular etiology of the disease is critical to identify and develop new therapeutic strategies and targets. Chromosome instability (CIN) is an abnormal phenotype, characterized by progressive numerical and/or structural chromosomal changes, which is observed in virtually all cancer types. CIN generates intratumoral heterogeneity, drives cancer development, and promotes metastatic progression, and thus, it is associated with highly aggressive, drug-resistant tumors and poor patient prognosis. As CIN is observed in both primary and metastatic lesions, innovative strategies that exploit CIN may offer therapeutic benefits and better outcomes for cancer patients. Unfortunately, exploiting CIN remains a significant challenge, as the aberrant mechanisms driving CIN and their causative roles in cancer have yet to be fully elucidated. The development and utilization of CIN-exploiting therapies is further complicated by the associated risks for off-target effects and secondary cancers. Accordingly, this review will assess the strengths and limitations of current CIN-exploiting therapies, and discuss emerging strategies designed to overcome these challenges to improve outcomes and survival for patients diagnosed with cancer.

Identifying DNase I hypersensitive sites as driver distal regulatory elements in breast cancer

Efforts to identify driver mutations in cancer have largely focused on genes, whereas non-coding sequences remain relatively unexplored. Here we develop a statistical method based on characteristics known to influence local mutation rate and a series of enrichment filters in order to identify distal regulatory elements harboring putative driver mutations in breast cancer. We identify ten DNase I hypersensitive sites that are significantly mutated in breast cancers and associated with the aberrant expression of neighboring genes. A pan-cancer analysis shows that three of these elements are significantly mutated across multiple cancer types and have mutation densities similar to protein-coding driver genes. Functional characterization of the most highly mutated DNase I hypersensitive sites in breast cancer (using in silico and experimental approaches) confirms that they are regulatory elements and affect the expression of cancer genes. Our study suggests that mutations of regulatory elements in tumors likely play an important role in cancer development.

Cancer driver mutations can occur within noncoding genomic sequences. Here, the authors develop a statistical approach to identify candidate noncoding driver mutations in DNase I hypersensitive sites in breast cancer and experimentally demonstrate they are regulatory elements of known cancer genes.

Evaluation of whole exome sequencing by targeted gene sequencing and Sanger sequencing

BACKGROUND

Targeted gene sequencing (TGS) and whole exome sequencing (WES) are being used in clinical testing in laboratories. We compared the performances of TGS and WES using the same DNA samples.

METHODS

DNA was extracted from 10 endometrial tumor tissue specimens. Sequencing were performed with an Illumina HiSeq 2000. We randomly selected variants to confirm through Sanger sequencing or mutant-enriched PCR with Sanger sequencing.

RESULTS

We found that the variants identified in both TGS and WES were true positives (47/47), regardless of the sequencing depth. Most variants found in TGS only were true positives (34/40), and most of the variants found by WES only were false positives (8/18). From these results, we suggest that the sequencing depth may not play important role in the accuracy of NGS-based methods. After analysis, we found that WES had a sensitivity of 72.70%, specificity of 96.27%, precision of 99.44%, and accuracy of 75.03%.

CONCLUSIONS

The results of NGS-based methods must currently be validated, especially for important reported variants regardless of the methods used, and for the use of WES in cancers a higher false negative rate must be considered. More sensitive methods should be used to confirm the NGS results in uneven cancer tissues.

Somatic copy number alterations detected by ultra-deep targeted sequencing predict prognosis in oral cavity squamous cell carcinoma

Background

Ultra-deep targeted sequencing (UDT-Seq) has advanced our knowledge on the incidence and functional significance of somatic mutations. However, the utility of UDT-Seq in detecting copy number alterations (CNAs) remains unclear. With the goal of improving molecular prognostication and identifying new therapeutic targets, we designed this study to assess whether UDT-Seq may be useful for detecting CNA in oral cavity squamous cell carcinoma (OSCC).

Methods

We sequenced a panel of clinically actionable cancer mutations in 310 formalin-fixed paraffin-embedded OSCC specimens. A linear model was developed to overcome uneven coverage across target regions and multiple samples. The 5-year rates of secondary primary tumors, local recurrence, neck recurrence, distant metastases, and survival served as the outcome measures. We confirmed the prognostic significance of the CNA signatures in an independent sample of 105 primary OSCC specimens.

Results

The CNA burden across 10 targeted genes was found to predict prognosis in two independent cohorts. FGFR1 and PIK3CAamplifications were associated with prognosis independent of clinical risk factors. Genes exhibiting CNA were clustered in the proteoglycan metabolism, the FOXO signaling, and the PI3K-AKT signaling pathways, for which targeted drugs are already available or currently under development.

Conclusions

UDT-Seq is clinically useful to identify CNA, which significantly improve the prognostic information provided by traditional clinicopathological risk factors in OSCC patients.

VarDict: a novel and versatile variant caller for next-generation sequencing in cancer research

Accurate variant calling in next generation sequencing (NGS) is critical to understand cancer genomes better. Here we present VarDict, a novel and versatile variant caller for both DNA- and RNA-sequencing data. VarDict simultaneously calls SNV, MNV, InDels, complex and structural variants, expanding the detected genetic driver landscape of tumors. It performs local realignments on the fly for more accurate allele frequency estimation. VarDict performance scales linearly to sequencing depth, enabling ultra-deep sequencing used to explore tumor evolution or detect tumor DNA circulating in blood. In addition, VarDict performs amplicon aware variant calling for polymerase chain reaction (PCR)-based targeted sequencing often used in diagnostic settings, and is able to detect PCR artifacts. Finally, VarDict also detects differences in somatic and loss of heterozygosity variants between paired samples. VarDict reprocessing of The Cancer Genome Atlas (TCGA) Lung Adenocarcinoma dataset called known driver mutations in KRAS, EGFR, BRAF, PIK3CA and MET in 16% more patients than previously published variant calls. We believe VarDict will greatly facilitate application of NGS in clinical cancer research.

A Systematic Assessment of Accuracy in Detecting Somatic Mosaic Variants by Deep Amplicon Sequencing: Application to NF2 Gene

The accurate detection of low-allelic variants is still challenging, particularly for the identification of somatic mosaicism, where matched control sample is not available. High throughput sequencing, by the simultaneous and independent analysis of thousands of different DNA fragments, might overcome many of the limits of traditional methods, greatly increasing the sensitivity. However, it is necessary to take into account the high number of false positives that may arise due to the lack of matched control samples. Here, we applied deep amplicon sequencing to the analysis of samples with known genotype and variant allele fraction (VAF) followed by a tailored statistical analysis. This method allowed to define a minimum value of VAF for detecting mosaic variants with high accuracy. Then, we exploited the estimated VAF to select candidate alterations in NF2 gene in 34 samples with unknown genotype (30 blood and 4 tumor DNAs), demonstrating the suitability of our method. The strategy we propose optimizes the use of deep amplicon sequencing for the identification of low abundance variants. Moreover, our method can be applied to different high throughput sequencing approaches to estimate the background noise and define the accuracy of the experimental design.

Biased estimates of clonal evolution and subclonal heterogeneity can arise from PCR duplicates in deep sequencing experiments

Accurate allele frequencies are important for measuring subclonal heterogeneity and clonal evolution. Deep-targeted sequencing data can contain PCR duplicates, inflating perceived read depth. Here we adapted the Illumina TruSeq Custom Amplicon kit to include single molecule tagging (SMT) and show that SMT-identified duplicates arise from PCR. We demonstrate that retention of PCR duplicate reads can imply clonal evolution when none exists, while their removal effectively controls the false positive rate. Additionally, PCR duplicates alter estimates of subclonal heterogeneity in tumor samples. Our method simplifies PCR duplicate identification and emphasizes their removal in studies of tumor heterogeneity and clonal evolution.

Mung bean nuclease treatment increases capture specificity of microdroplet-PCR based targeted DNA enrichment

Targeted DNA enrichment coupled with next generation sequencing has been increasingly used for interrogation of select sub-genomic regions at high depth of coverage in a cost effective manner. Specificity measured by on-target efficiency is a key performance metric for target enrichment. Non-specific capture leads to off-target reads, resulting in waste of sequencing throughput on irrelevant regions. Microdroplet-PCR allows simultaneous amplification of up to thousands of regions in the genome and is among the most commonly used strategies for target enrichment. Here we show that carryover of single-stranded template genomic DNA from microdroplet-PCR constitutes a major contributing factor for off-target reads in the resultant libraries. Moreover, treatment of microdroplet-PCR enrichment products with a nuclease specific to single-stranded DNA alleviates off-target load and improves enrichment specificity. We propose that nuclease treatment of enrichment products should be incorporated in the workflow of targeted sequencing using microdroplet-PCR for target capture. These findings may have a broad impact on other PCR based applications for which removal of template DNA is beneficial.