Discrepancies in tumor mutation burden reporting from sequential endobronchial ultrasound transbronchial needle aspiration samples within single lymph node stations - brief report.

Introduction: Tumour Mutation Burden (TMB) is a potential biomarker for immune cancer therapies. Here we investigated parameters that might affect TMB using duplicate cytology smears obtained from endobronchial ultrasound transbronchial needle aspiration (EBUS TBNA)-sampled malignant lymph nodes. Methods: Individual Diff-Quik cytology smears were prepared for each needle pass. DNA extracted from each smear underwent sequencing using large gene panel (TruSight Oncology 500 (TSO500 - Illumina)). TMB was estimated using the TSO500 Local App v. 2.0 (Illumina). Results: Twenty patients had two or more Diff-Quik smears (total 45 smears) which passed sequencing quality control. Average smear TMB was 8.7 ± 5.0 mutations per megabase (Mb). Sixteen of the 20 patients had paired samples with minimal differences in TMB score (average difference 1.3 ± 0.85). Paired samples from 13 patients had concordant TMB (scores below or above a threshold of 10 mutations/Mb). Markedly discrepant TMB was observed in four cases, with an average difference of 11.3 ± 2.7 mutations/Mb. Factors affecting TMB calling included sample tumour content, the amount of DNA used in sequencing, and bone fide heterogeneity of node tumour between paired samples. Conclusion: TMB assessment is feasible from EBUS-TBNA smears from a single needle pass. Repeated samples of a lymph node station have minimal variation in TMB in most cases. However, this novel data shows how tumour content and minor change in site of node sampling can impact TMB. Further study is needed on whether all node aspirates should be combined in 1 sample, or whether testing independent nodes using smears is needed.

Evaluating Diff-Quik cytology smears for large-panel mutation testing in lung cancer-Predicting DNA content and success with low-malignant-cellularity samples.

BACKGROUND: Cytology smears are commonly collected during endobronchial ultrasound-guided transbronchial needle aspiration (EBUS TBNA) procedures but are rarely used for molecular testing. Studies are needed to demonstrate their great potential, in particular for the prediction of malignant cell DNA content and for utility in molecular diagnostics using large gene panels. METHODS: A prospective study was performed on samples from 66 patients with malignant lymph nodes who underwent EBUS TBNA. All patients had air-dried, Diff-Quik cytology smears and formalin-fixed, paraffin-embedded cell blocks collected for cytopathology and molecular testing. One hundred eighty-five smears were evaluated by microscopy to estimate malignant cell percentage and abundance and to calculate smear size and were subjected to DNA extraction. DNA from 56 smears from 27 patients was sequenced with the TruSight Oncology 500 assay (Illumina). RESULTS: Each microscopy parameter had a significant effect on the DNA yield. An algorithm was developed that predicted a >50-ng DNA yield of a smear with an area under the curve of 0.86. Fifty DNA samples (89%) with varying malignant yields were successfully sequenced. Low-malignant-cell content (<25%) and smear area (<15%) were the main reasons for failure. All standard-of-care mutations were detected in replicate smears from individual patients, regardless of malignant cell content. Tier 1/2 mutations were discovered in two cases where standard-of-care specimens were inadequate for sequencing. Smears were scored for tumor mutation burden. CONCLUSIONS: Microscopy of Diff-Quik smears can triage samples for comprehensive panel sequencing, which highlights smears as an excellent alternative to traditional testing with cell blocks.

Comparative Genomics Provides Etiologic and Biological Insight into Melanoma Subtypes.

Melanoma is a cancer of melanocytes, with multiple subtypes based on body site location. Cutaneous melanoma is associated with skin exposed to ultraviolet radiation; uveal melanoma occurs in the eyes; mucosal melanoma occurs in internal mucous membranes; and acral melanoma occurs on the palms, soles, and nail beds. Here, we present the largest whole-genome sequencing study of melanoma to date, with 570 tumors profiled, as well as methylation and RNA sequencing for subsets of tumors. Uveal melanoma is genomically distinct from other melanoma subtypes, harboring the lowest tumor mutation burden and with significantly mutated genes in the G-protein signaling pathway. Most cutaneous, acral, and mucosal melanomas share alterations in components of the MAPK, PI3K, p53, p16, and telomere pathways. However, the mechanism by which these pathways are activated or inactivated varies between melanoma subtypes. Additionally, we identify potential novel germline predisposition genes for some of the less common melanoma subtypes. SIGNIFICANCE: This is the largest whole-genome analysis of melanoma to date, comprehensively comparing the genomics of the four major melanoma subtypes. This study highlights both similarities and differences between the subtypes, providing insights into the etiology and biology of melanoma. This article is highlighted in the In This Issue feature, p. 2711.

