Systematic analysis of noncoding somatic mutations and gene expression alterations across 14 tumor types. Nat Genet. 2014;46:1258-1263. [PMID: 25383969 DOI: 10.1038/ng.3141]

Benefits for children with suspected cancer from routine whole-genome sequencing

Clinical whole-genome sequencing (WGS) has been shown to deliver potential benefits to children with cancer and to alter treatment in high-risk patient groups. It remains unknown whether offering WGS to every child with suspected cancer can change patient management. We collected WGS variant calls and clinical and diagnostic information from 281 children (282 tumors) across two English units (n = 152 from a hematology center, n = 130 from a solid tumor center) where WGS had become a routine test. Our key finding was that variants uniquely attributable to WGS changed the management in ~7% (20 out of 282) of cases while providing additional disease-relevant findings, beyond standard-of-care molecular tests, in 108 instances for 83 (29%) cases. Furthermore, WGS faithfully reproduced every standard-of-care molecular test (n = 738) and revealed several previously unknown genomic features of childhood tumors. We show that WGS can be delivered as part of routine clinical care to children with suspected cancer and can change clinical management by delivering unexpected genomic insights. Our experience portrays WGS as a clinically impactful assay for routine practice, providing opportunities for assay consolidation and for delivery of molecularly informed patient care.

Hidden secrets of the cancer genome: unlocking the impact of non-coding mutations in gene regulatory elements

Discoveries in the field of genomics have revealed that non-coding genomic regions are not merely "junk DNA", but rather comprise critical elements involved in gene expression. These gene regulatory elements (GREs) include enhancers, insulators, silencers, and gene promoters. Notably, new evidence shows how mutations within these regions substantially influence gene expression programs, especially in the context of cancer. Advances in high-throughput sequencing technologies have accelerated the identification of somatic and germline single nucleotide mutations in non-coding genomic regions. This review provides an overview of somatic and germline non-coding single nucleotide alterations affecting transcription factor binding sites in GREs, specifically involved in cancer biology. It also summarizes the technologies available for exploring GREs and the challenges associated with studying and characterizing non-coding single nucleotide mutations. Understanding the role of GRE alterations in cancer is essential for improving diagnostic and prognostic capabilities in the precision medicine era, leading to enhanced patient-centered clinical outcomes.

Pan-cancer experimental characteristic of human transcriptional patterns connected with telomerase reverse transcriptase (TERT) gene expression status

The TERT gene encodes the reverse transcriptase subunit of telomerase and is normally transcriptionally suppressed in differentiated human cells but reactivated in cancers where its expression is frequently associated with poor survival prognosis. Here we experimentally assessed the RNA sequencing expression patterns associated with TERT transcription in 1039 human cancer samples of 27 tumor types. We observed a bimodal distribution of TERT expression where ∼27% of cancer samples did not express TERT and the rest showed a bell-shaped distribution. Expression of TERT strongly correlated with 1443 human genes including 103 encoding transcriptional factor proteins. Comparison of TERT- positive and negative cancers showed the differential activation of 496 genes and 1975 molecular pathways. Therein, 32/38 (84%) of DNA repair pathways were hyperactivated in TERT+ cancers which was also connected with accelerated replication, transcription, translation, and cell cycle progression. In contrast, the level of 40 positive cell cycle regulator proteins and a set of epithelial-to-mesenchymal transition pathways was specific for the TERT- group suggesting different proliferation strategies for both groups of cancer. Our pilot study showed that the TERT+ group had ∼13% of cancers with C228T or C250T mutated TERT promoter. However, the presence of promoter mutations was not associated with greater TERT expression compared with other TERT+ cancers, suggesting parallel mechanisms of its transcriptional activation in cancers. In addition, we detected a decreased expression of L1 retrotransposons in the TERT+ group, and further decreased L1 expression in promoter mutated TERT+ cancers. TERT expression was correlated with 17 genes encoding molecular targets of cancer therapeutics and may relate to differential survival patterns of TERT- positive and negative cancers.

A comparative study of structural variant calling in WGS from Alzheimer's disease families

Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.

Stroma-infiltrating T cell spatiotypes define immunotherapy outcomes in adolescent and young adult patients with melanoma

The biological underpinnings of therapeutic resistance to immune checkpoint inhibitors (ICI) in adolescent and young adult (AYA) melanoma patients are incompletely understood. Here, we characterize the immunogenomic profile and spatial architecture of the tumor microenvironment (TME) in AYA (aged ≤ 30 years) and older adult (aged 31-84 years) patients with melanoma, to determine the AYA-specific features associated with ICI treatment outcomes. We identify two ICI-resistant spatiotypes in AYA patients with melanoma showing stroma-infiltrating lymphocytes (SILs) that are distinct from the adult TME. The SILhigh subtype was enriched in regulatory T cells in the peritumoral space and showed upregulated expression of immune checkpoint molecules, while the SILlow subtype showed a lack of immune activation. We establish a young immunosuppressive melanoma score that can predict ICI responsiveness in AYA patients and propose personalized therapeutic strategies for the ICI-resistant subgroups. These findings highlight the distinct immunogenomic profile of AYA patients, and individualized TME features in ICI-resistant AYA melanoma that require patient-specific treatment strategies.

Diphthamide - a conserved modification of eEF2 with clinical relevance

Diphthamide, a complex modification on eukaryotic translation elongation factor 2 (eEF2), assures reading-frame fidelity during translation. Diphthamide and enzymes for its synthesis are conserved in eukaryotes and archaea. Originally identified as target for diphtheria toxin (DT) in humans, its clinical relevance now proves to be broader than the link to pathogenic bacteria. Diphthamide synthesis enzymes (DPH1 and DPH3) are associated with cancer, and DPH gene mutations can cause diphthamide deficiency syndrome (DDS). Finally, new analyses provide evidence that diphthamide may restrict propagation of viruses including SARS-CoV-2 and HIV-1, and that DPH enzymes are targeted by viruses for degradation to overcome this restriction. This review describes how diphthamide is synthesized and functions in translation, and covers its clinical relevance in human development, cancer, and infectious diseases.

Genomic alterations in thyroid cancer: biological and clinical insights

Tumours can arise from thyroid follicular cells if they acquire driver mutations that constitutively activate the MAPK signalling pathway. In addition, a limited set of additional mutations in key genes drive tumour progression towards more aggressive and less differentiated disease. Unprecedented insights into thyroid tumour biology have come from the breadth of thyroid tumour sequencing data from patients and the wide range of mutation-specific mechanisms identified in experimental models, in combination with the genomic simplicity of thyroid cancers. This knowledge is gradually being translated into refined strategies to stratify, manage and treat patients with thyroid cancer. This Review summarizes the biological underpinnings of the genetic alterations involved in thyroid cancer initiation and progression. We also provide a rationale for and discuss specific examples of how to implement genomic information to inform both recommended and investigational approaches to improve thyroid cancer prognosis, redifferentiation strategies and targeted therapies.

Targeting oncogenic TERT promoter variants by allele-specific epigenome editing

BACKGROUND

Activation of dominant oncogenes by small or structural genomic alterations is a common driver mechanism in many cancers. Silencing of such dominantly activated oncogenic alleles, thus, is a promising strategy to treat cancer. Recently, allele-specific epigenome editing (ASEE) has been described as a means to reduce transcription of genes in an allele-specific manner. In cancer, specificity to an oncogenic allele can be reached by either targeting directly a pathogenic single-nucleotide variant or a polymorphic single-nucleotide variant linked to the oncogenic allele. To investigate the potential of ASEE in cancer, we here explored this approach by targeting variants at the TERT promoter region. The TERT promoter region has been described as one of the most frequently mutated non-coding cancer drivers.

RESULTS

Sequencing of the TERT promoter in cancer cell lines showed 53% (41/77) to contain at least one heterozygous sequence variant allowing allele distinction. We chose the hepatoblastoma cell line Hep-G2 and the lung cancer cell line A-549 for this proof-of-principle study, as they contained two different kinds of variants, namely the activating mutation C228T in the TERT core promoter and the common SNP rs2853669 in the THOR region, respectively. These variants were targeted in an allele-specific manner using sgRNA-guided dCas9-DNMT3A-3L complexes. In both cell lines, we successfully introduced DNA methylation specifically to the on-target allele of the TERT promoter with limited background methylation on the off-target allele or an off-target locus (VEGFA), respectively. We observed a maximum CpG methylation gain of 39% and 76% on the target allele when targeting the activating mutation and the common SNP, respectively. The epigenome editing translated into reduced TERT RNA expression in Hep-G2.

CONCLUSIONS

We applied an ASEE-mediated approach to silence TERT allele specifically. Our results show that the concept of dominant oncogene inactivation by allele-specific epigenome editing can be successfully translated into cancer models. This new strategy may have important advantages in comparison with existing therapeutic approaches, e.g., targeting telomerase, especially with regard to reducing adverse side effects.

Reactivation of telomerase reverse transcriptase expression in cancer: the role of TERT promoter mutations

Telomerase activity and telomere elongation are essential conditions for the unlimited proliferation of neoplastic cells. Point mutations in the core promoter region of the telomerase reverse transcriptase (TERT) gene have been found to occur at high frequencies in several tumour types and considered a primary cause of telomerase reactivation in cancer cells. These mutations promote TERT gene expression by multiple mechanisms, including the generation of novel binding sites for nuclear transcription factors, displacement of negative regulators from DNA G-quadruplexes, recruitment of epigenetic activators and disruption of long-range interactions between TERT locus and telomeres. Furthermore, TERT promoter mutations cooperate with TPP1 promoter nucleotide changes to lengthen telomeres and with mutated BRAF and FGFR3 oncoproteins to enhance oncogenic signalling in cancer cells. TERT promoter mutations have been recognized as an early marker of tumour development or a major indicator of poor outcome and reduced patients survival in several cancer types. In this review, we summarize recent findings on the role of TERT promoter mutations, telomerase expression and telomeres elongation in cancer development, their clinical significance and therapeutic opportunities.

Telomerase Upregulation Induces Progression of Mouse BrafV600E-Driven Thyroid Cancers and Triggers Nontelomeric Effects

Mutations in the promoter of the telomerase reverse transcriptase (TERT) gene are the paradigm of a cross-cancer alteration in a noncoding region. TERT promoter mutations (TPM) are biomarkers of poor prognosis in cancer, including thyroid tumors. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert-123C>T) and crossed it with thyroid-specific BrafV600E-mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all BrafV600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of BrafV600E+Tert-123C>T and BrafV600E+K5-Tert mice progressed to poorly differentiated cancers at week 20, respectively. Tert-upregulated tumors showed increased mitosis and necrosis in areas of solid growth, and older animals displayed anaplastic-like features, that is, spindle cells and macrophage infiltration. Murine TPM increased Tert transcription in vitro and in vivo, but temporal and intratumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine, and chemokine signaling, were overactivated. These models constitute useful preclinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs.

IMPLICATIONS

Telomerase-driven cancer progression activates pathways that can be dissected and perhaps therapeutically exploited.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

Identifying functional regulatory mutation blocks by integrating genome sequencing and transcriptome data

Millions of single nucleotide variants (SNVs) exist in the human genome; however, it remains challenging to identify functional SNVs associated with diseases. We propose a non-encoding SNVs analysis tool bpb3, BayesPI-BAR version 3, aiming to identify the functional mutation blocks (FMBs) by integrating genome sequencing and transcriptome data. The identified FMBs display high frequency SNVs, significant changes in transcription factors (TFs) binding affinity and are nearby the regulatory regions of differentially expressed genes. A two-level Bayesian approach with a biophysical model for protein-DNA interactions is implemented, to compute TF-DNA binding affinity changes based on clustered position weight matrices (PWMs) from over 1700 TF-motifs. The epigenetic data, such as the DNA methylome can also be integrated to scan FMBs. By testing the datasets from follicular lymphoma and melanoma, bpb3 automatically and robustly identifies FMBs, demonstrating that bpb3 can provide insight into patho-mechanisms, and therapeutic targets from transcriptomic and genomic data.

PLEKHS1 drives PI3Ks and remodels pathway homeostasis in PTEN-null prostate

