PatternCNV: a versatile tool for detecting copy number changes from exome sequencing data

Breast Cancer Risk Modification in Women with Pathogenic Variants in BRCA1, BRCA2, ATM, CHEK2, and PALB2

SIGNIFICANCE

There is limited information on whether established risk factors increase breast cancer risk from PVs. In the WHI, PV carriers had no substantial (≥2-fold) increase with most risk factors, except potentially MHT in ATM or CHEK2 carriers. The results may inform counseling and research on MHT.

Susceptibility gene mutations in germline and tumors of patients with HER2-negative advanced breast cancer

Germline mutations in BRCA1 and BRCA2 (gBRCA1/2) are required for a PARP inhibitor therapy in patients with HER2-negative (HER2-) advanced breast cancer (aBC). However, little is known about the prognostic impact of gBRCA1/2 mutations in aBC patients treated with chemotherapy. This study aimed to investigate the frequencies and prognosis of germline and somatic BRCA1/2 mutations in HER2- aBC patients receiving the first chemotherapy in the advanced setting. Patients receiving their first chemotherapy for HER2- aBC were retrospectively selected from the prospective PRAEGNANT registry (NCT02338167). Genotyping of 26 cancer predisposition genes was performed with germline DNA of 471 patients and somatic tumor DNA of 94 patients. Mutation frequencies, progression-free and overall survival (PFS, OS) according to germline mutation status were assessed. gBRCA1/2 mutations were present in 23 patients (4.9%), and 33 patients (7.0%) had mutations in other cancer risk genes. Patients with a gBRCA1/2 mutation had a better OS compared to non-mutation carriers (HR: 0.38; 95%CI: 0.17-0.86). PFS comparison was not statistically significant. Mutations in other risk genes did not affect prognosis. Two somatic BRCA2 mutations were found in 94 patients without gBRCA1/2 mutations. Most frequently somatic mutated genes were TP53 (44.7%), CDH1 (10.6%) and PTEN (6.4%). In conclusion, aBC patients with gBRCA1/2 mutations had a more favorable prognosis under chemotherapy compared to non-mutation carriers. The mutation frequency of ~5% with gBRCA1/2 mutations together with improved outcome indicates that germline genotyping of all metastatic patients for whom a PARP inhibitor therapy is indicated should be considered.

Utility of Targeted Sequencing Compared to FISH for Detection of Chronic Lymphocytic Leukemia Copy Number Alterations

Chronic lymphocytic leukemia (CLL) is characterized by multiple copy number alterations (CNAs) and somatic mutations that are central to disease prognosis, risk stratification, and mechanisms of therapy resistance. Fluorescence in situ hybridization (FISH) panels are widely used in clinical applications as the gold standard for screening prognostic chromosomal abnormalities in CLL. DNA sequencing is an alternative approach to identifying CNAs but is not an established method for clinical CNA screening. We sequenced DNA from 509 individuals with CLL or monoclonal B-cell lymphocytosis (MBL), the precursor to CLL, using a targeted sequencing panel of 59 recurrently mutated genes in CLL and additional amplicons across regions affected by clinically relevant CNAs [i.e., del(17p), del(11q), del(13q), and trisomy 12]. We used the PatternCNV algorithm to call CNA and compared the concordance of calling clinically relevant CNAs by targeted sequencing to that of FISH. We found a high accuracy of calling CNAs via sequencing compared to FISH. With FISH as the gold standard, the specificity of targeted sequencing was >95%, sensitivity was >86%, positive predictive value was >90%, and negative predictive value was >84% across the clinically relevant CNAs. Using targeted sequencing, we were also able to identify other common CLL-associated CNAs, including del(6q), del(14q), and gain 8q, as well as complex karyotype, defined as the presence of 3 or more chromosomal abnormalities, in 26 patients. In a single and cost-effective assay that can be performed on stored DNA samples, targeted sequencing can simultaneously detect CNAs, somatic mutations, and complex karyotypes, which are all important prognostic features in CLL.

Follicular lymphoma B cells exhibit heterogeneous transcriptional states with associated somatic alterations and tumor microenvironments

Follicular lymphoma (FL) is an indolent non-Hodgkin lymphoma of germinal center origin, which presents with significant biologic and clinical heterogeneity. Using RNA-seq on B cells sorted from 87 FL biopsies, combined with machine-learning approaches, we identify 3 transcriptional states that divide the biological ontology of FL B cells into inflamed, proliferative, and chromatin-modifying states, with relationship to prior GC B cell phenotypes. When integrated with whole-exome sequencing and immune profiling, we find that each state was associated with a combination of mutations in chromatin modifiers, copy-number alterations to TNFAIP3, and T follicular helper cells (Tfh) cell interactions, or primarily by a microenvironment rich in activated T cells. Altogether, these data define FL B cell transcriptional states across a large cohort of patients, contribute to our understanding of FL heterogeneity at the tumor cell level, and provide a foundation for guiding therapeutic intervention.

Differential prognosis of single and multiple TP53 abnormalities in high-count MBL and untreated CLL

TP53 aberrations, including mutations and deletion of 17p13, are important adverse prognostic markers in chronic lymphocytic leukemia (CLL) but are less studied in high count monoclonal B-cell lymphocytosis (HCMBL), an asymptomatic pre-malignant stage of CLL. Here we estimated the prevalence and impact of TP53 aberrations in 1,230 newly diagnosed treatment-naïve individuals (849 CLL, 381 HCMBL). We defined TP53 state as: wild-type (no TP53 mutations and normal 17p), single-hit (del(17p) or one TP53 mutation), or multi-hit (TP53 mutation and del(17p), TP53 mutation and loss of heterozygosity, or multiple TP53 mutations). Cox regression was used to estimate hazard ratios (HR) and 95% confidence intervals (CI) for time to first treatment and overall survival by TP53 state. We found 64 (7.5%) CLL patients and 17 (4.5%) HCMBL individuals had TP53 mutations with variant allele fraction >10%. Del(17p) was present in 58 (6.8%) of CLL and 11 (2.9%) of HCMBL cases. Most individuals had wild-type (N=1,128, 91.7%) TP53 state, followed by multi-hit (N=55, 4.5%) and then single-hit (N=47, 3.8%) TP53 state. The risk of shorter time to therapy and death increased with the number of TP53 abnormalities. Compared to wild-type patients, multi-hit patients had 3-fold and single-hit patients had 1.5-fold increased risk of requiring therapy. Multi-hit patients also had 2.9-fold increased risk of death compared to wild-type. These results remained stable after accounting for other known poor prognostic factors. Both TP53 mutations and del(17p) may provide important prognostic information for HCMBL and CLL that would be missed if only one were measured.

Integration of multiomics data shows down regulation of mismatch repair and tubulin pathways in triple-negative chemotherapy-resistant breast tumors

BACKGROUND

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype. Patients with TNBC are primarily treated with neoadjuvant chemotherapy (NAC). The response to NAC is prognostic, with reductions in overall survival and disease-free survival rates in those patients who do not achieve a pathological complete response (pCR). Based on this premise, we hypothesized that paired analysis of primary and residual TNBC tumors following NAC could identify unique biomarkers associated with post-NAC recurrence.

METHODS AND RESULTS

