COSMIC Database Research Articles

Abstract 4982: The somatic mutation signature detected in gynecological mucosal melanoma

Abstract Background: Gynecological Mucosal Melanoma(GMM) is aggressive disease with a poor prognosis, where few studies focusing on its target drug and mutation signatures. Cancer as outcome of DNA mutagenesis process and can be inferred from somatic mutations signature by analyzing the genome sequence. Due to the rapid development of next-generation sequencing (NGS) and molecularly driven targeted therapy, discovering different mutation signature in cancer plays an important role in precision medicine. Methods: The study was conducted at a high-volume melanoma center, the Peking University Cancer Hospital. We enrolled female patients with melanoma (n = 22) from 2017 to 2021. Subsequently, whole-exome sequencing performed on the Formalin-fixed paraffin-embedded (FFPE) specimens from 22 GMM patients and match normal tissues of patients. Tumor mutational burden (TMB) was measured as the number of nonsynonymous somatic mutations of whole-exome sequence coding region. Mutation signature detection base pair of six substitution subtypes: C>A, C>G, C>T, T>A, T>C, and T>G. NMF can be used to decompose matrix A into two nonnegative matrices. Results: In all mutations of 22 GMM patients, we screen two mutation signatures. Signature 1 showed a strong correlation [0.77;cosine similarity value (CSV)] to the signature 3 from the COSMIC database, which represents failure of DNA double-strand break-repair by homologous recombination, which predominated in 55%(12/22) GMM patients. Signature 2 showed a strong correlation to Cosmic signature 1 (CSV = 0.791), which is has been found in all cancer types correlated with age cancer diagnosis. Patient-18 did not any treatment in our hospital, follow-up was lost. Kaplan-Meier revealed that Signature 2 was associated with longer overall survival in GMM patients and later recurrence or progression. Conclusions: This study aimed to identify prognostic markers and the development of targeted therapy for the treatment of GMM. Citation Format: Zhongwu Li, Lianyi Shen, Siyao Liu. The somatic mutation signature detected in gynecological mucosal melanoma. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4982.

Comprehensive genomic profiling of upper tract urothelial carcinoma and urothelial carcinoma of the bladder identifies distinct molecular characterizations with potential implications for targeted therapy & immunotherapy.

Despite the genomic landscape of urothelial carcinomas (UC) patients, especially those with UC of bladder (UCB), has been comprehensively delineated and associated with pathogenetic mechanisms and treatment preferences, the genomic characterization of upper tract UC (UTUC) has yet to be fully elucidated. A total of 131 Chinese UTUC (74 renal pelvis & 57 ureter) and 118 UCB patients were enrolled in the present study, and targeted next-generation sequencing (NGS) of 618 cancer-associated genes were conducted to exhibit the profile of somatic and germline alterations. The COSMIC database, including 30 mutational signatures, were utilized to evaluate the mutational spectrums. Moreover, TCGA-UCB, MSKCC-UCB, and MSKCC-UTUC datasets were retrieved for preforming genomic alterations (GAs) comparison analysis between Western and Chinese UC patients. In our cohort, 93.98% and 56.63% of UC patients were identified with oncogenic and actionable somatic alterations, respectively. Meanwhile, 11.24% of Chinese UC patients (of 14.50% and 7.63% of UTUC and UCB cases, respectively) were identified to harbor a total of 32 pathogenic/likely-pathogenic germline variants in 22 genes, with DNA damage repair (DDR)-associated BRCA1 (1.20%) and CHEK2 (1.20%) being the most prevalent. Chinese UTUC and UCB patients possessed distinct somatic genomic characteristics, especially with significantly different prevalence in KMT2D/C/A, GNAQ, ERCC2, RB1, and PPM1D. In addition, we also found notable differences in the prevalence of ELF3, TP53, PMS2, and FAT4 between renal pelvis and ureter carcinomas. Moreover, 22.90% and 33.90% of UTUC and UCB patients, respectively, had at least one deleterious/likely deleterious alteration in DDR related genes/pathways. Subsequently, mutational signature analysis revealed that UC patients with mutational signature 22, irrespective of UTUC or UCB, consistently had the markedly higher level of tumor mutational burden (TMB), which was proved to be positively correlated with the objective complete/partial response rate in the IMvigor210 cohort. By comparison, Chinese and Western UTUC patients also differed regrading GAs in oncogenic-related genes/pathways, especially in TP53, RTK/RAS, and PI3K pathways; besides, more alterations in WNT pathway but less TP53, RTK/RAS, HIPPO, and PI3K pathways were identified in Chinese UCB. The in-depth analysis of genomic mutational landscapes revealed distinct pathogenetic mechanisms between Chinese UTUC and UCB, and specific genomic characterizations could identify high risk population of UTUC/UCB and provided information regarding the selection of alternative therapeutic regimens.

Whole-exome sequencing reveals genetic variants in low-risk and high-risk neuroblastoma

