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.
Uveal melanoma (UM) is the most common intraocular tumour in adults and despite surgical or radiation treatment of primary tumours, ~50% of patients progress to metastatic disease. Therapeutic ...options for metastatic UM are limited, with clinical trials having little impact. Here we perform whole-genome sequencing (WGS) of 103 UM from all sites of the uveal tract (choroid, ciliary body, iris). While most UM have low tumour mutation burden (TMB), two subsets with high TMB are seen; one driven by germline MBD4 mutation, and another by ultraviolet radiation (UVR) exposure, which is restricted to iris UM. All but one tumour have a known UM driver gene mutation (GNAQ, GNA11, BAP1, PLCB4, CYSLTR2, SF3B1, EIF1AX). We identify three other significantly mutated genes (TP53, RPL5 and CENPE).
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.
Abstract
Summary
Identification of duplicate templates is a common preprocessing step in bulk sequence analysis; for large libraries, this can be resource intensive. Here, we present streammd: a ...fast, memory-efficient, single-pass duplicate marker operating on the principle of a Bloom filter. streammd closely reproduces outputs from Picard MarkDuplicates while being substantially faster, and requires much less memory than SAMBLASTER.
Availability and implementation
streammd is a C++ program available from GitHub https://github.com/delocalizer/streammd under the MIT license.
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.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The clinical classification of variants may change with new information, however, there is limited guidance on how often significant changes in variant classification occur. We used ClinVar to ...examine how variant classification changes over time. We developed a custom parser and accessed variant data from ClinVar between January 2015 and July 2021. The ClinVar‐assigned “aggregate” classification of variants in 121 hereditary cancer genes was harmonized across releases to align to the American College of Medical Genetics and Genomics and the Association for Molecular Pathology terms. Aggregate classification categories were grouped as: benign/likely benign (B/LB); likely pathogenic/pathogenic (LP/P); variant of uncertain significance (VUS); conflicting interpretations of pathogenicity (Conflicting); or Other. We profiled changes in aggregate variant classification between consecutive semi‐annual ClinVar releases. The proportion of variants that changed aggregate classification between semi‐annual ClinVar releases ranged from 0.6% to 6.4%. The most frequent changes were “VUS to conflicting,” “other to LP/P,” and “B/LB to Conflicting.” A limited number of variants changed aggregate classification from “LP/P to B/LB,” or vice versa. Our analysis indicates need for regular reassessment of clinical variant interpretations. The parser developed for this project will facilitate extraction of relevant interpretation data from ClinVar.
Modelling changes in classification of variants in cancer genes over time
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
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.
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.
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.
We propose accounting for genomic and immune microenvironment status may influence therapeutic planning in the future.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Tumour cellularity, the relative proportion of tumour and normal cells in a sample, affects the sensitivity of mutation detection, copy number analysis, cancer gene expression and methylation ...profiling. Tumour cellularity is traditionally estimated by pathological review of sectioned specimens; however this method is both subjective and prone to error due to heterogeneity within lesions and cellularity differences between the sample viewed during pathological review and tissue used for research purposes. In this paper we describe a statistical model to estimate tumour cellularity from SNP array profiles of paired tumour and normal samples using shifts in SNP allele frequency at regions of loss of heterozygosity (LOH) in the tumour. We also provide qpure, a software implementation of the method. Our experiments showed that there is a medium correlation 0.42 (Formula: see text-value=0.0001) between tumor cellularity estimated by qpure and pathology review. Interestingly there is a high correlation 0.87 (Formula: see text-value Formula: see text 2.2e-16) between cellularity estimates by qpure and deep Ion Torrent sequencing of known somatic KRAS mutations; and a weaker correlation 0.32 (Formula: see text-value=0.004) between IonTorrent sequencing and pathology review. This suggests that qpure may be a more accurate predictor of tumour cellularity than pathology review. qpure can be downloaded from https://sourceforge.net/projects/qpure/.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The importance of epigenetic modifications such as DNA methylation in tumorigenesis is increasingly being appreciated. To define the genome‐wide pattern of DNA methylation in pancreatic ductal ...adenocarcinomas (PDAC), we captured the methylation profiles of 167 untreated resected PDACs and compared them to a panel of 29 adjacent nontransformed pancreata using high‐density arrays. A total of 11,634 CpG sites associated with 3,522 genes were significantly differentially methylated (DM) in PDAC and were capable of segregating PDAC from non‐malignant pancreas, regardless of tumor cellularity. As expected, PDAC hypermethylation was most prevalent in the 5′ region of genes (including the proximal promoter, 5′UTR and CpG islands). Approximately 33% DM genes showed significant inverse correlation with mRNA expression levels. Pathway analysis revealed an enrichment of aberrantly methylated genes involved in key molecular mechanisms important to PDAC: TGF‐β, WNT, integrin signaling, cell adhesion, stellate cell activation and axon guidance. Given the recent discovery that SLIT‐ROBO mutations play a clinically important role in PDAC, the role of epigenetic perturbation of axon guidance was pursued in more detail. Bisulfite amplicon deep sequencing and qRT‐PCR expression analyses confirmed recurrent perturbation of axon guidance pathway genes SLIT2, SLIT3, ROBO1, ROBO3, ITGA2 and MET and suggests epigenetic suppression of SLIT‐ROBO signaling and up‐regulation of MET and ITGA2 expression. Hypomethylation of MET and ITGA2 correlated with high gene expression, which was associated with poor survival. These data suggest that aberrant methylation plays an important role in pancreatic carcinogenesis affecting core signaling pathways with potential implications for the disease pathophysiology and therapy.
What's new?
Based on a large genome‐wide scan of DNA methylation, this study reports that global DNA methylation patterns can robustly segregate tumor and non‐malignant pancreata. Cancer methylation also affects key pathways in pancreatic carcinogenesis, including TGF‐β, WNT, and axon guidance signaling. This study confirms that methylation plays an important role in the development and progression of pancreatic cancer, with implications for both ongoing research and therapeutic development.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK