Genetic diagnosis of familial breast cancer using clonal sequencing

Assessing the utility of long-read nanopore sequencing for rapid and efficient characterization of mobile element insertions

Short-read next generation sequencing (NGS) has become the predominant first-line technique used to diagnose patients with rare genetic conditions. Inherent limitations of short-read technology, notably for the detection and characterization of complex insertion-containing variants, are offset by the ability to concurrently screen many disease genes. "Third-generation" long-read sequencers are increasingly being deployed as an orthogonal adjunct technology, but their full potential for molecular genetic diagnosis has yet to be exploited. Here, we describe three diagnostic cases in which pathogenic mobile element insertions were refractory to characterization by short-read sequencing. To validate the accuracy of the long-read technology, we first used Sanger sequencing to confirm the integration sites and derive curated benchmark sequences of the variant-containing alleles. Long-read nanopore sequencing was then performed on locus-specific amplicons. Pairwise comparison between these data and the previously determined benchmark alleles revealed 100% identity of the variant-containing sequences. We demonstrate a number of technical advantages over existing wet-laboratory approaches, including in silico size selection of a mixed pool of amplification products, and the relative ease with which an automated informatics workflow can be established. Our findings add to a growing body of literature describing the diagnostic utility of long-read sequencing.

Simultaneous detection of BRCA mutations and large genomic rearrangements in germline DNA and FFPE tumor samples

The development of breast and ovarian cancer is strongly connected to the inactivation of the BRCA1 and BRCA2 genes by different germline and somatic alterations, and their diagnosis has great significance in targeted tumor therapy, since recently approved PARP inhibitors show high efficiency in the treatment of BRCA-deficient tumors. This raises the need for new diagnostic methods that are capable of performing an integrative mutation analysis of the BRCA genes not only from germline DNA but also from formalin-fixed and paraffin-embedded (FFPE) tumor samples. Here we describe the development of such a methodology based on next-generation sequencing and a new bioinformatics software for data analysis. The diagnostic method was initially developed on an Illumina MiSeq NGS platform using germline-mutated stem cell lines and then adapted for the Ion Torrent PGM NGS platform as well. We also investigated the usability of NGS coverage data for the detection of copy number variations and exon deletions as a replacement of the conventional MLPA technique. Finally, we tested the developed workflow on FFPE samples from breast and ovarian cancer patients. Our method meets the sensitivity and specificity requirements for the genetic diagnosis of breast and ovarian cancers both from germline and FFPE samples.

Characterization and Genomic Localization of a SMAD4 Processed Pseudogene

