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.

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.

High-Throughput Amplicon-Based Copy Number Detection of 11 Genes in Formalin-Fixed Paraffin-Embedded Ovarian Tumour Samples by MLPA-Seq.

Whilst next generation sequencing can report point mutations in fixed tissue tumour samples reliably, the accurate determination of copy number is more challenging. The conventional Multiplex Ligation-dependent Probe Amplification (MLPA) assay is an effective tool for measurement of gene dosage, but is restricted to around 50 targets due to size resolution of the MLPA probes. By switching from a size-resolved format, to a sequence-resolved format we developed a scalable, high-throughput, quantitative assay. MLPA-seq is capable of detecting deletions, duplications, and amplifications in as little as 5ng of genomic DNA, including from formalin-fixed paraffin-embedded (FFPE) tumour samples. We show that this method can detect BRCA1, BRCA2, ERBB2 and CCNE1 copy number changes in DNA extracted from snap-frozen and FFPE tumour tissue, with 100% sensitivity and >99.5% specificity.

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.

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.

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.

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.

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.

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.

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.

Mutations in NMNAT1 cause Leber congenital amaurosis and identify a new disease pathway for retinal degeneration.

Leber congenital amaurosis (LCA) is a blinding retinal disease that presents within the first year after birth. Using exome sequencing, we identified mutations in the nicotinamide adenine dinucleotide (NAD) synthase gene NMNAT1 encoding nicotinamide mononucleotide adenylyltransferase 1 in eight families with LCA, including the family in which LCA was originally linked to the LCA9 locus. Notably, all individuals with NMNAT1 mutations also have macular colobomas, which are severe degenerative entities of the central retina (fovea) devoid of tissue and photoreceptors. Functional assays of the proteins encoded by the mutant alleles identified in our study showed that the mutations reduce the enzymatic activity of NMNAT1 in NAD biosynthesis and affect protein folding. Of note, recent characterization of the slow Wallerian degeneration (Wld(s)) mouse model, in which prolonged axonal survival after injury is observed, identified NMNAT1 as a neuroprotective protein when ectopically expressed. Our findings identify a new disease mechanism underlying LCA and provide the first link between endogenous NMNAT1 dysfunction and a human nervous system disorder.

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.

Next-generation sequencing for simultaneous determination of human papillomavirus load, subtype, and associated genomic copy number changes in tumors.

Human papillomavirus (HPV) infection in cases of squamous cell carcinoma of the oropharynx is a powerful predictive and prognostic biomarker. We describe how the use of next-generation sequencing can provide a novel method for the detection of HPV in DNA isolated from formalin-fixed paraffin-embedded tissues. Using this methodology in a cohort of 44 head and neck tumors, we identified the samples that contained HPV sequences, the viral subtype involved, and a direct readout of viral load. Specificity of HPV detection by sequencing compared to traditional detection methods using either PCR or p16 immunohistochemistry was 100%. Sensitivity was 50% when either compared to PCR [confidence interval (CI) = 29% to 71%] or 75% when compared to p16 (CI = 47% to 91%). In addition, we demonstrate the ability of next-generation sequencing to detect other HPV subtypes that would not have been detected by traditional methods, and we demonstrated the ability to apply this method to any tumor and any virus in a panel of eight human cancer cell lines. This methodology also provides a tumor genomic copy number karyogram, and in the samples analyzed here, a lower level of chromosome instability was detected in HPV-positive tumors compared to HPV-negative tumors, as observed in previous studies. Thus, the use of next-generation sequencing for the detection of HPV provides a multiplicity of data with clinical significance in a single test.

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 in MEGF10, a regulator of satellite cell myogenesis, cause early onset myopathy, areflexia, respiratory distress and dysphagia (EMARDD).

Infantile myopathies with diaphragmatic paralysis are genetically heterogeneous, and clinical symptoms do not assist in differentiating between them. We used phased haplotype analysis with subsequent targeted exome sequencing to identify MEGF10 mutations in a previously unidentified type of infantile myopathy with diaphragmatic weakness, areflexia, respiratory distress and dysphagia. MEGF10 is highly expressed in activated satellite cells and regulates their proliferation as well as their differentiation and fusion into multinucleated myofibers, which are greatly reduced in muscle from individuals with early onset myopathy, areflexia, respiratory distress and dysphagia.

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.

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.