Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal.

Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5' UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor's most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

Tracking Cancer Evolution Reveals Constrained Routes to Metastases: TRACERx Renal.

Clear-cell renal cell carcinoma (ccRCC) exhibits a broad range of metastatic phenotypes that have not been systematically studied to date. Here, we analyzed 575 primary and 335 metastatic biopsies across 100 patients with metastatic ccRCC, including two cases sampledat post-mortem. Metastatic competence was afforded by chromosome complexity, and we identify 9p loss as a highly selected event driving metastasis and ccRCC-related mortality (p = 0.0014). Distinct patterns of metastatic dissemination were observed, including rapid progression to multiple tissue sites seeded by primary tumors of monoclonal structure. By contrast, we observed attenuated progression in cases characterized by high primary tumor heterogeneity, with metastatic competence acquired gradually and initial progression to solitary metastasis. Finally, we observed early divergence of primitive ancestral clones and protracted latency of up to two decades as a feature of pancreatic metastases.

Neuronal programming by microbiota regulates intestinal physiology.

Neural control of the function of visceral organs is essential for homeostasis and health. Intestinal peristalsis is critical for digestive physiology and host defence, and is often dysregulated in gastrointestinal disorders1. Luminal factors, such as diet and microbiota, regulate neurogenic programs of gut motility2-5, but the underlying molecular mechanisms remain unclear. Here we show that the transcription factor aryl hydrocarbon receptor (AHR) functions as a biosensor in intestinal neural circuits, linking their functional output to the microbial environment of the gut lumen. Using nuclear RNA sequencing of mouse enteric neurons that represent distinct intestinal segments and microbiota states, we demonstrate that the intrinsic neural networks of the colon exhibit unique transcriptional profiles that are controlled by the combined effects of host genetic programs and microbial colonization. Microbiota-induced expression of AHR in neurons of the distal gastrointestinal tract enables these neurons to respond to the luminal environment and to induce expression of neuron-specific effector mechanisms. Neuron-specific deletion of Ahr, or constitutive overexpression of its negative feedback regulator CYP1A1, results in reduced peristaltic activity of the colon, similar to that observed in microbiota-depleted mice. Finally, expression of Ahr in the enteric neurons of mice treated with antibiotics partially restores intestinal motility. Together, our experiments identify AHR signalling in enteric neurons as a regulatory node that integrates the luminal environment with the physiological output of intestinal neural circuits to maintain gut homeostasis and health.

Pervasive chromosomal instability and karyotype order in tumour evolution.

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes1,2. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution1,3,4. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2+ breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Author Correction: Metastatic-niche labelling reveals parenchymal cells with stem features.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Metastatic-niche labelling reveals parenchymal cells with stem features.

Direct investigation of the early cellular changes induced by metastatic cells within the surrounding tissue remains a challenge. Here we present a system in which metastatic cancer cells release a cell-penetrating fluorescent protein, which is taken up by neighbouring cells and enables spatial identification of the local metastatic cellular environment. Using this system, tissue cells with low representation in the metastatic niche can be identified and characterized within the bulk tissue. To highlight its potential, we applied this strategy to study the cellular environment of metastatic breast cancer cells in the lung. We report the presence of cancer-associated parenchymal cells, which exhibit stem-cell-like features, expression of lung progenitor markers, multi-lineage differentiation potential and self-renewal activity. In ex vivo assays, lung epithelial cells acquire a cancer-associated parenchymal-cell-like phenotype when co-cultured with cancer cells and support their growth. These results highlight the potential of this method as a platform for new discoveries.

ERK-Induced Activation of TCF Family of SRF Cofactors Initiates a Chromatin Modification Cascade Associated with Transcription.

We investigated the relationship among ERK signaling, histone modifications, and transcription factor activity, focusing on the ERK-regulated ternary complex factor family of SRF partner proteins. In MEFs, activation of ERK by TPA stimulation induced a common pattern of H3K9acS10ph, H4K16ac, H3K27ac, H3K9acK14ac, and H3K4me3 at hundreds of transcription start site (TSS) regions and remote regulatory sites. The magnitude of the increase in histone modification correlated well with changes in transcription. H3K9acS10ph preceded the other modifications. Most induced changes were TCF dependent, but TCF-independent TSSs exhibited the same hierarchy, indicating that it reflects gene activation per se. Studies with TCF Elk-1 mutants showed that TCF-dependent ERK-induced histone modifications required Elk-1 to be phosphorylated and competent to activate transcription. Analysis of direct TCF-SRF target genes and chromatin modifiers confirmed this and showed that H3S10ph required only Elk-1 phosphorylation. Induction of histone modifications following ERK stimulation is thus directed by transcription factor activation and transcription.