Like many clinical diagnostic laboratories, the Yorkshire Regional Genetics Service undertakes routine investigation of cancer-predisposed individuals by high-throughput sequencing of patient DNA that has been target-enriched for genes associated with hereditary cancer. Accurate diagnosis using such reagents requires alertness regarding rare nonpathogenic variants that may interfere with variant calling. In a cohort of 2042 such cases, we identified 5 that initially appeared to be carriers of a 95-bp deletion of SMAD4 intron 6. More detailed analysis indicated that these individuals all carried one copy of a SMAD4 processed gene. Because of its interference with diagnostic analysis, we characterized this processed gene in detail. Whole-genome sequencing and confirmatory Sanger sequencing of junction PCR products were used to show that in each of the 5 cases, the SMAD4 processed gene was integrated at the same position on chromosome 9, located within the last intron of the SCAI gene. This rare polymorphic processed gene therefore reflects the occurrence of a single ancestral retrotransposition event. Compared to the reference SMAD4 mRNA sequence NM_005359.5 (https://www.ncbi.nlm.nih.gov/nucleotide), the 5' and 3' untranslated regions of the processed gene are both truncated, but its open reading frame is unaltered. Our experience leads us to advocate the use of an RNA-seq aligner as part of diagnostic assay quality assurance, since this allows recognition of processed pseudogenes in a comparatively facile automated fashion.

New challenges for BRCA testing: a view from the diagnostic laboratory

Increased demand for BRCA testing is placing pressures on diagnostic laboratories to raise their mutation screening capacity and handle the challenges associated with classifying BRCA sequence variants for clinical significance, for example interpretation of pathogenic mutations or variants of unknown significance, accurate determination of large genomic rearrangements and detection of somatic mutations in DNA extracted from formalin-fixed, paraffin-embedded tumour samples. Many diagnostic laboratories are adopting next-generation sequencing (NGS) technology to increase their screening capacity and reduce processing time and unit costs. However, migration to NGS introduces complexities arising from choice of components of the BRCA testing workflow, such as NGS platform, enrichment method and bioinformatics analysis process. An efficient, cost-effective accurate mutation detection strategy and a standardised, systematic approach to the reporting of BRCA test results is imperative for diagnostic laboratories. This review covers the challenges of BRCA testing from the perspective of a diagnostics laboratory.

Clinical Applications of Next-Generation Sequencing in Cancer Diagnosis

With the advancement and improvement of new sequencing technology, next-generation sequencing (NGS) has been applied increasingly in cancer genomics research fields. More recently, NGS has been adopted in clinical oncology to advance personalized treatment of cancer. NGS is utilized to novel diagnostic and rare cancer mutations, detection of translocations, inversions, insertions and deletions, detection of copy number variants, detect familial cancer mutation carriers, provide the molecular rationale for appropriate targeted, therapeutic and prognostic. NGS holds many advantages, such as the ability to fully sequence all types of mutations for a large number of genes (hundreds to thousands) and the sensitivity, speed in a single test at a relatively low cost compared to be other sequencing modalities. Here we described the technology, methods and applications that can be immediately considered and some of the challenges that lie ahead.

Evaluation of an amplicon-based next-generation sequencing panel for detection of BRCA1 and BRCA2 genetic variants

The recent advances in the next-generation sequencing (NGS) technology have enabled fast, accurate, and cost-effective genetic testing. Here, we evaluated the performance of a targeted NGS panel for BRCA1/2 sequencing and confirmed its applicability in routine clinical diagnostics. We tested samples from 88 patients using the TruSeq custom panel (Illumina Inc, USA) and a MiSeq sequencer (Illumina) and compared the results to the outcomes of conventional Sanger sequencing. All 1015 sequence variations identified by Sanger sequencing were detected by NGS, except for one missense variant that might have been missed due to a rare mutation on a primer-binding site. One deletion variation, c.1909 + 12delT of BRCA2, was falsely called in all samples due to a homopolymer error. In addition, seven different single-nucleotide substitutions with low variant frequencies (range: 16.2-33.3 %) were falsely called by NGS. In a separate batch, 10 different false-positive variations were found in five samples. The overall sensitivity and positive predictive value of NGS were estimated to be 99.9 and 87.5 %, respectively. The false-positive results could be excluded by setting quality and alternative allele ratio filters and/or by visual inspection using the IGV software. Targeted NGS panel for BRCA1 and BRCA2 showed an excellent agreement with Sanger sequencing results. We therefore conclude that this NGS panel can be used for routine diagnostic method in a clinical genetic laboratory.

A Chromosome 7 Pericentric Inversion Defined at Single-Nucleotide Resolution Using Diagnostic Whole Genome Sequencing in a Patient with Hand-Foot-Genital Syndrome

Next generation sequencing methodologies are facilitating the rapid characterisation of novel structural variants at nucleotide resolution. These approaches are particularly applicable to variants initially identified using alternative molecular methods. We report a child born with bilateral postaxial syndactyly of the feet and bilateral fifth finger clinodactyly. This was presumed to be an autosomal recessive syndrome, due to the family history of consanguinity. Karyotype analysis revealed a homozygous pericentric inversion of chromosome 7 (46,XX,inv(7)(p15q21)x2) which was confirmed to be heterozygous in both unaffected parents. Since the resolution of the karyotype was insufficient to identify any putatively causative gene, we undertook medium-coverage whole genome sequencing using paired-end reads, in order to elucidate the molecular breakpoints. In a two-step analysis, we first narrowed down the region by identifying discordant read-pairs, and then determined the precise molecular breakpoint by analysing the mapping locations of “soft-clipped” breakpoint-spanning reads. PCR and Sanger sequencing confirmed the identified breakpoints, both of which were located in intergenic regions. Significantly, the 7p15 breakpoint was located 523 kb upstream of HOXA13, the locus for hand-foot-genital syndrome. By inference from studies of HOXA locus control in the mouse, we suggest that the inversion has delocalised a HOXA13 enhancer to produce the phenotype observed in our patient. This study demonstrates how modern genetic diagnostic approach can characterise structural variants at nucleotide resolution and provide potential insights into functional regulation.

Using Genetics to Identify Hereditary Colorectal Polyposis and Cancer Syndromes in Your Patient

The majority of patients with colorectal polyps and cancer do not have a Mendelian cause of the disease. Age, lifestyle, and environmental factors interact with complex genetic traits to contribute to the etiology. However, approximately 5-10 % of patients with colorectal cancer (CRC) and more than 40 % of patients meeting specific clinical features of the hereditary polyposis syndromes have a discoverable, actionable genetic cause which will significantly alter their medical management.

Enhanced diagnostic yield in Meckel-Gruber and Joubert syndrome through exome sequencing supplemented with split-read mapping

Background

The widespread adoption of high-throughput sequencing technologies by genetic diagnostic laboratories has enabled significant expansion of their testing portfolios. Rare autosomal recessive conditions have been a particular focus of many new services. Here we report a cohort of 26 patients referred for genetic analysis of Joubert (JBTS) and Meckel-Gruber (MKS) syndromes, two clinically and genetically heterogeneous neurodevelopmental conditions that define a phenotypic spectrum, with MKS at the severe end.

Methods

Exome sequencing was performed for all cases, using Agilent SureSelect v5 reagents and Illumina paired-end sequencing. For two cases medium-coverage (9×) whole genome sequencing was subsequently undertaken.

Results

Using a standard analysis pipeline for the detection of single nucleotide and small insertion or deletion variants, molecular diagnoses were confirmed in 12 cases (4 %). Seeking to determine whether our cohort harboured pathogenic copy number variants (CNV), in JBTS- or MKS-associated genes, targeted comparative read-depth analysis was performed using FishingCNV. These analyses identified a putative intragenic AHI1 deletion that included three exons spanning at least 3.4 kb and an intergenic MPP4 to TMEM237 deletion that included exons spanning at least 21.5 kb. Whole genome sequencing enabled confirmation of the deletion-containing alleles and precise characterisation of the mutation breakpoints at nucleotide resolution. These data were validated following development of PCR-based assays that could be subsequently used for “cascade” screening and/or prenatal diagnosis.

Conclusions

Our investigations expand the AHI1 and TMEM237 mutation spectrum and highlight the importance of performing CNV screening of disease-associated genes. We demonstrate a robust increasingly cost-effective CNV detection workflow that is applicable to all MKS/JBTS referrals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-015-0265-z) contains supplementary material, which is available to authorized users.

Diagnostic Approach to Hereditary Colorectal Cancer Syndromes

Approximately 5 to 10% of colorectal cancers develop within a known hereditary syndrome. Specific underlying genetic mutations drive the clinical phenotype and it is imperative to determine the genetic etiology to provide meaningful surveillance and intervention. Recognizing potential patients and families with a hereditary predisposition is the first step in management. Syndromes can be categorized according to polyp burden as polyposis or nonpolyposis. Clinical assessment should start with a personal and family medical history, physical examination, and evaluation for the presence and type of colorectal polyps or cancers. Key information is gained from these simple steps and should guide the specific genetic analysis for diagnosis. Genetic counseling is a critical component to any hereditary colorectal cancer program and should be conducted before genetic testing to provide education about the implications of test results. This review focuses on the thought process that drives initial clinical evaluation and guides genetic testing for patients with suspected hereditary colorectal cancer syndromes.

Addressing health disparities in Hispanic breast cancer: accurate and inexpensive sequencing of BRCA1 and BRCA2

Background

