Driver mutations in GNAQ and GNA11 genes as potential targets for precision immunotherapy in uveal melanoma patients.

Uveal melanoma (UM) is the most common ocular malignancy in adults. Nearly 95% of UM patients carry the mutually exclusive mutations in the homologous genes GNAQ (amino acid change Q209L/Q209P) and GNA11 (aminoacid change Q209L). UM is located in an immunosuppressed organ and does not suffer immunoediting. Therefore, we hypothesize that driver mutations in GNAQ/11 genes could be recognized by the immune system. Genomic and transcriptomic data from primary uveal tumors were collected from the TCGA-UM dataset (n = 80) and used to assess the immunogenic potential for GNAQ/GNA11 Q209L/Q209P mutations using a variety of tools and HLA type information. All prediction tools showed stronger GNAQ/11 Q209L binding to HLA than GNAQ/11 Q209P. The immunogenicity analysis revealed that Q209L is likely to be presented by more than 73% of individuals in 1000 G databases whereas Q209P is only predicted to be presented in 24% of individuals. GNAQ/11 Q209L showed a higher likelihood to be presented by HLA-I molecules than almost all driver mutations analyzed. Finally, samples carrying Q209L had a higher immune-reactive phenotype. Regarding cancer risk, seven HLA genotypes with low Q209L affinity show higher frequency in uveal melanoma patients than in the general population. However, no clear association was found between any HLA genotype and survival. Results suggest a high potential immunogenicity of the GNAQ/11 Q209L variant that could allow the generation of novel therapeutic tools to treat UM like neoantigen vaccinations.

Copy number architectures define treatment-mediated selection of lethal prostate cancer clones.

Despite initial responses to hormone treatment, metastatic prostate cancer invariably evolves to a lethal state. To characterize the intra-patient evolutionary relationships of metastases that evade treatment, we perform genome-wide copy number profiling and bespoke approaches targeting the androgen receptor (AR) on 167 metastatic regions from 11 organs harvested post-mortem from 10 men who died from prostate cancer. We identify diverse and patient-unique alterations clustering around the AR in metastases from every patient with evidence of independent acquisition of related genomic changes within an individual and, in some patients, the co-existence of AR-neutral clones. Using the genomic boundaries of pan-autosome copy number changes, we confirm a common clone of origin across metastases and diagnostic biopsies, and identified in individual patients, clusters of metastases occupied by dominant clones with diverged autosomal copy number alterations. These autosome-defined clusters are characterized by cluster-specific AR gene architectures, and in two index cases are topologically more congruent than by chance (p-values 3.07 × 10-8 and 6.4 × 10-4). Integration with anatomical sites suggests patterns of spread and points of genomic divergence. Here, we show that copy number boundaries identify treatment-selected clones with putatively distinct lethal trajectories.

Accumulation of copy number alterations and clinical progression across advanced prostate cancer.

BACKGROUND: Genomic copy number alterations commonly occur in prostate cancer and are one measure of genomic instability. The clinical implication of copy number change in advanced prostate cancer, which defines a wide spectrum of disease from high-risk localised to metastatic, is unknown. METHODS: We performed copy number profiling on 688 tumour regions from 300 patients, who presented with advanced prostate cancer prior to the start of long-term androgen deprivation therapy (ADT), in the control arm of the prospective randomised STAMPEDE trial. Patients were categorised into metastatic states as follows; high-risk non-metastatic with or without local lymph node involvement, or metastatic low/high volume. We followed up patients for a median of 7 years. Univariable and multivariable Cox survival models were fitted to estimate the association between the burden of copy number alteration as a continuous variable and the hazard of death or disease progression. RESULTS: The burden of copy number alterations positively associated with radiologically evident distant metastases at diagnosis (P=0.00006) and showed a non-linear relationship with clinical outcome on univariable and multivariable analysis, characterised by a sharp increase in the relative risk of progression (P=0.003) and death (P=0.045) for each unit increase, stabilising into more modest increases with higher copy number burdens. This association between copy number burden and outcome was similar in each metastatic state. Copy number loss occurred significantly more frequently than gain at the lowest copy number burden quartile (q=4.1 × 10-6). Loss of segments in chromosome 5q21-22 and gains at 8q21-24, respectively including CHD1 and cMYC occurred more frequently in cases with higher copy number alteration (for either region: Kolmogorov-Smirnov distance, 0.5; adjusted P<0.0001). Copy number alterations showed variability across tumour regions in the same prostate. This variance associated with increased risk of distant metastases (Kruskal-Wallis test P=0.037). CONCLUSIONS: Copy number alteration in advanced prostate cancer associates with increased risk of metastases at diagnosis. Accumulation of a limited number of copy number alterations associates with most of the increased risk of disease progression and death. The increased likelihood of involvement of specific segments in high copy number alteration burden cancers may suggest an order underlying the accumulation of copy number changes. TRIAL REGISTRATION: ClinicalTrials.gov NCT00268476 , registered on December 22, 2005. EudraCT  2004-000193-31 , registered on October 4, 2004.

