Preclinical evaluation of CDK4 phosphorylation predicts high sensitivity of pleural mesotheliomas to CDK4/6 inhibition

Malignant pleural mesothelioma (MPM) is an aggressive cancer with limited therapeutic options. We evaluated the impact of CDK4/6 inhibition by palbociclib in 28 MPM cell lines including 19 patient-derived ones, using various approaches including RNA-sequencing. Palbociclib strongly and durably inhibited the proliferation of 23 cell lines, indicating a unique sensitivity of MPM to CDK4/6 inhibition. When observed, insensitivity to palbociclib was mostly explained by the lack of active T172-phosphorylated CDK4. This was associated with high p16INK4A (CDKN2A) levels that accompany RB1 defects or inactivation, or (unexpectedly) CCNE1 overexpression in the presence of wild-type RB1. Prolonged palbociclib treatment irreversibly inhibited proliferation despite re-induction of cell cycle genes upon drug washout. A senescence-associated secretory phenotype including various potentially immunogenic components was irreversibly induced. Phosphorylated CDK4 was detected in 80% of 47 MPMs indicating their sensitivity to CDK4/6 inhibitors. Its absence in some highly proliferative MPMs was linked to very high p16 (CDKN2A) expression, which was also observed in public datasets in tumours from short-survival patients. Our study supports the evaluation of CDK4/6 inhibitors for MPM treatment, in monotherapy or combination therapy.

Genomic evolution shapes prostate cancer disease type

The development of cancer is an evolutionary process involving the sequential acquisition of genetic alterations that disrupt normal biological processes, enabling tumor cells to rapidly proliferate and eventually invade and metastasize to other tissues. We investigated the genomic evolution of prostate cancer through the application of three separate classification methods, each designed to investigate a different aspect of tumor evolution. Integrating the results revealed the existence of two distinct types of prostate cancer that arise from divergent evolutionary trajectories, designated as the Canonical and Alternative evolutionary disease types. We therefore propose the evotype model for prostate cancer evolution wherein Alternative-evotype tumors diverge from those of the Canonical-evotype through the stochastic accumulation of genetic alterations associated with disruptions to androgen receptor DNA binding. Our model unifies many previous molecular observations, providing a powerful new framework to investigate prostate cancer disease progression.

Aneuploidy and complex genomic rearrangements in cancer evolution

Mutational processes that alter large genomic regions occur frequently in developing tumors. They range from simple copy number gains and losses to the shattering and reassembly of entire chromosomes. These catastrophic events, such as chromothripsis, chromoplexy and the formation of extrachromosomal DNA, affect the expression of many genes and therefore have a substantial effect on the fitness of the cells in which they arise. In this review, we cover large genomic alterations, the mechanisms that cause them and their effect on tumor development and evolution.

CDK4 phosphorylation status and rational use for combining CDK4/6 and BRAF/MEK inhibition in advanced thyroid carcinomas

Background

CDK4/6 inhibitors (CDK4/6i) have been established as standard treatment against advanced Estrogen Receptor-positive breast cancers. These drugs are being tested against several cancers, including in combinations with other therapies. We identified the T172-phosphorylation of CDK4 as the step determining its activity, retinoblastoma protein (RB) inactivation, cell cycle commitment and sensitivity to CDK4/6i. Poorly differentiated (PDTC) and anaplastic (ATC) thyroid carcinomas, the latter considered one of the most lethal human malignancies, represent major clinical challenges. Several molecular evidence suggest that CDK4/6i could be considered for treating these advanced thyroid cancers.

Methods

We analyzed by two-dimensional gel electrophoresis the CDK4 modification profile and the presence of T172-phosphorylated CDK4 in a collection of 98 fresh-frozen tissues and in 21 cell lines. A sub-cohort of samples was characterized by RNA sequencing and immunohistochemistry. Sensitivity to CDK4/6i (palbociclib and abemaciclib) was assessed by BrdU incorporation/viability assays. Treatment of cell lines with CDK4/6i and combination with BRAF/MEK inhibitors (dabrafenib/trametinib) was comprehensively evaluated by western blot, characterization of immunoprecipitated CDK4 and CDK2 complexes and clonogenic assays.

Results

CDK4 phosphorylation was detected in all well-differentiated thyroid carcinomas (n=29), 19/20 PDTC, 16/23 ATC and 18/21 thyroid cancer cell lines, including 11 ATC-derived ones. Tumors and cell lines without phosphorylated CDK4 presented very high p16CDKN2A levels, which were associated with proliferative activity. Absence of CDK4 phosphorylation in cell lines was associated with CDK4/6i insensitivity. RB1 defects (the primary cause of intrinsic CDK4/6i resistance) were not found in 5/7 tumors without detectable phosphorylated CDK4. A previously developed 11-gene expression signature identified the likely unresponsive tumors, lacking CDK4 phosphorylation. In cell lines, palbociclib synergized with dabrafenib/trametinib by completely and permanently arresting proliferation. These combinations prevented resistance mechanisms induced by palbociclib, most notably Cyclin E1-CDK2 activation and a paradoxical stabilization of phosphorylated CDK4 complexes.

Conclusion

Our study supports further clinical evaluation of CDK4/6i and their combination with anti-BRAF/MEK therapies as a novel effective treatment against advanced thyroid tumors. Moreover, the complementary use of our 11 genes predictor with p16/KI67 evaluation could represent a prompt tool for recognizing the intrinsically CDK4/6i insensitive patients, who are potentially better candidates to immediate chemotherapy.

Abstract PR008: Timing complex copy number gains in whole genome duplicated tumors

The accumulation of genetic copy number aberrations is a process common to the development of many cancers. Tumors that have undergone a whole genome duplication (WGD) are known to have high degrees of genomic instability and often contain genomic regions that have undergone a series of copy number gains resulting in multiple copies of both alleles. Metastases are known to have an even greater rate of WGDs and general genomic instability. WGDs provide a powerful handle to time genetic events, such as single nucleotide variants or copy-number aberrations, in tumor evolution. However, the relative timing of copy number gains in whole genome duplicated tumors is complicated by the fact that there are multiple plausible route histories that could give rise to them. Typically, the most parsimonious event history is assumed to take place but this has never been validated. Here we describe a method, GRITIC, that overcomes this problem by evaluating all possible histories that could result in a given copy number state and determining the timing of the most likely series of gains. GRITIC allows the gains leading to any copy number state to be inferred. On a representative simulated cohort, it accurately measures the timing of complex copy number gains. By applying GRITIC to 6,010 primary and metastatic tumor samples from the Pan-Cancer Analysis of Whole Genomes and Hartwig Medical Foundation datasets, we find that the principle of maximum parsimony is violated in at least 25% of all copy number gains in whole genome duplicated tumors, with gains occurring both much earlier and later than thought under this assumption. We also find evidence for punctuated bursts of gains in WGD tumors, independent of the duplication itself. Finally, GRITIC allows us to look at the influence of WGD on copy-number events by comparing pre- and post-WGD copy number aberration landscapes. We find that the frequency with which different chromosome arms are gained before a WGD is highly correlated with the frequency that they are gained post-WGD. In contrast, the correlation between pre- and post-WGD losses is much smaller. We suggest that this indicates that a genome duplication changes the evolutionary trajectories of copy number losses but not gains. As an exception to this trend, chromosome 8 was disproportionately likely to be gained pre-WGD in metastatic prostate cancer as well as 1q in both breast and liver cancers. Altogether, GRITIC allows for a more accurate inference of evolutionary histories in different cancer types and better insights into the early genomic events that occur in tumors that underwent a WGD.

Citation Format: Toby M. Baker, Siqi Lai, Stefan Dentro, Maxime Tarabichi, Peter Van Loo. Timing complex copy number gains in whole genome duplicated tumors [abstract]. In: Proceedings of the AACR Special Conference: Cancer Metastasis; 2022 Nov 14-17; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;83(2 Suppl_2):Abstract nr PR008.

MEDICC2: whole-genome doubling aware copy-number phylogenies for cancer evolution

Aneuploidy, chromosomal instability, somatic copy-number alterations, and whole-genome doubling (WGD) play key roles in cancer evolution and provide information for the complex task of phylogenetic inference. We present MEDICC2, a method for inferring evolutionary trees and WGD using haplotype-specific somatic copy-number alterations from single-cell or bulk data. MEDICC2 eschews simplifications such as the infinite sites assumption, allowing multiple mutations and parallel evolution, and does not treat adjacent loci as independent, allowing overlapping copy-number events. Using simulations and multiple data types from 2780 tumors, we use MEDICC2 to demonstrate accurate inference of phylogenies, clonal and subclonal WGD, and ancestral copy-number states.

Abstract 2286: Profiling the complex rearrangement landscape of sarcoma

Genomic rearrangements are key mutational processes in a number of bone and soft tissue tumors (sarcomas). These events are used for both disease classification and as prognostic and predictive biomarkers. However, the rearrangement architecture and mutational processes giving rise to many of these events remain poorly characterized. Recent data indicate that sarcomas show particularly high frequencies of complex structural rearrangements, including patterns which do not fit those of known mutational mechanisms. As the largest whole genome sequencing (WGS) dataset of sarcomas to date, the Genomics England (GE) 100,000 Genomes Project represents a unique cohort in which to examine these events.

Structural variants (SVs) were identified using five SV callers (Manta, Delly, LUMPY, SvABA and GRIDSS) in WGS data from 1330 GE samples. SVs were filtered to retain only those passing the quality filter, and to remove variants present in either matched germline samples or >1% of a panel of normal variants. Caller performance was evaluated by comparison of distances between SV and copy number call breakpoints generated by Battenberg. SVs identified by at least two callers were taken forward as high-quality calls. Chromothriptic events, extrachromosomal DNA (ecDNA) and breakage-fusion-bridges were identified using established algorithms.

In keeping with previous data, the prevalence of structural variants varies by tumor type, with particularly high numbers of SVs observed in liposarcoma and osteosarcoma. Many tumor types included notable outlier samples; these were more likely to show evidence of chromothripsis, as identified by a dedicated chromothripsis caller. Liposarcomas show the highest frequency of chromothripsis, as described previously. In addition, we have examined some tumor types at higher granularity than previously available, demonstrating directly that the prevalence of chromothripsis varies by histological subtype. For example, chromothriptic events are identified in 100% of well-differentiated liposarcomas, but less than 10% of myxoid liposarcomas. Similarly, the presence of ecDNA varies by tumor type, with the novel observation of particularly high rates in angiosarcoma.

Despite recent advances in histological classification, the survival for patients with sarcoma has remained largely static for over 40 years. Determination of the mutational processes generating these rearrangements will shed light on the pathogenesis of these tumors.

Citation Format: Sara Waise, Tom Lesluyes, Jonas Demeulemeester, Maxime Tarabichi, Sarcoma Genomics England Clinical Interpretation Partnership, Nischalan Pillay, Adrienne M. Flanagan, Peter Van Loo. Profiling the complex rearrangement landscape of sarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2286.

Abstract 1606: Chromosomal instability drives spatial and temporal phenotypic diversity in Schwann cancer cells

Background Malignant peripheral nerve sheath tumors (MPNSTs) are rare, highly aggressive soft tissue sarcomas derived from Schwann cells with poor prognosis. Previous studies, although limited in scale and depth, have shown that these tumors are typically genomically complex, with a moderate burden of single nucleotide variants but extensive copy number aberrations (CNAs).

Methods We applied a multi-omics approach to scrutinize the evolution and heterogeneity of a primary MPNST and five recurrence regions from one patient with an NF1 phenotype. Techniques utilized included bulk whole genome sequencing (WGS), single-cell DNA sequencing (scDNA-seq), single-cell RNA sequencing (scRNA-seq), simultaneous genome and transcriptome sequencing (G&T-seq), spatial transcriptomics, and laser capture microdissection (LCM)-based spatial genomics.

Results We find significant heterogeneity in copy number profiles, suggestive of ongoing chromosomal instability and evolution. We perform in-depth tumor phylogenetic reconstruction from bulk WGS data. scDNA-seq revealed further CNA heterogeneity across and within regions, allowing us to refine the tree down to single-cell resolution. De novo somatic variant calling in scDNA-seq data confirms the reconstructed tree and identifies additional subclonal mutations.

We also profile different populations of tumor and non-tumor cells and confirmed this with a genotyping approach. CNA profiles inferred from scRNA-seq reflect the within-region heterogeneity seen in scDNA-seq. We find tumor cells with similar CNA profiles typically cluster together by their gene expression profiles as well, implying that the majority of gene expression heterogeneity in this tumor is underpinned by copy number changes. By using CNA profiles as a native barcoding system, we link genetic subclones to clusters of cells and their transcriptomes to characterize the effects of gene dosage.

Finally, we explore the spatial relationships between tumor subclones using LCM. We detect heterogeneity within tissue sections evidencing local clonal expansions and CNA events that follow spatial distributions. Using spatial transcriptomics, we show tumor cells are homogeneously admixed with tumor microenvironment populations.

Conclusions Our work demonstrates the power of combining spatial multi-omics at the single-cell and bulk levels to study cancer evolution. This enables a spatio-temporal representation of a tumor’s development with annotation of genetic and transcriptomic events.

Citation Format: Haixi Yan, Jonas Demeulemeester, Annelien Verfaillie, Cristina Cotobal Martin, Tom Kaufmann, Roland Schwarz, Thierry Voet, Simone Zaccaria, Adrienne Flanagan, Maxime Tarabichi, Peter Van Loo. Chromosomal instability drives spatial and temporal phenotypic diversity in Schwann cancer cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1606.

Abstract 1698: Copy number diversity within and across tumor types

Introduction Cancers commonly accrue copy number gains and losses during their development. An improved understanding of their contribution to tumorigenesis can be gained by studying their distribution and evolution within and across tumor types.

Materials and Methods We analyzed copy number profiles of 16,765 tumours from the Hartwig Medical Foundation, Pan-Cancer Analysis of Whole Genomes (PCAWG) and The Cancer Genome Atlas (TCGA) cohorts. We used a novel timing method, ComplexTiming, to time the occurrence of copy number gains, and a permutation-based approach to assess the distribution of copy number events relative to a null distribution.

Results Complex copy number states in cancers that have undergone a whole-genome duplication (WGD) are commonly assumed to develop through a most parsimonious route, where the series of events leading to the complex copy number state requiring the least number of steps, is proposed to be the most likely one. Our newly developed ComplexTiming approach allows unique insights into the actual routes that led to these complex copy number states. In our large pan-cancer cohort, we found that the most parsimonious set of events does not occur in approximately 40% of cases. Additional gains more frequently occurred after WGD than before, reflecting a previously observed increase in instability. In squamous cell lung cancer, gains of 3q were frequently found before WGD. Additionally, 5q and 1q were significantly more likely to be gained before WGD in colorectal and melanoma metastases, respectively.

Evaluating the distribution of copy number events, we found that the majority of copy number events are distributed across tumors consistent with a permutation model that assigns copy number events randomly. However, interestingly, in the majority of cancer types, we found copy number events that occurred more often than expected in individual tumors with low numbers of events, such as loss of chromosome 9 in bladder cancer, loss of 20p in colorectal adenocarcinomas, and gain of 1q in breast adenocarcinomas. These results suggest that the high-frequency of these specific copy number events is due to selection, while other copy number changes may occur neutrally due to innate chromosome instability.

Conclusion The landscape of copy number gains in tumors is diverse. By analyzing the timing of these events as well as their distribution, we identify a subset of frequent copy number changes in different cancer types that are likely under selection.

Citation Format: Toby Baker, Stefan Dentro, Paul Spellman, Maxime Tarabichi, Peter Van Loo. Copy number diversity within and across tumor types [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1698.

A pan-cancer compendium of chromosomal instability

Chromosomal instability (CIN) results in the accumulation of large-scale losses, gains and rearrangements of DNA¹. The broad genomic complexity caused by CIN is a hallmark of cancer²; however, there is no systematic framework to measure different types of CIN and their effect on clinical phenotypes pan-cancer. Here we evaluate the extent, diversity and origin of CIN across 7,880 tumours representing 33 cancer types. We present a compendium of 17 copy number signatures that characterize specific types of CIN, with putative aetiologies supported by multiple independent data sources. The signatures predict drug response and identify new drug targets. Our framework refines the understanding of impaired homologous recombination, which is one of the most therapeutically targetable types of CIN. Our results illuminate a fundamental structure underlying genomic complexity in human cancers and provide a resource to guide future CIN research.

Signatures of copy number alterations in human cancer

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

Leveraging single cell sequencing to unravel intra-tumour heterogeneity and tumour evolution in human cancers

Intra-tumour heterogeneity and tumour evolution are well-documented phenomena in human cancers. While the advent of next-generation sequencing technologies has facilitated the large-scale capture of genomic data, the field of single cell genomics is nascent but rapidly advancing and generating many new insights into the complex molecular mechanisms of tumour biology. In this review, we provide an overview of current single cell DNA sequencing technologies, exploring how recent methodological advancements have enumerated new insights into intra-tumour heterogeneity and tumour evolution. Areas highlighted include the potential power of single cell genome sequencing studies to explore evolutionary dynamics contributing to tumourigenesis through to progression, metastasis and therapy resistance. We also explore the use of in-situ sequencing technologies to study intra-tumour heterogeneity in a spatial context, as well as examining the use of single cell genomics to perform lineage tracing in both normal and malignant tissues. Finally, we consider the use of multi-modal single cell sequencing technologies. Taken together, it is hoped that these many facets of single cell genome sequencing will improve our understanding of tumourigenesis, progression and lethality in cancer leading to the development of novel therapies. This article is protected by copyright. All rights reserved.

Thyroid cancer under the scope of emerging technologies

The vast majority of thyroid cancers originate from follicular cells. We outline outstanding issues at each step along the path of cancer patient care, from prevention to post-treatment follow-up and highlight how emerging technologies will help address them in the coming years. Three directions will dominate the coming technological landscape. Genomics will reveal tumoral evolutionary history and shed light on how these cancers arise from the normal epithelium and the genomics alteration driving their progression. Transcriptomics will gain cellular and spatial resolution providing a full account of intra-tumor heterogeneity and opening a window on the microenvironment supporting thyroid tumor growth. Artificial intelligence will set morphological analysis on an objective quantitative ground laying the foundations of a systematic thyroid tumor classification system. It will also integrate into unified representations the molecular and morphological perspectives on thyroid cancer.

Therapeutic vulnerability to PARP1,2 inhibition in RB1-mutant osteosarcoma

Loss-of-function mutations in the RB1 tumour suppressor are key drivers in cancer, including osteosarcoma. RB1 loss-of-function compromises genome-maintenance and hence could yield vulnerability to therapeutics targeting such processes. Here we demonstrate selective hypersensitivity to clinically-approved inhibitors of Poly-ADP-Polymerase1,2 inhibitors (PARPi) in RB1-defective cancer cells, including an extended panel of osteosarcoma-derived lines. PARPi treatment results in extensive cell death in RB1-defective backgrounds and prolongs survival of mice carrying human RB1-defective osteosarcoma grafts. PARPi sensitivity is not associated with canonical homologous recombination defect (HRd) signatures that predict PARPi sensitivity in cancers with BRCA1,2 loss, but is accompanied by rapid activation of DNA replication checkpoint signalling, and active DNA replication is a prerequisite for sensitivity. Importantly, sensitivity in backgrounds with natural or engineered RB1 loss surpasses that seen in BRCA-mutated backgrounds where PARPi have established clinical benefit. Our work provides evidence that PARPi sensitivity extends beyond cancers identifiable by HRd and advocates PARP1,2 inhibition as a personalised strategy for RB1-mutated osteosarcoma and other cancers.

A pan-cancer landscape of somatic mutations in non-unique regions of the human genome

A substantial fraction of the human genome displays high sequence similarity with at least one other genomic sequence, posing a challenge for the identification of somatic mutations from short-read sequencing data. Here we annotate genomic variants in 2,658 cancers from the Pan-Cancer Analysis of Whole Genomes (PCAWG) cohort with links to similar sites across the human genome. We train a machine learning model to use signals distributed over multiple genomic sites to call somatic events in non-unique regions and validate the data against linked-read sequencing in an independent dataset. Using this approach, we uncover previously hidden mutations in ~1,700 coding sequences and in thousands of regulatory elements, including in known cancer genes, immunoglobulins and highly mutated gene families. Mutations in non-unique regions are consistent with mutations in unique regions in terms of mutation burden and substitution profiles. The analysis provides a systematic summary of the mutation events in non-unique regions at a genome-wide scale across multiple human cancers.

DeCiFering the subclonal composition of tumors

Accurately identifying the subclones that make up tumors is critical for understanding cancer biology. In an article in this issue of Cell Systems, Satas et al. examine mutations with an evolutionary perspective to decipher the composition of tumors.

Chemotherapy induces canalization of cell state in childhood B-cell precursor acute lymphoblastic leukemia

Comparison of intratumor genetic heterogeneity in cancer at diagnosis and relapse suggests that chemotherapy induces bottleneck selection of subclonal genotypes. However, evolutionary events subsequent to chemotherapy could also explain changes in clonal dominance seen at relapse. We, therefore, investigated the mechanisms of selection in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) during induction chemotherapy where maximal cytoreduction occurs. To distinguish stochastic versus deterministic events, individual leukemias were transplanted into multiple xenografts and chemotherapy administered. Analyses of the immediate post-treatment leukemic residuum at single-cell resolution revealed that chemotherapy has little impact on genetic heterogeneity. Rather, it acts on extensive, previously unappreciated, transcriptional and epigenetic heterogeneity in BCP-ALL, dramatically reducing the spectrum of cell states represented, leaving a genetically polyclonal but phenotypically uniform population with hallmark signatures relating to developmental stage, cell cycle and metabolism. Hence, canalization of cell state accounts for a significant component of bottleneck selection during induction chemotherapy.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Breast tumours maintain a reservoir of subclonal diversity during expansion

Our knowledge of copy number evolution during the expansion of primary breast tumours is limited^1,2. Here, to investigate this process, we developed a single-cell, single-molecule DNA-sequencing method and performed copy number analysis of 16,178 single cells from 8 human triple-negative breast cancers and 4 cell lines. The results show that breast tumours and cell lines comprise a large milieu of subclones (7-22) that are organized into a few (3-5) major superclones. Evolutionary analysis suggests that after clonal TP53 mutations, multiple loss-of-heterozygosity events and genome doubling, there was a period of transient genomic instability followed by ongoing copy number evolution during the primary tumour expansion. By subcloning single daughter cells in culture, we show that tumour cells rediversify their genomes and do not retain isogenic properties. These data show that triple-negative breast cancers continue to evolve chromosome aberrations and maintain a reservoir of subclonal diversity during primary tumour growth.

Common clonal origin of conventional T cells and induced regulatory T cells in breast cancer patients

Regulatory CD4+ T cells (Treg) prevent tumor clearance by conventional T cells (Tconv) comprising a major obstacle of cancer immune-surveillance. Hitherto, the mechanisms of Treg repertoire formation in human cancers remain largely unclear. Here, we analyze Treg clonal origin in breast cancer patients using T-Cell Receptor and single-cell transcriptome sequencing. While Treg in peripheral blood and breast tumors are clonally distinct, Tconv clones, including tumor-antigen reactive effectors (Teff), are detected in both compartments. Tumor-infiltrating CD4+ cells accumulate into distinct transcriptome clusters, including early activated Tconv, uncommitted Teff, Th1 Teff, suppressive Treg and pro-tumorigenic Treg. Trajectory analysis suggests early activated Tconv differentiation either into Th1 Teff or into suppressive and pro-tumorigenic Treg. Importantly, Tconv, activated Tconv and Treg share highly-expanded clones contributing up to 65% of intratumoral Treg. Here we show that Treg in human breast cancer may considerably stem from antigen-experienced Tconv converting into secondary induced Treg through intratumoral activation.

A practical guide to cancer subclonal reconstruction from DNA sequencing

Subclonal reconstruction from bulk tumor DNA sequencing has become a pillar of cancer evolution studies, providing insight into the clonality and relative ordering of mutations and mutational processes. We provide an outline of the complex computational approaches used for subclonal reconstruction from single and multiple tumor samples. We identify the underlying assumptions and uncertainties in each step and suggest best practices for analysis and quality assessment. This guide provides a pragmatic resource for the growing user community of subclonal reconstruction methods.

Drivers underpinning the malignant transformation of giant cell tumour of bone

The rare benign giant cell tumour of bone (GCTB) is defined by an almost unique mutation in the H3.3 family of histone genes H3‐3A or H3‐3B; however, the same mutation is occasionally found in primary malignant bone tumours which share many features with the benign variant. Moreover, lung metastases can occur despite the absence of malignant histological features in either the primary or metastatic lesions. Herein we investigated the genetic events of 17 GCTBs including benign and malignant variants and the methylation profiles of 122 bone tumour samples including GCTBs. Benign GCTBs possessed few somatic alterations and no other known drivers besides the H3.3 mutation, whereas all malignant tumours harboured at least one additional driver mutation and exhibited genomic features resembling osteosarcomas, including high mutational burden, additional driver event(s), and a high degree of aneuploidy. The H3.3 mutation was found to predate the development of aneuploidy. In contrast to osteosarcomas, malignant H3.3‐mutated tumours were enriched for a variety of alterations involving TERT, other than amplification, suggesting telomere dysfunction in the transformation of benign to malignant GCTB. DNA sequencing of the benign metastasising GCTB revealed no additional driver alterations; polyclonal seeding in the lung was identified, implying that the metastatic lesions represent an embolic event. Unsupervised clustering of DNA methylation profiles revealed that malignant H3.3‐mutated tumours are distinct from their benign counterpart, and other bone tumours. Differential methylation analysis identified CCND1, encoding cyclin D1, as a plausible cancer driver gene in these tumours because hypermethylation of the CCND1 promoter was specific for GCTBs. We report here the genomic and methylation patterns underlying the rare clinical phenomena of benign metastasising and malignant transformation of GCTB and show how the combination of genomic and epigenomic findings could potentially distinguish benign from malignant GCTBs, thereby predicting aggressive behaviour in challenging diagnostic cases. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Abstract 1320: The landscape of somatic substitutions in the repetitive genome across cancer types

Somatic mutation calling from bulk DNA sequencing is a complex problem susceptible to elevated false positive rates. High mapping quality is considered an important feature of reliable variant calls. At the resolution of short reads, ~10% of the genome displays high sequence similarity with at least one other genomic region and is assigned low mapping quality by alignment algorithms. These low-mapping-quality regions represent recurrent blind spots for mutation callers, which discard many of the variants they harbor, overlooking true biological variation.

Here, we developed a pipeline to call substitutions in the low-mapping-quality genome. We used a published thesaurus approach to annotate the variant positions with their high-similarity links. We trained a classifier to emulate high-quality consensus calls made in unique regions using 20 features unrelated to mapping quality, reaching ~95% accuracy in those regions. In an independent sample more than 90% of the thesaurus calls were validated through linked-read sequencing. We then applied the classifier to all candidate substitutions of 2,658 cancer whole genomes from the PCAWG/ICGC consortium including variants falling in low-mapability regions.

We retrieve hidden thesaurus variants genome-wide in ~6% of the genome, including genic, coding, and promoter regions. Thesaurus calls are directly proportional in numbers to somatic calls falling in the low-mapping-quality genome and share a similar trinucleotide context spectrum. Rescuing these mutations reveal hidden signal in known cancer genes, including PIK3CA, and excess of mutations genome-wide in promoter, untranslated, and coding regions of many other genes. We also find potential excess of non-synonymous mutations, including in genes from the TRIM and POTE families, having been previously implicated in multiple cancer types.

