panelcn.MOPS: Copy-number detection in targeted NGS panel data for clinical diagnostics

Targeted capture-based NGS is superior to multiplex PCR-based NGS for hereditary BRCA1 and BRCA2 gene analysis in FFPE tumor samples

BACKGROUND

With the introduction of Olaparib treatment for BRCA-deficient recurrent ovarian cancer, testing for somatic and/or germline mutations in BRCA1/2 genes in tumor tissues became essential for treatment decisions. In most cases only formalin-fixed paraffin-embedded (FFPE) samples, containing fragmented and chemically modified DNA of minor quality, are available. Thus, multiplex PCR-based sequencing is most commonly applied in routine molecular testing, which is predominantly focused on the identification of known hot spot mutations in oncogenes.

METHODS

We compared the overall performance of an adjusted targeted capture-based enrichment protocol and a multiplex PCR-based approach for calling of pathogenic SNVs and InDels using DNA extracted from 13 FFPE tissue samples. We further applied both strategies to seven blood samples and five matched FFPE tumor tissues of patients with known germline exon-spanning deletions and gene-wide duplications in BRCA1/2 to evaluate CNV detection based solely on panel NGS data. Finally, we analyzed DNA from FFPE tissues of 11 index patients from families suspected of having hereditary breast and ovarian cancer, of whom no blood samples were available for testing, in order to identify underlying pathogenic germline BRCA1/2 mutations.

RESULTS

The multiplex PCR-based protocol produced inhomogeneous coverage among targets of each sample and between samples as well as sporadic amplicon drop out, leading to insufficiently or non-covered nucleotides, which subsequently hindered variant detection. This protocol further led to detection of PCR-artifacts that could easily have been misinterpreted as pathogenic mutations. No such limitations were observed by application of an adjusted targeted capture-based protocol, which allowed for CNV calling with 86% sensitivity and 100% specificity. All pathogenic CNVs were confirmed in the five matched FFPE tumor samples from patients carrying known pathogenic germline mutations and we additionally identified somatic loss of the second allele in BRCA1/2. Furthermore we detected pathogenic BRCA1/2 variants in four the eleven FFPE samples from patients of whom no blood was available for analysis.

CONCLUSIONS

We demonstrate that an adjusted targeted capture-based enrichment protocol is superior to commonly applied multiplex PCR-based protocols for reliable BRCA1/2 variant detection, including CNV-detection, using FFPE tumor samples.

Detection of chromosome structural variation by targeted next-generation sequencing and a deep learning application

Molecular testing is increasingly important in cancer diagnosis. Targeted next generation sequencing (NGS) is widely accepted method but structural variation (SV) detection by targeted NGS remains challenging. In the brain tumor, identification of molecular alterations, including 1p/19q co-deletion, is essential for accurate glial tumor classification. Hence, we used targeted NGS to detect 1p/19q co-deletion using a newly developed deep learning (DL) model in 61 tumors, including 19 oligodendroglial tumors. An ensemble 1-dimentional convolution neural network was developed and used to detect the 1p/19q co-deletion. External validation was performed using 427 low-grade glial tumors from The Cancer Genome Atlas (TCGA). Manual review of the copy number plot from the targeted NGS identified the 1p/19q co-deletion in all 19 oligodendroglial tumors. Our DL model also perfectly detected the 1p/19q co-deletion (area under the curve, AUC = 1) in the testing set, and yielded reproducible results (AUC = 0.9652) in the validation set (n = 427), although the validation data were generated on a completely different platform (SNP Array 6.0 platform). In conclusion, targeted NGS using a cancer gene panel is a promising approach for classifying glial tumors, and DL can be successfully integrated for the SV detection in NGS data.

No Overt Clinical Immunodeficiency Despite Immune Biological Abnormalities in Patients With Constitutional Mismatch Repair Deficiency