Comprehensive genomic and tumour immune profiling reveals potential therapeutic targets in malignant pleural mesothelioma.

BACKGROUND: Malignant pleural mesothelioma (MPM) has a poor overall survival with few treatment options. Whole genome sequencing (WGS) combined with the immune features of MPM offers the prospect of identifying changes that could inform future clinical trials. METHODS: We analysed somatic mutations from 229 MPM samples, including previously published data and 58 samples that had undergone WGS within this study. This was combined with RNA-seq analysis to characterize the tumour immune environment. RESULTS: The comprehensive genome analysis identified 12 driver genes, including new candidate genes. Whole genome doubling was a frequent event that correlated with shorter survival. Mutational signature analysis revealed SBS5/40 were dominant in 93% of samples, and defects in homologous recombination repair were infrequent in our cohort. The tumour immune environment contained high M2 macrophage infiltrate linked with MMP2, MMP14, TGFB1 and CCL2 expression, representing an immune suppressive environment. The expression of TGFB1 was associated with overall survival. A small subset of samples (less than 10%) had a higher proportion of CD8 T cells and a high cytolytic score, suggesting a 'hot' immune environment independent of the somatic mutations. CONCLUSIONS: We propose accounting for genomic and immune microenvironment status may influence therapeutic planning in the future.

Combined Inhibition of G9a and EZH2 Suppresses Tumor Growth via Synergistic Induction of IL24-Mediated Apoptosis.

G9a and EZH2 are two histone methyltransferases commonly upregulated in several cancer types, yet the precise roles that these enzymes play cooperatively in cancer is unclear. We demonstrate here that frequent concurrent upregulation of both G9a and EZH2 occurs in several human tumors. These methyltransferases cooperatively repressed molecular pathways responsible for tumor cell death. In genetically distinct tumor subtypes, concomitant inhibition of G9a and EZH2 potently induced tumor cell death, highlighting the existence of tumor cell survival dependency at the epigenetic level. G9a and EZH2 synergistically repressed expression of genes involved in the induction of endoplasmic reticulum (ER) stress and the production of reactive oxygen species. IL24 was essential for the induction of tumor cell death and was identified as a common target of G9a and EZH2. Loss of function of G9a and EZH2 activated the IL24-ER stress axis and increased apoptosis in cancer cells while not affecting normal cells. These results indicate that G9a and EZH2 promotes the evasion of ER stress-mediated apoptosis by repressing IL24 transcription, therefore suggesting that their inhibition may represent a potential therapeutic strategy for solid cancers. SIGNIFICANCE: These findings demonstrate a novel role for G9a and EZH2 histone methyltransferases in suppressing apoptosis, which can be targeted with small molecule inhibitors as a potential approach to improve solid cancer treatment.

Patient-derived xenograft models capture genomic heterogeneity in endometrial cancer.

BACKGROUND: Endometrial cancer (EC) is a major gynecological cancer with increasing incidence. It comprises four molecular subtypes with differing etiology, prognoses, and responses to chemotherapy. In the future, clinical trials testing new single agents or combination therapies will be targeted to the molecular subtype most likely to respond. As pre-clinical models that faithfully represent the molecular subtypes of EC are urgently needed, we sought to develop and characterize a panel of novel EC patient-derived xenograft (PDX) models. METHODS: Here, we report whole exome or whole genome sequencing of 11 PDX models and their matched primary tumor. Analysis of multiple PDX lineages and passages was performed to study tumor heterogeneity across lineages and/or passages. Based on recent reports of frequent defects in the homologous recombination (HR) pathway in EC, we assessed mutational signatures and HR deficiency scores and correlated these with in vivo responses to the PARP inhibitor (PARPi) talazoparib in six PDXs representing the copy number high/p53-mutant and mismatch-repair deficient molecular subtypes of EC. RESULTS: PDX models were successfully generated from grade 2/3 tumors, including three uterine carcinosarcomas. The models showed similar histomorphology to the primary tumors and represented all four molecular subtypes of EC, including five mismatch-repair deficient models. The different PDX lineages showed a wide range of inter-tumor and intra-tumor heterogeneity. However, for most PDX models, one arm recapitulated the molecular landscape of the primary tumor without major genomic drift. An in vivo response to talazoparib was detected in four copy number high models. Two models (carcinosarcomas) showed a response consistent with stable disease and two models (one copy number high serous EC and another carcinosarcoma) showed significant tumor growth inhibition, albeit one consistent with progressive disease; however, all lacked the HR deficiency genomic signature. CONCLUSIONS: EC PDX models represent the four molecular subtypes of disease and can capture intra-tumor heterogeneity of the original primary tumor. PDXs of the copy number high molecular subtype showed sensitivity to PARPi; however, deeper and more durable responses will likely require combination of PARPi with other agents.