This study aimed to investigate the genetic aberrations in neuroblastoma (NB) by comparing high and low-risk NB patients by whole-exome sequencing (WES) and to reveal the heterogeneity and association between somatic variants and clinical features. Seven NB patients with available clinical data were included in the study (4 in the low-risk group and 3 in the high-risk group). WES was performed and somatic variants associated with NB genes in the COSMIC database were selected through bioinformatics pipeline analysis. Variants were determined using the Integrative Genomics Viewer (IGV). Some gene variations were found in both groups, including variations in oncogene and tumor suppressor genes. In general, candidate gene variations were associated with chromatin remodeling complexes, the RAS pathway, cell proliferation, and DNA repair mechanism. Some variations in CSF1R, MSH6, PTPN11, SOX9, RET, TSC1, and DNMT1 genes were detected only in high-risk patients, while EP300, TET2, MYCN, PRDM1, and ARID2 gene variations were detected only in low-risk patients. When high-risk gene variants were compared with the cBioportal cancer genomic database, two common gene variants (ARID1A and NCOR2) were identified. However, when low-risk gene variants were compared with the cBioportal cancer genomic database, no common genes were found. GO/KEGG enrichment analysis was performed to find relevant biological processes and molecular pathways related to gene variants, which will help to decipher the molecular mechanisms of NB tumorigenesis and the phenotypic differences between high-risk and low-risk patients.

Analysis Of Common Somatic Mutations In Colorectal Carcinoma And Associated Dysregulated Pathwaysarts.

Identification of gene targets and biological pathways involved in colorectal carcinoma (CRC) is essential for better management of patients. Our study aims to highlight common somatic mutations in colorectal carcinoma and to identify dysregulated pathways and gene enrichment based on KRAS and BRAF interaction network analysis. By using cancer browser tool in COSMIC database, mutation frequencies of the top 20 mutated genes listed for colorectal adenocarcinoma were identified. The most frequent variants of selected genes were explored with ClinVar database which led to identification of protein change along with its cytogenic location, variant type, variant length and the associated single nucleotide polymorphism (SNP). These identified SNPs were searched in Pakistani database using 1000genome in an attempt to identify common polymorphisms. Using the database ClinicalTrial.gov the number of clinical trials based upon these selected mutations was explored. Enrichment and protein interaction (PI) analysis of KRAS and BRAF was carried out to reveal significant biological pathways associated with these genes. In cumulative data, among all variants about 57% of substitution mutations are observed to be G>A including mutations in KRAS, Tp53, SMAD4, PI3K and NRAS. The mutations of KRAS (c.35G>A), TP53 (c.524G>A) and APC (c.4348C>T) were found to be pathogenic with single nucleotide variation and variant length of 1bp. Searching 1000genome database revealed that 100 % of alleles found in East Asian population studied are 'C'(frequency=1). Significant biological pathways (<0.05) identified by our search include Trk receptor signalling mediated by the MAPK pathway, signalling to p38 via RIT and RIN, signalling to ERKs, Frs2-mediated activation, ARMS-mediated activation and prolonged ERK activation events. Our study highlights the role of genetic profiling in CRC, with emphasis on mutations which may define treatment outcome. Targeting several collateral pathways simultaneously may be further explored to improve colorectal cancer therapeutics.

Abstract PR005: TP53 somatic evolution in healthy fallopian tube across the full spectrum of human lifespan

Abstract Introduction: TP53 is the most mutated gene in human cancers but research in recent years has demonstrated that TP53 mutations are also prevalent in healthy tissue. Little is known, however, about how these mutations accumulate with age and how this process might differ between and within individuals. This gap in knowledge is largely due to the difficulty of measuring clonal expansions in normal tissue since mutant clones are typically small and escape the detection of standard sequencing technologies. We previously demonstrated the high sensitivity of duplex sequencing (duplex-seq) to detect low frequency mutations in normal tissue. Here we use duplex-seq to characterize TP53 mutations across the full range of human lifespan in fallopian tube epithelium. Fallopian tube is the site of origin of high-grade serous ovarian carcinoma, a very aggressive cancer type that is almost exclusively driven by TP53 mutations. We hypothesize that TP53 clonal expansions increase in number, size, and pathogenicity during human lifespan, especially in the fimbria, which is the distal part of the tube and the predominant site for cancer formation. Methods: Frozen fallopian tube was collected at autopsy from 84 individuals spanning the full range of human lifespan (from newborn to &amp;gt;90-year-old). For most individuals, collection included a sample from the fallopian tube fimbria (distal part) and ampulla (proximal part). For a subset of patients, both right and left fallopian tubes were sampled. Frozen OCT embedded tissue was sectioned, pathologically examined to determine the presence of fallopian tube epithelium, and macrodissected. TP53 mutations were identified by ultra-deep duplex-seq (12,000x) and were compared between fallopian tube fimbria and ampulla, between right and left tubes, across decades of age, and with cancer mutations reported in the COSMIC database. Results: About 11% of samples have been analyzed, and the full cohort is expected to be completed within the next two months. TP53 mutations were identified in all individuals at an average of 6 mutations per sample (min=0, max=15). The number of mutations significantly increased with age and, in this smaller subset, did not significantly differ between bilateral samples or between fimbria and ampulla. The size of the mutant clones also increased with age. Most mutations were missense, but nonsense, indels, and splice mutations were also detected, especially at older ages. More than half of the mutations were predicted to be pathogenic and coincided with mutations commonly found in ovarian cancer. Conclusions: Ultra-deep sequence has enabled unprecedented resolution to characterize TP53 somatic evolution in normal tissue. TP53 mutant clones accumulate in the fallopian tube with age and are enriched for pathogenic mutations, indicating clonal expansion upon selective advantage. Deep characterization of these mutant clones will inform our knowledge of the normal physiological process of somatic evolution during healthy aging and its potential contribution to cancer formation later in life. Citation Format: Rosana Risques, Coohleen Coombes, Jeanne Fredrickson, Elizabeth U. Parker, Brendan Kohrn, Elena S. Latorre-Esteves, Emma Hazard, Eden M. Milan, Ronit Katz. TP53 somatic evolution in healthy fallopian tube across the full spectrum of human lifespan [abstract]. In: Proceedings of the AACR Special Conference: Aging and Cancer; 2022 Nov 17-20; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2022;83(2 Suppl_1):Abstract nr PR005.