We investigated 24 samples from 12 non-LAR TNBC patients with paired pre- and post-NAC data, including four patients with recurrence shortly after surgery (< 24 months) and eight who remained recurrence-free (> 48 months). These tumors were collected from a prospective NAC breast cancer study (BEAUTY) conducted at the Mayo Clinic. Differential expression analysis of pre-NAC biopsies showed minimal gene expression differences between early recurrent and nonrecurrent TNBC tumors; however, post-NAC samples demonstrated significant alterations in expression patterns in response to intervention. Topological-level differences associated with early recurrence were implicated in 251 gene sets, and an independent assessment of microarray gene expression data from the 9 paired non-LAR samples available in the NAC I-SPY1 trial confirmed 56 gene sets. Within these 56 gene sets, 113 genes were observed to be differentially expressed in the I-SPY1 and BEAUTY post-NAC studies. An independent (n = 392) breast cancer dataset with relapse-free survival (RFS) data was used to refine our gene list to a 17-gene signature. A threefold cross-validation analysis of the gene signature with the combined BEAUTY and I-SPY1 data yielded an average AUC of 0.88 for six machine-learning models. Due to the limited number of studies with pre- and post-NAC TNBC tumor data, further validation of the signature is needed.

CONCLUSION

Analysis of multiomics data from post-NAC TNBC chemoresistant tumors showed down regulation of mismatch repair and tubulin pathways. Additionally, we identified a 17-gene signature in TNBC associated with post-NAC recurrence enriched with down-regulated immune genes.

Comparison of biochemical, microbial and mucosal mRNA expression in bile acid diarrhoea and irritable bowel syndrome with diarrhoea

OBJECTIVE

There are altered mucosal functions in irritable bowel syndrome with diarrhoea (IBS-D); ~30% of patients with IBS-D have abnormal bile acid (BA) metabolism (ABAM) and diarrhoea (summarised as BAD).

AIM

To compare biochemical parameters, gastrointestinal and colonic transit, rectal sensation and pathobiological mechanisms in IBS-D without ABAM and in BAD (serum 7C4>52 ng/mL).

DESIGN

In patients with Rome III criteria of IBS-D, we compared biochemical features, colonic transit, rectal sensation, deep genotype of five BA-related genes, ileal and colonic mucosal mRNA (differential expression (DE) analysis) and stool dysbiosis (including functional analysis of microbiome). Results in BAD were compared with IBS-D without ABAM.

RESULTS

Compared with 161 patients with IBS-D without ABAM, 44 patients with BAD had significantly faster colonic transit, lower microbial alpha diversity, different compositional profile (beta diversity) and higher Firmicutes to Bacteroidetes ratio with evidence of decreased expression of bile acid thiol ligase (involved in transformation of primary to secondary BAs) and decreased sulfatases. In BAD (compared with IBS-D without ABAM), terminal ileal biopsies showed downregulation of SLC44A5 (a BA transporter), and ascending colon biopsies showed upregulation in barrier-weakening genes (CLDN2), serine protease inhibitors, immune activation, cellular differentiation and a cellular transporter (FABP6; BA binding). No DE of genes was documented in descending colon biopsies. The two groups had similar rectal sensation.

CONCLUSION

Though sharing clinical symptoms with IBS-D, BAD is associated with biological differences and mechanisms that have potential to enhance diagnosis and treatment targeting barrier dysfunction, inflammatory and microbial changes.

Rosette-forming glioneuronal tumours are midline, FGFR1-mutated tumours

Rosette-forming glioneuronal tumour (RGNT) is a rare CNS WHO grade 1 brain neoplasm. According to WHO 2021, essential diagnostic criteria are a "biphasic histomorphology with neurocytic and a glial component, and uniform neurocytes forming rosettes and/or perivascular pseudorosettes associated with synaptophysin expression" and/or DNA methylation profile of RGNT whereas "FGFR1 mutation with co-occurring PIK3CA and/or NF1 mutation" are desirable criteria.

MATERIAL AND METHODS

We report a series of 46 cases fulfilling the essential pathological diagnostic criteria for RGNT. FGFR1 and PIK3CA hotspot mutations were searched for by multiplex digital PCR in all cases whereas DNA methylation profiling and/or PIK3R1 and NF1 alterations were analysed in a subset of cases.

RESULTS

Three groups were observed. The first one included 21 intracranial midline tumours demonstrating FGFR1 mutation associated with PIK3CA or PIK3R1 (n=19) or NF1 (n=1) or PIK3CA and NF1 (n=1) mutation. By DNA methylation profiling, 8 cases were classified as RGNT (they demonstrated FGFR1 and PIK3CA or PIK3R1 mutations). Group 2 comprised 11 cases associated with one single FGFR1 mutation. Group 3 included 6 cases classified as LGG other than RGNT (1/6 showed FGFR1 mutation and one a FGFR1 and NF1 mutation) and 8 cases without FGFR1 mutation. Groups 2 and 3 were enriched in lateral and spinal cases.

CONCLUSIONS

We suggest adding FGFR1 mutation and intracranial midline location as essential diagnostic criteria. When DNA methylation profiling is not available, a RGNT diagnosis remains certain in cases demonstrating characteristic pathological features and FGFR1 mutation associated with either PIK3CA or PIK3R1 mutation.

CoverageMaster: comprehensive CNV detection and visualization from NGS short reads for genetic medicine applications

CoverageMaster (CoM) is a copy number variation (CNV) calling algorithm based on depth-of-coverage maps designed to detect CNVs of any size in exome [whole exome sequencing (WES)] and genome [whole genome sequencing (WGS)] data. The core of the algorithm is the compression of sequencing coverage data in a multiscale Wavelet space and the analysis through an iterative Hidden Markov Model. CoM processes WES and WGS data at nucleotide scale resolution and accurately detects and visualizes full size range CNVs, including single or partial exon deletions and duplications. The results obtained with this approach support the possibility for coverage-based CNV callers to replace probe-based methods such as array comparative genomic hybridization and multiplex ligation-dependent probe amplification in the near future.

Development and Validation of a Next-Generation Sequencing Panel for Syndromic and Nonsyndromic Hearing Loss

BACKGROUND

Deafness and hearing loss are common conditions that can be seen independently or as part of a syndrome and are often mediated by genetic causes. We sought to develop and validate a hereditary hearing loss panel (HHLP) to detect single nucleotide variants (SNVs), insertions and deletions (indels), and copy number variants (CNVs) in 166 genes related to nonsyndromic and syndromic hearing loss.

METHODS

We developed a custom-capture next-generation sequencing (NGS) reagent to detect all coding regions, ±10 flanking bp, for the 166 genes related to nonsyndromic and syndromic hearing loss. Our validation consisted of testing 52 samples to establish accuracy, reproducibility, and analytical sensitivity. In addition to NGS, supplementary methods, including multiplex ligation-dependent probe amplification, long-range PCR, and Sanger sequencing, were used to ensure coverage of regions that had high complexity or homology.

RESULTS

We observed 100% positive and negative percentage agreement for detection of SNVs (n = 362), small indels (1-22 bp, n = 25), and CNVs (gains, n = 8; losses, n = 17). Finally, we showed that this assay was able to detect variants with a variant allele frequency ≥20% for SNVs and indels and ≥30% to 35% for CNVs.

CONCLUSIONS

We validated an HHLP that detects SNVs, indels, and CNVs in 166 genes related to syndromic and nonsyndromic hearing loss. The results of this assay can be utilized to confirm a diagnosis of hearing loss and related syndromic disorders associated with known causal genes.

Benchmarking germline CNV calling tools from exome sequencing data