The PIP3/PI3K network is a central regulator of metabolism and is frequently activated in cancer, commonly by loss of the PIP3/PI(3,4)P2 phosphatase, PTEN. Despite huge research investment, the drivers of the PI3K network in normal tissues and how they adapt to overactivation are unclear. We find that in healthy mouse prostate PI3K activity is driven by RTK/IRS signaling and constrained by pathway feedback. In the absence of PTEN, the network is dramatically remodeled. A poorly understood YXXM- and PIP3/PI(3,4)P2-binding PH domain-containing adaptor, PLEKHS1, became the dominant activator and was required to sustain PIP3, AKT phosphorylation, and growth in PTEN-null prostate. This was because PLEKHS1 evaded pathway-feedback and experienced enhanced PI3K- and Src-family kinase-dependent phosphorylation of Y258XXM, eliciting PI3K activation. hPLEKHS1 mRNA and activating Y419 phosphorylation of hSrc correlated with PI3K pathway activity in human prostate cancers. We propose that in PTEN-null cells receptor-independent, Src-dependent tyrosine phosphorylation of PLEKHS1 creates positive feedback that escapes homeostasis, drives PIP3 signaling, and supports tumor progression.

Mechanistic Insights of Thyroid Cancer Progression

Differentiated thyroid cancers (DTCs) are primarily initiated by mutations that activate the MAPK signaling cascade, typically at BRAF or RAS oncoproteins. DTCs can evolve to more aggressive forms, specifically, poorly differentiated (PDTC) and anaplastic thyroid cancers (ATC), by acquiring additional genetic alterations which deregulate key pathways. In this review, we focused on bona fide mutations involved in thyroid cancer progression for which consistent mechanistic data exist. Here we summarized the relevant literature, spanning approximately 2 decades, highlighting genetic alterations that are unquestionably enriched in PDTC/ATC. We describe the relevant functional data obtained in multiple in vitro and in vivo thyroid cancer models employed to study genetic alterations in the following genes and functional groups: TP53, effectors of the PI3K/AKT pathway, TERT promoter, members of the SWI/SNF chromatin remodeling complex, NF2, and EIF1AX. In addition, we briefly discuss other genetic alterations that are selected in aggressive thyroid tumors but for which mechanistic data is still either limited or nonexistent. Overall, we argue for the importance conveyed by preclinical studies for the clinical translation of genomic knowledge of thyroid cancers.

Performance Comparison of Computational Methods for the Prediction of the Function and Pathogenicity of Non-coding Variants

Non-coding variants in the human genome significantly influence human traits and complex diseases via their regulation and modification effects. Hence, an increasing number of computational methods are developed to predict the effects of variants in human non-coding sequences. However, it is difficult for inexperienced users to select appropriate computational methods from dozens of available methods. To solve this issue, we assessed 12 performance metrics of 24 methods on four independent non-coding variant benchmark datasets: (1) rare germline variants from clinical relevant sequence variants (ClinVar), (2) rare somatic variants from Catalogue Of Somatic Mutations In Cancer (COSMIC), (3) common regulatory variants from curated expression quantitative trait locus (eQTL) data, and (4) disease-associated common variants from curated genome-wide association studies (GWAS). All 24 tested methods performed differently under various conditions, indicating varying strengths and weaknesses under different scenarios. Importantly, the performance of existing methods was acceptable for rare germline variants from ClinVar with the area under the receiver operating characteristic curve (AUROC) of 0.4481-0.8033 and poor for rare somatic variants from COSMIC (AUROC = 0.4984-0.7131), common regulatory variants from curated eQTL data (AUROC = 0.4837-0.6472), and disease-associated common variants from curated GWAS (AUROC = 0.4766-0.5188). We also compared the prediction performance of 24 methods for non-coding de novo mutations in autism spectrum disorder, and found that the combined annotation-dependent depletion (CADD) and context-dependent tolerance score (CDTS) methods showed better performance. Summarily, we assessed the performance of 24 computational methods under diverse scenarios, providing preliminary advice for proper tool selection and guiding the development of new techniques in interpreting non-coding variants.

InTERTwined: how TERT promoter mutations impact BRAFV600E-driven thyroid cancers

Thyroid cancers are often initiated by the acquisition of a BRAFV600E mutation. BRAFV600E-driven thyroid tumors display a wide range of behaviors, from the slow-growing papillary carcinomas to the highly aggressive anaplastic. Mutations in the promoter of TERT (telomerase reverse transcriptase) gene were discovered a decade ago and identified as prevalent events in thyroid cancers. Multiple studies showed that TERT promoter mutations, particularly when co-occurring with BRAFV600E, are markers of poor prognosis across thyroid cancer subtypes, and can be implemented for routine clinical stratification. Mechanistically, TERT promoter mutations reactivate telomerase expression via the differential recruitment of transcriptional complexes. Re-expression of TERT impacts tumor biology, plausibly via both the well-known function of telomerase maintaining telomeres and by affecting other cancer-relevant processes.

Base-resolution UV footprinting by sequencing reveals distinctive damage signatures for DNA-binding proteins

Decades ago, it was shown that proteins binding to DNA can quantitatively alter the formation of DNA damage by UV light. This established the principle of UV footprinting for non-intrusive study of protein-DNA contacts in living cells, albeit at limited scale and precision. Here, we perform deep base-resolution quantification of the principal UV damage lesion, the cyclobutane pyrimidine dimer (CPD), at select human promoter regions using targeted CPD sequencing. Several transcription factors exhibited distinctive and repeatable damage signatures indicative of site occupancy, involving strong (up to 17-fold) position-specific elevations and reductions in CPD formation frequency relative to naked DNA. Positive damage modulation at some ETS transcription factor binding sites coincided at base level with melanoma somatic mutation hotspots. Our work provides proof of concept for the study of protein-DNA interactions at individual loci using light and sequencing, and reveals widespread and potent modulation of UV damage in regulatory regions.

Detecting recurrent passenger mutations in melanoma by targeted UV damage sequencing

Sequencing of melanomas has identified hundreds of recurrent mutations in both coding and non-coding DNA. These include a number of well-characterized oncogenic driver mutations, such as coding mutations in the BRAF and NRAS oncogenes, and non-coding mutations in the promoter of telomerase reverse transcriptase (TERT). However, the molecular etiology and significance of most of these mutations is unknown. Here, we use a new method known as CPD-capture-seq to map UV-induced cyclobutane pyrimidine dimers (CPDs) with high sequencing depth and single nucleotide resolution at sites of recurrent mutations in melanoma. Our data reveal that many previously identified drivers and other recurrent mutations in melanoma occur at CPD hotspots in UV-irradiated melanocytes, often associated with an overlapping binding site of an E26 transformation-specific (ETS) transcription factor. In contrast, recurrent mutations in the promoters of a number of known or suspected cancer genes are not associated with elevated CPD levels. Our data indicate that a subset of recurrent protein-coding mutations are also likely caused by ETS-induced CPD hotspots. This analysis indicates that ETS proteins profoundly shape the mutation landscape of melanoma and reveals a method for distinguishing potential driver mutations from passenger mutations whose recurrence is due to elevated UV damage.