Altogether, we developed a pipeline to call somatic substitutions in the low-mapping-quality genome and uncovered hidden somatic changes along the genomes of human cancers. In the future, this pipeline could be extended to indels and structural variants, and applied to the study of de novo germline variants.

Citation Format: Maxime Tarabichi, Jonas Demeulemeester, Annelien Verfaillie, Peter Van Loo, Tomasz Konopka. The landscape of somatic substitutions in the repetitive genome across cancer types [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1320.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

The evolutionary history of 2,658 cancers

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

A community effort to create standards for evaluating tumor subclonal reconstruction

Tumor DNA sequencing data can be interpreted by computational methods that analyze genomic heterogeneity to infer evolutionary dynamics. A growing number of studies have used these approaches to link cancer evolution with clinical progression and response to therapy. Although the inference of tumor phylogenies is rapidly becoming standard practice in cancer genome analyses, standards for evaluating them are lacking. To address this need, we systematically assess methods for reconstructing tumor subclonality. First, we elucidate the main algorithmic problems in subclonal reconstruction and develop quantitative metrics for evaluating them. Then we simulate realistic tumor genomes that harbor all known clonal and subclonal mutation types and processes. Finally, we benchmark 580 tumor reconstructions, varying tumor read depth, tumor type and somatic variant detection. Our analysis provides a baseline for the establishment of gold-standard methods to analyze tumor heterogeneity.

Investigation of somatic CNVs in brains of synucleinopathy cases using targeted SNCA analysis and single cell sequencing

Synucleinopathies are mostly sporadic neurodegenerative disorders of partly unexplained aetiology, and include Parkinson's disease (PD) and multiple system atrophy (MSA). We have further investigated our recent finding of somatic SNCA (α-synuclein) copy number variants (CNVs, specifically gains) in synucleinopathies, using Fluorescent in-situ Hybridisation for SNCA, and single-cell whole genome sequencing for the first time in a synucleinopathy. In the cingulate cortex, mosaicism levels for SNCA gains were higher in MSA and PD than controls in neurons (> 2% in both diseases), and for MSA also in non-neurons. In MSA substantia nigra (SN), we noted SNCA gains in > 3% of dopaminergic (DA) neurons (identified by neuromelanin) and neuromelanin-negative cells, including olig2-positive oligodendroglia. Cells with CNVs were more likely to have α-synuclein inclusions, in a pattern corresponding to cell categories mostly relevant to the disease: DA neurons in Lewy-body cases, and other cells in the striatonigral degeneration-dominant MSA variant (MSA-SND). Higher mosaicism levels in SN neuromelanin-negative cells may correlate with younger onset in typical MSA-SND, and in cingulate neurons with younger death in PD. Larger sample sizes will, however, be required to confirm these putative findings. We obtained genome-wide somatic CNV profiles from 169 cells from the substantia nigra of two MSA cases, and pons and putamen of one. These showed somatic CNVs in ~ 30% of cells, with clonality and origins in segmental duplications for some. CNVs had distinct profiles based on cell type, with neurons having a mix of gains and losses, and other cells having almost exclusively gains, although control data sets will be required to determine possible disease relevance. We propose that somatic SNCA CNVs may contribute to the aetiology and pathogenesis of synucleinopathies, and that genome-wide somatic CNVs in MSA brain merit further study.

2019 John Insall Award: Fructosamine is a better glycaemic marker compared with glycated haemoglobin (HbA1C) in predicting adverse outcomes following total knee arthroplasty: a prospective multicentre study

AIMS

The best marker for assessing glycaemic control prior to total knee arthroplasty (TKA) remains unknown. The purpose of this study was to assess the utility of fructosamine compared with glycated haemoglobin (HbA1c) in predicting early complications following TKA, and to determine the threshold above which the risk of complications increased markedly.

PATIENTS AND METHODS

This prospective multi-institutional study evaluated primary TKA patients from four academic institutions. Patients (both diabetics and non-diabetics) were assessed using fructosamine and HbA1c levels within 30 days of surgery. Complications were assessed for 12 weeks from surgery and included prosthetic joint infection (PJI), wound complication, re-admission, re-operation, and death. The Youden's index was used to determine the cut-off for fructosamine and HbA1c associated with complications. Two additional cut-offs for HbA1c were examined: 7% and 7.5% and compared with fructosamine as a predictor for complications.

RESULTS

Overall, 1119 patients (441 men, 678 women) were included in the study. Fructosamine level of 293 µmol/l was identified as the optimal cut-off associated with complications. Patients with high fructosamine (> 293 µmol/l) were 11.2 times more likely to develop PJI compared with patients with low fructosamine (p = 0.001). Re-admission and re-operation rates were 4.2 and 4.5 times higher in patients with fructosamine above the threshold (p = 0.005 and p = 0.019, respectively). One patient (1.7%) from the elevated fructosamine group died compared with one patient (0.1%) in the normal fructosamine group (p = 0.10). These complications remained statistically significant in multiple regression analysis. Unlike fructosamine, all three cut-offs for HbA1c failed to show a significant association with complications.

CONCLUSION

Fructosamine is a valid and an excellent predictor of complications following TKA. It better reflects the glycaemic control, has greater predictive power for adverse events, and responds quicker to treatment compared with HbA1c. These findings support the screening of all patients undergoing TKA using fructosamine and in those with a level above 293 µmol/l, the risk of surgery should be carefully weighed against its benefit. Cite this article: Bone Joint J 2019;101-B(7 Supple C):3-9.

Undifferentiated Sarcomas Develop through Distinct Evolutionary Pathways

Undifferentiated sarcomas (USARCs) of adults are diverse, rare, and aggressive soft tissue cancers. Recent sequencing efforts have confirmed that USARCs exhibit one of the highest burdens of structural aberrations across human cancer. Here, we sought to unravel the molecular basis of the structural complexity in USARCs by integrating DNA sequencing, ploidy analysis, gene expression, and methylation profiling. We identified whole genome duplication as a prevalent and pernicious force in USARC tumorigenesis. Using mathematical deconvolution strategies to unravel the complex copy-number profiles and mutational timing models we infer distinct evolutionary pathways of these rare cancers. In addition, 15% of tumors exhibited raised mutational burdens that correlated with gene expression signatures of immune infiltration, and good prognosis.

Abstract 218: The evolutionary history of 2,658 cancers

Cancer develops through a continuous process of somatic evolution. Whole genome sequencing provides a snapshot of the tumor genome at the point of sampling, however, the data can contain information that permits the reconstruction of a tumor's evolutionary past.

Here, we apply such life history analyses on an unprecedented scale, to a set of 2,658 tumors spanning 39 cancer types. We estimated the timing of large chromosomal gains during tumor evolution, by comparing the rates of doubled to non-doubled point mutations within gained regions. Although we find that such events typically occur in the second half of clonal evolution, we also observe distinctive and early chromosomal gains in some cancer types, such as gains of chromosomes 7, 19 and 20 in glioblastoma, and isochromosome 17q in medulloblastoma. By integrating these results with the qualitative timing of individual driver mutations, we obtained an overall ranking, from early to late, of frequent somatic events per cancer type, which both identified novel patterns of tumor evolution, and incorporated additional detail into known models, such as the progression of APC-KRAS-TP53 in colorectal cancer proposed by Vogelstein and Fearon.

To estimate how mutational processes acting on the tumor genome change over time, we classified mutations in each sample according to three broad time periods (early clonal, late clonal, and subclonal), and quantified the activity of mutational signatures in each period. Most mutational processes appear to remain remarkably constant, however, certain signatures show clear and consistent changes during clonal evolution. Particularly, mutational signatures associated with exposure to carcinogens, such as smoking and UV light, tend to decrease over time. In contrast, signatures associated with defective endogenous processes, such as APOBEC mutagenesis and defective double strand break repair, show an increase between early and late phases of tumor evolution.

Making use of clock-like mutational signatures, we converted mutational time estimates for large events, such as whole genome duplication (WGD), and the emergence of the most recent common ancestor (MRCA), into real time estimates, which allowed us to combine our analyses into overall timelines of cancer evolution, per tumor type. For example, the typical timeline of ovarian adenocarcinoma development shows that early tumor evolution is characterized by mutations in TP53, and widespread genome instability, with WGD events taking place on average 8 years prior to diagnosis. In later stages of evolution, signatures of defective repair processes increase, and the MRCA emerges on average 1 year before diagnosis.

Taken together, these data reveal the common and divergent evolutionary trajectories available to a cancer, which might be crucial in understanding specific tumor biology, and in providing new opportunities for early detection and cancer prevention.

Citation Format: Clemency Jolly, Moritz Gerstung, Ignaty Leshchiner, Stefan C. Dentro, Santiago Gonzalez, Thomas J. Mitchell, Yulia Rubanova, Pavana Anur, Daniel Rosebrock, Kaixian Yu, Maxime Tarabichi, Amit Deshwar, Jeff Wintersinger, Kortine Kleinheinz, Ignacio Vásquez-García, Kerstin Haase, Subhajit Sengupta, Geoff Macintyre, Salem Malikic, Nilgun Donmez, Dimitri G. Livitz, Mark Cmero, Jonas Demeulemeester, Steve Schumacher, Yu Fan, Xiaotong Yao, Juhee Lee, Matthias Schlesner, Paul C. Boutros, David D. Bowtell, Hongtu Zhu, Gad Getz, Marcin Imielinski, Rameen Beroukhim, S Cenk Sahinalp, Yuan Ji, Martin Peifer, Florian Markowetz, Ville Mustonen, Ke Juan, Wenyi Wang, Quaid D. Morris, Paul T. Spellman, David C. Wedge, Peter Van Loo, PCAWG Evolution and Heterogeneity Working Group. The evolutionary history of 2,658 cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 218.

Recurrent rearrangements of FOS and FOSB define osteoblastoma

The transcription factor FOS has long been implicated in the pathogenesis of bone tumours, following the discovery that the viral homologue, v-fos, caused osteosarcoma in laboratory mice. However, mutations of FOS have not been found in human bone-forming tumours. Here, we report recurrent rearrangement of FOS and its paralogue, FOSB, in the most common benign tumours of bone, osteoblastoma and osteoid osteoma. Combining whole-genome DNA and RNA sequences, we find rearrangement of FOS in five tumours and of FOSB in one tumour. Extending our findings into a cohort of 55 cases, using FISH and immunohistochemistry, provide evidence of ubiquitous mutation of FOS or FOSB in osteoblastoma and osteoid osteoma. Overall, our findings reveal a human bone tumour defined by mutations of FOS and FOSB.

Distinctive Desmoplastic 3D Morphology Associated With BRAFV600E in Papillary Thyroid Cancers

CONTEXT

Although 60% of papillary thyroid carcinomas are BRAFV600E mutant (PTCV600E), the increased aggressiveness of these cancers is still debated.

OBJECTIVE

For PTCV600E we aimed to further characterize the extent of the stroma and its activation, the three-dimensional (3D) tumor-stroma interface, and the proliferation rates of tumor and stromal fibroblasts.

DESIGN

We analyzed exomes, transcriptomes, and images of 364 papillary thyroid carcinoma (PTCs) from The Cancer Genome Atlas (TCGA), including 211 PTCV600E; stained 22 independent PTCs for BRAFV600E and Ki67; sequenced the exomes and stained BRAFV600E in 5 primary tumor blocks and 4 nodal metastases from one patient with PTCV600E; and reconstructed the 3D volumes of one tumor and one metastatic block at histological resolution.

RESULTS

In TCGA, BRAFV600E was associated with higher expression of proliferation markers and lower expression of thyroid differentiation markers, independently of tumor purity. Moreover, PTCV600E, in line with their overall lower purity, also had higher expression of fibroblast- and T cell-associated genes and presented more fibrosis. Tumor cells that appeared disconnected on two-dimensional histological slices were revealed to be part of a unique tumor component in the 3D reconstructed microvolumes, and they formed a surprisingly complex connected space, infiltrating a proliferative stroma. Finally, in our PTC set, both stromal fibroblasts and tumor cells presented higher proliferation rates in PTCV600E.

CONCLUSIONS

Our results support the increased aggressiveness associated with BRAFV600E in PTC and shed light on the important role of the stroma in tumor expansion. The greater and more active fibrotic component predicts better efficiency of combined targeted treatments, as previously proposed for melanomaV600E.

Can next generation sequencing play a role in detecting pathogens in synovial fluid?

Aims

The diagnosis of periprosthetic joint infection can be difficult due to the high rate of culture-negative infections. The aim of this study was to assess the use of next-generation sequencing for detecting organisms in synovial fluid.

Materials and Methods

In this prospective, single-blinded study, 86 anonymized samples of synovial fluid were obtained from patients undergoing aspiration of the hip or knee as part of the investigation of a periprosthetic infection. A panel of synovial fluid tests, including levels of C-reactive protein, human neutrophil elastase, total neutrophil count, alpha-defensin, and culture were performed prior to next-generation sequencing.

Results

Of these 86 samples, 30 were alpha-defensin-positive and culture-positive (Group I), 24 were alpha-defensin-positive and culture-negative (Group II) and 32 were alpha-defensin-negative and culture-negative (Group III). Next-generation sequencing was concordant with 25 results for Group I. In four of these, it detected antibiotic resistant bacteria whereas culture did not. In another four samples with relatively low levels of inflammatory biomarkers, culture was positive but next-generation sequencing was negative. A total of ten samples had a positive next-generation sequencing result and a negative culture. In five of these, alpha-defensin was positive and the levels of inflammatory markers were high. In the other five, alpha-defensin was negative and the levels of inflammatory markers were low. While next-generation sequencing detected several organisms in each sample, in most samples with a higher probability of infection, there was a predominant organism present, while in those presumed not to be infected, many organisms were identified with no predominant organism.

Conclusion

Pathogens causing periprosthetic infection in both culture-positive and culture-negative samples of synovial fluid could be identified by next-generation sequencing. Cite this article: Bone Joint J 2018;100-B:127–33.

miRNA expression and function in thyroid carcinomas: a comparative and critical analysis and a model for other cancers

As in many cancer types, miRNA expression profiles and functions have become an important field of research on non-medullary thyroid carcinomas, the most common endocrine cancers. This could lead to the establishment of new diagnostic tests and new cancer therapies. However, different studies showed important variations in their research strategies and results. In addition, the action of miRNAs is poorly considered as a whole because of the use of underlying dogmatic truncated concepts. These lead to discrepancies and limits rarely considered. Recently, this field has been enlarged by new miRNA functional and expression studies. Moreover, studies using next generation sequencing give a new view of general miRNA differential expression profiles of papillary thyroid carcinoma. We analyzed in detail this literature from both physiological and differential expression points of view. Based on explicit examples, we reviewed the progresses but also the discrepancies and limits trying to provide a critical approach of where this literature may lead. We also provide recommendations for future studies. The conclusions of this systematic analysis could be extended to other cancer types.

The human bitter taste receptor T2R38 is broadly tuned for bacterial compounds

T2R38 has been shown to be a specific bacterial detector implicated in innate immune defense mechanism of human upper airway. Several clinical studies have demonstrated that this receptor is associated with the development of chronic rhinosinusitis (CRS). T2R38 was previously reported to bind to homoserine lactones (HSL), quorum sensing molecules specific of Pseudomonas Aeruginosa and other gram negative species. Nevertheless, these bacteria are not the major pathogens found in CRS. Here we report on the identification of bacterial metabolites acting as new agonists of T2R38 based on a single cell calcium imaging study. Two quorum sensing molecules (Agr D1 thiolactone from Staphylococcus Aureus and CSP-1 from Streptococcus Pneumoniae) and a list of 32 bacterial metabolites from pathogens frequently implicated in CRS were tested. First, we observed that HSL failed to activate T2R38 in our experimental system, but that the dimethylsulfoxide (DMSO), used as a solvent for these lactones may, by itself, account for the agonistic effect previously described. Secondly, we showed that both Agr D1 thiolactone and CSP-1 are inactive but that at least 7 bacterial metabolites (acetone, 2-butanone, 2-pentanone, 2-methylpropanal, dimethyl disulfide, methylmercaptan, γ-butyrolactone) are able to specifically trigger this receptor. T2R38 is thus much more broadly tuned for bacterial compounds than previously thought.

A research note regarding "Variation in cancer risk among tissues can be explained by the number of stem cell divisions"

Tomasetti and Vogelstein argued that 2/3 of human cancers are due to 'bad luck' and that "primary prevention measures [against cancer] are not likely to be very effective". We demonstrate that their calculations for hepatocellular carcinomas overlooked a major subset of these cancers proven to be preventable through vaccination. The problem, which is not limited to hepatocellular carcinoma, arises from the general reliance of their analysis on average incidences in the United States and the omission of incidences in specific risk groups.

Intratumor heterogeneity and clonal evolution in an aggressive papillary thyroid cancer and matched metastases

The contribution of intratumor heterogeneity to thyroid metastatic cancers is still unknown. The clonal relationships between the primary thyroid tumors and lymph nodes (LN) or distant metastases are also poorly understood. The objective of this study was to determine the phylogenetic relationships between matched primary thyroid tumors and metastases. We searched for non-synonymous single-nucleotide variants (nsSNVs), gene fusions, alternative transcripts, and loss of heterozygosity (LOH) by paired-end massively parallel sequencing of cDNA (RNA-Seq) in a patient diagnosed with an aggressive papillary thyroid cancer (PTC). Seven tumor samples from a stage IVc PTC patient were analyzed by RNA-Seq: two areas from the primary tumor, four areas from two LN metastases, and one area from a pleural metastasis (PLM). A large panel of other thyroid tumors was used for Sanger sequencing screening. We identified seven new nsSNVs. Some of these were early events clonally present in both the primary PTC and the three matched metastases. Other nsSNVs were private to the primary tumor, the LN metastases and/or the PLM. Three new gene fusions were identified. A novel cancer-specific KAZN alternative transcript was detected in this aggressive PTC and in dozens of additional thyroid tumors. The PLM harbored an exclusive whole-chromosome 19 LOH. We have presented the first, to our knowledge, deep sequencing study comparing the mutational spectra in a PTC and both LN and distant metastases. This study has yielded novel findings concerning intra-tumor heterogeneity, clonal evolution and metastases dissemination in thyroid cancer.

Systems biology of cancer: entropy, disorder, and selection-driven evolution to independence, invasion and "swarm intelligence"

Our knowledge of the biology of solid cancer has greatly progressed during the last few years, and many excellent reviews dealing with the various aspects of this biology have appeared. In the present review, we attempt to bring together these subjects in a general systems biology narrative. It starts from the roles of what we term entropy of signaling and noise in the initial oncogenic events, to the first major transition of tumorigenesis: the independence of the tumor cell and the switch in its physiology, i.e., from subservience to the organism to its own independent Darwinian evolution. The development after independence involves a constant dynamic reprogramming of the cells and the emergence of a sort of collective intelligence leading to invasion and metastasis and seldom to the ultimate acquisition of immortality through inter-individual infection. At each step, the probability of success is minimal to infinitesimal, but the number of cells possibly involved and the time scale account for the relatively high occurrence of tumorigenesis and metastasis in multicellular organisms.

Profiling of olfactory receptor gene expression in whole human olfactory mucosa

Olfactory perception is mediated by a large array of olfactory receptor genes. The human genome contains 851 olfactory receptor gene loci. More than 50% of the loci are annotated as nonfunctional due to frame-disrupting mutations. Furthermore haplotypic missense alleles can be nonfunctional resulting from substitution of key amino acids governing protein folding or interactions with signal transduction components. Beyond their role in odor recognition, functional olfactory receptors are also required for a proper targeting of olfactory neuron axons to their corresponding glomeruli in the olfactory bulb. Therefore, we anticipate that profiling of olfactory receptor gene expression in whole human olfactory mucosa and analysis in the human population of their expression should provide an opportunity to select the frequently expressed and potentially functional olfactory receptors in view of a systematic deorphanization. To address this issue, we designed a TaqMan Low Density Array (Applied Biosystems), containing probes for 356 predicted human olfactory receptor loci to investigate their expression in whole human olfactory mucosa tissues from 26 individuals (13 women, 13 men; aged from 39 to 81 years, with an average of 67±11 years for women and 63±12 years for men). Total RNA isolation, DNase treatment, RNA integrity evaluation and reverse transcription were performed for these 26 samples. Then 384 targeted genes (including endogenous control genes and reference genes specifically expressed in olfactory epithelium for normalization purpose) were analyzed using the same real-time reverse transcription PCR platform. On average, the expression of 273 human olfactory receptor genes was observed in the 26 selected whole human olfactory mucosa analyzed, of which 90 were expressed in all 26 individuals. Most of the olfactory receptors deorphanized to date on the basis of sensitivity to known odorant molecules, which are described in the literature, were found in the expressed olfactory receptors gene set.

Comparative analysis of the thyrocytes and T cells: responses to H2O2 and radiation reveals an H2O2-induced antioxidant transcriptional program in thyrocytes

CONTEXT

Radiation is an established cause of thyroid cancer, and growing evidence supports a role for hydrogen peroxide (H2O2) in spontaneous thyroid carcinogenesis. Little is known about the molecular programs activated by these agents in thyrocytes.

OBJECTIVE

The purpose of this study was to compare the responses of thyrocytes and T cells to H2O2 and radiation.

METHODS

We profiled the DNA damage and cell death induced by γ-radiation (0.1-5 Gy) and H2O2 (0.0025-0.3 mM) in primary human thyrocytes and T cells. We next prepared thyroid and T-cell primary cultures from 8 donors operated for noncancerous thyroid pathological conditions and profiled their genome-wide transcriptional response 4 hours after (1) exposure to 1-Gy radiation, (2) treatment with H2O2 and (3) no treatment. Two H2O2 concentrations were investigated, calibrated in each cell type to elicit levels of single- and double-strand breaks equivalent to 1-Gy γ-radiation.

RESULTS

Although thyrocytes and T cells had comparable radiation responses, 3- to 10-fold more H2O2 was needed to induce detectable DNA damage in thyrocytes. At H2O2 and radiation doses inducing double-strand breaks, cell death occurred after 24 hours in T cells but not in thyrocytes. The transcriptional responses of thyrocytes and T cells to radiation were similar, involving DNA repair and cell death genes. In addition to this transcriptional program, H2O2 also up-regulated antioxidant genes in thyrocytes, including glutathione peroxidases and heme oxygenase at the double-strand breaks-inducing concentration. In contrast, a transcriptional storm involving thousands of genes was raised in T cells. Finally, we showed that inhibiting glutathione peroxidases activity increased the DNA damaging effect of H2O2 in thyrocytes.

CONCLUSION

We propose that high H2O2 production in thyrocytes is matched with specific transcriptionally regulated antioxidant protection.

5-aza-2'-deoxycytidine has minor effects on differentiation in human thyroid cancer cell lines, but modulates genes that are involved in adaptation in vitro

BACKGROUND

In thyroid cancer, the lack of response to specific treatment, for example, radioactive iodine, can be caused by a loss of differentiation characteristics of tumor cells. It is hypothesized that this loss is due to epigenetic modifications. Therefore, drugs releasing epigenetic repression have been proposed to reverse this silencing.

METHODS

We investigated which genes were reinduced in dedifferentiated human thyroid cancer cell lines when treated with the demethylating agent 5-aza-2'-deoxycytidine (5-AzadC) and the histone deacetylase inhibitors trichostatin A (TSA) and suberoylanilide hydroxamic acid, by using reverse transcriptase-polymerase chain reaction and microarrays. These results were compared to the expression patterns in in vitro human differentiated thyrocytes and in in vivo dedifferentiated thyroid cancers. In addition, the effects of 5-AzadC on DNA quantities and cell viability were investigated.

RESULTS

Among the canonical thyroid differentiation markers, most were not, or only to a minor extent, re-expressed by 5-AzadC, whether or not combined with TSA or forskolin, an inducer of differentiation in normal thyrocytes. Furthermore, 5-AzadC-modulated overall mRNA expression profiles showed only few commonly regulated genes compared to differentiated cultured primary thyrocytes. In addition, most of the commonly strongly 5-AzadC-induced genes in cell lines were either not regulated or upregulated in anaplastic thyroid carcinomas. Further analysis of which genes were induced by 5-AzadC showed that they were involved in pathways such as apoptosis, antigen presentation, defense response, and cell migration. A number of these genes had similar expression responses in 5-AzadC-treated nonthyroid cell lines.

CONCLUSIONS

Our results suggest that 5-AzadC is not a strong inducer of differentiation in thyroid cancer cell lines. Under the studied conditions and with the model used, 5-AzadC treatment does not appear to be a potential redifferentiation treatment for dedifferentiated thyroid cancer. However, this may reflect primarily the inadequacy of the model rather than that of the treatment. Moreover, the observation that 5-AzadC negatively affected cell viability in cell lines could still suggest a therapeutic opportunity. Some of the genes that were modulated by 5-AzadC were also induced in nonthyroid cancer cell lines, which might be explained by an epigenetic modification resulting in the adaptation of the cell lines to their culture conditions.

Piecewise polynomial representations of genomic tracks

Genomic data from micro-array and sequencing projects consist of associations of measured values to chromosomal coordinates. These associations can be thought of as functions in one dimension and can thus be stored, analyzed, and interpreted as piecewise-polynomial curves. We present a general framework for building piecewise polynomial representations of genome-scale signals and illustrate some of its applications via examples. We show that piecewise constant segmentation, a typical step in copy-number analyses, can be carried out within this framework for both array and (DNA) sequencing data offering advantages over existing methods in each case. Higher-order polynomial curves can be used, for example, to detect trends and/or discontinuities in transcription levels from RNA-seq data. We give a concrete application of piecewise linear functions to diagnose and quantify alignment quality at exon borders (splice sites). Our software (source and object code) for building piecewise polynomial models is available at http://sourceforge.net/projects/locsmoc/.

Thyroid cancer cell lines: an overview

Human thyroid cancer cell lines are the most used models for thyroid cancer studies. They must be used with detailed knowledge of their characteristics. These in vitro cell lines originate from differentiated and dedifferentiated in vivo human thyroid tumors. However, it has been shown that mRNA expression profiles of these cell lines were closer to dedifferentiated in vivo thyroid tumors (anaplastic thyroid carcinoma, ATC) than to differentiated ones. Here an overview of the knowledge of these models was made. The mutational status of six human thyroid cancer cell lines (WRO, FTC133, BCPAP, TPC1, K1, and 8505C) was in line with previously reported findings for 10 genes frequently mutated in thyroid cancer. However, the presence of a BRAF mutation (T1799A: V600E) in WRO questions the use of this cell line as a model for follicular thyroid carcinoma (FTC). Next, to investigate the biological meaning of the modulated mRNAs in these cells, a pathway analysis on previously obtained mRNA profiles was performed on five cell lines. In five cell lines, the MHC class II pathway was down-regulated and in four of them, ribosome biosynthesis and translation pathways were up-regulated. mRNA expression profiles of the cell lines were also compared to those of the different types of thyroid cancers. Three datasets originating from different microarray platforms and derived from distinct laboratories were used. This meta-analysis showed a significant higher correlation between the profiles of the thyroid cancer cell lines and ATC, than to differentiated thyroid tumors (i.e., PTC or FTC) specifically for DNA replication. This already observed higher correlation was obtained here with an increased number of in vivo tumors and using different platforms. In summary, this would suggest that some papillary thyroid carcinoma or follicular thyroid carcinoma (PTC or FTC) cell lines (i.e., TPC-1) might have partially lost their original DNA synthesis/replication regulation mechanisms during their in vitro cell adaptation/evolution.