Germline mutations in the BRCA1 and BRCA2 genes account for 20–25 % of inherited breast cancers and about 10 % of all breast cancer cases. Detection of BRCA mutation carriers can lead to therapeutic interventions such as mastectomy, oophorectomy, hormonal prevention therapy, improved screening, and targeted therapies such as PARP-inhibition. We estimate that African Americans and Hispanics are 4–5 times less likely to receive BRCA screening, despite having similar mutation frequencies as non-Jewish Caucasians, who have higher breast cancer mortality. To begin addressing this health disparity, we initiated a nationwide trial of BRCA testing of Latin American women with breast cancer. Patients were recruited through community organizations, clinics, public events, and by mail and Internet. Subjects completed the consent process and questionnaire, and provided a saliva sample by mail or in person. DNA from 120 subjects was used to sequence the entirety of BRCA1 and BRCA2 coding regions and splice sites, and validate pathogenic mutations, with a total material cost of $85/subject. Subjects ranged in age from 23 to 81 years (mean age, 51 years), 6 % had bilateral disease, 57 % were ER/PR+, 23 % HER2+, and 17 % had triple-negative disease.

Results

A total of seven different predicted deleterious mutations were identified, one newly described and the rest rare. In addition, four variants of unknown effect were found.

Conclusions

Application of this strategy on a larger scale could lead to improved cancer care of minority and underserved populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0088-z) contains supplementary material, which is available to authorized users.

Towards a Next-Generation Sequencing Diagnostic Service for Tumour Genotyping: A Comparison of Panels and Platforms

Detection of clinically actionable mutations in diagnostic tumour specimens aids in the selection of targeted therapeutics. With an ever increasing number of clinically significant mutations identified, tumour genetic diagnostics is moving from single to multigene analysis. As it is still not feasible for routine diagnostic laboratories to perform sequencing of the entire cancer genome, our approach was to undertake targeted mutation detection. To optimise our diagnostic workflow, we evaluated three target enrichment strategies using two next-generation sequencing (NGS) platforms (Illumina MiSeq and Ion PGM). The target enrichment strategies were Fluidigm Access Array custom amplicon panel including 13 genes (MiSeq sequencing), the Oxford Gene Technologies (OGT) SureSeq Solid Tumour hybridisation panel including 60 genes (MiSeq sequencing), and an Ion AmpliSeq Cancer Hotspot Panel including 50 genes (Ion PGM sequencing). DNA extracted from formalin-fixed paraffin-embedded (FFPE) blocks of eight previously characterised cancer cell lines was tested using the three panels. Matching genomic DNA from fresh cultures of these cell lines was also tested using the custom Fluidigm panel and the OGT SureSeq Solid Tumour panel. Each panel allowed mutation detection of core cancer genes including KRAS, BRAF, and EGFR. Our results indicate that the panels enable accurate variant detection despite sequencing from FFPE DNA.

Genome, Exome, and Targeted Next-Generation Sequencing in Neonatal Diabetes

The use of targeted gene panels now allows the analysis of all the genes known to cause a disease in a single test. For neonatal diabetes, this has resulted in a paradigm shift with patients receiving a genetic diagnosis early and the genetic results guiding their clinical management. Exome and genome sequencing are powerful tools to identify novel genetic causes of known diseases. For neonatal diabetes, the use of these technologies has resulted in the identification of 2 novel disease genes (GATA6 and STAT3) and a novel regulatory element of PTF1A, in which mutations cause pancreatic agenesis.

Clinical applications of next generation sequencing in cancer: from panels, to exomes, to genomes

This article will review recent impact of massively parallel next-generation sequencing (NGS) in our understanding and treatment of cancer. While whole exome sequencing (WES) remains popular and effective as a method of genetically profiling different cancers, advances in sequencing technology has enabled an increasing number of whole-genome based studies. Clinically, NGS has been used or is being developed for genetic screening, diagnostics, and clinical assessment. Though challenges remain, clinicians are in the early stages of using genetic data to make treatment decisions for cancer patients. As the integration of NGS in the study and treatment of cancer continues to mature, we believe that the field of cancer genomics will need to move toward more complete 100% genome sequencing. Current technologies and methods are largely limited to coding regions of the genome. A number of recent studies have demonstrated that mutations in non-coding regions may have direct tumorigenic effects or lead to genetic instability. Non-coding regions represent an important frontier in cancer genomics.

Next-generation Sequencing of RYR1 and CACNA1S in Malignant Hyperthermia and Exertional Heat Illness

BACKGROUND

Variants in RYR1 are associated with the majority of cases of malignant hyperthermia (MH), a form of heat illness pharmacogenetically triggered by general anesthetics, and they have also been associated with exertional heat illness (EHI). CACNA1S has also been implicated in MH. The authors applied a targeted next-generation sequencing approach to identify variants in RYR1 and CACNA1S in a cohort of unrelated patients diagnosed with MH susceptibility. They also provide the first comprehensive report of sequencing of these two genes in a cohort of survivors of EHI.

METHODS

DNA extracted from blood was genotyped using a "long" polymerase chain reaction technique, with sequencing on the Illumina GAII or MiSeq platforms (Illumina Inc., USA). Variants were assessed for pathogenicity using bioinformatic approaches. For further follow-up, DNA from additional family members and up to 211 MH normal and 556 MH-susceptible unrelated individuals was tested.

RESULTS

In 29 MH patients, the authors identified three pathogenic and four novel RYR1 variants, with a further five RYR1 variants previously reported in association with MH. Three novel RYR1 variants were found in the EHI cohort (n = 28) along with two more previously reported in association with MH. Two other variants were reported previously associated with centronuclear myopathy. The authors found one and three rare variants of unknown significance in CACNA1S in the MH and EHI cohorts, respectively.

CONCLUSIONS

Targeted next-generation sequencing proved efficient at identifying diagnostically useful and potentially implicated variants in RYR1 and CACNA1S in MH and EHI.

Robust diagnostic genetic testing using solution capture enrichment and a novel variant-filtering interface