Evaluation of Breast Implant-Associated Anaplastic Large Cell Lymphoma With Whole Exome and Genome Sequencing

BACKGROUND

Breast implant-associated anaplastic large cell lymphoma (BIA-ALCL) is a rare malignancy originating from the periprosthetic capsule of a textured, most often macrotextured, breast implant. Identified in women whose indications for breast implants can be either aesthetic or reconstructive, the genomic underpinnings of this disease are only beginning to be elucidated.

OBJECTIVES

The aim of this study was to evaluate the exomes, and in some cases the entire genome, of patients with BIA-ALCL. Specific attention was paid to copy number alterations, chromosomal translocations, and other genomic abnormalities overrepresented in patients with BIA-ALCL.

METHODS

Whole-exome sequencing was performed on 6 patients, and whole-genome sequencing on 3 patients, with the Illumina NovaSeq 6000 sequencer. Data were analyzed with the Illumina DRAGEN Bio-IT Platform and the ChromoSeq pipeline. The Pathseq Genome Analysis Toolkit pipeline was used to detect the presence of microbial genomes in the sequenced samples.

RESULTS

Two cases with STAT3 mutations and 2 cases with NRAS mutations were noted. A critically deleted 7-Mb region was identified at the 11q22.3 region of chromosome 11, and multiple nonrecurrent chromosomal rearrangements were identified by whole-genome sequencing. Recurrent gene-level rearrangements, however, were not identified. None of the samples showed evidence of potential microbial pathogens.

CONCLUSIONS

Although no recurrent mutations were identified, this study identified mutations in genes not previously reported with BIA-ALCL or other forms of ALCL. Furthermore, not previously reported with BIA-ALCL, 11q22.3 deletions were consistent across whole-genome sequencing cases and present in some exomes.

LEVEL OF EVIDENCE: 5

Telomerase reactivation induces progression of mouse Braf V600E -driven thyroid cancers without telomere lengthening

Mutations in the promoter of the telomerase reverse transcriptase ( TERT ) gene are the paradigm of a cross-cancer alteration in a non-coding region. TERT promoter mutations (TPMs) are biomarkers of poor prognosis in several tumors, including thyroid cancers. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert -123C>T ) and crossed it with thyroid-specific Braf V600E -mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all Braf V600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of Braf V600E +Tert -123C>T and Braf V600E +K5-Tert mice progressed to poorly differentiated thyroid cancers at week 20, respectively. Braf+Tert tumors showed increased mitosis and necrosis in areas of solid growth, and older animals from these cohorts displayed anaplastic-like features, i.e., spindle cells and macrophage infiltration. Murine Tert promoter mutation increased Tert transcription in vitro and in vivo , but temporal and intra-tumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine and chemokine signaling, were overactivated. Braf+Tert animals remained responsive to MAPK pathway inhibitors. These models constitute useful pre-clinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs.

Genomic and epigenomic profile of thyroid cancer

Thyroid cancer is the most common malignancy of the endocrine system, and its incidence has been steadily increasing. Advances in sequencing have allowed analysis of the entire cancer genome, and has provided new information on the genetic lesions and modifications responsible for the onset, progression, dedifferentiation and metastasis of thyroid carcinomas. Moreover, integrated genomics has advanced our understanding of the development of cancer and its behavior, and has facilitated the identification of new genetic mutations and molecular pathways. The functional analysis of epigenetic modifications, such as DNA methylation, histone acetylation and non-coding RNAs, have contributed to define new regulatory mechanisms that control cell malignancy in thyroid cancer, especially aggressive forms. Here we review the most recent advances in genomics and epigenomics of thyroid cancer, which have resulted in a new classification and interpretation of the initiation and progression of thyroid tumors, providing new tools and opportunities for further investigation and for the clinical development of new treatment strategies.

Cis-regulatory mutations associate with transcriptional and post-transcriptional deregulation of gene regulatory programs in cancers

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

Integrated Analysis of Altered lncRNA, circRNA, microRNA, and mRNA Expression in Hepatocellular Carcinoma Carrying TERT Promoter Mutations

BACKGROUND AND OBJECTIVES

Telomerase reverse transcriptase (TERT) promoter mutations are one of the most common mutations responsible for the development of hepatocellular carcinoma (HCC). Noncoding RNAs (ncRNAs) play a regulatory role in different cancers through the long noncoding RNA (lncRNA)/circular RNA (circRNA)-microRNA (miRNA)-mRNA axis. The aim of the study was to explore the influence of TERT promoter mutations on the ncRNA regulatory network in HCC.

METHODS

Four tumor samples with a wildtype TERT promoter and four tumor samples with TERT promoter mutations (sequencing cohort) were collected from HCC patients for high-throughput next-generation sequencing. Selected ncRNAs and mRNAs were validated by qPCR in 15 HCC tumors with a wildtype TERT promoter and seven HCC tumors with TERT promoter mutations (validation cohort, including the sequencing cohort).

RESULTS

In the mutant TERT promoter group, 536 lncRNAs, 21 circRNAs, 41 miRNAs, and 266 mRNAs were significantly up-regulated, while 1745 lncRNAs, 23 circRNAs, 32 miRNAs, and 1117 mRNAs were significantly down-regulated (P < 0.05) compared with the findings in wildtype group. AL360169.3-201, LINC02672-203, hsa_circ_0021412, hsa-miR-29b-1-5p, hsa-miR-4699-5p, hsa-miR-199a-5p, REG3A, SFRP5, and GSTM1 were verified at the RNA expression level to validate the sequencing results. A differentially expressed lncRNA/circRNA-miRNA-mRNA network was constructed to explore the effects of TERT promoter mutations on ncRNA regulation. Two ncRNA regulatory axes associated with TERT promoter mutations (hsa_circ_0003154/hsa_circ_0008952/IGLL5-AS1/LINC576/LINC575-hsa-miR-1260b -CLPTM1L/GSTM1 and hsa_circ_0031584/LINC2101-hsa-miR-214-3p-CD151) had carcinogenic potential.

