Identification of liver cancer progenitors whose malignant progression depends on autocrine IL-6 signaling.

Hepatocellular carcinoma (HCC) is a slowly developing malignancy postulated to evolve from premalignant lesions in chronically damaged livers. However, it was never established that premalignant lesions actually contain tumor progenitors that give rise to cancer. Here, we describe isolation and characterization of HCC progenitor cells (HcPCs) from different mouse HCC models. Unlike fully malignant HCC, HcPCs give rise to cancer only when introduced into a liver undergoing chronic damage and compensatory proliferation. Although HcPCs exhibit a similar transcriptomic profile to bipotential hepatobiliary progenitors, the latter do not give rise to tumors. Cells resembling HcPCs reside within dysplastic lesions that appear several months before HCC nodules. Unlike early hepatocarcinogenesis, which depends on paracrine IL-6 production by inflammatory cells, due to upregulation of LIN28 expression, HcPCs had acquired autocrine IL-6 signaling that stimulates their in vivo growth and malignant progression. This may be a general mechanism that drives other IL-6-producing malignancies.

Chromosomal instability by mutations in the novel minor spliceosome component CENATAC.

Aneuploidy is the leading cause of miscarriage and congenital birth defects, and a hallmark of cancer. Despite this strong association with human disease, the genetic causes of aneuploidy remain largely unknown. Through exome sequencing of patients with constitutional mosaic aneuploidy, we identified biallelic truncating mutations in CENATAC (CCDC84). We show that CENATAC is a novel component of the minor (U12-dependent) spliceosome that promotes splicing of a specific, rare minor intron subtype. This subtype is characterized by AT-AN splice sites and relatively high basal levels of intron retention. CENATAC depletion or expression of disease mutants resulted in excessive retention of AT-AN minor introns in ˜ 100 genes enriched for nucleocytoplasmic transport and cell cycle regulators, and caused chromosome segregation errors. Our findings reveal selectivity in minor intron splicing and suggest a link between minor spliceosome defects and constitutional aneuploidy in humans.

Mutations in Epigenetic Regulation Genes Are a Major Cause of Overgrowth with Intellectual Disability.

To explore the genetic architecture of human overgrowth syndromes and human growth control, we performed experimental and bioinformatic analyses of 710 individuals with overgrowth (height and/or head circumference ≥+2 SD) and intellectual disability (OGID). We identified a causal mutation in 1 of 14 genes in 50% (353/710). This includes HIST1H1E, encoding histone H1.4, which has not been associated with a developmental disorder previously. The pathogenic HIST1H1E mutations are predicted to result in a product that is less effective in neutralizing negatively charged linker DNA because it has a reduced net charge, and in DNA binding and protein-protein interactions because key residues are truncated. Functional network analyses demonstrated that epigenetic regulation is a prominent biological process dysregulated in individuals with OGID. Mutations in six epigenetic regulation genes-NSD1, EZH2, DNMT3A, CHD8, HIST1H1E, and EED-accounted for 44% of individuals (311/710). There was significant overlap between the 14 genes involved in OGID and 611 genes in regions identified in GWASs to be associated with height (p = 6.84 × 10-8), suggesting that a common variation impacting function of genes involved in OGID influences height at a population level. Increased cellular growth is a hallmark of cancer and there was striking overlap between the genes involved in OGID and 260 somatically mutated cancer driver genes (p = 1.75 × 10-14). However, the mutation spectra of genes involved in OGID and cancer differ, suggesting complex genotype-phenotype relationships. These data reveal insights into the genetic control of human growth and demonstrate that exome sequencing in OGID has a high diagnostic yield.

BAP1 mutation is a frequent somatic event in peritoneal malignant mesothelioma.