Whole-exome sequencing is an attractive alternative to microarray analysis because of the low cost and potential ability to detect copy number variations (CNV) of various sizes (from 1-2 exons to several Mb). Previous comparison of the most popular CNV calling tools showed a high portion of false-positive calls. Moreover, due to a lack of a gold standard CNV set, the results are limited and incomparable. Here, we aimed to perform a comprehensive analysis of tools capable of germline CNV calling available at the moment using a single CNV standard and reference sample set. Compiling variants from previous studies with Bayesian estimation approach, we constructed an internal standard for NA12878 sample (pilot National Institute of Standards and Technology Reference Material) including 110,050 CNV or non-CNV exons. The standard was used to evaluate the performance of 16 germline CNV calling tools on the NA12878 sample and 10 correlated exomes as a reference set with respect to length distribution, concordance, and efficiency. Each algorithm had a certain range of detected lengths and showed low concordance with other tools. Most tools are focused on detection of a limited number of CNVs one to seven exons long with a false-positive rate below 50%. EXCAVATOR2, exomeCopy, and FishingCNV focused on detection of a wide range of variations but showed low precision. Upon unified comparison, the tools were not equivalent. The analysis performed allows choosing algorithms or ensembles of algorithms most suitable for a specific goal, e.g. population studies or medical genetics.

HBOS-CNV: A New Approach to Detect Copy Number Variations From Next-Generation Sequencing Data

Copy number variation (CNV) is a genomic mutation that plays an important role in tumor evolution and tumor genesis. Accurate detection of CNVs from next-generation sequencing (NGS) data is still a challenging task due to artifacts such as uneven mapped reads and unbalanced amplitudes of gains and losses. This study proposes a new approach called HBOS-CNV to detect CNVs from NGS data. The central point of HBOS-CNV is that it uses a new statistic, the histogram-based outlier score (HBOS), to evaluate the fluctuation of genome bins to determine those of changed copy numbers. In comparison with existing statistics in the evaluation of CNVs, HBOS is a non-linearly transformed value from the observed read depth (RD) value of each genome bin, having the potential ability to relieve the effects resulted from the above artifacts. In the calculation of HBOS values, a dynamic width histogram is utilized to depict the density of bins on the genome being analyzed, which can reduce the effects of noises partially contributed by mapping and sequencing errors. The evaluation of genome bins using such a new statistic can lead to less extremely significant CNVs having a high probability of detection. We evaluated this method using a large number of simulation datasets and compared it with four existing methods (CNVnator, CNV-IFTV, CNV-LOF, and iCopyDav). The results demonstrated that our proposed method outperforms the others in terms of sensitivity, precision, and F1-measure. Furthermore, we applied the proposed method to a set of real sequencing samples from the 1000 Genomes Project and determined a number of CNVs with biological meanings. Thus, the proposed method can be regarded as a routine approach in the field of genome mutation analysis for cancer samples.

Effect of Germline Mutations in Homologous Recombination Repair Genes on Overall Survival of Patients with Pancreatic Adenocarcinoma

PURPOSE

To compare the clinical characteristics and overall survival (OS) of germline mutation carriers in homologous recombination repair (HRR) genes and noncarriers with pancreatic ductal adenocarcinoma (PDAC).

EXPERIMENTAL DESIGN

Germline DNA from 3,078 patients with PDAC enrolled in a prospective registry at Mayo Clinic between 2000 and 2017 was analyzed for mutations in 37 cancer predisposition genes. Characteristics and OS of patients with mutations in eight genes (ATM, BARD1, BRCA1, BRCA2, BRIP1, PALB2, RAD51C, and RAD51D) involved in HRR were compared with patients testing negative for mutations in all 37 genes.

RESULTS

The 175 HRR mutation carriers and 2,730 noncarriers in the study had a median duration of follow-up of 9.9 years. HRR mutation carriers were younger (median age at diagnosis: 63 vs. 66 years, P < 0.001) and more likely to have metastatic disease at diagnosis (46% vs. 36%, P = 0.004). In a multivariable model adjusting for sex, age at diagnosis, and tumor staging, patients with germline HRR mutations had a significantly longer OS compared with noncarriers [HR, 0.83; 95% confidence interval (CI), 0.70-0.97; P = 0.02]. Further gene-level analysis demonstrated that germline ATM mutation carriers had longer OS compared with patients without germline mutations in any of the 37 genes (HR, 0.72; 95% CI, 0.55-0.94; P = 0.01).

CONCLUSIONS

This study demonstrates that germline mutation carrier status in PDAC is associated with longer OS compared with noncarriers. Further research into tumor biology and response to platinum-based chemotherapy in germline mutation carriers with PDAC are needed to better understand the association with longer OS.

Whole Exome Sequencing Among 26 Patients With Indeterminate Acute Liver Failure: A Pilot Study

INTRODUCTION

The etiology of acute liver failure (ALF) remains an important prognostic factor. The Acute Liver Failure Study Group recently reported that 150 of 2,718 adult patients with ALF (5.5%) had an indeterminate etiology. Our aim was to use whole exome sequencing to identify genetic variants associated with phenotypic, biochemical, and histologic features among patients with indeterminate ALF.

METHODS

This effort has defined a cohort of well-pedigreed patients with indeterminate ALF; DNA samples extracted from whole blood samples were obtained from 26 respective patients with indeterminate ALF. These samples were kept at the Acute Liver Failure Study Group repository at the NIDDK, Bethesda. Whole exome sequencing and bioinformatics analysis were performed at the Mayo Clinic Center of Individualized Medicine in Rochester, MN.

RESULTS

Of the 26 patients, 8 survived spontaneously, 6 died, and 12 underwent a liver transplantation; all those transplanted were alive at 21 days after enrollment in the study. Twenty-two of the 26 patients presented as ALF. We found 12 variants associated with 11 genes. The most common variant was rs4940595 in the SERPINB11 gene which was found in 23 of the 26 patients. This variant had a stop codon; no reports of disorders have been associated with this variant. The next most commonly found variant was rs1135840 in the CYP2D6 gene; this mutation is a missense_variant and has been reported to be associated with hepatotoxicity of antituberculous therapy. None of our patients were receiving this therapy. We also found a significant asymmetric distribution of rs1800754 of the CYP2D7 gene and rs1135840 of the CYP2D6 gene between patients who survived spontaneously (75%) and those who died or underwent liver transplantation (30.5% and 25%, respectively).

DISCUSSION

We found 12 variants of 11 genes significantly associated with ALF among adults with indeterminate etiology. We also found a significant asymmetric distribution of 2 variants belonging to the CYP2D7 and CYP2D6 genes, respectively, between those who survived spontaneously and those who died or underwent liver transplantation. The 2 most common variants, rs4940595 and rs1135840, of the SERPINB11 and CYP2D6 genes, respectively, found in our patients with ALF have been described as potentially important in the adaptive response combating the emergence of infectious diseases and associated with hepatotoxicity of antituberculous therapy, respectively. Our findings need to be expanded to include more patients with indeterminate ALF as well as viral, drug toxicity, and autoimmune etiologies to determine whether our findings are associated with the specific etiology, indeterminate, or with the overall ALF syndrome itself.

Genomic and Phenotypic Characterization of a Broad Panel of Patient-Derived Xenografts Reflects the Diversity of Glioblastoma

PURPOSE

Glioblastoma is the most frequent and lethal primary brain tumor. Development of novel therapies relies on the availability of relevant preclinical models. We have established a panel of 96 glioblastoma patient-derived xenografts (PDX) and undertaken its genomic and phenotypic characterization.

EXPERIMENTAL DESIGN

PDXs were established from glioblastoma, IDH-wildtype (n = 93), glioblastoma, IDH-mutant (n = 2), diffuse midline glioma, H3 K27M-mutant (n = 1), and both primary (n = 60) and recurrent (n = 34) tumors. Tumor growth rates, histopathology, and treatment response were characterized. Integrated molecular profiling was performed by whole-exome sequencing (WES, n = 83), RNA-sequencing (n = 68), and genome-wide methylation profiling (n = 76). WES data from 24 patient tumors was compared with derivative models.