CONCLUSION

This study provides novel insights into the role of TERT promoter mutations on ncRNAs regulatory network in HCC progression.

Somatic mutation distribution across tumour cohorts provides a signal for positive selection in cancer

Cancer gene discovery is reliant on distinguishing driver mutations from a multitude of passenger mutations in tumour genomes. While driver genes may be revealed based on excess mutation recurrence or clustering, there is a need for orthogonal principles. Here, we take advantage of the fact that non-cancer genes, containing only passenger mutations under neutral selection, exhibit a likelihood of mutagenesis in a given tumour determined by the tumour's mutational signature and burden. This relationship can be disrupted by positive selection, leading to a difference in the distribution of mutated cases across a cohort for driver and passenger genes. We apply this principle to detect cancer drivers independently of recurrence in large pan-cancer cohorts, and show that our method (SEISMIC) performs comparably to traditional approaches and can provide resistance to known confounding mutational phenomena. Being based on a different principle, the approach provides a much-needed complement to existing methods for detecting signals of selection.

Genome-wide mapping of somatic mutation rates uncovers drivers of cancer

Identification of cancer driver mutations that confer a proliferative advantage is central to understanding cancer; however, searches have often been limited to protein-coding sequences and specific non-coding elements (for example, promoters) because of the challenge of modeling the highly variable somatic mutation rates observed across tumor genomes. Here we present Dig, a method to search for driver elements and mutations anywhere in the genome. We use deep neural networks to map cancer-specific mutation rates genome-wide at kilobase-scale resolution. These estimates are then refined to search for evidence of driver mutations under positive selection throughout the genome by comparing observed to expected mutation counts. We mapped mutation rates for 37 cancer types and applied these maps to identify putative drivers within intronic cryptic splice regions, 5' untranslated regions and infrequently mutated genes. Our high-resolution mutation rate maps, available for web-based exploration, are a resource to enable driver discovery genome-wide.

The Role of Telomerase in Breast Cancer's Response to Therapy

Currently, breast cancer appears to be the most widespread cancer in the world and the most common cause of cancer deaths. This specific type of cancer affects women in both developed and developing countries. Prevention and early diagnosis are very important factors for good prognosis. A characteristic feature of cancer cells is the ability of unlimited cell division, which makes them immortal. Telomeres, which are shortened with each cell division in normal cells, are rebuilt in cancer cells by the enzyme telomerase, which is expressed in more than 85% of cancers (up to 100% of adenocarcinomas, including breast cancer). Telomerase may have different functions that are related to telomeres or unrelated. It has been shown that high activity of the enzyme in cancer cells is associated with poor cell sensitivity to therapies. Therefore, telomerase has become a potential target for cancer therapies. The low efficacy of therapies has resulted in the search for new combined and more effective therapeutic methods, including the involvement of telomerase inhibitors and telomerase-targeted immunotherapy.

Multiplexed functional genomic assays to decipher the noncoding genome

Linkage disequilibrium and the incomplete regulatory annotation of the noncoding genome complicates the identification of functional noncoding genetic variants and their causal association with disease. Current computational methods for variant prioritization have limited predictive value, necessitating the application of highly parallelized experimental assays to efficiently identify functional noncoding variation. Here, we summarize two distinct approaches, massively parallel reporter assays and CRISPR-based pooled screens and describe their flexible implementation to characterize human noncoding genetic variation at unprecedented scale. Each approach provides unique advantages and limitations, highlighting the importance of multimodal methodological integration. These multiplexed assays of variant effects are undoubtedly poised to play a key role in the experimental characterization of noncoding genetic risk, informing our understanding of the underlying mechanisms of disease-associated loci and the development of more robust predictive classification algorithms.

Clinical utility of whole-genome sequencing in precision oncology

Precision diagnostics is one of the two pillars of precision medicine. Sequencing efforts in the past decade have firmly established cancer as a primarily genetically driven disease. This concept is supported by therapeutic successes aimed at particular pathways that are perturbed by specific driver mutations in protein-coding domains and reflected in three recent FDA tissue agnostic cancer drug approvals. In addition, there is increasing evidence from studies that interrogate the entire genome by whole-genome sequencing that acquired global and complex genomic aberrations including those in non-coding regions of the genome might also reflect clinical outcome. After addressing technical, logistical, financial and ethical challenges, national initiatives now aim to introduce clinical whole-genome sequencing into real-world diagnostics as a rational and potentially cost-effective tool for response prediction in cancer and to identify patients who would benefit most from 'expensive' targeted therapies and recruitment into clinical trials. However, so far, this has not been accompanied by a systematic and prospective evaluation of the clinical utility of whole-genome sequencing within clinical trials of uniformly treated patients of defined clinical outcome. This approach would also greatly facilitate novel predictive biomarker discovery and validation, ultimately reducing size and duration of clinical trials and cost of drug development. This manuscript is the third in a series of three to review and critically appraise the potential and challenges of clinical whole-genome sequencing in solid tumors and hematological malignancies.

High-throughput techniques enable advances in the roles of DNA and RNA secondary structures in transcriptional and post-transcriptional gene regulation

The most stable structure of DNA is the canonical right-handed double helix termed B DNA. However, certain environments and sequence motifs favor alternative conformations, termed non-canonical secondary structures. The roles of DNA and RNA secondary structures in transcriptional regulation remain incompletely understood. However, advances in high-throughput assays have enabled genome wide characterization of some secondary structures. Here, we describe their regulatory functions in promoters and 3’UTRs, providing insights into key mechanisms through which they regulate gene expression. We discuss their implication in human disease, and how advances in molecular technologies and emerging high-throughput experimental methods could provide additional insights.

Epidemiology of Glioblastoma Multiforme–Literature Review

Simple Summary

Glioblastoma multiforme (GBM) is one of the most aggressive malignancies, accounting for 14.5% of all central nervous system tumors and 48.6% of malignant central nervous system tumors. The median overall survival (OS) of GBM patients is only 15 months. The aim of this review was to provide an overview of the epidemiology of GBM and factors that may have a significant impact on the risk of GBM.

