Spatial Transcriptomics Reveal Pitfalls and Opportunities for the Detection of Rare High-Plasticity Breast Cancer Subtypes.

Breast cancer is one of the most prominent types of cancers, in which therapeutic resistance is a major clinical concern. Specific subtypes, such as claudin-low and metaplastic breast carcinoma (MpBC), have been associated with high nongenetic plasticity, which can facilitate resistance. The similarities and differences between these orthogonal subtypes, identified by molecular and histopathological analyses, respectively, remain insufficiently characterized. Furthermore, adequate methods to identify high-plasticity tumors to better anticipate resistance are lacking. Here, we analyzed 11 triple-negative breast tumors, including 3 claudin-low and 4 MpBC, via high-resolution spatial transcriptomics. We combined pathological annotations and deconvolution approaches to precisely identify tumor spots, on which we performed signature enrichment, differential expression, and copy number analyses. We used The Cancer Genome Atlas and Cancer Cell Line Encyclopedia public databases for external validation of expression markers. By focusing our spatial transcriptomic analyses on tumor cells in MpBC samples, we bypassed the negative impact of stromal contamination and identified specific markers that are neither expressed in other breast cancer subtypes nor expressed in stromal cells. Three markers (BMPER, POPDC3, and SH3RF3) were validated in external expression databases encompassing bulk tumor material and stroma-free cell lines. We unveiled that existing bulk expression signatures of high-plasticity breast cancers are relevant in mesenchymal transdifferentiated compartments but can be hindered by abundant stromal cells in tumor samples, negatively impacting their clinical applicability. Spatial transcriptomic analyses constitute powerful tools to identify specific expression markers and could thus enhance diagnosis and clinical care of rare high-plasticity breast cancers.

Evolutionary history of human colitis-associated colorectal cancer.