RESULTS

PDXs recapitulate many key phenotypic and molecular features of patient tumors. Orthotopic PDXs show characteristic tumor morphology and invasion patterns, but largely lack microvascular proliferation and necrosis. PDXs capture common and rare molecular drivers, including alterations of TERT, EGFR, PTEN, TP53, BRAF, and IDH1, most at frequencies comparable with human glioblastoma. However, PDGFRA amplification was absent. RNA-sequencing and genome-wide methylation profiling demonstrated broad representation of glioblastoma molecular subtypes. MGMT promoter methylation correlated with increased survival in response to temozolomide. WES of 24 matched patient tumors showed preservation of most genetic driver alterations, including EGFR amplification. However, in four patient-PDX pairs, driver alterations were gained or lost on engraftment, consistent with clonal selection.

CONCLUSIONS

Our PDX panel captures the molecular heterogeneity of glioblastoma and recapitulates many salient genetic and phenotypic features. All models and genomic data are openly available to investigators.

Contribution of Inherited DNA-Repair Gene Mutations to Hormone-Sensitive and Castrate-Resistant Metastatic Prostate Cancer and Implications for Clinical Outcome

PURPOSE

To compare the prevalence of germline mutations in metastatic hormone-sensitive prostate cancer (mHSPC) and metastatic castrate-resistant prostate cancer (mCRPC) and assess the impact of mutations on progression to castration resistance and overall survival.

METHODS

Targeted sequencing of germline DNA from 704 men (221 at the time of mHSPC and 483 at the time of mCRPC) enrolled in two advanced prostate cancer registries at Mayo Clinic between 2003 and 2013 was performed for 21 predisposition genes. Frequencies of pathogenic mutations were compared in patients and reference controls to identify genes enriched in metastatic prostate cancer. Multivariable Cox proportional hazards regression was used to identify predictors of progression to mCRPC and overall survival.

RESULTS

Sixty-eight germline mutations in 12 genes were identified in 66 men (9.4%). Mutations in ATM, BRCA2, CHEK2, FANCM, and TP53 were significantly enriched (odds ratio greater than 2.0) in the metastatic cohorts compared with reference controls. The frequency of germline mutations was similar for patients with mHSPC and mCRPC (11.8% v 8.3%; P = .16). The median time to progression from mHSPC to mCRPC was 23.1 and 32.5 months for patients with and without mutations, respectively (P = .96). Although older age at diagnosis, Gleason score greater than 7, elevated alkaline phosphatase level, and high volume of disease were associated with shorter duration of progression to mCRPC and poor overall survival, mutation status was not (progression to mCRPC hazard ratio, 0.81; 95% CI, 0.61 to 1.09; P = .17; overall survival hazard ratio, 1.00; 95% CI, 0.75 to 1.34; P = .98).

CONCLUSION

Similarly elevated rates of germline predisposition gene mutations in mHSPC and mCRPC suggest that germline genetic testing may help to guide medical management for all patients with advanced metastatic prostate cancer. Mutation status was not associated with shorter progression to mCRPC or poor overall survival.

Copy Number Variation Detection Using Total Variation

Next-generation sequencing (NGS) technologies offer new opportunities for precise and accurate identification of genomic aberrations, including copy number variations (CNVs). For high-throughput NGS data, using depth of coverage has become a major approach to identify CNVs, especially for whole exome sequencing (WES) data. Due to the high level of noise and biases of read-count data and complexity of the WES data, existing CNV detection tools identify many false CNV segments. Besides, NGS generates a huge amount of data, requiring to use effective and efficient methods. In this work, we propose a novel segmentation algorithm based on the total variation approach to detect CNVs more precisely and efficiently using WES data. The proposed method also filters out outlier read-counts and identifies significant change points to reduce false positives. We used real and simulated data to evaluate the performance of the proposed method and compare its performance with those of other commonly used CNV detection methods. Using simulated and real data, we show that the proposed method outperforms the existing CNV detection methods in terms of accuracy and false discovery rate and has a faster runtime compared to the circular binary segmentation method.

Molecular Profiling Reclassifies Adult Astroblastoma into Known and Clinically Distinct Tumor Entities with Frequent Mitogen-Activated Protein Kinase Pathway Alterations

BACKGROUND

Astroblastoma (ABM) is a rare glial brain tumor. Recurrent meningioma 1 (MN1) alterations have been recently identified in most pediatric cases. Adolescent and adult cases, however, remain molecularly poorly defined.

MATERIALS AND METHODS

We performed clinical and molecular characterization of a retrospective cohort of 14 adult and 1 adolescent ABM.

RESULTS

Strikingly, we found that MN1 fusions are a rare event in this age group (1/15). Using methylation profiling and targeted sequencing, most cases were reclassified as either pleomorphic xanthoastrocytomas (PXA)-like or high-grade glioma (HGG)-like. PXA-like ABM show BRAF mutation (6/7 with V600E mutation and 1/7 with G466E mutation) and CD34 expression. Conversely, HGG-like ABM harbored specific alterations of diffuse midline glioma (2/5) or glioblastoma (GBM; 3/5). These latter patients showed an unfavorable clinical course with significantly shorter overall survival (p = .021). Mitogen-activated protein kinase pathway alterations (including FGFR fusion, BRAF and NF1 mutations) were present in 10 of 15 patients and overrepresented in the HGG-like group (3/5) compared with previously reported prevalence of these alterations in GBM and diffuse midline glioma.

CONCLUSION

We suggest that gliomas with astroblastic features include a variety of molecularly sharply defined entities. Adult ABM harboring molecular features of PXA and HGG should be reclassified. Central nervous system high-grade neuroepithelial tumors with MN1 alterations and histology of ABM appear to be uncommon in adults. Astroblastic morphology in adults should thus prompt thorough molecular investigation aiming at a clear histomolecular diagnosis and identifying actionable drug targets, especially in the mitogen-activated protein kinase pathway.

IMPLICATIONS FOR PRACTICE

Astroblastoma (ABM) remains a poorly defined and controversial entity. Although meningioma 1 alterations seem to define a large subset of pediatric cases, adult cases remain molecularly poorly defined. This comprehensive molecular characterization of 1 adolescent and 14 adult ABM revealed that adult ABM histology comprises several molecularly defined entities, which explains clinical diversity and identifies actionable targets. Namely, pleomorphic xanthoastrocytoma-like ABM cases show a favorable prognosis whereas high-grade glioma (glioblastoma and diffuse midline gliome)-like ABM show significantly worse clinical courses. These results call for in-depth molecular analysis of adult gliomas with astroblastic features for diagnostic and therapeutic purposes.

A Prospective Correlation of Tissue Histopathology With Nucleic Acid Yield in Metastatic Castration-Resistant Prostate Cancer Biopsy Specimens

Objective

To determine histopathologic, exome, and transcriptome nucleic acid material yield from prospectively collected metastatic tissue biopsy specimens in patients with metastatic castration-resistant prostate cancer (mCRPC).

Patients and Methods

Patients with mCRPC initiating abiraterone acetate therapy underwent 2 serial metastatic site core needle biopsies after study activation on May 17, 2013. Multiple cores were obtained, and from each core, 1- to 2-mm segments were separated and formalin fixed for histopathologic examination. Tumor purity was determined for DNA and RNA from the rest of the biopsy specimen. RNA quality was assessed by calculation of an RNA integrity number and a DV200 score.