The Immunogenic Potential of Recurrent Cancer Drug Resistance Mutations: An In Silico Study.

Cancer somatic mutations have been identified as a source of antigens that can be targeted by cancer immunotherapy. In this work, expanding on previous studies, we analyze the HLA-presentation properties of mutations that are known to drive resistance to cancer targeted-therapies. We survey a large dataset of mutations that confer resistance to different drugs and occur in numerous genes and tumor types. We show that a significant number of them are predicted in silico to be potentially immunogenic across a large proportion of the human population. Further, by analyzing a cohort of patients carrying a small subset of these resistance mutations, we provide evidence that what is observed in the general population may be indicative of the mutations' immunogenic potential in resistant patients. Two of the mutations in our dataset had previously been experimentally validated by others and it was confirmed that some of their associated neopeptides elicit T-cell responses in vitro. The identification of potent cancer-specific antigens can be instrumental for developing more effective immunotherapies. In this work, we propose a novel list of drug-resistance mutations, several of which are recurrent, that could be of particular interest in the context of off-the-shelf precision immunotherapies such as therapeutic cancer vaccines.

Clinical BRCA1/2 Reversion Analysis Identifies Hotspot Mutations and Predicted Neoantigens Associated with Therapy Resistance.

Reversion mutations in BRCA1 or BRCA2 are associated with resistance to PARP inhibitors and platinum. To better understand the nature of these mutations, we collated, codified, and analyzed more than 300 reversions. This identified reversion "hotspots" and "deserts" in regions encoding the N and C terminus, respectively, of BRCA2, suggesting that pathogenic mutations in these regions may be at higher or lower risk of reversion. Missense and splice-site pathogenic mutations in BRCA1/2 also appeared less likely to revert than truncating mutations. Most reversions were <100 bp deletions. Although many deletions exhibited microhomology, this was not universal, suggesting that multiple DNA-repair processes cause reversion. Finally, we found that many reversions were predicted to encode immunogenic neopeptides, suggesting a route to the treatment of reverted disease. As well as providing a freely available database for the collation of future reversion cases, these observations have implications for how drug resistance might be managed in BRCA-mutant cancers. SIGNIFICANCE: Reversion mutations in BRCA genes are a major cause of clinical platinum and PARP inhibitor resistance. This analysis of all reported clinical reversions suggests that the position of BRCA2 mutations affects the risk of reversion. Many reversions are also predicted to encode tumor neoantigens, providing a potential route to targeting resistance.This article is highlighted in the In This Issue feature, p. 1426.

Genome-wide plasma DNA methylation features of metastatic prostate cancer.

Tumor DNA circulates in the plasma of cancer patients admixed with DNA from noncancerous cells. The genomic landscape of plasma DNA has been characterized in metastatic castration-resistant prostate cancer (mCRPC) but the plasma methylome has not been extensively explored. Here, we performed next-generation sequencing (NGS) on plasma DNA with and without bisulfite treatment from mCRPC patients receiving either abiraterone or enzalutamide in the pre- or post-chemotherapy setting. Principal component analysis on the mCRPC plasma methylome indicated that the main contributor to methylation variance (principal component one, or PC1) was strongly correlated with genomically determined tumor fraction (r = -0.96; P < 10-8) and characterized by hypermethylation of targets of the polycomb repressor complex 2 components. Further deconvolution of the PC1 top-correlated segments revealed that these segments are comprised of methylation patterns specific to either prostate cancer or prostate normal epithelium. To extract information specific to an individual's cancer, we then focused on an orthogonal methylation signature, which revealed enrichment for androgen receptor binding sequences and hypomethylation of these segments associated with AR copy number gain. Individuals harboring this methylation pattern had a more aggressive clinical course. Plasma methylome analysis can accurately quantitate tumor fraction and identify distinct biologically relevant mCRPC phenotypes.