BACKGROUND: Malignant mesothelioma (MM) arises from mesothelial cells that line the pleural, peritoneal and pericardial surfaces. The majority of MMs are pleural and have been associated with asbestos exposure. Previously, pleural MMs have been genetically characterized by the loss of BAP1 (40-60%) as well as loss of NF2 (75%) and CDKN2A (60%). The rare peritoneal form of MM occurs in ~10% cases. With only ~300 cases diagnosed in the US per year, its link to asbestos exposure is not clear and its mutational landscape unknown. METHODS: We analyzed the somatic mutational landscape of 12 peritoneal MM of epitheloid subtype using copy number analysis (N = 9), whole exome sequencing (N = 7) and targeted sequencing (N = 12). RESULTS: Peritoneal MM display few copy number alterations, with most samples having less than 10 Mbp total changes, mostly through deletions and no high copy number amplification. Chromosome band 3p21 encoding BAP1 is the most recurrently deleted region (5/9), while, in contrast to pleural MM, NF2 and CDKN2A are not affected. We further identified 87 non-silent mutations across 7 sequenced tumors, with a median of 8 mutated genes per tumor, resulting in a very low mutation rate (median 1.3 10(-6)). BAP1 was the only recurrently mutated gene (N = 3/7). In one additional case, loss of the entire chromosome 3 leaves a non-functional copy of BAP1 carrying a rare nonsense germline variant, thus suggesting a potential genetic predisposition in this patient. Finally, with targeted sequencing of BAP1 in 3 additional cases, we conclude that BAP1 is frequently altered through copy number losses (N = 3/12), mutations (N = 3/12) or both (N = 2/12) sometimes at a sub-clonal level. CONCLUSION: Our findings suggest a major role for BAP1 in peritoneal MM susceptibility and oncogenesis and indicate important molecular differences to pleural MM as well as potential strategies for therapy and prevention.

Mutascope: sensitive detection of somatic mutations from deep amplicon sequencing.

SUMMARY: We present Mutascope, a sequencing analysis pipeline specifically developed for the identification of somatic variants present at low-allelic fraction from high-throughput sequencing of amplicons from matched tumor-normal specimen. Using datasets reproducing tumor genetic heterogeneity, we demonstrate that Mutascope has a higher sensitivity and generates fewer false-positive calls than tools designed for shotgun sequencing or diploid genomes. AVAILABILITY: Freely available on the web at http://sourceforge.net/projects/mutascope/. CONTACT: oharismendy@ucsd.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of high-confidence somatic mutations in whole genome sequence of formalin-fixed breast cancer specimens.

The utilization of archived, formalin-fixed paraffin-embedded (FFPE) tumor samples for massive parallel sequencing has been challenging due to DNA damage and contamination with normal stroma. Here, we perform whole genome sequencing of DNA isolated from two triple-negative breast cancer tumors archived for >11 years as 5 µm FFPE sections and matched germline DNA. The tumor samples show differing amounts of FFPE damaged DNA sequencing reads revealed as relatively high alignment mismatch rates enriched for C · G > T · A substitutions compared to germline samples. This increase in mismatch rate is observable with as few as one million reads, allowing for an upfront evaluation of the sample integrity before whole genome sequencing. By applying innovative quality filters incorporating global nucleotide mismatch rates and local mismatch rates, we present a method to identify high-confidence somatic mutations even in the presence of FFPE induced DNA damage. This results in a breast cancer mutational profile consistent with previous studies and revealing potentially important functional mutations. Our study demonstrates the feasibility of performing genome-wide deep sequencing analysis of FFPE archived tumors of limited sample size such as residual cancer after treatment or metastatic biopsies.

Evaluation of ultra-deep targeted sequencing for personalized breast cancer care.

INTRODUCTION: The increasing number of targeted therapies, together with a deeper understanding of cancer genetics and drug response, have prompted major healthcare centers to implement personalized treatment approaches relying on high-throughput tumor DNA sequencing. However, the optimal way to implement this transformative methodology is not yet clear. Current assays may miss important clinical information such as the mutation allelic fraction, the presence of sub-clones or chromosomal rearrangements, or the distinction between inherited variants and somatic mutations. Here, we present the evaluation of ultra-deep targeted sequencing (UDT-Seq) to generate and interpret the molecular profile of 38 breast cancer patients from two academic medical centers. METHODS: We sequenced 47 genes in matched germline and tumor DNA samples from 38 breast cancer patients. The selected genes, or the pathways they belong to, can be targeted by drugs or are important in familial cancer risk or drug metabolism. RESULTS: Relying on the added value of sequencing matched tumor and germline DNA and using a dedicated analysis, UDT-Seq has a high sensitivity to identify mutations in tumors with low malignant cell content. Applying UDT-Seq to matched tumor and germline specimens from the 38 patients resulted in a proposal for at least one targeted therapy for 22 patients, the identification of tumor sub-clones in 3 patients, the suggestion of potential adverse drug effects in 3 patients and a recommendation for genetic counseling for 2 patients. CONCLUSION: Overall our study highlights the additional benefits of a sequencing strategy, which includes germline DNA and is optimized for heterogeneous tumor tissues.

Insights into BRCA Cancer Predisposition from Integrated Germline and Somatic Analyses in 7632 Cancers.

BACKGROUND: It is often assumed any cancer in a germline BRCA1 or BRCA2 (collectively termed BRCA) mutation carrier was caused by that mutation. It is also often assumed the occurrence of breast or ovarian cancer in an individual with a variant of uncertain significance (VUS) suggests the VUS is pathogenic. These assumptions have profound management implications for cancer patients and healthy individuals. METHODS: We compared the frequency of BRCA mutations, allele loss, and Signature 3 in 7632 individuals with 28 cancers and 1000 population controls. Because only increased frequency was the focus of the study, all statistical tests were one-sided. RESULTS: Individuals with breast or ovarian cancer had increased germline BRCA pathogenic mutation frequencies compared to controls (P = 1.0x10-10 and P = 1.4x10-34, respectively). There was no increase in other cancer types. Wild-type allele loss and Signature 3 were statistically significantly higher in breast and ovarian cancers with BRCA mutations compared with other cancers with BRCA mutations (P = 5.1x10-10 and P = 3.7x10-9) and cancers without BRCA mutations (P = 2.8x10-53 and P = 1.0x10-134). There was no difference between non-breast and non-ovarian cancers with BRCA mutations and cancers without BRCA mutations. Allele loss and Signature 3 were statistically significantly higher in breast and ovarian cancers in individuals with BRCA pathogenic mutations compared to those with VUS (P = 3.8x10-17 and P = 1.6x10-8) or benign variants (P = 1.2x10-28 and P = 2.2x10-10). There was no difference between individuals with BRCA VUS and those with benign variants. CONCLUSIONS: These data show that non-breast and non-ovarian cancers in individuals with germline BRCA pathogenic mutations are often not causally related to the mutation and that BRCA VUS are highly unlikely to be pathogenic. These results should reduce inappropriate management of germline BRCA information.

Evaluation of Cancer-Based Criteria for Use in Mainstream BRCA1 and BRCA2 Genetic Testing in Patients With Breast Cancer.

Importance: Increasing BRCA1 and BRCA2 (collectively termed herein as BRCA) gene testing is required to improve cancer management and prevent BRCA-related cancers. Objective: To evaluate mainstream genetic testing using cancer-based criteria in patients with cancer. Design, Setting, and Participants: A quality improvement study and cost-effectiveness analysis of different BRCA testing selection criteria and access procedures to evaluate feasibility, acceptability, and mutation detection performance was conducted at the Royal Marsden National Health Service Foundation Trust as part of the Mainstreaming Cancer Genetics (MCG) Programme. Participants included 1184 patients with cancer who were undergoing genetic testing between September 1, 2013, and February 28, 2017. Main Outcomes and Measures: Mutation rates, quality-adjusted life-years (QALYs), and incremental cost-effectiveness ratios were the primary outcomes. Results: Of the 1184 patients (1158 women [97.8%]) meeting simple cancer-based criteria, 117 had a BRCA mutation (9.9%). The mutation rate was similar in retrospective United Kingdom (10.2% [235 of 2294]) and prospective Malaysian (9.7% [103 of 1061]) breast cancer studies. If traditional family history criteria had been used, more than 50% of the mutation-positive individuals would have been missed. Of the 117 mutation-positive individuals, 115 people (98.3%) attended their genetics appointment and cascade to relatives is underway in all appropriate families (85 of 85). Combining with the equivalent ovarian cancer study provides 5 simple cancer-based criteria for BRCA testing with a 10% mutation rate: (1) ovarian cancer; (2) breast cancer diagnosed when patients are 45 years or younger; (3) 2 primary breast cancers, both diagnosed when patients are 60 years or younger; (4) triple-negative breast cancer; and (5) male breast cancer. A sixth criterion-breast cancer plus a parent, sibling, or child with any of the other criteria-can be added to address family history. Criteria 1 through 5 are considered the MCG criteria, and criteria 1 through 6 are considered the MCGplus criteria. Testing using MCG or MCGplus criteria is cost-effective with cost-effectiveness ratios of $1330 per discounted QALYs and $1225 per discounted QALYs, respectively, and appears to lead to cancer and mortality reductions (MCG: 804 cancers, 161 deaths; MCGplus: 1020 cancers, 204 deaths per year over 50 years). Use of MCG or MCGplus criteria might allow detection of all BRCA mutations in patients with breast cancer in the United Kingdom through testing one-third of patients. Feedback questionnaires from 259 patients and 23 cancer team members (12 oncologists, 8 surgeons, and 3 nurse specialists) showed acceptability of the process with 100% of patients pleased they had genetic testing and 100% of cancer team members confident to approve patients for genetic testing. Use of MCGplus criteria also appeared to be time and resource efficient, requiring 95% fewer genetic consultations than the traditional process. Conclusions and Relevance: This study suggests that mainstream testing using simple, cancer-based criteria might be able to efficiently deliver consistent, cost-effective, patient-centered BRCA testing.

