Exploring the feasibility and utility of exome-scale tumour sequencing in a clinical setting.

BACKGROUND: Technology has progressed from single gene panel to large-scale genomic sequencing. This is raising expectations from clinicians and patients alike. The utility and performance of this technology in a clinical setting needs to be evaluated. AIM: This pilot study investigated the feasibility of using exome-scale sequencing (ESS) to identify molecular drivers within cancers in real-time for Precision Oncology in the clinic. METHODS: Between March 2014 and March 2015, the Victorian Comprehensive Cancer Centre Alliance explored the feasibility and utility of ESS in a pilot study. DNA extracted from the tumour specimens underwent both ESS and targeted 'hotspot' sequencing (TS). Blood was taken for germline analysis. A multi-disciplinary molecular tumour board determined the clinical relevance of identified mutations; in particular, whether they were 'actionable' and/or 'druggable'. RESULTS: Of 23 patients screened, 15 (65%) met the tissue requirements for genomic analysis. TS and ESS were successful in all cases. ESS identified pathogenic somatic variants in 73% (11/15 cases) versus 53% (8/15 cases) using TS. Clinically focused ESS identified 63 variants, consisting of 30 somatic variants (including all 13 identified by TS) and 33 germline variants. Overall, there were 48 unique variants. ESS had a clinical impact in 53% (8/15 cases); 47% (7/15 cases) were referred to the familial cancer clinic, and 'druggable' targets were identified in 53% (8/15 cases). CONCLUSION: ESS of tumour DNA impacted clinical decision-making in 53%, with 20% more pathogenic variants identified through ESS than TS. The identification of germline variants in 47% was an unexpected finding.

BRAF Inhibition in BRAFV600E-Positive Anaplastic Thyroid Carcinoma.

The efficacy of targeted monotherapy for BRAF(V600E)-positive anaplastic thyroid carcinomas (ATC) is not established. We report 2 cases of BRAF(V600E)-positive ATC treated with a BRAF inhibitor. A 49-year-old woman with a T4bN1bM0 ATC manifested symptomatic metastatic disease 8 weeks after radical chemoradiotherapy. Within 1 month of BRAF inhibitor monotherapy, a complete symptomatic response was observed, with FDG-PET scan confirming metabolic and radiologic response. Treatment was terminated after 3 months because of disease progression. The patient died 11 months after primary diagnosis. A 67-year-old man received first-line BRAF inhibitor for a T4aN1bM0 ATC. Within 10 days of treatment his pain had stabilized and his tumor had clinically halved in size. Stable disease was achieved for 11 weeks but the patient died 11 months after diagnosis because of disease progression. BRAF inhibitor monotherapy in ATC may obtain clinical benefit of short duration. Upfront combination therapy should be investigated in this patient subgroup.

AmpliVar: mutation detection in high-throughput sequence from amplicon-based libraries.

Conventional means of identifying variants in high-throughput sequencing align each read against a reference sequence, and then call variants at each position. Here, we demonstrate an orthogonal means of identifying sequence variation by grouping the reads as amplicons prior to any alignment. We used AmpliVar to make key-value hashes of sequence reads and group reads as individual amplicons using a table of flanking sequences. Low-abundance reads were removed according to a selectable threshold, and reads above this threshold were aligned as groups, rather than as individual reads, permitting the use of sensitive alignment tools. We show that this approach is more sensitive, more specific, and more computationally efficient than comparable methods for the analysis of amplicon-based high-throughput sequencing data. The method can be extended to enable alignment-free confirmation of variants seen in hybridization capture target-enrichment data.

Mutations in C4orf26, encoding a peptide with in vitro hydroxyapatite crystal nucleation and growth activity, cause amelogenesis imperfecta.

Autozygosity mapping and clonal sequencing of an Omani family identified mutations in the uncharacterized gene, C4orf26, as a cause of recessive hypomineralized amelogenesis imperfecta (AI), a disease in which the formation of tooth enamel fails. Screening of a panel of 57 autosomal-recessive AI-affected families identified eight further families with loss-of-function mutations in C4orf26. C4orf26 encodes a putative extracellular matrix acidic phosphoprotein expressed in the enamel organ. A mineral nucleation assay showed that the protein's phosphorylated C terminus has the capacity to promote nucleation of hydroxyapatite, suggesting a possible function in enamel mineralization during amelogenesis.

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.