Multiomic profiling of checkpoint inhibitor-treated melanoma: Identifying predictors of response and resistance, and markers of biological discordance.

We concurrently examine the whole genome, transcriptome, methylome, and immune cell infiltrates in baseline tumors from 77 patients with advanced cutaneous melanoma treated with anti-PD-1 with or without anti-CTLA-4. We show that high tumor mutation burden (TMB), neoantigen load, expression of IFNγ-related genes, programmed death ligand expression, low PSMB8 methylation (therefore high expression), and T cells in the tumor microenvironment are associated with response to immunotherapy. No specific mutation correlates with therapy response. A multivariable model combining the TMB and IFNγ-related gene expression robustly predicts response (89% sensitivity, 53% specificity, area under the curve [AUC], 0.84); tumors with high TMB and a high IFNγ signature show the best response to immunotherapy. This model validates in an independent cohort (80% sensitivity, 59% specificity, AUC, 0.79). Except for a JAK3 loss-of-function mutation, for patients who did not respond as predicted there is no obvious biological mechanism that clearly explained their outlier status, consistent with intratumor and intertumor heterogeneity in response to immunotherapy.

IFN-λ therapy prevents severe gastrointestinal graft-versus-host disease.

Immunopathology and intestinal stem cell (ISC) loss in the gastrointestinal (GI) tract is the prima facie manifestation of graft-versus-host disease (GVHD) and is responsible for significant mortality after allogeneic bone marrow transplantation (BMT). Approaches to prevent GVHD to date focus on immune suppression. Here, we identify interferon-λ (IFN-λ; interleukin-28 [IL-28]/IL-29) as a key protector of GI GVHD immunopathology, notably within the ISC compartment. Ifnlr1-/- mice displayed exaggerated GI GVHD and mortality independent of Paneth cells and alterations to the microbiome. Ifnlr1-/- intestinal organoid growth was significantly impaired, and targeted Ifnlr1 deficiency exhibited effects intrinsic to recipient Lgr5+ ISCs and natural killer cells. PEGylated recombinant IL-29 (PEG-rIL-29) treatment of naive mice enhanced Lgr5+ ISC numbers and organoid growth independent of both IL-22 and type I IFN and modulated proliferative and apoptosis gene sets in Lgr5+ ISCs. PEG-rIL-29 treatment improved survival, reduced GVHD severity, and enhanced epithelial proliferation and ISC-derived organoid growth after BMT. The preservation of ISC numbers in response to PEG-rIL-29 after BMT occurred both in the presence and absence of IFN-λ-signaling in recipient natural killer cells. IFN-λ is therefore an attractive and rapidly testable approach to prevent ISC loss and immunopathology during GVHD.

Injection site vaccinology of a recombinant vaccinia-based vector reveals diverse innate immune signatures.

Poxvirus systems have been extensively used as vaccine vectors. Herein a RNA-Seq analysis of intramuscular injection sites provided detailed insights into host innate immune responses, as well as expression of vector and recombinant immunogen genes, after vaccination with a new multiplication defective, vaccinia-based vector, Sementis Copenhagen Vector. Chikungunya and Zika virus immunogen mRNA and protein expression was associated with necrosing skeletal muscle cells surrounded by mixed cellular infiltrates. The multiple adjuvant signatures at 12 hours post-vaccination were dominated by TLR3, 4 and 9, STING, MAVS, PKR and the inflammasome. Th1 cytokine signatures were dominated by IFNγ, TNF and IL1ß, and chemokine signatures by CCL5 and CXCL12. Multiple signatures associated with dendritic cell stimulation were evident. By day seven, vaccine transcripts were absent, and cell death, neutrophil, macrophage and inflammation annotations had abated. No compelling arthritis signatures were identified. Such injection site vaccinology approaches should inform refinements in poxvirus-based vector design.