Results

A total of 89 patients underwent 172 uniformly processed core needle biopsies (89 on visit 1 and 83 on visit 2) between May 30, 2013, and September 10, 2015. Metastatic sites biopsied included bone (131), lymph nodes (31), liver (5), lung (3), and pelvic soft tissues (2). Of the 172 biopsy specimens, 85 (49%) had at least one of the multiple cores positive for tumor on histopathologic examination (53 of 88 [60%] from visit 1 and 32 of 83 [39%] from visit 2; P=.006). Metastatic carcinoma was observed in 50 of 130 bone lesion specimens (38%), compared to 35 of 41 nonbone specimens (85%) (P<.001). More than 10% tumoral DNA purity was observed in 89% and 80% of visit 1 and visit 2 biopsy specimens, respectively. Similarly, more than 10% tumor RNA purity was observed in 79% of visit 1 vs 59% for visit 2 (P=.008). In all, 134 of 172 procedures (78%) yielded tumor material either by histopathologic or nucleic acid purity analysis.

Conclusion

This study found that biopsy specimens from mCRPC sites yield adequate histopathologic, exome, and transcriptome material in most, but not all, cases. This finding has relevance for future genome sequencing studies on the introduction of targeted therapeutic agents.

Trial Registration

clinicaltrials.gov Identifier: 01953640.

RNA-Seq detects a SAMD12-EXT1 fusion transcript and leads to the discovery of an EXT1 deletion in a child with multiple osteochondromas

Background

We describe a patient presenting with pachygyria, epilepsy, developmental delay, short stature, failure to thrive, facial dysmorphisms, and multiple osteochondromas.

Methods

The patient underwent extensive genetic testing and analysis in an attempt to diagnose the cause of his condition. Clinical testing included metaphase karyotyping, array comparative genomic hybridization, direct sequencing and multiplex ligation‐dependent probe amplification and trio‐based exome sequencing. Subsequently, research‐based whole transcriptome sequencing was conducted to determine whether it might shed light on the undiagnosed phenotype.

Results

Clinical exome sequencing of patient and parent samples revealed a maternally inherited splice‐site variant in the doublecortin (DCX) gene that was classified as likely pathogenic and diagnostic of the patient's neurological phenotype. Clinical array comparative genome hybridization analysis revealed a 16p13.3 deletion that could not be linked to the patient phenotype based on affected genes. Further clinical testing to determine the cause of the patient's multiple osteochondromas was unrevealing despite extensive profiling of the most likely causative genes, EXT1 and EXT2, including mutation screening by direct sequence analysis and multiplex ligation‐dependent probe amplification. Whole transcriptome sequencing identified a SAMD12‐EXT1 fusion transcript that could have resulted from a chromosomal deletion, leading to the loss of EXT1 function. Re‐review of the clinical array comparative genomic hybridization results indicated a possible unreported mosaic deletion affecting the SAMD12 and EXT1 genes that corresponded precisely to the introns predicted to be affected by a fusion‐causing deletion. The existence of the mosaic deletion was subsequently confirmed clinically by an increased density copy number array and orthogonal methodologies

Conclusions

While mosaic mutations and deletions of EXT1 and EXT2 have been reported in the context of multiple osteochondromas, to our knowledge, this is the first time that transcriptomics technologies have been used to diagnose a patient via fusion transcript analysis in the congenital disease setting.

Genetic burden and associations with adverse neurodevelopment in neonates with congenital heart disease

BACKGROUND

Up to 20% of children with congenital heart disease (CHD) undergoing cardiac surgery develop neurodevelopmental disabilities (NDD), with some studies reporting persistent impairment. Recent large-scale studies have demonstrated shared genetic mechanisms contributing to CHD and NDD. In this study, a targeted approach was applied to assess direct clinical applicability of this information.

METHODS

A gene panel comprising 148 known CHD and/or NDD genes was used to sequence 15 patients with CHD + NDD, 15 patients with CHD, and 15 healthy controls. The number and types of variants between the 3 groups were compared using Poisson log-linear regression, and the SNP-set (Sequence) Kernel Association Test-Optimized was used to conduct single-gene and gene-pathway burden analyses.

RESULTS

A significant increase in rare (minor allele frequency < 0.01) and novel variants was identified between the CHD + NDD cohort and controls, P < .001 and P = .001, respectively. There was also a significant increase in rare variants in the CHD cohort compared with controls (P = .04). Rare variant burden analyses implicated pathways associated with "neurotransmitters," "axon guidance," and those incorporating "RASopathy" genes in the development of NDD in CHD patients.

CONCLUSIONS

These findings suggest that an increase in novel and rare variants in known CHD and/or NDD genes is associated with the development of NDD in patients with CHD. Furthermore, burden analyses point toward rare variant burden specifically in pathways related to brain development and function as contributors to NDD. Although promising variants and pathways were identified, further research, utilizing whole-genome approaches, is required prior to demonstrating clinical utility in this patient group.

Genomic approaches to the assessment of human spina bifida risk

Structural birth defects are a leading cause of mortality and morbidity in children world-wide, affecting as much as 6% of all live births. Among these conditions, neural tube defects (NTDs), including spina bifida and anencephaly, arise from a combination of complex gene and environment interactions that are as yet poorly understood within human populations. Rapid advances in massively parallel DNA sequencing and bioinformatics allow for analyses of the entire genome beyond the 2% of the genomic sequence covering protein coding regions. Efforts to collect and analyze these large datasets hold promise for illuminating gene network variations and eventually epigenetic events that increase individual risk for failure to close the neural tube. In this review, we discuss current challenges for DNA genome sequence analysis of NTD affected populations, and compare experience in the field with other complex genetic disorders for which large datasets are accumulating. The ultimate goal of this research is to find strategies for optimizing conditions that promote healthy birth outcomes for individual couples. Birth Defects Research 109:120-128, 2017. © 2016 Wiley Periodicals, Inc.

Myopathy With SQSTM1 and TIA1 Variants: Clinical and Pathological Features

Objective

The aim of this study is to identify the molecular defect of three unrelated individuals with late-onset predominant distal myopathy; to describe the spectrum of phenotype resulting from the contributing role of two variants in genes located on two different chromosomes; and to highlight the underappreciated complex forms of genetic myopathies.

Patients and methods

Clinical and laboratory data of three unrelated probands with predominantly distal weakness manifesting in the sixth-seventh decade of life, and available affected and unaffected family members were reviewed. Next-generation sequencing panel, whole exome sequencing, and targeted analyses of family members were performed to elucidate the genetic etiology of the myopathy.

Results

Genetic analyses detected two contributing variants located on different chromosomes in three unrelated probands: a heterozygous pathogenic mutation in SQSTM1 (c.1175C>T, p.Pro392Leu) and a heterozygous variant in TIA1 (c.1070A>G, p.Asn357Ser). The affected fraternal twin of one proband also carries both variants, while the unaffected family members harbor one or none. Two unrelated probands (family 1, II.3, and family 3, II.1) have a distal myopathy with rimmed vacuoles that manifested with index extensor weakness; the other proband (family 2, I.1) has myofibrillar myopathy manifesting with hypercapnic respiratory insufficiency and distal weakness.

Conclusion

The findings indicate that all the affected individuals have a myopathy associated with both variants in SQSTM1 and TIA1, respectively, suggesting that the two variants determine the phenotype and likely functionally interact. We speculate that the TIA1 variant is a modifier of the SQSTM1 mutation. We identify the combination of SQSTM1 and TIA1 variants as a novel genetic defect associated with myofibrillar myopathy and suggest to consider sequencing both genes in the molecular investigation of myopathy with rimmed vacuoles and myofibrillar myopathy although additional studies are needed to investigate the digenic nature of the disease.