Identification of new Wilms tumour predisposition genes: an exome sequencing study.

BACKGROUND: Wilms tumour is the most common childhood renal cancer and is genetically heterogeneous. While several Wilms tumour predisposition genes have been identified, there is strong evidence that further predisposition genes are likely to exist. Our study aim was to identify new predisposition genes for Wilms tumour. METHODS: In this exome sequencing study, we analysed lymphocyte DNA from 890 individuals with Wilms tumour, including 91 affected individuals from 49 familial Wilms tumour pedigrees. We used the protein-truncating variant prioritisation method to prioritise potential disease-associated genes for further assessment. We evaluated new predisposition genes in exome sequencing data that we generated in 334 individuals with 27 other childhood cancers and in exome data from The Cancer Genome Atlas obtained from 7632 individuals with 28 adult cancers. FINDINGS: We identified constitutional cancer-predisposing mutations in 33 individuals with childhood cancer. The three identified genes with the strongest signal in the protein-truncating variant prioritisation analyses were TRIM28, FBXW7, and NYNRIN. 21 of 33 individuals had a mutation in TRIM28; there was a strong parent-of-origin effect, with all ten inherited mutations being maternally transmitted (p=0·00098). We also found a strong association with the rare epithelial subtype of Wilms tumour, with 14 of 16 tumours being epithelial or epithelial predominant. There were no TRIM28 mutations in individuals with other childhood or adult cancers. We identified truncating FBXW7 mutations in four individuals with Wilms tumour and a de-novo non-synonymous FBXW7 mutation in a child with a rhabdoid tumour. Biallelic truncating mutations in NYNRIN were identified in three individuals with Wilms tumour, which is highly unlikely to have occurred by chance (p<0·0001). Finally, we identified two de-novo KDM3B mutations, supporting the role of KDM3B as a childhood cancer predisposition gene. INTERPRETATION: The four new Wilms tumour predisposition genes identified-TRIM28, FBXW7, NYNRIN, and KDM3B-are involved in diverse biological processes and, together with the other 17 known Wilms tumour predisposition genes, account for about 10% of Wilms tumour cases. The overlap between these 21 constitutionally mutated predisposition genes and 20 genes somatically mutated in Wilms tumour is limited, consisting of only four genes. We recommend that all individuals with Wilms tumour should be offered genetic testing and particularly, those with epithelial Wilms tumour should be offered TRIM28 genetic testing. Only a third of the familial Wilms tumour clusters we analysed were attributable to known genes, indicating that further Wilms tumour predisposition factors await discovery. FUNDING: Wellcome Trust.

Mutations in topoisomerase IIß result in a B cell immunodeficiency.