Abstract

Glioblastoma multiforme (GBM) is one of the most aggressive malignancies, with a median overall survival of approximately 15 months. In this review, we analyze the pathogenesis of GBM, as well as epidemiological data, by age, gender, and tumor location. The data indicate that GBM is the higher-grade primary brain tumor and is significantly more common in men. The risk of being diagnosed with glioma increases with age, and median survival remains low, despite medical advances. In addition, it is difficult to determine clearly how GBM is influenced by stimulants, certain medications (e.g., NSAIDs), cell phone use, and exposure to heavy metals.

Micropeptides translated from putative long non-coding RNAs

Long non-coding RNAs (lncRNAs) transcribed in mammals and eukaryotes were thought to have no protein coding capability. However, recent studies have suggested that plenty of lncRNAs are mis-annotated and virtually contain coding sequences which are translated into functional peptides by ribosomal machinery, and these functional peptides are called micropeptides or small peptides. Here we review the rapidly advancing field of micropeptides translated from putative lncRNAs, describe the strategies for their identification, and elucidate their critical roles in many fundamental biological processes. We also discuss the prospects of research in micropeptides and the potential applications of micropeptides.

Core promoter mutation contributes to abnormal gene expression in bladder cancer

Background

Bladder cancer is one of the most mortal cancers. Bladder cancer has distinct gene expression signature, highlighting altered gene expression plays important roles in bladder cancer etiology. However, the mechanism for how the regulatory disorder causes the altered expression in bladder cancer remains elusive. Core promoter controls transcriptional initiation. We hypothesized that mutation in core promoter abnormality could cause abnormal transcriptional initiation thereby the altered gene expression in bladder cancer.

Methods

In this study, we performed a genome-wide characterization of core promoter mutation in 77 Spanish bladder cancer cases.

Results

We identified 69 recurrent somatic mutations in 61 core promoters of 62 genes and 28 recurrent germline mutations in 20 core promoters of 21 genes, including TERT, the only gene known with core promoter mutation in bladder cancer, and many oncogenes and tumor suppressors. From the RNA-seq data from bladder cancer, we observed  altered expression of the core promoter-mutated genes. We further validated the effects of core promoter mutation on gene expression by using luciferase reporter gene assay. We also identified potential drugs targeting the core promoter-mutated genes.

Conclusions

Data from our study highlights that core promoter mutation contributes to bladder cancer development through altering gene expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09178-z.

Mutations in Noncoding Cis-Regulatory Elements Reveal Cancer Driver Cistromes in Luminal Breast Cancer

Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and noncoding cancer driver plexuses from somatic mutations. However, differentiating driver from passenger events among noncoding genetic variants remains a challenge. Herein, we reveal cancer-driver cis-regulatory elements linked to transcription factors previously shown to be involved in development of luminal breast cancers by defining a tumor-enriched catalogue of approximately 100,000 unique cis-regulatory elements from 26 primary luminal estrogen receptor (ER)⁺ progesterone receptor (PR)⁺ breast tumors. Integrating this catalog with somatic mutations from 350 publicly available breast tumor whole genomes, we uncovered cancer driver cistromes, defined as the sum of binding sites for a transcription factor, for ten transcription factors in luminal breast cancer such as FOXA1 and ER, nine of which are essential for growth in breast cancer with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of noncoding mutations over cis-regulatory elements concatenated into a functional unit. IMPLICATIONS: Mapping the accessible chromatin of luminal breast cancer led to discovery of an accumulation of mutations within cistromes of transcription factors essential to luminal breast cancer. This demonstrates coopting of regulatory networks to drive cancer and provides a framework to derive insight into the noncoding space of cancer.

Spectrum of mitochondrial genomic variation in parathyroid neoplasms

PURPOSE

Mitochondrial DNA (mtDNA) variations have been implicated in various cancer types. Several attempts have been made in benign parathyroid adenoma (PA); however, mtDNA variants and copy number alterations have not been investigated in parathyroid carcinoma (PC). The present study aimed to explore the variations in the mitochondrial genome in parathyroid neoplasms.

METHODS

In this study, ultradeep targeted sequencing of mtDNA was performed in 12 cases with PC and 25 cases with PA. Germline and somatic variants in the mitochondrial genome were analysed according to clinical features. The mtDNA copy number relative to nuclear DNA was measured.

RESULTS

A total of 821 germline variants and 23 somatic mutations were identified. All 160 germline nonsynonymous variants in the coding sequence (CDS) were predicted to be likely benign, and germline variants frequently occurred (26.3%) in the displacement loop region. Tumour somatic mutation in a CDS with heteroplasmic factor (HF) >0.8 showed a higher mtDNA copy number (p = 0.039). Using Spearman's rank test, tumour mtDNA copy number revealed a positive correlation with serum levels of intact parathyroid hormone and calcium.

CONCLUSION

Our results demonstrate that null recurrent mtDNA mutational hotspots were PC specific. An increased mtDNA copy number in parathyroid neoplasms was associated with somatic CDS mutations with a high HF and major clinical features. Larger cohort investigations should be performed in future analyses.

ISOTOPE: ISOform-guided prediction of epiTOPEs in cancer

Immunotherapies provide effective treatments for previously untreatable tumors and identifying tumor-specific epitopes can help elucidate the molecular determinants of therapy response. Here, we describe a pipeline, ISOTOPE (ISOform-guided prediction of epiTOPEs In Cancer), for the comprehensive identification of tumor-specific splicing-derived epitopes. Using RNA sequencing and mass spectrometry for MHC-I associated proteins, ISOTOPE identified neoepitopes from tumor-specific splicing events that are potentially presented by MHC-I complexes. Analysis of multiple samples indicates that splicing alterations may affect the production of self-epitopes and generate more candidate neoepitopes than somatic mutations. Although there was no difference in the number of splicing-derived neoepitopes between responders and non-responders to immune therapy, higher MHC-I binding affinity was associated with a positive response. Our analyses highlight the diversity of the immunogenic impacts of tumor-specific splicing alterations and the importance of studying splicing alterations to fully characterize tumors in the context of immunotherapies. ISOTOPE is available at https://github.com/comprna/ISOTOPE.

Naturally occurring Telomerase-specific CD4 T cell Immunity in Melanoma