Recessive TAF1A mutations reveal ribosomopathy in siblings with end-stage pediatric dilated cardiomyopathy

Non-ischemic dilated cardiomyopathy (DCM) has been recognized as a heritable disorder for over 25 years, yet clinical genetic testing is non-diagnostic in >50% of patients, underscoring the ongoing need for DCM gene discovery. Here, whole exome sequencing uncovered a novel molecular basis for idiopathic end-stage heart failure in two sisters who underwent cardiac transplantation at three years of age. Compound heterozygous recessive mutations in TAF1A, encoding an RNA polymerase I complex protein, were associated with marked fibrosis of explanted hearts and gene-specific nucleolar segregation defects in cardiomyocytes, indicative of impaired ribosomal RNA synthesis. Knockout of the homologous gene in zebrafish recapitulated a heart failure phenotype with pericardial edema, decreased ventricular systolic function, and embryonic mortality. These findings expand the clinical spectrum of ribosomopathies to include pediatric DCM.

A short review of variants calling for single-cell-sequencing data with applications

The field of single-cell sequencing is fleetly expanding, and many techniques have been developed in the past decade. With this technology, biologists can study not only the heterogeneity between two adjacent cells in the same tissue or organ, but also the evolutionary relationships and degenerative processes in a single cell. Calling variants is the main purpose in analyzing single cell sequencing (SCS) data. Currently, some popular methods used for bulk-cell-sequencing data analysis are tailored directly to be applied in dealing with SCS data. However, SCS requires an extra step of genome amplification to accumulate enough quantity for satisfying sequencing needs. The amplification yields large biases and thus raises challenge for using the bulk-cell-sequencing methods. In order to provide guidance for the development of specialized analyzed methods as well as using currently developed tools for SNS, this paper aims to bridge the gap. In this paper, we firstly introduced two popular genome amplification methods and compared their capabilities. Then we introduced a few popular models for calling single-nucleotide polymorphisms and copy-number variations. Finally, break-through applications of SNS were summarized to demonstrate its potential in researching cell evolution.

Targeted sequencing of 36 known or putative colorectal cancer susceptibility genes

BACKGROUND

Mutations in several genes predispose to colorectal cancer. Genetic testing for hereditary colorectal cancer syndromes was previously limited to single gene tests; thus, only a very limited number of genes were tested, and rarely those infrequently mutated in colorectal cancer. Next-generation sequencing technologies have made it possible to sequencing panels of genes known and suspected to influence colorectal cancer susceptibility.

METHODS

Targeted sequencing of 36 known or putative CRC susceptibility genes was conducted for 1231 CRC cases from five subsets: (1) Familial Colorectal Cancer Type X (n = 153); (2) CRC unselected by tumor immunohistochemical or microsatellite stability testing (n = 548); (3) young onset (age <50 years) (n = 333); (4) proficient mismatch repair (MMR) in cases diagnosed at ≥50 years (n = 68); and (5) deficient MMR CRCs with no germline mutations in MLH1, MSH2, MSH6, or PMS2 (n = 129). Ninety-three unaffected controls were also sequenced.

RESULTS

Overall, 29 nonsense, 43 frame-shift, 13 splice site, six initiator codon variants, one stop codon, 12 exonic deletions, 658 missense, and 17 indels were identified. Missense variants were reviewed by genetic counselors to determine pathogenicity; 13 were pathogenic, 61 were not pathogenic, and 584 were variants of uncertain significance. Overall, we identified 92 cases with pathogenic mutations in APC,MLH1,MSH2,MSH6, or multiple pathogenic MUTYH mutations (7.5%). Four cases with intact MMR protein expression by immunohistochemistry carried pathogenic MMR mutations.

CONCLUSIONS

Results across case subsets may help prioritize genes for inclusion in clinical gene panel tests and underscore the issue of variants of uncertain significance both in well-characterized genes and those for which limited experience has accumulated.

An evaluation of copy number variation detection tools for cancer using whole exome sequencing data

BACKGROUND

Recently copy number variation (CNV) has gained considerable interest as a type of genomic/genetic variation that plays an important role in disease susceptibility. Advances in sequencing technology have created an opportunity for detecting CNVs more accurately. Recently whole exome sequencing (WES) has become primary strategy for sequencing patient samples and study their genomics aberrations. However, compared to whole genome sequencing, WES introduces more biases and noise that make CNV detection very challenging. Additionally, tumors' complexity makes the detection of cancer specific CNVs even more difficult. Although many CNV detection tools have been developed since introducing NGS data, there are few tools for somatic CNV detection for WES data in cancer.

RESULTS

In this study, we evaluated the performance of the most recent and commonly used CNV detection tools for WES data in cancer to address their limitations and provide guidelines for developing new ones. We focused on the tools that have been designed or have the ability to detect cancer somatic aberrations. We compared the performance of the tools in terms of sensitivity and false discovery rate (FDR) using real data and simulated data. Comparative analysis of the results of the tools showed that there is a low consensus among the tools in calling CNVs. Using real data, tools show moderate sensitivity (~50% - ~80%), fair specificity (~70% - ~94%) and poor FDRs (~27% - ~60%). Also, using simulated data we observed that increasing the coverage more than 10× in exonic regions does not improve the detection power of the tools significantly.

CONCLUSIONS

The limited performance of the current CNV detection tools for WES data in cancer indicates the need for developing more efficient and precise CNV detection methods. Due to the complexity of tumors and high level of noise and biases in WES data, employing advanced novel segmentation, normalization and de-noising techniques that are designed specifically for cancer data is necessary. Also, CNV detection development suffers from the lack of a gold standard for performance evaluation. Finally, developing tools with user-friendly user interfaces and visualization features can enhance CNV studies for a broader range of users.

High prevalence of TP53 mutations is associated with poor survival and an EMT signature in gliosarcoma patients

Gliosarcoma (GS) is a rare variant (2%) of glioblastoma (GBM) that poses clinical genomic challenges because of its poor prognosis and limited genomic information. To gain a comprehensive view of the genomic alterations in GS and to understand the molecular etiology of GS, we applied whole-exome sequencing analyses for 28 GS cases (6 blood-matched fresh-frozen tissues for the discovery set, 22 formalin-fixed paraffin-embedded tissues for the validation set) and copy-number variation microarrays for 5 blood-matched fresh-frozen tissues. TP53 mutations were more prevalent in the GS cases (20/28, 70%) compared to the GBM cases (29/90, 32%), and the GS patients with TP53 mutations showed a significantly shorter survival (multivariate Cox analysis, hazard ratio=23.9, 95% confidence interval, 2.87-199.63, P=0.003). A pathway analysis showed recurrent alterations in MAPK signaling (EGFR, RASGRF2 and TP53), phosphatidylinositol/calcium signaling (CACNA1s, PLCs and ITPRs) and focal adhesion/tight junction (PTEN and PAK3) pathways. Genomic profiling of the matched recurrent GS cases detected the occurrence of TP53 mutations in two recurrent GS cases, which suggests that TP53 mutations play a role in treatment resistance. Functionally, we found that TP53 mutations are associated with the epithelial-mesenchymal transition (EMT) process of sarcomatous components of GS. We provide the first comprehensive genome-wide genetic alternation profiling of GS, which suggests novel prognostic subgroups in GS patients based on their TP53 mutation status and provides new insight in the pathogenesis and targeted treatment of GS.