Whole-genome sequencing of acral melanoma reveals genomic complexity and diversity.

To increase understanding of the genomic landscape of acral melanoma, a rare form of melanoma occurring on palms, soles or nail beds, whole genome sequencing of 87 tumors with matching transcriptome sequencing for 63 tumors was performed. Here we report that mutational signature analysis reveals a subset of tumors, mostly subungual, with an ultraviolet radiation signature. Significantly mutated genes are BRAF, NRAS, NF1, NOTCH2, PTEN and TYRP1. Mutations and amplification of KIT are also common. Structural rearrangement and copy number signatures show that whole genome duplication, aneuploidy and complex rearrangements are common. Complex rearrangements occur recurrently and are associated with amplification of TERT, CDK4, MDM2, CCND1, PAK1 and GAB2, indicating potential therapeutic options.

Using whole-genome sequencing data to derive the homologous recombination deficiency scores.

The homologous recombination deficiency (HRD) score was developed using whole-genome copy number data derived from arrays as a way to infer deficiency in the homologous recombination DNA damage repair pathway (in particular BRCA1 or BRCA2 deficiency) in breast cancer samples. The score has utility in understanding tumour biology and may be indicative of response to certain therapeutic strategies. Studies have used whole-exome sequencing to derive the HRD score, however, with increasing use of whole-genome sequencing (WGS) to characterise tumour genomes, there has yet to be a comprehensive comparison between HRD scores derived by array versus WGS. Here we demonstrate that there is both a high correlation and a good agreement between array- and WGS-derived HRD scores and between the scores derived from WGS and downsampled WGS to represent shallow WGS. For samples with an HRD score close to threshold for stratifying HR proficiency or deficiency there was however some disagreement in the HR status between array and WGS data, highlighting the importance of not relying on a single method of ascertaining the homologous recombination status of a tumour.

APC Mutation Marks an Aggressive Subtype of BRAF Mutant Colorectal Cancers.

BACKGROUND: WNT activation is a hallmark of colorectal cancer. BRAF mutation is present in 15% of colorectal cancers, and the role of mutations in WNT signaling regulators in this context is unclear. Here, we evaluate the mutational landscape of WNT signaling regulators in BRAF mutant cancers. METHODS: we performed exome-sequencing on 24 BRAF mutant colorectal cancers and analyzed these data in combination with 175 publicly available BRAF mutant colorectal cancer exomes. We assessed the somatic mutational landscape of WNT signaling regulators, and performed hotspot and driver mutation analyses to identify potential drivers of WNT signaling. The effects of Apc and Braf mutation were modelled, in vivo, using the Apcmin/+ and BrafV637/Villin-CreERT2/+ mouse, respectively. RESULTS: RNF43 was the most frequently mutated WNT signaling regulator (41%). Mutations in the beta-catenin destruction complex occurred in 48% of cancers. Hotspot analyses identified potential cancer driver genes in the WNT signaling cascade, including MEN1, GNG12 and WNT16. Truncating APC mutation was identified in 20.8% of cancers. Truncating APC mutation was associated with early age at diagnosis (p < 2 × 10-5), advanced stage (p < 0.01), and poor survival (p = 0.026). Apcmin/+/BrafV637 animals had more numerous and larger SI and colonic lesions (p < 0.0001 and p < 0.05, respectively), and a markedly reduced survival (median survival: 3.2 months, p = 8.8 × 10-21), compared to animals with Apc or Braf mutation alone. CONCLUSIONS: the WNT signaling axis is frequently mutated in BRAF mutant colorectal cancers. WNT16 and MEN1 may be novel drivers of aberrant WNT signaling in colorectal cancer. Co-mutation of BRAF and APC generates an extremely aggressive neoplastic phenotype that is associated with poor patient outcome.

Whole-genome landscape of mucosal melanoma reveals diverse drivers and therapeutic targets.