Extreme intratumour heterogeneity and driver evolution in mismatch repair deficient gastro-oesophageal cancer.

Mismatch repair deficient (dMMR) gastro-oesophageal adenocarcinomas (GOAs) show better outcomes than their MMR-proficient counterparts and high immunotherapy sensitivity. The hypermutator-phenotype of dMMR tumours theoretically enables high evolvability but their evolution has not been investigated. Here we apply multi-region exome sequencing (MSeq) to four treatment-naive dMMR GOAs. This reveals extreme intratumour heterogeneity (ITH), exceeding ITH in other cancer types >20-fold, but also long phylogenetic trunks which may explain the exquisite immunotherapy sensitivity of dMMR tumours. Subclonal driver mutations are common and parallel evolution occurs in RAS, PIK3CA, SWI/SNF-complex genes and in immune evasion regulators. MSeq data and evolution analysis of single region-data from 64 MSI GOAs show that chromosome 8 gains are early genetic events and that the hypermutator-phenotype remains active during progression. MSeq may be necessary for biomarker development in these heterogeneous cancers. Comparison with other MSeq-analysed tumour types reveals mutation rates and their timing to determine phylogenetic tree morphologies.

Author Correction: Extreme intratumour heterogeneity and driver evolution in mismatch repair deficient gastro-oesophageal cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Identification of single nucleotide variants using position-specific error estimation in deep sequencing data.

BACKGROUND: Targeted deep sequencing is a highly effective technology to identify known and novel single nucleotide variants (SNVs) with many applications in translational medicine, disease monitoring and cancer profiling. However, identification of SNVs using deep sequencing data is a challenging computational problem as different sequencing artifacts limit the analytical sensitivity of SNV detection, especially at low variant allele frequencies (VAFs). METHODS: To address the problem of relatively high noise levels in amplicon-based deep sequencing data (e.g. with the Ion AmpliSeq technology) in the context of SNV calling, we have developed a new bioinformatics tool called AmpliSolve. AmpliSolve uses a set of normal samples to model position-specific, strand-specific and nucleotide-specific background artifacts (noise), and deploys a Poisson model-based statistical framework for SNV detection. RESULTS: Our tests on both synthetic and real data indicate that AmpliSolve achieves a good trade-off between precision and sensitivity, even at VAF below 5% and as low as 1%. We further validate AmpliSolve by applying it to the detection of SNVs in 96 circulating tumor DNA samples at three clinically relevant genomic positions and compare the results to digital droplet PCR experiments. CONCLUSIONS: AmpliSolve is a new tool for in-silico estimation of background noise and for detection of low frequency SNVs in targeted deep sequencing data. Although AmpliSolve has been specifically designed for and tested on amplicon-based libraries sequenced with the Ion Torrent platform it can, in principle, be applied to other sequencing platforms as well. AmpliSolve is freely available at https://github.com/dkleftogi/AmpliSolve .

In Vivo Modeling of Chemoresistant Neuroblastoma Provides New Insights into Chemorefractory Disease and Metastasis.