Whole-exome sequencing and bioinformatic analyses revealed differences in gene mutation profiles in papillary thyroid cancer patients with and without benign thyroid goitre background.

Papillary thyroid cancer (PTC) is the most common thyroid malignancy. Concurrent presence of cytomorphological benign thyroid goitre (BTG) and PTC lesion is often detected. Aberrant protein profiles were previously reported in patients with and without BTG cytomorphological background. This study aimed to evaluate gene mutation profiles to further understand the molecular mechanism underlying BTG, PTC without BTG background and PTC with BTG background. Patients were grouped according to the histopathological examination results: (i) BTG patients (n = 9), (ii) PTC patients without BTG background (PTCa, n = 8), and (iii) PTC patients with BTG background (PTCb, n = 5). Whole-exome sequencing (WES) was performed on genomic DNA extracted from thyroid tissue specimens. Nonsynonymous and splice-site variants with MAF of ≤ 1% in the 1000 Genomes Project were subjected to principal component analysis (PCA). PTC-specific SNVs were filtered against OncoKB and COSMIC while novel SNVs were screened through dbSNP and COSMIC databases. Functional impacts of the SNVs were predicted using PolyPhen-2 and SIFT. Protein-protein interaction (PPI) enrichment of the tumour-related genes was analysed using Metascape and MCODE algorithm. PCA plots showed distinctive SNV profiles among the three groups. OncoKB and COSMIC database screening identified 36 tumour-related genes including BRCA2 and FANCD2 in all groups. BRAF and 19 additional genes were found only in PTCa and PTCb. "Pathways in cancer", "DNA repair" and "Fanconi anaemia pathway" were among the top networks shared by all groups. However, signalling pathways related to tyrosine kinases were the most significantly enriched in PTCa while "Jak-STAT signalling pathway" and "Notch signalling pathway" were the only significantly enriched in PTCb. Ten SNVs were PTC-specific of which two were novel; DCTN1 c.2786C>G (p.Ala929Gly) and TRRAP c.8735G>C (p.Ser2912Thr). Four out of the ten SNVs were unique to PTCa. Distinctive gene mutation patterns detected in this study corroborated the previous protein profile findings. We hypothesised that the PTCa and PTCb subtypes differed in the underlying molecular mechanisms involving tyrosine kinase, Jak-STAT and Notch signalling pathways. The potential applications of the SNVs in differentiating the benign from the PTC subtypes requires further validation in a larger sample size.

Investigation of the Potential Correlation Between RNA-Binding Proteins in the Evolutionarily Conserved MEX3 Family and Non-small-Cell Lung Cancer

Members of the MEX3 (muscle excess 3) family, uniquely characterised as mRNA binding proteins, play emerging roles in the post-transcriptional regulation of programmed biological processes, including tumour cell death and immune mechanisms, and have been shown to be involved in a variety of diseases. However, the role of MEX3 in non-small cell lung cancer (NSCLC) has not been fully elucidated. In this study, we found no significant changes in the sequence and copy number of the MEX3 gene through analysis using the COSMIC database, revealing its stability during malignancy development. Its expression in NSCLC was examined using the Oncomine™ database, and the prognosis of each member gene was analysed by Kaplan–Meier. The results showed that overexpression of MEX3A, MEX3B, MEX3C and MEX3D was associated with significantly worse OS in patients with LUAD, while overexpression of MEX3D was also associated with significantly worse OS in patients with LUSC. Afterwards, we applied the Tumour Immunology Estimation Resource (TIMER) tool to assess the correlation between different MEX3 and infiltrative immune cell infiltration. Ultimately, we found that most MEX3 members were highly expressed in NSCLC, with high expression suggesting poor prognosis and correlating with immune cell infiltration. The complexity and heterogeneity of NSCLC was understood through MEX3, setting the framework for the prognostic impact of MEX3 in NSCLC patients and the development of new targeted therapeutic strategies in the future.

Abstract A021: Clinical and molecular characterization of naturally-occurring colorectal cancer in rhesus macaques reveals mismatch repair deficiency driven by epigenetic silencing of MLH1