Knowledge of key drivers and therapeutic targets in mucosal melanoma is limited due to the paucity of comprehensive mutation data on this rare tumor type. To better understand the genomic landscape of mucosal melanoma, here we describe whole genome sequencing analysis of 67 tumors and validation of driver gene mutations by exome sequencing of 45 tumors. Tumors have a low point mutation burden and high numbers of structural variants, including recurrent structural rearrangements targeting TERT, CDK4 and MDM2. Significantly mutated genes are NRAS, BRAF, NF1, KIT, SF3B1, TP53, SPRED1, ATRX, HLA-A and CHD8. SF3B1 mutations occur more commonly in female genital and anorectal melanomas and CTNNB1 mutations implicate a role for WNT signaling defects in the genesis of some mucosal melanomas. TERT aberrations and ATRX mutations are associated with alterations in telomere length. Mutation profiles of the majority of mucosal melanomas suggest potential susceptibility to CDK4/6 and/or MEK inhibitors.

Molecular Genomic Profiling of Melanocytic Nevi.

The benign melanocytic nevus is the most common tumor in humans and rarely transforms into cutaneous melanoma. Elucidation of the nevus genome is required to better understand the molecular steps of progression to melanoma. We performed whole genome sequencing on a series of 14 benign melanocytic nevi consisting of both congenital and acquired types. All nevi had driver mutations in the MAPK signaling pathway, either BRAF V600E or NRAS Q61R/L. No additional definite driver mutations were identified. Somatic mutations in nevi with higher mutation loads showed a predominance of mutational signatures 7a and 7b, consistent with UVR exposure, whereas nevi with lower mutation loads (including all three congenital nevi) had a predominance of the ubiquitous signatures 1 and 5. Two nevi had mutations in promoter regions predicted to bind E26 transformation-specific family transcription factors, as well as subclonal mutations in the TERT promoter. This paper presents whole genome data from melanocytic nevi. We confirm that UVR is involved in the etiology of a subset of nevi. This study also establishes that TERT promoter mutations are present in morphologically benign skin nevi in subclonal populations, which has implications regarding the interpretation of this emerging biomarker in sensitive assays.

Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma.

Glioblastomas are the most common and lethal neoplasms of the central nervous system. Neighbouring glioma cells maintain extreme degrees of genetic and phenotypic variation that form intratumoural heterogeneity. This genetic diversity allows the most adaptive tumour clones to develop treatment resistance, ultimately leading to disease recurrence. We aimed to model this phenomenon and test the effectiveness of several targeted therapeutic interventions to overcome therapy resistance. Heterogeneous tumour masses were first deconstructed into single tumour cells, which were expanded independently as single-cell clones. Single nucleotide polymorphism arrays, whole-genome and RNA sequencing, and CpG methylation analysis validated the unique molecular profile of each tumour clone, which displayed distinct pathologic features, including cell morphology, growth rate, and resistance to temozolomide and ionizing radiation. We also identified variable sensitivities to AURK, CDK, and EGFR inhibitors which were consistent with the heterogeneous molecular alterations that each clone harboured. These targeted therapies effectively eliminated the temozolomide- and/or irradiation-resistant clones and also parental polyclonal cells. Our findings indicate that polyclonal tumours create a dynamic environment that consists of diverse tumour elements and treatment responses. Designing targeted therapies based on a range of molecular profiles can be a more effective strategy to eradicate treatment resistance, recurrence, and metastasis.

Complex structural rearrangements are present in high-grade dysplastic Barrett's oesophagus samples.

BACKGROUND: Oesophageal adenocarcinoma (EAC) incidence is increasing and has a poor survival rate. Barrett's oesophagus (BE) is a precursor condition that is associated with EAC and often occurs in conjunction with chronic gastro-oesophageal reflux, however many individuals diagnosed with BE never progress to cancer. An understanding of the genomic features of BE and EAC may help with the early identification of at-risk individuals. METHODS: In this study, we assessed the genomic features of 16 BE samples using whole-genome sequencing. These included non-dysplastic samples collected at two time-points from two BE patients who had not progressed to EAC over several years. Seven other non-dysplastic samples and five dysplastic BE samples with high-grade dysplasia were also examined. We compared the genome profiles of these 16 BE samples with 22 EAC samples. RESULTS: We observed that samples from the two non-progressor individuals had low numbers of somatic single nucleotide variants, indels and structural variation events compared to dysplastic and the remaining non-dysplastic BE. EAC had the highest level of somatic genomic variations. Mutational signature 17, which is common in EAC, was also present in non-dysplastic and dysplastic BE, but was not present in the non-progressors. Many dysplastic samples had mutations in genes previously reported in EAC, whereas only mutations in CDKN2A or in the fragile site genes appeared common in non-dysplastic samples. Rearrangement signatures were used to identify a signature associated with localised complex events such as chromothripsis and breakage fusion-bridge that are characteristic of EACs. Two dysplastic BE samples had a high contribution of this signature and contained evidence of localised rearrangements. Two other dysplastic samples also had regions of localised structural rearrangements. There was no evidence for complex events in non-dysplastic samples. CONCLUSIONS: The presence of complex localised rearrangements in dysplastic samples indicates a need for further investigations into the role such events play in the progression from BE to EAC.