Next generation sequencing identifies mutations in Atonal homolog 7 (ATOH7) in families with global eye developmental defects.

The atonal homolog 7 (ATOH7) gene encodes a transcription factor involved in determining the fate of retinal progenitor cells and is particularly required for optic nerve and ganglion cell development. Using a combination of autozygosity mapping and next generation sequencing, we have identified homozygous mutations in this gene, p.E49V and p.P18RfsX69, in two consanguineous families diagnosed with multiple ocular developmental defects, including severe vitreoretinal dysplasia, optic nerve hypoplasia, persistent fetal vasculature, microphthalmia, congenital cataracts, microcornea, corneal opacity and nystagmus. Most of these clinical features overlap with defects in the Norrin/ß-catenin signalling pathway that is characterized by dysgenesis of the retinal and hyaloid vasculature. Our findings document Mendelian mutations within ATOH7 and imply a role for this molecule in the development of structures at the front as well as the back of the eye. This work also provides further insights into the function of ATOH7, especially its importance in retinal vascular development and hyaloid regression.

Mutations causing familial biparental hydatidiform mole implicate c6orf221 as a possible regulator of genomic imprinting in the human oocyte.

Familial biparental hydatidiform mole (FBHM) is the only known pure maternal-effect recessive inherited disorder in humans. Affected women, although developmentally normal themselves, suffer repeated pregnancy loss because of the development of the conceptus into a complete hydatidiform mole in which extraembryonic trophoblastic tissue develops but the embryo itself suffers early demise. This developmental phenotype results from a genome-wide failure to correctly specify or maintain a maternal epigenotype at imprinted loci. Most cases of FBHM result from mutations of NLRP7, but genetic heterogeneity has been demonstrated. Here, we report biallelic mutations of C6orf221 in three families with FBHM. The previously described biological properties of their respective gene families suggest that NLRP7 and C6orf221 may interact as components of an oocyte complex that is directly or indirectly required for determination of epigenetic status on the oocyte genome.

Homozygous mutations in PXDN cause congenital cataract, corneal opacity, and developmental glaucoma.

Anterior segment dysgenesis describes a group of heterogeneous developmental disorders that affect the anterior chamber of the eye and are associated with an increased risk of glaucoma. Here, we report homozygous mutations in peroxidasin (PXDN) in two consanguineous Pakistani families with congenital cataract-microcornea with mild to moderate corneal opacity and in a consanguineous Cambodian family with developmental glaucoma and severe corneal opacification. These results highlight the diverse ocular phenotypes caused by PXDN mutations, which are likely due to differences in genetic background and environmental factors. Peroxidasin is an extracellular matrix-associated protein with peroxidase catalytic activity, and we confirmed localization of the protein to the cornea and lens epithelial layers. Our findings imply that peroxidasin is essential for normal development of the anterior chamber of the eye, where it may have a structural role in supporting cornea and lens architecture as well as an enzymatic role as an antioxidant enzyme in protecting the lens, trabecular meshwork, and cornea against oxidative damage.

Accurately identifying low-allelic fraction variants in single samples with next-generation sequencing: applications in tumor subclone resolution.

Current methods for resolving genetically distinct subclones in tumor samples require somatic mutations to be clustered by allelic frequencies, which are determined by applying a variant calling program to next-generation sequencing data. Such programs were developed to accurately distinguish true polymorphisms and somatic mutations from the artifactual nonreference alleles introduced during library preparation and sequencing. However, numerous variant callers exist with no clear indication of the best performer for subclonal analysis, in which the accuracy of the assigned variant frequency is as important as correctly indicating whether the variant is present or not. Furthermore, sequencing depth (the number of times that a genomic position is sequenced) affects the ability to detect low-allelic fraction variants and accurately assign their allele frequencies. We created two synthetic sequencing datasets, and sequenced real KRAS amplicons, with variants spiked in at specific ratios, to assess which caller performs best in terms of both variant detection and assignment of allelic frequencies. We also assessed the sequencing depths required to detect low-allelic fraction variants. We found that VarScan2 performed best overall with sequencing depths of 100×, 250×, 500×, and 1,000× required to accurately identify variants present at 10%, 5%, 2.5%, and 1%, respectively.

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.