Neuroblastoma is a pediatric cancer that is frequently metastatic and resistant to conventional treatment. In part, a lack of natively metastatic, chemoresistant in vivo models has limited our insight into the development of aggressive disease. The Th-MYCN genetically engineered mouse model develops rapidly progressive chemosensitive neuroblastoma and lacks clinically relevant metastases. To study tumor progression in a context more reflective of clinical therapy, we delivered multicycle treatment with cyclophosphamide to Th-MYCN mice, individualizing therapy using MRI, to generate the Th-MYCN CPM32 model. These mice developed chemoresistance and spontaneous bone marrow metastases. Tumors exhibited an altered immune microenvironment with increased stroma and tumor-associated fibroblasts. Analysis of copy number aberrations revealed genomic changes characteristic of human MYCN-amplified neuroblastoma, specifically copy number gains at mouse chromosome 11, syntenic with gains on human chromosome 17q. RNA sequencing revealed enriched expression of genes associated with 17q gain and upregulation of genes associated with high-risk neuroblastoma, such as the cell-cycle regulator cyclin B1-interacting protein 1 (Ccnb1ip1) and thymidine kinase (TK1). The antiapoptotic, prometastatic JAK-STAT3 pathway was activated in chemoresistant tumors, and treatment with the JAK1/JAK2 inhibitor CYT387 reduced progression of chemoresistant tumors and increased survival. Our results highlight that under treatment conditions that mimic chemotherapy in human patients, Th-MYCN mice develop genomic, microenvironmental, and clinical features reminiscent of human chemorefractory disease. The Th-MYCN CPM32 model therefore is a useful tool to dissect in detail mechanisms that drive metastasis and chemoresistance, and highlights dysregulation of signaling pathways such as JAK-STAT3 that could be targeted to improve treatment of aggressive disease. SIGNIFICANCE: An in vivo mouse model of high-risk treatment-resistant neuroblastoma exhibits changes in the tumor microenvironment, widespread metastases, and sensitivity to JAK1/2 inhibition.

Genomic and Transcriptomic Determinants of Therapy Resistance and Immune Landscape Evolution during Anti-EGFR Treatment in Colorectal Cancer.

Despite biomarker stratification, the anti-EGFR antibody cetuximab is only effective against a subgroup of colorectal cancers (CRCs). This genomic and transcriptomic analysis of the cetuximab resistance landscape in 35 RAS wild-type CRCs identified associations of NF1 and non-canonical RAS/RAF aberrations with primary resistance and validated transcriptomic CRC subtypes as non-genetic predictors of benefit. Sixty-four percent of biopsies with acquired resistance harbored no genetic resistance drivers. Most of these had switched from a cetuximab-sensitive transcriptomic subtype at baseline to a fibroblast- and growth factor-rich subtype at progression. Fibroblast-supernatant conferred cetuximab resistance in vitro, confirming a major role for non-genetic resistance through stromal remodeling. Cetuximab treatment increased cytotoxic immune infiltrates and PD-L1 and LAG3 immune checkpoint expression, potentially providing opportunities to treat cetuximab-resistant CRCs with immunotherapy.

Mutations in MAST1 Cause Mega-Corpus-Callosum Syndrome with Cerebellar Hypoplasia and Cortical Malformations.

Corpus callosum malformations are associated with a broad range of neurodevelopmental diseases. We report that de novo mutations in MAST1 cause mega-corpus-callosum syndrome with cerebellar hypoplasia and cortical malformations (MCC-CH-CM) in the absence of megalencephaly. We show that MAST1 is a microtubule-associated protein that is predominantly expressed in post-mitotic neurons and is present in both dendritic and axonal compartments. We further show that Mast1 null animals are phenotypically normal, whereas the deletion of a single amino acid (L278del) recapitulates the distinct neurological phenotype observed in patients. In animals harboring Mast1 microdeletions, we find that the PI3K/AKT3/mTOR pathway is unperturbed, whereas Mast2 and Mast3 levels are diminished, indicative of a dominant-negative mode of action. Finally, we report that de novo MAST1 substitutions are present in patients with autism and microcephaly, raising the prospect that mutations in this gene give rise to a spectrum of neurodevelopmental diseases.

Microenvironmental niche divergence shapes BRCA1-dysregulated ovarian cancer morphological plasticity.

How tumor microenvironmental forces shape plasticity of cancer cell morphology is poorly understood. Here, we conduct automated histology image and spatial statistical analyses in 514 high grade serous ovarian samples to define cancer morphological diversification within the spatial context of the microenvironment. Tumor spatial zones, where cancer cell nuclei diversify in shape, are mapped in each tumor. Integration of this spatially explicit analysis with omics and clinical data reveals a relationship between morphological diversification and the dysregulation of DNA repair, loss of nuclear integrity, and increased disease mortality. Within the Immunoreactive subtype, spatial analysis further reveals significantly lower lymphocytic infiltration within diversified zones compared with other tumor zones, suggesting that even immune-hot tumors contain cells capable of immune escape. Our findings support a model whereby a subpopulation of morphologically plastic cancer cells with dysregulated DNA repair promotes ovarian cancer progression through positive selection by immune evasion.