The next-generation Open Targets Platform: reimagined, redesigned, rebuilt.

The Open Targets Platform (https://platform.opentargets.org/) is an open source resource to systematically assist drug target identification and prioritisation using publicly available data. Since our last update, we have reimagined, redesigned, and rebuilt the Platform in order to streamline data integration and harmonisation, expand the ways in which users can explore the data, and improve the user experience. The gene-disease causal evidence has been enhanced and expanded to better capture disease causality across rare, common, and somatic diseases. For target and drug annotations, we have incorporated new features that help assess target safety and tractability, including genetic constraint, PROTACtability assessments, and AlphaFold structure predictions. We have also introduced new machine learning applications for knowledge extraction from the published literature, clinical trial information, and drug labels. The new technologies and frameworks introduced since the last update will ease the introduction of new features and the creation of separate instances of the Platform adapted to user requirements. Our new Community forum, expanded training materials, and outreach programme support our users in a range of use cases.

Validation of epidermal AMBRA1 and loricrin (AMBLor) as a prognostic biomarker for nonulcerated American Joint Committee on Cancer stage I/II cutaneous melanoma.

BACKGROUND: Combined expression of the autophagy-regulatory protein AMBRA1 (activating molecule in Beclin1-regulated autophagy) and the terminal differentiation marker loricrin in the peritumoral epidermis of stage I melanomas can identify tumour subsets at low risk of -metastasis. OBJECTIVES: To validate the combined expression of peritumoral AMBRA1 and loricrin (AMBLor) as a prognostic biomarker able to identify both stage I and II melanomas at low risk of tumour recurrence. METHODS: Automated immunohistochemistry was used to analyse peritumoral AMBRA1 and loricrin expression in geographically distinct discovery (n = 540) and validation (n = 300) cohorts of nonulcerated American Joint Committee on Cancer (AJCC) stage I and II melanomas. AMBLor status was correlated with clinical outcomes in the discovery and validation cohorts separately and combined. RESULTS: Analysis of AMBLor in the discovery cohort revealed a recurrence-free survival (RFS) rate of 95.5% in the AMBLor low-risk group vs. 81.7% in the AMBLor at-risk group (multivariate log-rank, P < 0.001) and a negative predictive value (NPV) of 96.0%. In the validation cohort, AMBLor analysis revealed a RFS rate of 97.6% in the AMBLor low-risk group vs. 78.3% in the at-risk group (multivariate log-rank, P < 0.001) and a NPV of 97.6%. In a multivariate model considering AMBLor, Breslow thickness, age and sex, analysis of the combined discovery and validation cohorts showed that the estimated effect of AMBLor was statistically significant, with a hazard ratio of 3.469 (95% confidence interval 1.403-8.580, P = 0.007) and an overall NPV of 96.5%. CONCLUSIONS: These data provide further evidence validating AMBLor as a prognostic biomarker to identify nonulcerated AJCC stage I and II melanoma tumours at low risk of disease recurrence.

SRF Co-factors Control the Balance between Cell Proliferation and Contractility.

The ERK-regulated ternary complex factors (TCFs) act with the transcription factor serum response factor (SRF) to activate mitogen-induced transcription. However, the extent of their involvement in the immediate-early transcriptional response, and their wider functional significance, has remained unclear. We show that, in MEFs, TCF inactivation significantly inhibits over 60% of TPA-inducible gene transcription and impairs cell proliferation. Using integrated SRF ChIP-seq and Hi-C data, we identified over 700 TCF-dependent SRF direct target genes involved in signaling, transcription, and proliferation. These also include a significant number of cytoskeletal gene targets for the Rho-regulated myocardin-related transcription factor (MRTF) SRF cofactor family. The TCFs act as general antagonists of MRTF-dependent SRF target gene expression, competing directly with the MRTFs for access to SRF. As a result, TCF-deficient MEFs exhibit hypercontractile and pro-invasive behavior. Thus, competition between TCFs and MRTFs for SRF determines the balance between antagonistic proliferative and contractile programs of gene expression.

Mutation of cancer driver MLL2 results in transcription stress and genome instability.

Genome instability is a recurring feature of tumorigenesis. Mutation in MLL2, encoding a histone methyltransferase, is a driver in numerous different cancer types, but the mechanism is unclear. Here, we present evidence that MLL2 mutation results in genome instability. Mouse cells in which MLL2 gene deletion can be induced display elevated levels of sister chromatid exchange, gross chromosomal aberrations, 53BP1 foci, and micronuclei. Human MLL2 knockout cells are characterized by genome instability as well. Interestingly, MLL2 interacts with RNA polymerase II (RNAPII) and RECQL5, and, although MLL2 mutated cells have normal overall H3K4me levels in genes, nucleosomes in the immediate vicinity of RNAPII are hypomethylated. Importantly, MLL2 mutated cells display signs of substantial transcription stress, and the most affected genes overlap with early replicating fragile sites, show elevated levels of γH2AX, and suffer frequent mutation. The requirement for MLL2 in the maintenance of genome stability in genes helps explain its widespread role in cancer and points to transcription stress as a strong driver in tumorigenesis.

Rho-actin signaling to the MRTF coactivators dominates the immediate transcriptional response to serum in fibroblasts.

The transcription factor SRF (serum response factor) recruits two families of coactivators, the MRTFs (myocardin-related transcription factors) and the TCFs (ternary complex factors), to couple gene transcription to growth factor signaling. Here we investigated the role of the SRF network in the immediate transcriptional response of fibroblasts to serum stimulation. SRF recruited its cofactors in a gene-specific manner, and virtually all MRTF binding was directed by SRF. Much of SRF DNA binding was serum-inducible, reflecting a requirement for MRTF-SRF complex formation in nucleosome displacement. We identified 960 serum-responsive SRF target genes, which were mostly MRTF-controlled, as assessed by MRTF chromatin immunoprecipitation (ChIP) combined with deep sequencing (ChIP-seq) and/or sensitivity to MRTF-linked signals. MRTF activation facilitates RNA polymerase II (Pol II) recruitment or promoter escape according to gene context. MRTF targets encode regulators of the cytoskeleton, transcription, and cell growth, underpinning the role of SRF in cytoskeletal dynamics and mechanosensing. Finally, we show that specific activation of either MRTFs or TCFs can reset the circadian clock.

Inherited duplications of PPP2R3B predispose to nevi and melanoma via a C21orf91-driven proliferative phenotype.

PURPOSE: Much of the heredity of melanoma remains unexplained. We sought predisposing germline copy-number variants using a rare disease approach. METHODS: Whole-genome copy-number findings in patients with melanoma predisposition syndrome congenital melanocytic nevus were extrapolated to a sporadic melanoma cohort. Functional effects of duplications in PPP2R3B were investigated using immunohistochemistry, transcriptomics, and stable inducible cellular models, themselves characterized using RNAseq, quantitative real-time polymerase chain reaction (qRT-PCR), reverse phase protein arrays, immunoblotting, RNA interference, immunocytochemistry, proliferation, and migration assays. RESULTS: We identify here a previously unreported genetic susceptibility to melanoma and melanocytic nevi, familial duplications of gene PPP2R3B. This encodes PR70, a regulatory unit of critical phosphatase PP2A. Duplications increase expression of PR70 in human nevus, and increased expression in melanoma tissue correlates with survival via a nonimmunological mechanism. PPP2R3B overexpression induces pigment cell switching toward proliferation and away from migration. Importantly, this is independent of the known microphthalmia-associated transcription factor (MITF)-controlled switch, instead driven by C21orf91. Finally, C21orf91 is demonstrated to be downstream of MITF as well as PR70. CONCLUSION: This work confirms the power of a rare disease approach, identifying a previously unreported copy-number change predisposing to melanocytic neoplasia, and discovers C21orf91 as a potentially targetable hub in the control of phenotype switching.

Complement pathway gene activation and rising circulating immune complexes characterize early disease in HIV-associated tuberculosis.

The transition between latent and active tuberculosis (TB) occurs before symptom onset. Better understanding of the early events in subclinical disease will facilitate the development of diagnostics and interventions that improve TB control. This is particularly relevant in the context of HIV-1 coinfection where progression of TB is more likely. In a recent study using [18F]-fluoro-2-deoxy-d-glucose positron emission/computed tomography (FDG-PET/CT) on 35 asymptomatic, HIV-1-infected adults, we identified 10 participants with radiographic evidence of subclinical disease, significantly more likely to progress than the 25 participants without. To gain insight into the biological events in early disease, we performed blood-based whole genome transcriptomic analysis on these participants and 15 active patients with TB. We found transcripts representing the classical complement pathway and Fcγ receptor 1 overabundant from subclinical stages of disease. Levels of circulating immune (antibody/antigen) complexes also increased in subclinical disease and were highly correlated with C1q transcript abundance. To validate our findings, we analyzed transcriptomic data from a publicly available dataset where samples were available in the 2 y before TB disease presentation. Transcripts representing the classical complement pathway and Fcγ receptor 1 were also differentially expressed in the 12 mo before disease presentation. Our results indicate that levels of antibody/antigen complexes increase early in disease, associated with increased gene expression of C1q and Fcγ receptors that bind them. Understanding the role this plays in disease progression may facilitate development of interventions that prevent this, leading to a more favorable outcome and may also be important to diagnostic development.

Mycobacterium tuberculosis requires glyoxylate shunt and reverse methylcitrate cycle for lactate and pyruvate metabolism.

Bacterial nutrition is an essential aspect of host-pathogen interaction. For the intracellular pathogen Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis in humans, fatty acids derived from lipid droplets are considered the major carbon source. However, many other soluble nutrients are available inside host cells and may be used as alternative carbon sources. Lactate and pyruvate are abundant in human cells and fluids, particularly during inflammation. In this work, we study Mtb metabolism of lactate and pyruvate combining classic microbial physiology with a 'multi-omics' approach consisting of transposon-directed insertion site sequencing (TraDIS), RNA-seq transcriptomics, proteomics and stable isotopic labelling coupled with mass spectrometry-based metabolomics. We discovered that Mtb is well adapted to use both lactate and pyruvate and that their metabolism requires gluconeogenesis, valine metabolism, the Krebs cycle, the GABA shunt, the glyoxylate shunt and the methylcitrate cycle. The last two pathways are traditionally associated with fatty acid metabolism and, unexpectedly, we found that in Mtb the methylcitrate cycle operates in reverse, to allow optimal metabolism of lactate and pyruvate. Our findings reveal a novel function for the methylcitrate cycle as a direct route for the biosynthesis of propionyl-CoA, the essential precursor for the biosynthesis of the odd-chain fatty acids.

Tracking the Evolution of Non-Small-Cell Lung Cancer.

BACKGROUND: Among patients with non-small-cell lung cancer (NSCLC), data on intratumor heterogeneity and cancer genome evolution have been limited to small retrospective cohorts. We wanted to prospectively investigate intratumor heterogeneity in relation to clinical outcome and to determine the clonal nature of driver events and evolutionary processes in early-stage NSCLC. METHODS: In this prospective cohort study, we performed multiregion whole-exome sequencing on 100 early-stage NSCLC tumors that had been resected before systemic therapy. We sequenced and analyzed 327 tumor regions to define evolutionary histories, obtain a census of clonal and subclonal events, and assess the relationship between intratumor heterogeneity and recurrence-free survival. RESULTS: We observed widespread intratumor heterogeneity for both somatic copy-number alterations and mutations. Driver mutations in EGFR, MET, BRAF, and TP53 were almost always clonal. However, heterogeneous driver alterations that occurred later in evolution were found in more than 75% of the tumors and were common in PIK3CA and NF1 and in genes that are involved in chromatin modification and DNA damage response and repair. Genome doubling and ongoing dynamic chromosomal instability were associated with intratumor heterogeneity and resulted in parallel evolution of driver somatic copy-number alterations, including amplifications in CDK4, FOXA1, and BCL11A. Elevated copy-number heterogeneity was associated with an increased risk of recurrence or death (hazard ratio, 4.9; P=4.4×10-4), which remained significant in multivariate analysis. CONCLUSIONS: Intratumor heterogeneity mediated through chromosome instability was associated with an increased risk of recurrence or death, a finding that supports the potential value of chromosome instability as a prognostic predictor. (Funded by Cancer Research UK and others; TRACERx ClinicalTrials.gov number, NCT01888601 .).

Small RNA profiling of Xenopus embryos reveals novel miRNAs and a new class of small RNAs derived from intronic transposable elements.

Small RNA control of gene expression is critical for developmental processes in vertebrate embryos. To determine the dynamics of small RNA expression and to uncover novel small RNAs in the early vertebrate embryo, we performed high-throughput sequencing of all small RNAs in Xenopus tropicalis embryos at three developmental time points and in dissected halves of gastrula embryos. This analysis allowed us to identify novel microRNAs and we show that microRNA expression is highly dynamic and spatially localized in early embryos. In addition, we have developed a microRNA prediction pipeline and demonstrate that it has the power to predict new miRNAs that are experimentally detectable in frogs, mice, and humans. By combining the small RNA sequencing with mRNA profiling at the different developmental stages, we identify a new class of small noncoding RNAs that we name siteRNAs, which align in clusters to introns of protein-coding genes. We show that siteRNAs are derived from remnants of transposable elements present in the introns. We find that genes containing clusters of siteRNAs are transcriptionally repressed as compared with all genes. Furthermore, we show that this is true for individual genes containing siteRNA clusters, and that these genes are enriched in specific repressive histone modifications. Our data thus suggest a new mechanism of siteRNA-mediated gene silencing in vertebrates, and provide an example of how mobile elements can affect gene regulation.

Tracking genomic cancer evolution for precision medicine: the lung TRACERx study.

The importance of intratumour genetic and functional heterogeneity is increasingly recognised as a driver of cancer progression and survival outcome. Understanding how tumour clonal heterogeneity impacts upon therapeutic outcome, however, is still an area of unmet clinical and scientific need. TRACERx (TRAcking non-small cell lung Cancer Evolution through therapy [Rx]), a prospective study of patients with primary non-small cell lung cancer (NSCLC), aims to define the evolutionary trajectories of lung cancer in both space and time through multiregion and longitudinal tumour sampling and genetic analysis. By following cancers from diagnosis to relapse, tracking the evolutionary trajectories of tumours in relation to therapeutic interventions, and determining the impact of clonal heterogeneity on clinical outcomes, TRACERx may help to identify novel therapeutic targets for NSCLC and may also serve as a model applicable to other cancer types.

The Wilms Tumor-1 (WT1) rs2234593 variant is a prognostic factor in normal karyotype acute myeloid leukemia.

The single-nucleotide polymorphism (SNP) within Wilms tumor-1 (WT1) exon 7, rs16754, has been arguably reported to be implicated in acute myeloid leukemia (AML) prognosis. We assessed the potential association of selected WT1 SNPs as well as WT1 mutations in normal karyotype (NK)-AML and evaluated the prognostic value of these normal gene variants. Diagnostic samples from a series of 474 young adult NK-AML patients were used to genotype five WT1 SNPs using TaqMan assays and to directly sequence WT1 exons 7 and 9. Analysis of five WT1 gene variants showed an association of rs2234593 allele C with WT1 Ex7 mutation. Prognostic study of the same variants identified rs2234593 significantly associated with relapse and overall survival (OS). Patients with rs2234593AA/AC showed significantly higher 10-year OS (50 vs 36 %, hazard ratio (HR) = 0.69 (0.52­0.90), p = 0.006) and lower cumulative incidence of relapse (CIR) (36 vs 51 %, HR = 0.62 (0.45­0.86), p = 0.004) compared to those with rs2234593CC. The effect of AA genotype on CIR remained significant after adjustment for basic covariates including FLT3 internal-tandem duplication (FLT3-ITD) and nucleophosmin 1 (NPM1) mutations (HR = 0.60 (0.41­0.89), p = 0.009), with some evidence of improved survival (HR = 0.75 (0.55­1.03), p = 0.07). A multivariate analysis showed WT1 Ex7-mutant as the major relapse predictor, with a tendency for rs2234593-A effect after allowing for Ex7 mutation (p = 0.07). No adjusted risk benefit was found for previously reported rs16754-G. In conclusion, WT1 normal gene variant rs2234593 is associated with mutational status of WT1 Ex7 and is a further prognostic marker independent from FLT3-ITD and NPM1 mutations in NK-AML.