Simple and efficient identification of rare recessive pathologically important sequence variants from next generation exome sequence data.

Massively parallel ("next generation") DNA sequencing (NGS) has quickly become the method of choice for seeking pathogenic mutations in rare uncharacterized monogenic diseases. Typically, before DNA sequencing, protein-coding regions are enriched from patient genomic DNA, representing either the entire genome ("exome sequencing") or selected mapped candidate loci. Sequence variants, identified as differences between the patient's and the human genome reference sequences, are then filtered according to various quality parameters. Changes are screened against datasets of known polymorphisms, such as dbSNP and the 1000 Genomes Project, in the effort to narrow the list of candidate causative variants. An increasing number of commercial services now offer to both generate and align NGS data to a reference genome. This potentially allows small groups with limited computing infrastructure and informatics skills to utilize this technology. However, the capability to effectively filter and assess sequence variants is still an important bottleneck in the identification of deleterious sequence variants in both research and diagnostic settings. We have developed an approach to this problem comprising a user-friendly suite of programs that can interactively analyze, filter and screen data from enrichment-capture NGS data. These programs ("Agile Suite") are particularly suitable for small-scale gene discovery or for diagnostic analysis.

Prostaglandin transporter mutations cause pachydermoperiostosis with myelofibrosis.

Pachydermoperiostosis, or primary hypertrophic osteoarthropathy (PHO), is an inherited multisystem disorder, whose features closely mimic the reactive osteoarthropathy that commonly accompanies neoplastic and inflammatory pathologies. We previously described deficiency of the prostaglandin-degrading enzyme 15-hydroxyprostaglandin dehydrogenase (HPGD) as a cause of this condition, implicating elevated circulating prostaglandin E(2) (PGE(2)) as causative of PHO, and perhaps also as the principal mediator of secondary HO. However, PHO is genetically heterogeneous. Here, we use whole-exome sequencing to identify recessive mutations of the prostaglandin transporter SLCO2A1, in individuals lacking HPGD mutations. We performed exome sequencing of four probands with severe PHO, followed by conventional mutation analysis of SLCO2A1 in nine others. Biallelic SLCO2A1 mutations were identified in 12 of the 13 families. Affected individuals had elevated urinary PGE(2), but unlike HPGD-deficient patients, also excreted considerable quantities of the PGE(2) metabolite, PGE-M. Clinical differences between the two groups were also identified, notably that SLCO2A1-deficient individuals have a high frequency of severe anemia due to myelofibrosis. These findings reinforce the key role of systemic or local prostaglandin excess as the stimulus to HO. They also suggest that the induction or maintenance of hematopoietic stem cells by prostaglandin may depend upon transporter activity.

Human Variome Project country nodes: documenting genetic information within a country.