Abstract Naturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and emerge at similar life stages as in human patients, but the carcinogenesis and molecular background remain elusive, limiting the value of the model for human disease and cancer therapy. In order to better understand this naturally occurring model, we established a diagnostic and staging pipeline including imaging (18F-FDG PET, plain and contrast-enhanced CT), histopathology, and clinical assessments. We observed cancer symptoms and co-morbidities in our cohort (n=16, all Indian-origin, 11 females, mean age at arrival 20.2y) such as hypoalbuminemia, fecal occult blood, and microcytic anemia, as frequently described in human CRC patients. Pathologically, all cancers were right-sided, involving the proximal colon and/or ileocecocolic junction, and most were densely fibrotic, restricting the colonic lumen. Most cancers appeared with glandular morphology and some (18.8%) had mucinous components. Immunohistochemistry revealed loss of MLH1 and PMS2 in 100% of investigated CRCs, indicating mismatch repair deficiency which furthermore resulted in microsatellite instability (PCR &amp; fragment analysis). Whole exome sequencing revealed the close genetic relatedness to human CRCs, particularly exemplified by mutations affecting KRAS (37.5%, e.g., p.G12D), APC (31.3%), TP53 (18.8%, e.g. p.R175H), ARID1A (56.3%), and ALK (43.8%), as similarly annotated in the human COSMIC database. Nonetheless, somatic mutations do not explain the loss of MLH1 in the entire CRC cohort. Transcriptomics on the other hand revealed the transcriptional suppression of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRC in comparison to adjacent healthy colon. Moreover, comparison of differentially expressed gene sets of rhesus CRC with a human annotated database (IPA) confirmed the disease similarities observed clinically, genetically, and histopathologically. Subsequently, we investigated DNA methylation of the promoter region of MLH1 and retranslated markers for the CpG island methylator phenotype (CIMP) as described in human CRC. While only 56.3% of CRCs were considered CIMP positive (≥3/5 markers hypermethylated), 100% of investigated CRCs exhibited MLH1 promoter hypermethylation. As a result, epigenetic silencing is suggested to suppress MLH1 transcription, cause the loss of MLH1 protein, and drive mismatch repair deficiency and genomic instability in naturally occurring CRC in rhesus macaques. We therefore consider spontaneous, uninduced CRC in rhesus macaques, their treatment-naïve nature, and their unaltered immune competence an outstanding model for human disease and in particular for human cancer immunotherapy. Citation Format: Simon Deycmar, Brendan Johnson, Karina Ray, David Caudell Caudell, John Olson, Greg Dugan, W. Shane Sills, Declan Ryan, Christopher Whitlow, Kiran K. Solingapuram Sai, Betsy Ferguson, Benjamin Bimber, Cassandra Cullin, Brandy Dozier, Emily Romero, Francois Villinger, Armando Burgos, Jeff Chou, Bruno Gomes, Michael Hettich, Maurizio Ceppi, Jehad Charo, J. Mark Cline. Clinical and molecular characterization of naturally-occurring colorectal cancer in rhesus macaques reveals mismatch repair deficiency driven by epigenetic silencing of MLH1 [abstract]. In: Proceedings of the AACR Special Conference on Colorectal Cancer; 2022 Oct 1-4; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_1):Abstract nr A021.

In Silico and Structure-Based Assessment of Similar Variants Discovered in Tandem Repeats of BRCT Domains of BRCA1 and BARD1 To Characterize the Folding Pattern.

BRCA1 and BARD1 are important proteins in the homologous DNA damage repair pathways. Different genetic variants identified in these proteins have been clinically correlated with the occurrence of hereditary breast and ovarian cancer (HBOC). Variants of unknown significance (VUS) reported in the BRCT domains of BRCA1 and BARD1 substantiate the importance of BRCT domain-containing proteins for genomic integrity. To classify the pathogenicity of variants, in silico, structural and molecular dynamics (MD)-based approaches were explored. Different variants reported in the BRCT region were retrieved from cBioPortal, LOVD3, BRCA Exchange, and COSMIC databases to evaluate the pathogenicity. Multiple sequence alignment and superimposition of the structures of BRCA1 BRCT and BARD1 BRCT domains were performed to compare alterations in folding patterns. From 11 in silico predictions servers, variants reported to be pathogenic by 70% of the servers were considered for structural analysis. To our observations, four residue pairs of both the proteins were reported, harboring 11 variants, H1686Y, W1718L, P1749L, P1749S, and W1837L variants for BRCA1 BRCT and H606D, H606N, W635L, P657L, P657S, and W762F for BARD1 BRCT. MD simulations of the BRCT repeat regions of these variants and wild-type proteins were performed to evaluate the differences of folding patterns. Root mean square deviation (RMSD), R g, solvent-accessible surface area (SASA), and root mean square fluctuation (RMSF) of variants showed slight differences in the folding patterns from the wild-type proteins. Furthermore, principal components analysis of H1686Y, P1749S, and W1718L variants of BRCA1 showed less flexibility than the wild type, whereas that of H606D, W635L, and W762F of BARD1 showed more flexibility than the wild type. Normal mode analysis of the energy minima from the simulation trajectories revealed that most of the variants do not show much differences in the flexibility compared to the wild-type proteins, except for the discrete regions in the BRCT repeats, most prominently in the 1798-1801 amino acid region of BRCA1 and at the residue 744 in BARD1.

Identification of potential molecular players and therapeutic drug molecules in Melphalan resistant Multiple myeloma by integrated bioinformatics analysis