Ultra-Sensitive Mutation Detection and Genome-Wide DNA Copy Number Reconstruction by Error-Corrected Circulating Tumor DNA Sequencing.

BACKGROUND: Circulating free DNA sequencing (cfDNA-Seq) can portray cancer genome landscapes, but highly sensitive and specific technologies are necessary to accurately detect mutations with often low variant frequencies. METHODS: We developed a customizable hybrid-capture cfDNA-Seq technology using off-the-shelf molecular barcodes and a novel duplex DNA molecule identification tool for enhanced error correction. RESULTS: Modeling based on cfDNA yields from 58 patients showed that this technology, requiring 25 ng of cfDNA, could be applied to >95% of patients with metastatic colorectal cancer (mCRC). cfDNA-Seq of a 32-gene, 163.3-kbp target region detected 100% of single-nucleotide variants, with 0.15% variant frequency in spike-in experiments. Molecular barcode error correction reduced false-positive mutation calls by 97.5%. In 28 consecutively analyzed patients with mCRC, 80 out of 91 mutations previously detected by tumor tissue sequencing were called in the cfDNA. Call rates were similar for point mutations and indels. cfDNA-Seq identified typical mCRC driver mutations in patients in whom biopsy sequencing had failed or did not include key mCRC driver genes. Mutations only called in cfDNA but undetectable in matched biopsies included a subclonal resistance driver mutation to anti-EGFR antibodies in KRAS, parallel evolution of multiple PIK3CA mutations in 2 cases, and TP53 mutations originating from clonal hematopoiesis. Furthermore, cfDNA-Seq off-target read analysis allowed simultaneous genome-wide copy number profile reconstruction in 20 of 28 cases. Copy number profiles were validated by low-coverage whole-genome sequencing. CONCLUSIONS: This error-corrected, ultradeep cfDNA-Seq technology with a customizable target region and publicly available bioinformatics tools enables broad insights into cancer genomes and evolution. CLINICALTRIALSGOV IDENTIFIER: NCT02112357.

A Novel Strategy Combining Array-CGH, Whole-exome Sequencing and In Utero Electroporation in Rodents to Identify Causative Genes for Brain Malformations.

Birth defects that involve the cerebral cortex - also known as malformations of cortical development (MCD) - are important causes of intellectual disability and account for 20-40% of drug-resistant epilepsy in childhood. High-resolution brain imaging has facilitated in vivo identification of a large group of MCD phenotypes. Despite the advances in brain imaging, genomic analysis and generation of animal models, a straightforward workflow to systematically prioritize candidate genes and to test functional effects of putative mutations is missing. To overcome this problem, an experimental strategy enabling the identification of novel causative genes for MCD was developed and validated. This strategy is based on identifying candidate genomic regions or genes via array-CGH or whole-exome sequencing and characterizing the effects of their inactivation or of overexpression of specific mutations in developing rodent brains via in utero electroporation. This approach led to the identification of the C6orf70 gene, encoding for a putative vesicular protein, to the pathogenesis of periventricular nodular heterotopia, a MCD caused by defective neuronal migration.

Making the most of RNA-seq: Pre-processing sequencing data with Opossum for reliable SNP variant detection.

Identifying variants from RNA-seq (transcriptome sequencing) data is a cost-effective and versatile alternative to whole-genome sequencing. However, current variant callers do not generally behave well with RNA-seq data due to reads encompassing intronic regions. We have developed a software programme called Opossum to address this problem. Opossum pre-processes RNA-seq reads prior to variant calling, and although it has been designed to work specifically with Platypus, it can be used equally well with other variant callers such as GATK HaplotypeCaller. In this work, we show that using Opossum in conjunction with either Platypus or GATK HaplotypeCaller maintains precision and improves the sensitivity for SNP detection compared to the GATK Best Practices pipeline. In addition, using it in combination with Platypus offers a substantial reduction in run times compared to the GATK pipeline so it is ideal when there are only limited time or computational resources available.

Combination of Whole Genome Sequencing, Linkage, and Functional Studies Implicates a Missense Mutation in Titin as a Cause of Autosomal Dominant Cardiomyopathy With Features of Left Ventricular Noncompaction.