The Human Variome Project (http://www.humanvariomeproject.org) is an international effort aiming to systematically collect and share information on all human genetic variation. The two main pillars of this effort are gene/disease-specific databases and a network of Human Variome Project Country Nodes. The latter are nationwide efforts to document the genomic variation reported within a specific population. The development and successful operation of the Human Variome Project Country Nodes are of utmost importance to the success of Human Variome Project's aims and goals because they not only allow the genetic burden of disease to be quantified in different countries, but also provide diagnosticians and researchers access to an up-to-date resource that will assist them in their daily clinical practice and biomedical research, respectively. Here, we report the discussions and recommendations that resulted from the inaugural meeting of the International Confederation of Countries Advisory Council, held on 12th December 2011, during the 2011 Human Variome Project Beijing Meeting. We discuss the steps necessary to maximize the impact of the Country Node effort for developing regional and country-specific clinical genetics resources and summarize a few well-coordinated genetic data collection initiatives that would serve as paradigms for similar projects.

Mutation of the variant alpha-tubulin TUBA8 results in polymicrogyria with optic nerve hypoplasia.

The critical importance of cytoskeletal function for correct neuronal migration during development of the cerebral cortex has been underscored by the identities of germline mutations underlying a number of human neurodevelopmental disorders. The proteins affected include TUBA1A, a major alpha-tubulin isoform, and microtubule-associated components such as doublecortin, and LIS1. Mutations in these genes are associated with the anatomical abnormality lissencephaly, which is believed to reflect failure of neuronal migration. An important recent observation has been the dependence of cortical neuronal migration upon acetylation of alpha-tubulin at lysine 40 by the histone acetyltransferase Elongator complex. Here, we describe a recognizable autosomal recessive syndrome, characterized by generalized polymicrogyria in association with optic nerve hypoplasia (PMGOH). By autozygosity mapping, we show that the molecular basis for this condition is mutation of the TUBA8 gene, encoding a variant alpha-tubulin of unknown function that is not susceptible to the lysine 40 acetylation that regulates microtubule function during cortical neuron migration. Together with the unique expression pattern of TUBA8 within the developing cerebral cortex, these observations suggest a role for this atypical microtubule component in regulating mammalian brain development.

GeneScreen: a program for high-throughput mutation detection in DNA sequence electropherograms.

BACKGROUND: While massively parallel DNA sequencing methods continue to evolve rapidly, the benchmark technique for detection and verification of rare (particularly disease-causing) sequence variants remains four-colour dye-terminator sequencing by capillary electrophoresis. The high throughput and long read lengths currently available have shifted the bottleneck in mutation detection away from data generation to data analysis. While excellent computational methods have been developed for quantifying sequence accuracy and detecting variants, either during de novo sequence assembly or for single-nucleotide polymorphism detection, the identification, verification and annotation of very rare sequence variants remains a rather labour-intensive process for which few software aids exist. AIM: To provide a freely available, intuitive software application for highly efficient mutation screening of large sequence batches. METHODS AND RESULTS: The authors developed GeneScreen, a desktop program that analyses capillary electropherograms and compares their sequences with a known reference for identification of mutations. The detected sequence variants are then made available for rapid assessment and annotation via a graphical user interface, allowing chosen variants to be exported for reporting and archiving. The program was validated using more than 16,000 diagnostic laboratory sequence traces. CONCLUSION: Using GeneScreen, a single user requires only a few minutes to identify rare mutations in hundreds of sequence traces, with comparable sensitivity to expensive commercial products.

Using next-generation sequencing for high resolution multiplex analysis of copy number variation from nanogram quantities of DNA from formalin-fixed paraffin-embedded specimens.

The use of next-generation sequencing technologies to produce genomic copy number data has recently been described. Most approaches, however, reply on optimal starting DNA, and are therefore unsuitable for the analysis of formalin-fixed paraffin-embedded (FFPE) samples, which largely precludes the analysis of many tumour series. We have sought to challenge the limits of this technique with regards to quality and quantity of starting material and the depth of sequencing required. We confirm that the technique can be used to interrogate DNA from cell lines, fresh frozen material and FFPE samples to assess copy number variation. We show that as little as 5 ng of DNA is needed to generate a copy number karyogram, and follow this up with data from a series of FFPE biopsies and surgical samples. We have used various levels of sample multiplexing to demonstrate the adjustable resolution of the methodology, depending on the number of samples and available resources. We also demonstrate reproducibility by use of replicate samples and comparison with microarray-based comparative genomic hybridization (aCGH) and digital PCR. This technique can be valuable in both the analysis of routine diagnostic samples and in examining large repositories of fixed archival material.

Illuminator, a desktop program for mutation detection using short-read clonal sequencing.

Current methods for sequencing clonal populations of DNA molecules yield several gigabases of data per day, typically comprising reads of < 100 nt. Such datasets permit widespread genome resequencing and transcriptome analysis or other quantitative tasks. However, this huge capacity can also be harnessed for the resequencing of smaller (gene-sized) target regions, through the simultaneous parallel analysis of multiple subjects, using sample "tagging" or "indexing". These methods promise to have a huge impact on diagnostic mutation analysis and candidate gene testing. Here we describe a software package developed for such studies, offering the ability to resolve pooled samples carrying barcode tags and to align reads to a reference sequence using a mutation-tolerant process. The program, Illuminator, can identify rare sequence variants, including insertions and deletions, and permits interactive data analysis on standard desktop computers. It facilitates the effective analysis of targeted clonal sequencer data without dedicated computational infrastructure or specialized training.

Genetic diagnosis of familial breast cancer using clonal sequencing.

Using conventional Sanger sequencing as a reference standard, we compared the sensitivity, specificity, and capacity of the Illumina GA II platform for the detection of TP53, BRCA1, and BRCA2 mutations in established tumor cell lines and DNA from patients with germline mutations. A total of 656 coding variants were identified in four cell lines and 65 patient DNAs. All of the known pathogenic mutations (including point mutations and insertions/deletions of up to 16 nucleotides) were identified, using a combination of the Illumina data analysis pipeline with custom and commercial sequence alignment software. In our configuration, clonal sequencing outperforms current diagnostic methods, providing a reduction in analysis times and in reagent costs compared with conventional sequencing. These improvements open the possibility of BRCA1/2 testing for a wider spectrum of at-risk women, and will allow the genetic classification of tumors prior to the use of novel PARP inhibitors to treat BRCA-deficient breast cancers.

Genetic and epigenetic analysis of recurrent hydatidiform mole.

Familial biparental hydatidiform mole (FBHM) is a maternal-effect autosomal recessive disorder in which recurrent pregnancy failure with molar degeneration occurs. The phenotype mimics molar pregnancy due to androgenesis, despite the normal genetic makeup of the conceptus. FBHM appears to result from a failure to establish correct maternal epigenetic identity at imprinted loci during oogenesis. Several women affected with FBHM have previously been shown to have biallelic mutations in the NLRP7 gene (NALP7). Here, we present the results of epigenetic and mutational analysis on FBHM patients from 11 families, 10 of them novel. We demonstrate a methylation defect at imprinted loci in tissue from four new FBHM cases. Biallelic NLRP7 mutations, including eight previously undescribed mutations, were found in all but one family. These results indicate for the first time that maternal imprints at some loci may be correctly specified in FBHM conceptions, since differential methylation of SGCE/PEG10 was preserved in all four cases.

Clinical Utility of Real-Time Targeted Molecular Profiling in the Clinical Management of Ovarian Cancer: The ALLOCATE Study.

PURPOSE: The ALLOCATE study was designed as a pilot to demonstrate the feasibility and clinical utility of real-time targeted molecular profiling of patients with recurrent or advanced ovarian cancer for identification of potential targeted therapies. PATIENTS AND METHODS: A total of 113 patients with ovarian cancer of varying histologies were recruited from two tertiary hospitals, with 99 patient cases suitable for prospective analysis. Targeted molecular and methylation profiling of fresh biopsy and archived tumor samples were performed by screening for mutations or copy-number variations in 44 genes and for promoter methylation of BRCA1 and RAD51C. RESULTS: Somatic genomic or methylation events were identified in 85% of all patient cases, with potentially actionable events with defined targeted therapies (including four resistance events) detected in 60% of all patient cases. On the basis of these findings, six patients received molecularly guided therapy, three patients had unsuspected germline cancer-associated BRCA1/2 mutations and were referred for genetic counseling, and two intermediate differentiated (grade 2) serous ovarian carcinomas were reclassified as low grade, leading to changes in clinical management. Additionally, secondary reversion mutations in BRCA1/2 were identified in fresh biopsy samples of two patients, consistent with clinical platinum/poly (ADP-ribose) polymerase inhibitor resistance. Timely reporting of results if molecular testing is done at disease recurrence, as well as early referral for patients with platinum-resistant cancers, were identified as factors that could improve the clinical utility of molecular profiling. CONCLUSION: ALLOCATE molecular profiling identified known genomic and methylation alterations of the different ovarian cancer subtypes and was deemed feasible and useful in routine clinical practice. Better patient selection and access to a wider range of targeted therapies or clinical trials will further enhance the clinical utility of molecular profiling.