CD4 T cells play a key role in anticancer immunity. Here, we investigate the clinical relevance of circulating CD4 Th1 response against telomerase (anti-TERT Th1 response) in melanoma patients. The spontaneous anti-TERT Th1 response was detected in 54.5% (85/156) of melanoma patients before treatment. The prevalence of this systemic response was inversely related to Breslow thickness above 1mm and AJCC stage ≥ II (P = 0.001 and 0.032). In contrast to patients treated by targeted therapies, the anti-TERT Th1 immunity was associated with objective response after immune checkpoint inhibitors (ICI) treatment. Hence 86% (18/21) of responder patients exhibited pre-existing anti-TERT Th1 versus 35% (6/19) in non-responders (P = 0.001). This response was also associated with increased progression free survival and overall survival in melanoma patients treated with ICI (P = 0.0008 and 0.012 respectively). Collectively, the presence of circulating anti-TERT Th1 response is inversely related to melanoma evolution and appears to be a predictive factor of response to immunotherapy. Our results highlight the interest of telomerase-specific CD4 Th1 response as a promising blood based biomarker of immune checkpoint inhibitors therapy in melanoma.

Non-coding driver mutations in human cancer

Tumour formation involves random mutagenic events and positive evolutionary selection acting on a subset of such events, referred to as driver mutations. A decade of careful surveying of tumour DNA using exome-based analyses has revealed a multitude of protein-coding somatic driver mutations, some of which are clinically actionable. Today, a transition towards whole-genome analysis is well under way, technically enabling the discovery of potential driver mutations occurring outside protein-coding sequences. Mutations are abundant in this vast non-coding space, which is more than 50 times larger than the coding exome, but reliable identification of selection signals in non-coding DNA remains a challenge. In this Review, we discuss recent findings in the field, where the emerging landscape is one in which non-coding driver mutations appear to be relatively infrequent. Nevertheless, we highlight several notable discoveries. We consider possible reasons for the relative absence of non-coding driver events, as well as the difficulties associated with detecting signals of positive selection in non-coding DNA.

Occurrence, functionality and abundance of the TERT promoter mutations

Telomere shortening at chromosomal ends due to the constraints of the DNA replication process acts as a tumor suppressor by restricting the replicative potential in primary cells. Cancers evade that limitation primarily through the reactivation of telomerase via different mechanisms. Mutations within the promoter of the telomerase reverse transcriptase (TERT) gene represent a definite mechanism for the ribonucleic enzyme regeneration predominantly in cancers that arise from tissues with low rates of self-renewal. The promoter mutations cause a moderate increase in TERT transcription and consequent telomerase upregulation to the levels sufficient to delay replicative senescence but not prevent bulk telomere shortening and genomic instability. Since the discovery, a staggering number of studies have resolved the discrete aspects, effects and clinical relevance of the TERT promoter mutations. The promoter mutations link transcription of TERT with oncogenic pathways, associate with markers of poor outcome and define patients with reduced survivals in several cancers. In this review, we discuss the occurrence and impact of the promoter mutations and highlight the mechanism of TERT activation. We further deliberate on the foundational question of the abundance of the TERT promoter mutations and a general dearth of functional mutations within noncoding sequences, as evident from pan-cancer analysis of the whole-genomes. We posit that the favorable genomic constellation within the TERT promoter may be less than a common occurrence in other noncoding functional elements. Besides, the evolutionary constraints limit the functional fraction within the human genome, hence the lack of abundant mutations outside the coding sequences.

Functional annotation of noncoding mutations in cancer

Recurrent regulatory mutations affecting transcription factor binding sites in 2,500 cancer samples.

In a cancer genome, the noncoding sequence contains the vast majority of somatic mutations. While very few are expected to be cancer drivers, those affecting regulatory elements have the potential to have downstream effects on gene regulation that may contribute to cancer progression. To prioritize regulatory mutations, we screened somatic mutations in the Pan-Cancer Analysis of Whole Genomes cohort of 2,515 cancer genomes on individual bases to assess their potential regulatory roles in their respective cancer types. We found a highly significant enrichment of regulatory mutations associated with the deamination signature overlapping a CpG site in the CCAAT/Enhancer Binding Protein β recognition sites in many cancer types. Overall, 5,749 mutated regulatory elements were identified in 1,844 tumor samples from 39 cohorts containing 11,962 candidate regulatory mutations. Our analysis indicated 20 or more regulatory mutations in 5.5% of the samples, and an overall average of six per tumor. Several recurrent elements were identified, and major cancer-related pathways were significantly enriched for genes nearby the mutated regulatory elements. Our results provide a detailed view of the role of regulatory elements in cancer genomes.

Analysis of gene expression and mutation data points on contribution of transcription to the mutagenesis by APOBEC enzymes

Since the discovery of the role of the APOBEC enzymes in human cancers, the mechanisms of this type of mutagenesis remain little understood. Theoretically, targeting of single-stranded DNA by the APOBEC enzymes could occur during cellular processes leading to the unwinding of DNA double-stranded structure. Some evidence points to the importance of replication in the APOBEC mutagenesis, while the role of transcription is still underexplored. Here, we analyzed gene expression and whole genome sequencing data from five types of human cancers with substantial APOBEC activity to estimate the involvement of transcription in the APOBEC mutagenesis and compare its impact with that of replication. Using the TCN motif as the mutation signature of the APOBEC enzymes, we observed a correlation of active APOBEC mutagenesis with gene expression, confirmed the increase of APOBEC-induced mutations in early-replicating regions and estimated the relative impact of transcription and replication on the APOBEC mutagenesis. We also found that the known effect of higher density of APOBEC-induced mutations on the lagging strand was highest in middle-replicating regions and observed higher APOBEC mutation density on the sense strand, the latter bias positively correlated with the gene expression level.

Enhanced C/EBP binding to G·T mismatches facilitates fixation of CpG mutations in cancer and adult stem cells