Multiple myeloma (MM), second most common haematological malignancy, still remains incurable with the acquiring of drug resistance. Melphalan therapy which has been used as one of the key therapies against MM, is impeded by the occurrence of Melphalan resistance. The precise underlying mechanism of this acquired Melphalan resistance in MM is not yet well deciphered. Therefore, this present study was aimed to identify the differentially expressed genes (DEGs) associated micro RNAs (miRNAs) and transcription factors (TFs) from the microarray dataset of Melphalan resistant and Melphalan sensitive MM cell lines, obtained from Gene Expression Omnibus (GEO) database. DEGs were analysed using GEO2R tool from the dataset. Then the gene ontology (GO) and pathway enrichment analysis were executed by using DAVID database. Protein-protein interaction (PPI) network of DEGs was constructed and analysed by using STRING database and Cytoscape tool. Genetic alterations in DEGs were studied through COSMIC database. Network of interaction of DEGs and miRNAs as well as TFs was constructed and analysed by using mirDIP, TRRUST and miRNet tools. Drug gene interaction was examined to identify potential drug molecules by DGIdb tool. Various DEGs that might play pivotal role in Melphalan resistant MM, were detected and selected for further analysis. miRNA analysis detected hsa-mir- 21-3p, hsa-mir-27a-5p, hsa-mir-129-2-3p that could interact with maximum target DEGs. One TF, STAT1 was found to regulate the expression of selected DEGs. The whole study may give a better understanding about the Melphalan resistance in MM.

Effect of naturally-occurring mutations on the stability and function of cancer-associated NQO1: Comparison of experiments and computation.

Recent advances in DNA sequencing technologies are revealing a large individual variability of the human genome. Our capacity to establish genotype-phenotype correlations in such large-scale is, however, limited. This task is particularly challenging due to the multifunctional nature of many proteins. Here we describe an extensive analysis of the stability and function of naturally-occurring variants (found in the COSMIC and gnomAD databases) of the cancer-associated human NAD(P)H:quinone oxidoreductase 1 (NQO1). First, we performed in silico saturation mutagenesis studies (>5,000 substitutions) aimed to identify regions in NQO1 important for stability and function. We then experimentally characterized twenty-two naturally-occurring variants in terms of protein levels during bacterial expression, solubility, thermal stability, and coenzyme binding. These studies showed a good overall correlation between experimental analysis and computational predictions; also the magnitude of the effects of the substitutions are similarly distributed in variants from the COSMIC and gnomAD databases. Outliers in these experimental-computational genotype-phenotype correlations remain, and we discuss these on the grounds and limitations of our approaches. Our work represents a further step to characterize the mutational landscape of NQO1 in the human genome and may help to improve high-throughput in silico tools for genotype-phenotype correlations in this multifunctional protein associated with disease.

High Prevalence of Epigenetic Mutations in Histiocytic and Dendritic Cell Disorders: Results from Molecular Analysis of a Large Cohort from Histiocytosis Working Group

Introduction Over the last decade, recurrent oncogenic alterations in the MAPK/ERK pathway have been discovered in various histiocytic and dendritic cell disorders, especially Erdheim-Chester disease (ECD), Langerhans cell histiocytosis (LCH), and Rosai-Dorfman disease (RDD). These discoveries have led to the successful use of targeted BRAF- and MEK-inhibitors. However, most data are derived from small institutional studies mostly focusing on ECD and LCH, which leaves room for a comprehensive evaluation of the molecular underpinnings in a large cohort. Methods Formalin-fixed paraffin embedded tumor biopsies from various histiocytic/dendritic cell disorders seen at a tertiary institution from 2001-2021 were subjected to whole exome sequencing (WES) using one of two assays: in-house SureSelectXT AllExon V5+UTR or commercially available Tempus xT assay. Identified variants were filtered to remove polymorphic, low coding impact and low read coverage variants. A variant was then selected for further analysis if it met one of the following criteria: 1) observed as somatic in ≥1 tumor-germline pairs or COSMIC database or 2) observed in at least three samples with an allele frequency of ≥10% in ≥1 sample. Cases were clustered into mutually exclusive groups based on their primary driver mutational profile using the following funneling order: BRAF-V600E (BRAF-V600E mutations), BRAF*/NRAS/KRAS (BRAF non-V600E, KRAS, or NRAS mutations), CSF1R (CSF1R mutations), MAPK (MAPK/MEK mutations), and PI3K (AKT3, INPP4B, MTOR, PIK3CA, PIK3R1, PIK3R2, PIK3R3, PPP2R1A, MTR mutations). A cancer specific pathway enrichment was performed on the unknown cases, which identified epigenetic mutations as major pathway alterations observed in the cohort. We also correlated the presence of epigenetic mutations with clinical features and outcomes. Results We included 114 cases with histiocytic and dendritic cell disorders in the study. Median age at diagnosis was 51y (IQR 37-64). The disease distribution was as follows: 46 LCH, 20 ECD, 17 RDD, 8 follicular dendritic cell sarcoma (FDCS), 1 adult xanthogranuloma, 12 malignant histiocytosis (MH; including histiocytic and Langerhans cell sarcoma), and 10 multicentric reticulohistiocytosis (MRH). Based on the known driver mutations in histiocytoses, the following 5 molecular groups were identified: BRAF-V600E (n= 19, 17%), BRAF*/KRAS/NRAS (n=17, 15%), MAPK (n=14, 12%), PI3K (n=8, 7%), CSF1R (n=2, 2%) (Figure 1). A large proportion (n=48, 42%) did not have any of the known driver mutations, while 6 cases (5%) had no alterations identified. Cancer-specific enrichment analysis in the entire cohort revealed mutations in epigenetic genes in 47 (41%) cases, including 30 (26%) cases with identifiable driver mutations and 17 (15%) cases without known driver mutations. Median number of epigenetic mutations was 1 (range 1-11); 21% cases had 2 and 10% cases had 3 epigenetic mutations, respectively. The most common epigenetic mutations included those in KMT2C (34%), KMT2D (23%), and ARID1A (19%). The mutational groups with highest prevalence of epigenetic mutations were CSF1R (100%), MAPK (71%), and PI3K (63%), while lowest prevalence of these alterations was found in BRAF-V600E group (16%; p=0.0003). The epigenetic gene mutations were seen more commonly in MRH (90%) and MH (50%) compared with ECD (45%), RDD (41%), FDCS (38%), and LCH (33%; p=0.037). There was no association between epigenetic gene mutations and age at diagnosis or organ involvement. The median overall survival (OS) was not significantly different between those with epigenetic mutations and those without such alterations (22y vs. 17y, p=0.89). Conclusions In this large study of patients with diverse histiocytic and dendritic cell disorders, we found a high prevalence of mutations in epigenetic genes (41%). These mutations were specifically enriched in tumors with non-BRAF-V600E mutations, and in MRH/MH. These data shed light on additional mechanisms of disease pathogenesis, and with available epigenetic therapies, warrant further investigation. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Whole exome sequencing identified a novel POT1 variant as a candidate pathogenic allele underlying a Li-Fraumeni-like family.