Neurology Individualized Medicine: When to Use Next-Generation Sequencing Panels

Next-generation sequencing (NGS) is increasingly being applied to clinical testing. This practice is predicted to grow especially in neurology clinics because many of their patients have monogenetic causes for their "diagnostic odyssey." The cost of sequencing has been steadily decreasing, but the cost of DNA sequencing is a minor part of the total cost. Downstream data analysis, storage, and interpretation account for most of the total expense. In patients with nonspecific neurologic disorders in which an extensive number of genetic differential diagnoses exist, whole-genome sequencing (WGS) or whole-exome sequencing (WES) has shown promise in the identification of genetic causes. However, both WGS and WES have incomplete coverage and produce a large number of rare variants of unknown importance. In addition, ethical dilemmas are often created by unexpected findings in genes unrelated to the initial sequencing indication. Targeted-panel NGS starts with the capture of a set of disease-focused genes, followed by massive parallel sequencing. For many genetically heterogeneous neurologic disorders, a genetic panel that is disease focused yet inclusive of a large genetic differential diagnosis can be defined to reduce cost, increase turnaround time, and optimize performance. Targeted-panel NGS is currently the preferred first-tier approach because it provides a reliable clinical application while eliminating unexpected ethical dilemmas. Targeted-panel NGS is leading to a paradigm shift in the diagnosis of many neurologic disorders, enabling individualized precision medicine. In this review, we provide an overview of WGS, WES, and targeted-panel NGS in consideration of their utility in clinical testing for neurologic diseases.

Sequencing Structural Variants in Cancer for Precision Therapeutics

The identification of mutations that guide therapy selection for patients with cancer is now routine in many clinical centres. The majority of assays used for solid tumour profiling use DNA sequencing to interrogate somatic point mutations because they are relatively easy to identify and interpret. Many cancers, however, including high-grade serous ovarian, oesophageal, and small-cell lung cancer, are driven by somatic structural variants that are not measured by these assays. Therefore, there is currently an unmet need for clinical assays that can cheaply and rapidly profile structural variants in solid tumours. In this review we survey the landscape of 'actionable' structural variants in cancer and identify promising detection strategies based on massively-parallel sequencing.

Assessing the reproducibility of exome copy number variations predictions

Background

Reproducibility is receiving increased attention across many domains of science and genomics is no exception. Efforts to identify copy number variations (CNVs) from exome sequence (ES) data have been increasing. Many algorithms have been published to discover CNVs from exomes and a major challenge is the reproducibility in other datasets. Here we test exome CNV calling reproducibility under three conditions: data generated by different sequencing centers; varying sample sizes; and varying capture methodology.

Methods

Four CNV tools were tested: eXome Hidden Markov Model (XHMM), Copy Number Inference From Exome Reads (CoNIFER), EXCAVATOR, and Copy Number Analysis for Targeted Resequencing (CONTRA). To examine the reproducibility, we ran the callers on four datasets, varying sample sizes of N = 10, 30, 75, 100, 300, and data with different capture methodology. We examined the false negative (FN) calls and false positive (FP) calls for potential limitations of the CNV callers. The positive predictive value (PPV) was measured by checking the CNV call concordance against single nucleotide polymorphism array.

Results

Using independently generated datasets, we examined the PPV for each dataset and observed wide range of PPVs. The PPV values were highly data dependent (p <0.001). For the sample sizes and capture method analyses, we tested the callers in triplicates. Both analyses resulted in wide ranges of PPVs, even for the same test. Interestingly, negative correlations between the PPV and the sample sizes were observed for CoNIFER (ρ = –0.80). Further examination of FN calls showed that 44 % of these were missed by all callers and were attributed to the CNV size (46 % spanned ≤3 exons). Overlap of the FP calls showed that FPs were unique to each caller, indicative of algorithm dependency.

Conclusions

Our results demonstrate that further improvements in CNV callers are necessary to improve reproducibility and to include wider spectrum of CNVs (including the small CNVs). These CNV callers should be evaluated on multiple independent, heterogeneously generated datasets of varying size to increase robustness and utility. These approaches to the evaluation of exome CNV are essential to support wide utility and applicability of CNV discovery in exome studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0336-6) contains supplementary material, which is available to authorized users.

CNVkit: Genome-Wide Copy Number Detection and Visualization from Targeted DNA Sequencing

Germline copy number variants (CNVs) and somatic copy number alterations (SCNAs) are of significant importance in syndromic conditions and cancer. Massively parallel sequencing is increasingly used to infer copy number information from variations in the read depth in sequencing data. However, this approach has limitations in the case of targeted re-sequencing, which leaves gaps in coverage between the regions chosen for enrichment and introduces biases related to the efficiency of target capture and library preparation. We present a method for copy number detection, implemented in the software package CNVkit, that uses both the targeted reads and the nonspecifically captured off-target reads to infer copy number evenly across the genome. This combination achieves both exon-level resolution in targeted regions and sufficient resolution in the larger intronic and intergenic regions to identify copy number changes. In particular, we successfully inferred copy number at equivalent to 100-kilobase resolution genome-wide from a platform targeting as few as 293 genes. After normalizing read counts to a pooled reference, we evaluated and corrected for three sources of bias that explain most of the extraneous variability in the sequencing read depth: GC content, target footprint size and spacing, and repetitive sequences. We compared the performance of CNVkit to copy number changes identified by array comparative genomic hybridization. We packaged the components of CNVkit so that it is straightforward to use and provides visualizations, detailed reporting of significant features, and export options for integration into existing analysis pipelines. CNVkit is freely available from https://github.com/etal/cnvkit.

Target-enrichment sequencing and copy number evaluation in inherited polyneuropathy

OBJECTIVE

To assess the efficiency of target-enrichment next-generation sequencing (NGS) with copy number assessment in inherited neuropathy diagnosis.

METHODS

A 197 polyneuropathy gene panel was designed to assess for mutations in 93 patients with inherited or idiopathic neuropathy without known genetic cause. We applied our novel copy number variation algorithm on NGS data, and validated the identified copy number mutations using CytoScan (Affymetrix). Cost and efficacy of this targeted NGS approach was compared to earlier evaluations.

RESULTS

Average coverage depth was ∼760× (median = 600, 99.4% > 100×). Among 93 patients, 18 mutations were identified in 17 cases (18%), including 3 copy number mutations: 2 PMP22 duplications and 1 MPZ duplication. The 2 patients with PMP22 duplication presented with bulbar and respiratory involvement and had absent extremity nerve conductions, leading to axonal diagnosis. Average onset age of these 17 patients was 25 years (2-61 years), vs 45 years for those without genetic discovery. Among those with onset age less than 40 years, the diagnostic yield of targeted NGS approach is high (27%) and cost savings is significant (∼20%). However, the cost savings for patients with late onset age and without family history is not demonstrated.

CONCLUSIONS

Incorporating copy number analysis in target-enrichment NGS approach improved the efficiency of mutation discovery for chronic, inherited, progressive length-dependent polyneuropathy diagnosis. The new technology is facilitating a simplified genetic diagnostic algorithm utilizing targeted NGS, clinical phenotypes, age at onset, and family history to improve diagnosis efficiency. Our findings prompt a need for updating the current practice parameters and payer guidelines.

FOXA1 defines cancer cell specificity