Somatic mutations in regulatory sites of human stem cells affect cell identity or cause malignant transformation. By mining the human genome for co-occurrence of mutations and transcription factor binding sites, we show that C/EBP binding sites are strongly enriched with [C > T]G mutations in cancer and adult stem cells, which is of special interest because C/EBPs regulate cell fate and differentiation. In vitro protein-DNA binding assay and structural modeling of the CEBPB-DNA complex show that the G·T mismatch in the core CG dinucleotide strongly enhances affinity of the binding site. We conclude that enhanced binding of C/EBPs shields CpG·TpG mismatches from DNA repair, leading to selective accumulation of [C > T]G mutations and consequent deterioration of the binding sites. This mechanism of targeted mutagenesis highlights the effect of a mutational process on certain regulatory sites and reveals the molecular basis of putative regulatory alterations in stem cells.

Selective requirement of MYB for oncogenic hyperactivation of a translocated enhancer in leukemia

In acute myeloid leukemia (AML) with inv(3)(q21;q26) or t(3;3)(q21;q26), a translocated GATA2 enhancer drives oncogenic expression of EVI1. We generated an EVI1-GFP AML model and applied an unbiased CRISPR/Cas9 enhancer scan to uncover sequence motifs essential for EVI1 transcription. Using this approach, we pinpointed a single regulatory element in the translocated GATA2 enhancer that is critically required for aberrant EVI1 expression. This element contained a DNA binding motif for the transcription factor MYB which specifically occupied this site at the translocated allele and was dispensable for GATA2 expression. MYB knockout as well as peptidomimetic blockade of CBP/p300-dependent MYB functions resulted in downregulation of EVI1 but not of GATA2. Targeting MYB or mutating its DNA-binding motif within the GATA2 enhancer resulted in myeloid differentiation and cell death, suggesting that interference with MYB-driven EVI1 transcription provides a potential entry point for therapy of inv(3)/t(3;3) AMLs.

Functional and genetic determinants of mutation rate variability in regulatory elements of cancer genomes

Background

Cancer genomes are shaped by mutational processes with complex spatial variation at multiple scales. Entire classes of regulatory elements are affected by local variations in mutation frequency. However, the underlying mechanisms with functional and genetic determinants remain poorly understood.

Results

We characterise the mutational landscape of 1.3 million gene-regulatory and chromatin architectural elements in 2419 whole cancer genomes with transcriptional and pathway activity, functional conservation and recurrent driver events. We develop RM2, a statistical model that quantifies mutational enrichment or depletion in classes of genomic elements through genetic, trinucleotide and megabase-scale effects. We report a map of localised mutational processes affecting CTCF binding sites, transcription start sites (TSS) and tissue-specific open-chromatin regions. Increased mutation frequency in TSSs associates with mRNA abundance in most cancer types, while open-chromatin regions are generally enriched in mutations. We identify ~ 10,000 CTCF binding sites with core DNA motifs and constitutive binding in 66 cell types that represent focal points of mutagenesis. We detect site-specific mutational signature enrichments, such as SBS40 in open-chromatin regions in prostate cancer and SBS17b in CTCF binding sites in gastrointestinal cancers. Candidate drivers of localised mutagenesis are also apparent: BRAF mutations associate with mutational enrichments at CTCF binding sites in melanoma, and ARID1A mutations with TSS-specific mutagenesis in pancreatic cancer.

Conclusions

Our method and catalogue of localised mutational processes provide novel perspectives to cancer genome evolution, mutagenesis, DNA repair and driver gene discovery. The functional and genetic correlates of mutational processes suggest mechanistic hypotheses for future studies.

Cell-Free DNA Detection of Tumor Mutations in Heterogeneous, Localized Prostate Cancer Via Targeted, Multiregion Sequencing

Cell-free DNA (cfDNA) may allow for minimally invasive identification of biologically relevant genomic alterations and genetically distinct tumor subclones. Although existing biomarkers may detect localized prostate cancer, additional strategies interrogating genomic heterogeneity are necessary for identifying and monitoring aggressive disease. In this study, we aimed to evaluate whether circulating tumor DNA can detect genomic alterations present in multiple regions of localized prostate tumor tissue.

METHODS

Low-pass whole-genome and targeted sequencing with a machine-learning guided 2.5-Mb targeted panel were used to identify single nucleotide variants, small insertions and deletions (indels), and copy-number alterations in cfDNA. The majority of this study focuses on the subset of 21 patients with localized disease, although 45 total individuals were evaluated, including 15 healthy controls and nine men with metastatic castration-resistant prostate cancer. Plasma cfDNA was barcoded with duplex unique molecular identifiers. For localized cases, matched tumor tissue was collected from multiple regions (one to nine samples per patient) for comparison.

RESULTS

Somatic tumor variants present in heterogeneous tumor foci from patients with localized disease were detected in cfDNA, and cfDNA mutational burden was found to track with disease severity. Somatic tissue alterations were identified in cfDNA, including nonsynonymous variants in FOXA1, PTEN, MED12, and ATM. Detection of these overlapping variants was associated with seminal vesicle invasion (P = .019) and with the number of variants initially found in the matched tumor tissue samples (P = .0005).

CONCLUSION

Our findings demonstrate the potential of targeted cfDNA sequencing to detect somatic tissue alterations in heterogeneous, localized prostate cancer, especially in a setting where matched tumor tissue may be unavailable (ie, active surveillance or treatment monitoring).

Clonal expansion in non-cancer tissues

Cancer is a clonal disorder derived from a single ancestor cell and its progenies that are positively selected by acquisition of 'driver mutations'. However, the evolution of positively selected clones does not necessarily imply the presence of cancer. On the contrary, it has become clear that expansion of these clones in phenotypically normal or non-cancer tissues is commonly seen in association with ageing and/or in response to environmental insults and chronic inflammation. Recent studies have reported expansion of clones harbouring mutations in cancer driver genes in the blood, skin, oesophagus, bronchus, liver, endometrium and bladder, where the expansion could be so extensive that tissues undergo remodelling of an almost entire tissue. The presence of common cancer driver mutations in normal tissues suggests a strong link to cancer development, providing an opportunity to understand early carcinogenic processes. Nevertheless, some driver mutations are unique to normal tissues or have a mutation frequency that is much higher in normal tissue than in cancer, indicating that the respective clones may not necessarily be destined for evolution to cancer but even negatively selected for carcinogenesis depending on the mutated gene. Moreover, tissues that are remodelled by genetically altered clones might define functionalities of aged tissues or modified inflammatory processes. In this Review, we provide an overview of major findings on clonal expansion in phenotypically normal or non-cancer tissues and discuss their biological significance not only in cancer development but also in ageing and inflammatory diseases.