Targeted hybridization enrichment prior to next-generation sequencing is a widespread method for characterizing sequence variation in a research setting, and is being adopted by diagnostic laboratories. However, the number of variants identified can overwhelm clinical laboratories with strict time constraints, the final interpretation of likely pathogenicity being a particular bottleneck. To address this, we have developed an approach in which, after automatic variant calling on a standard unix pipeline, subsequent variant filtering is performed interactively, using AgileExomeFilter and AgilePindelFilter (http://dna.leeds.ac.uk/agile), tools designed for clinical scientists with standard desktop computers. To demonstrate the method's diagnostic efficacy, we tested 128 patients using (1) a targeted capture of 36 cancer-predisposing genes or (2) whole-exome capture for diagnosis of the genetically heterogeneous disorder primary ciliary dyskinesia (PCD). In the cancer cohort, complete concordance with previous diagnostic data was achieved across 793 variant genotypes. A high yield (42%) was also achieved for exome-based PCD diagnosis, underscoring the scalability of our method. Simple adjustments to the variant filtering parameters further allowed the identification of a homozygous truncating mutation in a presumptive new PCD gene, DNAH8. These tools should allow diagnostic laboratories to expand their testing portfolios flexibly, using a standard set of reagents and techniques.

A high-plex PCR approach for massively parallel sequencing

Current methods for targeted massively parallel sequencing (MPS) have several drawbacks, including limited design flexibility, expense, and protocol complexity, which restrict their application to settings involving modest target size and requiring low cost and high throughput. To address this, we have developed Hi-Plex, a PCR-MPS strategy intended for high-throughput screening of multiple genomic target regions that integrates simple, automated primer design software to control product size. Featuring permissive thermocycling conditions and clamp bias reduction, our protocol is simple, cost- and time-effective, uses readily available reagents, does not require expensive instrumentation, and requires minimal optimization. In a 60-plex assay targeting the breast cancer predisposition genes PALB2 and XRCC2, we applied Hi-Plex to 100 ng LCL-derived DNA, and 100 ng and 25 ng FFPE tumor-derived DNA. Altogether, at least 86.94% of the human genome-mapped reads were on target, and 100% of targeted amplicons were represented within 25-fold of the mean. Using 25 ng FFPE-derived DNA, 95.14% of mapped reads were on-target and relative representation ranged from 10.1-fold lower to 5.8-fold higher than the mean. These results were obtained using only the initial automatically-designed primers present in equal concentration. Hi-Plex represents a powerful new approach for screening panels of genomic target regions.

Next-generation sequencing of the BRCA1 and BRCA2 genes for the genetic diagnostics of hereditary breast and/or ovarian cancer

Genetic testing for hereditary breast and/or ovarian cancer mostly relies on laborious molecular tools that use Sanger sequencing to scan for mutations in the BRCA1 and BRCA2 genes. We explored a more efficient genetic screening strategy based on next-generation sequencing of the BRCA1 and BRCA2 genes in 210 hereditary breast and/or ovarian cancer patients. We first validated this approach in a cohort of 115 samples with previously known BRCA1 and BRCA2 mutations and polymorphisms. Genomic DNA was amplified using the Ion AmpliSeq BRCA1 and BRCA2 panel. The DNA Libraries were pooled, barcoded, and sequenced using an Ion Torrent Personal Genome Machine sequencer. The combination of different robust bioinformatics tools allowed detection of all previously known pathogenic mutations and polymorphisms in the 115 samples, without detecting spurious pathogenic calls. We then used the same assay in a discovery cohort of 95 uncharacterized hereditary breast and/or ovarian cancer patients for BRCA1 and BRCA2. In addition, we describe the allelic frequencies across 210 hereditary breast and/or ovarian cancer patients of 74 unique definitely and likely pathogenic and uncertain BRCA1 and BRCA2 variants, some of which have not been previously annotated in the public databases. Targeted next-generation sequencing is ready to substitute classic molecular methods to perform genetic testing on the BRCA1 and BRCA2 genes and provides a greater opportunity for more comprehensive testing of at-risk patients.

Genetic testing in hereditary breast and ovarian cancer using massive parallel sequencing

High throughput methods such as next generation sequencing are increasingly used in molecular diagnosis. The aim of this study was to develop a workflow for the detection of BRCA1 and BRCA2 mutations using massive parallel sequencing in a 454 GS Junior bench top sequencer. Our approach was first validated in a panel of 23 patients containing 62 unique variants that had been previously Sanger sequenced. Subsequently, 101 patients with familial breast and ovarian cancer were studied. BRCA1 and BRCA2 exon enrichment has been performed by PCR amplification using the BRCA MASTR kit (Multiplicom). Bioinformatic analysis of reads is performed with the AVA software v2.7 (Roche). In total, all 62 variants were detected resulting in a sensitivity of 100%. 71 false positives were called resulting in a specificity of 97.35%. All of them correspond to deletions located in homopolymeric stretches. The analysis of the homopolymers stretches of 6 bp or longer using the BRCA HP kit (Multiplicom) increased the specificity of the detection of BRCA1 and BRCA2 mutations to 99.99%. We show here that massive parallel pyrosequencing can be used as a diagnostic strategy to test for BRCA1 and BRCA2 mutations meeting very stringent sensitivity and specificity parameters replacing traditional Sanger sequencing with a lower cost.

Utilization of multigene panels in hereditary cancer predisposition testing: analysis of more than 2,000 patients

PURPOSE

The aim of this study was to determine the clinical and molecular characteristics of 2,079 patients who underwent hereditary cancer multigene panel testing.

METHODS

Panels included comprehensive analysis of 14-22 cancer susceptibility genes (BRCA1 and BRCA2 not included), depending on the panel ordered (BreastNext, OvaNext, ColoNext, or CancerNext). Next-generation sequencing and deletion/duplication analyses were performed for all genes except EPCAM (deletion/duplication analysis only). Clinical histories of ColoNext patients harboring mutations in genes with well-established diagnostic criteria were assessed to determine whether diagnostic/testing criteria were met.

RESULTS

Positive rates were defined as the proportion of patients with a pathogenic mutation/likely pathogenic variant(s) and were as follows: 7.4% for BreastNext, 7.2% for OvaNext, 9.2% for ColoNext, and 9.6% for CancerNext. Inconclusive results were found in 19.8% of BreastNext, 25.6% of OvaNext, 15.1% of ColoNext, and 23.5% of CancerNext tests. Based on information submitted by clinicians, 30% of ColoNext patients with mutations in genes with well-established diagnostic criteria did not meet corresponding criteria.

CONCLUSION

Our data point to an important role for targeted multigene panels in diagnosing hereditary cancer predisposition, particularly for patients with clinical histories spanning several possible diagnoses and for patients with suspicious clinical histories not meeting diagnostic criteria for a specific hereditary cancer syndrome.

Next-generation sequencing identifies germline MRE11A variants as markers of radiotherapy outcomes in muscle-invasive bladder cancer

BACKGROUND

Muscle-invasive bladder cancer (MIBC) can be cured by radical radiotherapy (RT). We previously found tumour MRE11 expression to be predictive of survival following RT in MIBC, and this was independently validated in a separate institute. Here, we investigated germline MRE11A variants as possible predictors of RT outcomes in MIBC, using next-generation sequencing (NGS).

PATIENTS AND METHODS

The MRE11A gene was amplified in germline DNA from 186 prospectively recruited MIBC patients treated with RT and sequenced using bar-coded multiplexed NGS. Germline variants were analysed for associations with cancer-specific survival (CSS). For validation as a prognostic or predictive marker, rs1805363 was then genotyped in a cystectomy-treated MIBC cohort of 256 individuals. MRE11A mRNA isoform expression was measured in bladder cancer cell lines and primary tumour samples.

RESULTS

Carriage of at least one of six (five novel) rare variants was associated with the worse RT outcome (hazard ratio [HR] 4.04, 95% confidence interval [95% CI] 1.42-11.51, P = 0.009). The single-nucleotide polymorphism (SNP), rs1805363 (minor allele frequency 11%), was also associated with worse CSS (per-allele HR 2.10, 95% CI 1.34-3.28, Ptrend = 0.001) following RT in MIBC, with a gene-dosage effect observed, but no effect seen on CSS in the cystectomy cohort (Ptrend = 0.89). Furthermore, rs1805363 influenced relative MRE11A isoform expression, with increased isoform 2 expression with carriage of the rs1805363 minor A allele.

CONCLUSIONS

Germline MRE11A SNP rs1805363 was predictive of RT, but not of cystectomy outcome in MIBC. If successfully validated in an independent RT-treated cohort, this SNP could be a useful clinical tool for selecting patients for bladder-conserving treatment.

ROVER variant caller: read-pair overlap considerate variant-calling software applied to PCR-based massively parallel sequencing datasets

Background

We recently described Hi-Plex, a highly multiplexed PCR-based target-enrichment system for massively parallel sequencing (MPS), which allows the uniform definition of library size so that subsequent paired-end sequencing can achieve complete overlap of read pairs. Variant calling from Hi-Plex-derived datasets can thus rely on the identification of variants appearing in both reads of read-pairs, permitting stringent filtering of sequencing chemistry-induced errors. These principles underly ROVER software (derived from Read Overlap PCR-MPS variant caller), which we have recently used to report the screening for genetic mutations in the breast cancer predisposition gene PALB2. Here, we describe the algorithms underlying ROVER and its usage.

Results

ROVER enables users to quickly and accurately identify genetic variants from PCR-targeted, overlapping paired-end MPS datasets. The open-source availability of the software and threshold tailorability enables broad access for a range of PCR-MPS users.

Methods

ROVER is implemented in Python and runs on all popular POSIX-like operating systems (Linux, OS X). The software accepts a tab-delimited text file listing the coordinates of the target-specific primers used for targeted enrichment based on a specified genome-build. It also accepts aligned sequence files resulting from mapping to the same genome-build. ROVER identifies the amplicon a given read-pair represents and removes the primer sequences by using the mapping co-ordinates and primer co-ordinates. It considers overlapping read-pairs with respect to primer-intervening sequence. Only when a variant is observed in both reads of a read-pair does the signal contribute to a tally of read-pairs containing or not containing the variant. A user-defined threshold informs the minimum number of, and proportion of, read-pairs a variant must be observed in for a ‘call’ to be made. ROVER also reports the depth of coverage across amplicons to facilitate the identification of any regions that may require further screening.

Conclusions

ROVER can facilitate rapid and accurate genetic variant calling for a broad range of PCR-MPS users.

A massive parallel sequencing workflow for diagnostic genetic testing of mismatch repair genes

The purpose of this study was to develop a massive parallel sequencing (MPS) workflow for diagnostic analysis of mismatch repair (MMR) genes using the GS Junior system (Roche). A pathogenic variant in one of four MMR genes, (MLH1, PMS2, MSH6, and MSH2), is the cause of Lynch Syndrome (LS), which mainly predispose to colorectal cancer. We used an amplicon-based sequencing method allowing specific and preferential amplification of the MMR genes including PMS2, of which several pseudogenes exist. The amplicons were pooled at different ratios to obtain coverage uniformity and maximize the throughput of a single-GS Junior run. In total, 60 previously identified and distinct variants (substitutions and indels), were sequenced by MPS and successfully detected. The heterozygote detection range was from 19% to 63% and dependent on sequence context and coverage. We were able to distinguish between false-positive and true-positive calls in homopolymeric regions by cross-sample comparison and evaluation of flow signal distributions. In addition, we filtered variants according to a predefined status, which facilitated variant annotation. Our study shows that implementation of MPS in routine diagnostics of LS can accelerate sample throughput and reduce costs without compromising sensitivity, compared to Sanger sequencing.

A comprehensive next generation sequencing-based genetic testing strategy to improve diagnosis of inherited pheochromocytoma and paraganglioma

CONTEXT

Pheochromocytomas and paragangliomas are notable for a high frequency of inherited cases, many of which present as apparently sporadic tumors.

OBJECTIVE

The objective of this study was to establish a comprehensive next generation sequencing (NGS)-based strategy for the diagnosis of patients with pheochromocytoma and paraganglioma by testing simultaneously for mutations in MAX, RET, SDHA, SDHB, SDHC, SDHD, SDHAF2, TMEM127, and VHL.

DESIGN

After the methodology for the assay was designed and established, it was validated on DNA samples with known genotype and then patients were studied prospectively.

SETTING

The study was performed in a diagnostic genetics laboratory.

PATIENTS

DNA samples from 205 individuals affected with adrenal or extraadrenal pheochromocytoma/head and neck paraganglioma (PPGL/HNPGL) were analyzed. A proof-of-principle study was performed using 85 samples known to contain a variant in 1 or more of the genes to be tested, followed by prospective analysis of an additional 120 samples.

MAIN OUTCOME MEASURES

We assessed the ability to use an NGS-based method to perform comprehensive analysis of genes implicated in inherited PPGL/HNPGL.

RESULTS

The proof-of-principle study showed that the NGS assay and analysis gave a sensitivity of 98.7%. A pathogenic mutation was identified in 16.6% of the prospective analysis cohort of 120 patients.

CONCLUSIONS

A comprehensive NGS-based strategy for the analysis of genes associated with predisposition to PPGL and HNPGL was established, validated, and introduced into diagnostic service. The new assay provides simultaneous analysis of 9 genes and allows more rapid and cost-effective mutation detection than the previously used conventional Sanger sequencing-based methodology.

Systematic detection of pathogenic alu element insertions in NGS-based diagnostic screens: the BRCA1/BRCA2 example

Pathogenic Alu element insertions are rarely reported, whereas their occurrence is expected to be much higher. Alu containing alleles are usually out-competed during the PCR process and consequently undetectable with the classical screening methods. However, with the introduction of the next generation sequencing (NGS) technology in the diagnostic field, new opportunities are emerging. NGS data for a particular genomic region can be seen as the summation of all the individual sequences (reads) obtained for that region and no longer as the mean of this sum as it is the case for traditional Sanger sequencing. Because each single read covering that region is expected to be generated from a different template molecule, the presence of one single mutant read must theoretically be sufficient to identify the mutation. However, generation and identification of mutant reads bearing Alu insertions remains challenging and several wet/dry bench parameters need to be optimized. Hereby we present the proof of principle of a NGS-based mutation screening procedure allowing the detection of inherited Alu insertions within any predefined sequence by investigating 2 cases: c.1739_1740insAlu in BRCA1 and c.156_157insAlu in BRCA2.

Application of massively parallel sequencing to genetic diagnosis in multiplex families with idiopathic sensorineural hearing impairment

Despite the clinical utility of genetic diagnosis to address idiopathic sensorineural hearing impairment (SNHI), the current strategy for screening mutations via Sanger sequencing suffers from the limitation that only a limited number of DNA fragments associated with common deafness mutations can be genotyped. Consequently, a definitive genetic diagnosis cannot be achieved in many families with discernible family history. To investigate the diagnostic utility of massively parallel sequencing (MPS), we applied the MPS technique to 12 multiplex families with idiopathic SNHI in which common deafness mutations had previously been ruled out. NimbleGen sequence capture array was designed to target all protein coding sequences (CDSs) and 100 bp of the flanking sequence of 80 common deafness genes. We performed MPS on the Illumina HiSeq2000, and applied BWA, SAMtools, Picard, GATK, Variant Tools, ANNOVAR, and IGV for bioinformatics analyses. Initial data filtering with allele frequencies (<5% in the 1000 Genomes Project and 5400 NHLBI exomes) and PolyPhen2/SIFT scores (>0.95) prioritized 5 indels (insertions/deletions) and 36 missense variants in the 12 multiplex families. After further validation by Sanger sequencing, segregation pattern, and evolutionary conservation of amino acid residues, we identified 4 variants in 4 different genes, which might lead to SNHI in 4 families compatible with autosomal dominant inheritance. These included GJB2 p.R75Q, MYO7A p.T381M, KCNQ4 p.S680F, and MYH9 p.E1256K. Among them, KCNQ4 p.S680F and MYH9 p.E1256K were novel. In conclusion, MPS allows genetic diagnosis in multiplex families with idiopathic SNHI by detecting mutations in relatively uncommon deafness genes.

Nonoptical massive parallel DNA sequencing of BRCA1 and BRCA2 genes in a diagnostic setting

The introduction of the benchtop massive parallel sequencers made it possible for the majority of clinical diagnostic laboratories to gain access to this fast evolving technology. In this study, using the Ion Torrent Personal Genome Machine, we present a strategy for the molecular diagnosis of hereditary breast and ovarian cancer and respective analytical validation. The methodology relies on a multiplex PCR amplification of the BRCA1 and BRCA2 genes combined with a variant prioritization pipeline, designed to minimize the number of false-positive calls without the introduction of false-negative results. A training set of samples was used to optimize the entire process, and a second set was used to validate and independently evaluate the performance of the workflow. Performing the study in a blind manner relative to the variants in the samples and using conventional Sanger sequencing as standard, the workflow resulted in a strategy with a maximum analytical sensitivity ≥98.6% with a confidence of 95% and a specificity of 96.9%. Importantly, no true variant was missed. This study presents a comprehensive massive parallel sequencing-Sanger sequencing based strategy, which results in a high analytical sensitivity assay that provides a time- and cost-effective strategy for the identification of mutations in the BRCA1 and BRCA2 genes.

Simultaneous sequencing of 24 genes associated with steroid-resistant nephrotic syndrome

BACKGROUND AND OBJECTIVES

Up to 95% of children presenting with steroid-resistant nephrotic syndrome in early life will have a pathogenic single-gene mutation in 1 of 24 genes currently associated with this disease. Others may be affected by polymorphic variants. There is currently no accepted diagnostic algorithm for clinical genetic testing. The hypothesis was that the increasing reliability of next generation sequencing allows comprehensive one-step genetic investigation of this group and similar patient groups.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

This study used next generation sequencing to screen 446 genes, including the 24 genes known to be associated with hereditary steroid-resistant nephrotic syndrome. The first 36 pediatric patients collected through a national United Kingdom Renal Registry were chosen with comprehensive phenotypic detail. Significant variants detected by next generation sequencing were confirmed by conventional Sanger sequencing.

RESULTS

Analysis revealed known and novel disease-associated variations in expected genes such as NPHS1, NPHS2, and PLCe1 in 19% of patients. Phenotypically unexpected mutations were also detected in COQ2 and COL4A4 in two patients with isolated nephropathy and associated sensorineural deafness, respectively. The presence of an additional heterozygous polymorphism in WT1 in a patient with NPHS1 mutation was associated with earlier-onset disease, supporting modification of phenotype through genetic epistasis.

CONCLUSIONS

This study shows that next generation sequencing analysis of pediatric steroid-resistant nephrotic syndrome patients is accurate and revealing. This analysis should be considered part of the routine genetic workup of diseases such as childhood steroid-resistant nephrotic syndrome, where the chance of genetic mutation is high but requires sequencing of multiple genes.

Next-generation sequencing meets genetic diagnostics: development of a comprehensive workflow for the analysis of BRCA1 and BRCA2 genes

Next-generation sequencing (NGS) is changing genetic diagnosis due to its huge sequencing capacity and cost-effectiveness. The aim of this study was to develop an NGS-based workflow for routine diagnostics for hereditary breast and ovarian cancer syndrome (HBOCS), to improve genetic testing for BRCA1 and BRCA2. A NGS-based workflow was designed using BRCA MASTR kit amplicon libraries followed by GS Junior pyrosequencing. Data analysis combined Variant Identification Pipeline freely available software and ad hoc R scripts, including a cascade of filters to generate coverage and variant calling reports. A BRCA homopolymer assay was performed in parallel. A research scheme was designed in two parts. A Training Set of 28 DNA samples containing 23 unique pathogenic mutations and 213 other variants (33 unique) was used. The workflow was validated in a set of 14 samples from HBOCS families in parallel with the current diagnostic workflow (Validation Set). The NGS-based workflow developed permitted the identification of all pathogenic mutations and genetic variants, including those located in or close to homopolymers. The use of NGS for detecting copy-number alterations was also investigated. The workflow meets the sensitivity and specificity requirements for the genetic diagnosis of HBOCS and improves on the cost-effectiveness of current approaches.

Integrating next-generation sequencing into the diagnostic testing of inherited cancer predisposition

The clinical application of next-generation sequencing (NGS) as a diagnostic tool has become increasingly evident. The coupling of NGS technologies with new genomic sequence enrichment methods has made the sequencing of panels of target genes technically feasible, at the same time as making such an approach cost-effective for diagnostic applications. In this article, we discuss recent studies that have applied NGS in the diagnostic setting in relation to hereditary cancer.

Mutation detection by clonal sequencing of PCR amplicons and grouped read typing is applicable to clinical diagnostics

We describe a sensitive technique for mutation detection using clonal sequencing. We analyzed DNA extracted from 13 cancer cell lines and 35 tumor samples and applied a novel approach to identify disease-associated somatic mutations. By matching reads against an index of known variants, noise can be dramatically reduced, enabling the detection and quantification of those variants, even when they are present at less than 1% of the total sequenced population; this is comparable to, or better than, current diagnostic methods. Following the identification or exclusion of known variants, unmatched reads are grouped for BLAST searching to identify novel variants or contaminants. Known variants, novel variants, and contaminants were readily identified in tumor tissue using this approach. Our approach also enables an estimation of the per-base sequencing error rate, providing a confidence threshold for interpretation of the results in the clinic. This novel approach has immediate applicability to clinical testing for disease-associated genetic variants.

Development of a next-generation sequencing method for BRCA mutation screening: a comparison between a high-throughput and a benchtop platform

In a clinical setting, next-generation sequencing (NGS) approaches for the enrichment and resequencing of DNA targets may have limitations in throughput, cost, or accuracy. We evaluated an NGS workflow for targeted DNA sequencing for mutation detection. Targeted sequence data of the BRCA1 and BRCA2 genes, generated using a PCR-based, multiplexed NGS approach using the SOLiD 4 (n = 24) and Ion Torrent PGM (n = 20) next-generation sequencers, were evaluated against sequence data obtained by Sanger sequencing. The overall sensitivity for SOLiD and PGM were 97.8% (95% CI = 94.7 to 100.0) and 98.9% (95% CI = 96.8 to 100.0) respectively. The specificity for the SOLiD platform was high, at 100.0% (95% CI = 99.3 to 100.0). PGM correctly identified all 3 indels, but 68 false-positive indels were also called. Equimolar normalization of amplicons was not necessary for successful NGS. Both platforms are highly amenable to scale-up, potentially reducing the reagent cost for BRCA testing to <US$200. Only 325 ng of DNA per patient is required, with similar coverage and accuracy obtained using DNA derived from peripheral blood or buccal wash samples. The strategy described is accurate and easy to incorporate into conventional workflow, and shows potential for mutation screening of clinically important gene targets in genetic disorders.

Long-range PCR and next-generation sequencing of BRCA1 and BRCA2 in breast cancer

Individuals and families carrying mutations in BRCA1 and BRCA2 (BRCA1/2) have a markedly elevated risk of developing breast and ovarian cancers. The first-generation of BRCA1/2 mutation analysis targeted only the coding exons and has implicated protein-truncating mutations (indel, nonsense) in BRCA1/2 inactivation. Recently, heritable breast cancers have also been attributed to other exonic mutations (missense, silent) and mutations in introns and untranslated regions. However, analysis of these alterations has been prohibitively laborious and cost intensive, and the proportion of cases carrying mutations in unscreened regions of BRCA1/2 and other predisposition genes is unknown. We have developed and validated a next-generation sequencing (NGS) approach for BRCA1/2 mutation analysis by applying long-range PCR and deep sequencing. Genomic DNA from familial breast cancer patients (N = 12) were screened and NGS successfully identified all 19 distinct (51 total) BRCA1 and 35 distinct (63 total) BRCA2 sequence alterations detectable by the Sanger sequencing, with no false-negative or positive results. In addition, we report the robust detection of variants from introns and untranslated regions. These results illustrate that NGS can provide comprehensive genetic information more quickly, accurately, and at a lower cost than conventional approaches, and we propose NGS to be a more effective method for BRCA1/2 mutational analysis. Advances in NGS will play an important role in enabling molecular diagnostics and personalized treatment of breast and ovarian cancers.

ColoSeq provides comprehensive lynch and polyposis syndrome mutational analysis using massively parallel sequencing

Lynch syndrome (hereditary nonpolyposis colon cancer) and adenomatous polyposis syndromes frequently have overlapping clinical features. Current approaches for molecular genetic testing are often stepwise, taking a best-candidate gene approach with testing of additional genes if initial results are negative. We report a comprehensive assay called ColoSeq that detects all classes of mutations in Lynch and polyposis syndrome genes using targeted capture and massively parallel next-generation sequencing on the Illumina HiSeq2000 instrument. In blinded specimens and colon cancer cell lines with defined mutations, ColoSeq correctly identified 28/28 (100%) pathogenic mutations in MLH1, MSH2, MSH6, PMS2, EPCAM, APC, and MUTYH, including single nucleotide variants (SNVs), small insertions and deletions, and large copy number variants. There was 100% reproducibility of detection mutation between independent runs. The assay correctly identified 222 of 224 heterozygous SNVs (99.4%) in HapMap samples, demonstrating high sensitivity of calling all variants across each captured gene. Average coverage was greater than 320 reads per base pair when the maximum of 96 index samples with barcodes were pooled. In a specificity study of 19 control patients without cancer from different ethnic backgrounds, we did not find any pathogenic mutations but detected two variants of uncertain significance. ColoSeq offers a powerful, cost-effective means of genetic testing for Lynch and polyposis syndromes that eliminates the need for stepwise testing and multiple follow-up clinical visits.

From morphologic to molecular: established and emerging molecular diagnostics for breast carcinoma

Diagnostics in the field of breast carcinoma are constantly evolving. The recent wave of molecular methodologies, both microscope and non-microscope based, have opened new ways to gain insight into this disease process and have moved clinical diagnostics closer to a 'personalized medicine' approach. In this review we highlight some of the advancements that laboratory medicine technology is making toward guiding the diagnosis, prognosis, and therapy selection for patients affected by breast carcinoma. The content of the article is largely structured by methodology, with a distinct emphasis on both microscope based and non-microscope based diagnostic formats. Where possible, we have attempted to emphasize the potential benefits as well as limitations to each of these technologies. Successful molecular diagnostics, applied in concert within the morphologic context of a patient's tumor, are what will lay the foundation for personalized therapy and allow a more sophisticated approach to clinical trial stratification. The future of breast cancer diagnostics looks challenging, but it is also a field of great opportunity. Never before have there been such a plethora of new tools available for disease investigation or candidate therapy selection.

Massively parallel sequencing for early molecular diagnosis in Leber congenital amaurosis

PURPOSE

Leber congenital amaurosis (LCA) is a rare congenital retinal dystrophy associated with 16 genes. Recent breakthroughs in LCA gene therapy offer the first prospect of treating inherited blindness, which requires an unequivocal and early molecular diagnosis. While present genetic tests do not address this due to a tremendous genetic heterogeneity, massively parallel sequencing (MPS) strategies might bring a solution. Here, we developed a comprehensive molecular test for LCA based on targeted MPS of all exons of 16 known LCA genes.

METHODS

We designed a unique and flexible workflow for targeted resequencing of all 236 exons from 16 LCA genes based on quantitative PCR (qPCR) amplicon ligation, shearing, and parallel sequencing of multiple patients on a single lane of a short-read sequencer. Twenty-two prescreened LCA patients were included, five of whom had a known molecular cause.

RESULTS

Validation of 107 variations was performed as proof of concept. In addition, the causal genetic defect and a single heterozygous mutation were identified in 3 and 5, respectively, of 17 patients without previously identified mutations.

CONCLUSION

We propose a novel targeted MPS-based approach that is suitable for accurate, fast, and cost-effective early molecular testing in LCA, and easily applicable in other genetic disorders.

Applications of second generation sequencing technologies in complex disorders

Second generation sequencing (2ndGS) technologies generate unprecedented amounts of sequence data very rapidly and at relatively limited costs, allowing the sequence of a human genome to be completed in a few weeks. The principle is on the basis of generating millions of relatively short reads from amplified single DNA fragments using iterative cycles of nucleotide extensions. However, the data generated on this scale present new challenges in interpretation, data analysis and data management. 2ndGS technologies are becoming widespread and are profoundly impacting biomedical research. Common applications include whole-genome sequencing, target resequencing, characterization of structural and copy number variation, profiling epigenetic modifications, transcriptome sequencing and identification of infectious agents. New methodologies and instruments that will enable to sequence the complete human genome in less than a day at a cost of less than $1,000 are currently in development.

A novel application of pattern recognition for accurate SNP and indel discovery from high-throughput data: targeted resequencing of the glucocorticoid receptor co-chaperone FKBP5 in a Caucasian population

The detection of single nucleotide polymorphisms (SNPs) and insertion/deletions (indels) with precision from high-throughput data remains a significant bioinformatics challenge. Accurate detection is necessary before next-generation sequencing can routinely be used in the clinic. In research, scientific advances are inhibited by gaps in data, exemplified by the underrepresented discovery of rare variants, variants in non-coding regions and indels. The continued presence of false positives and false negatives prevents full automation and requires additional manual verification steps. Our methodology presents applications of both pattern recognition and sensitivity analysis to eliminate false positives and aid in the detection of SNP/indel loci and genotypes from high-throughput data. We chose FK506-binding protein 51(FKBP5) (6p21.31) for our clinical target because of its role in modulating pharmacological responses to physiological and synthetic glucocorticoids and because of the complexity of the genomic region. We detected genetic variation across a 160 kb region encompassing FKBP5. 613 SNPs and 57 indels, including a 3.3 kb deletion were discovered. We validated our method using three independent data sets and, with Sanger sequencing and Affymetrix and Illumina microarrays, achieved 99% concordance. Furthermore we were able to detect 267 novel rare variants and assess linkage disequilibrium. Our results showed both a sensitivity and specificity of 98%, indicating near perfect classification between true and false variants. The process is scalable and amenable to automation, with the downstream filters taking only 1.5h to analyze 96 individuals simultaneously. We provide examples of how our level of precision uncovered the interactions of multiple loci, their predicted influences on mRNA stability, perturbations of the hsp90 binding site, and individual variation in FKBP5 expression. Finally we show how our discovery of rare variants may change current conceptions of evolution at this locus.

Review of massively parallel DNA sequencing technologies

Since the development of technologies that can determine the base-pair sequence of DNA, the ability to sequence genes has contributed much to science and medicine. However, it has remained a relatively costly and laborious process, hindering its use as a routine biomedical tool. Recent times are seeing rapid developments in this field, both in the availability of novel sequencing platforms, as well as supporting technologies involved in processes such as targeting and data analysis. This is leading to significant reductions in the cost of sequencing a human genome and the potential for its use as a routine biomedical tool. This review is a snapshot of this rapidly moving field examining the current state of the art, forthcoming developments and some of the issues still to be resolved prior to the use of new sequencing technologies in routine clinical diagnosis.