OBJECTIVE: IBD confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multiregion sequencing. DESIGN: Exome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 patients with CA-CRC (n=55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using single nucleotide polymorphism arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n=9), low-grade dysplasia (LGD; n=30), high-grade dysplasia (HGD; n=13), mixed LGD/HGD (n=7) and CA-CRC (n=19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC. RESULTS: 10/12 tumours were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNA) per Mb, ~20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated ageing-associated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighbouring CA-CRC, indicating a precancer mutational field. CA-CRCs were often near triploid (40%) or near tetraploid (20%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated 'catastrophic' CNA increase. CONCLUSIONS: Evolutionary genomic analysis revealed precancer clones bearing extensive SNAs and CNAs, with progression to cancer involving a dramatic accrual of CNAs at HGD. Detection of the cancerised field is an encouraging prospect for surveillance, but punctuated evolution may limit the window for early detection.

Evolution of oesophageal adenocarcinoma from metaplastic columnar epithelium without goblet cells in Barrett's oesophagus.

OBJECTIVE: Barrett's oesophagus commonly presents as a patchwork of columnar metaplasia with and without goblet cells in the distal oesophagus. The presence of metaplastic columnar epithelium with goblet cells on oesophageal biopsy is a marker of cancer progression risk, but it is unclear whether clonal expansion and progression in Barrett's oesophagus is exclusive to columnar epithelium with goblet cells. DESIGN: We developed a novel method to trace the clonal ancestry of an oesophageal adenocarcinoma across an entire Barrett's segment. Clonal expansions in Barrett's mucosa were identified using cytochrome c oxidase enzyme histochemistry. Somatic mutations were identified through mitochondrial DNA sequencing and single gland whole exome sequencing. RESULTS: By tracing the clonal origin of an oesophageal adenocarcinoma across an entire Barrett's segment through a combination of histopathological spatial mapping and clonal ordering, we find that this cancer developed from a premalignant clonal expansion in non-dysplastic ('cardia-type') columnar metaplasia without goblet cells. CONCLUSION: Our data demonstrate the premalignant potential of metaplastic columnar epithelium without goblet cells in the context of Barrett's oesophagus.

Derivation of genetic biomarkers for cancer risk stratification in Barrett's oesophagus: a prospective cohort study.

OBJECTIVE: The risk of developing adenocarcinoma in non-dysplastic Barrett's oesophagus is low and difficult to predict. Accurate tools for risk stratification are needed to increase the efficiency of surveillance. We aimed to develop a prediction model for progression using clinical variables and genetic markers. METHODS: In a prospective cohort of patients with non-dysplastic Barrett's oesophagus, we evaluated six molecular markers: p16, p53, Her-2/neu, 20q, MYC and aneusomy by DNA fluorescence in situ hybridisation on brush cytology specimens. Primary study outcomes were the development of high-grade dysplasia or oesophageal adenocarcinoma. The most predictive clinical variables and markers were determined using Cox proportional-hazards models, receiver operating characteristic curves and a leave-one-out analysis. RESULTS: A total of 428 patients participated (345 men; median age 60 years) with a cumulative follow-up of 2019 patient-years (median 45 months per patient). Of these patients, 22 progressed; nine developed high-grade dysplasia and 13 oesophageal adenocarcinoma. The clinical variables, age and circumferential Barrett's length, and the markers, p16 loss, MYC gain and aneusomy, were significantly associated with progression on univariate analysis. We defined an 'Abnormal Marker Count' that counted abnormalities in p16, MYC and aneusomy, which significantly improved risk prediction beyond using just age and Barrett's length. In multivariate analysis, these three factors identified a high-risk group with an 8.7-fold (95% CI 2.6 to 29.8) increased HR when compared with the low-risk group, with an area under the curve of 0.76 (95% CI 0.66 to 0.86). CONCLUSIONS: A prediction model based on age, Barrett's length and the markers p16, MYC and aneusomy determines progression risk in non-dysplastic Barrett's oesophagus.

Regulation of the Expression of Heparan Sulfate 3-O-Sulfotransferase 3B (HS3ST3B) by Inflammatory Stimuli in Human Monocytes.

Heparan sulfate (HS) is recognized as an important player in a wide range of dynamic steps of inflammatory reactions. Thereby, structural HS remodeling is likely to play an important role in the regulation of inflammatory and immune responses; however, little is known about underlying mechanism. In this study, we analyzed the regulation of expression of HS 3-O-sulfotransferases (HS3STs) in response to inflammatory stimuli. We found that among the seven HS3ST isoenzymes, only the expression of HS3ST3B was markedly up-regulated in human primary monocytes and the related cell line THP1 after exposure to TLR agonists. TNF-α was also efficient, to a lesser extent, to increase HS3ST3B expression, while IL-6, IL-4, and IFN-γ were poor inducers. We then analyzed the molecular mechanisms that regulate the high expression of HS3ST3B in response to LPS. Based on the expression of HS3ST3B transcripts and on the response of a reporter gene containing the HS3ST3B1 promoter, we provide evidence that LPS induces a rapid and strong transcription of HS3ST3B1 gene, which was mainly dependent on the activation of NF-κB and JNK signaling pathways. Additionally, active p38 MAPK and de novo synthesized proteins are involved in post-transcriptional mechanisms to maintain a high level of HS3ST3B mRNA to a steady state. Altogether, our findings indicate that HS3ST3B1 gene behaves as a primary response gene, suggesting that it may play an important role in making 3-O-sulfated HS with specific functions in the regulation of inflammatory and immune responses. J. Cell. Biochem. 117: 1529-1542, 2016. © 2015 Wiley Periodicals, Inc.

Characterization of hepatitis C virus interaction with heparan sulfate proteoglycans.

UNLABELLED: Hepatitis C virus (HCV) entry involves binding to cell surface heparan sulfate (HS) structures. However, due to the lipoprotein-like structure of HCV, the exact contribution of virion components to this interaction remains controversial. Here, we investigated the relative contribution of HCV envelope proteins and apolipoprotein E in the HS-binding step. Deletion of hypervariable region 1, a region previously proposed to be involved in HS binding, did not alter HCV virion binding to HS, indicating that this region is not involved in this interaction in the context of a viral infection. Patient sera and monoclonal antibodies recognizing different regions of HCV envelope glycoproteins were also used in a pulldown assay with beads coated with heparin, a close HS structural homologue. Although isolated HCV envelope glycoproteins could interact with heparin, none of these antibodies was able to interfere with the virion-heparin interaction, strongly suggesting that at the virion surface, HCV envelope glycoproteins are not accessible for HS binding. In contrast, results from kinetic studies, heparin pulldown experiments, and inhibition experiments with anti-apolipoprotein E antibodies indicated that this apolipoprotein plays a major role in HCV-HS interaction. Finally, characterization of the HS structural determinants required for HCV infection by silencing of the enzymes involved in the HS biosynthesis pathway and by competition with modified heparin indicated that N- and 6-O-sulfation but not 2-O-sulfation is required for HCV infection and that the minimum HS oligosaccharide length required for HCV infection is a decasaccharide. Together, these data indicate that HCV hijacks apolipoprotein E to initiate its interaction with specific HS structures. IMPORTANCE: Hepatitis C is a global health problem. Hepatitis C virus (HCV) infects approximately 130 million individuals worldwide, with the majority of cases remaining undiagnosed and untreated. In most infected individuals, the virus evades the immune system and establishes a chronic infection. As a consequence, hepatitis C is the leading cause of cirrhosis, end-stage liver disease, hepatocellular carcinoma, and liver transplantation. Virus infection is initiated by entry of the virus into the host cell. In this study, we provide new insights into the viral and cellular determinants involved in the first step of HCV entry, the binding of the virus to host cells. We show that apolipoprotein E is likely responsible for virus binding to heparan sulfate and that N- and 6-O-sulfation of the heparan sulfate proteoglycans is required for HCV infection. In addition, the minimal HS length unit required for HCV infection is a decasaccharide.

Macrophage polarization alters the expression and sulfation pattern of glycosaminoglycans.

Macrophages are major cells of inflammatory process and take part in a large number of physiological and pathological processes. According to tissue environment, they can polarize into pro-inflammatory (M1) or alternative (M2) cells. Although many evidences have hinted to a potential role of cell-surface glycosaminoglycans (GAGs) in the functions of macrophages, the effect of M1 or M2 polarization on the biosynthesis of these polysaccharides has not been investigated so far. GAGs are composed of repeat sulfated disaccharide units. Heparan (HS) and chondroitin/dermatan sulfates (CS/DS) are the major GAGs expressed at the cell membrane. They are involved in numerous biological processes, which rely on their ability to selectively interact with a large panel of proteins. More than 20 genes encoding sulfotransferases have been implicated in HS and CS/DS biosynthesis, and the functional repertoire of HS and CS/DS has been related to the expression of these isoenzymes. In this study, we analyzed the expression of sulfotransferases as a response to macrophage polarization. We found that M1 and M2 activation drastically modified the profiles of expression of numerous HS and CS/DS sulfotransferases. This was accompanied by the expression of GAGs with distinct structural features. We then demonstrated that GAGs of M2 macrophages were efficient to present fibroblast growth factor-2 in an assay of tumor cell proliferation, thus indicating that changes in GAG structure may contribute to the functions of polarized macrophages. Altogether, our findings suggest a regulatory mechanism in which fine modifications in GAG biosynthesis may participate to the plasticity of macrophage functions.

Intratumor heterogeneity and branched evolution revealed by multiregion sequencing.

BACKGROUND: Intratumor heterogeneity may foster tumor evolution and adaptation and hinder personalized-medicine strategies that depend on results from single tumor-biopsy samples. METHODS: To examine intratumor heterogeneity, we performed exome sequencing, chromosome aberration analysis, and ploidy profiling on multiple spatially separated samples obtained from primary renal carcinomas and associated metastatic sites. We characterized the consequences of intratumor heterogeneity using immunohistochemical analysis, mutation functional analysis, and profiling of messenger RNA expression. RESULTS: Phylogenetic reconstruction revealed branched evolutionary tumor growth, with 63 to 69% of all somatic mutations not detectable across every tumor region. Intratumor heterogeneity was observed for a mutation within an autoinhibitory domain of the mammalian target of rapamycin (mTOR) kinase, correlating with S6 and 4EBP phosphorylation in vivo and constitutive activation of mTOR kinase activity in vitro. Mutational intratumor heterogeneity was seen for multiple tumor-suppressor genes converging on loss of function; SETD2, PTEN, and KDM5C underwent multiple distinct and spatially separated inactivating mutations within a single tumor, suggesting convergent phenotypic evolution. Gene-expression signatures of good and poor prognosis were detected in different regions of the same tumor. Allelic composition and ploidy profiling analysis revealed extensive intratumor heterogeneity, with 26 of 30 tumor samples from four tumors harboring divergent allelic-imbalance profiles and with ploidy heterogeneity in two of four tumors. CONCLUSIONS: Intratumor heterogeneity can lead to underestimation of the tumor genomics landscape portrayed from single tumor-biopsy samples and may present major challenges to personalized-medicine and biomarker development. Intratumor heterogeneity, associated with heterogeneous protein function, may foster tumor adaptation and therapeutic failure through Darwinian selection. (Funded by the Medical Research Council and others.).

Parallel evolution of tumour subclones mimics diversity between tumours.

Intratumour heterogeneity (ITH) may foster tumour adaptation and compromise the efficacy of personalized medicine approaches. The scale of heterogeneity within a tumour (intratumour heterogeneity) relative to genetic differences between tumours (intertumour heterogeneity) is unknown. To address this, we obtained 48 biopsies from eight stage III and IV clear cell renal cell carcinomas (ccRCCs) and used DNA copy-number analyses to compare biopsies from the same tumour with 440 single tumour biopsies from the Cancer Genome Atlas (TCGA). Unsupervised hierarchical clustering of TCGA and multi-region ccRCC samples revealed segregation of samples from the same tumour into unrelated clusters; 25% of multi-region samples appeared more similar to unrelated samples than to any other sample originating from the same tumour. We found that the majority of recurrent DNA copy number driver aberrations in single biopsies were not present ubiquitously in late-stage ccRCCs and were likely to represent subclonal events acquired during tumour progression. Such heterogeneous subclonal genetic alterations within individual tumours may impair the identification of robust ccRCC molecular subtypes classified by distinct copy number alterations and clinical outcomes. The co-existence of distinct subclonal copy number events in different regions of individual tumours reflects the diversification of individual ccRCCs through multiple evolutionary routes and may contribute to tumour sampling bias and impact upon tumour progression and clinical outcome.

Cyclophilin B attenuates the expression of TNF-α in lipopolysaccharide-stimulated macrophages through the induction of B cell lymphoma-3.

Extracellular cyclophilin A (CyPA) and CyPB have been well described as chemotactic factors for various leukocyte subsets, suggesting their contribution to inflammatory responses. Unlike CyPA, CyPB accumulates in extracellular matrixes, from which it is released by inflammatory proteases. Hence, we hypothesized that it could participate in tissue inflammation by regulating the activity of macrophages. In the current study, we confirmed that CyPB initiated in vitro migration of macrophages, but it did not induce production of proinflammatory cytokines. In contrast, pretreatment of macrophages with CyPB attenuated the expression of inflammatory mediators induced by LPS stimulation. The expression of TNF-α mRNA was strongly reduced after exposure to CyPB, but it was not accompanied by significant modification in LPS-induced activation of MAPK and NF-κB pathways. LPS activation of a reporter gene under the control of TNF-α gene promoter was also markedly decreased in cells treated with CyPB, suggesting a transcriptional mechanism of inhibition. Consistent with this hypothesis, we found that CyPB induced the expression of B cell lymphoma-3 (Bcl-3), which was accompanied by a decrease in the binding of NF-κB p65 to the TNF-α promoter. As expected, interfering with the expression of Bcl-3 restored cell responsiveness to LPS, thus confirming that CyPB acted by inhibiting initiation of TNF-α gene transcription. Finally, we found that CyPA was not efficient in attenuating the production of TNF-α from LPS-stimulated macrophages, which seemed to be due to a modest induction of Bcl-3 expression. Collectively, these findings suggest an unexpected role for CyPB in attenuation of the responses of proinflammatory macrophages.

Genome-wide RNA interference analysis of renal carcinoma survival regulators identifies MCT4 as a Warburg effect metabolic target.

Clear cell renal cell carcinoma (ccRCC) is the most common pathological subtype of kidney cancer. Here, we integrated an unbiased genome-wide RNA interference screen for ccRCC survival regulators with an analysis of recurrently overexpressed genes in ccRCC to identify new therapeutic targets in this disease. One of the most potent survival regulators, the monocarboxylate transporter MCT4 (SLC16A3), impaired ccRCC viability in all eight ccRCC lines tested and was the seventh most overexpressed gene in a meta-analysis of five ccRCC expression datasets. MCT4 silencing impaired secretion of lactate generated through glycolysis and induced cell cycle arrest and apoptosis. Silencing MCT4 resulted in intracellular acidosis, and reduction in intracellular ATP production together with partial reversion of the Warburg effect in ccRCC cell lines. Intra-tumoural heterogeneity in the intensity of MCT4 protein expression was observed in primary ccRCCs. MCT4 protein expression analysis based on the highest intensity of expression in primary ccRCCs was associated with poorer relapse-free survival, whereas modal intensity correlated with Fuhrman nuclear grade. Consistent with the potential selection of subclones enriched for MCT4 expression during disease progression, MCT4 expression was greater at sites of metastatic disease. These data suggest that MCT4 may serve as a novel metabolic target to reverse the Warburg effect and limit disease progression in ccRCC.

Identification of a Gene-Expression-Based Surrogate of Genomic Instability during Oral Carcinogenesis.

BACKGROUND: Our goal was to identify a gene-expression-based surrogate of genomic instability (GI) associated with the transformation of oral potentially malignant disorder (OPMD) into oral squamous cell carcinoma (OSCC). METHODS: GI was defined as the fraction of genome altered (FGA). Training sets included the CCLE and TCGA databases. The relevance of the enrichment score of the top correlated genes, referred to as the GIN score, was evaluated in eight independent public datasets from the GEO repository, including a cohort of patients with OPMD with available outcome. RESULTS: A set of 20 genes correlated with FGA in head and neck SCC were identified. A significant correlation was found between the 20-gene based GIN score and FGA in 95 esophagus SCC (r = 0.59) and 501 lung SCC (r = 0.63), and in 33 OPMD/OSCC (r = 0.38). A significantly increased GIN score was observed at different stages of oral carcinogenesis (normal-dysplasia -OSCC) in five independent datasets. The GIN score was higher in 10 OPMD that transformed into oral cancer compared to 10 nontransforming OPMD (p = 0.0288), and was associated with oral-cancer-free survival in 86 patients with OPMD (p = 0.0081). CONCLUSIONS: The GIN score is a gene-expression surrogate of GI, and is associated with oral carcinogenesis and OPMD malignant transformation.

Unmet Needs and Perspectives in Oral Cancer Prevention.

Oral potentially malignant disorders (OPMD) may precede oral squamous cell carcinoma (OSCC). Reported rates of malignant transformation of OPMD range from 3 to 50%. While some clinical, histological, and molecular factors have been associated with a high-risk OPMD, they are, to date, insufficiently accurate for treatment decision-making. Moreover, this range highlights differences in the clinical definition of OPMD, variation in follow-up periods, and molecular and biological heterogeneity of OPMD. Finally, while treatment of OPMD may improve outcome, standard therapy has been shown to be ineffective to prevent OSCC development in patients with OPMD. In this perspective paper, several experts discuss the main challenges in oral cancer prevention, in particular the need to (i) to define an OPMD classification system by integrating new pathological and molecular characteristics, aiming (ii) to better identify OPMD at high risk of malignant transformation, and (iii) to develop treatment strategies to eradicate OPMD or prevent malignant transformation.

Immunosuppressive niche engineering at the onset of human colorectal cancer.

The evolutionary dynamics of tumor initiation remain undetermined, and the interplay between neoplastic cells and the immune system is hypothesized to be critical in transformation. Colorectal cancer (CRC) presents a unique opportunity to study the transition to malignancy as pre-cancers (adenomas) and early-stage cancers are frequently resected. Here, we examine tumor-immune eco-evolutionary dynamics from pre-cancer to carcinoma using a computational model, ecological analysis of digital pathology data, and neoantigen prediction in 62 patient samples. Modeling predicted recruitment of immunosuppressive cells would be the most common driver of transformation. As predicted, ecological analysis reveals that progressed adenomas co-localized with immunosuppressive cells and cytokines, while benign adenomas co-localized with a mixed immune response. Carcinomas converge to a common immune "cold" ecology, relaxing selection against immunogenicity and high neoantigen burdens, with little evidence for PD-L1 overexpression driving tumor initiation. These findings suggest re-engineering the immunosuppressive niche may prove an effective immunotherapy in CRC.

Assessing Cell Activities rather than Identities to Interpret Intra-Tumor Phenotypic Diversity and Its Dynamics.

Despite advances in single-cell and molecular techniques, it is still unclear how to best quantify phenotypic heterogeneity in cancer cells that evolved beyond normal, known classifications. We present an approach to phenotypically characterize cells based on their activities rather than static classifications. We validated the detectability of specific activities (epithelial-mesenchymal transition, glycolysis) in single cells, using targeted RT-qPCR analyses and in vitro inductions. We analyzed 50 established activity signatures as a basis for phenotypic description in public data and computed cell-cell distances in 28,513 cells from 85 patients and 8 public datasets. Despite not relying on any classification, our measure correlated with standard diversity indices in populations of known structure. We identified bottlenecks as phenotypic diversity reduced upon colorectal cancer initiation. This suggests that focusing on what cancer cells do rather than what they are can quantify phenotypic diversity in universal fashion, to better understand and predict intra-tumor heterogeneity dynamics.

Comprehensive characterization of claudin-low breast tumors reflects the impact of the cell-of-origin on cancer evolution.

Claudin-low breast cancers are aggressive tumors defined by the low expression of key components of cellular junctions, associated with mesenchymal and stemness features. Although they are generally considered as the most primitive breast malignancies, their histogenesis remains elusive. Here we show that this molecular subtype of breast cancers exhibits a significant diversity, comprising three main subgroups that emerge from unique evolutionary processes. Genetic, gene methylation and gene expression analyses reveal that two of the subgroups relate, respectively, to luminal breast cancers and basal-like breast cancers through the activation of an EMT process over the course of tumor progression. The third subgroup is closely related to normal human mammary stem cells. This unique subgroup of breast cancers shows a paucity of genomic aberrations and a low frequency of TP53 mutations, supporting the emerging notion that the intrinsic properties of the cell-of-origin constitute a major determinant of the genetic history of tumorigenesis.

Separating the Local and Malignant Dimensions of Cancer Adaptation.

The repeatability observed across cancers arising in the same tissue can help understand the evolutionary process of tumour initiation. We recently developed a framework to quantify the local malignant adaptation of genetic clones in tissue-specific environments. In this Commentary, we argue that such a 1-dimensional model can be improved by separating its 2 components to obtain a dual scale: local adaptation, dictating proliferation rates in the local environment, and malignant adaptation, influencing the likelihood that a clone becomes cancerous and invasive. Such a change could strengthen our understanding of the population dynamics underlying cancer initiation and assess different evolutionary scenarios.

Quantifying local malignant adaptation in tissue-specific evolutionary trajectories by harnessing cancer's repeatability at the genetic level.

Cancer is a potentially lethal disease, in which patients with nearly identical genetic backgrounds can develop a similar pathology through distinct combinations of genetic alterations. We aimed to reconstruct the evolutionary process underlying tumour initiation, using the combination of convergence and discrepancies observed across 2,742 cancer genomes from nine tumour types. We developed a framework using the repeatability of cancer development to score the local malignant adaptation (LMA) of genetic clones, as their potential to malignantly progress and invade their environment of origin. Using this framework, we found that premalignant skin and colorectal lesions appeared specifically adapted to their local environment, yet insufficiently for full cancerous transformation. We found that metastatic clones were more adapted to the site of origin than to the invaded tissue, suggesting that genetics may be more important for local progression than for the invasion of distant organs. In addition, we used network analyses to investigate evolutionary properties at the system-level, highlighting that different dynamics of malignant progression can be modelled by such a framework in tumour-type-specific fashion. We find that occurrence-based methods can be used to specifically recapitulate the process of cancer initiation and progression, as well as to evaluate the adaptation of genetic clones to given environments. The repeatability observed in the evolution of most tumour types could therefore be harnessed to better predict the trajectories likely to be taken by tumours and preneoplastic lesions in the future.

Evolution of Barrett's esophagus through space and time at single-crypt and whole-biopsy levels.

The low risk of progression of Barrett's esophagus (BE) to esophageal adenocarcinoma can lead to over-diagnosis and over-treatment of BE patients. This may be addressed through a better understanding of the dynamics surrounding BE malignant progression. Although genetic diversity has been characterized as a marker of malignant development, it is still unclear how BE arises and develops. Here we uncover the evolutionary dynamics of BE at crypt and biopsy levels in eight individuals, including four patients that experienced malignant progression. We assay eight individual crypts and the remaining epithelium by SNP array for each of 6-11 biopsies over 2 time points per patient (358 samples in total). Our results indicate that most Barrett's segments are clonal, with similar number and inferred rates of alterations observed for crypts and biopsies. Divergence correlates with geographical location, being higher near the gastro-esophageal junction. Relaxed clock analyses show that genomic instability precedes and is enhanced by genome doubling. These results shed light on the clinically relevant evolutionary dynamics of BE.

Increased expression of the thyroid hormone nuclear receptor TRα1 characterizes intestinal tumors with high Wnt activity.

Our previous work demonstrated a key function of the thyroid hormone nuclear receptor TRα1, a T3-modulated transcription factor, in controlling intestinal development and homeostasis via the Wnt and Notch pathways. Importantly, increased expression of TRα1 in the intestinal epithelium in a mutated Apc genetic background (vil-TRα1/Apc+/1638N mice) accelerated tumorigenesis and contributed to a more aggressive tumor phenotype compared to that of the Apc mutants alone. Therefore, the aim of this study was to determine the relevance of this synergistic effect in human colorectal cancers and to gain insights into the mechanisms involved. We analyzed cohorts of patients by in silico and experimental approaches and observed increased TRα1 expression and a significant correlation between TRα1 levels and Wnt activity. TRα1 loss-of-function and gain-of-function in Caco2 cell lines not only confirmed that TRα1 levels control Wnt activity but also demonstrated the role of TRα1 in regulating cell proliferation and migration. Finally, upon investigation of the molecular mechanisms responsible for the Wnt-TRα1 association, we described the repression by TRα1 of several Wnt inhibitors, including Frzb, Sox17 and Wif1. In conclusion, our results underline an important functional interplay between the thyroid hormone nuclear receptor TRα1 and the canonical Wnt pathway in intestinal cancer initiation and progression. More importantly, we show for the first time that the expression of TRα1 is induced in human colorectal cancers.