Li-Fraumeni syndrome (LFS) and Li-Fraumeni-like (LFL) syndrome are rare hereditary diseases characterized by predisposition to a diverse spectrum of cancer types, primarily sarcoma. The pathogenic variants underlying the majority of LFL cases remain to be explored. We performed whole-exome sequencing (WES) on 13 core members of a large LFL family with highly aggregated incidences of cancers, including cases with sarcoma, non-small cell lung cancer and cardiac angiosarcoma, and conducted a comprehensive literature review of candidate gene associations in LFS/LFL syndromes or sarcoma to identify potential pathogenic germline variants. No germline variants in the best-known LFL/LFS-associated gene TP53 were detected. Of all the genes associated with LFS/LFL or sarcoma that we have surveyed, we identified a novel p.P35L germline variant in POT1 (protection of telomeres 1). Germline and somatic alterations in POT1 have been implicated in a series of familial cancers, including angiosarcoma, glioma, melanoma and colorectal cancer. This particular variant is located in the telomere-binding OB1 domain, which is important in maintaining the proper telomere length, and showed high conservation across different POT1 orthologues. No record of the variant was found in any of the 1000 genomes, ExAC, gnomAD, dpSNP and COSMIC databases. Prediction algorithms and in silico structural analysis suggested completely disrupted protein structure and function of POT1 in the presence of this mutation. Leveraging WES, we identified a novel germline risk allele, p.P35L in POT1, that likely predisposes to LFL syndrome. Our results support the routine testing of POT1 and other LFL/LFS-associated genes in the risk populations to enable early cancer diagnosis, prevention and intervention.

Cancer-associated Notch receptor variants lead to O-fucosylation defects that deregulate Notch signaling

NOTCH1 is a transmembrane receptor that initiates a signaling pathway involved in embryonic development of adult tissue homeostasis. The extracellular domain of NOTCH1 is composed largely of epidermal growth factor-like repeats (EGFs), many of which can be O-fucosylated at a specific consensus sequence by protein O-fucosyltransferase 1 (POFUT1). O-fucosylation of NOTCH1 is necessary for its function. The Notch pathway is deregulated in many cancers, and alteration of POFUT1 has been reported in several cancers, but further investigation is needed to assess whether there is deregulation of the Notch pathway associated with mutations that affect O-fucosylation in cancers. Using Biomuta and COSMIC databases, we selected nine NOTCH1 variants that could cause a change in O-fucosylation of key EGFs. Mass spectral glycoproteomic site mapping was used to identify alterations in O-fucosylation of EGFs containing the mutations. Cell-based NOTCH-1 signaling assays, ligand-binding assays, and cellsurface analysis were used to determine the effect of each mutation on Notch activation. Two variants led to a gain of function (GOF), six to a loss of function (LOF), and one had minimal effects. Most GOF and LOF were associated with a change in O-fucosylation. Finally, by comparing our results with known NOTCH1 alterations in cancers from which our mutations originated, we were able to establish a correlation between our results and the known GOF or LOF of NOTCH1 in these cancers. This study shows that point mutations in N1 can lead to alterations in O-fucosylation that deregulate the Notch pathway and be associated with cancer processes.

Towards Precision Oncology: The Role of Smoothened and Its Variants in Cancer.