Immunoglobulin class-switch recombination (CSR) and somatic hypermutations (SHMs) are prerequisites for antibody and immunoglobulin receptor maturation and adaptive immune diversity. The mismatch repair (MMR) machinery, consisting of homologs of MutSα, MutLα, and MutSβ (MSH2/MSH6, MLH1/PMS2, and MSH2/MSH3, respectively) and other proteins, is involved in CSR, primarily acting as a backup for nonhomologous end-joining repair of activation-induced cytidine deaminase-induced DNA mismatches and, furthermore, in addition to error-prone polymerases, in the repair of SHM-induced DNA breaks. A varying degree of antibody formation defect, from IgA or selective IgG subclass deficiency to common variable immunodeficiency and hyper-IgM syndrome, has been detected in a small number of patients with constitutional mismatch repair deficiency (CMMRD) due to biallelic loss-of-function mutations in one of the MMR genes (PMS2, MSH6, MLH1, or MSH2). To elucidate the clinical relevance of a presumed primary immunodeficiency (PID) in CMMRD, we systematically collected clinical history and laboratory data of a cohort of 15 consecutive, unrelated patients (10 not previously reported) with homozygous/compound heterozygous mutations in PMS2 (n = 8), MSH6 (n = 5), and MLH1 (n = 2), most of whom manifested with typical malignancies during childhood. Detailed descriptions of their genotypes, phenotypes, and family histories are provided. Importantly, none of the patients showed any clinical warning signs of PID (infections, immune dysregulation, inflammation, failure to thrive, etc.). Furthermore, we could not detect uniform or specific patterns of laboratory abnormalities. The concentration of IgM was increased in 3 out of 12, reduced in 3 out of 12, and normal in 6 out of 12 patients, while concentrations of IgG and IgG subclasses, except IgG4, and of IgA, and specific antibody formation were normal in most. Class-switched B memory cells were reduced in 5 out of 12 patients, and in 9 out of 12 also the CD38^hiIgM⁻ plasmablasts were reduced. Furthermore, results of next generation sequencing-based analyses of antigen-selected B-cell receptor rearrangements showed a significantly reduced frequency of SHM and an increased number of rearranged immunoglobulin heavy chain (IGH) transcripts that use IGHG3, IGHG1, and IGHA1 subclasses. T cell subsets and receptor repertoires were unaffected. Together, neither clinical nor routine immunological laboratory parameters were consistently suggestive of PID in these CMMRD patients, but previously shown abnormalities in SHM and rearranged heavy chain transcripts were confirmed.

Genetic analyses in a cohort of 191 pulmonary arterial hypertension patients

Background

Pulmonary arterial hypertension (PAH) is a progressive and fatal disorder associated with high pulmonary artery pressure. Genetic testing enables early diagnosis and offers an opportunity for family screening. To identify genetic mutations and help make a precise diagnosis, we performed genetic testing in 191 probands with PAH and tried to analyze the genotype-phenotype correlation.

Methods

Initially, PAH samples (n = 119) were submitted to BMPR2 screening using Sanger sequencing. Later, we developed a PAH panel test to identify causal mutations in 13 genes related to PAH and tried to call BMPR2 copy number variations (CNVs) with the panel data. Multiplex ligation-dependent probe amplification (MLPA) was used to search for CNVs in BMPR2, ACVRL1 and ENG. Notably, EIF2AK4 gene was also involved in the panel, which allowed to distinguish pulmonary veno-occlusive disease (PVOD)/pulmonary capillary hemangiomatosis (PCH) patients from idiopathic PAH (IPAH). Characteristics of patients were compared using t test for continuous variables.

Results

Pathogenic BMPR2 mutations were detected most frequently in 32 (17.9%) IPAH and 5 (41.7%) heritable PAH (HPAH) patients by sequencing, and 12 BMPR2 CNVs called from the panel data were all successfully confirmed by MLPA analysis. In addition, homozygous or compound heterozygous EIF2AK4 mutations were identified in 6 patients, who should be corrected to a diagnosis of PVOD/PCH. Genotype-phenotype correlation analysis revealed that PAH patients with BMPR2 mutations were younger at diagnosis (27.2y vs. 31.6y, p = 0.0003) and exhibited more severe pulmonary hemodynamic impairment and a worse cardiac index compared with those without BMPR2 mutations.

Conclusions

The panel assay represented a highly valuable tool in PAH genetic testing, not only for the detection of small sequence alterations, but also for an indication of BMPR2 CNVs, which had implications for the specific samples to perform further MLPA assay. Analyses of PAH causal genes have a great help to clinical diagnosis and deep implications in disease treatment.

Electronic supplementary material

The online version of this article (10.1186/s12931-018-0789-9) contains supplementary material, which is available to authorized users.