A plugin for the Ensembl Variant Effect Predictor that uses MaxEntScan to predict variant spliceogenicity.

SUMMARY: Assessing the pathogenicity of genetic variants can be a complex and challenging task. Spliceogenic variants, which alter mRNA splicing, may yield mature transcripts that encode non-functional protein products, an important predictor of Mendelian disease risk. However, most variant annotation tools do not adequately assess spliceogenicity outside the native splice site and thus the disease-causing potential of variants in other intronic and exonic regions is often overlooked. Here, we present a plugin for the Ensembl Variant Effect Predictor that packages MaxEntScan and extends its functionality to provide splice site predictions using a maximum entropy model. The plugin incorporates a sliding window algorithm to predict splice site loss or gain for any variant that overlaps a transcript feature. We also demonstrate the utility of the plugin by comparing our predictions to two mRNA splicing datasets containing several cancer-susceptibility genes. AVAILABILITY AND IMPLEMENTATION: Source code is freely available under the Apache License, Version 2.0: https://github.com/Ensembl/VEP_plugins. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization of a novel breast cancer cell line derived from a metastatic bone lesion of a breast cancer patient.

PURPOSE: We aimed to generate and characterize a novel cell line from a breast cancer bone metastasis to better study the progression of the disease. METHODS: The cell line, P7731, was derived from a metastatic bone lesion of a breast cancer patient and assessed for marker expression. P7731 was analyzed for DNA copy number variation, somatic mutations, and gene expression and was compared with the primary tumor. RESULTS: P7731 cells are negative for estrogen receptor alpha (ERα), progesterone receptor (PR), and HER2 (triple-negative); strongly express vimentin (100% of cells positive) and also express cytokeratins 8/18 and 19 but at lower frequencies. Flow cytometry indicates P7731 cells are predominantly CD44+/CD49f+/EpCAM-, consistent with a primitive, mesenchymal-like phenotype. The cell line is tumorigenic in immunocompromised mice. Exome sequencing identified a total of 45 and 76 somatic mutations in the primary tumor and cell line, respectively, of which 32 were identified in both samples and included mutations in known driver genes PIK3CA, TP53, and ARID1A. P7731 retains the DNA copy number alterations present in the matching primary tumor. Homozygous deletions detected in the cell line and in the primary tumor were found in regions containing three known (CDKN2A, CDKN2B, and CDKN1B) and 23 putative tumor suppressor genes. Cell line-specific gene amplification coupled with mRNA expression analysis revealed genes and pathways with potential pro-metastatic functions. CONCLUSION: This novel human breast cancer-bone metastasis cell line will be a useful model to study aspects of breast cancer biology, particularly metastasis-related changes from breast to bone.

Germline and somatic variant identification using BGISEQ-500 and HiSeq X Ten whole genome sequencing.

Technological innovation and increased affordability have contributed to the widespread adoption of genome sequencing technologies in biomedical research. In particular large cancer research consortia have embraced next generation sequencing, and have used the technology to define the somatic mutation landscape of multiple cancer types. These studies have primarily utilised the Illumina HiSeq platforms. In this study we performed whole genome sequencing of three malignant pleural mesothelioma and matched normal samples using a new platform, the BGISEQ-500, and compared the results obtained with Illumina HiSeq X Ten. Germline and somatic, single nucleotide variants and small insertions or deletions were independently identified from data aligned human genome reference. The BGISEQ-500 and HiSeq X Ten platforms showed high concordance for germline calls with genotypes from SNP arrays (>99%). The germline and somatic single nucleotide variants identified in both sequencing platforms were highly concordant (86% and 72% respectively). These results indicate the potential applicability of the BGISEQ-500 platform for the identification of somatic and germline single nucleotide variants by whole genome sequencing. The BGISEQ-500 datasets described here represent the first publicly-available cancer genome sequencing performed using this platform.

Corrigendum: Whole-genome landscape of pancreatic neuroendocrine tumours.

This corrects the article DOI: 10.1038/nature21063.