B cell development is a highly regulated process involving multiple differentiation steps, yet many details regarding this pathway remain unknown. Sequencing of patients with B cell-restricted immunodeficiency reveals autosomal dominant mutations in TOP2B. TOP2B encodes a type II topoisomerase, an essential gene required to alleviate topological stress during DNA replication and gene transcription, with no previously known role in B cell development. We use Saccharomyces cerevisiae, and knockin and knockout murine models, to demonstrate that patient mutations in TOP2B have a dominant negative effect on enzyme function, resulting in defective proliferation, survival of B-2 cells, causing a block in B cell development, and impair humoral function in response to immunization.

ICR142 Benchmarker: evaluating, optimising and benchmarking variant calling performance using the ICR142 NGS validation series.

Evaluating, optimising and benchmarking of next generation sequencing (NGS) variant calling performance are essential requirements for clinical, commercial and academic NGS pipelines. Such assessments should be performed in a consistent, transparent and reproducible fashion, using independently, orthogonally generated data. Here we present ICR142 Benchmarker, a tool to generate outputs for assessing germline base substitution and indel calling performance using the ICR142 NGS validation series, a dataset of Illumina platform-based exome sequence data from 142 samples together with Sanger sequence data at 704 sites. ICR142 Benchmarker provides summary and detailed information on the sensitivity, specificity and false detection rates of variant callers. ICR142 Benchmarker also automatically generates a single page report highlighting key performance metrics and how performance compares to widely-used open-source tools. We used ICR142 Benchmarker with VCF files outputted by GATK, OpEx and DeepVariant to create a benchmark for variant calling performance. This evaluation revealed pipeline-specific differences and shared challenges in variant calling, for example in detecting indels in short repeating sequence motifs. We next used ICR142 Benchmarker to perform regression testing with DeepVariant versions 0.5.2 and 0.6.1. This showed that v0.6.1 improves variant calling performance, but there was evidence of minor changes in indel calling behaviour that may benefit from attention. The data also allowed us to evaluate filters to optimise DeepVariant calling, and we recommend using 30 as the QUAL threshold for base substitution calls when using DeepVariant v0.6.1. Finally, we used ICR142 Benchmarker with VCF files from two commercial variant calling providers to facilitate optimisation of their in-house pipelines and to provide transparent benchmarking of their performance. ICR142 Benchmarker consistently and transparently analyses variant calling performance based on the ICR142 NGS validation series, using the standard VCF input and outputting informative metrics to enable user understanding of pipeline performance. ICR142 Benchmarker is freely available at https://github.com/RahmanTeamDevelopment/ICR142_Benchmarker/releases.

The Quality Sequencing Minimum (QSM): providing comprehensive, consistent, transparent next generation sequencing data quality assurance.

Next generation sequencing (NGS) is routinely used in clinical genetic testing. Quality management of NGS testing is essential to ensure performance is consistently and rigorously evaluated. Three primary metrics are used in NGS quality evaluation: depth of coverage, base quality and mapping quality. To provide consistency and transparency in the utilisation of these metrics we present the Quality Sequencing Minimum (QSM). The QSM defines the minimum quality requirement a laboratory has selected for depth of coverage (C), base quality (B) and mapping quality (M) and can be applied per base, exon, gene or other genomic region, as appropriate. The QSM format is CX_BY(P Y)_MZ(P Z). X is the parameter threshold for C, Y the parameter threshold for B, P Y the percentage of reads that must reach Y, Z the parameter threshold for M, P Z the percentage of reads that must reach Z. The data underlying the QSM is in the BAM file, so a QSM can be easily and automatically calculated in any NGS pipeline. We used the QSM to optimise cancer predisposition gene testing using the TruSight Cancer Panel (TSCP). We set the QSM as C50_B10(85)_M20(95). Test regions falling below the QSM were automatically flagged for review, with 100/1471 test regions QSM-flagged in multiple individuals. Supplementing these regions with 132 additional probes improved performance in 85/100. We also used the QSM to optimise testing of genes with pseudogenes such as PTEN and PMS2. In TSCP data from 960 individuals the median number of regions that passed QSM per sample was 1429 (97%).  Importantly, the QSM can be used at an individual report level to provide succinct, comprehensive quality assurance information about individual test performance. We believe many laboratories would find the QSM useful. Furthermore, widespread adoption of the QSM would facilitate consistent, transparent reporting of genetic test performance by different laboratories.