A transcription factor functions differentially and/or identically in multiple cell types. However, the mechanism for cell-specific regulation of a transcription factor remains to be elucidated. We address how a single transcription factor, forkhead box protein A1 (FOXA1), forms cell-specific genomic signatures and differentially regulates gene expression in four human cancer cell lines (HepG2, LNCaP, MCF7, and T47D). FOXA1 is a pioneer transcription factor in organogenesis and cancer progression. Genomewide mapping of FOXA1 by chromatin immunoprecipitation sequencing annotates that target genes associated with FOXA1 binding are mostly common to these cancer cells. However, most of the functional FOXA1 target genes are specific to each cancer cell type. Further investigations using CRISPR-Cas9 genome editing technology indicate that cell-specific FOXA1 regulation is attributable to unique FOXA1 binding, genetic variations, and/or potential epigenetic regulation. Thus, FOXA1 controls the specificity of cancer cell types. We raise a "flower-blooming" hypothesis for cell-specific transcriptional regulation based on these observations.

Mutational Landscapes of Sequential Prostate Metastases and Matched Patient Derived Xenografts during Enzalutamide Therapy

Developing patient derived models from individual tumors that capture the biological heterogeneity and mutation landscape in advanced prostate cancer is challenging, but essential for understanding tumor progression and delivery of personalized therapy in metastatic castrate resistant prostate cancer stage. To demonstrate the feasibility of developing patient derived xenograft models in this stage, we present a case study wherein xenografts were derived from cancer metastases in a patient progressing on androgen deprivation therapy and prior to initiating pre-chemotherapy enzalutamide treatment. Tissue biopsies from a metastatic rib lesion were obtained for sequencing before and after initiating enzalutamide treatment over a twelve-week period and also implanted subcutaneously as well as under the renal capsule in immuno-deficient mice. The genome and transcriptome landscapes of xenografts and the original patient tumor tissues were compared by performing whole exome and transcriptome sequencing of the metastatic tumor tissues and the xenografts at both time points. After comparing the somatic mutations, copy number variations, gene fusions and gene expression we found that the patient's genomic and transcriptomic alterations were preserved in the patient derived xenografts with high fidelity. These xenograft models provide an opportunity for predicting efficacy of existing and potentially novel drugs that is based on individual metastatic tumor expression signature and molecular pharmacology for delivery of precision medicine.

SAAS-CNV: A Joint Segmentation Approach on Aggregated and Allele Specific Signals for the Identification of Somatic Copy Number Alterations with Next-Generation Sequencing Data

Cancer genomes exhibit profound somatic copy number alterations (SCNAs). Studying tumor SCNAs using massively parallel sequencing provides unprecedented resolution and meanwhile gives rise to new challenges in data analysis, complicated by tumor aneuploidy and heterogeneity as well as normal cell contamination. While the majority of read depth based methods utilize total sequencing depth alone for SCNA inference, the allele specific signals are undervalued. We proposed a joint segmentation and inference approach using both signals to meet some of the challenges. Our method consists of four major steps: 1) extracting read depth supporting reference and alternative alleles at each SNP/Indel locus and comparing the total read depth and alternative allele proportion between tumor and matched normal sample; 2) performing joint segmentation on the two signal dimensions; 3) correcting the copy number baseline from which the SCNA state is determined; 4) calling SCNA state for each segment based on both signal dimensions. The method is applicable to whole exome/genome sequencing (WES/WGS) as well as SNP array data in a tumor-control study. We applied the method to a dataset containing no SCNAs to test the specificity, created by pairing sequencing replicates of a single HapMap sample as normal/tumor pairs, as well as a large-scale WGS dataset consisting of 88 liver tumors along with adjacent normal tissues. Compared with representative methods, our method demonstrated improved accuracy, scalability to large cancer studies, capability in handling both sequencing and SNP array data, and the potential to improve the estimation of tumor ploidy and purity.

Concurrent Mutations in ATM and Genes Associated with Common γ Chain Signaling in Peripheral T Cell Lymphoma

Peripheral T cell lymphoma (PTCL) is a heterogeneous malignancy with poor response to current therapeutic strategies and incompletely characterized genetics. We conducted whole exome sequencing of matched PTCL and non-malignant samples from 12 patients, spanning 8 subtypes, to identify potential oncogenic mutations in PTCL. Analysis of the mutations identified using computational algorithms, CHASM, PolyPhen2, PROVEAN, and MutationAssessor to predict the impact of these mutations on protein function and PTCL tumorigenesis, revealed 104 somatic mutations that were selected as high impact by all four algorithms. Our analysis identified recurrent somatic missense or nonsense mutations in 70 genes, 9 of which contained mutations predicted significant by all 4 algorithms: ATM, RUNX1T1, WDR17, NTRK3, TP53, TRMT12, CACNA2D1, INTS8, and KCNH8. We observed somatic mutations in ATM (ataxia telangiectasia-mutated) in 5 out of the 12 samples and mutations in the common gamma chain (γc) signaling pathway (JAK3, IL2RG, STAT5B) in 3 samples, all of which also harbored mutations in ATM. Our findings contribute insights into the genetics of PTCL and suggest a relationship between γc signaling and ATM in T cell malignancy.

Prevalence of Pathogenic Mutations in Cancer Predisposition Genes among Pancreatic Cancer Patients

The prevalence of germline pathogenic mutations in a comprehensive panel of cancer predisposition genes is not well-defined for patients with pancreatic ductal adenocarcinoma (PDAC). To estimate the frequency of mutations in a panel of 22 cancer predisposition genes, 96 patients unselected for a family history of cancer who were recruited to the Mayo Clinic Pancreatic Cancer patient registry over a 12-month period were screened by next-generation sequencing. Fourteen pathogenic mutations in 13 patients (13.5%) were identified in eight genes: four in ATM, two in BRCA2, CHEK2, and MSH6, and one in BARD1, BRCA1, FANCM, and NBN. These included nine mutations (9.4%) in established pancreatic cancer genes. Three mutations were found in patients with a first-degree relative with PDAC, and 10 mutations were found in patients with first- or second-degree relatives with breast, pancreas, colorectal, ovarian, or endometrial cancers. These results suggest that a substantial proportion of patients with PDAC carry germline mutations in predisposition genes associated with other cancers and that a better understanding of pancreatic cancer risk will depend on evaluation of families with broad constellations of tumors. These findings highlight the need for recommendations governing germline gene-panel testing of patients with pancreatic cancer.

TP53 mutations, tetraploidy and homologous recombination repair defects in early stage high-grade serous ovarian cancer

To determine early somatic changes in high-grade serous ovarian cancer (HGSOC), we performed whole genome sequencing on a rare collection of 16 low stage HGSOCs. The majority showed extensive structural alterations (one had an ultramutated profile), exhibited high levels of p53 immunoreactivity, and harboured a TP53 mutation, deletion or inactivation. BRCA1 and BRCA2 mutations were observed in two tumors, with nine showing evidence of a homologous recombination (HR) defect. Combined Analysis with The Cancer Genome Atlas (TCGA) indicated that low and late stage HGSOCs have similar mutation and copy number profiles. We also found evidence that deleterious TP53 mutations are the earliest events, followed by deletions or loss of heterozygosity (LOH) of chromosomes carrying TP53, BRCA1 or BRCA2. Inactivation of HR appears to be an early event, as 62.5% of tumours showed a LOH pattern suggestive of HR defects. Three tumours with the highest ploidy had little genome-wide LOH, yet one of these had a homozygous somatic frame-shift BRCA2 mutation, suggesting that some carcinomas begin as tetraploid then descend into diploidy accompanied by genome-wide LOH. Lastly, we found evidence that structural variants (SV) cluster in HGSOC, but are absent in one ultramutated tumor, providing insights into the pathogenesis of low stage HGSOC.