Targeted Next-generation Sequencing and Bioinformatics Pipeline to Evaluate Genetic Determinants of Constitutional Disease

Next-generation sequencing (NGS) is quickly revolutionizing how research into the genetic determinants of constitutional disease is performed. The technique is highly efficient with millions of sequencing reads being produced in a short time span and at relatively low cost. Specifically, targeted NGS is able to focus investigations to genomic regions of particular interest based on the disease of study. Not only does this further reduce costs and increase the speed of the process, but it lessens the computational burden that often accompanies NGS. Although targeted NGS is restricted to certain regions of the genome, preventing identification of potential novel loci of interest, it can be an excellent technique when faced with a phenotypically and genetically heterogeneous disease, for which there are previously known genetic associations. Because of the complex nature of the sequencing technique, it is important to closely adhere to protocols and methodologies in order to achieve sequencing reads of high coverage and quality. Further, once sequencing reads are obtained, a sophisticated bioinformatics workflow is utilized to accurately map reads to a reference genome, to call variants, and to ensure the variants pass quality metrics. Variants must also be annotated and curated based on their clinical significance, which can be standardized by applying the American College of Medical Genetics and Genomics Pathogenicity Guidelines. The methods presented herein will display the steps involved in generating and analyzing NGS data from a targeted sequencing panel, using the ONDRISeq neurodegenerative disease panel as a model, to identify variants that may be of clinical significance.

The ICR96 exon CNV validation series: a resource for orthogonal assessment of exon CNV calling in NGS data

Detection of deletions and duplications of whole exons (exon CNVs) is a key requirement of genetic testing. Accurate detection of this variant type has proved very challenging in targeted next-generation sequencing (NGS) data, particularly if only a single exon is involved. Many different NGS exon CNV calling methods have been developed over the last five years. Such methods are usually evaluated using simulated and/or in-house data due to a lack of publicly-available datasets with orthogonally generated results. This hinders tool comparisons, transparency and reproducibility. To provide a community resource for assessment of exon CNV calling methods in targeted NGS data, we here present the ICR96 exon CNV validation series. The dataset includes high-quality sequencing data from a targeted NGS assay (the TruSight Cancer Panel) together with Multiplex Ligation-dependent Probe Amplification (MLPA) results for 96 independent samples. 66 samples contain at least one validated exon CNV and 30 samples have validated negative results for exon CNVs in 26 genes. The dataset includes 46 exon CNVs in BRCA1, BRCA2, TP53, MLH1, MSH2, MSH6, PMS2, EPCAM or PTEN, giving excellent representation of the cancer predisposition genes most frequently tested in clinical practice. Moreover, the validated exon CNVs include 25 single exon CNVs, the most difficult type of exon CNV to detect. The FASTQ files for the ICR96 exon CNV validation series can be accessed through the European-Genome phenome Archive (EGA) under the accession number EGAS00001002428.

panelcn.MOPS: Copy-number detection in targeted NGS panel data for clinical diagnostics

Targeted next‐generation‐sequencing (NGS) panels have largely replaced Sanger sequencing in clinical diagnostics. They allow for the detection of copy‐number variations (CNVs) in addition to single‐nucleotide variants and small insertions/deletions. However, existing computational CNV detection methods have shortcomings regarding accuracy, quality control (QC), incidental findings, and user‐friendliness. We developed panelcn.MOPS, a novel pipeline for detecting CNVs in targeted NGS panel data. Using data from 180 samples, we compared panelcn.MOPS with five state‐of‐the‐art methods. With panelcn.MOPS leading the field, most methods achieved comparably high accuracy. panelcn.MOPS reliably detected CNVs ranging in size from part of a region of interest (ROI), to whole genes, which may comprise all ROIs investigated in a given sample. The latter is enabled by analyzing reads from all ROIs of the panel, but presenting results exclusively for user‐selected genes, thus avoiding incidental findings. Additionally, panelcn.MOPS offers QC criteria not only for samples, but also for individual ROIs within a sample, which increases the confidence in called CNVs. panelcn.MOPS is freely available both as R package and standalone software with graphical user interface that is easy to use for clinical geneticists without any programming experience. panelcn.MOPS combines high sensitivity and specificity with user‐friendliness rendering it highly suitable for routine clinical diagnostics.