CoverView: a sequence quality evaluation tool for next generation sequencing data.

Quality assurance and quality control are essential for robust next generation sequencing (NGS). Here we present CoverView, a fast, flexible, user-friendly quality evaluation tool for NGS data. CoverView processes mapped sequencing reads and user-specified regions to report depth of coverage, base and mapping quality metrics with increasing levels of detail from a chromosome-level summary to per-base profiles. CoverView can flag regions that do not fulfil user-specified quality requirements, allowing suboptimal data to be systematically and automatically presented for review. It also provides an interactive graphical user interface (GUI) that can be opened in a web browser and allows intuitive exploration of results. We have integrated CoverView into our accredited clinical cancer predisposition gene testing laboratory that uses the TruSight Cancer Panel (TSCP). CoverView has been invaluable for optimisation and quality control of our testing pipeline, providing transparent, consistent quality metric information and automatic flagging of regions that fall below quality thresholds. We demonstrate this utility with TSCP data from the Genome in a Bottle reference sample, which CoverView analysed in 13 seconds. CoverView uses data routinely generated by NGS pipelines, reads standard input formats, and rapidly creates easy-to-parse output text (.txt) files that are customised by a simple configuration file. CoverView can therefore be easily integrated into any NGS pipeline. CoverView and detailed documentation for its use are freely available at github.com/RahmanTeamDevelopment/CoverView/releases and www.icr.ac.uk/CoverView.

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

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

Biallelic TRIP13 mutations predispose to Wilms tumor and chromosome missegregation.

Through exome sequencing, we identified six individuals with biallelic loss-of-function mutations in TRIP13. All six developed Wilms tumor. Constitutional mosaic aneuploidies, microcephaly, developmental delay and seizures, which are features of mosaic variegated aneuploidy (MVA) syndrome, were more variably present. Through functional studies, we show that TRIP13-mutant patient cells have no detectable TRIP13 and have substantial impairment of the spindle assembly checkpoint (SAC), leading to a high rate of chromosome missegregation. Accurate segregation, as well as SAC proficiency, is rescued by restoring TRIP13 function. Individuals with biallelic TRIP13 or BUB1B mutations have a high risk of embryonal tumors, and here we show that their cells display severe SAC impairment. MVA due to biallelic CEP57 mutations, or of unknown cause, is not associated with embryonal tumors and cells from these individuals show minimal SAC deficiency. These data provide insights into the complex relationships between aneuploidy and carcinogenesis.

Mutational Profiling Can Establish Clonal or Independent Origin in Synchronous Bilateral Breast and Other Tumors.