The G protein-coupled receptor Smoothened (Smo) is a central signal transducer of the Hedgehog (Hh) pathway which has been linked to diverse forms of tumours. Stimulated by advancements in structural and functional characterisation, the Smo receptor has been recognised as an important therapeutic target in Hh-driven cancers, and several Smo inhibitors have now been approved for cancer therapy. This receptor is also known to be an oncoprotein itself and its gain-of-function variants have been associated with skin, brain, and liver cancers. According to the COSMIC database, oncogenic mutations of Smo have been identified in various other tumours, although their oncogenic effect remains unknown in these tissues. Drug resistance is a common challenge in cancer therapies targeting Smo, and data analysis shows that healthy individuals also harbour resistance mutations. Based on the importance of Smo in cancer progression and the high incidence of resistance towards Smo inhibitors, this review suggests that detection of Smo variants through tumour profiling could lead to increased precision and improved outcomes of anti-cancer treatments.

Comparative study of the evolution of cancer gene duplications across fish.

Comparative studies of cancer-related genes not only provide novel information about their evolution and function but also an understanding of cancer as a driving force in biological systems and species' life histories. So far, these studies have focused on mammals. Here, we provide the first comparative study of cancer-related gene copy number variation in fish. Fishes are a paraphyletic group whose last common ancestor is also an ancestor of the tetrapods, and accordingly, their tumour suppression mechanisms should include most of the mammalian mechanisms and also reveal novel (but potentially phylogenetically older) previously undetected mechanisms. We have matched the sequenced genomes of 65 fish species from the Ensemble database with the cancer gene information from the COSMIC database. By calculating the number of gene copies across species using the Ensembl CAFE data (providing species trees for gene copy number counts), we used a less resource-demanding method for homolog identification. Our analysis demonstrates a masked relationship between cancer-related gene copy number variation (CNV) and maximum lifespan in fish species, suggesting that a higher number of copies of tumour suppressor genes lengthens and the number of copies of oncogenes shortens lifespan. Based on the positive correlation between the number of copies of tumour suppressors and oncogenes, we show which species have more tumour suppressors in relation to oncogenes. It could be suggested that these species have stronger genetic defences against oncogenic processes. Fish studies could be a largely unexplored treasure trove for understanding the evolution and ecology of cancer, providing novel insights into the study of cancer and tumour suppression, in addition to fish evolution, life-history trade-offs, and ecology.

RAB6B is a potential prognostic marker and correlated with the remolding of tumor immune microenvironment in hepatocellular carcinoma.

Backgrounds: Hepatocellular carcinoma (HCC) is the most common type of primary liver cancer and the second leading cause of death among all cancers. The Ras-associated binding (Rab) proteins constitute the largest family of the Ras superfamily of small GTPases, which mainly mediate membrane trafficking processes. RAB6B is a member of Rab GTPases, and it has been found to be dysregulated in various tumors. However, the clinical significance, correlations with immune cells, and stroma infiltration of RAB6B in HCC remain unclear. Methods: RAB6B mRNA and protein expression in HCC were examined using the TIMER, HCCDB, UALCAN, and HPA databases. The genetic alterations of RAB6B were analyzed by cBioPortal and COSMIC databases. The correlations between RAB6B and tumor-infiltrating immune cells and cancer-associated fibroblasts were explored by using TIMER, TISIDB, and GEPIA databases. Co-expression networks of RAB6B were investigated based on LinkedOmics. Drug sensitivity was analyzed through the GDSC and CTRP databases. RAB6B was knocked down with siRNA in HCC cell lines. EdU assay was performed to detect the cell proliferation ability, flow cytometry was used to compare the differences in the ability of apoptosis, and MTT was used to evaluate the drug sensitivity in vitro. Results: RAB6B mRNA and protein expression were upregulated in the HCC tissues. Kaplan–Meier and Cox regression analyses suggested that highly expressed RAB6B was an independent prognostic factor for poor survival in HCC patients. Moreover, we found that RAB6B expression was positively correlated with the infiltration of immune cells in HCC, including some immunosuppressive cells, chemokines, and receptors, meanwhile RAB6B expression was associated with CD8+T cells exhaustion, resulting in an immunosuppressive microenvironment. Additionally, functional enrichment analysis indicated that RAB6B may be involved in ECM remodeling in the TME, and RAB6B expression was positively associated with CAFs infiltration. Furthermore, RAB6B presented a positive association with sensitivity to GDSC and CTRP drugs. RAB6B knockdown inhibited the cell proliferation and promoted apoptosis and sensitivity to cisplatin of HCC cells in vitro. Conclusion: Our study revealed that RAB6B is a potential biomarker for poor prognosis in HCC patients and correlates with the formation of the immunosuppressive microenvironment in HCC.

In-silico probing of AML related RUNX1 cancer-associated missense mutations: Predicted relationships to DNA binding and drug interactions.