BACKGROUND: High throughput next-generation sequencing techniques have made whole genome sequencing accessible in clinical practice; however, the abundance of variation in the human genomes makes the identification of a disease-causing mutation on a background of benign rare variants challenging. METHODS AND RESULTS: Here we combine whole genome sequencing with linkage analysis in a 3-generation family affected by cardiomyopathy with features of autosomal dominant left ventricular noncompaction cardiomyopathy. A missense mutation in the giant protein titin is the only plausible disease-causing variant that segregates with disease among the 7 surviving affected individuals, with interrogation of the entire genome excluding other potential causes. This A178D missense mutation, affecting a conserved residue in the second immunoglobulin-like domain of titin, was introduced in a bacterially expressed recombinant protein fragment and biophysically characterized in comparison to its wild-type counterpart. Multiple experiments, including size exclusion chromatography, small-angle x ray scattering, and circular dichroism spectroscopy suggest partial unfolding and domain destabilization in the presence of the mutation. Moreover, binding experiments in mammalian cells show that the mutation markedly impairs binding to the titin ligand telethonin. CONCLUSIONS: Here we present genetic and functional evidence implicating the novel A178D missense mutation in titin as the cause of a highly penetrant familial cardiomyopathy with features of left ventricular noncompaction. This expands the spectrum of titin's roles in cardiomyopathies. It furthermore highlights that rare titin missense variants, currently often ignored or left uninterpreted, should be considered to be relevant for cardiomyopathies and can be identified by the approach presented here.

Autosomal dominant osteopetrosis associated with renal tubular acidosis is due to a CLCN7 mutation.

The aim of this study was to identify the causative mutation in a family with an unusual presentation of autosomal dominant osteopetrosis (OPT), proximal renal tubular acidosis (RTA), renal stones, epilepsy, and blindness, a combination of features not previously reported. We undertook exome sequencing of one affected and one unaffected family member, followed by targeted analysis of known candidate genes to identify the causative mutation. This identified a missense mutation (c.643G>A; p.Gly215Arg) in the gene encoding the chloride/proton antiporter 7 (gene CLCN7, protein CLC-7), which was confirmed by amplification refractory mutation system (ARMS)-PCR, and to be present in the three available patients. CLC-7 mutations are known to cause autosomal dominant OPT type 2, also called Albers-Schonberg disease, which is characterized by osteosclerosis, predominantly of the spine, pelvis and skull base, resulting in bone fragility and fractures. Albers-Schonberg disease is not reported to be associated with RTA, but autosomal recessive OPT type 3 (OPTB3) with RTA is associated with carbonic anhydrase type 2 (CA2) mutations. No mutations were detected in CA2 or any other genes known to cause proximal RTA. Neither CLCN7 nor CA2 mutations have previously been reported to be associated with renal stones or epilepsy. Thus, we identified a CLCN7 mutation in a family with autosomal dominant osteopetrosis, RTA, renal stones, epilepsy, and blindness. © 2016 Wiley Periodicals, Inc.

Premalignant SOX2 overexpression in the fallopian tubes of ovarian cancer patients: Discovery and validation studies.

Current screening methods for ovarian cancer can only detect advanced disease. Earlier detection has proved difficult because the molecular precursors involved in the natural history of the disease are unknown. To identify early driver mutations in ovarian cancer cells, we used dense whole genome sequencing of micrometastases and microscopic residual disease collected at three time points over three years from a single patient during treatment for high-grade serous ovarian cancer (HGSOC). The functional and clinical significance of the identified mutations was examined using a combination of population-based whole genome sequencing, targeted deep sequencing, multi-center analysis of protein expression, loss of function experiments in an in-vivo reporter assay and mammalian models, and gain of function experiments in primary cultured fallopian tube epithelial (FTE) cells. We identified frequent mutations involving a 40kb distal repressor region for the key stem cell differentiation gene SOX2. In the apparently normal FTE, the region was also mutated. This was associated with a profound increase in SOX2 expression (p<2(-16)), which was not found in patients without cancer (n=108). Importantly, we show that SOX2 overexpression in FTE is nearly ubiquitous in patients with HGSOCs (n=100), and common in BRCA1-BRCA2 mutation carriers (n=71) who underwent prophylactic salpingo-oophorectomy. We propose that the finding of SOX2 overexpression in FTE could be exploited to develop biomarkers for detecting disease at a premalignant stage, which would reduce mortality from this devastating disease.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders.

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.