BACKGROUND: Synchronous tumors can be independent primary tumors or a primary-metastatic (clonal) pair, which may have clinical implications. Mutational profiling of tumor DNA is increasingly common in the clinic. We investigated whether mutational profiling can distinguish independent from clonal tumors in breast and other cancers, using a carefully defined test based on the Clonal Likelihood Score (CLS = 100 x # shared high confidence (HC) mutations/ # total HC mutations). METHODS: Statistical properties of a formal test using the CLS were investigated. A high CLS is evidence in favor of clonality; the test is implemented as a one-sided binomial test of proportions. Test parameters were empirically determined using 16,422 independent breast tumor pairs and 15 primary-metastatic tumor pairs from 10 cancer types using The Cancer Genome Atlas. RESULTS: We validated performance of the test with its established parameters, using five published data sets comprising 15,758 known independent tumor pairs (maximum CLS = 4.1%, minimum p-value = 0.48) and 283 known tumor clonal pairs (minimum CLS 13%, maximum p-value <0.01), across renal cell, testicular, and colorectal cancer. The CLS test correctly classified all validation samples but one, which it appears may have been incorrectly classified in the published data. As proof-of-concept we then applied the CLS test to two new cases of invasive synchronous bilateral breast cancer at our institution, each with one hormone receptor positive (ER+/PR+/HER2-) lobular and one triple negative ductal carcinoma. High confidence mutations were identified by exome sequencing and results were validated using deep targeted sequencing. The first tumor pair had CLS of 81% (p-value < 10-15), supporting clonality. In the second pair, no common mutations of 184 variants were validated (p-value >0.99), supporting independence. A plausible molecular mechanism for the shift from hormone receptor positive to triple negative was identified in the clonal pair. CONCLUSION: We have developed the statistical properties of a carefully defined Clonal Likelihood Score test from mutational profiling of tumor DNA. Under identified conditions, the test appears to reliably distinguish between synchronous tumors of clonal and of independent origin in several cancer types. This approach may have scientific and clinical utility.

Whole-exome sequencing reveals the mutational spectrum of testicular germ cell tumours.

Testicular germ cell tumours (TGCTs) are the most common cancer in young men. Here we perform whole-exome sequencing (WES) of 42 TGCTs to comprehensively study the cancer's mutational profile. The mutation rate is uniformly low in all of the tumours (mean 0.5 mutations per Mb) as compared with common cancers, consistent with the embryological origin of TGCT. In addition to expected copy number gain of chromosome 12p and mutation of KIT, we identify recurrent mutations in the tumour suppressor gene CDC27 (11.9%). Copy number analysis reveals recurring amplification of the spermatocyte development gene FSIP2 (15.3%) and a 0.4 Mb region at Xq28 (15.3%). Two treatment-refractory patients are shown to harbour XRCC2 mutations, a gene strongly implicated in defining cisplatin resistance. Our findings provide further insights into genes involved in the development and progression of TGCT.

Mutations in the transcriptional repressor REST predispose to Wilms tumor.

Wilms tumor is the most common childhood renal cancer. To identify mutations that predispose to Wilms tumor, we are conducting exome sequencing studies. Here we describe 11 different inactivating mutations in the REST gene (encoding RE1-silencing transcription factor) in four familial Wilms tumor pedigrees and nine non-familial cases. Notably, no similar mutations were identified in the ICR1000 control series (13/558 versus 0/993; P < 0.0001) or in the ExAC series (13/558 versus 0/61,312; P < 0.0001). We identified a second mutational event in two tumors, suggesting that REST may act as a tumor-suppressor gene in Wilms tumor pathogenesis. REST is a zinc-finger transcription factor that functions in cellular differentiation and embryonic development. Notably, ten of 11 mutations clustered within the portion of REST encoding the DNA-binding domain, and functional analyses showed that these mutations compromise REST transcriptional repression. These data establish REST as a Wilms tumor predisposition gene accounting for ∼2% of Wilms tumor.

Genome-wide mutational landscape of mucinous carcinomatosis peritonei of appendiceal origin.

BACKGROUND: Mucinous neoplasms of the appendix (MNA) are rare tumors which may progress from benign to malignant disease with an aggressive biological behavior. MNA is often diagnosed after metastasis to the peritoneal surfaces resulting in mucinous carcinomatosis peritonei (MCP). Genetic alterations in MNA are poorly characterized due to its low incidence, the hypo-cellularity of MCPs, and a lack of relevant pre-clinical models. As such, application of targeted therapies to this disease is limited to those developed for colorectal cancer and not based on molecular rationale. METHODS: We sequenced the whole exomes of 10 MCPs of appendiceal origin to identify genome-wide somatic mutations and copy number aberrations and validated significant findings in 19 additional cases. RESULTS: Our study demonstrates that MNA has a different molecular makeup than colorectal cancer. Most tumors have co-existing oncogenic mutations in KRAS (26/29) and GNAS (20/29) and are characterized by downstream PKA activation. High-grade tumors are GNAS wild-type (5/6), suggesting they do not progress from low-grade tumors. MNAs do share some genetic alterations with colorectal cancer including gain of 1q (5/10), Wnt, and TGFß pathway alterations. In contrast, mutations in TP53 (1/10) and APC (0/10), common in colorectal cancer, are rare in MNA. Concurrent activation of the KRAS and GNAS mediated signaling pathways appears to be shared with pancreatic intraductal papillary mucinous neoplasm. CONCLUSIONS: MNA genome-wide mutational analysis reveals genetic alterations distinct from colorectal cancer, in support of its unique pathophysiology and suggests new targeted therapeutic opportunities.