The molecular consequences of cancer associated mutations in Acute myeloid leukemia (AML) linked factors are not very well understood. Here, we interrogated the COSMIC database for missense mutations associated with the RUNX1 protein, that is frequently mis-regulated in AML, where we sought to identify recurrently mutated positions at the DNA-interacting interface. Indeed, six of the mutated residues, out of a total 417 residues examined within the DNA binding domain, evidenced reduced DNA association in in silico predictions. Further, given the prominence of RUNX1’s compromised function in AML, we asked the question if the mutations themselves might alter RUNX1’s interaction (off-target) with known FDA-approved drug molecules, including three currently used in treating AML. We identified several AML-associated mutations in RUNX1 that were calculated to enhance RUNX1’s interaction with specific drugs. Specifically, we retrieved data from the COSMIC database for cancer-associated mutations of RUNX1 by using R package “data.table” and “ggplot2” modules. In the presence of DNA and/or drug, we used docking scores and energetics of the complexes as tools to evaluate predicted interaction strengths with RUNX1. For example, we performed predictions of drug binding pockets involving Enasidenib, Giltertinib, and Midostaurin (AML associated), as well as ten different published cancer associated drug compounds. Docking of wild type RUNX1 with these 13 different cancer-associated drugs indicates that wild-type RUNX1 has a lower efficiency of binding while RUNX1 mutants R142K, D171N, R174Q, P176H, and R177Q suggested higher affinity of drug association. Literature evidence support our prediction and suggests the mutation R174Q affects RUNX1 DNA binding and could lead to compromised function. We conclude that specific RUNX1 mutations that lessen DNA binding facilitate the binding of a number of tested drug molecules. Further, we propose that molecular modeling and docking studies for RUNX1 in the presence of DNA and/or drugs enables evaluation of the potential impact of RUNX1 cancer associated mutations in AML.

In silico mutational analysis to identify the role and pathogenicity of BCL-w missense variants

BackgroundIntrinsic pathway of apoptosis is generally mediated by BCL-2 (B cell lymphoma 2) family of proteins; they either induce or inhibit the apoptosis. Overexpression of BCL-2 in cancer cell may lead to delay in apoptosis. BCL-w is the pro-survival member of the BCL-2 family. BCL2L2 gene is present on chromosome number 14 in humans, and it encodes BCL-w protein; BCL-w protein is 193 amino acids residues in length. Interactions among the BCL-2 proteins are very specific. The fate of cell is determined by the ratio of pro-apoptotic proteins to pro-survival proteins. BCL-w promotes cell survival. Studies suggested that overexpression of BCL-w protein is associated with many cancers including DLBCL, BL, colorectal cancers, gastric cancers, and many more. The cause of overexpression is translocations or gene amplification which will subsequently result in cancerous activity. ProcessFor in-silico analysis, BCL2L2 gene was retrieved from UniProt (UniProt ID: Q92843). 54 missense variants have been collected in BCL-w proteins from COSMIC database. Different tools were used to detect the deleteriousness of the variants. ResultIn silico mutational study reveals how the non-synonymous mutations directly affect the protein’s native structure and its function. Variant mutational analysis with PolyPhen-2 revealed that out of 55 variants, 28 of the missense mutations was probably damaging with a score ranging from 0.9 to 1, while 24 variants were benign with a score ranging from 0 to 0.4. ConclusionsThis in silico work aims to determine how missense mutations in BCL-w protein affect the activity of the protein, the stability of the protein, and to determine the pathogenicity of the variants. Prediction of pathogenicity of variants will reveal if the missense mutation has a damaging effect on the native structure of protein or not. Prediction of protein stability will reveal whether the mutation has a stabilizing or destabilizing effect on the protein.

The Genomic Landscape in Philadelphia-Negative Myeloproliferative Neoplasm Patients with Second Cancers.

Simple SummaryPhiladelphia-chromosome-negative myeloproliferative neoplasms (MPNs) are a group of clonal hematopoietic stem cell disorders characterized by the excessive production of differentiated myeloid cells. Other than the well-known propensity for leukemia transformation, MPN patients are also prone to developing second cancers (SCs) that arise from different embryonic dermal origins. This brings up an intriguing question: what is the molecular background leading to the development of SCs in MPN patients? To explore further, we used whole exome sequencing to characterize the genomic landscapes in 26 paired MPN samples stratified by the presence or absence of SCs. We found that mutated genes in MPN–SC samples were enriched in some critical immune-related pathways and inflammatory networks, an observation further supported by the increased plasma levels of cytokines. Our work provides the genomic landscape that profiles the genetic basis for SC tumorigenesis in MPN patients, and also demonstrates that inflammation could be indispensable in MPN–SC pathogenesis.Patients with myeloproliferative neoplasms (MPNs) are characterized by systemic inflammation. With the indolent nature of the diseases, second cancers (SCs) have emerged as a challenging issue in afflicted patients. Epidemiological studies have confirmed the excessive risk of SCs in MPNs, but little is known about their molecular basis. To explore further, we used whole exome sequencing to explore the genetic changes in the granulocytes of 26 paired MPN patients with or without SC. We noticed that MPN–SC patients harbor genomic variants of distinct genes, among which a unique pattern of co-occurrence or mutual exclusiveness could be identified. We also found that mutated genes in MPN–SC samples were enriched in immune-related pathways and inflammatory networks, an observation further supported by their increased plasma levels of TGF-β and IL-23. Noteworthily, variants of KRT6A, a gene capable of mediating tumor-associate macrophage activity, were more commonly detected in MPN–SC patients. Analysis through OncodriveCLUST disclosed that KRT6A replaces JAK2V617F as the more prominent disease driver in MPN–SC, whereas a major mutation in this gene (KRT6A c.745T>C) in our patients is linked to human carcinoma and predicted to be pathogenic in COSMIC database. Overall, we demonstrate that inflammation could be indispensable in MPN–SC pathogenesis.